Large-scale plasma proteomics comparisons through genetics and disease associations.

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Actionable Genotypes and Their Association with Life Span in Iceland.

BACKGROUND: In 2021, the American College of Medical Genetics and Genomics (ACMG) recommended reporting actionable genotypes in 73 genes associated with diseases for which preventive or therapeutic measures are available. Evaluations of the association of actionable genotypes in these genes with life span are currently lacking. METHODS: We assessed the prevalence of coding and splice variants in genes on the ACMG Secondary Findings, version 3.0 (ACMG SF v3.0), list in the genomes of 57,933 Icelanders. We assigned pathogenicity to all reviewed variants using reported evidence in the ClinVar database, the frequency of variants, and their associations with disease to create a manually curated set of actionable genotypes (variants). We assessed the relationship between these genotypes and life span and further examined the specific causes of death among carriers. RESULTS: Through manual curation of 4405 sequence variants in the ACMG SF v3.0 genes, we identified 235 actionable genotypes in 53 genes. Of the 57,933 participants, 2306 (4.0%) carried at least one actionable genotype. We found shorter median survival among persons carrying actionable genotypes than among noncarriers. Specifically, we found that carrying an actionable genotype in a cancer gene was associated with survival that was 3 years shorter than that among noncarriers, with causes of death among carriers attributed primarily to cancer-related conditions. Furthermore, we found evidence of association between carrying an actionable genotype in certain genes in the cardiovascular disease group and a reduced life span. CONCLUSIONS: On the basis of the ACMG SF v3.0 guidelines, we found that approximately 1 in 25 Icelanders carried an actionable genotype and that carrying such a genotype was associated with a reduced life span. (Funded by deCODE Genetics-Amgen.).

Parental influence on human germline de novo mutations in 1,548 trios from Iceland.

The characterization of mutational processes that generate sequence diversity in the human genome is of paramount importance both to medical genetics and to evolutionary studies. To understand how the age and sex of transmitting parents affect de novo mutations, here we sequence 1,548 Icelanders, their parents, and, for a subset of 225, at least one child, to 35× genome-wide coverage. We find 108,778 de novo mutations, both single nucleotide polymorphisms and indels, and determine the parent of origin of 42,961. The number of de novo mutations from mothers increases by 0.37 per year of age (95% CI 0.32-0.43), a quarter of the 1.51 per year from fathers (95% CI 1.45-1.57). The number of clustered mutations increases faster with the mother's age than with the father's, and the genomic span of maternal de novo mutation clusters is greater than that of paternal ones. The types of de novo mutation from mothers change substantially with age, with a 0.26% (95% CI 0.19-0.33%) decrease in cytosine-phosphate-guanine to thymine-phosphate-guanine (CpG>TpG) de novo mutations and a 0.33% (95% CI 0.28-0.38%) increase in C>G de novo mutations per year, respectively. Remarkably, these age-related changes are not distributed uniformly across the genome. A striking example is a 20 megabase region on chromosome 8p, with a maternal C>G mutation rate that is up to 50-fold greater than the rest of the genome. The age-related accumulation of maternal non-crossover gene conversions also mostly occurs within these regions. Increased sequence diversity and linkage disequilibrium of C>G variants within regions affected by excess maternal mutations indicate that the underlying mutational process has persisted in humans for thousands of years. Moreover, the regional excess of C>G variation in humans is largely shared by chimpanzees, less by gorillas, and is almost absent from orangutans. This demonstrates that sequence diversity in humans results from evolving interactions between age, sex, mutation type, and genomic location.

Evaluation of Large-Scale Proteomics for Prediction of Cardiovascular Events.

