Neuroimaging of the Philadelphia neurodevelopmental cohort

Longitudinal Development of Neurocognitive Functioning and Gray Matter Volume in Youths With Recurrent Psychosis Spectrum Symptoms

BACKGROUND AND HYPOTHESIS

Neurodevelopmental risk-factor models of psychosis highlight the importance of early developmental deviations in the emergence of psychosis. However, few longitudinal studies map neurodevelopment and neurocognitive trajectories across age in preclinical psychosis. We investigated longitudinal trajectories in neurocognition and brain volume in a community cohort of adolescents with recurrent psychosis spectrum (PS) symptoms, tracking their development into young adulthood compared to their typically developing (TD) peers.

STUDY DESIGN

Utilizing the Philadelphia Neurodevelopmental Cohort, we analyzed data of 231 youths aged 8-30 with at least one follow-up assessment, including 88 with PS.

STUDY RESULTS

Individuals with PS showed similar developmental trajectories but demonstrated significant impairments in executive functioning (t = -2.81, q = 0.010), memory (t = -2.34, q = 0.019), complex cognition (t = -3.72, q = 0.001), social cognition (t = -2.73, q = 0.010), motor (t = -2.50, q = 0.015), and general cognition (t = -3.20, q = 0.004). Lower cortical (t = -2.46, P = .014) and subcortical (t = -2.41, P = .016) gray matter volume in the recurrent PS group compared to the TD group were documented with age-related group differences becoming less pronounced by young adulthood. Further analyses revealed age-by-group interactions (qs < 0.05) observed in a few temporal and frontal regions, with differences between groups at earlier ages.

CONCLUSIONS

These findings suggest that recurrent PS symptoms are linked to early neurocognitive and brain structure deficits, highlighting the need for interventions to reduce psychosis risk and support healthy neurodevelopment.

Diverging Effects of Violence Exposure and Psychiatric Symptoms on Amygdala-Prefrontal Maturation During Childhood and Adolescence

BACKGROUND

Violence exposure during childhood and adolescence is associated with increased prevalence and severity of psychopathology. Neurobiological correlates suggest that abnormal maturation of emotion-related brain circuitry, such as the amygdala-prefrontal cortex (PFC) circuit, may underlie the development of psychiatric symptoms after exposure. However, it remains unclear how amygdala-PFC circuit maturation is related to psychiatric risk in the context of violence.

METHODS

In this study, we analyzed individual differences in amygdala-PFC circuit maturity using data collected from the PNC (Philadelphia Neurodevelopmental Cohort) (n = 1133 youths). Neurodevelopment models of amygdala-PFC resting-state functional connectivity were built using deep learning and trained to predict chronological age in typically developing youths (not violence exposed and without a psychiatric diagnosis). Using the brain age gap estimate, an index of relative circuit maturation, patterns of atypical neurodevelopment were investigated.

RESULTS

Violence exposure was associated with delayed maturation of basolateral amygdala (BLA)-PFC circuits, driven by increased BLA-medial orbitofrontal cortex functional connectivity. In contrast, increased psychiatric symptoms were associated with advanced maturation of BLA-PFC functional connectivity, driven by decreased BLA-dorsolateral PFC functional connectivity.

CONCLUSIONS

Delayed frontoamygdala maturation after exposure to violence suggests atypical, but adaptive, development of threat appraisal processes, potentially reflecting a greater threat generalization characteristic of younger children. Advanced circuit maturation with increasing symptoms suggests divergent neurodevelopmental mechanisms underlying illness after emotion circuits have adapted to adversity, exacerbated by preexisting vulnerabilities to early maturation. Disentangling the effects of adversity and psychopathology on neurodevelopment is crucial for helping youths recover from violence and preventing illness from continuing into adulthood.

Presence, severity, and functional associations of incomplete hippocampal inversion in 22q11.2 deletion syndrome

BACKGROUND

The hippocampus is smaller and functionally disrupted in individuals with 22q11.2 deletion syndrome (22q11DS), though the cause remains unclear. During gestational weeks 20-30 an inversion in the dentate gyrus and cornu ammonis occurs. This process can go awry resulting in incomplete hippocampal inversion (IHI). In the general population, IHI is more common in the left hemisphere than the right; yet its prevalence, severity, and functional impact in 22q11DS remain unexplored. Investigating IHI in 22q11DS could uncover morphological hippocampal abnormalities linked to neuropsychiatric and neurocognitive symptoms.

METHODS

Using 3T structural MRI data, the presence and severity of IHI were assessed in 22q11DS (n=108) and healthy comparison (HC; n=633) individuals. Total and subregional hippocampal volume, psychopathology, and hippocampal-based memory were evaluated.

RESULTS

IHI prevalence was significantly higher in 22q11DS compared to HC in both the left (63% vs. 30%, p<0.001) and right hemispheres (29% vs. 8%, p<0.001). IHI severity was also greater in 22q11DS (p<0.001) bilaterally. IHI influenced hippocampal volume differences, with left IHI primarily affecting the head (p<0.01) and tail (p<0.001) and right IHI affecting only the tail (p<0.001). In exploratory analyses within 22q11DS, left IHI presence was linked to poorer face memory (p<0.05) but not psychopathology.

CONCLUSIONS

These findings highlight a high prevalence of hippocampal morphological alterations in 22q11DS, which are associated with memory performance. Earlier developmental and longitudinal studies are needed to clarify the role of IHI in 22q11DS sequelae.

Optimizing Biophysical Large-Scale Brain Circuit Models With Deep Neural Networks

Biophysical modeling provides mechanistic insights into brain function, from single-neuron dynamics to large-scale circuit models bridging macro-scale brain activity with microscale measurements. Biophysical models are governed by biologically meaningful parameters, many of which can be experimentally measured. Some parameters are unknown, and optimizing their values can dramatically improve adherence to experimental data, significantly enhancing biological plausibility. Previous optimization methods - such as exhaustive search, gradient descent, evolutionary strategies and Bayesian optimization - require repeated, computationally expensive numerical integration of biophysical differential equations, limiting scalability to population-level datasets. Here, we introduce DELSSOME (DEep Learning for Surrogate Statistics Optimization in MEan field modeling), a framework that bypasses numerical integration by directly predicting whether model parameters produce realistic brain dynamics. When applied to the widely used feedback inhibition control (FIC) mean field model, DELSSOME achieves a 2000× speedup over Euler integration. By embedding DELSSOME within an evolutionary optimization strategy, trained models generalize to new datasets without additional tuning, enabling a 50× speedup in FIC model estimation while preserving neurobiological insights. The massive acceleration facilitates large-scale mechanistic modeling in population-level neuroscience, unlocking new opportunities for understanding brain function.

Dev-Atlas: A reference atlas of functional brain networks for typically developing adolescents

