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 the neurodevelopmental predictors of psychopathology.

Neuroimaging research has uncovered a multitude of neural abnormalities associated with psychopathology, but few prediction-based studies have been conducted during adolescence, and even fewer used neurobiological features that were extracted across multiple neuroimaging modalities. This gap in the literature is critical, as deriving accurate brain-based models of psychopathology is an essential step towards understanding key neural mechanisms and identifying high-risk individuals. As such, we trained adaptive tree-boosting algorithms on multimodal neuroimaging features from the Lifespan Human Connectome Developmental (HCP-D) sample that contained 956 participants between the ages of 8 to 22 years old. Our feature space consisted of 1037 anatomical, 1090 functional, and 192 diffusion MRI features, which were used to derive models that separately predicted internalizing symptoms, externalizing symptoms, and the general psychopathology factor. We found that multimodal models were the most accurate, but all brain-based models of psychopathology yielded out-of-sample predictions that were weakly correlated with actual symptoms (r2 < 0.15). White matter microstructural properties, including orientation dispersion indices and intracellular volume fractions, were the most predictive of general psychopathology, followed by cortical thickness and functional connectivity. Spatially, the most predictive features of general psychopathology were primarily localized within the default mode and dorsal attention networks. These results were mostly consistent across all dimensions of psychopathology, except orientation dispersion indices and the default mode network were not as heavily weighted in the prediction of internalizing and externalizing symptoms. Taken with prior literature, it appears that neurobiological features are an important part of the equation for predicting psychopathology but relying exclusively on neural markers is clearly not sufficient, especially among adolescent samples with subclinical symptoms. Consequently, risk factor models of psychopathology may benefit from incorporating additional sources of information that have also been shown to explain individual differences, such as psychosocial factors, environmental stressors, and genetic vulnerabilities.

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.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

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.

Caregiver monitoring, but not caregiver warmth, is associated with general cognition in two large sub-samples of youth.

Individual differences in cognitive abilities emerge early during development, and children with poorer cognition are at increased risk for adverse outcomes as they enter adolescence. Caregiving plays an important role in supporting cognitive development, yet it remains unclear how specific types of caregiving behaviors may shape cognition, highlighting the need for large-scale studies. In the present study, we characterized replicable yet specific associations between caregiving behaviors and cognition in two large sub-samples of children ages 9-10 years old from the Adolescent Brain Cognitive Development Study® (ABCD). Across both discovery and replication sub-samples, we found that child reports of caregiver monitoring (supervision or regular knowledge of the child's whereabouts) were positively associated with general cognition abilities, after covarying for age, sex, household income, neighborhood deprivation, and parental education. This association was specific to the type of caregiving behavior (caregiver monitoring, but not caregiver warmth), and was most strongly associated with a broad domain of general cognition (but not executive function or learning/memory). Additionally, we found that caregiver monitoring partially mediated the association between household income and cognition, furthering our understanding of how socioeconomic disparities may contribute to disadvantages in cognitive development. Together, these findings underscore the influence of differences in caregiving behavior in shaping youth cognition. RESEARCH HIGHLIGHTS: Caregiver monitoring, but not caregiver warmth, is associated with cognitive performance in youth Caregiver monitoring partially mediates the association between household income and cognition Results replicated across two large matched samples from the Adolescent Brain Cognitive Development Study® (ABCD).

Characterizing the role of the structural connectome in seizure dynamics.

How does the human brain's structural scaffold give rise to its intricate functional dynamics? This is a central question in translational neuroscience that is particularly relevant to epilepsy, a disorder affecting over 50 million subjects worldwide. Treatment for medication-resistant focal epilepsy is often structural-through surgery or laser ablation-but structural targets, particularly in patients without clear lesions, are largely based on functional mapping via intracranial EEG. Unfortunately, the relationship between structural and functional connectivity in the seizing brain is poorly understood. In this study, we quantify structure-function coupling, specifically between white matter connections and intracranial EEG, across pre-ictal and ictal periods in 45 seizures from nine patients with unilateral drug-resistant focal epilepsy. We use high angular resolution diffusion imaging (HARDI) tractography to construct structural connectivity networks and correlate these networks with time-varying broadband and frequency-specific functional networks derived from coregistered intracranial EEG. Across all frequency bands, we find significant increases in structure-function coupling from pre-ictal to ictal periods. We demonstrate that short-range structural connections are primarily responsible for this increase in coupling. Finally, we find that spatiotemporal patterns of structure-function coupling are highly stereotyped for each patient. These results suggest that seizures harness the underlying structural connectome as they propagate. Mapping the relationship between structural and functional connectivity in epilepsy may inform new therapies to halt seizure spread, and pave the way for targeted patient-specific interventions.

Psychiatric Symptoms, Cognition, and Symptom Severity in Children.

Importance: Mental illnesses are a leading cause of disability globally, and functional disability is often in part caused by cognitive impairments across psychiatric disorders. However, studies have consistently reported seemingly opposite findings regarding the association between cognition and psychiatric symptoms. Objective: To determine if the association between general cognition and mental health symptoms diverges at different symptom severities in children. Design, Setting, and Participants: A total of 5175 children with complete data at 2 time points assessed 2 years apart (aged 9 to 11 years at the first assessment) from the ongoing Adolescent Brain and Cognitive Development (ABCD) study were evaluated for a general cognition factor and mental health symptoms from September 2016 to August 2020 at 21 sites across the US. Polynomial and generalized additive models afforded derivation of continuous associations between cognition and psychiatric symptoms across different ranges of symptom severity. Data were analyzed from December 2022 to April 2024. Main Outcomes and Measures: Aggregate cognitive test scores (general cognition) were primarily evaluated in relation to total and subscale-specific symptoms reported from the Child Behavioral Checklist. Results: The sample included 5175 children (2713 male [52.4%] and 2462 female [47.6%]; mean [SD] age, 10.9 [1.18] years). Previously reported mixed findings regarding the association between general cognition and symptoms may consist of several underlying, opposed associations that depend on the class and severity of symptoms. Linear models recovered differing associations between general cognition and mental health symptoms, depending on the range of symptom severities queried. Nonlinear models confirm that internalizing symptoms were significantly positively associated with cognition at low symptom burdens higher cognition = more symptoms) and significantly negatively associated with cognition at high symptom burdens. Conclusions and Relevance: The association between mental health symptoms and general cognition in this study was nonlinear. Internalizing symptoms were both positively and negatively associated with general cognition at a significant level, depending on the range of symptom severities queried in the analysis sample. These results appear to reconcile mixed findings in prior studies, which implicitly assume that symptom severity tracks linearly with cognitive ability across the entire spectrum of mental health. As the association between cognition and symptoms may be opposite in low vs high symptom severity samples, these results reveal the necessity of clinical enrichment in studies of cognitive impairment.

