Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity

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.

Polygenic risk score phenome-wide association study reveals an association between endometriosis and testosterone

BACKGROUND

Endometriosis affects 1 in 9 women, yet it is poorly understood with long diagnostic delays, invasive diagnoses, and poor treatment outcomes. Characterised by the presence of endometrial-like tissue outside of the uterus, its main symptoms are pain and infertility. Endometriosis often co-occurs with other conditions, which may provide insights into the origins of endometriosis.

METHODS

Here a polygenic risk score phenome-wide association study of endometriosis was conducted in the UK Biobank to investigate the pleiotropic effects of a genetic liability to endometriosis. The relationship between the polygenic risk score for endometriosis and health conditions, blood and urine biomarkers and reproductive factors were investigated separately in females, males and females without an endometriosis diagnosis. The relationship between endometriosis and the blood and urine biomarkers was further investigated using genetic correlation and Mendelian randomisation approaches to identify causal relationships.

RESULTS

Multiple health conditions, blood and urine biomarkers and reproductive factors were associated with genetic liability to endometriosis in each group, indicating many endometriosis comorbidities are not dependent on the physical manifestation of endometriosis. Differences in the associated traits between males and females highlighted the importance of sex-specific pathways in the overlap of endometriosis with many other traits. Notably, an association of genetic liability to endometriosis with lower testosterone levels was identified. Follow-up analysis utilising Mendelian randomisation approaches suggested lower testosterone may be causal for both endometriosis and clear cell ovarian cancer.

CONCLUSIONS

This study highlights the diversity of the pleiotropic effects of genetic risk to endometriosis irrespective of a diagnosis of endometriosis. A key finding was the identification of a causal effect of the genetic liability to lower testosterone on endometriosis using Mendelian randomisation.

Gene expression profiles separate endometriosis lesion subtypes and indicate a sensitivity of endometrioma to estrogen suppressive treatments through elevated ESR2 expression

BACKGROUND

Endometriosis is a common, gynaecological disease characterised by the presence of endometrial-like cells growing outside the uterus. Lesions appear at multiple locations, present with variation in appearance, size and depth of invasion. Despite hormones being the recommended first-line treatment, their efficacy, success and side effects vary widely amongst study populations. Current, hormonal medication for endometriosis is designed to suppress systemic oestrogen. Whether these hormones can influence the lesions themselves is not yet clear. Evidence of hormone receptor expression in endometriotic lesions and their ability to respond is conflicting. A variation in their expression, activation of transcriptional co-regulators and the potential to respond may contribute to their variation in patient outcomes. Identifying patients who would benefit from hormonal treatments remain an important goal in endometriosis research.

METHODS

Using gene expression data from endometriosis lesions including endometrioma (OMA, n = 28), superficial peritoneal lesions (SUP, n = 72) and deeply infiltrating lesions (DIE, n = 78), we performed principal component analysis, differential gene expression and gene correlation analyses to assess the impact of menstrual stage, lesion subtype and hormonal treatment on the gene expression.

RESULTS

The gene expression profiles did not vary based on menstrual stage, but could distinguish lesion subtypes with OMA significantly differentiating from both SUP and DIE. Additionally, the effect of oestrogen suppression medication altered the gene expression profile in OMA, while such effect was not observed in SUP or DIE. Analysis of the target receptors for hormonal medication indicated ESR2 was differentially expressed in OMA and that genes that correlated with ESR2 varied significantly between medicated and non-medicated OMA samples.

CONCLUSIONS

Our results demonstrate of the different lesion types OMA present with strongest response to hormonal treatment directly through ESR2. The data suggests that there may be the potential to target treatment options to individual patients based on pre-surgical diagnoses.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

Evidence of shared genetic factors in the etiology of gastrointestinal disorders and endometriosis and clinical implications for disease management

In clinical practice, the co-existence of endometriosis and gastrointestinal symptoms is often observed. Using large-scale datasets, we report a genetic correlation between endometriosis and irritable bowel syndrome (IBS), peptic ulcer disease (PUD), gastro-esophageal reflux disease (GORD), and a combined GORD/PUD medicated (GPM) phenotype. Mendelian randomization analyses support a causal relationship between genetic predisposition to endometriosis and IBS and GPM. Identification of shared risk loci highlights biological pathways that may contribute to the pathogenesis of both diseases, including estrogen regulation and inflammation, and potential therapeutic drug targets (CCKBR; PDE4B). Higher use of IBS, GORD, and PUD medications in women with endometriosis and higher use of hormone therapies in women with IBS, GORD, and PUD, support the co-occurrence of these conditions and highlight the potential for drug repositioning and drug contraindications. Our results provide evidence of shared disease etiology and have important clinical implications for diagnostic and treatment decisions for both diseases.

A molecular staging model for accurately dating the endometrial biopsy

Natural variability in menstrual cycle length, coupled with rapid changes in endometrial gene expression, makes it difficult to accurately define and compare different stages of the endometrial cycle. Here we develop and validate a method for precisely determining endometrial cycle stage based on global gene expression. Our 'molecular staging model' reveals significant and remarkably synchronised daily changes in expression for over 3400 endometrial genes throughout the cycle, with the most dramatic changes occurring during the secretory phase. Our study significantly extends existing data on the endometrial transcriptome, and for the first time enables identification of differentially expressed endometrial genes with increasing age and different ethnicities. It also allows reinterpretation of all endometrial RNA-seq and array data that has been published to date. Our molecular staging model will significantly advance understanding of endometrial-related disorders that affect nearly all women at some stage of their lives, such as heavy menstrual bleeding, endometriosis, adenomyosis, and recurrent implantation failure.

Insights from Mendelian randomization and genetic correlation analyses into the relationship between endometriosis and its comorbidities

BACKGROUND

Endometriosis remains a poorly understood disease, despite its high prevalence and debilitating symptoms. The overlap in symptoms and the increased risk of multiple other traits in women with endometriosis is becoming increasingly apparent through epidemiological data. Genetic studies offer a method of investigating these comorbid relationships through the assessment of causal relationships with Mendelian randomization (MR), as well as identification of shared genetic variants and genes involved across traits. This has the capacity to identify risk factors for endometriosis as well as provide insight into the aetiology of disease.

