STUDY QUESTION Can additional genetic variants for circulating anti-Müllerian hormone (AMH) levels be identified through a genome-wide association study (GWAS) meta-analysis including a large sample of premenopausal women? SUMMARY ANSWER We identified four loci associated with AMH levels at P < 5 × 10−8: the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 and CDCA7. WHAT IS KNOWN ALREADY AMH is expressed by antral stage ovarian follicles in women, and variation in age-specific circulating AMH levels has been associated with disease outcomes. However, the physiological mechanisms underlying these AMH-disease associations are largely unknown. STUDY DESIGN, SIZE, DURATION We performed a GWAS meta-analysis in which we combined summary statistics of a previous AMH GWAS with GWAS data from 3705 additional women from three different cohorts. PARTICIPANTS/MATERIALS, SETTING, METHODS In total, we included data from 7049 premenopausal female participants of European ancestry. The median age of study participants ranged from 15.3 to 48 years across cohorts. Circulating AMH levels were measured in either serum or plasma samples using different ELISA assays. Study-specific analyses were adjusted for age at blood collection and population stratification, and summary statistics were meta-analysed using a standard error-weighted approach. Subsequently, we functionally annotated GWAS variants that reached genome-wide significance (P < 5 × 10−8). We also performed a gene-based GWAS, pathway analysis and linkage disequilibrium score regression and Mendelian randomization (MR) analyses. MAIN RESULTS AND THE ROLE OF CHANCE We identified four loci associated with AMH levels at P < 5 × 10−8: the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 and CDCA7. The strongest signal was a missense variant in the AMH gene (rs10417628). Most prioritized genes at the other three identified loci were involved in cell cycle regulation. Genetic correlation analyses indicated a strong positive correlation among single nucleotide polymorphisms for AMH levels and for age at menopause (rg = 0.82, FDR = 0.003). Exploratory two-sample MR analyses did not support causal effects of AMH on breast cancer or polycystic ovary syndrome risk, but should be interpreted with caution as they may be underpowered and the validity of genetic instruments could not be extensively explored. LARGE SCALE DATA The full AMH GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION Whilst this study doubled the sample size of the most recent GWAS, the statistical power is still relatively low. As a result, we may still lack power to identify more genetic variants for AMH and to determine causal effects of AMH on, for example, breast cancer. Also, follow-up studies are needed to investigate whether the signal for the AMH gene is caused by reduced AMH detection by certain assays instead of actual lower circulating AMH levels. WIDER IMPLICATIONS OF THE FINDINGS Genes mapped to the MCM8, TEX41 and CDCA7 loci are involved in the cell cycle and processes such as DNA replication and apoptosis. The mechanism underlying their associations with AMH may affect the size of the ovarian follicle pool. Altogether, our results provide more insight into the biology of AMH and, accordingly, the biological processes involved in ovarian ageing. STUDY FUNDING/COMPETING INTEREST(S) Nurses’ Health Study and Nurses’ Health Study II were supported by research grants from the National Institutes of Health (CA172726, CA186107, CA50385, CA87969, CA49449, CA67262, CA178949). The UK Medical Research Council and Wellcome (217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This publication is the work of the listed authors, who will serve as guarantors for the contents of this article. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). Funding for the collection of genotype and phenotype data used here was provided by the British Heart Foundation (SP/07/008/24066), Wellcome (WT092830M and WT08806) and UK Medical Research Council (G1001357). M.C.B., A.L.G.S. and D.A.L. work in a unit that is funded by the University of Bristol and UK Medical Research Council (MC_UU_00011/6). M.C.B.’s contribution to this work was funded by a UK Medical Research Council Skills Development Fellowship (MR/P014054/1) and D.A.L. is a National Institute of Health Research Senior Investigator (NF-0616-10102). A.L.G.S. was supported by the study of Dynamic longitudinal exposome trajectories in cardiovascular and metabolic non-communicable diseases (H2020-SC1-2019-Single-Stage-RTD, project ID 874739). The Doetinchem Cohort Study was financially supported by the Ministry of Health, Welfare and Sports of the Netherlands. The funder had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. Ansh Labs performed the AMH measurements for the Doetinchem Cohort Study free of charge. Ansh Labs was not involved in the data analysis, interpretation or reporting, nor was it financially involved in any aspect of the study. R.M.G.V. was funded by the Honours Track of MSc Epidemiology, University Medical Center Utrecht with a grant from the Netherlands Organization for Scientific Research (NWO) (022.005.021). The Study of Women's Health Across the Nation (SWAN) has grant support from the National Institutes of Health (NIH), DHHS, through the National Institute on Aging (NIA), the National Institute of Nursing Research (NINR) and the NIH Office of Research on Women’s Health (ORWH) (U01NR004061; U01AG012505, U01AG012535, U01AG012531, U01AG012539, U01AG012546, U01AG012553, U01AG012554, U01AG012495). The SWAN Genomic Analyses and SWAN Legacy have grant support from the NIA (U01AG017719). The Generations Study was funded by Breast Cancer Now and the Institute of Cancer Research (ICR). The ICR acknowledges NHS funding to the NIHR Biomedical Research Centre. The content of this manuscript is solely the responsibility of the authors and does not necessarily represent official views of the funders. The Sister Study was funded by the Intramural Research Program of the National Institutes of Health (NIH), National Institute of Environmental Health Sciences (Z01-ES044005 to D.P.S.); the AMH assays were supported by the Avon Foundation (02-2012-065 to H.B. Nichols and D.P.S.). The breast cancer genome-wide association analyses were supported by the Government of Canada through Genome Canada and the Canadian Institutes of Health Research, the ‘Ministère de l’Économie, de la Science et de l’Innovation du Québec’ through Genome Québec and grant PSR-SIIRI-701, The National Institutes of Health (U19 CA148065, X01HG007492), Cancer Research UK (C1287/A10118, C1287/A16563, C1287/A10710) and The European Union (HEALTH-F2-2009-223175 and H2020 633784 and 634935). All studies and funders are listed in Michailidou et al. (Nature, 2017). F.J.M.B. has received fees and grant support from Merck Serono and Ferring BV. D.A.L. has received financial support from several national and international government and charitable funders as well as from Medtronic Ltd and Roche Diagnostics for research that is unrelated to this study. N.S. is scientific consultant for Ansh Laboratories. The other authors declare no competing interests.

