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.

Study design features that improve effect sizes in cross-sectional and longitudinal brain-wide association studies

Brain-wide association studies (BWAS) are a fundamental tool in discovering brain-behavior associations. Several recent studies showed that thousands of study participants are required to improve the replicability of BWAS because actual effect sizes are much smaller than those reported in smaller studies. Here, we perform analyses and meta-analyses of a robust effect size index (RESI) using 63 longitudinal and cross-sectional magnetic resonance imaging studies from the Lifespan Brain Chart Consortium (77,695 total scans) to demonstrate that optimizing study design is critical for improving standardized effect sizes and replicability in BWAS. A meta-analysis of brain volume associations with age indicates that BWAS with larger covariate variance have larger effect size estimates and that the longitudinal studies we examined have systematically larger standardized effect sizes than cross-sectional studies. We propose a cross-sectional RESI to adjust for the systematic difference in effect sizes between cross-sectional and longitudinal studies that allows investigators to quantify the benefit of conducting their study longitudinally. Analyzing age effects on global and regional brain measures from the United Kingdom Biobank and the Alzheimer's Disease Neuroimaging Initiative, we show that modifying longitudinal study design through sampling schemes to increase between-subject variability and adding a single additional longitudinal measurement per subject can improve effect sizes. However, evaluating these longitudinal sampling schemes on cognitive, psychopathology, and demographic associations with structural and functional brain outcome measures in the Adolescent Brain and Cognitive Development dataset shows that commonly used longitudinal models can, counterintuitively, reduce effect sizes. We demonstrate that the benefit of conducting longitudinal studies depends on the strengths of the between- and within-subject associations of the brain and non-brain measures. Explicitly modeling between- and within-subject effects avoids conflating the effects and allows optimizing effect sizes for them separately. These findings underscore the importance of considering study design features to improve the replicability of BWAS.

A general exposome factor explains individual differences in functional brain network topography and cognition in youth

Childhood environments are critical in shaping cognitive neurodevelopment. With the increasing availability of large-scale neuroimaging datasets with deep phenotyping of childhood environments, we can now build upon prior studies that have considered relationships between one or a handful of environmental and neuroimaging features at a time. Here, we characterize the combined effects of hundreds of inter-connected and co-occurring features of a child's environment ("exposome") and investigate associations with each child's unique, multidimensional pattern of functional brain network organization ("functional topography") and cognition. We apply data-driven computational models to measure the exposome and define personalized functional brain networks in pre-registered analyses. Across matched discovery (n=5139, 48.5% female) and replication (n=5137, 47.1% female) samples from the Adolescent Brain Cognitive Development study, the exposome was associated with current (ages 9-10) and future (ages 11-12) cognition. Changes in the exposome were also associated with changes in cognition after accounting for baseline scores. Cross-validated ridge regressions revealed that the exposome is reflected in functional topography and can predict performance across cognitive domains. Importantly, a single measure capturing a child's exposome could more accurately and parsimoniously predict cognition than a wealth of personalized neuroimaging data, highlighting the importance of children's complex, multidimensional environments in cognitive neurodevelopment.

A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

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.

A precision functional atlas of personalized network topography and probabilities

Although the general location of functional neural networks is similar across individuals, there is vast person-to-person topographic variability. To capture this, we implemented precision brain mapping functional magnetic resonance imaging methods to establish an open-source, method-flexible set of precision functional network atlases-the Masonic Institute for the Developing Brain (MIDB) Precision Brain Atlas. This atlas is an evolving resource comprising 53,273 individual-specific network maps, from more than 9,900 individuals, across ages and cohorts, including the Adolescent Brain Cognitive Development study, the Developmental Human Connectome Project and others. We also generated probabilistic network maps across multiple ages and integration zones (using a new overlapping mapping technique, Overlapping MultiNetwork Imaging). Using regions of high network invariance improved the reproducibility of executive function statistical maps in brain-wide associations compared to group average-based parcellations. Finally, we provide a potential use case for probabilistic maps for targeted neuromodulation. The atlas is expandable to alternative datasets with an online interface encouraging the scientific community to explore and contribute to understanding the human brain function more precisely.

Generalizable Links Between Borderline Personality Traits and Functional Connectivity

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Delineating Transdiagnostic Subtypes in Neurodevelopmental Disorders via Contrastive Graph Machine Learning of Brain Connectivity Patterns