OBJECTIVE AND RATIONALE

We aim to review the current literature assessing the relationship between endometriosis and other traits using genomic data, primarily through the methods of MR and genetic correlation. We critically examine the limitations of these studies in accordance with the assumptions of the utilized methods.

SEARCH METHODS

The PubMed database was used to search for peer-reviewed original research articles using the terms 'Mendelian randomization endometriosis' and '"genetic correlation" endometriosis'. Additionally, a Google Scholar search using the terms '"endometriosis" "mendelian randomization" "genetic correlation"' was performed. All relevant publications (n = 21) published up until 7 October 2022 were included in this review. Upon compilation of all traits with published MR and/or genetic correlation with endometriosis, additional epidemiological and genetic information on their comorbidity with endometriosis was sourced by searching for the trait in conjunction with 'endometriosis' on Google Scholar.

OUTCOMES

The association between endometriosis and multiple pain, gynaecological, cancer, inflammatory, gastrointestinal, psychological, and anthropometric traits has been assessed using MR analysis and genetic correlation analysis. Genetic correlation analyses provide evidence that genetic factors contributing to endometriosis are shared with multiple traits: migraine, uterine fibroids, subtypes of ovarian cancer, melanoma, asthma, gastro-oesophageal reflux disease, gastritis/duodenitis, and depression, suggesting the involvement of multiple biological mechanisms in endometriosis. The assessment of causality with MR has revealed several potential causes (e.g. depression) and outcomes (e.g. ovarian cancer and uterine fibroids) of a genetic predisposition to endometriosis; however, interpretation of these results requires consideration of potential violations of the MR assumptions.

WIDER IMPLICATIONS

Genomic studies have demonstrated that there is a molecular basis for the co-occurrence of endometriosis with other traits. Dissection of this overlap has identified shared genes and pathways, which provide insight into the biology of endometriosis. Thoughtful MR studies are necessary to ascertain causality of the comorbidities of endometriosis. Given the significant diagnostic delay of endometriosis of 7-11 years, determining risk factors is necessary to aid diagnosis and reduce the disease burden. Identification of traits for which endometriosis is a risk factor is important for holistic treatment and counselling of the patient. The use of genomic data to disentangle the overlap of endometriosis with other traits has provided insights into the aetiology of endometriosis.

Genomic characterisation of the overlap of endometriosis with 76 comorbidities identifies pleiotropic and causal mechanisms underlying disease risk

Comorbid conditions can be driven by underlying pleiotropic and causal mechanisms that can provide insights into shared molecular and biological processes contributing to disease risk. Endometriosis is a chronic condition affecting one in nine women of reproductive age and poses many challenges including lengthy diagnostic delays and limited treatment efficacy owing to poor understanding of disease aetiology. To shed light on the underlying biological mechanisms and to identify potential risk factors, we examine the epidemiological and genomic relationship between endometriosis and its comorbidities. In the UK Biobank 292 ICD10 codes were epidemiologically correlated with endometriosis diagnosis, including gynaecological, immune, infection, pain, psychiatric, cancer, gastrointestinal, urinary, bone and cardiovascular traits. A subset of the identified comorbidities (n = 76) underwent follow-up genetic analysis. Whilst Mendelian randomisation suggested causality was not responsible for most comorbid relationships, 22 traits were genetically correlated with endometriosis, including pain, gynaecological and gastrointestinal traits, suggestive of a shared genetic background. Pleiotropic genetic variants and genes were identified using gene-based and colocalisation analysis. Shared genetic risk factors and potential target genes suggest a diverse collection of biological systems are involved in these comorbid relationships including coagulation factors, development of the female reproductive tract and cell proliferation. These findings highlight the diversity of traits with epidemiological and genomic overlap with endometriosis and implicate a key role for pleiotropy in the comorbid relationships.

Global endometrial DNA methylation analysis reveals insights into mQTL regulation and associated endometriosis disease risk and endometrial function

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.

Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis

Identifying tissue-specific molecular signatures of active regulatory elements is critical to understanding gene regulatory mechanisms. In this study, transcription start sites (TSS) and enhancers were identified using Cap analysis of gene expression (CAGE) across endometrial stromal cell (ESC) samples obtained from women with (n = 4) and without endometriosis (n = 4). ESC TSSs and enhancers were compared to those reported in other tissue and cell types in FANTOM5 and were integrated with RNA-seq and ATAC-seq data from the same samples for regulatory activity and network analyses. CAGE tag count differences between women with and without endometriosis were statistically tested and tags within close proximity to genetic variants associated with endometriosis risk were identified. Over 90% of tag clusters mapping to promoters were observed in cells and tissues in FANTOM5. However, some potential cell-type-specific promoters and enhancers were also observed. Regions of open chromatin identified using ATAC-seq provided further evidence of the active transcriptional regions identified by CAGE. Despite the small sample number, there was evidence of differences associated with endometriosis at 210 consensus clusters, including IGFBP5, CALD1 and OXTR. ESC TSSs were also located within loci associated with endometriosis risk from genome-wide association studies. This study provides novel evidence of transcriptional differences in endometrial stromal cells associated with endometriosis and provides a valuable cell-type specific resource of active TSSs and enhancers in endometrial stromal cells.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Rare genetic variants underlie outlying levels of DNA methylation and gene-expression

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on DNA methylation (DNAm) as effect sizes are expected to be larger for molecular traits compared with complex traits. Here, we investigate DNAm in healthy ageing populations-the Lothian Birth Cohorts of 1921 and 1936-and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic single nucleotide polymorphisms (SNPs) within 1 kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels toward the tails of the population distribution. This study demonstrates the role of rare SNPs in the phenotypic variation of DNAm and the effect of extreme levels of DNAm on gene expression.

Establishing the Australian National Endometriosis Clinical and Scientific Trials (NECST) Registry: A protocol paper

