A comparison of methods to harmonize cortical thickness measurements across scanners and sites

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize ( https://github.com/ahzhu/eharmonize ).

Site effects in multisite fetal brain MRI: morphological insights into early brain development

OBJECTIVE

To evaluate multisite effects on fetal brain MRI. Specifically, to identify crucial acquisition factors affecting fetal brain structural measurements and developmental patterns, while assessing the effectiveness of existing harmonization methods in mitigating site effects.

MATERIALS AND METHODS

Between May 2017 and March 2022, T2-weighted fast spin-echo sequences in-utero MRI were performed on healthy fetuses from retrospectively recruited pregnant volunteers on four different scanners at four sites. A generalized additive model (GAM) was used to quantitatively assess site effects, including field strength (FS), manufacturer (M), in-plane resolution (R), and slice thickness (ST), on subcortical volume and cortical morphological measurements, including cortical thickness, curvature, and sulcal depth. Growth models were selected to elucidate the developmental trajectories of these morphological measurements. Welch's test was performed to evaluate the influence of site effects on developmental trajectories. The comBat-GAM harmonization method was applied to mitigate site-related biases.

RESULTS

The final analytic sample consisted of 340 MRI scans from 218 fetuses (mean GA, 30.1 weeks ± 4.4 [range, 21.7-40 weeks]). GAM results showed that lower FS and lower spatial resolution led to overestimations in selected brain regions of subcortical volumes and cortical morphological measurements. Only the peak cortical thickness in developmental trajectories was significantly influenced by the effects of FS and R. Notably, ComBat-GAM harmonization effectively removed site effects while preserving developmental patterns.

CONCLUSION

Our findings pinpointed the key acquisition factors in in-utero fetal brain MRI and underscored the necessity of data harmonization when pooling multisite data for fetal brain morphology investigations.

KEY POINTS

Question How do specific site MRI acquisition factors affect fetal brain imaging? Finding Lower FS and spatial resolution overestimated subcortical volumes and cortical measurements. Cortical thickness in developmental trajectories was influenced by FS and in-plane resolution. Clinical relevance This study provides important guidelines for the fetal MRI community when scanning fetal brains and underscores the necessity of data harmonization of cross-center fetal studies.

Evaluating the effects of high-throughput structural neuroimaging predictors on whole-brain functional connectome outcomes via network-based matrix-on-vector regression

The joint analysis of multimodal neuroimaging data is vital in brain research, revealing complex interactions between brain structures and functions. Our study is motivated by the analysis of a vast dataset of brain functional connectivity (FC) and multimodal structural imaging (SI) features from the UK Biobank. Specifically, we aim to investigate the effects of SI features, such as white matter microstructure integrity (WMMI) and cortical thickness, on the whole-brain functional connectome network. This analysis is inherently challenging due to the extensive structural-functional associations and the intricate network patterns present in multimodal high-dimensional neuroimaging data. To bridge methodological gaps, we developed a novel multi-level sub-graph extraction method (dense bipartite with nested unipartite graph) within a matrix(network)-on-vector regression model. This method identifies subsets of spatially specific SI features that intensely and systematically influence FC sub-networks, while effectively suppressing false positives in large-scale datasets. Applying our method to a multimodal neuroimaging dataset of 4242 participants ffrom the UK Biobank, we evaluated the effects of whole-brain WMMI and cortical thickness on resting-state FC. Our findings indicate that the WMMI in corticospinal tracts and inferior cerebellar peduncle significantly affect functional connections of sensorimotor, salience, and executive sub-networks, with an average correlation of 0.81 ($p < 0.001$).

Evaluation of ComBat Harmonization for Reducing Across-Tracer Differences in Regional Amyloid PET Analyses

Differences in amyloid positron emission tomography (PET) radiotracer pharmacokinetics and binding properties lead to discrepancies in amyloid-β uptake estimates. Harmonization of tracer-specific biases is crucial for optimal performance of downstream tasks. Here, we investigated the efficacy of ComBat, a data-driven harmonization model, for reducing tracer-specific biases in regional amyloid PET measurements from [18F]-florbetapir (FBP) and [11C]-Pittsburgh compound-B (PiB). One hundred thirteen head-to-head FBP-PiB scan pairs, scanned from the same subject within 90 days, were selected from the Open Access Series of Imaging Studies 3 (OASIS-3) dataset. The Centiloid scale, ComBat with no covariates, ComBat with biological covariates, and GAM-ComBat with biological covariates were used to harmonize both global and regional amyloid standardized uptake value ratios (SUVR). Variants of ComBat, including longitudinal ComBat and PEACE, were also tested. Intraclass correlation coefficient (ICC) and mean absolute error (MAE) were computed to measure the absolute agreement between tracers. Additionally, longitudinal amyloid SUVRs from an anti-amyloid drug trial were simulated using linear mixed effects modeling. Differences in rates-of-change between simulated treatment and placebo groups were tested, and change in statistical power/Type-I error after harmonization was quantified. In the head-to-head tracer comparison, ComBat with no covariates was the best at increasing ICC and decreasing MAE of both global summary and regional amyloid PET SUVRs between scan pairs of the same group of subjects. In the clinical trial simulation, harmonization with both Centiloid and ComBat increased statistical power of detecting true rate-of-change differences between groups and decreased false discovery rate in the absence of a treatment effect. The greatest benefit of harmonization was observed when groups exhibited differing FBP-to-PiB proportions. ComBat outperformed the Centiloid scale in harmonizing both global and regional amyloid estimates. Additionally, ComBat improved the detection of rate-of-change differences between clinical trial groups. Our findings suggest that ComBat is a viable alternative to Centiloid for harmonizing regional amyloid PET analyses.

Current best practices and future opportunities for reproducible findings using large-scale neuroimaging in psychiatry

Research into the brain basis of psychopathology is challenging due to the heterogeneity of psychiatric disorders, extensive comorbidities, underdiagnosis or overdiagnosis, multifaceted interactions with genetics and life experiences, and the highly multivariate nature of neural correlates. Therefore, increasingly larger datasets that measure more variables in larger cohorts are needed to gain insights. In this review, we present current "best practice" approaches for using existing databases, collecting and sharing new repositories for big data analyses, and future directions for big data in neuroimaging and psychiatry with an emphasis on contributing to collaborative efforts and the challenges of multi-study data analysis.

Evaluation of ComBat harmonization for reducing across-tracer differences in regional amyloid PET analyses

Introduction

Differences in amyloid positron emission tomography (PET) radiotracer pharmacokinetics and binding properties lead to discrepancies in amyloid-β uptake estimates. Harmonization of tracer-specific biases is crucial for optimal performance of downstream tasks. Here, we investigated the efficacy of ComBat, a data-driven harmonization model, for reducing tracer-specific biases in regional amyloid PET measurements from [18F]-florbetapir (FBP) and [11C]-Pittsburgh Compound-B (PiB).

Methods

One-hundred-thirteen head-to-head FBP-PiB scan pairs, scanned from the same subject within ninety days, were selected from the Open Access Series of Imaging Studies 3 (OASIS-3) dataset. The Centiloid scale, ComBat with no covariates, ComBat with biological covariates, and GAM-ComBat with biological covariates were used to harmonize both global and regional amyloid standardized uptake value ratios (SUVR). Variants of ComBat, including longitudinal ComBat and PEACE, were also tested. Intraclass correlation coefficient (ICC) and mean absolute error (MAE) were computed to measure the absolute agreement between tracers. Additionally, longitudinal amyloid SUVRs from an anti-amyloid drug trial were simulated using linear mixed effects modeling. Differences in rates-of-change between simulated treatment and placebo groups were tested, and change in statistical power/Type-I error after harmonization was quantified.

Results

In the head-to-head tracer comparison, ComBat with no covariates was the best at increasing ICC and decreasing MAE of both global summary and regional amyloid PET SUVRs between scan pairs of the same group of subjects. In the clinical trial simulation, harmonization with both Centiloid and ComBat increased statistical power of detecting true rate-of-change differences between groups and decreased false discovery rate in the absence of a treatment effect. The greatest benefit of harmonization was observed when groups exhibited differing FBP-to-PiB proportions.

Conclusion

ComBat outperformed the Centiloid scale in harmonizing both global and regional amyloid estimates. Additionally, ComBat improved the detection of rate-of-change differences between clinical trial groups. Our findings suggest that ComBat is a viable alternative to Centiloid for harmonizing regional amyloid PET analyses.

Modeling the impact of MRI acquisition bias on structural connectomes: Harmonizing structural connectomes

One way to increase the statistical power and generalizability of neuroimaging studies is to collect data at multiple sites or merge multiple cohorts. However, this usually comes with site-related biases due to the heterogeneity of scanners and acquisition parameters, negatively impacting sensitivity. Brain structural connectomes are not an exception: Being derived from T1-weighted and diffusion-weighted magnetic resonance images, structural connectivity is impacted by differences in imaging protocol. Beyond minimizing acquisition parameter differences, removing bias with postprocessing is essential. In this work we create, from the exhaustive Human Connectome Project Young Adult dataset, a resampled dataset of different b-values and spatial resolutions, modeling a cohort scanned across multiple sites. After demonstrating the statistical impact of acquisition parameters on connectivity, we propose a linear regression with explicit modeling of b-value and spatial resolution, and validate its performance on separate datasets. We show that b-value and spatial resolution affect connectivity in different ways and that acquisition bias can be reduced using a linear regression informed by the acquisition parameters while retaining interindividual differences and hence boosting fingerprinting performance. We also demonstrate the generative potential of our model, and its generalization capability in an independent dataset reflective of typical acquisition practices in clinical settings.

Longitudinal development of resting-state functional connectivity in adolescents with and without internalizing disorders

Longitudinal studies using resting-state functional magnetic resonance imaging (rs-fMRI) focused on adolescent internalizing psychopathology are scarce and have mostly investigated standardized treatment effects on functional connectivity (FC) of the full amygdala. The role of amygdala subregions and large resting-state networks had yet to be elucidated, and treatment is in practice often personalized. Here, longitudinal FC development of amygdala subregions and whole-brain networks are investigated in a clinically representative sample. Treatment-naïve adolescents with clinical depression and comorbid anxiety who started care-as-usual (n = 23; INT) and healthy controls (n = 24; HC) participated in rs-fMRI scans and questionnaires at baseline (before treatment) and after three months. Changes between and within groups over time in FC of the laterobasal amygdala (LBA), centromedial amygdala (CMA) and whole-brain networks derived from independent component analysis (ICA) were investigated. Groups differed significantly in FC development of the right LBA to the postcentral gyrus and the left LBA to the frontal pole. Within INT, FC to the frontal pole and postcentral gyrus changed over time while changes in FC of the right LBA were also linked to symptom change. No significant interactions were observed when considering FC from CMA bilateral seeds or within ICA-derived networks. Results in this cohort suggest divergent longitudinal development of FC from bilateral LBA subregions in adolescents with internalizing disorders compared to healthy peers, possibly reflecting nonspecific treatment effects. Moreover, associations were found with symptom change. These results highlight the importance of differentiation of amygdala subregions in neuroimaging research in adolescents.

ComBatLS: A location- and scale-preserving method for multi-site image harmonization

Recent work has leveraged massive datasets and advanced harmonization methods to construct normative models of neuroanatomical features and benchmark individuals' morphology. However, current harmonization tools do not preserve the effects of biological covariates including sex and age on features' variances; this failure may induce error in normative scores, particularly when such factors are distributed unequally across sites. Here, we introduce a new extension of the popular ComBat harmonization method, ComBatLS, that preserves biological variance in features' locations and scales. We use UK Biobank data to show that ComBatLS robustly replicates individuals' normative scores better than other ComBat methods when subjects are assigned to sex-imbalanced synthetic "sites". Additionally, we demonstrate that ComBatLS significantly reduces sex biases in normative scores compared to traditional methods. Finally, we show that ComBatLS successfully harmonizes consortium data collected across over 50 studies. R implementation of ComBatLS is available at https://github.com/andy1764/ComBatFamily.

White matter trajectories over the lifespan

White matter (WM) changes occur throughout the lifespan at a different rate for each developmental period. We aggregated 10879 structural MRIs and 6186 diffusion-weighted MRIs from participants between 2 weeks to 100 years of age. Age-related changes in gray matter and WM partial volumes and microstructural WM properties, both brain-wide and on 29 reconstructed tracts, were investigated as a function of biological sex and hemisphere, when appropriate. We investigated the curve fit that would best explain age-related differences by fitting linear, cubic, quadratic, and exponential models to macro and microstructural WM properties. Following the first steep increase in WM volume during infancy and childhood, the rate of development slows down in adulthood and decreases with aging. Similarly, microstructural properties of WM, particularly fractional anisotropy (FA) and mean diffusivity (MD), follow independent rates of change across the lifespan. The overall increase in FA and decrease in MD are modulated by demographic factors, such as the participant's age, and show different hemispheric asymmetries in some association tracts reconstructed via probabilistic tractography. All changes in WM macro and microstructure seem to follow nonlinear trajectories, which also differ based on the considered metric. Exponential changes occurred for the WM volume and FA and MD values in the first five years of life. Collectively, these results provide novel insight into how changes in different metrics of WM occur when a lifespan approach is considered.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

Intrusive Traumatic Re-Experiencing Domain: Functional Connectivity Feature Classification by the ENIGMA PTSD Consortium

Background

Intrusive traumatic re-experiencing domain (ITRED) was recently introduced as a novel perspective on posttraumatic psychopathology, proposing to focus research of posttraumatic stress disorder (PTSD) on the unique symptoms of intrusive and involuntary re-experiencing of the trauma, namely, intrusive memories, nightmares, and flashbacks. The aim of the present study was to explore ITRED from a neural network connectivity perspective.

Methods

Data were collected from 9 sites taking part in the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) PTSD Consortium (n= 584) and included itemized PTSD symptom scores and resting-state functional connectivity (rsFC) data. We assessed the utility of rsFC in classifying PTSD, ITRED-only (no PTSD diagnosis), and trauma-exposed (TE)-only (no PTSD or ITRED) groups using a machine learning approach, examining well-known networks implicated in PTSD. A random forest classification model was built on a training set using cross-validation, and the averaged cross-validation model performance for classification was evaluated using the area under the curve. The model was tested using a fully independent portion of the data (test dataset), and the test area under the curve was evaluated.

