ENIGMA and the individual: Predicting factors that affect the brain in 35 countries worldwide

The power of many brains: Catalyzing neuropsychiatric discovery through open neuroimaging data and large-scale collaboration

Recent advances in open neuroimaging data are enhancing our comprehension of neuropsychiatric disorders. By pooling images from various cohorts, statistical power has increased, enabling the detection of subtle abnormalities and robust associations, and fostering new research methods. Global collaborations in imaging have furthered our knowledge of the neurobiological foundations of brain disorders and aided in imaging-based prediction for more targeted treatment. Large-scale magnetic resonance imaging initiatives are driving innovation in analytics and supporting generalizable psychiatric studies. We also emphasize the significant role of big data in understanding neural mechanisms and in the early identification and precise treatment of neuropsychiatric disorders. However, challenges such as data harmonization across different sites, privacy protection, and effective data sharing must be addressed. With proper governance and open science practices, we conclude with a projection of how large-scale imaging resources and collaborations could revolutionize diagnosis, treatment selection, and outcome prediction, contributing to optimal brain health.

Functional magnetic resonance imaging in schizophrenia: current evidence, methodological advances, limitations and future directions

Functional neuroimaging emerged with great promise and has provided fundamental insights into the neurobiology of schizophrenia. However, it has faced challenges and criticisms, most notably a lack of clinical translation. This paper provides a comprehensive review and critical summary of the literature on functional neuroimaging, in particular functional magnetic resonance imaging (fMRI), in schizophrenia. We begin by reviewing research on fMRI biomarkers in schizophrenia and the clinical high risk phase through a historical lens, moving from case-control regional brain activation to global connectivity and advanced analytical approaches, and more recent machine learning algorithms to identify predictive neuroimaging features. Findings from fMRI studies of negative symptoms as well as of neurocognitive and social cognitive deficits are then reviewed. Functional neural markers of these symptoms and deficits may represent promising treatment targets in schizophrenia. Next, we summarize fMRI research related to antipsychotic medication, psychotherapy and psychosocial interventions, and neurostimulation, including treatment response and resistance, therapeutic mechanisms, and treatment targeting. We also review the utility of fMRI and data-driven approaches to dissect the heterogeneity of schizophrenia, moving beyond case-control comparisons, as well as methodological considerations and advances, including consortia and precision fMRI. Lastly, limitations and future directions of research in the field are discussed. Our comprehensive review suggests that, in order for fMRI to be clinically useful in the care of patients with schizophrenia, research should address potentially actionable clinical decisions that are routine in schizophrenia treatment, such as which antipsychotic should be prescribed or whether a given patient is likely to have persistent functional impairment. The potential clinical utility of fMRI is influenced by and must be weighed against cost and accessibility factors. Future evaluations of the utility of fMRI in prognostic and treatment response studies may consider including a health economics analysis.

DeepResBat: deep residual batch harmonization accounting for covariate distribution differences

Pooling MRI data from multiple datasets requires harmonization to reduce undesired inter-site variabilities, while preserving effects of biological variables (or covariates). The popular harmonization approach ComBat uses a mixed effect regression framework that explicitly accounts for covariate distribution differences across datasets. There is also significant interest in developing harmonization approaches based on deep neural networks (DNNs), such as conditional variational autoencoder (cVAE). However, current DNN approaches do not explicitly account for covariate distribution differences across datasets. Here, we provide mathematical results, suggesting that not accounting for covariates can lead to suboptimal harmonization outcomes. We propose two DNN-based harmonization approaches that explicitly account for covariate distribution differences across datasets: covariate VAE (coVAE) and DeepResBat. The coVAE approach is a natural extension of cVAE by concatenating covariates and site information with site- and covariate-invariant latent representations. DeepResBat adopts a residual framework inspired by ComBat. DeepResBat first removes the effects of covariates with nonlinear regression trees, followed by eliminating site differences with cVAE. Finally, covariate effects are added back to the harmonized residuals. Using three datasets from three different continents with a total of 2787 participants and 10085 anatomical T1 scans, we find that DeepResBat and coVAE outperformed ComBat, CovBat and cVAE in terms of removing dataset differences, while enhancing biological effects of interest. However, coVAE hallucinates spurious associations between anatomical MRI and covariates even when no association exists. Therefore, future studies proposing DNN-based harmonization approaches should be aware of this false positive pitfall. Overall, our results suggest that DeepResBat is an effective deep learning alternative to ComBat.

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.

A neurocognitive model of early onset persistent and desistant antisocial behavior in early adulthood

It remains unclear which functional and neurobiological mechanisms are associated with persistent and desistant antisocial behavior in early adulthood. We reviewed the empirical literature and propose a neurocognitive social information processing model for early onset persistent and desistant antisocial behavior in early adulthood, focusing on how young adults evaluate, act upon, monitor, and learn about their goals and self traits. Based on the reviewed literature, we propose that persistent antisocial behavior is characterized by domain-general impairments in self-relevant and goal-related information processing, regulation, and learning, which is accompanied by altered activity in fronto-limbic brain areas. We propose that desistant antisocial development is associated with more effortful information processing, regulation and learning, that possibly balances self-relevant goals and specific situational characteristics. The proposed framework advances insights by considering individual differences such as psychopathic personality traits, and specific emotional characteristics (e.g., valence of social cues), to further illuminate functional and neural mechanisms underlying heterogenous developmental pathways. Finally, we address important open questions and offer suggestions for future research to improve scientific knowledge on general and context-specific expression and development of antisocial behavior in early adulthood.

A systematic review of neuroimaging epigenetic research: calling for an increased focus on development

Epigenetic mechanisms, such as DNA methylation (DNAm), have gained increasing attention as potential biomarkers and mechanisms underlying risk for neurodevelopmental, psychiatric and other brain-based disorders. Yet, surprisingly little is known about the extent to which DNAm is linked to individual differences in the brain itself, and how these associations may unfold across development - a time of life when many of these disorders emerge. Here, we systematically review evidence from the nascent field of Neuroimaging Epigenetics, combining structural or functional neuroimaging measures with DNAm, and the extent to which the developmental period (birth to adolescence) is represented in these studies. We identified 111 articles published between 2011-2021, out of which only a minority (21%) included samples under 18 years of age. Most studies were cross-sectional (85%), employed a candidate-gene approach (67%), and examined DNAm-brain associations in the context of health and behavioral outcomes (75%). Nearly half incorporated genetic data, and a fourth investigated environmental influences. Overall, studies support a link between peripheral DNAm and brain imaging measures, but there is little consistency in specific findings and it remains unclear whether DNAm markers present a cause, correlate or consequence of brain alterations. Overall, there is large heterogeneity in sample characteristics, peripheral tissue and brain outcome examined as well as the methods used. Sample sizes were generally low to moderate (median nall = 98, ndevelopmental = 80), and attempts at replication or meta-analysis were rare. Based on the strengths and weaknesses of existing studies, we propose three recommendations on how advance the field of Neuroimaging Epigenetics. We advocate for: (1) a greater focus on developmentally oriented research (i.e. pre-birth to adolescence); (2) the analysis of large, prospective, pediatric cohorts with repeated measures of DNAm and imaging to assess directionality; and (3) collaborative, interdisciplinary science to identify robust signals, triangulate findings and enhance translational potential.

Enhancing collaborative neuroimaging research: introducing COINSTAC Vaults for federated analysis and reproducibility

Collaborative neuroimaging research is often hindered by technological, policy, administrative, and methodological barriers, despite the abundance of available data. COINSTAC (The Collaborative Informatics and Neuroimaging Suite Toolkit for Anonymous Computation) is a platform that successfully tackles these challenges through federated analysis, allowing researchers to analyze datasets without publicly sharing their data. This paper presents a significant enhancement to the COINSTAC platform: COINSTAC Vaults (CVs). CVs are designed to further reduce barriers by hosting standardized, persistent, and highly-available datasets, while seamlessly integrating with COINSTAC's federated analysis capabilities. CVs offer a user-friendly interface for self-service analysis, streamlining collaboration, and eliminating the need for manual coordination with data owners. Importantly, CVs can also be used in conjunction with open data as well, by simply creating a CV hosting the open data one would like to include in the analysis, thus filling an important gap in the data sharing ecosystem. We demonstrate the impact of CVs through several functional and structural neuroimaging studies utilizing federated analysis showcasing their potential to improve the reproducibility of research and increase sample sizes in neuroimaging studies.

