A unified framework for estimating country-specific cumulative incidence for 18 diseases stratified by polygenic risk.

Polygenic scores (PGSs) offer the ability to predict genetic risk for complex diseases across the life course; a key benefit over short-term prediction models. To produce risk estimates relevant to clinical and public health decision-making, it is important to account for varying effects due to age and sex. Here, we develop a novel framework to estimate country-, age-, and sex-specific estimates of cumulative incidence stratified by PGS for 18 high-burden diseases. We integrate PGS associations from seven studies in four countries (N = 1,197,129) with disease incidences from the Global Burden of Disease. PGS has a significant sex-specific effect for asthma, hip osteoarthritis, gout, coronary heart disease and type 2 diabetes (T2D), with all but T2D exhibiting a larger effect in men. PGS has a larger effect in younger individuals for 13 diseases, with effects decreasing linearly with age. We show for breast cancer that, relative to individuals in the bottom 20% of polygenic risk, the top 5% attain an absolute risk for screening eligibility 16.3 years earlier. Our framework increases the generalizability of results from biobank studies and the accuracy of absolute risk estimates by appropriately accounting for age- and sex-specific PGS effects. Our results highlight the potential of PGS as a screening tool which may assist in the early prevention of common diseases.

Generation Scotland: an update on Scotland's longitudinal family health study.

PURPOSE: Generation Scotland (GS) is a large family-based cohort study established as a longitudinal resource for research into the genetic, lifestyle and environmental determinants of physical and mental health. It comprises extensive genetic, sociodemographic and clinical data from volunteers in Scotland. PARTICIPANTS: A total of 24 084 adult participants, including 5501 families, were recruited between 2006 and 2011. Within the cohort, 59% (approximately 14 209) are women, with an average age at recruitment of 49 years. Participants completed a health questionnaire and attended an in-person clinic visit, where detailed baseline data were collected on lifestyle information, cognitive function, personality traits and mental and physical health. Genotype array data are available for 20 026 (83%) participants, and blood-based DNA methylation (DNAm) data for 18 869 (78%) participants. Linkage to routine National Health Service datasets has been possible for 93% (n=22 402) of the cohort, creating a longitudinal resource that includes primary care, hospital attendance, prescription and mortality records. Multimodal brain imaging is available in 1069 individuals. FINDINGS TO DATE: GS has been widely used by researchers across the world to study the genetic and environmental basis of common complex diseases. Over 350 peer-reviewed papers have been published using GS data, contributing to research areas such as ageing, cancer, cardiovascular disease and mental health. Recontact studies have built on the GS cohort to collect additional prospective data to study chronic pain, major depressive disorder and COVID-19. FUTURE PLANS: To create a larger, richer, longitudinal resource, 'Next Generation Scotland' launched in May 2022 to expand the existing cohort by a target of 20 000 additional volunteers, now including anyone aged 12+ years. New participants complete online consent and questionnaires and provide postal saliva samples, from which genotype and salivary DNAm array data will be generated. The latest cohort information and how to access data can be found on the GS website (www.generationscotland.org).

The Effect of Mediterranean Diet on Thyroid Gland Activity.

The main goal of this research was to determine whether there is a correlation between adherence to the Mediterranean diet (assessed by the Mediterranean Diet Serving Score (MDSS)) and parameters indicating thyroid gland activity, such as concentration of thyroid-stimulating hormone (TSH), thyroid hormones (free triiodothyronine (fT3), free thyroxine (fT4)), thyroglobulin (Tg), antibodies to thyroid proteins (thyroglobulin antibodies (TgAb) and thyroid peroxidase antibodies (TPOAb)), and calcitonin (CT) in plasma and serum samples. An additional objective was to investigate whether there are differences in the values of the MDSS among clinical groups (euthyroid individuals, euthyroid individuals with positive TgAb and/or TPOAb, and hypothyroid and hyperthyroid participants). This cross-sectional study included 4620 participants over 18 years of age from the islands of Korcula and Vis, and the mainland city of Split. The MDSS was assessed from a food frequency questionnaire (FFQ). MDSS values were significantly higher in females compared to males and showed a positive association with the age of the participants. There was no significant difference in the MDSS values among the examined clinical groups. In the group of subjects with euthyroidism, a significant positive association was found between fT3 and the MDSS, while in the group of subjects with subclinical hypothyroidism, a significant positive association was observed between the MDSS and both fT3 and fT4. CT levels were also positively associated with the MDSS. Considering the significant positive association of the MDSS and both fT3 and fT4 levels in patients with subclinical hypothyroidism, the results of this study could be used to create guidelines for selecting an appropriate, potentially protective diet for these patients.