It is well accepted that the brain is functionally organized into multiple networks and extensive literature has demonstrated that the organization of these networks shows major changes during adolescence. Yet, there is limited option for a reference functional brain atlas derived from typically-developing adolescents, which is problematic as the reliable identification of functional brain networks crucially depends on the use of such reference functional atlases. In this context, we utilized resting-state functional MRI data from 1391 typically-developing youth aged 8-17 years to create an adolescent-specific reference atlas of functional brain networks. We further investigated the impact of age and sex on these networks. Using a multiscale individual component clustering algorithm, we identified 24 reliable functional brain networks, classified within six domains: Default-Mode (5 networks), Control (4 networks), Salience (3 networks), Attention (4 networks), Somatomotor (5 networks), and Visual (3 networks). We identified reliable and large effects of age on the spatial topography of these majority of networks, as well as on the functional network connectivity. Sex effects were not as widespread. We created a novel brain atlas, named Dev-Atlas, focused on a typically-developing sample, with the hope that this atlas can be used in future developmental neuroscience studies.

A lightweight generative model for interpretable subject-level prediction

Recent years have seen a growing interest in methods for predicting an unknown variable of interest, such as a subject's diagnosis, from medical images depicting its anatomical-functional effects. Methods based on discriminative modeling excel at making accurate predictions, but are challenged in their ability to explain their decisions in anatomically meaningful terms. In this paper, we propose a simple technique for single-subject prediction that is inherently interpretable. It augments the generative models used in classical human brain mapping techniques, in which the underlying cause-effect relations can be encoded, with a multivariate noise model that captures dominant spatial correlations. Experiments demonstrate that the resulting model can be efficiently inverted to make accurate subject-level predictions, while at the same time offering intuitive visual explanations of its inner workings. The method is easy to use: training is fast for typical training set sizes, and only a single hyperparameter needs to be set by the user. Our code is available at https://github.com/chiara-mauri/Interpretable-subject-level-prediction.

Two Axes of White Matter Development

Despite decades of neuroimaging research, how white matter develops along the length of major tracts in humans remains unknown. Here, we identify fundamental patterns of white matter maturation by examining developmental variation along major, long-range cortico-cortical tracts in youth ages 5-23 years using diffusion MRI from three large-scale, cross-sectional datasets (total N = 2,710). Across datasets, we delineate two replicable axes of human white matter development. First, we find a deep-to-superficial axis, in which superficial tract regions near the cortical surface exhibit greater age-related change than deep tract regions. Second, we demonstrate that the development of superficial tract regions aligns with the cortical hierarchy defined by the sensorimotor-association axis, with tract ends adjacent to sensorimotor cortices maturing earlier than those adjacent to association cortices. These results reveal developmental variation along tracts that conventional tract-average analyses have previously obscured, challenging the implicit assumption that white matter tracts mature uniformly along their length. Such developmental variation along tracts may have functional implications, including mitigating ephaptic coupling in densely packed deep tract regions and tuning neural synchrony through hierarchical development in superficial tract regions - ultimately refining neural transmission in youth.

Tensor dictionary-based heterogeneous transfer learning to study emotion-related gender differences in brain

In practice, collecting auxiliary labeled data with same feature space from multiple domains is difficult. Thus, we focus on the heterogeneous transfer learning to address the problem of insufficient sample sizes in neuroimaging. Viewing subjects, time, and features as dimensions, brain activation and dynamic functional connectivity data can be treated as high-order heterogeneous data with heterogeneity arising from distinct feature space. To use the heterogeneous priori knowledge from the low-dimensional brain activation data to improve the classification performance of high-dimensional dynamic functional connectivity data, we propose a tensor dictionary-based heterogeneous transfer learning framework. It combines supervised tensor dictionary learning with heterogeneous transfer learning for enhance high-order heterogeneous knowledge sharing. The former can encode the underlying discriminative features in high-order data into dictionaries, while the latter can transfer heterogeneous knowledge encoded in dictionaries through feature transformation derived from mathematical relationship between domains. The primary focus of this paper is gender classification using fMRI data to identify emotion-related brain gender differences during adolescence. Additionally, experiments on simulated data and EEG data are included to demonstrate the generalizability of the proposed method. Experimental results indicate that incorporating prior knowledge significantly enhances classification performance. Further analysis of brain gender differences suggests that temporal variability in brain activity explains differences in emotion regulation strategies between genders. By adopting the heterogeneous knowledge sharing strategy, the proposed framework can capture the multifaceted characteristics of the brain, improve the generalization of the model, and reduce training costs. Understanding the gender specific neural mechanisms of emotional cognition helps to develop the gender-specific treatments for neurological diseases.

Reproducible Brain Charts: An open data resource for mapping brain development and its associations with mental health

Major mental disorders are increasingly understood as disorders of brain development. Large and heterogeneous samples are required to define generalizable links between brain development and psychopathology. To this end, we introduce the Reproducible Brain Charts (RBC), an open data resource that integrates data from 5 large studies of brain development in youth from three continents (N=6,346; 45% Female). Confirmatory bifactor models were used to create harmonized psychiatric phenotypes that capture major dimensions of psychopathology. Following rigorous quality assurance, neuroimaging data were carefully curated and processed using consistent pipelines in a reproducible manner with DataLad, the Configurable Pipeline for the Analysis of Connectomes (C-PAC), and FreeSurfer. Initial analyses of RBC data emphasize the benefit of careful quality assurance and data harmonization in delineating developmental effects and associations with psychopathology. Critically, all RBC data - including harmonized psychiatric phenotypes, unprocessed images, and fully processed imaging derivatives - are openly shared without a data use agreement via the International Neuroimaging Data-sharing Initiative. Together, RBC facilitates large-scale, reproducible, and generalizable research in developmental and psychiatric neuroscience.

Regression and alignment for functional data and network topology

In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of preprocessing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.

Approximating R1 and R2: A Quantitative Approach to Clinical Weighted MRI

Weighted MRI images are widely used in clinical as well as open-source neuroimaging databases. Weighted images such as T1-weighted, T2-weighted, and proton density-weighted (T1w, T2w, and PDw, respectively) are used for evaluating the brain's macrostructure; however, their values cannot be used for microstructural analysis, as they lack physical meaning. Quantitative MRI (qMRI) relaxation rate parameters (e.g., R1 and R2) do contain microstructural physical meaning. Nevertheless, qMRI is rarely done in large-scale clinical databases. Currently, the weighted images ratio T1w/T2w is used as a quantifier to approximate the brain's microstructure. In this paper, we test three additional quantifiers that approximate quantitative maps, which can help bring quantitative MRI to the clinic for easy use. Following the signal equations and using simple mathematical operations, we combine the T1w, T2w, and PDw images to estimate the R1 and R2 maps. We find that two of these quantifiers (T1w/PDw and T1w/ln(T2w)) can approximate R1, and that (ln(T2w/PDw)) can approximate R2, in 3 datasets that were tested. We find that this approach also can be applied to T2w scans taken from widely available DTI (Diffusion Tensor Imaging) datasets. We tested these quantifiers on both in vitro phantom and in vivo human datasets. We found that the quantifiers accurately represent the quantitative parameters across datasets. Finally, we tested the quantifiers within a clinical context, and found that they are robust across datasets. Our work provides a simple pipeline to enhance the usability and quantitative accuracy of MRI weighted images.

