Joshua Faskowitz

Subcortical volumetric abnormalities in bipolar disorder.

Considerable uncertainty exists about the defining brain changes associated with bipolar disorder (BD). Understanding and quantifying the sources of uncertainty can help generate novel clinical hypotheses about etiology and assist in the development of biomarkers for indexing disease progression and prognosis. Here we were interested in quantifying case–control differences in intracranial volume (ICV) and each of eight subcortical brain measures: nucleus accumbens, amygdala, caudate, hippocampus, globus pallidus, putamen, thalamus, lateral ventricles. In a large study of 1710 BD patients and 2594 healthy controls, we found consistent volumetric reductions in BD patients for mean hippocampus (Cohen’s d=−0.232; P=3.50 × 10−7) and thalamus (d=−0.148; P=4.27 × 10−3) and enlarged lateral ventricles (d=−0.260; P=3.93 × 10−5) in patients. No significant effect of age at illness onset was detected. Stratifying patients based on clinical subtype (BD type I or type II) revealed that BDI patients had significantly larger lateral ventricles and smaller hippocampus and amygdala than controls. However, when comparing BDI and BDII patients directly, we did not detect any significant differences in brain volume. This likely represents similar etiology between BD subtype classifications. Exploratory analyses revealed significantly larger thalamic volumes in patients taking lithium compared with patients not taking lithium. We detected no significant differences between BDII patients and controls in the largest such comparison to date. Findings in this study should be interpreted with caution and with careful consideration of the limitations inherent to meta-analyzed neuroimaging comparisons.

Heritability and reliability of automatically segmented human hippocampal formation subregions.

The human hippocampal formation can be divided into a set of cytoarchitecturally and functionally distinct subregions, involved in different aspects of memory formation. Neuroanatomical disruptions within these subregions are associated with several debilitating brain disorders including Alzheimer’s disease, major depression, schizophrenia, and bipolar disorder. Multi-center brain imaging consortia, such as the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) consortium, are interested in studying disease effects on these subregions, and in the genetic factors that affect them. For large-scale studies, automated extraction and subsequent genomic association studies of these hippocampal subregion measures may provide additional insight. Here, we evaluated the test–retest reliability and transplatform reliability (1.5 T versus 3 T) of the subregion segmentation module in the FreeSurfer software package using three independent cohorts of healthy adults, one young (Queensland Twins Imaging Study, N = 39), another elderly (Alzheimer’s Disease Neuroimaging Initiative, ADNI-2, N = 163) and another mixed cohort of healthy and depressed participants (Max Planck Institute, MPIP, N = 598). We also investigated agreement between the most recent version of this algorithm (v6.0) and an older version (v5.3), again using the ADNI-2 and MPIP cohorts in addition to a sample from the Netherlands Study for Depression and Anxiety (NESDA) (N = 221). Finally, we estimated the heritability (h2) of the segmented subregion volumes using the full sample of young, healthy QTIM twins (N = 728). Test–retest reliability was high for all twelve subregions in the 3 T ADNI-2 sample (intraclass correlation coefficient (ICC) = 0.70–0.97) and moderate-to-high in the 4 T QTIM sample (ICC = 0.5–0.89). Transplatform reliability was strong for eleven of the twelve subregions (ICC = 0.66–0.96); however, the hippocampal fissure was not consistently reconstructed across 1.5 T and 3 T field strengths (ICC = 0.47–0.57). Between-version agreement was moderate for the hippocampal tail, subiculum and presubiculum (ICC = 0.78–0.84; Dice Similarity Coefficient (DSC) = 0.55–0.70), and poor for all other subregions (ICC = 0.34–0.81; DSC = 0.28–0.51). All hippocampal subregion volumes were highly heritable (h2 = 0.67–0.91). Our findings indicate that eleven of the twelve human hippocampal subregions segmented using FreeSurfer version 6.0 may serve as reliable and informative quantitative phenotypes for future multi-site imaging genetics initiatives such as those of the ENIGMA consortium.

Partitioning heritability analysis reveals a shared genetic basis of brain anatomy and schizophrenia

