The Philadelphia Neurodevelopmental Cohort: A publicly available resource for the study of normal and abnormal brain development in youth

MULTIMODAL NEURAL CORRELATES OF CHILDHOOD PSYCHOPATHOLOGY

Complex structural and functional changes occurring in typical and atypical development necessitate multidimensional approaches to better understand the risk of developing psychopathology. Here, we simultaneously examined structural and functional brain network patterns in relation to dimensions of psychopathology in the Adolescent Brain Cognitive Development dataset. Several components were identified, recapitulating the psychopathology hierarchy, with the general psychopathology (p) factor explaining most covariance with multimodal imaging features, while the internalizing, externalizing, and neurodevelopmental dimensions were each associated with distinct morphological and functional connectivity signatures. Connectivity signatures associated with the p factor and neurodevelopmental dimensions followed the sensory-to-transmodal axis of cortical organization, which is related to the emergence of complex cognition and risk for psychopathology. Results were consistent in two separate data subsamples, supporting generalizability, and robust to variations in analytical parameters. Our findings help in better understanding biological mechanisms underpinning dimensions of psychopathology, and could provide brain-based vulnerability markers.

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.

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 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.

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 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.

An Automated Tool to Classify and Transform Unstructured MRI Data into BIDS Datasets

The increasing use of neuroimaging in clinical research has driven the creation of many large imaging datasets. However, these datasets often rely on inconsistent naming conventions in image file headers to describe acquisition, and time-consuming manual curation is necessary. Therefore, we sought to automate the process of classifying and organizing magnetic resonance imaging (MRI) data according to acquisition types common to the clinical routine, as well as automate the transformation of raw, unstructured images into Brain Imaging Data Structure (BIDS) datasets. To do this, we trained an XGBoost model to classify MRI acquisition types using relatively few acquisition parameters that are automatically stored by the MRI scanner in image file metadata, which are then mapped to the naming conventions prescribed by BIDS to transform the input images to the BIDS structure. The model recognizes MRI types with 99.475% accuracy, as well as a micro/macro-averaged precision of 0.9995/0.994, a micro/macro-averaged recall of 0.9995/0.989, and a micro/macro-averaged F1 of 0.9995/0.991. Our approach accurately and quickly classifies MRI types and transforms unstructured data into standardized structures with little-to-no user intervention, reducing the barrier of entry for clinical scientists and increasing the accessibility of existing neuroimaging data.

Mapping Individual Differences in Intermodal Coupling in Neurodevelopment

Within-individual coupling between measures of brain structure and function evolves in development and may underlie differential risk for neuropsychiatric disorders. Despite increasing interest in the development of structure-function relationships, rigorous methods to quantify and test individual differences in coupling remain nascent. In this article, we explore and address gaps in approaches for testing and spatially localizing individual differences in intermodal coupling. We propose a new method, called CIDeR, which is designed to simultaneously perform hypothesis testing in a way that limits false positive results and improve detection of true positive results. Through a comparison across different approaches to testing individual differences in intermodal coupling, we delineate subtle differences in the hypotheses they test, which may ultimately lead researchers to arrive at different results. Finally, we illustrate the utility of CIDeR in two applications to brain development using data from the Philadelphia Neurodevelopmental Cohort.

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.

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.

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

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 Developmental Perspective on Early and Current Motor Abnormalities and Psychotic-Like Symptoms

BACKGROUND AND HYPOTHESIS

Psychotic-like experiences (PLEs) are prevalent in the general population and, because they represent a lower end of the psychosis vulnerability spectrum, may be useful in informing mechanistic understanding. Although it is well-understood that motor signs characterize formal psychotic disorders, the developmental trajectory of these features and their relationships with PLEs are less well-understood.

STUDY DESIGN

Data from 7559 adolescents and young adults (age 11-21) in the Philadelphia Neurodevelopmental Cohort were used to investigate whether early-life milestone-attainment delays relate to current adolescent sensorimotor functioning and positive and negative PLEs. Current sensorimotor functioning was assessed using the Computerized Finger Tapping task (assessing motor slowing) and Mouse Practice task (assessing sensorimotor planning).

STUDY RESULTS

Early developmental abnormalities were related to current adolescent-aged motor slowing (t(7415.3) = -7.74, corrected-P < .001) and impaired sensorimotor planning (t(7502.5) = 5.57, corrected-P < .001). There was a significant interaction between developmental delays and current sensorimotor functioning on positive and negative PLEs (t = 1.67-4.51), such that individuals with early developmental delays had a stronger positive relationship between sensorimotor dysfunction and PLEs. Importantly, interaction models were significantly better at explaining current PLEs than those treating early and current sensorimotor dysfunction independently (χ2 = 4.89-20.34).

CONCLUSIONS

These findings suggest a relationship between early developmental delays and current sensorimotor functioning in psychosis proneness and inform an understanding of heterotypic continuity as well as a neurodevelopmental perspective of motor circuits. Furthermore, results indicate that motor signs are a clear factor in the psychosis continuum, suggesting that they may represent a core feature of psychosis vulnerability.

Reliability and predictability of phenotype information from functional connectivity in large imaging datasets

