Menopause age, reproductive span and hormone therapy duration predict the volume of medial temporal lobe brain structures in postmenopausal women

Medial temporal lobe (MTL) atrophy is correlated with risk and severity of Alzheimer disease (AD) pathology and cognitive decline. Increasing evidence suggest that oestrogens affect the aging of MTL structures. Here we investigate the relationship between reproductive hormone exposure, polygenic scores for AD risk and oestradiol concentration, MTL anatomy and cognitive performance in postmenopausal women. To this end, we used data from 10,924 female participants in the UK Biobank from whom brain MRI and genetic data were available. We fitted linear regression models to test whether the volume of structures comprising the MTL were predicted by a) timing related to menopause, b) the use and timing of hormone replacement therapy (HRT) and c) polygenic scores for AD risk and oestradiol concentration. Results showed that longer use of HRT was associated with larger parahippocampal volumes (2.53 mm3/year, p = 0.042). A later age of natural menopause, and a longer reproductive span, was associated with larger hippocampal (6.08 and 5.72 mm3/year, p = 0.0006 and 0.0005), parahippocampal (4.17 mm3 and 4.19 mm3/year, p = 0.00006 and 0.00001), amygdala (2.10 and 2.22 mm3/year, p = 0.028 and 0.01) and perirhinal cortical (2.56 and 2.95 mm3/year, p = 0.028 and 0.008) volumes. Superior prospective memory performance was associated with later age at natural menopause, and a longer reproductive span (ß = 0.05 and 0.05 respectively, p = 0.019 and 0.019). Polygenic scores for AD risk and for oestradiol concentration were not associated with MTL volume and did not interact with menopause-related factors to affect MTL structure. Our results suggest that HRT use did not have any detrimental effects on cognition or brain structure, whilst greater exposure to reproductive hormones across time is associated both with slightly larger volumes of specific MTL structures and marginally superior memory performance, independent of genetic risk for AD and genetic predisposition for higher oestradiol levels. However, the clinical utility of maintenance of oestrogens post-menopause for brain health and protection against cognitive decline is curtailed by the small effect sizes observed.

Insights into the genetics of menopausal vasomotor symptoms: genome-wide analyses of routinely-collected primary care health records

BACKGROUND

Vasomotor symptoms (VMS) can often significantly impact women's quality of life at menopause. In vivo studies have shown that increased neurokinin B (NKB) / neurokinin 3 receptor (NK3R) signalling contributes to VMS, with previous genetic studies implicating the TACR3 gene locus that encodes NK3R. Large-scale genomic analyses offer the possibility of biological insights but few such studies have collected data on VMS, while proxy phenotypes such as hormone replacement therapy (HRT) use are likely to be affected by changes in clinical practice. We investigated the genetic basis of VMS by analysing routinely-collected health records.

METHODS

We performed a GWAS of VMS derived from linked primary-care records and cross-sectional self-reported HRT use in up to 153,152 women from UK Biobank, a population-based cohort. In a subset of this cohort (n = 39,356), we analysed exome-sequencing data to test the association with VMS of rare deleterious genetic variants. Finally, we used Mendelian randomisation analysis to investigate the reasons for HRT use over time.

RESULTS

Our GWAS of health-records derived VMS identified a genetic signal near TACR3 associated with a lower risk of VMS (OR=0.76 (95% CI 0.72,0.80) per A allele, P=3.7x10-27), which was consistent with previous studies, validating this approach. Conditional analyses demonstrated independence of genetic signals for puberty timing and VMS at the TACR3 locus, including a rare variant predicted to reduce functional NK3R levels that was associated with later menarche (P = 5 × 10-9) but showed no association with VMS (P = 0.6). Younger menopause age was causally-associated with greater HRT use before 2002 but not after.

CONCLUSIONS

We provide support for TACR3 in the genetic basis of VMS but unexpectedly find that rare genomic variants predicted to lower NK3R levels did not modify VMS, despite the proven efficacy of NK3R antagonists. Using genomics we demonstrate changes in genetic associations with HRT use over time, arising from a change in clinical practice since the early 2000s, which is likely to reflect a switch from preventing post-menopausal complications in women with earlier menopause to primarily treating VMS. Our study demonstrates that integrating routinely-collected primary care health records and genomic data offers great potential for exploring the genetic basis of symptoms.

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency

Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI.

Genetically proxied therapeutic prolyl-hydroxylase inhibition and cardiovascular risk

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

Detection and characterization of male sex chromosome abnormalities in the UK Biobank study

PURPOSE

The study aimed to systematically ascertain male sex chromosome abnormalities, 47,XXY (Klinefelter syndrome [KS]) and 47,XYY, and characterize their risks of adverse health outcomes.

METHODS

We analyzed genotyping array or exome sequence data in 207,067 men of European ancestry aged 40 to 70 years from the UK Biobank and related these to extensive routine health record data.

RESULTS

Only 49 of 213 (23%) of men whom we identified with KS and only 1 of 143 (0.7%) with 47,XYY had a diagnosis of abnormal karyotype on their medical records or self-report. We observed expected associations for KS with reproductive dysfunction (late puberty: risk ratio [RR] = 2.7; childlessness: RR = 4.2; testosterone concentration: RR = -3.8 nmol/L, all P < 2 × 10-8), whereas XYY men appeared to have normal reproductive function. Despite this difference, we identified several higher disease risks shared across both KS and 47,XYY, including type 2 diabetes (RR = 3.0 and 2.6, respectively), venous thrombosis (RR = 6.4 and 7.4, respectively), pulmonary embolism (RR = 3.3 and 3.7, respectively), and chronic obstructive pulmonary disease (RR = 4.4 and 4.6, respectively) (all P < 7 × 10-6).

CONCLUSION

