Large-scale analysis of structural brain asymmetries in schizophrenia via the ENIGMA consortium.

Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

Identifying covariate-related subnetworks for whole-brain connectome analysis.

Whole-brain connectome data characterize the connections among distributed neural populations as a set of edges in a large network, and neuroscience research aims to systematically investigate associations between brain connectome and clinical or experimental conditions as covariates. A covariate is often related to a number of edges connecting multiple brain areas in an organized structure. However, in practice, neither the covariate-related edges nor the structure is known. Therefore, the understanding of underlying neural mechanisms relies on statistical methods that are capable of simultaneously identifying covariate-related connections and recognizing their network topological structures. The task can be challenging because of false-positive noise and almost infinite possibilities of edges combining into subnetworks. To address these challenges, we propose a new statistical approach to handle multivariate edge variables as outcomes and output covariate-related subnetworks. We first study the graph properties of covariate-related subnetworks from a graph and combinatorics perspective and accordingly bridge the inference for individual connectome edges and covariate-related subnetworks. Next, we develop efficient algorithms to exact covariate-related subnetworks from the whole-brain connectome data with an $\ell_0$ norm penalty. We validate the proposed methods based on an extensive simulation study, and we benchmark our performance against existing methods. Using our proposed method, we analyze two separate resting-state functional magnetic resonance imaging data sets for schizophrenia research and obtain highly replicable disease-related subnetworks.

Connectome architecture shapes large-scale cortical alterations in schizophrenia: a worldwide ENIGMA study.

Schizophrenia is a prototypical network disorder with widespread brain-morphological alterations, yet it remains unclear whether these distributed alterations robustly reflect the underlying network layout. We tested whether large-scale structural alterations in schizophrenia relate to normative structural and functional connectome architecture, and systematically evaluated robustness and generalizability of these network-level alterations. Leveraging anatomical MRI scans from 2439 adults with schizophrenia and 2867 healthy controls from 26 ENIGMA sites and normative data from the Human Connectome Project (n = 207), we evaluated structural alterations of schizophrenia against two network susceptibility models: (i) hub vulnerability, which examines associations between regional network centrality and magnitude of disease-related alterations; (ii) epicenter mapping, which identifies regions whose typical connectivity profile most closely resembles the disease-related morphological alterations. To assess generalizability and specificity, we contextualized the influence of site, disease stages, and individual clinical factors and compared network associations of schizophrenia with that found in affective disorders. Our findings show schizophrenia-related cortical thinning is spatially associated with functional and structural hubs, suggesting that highly interconnected regions are more vulnerable to morphological alterations. Predominantly temporo-paralimbic and frontal regions emerged as epicenters with connectivity profiles linked to schizophrenia's alteration patterns. Findings were robust across sites, disease stages, and related to individual symptoms. Moreover, transdiagnostic comparisons revealed overlapping epicenters in schizophrenia and bipolar, but not major depressive disorder, suggestive of a pathophysiological continuity within the schizophrenia-bipolar-spectrum. In sum, cortical alterations over the course of schizophrenia robustly follow brain network architecture, emphasizing marked hub susceptibility and temporo-frontal epicenters at both the level of the group and the individual. Subtle variations of epicenters across disease stages suggest interacting pathological processes, while associations with patient-specific symptoms support additional inter-individual variability of hub vulnerability and epicenters in schizophrenia. Our work outlines potential pathways to better understand macroscale structural alterations, and inter- individual variability in schizophrenia.

Cross-cohort replicable resting-state functional connectivity in predicting symptoms and cognition of schizophrenia.