Copy-number variants and polygenic risk for intelligence confer risk for autism spectrum disorder irrespective of their effects on cognitive ability.

Introduction: Rare copy number variants (CNVs) and polygenic risk for intelligence (PRS-IQ) both confer susceptibility for autism spectrum disorder (ASD) but have opposing effects on cognitive ability. The field has struggled to disentangle the effects of these two classes of genomic variants on cognitive ability from their effects on ASD susceptibility, in part because previous studies did not include controls with cognitive measures. We aim to investigate the impact of these genomic variants on ASD risk while adjusting for their known effects on cognitive ability. Methods: In a cohort of 8,426 subjects with ASD and 169,804 controls with cognitive assessments, we found that rare coding CNVs and PRS-IQ increased ASD risk, even after adjusting for their effects on cognitive ability. Results: Bottom decile PRS-IQ and CNVs both decreased cognitive ability but had opposing effects on ASD risk. Models combining both classes of variants showed that the effects of rare CNVs and PRS-IQ on ASD risk and cognitive ability were largely additive, further suggesting that susceptibility for ASD is conferred independently from its effects on cognitive ability. Despite imparting mostly additive effects on ASD risk, rare CNVs and PRS-IQ showed opposing effects on core and associated features and developmental history among subjects with ASD. Discussion: Our findings suggest that cognitive ability itself may not be the factor driving the underlying liability for ASD conferred by these two classes of genomic variants. In other words, ASD risk and cognitive ability may be two distinct manifestations of CNVs and PRS-IQ. This study also highlights the challenge of understanding how genetic risk for ASD maps onto its dimensional traits.

Blood-based epigenome-wide analyses of chronic low-grade inflammation across diverse population cohorts.

Chronic inflammation is a hallmark of age-related disease states. The effectiveness of inflammatory proteins including C-reactive protein (CRP) in assessing long-term inflammation is hindered by their phasic nature. DNA methylation (DNAm) signatures of CRP may act as more reliable markers of chronic inflammation. We show that inter-individual differences in DNAm capture 50% of the variance in circulating CRP (N = 17,936, Generation Scotland). We develop a series of DNAm predictors of CRP using state-of-the-art algorithms. An elastic-net-regression-based predictor outperformed competing methods and explained 18% of phenotypic variance in the Lothian Birth Cohort of 1936 (LBC1936) cohort, doubling that of existing DNAm predictors. DNAm predictors performed comparably in four additional test cohorts (Avon Longitudinal Study of Parents and Children, Health for Life in Singapore, Southall and Brent Revisited, and LBC1921), including for individuals of diverse genetic ancestry and different age groups. The best-performing predictor surpassed assay-measured CRP and a genetic score in its associations with 26 health outcomes. Our findings forge new avenues for assessing chronic low-grade inflammation in diverse populations.

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

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

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Multiple Cardiac Biomarkers to Improve Prediction of Cardiovascular Events: Findings from the Generation Scotland Scottish Family Health Study.