Adaptive cognitive control circuit changes associated with problem-solving ability and depression symptom outcomes over 24 months.

Mechanistically targeted behavioral interventions are a much-needed strategy for improving outcomes in depression, especially for vulnerable populations with comorbidities such as obesity. Such interventions may change behavior and outcome by changing underlying neural circuit function. However, it is unknown how these circuit-level modifications unfold over intervention and how individual differences in early circuit-level modifications may explain the heterogeneity of treatment effects. We addressed this need within a clinical trial of problem-solving therapy for participants with depression symptoms and comorbid obesity, focusing on the cognitive control circuit as a putative neural mechanism of action. Functional magnetic resonance imaging was applied to measure the cognitive control circuit activity at five time points over 24 months. Compared with participants who received usual care, those receiving problem-solving therapy showed that attenuations in cognitive control circuit activity were associated with enhanced problem-solving ability, which suggests that this circuit plays a key role in the mechanisms of problem-solving therapy. Attenuations in circuit activity were also associated with improved depression symptoms. Changes in cognitive control circuit activity at 2 months better predicted changes in problem-solving ability and depression symptoms at 6, 12, and 24 months, with predictive improvements ranging from 17.8 to 104.0%, exceeding baseline demographic and symptom characteristics. Our findings suggest that targeting the circuit mechanism of action could enhance the prediction of treatment outcomes, warranting future model refinement and improvement to pave the way for its clinical application.

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.

Personalized brain circuit scores identify clinically distinct biotypes in depression and anxiety.

There is an urgent need to derive quantitative measures based on coherent neurobiological dysfunctions or 'biotypes' to enable stratification of patients with depression and anxiety. We used task-free and task-evoked data from a standardized functional magnetic resonance imaging protocol conducted across multiple studies in patients with depression and anxiety when treatment free (n = 801) and after randomization to pharmacotherapy or behavioral therapy (n = 250). From these patients, we derived personalized and interpretable scores of brain circuit dysfunction grounded in a theoretical taxonomy. Participants were subdivided into six biotypes defined by distinct profiles of intrinsic task-free functional connectivity within the default mode, salience and frontoparietal attention circuits, and of activation and connectivity within frontal and subcortical regions elicited by emotional and cognitive tasks. The six biotypes showed consistency with our theoretical taxonomy and were distinguished by symptoms, behavioral performance on general and emotional cognitive computerized tests, and response to pharmacotherapy as well as behavioral therapy. Our results provide a new, theory-driven, clinically validated and interpretable quantitative method to parse the biological heterogeneity of depression and anxiety. Thus, they represent a promising approach to advance precision clinical care in psychiatry.

Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy.

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

Reproducible Sex Differences in Personalized Functional Network Topography in Youth.

