Association between endometriosis and the interleukin 1A (IL1A) locus

STUDY QUESTION Are single-nucleotide polymorphisms (SNPs) at the interleukin 1A (IL1A) gene locus associated with endometriosis risk? SUMMARY ANSWER We found evidence for strong association between IL1A SNPs and endometriosis risk. WHAT IS KNOWN ALREADY Genetic factors contribute substantially to the complex aetiology of endometriosis and the disease has an estimated heritability of ∼51%. We, and others, have conducted genome-wide association (GWA) studies for endometriosis, which identified a total of nine independent risk loci. Recently, two small Japanese studies reported eight SNPs (rs6542095, rs11677416, rs3783550, rs3783525, rs3783553, rs2856836, rs1304037 and rs17561) at the IL1A gene locus as suggestively associated with endometriosis risk. There is also evidence of a link between inflammation and endometriosis. STUDY DESIGN, SIZE, DURATION We sought to further investigate the eight IL1A SNPs for association with endometriosis using an independent sample of 3908 endometriosis cases and 8568 controls of European and Japanese ancestry. The study was conducted between October 2013 and July 2014. PARTICIPANTS/MATERIALS, SETTING, METHODS By leveraging GWA data from our previous multi-ethnic GWA meta-analysis for endometriosis, we imputed variants in the IL1A region, using a recent 1000 Genomes reference panel. After combining summary statistics for the eight SNPs from our European and Japanese imputed data with the published results, a fixed-effect meta-analysis was performed. An additional meta-analysis restricted to endometriosis cases with moderate-to-severe (revised American Fertility Society stage 3 or 4) disease versus controls was also performed. MAIN RESULTS AND THE ROLE OF CHANCE All eight IL1A SNPs successfully replicated at P < 0.014 in the European imputed data with concordant direction and similar size to the effects reported in the original Japanese studies. Of these, three SNPs (rs6542095, rs3783550 and rs3783525) also showed association with endometriosis at a nominal P < 0.05 in our independent Japanese sample. Fixed-effect meta-analysis of the eight SNPs for moderate-to-severe endometriosis produced a genome-wide significant association for rs6542095 (odds ratio = 1.21; 95% confidence interval = 1.13–1.29; P = 3.43 × 10−8). LIMITATIONS, REASONS FOR CAUTION The meta-analysis for moderate-to-severe endometriosis included results of moderate-to-severe endometriosis cases from our European data sets and all endometriosis cases from the Japanese data sets, as disease stage information was not available for endometriosis cases in the Japanese data sets. WIDER IMPLICATIONS OF THE FINDINGS SNP rs6542095 is located ∼2.3 kb downstream of the IL1A gene and ∼6.9 kb upstream of cytoskeleton-associated protein 2-like (CKAP2L) gene. The IL1A gene encodes the IL1a protein, a member of the interleukin 1 cytokine family which is involved in various immune responses and inflammatory processes. These results provide important replication in an independent Japanese sample and, for the first time, association of the IL1A locus in endometriosis patients of European ancestry. SNPs within the IL1A locus may regulate other genes, but if IL1A is the target, our results provide supporting evidence for a link between inflammatory responses and the pathogenesis of endometriosis. STUDY FUNDING/COMPETING INTEREST(S) The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study.

Genetic burden associated with varying degrees of disease severity in endometriosis

