Genetic overlap analysis of endometriosis and asthma identifies shared loci implicating sex hormones and thyroid signalling pathways.

STUDY QUESTION: Is there a shared genetic or causal association of endometriosis with asthma or what biological mechanisms may underlie their potential relationships? SUMMARY ANSWER: Our results confirm a significant but non-causal association of endometriosis with asthma implicating shared genetic susceptibility and biological pathways in the mechanisms of the disorders, and potentially, their co-occurrence. WHAT IS KNOWN ALREADY: Some observational studies have reported a pattern of co-occurring relationship between endometriosis and asthma; however, there is conflicting evidence and the aetiology, as well as the underlying mechanisms of the relationship, remain unclear. STUDY DESIGN, SIZE, DURATION: We applied multiple statistical genetic approaches in the analysis of well-powered, genome-wide association study (GWAS) summary data to comprehensively assess the relationship of endometriosis with asthma. Endometriosis GWAS from the International Endogene Consortium (IEC, 17 054 cases and 191 858 controls) and asthma GWAS from the United Kingdom Biobank (UKB, 26 332 cases and 375 505 controls) were analysed. Additional asthma data from the Trans-National Asthma Genetic Consortium (TAGC, 19 954 cases and 107 715 controls) were utilized for replication testing. PARTICIPANTS/MATERIALS, SETTING, METHODS: We assessed single-nucleotide polymorphism (SNP)-level genetic overlap and correlation between endometriosis and asthma using SNP effect concordance analysis (SECA) and linkage disequilibrium score regression analysis (LDSC) methods, respectively. GWAS meta-analysis, colocalization (GWAS-PW), gene-based and pathway-based functional enrichment analysis methods were applied, respectively, to identify SNP loci, genomic regions, genes and biological pathways shared by endometriosis and asthma. Potential causal associations between endometriosis and asthma were assessed using Mendelian randomization (MR) methods. MAIN RESULTS AND THE ROLE OF CHANCE: SECA revealed significant concordance of SNP risk effects across the IEC endometriosis and the UKB asthma GWAS. Also, LDSC analysis found a positive and significant genetic correlation (rG = 0.16, P = 2.01 × 10-6) between the two traits. GWAS meta-analysis of the IEC endometriosis and UKB asthma GWAS identified 14 genome-wide significant (Pmeta-analysis < 5.0 × 10-8) independent loci, five of which are putatively novel. Three of these loci were consistently replicated using TAGC asthma GWAS and reinforced in colocalization and gene-based analyses. Additional shared genomic regions were identified in the colocalization analysis. MR found no evidence of a significant causal association between endometriosis and asthma. However, combining gene-based association results across the GWAS for endometriosis and asthma, we identified 17 shared genes with a genome-wide significant Fisher's combined P-value (FCPgene) <2.73 × 10-6. Additional analyses (independent gene-based analysis) replicated evidence of gene-level genetic overlap between endometriosis and asthma. Biological mechanisms including 'thyroid hormone signalling', 'abnormality of immune system physiology', 'androgen biosynthetic process' and 'brain-derived neurotrophic factor signalling pathway', among others, were significantly enriched for endometriosis and asthma in a pathway-based analysis. LARGE SCALE DATA: The GWAS for endometriosis data were sourced from the International Endogen Consortium (IEC) and can be accessed by contacting the consortium. The GWAS data for asthma are freely available online at Lee Lab (https://www.leelabsg.org/resources) and from the Trans-National Asthma Genetic Consortium (TAGC). LIMITATIONS, REASONS FOR CAUTION: Given we analysed GWAS datasets from mainly European populations, our results may not be generalizable to other ancestries. WIDER IMPLICATIONS OF THE FINDINGS: This study provides novel insights into mechanisms underpinning endometriosis and asthma, and potentially their observed relationship. Findings support a co-occurring relationship of endometriosis with asthma largely due to shared genetic components. Agents targeting 'selective androgen receptor modulators' may be therapeutically relevant in both disorders. Moreover, SNPs, loci, genes and biological pathways identified in our study provide potential targets for further investigation in endometriosis and asthma. STUDY FUNDING/COMPETING INTEREST(S): National Health and Medical Research Council (NHMRC) of Australia (241,944, 339,462, 389,927, 389,875, 389,891, 389,892, 389,938, 443,036, 442,915, 442,981, 496,610, 496,739, 552,485, 552,498, 1,026,033 and 1,050,208), Wellcome Trust (awards 076113 and 085475) and the Lundbeck Foundation (R102-A9118 and R155-2014-1724). All researchers had full independence from the funders. Authors do not have any conflict of interest.