Neurodevelopmental disorders, such as Attention Deficit/Hyperactivity Disorder (ADHD) and Autism Spectrum Disorder (ASD), are characterized by comorbidity and heterogeneity. Identifying distinct subtypes within these disorders can illuminate the underlying neurobiological and clinical characteristics, paving the way for more tailored treatments. We adopted a novel transdiagnostic approach across ADHD and ASD, using cutting-edge contrastive graph machine learning to determine subtypes based on brain network connectivity as revealed by resting-state functional magnetic resonance imaging. Our approach identified two generalizable subtypes characterized by robust and distinct functional connectivity patterns, prominently within the frontoparietal control network and the somatomotor network. These subtypes exhibited pronounced differences in major cognitive and behavioural measures. We further demonstrated the generalizability of these subtypes using data collected from independent study sites. Our data-driven approach provides a novel solution for parsing biological heterogeneity in neurodevelopmental disorders.

Q&A with Ted Satterthwaite and Joëlle Bagautdinova

We at Cell Reports discuss the work, interests, and mentoring experiences of Theodore Satterthwaite (TS) and his graduate student and co-author, Joëlle Bagautdinova (JB). They share with us their recent findings highlighting the relationship between the development of cognitive function and white matter and also talk about the challenges and technical advances in cognitive neuroscience and neuroimaging.

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.

Connectome architecture shapes large-scale cortical alterations in schizophrenia: a worldwide ENIGMA study

Schizophrenia is a prototypical network disorder with widespread brain-morphological alterations, yet it remains unclear whether these distributed alterations robustly reflect the underlying network layout. We tested whether large-scale structural alterations in schizophrenia relate to normative structural and functional connectome architecture, and systematically evaluated robustness and generalizability of these network-level alterations. Leveraging anatomical MRI scans from 2439 adults with schizophrenia and 2867 healthy controls from 26 ENIGMA sites and normative data from the Human Connectome Project (n = 207), we evaluated structural alterations of schizophrenia against two network susceptibility models: (i) hub vulnerability, which examines associations between regional network centrality and magnitude of disease-related alterations; (ii) epicenter mapping, which identifies regions whose typical connectivity profile most closely resembles the disease-related morphological alterations. To assess generalizability and specificity, we contextualized the influence of site, disease stages, and individual clinical factors and compared network associations of schizophrenia with that found in affective disorders. Our findings show schizophrenia-related cortical thinning is spatially associated with functional and structural hubs, suggesting that highly interconnected regions are more vulnerable to morphological alterations. Predominantly temporo-paralimbic and frontal regions emerged as epicenters with connectivity profiles linked to schizophrenia's alteration patterns. Findings were robust across sites, disease stages, and related to individual symptoms. Moreover, transdiagnostic comparisons revealed overlapping epicenters in schizophrenia and bipolar, but not major depressive disorder, suggestive of a pathophysiological continuity within the schizophrenia-bipolar-spectrum. In sum, cortical alterations over the course of schizophrenia robustly follow brain network architecture, emphasizing marked hub susceptibility and temporo-frontal epicenters at both the level of the group and the individual. Subtle variations of epicenters across disease stages suggest interacting pathological processes, while associations with patient-specific symptoms support additional inter-individual variability of hub vulnerability and epicenters in schizophrenia. Our work outlines potential pathways to better understand macroscale structural alterations, and inter- individual variability in schizophrenia.

Transcriptional cartography integrates multiscale biology of the human cortex

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Diffusion MRI head motion correction methods are highly accurate but impacted by denoising and sampling scheme

Head motion correction is particularly challenging in diffusion-weighted MRI (dMRI) scans due to the dramatic changes in image contrast at different gradient strengths and directions. Head motion correction is typically performed using a Gaussian Process model implemented in FSL's Eddy. Recently, the 3dSHORE-based SHORELine method was introduced that does not require shell-based acquisitions, but it has not been previously benchmarked. Here we perform a comprehensive evaluation of both methods on realistic simulations of a software fiber phantom that provides known ground-truth head motion. We demonstrate that both methods perform remarkably well, but that performance can be impacted by sampling scheme and the extent of head motion and the denoising strategy applied before head motion correction. Furthermore, we find Eddy benefits from denoising the data first with MP-PCA. In sum, we provide the most extensive known benchmarking of dMRI head motion correction, together with extensive simulation data and a reproducible workflow. PRACTITIONER POINTS: Both Eddy and SHORELine head motion correction methods performed quite well on a large variety of simulated data. Denoising with MP-PCA can improve head motion correction performance when Eddy is used. SHORELine effectively corrects motion in non-shelled diffusion spectrum imaging data.

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.

Genomic loci influence patterns of structural covariance in the human brain

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

