Smoller Lab at PNGU

Publications

Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects

Marshall CR, Howrigan DP, Merico D, Thiruvahindrapuram B, Wu W, Greer DS, et al. - Nature Genetics. doi: 10.1038/ng.3725

Copy number variants (CNVs) have been strongly implicated in the genetic etiology of schizophrenia (SCZ). However, genome-wide investigation of the contribution of CNV to risk has been hampered by limited sample sizes. We sought to address this obstacle by applying a centralized analysis pipeline to a SCZ cohort of 21,094 cases and 20,227 controls. A global enrichment of CNV burden was observed in cases (odds ratio (OR) = 1.11, P = 5.7 × 10-15), which persisted after excluding loci implicated in previous studies (OR = 1.07, P = 1.7 × 10-6). CNV burden was enriched for genes associated with synaptic function (OR = 1.68, P = 2.8 × 10-11) and neurobehavioral phenotypes in mouse (OR = 1.18, P = 7.3 × 10-5). Genome-wide significant evidence was obtained for eight loci, including 1q21.1, 2p16.3 (NRXN1), 3q29, 7q11.2, 15q13.3, distal 16p11.2, proximal 16p11.2 and 22q11.2. Suggestive support was found for eight additional candidate susceptibility and protective loci, which consisted predominantly of CNVs mediated by nonallelic homologous recombination.

Multidimensional heritability analysis of neuroanatomical shape

Ge T, Reuter M, Winkler AM, Holmes AJ, Lee PH, Tirrell LS, et al. - Nature Communications. doi: 10.1038/ncomms13291

In the dawning era of large-scale biomedical data, multidimensional phenotype vectors will play an increasing role in examining the genetic underpinnings of brain features, behaviour and disease. For example, shape measurements derived from brain MRI scans are multidimensional geometric descriptions of brain structure and provide an alternate class of phenotypes that remains largely unexplored in genetic studies. Here we extend the concept of heritability to multidimensional traits, and present the first comprehensive analysis of the heritability of neuroanatomical shape measurements across an ensemble of brain structures based on genome-wide SNP and MRI data from 1,320 unrelated, young and healthy individuals. We replicate our findings in an extended twin sample from the Human Connectome Project (HCP). Our results demonstrate that neuroanatomical shape can be significantly heritable, above and beyond volume, and can serve as a complementary phenotype to study the genetic determinants and clinical relevance of brain structure.

Morphometricity as a measure of the neuroanatomical signature of a trait

Sabuncu MR, Ge T, Holmes AJ, Smoller JW, Buckner RL, Fischl B, et al - Proceedings of the National Academy of Science of the United States of America (PNAS). doi: 10.1073/pnas.1604378113

Complex physiological and behavioral traits, including neurological and psychiatric disorders, often associate with distributed anatomical variation. This paper introduces a global metric, called morphometricity, as a measure of the anatomical signature of different traits. Morphometricity is defined as the proportion of phenotypic variation that can be explained by macroscopic brain morphology. We estimate morphometricity via a linear mixed-effects model that uses an anatomical similarity matrix computed based on measurements derived from structural brain MRI scans. We examined over 3,800 unique MRI scans from nine large-scale studies to estimate the morphometricity of a range of phenotypes, including clinical diagnoses such as Alzheimer's disease, and nonclinical traits such as measures of cognition. Our results demonstrate that morphometricity can provide novel insights about the neuroanatomical correlates of a diverse set of traits, revealing associations that might not be detectable through traditional statistical techniques.

Polygenic risk of Alzheimer disease is associated with early- and late-life processes

Mormino EC, Sperling RA, Holmes AJ, Buckner RL, De Jager PL, Smoller JW, et al. - Neurology. doi: 10.1212/WNL.0000000000002922

OBJECTIVE: To examine associations between aggregate genetic risk and Alzheimer disease (AD) markers in stages preceding the clinical symptoms of dementia using data from 2 large observational cohort studies.

METHODS: We computed polygenic risk scores (PGRS) using summary statistics from the International Genomics of Alzheimer's Project genome-wide association study of AD. Associations between PGRS and AD markers (cognitive decline, clinical progression, hippocampus volume, and β-amyloid) were assessed within older participants with dementia. Associations between PGRS and hippocampus volume were additionally examined within healthy younger participants (age 18-35 years).