Distinguishing deception from its confounds by improving the validity of fMRI-based neural prediction

Deception is a universal human behavior. Yet longstanding skepticism about the validity of measures used to characterize the biological mechanisms underlying deceptive behavior has relegated such studies to the scientific periphery. Here, we address these fundamental questions by applying machine learning methods and functional magnetic resonance imaging (fMRI) to signaling games capturing motivated deception in human participants. First, we develop an approach to test for the presence of confounding processes and validate past skepticism by showing that much of the predictive power of neural predictors trained on deception data comes from processes other than deception. Specifically, we demonstrate that discriminant validity is compromised by the predictor's ability to predict behavior in a control task that does not involve deception. Second, we show that the presence of confounding signals need not be fatal and that the validity of the neural predictor can be improved by removing confounding signals while retaining those associated with the task of interest. To this end, we develop a "dual-goal tuning" approach in which, beyond the typical goal of predicting the behavior of interest, the predictor also incorporates a second compulsory goal that enforces chance performance in the control task. Together, these findings provide a firmer scientific foundation for understanding the neural basis of a neglected class of behavior, and they suggest an approach for improving validity of neural predictors.

Tbx1 haploinsufficiency leads to local skull deformity, paraflocculus and flocculus dysplasia, and motor-learning deficit in 22q11.2 deletion syndrome

Neurodevelopmental disorders are thought to arise from intrinsic brain abnormalities. Alternatively, they may arise from disrupted crosstalk among tissues. Here we show the local reduction of two vestibulo-cerebellar lobules, the paraflocculus and flocculus, in mouse models and humans with 22q11.2 deletion syndrome (22q11DS). In mice, this paraflocculus/flocculus dysplasia is associated with haploinsufficiency of the Tbx1 gene. Tbx1 haploinsufficiency also leads to impaired cerebellar synaptic plasticity and motor learning. However, neural cell compositions and neurogenesis are not altered in the dysplastic paraflocculus/flocculus. Interestingly, 22q11DS and Tbx1+/- mice have malformations of the subarcuate fossa, a part of the petrous temporal bone, which encapsulates the paraflocculus/flocculus. Single-nuclei RNA sequencing reveals that Tbx1 haploinsufficiency leads to precocious differentiation of chondrocytes to osteoblasts in the petrous temporal bone autonomous to paraflocculus/flocculus cell populations. These findings suggest a previously unrecognized pathogenic structure/function relation in 22q11DS in which local skeletal deformity and cerebellar dysplasia result in behavioral deficiencies.

A Deep Dynamic Causal Learning Model to Study Changes in Dynamic Effective Connectivity During Brain Development

OBJECTIVE

Brain dynamic effective connectivity (dEC), characterizes the information transmission patterns between brain regions that change over time, which provides insight into the biological mechanism underlying brain development. However, most existing methods predominantly capture fixed or temporally invariant EC, leaving dEC largely unexplored.

METHODS

Herein we propose a deep dynamic causal learning model specifically designed to capture dEC. It includes a dynamic causal learner to detect time-varying causal relationships from spatio-temporal data, and a dynamic causal discriminator to validate these findings by comparing original and reconstructed data.

RESULTS

Our model outperforms established baselines in the accuracy of identifying dynamic causalities when tested on the simulated data. When applied to the Philadelphia Neurodevelopmental Cohort, the model uncovers distinct patterns in dEC networks across different age groups. Specifically, the evolution process of brain dEC networks in young adults is more stable than in children, and significant differences in information transfer patterns exist between them.

CONCLUSION

This study highlights the brain's developmental trajectory, where networks transition from undifferentiated to specialized structures with age, in accordance with the improvement of an individual's cognitive and information processing capability.

SIGNIFICANCE

The proposed model consists of the identification and verification of dynamic causality, utilizing the spatio-temporal fusing information from fMRI. As a result, it can accurately detect dEC and characterize its evolution over age.

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.

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Normative Neuroimaging Library: Designing a Comprehensive and Demographically Diverse Dataset of Healthy Controls to Support Traumatic Brain Injury Diagnostic and Therapeutic Development

The past decade has seen impressive advances in neuroimaging, moving from qualitative to quantitative outputs. Available techniques now allow for the inference of microscopic changes occurring in white and gray matter, along with alterations in physiology and function. These existing and emerging techniques hold the potential of providing unprecedented capabilities in achieving a diagnosis and predicting outcomes for traumatic brain injury (TBI) and a variety of other neurological diseases. To see this promise move from the research lab into clinical care, an understanding is needed of what normal data look like for all age ranges, sex, and other demographic and socioeconomic categories. Clinicians can only use the results of imaging scans to support their decision-making if they know how the results for their patient compare with a normative standard. This potential for utilizing magnetic resonance imaging (MRI) in TBI diagnosis motivated the American College of Radiology and Cohen Veterans Bioscience to create a reference database of healthy individuals with neuroimaging, demographic data, and characterization of psychological functioning and neurocognitive data that will serve as a normative resource for clinicians and researchers for development of diagnostics and therapeutics for TBI and other brain disorders. The goal of this article is to introduce the large, well-curated Normative Neuroimaging Library (NNL) to the research community. NNL consists of data collected from ∼1900 healthy participants. The highlights of NNL are (1) data are collected across a diverse population, including civilians, veterans, and active-duty service members with an age range (18-64 years) not well represented in existing datasets; (2) comprehensive structural and functional neuroimaging acquisition with state-of-the-art sequences (including structural, diffusion, and functional MRI; raw scanner data are preserved, allowing higher quality data to be derived in the future; standardized imaging acquisition protocols across sites reflect sequences and parameters often recommended for use with various neurological and psychiatric conditions, including TBI, post-traumatic stress disorder, stroke, neurodegenerative disorders, and neoplastic disease); and (3) the collection of comprehensive demographic details, medical history, and a broad structured clinical assessment, including cognition and psychological scales, relevant to multiple neurological conditions with functional sequelae. Thus, NNL provides a demographically diverse population of healthy individuals who can serve as a comparison group for brain injury study and clinical samples, providing a strong foundation for precision medicine. Use cases include the creation of imaging-derived phenotypes (IDPs), derivation of reference ranges of imaging measures, and use of IDPs as training samples for artificial intelligence-based biomarker development and for normative modeling to help identify injury-induced changes as outliers for precision diagnosis and targeted therapeutic development. On its release, NNL is poised to support the use of advanced imaging in clinician decision support tools, the validation of imaging biomarkers, and the investigation of brain-behavior anomalies, moving the field toward precision medicine.