Endometriosis is primarily characterized by the presence of tissue resembling endometrium outside the uterine cavity and is usually diagnosed by laparoscopy. The most commonly used classification of disease, the revised American Fertility Society (rAFS) system to grade endometriosis into different stages based on disease severity (I to IV), has been questioned as it does not correlate well with underlying symptoms, posing issues in diagnosis and choice of treatment. Using two independent European genome-wide association (GWA) datasets and top-level classification of the endometriosis cases based on rAFS [minimal or mild (Stage A) and moderate-to-severe (Stage B) disease], we previously showed that Stage B endometriosis has greater contribution of common genetic variation to its aetiology than Stage A disease. Herein, we extend our previous analysis to four endometriosis stages [minimal (Stage I), mild (Stage II), moderate (Stage III) and severe (Stage IV) disease] based on the rAFS classification system and compared the genetic burden across stages. Our results indicate that genetic burden increases from minimal to severe endometriosis. For the minimal disease, genetic factors may contribute to a lesser extent than other disease categories. Mild and moderate endometriosis appeared genetically similar, making it difficult to tease them apart. Consistent with our previous reports, moderate and severe endometriosis showed greater genetic burden than minimal or mild disease. Overall, our results provide new insights into the genetic architecture of endometriosis and further investigation in larger samples may help to understand better the aetiology of varying degrees of endometriosis, enabling improved diagnostic and treatment modalities.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around 1 in 7 people worldwide, but its molecular mechanisms remain poorly understood. Some debate exists over whether migraine is a disease of vascular dysfunction, or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we performed the largest genetic study of migraine to date, comprising 59,674 cases and 316,078 controls from 22 GWA studies. We identified 45 independent single nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 x 10-8) that map to 38 distinct genomic loci, including 28 loci not previously reported and the first locus identified on chromosome X. Furthermore, a subset analysis for migraine without aura (MO) identified seven of the same loci as from the full sample, whereas no loci reached genome-wide significance in the migraine with aura (MA) subset. In subsequent computational analyzes, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

A mega-analysis of genome-wide association studies for major depressive disorder

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 x 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 x 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs

Most psychiatric disorders are moderately to highly heritable. The degree to which genetic variation is unique to individual disorders or shared across disorders is unclear. To examine shared genetic etiology, we use genome-wide genotype data from the Psychiatric Genomics Consortium (PGC) for cases and controls in schizophrenia, bipolar disorder, major depressive disorder, autism spectrum disorders (ASD) and attention-deficit/hyperactivity disorder (ADHD). We apply univariate and bivariate methods for the estimation of genetic variation within and covariation between disorders. SNPs explained 17-29% of the variance in liability. The genetic correlation calculated using common SNPs was high between schizophrenia and bipolar disorder (0.68 +/- 0.04 s.e.), moderate between schizophrenia and major depressive disorder (0.43 +/- 0.06 s.e.), bipolar disorder and major depressive disorder (0.47 +/- 0.06 s.e.), and ADHD and major depressive disorder (0.32 +/- 0.07 s.e.), low between schizophrenia and ASD (0.16 +/- 0.06 s.e.) and non-significant for other pairs of disorders as well as between psychiatric disorders and the negative control of Crohn's disease. This empirical evidence of shared genetic etiology for psychiatric disorders can inform nosology and encourages the investigation of common pathophysiologies for related disorders.

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.

Genome-wide association meta-analysis identifies new endometriosis risk loci

We conducted a genome-wide association meta-analysis of 4,604 endometriosis cases and 9,393 controls of Japanese and European ancestry. We show that rs12700667 on chromosome 7p15.2, previously found to associate with disease in Europeans, replicates in Japanese (P = 3.6 x 10(-3)), and we confirm association of rs7521902 at 1p36.12 near WNT4. In addition, we establish an association of rs13394619 in GREB1 at 2p25.1 with endometriosis and identify a newly associated locus at 12q22 near VEZT (rs10859871). Excluding cases of European ancestry of minimal or unknown severity, we identified additional previously unknown loci at 2p14 (rs4141819), 6p22.3 (rs7739264) and 9p21.3 (rs1537377). All seven SNP effects were replicated in an independent cohort and associated at P <5 x 10(-8) in a combined analysis. Finally, we found a significant overlap in polygenic risk for endometriosis between the genome-wide association cohorts of European and Japanese descent (P = 8.8 x 10(-11)), indicating that many weakly associated SNPs represent true endometriosis risk loci and that risk prediction and future targeted disease therapy may be transferred across these populations.

Loci affecting gamma-glutamyl transferase in adults and adolescents show age x SNP interaction and cardiometabolic disease associations

