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

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

Tackling the dimensions in imaging genetics with CLUB-PLS

A major challenge in imaging genetics and similar fields is to link high-dimensional data in one domain, e.g., genetic data, to high dimensional data in a second domain, e.g., brain imaging data. The standard approach in the area are mass univariate analyses across genetic factors and imaging phenotypes. That entails executing one genome-wide association study (GWAS) for each pre-defined imaging measure. Although this approach has been tremendously successful, one shortcoming is that phenotypes must be pre-defined. Consequently, effects that are not confined to pre-selected regions of interest or that reflect larger brain-wide patterns can easily be missed. In this work we introduce a Partial Least Squares (PLS)-based framework, which we term Cluster-Bootstrap PLS (CLUB-PLS), that can work with large input dimensions in both domains as well as with large sample sizes. One key factor of the framework is to use cluster bootstrap to provide robust statistics for single input features in both domains. We applied CLUB-PLS to investigating the genetic basis of surface area and cortical thickness in a sample of 33,000 subjects from the UK Biobank. We found 107 genome-wide significant locus-phenotype pairs that are linked to 386 different genes. We found that a vast majority of these loci could be technically validated at a high rate: using classic GWAS or Genome-Wide Inferred Statistics (GWIS) we found that 85 locus-phenotype pairs exceeded the genome-wide suggestive (P<1e-05) threshold.Comment: 12 pages, 4 Figures, 2 Table

Novel genetic loci associated with hippocampal volume

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer’s disease (rg=−0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness

Diverging volumetric trajectories following pediatric traumatic brain injury.

Traumatic brain injury (TBI) is a significant public health concern, and can be especially disruptive in children, derailing on-going neuronal maturation in periods critical for cognitive development. There is considerable heterogeneity in post-injury outcomes, only partially explained by injury severity. Understanding the time course of recovery, and what factors may delay or promote recovery, will aid clinicians in decision-making and provide avenues for future mechanism-based therapeutics. We examined regional changes in brain volume in a pediatric/adolescent moderate-severe TBI (msTBI) cohort, assessed at two time points. Children were first assessed 2-5&nbsp;months post-injury, and again 12&nbsp;months later. We used tensor-based morphometry (TBM) to localize longitudinal volume expansion and reduction. We studied 21 msTBI patients (5 F, 8-18&nbsp;years old) and 26 well-matched healthy control children, also assessed twice over the same interval. In a prior paper, we identified a subgroup of msTBI patients, based on interhemispheric transfer time (IHTT), with significant structural disruption of the white matter (WM) at 2-5&nbsp;months post injury. We investigated how this subgroup (TBI-slow, N&nbsp;=&nbsp;11) differed in longitudinal regional volume changes from msTBI patients (TBI-normal, N&nbsp;=&nbsp;10) with normal WM structure and function. The TBI-slow group had longitudinal decreases in brain volume in several WM clusters, including the corpus callosum and hypothalamus, while the TBI-normal group showed increased volume in WM areas. Our results show prolonged atrophy of the WM over the first 18&nbsp;months post-injury in the TBI-slow group. The TBI-normal group shows a different pattern that could indicate a return to a healthy trajectory

Neuroimaging genomics in psychiatry—a translational approach

Neuroimaging genomics is a relatively new field focused on integrating genomic and imaging data in order to investigate the mechanisms underlying brain phenotypes and neuropsychiatric disorders. While early work in neuroimaging genomics focused on mapping the associations of candidate gene variants with neuroimaging measures in small cohorts, the lack of reproducible results inspired better-powered and unbiased large-scale approaches. Notably, genome-wide association studies (GWAS) of brain imaging in thousands of individuals around the world have led to a range of promising findings. Extensions of such approaches are now addressing epigenetics, gene-gene epistasis, and gene-environment interactions, not only in brain structure, but also in brain function. Complementary developments in systems biology might facilitate the translation of findings from basic neuroscience and neuroimaging genomics to clinical practice. Here, we review recent approaches in neuroimaging genomics-we highlight the latest discoveries, discuss advantages and limitations of current approaches, and consider directions by which the field can move forward to shed light on brain disorders

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

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

Brain structural covariance networks in obsessive-compulsive disorder: a graph analysis from the ENIGMA Consortium

Brain structural covariance networks reflect covariation in morphology of different brain areas and are thought to reflect common trajectories in brain development and maturation. Large-scale investigation of structural covariance networks in obsessive-compulsive disorder (OCD) may provide clues to the pathophysiology of this neurodevelopmental disorder. Using T1-weighted MRI scans acquired from 1616 individuals with OCD and 1463 healthy controls across 37 datasets participating in the ENIGMA-OCD Working Group, we calculated intra-individual brain structural covariance networks (using the bilaterally-averaged values of 33 cortical surface areas, 33 cortical thickness values, and six subcortical volumes), in which edge weights were proportional to the similarity between two brain morphological features in terms of deviation from healthy controls (i.e. z-score transformed). Global networks were characterized using measures of network segregation (clustering and modularity), network integration (global efficiency), and their balance (small-worldness), and their community membership was assessed. Hub profiling of regional networks was undertaken using measures of betweenness, closeness, and eigenvector centrality. Individually calculated network measures were integrated across the 37 datasets using a meta-analytical approach. These network measures were summated across the network density range of K = 0.10-0.25 per participant, and were integrated across the 37 datasets using a meta-analytical approach. Compared with healthy controls, at a global level, the structural covariance networks of OCD showed lower clustering (P < 0.0001), lower modularity (P < 0.0001), and lower small-worldness (P = 0.017). Detection of community membership emphasized lower network segregation in OCD compared to healthy controls. At the regional level, there were lower (rank-transformed) centrality values in OCD for volume of caudate nucleus and thalamus, and surface area of paracentral cortex, indicative of altered distribution of brain hubs. Centrality of cingulate and orbito-frontal as well as other brain areas was associated with OCD illness duration, suggesting greater involvement of these brain areas with illness chronicity. In summary, the findings of this study, the largest brain structural covariance study of OCD to date, point to a less segregated organization of structural covariance networks in OCD, and reorganization of brain hubs. The segregation findings suggest a possible signature of altered brain morphometry in OCD, while the hub findings point to OCD-related alterations in trajectories of brain development and maturation, particularly in cingulate and orbitofrontal regions

Workflow reuse in practice: a study of neuroimaging pipeline users

Workflow reuse is a major benefit of workflow systems and shared workflow repositories, but there are barely any studies that quantify the degree of reuse of workflows or the practical barriers that may stand in the way of successful reuse. In our own work, we hypothesize that defining workflow fragments improves reuse, since end-to-end workflows may be very specific and only partially reusable by others. This paper reports on a study of the current use of workflows and workflow fragments in labs that use the LONI Pipeline, a popular workflow system used mainly for neuroimaging research that enables users to define and reuse workflow fragments. We present an overview of the benefits of workflows and workflow fragments reported by users in informal discussions. We also report on a survey of researchers in a lab that has the LONI Pipeline installed, asking them about their experiences with reuse of workflow fragments and the actual benefits they perceive. This leads to quantifiable indicators of the reuse of workflows and workflow fragments in practice. Finally, we discuss barriers to further adoption of workflow fragments and workflow reuse that motivate further work