Personalized functional brain network topography is associated with individual differences in youth cognition

Individual differences in cognition during childhood are associated with important social, physical, and mental health outcomes in adolescence and adulthood. Given that cortical surface arealization during development reflects the brain's functional prioritization, quantifying variation in the topography of functional brain networks across the developing cortex may provide insight regarding individual differences in cognition. We test this idea by defining personalized functional networks (PFNs) that account for interindividual heterogeneity in functional brain network topography in 9-10 year olds from the Adolescent Brain Cognitive Development℠ Study. Across matched discovery (n = 3525) and replication (n = 3447) samples, the total cortical representation of fronto-parietal PFNs positively correlates with general cognition. Cross-validated ridge regressions trained on PFN topography predict cognition in unseen data across domains, with prediction accuracy increasing along the cortex's sensorimotor-association organizational axis. These results establish that functional network topography heterogeneity is associated with individual differences in cognition before the critical transition into adolescence.

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.

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.

Normative Modeling of Brain Morphometry Across the Lifespan Using CentileBrain: Algorithm Benchmarking and Model Optimization

We present an empirically benchmarked framework for sex-specific normative modeling of brain morphometry that can inform about the biological and behavioral significance of deviations from typical age-related neuroanatomical changes and support future study designs. This framework was developed using regional morphometric data from 37,407 healthy individuals (53% female; aged 3-90 years) following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The Multivariate Factorial Polynomial Regression (MFPR) emerged as the preferred algorithm optimized using nonlinear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins, and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3,000 study participants. The model and scripts described here are freely available through CentileBrain (https://centilebrain.org/).

Framewise multi-echo distortion correction for superior functional MRI

Functional MRI (fMRI) data are severely distorted by magnetic field (B0) inhomogeneities which currently must be corrected using separately acquired field map data. However, changes in the head position of a scanning participant across fMRI frames can cause changes in the B0 field, preventing accurate correction of geometric distortions. Additionally, field maps can be corrupted by movement during their acquisition, preventing distortion correction altogether. In this study, we use phase information from multi-echo (ME) fMRI data to dynamically sample distortion due to fluctuating B0 field inhomogeneity across frames by acquiring multiple echoes during a single EPI readout. Our distortion correction approach, MEDIC (Multi-Echo DIstortion Correction), accurately estimates B0 related distortions for each frame of multi-echo fMRI data. Here, we demonstrate that MEDIC's framewise distortion correction produces improved alignment to anatomy and decreases the impact of head motion on resting-state functional connectivity (RSFC) maps, in higher motion data, when compared to the prior gold standard approach (i.e., TOPUP). Enhanced framewise distortion correction with MEDIC, without the requirement for field map collection, furthers the advantage of multi-echo over single-echo fMRI.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

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.

Mapping the Relationship of White Matter Lesions to Depression in Multiple Sclerosis

BACKGROUND

Multiple sclerosis (MS) is an immune-mediated neurological disorder, and up to 50% of patients experience depression. We investigated how white matter network disruption is related to depression in MS.

METHODS

Using electronic health records, 380 participants with MS were identified. Depressed individuals (MS+Depression group; n = 232) included persons who had an ICD-10 depression diagnosis, had a prescription for antidepressant medication, or screened positive via Patient Health Questionnaire (PHQ)-2 or PHQ-9. Age- and sex-matched nondepressed individuals with MS (MS-Depression group; n = 148) included persons who had no prior depression diagnosis, had no psychiatric medication prescriptions, and were asymptomatic on PHQ-2 or PHQ-9. Research-quality 3T structural magnetic resonance imaging was obtained as part of routine care. We first evaluated whether lesions were preferentially located within the depression network compared with other brain regions. Next, we examined if MS+Depression patients had greater lesion burden and if this was driven by lesions in the depression network. Primary outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain.

RESULTS

MS lesions preferentially affected fascicles within versus outside the depression network (β = 0.09, 95% CI = 0.08 to 0.10, p < .001). MS+Depression patients had more lesion burden (β = 0.06, 95% CI = 0.01 to 0.10, p = .015); this was driven by lesions within the depression network (β = 0.02, 95% CI = 0.003 to 0.040, p = .020).

CONCLUSIONS

We demonstrated that lesion location and burden may contribute to depression comorbidity in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression patients had more disease than MS-Depression patients, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

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 the spatial structure 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 genomics 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 phenotype associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

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

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

ALE meta-analyses of voxel-based morphometry studies: Parameter validation via large-scale simulations

Activation likelihood estimation (ALE) meta-analysis has been applied to structural neuroimaging data since long, but up to now, any systematic assessment of the algorithm's behavior, power and sensitivity has been based on simulations using functional neuroimaging databases as their foundation. Here, we aimed to determine whether the guidelines offered by previous evaluations can be generalized to ALE meta-analyses of voxel-based morphometry (VBM) studies. We ran 365000 distinct ALE analyses filled with simulated experiments, randomly sampling parameters from BrainMap's VBM experiment database. We then examined the algorithm's sensitivity, its susceptibility to spurious convergence, and its susceptibility to excessive contributions by individual experiments. In general, the performance of the ALE algorithm was highly comparable between imaging modalities, with the algorithm's sensitivity and specificity reaching similar levels with structural data as previously observed with functional data. Because of the lower number of foci reported and the higher number of participants usually included in structural experiments, individual studies had, on average, a higher impact towards significant clusters. To prevent significant clusters from being driven by single experiments, we recommend that researchers include at least 23 experiments in a VBM ALE dataset, instead of the previously recommended minimum of n = 17. While these recommendations do not constitute hard borders, running ALE analyses on smaller datasets would require special diligence in assessing and reporting the contributions of experiments to individual clusters.

Development and application of novel performance validity metrics for computerized neurocognitive batteries

OBJECTIVES

Data from neurocognitive assessments may not be accurate in the context of factors impacting validity, such as disengagement, unmotivated responding, or intentional underperformance. Performance validity tests (PVTs) were developed to address these phenomena and assess underperformance on neurocognitive tests. However, PVTs can be burdensome, rely on cutoff scores that reduce information, do not examine potential variations in task engagement across a battery, and are typically not well-suited to acquisition of large cognitive datasets. Here we describe the development of novel performance validity measures that could address some of these limitations by leveraging psychometric concepts using data embedded within the Penn Computerized Neurocognitive Battery (PennCNB).

METHODS

We first developed these validity measures using simulations of invalid response patterns with parameters drawn from real data. Next, we examined their application in two large, independent samples: 1) children and adolescents from the Philadelphia Neurodevelopmental Cohort (n = 9498); and 2) adult servicemembers from the Marine Resiliency Study-II (n = 1444).

