ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion.

Arterial spin labeled (ASL) magnetic resonance imaging (MRI) is the primary method for noninvasively measuring regional brain perfusion in humans. We introduce ASLPrep, a suite of software pipelines that ensure the reproducible and generalizable processing of ASL MRI data.

Sex differences in the functional topography of association networks in youth.

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

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.

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

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

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.

Associations between neighborhood socioeconomic status, parental education, and executive system activation in youth.

Socioeconomic status (SES) can impact cognitive performance, including working memory (WM). As executive systems that support WM undergo functional neurodevelopment during adolescence, environmental stressors at both individual and community levels may influence cognitive outcomes. Here, we sought to examine how SES at the neighborhood and family level impacts task-related activation of the executive system during adolescence and determine whether this effect mediates the relationship between SES and WM performance. To address these questions, we studied 1,150 youths (age 8-23) that completed a fractal n-back WM task during functional magnetic resonance imaging at 3T as part of the Philadelphia Neurodevelopmental Cohort. We found that both higher neighborhood SES and parental education were associated with greater activation of the executive system to WM load, including the bilateral dorsolateral prefrontal cortex, posterior parietal cortex, and precuneus. The association of neighborhood SES remained significant when controlling for task performance, or related factors like exposure to traumatic events. Furthermore, high-dimensional multivariate mediation analysis identified distinct patterns of brain activity within the executive system that significantly mediated the relationship between measures of SES and task performance. These findings underscore the importance of multilevel environmental factors in shaping executive system function and WM in youth.

Spatially-enhanced clusterwise inference for testing and localizing intermodal correspondence.

With the increasing availability of neuroimaging data from multiple modalities-each providing a different lens through which to study brain structure or function-new techniques for comparing, integrating, and interpreting information within and across modalities have emerged. Recent developments include hypothesis tests of associations between neuroimaging modalities, which can be used to determine the statistical significance of intermodal associations either throughout the entire brain or within anatomical subregions or functional networks. While these methods provide a crucial foundation for inference on intermodal relationships, they cannot be used to answer questions about where in the brain these associations are most pronounced. In this paper, we introduce a new method, called CLEAN-R, that can be used both to test intermodal correspondence throughout the brain and also to localize this correspondence. Our method involves first adjusting for the underlying spatial autocorrelation structure within each modality before aggregating information within small clusters to construct a map of enhanced test statistics. Using structural and functional magnetic resonance imaging data from a subsample of children and adolescents from the Philadelphia Neurodevelopmental Cohort, we conduct simulations and data analyses where we illustrate the high statistical power and nominal type I error levels of our method. By constructing an interpretable map of group-level correspondence using spatially-enhanced test statistics, our method offers insights beyond those provided by earlier methods.

Curation of BIDS (CuBIDS): A workflow and software package for streamlining reproducible curation of large BIDS datasets.

The Brain Imaging Data Structure (BIDS) is a specification accompanied by a software ecosystem that was designed to create reproducible and automated workflows for processing neuroimaging data. BIDS Apps flexibly build workflows based on the metadata detected in a dataset. However, even BIDS valid metadata can include incorrect values or omissions that result in inconsistent processing across sessions. Additionally, in large-scale, heterogeneous neuroimaging datasets, hidden variability in metadata is difficult to detect and classify. To address these challenges, we created a Python-based software package titled "Curation of BIDS" (CuBIDS), which provides an intuitive workflow that helps users validate and manage the curation of their neuroimaging datasets. CuBIDS includes a robust implementation of BIDS validation that scales to large samples and incorporates DataLad--a version control software package for data--as an optional dependency to ensure reproducibility and provenance tracking throughout the entire curation process. CuBIDS provides tools to help users perform quality control on their images' metadata and identify unique combinations of imaging parameters. Users can then execute BIDS Apps on a subset of participants that represent the full range of acquisition parameters that are present, accelerating pipeline testing on large datasets.

