The Lifespan Human Connectome Project in Development: A large-scale study of brain connectivity development in 5-21 year olds

Structural deviations of the posterior fossa and the cerebellum and their cognitive links in a neurodevelopmental deletion syndrome

High-impact genetic variants associated with neurodevelopmental disorders provide biologically-defined entry points for mechanistic investigation. The 3q29 deletion (3q29Del) is one such variant, conferring a 40-100-fold increased risk for schizophrenia, as well as high risk for autism and intellectual disability. However, the mechanisms leading to neurodevelopmental disability remain largely unknown. Here, we report the first in vivo quantitative neuroimaging study in individuals with 3q29Del (N = 24) and neurotypical controls (N = 1608) using structural MRI. Given prior radiology reports of posterior fossa abnormalities in 3q29Del, we focused our investigation on the cerebellum and its tissue-types and lobules. Additionally, we compared the prevalence of cystic/cyst-like malformations of the posterior fossa between 3q29Del and controls and examined the association between neuroanatomical findings and quantitative traits to probe gene-brain-behavior relationships. 3q29Del participants had smaller cerebellar cortex volumes than controls, before and after correction for intracranial volume (ICV). An anterior-posterior gradient emerged in finer grained lobule-based and voxel-wise analyses. 3q29Del participants also had larger cerebellar white matter volumes than controls following ICV-correction and displayed elevated rates of posterior fossa arachnoid cysts and mega cisterna magna findings independent of cerebellar volume. Cerebellar white matter and subregional gray matter volumes were associated with visual-perception and visual-motor integration skills as well as IQ, while cystic/cyst-like malformations yielded no behavioral link. In summary, we find that abnormal development of cerebellar structures may represent neuroimaging-based biomarkers of cognitive and sensorimotor function in 3q29Del, adding to the growing evidence identifying cerebellar pathology as an intersection point between syndromic and idiopathic forms of neurodevelopmental disabilities.

Longitudinal hippocampal axis in large-scale cortical systems underlying development and episodic memory

The hippocampus is functionally specialized along its longitudinal axis with intricate interactions with cortical systems, which is crucial for understanding development and cognition. Using a well-established connectopic mapping technique on two large resting-state functional MRI datasets, we systematically quantified topographic organization of the hippocampal functional connectivity (hippocampal gradient) and its cortical interaction in developing brains. We revealed hippocampal functional hierarchy within the large-scale cortical brain systems, with the anterior hippocampus preferentially connected to an anterior temporal (AT) pathway and the posterior hippocampus embedded in a posterior medial (PM) pathway. We examined the developmental effects of the primary gradient and its whole-brain functional interaction. We observed increased functional specialization along the hippocampal long axis and found a general whole-brain connectivity shift from the posterior to the anterior hippocampus during development. Using phenotypic predictive modeling, we further delineated how the hippocampus is differentially integrated into the whole-brain cortical hierarchy underlying episodic memory and identified several key nodes within PM/AT systems. Our results highlight the importance of hippocampal gradient and its cortical interaction in development and for supporting episodic memory.

A review of artificial intelligence-based brain age estimation and its applications for related diseases

The study of brain age has emerged over the past decade, aiming to estimate a person's age based on brain imaging scans. Ideally, predicted brain age should match chronological age in healthy individuals. However, brain structure and function change in the presence of brain-related diseases. Consequently, brain age also changes in affected individuals, making the brain age gap (BAG)-the difference between brain age and chronological age-a potential biomarker for brain health, early screening, and identifying age-related cognitive decline and disorders. With the recent successes of artificial intelligence in healthcare, it is essential to track the latest advancements and highlight promising directions. This review paper presents recent machine learning techniques used in brain age estimation (BAE) studies. Typically, BAE models involve developing a machine learning regression model to capture age-related variations in brain structure from imaging scans of healthy individuals and automatically predict brain age for new subjects. The process also involves estimating BAG as a measure of brain health. While we discuss recent clinical applications of BAE methods, we also review studies of biological age that can be integrated into BAE research. Finally, we point out the current limitations of BAE's studies.

Brain-phenotype predictions of language and executive function can survive across diverse real-world data: Dataset shifts in developmental populations

Predictive modeling potentially increases the reproducibility and generalizability of neuroimaging brain-phenotype associations. Yet, the evaluation of a model in another dataset is underutilized. Among studies that undertake external validation, there is a notable lack of attention to generalization across dataset-specific idiosyncrasies (i.e., dataset shifts). Research settings, by design, remove the between-site variations that real-world and, eventually, clinical applications demand. Here, we rigorously test the ability of a range of predictive models to generalize across three diverse, unharmonized developmental samples: the Philadelphia Neurodevelopmental Cohort (n=1291), the Healthy Brain Network (n=1110), and the Human Connectome Project in Development (n=428). These datasets have high inter-dataset heterogeneity, encompassing substantial variations in age distribution, sex, racial and ethnic minority representation, recruitment geography, clinical symptom burdens, fMRI tasks, sequences, and behavioral measures. Through advanced methodological approaches, we demonstrate that reproducible and generalizable brain-behavior associations can be realized across diverse dataset features. Results indicate the potential of functional connectome-based predictive models to be robust despite substantial inter-dataset variability. Notably, for the HCPD and HBN datasets, the best predictions were not from training and testing in the same dataset (i.e., cross-validation) but across datasets. This result suggests that training on diverse data may improve prediction in specific cases. Overall, this work provides a critical foundation for future work evaluating the generalizability of brain-phenotype associations in real-world scenarios and clinical settings.

Release the Krakencoder: A unified brain connectome translation and fusion tool

Brain connectivity can be estimated in many ways, depending on modality and processing strategy. Here we present the Krakencoder, a joint connectome mapping tool that simultaneously, bidirectionally translates between structural (SC) and functional connectivity (FC), and across different atlases and processing choices via a common latent representation. These mappings demonstrate unprecedented accuracy and individual-level identifiability; the mapping between SC and FC has identifiability 42-54% higher than existing models. The Krakencoder combines all connectome flavors via a shared low-dimensional latent space. This "fusion" representation i) better reflects familial relatedness, ii) preserves age- and sex-relevant information and iii) enhances cognition-relevant information. The Krakencoder can be applied without retraining to new, out-of-age-distribution data while still preserving inter-individual differences in the connectome predictions and familial relationships in the latent representations. The Krakencoder is a significant leap forward in capturing the relationship between multi-modal brain connectomes in an individualized, behaviorally- and demographically-relevant way.

Development of Higher-Level Vision: A Network Perspective

Most studies on the development of the visual system have focused on the mechanisms shaping early visual stages up to the level of primary visual cortex (V1). Much less is known about the development of the stages after V1 that handle the higher visual functions fundamental to everyday life. The standard model for the maturation of these areas is that it occurs sequentially, according to the positions of areas in the adult hierarchy. Yet, the existing literature reviewed here paints a different picture, one in which the adult configuration emerges through a sequence of unique network configurations that are not mere partial versions of the adult hierarchy. In addition to studying higher visual development per se to fill major gaps in knowledge, it will be crucial to adopt a network-level perspective in future investigations to unravel normal developmental mechanisms, identify vulnerabilities to developmental disorders, and eventually devise treatments for these disorders.