RESULTS

Our performance validity metrics detected patterns of invalid responding in simulated data, even at subtle levels. Furthermore, a combination of these metrics significantly predicted previously established validity rules for these tests in both developmental and adult datasets. Moreover, most clinical diagnostic groups did not show reduced validity estimates.

CONCLUSIONS

These results provide proof-of-concept evidence for multivariate, data-driven performance validity metrics. These metrics offer a novel method for determining the performance validity for individual neurocognitive tests that is scalable, applicable across different tests, less burdensome, and dimensional. However, more research is needed into their application.

A critical period plasticity framework for the sensorimotor-association axis of cortical neurodevelopment

To understand human brain development it is necessary to describe not only the spatiotemporal patterns of neurodevelopment but also the neurobiological mechanisms that underlie them. Human neuroimaging studies have provided evidence for a hierarchical sensorimotor-to-association (S-A) axis of cortical neurodevelopment. Understanding the biological mechanisms that underlie this program of development using traditional neuroimaging approaches has been challenging. Animal models have been used to identify periods of enhanced experience-dependent plasticity - 'critical periods' - that progress along cortical hierarchies and are governed by a conserved set of neurobiological mechanisms that promote and then restrict plasticity. In this review we hypothesize that the S-A axis of cortical development in humans is partly driven by the cascading maturation of critical period plasticity mechanisms. We then describe how recent advances in in vivo neuroimaging approaches provide a promising path toward testing this hypothesis by linking signals derived from non-invasive imaging to critical period mechanisms.

Brain Growth Charts for Quantitative Analysis of Pediatric Clinical Brain MRI Scans with Limited Imaging Pathology

