Population-Based Study of Rare Coding Variants in NR5A1/SF-1.

Background: Steroidogenic Factor 1/Nuclear Receptor Subfamily 5 Group A Member 1 (SF-1/NR5A1) is critical for the development and function of sex organs, influencing steroidogenesis and reproduction. While rare deleterious NR5A1/SF-1 variants have been identified in individuals with various differences of sex development (DSD), primary ovarian insufficiency, and infertility, their impact on the general population remains unclear. Methods: We analyzed health records and exome sequencing data from up to 420 162 individuals (227 858 women) from the UK Biobank study to assess the impact of rare (frequency < 0.1%) predicted deleterious NR5A1/SF-1 variants on age at menopause and 26 other traits. Results: No carriers of rare protein truncating variants in NR5A1/SF-1 were identified. We found that the previously reported association of rare deleterious missense NR5A1/SF-1 variants with earlier age at menopause is driven by variants in the DNA binding domain (DBD) and ligand binding domain (LBD) (combined test: beta = -2.36 years/allele, [95% CI: 3.21, -1.51], N = 107 carriers, P = 4.6 × 10-8). Carriers also had a higher risk of adult obesity (OR = 1.061, [95% CI: 1.003, 1.104], N = 344, P = .015), particularly among women (OR = 1.095 [95% CI: 1.034, 1.163, P = 3.87 × 10-3], N = 176), but not men (OR = 1.019, [95% CI: 0.955, 1.088], P = .57, N = 168). Conclusion: Deleterious missense variants in the DBD and LBD likely disrupt NR5A1/SF-1 function. This study broadens the relevance of deleterious NR5A1/SF-1 variants beyond rare DSDs, suggesting the need for extended phenotyping and monitoring of affected individuals.

Genetic links between ovarian ageing, cancer risk and de novo mutation rates.

Human genetic studies of common variants have provided substantial insight into the biological mechanisms that govern ovarian ageing1. Here we report analyses of rare protein-coding variants in 106,973 women from the UK Biobank study, implicating genes with effects around five times larger than previously found for common variants (ETAA1, ZNF518A, PNPLA8, PALB2 and SAMHD1). The SAMHD1 association reinforces the link between ovarian ageing and cancer susceptibility1, with damaging germline variants being associated with extended reproductive lifespan and increased all-cause cancer risk in both men and women. Protein-truncating variants in ZNF518A are associated with shorter reproductive lifespan-that is, earlier age at menopause (by 5.61 years) and later age at menarche (by 0.56 years). Finally, using 8,089 sequenced trios from the 100,000 Genomes Project (100kGP), we observe that common genetic variants associated with earlier ovarian ageing associate with an increased rate of maternally derived de novo mutations. Although we were unable to replicate the finding in independent samples from the deCODE study, it is consistent with the expected role of DNA damage response genes in maintaining the genetic integrity of germ cells. This study provides evidence of genetic links between age of menopause and cancer risk.

High-resolution functional mapping of RAD51C by saturation genome editing.

Pathogenic variants in RAD51C confer an elevated risk of breast and ovarian cancer, while individuals homozygous for specific RAD51C alleles may develop Fanconi anemia. Using saturation genome editing (SGE), we functionally assess 9,188 unique variants, including >99.5% of all possible coding sequence single-nucleotide alterations. By computing changes in variant abundance and Gaussian mixture modeling (GMM), we functionally classify 3,094 variants to be disruptive and use clinical truth sets to reveal an accuracy/concordance of variant classification >99.9%. Cell fitness was the primary assay readout allowing us to observe a phenomenon where specific missense variants exhibit distinct depletion kinetics potentially suggesting that they represent hypomorphic alleles. We further explored our exhaustive functional map, revealing critical residues on the RAD51C structure and resolving variants found in cancer-segregating kindred. Furthermore, through interrogation of UK Biobank and a large multi-center ovarian cancer cohort, we find significant associations between SGE-depleted variants and cancer diagnoses.

QSOX2 Deficiency-induced short stature, gastrointestinal dysmotility and immune dysfunction.

