Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain.

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

Mapping human brain charts cross-sectionally and longitudinally.

Brain scans acquired across large, age-diverse cohorts have facilitated recent progress in establishing normative brain aging charts. Here, we ask the critical question of whether cross-sectional estimates of age-related brain trajectories resemble those directly measured from longitudinal data. We show that age-related brain changes inferred from cross-sectionally mapped brain charts can substantially underestimate actual changes measured longitudinally. We further find that brain aging trajectories vary markedly between individuals and are difficult to predict with population-level age trends estimated cross-sectionally. Prediction errors relate modestly to neuroimaging confounds and lifestyle factors. Our findings provide explicit evidence for the importance of longitudinal measurements in ascertaining brain development and aging trajectories.

Sex differences in the functional topography of association networks in youth.

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Traumatic stress load and stressor reactivity score associated with accelerated gray matter maturation in youths indexed by normative models.

Understanding how traumatic stress affects typical brain development during adolescence is critical to elucidate underlying mechanisms related to both maladaptive functioning and resilience after traumatic exposures. The current study aimed to map deviations from normative ranges of brain gray matter for youths with traumatic exposures. For each cortical and subcortical gray matter region, normative percentiles of variations were established using structural MRI from typically developing youths without any traumatic exposure (n = 245; age range = 8-23) from the Philadelphia Neurodevelopmental Cohort (PNC). The remaining PNC participants with neuroimaging data (n = 1129) were classified as either within the normative range (5-95%), delayed (>95%) or accelerated (<5%) maturational ranges for each region using the normative model. An averaged quantile regression index was calculated across all regions. Mediation models revealed that high traumatic stress load was positively associated with poorer cognitive functioning and greater psychopathology, and these associations were mediated by accelerated gray matter maturation. Furthermore, higher stressor reactivity scores, which represent a less resilient response under traumatic stress, were positively correlated with greater acceleration of gray matter maturation (r = 0.224, 95% CI = [0.17, 0.28], p < 0.001), suggesting that more accelerated maturation was linked to greater stressor response regardless of traumatic stress load. We conclude that traumatic stress is a source of deviation from normative brain development associated with poorer cognitive functioning and more psychopathology in the long run.

Associations between neighborhood socioeconomic status, parental education, and executive system activation in youth.

Socioeconomic status (SES) can impact cognitive performance, including working memory (WM). As executive systems that support WM undergo functional neurodevelopment during adolescence, environmental stressors at both individual and community levels may influence cognitive outcomes. Here, we sought to examine how SES at the neighborhood and family level impacts task-related activation of the executive system during adolescence and determine whether this effect mediates the relationship between SES and WM performance. To address these questions, we studied 1,150 youths (age 8-23) that completed a fractal n-back WM task during functional magnetic resonance imaging at 3T as part of the Philadelphia Neurodevelopmental Cohort. We found that both higher neighborhood SES and parental education were associated with greater activation of the executive system to WM load, including the bilateral dorsolateral prefrontal cortex, posterior parietal cortex, and precuneus. The association of neighborhood SES remained significant when controlling for task performance, or related factors like exposure to traumatic events. Furthermore, high-dimensional multivariate mediation analysis identified distinct patterns of brain activity within the executive system that significantly mediated the relationship between measures of SES and task performance. These findings underscore the importance of multilevel environmental factors in shaping executive system function and WM in youth.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization.

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

CoCoA: conditional correlation models with association size.

