Comparison of traveling-subject and ComBat harmonization methods for assessing structural brain characteristics

Noninvasive Magnetic Resonance Imaging Measures of Glymphatic System Activity

The comprehension of the glymphatic system, a postulated mechanism responsible for the removal of interstitial solutes within the central nervous system (CNS), has witnessed substantial progress recently. While direct measurement techniques involving fluorescence and contrast agent tracers have demonstrated success in animal studies, their application in humans is invasive and presents challenges. Hence, exploring alternative noninvasive approaches that enable glymphatic research in humans is imperative. This review primarily focuses on several noninvasive magnetic resonance imaging (MRI) techniques, encompassing perivascular space (PVS) imaging, diffusion tensor image analysis along the PVS, arterial spin labeling, chemical exchange saturation transfer, and intravoxel incoherent motion. These methodologies provide valuable insights into the dynamics of interstitial fluid, water permeability across the blood-brain barrier, and cerebrospinal fluid flow within the cerebral parenchyma. Furthermore, the review elucidates the underlying concept and clinical applications of these noninvasive MRI techniques, highlighting their strengths and limitations. It addresses concerns about the relationship between glymphatic system activity and pathological alterations, emphasizing the necessity for further studies to establish correlations between noninvasive MRI measurements and pathological findings. Additionally, the challenges associated with conducting multisite studies, such as variability in MRI systems and acquisition parameters, are addressed, with a suggestion for the use of harmonization methods, such as the combined association test (COMBAT), to enhance standardization and statistical power. Current research gaps and future directions in noninvasive MRI techniques for assessing the glymphatic system are discussed, emphasizing the need for larger sample sizes, harmonization studies, and combined approaches. In conclusion, this review provides invaluable insights into the application of noninvasive MRI methods for monitoring glymphatic system activity in the CNS. It highlights their potential in advancing our understanding of the glymphatic system, facilitating clinical applications, and paving the way for future research endeavors in this field. EVIDENCE LEVEL: 3 TECHNICAL EFFICACY: Stage 5.

Combining Clinical With Cognitive or Magnetic Resonance Imaging Data for Predicting Transition to Psychosis in Ultra High-Risk Patients: Data From the PACE 400 Cohort

BACKGROUND

Multimodal modeling that combines biological and clinical data shows promise in predicting transition to psychosis in individuals who are at ultra-high risk. Individuals who transition to psychosis are known to have deficits at baseline in cognitive function and reductions in gray matter volume in multiple brain regions identified by magnetic resonance imaging.

METHODS

In this study, we used Cox proportional hazards regression models to assess the additive predictive value of each modality-cognition, cortical structure information, and the neuroanatomical measure of brain age gap-to a previously developed clinical model using functioning and duration of symptoms prior to service entry as predictors in the Personal Assessment and Crisis Evaluation (PACE) 400 cohort. The PACE 400 study is a well-characterized cohort of Australian youths who were identified as ultra-high risk of transitioning to psychosis using the Comprehensive Assessment of At Risk Mental States (CAARMS) and followed for up to 18 years; it contains clinical data (from N = 416 participants), cognitive data (n = 213), and magnetic resonance imaging cortical parameters extracted using FreeSurfer (n = 231).

RESULTS

The results showed that neuroimaging, brain age gap, and cognition added marginal predictive information to the previously developed clinical model (fraction of new information: neuroimaging 0%-12%, brain age gap 7%, cognition 0%-16%).

CONCLUSIONS

In summary, adding a second modality to a clinical risk model predicting the onset of a psychotic disorder in the PACE 400 cohort showed little improvement in the fit of the model for long-term prediction of transition to psychosis.

Harmonizing T1-Weighted Images to Improve Consistency of Brain Morphology Among Different Scanner Manufacturers in Alzheimer's disease

BACKGROUND

Brain MRI scanner variability can introduce bias in measurements. Harmonizing scanner variability is crucial.

PURPOSE

To develop a harmonization method aimed at removing scanner variability, and to evaluate the consistency of results in multicenter studies.

STUDY TYPE

Retrospective.

POPULATION

Multicenter data from 170 healthy participants (males/females = 98/72; age = 73.8 ± 7.3) and 170 Alzheimer's disease patients (males/females = 98/72; age = 76.2 ± 8.5) were compared with reference data from another 340 participants.

FIELD STRENGTH/SEQUENCE

3-T, magnetization prepared rapid gradient echo and turbo field echo; 1.5-T, inversion recovery prepared fast spoiled gradient echo T1-weighted sequences.

ASSESSMENT

Gray matter (GM) brain images, obtained through segmentation of T1-weighted images, were utilized to evaluate the performance of the harmonization method using common orthogonal basis extraction (HCOBE) and four other methods (removal of artificial voxel effect by linear regression, RAVEL; Z_score; general linear model, GLM; ComBat). Linear discriminant analysis (LDA) was used to access the effectiveness of different methods in reducing scanner variability. The performance of harmonization methods in preserving GM volumes heterogeneity was evaluated by the similarity of the relationship between GM proportion and age in the reference and multicenter data. Furthermore, the consistency of the harmonized multicenter data with the reference data were evaluated based on classification results (train/test = 7/3) and brain atrophy.

STATISTICAL TESTS

Two-sample t-tests, area under the curve (AUC), and Dice coefficients were used to analyze the consistency of results from the reference and harmonized multicenter data. A P-value <0.01 was considered statistically significant.

RESULTS

HCOBE reduced the scanner variability from 0.09 before harmonization to 0.003 (ideal: 0, RAVEL/Z_score/GLM/ComBat = 0.087/0.003/0.006/0.013). GM volumes showed no significant difference (P = 0.52) between the reference and HCOBE-harmonized multicenter data. Consistency evaluation showed that AUC values of 0.95 for both reference and HCOBE-harmonized multicenter data (RAVEL/Z_score/GLM/ComBat = 0.86/0.86/0.84/0.89), and the Dice coefficient increased from 0.73 before harmonization to 0.82 (ideal: 1, RAVEL/Z_score/GLM/ComBat = 0.39/0.64/0.59/0.74).

DATA CONCLUSION

HCOBE may help to remove scanner variability and could improve the consistency of results in multicenter studies.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Harmonizing multisite data with the ComBat method for enhanced Parkinson's disease diagnosis via DAT-SPECT

Background

Dopamine transporter single-photon emission computed tomography (DAT-SPECT) is a crucial tool for evaluating patients with Parkinson's disease (PD). However, its implication is limited by inter-site variability in large multisite clinical trials. To overcome the limitation, a conventional prospective correction method employs linear regression with phantom scanning, which is effective yet available only in a prospective manner. An alternative, although relatively underexplored, involves retrospective modeling using a statistical method known as "combatting batch effects when combining batches of gene expression microarray data" (ComBat).