RESULTS: Within participants without dementia, elevated PGRS was associated with worse memory (p = 0.002) and smaller hippocampus (p = 0.002) at baseline, as well as greater longitudinal cognitive decline (memory: p = 0.0005, executive function: p = 0.01) and clinical progression (p < 0.00001). High PGRS was associated with AD-like levels of β-amyloid burden as measured with florbetapir PET (p = 0.03) but did not reach statistical significance for CSF β-amyloid (p = 0.11). Within the younger group, higher PGRS was associated with smaller hippocampus volume (p = 0.05). This pattern was evident when examining a PGRS that included many loci below the genome-wide association study (GWAS)-level significance threshold (16,123 single nucleotide polymorphisms), but not when PGRS was restricted to GWAS-level significant loci (18 single nucleotide polymorphisms).

CONCLUSIONS: Effects related to common genetic risk loci distributed throughout the genome are detectable among individuals without dementia. The influence of this genetic risk may begin in early life and make an individual more susceptible to cognitive impairment in late life. Future refinement of polygenic risk scores may help identify individuals at risk for AD dementia.

Genome-wide Association Studies of Posttraumatic Stress Disorder in 2 Cohorts of US Army Soldiers

Stein MB, Chen CY, Ursano RJ, Cai T, Gelernter J, Heeringa SG, et al. - JAMA Psychiatry. doi: 10.1001/jamapsychiatry.2016.0350

IMPORTANCE: Posttraumatic stress disorder (PTSD) is a prevalent, serious public health concern, particularly in the military. The identification of genetic risk factors for PTSD may provide important insights into the biological foundation of vulnerability and comorbidity.

OBJECTIVE: To discover genetic loci associated with the lifetime risk for PTSD in 2 cohorts from the Army Study to Assess Risk and Resilience in Servicemembers (Army STARRS).

DESIGN, SETTING, AND PARTICIPANTS: Two coordinated genome-wide association studies of mental health in the US military contributed participants. The New Soldier Study (NSS) included 3167 unique participants with PTSD and 4607 trauma-exposed control individuals; the Pre/Post Deployment Study (PPDS) included 947 unique participants with PTSD and 4969 trauma-exposed controls. The NSS data were collected from February 1, 2011, to November 30, 2012; the PDDS data, from January 9 to April 30, 2012. The primary analysis compared lifetime DSM-IV PTSD cases with trauma-exposed controls without lifetime PTSD. Data were analyzed from March 18 to December 27, 2015.

MAIN OUTCOMES AND MEASURES: Association analyses for PTSD used logistic regression models within each of 3 ancestral groups (European, African, and Latino American) by study, followed by meta-analysis. Heritability and genetic correlation and pleiotropy with other psychiatric and immune-related disorders were estimated.

RESULTS: The NSS population was 80.7% male (6277 of 7774 participants; mean [SD] age, 20.9 [3.3] years); the PPDS population, 94.4% male (5583 of 5916 participants; mean [SD] age, 26.5 [6.0] years). A genome-wide significant locus was found in ANKRD55 on chromosome 5 (rs159572; odds ratio [OR], 1.62; 95% CI, 1.37-1.92; P = 2.34 × 10-8) and persisted after adjustment for cumulative trauma exposure (adjusted OR, 1.64; 95% CI, 1.39-1.95; P = 1.18 × 10-8) in the African American samples from the NSS. A genome-wide significant locus was also found in or near ZNF626 on chromosome 19 (rs11085374; OR, 0.77; 95% CI, 0.70-0.85; P = 4.59 × 10-8) in the European American samples from the NSS. Similar results were not found for either single-nucleotide polymorphism in the corresponding ancestry group from the PPDS sample, in other ancestral groups, or in transancestral meta-analyses. Single-nucleotide polymorphism-based heritability was nonsignificant, and no significant genetic correlations were observed between PTSD and 6 mental disorders or 9 immune-related disorders. Significant evidence of pleiotropy was observed between PTSD and rheumatoid arthritis and, to a lesser extent, psoriasis.