Endometriosis is a common yet under-recognised chronic inflammatory disease, affecting 176 million women, trans and gender diverse people globally. The National Endometriosis Clinical and Scientific Trials (NECST) Registry is a new clinical registry, collecting and tracking diagnostic and treatment data, and patient-reported outcomes on people with endometriosis. The registry is a research priority action item from the 2018 National Action Plan for Endometriosis and aims to provide, large-scale, national and longitudinal population-based data on endometriosis. Working groups (consisting of patients with endometriosis, clinicians and researchers) developing the NECST Registry data dictionary and data collection platform started in 2019. Our data dictionary was developed based on existing and validated questionnaires, tools, meta-data and data cubes - World Endometriosis Research Foundation (WERF) Endometriosis Phenome and Biobanking Harmonisation Project (EPHect), endometriosis CORE outcomes set, patient-reported outcome measures, the International Statistical Classification of Diseases-10th Revision Australian Modification diagnosis codes, and Australian Government datasets: Australian Institute for Health and Welfare (for sociodemographic data), Medicare Benefits Schedule (MBS; for medical procedures) and the Pharmaceutical Benefits Scheme (PBS; for medical therapies). The resulting NECST Registry is an online, secure cloud-based database; prospectively collecting minimum core clinical and health data across eight patient and clinician modules and longitudinal data tracking disease life course. The NECST Registry has ethics approval (HREC/62508/MonH-2020) and is registered on the Australian New Zealand Clinical Trials Registry (ACTRN12622000987763).

New concepts on the etiology of endometriosis

Endometriosis is a serious, chronic disorder where endometrial tissue grows outside the uterus, causing severe pelvic pain and infertility. It affects 11% of women. Endometriosis is a multifactorial disorder of unclear etiology, although retrograde menstruation plays a major role. It has a genetic component with over 40 genetic risk factors mapped, although their mechanism of action is still emerging. New evidence suggests a role for retrograde menstruation of endometrial stem/progenitor cells, now that identifying markers of these cells are available. Recent lineage tracing and tissue clearing microscopy and 3D reconstruction has provided new understanding of endometrial glandular structure, particularly the horizontal orientation and interconnection of basalis glands. New sequencing technologies, particularly whole genome DNA sequencing are revealing somatic mutations, including in cancer driver genes, in normal and eutopic endometrium of patients with endometriosis, as well as ectopic endometriotic lesions. Methylome sequencing is offering insight into the regulation of genes and the role of the environmental factors. Single cell RNA sequencing reveals the transcriptome of individual endometrial cells, shedding new light on the diversity and range of cellular subpopulations of the major cell types present in the endometrium and in endometriotic lesions. New endometrial epithelial organoid cultures replicating glandular epithelium are providing tractable models for studying endometriosis. Organoids derived from menstrual fluid offer a non-invasive source of endometrial tissue and a new avenue for testing drugs and developing personalized medicine for treating endometriosis. These new approaches are rapidly advancing our understanding of endometriosis etiology.

The genetic basis of endometriosis and comorbidity with other pain and inflammatory conditions

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.

Multi-trait genome-wide association study of opioid addiction: OPRM1 and beyond

Opioid addiction (OA) is moderately heritable, yet only rs1799971, the A118G variant in OPRM1, has been identified as a genome-wide significant association with OA and independently replicated. We applied genomic structural equation modeling to conduct a GWAS of the new Genetics of Opioid Addiction Consortium (GENOA) data together with published studies (Psychiatric Genomics Consortium, Million Veteran Program, and Partners Health), comprising 23,367 cases and effective sample size of 88,114 individuals of European ancestry. Genetic correlations among the various OA phenotypes were uniformly high (rg > 0.9). We observed the strongest evidence to date for OPRM1: lead SNP rs9478500 (p = 2.56 × 10-9). Gene-based analyses identified novel genome-wide significant associations with PPP6C and FURIN. Variants within these loci appear to be pleiotropic for addiction and related traits.

Clinical characteristics and surgical management of endometriosis-associated infertility: A multicenter prospective cohort study

OBJECTIVE

To study clinical, surgical characteristics and the relationship between endometriosis lesion types and conception rate after surgery in infertile women with endometriosis.

METHODS

A prospective, multicenter cohort of 204 women (age 20-35 years) with endometriosis was followed up post-surgery between November 2017 and February 2020 at three tertiary-care hospitals.

RESULTS

Based on the severity of endometriosis lesion type, deep infiltrating endometriosis (DIE) (81/204, 39.7%) was the most common lesion; followed by ovarian endometriosis (OMA) (64/204, 31.4%), and superficial peritoneal endometriosis (SUP) (59/204, 28.9%). Endometriosis patients had a single lesion type (94/204, 46.1%), two lesion types (77/204, 37.7%), or three lesion types (33/204, 16.2%) with significant differences between regions (P < 0.001). Around 40% (37/95) of obese women had SUP (P = 0.003) whereas 78% (14/18) of underweight women had DIE (P < 0.001). Significant differences in mean Endometriosis Fertility Index scores between endometriosis lesion types and patients with one, two, and three types of lesions were observed (P < 0.001). The majority (22/32, 68.8%) of the women conceived naturally after the surgery. Half (16/32; 50%) of the women with a single lesion type conceived after the surgery; of which most (13/16, 81.2%) had SUP, followed by OMA (2/16, 12.5%), and DIE (1/16, 6.3%).

CONCLUSION

Women with SUP and only one type of endometriotic lesion were more likely to conceive post-surgery.

Changing genetic architecture of body mass index from infancy to early adulthood: an individual based pooled analysis of 25 twin cohorts

Background

Body mass index (BMI) shows strong continuity over childhood and adolescence and high childhood BMI is the strongest predictor of adult obesity. Genetic factors strongly contribute to this continuity, but it is still poorly known how their contribution changes over childhood and adolescence. Thus, we used the genetic twin design to estimate the genetic correlations of BMI from infancy to adulthood and compared them to the genetic correlations of height.

Methods

We pooled individual level data from 25 longitudinal twin cohorts including 38,530 complete twin pairs and having 283,766 longitudinal height and weight measures. The data were analyzed using Cholesky decomposition offering genetic and environmental correlations of BMI and height between all age combinations from 1 to 19 years of age.