BACKGROUND: Many studies have investigated whether single cardiac biomarkers improve cardiovascular risk prediction for primary prevention but whether a combined approach could further improve risk prediction is unclear. We aimed to test a sex-specific, combined cardiac biomarker approach for cardiovascular risk prediction. METHODS: In the Generation Scotland Scottish Family Health Study, N-terminal pro-B-type natriuretic peptide (NT-proBNP), growth differentiation factor-15 (GDF-15), cardiac troponin I (cTnI), cardiac troponin T (cTnT), and C-reactive protein (CRP) were measured in stored serum using automated immunoassays. Sex-specific Cox models that included SCORE2 risk factors evaluated addition of single and combined biomarkers for prediction of major adverse cardiovascular events (MACE). Combined biomarker models were compared to a baseline model that included SCORE2 risk factors. RESULTS: The study population comprised 18 383 individuals (58.9% women, median age of 48 years [25th-75th percentile, 35-58 years]). During the median follow up of 11.6 (25th-75th percentile, 10.8-13.0) years, MACE occurred in 942 (5.1%) individuals. The greatest increase in discrimination with addition of individual biomarkers to the base model was for women GDF-15 and for men NT-proBNP (change in c-index: + 0.010 for women and +0.005 for men). For women, combined biomarker models that included GDF-15 and NT-proBNP (+0.012) or GDF-15 and cTnI (+0.013), but not CRP or cTnT, further improved discrimination. For men, combined biomarker models that included NT-proBNP and GDF-15 (+0.007), NT-proBNP and cTnI (+0.006), or NT-proBNP and CRP (+0.008), but not cTnT, further improved discrimination. CONCLUSIONS: A combined biomarker approach, particularly the use of GDF-15, NT-proBNP and cTnI, further refined cardiovascular risk estimates.

Copy-number variants and polygenic risk for intelligence confer risk for autism spectrum disorder irrespective of their effects on cognitive ability.

Rare copy number variants (CNVs) and polygenic risk for intelligence (PRS-IQ) both confer risk for autism spectrum disorder (ASD) but have opposing effects on cognitive ability. The field has struggled to disentangle the effects of these two classes of genomic variants on cognitive ability from their effects on ASD risk, in part because previous studies did not include controls with cognitive measures. We aim to investigate the impact of these genomic variants on ASD risk while adjusting for their known effects on cognitive ability. In a cohort of 8,426 subjects with ASD and 169,804 controls with cognitive assessments, we found that rare coding CNVs and PRS-IQ increased ASD risk, even after adjusting for their effects on cognitive ability. Bottom decile PRS-IQ and CNVs both decreased cognitive ability but had opposing effects on ASD risk. Models combining both classes of variants showed that the effects of rare CNVs and PRS-IQ on ASD risk and cognitive ability were largely additive, further suggesting that risk for ASD is conferred independently from its effects on cognitive ability. Despite imparting mostly additive effects on ASD risk, rare CNVs and PRS-IQ showed opposing effects on core and associated features and developmental history among subjects with ASD. Our findings suggest that cognitive ability itself may not be the factor driving the underlying risk for ASD conferred by these two classes of genomic variants. In other words, ASD risk and cognitive ability may be two distinct manifestations of CNVs and PRS-IQ. This study also highlights the challenge of understanding how genetic risk for ASD maps onto its dimensional traits.

Mapping of the gene network that regulates glycan clock of ageing.

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

Blood-based epigenome-wide analyses of 19 common disease states: A longitudinal, population-based linked cohort study of 18,413 Scottish individuals.