KS and 47,XYY were mostly unrecognized but conferred substantially higher risks for metabolic, vascular, and respiratory diseases, which were only partially explained by higher levels of body mass index, deprivation, and smoking.

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 (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.

Identifying molecular mediators of the relationship between body mass index and endometrial cancer risk: a Mendelian randomization analysis

BACKGROUND

Endometrial cancer is the most common gynaecological cancer in high-income countries. Elevated body mass index (BMI) is an established modifiable risk factor for this condition and is estimated to confer a larger effect on endometrial cancer risk than any other cancer site. However, the molecular mechanisms underpinning this association remain unclear. We used Mendelian randomization (MR) to evaluate the causal role of 14 molecular risk factors (hormonal, metabolic and inflammatory markers) in endometrial cancer risk. We then evaluated and quantified the potential mediating role of these molecular traits in the relationship between BMI and endometrial cancer using multivariable MR.

METHODS

Genetic instruments to proxy 14 molecular risk factors and BMI were constructed by identifying single-nucleotide polymorphisms (SNPs) reliably associated (P < 5.0 × 10-8) with each respective risk factor in previous genome-wide association studies (GWAS). Summary statistics for the association of these SNPs with overall and subtype-specific endometrial cancer risk (12,906 cases and 108,979 controls) were obtained from a GWAS meta-analysis of the Endometrial Cancer Association Consortium (ECAC), Epidemiology of Endometrial Cancer Consortium (E2C2) and UK Biobank. SNPs were combined into multi-allelic models and odds ratios (ORs) and 95% confidence intervals (95% CIs) were generated using inverse-variance weighted random-effects models. The mediating roles of the molecular risk factors in the relationship between BMI and endometrial cancer were then estimated using multivariable MR.

RESULTS

In MR analyses, there was strong evidence that BMI (OR per standard deviation (SD) increase 1.88, 95% CI 1.69 to 2.09, P = 3.87 × 10-31), total testosterone (OR per inverse-normal transformed nmol/L increase 1.64, 95% CI 1.43 to 1.88, P = 1.71 × 10-12), bioavailable testosterone (OR per natural log transformed nmol/L increase: 1.46, 95% CI 1.29 to 1.65, P = 3.48 × 10-9), fasting insulin (OR per natural log transformed pmol/L increase: 3.93, 95% CI 2.29 to 6.74, P = 7.18 × 10-7) and sex hormone-binding globulin (SHBG, OR per inverse-normal transformed nmol/L increase 0.71, 95% CI 0.59 to 0.85, P = 2.07 × 10-4) had a causal effect on endometrial cancer risk. Additionally, there was suggestive evidence that total serum cholesterol (OR per mg/dL increase 0.90, 95% CI 0.81 to 1.00, P = 4.01 × 10-2) had an effect on endometrial cancer risk. In mediation analysis, we found evidence for a mediating role of fasting insulin (19% total effect mediated, 95% CI 5 to 34%, P = 9.17 × 10-3), bioavailable testosterone (15% mediated, 95% CI 10 to 20%, P = 1.43 × 10-8) and SHBG (7% mediated, 95% CI 1 to 12%, P = 1.81 × 10-2) in the relationship between BMI and endometrial cancer risk.

CONCLUSIONS

Our comprehensive MR analysis provides insight into potential causal mechanisms linking BMI with endometrial cancer risk and suggests targeting of insulinemic and hormonal traits as a potential strategy for the prevention of endometrial cancer.

Do sex hormones confound or mediate the effect of chronotype on breast and prostate cancer? A Mendelian randomization study

Morning-preference chronotype has been found to be protective against breast and prostate cancer. Sex hormones have been implicated in relation to chronotype and the development of both cancers. This study aimed to assess whether sex hormones confound or mediate the effect of chronotype on breast and prostate cancer using a Mendelian Randomization (MR) framework. Genetic variants associated with chronotype and sex hormones (total testosterone, bioavailable testosterone, sex hormone binding globulin, and oestradiol) (p<5×10-8) were obtained from published genome-wide association studies (n≤244,207 females and n≤205,527 males). These variants were used to investigate causal relationships with breast (nCases/nControls = 133,384/113,789) and prostate (nCases/nControls = 79,148/61,106) cancer using univariable, bidirectional and multivariable MR. In females, we found evidence for: I) Reduced risk of breast cancer per category increase in morning-preference (OR = 0.93, 95% CI:0. 88, 1.00); II) Increased risk of breast cancer per SD increase in bioavailable testosterone (OR = 1.10, 95% CI: 1.01, 1.19) and total testosterone (OR = 1.15, 95% CI:1.07, 1.23); III) Bidirectional effects between morning-preference and both bioavailable and total testosterone (e.g. mean SD difference in bioavailable testosterone = -0.08, 95% CI:-0.12, -0.05 per category increase in morning-preference vs difference in morning-preference category = -0.04, 95% CI: -0.08, 0.00 per SD increase in bioavailable testosterone). In males, we found evidence for: I) Reduced risk of prostate cancer per category increase in morning-preference (OR = 0.90, 95% CI: 0.83, 0.97) and II) Increased risk of prostate cancer per SD increase in bioavailable testosterone (OR = 1.22, 95% CI: 1.08, 1.37). No bidirectional effects were found between morning-preference and testosterone in males. While testosterone levels were causally implicated with both chronotype and cancer, there was inconsistent evidence for testosterone as a mediator of the relationship. The protective effect of morning-preference on both breast and prostate cancer is clinically interesting, although it may be difficult to effectively modify chronotype. Further studies are needed to investigate other potentially modifiable intermediates.

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.

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

Anti-Müllerian hormone (AMH) is expressed by antral stage ovarian follicles in women. Consequently, circulating AMH levels are detectable until menopause. Variation in age-specific AMH levels has been associated with breast cancer and polycystic ovary syndrome (PCOS), amongst other diseases. Identification of genetic variants underlying variation in AMH levels could provide clues about the physiological mechanisms that explain these AMH-disease associations. To date, only one variant in MCM8 has been identified to be associated with circulating AMH levels in women. We aimed to identify additional variants for AMH through a GWAS meta-analysis including data from 7049 premenopausal women of European ancestry, which more than doubles the sample size of the largest previous GWAS. 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 SNPs for AMH levels and for age at menopause (r g = 0.82, FDR=0.003). Exploratory Mendelian randomization analyses did not support a causal effect of AMH on breast cancer or PCOS risk, but should be interpreted with caution as they may be underpowered and the validity of genetic instruments could not be extensively explored. In conclusion, we identified a variant in the AMH gene and three other loci that may affect circulating AMH levels in women.

Large Copy-Number Variants in UK Biobank Caused by Clonal Hematopoiesis May Confound Penetrance Estimates

