Multi-site study of additive genetic effects on fractional anisotropy of cerebral white matter: Comparing meta and megaanalytical approaches for data pooling

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.

Unraveling multi-scale neuroimaging biomarkers and molecular foundations for schizophrenia: A combined multivariate pattern analysis and transcriptome-neuroimaging association study

AIMS

Schizophrenia is characterized by alterations in resting-state spontaneous brain activity; however, it remains uncertain whether variations at diverse spatial scales are capable of effectively distinguishing patients from healthy controls. Additionally, the genetic underpinnings of these alterations remain poorly elucidated. We aimed to address these questions in this study to gain better understanding of brain alterations and their underlying genetic factors in schizophrenia.

METHODS

A cohort of 103 individuals with diagnosed schizophrenia and 110 healthy controls underwent resting-state functional MRI scans. Spontaneous brain activity was assessed using the regional homogeneity (ReHo) metric at four spatial scales: voxel-level (Scale 1) and regional-level (Scales 2-4: 272, 53, 17 regions, respectively). For each spatial scale, multivariate pattern analysis was performed to classify schizophrenia patients from healthy controls, and a transcriptome-neuroimaging association analysis was performed to establish connections between gene expression data and ReHo alterations in schizophrenia.

RESULTS

The ReHo metrics at all spatial scales effectively discriminated schizophrenia from healthy controls. Scale 2 showed the highest classification accuracy at 84.6%, followed by Scale 1 (83.1%) and Scale 3 (78.5%), while Scale 4 exhibited the lowest accuracy (74.2%). Furthermore, the transcriptome-neuroimaging association analysis showed that there were not only shared but also unique enriched biological processes across the four spatial scales. These related biological processes were mainly linked to immune responses, inflammation, synaptic signaling, ion channels, cellular development, myelination, and transporter activity.

CONCLUSIONS

This study highlights the potential of multi-scale ReHo as a valuable neuroimaging biomarker in the diagnosis of schizophrenia. By elucidating the complex molecular basis underlying the ReHo alterations of this disorder, this study not only enhances our understanding of its pathophysiology, but also pave the way for future advancements in genetic diagnosis and treatment of schizophrenia.

Eye-brain connections revealed by multimodal retinal and brain imaging genetics

The retina, an anatomical extension of the brain, forms physiological connections with the visual cortex of the brain. Although retinal structures offer a unique opportunity to assess brain disorders, their relationship to brain structure and function is not well understood. In this study, we conducted a systematic cross-organ genetic architecture analysis of eye-brain connections using retinal and brain imaging endophenotypes. We identified novel phenotypic and genetic links between retinal imaging biomarkers and brain structure and function measures from multimodal magnetic resonance imaging (MRI), with many associations involving the primary visual cortex and visual pathways. Retinal imaging biomarkers shared genetic influences with brain diseases and complex traits in 65 genomic regions, with 18 showing genetic overlap with brain MRI traits. Mendelian randomization suggests bidirectional genetic causal links between retinal structures and neurological and neuropsychiatric disorders, such as Alzheimer's disease. Overall, our findings reveal the genetic basis for eye-brain connections, suggesting that retinal images can help uncover genetic risk factors for brain disorders and disease-related changes in intracranial structure and function.

A diffusion MRI tractography atlas for concurrent white matter mapping across Eastern and Western populations

The study of brain differences across Eastern and Western populations provides vital insights for understanding potential cultural and genetic influences on cognition and mental health. Diffusion MRI (dMRI) tractography is an important tool in assessing white matter (WM) connectivity and brain tissue microstructure across different populations. However, a comprehensive investigation into WM fiber tracts between Eastern and Western populations is challenged due to the lack of a cross-population WM atlas and the large site-specific variability of dMRI data. This study presents a dMRI tractography atlas, namely the East-West WM Atlas, for concurrent WM mapping between Eastern and Western populations and creates a large, harmonized dMRI dataset (n=306) based on the Human Connectome Project and the Chinese Human Connectome Project. The curated WM atlas, as well as subject-specific data including the harmonized dMRI data, the whole brain tractography data, and parcellated WM fiber tracts and their diffusion measures, are publicly released. This resource is a valuable addition to facilitating the exploration of brain commonalities and differences across diverse cultural backgrounds.

Multi-organ imaging-derived polygenic indexes for brain and body health

The UK Biobank (UKB) imaging project is a crucial resource for biomedical research, but is limited to 100,000 participants due to cost and accessibility barriers. Here we used genetic data to predict heritable imaging-derived phenotypes (IDPs) for a larger cohort. We developed and evaluated 4,375 IDP genetic scores (IGS) derived from UKB brain and body images. When applied to UKB participants who were not imaged, IGS revealed links to numerous phenotypes and stratified participants at increased risk for both brain and somatic diseases. For example, IGS identified individuals at higher risk for Alzheimer's disease and multiple sclerosis, offering additional insights beyond traditional polygenic risk scores of these diseases. When applied to independent external cohorts, IGS also stratified those at high disease risk in the All of Us Research Program and the Alzheimer's Disease Neuroimaging Initiative study. Our results demonstrate that, while the UKB imaging cohort is largely healthy and may not be the most enriched for disease risk management, it holds immense potential for stratifying the risk of various brain and body diseases in broader external genetic cohorts.

Chronotype and subjective sleep quality predict white matter integrity in young people with emerging mental disorders

Protecting brain health is a goal of early intervention. We explored whether sleep quality or chronotype could predict white matter (WM) integrity in emerging mental disorders. Young people (N = 364) accessing early-intervention clinics underwent assessments for chronotype, subjective sleep quality, and diffusion tensor imaging. Using machine learning, we examined whether chronotype or sleep quality (alongside diagnostic and demographic factors) could predict four measures of WM integrity: fractional anisotropy (FA), and radial, axial, and mean diffusivities (RD, AD and MD). We prioritised tracts that showed a univariate association with sleep quality or chronotype and considered predictors identified by ≥80% of machine learning (ML) models as 'important'. The most important predictors of WM integrity were demographics (age, sex and education) and diagnosis (depressive and bipolar disorders). Subjective sleep quality only predicted FA in the perihippocampal cingulum tract, whereas chronotype had limited predictive importance for WM integrity. To further examine links with mood disorders, we conducted a subgroup analysis. In youth with depressive and bipolar disorders, chronotype emerged as an important (often top-ranking) feature, predicting FA in the cingulum (cingulate gyrus), AD in the anterior corona radiata and genu of the corpus callosum, and RD in the corona radiata, anterior corona radiata, and genu of corpus callosum. Subjective quality was not important in this subgroup analysis. In summary, chronotype predicted altered WM integrity in the corona radiata and corpus callosum, whereas subjective sleep quality had a less significant role, suggesting that circadian factors may play a more prominent role in WM integrity in emerging mood disorders.

Heritability of cerebellar subregion volumes in adolescent and young adult twins

Twin studies have found gross cerebellar volume to be highly heritable. However, whether fine-grained regional volumes within the cerebellum are similarly heritable is still being determined. Anatomical MRI scans from two independent datasets (QTIM: Queensland Twin IMaging, N = 798, mean age 22.1 years; QTAB: Queensland Twin Adolescent Brain, N = 396, mean age 11.3 years) were combined with an optimised and automated cerebellum parcellation algorithm to segment and measure 28 cerebellar regions. We show that the heritability of regional volumetric measures varies widely across the cerebellum (h 2 $$ {h}^2 $$ 47%-91%). Additionally, the good to excellent test-retest reliability for a subsample of QTIM participants suggests that non-genetic variance in cerebellar volumes is due primarily to unique environmental influences rather than measurement error. We also show a consistent pattern of strong associations between the volumes of homologous left and right hemisphere regions. Associations were predominantly driven by genetic effects shared between lobules, with only sparse contributions from environmental effects. These findings are consistent with similar studies of the cerebrum and provide a first approximation of the upper bound of heritability detectable by genome-wide association studies.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) Study® has collected data from over 10,000 children across 21 sites, providing insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a dataset of harmonized and processed ABCD dMRI data (from release 3) has been created, comprising quality-controlled imaging data from 9,345 subjects, focusing exclusively on the baseline session, i.e., the first time point of the study. This resource required substantial computational time (approx. 50,000 CPU hours) for harmonization, whole-brain tractography, and white matter parcellation. The dataset includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts in full and low resolution, and 804 different dMRI-derived measures per subject (72.3 TB total size). Accessible via the NIMH Data Archive, it offers a large-scale dMRI dataset for studying structural connectivity in child and adolescent neurodevelopment. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Identification and validation of supervariants reveal novel loci associated with human white matter microstructure