Methods

We analyzed DAT-SPECT-specific binding ratios (SBRs) derived from 72 healthy older adults and 81 patients with PD registered in four clinical sites. We applied both the prospective correction and the retrospective ComBat correction to the original SBRs. Next, we compared the performance of the original and two corrected SBRs to differentiate the PD patients from the healthy controls. Diagnostic accuracy was assessed using the area under the receiver operating characteristic curve (AUC-ROC).

Results

The original SBRs were 6.13 ± 1.54 (mean ± standard deviation) and 2.03 ± 1.41 in the control and PD groups, respectively. After the prospective correction, the mean SBRs were 6.52 ± 1.06 and 2.40 ± 0.99 in the control and PD groups, respectively. After the retrospective ComBat correction, the SBRs were 5.25 ± 0.89 and 2.01 ± 0.73 in the control and PD groups, respectively, resulting in substantial changes in mean values with fewer variances. The original SBRs demonstrated fair performance in differentiating PD from controls (Hedges's g = 2.76; AUC-ROC = 0.936). Both correction methods improved discrimination performance. The ComBat-corrected SBR demonstrated comparable performance (g = 3.99 and AUC-ROC = 0.987) to the prospectively corrected SBR (g = 4.32 and AUC-ROC = 0.992) for discrimination.

Conclusion

Although we confirmed that SBRs fairly discriminated PD from healthy older adults without any correction, the correction methods improved their discrimination performance in a multisite setting. Our results support the utility of harmonization methods with ComBat for consolidating SBR-based diagnosis or stratification of PD in multisite studies. Nonetheless, given the substantial changes in the mean values of ComBat-corrected SBRs, caution is advised when interpreting them.

Abnormal functional connectivity of the reward network is associated with social communication impairments in autism spectrum disorder: A large-scale multi-site resting-state fMRI study

BACKGROUND

The social motivation hypothesis proposes that the social deficits of autism spectrum disorder (ASD) are related to reward system dysfunction. However, functional connectivity (FC) patterns of the reward network in ASD have not been systematically explored yet.

METHODS

The reward network was defined as eight regions of interest (ROIs) per hemisphere, including the nucleus accumbens (NAc), caudate, putamen, anterior cingulate cortex (ACC), ventromedial prefrontal cortex (vmPFC), orbitofrontal cortex (OFC), amygdala, and insula. We computed both the ROI-wise resting-state FC and seed-based whole-brain FC in 298 ASD participants and 348 typically developing (TD) controls from the Autism Brain Imaging Data Exchange I dataset. Two-sample t-tests were applied to obtain the aberrant FCs. Then, the association between aberrant FCs and clinical symptoms was assessed with Pearson's correlation or Spearman's correlation. In addition, Neurosynth Image Decoder was used to generate word clouds verifying the cognitive functions of the aberrant pathways. Furthermore, a three-way multivariate analysis of variance (MANOVA) was conducted to examine the effects of gender, subtype and age on the atypical FCs.

RESULTS

For the within network analysis, the left ACC showed weaker FCs with both the right amygdala and left NAc in ASD compared with TD, which were negatively correlated with the Autism Diagnostic Observation Schedule (ADOS) total scores and Social Responsiveness Scale (SRS) total scores respectively. For the whole-brain analysis, weaker FC (i.e., FC between the left vmPFC and left calcarine gyrus, and between the right vmPFC and left precuneus) accompanied by stronger FC (i.e., FC between the left caudate and right insula) were exhibited in ASD relative to TD, which were positively associated with the SRS motivation scores. Additionally, we detected the main effect of age on FC between the left vmPFC and left calcarine gyrus, of subtype on FC between the right vmPFC and left precuneus, of age and age-by-gender interaction on FC between the left caudate and right insula.

CONCLUSIONS

Our findings highlight the crucial role of abnormal FC patterns of the reward network in the core social deficits of ASD, which have the potential to reveal new biomarkers for ASD.

Using brain structural neuroimaging measures to predict psychosis onset for individuals at clinical high-risk

Machine learning approaches using structural magnetic resonance imaging (sMRI) can be informative for disease classification, although their ability to predict psychosis is largely unknown. We created a model with individuals at CHR who developed psychosis later (CHR-PS+) from healthy controls (HCs) that can differentiate each other. We also evaluated whether we could distinguish CHR-PS+ individuals from those who did not develop psychosis later (CHR-PS-) and those with uncertain follow-up status (CHR-UNK). T1-weighted structural brain MRI scans from 1165 individuals at CHR (CHR-PS+, n = 144; CHR-PS-, n = 793; and CHR-UNK, n = 228), and 1029 HCs, were obtained from 21 sites. We used ComBat to harmonize measures of subcortical volume, cortical thickness and surface area data and corrected for non-linear effects of age and sex using a general additive model. CHR-PS+ (n = 120) and HC (n = 799) data from 20 sites served as a training dataset, which we used to build a classifier. The remaining samples were used external validation datasets to evaluate classifier performance (test, independent confirmatory, and independent group [CHR-PS- and CHR-UNK] datasets). The accuracy of the classifier on the training and independent confirmatory datasets was 85% and 73% respectively. Regional cortical surface area measures-including those from the right superior frontal, right superior temporal, and bilateral insular cortices strongly contributed to classifying CHR-PS+ from HC. CHR-PS- and CHR-UNK individuals were more likely to be classified as HC compared to CHR-PS+ (classification rate to HC: CHR-PS+, 30%; CHR-PS-, 73%; CHR-UNK, 80%). We used multisite sMRI to train a classifier to predict psychosis onset in CHR individuals, and it showed promise predicting CHR-PS+ in an independent sample. The results suggest that when considering adolescent brain development, baseline MRI scans for CHR individuals may be helpful to identify their prognosis. Future prospective studies are required about whether the classifier could be actually helpful in the clinical settings.

Impaired segregation of the attention deficit hyperactivity disorder related pattern in children

BACKGROUND

Inattention is a key characteristic of attention deficit hyperactivity disorder (ADHD). Specific brain abnormalities associated with this symptom form a discernible pattern related with ADHD in children (i.e., ADHD related pattern) in our earlier research. The developmental processes of segregation and integration may be crucial to ADHD. However, how brains reconfigure these processes of the ADHD related pattern in different subtypes of ADHD and across sexes remain unclear.

METHODS

Nested-spectral partition method was applied to identify effects of subtype and sex on segregation and integration of the ADHD related pattern, using 145 ADHD patients and 135 typically developing controls (TDC) aged 7-14. Relationships between the measures and inattention symptoms were also investigated.

RESULTS

Children with ADHD exhibited lower segregation of the ADHD related pattern (p = 1.17 × 10-8) than TDCs. Only the main effect of subtype was significant (p = 1.14 × 10-5). Both ADHD-C (p = 2.16 × 10-6) and ADHD-I (p = 2.87 × 10-6) patients had lower segregation components relative to the TDC. Moreover, segregation components were negatively correlated with inattention scores.

CONCLUSIONS

This study identified impaired segregation in the ADHD related pattern of children with ADHD and found shared neural bases among different subtypes and sexes.

Neuroimaging Biomarkers for Drug Discovery and Development in Schizophrenia

Schizophrenia is a chronic mental illness that affects up to 1% of the population. While efficacious therapies are available for positive symptoms, effective treatment of cognitive and negative symptoms remains an unmet need after decades of research. New developments in the field of neuroimaging are accelerating our knowledge gain regarding the underlying pathophysiology of symptoms in schizophrenia and psychosis spectrum disorders, inspiring new targets for drug development. However, no validated and qualified biomarkers are currently available to support the development of new therapeutics. This review summarizes the current use of neuroimaging technology in clinical drug development for psychotic disorders. As exemplified by drug development programs that target NMDA receptor hypofunction, neuroimaging results play a critical role in target discovery and establishing target engagement and dose selection. Furthermore, pharmacological neuroimaging may provide response biomarkers that allow for early decision making in proof-of-concept studies that leverage pharmacological challenge models in healthy volunteers. That said, while response and predictive biomarkers are starting to be evaluated in patient populations, they continue to play a limited role. Novel approaches to neuroimaging data acquisition and analysis may aid the establishment of biomarkers that are predictive at the individual level in the future. Nevertheless, various gaps in knowledge need to be addressed and biomarkers need to be validated to establish them as "fit for purpose" in drug development.

Efficacy of MRI data harmonization in the age of machine learning: a multicenter study across 36 datasets

Pooling publicly-available MRI data from multiple sites allows to assemble extensive groups of subjects, increase statistical power, and promote data reuse with machine learning techniques. The harmonization of multicenter data is necessary to reduce the confounding effect associated with non-biological sources of variability in the data. However, when applied to the entire dataset before machine learning, the harmonization leads to data leakage, because information outside the training set may affect model building, and potentially falsely overestimate performance. We propose a 1) measurement of the efficacy of data harmonization; 2) harmonizer transformer, i.e., an implementation of the ComBat harmonization allowing its encapsulation among the preprocessing steps of a machine learning pipeline, avoiding data leakage by design. We tested these tools using brain T1-weighted MRI data from 1740 healthy subjects acquired at 36 sites. After harmonization, the site effect was removed or reduced, and we showed the data leakage effect in predicting individual age from MRI data, highlighting that introducing the harmonizer transformer into a machine learning pipeline allows for avoiding data leakage by design.