Postnatal growth failure is often attributed to dysregulated somatotropin action, however marked genetic and phenotypic heterogeneity exist. We report five patients from three families who present with short stature, immune dysfunction, atopic eczema and gastrointestinal pathology associated with recessive variants in QSOX2. QSOX2 encodes a nuclear membrane protein linked to disulphide isomerase and oxidoreductase activity. Loss of QSOX2 disrupts Growth hormone-mediated STAT5B nuclear translocation despite enhanced Growth hormone-induced STAT5B phosphorylation. Moreover, patient-derived dermal fibroblasts demonstrate Growth hormone-induced mitochondriopathy and reduced mitochondrial membrane potential. Located at the nuclear membrane, QSOX2 acts as a gatekeeper for regulating stabilisation and import of phosphorylated-STAT5B. Altogether, QSOX2 deficiency modulates human growth by impairing Growth hormone-STAT5B downstream activities and mitochondrial dynamics, which contribute to multi-system dysfunction. Furthermore, our work suggests that therapeutic recombinant insulin-like growth factor-1 may circumvent the Growth hormone-STAT5B dysregulation induced by pathological QSOX2 variants and potentially alleviate organ specific disease.

Damaging mutations in liver X receptor-α are hepatotoxic and implicate cholesterol sensing in liver health.

Liver X receptor-α (LXRα) regulates cellular cholesterol abundance and potently activates hepatic lipogenesis. Here we show that at least 1 in 450 people in the UK Biobank carry functionally impaired mutations in LXRα, which is associated with biochemical evidence of hepatic dysfunction. On a western diet, male and female mice homozygous for a dominant negative mutation in LXRα have elevated liver cholesterol, diffuse cholesterol crystal accumulation and develop severe hepatitis and fibrosis, despite reduced liver triglyceride and no steatosis. This phenotype does not occur on low-cholesterol diets and can be prevented by hepatocyte-specific overexpression of LXRα. LXRα knockout mice exhibit a milder phenotype with regional variation in cholesterol crystal deposition and inflammation inversely correlating with steatosis. In summary, LXRα is necessary for the maintenance of hepatocyte health, likely due to regulation of cellular cholesterol content. The inverse association between steatosis and both inflammation and cholesterol crystallization may represent a protective action of hepatic lipogenesis in the context of excess hepatic cholesterol.

Protein-truncating and rare missense variants in ATM and CHEK2 and associations with cancer in UK Biobank whole-exome sequence data.

BACKGROUND: Deleterious germline variants in ATM and CHEK2 have been associated with a moderately increased risk of breast cancer. Risks for other cancers remain unclear. METHODS: Cancer associations for coding variants in ATM and CHEK2 were evaluated using whole-exome sequence data from UK Biobank linked to cancer registration data (348 488 participants), and analysed both as a retrospective case-control and a prospective cohort study. Odds ratios, hazard ratios, and combined relative risks (RRs) were estimated by cancer type and gene. Separate analyses were performed for protein-truncating variants (PTVs) and rare missense variants (rMSVs; allele frequency <0.1%). RESULTS: PTVs in ATM were associated with increased risks of nine cancers at p<0.001 (pancreas, oesophagus, lung, melanoma, breast, ovary, prostate, bladder, lymphoid leukaemia (LL)), and three at p<0.05 (colon, diffuse non-Hodgkin's lymphoma (DNHL), rectosigmoid junction). Carriers of rMSVs had increased risks of four cancers (p<0.05: stomach, pancreas, prostate, Hodgkin's disease (HD)). RRs were highest for breast, prostate, and any cancer where rMSVs lay in the FAT or PIK domains, and had a Combined Annotation Dependent Depletion score in the highest quintile.PTVs in CHEK2 were associated with three cancers at p<0.001 (breast, prostate, HD) and six at p<0.05 (oesophagus, melanoma, ovary, kidney, DNHL, myeloid leukaemia). Carriers of rMSVs had increased risks of five cancers (p<0.001: breast, prostate, LL; p<0.05: melanoma, multiple myeloma). CONCLUSION: PTVs in ATM and CHEK2 are associated with a wide range of cancers, with the highest RR for pancreatic cancer in ATM PTV carriers. These findings can inform genetic counselling of carriers.