Results

rsFC signatures differentiated TE-only participants from PTSD and ITRED-only participants at about 60% accuracy. Conversely, rsFC signatures did not differentiate PTSD from ITRED-only individuals (45% accuracy). Common features differentiating TE-only participants from PTSD and ITRED-only participants mainly involved default mode network-related pathways. Some unique features, such as connectivity within the frontoparietal network, differentiated TE-only participants from one group (PTSD or ITRED-only) but to a lesser extent from the other group.

Conclusions

Neural network connectivity supports ITRED as a novel neurobiologically based approach to classifying posttrauma psychopathology.

Source-based morphometry reveals structural brain pattern abnormalities in 22q11.2 deletion syndrome

22q11.2 deletion syndrome (22q11DS) is the most frequently occurring microdeletion in humans. It is associated with a significant impact on brain structure, including prominent reductions in gray matter volume (GMV), and neuropsychiatric manifestations, including cognitive impairment and psychosis. It is unclear whether GMV alterations in 22q11DS occur according to distinct structural patterns. Then, 783 participants (470 with 22q11DS: 51% females, mean age [SD] 18.2 [9.2]; and 313 typically developing [TD] controls: 46% females, mean age 18.0 [8.6]) from 13 datasets were included in the present study. We segmented structural T1-weighted brain MRI scans and extracted GMV images, which were then utilized in a novel source-based morphometry (SBM) pipeline (SS-Detect) to generate structural brain patterns (SBPs) that capture co-varying GMV. We investigated the impact of the 22q11.2 deletion, deletion size, intelligence quotient, and psychosis on the SBPs. Seventeen GMV-SBPs were derived, which provided spatial patterns of GMV covariance associated with a quantitative metric (i.e., loading score) for analysis. Patterns of topographically widespread differences in GMV covariance, including the cerebellum, discriminated individuals with 22q11DS from healthy controls. The spatial extents of the SBPs that revealed disparities between individuals with 22q11DS and controls were consistent with the findings of the univariate voxel-based morphometry analysis. Larger deletion size was associated with significantly lower GMV in frontal and occipital SBPs; however, history of psychosis did not show a strong relationship with these covariance patterns. 22q11DS is associated with distinct structural abnormalities captured by topographical GMV covariance patterns that include the cerebellum. Findings indicate that structural anomalies in 22q11DS manifest in a nonrandom manner and in distinct covarying anatomical patterns, rather than a diffuse global process. These SBP abnormalities converge with previously reported cortical surface area abnormalities, suggesting disturbances of early neurodevelopment as the most likely underlying mechanism.