Schizophrenia (SZ) is a debilitating mental illness characterized by adolescence or early adulthood onset of psychosis, positive and negative symptoms, as well as cognitive impairments. Despite a plethora of studies leveraging functional connectivity (FC) from functional magnetic resonance imaging (fMRI) to predict symptoms and cognitive impairments of SZ, the findings have exhibited great heterogeneity. We aimed to identify congruous and replicable connectivity patterns capable of predicting positive and negative symptoms as well as cognitive impairments in SZ. Predictable functional connections (FCs) were identified by employing an individualized prediction model, whose replicability was further evaluated across three independent cohorts (BSNIP, SZ = 174; COBRE, SZ = 100; FBIRN, SZ = 161). Across cohorts, we observed that altered FCs in frontal-temporal-cingulate-thalamic network were replicable in prediction of positive symptoms, while sensorimotor network was predictive of negative symptoms. Temporal-parahippocampal network was consistently identified to be associated with reduced cognitive function. These replicable 23 FCs effectively distinguished SZ from healthy controls (HC) across three cohorts (82.7%, 90.2%, and 86.1%). Furthermore, models built using these replicable FCs showed comparable accuracies to those built using the whole-brain features in predicting symptoms/cognition of SZ across the three cohorts (r = .17-.33, p < .05). Overall, our findings provide new insights into the neural underpinnings of SZ symptoms/cognition and offer potential targets for further research and possible clinical interventions.

Reciprocal relationships between stress and depressive symptoms: the essential role of the nucleus accumbens.

BACKGROUND: Stress and depression have a reciprocal relationship, but the neural underpinnings of this reciprocity are unclear. We investigated neuroimaging phenotypes that facilitate the reciprocity between stress and depressive symptoms. METHODS: In total, 22 195 participants (52.0% females) from the population-based UK Biobank study completed two visits (initial visit: 2006-2010, age = 55.0 ± 7.5 [40-70] years; second visit: 2014-2019; age = 62.7 ± 7.5 [44-80] years). Structural equation modeling was used to examine the longitudinal relationship between self-report stressful life events (SLEs) and depressive symptoms. Cross-sectional data were used to examine the overlap between neuroimaging correlates of SLEs and depressive symptoms on the second visit among 138 multimodal imaging phenotypes. RESULTS: Longitudinal data were consistent with significant bidirectional causal relationship between SLEs and depressive symptoms. In cross-sectional analyses, SLEs were significantly associated with lower bilateral nucleus accumbal volume and lower fractional anisotropy of the forceps major. Depressive symptoms were significantly associated with extensive white matter hyperintensities, thinner cortex, lower subcortical volume, and white matter microstructural deficits, mainly in corticostriatal-limbic structures. Lower bilateral nucleus accumbal volume were the only imaging phenotypes with overlapping effects of depressive symptoms and SLEs (B = -0.032 to -0.023, p = 0.006-0.034). Depressive symptoms and SLEs significantly partially mediated the effects of each other on left and right nucleus accumbens volume (proportion of effects mediated = 12.7-14.3%, p < 0.001-p = 0.008). For the left nucleus accumbens, post-hoc seed-based analysis showed lower resting-state functional connectivity with the left orbitofrontal cortex (cluster size = 83 voxels, p = 5.4 × 10-5) in participants with high v. no SLEs. CONCLUSIONS: The nucleus accumbens may play a key role in the reciprocity between stress and depressive symptoms.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Brain ageing in schizophrenia: evidence from 26 international cohorts via the ENIGMA Schizophrenia consortium.

