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.

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.

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.

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.

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.

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.

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.

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.

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.

Large-scale exome sequence analysis identifies sex- and age-specific determinants of obesity.

Obesity contributes substantially to the global burden of disease and has a significant heritable component. Recent large-scale exome sequencing studies identified several genes in which rare, protein-coding variants have large effects on adult body mass index (BMI). Here we extended such work by performing sex-stratified associations in the UK Biobank study (N∼420,000). We identified genes in which rare heterozygous loss-of-function increases adult BMI in women (DIDO1, PTPRG, and SLC12A5) and in men (SLTM), with effect sizes up to ∼8 kg/m2. This is complemented by analyses implicating rare variants in OBSCN and MADD for recalled childhood adiposity. The known functions of these genes, as well as findings of common variant genome-wide pathway enrichment analyses, suggest a role for neuron death, apoptosis, and DNA damage response mechanisms in the susceptibility to obesity across the life-course. These findings highlight the importance of considering sex-specific and life-course effects in the genetic regulation of obesity.

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.

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency.

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

Population analyses of mosaic X chromosome loss identify genetic drivers and widespread signatures of cellular selection.

Mosaic loss of the X chromosome (mLOX) is the most commonly occurring clonal somatic alteration detected in the leukocytes of women, yet little is known about its genetic determinants or phenotypic consequences. To address this, we estimated mLOX in >900,000 women across eight biobanks, identifying 10% of women with detectable X loss in approximately 2% of their leukocytes. Out of 1,253 diseases examined, women with mLOX had an elevated risk of myeloid and lymphoid leukemias and pneumonia. Genetic analyses identified 49 common variants influencing mLOX, implicating genes with established roles in chromosomal missegregation, cancer predisposition, and autoimmune diseases. Complementary exome-sequence analyses identified rare missense variants in FBXO10 which confer a two-fold increased risk of mLOX. A small fraction of these associations were shared with mosaic Y chromosome loss in men, suggesting different biological processes drive the formation and clonal expansion of sex chromosome missegregation events. Allelic shift analyses identified alleles on the X chromosome which are preferentially retained, demonstrating that variation at many loci across the X chromosome is under cellular selection. A novel polygenic score including 44 independent X chromosome allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Collectively our results support a model where germline variants predispose women to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of subsequent clonal expansion.

Large scale phenotype imputation and in vivo functional validation implicate ADAMTS14 as an adiposity gene.

Obesity remains an unmet global health burden. Detrimental anatomical distribution of body fat is a major driver of obesity-mediated mortality risk and is demonstrably heritable. However, our understanding of the full genetic contribution to human adiposity is incomplete, as few studies measure adiposity directly. To address this, we impute whole-body imaging adiposity phenotypes in UK Biobank from the 4,366 directly measured participants onto the rest of the cohort, greatly increasing our discovery power. Using these imputed phenotypes in 392,535 participants yielded hundreds of genome-wide significant associations, six of which replicate in independent cohorts. The leading causal gene candidate, ADAMTS14, is further investigated in a mouse knockout model. Concordant with the human association data, the Adamts14-/- mice exhibit reduced adiposity and weight-gain under obesogenic conditions, alongside an improved metabolic rate and health. Thus, we show that phenotypic imputation at scale offers deeper biological insights into the genetics of human adiposity that could lead to therapeutic targets.

Damaging missense variants in IGF1R implicate a role for IGF-1 resistance in the etiology of type 2 diabetes.

Type 2 diabetes (T2D) is a heritable metabolic disorder. While population studies have identified hundreds of common genetic variants associated with T2D, the role of rare (frequency < 0.1%) protein-coding variation is less clear. We performed exome sequence analysis in 418,436 (n = 32,374 T2D cases) individuals in the UK Biobank. We identified previously reported genes (GCK, GIGYF1, HNF1A) in addition to missense variants in ZEB2 (n = 31 carriers; odds ratio [OR] = 5.5 [95% confidence interval = 2.5-12.0]; p = 6.4 × 10-7), MLXIPL (n = 245; OR = 2.3 [1.6-3.2]; p = 3.2 × 10-7), and IGF1R (n = 394; OR = 2.4 [1.8-3.2]; p = 1.3 × 10-10). Carriers of damaging missense variants within IGF1R were also shorter (-2.2 cm [-1.8 to -2.7]; p = 1.2 × 10-19) and had higher circulating insulin-like growth factor-1 (IGF-1) protein levels (2.3 nmol/L [1.7-2.9]; p = 2.8 × 10-14), indicating relative IGF-1 resistance. A likely causal role of IGF-1 resistance was supported by Mendelian randomization analyses using common variants. These results increase understanding of the genetic architecture of T2D and highlight the growth hormone/IGF-1 axis as a potential therapeutic target.

HOTAIR interacts with PRC2 complex regulating the regional preadipocyte transcriptome and human fat distribution.

Mechanisms governing regional human adipose tissue (AT) development remain undefined. Here, we show that the long non-coding RNA HOTAIR (HOX transcript antisense RNA) is exclusively expressed in gluteofemoral AT, where it is essential for adipocyte development. We find that HOTAIR interacts with polycomb repressive complex 2 (PRC2) and we identify core HOTAIR-PRC2 target genes involved in adipocyte lineage determination. Repression of target genes coincides with PRC2 promoter occupancy and H3K27 trimethylation. HOTAIR is also involved in modifying the gluteal adipocyte transcriptome through alternative splicing. Gluteal-specific expression of HOTAIR is maintained by defined regions of open chromatin across the HOTAIR promoter. HOTAIR expression levels can be modified by hormonal (estrogen, glucocorticoids) and genetic variation (rs1443512 is a HOTAIR eQTL associated with reduced gynoid fat mass). These data identify HOTAIR as a dynamic regulator of the gluteal adipocyte transcriptome and epigenome with functional importance for human regional AT development.