Search

Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention.

Endometriosis is a leading cause of pain and infertility affecting millions of women globally. Herein, we characterize variation in DNA methylation (DNAm) and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genotype data and methylation in endometrial samples from 984 deeply-phenotyped participants. We estimate that 15.4% of the variation in endometriosis is captured by DNAm and identify significant differences in DNAm profiles associated with stage III/IV endometriosis, endometriosis sub-phenotypes and menstrual cycle phase, including opening of the window for embryo implantation. Menstrual cycle phase was a major source of DNAm variation suggesting cellular and hormonally-driven changes across the cycle can regulate genes and pathways responsible for endometrial physiology and function. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independent cis-mQTLs including 51 associated with risk of endometriosis, highlighting candidate genes contributing to disease risk. Our work provides functional evidence for epigenetic targets contributing to endometriosis risk and pathogenesis. Data generated serve as a valuable resource for understanding tissue-specific effects of methylation on endometrial biology in health and disease.

Endometriosis is associated with increased risk of epithelial ovarian cancers (EOCs). Using data from large endometriosis and EOC genome-wide association meta-analyses, we estimate the genetic correlation and evaluate the causal relationship between genetic liability to endometriosis and EOC histotypes, and identify shared susceptibility loci. We estimate a significant genetic correlation (rg) between endometriosis and clear cell (rg = 0.71), endometrioid (rg = 0.48), and high-grade serous (rg = 0.19) ovarian cancer, associations supported by Mendelian randomization analyses. Bivariate meta-analysis identified 28 loci associated with both endometriosis and EOC, including 19 with evidence for a shared underlying association signal. Differences in the shared risk suggest different underlying pathways may contribute to the relationship between endometriosis and the different histotypes. Functional annotation using transcriptomic and epigenomic profiles of relevant tissues/cells highlights several target genes. This comprehensive analysis reveals profound genetic overlap between endometriosis and EOC histotypes with valuable genomic targets for understanding the biological mechanisms linking the diseases.

STUDY QUESTION Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER We identified four new loci, GNRH1 , FSHR , ZFPM1 , and IPO8 , in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR , have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r  = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S) Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE—European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Skłodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER N/A. [ABSTRACT FROM AUTHOR]

Introduction: Sex steroid hormone fluctuations may underlie both reproductive disorders and sex differences in lifetime depression prevalence. Previous studies report high comorbidity among reproductive disorders and between reproductive disorders and depression. This study sought to assess the multivariate genetic architecture of reproductive disorders and their loading onto a common genetic factor and investigated whether this latent factor shares a common genetic architecture with female depression, including perinatal depression (PND). Method: Using UK Biobank and FinnGen data, genome-wide association meta-analyses were conducted for nine reproductive disorders, and genetic correlation between disorders was estimated. Genomic Structural Equation Modelling identified a latent genetic factor underlying disorders, accounting for their significant genetic correlations. SNPs significantly associated with both latent factor and depression were identified. Results: Excellent model fit existed between a latent factor underlying five reproductive disorders (χ2 (5) = 6.4; AIC = 26.4; CFI = 1.00; SRMR = 0.03) with high standardised loadings for menorrhagia (0.96, SE = 0.05); ovarian cysts (0.94, SE = 0.05); endometriosis (0.83, SE = 0.05); menopausal symptoms (0.77, SE = 0.10); and uterine fibroids (0.65, SE = 0.05). This latent factor was genetically correlated with PND (rG = 0.37, SE = 0.15, p = 1.4e−03), depression in females only (rG = 0.48, SE = 0.06, p = 7.2e−11), and depression in both males and females (MD) (rG = 0.35, SE = 0.03, p = 1.8e−30), with its top locus associated with FSHB/ARL14EP (rs11031006; p = 9.1e−33). SNPs intronic to ESR1, significantly associated with the latent factor, were also associated with PND, female depression, and MD. Conclusion: A common genetic factor, correlated with depression, underlies risk of reproductive disorders, with implications for aetiology and treatment. Genetic variation in ESR1 is associated with reproductive disorders and depression, highlighting the importance of oestrogen signalling for both reproductive and mental health. [ABSTRACT FROM AUTHOR]

Background and Aims Ulcerative colitis [UC] is a major form of inflammatory bowel disease globally. Phenotypic heterogeneity is defined by several variables including age of onset and disease extent. The genetics of disease severity remains poorly understood. To further investigate this, we performed a genome wide association [GWA] study using an extremes of phenotype strategy. Methods We conducted GWA analyses in 311 patients with medically refractory UC [MRUC], 287 with non-medically refractory UC [non-MRUC] and 583 controls. Odds ratios [ORs] were calculated for known risk variants comparing MRUC and non-MRUC, and controls. Results MRUC–control analysis had the greatest yield of genome-wide significant single nucleotide polymorphisms [SNPs] [2018], including lead SNP = rs111838972 [OR = 1.82, p = 6.28 × 10−9] near MMEL1 and a locus in the human leukocyte antigen [HLA] region [lead SNP = rs144717024, OR = 12.23, p = 1.7 × 10−19]. ORs for the lead SNPs were significantly higher in MRUC compared to non-MRUC [ p < 9.0 × 10−6]. No SNPs reached significance in the non-MRUC–control analysis (top SNP, rs7680780 [OR 2.70, p = 5.56 × 10−8). We replicate findings for rs4151651 in the Complement Factor B [CFB] gene and demonstrate significant changes in CFB gene expression in active UC. Detailed HLA analyses support the strong associations with MHC II genes, particularly HLA-DQA1 , HLA-DQB1 and HLA-DRB1 in MRUC. Conclusions Our MRUC subgroup replicates multiple known UC risk variants in contrast to non-MRUC and demonstrates significant differences in effect sizes compared to those published. Non-MRUC cases demonstrate lower ORs similar to those published. Additional risk and prognostic loci may be identified by targeted recruitment of individuals with severe disease. [ABSTRACT FROM AUTHOR]