GWAS meta-analyses clarify the genetics of cervical phenotypes and inform risk stratification for cervical cancer.

Genome-wide association studies (GWAS) have successfully identified associations for cervical cancer, but the underlying mechanisms of cervical biology and pathology remain uncharacterised. Our GWAS meta-analyses fill this gap, as we characterise the genetic architecture of cervical phenotypes, including cervical ectropion, cervicitis, cervical dysplasia, as well as up to 9229 cases and 490 304 controls for cervical cancer from diverse ancestries. Leveraging the latest computational methods and gene expression data, we refine the association signals for cervical cancer and propose potential causal variants and genes at each locus. We prioritise PAX8/PAX8-AS1, LINC00339, CDC42, CLPTM1L, HLA-DRB1 and GSDMB as the most likely candidate genes for cervical cancer signals, providing insights into cervical cancer pathogenesis and supporting the involvement of reproductive tract development, immune response and cellular proliferation/apoptosis. We construct a genetic risk score (GRS) that is associated with cervical cancer [hazard ratios (HR) = 3.1 (1.7-5.6) for the top 15% vs lowest 15% of individuals], and with other HPV- and immune-system-related diagnoses in a phenome-wide association study analysis. Our results propose valuable leads for further functional studies and present a GRS for cervical cancer that allows additional risk stratification and could potentially be used to personalise the conventional screening strategies for groups more susceptible to cervical cancer.

Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry.

BACKGROUND: Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes. RESULTS: The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata. CONCLUSIONS: Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.

Proteome-Wide Mendelian Randomization Identifies Causal Links Between Blood Proteins and Acute Pancreatitis.

BACKGROUND & AIMS: Acute pancreatitis (AP) is a complex disease and the leading cause of gastrointestinal disease-related hospital admissions. Few therapeutic options exist for AP prevention. Blood proteins with causal evidence may represent promising drug targets, but few have been causally linked with AP. Our objective was to identify blood proteins linked with AP by combining genome-wide association meta-analysis and proteome-wide Mendelian randomization (MR) studies. METHODS: We performed a genome-wide association meta-analysis totalling 10,630 patients with AP and 844,679 controls and a series of inverse-variance weighted MR analyses using cis-acting variants on 4719 blood proteins from the deCODE study (N = 35,559) and 4979 blood proteins from the Fenland study (N = 10,708). RESULTS: The meta-analysis identified genome-wide significant variants (P <5 × 10-8) at 5 loci (ABCG5/8, TWIST2, SPINK1, PRSS2 and MORC4). The proteome-wide MR analyses identified 68 unique blood proteins that may causally be associated with AP, including 29 proteins validated in both data sets. Functional annotation of these proteins confirmed expression of many proteins in metabolic tissues responsible for digestion and energy metabolism, such as the esophagus, adipose tissue, and liver as well as acinar cells of the pancreas. Genetic colocalization and investigations into the druggable genome also identified potential drug targets for AP. CONCLUSIONS: This large genome-wide association study meta-analysis for AP identified new variants linked with AP as well as several blood proteins that may be causally associated with AP. This study provides new information on the genetic architecture of this disease and identified pathways related to AP, which may be further explored as possible therapeutic targets for AP.

Circulating anti-Müllerian hormone levels in pre-menopausal women: novel genetic insights from a genome-wide association meta-analysis.