Importance: Whether protein risk scores derived from a single plasma sample could be useful for risk assessment for atherosclerotic cardiovascular disease (ASCVD), in conjunction with clinical risk factors and polygenic risk scores, is uncertain. Objective: To develop protein risk scores for ASCVD risk prediction and compare them to clinical risk factors and polygenic risk scores in primary and secondary event populations. Design, Setting, and Participants: The primary analysis was a retrospective study of primary events among 13 540 individuals in Iceland (aged 40-75 years) with proteomics data and no history of major ASCVD events at recruitment (study duration, August 23, 2000 until October 26, 2006; follow-up through 2018). We also analyzed a secondary event population from a randomized, double-blind lipid-lowering clinical trial (2013-2016), consisting of individuals with stable ASCVD receiving statin therapy and for whom proteomic data were available for 6791 individuals. Exposures: Protein risk scores (based on 4963 plasma protein levels and developed in a training set in the primary event population); polygenic risk scores for coronary artery disease and stroke; and clinical risk factors that included age, sex, statin use, hypertension treatment, type 2 diabetes, body mass index, and smoking status at the time of plasma sampling. Main Outcomes and Measures: Outcomes were composites of myocardial infarction, stroke, and coronary heart disease death or cardiovascular death. Performance was evaluated using Cox survival models and measures of discrimination and reclassification that accounted for the competing risk of non-ASCVD death. Results: In the primary event population test set (4018 individuals [59.0% women]; 465 events; median follow-up, 15.8 years), the protein risk score had a hazard ratio (HR) of 1.93 per SD (95% CI, 1.75 to 2.13). Addition of protein risk score and polygenic risk scores significantly increased the C index when added to a clinical risk factor model (C index change, 0.022 [95% CI, 0.007 to 0.038]). Addition of the protein risk score alone to a clinical risk factor model also led to a significantly increased C index (difference, 0.014 [95% CI, 0.002 to 0.028]). Among White individuals in the secondary event population (6307 participants; 432 events; median follow-up, 2.2 years), the protein risk score had an HR of 1.62 per SD (95% CI, 1.48 to 1.79) and significantly increased C index when added to a clinical risk factor model (C index change, 0.026 [95% CI, 0.011 to 0.042]). The protein risk score was significantly associated with major adverse cardiovascular events among individuals of African and Asian ancestries in the secondary event population. Conclusions and Relevance: A protein risk score was significantly associated with ASCVD events in primary and secondary event populations. When added to clinical risk factors, the protein risk score and polygenic risk score both provided statistically significant but modest improvement in discrimination.

A rare IL33 loss-of-function mutation reduces blood eosinophil counts and protects from asthma.

IL-33 is a tissue-derived cytokine that induces and amplifies eosinophilic inflammation and has emerged as a promising new drug target for asthma and allergic disease. Common variants at IL33 and IL1RL1, encoding the IL-33 receptor ST2, associate with eosinophil counts and asthma. Through whole-genome sequencing and imputation into the Icelandic population, we found a rare variant in IL33 (NM_001199640:exon7:c.487-1G>C (rs146597587-C), allele frequency = 0.65%) that disrupts a canonical splice acceptor site before the last coding exon. It is also found at low frequency in European populations. rs146597587-C associates with lower eosinophil counts (ß = -0.21 SD, P = 2.5×10-16, N = 103,104), and reduced risk of asthma in Europeans (OR = 0.47; 95%CI: 0.32, 0.70, P = 1.8×10-4, N cases = 6,465, N controls = 302,977). Heterozygotes have about 40% lower total IL33 mRNA expression than non-carriers and allele-specific analysis based on RNA sequencing and phased genotypes shows that only 20% of the total expression is from the mutated chromosome. In half of those transcripts the mutation causes retention of the last intron, predicted to result in a premature stop codon that leads to truncation of 66 amino acids. The truncated IL-33 has normal intracellular localization but neither binds IL-33R/ST2 nor activates ST2-expressing cells. Together these data demonstrate that rs146597587-C is a loss of function mutation and support the hypothesis that IL-33 haploinsufficiency protects against asthma.

A rare splice donor mutation in the haptoglobin gene associates with blood lipid levels and coronary artery disease.

Common sequence variants at the haptoglobin gene (HP) have been associated with blood lipid levels. Through whole-genome sequencing of 8,453 Icelanders, we discovered a splice donor founder mutation in HP (NM_001126102.1:c.190 + 1G > C, minor allele frequency = 0.56%). This mutation occurs on the HP1 allele of the common copy number variant in HP and leads to a loss of function of HP1. It associates with lower levels of haptoglobin (P = 2.1 × 10-54), higher levels of non-high density lipoprotein cholesterol (ß = 0.26 mmol/l, P = 2.6 × 10-9) and greater risk of coronary artery disease (odds ratio = 1.30, 95% confidence interval: 1.10-1.54, P = 0.0024). Through haplotype analysis and with RNA sequencing, we provide evidence of a causal relationship between one of the two haptoglobin isoforms, namely Hp1, and lower levels of non-HDL cholesterol. Furthermore, we show that the HP1 allele associates with various other quantitative biological traits.

Adiposity-Dependent Regulatory Effects on Multi-tissue Transcriptomes.