Results

The genetic correlations of BMI and height were stronger than the trait correlations. For BMI, we found that genetic correlations decreased as the age between the assessments increased, a trend that was especially visible from early to middle childhood. In contrast, for height, the genetic correlations were strong between all ages. Age-to-age correlations between environmental factors shared by co-twins were found for BMI in early childhood but disappeared altogether by middle childhood. For height, shared environmental correlations persisted from infancy to adulthood.

Conclusions

Our results suggest that the genes affecting BMI change over childhood and adolescence leading to decreasing age-to-age genetic correlations. This change is especially visible from early to middle childhood indicating that new genetic factors start to affect BMI in middle childhood. Identifying mediating pathways of these genetic factors can open possibilities for interventions, especially for those children with high genetic predisposition to adult obesity.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

Identifying the Common Genetic Basis of Antidepressant Response

Background

Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction.

Methods

Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA.

Results

Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response.

Conclusions

This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

A multi-level investigation of the genetic relationship between endometriosis and ovarian cancer histotypes

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 (r_g) between endometriosis and clear cell (r_g = 0.71), endometrioid (r_g = 0.48), and high-grade serous (r_g = 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.

•

Endometriosis is genetically correlated with CCOC, ENOC, and HGSOC

•

Genetic liability to endometriosis confers risk of these EOC histotypes

•

Profound colocalization of genetic associations at endometriosis and EOC risk loci

•

Functional annotation highlights shared target genes elucidating the genetic link

Genetic Regulation of Transcription in the Endometrium in Health and Disease

The endometrium is a complex and dynamic tissue essential for fertility and implicated in many reproductive disorders. The tissue consists of glandular epithelium and vascularised stroma and is unique because it is constantly shed and regrown with each menstrual cycle, generating up to 10 mm of new mucosa. Consequently, there are marked changes in cell composition and gene expression across the menstrual cycle. Recent evidence shows expression of many genes is influenced by genetic variation between individuals. We and others have reported evidence for genetic effects on hundreds of genes in endometrium. The genetic factors influencing endometrial gene expression are highly correlated with the genetic effects on expression in other reproductive (e.g., in uterus and ovary) and digestive tissues (e.g., salivary gland and stomach), supporting a shared genetic regulation of gene expression in biologically similar tissues. There is also increasing evidence for cell specific genetic effects for some genes. Sample size for studies in endometrium are modest and results from the larger studies of gene expression in blood report genetic effects for a much higher proportion of genes than currently reported for endometrium. There is also emerging evidence for the importance of genetic variation on RNA splicing. Gene mapping studies for common disease, including diseases associated with endometrium, show most variation maps to intergenic regulatory regions. It is likely that genetic risk factors for disease function through modifying the program of cell specific gene expression. The emerging evidence from our gene mapping studies coupled with tissue specific studies, and the GTEx, eQTLGen and EpiMap projects, show we need to expand our understanding of the complex regulation of gene expression. These data also help to link disease genetic risk factors to specific target genes. Combining our data on genetic regulation of gene expression in endometrium, and cell types within the endometrium with gene mapping data for endometriosis and related diseases is beginning to uncover the specific genes and pathways responsible for increased risk of these diseases.

Educational attainment of same-sex and opposite-sex dizygotic twins: An individual-level pooled study of 19 twin cohorts

Comparing twins from same- and opposite-sex pairs can provide information on potential sex differences in a variety of outcomes, including socioeconomic-related outcomes such as educational attainment. It has been suggested that this design can be applied to examine the putative role of intrauterine exposure to testosterone for educational attainment, but the evidence is still disputed. Thus, we established an international database of twin data from 11 countries with 88,290 individual dizygotic twins born over 100 years and tested for differences between twins from same- and opposite-sex dizygotic pairs in educational attainment. Effect sizes with 95% confidence intervals (CI) were estimated by linear regression models after adjusting for birth year and twin study cohort. In contrast to the hypothesis, no difference was found in women (β = -0.05 educational years, 95% CI -0.11, 0.02). However, men with a same-sex co-twin were slightly more educated than men having an opposite-sex co-twin (β = 0.14 educational years, 95% CI 0.07, 0.21). No consistent differences in effect sizes were found between individual twin study cohorts representing Europe, the USA, and Australia or over the cohorts born during the 20th century, during which period the sex differences in education reversed favoring women in the latest birth cohorts. Further, no interaction was found with maternal or paternal education. Our results contradict the hypothesis that there would be differences in the intrauterine testosterone levels between same-sex and opposite-sex female twins affecting education. Our findings in men may point to social dynamics within same-sex twin pairs that may benefit men in their educational careers.

Longitudinal changes in employment following a diagnosis of endometriosis: Findings from an Australian cohort study

PURPOSE

Endometriosis is a chronic inflammatory disease affecting the reproductive, gastrointestinal and urinary systems. We examined changes in labour force participation among women with endometriosis following diagnosis.

METHODS

We analysed data from 4494 women born in 1973-78 from the Australian Longitudinal Study on Women's Health. We used multinomial logistic regression models with generalised estimating equations to examine changes in labour force participation among 468 women with surgically confirmed endometriosis, and 375 women with clinically suspected endometriosis, relative to a comparison group of 4151 women without endometriosis.

RESULTS

At diagnosis, women with surgically confirmed endometriosis were somewhat more likely to be working part-time (OR 1.26, 95% CI 0.94-1.68) or unemployed (OR 1.46, 95% CI 0.96-2.23) than before diagnosis. After diagnosis, women with surgically confirmed endometriosis remained somewhat more likely to be working part-time (OR 1.26, 95% CI 0.88-1.80) but were significantly more likely to be unemployed (OR 1.85, 95% CI 1.16-2.96) than before diagnosis. Labour force participation for women with clinically suspected endometriosis did not differ from women without endometriosis at diagnosis and did not change over time.

CONCLUSIONS

Women with surgically confirmed endometriosis transitioned out the labour force following diagnosis. Supportive workplace practices may help women remain in the labour force.

Identical twins carry a persistent epigenetic signature of early genome programming

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.

Neuropeptide S receptor 1 is a nonhormonal treatment target in endometriosis

Endometriosis is a common chronic inflammatory condition causing pelvic pain and infertility in women, with limited treatment options and 50% heritability. We leveraged genetic analyses in two species with spontaneous endometriosis, humans and the rhesus macaque, to uncover treatment targets. We sequenced DNA from 32 human families contributing to a genetic linkage signal on chromosome 7p13-15 and observed significant overrepresentation of predicted deleterious low-frequency coding variants in NPSR1, the gene encoding neuropeptide S receptor 1, in cases (predominantly stage III/IV) versus controls (P = 7.8 × 10^-4). Significant linkage to the region orthologous to human 7p13-15 was replicated in a pedigree of 849 rhesus macaques (P = 0.0095). Targeted association analyses in 3194 surgically confirmed, unrelated cases and 7060 controls revealed that a common insertion/deletion variant, rs142885915, was significantly associated with stage III/IV endometriosis (P = 5.2 × 10^-5; odds ratio, 1.23; 95% CI, 1.09 to 1.39). Immunohistochemistry, qRT-PCR, and flow cytometry experiments demonstrated that NPSR1 was expressed in glandular epithelium from eutopic and ectopic endometrium, and on monocytes in peritoneal fluid. The NPSR1 inhibitor SHA 68R blocked NPSR1-mediated signaling, proinflammatory TNF-α release, and monocyte chemotaxis in vitro (P < 0.01), and led to a significant reduction of inflammatory cell infiltrate and abdominal pain (P < 0.05) in a mouse model of peritoneal inflammation as well as in a mouse model of endometriosis. We conclude that the NPSR1/NPS system is a genetically validated, nonhormonal target for the treatment of endometriosis with likely increased relevance to stage III/IV disease.

Genetic Relationship Between Endometriosis and Melanoma

Epidemiological studies have observed that risk of endometriosis is associated with history of cutaneous melanoma and vice versa. Evidence for shared biological mechanisms between the two traits is limited. The aim of this study was to investigate the genetic correlation and causal relationship between endometriosis and melanoma. Summary statistics from genome-wide association meta-analyses (GWAS) for endometriosis and melanoma were used to estimate the genetic correlation between the traits and Mendelian randomization was used to test for a causal association. When using summary statistics from separate female and male melanoma cohorts we identified a significant positive genetic correlation between melanoma in females and endometriosis (rg = 0.144, se = 0.065, p = 0.025). However, we find no evidence of a correlation between endometriosis and melanoma in males or a combined melanoma dataset. Endometriosis was not genetically correlated with skin color, red hair, childhood sunburn occasions, ease of skin tanning, or nevus count suggesting that the correlation between endometriosis and melanoma in females is unlikely to be influenced by pigmentary traits. Mendelian Randomization analyses also provided evidence for a relationship between the genetic risk of melanoma in females and endometriosis. Colocalization analysis identified 27 genomic loci jointly associated with the two diseases regions that contain different causal variants influencing each trait independently. This study provides evidence of a small genetic correlation and relationship between the genetic risk of melanoma in females and endometriosis. Genetic risk does not equate to disease occurrence and differences in the pathogenesis and age of onset of both diseases means it is unlikely that occurrence of melanoma causes endometriosis. This study instead provides evidence that having an increased genetic risk for melanoma in females is related to increased risk of endometriosis. Larger GWAS studies with increased power will be required to further investigate these associations.

Genetic insights into biological mechanisms governing human ovarian ageing

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Genetic analyses of gynecological disease identify genetic relationships between uterine fibroids and endometrial cancer, and a novel endometrial cancer genetic risk region at the WNT4 1p36.12 locus

Endometriosis, polycystic ovary syndrome (PCOS) and uterine fibroids have been proposed as endometrial cancer risk factors; however, disentangling their relationships with endometrial cancer is complicated due to shared risk factors and comorbidities. Using genome-wide association study (GWAS) data, we explored the relationships between these non-cancerous gynecological diseases and endometrial cancer risk by assessing genetic correlation, causal relationships and shared risk loci. We found significant genetic correlation between endometrial cancer and PCOS, and uterine fibroids. Adjustment for genetically predicted body mass index (a risk factor for PCOS, uterine fibroids and endometrial cancer) substantially attenuated the genetic correlation between endometrial cancer and PCOS but did not affect the correlation with uterine fibroids. Mendelian randomization analyses suggested a causal relationship between only uterine fibroids and endometrial cancer. Gene-based analyses revealed risk regions shared between endometrial cancer and endometriosis, and uterine fibroids. Multi-trait GWAS analysis of endometrial cancer and the genetically correlated gynecological diseases identified a novel genome-wide significant endometrial cancer risk locus at 1p36.12, which replicated in an independent endometrial cancer dataset. Interrogation of functional genomic data at 1p36.12 revealed biologically relevant genes, including WNT4 which is necessary for the development of the female reproductive system. In summary, our study provides genetic evidence for a causal relationship between uterine fibroids and endometrial cancer. It further provides evidence that the comorbidity of endometrial cancer, PCOS and uterine fibroids may partly be due to shared genetic architecture. Notably, this shared architecture has revealed a novel genome-wide risk locus for endometrial cancer.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Meta-analysis of genome-wide DNA methylation identifies shared associations across neurodegenerative disorders

BACKGROUND

People with neurodegenerative disorders show diverse clinical syndromes, genetic heterogeneity, and distinct brain pathological changes, but studies report overlap between these features. DNA methylation (DNAm) provides a way to explore this overlap and heterogeneity as it is determined by the combined effects of genetic variation and the environment. In this study, we aim to identify shared blood DNAm differences between controls and people with Alzheimer's disease, amyotrophic lateral sclerosis, and Parkinson's disease.

RESULTS

We use a mixed-linear model method (MOMENT) that accounts for the effect of (un)known confounders, to test for the association of each DNAm site with each disorder. While only three probes are found to be genome-wide significant in each MOMENT association analysis of amyotrophic lateral sclerosis and Parkinson's disease (and none with Alzheimer's disease), a fixed-effects meta-analysis of the three disorders results in 12 genome-wide significant differentially methylated positions. Predicted immune cell-type proportions are disrupted across all neurodegenerative disorders. Protein inflammatory markers are correlated with profile sum-scores derived from disease-associated immune cell-type proportions in a healthy aging cohort. In contrast, they are not correlated with MOMENT DNAm-derived profile sum-scores, calculated using effect sizes of the 12 differentially methylated positions as weights.