Enhancing Collaborative Neuroimaging Research: Introducing COINSTAC Vaults for Federated Analysis and Reproducibility

Collaborative neuroimaging research is often hindered by technological, policy, administrative, and methodological barriers, despite the abundance of available data. COINSTAC is a platform that successfully tackles these challenges through federated analysis, allowing researchers to analyze datasets without publicly sharing their data. This paper presents a significant enhancement to the COINSTAC platform: COINSTAC Vaults (CVs). CVs are designed to further reduce barriers by hosting standardized, persistent, and highly-available datasets, while seamlessly integrating with COINSTAC's federated analysis capabilities. CVs offer a user-friendly interface for self-service analysis, streamlining collaboration and eliminating the need for manual coordination with data owners. Importantly, CVs can also be used in conjunction with open data as well, by simply creating a CV hosting the open data one would like to include in the analysis, thus filling an important gap in the data sharing ecosystem. We demonstrate the impact of CVs through several functional and structural neuroimaging studies utilizing federated analysis showcasing their potential to improve the reproducibility of research and increase sample sizes in neuroimaging studies.

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.

Diffusion-Weighted Imaging in Mild Traumatic Brain Injury: A Systematic Review of the Literature

There is evidence that diffusion-weighted imaging (DWI) is able to detect tissue alterations following mild traumatic brain injury (mTBI) that may not be observed on conventional neuroimaging; however, findings are often inconsistent between studies. This systematic review assesses patterns of differences in DWI metrics between those with and without a history of mTBI. A PubMed literature search was performed using relevant indexing terms for articles published prior to May 14, 2020. Findings were limited to human studies using DWI in mTBI. Articles were excluded if they were not full-length, did not contain original data, if they were case studies, pertained to military populations, had inadequate injury severity classification, or did not report post-injury interval. Findings were reported independently for four subgroups: acute/subacute pediatric mTBI, acute/subacute adult mTBI, chronic adult mTBI, and sport-related concussion, and all DWI acquisition and analysis methods used were included. Patterns of findings between studies were reported, along with strengths and weaknesses of the current state of the literature. Although heterogeneity of sample characteristics and study methods limited the consistency of findings, alterations in DWI metrics were most commonly reported in the corpus callosum, corona radiata, internal capsule, and long association pathways. Many acute/subacute pediatric studies reported higher FA and lower ADC or MD in various regions. In contrast, acute/subacute adult studies most commonly indicate lower FA within the context of higher MD and RD. In the chronic phase of recovery, FA may remain low, possibly indicating overall demyelination or Wallerian degeneration over time. Longitudinal studies, though limited, generally indicate at least a partial normalization of DWI metrics over time, which is often associated with functional improvement. We conclude that DWI is able to detect structural mTBI-related abnormalities that may persist over time, although future DWI research will benefit from larger samples, improved data analysis methods, standardized reporting, and increasing transparency.

Connectome-based predictive modeling of compulsion in obsessive-compulsive disorder

Compulsion is one of core symptoms of obsessive-compulsive disorder (OCD). Although many studies have investigated the neural mechanism of compulsion, no study has used brain-based measures to predict compulsion. Here, we used connectome-based predictive modeling (CPM) to identify networks that could predict the levels of compulsion based on whole-brain functional connectivity in 57 OCD patients. We then applied a computational lesion version of CPM to examine the importance of specific brain areas. We also compared the predictive network strength in OCD with unaffected first-degree relatives (UFDR) of patients and healthy controls. CPM successfully predicted individual level of compulsion and identified networks positively (primarily subcortical areas of the striatum and limbic regions of the hippocampus) and negatively (primarily frontoparietal regions) correlated with compulsion. The prediction power of the negative model significantly decreased when simulating lesions to the prefrontal cortex and cerebellum, supporting the importance of these regions for compulsion prediction. We found a similar pattern of network strength in the negative predictive network for OCD patients and their UFDR, demonstrating the potential of CPM to identify vulnerability markers for psychopathology.

Machine Learning and MRI-based Diagnostic Models for ADHD: Are We There Yet?

OBJECTIVE

Machine learning (ML) has been applied to develop magnetic resonance imaging (MRI)-based diagnostic classifiers for attention-deficit/hyperactivity disorder (ADHD). This systematic review examines this literature to clarify its clinical significance and to assess the implications of the various analytic methods applied.

METHODS

A comprehensive literature search on MRI-based diagnostic classifiers for ADHD was performed and data regarding the utilized models and samples were gathered.

RESULTS

We found that, although most studies reported the classification accuracies, they varied in choice of MRI modalities, ML models, cross-validation and testing methods, and sample sizes. We found that the accuracies of cross-validation methods inflated the performance estimation compared with those of a held-out test, compromising the model generalizability. Test accuracies have increased with publication year but were not associated with training sample sizes. Improved test accuracy over time was likely due to the use of better ML methods along with strategies to deal with data imbalances.

CONCLUSION

Ultimately, large multi-modal imaging datasets, and potentially the combination with other types of data, like cognitive data and/or genetics, will be essential to achieve the goal of developing clinically useful imaging classification tools for ADHD in the future.

The neurobiology of obsessive-compulsive personality disorder: a systematic review

BACKGROUND

A number of recent investigations have focused on the neurobiology of obsessive-compulsive personality disorder (OCPD). However, there have been few reviews of this literature with no detailed model proposed. We therefore undertook a systematic review of these investigations, aiming to map the available evidence and investigate whether it is possible to formulate a detailed model of the neurobiology of OCPD.

METHODS

OCPD can be considered from both categorical and dimensional perspectives. An electronic search was therefore conducted using terms that would address not only OCPD as a category, but also related constructs, such as perfectionism, that would capture research on neuropsychology, neuroimaging, neurochemistry, and neurogenetics.

RESULTS

A total of 1059 articles were retrieved, with 87 ultimately selected for abstract screening, resulting in a final selection of 49 articles focusing on neurobiological investigations relevant to OCPD. Impaired executive function and cognitive inflexibility are common neuropsychological traits in this condition, and suggest that obsessive-compulsive disorder (OCD) and OCPD may lie on a continuum. However, neuroimaging studies in OCPD indicate the involvement of specific neurocircuitry, including the precuneus and amygdala, and so suggest that OCD and OCPD may have important differences. Although OCPD has a heritable component, we found no well-powered genetic studies of OCPD.

CONCLUSION

Although knowledge in this area has advanced, there are insufficient data on which to base a comprehensive model of the neurobiology of OCPD. Given the clinical importance of OCPD, further work to understand the mechanisms that underpin this condition is warranted.

Association between CILP and IL-1α polymorphisms and phenotype-dependent intervertebral disc degeneration susceptibility: A meta-analysis

Background: The relationship between CILP (1184T>C) and IL-1α(+889C/T) polymorphisms and intervertebral disc degeneration (IDD) have been explored in several studies but the results were conflicting. The aim of the study was to evaluate and synthesize the currently available data on the association between CILP (1184T>C) and IL-1α(+889C/T) polymorphisms and susceptibility of phenotype-dependent radiologic IDD (RIDD) and symptomatic intervertebral disk herniation (SIDH).

Methods: A computerized literature search was in PubMed, Cochrane Library, Embase, China National Knowledge Infrastructure database, and Web of Science. The pooled results were presented as odds ratios (ORs) with 95% confidence intervals (CIs). Moreover, the false-positive report probability (FPRP) test and trial sequential analysis (TSA) were applied to estimate the significant results.

Results: Our evidence demonstrated that IL-1α(+889C/T) was significant associated with RIDD (allele model: OR = 1.34, 95%CI 1.03–1.74, p = 0.029) and SIDH (allele model: OR = 1.28, 95% CI 1.03–1.60, p = 0.028). However, the results were not noteworthy under the FPRP test and TSA analysis. Additionally, CILP (1184T>C) polymorphism was significantly associated with RIDD with adequate evidence (allele model: OR = 1.27, 95% CI 1.09–1.48, p = 0.002) instead of SIDH.

