Medial Demons Registration Localizes The Degree of Genetic Influence Over Subcortical Shape Variability: An N= 1480 Meta-Analysis

A worldwide study of subcortical shape as a marker for clinical staging in Parkinson's disease

Alterations in subcortical brain regions are linked to motor and non-motor symptoms in Parkinson's disease (PD). However, associations between clinical expression and regional morphological abnormalities of the basal ganglia, thalamus, amygdala and hippocampus are not well established. We analyzed 3D T1-weighted brain MRI and clinical data from 2525 individuals with PD and 1326 controls from 22 global sources in the ENIGMA-PD consortium. We investigated disease effects using mass univariate and multivariate models on the medial thickness of 27,120 vertices of seven bilateral subcortical structures. Shape differences were observed across all Hoehn and Yahr (HY) stages, as well as correlations with motor and cognitive symptoms. Notably, we observed incrementally thinner putamen from HY1, caudate nucleus and amygdala from HY2, hippocampus, nucleus accumbens, and thalamus from HY3, and globus pallidus from HY4-5. Subregions of the thalami were thicker in HY1 and HY2. Largely congruent patterns were associated with a longer time since diagnosis and worse motor symptoms and cognitive performance. Multivariate regression revealed patterns predictive of disease stage. These cross-sectional findings provide new insights into PD subcortical degeneration by demonstrating patterns of disease stage-specific morphology, largely consistent with ongoing degeneration.

ENIGMA-anxiety working group: Rationale for and organization of large-scale neuroimaging studies of anxiety disorders

Anxiety disorders are highly prevalent and disabling but seem particularly tractable to investigation with translational neuroscience methodologies. Neuroimaging has informed our understanding of the neurobiology of anxiety disorders, but research has been limited by small sample sizes and low statistical power, as well as heterogenous imaging methodology. The ENIGMA-Anxiety Working Group has brought together researchers from around the world, in a harmonized and coordinated effort to address these challenges and generate more robust and reproducible findings. This paper elaborates on the concepts and methods informing the work of the working group to date, and describes the initial approach of the four subgroups studying generalized anxiety disorder, panic disorder, social anxiety disorder, and specific phobia. At present, the ENIGMA-Anxiety database contains information about more than 100 unique samples, from 16 countries and 59 institutes. Future directions include examining additional imaging modalities, integrating imaging and genetic data, and collaborating with other ENIGMA working groups. The ENIGMA consortium creates synergy at the intersection of global mental health and clinical neuroscience, and the ENIGMA-Anxiety Working Group extends the promise of this approach to neuroimaging research on anxiety disorders.

GraphNet-based imaging biomarker model to explain levodopa-induced dyskinesia in Parkinson's disease

BACKGROUND AND OBJECTIVE

Levodopa-induced dyskinesia (LID) is a disabling complication of Parkinson's disease (PD). Imaging-based measurements, especially those related to the surface shape of the basal ganglia, have shown potential for explaining the severity of LID in PD. Here, we aimed to explore a novel application of the methodology to find biomarkers of LID severity in PD using regularization.

METHODS

We proposed an application of graph-constrained elastic net (GraphNet) regularization to detect surface-based shape biomarkers explaining the severity of LID and compared the approach with other conventional regularization methods. To examine the methods, we used two independent datasets, one as a training dataset to build the model, and the other dataset was used to validate the constructed model.

RESULTS

We found that the left striatum (putamen was the greatest and the caudate was second) was the most significant surface-based biomarker related to the severity of LID. Our results improved the interpretability of identified surface-based biomarkers compared to competing methods. We also found that GraphNet regularization improved prediction of the severity of LID better than the conventional regularization methods. Our model performed better in terms of root-mean-squared error and correlation coefficient between predicted and actual clinical scores.

CONCLUSION

The proposed algorithm offers an advantage of interpretable anatomical variations related to the deformation of the cortical surface. The experimental results showed that GraphNet regularization was robust to identify surface-based shape biomarkers related to both hypokinetic and hyperkinetic movement disorders.

Estimating the heritability of psychological measures in the Human Connectome Project dataset

