A subset of brain regions within adult functional connectivity networks demonstrate high reliability across early development.

The human cerebral cortex contains groups of areas that support sensory, motor, cognitive, and affective functions, often categorized as functional networks. These areas show stronger internal and weaker external functional connectivity (FC) and exhibit similar FC profiles within rather than between networks. Previous studies have demonstrated the development of these networks from nascent forms present before birth to their mature, adult-like topography in childhood. However, analyses often still use definitions based on adult functional networks. We aim to assess how this might lead to the misidentification of functional networks and explore potential consequences and solutions. Our findings suggest that even though adult networks provide only a marginally better than-chance description of the infant FC organization, misidentification was largely driven by specific areas. By restricting functional networks to areas showing adult-like network clustering, we observed consistent within-network FC both within and across scans and throughout development. Additionally, these areas were spatially closer to locations with low variability in network identity among adults. Our analysis aids in understanding the potential consequences of using adult networks "as is" and provides guidance for future research on selecting and utilizing functional network models based on the research question and scenario. Highlights: Specialized functional networks in the human cerebral cortex, evident in resting-state fMRI, support sensory, motor, cognitive, and affective functions and evolve throughout the lifespan.Existing studies have focused on age-specific networks for infants, but less on to what extent adult networks can describe infant functional connectivity (FC).Analysis revealed a subset of areas in infants showing adult-like network organization, with within-network FC exhibiting less variation across age and higher reliability across scans.These areas are posited near locations with low variability in functional network identity in adults, suggestive of the relationship between developmental sequence and interindividual variability in functional network organization.

Familial risk for depression moderates neural circuitry in healthy preadolescents to predict adolescent depression symptoms in the Adolescent Brain Cognitive Development (ABCD) Study.

BACKGROUND: There is an imminent need to identify neural markers during preadolescence that are linked to developing depression during adolescence, especially among youth at elevated familial risk. However, longitudinal studies remain scarce and exhibit mixed findings. Here we aimed to elucidate functional connectivity (FC) patterns among preadolescents that interact with familial depression risk to predict depression two years later. METHODS: 9-10 year-olds in the Adolescent Brain Cognitive Development (ABCD) Study were classified as healthy (i.e., no lifetime psychiatric diagnoses) at high familial risk for depression (HR; n=559) or at low familial risk for psychopathology (LR; n=1203). Whole-brain seed-to-voxel resting-state FC patterns with the amygdala, putamen, nucleus accumbens, and caudate were calculated. Multi-level, mixed-effects regression analyses were conducted to test whether FC at ages 9-10 interacted with familial risk to predict depression symptoms at ages 11-12. RESULTS: HR youth demonstrated stronger associations between preadolescent FC and adolescent depression symptoms (ps<0.001) as compared to LR youth (ps>0.001), primarily among amygdala/striatal FC with visual and sensory/somatomotor networks. CONCLUSIONS: Preadolescent amygdala and striatal FC may be useful biomarkers of adolescent-onset depression, particularly for youth with family histories of depression. This research may point to neurobiologically-informed approaches to prevention and intervention for depression in adolescents.

A general exposome factor explains individual differences in functional brain network topography and cognition in youth.

Childhood environments are critical in shaping cognitive neurodevelopment. With the increasing availability of large-scale neuroimaging datasets with deep phenotyping of childhood environments, we can now build upon prior studies that have considered relationships between one or a handful of environmental and neuroimaging features at a time. Here, we characterize the combined effects of hundreds of inter-connected and co-occurring features of a child's environment ("exposome") and investigate associations with each child's unique, multidimensional pattern of functional brain network organization ("functional topography") and cognition. We apply data-driven computational models to measure the exposome and define personalized functional brain networks in pre-registered analyses. Across matched discovery (n=5139, 48.5% female) and replication (n=5137, 47.1% female) samples from the Adolescent Brain Cognitive Development study, the exposome was associated with current (ages 9-10) and future (ages 11-12) cognition. Changes in the exposome were also associated with changes in cognition after accounting for baseline scores. Cross-validated ridge regressions revealed that the exposome is reflected in functional topography and can predict performance across cognitive domains. Importantly, a single measure capturing a child's exposome could more accurately and parsimoniously predict cognition than a wealth of personalized neuroimaging data, highlighting the importance of children's complex, multidimensional environments in cognitive neurodevelopment.