Comparing single- and multi-post labeling delays for the measurements of resting cerebral and hippocampal blood flow for cerebrovascular testing in midlife adults

Objectives

To assess the reliability and validity of measuring resting cerebral blood flow (CBF) and hippocampal CBF using a single-post-labeling delay (PLD) and a multi-PLD pseudo-continuous arterial spin labeling (pCASL) protocol for cerebrovascular reactivity (CVR) testing.

Methods

25 healthy, midlife adults (57 ± 4 years old) were imaged in a Siemens Prisma 3T magnetic resonance imaging (MRI) scanner. Resting CBF and hippocampal CBF were assessed using two pCASL protocols, our modified single-PLD protocol (pCASL-MOD) to accommodate the needs for CVR testing and the multi-PLD Human Connectome Project (HCP) Lifespan protocol to serve as the reference control (pCASL-HCP). During pCASL-MOD, CVR was calculated as the change in CBF from rest to hypercapnia (+9 mmHg increase in end-tidal partial pressure of carbon dioxide [PETCO2]) and then normalized for PETCO2. The reliability and validity in resting gray matter (GM) CBF, white matter (WM) CBF, and hippocampal CBF between pCASL-MOD and pCASL-HCP protocols were examined using correlation analyses, paired t-tests, and Bland Altman plots.

Results

The pCASL-MOD and pCASL-HCP protocols were significantly correlated for resting GM CBF [r = 0.72; F (1, 23) = 25.24, p < 0.0001], WM CBF [r = 0.57; F (1, 23) = 10.83, p = 0.003], and hippocampal CBF [r = 0.77; F (1, 23) = 32.65, p < 0.0001]. However, pCASL-MOD underestimated resting GM CBF (pCASL-MOD: 53.7 ± 11.1 v. pCASL-HCP: 69.1 ± 13.1 mL/100 g/min; p < 0.0001), WM CBF (pCASL-MOD: 32.4 ± 4.8 v. pCASL-HCP: 35.5 ± 6.9 mL/100 g/min; p = 0.01), and hippocampal CBF (pCASL-MOD: 50.5 ± 9.0 v. pCASL-HCP: 68.1 ± 12.5 mL/100 g/min; p < 0.0001). PETCO2 increased by 8.0 ± 0.7 mmHg to induce CVR (GM CBF: 4.8% ± 2.6%; WM CBF 2.9% ± 2.5%; and hippocampal CBF: 3.4% ± 3.8%).

Conclusion

Our single-PLD pCASL-MOD protocol reliably measured CBF and hippocampal CBF at rest given the significant correlation with the multi-PLD pCASL-HCP protocol. Despite the lower magnitude relative to pCASL-HCP, we recommend using our pCASL-MOD protocol for CVR testing in which an exact estimate of CBF is not required such as the assessment of relative change in CBF to hypercapnia.

Power and reproducibility in the external validation of brain-phenotype predictions

