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.

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.

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.

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.

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.

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.

Sex-chromosome dosage effects on gene expression in humans.

A fundamental question in the biology of sex differences has eluded direct study in humans: How does sex-chromosome dosage (SCD) shape genome function? To address this, we developed a systematic map of SCD effects on gene function by analyzing genome-wide expression data in humans with diverse sex-chromosome aneuploidies (XO, XXX, XXY, XYY, and XXYY). For sex chromosomes, we demonstrate a pattern of obligate dosage sensitivity among evolutionarily preserved X-Y homologs and update prevailing theoretical models for SCD compensation by detecting X-linked genes that increase expression with decreasing X- and/or Y-chromosome dosage. We further show that SCD-sensitive sex-chromosome genes regulate specific coexpression networks of SCD-sensitive autosomal genes with critical cellular functions and a demonstrable potential to mediate previously documented SCD effects on disease. These gene coexpression results converge with analysis of transcription factor binding site enrichment and measures of gene expression in murine knockout models to spotlight the dosage-sensitive X-linked transcription factor ZFX as a key mediator of SCD effects on wider genome expression. Our findings characterize the effects of SCD broadly across the genome, with potential implications for human phenotypic variation.

Human Cortical Thickness Organized into Genetically-determined Communities across Spatial Resolutions.

The cerebral cortex may be organized into anatomical genetic modules, communities of brain regions with shared genetic influences via pleiotropy. Such modules could represent novel phenotypes amenable to large-scale gene discovery. This modular structure was investigated with network analysis of in vivo MRI of extended pedigrees, revealing a "multiscale" structure where smaller and larger modules exist simultaneously and in partially overlapping fashion across spatial scales, in contrast to prior work suggesting a specific number of cortical thickness modules. Inter-regional genetic correlations, gene co-expression patterns and computational models indicate that two simple organizational principles account for a large proportion of the apparent complexity in the network of genetic correlations. First, regions are strongly genetically correlated with their homologs in the opposite cerebral hemisphere. Second, regions are strongly genetically correlated with nearby regions in the same hemisphere, with an initial steep decrease in genetic correlation with anatomical distance, followed by a more gradual decline. Understanding underlying organizational principles of genetic influence is a critical step towards a mechanistic model of how specific genes influence brain anatomy and mediate neuropsychiatric risk.

Waves of Maturation and Senescence in Micro-structural MRI Markers of Human Cortical Myelination over the Lifespan.

Seminal human brain histology work has demonstrated developmental waves of myelination. Here, using a micro-structural magnetic resonance imaging (MRI) marker linked to myelin, we studied fine-grained age differences to deduce waves of growth, stability, and decline of cortical myelination over the life-cycle. In 484 participants, aged 8-85 years, we fitted smooth growth curves to T1- to T2-weighted ratio in each of 360 regions from one of seven cytoarchitectonic classes. From the first derivatives of these generally inverted-U trajectories, we defined three milestones: the age at peak growth; the age at onset of a stable plateau; and the age at the onset of decline. Age at peak growth had a bimodal distribution comprising an early (pre-pubertal) wave of primary sensory and motor cortices and a later (post-pubertal) wave of association, insular and limbic cortices. Most regions reached stability in the 30-s but there was a second wave reaching stability in the 50-s. Age at onset of decline was also bimodal: in some right hemisphere regions, the curve declined from the 60-s, but in other left hemisphere regions, there was no significant decline from the stable plateau. These results are consistent with regionally heterogeneous waves of intracortical myelinogenesis and age-related demyelination.

On testing for spatial correspondence between maps of human brain structure and function.