Voxel-wise intermodal coupling analysis of two or more modalities using local covariance decomposition.

When individual subjects are imaged with multiple modalities, biological information is present not only within each modality, but also between modalities - that is, in how modalities covary at the voxel level. Previous studies have shown that local covariance structures between modalities, or intermodal coupling (IMCo), can be summarized for two modalities, and that two-modality IMCo reveals otherwise undiscovered patterns in neurodevelopment and certain diseases. However, previous IMCo methods are based on the slopes of local weighted linear regression lines, which are inherently asymmetric and limited to the two-modality setting. Here, we present a generalization of IMCo estimation which uses local covariance decompositions to define a symmetric, voxel-wise coupling coefficient that is valid for two or more modalities. We use this method to study coupling between cerebral blood flow, amplitude of low frequency fluctuations, and local connectivity in 803 subjects ages 8 through 22. We demonstrate that coupling is spatially heterogeneous, varies with respect to age and sex in neurodevelopment, and reveals patterns that are not present in individual modalities. As availability of multi-modal data continues to increase, principal-component-based IMCo (pIMCo) offers a powerful approach for summarizing relationships between multiple aspects of brain structure and function. An R package for estimating pIMCo is available at: https://github.com/hufengling/pIMCo.

Longitudinal Development of Brain Iron Is Linked to Cognition in Youth.

Brain iron is vital to multiple aspects of brain function, including oxidative metabolism, myelination, and neurotransmitter synthesis. Atypical iron concentration in the basal ganglia is associated with neurodegenerative disorders in aging and cognitive deficits. However, the normative development of brain iron concentration in adolescence and its relationship to cognition are less well understood. Here, we address this gap in a longitudinal sample of 922 humans aged 8-26 years at the first visit (M = 15.1, SD = 3.72; 336 males, 486 females) with up to four multiecho T2* scans each. Using this sample of 1236 imaging sessions, we assessed the longitudinal developmental trajectories of tissue iron in the basal ganglia. We quantified tissue iron concentration using R2* relaxometry within four basal ganglia regions, including the caudate, putamen, nucleus accumbens, and globus pallidus. The longitudinal development of R2* was modeled using generalized additive mixed models (GAMMs) with splines to capture linear and nonlinear developmental processes. We observed significant increases in R2* across all regions, with the greatest and most prolonged increases occurring in the globus pallidus and putamen. Further, we found that the developmental trajectory of R2* in the putamen is significantly related to individual differences in cognitive ability, such that greater cognitive ability is increasingly associated with greater iron concentration through late adolescence and young-adulthood. Together, our results suggest a prolonged period of basal ganglia iron enrichment that extends into the mid-twenties, with diminished iron concentration associated with poorer cognitive ability during late adolescence.SIGNIFICANCE STATEMENT Brain tissue iron is essential to healthy brain function. Atypical basal ganglia tissue iron levels have been linked to impaired cognition in iron deficient children and adults with neurodegenerative disorders. However, the normative developmental trajectory of basal ganglia iron concentration during adolescence and its association with cognition are less well understood. In the largest study of tissue iron development yet reported, we characterize the developmental trajectory of tissue iron concentration across the basal ganglia during adolescence and provide evidence that diminished iron content is associated with poorer cognitive performance even in healthy youth. These results highlight the transition from adolescence to adulthood as a period of dynamic maturation of tissue iron concentration in the basal ganglia.

A simple permutation-based test of intermodal correspondence.

