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.

Imaging chronic active lesions in multiple sclerosis: a consensus statement.

Chronic active lesions (CAL) are an important manifestation of chronic inflammation in multiple sclerosis (MS) and have implications for non-relapsing biological progression. In recent years, the discovery of innovative magnetic resonance imaging (MRI) and PET derived biomarkers has made it possible to detect CAL, and to some extent quantify them, in the brain of persons with MS, in vivo. Paramagnetic rim lesions on susceptibility-sensitive MRI sequences, MRI-defined slowly expanding lesions on T1-weighted (T1-w) and T2-w scans, and 18-kDa translocator protein-positive lesions on PET are promising candidate biomarkers of CAL. While partially overlapping, these biomarkers do not have equivalent sensitivity and specificity to histopathological CAL. Standardization in the use of available imaging measures for CAL identification, quantification, and monitoring is lacking. To fast-forward clinical translation of CAL, the North American Imaging in Multiple Sclerosis Cooperative developed a Consensus Statement, which provides guidance for the radiological definition and measurement of CAL. The proposed manuscript presents this Consensus Statement, summarizes the multistep process leading to it, and identifies the remaining major gaps in knowledge.

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.

Network enrichment significance testing in brain-phenotype association studies.

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

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

Penalized decomposition using residuals (PeDecURe) for feature extraction in the presence of nuisance variables.

Neuroimaging data are an increasingly important part of etiological studies of neurological and psychiatric disorders. However, mitigating the influence of nuisance variables, including confounders, remains a challenge in image analysis. In studies of Alzheimer's disease, for example, an imbalance in disease rates by age and sex may make it difficult to distinguish between structural patterns in the brain (as measured by neuroimaging scans) attributable to disease progression and those characteristic of typical human aging or sex differences. Concerningly, when not properly accounted for, nuisance variables pose threats to the generalizability and interpretability of findings from these studies. Motivated by this critical issue, in this work, we examine the impact of nuisance variables on feature extraction methods and propose Penalized Decomposition Using Residuals (PeDecURe), a new method for obtaining nuisance variable-adjusted features. PeDecURe estimates primary directions of variation which maximize covariance between partially residualized imaging features and a variable of interest (e.g., Alzheimer's diagnosis) while simultaneously mitigating the influence of nuisance variation through a penalty on the covariance between partially residualized imaging features and those variables. Using features derived using PeDecURe's first direction of variation, we train a highly accurate and generalizable predictive model, as evidenced by its robustness in testing samples with different underlying nuisance variable distributions. We compare PeDecURe to commonly used decomposition methods (principal component analysis (PCA) and partial least squares) as well as a confounder-adjusted variation of PCA. We find that features derived from PeDecURe offer greater accuracy and generalizability and lower correlations with nuisance variables compared with the other methods. While PeDecURe is primarily motivated by challenges that arise in the analysis of neuroimaging data, it is broadly applicable to data sets with highly correlated features, where novel methods to handle nuisance variables are warranted.

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.

Effect of patient factors, center, and era on Fontan timing: An observational study using the Pediatric Health Information Systems Database.

BACKGROUND: There are no consensus guidelines defining optimal timing for the Fontan operation, the last planned surgery in staged palliation for single-ventricle heart disease. OBJECTIVES: Identify patient-level characteristics, center-level variation, and secular trends driving Fontan timing. METHODS: A retrospective observational study of subjects who underwent Fontan from 2007 to 2021 at centers in the Pediatric Health Information Systems database was performed using linear mixed-effects modeling in which age at Fontan was regressed on patient characteristics and date of operation with center as random effect. RESULTS: We included 10,305 subjects (40.4% female, 44% non-white) at 47 centers. Median age at Fontan was 3.4 years (IQR 2.6-4.4). Hypoplastic left heart syndrome (-4.4 months, 95%CI -5.5 to -3.3) and concomitant conditions (-2.6 months, 95%CI -4.1 to -1.1) were associated with younger age at Fontan. Subjects with technology-dependence (+4.6 months, 95%CI 3.1-6.1) were older at Fontan. Black (+4.1 months, 95%CI 2.5-5.7) and Asian (+8.3 months, 95%CI 5.4-11.2) race were associated with older age at Fontan. There was significant variation in Fontan timing between centers. Center accounted for 10% of variation (ICC 0.10, 95%CI 0.07-0.14). Center surgical volume was not associated with Fontan timing (P = .21). Operation year was associated with age at Fontan, with a 3.1 month increase in age for every 5 years (+0.61 months, 95%CI 0.48-0.75). CONCLUSIONS: After adjusting for patient-level characteristics there remains significant inter-center variation in Fontan timing. Age at Fontan has increased. Future studies addressing optimal Fontan timing are warranted.