Conclusion: The current meta-analysis illustrated firm evidence that CILP (1184T>C) polymorphism was significantly associated with the susceptibility of RIDD. However, the significant associations between IL-1α(+889C/T) and RIDD and SIDH were less credible. Thus, more multi-center studies with diverse populations were required to verify the results.

The relationship between fat mass and obesity associated gene polymorphism rs9939609 and resting cerebral blood flow in a midlife sample with overweight and obesity

Background

The single nucleotide polymorphism (SNP) rs9939609 in the fat mass and obesity associated fat mass and obesity associated gene (FTO) gene has been linked with increased BMI in adults. Higher BMI has been associated with poor brain health and may exert deleterious effects on neurocognitive health through cerebral hypoperfusion. However, it is unclear if there is a relationship between the FTO genotype and cerebral perfusion, or whether FTO genotype moderates the effects of weight loss on cerebral perfusion. Using data from a randomized controlled behavioral weight loss trial in adults with overweight and obesity, we tested (1) whether carriers of the A allele for FTO rs9939609 demonstrate different patterns of resting cerebral blood flow (rCBF) compared to T carriers, and (2) whether the FTO genotype moderates the effects of weight loss on rCBF. We hypothesized that carriers of the A allele would exhibit lower resting CBF in frontal brain areas compared to T/T homozygotes at baseline, and that intervention-induced weight loss may partially remediate these differences.

Methods and results

One hundred and five adults (75.2% female, mean age 44.9 years) with overweight or obesity were included in the analyses. These participants represent a subsample of participants in a larger randomized controlled trial (NCT01500356). A resting pseudo-continuous arterial spin labeling (pCASL) scan was acquired to examine rCBF. Age, sex, and BMI were included as covariates. At baseline, A carriers had greater rCBF in a diffuse cluster extending into the brainstem, motor cortex, and occipital lobe, but lower perfusion in the temporal lobe. We found no evidence that FTO moderated the effect of the intervention group assignment on rCBF changes.

Conclusion

Overall, these results indicate that (a) individual variation in rCBF within a sample with overweight and obesity may be attributed to a common FTO variant, but (b) a weight loss intervention is effective at increasing rCBF, regardless of FTO genotype.

Goal-specific brain MRI harmonization

There is significant interest in pooling magnetic resonance image (MRI) data from multiple datasets to enable mega-analysis. Harmonization is typically performed to reduce heterogeneity when pooling MRI data across datasets. Most MRI harmonization algorithms do not explicitly consider downstream application performance during harmonization. However, the choice of downstream application might influence what might be considered as study-specific confounds. Therefore, ignoring downstream applications during harmonization might potentially limit downstream performance. Here we propose a goal-specific harmonization framework that utilizes downstream application performance to regularize the harmonization procedure. Our framework can be integrated with a wide variety of harmonization models based on deep neural networks, such as the recently proposed conditional variational autoencoder (cVAE) harmonization model. Three datasets from three different continents with a total of 2787 participants and 10,085 anatomical T1 scans were used for evaluation. We found that cVAE removed more dataset differences than the widely used ComBat model, but at the expense of removing desirable biological information as measured by downstream prediction of mini mental state examination (MMSE) scores and clinical diagnoses. On the other hand, our goal-specific cVAE (gcVAE) was able to remove as much dataset differences as cVAE, while improving downstream cross-sectional prediction of MMSE scores and clinical diagnoses.

Cross disorder comparisons of brain structure in schizophrenia, bipolar disorder, major depressive disorder, and 22q11.2 deletion syndrome: A review of ENIGMA findings

This review compares the main brain abnormalities in schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and 22q11.2 Deletion Syndrome (22q11DS) determined by ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) consortium investigations. We obtained ranked effect sizes for subcortical volumes, regional cortical thickness, cortical surface area, and diffusion tensor imaging abnormalities, comparing each of these disorders relative to healthy controls. In addition, the studies report on significant associations between brain imaging metrics and disorder-related factors such as symptom severity and treatments. Visual comparison of effect size profiles shows that effect sizes are generally in the same direction and scale in severity with the disorders (in the order SZ > BD > MDD). The effect sizes for 22q11DS, a rare genetic syndrome that increases the risk for psychiatric disorders, appear to be much larger than for either of the complex psychiatric disorders. This is consistent with the idea of generally larger effects on the brain of rare compared to common genetic variants. Cortical thickness and surface area effect sizes for 22q11DS with psychosis compared to 22q11DS without psychosis are more similar to those of SZ and BD than those of MDD; a pattern not observed for subcortical brain structures and fractional anisotropy effect sizes. The observed similarities in effect size profiles for cortical measures across the psychiatric disorders mimic those observed for shared genetic variance between these disorders reported based on family and genetic studies and are consistent with shared genetic risk for SZ and BD and structural brain phenotypes.

Quantitative mapping of the brain's structural connectivity using diffusion MRI tractography: A review

Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.

A Decentralized ComBat Algorithm and Applications to Functional Network Connectivity

Recent studies showed that working with neuroimage data collected from different research facilities or locations may incur additional source dependency, affecting the overall statistical power. This problem can be mitigated with data harmonization approaches. Recently, the ComBat method has become commonly adopted for various neuroimage modalities. While open neuroimaging datasets are becoming more common, a substantial amount of data is still unable to be shared for various reasons. In addition, current approaches require moving all the data to a central location, which requires additional resources and creates redundant copies of the same datasets. To address these issues, we propose a decentralized harmonization approach that does not create redundant copies of the original datasets and performs remote operations on the datasets separately without sharing any individual subject data, ensuring a certain level of privacy and reducing regulatory hurdles. We proposed a novel approach called "Decentralized ComBat" which can harmonize datasets separately without combining the datasets. We tested our model by harmonizing functional network connectivity datasets from two traumatic brain injury studies in a decentralized way. Also, we used simulations to analyze the performance and scalability of our model when the number of data collection sites increases. We compare the output with centralized ComBat and show that the proposed approach produces similar results, increasing the sensitivity of the functional network connectivity analysis and validating our approach. Simulations show that our model can be easily scaled to many more datasets based on the requirement. In sum, we believe this provides a powerful tool, further complementing open data and allowing for integrating public and private datasets.

Shared Genetic Etiology between Cortical Brain Morphology and Tobacco, Alcohol, and Cannabis Use

Genome-wide association studies (GWAS) have identified genetic variants associated with brain morphology and substance use behaviors (SUB). However, the genetic overlap between brain structure and SUB has not been well characterized. We leveraged GWAS summary data of 71 brain imaging measures and alcohol, tobacco, and cannabis use to investigate their genetic overlap using linkage disequilibrium score regression. We used genomic structural equation modeling to model a "common SUB genetic factor" and investigated its genetic overlap with brain structure. Furthermore, we estimated SUB polygenic risk scores (PRS) and examined whether they predicted brain imaging traits using the Adolescent Behavior and Cognitive Development (ABCD) study. We identified 8 significant negative genetic correlations, including between (1) alcoholic drinks per week and average cortical thickness, and (2) intracranial volume with age of smoking initiation. We observed 5 positive genetic correlations, including those between (1) insula surface area and lifetime cannabis use, and (2) the common SUB genetic factor and pericalcarine surface area. SUB PRS were associated with brain structure variation in ABCD. Our findings highlight a shared genetic etiology between cortical brain morphology and SUB and suggest that genetic variants associated with SUB may be causally related to brain structure differences.

Subcortical volumes across the lifespan: Data from 18,605 healthy individuals aged 3-90 years

Age has a major effect on brain volume. However, the normative studies available are constrained by small sample sizes, restricted age coverage and significant methodological variability. These limitations introduce inconsistencies and may obscure or distort the lifespan trajectories of brain morphometry. In response, we capitalized on the resources of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to examine age-related trajectories inferred from cross-sectional measures of the ventricles, the basal ganglia (caudate, putamen, pallidum, and nucleus accumbens), the thalamus, hippocampus and amygdala using magnetic resonance imaging data obtained from 18,605 individuals aged 3-90 years. All subcortical structure volumes were at their maximum value early in life. The volume of the basal ganglia showed a monotonic negative association with age thereafter; there was no significant association between age and the volumes of the thalamus, amygdala and the hippocampus (with some degree of decline in thalamus) until the sixth decade of life after which they also showed a steep negative association with age. The lateral ventricles showed continuous enlargement throughout the lifespan. Age was positively associated with inter-individual variability in the hippocampus and amygdala and the lateral ventricles. These results were robust to potential confounders and could be used to examine the functional significance of deviations from typical age-related morphometric patterns.