Many key findings in neuroimaging studies involve similarities between brain maps, but statistical methods used to measure these findings have varied. Current state-of-the-art methods involve comparing observed group-level brain maps (after averaging intensities at each image location across multiple subjects) against spatial null models of these group-level maps. However, these methods typically make strong and potentially unrealistic statistical assumptions, such as covariance stationarity. To address these issues, in this article we propose using subject-level data and a classical permutation testing framework to test and assess similarities between brain maps. Our method is comparable to traditional permutation tests in that it involves randomly permuting subjects to generate a null distribution of intermodal correspondence statistics, which we compare to an observed statistic to estimate a p-value. We apply and compare our method in simulated and real neuroimaging data from the Philadelphia Neurodevelopmental Cohort. We show that our method performs well for detecting relationships between modalities known to be strongly related (cortical thickness and sulcal depth), and it is conservative when an association would not be expected (cortical thickness and activation on the n-back working memory task). Notably, our method is the most flexible and reliable for localizing intermodal relationships within subregions of the brain and allows for generalizable statistical inference.

Sex differences in network controllability as a predictor of executive function in youth.

Executive function is a quintessential human capacity that emerges late in development and displays different developmental trends in males and females. Sex differences in executive function in youth have been linked to vulnerability to psychopathology as well as to behaviors that impinge on health, wellbeing, and longevity. Yet, the neurobiological basis of these differences is not well understood, in part due to the spatiotemporal complexity inherent in patterns of brain network maturation supporting executive function. Here we test the hypothesis that sex differences in impulsivity in youth stem from sex differences in the controllability of structural brain networks as they rewire over development. Combining methods from network neuroscience and network control theory, we characterize the network control properties of structural brain networks estimated from diffusion imaging data acquired in males and females in a sample of 879 youth aged 8-22 years. We summarize the control properties of these networks by estimating average and modal controllability, two statistics that probe the ease with which brain areas can drive the network towards easy versus difficult-to-reach states. We find that females have higher modal controllability in frontal, parietal, and subcortical regions while males have higher average controllability in frontal and subcortical regions. Furthermore, controllability profiles in males are negatively related to the false positive rate on a continuous performance task, a common measure of impulsivity. Finally, we find associations between average controllability and individual differences in activation during an n-back working memory task. Taken together, our findings support the notion that sex differences in the controllability of structural brain networks can partially explain sex differences in executive function. Controllability of structural brain networks also predicts features of task-relevant activation, suggesting the potential for controllability to represent context-specific constraints on network state more generally.

Robust spatial extent inference with a semiparametric bootstrap joint inference procedure.

Spatial extent inference (SEI) is widely used across neuroimaging modalities to adjust for multiple comparisons when studying brain-phenotype associations that inform our understanding of disease. Recent studies have shown that Gaussian random field (GRF)-based tools can have inflated family-wise error rates (FWERs). This has led to substantial controversy as to which processing choices are necessary to control the FWER using GRF-based SEI. The failure of GRF-based methods is due to unrealistic assumptions about the spatial covariance function of the imaging data. A permutation procedure is the most robust SEI tool because it estimates the spatial covariance function from the imaging data. However, the permutation procedure can fail because its assumption of exchangeability is violated in many imaging modalities. Here, we propose the (semi-) parametric bootstrap joint (PBJ; sPBJ) testing procedures that are designed for SEI of multilevel imaging data. The sPBJ procedure uses a robust estimate of the spatial covariance function, which yields consistent estimates of standard errors, even if the covariance model is misspecified. We use the methods to study the association between performance and executive functioning in a working memory functional magnetic resonance imaging study. The sPBJ has similar or greater power to the PBJ and permutation procedures while maintaining the nominal type 1 error rate in reasonable sample sizes. We provide an R package to perform inference using the PBJ and sPBJ procedures.

Preterm Modulation of Connectivity by Endogenous Generators: The Theta Temporal Activities in Coalescence with Slow Waves.