STUDY QUESTION: Can a genome-wide association study (GWAS) meta-analysis, including a large sample of young premenopausal women from a founder population from Northern Finland, identify novel genetic variants for circulating anti-Müllerian hormone (AMH) levels and provide insights into single-nucleotide polymorphism enrichment in different biological pathways and tissues involved in AMH regulation? SUMMARY ANSWER: The meta-analysis identified a total of six loci associated with AMH levels at P < 5 × 10-8, three of which were novel in or near CHEK2, BMP4, and EIF4EBP1, as well as highlighted significant enrichment in renal system vasculature morphogenesis, and the pituitary gland as the top associated tissue in tissue enrichment analysis. WHAT IS KNOWN ALREADY: AMH is expressed by preantral and small antral stage ovarian follicles in women, and variation in age-specific circulating AMH levels has been associated with several health conditions. However, the biological mechanisms underlying the association between health conditions and AMH levels are not yet fully understood. Previous GWAS have identified loci associated with AMH levels in pre-menopausal women, in or near MCM8, AMH, TEX41, and CDCA7. STUDY DESIGN, SIZE, DURATION: We performed a GWAS meta-analysis for circulating AMH level measurements in 9668 pre-menopausal women. PARTICIPANTS/MATERIALS, SETTING, METHODS: We performed a GWAS meta-analysis in which we combined 2619 AMH measurements (at age 31 years) from a prospective founder population cohort (Northern Finland Birth Cohort 1966, NFBC1966) with a previous GWAS meta-analysis that included 7049 pre-menopausal women (age range 15-48 years) (N = 9668). NFBC1966 AMH measurements were quantified using an automated assay. We annotated the genetic variants, combined different data layers to prioritize potential candidate genes, described significant pathways and tissues enriched by the GWAS signals, identified plausible regulatory roles using colocalization analysis, and leveraged publicly available summary statistics to assess genetic and phenotypic correlations with multiple traits. MAIN RESULTS AND THE ROLE OF CHANCE: Three novel genome-wide significant loci were identified. One of these is in complete linkage disequilibrium with c.1100delC in CHEK2, which is found to be 4-fold enriched in the Finnish population compared to other European populations. We propose a plausible regulatory effect of some of the GWAS variants linked to AMH, as they colocalize with GWAS signals associated with gene expression levels of BMP4, TEX41, and EIFBP41. Gene set analysis highlighted significant enrichment in renal system vasculature morphogenesis, and tissue enrichment analysis ranked the pituitary gland as the top association. LARGE SCALE DATA: The GWAS meta-analysis summary statistics are available for download from the GWAS Catalogue with accession number GCST90428625. LIMITATIONS, REASONS FOR CAUTION: This study only included women of European ancestry and the lack of sufficiently sized relevant tissue data in gene expression datasets hinders the assessment of potential regulatory effects in reproductive tissues. WIDER IMPLICATIONS OF THE FINDINGS: Our results highlight the increased power of founder populations and larger sample sizes to boost the discovery of novel trait-associated variants underlying variation in AMH levels, which aided the characterization of GWAS signals enrichment in different biological pathways and plausible genetic regulatory effects linked with AMH level variation for the first time. STUDY FUNDING/COMPETING INTEREST(S): This work has received funding from the European Union's Horizon 2020 Research and Innovation Programme under the MATER Marie Sklodowska-Curie Grant Agreement No. 813707 and Oulu University Scholarship Foundation and Paulon Säätiö Foundation. (N.P.-G.), Academy of Finland, Sigrid Jusélius Foundation, Novo Nordisk, University of Oulu, Roche Diagnostics (T.T.P.). This work was supported by the Estonian Research Council Grant 1911 (R.M.). J.R. was supported by the European Union's Horizon 2020 Research and Innovation Program under Grant Agreements No. 874739 (LongITools), 824989 (EUCAN-Connect), 848158 (EarlyCause), and 733206 (LifeCycle). U.V. was supported by the Estonian Research Council grant PRG (PRG1291). The NFBC1966 received financial support from University of Oulu Grant No. 24000692, Oulu University Hospital Grant No. 24301140, and ERDF European Regional Development Fund Grant No. 539/2010 A31592. T.T.P. has received grants from Roche, Perkin Elmer, and honoraria for scientific presentations from Gedeon Richter, Exeltis, Astellas, Roche, Stragen, Astra Zeneca, Merck, MSD, Ferring, Duodecim, and Ajaton Terveys. For all other authors, there are no competing interests.

The endometrial transcriptome of infertile women with and without implantation failure.