A critical issue in many neuroimaging studies is the comparison between brain maps. Nonetheless, it remains unclear how one should test hypotheses focused on the overlap or spatial correspondence between two or more brain maps. This "correspondence problem" affects, for example, the interpretation of comparisons between task-based patterns of functional activation, resting-state networks or modules, and neuroanatomical landmarks. To date, this problem has been addressed with remarkable variability in terms of methodological approaches and statistical rigor. In this paper, we address the correspondence problem using a spatial permutation framework to generate null models of overlap by applying random rotations to spherical representations of the cortical surface, an approach for which we also provide a theoretical statistical foundation. We use this method to derive clusters of cognitive functions that are correlated in terms of their functional neuroatomical substrates. In addition, using publicly available data, we formally demonstrate the correspondence between maps of task-based functional activity, resting-state fMRI networks and gyral-based anatomical landmarks. We provide open-access code to implement the methods presented for two commonly-used tools for surface based cortical analysis (https://www.github.com/spin-test). This spatial permutation approach constitutes a useful advance over widely-used methods for the comparison of cortical maps, thereby opening new possibilities for the integration of diverse neuroimaging data.

Impaired long distance functional connectivity and weighted network architecture in Alzheimer's disease.

Alzheimer's disease (AD) is increasingly recognized as a disconnection syndrome, which leads to cognitive impairment due to the disruption of functional activity across large networks or systems of interconnected brain regions. We explored abnormal functional magnetic resonance imaging (fMRI) resting-state dynamics, functional connectivity, and weighted functional networks, in a sample of patients with severe AD (N = 18) and age-matched healthy volunteers (N = 21). We found that patients had reduced amplitude and regional homogeneity of low-frequency fMRI oscillations, and reduced the strength of functional connectivity, in several regions previously described as components of the default mode network, for example, medial posterior parietal cortex and dorsal medial prefrontal cortex. In patients with severe AD, functional connectivity was particularly attenuated between regions that were separated by a greater physical distance; and loss of long distance connectivity was associated with less efficient global and nodal network topology. This profile of functional abnormality in severe AD was consistent with the results of a comparable analysis of data on 2 additional groups of patients with mild AD (N = 17) and amnestic mild cognitive impairment (MCI; N = 18). A greater degree of cognitive impairment, measured by the mini-mental state examination across all patient groups, was correlated with greater attenuation of functional connectivity, particularly over long connection distances, for example, between anterior and posterior components of the default mode network, and greater reduction of global and nodal network efficiency. These results indicate that neurodegenerative disruption of fMRI oscillations and connectivity in AD affects long-distance connections to hub nodes, with the consequent loss of network efficiency. This profile was evident also to a lesser degree in the patients with less severe cognitive impairment, indicating that the potential of resting-state fMRI measures as biomarkers or predictors of disease progression in AD.

Simple models of human brain functional networks.

Human brain functional networks are embedded in anatomical space and have topological properties--small-worldness, modularity, fat-tailed degree distributions--that are comparable to many other complex networks. Although a sophisticated set of measures is available to describe the topology of brain networks, the selection pressures that drive their formation remain largely unknown. Here we consider generative models for the probability of a functional connection (an edge) between two cortical regions (nodes) separated by some Euclidean distance in anatomical space. In particular, we propose a model in which the embedded topology of brain networks emerges from two competing factors: a distance penalty based on the cost of maintaining long-range connections; and a topological term that favors links between regions sharing similar input. We show that, together, these two biologically plausible factors are sufficient to capture an impressive range of topological properties of functional brain networks. Model parameters estimated in one set of functional MRI (fMRI) data on normal volunteers provided a good fit to networks estimated in a second independent sample of fMRI data. Furthermore, slightly detuned model parameters also generated a reasonable simulation of the abnormal properties of brain functional networks in people with schizophrenia. We therefore anticipate that many aspects of brain network organization, in health and disease, may be parsimoniously explained by an economical clustering rule for the probability of functional connectivity between different brain areas.

Human brain functional network changes associated with enhanced and impaired attentional task performance.

How is the cognitive performance of the human brain related to its topological and spatial organization as a complex network embedded in anatomical space? To address this question, we used nicotine replacement and duration of attentionally demanding task performance (time-on-task), as experimental factors expected, respectively, to enhance and impair cognitive function. We measured resting-state fMRI data, performance and brain activation on a go/no-go task demanding sustained attention, and subjective fatigue in n = 18 healthy, briefly abstinent, cigarette smokers scanned repeatedly in a placebo-controlled, crossover design. We tested the main effects of drug (placebo vs Nicorette gum) and time-on-task on behavioral performance and brain functional network metrics measured in binary graphs of 477 regional nodes (efficiency, measure of integrative topology; clustering, a measure of segregated topology; and the Euclidean physical distance between connected nodes, a proxy marker of wiring cost). Nicotine enhanced attentional task performance behaviorally and increased efficiency, decreased clustering, and increased connection distance of brain networks. Greater behavioral benefits of nicotine were correlated with stronger drug effects on integrative and distributed network configuration and with greater frequency of cigarette smoking. Greater time-on-task had opposite effects: it impaired attentional accuracy, decreased efficiency, increased clustering, and decreased connection distance of networks. These results are consistent with hypothetical predictions that superior cognitive performance should be supported by more efficient, integrated (high capacity) brain network topology at greater connection distance (high cost). They also demonstrate that brain network analysis can provide novel and theoretically principled pharmacodynamic biomarkers of pro-cognitive drug effects in humans.