The neuronal activity of the preterm brain is characterized by various endogenous activities whose roles in neurodevelopmental maturation processes have not been fully elucidated. The preterm EEG is characterized by discontinuities composed of short bursts of activity with dominant low frequencies. One of the earliest endogenous activities is the theta temporal activity in coalescence with slow waves (TTA-SW), which appears at 24 to 32 weeks of gestational age (wGA). The present study investigated the influence of TTA-SW on the spatial organization of the early preterm brain network. To achieve this objective, High-Density EEG data were recorded from preterm infants (29-32 wGA) and functional connectivity (FC) was estimated from the scalp EEG. TTA-SW, particularly in the theta band, induced increased FC between left temporal and left frontal areas and between left temporal and parietal areas with TTA-SW at the left temporal region, while FC was limited to the right temporal regions in the case of TTA-SW at the right temporal region. Regardless of the lateralization of TTA-SW, long-range FCs were observed between left frontal to left parietal areas, suggesting that these regions, together with the temporal region, provide a basis for coherent neuronal activation across distal cortical regions. TTA-SW dynamic features showed that brief phases of TTA-SW had an impact on both local and whole brain network organization, supporting the importance of TTA-SW as a biomarker of brain development.

Structural imaging studies of patients with chronic pain: an anatomical likelihood estimate meta-analysis.

ABSTRACT: Neuroimaging is a powerful tool to investigate potential associations between chronic pain and brain structure. However, the proliferation of studies across diverse chronic pain syndromes and heterogeneous results challenges data integration and interpretation. We conducted a preregistered anatomical likelihood estimate meta-analysis on structural magnetic imaging studies comparing patients with chronic pain and healthy controls. Specifically, we investigated a broad range of measures of brain structure as well as specific alterations in gray matter and cortical thickness. A total of 7849 abstracts of experiments published between January 1, 1990, and April 26, 2021, were identified from 8 databases and evaluated by 2 independent reviewers. Overall, 103 experiments with a total of 5075 participants met the preregistered inclusion criteria. After correction for multiple comparisons using the gold-standard family-wise error correction ( P < 0.05), no significant differences associated with chronic pain were found. However, exploratory analyses using threshold-free cluster enhancement revealed several spatially distributed clusters showing structural alterations in chronic pain. Most of the clusters coincided with regions implicated in nociceptive processing including the amygdala, thalamus, hippocampus, insula, anterior cingulate cortex, and inferior frontal gyrus. Taken together, these results suggest that chronic pain is associated with subtle, spatially distributed alterations of brain structure.

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

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

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

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

Intrinsic activity development unfolds along a sensorimotor-association cortical axis in youth.

Animal studies of neurodevelopment have shown that recordings of intrinsic cortical activity evolve from synchronized and high amplitude to sparse and low amplitude as plasticity declines and the cortex matures. Leveraging resting-state functional MRI (fMRI) data from 1,033 youths (ages 8-23 years), we find that this stereotyped refinement of intrinsic activity occurs during human development and provides evidence for a cortical gradient of neurodevelopmental change. Declines in the amplitude of intrinsic fMRI activity were initiated heterochronously across regions and were coupled to the maturation of intracortical myelin, a developmental plasticity regulator. Spatiotemporal variability in regional developmental trajectories was organized along a hierarchical, sensorimotor-association cortical axis from ages 8 to 18. The sensorimotor-association axis furthermore captured variation in associations between youths' neighborhood environments and intrinsic fMRI activity; associations suggest that the effects of environmental disadvantage on the maturing brain diverge most across this axis during midadolescence. These results uncover a hierarchical neurodevelopmental axis and offer insight into the progression of cortical plasticity in humans.

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.

A developmental reduction of the excitation:inhibition ratio in association cortex during adolescence.

Adolescence is hypothesized to be a critical period for the development of association cortex. A reduction of the excitation:inhibition (E:I) ratio is a hallmark of critical period development; however, it has been unclear how to assess the development of the E:I ratio using noninvasive neuroimaging techniques. Here, we used pharmacological fMRI with a GABAergic benzodiazepine challenge to empirically generate a model of E:I ratio based on multivariate patterns of functional connectivity. In an independent sample of 879 youth (ages 8 to 22 years), this model predicted reductions in the E:I ratio during adolescence, which were specific to association cortex and related to psychopathology. These findings support hypothesized shifts in E:I balance of association cortices during a neurodevelopmental critical period in adolescence.