Likely causal effects of insulin resistance and IGF-1 bioaction on childhood and adult adiposity: a Mendelian randomization study.

BACKGROUND: Circulating insulin and insulin-like growth factor-1 (IGF-1) concentrations are positively correlated with adiposity. However, the causal effects of insulin and IGF-1 on adiposity are unclear. METHODS: We performed two-sample Mendelian randomization analyses to estimate the likely causal effects of fasting insulin and IGF-1 on relative childhood adiposity and adult body mass index (BMI). To improve accuracy and biological interpretation, we applied Steiger filtering (to avoid reverse causality) and 'biological effect' filtering of fasting insulin and IGF-1 associated variants. RESULTS: Fasting insulin-increasing alleles (35 variants also associated with higher fasting glucose, indicative of insulin resistance) were associated with lower relative childhood adiposity (P = 3.8 × 10-3) and lower adult BMI (P = 1.4 × 10-5). IGF-1-increasing alleles also associated with taller childhood height (351 variants indicative of greater IGF-1 bioaction) showed no association with relative childhood adiposity (P = 0.077) or adult BMI (P = 0.562). Conversely, IGF-1-increasing alleles also associated with shorter childhood height (306 variants indicative of IGF-1 resistance) were associated with lower relative childhood adiposity (P = 6.7 × 10-3), but effects on adult BMI were inconclusive. CONCLUSIONS: Genetic causal modelling indicates negative effects of insulin resistance on childhood and adult adiposity, and negative effects of IGF-1 resistance on childhood adiposity. Our findings demonstrate the need to distinguish between bioaction and resistance when modelling variants associated with biomarker concentrations.

Loss of transient receptor potential channel 5 causes obesity and postpartum depression.

Hypothalamic neural circuits regulate instinctive behaviors such as food seeking, the fight/flight response, socialization, and maternal care. Here, we identified microdeletions on chromosome Xq23 disrupting the brain-expressed transient receptor potential (TRP) channel 5 (TRPC5). This family of channels detects sensory stimuli and converts them into electrical signals interpretable by the brain. Male TRPC5 deletion carriers exhibited food seeking, obesity, anxiety, and autism, which were recapitulated in knockin male mice harboring a human loss-of-function TRPC5 mutation. Women carrying TRPC5 deletions had severe postpartum depression. As mothers, female knockin mice exhibited anhedonia and depression-like behavior with impaired care of offspring. Deletion of Trpc5 from oxytocin neurons in the hypothalamic paraventricular nucleus caused obesity in both sexes and postpartum depressive behavior in females, while Trpc5 overexpression in oxytocin neurons in knock-in mice reversed these phenotypes. We demonstrate that TRPC5 plays a pivotal role in mediating innate human behaviors fundamental to survival, including food seeking and maternal care.

Saturation genome editing of BAP1 functionally classifies somatic and germline variants.

Many variants that we inherit from our parents or acquire de novo or somatically are rare, limiting the precision with which we can associate them with disease. We performed exhaustive saturation genome editing (SGE) of BAP1, the disruption of which is linked to tumorigenesis and altered neurodevelopment. We experimentally characterized 18,108 unique variants, of which 6,196 were found to have abnormal functions, and then used these data to evaluate phenotypic associations in the UK Biobank. We also characterized variants in a large population-ascertained tumor collection, in cancer pedigrees and ClinVar, and explored the behavior of cancer-associated variants compared to that of variants linked to neurodevelopmental phenotypes. Our analyses demonstrated that disruptive germline BAP1 variants were significantly associated with higher circulating levels of the mitogen IGF-1, suggesting a possible pathological mechanism and therapeutic target. Furthermore, we built a variant classifier with >98% sensitivity and specificity and quantify evidence strengths to aid precision variant interpretation.