CONCLUSIONS

We identify shared differentially methylated positions in whole blood between neurodegenerative disorders that point to shared pathogenic mechanisms. These shared differentially methylated positions may reflect causes or consequences of disease, but they are unlikely to reflect cell-type proportion differences.

Genetic Regulation of Physiological Reproductive Lifespan and Female Fertility

There is substantial genetic variation for common traits associated with reproductive lifespan and for common diseases influencing female fertility. Progress in high-throughput sequencing and genome-wide association studies (GWAS) have transformed our understanding of common genetic risk factors for complex traits and diseases influencing reproductive lifespan and fertility. The data emerging from GWAS demonstrate the utility of genetics to explain epidemiological observations, revealing shared biological pathways linking puberty timing, fertility, reproductive ageing and health outcomes. The observations also identify unique genetic risk factors specific to different reproductive diseases impacting on female fertility. Sequencing in patients with primary ovarian insufficiency (POI) have identified mutations in a large number of genes while GWAS have revealed shared genetic risk factors for POI and ovarian ageing. Studies on age at menopause implicate DNA damage/repair genes with implications for follicle health and ageing. In addition to the discovery of individual genes and pathways, the increasingly powerful studies on common genetic risk factors help interpret the underlying relationships and direction of causation in the regulation of reproductive lifespan, fertility and related traits.

