Martin, Nicholas G.

ABSTRACT Bipolar disorder’s (BD) clinical heterogeneity has an unresolved genetic basis. We meta-analyzed genome-wide association studies (GWAS) of 16 BD subphenotypes in 226,032 individuals from 57 cohorts (38,022 cases); 10 advanced to multivariate and multi-trait analyses. Four factors (compulsive, psychotic, dysregulated, internalizing) explained 82.8% of shared genetic variance. BD1 and BD2 loaded on distinct factors despite a high genetic correlation; 87.0% of common-factor loci were significant in neither subtype. Unipolar mania aligned with psychosis over internalizing, and was distinguishable from BD1, and rapid cycling showed heritable cross-domain liability. We identified 356 risk loci, 158 novel, including the first univariate-GWAS associations for psychosis, unipolar mania, rapid cycling and schizoaffective disorder—and 249 credible genes (89 high-confidence), 12 with approved-drug or clinical-phase annotations. Cell-type association showed a midbrain dopaminergic–GABAergic gradient along the psychotic factor. BD’s genetic architecture appears hierarchical—a general liability resolving into dimensions of course and comorbidity beyond diagnostic subtypes.

INTRODUCTION The cerebral cortex underlies our complex cognitive capabilities. Variations in human cortical surface area and thickness are associated with neurological, psychological, and behavioral traits and can be measured in vivo by magnetic resonance imaging (MRI). Studies in model organisms have identified genes that influence cortical structure, but little is known about common genetic variants that affect human cortical structure. RATIONALE To identify genetic variants associated with human cortical structure at both global and regional levels, we conducted a genome-wide association meta-analysis of brain MRI data from 51,665 individuals across 60 cohorts. We analyzed the surface area and average thickness of the whole cortex and 34 cortical regions with known functional specializations. RESULTS We identified 369 nominally genome-wide significant loci ( P < 5 × 10 −8 ) associated with cortical structure in a discovery sample of 33,992 participants of European ancestry. Of the 360 loci for which replication data were available, 241 loci influencing surface area and 66 influencing thickness remained significant after replication, with 237 loci passing multiple testing correction ( P < 8.3 × 10 −10 ; 187 influencing surface area and 50 influencing thickness). Common genetic variants explained 34% (SE = 3%) of the variation in total surface area and 26% (SE = 2%) in average thickness; surface area and thickness showed a negative genetic correlation ( r G = −0.32, SE = 0.05, P = 6.5 × 10 −12 ), which suggests that genetic influences have opposing effects on surface area and thickness. Bioinformatic analyses showed that total surface area is influenced by genetic variants that alter gene regulatory activity in neural progenitor cells during fetal development. By contrast, average thickness is influenced by active regulatory elements in adult brain samples, which may reflect processes that occur after mid-fetal development, such as myelination, branching, or pruning. When considered together, these results support the radial unit hypothesis that different developmental mechanisms promote surface area expansion and increases in thickness. To identify specific genetic influences on individual cortical regions, we controlled for global measures (total surface area or average thickness) in the regional analyses. After multiple testing correction, we identified 175 loci that influence regional surface area and 46 that influence regional thickness. Loci that affect regional surface area cluster near genes involved in the Wnt signaling pathway, which is known to influence areal identity. We observed significant positive genetic correlations and evidence of bidirectional causation of total surface area with both general cognitive functioning and educational attainment. We found additional positive genetic correlations between total surface area and Parkinson’s disease but did not find evidence of causation. Negative genetic correlations were evident between total surface area and insomnia, attention deficit hyperactivity disorder, depressive symptoms, major depressive disorder, and neuroticism. CONCLUSION This large-scale collaborative work enhances our understanding of the genetic architecture of the human cerebral cortex and its regional patterning. The highly polygenic architecture of the cortex suggests that distinct genes are involved in the development of specific cortical areas. Moreover, we find evidence that brain structure is a key phenotype along the causal pathway that leads from genetic variation to differences in general cognitive function. Identifying genetic influences on human cortical structure. ( A ) Measurement of cortical surface area and thickness from MRI. ( B ) Genomic locations of common genetic variants that influence global and regional cortical structure. ( C ) Our results support the radial unit hypothesis that the expansion of cortical surface area is driven by proliferating neural progenitor cells. ( D ) Cortical surface area shows genetic correlation with psychiatric and cognitive traits. Error bars indicate SE. IMAGE CREDITS: (A) K. COURTNEY; (C) M. R. GLASS

Abstract Eating disorders and substance use disorders frequently co‐occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa and problem alcohol use (genetic correlation [ r g ], twin‐based = 0.23‐0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome‐wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge eating, AN without binge eating, and a bulimia nervosa factor score), and eight substance‐use‐related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder and schizophrenia. Total study sample sizes per phenotype ranged from ~2400 to ~537 000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism‐based genetic correlations between eating disorder‐ and substance‐use‐related phenotypes. Significant positive genetic associations emerged between AUD and AN ( r g = 0.18; false discovery rate q = 0.0006), cannabis initiation and AN ( r g = 0.23; q < 0.0001), and cannabis initiation and AN with binge eating ( r g = 0.27; q = 0.0016). Conversely, significant negative genetic correlations were observed between three nondiagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge eating ( r gs = −0.19 to −0.23; qs < 0.04). The genetic correlation between AUD and AN was no longer significant after co‐varying for major depressive disorder loci. The patterns of association between eating disorder‐ and substance‐use‐related phenotypes highlights the potentially complex and substance‐specific relationships among these behaviors.

AbstractEpidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg = 0.23, P = 9.3 × 10−3), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10−3) and concordance in effect direction (P = 2.0 × 10−3) between the two diseases. Cross‐disease GWAS meta‐analysis highlighted 13 distinct loci associated at P ≤ 10−5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10−8, OR = 1.11, 95% CI = 1.07–1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross‐disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

Abstract Recent studies have shown that associations between common genetic variation and gene expression in breast tumors provide insight into the functional mechanisms underlying breast cancer risk loci. However, the relationship between genetic variation, gene expression, and risk of triple negative (TN) breast cancer remains largely unexplored. We performed a genome-wide expression quantitative trait locus (eQTL) analysis of TN breast cancer using 668 formalin-fixed paraffin embedded TN tumors from the Triple Negative Breast Cancer Consortium (TNBCC) to identify novel genes relevant to TN breast cancer risk and to further explore the biology underlying known TN risk loci. Cis-eQTLs were defined as correlations between single nucleotide polymorphisms (SNP) and expression of genes within 1Mb, and trans-eQTLs were defined as the remainder of all SNP-gene correlations. In a genome-wide analysis, we identified 68,012 cis-eQTLs, representing correlations between 42,193 SNPs and 2,092 genes that were significant at a 10% false discovery rate (FDR). We also identified 70 trans-eQTLs across two unique genes at the same FDR threshold. Five cis-eQTLs involved SNPs previously associated with TN breast cancer risk (p&lt;5.0x10-5) in 3,700 TN cases and 4,700 controls from the TNBCC. Three of these, located on chromosomes 5p (cis-eQTL p=4.8x10-34; risk p=1.1x10-7), 7p (cis-eQTL p=6.1x10-10; risk p=1.3x10-7), and 17q (cis-eQTL p=1.8x10-15; risk p=2.9x10-7), have previously been reported in public databases of normal tissue. The remaining two TN breast cancer cis-eQTLs, located on chromosomes 8q (cis-eQTL p=6.3x10-6, risk p=8.3x10-6) and 14q (cis-eQTL p=3.1x10-6, risk p=1.2x10-5), have not been reported as eQTLs in other tissues. In addition, 12 of 25 known TN breast cancer risk loci (PEX14, MDM4, 2q31.1, ESR1, 11q13.1, 11q24.3, 12p13.1, PTHLH, NTN4, 12q24, 19p13.1, MLK1), contained at least one significant cis-eQTL. In summary, we have identified five cis-eQTLs that may identify novel risk loci for TN breast cancer. Citation Format: Kristen S. Purrington, Drakoulis Yannoukakos, Jane Carpenter, Heli Nevanlinna, Angela Cox, Gianluca Severi, Christine Ambrosone, Amanda Ewart Toland, Andrew K. Godwin, Hiltrud Brauch, Peter A. Fasching, Penelope Miron, Jenny Chang-Claude, Nicholas G. Martin, Grant W. Montgomery, Vessela Kristensen, Hoda Anton-Culver, Paul Goodfellow, Janet E. Olson, Hugues Sicotte, Naresh Prodduturi, Daniel W. Visscher, Jeanette E. Eckel-Passow, S. Keith Anderson, Seth Slettedahl, Curtis Olswold, Xianshu Wang, V. Shane Pankratz, Susan Slager, Wei Zheng, Arto Mannermaa, Ute Hamann, Diana M. Eccles, Celine M. Vachon, Fergus J. Couch. Expression quantitative trait locus analysis of triple negative breast cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3266. doi:10.1158/1538-7445.AM2014-3266

Genetic College Many genomic elements in humans are associated with behavior, including educational attainment. In a genome-wide association study including more than 100,000 samples, Rietveld et al. (p. 1467 , published online 30 May; see the Perspective by Flint and Munafò ) looked for genes related to educational attainment in Caucasians. Small genetic effects at three loci appeared to impact educational attainment.

Aberrant connectivity is implicated in many neurological and psychiatric disorders, including Alzheimer’s disease and schizophrenia. However, other than a few disease-associated candidate genes, we know little about the degree to which genetics play a role in the brain networks; we know even less about specific genes that influence brain connections. Twin and family-based studies can generate estimates of overall genetic influences on a trait, but genome-wide association scans (GWASs) can screen the genome for specific variants influencing the brain or risk for disease. To identify the heritability of various brain connections, we scanned healthy young adult twins with high-field, high-angular resolution diffusion MRI. We adapted GWASs to screen the brain’s connectivity pattern, allowing us to discover genetic variants that affect the human brain’s wiring. The association of connectivity with the SPON1 variant at rs2618516 on chromosome 11 (11p15.2) reached connectome-wide, genome-wide significance after stringent statistical corrections were enforced, and it was replicated in an independent subsample. rs2618516 was shown to affect brain structure in an elderly population with varying degrees of dementia. Older people who carried the connectivity variant had significantly milder clinical dementia scores and lower risk of Alzheimer’s disease. As a posthoc analysis, we conducted GWASs on several organizational and topological network measures derived from the matrices to discover variants in and around genes associated with autism ( MACROD2 ), development ( NEDD4 ), and mental retardation ( UBE2A ) significantly associated with connectivity. Connectome-wide, genome-wide screening offers substantial promise to discover genes affecting brain connectivity and risk for brain diseases.