Schizophrenia (SZ) is associated with an increased risk of life-long cognitive impairments, age-related chronic disease, and premature mortality. We investigated evidence for advanced brain ageing in adult SZ patients, and whether this was associated with clinical characteristics in a prospective meta-analytic study conducted by the ENIGMA Schizophrenia Working Group. The study included data from 26 cohorts worldwide, with a total of 2803 SZ patients (mean age 34.2 years; range 18-72 years; 67% male) and 2598 healthy controls (mean age 33.8 years, range 18-73 years, 55% male). Brain-predicted age was individually estimated using a model trained on independent data based on 68 measures of cortical thickness and surface area, 7 subcortical volumes, lateral ventricular volumes and total intracranial volume, all derived from T1-weighted brain magnetic resonance imaging (MRI) scans. Deviations from a healthy brain ageing trajectory were assessed by the difference between brain-predicted age and chronological age (brain-predicted age difference [brain-PAD]). On average, SZ patients showed a higher brain-PAD of +3.55 years (95% CI: 2.91, 4.19; I2 = 57.53%) compared to controls, after adjusting for age, sex and site (Cohen's d = 0.48). Among SZ patients, brain-PAD was not associated with specific clinical characteristics (age of onset, duration of illness, symptom severity, or antipsychotic use and dose). This large-scale collaborative study suggests advanced structural brain ageing in SZ. Longitudinal studies of SZ and a range of mental and somatic health outcomes will help to further evaluate the clinical implications of increased brain-PAD and its ability to be influenced by interventions.

In vivo white matter microstructure in adolescents with early-onset psychosis: a multi-site mega-analysis.

Emerging evidence suggests brain white matter alterations in adolescents with early-onset psychosis (EOP; age of onset <18 years). However, as neuroimaging methods vary and sample sizes are modest, results remain inconclusive. Using harmonized data processing protocols and a mega-analytic approach, we compared white matter microstructure in EOP and healthy controls using diffusion tensor imaging (DTI). Our sample included 321 adolescents with EOP (median age = 16.6 years, interquartile range (IQR) = 2.14, 46.4% females) and 265 adolescent healthy controls (median age = 16.2 years, IQR = 2.43, 57.7% females) pooled from nine sites. All sites extracted mean fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) for 25 white matter regions of interest per participant. ComBat harmonization was performed for all DTI measures to adjust for scanner differences. Multiple linear regression models were fitted to investigate case-control differences and associations with clinical variables in regional DTI measures. We found widespread lower FA in EOP compared to healthy controls, with the largest effect sizes in the superior longitudinal fasciculus (Cohen's d = 0.37), posterior corona radiata (d = 0.32), and superior fronto-occipital fasciculus (d = 0.31). We also found widespread higher RD and more localized higher MD and AD. We detected significant effects of diagnostic subgroup, sex, and duration of illness, but not medication status. Using the largest EOP DTI sample to date, our findings suggest a profile of widespread white matter microstructure alterations in adolescents with EOP, most prominently in male individuals with early-onset schizophrenia and individuals with a shorter duration of illness.

A multivariate to multivariate approach for voxel-wise genome-wide association analysis.

The joint analysis of imaging-genetics data facilitates the systematic investigation of genetic effects on brain structures and functions with spatial specificity. We focus on voxel-wise genome-wide association analysis, which may involve trillions of single nucleotide polymorphism (SNP)-voxel pairs. We attempt to identify underlying organized association patterns of SNP-voxel pairs and understand the polygenic and pleiotropic networks on brain imaging traits. We propose a bi-clique graph structure (ie, a set of SNPs highly correlated with a cluster of voxels) for the systematic association pattern. Next, we develop computational strategies to detect latent SNP-voxel bi-cliques and an inference model for statistical testing. We further provide theoretical results to guarantee the accuracy of our computational algorithms and statistical inference. We validate our method by extensive simulation studies, and then apply it to the whole genome genetic and voxel-level white matter integrity data collected from 1052 participants of the human connectome project. The results demonstrate multiple genetic loci influencing white matter integrity measures on splenium and genu of the corpus callosum.

An in-depth association analysis of genetic variants within nicotine-related loci: Meeting in middle of GWAS and genetic fine-mapping.

In the last two decades of Genome-wide association studies (GWAS), nicotine-dependence-related genetic loci (e.g., nicotinic acetylcholine receptor - nAChR subunit genes) are among the most replicable genetic findings. Although GWAS results have reported tens of thousands of SNPs within these loci, further analysis (e.g., fine-mapping) is required to identify the causal variants. However, it is computationally challenging for existing fine-mapping methods to reliably identify causal variants from thousands of candidate SNPs based on the posterior inclusion probability. To address this challenge, we propose a new method to select SNPs by jointly modeling the SNP-wise inference results and the underlying structured network patterns of the linkage disequilibrium (LD) matrix. We use adaptive dense subgraph extraction method to recognize the latent network patterns of the LD matrix and then apply group LASSO to select causal variant candidates. We applied this new method to the UK biobank data to identify the causal variant candidates for nicotine addiction. Eighty-one nicotine addiction-related SNPs (i.e.,-log(p) > 50) of nAChR were selected, which are highly correlated (average r2>0.8) although they are physically distant (e.g., >200 kilobase away) and from various genes. These findings revealed that distant SNPs from different genes can show higher LD r2 than their neighboring SNPs, and jointly contribute to a complex trait like nicotine addiction.

Association between brain similarity to severe mental illnesses and comorbid cerebral, physical, and cognitive impairments.

Severe mental illnesses (SMIs) are often associated with compromised brain health, physical comorbidities, and cognitive deficits, but it is incompletely understood whether these comorbidities are intrinsic to SMI pathophysiology or secondary to having SMIs. We tested the hypothesis that cerebral, cardiometabolic, and cognitive impairments commonly observed in SMIs can be observed in non-psychiatric individuals with SMI-like brain patterns of deviation as seen on magnetic resonance imaging. 22,883 participants free of common neuropsychiatric conditions from the UK Biobank (age = 63.4 ± 7.5 years, range = 45-82 years, 50.9% female) were split into discovery and replication samples. The regional vulnerability index (RVI) was used to quantify each participant's respective brain similarity to meta-analytical patterns of schizophrenia spectrum disorder, bipolar disorder, and major depressive disorder in gray matter thickness, subcortical gray matter volume, and white matter integrity. Cluster analysis revealed five clusters with distinct RVI profiles. Compared with a cluster with no RVI elevation, a cluster with RVI elevation across all SMIs and brain structures showed significantly higher volume of white matter hyperintensities (Cohen's d = 0.59, pFDR < 10-16), poorer cardiovascular (Cohen's d = 0.30, pFDR < 10-16) and metabolic (Cohen's d = 0.12, pFDR = 1.3 × 10-4) health, and slower speed of information processing (|Cohen's d| = 0.11-0.17, pFDR = 1.6 × 10-3-4.6 × 10-8). This cluster also had significantly higher level of C-reactive protein and alcohol use (Cohen's d = 0.11 and 0.28, pFDR = 4.1 × 10-3 and 1.1 × 10-11). Three other clusters with respective RVI elevation in gray matter thickness, subcortical gray matter volume, and white matter integrity showed intermediate level of white matter hyperintensities, cardiometabolic health, and alcohol use. Our results suggest that cerebral, physical, and cognitive impairments in SMIs may be partly intrinsic via shared pathophysiological pathways with SMI-related brain anatomical changes.

