ObjectiveThe objective of this study was to aggregate data for the first genomewide association study meta-analysis of cluster headache, to identify genetic risk variants, and gain biological... Show moreObjectiveThe objective of this study was to aggregate data for the first genomewide association study meta-analysis of cluster headache, to identify genetic risk variants, and gain biological insights.MethodsA total of 4,777 cases (3,348 men and 1,429 women) with clinically diagnosed cluster headache were recruited from 10 European and 1 East Asian cohorts. We first performed an inverse-variance genomewide association meta-analysis of 4,043 cases and 21,729 controls of European ancestry. In a secondary trans-ancestry meta-analysis, we included 734 cases and 9,846 controls of East Asian ancestry. Candidate causal genes were prioritized by 5 complementary methods: expression quantitative trait loci, transcriptome-wide association, fine-mapping of causal gene sets, genetically driven DNA methylation, and effects on protein structure. Gene set and tissue enrichment analyses, genetic correlation, genetic risk score analysis, and Mendelian randomization were part of the downstream analyses.ResultsThe estimated single nucleotide polymorphism (SNP)-based heritability of cluster headache was 14.5%. We identified 9 independent signals in 7 genomewide significant loci in the primary meta-analysis, and one additional locus in the trans-ethnic meta-analysis. Five of the loci were previously known. The 20 genes prioritized as potentially causal for cluster headache showed enrichment to artery and brain tissue. Cluster headache was genetically correlated with cigarette smoking, risk-taking behavior, attention deficit hyperactivity disorder (ADHD), depression, and musculoskeletal pain. Mendelian randomization analysis indicated a causal effect of cigarette smoking intensity on cluster headache. Three of the identified loci were shared with migraine. Show less
Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE... Show moreHuman populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics. Show less
Genome-wide association analyses identify 123 susceptibility loci for migraine and implicate neurovascular mechanisms in its pathophysiology. Subtype analyses highlight risk loci specific for... Show moreGenome-wide association analyses identify 123 susceptibility loci for migraine and implicate neurovascular mechanisms in its pathophysiology. Subtype analyses highlight risk loci specific for migraine with or without aura in addition to shared risk variants.Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 controls and identified 123 loci, of which 86 are previously unknown. These loci provide an opportunity to evaluate shared and distinct genetic components in the two main migraine subtypes: migraine with aura and migraine without aura. Stratification of the risk loci using 29,679 cases with subtype information indicated three risk variants that seem specific for migraine with aura (in HMOX2, CACNA1A and MPPED2), two that seem specific for migraine without aura (near SPINK2 and near FECH) and nine that increase susceptibility for migraine regardless of subtype. The new risk loci include genes encoding recent migraine-specific drug targets, namely calcitonin gene-related peptide (CALCA/CALCB) and serotonin 1F receptor (HTR1F). Overall, genomic annotations among migraine-associated variants were enriched in both vascular and central nervous system tissue/cell types, supporting unequivocally that neurovascular mechanisms underlie migraine pathophysiology. Show less
Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions... Show moreCommon single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries. Show less
Objective Identifying common genetic variants that confergenetic risk for cluster headache. Methods We conducted a case-control study in the Dutch Leiden University Cluster headache neuro-Analysis... Show moreObjective Identifying common genetic variants that confergenetic risk for cluster headache. Methods We conducted a case-control study in the Dutch Leiden University Cluster headache neuro-Analysis program (LUCA) study population (n = 840) and unselected controls from the Netherlands Epidemiology of Obesity Study (NEO; n = 1,457). Replication was performed in a Norwegian sample of 144 cases from the Trondheim Cluster headache sample and 1,800 controls from the Nord-Trondelag Health Survey (HUNT). Gene set and tissue enrichment analyses, blood cell-derived RNA-sequencing of genes around the risk loci and linkage disequilibrium score regression were part of the downstream analyses. Results An association was found with cluster headache for 4 independent loci (r(2) < 0.1) with genomewide significance (p < 5 x 10(-8)), rs11579212 (odds ratio [OR] = 1.51, 95% confidence interval [CI] = 1.33-1.72 near RP11-815 M8.1), rs6541998 (OR = 1.53, 95% CI = 1.37-1.74 near MERTK), rs10184573 (OR = 1.43, 95% CI = 1.26-1.61 near AC093590.1), and rs2499799 (OR = 0.62, 95% CI = 0.54-0.73 near UFL1/FHL5), collectively explaining 7.2% of the variance of cluster headache. SNPs rs11579212, rs10184573, and rs976357, as proxy SNP for rs2499799 (r(2) = 1.0), replicated in the Norwegian sample (p < 0.05). Gene-based mapping yielded ASZ1 as possible fifth locus. RNA-sequencing indicated differential expression of POLR1B and TMEM87B in cluster headache patients. Interpretation This genomewide association study (GWAS) identified and replicated genetic risk loci for cluster headaches with effect sizes larger than those typically seen in complex genetic disorders. ANN NEUROL 2021 Show less
BACKGROUND: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a... Show moreBACKGROUND: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression.METHODS: We fine-mapped the classical MHC (chr6: 29.6-33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 x 10(-6) ) and a candidate threshold (1.6 x 10(-4) ).RESULTS: No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLAB*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97-0.99).CONCLUSIONS: We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes. Show less
Depression is more frequent among individuals exposed to traumatic events. Both trauma exposure and depression are heritable. However, the relationship between these traits, including the role of... Show moreDepression is more frequent among individuals exposed to traumatic events. Both trauma exposure and depression are heritable. However, the relationship between these traits, including the role of genetic risk factors, is complex and poorly understood. When modelling trauma exposure as an environmental influence on depression, both gene-environment correlations and gene-environment interactions have been observed. The UK Biobank concurrently assessed Major Depressive Disorder (MDD) and self-reported lifetime exposure to traumatic events in 126,522 genotyped individuals of European ancestry. We contrasted genetic influences on MDD stratified by reported trauma exposure (final sample size range: 24,094-92,957). The SNP-based heritability of MDD with reported trauma exposure (24%) was greater than MDD without reported trauma exposure (12%). Simulations showed that this is not confounded by the strong, positive genetic correlation observed between MDD and reported trauma exposure. We also observed that the genetic correlation between MDD and waist circumference was only significant in individuals reporting trauma exposure (r(g) = 0.24, p = 1.8 x 10(-7) versus r(g) = -0.05, p = 0.39 in individuals not reporting trauma exposure, difference p = 2.3 x 10(-4)). Our results suggest that the genetic contribution to MDD is greater when reported trauma is present, and that a complex relationship exists between reported trauma exposure, body composition, and MDD. Show less
ObjectiveTo assess whether the polygenic risk score (PRS) for migraine is associated with acute and/or prophylactic migraine treatment response.MethodsWe interviewed 2,219 unrelated patients at the... Show moreObjectiveTo assess whether the polygenic risk score (PRS) for migraine is associated with acute and/or prophylactic migraine treatment response.MethodsWe interviewed 2,219 unrelated patients at the Danish Headache Center using a semistructured interview to diagnose migraine and assess acute and prophylactic drug response. All patients were genotyped. A PRS was calculated with the linkage disequilibrium pred algorithm using summary statistics from the most recent migraine genome-wide association study comprising similar to 375,000 cases and controls. The PRS was scaled to a unit corresponding to a twofold increase in migraine risk, using 929 unrelated Danish controls as reference. The association of the PRS with treatment response was assessed by logistic regression, and the predictive power of the model by area under the curve using a case-control design with treatment response as outcome.ResultsA twofold increase in migraine risk associates with positive response to migraine-specific acute treatment (odds ratio [OR] = 1.25 [95% confidence interval (CI) = 1.05-1.49]). The association between migraine risk and migraine-specific acute treatment was replicated in an independent cohort consisting of 5,616 triptan users with prescription history (OR = 3.20 [95% CI = 1.26-8.14]). No association was found for acute treatment with non-migraine-specific weak analgesics and prophylactic treatment response.ConclusionsThe migraine PRS can significantly identify subgroups of patients with a higher-than-average likelihood of a positive response to triptans, which provides a first step toward genetics-based precision medicine in migraine. Show less
Choi, K.W.; Chen, C.Y.; Stein, M.B.; Klimentidis, Y.C.; Wang, M.J.; Koenen, K.C.; ... ; Major Depressive Disorder Working 2019
IMPORTANCE Increasing evidence shows that physical activity is associated with reduced risk for depression, pointing to a potential modifiable target for prevention. However, the causality and... Show moreIMPORTANCE Increasing evidence shows that physical activity is associated with reduced risk for depression, pointing to a potential modifiable target for prevention. However, the causality and direction of this association are not clear; physical activity may protect against depression, and/or depression may result in decreased physical activity.OBJECTIVE To examine bidirectional relationships between physical activity and depression using a genetically informed method for assessing potential causal inference.DESIGN, SETTING, AND PARTICIPANTS This 2-sample mendelian randomization (MR) used independent top genetic variants associated with 2 physical activity phenotypes-self- reported (n = 377 234) and objective accelerometer-based (n = 91 084)-and with major depressive disorder (MDD) (n = 143 265) as genetic instruments from the largest available, nonoverlapping genome-wide association studies (GWAS). GWAS were previously conducted in diverse observational cohorts, including the UK Biobank (for physical activity) and participating studies in the Psychiatric Genomics Consortium (for MDD) among adults of European ancestry. Mendelian randomization estimates from each genetic instrument were combined using inverse variance weighted meta-analysis, with alternate methods (eg, weighted median, MR Egger, MR-Pleiotropy Residual Sum and Outlier [PRESSO]) and multiple sensitivity analyses to assess horizontal pleiotropy and remove outliers. Data were analyzed from May 10 through July 31, 2018.MAIN OUTCOMES AND MEASURES MDD and physical activity.RESULTS GWAS summary data were available for a combined sample size of 611 583 adult participants. Mendelian randomization evidence suggested a protective relationship between accelerometer-based activity and MDD (odds ratio [ OR], 0.74 for MDD per 1-SD increase in mean acceleration; 95% CI, 0.59-0.92; P = .006). In contrast, there was no statistically significant relationship between MDD and accelerometer-based activity (beta = -0.08 in mean acceleration per MDD vs control status; 95% CI, -0.47 to 0.32; P = .70). Furthermore, there was no significant relationship between self-reported activity and MDD (OR, 1.28 for MDD per 1-SD increase in metabolic-equivalent minutes of reported moderate-to-vigorous activity; 95% CI, 0.57-3.37; P = .48), or between MDD and self-reported activity (beta = 0.02 per MDD in standardized metabolic-equivalent minutes of reported moderate-to-vigorous activity per MDD vs control status; 95% CI, -0.008 to 0.05; P = .15).CONCLUSIONS AND RELEVANCE Using genetic instruments identified from large-scale GWAS, robust evidence supports a protective relationship between objectively assessed-but not self-reported-physical activity and the risk for MDD. Findings point to the importance of objective measurement of physical activity in epidemiologic studies of mental health and support the hypothesis that enhancing physical activity may be an effective prevention strategy for depression. Show less
Grove, J.; Ripke, S.; Als, T.D.; Mattheisen, M.; Walters, R.K.; Won, H.; ... ; 23andMe Re 2019
Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially... Show moreAutism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD. Show less
Stress is associated with poorer physical and mental health. To improve our understanding of this link, we performed genome-wide association studies (GWAS) of depressive symptoms and genome-wide by... Show moreStress is associated with poorer physical and mental health. To improve our understanding of this link, we performed genome-wide association studies (GWAS) of depressive symptoms and genome-wide by environment interaction studies (GWEIS) of depressive symptoms and stressful life events (SLE) in two UK population-based cohorts (Generation Scotland and UK Biobank). No SNP was individually significant in either GWAS, but gene-based tests identified six genes associated with depressive symptoms in UK Biobank (DCC, ACSS3, DRD2, STAG1, FOXP2 and KYNU; p < 2.77 x 10(-6)). Two SNPs with genome-wide significant GxE effects were identified by GWEIS in Generation Scotland: rs12789145 (53-kb downstream PIWIL4; p = 4.95 x 10(-9); total SLE) and rs17070072 (intronic to ZCCHC2; p = 1.46 x 10(-8); dependent SLE). A third locus upstream CYLC2 (rs12000047 and rs12005200, p < 2.00 x 10(-8); dependent SLE) when the joint effect of the SNP main and GxE effects was considered. GWEIS gene-based tests identified: MTNR1B with GxE effect with dependent SLE in Generation Scotland; and PHF2 with the joint effect in UK Biobank (p < 2.77 x 10(-6)). Polygenic risk scores (PRSs) analyses incorporating GxE effects improved the prediction of depressive symptom scores, when using weights derived from either the UK Biobank GWAS of depressive symptoms (p = 0.01) or the PGC GWAS of major depressive disorder (p = 5.91 x 10(-3)). Using an independent sample, PRS derived using GWEIS GxE effects provided evidence of shared aetiologies between depressive symptoms and schizotypal personality, heart disease and COPD. Further such studies are required and may result in improved treatments for depression and other stress-related conditions. Show less
AbstractElectroconvulsive therapy (ECT) is the treatment of choice for severe and treatment-resistantdepression; disorder severity and unfavorable treatment outcomes are shown to be influencedby an... Show moreAbstractElectroconvulsive therapy (ECT) is the treatment of choice for severe and treatment-resistantdepression; disorder severity and unfavorable treatment outcomes are shown to be influencedby an increased genetic burden for major depression (MD). Here, we tested whether ECT assign-ment and response/nonresponse are associated with an increased genetic burden for majordepression (MD) using polygenic risk score (PRS), which summarize the contribution of disease-related common risk variants. Fifty-one psychiatric inpatients suffering from a major depressiveepisode underwent ECT. MD-PRS were calculated for these inpatients and a separatepopulation-based sample (n = 3,547 healthy; n = 426 self-reported depression) based on sum-mary statistics from the Psychiatric Genomics Consortium MDD-working group (Cases:n = 59,851; Controls: n = 113,154). MD-PRS explained a significant proportion of disease statusbetween ECT patients and healthy controls (p = .022, R2 = 1.173%); patients showed higherMD-PRS. MD-PRS in population-based depression self-reporters were intermediate betweenECT patients and controls (n.s.). Significant associations between MD-PRS and ECT response(50% reduction in Hamilton depression rating scale scores) were not observed. Our findings indi-cate that ECT cohorts show an increased genetic burden for MD and are consistent with thehypothesis that treatment-resistant MD patients represent a subgroup with an increased geneticrisk for MD. Larger samples are needed to better substantiate these findings. Show less