SpaGCN: Integrating gene expression, spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network.

Recent advances in spatially resolved transcriptomics (SRT) technologies have enabled comprehensive characterization of gene expression patterns in the context of tissue microenvironment. To elucidate spatial gene expression variation, we present SpaGCN, a graph convolutional network approach that integrates gene expression, spatial location and histology in SRT data analysis. Through graph convolution, SpaGCN aggregates gene expression of each spot from its neighboring spots, which enables the identification of spatial domains with coherent expression and histology. The subsequent domain guided differential expression (DE) analysis then detects genes with enriched expression patterns in the identified domains. Analyzing seven SRT datasets using SpaGCN, we show it can detect genes with much more enriched spatial expression patterns than competing methods. Furthermore, genes detected by SpaGCN are transferrable and can be utilized to study spatial variation of gene expression in other datasets. SpaGCN is computationally fast, platform independent, making it a desirable tool for diverse SRT studies.

Early Magnetic Resonance Imaging Features of New Paramagnetic Rim Lesions in Multiple Sclerosis.

OBJECTIVE: To determine early magnetic resonance imaging (MRI) features of new multiple sclerosis (MS) lesions that will develop into paramagnetic rim lesions (PRLs), which have been associated with progressive tissue injury in MS. METHODS: New contrast-enhancing lesions observed on routine clinical MRI were imaged at 7 T within 4 weeks of observation, and 3 and 6 months later. The 6-month MRI was used to classify PRL status (PRL or non-PRL). The relationship between early lesion characteristics and subsequent PRL status was assessed using generalized linear mixed effects models. Random forest classification was performed to classify early predictors of subsequent PRL status. RESULTS: From 93 contrast-enhancing lesions in 23 MS patients, 37 lesions developed into a PRL. In lesions that developed into PRLs compared with those that did not, the average lesion T1 on the initial 7 T MRI was 1994 ms compared with 1,670 ms (p-value <0.001), and the average volume was 168.7 mL compared with 44 mL (p-value <0.001) in lesions that did not. These volume differences were also found on 3 T scans (p-value <0.001), and for intensity-normalized T1 -w (p-value = 0.011) and fluid-attenuated inversion recovery (p-value = 0.005). The area under the receiver operating characteristic curve for the random forest classification with leave-one-out cross-validation was found to be 0.86 using initial 7 T features. INTERPRETATION: New MS lesions that evolve into PRLs can be identified early in lesion evolution. These findings suggest that biological mechanisms underlying PRL development begin early, which has important implications for clinical trials targeting PRLs development and subsequent therapeutics. ANN NEUROL 2023;94:736-744.

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.

Choroid plexus volume differentiates MS from its mimics.

This study aimed to determine whether choroid plexus volume (CPV) could differentiate multiple sclerosis (MS) from its mimics. A secondary analysis of two previously enrolled studies, 50 participants with MS and 64 with alternative diagnoses were included. CPV was automatically segmented from 3T magnetic resonance imaging (MRI), followed by manual review to remove misclassified tissue. Mean normalized choroid plexus volume (nCPV) to intracranial volume demonstrated relatively high specificity for MS participants in each cohort (0.80 and 0.76) with an area under the receiver-operator characteristic curve of 0.71 (95% confidence interval (CI) = 0.55-0.87) and 0.65 (95% CI = 0.52-0.77). In this preliminary study, nCPV differentiated MS from its mimics.

A multicenter pilot study evaluating simplified central vein assessment for the diagnosis of multiple sclerosis.

BACKGROUND: The central vein sign (CVS) is a proposed magnetic resonance imaging (MRI) biomarker for multiple sclerosis (MS); the optimal method for abbreviated CVS scoring is not yet established. OBJECTIVE: The aim of this study was to evaluate the performance of a simplified approach to CVS assessment in a multicenter study of patients being evaluated for suspected MS. METHODS: Adults referred for possible MS to 10 sites were recruited. A post-Gd 3D T2*-weighted MRI sequence (FLAIR*) was obtained in each subject. Trained raters at each site identified up to six CVS-positive lesions per FLAIR* scan. Diagnostic performance of CVS was evaluated for a diagnosis of MS which had been confirmed using the 2017 McDonald criteria at thresholds including three positive lesions (Select-3*) and six positive lesions (Select-6*). Inter-rater reliability assessments were performed. RESULTS: Overall, 78 participants were analyzed; 37 (47%) were diagnosed with MS, and 41 (53%) were not. The mean age of participants was 45 (range: 19-64) years, and most were female (n = 55, 71%). The area under the receiver operating characteristic curve (AUROC) for the simplified counting method was 0.83 (95% CI: 0.73-0.93). Select-3* and Select-6* had sensitivity of 81% and 65% and specificity of 68% and 98%, respectively. Inter-rater agreement was 78% for Select-3* and 83% for Select-6*. CONCLUSION: A simplified method for CVS assessment in patients referred for suspected MS demonstrated good diagnostic performance and inter-rater agreement.