Serum gamma-glutamyl transferase (GGT) activity is a marker of liver disease which is also prospectively associated with the risk of all-cause mortality, cardiovascular disease, type 2 diabetes and cancers. We have discovered novel loci affecting GGT in a genome-wide association study (rs1497406 in an intergenic region of chromosome 1, P = 3.9 x 10(-8); rs944002 in C14orf73 on chromosome 14, P = 4.7 x 10(-13); rs340005 in RORA on chromosome 15, P = 2.4 x 10(-8)), and a highly significant heterogeneity between adult and adolescent results at the GGT1 locus on chromosome 22 (maximum P(HET) = 5.6 x 10(-12) at rs6519520). Pathway analysis of significant and suggestive single-nucleotide polymorphism associations showed significant overlap between genes affecting GGT and those affecting common metabolic and inflammatory diseases, and identified the hepatic nuclear factor (HNF) family as controllers of a network of genes affecting GGT. Our results reinforce the disease associations of GGT and demonstrate that control by the GGT1 locus varies with age.

Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis

Endometriosis is a common gynecological disease associated with pelvic pain and subfertility. We conducted a genome-wide association study (GWAS) in 3,194 individuals with surgically confirmed endometriosis (cases) and 7,060 controls from Australia and the UK. Polygenic predictive modeling showed significantly increased genetic loading among 1,364 cases with moderate to severe endometriosis. The strongest association signal was on 7p15.2 (rs12700667) for 'all' endometriosis (P = 2.6 x 10(-)(7), odds ratio (OR) = 1.22, 95% CI 1.13-1.32) and for moderate to severe disease (P = 1.5 x 10(-)(9), OR = 1.38, 95% CI 1.24-1.53). We replicated rs12700667 in an independent cohort from the United States of 2,392 self-reported, surgically confirmed endometriosis cases and 2,271 controls (P = 1.2 x 10(-)(3), OR = 1.17, 95% CI 1.06-1.28), resulting in a genome-wide significant P value of 1.4 x 10(-)(9) (OR = 1.20, 95% CI 1.13-1.27) for 'all' endometriosis in our combined datasets of 5,586 cases and 9,331 controls. rs12700667 is located in an intergenic region upstream of the plausible candidate genes NFE2L3 and HOXA10.

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, approximately 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.

A quantitative-trait genome-wide association study of alcoholism risk in the community: findings and implications

BACKGROUND Given moderately strong genetic contributions to variation in alcoholism and heaviness of drinking (50% to 60% heritability) with high correlation of genetic influences, we have conducted a quantitative trait genome-wide association study (GWAS) for phenotypes related to alcohol use and dependence. METHODS Diagnostic interview and blood/buccal samples were obtained from sibships ascertained through the Australian Twin Registry. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed with 8754 individuals (2062 alcohol-dependent cases) selected for informativeness for alcohol use disorder and associated quantitative traits. Family-based association tests were performed for alcohol dependence, dependence factor score, and heaviness of drinking factor score, with confirmatory case-population control comparisons using an unassessed population control series of 3393 Australians with genome-wide SNP data. RESULTS No findings reached genome-wide significance (p = 8.4 x 10(-8) for this study), with lowest p value for primary phenotypes of 1.2 x 10(-7). Convergent findings for quantitative consumption and diagnostic and quantitative dependence measures suggest possible roles for a transmembrane protein gene (TMEM108) and for ANKS1A. The major finding, however, was small effect sizes estimated for individual SNPs, suggesting that hundreds of genetic variants make modest contributions (1/4% of variance or less) to alcohol dependence risk. CONCLUSIONS We conclude that: - 1) meta-analyses of consumption data may contribute usefully to gene discovery; - 2) translation of human alcoholism GWAS results to drug discovery or clinically useful prediction of risk will be challenging, and; - 3) through accumulation across studies, GWAS data may become valuable for improved genetic risk differentiation in research in biological psychiatry (e.g., prospective high-risk or resilience studies).

Association between ORMDL3, IL1RL1 and a deletion on chromosome 17q21 with asthma risk in Australia