As an essential part of the central nervous system, white matter coordinates communications between different brain regions and is related to a wide range of neurodegenerative and neuropsychiatric disorders. Previous genome-wide association studies (GWASs) have uncovered loci associated with white matter microstructure. However, GWASs suffer from limited reproducibility and difficulties in detecting multi-single-nucleotide polymorphism (multi-SNP) and epistatic effects. In this study, we adopt the concept of supervariants, a combination of alleles in multiple loci, to account for potential multi-SNP effects. We perform supervariant identification and validation to identify loci associated with 22 white matter fractional anisotropy phenotypes derived from diffusion tensor imaging. To increase reproducibility, we use United Kingdom (UK) Biobank White British (n = 30,842) data for discovery and internal validation, and UK Biobank White but non-British (n = 1927) data, Europeans from the Adolescent Brain Cognitive Development study (n = 4399) data, and Europeans from the Human Connectome Project (n = 319) data for external validation. We identify 23 novel loci on the discovery set that have not been reported in the previous GWASs on white matter microstructure. Among them, three supervariants on genomic regions 5q35.1, 8p21.2, and 19q13.32 have P-values lower than 0.05 in the meta-analysis of the three independent validation data sets. These supervariants contain genetic variants located in genes that have been related to brain structures, cognitive functions, and neuropsychiatric diseases. Our findings provide a better understanding of the genetic architecture underlying white matter microstructure.

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

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

Microstructural Organization of Distributed White Matter Associated With Fine Motor Control in US Service Members With Mild Traumatic Brain Injury

Mild traumatic brain injury (mTBI) is the most common form of brain injury. While most individuals recover from mTBI, roughly 20% experience persistent symptoms, potentially including reduced fine motor control. We investigate relationships between regional white matter organization and subcortical volumes associated with performance on the Grooved Pegboard (GPB) test in a large cohort of military Service Members and Veterans (SM&Vs) with and without a history of mTBI(s). Participants were enrolled in the Long-term Impact of Military-relevant Brain Injury Consortium-Chronic Effects of Neurotrauma Consortium. SM&Vs with a history of mTBI(s) (n = 847) and without mTBI (n = 190) underwent magnetic resonance imaging and the GPB test. We first examined between-group differences in GPB completion time. We then investigated associations between GPB performance and regional structural imaging measures (tractwise diffusivity, subcortical volumes, and cortical thickness) in SM&Vs with a history of mTBI(s). Lastly, we explored whether mTBI history moderated associations between imaging measures and GPB performance. SM&Vs with mTBI(s) performed worse than those without mTBI(s) on the non-dominant hand GPB test at a trend level (p < 0.1). Higher fractional anisotropy (FA) of tracts including the posterior corona radiata, superior longitudinal fasciculus, and uncinate fasciculus were associated with better GPB performance in the dominant hand in SM&Vs with mTBI(s). These findings support that the organization of several white matter bundles are associated with fine motor performance in SM&Vs. We did not observe that mTBI history moderated associations between regional FA and GPB test completion time, suggesting that chronic mTBI may not significantly influence fine motor control.

Structural impairment in superficial and deep white matter in schizophrenia

OBJECTIVE

Although disconnectivity among brain regions has been one of the main hypotheses for schizophrenia, the superficial white matter (SWM) has received less attention in schizophrenia research than the deep white matter (DWM) owing to the challenge of consistent reconstruction across subjects.

METHODS

We obtained the diffusion magnetic resonance imaging (dMRI) data of 223 healthy controls and 143 patients with schizophrenia. After harmonising the raw dMRIs from three different studies, we performed whole-brain two-tensor tractography and fibre clustering on the tractography data. We compared the fractional anisotropy (FA) of white matter tracts between healthy controls and patients with schizophrenia. Spearman's rho was adopted for the associations with clinical symptoms measured by the Positive and Negative Syndrome Scale (PANSS). The Bonferroni correction was used to adjust multiple testing.

RESULTS

Among the 33 DWM and 8 SWM tracts, patients with schizophrenia had a lower FA in 14 DWM and 4 SWM tracts than healthy controls, with small effect sizes. In the patient group, the FA deviations of the corticospinal and superficial-occipital tracts were negatively correlated with the PANSS negative score; however, this correlation was not evident after adjusting for multiple testing.

CONCLUSION

We observed the structural impairments of both the DWM and SWM tracts in patients with schizophrenia. The SWM could be a potential target of interest in future research on neural biomarkers for schizophrenia.

Ancestral, Pregnancy, and Negative Early-Life Risks Shape Children's Brain (Dis)similarity to Schizophrenia

BACKGROUND

Familial, obstetric, and early-life environmental risks for schizophrenia spectrum disorder (SSD) alter normal cerebral development, leading to the formation of characteristic brain deficit patterns prior to onset of symptoms. We hypothesized that the insidious effects of these risks may increase brain similarity to adult SSD deficit patterns in prepubescent children.

METHODS

We used data collected by the Adolescent Brain Cognitive Development (ABCD) Study (N = 8940, age = 9.9 ± 0.1 years, 4307/4633 female/male), including 727 (age = 9.9 ± 0.1 years, 351/376 female/male) children with family history of SSD, to evaluate unfavorable cerebral effects of ancestral SSD history, pre/perinatal environment, and negative early-life environment. We used a regional vulnerability index to measure the alignment of a child's cerebral patterns with the adult SSD pattern derived from a large meta-analysis of case-control differences.

RESULTS

In children with a family history of SSD, the regional vulnerability index captured significantly more variance in ancestral history than traditional whole-brain and regional brain measurements. In children with and without family history of SSD, the regional vulnerability index also captured more variance associated with negative pre/perinatal environment and early-life experiences than traditional brain measurements.

CONCLUSIONS

In summary, in a cohort in which most children will not develop SSD, familial, pre/perinatal, and early developmental risks can alter brain patterns in the direction observed in adult patients with SSD. Individual similarity to adult SSD patterns may provide an early biomarker of the effects of genetic and developmental risks on the brain prior to psychotic or prodromal symptom onset.

Heart-brain connections: Phenotypic and genetic insights from magnetic resonance images

Cardiovascular health interacts with cognitive and mental health in complex ways, yet little is known about the phenotypic and genetic links of heart-brain systems. We quantified heart-brain connections using multiorgan magnetic resonance imaging (MRI) data from more than 40,000 subjects. Heart MRI traits displayed numerous association patterns with brain gray matter morphometry, white matter microstructure, and functional networks. We identified 80 associated genomic loci (P < 6.09 × 10-10) for heart MRI traits, which shared genetic influences with cardiovascular and brain diseases. Genetic correlations were observed between heart MRI traits and brain-related traits and disorders. Mendelian randomization suggests that heart conditions may causally contribute to brain disorders. Our results advance a multiorgan perspective on human health by revealing heart-brain connections and shared genetic influences.