Genome-wide association and epidemiological analyses reveal common genetic origins between uterine leiomyomata and endometriosis.

Uterine leiomyomata (UL) are the most common neoplasms of the female reproductive tract and primary cause for hysterectomy, leading to considerable morbidity and high economic burden. Here we conduct a GWAS meta-analysis in 35,474 cases and 267,505 female controls of European ancestry, identifying eight novel genome-wide significant (P < 5 × 10-8) loci, in addition to confirming 21 previously reported loci, including multiple independent signals at 10 loci. Phenotypic stratification of UL by heavy menstrual bleeding in 3409 cases and 199,171 female controls reveals genome-wide significant associations at three of the 29 UL loci: 5p15.33 (TERT), 5q35.2 (FGFR4) and 11q22.3 (ATM). Four loci identified in the meta-analysis are also associated with endometriosis risk; an epidemiological meta-analysis across 402,868 women suggests at least a doubling of risk for UL diagnosis among those with a history of endometriosis. These findings increase our understanding of genetic contribution and biology underlying UL development, and suggest overlapping genetic origins with endometriosis.

A continuum of genetic liability for minor and major depression.

The recent success of a large genome-wide association (GWA) study-analysing 130 620 major depression cases and 347 620 controls-in identifying the first single-nucleotide polymorphism (SNP) loci robustly associated with major depression in Europeans confirms that immense sample sizes are required to identify risk loci for depression. Given the phenotypic similarity between major depressive disorder (MDD) and the less severe minor depressive disorder (MiDD), we hypothesised that broadening the case definition to include MiDD may be an efficient approach to increase sample sizes in GWA studies of depression. By analysing two large twin pair cohorts, we show that minor depression and major depression lie on a single genetic continuum, with major depression being more severe but not aetiologically distinct from minor depression. Furthermore, we estimate heritabilities of 37% for minor depression, 46% for major depression and 48% for minor or major depression in a cohort of older adults (aged 50-92). However, the heritability of minor or major depression was estimated at 40% in a cohort of younger adults (aged 23-38). Moreover, two robust major depression-risk SNPs nominally associated with major depression in our Australian GWA data set produced more significant evidence for association with minor or major depression. Hence, broadening the case phenotype in GWA studies to include subthreshold definitions, such as MiDD, should facilitate the identification of additional genetic risk loci for depression.

Short telomere length is associated with impaired cognitive performance in European ancestry cohorts.

The association between telomere length (TL) dynamics on cognitive performance over the life-course is not well understood. This study meta-analyses observational and causal associations between TL and six cognitive traits, with stratifications on APOE genotype, in a Mendelian Randomization (MR) framework. Twelve European cohorts (N=17 052; mean age=59.2±8.8 years) provided results for associations between qPCR-measured TL (T/S-ratio scale) and general cognitive function, mini-mental state exam (MMSE), processing speed by digit symbol substitution test (DSST), visuospatial functioning, memory and executive functioning (STROOP). In addition, a genetic risk score (GRS) for TL including seven known genetic variants for TL was calculated, and used in associations with cognitive traits as outcomes in all cohorts. Observational analyses showed that longer telomeres were associated with better scores on DSST (ß=0.051 per s.d.-increase of TL; 95% confidence interval (CI): 0.024, 0.077; P=0.0002), and MMSE (ß=0.025; 95% CI: 0.002, 0.047; P=0.03), and faster STROOP (ß=-0.053; 95% CI: -0.087, -0.018; P=0.003). Effects for DSST were stronger in APOE ɛ4 non-carriers (ß=0.081; 95% CI: 0.045, 0.117; P=1.0 × 10-5), whereas carriers performed better in STROOP (ß=-0.074; 95% CI: -0.140, -0.009; P=0.03). Causal associations were found for STROOP only (ß=-0.598 per s.d.-increase of TL; 95% CI: -1.125, -0.072; P=0.026), with a larger effect in ɛ4-carriers (ß=-0.699; 95% CI: -1.330, -0.069; P=0.03). Two-sample replication analyses using CHARGE summary statistics showed causal effects between TL and general cognitive function and DSST, but not with STROOP. In conclusion, we suggest causal effects from longer TL on better cognitive performance, where APOE ɛ4-carriers might be at differential risk.

The future for genetic studies in reproduction.