Many scientific questions can be formulated as hypotheses about conditional correlations. For instance, in tests of cognitive and physical performance, the trade-off between speed and accuracy motivates study of the two variables together. A natural question is whether speed-accuracy coupling depends on other variables, such as sustained attention. Classical regression techniques, which posit models in terms of covariates and outcomes, are insufficient to investigate the effect of a third variable on the symmetric relationship between speed and accuracy. In response, we propose a conditional correlation model with association size, a likelihood-based statistical framework to estimate the conditional correlation between speed and accuracy as a function of additional variables. We propose novel measures of the association size, which are analogous to effect sizes on the correlation scale while adjusting for confound variables. In simulation studies, we compare likelihood-based estimators of conditional correlation to semiparametric estimators adapted from genomic studies and find that the former achieves lower bias and variance under both ideal settings and model assumption misspecification. Using neurocognitive data from the Philadelphia Neurodevelopmental Cohort, we demonstrate that greater sustained attention is associated with stronger speed-accuracy coupling in a complex reasoning task while controlling for age. By highlighting conditional correlations as the outcome of interest, our model provides complementary insights to traditional regression modeling and partitioned correlation analyses.

Brain Growth Charts for Quantitative Analysis of Pediatric Clinical Brain MRI Scans with Limited Imaging Pathology.

Background Clinically acquired brain MRI scans represent a valuable but underused resource for investigating neurodevelopment due to their technical heterogeneity and lack of appropriate controls. These barriers have curtailed retrospective studies of clinical brain MRI scans compared with more costly prospectively acquired research-quality brain MRI scans. Purpose To provide a benchmark for neuroanatomic variability in clinically acquired brain MRI scans with limited imaging pathology (SLIPs) and to evaluate if growth charts from curated clinical MRI scans differed from research-quality MRI scans or were influenced by clinical indication for the scan. Materials and Methods In this secondary analysis of preexisting data, clinical brain MRI SLIPs from an urban pediatric health care system (individuals aged ≤22 years) were scanned across nine 3.0-T MRI scanners. The curation process included manual review of signed radiology reports and automated and manual quality review of images without gross pathology. Global and regional volumetric imaging phenotypes were measured using two image segmentation pipelines, and clinical brain growth charts were quantitatively compared with charts derived from a large set of research controls in the same age range by means of Pearson correlation and age at peak volume. Results The curated clinical data set included 532 patients (277 male; median age, 10 years [IQR, 5-14 years]; age range, 28 days after birth to 22 years) scanned between 2005 and 2020. Clinical brain growth charts were highly correlated with growth charts derived from research data sets (22 studies, 8346 individuals [4947 male]; age range, 152 days after birth to 22 years) in terms of normative developmental trajectories predicted by the models (median r = 0.979). Conclusion The clinical indication of the scans did not significantly bias the output of clinical brain charts. Brain growth charts derived from clinical controls with limited imaging pathology were highly correlated with brain charts from research controls, suggesting the potential of curated clinical MRI scans to supplement research data sets. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Ertl-Wagner and Pai in this issue.

Development of structure-function coupling in human brain networks during youth.

The protracted development of structural and functional brain connectivity within distributed association networks coincides with improvements in higher-order cognitive processes such as executive function. However, it remains unclear how white-matter architecture develops during youth to directly support coordinated neural activity. Here, we characterize the development of structure-function coupling using diffusion-weighted imaging and n-back functional MRI data in a sample of 727 individuals (ages 8 to 23 y). We found that spatial variability in structure-function coupling aligned with cortical hierarchies of functional specialization and evolutionary expansion. Furthermore, hierarchy-dependent age effects on structure-function coupling localized to transmodal cortex in both cross-sectional data and a subset of participants with longitudinal data (n = 294). Moreover, structure-function coupling in rostrolateral prefrontal cortex was associated with executive performance and partially mediated age-related improvements in executive function. Together, these findings delineate a critical dimension of adolescent brain development, whereby the coupling between structural and functional connectivity remodels to support functional specialization and cognition.

Imaging local genetic influences on cortical folding.

