- 1Biomedical Engineering Research Centre, Division of Biomedical Engineering, Department of Human Biology, University of Cape Town, Cape Town, South Africa
- 2Neurosciences Institute, University of Cape Town, Cape Town, South Africa
- 3Cape Universities Body Imaging Centre, University of Cape Town, Cape Town, South Africa
- 4School of Speech, Language, and Hearing Sciences, College of Health and Human Services, San Diego, CA, United States
- 5Family Centre for Research with Ubuntu, Department of Pediatrics and Child Health, Stellenbosch University, Cape Town, South Africa
- 6Department of Biomedical Physiology and Kinesiology, Simon Fraser University, Burnaby, BC, Canada
- 7ImageTech, Simon Fraser University, Surrey, BC, Canada
Introduction: Children living with perinatally acquired HIV (CPHIV) demonstrate hearing impairments and language processing delays even in the presence of combination antiretroviral therapy (cART). Investigations on the effect of HIV on the auditory system have predominantly focused on the peripheral auditory system. Additionally, language processing requires the efficient interaction between central auditory system (CAS) brain regions and non-auditory regions. Investigating the functional connectivity (FC) within the CAS and between the CAS and non-auditory regions may reveal the influence of HIV on regions involved in auditory function.
Methods: Within a Bayesian statistical framework, we used resting-state functional magnetic resonance imaging to map FC in the CAS as well as between CAS regions and non-auditory regions of 11-year-old CPHIV. Graph theory was used to investigate the regional effects of HIV on brain network properties. We explored the relationships between FC and neurocognitive outcomes. We hypothesized that CPHIV would show disruptions in FC between CAS regions as well as between CAS and non-auditory regions. Secondly, we hypothesized that in CPHIV, regional brain network properties would be altered compared to their uninfected peers (CHUU). Finally we hypothesized that FC and functional network regional outcomes would be related to neurocognitive outcomes.
Results: Our investigation revealed lower FC of the primary auditory cortex (PAC) in CPHIV as well as disruptions in FC between CAS regions and non-auditory regions including hippocampal sub-regions, the lingual gyri and basal ganglia. Functional network analysis revealed lower nodal degree and efficiency in CAS regions including the cochlear nucleus/superior olivary complex and the inferior colliculus. We also report associations between the nodal efficiency of middle temporal and superior frontal regions and delayed recall, a neurocognitive marker of working memory, present in CHUU but not in CPHIV.
Discussion: Our results demonstrate FC alterations in the PAC and between CAS regions and non-auditory regions involved in limbic, visual and motor processing, as well as disruptions to the regional properties of the CAS regions in the functional brain network. These results provide insight into the state of the CAS FC in the presence of HIV and its possible role in the hearing and language impairments seen in this population.
Introduction
The efficient processing of sound involves both the perception of sound by the peripheral nervous system and the processing of auditory signals in the central nervous system (CNS) (Litovsky, 2015). While basic auditory abilities are developed at birth, continuous exposure to auditory stimulation in childhood results in functional changes in the central auditory system (CAS) (Fitzroy et al., 2015; Skoe et al., 2015). The CAS, which comprises the cochlear nucleus (CN) and inferior colliculus (IC) in the brainstem, medial geniculate nucleus (MGN) in the thalamus and primary auditory cortex (PAC) in the superior temporal gyrus (STG), interacts heavily with other brain areas, including the temporal and frontal lobe regions, for efficient language processing and communication (Maffei et al., 2015; Middlebrooks et al., 2016).
Children perinatally infected with the human immunodeficiency virus (CPHIV) are at risk of hearing impairments and disruptions in language development (Chao et al., 2012; Maro et al., 2017; Rice et al., 2012; Torre et al., 2012, 2015). Associations between auditory function and HIV status have been reported (Niemczak et al., 2021). Further, CPHIV are at risk of lower working memory and verbal fluency outcomes (Lima et al., 2023; Morgan et al., 2015). Investigating the function of the CAS may provide a window into the relationship between HIV, language and the brain in these at-risk populations.
Resting-state functional magnetic resonance imaging (RS-fMRI) is a tool to investigate co-activation of anatomically separate brain areas, which reflects functional communication between areas (van den Heuvel and Hulshoff Pol, 2010). In the literature, functional co-activation is often referred to as functional connectivity (FC). Coupled with a graph theory analysis, RS-fMRI can provide insight into information transfer and integration in the functional brain network. In a RS-fMRI based graph-theoretical framework, a graph is a set of brain regions (nodes) functionally connected by edges. Using this approach, we can investigate network measures of integration and segregation at the regional/nodal level as well as at a global level revealing the arrangement or topology of the network and its regions (Bullmore et al., 2016).
To date, there are only a handful of RS-fMRI studies in the context of children living with HIV and none focused on the CAS. Yadav et al. (2018) reported lower FC in the middle temporal gyrus MTG of CPHIV, which could potentially affect auditory processing (Yadav et al., 2018). Another study showed HIV-related altered FC between brain regions in the DMN and regions belonging to resting-state networks (RSNs) that work with the auditory network, including the visual and somatosensory networks (Herting et al., 2015). Our group has used neuroimaging techniques to study brain development since age 5 years in a cohort of CPHIV who initiated antiretroviral therapy (ART) in infancy. At 7 years old, we observed reduced FC in the inferior frontal gyrus of CPHIV, a region that is associated with the auditory RSN (Toich et al., 2018). Although these findings allude to potential alterations in auditory functional regions in CPHIV, investigating the specific regions involved in auditory function as well as their interactions with other brain regions that participate in language processing and communication is necessary.
We used RS-fMRI and a graph-theoretical framework to assess FC of the CAS in a cohort of 11-year old children. We hypothesized significant differences in FC in CPHIV as compared to CHUU children. Furthermore, we posited that the CAS in CPHIV would have reduced functional strength and efficiency compared to CHUU. Lastly, to better understand the possible relationship to auditory outcomes, we explored associations between functional deficits and cognitive outcomes related to language and hearing.
Methods
Study cohort
Participants were 84 (11–12-year-old) children (25 CHUU, 59 CPHIV) who completed structural and RS-fMRI scans. The children living with HIV were recruited from the randomized Children with HIV Early Antiretroviral Therapy (CHER) trial and have been monitored since birth at the Family Centre for Research with Ubuntu (FAMCRU), Tygerberg Children’s Hospital, in Cape Town, South Africa (Cotton et al., 2014; Violari et al., 2008). For the CHER trial, infants were randomized into one of three treatment groups; two arms involved early (before 12 weeks of age) ART and one deferred ART until CD4% < 25% in the first year or CD4 < 20% thereafter, or if clinical disease progression criteria presented (Cotton et al., 2014). The uninfected children were recruited from a parallel vaccine study (Madhi et al., 2010).
T1-weighted (T1w) structural and RS-fMRI image acquisition
T1w structural and RS-fMRI acquisitions were done in the same session using a 3 Tesla Siemens Skyra MRI scanner (Erlangen, Germany) at University of Cape Town’s Cape Universities Body Imaging Centre (CUBIC). Children were first familiarized with the scanning procedures on a mock scanner. The T1w structural was acquired for 6 min using a multiecho magnetization prepared rapid gradient echo (ME-MPRAGE) sequence with TR = 2,530 ms, TEs = 1.69/3.54/5.39/7.24 ms, inversion time (TI) 1,100 ms, flip angle = 7o, and FOV = 224 × 224 × 176 mm3. RS-fMRI data were acquired for 6 min using an interleaved multi-slice 2D gradient echo, echo planar imaging (EPI) sequence with the following parameters: 33 interleaved slices, slice thickness 4 mm, slice gap 1 mm, voxel size 3 × 3 × 4 mm3, FOV = 250 × 250 × 164 mm3, TR/TE = 2,000/30 ms, flip angle 90°, 180 EPI volumes.
Neurodevelopmental measures: KABC-II auditory/language-related scores
Participants were administered the Kaufman Assessment Battery for Children – Second Edition (KABC-II) (Kaufman and Kaufman, 2014) by trained research assistants supervised by a licensed psychologist. The neurodevelopmental scores were part of the Kaufman Assessment Battery for Children – Second Edition (KABC-II) (Kaufman and Kaufman, 2014) and were measured in the same group of children at 11 years old. Subtests of the battery include Learning, (Atlantis, Atlantis delayed, Rebus, Rebus delayed, Delayed Recall) which assesses the ability to store and retrieve information, Planning (Story completion and Pattern reasoning) which assesses the ability to solve nonverbal problems, Sequential processing (Number recall, Word order and Hand movement) which assesses ability to solve problems using auditory working memory and visual information presented serially, and Simultaneous processing (Rover, Block Counting, Conceptual Thinking, Triangles, Face Recognition, and Gestalt Closure) which assesses the ability to solve spatial or logistical problems that require simultaneous processing of multiple related stimuli (Kaufman and Kaufman, 2014). The sum of all scores generates a global cognitive score called the Mental Processing index. Moreover, summing the subtests that do not require verbal expression (Hand Movements, Block Counting, Triangles, Pattern Reasoning, Story Completion, Conceptual Thinking, and Face Recognition) generates the Nonverbal Index score (Kaufman and Kaufman, 2014). For the current study we chose the tests related to auditory and language processing: Number recall, Word order, Learning, Atlantis, Atlantis delayed, Rebus, Rebus delayed and Delayed Recall, Story completion, Sequential processing, Mental processing and Nonverbal index. Moreover, we included Animal naming, a neurocognitive score relating to word generation and categorization (Strauss et al., 2010).
Ethical considerations
The study protocol was approved by the Human Research Ethics Committees of the University of Cape Town and Stellenbosch University. Parents/guardians provided written informed consent and children oral assent.
RS-fMRI pre-processing
The process described here was done for each subject in the Analysis of Functional Neuroimages (AFNI) software. The first 5 EPI volumes were removed to account for the time it takes for signal stabilization. Then we ‘despiked’ each volume to truncate the signal in voxels with artificially high signal, and performed temporal interpolation to adjust the acquisition time of each slice to that of the first slice of the volume. Next each EPI volume was co-registered to the reference EPI volume (6th EPI volume in the current study), generating 6 motion parameters (3 rotations and 3 translations) per volume. The motion estimates were de-meaned by removing the mean per run, and per-TR motion derivates were computed for each volume as the backward difference between the motion estimates of said volume and the prior volume [i.e., Δx(t) = x(t)-x(t-1), where x is any one of the 6 motion parameters]. The EPI volumes were spatially smoothed with a Gaussian kernel of 6 mm full width at half maximum (FWHM), and a brain mask of the EPI dataset created. The T1w structural volume was next segmented into white matter (WM), gray matter (GM) and cerebrospinal fluid (CSF). We then aligned the subject’s EPI volumes to his/her T1w structural volume, and used the resulting co-registered CSF and WM masks to extract WM and CSF signals from the EPI data. Finally, signals from non-neuronal sources were removed from the EPI data by regressing out the mean WM and CSF signals, as well as the 6 de-meaned motion parameters and their per-TR derivatives. We censored pairs of TRs where the Euclidean norm (i.e., the square root of the sum of squares) of the per-TR motion derivatives exceeded 0.3 mm, as well as TRs for which 5% or more of the brain mask was classified as outliers by afni_proc. Notably, a per-TR Euclidean norm of 0.3 mm corresponds to 0.3 mm motion in a single direction, or translations of 0.12 mm and rotations of 0.11° in each direction. Participants with fewer than 160 EPI volumes after truncation were excluded from further analyses. EPI time series data of remaining participants were truncated to 160 volumes.
GM regions of interest
To fully examine the CAS, we included both manually traced and automatically segmented regions of interest (ROI). ROIs included 4 bilateral manually traced GM regions that form part of the auditory pathway and 126 automatically segmented cortical and subcortical GM regions. The manual tracings were performed on the T1w structural volumes by an expert neuroanatomist using Freeview software FreeSurfer v7.1.0 image analysis suite: (http://surfer.nmr.mgh.harvard.edu/) on a Lenovo ThinkPad Yoga370 tablet. The 4 manually traced auditory ROIs included a region at the junction between brainstem pons and medulla and consists of the cochlear nucleus and superior olivary complex (CN/SOC), the IC in the midbrain, the MGN in the thalamus and the PAC. For automatic parcellation we used Connectome mapper (CMP) 3 (v3.0.0-RC4), an imaging processing pipeline to produce multi-resolution morphological subject atlases (Daducci et al., 2012; Tourbier et al., 2022). We used the lowest resolution which consisted of a set of 126 subcortical and cortical ROIs. A list of all the ROIs in the atlas and their abbreviations can be found in Supplementary Appendix A. We then merged the manually traced and automatically segmented ROIs. Automatically segmented ROIs that overlapped with manually traced ROIs were excluded. The number of excluded ROIs varied between 81 and 122 across subjects. Finally, 80 ROIs were common to all participants. Masks of these ROIs were mapped to the RS-fMRI mask and used as seeds for the FC analyses.
FC and functional network definition
The process described here was done for each subject. The time series signals within each seed were extracted to compute the mean time series for each ROI. We calculated a Pearson correlation coefficient between the mean time series of each ROI and every other ROI (3dNetCorr) (Taylor and Saad, 2013). The correlation coefficients were Fisher Z-transformed to obtain approximately normally distributed correlations. This resulted in an 80×80 correlation matrix for each subject with 3,160 connections. Of those connections, 604 pairs included at least one auditory ROI. This matrix was the input for our FC analysis.
To define each subject’s functional network, we discarded weak and/or negative edge weights (Fisher Z-transformed r values < 0.1003 equivalent to Pearson r values < 0.1000) from each subject’s 80×80 Fisher Z-transformed matrix (Figure 1). Removing negative edge weights allows for an intuitive interpretation of the relationship between nodes, i.e., strong positive weights mean nodes work together in a correlated way (Fornito et al., 2016; Rubinov and Sporns, 2010). We then computed nodal graph measures using iGraph and braingraph packages in R (v 4.1.3). We computed nodal degree and nodal strength – measures of the number and weightings of connections to a node, respectively. Nodes with high degree or strength tend to be identified as network hubs as they are generally positioned where they can make strong contributions to the global network (Hoff et al., 2013; van den Heuvel and Sporns, 2013). In addition we computed nodal transitivity – a measure of densely connected groups of nodes. Nodes with high transitivity in the brain are surrounded by highly interconnected nodes forming clusters that may translate as functional coordination (Fornito et al., 2016). Finally we computed nodal efficiency and nodal local efficiency. Nodal efficiency is a measure of the efficiency of communication between a node and all the other nodes in the network and nodes that have high degree tend to have high efficiency (Fornito et al., 2016; Sporns, 2018). On the other hand, nodal local efficiency is a measure of the efficiency of the graph formed in the absence of the node. This measure shows the graph’s level of fault tolerance and nodes with high nodal local efficiency reflect efficient integration with their neighbors (Fornito et al., 2016).

Figure 1. Functional network of a random control subject: Fisher Z-transformed r values of weak and negative connections (Fisher Z-transformed r value < 0.1003) were discarded. A – Full 80×80 RS connectivity matrix, B – the columns containing auditory ROIs. From left column to right: (R) PAC, (R) MGN, (R) IC, (R) SOC/CN, (L)_ SOC/CN, (L) IC, (L) MGN, (L) PAC. (R/L) – Right/Left.
Statistical analyses
To investigate the effects of living with HIV on both FC and functional network organization we employed Bayesian multilevel modeling (BML) (Chen et al., 2019a; Chen et al., 2019b). An advantage of the Bayesian framework is its non-dichotomous interpretation of results, unlike binary significance of the p-value statistic in conventional null-hypothesis significance testing (NHST) (Kruschke and Liddell, 2018). Briefly, consider a scenario where we are investigating the difference between two group means (effect), CHUU-CPHIV. Under the null hypothesis (H0), the effect equals 0. A Bayesian analysis estimates a probability distribution (posterior distribution) of the alternative hypothesis (H1: effect > 0), based on the data, i.e., P(effect>0|data). The data is assigned a prior probability distribution based on an assumption, such as that the ROIs are connected to each other through Gaussian or Student’s t-distribution. The area under the posterior distribution curve that is greater than 0, the posterior probability (P+), is interpreted as the strength of evidence that the effect is positive, given the data. That is, the strength of evidence that the reference group (CHUU) has a higher mean than the test group (CPHIV), based on the data. A very high P + represents a strong probability that the effect is greater than 0, while a very low P + means there is a weak probability the effect is greater than 0. Note, a very low P + corresponds with a strong probability of the effect being less than 0. Figure 2 shows the posterior probability distribution for the nodal strength of a random brain ROI, for the comparison between CHUU and CPHIV.

Figure 2. An example of posterior probability distribution of the nodal strength of a randomly chosen brain ROI for the comparison CHUU-CPHIV. The horizontal axes indicates the effect size, the vertical axes (black) indicates zero effect or difference between the two group means and the vertical green line indicates the mean difference. Here, the area under the curve greater than 0 (red), P + = 0.60, indicates we found with 60% probability that the CHUU group has a bigger mean than CPHIV given the data.
Furthermore, using the framework, a multilevel model is formed which integrates variability estimates at the group and subject levels in estimating the uncertainty of effect at each brain region or ROI. This integrative model is done via Monte-Carlo Markov chains (MCMC) at a predetermined number of iterations resulting in the partial pooling of the effect of interest toward the centre (Chen et al., 2019a; Chen et al., 2019b). BML for investigating effects is particularly well suited for neuroimaging studies as we know brain regions are not independent, i.e., one ROI could be involved in multiple pairwise connections. In addition, unlike the NHST assumption that the likelihood of an effect is equally anywhere between negative infinity and positive infinity for every brain region, BML assumes the likelihood of an effect more likely follows a Gaussian distribution (Chen et al., 2021). In NHST models this is corrected for via multiple comparisons correction but the BML approach not only negates the need for multiple comparisons correction, it also captures the covariant structure of the brain, allowing for a more robust and efficient inference.
To investigate potential FC group differences, we ran the matrix-based analysis (MBA) program which uses BML to pool and share variability estimates of ROIs in each subject’s Fisher Z-transformed correlation matrix. The contribution of each ROI in the matrix is captured in the model, showing its overall contribution to connections in the matrix. MBA offers effect estimations of both region and region pairs or connections (Chen et al., 2019a). The reporting of ROI effects is a unique feature of MBA, outputting information about regional contributions. Region effects characterize regional contributions within a multilevel framework, complementing paired FC. We included sex as a confounding variable in all our models because studies find biological sex is related to different patterns of functional activation, which may be confounded by pubertal hormones during adolescence (Gur and Gur, 2016; Lenroot and Giedd, 2010). Our final MBA model was of the form:
where ROI1:ROI2 represent a connection. As such the resulting effect estimations were at the ROI level, i.e., 80 ROIs and region pair level, i.e., 3,160 region pairs.
To investigate functional network organization between CHUU and CPHIV, we ran the region-based analysis (RBA) program for each graph theory measure (Chen et al., 2019b). RBA uses BML to pool and share information among brain regions or nodes. Here we also included sex as a confounding variable and our final RBA model was of the form:
The resulting effect estimates were for every node.
To assess the quality of MCMC we used the split Ř statistic which gauges the consistency of an ensemble of Markov chains. Fully converged chains correspond to Ř = 1. In the present study we used an arbitrary threshold of Ř < 1.1 for chains that converged well in both the MBA and RBA analyses.
Associations between functional outcomes and auditory/language-related scores
We examined the association between functional outcomes (FC between ROIs/graph measure outcomes for ROIs) and KABC-II auditory/language-related neurodevelopmental scores. Associations were investigated among connections or ROI pairs that included an auditory ROI, i.e., 604 connections. Graph measure associations were investigated with all 5 graph measures (i.e., node degree, strength, transitivity, nodal and local efficiency) for each of the 80 nodes.
To investigate the associations we ran Pearson correlation statistical tests. Where a variable was not normally distributed (Shapiro–Wilk p < 0.05), we investigated the relationship via Spearman correlation statistical test as it is resistant to outliers. We included sex as a confounding variable and investigated handedness as potential confounder due to evidence of functional asymmetry in language processing areas that develops with age (Olulade et al., 2020; Vassal et al., 2016). We corrected for multiple comparisons using the False Discovery Rate (FDR) method and considered FDR q < 0.05 as our significance level.
Associations between functional outcomes and immune health markers
We ran analyses to investigate the association between functional outcomes and immune health markers within the CPHIV group. Correlation tests were either Pearson or Spearman depending on normality. We investigated early immune health (CD4% at study enrolment), immune health of participant at the time of 11 year scan and treatment related health markers (age at cART initiation and first VL suppression). We corrected for multiple comparisons using the FDR method and considered FDR q < 0.05 as our significance level.
Results
Out of 84 participants who underwent MRI scanning, a total of 6 participants (2 CHUU; 4 CPHIV) were excluded due to poor image quality and/or not meeting motion inclusion criteria. We present results from a sample of 78 children (23 CHUU; 55 CPHIV), demographics of which are presented in Table 1.

Table 1. Sample characteristics of children perinatally infected and living with (CPHIV) and uninfected controls (CHUU) (N = 78).
We grouped our results based on strength of evidence as follows: very strong evidence of an effect – P + > = 0.975 or P + < = 0.025, strong evidence of an effect – P + > = 0.95 or P + < = 0.05 and evidence of an effect – P + > = 0.90 or P + < = 0.10. Our discussion mostly focuses on the connections that showed strong evidence of an effect.
Functional connectivity results
The Ř values were less than 1.1 for the CHUU vs. CPHIV FC analysis indicating all MCMC chains converged well. For a summary of model results please refer to Supplementary Appendix B.
ROI effects
Among the 80 ROIs, 14 ROIs showed strong evidence of lower FC in the CPHIV group. Regions include the left PAC, bilateral hippocampus, bilateral hippocampus CA1, bilateral hippocampus CA3, bilateral hippocampus CA4, left hippocampus molecular layer, left hippocampus presubiculum, right hippocampus GCDG, and bilateral lingual gyrus. In addition, 10 ROIs showed strong evidence of higher FC (P + < = 0.05) in the CPHIV group including the left mediodorsal thalamic nucleus, right anterior thalamic nucleus, left caudate, left postcentral gyrus, right precentral gyrus, right rMFC, right parsopecularis, right pallidum and bilateral putamen (Figure 3). Supplementary Appendix C shows the ROI effect estimates including uncertainty intervals and posterior probabilities of the effect being positive.

Figure 3. We found 73 connections showing evidence (P+ > = 0.90 | P + <=0.10) of altered FC between the CPHIV and CHUU groups (right). 43 connections showed lower FC (blue links) while 30 connections showed higher FC (red links) in the CPHIV group compared to CHUU. We found 29 connections showing strong evidence (P+ > =0.95 | P + <=0.05) (middle) and 6 connections showing very strong evidence (P+ > =0.975 | P + <=0.025) of altered FC (left) between groups. ROIs are colored according to direction of ROI effect with red and blue ROIs indicating higher and lower ROI FC in the CPHIV compared to CHUU, respectively. ROIs are sized according to the magnitude of ROI effect. Lower FC in the CPHIV group was predominantly from the left PAC and right CN/SOC to hippocampal subregions. Higher FC was predominantly from the bilateral MGN, left IC, right PAC to basal ganglia regions (putamen, pallidum and caudate) as well as to frontal lobe regions. (L/R)_parstr – parstriangularis, (L/R)_ pars – parsorbitalis, (L/R)_parso – parsopercularis, (L/R)_rMFC – rostral middle-frontal cortex, (L/R)_SFC – superior frontal cortex, (L/R)_cMFC – caudal middle-frontal cortex, (L/R)_precentral - precentral gyrus, (L/R)_paracentral – paracentral gyrus, (L/R)_rACC - rostral anterior cingulate cortex, (L/R)_cACC – caudal anterior cingulate cortex, (L/R)_PCC – posterior cingulate cortex, (L/R)_Isth – isthmus, (L/R)_postcentral – postcentral gyrus, (L/R)_SPC – superior parietal cortex, (L/R)_IPC – inferior parietal cortex, (L/R)_pericalc – pericalcarine, (L/R)_LOC – lateral occipital cortex, (L/R)_lingual – lingual gyrus, (L/R)_bankssts – banks of superior temporal sulcus, (L/R)_PAC – primary auditory cortex, (L/R)_Anterior – anterior thalamic nucleus, (L/R)_Medio_Dorsal – medio-dorsal thalamic nucleus, (L/R)_MGN – medial geniculate nucleus, (L/R)_VLD – ventral latero-dorsal thalamic nucleus, (L/R)_CLLPM_Pulvinar – central latero-lateral posterior medial pulvinar, (L/R)_V_Anterior – ventral anterior thalamic nucleus, (L/R)_Hippo – hippocampus, (L/R)_Hippo_Pres – hippocampus presubiculum, (L/R)_Hippo_CA1,3,4 – hippocampus cornu ammonis 1,3,4, (L/R)_Hippo_Mol_layer_HP – hippocampus molecular layer, (L/R)_Hippo_Tail – hippocampus tail, (L/R)_IC – inferior colliculus, (L/R)_SOC_CN – cochlear nucleus/superior olivary complex. L/R – left/right.
ROI pair effects
Compared to CHUU, the CPHIV showed strong evidence of lower FC between the left PAC and bilateral lingual gyrus, bilateral hippocampus CA1, bilateral hippocampus CA3, left hippocampus, left hippocampus presubiculum, and the left hippocampus molecular layer. The right CN/SOC also showed lower FC to the same regions as the left PAC in the CPHIV group. Further, we found strong evidence of higher FC (P + < = 0.05) in the CPHIV group between the right PAC and right putamen, bilateral MGN and right putamen, right pallidum, right parsopercularis, bilateral IC to the right putamen, and left CN/SOC to the right putamen (Figure 3). ROI pair effect estimates for all connections to auditory regions can be found in Supplementary Appendix C.
Functional network results
RBA showed strong evidence (P + > = 0.95) of lower nodal degree of the left IC in the CPHIV group compared to CHUU. The right IC showed evidence of lower nodal degree in the same group (P + = 0.94). Additionally, we found strong evidence (P+ > = 0.95) of lower nodal local efficiency of the right CN/SOC, left pulvinar and right IC in the CPHIV group compared to CHUU. The left IC showed evidence of lower nodal local efficiency in the same group (P + = 0.91) (Figure 4). The nodal effect estimates for all graph measures can be found here Supplementary Appendix D.

Figure 4. Posterior probability distribution plots for the effect magnitude of the comparison CHUU-CPHIV difference between mean graph measures. The graph measures include nodal degree, strength and transitivity, and both nodal and local efficiency. The green vertical line marks the location of no effect. (L/R)_parstr – parstriangularis, (L/R)_ pars – parsorbitalis, (L/R)_parso – parsopercularis, (L/R)_rMFC – rostral middle-frontal cortex, (L/R)_SFC – superior frontal cortex, (L/R)_cMFC – caudal middle-frontal cortex, (L/R)_precentral – precentral gyrus, (L/R)_paracentral – paracentral gyrus, (L/R)_rACC – rostral anterior cingulate cortex, (L/R)_cACC – caudal anterior cingulate cortex, (L/R)_PCC – posterior cingulate cortex, (L/R)_Isth – isthmus, (L/R)_postcentral – postcentral gyrus, (L/R)_SPC – superior parietal cortex, (L/R)_IPC – inferior parietal cortex, (L/R)_pericalc – pericalcarine, (L/R)_LOC – lateral occipital cortex, (L/R)_lingual – lingual gyrus, (L/R)_bankssts – banks of superior temporal sulcus, (L/R)_PAC – primary auditory cortex, (L/R)_Anterior – anterior thalamic nucleus, (L/R)_Medio_Dorsal – medio-dorsal thalamic nucleus, (L/R)_MGN – medial geniculate nucleus, (L/R)_VLD – ventral latero-dorsal thalamic nucleus, (L/R)_CLLPM_Pulvinar – central latero-lateral posterior medial pulvinar, (L/R)_V_Anterior – ventral anterior thalamic nucleus, (L/R)_Hippo – hippocampus, (L/R)_Hippo_Pres – hippocampus presubiculum, (L/R)_Hippo_CA1,3,4 – hippocampus cornu ammonis 1,3,4, (L/R)_Hippo_Mol_layer_HP – hippocampus molecular layer, (L/R)_Hippo_Tail – hippocampus tail, (L/R)_IC – inferior colliculus, (L/R)_SOC_CN – cochlear nucleus/superior olivary complex. L/R – left/right.
Associations between functional outcomes and auditory/language-related scores
Among the sample as a whole, we found no associations between FC and the auditory/language-related scores after correcting for multiple comparisons. However, among CHUU, we found negative associations between standardized Delayed Recall scores and nodal efficiency in the right middle temporal cortex (MTC), as well as a significantly positive association between standardized Delayed Recall and nodal efficiency in the left superior frontal cortex (SFC) (FDR q < 0.05) (Figure 5). These associations were not seen in CPHIV.

Figure 5. Associations of nodal efficiency in the left superior frontal cortex (SFC) and right middle temporal cortex (MTC) with delayed recall test scores assessed on the Kaufman assessment battery for children – second edition (KABC-II). Std = standardized. L – left, R –right.
Associations between functional outcomes and immune health markers
We found no significant associations between functional outcomes and immune health makers after multiple comparison correction.
Discussion
The current study investigated the functional integrity of the CAS in relation to HIV in an adolescent cohort. CPHIV demonstrated reduced FC to the PAC, from hippocampal areas and the lingual gyrus which are involved in limbic and visual functions, respectively. In addition, we reported evidence of increased FC to the bilateral MGN, with connections to the parsopercularis, an important region in speech production, most affected. Our analysis shows HIV infection influences the relationships between auditory structures and non-auditory systems, including the auditory–visual and auditory-limbic, in 11-year old children.
Auditory-limbic connectivity
Both ROI and region pair analysis identified lower FC related to HIV infection in the left PAC, bilateral lingual gyrus and bilateral hippocampal regions. Altered FC of the left PAC is noteworthy as it is responsible for language processing, which displays hemisphere lateralisation (Bernal and Altman, 2012; Binder et al., 1996; Olulade et al., 2020). In fact, language lateralisation increases significantly between childhood and adolescence (Szaflarski et al., 2006). The majority of our cohort was right-handed, suggesting left hemisphere dominance for language processing. However, even though alterations to the left PAC suggest direct effects on language processing we did not find significant associations between FC to the region and language related neurodevelopmental outcomes. Instead, we found evidence of reduced FC between the left PAC and non-auditory regions including the lingual gyri and hippocampal regions.
The link between the hippocampus, a limbic system region located in the medial part of the temporal lobe, and the auditory system has been previously reported (Kraus and Canlon, 2012; Zhang et al., 2022). The hippocampus is involved in the formation of long-term auditory memories (Nyhus and Curran, 2010; Watanabe et al., 2008) as well as the comprehension of spoken language (Davis and Johnsrude, 2003; van de Ven et al., 2020). Cortico-limbic FC is particularly vulnerable in the adolescent brain as cognition-emotion interactions are being refined (Simmonds et al., 2014). Lower FC between hippocampal regions and PAC point to possible interruptions of these relationships in the presence of HIV. Further, hippocampal neuropathology has been linked to hearing impairment (Park et al., 2018; Uchida et al., 2018) and associated with poor memory (Squire et al., 2001). Notably, CPHIV in this cohort at age 7 were found to have lower hippocampal volumes (Nwosu et al., 2018).
Auditory–visual connectivity
On the other hand, the lingual gyrus is part of the visual network and involved in visual memory and visual language processing (Palejwala et al., 2021). We have previously reported reduced FC to the lingual gyrus in the same cohort at age 7 (Toich et al., 2018). The visual system relies heavily on auditory input for efficient and precise perception of the environment around us (Gori et al., 2012; Yu et al., 2018). Additionally, adult-like visual–auditory integration only reaches its peak in late-childhood making it at risk to disruptions in young children (Gori et al., 2012). Our results point to a vulnerability of the auditory and visual memory system in the presence of HIV.
The right CN/SOC did not show evidence of altered region based FC, however its connections to hippocampal regions and the bilateral lingual gyrus demonstrated lower FC in the CPHIV group compared to CHUU. The CN and SOC are crucial auditory processing structures. The CN is the initial site of subcortical auditory processing (Hackett, 2015), where auditory signals are enhanced and sent to the SOC. In the SOC, the signals are filtered and a crossing of signals from the contralateral SOC takes place (Peterson et al., 2021). Even at this level, the CAS is interacting with other systems for cognitive processing of attentional and motor actions (Stein and Stanford, 2008) to facilitate one’s behavioral or emotional response to auditory stimulation. These results echo alterations in auditory-limbic and auditory–visual FC. In addition, graph theory analysis revealed lower nodal degree and nodal local efficiency in the CN/SOC. The local efficiency of a node is the efficiency of the graph formed by the node’s neighbors (direct edges) in its absence. Unlike nodal efficiency, which reports the integration of a node within the entire network, local efficiency measures how efficiently the neighbors can communicate in the event of damage to the node (Fornito et al., 2016). This points to reduced integration among the immediate neighbors of both the CN and SOC in the presence of HIV.
Auditory-motor connectivity
We observed several connections between CAS regions and non-auditory regions showing evidence of higher FC in the CPHIV group compared to CHUU. The non-auditory regions were distributed among basal ganglia, thalamic nuclei and frontal areas. The basal ganglia are of particular interest as they show alterations in CPHIV in childhood across imaging modalities (Hoare et al., 2014; Mbugua et al., 2016; Pajkrt et al., 2016; Robertson et al., 2018). Within this cohort, we have reported HIV-related alterations and associations in the basal ganglia at earlier ages. At ages 5 and 7, putamen volumes were larger and smaller in the CPHIV group, respectively, compared to controls (Nwosu et al., 2018; Randall et al., 2017). At age 7 years, we reported a negative association between FC and CD4 count in infancy among CPHIV pointing to higher FC linked with a weaker immune system (Toich et al., 2018). A longitudinal study from ages 5 to 11 years of neurometabolism in our cohort reported elevated choline levels in the basal ganglia, suggesting persistent inflammation (van Biljon et al., 2021). Our current results at age 11 suggest overcompensation from all levels of the CAS to the basal ganglia, which may be related to the volumetric and neurometabolic differences observed at earlier ages. Basal ganglia are broadly involved in motor control and the putamen interacts with temporal regions to initiate phonological processes (Lanciego et al., 2012).
Auditory-language connectivity
The bilateral MGN also showed higher FC to the parsopercularis in the presence of HIV. The parsopercularis, a region part of the inferior frontal gyrus, plays a role in language and phonological processing. Interestingly altered FC in the inferior frontal gyrus has been reported previously in the same cohort, although FC was lower in CPHIV in the analysis at age 7 compared to CHUU (Toich et al., 2018). Altogether our results suggest increased interaction between the CAS and regions in other brain systems including limbic-motor, visual and language areas in the presence of HIV. Interestingly children with hearing impairments demonstrate a similar pattern of overcompensation in the auditory network, motor-limbic and visual areas (Luan et al., 2019; Shi et al., 2016; Wang et al., 2021).
Besides the reduced nodal integrity of the CN/SOC, our functional network analysis also revealed lower nodal degree and local efficiency in the IC of the CPHIV group as compared to CHUU. This result is unexpected as the IC did not show any evidence of FC alterations between groups. The reduced nodal integrity observed in the current analysis may be related to the lower degree observed in this node as well as the previously discussed CN/SOC. The nodal degree reflects the density of connections to the node, with a low degree reflecting nodes that are not as highly connected with other regions. It is worth noting that both these regions are in the brainstem where the measured BOLD signal is not high. Indeed, they were the structures with the lowest number of connections in Figure 1B.
Additionally, we found a positive association between nodal efficiency of the SFC and the Delayed Recall test as well as a negative association between nodal efficiency of the MTG and Delayed Recall. Both associations were seen in the CHUU group but not in CPHIV. Delayed Recall is part of the KABC Learning subscale (Kaufman and Kaufman, 2014). More specifically Delayed Recall is a composite score of Atlantis delayed, where children are asked to point to pictures they would have learnt 20 min earlier, and Rebus delayed where the children are taught pictures that have a word that they have to put in a sentence and recall 20 min later. Though both the tests have visual cues they assess the child’s language processing ability in coming up with a sentence as well as their ability to store and retrieve memories. The SFC is recruited during tasks of delayed recognition (Yee et al., 2010) while the MTG is involved in phonologic and sematic language processing (Middlebrooks et al., 2016), both processes of which interact with the CAS to facilitate speech production for communicating (Maffei et al., 2015). In addition temporo-frontal regions show altered FC in CPHIV (Toich et al., 2018; Yadav et al., 2018). Our current results point to possible disruption in the relationships between temporo-frontal nodal integrity and neurodevelopment relating to working memory in CPHIV.
Limitations and future work
The origin and interpretation of negative FC has been an area of debate (Chen et al., 2011; Goelman et al., 2014; Mohanty et al., 2020) and the role of negative correlations as edge weights is not yet fully understood. Nodes with strong positive correlations suggest that the nodes work together whereas strong negative correlations imply the nodes have an antagonistic relationship (Fornito et al., 2016). In the present study we used only positive weights (Fisher Z-value >0.1) as the measure of pairwise connectivity in the functional network. This not only removes weak edge weights that represent spurious connections and may alter network topology, it also simplifies the interpretation of FC to cooperation between nodes. However, this approach discards information about other functional interactions in the network, particularly the very high but negative correlations. Negative correlations may imply that the pairs of regions work together in an anti-correlated way (Keller et al., 2015).
In addition, the choice of a threshold of Fisher Z-value > 0.1 was arbitrary. The threshold choice was used to improve the number of pairwise connections between subcortical CAS regions deep in the brainstem and non-auditory ROIs. As seen in Figure 1B of the functional network matrix, most of the FC to the CAS was to the PAC, then to the MGN, less to the IC and even less to the CN/SOC in the brainstem. Using a range of adjacency matrices during network definition, for example Fisher Z-value > 0.1 with 0.05 increments, as opposed to the current network definition of Fisher Z-value > 0.1 may help observe if the current results persist.
While this work examines the functional auditory system and its associations with language measures, it does not include measures of hearing. Future work is needed to link the identified alterations with hearing-related outcomes. This study examined associations with ART initiation, however it did not look at the potential influence of early ART interruptions or ART initiation group differences within the first 18 months of life. Future work may explore the question of variation in early ART initiation and HIV-related group auditory differences in more statistical depth. In addition, given the dynamic nature of the brain in adolescence, longitudinal analysis is needed to examine the long-term impact of the results presented.
Conclusion
To our knowledge this is the first study to investigate the functional integrity of the CAS and its interactions with other systems in CPHIV. Our results show altered FC between the CAS and non-auditory regions involved in limbic, visual and motor processes in CPHIV. Furthermore, CPHIV show altered nodal topology of the auditory brainstem and midbrain, with nodal measures related to cognitive scores in the CHUU group only. These results provide insight on FC in CPHIV and point to risk of auditory and language disruptions based on CAS alterations.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The studies involving humans were approved by the Human Research Ethics Committees of the Universities of Cape Town and Stellenbosch. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
JM: Formal analysis, Methodology, Visualization, Writing – original draft. MJ: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing. EM: Conceptualization, Funding acquisition, Resources, Writing – review & editing. PT: Conceptualization, Funding acquisition, Resources, Writing – review & editing. BL: Conceptualization, Funding acquisition, Writing – review & editing. MH: Conceptualization, Funding acquisition, Resources, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Institute On Deafness and Other Communication Disorders of the National Institutes of Health under Award Number R01DC015984; South African Medical Research Council (SAMRC); South African National Research Foundation (NRF) grants CPR20110614000019421 and CPRR150723129691; and the NRF/DST South African Research Chairs Initiative; Support for the CHER study, which provided the infrastructure for the neurodevelopmental study, was provided by the US National Institute of Allergy and Infectious Diseases through the CIPRA network, Grant U19 AI53217; the Departments of Health of the Western Cape and Gauteng, South Africa; and GlaxoSmithKline/Viiv Healthcare. Additional support was provided with Federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, United States Department of Health and Human Services, under Contract No. HHSN272200800014C.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsys.2025.1508516/full#supplementary-material
References
Bernal, B., and Altman, N. R. (2012). Auditory functional MR Imaging. Amer. J. Roentgenol. 176, 1009–1015. doi: 10.2214/AJR.176.4.1761009
Binder, J. R., Swanson, S. J., Hammeke, T. A., Morris, G. L., Mueller, W. M., Fischer, M., et al. (1996). Determination of language dominance using functional MRI. Neurology 46, 978–984. doi: 10.1212/WNL.46.4.978
Bullmore, E. T., Fornito, A., and Zalesky, A. (2016). “Fundamentals of brain network analysis” in Fundamentals of brain network analysis.
Chao, C. K., Czechowicz, J. A., Messner, A. H., Alarcoń, J., Roca, L. K., Larragán Rodriguez, M. M., et al. (2012). High prevalence of hearing impairment in HIV-infected Peruvian children. Otolaryngol. Head Neck Surg. 146, 259–265. doi: 10.1177/0194599811429271
Chen, G., Bürkner, P. C., Taylor, P. A., Li, Z., Yin, L., Glen, D. R., et al. (2019a). An integrative Bayesian approach to matrix-based analysis in neuroimaging. Hum. Brain Mapp. 40, 4072–4090. doi: 10.1002/HBM.24686
Chen, G., Chen, G., Xie, C., and Li, S. J. (2011). Negative functional connectivity and its dependence on the shortest path length of positive network in the resting-state human brain. Brain Connect. 1, 195–206. doi: 10.1089/BRAIN.2011.0025
Chen, G., Taylor, P. A., Stoddard, J., Cox, R. W., Bandettini, P. A., and Pessoa, L. (2021). Sources of information waste in neuroimaging: mishandling structures, thinking dichotomously, and over-reducing data. BioRxiv 2021:3246. doi: 10.1101/2021.05.09.443246
Chen, G., Xiao, Y., Taylor, P. A., Rajendra, J. K., Riggins, T., Geng, F., et al. (2019b). Handling multiplicity in neuroimaging through Bayesian lenses with multilevel modeling. Neuroinformatics 17, 515–545. doi: 10.1007/s12021-018-9409-6
Cotton, M. F., Violari, A., Otwombe, K., Panchia, R., Dobbels, E., Rabie, H., et al. (2014). Early limited antiretroviral therapy is superior to deferred therapy in HIV-infected south African infants: Results from the CHER (children with HIV early antiRetroviral) randomized trial mark. Lancet 382, 1555–1563. doi: 10.1016/S0140-6736(13)61409-9.Early
Daducci, A., Gerhard, S., Griffa, A., Lemkaddem, A., Cammoun, L., Gigandet, X., et al. (2012). The connectome mapper: an open-source processing pipeline to map connectomes with MRI. PLoS One 7:e48121. doi: 10.1371/JOURNAL.PONE.0048121
Davis, M. H., and Johnsrude, I. S. (2003). Hierarchical processing in spoken language comprehension. J. Neurosci. 23, 3423–3431. doi: 10.1523/JNEUROSCI.23-08-03423.2003
Fitzroy, A. B., Krizman, J., Tierney, A., Agouridou, M., and Kraus, N. (2015). Longitudinal maturation of auditory cortical function during adolescence. Front. Hum. Neurosci. 9. doi: 10.3389/fnhum.2015.00530
Goelman, G., Gordon, N., and Bonne, O. (2014). Maximizing negative correlations in resting-state functional connectivity MRI by time-lag. PLoS One 9:e111554. doi: 10.1371/JOURNAL.PONE.0111554
Gori, M., Sandini, G., and Burr, D. (2012). Development of visuo-auditory integration in space and time. Front. Integrative Neurosci. 6:77. doi: 10.3389/fnint.2012.00077
Gur, R. E., and Gur, R. C. (2016). Sex differences in brain and behavior in adolescence: findings from the Philadelphia neurodevelopmental cohort. Neurosci. Biobehav. Rev. 70, 159–170. doi: 10.1016/J.NEUBIOREV.2016.07.035
Hackett, T. A. (2015). Anatomic organization of the auditory cortex. Handb. Clin. Neurol. 129, 27–53. doi: 10.1016/B978-0-444-62630-1.00002-0
Herting, M. M., Uban, K. A., Williams, P. L., Gautam, P., Huo, Y., Malee, K., et al. (2015). Default mode connectivity in youth with perinatally acquired HIV. Medicine 94:e1417. doi: 10.1097/MD.0000000000001417
Hoare, J., Ransford, G. L., Phillips, N., Amos, T., Donald, K., and Stein, D. J. (2014). Systematic review of neuroimaging studies in vertically transmitted HIV positive children and adolescents. Metab. Brain Dis. 29, 221–229. doi: 10.1007/s11011-013-9456-5
Hoff, G. E. A. J., Van den Heuvel, M. P., Benders, M. J. N. L., Kersbergen, K. J., and De Vries, L. S. (2013). On development of functional brain connectivity in the young brain. Front. Hum. Neurosci. 7, 1–7. doi: 10.3389/fnhum.2013.00650
Kaufman, A. S., and Kaufman, N. L. (2014). “Kaufman assessment battery for children, second edition” in Encyclopedia of special education (John Wiley & Sons, Inc). doi: 10.1002/9781118660584.ese1324
Keller, J. B., Hedden, T., Thompson, T. W., Anteraper, S. A., Gabrieli, J. D. E., and Whitfield-Gabrieli, S. (2015). Resting-state anticorrelations between medial and lateral prefrontal cortex: association with working memory, aging, and individual differences. Cortex 64, 271–280. doi: 10.1016/J.CORTEX.2014.12.001
Kraus, K. S., and Canlon, B. (2012). Neuronal connectivity and interactions between the auditory and limbic systems. Effects of noise and tinnitus. Hear. Res. 288, 34–46. doi: 10.1016/J.HEARES.2012.02.009
Kruschke, J. K., and Liddell, T. M. (2018). The Bayesian new statistics: hypothesis testing, estimation, meta-analysis, and power analysis from a Bayesian perspective. Psychon. Bull. Rev. 25, 178–206. doi: 10.3758/S13423-016-1221-4/FIGURES/12
Lanciego, J. L., Luquin, N., and Obeso, J. A. (2012). Functional neuroanatomy of the basal ganglia. Cold Spring Harb. Perspect. Med. 2. doi: 10.1101/cshperspect.a009621
Lenroot, R. K., and Giedd, J. N. (2010). Sex differences in the adolescent brain. Brain Cogn. 72, 46–55. doi: 10.1016/j.bandc.2009.10.008
Lima, J. V. D. S., Morais, C. F. M. de, Zamberlan-Amorim, N. E., Mandrá, P. P., and Reis, A. C. M. B. (2023). Neurocognitive function in children with cochlear implants and hearing aids: a systematic review. Front. Neurosci. 17:1242949. doi: 10.3389/FNINS.2023.1242949/BIBTEX
Luan, Y., Wang, C., Jiao, Y., Tang, T., Zhang, J., and Teng, G. J. (2019). Dysconnectivity of multiple resting-state networks associated with higher-order functions in sensorineural hearing loss. Front. Neurosci. 13. doi: 10.3389/fnins.2019.00055
Madhi, S. A., Adrian, P., Cotton, M. F., McIntyre, J. A., Jean-Philippe, P., Meadows, S., et al. (2010). Effect of HIV infection status and anti-retroviral treatment on quantitative and qualitative antibody responses to pneumococcal conjugate vaccine in infants. J. Infect. Dis. 202, 355–361. doi: 10.1086/653704
Maffei, C., Soria, G., Prats-Galino, A., and Catani, M. (2015). Imaging white-matter pathways of the auditory system with diffusion imaging tractography.
Maro, I., Fellows, A., Clavier, O., Gui, J., Rieke, C., Wilbur, J., et al. (2017). Auditory impairments in HIV-infected children. Physiol. Behav. 176, 139–148. doi: 10.1016/j.physbeh.2017.03.040
Mbugua, K. K., Holmes, M. J., Cotton, M. F., Ratai, E.-M., Little, F., Hess, A. T., et al. (2016). HIV-associated CD4/8 depletion in infancy is associated with neurometabolic reductions in the basal ganglia at age 5 years despite early antiretroviral therapy HHS public access. AIDS 30, 1353–1362. doi: 10.1097/QAD.0000000000001082
Middlebrooks, E. H., Yagmurlu, K., Szaflarski, J. P., Rahman, M., and Bozkurt, B. (2016). A contemporary framework of language processing in the human brain in the context of preoperative and intraoperative language mapping. Neuroradiology 59, 69–87. doi: 10.1007/S00234-016-1772-0
Mohanty, R., Sethares, W. A., Nair, V. A., and Prabhakaran, V. (2020). Rethinking measures of functional connectivity via feature extraction. Sci. Rep. 10:1298. doi: 10.1038/s41598-020-57915-w
Morgan, E. E., Iudicello, J. E., Cattie, J. E., Blackstone, K., Grant, I., and Woods, S. P. (2015). Neurocognitive impairment is associated with lower health literacy among persons living with HIV infection. AIDS Behav. 19, 166–177. doi: 10.1007/s10461-014-0851-7
Nwosu, E. C., Robertson, F. C., Holmes, M. J., Cotton, M. F., Dobbels, E., Little, F., et al. (2018). Altered brain morphometry in 7-year old HIV-infected children on early ART. Metab. Brain Dis. 33, 523–535. doi: 10.1007/s11011-017-0162-6
Nyhus, E., and Curran, T. (2010). Functional role of gamma and theta oscillations in episodic memory. Neurosci. Biobehav. Rev. 34, 1023–1035. doi: 10.1016/J.NEUBIOREV.2009.12.014
Olulade, O. A., Seydell-Greenwald, A., Chambers, C. E., Turkeltaub, P. E., Dromerick, A. W., Berl, M. M., et al. (2020). The neural basis of language development: changes in lateralization over age. Proc. Natl. Acad. Sci. USA 117, 23477–23483. doi: 10.1073/pnas.1905590117
Pajkrt, D., Blokhuis, C., Caan, M., and Kootstra, N. (2016). Neurodevelopmental delay in pediatric HIV/AIDS: current perspectives. Neurobehavioral HIV Medicine 7:1. doi: 10.2147/NBHIV.S68954
Palejwala, A. H., Dadario, N. B., Young, I. M., O’Connor, K., Briggs, R. G., Conner, A. K., et al. (2021). Anatomy and white matter connections of the lingual gyrus and cuneus. World Neurosurg. 151, e426–e437. doi: 10.1016/J.WNEU.2021.04.050
Park, S. Y., Kim, M. J., Kim, H. L., Kim, D. K., Yeo, S. W., and Park, S. N. (2018). Cognitive decline and increased hippocampal p-tau expression in mice with hearing loss. Behav. Brain Res. 342, 19–26. doi: 10.1016/J.BBR.2018.01.003
Peterson, D. C., Reddy, V., and Hamel, R. N. (2021). Neuroanatomy. StatPearls: Auditory Pathway. Available online at: https://www.ncbi.nlm.nih.gov/books/NBK532311/
Randall, S. R., Warton, C. M. R., Holmes, M. J., Cotton, M. F., Laughton, B., van der Kouwe, A. J. W., et al. (2017). Larger subcortical gray matter structures and smaller corpora Callosa at age 5 years in HIV infected children on early ART. Front. Neuroanat. 11:95. doi: 10.3389/fnana.2017.00095
Rice, M. L., Buchanan, A. L., Siberry, G. K., Malee, K. M., Zeldow, B., Frederick, T., et al. (2012). Language impairment in children perinatally infected with hiv compared to children who were hiv-exposed and uninfected. DBP. 33, 112–123. doi: 10.1097/DBP.0b013e318241ed23
Robertson, F. C., Holmes, M. J., Cotton, M. F., Dobbels, E., Little, F., Laughton, B., et al. (2018). Perinatal HIV infection or exposure is associated with low n-acetylaspartate and glutamate in basal ganglia at age 9 but not 7 years. Front. Hum. Neurosci. 12:145. doi: 10.3389/fnhum.2018.00145
Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. NeuroImage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003
Shi, B., Yang, L.-Z., Liu, Y., Zhao, S.-L., Wang, Y., Gu, F., et al. (2016). Early-onset hearing loss reorganizes the visual and auditory network in children without cochlear implantation. Neuroreport 27, 197–202. doi: 10.1097/WNR.0000000000000524
Simmonds, D. J., Hallquist, M. N., Asato, M., and Luna, B. (2014). Developmental stages and sex differences of white matter and behavioral development through adolescence: a longitudinal diffusion tensor imaging (DTI) study. NeuroImage 92, 356–368. doi: 10.1016/J.NEUROIMAGE.2013.12.044
Skoe, E., Krizman, J., Anderson, S., and Kraus, N. (2015). Stability and plasticity of auditory brainstem function across the lifespan. Cerebral Cortex 25, 1415–1426. doi: 10.1093/CERCOR/BHT311
Sporns, O. (2018). Graph theory methods: applications in brain networks. Dialogues Clin. Neurosci. 20, 111–121. doi: 10.31887/DCNS.2018.20.2/OSPORNS
Squire, L. R., Schmolck, H., and Stark, S. M. (2001). Impaired Auditory Recognition Memory in Amnesic Patients with Medial Temporal Lobe Lesions. 8, 252–256. doi: 10.1101/lm.42001
Stein, B. E., and Stanford, T. R. (2008). Multisensory integration: current issues from the perspective of the single neuron [nature reviews neuroscience (2008) 9, (255–266)]. Nat. Rev. Neurosci. 9:406. doi: 10.1038/nrn2377
Strauss, E., Sherman, E. M. S., and Spreen, O. (2010). E. Strauss, E. M. S. Sherman, & O. Spreen, A Compendium of Neuropsychological Tests: Administration, Norms, and Commentary. Applied Neuropsychol. 14, 62–63. doi: 10.1080/09084280701280502
Szaflarski, J. P., Holland, S. K., Schmithorst, V. J., and Byars, A. W. (2006). fMRI study of language lateralization in children and adults. Hum. Brain Mapp. 27, 202–212. doi: 10.1002/hbm.20177
Taylor, P. A., and Saad, Z. S. (2013). FATCAT: (an efficient) functional and Tractographic connectivity analysis toolbox. Brain Connect. 3, 523–535. doi: 10.1089/brain.2013.0154
Toich, J. T. F., Taylor, P. A., Holmes, M. J., Gohel, S., Cotton, M. F., Dobbels, E., et al. (2018). Functional connectivity alterations between networks and associations with infant immune health within networks in HIV infected children on early treatment: a study at 7 years. Front. Hum. Neurosci. 11, 1–15. doi: 10.3389/fnhum.2017.00635
Torre, P., Cook, A., Elliott, H., Dawood, G., and Laughton, B. (2015). Hearing assessment data in HIV-infected and uninfected children of Cape Town, South Africa. AIDS Care Psychological Socio Medical Aspects of AIDS/HIV 27, 1037–1041. doi: 10.1080/09540121.2015.1021746
Torre, P., Zeldow, B., Hoffman, H. J., Buchanan, A., Siberry, G. K., Rice, M., et al. (2012). Hearing loss in perinatally HIV-infected and HIV-exposed but uninfected children and adolescents. Pediatr. Infect. Dis. J. 31, 835–841. doi: 10.1097/INF.0b013e31825b9524
Tourbier, S., Rue-Queralt, J., Glomb, K., Aleman-Gomez, Y., Mullier, E., Griffa, A., et al. (2022). Connectome mapper 3: a flexible and open-source pipeline software for multiscale multimodal human connectome mapping. J. Open Source Software 7:4248. doi: 10.21105/JOSS.04248
Uchida, Y., Nishita, Y., Kato, T., Iwata, K., Sugiura, S., Suzuki, H., et al. (2018). Smaller hippocampal volume and degraded peripheral hearing among Japanese community dwellers. Front. Aging Neurosci. 10:319. doi: 10.3389/FNAGI.2018.00319/FULL
van Biljon, N., Robertson, F., Holmes, M., Cotton, M. F., Laughton, B., van der Kouwe, A., et al. (2021). Multivariate approach for longitudinal analysis of brain metabolite levels from ages 5-11 years in children with perinatal HIV infection. NeuroImage 237:118101. doi: 10.1016/j.neuroimage.2021.118101
van de Ven, V., Waldorp, L., and Christoffels, I. (2020). Hippocampus plays a role in speech feedback processing. NeuroImage 223:117319. doi: 10.1016/J.NEUROIMAGE.2020.117319
van den Heuvel, M. P., and Hulshoff Pol, H. E. (2010). Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 20, 519–534. doi: 10.1016/j.euroneuro.2010.03.008
van den Heuvel, M. P., and Sporns, O. (2013). Network hubs in the human brain. Trends Cogn. Sci. 17, 683–696. doi: 10.1016/J.TICS.2013.09.012
Vassal, F., Schneider, F., Boutet, C., Jean, B., Sontheimer, A., and Lemaire, J.-J. (2016). Combined DTI Tractography and functional MRI study of the language connectome in healthy volunteers: extensive mapping of white matter fascicles and cortical activations. PLoS One 11:152614:e0152614. doi: 10.1371/journal.pone.0152614
Violari, A., Cotton, M. F., Gibb, D. M., Babiker, A. G., Steyn, J., Madhi, S. A., et al. (2008). Early antiretroviral therapy and mortality among HIV-infected infants. N. Engl. J. Med. 359, 2233–2244. doi: 10.1056/NEJMoa0800971
Wang, S., Chen, B., Yu, Y., Yang, H., Cui, W., Fan, G., et al. (2021). Altered resting-state functional network connectivity in profound sensorineural hearing loss infants within an early sensitive period: a group ICA study. Hum. Brain Mapp. 42, 4314–4326. doi: 10.1002/HBM.25548
Watanabe, T., Yagishita, S., and Kikyo, H. (2008). Memory of music: roles of right hippocampus and left inferior frontal gyrus. NeuroImage 39, 483–491. doi: 10.1016/J.NEUROIMAGE.2007.08.024
Yadav, S. K., Gupta, R. K., Hashem, S., Bhat, A. A., Garg, R. K., Venkatesh, V., et al. (2018). Changes in resting-state functional brain activity are associated with waning cognitive functions in HIV-infected children. NeuroImage 20, 1204–1210. doi: 10.1016/j.nicl.2018.10.028
Yee, L. T. S., Roe, K., and Courtney, S. M. (2010). Selective involvement of superior frontal cortex during working memory for shapes. J. Neurophysiol. 103, 557–563. doi: 10.1152/JN.91299.2008
Yu, Z., Li, L., Song, J., and Lv, H. (2018). The study of visual-auditory interactions on lower limb motor imagery. Front. Neurosci. 12:509. doi: 10.3389/fnins.2018.00509
Keywords: central auditory system, HIV, resting-state functional connectivity, functional brain network, graph theory, Bayesian analysis
Citation: Madzime J, Jankiewicz M, Meintjes EM, Torre P III, Laughton B and Holmes MJ (2025) Evidence of functional connectivity disruptions between auditory and non-auditory regions in adolescents living with HIV. Front. Syst. Neurosci. 19:1508516. doi: 10.3389/fnsys.2025.1508516
Edited by:
Paolo Taurisano, University of Bari Aldo Moro, ItalyReviewed by:
Charalampos D. Moschopoulos, University General Hospital Attikon, GreeceJoga Chaganti, Thomas Jefferson University, United States
Copyright © 2025 Madzime, Jankiewicz, Meintjes, Torre, Laughton and Holmes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joanah Madzime, anMubWFkemltZUBnbWFpbC5jb20=