Large copy-number variants (CNVs) are strongly associated with both developmental delay and cancer, but the type of disease depends strongly on when and where the mutation occurred, i.e., germline versus somatic. We used microarray data from UK Biobank to investigate the prevalence and penetrance of large autosomal CNVs and chromosomal aneuploidies using a standard CNV detection algorithm not designed for detecting mosaic variants. We found 160 individuals that carry >10 Mb copy number changes, including 56 with whole chromosome aneuploidies. Nineteen (12%) individuals had a diagnosis of Down syndrome or other developmental disorder, while 84 (52.5%) individuals had a diagnosis of hematological malignancies or chronic myeloproliferative disorders. Notably, there was no evidence of mosaicism in the blood for many of these large CNVs, so they could easily be mistaken for germline alleles even when caused by somatic mutations. We therefore suggest that somatic mutations associated with blood cancers may result in false estimates of rare variant penetrance from population biobanks.

Using human genetics to understand the disease impacts of testosterone in men and women

Testosterone supplementation is commonly used for its effects on sexual function, bone health and body composition, yet its effects on disease outcomes are unknown. To better understand this, we identified genetic determinants of testosterone levels and related sex hormone traits in 425,097 UK Biobank study participants. Using 2,571 genome-wide significant associations, we demonstrate that the genetic determinants of testosterone levels are substantially different between sexes and that genetically higher testosterone is harmful for metabolic diseases in women but beneficial in men. For example, a genetically determined 1 s.d. higher testosterone increases the risks of type 2 diabetes (odds ratio (OR) = 1.37 (95% confidence interval (95% CI): 1.22-1.53)) and polycystic ovary syndrome (OR = 1.51 (95% CI: 1.33-1.72)) in women, but reduces type 2 diabetes risk in men (OR = 0.86 (95% CI: 0.76-0.98)). We also show adverse effects of higher testosterone on breast and endometrial cancers in women and prostate cancer in men. Our findings provide insights into the disease impacts of testosterone and highlight the importance of sex-specific genetic analyses.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Using genetics to understand the causal influence of higher BMI on depression

BACKGROUND

Depression is more common in obese than non-obese individuals, especially in women, but the causal relationship between obesity and depression is complex and uncertain. Previous studies have used genetic variants associated with BMI to provide evidence that higher body mass index (BMI) causes depression, but have not tested whether this relationship is driven by the metabolic consequences of BMI nor for differences between men and women.

METHODS

We performed a Mendelian randomization study using 48 791 individuals with depression and 291 995 controls in the UK Biobank, to test for causal effects of higher BMI on depression (defined using self-report and Hospital Episode data). We used two genetic instruments, both representing higher BMI, but one with and one without its adverse metabolic consequences, in an attempt to 'uncouple' the psychological component of obesity from the metabolic consequences. We further tested causal relationships in men and women separately, and using subsets of BMI variants from known physiological pathways.

RESULTS

Higher BMI was strongly associated with higher odds of depression, especially in women. Mendelian randomization provided evidence that higher BMI partly causes depression. Using a 73-variant BMI genetic risk score, a genetically determined one standard deviation (1 SD) higher BMI (4.9 kg/m2) was associated with higher odds of depression in all individuals [odds ratio (OR): 1.18, 95% confidence interval (CI): 1.09, 1.28, P = 0.00007) and women only (OR: 1.24, 95% CI: 1.11, 1.39, P = 0.0001). Meta-analysis with 45 591 depression cases and 97 647 controls from the Psychiatric Genomics Consortium (PGC) strengthened the statistical confidence of the findings in all individuals. Similar effect size estimates were obtained using different Mendelian randomization methods, although not all reached P < 0.05. Using a metabolically favourable adiposity genetic risk score, and meta-analysing data from the UK biobank and PGC, a genetically determined 1 SD higher BMI (4.9 kg/m2) was associated with higher odds of depression in all individuals (OR: 1.26, 95% CI: 1.06, 1.50], P = 0.010), but with weaker statistical confidence.

CONCLUSIONS

Higher BMI, with and without its adverse metabolic consequences, is likely to have a causal role in determining the likelihood of an individual developing depression.

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Genome-wide association study of anti-Müllerian hormone levels in pre-menopausal women of late reproductive age and relationship with genetic determinants of reproductive lifespan

Anti-Müllerian hormone (AMH) is required for sexual differentiation in the fetus, and in adult females AMH is produced by growing ovarian follicles. Consequently, AMH levels are correlated with ovarian reserve, declining towards menopause when the oocyte pool is exhausted. A previous genome-wide association study identified three genetic variants in and around the AMH gene that explained 25% of variation in AMH levels in adolescent males but did not identify any genetic associations reaching genome-wide significance in adolescent females. To explore the role of genetic variation in determining AMH levels in women of late reproductive age, we carried out a genome-wide meta-analysis in 3344 pre-menopausal women from five cohorts (median age 44-48 years at blood draw). A single genetic variant, rs16991615, previously associated with age at menopause, reached genome-wide significance at P = 3.48 × 10-10, with a per allele difference in age-adjusted inverse normal AMH of 0.26 standard deviations (SD) (95% confidence interval (CI) [0.18,0.34]). We investigated whether genetic determinants of female reproductive lifespan were more generally associated with pre-menopausal AMH levels. Genetically-predicted age at menarche had no robust association but genetically-predicted age at menopause was associated with lower AMH levels by 0.18 SD (95% CI [0.14,0.21]) in age-adjusted inverse normal AMH per one-year earlier age at menopause. Our findings provide genetic support for the well-established use of AMH as a marker of ovarian reserve.

Genetic studies of accelerometer-based sleep measures yield new insights into human sleep behaviour

Sleep is an essential human function but its regulation is poorly understood. Using accelerometer data from 85,670 UK Biobank participants, we perform a genome-wide association study of 8 derived sleep traits representing sleep quality, quantity and timing, and validate our findings in 5,819 individuals. We identify 47 genetic associations at P < 5 × 10-8, of which 20 reach a stricter threshold of P < 8 × 10-10. These include 26 novel associations with measures of sleep quality and 10 with nocturnal sleep duration. The majority of identified variants associate with a single sleep trait, except for variants previously associated with restless legs syndrome. For sleep duration we identify a missense variant (p.Tyr727Cys) in PDE11A as the likely causal variant. As a group, sleep quality loci are enriched for serotonin processing genes. Although accelerometer-derived measures of sleep are imperfect and may be affected by restless legs syndrome, these findings provide new biological insights into sleep compared to previous efforts based on self-report sleep measures.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms

Being a morning person is a behavioural indicator of a person’s underlying circadian rhythm. Using genome-wide data from 697,828 UK Biobank and 23andMe participants we increase the number of genetic loci associated with being a morning person from 24 to 351. Using data from 85,760 individuals with activity-monitor derived measures of sleep timing we find that the chronotype loci associate with sleep timing: the mean sleep timing of the 5% of individuals carrying the most morningness alleles is 25 min earlier than the 5% carrying the fewest. The loci are enriched for genes involved in circadian regulation, cAMP, glutamate and insulin signalling pathways, and those expressed in the retina, hindbrain, hypothalamus, and pituitary. Using Mendelian Randomisation, we show that being a morning person is causally associated with better mental health but does not affect BMI or risk of Type 2 diabetes. This study offers insights into circadian biology and its links to disease in humans.

GWAS have previously found 24 genomic loci associated with chronotype, an individual’s preference for early or late sleep timing. Here, the authors identify 327 additional loci in a sample of 697,828 individuals and further explore the relationships of chronotype with metabolic and psychiatric diseases.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother–child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10⁻⁸. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

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⁻¹⁴), 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.

Genome-Wide and Abdominal MRI Data Provide Evidence That a Genetically Determined Favorable Adiposity Phenotype Is Characterized by Lower Ectopic Liver Fat and Lower Risk of Type 2 Diabetes, Heart Disease, and Hypertension

Recent genetic studies have identified alleles associated with opposite effects on adiposity and risk of type 2 diabetes. We aimed to identify more of these variants and test the hypothesis that such favorable adiposity alleles are associated with higher subcutaneous fat and lower ectopic fat. We combined MRI data with genome-wide association studies of body fat percentage (%) and metabolic traits. We report 14 alleles, including 7 newly characterized alleles, associated with higher adiposity but a favorable metabolic profile. Consistent with previous studies, individuals carrying more favorable adiposity alleles had higher body fat % and higher BMI but lower risk of type 2 diabetes, heart disease, and hypertension. These individuals also had higher subcutaneous fat but lower liver fat and a lower visceral-to-subcutaneous adipose tissue ratio. Individual alleles associated with higher body fat % but lower liver fat and lower risk of type 2 diabetes included those in PPARG, GRB14, and IRS1, whereas the allele in ANKRD55 was paradoxically associated with higher visceral fat but lower risk of type 2 diabetes. Most identified favorable adiposity alleles are associated with higher subcutaneous and lower liver fat, a mechanism consistent with the beneficial effects of storing excess triglycerides in metabolically low-risk depots.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

A Common Allele in FGF21 Associated with Sugar Intake Is Associated with Body Shape, Lower Total Body-Fat Percentage, and Higher Blood Pressure

Fibroblast growth factor 21 (FGF21) is a hormone that has insulin-sensitizing properties. Some trials of FGF21 analogs show weight loss and lipid-lowering effects. Recent studies have shown that a common allele in the FGF21 gene alters the balance of macronutrients consumed, but there was little evidence of an effect on metabolic traits. We studied a common FGF21 allele (A:rs838133) in 451,099 people from the UK Biobank study, aiming to use the human allele to inform potential adverse and beneficial effects of targeting FGF21. We replicated the association between the A allele and higher percentage carbohydrate intake. We then showed that this allele is more strongly associated with higher blood pressure and waist-hip ratio, despite an association with lower total body-fat percentage, than it is with BMI or type 2 diabetes. These human phenotypes of variation in the FGF21 gene will inform research into FGF21's mechanisms and therapeutic potential.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

Red blood cell distribution width: Genetic evidence for aging pathways in 116,666 volunteers

INTRODUCTION

Variability in red blood cell volumes (distribution width, RDW) increases with age and is strongly predictive of mortality, incident coronary heart disease and cancer. We investigated inherited genetic variation associated with RDW in 116,666 UK Biobank human volunteers.

RESULTS

A large proportion RDW is explained by genetic variants (29%), especially in the older group (60+ year olds, 33.8%, <50 year olds, 28.4%). RDW was associated with 194 independent genetic signals; 71 are known for conditions including autoimmune disease, certain cancers, BMI, Alzheimer's disease, longevity, age at menopause, bone density, myositis, Parkinson's disease, and age-related macular degeneration. Exclusion of anemic participants did not affect the overall findings. Pathways analysis showed enrichment for telomere maintenance, ribosomal RNA, and apoptosis. The majority of RDW-associated signals were intronic (119 of 194), including SNP rs6602909 located in an intron of oncogene GAS6, an eQTL in whole blood.

CONCLUSIONS

Although increased RDW is predictive of cardiovascular outcomes, this was not explained by known CVD or related lipid genetic risks, and a RDW genetic score was not predictive of incident disease. The predictive value of RDW for a range of negative health outcomes may in part be due to variants influencing fundamental pathways of aging.

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to ∼370,000 women, we identify 389 independent signals (P < 5 × 10-8) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain ∼7.4% of the population variance in age at menarche, corresponding to ∼25% of the estimated heritability. We implicate ∼250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility.

Quantifying the extent to which index event biases influence large genetic association studies

As genetic association studies increase in size to 100 000s of individuals, subtle biases may influence conclusions. One possible bias is 'index event bias' (IEB) that appears due to the stratification by, or enrichment for, disease status when testing associations between genetic variants and a disease-associated trait. We aimed to test the extent to which IEB influences some known trait associations in a range of study designs and provide a statistical framework for assessing future associations. Analyzing data from 113 203 non-diabetic UK Biobank participants, we observed three (near TCF7L2, CDKN2AB and CDKAL1) overestimated (body mass index (BMI) decreasing) and one (near MTNR1B) underestimated (BMI increasing) associations among 11 type 2 diabetes risk alleles (at P  <  0.05). IEB became even stronger when we tested a type 2 diabetes genetic risk score composed of these 11 variants (-0.010 standard deviations BMI per allele, P  =  5 × 10- 4), which was confirmed in four additional independent studies. Similar results emerged when examining the effect of blood pressure increasing alleles on BMI in normotensive UK Biobank samples. Furthermore, we demonstrated that, under realistic scenarios, common disease alleles would become associated at P <  5 × 10- 8 with disease-related traits through IEB alone, if disease prevalence in the sample differs appreciably from the background population prevalence. For example, some hypertension and type 2 diabetes alleles will be associated with BMI in sample sizes of  >500 000 if the prevalence of those diseases differs by >10% from the background population. In conclusion, IEB may result in false positive or negative genetic associations in very large studies stratified or strongly enriched for/against disease cases.