Genetic factors contribute to risk of many common diseases affecting reproduction and fertility. In recent years, methods for genome-wide association studies (GWAS) have revolutionized gene discovery for common traits and diseases. Results of GWAS are documented in the Catalog of Published Genome-Wide Association Studies at the National Human Genome Research Institute and report over 70 publications for 32 traits and diseases associated with reproduction. These include endometriosis, uterine fibroids, age at menarche and age at menopause. Results that pass appropriate stringent levels of significance are generally well replicated in independent studies. Examples of genetic variation affecting twinning rate, infertility, endometriosis and age at menarche demonstrate that the spectrum of disease-related variants for reproductive traits is similar to most other common diseases. GWAS 'hits' provide novel insights into biological pathways and the translational value of these studies lies in discovery of novel gene targets for biomarkers, drug development and greater understanding of environmental factors contributing to disease risk. Results also show that genetic data can help define sub-types of disease and co-morbidity with other traits and diseases. To date, many studies on reproductive traits have used relatively small samples. Future genetic marker studies in large samples with detailed phenotypic and clinical information will yield new insights into disease risk, disease classification and co-morbidity for many diseases associated with reproduction and infertility.

Using genomic data to make indirect (and unauthorized) estimates of disease risk.

The number of genetic factors associated with common human traits and disease is increasing rapidly, and the general public is utilizing affordable, direct-to-consumer genetic tests. The results of these tests are often in the public domain. A combination of factors has increased the potential for the indirect estimation of an individual's risk for a particular trait. Here we explain the basic principals underlying risk estimation which allowed us to test the ability to make an indirect risk estimation from genetic data by imputing Dr. James Watson's redacted apolipoprotein E gene (APOE) information. The principles underlying risk prediction from genetic data have been well known and applied for many decades, however, the recent increase in genomic knowledge, and advances in mathematical and statistical techniques and computational power, make it relatively easy to make an accurate but indirect estimation of risk. There is a current hazard for indirect risk estimation that is relevant not only to the subject but also to individuals related to the subject; this risk will likely increase as more detailed genomic data and better computational tools become available.

Confirmation that Xq27 and Xq28 are susceptibility loci for migraine in independent pedigrees and a case-control cohort.

Investigations into migraine genetics have suggested that susceptibility loci exist on the X chromosome. These reports are supported by evidence that demonstrates male probands as having a higher proportion of affected first-degree relatives as well as the female preponderance of 3:1 that the disorder displays. We have previously implicated the Xq24-28 locus in migraine using two independent multigenerational Australian pedigrees that demonstrated excess allele sharing at the Xq24, Xq27 and Xq28 loci. Here, we expand this work to investigate a further six independent migraine pedigrees using 11 microsatellite markers spanning the Xq27­28 region. Furthermore, 11 candidate genes are investigated in an Australian case-control cohort consisting of 500 cases and 500 controls. Microsatellite analysis showed evidence of excess allele sharing to the Xq27 marker DXS8043 (LOD* 1.38 P00.005) in MF879 whilst a second independent pedigree showed excess allele sharing to DXS8061 at Xq28 (LOD* 1.5 P00.004). Furthermore, analysis of these key markers in a case control cohort showed significant association to migraine in females at the DXS8043 marker (T1 P00.009) and association with MO at DXS8061 (T1 P00.05). Further analysis of 11 key genes across these regions showed significant association of a three-marker risk haplotype in the NSDHL gene at Xq28 (P00.0082). The results of this study add further support to the presence of migraine susceptibility loci on chromosome Xq27 and Xq28 as well as point to potential candidate genes in the regions that warrant further investigation.

