ABSTRACT
Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor two conditionally distinct, non-coding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR=0.180; MAFEUR=0.008). Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation.
ABSTRACT
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.
Subject(s)
Gene Frequency , Menarche , Puberty , Humans , Female , Menarche/genetics , Puberty/genetics , Animals , Multifactorial Inheritance/genetics , Mice , Genome-Wide Association Study , Adolescent , Puberty, Precocious/genetics , Polymorphism, Single Nucleotide , Receptors, G-Protein-Coupled/genetics , Puberty, Delayed/genetics , ChildABSTRACT
Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to 23 genes. Investigating these genes' functional implications shed light on neurological, thyroidal, bone metabolism, and hematopoietic pathways that necessitate future investigation for blood pressure management that caters to sleep health lifestyle. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausible nature of distinct influences of both sleep duration extremes in cardiovascular health. Several of our loci are specific towards a particular population background or sex, emphasizing the importance of addressing heterogeneity entangled in gene-environment interactions, when considering precision medicine design approaches for blood pressure management.
ABSTRACT
BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.
Subject(s)
CpG Islands , DNA Methylation , Diabetes Mellitus, Type 2 , Humans , Diabetes Mellitus, Type 2/genetics , Female , Male , Genome-Wide Association Study , Genetic Predisposition to Disease , Middle Aged , Epigenesis, Genetic , Otx Transcription Factors/genetics , Otx Transcription Factors/metabolism , Diabetes Complications/genetics , Gene Expression ProfilingABSTRACT
Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.
Subject(s)
Biomarkers , Genome-Wide Association Study , Metabolomics , Female , Humans , Pregnancy , Acetone/blood , Acetone/metabolism , Biomarkers/blood , Biomarkers/metabolism , Cholestasis, Intrahepatic/blood , Cholestasis, Intrahepatic/genetics , Cholestasis, Intrahepatic/metabolism , Cohort Studies , Genome-Wide Association Study/methods , Hypertension/blood , Hypertension/genetics , Hypertension/metabolism , Lipoproteins/genetics , Lipoproteins/metabolism , Magnetic Resonance Spectroscopy , Mendelian Randomization Analysis , Metabolic Networks and Pathways/genetics , Phenotype , Polymorphism, Single Nucleotide/genetics , Pregnancy Complications/blood , Pregnancy Complications/genetics , Pregnancy Complications/metabolismABSTRACT
Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.
ABSTRACT
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.
Subject(s)
Fertility , Genome-Wide Association Study , Twinning, Dizygotic , Animals , Female , Humans , Pregnancy , Carrier Proteins/genetics , Fertility/genetics , Hormones , Proteins/genetics , United States , Zebrafish/geneticsABSTRACT
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.
Subject(s)
Galactose , Genome-Wide Association Study , Gene Regulatory Networks , Immunoglobulin G/genetics , Polysaccharides/metabolismABSTRACT
High-dimensional flow cytometry is the gold standard to study the human immune system in large cohorts. However, large sample sizes increase inter-experimental variation because of technical and experimental inaccuracies introduced by batch variability. Our high-throughput sample processing pipeline in combination with 28-color flow cytometry focuses on increased throughput (192 samples/experiment) and high reproducibility. We implemented quality control checkpoints to reduce technical and experimental variation. Finally, we integrated FlowSOM clustering to facilitate automated data analysis and demonstrate the reproducibility of our pipeline in a study with 3,357 samples. We reveal age-associated immune dynamics in 2,300 individuals, signified by decreasing T and B cell subsets with age. In addition, by combining genetic analyses, our approach revealed unique immune signatures associated with a single nucleotide polymorphism (SNP) that abrogates CD45 isoform splicing. In summary, we provide a versatile and reliable high-throughput, flow cytometry-based pipeline for immune discovery and exploration in large cohorts.
Subject(s)
B-Lymphocyte Subsets , Leukocytes , Humans , Immunophenotyping , Reproducibility of Results , Flow Cytometry/methodsABSTRACT
BACKGROUND: B vitamins such as folate (B9), B6, and B12 are key in one carbon metabolism, which generates methyl donors for DNA methylation. Several studies have linked differential methylation to self-reported intakes of folate and B12, but these estimates can be imprecise, while metabolomic biomarkers can offer an objective assessment of dietary intakes. We explored blood metabolomic biomarkers of folate and vitamins B6 and B12, to carry out epigenome-wide analyses across up to three European cohorts. Associations between self-reported habitual daily B vitamin intakes and 756 metabolites (Metabolon Inc.) were assessed in serum samples from 1064 UK participants from the TwinsUK cohort. The identified B vitamin metabolomic biomarkers were then used in epigenome-wide association tests with fasting blood DNA methylation levels at 430,768 sites from the Infinium HumanMethylation450 BeadChip in blood samples from 2182 European participants from the TwinsUK and KORA cohorts. Candidate signals were explored for metabolite associations with gene expression levels in a subset of the TwinsUK sample (n = 297). Metabolomic biomarker epigenetic associations were also compared with epigenetic associations of self-reported habitual B vitamin intakes in samples from 2294 European participants. RESULTS: Eighteen metabolites were associated with B vitamin intakes after correction for multiple testing (Bonferroni-adj. p < 0.05), of which 7 metabolites were available in both cohorts and tested for epigenome-wide association. Three metabolites - pipecolate (metabolomic biomarker of B6 and folate intakes), pyridoxate (marker of B6 and folate) and docosahexaenoate (DHA, marker of B6) - were associated with 10, 3 and 1 differentially methylated positions (DMPs), respectively. The strongest association was observed between DHA and DMP cg03440556 in the SCD gene (effect = 0.093 ± 0.016, p = 4.07E-09). Pyridoxate, a catabolic product of vitamin B6, was inversely associated with CpG methylation near the SLC1A5 gene promoter region (cg02711608 and cg22304262) and with SLC7A11 (cg06690548), but not with corresponding changes in gene expression levels. The self-reported intake of folate and vitamin B6 had consistent but non-significant associations with the epigenetic signals. CONCLUSION: Metabolomic biomarkers are a valuable approach to investigate the effects of dietary B vitamin intake on the human epigenome.
Subject(s)
Vitamin B Complex , Humans , Vitamin B 12 , Epigenome , DNA Methylation , Folic Acid , Vitamin B 6 , Biomarkers , Minor Histocompatibility Antigens , Amino Acid Transport System ASCABSTRACT
Short-chain fatty acids (SCFA) are involved in immune system and inflammatory responses. We comprehensively assessed the host genetic and gut microbial contribution to a panel of eight serum and stool SCFAs in two cohorts (TwinsUK, n = 2507; ZOE PREDICT-1, n = 328), examined their postprandial changes and explored their links with chronic and acute inflammatory responses in healthy individuals and trauma patients. We report low concordance between circulating and fecal SCFAs, significant postprandial changes in most circulating SCFAs, and a heritable genetic component (average h2: serum = 14%(SD = 14%); stool = 12%(SD = 6%)). Furthermore, we find that gut microbiome can accurately predict their fecal levels (AUC>0.71) while presenting weaker associations with serum. Finally, we report different correlation patterns with inflammatory markers depending on the type of inflammatory response (chronic or acute trauma). Our results illustrate the breadth of the physiological relevance of SCFAs on human inflammatory and metabolic responses highlighting the need for a deeper understanding of this important class of molecules.
Subject(s)
Gastrointestinal Microbiome , Humans , Gastrointestinal Microbiome/genetics , Fatty Acids, Volatile/metabolism , Feces , InflammationABSTRACT
BACKGROUND: A dysregulated postprandial metabolic response is a risk factor for chronic diseases, including type 2 diabetes mellitus (T2DM). The plasma protein N-glycome is implicated in both lipid metabolism and T2DM risk. Hence, we first investigate the relationship between the N-glycome and postprandial metabolism and then explore the mediatory role of the plasma N-glycome in the relationship between postprandial lipaemia and T2DM. METHODS: We included 995 individuals from the ZOE-PREDICT 1 study with plasma N-glycans measured by ultra-performance liquid chromatography at fasting and triglyceride, insulin, and glucose levels measured at fasting and following a mixed-meal challenge. Linear mixed models were used to investigate the associations between plasma protein N-glycosylation and metabolic response (fasting, postprandial (Cmax), or change from fasting). A mediation analysis was used to further explore the relationship of the N-glycome in the prediabetes (HbA1c = 39-47 mmol/mol (5.7-6.5%))-postprandial lipaemia association. RESULTS: We identified 36 out of 55 glycans significantly associated with postprandial triglycerides (Cmax ß ranging from -0.28 for low-branched glycans to 0.30 for GP26) after adjusting for covariates and multiple testing (padjusted < 0.05). N-glycome composition explained 12.6% of the variance in postprandial triglycerides not already explained by traditional risk factors. Twenty-seven glycans were also associated with postprandial glucose and 12 with postprandial insulin. Additionally, 3 of the postprandial triglyceride-associated glycans (GP9, GP11, and GP32) also correlate with prediabetes and partially mediate the relationship between prediabetes and postprandial triglycerides. CONCLUSIONS: This study provides a comprehensive overview of the interconnections between plasma protein N-glycosylation and postprandial responses, demonstrating the incremental predictive benefit of N-glycans. We also suggest a considerable proportion of the effect of prediabetes on postprandial triglycerides is mediated by some plasma N-glycans.
Subject(s)
Diabetes Mellitus, Type 2 , Hyperlipidemias , Prediabetic State , Humans , Blood Glucose/metabolism , Triglycerides , Insulin , Polysaccharides , Blood ProteinsABSTRACT
Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( SERPINA1, ZFP36L2 , and TLR10) signals contain predicted deleterious missense variants. Two loci, SOCS3 and HPN , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( FGG;FGB;FGA ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation. Key Points: Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.
ABSTRACT
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.
ABSTRACT
Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.
Subject(s)
Depression , Tandem Mass Spectrometry , Humans , Depression/metabolism , Diet , Metabolome/genetics , Vitamin A/metabolism , Hippurates , Metabolomics/methodsABSTRACT
Primary and secondary bile acids (BAs) influence metabolism and inflammation, and the gut microbiome modulates levels of BAs. We systematically explore the host genetic, gut microbial, and habitual dietary contribution to a panel of 19 serum and 15 stool BAs in two population-based cohorts (TwinsUK, n = 2,382; ZOE PREDICT-1, n = 327) and assess changes post-bariatric surgery and after nutritional interventions. We report that BAs have a moderately heritable genetic component, and the gut microbiome accurately predicts their levels in serum and stool. The secondary BA isoursodeoxycholate (isoUDCA) can be explained mostly by gut microbes (area under the receiver operating characteristic curve [AUC] = â¼80%) and associates with post-prandial lipemia and inflammation (GlycA). Furthermore, circulating isoUDCA decreases significantly 1 year after bariatric surgery (ß = -0.72, p = 1 × 10-5) and in response to fiber supplementation (ß = -0.37, p < 0.03) but not omega-3 supplementation. In healthy individuals, isoUDCA fasting levels correlate with pre-meal appetite (p < 1 × 10-4). Our findings indicate an important role for isoUDCA in lipid metabolism, appetite, and, potentially, cardiometabolic risk.
Subject(s)
Bariatric Surgery , Bile Acids and Salts , Humans , Appetite , Bariatric Surgery/adverse effects , Feces , InflammationABSTRACT
BACKGROUND: An association between blood pressure and aortic stiffness is well known, but ambiguity remains as to whether one precedes the other. This study aimed to investigate the association of aortic stiffness with contemporaneous versus historic blood pressure and direction of causality between aortic stiffening and hypertension in female twins. METHODS: Aortic stiffness, measured by carotid-femoral pulse wave velocity (PWV), and mean arterial pressure (MAP) was recorded in 2037 female TwinsUK participants (mean age: 62.4±9.7 years) at a single time point. A subset of 947 participants had repeat PWV and MAP measures (mean interval 5.5±1.7 years) with additional historic MAP (mean interval 6.6±3.3 years before baseline). RESULTS: Cross-sectional multivariable linear regression analysis confirmed PWV significantly associated with age and MAP. In longitudinal analysis, annual progression of PWV was not associated with historic MAP (standardized beta coefficient [ß]=-0.02, P=0.698), weakly associated with baseline MAP (ß=0.09, P=0.049) but strongly associated with progression (from baseline to most recent measurement) of MAP (ß= 0.26, P<0.001). Progression of MAP associated with both baseline and progression of PWV (ß=0.13, P=0.003 and ß=0.24, P<0.001, respectively). CONCLUSIONS: Progression of aortic stiffness associates more strongly with contemporaneous MAP compared with historic MAP. In contrast, progression of MAP is associated with prior arterial stiffness. These findings suggest a bidirectional relationship between arterial stiffness and blood pressure, and that lowering blood pressure may prevent a cycle of arterial stiffening and hypertension.
Subject(s)
Hypertension , Vascular Stiffness , Humans , Female , Middle Aged , Aged , Blood Pressure/physiology , Pulse Wave Analysis , Cross-Sectional Studies , Arterial Pressure/physiology , Vascular Stiffness/physiologyABSTRACT
Background: Ageing is a heterogenous process characterised by cellular and molecular hallmarks, including changes to haematopoietic stem cells and is a primary risk factor for chronic diseases. X chromosome inactivation (XCI) randomly transcriptionally silences either the maternal or paternal X in each cell of 46, XX females to balance the gene expression with 46, XY males. Age acquired XCI-skew describes the preferential selection of cells across a tissue resulting in an imbalance of XCI, which is particularly prevalent in blood tissues of ageing females, and yet its clinical consequences are unknown. Methods: We assayed XCI in 1575 females from the TwinsUK population cohort using DNA extracted from whole blood. We employed prospective, cross-sectional, and intra-twin study designs to characterise the relationship of XCI-skew with molecular and cellular measures of ageing, cardiovascular disease risk, and cancer diagnosis. Results: We demonstrate that XCI-skew is independent of traditional markers of biological ageing and is associated with a haematopoietic bias towards the myeloid lineage. Using an atherosclerotic cardiovascular disease risk score, which captures traditional risk factors, XCI-skew is associated with an increased cardiovascular disease risk both cross-sectionally and within XCI-skew discordant twin pairs. In a prospective 10 year follow-up study, XCI-skew is predictive of future cancer incidence. Conclusions: Our study demonstrates that age acquired XCI-skew captures changes to the haematopoietic stem cell population and has clinical potential as a unique biomarker of chronic disease risk. Funding: KSS acknowledges funding from the Medical Research Council [MR/M004422/1 and MR/R023131/1]. JTB acknowledges funding from the ESRC [ES/N000404/1]. MM acknowledges funding from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. TwinsUK is funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London.
Subject(s)
Cardiovascular Diseases , X Chromosome Inactivation , Female , Humans , Male , Cardiovascular Diseases/genetics , Cross-Sectional Studies , Follow-Up Studies , Outcome Assessment, Health Care , Prospective StudiesABSTRACT
Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.