Psychiatric neuroimaging at a crossroads: Insights from psychiatric genetics

Thanks to methodological advances, large-scale data collections, and longitudinal designs, psychiatric neuroimaging is better equipped than ever to identify the neurobiological underpinnings of youth mental health problems. However, the complexity of such endeavors has become increasingly evident, as the field has been confronted by limited clinical relevance, inconsistent results, and small effect sizes. Some of these challenges parallel those historically encountered by psychiatric genetics. In past genetic research, robust findings were historically undermined by oversimplified biological hypotheses, mistaken assumptions about expectable effect sizes, replication problems, confounding by population structure, and shared biological patterns across disorders. Overcoming these challenges has contributed to current successes in the field. Drawing parallels across psychiatric genetics and neuroimaging, we identify key shared challenges as well as pinpoint relevant insights that could be gained in psychiatric neuroimaging from the transition that occurred from the candidate gene to (post) genome-wide "eras" of psychiatric genetics. Finally, we discuss the prominent developmental component of psychiatric neuroimaging and how that might be informed by epidemiological and omics approaches. The evolution of psychiatric genetic research offers valuable insights that may expedite the resolution of key challenges in psychiatric neuroimaging, thus potentially moving our understanding of psychiatric pathophysiology forward.

The General Psychopathology 'p' Factor in Adolescence: Multi-Informant Assessment and Computerized Adaptive Testing

Accumulating evidence supports the presence of a general psychopathology dimension, the p factor ('p'). Despite growing interest in the p factor, questions remain about how p is assessed. Although multi-informant assessment of psychopathology is commonplace in clinical research and practice with children and adolescents, almost no research has taken a multi-informant approach to studying youth p or has examined the degree of concordance between parent and youth reports. Further, estimating p requires assessment of a large number of symptoms, resulting in high reporter burden that may not be feasible in many clinical and research settings. In the present study, we used bifactor multidimensional item response theory models to estimate parent- and adolescent-reported p in a large community sample of youth (11-17 years) and parents (N = 5,060 dyads). We examined agreement between parent and youth p scores and associations with assessor-rated youth global functioning. We also applied computerized adaptive testing (CAT) simulations to parent and youth reports to determine whether adaptive testing substantially alters agreement on p or associations with youth global functioning. Parent-youth agreement on p was moderate (r =.44) and both reports were negatively associated with youth global functioning. Notably, 7 out of 10 of the highest loading items were common across reporters. CAT reduced the average number of items administered by 57%. Agreement between CAT-derived p scores was similar to the full form (r =.40) and CAT scores were negatively correlated with youth functioning. These novel results highlight the promise and potential clinical utility of a multi-informant p factor approach.

Similarity-based multimodal regression

To better understand complex human phenotypes, large-scale studies have increasingly collected multiple data modalities across domains such as imaging, mobile health, and physical activity. The properties of each data type often differ substantially and require either separate analyses or extensive processing to obtain comparable features for a combined analysis. Multimodal data fusion enables certain analyses on matrix-valued and vector-valued data, but it generally cannot integrate modalities of different dimensions and data structures. For a single data modality, multivariate distance matrix regression provides a distance-based framework for regression accommodating a wide range of data types. However, no distance-based method exists to handle multiple complementary types of data. We propose a novel distance-based regression model, which we refer to as Similarity-based Multimodal Regression (SiMMR), that enables simultaneous regression of multiple modalities through their distance profiles. We demonstrate through simulation, imaging studies, and longitudinal mobile health analyses that our proposed method can detect associations between clinical variables and multimodal data of differing properties and dimensionalities, even with modest sample sizes. We perform experiments to evaluate several different test statistics and provide recommendations for applying our method across a broad range of scenarios.

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.

A Bayesian incorporated linear non-Gaussian acyclic model for multiple directed graph estimation to study brain emotion circuit development in adolescence

Emotion perception is essential to affective and cognitive development which involves distributed brain circuits. Emotion identification skills emerge in infancy and continue to develop throughout childhood and adolescence. Understanding the development of the brain's emotion circuitry may help us explain the emotional changes during adolescence. In this work, we aim to deepen our understanding of emotion-related functional connectivity (FC) from association to causation. We proposed a Bayesian incorporated linear non-Gaussian acyclic model (BiLiNGAM), which incorporated association model into the estimation pipeline. Simulation results indicated stable and accurate performance over various settings, especially when the sample size was small. We used fMRI data from the Philadelphia Neurodevelopmental Cohort (PNC) to validate the approach. It included 855 individuals aged 8-22 years who were divided into five different adolescent stages. Our network analysis revealed the development of emotion-related intra- and intermodular connectivity and pinpointed several emotion-related hubs. We further categorized the hubs into two types: in-hubs and out-hubs, as the center of receiving and distributing information, respectively. In addition, several unique developmental hub structures and group-specific patterns were discovered. Our findings help provide a directed FC template of brain network organization underlying emotion processing during adolescence.

Explainable spatio-temporal graph evolution learning with applications to dynamic brain network analysis during development

Modeling dynamic interactions among network components is crucial to uncovering the evolution mechanisms of complex networks. Recently, spatio-temporal graph learning methods have achieved noteworthy results in characterizing the dynamic changes of inter-node relations (INRs). However, challenges remain: The spatial neighborhood of an INR is underexploited, and the spatio-temporal dependencies in INRs' dynamic changes are overlooked, ignoring the influence of historical states and local information. In addition, the model's explainability has been understudied. To address these issues, we propose an explainable spatio-temporal graph evolution learning (ESTGEL) model to model the dynamic evolution of INRs. Specifically, an edge attention module is proposed to utilize the spatial neighborhood of an INR at multi-level, i.e., a hierarchy of nested subgraphs derived from decomposing the initial node-relation graph. Subsequently, a dynamic relation learning module is proposed to capture the spatio-temporal dependencies of INRs. The INRs are then used as adjacent information to improve the node representation, resulting in comprehensive delineation of dynamic evolution of the network. Finally, the approach is validated with real data on brain development study. Experimental results on dynamic brain networks analysis reveal that brain functional networks transition from dispersed to more convergent and modular structures throughout development. Significant changes are observed in the dynamic functional connectivity (dFC) associated with functions including emotional control, decision-making, and language processing.

The role of functional emotion circuits in distinct dimensions of psychopathology in youth