Background Clinically acquired brain MRI scans represent a valuable but underused resource for investigating neurodevelopment due to their technical heterogeneity and lack of appropriate controls. These barriers have curtailed retrospective studies of clinical brain MRI scans compared with more costly prospectively acquired research-quality brain MRI scans. Purpose To provide a benchmark for neuroanatomic variability in clinically acquired brain MRI scans with limited imaging pathology (SLIPs) and to evaluate if growth charts from curated clinical MRI scans differed from research-quality MRI scans or were influenced by clinical indication for the scan. Materials and Methods In this secondary analysis of preexisting data, clinical brain MRI SLIPs from an urban pediatric health care system (individuals aged ≤22 years) were scanned across nine 3.0-T MRI scanners. The curation process included manual review of signed radiology reports and automated and manual quality review of images without gross pathology. Global and regional volumetric imaging phenotypes were measured using two image segmentation pipelines, and clinical brain growth charts were quantitatively compared with charts derived from a large set of research controls in the same age range by means of Pearson correlation and age at peak volume. Results The curated clinical data set included 532 patients (277 male; median age, 10 years [IQR, 5-14 years]; age range, 28 days after birth to 22 years) scanned between 2005 and 2020. Clinical brain growth charts were highly correlated with growth charts derived from research data sets (22 studies, 8346 individuals [4947 male]; age range, 152 days after birth to 22 years) in terms of normative developmental trajectories predicted by the models (median r = 0.979). Conclusion The clinical indication of the scans did not significantly bias the output of clinical brain charts. Brain growth charts derived from clinical controls with limited imaging pathology were highly correlated with brain charts from research controls, suggesting the potential of curated clinical MRI scans to supplement research data sets. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Ertl-Wagner and Pai in this issue.

Myelination and excitation-inhibition balance synergistically shape structure-function coupling across the human cortex

Recent work has demonstrated that the relationship between structural and functional connectivity varies regionally across the human brain, with reduced coupling emerging along the sensory-association cortical hierarchy. The biological underpinnings driving this expression, however, remain largely unknown. Here, we postulate that intracortical myelination and excitation-inhibition (EI) balance mediate the heterogeneous expression of structure-function coupling (SFC) and its temporal variance across the cortical hierarchy. We employ atlas- and voxel-based connectivity approaches to analyze neuroimaging data acquired from two groups of healthy participants. Our findings are consistent across six complementary processing pipelines: 1) SFC and its temporal variance respectively decrease and increase across the unimodal-transmodal and granular-agranular gradients; 2) increased myelination and lower EI-ratio are associated with more rigid SFC and restricted moment-to-moment SFC fluctuations; 3) a gradual shift from EI-ratio to myelination as the principal predictor of SFC occurs when traversing from granular to agranular cortical regions. Collectively, our work delivers a framework to conceptualize structure-function relationships in the human brain, paving the way for an improved understanding of how demyelination and/or EI-imbalances induce reorganization in brain disorders.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

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.

Deviations from a typical development of the cerebellum in youth are associated with psychopathology, executive functions and educational outcomes

BACKGROUND

Understanding deviations from typical brain development is a promising approach to comprehend pathophysiology in childhood and adolescence. We investigated if cerebellar volumes different than expected for age and sex could predict psychopathology, executive functions and academic achievement.

METHODS

Children and adolescents aged 6-17 years from the Brazilian High-Risk Cohort Study for Mental Conditions had their cerebellar volume estimated using Multiple Automatically Generated Templates from T1-weighted images at baseline (n = 677) and at 3-year follow-up (n = 447). Outcomes were assessed using the Child Behavior Checklist and standardized measures of executive functions and school achievement. Models of typically developing cerebellum were based on a subsample not exposed to risk factors and without mental-health conditions (n = 216). Deviations from this model were constructed for the remaining individuals (n = 461) and standardized variation from age and sex trajectory model was used to predict outcomes in cross-sectional, longitudinal and mediation analyses.

RESULTS

Cerebellar volumes higher than expected for age and sex were associated with lower externalizing specific factor and higher executive functions. In a longitudinal analysis, deviations from typical development at baseline predicted inhibitory control at follow-up, and cerebellar deviation changes from baseline to follow-up predicted changes in reading and writing abilities. The association between deviations in cerebellar volume and academic achievement was mediated by inhibitory control.

CONCLUSIONS

Deviations in the cerebellar typical development are associated with outcomes in youth that have long-lasting consequences. This study highlights both the potential of typical developing models and the important role of the cerebellum in mental health, cognition and education.

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

Generalizable links between symptoms of borderline personality disorder and functional connectivity

Background |

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

Methods |