INTRODUCTION: Implantation failure after transferring morphologically "good-quality" embryos in in vitro fertilization/intracytoplasmic sperm injection (IVF/ICSI) may be explained by impaired endometrial receptivity. Analyzing the endometrial transcriptome analysis may reveal the underlying processes and could help in guiding prognosis and using targeted interventions for infertility. This exploratory study investigated whether the endometrial transcriptome profile was associated with short-term or long-term implantation outcomes (ie success or failure). MATERIAL AND METHODS: Mid-luteal phase endometrial biopsies of 107 infertile women with one full failed IVF/ICSI cycle, obtained within an endometrial scratching trial, were subjected to RNA-sequencing and differentially expressed genes analysis with covariate adjustment (age, body mass index, luteinizing hormone [LH]-day). Endometrial transcriptomes were compared between implantation failure and success groups in the short term (after the second fresh IVF/ICSI cycle) and long term (including all fresh and frozen cycles within 12 months). The short-term analysis included 85/107 women (33 ongoing pregnancy vs 52 no pregnancy), excluding 22/107 women. The long-term analysis included 46/107 women (23 'fertile' group, ie infertile women with a live birth after ≤3 embryos transferred vs 23 recurrent implantation failure group, ie no live birth after ≥3 good quality embryos transferred), excluding 61/107 women not fitting these categories. As both analyses drew from the same pool of 107 samples, there was some sample overlap. Additionally, cell type enrichment scores and endometrial receptivity were analyzed, and an endometrial development pseudo-timeline was constructed to estimate transcriptomic deviations from the optimum receptivity day (LH + 7), denoted as ΔWOI (window of implantation). RESULTS: There were no significantly differentially expressed genes between implantation failure and success groups in either the short-term or long-term analyses. Principal component analysis initially showed two clusters in the long-term analysis, unrelated to clinical phenotype and no longer distinct following covariate adjustment. Cell type enrichment scores did not differ significantly between groups in both analyses. However, endometrial receptivity analysis demonstrated a potentially significant displacement of the WOI in the non-pregnant group compared with the ongoing pregnant group in the short-term analysis. CONCLUSIONS: No distinct endometrial transcriptome profile was associated with either implantation failure or success in infertile women. However, there may be differences in the extent to which the WOI is displaced.

Genome-wide Study Identifies Association between HLA-B∗55:01 and Self-Reported Penicillin Allergy.

Hypersensitivity reactions to drugs are often unpredictable and can be life threatening, underscoring a need for understanding their underlying mechanisms and risk factors. The extent to which germline genetic variation influences the risk of commonly reported drug allergies such as penicillin allergy remains largely unknown. We extracted data from the electronic health records of more than 600,000 participants from the UK, Estonian, and Vanderbilt University Medical Center's BioVU biobanks to study the role of genetic variation in the occurrence of self-reported penicillin hypersensitivity reactions. We used imputed SNP to HLA typing data from these cohorts to further fine map the human leukocyte antigen (HLA) association and replicated our results in 23andMe's research cohort involving a total of 1.12 million individuals. Genome-wide meta-analysis of penicillin allergy revealed two loci, including one located in the HLA region on chromosome 6. This signal was further fine-mapped to the HLA-B∗55:01 allele (OR 1.41 95% CI 1.33-1.49, p value 2.04 × 10-31) and confirmed by independent replication in 23andMe's research cohort (OR 1.30 95% CI 1.25-1.34, p value 1.00 × 10-47). The lead SNP was also associated with lower lymphocyte counts and in silico follow-up suggests a potential effect on T-lymphocytes at HLA-B∗55:01. We also observed a significant hit in PTPN22 and the GWAS results correlated with the genetics of rheumatoid arthritis and psoriasis. We present robust evidence for the role of an allele of the major histocompatibility complex (MHC) I gene HLA-B in the occurrence of penicillin allergy.

Genome-wide association study meta-analysis supports association between MUC1 and ectopic pregnancy.