A time-course study of actively stained mouse brains: Diffusion tensor imaging parameters and connectomic stability over 1 year

While the application of diffusion tensor imaging (DTI), tractography, and connectomics to fixed tissue is a common practice today, there have been limited studies examining the effects of fixation on brain microstructure over extended periods. This mouse model time-course study reports the changes of regional brain volumes and diffusion scalar parameters, such as fractional anisotropy, across 12 representative brain regions as measures of brain structural stability. The scalar DTI parameters and regional volumes were highly variable over the first 2 weeks after fixation. The same parameters were consistent over a 2-8-week window after fixation, which means confounds from tissue stability over that scanning window were minimal. Quantitative connectomes were analyzed over the same time with extension out to 1 year. While there was some change in the scalar metrics at 1 year after fixation, these changes were sufficiently small, particularly in white matter, to support reproducible connectomes over a period ranging from 2-weeks to 1-year post-fixation. These findings delineate a scanning period, during which brain volumes, diffusion scalar metrics, and connectomes are remarkably consistent.

Cortical thickness across the lifespan: Data from 17,075 healthy individuals aged 3-90 years

Delineating the association of age and cortical thickness in healthy individuals is critical given the association of cortical thickness with cognition and behavior. Previous research has shown that robust estimates of the association between age and brain morphometry require large-scale studies. In response, we used cross-sectional data from 17,075 individuals aged 3-90 years from the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to infer age-related changes in cortical thickness. We used fractional polynomial (FP) regression to quantify the association between age and cortical thickness, and we computed normalized growth centiles using the parametric Lambda, Mu, and Sigma method. Interindividual variability was estimated using meta-analysis and one-way analysis of variance. For most regions, their highest cortical thickness value was observed in childhood. Age and cortical thickness showed a negative association; the slope was steeper up to the third decade of life and more gradual thereafter; notable exceptions to this general pattern were entorhinal, temporopolar, and anterior cingulate cortices. Interindividual variability was largest in temporal and frontal regions across the lifespan. Age and its FP combinations explained up to 59% variance in cortical thickness. These results may form the basis of further investigation on normative deviation in cortical thickness and its significance for behavioral and cognitive outcomes.

SHIP-MR and Radiology: 12 Years of Whole-Body Magnetic Resonance Imaging in a Single Center

The Study of Health in Pomerania (SHIP), a population-based study from a rural state in northeastern Germany with a relatively poor life expectancy, supplemented its comprehensive examination program in 2008 with whole-body MR imaging at 1.5 T (SHIP-MR). We reviewed more than 100 publications that used the SHIP-MR data and analyzed which sequences already produced fruitful scientific outputs and which manuscripts have been referenced frequently. Upon reviewing the publications about imaging sequences, those that used T1-weighted structured imaging of the brain and a gradient-echo sequence for R2* mapping obtained the highest scientific output; regarding specific body parts examined, most scientific publications focused on MR sequences involving the brain and the (upper) abdomen. We conclude that population-based MR imaging in cohort studies should define more precise goals when allocating imaging time. In addition, quality control measures might include recording the number and impact of published work, preferably on a bi-annual basis and starting 2 years after initiation of the study. Structured teaching courses may enhance the desired output in areas that appear underrepresented.

Neuroimaging Association Scores: reliability and validity of aggregate measures of brain structural features linked to mental disorders in youth

In genetics, aggregation of many loci with small effect sizes into a single score improved prediction. Nevertheless, studies applying easily replicable weighted scores to neuroimaging data are lacking. Our aim was to assess the reliability and validity of the Neuroimaging Association Score (NAS), which combines information from structural brain features previously linked to mental disorders. Participants were 726 youth (aged 6-14) from two cities in Brazil who underwent MRI and psychopathology assessment at baseline and 387 at 3-year follow-up. Results were replicated in two samples: IMAGEN (n = 1627) and the Healthy Brain Network (n = 843). NAS were derived by summing the product of each standardized brain feature by the effect size of the association of that brain feature with seven psychiatric disorders documented by previous meta-analyses. NAS were calculated for surface area, cortical thickness and subcortical volumes using T1-weighted scans. NAS reliability, temporal stability and psychopathology and cognition prediction were analyzed. NAS for surface area showed high internal consistency and 3-year stability and predicted general psychopathology and cognition with higher replicability than specific symptomatic domains for all samples. They also predicted general psychopathology with higher replicability than single structures alone, accounting for 1-3% of the variance, but without directionality. The NAS for cortical thickness and subcortical volumes showed lower internal consistency and less replicable associations with behavioural phenotypes. These findings indicate the NAS based on surface area might be replicable markers of general psychopathology, but these links are unlikely to be causal or clinically useful yet.

Federated Analysis of Neuroimaging Data: A Review of the Field

The field of neuroimaging has embraced sharing data to collaboratively advance our understanding of the brain. However, data sharing, especially across sites with large amounts of protected health information (PHI), can be cumbersome and time intensive. Recently, there has been a greater push towards collaborative frameworks that enable large-scale federated analysis of neuroimaging data without the data having to leave its original location. However, there still remains a need for a standardized federated approach that not only allows for data sharing adhering to the FAIR (Findability, Accessibility, Interoperability, Reusability) data principles, but also streamlines analyses and communication while maintaining subject privacy. In this paper, we review a non-exhaustive list of neuroimaging analytic tools and frameworks currently in use. We then provide an update on our federated neuroimaging analysis software system, the Collaborative Informatics and Neuroimaging Suite Toolkit for Anonymous Computation (COINSTAC). In the end, we share insights on future research directions for federated analysis of neuroimaging data.

A Correlated Noise-assisted Decentralized Differentially Private Estimation Protocol, and its application to fMRI Source Separation

Blind source separation algorithms such as independent component analysis (ICA) are widely used in the analysis of neuroimaging data. To leverage larger sample sizes, different data holders/sites may wish to collaboratively learn feature representations. However, such datasets are often privacy-sensitive, precluding centralized analyses that pool the data at one site. In this work, we propose a differentially private algorithm for performing ICA in a decentralized data setting. Due to the high dimension and small sample size, conventional approaches to decentralized differentially private algorithms suffer in terms of utility. When centralizing the data is not possible, we investigate the benefit of enabling limited collaboration in the form of generating jointly distributed random noise. We show that such (anti) correlated noise improves the privacy-utility trade-off, and can reach the same level of utility as the corresponding non-private algorithm for certain parameter choices. We validate this benefit using synthetic and real neuroimaging datasets. We conclude that it is possible to achieve meaningful utility while preserving privacy, even in complex signal processing systems.

Genomics of Obsessive-Compulsive Disorder-Toward Personalized Medicine in the Era of Big Data

Pathogenesis of obsessive-compulsive disorder (OCD) mainly involves dysregulation of serotonergic neurotransmission, but a number of other factors are involved. Genetic underprints of OCD fall under the category of "common disease common variant hypothesis," that suggests that if a disease that is heritable is common in the population (a prevalence >1-5%), then the genetic contributors-specific variations in the genetic code-will also be common in the population. Therefore, the genetic contribution in OCD is believed to come from multiple genes simultaneously and it is considered a polygenic disorder. Genomics offers a number of advanced tools to determine causal relationship between the exposure and the outcome of interest. Particularly, methods such as polygenic risk score (PRS) or Mendelian Randomization (MR) enable investigation of new pathways involved in OCD pathogenesis. This premise is also facilitated by the existence of publicly available databases that include vast study samples. Examples include population-based studies such as UK Biobank, China Kadoorie Biobank, Qatar Biobank, All of US Program sponsored by National Institute of Health or Generations launched by Yale University, as well as disease-specific databases, that include patients with OCD and co-existing pathologies, with the following examples: Psychiatric Genomics Consortium (PGC), ENIGMA OCD, The International OCD Foundation Genetics Collaborative (IOCDF-GC) or OCD Collaborative Genetic Association Study. The aim of this review is to present a comprehensive overview of the available Big Data resources for the study of OCD pathogenesis in the context of genomics and demonstrate that OCD should be considered a disorder which requires the approaches offered by personalized medicine.