The Human Connectome Project (HCP) is a large structural and functional MRI dataset with a rich array of behavioral and genotypic measures, as well as a biologically verified family structure. This makes it a valuable resource for investigating questions about individual differences, including questions about heritability. While its MRI data have been analyzed extensively in this regard, to our knowledge a comprehensive estimation of the heritability of the behavioral dataset has never been conducted. Using a set of behavioral measures of personality, emotion and cognition, we show that it is possible to re-identify the same individual across two testing times (fingerprinting), and to identify identical twins significantly above chance. Standard heritability estimates of 37 behavioral measures were derived from twin correlations, and machine-learning models (univariate linear model, Ridge classifier and Random Forest model) were trained to classify monozygotic twins and dizygotic twins. Correlations between the standard heritability metric and each set of model weights ranged from 0.36 to 0.7, and questionnaire-based and task-based measures did not differ significantly in their heritability. We further explored the heritability of a smaller number of latent factors extracted from the 37 measures and repeated the heritability estimation; in this case, the correlations between the standard heritability and each set of model weights were lower, ranging from 0.05 to 0.43. One specific discrepancy arose for the general intelligence factor, which all models assigned high importance, but the standard heritability calculation did not. We present a thorough investigation of the heritabilities of the behavioral measures in the HCP as a resource for other investigators, and illustrate the utility of machine-learning methods for qualitative characterization of the differential heritability across diverse measures.

Subcortical shape alterations in major depressive disorder: Findings from the ENIGMA major depressive disorder working group

Alterations in regional subcortical brain volumes have been investigated as part of the efforts of an international consortium, ENIGMA, to identify reliable neural correlates of major depressive disorder (MDD). Given that subcortical structures are comprised of distinct subfields, we sought to build significantly from prior work by precisely mapping localized MDD‐related differences in subcortical regions using shape analysis. In this meta‐analysis of subcortical shape from the ENIGMA‐MDD working group, we compared 1,781 patients with MDD and 2,953 healthy controls (CTL) on individual measures of shape metrics (thickness and surface area) on the surface of seven bilateral subcortical structures: nucleus accumbens, amygdala, caudate, hippocampus, pallidum, putamen, and thalamus. Harmonized data processing and statistical analyses were conducted locally at each site, and findings were aggregated by meta‐analysis. Relative to CTL, patients with adolescent‐onset MDD (≤ 21 years) had lower thickness and surface area of the subiculum, cornu ammonis (CA) 1 of the hippocampus and basolateral amygdala (Cohen's d = −0.164 to −0.180). Relative to first‐episode MDD, recurrent MDD patients had lower thickness and surface area in the CA1 of the hippocampus and the basolateral amygdala (Cohen's d = −0.173 to −0.184). Our results suggest that previously reported MDD‐associated volumetric differences may be localized to specific subfields of these structures that have been shown to be sensitive to the effects of stress, with important implications for mapping treatments to patients based on specific neural targets and key clinical features.

Amygdala Nuclei Volume and Shape in Military Veterans With Posttraumatic Stress Disorder

BACKGROUND

The amygdala is a subcortical structure involved in socioemotional and associative fear learning processes relevant for understanding the mechanisms of posttraumatic stress disorder (PTSD). Research in animals indicates that the amygdala is a heterogeneous structure in which the basolateral and centromedial divisions are susceptible to stress. While the amygdala complex is implicated in the pathophysiology of PTSD, little is known about the specific contributions of the individual nuclei that constitute the amygdala complex.

METHODS

Military veterans (n = 355), including military veterans with PTSD (n = 149) and trauma-exposed control subjects without PTSD (n = 206), underwent high-resolution T1-weighted anatomical scans. Automated FreeSurfer segmentation of the amygdala yielded 9 structures: basal, lateral, accessory basal, anterior amygdaloid, and central, medial, cortical, and paralaminar nuclei, along with the corticoamygdaloid transition zone. Subregional volumes were compared between groups using ordinary-least-squares regression with relevant demographic and clinical regressors followed by 3-dimensional shape analysis of whole amygdala.

RESULTS

PTSD was associated with smaller left and right lateral and paralaminar nuclei, but with larger left and right central, medial, and cortical nuclei (p < .05, false discovery rate corrected). Shape analyses revealed lower radial distance in anterior bilateral amygdala and lower Jacobian determinant in posterior bilateral amygdala in PTSD compared with control subjects.

CONCLUSIONS