Background: A key step towards understanding psychiatric disorders that disproportionately impact female mental health is delineating the emergence of sex-specific patterns of brain organization at the critical transition from childhood to adolescence. Prior work suggests that individual differences in the spatial organization of functional brain networks across the cortex are associated with psychopathology and differ systematically by sex. Aims: We aimed to evaluate the impact of sex on the spatial organization of person-specific functional brain networks. Method: We leveraged person-specific atlases of functional brain networks defined using nonnegative matrix factorization in a sample of n = 6437 youths from the Adolescent Brain Cognitive Development Study. Across independent discovery and replication samples, we used generalized additive models to uncover associations between sex and the spatial layout ("topography") of personalized functional networks (PFNs). Next, we trained support vector machines to classify participants' sex from multivariate patterns of PFN topography. Finally, we leveraged transcriptomic data from the Allen Human Brain Atlas to evaluate spatial correlations between sex differences in PFN topography and gene expression. Results: Sex differences in PFN topography were greatest in association networks including the fronto-parietal, ventral attention, and default mode networks. Machine learning models trained on participants' PFNs were able to classify participant sex with high accuracy. Brain regions with the greatest sex differences in PFN topography were enriched in expression of X-linked genes as well as genes expressed in astrocytes and excitatory neurons. Conclusions: Sex differences in PFN topography are robust, replicate across large-scale samples of youth, and are associated with expression patterns of X-linked genes. These results suggest a potential contributor to the female-biased risk in depressive and anxiety disorders that emerge at the transition from childhood to adolescence.

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.

Hierarchical functional system development supports executive function.

In this perspective, we describe how developmental improvements in youth executive function (EF) are supported by hierarchically organized maturational changes in functional brain systems. We first highlight evidence that functional brain systems are embedded within a hierarchical sensorimotor-association axis of cortical organization. We then review data showing that functional system developmental profiles vary along this axis: systems near the associative end become more functionally segregated, while those in the middle become more integrative. Developmental changes that strengthen the hierarchical organization of the cortex may support EF by facilitating top-down information flow and balancing within- and between-system communication. We propose a central role for attention and frontoparietal control systems in the maturation of healthy EF and suggest that reduced functional system differentiation across the sensorimotor-association axis contributes to transdiagnostic EF deficits.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges.

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

Functional imaging studies of acute administration of classic psychedelics, ketamine, and MDMA: Methodological limitations and convergent results.

Functional magnetic resonance imaging (fMRI) is increasingly used to non-invasively study the acute impact of psychedelics on the human brain. While fMRI is a promising tool for measuring brain function in response to psychedelics, it also has known methodological challenges. We conducted a systematic review of fMRI studies examining acute responses to experimentally administered psychedelics in order to identify convergent findings and characterize heterogeneity in the literature. We reviewed 91 full-text papers; these studies were notable for substantial heterogeneity in design, task, dosage, drug timing, and statistical approach. Data recycling was common, with 51 unique samples across 91 studies. Fifty-seven studies (54%) did not meet contemporary standards for Type I error correction or control of motion artifact. Psilocybin and LSD were consistently reported to moderate the connectivity architecture of the sensorimotor-association cortical axis. Studies also consistently reported that ketamine administration increased activation in the dorsomedial prefrontal cortex. Moving forward, use of best practices such as pre-registration, standardized image processing and statistical testing, and data sharing will be important in this rapidly developing field.

Individual Differences in Delay Discounting are Associated with Dorsal Prefrontal Cortex Connectivity in Youth.

Delay discounting is a measure of impulsive choice relevant in adolescence as it predicts many real-life outcomes, including substance use disorders, obesity, and academic achievement. However, the functional networks underlying individual differences in delay discounting during youth remain incompletely described. Here we investigate the association between multivariate patterns of functional connectivity and individual differences in impulsive choice in a large sample of youth. A total of 293 youth (9-23 years) completed a delay discounting task and underwent resting-state fMRI at 3T. A connectome-wide analysis using multivariate distance-based matrix regression was used to examine whole-brain relationships between delay discounting and functional connectivity was then performed. These analyses revealed that individual differences in delay discounting were associated with patterns of connectivity emanating from the left dorsal prefrontal cortex, a hub of the default mode network. Delay discounting was associated with greater functional connectivity between the dorsal prefrontal cortex and other parts of the default mode network, and reduced connectivity with regions in the dorsal and ventral attention networks. These results suggest that delay discounting in youth is associated with individual differences in relationships both within the default mode network and between the default mode and networks involved in attentional and cognitive control.

Individual differences in delay discounting are associated with dorsal prefrontal cortex connectivity in children, adolescents, and adults.

Delay discounting is a measure of impulsive choice relevant in adolescence as it predicts many real-life outcomes, including obesity and academic achievement. However, resting-state functional networks underlying individual differences in delay discounting during youth remain incompletely described. Here we investigate the association between multivariate patterns of functional connectivity and individual differences in impulsive choice in a large sample of children, adolescents, and adults. A total of 293 participants (9-23 years) completed a delay discounting task and underwent 3T resting-state fMRI. A connectome-wide analysis using multivariate distance-based matrix regression was used to examine whole-brain relationships between delay discounting and functional connectivity. These analyses revealed that individual differences in delay discounting were associated with patterns of connectivity emanating from the left dorsal prefrontal cortex, a default mode network hub. Greater delay discounting was associated with greater functional connectivity between the dorsal prefrontal cortex and other default mode network regions, but reduced connectivity with regions in the dorsal and ventral attention networks. These results suggest delay discounting in children, adolescents, and adults is associated with individual differences in relationships both within the default mode network and between the default mode and networks involved in attentional and cognitive control.

Development of top-down cortical propagations in youth.

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.