One of the central objectives of contemporary neuroimaging research is to create predictive models that can disentangle the connection between patterns of functional connectivity across the entire brain and various behavioral traits. Previous studies have shown that models trained to predict behavioral features from the individual's functional connectivity have modest to poor performance. In this study, we trained models that predict observable individual traits (phenotypes) and their corresponding singular value decomposition (SVD) representations - herein referred to as latent phenotypes from resting state functional connectivity. For this task, we predicted phenotypes in two large neuroimaging datasets: the Human Connectome Project (HCP) and the Philadelphia Neurodevelopmental Cohort (PNC). We illustrate the importance of regressing out confounds, which could significantly influence phenotype prediction. Our findings reveal that both phenotypes and their corresponding latent phenotypes yield similar predictive performance. Interestingly, only the first five latent phenotypes were reliably identified, and using just these reliable phenotypes for predicting phenotypes yielded a similar performance to using all latent phenotypes. This suggests that the predictable information is present in the first latent phenotypes, allowing the remainder to be filtered out without any harm in performance. This study sheds light on the intricate relationship between functional connectivity and the predictability and reliability of phenotypic information, with potential implications for enhancing predictive modeling in the realm of neuroimaging research.

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.

Menarche, pubertal timing and the brain: female-specific patterns of brain maturation beyond age-related development

BACKGROUND

Puberty depicts a period of profound and multifactorial changes ranging from social to biological factors. While brain development in youths has been studied mostly from an age perspective, recent evidence suggests that pubertal measures may be more sensitive to study adolescent neurodevelopment, however, studies on pubertal timing in relation to brain development are still scarce.

METHODS

We investigated if pre- vs. post-menarche status can be classified using machine learning on cortical and subcortical structural magnetic resonance imaging (MRI) data from strictly age-matched adolescent females from the Adolescent Brain Cognitive Development (ABCD) cohort. For comparison of the identified menarche-related patterns to age-related patterns of neurodevelopment, we trained a brain age prediction model on data from the Philadelphia Neurodevelopmental Cohort and applied it to the same ABCD data, yielding differences between predicted and chronological age referred to as brain age gaps. We tested the sensitivity of both these frameworks to measures of pubertal maturation, specifically age at menarche and puberty status.

RESULTS

The machine learning model achieved moderate but statistically significant accuracy in the menarche classification task, yielding for each subject a class probability ranging from 0 (pre-) to 1 (post- menarche). Comparison to brain age predictions revealed shared and distinct patterns of neurodevelopment captured by both approaches. Continuous menarche class probabilities were positively associated with brain age gaps, but only the menarche class probabilities-not the brain age gaps-were associated with age at menarche.

CONCLUSIONS

This study demonstrates the use of a machine learning model to classify menarche status from structural MRI data while accounting for age-related neurodevelopment. Given its sensitivity towards measures of puberty timing, our work suggests that menarche class probabilities may be developed toward an objective brain-based marker of pubertal development.

The future of neuroscience in developmental psychopathology

Developmental psychopathology started as an intersection of fields and is now a field itself. As we contemplate the future of this field, we consider the ways in which a newer, interdisciplinary field - human developmental neuroscience - can inform, and be informed by, developmental psychopathology. To do so, we outline principles of developmental psychopathology and how they are and/or can be implemented in developmental neuroscience. In turn, we highlight how the collaboration between these fields can lead to richer models and more impactful translation. In doing so, we describe the ways in which models from developmental psychopathology can enrich developmental neuroscience and future directions for developmental psychopathology.

Neuroimaging and Biosample Collection in the Toronto Adolescent and Youth Cohort Study: Rationale, Methods, and Early Data

BACKGROUND

The Toronto Adolescent and Youth (TAY) Cohort Study will characterize the neurobiological trajectories of psychosis spectrum symptoms, functioning, and suicidality (i.e., suicidal thoughts and behaviors) in youth seeking mental health care. Here, we present the neuroimaging and biosample component of the protocol. We also present feasibility and quality control metrics for the baseline sample collected thus far.

METHODS

The current study includes youths (ages 11-24 years) who were referred to child and youth mental health services within a large tertiary care center in Toronto, Ontario, Canada, with target recruitment of 1500 participants. Participants were offered the opportunity to provide any or all of the following: 1) 1-hour magnetic resonance imaging (MRI) scan (electroencephalography if ineligible for or declined MRI), 2) blood sample for genomic and proteomic data (or saliva if blood collection was declined or not feasible) and urine sample, and 3) heart rate recording to assess respiratory sinus arrhythmia.

RESULTS

Of the first 417 participants who consented to participate between May 4, 2021, and February 2, 2023, 412 agreed to participate in the imaging and biosample protocol. Of these, 334 completed imaging, 341 provided a biosample, 338 completed respiratory sinus arrhythmia, and 316 completed all 3. Following quality control, data usability was high (MRI: T1-weighted 99%, diffusion-weighted imaging 99%, arterial spin labeling 90%, resting-state functional MRI 95%, task functional MRI 90%; electroencephalography: 83%; respiratory sinus arrhythmia: 99%).