Brain deficit patterns of metabolic illnesses overlap with those for major depressive disorder: A new metric of brain metabolic disease.

Metabolic illnesses (MET) are detrimental to brain integrity and are common comorbidities in patients with mental illnesses, including major depressive disorder (MDD). We quantified effects of MET on standard regional brain morphometric measures from 3D brain MRI as well as diffusion MRI in a large sample of UK BioBank participants. The pattern of regional effect sizes of MET in non-psychiatric UKBB subjects was significantly correlated with the spatial profile of regional effects reported by the largest meta-analyses in MDD but not in bipolar disorder, schizophrenia or Alzheimer's disease. We used a regional vulnerability index (RVI) for MET (RVI-MET) to measure individual's brain similarity to the expected patterns in MET in the UK Biobank sample. Subjects with MET showed a higher effect size for RVI-MET than for any of the individual brain measures. We replicated elevation of RVI-MET in a sample of MDD participants with MET versus non-MET. RVI-MET scores were significantly correlated with the volume of white matter hyperintensities, a neurological consequence of MET and age, in both groups. Higher RVI-MET in both samples was associated with obesity, tobacco smoking and frequent alcohol use but was unrelated to antidepressant use. In summary, MET effects on the brain were regionally specific and individual similarity to the pattern was more strongly associated with MET than any regional brain structural metric. Effects of MET overlapped with the reported brain differences in MDD, likely due to higher incidence of MET, smoking and alcohol use in subjects with MDD.