Blockchain-enabled immutable, distributed, and highly available clinical research activity logging system for federated COVID-19 data analysis from multiple institutions

OBJECTIVE

We aimed to develop a distributed, immutable, and highly available cross-cloud blockchain system to facilitate federated data analysis activities among multiple institutions.

MATERIALS AND METHODS

We preprocessed 9166 COVID-19 Structured Query Language (SQL) code, summary statistics, and user activity logs, from the GitHub repository of the Reliable Response Data Discovery for COVID-19 (R2D2) Consortium. The repository collected local summary statistics from participating institutions and aggregated the global result to a COVID-19-related clinical query, previously posted by clinicians on a website. We developed both on-chain and off-chain components to store/query these activity logs and their associated queries/results on a blockchain for immutability, transparency, and high availability of research communication. We measured run-time efficiency of contract deployment, network transactions, and confirmed the accuracy of recorded logs compared to a centralized baseline solution.

RESULTS

The smart contract deployment took 4.5 s on an average. The time to record an activity log on blockchain was slightly over 2 s, versus 5-9 s for baseline. For querying, each query took on an average less than 0.4 s on blockchain, versus around 2.1 s for baseline.

DISCUSSION

The low deployment, recording, and querying times confirm the feasibility of our cross-cloud, blockchain-based federated data analysis system. We have yet to evaluate the system on a larger network with multiple nodes per cloud, to consider how to accommodate a surge in activities, and to investigate methods to lower querying time as the blockchain grows.

CONCLUSION

Blockchain technology can be used to support federated data analysis among multiple institutions.

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) study has collected data from over 10,000 children across 21 sites, providing valuable insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a database of harmonized and processed ABCD dMRI data has been created, comprising quality-controlled imaging data from 9345 subjects. This resource required significant computational effort, taking ~50,000 CPU hours to harmonize the data, perform white matter parcellation, and run whole brain tractography. The database includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts both in full-resolution (for analysis) and low-resolution (for visualization), and 804 different dMRI-derived measures per subject. It is available via the NIMH Data Archive and offers tremendous potential for scientific discoveries in structural connectivity studies of neurodevelopment in children and adolescents. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Sex-specific differences in the genetic and environmental effects on cardiac phenotypic variation assessed by echocardiography

The drivers of sexual dimorphism in heart failure phenotypes are currently poorly understood. Divergent phenotypes may result from differences in heritability and genetic versus environmental influences on the interplay of cardiac structure and function. To assess sex-specific heritability and genetic versus environmental contributions to variation and inter-relations between echocardiography traits in a large community-based cohort. We studied Framingham Heart Study participants of Offspring Cohort examination 8 (2005-2008) and Third Generation Cohort examination 1 (2002-2005). Five cardiac traits and six functional traits were measured using standardized echocardiography. Sequential Oligogenic Linkage Analysis Routines (SOLAR) software was used to perform singular and bivariate quantitative trait linkage analysis. In our study of 5674 participants (age 49 ± 15 years; 54% women), heritability for all traits was significant for both men and women. There were no significant differences in traits between men and women. Within inter-trait correlations, there were two genetic, and four environmental trait pairs with sex-based differences. Within both significant genetic trait pairs, men had a positive relation, and women had no significant relation. We observed significant sex-based differences in inter-trait genetic and environmental correlations between cardiac structure and function. These findings highlight potential pathways of sex-based divergent heart failure phenotypes.

Eye-brain connections revealed by multimodal retinal and brain imaging genetics in the UK Biobank

As an anatomical extension of the brain, the retina of the eye is synaptically connected to the visual cortex, establishing physiological connections between the eye and the brain. Despite the unique opportunity retinal structures offer for assessing brain disorders, less is known about their relationship to brain structure and function. Here we present a systematic cross-organ genetic architecture analysis of eye-brain connections using retina and brain imaging endophenotypes. Novel phenotypic and genetic links were identified between retinal imaging biomarkers and brain structure and function measures derived from multimodal magnetic resonance imaging (MRI), many of which were involved in the visual pathways, including the primary visual cortex. In 65 genomic regions, retinal imaging biomarkers shared genetic influences with brain diseases and complex traits, 18 showing more genetic overlaps with brain MRI traits. Mendelian randomization suggests that retinal structures have bidirectional genetic causal links with neurological and neuropsychiatric disorders, such as Alzheimer's disease. Overall, cross-organ imaging genetics reveals a genetic basis for eye-brain connections, suggesting that the retinal images can elucidate genetic risk factors for brain disorders and disease-related changes in intracranial structure and function.

Functional hybrid factor regression model for handling heterogeneity in imaging studies

This paper develops a functional hybrid factor regression modelling framework to handle the heterogeneity of many large-scale imaging studies, such as the Alzheimer's disease neuroimaging initiative study. Despite the numerous successes of those imaging studies, such heterogeneity may be caused by the differences in study environment, population, design, protocols or other hidden factors, and it has posed major challenges in integrative analysis of imaging data collected from multicentres or multistudies. We propose both estimation and inference procedures for estimating unknown parameters and detecting unknown factors under our new model. The asymptotic properties of both estimation and inference procedures are systematically investigated. The finite-sample performance of our proposed procedures is assessed by using Monte Carlo simulations and a real data example on hippocampal surface data from the Alzheimer's disease study.

Comparison of Diffusion Tensor Imaging Metrics in Normal-Appearing White Matter to Cerebrovascular Lesions and Correlation with Cerebrovascular Disease Risk Factors and Severity

Alterations in tissue microstructure in normal-appearing white matter (NAWM), specifically measured by diffusion tensor imaging (DTI) fractional anisotropy (FA), have been associated with cognitive outcomes following stroke. The purpose of this study was to comprehensively compare conventional DTI measures of tissue microstructure in NAWM to diverse vascular brain lesions in people with cerebrovascular disease (CVD) and to examine associations between FA in NAWM and cerebrovascular risk factors. DTI metrics including fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were measured in cerebral tissues and cerebrovascular anomalies from 152 people with CVD participating in the Ontario Neurodegenerative Disease Research Initiative (ONDRI). Ten cerebral tissue types were segmented including NAWM, and vascular lesions including stroke, periventricular and deep white matter hyperintensities, periventricular and deep lacunar infarcts, and perivascular spaces (PVS) using T1-weighted, proton density-weighted, T2-weighted, and fluid attenuated inversion recovery MRI scans. Mean DTI metrics were measured in each tissue region using a previously developed DTI processing pipeline and compared between tissues using multivariate analysis of covariance. Associations between FA in NAWM and several CVD risk factors were also examined. DTI metrics in vascular lesions differed significantly from healthy tissue. Specifically, all tissue types had significantly different MD values, while FA was also found to be different in most tissue types. FA in NAWM was inversely related to hypertension and modified Rankin scale (mRS). This study demonstrated the differences between conventional DTI metrics, FA, MD, AD, and RD, in cerebral vascular lesions and healthy tissue types. Therefore, incorporating DTI to characterize the integrity of the tissue microstructure could help to define the extent and severity of various brain vascular anomalies. The association between FA within NAWM and clinical evaluation of hypertension and disability provides further evidence that white matter microstructural integrity is impacted by cerebrovascular function.

dsMTL - a computational framework for privacy-preserving, distributed multi-task machine learning

MOTIVATION

In multi-cohort machine learning studies, it is critical to differentiate between effects that are reproducible across cohorts and those that are cohort-specific. Multi-task learning (MTL) is a machine learning approach that facilitates this differentiation through the simultaneous learning of prediction tasks across cohorts. Since multi-cohort data can often not be combined into a single storage solution, there would be the substantial utility of an MTL application for geographically distributed data sources.

RESULTS