Several mental disorders emerge during childhood or adolescence and are often characterized by socioemotional difficulties, including alterations in emotion perception. Emotional facial expressions are processed in discrete functional brain modules whose connectivity patterns encode emotion categories, but the involvement of these neural circuits in psychopathology in youth is poorly understood. This study examined the associations between activation and functional connectivity patterns in emotion circuits and psychopathology during development. We used task-based fMRI data from the Philadelphia Neurodevelopmental Cohort (PNC, N = 1221, 8-23 years) and conducted generalized psycho-physiological interaction (gPPI) analyses. Measures of psychopathology were derived from an independent component analysis of questionnaire data. The results showed positive associations between identifying fearful, sad, and angry faces and depressive symptoms, and a negative relationship between sadness recognition and positive psychosis symptoms. We found a positive main effect of depressive symptoms on BOLD activation in regions overlapping with the default mode network, while individuals reporting higher levels of norm-violating behavior exhibited emotion-specific lower functional connectivity within regions of the salience network and between modules that overlapped with the salience and default mode network. Our findings illustrate the relevance of functional connectivity patterns underlying emotion processing for behavioral problems in children and adolescents.

Large-scale analysis of structural brain asymmetries during neurodevelopment: Associations with age and sex in 4265 children and adolescents

Only a small number of studies have assessed structural differences between the two hemispheres during childhood and adolescence. However, the existing findings lack consistency or are restricted to a particular brain region, a specific brain feature, or a relatively narrow age range. Here, we investigated associations between brain asymmetry and age as well as sex in one of the largest pediatric samples to date (n = 4265), aged 1-18 years, scanned at 69 sites participating in the ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) consortium. Our study revealed that significant brain asymmetries already exist in childhood, but their magnitude and direction depend on the brain region examined and the morphometric measurement used (cortical volume or thickness, regional surface area, or subcortical volume). With respect to effects of age, some asymmetries became weaker over time while others became stronger; sometimes they even reversed direction. With respect to sex differences, the total number of regions exhibiting significant asymmetries was larger in females than in males, while the total number of measurements indicating significant asymmetries was larger in males (as we obtained more than one measurement per cortical region). The magnitude of the significant asymmetries was also greater in males. However, effect sizes for both age effects and sex differences were small. Taken together, these findings suggest that cerebral asymmetries are an inherent organizational pattern of the brain that manifests early in life. Overall, brain asymmetry appears to be relatively stable throughout childhood and adolescence, with some differential effects in males and females.

Domain-specific associations between psychopathology and neurocognitive functioning

BACKGROUND

Neurocognitive dysfunction is a transdiagnostic finding in psychopathology, but relationships among cognitive domains and general and specific psychopathology dimensions remain unclear. This study aimed to examine associations between cognition and psychopathology dimensions in a large youth cohort.

METHOD

The sample (N = 9350; age 8-21 years) was drawn from the Philadelphia Neurodevelopmental Cohort. Data from structured clinical interviews were modeled using bifactor confirmatory factor analysis (CFA), resulting in an overall psychopathology ('p') factor score and six orthogonal psychopathology dimensions: dysphoria/distress, obsessive-compulsive, behavioral/externalizing, attention-deficit/hyperactivity, phobias, and psychosis. Neurocognitive data were aggregated using correlated-traits CFA into five factors: executive functioning, memory, complex cognition, social cognition, and sensorimotor speed. We examined relationships among specific and general psychopathology dimensions and neurocognitive factors.

RESULTS

The final model showed both overall and specific associations between cognitive functioning and psychopathology, with acceptable fit (CFI = 0.91; TLI = 0.90; RMSEA = 0.024; SRMR = 0.054). Overall psychopathology and most psychopathology dimensions were negatively associated with neurocognitive functioning (phobias [p < 0.0005], behavioral/externalizing [p < 0.0005], attention-deficit/hyperactivity [p < 0.0005], psychosis [p < 0.0005 to p < 0.05]), except for dysphoria/distress and obsessive-compulsive symptoms, which were positively associated with complex cognition (p < 0.05 and p < 0.01, respectively).

CONCLUSION

By modeling a broad range of cognitive and psychopathology domains in a large, diverse sample of youth, we found aspects of neurocognitive functioning shared across clinical phenotypes, as well as domain-specific patterns. Findings support transdiagnostic examination of cognitive performance to parse variability in the link between neurocognitive functioning and clinical phenotypes.

Adolescent maturation of cortical excitation-inhibition balance based on individualized biophysical network modeling

The balance of excitation and inhibition is a key functional property of cortical microcircuits which changes through the lifespan. Adolescence is considered a crucial period for the maturation of excitation-inhibition balance. This has been primarily observed in animal studies, yet human in vivo evidence on adolescent maturation of the excitation-inhibition balance at the individual level is limited. Here, we developed an individualized in vivo marker of regional excitation-inhibition balance in human adolescents, estimated using large-scale simulations of biophysical network models fitted to resting-state functional magnetic resonance imaging data from two independent cross-sectional (N = 752) and longitudinal (N = 149) cohorts. We found a widespread relative increase of inhibition in association cortices paralleled by a relative age-related increase of excitation, or lack of change, in sensorimotor areas across both datasets. This developmental pattern co-aligned with multiscale markers of sensorimotor-association differentiation. The spatial pattern of excitation-inhibition development in adolescence was robust to inter-individual variability of structural connectomes and modeling configurations. Notably, we found that alternative simulation-based markers of excitation-inhibition balance show a variable sensitivity to maturational change. Taken together, our study highlights an increase of inhibition during adolescence in association areas using cross sectional and longitudinal data, and provides a robust computational framework to estimate microcircuit maturation in vivo at the individual level.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Network enrichment significance testing in brain-phenotype association studies

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

A Demographic-Conditioned Variational Autoencoder for fMRI Distribution Sampling and Removal of Confounds

Objective

fMRI and derived measures such as functional connectivity (FC) have been used to predict brain age, general fluid intelligence, psychiatric disease status, and preclinical neurodegenerative disease. However, it is not always clear that all demographic confounds, such as age, sex, and race, have been removed from fMRI data. Additionally, many fMRI datasets are restricted to authorized researchers, making dissemination of these valuable data sources challenging.

Methods

We create a variational autoencoder (VAE)-based model, DemoVAE, to decorrelate fMRI features from demographics and generate high-quality synthetic fMRI data based on user-supplied demographics. We train and validate our model using two large, widely used datasets, the Philadelphia Neurodevel-opmental Cohort (PNC) and Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP).

Results

We find that DemoVAE recapitulates group differences in fMRI data while capturing the full breadth of individual variations. Significantly, we also find that most clinical and computerized battery fields that are correlated with fMRI data are not correlated with DemoVAE latents. An exception are several fields related to schizophrenia medication and symptom severity.

Conclusion