Multi-feature concatenation and multi-classifier stacking: An interpretable and generalizable machine learning method for MDD discrimination with rsfMRI

Major depressive disorder (MDD) is a serious and heterogeneous psychiatric disorder that needs accurate diagnosis. Resting-state functional MRI (rsfMRI), which captures multiple perspectives on brain structure, function, and connectivity, is increasingly applied in the diagnosis and pathological research of MDD. Different machine learning algorithms are then developed to exploit the rich information in rsfMRI and discriminate MDD patients from normal controls. Despite recent advances reported, the MDD discrimination accuracy has room for further improvement. The generalizability and interpretability of the discrimination method are not sufficiently addressed either. Here, we propose a machine learning method (MFMC) for MDD discrimination by concatenating multiple features and stacking multiple classifiers. MFMC is tested on the REST-meta-MDD data set that contains 2428 subjects collected from 25 different sites. MFMC yields 96.9% MDD discrimination accuracy, demonstrating a significant improvement over existing methods. In addition, the generalizability of MFMC is validated by the good performance when the training and testing subjects are from independent sites. The use of XGBoost as the meta classifier allows us to probe the decision process of MFMC. We identify 13 feature values related to 9 brain regions including the posterior cingulate gyrus, superior frontal gyrus orbital part, and angular gyrus, which contribute most to the classification and also demonstrate significant differences at the group level. The use of these 13 feature values alone can reach 87% of MFMC's full performance when taking all feature values. These features may serve as clinically useful diagnostic and prognostic biomarkers for MDD in the future.

Traveling Subject-Informed Harmonization Increases Reliability of Brain Diffusion Tensor and Neurite Mapping

Diffusion-weighted magnetic resonance imaging (dMRI) of brain has helped elucidate the microstructural changes of psychiatric and neurodegenerative disorders. Inconsistency between MRI models has hampered clinical application of dMRI-based metrics. Using harmonized dMRI data of 300 scans from 69 traveling subjects (TS) scanning the same individuals at multiple conditions with 13 MRI models and 2 protocols, the widely-used metrics such as diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) were evaluated before and after harmonization with a combined association test (ComBat) or TS-based general linear model (TS-GLM). Results showed that both ComBat and TS-GLM significantly reduced the effects of the MRI site, model, and protocol for diffusion metrics while maintaining the intersubject biological effects. The harmonization power of TS-GLM based on TS data model is more powerful than that of ComBat. In conclusion, our research demonstrated that although ComBat and TS-GLM harmonization approaches were effective at reducing the scanner effects of the site, model, and protocol for DTI and NODDI metrics in WM, they exhibited high retainability of biological effects. Therefore, we suggest that, after harmonizing DTI and NODDI metrics, a multisite study with large cohorts can accurately detect small pathological changes by retaining pathological effects.

Effects of MRI scanner manufacturers in classification tasks with deep learning models

Deep learning has become a leading subset of machine learning and has been successfully employed in diverse areas, ranging from natural language processing to medical image analysis. In medical imaging, researchers have progressively turned towards multi-center neuroimaging studies to address complex questions in neuroscience, leveraging larger sample sizes and aiming to enhance the accuracy of deep learning models. However, variations in image pixel/voxel characteristics can arise between centers due to factors including differences in magnetic resonance imaging scanners. Such variations create challenges, particularly inconsistent performance in machine learning-based approaches, often referred to as domain shift, where the trained models fail to achieve satisfactory or improved results when confronted with dissimilar test data. This study analyzes the performance of multiple disease classification tasks using multi-center MRI data obtained from three widely used scanner manufacturers (GE, Philips, and Siemens) across several deep learning-based networks. Furthermore, we investigate the efficacy of mitigating scanner vendor effects using ComBat-based harmonization techniques when applied to multi-center datasets of 3D structural MR images. Our experimental results reveal a substantial decline in classification performance when models trained on one type of scanner manufacturer are tested with data from different manufacturers. Moreover, despite applying ComBat-based harmonization, the harmonized images do not demonstrate any noticeable performance enhancement for disease classification tasks.

