Genome-wide association study meta-analysis identifies three novel loci for circulating anti-Müllerian hormone levels in women

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 (r(g) = 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

FTO genetic variants, dietary intake and body mass index: insights from 177 330 individuals

FTO is the strongest known genetic susceptibility locus for obesity. Experimental studies in animals suggest the potential roles of FTO in regulating food intake. The interactive relation among FTO variants, dietary intake and body mass index (BMI) is complex and results from previous often small-scale studies in humans are highly inconsistent. We performed large-scale analyses based on data from 177 330 adults (154 439 Whites, 5776 African Americans and 17 115 Asians) from 40 studies to examine: (i) the association between the FTO-rs9939609 variant (or a proxy single-nucleotide polymorphism) and total energy and macronutrient intake and (ii) the interaction between the FTO variant and dietary intake on BMI. The minor allele (A-allele) of the FTO-rs9939609 variant was associated with higher BMI in Whites (effect per allele = 0.34 [0.31, 0.37] kg/m2, P = 1.9 × 10−105), and all participants (0.30 [0.30, 0.35] kg/m2, P = 3.6 × 10−107). The BMI-increasing allele of the FTO variant showed a significant association with higher dietary protein intake (effect per allele = 0.08 [0.06, 0.10] %, P = 2.4 × 10−16), and relative weak associations with lower total energy intake (−6.4 [−10.1, −2.6] kcal/day, P = 0.001) and lower dietary carbohydrate intake (−0.07 [−0.11, −0.02] %, P = 0.004). The associations with protein (P = 7.5 × 10−9) and total energy (P = 0.002) were attenuated but remained significant after adjustment for BMI. We did not find significant interactions between the FTO variant and dietary intake of total energy, protein, carbohydrate or fat on BMI. Our findings suggest a positive association between the BMI-increasing allele of FTO variant and higher dietary protein intake and offer insight into potential link between FTO, dietary protein intake and adiposit

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10−8) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10−8). The top IBC association for SBP was rs2012318 (P= 6.4 × 10−6) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10−6) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexit

A Meta-analysis of Gene Expression Signatures of Blood Pressure and Hypertension

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p&lt;0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%–9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk.

Blood pressure is a heritable trait influenced by several biological pathways and responsive to environmental stimuli. Over one billion people worldwide have hypertension (≥140 mm Hg systolic blood pressure or  ≥90 mm Hg diastolic blood pressure). Even small increments in blood pressure are associated with an increased risk of cardiovascular events. This genome-wide association study of systolic and diastolic blood pressure, which used a multi-stage design in 200,000 individuals of European descent, identified sixteen novel loci: six of these loci contain genes previously known or suspected to regulate blood pressure (GUCY1A3-GUCY1B3, NPR3-C5orf23, ADM, FURIN-FES, GOSR2, GNAS-EDN3); the other ten provide new clues to blood pressure physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke and coronary artery disease, but not kidney disease or kidney function. We also observed associations with blood pressure in East Asian, South Asian and African ancestry individuals. Our findings provide new insights into the genetics and biology of blood pressure, and suggest potential novel therapeutic pathways for cardiovascular disease prevention

Causal effect of plasminogen activator inhibitor type 1 on coronary heart disease

Background--Plasminogen activator inhibitor type 1 (PAI-1) plays an essential role in the fibrinolysis system and thrombosis. Population studies have reported that blood PAI-1 levels are associated with increased risk of coronary heart disease (CHD). However, it is unclear whether the association reflects a causal influence of PAI-1 on CHD risk. Methods and Results--To evaluate the association between PAI-1 and CHD, we applied a 3-step strategy. First, we investigated the observational association between PAI-1 and CHD incidence using a systematic review based on a literature search for PAI-1 and CHD studies. Second, we explored the causal association between PAI-1 and CHD using a Mendelian randomization approach using summary statistics from large genome-wide association studies. Finally, we explored the causal effect of PAI-1 on cardiovascular risk factors including metabolic and subclinical atherosclerosis measures. In the systematic meta-analysis, the highest quantile of blood PAI-1 level was associated with higher CHD risk comparing with the lowest quantile (odds ratio=2.17; 95% CI: 1.53, 3.07) in an age- and sex-adjusted model. The effect size was reduced in studies using a multivariable-adjusted model (odds ratio=1.46; 95% CI: 1.13, 1.88). The Mendelian randomization analyses suggested a causal effect of increased PAI-1 level on CHD risk (odds ratio=1.22 per unit increase of log-transformed PAI-1; 95% CI: 1.01, 1.47). In addition, we also detected a causal effect of PAI-1 on elevating blood glucose and high-density lipoprotein cholesterol. Conclusions--Our study indicates a causal effect of elevated PAI-1 level on CHD risk, which may be mediated by glucose dysfunction

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function.

Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, 19 associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation of glucose metabolism. Chromatin state mapping and DNase I hypersensitivity analyses across adult tissues demonstrate preferential mapping of associated variants to regulatory regions in kidney but not extra-renal tissues. These findings suggest that genetic determinants of eGFR are mediated largely through direct effects within the kidney and highlight important cell types and biological pathways