Meta-analysis of genome-wide DNA methylation and integrative omics of age in human skeletal muscle

International audienceBackground: Knowledge of age-related DNA methylation changes in skeletal muscle is limited, yet this tissue is severely affected by ageing in humans.Methods: We conducted a large-scale epigenome-wide association study meta-analysis of age in human skeletal muscle from 10 studies (total n = 908 muscle methylomes from men and women aged 18-89 years old). We explored the genomic context of age-related DNA methylation changes in chromatin states, CpG islands, and transcription factor binding sites and performed gene set enrichment analysis. We then integrated the DNA methylation data with known transcriptomic and proteomic age-related changes in skeletal muscle. Finally, we updated our recently developed muscle epigenetic clock (https://bioconductor.org/packages/release/bioc/html/MEAT.html).Results: We identified 6710 differentially methylated regions at a stringent false discovery rate <0.005, spanning 6367 unique genes, many of which related to skeletal muscle structure and development. We found a strong increase in DNA methylation at Polycomb target genes and bivalent chromatin domains and a concomitant decrease in DNA methylation at enhancers. Most differentially methylated genes were not altered at the mRNA or protein level, but they were nonetheless strongly enriched for genes showing age-related differential mRNA and protein expression. After adding a substantial number of samples from five datasets (+371), the updated version of the muscle clock (MEAT 2.0, total n = 1053 samples) performed similarly to the original version of the muscle clock (median of 4.4 vs. 4.6 years in age prediction error), suggesting that the original version of the muscle clock was very accurate.Conclusions: We provide here the most comprehensive picture of DNA methylation ageing in human skeletal muscle and reveal widespread alterations of genes involved in skeletal muscle structure, development, and differentiation. We have made our results available as an open-access, user-friendly, web-based tool called MetaMeth (https://sarah-voisin.shinyapps.io/MetaMeth/)

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes

Genome Wide Analysis of Drug-Induced Torsades de Pointes: Lack of Common Variants with Large Effect Sizes

Marked prolongation of the QT interval on the electrocardiogram associated with the polymorphic ventricular tachycardia Torsades de Pointes is a serious adverse event during treatment with antiarrhythmic drugs and other culprit medications, and is a common cause for drug relabeling and withdrawal. Although clinical risk factors have been identified, the syndrome remains unpredictable in an individual patient. Here we used genome-wide association analysis to search for common predisposing genetic variants. Cases of drug-induced Torsades de Pointes (diTdP), treatment tolerant controls, and general population controls were ascertained across multiple sites using common definitions, and genotyped on the Illumina 610k or 1M-Duo BeadChips. Principal Components Analysis was used to select 216 Northwestern European diTdP cases and 771 ancestry-matched controls, including treatment-tolerant and general population subjects. With these sample sizes, there is 80% power to detect a variant at genome-wide significance with minor allele frequency of 10% and conferring an odds ratio of ≥2.7. Tests of association were carried out for each single nucleotide polymorphism (SNP) by logistic regression adjusting for gender and population structure. No SNP reached genome wide-significance; the variant with the lowest P value was rs2276314, a non-synonymous coding variant in C18orf21 (p  =  3×10(-7), odds ratio = 2, 95% confidence intervals: 1.5-2.6). The haplotype formed by rs2276314 and a second SNP, rs767531, was significantly more frequent in controls than cases (p  =  3×10(-9)). Expanding the number of controls and a gene-based analysis did not yield significant associations. This study argues that common genomic variants do not contribute importantly to risk for drug-induced Torsades de Pointes across multiple drugs

Coding Variation in ANGPTL4, LPL, and SVEP1 and the Risk of Coronary Disease.

BACKGROUND: The discovery of low-frequency coding variants affecting the risk of coronary artery disease has facilitated the identification of therapeutic targets. METHODS: Through DNA genotyping, we tested 54,003 coding-sequence variants covering 13,715 human genes in up to 72,868 patients with coronary artery disease and 120,770 controls who did not have coronary artery disease. Through DNA sequencing, we studied the effects of loss-of-function mutations in selected genes. RESULTS: We confirmed previously observed significant associations between coronary artery disease and low-frequency missense variants in the genes LPA and PCSK9. We also found significant associations between coronary artery disease and low-frequency missense variants in the genes SVEP1 (p.D2702G; minor-allele frequency, 3.60%; odds ratio for disease, 1.14; P=4.2×10(-10)) and ANGPTL4 (p.E40K; minor-allele frequency, 2.01%; odds ratio, 0.86; P=4.0×10(-8)), which encodes angiopoietin-like 4. Through sequencing of ANGPTL4, we identified 9 carriers of loss-of-function mutations among 6924 patients with myocardial infarction, as compared with 19 carriers among 6834 controls (odds ratio, 0.47; P=0.04); carriers of ANGPTL4 loss-of-function alleles had triglyceride levels that were 35% lower than the levels among persons who did not carry a loss-of-function allele (P=0.003). ANGPTL4 inhibits lipoprotein lipase; we therefore searched for mutations in LPL and identified a loss-of-function variant that was associated with an increased risk of coronary artery disease (p.D36N; minor-allele frequency, 1.9%; odds ratio, 1.13; P=2.0×10(-4)) and a gain-of-function variant that was associated with protection from coronary artery disease (p.S447*; minor-allele frequency, 9.9%; odds ratio, 0.94; P=2.5×10(-7)). CONCLUSIONS: We found that carriers of loss-of-function mutations in ANGPTL4 had triglyceride levels that were lower than those among noncarriers; these mutations were also associated with protection from coronary artery disease. (Funded by the National Institutes of Health and others.).Supported by a career development award from the National Heart, Lung, and Blood Institute, National Institutes of Health (NIH) (K08HL114642 to Dr. Stitziel) and by the Foundation for Barnes–Jewish Hospital. Dr. Peloso is supported by the National Heart, Lung, and Blood Institute of the NIH (award number K01HL125751). Dr. Kathiresan is supported by a Research Scholar award from the Massachusetts General Hospital, the Donovan Family Foundation, grants from the NIH (R01HL107816 and R01HL127564), a grant from Fondation Leducq, and an investigator-initiated grant from Merck. Dr. Merlini was supported by a grant from the Italian Ministry of Health (RFPS-2007-3-644382). Drs. Ardissino and Marziliano were supported by Regione Emilia Romagna Area 1 Grants. Drs. Farrall and Watkins acknowledge the support of the Wellcome Trust core award (090532/Z/09/Z), the British Heart Foundation (BHF) Centre of Research Excellence. Dr. Schick is supported in part by a grant from the National Cancer Institute (R25CA094880). Dr. Goel acknowledges EU FP7 & Wellcome Trust Institutional strategic support fund. Dr. Deloukas’s work forms part of the research themes contributing to the translational research portfolio of Barts Cardiovascular Biomedical Research Unit, which is supported and funded by the National Institute for Health Research (NIHR). Drs. Webb and Samani are funded by the British Heart Foundation, and Dr. Samani is an NIHR Senior Investigator. Dr. Masca was supported by the NIHR Leicester Cardiovascular Biomedical Research Unit (BRU), and this work forms part of the portfolio of research supported by the BRU. Dr. Won was supported by a postdoctoral award from the American Heart Association (15POST23280019). Dr. McCarthy is a Wellcome Trust Senior Investigator (098381) and an NIHR Senior Investigator. Dr. Danesh is a British Heart Foundation Professor, European Research Council Senior Investigator, and NIHR Senior Investigator. Drs. Erdmann, Webb, Samani, and Schunkert are supported by the FP7 European Union project CVgenes@ target (261123) and the Fondation Leducq (CADgenomics, 12CVD02). Drs. Erdmann and Schunkert are also supported by the German Federal Ministry of Education and Research e:Med program (e:AtheroSysMed and sysINFLAME), and Deutsche Forschungsgemeinschaft cluster of excellence “Inflammation at Interfaces” and SFB 1123. Dr. Kessler received a DZHK Rotation Grant. The analysis was funded, in part, by a Programme Grant from the BHF (RG/14/5/30893 to Dr. Deloukas). Additional funding is listed in the Supplementary Appendix.This is the author accepted manuscript. The final version is available from the Massachusetts Medical Society via http://dx.doi.org/10.1056/NEJMoa150765

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value OBFC1indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes

Validation of Human Telomere Length Multi-Ancestry Meta-Analysis Association Signals Identifies POP5 and KBTBD6 as Human Telomere Length Regulation Genes

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation

Phenome-wide association analysis of LDL-cholesterol lowering genetic variants in PCSK9

Abstract: Background: We characterised the phenotypic consequence of genetic variation at the PCSK9 locus and compared findings with recent trials of pharmacological inhibitors of PCSK9. Methods: Published and individual participant level data (300,000+ participants) were combined to construct a weighted PCSK9 gene-centric score (GS). Seventeen randomized placebo controlled PCSK9 inhibitor trials were included, providing data on 79,578 participants. Results were scaled to a one mmol/L lower LDL-C concentration. Results: The PCSK9 GS (comprising 4 SNPs) associations with plasma lipid and apolipoprotein levels were consistent in direction with treatment effects. The GS odds ratio (OR) for myocardial infarction (MI) was 0.53 (95% CI 0.42; 0.68), compared to a PCSK9 inhibitor effect of 0.90 (95% CI 0.86; 0.93). For ischemic stroke ORs were 0.84 (95% CI 0.57; 1.22) for the GS, compared to 0.85 (95% CI 0.78; 0.93) in the drug trials. ORs with type 2 diabetes mellitus (T2DM) were 1.29 (95% CI 1.11; 1.50) for the GS, as compared to 1.00 (95% CI 0.96; 1.04) for incident T2DM in PCSK9 inhibitor trials. No genetic associations were observed for cancer, heart failure, atrial fibrillation, chronic obstructive pulmonary disease, or Alzheimer’s disease – outcomes for which large-scale trial data were unavailable. Conclusions: Genetic variation at the PCSK9 locus recapitulates the effects of therapeutic inhibition of PCSK9 on major blood lipid fractions and MI. While indicating an increased risk of T2DM, no other possible safety concerns were shown; although precision was moderate