BACKGROUND: DNA methylation is a dynamic epigenetic mechanism that occurs at cytosine-phosphate-guanine dinucleotide (CpG) sites. Epigenome-wide association studies (EWAS) investigate the strength of association between methylation at individual CpG sites and health outcomes. Although blood methylation may act as a peripheral marker of common disease states, previous EWAS have typically focused only on individual conditions and have had limited power to discover disease-associated loci. This study examined the association of blood DNA methylation with the prevalence of 14 disease states and the incidence of 19 disease states in a single population of over 18,000 Scottish individuals. METHODS AND FINDINGS: DNA methylation was assayed at 752,722 CpG sites in whole-blood samples from 18,413 volunteers in the family-structured, population-based cohort study Generation Scotland (age range 18 to 99 years). EWAS tested for cross-sectional associations between baseline CpG methylation and 14 prevalent disease states, and for longitudinal associations between baseline CpG methylation and 19 incident disease states. Prevalent cases were self-reported on health questionnaires at the baseline. Incident cases were identified using linkage to Scottish primary (Read 2) and secondary (ICD-10) care records, and the censoring date was set to October 2020. The mean time-to-diagnosis ranged from 5.0 years (for chronic pain) to 11.7 years (for Coronavirus Disease 2019 (COVID-19) hospitalisation). The 19 disease states considered in this study were selected if they were present on the World Health Organisation's 10 leading causes of death and disease burden or included in baseline self-report questionnaires. EWAS models were adjusted for age at methylation typing, sex, estimated white blood cell composition, population structure, and 5 common lifestyle risk factors. A structured literature review was also conducted to identify existing EWAS for all 19 disease states tested. The MEDLINE, Embase, Web of Science, and preprint servers were searched to retrieve relevant articles indexed as of March 27, 2023. Fifty-four of approximately 2,000 indexed articles met our inclusion criteria: assayed blood-based DNA methylation, had >20 individuals in each comparison group, and examined one of the 19 conditions considered. First, we assessed whether the associations identified in our study were reported in previous studies. We identified 69 associations between CpGs and the prevalence of 4 conditions, of which 58 were newly described. The conditions were breast cancer, chronic kidney disease, ischemic heart disease, and type 2 diabetes mellitus. We also uncovered 64 CpGs that associated with the incidence of 2 disease states (COPD and type 2 diabetes), of which 56 were not reported in the surveyed literature. Second, we assessed replication across existing studies, which was defined as the reporting of at least 1 common site in >2 studies that examined the same condition. Only 6/19 disease states had evidence of such replication. The limitations of this study include the nonconsideration of medication data and a potential lack of generalizability to individuals that are not of Scottish and European ancestry. CONCLUSIONS: We discovered over 100 associations between blood methylation sites and common disease states, independently of major confounding risk factors, and a need for greater standardisation among EWAS on human disease.

Educational attainment, health outcomes and mortality: a within-sibship Mendelian randomization study.

BACKGROUND: Previous Mendelian randomization (MR) studies using population samples (population MR) have provided evidence for beneficial effects of educational attainment on health outcomes in adulthood. However, estimates from these studies may have been susceptible to bias from population stratification, assortative mating and indirect genetic effects due to unadjusted parental genotypes. MR using genetic association estimates derived from within-sibship models (within-sibship MR) can avoid these potential biases because genetic differences between siblings are due to random segregation at meiosis. METHODS: Applying both population and within-sibship MR, we estimated the effects of genetic liability to educational attainment on body mass index (BMI), cigarette smoking, systolic blood pressure (SBP) and all-cause mortality. MR analyses used individual-level data on 72 932 siblings from UK Biobank and the Norwegian HUNT study, and summary-level data from a within-sibship Genome-wide Association Study including >140 000 individuals. RESULTS: Both population and within-sibship MR estimates provided evidence that educational attainment decreased BMI, cigarette smoking and SBP. Genetic variant-outcome associations attenuated in the within-sibship model, but genetic variant-educational attainment associations also attenuated to a similar extent. Thus, within-sibship and population MR estimates were largely consistent. The within-sibship MR estimate of education on mortality was imprecise but consistent with a putative effect. CONCLUSIONS: These results provide evidence of beneficial individual-level effects of education (or liability to education) on adulthood health, independently of potential demographic and family-level confounders.

Genome-wide association study of chronic sputum production implicates loci involved in mucus production and infection.