MISPEL: A supervised deep learning harmonization method for multi-scanner neuroimaging data

Large-scale data obtained from aggregation of already collected multi-site neuroimaging datasets has brought benefits such as higher statistical power, reliability, and robustness to the studies. Despite these promises from growth in sample size, substantial technical variability stemming from differences in scanner specifications exists in the aggregated data and could inadvertently bias any downstream analyses on it. Such a challenge calls for data normalization and/or harmonization frameworks, in addition to comprehensive criteria to estimate the scanner-related variability and evaluate the harmonization frameworks. In this study, we propose MISPEL (Multi-scanner Image harmonization via Structure Preserving Embedding Learning), a supervised multi-scanner harmonization method that is naturally extendable to more than two scanners. We also designed a set of criteria to investigate the scanner-related technical variability and evaluate the harmonization techniques. As an essential requirement of our criteria, we introduced a multi-scanner matched dataset of 3T T1 images across four scanners, which, to the best of our knowledge is one of the few datasets of this kind. We also investigated our evaluations using two popular segmentation frameworks: FSL and segmentation in statistical parametric mapping (SPM). Lastly, we compared MISPEL to popular methods of normalization and harmonization, namely White Stripe, RAVEL, and CALAMITI. MISPEL outperformed these methods and is promising for many other neuroimaging modalities.

Removal of site effects and enhancement of signal using dual projection independent component analysis for pooling multi-site MRI data

Combining magnetic resonance imaging (MRI) data from multi-site studies is a popular approach for constructing larger datasets to greatly enhance the reliability and reproducibility of neuroscience research. However, the scanner/site variability is a significant confound that complicates the interpretation of the results, so effective and complete removal of the scanner/site variability is necessary to realise the full advantages of pooling multi-site datasets. Independent component analysis (ICA) and general linear model (GLM) based harmonisation methods are the two primary methods used to eliminate scanner/site effects. Unfortunately, there are challenges with both ICA-based and GLM-based harmonisation methods to remove site effects completely when the signals of interest and scanner/site effects-related variables are correlated, which may occur in neuroscience studies. In this study, we propose an effective and powerful harmonisation strategy that implements dual projection (DP) theory based on ICA to remove the scanner/site effects more completely. This method can separate the signal effects correlated with site variables from the identified site effects for removal without losing signals of interest. Both simulations and vivo structural MRI datasets, including a dataset from Autism Brain Imaging Data Exchange II and a travelling subject dataset from the Strategic Research Program for Brain Sciences, were used to test the performance of a DP-based ICA harmonisation method. Results show that DP-based ICA harmonisation has superior performance for removing site effects and enhancing the sensitivity to detect signals of interest as compared with GLM-based and conventional ICA harmonisation methods.

Harmonization of multi-site functional MRI data with dual-projection based ICA model

Modern neuroimaging studies frequently merge magnetic resonance imaging (MRI) data from multiple sites. A larger and more diverse group of participants can increase the statistical power, enhance the reliability and reproducibility of neuroimaging research, and obtain findings more representative of the general population. However, measurement biases caused by site differences in scanners represent a barrier when pooling data collected from different sites. The existence of site effects can mask biological effects and lead to spurious findings. We recently proposed a powerful denoising strategy that implements dual-projection (DP) theory based on ICA to remove site-related effects from pooled data, demonstrating the method for simulated and in vivo structural MRI data. This study investigates the use of our DP-based ICA denoising method for harmonizing functional MRI (fMRI) data collected from the Autism Brain Imaging Data Exchange II. After frequency-domain and regional homogeneity analyses, two modalities, including amplitude of low frequency fluctuation (ALFF) and regional homogeneity (ReHo), were used to validate our method. The results indicate that DP-based ICA denoising method removes unwanted site effects for both two fMRI modalities, with increases in the significance of the associations between non-imaging variables (age, sex, etc.) and fMRI measures. In conclusion, our DP method can be applied to fMRI data in multi-site studies, enabling more accurate and reliable neuroimaging research findings.

ComBat Harmonization: Empirical Bayes versus fully Bayes approaches

Studying small effects or subtle neuroanatomical variation requires large-scale sample size data. As a result, combining neuroimaging data from multiple datasets is necessary. Variation in acquisition protocols, magnetic field strength, scanner build, and many other non-biologically related factors can introduce undesirable bias into studies. Hence, harmonization is required to remove the bias-inducing factors from the data. ComBat is one of the most common methods applied to features from structural images. ComBat models the data using a hierarchical Bayesian model and uses the empirical Bayes approach to infer the distribution of the unknown factors. The empirical Bayes harmonization method is computationally efficient and provides valid point estimates. However, it tends to underestimate uncertainty. This paper investigates a new approach, fully Bayesian ComBat, where Monte Carlo sampling is used for statistical inference. When comparing fully Bayesian and empirical Bayesian ComBat, we found Empirical Bayesian ComBat more effectively removed scanner strength information and was much more computationally efficient. Conversely, fully Bayesian ComBat better preserved biological disease and age-related information while performing more accurate harmonization on traveling subjects. The fully Bayesian approach generates a rich posterior distribution, which is useful for generating simulated imaging features for improving classifier performance in a limited data setting. We show the generative capacity of our model for augmenting and improving the detection of patients with Alzheimer's disease. Posterior distributions for harmonized imaging measures can also be used for brain-wide uncertainty comparison and more principled downstream statistical analysis.Code for our new fully Bayesian ComBat extension is available at https://github.com/batmanlab/BayesComBat.

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.

Child Developmental MRI (CDM) project: protocol for a multi-centre, cross-sectional study on elucidating the pathophysiology of attention-deficit/hyperactivity disorder and autism spectrum disorder through a multi-dimensional approach

INTRODUCTION

Neuroimaging studies on attention-deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) have demonstrated differences in extensive brain structure, activity and network. However, there remains heterogeneity and inconsistency across these findings, presumably because of the diversity of the disorders themselves, small sample sizes, and site and parameter differences in MRI scanners, and their overall pathogenesis remains unclear. To address these gaps in the literature, we will apply the travelling-subject approach to correct site differences in MRI scanners and clarify brain structure and network characteristics of children with ADHD and ASD using large samples collected in a multi-centre collaboration. In addition, we will investigate the relationship between these characteristics and genetic, epigenetic, biochemical markers, and behavioural and psychological measures.