Gene-obesogenic environment interactions in the UK Biobank study

Background

Previous studies have suggested that modern obesogenic environments accentuate the genetic risk of obesity. However, these studies have proven controversial as to which, if any, measures of the environment accentuate genetic susceptibility to high body mass index (BMI).

Methods

We used up to 120 000 adults from the UK Biobank study to test the hypothesis that high-risk obesogenic environments and behaviours accentuate genetic susceptibility to obesity. We used BMI as the outcome and a 69-variant genetic risk score (GRS) for obesity and 12 measures of the obesogenic environment as exposures. These measures included Townsend deprivation index (TDI) as a measure of socio-economic position, TV watching, a 'Westernized' diet and physical activity. We performed several negative control tests, including randomly selecting groups of different average BMIs, using a simulated environment and including sun-protection use as an environment.

Results

We found gene-environment interactions with TDI (Pinteraction = 3 × 10 -10 ), self-reported TV watching (Pinteraction = 7 × 10 -5 ) and self-reported physical activity (Pinteraction = 5 × 10 -6 ). Within the group of 50% living in the most relatively deprived situations, carrying 10 additional BMI-raising alleles was associated with approximately 3.8 kg extra weight in someone 1.73 m tall. In contrast, within the group of 50% living in the least deprivation, carrying 10 additional BMI-raising alleles was associated with approximately 2.9 kg extra weight. The interactions were weaker, but present, with the negative controls, including sun-protection use, indicating that residual confounding is likely.

Conclusions

Our findings suggest that the obesogenic environment accentuates the risk of obesity in genetically susceptible adults. Of the factors we tested, relative social deprivation best captures the aspects of the obesogenic environment responsible.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

Genome-wide associations for birth weight and correlations with adult disease

Birth weight (BW) has been shown to be influenced by both fetal and maternal factors and in observational studies is reproducibly associated with future risk of adult metabolic diseases including type 2 diabetes (T2D) and cardiovascular disease. These life-course associations have often been attributed to the impact of an adverse early life environment. Here, we performed a multi-ancestry genome-wide association study (GWAS) meta-analysis of BW in 153,781 individuals, identifying 60 loci where fetal genotype was associated with BW (P < 5 × 10-8). Overall, approximately 15% of variance in BW was captured by assays of fetal genetic variation. Using genetic association alone, we found strong inverse genetic correlations between BW and systolic blood pressure (Rg = -0.22, P = 5.5 × 10-13), T2D (Rg = -0.27, P = 1.1 × 10-6) and coronary artery disease (Rg = -0.30, P = 6.5 × 10-9). In addition, using large -cohort datasets, we demonstrated that genetic factors were the major contributor to the negative covariance between BW and future cardiometabolic risk. Pathway analyses indicated that the protein products of genes within BW-associated regions were enriched for diverse processes including insulin signalling, glucose homeostasis, glycogen biosynthesis and chromatin remodelling. There was also enrichment of associations with BW in known imprinted regions (P = 1.9 × 10-4). We demonstrate that life-course associations between early growth phenotypes and adult cardiometabolic disease are in part the result of shared genetic effects and identify some of the pathways through which these causal genetic effects are mediated.

Length of FMR1 repeat alleles within the normal range does not substantially affect the risk of early menopause

STUDY QUESTION

Is the length of FMR1 repeat alleles within the normal range associated with the risk of early menopause?

SUMMARY ANSWER

The length of repeat alleles within the normal range does not substantially affect risk of early menopause.

WHAT IS KNOWN ALREADY

There is a strong, well-established relationship between length of premutation FMR1 alleles and age at menopause, suggesting that this relationship could continue into the normal range. Within the normal range, there is conflicting evidence; differences in ovarian reserve have been identified with FMR1 repeat allele length, but a recent population-based study did not find any association with age at menopause as a quantitative trait.

STUDY DESIGN, SIZE, DURATION

We analysed cross-sectional baseline survey data collected at recruitment from 2004 to 2010 from a population-based, prospective epidemiological cohort study of >110 000 women to investigate whether repeat allele length was associated with early menopause.

PARTICIPANTS/MATERIALS, SETTING, METHOD

We included 4333 women from the Breakthrough Generations Study (BGS), of whom 2118 were early menopause cases (menopause under 46 years) and 2215 were controls. We analysed the relationship between length of FMR1 alleles and early menopause using logistic regression with allele length as continuous and categorical variables. We also conducted analyses with the outcome age at menopause as a quantitative trait as well as appropriate sensitivity and exploratory analyses.

MAIN RESULTS AND THE ROLE OF CHANCE

There was no association of the shorter or longer FMR1 allele or their combined genotype with the clinically relevant end point of early menopause in our main analysis. Likewise, there were no associations with age at menopause as a quantitative trait in our secondary analysis.

LIMITATIONS, REASONS FOR CAUTION

Women with homozygous alleles in the normal range may have undetected FMR1 premutation alleles, although there was no evidence to suggest this. We estimate minor dilution of risk of early menopause from the likely inclusion of some women with menopause at over 45 years in the early menopause cases due to age-rounding bias in self-reports.

WIDER IMPLICATIONS OF THE FINDINGS

There is no robust evidence in this large study that variation within the normal range of FMR1 repeat alleles influences timing of menopause in the general population, which contradicts findings from some earlier, mainly smaller studies. The FMR1 CGG repeat polymorphism in the normal range is unlikely to contribute to genetic susceptibility to early menopause.

STUDY FUNDING/COMPETING INTERESTS

We thank Breast Cancer Now and The Institute of Cancer Research for funding the BGS. The Institute of Cancer Research acknowledges NHS funding to the NIHR Biomedical Research Centre. The study was funded by the Wellcome Trust (grant number 085943). There are no competing interests.

TRIAL REGISTRATION NUMBER

Not applicable.

Lessons from Genome-Wide Association Studies in Reproductive Medicine: Menopause