CONCLUSIONS AND RELEVANCE: In the largest genome-wide association study of PTSD to date, involving a US military sample, limited evidence of association for specific loci was found. Further efforts are needed to replicate the genome-wide significant association with ANKRD55-associated in prior research with several autoimmune and inflammatory disorders-and to clarify the nature of the genetic overlap observed between PTSD and rheumatoid arthritis and psoriasis.

Genome-wide association study of generalized anxiety symptoms in the Hispanic Community Health Study/Study of Latinos

Dunn EC, Sofer T, Gallo LC, Gogarten SM, Kerr KF, Chen CY, et al. - American Journal of Medicine Genetics Part B: Neuropsychiatric Genetics. doi: 10.1002/ajmg.b.32448

Although generalized anxiety disorder (GAD) is heritable and aggregates in families, no genomic loci associated with GAD have been reported. We aimed to discover potential loci by conducting a genome‐wide analysis of GAD symptoms in a large, population‐based sample of Hispanic/Latino adults. Data came from 12,282 participants (aged 18–74) in the Hispanic Community Health Study/Study of Latinos. Using a shortened Spielberger Trait Anxiety measure, we analyzed the following: (i) a GAD symptoms score restricted to the three items tapping diagnostic features of GAD as defined by DSM‐V; and (ii) a total trait anxiety score based on summing responses to all ten items. We first calculated the heritability due to common variants (h2SNP) and then conducted a genome‐wide association study (GWAS) of GAD symptoms. Replication was attempted in three independent Hispanic cohorts (Multi‐Ethnic Study of Atherosclerosis, Women's Health Initiative, Army STARRS). The GAD symptoms score showed evidence of modest heritability (7.2%; P = 0.03), while the total trait anxiety score did not (4.97%; P = 0.20). One genotyped SNP (rs78602344) intronic to thrombospondin 2 (THBS2) was nominally associated (P = 5.28 × 10−8) in the primary analysis adjusting for psychiatric medication use and significantly associated with the GAD symptoms score in the analysis excluding medication users (P = 4.18 × 10−8). However, meta‐analysis of the replication samples did not support this association. Although we identified a genome‐wide significant locus in this sample, we were unable to replicate this finding. Evidence for heritability was also only detected for GAD symptoms, and not the trait anxiety measure, suggesting differential genetic influences within the domain of trait anxiety.

Individual Differences in Cognitive Control Circuit Anatomy Link Sensation Seeking, Impulsivity, and Substance Use

Holmes AJ, Hollinshead MO, Roffman JL, Smoller JW, Buckner RL. - Journal of Neuroscience. doi: 10.1523/JNEUROSCI.3206-15.2016

Individuals vary widely in their tendency to seek stimulation and act impulsively, early developing traits with genetic origins. Failures to regulate these behaviors increase risk for maladaptive outcomes including substance abuse. Here, we explored the neuroanatomical correlates of sensation seeking and impulsivity in healthy young adults. Our analyses revealed links between sensation seeking and reduced cortical thickness that were preferentially localized to regions implicated in cognitive control, including anterior cingulate and middle frontal gyrus (n = 1015). These associations generalized to self-reported motor impulsivity, replicated in an independent group (n = 219), and correlated with heightened alcohol, tobacco, and caffeine use. Critically, the relations between sensation seeking and brain structure were evident in participants without a history of alcohol or tobacco use, suggesting that observed associations with anatomy are not solely a consequence of substance use. These results demonstrate that individual differences in the tendency to seek stimulation, act on impulse, and engage in substance use are correlated with the anatomical structure of cognitive control circuitry. Our findings suggest that, in healthy populations, covariation across these complex multidimensional behaviors may in part originate from a common underlying biology.

SIGNIFICANCE STATEMENT: Impaired cognitive control may result in a tendency to seek stimulation impulsively and an increased risk for maladaptive outcomes, including substance abuse. Here, we examined the structural correlates of sensation seeking and impulsivity in a large cohort of healthy young adults. Our analyses revealed links between sensation seeking and reduced cortical thickness that were preferentially localized to regions implicated in cognitive control, including anterior cingulate and middle frontal gyrus. The observed associations generalized to motor impulsivity, replicated in an independent group, and predicted heightened alcohol, tobacco, and caffeine use. These data indicate that normal variability in cognitive control system anatomy predicts sensation seeking and motor impulsivity in the healthy populations, potentially increasing risk for substance use disorders.

Genome-wide Association Study (GWAS) and Genome-Wide Environment Interaction Study (GWEIS) of Depressive Symptoms in African American and Hispanic/Latina Women

Dunn EC, Wiste A, Radmanesh F, Almli LM, Gogarten SM, Sofer T, et al. - Depression and Anxiety. doi: 10.1002/da.22484

BACKGROUND: Genome-wide association studies (GWAS) have made little progress in identifying variants linked to depression. We hypothesized that examining depressive symptoms and considering gene-environment interaction (GxE) might improve efficiency for gene discovery. We therefore conducted a GWAS and genome-wide by environment interaction study (GWEIS) of depressive symptoms.

METHODS: Using data from the SHARe cohort of the Women's Health Initiative, comprising African Americans (n = 7,179) and Hispanics/Latinas (n = 3,138), we examined genetic main effects and GxE with stressful life events and social support. We also conducted a heritability analysis using genome-wide complex trait analysis (GCTA). Replication was attempted in four independent cohorts.

RESULTS: No SNPs achieved genome-wide significance for main effects in either discovery sample. The top signals in African Americans were rs73531535 (located 20 kb from GPR139, P = 5.75 × 10(-8) ) and rs75407252 (intronic to CACNA2D3, P = 6.99 × 10(-7) ). In Hispanics/Latinas, the top signals were rs2532087 (located 27 kb from CD38, P = 2.44 × 10(-7) ) and rs4542757 (intronic to DCC, P = 7.31 × 10(-7) ). In the GEWIS with stressful life events, one interaction signal was genome-wide significant in African Americans (rs4652467; P = 4.10 × 10(-10) ; located 14 kb from CEP350). This interaction was not observed in a smaller replication cohort. Although heritability estimates for depressive symptoms and stressful life events were each less than 10%, they were strongly genetically correlated (rG = 0.95), suggesting that common variation underlying self-reported depressive symptoms and stressful life event exposure, though modest on their own, were highly overlapping in this sample.

CONCLUSIONS: Our results underscore the need for larger samples, more GEWIS, and greater investigation into genetic and environmental determinants of depressive symptoms in minorities.

Genetic variants associated with response to lithium treatment in bipolar disorder: a genome-wide association study

Hou L, Heilbronner U, Degenhardt F, Adli M, Akiyama K, Akula N, et al. - Lancet. doi: 10.1016/S0140-6736(16)00143-4

BACKGROUND: Lithium is a first-line treatment in bipolar disorder, but individual response is variable. Previous studies have suggested that lithium response is a heritable trait. However, no genetic markers of treatment response have been reproducibly identified.

METHODS: Here, we report the results of a genome-wide association study of lithium response in 2563 patients collected by 22 participating sites from the International Consortium on Lithium Genetics (ConLiGen). Data from common single nucleotide polymorphisms (SNPs) were tested for association with categorical and continuous ratings of lithium response. Lithium response was measured using a well established scale (Alda scale). Genotyped SNPs were used to generate data at more than 6 million sites, using standard genomic imputation methods. Traits were regressed against genotype dosage. Results were combined across two batches by meta-analysis.

FINDINGS: A single locus of four linked SNPs on chromosome 21 met genome-wide significance criteria for association with lithium response (rs79663003, p=1·37 × 10(-8); rs78015114, p=1·31 × 10(-8); rs74795342, p=3·31 × 10(-9); and rs75222709, p=3·50 × 10(-9)). In an independent, prospective study of 73 patients treated with lithium monotherapy for a period of up to 2 years, carriers of the response-associated alleles had a significantly lower rate of relapse than carriers of the alternate alleles (p=0·03268, hazard ratio 3·8, 95% CI 1·1-13·0).

INTERPRETATION: The response-associated region contains two genes for long, non-coding RNAs (lncRNAs), AL157359.3 and AL157359.4. LncRNAs are increasingly appreciated as important regulators of gene expression, particularly in the CNS. Confirmed biomarkers of lithium response would constitute an important step forward in the clinical management of bipolar disorder. Further studies are needed to establish the biological context and potential clinical utility of these findings.

LOCATION

Massachusetts General Hospital
Simches Research Building
185 Cambridge Street
Boston, MA 02114