Linking the genomic signatures of human beat synchronization and learned song in birds

The development of rhythmicity is foundational to communicative and social behaviours in humans and many other species, and mechanisms of synchrony could be conserved across species. The goal of the current paper is to explore evolutionary hypotheses linking vocal learning and beat synchronization through genomic approaches, testing the prediction that genetic underpinnings of birdsong also contribute to the aetiology of human interactions with musical beat structure. We combined state-of-the-art-genomic datasets that account for underlying polygenicity of these traits: birdsong genome-wide transcriptomics linked to singing in zebra finches, and a human genome-wide association study of beat synchronization. Results of competitive gene set analysis revealed that the genetic architecture of human beat synchronization is significantly enriched for birdsong genes expressed in songbird Area X (a key nucleus for vocal learning, and homologous to human basal ganglia). These findings complement ethological and neural evidence of the relationship between vocal learning and beat synchronization, supporting a framework of some degree of common genomic substrates underlying rhythm-related behaviours in two clades, humans and songbirds (the largest evolutionary radiation of vocal learners). Future cross-species approaches investigating the genetic underpinnings of beat synchronization in a broad evolutionary context are discussed. This article is part of the theme issue 'Synchrony and rhythm interaction: from the brain to behavioural ecology'.

A single mode of population covariation associates brain networks structure and behavior and predicts individual subjects' age

Multiple human behaviors improve early in life, peaking in young adulthood, and declining thereafter. Several properties of brain structure and function progress similarly across the lifespan. Cognitive and neuroscience research has approached aging primarily using associations between a few behaviors, brain functions, and structures. Because of this, the multivariate, global factors relating brain and behavior across the lifespan are not well understood. We investigated the global patterns of associations between 334 behavioral and clinical measures and 376 brain structural connections in 594 individuals across the lifespan. A single-axis associated changes in multiple behavioral domains and brain structural connections (r = 0.5808). Individual variability within the single association axis well predicted the age of the subject (r = 0.6275). Representational similarity analysis evidenced global patterns of interactions across multiple brain network systems and behavioral domains. Results show that global processes of human aging can be well captured by a multivariate data fusion approach.

Large-scale collaboration in ENIGMA-EEG: A perspective on the meta-analytic approach to link neurological and psychiatric liability genes to electrophysiological brain activity

BACKGROUND AND PURPOSE

The ENIGMA-EEG working group was established to enable large-scale international collaborations among cohorts that investigate the genetics of brain function measured with electroencephalography (EEG). In this perspective, we will discuss why analyzing the genetics of functional brain activity may be crucial for understanding how neurological and psychiatric liability genes affect the brain.

METHODS

We summarize how we have performed our currently largest genome-wide association study of oscillatory brain activity in EEG recordings by meta-analyzing the results across five participating cohorts, resulting in the first genome-wide significant hits for oscillatory brain function located in/near genes that were previously associated with psychiatric disorders. We describe how we have tackled methodological issues surrounding genetic meta-analysis of EEG features. We discuss the importance of harmonizing EEG signal processing, cleaning, and feature extraction. Finally, we explain our selection of EEG features currently being investigated, including the temporal dynamics of oscillations and the connectivity network based on synchronization of oscillations.

RESULTS

We present data that show how to perform systematic quality control and evaluate how choices in reference electrode and montage affect individual differences in EEG parameters.

CONCLUSION

The long list of potential challenges to our large-scale meta-analytic approach requires extensive effort and organization between participating cohorts; however, our perspective shows that these challenges are surmountable. Our perspective argues that elucidating the genetic of EEG oscillatory activity is a worthwhile effort in order to elucidate the pathway from gene to disease liability.

Exposure to methylphenidate during peri-adolescence decouples the prefrontal cortex: a multimodal MRI study

This study was designed to assess the effects of daily psychostimulant exposure during juvenility and peri-adolescence on brain morphology and functional connectivity using multimodal magnetic resonance imaging. We hypothesized that long-term exposure to methylphenidate would enhance connectivity with the prefrontal cortex. Male rats were given daily injections of either methylphenidate (n=10), dextroamphetamine (n=10) or saline vehicle (n=10) from postnatal day 21 to 42. They were imaged between postnatal day 43 and 48. Voxel-based morphometry, diffusion weighted imaging, and resting state functional connectivity were used to quantify brain structure and function. Images from each modality were registered and analyzed, using a 3D MRI rat atlas providing site-specific data over 171 different brain areas. Following imaging, rats were tested for cognitive function using novel object preference. Long-lasting psychostimulant treatment was associated with only a few significant changes in brain volume and measures of anisotropy compared to vehicle. Resting state functional connectivity imaging revealed decreased coupling between the prefrontal cortex, basal ganglia and sensory motor cortices. There were no significant differences between experimental groups for cognitive behavior. In this exploratory study, we showed that chronic psychostimulant treatment throughout juvenility and preadolescence has a minimal effect on brain volume and gray matter microarchitecture, but significantly uncouples the connectivity in the cerebral/basal ganglia circuitry.

Maturational trajectories of pericortical contrast in typical brain development

In the last few years, a significant amount of work has aimed to characterize maturational trajectories of cortical development. The role of pericortical microstructure putatively characterized as the gray-white matter contrast (GWC) at the pericortical gray-white matter boundary and its relationship to more traditional morphological measures of cortical morphometry has emerged as a means to examine finer grained neuroanatomical underpinnings of cortical changes. In this work, we characterize the GWC developmental trajectories in a representative sample (n = 394) of children and adolescents (~4 to ~22 years of age), with repeated scans (1-3 scans per subject, total scans n = 819). We tested whether linear, quadratic, or cubic trajectories of contrast development best described changes in GWC. A best-fit model was identified vertex-wise across the whole cortex via the Akaike Information Criterion (AIC). GWC across nearly the whole brain was found to significantly change with age. Cubic trajectories were likeliest for 63% of vertices, quadratic trajectories were likeliest for 20% of vertices, and linear trajectories were likeliest for 16% of vertices. A main effect of sex was observed in some regions, where males had a higher GWC than females. However, no sex by age interactions were found on GWC. In summary, our results suggest a progressive decrease in GWC at the pericortical boundary throughout childhood and adolescence. This work contributes to efforts seeking to characterize typical, healthy brain development and, by extension, can help elucidate aberrant developmental trajectories.

Prioritizing Genetic Contributors to Cortical Alterations in 22q11.2 Deletion Syndrome Using Imaging Transcriptomics

22q11.2 deletion syndrome (22q11DS) results from a hemizygous deletion that typically spans 46 protein-coding genes and is associated with widespread alterations in brain morphology. The specific genetic mechanisms underlying these alterations remain unclear. In the 22q11.2 ENIGMA Working Group, we characterized cortical alterations in individuals with 22q11DS (n = 232) versus healthy individuals (n = 290) and conducted spatial convergence analyses using gene expression data from the Allen Human Brain Atlas to prioritize individual genes that may contribute to altered surface area (SA) and cortical thickness (CT) in 22q11DS. Total SA was reduced in 22q11DS (Z-score deviance = -1.04), with prominent reductions in midline posterior and lateral association regions. Mean CT was thicker in 22q11DS (Z-score deviance = +0.64), with focal thinning in a subset of regions. Regional expression of DGCR8 was robustly associated with regional severity of SA deviance in 22q11DS; AIFM3 was also associated with SA deviance. Conversely, P2RX6 was associated with CT deviance. Exploratory analysis of gene targets of microRNAs previously identified as down-regulated due to DGCR8 deficiency suggested that DGCR8 haploinsufficiency may contribute to altered corticogenesis in 22q11DS by disrupting cell cycle modulation. These findings demonstrate the utility of combining neuroanatomic and transcriptomic datasets to derive molecular insights into complex, multigene copy number variants.