In recent years, common genetic variants have been identified by genome-wide association studies (GWASs) that have led to the detection of 44 genetic loci associated with approximately 6% of common variation in age at natural menopause. In the latest GWAS, doubling the sample size to approximately 70,000 women more than doubled the number of signals identified, from 17 to 56. In addition, low-frequency coding variants (< 5% minor allele frequency), with relatively large effect sizes, have been identified in two genes, by analyzing genome-wide exome data. GWAS has been very successful in identifying novel biological pathways involved in reproductive aging. Approximately two-thirds of the loci reported so far include genes involved in DNA damage response (DDR), highlighting the importance of this pathway in determining oocyte reserve. In addition, GWAS demonstrates that the hypothalamic-pituitary axis is involved in menopause timing as well as puberty timing, showing the first genetic link between timing of the start and end of reproductive life. Genetic variants have been used to explore the causal relationships between menopause timing and breast cancer. These studies demonstrate that for a 1 year increase in menopause age, there is a 6% increase in breast cancer risk, a value approximately double the estimate from epidemiological studies. Prolonged exposure to estrogen during reproductive life is the likely mechanism, rather than a direct effect of DDR variants on cancer risk. Further work is needed to determine the mechanism for the effect of each variant identified by GWAS and more variants will undoubtedly be discovered as sample sizes increase, denser single nucleotide polymorphism arrays and reference genomes are used, and populations from diverse ethnic groups are studied.

Genome-Wide Association Analyses in 128,266 Individuals Identifies New Morningness and Sleep Duration Loci

Disrupted circadian rhythms and reduced sleep duration are associated with several human diseases, particularly obesity and type 2 diabetes, but until recently, little was known about the genetic factors influencing these heritable traits. We performed genome-wide association studies of self-reported chronotype (morning/evening person) and self-reported sleep duration in 128,266 white British individuals from the UK Biobank study. Sixteen variants were associated with chronotype (P<5x10^-8), including variants near the known circadian rhythm genes RGS16 (1.21 odds of morningness, 95% CI [1.15, 1.27], P = 3x10^-12) and PER2 (1.09 odds of morningness, 95% CI [1.06, 1.12], P = 4x10^-10). The PER2 signal has previously been associated with iris function. We sought replication using self-reported data from 89,283 23andMe participants; thirteen of the chronotype signals remained associated at P<5x10^-8 on meta-analysis and eleven of these reached P<0.05 in the same direction in the 23andMe study. We also replicated 9 additional variants identified when the 23andMe study was used as a discovery GWAS of chronotype (all P<0.05 and meta-analysis P<5x10^-8). For sleep duration, we replicated one known signal in PAX8 (2.6 minutes per allele, 95% CI [1.9, 3.2], P = 5.7x10^-16) and identified and replicated two novel associations at VRK2 (2.0 minutes per allele, 95% CI [1.3, 2.7], P = 1.2x10^-9; and 1.6 minutes per allele, 95% CI [1.1, 2.2], P = 7.6x10^-9). Although we found genetic correlation between chronotype and BMI (rG = 0.056, P = 0.05); undersleeping and BMI (rG = 0.147, P = 1x10^-5) and oversleeping and BMI (rG = 0.097, P = 0.04), Mendelian Randomisation analyses, with limited power, provided no consistent evidence of causal associations between BMI or type 2 diabetes and chronotype or sleep duration. Our study brings the total number of loci associated with chronotype to 22 and with sleep duration to three, and provides new insights into the biology of sleep and circadian rhythms in humans.

Numerous studies have identified links between too little or too much sleep and circadian misalignment with metabolic disorders such as obesity and type 2 diabetes. However, cause-and-effect is not easily determined, because of multiple confounding factors affecting both sleep patterns and disease risk. Using the first release of the UK Biobank study, which combines detailed measurements and questionnaire data with genetic data, we investigate the genetics of two self-report sleep measures, chronotype and average sleep duration, in 128,266 white British individuals. We replicate previous genetic associations and identify seven and two novel genetic variants influencing chronotype and sleep duration, respectively. Associated variants are located near genes implicated in circadian rhythm regulation (RGS16, PER2), near a serotonin receptor gene (HTR6) and another gene (INADL) encoding a protein thought to be important in photosensitive retinal cells, cells known to communicate with the brain’s primary circadian pacemaker. Using the genetic risk factors, we estimate the unconfounded causal associations of BMI and type 2 diabetes on sleep patterns (and vice versa) through Mendelian Randomisation. However, we find no evidence for causal associations in either direction. The full UK Biobank release of 500,000 individuals will boost our power to detect causal associations.

Genetic Evidence for a Link Between Favorable Adiposity and Lower Risk of Type 2 Diabetes, Hypertension, and Heart Disease

Recent genetic studies have identified some alleles that are associated with higher BMI but lower risk of type 2 diabetes, hypertension, and heart disease. These "favorable adiposity" alleles are collectively associated with lower insulin levels and higher subcutaneous-to-visceral adipose tissue ratio and may protect from disease through higher adipose storage capacity. We aimed to use data from 164,609 individuals from the UK Biobank and five other studies to replicate associations between a genetic score of 11 favorable adiposity variants and adiposity and risk of disease, to test for interactions between BMI and favorable adiposity genetics, and to test effects separately in men and women. In the UK Biobank, the 50% of individuals carrying the most favorable adiposity alleles had higher BMIs (0.120 kg/m(2) [95% CI 0.066, 0.174]; P = 1E-5) and higher body fat percentage (0.301% [0.230, 0.372]; P = 1E-16) compared with the 50% of individuals carrying the fewest alleles. For a given BMI, the 50% of individuals carrying the most favorable adiposity alleles were at lower risk of type 2 diabetes (odds ratio [OR] 0.837 [0.784, 0.894]; P = 1E-7), hypertension (OR 0.935 [0.911, 0.958]; P = 1E-7), and heart disease (OR 0.921 [0.872, 0.973]; P = 0.003) and had lower blood pressure (systolic -0.859 mmHg [-1.099, -0.618]; P = 3E-12 and diastolic -0.394 mmHg [-0.534, -0.254]; P = 4E-8). In women, these associations could be explained by the observation that the alleles associated with higher BMI but lower risk of disease were also associated with a favorable body fat distribution, with a lower waist-to-hip ratio (-0.004 cm [95% CI -0.005, -0.003] 50% vs. 50%; P = 3E-14), but in men, the favorable adiposity alleles were associated with higher waist circumference (0.454 cm [0.267, 0.641] 50% vs. 50%; P = 2E-6) and higher waist-to-hip ratio (0.0013 [0.0003, 0.0024] 50% vs. 50%; P = 0.01). Results were strengthened when a meta-analysis with five additional studies was conducted. There was no evidence of interaction between a genetic score consisting of known BMI variants and the favorable adiposity genetic score. In conclusion, different molecular mechanisms that lead to higher body fat percentage (with greater subcutaneous storage capacity) can have different impacts on cardiometabolic disease risk. Although higher BMI is associated with higher risk of diseases, better fat storage capacity could reduce the risk.