Our model generates fMRI data that captures the full distribution of FC better than traditional VAE or GAN models. We also find that most prediction using fMRI data is dependent on correlation with, and prediction of, demographics.

Significance

Our DemoVAE model allows for generation of high quality synthetic data conditioned on subject demographics as well as the removal of the confounding effects of demographics. We identify that FC-based prediction tasks are highly influenced by demographic confounds.

A Demographic-Conditioned Variational Autoencoder for fMRI Distribution Sampling and Removal of Confounds

Objective

fMRI and derived measures such as functional connectivity (FC) have been used to predict brain age, general fluid intelligence, psychiatric disease status, and preclinical neurodegenerative disease. However, it is not always clear that all demographic confounds, such as age, sex, and race, have been removed from fMRI data. Additionally, many fMRI datasets are restricted to authorized researchers, making dissemination of these valuable data sources challenging.

Methods

We create a variational autoencoder (VAE)-based model, DemoVAE, to decorrelate fMRI features from demographics and generate high-quality synthetic fMRI data based on user-supplied demographics. We train and validate our model using two large, widely used datasets, the Philadelphia Neurodevelopmental Cohort (PNC) and Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP).

Results

We find that DemoVAE recapitulates group differences in fMRI data while capturing the full breadth of individual variations. Significantly, we also find that most clinical and computerized battery fields that are correlated with fMRI data are not correlated with DemoVAE latents. An exception are several fields related to schizophrenia medication and symptom severity.

Conclusion

Our model generates fMRI data that captures the full distribution of FC better than traditional VAE or GAN models. We also find that most prediction using fMRI data is dependent on correlation with, and prediction of, demographics.

Significance

Our DemoVAE model allows for generation of high quality synthetic data conditioned on subject demographics as well as the removal of the confounding effects of demographics. We identify that FC-based prediction tasks are highly influenced by demographic confounds.

pNet: A toolbox for personalized functional networks modeling

Personalized functional networks (FNs) derived from functional magnetic resonance imaging (fMRI) data are useful for characterizing individual variations in the brain functional topography associated with the brain development, aging, and disorders. To facilitate applications of the personalized FNs with enhanced reliability and reproducibility, we develop an open-source toolbox that is user-friendly, extendable, and includes rigorous quality control (QC), featuring multiple user interfaces (graphics, command line, and a step-by-step guideline) and job-scheduling for high performance computing (HPC) clusters. Particularly, the toolbox, named personalized functional network modeling (pNet), takes fMRI inputs in either volumetric or surface type, ensuring compatibility with multiple fMRI data formats, and computes personalized FNs using two distinct modeling methods: one method optimizes the functional coherence of FNs, while the other enhances their independence. Additionally, the toolbox provides HTML-based reports for QC and visualization of personalized FNs. The toolbox is developed in both MATLAB and Python platforms with a modular design to facilitate extension and modification by users familiar with either programming language. We have evaluated the toolbox on two fMRI datasets and demonstrated its effectiveness and user-friendliness with interactive and scripting examples. pNet is publicly available at https://github.com/MLDataAnalytics/pNet.

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.

Interpretable Cognitive Ability Prediction: A Comprehensive Gated Graph Transformer Framework for Analyzing Functional Brain Networks

Graph convolutional deep learning has emerged as a promising method to explore the functional organization of the human brain in neuroscience research. This paper presents a novel framework that utilizes the gated graph transformer (GGT) model to predict individuals' cognitive ability based on functional connectivity (FC) derived from fMRI. Our framework incorporates prior spatial knowledge and uses a random-walk diffusion strategy that captures the intricate structural and functional relationships between different brain regions. Specifically, our approach employs learnable structural and positional encodings (LSPE) in conjunction with a gating mechanism to efficiently disentangle the learning of positional encoding (PE) and graph embeddings. Additionally, we utilize the attention mechanism to derive multi-view node feature embeddings and dynamically distribute propagation weights between each node and its neighbors, which facilitates the identification of significant biomarkers from functional brain networks and thus enhances the interpretability of the findings. To evaluate our proposed model in cognitive ability prediction, we conduct experiments on two large-scale brain imaging datasets: the Philadelphia Neurodevelopmental Cohort (PNC) and the Human Connectome Project (HCP). The results show that our approach not only outperforms existing methods in prediction accuracy but also provides superior explainability, which can be used to identify important FCs underlying cognitive behaviors.

Age Prediction Using Resting-State Functional MRI

The increasing lifespan and large individual differences in cognitive capability highlight the importance of comprehending the aging process of the brain. Contrary to visible signs of bodily ageing, like greying of hair and loss of muscle mass, the internal changes that occur within our brains remain less apparent until they impair function. Brain age, distinct from chronological age, reflects our brain's health status and may deviate from our actual chronological age. Notably, brain age has been associated with mortality and depression. The brain is plastic and can compensate even for severe structural damage by rewiring. Functional characterization offers insights that structural cannot provide. Contrary to the multitude of studies relying on structural magnetic resonance imaging (MRI), we utilize resting-state functional MRI (rsfMRI). We also address the issue of inclusion of subjects with abnormal brain ageing through outlier removal. In this study, we employ the Least Absolute Shrinkage and Selection Operator (LASSO) to identify the 39 most predictive correlations derived from the rsfMRI data. The data is from a cohort of 176 healthy right-handed volunteers, aged 18-78 years (95/81 male/female, mean age 48, SD 17) collected at the Mind Research Imaging Center at the National Cheng Kung University. We establish a normal reference model by excluding 68 outliers, which achieves a leave-one-out mean absolute error of 2.48 years. By asking which additional features that are needed to predict the chronological age of the outliers with a smaller error, we identify correlations predictive of abnormal aging. These are associated with the Default Mode Network (DMN). Our normal reference model has the lowest prediction error among published models evaluated on adult subjects of almost all ages and is thus a candidate for screening for abnormal brain aging that has not yet manifested in cognitive decline. This study advances our ability to predict brain aging and provides insights into potential biomarkers for assessing brain age, suggesting that the role of DMN in brain aging should be studied further.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Somatomotor-visual resting state functional connectivity increases after 2 years in the UK Biobank longitudinal cohort

Purpose

Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, in which high connectivity among all brain regions changes to a more modular structure with maturation. We examine FC changes in older adults after 2 years of aging in the UK Biobank (UKB) longitudinal cohort.

Approach

We process fMRI connectivity data using the Power264 atlas and then test whether the average internetwork FC changes in the 2722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t -test. We also compare the ability of Power264 and UKB-provided, independent component analysis (ICA)-based FC to determine which of a longitudinal scan pair is older. Finally, we investigate cross-sectional FC changes as well as differences due to differing scanner tasks in the UKB, Philadelphia Neurodevelopmental Cohort, and Alzheimer's Disease Neuroimaging Initiative datasets.

Results