We used high-quality functional Magnetic Resonance Imaging (fMRI) data from young adults from the Human Connectome Project: Young Adults (HCP-YA; N = 870, ages 22-37 years, 457 female) and youth from the Human Connectome Project: Development (HCP-D; N = 223, age range 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five Factor Inventory (NEO-FFI). A ridge regression model with 10-fold cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity, while controlling for in-scanner motion, age, and sex. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity.

Results |

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

Conclusion |

Individual differences in functional connectivity in developmentally-sensitive regions are associated with the symptoms of BPD.

A reproducible and generalizable software workflow for analysis of large-scale neuroimaging data collections using BIDS Apps

Neuroimaging research faces a crisis of reproducibility. With massive sample sizes and greater data complexity, this problem becomes more acute. Software that operates on imaging data defined using the Brain Imaging Data Structure (BIDS) - BIDS Apps - have provided a substantial advance. However, even using BIDS Apps, a full audit trail of data processing is a necessary prerequisite for fully reproducible research. Obtaining a faithful record of the audit trail is challenging - especially for large datasets. Recently, the FAIRly big framework was introduced as a way to facilitate reproducible processing of large-scale data by leveraging DataLad - a version control system for data management. However, the current implementation of this framework was more of a proof of concept, and could not be immediately reused by other investigators for different use cases. Here we introduce the BIDS App Bootstrap (BABS), a user-friendly and generalizable Python package for reproducible image processing at scale. BABS facilitates the reproducible application of BIDS Apps to large-scale datasets. Leveraging DataLad and the FAIRly big framework, BABS tracks the full audit trail of data processing in a scalable way by automatically preparing all scripts necessary for data processing and version tracking on high performance computing (HPC) systems. Currently, BABS supports jobs submissions and audits on Sun Grid Engine (SGE) and Slurm HPCs with a parsimonious set of programs. To demonstrate its scalability, we applied BABS to data from the Healthy Brain Network (HBN; n=2,565). Taken together, BABS allows reproducible and scalable image processing and is broadly extensible via an open-source development model.

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

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

Alprazolam modulates persistence energy during emotion processing in first-degree relatives of individuals with schizophrenia: a network control study

Schizophrenia is marked by deficits in facial affect processing associated with abnormalities in GABAergic circuitry, deficits also found in first-degree relatives. Facial affect processing involves a distributed network of brain regions including limbic regions like amygdala and visual processing areas like fusiform cortex. Pharmacological modulation of GABAergic circuitry using benzodiazepines like alprazolam can be useful for studying this facial affect processing network and associated GABAergic abnormalities in schizophrenia. Here, we use pharmacological modulation and computational modeling to study the contribution of GABAergic abnormalities toward emotion processing deficits in schizophrenia. Specifically, we apply principles from network control theory to model persistence energy - the control energy required to maintain brain activation states - during emotion identification and recall tasks, with and without administration of alprazolam, in a sample of first-degree relatives and healthy controls. Here, persistence energy quantifies the magnitude of theoretical external inputs during the task. We find that alprazolam increases persistence energy in relatives but not in controls during threatening face processing, suggesting a compensatory mechanism given the relative absence of behavioral abnormalities in this sample of unaffected relatives. Further, we demonstrate that regions in the fusiform and occipital cortices are important for facilitating state transitions during facial affect processing. Finally, we uncover spatial relationships (i) between regional variation in differential control energy (alprazolam versus placebo) and (ii) both serotonin and dopamine neurotransmitter systems, indicating that alprazolam may exert its effects by altering neuromodulatory systems. Together, these findings provide a new perspective on the distributed emotion processing network and the effect of GABAergic modulation on this network, in addition to identifying an association between schizophrenia risk and abnormal GABAergic effects on persistence energy during threat processing.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

Development of Iron Status Measures during Youth: Associations with Sex, Neighborhood Socioeconomic Status, Cognitive Performance, and Brain Structure

BACKGROUND

Iron is essential to brain function, and iron deficiency during youth may adversely impact neurodevelopment. Understanding the developmental time course of iron status and its association with neurocognitive functioning is important for identifying windows for intervention.

OBJECTIVES

This study aimed to characterize developmental change in iron status and understand its association with cognitive performance and brain structure during adolescence using data from a large pediatric health network.

METHODS

This study included a cross-sectional sample of 4899 participants (2178 males; aged 8-22 y at the time of participation, M [SD] = 14.24 [3.7]) who were recruited from the Children's Hospital of Philadelphia network. Prospectively collected research data were enriched with electronic medical record data that included hematological measures related to iron status, including serum hemoglobin, ferritin, and transferrin (33,015 total samples). At the time of participation, cognitive performance was assessed using the Penn Computerized Neurocognitive Battery, and brain white matter integrity was assessed using diffusion-weighted MRI in a subset of individuals.

RESULTS

Developmental trajectories were characterized for all metrics and revealed that sex differences emerged after menarche such that females had reduced iron status relative to males [all R2partial > 0.008; all false discovery rates (FDRs) < 0.05]. Higher socioeconomic status was associated with higher hemoglobin concentrations throughout development (R2partial = 0.005; FDR < 0.001), and the association was greatest during adolescence. Higher hemoglobin concentrations were associated with better cognitive performance during adolescence (R2partial = 0.02; FDR < 0.001) and mediated the association between sex and cognition (mediation effect = -0.107; 95% CI: -0.191, -0.02). Higher hemoglobin concentration was also associated with greater brain white matter integrity in the neuroimaging subsample (R2partial = 0.06, FDR = 0.028).

CONCLUSIONS

Iron status evolves during youth and is lowest in females and individuals of low socioeconomic status during adolescence. Diminished iron status during adolescence has consequences for neurocognition, suggesting that this critical period of neurodevelopment may be an important window for intervention that has the potential to reduce health disparities in at-risk populations.

Mapping the relationship of white matter lesions to depression in multiple sclerosis

Importance

Multiple sclerosis (MS) is an immune-mediated neurological disorder that affects nearly one million people in the United States. Up to 50% of patients with MS experience depression.

Objective

To investigate how white matter network disruption is related to depression in MS.

Design

Retrospective case-control study of participants who received research-quality 3-tesla neuroimaging as part of MS clinical care from 2010-2018. Analyses were performed from May 1 to September 30, 2022.

Setting

Single-center academic medical specialty MS clinic.

Participants

Participants with MS were identified via the electronic health record (EHR). All participants were diagnosed by an MS specialist and completed research-quality MRI at 3T. After excluding participants with poor image quality, 783 were included. Inclusion in the depression group (MS+Depression) required either: 1) ICD-10 depression diagnosis (F32-F34.*); 2) prescription of antidepressant medication; or 3) screening positive via Patient Health Questionnaire-2 (PHQ-2) or -9 (PHQ-9). Age- and sex-matched nondepressed comparators (MS-Depression) included persons with no depression diagnosis, no psychiatric medications, and were asymptomatic on PHQ-2/9.