Variants in the FTO and CDKAL1 loci have recessive effects on risk of obesity and type 2 diabetes, respectively

AIMS/HYPOTHESIS

Genome-wide association (GWA) studies have identified hundreds of common genetic variants associated with obesity and type 2 diabetes. These studies have usually focused on additive association tests. Identifying deviations from additivity may provide new biological insights and explain some of the missing heritability for these diseases.

METHODS

We performed a GWA study using a dominance deviation model for BMI, obesity (29,925 cases) and type 2 diabetes (4,040 cases) in 120,286 individuals of British ancestry from the UK Biobank study. We also investigated whether single nucleotide polymorphisms previously shown to be associated with these traits showed any enrichment for departures from additivity.

RESULTS

Known obesity-associated variants in FTO showed strong evidence of deviation from additivity (p DOMDEV = 3 × 10(-5)) through a recessive effect of the allele associated with higher BMI. The average BMI of individuals carrying zero, one or two BMI-raising alleles was 27.27 (95% CI 27.22, 27.31) kg/m(2), 27.54 (95% CI 27.50, 27.58) kg/m(2) and 28.07 (95% CI 28.00, 28.14) kg/m(2), respectively. A similar effect was observed in 105,643 individuals from the GIANT Consortium (p DOMDEV = 0.003; meta-analysis p DOMDEV = 1 × 10(-7)). For type 2 diabetes, we detected a recessive effect (p DOMDEV = 5 × 10(-4)) at CDKAL1. Relative to homozygous non-risk allele carriers, homozygous risk allele carriers had an OR of 1.48 (95% CI 1.32, 1.65), while the heterozygous group had an OR of 1.06 (95% CI 0.99, 1.14), a result consistent with that of a previous study. We did not identify any novel associations at genome-wide significance.

CONCLUSIONS/INTERPRETATION

Although we found no evidence of widespread non-additive genetic effects contributing to obesity and type 2 diabetes risk, we did find robust examples of recessive effects at the FTO and CDKAL1 loci.

ACCESS TO RESEARCH MATERIALS

Summary statistics are available at www.t2diabetesgenes.org and by request (a.r.wood@exeter.ac.uk). All underlying data are available on application from the UK Biobank.

Height, body mass index, and socioeconomic status: mendelian randomisation study in UK Biobank

OBJECTIVE

To determine whether height and body mass index (BMI) have a causal role in five measures of socioeconomic status.

DESIGN

Mendelian randomisation study to test for causal effects of differences in stature and BMI on five measures of socioeconomic status. Mendelian randomisation exploits the fact that genotypes are randomly assigned at conception and thus not confounded by non-genetic factors.

SETTING

UK Biobank.

PARTICIPANTS

119,669 men and women of British ancestry, aged between 37 and 73 years.

MAIN OUTCOME MEASURES

Age completed full time education, degree level education, job class, annual household income, and Townsend deprivation index.

RESULTS

In the UK Biobank study, shorter stature and higher BMI were observationally associated with several measures of lower socioeconomic status. The associations between shorter stature and lower socioeconomic status tended to be stronger in men, and the associations between higher BMI and lower socioeconomic status tended to be stronger in women. For example, a 1 standard deviation (SD) higher BMI was associated with a £210 (€276; $300; 95% confidence interval £84 to £420; P=6 × 10(-3)) lower annual household income in men and a £1890 (£1680 to £2100; P=6 × 10(-15)) lower annual household income in women. Genetic analysis provided evidence that these associations were partly causal. A genetically determined 1 SD (6.3 cm) taller stature caused a 0.06 (0.02 to 0.09) year older age of completing full time education (P=0.01), a 1.12 (1.07 to 1.18) times higher odds of working in a skilled profession (P=6 × 10(-7)), and a £1130 (£680 to £1580) higher annual household income (P=4 × 10(-8)). Associations were stronger in men. A genetically determined 1 SD higher BMI (4.6 kg/m(2)) caused a £2940 (£1680 to £4200; P=1 × 10(-5)) lower annual household income and a 0.10 (0.04 to 0.16) SD (P=0.001) higher level of deprivation in women only.

CONCLUSIONS

These data support evidence that height and BMI play an important partial role in determining several aspects of a person's socioeconomic status, especially women's BMI for income and deprivation and men's height for education, income, and job class. These findings have important social and health implications, supporting evidence that overweight people, especially women, are at a disadvantage and that taller people, especially men, are at an advantage.

Human longevity is influenced by many genetic variants: evidence from 75,000 UK Biobank participants

Variation in human lifespan is 20 to 30% heritable in twins but few genetic variants have been identified. We undertook a Genome Wide Association Study (GWAS) using age at death of parents of middle-aged UK Biobank participants of European decent (n=75,244 with father's and/or mother's data, excluding early deaths). Genetic risk scores for 19 phenotypes (n=777 proven variants) were also tested. In GWAS, a nicotine receptor locus(CHRNA3, previously associated with increased smoking and lung cancer) was associated with fathers' survival. Less common variants requiring further confirmation were also identified. Offspring of longer lived parents had more protective alleles for coronary artery disease, systolic blood pressure, body mass index, cholesterol and triglyceride levels, type-1 diabetes, inflammatory bowel disease and Alzheimer's disease. In candidate analyses, variants in the TOMM40/APOE locus were associated with longevity, but FOXO variants were not. Associations between extreme longevity (mother >=98 years, fathers >=95 years, n=1,339) and disease alleles were similar, with an additional association with HDL cholesterol (p=5.7x10-3). These results support a multiple protective factors model influencing lifespan and longevity (top 1% survival) in humans, with prominent roles for cardiovascular-related pathways. Several of these genetically influenced risks, including blood pressure and tobacco exposure, are potentially modifiable.

Genetic evidence that lower circulating FSH levels lengthen menstrual cycle, increase age at menopause and impact female reproductive health

STUDY QUESTION

How does a genetic variant in the FSHB promoter, known to alter FSH levels, impact female reproductive health?

SUMMARY ANSWER

The T allele of the FSHB promoter polymorphism (rs10835638; c.-211G>T) results in longer menstrual cycles and later menopause and, while having detrimental effects on fertility, is protective against endometriosis.