Genome-wide association study of major depressive disorder: new results, meta-analysis, and lessons learned.

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10); and (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

The genetic association between personality and major depression or bipolar disorder. A polygenic score analysis using genome-wide association data.

The relationship between major depressive disorder (MDD) and bipolar disorder (BD) remains controversial. Previous research has reported differences and similarities in risk factors for MDD and BD, such as predisposing personality traits. For example, high neuroticism is related to both disorders, whereas openness to experience is specific for BD. This study examined the genetic association between personality and MDD and BD by applying polygenic scores for neuroticism, extraversion, openness to experience, agreeableness and conscientiousness to both disorders. Polygenic scores reflect the weighted sum of multiple single-nucleotide polymorphism alleles associated with the trait for an individual and were based on a meta-analysis of genome-wide association studies for personality traits including 13,835 subjects. Polygenic scores were tested for MDD in the combined Genetic Association Information Network (GAIN-MDD) and MDD2000+ samples (N=8921) and for BD in the combined Systematic Treatment Enhancement Program for Bipolar Disorder and Wellcome Trust Case-Control Consortium samples (N=6329) using logistic regression analyses. At the phenotypic level, personality dimensions were associated with MDD and BD. Polygenic neuroticism scores were significantly positively associated with MDD, whereas polygenic extraversion scores were significantly positively associated with BD. The explained variance of MDD and BD, ∼0.1%, was highly comparable to the variance explained by the polygenic personality scores in the corresponding personality traits themselves (between 0.1 and 0.4%). This indicates that the proportions of variance explained in mood disorders are at the upper limit of what could have been expected. This study suggests shared genetic risk factors for neuroticism and MDD on the one hand and for extraversion and BD on the other.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction.

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

A reanalysis of 409 European-Ancestry and African American schizophrenia pedigrees reveals significant linkage to 8p23.3 with evidence of locus heterogeneity.

The detection and replication of schizophrenia risk loci can require substantial sample sizes, which has prompted various collaborative efforts for combining multiple samples. However, pooled samples may comprise sub-samples with substantial population genetic differences, including allele frequency differences. We investigated the impact of population differences via linkage reanalysis of Molecular Genetics of Schizophrenia 1 (MGS1) affected sibling-pair data, comprising two samples of distinct ancestral origin: European (EA: 263 pedigrees) and African-American (AA: 146 pedigrees). To exploit the linkage information contained within these distinct continental samples, we performed separate analyses of the individual samples, allowing for within-sample locus heterogeneity, and the pooled sample, allowing for both within-sample and between-sample heterogeneity. Significance levels, corrected for the multiple tests, were determined empirically. For all suggestive peaks, stronger linkage evidence was obtained in either the EA or AA sample than the combined sample, regardless of how heterogeneity was modeled for the latter. Notably, we report genomewide significant linkage of schizophrenia to 8p23.3 and evidence for a second, independent susceptibility locus, reaching suggestive linkage, 29 cM away on 8p21.3. We also detected suggestive linkage on chromosomes 5p13.3 and 7q36.2. Many regions showed pronounced differences in the extent of linkage between the EA and AA samples. This reanalysis highlights the potential impact of population differences upon linkage evidence in pooled data and demonstrates a useful approach for the analysis of samples drawn from distinct continental groups.

Trait components provide tools to dissect the genetic susceptibility of migraine.

The commonly used "end diagnosis" phenotype that is adopted in linkage and association studies of complex traits is likely to represent an oversimplified model of the genetic background of a disease. This is also likely to be the case for common types of migraine, for which no convincingly associated genetic variants have been reported. In headache disorders, most genetic studies have used end diagnoses of the International Headache Society (IHS) classification as phenotypes. Here, we introduce an alternative strategy; we use trait components--individual clinical symptoms of migraine--to determine affection status in genomewide linkage analyses of migraine-affected families. We identified linkage between several traits and markers on chromosome 4q24 (highest LOD score under locus heterogeneity [HLOD] 4.52), a locus we previously reported to be linked to the end diagnosis migraine with aura. The pulsation trait identified a novel locus on 17p13 (HLOD 4.65). Additionally, a trait combination phenotype (IHS full criteria) revealed a locus on 18q12 (HLOD 3.29), and the age at onset trait revealed a locus on 4q28 (HLOD 2.99). Furthermore, suggestive or nearly suggestive evidence of linkage to four additional loci was observed with the traits phonophobia (10q22) and aggravation by physical exercise (12q21, 15q14, and Xp21), and, interestingly, these loci have been linked to migraine in previous studies. Our findings suggest that the use of symptom components of migraine instead of the end diagnosis provides a useful tool in stratifying the sample for genetic studies.

Familial isolated hyperparathyroidism is linked to a 1.7 Mb region on chromosome 2p13.3-14.

BACKGROUND: Familial isolated hyperparathyroidism (FIHP) is an autosomal dominantly inherited form of primary hyperparathyroidism. Although comprising only about 1% of cases of primary hyperparathyroidism, identification and functional analysis of a causative gene for FIHP is likely to advance our understanding of parathyroid physiology and pathophysiology. METHODS: A genome-wide screen of DNA from seven pedigrees with FIHP was undertaken in order to identify a region of genetic linkage with the disorder. RESULTS: Multipoint linkage analysis identified a region of suggestive linkage (LOD score 2.68) on chromosome 2. Fine mapping with the addition of three other families revealed significant linkage adjacent to D2S2368 (maximum multipoint LOD score 3.43). Recombination events defined a 1.7 Mb region of linkage between D2S2368 and D2S358 in nine pedigrees. Sequencing of the two most likely candidate genes in this region, however, did not identify a gene for FIHP. CONCLUSIONS: We conclude that a causative gene for FIHP lies within this interval on chromosome 2. This is a major step towards eventual precise identification of a gene for FIHP, likely to be a key component in the genetic regulation of calcium homeostasis.

Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia.

The discovery of genetic factors that contribute to schizophrenia susceptibility is a key challenge in understanding the etiology of this disease. Here, we report the identification of a novel schizophrenia candidate gene on chromosome 1q32, plexin A2 (PLXNA2), in a genome-wide association study using 320 patients with schizophrenia of European descent and 325 matched controls. Over 25,000 single-nucleotide polymorphisms (SNPs) located within approximately 14,000 genes were tested. Out of 62 markers found to be associated with disease status, the most consistent finding was observed for a candidate locus on chromosome 1q32. The marker SNP rs752016 showed suggestive association with schizophrenia (odds ratio (OR) = 1.49, P = 0.006). This result was confirmed in an independent case-control sample of European Americans (combined OR = 1.38, P = 0.035) and similar genetic effects were observed in smaller subsets of Latin Americans (OR = 1.26) and Asian Americans (OR = 1.37). Supporting evidence was also obtained from two family-based collections, one of which reached statistical significance (OR = 2.2, P = 0.02). High-density SNP mapping showed that the region of association spans approximately 60 kb of the PLXNA2 gene. Eight out of 14 SNPs genotyped showed statistically significant differences between cases and controls. These results are in accordance with previous genetic findings that identified chromosome 1q32 as a candidate region for schizophrenia. PLXNA2 is a member of the transmembrane semaphorin receptor family that is involved in axonal guidance during development and may modulate neuronal plasticity and regeneration. The PLXNA2 ligand semaphorin 3A has been shown to be upregulated in the cerebellum of individuals with schizophrenia. These observations, together with the genetic results, make PLXNA2 a likely candidate for the 1q32 schizophrenia susceptibility locus.

A genome-wide scan provides evidence for loci influencing a severe heritable form of common migraine.

Migraine is a prevalent neurovascular disease with a significant genetic component. Linkage studies have so far identified migraine susceptibility loci on chromosomes 1, 4, 6, 11, 14, 19 and X. We performed a genome-wide scan of 92 Australian pedigrees phenotyped for migraine with and without aura and for a more heritable form of "severe" migraine. Multipoint non-parametric linkage analysis revealed suggestive linkage on chromosome 18p11 for the severe migraine phenotype (LOD*=2.32, P=0.0006) and chromosome 3q (LOD*=2.28, P=0.0006). Excess allele sharing was also observed at multiple different chromosomal regions, some of which overlap with, or are directly adjacent to, previously implicated migraine susceptibility regions. We have provided evidence for two loci involved in severe migraine susceptibility and conclude that dissection of the "migraine" phenotype may be helpful for identifying susceptibility genes that influence the more heritable clinical (symptom) profiles in affected pedigrees. Also, we concluded that the genetic aetiology of the common (International Headache Society) forms of the disease is probably comprised of a number of low to moderate effect susceptibility genes, perhaps acting synergistically, and this effect is not easily detected by traditional single-locus linkage analyses of large samples of affected pedigrees.

Separate and interacting effects within the catechol-O-methyltransferase (COMT) are associated with schizophrenia.

Several lines of evidence have implicated the catechol-O-methyltransferase (COMT) gene as a candidate for schizophrenia (SZ) susceptibility, not only because it encodes a key dopamine catabolic enzyme but also because it maps to the velocardiofacial syndrome region of chromosome 22q11 which has long been associated with SZ predisposition. The interest in COMT as a candidate SZ risk factor has led to numerous case-control and family-based studies, with the majority placing emphasis on examining a functional Val/Met polymorphism within this enzyme. Unfortunately, these studies have continually produced conflicting results. To assess the genetic contribution of other COMT variants to SZ susceptibility, we investigated three single-nucleotide polymorphisms (SNPs) (rs737865, rs4633, rs165599) in addition to the Val/Met variant (rs4680) in a highly selected sample of Australian Caucasian families containing 107 patients with SZ. The Val/Met and rs4633 variants showed nominally significant associations with SZ (P<0.05), although neither of the individual SNPs remained significant after adjusting for multiple testing (most significant P=0.1174). However, haplotype analyses showed strong evidence of an association; the most significant being the three-marker haplotype rs737865-rs4680-rs165599 (global P=0.0022), which spans more than 26 kb. Importantly, conditional analyses indicated the presence of two separate and interacting effects within this haplotype, irrespective of gender. In addition, our results indicate the Val/Met polymorphism is not disease-causing and is simply in strong linkage disequilibrium with a causative effect, which interacts with another as yet unidentified variant approximately 20 kb away. These results may help explain the inconsistent results reported on the Val/Met polymorphism and have important implications for future investigations into the role of COMT in SZ susceptibility.