CONCLUSIONS

The high consent rates, good completion rates, and high data usability reported here demonstrate the feasibility of collecting and using brain imaging and biosamples in a large clinical cohort of youths seeking mental health care.

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.

Association between suicidal ideation and tandem repeats in contactins

Background

Death by suicide is one of the leading causes of death among adolescents. Genome-wide association studies (GWAS) have identified loci that associate with suicidal ideation and related behaviours. One such group of loci are the six contactin genes (CNTN1-6) that are critical to neurodevelopment through regulating neurite structure. Because single nucleotide polymorphisms (SNPs) detected by GWAS often map to non-coding intergenic regions, we investigated whether repetitive variants in CNTNs associated with suicidality in a young cohort aged 8 to 21. Understanding the genetic liability of suicidal thought and behavior in this age group will promote early intervention and treatment.

Methods

Genotypic and phenotypic data were obtained from the Philadelphia Neurodevelopment Cohort (PNC). Across six CNTNs, 232 short tandem repeats (STRs) were analyzed in up to 4,595 individuals of European ancestry who expressed current, previous, or no suicidal ideation. STRs were imputed into SNP arrays using a phased SNP-STR haplotype reference panel from the 1000 Genomes Project. We tested several additive and interactive models of locus-level burden (i.e., sum of STR alleles) with respect to suicidal ideation. Additive models included sex, birth year, developmental stage ("DevStage"), and the first 10 principal components of ancestry as covariates; interactive models assessed the effect of STR-by-DevStage considering all other covariates.

Results

CNTN1-[T]N interacted with DevStage to increase risk for current suicidal ideation (CNTN1-[T]N-by-DevStage; p = 0.00035). Compared to the youngest age group, the middle (OR = 1.80, p = 0.0514) and oldest (OR = 3.82, p = 0.0002) participant groups had significantly higher odds of suicidal ideation as their STR length expanded; this result was independent of polygenic scores for suicidal ideation.

Discussion

These findings highlight diversity in the genetic effects (i.e., SNP and STR) acting on suicidal thoughts and behavior and advance our understanding of suicidal ideation across childhood and adolescence.

The general factor of psychopathology (p): Choosing among competing models and interpreting p

Over the past 10 years, the general factor of psychopathology, p, has attracted interest and scrutiny. We review the history of the idea that all mental disorders share something in common, p; how we arrived at this idea; and how it became conflated with a statistical representation, the Bi-Factor Model. We then leverage the Environmental Risk (E-Risk) longitudinal twin study to examine the properties and nomological network of different statistical representations of p. We find that p performed similarly regardless of how it was modelled, suggesting that if the sample and content are the same the resulting p factor will be similar. We suggest that the meaning of p is not to be found by dueling over statistical models but by conducting well-specified criterion-validation studies and developing new measurement approaches. We outline new directions to refresh research efforts to uncover what all mental disorders have in common.

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.

CoCoA: conditional correlation models with association size

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

Transdiagnostic biomarkers of mental illness across the lifespan: A systematic review examining the genetic and neural correlates of latent transdiagnostic dimensions of psychopathology in the general population

This systematic review synthesizes evidence from research investigating the biological correlates of latent transdiagnostic dimensions of psychopathology (e.g., the p-factor, internalizing, externalizing) across the lifespan. Eligibility criteria captured genomic and neuroimaging studies investigating general and/or specific dimensions in general population samples across all age groups. MEDLINE, Embase, and PsycINFO were searched for relevant studies published up to March 2023 and 46 studies were selected for inclusion. The results revealed several biological correlates consistently associated with transdiagnostic dimensions of psychopathology, including polygenic scores for ADHD and neuroticism, global surface area and global gray matter volume. Shared and unique associations between symptom dimensions are highlighted, as are potential age-specific differences in biological associations. Findings are interpreted with reference to key methodological differences across studies. The included studies provide compelling evidence that the general dimension of psychopathology reflects common underlying genetic and neurobiological vulnerabilities that are shared across diverse manifestations of mental illness. Substantive interpretations of general psychopathology in the context of genetic and neurobiological evidence are discussed.

XCP-D: A Robust Pipeline for the post-processing of fMRI data

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

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.

Domain adapted brain network fusion captures variance related to pubertal brain development and mental health

Puberty demarks a period of profound brain dynamics that orchestrates changes to a multitude of neuroimaging-derived phenotypes. This complexity poses a dimensionality problem when attempting to chart an individual's brain development over time. Here, we illustrate that shifts in subject similarity of brain imaging data relate to pubertal maturation in the longitudinal ABCD study. Given that puberty depicts a critical window for emerging mental health issues, we additionally show that our model is capable of capturing variance in the adolescent brain related to psychopathology in a population-based and a clinical cohort. These results suggest that low-dimensional reference spaces based on subject similarities render useful to chart variance in brain development in youths.

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.

Automated analysis of low-field brain MRI in cerebral malaria