Maternal obesity increases hypothalamic miR-505-5p expression in mouse offspring leading to altered fatty acid sensing and increased intake of high-fat food.

In utero exposure to maternal obesity programs increased obesity risk. Animal models show that programmed offspring obesity is preceded by hyperphagia, but the mechanisms that mediate these changes are unknown. Using a mouse model of maternal obesity, we observed increased intake of a high-fat diet (HFD) in offspring of obese mothers that precedes the development of obesity. Through small RNA sequencing, we identified programmed overexpression of hypothalamic miR-505-5p that is established in the fetus, lasts to adulthood and is maintained in hypothalamic neural progenitor cells cultured in vitro. Metabolic hormones and long-chain fatty acids associated with obesity increase miR-505-5p expression in hypothalamic neurons in vitro. We demonstrate that targets of miR-505-5p are enriched in fatty acid metabolism pathways and overexpression of miR-505-5p decreased neuronal fatty acid metabolism in vitro. miR-505-5p targets are associated with increased BMI in human genetic studies. Intra-cerebroventricular injection of miR-505-5p in wild-type mice increased HFD intake, mimicking the phenotype observed in offspring exposed to maternal obesity. Conversely, maternal exercise intervention in an obese mouse pregnancy rescued the programmed increase of hypothalamic miR-505-5p in offspring of obese dams and reduced HFD intake to control offspring levels. This study identifies a novel mechanism by which maternal obesity programs obesity in offspring via increased intake of high-fat foods.

Association between mosaic loss of chromosome Y and pulmonary fibrosis susceptibility and severity.

Background Pulmonary fibrosis (PF) is a rare lung disease with diverse pathogenesis and multiple interconnected underlying biological mechanisms. Mosaic loss of chromosome Y (mLOY) is one of the most common forms of acquired chromosome abnormality in men, which has been reported to be associated with increased risk of various chronic progressive diseases including fibrotic diseases. However, the exact role of mLOY in the development of PF remains elusive and to be elucidated. Methods: We adopted three complementary approaches to explore the role of mLOY in the pathogenesis of PF. We used copy number on chromosome Y to estimate mLOY comparing patients in PROFILE and gnomAD cohorts and between cases and control patients from the GE100KGP cohort. Correlation of mLOY with demographic and clinical variables was tested using patients from PROFILE cohort. Lung single-cell transcriptomic data were analysed to assess the cell types implicated in mLOY. We performed Mendelian randomisation to examine the causal relationship between mLOY, IPF, and telomere length. Results: The genetic analysis suggests that mLOY is found in PF from both case cohorts but when compared with an age matched population the effect is minimal (P = 0.0032). mLOY is related to age (P = 0.00021) and shorter telomere length (P = 0.0081) rather than PF severity or progression. Single-cell analysis indicates that mLOY appears to be found primarily in immune cells and appears to be related to presence and severity of fibrosis. Mendelian randomisation demonstrates that mLOY is not on the causal pathway for IPF, but partial evidence supports that telomere shortening is on the causal pathway for mLOY. Conclusion: Our study confirms the existence of mLOY in PF patients and suggests that mLOY is not a major driver of IPF. The combined evidence suggests a triangulation model where telomere shortening leads to both IPF and mLOY.

Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Metabolome-wide Mendelian randomization for age at menarche and age at natural menopause.