STUDY QUESTION: Can we identify genetic variants associated with ectopic pregnancy by undertaking the first genome-wide association study (GWAS) leveraging two large-scale biobank initiatives? SUMMARY ANSWER: We identified two novel genome-wide significant associations with ectopic pregnancy, highlighting MUC1 (mucin 1) as the most plausible affected gene. WHAT IS KNOWN ALREADY: Ectopic pregnancy is an important cause of maternal morbidity and mortality worldwide. Despite being a common early pregnancy complication, the genetic predisposition to this condition remains understudied and no large scale genetic studies have been performed so far. STUDY DESIGN, SIZE, DURATION: A GWAS meta-analysis including 7070 women with ectopic pregnancy and 248 810 controls from Estonian Biobank and the FinnGen study. PARTICIPANTS/MATERIALS, SETTING, METHODS: We identified ectopic pregnancy cases from national registers by ICD (International Classification of Disease) codes (ICD-10 O00), and all remaining women were considered controls. We carried out standard GWAS meta-analysis and additionally annotated GWAS signals, analysed co-localization with quantitative trait loci, estimated genetic correlations and identified associated phenotypes to characterize the genetic signals, as well as to analyse the genetic and phenotypic relationships with the condition. MAIN RESULTS AND THE ROLE OF CHANCE: We identified two genome-wide significant loci on chromosomes 1 (rs4971091, P = 5.32×10-9) and 10 (rs11598956, P = 2.41×10-8) potentially associated with ectopic pregnancy. Follow-up analyses propose MUC1, which codes for an epithelial glycoprotein with an important role in barrier function, as the most likely candidate gene for the association on chromosome 1. We also characterize the phenotypic and genetic correlations with other phenotypes, identifying a genetic correlation with smoking and diseases of the (genito)urinary and gastrointestinal system, and phenotypic correlations with various reproductive health diagnoses, reflecting the previously known epidemiological associations. LARGE SCALE DATA: The GWAS meta-analysis summary statistics are available from the GWAS Catalogue (GCST90272883). LIMITATIONS, REASONS FOR CAUTION: The main limitation is that the findings are based on European-based ancestry populations, with limited data on other populations, and we only captured maternal genomes. Additionally, further larger meta-analysis or independent studies are needed to validate these findings. WIDER IMPLICATIONS OF THE FINDINGS: This study encourages the use of large-scale genetic datasets to unravel genetic factors linked to ectopic pregnancy, which is difficult to study in experimental settings. Increased sample size might bring additional genetic factors associating with ectopic pregnancy and inform its heritability. Altogether, our results provide more insight into the biology of ectopic pregnancy and, accordingly, the biological processes governing embryo implantation. STUDY FUNDING/COMPETING INTEREST(S): N.P.G. was supported by MATER Marie Sklodowska-Curie which received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No. 813707. This study was funded by European Union through the European Regional Development Fund Project No. 2014-2020.4.01.15-0012 GENTRANSMED. Computations were performed in the High-Performance Computing Center of University of Tartu. The authors declare no competing interests.

Uniting biobank resources reveals novel genetic pathways modulating susceptibility for atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is a common chronic inflammatory skin disease with high heritability. Previous genome-wide association studies have identified several loci predisposing to AD. These findings explain approximately 30% of the variance in AD susceptibility, suggesting that further work is required to fully understand the genetic underpinnings. OBJECTIVE: We sought to gain additional understanding of the genetic contribution to AD risk by using biobank resources. METHODS: We completed a genome-wide meta-analysis of AD in 796,661 individuals (Ncases = 22,474) from the FinnGen study, the Estonian Biobank, and the UK Biobank. We further performed downstream in silico analyses to characterize the risk variants at the novel loci. RESULTS: We report 30 loci associating with AD (P < 5 × 10-8), 5 of which are novel. In 2 of the novel loci, we identified missense mutations with deleterious predictions in desmocollin 1 and serpin family B member 7, genes encoding proteins crucial to epidermal strength and integrity. CONCLUSIONS: These findings elucidate novel genetic pathways involved in AD pathophysiology. The likely involvement of desmocollin 1 and serpin family B member 7 in AD pathogenesis may offer opportunities for the development of novel treatment strategies for AD in the future.

GWAS Identifies Risk Locus for Erectile Dysfunction and Implicates Hypothalamic Neurobiology and Diabetes in Etiology.