Recent progress in deciphering mechanisms of human brain cortical folding leave unexplained whether spatially patterned genetic influences contribute to this folding. High-resolution in vivo brain MRI can be used to estimate genetic correlations (covariability due to shared genetic factors) in interregional cortical thickness, and biomechanical studies predict an influence of cortical thickness on folding patterns. However, progress has been hampered because shared genetic influences related to folding patterns likely operate at a scale that is much more local (<1 cm) than that addressed in prior imaging studies. Here, we develop methodological approaches to examine local genetic influences on cortical thickness and apply these methods to two large, independent samples. We find that such influences are markedly heterogeneous in strength, and in some cortical areas are notably stronger in specific orientations relative to gyri or sulci. The overall, phenotypic local correlation has a significant basis in shared genetic factors and is highly symmetric between left and right cortical hemispheres. Furthermore, the degree of local cortical folding relates systematically with the strength of local correlations, which tends to be higher in gyral crests and lower in sulcal fundi. The relationship between folding and local correlations is stronger in primary sensorimotor areas and weaker in association areas such as prefrontal cortex, consistent with reduced genetic constraints on the structural topology of association cortex. Collectively, our results suggest that patterned genetic influences on cortical thickness, measurable at the scale of in vivo MRI, may be a causal factor in the development of cortical folding.

Spatially-enhanced clusterwise inference for testing and localizing intermodal correspondence.

With the increasing availability of neuroimaging data from multiple modalities-each providing a different lens through which to study brain structure or function-new techniques for comparing, integrating, and interpreting information within and across modalities have emerged. Recent developments include hypothesis tests of associations between neuroimaging modalities, which can be used to determine the statistical significance of intermodal associations either throughout the entire brain or within anatomical subregions or functional networks. While these methods provide a crucial foundation for inference on intermodal relationships, they cannot be used to answer questions about where in the brain these associations are most pronounced. In this paper, we introduce a new method, called CLEAN-R, that can be used both to test intermodal correspondence throughout the brain and also to localize this correspondence. Our method involves first adjusting for the underlying spatial autocorrelation structure within each modality before aggregating information within small clusters to construct a map of enhanced test statistics. Using structural and functional magnetic resonance imaging data from a subsample of children and adolescents from the Philadelphia Neurodevelopmental Cohort, we conduct simulations and data analyses where we illustrate the high statistical power and nominal type I error levels of our method. By constructing an interpretable map of group-level correspondence using spatially-enhanced test statistics, our method offers insights beyond those provided by earlier methods.

Curation of BIDS (CuBIDS): A workflow and software package for streamlining reproducible curation of large BIDS datasets.

The Brain Imaging Data Structure (BIDS) is a specification accompanied by a software ecosystem that was designed to create reproducible and automated workflows for processing neuroimaging data. BIDS Apps flexibly build workflows based on the metadata detected in a dataset. However, even BIDS valid metadata can include incorrect values or omissions that result in inconsistent processing across sessions. Additionally, in large-scale, heterogeneous neuroimaging datasets, hidden variability in metadata is difficult to detect and classify. To address these challenges, we created a Python-based software package titled "Curation of BIDS" (CuBIDS), which provides an intuitive workflow that helps users validate and manage the curation of their neuroimaging datasets. CuBIDS includes a robust implementation of BIDS validation that scales to large samples and incorporates DataLad--a version control software package for data--as an optional dependency to ensure reproducibility and provenance tracking throughout the entire curation process. CuBIDS provides tools to help users perform quality control on their images' metadata and identify unique combinations of imaging parameters. Users can then execute BIDS Apps on a subset of participants that represent the full range of acquisition parameters that are present, accelerating pipeline testing on large datasets.

Voxel-wise intermodal coupling analysis of two or more modalities using local covariance decomposition.