Obesity is a global epidemic that is causally associated with a range of diseases, including type 2 diabetes and cardiovascular disease, at the population-level. However, there is marked heterogeneity in obesity-related outcomes among individuals. This might reflect genotype-dependent responses to adiposity. Given that adiposity, measured by BMI, is associated with widespread changes in gene expression and regulatory variants mediate the majority of known complex trait loci, we sought to identify gene-by-BMI (G × BMI) interactions on the regulation of gene expression in a multi-tissue RNA-sequencing (RNA-seq) dataset from the TwinsUK cohort (n = 856). At a false discovery rate of 5%, we identified 16 cis G × BMI interactions (top cis interaction: CHURC1, rs7143432, p = 2.0 × 10(-12)) and one variant regulating 53 genes in trans (top trans interaction: ZNF423, rs3851570, p = 8.2 × 10(-13)), all in adipose tissue. The interactions were adipose-specific and enriched for variants overlapping adipocyte enhancers, and regulated genes were enriched for metabolic and inflammatory processes. We replicated a subset of the interactions in an independent adipose RNA-seq dataset (deCODE genetics, n = 754). We also confirmed the interactions with an alternate measure of obesity, dual-energy X-ray absorptiometry (DXA)-derived visceral-fat-volume measurements, in a subset of TwinsUK individuals (n = 682). The identified G × BMI regulatory effects demonstrate the dynamic nature of gene regulation and reveal a functional mechanism underlying the heterogeneous response to obesity. Additionally, we have provided a web browser allowing interactive exploration of the dataset, including of association between expression, BMI, and G × BMI regulatory effects in four tissues.

Variant ASGR1 Associated with a Reduced Risk of Coronary Artery Disease.

BACKGROUND: Several sequence variants are known to have effects on serum levels of non-high-density lipoprotein (HDL) cholesterol that alter the risk of coronary artery disease. METHODS: We sequenced the genomes of 2636 Icelanders and found variants that we then imputed into the genomes of approximately 398,000 Icelanders. We tested for association between these imputed variants and non-HDL cholesterol levels in 119,146 samples. We then performed replication testing in two populations of European descent. We assessed the effects of an implicated loss-of-function variant on the risk of coronary artery disease in 42,524 case patients and 249,414 controls from five European ancestry populations. An augmented set of genomes was screened for additional loss-of-function variants in a target gene. We evaluated the effect of an implicated variant on protein stability. RESULTS: We found a rare noncoding 12-base-pair (bp) deletion (del12) in intron 4 of ASGR1, which encodes a subunit of the asialoglycoprotein receptor, a lectin that plays a role in the homeostasis of circulating glycoproteins. The del12 mutation activates a cryptic splice site, leading to a frameshift mutation and a premature stop codon that renders a truncated protein prone to degradation. Heterozygous carriers of the mutation (1 in 120 persons in our study population) had a lower level of non-HDL cholesterol than noncarriers, a difference of 15.3 mg per deciliter (0.40 mmol per liter) (P=1.0×10(-16)), and a lower risk of coronary artery disease (by 34%; 95% confidence interval, 21 to 45; P=4.0×10(-6)). In a larger set of sequenced samples from Icelanders, we found another loss-of-function ASGR1 variant (p.W158X, carried by 1 in 1850 persons) that was also associated with lower levels of non-HDL cholesterol (P=1.8×10(-3)). CONCLUSIONS: ASGR1 haploinsufficiency was associated with reduced levels of non-HDL cholesterol and a reduced risk of coronary artery disease. (Funded by the National Institutes of Health and others.).

Multi-nucleotide de novo Mutations in Humans.

Mutation of the DNA molecule is one of the most fundamental processes in biology. In this study, we use 283 parent-offspring trios to estimate the rate of mutation for both single nucleotide variants (SNVs) and short length variants (indels) in humans and examine the mutation process. We found 17812 SNVs, corresponding to a mutation rate of 1.29 × 10-8 per position per generation (PPPG) and 1282 indels corresponding to a rate of 9.29 × 10-10 PPPG. We estimate that around 3% of human de novo SNVs are part of a multi-nucleotide mutation (MNM), with 558 (3.1%) of mutations positioned less than 20kb from another mutation in the same individual (median distance of 525bp). The rate of de novo mutations is greater in late replicating regions (p = 8.29 × 10-19) and nearer recombination events (p = 0.0038) than elsewhere in the genome.