Genetic risk factors for endometriosis near estrogen receptor 1 and coexpression of genes in this region in endometrium

The etiology and pathogenesis of endometriosis are complex with both genetic and environmental factors contributing to disease risk. Genome-wide association studies (GWAS) have identified multiple signals in the estrogen receptor 1 (ESR1) region associated with endometriosis and other reproductive traits and diseases. In addition, candidate gene association studies identified signals in the ESR1 region associated with endometriosis risk suggesting genetic regulation of genes in this region may be important for reproductive health. This study aimed to investigate hormonal and genetic regulation of genes in the ESR1 region in human endometrium. Changes in serum oestradiol and progesterone concentrations and expression of hormone receptors ESR1 and progesterone receptor (PGR) were assessed in endometrial samples from 135 women collected at various stages of the menstrual cycle. Correlation between hormone concentrations, receptor expression and expression of genes in the ESR1 locus was investigated. The effect of endometriosis risk variants on expression of genes in the region was analyzed to identify gene targets. Hormone concentrations and receptor expression varied significantly across the menstrual cycle. Expression of genes in the ESR1 region correlated with progesterone concentration; however, they were more strongly correlated with expression of ESR1 and PGR suggesting coregulation of genes. There was no evidence that endometriosis risk variants directly regulated expression of genes in the region. Limited sample size and cellular heterogeneity in endometrial tissue may impact the ability to detect significant genetic effects on gene expression. Effects of these variants should be validated in a larger dataset and in relevant individual cell types.

A large-scale genome-wide association study meta-analysis of cannabis use disorder

BACKGROUND

Variation in liability to cannabis use disorder has a strong genetic component (estimated twin and family heritability about 50-70%) and is associated with negative outcomes, including increased risk of psychopathology. The aim of the study was to conduct a large genome-wide association study (GWAS) to identify novel genetic variants associated with cannabis use disorder.

METHODS

To conduct this GWAS meta-analysis of cannabis use disorder and identify associations with genetic loci, we used samples from the Psychiatric Genomics Consortium Substance Use Disorders working group, iPSYCH, and deCODE (20 916 case samples, 363 116 control samples in total), contrasting cannabis use disorder cases with controls. To examine the genetic overlap between cannabis use disorder and 22 traits of interest (chosen because of previously published phenotypic correlations [eg, psychiatric disorders] or hypothesised associations [eg, chronotype] with cannabis use disorder), we used linkage disequilibrium score regression to calculate genetic correlations.

FINDINGS

We identified two genome-wide significant loci: a novel chromosome 7 locus (FOXP2, lead single-nucleotide polymorphism [SNP] rs7783012; odds ratio [OR] 1·11, 95% CI 1·07-1·15, p=1·84 × 10-9) and the previously identified chromosome 8 locus (near CHRNA2 and EPHX2, lead SNP rs4732724; OR 0·89, 95% CI 0·86-0·93, p=6·46 × 10-9). Cannabis use disorder and cannabis use were genetically correlated (rg 0·50, p=1·50 × 10-21), but they showed significantly different genetic correlations with 12 of the 22 traits we tested, suggesting at least partially different genetic underpinnings of cannabis use and cannabis use disorder. Cannabis use disorder was positively genetically correlated with other psychopathology, including ADHD, major depression, and schizophrenia.

INTERPRETATION

These findings support the theory that cannabis use disorder has shared genetic liability with other psychopathology, and there is a distinction between genetic liability to cannabis use and cannabis use disorder.

FUNDING