When individual subjects are imaged with multiple modalities, biological information is present not only within each modality, but also between modalities - that is, in how modalities covary at the voxel level. Previous studies have shown that local covariance structures between modalities, or intermodal coupling (IMCo), can be summarized for two modalities, and that two-modality IMCo reveals otherwise undiscovered patterns in neurodevelopment and certain diseases. However, previous IMCo methods are based on the slopes of local weighted linear regression lines, which are inherently asymmetric and limited to the two-modality setting. Here, we present a generalization of IMCo estimation which uses local covariance decompositions to define a symmetric, voxel-wise coupling coefficient that is valid for two or more modalities. We use this method to study coupling between cerebral blood flow, amplitude of low frequency fluctuations, and local connectivity in 803 subjects ages 8 through 22. We demonstrate that coupling is spatially heterogeneous, varies with respect to age and sex in neurodevelopment, and reveals patterns that are not present in individual modalities. As availability of multi-modal data continues to increase, principal-component-based IMCo (pIMCo) offers a powerful approach for summarizing relationships between multiple aspects of brain structure and function. An R package for estimating pIMCo is available at: https://github.com/hufengling/pIMCo.

A simple permutation-based test of intermodal correspondence.

Many key findings in neuroimaging studies involve similarities between brain maps, but statistical methods used to measure these findings have varied. Current state-of-the-art methods involve comparing observed group-level brain maps (after averaging intensities at each image location across multiple subjects) against spatial null models of these group-level maps. However, these methods typically make strong and potentially unrealistic statistical assumptions, such as covariance stationarity. To address these issues, in this article we propose using subject-level data and a classical permutation testing framework to test and assess similarities between brain maps. Our method is comparable to traditional permutation tests in that it involves randomly permuting subjects to generate a null distribution of intermodal correspondence statistics, which we compare to an observed statistic to estimate a p-value. We apply and compare our method in simulated and real neuroimaging data from the Philadelphia Neurodevelopmental Cohort. We show that our method performs well for detecting relationships between modalities known to be strongly related (cortical thickness and sulcal depth), and it is conservative when an association would not be expected (cortical thickness and activation on the n-back working memory task). Notably, our method is the most flexible and reliable for localizing intermodal relationships within subregions of the brain and allows for generalizable statistical inference.

Associations of cannabis use disorder with cognition, brain structure, and brain function in African Americans.

Although previous studies have highlighted associations of cannabis use with cognition and brain morphometry, critical questions remain with regard to the association between cannabis use and brain structural and functional connectivity. In a cross-sectional community sample of 205 African Americans (age 18-70) we tested for associations of cannabis use disorder (CUD, n = 57) with multi-domain cognitive measures and structural, diffusion, and resting state brain-imaging phenotypes. Post hoc model evidence was computed with Bayes factors (BF) and posterior probabilities of association (PPA) to account for multiple testing. General cognitive functioning, verbal intelligence, verbal memory, working memory, and motor speed were lower in the CUD group compared with non-users (p < .011; 1.9 < BF < 3,217). CUD was associated with altered functional connectivity in a network comprising the motor-hand region in the superior parietal gyri and the anterior insula (p < .04). These differences were not explained by alcohol, other drug use, or education. No associations with CUD were observed in cortical thickness, cortical surface area, subcortical or cerebellar volumes (0.12 < BF < 1.5), or graph-theoretical metrics of resting state connectivity (PPA < 0.01). In a large sample collected irrespective of cannabis used to minimize recruitment bias, we confirm the literature on poorer cognitive functioning in CUD, and an absence of volumetric brain differences between CUD and non-CUD. We did not find evidence for or against a disruption of structural connectivity, whereas we did find localized resting state functional dysconnectivity in CUD. There was sufficient proof, however, that organization of functional connectivity as determined via graph metrics does not differ between CUD and non-user group.

Childhood Obesity, Cortical Structure, and Executive Function in Healthy Children.

The development of executive function is linked to maturation of prefrontal cortex (PFC) in childhood. Childhood obesity has been associated with changes in brain structure, particularly in PFC, as well as deficits in executive functions. We aimed to determine whether differences in cortical structure mediate the relationship between executive function and childhood obesity. We analyzed MR-derived measures of cortical thickness for 2700 children between the ages of 9 and 11 years, recruited as part of the NIH Adolescent Brain and Cognitive Development (ABCD) study. We related our findings to measures of executive function and body mass index (BMI). In our analysis, increased BMI was associated with significantly reduced mean cortical thickness, as well as specific bilateral reduced cortical thickness in prefrontal cortical regions. This relationship remained after accounting for age, sex, race, parental education, household income, birth-weight, and in-scanner motion. Increased BMI was also associated with lower executive function. Reduced thickness in the rostral medial and superior frontal cortex, the inferior frontal gyrus, and the lateral orbitofrontal cortex partially accounted for reductions in executive function. These results suggest that childhood obesity is associated with compromised executive function. This relationship may be partly explained by BMI-associated reduced cortical thickness in the PFC.

