1
|
Lu Y, Xu K, Maydanchik N, Kang B, Pierce BL, Yang F, Chen LS. An integrative multi-context Mendelian randomization method for identifying risk genes across human tissues. Am J Hum Genet 2024; 111:1736-1749. [PMID: 39053459 PMCID: PMC11339623 DOI: 10.1016/j.ajhg.2024.06.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2024] [Revised: 06/11/2024] [Accepted: 06/24/2024] [Indexed: 07/27/2024] Open
Abstract
Mendelian randomization (MR) provides valuable assessments of the causal effect of exposure on outcome, yet the application of conventional MR methods for mapping risk genes encounters new challenges. One of the issues is the limited availability of expression quantitative trait loci (eQTLs) as instrumental variables (IVs), hampering the estimation of sparse causal effects. Additionally, the often context- or tissue-specific eQTL effects challenge the MR assumption of consistent IV effects across eQTL and GWAS data. To address these challenges, we propose a multi-context multivariable integrative MR framework, mintMR, for mapping expression and molecular traits as joint exposures. It models the effects of molecular exposures across multiple tissues in each gene region, while simultaneously estimating across multiple gene regions. It uses eQTLs with consistent effects across more than one tissue type as IVs, improving IV consistency. A major innovation of mintMR involves employing multi-view learning methods to collectively model latent indicators of disease relevance across multiple tissues, molecular traits, and gene regions. The multi-view learning captures the major patterns of disease relevance and uses these patterns to update the estimated tissue relevance probabilities. The proposed mintMR iterates between performing a multi-tissue MR for each gene region and joint learning the disease-relevant tissue probabilities across gene regions, improving the estimation of sparse effects across genes. We apply mintMR to evaluate the causal effects of gene expression and DNA methylation for 35 complex traits using multi-tissue QTLs as IVs. The proposed mintMR controls genome-wide inflation and offers insights into disease mechanisms.
Collapse
Affiliation(s)
- Yihao Lu
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
| | - Ke Xu
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, USA
| | - Nathaniel Maydanchik
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
| | - Bowei Kang
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
| | - Brandon L Pierce
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
| | - Fan Yang
- Yau Mathematical Sciences Center, Tsinghua University, Beijing, China; Yanqi Lake Beijing Institute of Mathematical Sciences and Applications, Beijing, China.
| | - Lin S Chen
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA.
| |
Collapse
|
2
|
Kazan HH, Kasakolu A, Koncagul S, Ergun MA, John G, Sultanov RI, Zhelankin AV, Semenova EA, Yusupov RA, Kulemin NA, Larin AK, Generozov EV, Bulgay C, Ahmetov II. Association analysis of indel variants and gene expression identifies MDM4 as a novel locus for skeletal muscle hypertrophy and power athlete status. Exp Physiol 2024. [PMID: 39041487 DOI: 10.1113/ep091992] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Accepted: 06/26/2024] [Indexed: 07/24/2024]
Abstract
Insertions and deletions (indels) are the second most common type of variation in the human genome. However, limited data on their associations with exercise-related phenotypes have been documented. The aim of the present study was to examine the association between 18,370 indel variants and power athlete status, followed by additional studies in 357,246 individuals. In the discovery phase, the D allele of the MDM4 gene rs35493922 I/D polymorphism was over-represented in power athletes compared with control subjects (P = 7.8 × 10-9) and endurance athletes (P = 0.0012). These findings were replicated in independent cohorts, showing a higher D allele frequency in power athletes compared with control subjects (P = 0.016) and endurance athletes (P = 0.031). Furthermore, the D allele was positively associated (P = 0.0013) with greater fat-free mass in the UK Biobank. MDM4 encodes a protein that inhibits the activity of p53, which induces muscle fibre atrophy. Accordingly, we found that MDM4 expression was significantly higher in the vastus lateralis of power athletes compared with endurance athletes (P = 0.0009) and was positively correlated with the percentage of fast-twitch muscle fibres (P = 0.0062) and the relative area occupied by fast-twitch muscle fibres (P = 0.0086). The association between MDM4 gene expression and an increased proportion of fast-twitch muscle fibres was confirmed in two additional cohorts. Finally, we found that the MDM4 DD genotype was associated with increased MDM4 gene expression in vastus lateralis and greater cross-sectional area of fast-twitch muscle fibres. In conclusion, MDM4 is suggested to be a potential regulator of muscle fibre specification and size, with its indel variant being associated with power athlete status. HIGHLIGHTS: What is the central question of this study? Which indel variants are functional and associated with sport- and exercise-related traits? What is the main finding and its importance? Out of 18,370 tested indels, the MDM4 gene rs35493922 I/D polymorphism was found to be the functional variant (affecting gene expression) and the most significant, with the deletion allele showing associations with power athlete status, fat-free mass and cross-sectional area of fast-twitch muscle fibres. Furthermore, the expression of MDM4 was positively correlated with the percentage of fast-twitch muscle fibres and the relative area occupied by fast-twitch muscle fibres.
Collapse
Affiliation(s)
- Hasan H Kazan
- Department of Medical Biology, Gulhane Faculty of Medicine, University of Health Sciences, Ankara, Türkiye
| | - Anıl Kasakolu
- Graduate School of Natural and Applied Sciences, Ankara University, Ankara, Türkiye
| | - Seyrani Koncagul
- Graduate School of Natural and Applied Sciences, Ankara University, Ankara, Türkiye
| | - Mehmet A Ergun
- Department of Medical Genetics, Faculty of Medicine, Gazi University, Ankara, Türkiye
| | - George John
- Transform Specialist Medical Centre, Dubai, UAE
| | - Rinat I Sultanov
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
| | - Andrey V Zhelankin
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
| | - Ekaterina A Semenova
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
- Research Institute of Physical Culture and Sport, Volga Region State University of Physical Culture, Sport and Tourism, Kazan, Russia
| | - Rinat A Yusupov
- Department of Physical Culture and Sport, Kazan National Research Technical University Named after A.N. Tupolev-KAI, Kazan, Russia
| | - Nikolay A Kulemin
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
| | - Andrey K Larin
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
| | - Edward V Generozov
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
| | - Celal Bulgay
- Sports Science Faculty, Bingol University, Bingol, Türkiye
| | - Ildus I Ahmetov
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
- Sports Genetics Laboratory, St Petersburg Research Institute of Physical Culture, St Petersburg, Russia
- Laboratory of Genetics of Aging and Longevity, Kazan State Medical University, Kazan, Russia
- Research Institute for Sport and Exercise Sciences, Liverpool John Moores University, Liverpool, UK
| |
Collapse
|
3
|
Grenko CM, Taylor HJ, Bonnycastle LL, Xue D, Lee BN, Weiss Z, Yan T, Swift AJ, Mansell EC, Lee A, Robertson CC, Narisu N, Erdos MR, Chen S, Collins FS, Taylor DL. Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 h. Diabetologia 2024:10.1007/s00125-024-06214-4. [PMID: 38967666 DOI: 10.1007/s00125-024-06214-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/04/2024] [Accepted: 05/08/2024] [Indexed: 07/06/2024]
Abstract
AIMS/HYPOTHESIS Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycaemia, beta cell glucotoxicity and subsequently type 2 diabetes. In this study, we explored the effects of in vitro hyperglycaemic conditions on human pancreatic islet gene expression across 24 h in six pancreatic cell types: alpha; beta; gamma; delta; ductal; and acinar. We hypothesised that genes associated with hyperglycaemic conditions may be relevant to the onset and progression of diabetes. METHODS We exposed human pancreatic islets from two donors to low (2.8 mmol/l) and high (15.0 mmol/l) glucose concentrations over 24 h in vitro. To assess the transcriptome, we performed single-cell RNA-seq (scRNA-seq) at seven time points. We modelled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Additionally, we integrated genomic features and genetic summary statistics to nominate candidate effector genes. For three of these genes, we functionally characterised the effect on insulin production and secretion using CRISPR interference to knock down gene expression in EndoC-βH1 cells, followed by a glucose-stimulated insulin secretion assay. RESULTS In the discrete time models, we identified 1344 genes associated with time and 668 genes associated with glucose exposure across all cell types and time points. In the continuous time models, we identified 1311 genes associated with time, 345 genes associated with glucose exposure and 418 genes associated with interaction effects between time and glucose across all cell types. By integrating these expression profiles with summary statistics from genetic association studies, we identified 2449 candidate effector genes for type 2 diabetes, HbA1c, random blood glucose and fasting blood glucose. Of these candidate effector genes, we showed that three (ERO1B, HNRNPA2B1 and RHOBTB3) exhibited an effect on glucose-stimulated insulin production and secretion in EndoC-βH1 cells. CONCLUSIONS/INTERPRETATION The findings of our study provide an in-depth characterisation of the 24 h transcriptomic response of human pancreatic islets to glucose exposure at a single-cell resolution. By integrating differentially expressed genes with genetic signals for type 2 diabetes and glucose-related traits, we provide insights into the molecular mechanisms underlying glucose homeostasis. Finally, we provide functional evidence to support the role of three candidate effector genes in insulin secretion and production. DATA AVAILABILITY The scRNA-seq data from the 24 h glucose exposure experiment performed in this study are available in the database of Genotypes and Phenotypes (dbGap; https://www.ncbi.nlm.nih.gov/gap/ ) with accession no. phs001188.v3.p1. Study metadata and summary statistics for the differential expression, gene set enrichment and candidate effector gene prediction analyses are available in the Zenodo data repository ( https://zenodo.org/ ) under accession number 11123248. The code used in this study is publicly available at https://github.com/CollinsLabBioComp/publication-islet_glucose_timecourse .
Collapse
Affiliation(s)
- Caleb M Grenko
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
- Graduate School of Biomedical Sciences, Mayo Clinic, Rochester, MN, USA
| | - Henry J Taylor
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA.
- British Heart Foundation Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK.
- Heart and Lung Research Institute, University of Cambridge, Cambridge, UK.
| | - Lori L Bonnycastle
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Dongxiang Xue
- Department of Surgery, Weill Cornell Medicine, New York, NY, USA
- Center for Genomic Health, Weill Cornell Medicine, New York, NY, USA
| | - Brian N Lee
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Zoe Weiss
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Tingfen Yan
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Amy J Swift
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Erin C Mansell
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Angela Lee
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Catherine C Robertson
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Narisu Narisu
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Michael R Erdos
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Shuibing Chen
- Department of Surgery, Weill Cornell Medicine, New York, NY, USA
- Center for Genomic Health, Weill Cornell Medicine, New York, NY, USA
| | - Francis S Collins
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA.
| | - D Leland Taylor
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| |
Collapse
|
4
|
Nair VD, Pincas H, Smith GR, Zaslavsky E, Ge Y, Amper MAS, Vasoya M, Chikina M, Sun Y, Raja AN, Mao W, Gay NR, Esser KA, Smith KS, Zhao B, Wiel L, Singh A, Lindholm ME, Amar D, Montgomery S, Snyder MP, Walsh MJ, Sealfon SC. Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures. CELL GENOMICS 2024; 4:100421. [PMID: 38697122 PMCID: PMC11228891 DOI: 10.1016/j.xgen.2023.100421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 07/07/2023] [Accepted: 09/12/2023] [Indexed: 05/04/2024]
Abstract
Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.
Collapse
Affiliation(s)
- Venugopalan D Nair
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA.
| | - Hanna Pincas
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Gregory R Smith
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Elena Zaslavsky
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Yongchao Ge
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Mary Anne S Amper
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Mital Vasoya
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Maria Chikina
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Yifei Sun
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | | | - Weiguang Mao
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Nicole R Gay
- Department of Genetics, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Karyn A Esser
- Department of Physiology and Aging, University of Florida, Gainesville, FL 32610, USA
| | - Kevin S Smith
- Departments of Pathology and Genetics, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Bingqing Zhao
- Department of Genetics, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Laurens Wiel
- Department of Medicine, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Aditya Singh
- Department of Medicine, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Malene E Lindholm
- Department of Medicine, Stanford School of Medicine, Stanford, CA 94305, USA
| | - David Amar
- Department of Medicine, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Stephen Montgomery
- Departments of Pathology and Genetics, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Michael P Snyder
- Department of Genetics, Stanford School of Medicine, Stanford, CA 94305, USA
| | - Martin J Walsh
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Stuart C Sealfon
- Department of Neurology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA.
| |
Collapse
|
5
|
Bittel AJ, Chen YW. DNA Methylation in the Adaptive Response to Exercise. Sports Med 2024; 54:1419-1458. [PMID: 38561436 DOI: 10.1007/s40279-024-02011-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/23/2024] [Indexed: 04/04/2024]
Abstract
Emerging evidence published over the past decade has highlighted the role of DNA methylation in skeletal muscle function and health, including as an epigenetic transducer of the adaptive response to exercise. In this review, we aim to synthesize the latest findings in this field to highlight: (1) the shifting understanding of the genomic localization of altered DNA methylation in response to acute and chronic aerobic and resistance exercise in skeletal muscle (e.g., promoter, gene bodies, enhancers, intergenic regions, un-annotated regions, and genome-wide methylation); (2) how these global/regional methylation changes relate to transcriptional activity following exercise; and (3) the factors (e.g., individual demographic or genetic features, dietary, training history, exercise parameters, local epigenetic characteristics, circulating hormones) demonstrated to alter both the pattern of DNA methylation after exercise, and the relationship between DNA methylation and gene expression. Finally, we discuss the changes in non-CpG methylation and 5-hydroxymethylation after exercise, as well as the importance of emerging single-cell analyses to future studies-areas of increasing focus in the field of epigenetics. We anticipate that this review will help generate a framework for clinicians and researchers to begin developing and testing exercise interventions designed to generate targeted changes in DNA methylation as part of a personalized exercise regimen.
Collapse
Affiliation(s)
- Adam J Bittel
- Research Center for Genetic Medicine, Children's National Hospital, 111 Michigan Ave NW, Washington, DC, 20010, USA.
| | - Yi-Wen Chen
- Research Center for Genetic Medicine, Children's National Hospital, 111 Michigan Ave NW, Washington, DC, 20010, USA
- Department of Genomics and Precision Medicine, The George Washington University School of Medicine and Health Science, 111 Michigan Ave NW, Washington, DC, 20010, USA
- Department of Integrative Systems Biology, Institute for Biomedical Sciences, The George Washington University, 2121 I St NW, Washington, DC, 20052, USA
| |
Collapse
|
6
|
Qin Q, Zhou Y, Guo J, Chen Q, Tang W, Li Y, You J, Li Q. Conserved methylation signatures associate with the tumor immune microenvironment and immunotherapy response. Genome Med 2024; 16:47. [PMID: 38566132 PMCID: PMC10985907 DOI: 10.1186/s13073-024-01318-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Accepted: 03/20/2024] [Indexed: 04/04/2024] Open
Abstract
BACKGROUND Aberrant DNA methylation is a major characteristic of cancer genomes. It remains unclear which biological processes determine epigenetic reprogramming and how these processes influence the variants in the cancer methylome, which can further impact cancer phenotypes. METHODS We performed pairwise permutations of 381,900 loci in 569 paired DNA methylation profiles of cancer tissue and matched normal tissue from The Cancer Genome Atlas (TCGA) and defined conserved differentially methylated positions (DMPs) based on the resulting null distribution. Then, we derived independent methylation signatures from 2,465 cancer-only methylation profiles from the TCGA and 241 cell line-based methylation profiles from the Genomics of Drug Sensitivity in Cancer (GDSC) cohort using nonnegative matrix factorization (NMF). We correlated DNA methylation signatures with various clinical and biological features, including age, survival, cancer stage, tumor immune microenvironment factors, and immunotherapy response. We inferred the determinant genes of these methylation signatures by integrating genomic and transcriptomic data and evaluated the impact of these signatures on cancer phenotypes in independent bulk and single-cell RNA/methylome cohorts. RESULTS We identified 7,364 differentially methylated positions (2,969 Hyper-DMPs and 4,395 Hypo-DMPs) in nine cancer types from the TCGA. We subsequently retrieved three highly conserved, independent methylation signatures (Hyper-MS1, Hypo-MS1, and Hypo-MS4) from cancer tissues and cell lines based on these Hyper and Hypo-DMPs. Our data suggested that Hypo-MS4 activity predicts poor survival and is associated with immunotherapy response and distant tumor metastasis, and Hypo-MS4 activity is related to TP53 mutation and FOXA1 binding specificity. In addition, we demonstrated a correlation between the activities of Hypo-MS4 in cancer cells and the fractions of regulatory CD4 + T cells with the expression levels of immunological genes in the tumor immune microenvironment. CONCLUSIONS Our findings demonstrated that the methylation signatures of distinct biological processes are associated with immune activity in the cancer microenvironment and predict immunotherapy response.
Collapse
Affiliation(s)
- Qingqing Qin
- Department of Hematology, The First Affiliated Hospital of Xiamen University and Institute of Hematology, School of Medicine, Xiamen University, Xiamen, 361003, China
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China
- Department of Pediatrics, Women and Children's Hospital, School of Medicine, Xiamen University, Xiamen, 361003, China
| | - Ying Zhou
- Department of Hematology, The First Affiliated Hospital of Xiamen University and Institute of Hematology, School of Medicine, Xiamen University, Xiamen, 361003, China
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China
- Department of Pediatrics, Women and Children's Hospital, School of Medicine, Xiamen University, Xiamen, 361003, China
| | - Jintao Guo
- Department of Hematology, The First Affiliated Hospital of Xiamen University and Institute of Hematology, School of Medicine, Xiamen University, Xiamen, 361003, China
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China
- Department of Pediatrics, Women and Children's Hospital, School of Medicine, Xiamen University, Xiamen, 361003, China
| | - Qinwei Chen
- Department of Hematology, The First Affiliated Hospital of Xiamen University and Institute of Hematology, School of Medicine, Xiamen University, Xiamen, 361003, China
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China
| | - Weiwei Tang
- Department of Medical Oncology, School of Medicine, The First Affiliated Hospital of Xiamen University and Institute of Hematology, Xiamen University, Xiamen, 361003, China
- Xiamen Key Laboratory of Antitumor Drug Transformation Research, The First Affiliated Hospital of Xiamen University, School of Medicine, Xiamen University, The School of Clinical Medicine of Fujian, Medical University, Xiamen, 361003, China
| | - Yuchen Li
- Department of Hematology, The First Affiliated Hospital of Xiamen University and Institute of Hematology, School of Medicine, Xiamen University, Xiamen, 361003, China
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China
- Department of Pediatrics, Women and Children's Hospital, School of Medicine, Xiamen University, Xiamen, 361003, China
| | - Jun You
- Department of Gastrointestinal Oncology Surgery, The First Affiliated Hospital of Xiamen University, Cancer Center, Xiamen, 361003, China
| | - Qiyuan Li
- Department of Hematology, The First Affiliated Hospital of Xiamen University and Institute of Hematology, School of Medicine, Xiamen University, Xiamen, 361003, China.
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China.
- Department of Pediatrics, Women and Children's Hospital, School of Medicine, Xiamen University, Xiamen, 361003, China.
| |
Collapse
|
7
|
Lu Y, Oliva M, Pierce BL, Liu J, Chen LS. Integrative cross-omics and cross-context analysis elucidates molecular links underlying genetic effects on complex traits. Nat Commun 2024; 15:2383. [PMID: 38493154 PMCID: PMC10944527 DOI: 10.1038/s41467-024-46675-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 03/06/2024] [Indexed: 03/18/2024] Open
Abstract
Genetic effects on functionally related 'omic' traits often co-occur in relevant cellular contexts, such as tissues. Motivated by the multi-tissue methylation quantitative trait loci (mQTLs) and expression QTLs (eQTLs) analysis, we propose X-ING (Cross-INtegrative Genomics) for cross-omics and cross-context integrative analysis. X-ING takes as input multiple matrices of association statistics, each obtained from different omics data types across multiple cellular contexts. It models the latent binary association status of each statistic, captures the major association patterns among omics data types and contexts, and outputs the posterior mean and probability for each input statistic. X-ING enables the integration of effects from different omics data with varying effect distributions. In the multi-tissue cis-association analysis, X-ING shows improved detection and replication of mQTLs by integrating eQTL maps. In the trans-association analysis, X-ING reveals an enrichment of trans-associations in many disease/trait-relevant tissues.
Collapse
Affiliation(s)
- Yihao Lu
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
| | - Meritxell Oliva
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
- Genomics Research Center, AbbVie, North Chicago, IL, USA
| | - Brandon L Pierce
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA
| | - Jin Liu
- School of Data Science, The Chinese University of Hong Kong-Shenzhen, Shenzhen, China.
| | - Lin S Chen
- Department of Public Health Sciences, The University of Chicago, Chicago, IL, USA.
| |
Collapse
|
8
|
Lu Y, Xu K, Kang B, Pierce BL, Yang F, Chen LS. An integrative multi-context Mendelian randomization method for identifying risk genes across human tissues. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2024:2024.03.04.24303731. [PMID: 38496462 PMCID: PMC10942526 DOI: 10.1101/2024.03.04.24303731] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
Mendelian randomization (MR) provides valuable assessments of the causal effect of exposure on outcome, yet the application of conventional MR methods for mapping risk genes encounters new challenges. One of the issues is the limited availability of expression quantitative trait loci (eQTLs) as instrumental variables (IVs), hampering the estimation of sparse causal effects. Additionally, the often context/tissue-specific eQTL effects challenge the MR assumption of consistent IV effects across eQTL and GWAS data. To address these challenges, we propose a multi-context multivariable integrative MR framework, mintMR, for mapping expression and molecular traits as joint exposures. It models the effects of molecular exposures across multiple tissues in each gene region, while simultaneously estimating across multiple gene regions. It uses eQTLs with consistent effects across more than one tissue type as IVs, improving IV consistency. A major innovation of mintMR involves employing multi-view learning methods to collectively model latent indicators of disease relevance across multiple tissues, molecular traits, and gene regions. The multi-view learning captures the major patterns of disease-relevance and uses these patterns to update the estimated tissue relevance probabilities. The proposed mintMR iterates between performing a multi-tissue MR for each gene region and joint learning the disease-relevant tissue probabilities across gene regions, improving the estimation of sparse effects across genes. We apply mintMR to evaluate the causal effects of gene expression and DNA methylation for 35 complex traits using multi-tissue QTLs as IVs. The proposed mintMR controls genome-wide inflation and offers new insights into disease mechanisms.
Collapse
|
9
|
Advani J, Mehta PA, Hamel AR, Mehrotra S, Kiel C, Strunz T, Corso-Díaz X, Kwicklis M, van Asten F, Ratnapriya R, Chew EY, Hernandez DG, Montezuma SR, Ferrington DA, Weber BHF, Segrè AV, Swaroop A. QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration. Nat Commun 2024; 15:1972. [PMID: 38438351 PMCID: PMC10912779 DOI: 10.1038/s41467-024-46063-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 02/12/2024] [Indexed: 03/06/2024] Open
Abstract
DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.
Collapse
Affiliation(s)
- Jayshree Advani
- Neurobiology, Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD, USA
| | - Puja A Mehta
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, MA, USA
- Department of Ophthalmology, Harvard Medical School, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Andrew R Hamel
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, MA, USA
- Department of Ophthalmology, Harvard Medical School, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Sudeep Mehrotra
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, MA, USA
- Department of Ophthalmology, Harvard Medical School, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Christina Kiel
- Institute of Human Genetics, University of Regensburg, Regensburg, Germany
| | - Tobias Strunz
- Institute of Human Genetics, University of Regensburg, Regensburg, Germany
| | - Ximena Corso-Díaz
- Neurobiology, Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD, USA
| | - Madeline Kwicklis
- Neurobiology, Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD, USA
| | - Freekje van Asten
- Neurobiology, Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD, USA
| | - Rinki Ratnapriya
- Neurobiology, Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD, USA
| | - Emily Y Chew
- Division of Epidemiology and Clinical Applications, Clinical Trials Branch, National Eye Institute, National Institutes of Health, Bethesda, MD, USA
| | - Dena G Hernandez
- Laboratory of Neurogenetics, National Institute of Aging, National Institutes of Health, Bethesda, MD, USA
| | - Sandra R Montezuma
- Department of Ophthalmology and Visual Neurosciences, University of Minnesota, Minneapolis, MN, USA
| | - Deborah A Ferrington
- Department of Ophthalmology and Visual Neurosciences, University of Minnesota, Minneapolis, MN, USA
- Doheny Eye Institute, Pasadena, CA, USA
| | - Bernhard H F Weber
- Institute of Human Genetics, University of Regensburg, Regensburg, Germany
- Institute of Clinical Human Genetics, University Hospital Regensburg, Regensburg, Germany
| | - Ayellet V Segrè
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, MA, USA.
- Department of Ophthalmology, Harvard Medical School, Boston, MA, USA.
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
| | - Anand Swaroop
- Neurobiology, Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD, USA.
| |
Collapse
|
10
|
Pahl MC, Liu L, Pippin JA, Wagley Y, Boehm K, Hankenson KD, Wells AD, Yang W, Grant SFA. Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements. EBioMedicine 2024; 101:105038. [PMID: 38417377 PMCID: PMC10909706 DOI: 10.1016/j.ebiom.2024.105038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Revised: 02/13/2024] [Accepted: 02/14/2024] [Indexed: 03/01/2024] Open
Abstract
BACKGROUND Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms. METHODS We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts. FINDINGS We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro. INTERPRETATION Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting. FUNDING This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.
Collapse
Affiliation(s)
- Matthew C Pahl
- Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Division of Human Genetics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Lin Liu
- Department of Medicine, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Cell and Developmental Biology, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Institute for Diabetes, Obesity & Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia PA19104, USA
| | - James A Pippin
- Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Division of Human Genetics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Yadav Wagley
- Orthopaedic Research Laboratories, Department of Orthopaedic Surgery, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Keith Boehm
- Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Division of Human Genetics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Kurt D Hankenson
- Orthopaedic Research Laboratories, Department of Orthopaedic Surgery, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Andrew D Wells
- Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, 3615 Civic Center Boulevard, Philadelphia, PA, USA; Institute for Immunology, Perelman School of Medicine, University of Pennsylvania, 3615 Civic Center Boulevard, Philadelphia, PA, USA
| | - Wenli Yang
- Department of Medicine, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Cell and Developmental Biology, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Institute for Diabetes, Obesity & Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia PA19104, USA.
| | - Struan F A Grant
- Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Division of Human Genetics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Department of Pediatrics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA.
| |
Collapse
|
11
|
Das SS, Das SK. Common and ethnic-specific derangements in skeletal muscle transcriptome associated with obesity. Int J Obes (Lond) 2024; 48:330-338. [PMID: 37993634 DOI: 10.1038/s41366-023-01417-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Revised: 10/25/2023] [Accepted: 11/07/2023] [Indexed: 11/24/2023]
Abstract
BACKGROUND Obesity is a common disease with a higher prevalence among African Americans. Obesity alters cellular function in many tissues, including skeletal muscle, and is a risk factor for many life-threatening diseases, including cardiovascular disease and diabetes. The similarities and differences in molecular mechanisms that may explain ethnic disparities in obesity between African and European ancestry individuals have not been studied. METHODS In this study, data from transcriptome-wide analyses on skeletal muscle tissues from well-powered human cohorts were used to compare genes and biological pathways affected by obesity in European and African ancestry populations. Data on obesity-induced differentially expressed transcripts and GWAS-identified SNPs were integrated to prioritize target genes for obesity-associated genetic variants. RESULTS Linear regression analysis in the FUSION (European, N = 301) and AAGMEx (African American, N = 256) cohorts identified a total of 2569 body mass index (BMI)-associated transcripts (q < 0.05), of which 970 genes (at p < 0.05) are associated in both cohorts, and the majority showed the same direction of effect on BMI. Biological pathway analyses, including over-representation and gene-set enrichment analyses, identified enrichment of protein synthesis pathways (e.g., ribosomal function) and the ceramide signaling pathway in both cohorts among BMI-associated down- and up-regulated transcripts, respectively. A comparison using the IPA-tool suggested the activation of inflammation pathways only in Europeans with obesity. Interestingly, these analyses suggested repression of the mitochondrial oxidative phosphorylation pathway in Europeans but showed its activation in African Americans. Integration of SNP-to-Gene analyses-predicted target genes for obesity-associated genetic variants (GWAS-identified SNPs) and BMI-associated transcripts suggested that these SNPs might cause obesity by altering the expression of 316 critical target genes (e.g., GRB14) in the muscle. CONCLUSIONS This study provides a replication of obesity-associated transcripts and biological pathways in skeletal muscle across ethnicities, but also identifies obesity-associated processes unique in either African or European ancestry populations.
Collapse
Affiliation(s)
- Sreejon S Das
- The School of Biotechnology at Atkins, Atkins Academic and Technology High, Winston-Salem, NC, 27101, USA
| | - Swapan K Das
- Department of Internal Medicine, Section of Endocrinology and Metabolism, Wake Forest University School of Medicine, Winston-Salem, NC, 27157, USA.
| |
Collapse
|
12
|
Zhang Z, Wang S, Jiang L, Wei J, Lu C, Li S, Diao Y, Fang Z, He S, Tan T, Yang Y, Zou K, Shi J, Lin J, Chen L, Bao C, Fei J, Fang H. Priority index for critical Covid-19 identifies clinically actionable targets and drugs. Commun Biol 2024; 7:189. [PMID: 38366110 PMCID: PMC10873402 DOI: 10.1038/s42003-024-05897-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2023] [Accepted: 02/07/2024] [Indexed: 02/18/2024] Open
Abstract
While genome-wide studies have identified genomic loci in hosts associated with life-threatening Covid-19 (critical Covid-19), the challenge of resolving these loci hinders further identification of clinically actionable targets and drugs. Building upon our previous success, we here present a priority index solution designed to address this challenge, generating the target and drug resource that consists of two indexes: the target index and the drug index. The primary purpose of the target index is to identify clinically actionable targets by prioritising genes associated with Covid-19. We illustrate the validity of the target index by demonstrating its ability to identify pre-existing Covid-19 phase-III drug targets, with the majority of these targets being found at the leading prioritisation (leading targets). These leading targets have their evolutionary origins in Amniota ('four-leg vertebrates') and are predominantly involved in cytokine-cytokine receptor interactions and JAK-STAT signaling. The drug index highlights opportunities for repurposing clinically approved JAK-STAT inhibitors, either individually or in combination. This proposed strategic focus on the JAK-STAT pathway is supported by the active pursuit of therapeutic agents targeting this pathway in ongoing phase-II/III clinical trials for Covid-19.
Collapse
Affiliation(s)
- Zhiqiang Zhang
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
- School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Shan Wang
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Lulu Jiang
- Translational Health Sciences, University of Bristol, Bristol, BS1 3NY, UK
| | - Jianwen Wei
- Network and Information Center, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Chang Lu
- MRC London Institute of Medical Sciences, Imperial College London, London, W12 0HS, UK
| | - Shengli Li
- Precision Research Center for Refractory Diseases, Institute for Clinical Research, Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 201620, China
| | - Yizhu Diao
- College of Finance and Statistics, Hunan University, Changsha, 410079, Hunan, China
| | - Zhongcheng Fang
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
- School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Shuo He
- College of Health Science and Technology, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Tingting Tan
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Yisheng Yang
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Kexin Zou
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
- School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jiantao Shi
- Key Laboratory of RNA Science and Engineering, Shanghai Institute of Biochemistry and Cell Biology, Center for Excellence in Molecular Cell Science, Chinese Academy of Sciences, Shanghai, 200031, China
| | - James Lin
- Network and Information Center, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Liye Chen
- Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Oxford, OX3 7LD, UK.
| | - Chaohui Bao
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China.
- Department of General Surgery, Ruijin Hospital Luwan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, 200020, China.
| | - Jian Fei
- Department of General Surgery, Ruijin Hospital Luwan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, 200020, China.
- Department of General Surgery, Pancreatic Disease Center, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China.
| | - Hai Fang
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China.
| |
Collapse
|
13
|
Kober KM, Berger L, Roy R, Olshen A. Torch-eCpG: a fast and scalable eQTM mapper for thousands of molecular phenotypes with graphical processing units. BMC Bioinformatics 2024; 25:71. [PMID: 38355413 PMCID: PMC10867984 DOI: 10.1186/s12859-024-05670-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Accepted: 01/22/2024] [Indexed: 02/16/2024] Open
Abstract
BACKGROUND Gene expression may be regulated by the DNA methylation of regulatory elements in cis, distal, and trans regions. One method to evaluate the relationship between DNA methylation and gene expression is the mapping of expression quantitative trait methylation (eQTM) loci (also called expression associated CpG loci, eCpG). However, no open-source tools are available to provide eQTM mapping. In addition, eQTM mapping can involve a large number of comparisons which may prevent the analyses due to limitations of computational resources. Here, we describe Torch-eCpG, an open-source tool to perform eQTM mapping that includes an optimized implementation that can use the graphical processing unit (GPU) to reduce runtime. RESULTS We demonstrate the analyses using the tool are reproducible, up to 18 × faster using the GPU, and scale linearly with increasing methylation loci. CONCLUSIONS Torch-eCpG is a fast, reliable, and scalable tool to perform eQTM mapping. Source code for Torch-eCpG is available at https://github.com/kordk/torch-ecpg .
Collapse
Affiliation(s)
- Kord M Kober
- School of Nursing, University of California San Francisco, San Francisco, CA, USA.
- Helen Diller Family Comprehensive Cancer Center, University of California San Francisco, San Francisco, CA, USA.
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA.
| | - Liam Berger
- School of Nursing, University of California San Francisco, San Francisco, CA, USA
| | - Ritu Roy
- Helen Diller Family Comprehensive Cancer Center, University of California San Francisco, San Francisco, CA, USA
| | - Adam Olshen
- Helen Diller Family Comprehensive Cancer Center, University of California San Francisco, San Francisco, CA, USA
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
| |
Collapse
|
14
|
Kwong A, Zawistowski M, Fritsche LG, Zhan X, Bragg-Gresham J, Branham KE, Advani J, Othman M, Ratnapriya R, Teslovich TM, Stambolian D, Chew EY, Abecasis GR, Swaroop A. Whole genome sequencing of 4,787 individuals identifies gene-based rare variants in age-related macular degeneration. Hum Mol Genet 2024; 33:374-385. [PMID: 37934784 PMCID: PMC10840384 DOI: 10.1093/hmg/ddad189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2023] [Revised: 10/12/2023] [Accepted: 10/31/2023] [Indexed: 11/09/2023] Open
Abstract
Genome-wide association studies have contributed extensively to the discovery of disease-associated common variants. However, the genetic contribution to complex traits is still largely difficult to interpret. We report a genome-wide association study of 2394 cases and 2393 controls for age-related macular degeneration (AMD) via whole-genome sequencing, with 46.9 million genetic variants. Our study reveals significant single-variant association signals at four loci and independent gene-based signals in CFH, C2, C3, and NRTN. Using data from the Exome Aggregation Consortium (ExAC) for a gene-based test, we demonstrate an enrichment of predicted rare loss-of-function variants in CFH, CFI, and an as-yet unreported gene in AMD, ORMDL2. Our method of using a large variant list without individual-level genotypes as an external reference provides a flexible and convenient approach to leverage the publicly available variant datasets to augment the search for rare variant associations, which can explain additional disease risk in AMD.
Collapse
Affiliation(s)
- Alan Kwong
- Department of Biostatistics and Center for Statistical Genetics, University of Michigan, 1415 Washington Heights, Ann Arbor, MI 48109, United States
| | - Matthew Zawistowski
- Department of Biostatistics and Center for Statistical Genetics, University of Michigan, 1415 Washington Heights, Ann Arbor, MI 48109, United States
| | - Lars G Fritsche
- Department of Biostatistics and Center for Statistical Genetics, University of Michigan, 1415 Washington Heights, Ann Arbor, MI 48109, United States
| | - Xiaowei Zhan
- Southwestern Medical Center, University of Texas, 5323 Harry Hines Blvd, Dallas, TX 75390, United States
| | - Jennifer Bragg-Gresham
- Kidney Epidemiology and Cost Center, Department of Internal Medicine-Nephrology, University of Michigan, 1415 Washington Heights, Ann Arbor, MI 48109, United States
| | - Kari E Branham
- Department of Ophthalmology and Visual Sciences, University of Michigan Kellogg Eye Center, 1000 Wall St, Ann Arbor, MI 48105, United States
| | - Jayshree Advani
- Neurobiology-Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, 6 Center Drive, MSC 0610, Bethesda, MD 20892, United States
| | - Mohammad Othman
- Department of Ophthalmology and Visual Sciences, University of Michigan Kellogg Eye Center, 1000 Wall St, Ann Arbor, MI 48105, United States
| | - Rinki Ratnapriya
- Neurobiology-Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, 6 Center Drive, MSC 0610, Bethesda, MD 20892, United States
| | - Tanya M Teslovich
- Regeneron Pharmaceuticals Inc., 777 Old Saw Mill River Rd, Tarrytown, NY 10591, United States
| | - Dwight Stambolian
- Department of Ophthalmology, Perelman School of Medicine, University of Pennsylvania Medical School, 51 N. 39th Street, Philadelphia, PA 19104, United States
| | - Emily Y Chew
- Division of Epidemiology and Clinical Application, National Eye Institute, National Institutes of Health, 10 Center Drive Building 10-CRC, Bethesda, MD 20892, United States
| | - Gonçalo R Abecasis
- Department of Biostatistics and Center for Statistical Genetics, University of Michigan, 1415 Washington Heights, Ann Arbor, MI 48109, United States
- Regeneron Pharmaceuticals Inc., 777 Old Saw Mill River Rd, Tarrytown, NY 10591, United States
| | - Anand Swaroop
- Neurobiology-Neurodegeneration and Repair Laboratory, National Eye Institute, National Institutes of Health, 6 Center Drive, MSC 0610, Bethesda, MD 20892, United States
| |
Collapse
|
15
|
Varshney A, Manickam N, Orchard P, Tovar A, Zhang Z, Feng F, Erdos MR, Narisu N, Ventresca C, Nishino K, Rai V, Stringham HM, Jackson AU, Tamsen T, Gao C, Yang M, Koues OI, Welch JD, Burant CF, Williams LK, Jenkinson C, DeFronzo RA, Norton L, Saramies J, Lakka TA, Laakso M, Tuomilehto J, Mohlke KL, Kitzman JO, Koistinen HA, Liu J, Boehnke M, Collins FS, Scott LJ, Parker SCJ. Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.15.571696. [PMID: 38168419 PMCID: PMC10760134 DOI: 10.1101/2023.12.15.571696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.
Collapse
Affiliation(s)
- Arushi Varshney
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Nandini Manickam
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Peter Orchard
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Adelaide Tovar
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Zhenhao Zhang
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Fan Feng
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Michael R Erdos
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Narisu Narisu
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Christa Ventresca
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
- Dept. of Human Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Kirsten Nishino
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Vivek Rai
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Heather M Stringham
- Department of Biostatistics, Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Anne U Jackson
- Department of Biostatistics, Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Tricia Tamsen
- Biomedical Research Core Facilities Advanced Genomics Core, University of Michigan, Ann Arbor, MI, USA
| | - Chao Gao
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Mao Yang
- Department of Internal Medicine, Center for Individualized and Genomic Medicine Research, Henry Ford Hospital, Detroit, MI, USA
| | - Olivia I Koues
- Biomedical Research Core Facilities Advanced Genomics Core, University of Michigan, Ann Arbor, MI, USA
| | - Joshua D Welch
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Charles F Burant
- Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, USA
| | - L Keoki Williams
- Department of Internal Medicine, Center for Individualized and Genomic Medicine Research, Henry Ford Hospital, Detroit, MI, USA
| | - Chris Jenkinson
- South Texas Diabetes and Obesity Research Institute, School of Medicine, University of Texas, Rio Grande Valley, TX, USA
| | - Ralph A DeFronzo
- Department of Medicine/Diabetes Division, University of Texas Health, San Antonio, TX, USA
| | - Luke Norton
- Department of Medicine/Diabetes Division, University of Texas Health, San Antonio, TX, USA
| | - Jouko Saramies
- Savitaipale Health Center, South Karelia Central Hospital, Lappeenranta, Finland
| | - Timo A Lakka
- Institute of Biomedicine, University of Eastern Finland, Kuopio, Finland
| | - Markku Laakso
- Institute of Clinical Medicine, University of Eastern Finland, Kuopio, Finland
| | - Jaakko Tuomilehto
- Dept. of Public Health and Welfare, Finnish Institute for Health and Welfare, Helsinki, Finland
- Dept. of Public Health, University of Helsinki, Helsinki, Finland
- Diabetes Research Group, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Karen L Mohlke
- Dept. of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | - Jacob O Kitzman
- Dept. of Human Genetics, University of Michigan, Ann Arbor, MI, USA
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Heikki A Koistinen
- Dept. of Public Health and Welfare, Finnish Institute for Health and Welfare, Helsinki, Finland
- Department of Medicine, University of Helsinki and Helsinki University Hospital, Helsinki, Finland
- Minerva Foundation Institute for Medical Research, Helsinki, Finland
| | - Jie Liu
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Michael Boehnke
- Department of Biostatistics, Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Francis S Collins
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Laura J Scott
- Department of Biostatistics, Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Stephen C J Parker
- Dept. of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
- Dept. of Human Genetics, University of Michigan, Ann Arbor, MI, USA
- Department of Biostatistics, Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
16
|
Kim S, Qin Y, Park HJ, Yue M, Xu Z, Forno E, Chen W, Celedón JC. Methyl-TWAS: A powerful method for in silico transcriptome-wide association studies (TWAS) using long-range DNA methylation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.10.566586. [PMID: 38014125 PMCID: PMC10680683 DOI: 10.1101/2023.11.10.566586] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
In silico transcriptome-wide association studies (TWAS) are commonly used to test whether expression of specific genes is linked to a complex trait. However, genotype-based in silico TWAS such as PrediXcan, exhibit low prediction accuracy for a majority of genes because genotypic data lack tissue- and disease-specificity and are not affected by the environment. Because methylation is tissue-specific and, like gene expression, can be modified by environment or disease status, methylation should predict gene expression with more accuracy than SNPs. Therefore, we propose Methyl-TWAS, the first approach that utilizes long-range methylation markers to impute gene expression for in silico TWAS through penalized regression. Methyl-TWAS 1) predicts epigenetically regulated/associated expression (eGReX), which incorporates tissue-specific expression and both genetically- (GReX) and environmentally-regulated expression to identify differentially expressed genes (DEGs) that could not be identified by genotype-based methods; and 2) incorporates both cis- and trans- CpGs, including various regulatory regions to identify DEGs that would be missed using cis- methylation only. Methyl-TWAS outperforms PrediXcan and two other methods in imputing gene expression in the nasal epithelium, particularly for immunity-related genes and DEGs in atopic asthma. Methyl-TWAS identified 3,681 (85.2%) of the 4,316 DEGs identified in a previous TWAS of atopic asthma using measured expression, while PrediXcan could not identify any gene. Methyl-TWAS also outperforms PrediXcan for expression imputation as well as in silico TWAS in white blood cells. Methyl-TWAS is a valuable tool for in silico TWAS, leveraging a growing body of publicly available genome-wide DNA methylation data for a variety of human tissues.
Collapse
Affiliation(s)
- Soyeon Kim
- Division of Pulmonary Medicine, Department of Pediatrics, UPMC Children’s Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, PA, USA
| | - Yidi Qin
- Department of Human Genetics, School of Public Health, University of Pittsburgh, Pittsburgh, PA, USA
| | - Hyun Jung Park
- Department of Human Genetics, School of Public Health, University of Pittsburgh, Pittsburgh, PA, USA
| | - Molin Yue
- Department of Biostatistics, School of Public Health, University of Pittsburgh, Pittsburgh, PA, USA
| | - Zhongli Xu
- Division of Pulmonary Medicine, Department of Pediatrics, UPMC Children’s Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, PA, USA
- School of Medicine, Tsinghua University, Beijing, China
| | - Erick Forno
- Division of Pulmonary Medicine, Department of Pediatrics, UPMC Children’s Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, PA, USA
| | - Wei Chen
- Division of Pulmonary Medicine, Department of Pediatrics, UPMC Children’s Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, PA, USA
| | - Juan C. Celedón
- Division of Pulmonary Medicine, Department of Pediatrics, UPMC Children’s Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
17
|
Du Z, Iyyanki T, Lessard S, Chao M, Asbrand C, Nassar D, Klinger K, de Rinaldis E, Khader S, Chatelain C. Genome-wide association study analysis of disease severity in Acne reveals novel biological insights. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2023:2023.11.13.23298473. [PMID: 38014089 PMCID: PMC10680891 DOI: 10.1101/2023.11.13.23298473] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Acne vulgaris is a common skin disease that affects >85% of teenage young adults among which >8% develop severe lesions that leaves permanent scars. Genetic heritability studies of acne in twin cohorts have estimated that the heritability for acne is 80%. Previous genome-wide association studies (GWAS) have identified 50 genetic loci associated with increased risk of developing acne when compared to healthy individuals. However only a few studies have investigated genetic association with disease severity. GWAS of disease progression may provide a more effective approach to unveil potential disease modifying therapeutic targets. Here, we performed a multi-ethnic GWAS analysis to capture disease severity in acne patients by using individuals with normal acne as a control. Our cohort consists of a total of 2,956 participants, including 290 severe acne cases and 930 normal acne controls from FinnGen, and 522 cases and 1,214 controls from BioVU. We also performed mendelian randomization (MR), colocalization analyses and transcriptome-wide association study (TWAS) to identify putative causal genes. Lastly, we performed gene-set enrichment analysis using MAGMA to implicate biological pathways that drive disease severity in Acne. We identified two new loci associated with acne severity at the genome-wide significance level, six novel associated genes by MR, colocalization and TWAS analyses, including genes CDC7, SLC7A1, ADAM23, TTLL10, CDK20 and DNAJA4 , and 5 novel pathways by MAGMA analyses. Our study suggests that the etiologies of acne susceptibility and severity have limited overlap, with only 26% of known acne risk loci presenting nominal association with acne severity and none of the novel severity associated genes reported as associated with acne risk in previous GWAS.
Collapse
|
18
|
Zhu Z, Chen X, Zhang S, Yu R, Qi C, Cheng L, Zhang X. Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective. Hum Genet 2023; 142:1543-1560. [PMID: 37755483 DOI: 10.1007/s00439-023-02602-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 09/14/2023] [Indexed: 09/28/2023]
Abstract
Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.
Collapse
Affiliation(s)
- Zijun Zhu
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, Heilongjiang, China
| | - Xinyu Chen
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, Heilongjiang, China
| | - Sainan Zhang
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, Heilongjiang, China
| | - Rui Yu
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, Heilongjiang, China
| | - Changlu Qi
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, Heilongjiang, China
| | - Liang Cheng
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, Heilongjiang, China.
- NHC Key Laboratory of Molecular Probe and Targeted Diagnosis and Therapy, Harbin Medical University, Harbin, 150028, Heilongjiang, China.
| | - Xue Zhang
- NHC Key Laboratory of Molecular Probe and Targeted Diagnosis and Therapy, Harbin Medical University, Harbin, 150028, Heilongjiang, China
- McKusick-Zhang Center for Genetic Medicine, State Key Laboratory of Medical Molecular Biology, Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, 100005, China
| |
Collapse
|
19
|
Li T, Ferraro N, Strober BJ, Aguet F, Kasela S, Arvanitis M, Ni B, Wiel L, Hershberg E, Ardlie K, Arking DE, Beer RL, Brody J, Blackwell TW, Clish C, Gabriel S, Gerszten R, Guo X, Gupta N, Johnson WC, Lappalainen T, Lin HJ, Liu Y, Nickerson DA, Papanicolaou G, Pritchard JK, Qasba P, Shojaie A, Smith J, Sotoodehnia N, Taylor KD, Tracy RP, Van Den Berg D, Wheeler MT, Rich SS, Rotter JI, Battle A, Montgomery SB. The functional impact of rare variation across the regulatory cascade. CELL GENOMICS 2023; 3:100401. [PMID: 37868038 PMCID: PMC10589633 DOI: 10.1016/j.xgen.2023.100401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Revised: 03/08/2023] [Accepted: 08/10/2023] [Indexed: 10/24/2023]
Abstract
Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.
Collapse
Affiliation(s)
- Taibo Li
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Nicole Ferraro
- Biomedical Informatics Training Program, Stanford University, Stanford, CA, USA
| | - Benjamin J. Strober
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Harvard School of Public Health, Epidemiology Department, Boston, MA, USA
| | | | - Silva Kasela
- New York Genome Center, New York, NY, USA
- Department of Systems Biology, Columbia University, New York, NY, USA
| | - Marios Arvanitis
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Medicine, Division of Cardiology, Johns Hopkins School of Medicine, Baltimore, MD, USA
| | - Bohan Ni
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA
| | - Laurens Wiel
- Division of Cardiovascular Medicine, Department of Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | | | | | - Dan E. Arking
- McKusick-Nathans Institute, Department of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Rebecca L. Beer
- National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, USA
| | - Jennifer Brody
- Cardiovascular Health Research Unit, Departments of Medicine and Epidemiology, University of Washington, Seattle, WA, USA
| | - Thomas W. Blackwell
- Department of Biostatistics, School of Public Health, University of Michigan, Ann Arbor, MI, USA
| | - Clary Clish
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | | | - Robert Gerszten
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Cardiovascular Institute, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, USA
| | - Xiuqing Guo
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Namrata Gupta
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - W. Craig Johnson
- Collaborative Health Studies Coordinating Center, University of Washington, Seattle, WA, USA
| | - Tuuli Lappalainen
- New York Genome Center, New York, NY, USA
- Department of Systems Biology, Columbia University, New York, NY, USA
| | - Henry J. Lin
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Yongmei Liu
- Department of Medicine, Duke University School of Medicine, Durham, NC, USA
| | | | - George Papanicolaou
- National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, USA
| | | | - Pankaj Qasba
- National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, USA
| | - Ali Shojaie
- Department of Biostatistics, University of Washington School of Public Health, Seattle, WA, USA
| | - Josh Smith
- Department of Genome Sciences, University of Washington, Seattle, WA, USA
| | - Nona Sotoodehnia
- Cardiovascular Health Research Unit, Departments of Medicine and Epidemiology, University of Washington, Seattle, WA, USA
| | - Kent D. Taylor
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Russell P. Tracy
- Laboratory for Clinical Biochemistry Research, University of Vermont, Burlington, VT, USA
| | - David Van Den Berg
- Department of Preventive Medicine, University of Southern California, Los Angeles, CA, USA
| | - Matthew T. Wheeler
- Division of Cardiovascular Medicine, Department of Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Stephen S. Rich
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA, USA
| | - Jerome I. Rotter
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Alexis Battle
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA
- McKusick-Nathans Institute, Department of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Malone Center for Engineering of Healthcare, Johns Hopkins University, Baltimore, MD, USA
| | - Stephen B. Montgomery
- Department of Genetics, Stanford University, Stanford, CA, USA
- Department of Pathology, Stanford University, Stanford, CA, USA
| |
Collapse
|
20
|
Kim S, Xu Z, Forno E, Qin Y, Park HJ, Yue M, Yan Q, Manni ML, Acosta-Pérez E, Canino G, Chen W, Celedón JC. Cis- and trans-eQTM analysis reveals novel epigenetic and transcriptomic immune markers of atopic asthma in airway epithelium. J Allergy Clin Immunol 2023; 152:887-898. [PMID: 37271320 PMCID: PMC10592527 DOI: 10.1016/j.jaci.2023.05.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Revised: 04/03/2023] [Accepted: 05/02/2023] [Indexed: 06/06/2023]
Abstract
BACKGROUND Expression quantitative trait methylation (eQTM) analyses uncover associations between DNA methylation markers and gene expression. Most eQTM analyses of complex diseases have focused on cis-eQTM pairs (within 1 megabase). OBJECTIVES This study sought to identify cis- and trans-methylation markers associated with gene expression in airway epithelium from youth with and without atopic asthma. METHODS In this study, the investigators conducted both cis- and trans-eQTM analyses in nasal (airway) epithelial samples from 158 Puerto Rican youth with atopic asthma and 100 control subjects without atopy or asthma. The investigators then attempted to replicate their findings in nasal epithelial samples from 2 studies of children, while also examining whether their results in nasal epithelium overlap with those from an eQTM analysis in white blood cells from the Puerto Rican subjects. RESULTS This study identified 9,108 cis-eQTM pairs and 2,131,500 trans-eQTM pairs. Trans-associations were significantly enriched for transcription factor and microRNA target genes. Furthermore, significant cytosine-phosphate-guanine sites (CpGs) were differentially methylated in atopic asthma and significant genes were enriched for genes differentially expressed in atopic asthma. In this study, 50.7% to 62.6% of cis- and trans-eQTM pairs identified in Puerto Rican youth were replicated in 2 smaller cohorts at false discovery rate-adjusted P < .1. Replicated genes in the trans-eQTM analysis included biologically plausible asthma-susceptibility genes (eg, HDC, NLRP3, ITGAE, CDH26, and CST1) and are enriched in immune pathways. CONCLUSIONS Studying both cis- and trans-epigenetic regulation of airway epithelial gene expression can identify potential causal and regulatory pathways or networks for childhood asthma. Trans-eQTM CpGs may regulate gene expression in airway epithelium through effects on transcription factor and microRNA target genes.
Collapse
Affiliation(s)
- Soyeon Kim
- Division of Pediatric Pulmonary Medicine, University of Pittsburgh Medical Center Children's Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, Pa; Department of Pediatrics, School of Medicine, University of Pittsburgh, Pittsburgh, Pa
| | - Zhongli Xu
- Division of Pediatric Pulmonary Medicine, University of Pittsburgh Medical Center Children's Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, Pa; School of Medicine, Tsinghua University, Beijing, China
| | - Erick Forno
- Division of Pediatric Pulmonary Medicine, University of Pittsburgh Medical Center Children's Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, Pa; Department of Pediatrics, School of Medicine, University of Pittsburgh, Pittsburgh, Pa
| | - Yidi Qin
- Department of Human Genetics, School of Public Health, University of Pittsburgh, Pittsburgh, Pa
| | - Hyun Jung Park
- Department of Human Genetics, School of Public Health, University of Pittsburgh, Pittsburgh, Pa
| | - Molin Yue
- Department of Biostatistics, School of Public Health, University of Pittsburgh, Pittsburgh, Pa
| | - Qi Yan
- Department of Obstetrics and Gynecology, Columbia University, New York, NY
| | - Michelle L Manni
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh
| | - Edna Acosta-Pérez
- Behavioral Sciences Research Institute, University of Puerto Rico, San Juan, Puerto Rico
| | - Glorisa Canino
- Behavioral Sciences Research Institute, University of Puerto Rico, San Juan, Puerto Rico
| | - Wei Chen
- Division of Pediatric Pulmonary Medicine, University of Pittsburgh Medical Center Children's Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, Pa; School of Medicine, Tsinghua University, Beijing, China
| | - Juan C Celedón
- Division of Pediatric Pulmonary Medicine, University of Pittsburgh Medical Center Children's Hospital of Pittsburgh, University of Pittsburgh, Pittsburgh, Pa; Department of Pediatrics, School of Medicine, University of Pittsburgh, Pittsburgh, Pa.
| |
Collapse
|
21
|
Maimaiti A, Turhon M, Abulaiti A, Dilixiati Y, Zhang F, Axieer A, Kadeer K, Zhang Y, Maimaitili A, Yang X. DNA methylation regulator-mediated modification patterns and risk of intracranial aneurysm: a multi-omics and epigenome-wide association study integrating machine learning, Mendelian randomization, eQTL and mQTL data. J Transl Med 2023; 21:660. [PMID: 37742034 PMCID: PMC10518114 DOI: 10.1186/s12967-023-04512-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 09/08/2023] [Indexed: 09/25/2023] Open
Abstract
BACKGROUND Intracranial aneurysms (IAs) pose a significant and intricate challenge. Elucidating the interplay between DNA methylation and IA pathogenesis is paramount to identify potential biomarkers and therapeutic interventions. METHODS We employed a comprehensive bioinformatics investigation of DNA methylation in IA, utilizing a transcriptomics-based methodology that encompassed 100 machine learning algorithms, genome-wide association studies (GWAS), Mendelian randomization (MR), and summary-data-based Mendelian randomization (SMR). Our sophisticated analytical strategy allowed for a systematic assessment of differentially methylated genes and their implications on the onset, progression, and rupture of IA. RESULTS We identified DNA methylation-related genes (MRGs) and associated molecular pathways, and the MR and SMR analyses provided evidence for potential causal links between the observed DNA methylation events and IA predisposition. CONCLUSION These insights not only augment our understanding of the molecular underpinnings of IA but also underscore potential novel biomarkers and therapeutic avenues. Although our study faces inherent limitations and hurdles, it represents a groundbreaking initiative in deciphering the intricate relationship between genetic, epigenetic, and environmental factors implicated in IA pathogenesis.
Collapse
Affiliation(s)
- Aierpati Maimaiti
- Department of Neurosurgery, Xinjiang Medical University Affiliated First Hospital, Urumqi, Xinjiang, 830017, People's Republic of China
| | - Mirzat Turhon
- Department of Interventional Neuroradiology, Beijing Neurosurgical Institute, Capital Medical University, 100070, Beijing, People's Republic of China
- Department of Interventional Neuroradiology, Beijing TianTan Hospital, Capital Medical University, Beijing, People's Republic of China
| | - Aimitaji Abulaiti
- Xinjiang Medical University, Urumqi, Xinjiang, People's Republic of China
| | | | - Fujunhui Zhang
- Department of Interventional Neuroradiology, Beijing Neurosurgical Institute, Capital Medical University, 100070, Beijing, People's Republic of China
- Department of Interventional Neuroradiology, Beijing TianTan Hospital, Capital Medical University, Beijing, People's Republic of China
| | - Aximujiang Axieer
- Department of Neurosurgery, Xinjiang Medical University Affiliated First Hospital, Urumqi, Xinjiang, 830017, People's Republic of China
| | - Kaheerman Kadeer
- Department of Neurosurgery, Xinjiang Medical University Affiliated First Hospital, Urumqi, Xinjiang, 830017, People's Republic of China
| | - Yisen Zhang
- Department of Interventional Neuroradiology, Beijing Neurosurgical Institute, Capital Medical University, 100070, Beijing, People's Republic of China.
- Department of Interventional Neuroradiology, Beijing TianTan Hospital, Capital Medical University, Beijing, People's Republic of China.
| | - Aisha Maimaitili
- Department of Neurosurgery, Xinjiang Medical University Affiliated First Hospital, Urumqi, Xinjiang, 830017, People's Republic of China.
| | - Xinjian Yang
- Department of Interventional Neuroradiology, Beijing Neurosurgical Institute, Capital Medical University, 100070, Beijing, People's Republic of China.
- Department of Interventional Neuroradiology, Beijing TianTan Hospital, Capital Medical University, Beijing, People's Republic of China.
| |
Collapse
|
22
|
Kerimov N, Tambets R, Hayhurst JD, Rahu I, Kolberg P, Raudvere U, Kuzmin I, Chowdhary A, Vija A, Teras HJ, Kanai M, Ulirsch J, Ryten M, Hardy J, Guelfi S, Trabzuni D, Kim-Hellmuth S, Rayner W, Finucane H, Peterson H, Mosaku A, Parkinson H, Alasoo K. eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs. PLoS Genet 2023; 19:e1010932. [PMID: 37721944 PMCID: PMC10538656 DOI: 10.1371/journal.pgen.1010932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 09/28/2023] [Accepted: 08/22/2023] [Indexed: 09/20/2023] Open
Abstract
The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.
Collapse
Affiliation(s)
- Nurlan Kerimov
- Institute of Computer Science, University of Tartu, Tartu, Estonia
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Ralf Tambets
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - James D. Hayhurst
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Ida Rahu
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Peep Kolberg
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Uku Raudvere
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Ivan Kuzmin
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Anshika Chowdhary
- Institute of Translational Genomics, Helmholtz Munich, Neuherberg, Germany
| | - Andreas Vija
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Hans J. Teras
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Masahiro Kanai
- Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, Massachusetts, United States of America
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
| | - Jacob Ulirsch
- Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, Massachusetts, United States of America
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
| | - Mina Ryten
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London, United Kingdom
| | - John Hardy
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London, United Kingdom
| | - Sebastian Guelfi
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London, United Kingdom
| | - Daniah Trabzuni
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London, United Kingdom
| | - Sarah Kim-Hellmuth
- Institute of Translational Genomics, Helmholtz Munich, Neuherberg, Germany
- Department of Pediatrics, Dr. von Hauner Children’s Hospital, University Hospital LMU Munich, Munich, Germany
| | - William Rayner
- Institute of Translational Genomics, Helmholtz Munich, Neuherberg, Germany
| | - Hilary Finucane
- Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, Massachusetts, United States of America
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
| | - Hedi Peterson
- Institute of Computer Science, University of Tartu, Tartu, Estonia
| | - Abayomi Mosaku
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Helen Parkinson
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Kaur Alasoo
- Institute of Computer Science, University of Tartu, Tartu, Estonia
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| |
Collapse
|
23
|
Mortlock S, Houshdaran S, Kosti I, Rahmioglu N, Nezhat C, Vitonis AF, Andrews SV, Grosjean P, Paranjpe M, Horne AW, Jacoby A, Lager J, Opoku-Anane J, Vo KC, Manvelyan E, Sen S, Ghukasyan Z, Collins F, Santamaria X, Saunders P, Kober K, McRae AF, Terry KL, Vallvé-Juanico J, Becker C, Rogers PAW, Irwin JC, Zondervan K, Montgomery GW, Missmer S, Sirota M, Giudice L. Global endometrial DNA methylation analysis reveals insights into mQTL regulation and associated endometriosis disease risk and endometrial function. Commun Biol 2023; 6:780. [PMID: 37587191 PMCID: PMC10432557 DOI: 10.1038/s42003-023-05070-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Accepted: 06/23/2023] [Indexed: 08/18/2023] Open
Abstract
Endometriosis is a leading cause of pain and infertility affecting millions of women globally. Herein, we characterize variation in DNA methylation (DNAm) and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genotype data and methylation in endometrial samples from 984 deeply-phenotyped participants. We estimate that 15.4% of the variation in endometriosis is captured by DNAm and identify significant differences in DNAm profiles associated with stage III/IV endometriosis, endometriosis sub-phenotypes and menstrual cycle phase, including opening of the window for embryo implantation. Menstrual cycle phase was a major source of DNAm variation suggesting cellular and hormonally-driven changes across the cycle can regulate genes and pathways responsible for endometrial physiology and function. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independent cis-mQTLs including 51 associated with risk of endometriosis, highlighting candidate genes contributing to disease risk. Our work provides functional evidence for epigenetic targets contributing to endometriosis risk and pathogenesis. Data generated serve as a valuable resource for understanding tissue-specific effects of methylation on endometrial biology in health and disease.
Collapse
Affiliation(s)
- Sally Mortlock
- The Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD, 4072, Australia.
| | - Sahar Houshdaran
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Idit Kosti
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
| | - Nilufer Rahmioglu
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
- Oxford Endometriosis CaRe Centre, Nuffield Department of Women's and Reproductive Health, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Camran Nezhat
- Stanford University Medical Center, Palo Alto, CA, USA
- University of California San Francisco, San Francisco, CA, USA
- Camran Nezhat Institute, Center for Special Minimally Invasive and Robotic Surgery, Woodside, CA, USA
| | - Allison F Vitonis
- Obstetrics and Gynecology Epidemiology Center, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Shan V Andrews
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
| | - Parker Grosjean
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
| | - Manish Paranjpe
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
| | - Andrew W Horne
- MRC Centre for Reproductive Health, University of Edinburgh, QMRI, Edinburgh, UK
| | - Alison Jacoby
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Jeannette Lager
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Jessica Opoku-Anane
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Kim Chi Vo
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Evelina Manvelyan
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Sushmita Sen
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Zhanna Ghukasyan
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Frances Collins
- MRC Centre for Reproductive Health, University of Edinburgh, QMRI, Edinburgh, UK
| | - Xavier Santamaria
- Carlos Simon Foundation, Health Research Institute, Valencia, Spain
- Group of Biomedical Research in Gynecology, Vall d'Hebron Research Institute, Barcelona, Spain
| | - Philippa Saunders
- Centre for Inflammation Research, Institute for Regeneration and Repair University of Edinburgh, Edinburgh, UK
| | - Kord Kober
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
- Department of Physiological Nursing, University of California San Francisco, San Francisco, CA, USA
| | - Allan F McRae
- The Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Kathryn L Terry
- Obstetrics and Gynecology Epidemiology Center, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Boston Center for Endometriosis, Boston Children's Hospital and Brigham and Women's Hospital, Boston, MA, USA
| | - Júlia Vallvé-Juanico
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
- Group of Biomedical Research in Gynecology, Vall d'Hebron Research Institute, Barcelona, Spain
| | - Christian Becker
- Oxford Endometriosis CaRe Centre, Nuffield Department of Women's and Reproductive Health, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Peter A W Rogers
- University of Melbourne Department of Obstetrics and Gynaecology, Royal Women's Hospital, Melbourne, Australia
| | - Juan C Irwin
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Krina Zondervan
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
- Oxford Endometriosis CaRe Centre, Nuffield Department of Women's and Reproductive Health, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Grant W Montgomery
- The Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Stacey Missmer
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Boston Center for Endometriosis, Boston Children's Hospital and Brigham and Women's Hospital, Boston, MA, USA
- Division of Adolescent and Young Adult Medicine, Department of Medicine, Boston Children's Hospital and Harvard Medical School, Boston, MA, USA
- Department of Obstetrics, Gynecology, and Reproductive Biology, College of Human Medicine, Michigan State University, Grand Rapids, MI, USA
| | - Marina Sirota
- Bakar Computational Health Sciences Institute, University of California San Francisco, San Francisco, CA, USA
- Department of Pediatrics, Division of Neonatology, University of California San Francisco, San Francisco, CA, USA
| | - Linda Giudice
- Center for Reproductive Sciences, Department of Obstetrics, Gynecology & Reproductive Sciences, University of California San Francisco, San Francisco, CA, USA.
| |
Collapse
|
24
|
Grenko CM, Bonnycastle LL, Taylor HJ, Yan T, Swift AJ, Robertson CC, Narisu N, Erdos MR, Collins FS, Taylor DL. Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 hours. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.06.543931. [PMID: 37333221 PMCID: PMC10274787 DOI: 10.1101/2023.06.06.543931] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycemia, beta cell glucotoxicity, and ultimately type 2 diabetes (T2D). In this study, we sought to explore the effects of hyperglycemia on human pancreatic islet (HPI) gene expression by exposing HPIs from two donors to low (2.8mM) and high (15.0mM) glucose concentrations over 24 hours, assaying the transcriptome at seven time points using single-cell RNA sequencing (scRNA-seq). We modeled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Across all cell types, we identified 1,528 genes associated with time, 1,185 genes associated with glucose exposure, and 845 genes associated with interaction effects between time and glucose. We clustered differentially expressed genes across cell types and found 347 modules of genes with similar expression patterns across time and glucose conditions, including two beta cell modules enriched in genes associated with T2D. Finally, by integrating genomic features from this study and genetic summary statistics for T2D and related traits, we nominate 363 candidate effector genes that may underlie genetic associations for T2D and related traits.
Collapse
Affiliation(s)
- Caleb M. Grenko
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - Lori L. Bonnycastle
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - Henry J. Taylor
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
- British Heart Foundation Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, Cambridge CB1 8RN, UK
| | - Tingfen Yan
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - Amy J. Swift
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - Catherine C. Robertson
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Narisu Narisu
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - Michael R. Erdos
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - Francis S. Collins
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| | - D. Leland Taylor
- Center for Precision Health Research, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD 20892, USA
| |
Collapse
|
25
|
Advani J, Corso-Diaz X, Kwicklis M, van Asten F, Ratnapriya R, Mehta P, Hamel A, Mahrotra S, Segrè A, Kiel C, Strunz T, Weber B, Chew E, Hernandez D, Montezuma S, Ferrington D, Swaroop A. QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration. RESEARCH SQUARE 2023:rs.3.rs-3011096. [PMID: 37398472 PMCID: PMC10312909 DOI: 10.21203/rs.3.rs-3011096/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/04/2023]
Abstract
DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.
Collapse
Affiliation(s)
| | | | | | | | | | - Puja Mehta
- Department of Ophthalmology, Massachusetts Eye and Ear, Harvard Medical School, Boston, MA, USA
| | - Andrew Hamel
- Department of Ophthalmology, Massachusetts Eye and Ear
| | | | | | | | | | | | - Emily Chew
- National Eye Institute/National Institutes of Health
| | | | | | | | - Anand Swaroop
- National Eye Institute, National Institutes of Health
| |
Collapse
|
26
|
Williamson A, Norris DM, Yin X, Broadaway KA, Moxley AH, Vadlamudi S, Wilson EP, Jackson AU, Ahuja V, Andersen MK, Arzumanyan Z, Bonnycastle LL, Bornstein SR, Bretschneider MP, Buchanan TA, Chang YC, Chuang LM, Chung RH, Clausen TD, Damm P, Delgado GE, de Mello VD, Dupuis J, Dwivedi OP, Erdos MR, Fernandes Silva L, Frayling TM, Gieger C, Goodarzi MO, Guo X, Gustafsson S, Hakaste L, Hammar U, Hatem G, Herrmann S, Højlund K, Horn K, Hsueh WA, Hung YJ, Hwu CM, Jonsson A, Kårhus LL, Kleber ME, Kovacs P, Lakka TA, Lauzon M, Lee IT, Lindgren CM, Lindström J, Linneberg A, Liu CT, Luan J, Aly DM, Mathiesen E, Moissl AP, Morris AP, Narisu N, Perakakis N, Peters A, Prasad RB, Rodionov RN, Roll K, Rundsten CF, Sarnowski C, Savonen K, Scholz M, Sharma S, Stinson SE, Suleman S, Tan J, Taylor KD, Uusitupa M, Vistisen D, Witte DR, Walther R, Wu P, Xiang AH, Zethelius B, Ahlqvist E, Bergman RN, Chen YDI, Collins FS, Fall T, Florez JC, Fritsche A, Grallert H, Groop L, Hansen T, Koistinen HA, Komulainen P, Laakso M, Lind L, Loeffler M, März W, Meigs JB, Raffel LJ, Rauramaa R, Rotter JI, Schwarz PEH, Stumvoll M, Sundström J, Tönjes A, Tuomi T, Tuomilehto J, Wagner R, Barroso I, Walker M, Grarup N, Boehnke M, Wareham NJ, Mohlke KL, Wheeler E, O'Rahilly S, Fazakerley DJ, Langenberg C. Genome-wide association study and functional characterization identifies candidate genes for insulin-stimulated glucose uptake. Nat Genet 2023; 55:973-983. [PMID: 37291194 PMCID: PMC7614755 DOI: 10.1038/s41588-023-01408-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 04/26/2023] [Indexed: 06/10/2023]
Abstract
Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P < 5 × 10-8) not previously associated with postchallenge insulin resistance, eight of which were shown to share their genetic architecture with type 2 diabetes in colocalization analyses. We investigated candidate genes at a subset of associated loci in cultured cells and identified nine candidate genes newly implicated in the expression or trafficking of GLUT4, the key glucose transporter in postprandial glucose uptake in muscle and fat. By focusing on postprandial insulin resistance, we highlighted the mechanisms of action at type 2 diabetes loci that are not adequately captured by studies of fasting glycemic traits.
Collapse
Affiliation(s)
- Alice Williamson
- MRC Epidemiology Unit Institute of Metabolic Science, University of Cambridge School of Clinical Medicine, Cambridge, UK
- Metabolic Research Laboratories Wellcome Trust-MRC Institute of Metabolic Science, Department of Clinical Biochemistry, University of Cambridge, Cambridge, UK
| | - Dougall M Norris
- Metabolic Research Laboratories Wellcome Trust-MRC Institute of Metabolic Science, Department of Clinical Biochemistry, University of Cambridge, Cambridge, UK
| | - Xianyong Yin
- Department of Biostatistics, University of Michigan, Ann Arbor, MI, USA
- Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
- Department of Epidemiology, School of Public Health, Nanjing Medical University, Nanjing, China
| | - K Alaine Broadaway
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | - Anne H Moxley
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | | | - Emma P Wilson
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | - Anne U Jackson
- Department of Biostatistics, University of Michigan, Ann Arbor, MI, USA
- Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Vasudha Ahuja
- Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland
| | - Mette K Andersen
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Zorayr Arzumanyan
- Department of Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Lori L Bonnycastle
- Center for Precision Health Research National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Stefan R Bornstein
- Department of Internal Medicine III, Metabolic and Vascular Medicine, Medical Faculty Carl Gustav Carus, Dresden, Germany
- Helmholtz Zentrum München Paul Langerhans Institute Dresden (PLID), University Hospital and Faculty of Medicine TU Dresden, Dresden, Germany
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
| | - Maxi P Bretschneider
- Department of Internal Medicine III, Metabolic and Vascular Medicine, Medical Faculty Carl Gustav Carus, Dresden, Germany
- Helmholtz Zentrum München Paul Langerhans Institute Dresden (PLID), University Hospital and Faculty of Medicine TU Dresden, Dresden, Germany
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
| | - Thomas A Buchanan
- Department of Medicine, Division of Endocrinology and Diabetes, Keck School of Medicine USC, Los Angeles, CA, USA
| | - Yi-Cheng Chang
- Graduate Institute of Medical Genomics and Proteomics, National Taiwan University, Taipei City, Taiwan
- Internal Medicine, National Taiwan University Hospital, Taipei City, Taiwan
- Institute of Biomedical Sciences, Academia Sinica, Taipei City, Taiwan
| | - Lee-Ming Chuang
- Department of Internal Medicine, Division of Endocrinology and Metabolism, National Taiwan University Hospital, Taipei City, Taiwan
| | - Ren-Hua Chung
- Institute of Population Health Sciences, National Health Research Institutes, Toufen, Taiwan
| | - Tine D Clausen
- Department of Gynecology and Obstetrics, Nordsjaellands Hospital, Hillerød, Denmark
- Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Peter Damm
- Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
- Center for Pregnant Women with Diabetes, Rigshospitalet, Copenhagen, Denmark
- Department of Obstetrics, Rigshospitalet, Copenhagen, Denmark
| | - Graciela E Delgado
- Vth Department of Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Vanessa D de Mello
- Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland
| | - Josée Dupuis
- Department of Biostatistics, Boston University School of Public Health, Boston, MA, USA
- Department of Epidemiology, Biostatistics and Occupational Health, McGill University, Montréal, Quebec, Canada
| | - Om P Dwivedi
- Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland
| | - Michael R Erdos
- Center for Precision Health Research National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | | | | | - Christian Gieger
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
- Institute of Epidemiology II, Helmholtz Zentrum München-German Research Center for Environmental Health, Neuherberg, Germany
| | - Mark O Goodarzi
- Department of Medicine, Division of Endocrinology, Diabetes and Metabolism Cedars-Sinai Medical Center, Los Angeles, CA, USA
| | - Xiuqing Guo
- Department of Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Stefan Gustafsson
- Department of Medical Sciences, Clinical Epidemiology, Uppsala University, Uppsala, Sweden
| | - Liisa Hakaste
- Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland
| | - Ulf Hammar
- Department of Medical Sciences, Molecular Epidemiology, Uppsala University, Uppsala, Sweden
| | - Gad Hatem
- Clinical Sciences Malmö, Genomics, Diabetes and Endocrinology, Lund University, Malmö, Sweden
| | - Sandra Herrmann
- Helmholtz Zentrum München Paul Langerhans Institute Dresden (PLID), University Hospital and Faculty of Medicine TU Dresden, Dresden, Germany
- Department of Internal Medicine III, Prevention and Care of Diabetes, Medical Faculty Carl Gustav Carus, Dresden, Germany
| | - Kurt Højlund
- Steno Diabetes Center Odense, Odense University Hospital, Odense, Denmark
| | - Katrin Horn
- Medical Faculty Institute for Medical Informatics, Statistics and Epidemiology, Leipzig, Germany
- LIFE Leipzig Research Center for Civilization Diseases, Medical Faculty, Leipzig, Germany
| | - Willa A Hsueh
- Internal Medicine, Endocrinology, Diabetes and Metabolism, The Ohio State University Wexner Medical Center, Columbus, OH, USA
| | - Yi-Jen Hung
- Institute of Preventive Medicine, National Defense Medical Center, New Taipei City, Taiwan
| | - Chii-Min Hwu
- Department of Medicine Section of Endocrinology and Metabolism, Taipei Veterans General Hospital, Taipei City, Taiwan
| | - Anna Jonsson
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Line L Kårhus
- Center for Clinical Research and Prevention, Copenhagen University Hospital - Bispebjerg and Frederiksberg, Copenhagen, Denmark
| | - Marcus E Kleber
- Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland
- SYNLAB MVZ Humangenetik Mannheim, Mannheim, Germany
| | - Peter Kovacs
- Medical Department III-Endocrinology, Nephrology, Rheumatology, University of Leipzig Medical Center, Leipzig, Germany
| | - Timo A Lakka
- Institute of Biomedicine, School of Medicine, University of Eastern Finland, Kuopio, Finland
- Department of Clinical Physiology and Nuclear Medicine, Kuopio University Hospital, Kuopio, Finland
- Foundation for Research in Health Exercise and Nutrition, Kuopio Research Institute of Exercise Medicine, Kuopio, Finland
| | - Marie Lauzon
- Department of Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - I-Te Lee
- Department of Internal Medicine Division of Endocrinology and Metabolism, Taichung Veterans General Hospital, Taichung City, Taiwan
- School of Medicine, National Yang Ming Chiao Tung University, Taipei City, Taiwan
- School of Medicine, Chung Shan Medical University, Taichung City, Taiwan
| | - Cecilia M Lindgren
- Big Data Institute Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, UK
- Nuffield Department of Population Health, University of Oxford, Oxford, UK
- Wellcome Trust Centre Human Genetics, University of Oxford, Oxford, UK
- Broad Institute, Cambridge, MA, USA
| | | | - Allan Linneberg
- Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
- Center for Clinical Research and Prevention, Copenhagen University Hospital - Bispebjerg and Frederiksberg, Copenhagen, Denmark
| | - Ching-Ti Liu
- Department of Biostatistics, Boston University School of Public Health, Boston, MA, USA
| | - Jian'an Luan
- MRC Epidemiology Unit Institute of Metabolic Science, University of Cambridge School of Clinical Medicine, Cambridge, UK
| | - Dina Mansour Aly
- Clinical Sciences Malmö, Genomics, Diabetes and Endocrinology, Lund University, Malmö, Sweden
| | - Elisabeth Mathiesen
- Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
- Center for Pregnant Women with Diabetes, Rigshospitalet, Copenhagen, Denmark
- Department of Endocrinology Rigshospitalet, Copenhagen, Denmark
| | - Angela P Moissl
- Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland
- Institute of Nutritional Sciences, Friedrich-Schiller-University, Jena, Germany
- Competence Cluster for Nutrition and Cardiovascular Health (nutriCARD) Halle-Jena, Jena, Germany
| | - Andrew P Morris
- Centre for Genetics and Genomics Versus Arthritis Centre for Musculoskeletal Research, The University of Manchester, Manchester, UK
| | - Narisu Narisu
- Center for Precision Health Research National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Nikolaos Perakakis
- Department of Internal Medicine III, Metabolic and Vascular Medicine, Medical Faculty Carl Gustav Carus, Dresden, Germany
- Helmholtz Zentrum München Paul Langerhans Institute Dresden (PLID), University Hospital and Faculty of Medicine TU Dresden, Dresden, Germany
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
| | - Annette Peters
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
- Institute of Epidemiology II, Helmholtz Zentrum München-German Research Center for Environmental Health, Neuherberg, Germany
| | - Rashmi B Prasad
- Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland
- Clinical Sciences Malmö, Genomics, Diabetes and Endocrinology, Lund University, Malmö, Sweden
| | - Roman N Rodionov
- Department of Internal Medicine III, University Center for Vascular Medicine, Medical Faculty Carl Gustav Carus, Dresden, Germany
- College of Medicine and Public Health, Flinders University and Flinders Medical Centre, Adelaide, Australia
| | - Kathryn Roll
- Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Carsten F Rundsten
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Chloé Sarnowski
- Department of Epidemiology, Human Genetics and Environmental Sciences, The University of Texas Health Science Center, Houston, TX, USA
| | - Kai Savonen
- Foundation for Research in Health Exercise and Nutrition, Kuopio Research Institute of Exercise Medicine, Kuopio, Finland
| | - Markus Scholz
- Medical Faculty Institute for Medical Informatics, Statistics and Epidemiology, Leipzig, Germany
- LIFE Leipzig Research Center for Civilization Diseases, Medical Faculty, Leipzig, Germany
| | - Sapna Sharma
- Research Unit of Molecular Epidemiology, Helmholtz Zentrum München-German Research Center for Environmental Health, Neuherberg, Germany
- Food Chemistry and Molecular and Sensory Science, Technical University of Munich, Freising-Weihenstephan, München, Germany
| | - Sara E Stinson
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Sufyan Suleman
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Jingyi Tan
- Department of Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Kent D Taylor
- Department of Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Matti Uusitupa
- Department of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland
| | - Dorte Vistisen
- Clinical Research, Steno Diabetes Center Copenhagen, Herlev, Denmark
- Department of Public Health, University of Copenhagen, Copenhagen, Denmark
| | - Daniel R Witte
- Steno Diabetes Center Aarhus, Aarhus, Denmark
- Department of Public Health, Aarhus University, Aarhus, Denmark
| | - Romy Walther
- Helmholtz Zentrum München Paul Langerhans Institute Dresden (PLID), University Hospital and Faculty of Medicine TU Dresden, Dresden, Germany
- Department of Internal Medicine III, Pathobiochemistry, Medical Faculty Carl Gustav Carus, Dresden, Germany
| | - Peitao Wu
- Department of Biostatistics, Boston University School of Public Health, Boston, MA, USA
| | - Anny H Xiang
- Research and Evaluation, Division of Biostatistics, Kaiser Permanente Southern California, Pasadena, CA, USA
| | - Björn Zethelius
- Department of Public Health and Caring Sciences, Geriatrics, Uppsala University, Uppsala, Sweden
| | - Emma Ahlqvist
- Clinical Sciences Malmö, Genomics, Diabetes and Endocrinology, Lund University, Malmö, Sweden
| | - Richard N Bergman
- Diabetes and Obesity Research Institute, Cedars-Sinai Medical Center, Los Angeles, CA, USA
| | - Yii-Der Ida Chen
- Department of Pediatrics, Genomic Outcomes, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Francis S Collins
- Center for Precision Health Research National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Tove Fall
- Department of Medical Sciences, Molecular Epidemiology, Uppsala University, Uppsala, Sweden
| | - Jose C Florez
- Diabetes Unit and Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA
- Programs in Metabolism and Medical and Population Genetics, The Broad Institute, Cambridge, MA, USA
- Department of Medicine, Harvard Medical School, Boston, MA, USA
| | - Andreas Fritsche
- Department of Internal Medicine IV, University Hospital Tübingen, Tübingen, Germany
| | - Harald Grallert
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
- Institute of Epidemiology II, Helmholtz Zentrum München-German Research Center for Environmental Health, Neuherberg, Germany
- Research Unit of Molecular Epidemiology, Helmholtz Zentrum München-German Research Center for Environmental Health, Neuherberg, Germany
| | - Leif Groop
- Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland
- Clinical Sciences Malmö, Genomics, Diabetes and Endocrinology, Lund University, Lund, Sweden
| | - Torben Hansen
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Heikki A Koistinen
- Department of Public Health and Welfare, Finnish Institute for Health and Welfare, Helsinki, Finland
- Department of Medicine, University of Helsinki and Helsinki University Hospital, Helsinki, Finland
- Minerva Foundation Institute for Medical Research, Helsinki, Finland
| | - Pirjo Komulainen
- Foundation for Research in Health Exercise and Nutrition, Kuopio Research Institute of Exercise Medicine, Kuopio, Finland
| | - Markku Laakso
- Institute of Clinical Medicine, University of Eastern Finland, Kuopio, Finland
| | - Lars Lind
- Department of Medical Sciences, Clinical Epidemiology, Uppsala University, Uppsala, Sweden
| | - Markus Loeffler
- Medical Faculty Institute for Medical Informatics, Statistics and Epidemiology, Leipzig, Germany
- LIFE Leipzig Research Center for Civilization Diseases, Medical Faculty, Leipzig, Germany
| | - Winfried März
- Vth Department of Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
- Synlab Academy, SYNLAB Holding Deutschland GmbH, Mannheim, Germany
| | - James B Meigs
- Department of Internal Medicine IV, University Hospital Tübingen, Tübingen, Germany
- Clinical Sciences Malmö, Genomics, Diabetes and Endocrinology, Lund University, Lund, Sweden
- Department of Medicine Division of General Internal Medicine, Massachusetts General Hospital, Boston, MA, USA
| | - Leslie J Raffel
- Department of Pediatrics, Genetic and Genomic Medicine, University of California, Irvine, CA, USA
| | - Rainer Rauramaa
- Foundation for Research in Health Exercise and Nutrition, Kuopio Research Institute of Exercise Medicine, Kuopio, Finland
| | - Jerome I Rotter
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Peter E H Schwarz
- Helmholtz Zentrum München Paul Langerhans Institute Dresden (PLID), University Hospital and Faculty of Medicine TU Dresden, Dresden, Germany
- German Center for Diabetes Research (DZD e.V.), Neuherberg, Germany
- Department of Internal Medicine III, Prevention and Care of Diabetes, Medical Faculty Carl Gustav Carus, Dresden, Germany
| | - Michael Stumvoll
- Medical Department III-Endocrinology, Nephrology, Rheumatology, University of Leipzig Medical Center, Leipzig, Germany
| | - Johan Sundström
- Department of Medical Sciences, Clinical Epidemiology, Uppsala University, Uppsala, Sweden
| | - Anke Tönjes
- Medical Department III-Endocrinology, Nephrology, Rheumatology, University of Leipzig Medical Center, Leipzig, Germany
| | - Tiinamaija Tuomi
- Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland
- Folkhälsan Research Center, Helsinki, Finland
| | - Jaakko Tuomilehto
- Department of Public Health, University of Helsinki, Helsinki, Finland
- Population Health Unit, Finnish Institute for Health and Welfare, Helsinki, Finland
- Diabetes Research Group, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Robert Wagner
- Department of Internal Medicine IV, University Hospital Tübingen, Tübingen, Germany
| | - Inês Barroso
- Exeter Centre of Excellence for Diabetes Research (EXCEED), Genetics of Complex Traits, University of Exeter Medical School, University of Exeter, Exeter, UK
| | - Mark Walker
- Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne, UK
| | - Niels Grarup
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Michael Boehnke
- Department of Biostatistics, University of Michigan, Ann Arbor, MI, USA
- Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA
| | - Nicholas J Wareham
- MRC Epidemiology Unit Institute of Metabolic Science, University of Cambridge School of Clinical Medicine, Cambridge, UK
| | - Karen L Mohlke
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA.
| | - Eleanor Wheeler
- MRC Epidemiology Unit Institute of Metabolic Science, University of Cambridge School of Clinical Medicine, Cambridge, UK.
| | - Stephen O'Rahilly
- Metabolic Research Laboratories Wellcome Trust-MRC Institute of Metabolic Science, Department of Clinical Biochemistry, University of Cambridge, Cambridge, UK.
| | - Daniel J Fazakerley
- Metabolic Research Laboratories Wellcome Trust-MRC Institute of Metabolic Science, Department of Clinical Biochemistry, University of Cambridge, Cambridge, UK.
| | - Claudia Langenberg
- MRC Epidemiology Unit Institute of Metabolic Science, University of Cambridge School of Clinical Medicine, Cambridge, UK.
- Computational Medicine, Berlin Institute of Health at Charité-Universitätsmedizin, Berlin, Germany.
- Precision Healthcare University Research Institute, Queen Mary University of London, London, UK.
| |
Collapse
|
27
|
Kerimov N, Tambets R, Hayhurst JD, Rahu I, Kolberg P, Raudvere U, Kuzmin I, Chowdhary A, Vija A, Teras HJ, Kanai M, Ulirsch J, Ryten M, Hardy J, Guelfi S, Trabzuni D, Kim-Hellmuth S, Rayner W, Finucane H, Peterson H, Mosaku A, Parkinson H, Alasoo K. Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.06.535816. [PMID: 37066341 PMCID: PMC10104061 DOI: 10.1101/2023.04.06.535816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/18/2023]
Abstract
Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.
Collapse
Affiliation(s)
- Nurlan Kerimov
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Ralf Tambets
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - James D Hayhurst
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Ida Rahu
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Peep Kolberg
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Uku Raudvere
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Ivan Kuzmin
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Anshika Chowdhary
- Institute of Translational Genomics, Helmholtz Munich, Neuherberg, Germany
| | - Andreas Vija
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Hans J Teras
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Masahiro Kanai
- Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, MA, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Jacob Ulirsch
- Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, MA, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Mina Ryten
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London
| | - John Hardy
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London
| | - Sebastian Guelfi
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London
| | - Daniah Trabzuni
- Department of Genetics and Genomic Medicine, Great Ormond Street Institute of Child Health, University College London, London
| | - Sarah Kim-Hellmuth
- Institute of Translational Genomics, Helmholtz Munich, Neuherberg, Germany
- Department of Pediatrics, Dr. von Hauner Children's Hospital, University Hospital LMU Munich, Munich, Germany
| | - Will Rayner
- Institute of Translational Genomics, Helmholtz Munich, Neuherberg, Germany
| | - Hilary Finucane
- Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, MA, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Hedi Peterson
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
| | - Abayomi Mosaku
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Helen Parkinson
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Kaur Alasoo
- Institute of Computer Science, University of Tartu, Tartu, 51009, Estonia
- Open Targets, South Building, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| |
Collapse
|
28
|
Auwerx C, Sadler MC, Woh T, Reymond A, Kutalik Z, Porcu E. Exploiting the mediating role of the metabolome to unravel transcript-to-phenotype associations. eLife 2023; 12:81097. [PMID: 36891970 PMCID: PMC9998083 DOI: 10.7554/elife.81097] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Accepted: 01/17/2023] [Indexed: 02/17/2023] Open
Abstract
Despite the success of genome-wide association studies (GWASs) in identifying genetic variants associated with complex traits, understanding the mechanisms behind these statistical associations remains challenging. Several methods that integrate methylation, gene expression, and protein quantitative trait loci (QTLs) with GWAS data to determine their causal role in the path from genotype to phenotype have been proposed. Here, we developed and applied a multi-omics Mendelian randomization (MR) framework to study how metabolites mediate the effect of gene expression on complex traits. We identified 216 transcript-metabolite-trait causal triplets involving 26 medically relevant phenotypes. Among these associations, 58% were missed by classical transcriptome-wide MR, which only uses gene expression and GWAS data. This allowed the identification of biologically relevant pathways, such as between ANKH and calcium levels mediated by citrate levels and SLC6A12 and serum creatinine through modulation of the levels of the renal osmolyte betaine. We show that the signals missed by transcriptome-wide MR are found, thanks to the increase in power conferred by integrating multiple omics layer. Simulation analyses show that with larger molecular QTL studies and in case of mediated effects, our multi-omics MR framework outperforms classical MR approaches designed to detect causal relationships between single molecular traits and complex phenotypes.
Collapse
Affiliation(s)
- Chiara Auwerx
- Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland.,University Center for Primary Care and Public Health, Lausanne, Switzerland.,Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| | - Marie C Sadler
- Swiss Institute of Bioinformatics, Lausanne, Switzerland.,University Center for Primary Care and Public Health, Lausanne, Switzerland.,Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| | - Tristan Woh
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| | - Alexandre Reymond
- Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland
| | - Zoltán Kutalik
- Swiss Institute of Bioinformatics, Lausanne, Switzerland.,University Center for Primary Care and Public Health, Lausanne, Switzerland.,Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| | - Eleonora Porcu
- Center for Integrative Genomics, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland.,University Center for Primary Care and Public Health, Lausanne, Switzerland
| |
Collapse
|
29
|
Bulgay C, Kasakolu A, Kazan HH, Mijaica R, Zorba E, Akman O, Bayraktar I, Ekmekci R, Koncagul S, Ulucan K, Semenova EA, Larin AK, Kulemin NA, Generozov EV, Balint L, Badicu G, Ahmetov II, Ergun MA. Exome-Wide Association Study of Competitive Performance in Elite Athletes. Genes (Basel) 2023; 14:genes14030660. [PMID: 36980932 PMCID: PMC10048216 DOI: 10.3390/genes14030660] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 02/27/2023] [Accepted: 03/04/2023] [Indexed: 03/09/2023] Open
Abstract
The aim of the study was to identify genetic variants associated with personal best scores in Turkish track and field athletes and to compare allelic frequencies between sprint/power and endurance athletes and controls using a whole-exome sequencing (WES) approach, followed by replication studies in independent cohorts. The discovery phase involved 60 elite Turkish athletes (31 sprint/power and 29 endurance) and 20 ethnically matched controls. The replication phase involved 1132 individuals (115 elite Russian sprinters, 373 elite Russian endurance athletes (of which 75 athletes were with VO2max measurements), 209 controls, 148 Russian and 287 Finnish individuals with muscle fiber composition and cross-sectional area (CSA) data). None of the single nucleotide polymorphisms (SNPs) reached an exome-wide significance level (p < 2.3 × 10−7) in genotype–phenotype and case–control studies of Turkish athletes. However, of the 53 nominally (p < 0.05) associated SNPs, four functional variants were replicated. The SIRT1 rs41299232 G allele was significantly over-represented in Turkish (p = 0.047) and Russian (p = 0.018) endurance athletes compared to sprint/power athletes and was associated with increased VO2max (p = 0.037) and a greater proportion of slow-twitch muscle fibers (p = 0.035). The NUP210 rs2280084 A allele was significantly over-represented in Turkish (p = 0.044) and Russian (p = 0.012) endurance athletes compared to sprint/power athletes. The TRPM2 rs1785440 G allele was significantly over-represented in Turkish endurance athletes compared to sprint/power athletes (p = 0.034) and was associated with increased VO2max (p = 0.008). The AGRN rs4074992 C allele was significantly over-represented in Turkish sprint/power athletes compared to endurance athletes (p = 0.037) and was associated with a greater CSA of fast-twitch muscle fibers (p = 0.024). In conclusion, we present the first WES study of athletes showing that this approach can be used to identify novel genetic markers associated with exercise- and sport-related phenotypes.
Collapse
Affiliation(s)
- Celal Bulgay
- Sports Science Faculty, Bingol University, 12000 Bingol, Turkey
| | - Anıl Kasakolu
- Faculty of Agriculture, Ankara University, 06000 Ankara, Turkey
| | - Hasan Hüseyin Kazan
- Medical Genetics Department, Faculty of Medicine, Near East University, 1010–1107 Nicosia, Cyprus
- DESAM Institute, Near East University, 1010–1107 Nicosia, Cyprus
| | - Raluca Mijaica
- Department of Physical Education and Special Motricity, Faculty of Physical Education and Mountain Sports, Transilvania University, 500068 Braşov, Romania
- Correspondence:
| | - Erdal Zorba
- Sports Science Faculty, Gazi University, 06560 Ankara, Turkey
| | - Onur Akman
- Sports Science Faculty, Bayburt University, 69000 Bayburt, Turkey
| | - Isık Bayraktar
- Sports Science Faculty, Alanya Alaaddin Keykubat University, 07450 Alanya, Turkey
| | - Rıdvan Ekmekci
- Sports Science Faculty, Pamukkale University, 20160 Denizli, Turkey
| | | | - Korkut Ulucan
- Sports Department of Medical Biology and Genetics, Marmara University, 34722 Istanbul, Turkey
| | - Ekaterina A. Semenova
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
- Research Institute of Physical Culture and Sport, Volga Region State University of Physical Culture, Sport and Tourism, 420138 Kazan, Russia
| | - Andrey K. Larin
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Nikolay A. Kulemin
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Edward V. Generozov
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Lorand Balint
- Department of Physical Education and Special Motricity, Faculty of Physical Education and Mountain Sports, Transilvania University, 500068 Braşov, Romania
| | - Georgian Badicu
- Department of Physical Education and Special Motricity, Faculty of Physical Education and Mountain Sports, Transilvania University, 500068 Braşov, Romania
| | - Ildus I. Ahmetov
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
- Laboratory of Genetics of Aging and Longevity, Kazan State Medical University, 420012 Kazan, Russia
- Department of Physical Education, Plekhanov Russian University of Economics, 115093 Moscow, Russia
- Research Institute for Sport and Exercise Sciences, Liverpool John Moores University, Liverpool L3 5AF, UK
| | - Mehmet Ali Ergun
- Department of Medical Genetics, Faculty of Medicine, Gazi University, 06560 Ankara, Turkey
| |
Collapse
|
30
|
Hemerich D, Smit RAJ, Preuss M, Stalbow L, van der Laan SW, Asselbergs FW, van Setten J, Tragante V. Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes. Sci Rep 2023; 13:3579. [PMID: 36864090 PMCID: PMC9981672 DOI: 10.1038/s41598-023-30369-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 02/21/2023] [Indexed: 03/04/2023] Open
Abstract
Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.
Collapse
Affiliation(s)
- Daiane Hemerich
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Roelof A J Smit
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Michael Preuss
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Lauren Stalbow
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Sander W van der Laan
- Central Diagnostics Laboratory, Division Laboratories, Pharmacy, and Biomedical Genetics, University Medical Center Utrecht, Utrecht University, Utrecht, The Netherlands
| | - Folkert W Asselbergs
- Department of Cardiology, Amsterdam University Medical Centers, University of Amsterdam, Amsterdam, The Netherlands
- Health Data Research UK and Institute of Health Informatics, University College London, London, UK
| | - Jessica van Setten
- Department of Cardiology, UMC Utrecht, Utrecht University, Utrecht, The Netherlands
| | - Vinicius Tragante
- Department of Cardiology, UMC Utrecht, Utrecht University, Utrecht, The Netherlands.
| |
Collapse
|
31
|
Yang G, Mishra M, Perera MA. Multi-Omics Studies in Historically Excluded Populations: The Road to Equity. Clin Pharmacol Ther 2023; 113:541-556. [PMID: 36495075 DOI: 10.1002/cpt.2818] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Accepted: 12/01/2022] [Indexed: 12/14/2022]
Abstract
Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.
Collapse
Affiliation(s)
- Guang Yang
- Department of Pharmacology, Center for Pharmacogenomics, Feinberg School of Medicine, Northwestern University, Chicago, Illinois, USA
| | - Mrinal Mishra
- Department of Pharmacology, Center for Pharmacogenomics, Feinberg School of Medicine, Northwestern University, Chicago, Illinois, USA
| | - Minoli A Perera
- Department of Pharmacology, Center for Pharmacogenomics, Feinberg School of Medicine, Northwestern University, Chicago, Illinois, USA
| |
Collapse
|
32
|
Zhelankin AV, Iulmetova LN, Ahmetov II, Generozov EV, Sharova EI. Diversity and Differential Expression of MicroRNAs in the Human Skeletal Muscle with Distinct Fiber Type Composition. Life (Basel) 2023; 13:659. [PMID: 36983815 PMCID: PMC10056610 DOI: 10.3390/life13030659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Revised: 02/22/2023] [Accepted: 02/23/2023] [Indexed: 03/06/2023] Open
Abstract
The ratio of fast- and slow-twitch fibers in human skeletal muscle is variable and largely determined by genetic factors. In this study, we investigated the contribution of microRNA (miRNA) in skeletal muscle fiber type composition. The study involved biopsy samples of the vastus lateralis muscle from 24 male participants with distinct fiber type ratios. The miRNA study included samples from five endurance athletes and five power athletes with the predominance of slow-twitch (61.6-72.8%) and fast-twitch (69.3-80.7%) fibers, respectively. Total and small RNA were extracted from tissue samples. Total RNA sequencing (N = 24) revealed 352 differentially expressed genes between the groups with the predominance of fast- and slow-twitch muscle fibers. Small RNA sequencing showed upregulation of miR-206, miR-501-3p and miR-185-5p, and downregulation of miR-499a-5p and miR-208-5p in the group of power athletes with fast-twitch fiber predominance. Two miRtronic miRNAs, miR-208b-3p and miR-499a-5p, had strong correlations in expression with their host genes (MYH7 and MYH7B, respectively). Correlations between the expression of miRNAs and their experimentally validated messenger RNA (mRNA) targets were calculated, and 11 miRNA-mRNA interactions with strong negative correlations were identified. Two of them belonged to miR-208b-3p and miR-499a-5p, indicating their regulatory links with the expression of CDKN1A and FOXO4, respectively.
Collapse
Affiliation(s)
- Andrey V. Zhelankin
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Liliia N. Iulmetova
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Ildus I. Ahmetov
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
- Research Institute for Sport and Exercise Sciences, Liverpool John Moores University, Liverpool L3 5AF, UK
| | - Eduard V. Generozov
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Elena I. Sharova
- Department of Molecular Biology and Genetics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| |
Collapse
|
33
|
Taylor HJ, Hung YH, Narisu N, Erdos MR, Kanke M, Yan T, Grenko CM, Swift AJ, Bonnycastle LL, Sethupathy P, Collins FS, Taylor DL. Human pancreatic islet microRNAs implicated in diabetes and related traits by large-scale genetic analysis. Proc Natl Acad Sci U S A 2023; 120:e2206797120. [PMID: 36757889 PMCID: PMC9963967 DOI: 10.1073/pnas.2206797120] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 01/11/2023] [Indexed: 02/10/2023] Open
Abstract
Genetic studies have identified ≥240 loci associated with the risk of type 2 diabetes (T2D), yet most of these loci lie in non-coding regions, masking the underlying molecular mechanisms. Recent studies investigating mRNA expression in human pancreatic islets have yielded important insights into the molecular drivers of normal islet function and T2D pathophysiology. However, similar studies investigating microRNA (miRNA) expression remain limited. Here, we present data from 63 individuals, the largest sequencing-based analysis of miRNA expression in human islets to date. We characterized the genetic regulation of miRNA expression by decomposing the expression of highly heritable miRNAs into cis- and trans-acting genetic components and mapping cis-acting loci associated with miRNA expression [miRNA-expression quantitative trait loci (eQTLs)]. We found i) 84 heritable miRNAs, primarily regulated by trans-acting genetic effects, and ii) 5 miRNA-eQTLs. We also used several different strategies to identify T2D-associated miRNAs. First, we colocalized miRNA-eQTLs with genetic loci associated with T2D and multiple glycemic traits, identifying one miRNA, miR-1908, that shares genetic signals for blood glucose and glycated hemoglobin (HbA1c). Next, we intersected miRNA seed regions and predicted target sites with credible set SNPs associated with T2D and glycemic traits and found 32 miRNAs that may have altered binding and function due to disrupted seed regions. Finally, we performed differential expression analysis and identified 14 miRNAs associated with T2D status-including miR-187-3p, miR-21-5p, miR-668, and miR-199b-5p-and 4 miRNAs associated with a polygenic score for HbA1c levels-miR-216a, miR-25, miR-30a-3p, and miR-30a-5p.
Collapse
Affiliation(s)
- Henry J. Taylor
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
- British Heart Foundation Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, CambridgeCB2 0BB, UK
- Heart and Lung Research Institute, University of Cambridge, CambridgeCB2 0BB, UK
| | - Yu-Han Hung
- Department of Biomedical Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY14853
| | - Narisu Narisu
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - Michael R. Erdos
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - Matthew Kanke
- Department of Biomedical Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY14853
| | - Tingfen Yan
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - Caleb M. Grenko
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - Amy J. Swift
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - Lori L. Bonnycastle
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - Praveen Sethupathy
- Department of Biomedical Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY14853
| | - Francis S. Collins
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| | - D. Leland Taylor
- Center for Precision Health Research, National Human Genome Research Institute, NIH, Bethesda, MD20892
| |
Collapse
|
34
|
Gunasekara CJ, MacKay H, Scott CA, Li S, Laritsky E, Baker MS, Grimm SL, Jun G, Li Y, Chen R, Wiemels JL, Coarfa C, Waterland RA. Systemic interindividual epigenetic variation in humans is associated with transposable elements and under strong genetic control. Genome Biol 2023; 24:2. [PMID: 36631879 PMCID: PMC9835319 DOI: 10.1186/s13059-022-02827-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 12/01/2022] [Indexed: 01/13/2023] Open
Abstract
BACKGROUND Genetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition. RESULTS We use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R2=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease. CONCLUSIONS A focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.
Collapse
Affiliation(s)
- Chathura J. Gunasekara
- grid.508989.50000 0004 6410 7501USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX USA
| | - Harry MacKay
- grid.508989.50000 0004 6410 7501USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX USA
| | - C. Anthony Scott
- grid.508989.50000 0004 6410 7501USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX USA
| | - Shaobo Li
- grid.42505.360000 0001 2156 6853Keck School of Medicine, University of Southern California, Los Angeles, CA USA
| | - Eleonora Laritsky
- grid.508989.50000 0004 6410 7501USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX USA
| | - Maria S. Baker
- grid.508989.50000 0004 6410 7501USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX USA
| | - Sandra L. Grimm
- grid.39382.330000 0001 2160 926XDepartment of Molecular and Cellular Biology, Baylor College of Medicine, Houston, TX USA
| | - Goo Jun
- grid.267308.80000 0000 9206 2401Human Genetics Center, University of Texas Health Science Center at Houston, Houston, TX USA
| | - Yumei Li
- grid.39382.330000 0001 2160 926XDepartment of Molecular & Human Genetics, Baylor College of Medicine, Houston, TX USA
| | - Rui Chen
- grid.39382.330000 0001 2160 926XDepartment of Molecular & Human Genetics, Baylor College of Medicine, Houston, TX USA
| | - Joseph L. Wiemels
- grid.42505.360000 0001 2156 6853Keck School of Medicine, University of Southern California, Los Angeles, CA USA
| | - Cristian Coarfa
- grid.39382.330000 0001 2160 926XDepartment of Molecular and Cellular Biology, Baylor College of Medicine, Houston, TX USA ,grid.39382.330000 0001 2160 926XDan L Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX USA
| | - Robert A. Waterland
- grid.508989.50000 0004 6410 7501USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX USA ,grid.39382.330000 0001 2160 926XDepartment of Molecular & Human Genetics, Baylor College of Medicine, Houston, TX USA
| |
Collapse
|
35
|
Oliva M, Demanelis K, Lu Y, Chernoff M, Jasmine F, Ahsan H, Kibriya MG, Chen LS, Pierce BL. DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits. Nat Genet 2023; 55:112-122. [PMID: 36510025 PMCID: PMC10249665 DOI: 10.1038/s41588-022-01248-z] [Citation(s) in RCA: 53] [Impact Index Per Article: 53.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Accepted: 10/26/2022] [Indexed: 12/14/2022]
Abstract
Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.
Collapse
Affiliation(s)
- Meritxell Oliva
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA.
| | - Kathryn Demanelis
- Department of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Yihao Lu
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA
| | - Meytal Chernoff
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA
| | - Farzana Jasmine
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA
| | - Habibul Ahsan
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Comprehensive Cancer Center, University of Chicago, Chicago, IL, USA
- Department of Medicine, University of Chicago, Chicago, IL, USA
- Institute for Population and Precision Health, University of Chicago, Chicago, IL, USA
| | - Muhammad G Kibriya
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA
- Institute for Population and Precision Health, University of Chicago, Chicago, IL, USA
| | - Lin S Chen
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA.
| | - Brandon L Pierce
- Department of Public Health Sciences, University of Chicago, Chicago, IL, USA.
- Department of Human Genetics, University of Chicago, Chicago, IL, USA.
- Comprehensive Cancer Center, University of Chicago, Chicago, IL, USA.
| |
Collapse
|
36
|
Semenova EA, Zempo H, Miyamoto-Mikami E, Kumagai H, Larin AK, Sultanov RI, Babalyan KA, Zhelankin AV, Tobina T, Shiose K, Kakigi R, Tsuzuki T, Ichinoseki-Sekine N, Kobayashi H, Naito H, Burniston J, Generozov EV, Fuku N, Ahmetov II. Genome-Wide Association Study Identifies CDKN1A as a Novel Locus Associated with Muscle Fiber Composition. Cells 2022; 11:cells11233910. [PMID: 36497168 PMCID: PMC9737696 DOI: 10.3390/cells11233910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Revised: 11/24/2022] [Accepted: 12/01/2022] [Indexed: 12/12/2022] Open
Abstract
Muscle fiber composition is associated with physical performance, with endurance athletes having a high proportion of slow-twitch muscle fibers compared to power athletes. Approximately 45% of muscle fiber composition is heritable, however, single nucleotide polymorphisms (SNP) underlying inter-individual differences in muscle fiber types remain largely unknown. Based on three whole genome SNP datasets, we have shown that the rs236448 A allele located near the cyclin-dependent kinase inhibitor 1A (CDKN1A) gene was associated with an increased proportion of slow-twitch muscle fibers in Russian (n = 151; p = 0.039), Finnish (n = 287; p = 0.03), and Japanese (n = 207; p = 0.008) cohorts (meta-analysis: p = 7.9 × 10−5. Furthermore, the frequency of the rs236448 A allele was significantly higher in Russian (p = 0.045) and Japanese (p = 0.038) elite endurance athletes compared to ethnically matched power athletes. On the contrary, the C allele was associated with a greater proportion of fast-twitch muscle fibers and a predisposition to power sports. CDKN1A participates in cell cycle regulation and is suppressed by the miR-208b, which has a prominent role in the activation of the slow myofiber gene program. Bioinformatic analysis revealed that the rs236448 C allele was associated with increased CDKN1A expression in whole blood (p = 8.5 × 10−15) and with greater appendicular lean mass (p = 1.2 × 10−5), whereas the A allele was associated with longer durations of exercise (p = 0.044) reported amongst the UK Biobank cohort. Furthermore, the expression of CDKN1A increased in response to strength (p < 0.0001) or sprint (p = 0.00035) training. Accordingly, we found that CDKN1A expression is significantly (p = 0.002) higher in the m. vastus lateralis of strength athletes compared to endurance athletes and is positively correlated with the percentage of fast-twitch muscle fibers (p = 0.018). In conclusion, our data suggest that the CDKN1A rs236448 SNP may be implicated in the determination of muscle fiber composition and may affect athletic performance.
Collapse
Affiliation(s)
- Ekaterina A. Semenova
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
- Research Institute of Physical Culture and Sport, Volga Region State University of Physical Culture, Sport and Tourism, 420138 Kazan, Russia
| | - Hirofumi Zempo
- Faculty of Health and Nutrition, Tokyo Seiei College, Tokyo 124-0025, Japan
| | - Eri Miyamoto-Mikami
- Graduate School of Health and Sports Science, Juntendo University, Chiba 270-1695, Japan
| | - Hiroshi Kumagai
- Graduate School of Health and Sports Science, Juntendo University, Chiba 270-1695, Japan
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089, USA
| | - Andrey K. Larin
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Rinat I. Sultanov
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Konstantin A. Babalyan
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Andrey V. Zhelankin
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Takuro Tobina
- Faculty of Nursing and Nutrition, University of Nagasaki, Nagasaki 851-2195, Japan
| | - Keisuke Shiose
- Faculty of Education, University of Miyazaki, Miyazaki 889-2192, Japan
| | - Ryo Kakigi
- Faculty of Management & Information Science, Josai International University, Chiba 283-8555, Japan
| | | | - Noriko Ichinoseki-Sekine
- Graduate School of Health and Sports Science, Juntendo University, Chiba 270-1695, Japan
- Faculty of Liberal Arts, The Open University of Japan, Chiba 261-8586, Japan
| | - Hiroyuki Kobayashi
- Department of General Medicine, Mito Medical Center, Tsukuba University Hospital, Ibaraki 310-0015, Japan
| | - Hisashi Naito
- Graduate School of Health and Sports Science, Juntendo University, Chiba 270-1695, Japan
| | - Jatin Burniston
- Research Institute for Sport and Exercise Sciences, Liverpool John Moores University, Liverpool L3 5AF, UK
| | - Edward V. Generozov
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
| | - Noriyuki Fuku
- Graduate School of Health and Sports Science, Juntendo University, Chiba 270-1695, Japan
| | - Ildus I. Ahmetov
- Department of Molecular Biology and Genetics, Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, 119435 Moscow, Russia
- Research Institute for Sport and Exercise Sciences, Liverpool John Moores University, Liverpool L3 5AF, UK
- Department of Physical Education, Plekhanov Russian University of Economics, 115093 Moscow, Russia
- Laboratory of Genetics of Aging and Longevity, Kazan State Medical University, 420012 Kazan, Russia
- Correspondence:
| |
Collapse
|
37
|
Elucidating the genetic architecture of DNA methylation to identify promising molecular mechanisms of disease. Sci Rep 2022; 12:19564. [PMID: 36380121 PMCID: PMC9664436 DOI: 10.1038/s41598-022-24100-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Accepted: 11/10/2022] [Indexed: 11/16/2022] Open
Abstract
DNA methylation commonly occurs at cytosine-phosphate-guanine sites (CpGs) that can serve as biomarkers for many diseases. We analyzed whole genome sequencing data to identify DNA methylation quantitative trait loci (mQTLs) in 4126 Framingham Heart Study participants. Our mQTL mapping identified 94,362,817 cis-mQTLvariant-CpG pairs (for 210,156 unique autosomal CpGs) at P < 1e-7 and 33,572,145 trans-mQTL variant-CpG pairs (for 213,606 unique autosomal CpGs) at P < 1e-14. Using cis-mQTL variants for 1258 CpGs associated with seven cardiovascular disease (CVD) risk factors, we found 104 unique CpGs that colocalized with at least one CVD trait. For example, cg11554650 (PPP1R18) colocalized with type 2 diabetes, and was driven by a single nucleotide polymorphism (rs2516396). We performed Mendelian randomization (MR) analysis and demonstrated 58 putatively causal relations of CVD risk factor-associated CpGs to one or more risk factors (e.g., cg05337441 [APOB] with LDL; MR P = 1.2e-99, and 17 causal associations with coronary artery disease (e.g. cg08129017 [SREBF1] with coronary artery disease; MR P = 5e-13). We also showed that three CpGs, e.g., cg14893161 (PM20D1), are putatively causally associated with COVID-19 severity. To assist in future analyses of the role of DNA methylation in disease pathogenesis, we have posted a comprehensive summary data set in the National Heart, Lung, and Blood Institute's BioData Catalyst.
Collapse
|
38
|
Borie R, Cardwell J, Konigsberg IR, Moore CM, Zhang W, Sasse SK, Gally F, Dobrinskikh E, Walts A, Powers J, Brancato J, Rojas M, Wolters PJ, Brown KK, Blackwell TS, Nakanishi T, Richards JB, Gerber AN, Fingerlin TE, Sachs N, Pulit SL, Zappala Z, Schwartz DA, Yang IV. Colocalization of Gene Expression and DNA Methylation with Genetic Risk Variants Supports Functional Roles of MUC5B and DSP in Idiopathic Pulmonary Fibrosis. Am J Respir Crit Care Med 2022; 206:1259-1270. [PMID: 35816432 PMCID: PMC9746850 DOI: 10.1164/rccm.202110-2308oc] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Accepted: 07/05/2022] [Indexed: 11/16/2022] Open
Abstract
Rationale: Common genetic variants have been associated with idiopathic pulmonary fibrosis (IPF). Objectives: To determine functional relevance of the 10 IPF-associated common genetic variants we previously identified. Methods: We performed expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL) mapping, followed by co-localization of eQTL and mQTL with genetic association signals and functional validation by luciferase reporter assays. Illumina multi-ethnic genotyping arrays, mRNA sequencing, and Illumina 850k methylation arrays were performed on lung tissue of participants with IPF (234 RNA and 345 DNA samples) and non-diseased controls (188 RNA and 202 DNA samples). Measurements and Main Results: Focusing on genetic variants within 10 IPF-associated genetic loci, we identified 27 eQTLs in controls and 24 eQTLs in cases (false-discovery-rate-adjusted P < 0.05). Among these signals, we identified associations of lead variants rs35705950 with expression of MUC5B and rs2076295 with expression of DSP in both cases and controls. mQTL analysis identified CpGs in gene bodies of MUC5B (cg17589883) and DSP (cg08964675) associated with the lead variants in these two loci. We also demonstrated strong co-localization of eQTL/mQTL and genetic signal in MUC5B (rs35705950) and DSP (rs2076295). Functional validation of the mQTL in MUC5B using luciferase reporter assays demonstrates that the CpG resides within a putative internal repressor element. Conclusions: We have established a relationship of the common IPF genetic risk variants rs35705950 and rs2076295 with respective changes in MUC5B and DSP expression and methylation. These results provide additional evidence that both MUC5B and DSP are involved in the etiology of IPF.
Collapse
Affiliation(s)
| | | | | | - Camille M. Moore
- Department of Biostatistics and Bioinformatics and
- Center for Genes, Environment, and Health
| | | | | | - Fabienne Gally
- Department of Medicine
- Department of Immunology and Genomic Medicine, National Jewish Health, Denver, Colorado
| | | | | | | | | | - Mauricio Rojas
- Department of Internal Medicine, Ohio State College of Medicine, The Ohio State University, Columbus, Ohio
| | - Paul J. Wolters
- Department of Medicine, University of California, San Francisco, California
| | | | - Timothy S. Blackwell
- Department of Medicine, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Tomoko Nakanishi
- Department of Human Genetics, Lady Davis Institute, Jewish General Hospital, McGill University, Montréal, Canada
| | - J. Brent Richards
- Department of Human Genetics, Lady Davis Institute, Jewish General Hospital, McGill University, Montréal, Canada
| | - Anthony N. Gerber
- Department of Medicine
- Department of Medicine, and
- Department of Immunology and Genomic Medicine, National Jewish Health, Denver, Colorado
| | - Tasha E. Fingerlin
- Department of Biostatistics and Bioinformatics and
- Center for Genes, Environment, and Health
- Department of Immunology and Genomic Medicine, National Jewish Health, Denver, Colorado
| | - Norman Sachs
- Cell Biology, Vertex Pharmaceuticals, San Diego, California; and
| | - Sara L. Pulit
- Computational Genomics, Vertex Pharmaceuticals, Boston, Massachusetts
| | - Zachary Zappala
- Computational Genomics, Vertex Pharmaceuticals, Boston, Massachusetts
| | - David A. Schwartz
- Department of Medicine
- Department of Microbiology and Immunology, University of Colorado Anschutz Medical Campus; Aurora, Colorado
| | - Ivana V. Yang
- Department of Medicine
- Department of Epidemiology, Colorado School of Public Health, Aurora, Colorado
| |
Collapse
|
39
|
Huang D, Feng X, Yang H, Wang J, Zhang W, Fan X, Dong X, Chen K, Yu Y, Ma X, Yi X, Li M. QTLbase2: an enhanced catalog of human quantitative trait loci on extensive molecular phenotypes. Nucleic Acids Res 2022; 51:D1122-D1128. [PMID: 36330927 PMCID: PMC9825467 DOI: 10.1093/nar/gkac1020] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Revised: 10/17/2022] [Accepted: 10/21/2022] [Indexed: 11/06/2022] Open
Abstract
Deciphering the fine-scale molecular mechanisms that shape the genetic effects at disease-associated loci from genome-wide association studies (GWAS) remains challenging. The key avenue is to identify the essential molecular phenotypes that mediate the causal variant and disease under particular biological conditions. Therefore, integrating GWAS signals with context-specific quantitative trait loci (QTLs) (such as different tissue/cell types, disease states, and perturbations) from extensive molecular phenotypes would present important strategies for full understanding of disease genetics. Via persistent curation and systematic data processing of large-scale human molecular trait QTLs (xQTLs), we updated our previous QTLbase database (now QTLbase2, http://mulinlab.org/qtlbase) to comprehensively analyze and visualize context-specific QTLs across 22 molecular phenotypes and over 95 tissue/cell types. Overall, the resource features the following major updates and novel functions: (i) 960 more genome-wide QTL summary statistics from 146 independent studies; (ii) new data for 10 previously uncompiled QTL types; (iii) variant query scope expanded to fit 195 QTL datasets based on whole-genome sequencing; (iv) supports filtering and comparison of QTLs for different biological conditions, such as stimulation types and disease states; (v) a new linkage disequilibrium viewer to facilitate variant prioritization across tissue/cell types and QTL types.
Collapse
Affiliation(s)
| | | | - Hongxi Yang
- Department of Pharmacology, The Province and Ministry Co-sponsored Collaborative Innovation Center for Medical Epigenetics, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
| | - Jianhua Wang
- Department of Pharmacology, The Province and Ministry Co-sponsored Collaborative Innovation Center for Medical Epigenetics, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
| | - Wenwen Zhang
- Department of Pharmacology, The Province and Ministry Co-sponsored Collaborative Innovation Center for Medical Epigenetics, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
| | - Xutong Fan
- Department of Pharmacology, The Province and Ministry Co-sponsored Collaborative Innovation Center for Medical Epigenetics, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
| | - Xiaobao Dong
- Department of Bioinformatics, Tianjin Key Laboratory of Molecular Cancer Epidemiology, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin, China
| | - Kexin Chen
- Department of Bioinformatics, Tianjin Key Laboratory of Molecular Cancer Epidemiology, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin, China
| | - Ying Yu
- Department of Pharmacology, The Province and Ministry Co-sponsored Collaborative Innovation Center for Medical Epigenetics, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
| | - Xin Ma
- Correspondence may also be addressed to Xin Ma.
| | - Xianfu Yi
- Correspondence may also be addressed to Xianfu Yi.
| | - Mulin Jun Li
- To whom correspondence should be addressed. Tel: +86 22 83336668; Fax: +86 22 83336668;
| |
Collapse
|
40
|
Bykova M, Hou Y, Eng C, Cheng F. Quantitative trait locus (xQTL) approaches identify risk genes and drug targets from human non-coding genomes. Hum Mol Genet 2022; 31:R105-R113. [PMID: 36018824 PMCID: PMC9989738 DOI: 10.1093/hmg/ddac208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 08/18/2022] [Accepted: 08/19/2022] [Indexed: 11/13/2022] Open
Abstract
Advances and reduction of costs in various sequencing technologies allow for a closer look at variations present in the non-coding regions of the human genome. Correlating non-coding variants with large-scale multi-omic data holds the promise not only of a better understanding of likely causal connections between non-coding DNA and expression of traits but also identifying potential disease-modifying medicines. Genome-phenome association studies have created large datasets of DNA variants that are associated with multiple traits or diseases, such as Alzheimer's disease; yet, the functional consequences of variants, in particular of non-coding variants, remain largely unknown. Recent advances in functional genomics and computational approaches have led to the identification of potential roles of DNA variants, such as various quantitative trait locus (xQTL) techniques. Multi-omic assays and analytic approaches toward xQTL have identified links between genetic loci and human transcriptomic, epigenomic, proteomic and metabolomic data. In this review, we first discuss the recent development of xQTL from multi-omic findings. We then highlight multimodal analysis of xQTL and genetic data for identification of risk genes and drug targets using Alzheimer's disease as an example. We finally discuss challenges and future research directions (e.g. artificial intelligence) for annotation of non-coding variants in complex diseases.
Collapse
Affiliation(s)
- Marina Bykova
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Yuan Hou
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Charis Eng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
- Department of Molecular Medicine, Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA
| | - Feixiong Cheng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
- Department of Molecular Medicine, Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA
| |
Collapse
|
41
|
Wang K, Dai R, Xia Y, Tian J, Jiao C, Mikhailova T, Zhang C, Chen C, Liu C. Spatiotemporal specificity of correlated DNA methylation and gene expression pairs across different human tissues and stages of brain development. Epigenetics 2022; 17:1110-1127. [PMID: 34652256 PMCID: PMC9543113 DOI: 10.1080/15592294.2021.1993607] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
DNA methylation (DNAm) that occurs on promoter regions is primarily considered to repress gene expression. Previous studies indicated that DNAm could also show positive correlations with gene expression. Both DNAm and gene expression profiles are known to be tissue- and development-specific. This study aims to investigate how DNAm and gene expression are coordinated across different human tissues and developmental stages, as well as the biological significance of such correlations. By analyzing 2,239 samples with both DNAm and gene expression data in the same human subjects obtained from six published datasets, we evaluated the correlations between gene and CpG pairs (GCPs) at cis-regions and compared significantly correlated GCPs (cGCPs) across different tissues and brains at different age groups. A total of 37,363 cGCPs was identified in the six datasets; approximately 38% of the cGCPs were positively correlated. The majority (>90%) of cGCPs was tissue- or development-specific. We also observed that the correlation direction can be opposite in different tissues and ages. Further analysis highlights the importance of cGCPs for their cellular functions and potential roles in complex traits and human diseases. For instance, the early developmental brain possessed a highly unique set of cGCPs that were associated with neurogenesis and psychiatric disorders. By assessing the epigenetic factors involved in cGCPs, we discovered novel regulatory mechanisms of positive cGCPs distinct from negative cGCPs, which were related to multiple factors, such as H3K27me3, CTCF, and JARD2. The catalogue of cGCPs compiled can be used to guide functional interpretation of genetic and epigenetic studies.
Collapse
Affiliation(s)
- Kangli Wang
- Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, and Department of Psychiatry, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China
| | - Rujia Dai
- Department of Psychiatry, State University of New York Upstate Medical University, Syracuse, NY, USA
| | - Yan Xia
- Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, and Department of Psychiatry, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China,Department of Psychiatry, State University of New York Upstate Medical University, Syracuse, NY, USA
| | - Jianghua Tian
- Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, and Department of Psychiatry, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China
| | - Chuan Jiao
- Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, and Department of Psychiatry, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China,Department of Psychiatry, State University of New York Upstate Medical University, Syracuse, NY, USA
| | - Tatiana Mikhailova
- Department of Psychiatry, State University of New York Upstate Medical University, Syracuse, NY, USA
| | - Chunling Zhang
- Department of Neuroscience and Physiology, State University of New York Upstate Medical University, Syracuse, Ny, USA
| | - Chao Chen
- Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, and Department of Psychiatry, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China,National Clinical Research Center on Mental Disorders, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, Hunan, China
| | - Chunyu Liu
- Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, and Department of Psychiatry, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China,Department of Psychiatry, State University of New York Upstate Medical University, Syracuse, NY, USA,Department of Neuroscience and Physiology, State University of New York Upstate Medical University, Syracuse, Ny, USA,CONTACT Chunyu Liu Department of Psychiatry, State University of New York Upstate Medical University, Syracuse, NY13210, USA; Chao Chen Center for Medical Genetics, Central South University, Changsha, Hunan 410005, China
| |
Collapse
|
42
|
Sleiman MB, Roy S, Gao AW, Sadler MC, von Alvensleben GVG, Li H, Sen S, Harrison DE, Nelson JF, Strong R, Miller RA, Kutalik Z, Williams RW, Auwerx J. Sex- and age-dependent genetics of longevity in a heterogeneous mouse population. Science 2022; 377:eabo3191. [PMID: 36173858 PMCID: PMC9905652 DOI: 10.1126/science.abo3191] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
DNA variants that modulate life span provide insight into determinants of health, disease, and aging. Through analyses in the UM-HET3 mice of the Interventions Testing Program (ITP), we detected a sex-independent quantitative trait locus (QTL) on chromosome 12 and identified sex-specific QTLs, some of which we detected only in older mice. Similar relations between life history and longevity were uncovered in mice and humans, underscoring the importance of early access to nutrients and early growth. We identified common age- and sex-specific genetic effects on gene expression that we integrated with model organism and human data to create a hypothesis-building interactive resource of prioritized longevity and body weight genes. Finally, we validated Hipk1, Ddost, Hspg2, Fgd6, and Pdk1 as conserved longevity genes using Caenorhabditis elegans life-span experiments.
Collapse
Affiliation(s)
- Maroun Bou Sleiman
- Laboratory of Integrative Systems Physiology, Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Suheeta Roy
- Department of Genetics, Genomics and Informatics, University of Tennessee Health Science Center (UTHSC), Memphis, TN 38163, USA
| | - Arwen W. Gao
- Laboratory of Integrative Systems Physiology, Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Marie C. Sadler
- Institute of Primary Care and Public Health (Unisante), University of Lausanne, Lausanne 1011, Switzerland
- Swiss Institute of Bioinformatics, Lausanne 1015, Switzerland
- Department of Computational Biology, University of Lausanne, Lausanne 1015, Switzerland
| | - Giacomo V. G. von Alvensleben
- Laboratory of Integrative Systems Physiology, Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Hao Li
- Laboratory of Integrative Systems Physiology, Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Saunak Sen
- Department of Preventive Medicine, University of Tennessee Health Science Center, Memphis, TN 38163, USA
| | | | - James F. Nelson
- Barshop Center for Longevity Studies at the University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA
| | - Randy Strong
- Barshop Center for Longevity Studies at the University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA
- South Texas Veterans Healthcare System, San Antonio, TX 78229, USA
| | - Richard A. Miller
- Department of Pathology, University of Michigan Geriatrics Center, Ann Arbor, MI 48109-2200, USA
| | - Zoltán Kutalik
- Institute of Primary Care and Public Health (Unisante), University of Lausanne, Lausanne 1011, Switzerland
- Swiss Institute of Bioinformatics, Lausanne 1015, Switzerland
- Department of Computational Biology, University of Lausanne, Lausanne 1015, Switzerland
| | - Robert W. Williams
- Department of Genetics, Genomics and Informatics, University of Tennessee Health Science Center (UTHSC), Memphis, TN 38163, USA
| | - Johan Auwerx
- Laboratory of Integrative Systems Physiology, Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| |
Collapse
|
43
|
Profiling Analysis of N6-Methyladenosine mRNA Methylation Reveals Differential m6A Patterns during the Embryonic Skeletal Muscle Development of Ducks. Animals (Basel) 2022; 12:ani12192593. [PMID: 36230334 PMCID: PMC9559603 DOI: 10.3390/ani12192593] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 09/03/2022] [Accepted: 09/23/2022] [Indexed: 12/03/2022] Open
Abstract
Simple Summary Recent studies show that N6-methyladenosine (m6A) modification, the most common RNA chemical modification, influences the modification, processing, transport, and translation of RNA. N6-methyladenosine is an epigenetic modification that influences skeletal myogenesis and skeletal muscle development. However, the N6-methyladenosine modification profile and its function during poultry muscle development is unclear, and there is only one report about m6A modification in ducks, which focuses on duck hepatitis A virus infection. Here, we compared the m6A modification profiles between E13 (embryonic day 13) and E19 (embryonic day 19) in duck breast muscle differentiation using MeRIP-seq, and evaluated the expression profile of the methyl transferase METTL14 and its cofactors during breast muscle development. This is the first study of N6-methyladenosine modification patterns in duck muscle tissue. The current study not only elucidates the regulation mechanisms of duck skeletal muscle development, but also lays the groundwork for studying the role of RNA modification in poultry muscle development. Abstract N6-Methyladenosine is a reversible epigenetic modification that influences muscle development. However, the m6A modification profile during poultry skeletal muscle development is poorly understood. Here, we utilized m6A-specific methylated RNA immunoprecipitation sequencing to identify m6A sites during two stages of breast muscle development in ducks: embryonic days 13 (E13) and E19. MeRIP-seq detected 19,024 and 18,081 m6A peaks in the E13 and E19 groups, respectively. Similarly to m6A distribution in mammalian transcripts, our results revealed GGACU as the main m6A motif in duck breast muscle; they also revealed that m6A peaks are mainly enriched near the stop codons. In addition, motif sequence analysis and gene expression analysis demonstrated that m6A modification in duck embryo skeletal muscles may be mediated by the methyltransferase-like 14. GO and KEGG analysis showed that m6A peaks containing genes at E19 were mainly enriched in muscle-differentiation- and muscle-growth-related pathways, whereas m6A peaks containing genes in E13 were mainly enriched in embryonic development and cell proliferation pathways. Combined analysis of MeRIP-seq and RNA-seq showed that the mRNA expression may be affected by m6A modification. Moreover, qRT-PCR analysis of the expression of METTL14 and its cofactors (WTAP, ZC3H13, RBM15 and VIRMA) during duck embryonic skeletal muscle development in breast and leg muscle samples revealed a significant downward trend as the developmental age progressed. Our results demonstrated that m6A mRNA methylation modifications control muscle development in ducks. This is the first study of m6A modification patterns in duck muscle tissue development, and it lays the foundation for the study of the effects of RNA modification on poultry skeletal muscle development.
Collapse
|
44
|
Fan Y, Ren C, Deng K, Zhang Z, Li J, Deng M, Zhang Y, Wang F. The regulation of LncRNA GTL2 expression by DNA methylation during sheep skeletal muscle development. Genomics 2022; 114:110453. [PMID: 36030023 DOI: 10.1016/j.ygeno.2022.110453] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2022] [Revised: 08/08/2022] [Accepted: 08/16/2022] [Indexed: 11/27/2022]
Abstract
DNA methylation has crucial roles in regulating the expression of genes involved in skeletal muscle development. However, the DNA methylation pattern of lncRNA during sheep skeletal muscle development remains unclear. This study investigated previous WGBS and LncRNA data in skeletal muscle of sheep (fetus and adult). We then focused on LncRNA GTL2, which is differentially expressed in skeletal muscle and has multiple DMRs. We found that the expression level of GTL2 decreased with age. GTL2 DMRs methylation levels were significantly higher in adult muscle than in fetal muscle. After 5AZA treatment, GTL2 expression was significantly increased in a dose-dependent manner.The dCas9-DNMT3A-sgRNA significantly reduced the expression level of GTL2 in cells, but increased GTL2 DMR methylation levels. The above studies indicate that dCas9-DNMT3A can effectively increase the methylation level in the DMR region of GTL2, the expression level of GTL2 is regulated by DNA methylation during muscle development.
Collapse
Affiliation(s)
- Yixuan Fan
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China
| | - Caifang Ren
- School of Medicine, Jiangsu University, Zhenjiang, Jiangsu Province, China
| | - Kaiping Deng
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China
| | - Zhen Zhang
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China
| | - Juan Li
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China
| | - Mingtian Deng
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China
| | - Yanli Zhang
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China
| | - Feng Wang
- Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu Province, China.
| |
Collapse
|
45
|
Bankier S, Michoel T. eQTLs as causal instruments for the reconstruction of hormone linked gene networks. Front Endocrinol (Lausanne) 2022; 13:949061. [PMID: 36060942 PMCID: PMC9428692 DOI: 10.3389/fendo.2022.949061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/20/2022] [Accepted: 07/25/2022] [Indexed: 11/17/2022] Open
Abstract
Hormones act within in highly dynamic systems and much of the phenotypic response to variation in hormone levels is mediated by changes in gene expression. The increase in the number and power of large genetic association studies has led to the identification of hormone linked genetic variants. However, the biological mechanisms underpinning the majority of these loci are poorly understood. The advent of affordable, high throughput next generation sequencing and readily available transcriptomic databases has shown that many of these genetic variants also associate with variation in gene expression levels as expression Quantitative Trait Loci (eQTLs). In addition to further dissecting complex genetic variation, eQTLs have been applied as tools for causal inference. Many hormone networks are driven by transcription factors, and many of these genes can be linked to eQTLs. In this mini-review, we demonstrate how causal inference and gene networks can be used to describe the impact of hormone linked genetic variation upon the transcriptome within an endocrinology context.
Collapse
Affiliation(s)
- Sean Bankier
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | | |
Collapse
|
46
|
Damer A, El Meniawy S, McPherson R, Wells G, Harper ME, Dent R. Association of muscle fiber type with measures of obesity: A systematic review. Obes Rev 2022; 23:e13444. [PMID: 35293095 DOI: 10.1111/obr.13444] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Revised: 02/22/2022] [Accepted: 02/23/2022] [Indexed: 12/18/2022]
Abstract
Obesity derives from an extended period of positive energy imbalance due to a complex interplay of environmental and biological factors. Muscle fiber type and physiology have been hypothesized to affect metabolism and energy expenditure and thus to affect an individual's propensity to gain weight. However, there have been conflicting reports regarding a relationship between muscle fiber type and obesity. Here, we systematically reviewed literature investigating this topic from PubMed, Web of Science, and EMBASE. Of these, 32 articles were included in our final review and analysis. Most studies (22/32) reported a significant relationship between muscle fiber-type proportion and a measure of obesity. Overall, there was a significant negative relationship between the proportion of type I fibers and body mass index (BMI) and a significant positive relationship between the proportion of type IIX fibers and BMI. Moreover, between-group comparisons indicated a greater prevalence of type IIX fibers and a lower prevalence of type I fibers in patients living with obesity relative to lean individuals. These significant relationships were confirmed in a meta-analysis of these data. The causal nature of these associations remains to be evaluated.
Collapse
Affiliation(s)
- Alameen Damer
- Faculty of Medicine, University of Ottawa, Ottawa, Ontario, Canada
| | | | - Ruth McPherson
- Atherogenomics Laboratory, University of Ottawa Heart Institute, Ottawa, Ontario, Canada
| | - George Wells
- School of Epidemiology and Public Health, Faculty of Medicine, University of Ottawa, Ottawa, Ontario, Canada.,Cardiovascular Research Methods Centre, University of Ottawa Heart Institute, Ottawa, Ontario, Canada
| | - Mary-Ellen Harper
- Department of Biochemistry, Microbiology and Immunology, Faculty of Medicine, Ottawa Institute of Systems Biology, University of Ottawa, Ottawa, Ontario, Canada
| | - Robert Dent
- Department of Medicine, Division of Endocrinology, University of Ottawa, Ottawa, Ontario, Canada
| |
Collapse
|
47
|
Robinson N, Casement J, Gunter MJ, Huybrechts I, Agudo A, Barranco MR, Eichelmann F, Johnson T, Kaaks R, Pala V, Panico S, Sandanger TM, Schultze MB, Travis RC, Tumino R, Vineis P, Weiderpass E, Skinner R, Sharp L, McKay JA, Strathdee G. Anti-cancer therapy is associated with long-term epigenomic changes in childhood cancer survivors. Br J Cancer 2022; 127:288-300. [PMID: 35354948 PMCID: PMC9296636 DOI: 10.1038/s41416-022-01792-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 02/24/2022] [Accepted: 03/11/2022] [Indexed: 11/09/2022] Open
Abstract
BACKGROUND Childhood cancer survivors (CCS) exhibit significantly increased chronic diseases and premature death. Abnormalities in DNA methylation are associated with development of chronic diseases and reduced life expectancy. We investigated the hypothesis that anti-cancer treatments are associated with long-term DNA methylation changes that could be key drivers of adverse late health effects. METHODS Genome-wide DNA methylation was assessed using MethylationEPIC arrays in paired samples (before/after therapy) from 32 childhood cancer patients. Separately, methylation was determined in 32 samples from different adult CCS (mean 22-years post-diagnosis) and compared with cancer-free controls (n = 284). RESULTS Widespread DNA methylation changes were identified post-treatment in childhood cancer patients, including 146 differentially methylated regions (DMRs), which were consistently altered in the 32 post-treatment samples. Analysis of adult CCS identified matching methylation changes at 107/146 of the DMRs, suggesting potential long-term retention of post-therapy changes. Adult survivors also exhibited epigenetic age acceleration, independent of DMR methylation. Furthermore, altered methylation at the DUSP6 DMR was significantly associated with early mortality, suggesting altered methylation may be prognostic for some late adverse health effects in CCS. CONCLUSIONS These novel methylation changes could serve as biomarkers for assessing normal cell toxicity in ongoing treatments and predicting long-term health outcomes in CCS.
Collapse
Affiliation(s)
- Natassia Robinson
- Newcastle University Centre for Cancer, Biosciences Institute, Newcastle University, Newcastle upon Tyne, UK
| | - John Casement
- Bioinformatic Support Unit, Newcastle University, Newcastle upon Tyne, UK
| | - Marc J Gunter
- Section of Nutrition and Metabolism, IARC, Lyon, France
| | | | - Antonio Agudo
- Unit of Nutrition and Cancer, Catalan Institute of Oncology - ICO, Nutrition and Cancer Group, Bellvitge Biomedical Research Institute - IDIBELL, L'Hospitalet de Llobregat, Barcelona, 08908, Spain
| | | | - Fabian Eichelmann
- German Institute of Human Nutrition Potsdam-Rehbrücke, Heidelberg, Germany
| | - Theron Johnson
- German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Rudolf Kaaks
- German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Valeria Pala
- Fondazione IRCCS Istituto Nazionale dei Tumori di Milano, Milan, Italy
| | - Salvatore Panico
- Dipartimento di Medicina Clinica e Chirurgia, Federico II University, Naples, Italy
| | | | | | - Ruth C Travis
- Cancer Epidemiology Unit, University of Oxford, Oxford, UK
| | - Rosario Tumino
- Cancer Registry and Histopathology Department Azienda Sanitaria Provinciale (ASP), Ragusa, Italy
| | - Paolo Vineis
- MRC Centre for Environment and Health, Imperial College London, London, UK
| | | | - Roderick Skinner
- Newcastle University Centre for Cancer, Translational and Clinical Research Institute, Newcastle University, Newcastle upon Tyne, UK
- Department of Paediatric and Adolescent Haematology/Oncology, Great North Children's Hospital, Newcastle, UK
| | - Linda Sharp
- Newcastle University Centre for Cancer, Population Health Sciences Institute, Newcastle upon Tyne, UK
| | - Jill A McKay
- Department of Applied Sciences, Northumbria University, Newcastle upon Tyne, UK
| | - Gordon Strathdee
- Newcastle University Centre for Cancer, Biosciences Institute, Newcastle University, Newcastle upon Tyne, UK.
| |
Collapse
|
48
|
Ling C, Bacos K, Rönn T. Epigenetics of type 2 diabetes mellitus and weight change - a tool for precision medicine? Nat Rev Endocrinol 2022; 18:433-448. [PMID: 35513492 DOI: 10.1038/s41574-022-00671-w] [Citation(s) in RCA: 40] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 03/29/2022] [Indexed: 12/12/2022]
Abstract
Pioneering studies performed over the past few decades demonstrate links between epigenetics and type 2 diabetes mellitus (T2DM), the metabolic disorder with the most rapidly increasing prevalence in the world. Importantly, these studies identified epigenetic modifications, including altered DNA methylation, in pancreatic islets, adipose tissue, skeletal muscle and the liver from individuals with T2DM. As non-genetic factors that affect the risk of T2DM, such as obesity, unhealthy diet, physical inactivity, ageing and the intrauterine environment, have been associated with epigenetic modifications in healthy individuals, epigenetics probably also contributes to T2DM development. In addition, genetic factors associated with T2DM and obesity affect the epigenome in human tissues. Notably, causal mediation analyses found DNA methylation to be a potential mediator of genetic associations with metabolic traits and disease. In the past few years, translational studies have identified blood-based epigenetic markers that might be further developed and used for precision medicine to help patients with T2DM receive optimal therapy and to identify patients at risk of complications. This Review focuses on epigenetic mechanisms in the development of T2DM and the regulation of body weight in humans, with a special focus on precision medicine.
Collapse
Affiliation(s)
- Charlotte Ling
- Epigenetics and Diabetes Unit, Department of Clinical Sciences, Lund University Diabetes Centre, Lund University, Scania University Hospital, Malmö, Sweden.
| | - Karl Bacos
- Epigenetics and Diabetes Unit, Department of Clinical Sciences, Lund University Diabetes Centre, Lund University, Scania University Hospital, Malmö, Sweden
| | - Tina Rönn
- Epigenetics and Diabetes Unit, Department of Clinical Sciences, Lund University Diabetes Centre, Lund University, Scania University Hospital, Malmö, Sweden
| |
Collapse
|
49
|
Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease. Nat Genet 2022; 54:950-962. [PMID: 35710981 DOI: 10.1038/s41588-022-01097-w] [Citation(s) in RCA: 75] [Impact Index Per Article: 37.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Accepted: 05/09/2022] [Indexed: 12/29/2022]
Abstract
More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.
Collapse
|
50
|
Lancaster EE, Vladimirov VI, Riley BP, Landry JW, Roberson-Nay R, York TP. LARGE-SCALE INTEGRATION OF DNA METHYLATION AND GENE EXPRESSION ARRAY PLATFORMS IDENTIFIES BOTH cis AND trans RELATIONSHIPS. Epigenetics 2022; 17:1753-1773. [PMID: 35608069 PMCID: PMC9621057 DOI: 10.1080/15592294.2022.2079293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Although epigenome-wide association studies (EWAS) have been successful in identifying DNA methylation (DNAm) patterns associated with disease states, any further characterization of etiologic mechanisms underlying disease remains elusive. This knowledge gap does not originate from a lack of DNAm–trait associations, but rather stems from study design issues that affect the interpretability of EWAS results. Despite known limitations in predicting the function of a particular CpG site, most EWAS maintain the broad assumption that altered DNAm results in a concomitant change of transcription at the most proximal gene. This study integrated DNAm and gene expression (GE) measurements in two cohorts, the Adolescent and Young Adult Twin Study (AYATS) and the Pregnancy, Race, Environment, Genes (PREG) study, to improve the understanding of epigenomic regulatory mechanisms. CpG sites associated with GE in cis were enriched in areas of transcription factor binding and areas of intermediate-to-low CpG density. CpG sites associated with trans GE were also enriched in areas of known regulatory significance, including enhancer regions. These results highlight issues with restricting DNAm-transcript annotations to small genomic intervals and question the validity of assuming a cis DNAm–GE pathway. Based on these findings, the interpretation of EWAS results is limited in studies without multi-omic support and further research should identify genomic regions in which GE-associated DNAm is overrepresented. An in-depth characterization of GE-associated CpG sites could improve predictions of the downstream functional impact of altered DNAm and inform best practices for interpreting DNAm–trait associations generated by EWAS.
Collapse
Affiliation(s)
- Eva E Lancaster
- Department of Psychiatry, Virginia Commonwealth University, Richmond, VA 23220
| | | | - Brien P Riley
- Department of Psychiatry, Virginia Commonwealth University, Richmond, VA 23220.,Department of Human and Molecular Genetics, Virginia Commonwealth University, Richmond, VA 23220
| | - Joseph W Landry
- Department of Human and Molecular Genetics, Virginia Commonwealth University, Richmond, VA 23220
| | - Roxann Roberson-Nay
- Department of Psychiatry, Virginia Commonwealth University, Richmond, VA 23220.,Department of Psychology, Virginia Commonwealth University, Richmond, VA 23220
| | - Timothy P York
- Department of Human and Molecular Genetics, Virginia Commonwealth University, Richmond, VA 23220.,Department of Obstetrics and Gynecology, Virginia Commonwealth University, Richmond, VA 23220
| |
Collapse
|