METHODS AND ANALYSIS

We will collect resting-state functional MRI (fMRI) and T1-weighted and diffusion-weighted MRI data from 15 healthy adults as travelling subjects and 300 children (ADHD, n=100; ASD, n=100; and typical development, n=100) with multi-dimensional assessments. We will also apply data from more than 1000 samples acquired in our previous neuroimaging studies on ADHD and ASD.

ETHICS AND DISSEMINATION

The study protocol has been approved by the Research Ethics Committee of the University of Fukui Hospital (approval no: 20220601). Our study findings will be submitted to scientific peer-reviewed journals and conferences.

DeepComBat: A Statistically Motivated, Hyperparameter-Robust, Deep Learning Approach to Harmonization of Neuroimaging Data

Neuroimaging data from multiple batches (i.e. acquisition sites, scanner manufacturer, datasets, etc.) are increasingly necessary to gain new insights into the human brain. However, multi-batch data, as well as extracted radiomic features, exhibit pronounced technical artifacts across batches. These batch effects introduce confounding into the data and can obscure biological effects of interest, decreasing the generalizability and reproducibility of findings. This is especially true when multi-batch data is used alongside complex downstream analysis models, such as machine learning methods. Image harmonization methods seeking to remove these batch effects are important for mitigating these issues; however, significant multivariate batch effects remain in the data following harmonization by current state-of-the-art statistical and deep learning methods. We present DeepCombat, a deep learning harmonization method based on a conditional variational autoencoder architecture and the ComBat harmonization model. DeepCombat learns and removes subject-level batch effects by accounting for the multivariate relationships between features. Additionally, DeepComBat relaxes a number of strong assumptions commonly made by previous deep learning harmonization methods and is empirically robust across a wide range of hyperparameter choices. We apply this method to neuroimaging data from a large cognitive-aging cohort and find that DeepCombat outperforms existing methods, as assessed by a battery of machine learning methods, in removing scanner effects from cortical thickness measurements while preserving biological heterogeneity. Additionally, DeepComBat provides a new perspective for statistically-motivated deep learning harmonization methods.

Machine Learning for Brain MRI Data Harmonisation: A Systematic Review

BACKGROUND

Magnetic Resonance Imaging (MRI) data collected from multiple centres can be heterogeneous due to factors such as the scanner used and the site location. To reduce this heterogeneity, the data needs to be harmonised. In recent years, machine learning (ML) has been used to solve different types of problems related to MRI data, showing great promise.

OBJECTIVE

This study explores how well various ML algorithms perform in harmonising MRI data, both implicitly and explicitly, by summarising the findings in relevant peer-reviewed articles. Furthermore, it provides guidelines for the use of current methods and identifies potential future research directions.

METHOD

This review covers articles published through PubMed, Web of Science, and IEEE databases through June 2022. Data from studies were analysed based on the criteria of Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). Quality assessment questions were derived to assess the quality of the included publications.

RESULTS

a total of 41 articles published between 2015 and 2022 were identified and analysed. In the review, MRI data has been found to be harmonised either in an implicit (n = 21) or an explicit (n = 20) way. Three MRI modalities were identified: structural MRI (n = 28), diffusion MRI (n = 7) and functional MRI (n = 6).

CONCLUSION

Various ML techniques have been employed to harmonise different types of MRI data. There is currently a lack of consistent evaluation methods and metrics used across studies, and it is recommended that the issue be addressed in future studies. Harmonisation of MRI data using ML shows promises in improving performance for ML downstream tasks, while caution should be exercised when using ML-harmonised data for direct interpretation.

Disrupted network integration and segregation involving the default mode network in autism spectrum disorder

Changes in the brain's default mode network (DMN) in the resting state are closely related to autism spectrum disorder (ASD). Module segmentation can effectively elucidate the neural mechanism of ASD and explore intra- and inter-network connections by means of the participation coefficient (PC). We used that resting-state fMRI data from 269 ASD patients and 340 healthy controls (HCs) in the current study. From the results, ASD subjects showed a significantly higher PC of the DMN than HC subjects. This difference was related to lower intra-module connections within the DMN and higher inter-network connections between the DMN and other networks. When the subjects were split into age groups, the results were verified in the 7-12- and 12-18-year-old age groups but not in the young adult group (18-25 years). When the subjects were divided according to different subtypes of ASD, the results were also observed in the classic autism and pervasive developmental disorder groups, but not in the Asperger disorder group. In conclusions, less developed network segregation in the DMN could be a valid biomarker for ASD. This provides network scientists with new insights into the intermodular connectivity configurations of complex networks from different dimensions in a systematic and comprehensive manner.

Reproducibility of cerebellar involvement as quantified by consensus structural MRI biomarkers in advanced essential tremor

Essential tremor (ET) is the most prevalent movement disorder with poorly understood etiology. Some neuroimaging studies report cerebellar involvement whereas others do not. This discrepancy may stem from underpowered studies, differences in statistical modeling or variation in magnetic resonance imaging (MRI) acquisition and processing. To resolve this, we investigated the cerebellar structural differences using a local advanced ET dataset augmented by matched controls from PPMI and ADNI. We tested the hypothesis of cerebellar involvement using three neuroimaging biomarkers: VBM, gray/white matter volumetry and lobular volumetry. Furthermore, we assessed the impacts of statistical models and segmentation pipelines on results. Results indicate that the detected cerebellar structural changes vary with methodology. Significant reduction of right cerebellar gray matter and increase of the left cerebellar white matter were the only two biomarkers consistently identified by multiple methods. Results also show substantial volumetric overestimation from SUIT-based segmentation-partially explaining previous literature discrepancies. This study suggests that current estimation of cerebellar involvement in ET may be overemphasized in MRI studies and highlights the importance of methods sensitivity analysis on results interpretation. ET datasets with large sample size and replication studies are required to improve our understanding of regional specificity of cerebellum involvement in ET. PROTOCOL REGISTRATION: The stage 1 protocol for this Registered Report was accepted in principle on 21 March 2022. The protocol, as accepted by the journal, can be found at: https://doi.org/10.6084/m9.figshare.19697776 .

Comparison of different approaches to manage multi-site magnetic resonance spectroscopy clinical data analysis

Introduction

The effects caused by differences in data acquisition can be substantial and may impact data interpretation in multi-site/scanner studies using magnetic resonance spectroscopy (MRS). Given the increasing use of multi-site studies, a better understanding of how to account for different scanners is needed. Using data from a concussion population, we compare ComBat harmonization with different statistical methods in controlling for site, vendor, and scanner as covariates to determine how to best control for multi-site data.

Methods