A Brain Morphometry Study with Across-Site Harmonization Using a ComBat-Generalized Additive Model in Children and Adolescents

Regional anatomical structures of the brain are intimately connected to functions corresponding to specific regions and the temporospatial pattern of genetic expression and their functions from the fetal period to old age. Therefore, quantitative brain morphometry has often been employed in neuroscience investigations, while controlling for the scanner effect of the scanner is a critical issue for ensuring accuracy in brain morphometric studies of rare orphan diseases due to the lack of normal reference values available for multicenter studies. This study aimed to provide across-site normal reference values of global and regional brain volumes for each sex and age group in children and adolescents. We collected magnetic resonance imaging (MRI) examinations of 846 neurotypical participants aged 6.0-17.9 years (339 male and 507 female participants) from 5 institutions comprising healthy volunteers or neurotypical patients without neurological disorders, neuropsychological disorders, or epilepsy. Regional-based analysis using the CIVET 2.1.0. pipeline provided regional brain volumes, and the measurements were across-site combined using ComBat-GAM harmonization. The normal reference values of global and regional brain volumes and lateral indices in our study could be helpful for evaluating the characteristics of the brain morphology of each individual in a clinical setting and investigating the brain morphology of ultra-rare diseases.

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.

Multiscale functional connectivity patterns of the aging brain learned from harmonized rsfMRI data of the multi-cohort iSTAGING study