A central challenge of medical imaging studies is to extract biomarkers that characterize disease pathology or outcomes. Modern automated approaches have found tremendous success in high-resolution, high-quality magnetic resonance images. These methods, however, may not translate to low-resolution images acquired on magnetic resonance imaging (MRI) scanners with lower magnetic field strength. In low-resource settings where low-field scanners are more common and there is a shortage of radiologists to manually interpret MRI scans, it is critical to develop automated methods that can augment or replace manual interpretation, while accommodating reduced image quality. We present a fully automated framework for translating radiological diagnostic criteria into image-based biomarkers, inspired by a project in which children with cerebral malaria (CM) were imaged using low-field 0.35 Tesla MRI. We integrate multiatlas label fusion, which leverages high-resolution images from another sample as prior spatial information, with parametric Gaussian hidden Markov models based on image intensities, to create a robust method for determining ventricular cerebrospinal fluid volume. We also propose normalized image intensity and texture measurements to determine the loss of gray-to-white matter tissue differentiation and sulcal effacement. These integrated biomarkers have excellent classification performance for determining severe brain swelling due to CM.

Harmonizing bifactor models of psychopathology between distinct assessment instruments: Reliability, measurement invariance, and authenticity

OBJECTIVES

Model configuration is important for mental health data harmonization. We provide a method to investigate the performance of different bifactor model configurations to harmonize different instruments.

METHODS

We used data from six samples from the Reproducible Brain Charts initiative (N = 8,606, ages 5-22 years, 41.0% females). We harmonized items from two psychopathology instruments, Child Behavior Checklist (CBCL) and GOASSESS, based on semantic content. We estimated bifactor models using confirmatory factor analysis, and calculated their model fit, factor reliability, between-instrument invariance, and authenticity (i.e., the correlation and factor score difference between the harmonized and original models).

RESULTS

Five out of 12 model configurations presented acceptable fit and were instrument-invariant. Correlations between the harmonized factor scores and the original full-item models were high for the p-factor (>0.89) and small to moderate (0.12-0.81) for the specific factors. 6.3%-50.9% of participants presented factor score differences between harmonized and original models higher than 0.5 z-score.

CONCLUSIONS

The CBCL-GOASSESS harmonization indicates that few models provide reliable specific factors and are instrument-invariant. Moreover, authenticity was high for the p-factor and moderate for specific factors. Future studies can use this framework to examine the impact of harmonizing instruments in psychiatric research.

Genome-wide analysis of a model-derived binge eating disorder phenotype identifies risk loci and implicates iron metabolism

Binge eating disorder (BED) is the most common eating disorder, yet its genetic architecture remains largely unknown. Studying BED is challenging because it is often comorbid with obesity, a common and highly polygenic trait, and it is underdiagnosed in biobank data sets. To address this limitation, we apply a supervised machine-learning approach (using 822 cases of individuals diagnosed with BED) to estimate the probability of each individual having BED based on electronic medical records from the Million Veteran Program. We perform a genome-wide association study of individuals of African (n = 77,574) and European (n = 285,138) ancestry while controlling for body mass index to identify three independent loci near the HFE, MCHR2 and LRP11 genes and suggest APOE as a risk gene for BED. We identify shared heritability between BED and several neuropsychiatric traits, and implicate iron metabolism in the pathophysiology of BED. Overall, our findings provide insights into the genetics underlying BED and suggest directions for future translational research.

Neurogenetic mechanisms of risk for ADHD: Examining associations of polygenic scores and brain volumes in a population cohort

BACKGROUND

ADHD polygenic scores (PGSs) have been previously shown to predict ADHD outcomes in several studies. However, ADHD PGSs are typically correlated with ADHD but not necessarily reflective of causal mechanisms. More research is needed to elucidate the neurobiological mechanisms underlying ADHD. We leveraged functional annotation information into an ADHD PGS to (1) improve the prediction performance over a non-annotated ADHD PGS and (2) test whether volumetric variation in brain regions putatively associated with ADHD mediate the association between PGSs and ADHD outcomes.

METHODS

Data were from the Philadelphia Neurodevelopmental Cohort (N = 555). Multiple mediation models were tested to examine the indirect effects of two ADHD PGSs-one using a traditional computation involving clumping and thresholding and another using a functionally annotated approach (i.e., AnnoPred)-on ADHD inattention (IA) and hyperactivity-impulsivity (HI) symptoms, via gray matter volumes in the cingulate gyrus, angular gyrus, caudate, dorsolateral prefrontal cortex (DLPFC), and inferior temporal lobe.

RESULTS

A direct effect was detected between the AnnoPred ADHD PGS and IA symptoms in adolescents. No indirect effects via brain volumes were detected for either IA or HI symptoms. However, both ADHD PGSs were negatively associated with the DLPFC.

CONCLUSIONS

The AnnoPred ADHD PGS was a more developmentally specific predictor of adolescent IA symptoms compared to the traditional ADHD PGS. However, brain volumes did not mediate the effects of either a traditional or AnnoPred ADHD PGS on ADHD symptoms, suggesting that we may still be underpowered in clarifying brain-based biomarkers for ADHD using genetic measures.

Predictors of transition in patients with clinical high risk for psychosis: an umbrella review