Aims/hypothesis Given its role in ovarian follicle development, circulating anti-Müllerian hormone (AMH) is considered to be a marker of reproductive ageing. Although accelerated reproductive ageing has been associated with a higher risk of type 2 diabetes, research on the relationship between AMH and type 2 diabetes risk is scarce. Therefore, we aimed to investigate whether age-specific AMH levels and age-related AMH trajectories are associated with type 2 diabetes risk in women. Methods We measured AMH in repeated plasma samples from 3293 female participants (12,460 samples in total), aged 20–59 years at recruitment, from the Doetinchem Cohort Study, a longitudinal study with follow-up visits every 5 years. We calculated age-specific AMH tertiles at baseline to account for the strong AMH–age correlation. Cox proportional hazards models adjusted for confounders were used to assess the association between baseline age-specific AMH tertiles and incident type 2 diabetes. We applied linear mixed models to compare age-related AMH trajectories for women who developed type 2 diabetes with trajectories for women who did not develop diabetes. Results During a median follow-up of 20 years, 163 women developed type 2 diabetes. Lower baseline age-specific AMH levels were associated with a higher type 2 diabetes risk (HRT2vsT3 1.24 [95% CI 0.81, 1.92]; HRT1vsT3 1.62 [95% CI 1.06, 2.48]; ptrend = 0.02). These findings seem to be supported by predicted AMH trajectories, which suggested that plasma AMH levels were lower at younger ages in women who developed type 2 diabetes compared with women who did not. The trajectories also suggested that AMH levels declined at a slower rate in women who developed type 2 diabetes, although differences in trajectories were not statistically significant. Conclusions/interpretation We observed that lower age-specific AMH levels were associated with a higher risk of type 2 diabetes in women. Longitudinal analyses did not show clear evidence of differing AMH trajectories between women who developed type 2 diabetes compared with women who did not, possibly because these analyses were underpowered. Further research is needed to investigate whether AMH is part of the biological mechanism explaining the association between reproductive ageing and type 2 diabetes.