We find a 6.8% average increase in somatomotor network (SMT)-visual network (VIS) connectivity from younger to older scans (corrected p < 10 - 15 ) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all internetwork connections, the average SMT-VIS connectivity is the best predictor of relative scan age. Using the full FC and a training set of 2000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC.

Conclusions

We conclude that SMT-VIS connectivity increases with age in the UKB longitudinal cohort and that resting state FC increases with age in the UKB cross-sectional cohort.

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.

Regression and Alignment for Functional Data and Network Topology

In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of pre-processing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.

Brain-phenotype predictions can survive across diverse real-world data

Recent work suggests that machine learning models predicting psychiatric treatment outcomes based on clinical data may fail when applied to unharmonized samples. Neuroimaging predictive models offer the opportunity to incorporate neurobiological information, which may be more robust to dataset shifts. Yet, among the minority of neuroimaging studies that undertake any form of external validation, there is a notable lack of attention to generalization across dataset-specific idiosyncrasies. 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 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. We demonstrate that reproducible and generalizable brain-behavior associations can be realized across diverse dataset features with sample sizes in the hundreds. Results indicate the potential of functional connectivity-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 neuroimaging predictive models in real-world scenarios and clinical settings.

Development and Maturation of the Human Brain, from Infancy to Adolescence

This chapter describes basic principles and key findings regarding the development and maturation of the human brain, the former referring to the pre-natal and early post-natal periods and the latter concerning childhood and adolescence. In both cases, we focus on brain structure as revealed in vivo with multi-modal magnetic resonance imaging (MRI). We begin with a few numbers about the human brain and its cellular composition and a brief overview of a number of MRI-based metrics used to characterize age-related variations in grey and white matter. We then proceed with synthesizing current knowledge about developmental and maturational changes in the cerebral cortex (its thickness, surface area, and intra-cortical myelination) and the underlying white matter (volume and structural properties). To facilitate biological interpretations of MRI-derived metrics, we introduce the concept of virtual histology. We conclude the chapter with a few notes about future directions in the study of factors shaping the human brain from conception onwards.

The effects of data leakage on connectome-based machine learning models

Predictive modeling has now become a central technique in neuroimaging to identify complex brain-behavior relationships and test their generalizability to unseen data. However, data leakage, which unintentionally breaches the separation between data used to train and test the model, undermines the validity of predictive models. Previous literature suggests that leakage is generally pervasive in machine learning, but few studies have empirically evaluated the effects of leakage in neuroimaging data. Although leakage is always an incorrect practice, understanding the effects of leakage on neuroimaging predictive models provides insight into the extent to which leakage may affect the literature. Here, we investigated the effects of leakage on machine learning models in two common neuroimaging modalities, functional and structural connectomes. Using over 400 different pipelines spanning four large datasets and three phenotypes, we evaluated five forms of leakage fitting into three broad categories: feature selection, covariate correction, and lack of independence between subjects. As expected, leakage via feature selection and repeated subjects drastically inflated prediction performance. Notably, other forms of leakage had only minor effects (e.g., leaky site correction) or even decreased prediction performance (e.g., leaky covariate regression). In some cases, leakage affected not only prediction performance, but also model coefficients, and thus neurobiological interpretations. Finally, we found that predictive models using small datasets were more sensitive to leakage. Overall, our results illustrate the variable effects of leakage on prediction pipelines and underscore the importance of avoiding data leakage to improve the validity and reproducibility of predictive modeling.

Development of white matter fiber covariance networks supports executive function in youth

During adolescence, the brain undergoes extensive changes in white matter structure that support cognition. Data-driven approaches applied to cortical surface properties have led the field to understand brain development as a spatially and temporally coordinated mechanism that follows hierarchically organized gradients of change. Although white matter development also appears asynchronous, previous studies have relied largely on anatomical tract-based atlases, precluding a direct assessment of how white matter structure is spatially and temporally coordinated. Harnessing advances in diffusion modeling and machine learning, we identified 14 data-driven patterns of covarying white matter structure in a large sample of youth. Fiber covariance networks aligned with known major tracts, while also capturing distinct patterns of spatial covariance across distributed white matter locations. Most networks showed age-related increases in fiber network properties, which were also related to developmental changes in executive function. This study delineates data-driven patterns of white matter development that support cognition.

Age-associated alterations in thalamocortical structural connectivity in youths with a psychosis-spectrum disorder

Psychotic symptoms typically emerge in adolescence. Age-associated thalamocortical connectivity differences in psychosis remain unclear. We analyzed diffusion-weighted imaging data from 1254 participants 8-23 years old (typically developing (TD):N = 626, psychosis-spectrum (PS): N = 329, other psychopathology (OP): N = 299) from the Philadelphia Neurodevelopmental Cohort. We modeled thalamocortical tracts using deterministic fiber tractography, extracted Q-Space Diffeomorphic Reconstruction (QSDR) and diffusion tensor imaging (DTI) measures, and then used generalized additive models to determine group and age-associated thalamocortical connectivity differences. Compared to other groups, PS exhibited thalamocortical reductions in QSDR global fractional anisotropy (GFA, p-values range = 3.0 × 10-6-0.05) and DTI fractional anisotropy (FA, p-values range = 4.2 × 10-4-0.03). Compared to TD, PS exhibited shallower thalamus-prefrontal age-associated increases in GFA and FA during mid-childhood, but steeper age-associated increases during adolescence. TD and OP exhibited decreases in thalamus-frontal mean and radial diffusivities during adolescence; PS did not. Altered developmental trajectories of thalamocortical connectivity may contribute to the disruptions observed in adults with psychosis.

ImageNomer: Description of a functional connectivity and omics analysis tool and case study identifying a race confound

Most packages for the analysis of fMRI-based functional connectivity (FC) and genomic data are used with a programming language interface, lacking an easy-to-navigate GUI frontend. This exacerbates two problems found in these types of data: demographic confounds and quality control in the face of high dimensionality of features. The reason is that it is too slow and cumbersome to use a programming interface to create all the necessary visualizations required to identify all correlations, confounding effects, or quality control problems in a dataset. FC in particular usually contains tens of thousands of features per subject, and can only be summarized and efficiently explored using visualizations. To remedy this situation, we have developed ImageNomer, a data visualization and analysis tool that allows inspection of both subject-level and cohort-level demographic, genomic, and imaging features. The software is Python-based, runs in a self-contained Docker image, and contains a browser-based GUI frontend. We demonstrate the usefulness of ImageNomer by identifying an unexpected race confound when predicting achievement scores in the Philadelphia Neurodevelopmental Cohort (PNC) dataset, which contains multitask fMRI and single nucleotide polymorphism (SNP) data of healthy adolescents. In the past, many studies have attempted to use FC to identify achievement-related features in fMRI. Using ImageNomer to visualize trends in achievement scores between races, we find a clear potential for confounding effects if race can be predicted using FC. Using correlation analysis in the ImageNomer software, we show that FCs correlated with Wide Range Achievement Test (WRAT) score are in fact more highly correlated with race. Investigating further, we find that whereas both FC and SNP (genomic) features can account for 10-15% of WRAT score variation, this predictive ability disappears when controlling for race. We also use ImageNomer to investigate race-FC correlation in the Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP) dataset. In this work, we demonstrate the advantage of our ImageNomer GUI tool in data exploration and confound detection. Additionally, this work identifies race as a strong confound in FC data and casts doubt on the possibility of finding unbiased achievement-related features in fMRI and SNP data of healthy adolescents.