The data for the current study included 545 MRS datasets to measure tNAA, tCr, tCho, Glx, and mI to study the pediatric concussion acquired across five sites, six scanners, and two different MRI vendors. For each metabolite, the site and vendor were accounted for in seven different models of general linear models (GLM) or mixed-effects models while testing for group differences between the concussion and orthopedic injury. Models 1 and 2 controlled for vendor and site. Models 3 and 4 controlled for scanner. Models 5 and 6 controlled for site applied to data harmonized by vendor using ComBat. Model 7 controlled for scanner applied to data harmonized by scanner using ComBat. All the models controlled for age and sex as covariates.

Results

Models 1 and 2, controlling for site and vendor, showed no significant group effect in any metabolites, but the vendor and site were significant factors in the GLM. Model 3, which included a scanner, showed a significant group effect for tNAA and tCho, and the scanner was a significant factor. Model 4, controlling for the scanner, did not show a group effect in the mixed model. The data harmonized by the vendor using ComBat (Models 5 and 6) had no significant group effect in both the GLM and mixed models. Lastly, the data harmonized by the scanner using ComBat (Model 7) showed no significant group effect. The individual site data suggest there were no group differences.

Conclusion

Using data from a large clinical concussion population, different analysis techniques to control for site, vendor, and scanner in MRS data yielded different results. The findings support the use of ComBat harmonization for clinical MRS data, as it removes the site and vendor effects.

Removing the effects of the site in brain imaging machine-learning - Measurement and extendable benchmark

Multisite machine-learning neuroimaging studies, such as those conducted by the ENIGMA Consortium, need to remove the differences between sites to avoid effects of the site (EoS) that may prevent or fraudulently help the creation of prediction models, leading to impoverished or inflated prediction accuracy. Unfortunately, we have shown earlier that current Methods Aiming to Remove the EoS (MAREoS, e.g., ComBat) cannot remove complex EoS (e.g., including interactions between regions). And complex EoS may bias the accuracy. To overcome this hurdle, groups worldwide are developing novel MAREoS. However, we cannot assess their effectiveness because EoS may either inflate or shrink the accuracy, and MAREoS may both remove the EoS and degrade the data. In this work, we propose a strategy to measure the effectiveness of a MAREoS in removing different types of EoS. FOR MAREOS DEVELOPERS, we provide two multisite MRI datasets with only simple true effects (i.e., detectable by most machine-learning algorithms) and two with only simple EoS (i.e., removable by most MAREoS). First, they should use these datasets to fit machine-learning algorithms after applying the MAREoS. Second, they should use the formulas we provide to calculate the relative accuracy change associated with the MAREoS in each dataset and derive an EoS-removal effectiveness statistic. We also offer similar datasets and formulas for complex true effects and EoS that include first-order interactions. FOR MACHINE-LEARNING RESEARCHERS, we provide an extendable benchmark website to show: a) the types of EoS they should remove for each given machine-learning algorithm and b) the effectiveness of each MAREoS for removing each type of EoS. Relevantly, a MAREoS only able to remove the simple EoS may suffice for simple machine-learning algorithms, whereas more complex algorithms need a MAREoS that can remove more complex EoS. For instance, ComBat removes all simple EoS as needed for predictions based on simple lasso algorithms, but it leaves residual complex EoS that may bias the predictions based on standard support vector machine algorithms.

Structural brain abnormalities in schizophrenia patients with a history and presence of auditory verbal hallucination

Although many studies have demonstrated structural brain abnormalities associated with auditory verbal hallucinations (AVH) in schizophrenia, the results remain inconsistent because of the small sample sizes and the reliability of clinical interviews. We compared brain morphometries in 204 participants, including 58 schizophrenia patients with a history of AVH (AVH + ), 29 without a history of AVH (AVH-), and 117 healthy controls (HCs) based on a detailed inspection of medical records. We further divided the AVH+ group into 37 patients with and 21 patients without hallucinations at the time of the MRI scans (AVH++ and AVH+-, respectively) via clinical interviews to explore the morphological differences according to the persistence of AVH. The AVH + group had a smaller surface area in the left caudal middle frontal gyrus (F = 7.28, FDR-corrected p = 0.0008) and precentral gyrus (F = 7.68, FDR-corrected p = 0.0006) compared to the AVH- group. The AVH+ patients had a smaller surface area in the left insula (F = 7.06, FDR-corrected p = 0.001) and a smaller subcortical volume in the bilateral hippocampus (right: F = 13.34, FDR-corrected p = 0.00003; left: F = 6.80, FDR-corrected p = 0.001) compared to the HC group. Of these significantly altered areas, the AVH++ group showed significantly smaller bilateral hippocampal volumes compared to the AVH+- group, and a smaller surface area in the left precentral gyrus and caudal middle frontal gyrus compared to the AVH- group. Our findings highlighted the distinct pattern of structural alteration between the history and presence of AVH in schizophrenia, and the importance of integrating multiple criteria to elucidate the neuroanatomical mechanisms.

Validation of cross-sectional and longitudinal ComBat harmonization methods for magnetic resonance imaging data on a travelling subject cohort

Background

The growth in multi-center neuroimaging studies generated a need for methods that mitigate the differences in hardware and acquisition protocols across sites i.e., scanner effects. ComBat harmonization methods have shown promise but have not yet been tested on all the data types commonly studied with magnetic resonance imaging (MRI). This study aimed to validate neuroCombat, longCombat and gamCombat on both structural and diffusion metrics in both cross-sectional and longitudinal data.

Methods

We used a travelling subject design whereby 73 healthy volunteers contributed 161 scans across two sites and four machines using one T1 and five diffusion MRI protocols. Scanner was defined as a composite of site, machine and protocol. A common pipeline extracted two structural metrics (volumes and cortical thickness) and two diffusion tensor imaging metrics (mean diffusivity and fractional anisotropy) for seven regions of interest including gray and (except for cortical thickness) white matter regions.

Results

Structural data exhibited no significant scanner effect and therefore did not benefit from harmonization in our particular cohort. Indeed, attempting harmonization obscured the true biological effect for some regions of interest. Diffusion data contained marked scanner effects and was successfully harmonized by all methods, resulting in smaller scanner effects and better detection of true biological effects. LongCombat less effectively reduced the scanner effect for cross-sectional white matter data but had a slightly lower probability of incorrectly finding group differences in simulations, compared to neuroCombat and gamCombat. False positive rates for all methods and all metrics did not significantly exceed 5%.

Conclusions

Statistical harmonization of structural data is not always necessary and harmonization in the absence of a scanner effect may be harmful. Harmonization of diffusion MRI data is highly recommended with neuroCombat, longCombat and gamCombat performing well in cross-sectional and longitudinal settings.