Exposure

Depression diagnosis.

Main Outcomes and Measures

We first evaluated if lesions were preferentially located within the depression network compared to other brain regions. Next, we examined if MS+Depression patients had greater lesion burden, and if this was driven by lesions specifically in the depression network. Outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain. Secondary measures included between-diagnosis lesion burden, stratified by brain network. Linear mixed-effects models were employed.

Results

Three hundred-eighty participants met inclusion criteria, (232 MS+Depression: age[SD]=49[12], %females=86; 148 MS-Depression: age[SD]=47[13], %females=79). MS lesions preferentially affected fascicles within versus outside the depression network (β=0.09, 95% CI=0.08-0.10, P<0.001). MS+Depression had more white matter lesion burden (β=0.06, 95% CI=0.01-0.10, P=0.015); this was driven by lesions within the depression network (β=0.02, 95% CI 0.003-0.040, P=0.020).

Conclusions and Relevance

We provide new evidence supporting a relationship between white matter lesions and depression in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression had more disease than MS-Depression, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

Dissecting Symptom-linked Dimensions of Resting-State Electroencephalographic Functional Connectivity in Autism with Contrastive Learning

Autism spectrum disorder (ASD) is a common neurodevelopmental disorder characterized by social interaction deficits, communication difficulties, and restricted/repetitive behaviors or fixated interests. Despite its high prevalence, development of effective therapy for ASD is hindered by its symptomatic and neurophysiological heterogeneities. To collectively dissect the ASD heterogeneity in neurophysiology and symptoms, we develop a new analytical framework combining contrastive learning and sparse canonical correlation analysis to identify resting-state EEG connectivity dimensions linked to ASD behavioral symptoms within 392 ASD samples. Two dimensions are successfully identified, showing significant correlations with social/communication deficits (r = 0.70) and restricted/repetitive behaviors (r = 0.45), respectively. We confirm the robustness of these dimensions through cross-validation and further demonstrate their generalizability using an independent dataset of 223 ASD samples. Our results reveal that the right inferior parietal lobe is the core region displaying EEG activity associated with restricted/repetitive behaviors, and functional connectivity between the left angular gyrus and the right middle temporal gyrus is a promising biomarker of social/communication deficits. Overall, these findings provide a promising avenue to parse ASD heterogeneity with high clinical translatability, paving the way for treatment development and precision medicine for ASD.

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

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

Activation of Internal Correctness Monitoring Circuitry in Youths With Psychosis Spectrum Symptoms

BACKGROUND