Altered Sex Chromosome Dosage Induces Coordinated Shifts in Cortical Anatomy and Anatomical Covariance.

Sex chromosome dosage (SCD) variation increases risk for neuropsychiatric impairment, which may reflect direct SCD effects on brain organization. Here, we 1) map cumulative X- and Y-chromosome dosage effects on regional cortical thickness (CT) and investigate potential functional implications of these effects using Neurosynth, 2) test if this map is organized by patterns of CT covariance that are evident in health, and 3) characterize SCD effects on CT covariance itself. We modeled SCD effects on CT and CT covariance for 308 equally sized regions of the cortical sheet using structural neuroimaging data from 301 individuals with varying numbers of sex chromosomes (169 euploid, 132 aneuploid). Mounting SCD increased CT in the rostral frontal cortex and decreased CT in the lateral temporal cortex, bilaterally. Regions targeted by SCD were associated with social functioning, language processing, and comprehension. Cortical regions with a similar degree of SCD-sensitivity showed heightened CT covariance in health. Finally, greater SCD also increased covariance among regions similarly affected by SCD. Our study both 1) develops novel methods for comparing typical and disease-related structural covariance networks in the brain and 2) uses these techniques to resolve and identify organizing principles for SCD effects on regional cortical anatomy and anatomical covariance.

Genetic Contributions to Multivariate Data-Driven Brain Networks Constructed via Source-Based Morphometry.

Identifying genetic factors underlying neuroanatomical variation has been difficult. Traditional methods have used brain regions from predetermined parcellation schemes as phenotypes for genetic analyses, although these parcellations often do not reflect brain function and/or do not account for covariance between regions. We proposed that network-based phenotypes derived via source-based morphometry (SBM) may provide additional insight into the genetic architecture of neuroanatomy given its data-driven approach and consideration of covariance between voxels. We found that anatomical SBM networks constructed on ~ 20 000 individuals from the UK Biobank were heritable and shared functionally meaningful genetic overlap with each other. We additionally identified 27 unique genetic loci that contributed to one or more SBM networks. Both GWA and genetic correlation results indicated complex patterns of pleiotropy and polygenicity similar to other complex traits. Lastly, we found genetic overlap between a network related to the default mode and schizophrenia, a disorder commonly associated with neuroanatomic alterations.

Minimal Relationship between Local Gyrification and General Cognitive Ability in Humans.

Previous studies suggest that gyrification is associated with superior cognitive abilities in humans, but the strength of this relationship remains unclear. Here, in two samples of related individuals (total N = 2882), we calculated an index of local gyrification (LGI) at thousands of cortical surface points using structural brain images and an index of general cognitive ability (g) using performance on cognitive tests. Replicating previous studies, we found that phenotypic and genetic LGI-g correlations were positive and statistically significant in many cortical regions. However, all LGI-g correlations in both samples were extremely weak, regardless of whether they were significant or nonsignificant. For example, the median phenotypic LGI-g correlation was 0.05 in one sample and 0.10 in the other. These correlations were even weaker after adjusting for confounding neuroanatomical variables (intracranial volume and local cortical surface area). Furthermore, when all LGIs were considered together, at least 89% of the phenotypic variance of g remained unaccounted for. We conclude that the association between LGI and g is too weak to have profound implications for our understanding of the neurobiology of intelligence. This study highlights potential issues when focusing heavily on statistical significance rather than effect sizes in large-scale observational neuroimaging studies.