Alterations in select amygdala subnuclear volumes and regional shape distortions are associated with PTSD in military veterans. Volume differences of the lateral nucleus and the centromedial complex associated with PTSD demonstrate a subregion-specific pattern that is consistent with their functional roles in fear learning and fear expression behaviors.

Mapping abnormal subcortical neurodevelopment in a cohort of Thai children with HIV

Alterations in subcortical brain structures have been reported in adults with HIV and, to a lesser extent, pediatric cohorts. The extent of longitudinal structural abnormalities in children with perinatal HIV infection (PaHIV) remains unclear. We modeled subcortical morphometry from whole brain structural magnetic resonance imaging (1.5 T) scans of 43 Thai children with PaHIV (baseline age = 11.09±2.36 years) and 50 HIV- children (11.26±2.80 years) using volumetric and surface-based shape analyses. The PaHIV sample were randomized to initiate combination antiretroviral treatment (cART) when CD4 counts were 15-24% (immediate: n = 22) or when CD4 < 15% (deferred: n = 21). Follow-up scans were acquired approximately 52 weeks after baseline. Volumetric and shape descriptors capturing local thickness and surface area dilation were defined for the bilateral accumbens, amygdala, putamen, pallidum, thalamus, caudate, and hippocampus. Regression models adjusting for clinical and demographic variables examined between and within group differences in morphometry associated with HIV. We assessed whether baseline CD4 count and cART status or timing associated with brain maturation within the PaHIV group. All models were adjusted for multiple comparisons using the false discovery rate. A pallidal subregion was significantly thinner in children with PaHIV. Regional thickness, surface area, and volume of the pallidum was associated with CD4 count in children with PaHIV. Longitudinal morphometry was not associated with HIV or cART status or timing, however, the trajectory of the left pallidum volume was positively associated with baseline CD4 count. Our findings corroborate reports in adult cohorts demonstrating a high predilection for HIV-mediated abnormalities in the basal ganglia, but suggest the effect of stable PaHIV infection on morphological aspects of brain development may be subtle.

Machine Learning for Large-Scale Quality Control of 3D Shape Models in Neuroimaging

As very large studies of complex neuroimaging phenotypes become more common, human quality assessment of MRI-derived data remains one of the last major bottlenecks. Few attempts have so far been made to address this issue with machine learning. In this work, we optimize predictive models of quality for meshes representing deep brain structure shapes. We use standard vertex-wise and global shape features computed homologously across 19 cohorts and over 7500 human-rated subjects, training kernelized Support Vector Machine and Gradient Boosted Decision Trees classifiers to detect meshes of failing quality. Our models generalize across datasets and diseases, reducing human workload by 30-70%, or equivalently hundreds of human rater hours for datasets of comparable size, with recall rates approaching inter-rater reliability.

Mapping 22q11.2 Gene Dosage Effects on Brain Morphometry

Reciprocal chromosomal rearrangements at the 22q11.2 locus are associated with elevated risk of neurodevelopmental disorders. The 22q11.2 deletion confers the highest known genetic risk for schizophrenia, but a duplication in the same region is strongly associated with autism and is less common in schizophrenia cases than in the general population. Here we conducted the first study of 22q11.2 gene dosage effects on brain structure in a sample of 143 human subjects: 66 with 22q11.2 deletions (22q-del; 32 males), 21 with 22q11.2 duplications (22q-dup; 14 males), and 56 age- and sex-matched controls (31 males). 22q11.2 gene dosage varied positively with intracranial volume, gray and white matter volume, and cortical surface area (deletion < control < duplication). In contrast, gene dosage varied negatively with mean cortical thickness (deletion > control > duplication). Widespread differences were observed for cortical surface area with more localized effects on cortical thickness. These diametric patterns extended into subcortical regions: 22q-dup carriers had a significantly larger right hippocampus, on average, but lower right caudate and corpus callosum volume, relative to 22q-del carriers. Novel subcortical shape analysis revealed greater radial distance (thickness) of the right amygdala and left thalamus, and localized increases and decreases in subregions of the caudate, putamen, and hippocampus in 22q-dup relative to 22q-del carriers. This study provides the first evidence that 22q11.2 is a genomic region associated with gene-dose-dependent brain phenotypes. Pervasive effects on cortical surface area imply that this copy number variant affects brain structure early in the course of development.SIGNIFICANCE STATEMENT Probing naturally occurring reciprocal copy number variation in the genome may help us understand mechanisms underlying deviations from typical brain and cognitive development. The 22q11.2 genomic region is particularly susceptible to chromosomal rearrangements and contains many genes crucial for neuronal development and migration. Not surprisingly, reciprocal genomic imbalances at this locus confer some of the highest known genetic risks for developmental neuropsychiatric disorders. Here we provide the first evidence that brain morphology differs meaningfully as a function of reciprocal genomic variation at the 22q11.2 locus. Cortical thickness and surface area were affected in opposite directions with more widespread effects of gene dosage on cortical surface area.