Neuroimaging PheWAS (Phenome-Wide Association Study): A Free Cloud-Computing Platform for Big-Data, Brain-Wide Imaging Association Studies

Large-scale, case-control genome-wide association studies (GWASs) have revealed genetic variations associated with diverse neurological and psychiatric disorders. Recent advances in neuroimaging and genomic databases of large healthy and diseased cohorts have empowered studies to characterize effects of the discovered genetic factors on brain structure and function, implicating neural pathways and genetic mechanisms in the underlying biology. However, the unprecedented scale and complexity of the imaging and genomic data requires new advanced biomedical data science tools to manage, process and analyze the data. In this work, we introduce Neuroimaging PheWAS (phenome-wide association study): a web-based system for searching over a wide variety of brain-wide imaging phenotypes to discover true system-level gene-brain relationships using a unified genotype-to-phenotype strategy. This design features a user-friendly graphical user interface (GUI) for anonymous data uploading, study definition and management, and interactive result visualizations as well as a cloud-based computational infrastructure and multiple state-of-art methods for statistical association analysis and multiple comparison correction. We demonstrated the potential of Neuroimaging PheWAS with a case study analyzing the influences of the apolipoprotein E (APOE) gene on various brain morphological properties across the brain in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. Benchmark tests were performed to evaluate the system's performance using data from UK Biobank. The Neuroimaging PheWAS system is freely available. It simplifies the execution of PheWAS on neuroimaging data and provides an opportunity for imaging genetics studies to elucidate routes at play for specific genetic variants on diseases in the context of detailed imaging phenotypic data.

Structural brain imaging studies offer clues about the effects of the shared genetic etiology among neuropsychiatric disorders

Genomewide association studies have found significant genetic correlations among many neuropsychiatric disorders. In contrast, we know much less about the degree to which structural brain alterations are similar among disorders and, if so, the degree to which such similarities have a genetic etiology. From the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium, we acquired standardized mean differences (SMDs) in regional brain volume and cortical thickness between cases and controls. We had data on 41 brain regions for: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), epilepsy, major depressive disorder (MDD), obsessive compulsive disorder (OCD), and schizophrenia (SCZ). These data had been derived from 24,360 patients and 37,425 controls. The SMDs were significantly correlated between SCZ and BD, OCD, MDD, and ASD. MDD was positively correlated with BD and OCD. BD was positively correlated with OCD and negatively correlated with ADHD. These pairwise correlations among disorders were correlated with the corresponding pairwise correlations among disorders derived from genomewide association studies (r = 0.494). Our results show substantial similarities in sMRI phenotypes among neuropsychiatric disorders and suggest that these similarities are accounted for, in part, by corresponding similarities in common genetic variant architectures.

Increased facial asymmetry in focal epilepsies associated with unilateral lesions

The epilepsies are now conceptualized as network disruptions: focal epilepsies are considered to have network alterations in the hemisphere of seizure onset, whilst generalized epilepsies are considered to have bi-hemispheric network changes. Increasingly, many epilepsies are also considered to be neurodevelopmental disorders, with early changes in the brain underpinning seizure biology. The development of the structure of the face is influenced by complex molecular interactions between surface ectoderm and underlying developing forebrain and neural crest cells. This influence is likely to continue postnatally, given the evidence of facial growth changes over time in humans until at least 18 years of age. In this case-control study, we hypothesized that people with lateralized focal epilepsies (i.e. unilateral network changes) have an increased degree of facial asymmetry, compared with people with generalized epilepsies or controls without epilepsy. We applied three-dimensional stereophotogrammetry and dense surface models to evaluate facial asymmetry in people with epilepsy, aiming to generate new tools to explore pathophysiological mechanisms in epilepsy. We analysed neuroimaging data to explore the correlation between face and brain asymmetry. We consecutively recruited 859 people with epilepsy attending the epilepsy clinics at a tertiary referral centre. We used dense surface modelling of the full face and signature analyses of three-dimensional facial photographs to analyse facial differences between 378 cases and 205 healthy controls. Neuroimaging around the time of the facial photograph was available for 234 cases. We computed the brain asymmetry index between contralateral regions. Cases with focal symptomatic epilepsy associated with unilateral lesions showed greater facial asymmetry compared to controls (P = 0.0001, two-sample t-test). This finding was confirmed by linear regression analysis after controlling for age and gender. We also found a significant correlation between duration of illness and the brain asymmetry index of total average cortical thickness (r = -0.19, P = 0.0075) but not for total average surface area (r = 0.06, P = 0.3968). There was no significant correlation between facial asymmetry and asymmetry of regional cortical thickness or surface area. We propose that the greater facial asymmetry in cases with focal epilepsy caused by unilateral abnormality might be explained by early unilateral network disruption, and that this is independent of underlying brain asymmetry. Three-dimensional stereophotogrammetry and dense surface modelling are a novel powerful phenotyping tool in epilepsy that may permit greater understanding of pathophysiology in epilepsy, and generate further insights into the development of cerebral networks underlying epilepsy, and the genetics of facial and neural development.

Case Report: Chemotherapy Indication in a Case of Neurofibromatosis Type 1 Presenting Optic Pathway Glioma: A One-Year Clinical Case Study Using Differential Tractography Approach

Neurofibromatosis type 1 (NF1) is associated with peripheral and central nervous system tumors. It is noteworthy that the regions in which these tumors frequently arise are the optic pathways (OPs) and the brainstem. Thus, we decided to trace the procedure of diffusion Magnetic Resonance Imaging (dMRI) alterations along with Short-Wavelength Automated Perimetry (SWAP) examinations of the OPs after surgery and chemotherapy over 1 year, which enabled us to evaluate chemotherapy's efficacy in an NF1 patient with an OP tumor. In this study, a 25-year-old woman with NF1 and left optic radiation (OR) glioma underwent surgery to remove the glioma. Immunohistochemistry (IHC) revealed a Pilocytic Astrocytoma (PA) WHO grade I. Post-operation chemotherapy done using nine treatment cycles of administering Temozolomide (TMZ) for 5 days every 4 weeks. Applying the region of interest (ROI) differential tractography method and SWAP four times every 3 months allowed us to follow the patient's visual acuity alterations longitudinally. The differential deterministic tractography method and statistical analyses enabled us to discover the white matter (WM) tracts anisotropy alterations over time. Furthermore, statistical analyses on the SWAP results along time illustrated possible alterations in visual acuity. Then, we could compare and associate the findings with the SWAP examinations and patient symptoms longitudinally. Statistical analyses of SWAP tests revealed a significant improvement in visual fields, and longitudinal differential tractography showed myelination and dense axonal packing in the left OR after 1 year of treatment. In this study, we examined an old hypothesis suggesting that chemotherapy is more effective than radiotherapy for NF1 patients with OP gliomas (OPGs) because of the radiation side effects on the visual field, cognition, and cerebrovascular complications. Our longitudinal clinical case study involving dMRI and SWAP on a single NF1-OPG patient showed that chemotherapy did not suppress the OP myelination over time. However, it should be noted that this is a clinical case study, and, therefore, the generalization of results is limited. Future investigations might focus on genetic-based imaging, particularly in more cases. Further, meta-analyses are recommended for giving a proper Field Of View (FOV) to researchers as a subtle clue regarding precision medicine.

NEURO-LEARN: a Solution for Collaborative Pattern Analysis of Neuroimaging Data

The development of neuroimaging instrumentation has boosted neuroscience researches. Consequently, both the fineness and the cost of data acquisition have profoundly increased, leading to the main bottleneck of this field: limited sample size and high dimensionality of neuroimaging data. Therefore, the emphasis of ideas of data pooling and research collaboration has increased over the past decade. Collaborative analysis techniques emerge as the idea developed. In this paper, we present NEURO-LEARN, a solution for collaborative pattern analysis of neuroimaging data. Its collaboration scheme consists of four parts: projects, data, analysis, and reports. While data preparation workflows defined in projects reduce the high dimensionality of neuroimaging data by collaborative computation, pooling of derived data and sharing of pattern analysis workflows along with generated reports on the Web enlarge the sample size and ensure the reliability and reproducibility of pattern analysis. Incorporating this scheme, NEURO-LEARN provides an easy-to-use Web application that allows users from different sites to share projects and processed data, perform pattern analysis, and obtain result reports. We anticipate that this solution will help neuroscientists to enlarge sample size, conquer the curse of dimensionality and conduct reproducible studies on neuroimaging data with efficiency and validity.