BACKGROUND: The role of metabolism in the variation of age at menarche (AAM) and age at natural menopause (ANM) in the female population is not entirely known. We aimed to investigate the causal role of circulating metabolites in AAM and ANM using Mendelian randomization (MR). METHODS: We combined MR with genetic colocalization to investigate potential causal associations between 658 metabolites and AAM and between 684 metabolites and ANM. We extracted genetic instruments for our exposures from four genome-wide association studies (GWAS) on circulating metabolites and queried the effects of these variants on the outcomes in two large GWAS from the ReproGen consortium. Additionally, we assessed the mediating role of the body mass index (BMI) in these associations, identified metabolic pathways implicated in AAM and ANM, and sought validation for selected metabolites in the Avon Longitudinal Study of Parents and Children (ALSPAC). RESULTS: Our analysis identified 10 candidate metabolites for AAM, but none of them colocalized with AAM. For ANM, 76 metabolites were prioritized (FDR-adjusted MR P-value ≤ 0.05), with 17 colocalizing, primarily in the glycerophosphocholines class, including the omega-3 fatty acid and phosphatidylcholine (PC) categories. Pathway analyses and validation in ALSPAC mothers also highlighted the role of omega and polyunsaturated fatty acids levels in delaying age at menopause. CONCLUSIONS: Our study suggests that metabolites from the glycerophosphocholine and fatty acid families play a causal role in the timing of both menarche and menopause. This underscores the significance of specific metabolic pathways in the biology of female reproductive longevity.

Rare variant associations with birth weight identify genes involved in adipose tissue regulation, placental function and insulin-like growth factor signalling.

Investigating the genetic factors influencing human birth weight may lead to biological insights into fetal growth and long-term health. Genome-wide association studies of birth weight have highlighted associated variants in more than 200 regions of the genome, but the causal genes are mostly unknown. Rare genetic variants with robust evidence of association are more likely to point to causal genes, but to date, only a few rare variants are known to influence birth weight. We aimed to identify genes that harbour rare variants that impact birth weight when carried by either the fetus or the mother, by analysing whole exome sequence data in UK Biobank participants. We annotated rare (minor allele frequency <0.1%) protein-truncating or high impact missense variants on whole exome sequence data in up to 234,675 participants with data on their own birth weight (fetal variants), and up to 181,883 mothers who reported the birth weight of their first child (maternal variants). Variants within each gene were collapsed to perform gene burden tests and for each associated gene, we compared the observed fetal and maternal effects. We identified 8 genes with evidence of rare fetal variant effects on birth weight, of which 2 also showed maternal effects. One additional gene showed evidence of maternal effects only. We observed 10/11 directionally concordant associations in an independent sample of up to 45,622 individuals (sign test P=0.01). Of the genes identified, IGF1R and PAPPA2 (fetal and maternal-acting) have known roles in insulin-like growth factor bioavailability and signalling. PPARG, INHBE and ACVR1C (all fetal-acting) have known roles in adipose tissue regulation and rare variants in the latter two also showed associations with favourable adiposity patterns in adults. We highlight the dual role of PPARG in both adipocyte differentiation and placental angiogenesis. NOS3, NRK, and ADAMTS8 (fetal and maternal-acting) have been implicated in both placental function and hypertension. Analysis of rare coding variants has identified regulators of fetal adipose tissue and fetoplacental angiogenesis as determinants of birth weight, as well as further evidence for the role of insulin-like growth factors.

Protein-truncating variants in BSN are associated with severe adult-onset obesity, type 2 diabetes and fatty liver disease.

Obesity is a major risk factor for many common diseases and has a substantial heritable component. To identify new genetic determinants, we performed exome-sequence analyses for adult body mass index (BMI) in up to 587,027 individuals. We identified rare loss-of-function variants in two genes (BSN and APBA1) with effects substantially larger than those of well-established obesity genes such as MC4R. In contrast to most other obesity-related genes, rare variants in BSN and APBA1 were not associated with normal variation in childhood adiposity. Furthermore, BSN protein-truncating variants (PTVs) magnified the influence of common genetic variants associated with BMI, with a common variant polygenic score exhibiting an effect twice as large in BSN PTV carriers than in noncarriers. Finally, we explored the plasma proteomic signatures of BSN PTV carriers as well as the functional consequences of BSN deletion in human induced pluripotent stem cell-derived hypothalamic neurons. Collectively, our findings implicate degenerative processes in synaptic function in the etiology of adult-onset obesity.

Loss of GPR75 protects against non-alcoholic fatty liver disease and body fat accumulation.