Self-directed performance monitoring is a critical contributor to cognitive performance and general functioning and is impacted by psychiatric symptoms and personality traits, but has been understudied in psychosis-risk states. We have shown that ventral striatum (VS) responds to correctness during cognitive tasks where no explicit feedback is required, and this intrinsic reinforcement response is reduced in schizophrenia.

METHODS

Here, we examined this phenomenon in youths (n = 796, age range 11-22 years) from the Philadelphia Neurodevelopmental Cohort (PNC) performing a working memory functional magnetic resonance imaging task. We hypothesized that VS would respond to internal correctness monitoring, while classic salience network regions, such as dorsal anterior cingulate cortex and anterior insular cortex, would reflect internal error monitoring and that these responses would increase with age. We expected that neurobehavioral measures of performance monitoring would be reduced in youths with subclinical psychosis spectrum features and would correlate with amotivation severity.

RESULTS

Supporting these hypotheses, we found correct>incorrect activation in VS and incorrect>correct activation in anterior cingulate cortex and anterior insular cortex. Furthermore, VS activation was positively correlated with age, reduced in youths with psychosis spectrum features, and inversely correlated with amotivation. However, these patterns were not significant in anterior cingulate cortex and anterior insular cortex.

CONCLUSIONS

These findings advance our understanding of the neural underpinnings of performance monitoring and its impairment in adolescents with psychosis spectrum features. Such understanding can facilitate investigation of the developmental trajectory of normative and aberrant performance monitoring; contribute to early identification of youths at elevated risk for poor academic, occupational, or psychiatric outcomes; and provide potential targets for therapeutic development.

ModelArray: An R package for statistical analysis of fixel-wise data

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

Development of top-down cortical propagations in youth

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

Assessment of Neuroanatomical Endophenotypes of Autism Spectrum Disorder and Association With Characteristics of Individuals With Schizophrenia and the General Population

Importance

Autism spectrum disorder (ASD) is associated with significant clinical, neuroanatomical, and genetic heterogeneity that limits precision diagnostics and treatment.

Objective

To assess distinct neuroanatomical dimensions of ASD using novel semisupervised machine learning methods and to test whether the dimensions can serve as endophenotypes also in non-ASD populations.

Design, Setting, and Participants

This cross-sectional study used imaging data from the publicly available Autism Brain Imaging Data Exchange (ABIDE) repositories as the discovery cohort. The ABIDE sample included individuals diagnosed with ASD aged between 16 and 64 years and age- and sex-match typically developing individuals. Validation cohorts included individuals with schizophrenia from the Psychosis Heterogeneity Evaluated via Dimensional Neuroimaging (PHENOM) consortium and individuals from the UK Biobank to represent the general population. The multisite discovery cohort included 16 internationally distributed imaging sites. Analyses were performed between March 2021 and March 2022.

Main Outcomes and Measures

The trained semisupervised heterogeneity through discriminative analysis models were tested for reproducibility using extensive cross-validations. It was then applied to individuals from the PHENOM and the UK Biobank. It was hypothesized that neuroanatomical dimensions of ASD would display distinct clinical and genetic profiles and would be prominent also in non-ASD populations.

Results

Heterogeneity through discriminative analysis models trained on T1-weighted brain magnetic resonance images of 307 individuals with ASD (mean [SD] age, 25.4 [9.8] years; 273 [88.9%] male) and 362 typically developing control individuals (mean [SD] age, 25.8 [8.9] years; 309 [85.4%] male) revealed that a 3-dimensional scheme was optimal to capture the ASD neuroanatomy. The first dimension (A1: aginglike) was associated with smaller brain volume, lower cognitive function, and aging-related genetic variants (FOXO3; Z = 4.65; P = 1.62 × 10-6). The second dimension (A2: schizophrenialike) was characterized by enlarged subcortical volumes, antipsychotic medication use (Cohen d = 0.65; false discovery rate-adjusted P = .048), partially overlapping genetic, neuroanatomical characteristics to schizophrenia (n = 307), and significant genetic heritability estimates in the general population (n = 14 786; mean [SD] h2, 0.71 [0.04]; P < 1 × 10-4). The third dimension (A3: typical ASD) was distinguished by enlarged cortical volumes, high nonverbal cognitive performance, and biological pathways implicating brain development and abnormal apoptosis (mean [SD] β, 0.83 [0.02]; P = 4.22 × 10-6).

Conclusions and Relevance

This cross-sectional study discovered 3-dimensional endophenotypic representation that may elucidate the heterogeneous neurobiological underpinnings of ASD to support precision diagnostics. The significant correspondence between A2 and schizophrenia indicates a possibility of identifying common biological mechanisms across the 2 mental health diagnoses.