Reproducible brain-wide association studies require thousands of individuals.

Magnetic resonance imaging (MRI) has transformed our understanding of the human brain through well-replicated mapping of abilities to specific structures (for example, lesion studies) and functions1-3 (for example, task functional MRI (fMRI)). Mental health research and care have yet to realize similar advances from MRI. A primary challenge has been replicating associations between inter-individual differences in brain structure or function and complex cognitive or mental health phenotypes (brain-wide association studies (BWAS)). Such BWAS have typically relied on sample sizes appropriate for classical brain mapping4 (the median neuroimaging study sample size is about 25), but potentially too small for capturing reproducible brain-behavioural phenotype associations5,6. Here we used three of the largest neuroimaging datasets currently available-with a total sample size of around 50,000 individuals-to quantify BWAS effect sizes and reproducibility as a function of sample size. BWAS associations were smaller than previously thought, resulting in statistically underpowered studies, inflated effect sizes and replication failures at typical sample sizes. As sample sizes grew into the thousands, replication rates began to improve and effect size inflation decreased. More robust BWAS effects were detected for functional MRI (versus structural), cognitive tests (versus mental health questionnaires) and multivariate methods (versus univariate). Smaller than expected brain-phenotype associations and variability across population subsamples can explain widespread BWAS replication failures. In contrast to non-BWAS approaches with larger effects (for example, lesions, interventions and within-person), BWAS reproducibility requires samples with thousands of individuals.

Cumulative Effects of Resting-State Connectivity Across All Brain Networks Significantly Correlate with Attention-Deficit Hyperactivity Disorder Symptoms.