Genome-wide association studies followed by replication provide a powerful approach to map genetic risk factors for asthma. We sought to search for new variants associated with asthma and attempt to replicate the association with four loci reported previously (ORMDL3, PDE4D, DENND1B and IL1RL1). Genome-wide association analyses of individual single nucleotide polymorphisms (SNPs), rare copy number variants (CNVs) and overall CNV burden were carried out in 986 asthma cases and 1846 asthma-free controls from Australia. The most-associated locus in the SNP analysis was ORMDL3 (rs6503525, P = 4.8 x 10(-)(7)). Five other loci were associated with P < 10(-)(5), most notably the chemokine CXC motif ligand 14 (CXCL14) gene (rs31263, P = 7.8 x 10(-)(6)). We found no evidence for association with the specific risk variants reported recently for PDE4D, DENND1B and ILR1L1. However, a variant in IL1RL1 that is in low linkage disequilibrium with that reported previously was associated with asthma risk after accounting for all variants tested (rs10197862, gene wide P = 0.01). This association replicated convincingly in an independent cohort (P = 2.4 x 10(-)(4)). A 300-kb deletion on chromosome 17q21 was associated with asthma risk, but this did not reach experiment-wide significance. Asthma cases and controls had comparable CNV rates, length and number of genes affected by deletions or duplications. In conclusion, we confirm the association between asthma risk and variants in ORMDL3 and identify a novel risk variant in IL1RL1. Follow-up of the 17q21 deletion in larger cohorts is warranted.

GWAS of butyrylcholinesterase activity identifies four novel loci, independent effects within BCHE and secondary associations with metabolic risk factors

Serum butyrylcholinesterase (BCHE) activity is associated with obesity, blood pressure and biomarkers of cardiovascular and diabetes risk. We have conducted a genome-wide association scan to discover genetic variants affecting BCHE activity, and to clarify whether the associations between BCHE activity and cardiometabolic risk factors are caused by variation in BCHE or whether BCHE variation is secondary to the metabolic abnormalities. We measured serum BCHE in adolescents and adults from three cohorts of Australian twin and family studies. The genotypes from approximately 2.4 million single-nucleotide polymorphisms (SNPs) were available in 8791 participants with BCHE measurements. We detected significant associations with BCHE activity at three independent groups of SNPs at the BCHE locus (P = 5.8 x 10(-262), 7.8 x 10(-47), 2.9 x 10(-12)) and at four other loci: RNPEP (P = 9.4 x 10(-16)), RAPH1-ABI2 (P = 4.1 x 10(-18)), UGT1A1 (P = 4.0 x 10(-8)) and an intergenic region on chromosome 8 (P = 1.4 x 10(-8)). These loci affecting BCHE activity were not associated with metabolic risk factors. On the other hand, SNPs in genes previously associated with metabolic risk had effects on BCHE activity more often than can be explained by chance. In particular, SNPs within FTO and GCKR were associated with BCHE activity, but their effects were partly mediated by body mass index and triglycerides, respectively. We conclude that variation in BCHE activity is due to multiple variants across the spectrum from uncommon/large effect to common/small effect, and partly results from (rather than causes) metabolic abnormalities.

Common SNPs explain a large proportion of the heritability for human height

SNPs discovered by genome-wide association studies (GWASs) account for only a small fraction of the genetic variation of complex traits in human populations. Where is the remaining heritability? We estimated the proportion of variance for human height explained by 294,831 SNPs genotyped on 3,925 unrelated individuals using a linear model analysis, and validated the estimation method with simulations based on the observed genotype data. We show that 45% of variance can be explained by considering all SNPs simultaneously. Thus, most of the heritability is not missing but has not previously been detected because the individual effects are too small to pass stringent significance tests. We provide evidence that the remaining heritability is due to incomplete linkage disequilibrium between causal variants and genotyped SNPs, exacerbated by causal variants having lower minor allele frequency than the SNPs explored to date.

A genome-wide association study of Cloninger's temperament scales: implications for the evolutionary genetics of personality

Variation in personality traits is 30-60% attributed to genetic influences. Attempts to unravel these genetic influences at the molecular level have, so far, been inconclusive. We performed the first genome-wide association study of Cloninger's temperament scales in a sample of 5117 individuals, in order to identify common genetic variants underlying variation in personality. Participants' scores on Harm Avoidance, Novelty Seeking, Reward Dependence, and Persistence were tested for association with 1,252,387 genetic markers. We also performed gene-based association tests and biological pathway analyses. No genetic variants that significantly contribute to personality variation were identified, while our sample provides over 90% power to detect variants that explain only 1% of the trait variance. This indicates that individual common genetic variants of this size or greater do not contribute to personality trait variation, which has important implications regarding the genetic architecture of personality and the evolutionary mechanisms by which heritable variation is maintained.