Diagnosis of a clinical high-risk (CHR) state enables timely treatment of individuals at risk for a psychotic disorder, thereby contributing to improving illness outcomes. However, only a minority of patients diagnosed with CHR will make the transition to overt psychosis. To identify patients most likely to benefit from early intervention, several studies have investigated characteristics that distinguish CHR patients who will later develop a psychotic disorder from those who will not. We aimed to summarize evidence from systematic reviews and meta-analyses on predictors of transition to psychosis in CHR patients, among characteristics and biomarkers assessed at baseline. A systematic search was conducted in Pubmed, Scopus, PsychInfo and Cochrane databases to identify reviews and meta-analyses of studies that investigated specific baseline predictors or biomarkers for transition to psychosis in CHR patients using a cross-sectional or longitudinal design. Non-peer-reviewed publications, gray literature, narrative reviews and publications not written in English were excluded from analyses. We provide a narrative synthesis of results from all included reviews and meta-analyses. For each included publication, we indicate the number of studies cited in each domain and its quality rating. A total of 40 publications (21 systematic reviews and 19 meta-analyses) that reviewed a total of 272 original studies qualified for inclusion. Baseline predictors most consistently associated with later transition included clinical characteristics such as attenuated psychotic and negative symptoms and functioning, verbal memory deficits and the electrophysiological marker of mismatch negativity. Few predictors reached a level of evidence sufficient to inform clinical practice, reflecting generalizability issues in a field characterized by studies with small, heterogeneous samples and relatively few transition events. Sample pooling and harmonization of methods across sites and projects are necessary to overcome these limitations.

Using network control theory to study the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. 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 dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Replication and Refinement of Brain Age Model for adolescent development

The discrepancy between chronological age and estimated brain age, known as the brain age gap, may serve as a biomarker to reveal brain development and neuropsychiatric problems. This has motivated many studies focusing on the accurate estimation of brain age using different features and models, of which the generalizability is yet to be tested. Our recent study has demonstrated that conventional machine learning models can achieve high accuracy on brain age prediction during development using only a small set of selected features from multimodal brain imaging data. In the current study, we tested the replicability of various brain age models on the Adolescent Brain Cognitive Development (ABCD) cohort. We proposed a new refined model to improve the robustness of brain age prediction. The direct replication test for existing brain age models derived from the age range of 8-22 years onto the ABCD participants at baseline (9 to 10 years old) and year-two follow-up (11 to 12 years old) indicate that pre-trained models could capture the overall mean age failed precisely estimating brain age variation within a narrow range. The refined model, which combined broad prediction of the pre-trained model and granular information with the narrow age range, achieved the best performance with a mean absolute error of 0.49 and 0.48 years on the baseline and year-two data, respectively. The brain age gap yielded by the refined model showed significant associations with the participants' information processing speed and verbal comprehension ability on baseline data.

Connectome-based machine learning models are vulnerable to subtle data manipulations

Neuroimaging-based predictive models continue to improve in performance, yet a widely overlooked aspect of these models is "trustworthiness," or robustness to data manipulations. High trustworthiness is imperative for researchers to have confidence in their findings and interpretations. In this work, we used functional connectomes to explore how minor data manipulations influence machine learning predictions. These manipulations included a method to falsely enhance prediction performance and adversarial noise attacks designed to degrade performance. Although these data manipulations drastically changed model performance, the original and manipulated data were extremely similar (r = 0.99) and did not affect other downstream analysis. Essentially, connectome data could be inconspicuously modified to achieve any desired prediction performance. Overall, our enhancement attacks and evaluation of existing adversarial noise attacks in connectome-based models highlight the need for counter-measures that improve the trustworthiness to preserve the integrity of academic research and any potential translational applications.

Novel effects of Ras-MAPK pathogenic variants on the developing human brain and their link to gene expression and inhibition abilities

The RASopathies are genetic syndromes associated with pathogenic variants causing dysregulation of the Ras/mitogen-activated protein kinase (Ras-MAPK) pathway, essential for brain development, and increased risk for neurodevelopmental disorders. Yet, the effects of most pathogenic variants on the human brain are unknown. We examined: (1) How Ras-MAPK activating variants of PTPN11/SOS1 protein-coding genes affect brain anatomy. (2) The relationship between PTPN11 gene expression levels and brain anatomy, and (3) The relevance of subcortical anatomy to attention and memory skills affected in the RASopathies. We collected structural brain MRI and cognitive-behavioral data from 40 pre-pubertal children with Noonan syndrome (NS), caused by PTPN11 (n = 30) or SOS1 (n = 10) variants (age 8.53 ± 2.15, 25 females), and compared them to 40 age- and sex-matched typically developing controls (9.24 ± 1.62, 27 females). We identified widespread effects of NS on cortical and subcortical volumes and on determinants of cortical gray matter volume, surface area (SA), and cortical thickness (CT). In NS, we observed smaller volumes of bilateral striatum, precentral gyri, and primary visual area (d's < -0.8), and extensive effects on SA (d's > |0.8|) and CT (d's > |0.5|) relative to controls. Further, SA effects were associated with increasing PTPN11 gene expression, most prominently in the temporal lobe. Lastly, PTPN11 variants disrupted normative relationships between the striatum and inhibition functioning. We provide evidence for the effects of Ras-MAPK pathogenic variants on striatal and cortical anatomy as well as links between PTPN11 gene expression and cortical SA increases, and striatal volume and inhibition skills. These findings provide essential translational information on the Ras-MAPK pathway's effect on human brain development and function.

Systematic comparisons of different quality control approaches applied to three large pediatric neuroimaging datasets

INTRODUCTION

Poor quality T1-weighted brain scans systematically affect the calculation of brain measures. Removing the influence of such scans requires identifying and excluding scans with noise and artefacts through a quality control (QC) procedure. While QC is critical for brain imaging analyses, it is not yet clear whether different QC approaches lead to the exclusion of the same participants. Further, the removal of poor-quality scans may unintentionally introduce a sampling bias by excluding the subset of participants who are younger and/or feature greater clinical impairment. This study had two aims: 1) examine whether different QC approaches applied to T1-weighted scans would exclude the same participants, and 2) examine how exclusion of poor-quality scans impacts specific demographic, clinical and brain measure characteristics between excluded and included participants in three large pediatric neuroimaging samples.