A precision functional atlas of personalized network topography and probabilities.

Although the general location of functional neural networks is similar across individuals, there is vast person-to-person topographic variability. To capture this, we implemented precision brain mapping functional magnetic resonance imaging methods to establish an open-source, method-flexible set of precision functional network atlases-the Masonic Institute for the Developing Brain (MIDB) Precision Brain Atlas. This atlas is an evolving resource comprising 53,273 individual-specific network maps, from more than 9,900 individuals, across ages and cohorts, including the Adolescent Brain Cognitive Development study, the Developmental Human Connectome Project and others. We also generated probabilistic network maps across multiple ages and integration zones (using a new overlapping mapping technique, Overlapping MultiNetwork Imaging). Using regions of high network invariance improved the reproducibility of executive function statistical maps in brain-wide associations compared to group average-based parcellations. Finally, we provide a potential use case for probabilistic maps for targeted neuromodulation. The atlas is expandable to alternative datasets with an online interface encouraging the scientific community to explore and contribute to understanding the human brain function more precisely.

Perinatal Western-style diet exposure associated with decreased microglial counts throughout the arcuate nucleus of the hypothalamus in Japanese macaques.

Perinatal exposure to a high-fat, high-sugar Western-style diet (WSD) is associated with altered neural circuitry in the melanocortin system. This association may have an underlying inflammatory component, as consumption of a WSD during pregnancy can lead to an elevated inflammatory environment. Our group previously demonstrated that prenatal WSD exposure was associated with increased markers of inflammation in the placenta and fetal hypothalamus in Japanese macaques. In this follow-up study, we sought to determine whether this heightened inflammatory state persisted into the postnatal period, as prenatal exposure to inflammation has been shown to reprogram offspring immune function and long-term neuroinflammation would present a potential means for prolonged disruptions to microglia-mediated neuronal circuit formation. Neuroinflammation was approximated in 1-yr-old offspring by counting resident microglia and peripherally derived macrophages in the region of the hypothalamus examined in the fetal study, the arcuate nucleus (ARC). Microglia and macrophages were immunofluorescently stained with their shared marker, ionized calcium-binding adapter molecule 1 (Iba1), and quantified in 11 regions along the rostral-caudal axis of the ARC. A mixed-effects model revealed main effects of perinatal diet (P = 0.011) and spatial location (P = 0.003) on Iba1-stained cell count. Perinatal WSD exposure was associated with a slight decrease in the number of Iba1-stained cells, and cells were more densely located in the center of the ARC. These findings suggest that the heightened inflammatory state experienced in utero does not persist postnatally. This inflammatory response trajectory could have important implications for understanding how neurodevelopmental disorders progress.NEW & NOTEWORTHY Prenatal Western-style diet exposure is associated with increased microglial activity in utero. However, we found a potentially neuroprotective reduction in microglia count during early postnatal development. This trajectory could inform the timing of disruptions to microglia-mediated neuronal circuit formation. Additionally, this is the first study in juvenile macaques to characterize the distribution of microglia along the rostral-caudal axis of the arcuate nucleus of the hypothalamus. Nearby neuronal populations may be greater targets during inflammatory insults.

Cumulative Effects of Resting-State Connectivity Across All Brain Networks Significantly Correlate with Attention-Deficit Hyperactivity Disorder Symptoms.