Using multivariate endophenotypes to identify psychophysiological mechanisms associated with polygenic scores for substance use, schizophrenia, and education attainment

BACKGROUND

To better characterize brain-based mechanisms of polygenic liability for psychopathology and psychological traits, we extended our previous report (Liu et al. Psychophysiological endophenotypes to characterize mechanisms of known schizophrenia genetic loci. Psychological Medicine, 2017), focused solely on schizophrenia, to test the association between multivariate psychophysiological candidate endophenotypes (including novel measures of θ/δ oscillatory activity) and a range of polygenic scores (PGSs), namely alcohol/cannabis/nicotine use, an updated schizophrenia PGS (containing 52 more genome-wide significant loci than the PGS used in our previous report) and educational attainment.

METHOD

A large community-based twin/family sample (N = 4893) was genome-wide genotyped and imputed. PGSs were constructed for alcohol use, regular smoking initiation, lifetime cannabis use, schizophrenia, and educational attainment. Eleven endophenotypes were assessed: visual oddball task event-related electroencephalogram (EEG) measures (target-related parietal P3 amplitude, frontal θ, and parietal δ energy/inter-trial phase clustering), band-limited resting-state EEG power, antisaccade error rate. Principal component analysis exploited covariation among endophenotypes to extract a smaller number of meaningful dimensions/components for statistical analysis.

RESULTS

Endophenotypes were heritable. PGSs showed expected intercorrelations (e.g. schizophrenia PGS correlated positively with alcohol/nicotine/cannabis PGSs). Schizophrenia PGS was negatively associated with an event-related P3/δ component [β = -0.032, nonparametric bootstrap 95% confidence interval (CI) -0.059 to -0.003]. A prefrontal control component (event-related θ/antisaccade errors) was negatively associated with alcohol (β = -0.034, 95% CI -0.063 to -0.006) and regular smoking PGSs (β = -0.032, 95% CI -0.061 to -0.005) and positively associated with educational attainment PGS (β = 0.031, 95% CI 0.003-0.058).

CONCLUSIONS

Evidence suggests that multivariate endophenotypes of decision-making (P3/δ) and cognitive/attentional control (θ/antisaccade error) relate to alcohol/nicotine, schizophrenia, and educational attainment PGSs and represent promising targets for future research.

White matter microstructure and its relation to clinical features of obsessive-compulsive disorder: findings from the ENIGMA OCD Working Group