METHODS

We used T1-weighted, resting-state fMRI, demographic and clinical data from the Province of Ontario Neurodevelopmental Disorders network (Aim 1: n=553, Aim 2: n=465), the Healthy Brain Network (Aim 1: n=1051, Aim 2: n=558), and the Philadelphia Neurodevelopmental Cohort (Aim 1: n=1087; Aim 2: n=619). Four different QC approaches were applied to T1-weighted MRI (visual QC, metric QC, automated QC, fMRI-derived QC). We used tetrachoric correlation and inter-rater reliability analyses to examine whether different QC approaches excluded the same participants. We examined differences in age, mental health symptoms, everyday/adaptive functioning, IQ and structural MRI-derived brain indices between participants that were included versus excluded following each QC approach.

RESULTS

Dataset-specific findings revealed mixed results with respect to overlap of QC exclusion. However, in POND and HBN, we found a moderate level of overlap between visual and automated QC approaches (r_tet=0.52-0.59). Implementation of QC excluded younger participants, and tended to exclude those with lower IQ, and lower everyday/adaptive functioning scores across several approaches in a dataset-specific manner. Across nearly all datasets and QC approaches examined, excluded participants had lower estimates of cortical thickness and subcortical volume, but this effect did not differ by QC approach.

CONCLUSION

The results of this study provide insight into the influence of QC decisions on structural pediatric imaging analyses. While different QC approaches exclude different subsets of participants, the variation of influence of different QC approaches on clinical and brain metrics is minimal in large datasets. Overall, implementation of QC tends to exclude participants who are younger, and those who have more cognitive and functional impairment. Given that automated QC is standardized and can reduce between-study differences, the results of this study support the potential to use automated QC for large pediatric neuroimaging datasets.

Novel effects of Ras-MAPK pathogenic variants on the developing human brain and their link to gene expression and inhibition abilities