Here, we describe the development of "dsMTL", a computational framework for privacy-preserving, distributed multi-task machine learning that includes three supervised and one unsupervised algorithms. First, we derive the theoretical properties of these methods and the relevant machine learning workflows to ensure the validity of the software implementation. Second, we implement dsMTL as a library for the R programming language, building on the DataSHIELD platform that supports the federated analysis of sensitive individual-level data. Third, we demonstrate the applicability of dsMTL for comorbidity modeling in distributed data. We show that comorbidity modeling using dsMTL outperformed conventional, federated machine learning, as well as the aggregation of multiple models built on the distributed datasets individually. The application of dsMTL was computationally efficient and highly scalable when applied to moderate-size (n < 500), real expression data given the actual network latency.

AVAILABILITY

dsMTL is freely available at https://github.com/transbioZI/dsMTLBase (server-side package) and https://github.com/transbioZI/dsMTLClient (client-side package).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Implicating Causal Brain Magnetic Resonance Imaging in Glaucoma Using Mendelian Randomization

Background

Glaucoma is hypothesized to originate in the brain but manifests as an eye disease as it possesses the common features of neurodegeneration diseases. But there is no evidence to demonstrate the primary brain changes in glaucoma patients. In the present study, we have used Mendelian randomization (MR) to understand the causal effect of brain alterations on glaucoma.

Methods

Our MR study was carried out using summary statistics from genome-wide associations for 110 diffusion tensor imaging (DTI) measurements of white matter (WM) tracts (17,706 individuals), 101 brain region-of-interest (ROI) volumes (19,629 individuals), and glaucoma (8,591 cases, 210,201 control subjects). The causal relationship was evaluated by multiplicative random effects inverse variance weighted (IVW) method and verified by two other MR methods, including MR Egger, weighted median, and extensive sensitivity analyses.

Results

Genetic liability to fornix fractional anisotropy (FX.FA) (OR = 0.71, 95%CI = 0.56–0.88, P = 2.44 × 10^–3), and uncinate fasciculus UNC.FA (OR = 0.65, 95%CI = 0.48–0.88, P = 5.57 × 10^–3) was associated with a low risk of glaucoma. Besides, the right ventral diencephalon (OR = 1.72, 95%CI = 1.17–2.52, P = 5.64 × 10^–3) and brain stem (OR = 1.35, 95%CI = 1.08–1.69, P = 8.94 × 10^–3) were associated with the increased risk of glaucoma. No heterogeneity and pleiotropy were detected.

Conclusion

Our study suggests that the fornix and uncinate fasciculus degenerations and injures of the right ventral diencephalon and brain stem potentially increase the occurrence of glaucoma and reveal the existence of the brain-eye axis.

ErbB Signaling Pathway Genes Are Differentially Expressed in Monozygotic Twins Discordant for Sports-Related Concussion

Transcriptional changes involved in neuronal recovery after sports-related concussion (SRC) may be obscured by inter-individual variation in mRNA expression and nonspecific changes related to physical exertion. Using a co-twin study, the objective of this study was to identify important differences in mRNA expression among a single pair of monozygotic (MZ) twins discordant for concussion. A pair of MZ twins were enrolled as part of a larger study of concussion biomarkers among collegiate athletes. During the study, Twin A sustained SRC, allowing comparison of mRNA expression to the nonconcussed Twin B. Twin A clinically recovered by Day 7. mRNA expression was measured pre-injury and at 6 h and 7 days postinjury using Affymetrix HG-U133 Plus 2.0 microarray. Changes in mRNA expression from pre-injury to each postinjury time point were compared between the twins; differences >1.5-fold were considered important. Kyoto Encyclopedia of Genes and Genomes identified biologic networks associated with important transcripts. Among 38,000 analyzed genes, important changes were identified in 153 genes. The ErbB (epidermal growth factor receptor) signaling pathway was identified as the top transcriptional network from pre-injury to 7 days postinjury. Genes in this pathway with important transcriptional changes included epidermal growth factor (2.41), epiregulin (1.73), neuregulin 1 (1.54) and mechanistic target of rapamycin (1.51). In conclusion, the ErbB signaling pathway was identified as a potential regulator of clinical recovery in a MZ twin pair discordant for SRC. A co-twin study design may be a useful method for identifying important gene pathways associated with concussion recovery.

MAGNIMS recommendations for harmonization of MRI data in MS multicenter studies

•

Sharing data from cooperative studies is essential to develop new biomarkers in MS.

•

Differences in MRI acquisition, analysis, storage represent a substantial constraint.

•

We review the state of the art and developments in the harmonization of MRI.

•

We provide recommendations to harmonize large MRI datasets in the MS field.

There is an increasing need of sharing harmonized data from large, cooperative studies as this is essential to develop new diagnostic and prognostic biomarkers. In the field of multiple sclerosis (MS), the issue has become of paramount importance due to the need to translate into the clinical setting some of the most recent MRI achievements. However, differences in MRI acquisition parameters, image analysis and data storage across sites, with their potential bias, represent a substantial constraint. This review focuses on the state of the art, recent technical advances, and desirable future developments of the harmonization of acquisition, analysis and storage of large-scale multicentre MRI data of MS cohorts. Huge efforts are currently being made to achieve all the requirements needed to provide harmonized MRI datasets in the MS field, as proper management of large imaging datasets is one of our greatest opportunities and challenges in the coming years. Recommendations based on these achievements will be provided here. Despite the advances that have been made, the complexity of these tasks requires further research by specialized academical centres, with dedicated technical and human resources. Such collective efforts involving different professional figures are of crucial importance to offer to MS patients a personalised management while minimizing consumption of resources.

Longitudinal assessment of early-life white matter development with quantitative relaxometry in nonhuman primates

Alterations in white matter (WM) development are associated with many neuropsychiatric and neurodevelopmental disorders. Most MRI studies examining WM development employ diffusion tensor imaging (DTI), which relies on estimating diffusion patterns of water molecules as a reflection of WM microstructure. Quantitative relaxometry, an alternative method for characterizing WM microstructural changes, is based on molecular interactions associated with the magnetic relaxation of protons. In a longitudinal study of 34 infant non-human primates (NHP) (Macaca mulatta) across the first year of life, we implement a novel, high-resolution, T1-weighted MPnRAGE sequence to examine WM trajectories of the longitudinal relaxation rate (qR1) in relation to DTI metrics and gestational age at scan. To the best of our knowledge, this is the first study to assess developmental WM trajectories in NHPs using quantitative relaxometry and the first to directly compare DTI and relaxometry metrics during infancy. We demonstrate that qR1 exhibits robust logarithmic growth, unfolding in a posterior-anterior and medial-lateral fashion, similar to DTI metrics. On a within-subject level, DTI metrics and qR1 are highly correlated, but are largely unrelated on a between-subject level. Unlike DTI metrics, gestational age at birth (time in utero) is a strong predictor of early postnatal qR1 levels. Whereas individual differences in DTI metrics are maintained across the first year of life, this is not the case for qR1. These results point to the similarities and differences in using quantitative relaxometry and DTI in developmental studies, providing a basis for future studies to characterize the unique processes that these measures reflect at the cellular and molecular level.

Common variants contribute to intrinsic human brain functional networks