National Institute of Mental Health; National Institute on Alcohol Abuse and Alcoholism; National Institute on Drug Abuse; Center for Genomics and Personalized Medicine and the Centre for Integrative Sequencing; The European Commission, Horizon 2020; National Institute of Child Health and Human Development; Health Research Council of New Zealand; National Institute on Aging; Wellcome Trust Case Control Consortium; UK Research and Innovation Medical Research Council (UKRI MRC); The Brain & Behavior Research Foundation; National Institute on Deafness and Other Communication Disorders; Substance Abuse and Mental Health Services Administration (SAMHSA); National Institute of Biomedical Imaging and Bioengineering; National Health and Medical Research Council (NHMRC) Australia; Tobacco-Related Disease Research Program of the University of California; Families for Borderline Personality Disorder Research (Beth and Rob Elliott) 2018 NARSAD Young Investigator Grant; The National Child Health Research Foundation (Cure Kids); The Canterbury Medical Research Foundation; The New Zealand Lottery Grants Board; The University of Otago; The Carney Centre for Pharmacogenomics; The James Hume Bequest Fund; National Institutes of Health: Genes, Environment and Health Initiative; National Institutes of Health; National Cancer Institute; The William T Grant Foundation; Australian Research Council; The Virginia Tobacco Settlement Foundation; The VISN 1 and VISN 4 Mental Illness Research, Education, and Clinical Centers of the US Department of Veterans Affairs; The 5th Framework Programme (FP-5) GenomEUtwin Project; The Lundbeck Foundation; NIH-funded Shared Instrumentation Grant S10RR025141; Clinical Translational Sciences Award grants; National Institute of Neurological Disorders and Stroke; National Heart, Lung, and Blood Institute; National Institute of General Medical Sciences.

Commentary: lessons from molecular genetic studies on reporting false-positive results

Poor replication of published research results is the subject of debate. A common problem is the failure to adequately account for multiple testing issues. In this regard, the evolution of mapping studies to identify genetic risk factors for common diseases has been instructive. Large genome-wide association studies (GWAS) reliably detect the genetic factors with small effects that contribute to risk for many common diseases. GWAS superseded candidate gene studies from the previous decade and looking back, almost no genetic risk factors reported from earlier candidate gene studies replicate in the GWAS results. Candidate gene studies often used small samples and failed to appreciate and adequately account for the multiple testing issues. The failure to replicate results from most candidate gene studies highlights the importance of study power and appropriate statistical analysis to prevent publication of false-positive results.

Cellular Origins of Endometriosis: Towards Novel Diagnostics and Therapeutics

Endometriosis remains an enigmatic disease of unknown etiology, with delayed diagnosis and poor therapeutic options. This review will discuss the cellular, physiological, and genomic evidence of Sampson's hypothesis of retrograde menstruation as a cause of pelvic endometriosis and as the basis of phenotypic heterogeneity of the disease. We postulate that collaborative research at the single cell level focused on unlocking the cellular, physiological, and genomic mechanisms of endometriosis will be accompanied by advances in personalized diagnosis and therapies that target unique subtypes of endometriosis disease. These advances will address the clinical conundrums of endometriosis clinical care—including diagnostic delay, suboptimal treatments, disease recurrence, infertility, chronic pelvic pain, and quality of life. There is an urgent need to improve outcomes for women with endometriosis. To achieve this, it is imperative that we understand which cells form the lesions, how they arrive at distant sites, and what factors govern their ability to survive and invade at ectopic locations. This review proposes new research avenues to address these basic questions of endometriosis pathobiology that will lay the foundations for new diagnostic tools and treatment pathways.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma.

Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness.

To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness.

The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry.

We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10⁻⁸) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Multiplex melanoma families are enriched for polygenic risk

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental etiological factors such as ultraviolet radiation, cutaneous melanoma has a strong genetic component. Genetic risks for cutaneous melanoma range from rare, high-penetrance mutations to common, low-penetrance variants. Known high-penetrance mutations account for only about half of all densely affected cutaneous melanoma families, and the causes of familial clustering in the remainder are unknown. We hypothesize that some clustering is due to the cumulative effect of a large number of variants of individually small effect. Common, low-penetrance genetic risk variants can be combined into polygenic risk scores. We used a polygenic risk score for cutaneous melanoma to compare families without known high-penetrance mutations with unrelated melanoma cases and melanoma-free controls. Family members had significantly higher mean polygenic load for cutaneous melanoma than unrelated cases or melanoma-free healthy controls (Bonferroni-corrected t-test P = 1.5 × 10-5 and 6.3 × 10-45, respectively). Whole genome sequencing of germline DNA from 51 members of 21 families with low polygenic risk for melanoma identified a CDKN2A p.G101W mutation in a single family but no other candidate high-penetrance melanoma susceptibility genes. This work provides further evidence that melanoma, like many other common complex disorders, can arise from the joint action of multiple predisposing factors, including rare high-penetrance mutations, as well as via a combination of large numbers of alleles of small effect.

Genetic analysis of endometriosis and depression identifies shared loci and implicates causal links with gastric mucosa abnormality

Evidence from observational studies indicates that endometriosis and depression often co-occur. However, conflicting evidence exists, and the etiology as well as biological mechanisms underlying their comorbidity remain unknown. Utilizing genome-wide association study (GWAS) data, we comprehensively assessed the relationship between endometriosis and depression. Single nucleotide polymorphism effect concordance analysis (SECA) found a significant genetic overlap between endometriosis and depression (P_{Fsig-permuted} = 9.99 × 10^-4). Linkage disequilibrium score regression (LDSC) analysis estimated a positive and highly significant genetic correlation between the two traits (r_G = 0.27, P = 8.85 × 10^-27). A meta-analysis of endometriosis and depression GWAS (sample size = 709,111), identified 20 independent genome-wide significant loci (P < 5 × 10^-8), of which eight are novel. Mendelian randomization analysis (MR) suggests a causal effect of depression on endometriosis. Combining gene-based association results across endometriosis and depression GWAS, we identified 22 genes with a genome-wide significant Fisher's combined P value (FCP_gene < 2.75 × 10^-6). Genes with a nominal gene-based association (P_gene < 0.05) were significantly enriched across endometriosis and depression (P_{binomial-test} = 2.90 × 10^-4). Also, genes overlapping the two traits at P_gene < 0.1 (P_{binomial-test} = 1.31 × 10^-5) were significantly enriched for the biological pathways 'cell-cell adhesion', 'inositol phosphate metabolism', 'Hippo-Merlin signaling dysregulation' and 'gastric mucosa abnormality'. These results reveal a shared genetic etiology for endometriosis and depression. Indeed, additional analyses found evidence of a causal association between each of endometriosis and depression and at least one abnormal condition of gastric mucosa. Our study confirms the comorbidity of endometriosis and depression, implicates links with gastric mucosa abnormalities in their causal pathways and reveals potential therapeutic targets for further investigation.