Insertion of an SVA-E retrotransposon into the CASP8 gene is associated with protection against prostate cancer.

Transcriptional and splicing anomalies have been observed in intron 8 of the CASP8 gene (encoding procaspase-8) in association with cutaneous basal-cell carcinoma (BCC) and linked to a germline SNP rs700635. Here, we show that the rs700635[C] allele, which is associated with increased risk of BCC and breast cancer, is protective against prostate cancer [odds ratio (OR) = 0.91, P = 1.0 × 10(-6)]. rs700635[C] is also associated with failures to correctly splice out CASP8 intron 8 in breast and prostate tumours and in corresponding normal tissues. Investigation of rs700635[C] carriers revealed that they have a human-specific short interspersed element-variable number of tandem repeat-Alu (SINE-VNTR-Alu), subfamily-E retrotransposon (SVA-E) inserted into CASP8 intron 8. The SVA-E shows evidence of prior activity, because it has transduced some CASP8 sequences during subsequent retrotransposition events. Whole-genome sequence (WGS) data were used to tag the SVA-E with a surrogate SNP rs1035142[T] (r(2) = 0.999), which showed associations with both the splicing anomalies (P = 6.5 × 10(-32)) and with protection against prostate cancer (OR = 0.91, P = 3.8 × 10(-7)).

Rate of de novo mutations and the importance of father's age to disease risk.

Mutations generate sequence diversity and provide a substrate for selection. The rate of de novo mutations is therefore of major importance to evolution. Here we conduct a study of genome-wide mutation rates by sequencing the entire genomes of 78 Icelandic parent-offspring trios at high coverage. We show that in our samples, with an average father's age of 29.7, the average de novo mutation rate is 1.20 × 10(-8) per nucleotide per generation. Most notably, the diversity in mutation rate of single nucleotide polymorphisms is dominated by the age of the father at conception of the child. The effect is an increase of about two mutations per year. An exponential model estimates paternal mutations doubling every 16.5 years. After accounting for random Poisson variation, father's age is estimated to explain nearly all of the remaining variation in the de novo mutation counts. These observations shed light on the importance of the father's age on the risk of diseases such as schizophrenia and autism.

A Splice Region Variant in LDLR Lowers Non-high Density Lipoprotein Cholesterol and Protects against Coronary Artery Disease.

Through high coverage whole-genome sequencing and imputation of the identified variants into a large fraction of the Icelandic population, we found four independent signals in the low density lipoprotein receptor gene (LDLR) that associate with levels of non-high density lipoprotein cholesterol (non-HDL-C) and coronary artery disease (CAD). Two signals are novel with respect to association with non-HDL-C and are represented by non-coding low frequency variants (between 2-4% frequency), the splice region variant rs72658867-A in intron 14 and rs17248748-T in intron one. These two novel associations were replicated in three additional populations. Both variants lower non-HDL-C levels (rs72658867-A, non-HDL-C effect = -0.44 mmol/l, Padj = 1.1 × 10⁻8° and rs17248748-T, non-HDL-C effect = -0.13 mmol/l, Padj = 1.3 × 10⁻¹²) and confer protection against CAD (rs72658867-A, OR = 0.76 and Padj = 2.7 × 10⁻8 and rs17248748-T, OR = 0.92 and Padj = 0.022). The LDLR splice region variant, rs72658867-A, located at position +5 in intron 14 (NM_000527:c.2140+5G>A), causes retention of intron 14 during transcription and is expected to produce a truncated LDL receptor lacking domains essential for function of the receptor. About half of the transcripts generated from chromosomes carrying rs72658867-A are characterized by this retention of the intron. The same variant also increases LDLR mRNA expression, however, the wild type transcripts do not exceed levels in non-carriers. This demonstrates that sequence variants that disrupt the LDL receptor can lower non-HDL-C and protect against CAD.

A frameshift deletion in the sarcomere gene MYL4 causes early-onset familial atrial fibrillation.