Erectile dysfunction (ED) is a common condition affecting more than 20% of men over 60 years, yet little is known about its genetic architecture. We performed a genome-wide association study of ED in 6,175 case subjects among 223,805 European men and identified one locus at 6q16.3 (lead variant rs57989773, OR 1.20 per C-allele; p = 5.71 × 10-14), located between MCHR2 and SIM1. In silico analysis suggests SIM1 to confer ED risk through hypothalamic dysregulation. Mendelian randomization provides evidence that genetic risk of type 2 diabetes mellitus is a cause of ED (OR 1.11 per 1-log unit higher risk of type 2 diabetes). These findings provide insights into the biological underpinnings and the causes of ED and may help prioritize the development of future therapies for this common disorder.

Leveraging Northern European population history: novel low-frequency variants for polycystic ovary syndrome.

STUDY QUESTION: Can we identify novel variants associated with polycystic ovary syndrome (PCOS) by leveraging the unique population history of Northern Europe? SUMMARY ANSWER: We identified three novel genome-wide significant associations with PCOS, with two putative independent causal variants in the checkpoint kinase 2 (CHEK2) gene and a third in myosin X (MYO10). WHAT IS KNOWN ALREADY: PCOS is a common, complex disorder with unknown aetiology. While previous genome-wide association studies (GWAS) have mapped several loci associated with PCOS, the analysis of populations with unique population history and genetic makeup has the potential to uncover new low-frequency variants with larger effects. STUDY DESIGN, SIZE, DURATION: A population-based case-control GWAS was carried out. PARTICIPANTS/MATERIALS, SETTING, METHODS: We identified PCOS cases from national registers by ICD codes (ICD-10 E28.2, ICD-9 256.4, or ICD-8 256.90), and all remaining women were considered controls. We then conducted a three-stage case-control GWAS: in the discovery phase, we had a total of 797 cases and 140 558 controls from the FinnGen study. For validation, we used an independent dataset from the Estonian Biobank, including 2812 cases and 89 230 controls. Finally, we performed a joint meta-analysis of 3609 cases and 229 788 controls from both cohorts. Additionally, we reran the association analyses including BMI as a covariate, with 2169 cases and 160 321 controls from both cohorts. MAIN RESULTS AND THE ROLE OF CHANCE: Two out of the three novel genome-wide significant variants associating with PCOS, rs145598156 (P = 3.6×10-8, odds ratio (OR) = 3.01 [2.02-4.50] minor allele frequency (MAF) = 0.005) and rs182075939 (P = 1.9×10-16, OR = 1.69 [1.49-1.91], MAF = 0.04), were found to be enriched in the Finnish and Estonian populations and are tightly linked to a deletion c.1100delC (r2 = 0.95) and a missense I157T (r2 = 0.83) in CHEK2. The third novel association is a common variant near MYO10 (rs9312937, P = 1.7 × 10-8, OR = 1.16 [1.10-1.23], MAF = 0.44). We also replicated four previous reported associations near the genes Erb-B2 Receptor Tyrosine Kinase 4 (ERBB4), DENN Domain Containing 1A (DENND1A), FSH Subunit Beta (FSHB) and Zinc Finger And BTB Domain Containing 16 (ZBTB16). When adding BMI as a covariate only one of the novel variants remained genome-wide significant in the meta-analysis (the EstBB lead signal in CHEK2 rs182075939, P = 1.9×10-16, OR = 1.74 [1.5-2.01]) possibly owing to reduced sample size. LARGE SCALE DATA: The age- and BMI-adjusted GWAS meta-analysis summary statistics are available for download from the GWAS Catalog with accession numbers GCST90044902 and GCST90044903. LIMITATIONS, REASONS FOR CAUTION: The main limitation was the low prevalence of PCOS in registers; however, the ones with the diagnosis most likely represent the most severe cases. Also, BMI data were not available for all (63% for FinnGen, 76% for EstBB), and the biobank setting limited the accessibility of PCOS phenotypes and laboratory values. WIDER IMPLICATIONS OF THE FINDINGS: This study encourages the use of isolated populations to perform genetic association studies for the identification of rare variants contributing to the genetic landscape of complex diseases such as PCOS. STUDY FUNDING/COMPETING INTEREST(S): This work has received funding from the European Union's Horizon 2020 research and innovation programme under the MATER Marie Sklodowska-Curie grant agreement No. 813707 (N.P.-G., T.L., T.P.), the Estonian Research Council grant (PRG687, T.L.), the Academy of Finland grants 315921 (T.P.), 321763 (T.P.), 297338 (J.K.), 307247 (J.K.), 344695 (H.L.), Novo Nordisk Foundation grant NNF17OC0026062 (J.K.), the Sigrid Juselius Foundation project grants (T.L., J.K., T.P.), Finska Läkaresällskapet (H.L.) and Jane and Aatos Erkko Foundation (H.L.). The funders had no role in study design, data collection and analysis, publishing or preparation of the manuscript. The authors declare no conflicts of interest.