Prevalence and incidence of endometriosis in Australian women: a data linkage cohort study

OBJECTIVE

To estimate the prevalence and incidence of endometriosis among Australian women.

DESIGN

Population-based cohort study linked to administrative health records.

SETTING

Secondary analysis of seven surveys collected between 2000 and 2018 from a population-based cohort study.

POPULATION

A total of 13 508 Australian women, born 1973-78, from a prospective cohort study of 14 247 women conducted between 1996 and 2018.

METHODS

During 2000 and 2018, self-reported longitudinal survey data were linked to three administrative health databases to separately identify women with clinically confirmed or suspected endometriosis across the multiple data sources.

MAIN OUTCOME MEASURES

Prevalence and incidence of clinically confirmed endometriosis in the cohort were first estimated using national hospital data. Data were then combined with other administrative health databases and the survey data to capture all clinically confirmed and suspected diagnoses of endometriosis.

RESULTS

The cumulative prevalence of clinically confirmed endometriosis was 6.0% (95% CI 5.8-6.2%) by age 40-44 years. The cumulative prevalence increased to 11.4% (95% CI 11.1-11.7%) when adding diagnoses of clinically suspected endometriosis. Age-specific incidence estimates peaked to 6 per 1000 person-years at age 30-34 years.

CONCLUSIONS

Among 13 508 Australian women followed for 20 years, one in nine women had clinically confirmed or suspected endometriosis by the age of 44, with most diagnosed during their early thirties. Endometriosis is a significant public health issue requiring increased surveillance, clinical awareness and management. Efforts to expand knowledge on the aetiology of the disease and optimal methods for disease management are crucial to women's health.

TWEETABLE ABSTRACT

In a national study of 13 508 Australian women, one in nine women were diagnosed with endometriosis by age 44.

Genetic and environmental variation in educational attainment: an individual-based analysis of 28 twin cohorts

We investigated the heritability of educational attainment and how it differed between birth cohorts and cultural-geographic regions. A classical twin design was applied to pooled data from 28 cohorts representing 16 countries and including 193,518 twins with information on educational attainment at 25 years of age or older. Genetic factors explained the major part of individual differences in educational attainment (heritability: a2 = 0.43; 0.41-0.44), but also environmental variation shared by co-twins was substantial (c2 = 0.31; 0.30-0.33). The proportions of educational variation explained by genetic and shared environmental factors did not differ between Europe, North America and Australia, and East Asia. When restricted to twins 30 years or older to confirm finalized education, the heritability was higher in the older cohorts born in 1900-1949 (a2 = 0.44; 0.41-0.46) than in the later cohorts born in 1950-1989 (a2 = 0.38; 0.36-0.40), with a corresponding lower influence of common environmental factors (c2 = 0.31; 0.29-0.33 and c2 = 0.34; 0.32-0.36, respectively). In conclusion, both genetic and environmental factors shared by co-twins have an important influence on individual differences in educational attainment. The effect of genetic factors on educational attainment has decreased from the cohorts born before to those born after the 1950s.

Generation of immortalized human endometrial stromal cell lines with different endometriosis risk genotypes

Endometriotic lesions are composed in part of endometrial-like stromal cells, however, there is a shortage of immortalized human endometrial stromal cultures available for research. As genetic factors play a role in endometriosis risk, it is important that genotype is also incorporated into analysis of pathological mechanisms. Human telomerase reverse transcriptase (hTERT) immortalization (using Lenti-hTERT-green fluorescent protein virus) took place following genotype selection; 13 patients homozygous for either the risk or non-risk 'other' allele for one or more important endometriosis risk single nucleotide polymorphism on chromosome 1p36.12 (rs3820282, rs56318008, rs55938609, rs12037376, rs7521902 or rs12061255). Short tandem repeat DNA profiling validated that donor tissue matched that of the immortalized cell lines and confirmed that cultures were genetically novel. Expression of morphological markers (vimentin and cytokeratin) and key genes of interest (telomerase, estrogen and progesterone receptors and LINC00339) were examined and functional assays for cell proliferation, steroid hormone and inflammatory responses were performed for 7/13 cultures. All endometrial stromal cell lines maintained their fibroblast-like morphology (vimentin-positive) and homozygous endometriosis-risk genotype following introduction of hTERT. Furthermore, the new stromal cultures demonstrated positive and diverse responses to hormones (proliferation and decidualisation changes) and inflammation (dose-dependent response), while maintaining hormone receptor expression. In conclusion, we successfully developed a range of human endometrial stromal cell lines that carry important endometriosis-risk alleles. The wider implications of this approach go beyond advancing endometriosis research; these cell lines will be valuable tools for multiple endometrial pathologies offering a level of genetic and phenotypic diversity not previously available.

Obesity does not alter endometrial gene expression in women with endometriosis

RESEARCH QUESTION

Does obesity affect endometrial gene expression in women with endometriosis, specifically women with stage I disease?

DESIGN

Differential gene expression analysis was conducted on endometrium from women with and without endometriosis (n = 169). Women were diagnosed after surgical visualization and staged according to the revised American Society for Reproductive Medicine (stage I-IV). Women were grouped by body mass index (BMI) (kg/m²) as underweight, normal, pre-obese or obese. After accounting for menstrual cycle stage, endometrial gene expression was analysed by BMI (continuous and grouped) in women with endometriosis, and in non-endometriosis controls.

RESULTS

No significant interaction effect was found between BMI and endometriosis status on endometrial gene expression. We have previously reported that obese women with endometriosis have a reduced incidence of stage I disease; however, stratifying our analysis into stage I endometriosis versus combined II, III and IV endometriosis failed to reveal any differentially expressed endometrial genes between normal, pre-obese and obese patients.

CONCLUSIONS

Despite obesity having deleterious effects on endometrial gene expression in other gynaecological pathologies, e.g. endometrial cancer and polycystic ovary syndrome, our results do not support an association between BMI and altered endometrial gene expression in women with or without endometriosis.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.