The RASopathies are genetic syndromes associated with pathogenic variants causing dysregulation of the Ras/mitogen-activated protein kinase (Ras-MAPK) pathway, essential for brain development, and increased risk for neurodevelopmental disorders. Yet, the effects of most pathogenic variants on the human brain are unknown. We examined: 1. How Ras-MAPK activating variants of PTPN11 / SOS1 protein-coding genes affect brain anatomy. 2. The relationship between PTPN11 gene expression levels and brain anatomy, and 3. The relevance of subcortical anatomy to attention and memory skills affected in the RASopathies. We collected structural brain MRI and cognitive-behavioral data from 40 pre-pubertal children with Noonan syndrome (NS), caused by PTPN11 ( n  = 30) or SOS1 ( n  = 10) variants (age 8.53 ± 2.15, 25 females), and compared them to 40 age- and sex-matched typically developing controls (9.24 ± 1.62, 27 females). We identified widespread effects of NS on cortical and subcortical volumes and on determinants of cortical gray matter volume, surface area (SA) and cortical thickness (CT). In NS, we observed smaller volumes of bilateral striatum, precentral gyri, and primary visual area ( d 's<-0.8), and extensive effects on SA ( d 's>|0.8|) and CT ( d 's>|0.5|) relative to controls. Further, SA effects were associated with increasing PTPN11 gene expression, most prominently in the temporal lobe. Lastly, PTPN11 variants disrupted normative relationships between the striatum and inhibition functioning. We provide evidence for effects of Ras-MAPK pathogenic variants on striatal and cortical anatomy as well as links between PTPN11 gene expression and cortical SA increases, and striatal volume and inhibition skills. These findings provide essential translational information on the Ras-MAPK pathway's effect on human brain development and function.

Shared pattern of impaired social communication and cognitive ability in the youth brain across diagnostic boundaries

BACKGROUND

Abnormalities in brain structure are shared across diagnostic categories. Given the high rate of comorbidity, the interplay of relevant behavioural factors may also cross these classic boundaries.

METHODS

We aimed to detect brain-based dimensions of behavioural factors using canonical correlation and independent component analysis in a clinical youth sample (n = 1732, 64 % male, age: 5-21 years).

RESULTS

We identified two correlated patterns of brain structure and behavioural factors. The first mode reflected physical and cognitive maturation (r = 0.92, p = .005). The second mode reflected lower cognitive ability, poorer social skills, and psychological difficulties (r = 0.92, p = .006). Elevated scores on the second mode were a common feature across all diagnostic boundaries and linked to the number of comorbid diagnoses independently of age. Critically, this brain pattern predicted normative cognitive deviations in an independent population-based sample (n = 1253, 54 % female, age: 8-21 years), supporting the generalisability and external validity of the reported brain-behaviour relationships.

CONCLUSIONS

These results reveal dimensions of brain-behaviour associations across diagnostic boundaries, highlighting potent disorder-general patterns as the most prominent. In addition to providing biologically informed patterns of relevant behavioural factors for mental illness, this contributes to a growing body of evidence in favour of transdiagnostic approaches to prevention and intervention.

General psychopathology factor (p-factor) prediction using resting-state functional connectivity and a scanner-generalization neural network

The general psychopathology factor (p-factor) represents shared variance across mental disorders based on psychopathologic symptoms. The Adolescent Brain Cognitive Development (ABCD) Study offers an unprecedented opportunity to investigate functional networks (FNs) from functional magnetic resonance imaging (fMRI) associated with the psychopathology of an adolescent cohort (n > 10,000). However, the heterogeneities associated with the use of multiple sites and multiple scanners in the ABCD Study need to be overcome to improve the prediction of the p-factor using fMRI. We proposed a scanner-generalization neural network (SGNN) to predict the individual p-factor by systematically reducing the scanner effect for resting-state functional connectivity (RSFC). We included 6905 adolescents from 18 sites whose fMRI data were collected using either Siemens or GE scanners. The p-factor was estimated based on the Child Behavior Checklist (CBCL) scores available in the ABCD study using exploratory factor analysis. We evaluated the Pearson's correlation coefficients (CCs) for p-factor prediction via leave-one/two-site-out cross-validation (LOSOCV/LTSOCV) and identified important FNs from the weight features (WFs) of the SGNN. The CCs were higher for the SGNN than for alternative models when using both LOSOCV (0.1631 ± 0.0673 for the SGNN vs. 0.1497 ± 0.0710 for kernel ridge regression [KRR]; p < 0.05 from a two-tailed paired t-test) and LTSOCV (0.1469 ± 0.0381 for the SGNN vs. 0.1394 ± 0.0359 for KRR; p = 0.01). It was found that (a) the default-mode and dorsal attention FNs were important for p-factor prediction, and (b) the intra-visual FN was important for scanner generalization. We demonstrated the efficacy of our novel SGNN model for p-factor prediction while simultaneously eliminating scanner-related confounding effects for RSFC.

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

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

A Review of Fusion Methods for Omics and Imaging Data

The development of omics data and biomedical images has greatly advanced the progress of precision medicine in diagnosis, treatment, and prognosis. The fusion of omics and imaging data, i.e., omics-imaging fusion, offers a new strategy for understanding complex diseases. However, due to a variety of issues such as the limited number of samples, high dimensionality of features, and heterogeneity of different data types, efficiently learning complementary or associated discriminative fusion information from omics and imaging data remains a challenge. Recently, numerous machine learning methods have been proposed to alleviate these problems. In this review, from the perspective of fusion levels and fusion methods, we first provide an overview of preprocessing and feature extraction methods for omics and imaging data, and comprehensively analyze and summarize the basic forms and variations of commonly used and newly emerging fusion methods, along with their advantages, disadvantages and the applicable scope. We then describe public datasets and compare experimental results of various fusion methods on the ADNI and TCGA datasets. Finally, we discuss future prospects and highlight remaining challenges in the field.

Transdiagnostic symptom networks in relation to mental health service use in community youth

OBJECTIVE

The objective of this study is to scrutinize whether psychopathology symptom networks differ between those with and without lifetime: treatment seeking, treatment and treatment of longer duration.

METHODS

We created non-exclusive groups of subjects with versus without lifetime treatment seeking, treatment and treatment of mid-long-term duration. We estimated Ising models and carried out network comparison tests (NCTs) to compare (a) overall connectivity and (b) network structure. Furthermore, we examined node strength. We used propensity score matching (PSM) to minimize potential confounding by indication for service use.

RESULTS

Based on data from 9,172 participants, there were no statistically significant differences in overall connectivity and network structure in those with versus without lifetime: treatment seeking (p = .75 and p = .82, respectively), treatment (p = .63 and p = .49, respectively) and treatment of mid-longterm duration (p = .15 and p = .62, respectively). Notably, comparing networks with versus without service use consistently revealed higher node strength in 'obsessions' and 'aggression' and lower node strength in 'elevated mood' in all networks with service use.

CONCLUSIONS

Findings suggest that after adjusting for potential confounding by indication for service use, there was no indication of an association in overall connectivity or network structure for lifetime treatment seeking, treatment and treatment of longer duration. However, selected structurally important symptoms differed consistently in all three comparisons. Our findings highlight the potential of network analysis methods to examine treatment mechanisms and outcomes. Specifically, more granular network characteristics on the node level may complement and enrich traditional outcomes in clinical research.

Asymmetric signaling across the hierarchy of cytoarchitecture within the human connectome

Cortical variations in cytoarchitecture form a sensory-fugal axis that shapes regional profiles of extrinsic connectivity and is thought to guide signal propagation and integration across the cortical hierarchy. While neuroimaging work has shown that this axis constrains local properties of the human connectome, it remains unclear whether it also shapes the asymmetric signaling that arises from higher-order topology. Here, we used network control theory to examine the amount of energy required to propagate dynamics across the sensory-fugal axis. Our results revealed an asymmetry in this energy, indicating that bottom-up transitions were easier to complete compared to top-down. Supporting analyses demonstrated that asymmetries were underpinned by a connectome topology that is wired to support efficient bottom-up signaling. Lastly, we found that asymmetries correlated with differences in communicability and intrinsic neuronal time scales and lessened throughout youth. Our results show that cortical variation in cytoarchitecture may guide the formation of macroscopic connectome topology.

Closing the life-cycle of normative modeling using federated hierarchical Bayesian regression

Clinical neuroimaging data availability has grown substantially in the last decade, providing the potential for studying heterogeneity in clinical cohorts on a previously unprecedented scale. Normative modeling is an emerging statistical tool for dissecting heterogeneity in complex brain disorders. However, its application remains technically challenging due to medical data privacy issues and difficulties in dealing with nuisance variation, such as the variability in the image acquisition process. Here, we approach the problem of estimating a reference normative model across a massive population using a massive multi-center neuroimaging dataset. To this end, we introduce a federated probabilistic framework using hierarchical Bayesian regression (HBR) to complete the life-cycle of normative modeling. The proposed model provides the possibilities to learn, update, and adapt the model parameters on decentralized neuroimaging data. Our experimental results confirm the superiority of HBR in deriving more accurate normative ranges on large multi-site neuroimaging datasets compared to the current standard methods. In addition, our approach provides the possibility to recalibrate and reuse the learned model on local datasets and even on datasets with very small sample sizes. The proposed method will facilitate applications of normative modeling as a medical tool for screening the biological deviations in individuals affected by complex illnesses such as mental disorders.

Towards mapping neuro-behavioral heterogeneity of psychedelic neurobiology in humans

Precision psychiatry aims to identify markers of inter-individual variability that allow predicting the right treatment for each patient. However, bridging the gap between molecular-level manipulations and neural systems-level functional alterations remains an unsolved problem in psychiatry. After decades of low success rates in pharmaceutical R&D for psychiatric drugs, multiple studies now point to the potential of psychedelics as a promising fast-acting and long-lasting treatment for some psychiatric symptoms. Yet, given the highly psychoactive nature of these substances, a precision medicine approach is essential to map the neural signals related to clinical efficacy in order to identify patients who can maximally benefit from this treatment. Recent studies have shown that bridging the gap between pharmacology, systems-level neural response in humans and individual experience is possible for psychedelic substances, therefore paving the way for a precision neuropsychiatric therapeutic development. Specifically, it has been shown that the integration of brain-wide PET or transcriptomic data, i.e. receptor distribution for the serotonin 2A receptor, with computational neuroimaging methods can simulate the effect of psychedelics on the human brain. These novel 'computational psychiatry' approaches allow for modeling inter-individual differences in neural as well as subjective effects of psychedelic substances. Collectively, this review provides a deep dive into psychedelic pharmaco-neuroimaging studies with a core focus on how recent computational psychiatry advances in biophysically based circuit modeling can be leveraged to predict individual responses. Finally, we emphasize the importance of human pharmacological neuroimaging for the continued precision therapeutic development of psychedelics.

Leveraging edge-centric networks complements existing network-level inference for functional connectomes

The human connectome is modular with distinct brain regions clustering together to form large-scale communities, or networks. This concept has recently been leveraged in novel inferencing procedures by averaging the edge-level statistics within networks to induce more powerful inferencing at the network level. However, these networks are constructed based on the similarity between pairs of nodes. Emerging work has described novel edge-centric networks, which instead use the similarity between pairs of edges to construct networks. In this work, we use these edge-centric networks in a network-level inferencing procedure and compare this novel method to traditional inferential procedures and the network-level procedure using node-centric networks. We use data from the Human Connectome Project, the Healthy Brain Network, and the Philadelphia Neurodevelopmental Cohort and use a resampling technique with various sample sizes (n=40, 80, 120) to probe the power and specificity of each method. Across datasets and sample sizes, using the edge-centric networks outperforms using node-centric networks for inference as well as edge-level FDR correction and NBS. Additionally, the edge-centric networks were found to be more consistent in clustering effect sizes of similar values as compared to node-centric networks, although node-centric networks often had a lower average within-network effect size variability. Together, these findings suggest that using edge-centric networks for network-level inference can procure relatively powerful results while remaining similarly accurate to the underlying edge-level effects across networks, complementing previous inferential methods.