Brain-phenotype predictive models seek to identify reproducible and generalizable brain-phenotype associations. External validation, or the evaluation of a model in external datasets, is the gold standard in evaluating the generalizability of models in neuroimaging. Unlike typical studies, external validation involves two sample sizes: the training and the external sample sizes. Thus, traditional power calculations may not be appropriate. Here we ran over 900 million resampling-based simulations in functional and structural connectivity data to investigate the relationship between training sample size, external sample size, phenotype effect size, theoretical power and simulated power. Our analysis included a wide range of datasets: the Healthy Brain Network, the Adolescent Brain Cognitive Development Study, the Human Connectome Project (Development and Young Adult), the Philadelphia Neurodevelopmental Cohort, the Queensland Twin Adolescent Brain Project, and the Chinese Human Connectome Project; and phenotypes: age, body mass index, matrix reasoning, working memory, attention problems, anxiety/depression symptoms and relational processing. High effect size predictions achieved adequate power with training and external sample sizes of a few hundred individuals, whereas low and medium effect size predictions required hundreds to thousands of training and external samples. In addition, most previous external validation studies used sample sizes prone to low power, and theoretical power curves should be adjusted for the training sample size. Furthermore, model performance in internal validation often informed subsequent external validation performance (Pearson's r difference <0.2), particularly for well-harmonized datasets. These results could help decide how to power future external validation studies.

FunMaps: a method for parcellating functional brain networks using resting-state functional MRI data

Parcellations of resting-state functional magnetic resonance imaging (rs-fMRI) data are widely used to create topographical maps of functional networks in the human brain. While such network maps are highly useful for studying brain organization and function, they usually require large sample sizes to make them, thus creating practical limitations for researchers that would like to carry out parcellations on data collected in their labs. Furthermore, it can be difficult to quantitatively evaluate the results of a parcellation since networks are usually identified using a clustering algorithm, like principal components analysis, on the results of a single group-averaged connectivity map. To address these challenges, we developed the FunMaps method: a parcellation routine that intrinsically incorporates stability and replicability of the parcellation by keeping only network distinctions that agree across halves of the data over multiple random iterations. Here, we demonstrate the efficacy and flexibility of FunMaps, while describing step-by-step instructions for running the program. The FunMaps method is publicly available on GitHub (https://github.com/persichetti-lab/FunMaps). It includes source code for running the parcellation and auxiliary code for preparing data, evaluating the parcellation, and displaying the results.

DeepResBat: Deep residual batch harmonization accounting for covariate distribution differences

Pooling MRI data from multiple datasets requires harmonization to reduce undesired inter-site variabilities, while preserving effects of biological variables (or covariates). The popular harmonization approach ComBat uses a mixed effect regression framework that explicitly accounts for covariate distribution differences across datasets. There is also significant interest in developing harmonization approaches based on deep neural networks (DNNs), such as conditional variational autoencoder (cVAE). However, current DNN approaches do not explicitly account for covariate distribution differences across datasets. Here, we provide mathematical results, suggesting that not accounting for covariates can lead to suboptimal harmonization. We propose two DNN-based covariate-aware harmonization approaches: covariate VAE (coVAE) and DeepResBat. The coVAE approach is a natural extension of cVAE by concatenating covariates and site information with site- and covariate-invariant latent representations. DeepResBat adopts a residual framework inspired by ComBat. DeepResBat first removes the effects of covariates with nonlinear regression trees, followed by eliminating site differences with cVAE. Finally, covariate effects are added back to the harmonized residuals. Using three datasets from three continents with a total of 2787 participants and 10,085 anatomical T1 scans, we find that DeepResBat and coVAE outperformed ComBat, CovBat and cVAE in terms of removing dataset differences, while enhancing biological effects of interest. However, coVAE hallucinates spurious associations between anatomical MRI and covariates even when no association exists. Future studies proposing DNN-based harmonization approaches should be aware of this false positive pitfall. Overall, our results suggest that DeepResBat is an effective deep learning alternative to ComBat. Code for DeepResBat can be found here: https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/harmonization/An2024_DeepResBat.

Quantifying brain development in the HEALthy Brain and Child Development (HBCD) Study: The magnetic resonance imaging and spectroscopy protocol

The HEALthy Brain and Child Development (HBCD) Study, a multi-site prospective longitudinal cohort study, will examine human brain, cognitive, behavioral, social, and emotional development beginning prenatally and planned through early childhood. The acquisition of multimodal magnetic resonance-based brain development data is central to the study's core protocol. However, application of Magnetic Resonance Imaging (MRI) methods in this population is complicated by technical challenges and difficulties of imaging in early life. Overcoming these challenges requires an innovative and harmonized approach, combining age-appropriate acquisition protocols together with specialized pediatric neuroimaging strategies. The HBCD MRI Working Group aimed to establish a core acquisition protocol for all 27 HBCD Study recruitment sites to measure brain structure, function, microstructure, and metabolites. Acquisition parameters of individual modalities have been matched across MRI scanner platforms for harmonized acquisitions and state-of-the-art technologies are employed to enable faster and motion-robust imaging. Here, we provide an overview of the HBCD MRI protocol, including decisions of individual modalities and preliminary data. The result will be an unparalleled resource for examining early neurodevelopment which enables the larger scientific community to assess normative trajectories from birth through childhood and to examine the genetic, biological, and environmental factors that help shape the developing brain.

Group-common and individual-specific effects of structure-function coupling in human brain networks with graph neural networks

The human cerebral cortex is organized into functionally segregated but synchronized regions bridged by the structural connectivity of white matter pathways. While structure-function coupling has been implicated in cognitive development and neuropsychiatric disorders, it remains unclear to what extent the structure-function coupling reflects a group-common characteristic or varies across individuals, at both the global and regional brain levels. By leveraging two independent, high-quality datasets, we found that the graph neural network accurately predicted unseen individuals' functional connectivity from structural connectivity, reflecting a strong structure-function coupling. This coupling was primarily driven by network topology and was substantially stronger than that of the correlation approaches. Moreover, we observed that structure-function coupling was dominated by group-common effects, with subtle yet significant individual-specific effects. The regional group and individual effects of coupling were hierarchically organized across the cortex along a sensorimotor-association axis, with lower group and higher individual effects in association cortices. These findings emphasize the importance of considering both group and individual effects in understanding cortical structure-function coupling, suggesting insights into interpreting individual differences of the coupling and informing connectivity-guided therapeutics.

COVID-19 lockdown effects on adolescent brain structure suggest accelerated maturation that is more pronounced in females than in males

Adolescence is a period of substantial social-emotional development, accompanied by dramatic changes to brain structure and function. Social isolation due to lockdowns that were imposed because of the COVID-19 pandemic had a detrimental impact on adolescent mental health, with the mental health of females more affected than males. We assessed the impact of the COVID-19 pandemic lockdowns on adolescent brain structure with a focus on sex differences. We collected MRI structural data longitudinally from adolescents prior to and after the pandemic lockdowns. The pre-COVID data were used to create a normative model of cortical thickness change with age during typical adolescent development. Cortical thickness values in the post-COVID data were compared to this normative model. The analysis revealed accelerated cortical thinning in the post-COVID brain, which was more widespread throughout the brain and greater in magnitude in females than in males. When measured in terms of equivalent years of development, the mean acceleration was found to be 4.2 y in females and 1.4 y in males. Accelerated brain maturation as a result of chronic stress or adversity during development has been well documented. These findings suggest that the lifestyle disruptions associated with the COVID-19 pandemic lockdowns caused changes in brain biology and had a more severe impact on the female than the male brain.

Testing a Conceptual Model of Early Adversity, Neural Function, and Psychopathology: Protocol for a Retrospective Observational Cohort Study

BACKGROUND

Early adversity, broadly defined as a set of negative exposures during childhood, is extremely common and increases risk for psychopathology across the life span. Previous research suggests that separate dimensions of adversity increase risk through developmental plasticity mechanisms shaping unique neurobiological pathways. Specifically, research suggests that deprivation is associated with deficits in higher order cognition, while threat is associated with atypicality in fear learning and emotion dysregulation. However, most of this research has been conducted in adolescent and adult samples, long after exposure to adversity occurs and far from periods of peak developmental plasticity.

OBJECTIVE

The Wellness Health and Life Experiences (WHALE) study examines the neurobiological and behavioral mechanisms by which deprivation, threat, and unpredictability increase risk for psychopathology in early childhood (age 4-7 years) directly following periods of peak developmental plasticity. The objective of this study is to describe the study rationale and aims, the research design and procedures, and the analytical plan to test the study hypotheses.

METHODS

This is a retrospective cohort study that examines associations between exposure to deprivation and threat and their hypothesized neurobiological mechanisms, how these neurobiological mechanisms link early adversity and psychopathology, and associations between unpredictability, reward learning, and psychopathology. The sample was a convenience sample of children (aged 4-7 years) and their families, identified through flyers, email blasts to listserves, school-based advertising, and involvement in community events. Data were collected during a home visit, a subsequent laboratory visit, and a final neuroimaging visit. Planned analyses include linear regression, path analyses, and functional magnetic resonance imaging analyses to explore the role of neural function in the association between early adversity and psychopathology.

RESULTS

Participants (N=301) have been recruited into the study, and data collection has commenced. The expected results will be available in 2024.

CONCLUSIONS

The findings of this study will help elucidate the neurobiological mechanisms by which early adversity increases risk for psychopathology in early childhood. This study represents the earliest test of an influential theory of biological embedding of early adversity.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/59636.

Generalizable Links Between Borderline Personality Traits and Functional Connectivity

BACKGROUND

Symptoms of borderline personality disorder (BPD) often manifest during adolescence, but the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here, we aimed to investigate how multivariate patterns of functional connectivity are associated with borderline personality traits in large samples of young adults and adolescents.

METHODS

We used functional magnetic resonance imaging data from young adults and adolescents from the HCP-YA (Human Connectome Project Young Adult) (n = 870, ages 22-37 years, 457 female) and the HCP-D (Human Connectome Project Development) (n = 223, ages 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five-Factor Inventory. A ridge regression model with cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity.

RESULTS

Multivariate functional connectivity patterns significantly predicted out-of-sample BPD scores in unseen data in young adults (HCP-YA ppermuted = .001) and older adolescents (HCP-D ppermuted = .001). Regional predictive capacity was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD scores aligned with those associated with development in youth.

CONCLUSIONS

Individual differences in functional connectivity in developmentally sensitive regions are associated with borderline personality traits.

Lifespan development of thalamic nuclei and characterizing thalamic nuclei abnormalities in schizophrenia using normative modeling

Thalamic abnormalities have been repeatedly implicated in the pathophysiology of schizophrenia and other neurodevelopmental disorders. Uncovering the etiology of thalamic abnormalities and how they may contribute to illness phenotypes faces at least two obstacles. First, the typical developmental trajectories of thalamic nuclei and their association with cognition across the lifespan are largely unknown. Second, modest effect sizes indicate marked individual differences and pose a significant challenge to personalized medicine. To address these knowledge gaps, we characterized the development of thalamic nuclei volumes using normative models generated from the Human Connectome Project Lifespan datasets (5-100+ years), then applied them to an independent clinical cohort to determine the frequency of thalamic volume deviations in people with schizophrenia (17-61 years). Normative models revealed diverse non-linear age effects across the lifespan. Association nuclei exhibited negative age effects during youth but stabilized in adulthood until turning negative again with older age. Sensorimotor nuclei volumes remained relatively stable through youth and adulthood until also turning negative with older age. Up to 18% of individuals with schizophrenia exhibited abnormally small (i.e., below the 5th centile) mediodorsal and pulvinar volumes, and the degree of deviation, but not raw volumes, correlated with the severity of cognitive impairment. While case-control differences are robust, only a minority of patients demonstrate unusually small thalamic nuclei volumes. Normative modeling enables the identification of these individuals, which is a necessary step toward precision medicine.

Unveiling the core functional networks of cognition: An ontology-guided machine learning approach

Deciphering the functional architecture that underpins diverse cognitive functions is fundamental quest in neuroscience. In this study, we employed an innovative machine learning framework that integrated cognitive ontology with functional connectivity analysis to identify brain networks essential for cognition. We identified a core assembly of functional connectomes, primarily located within the association cortex, which showed superior predictive performance compared to two conventional methods widely employed in previous research across various cognitive domains. Our approach achieved a mean prediction accuracy of 0.13 across 16 cognitive tasks, including working memory, reading comprehension, and sustained attention, outperforming the traditional methods' accuracy of 0.08. In contrast, our method showed limited predictive power for sensory, motor, and emotional functions, with a mean prediction accuracy of 0.03 across 9 relevant tasks, slightly lower than the traditional methods' accuracy of 0.04. These cognitive connectomes were further characterized by distinctive patterns of resting-state functional connectivity, structural connectivity via white matter tracts, and gene expression, highlighting their neurogenetic underpinnings. Our findings reveal a domain-general functional network fingerprint that pivotal to cognition, offering a novel computational approach to explore the neural foundations of cognitive abilities.

Impact of Stimulation Duration in taVNS-Exploring Multiple Physiological and Cognitive Outcomes

Transcutaneous auricular vagus nerve stimulation (taVNS) is a non-invasive neuromodulation technique that modulates the noradrenergic activity of the locus coeruleus (LC). Yet, there is still uncertainty about the most effective stimulation and reliable outcome parameters. In a double blind, sham-controlled study including a sample of healthy young individuals (N = 29), we compared a shorter (3.4 s) and a longer (30 s) stimulation duration and investigated the effects of taVNS (real vs. sham) on saliva samples (alpha amylase and cortisol concentration), pupil (pupillary light reflex and pupil size at rest) and EEG data (alpha and theta activity at rest, ERPs for No-Go signals), and cognitive tasks (Go/No-Go and Stop Signal Tasks). Salivary alpha amylase concentration was significantly increased in the real as compared to sham stimulation for the 30 s stimulation condition. In the 3.4 s stimulation condition, we found prolonged reaction times and increased error rates in the Go/No-Go task and increased maximum acceleration in the pupillary light reflex. For the other outcomes, no significant differences were found. Our results show that prolonged stimulation increases salivary alpha-amylase, which was expected from the functional properties of the LC. The finding of longer response times to short taVNS stimulation was not expected and cannot be explained by an increase in LC activity. We also discuss the difficulties in assessing pupil size as an expression of taVNS-mediated LC functional changes.

Integrated Brain Connectivity Analysis with fMRI, DTI, and sMRI Powered by Interpretable Graph Neural Networks

Multimodal neuroimaging modeling has become a widely used approach but confronts considerable challenges due to heterogeneity, which encompasses variability in data types, scales, and formats across modalities. This variability necessitates the deployment of advanced computational methods to integrate and interpret these diverse datasets within a cohesive analytical framework. In our research, we amalgamate functional magnetic resonance imaging (fMRI), diffusion tensor imaging (DTI), and structural MRI (sMRI) into a cohesive framework. This integration capitalizes on the unique strengths of each modality and their inherent interconnections, aiming for a comprehensive understanding of the brain's connectivity and anatomical characteristics. Utilizing the Glasser atlas for parcellation, we integrate imaging-derived features from various modalities-functional connectivity from fMRI, structural connectivity from DTI, and anatomical features from sMRI-within consistent regions. Our approach incorporates a masking strategy to differentially weight neural connections, thereby facilitating a holistic amalgamation of multimodal imaging data. This technique enhances interpretability at connectivity level, transcending traditional analyses centered on singular regional attributes. The model is applied to the Human Connectome Project's Development study to elucidate the associations between multimodal imaging and cognitive functions throughout youth. The analysis demonstrates improved predictive accuracy and uncovers crucial anatomical features and essential neural connections, deepening our understanding of brain structure and function. This study not only advances multi-modal neuroimaging analytics by offering a novel method for the integrated analysis of diverse imaging modalities but also improves the understanding of intricate relationship between the brain's structural and functional networks and cognitive development.

Human adolescent brain similarity development is different for paralimbic versus neocortical zones

Adolescent development of human brain structural and functional networks is increasingly recognized as fundamental to emergence of typical and atypical adult cognitive and emotional proodal magnetic resonance imaging (MRI) data collected from N [Formula: see text] 300 healthy adolescents (51%; female; 14 to 26 y) each scanned repeatedly in an accelerated longitudinal design, to provide an analyzable dataset of 469 structural scans and 448 functional MRI scans. We estimated the morphometric similarity between each possible pair of 358 cortical areas on a feature vector comprising six macro- and microstructural MRI metrics, resulting in a morphometric similarity network (MSN) for each scan. Over the course of adolescence, we found that morphometric similarity increased in paralimbic cortical areas, e.g., insula and cingulate cortex, but generally decreased in neocortical areas, and these results were replicated in an independent developmental MRI cohort (N [Formula: see text] 304). Increasing hubness of paralimbic nodes in MSNs was associated with increased strength of coupling between their morphometric similarity and functional connectivity. Decreasing hubness of neocortical nodes in MSNs was associated with reduced strength of structure-function coupling and increasingly diverse functional connections in the corresponding fMRI networks. Neocortical areas became more structurally differentiated and more functionally integrative in a metabolically expensive process linked to cortical thinning and myelination, whereas paralimbic areas specialized for affective and interoceptive functions became less differentiated, as hypothetically predicted by a developmental transition from periallocortical to proisocortical organization of the cortex. Cytoarchitectonically distinct zones of the human cortex undergo distinct neurodevelopmental programs during typical adolescence.

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.

The Human Connectome Project of adolescent anxiety and depression dataset

This article describes primary data and resources available from the Boston Adolescent Neuroimaging of Depression and Anxiety (BANDA) study, a novel arm of the Human Connectome Project (HCP). Data were collected from 215 adolescents (14-17 years old), 152 of whom had current diagnoses of anxiety and/or depressive disorders at study intake. Data include cross-sectional structural (T1- and T2-weighted), functional (resting state and three tasks), and diffusion-weighted magnetic resonance images. Both unprocessed and HCP minimally-preprocessed imaging data are available within the data release packages. Adolescent and parent clinical interview data, as well as cognitive and neuropsychological data are also included within these packages. Release packages additionally provide data collected from self-report measures assessing key features of adolescent psychopathology, including: anxious and depressive symptom dimensions, behavioral inhibition/activation, exposure to stressful life events, and risk behaviors. Finally, the release packages include 6- and 12-month longitudinal data acquired from clinical measures. Data are publicly accessible through the National Institute of Mental Health Data Archive (ID: #2505).

The associations of socioeconomic position with structural brain damage and connectivity and cognitive functioning: The Maastricht Study

BACKGROUND

Socioeconomic inequalities in cognitive impairment may partly act through structural brain damage and reduced connectivity. This study investigated the extent to which the association of early-life socioeconomic position (SEP) with later-life cognitive functioning is mediated by later-life SEP, and whether the associations of SEP with later-life cognitive functioning can be explained by structural brain damage and connectivity.

METHODS

We used cross-sectional data from the Dutch population-based Maastricht Study (n = 4,839; mean age 59.2 ± 8.7 years, 49.8% women). Early-life SEP was assessed by self-reported poverty during childhood and parental education. Later-life SEP included education, occupation, and current household income. Participants underwent cognitive testing and 3-T magnetic resonance imaging to measure volumes of white matter hyperintensities, grey matter, white matter, cerebrospinal fluid, and structural connectivity. Multiple linear regression analyses tested the associations between SEP, markers of structural brain damage and connectivity, and cognitive functioning. Mediation was tested using structural equation modeling.

RESULTS

Although there were direct associations between both indicators of SEP and later-life cognitive functioning, a large part of the association between early-life SEP and later-life cognitive functioning was explained by later-life SEP (72.2%). The extent to which structural brain damage or connectivity acted as mediators between SEP and cognitive functioning was small (up to 5.9%).

CONCLUSIONS

We observed substantial SEP differences in later-life cognitive functioning. Associations of structural brain damage and connectivity with cognitive functioning were relatively small, and only marginally explained the SEP gradients in cognitive functioning.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

Dissecting task-based fMRI activity using normative modelling: an application to the Emotional Face Matching Task

Functional neuroimaging has contributed substantially to understanding brain function but is dominated by group analyses that index only a fraction of the variation in these data. It is increasingly clear that parsing the underlying heterogeneity is crucial to understand individual differences and the impact of different task manipulations. We estimate large-scale (N = 7728) normative models of task-evoked activation during the Emotional Face Matching Task, which enables us to bind heterogeneous datasets to a common reference and dissect heterogeneity underlying group-level analyses. We apply this model to a heterogenous patient cohort, to map individual differences between patients with one or more mental health diagnoses relative to the reference cohort and determine multivariate associations with transdiagnostic symptom domains. For the face>shapes contrast, patients have a higher frequency of extreme deviations which are spatially heterogeneous. In contrast, normative models for faces>baseline have greater predictive value for individuals' transdiagnostic functioning. Taken together, we demonstrate that normative modelling of fMRI task-activation can be used to illustrate the influence of different task choices and map replicable individual differences, and we encourage its application to other neuroimaging tasks in future studies.

Development of human visual cortical function: A scoping review of task- and naturalistic-fMRI studies through the interactive specialization and maturational frameworks

Overarching theories such as the interactive specialization and maturational frameworks have been proposed to describe human functional brain development. However, these frameworks have not yet been systematically examined across the fMRI literature. Visual processing is one of the most well-studied fields in neuroimaging, and research in this area has recently expanded to include naturalistic paradigms that facilitate study in younger age ranges, allowing for an in-depth critical appraisal of these frameworks across childhood. To this end, we conducted a scoping review of 94 developmental visual fMRI studies, including both traditional experimental task and naturalistic studies, across multiple sub-domains (early visual processing, category-specific higher order processing, naturalistic visual processing). We found that across domains, many studies reported progressive development, but few studies describe regressive or emergent changes necessary to fit the maturational or interactive specialization frameworks. Our findings suggest a need for the expansion of developmental frameworks and clearer reporting of both progressive and regressive changes, along with well-powered, longitudinal studies.

Neural underpinnings of individual differences in emotion regulation: A systematic review

This review synthesises individual differences in neural processes related to emotion regulation (ER). It comprises individual differences in self-reported and physiological regulation success, self-reported ER-related traits, and demographic variables, to assess their correlation with brain activation during ER tasks. Considering region-of-interest (ROI) and whole-brain analyses, the review incorporated data from 52 functional magnetic resonance imaging studies. Results can be summarized as follows: (1) Self-reported regulation success (assessed by emotional state ratings after regulation) and self-reported ER-related traits (assessed by questionnaires) correlated with brain activity in the lateral prefrontal cortex. (2) Amygdala activation correlated with ER-related traits only in ROI analyses, while it was associated with regulation success in whole-brain analyses. (3) For demographic and physiological measures, there was no systematic overlap in effects reported across studies. In showing that individual differences in regulation success and ER-related traits can be traced back to differences in the neural activity of brain regions associated with emotional reactivity (amygdala) and cognitive control (lateral prefrontal cortex), our findings can inform prospective personalised intervention models.

Anatomy-aware and acquisition-agnostic joint registration with SynthMorph

Affine image registration is a cornerstone of medical-image analysis. While classical algorithms can achieve excellent accuracy, they solve a time-consuming optimization for every image pair. Deep-learning (DL) methods learn a function that maps an image pair to an output transform. Evaluating the function is fast, but capturing large transforms can be challenging, and networks tend to struggle if a test-image characteristic shifts from the training domain, such as the resolution. Most affine methods are agnostic to the anatomy the user wishes to align, meaning the registration will be inaccurate if algorithms consider all structures in the image. We address these shortcomings with SynthMorph, a fast, symmetric, diffeomorphic, and easy-to-use DL tool for joint affine-deformable registration of any brain image without preprocessing. First, we leverage a strategy that trains networks with widely varying images synthesized from label maps, yielding robust performance across acquisition specifics unseen at training. Second, we optimize the spatial overlap of select anatomical labels. This enables networks to distinguish anatomy of interest from irrelevant structures, removing the need for preprocessing that excludes content which would impinge on anatomy-specific registration. Third, we combine the affine model with a deformable hypernetwork that lets users choose the optimal deformation-field regularity for their specific data, at registration time, in a fraction of the time required by classical methods. This framework is applicable to learning anatomy-aware, acquisition-agnostic registration of any anatomy with any architecture, as long as label maps are available for training. We analyze how competing architectures learn affine transforms and compare state-of-the-art registration tools across an extremely diverse set of neuroimaging data, aiming to truly capture the behavior of methods in the real world. SynthMorph demonstrates high accuracy and is available at https://w3id.org/synthmorph, as a single complete end-to-end solution for registration of brain magnetic resonance imaging (MRI) data.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7-21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Structural connectivity matures along a sensorimotor-association connectional axis in youth

Childhood and adolescence are associated with protracted developmental remodeling of cortico-cortical structural connectivity. However, how heterochronous development in white matter structural connectivity spatially and temporally unfolds across the macroscale human connectome remains unknown. Leveraging non-invasive diffusion MRI data from both cross-sectional (N = 590) and longitudinal (baseline: N = 3,949; two-year follow-up: N = 3,155) developmental datasets, we found that structural connectivity development diverges along a pre-defined sensorimotor-association (S-A) connectional axis from ages 8.1 to 21.9 years. Specifically, we observed a continuum of developmental profiles that spans from an early childhood increase in connectivity strength in sensorimotor-sensorimotor connections to a late adolescent increase in association-association connectional strength. The S-A connectional axis also captured spatial variations in associations between structural connectivity and both higher-order cognition and general psychopathology. Together, our findings reveal a hierarchical axis in the development of structural connectivity across the human connectome.

A sensorimotor-association axis of thalamocortical connection development

Human cortical development follows a sensorimotor-to-association sequence during childhood and adolescence1-6. The brain's capacity to enact this sequence over decades indicates that it relies on intrinsic mechanisms to regulate inter-regional differences in the timing of cortical maturation, yet regulators of human developmental chronology are not well understood. Given evidence from animal models that thalamic axons modulate windows of cortical plasticity7-12, here we evaluate the overarching hypothesis that structural connections between the thalamus and cortex help to coordinate cortical maturational heterochronicity during youth. We first introduce, cortically annotate, and anatomically validate a new atlas of human thalamocortical connections using diffusion tractography. By applying this atlas to three independent youth datasets (ages 8-23 years; total N = 2,676), we reproducibly demonstrate that thalamocortical connections develop along a maturational gradient that aligns with the cortex's sensorimotor-association axis. Associative cortical regions with thalamic connections that take longest to mature exhibit protracted expression of neurochemical, structural, and functional markers indicative of higher circuit plasticity as well as heightened environmental sensitivity. This work highlights a central role for the thalamus in the orchestration of hierarchically organized and environmentally sensitive windows of cortical developmental malleability.

Connectome-based Brain Marker Moderates the Relationship between Childhood Adversity and Transdiagnostic Psychopathology during Early Adolescence

Importance

Identifying brain-based markers of resiliency that reliably predict who is and is not at elevated risk for developing psychopathology among children who experience adverse childhood experiences (ACEs) is important for improving our mechanistic understanding of these etiological links between child adversity and psychopathology and guiding precision medicine and prevention efforts for reducing psychiatric impact of ACEs.

Objective

To examine associations between ACEs and transdiagnostic psychopathology during the transition from preadolescence to early adolescence and test whether these associations are moderated by a hypothesized resilience factor, a previously identified connectome variate (CV) that is associated with higher cognitive function and lower psychopathology.

Design Setting and Participants

This study was conducted in a longitudinal design based on multicenter data from a community cohort of U.S. youth aged of 9-11 at baseline, who participated in the Adolescent Brain Cognitive Development (ABCD) study (N=7,382 at baseline and 6,813 at 2-year follow-up). Linear regression models and moderation analyses were used to characterize concurrent and prospective associations between lifetime ACEs and number of DSM-5 psychiatric disorders (indexing transdiagnostic psychopathology) and to determine if individual variations in these associations were moderated by the CV derived from resting-state fMRI at baseline.

Main Outcomes and Measures

Cumulative number of current DSM-5 psychiatric disorders assessed using the computerized self-admin version Kiddie Schedule for Affective Disorders and Schizophrenia (KSADS-5) and lifetime ACEs assessed from child and parent reports at baseline (9-10 years) and 2-year-follow-up (11-12 years).

Results

ACE total scores correlated positively with the cumulative number of current DSM-5 psychiatric disorders at both baseline (r =.258, p < .001) and 2-year follow-up (r =.257, p < .001). The baseline CV score moderated the ACE-disorder associations at baseline (B = -0.021, p < .001) and at 2-year follow-up (B = -0.018, p = .008), as well as the association between the changes in ACE and in the number of disorders from baseline to year 2 (B = -0.012, p = .045). Post-hoc analyses further showed that the moderation effect of CV on ACE-psychopathology associations was specific to the threat-related ACEs and to female youth.

Conclusions and Relevance

These findings provide preliminary evidence for a connectome-based resiliency marker and suggest that functional connectivity strength in a broad system including frontal-parietal cortices and subcortical nuclei relevant to cognitive control may protect preadolescents who have experienced lifetime ACEs--especially females and those experiencing threat-related ACEs--from developing transdiagnostic psychopathology.

Brain age has limited utility as a biomarker for capturing fluid cognition in older individuals

One well-known biomarker candidate that supposedly helps capture fluid cognition is Brain Age, or a predicted value based on machine-learning models built to predict chronological age from brain MRI. To formally evaluate the utility of Brain Age for capturing fluid cognition, we built 26 age-prediction models for Brain Age based on different combinations of MRI modalities, using the Human Connectome Project in Aging (n=504, 36-100 years old). First, based on commonality analyses, we found a large overlap between Brain Age and chronological age: Brain Age could uniquely add only around 1.6% in explaining variation in fluid cognition over and above chronological age. Second, the age-prediction models that performed better at predicting chronological age did NOT necessarily create better Brain Age for capturing fluid cognition over and above chronological age. Instead, better-performing age-prediction models created Brain Age that overlapped larger with chronological age, up to around 29% out of 32%, in explaining fluid cognition. Third, Brain Age missed around 11% of the total variation in fluid cognition that could have been explained by the brain variation. That is, directly predicting fluid cognition from brain MRI data (instead of relying on Brain Age and chronological age) could lead to around a 1/3-time improvement of the total variation explained. Accordingly, we demonstrated the limited utility of Brain Age as a biomarker for fluid cognition and made some suggestions to ensure the utility of Brain Age in explaining fluid cognition and other phenotypes of interest.

Examining the Association between Punishment and Reward Sensitivity and Response Inhibition to Previously-Incentivized Cues across Development

Processing and learning from affective cues to guide goal-directed behavior may be particularly important during adolescence; yet the factors that promote and/or disrupt the ability to integrate value in order to guide decision making across development remain unclear. The present study (N = 1046) assessed individual difference factors (self-reported punishment and reward sensitivity) related to whether previously-rewarded and previously-punished cues differentially impact goal-directed behavior (response inhibition) in a large developmental sample. Participants were between the ages of 8-21 years (Mage = 14.29, SD = 3.97, 50.38% female). Previously-rewarded cues improved response inhibition among participants age 14 and older. Further, punishment sensitivity predicted overall improved response inhibition among participants aged 10 to 18. The results highlight two main factors that are associated with improvements in the ability to integrate value to guide goal-directed behaviour - cues in the environment (e.g., reward-laden cues) and individual differences in punishment sensitivity. These findings have implications for both educational and social policies aimed at characterizing the ways in which youth integrate value to guide decision making.

Subcortico-Cortical Dysconnectivity in ADHD: A Voxel-Wise Mega-Analysis Across Multiple Cohorts

OBJECTIVE

A large body of functional MRI research has examined a potential role for subcortico-cortical loops in the pathogenesis of attention deficit hyperactivity disorder (ADHD), but has produced inconsistent findings. The authors performed a mega-analysis of six neuroimaging data sets to examine associations between ADHD diagnosis and traits and subcortico-cortical connectivity.

METHODS

Group differences were examined in the functional connectivity of four subcortical seeds in 1,696 youths with ADHD diagnoses (66.39% males; mean age, 10.83 years [SD=2.17]) and 6,737 unaffected control subjects (47.05% males; mean age, 10.33 years [SD=1.30]). The authors examined associations between functional connectivity and ADHD traits (total N=9,890; 50.3% males; mean age, 10.77 years [SD=1.96]). Sensitivity analyses were used to examine specificity relative to commonly comorbid internalizing and non-ADHD externalizing problems. The authors further examined results within motion-matched subsamples, and after adjusting for estimated intelligence.

RESULTS

In the group comparison, youths with ADHD showed greater connectivity between striatal seeds and temporal, fronto-insular, and supplementary motor regions, as well as between the amygdala and dorsal anterior cingulate cortex, compared with control subjects. Similar findings emerged when ADHD traits were considered and when alternative seed definitions were adopted. Dominant associations centered on the connectivity of the caudate bilaterally. Findings were not driven by in-scanner motion and were not shared with commonly comorbid internalizing and externalizing problems. Effect sizes were small (largest peak d, 0.15).

CONCLUSIONS

The findings from this large-scale mega-analysis support established links with subcortico-cortical circuits, which were robust to potential confounders. However, effect sizes were small, and it seems likely that resting-state subcortico-cortical connectivity can capture only a fraction of the complex pathophysiology of ADHD.

Morphometric brain organization across the human lifespan reveals increased dispersion linked to cognitive performance

The human brain is organized as segregation and integration units and follows complex developmental trajectories throughout life. The cortical manifold provides a new means of studying the brain's organization in a multidimensional connectivity gradient space. However, how the brain's morphometric organization changes across the human lifespan remains unclear. Here, leveraging structural magnetic resonance imaging scans from 1,790 healthy individuals aged 8 to 89 years, we investigated age-related global, within- and between-network dispersions to reveal the segregation and integration of brain networks from 3D manifolds based on morphometric similarity network (MSN), combining multiple features conceptualized as a "fingerprint" of an individual's brain. Developmental trajectories of global dispersion unfolded along patterns of molecular brain organization, such as acetylcholine receptor. Communities were increasingly dispersed with age, reflecting more disassortative morphometric similarity profiles within a community. Increasing within-network dispersion of primary motor and association cortices mediated the influence of age on the cognitive flexibility of executive functions. We also found that the secondary sensory cortices were decreasingly dispersed with the rest of the cortices during aging, possibly indicating a shift of secondary sensory cortices across the human lifespan from an extreme to a more central position in 3D manifolds. Together, our results reveal the age-related segregation and integration of MSN from the perspective of a multidimensional gradient space, providing new insights into lifespan changes in multiple morphometric features of the brain, as well as the influence of such changes on cognitive performance.

White matter trajectories over the lifespan

White matter (WM) changes occur throughout the lifespan at a different rate for each developmental period. We aggregated 10879 structural MRIs and 6186 diffusion-weighted MRIs from participants between 2 weeks to 100 years of age. Age-related changes in gray matter and WM partial volumes and microstructural WM properties, both brain-wide and on 29 reconstructed tracts, were investigated as a function of biological sex and hemisphere, when appropriate. We investigated the curve fit that would best explain age-related differences by fitting linear, cubic, quadratic, and exponential models to macro and microstructural WM properties. Following the first steep increase in WM volume during infancy and childhood, the rate of development slows down in adulthood and decreases with aging. Similarly, microstructural properties of WM, particularly fractional anisotropy (FA) and mean diffusivity (MD), follow independent rates of change across the lifespan. The overall increase in FA and decrease in MD are modulated by demographic factors, such as the participant's age, and show different hemispheric asymmetries in some association tracts reconstructed via probabilistic tractography. All changes in WM macro and microstructure seem to follow nonlinear trajectories, which also differ based on the considered metric. Exponential changes occurred for the WM volume and FA and MD values in the first five years of life. Collectively, these results provide novel insight into how changes in different metrics of WM occur when a lifespan approach is considered.

Lifetime residential history collection and processing for environmental data linkages in the ABCD study

By using geospatial information such as participants' residential history along with external datasets of environmental exposures, ongoing studies can enrich their cohorts to investigate the role of the environment on brain-behavior health outcomes. However, challenges may arise if clear guidance and key quality control steps are not taken at the outset of data collection of residential information. Here, we detail the protocol development aimed at improving the collection of lifetime residential address information from the Adolescent Brain Cognitive Development (ABCD) Study. This protocol generates a workflow for minimizing gaps in residential information, improving data collection processes, and reducing misclassification error in exposure estimates.

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.

Cerebral perivascular spaces as predictors of dementia risk and accelerated brain atrophy

Cerebral small vessel disease, an important risk factor for dementia, lacks robust, in vivo measurement methods. Perivascular spaces (PVS) on brain MRI are surrogates for small parenchymal blood vessels and their perivascular compartment, and may relate to brain health. We developed a novel, robust algorithm to automatically assess PVS count and size on MRI, and investigated their relationship with dementia risk and brain atrophy. We analyzed 46,478 clinical measurements of cognitive functioning and 20,845 brain MRI scans from 10,004 participants (71.1±9.7 years-old, 56.6% women). Fewer PVS and larger PVS diameter at baseline were associated with higher dementia risk and accelerated brain atrophy. Longitudinal trajectories of PVS markers were significantly different in non-demented individuals who converted to dementia compared with non-converters. In simulated placebo-controlled trials for treatments targeting cognitive decline, screening out participants less likely to develop dementia based on our PVS markers enhanced the power of the trial. These novel radiographic cerebrovascular markers may improve risk-stratification of individuals, potentially reducing cost and increasing throughput of clinical trials to combat dementia.

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.

Unveiling the Core Functional Networks of Cognition: An Ontology-Guided Machine Learning Approach

Deciphering the functional architecture that underpins diverse cognitive functions is fundamental quest in neuroscience. In this study, we employed an innovative machine learning framework that integrated cognitive ontology with functional connectivity analysis to identify brain networks essential for cognition. We identified a core assembly of functional connectomes, primarily located within the association cortex, which showed superior predictive performance compared to two conventional methods widely employed in previous research across various cognitive domains. Our approach achieved a mean prediction accuracy of 0.13 across 16 cognitive tasks, including working memory, reading comprehension, and sustained attention, outperforming the traditional methods' accuracy of 0.08. In contrast, our method showed limited predictive power for sensory, motor, and emotional functions, with a mean prediction accuracy of 0.03 across 9 relevant tasks, slightly lower than the traditional methods' accuracy of 0.04. These cognitive connectomes were further characterized by distinctive patterns of resting-state functional connectivity, structural connectivity via white matter tracts, and gene expression, highlighting their neurogenetic underpinnings. Our findings reveal a domain-general functional network fingerprint that pivotal to cognition, offering a novel computational approach to explore the neural foundations of cognitive abilities.

Negative association between neurovascular coupling and cortical gray matter volume during the lifespan

Recent studies have established the moment-to-moment turnover of the blood-oxygen-level-dependent signal (TBOLD) at resting state as a key measure of local cortical brain function. Here, we sought to extend that line of research by evaluating TBOLD in 70 cortical areas with respect to corresponding brain volume, age, and sex across the lifespan in 1,344 healthy participants including 633 from the Human Connectome Project (HCP)-Development cohort (294 males and 339 females, age range 8-21 yr) and 711 healthy participants from HCP-Aging cohort (316 males and 395 females, 36-90 yr old). In both groups, we found that 1) TBOLD increased with age, 2) volume decreased with age, and 3) TBOLD and volume were highly significantly negatively correlated, independent of age. The inverse association between TBOLD and volume was documented in nearly all 70 brain areas and for both sexes, with slightly stronger associations documented for males. The strong correspondence between TBOLD and volume across age and sex suggests a common influence such as chronic neuroinflammation contributing to reduced cortical volume and increased TBOLD across the lifespan.NEW & NOTEWORTHY We report a significant negative association between resting functional magnetic resonance imaging (fMRI) blood-oxygen-level-dependent (BOLD) signal turnover (TBOLD) and cortical gray matter volume across the lifespan, such that TBOLD increased whereas volume decreased. We attribute this association to a hypothesized chronic, low-grade neuroinflammation, probably induced by various neurotropic pathogens, including human herpes viruses known to be dormant in the brain in a latent state and reactivated by stress, fever, and various environmental exposures, such as ultraviolet light.

What is developmental about developmental prosopagnosia?

Developmental prosopagnosia (DP) is characterised by difficulties recognising face identities and is associated with diverse co-occurring object recognition difficulties. The high co-occurrence rate and heterogeneity of associated difficulties in DP is an intrinsic feature of developmental conditions, where co-occurrence of difficulties is the rule, rather than the exception. However, despite its name, cognitive and neural theories of DP rarely consider the developmental context in which these difficulties occur. This leaves a large gap in our understanding of how DP emerges in light of the developmental trajectory of face recognition. Here, we argue that progress in the field requires re-considering the developmental origins of differences in face recognition abilities, rather than studying the end-state alone. In practice, considering development in DP necessitates a re-evaluation of current approaches in recruitment, design, and analyses.

Co-representation of Functional Brain Networks Is Shaped by Cortical Myeloarchitecture and Reveals Individual Behavioral Ability

Large-scale functional networks are spatially distributed in the human brain. Despite recent progress in differentiating their functional roles, how the brain navigates the spatial coordination among them and the biological relevance of this coordination is still not fully understood. Capitalizing on canonical individualized networks derived from functional MRI data, we proposed a new concept, that is, co-representation of functional brain networks, to delineate the spatial coordination among them. To further quantify the co-representation pattern, we defined two indexes, that is, the co-representation specificity (CoRS) and intensity (CoRI), for separately measuring the extent of specific and average expression of functional networks at each brain location by using the data from both sexes. We found that the identified pattern of co-representation was anchored by cortical regions with three types of cytoarchitectural classes along a sensory-fugal axis, including, at the first end, primary (idiotypic) regions showing high CoRS, at the second end, heteromodal regions showing low CoRS and high CoRI, at the third end, paralimbic regions showing low CoRI. Importantly, we demonstrated the critical role of myeloarchitecture in sculpting the spatial distribution of co-representation by assessing the association with the myelin-related neuroanatomical and transcriptomic profiles. Furthermore, the significance of manifesting the co-representation was revealed in its prediction of individual behavioral ability. Our findings indicated that the spatial coordination among functional networks was built upon an anatomically configured blueprint to facilitate neural information processing, while advancing our understanding of the topographical organization of the brain by emphasizing the assembly of functional networks.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7 to 21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Heterogeneity Analysis on Multi-state Brain Functional Connectivity and Adolescent Neurocognition

Brain functional connectivity or connectome, a unique measure for brain functional organization, provides a great potential to explain the neurobiological underpinning of behavioral profiles. Existing connectome-based analyses highly concentrate on brain activities under a single cognitive state, and fail to consider heterogeneity when attempting to characterize brain-to-behavior relationships. In this work, we study the complex impact of multi-state functional connectivity on behaviors by analyzing the data from a recent landmark brain development and child health study. We propose a nonparametric, Bayesian supervised heterogeneity analysis to uncover neurodevelopmental subtypes with distinct effect mechanisms. We impose stochastic block structures to identify network-based functional phenotypes and develop a variational expectation-maximization algorithm to facilitate an efficient posterior computation. Through integrating resting-state and task-related functional connectomes, we dissect heterogeneous effect mechanisms on children's fluid intelligence from the functional network phenotypes including Fronto-parietal Network and Default Mode Network under different cognitive states. Based on extensive simulations, we further confirm the superior performance of our method on uncovering brain-to-behavior relationships.

Normative modeling of thalamic nuclear volumes

Thalamic nuclei have been implicated in neurodegenerative and neuropsychiatric disorders. Normative models for thalamic nuclear volumes have not been proposed thus far. The aim of this work was to establish normative models of thalamic nuclear volumes and subsequently investigate changes in thalamic nuclei in cognitive and psychiatric disorders. Volumes of the bilateral thalami and 12 nuclear regions were generated from T1 MRI data using a novel segmentation method (HIPS-THOMAS) in healthy control subjects (n=2374) and non-control subjects (n=695) with early and late mild cognitive impairment (EMCI, LMCI), Alzheimer's disease (AD), Early psychosis and Schizophrenia, Bipolar disorder, and Attention deficit hyperactivity disorder. Three different normative modelling methods were evaluated while controlling for sex, intracranial volume, and site. Z-scores and extreme z-score deviations were calculated and compared across phenotypes. GAMLSS models performed the best. Statistically significant shifts in z-score distributions consistent with atrophy were observed for most phenotypes. Shifts of progressively increasing magnitude were observed bilaterally from EMCI to AD with larger shifts in the left thalamic regions. Heterogeneous shifts were observed in psychiatric diagnoses with a predilection for the right thalamic regions. Here we present the first normative models of thalamic nuclear volumes and highlight their utility in evaluating heterogenous disorders such as Schizophrenia.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

Data leakage inflates prediction performance in connectome-based machine learning models

Predictive modeling is a central technique in neuroimaging to identify brain-behavior relationships and test their generalizability to unseen data. However, data leakage undermines the validity of predictive models by breaching the separation between training and test data. Leakage is always an incorrect practice but still pervasive in machine learning. Understanding its effects on neuroimaging predictive models can inform how leakage affects existing literature. Here, we investigate the effects of five forms of leakage-involving feature selection, covariate correction, and dependence between subjects-on functional and structural connectome-based machine learning models across four datasets and three phenotypes. Leakage via feature selection and repeated subjects drastically inflates prediction performance, whereas other forms of leakage have minor effects. Furthermore, small datasets exacerbate the effects of leakage. Overall, our results illustrate the variable effects of leakage and underscore the importance of avoiding data leakage to improve the validity and reproducibility of predictive modeling.