Approximately 1 in 4 people worldwide have non-alcoholic fatty liver disease (NAFLD); however, there are currently no medications to treat this condition. This study investigated the role of adiposity-associated orphan G protein-coupled receptor 75 (GPR75) in liver lipid accumulation. We profiled Gpr75 expression and report that it is most abundant in the brain. Next, we generated the first single-cell-level analysis of Gpr75 and identified a subpopulation co-expressed with key appetite-regulating hypothalamic neurons. CRISPR-Cas9-deleted Gpr75 mice fed a palatable western diet high in fat adjusted caloric intake to remain in energy balance, thereby preventing NAFLD. Consistent with mouse results, analysis of whole-exome sequencing data from 428,719 individuals (UK Biobank) revealed that variants in GPR75 are associated with a reduced likelihood of hepatic steatosis. Here, we provide a significant advance in understanding of the expression and function of GPR75, demonstrating that it is a promising pharmaceutical target for NAFLD treatment.

Mendelian randomization identifies circulating proteins as biomarkers for age at menarche and age at natural menopause.

Age at menarche (AAM) and age at natural menopause (ANM) are highly heritable traits and have been linked to various health outcomes. We aimed to identify circulating proteins associated with altered ANM and AAM using an unbiased two-sample Mendelian randomization (MR) and colocalization approach. By testing causal effects of 1,271 proteins on AAM, we identified 22 proteins causally associated with AAM in MR, among which 13 proteins (GCKR, FOXO3, SEMA3G, PATE4, AZGP1, NEGR1, LHB, DLK1, ANXA2, YWHAB, DNAJB12, RMDN1 and HPGDS) colocalized. Among 1,349 proteins tested for causal association with ANM using MR, we identified 19 causal proteins among which 7 proteins (CPNE1, TYMP, DNER, ADAMTS13, LCT, ARL and PLXNA1) colocalized. Follow-up pathway and gene enrichment analyses demonstrated links between AAM-related proteins and obesity and diabetes, and between AAM and ANM-related proteins and various types of cancer. In conclusion, we identified proteomic signatures of reproductive ageing in women, highlighting biological processes at both ends of the reproductive lifespan.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Metabolic control of puberty: 60 years in the footsteps of Kennedy and Mitra's seminal work.

An individual's nutritional status has a powerful effect on sexual maturation. Puberty onset is delayed in response to chronic energy insufficiency and is advanced under energy abundance. The consequences of altered pubertal timing for human health are profound. Late puberty increases the chances of cardiometabolic, musculoskeletal and neurocognitive disorders, whereas early puberty is associated with increased risks of adult obesity, type 2 diabetes mellitus, cardiovascular diseases and various cancers, such as breast, endometrial and prostate cancer. Kennedy and Mitra's trailblazing studies, published in 1963 and using experimental models, were the first to demonstrate that nutrition is a key factor in puberty onset. Building on this work, the field has advanced substantially in the past decade, which is largely due to the impressive development of molecular tools for experimentation and population genetics. In this Review, we discuss the latest advances in basic and translational sciences underlying the nutritional and metabolic control of pubertal development, with a focus on perspectives and future directions.

Saturation genome editing of DDX3X clarifies pathogenicity of germline and somatic variation.

Loss-of-function of DDX3X is a leading cause of neurodevelopmental disorders (NDD) in females. DDX3X is also a somatically mutated cancer driver gene proposed to have tumour promoting and suppressing effects. We perform saturation genome editing of DDX3X, testing in vitro the functional impact of 12,776 nucleotide variants. We identify 3432 functionally abnormal variants, in three distinct classes. We train a machine learning classifier to identify functionally abnormal variants of NDD-relevance. This classifier has at least 97% sensitivity and 99% specificity to detect variants pathogenic for NDD, substantially out-performing in silico predictors, and resolving up to 93% of variants of uncertain significance. Moreover, functionally-abnormal variants can account for almost all of the excess nonsynonymous DDX3X somatic mutations seen in DDX3X-driven cancers. Systematic maps of variant effects generated in experimentally tractable cell types have the potential to transform clinical interpretation of both germline and somatic disease-associated variation.