BACKGROUND: Chronic sputum production impacts on quality of life and is a feature of many respiratory diseases. Identification of the genetic variants associated with chronic sputum production in a disease agnostic sample could improve understanding of its causes and identify new molecular targets for treatment. METHODS: We conducted a genome-wide association study (GWAS) of chronic sputum production in UK Biobank. Signals meeting genome-wide significance (p<5×10-8) were investigated in additional independent studies, were fine-mapped and putative causal genes identified by gene expression analysis. GWASs of respiratory traits were interrogated to identify whether the signals were driven by existing respiratory disease among the cases and variants were further investigated for wider pleiotropic effects using phenome-wide association studies (PheWASs). RESULTS: From a GWAS of 9714 cases and 48 471 controls, we identified six novel genome-wide significant signals for chronic sputum production including signals in the human leukocyte antigen (HLA) locus, chromosome 11 mucin locus (containing MUC2, MUC5AC and MUC5B) and FUT2 locus. The four common variant associations were supported by independent studies with a combined sample size of up to 2203 cases and 17 627 controls. The mucin locus signal had previously been reported for association with moderate-to-severe asthma. The HLA signal was fine-mapped to an amino acid change of threonine to arginine (frequency 36.8%) in HLA-DRB1 (HLA-DRB1*03:147). The signal near FUT2 was associated with expression of several genes including FUT2, for which the direction of effect was tissue dependent. Our PheWAS identified a wide range of associations including blood cell traits, liver biomarkers, infections, gastrointestinal and thyroid-associated diseases, and respiratory disease. CONCLUSIONS: Novel signals at the FUT2 and mucin loci suggest that mucin fucosylation may be a driver of chronic sputum production even in the absence of diagnosed respiratory disease and provide genetic support for this pathway as a target for therapeutic intervention.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans.

BACKGROUND: The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging. RESULTS: We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks. CONCLUSIONS: These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Genome-wide aggregated trans-effects on risk of type 1 diabetes: A test of the "omnigenic" sparse effector hypothesis of complex trait genetics.

The "omnigenic" hypothesis postulates that the polygenic effects of common SNPs on a typical complex trait are mediated through trans-effects on expression of a relatively sparse set of effector ("core") genes. We tested this hypothesis in a study of 4,964 cases of type 1 diabetes (T1D) and 7,497 controls by using summary statistics to calculate aggregated (excluding the HLA region) trans-scores for gene expression in blood. From associations of T1D with aggregated trans-scores, nine putative core genes were identified, of which three-STAT1, CTLA4 and FOXP3-are genes in which variants cause monogenic forms of autoimmune diabetes. Seven of these genes affect the activity of regulatory T cells, and two are involved in immune responses to microbial lipids. Four T1D-associated genomic regions could be identified as master regulators via trans-effects on gene expression. These results support the sparse effector hypothesis and reshape our understanding of the genetic architecture of T1D.

Genome-Wide Meta-Analysis Identifies Multiple Novel Rare Variants to Predict Common Human Infectious Diseases Risk.

Infectious diseases still threaten global human health, and host genetic factors have been indicated as determining risk factors for observed variations in disease susceptibility, severity, and outcome. We performed a genome-wide meta-analysis on 4624 subjects from the 10,001 Dalmatians cohort, with 14 infection-related traits. Despite a rather small number of cases in some instances, we detected 29 infection-related genetic associations, mostly belonging to rare variants. Notably, the list included the genes CD28, INPP5D, ITPKB, MACROD2, and RSF1, all of which have known roles in the immune response. Expanding our knowledge on rare variants could contribute to the development of genetic panels that could assist in predicting an individual's life-long susceptibility to major infectious diseases. In addition, longitudinal biobanks are an interesting source of information for identifying the host genetic variants involved in infectious disease susceptibility and severity. Since infectious diseases continue to act as a selective pressure on our genomes, there is a constant need for a large consortium of biobanks with access to genetic and environmental data to further elucidate the complex mechanisms behind host-pathogen interactions and infectious disease susceptibility.

Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus.

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.