Identification of replicable neuroimaging correlates of attention-deficit hyperactivity disorder (ADHD) has been hindered by small sample sizes, small effects, and heterogeneity of methods. Given evidence that ADHD is associated with alterations in widely distributed brain networks and the small effects of individual brain features, a whole-brain perspective focusing on cumulative effects is warranted. The use of large, multisite samples is crucial for improving reproducibility and clinical utility of brain-wide MRI association studies. To address this, a polyneuro risk score (PNRS) representing cumulative, brain-wide, ADHD-associated resting-state functional connectivity was constructed and validated using data from the Adolescent Brain Cognitive Development (ABCD, N = 5,543, 51.5% female) study, and was further tested in the independent Oregon-ADHD-1000 case-control cohort (N = 553, 37.4% female). The ADHD PNRS was significantly associated with ADHD symptoms in both cohorts after accounting for relevant covariates (p < 0.001). The most predictive PNRS involved all brain networks, though the strongest effects were concentrated among the default mode and cingulo-opercular networks. In the longitudinal Oregon-ADHD-1000, non-ADHD youth had significantly lower PNRS (Cohen's d = -0.318, robust p = 5.5 × 10-4) than those with persistent ADHD (age 7-19). The PNRS, however, did not mediate polygenic risk for ADHD. Brain-wide connectivity was robustly associated with ADHD symptoms in two independent cohorts, providing further evidence of widespread dysconnectivity in ADHD. Evaluation in enriched samples demonstrates the promise of the PNRS approach for improving reproducibility in neuroimaging studies and unraveling the complex relationships between brain connectivity and behavioral disorders.

Study design features that improve effect sizes in cross-sectional and longitudinal brain-wide association studies.

Brain-wide association studies (BWAS) are a fundamental tool in discovering brain-behavior associations. Several recent studies showed that thousands of study participants are required to improve the replicability of BWAS because actual effect sizes are much smaller than those reported in smaller studies. Here, we perform analyses and meta-analyses of a robust effect size index (RESI) using 63 longitudinal and cross-sectional magnetic resonance imaging studies from the Lifespan Brain Chart Consortium (77,695 total scans) to demonstrate that optimizing study design is critical for improving standardized effect sizes and replicability in BWAS. A meta-analysis of brain volume associations with age indicates that BWAS with larger covariate variance have larger effect size estimates and that the longitudinal studies we examined have systematically larger standardized effect sizes than cross-sectional studies. We propose a cross-sectional RESI to adjust for the systematic difference in effect sizes between cross-sectional and longitudinal studies that allows investigators to quantify the benefit of conducting their study longitudinally. Analyzing age effects on global and regional brain measures from the United Kingdom Biobank and the Alzheimer's Disease Neuroimaging Initiative, we show that modifying longitudinal study design through sampling schemes to increase between-subject variability and adding a single additional longitudinal measurement per subject can improve effect sizes. However, evaluating these longitudinal sampling schemes on cognitive, psychopathology, and demographic associations with structural and functional brain outcome measures in the Adolescent Brain and Cognitive Development dataset shows that commonly used longitudinal models can, counterintuitively, reduce effect sizes. We demonstrate that the benefit of conducting longitudinal studies depends on the strengths of the between- and within-subject associations of the brain and non-brain measures. Explicitly modeling between- and within-subject effects avoids conflating the effects and allows optimizing effect sizes for them separately. These findings underscore the importance of considering study design features to improve the replicability of BWAS.

Neural Circuit Markers of Familial Risk for Depression Among Healthy Youths in the Adolescent Brain Cognitive Development Study.

BACKGROUND: Family history of depression is a robust predictor of early-onset depression, which may confer risk through alterations in neural circuits that have been implicated in reward and emotional processing. These alterations may be evident in youths who are at familial risk for depression but who do not currently have depression. However, the identification of robust and replicable findings has been hindered by few studies and small sample sizes. In the current study, we sought to identify functional connectivity (FC) patterns associated with familial risk for depression. METHODS: Participants included healthy (i.e., no lifetime psychiatric diagnoses) youths at high familial risk for depression (HR) (n = 754; at least one parent with a history of depression) and healthy youths at low familial risk for psychiatric problems (LR) (n = 1745; no parental history of psychopathology) who were 9 to 10 years of age and from the Adolescent Brain Cognitive Development (ABCD) Study sample. We conducted whole-brain seed-to-voxel analyses to examine group differences in resting-state FC with the amygdala, caudate, nucleus accumbens, and putamen. We hypothesized that HR youths would exhibit global amygdala hyperconnectivity and striatal hypoconnectivity patterns primarily driven by maternal risk. RESULTS: HR youths exhibited weaker caudate-angular gyrus FC than LR youths (α = 0.04, Cohen's d = 0.17). HR youths with a history of maternal depression specifically exhibited weaker caudate-angular gyrus FC (α = 0.03, Cohen's d = 0.19) as well as weaker caudate-dorsolateral prefrontal cortex FC (α = 0.04, Cohen's d = 0.21) than LR youths. CONCLUSIONS: Weaker striatal connectivity may be related to heightened familial risk for depression, primarily driven by maternal history. Identifying brain-based markers of depression risk in youths can inform approaches to improving early detection, diagnosis, and treatment.

XCP-D: A Robust Pipeline for the post-processing of fMRI data.

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

Personalized functional brain network topography is associated with individual differences in youth cognition.

Individual differences in cognition during childhood are associated with important social, physical, and mental health outcomes in adolescence and adulthood. Given that cortical surface arealization during development reflects the brain's functional prioritization, quantifying variation in the topography of functional brain networks across the developing cortex may provide insight regarding individual differences in cognition. We test this idea by defining personalized functional networks (PFNs) that account for interindividual heterogeneity in functional brain network topography in 9-10 year olds from the Adolescent Brain Cognitive Development℠ Study. Across matched discovery (n = 3525) and replication (n = 3447) samples, the total cortical representation of fronto-parietal PFNs positively correlates with general cognition. Cross-validated ridge regressions trained on PFN topography predict cognition in unseen data across domains, with prediction accuracy increasing along the cortex's sensorimotor-association organizational axis. These results establish that functional network topography heterogeneity is associated with individual differences in cognition before the critical transition into adolescence.

Motor networks, but also non-motor networks predict motor signs in Parkinson's disease.

OBJECTIVE: Investigate the brain functional networks associated with motor impairment in people with Parkinson's disease (PD). BACKGROUND: PD is primarily characterized by motor dysfunction. Resting-state functional connectivity (RsFC) offers a unique opportunity to non-invasively characterize brain function. In this study, we hypothesized that the motor dysfunction observed in people with PD involves atypical connectivity not only in motor but also in higher-level attention networks. Understanding the interaction between motor and non-motor RsFC that are related to the motor signs could provide insights into PD pathophysiology. METHODS: We used data from 88 people with PD (mean age: 68.2(SD:10), 55 M/33F) coming from 2 cohorts. Motor severity was assessed in practical OFF-medication state, using MDS-UPDRS Part-III motor scores (mean: 49 (SD:10)). RsFC was characterized using an atlas of 384 regions that were grouped into 13 functional networks. Associations between RsFC and motor severity were assessed independently for each RsFC using predictive modeling. RESULTS: The top 5 % models that predicted the MDS-UPDRS-III motor scores with effect size >0.5 were the connectivity between (1) the somatomotor and Subcortical-Basal-ganglia, (2) somatomotor and Visual and (3) CinguloOpercular (CiO) and language/Ventral attention (Lan/VeA) network pairs. DISCUSSION: Our findings suggest that, along with motor networks, visual- and attention-related cortical networks are also associated with the motor symptoms of PD. Non-motor networks may be involved indirectly in motor-coordination. When people with PD have deficits in motor networks, more attention may be needed to carry out formerly automatic motor functions, consistent with compensatory mechanisms in parkinsonian movement disorders.

Early-in-life isoflurane exposure alters resting-state functional connectivity in juvenile non-human primates.

BACKGROUND: Clinical studies suggest that anaesthesia exposure early in life affects neurobehavioural development. We designed a non-human primate (NHP) study to evaluate cognitive, behavioural, and brain functional and structural alterations after isoflurane exposure during infancy. These NHPs displayed decreased close social behaviour and increased astrogliosis in specific brain regions, most notably in the amygdala. Here we hypothesise that resting-state functional connectivity MRI can detect alterations in connectivity of brain areas that relate to these social behaviours and astrogliosis. METHODS: Imaging was performed in 2-yr-old NHPs under light anaesthesia, after early-in-life (postnatal days 6-12) exposure to 5 h of isoflurane either one or three times, or to room air. Brain images were segmented into 82 regions of interest; the amygdala and the posterior cingulate cortex were chosen for a seed-based resting-state functional connectivity MRI analysis. RESULTS: We found differences between groups in resting-state functional connectivity of the amygdala and the auditory cortices, medial premotor cortex, and posterior cingulate cortex. There were also alterations in resting-state functional connectivity between the posterior cingulate cortex and secondary auditory, polar prefrontal, and temporal cortices, and the anterior insula. Relationships were identified between resting-state functional connectivity alterations and the decrease in close social behaviour and increased astrogliosis. CONCLUSIONS: Early-in-life anaesthesia exposure in NHPs is associated with resting-state functional connectivity alterations of the amygdala and the posterior cingulate cortex with other brain regions, evident at the juvenile age of 2 yr. These changes in resting-state functional connectivity correlate with the decrease in close social behaviour and increased astrogliosis. Using resting-state functional connectivity MRI to study the neuronal underpinnings of early-in-life anaesthesia-induced behavioural alterations could facilitate development of a biomarker for anaesthesia-induced developmental neurotoxicity.

Using synthetic MR images for distortion correction.

Functional MRI (fMRI) data acquired using echo-planar imaging (EPI) are highly distorted by magnetic field inhomogeneities. Distortion and differences in image contrast between EPI and T1-weighted and T2-weighted (T1w/T2w) images makes their alignment a challenge. Typically, field map data are used to correct EPI distortions. Alignments achieved with field maps can vary greatly and depends on the quality of field map data. However, many public datasets lack field map data entirely. Additionally, reliable field map data is often difficult to acquire in high-motion pediatric or developmental cohorts. To address this, we developed Synth, a software package for distortion correction and cross-modal image registration that does not require field map data. Synth combines information from T1w and T2w anatomical images to construct an idealized undistorted synthetic image with similar contrast properties to EPI data. This synthetic image acts as an effective reference for individual-specific distortion correction. Using pediatric (ABCD: Adolescent Brain Cognitive Development) and adult (MSC: Midnight Scan Club; HCP: Human Connectome Project) data, we demonstrate that Synth performs comparably to field map distortion correction approaches, and often outperforms them. Field map-less distortion correction with Synth allows accurate and precise registration of fMRI data with missing or corrupted field map information.

Polyneuro risk scores capture widely distributed connectivity patterns of cognition.

Resting-state functional connectivity (RSFC) is a powerful tool for characterizing brain changes, but it has yet to reliably predict higher-order cognition. This may be attributed to small effect sizes of such brain-behavior relationships, which can lead to underpowered, variable results when utilizing typical sample sizes (N∼25). Inspired by techniques in genomics, we implement the polyneuro risk score (PNRS) framework - the application of multivariate techniques to RSFC data and validation in an independent sample. Utilizing the Adolescent Brain Cognitive Development® cohort split into two datasets, we explore the framework's ability to reliably capture brain-behavior relationships across 3 cognitive scores - general ability, executive function, learning & memory. The weight and significance of each connection is assessed in the first dataset, and a PNRS is calculated for each participant in the second. Results support the PNRS framework as a suitable methodology to inspect the distribution of connections contributing towards behavior, with explained variance ranging from 1.0 % to 21.4 %. For the outcomes assessed, the framework reveals globally distributed, rather than localized, patterns of predictive connections. Larger samples are likely necessary to systematically identify the specific connections contributing towards complex outcomes. The PNRS framework could be applied translationally to identify neurologically distinct subtypes of neurodevelopmental disorders.

Functional maturation in visual pathways predicts attention to the eyes in infant rhesus macaques: Effects of social status.

Differences in looking at the eyes of others are one of the earliest behavioral markers for social difficulties in neurodevelopmental disabilities, including autism. However, it is unknown how early visuo-social experiences relate to the maturation of infant brain networks that process visual social stimuli. We investigated functional connectivity (FC) within the ventral visual object pathway as a contributing neural system. Densely sampled, longitudinal eye-tracking and resting state fMRI (rs-fMRI) data were collected from infant rhesus macaques, an important model of human social development, from birth through 6 months of age. Mean trajectories were fit for both datasets and individual trajectories from subjects with both eye-tracking and rs-fMRI data were used to test for brain-behavior relationships. Exploratory findings showed infants with greater increases in FC between left V1 to V3 visual areas have an earlier increase in eye-looking before 2 months. This relationship was moderated by social status such that infants with low social status had a stronger association between left V1 to V3 connectivity and eye-looking than high status infants. Results indicated that maturation of the visual object pathway may provide an important neural substrate supporting adaptive transitions in social visual attention during infancy.

Does Cueing Need Attention? A Pilot Study in People with Parkinson's Disease.

We previously showed that both open-loop (beat of a metronome) and closed-loop (phase-dependent tactile feedback) cueing may be similarly effective in reducing Freezing of Gait (FoG), assessed with a quantitative FoG Index, while turning in place in the laboratory in a group of people with Parkinson's disease (PD). Despite the similar changes on the FoG Index, it is not known whether both cueing responses require attentional control, which would explain FoG Index improvement. The mechanisms underlying cueing responses are poorly understood. Here, we tested the hypothesis that the salience network would predict responsiveness (i.e., FoG Index improvement) to open-loop and closed-loop cueing in people with and without FoG of PD, as salience network contributes to tasks requiring attention to external stimuli in healthy adults. Thirteen people with PD with high-quality imaging data were analyzed to characterize relationships between resting-state MRI functional connectivity and responses to cues. The interaction of the salience network and retrosplenial-temporal networks was the best predictor of responsiveness to open-loop cueing, presenting the largest effect size (d = 1.16). The interaction between the salience network and subcortical as well as cingulo-parietal and subcortical networks were the strongest predictors of responsiveness to closed-loop cueing, presenting the largest effect sizes (d = 1.06 and d = 0.84, respectively). Salience network activity was a common predictor of responsiveness to both cueing, which suggests that auditory and proprioceptive stimuli during turning may require some level of cognitive and insular activity, anchored within the salience network, which explain FoG Index improvements in people with PD.

Use of connectotyping on task functional MRI data reveals dynamic network level cross talking during task performance.

Task-based functional MRI (fMRI) has greatly improved understanding of brain functioning, enabling the identification of brain areas associated with specific cognitive operations. Traditional analyses are limited to associating activation patterns in particular regions with specific cognitive operation, largely ignoring regional cross-talk or dynamic connectivity, which we propose is crucial for characterization of brain function in the context of task fMRI. We use connectotyping, which efficiently models functional brain connectivity to reveal the progression of temporal brain connectivity patterns in task fMRI. Connectotyping was employed on data from twenty-four participants (12 male, mean age 24.8 years, 2.57 std. dev) who performed a widely spaced event-related fMRI word vs. pseudoword decision task, where stimuli were presented every 20 s. After filtering for movement, we ended up with 15 participants that completed each trial and had enough usable data for our analyses. Connectivity matrices were calculated per participant across time for each stimuli type. A Repeated Measures ANOVA applied on the connectotypes was used to characterize differences across time for words and pseudowords. Our group level analyses found significantly different dynamic connectivity patterns during word vs. pseudoword processing between the Fronto-Parietal and Cingulo-Parietal Systems, areas involved in cognitive task control, memory retrieval, and semantic processing. Our findings support the presence of dynamic changes in functional connectivity during task execution and that such changes can be characterized using connectotyping but not with traditional Pearson's correlations.

Carotenoids improve the development of cerebral cortical networks in formula-fed infant macaques.

Nutrition during the first years of life has a significant impact on brain development. This study characterized differences in brain maturation from birth to 6 months of life in infant macaques fed formulas differing in content of lutein, ß-carotene, and other carotenoids using Magnetic Resonance Imaging to measure functional connectivity. We observed differences in functional connectivity based on the interaction of diet, age and brain networks. Post hoc analysis revealed significant diet-specific differences between insular-opercular and somatomotor networks at 2 months of age, dorsal attention and somatomotor at 4 months of age, and within somatomotor and between somatomotor-visual and auditory-dorsal attention networks at 6 months of age. Overall, we found a larger divergence in connectivity from the breastfeeding group in infant macaques fed formula containing no supplemental carotenoids in comparison to those fed formula supplemented with carotenoids. These findings suggest that carotenoid formula supplementation influences functional brain development.