The human brain forms functional networks of correlated activity, which have been linked with both cognitive and clinical outcomes. However, the genetic variants affecting brain function are largely unknown. Here, we used resting-state functional magnetic resonance images from 47,276 individuals to discover and validate common genetic variants influencing intrinsic brain activity. We identified 45 new genetic regions associated with brain functional signatures (P < 2.8 × 10^-11), including associations to the central executive, default mode, and salience networks involved in the triple-network model of psychopathology. A number of brain activity-associated loci colocalized with brain disorders (e.g., the APOE ε4 locus with Alzheimer's disease). Variation in brain function was genetically correlated with brain disorders, such as major depressive disorder and schizophrenia. Together, our study provides a step forward in understanding the genetic architecture of brain functional networks and their genetic links to brain-related complex traits and disorders.

ENIGMA-DTI: Translating reproducible white matter deficits into personalized vulnerability metrics in cross-diagnostic psychiatric research

The ENIGMA-DTI (diffusion tensor imaging) workgroup supports analyses that examine the effects of psychiatric, neurological, and developmental disorders on the white matter pathways of the human brain, as well as the effects of normal variation and its genetic associations. The seven ENIGMA disorder-oriented working groups used the ENIGMA-DTI workflow to derive patterns of deficits using coherent and coordinated analyses that model the disease effects across cohorts worldwide. This yielded the largest studies detailing patterns of white matter deficits in schizophrenia spectrum disorder (SSD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), posttraumatic stress disorder (PTSD), traumatic brain injury (TBI), and 22q11 deletion syndrome. These deficit patterns are informative of the underlying neurobiology and reproducible in independent cohorts. We reviewed these findings, demonstrated their reproducibility in independent cohorts, and compared the deficit patterns across illnesses. We discussed translating ENIGMA-defined deficit patterns on the level of individual subjects using a metric called the regional vulnerability index (RVI), a correlation of an individual's brain metrics with the expected pattern for a disorder. We discussed the similarity in white matter deficit patterns among SSD, BD, MDD, and OCD and provided a rationale for using this index in cross-diagnostic neuropsychiatric research. We also discussed the difference in deficit patterns between idiopathic schizophrenia and 22q11 deletion syndrome, which is used as a developmental and genetic model of schizophrenia. Together, these findings highlight the importance of collaborative large-scale research to provide robust and reproducible effects that offer insights into individual vulnerability and cross-diagnosis features.

Mega-analysis methods in ENIGMA: The experience of the generalized anxiety disorder working group

The ENIGMA group on Generalized Anxiety Disorder (ENIGMA-Anxiety/GAD) is part of a broader effort to investigate anxiety disorders using imaging and genetic data across multiple sites worldwide. The group is actively conducting a mega-analysis of a large number of brain structural scans. In this process, the group was confronted with many methodological challenges related to study planning and implementation, between-country transfer of subject-level data, quality control of a considerable amount of imaging data, and choices related to statistical methods and efficient use of resources. This report summarizes the background information and rationale for the various methodological decisions, as well as the approach taken to implement them. The goal is to document the approach and help guide other research groups working with large brain imaging data sets as they develop their own analytic pipelines for mega-analyses.

White matter impairments in patients with schizophrenia: A multisite diffusion MRI study

There is a lack of convincing and replicative findings regarding white matter abnormalities in schizophrenia. Several multisite diffusion magnetic resonance imaging (dMRI) studies have been conducted to increase statistical power and reveal subtle white matter changes. Data pooling methods are crucial in joint analysis to compensate for the use of different scanners and image acquisition parameters. A harmonization method using raw dMRI data was developed to overcome the limited generalizability of previous data pooling methods. We obtained dMRI data of 242 healthy controls and 190 patients with schizophrenia from four different study sites. After applying the harmonization method to the raw dMRI data, a two-tensor whole-brain tractography was performed, and diffusion measures were compared between the two groups. The correlation of fractional anisotropy (FA) with the positive and negative symptoms was evaluated, and the interaction effect of diagnosis-by-age, age-squared, and sex was examined. The following white matter tracts showed significant group differences in the FA: the right superior longitudinal fascicle (SLF), the left-to-right lateral orbitofrontal commissural tract, pars orbitalis (pOr-pOr) commissural tract, and pars triangularis (pTr-pTr) commissural tract. The FA of the right SLF and pTr-pTr commissural tract were significantly associated with the Positive and Negative Syndrome Scale (PANSS) positive and negative scores. No significant interaction effect was observed. These findings add to the evidence on structural brain abnormalities in schizophrenia and can aid in obtaining a better understanding of the biological foundations of schizophrenia.

Comparing empirical kinship derived heritability for imaging genetics traits in the UK biobank and human connectome project

Imaging genetics analyses use neuroimaging traits as intermediate phenotypes to infer the degree of genetic contribution to brain structure and function in health and/or illness. Coefficients of relatedness (CR) summarize the degree of genetic similarity among subjects and are used to estimate the heritability – the proportion of phenotypic variance explained by genetic factors. The CR can be inferred directly from genome-wide genotype data to explain the degree of shared variation in common genetic polymorphisms (SNP-heritability) among related or unrelated subjects. We developed a central processing and graphics processing unit (CPU and GPU) accelerated Fast and Powerful Heritability Inference (FPHI) approach that linearizes likelihood calculations to overcome the ~N^2–3 computational effort dependency on sample size of classical likelihood approaches. We calculated for 60 regional and 1.3 × 10⁵ voxel-wise traits in N = 1,206 twin and sibling participants from the Human Connectome Project (HCP) (550 M/656 F, age = 28.8 ± 3.7 years) and N = 37,432 (17,531 M/19,901 F; age = 63.7 ± 7.5 years) participants from the UK Biobank (UKBB). The FPHI estimates were in excellent agreement with heritability values calculated using Genome-wide Complex Trait Analysis software (r = 0.96 and 0.98 in HCP and UKBB sample) while significantly reducing computational (10^2–4 times). The regional and voxel-wise traits heritability estimates for the HCP and UKBB were likewise in excellent agreement (r = 0.63–0.76, p < 10⁻¹⁰). In summary, the hardware-accelerated FPHI made it practical to calculate heritability values for voxel-wise neuroimaging traits, even in very large samples such as the UKBB. The patterns of additive genetic variance in neuroimaging traits measured in a large sample of related and unrelated individuals showed excellent agreement regardless of the estimation method. The code and instruction to execute these analyses are available at www.solar-eclipse-genetics.org.

Realistic image normalization for multi-Domain segmentation

Image normalization is a building block in medical image analysis. Conventional approaches are customarily employed on a per-dataset basis. This strategy, however, prevents the current normalization algorithms from fully exploiting the complex joint information available across multiple datasets. Consequently, ignoring such joint information has a direct impact on the processing of segmentation algorithms. This paper proposes to revisit the conventional image normalization approach by, instead, learning a common normalizing function across multiple datasets. Jointly normalizing multiple datasets is shown to yield consistent normalized images as well as an improved image segmentation when intensity shifts are large. To do so, a fully automated adversarial and task-driven normalization approach is employed as it facilitates the training of realistic and interpretable images while keeping performance on par with the state-of-the-art. The adversarial training of our network aims at finding the optimal transfer function to improve both, jointly, the segmentation accuracy and the generation of realistic images. We have evaluated the performance of our normalizer on both infant and adult brain images from the iSEG, MRBrainS and ABIDE datasets. The results indicate that our contribution does provide an improved realism to the normalized images, while retaining a segmentation accuracy at par with the state-of-the-art learnable normalization approaches.

Maternal serotonin transporter genotype and offsprings' clinical and cognitive measures of ADHD and ASD

Serotonin (5-HT) is an important factor for prenatal neurodevelopment whereby its neurotrophic actions can be regulated through maternal-fetal interactions. We explored if maternal 5-HTTLPR genotype is associated with clinical and cognitive measures of attention-deficit/hyperactivity disorder (ADHD) and comorbid autism spectrum disorder (ASD) in typically-developing and ADHD-diagnosed offspring, beyond classical inheritance and environmental- and comorbidity-mediators/confounders. Family-based variance decomposition analyses were performed incorporating 6-31 year-old offsprings' as well as parental genotypes of 462 ADHD and control families from the NeuroIMAGE cohort. Dependent measures were offsprings' ADHD symptom- and ASD trait-scores and cognitive measures including executive functioning (including response inhibition and cognitive flexibility), sustained attention, reward processing, motor control, and emotion recognition. Offsprings' stereotyped behavior was predicted by an interaction between maternal 5-HTTLPR genotype and offsprings' sex. Furthermore, offspring of mothers with low-expressing genotypes demonstrated larger reward-related reductions in reaction time. While specifically adult male offspring of these mothers reported a faster reversal learning with less errors, specifically young female offspring of these mothers were more accurate in identifying happy faces. Adult offspring from the mothers with low-expressing 5-HTTLPR genotypes were also slower in identifying happy faces. However, this association seemed to be mediated by offsprings' high anxiety levels. In sum, we found some support for a role of the maternal 5-HT system in modulating fetal brain development and behavior. Offsprings' cognitive measures might be more sensitive to small alterations within the maternal 5-HT system than their ADHD and ASD clinical phenotypes. Further studies are needed to specify the association between maternal genotype and risk for neurodevelopmental disorders.

Altered white matter microstructural organization in posttraumatic stress disorder across 3047 adults: results from the PGC-ENIGMA PTSD consortium

A growing number of studies have examined alterations in white matter organization in people with posttraumatic stress disorder (PTSD) using diffusion MRI (dMRI), but the results have been mixed which may be partially due to relatively small sample sizes among studies. Altered structural connectivity may be both a neurobiological vulnerability for, and a result of, PTSD. In an effort to find reliable effects, we present a multi-cohort analysis of dMRI metrics across 3047 individuals from 28 cohorts currently participating in the PGC-ENIGMA PTSD working group (a joint partnership between the Psychiatric Genomics Consortium and the Enhancing NeuroImaging Genetics through Meta-Analysis consortium). Comparing regional white matter metrics across the full brain in 1426 individuals with PTSD and 1621 controls (2174 males/873 females) between ages 18-83, 92% of whom were trauma-exposed, we report associations between PTSD and disrupted white matter organization measured by lower fractional anisotropy (FA) in the tapetum region of the corpus callosum (Cohen's d = -0.11, p = 0.0055). The tapetum connects the left and right hippocampus, for which structure and function have been consistently implicated in PTSD. Results were consistent even after accounting for the effects of multiple potentially confounding variables: childhood trauma exposure, comorbid depression, history of traumatic brain injury, current alcohol abuse or dependence, and current use of psychotropic medications. Our results show that PTSD may be associated with alterations in the broader hippocampal network.

Large-scale GWAS reveals genetic architecture of brain white matter microstructure and genetic overlap with cognitive and mental health traits (n = 17,706)

Individual variations of white matter (WM) tracts are known to be associated with various cognitive and neuropsychiatric traits. Diffusion tensor imaging (DTI) and genome-wide single-nucleotide polymorphism (SNP) data from 17,706 UK Biobank participants offer the opportunity to identify novel genetic variants of WM tracts and explore the genetic overlap with other brain-related complex traits. We analyzed the genetic architecture of 110 tract-based DTI parameters, carried out genome-wide association studies (GWAS), and performed post-GWAS analyses, including association lookups, gene-based association analysis, functional gene mapping, and genetic correlation estimation. We found that DTI parameters are substantially heritable for all WM tracts (mean heritability 48.7%). We observed a highly polygenic architecture of genetic influence across the genome (p value = 1.67 × 10^-05) as well as the enrichment of genetic effects for active SNPs annotated by central nervous system cells (p value = 8.95 × 10^-12). GWAS identified 213 independent significant SNPs associated with 90 DTI parameters (696 SNP-level and 205 locus-level associations; p value < 4.5 × 10^-10, adjusted for testing multiple phenotypes). Gene-based association study prioritized 112 significant genes, most of which are novel. More importantly, association lookups found that many of the novel SNPs and genes of DTI parameters have previously been implicated with cognitive and mental health traits. In conclusion, the present study identifies many new genetic variants at SNP, locus and gene levels for integrity of brain WM tracts and provides the overview of pleiotropy with cognitive and mental health traits.

Common genetic variation influencing human white matter microstructure

Brain regions communicate with each other through tracts of myelinated axons, commonly referred to as white matter. We identified common genetic variants influencing white matter microstructure using diffusion magnetic resonance imaging of 43,802 individuals. Genome-wide association analysis identified 109 associated loci, 30 of which were detected by tract-specific functional principal components analysis. A number of loci colocalized with brain diseases, such as glioma and stroke. Genetic correlations were observed between white matter microstructure and 57 complex traits and diseases. Common variants associated with white matter microstructure altered the function of regulatory elements in glial cells, particularly oligodendrocytes. This large-scale tract-specific study advances the understanding of the genetic architecture of white matter and its genetic links to a wide spectrum of clinical outcomes.

Abnormal hemispheric asymmetry of both brain function and structure in attention deficit/hyperactivity disorder: a meta-analysis of individual participant data

Aberration in the asymmetric nature of the human brain is associated with several mental disorders, including attention deficit/hyperactivity disorder (ADHD). In ADHD, these aberrations are thought to reflect key hemispheric differences in the functioning of attention, although the structural and functional bases of these defects are yet to be fully characterized. In this study, we applied a comprehensive meta-analysis to multimodal imaging datasets from 627 subjects (303 typically developing control [TDCs] and 324 patients with ADHD) with both resting-state functional and structural magnetic resonance imaging (MRI), from seven independent publicly available datasets of the ADHD-200 sample. We performed lateralization analysis and calculated the combined effects of ADHD on each of three cortical regional measures (grey matter volume - GMV, fractional amplitude of low frequency fluctuations at rest -fALFF, and regional homogeneity -ReHo). We found that compared with TDC, 68%,73% and 66% of regions showed statistically significant ADHD disorder effects on the asymmetry of GMV, fALFF, and ReHo, respectively, (false discovery rate corrected, q = 0.05). Forty-one percent (41%) of regions had both structural and functional abnormalities in asymmetry, located in the prefrontal, frontal, and subcortical cortices, and the cerebellum. Furthermore, brain asymmetry indices in these regions were higher in children with more severe ADHD symptoms, indicating a crucial pathoplastic role for asymmetry. Our findings highlight the functional asymmetry in ADHD which has (1) a strong structural basis, and thus is likely to be developmental in nature; and (2) is strongly linked to symptom burden and IQ and may carry a possible prognostic value for grading the severity of ADHD.

Evidence of Genetic Overlap Between Circadian Preference and Brain White Matter Microstructure

Several neuroimaging studies have reported associations between brain white matter microstructure and chronotype. However, it is unclear whether those phenotypic relationships are causal or underlined by genetic factors. In the present study, we use genetic data to examine the genetic overlap and infer causal relationships between chronotype and diffusion tensor imaging (DTI) measures. We identify 29 significant pairwise genetic correlations, of which 13 also show evidence for a causal association. Genetic correlations were identified between chronotype and brain-wide mean, axial and radial diffusivities. When exploring individual tracts, 10 genetic correlations were observed with mean diffusivity, 10 with axial diffusivity, 4 with radial diffusivity and 2 with mode of anisotropy. We found evidence for a possible causal association of eveningness with white matter microstructure measures in individual tracts including the posterior limb and the retrolenticular part of the internal capsule; the genu and splenium of the corpus callosum and the posterior, superior and anterior regions of the corona radiata. Our findings contribute to the understanding of how genes influence circadian preference and brain white matter and provide a new avenue for investigating the role of chronotype in health and disease.

Genome-Wide Association Analysis of Neonatal White Matter Microstructure

A better understanding of genetic influences on early white matter development could significantly advance our understanding of neurological and psychiatric conditions characterized by altered integrity of axonal pathways. We conducted a genome-wide association study (GWAS) of diffusion tensor imaging (DTI) phenotypes in 471 neonates. We used a hierarchical functional principal regression model (HFPRM) to perform joint analysis of 44 fiber bundles. HFPRM revealed a latent measure of white matter microstructure that explained approximately 50% of variation in our tractography-based measures and accounted for a large proportion of heritable variation in each individual bundle. An intronic SNP in PSMF1 on chromosome 20 exceeded the conventional GWAS threshold of 5 x 10-8 (p = 4.61 x 10-8). Additional loci nearing genome-wide significance were located near genes with known roles in axon growth and guidance, fasciculation, and myelination.

[The relationship between inflammation, cognitive disorders and neuroimaging data in schizophrenia]

OBJECTIVE

To search for the relationship between the results of functional imaging, immunological parameters and laboratory markers of inflammation in schizophrenia, taking into account cognitive impairment in patients, and to consider the possibility of using a multidisciplinary approach to diagnosis, treatment and prognosis of schizophrenia.

MATERIAL AND METHODS

The study included 25 patients with schizophrenia and 13 healthy volunteers. Psychiatric scales were administered to evaluate the patient's condition. The main indicators of humoral immunity, the level of markers of inflammation, key pro-inflammatory and anti-inflammatory cytokines, and growth factor VEGF were determined by ELISA. Brain MRI was performed. All calculated tractographic data are included in the connection database to study the effect of immunological markers and the degree of severity of cognitive impairment.

RESULTS AND CONCLUSION

Levels of markers of systemic inflammation and growth factor VEGF-A as well as the activation of humoral immunity are increased in patients with schizophrenia compared with controls. For the first time, the relationship of immunological parameters with the coefficient of quantitative anisotropy in the area of the corpus callosum in schizophrenia was revealed. The results indicate the possible value of indicators of the activation of the humoral immune response and systemic inflammation as markers of neurophysiological changes and cognitive dysfunction in schizophrenia.

Genes associated with gray matter volume alterations in schizophrenia

Although both schizophrenia and gray matter volume (GMV) show high heritability, however, genes accounting for GMV alterations in schizophrenia remain largely unknown. Based on risk genes identified in schizophrenia by the genome-wide association study of the Schizophrenia Working Group of the Psychiatric Genomics Consortium, we used transcription-neuroimaging association analysis to test that which of these genes are associated with GMV changes in schizophrenia. For each brain tissue sample, the expression profiles of 196 schizophrenia risk genes were extracted from six donated normal brains of the Allen Human Brain Atlas, and GMV differences between patients with schizophrenia and healthy controls were calculated based on five independent case-control structural MRI datasets (276 patients and 284 controls). Genes associated with GMV changes in schizophrenia were identified by performing cross-sample spatial correlations between expression levels of each gene and case-control GMV difference derived from the five MRI datasets integrated by harmonization and meta-analysis. We found that expression levels of 98 genes consistently showed significant cross-sample spatial correlations with GMV changes in schizophrenia. These genes were functionally enriched for chemical synaptic transmission, central nervous system development, and cell projection. Overall, this study provides a set of genes possibly associated with GMV changes in schizophrenia, which could be used as candidate genes to explore biological mechanisms underlying the structural impairments in schizophrenia.

The reliability and heritability of cortical folds and their genetic correlations across hemispheres

Cortical folds help drive the parcellation of the human cortex into functionally specific regions. Variations in the length, depth, width, and surface area of these sulcal landmarks have been associated with disease, and may be genetically mediated. Before estimating the heritability of sulcal variation, the extent to which these metrics can be reliably extracted from in-vivo MRI must be established. Using four independent test-retest datasets, we found high reliability across the brain (intraclass correlation interquartile range: 0.65-0.85). Heritability estimates were derived for three family-based cohorts using variance components analysis and pooled (total N > 3000); the overall sulcal heritability pattern was correlated to that derived for a large population cohort (N > 9000) calculated using genomic complex trait analysis. Overall, sulcal width was the most heritable metric, and earlier forming sulci showed higher heritability. The inter-hemispheric genetic correlations were high, yet select sulci showed incomplete pleiotropy, suggesting hemisphere-specific genetic influences.

Interhemispheric Relationship of Genetic Influence on Human Brain Connectivity

To understand the origins of interhemispheric differences and commonalities/coupling in human brain wiring, it is crucial to determine how homologous interregional connectivities of the left and right hemispheres are genetically determined and related. To address this, in the present study, we analyzed human twin and pedigree samples with high-quality diffusion magnetic resonance imaging tractography and estimated the heritability and genetic correlation of homologous left and right white matter (WM) connections. The results showed that the heritability of WM connectivity was similar and coupled between the 2 hemispheres and that the degree of overlap in genetic factors underlying homologous WM connectivity (i.e., interhemispheric genetic correlation) varied substantially across the human brain: from complete overlap to complete nonoverlap. Particularly, the heritability was significantly stronger and the chance of interhemispheric complete overlap in genetic factors was higher in subcortical WM connections than in cortical WM connections. In addition, the heritability and interhemispheric genetic correlations were stronger for long-range connections than for short-range connections. These findings highlight the determinants of the genetics underlying WM connectivity and its interhemispheric relationships, and provide insight into genetic basis of WM connectivity asymmetries in both healthy and disease states.

Effect of DAOA genetic variation on white matter alteration in corpus callosum in patients with first-episode schizophrenia

D-amino acid oxidase activator (DAOA) gene, which plays a crucial role in the process of glutamatergic transmission and mitochondrial function, is frequently linked with the liability for schizophrenia. We aimed to investigate whether the variation of DAOA rs2391191 is associated with alterations in white matter integrity of first-episode schizophrenia (FES) patients; and whether it influences the association between white matter integrity, cognitive function and clinical symptoms of schizophrenia. Forty-six patients with FES and forty-nine healthy controls underwent DTI and were genotyped for DAOA rs2391191. Psychopathological assessments were performed by Brief Psychiatric Rating Scale (BPRS) and Scale for Assessment of Negative Symptoms (SANS). Cognitive function was assessed by MATRICS Consensus Cognitive Battery (MCCB). Schizophrenia patients presented lower fractional anisotropy (FA) and higher radial diffusivity (RD), mainly spreading over the corpus callosum and corona radiata compared with healthy controls. Compared with patients carrying G allele, patients with AA showed lower FA in the body of corpus callosum, and higher RD in the genu of corpus callosum, right superior and anterior corona radiata, and left posterior corona radiata. In patients carrying G allele, FA in body of corpus callosum was positively correlated with working memory, RD in genu of corpus callosum was negatively associated with the speed of processing, working memory, and the composite score of MCCB, while no significant correlations were found in AA homozygotes. In our study, patients with FES presented abnormal white matter integrity in corpus callosum and corona radiata. Furthermore, this abnormality was associated with the genetic variation of DAOA rs2391191, with AA homozygotes showing less white matter integrity in the corpus callosum. Our findings possibly provide further support to the evidence that DAOA regulates the process of glutamatergic neurotransmission and mitochondrial function in the pathophysiological mechanism of schizophrenia.

White matter abnormalities across different epilepsy syndromes in adults: an ENIGMA-Epilepsy study

The epilepsies are commonly accompanied by widespread abnormalities in cerebral white matter. ENIGMA-Epilepsy is a large quantitative brain imaging consortium, aggregating data to investigate patterns of neuroimaging abnormalities in common epilepsy syndromes, including temporal lobe epilepsy, extratemporal epilepsy, and genetic generalized epilepsy. Our goal was to rank the most robust white matter microstructural differences across and within syndromes in a multicentre sample of adult epilepsy patients. Diffusion-weighted MRI data were analysed from 1069 healthy controls and 1249 patients: temporal lobe epilepsy with hippocampal sclerosis (n = 599), temporal lobe epilepsy with normal MRI (n = 275), genetic generalized epilepsy (n = 182) and non-lesional extratemporal epilepsy (n = 193). A harmonized protocol using tract-based spatial statistics was used to derive skeletonized maps of fractional anisotropy and mean diffusivity for each participant, and fibre tracts were segmented using a diffusion MRI atlas. Data were harmonized to correct for scanner-specific variations in diffusion measures using a batch-effect correction tool (ComBat). Analyses of covariance, adjusting for age and sex, examined differences between each epilepsy syndrome and controls for each white matter tract (Bonferroni corrected at P < 0.001). Across 'all epilepsies' lower fractional anisotropy was observed in most fibre tracts with small to medium effect sizes, especially in the corpus callosum, cingulum and external capsule. There were also less robust increases in mean diffusivity. Syndrome-specific fractional anisotropy and mean diffusivity differences were most pronounced in patients with hippocampal sclerosis in the ipsilateral parahippocampal cingulum and external capsule, with smaller effects across most other tracts. Individuals with temporal lobe epilepsy and normal MRI showed a similar pattern of greater ipsilateral than contralateral abnormalities, but less marked than those in patients with hippocampal sclerosis. Patients with generalized and extratemporal epilepsies had pronounced reductions in fractional anisotropy in the corpus callosum, corona radiata and external capsule, and increased mean diffusivity of the anterior corona radiata. Earlier age of seizure onset and longer disease duration were associated with a greater extent of diffusion abnormalities in patients with hippocampal sclerosis. We demonstrate microstructural abnormalities across major association, commissural, and projection fibres in a large multicentre study of epilepsy. Overall, patients with epilepsy showed white matter abnormalities in the corpus callosum, cingulum and external capsule, with differing severity across epilepsy syndromes. These data further define the spectrum of white matter abnormalities in common epilepsy syndromes, yielding more detailed insights into pathological substrates that may explain cognitive and psychiatric co-morbidities and be used to guide biomarker studies of treatment outcomes and/or genetic research.

Neuroimaging studies within Cognitive Genetics Collaborative Research Organization aiming to replicate and extend works of ENIGMA

Reproducibility is one of the most important issues for generalizing the results of clinical research; however, low reproducibility in neuroimaging studies is well known. To overcome this problem, the Enhancing Neuroimaging Genetics through Meta‐Analysis (ENIGMA) consortium, an international neuroimaging consortium, established standard protocols for imaging analysis and employs either meta‐ and mega‐analyses of psychiatric disorders with large sample sizes. The Cognitive Genetics Collaborative Research Organization (COCORO) in Japan promotes neurobiological studies in psychiatry and has successfully replicated and extended works of ENIGMA especially for neuroimaging studies. For example, (a) the ENIGMA consortium showed subcortical regional volume alterations in patients with schizophrenia (n = 2,028) compared to controls (n = 2,540) across 15 cohorts using meta‐analysis. COCORO replicated the volumetric changes in patients with schizophrenia (n = 884) compared to controls (n = 1,680) using the ENIGMA imaging analysis protocol and mega‐analysis. Furthermore, a schizophrenia‐specific leftward asymmetry for the pallidum volume was demonstrated; and (b) the ENIGMA consortium identified white matter microstructural alterations in patients with schizophrenia (n = 1,963) compared to controls (n = 2,359) across 29 cohorts. Using the ENIGMA protocol, a study from COCORO showed similar results in patients with schizophrenia (n = 696) compared to controls (n = 1,506) from 12 sites using mega‐analysis. Moreover, the COCORO study found that schizophrenia, bipolar disorder (n = 211) and autism spectrum disorder (n = 126), but not major depressive disorder (n = 398), share similar white matter microstructural alterations, compared to controls. Further replication and harmonization of the ENIGMA consortium and COCORO will contribute to the generalization of their research findings.

Harmonization of Brain Diffusion MRI: Concepts and Methods

MRI diffusion data suffers from significant inter- and intra-site variability, which hinders multi-site and/or longitudinal diffusion studies. This variability may arise from a range of factors, such as hardware, reconstruction algorithms and acquisition settings. To allow a reliable comparison and joint analysis of diffusion data across sites and over time, there is a clear need for robust data harmonization methods. This review article provides a comprehensive overview of diffusion data harmonization concepts and methods, and their limitations. Overall, the methods for the harmonization of multi-site diffusion images can be categorized in two main groups: diffusion parametric map harmonization (DPMH) and diffusion weighted image harmonization (DWIH). Whereas DPMH harmonizes the diffusion parametric maps (e.g., FA, MD, and MK), DWIH harmonizes the diffusion-weighted images. Defining a gold standard harmonization technique for dMRI data is still an ongoing challenge. Nevertheless, in this paper we provide two classification tools, namely a feature table and a flowchart, which aim to guide the readers in selecting an appropriate harmonization method for their study.

Cingulum and abnormal psychological stress response in schizophrenia

Stress is implicated in many aspects of schizophrenia, including heightened distress intolerance. We examined how affect and microstructure of major brain tracts involved in regulating affect may contribute to distress intolerance in schizophrenia. Patients with schizophrenia spectrum disorders (n = 78) and community controls (n = 95) completed diffusion weighted imaging and performed psychological stress tasks. Subjective affect was collected pre and post stressors. Individuals who did not persist during one or both stress tasks were considered distress intolerant (DI), and otherwise distress tolerant (DT). Fractional anisotropy (FA) of the dorsal cingulum showed a significant diagnosis x DT/DI phenotype interaction (p = 0.003). Post-hoc tests showed dorsal cingulum FA was significantly lower in DI patients compared with DI controls (p < 0.001), but not different between DT groups (p = 0.27). Regarding affect responses to stress, irritability showed the largest stress-related change (p < 0.001), but irritability changes were significantly reduced in DI patients compared to DI controls (p = 0.006). The relationship between irritability change and performance errors also differed among patients (ρ = -0.29, p = 0.011) and controls (ρ = 0.21, p = 0.042). Further modeling highlighted the explanatory power of dorsal cingulum for predicting DI even after performance and irritability were taken into account. Distress intolerance during psychological stress exposure is related to microstructural properties of the dorsal cingulum, a key structure for cognitive control and emotion regulation. In schizophrenia, the affective response to psychological stressors is abnormal, and distress intolerant patients had significantly reduced dorsal cingulum FA compared to distress intolerant controls. The findings provide new insight regarding distress intolerance in schizophrenia.

Inter-site harmonization based on dual generative adversarial networks for diffusion tensor imaging: application to neonatal white matter development

Background

Site-specific variations are challenges for pooling analyses in multi-center studies. This work aims to propose an inter-site harmonization method based on dual generative adversarial networks (GANs) for diffusion tensor imaging (DTI) derived metrics on neonatal brains.

Results

DTI-derived metrics (fractional anisotropy, FA; mean diffusivity, MD) are obtained on age-matched neonates without magnetic resonance imaging (MRI) abnormalities: 42 neonates from site 1 and 42 neonates from site 2. Significant inter-site differences of FA can be observed. The proposed harmonization approach and three conventional methods (the global-wise scaling, the voxel-wise scaling, and the ComBat) are performed on DTI-derived metrics from two sites. During the tract-based spatial statistics, inter-site differences can be removed by the proposed dual GANs method, the voxel-wise scaling, and the ComBat. Among these methods, the proposed method holds the lowest median values in absolute errors and root mean square errors. During the pooling analysis of two sites, Pearson correlation coefficients between FA and the postmenstrual age after harmonization are larger than those before harmonization. The effect sizes (Cohen’s d between males and females) are also maintained by the harmonization procedure.

Conclusions

The proposed dual GANs-based harmonization method is effective to harmonize neonatal DTI-derived metrics from different sites. Results in this study further suggest that the GANs-based harmonization is a feasible pre-processing method for pooling analyses in multi-center studies.