AIMS: Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia in man, causing substantial morbidity and mortality with a major worldwide public health impact. It is increasingly recognized as a highly heritable condition. This study aimed to determine genetic risk factors for early-onset AF. METHODS AND RESULTS: We sequenced the whole genomes of 8453 Icelanders and imputed genotypes of the 25.5 million sequence variants we discovered into 1799 Icelanders with early-onset AF (diagnosed before 60 years of age) and 337 453 controls. Each sequence variant was tested for association based on multiplicative and recessive inheritance models. We discovered a rare frameshift deletion in the myosin MYL4 gene (c.234delC) that associates with early-onset AF under a recessive mode of inheritance (allelic frequency = 0.58%). We found eight homozygous carriers of the mutation, all of whom had early-onset AF. Six of the homozygotes were diagnosed by the age of 30 and the remaining two in their 50s. Three of the homozygotes had received pacemaker implantations due to sick sinus syndrome, three had suffered an ischemic stroke, and one suffered sudden cardiac death. CONCLUSIONS: Through a population approach we found a loss of function mutation in the myosin gene MYL4 that, in the homozygous state, is completely penetrant for early-onset AF. The finding may provide novel mechanistic insight into the pathophysiology of this complex arrhythmia.

Genome-wide association study yields variants at 20p12.2 that associate with urinary bladder cancer.

Genome-wide association studies (GWAS) of urinary bladder cancer (UBC) have yielded common variants at 12 loci that associate with risk of the disease. We report here the results of a GWAS of UBC including 1670 UBC cases and 90 180 controls, followed by replication analysis in additional 5266 UBC cases and 10 456 controls. We tested a dataset containing 34.2 million variants, generated by imputation based on whole-genome sequencing of 2230 Icelanders. Several correlated variants at 20p12, represented by rs62185668, show genome-wide significant association with UBC after combining discovery and replication results (OR = 1.19, P = 1.5 × 10(-11) for rs62185668-A, minor allele frequency = 23.6%). The variants are located in a non-coding region approximately 300 kb upstream from the JAG1 gene, an important component of the Notch signaling pathways that may be oncogenic or tumor suppressive in several forms of cancer. Our results add to the growing number of UBC risk variants discovered through GWAS.

Variants at the Interleukin 1 Gene Locus and Pericarditis.

Importance: Recurrent pericarditis is a treatment challenge and often a debilitating condition. Drugs inhibiting interleukin 1 cytokines are a promising new treatment option, but their use is based on scarce biological evidence and clinical trials of modest sizes, and the contributions of innate and adaptive immune processes to the pathophysiology are incompletely understood. Objective: To use human genomics, transcriptomics, and proteomics to shed light on the pathogenesis of pericarditis. Design, Setting, and Participants: This was a meta-analysis of genome-wide association studies of pericarditis from 5 countries. Associations were examined between the pericarditis-associated variants and pericarditis subtypes (including recurrent pericarditis) and secondary phenotypes. To explore mechanisms, associations with messenger RNA expression (cis-eQTL), plasma protein levels (pQTL), and CpG methylation of DNA (ASM-QTL) were assessed. Data from Iceland (deCODE genetics, 1983-2020), Denmark (Copenhagen Hospital Biobank/Danish Blood Donor Study, 1977-2022), the UK (UK Biobank, 1953-2021), the US (Intermountain, 1996-2022), and Finland (FinnGen, 1970-2022) were included. Data were analyzed from September 2022 to August 2023. Exposure: Genotype. Main Outcomes and Measures: Pericarditis. Results: In this genome-wide association study of 4894 individuals with pericarditis (mean [SD] age at diagnosis, 51.4 [17.9] years, 2734 [67.6%] male, excluding the FinnGen cohort), associations were identified with 2 independent common intergenic variants at the interleukin 1 locus on chromosome 2q14. The lead variant was rs12992780 (T) (effect allele frequency [EAF], 31%-40%; odds ratio [OR], 0.83; 95% CI, 0.79-0.87; P = 6.67 × 10-16), downstream of IL1B and the secondary variant rs7575402 (A or T) (EAF, 45%-55%; adjusted OR, 0.89; 95% CI, 0.85-0.93; adjusted P = 9.6 × 10-8). The lead variant rs12992780 had a smaller odds ratio for recurrent pericarditis (0.76) than the acute form (0.86) (P for heterogeneity = .03) and rs7575402 was associated with CpG methylation overlapping binding sites of 4 transcription factors known to regulate interleukin 1 production: PU.1 (encoded by SPI1), STAT1, STAT3, and CCAAT/enhancer-binding protein ß (encoded by CEBPB). Conclusions and Relevance: This study found an association between pericarditis and 2 independent sequence variants at the interleukin 1 gene locus. This finding has the potential to contribute to development of more targeted and personalized therapy of pericarditis with interleukin 1-blocking drugs.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Sequence variants affecting voice pitch in humans.

The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.