WHAT IS KNOWN ALREADY

The FSHB promoter polymorphism (rs10835638; c.-211G>T) affects levels of FSHB transcription and, as a result, circulating levels of FSH. FSH is required for normal fertility and genetic variants at the FSHB locus are associated with age at menopause and polycystic ovary syndrome (PCOS).

STUDY DESIGN, SIZE, DURATION

We used cross-sectional data from the UK Biobank to look at associations between the FSHB promoter polymorphism and reproductive traits, and performed a genome-wide association study (GWAS) for length of menstrual cycle.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We included white British individuals aged 40-69 years in 2006-2010, in the May 2015 release of genetic data from UK Biobank. We tested the FSH-lowering T allele of the FSHB promoter polymorphism (rs10835638; c.-211G>T) for associations with 29, mainly female, reproductive phenotypes in up to 63 350 women and 56 608 men. We conducted a GWAS in 9534 individuals to identify genetic variants associated with length of menstrual cycle.

MAIN RESULTS AND THE ROLE OF CHANCE

The FSH-lowering T allele of the FSHB promoter polymorphism (rs10835638; MAF 0.16) was associated with longer menstrual cycles [0.16 SD (c. 1 day) per minor allele; 95% confidence interval (CI) 0.12-0.20; P = 6 × 10(-16)], later age at menopause (0.13 years per minor allele; 95% CI 0.04-0.22; P = 5.7 × 10(-3)), greater female nulliparity [odds ratio (OR) = 1.06; 95% CI 1.02-1.11; P = 4.8 × 10(-3)] and lower risk of endometriosis (OR = 0.79; 95% CI 0.69-0.90; P = 4.1 × 10(-4)). The FSH-lowering T allele was not associated with other female reproductive illnesses or conditions in our study and we did not replicate associations with male infertility or PCOS. In the GWAS for menstrual cycle length, only variants near the FSHB gene reached genome-wide significance (P < 5 × 10(-9)).

LIMITATIONS, REASONS FOR CAUTION

The data included might be affected by recall bias. Cycle length was not available for 25% of women still cycling (1% did not answer, 6% did not know and for 18% cycle length was recorded as 'irregular'). Women with a cycle length recorded were aged over 40 and were approaching menopause; however, we did not find evidence that this affected the results. Many of the groups with illnesses had relatively small sample sizes and so the study may have been under-powered to detect an effect.

WIDER IMPLICATIONS OF THE FINDINGS

We found a strong novel association between a genetic variant that lowers FSH levels and longer menstrual cycles, at a locus previously robustly associated with age at menopause. The variant was also associated with nulliparity and endometriosis risk. These findings should now be verified in a second independent group of patients. We conclude that lifetime differences in circulating levels of FSH between individuals can influence menstrual cycle length and a range of reproductive outcomes, including menopause timing, infertility, endometriosis and PCOS.

STUDY FUNDING/COMPETING INTERESTS

None.

TRIAL REGISTRATION NUMBER

Not applicable.

Corrigendum: Rare coding variants and X-linked loci associated with age at menarche

Nature Communications 6, Article number: 7756 (2015); Published 4 August 2015; Updated 17 December 2015 In the Results section and in the legend of Table 1 of this Article, the company deCODE genetics, Inc. is incorrectly referred to as ‘Diabetes Epidemiology: collaborative analysis of Diagnostic criteria in Europe’.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Rare coding variants and X-linked loci associated with age at menarche

More than 100 loci have been identified for age at menarche by genome-wide association studies; however, collectively these explain only ∼3% of the trait variance. Here we test two overlooked sources of variation in 192,974 European ancestry women: low-frequency protein-coding variants and X-chromosome variants. Five missense/nonsense variants (in ALMS1/LAMB2/TNRC6A/TACR3/PRKAG1) are associated with age at menarche (minor allele frequencies 0.08-4.6%; effect sizes 0.08-1.25 years per allele; P<5 × 10(-8)). In addition, we identify common X-chromosome loci at IGSF1 (rs762080, P=9.4 × 10(-13)) and FAAH2 (rs5914101, P=4.9 × 10(-10)). Highlighted genes implicate cellular energy homeostasis, post-transcriptional gene silencing and fatty-acid amide signalling. A frequently reported mutation in TACR3 for idiopathic hypogonatrophic hypogonadism (p.W275X) is associated with 1.25-year-later menarche (P=2.8 × 10(-11)), illustrating the utility of population studies to estimate the penetrance of reportedly pathogenic mutations. Collectively, these novel variants explain ∼0.5% variance, indicating that these overlooked sources of variation do not substantially explain the 'missing heritability' of this complex trait.

Genome-wide association study with 1000 genomes imputation identifies signals for nine sex hormone-related phenotypes

Genetic factors contribute strongly to sex hormone levels, yet knowledge of the regulatory mechanisms remains incomplete. Genome-wide association studies (GWAS) have identified only a small number of loci associated with sex hormone levels, with several reproductive hormones yet to be assessed. The aim of the study was to identify novel genetic variants contributing to the regulation of sex hormones. We performed GWAS using genotypes imputed from the 1000 Genomes reference panel. The study used genotype and phenotype data from a UK twin register. We included 2913 individuals (up to 294 males) from the Twins UK study, excluding individuals receiving hormone treatment. Phenotypes were standardised for age, sex, BMI, stage of menstrual cycle and menopausal status. We tested 7 879 351 autosomal SNPs for association with levels of dehydroepiandrosterone sulphate (DHEAS), oestradiol, free androgen index (FAI), follicle-stimulating hormone (FSH), luteinizing hormone (LH), prolactin, progesterone, sex hormone-binding globulin and testosterone. Eight independent genetic variants reached genome-wide significance (P<5 × 10⁻⁸), with minor allele frequencies of 1.3–23.9%. Novel signals included variants for progesterone (P=7.68 × 10⁻¹²), oestradiol (P=1.63 × 10⁻⁸) and FAI (P=1.50 × 10⁻⁸). A genetic variant near the FSHB gene was identified which influenced both FSH (P=1.74 × 10⁻⁸) and LH (P=3.94 × 10⁻⁹) levels. A separate locus on chromosome 7 was associated with both DHEAS (P=1.82 × 10⁻¹⁴) and progesterone (P=6.09 × 10⁻¹⁴). This study highlights loci that are relevant to reproductive function and suggests overlap in the genetic basis of hormone regulation.