Deciphering the causal relationship between blood pressure and regional white matter integrity: A two-sample Mendelian randomization study.

Elevated arterial blood pressure (BP) is a common risk factor for cerebrovascular and cardiovascular diseases, but no causal relationship has been established between BP and cerebral white matter (WM) integrity. In this study, we performed a two-sample Mendelian randomization (MR) analysis with individual-level data by defining two nonoverlapping sets of European ancestry individuals (genetics-exposure set: N = 203,111; mean age = 56.71 years, genetics-outcome set: N = 16,156; mean age = 54.61 years) from UK Biobank to evaluate the causal effects of BP on regional WM integrity, measured by fractional anisotropy of diffusion tensor imaging. Two BP traits: systolic and diastolic blood pressure were used as exposures. Genetic variant was carefully selected as instrumental variable (IV) under the MR analysis assumptions. We existing large-scale genome-wide association study summary data for validation. The main method used was a generalized version of inverse-variance weight method while other MR methods were also applied for consistent findings. Two additional MR analyses were performed to exclude the possibility of reverse causality. We found significantly negative causal effects (FDR-adjusted p < .05; every 10 mmHg increase in BP leads to a decrease in FA value by .4% ~ 2%) of BP traits on a union set of 17 WM tracts, including brain regions related to cognitive function and memory. Our study extended the previous findings of association to causation for regional WM integrity, providing insights into the pathological processes of elevated BP that might chronically alter the brain microstructure in different regions.

Frontal white matter association with sleep quality and the role of stress.

An important measure of brain health is the integrity of white matter connectivity structures that link brain regions. Studies have found an association between poorer sleep quality and decreased white matter integrity. Stress is among the strongest predictors of sleep quality. This study aimed to evaluate the association between sleep quality and white matter and to test if the relationship persisted after accounting for stress. White matter microstructures were measured by diffusion tensor imaging in a population of Old Order Amish/Mennonite (N = 240). Sleep quality was determined by the Pittsburgh Sleep Quality Index. Current stress levels were measured by the perceived stress scale. Exposure to lifetime stress was measured by the lifetime stressor inventory. Microstructures of four white matter tracts: left and right anterior limbs of internal capsule, left anterior corona radiata, and genu of corpus callosum were significantly correlated with sleep quality (all p ≤ 0.001). The current stress level was a significant predictor of sleep quality (p ≤ 0.001) while lifetime stress was not. PSQI remained significantly associated with white matter integrity in these frontal tracts (all p < 0.01) after accounting for current stress and lifetime stress, while current and lifetime stress were not significant predictors of white matter in any of the four models. Sleep quality did not have any substantial mediation role between stress and white matter integrity. Sleep quality was significantly associated with several frontal white matter tracts that connect brain structures important for sleep regulation regardless of current or past stress levels.

History of suicide attempt associated with amygdala and hippocampus changes among individuals with schizophrenia.