To learn multiscale functional connectivity patterns of the aging brain, we built a brain age prediction model of functional connectivity measures at seven scales on a large fMRI dataset, consisting of resting-state fMRI scans of 4186 individuals with a wide age range (22 to 97 years, with an average of 63) from five cohorts. We computed multiscale functional connectivity measures of individual subjects using a personalized functional network computational method, harmonized the functional connectivity measures of subjects from multiple datasets in order to build a functional brain age model, and finally evaluated how functional brain age gap correlated with cognitive measures of individual subjects. Our study has revealed that functional connectivity measures at multiple scales were more informative than those at any single scale for the brain age prediction, the data harmonization significantly improved the brain age prediction performance, and the data harmonization in the functional connectivity measures' tangent space worked better than in their original space. Moreover, brain age gap scores of individual subjects derived from the brain age prediction model were significantly correlated with clinical and cognitive measures. Overall, these results demonstrated that multiscale functional connectivity patterns learned from a large-scale multi-site rsfMRI dataset were informative for characterizing the aging brain and the derived brain age gap was associated with cognitive and clinical measures.

Site effects how-to and when: An overview of retrospective techniques to accommodate site effects in multi-site neuroimaging analyses

Site differences, or systematic differences in feature distributions across multiple data-acquisition sites, are a known source of heterogeneity that may adversely affect large-scale meta- and mega-analyses of independently collected neuroimaging data. They influence nearly all multi-site imaging modalities and biomarkers, and methods to compensate for them can improve reliability and generalizability in the analysis of genetics, omics, and clinical data. The origins of statistical site effects are complex and involve both technical differences (scanner vendor, head coil, acquisition parameters, imaging processing) and differences in sample characteristics (inclusion/exclusion criteria, sample size, ancestry) between sites. In an age of expanding international consortium research, there is a growing need to disentangle technical site effects from sample characteristics of interest. Numerous statistical and machine learning methods have been developed to control for, model, or attenuate site effects - yet to date, no comprehensive review has discussed the benefits and drawbacks of each for different use cases. Here, we provide an overview of the different existing statistical and machine learning methods developed to remove unwanted site effects from independently collected neuroimaging samples. We focus on linear mixed effect models, the ComBat technique and its variants, adjustments based on image quality metrics, normative modeling, and deep learning approaches such as generative adversarial networks. For each method, we outline the statistical foundation and summarize strengths and weaknesses, including their assumptions and conditions of use. We provide information on software availability and comment on the ease of use and the applicability of these methods to different types of data. We discuss validation and comparative reports, mention caveats and provide guidance on when to use each method, depending on context and specific research questions.