Pathways from threat exposure to psychotic symptoms in youth: The role of emotion recognition bias and brain structure

BACKGROUND

Research supports an association between threatening experiences in childhood and psychosis. It is possible that early threat exposure disrupts the development of emotion recognition (specifically, producing a bias for facial expressions relating to threat) and the brain structures subserving it, contributing to psychosis development.

METHODS

Using data from the Philadelphia Neurodevelopmental Cohort, we examined associations between threat exposure and both the misattribution of facial expressions to fear/anger in an emotion recognition task, and gray matter volumes in key emotion processing regions. Our sample comprised youth with psychosis spectrum symptoms (N = 304), control youth (N = 787), and to evaluate specificity, youth with internalizing symptoms (N = 92). The moderating effects of group and sex were examined.

RESULTS

Both the psychosis spectrum and internalizing groups had higher levels of threat exposure than controls. In the total sample, threat exposure was associated with lower left medial prefrontal cortex (mPFC) volume but not misattributions to fear/anger. The effects of threat exposure did not significantly differ by group or sex.

CONCLUSIONS

The findings of this study provide evidence for an effect of threat exposure on mPFC morphology, but do not support an association between threat exposure and a recognition bias for threat-related expressions, that is particularly pronounced in psychosis. Future research should investigate factors linking transdiagnostic alterations related to threat exposure with psychotic symptoms, and attempt to clarify the mechanisms underpinning emotion recognition misattributions in threat-exposed youth.

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

Identifying reproducible and generalizable brain-phenotype associations is a central goal of neuroimaging. Consistent with this goal, prediction frameworks evaluate brain-phenotype models in unseen data. Most prediction studies train and evaluate a model in the same dataset. However, external validation, or the evaluation of a model in an external dataset, provides a better assessment of robustness and generalizability. Despite the promise of external validation and calls for its usage, the statistical power of such studies has yet to be investigated. In this work, we ran over 60 million simulations across several datasets, phenotypes, and sample sizes to better understand how the sizes of the training and external datasets affect statistical power. We found that prior external validation studies used sample sizes prone to low power, which may lead to false negatives and effect size inflation. Furthermore, increases in the external sample size led to increased simulated power directly following theoretical power curves, whereas changes in the training dataset size offset the simulated power curves. Finally, we compared the performance of a model within a dataset to the external performance. The within-dataset performance was typically within r=0.2 of the cross-dataset performance, which could help decide how to power future external validation studies. Overall, our results illustrate the importance of considering the sample sizes of both the training and external datasets when performing external validation.

White matter functional gradients and their formation in adolescence

It is well known that functional magnetic resonance imaging (fMRI) is a widely used tool for studying brain activity. Recent research has shown that fluctuations in fMRI data can reflect functionally meaningful patterns of brain activity within the white matter. We leveraged resting-state fMRI from an adolescent population to characterize large-scale white matter functional gradients and their formation during adolescence. The white matter showed gray-matter-like unimodal-to-transmodal and sensorimotor-to-visual gradients with specific cognitive associations and a unique superficial-to-deep gradient with nonspecific cognitive associations. We propose two mechanisms for their formation in adolescence. One is a "function-molded" mechanism that may mediate the maturation of the transmodal white matter via the transmodal gray matter. The other is a "structure-root" mechanism that may support the mutual mediation roles of the unimodal and transmodal white matter maturation during adolescence. Thus, the spatial layout of the white matter functional gradients is in concert with the gray matter functional organization. The formation of the white matter functional gradients may be driven by brain anatomical wiring and functional needs.

A systematic review and meta-analysis of resting-state fMRI in anxiety disorders: Need for data sharing to move the field forward

Anxiety disorders are among the most prevalent psychiatric disorders. Neuroimaging findings remain uncertain, and resting state functional magnetic resonance (rs-fMRI) connectivity is of particular interest since it is a scalable functional imaging modality. Given heterogeneous past findings for rs-fMRI in anxious individuals, we characterize patterns across anxiety disorders by conducting a systematic review and meta-analysis. Studies were included if they contained at the time of scanning both a healthy group and a patient group. Due to insufficient study numbers, the quantitative meta-analysis only included seed-based studies. We performed an activation likelihood estimation (ALE) analysis that compared patients and healthy volunteers. All analyses were corrected for family-wise error with a cluster-level threshold of p < .05. Patients exhibited hypo-connectivity between the amygdala and the medial frontal gyrus, anterior cingulate cortex, and cingulate gyrus. This finding, however, was not robust to potential file-drawer effects. Though limited by strict inclusion criteria, our results highlight the heterogeneous nature of reported findings. This underscores the need for data sharing when attempting to detect reliable patterns of disruption in brain activity across anxiety disorders.

Transmit Radiant Individuality to Offspring (TRIO) study: investigating intergenerational transmission effects on brain development

Intergenerational transmission is a crucial aspect of human development. Although prior studies have demonstrated the continuity of psychopathology and maladaptive upbringing environments between parents and offspring, the underlying neurobiological mechanisms remain unclear. We have begun a novel neuroimaging research project, the Transmit Radiant Individuality to Offspring (TRIO) study, which focuses on biological parent-offspring trios. The participants of the TRIO study were Japanese parent-offspring trios consisting of offspring aged 10-40 and their biological mother and father. Structural and functional brain images of all participants were acquired using magnetic resonance imaging (MRI). Saliva samples were collected for DNA analysis. We obtained psychosocial information, such as intelligence, mental health problems, personality traits, and experiences during the developmental period from each parent and offspring in the same manner as much as possible. By April 2023, we completed data acquisition from 174 trios consisting of fathers, mothers, and offspring. The target sample size was 310 trios. However, we plan to conduct genetic and epigenetic analyses, and the sample size is expected to be expanded further while developing this project into a multi-site collaborative study in the future. The TRIO study can challenge the elucidation of the mechanism of intergenerational transmission effects on human development by collecting diverse information from parents and offspring at the molecular, neural, and behavioral levels. Our study provides interdisciplinary insights into how individuals' lives are involved in the construction of the lives of their descendants in the subsequent generation.