Microstructural alterations in cortico-subcortical connections are thought to be present in obsessive-compulsive disorder (OCD). However, prior studies have yielded inconsistent findings, perhaps because small sample sizes provided insufficient power to detect subtle abnormalities. Here we investigated microstructural white matter alterations and their relation to clinical features in the largest dataset of adult and pediatric OCD to date. We analyzed diffusion tensor imaging metrics from 700 adult patients and 645 adult controls, as well as 174 pediatric patients and 144 pediatric controls across 19 sites participating in the ENIGMA OCD Working Group, in a cross-sectional case-control magnetic resonance study. We extracted measures of fractional anisotropy (FA) as main outcome, and mean diffusivity, radial diffusivity, and axial diffusivity as secondary outcomes for 25 white matter regions. We meta-analyzed patient-control group differences (Cohen's d) across sites, after adjusting for age and sex, and investigated associations with clinical characteristics. Adult OCD patients showed significant FA reduction in the sagittal stratum (d = -0.21, z = -3.21, p = 0.001) and posterior thalamic radiation (d = -0.26, z = -4.57, p < 0.0001). In the sagittal stratum, lower FA was associated with a younger age of onset (z = 2.71, p = 0.006), longer duration of illness (z = -2.086, p = 0.036), and a higher percentage of medicated patients in the cohorts studied (z = -1.98, p = 0.047). No significant association with symptom severity was found. Pediatric OCD patients did not show any detectable microstructural abnormalities compared to controls. Our findings of microstructural alterations in projection and association fibers to posterior brain regions in OCD are consistent with models emphasizing deficits in connectivity as an important feature of this disorder.

Comparison of regional brain deficit patterns in common psychiatric and neurological disorders as revealed by big data

•

RVI for MDD and AD was derived based on large meta-analytical findings.

•

RVI-MDD and AD were significantly elevated in UKBB subjects with respective illnesses.

•

There was no elevation of RVI-MDD in subjects with AD or RVI-AD in subjects with MDD.

•

RVI captures neuroanatomic deviation patterns.

•

RVI is a useful biomarker for assessing similarity to neuropsychiatric illnesses.

Neurological and psychiatric illnesses are associated with regional brain deficit patterns that bear unique signatures and capture illness-specific characteristics. The Regional Vulnerability Index (RVI) was developed to quantify brain similarity by comparing individual white matter microstructure, cortical gray matter thickness and subcortical gray matter structural volume measures with neuroanatomical deficit patterns derived from large-scale meta-analytic studies. We tested the specificity of the RVI approach for major depressive disorder (MDD) and Alzheimer’s disease (AD) in a large epidemiological sample of UK Biobank (UKBB) participants (N = 19,393; 9138 M/10,255F; age = 64.8 ± 7.4 years). Compared to controls free of neuropsychiatric disorders, participants with MDD (N = 2,248; 805 M/1443F; age = 63.4 ± 7.4) had significantly higher RVI-MDD values (t = 5.6, p = 1·10⁻⁸), but showed no detectable difference in RVI-AD (t = 2.0, p = 0.10). Subjects with dementia (N = 7; 4 M/3F; age = 68.6 ± 8.6 years) showed significant elevation in RVI-AD (t = 4.2, p = 3·10⁻⁵) but not RVI-MDD (t = 2.1, p = 0.10) compared to controls. Even within affective illnesses, participants with bipolar disorder (N = 54) and anxiety disorder (N = 773) showed no significant elevation in whole-brain RVI-MDD. Participants with Parkinson’s disease (N = 37) showed elevation in RVI-AD (t = 2.4, p = 0.01) while subjects with stroke (N = 247) showed no such elevation (t = 1.1, p = 0.3). In summary, we demonstrated elevation in RVI-MDD and RVI-AD measures in the respective illnesses with strong replicability that is relatively specific to the respective diagnoses. These neuroanatomic deviation patterns offer a useful biomarker for population-wide assessments of similarity to neuropsychiatric illnesses.

Structural Imaging in OCD

Anatomical imaging in OCD using magnetic resonance imaging (MRI) has been performed since the late 1980s. MRI research was further stimulated with the advent of automated image processing techniques such as voxel-based morphometry (VBM) and surface-based methods (e.g., FreeSurfer) which allow for detailed whole-brain data analyses. Early studies suggesting involvement of corticostriatal circuitry (particularly orbitofrontal cortex and ventral striatum) have been complemented by meta-analyses and pooled analyses indicating additional involvement of posterior brain regions, in particular parietal cortex. Recent large-scale meta-analyses from the ENIGMA consortium have revealed greater pallidum and smaller hippocampus volume in adult OCD, coupled with parietal cortical thinning. Frontal cortical thinning was only observed in medicated patients. Previous reports of symptom dimension-specific alterations were not confirmed. In paediatric OCD, thalamus enlargement has been a consistent finding. Studies investigating white matter volume (VBM) or integrity (using diffusion tensor imaging (DTI)) have shown mixed results, with recent DTI meta-analyses mainly showing involvement of posterior cortical-subcortical tracts in addition to subcortical-prefrontal connections. To which extent these abnormalities are unique to OCD or common to other psychiatric disorders is unclear, as few comparative studies have been performed. Overall, neuroanatomical alterations in OCD appear to be subtle and may vary with time, stressing the need for adequately powered longitudinal studies. Although multivariate approaches using machine learning methodologies have so far been disappointing in distinguishing individual OCD patients from healthy controls, including multimodal data in such analyses may aid in further establishing a neurobiological profile of OCD.

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

A recipe for accurate estimation of lifespan brain trajectories, distinguishing longitudinal and cohort effects

We address the problem of estimating how different parts of the brain develop and change throughout the lifespan, and how these trajectories are affected by genetic and environmental factors. Estimation of these lifespan trajectories is statistically challenging, since their shapes are typically highly nonlinear, and although true change can only be quantified by longitudinal examinations, as follow-up intervals in neuroimaging studies typically cover less than 10% of the lifespan, use of cross-sectional information is necessary. Linear mixed models (LMMs) and structural equation models (SEMs) commonly used in longitudinal analysis rely on assumptions which are typically not met with lifespan data, in particular when the data consist of observations combined from multiple studies. While LMMs require a priori specification of a polynomial functional form, SEMs do not easily handle data with unstructured time intervals between measurements. Generalized additive mixed models (GAMMs) offer an attractive alternative, and in this paper we propose various ways of formulating GAMMs for estimation of lifespan trajectories of 12 brain regions, using a large longitudinal dataset and realistic simulation experiments. We show that GAMMs are able to more accurately fit lifespan trajectories, distinguish longitudinal and cross-sectional effects, and estimate effects of genetic and environmental exposures. Finally, we discuss and contrast questions related to lifespan research which strictly require repeated measures data and questions which can be answered with a single measurement per participant, and in the latter case, which simplifying assumptions that need to be made. The examples are accompanied with R code, providing a tutorial for researchers interested in using GAMMs.

Neuroharmony: A new tool for harmonizing volumetric MRI data from unseen scanners

•We present Neuroharmony, a harmonization tool for images from unseen scanners. •We developed Neuroharmony using a total of 15,026 sMRI images. •The tool was able to reduce scanner-related bias from unseen scans. •Neuroharmony represents a significant step towards imaging-based clinical tools. •Neuroharmony is available at https://github.com/garciadias/Neuroharmony .

Meta-analysis of generalized additive models in neuroimaging studies

Analyzing data from multiple neuroimaging studies has great potential in terms of increasing statistical power, enabling detection of effects of smaller magnitude than would be possible when analyzing each study separately and also allowing to systematically investigate between-study differences. Restrictions due to privacy or proprietary data as well as more practical concerns can make it hard to share neuroimaging datasets, such that analyzing all data in a common location might be impractical or impossible. Meta-analytic methods provide a way to overcome this issue, by combining aggregated quantities like model parameters or risk ratios. Most meta-analytic tools focus on parametric statistical models, and methods for meta-analyzing semi-parametric models like generalized additive models have not been well developed. Parametric models are often not appropriate in neuroimaging, where for instance age-brain relationships may take forms that are difficult to accurately describe using such models. In this paper we introduce meta-GAM, a method for meta-analysis of generalized additive models which does not require individual participant data, and hence is suitable for increasing statistical power while upholding privacy and other regulatory concerns. We extend previous works by enabling the analysis of multiple model terms as well as multivariate smooth functions. In addition, we show how meta-analytic p-values can be computed for smooth terms. The proposed methods are shown to perform well in simulation experiments, and are demonstrated in a real data analysis on hippocampal volume and self-reported sleep quality data from the Lifebrain consortium. We argue that application of meta-GAM is especially beneficial in lifespan neuroscience and imaging genetics. The methods are implemented in an accompanying R package metagam, which is also demonstrated.

Using machine learning to quantify structural MRI neurodegeneration patterns of Alzheimer's disease into dementia score: Independent validation on 8,834 images from ADNI, AIBL, OASIS, and MIRIAD databases

Biomarkers for dementia of Alzheimer's type (DAT) are sought to facilitate accurate prediction of the disease onset, ideally predating the onset of cognitive deterioration. T1-weighted magnetic resonance imaging (MRI) is a commonly used neuroimaging modality for measuring brain structure in vivo, potentially providing information enabling the design of biomarkers for DAT. We propose a novel biomarker using structural MRI volume-based features to compute a similarity score for the individual's structural patterns relative to those observed in the DAT group. We employed ensemble-learning framework that combines structural features in most discriminative ROIs to create an aggregate measure of neurodegeneration in the brain. This classifier is trained on 423 stable normal control (NC) and 330 DAT subjects, where clinical diagnosis is likely to have the highest certainty. Independent validation on 8,834 unseen images from ADNI, AIBL, OASIS, and MIRIAD Alzheimer's disease (AD) databases showed promising potential to predict the development of DAT depending on the time-to-conversion (TTC). Classification performance on stable versus progressive mild cognitive impairment (MCI) groups achieved an AUC of 0.81 for TTC of 6 months and 0.73 for TTC of up to 7 years, achieving state-of-the-art results. The output score, indicating similarity to patterns seen in DAT, provides an intuitive measure of how closely the individual's brain features resemble the DAT group. This score can be used for assessing the presence of AD structural atrophy patterns in normal aging and MCI stages, as well as monitoring the progression of the individual's brain along with the disease course.

Deep white matter analysis (DeepWMA): Fast and consistent tractography segmentation

White matter tract segmentation, i.e. identifying tractography fibers (streamline trajectories) belonging to anatomically meaningful fiber tracts, is an essential step to enable tract quantification and visualization. In this study, we present a deep learning tractography segmentation method (DeepWMA) that allows fast and consistent identification of 54 major deep white matter fiber tracts from the whole brain. We create a large-scale training tractography dataset of 1 million labeled fiber samples, and we propose a novel 2D multi-channel feature descriptor (FiberMap) that encodes spatial coordinates of points along each fiber. We learn a convolutional neural network (CNN) fiber classification model based on FiberMap and obtain a high fiber classification accuracy of 90.99% on the training tractography data with ground truth fiber labels. Then, the method is evaluated on a test dataset of 597 diffusion MRI scans from six independently acquired populations across genders, the lifespan (1 day - 82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). We perform comparisons with two state-of-the-art tract segmentation methods. Experimental results show that our method obtains a highly consistent tract segmentation result, where on average over 99% of the fiber tracts are successfully identified across all subjects under study, most importantly, including neonates and patients with space-occupying brain tumors. We also demonstrate good generalization of the method to tractography data from multiple different fiber tracking methods. The proposed method leverages deep learning techniques and provides a fast and efficient tool for brain white matter segmentation in large diffusion MRI tractography datasets.

Large-scale mapping of cortical alterations in 22q11.2 deletion syndrome: Convergence with idiopathic psychosis and effects of deletion size

The 22q11.2 deletion (22q11DS) is a common chromosomal microdeletion and a potent risk factor for psychotic illness. Prior studies reported widespread cortical changes in 22q11DS, but were generally underpowered to characterize neuroanatomic abnormalities associated with psychosis in 22q11DS, and/or neuroanatomic effects of variability in deletion size. To address these issues, we developed the ENIGMA (Enhancing Neuro Imaging Genetics Through Meta-Analysis) 22q11.2 Working Group, representing the largest analysis of brain structural alterations in 22q11DS to date. The imaging data were collected from 10 centers worldwide, including 474 subjects with 22q11DS (age = 18.2 ± 8.6; 46.9% female) and 315 typically developing, matched controls (age = 18.0 ± 9.2; 45.9% female). Compared to controls, 22q11DS individuals showed thicker cortical gray matter overall (left/right hemispheres: Cohen's d = 0.61/0.65), but focal thickness reduction in temporal and cingulate cortex. Cortical surface area (SA), however, showed pervasive reductions in 22q11DS (left/right hemispheres: d = -1.01/-1.02). 22q11DS cases vs. controls were classified with 93.8% accuracy based on these neuroanatomic patterns. Comparison of 22q11DS-psychosis to idiopathic schizophrenia (ENIGMA-Schizophrenia Working Group) revealed significant convergence of affected brain regions, particularly in fronto-temporal cortex. Finally, cortical SA was significantly greater in 22q11DS cases with smaller 1.5 Mb deletions, relative to those with typical 3 Mb deletions. We found a robust neuroanatomic signature of 22q11DS, and the first evidence that deletion size impacts brain structure. Psychotic illness in this highly penetrant deletion was associated with similar neuroanatomic abnormalities to idiopathic schizophrenia. These consistent cross-site findings highlight the homogeneity of this single genetic etiology, and support the suitability of 22q11DS as a biological model of schizophrenia.