Abnormalities in subcortical brain structures may reflect higher suicide risk in mood disorders, but less is known about its associations for schizophrenia. This cross-sectional imaging study aimed to explore whether the history of suicide attempts was associated with subcortical changes among individuals with schizophrenia. We recruited 44 individuals with schizophrenia and a history of suicide attempts (SZ-SA) and 44 individuals with schizophrenia but without a history of suicide attempts (SZ-NSA) and 44 healthy controls. Linear regression showed that SZ-SA had smaller volumes of the hippocampus (Cohen's d = -0.72), the amygdala (Cohen's d = -0.69), and some nuclei of the amygdala (Cohen's d, -0.57 to -0.72) than SZ-NSA after adjusting for age, sex, illness phase, and intracranial volume. There was no difference in the volume of the subfields of the hippocampus. It suggests the history of suicide attempts is associated with subcortical volume alterations in schizophrenia.

Association of cortical thickness and cognition with schizophrenia treatment resistance.

AIM: Approximately a third of patients with schizophrenia fail to adequately respond to antipsychotic medications, a condition known as treatment resistance (TR). We aimed to assess cognitive and cortical thickness deficits and their relationship to TR in schizophrenia. METHOD: We recruited patients with schizophrenia (n = 127), including patients at treatment initiation (n = 45), treatment-responsive patients (n = 40) and TR patients (n = 42), and healthy controls (n = 83). Clinical symptoms, neurocognitive function, and structural images were assessed. We performed group comparisons, and explored association of cortical thickness and cognition with TR. RESULTS: The TR patients showed significantly more severe clinical symptoms and cognitive impairment relative to the treatment-responsive group. Compared to healthy controls, 56 of 68 brain regions showed significantly reduced cortical thickness in patients with schizophrenia. Reductions in five regions were significantly associated with TR (reduction in TR relative to treatment-responsive patients), i.e. in the right caudal middle frontal gyrus, superior frontal cortex, fusiform gyrus, pars opercularis of the inferior frontal cortex, and supramarginal cortex. Cognition deficits were also significantly correlated with cortical thickness in these five regions in patients with schizophrenia. Cortical thickness of the right caudal middle frontal gyrus, superior frontal cortex and pars opercularis of the inferior frontal cortex also significantly mediated effects of cognitive deficits on TR. CONCLUSION: Treatment resistance in schizophrenia was associated with reduced thickness in the right caudal middle frontal gyrus, superior frontal cortex, fusiform gyrus, pars opercularis of the inferior frontal cortex, and supramarginal cortex. Cortical abnormalities further mediate cognitive deficits known to be associated with TR.

Mediation Analysis for High-Dimensional Mediators and Outcomes with an Application to Multimodal Imaging Data.

Multimodal neuroimaging data have attracted increasing attention for brain research. An integrated analysis of multimodal neuroimaging data and behavioral or clinical measurements provides a promising approach for comprehensively and systematically investigating the underlying neural mechanisms of different phenotypes. However, such an integrated data analysis is intrinsically challenging due to the complex interactive relationships between the multimodal multivariate imaging variables. To address this challenge, a novel multivariate-mediator and multivariate-outcome mediation model (MMO) is proposed to simultaneously extract the latent systematic mediation patterns and estimate the mediation effects based on a dense bi-cluster graph approach. A computationally efficient algorithm is developed for dense bicluster structure estimation and inference to identify the mediation patterns with multiple testing correction. The performance of the proposed method is evaluated by an extensive simulation analysis with comparison to the existing methods. The results show that MMO performs better in terms of both the false discovery rate and sensitivity compared to existing models. The MMO is applied to a multimodal imaging dataset from the Human Connectome Project to investigate the effect of systolic blood pressure on whole-brain imaging measures for the regional homogeneity of the blood oxygenation level-dependent signal through the cerebral blood flow.

Cerebral blood flow and cardiovascular risk effects on resting brain regional homogeneity.

Regional homogeneity (ReHo) is a measure of local functional brain connectivity that has been reported to be altered in a wide range of neuropsychiatric disorders. Computed from brain resting-state functional MRI time series, ReHo is also sensitive to fluctuations in cerebral blood flow (CBF) that in turn may be influenced by cerebrovascular health. We accessed cerebrovascular health with Framingham cardiovascular risk score (FCVRS). We hypothesize that ReHo signal may be influenced by regional CBF; and that these associations can be summarized as FCVRS→CBF→ReHo. We used three independent samples to test this hypothesis. A test-retest sample of N = 30 healthy volunteers was used for test-retest evaluation of CBF effects on ReHo. Amish Connectome Project (ACP) sample (N = 204, healthy individuals) was used to evaluate association between FCVRS and ReHo and testing if the association diminishes given CBF. The UKBB sample (N = 6,285, healthy participants) was used to replicate the effects of FCVRS on ReHo. We observed strong CBF→ReHo links (p<2.5 × 10-3) using a three-point longitudinal sample. In ACP sample, marginal and partial correlations analyses demonstrated that both CBF and FCVRS were significantly correlated with the whole-brain average (p<10-6) and regional ReHo values, with the strongest correlations observed in frontal, parietal, and temporal areas. Yet, the association between ReHo and FCVRS became insignificant once the effect of CBF was accounted for. In contrast, CBF→ReHo remained significantly linked after adjusting for FCVRS and demographic covariates (p<10-6). Analysis in N = 6,285 replicated the FCVRS→ReHo effect (p = 2.7 × 10-27). In summary, ReHo alterations in health and neuropsychiatric illnesses may be partially driven by region-specific variability in CBF, which is, in turn, influenced by cardiovascular factors.

An integrated cluster-wise significance measure for fMRI analysis.

Cluster-wise inference is widely used in fMRI analysis. The cluster-level statistic is often obtained by counting the number of intra-cluster voxels which surpass a voxel-level statistical significance threshold. This measure can be sub-optimal regarding the power and false-positive error rate because the suprathreshold voxel count neglects the voxel-wise significance levels and ignores the dependence between voxels. This article aims to provide a new Integrated Cluster-wise significance Measure (ICM) for cluster-level significance determination in cluster-wise fMRI analysis by integrating cluster extent, voxel-level significance (e.g., p values), and activation dependence between within-cluster voxels. We develop a computationally efficient strategy for ICM based on probabilistic approximation theories. Consequently, the computational load for ICM-based cluster-wise inference (e.g., permutation tests) is affordable. We validate the proposed method via extensive simulations and then apply it to two fMRI data sets. The results demonstrate that ICM can improve the power with well-controlled family-wise error (FWE).

A new multimodality fusion classification approach to explore the uniqueness of schizophrenia and autism spectrum disorder.

Schizophrenia (SZ) and autism spectrum disorder (ASD) sharing overlapping symptoms have a long history of diagnostic confusion. It is unclear what their differences at a brain level are. Here, we propose a multimodality fusion classification approach to investigate their divergence in brain function and structure. Using brain functional network connectivity (FNC) calculated from resting-state fMRI data and gray matter volume (GMV) estimated from sMRI data, we classify the two disorders using the main data (335 SZ and 380 ASD patients) via an unbiased 10-fold cross-validation pipeline, and also validate the classification generalization ability on an independent cohort (120 SZ and 349 ASD patients). The classification accuracy reached up to 83.08% for the testing data and 72.10% for the independent data, significantly better than the results from using the single-modality features. The discriminative FNCs that were automatically selected primarily involved the sub-cortical, default mode, and visual domains. Interestingly, all discriminative FNCs relating to the default mode network showed an intermediate strength in healthy controls (HCs) between SZ and ASD patients. Their GMV differences were mainly driven by the frontal gyrus, temporal gyrus, and insula. Regarding these regions, the mean GMV of HC fell intermediate between that of SZ and ASD, and ASD showed the highest GMV. The middle frontal gyrus was associated with both functional and structural differences. In summary, our work reveals the unique neuroimaging characteristics of SZ and ASD that can achieve high and generalizable classification accuracy, supporting their potential as disorder-specific neural substrates of the two entwined disorders.