Schizophrenia is a devastating neurodevelopmental disorder with a complex genetic etiology. Widespread cortical gray matter loss has been observed in patients and prodromal samples. However, it remains unresolved whether schizophrenia-associated cortical structure variations arise due to disease etiology or secondary to the illness. Here we address this question using a partitioning-based heritability analysis of genome-wide single-nucleotide polymorphism (SNP) and neuroimaging data from 1750 healthy individuals. We find that schizophrenia-associated genetic variants explain a significantly enriched proportion of trait heritability in eight brain phenotypes (false discovery rate=10%). In particular, intracranial volume and left superior frontal gyrus thickness exhibit significant and robust associations with schizophrenia genetic risk under varying SNP selection conditions. Cross-disorder comparison suggests that the neurogenetic architecture of schizophrenia-associated brain regions is, at least in part, shared with other psychiatric disorders. Our study highlights key neuroanatomical correlates of schizophrenia genetic risk in the general population. These may provide fundamental insights into the complex pathophysiology of the illness, and a potential link to neurocognitive deficits shaping the disorder.

Relationship of a common OXTR gene variant to brain structure and default mode network function in healthy humans

ABSTRACT A large body of research suggests that oxytocin receptor (OXTR) gene polymorphisms may influence both social behaviors and psychiatric conditions related to social deficits, such as autism spectrum disorders (ASDs), schizophrenia, and mood and anxiety disorders. However, the neural mechanism underlying these associations is still unclear. Relative to controls, patients with these psychiatric conditions show differences in brain structure, and in resting state fMRI (rs‐fMRI) signal synchronicity among default mode network (DMN) regions (also known as functional connectivity). We used a stepwise imaging genetics approach in 328 healthy young adults to test the hypothesis that 10 SNPs in OXTR are associated with differences in DMN synchronicity and structure of some of the associated brain regions. As OXTR effects may be sex‐dependent, we also tested whether our findings were modulated by sex. OXTR rs2254298 A allele carriers had significantly lower rsFC with PCC in a cluster extending from the right fronto‐insular cortex to the putamen and globus pallidus, and in bilateral dorsal anterior cingulate cortex (dACC) compared to individuals with the GG genotype; all observed effects were found only in males. Moreover, compared to the male individuals with GG genotype ofrs2254298, the male A allele carriers demonstrated significantly thinner cortical gray matter in the bilateral dACC. Our findings suggest that there may be sexually dimorphic mechanisms by which a naturally occurring variation of the OXTR gene may influence brain structure and function in DMN‐related regions implicated in neuropsychiatric disorders. HIGHLIGHTSA sex‐dependent impact of OXTR variants on the structure and function in DMN.A neural mechanism for genetically increased risk for social impairments in males.

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 months post-injury, and again 12 months later. We used tensor-based morphometry (TBM) to localize longitudinal volume expansion and reduction. We studied 21 msTBI patients (5 F, 8–18 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 months post injury. We investigated how this subgroup (TBI-slow, N = 11) differed in longitudinal regional volume changes from msTBI patients (TBI-normal, N = 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 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.

Continuous representations of brain connectivity using spatial point processes

HighlightsGeneralizes traditional connectome count matrices to spatial process of tracts.Provides fast estimator, with efficient hyper parameter tuning.Provides results showing improved reliability (as measured by ICC score).Includes demonstration analysis using analogous “degree” function.Significant differences in example analysis between sexes. Graphical abstract No Caption available. Abstract We present a continuous model for structural brain connectivity based on the Poisson point process. The model treats each streamline curve in a tractography as an observed event in connectome space, here the product space of the gray matter/white matter interfaces. We approximate the model parameter via kernel density estimation. To deal with the heavy computational burden, we develop a fast parameter estimation method by pre‐computing associated Legendre products of the data, leveraging properties of the spherical heat kernel. We show how our approach can be used to assess the quality of cortical parcellations with respect to connectivity. We further present empirical results that suggest that “discrete” connectomes derived from our model have substantially higher test‐retest reliability compared to standard methods. In this, the expanded form of this paper for journal publication, we also explore parcellation free analysis techniques that avoid the use of explicit partitions of the cortical surface altogether. We provide an analysis of sex effects on our proposed continuous representation, demonstrating the utility of this approach.

Systemic inflammation as a predictor of brain aging: Contributions of physical activity, metabolic risk, and genetic risk

&NA; Inflammatory processes may contribute to risk for Alzheimers disease (AD) and age‐related brain degeneration. Metabolic and genetic risk factors, and physical activity may, in turn, influence these inflammatory processes. Some of these risk factors are modifiable, and interact with each other. Understanding how these processes together relate to brain aging will help to inform future interventions to treat or prevent cognitive decline. We used brain magnetic resonance imaging (MRI) to scan 335 older adult humans (mean age 77.3 ± 3.4 years) who remained non‐demented for the duration of the 9‐year longitudinal study. We used structural equation modeling (SEM) in a subset of 226 adults to evaluate whether measures of baseline peripheral inflammation (serum C‐reactive protein levels; CRP), mediated the baseline contributions of genetic and metabolic risk, and physical activity, to regional cortical thickness in AD‐relevant brain regions at study year 9. We found that both baseline metabolic risk and AD risk variant apolipoprotein E &egr;4 (APOE4), modulated baseline serum CRP. Higher baseline CRP levels, in turn, predicted thinner regional cortex at year 9, and mediated an effect between higher metabolic risk and thinner cortex in those regions. A higher polygenic risk score composed of variants in immune‐associated AD risk genes (other than APOE) was associated with thinner regional cortex. However, CRP levels did not mediate this effect, suggesting that other mechanisms may be responsible for the elevated AD risk. We found interactions between genetic and environmental factors and structural brain health. Our findings support the role of metabolic risk and peripheral inflammation in age‐related brain decline. Graphical abstract Figure. No caption available. HighlightsHigher baseline blood C‐reactive protein levels were correlated with thinner posterior cingulate cortex, at study year 9.Higher metabolic risk was strongly associated with higher blood CRP levels.Higher blood CRP levels mediated an effect between higher metabolic risk and thinner cortex.A polygenic score of AD risk variants in immune‐related genes was associated with thinner cortex but not with CRP levels.

Evaluating 35 Methods to Generate Structural Connectomes Using Pairwise Classification

There is no consensus on how to construct structural brain networks from diffusion MRI. How variations in pre-processing steps affect network reliability and its ability to distinguish subjects remains opaque. In this work, we address this issue by comparing 35 structural connectome-building pipelines. We vary diffusion reconstruction models, tractography algorithms and parcellations. Next, we classify structural connectome pairs as either belonging to the same individual or not. Connectome weights and eight topological derivative measures form our feature set. For experiments, we use three test-retest datasets from the Consortium for Reliability and Reproducibility (CoRR) comprised of a total of 105 individuals. We also compare pairwise classification results to a commonly used parametric test-retest measure, Intraclass Correlation Coefficient (ICC) (Code and results are available at https://github.com/lodurality/35_methods_MICCAI_2017).

Classifying Phenotypes Based on the Community Structure of Human Brain Networks.

Human anatomical brain networks derived from the analysis of neuroimaging data are known to demonstrate modular organization. Modules, or communities, of cortical brain regions capture information about the structure of connections in the entire network. Hence, anatomical changes in network connectivity (e.g., caused by a certain disease) should translate into changes in the community structure of brain regions. This means that essential structural differences between phenotypes (e.g., healthy and diseased) should be reflected in how brain networks cluster into communities. To test this hypothesis, we propose a pipeline to classify brain networks based on their underlying community structure. We consider network partitionings into both non-overlapping and overlapping communities and introduce a distance between connectomes based on whether or not they cluster into modules similarly. We next construct a classifier that uses partitioning-based kernels to predict a phenotype from brain networks. We demonstrate the performance of the proposed approach in a task of classifying structural connectomes of healthy subjects and those with mild cognitive impairment and Alzheimer’s disease.

Blockmodels for connectome analysis

In the present work we study a family of generative network model and its applications for modeling the human connectome. We introduce a minor but novel variant of the Mixed Membership Stochastic Blockmodel and apply it and two other related model to two human connectome datasets (ADNI and a Bipolar Disorder dataset) with both control and diseased subjects. We further provide a simple generative classifier that, alongside more discriminating methods, provides evidence that blockmodels accurately summarize tractography count networks with respect to a disease classification task.