Science Forum: How failure to falsify in high-volume science contributes to the replication crisis

The number of scientific papers published every year continues to increase, but scientific knowledge is not progressing at the same rate. Here we argue that a greater emphasis on falsification – the direct testing of strong hypotheses – would lead to faster progress by allowing well-specified hypotheses to be eliminated. We describe an example from neuroscience where there has been little work to directly test two prominent but incompatible hypotheses related to traumatic brain injury. Based on this example, we discuss how building strong hypotheses and then setting out to falsify them can bring greater precision to the clinical neurosciences, and argue that this approach could be beneficial to all areas of science.

Harmonization of multi-site MRS data with ComBat

Magnetic resonance spectroscopy (MRS) is a non-invasive neuroimaging technique used to measure brain chemistry in vivo and has been used to study the healthy brain as well as neuropathology in numerous neurological disorders. The number of multi-site studies using MRS are increasing; however, non-biological variability introduced during data collection across multiple sites, such as differences in scanner vendors and site-specific acquisition implementations for MRS, can obscure detection of biological effects of interest. ComBat is a data harmonization technique that can remove non-biological sources of variance in multisite studies. It has been validated for use with structural and functional MRI metrics but not for MRS metabolites. This study investigated the validity of using ComBat to harmonize MRS metabolites for vendor and site differences. Analyses were performed using data acquired across 20 sites and included edited MRS for GABA+ (N=218) and macromolecule-suppressed GABA data (N=209), as well as standard PRESS data to quantify NAA, creatine, choline, and glutamate (N=190). ComBat harmonization successfully mitigated vendor and site differences for all metabolites of interest. Moreover, significant associations were detected between sex and choline levels and between age and glutamate and GABA+ levels that were not detectable prior to harmonization, confirming the importance of removing site and vendor effects in multi-site data. In conclusion, ComBat harmonization can be successfully applied to MRS data in multi-site MRS studies.

A deep learning-based multisite neuroimage harmonization framework established with a traveling-subject dataset

The accumulation of multisite large-sample MRI datasets collected during large brain research projects in the last decade has provided critical resources for understanding the neurobiological mechanisms underlying cognitive functions and brain disorders. However, the significant site effects observed in imaging data and their derived structural and functional features have prevented the derivation of consistent findings across multiple studies. The development of harmonization methods that can effectively eliminate complex site effects while maintaining biological characteristics in neuroimaging data has become a vital and urgent requirement for multisite imaging studies. Here, we propose a deep learning-based framework to harmonize imaging data obtained from pairs of sites, in which site factors and brain features can be disentangled and encoded. We trained the proposed framework with a publicly available traveling subject dataset from the Strategic Research Program for Brain Sciences (SRPBS) and harmonized the gray matter volume maps derived from eight source sites to a target site. The proposed framework significantly eliminated intersite differences in gray matter volumes. The embedded encoders successfully captured both the abstract textures of site factors and the concrete brain features. Moreover, the proposed framework exhibited outstanding performance relative to conventional statistical harmonization methods in terms of site effect removal, data distribution homogenization, and intrasubject similarity improvement. Finally, the proposed harmonization network provided fixable expandability, through which new sites could be linked to the target site via indirect schema without retraining the whole model. Together, the proposed method offers a powerful and interpretable deep learning-based harmonization framework for multisite neuroimaging data that can enhance reliability and reproducibility in multisite studies regarding brain development and brain disorders.

Lifespan Volume Trajectories From Non–harmonized T1–Weighted MRI Do Not Differ After Site Correction Based on Traveling Human Phantoms

Multi–site imaging consortiums strive to increase participant numbers by pooling data across sites, but scanner related differences can bias results. This study combines data from three research MRI centers, including three different scanner models from two vendors, to examine non–harmonized T1–weighted brain imaging protocols in two cohorts. First, 23 human traveling phantoms were scanned twice each at all three sites (six scans per person; 138 scans total) to quantify within–participant variability of brain volumes (total brain, white matter, gray matter, lateral ventricles, thalamus, caudate, putamen and globus pallidus), and to calculate site–specific correction factors for each structure. Sample size calculations were used to determine the number of traveling phantoms needed to achieve effect sizes for observed differences to help guide future studies. Next, cross–sectional lifespan volume trajectories were examined in 856 healthy participants (5—91 years of age) scanned at these sites. Cross–sectional trajectories of volume versus age for each structure were then compared before and after application of traveling phantom based site–specific correction factors, as well as correction using the open–source method ComBat. Although small systematic differences between sites were observed in the traveling phantom analysis, correction for site using either method had little impact on the lifespan trajectories. Only white matter had small but significant differences in the intercept parameter after ComBat correction (but not traveling phantom based correction), while no other fits differed. This suggests that age–related changes over the lifespan outweigh systematic differences between scanners for volumetric analysis. This work will help guide pooling of multisite datasets as well as meta–analyses of data from non–harmonized protocols.

Quantitative MRI Harmonization to Maximize Clinical Impact: The RIN-Neuroimaging Network

Neuroimaging studies often lack reproducibility, one of the cardinal features of the scientific method. Multisite collaboration initiatives increase sample size and limit methodological flexibility, therefore providing the foundation for increased statistical power and generalizable results. However, multisite collaborative initiatives are inherently limited by hardware, software, and pulse and sequence design heterogeneities of both clinical and preclinical MRI scanners and the lack of benchmark for acquisition protocols, data analysis, and data sharing. We present the overarching vision that yielded to the constitution of RIN-Neuroimaging Network, a national consortium dedicated to identifying disease and subject-specific in-vivo neuroimaging biomarkers of diverse neurological and neuropsychiatric conditions. This ambitious goal needs efforts toward increasing the diagnostic and prognostic power of advanced MRI data. To this aim, 23 Italian Scientific Institutes of Hospitalization and Care (IRCCS), with technological and clinical specialization in the neurological and neuroimaging field, have gathered together. Each IRCCS is equipped with high- or ultra-high field MRI scanners (i.e., ≥3T) for clinical or preclinical research or has established expertise in MRI data analysis and infrastructure. The actions of this Network were defined across several work packages (WP). A clinical work package (WP1) defined the guidelines for a minimum standard clinical qualitative MRI assessment for the main neurological diseases. Two neuroimaging technical work packages (WP2 and WP3, for clinical and preclinical scanners) established Standard Operative Procedures for quality controls on phantoms as well as advanced harmonized quantitative MRI protocols for studying the brain of healthy human participants and wild type mice. Under FAIR principles, a web-based e-infrastructure to store and share data across sites was also implemented (WP4). Finally, the RIN translated all these efforts into a large-scale multimodal data collection in patients and animal models with dementia (i.e., case study). The RIN-Neuroimaging Network can maximize the impact of public investments in research and clinical practice acquiring data across institutes and pathologies with high-quality and highly-consistent acquisition protocols, optimizing the analysis pipeline and data sharing procedures.