Identification of replicable neuroimaging correlates of attention-deficit hyperactivity disorder (ADHD) has been hindered by small sample sizes, small effects, and heterogeneity of methods. Given evidence that ADHD is associated with alterations in widely distributed brain networks and the small effects of individual brain features, a whole-brain perspective focusing on cumulative effects is warranted. The use of large, multisite samples is crucial for improving reproducibility and clinical utility of brain-wide MRI association studies. To address this, a polyneuro risk score (PNRS) representing cumulative, brain-wide, ADHD-associated resting-state functional connectivity was constructed and validated using data from the Adolescent Brain Cognitive Development (ABCD, N = 5,543, 51.5% female) study, and was further tested in the independent Oregon-ADHD-1000 case-control cohort (N = 553, 37.4% female). The ADHD PNRS was significantly associated with ADHD symptoms in both cohorts after accounting for relevant covariates (p < 0.001). The most predictive PNRS involved all brain networks, though the strongest effects were concentrated among the default mode and cingulo-opercular networks. In the longitudinal Oregon-ADHD-1000, non-ADHD youth had significantly lower PNRS (Cohen's d = -0.318, robust p = 5.5 × 10-4) than those with persistent ADHD (age 7-19). The PNRS, however, did not mediate polygenic risk for ADHD. Brain-wide connectivity was robustly associated with ADHD symptoms in two independent cohorts, providing further evidence of widespread dysconnectivity in ADHD. Evaluation in enriched samples demonstrates the promise of the PNRS approach for improving reproducibility in neuroimaging studies and unraveling the complex relationships between brain connectivity and behavioral disorders.

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.

Early-in-life isoflurane exposure alters resting-state functional connectivity in juvenile non-human primates.

BACKGROUND: Clinical studies suggest that anaesthesia exposure early in life affects neurobehavioural development. We designed a non-human primate (NHP) study to evaluate cognitive, behavioural, and brain functional and structural alterations after isoflurane exposure during infancy. These NHPs displayed decreased close social behaviour and increased astrogliosis in specific brain regions, most notably in the amygdala. Here we hypothesise that resting-state functional connectivity MRI can detect alterations in connectivity of brain areas that relate to these social behaviours and astrogliosis. METHODS: Imaging was performed in 2-yr-old NHPs under light anaesthesia, after early-in-life (postnatal days 6-12) exposure to 5 h of isoflurane either one or three times, or to room air. Brain images were segmented into 82 regions of interest; the amygdala and the posterior cingulate cortex were chosen for a seed-based resting-state functional connectivity MRI analysis. RESULTS: We found differences between groups in resting-state functional connectivity of the amygdala and the auditory cortices, medial premotor cortex, and posterior cingulate cortex. There were also alterations in resting-state functional connectivity between the posterior cingulate cortex and secondary auditory, polar prefrontal, and temporal cortices, and the anterior insula. Relationships were identified between resting-state functional connectivity alterations and the decrease in close social behaviour and increased astrogliosis. CONCLUSIONS: Early-in-life anaesthesia exposure in NHPs is associated with resting-state functional connectivity alterations of the amygdala and the posterior cingulate cortex with other brain regions, evident at the juvenile age of 2 yr. These changes in resting-state functional connectivity correlate with the decrease in close social behaviour and increased astrogliosis. Using resting-state functional connectivity MRI to study the neuronal underpinnings of early-in-life anaesthesia-induced behavioural alterations could facilitate development of a biomarker for anaesthesia-induced developmental neurotoxicity.

An open-access accelerated adult equivalent of the ABCD Study neuroimaging dataset (a-ABCD).

As public access to longitudinal developmental datasets like the Adolescent Brain Cognitive Development StudySM (ABCD Study®) increases, so too does the need for resources to benchmark time-dependent effects. Scan-to-scan changes observed with repeated imaging may reflect development but may also reflect practice effects, day-to-day variability in psychological states, and/or measurement noise. Resources that allow disentangling these time-dependent effects will be useful in quantifying actual developmental change. We present an accelerated adult equivalent of the ABCD Study dataset (a-ABCD) using an identical imaging protocol to acquire magnetic resonance imaging (MRI) structural, diffusion-weighted, resting-state and task-based data from eight adults scanned five times over five weeks. We report on the task-based imaging data (n = 7). In-scanner stop-signal (SST), monetary incentive delay (MID), and emotional n-back (EN-back) task behavioral performance did not change across sessions. Post-scan recognition memory for emotional n-back stimuli, however, did improve as participants became more familiar with the stimuli. Functional MRI analyses revealed that patterns of task-based activation reflecting inhibitory control in the SST, reward success in the MID task, and working memory in the EN-back task were more similar within individuals across repeated scan sessions than between individuals. Within-subject, activity was more consistent across sessions during the EN-back task than in the SST and MID task, demonstrating differences in fMRI data reliability as a function of task. The a-ABCD dataset provides a unique testbed for characterizing the reliability of brain function, structure, and behavior across imaging modalities in adulthood and benchmarking neurodevelopmental change observed in the open-access ABCD Study.

Synthesizing pseudo-T2w images to recapture missing data in neonatal neuroimaging with applications in rs-fMRI.

T1- and T2-weighted (T1w and T2w) images are essential for tissue classification and anatomical localization in Magnetic Resonance Imaging (MRI) analyses. However, these anatomical data can be challenging to acquire in non-sedated neonatal cohorts, which are prone to high amplitude movement and display lower tissue contrast than adults. As a result, one of these modalities may be missing or of such poor quality that they cannot be used for accurate image processing, resulting in subject loss. While recent literature attempts to overcome these issues in adult populations using synthetic imaging approaches, evaluation of the efficacy of these methods in pediatric populations and the impact of these techniques in conventional MR analyses has not been performed. In this work, we present two novel methods to generate pseudo-T2w images: the first is based in deep learning and expands upon previous models to 3D imaging without the requirement of paired data, the second is based in nonlinear multi-atlas registration providing a computationally lightweight alternative. We demonstrate the anatomical accuracy of pseudo-T2w images and their efficacy in existing MR processing pipelines in two independent neonatal cohorts. Critically, we show that implementing these pseudo-T2w methods in resting-state functional MRI analyses produces virtually identical functional connectivity results when compared to those resulting from T2w images, confirming their utility in infant MRI studies for salvaging otherwise lost subject data.

Filtering respiratory motion artifact from resting state fMRI data in infant and toddler populations.

The importance of motion correction when processing resting state functional magnetic resonance imaging (rs-fMRI) data is well-established in adult cohorts. This includes adjustments based on self-limited, large amplitude subject head motion, as well as factitious rhythmic motion induced by respiration. In adults, such respiration artifact can be effectively removed by applying a notch filter to the motion trace, resulting in higher amounts of data retained after frame censoring (e.g., "scrubbing") and more reliable correlation values. Due to the unique physiological and behavioral characteristics of infants and toddlers, rs-fMRI processing pipelines, including methods to identify and remove colored noise due to subject motion, must be appropriately modified to accurately reflect true neuronal signal. These younger cohorts are characterized by higher respiration rates and lower-amplitude head movements than adults; thus, the presence and significance of comparable respiratory artifact and the subsequent necessity of applying similar techniques remain unknown. Herein, we identify and characterize the consistent presence of respiratory artifact in rs-fMRI data collected during natural sleep in infants and toddlers across two independent cohorts (aged 8-24 months) analyzed using different pipelines. We further demonstrate how removing this artifact using an age-specific notch filter allows for both improved data quality and data retention in measured results. Importantly, this work reveals the critical need to identify and address respiratory-driven head motion in fMRI data acquired in young populations through the use of age-specific motion filters as a mechanism to optimize the accuracy of measured results in this population.

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.

Maternal Interleukin-6 Is Associated With Macaque Offspring Amygdala Development and Behavior.

Human and animal cross-sectional studies have shown that maternal levels of the inflammatory cytokine interleukin-6 (IL-6) may compromise brain phenotypes assessed at single time points. However, how maternal IL-6 associates with the trajectory of brain development remains unclear. We investigated whether maternal IL-6 levels during pregnancy relate to offspring amygdala volume development and anxiety-like behavior in Japanese macaques. Magnetic resonance imaging (MRI) was administered to 39 Japanese macaque offspring (Female: 18), providing at least one or more time points at 4, 11, 21, and 36 months of age with a behavioral assessment at 11 months of age. Increased maternal third trimester plasma IL-6 levels were associated with offspring's smaller left amygdala volume at 4 months, but with more rapid amygdala growth from 4 to 36 months. Maternal IL-6 predicted offspring anxiety-like behavior at 11 months, which was mediated by reduced amygdala volumes in the model's intercept (i.e., 4 months). The results increase our understanding of the role of maternal inflammation in the development of neurobehavioral disorders by detailing the associations of a commonly examined inflammatory indicator, IL-6, on amygdala volume growth over time, and anxiety-like behavior.

Correction of respiratory artifacts in MRI head motion estimates.

Head motion represents one of the greatest technical obstacles in magnetic resonance imaging (MRI) of the human brain. Accurate detection of artifacts induced by head motion requires precise estimation of movement. However, head motion estimates may be corrupted by artifacts due to magnetic main field fluctuations generated by body motion. In the current report, we examine head motion estimation in multiband resting state functional connectivity MRI (rs-fcMRI) data from the Adolescent Brain and Cognitive Development (ABCD) Study and comparison 'single-shot' datasets. We show that respirations contaminate movement estimates in functional MRI and that respiration generates apparent head motion not associated with functional MRI quality reductions. We have developed a novel approach using a band-stop filter that accurately removes these respiratory effects from motion estimates. Subsequently, we demonstrate that utilizing a band-stop filter improves post-processing fMRI data quality. Lastly, we demonstrate the real-time implementation of motion estimate filtering in our FIRMM (Framewise Integrated Real-Time MRI Monitoring) software package.

Neural Circuit Markers of Familial Risk for Depression Among Healthy Youths in the Adolescent Brain Cognitive Development Study.

BACKGROUND: Family history of depression is a robust predictor of early-onset depression, which may confer risk through alterations in neural circuits that have been implicated in reward and emotional processing. These alterations may be evident in youths who are at familial risk for depression but who do not currently have depression. However, the identification of robust and replicable findings has been hindered by few studies and small sample sizes. In the current study, we sought to identify functional connectivity (FC) patterns associated with familial risk for depression. METHODS: Participants included healthy (i.e., no lifetime psychiatric diagnoses) youths at high familial risk for depression (HR) (n = 754; at least one parent with a history of depression) and healthy youths at low familial risk for psychiatric problems (LR) (n = 1745; no parental history of psychopathology) who were 9 to 10 years of age and from the Adolescent Brain Cognitive Development (ABCD) Study sample. We conducted whole-brain seed-to-voxel analyses to examine group differences in resting-state FC with the amygdala, caudate, nucleus accumbens, and putamen. We hypothesized that HR youths would exhibit global amygdala hyperconnectivity and striatal hypoconnectivity patterns primarily driven by maternal risk. RESULTS: HR youths exhibited weaker caudate-angular gyrus FC than LR youths (α = 0.04, Cohen's d = 0.17). HR youths with a history of maternal depression specifically exhibited weaker caudate-angular gyrus FC (α = 0.03, Cohen's d = 0.19) as well as weaker caudate-dorsolateral prefrontal cortex FC (α = 0.04, Cohen's d = 0.21) than LR youths. CONCLUSIONS: Weaker striatal connectivity may be related to heightened familial risk for depression, primarily driven by maternal history. Identifying brain-based markers of depression risk in youths can inform approaches to improving early detection, diagnosis, and treatment.

Familial risk for depression moderates neural circuitry in healthy preadolescents to predict adolescent depression symptoms in the Adolescent Brain Cognitive Development (ABCD) Study.

BACKGROUND: There is an imminent need to identify neural markers during preadolescence that are linked to developing depression during adolescence, especially among youth at elevated familial risk. However, longitudinal studies remain scarce and exhibit mixed findings. Here we aimed to elucidate functional connectivity (FC) patterns among preadolescents that interact with familial depression risk to predict depression two years later. METHODS: 9-10 year-olds in the Adolescent Brain Cognitive Development (ABCD) Study were classified as healthy (i.e., no lifetime psychiatric diagnoses) at high familial risk for depression (HR; n=559) or at low familial risk for psychopathology (LR; n=1203). Whole-brain seed-to-voxel resting-state FC patterns with the amygdala, putamen, nucleus accumbens, and caudate were calculated. Multi-level, mixed-effects regression analyses were conducted to test whether FC at ages 9-10 interacted with familial risk to predict depression symptoms at ages 11-12. RESULTS: HR youth demonstrated stronger associations between preadolescent FC and adolescent depression symptoms (ps<0.001) as compared to LR youth (ps>0.001), primarily among amygdala/striatal FC with visual and sensory/somatomotor networks. CONCLUSIONS: Preadolescent amygdala and striatal FC may be useful biomarkers of adolescent-onset depression, particularly for youth with family histories of depression. This research may point to neurobiologically-informed approaches to prevention and intervention for depression in adolescents.

Multi-level fMRI analysis applied to hemispheric specialization in the language network, functional areas, and their behavioral correlations in the ABCD sample.

Prior research suggests that the organization of the language network in the brain is left-dominant and becomes more lateralized with age and increasing language skill. The age at which specific components of the language network become adult-like varies depending on the abilities they subserve. So far, a large, developmental study has not included a language task paradigm, so we introduce a method to study resting-state laterality in the Adolescent Brain Cognitive Development (ABCD) study. Our approach mixes source timeseries between left and right homotopes of the (1) inferior frontal and (2) middle temporal gyri and (3) a region we term "Wernicke's area" near the supramarginal gyrus. Our large subset sample size of ABCD (n = 6153) allows improved reliability and validity compared to previous, smaller studies of brain-behavior associations. We show that behavioral metrics from the NIH Youth Toolbox and other resources are differentially related to tasks with a larger linguistic component over ones with less (e.g., executive function-dominant tasks). These baseline characteristics of hemispheric specialization in youth are critical for future work determining the correspondence of lateralization with language onset in earlier stages of development.

Heritability estimation of subcortical volumes in a multi-ethnic multi-site cohort study.

Heritability of regional subcortical brain volumes (rSBVs) describes the role of genetics in middle and inner brain development. rSBVs are highly heritable in adults but are not characterized well in adolescents. The Adolescent Brain Cognitive Development study (ABCD), taken over 22 US sites, provides data to characterize the heritability of subcortical structures in adolescence. In ABCD, site-specific effects co-occur with genetic effects which can bias heritability estimates. Existing methods adjusting for site effects require additional steps to adjust for site effects and can lead to inconsistent estimation. We propose a random-effect model-based method of moments approach that is a single step estimator and is a theoretically consistent estimator even when sites are imbalanced and performs well under simulations. We compare methods on rSBVs from ABCD. The proposed approach yielded heritability estimates similar to previous results derived from single-site studies. The cerebellum cortex and hippocampus were the most heritable regions (> 50%).

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.

Polyneuro risk scores capture widely distributed connectivity patterns of cognition.

Resting-state functional connectivity (RSFC) is a powerful tool for characterizing brain changes, but it has yet to reliably predict higher-order cognition. This may be attributed to small effect sizes of such brain-behavior relationships, which can lead to underpowered, variable results when utilizing typical sample sizes (N∼25). Inspired by techniques in genomics, we implement the polyneuro risk score (PNRS) framework - the application of multivariate techniques to RSFC data and validation in an independent sample. Utilizing the Adolescent Brain Cognitive Development® cohort split into two datasets, we explore the framework's ability to reliably capture brain-behavior relationships across 3 cognitive scores - general ability, executive function, learning & memory. The weight and significance of each connection is assessed in the first dataset, and a PNRS is calculated for each participant in the second. Results support the PNRS framework as a suitable methodology to inspect the distribution of connections contributing towards behavior, with explained variance ranging from 1.0 % to 21.4 %. For the outcomes assessed, the framework reveals globally distributed, rather than localized, patterns of predictive connections. Larger samples are likely necessary to systematically identify the specific connections contributing towards complex outcomes. The PNRS framework could be applied translationally to identify neurologically distinct subtypes of neurodevelopmental disorders.

BIBSNet: A Deep Learning Baby Image Brain Segmentation Network for MRI Scans.

Objectives: Brain segmentation of infant magnetic resonance (MR) images is vitally important in studying developmental mental health and disease. The infant brain undergoes many changes throughout the first years of postnatal life, making tissue segmentation difficult for most existing algorithms. Here, we introduce a deep neural network BIBSNet (Baby and Infant Brain Segmentation Neural Network), an open-source, community-driven model that relies on data augmentation and a large sample size of manually annotated images to facilitate the production of robust and generalizable brain segmentations. Experimental Design: Included in model training and testing were MR brain images on 84 participants with an age range of 0-8 months (median postmenstrual ages of 13.57 months). Using manually annotated real and synthetic segmentation images, the model was trained using a 10-fold cross-validation procedure. Testing occurred on MRI data processed with the DCAN labs infant-ABCD-BIDS processing pipeline using segmentations produced from gold standard manual annotation, joint-label fusion (JLF), and BIBSNet to assess model performance. Principal Observations: Using group analyses, results suggest that cortical metrics produced using BIBSNet segmentations outperforms JLF segmentations. Additionally, when analyzing individual differences, BIBSNet segmentations perform even better. Conclusions: BIBSNet segmentation shows marked improvement over JLF segmentations across all age groups analyzed. The BIBSNet model is 600x faster compared to JLF and can be easily included in other processing pipelines.

Using synthetic MR images for distortion correction.

Functional MRI (fMRI) data acquired using echo-planar imaging (EPI) are highly distorted by magnetic field inhomogeneities. Distortion and differences in image contrast between EPI and T1-weighted and T2-weighted (T1w/T2w) images makes their alignment a challenge. Typically, field map data are used to correct EPI distortions. Alignments achieved with field maps can vary greatly and depends on the quality of field map data. However, many public datasets lack field map data entirely. Additionally, reliable field map data is often difficult to acquire in high-motion pediatric or developmental cohorts. To address this, we developed Synth, a software package for distortion correction and cross-modal image registration that does not require field map data. Synth combines information from T1w and T2w anatomical images to construct an idealized undistorted synthetic image with similar contrast properties to EPI data. This synthetic image acts as an effective reference for individual-specific distortion correction. Using pediatric (ABCD: Adolescent Brain Cognitive Development) and adult (MSC: Midnight Scan Club; HCP: Human Connectome Project) data, we demonstrate that Synth performs comparably to field map distortion correction approaches, and often outperforms them. Field map-less distortion correction with Synth allows accurate and precise registration of fMRI data with missing or corrupted field map information.

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.