Correction: Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

[This corrects the article DOI: 10.1371/journal.pgen.1007813.].

Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

Polycystic ovary syndrome (PCOS) is a disorder characterized by hyperandrogenism, ovulatory dysfunction and polycystic ovarian morphology. Affected women frequently have metabolic disturbances including insulin resistance and dysregulation of glucose homeostasis. PCOS is diagnosed with two different sets of diagnostic criteria, resulting in a phenotypic spectrum of PCOS cases. The genetic similarities between cases diagnosed based on the two criteria have been largely unknown. Previous studies in Chinese and European subjects have identified 16 loci associated with risk of PCOS. We report a fixed-effect, inverse-weighted-variance meta-analysis from 10,074 PCOS cases and 103,164 controls of European ancestry and characterisation of PCOS related traits. We identified 3 novel loci (near PLGRKT, ZBTB16 and MAPRE1), and provide replication of 11 previously reported loci. Only one locus differed significantly in its association by diagnostic criteria; otherwise the genetic architecture was similar between PCOS diagnosed by self-report and PCOS diagnosed by NIH or non-NIH Rotterdam criteria across common variants at 13 loci. Identified variants were associated with hyperandrogenism, gonadotropin regulation and testosterone levels in affected women. Linkage disequilibrium score regression analysis revealed genetic correlations with obesity, fasting insulin, type 2 diabetes, lipid levels and coronary artery disease, indicating shared genetic architecture between metabolic traits and PCOS. Mendelian randomization analyses suggested variants associated with body mass index, fasting insulin, menopause timing, depression and male-pattern balding play a causal role in PCOS. The data thus demonstrate 3 novel loci associated with PCOS and similar genetic architecture for all diagnostic criteria. The data also provide the first genetic evidence for a male phenotype for PCOS and a causal link to depression, a previously hypothesized comorbid disease. Thus, the genetics provide a comprehensive view of PCOS that encompasses multiple diagnostic criteria, gender, reproductive potential and mental health.

Large-scale meta-analysis highlights the hypothalamic-pituitary-gonadal axis in the genetic regulation of menstrual cycle length.

The normal menstrual cycle requires a delicate interplay between the hypothalamus, pituitary and ovary. Therefore, its length is an important indicator of female reproductive health. Menstrual cycle length has been shown to be partially controlled by genetic factors, especially in the follicle-stimulating hormone beta-subunit (FSHB) locus. A genome-wide association study meta-analysis of menstrual cycle length in 44 871 women of European ancestry confirmed the previously observed association with the FSHB locus and identified four additional novel signals in, or near, the GNRH1, PGR, NR5A2 and INS-IGF2 genes. These findings not only confirm the role of the hypothalamic-pituitary-gonadal axis in the genetic regulation of menstrual cycle length but also highlight potential novel local regulatory mechanisms, such as those mediated by IGF2.

Endometrial receptivity revisited: endometrial transcriptome adjusted for tissue cellular heterogeneity.