Alcohol for seizure induction in the epilepsy monitoring unit.

RATIONALE: Seizure induction techniques are used in the epilepsy monitoring unit (EMU) to increase diagnostic yield and reduce length of stay. There are insufficient data on the efficacy of alcohol as an induction technique. METHODS: We performed a retrospective cohort study using six years of EMU data at our institution. We compared cases who received alcohol for seizure induction to matched controls who did not. The groups were matched on the following variables: age, reason for admission, length of stay, number of antiseizure medications (ASM) at admission, whether ASMs were tapered during admission, and presence of interictal epileptiform discharges. We used both propensity score and exact matching strategies. We compared the likelihood of epileptic seizures and nonepileptic events in cases versus controls using Kaplan-Meier time-to-event analysis, as well as odds ratios for these outcomes occurring at any time during the admission. RESULTS: We analyzed 256 cases who received alcohol (median dose 2.5 standard drinks) and 256 propensity score-matched controls. Cases who received alcohol were no more likely than controls to have an epileptic seizure (X2(1) = 0.01, p = 0.93) or nonepileptic event (X2(1) = 2.1, p = 0.14) in the first 48 h after alcohol administration. For the admission overall, cases were no more likely to have an epileptic seizure (OR 0.89, 95 % CI 0.61-1.28, p = 0.58), nonepileptic event (OR 0.97, CI 0.62-1.53, p = 1.00), nor require rescue benzodiazepine (OR 0.63, CI 0.35-1.12, p = 0.15). Stratified analyses revealed no increased risk of epileptic seizure in any subgroups. Sensitivity analysis using exact matching showed that results were robust to matching strategy. CONCLUSIONS: Alcohol was not an effective induction technique in the EMU. This finding has implications for counseling patients with epilepsy about the risks of drinking alcohol in moderation in their daily lives.

Neural Responses to Auditory Food Stimuli Following Cognitive Behavioral Therapy for Binge-Eating Disorder.

OBJECTIVE: Adults with binge-eating disorder (BED), compared with those without BED, demonstrate higher blood-oxygen-level-dependent (BOLD) response to food cues in reward-related regions of the brain. It is not known whether cognitive behavioral therapy (CBT) can reverse this reward system hyperactivation. This randomized controlled trial (RCT) assessed changes in BOLD response to binge-eating cues following CBT versus wait-list control (WLC). METHOD: Females with BED (N = 40) were randomized to CBT or WLC. Participants completed assessments at baseline and 16 weeks including measures of eating and appetite and functional magnetic resonance imaging (fMRI) to measure BOLD response while listening to personalized scripts of binge-eating and neutral-relaxing cues. Data were analyzed using general linear models with mixed effects. RESULTS: Overall retention rate was 87.5%. CBT achieved significantly greater reductions in binge-eating episodes than WLC (mean ± standard error decline of 14.6 ± 2.7 vs. 5.7 ± 2.8 episodes in the past 28 days, respectively; p = 0.03). CBT and WLC did not differ significantly in changes in neural responses to binge-eating stimuli during the fMRI sessions. Compared with WLC, CBT had significantly greater improvements in reward-based eating drive, disinhibition, and hunger as assessed by questionnaires (ps < 0.05). DISCUSSION: CBT was effective in reducing binge eating, but, contrary to our hypothesis, CBT did not improve BOLD response to auditory binge-eating stimuli in reward regions of the brain. Further studies are needed to assess mechanisms underlying improvements with CBT for BED. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT03604172.

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.

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.

Matrix decomposition for modeling lesion development processes in multiple sclerosis.

Our main goal is to study and quantify the evolution of multiple sclerosis lesions observed longitudinally over many years in multi-sequence structural magnetic resonance imaging (sMRI). To achieve that, we propose a class of functional models for capturing the temporal dynamics and spatial distribution of the voxel-specific intensity trajectories in all sMRI sequences. To accommodate the hierarchical data structure (observations nested within voxels, which are nested within lesions, which, in turn, are nested within study participants), we use structured functional principal component analysis. We propose and evaluate the finite sample properties of hypothesis tests of therapeutic intervention effects on lesion evolution while accounting for the multilevel structure of the data. Using this novel testing strategy, we found statistically significant differences in lesion evolution between treatment groups.

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.

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.