STUDY QUESTION Do polymorphisms in the anti-Müllerian hormone (AMH) promoter have an effect on AMH levels in patients with polycystic ovary syndrome (PCOS)? SUMMARY ANSWER We have identified a novel AMH promoter polymorphism rs10406324 that is associated with lower serum AMH levels and is suggested to play a role in the mechanism of regulation of AMH gene expression in women. WHAT IS KNOWN ALREADY Follicle number is positively correlated with serum AMH levels, reflected by elevated AMH levels in women with PCOS. In addition, it is suggested that AMH production per follicle is higher in women with PCOS than in normo-ovulatory women, implying an altered regulation of AMH in PCOS. STUDY DESIGN, SIZE, DURATION A discovery cohort of 655 PCOS women of Northern European ancestry and both an internal and external validation PCOS cohort (n = 458 and n = 321, respectively) were included in this study. Summary-level data of an AMH genome-wide association study meta-analysis including 7049 normo-ovulatory women was included as a control cohort. A genetic approach was taken through association analysis and in silico analysis of the associated variants in the AMH promoter. In vitro analysis was performed to investigate the functional mechanisms. PARTICIPANTS/MATERIALS, SETTING, METHODS All common two-allelic single-nucleotide polymorphisms (SNPs) in the region Chr19:2 245 353–2 250 827 bp (Build 37) were selected for the analysis. Linear regression analyses were performed to determine the association between SNPs in the AMH promoter region and serum AMH levels. For the in silico analysis, the webtools ‘HaploReg’ v4.1 for ENCODE prediction weight matrices and ‘atSNP’ were used. In vitro analysis was performed using KK1 cells, a mouse granulosa cell line and COV434 cells, a human granulosa tumor cell line. Cells were transfected with the reference or the variant human AMH promoter reporter construct together with several transcription factors (TFs). Dual-Glo® Luciferase Assay was performed to measure the luciferase activity. MAIN RESULTS AND THE ROLE OF CHANCE Polymorphism rs10406324 was significantly associated with serum AMH levels in all three PCOS cohorts. Carriers of the minor allele G had significantly lower log-transformed serum AMH levels compared to non-carriers (P = 8.58 × 10−8, P = 1.35 × 10−3 and P = 1.24 × 10−3, respectively). This result was validated in a subsequent meta-analysis (P = 3.24 × 10−12). Interestingly, rs10406324 was not associated with follicle count, nor with other clinical traits. Also, in normo-ovulatory women, the minor allele of this variant was associated with lower serum AMH levels (P = 1.04 × 10−5). These findings suggest that polymorphism rs10406324 plays a role in the regulation of AMH expression, irrespective of clinical background. In silico analysis suggested a decreased binding affinity of the TFs steroidogenenic factor 1, estrogen-related receptor alpha and glucocorticoid receptor to the minor allele G variant, however in vitro analysis did not show a difference in promoter activity between the A and G allele. LIMITATIONS, REASONS FOR CAUTION Functional analyses were performed in a mouse and a human granulosa cell line using an AMH promoter reporter construct. This may have limited assessment of the impact of the polymorphism on higher order chromatin structures. Human granulosa cells generated from induced pluripotent stem cells, combined with gene editing, may provide a method to elucidate the exact mechanism behind the decrease in serum AMH levels in carriers of the −210 G allele. We acknowledge that the lack of follicle number in the external validation and the control cohort is a limitation of the paper. Although we observed that the association between rs10406324 and AMH levels was independent of follicle number in our discovery and internal validation PCOS cohorts, we cannot fully rule out that the observed effects on serum AMH levels are, in part, caused by differences in follicle number. WIDER IMPLICATIONS OF THE FINDINGS These results suggest that variations in serum AMH levels are not only caused by differences in follicle number but also by genetic factors. Therefore, the genetic context should be taken into consideration when assessing serum AMH levels in women. This may have clinical consequences when serum AMH levels are used as a marker for the polycystic ovarian morphology phenotype. STUDY FUNDING/COMPETING INTEREST(S) No external funding was used. J.S.E.L. has received consultancy fees from the following companies: Ferring, Roche Diagnostics and Ansh Labs and has received travel reimbursement from Ferring. J.A.V. has received royalties from AMH assays, paid to the institute/lab with no personal financial gain. The other authors declare no competing interests. TRIAL REGISTRATION NUMBER N/A.