STUDY QUESTION: Does cellular composition of the endometrial biopsy affect the gene expression profile of endometrial whole-tissue samples? SUMMARY ANSWER: The differences in epithelial and stromal cell proportions in endometrial biopsies modify the whole-tissue gene expression profiles and affect the results of differential expression analyses. WHAT IS ALREADY KNOWN: Each cell type has its unique gene expression profile. The proportions of epithelial and stromal cells vary in endometrial tissue during the menstrual cycle, along with individual and technical variation due to the method and tools used to obtain the tissue biopsy. STUDY DESIGN, SIZE, DURATION: Using cell-population specific transcriptome data and computational deconvolution approach, we estimated the epithelial and stromal cell proportions in whole-tissue biopsies taken during early secretory and mid-secretory phases. The estimated cellular proportions were used as covariates in whole-tissue differential gene expression analysis. Endometrial transcriptomes before and after deconvolution were compared and analysed in biological context. PARTICIPANTS/MATERIAL, SETTING, METHODS: Paired early- and mid-secretory endometrial biopsies were obtained from 35 healthy, regularly cycling, fertile volunteers, aged 23-36 years, and analysed by RNA sequencing. Differential gene expression analysis was performed using two approaches. In one of them, computational deconvolution was applied as an intermediate step to adjust for the proportions of epithelial and stromal cells in the endometrial biopsy. The results were then compared to conventional differential expression analysis. Ten paired endometrial samples were analysed with qPCR to validate the results. MAIN RESULTS AND THE ROLE OF CHANCE: The estimated average proportions of stromal and epithelial cells in early secretory phase were 65% and 35%, and during mid-secretory phase, 46% and 54%, respectively, correlating well with the results of histological evaluation (r = 0.88, P = 1.1 × 10-6). Endometrial tissue transcriptomic analysis showed that approximately 26% of transcripts (n = 946) differentially expressed in receptive endometrium in cell-type unadjusted analysis also remain differentially expressed after adjustment for biopsy cellular composition. However, the other 74% (n = 2645) become statistically non-significant after adjustment for biopsy cellular composition, underlining the impact of tissue heterogeneity on differential expression analysis. The results suggest new mechanisms involved in endometrial maturation, involving genes like LINC01320, SLC8A1 and GGTA1P, described for the first time in context of endometrial receptivity. LARGE-SCALE DATA: The RNA-seq data presented in this study is deposited in the Gene Expression Omnibus database with accession number GSE98386. LIMITATIONS REASONS FOR CAUTION: Only dominant endometrial cell types were considered in gene expression profile deconvolution; however, other less frequent endometrial cell types also contribute to the whole-tissue gene expression profile. WIDER IMPLICATIONS OF THE FINDINGS: The better understanding of molecular processes during transition from pre-receptive to receptive endometrium serves to improve the effectiveness and personalization of assisted reproduction protocols. Biopsy cellular composition should be taken into account in future endometrial 'omics' studies, where tissue heterogeneity could potentially influence the results. STUDY FUNDING/COMPETING INTEREST(S): This study was funded by: Estonian Ministry of Education and Research (grant IUT34-16); Enterprise Estonia (EU48695); the EU-FP7 Eurostars program (NOTED, EU41564); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (SARM, EU324509); Horizon 2020 innovation program (WIDENLIFE, EU692065); MSCA-RISE-2015 project MOMENDO (No 691058) and the Miguel Servet Program Type I of Instituto de Salud Carlos III (CP13/00038); Spanish Ministry of Economy, Industry and Competitiveness (MINECO) and European Regional Development Fund (FEDER): grants RYC-2016-21199 and ENDORE SAF2017-87526. Authors confirm no competing interests.

Identifying Atopic Dermatitis Risk Loci in 1,094,060 Individuals with Subanalysis of Disease Severity and Onset.

Atopic dermatitis (AD) is a common inflammatory skin disease highly attributable to genetic factors. In this study, we report results from a genome-wide meta-analysis of AD in 37,541 cases and 1,056,519 controls with data from the FinnGen project, the Estonian Biobank, the UK Biobank, the EAGLE Consortium, and the BioBank Japan. We detected 77 independent AD-associated loci, of which 10 were, to our knowledge, previously unreported. The associated loci showed enrichment in various immune regulatory processes. We further performed subgroup analyses of mild and severe AD and of early- and late-onset AD, with data from the FinnGen project. Fifty-five of the 79 tested variants in the associated loci showed larger effect estimates for severe than for mild AD as determined through administered treatment. The age of onset, as determined by the first hospital visit with AD diagnosis, was lower in patients with particular AD-risk alleles. Our findings add to the knowledge of the genetic background of AD and may underlie the development of new therapeutic strategies.