Quantitative Structural Brain Magnetic Resonance Imaging Analyses: Methodological Overview and Application to Rett Syndrome

Congenital genetic disorders often present with neurological manifestations such as neurodevelopmental disorders, motor developmental retardation, epilepsy, and involuntary movement. Through qualitative morphometric evaluation of neuroimaging studies, remarkable structural abnormalities, such as lissencephaly, polymicrogyria, white matter lesions, and cortical tubers, have been identified in these disorders, while no structural abnormalities were identified in clinical settings in a large population. Recent advances in data analysis programs have led to significant progress in the quantitative analysis of anatomical structural magnetic resonance imaging (MRI) and diffusion-weighted MRI tractography, and these approaches have been used to investigate psychological and congenital genetic disorders. Evaluation of morphometric brain characteristics may contribute to the identification of neuroimaging biomarkers for early diagnosis and response evaluation in patients with congenital genetic diseases. This mini-review focuses on the methodologies and attempts employed to study Rett syndrome using quantitative structural brain MRI analyses, including voxel- and surface-based morphometry and diffusion-weighted MRI tractography. The mini-review aims to deepen our understanding of how neuroimaging studies are used to examine congenital genetic disorders.

Application of a Machine Learning Algorithm for Structural Brain Images in Chronic Schizophrenia to Earlier Clinical Stages of Psychosis and Autism Spectrum Disorder: A Multiprotocol Imaging Dataset Study

BACKGROUND AND HYPOTHESIS

Machine learning approaches using structural magnetic resonance imaging (MRI) can be informative for disease classification; however, their applicability to earlier clinical stages of psychosis and other disease spectra is unknown. We evaluated whether a model differentiating patients with chronic schizophrenia (ChSZ) from healthy controls (HCs) could be applied to earlier clinical stages such as first-episode psychosis (FEP), ultra-high risk for psychosis (UHR), and autism spectrum disorders (ASDs).

STUDY DESIGN

Total 359 T1-weighted MRI scans, including 154 individuals with schizophrenia spectrum (UHR, n = 37; FEP, n = 24; and ChSZ, n = 93), 64 with ASD, and 141 HCs, were obtained using three acquisition protocols. Of these, data regarding ChSZ (n = 75) and HC (n = 101) from two protocols were used to build a classifier (training dataset). The remainder was used to evaluate the classifier (test, independent confirmatory, and independent group datasets). Scanner and protocol effects were diminished using ComBat.

STUDY RESULTS

The accuracy of the classifier for the test and independent confirmatory datasets were 75% and 76%, respectively. The bilateral pallidum and inferior frontal gyrus pars triangularis strongly contributed to classifying ChSZ. Schizophrenia spectrum individuals were more likely to be classified as ChSZ compared to ASD (classification rate to ChSZ: UHR, 41%; FEP, 54%; ChSZ, 70%; ASD, 19%; HC, 21%).

CONCLUSION

We built a classifier from multiple protocol structural brain images applicable to independent samples from different clinical stages and spectra. The predictive information of the classifier could be useful for applying neuroimaging techniques to clinical differential diagnosis and predicting disease onset earlier.

Gray Matter Network Associated With Attention in Children With Attention Deficit Hyperactivity Disorder

BACKGROUND

Attention deficit hyperactivity disorder (ADHD) is one of the most prevalent childhood-onset neurodevelopmental disorders; however, the underlying neural mechanisms for the inattention symptom remain elusive for children with ADHD. At present, the majority of studies have analyzed the structural MRI (sMRI) with the univariate method, which fails to demonstrate the interregional covarying relationship of gray matter (GM) volumes among brain regions. The scaled subprofile model of principal component analysis (SSM-PCA) is a multivariate method, which can detect more robust brain-behavioral phenotype association compared to the univariate analysis method. This study aims to identify the GM network associated with attention in children with ADHD by applying SSM-PCA to the sMRI.

METHODS

The sMRI of 209 children with ADHD and 209 typically developing controls (TDCs) aged 7-14 years from the ADHD-200 dataset was used for anatomical computation, and the GM volume in each brain region was acquired. Then, SSM-PCA was applied to the GM volumes of all the subjects to capture the GM network of children with ADHD (i.e., ADHD-related pattern). The relationship between the expression of ADHD-related pattern and inattention symptom was further investigated. Finally, the influence of sample size on the analysis of this study was explored.

RESULTS

The ADHD-related pattern mainly included putamen, pallium, caudate, thalamus, right accumbens, superior/middle/inferior frontal cortex, superior occipital cortex, superior parietal cortex, and left middle occipital cortex. In addition, the expression of the ADHD-related pattern was related to inattention scores measured by the Conners' Parent Rating Scale long version (CPRS-LV; r = 0.25, p = 0.0004) and the DuPaul ADHD Rating Scale IV (ADHD-RS; r = 0.18, p = 0.03). Finally, we found that when the sample size was 252, the results of ADHD-related pattern were relatively reliable. Similarly, the sample size needed to be 162 when exploring the relationship between ADHD-related pattern and behavioral indicator measured by CPRS-LV.

CONCLUSION

We captured a GM network associated with attention in children with ADHD, which is different from that in adolescents and adults with ADHD. Our findings may shed light on the diverse neural mechanisms of inattention and provide treatment targets for children with ADHD.

Neuroimaging brain growth charts: A road to mental health

Mental disorders are common health concerns and contribute to a heavy global burden on our modern society. It is challenging to identify and treat them timely. Neuroimaging evidence suggests the incidence of various psychiatric and behavioral disorders is closely related to the atypical development of brain structure and function. The identification and understanding of atypical brain development provide chances for clinicians to detect mental disorders earlier, perhaps even prior to onset, and treat them more precisely. An invaluable and necessary method in identifying and monitoring atypical brain development are growth charts of typically developing individuals in the population. The brain growth charts can offer a series of standard references on typical neurodevelopment, representing an important resource for the scientific and medical communities. In the present paper, we review the relationship between mental disorders and atypical brain development from a perspective of normative brain development by surveying the recent progress in the development of brain growth charts, including four aspects on growth chart utility: 1) cohorts, 2) measures, 3) mechanisms, and 4) clinical translations. In doing so, we seek to clarify the challenges and opportunities in charting brain growth, and to promote the application of brain growth charts in clinical practice.