Identifying and revealing different brain neural activities of cognitive subtypes in early course schizophrenia

Background

Cognitive subtypes of schizophrenia may exhibit different neurobiological characteristics. This study aimed to reveal the underlying neurobiological features between cognitive subtypes in the early course of schizophrenia (ECS). According to prior studies, we hypothesized to identify 2–4 distinct cognitive subtypes. We further hypothesized that the subtype with relatively poorer cognitive function might have lower brain spontaneous neural activity than the subtype with relatively better cognitive function.

Method

Cognitive function was assessed by the MATRICS Consensus Cognitive Battery (MCCB). Resting-state functional magnetic resonance imaging scanning was conducted for each individual. There were 155 ECS individuals and 97 healthy controls (HCs) included in the subsequent analysis. Latent profile analysis (LPA) was used to identify the cognitive subtypes in ECS individuals, and amplitude of low-frequency fluctuations (ALFFs) was used to measure brain spontaneous neural activity in ECS individuals and HCs.

Results

LPA identified two cognitive subtypes in ECS individuals, containing a severely impaired subtype (SI, n = 63) and a moderately impaired subtype (MI, n = 92). Compared to HCs, ECS individuals exhibited significantly increased ALFF in the left caudate and bilateral thalamus and decreased ALFF in the bilateral medial prefrontal cortex and bilateral posterior cingulate cortex/precuneus (PCC/PCu). In ECS cognitive subtypes, SI showed significantly higher ALFF in the left precentral gyrus (PreCG) and lower ALFF in the left PCC/PCu than MI. Furthermore, ALFFs of left PreCG were negatively correlated with several MCCB cognitive domains in ECS individuals, while ALFF of left PCC/PCu presented opposite correlations.

Conclusion

Our findings suggest that differences in the brain spontaneous neural activity of PreCG and PCC/PCu might be the potential neurobiological features of the cognitive subtypes in ECS, which may deepen our understanding of the role of PreCG and PCC/PCu in the pathogenesis of cognitive impairment in schizophrenia.