Uncovering the shared genetic components of thyroid disorders and reproductive health.

OBJECTIVE: The aim of the study is to map the shared genetic component and relationships between thyroid and reproductive health traits to improve the understanding of the interplay between those domains. DESIGN: A large-scale genetic analysis of thyroid traits (hyper- and hypothyroidism, and thyroid-stimulating hormone levels) was conducted in up to 743 088 individuals of European ancestry from various cohorts. METHODS: We evaluated genetic associations using genome-wide association study (GWAS) meta-analysis, GWAS Catalog lookup, gene prioritization, mouse phenotype lookup, and genetic correlation analysis. RESULTS: GWAS meta-analysis results for thyroid phenotypes showed that 50 lead variants out of 253 (including 5/52 of the novel hits) were linked to reproductive health in previous literature. Genetic correlation analyses revealed significant correlations between hypothyroidism and reproductive phenotypes. The results showed that 31.9% of thyroid-associated genes also had an impact on reproductive phenotypes, with the most affected functions being related to genitourinary tract issues. CONCLUSIONS: The study discovers novel genetic loci linked to thyroid phenotypes and highlights the shared genetic determinants between thyroid function and reproductive health, providing evidence for the genetic pleiotropy and shared biological mechanisms between these traits in both sexes.

Genome-wide association study in Estonia reveals importance of vaginal epithelium associated genes in case of recurrent vaginitis.

Recurrent vaginitis is a leading reason for visiting a gynaecologist, with bacterial vaginosis (BV) and vulvovaginal candidiasis (VVC) being the most common diagnoses. Reasons and mechanisms behind their recurrent nature are poorly understood. We conducted a genome-wide association study (GWAS) to find possible genetic risk factors for recurrent vaginitis using data from a large population-based biobank, the Estonian Biobank. The study included 6870 cases (at least two episodes of vaginitis) and 5945 controls (no vaginitis episodes). GWAS approach included single marker and gene-based analyses, followed by functional annotation of associated variants and candidate gene mapping.In single marker analysis, one statistically significant (P = 7.8 × 10-9) variant rs1036732378 was identified on chromosome 10. The gene-based association analysis identified one gene, KRT6A, that exceeded the recommended significance threshold (P = 2.6 × 10-6). This is a member of the keratin protein family and is expressed during differentiation in epithelial tissues.Functional mapping and annotation of genetic associations by using adjusted significance level identified 22 potential risk loci that may be associated with recurrent vaginitis phenotype. Comparison of our results with previous studies provided nominal support for LBP (associated with immune response to vaginal bacteria) and PRKCH genes (possible role in keratinocyte differentiation and susceptibility to candidiasis).In conclusion, this study is the first highlighting a potential role of the vaginal epithelium in recurrent vaginitis.

Genetic Study of Psoriasis Highlights its Close Link with Socioeconomic Status and Affective Symptoms.

Psoriasis is an inflammatory skin disease with an estimated heritability of around 70%. Previous GWASs have detected several risk loci for psoriasis. To further improve the understanding of the genetic risk factors impacting the disease, we conducted a discovery GWAS in FinnGen and a subsequent replication and meta-analysis with data from the Estonian Biobank and the UK Biobank; the study sample included 925,649 individuals (22,659 cases and 902,990 controls), the largest sample for psoriasis yet. In addition, we conducted downstream analyses to find out more about psoriasis' cross-trait genetic correlations and causal relationships. We report 6 risk loci, which, to our knowledge, are previously unreported, most of which harbor genes related to NF-κB signaling pathway and overall immunity. Genetic correlations highlight the relationship between psoriasis and smoking, higher body weight, and lower education level. In addition, we report causal relationships between psoriasis and mood symptoms as well as 2-directioned causal relationship between psoriasis and lower education level. Our results provide further knowledge on psoriasis risk factors, which may be useful in the development of future treatment strategies.