APPROXIMATING PRINCIPAL GENETIC COMPONENTS OF SUBCORTICAL SHAPE

Optimal representations of the genetic structure underlying complex neuroimaging phenotypes lie at the heart of our quest to discover the genetic code of the brain. Here, we suggest a strategy for achieving such a representation by decomposing the genetic covariance matrix of complex phenotypes into maximally heritable and genetically independent components. We show that such a representation can be approximated well with eigenvectors of the genetic covariance based on a large family study. Using 520 twin pairs from the QTIM dataset, we estimate 500 principal genetic components of 54,000 vertex-wise shape features representing fourteen subcortical regions. We show that our features maintain their desired properties in practice. Further, the genetic components are found to be significantly associated with the CLU and PICALM genes in an unrelated Alzheimer's Disease (AD) dataset. The same genes are not significantly associated with other volume and shape measures in this dataset.

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

In this review, we discuss recent work by the ENIGMA Consortium (http://enigma.ini.usc.edu) - a global alliance of over 500 scientists spread across 200 institutions in 35 countries collectively analyzing brain imaging, clinical, and genetic data. Initially formed to detect genetic influences on brain measures, ENIGMA has grown to over 30 working groups studying 12 major brain diseases by pooling and comparing brain data. In some of the largest neuroimaging studies to date - of schizophrenia and major depression - ENIGMA has found replicable disease effects on the brain that are consistent worldwide, as well as factors that modulate disease effects. In partnership with other consortia including ADNI, CHARGE, IMAGEN and others1, ENIGMA's genomic screens - now numbering over 30,000 MRI scans - have revealed at least 8 genetic loci that affect brain volumes. Downstream of gene findings, ENIGMA has revealed how these individual variants - and genetic variants in general - may affect both the brain and risk for a range of diseases. The ENIGMA consortium is discovering factors that consistently affect brain structure and function that will serve as future predictors linking individual brain scans and genomic data. It is generating vast pools of normative data on brain measures - from tens of thousands of people - that may help detect deviations from normal development or aging in specific groups of subjects. We discuss challenges and opportunities in applying these predictors to individual subjects and new cohorts, as well as lessons we have learned in ENIGMA's efforts so far.

Heritability of the shape of subcortical brain structures in the general population

The volumes of subcortical brain structures are highly heritable, but genetic underpinnings of their shape remain relatively obscure. Here we determine the relative contribution of genetic factors to individual variation in the shape of seven bilateral subcortical structures: the nucleus accumbens, amygdala, caudate, hippocampus, pallidum, putamen and thalamus. In 3,686 unrelated individuals aged between 45 and 98 years, brain magnetic resonance imaging and genotyping was performed. The maximal heritability of shape varies from 32.7 to 53.3% across the subcortical structures. Genetic contributions to shape extend beyond influences on intracranial volume and the gross volume of the respective structure. The regional variance in heritability was related to the reliability of the measurements, but could not be accounted for by technical factors only. These findings could be replicated in an independent sample of 1,040 twins. Differences in genetic contributions within a single region reveal the value of refined brain maps to appreciate the genetic complexity of brain structures.

Volumetric and shape analyses of subcortical structures in United States service members with mild traumatic brain injury

Mild traumatic brain injury (mTBI) is a significant health concern. The majority who sustain mTBI recover, although ~20 % continue to experience symptoms that can interfere with quality of life. Accordingly, there is a critical need to improve diagnosis, prognostic accuracy, and monitoring (recovery trajectory over time) of mTBI. Volumetric magnetic resonance imaging (MRI) has been successfully utilized to examine TBI. One promising improvement over standard volumetric approaches is to analyze high-dimensional shape characteristics of brain structures. In this study, subcortical shape and volume in 76 Service Members with mTBI was compared to 59 Service Members with orthopedic injury (OI) and 17 with post-traumatic stress disorder (PTSD) only. FreeSurfer was used to quantify structures from T1-weighted 3 T MRI data. Radial distance (RD) and Jacobian determinant (JD) were defined vertex-wise on parametric mesh-representations of subcortical structures. Linear regression was used to model associations between morphometry (volume and shape), TBI status, and time since injury (TSI) correcting for age, sex, intracranial volume, and level of education. Volumetric data was not significantly different between the groups. JD was significantly increased in the accumbens and caudate and significantly reduced in the thalamus of mTBI participants. Additional significant associations were noted between RD of the amygdala and TSI. Positive trend-level associations between TSI and the amygdala and accumbens were observed, while a negative association was observed for third ventricle. Our findings may aid in the initial diagnosis of mTBI, provide biological targets for functional examination, and elucidate regions that may continue remodeling after injury.

Effect of Electroconvulsive Therapy on Striatal Morphometry in Major Depressive Disorder

Patients with major depression show reductions in striatal and paleostriatal volumes. The functional integrity and connectivity of these regions are also shown to change with antidepressant response. Electroconvulsive therapy (ECT) is a robust and rapidly acting treatment for severe depression. However, whether morphological changes in the dorsal and ventral striatum/pallidum relate to or predict therapeutic response to ECT is unknown. Using structural MRI, we assessed cross-sectional effects of diagnosis and longitudinal effects of ECT for volume and surface-based shape metrics of the caudate, putamen, pallidum, and nucleus accumbens in 53 depressed patients (mean age: 44.1 years, 13.8 SD; 52% female) and 33 healthy controls (mean age: 39.3 years, 12.4 SD; 57% female). Patients were assessed before ECT, after their second ECT, and after completing an ECT treatment index. Controls were evaluated at two time points. Support vector machines determined whether morphometric measures at baseline predicted ECT-related clinical response. Patients showed smaller baseline accumbens and pallidal volumes than controls (P<0.05). Increases in left putamen volume (P<0.03) occurred with ECT. Global increases in accumbens volume and local changes in pallidum and caudate volume occurred in patients defined as treatment responders. Morphometric changes were absent across time in controls. Baseline volume and shape metrics predicted overall response to ECT with up to 89% accuracy. Results support that ECT elicits structural plasticity in the dorsal and ventral striatum/pallidum. The morphometry of these structures, forming key components of limbic-cortical-striatal-pallidal-thalamic circuitry involved in mood and emotional regulation, may determine patients likely to benefit from treatment.

Mapping abnormal subcortical brain morphometry in an elderly HIV+ cohort

Over 50% of HIV + individuals exhibit neurocognitive impairment and subcortical atrophy, but the profile of brain abnormalities associated with HIV is still poorly understood. Using surface-based shape analyses, we mapped the 3D profile of subcortical morphometry in 63 elderly HIV + participants and 31 uninfected controls. The thalamus, caudate, putamen, pallidum, hippocampus, amygdala, brainstem, accumbens, callosum and ventricles were segmented from high-resolution MRIs. To investigate shape-based morphometry, we analyzed the Jacobian determinant (JD) and radial distances (RD) defined on each region's surfaces. We also investigated effects of nadir CD4 + T-cell counts, viral load, time since diagnosis (TSD) and cognition on subcortical morphology. Lastly, we explored whether HIV + participants were distinguishable from unaffected controls in a machine learning context. All shape and volume features were included in a random forest (RF) model. The model was validated with 2-fold cross-validation. Volumes of HIV + participants' bilateral thalamus, left pallidum, left putamen and callosum were significantly reduced while ventricular spaces were enlarged. Significant shape variation was associated with HIV status, TSD and the Wechsler adult intelligence scale. HIV + people had diffuse atrophy, particularly in the caudate, putamen, hippocampus and thalamus. Unexpectedly, extended TSD was associated with increased thickness of the anterior right pallidum. In the classification of HIV + participants vs. controls, our RF model attained an area under the curve of 72%.

•

We model subcortical morphometry of elderly HIV + participants.

•

We explore classifying HIV status based on shape and volume of brain regions.

•

Morphometry of brain regions was associated with infection status and duration.

•

HIV status was classifiable with 72% accuracy in morphometry-based classifiers.