1
|
Damarov IS, Korbolina EE, Rykova EY, Merkulova TI. Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response. Int J Mol Sci 2024; 25:9297. [PMID: 39273245 PMCID: PMC11394919 DOI: 10.3390/ijms25179297] [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/11/2024] [Revised: 08/14/2024] [Accepted: 08/26/2024] [Indexed: 09/15/2024] Open
Abstract
The goal of our study was to identify and assess the functionally significant SNPs with potentially important roles in the development of type 2 diabetes mellitus (T2DM) and/or their effect on individual response to antihyperglycemic medication with metformin. We applied a bioinformatics approach to identify the regulatory SNPs (rSNPs) associated with allele-asymmetric binding and expression events in our paired ChIP-seq and RNA-seq data for peripheral blood mononuclear cells (PBMCs) of nine healthy individuals. The rSNP outcomes were analyzed using public data from the GWAS (Genome-Wide Association Studies) and Genotype-Tissue Expression (GTEx). The differentially expressed genes (DEGs) between healthy and T2DM individuals (GSE221521), including metformin responders and non-responders (GSE153315), were searched for in GEO RNA-seq data. The DEGs harboring rSNPs were analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). We identified 14,796 rSNPs in the promoters of 5132 genes of human PBMCs. We found 4280 rSNPs to associate with both phenotypic traits (GWAS) and expression quantitative trait loci (eQTLs) from GTEx. Between T2DM patients and controls, 3810 rSNPs were detected in the promoters of 1284 DEGs. Based on the protein-protein interaction (PPI) network, we identified 31 upregulated hub genes, including the genes involved in inflammation, obesity, and insulin resistance. The top-ranked 10 enriched KEGG pathways for these hubs included insulin, AMPK, and FoxO signaling pathways. Between metformin responders and non-responders, 367 rSNPs were found in the promoters of 131 DEGs. Genes encoding transcription factors and transcription regulators were the most widely represented group and many were shown to be involved in the T2DM pathogenesis. We have formed a list of human rSNPs that add functional interpretation to the T2DM-association signals identified in GWAS. The results suggest candidate causal regulatory variants for T2DM, with strong enrichment in the pathways related to glucose metabolism, inflammation, and the effects of metformin.
Collapse
Affiliation(s)
- Igor S Damarov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| | - Elena E Korbolina
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| | - Elena Y Rykova
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
- Department of Engineering Problems of Ecology, Novosibirsk State Technical University, 630087 Novosibirsk, Russia
| | - Tatiana I Merkulova
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| |
Collapse
|
2
|
Zou X, Gomez ZW, Reddy TE, Allen AS, Majoros WH. Bayesian Estimation of Allele-Specific Expression in the Presence of Phasing Uncertainty. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.08.09.607371. [PMID: 39211106 PMCID: PMC11361064 DOI: 10.1101/2024.08.09.607371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Motivation Allele-specific expression (ASE) analyses aim to detect imbalanced expression of maternal versus paternal copies of an autosomal gene. Such allelic imbalance can result from a variety of cis-acting causes, including disruptive mutations within one copy of a gene that impact the stability of transcripts, as well as regulatory variants outside the gene that impact transcription initiation. Current methods for ASE estimation suffer from a number of shortcomings, such as relying on only one variant within a gene, assuming perfect phasing information across multiple variants within a gene, or failing to account for alignment biases and possible genotyping errors. Results We developed BEASTIE, a Bayesian hierarchical model designed for precise ASE quantification at the gene level, based on given genotypes and RNA-Seq data. BEASTIE addresses the complexities of allelic mapping bias, genotyping error, and phasing errors by incorporating empirical phasing error rates derived from Genome-in-a-Bottle individual NA12878. BEASTIE surpasses existing methods in accuracy, especially in scenarios with high phasing errors. This improvement is critical for identifying rare genetic variants often obscured by such errors. Through rigorous validation on simulated data and application to real data from the 1000 Genomes Project, we establish the robustness of BEASTIE. These findings underscore the value of BEASTIE in revealing patterns of ASE across gene sets and pathways. Availability and Implementation The software is freely available from https://github.com/x811zou/BEASTIE . BEASTIE is available as Python source code and as a Docker image. Supplementary information Additional information is available online.
Collapse
|
3
|
Qi G, Battle A. Computational methods for allele-specific expression in single cells. Trends Genet 2024:S0168-9525(24)00169-0. [PMID: 39127549 DOI: 10.1016/j.tig.2024.07.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2024] [Revised: 07/16/2024] [Accepted: 07/17/2024] [Indexed: 08/12/2024]
Abstract
Allele-specific expression (ASE) is a powerful signal that can be used to investigate multiple molecular mechanisms, such as cis-regulatory effects and imprinting. Single-cell RNA-sequencing (scRNA-seq) enables ASE characterization at the resolution of individual cells. In this review, we highlight the computational methods for processing and analyzing single-cell ASE data. We first describe a bioinformatics pipeline to obtain ASE counts from raw reads synthesized from previous literature. We then discuss statistical methods for detecting allelic imbalance and its variability across conditions using scRNA-seq data. In addition, we describe other methods that use single-cell ASE to address specific biological questions. Finally, we discuss future directions and emphasize the need for an integrated, optimized bioinformatics pipeline, and further development of statistical methods for different technologies.
Collapse
Affiliation(s)
- Guanghao Qi
- Department of Biostatistics, University of Washington, Seattle, WA 98195, USA.
| | - Alexis Battle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA; Department of Computer Science, Johns Hopkins University, Baltimore, MD 21218, USA; Department of Genetic Medicine, Johns Hopkins University, Baltimore, MD 21205, USA.
| |
Collapse
|
4
|
Tsouris A, Brach G, Schacherer J, Hou J. Non-additive genetic components contribute significantly to population-wide gene expression variation. CELL GENOMICS 2024; 4:100459. [PMID: 38190102 PMCID: PMC10794783 DOI: 10.1016/j.xgen.2023.100459] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 09/19/2023] [Accepted: 11/09/2023] [Indexed: 01/09/2024]
Abstract
Gene expression variation, an essential step between genotype and phenotype, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence gene expression variation remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originating from genetically divergent Saccharomyces cerevisiae isolates. Our analysis across 5,087 transcript abundance traits showed that non-additive components account for 36% of the gene expression variance on average. By comparing allele-specific read counts in parent-hybrid trios, we found that trans-regulatory changes underlie the majority of gene expression variation in the population. Remarkably, most cis-regulatory variations are also exaggerated or attenuated by additional trans effects. Overall, we showed that the transcriptome is globally buffered at the genetic level mainly due to trans-regulatory variation in the population.
Collapse
Affiliation(s)
- Andreas Tsouris
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France
| | - Gauthier Brach
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France; Institut Universitaire de France (IUF), Paris, France.
| | - Jing Hou
- Université de Strasbourg, CNRS, GMGM UMR, 7156 Strasbourg, France.
| |
Collapse
|
5
|
Beiki H, Murdoch BM, Park CA, Kern C, Kontechy D, Becker G, Rincon G, Jiang H, Zhou H, Thorne J, Koltes JE, Michal JJ, Davenport K, Rijnkels M, Ross PJ, Hu R, Corum S, McKay S, Smith TPL, Liu W, Ma W, Zhang X, Xu X, Han X, Jiang Z, Hu ZL, Reecy JM. Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology. Gigascience 2024; 13:giae019. [PMID: 38626724 PMCID: PMC11020238 DOI: 10.1093/gigascience/giae019] [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/11/2023] [Revised: 07/29/2023] [Accepted: 03/27/2024] [Indexed: 04/18/2024] Open
Abstract
BACKGROUND The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS These validated results show significant improvement over current bovine genome annotations.
Collapse
Affiliation(s)
- Hamid Beiki
- Department of Animal Science, Iowa State University, Ames, IA 50011, USA
| | - Brenda M Murdoch
- Department of Animal and Veterinary and Food Science, University of Idaho, ID 83844, USA
| | - Carissa A Park
- Department of Animal Science, Iowa State University, Ames, IA 50011, USA
| | - Chandlar Kern
- Department of Animal Science, Pennsylvania State University, PA 16802, USA
| | - Denise Kontechy
- Department of Animal and Veterinary and Food Science, University of Idaho, ID 83844, USA
| | - Gabrielle Becker
- Department of Animal and Veterinary and Food Science, University of Idaho, ID 83844, USA
| | | | - Honglin Jiang
- Department of Animal and Poultry Sciences, Virginia Tech, VA 24060, USA
| | - Huaijun Zhou
- Department of Animal Science, University of California, Davis, CA 95616, USA
| | - Jacob Thorne
- Department of Animal and Veterinary and Food Science, University of Idaho, ID 83844, USA
| | - James E Koltes
- Department of Animal Science, Iowa State University, Ames, IA 50011, USA
| | - Jennifer J Michal
- Department of Animal Science, Washington State University, WA 99164, USA
| | - Kimberly Davenport
- Department of Animal and Veterinary and Food Science, University of Idaho, ID 83844, USA
| | - Monique Rijnkels
- Department of Veterinary Integrative Biosciences, Texas A&M University, TX 77843, USA
| | - Pablo J Ross
- Department of Animal Science, University of California, Davis, CA 95616, USA
| | - Rui Hu
- Department of Animal and Poultry Sciences, Virginia Tech, VA 24060, USA
| | - Sarah Corum
- Zoetis, Parsippany-Troy Hills, NJ 07054, USA
| | | | | | - Wansheng Liu
- Department of Animal Science, Pennsylvania State University, PA 16802, USA
| | - Wenzhi Ma
- Department of Animal Science, Pennsylvania State University, PA 16802, USA
| | - Xiaohui Zhang
- Department of Animal Science, Washington State University, WA 99164, USA
| | - Xiaoqing Xu
- Department of Animal Science, University of California, Davis, CA 95616, USA
| | - Xuelei Han
- Department of Animal Science, Washington State University, WA 99164, USA
| | - Zhihua Jiang
- Department of Animal Science, Washington State University, WA 99164, USA
| | - Zhi-Liang Hu
- Department of Animal Science, Iowa State University, Ames, IA 50011, USA
| | - James M Reecy
- Department of Animal Science, Iowa State University, Ames, IA 50011, USA
| |
Collapse
|
6
|
Dyer NA, Lucas ER, Nagi SC, McDermott DP, Brenas JH, Miles A, Clarkson CS, Mawejje HD, Wilding CS, Halfon MS, Asma H, Heinz E, Donnelly MJ. Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.22.568226. [PMID: 38045426 PMCID: PMC10690255 DOI: 10.1101/2023.11.22.568226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2023]
Abstract
Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis-regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.
Collapse
Affiliation(s)
- Naomi A Dyer
- Department of Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
| | - Eric R Lucas
- Department of Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
| | - Sanjay C Nagi
- Department of Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
| | - Daniel P McDermott
- Department of Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
| | - Jon H Brenas
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Alistair Miles
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Chris S Clarkson
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Henry D Mawejje
- Infectious Diseases Research Collaboration (IDRC), Plot 2C Nakasero Hill Road, P.O.Box 7475, Kampala, Uganda
| | - Craig S Wilding
- School of Biological and Environmental Sciences, Liverpool John Moores University, Byrom Street, Liverpool, L3 3AF, UK
| | - Marc S Halfon
- Department of Biochemistry, Jacobs School of Medicine & Biomedical Sciences, University at Buffalo-State University of New York, 955 Main Street, Buffalo, New York 14203, USA
| | - Hasiba Asma
- Department of Biochemistry, Jacobs School of Medicine & Biomedical Sciences, University at Buffalo-State University of New York, 955 Main Street, Buffalo, New York 14203, USA
| | - Eva Heinz
- Department of Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
- Department of Clinical Sciences, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
| | - Martin J Donnelly
- Department of Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, L3 5QA, UK
| |
Collapse
|
7
|
Qi G, Strober BJ, Popp JM, Keener R, Ji H, Battle A. Single-cell allele-specific expression analysis reveals dynamic and cell-type-specific regulatory effects. Nat Commun 2023; 14:6317. [PMID: 37813843 PMCID: PMC10562474 DOI: 10.1038/s41467-023-42016-9] [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: 01/06/2023] [Accepted: 09/27/2023] [Indexed: 10/11/2023] Open
Abstract
Differential allele-specific expression (ASE) is a powerful tool to study context-specific cis-regulation of gene expression. Such effects can reflect the interaction between genetic or epigenetic factors and a measured context or condition. Single-cell RNA sequencing (scRNA-seq) allows the measurement of ASE at individual-cell resolution, but there is a lack of statistical methods to analyze such data. We present Differential Allelic Expression using Single-Cell data (DAESC), a powerful method for differential ASE analysis using scRNA-seq from multiple individuals, with statistical behavior confirmed through simulation. DAESC accounts for non-independence between cells from the same individual and incorporates implicit haplotype phasing. Application to data from 105 induced pluripotent stem cell (iPSC) lines identifies 657 genes dynamically regulated during endoderm differentiation, with enrichment for changes in chromatin state. Application to a type-2 diabetes dataset identifies several differentially regulated genes between patients and controls in pancreatic endocrine cells. DAESC is a powerful method for single-cell ASE analysis and can uncover novel insights on gene regulation.
Collapse
Affiliation(s)
- Guanghao Qi
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, USA
- Department of Biostatistics, University of Washington, Seattle, WA, 98195, USA
| | - Benjamin J Strober
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, 02115, USA
| | - Joshua M Popp
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, USA
| | - Rebecca Keener
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, USA
| | - Hongkai Ji
- Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, 21205, USA
| | - Alexis Battle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, USA.
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, 21218, USA.
- Department of Genetic Medicine, Johns Hopkins University, Baltimore, MD, 21205, USA.
| |
Collapse
|
8
|
Diaz-Ortiz ME, Jain N, Gallagher MD, Posavi M, Unger TL, Chen-Plotkin AS. Testing for Allele-specific Expression from Human Brain Samples. Bio Protoc 2023; 13:e4832. [PMID: 37817908 PMCID: PMC10560631 DOI: 10.21769/bioprotoc.4832] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Revised: 07/16/2023] [Accepted: 07/17/2023] [Indexed: 10/12/2023] Open
Abstract
Many single nucleotide polymorphisms (SNPs) identified by genome-wide association studies exert their effects on disease risk as expression quantitative trait loci (eQTL) via allele-specific expression (ASE). While databases for probing eQTLs in tissues from normal individuals exist, one may wish to ascertain eQTLs or ASE in specific tissues or disease-states not characterized in these databases. Here, we present a protocol to assess ASE of two possible target genes (GPNMB and KLHL7) of a known genome-wide association study (GWAS) Parkinson's disease (PD) risk locus in postmortem human brain tissue from PD and neurologically normal individuals. This was done using a sequence of RNA isolation, cDNA library generation, enrichment for transcripts of interest using customizable cDNA capture probes, paired-end RNA sequencing, and subsequent analysis. This method provides increased sensitivity relative to traditional bulk RNAseq-based and a blueprint that can be extended to the study of other genes, tissues, and disease states. Key features • Analysis of GPNMB allele-specific expression (ASE) in brain lysates from cognitively normal controls (NC) and Parkinson's disease (PD) individuals. • Builds on the ASE protocol of Mayba et al. (2014) and extends application from cells to human tissue. • Increased sensitivity by enrichment for desired transcript via RNA CaptureSeq (Mercer et al., 2014). • Optimized for human brain lysates from cingulate gyrus, caudate nucleus, and cerebellum.
Collapse
Affiliation(s)
- Maria E Diaz-Ortiz
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Nimansha Jain
- Department of Neurology, Hope Center for Neurological Disorders, Knight Alzheimer Disease, Research Center, Washington University, St. Louis, MO, USA
| | | | - Marijan Posavi
- Department of Biology, University of New Mexico, Albuquerque, NM, USA
| | - Travis L Unger
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Alice S Chen-Plotkin
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| |
Collapse
|
9
|
Durham M, Vadde S, Brooks AN. MET exon 14 skipping is overexpressed in an allele-specific manner in lung adenocarcinoma primary samples. MICROPUBLICATION BIOLOGY 2023; 2023:10.17912/micropub.biology.000957. [PMID: 37829573 PMCID: PMC10565573 DOI: 10.17912/micropub.biology.000957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Figures] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Revised: 09/20/2023] [Accepted: 09/22/2023] [Indexed: 10/14/2023]
Abstract
MET exon 14 skipping ( METΔ14 ) is a well-characterized oncogene in the Ras-MAPK pathway driving lung adenocarcinoma (LUAD). Previous studies on METΔ14 revealed this aberrantly spliced oncogene is expressed in LUAD primary samples and is associated with heterozygous somatic mutations and deletions near exon 14 splice sites. Upon further examination of DNA and RNA sequencing data from primary samples, we highlight that METΔ14 is overexpressed in an allele-specific manner. These data suggest that dose-dependence of METΔ14 may be critical to oncogenesis.
Collapse
Affiliation(s)
- Megan Durham
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Cruz, Santa Cruz, California, United States
| | - Swetha Vadde
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Cruz, Santa Cruz, California, United States
| | - Angela N Brooks
- Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, California, United States
| |
Collapse
|
10
|
Tsouris A, Brach G, Schacherer J, Hou J. Non-additive genetic components contribute significantly to population-wide gene expression variation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.21.550013. [PMID: 37546809 PMCID: PMC10401925 DOI: 10.1101/2023.07.21.550013] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Gene expression variation, an essential step between genomic variation and phenotypic landscape, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence the heritability of expression traits remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originated from genetically divergent yeast isolates. We estimated the broad- and narrow-sense heritability across 5,087 transcript abundance traits and showed that non-additive components account for 36% of the phenotypic variance on average. By comparing allelic expression ratios in the hybrid and the corresponding parental pair, we identified regulatory changes in 25% of all cases, with a majority acting in trans. We further showed that trans-regulation could underlie coordinated expression variation across highly connected genes, resulting in significantly higher non-additive variance and most likely in some of the missing heritability of gene expression traits.
Collapse
Affiliation(s)
- Andreas Tsouris
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Gauthier Brach
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
- Institut Universitaire de France (IUF), Paris, France
| | - Jing Hou
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| |
Collapse
|
11
|
Murani E, Hadlich F. Exploration of genotype-by-environment interactions affecting gene expression responses in porcine immune cells. Front Genet 2023; 14:1157267. [PMID: 37007953 PMCID: PMC10061014 DOI: 10.3389/fgene.2023.1157267] [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: 02/02/2023] [Accepted: 03/06/2023] [Indexed: 03/18/2023] Open
Abstract
As one of the keys to healthy performance, robustness of farm animals is gaining importance, and with this comes increasing interest in genetic dissection of genotype-by-environment interactions (G×E). Changes in gene expression are among the most sensitive responses conveying adaptation to environmental stimuli. Environmentally responsive regulatory variation thus likely plays a central role in G×E. In the present study, we set out to detect action of environmentally responsive cis-regulatory variation by the analysis of condition-dependent allele specific expression (cd-ASE) in porcine immune cells. For this, we harnessed mRNA-sequencing data of peripheral blood mononuclear cells (PBMCs) stimulated in vitro with lipopolysaccharide, dexamethasone, or their combination. These treatments mimic common challenges such as bacterial infection or stress, and induce vast transcriptome changes. About two thirds of the examined loci showed significant ASE in at least one treatment, and out of those about ten percent exhibited cd-ASE. Most of the ASE variants were not yet reported in the PigGTEx Atlas. Genes showing cd-ASE were enriched in cytokine signaling in immune system and include several key candidates for animal health. In contrast, genes showing no ASE featured cell-cycle related functions. We confirmed LPS-dependent ASE for one of the top candidates, SOD2, which ranks among the major response genes in LPS-stimulated monocytes. The results of the present study demonstrate the potential of in vitro cell models coupled with cd-ASE analysis for the investigation of G×E in farm animals. The identified loci may benefit efforts to unravel the genetic basis of robustness and improvement of health and welfare in pigs.
Collapse
Affiliation(s)
- Eduard Murani
- Institute of Genome Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | | |
Collapse
|
12
|
Diaz-Ortiz ME, Seo Y, Posavi M, Carceles Cordon M, Clark E, Jain N, Charan R, Gallagher MD, Unger TL, Amari N, Skrinak RT, Davila-Rivera R, Brody EM, Han N, Zack R, Van Deerlin VM, Tropea TF, Luk KC, Lee EB, Weintraub D, Chen-Plotkin AS. GPNMB confers risk for Parkinson's disease through interaction with α-synuclein. Science 2022; 377:eabk0637. [PMID: 35981040 PMCID: PMC9870036 DOI: 10.1126/science.abk0637] [Citation(s) in RCA: 59] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Many risk loci for Parkinson's disease (PD) have been identified by genome-wide association studies (GWASs), but target genes and mechanisms remain largely unknown. We linked the GWAS-derived chromosome 7 locus (sentinel single-nucleotide polymorphism rs199347) to GPNMB through colocalization analyses of expression quantitative trait locus and PD risk signals, confirmed by allele-specific expression studies in the human brain. In cells, glycoprotein nonmetastatic melanoma protein B (GPNMB) coimmunoprecipitated and colocalized with α-synuclein (aSyn). In induced pluripotent stem cell-derived neurons, loss of GPNMB resulted in loss of ability to internalize aSyn fibrils and develop aSyn pathology. In 731 PD and 59 control biosamples, GPNMB was elevated in PD plasma, associating with disease severity. Thus, GPNMB represents a PD risk gene with potential for biomarker development and therapeutic targeting.
Collapse
Affiliation(s)
- Maria E. Diaz-Ortiz
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Department of Bioengineering, School of Engineering and Applied Sciences, University of Pennsylvania, Philadelphia, PA, USA
| | - Yunji Seo
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Marijan Posavi
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Marc Carceles Cordon
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Elisia Clark
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Nimansha Jain
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Department of Neurology, Hope Center for Neurological Disorders, Knight Alzheimer Disease, Research Center, Washington University, St. Louis, MO, USA
| | - Rakshita Charan
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Flagship Pioneering, Cambridge, MA, USA
| | - Michael D. Gallagher
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Whitehead Institute for Biomedical Research, Cambridge, MA, USA
| | - Travis L. Unger
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Noor Amari
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - R. Tyler Skrinak
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Roseanne Davila-Rivera
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Eliza M. Brody
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Noah Han
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Rebecca Zack
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Vivianna M. Van Deerlin
- Center for Neurodegenerative Disease Research, Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Thomas F. Tropea
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Kelvin C. Luk
- Center for Neurodegenerative Disease Research, Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Edward B. Lee
- Center for Neurodegenerative Disease Research, Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Daniel Weintraub
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Parkinson’s Disease Research, Education and Clinical Center (PADRECC), Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA, USA
| | - Alice S. Chen-Plotkin
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| |
Collapse
|
13
|
Deng W, Mou T, Pawitan Y, Vu TN. Quantification of mutant–allele expression at isoform level in cancer from RNA-seq data. NAR Genom Bioinform 2022; 4:lqac052. [PMID: 35855322 PMCID: PMC9278039 DOI: 10.1093/nargab/lqac052] [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] [Received: 05/10/2021] [Revised: 06/26/2022] [Accepted: 07/04/2022] [Indexed: 11/13/2022] Open
Abstract
Even though the role of DNA mutations in cancer is well recognized, current quantification of the RNA expression, performed either at gene or isoform level, typically ignores the mutation status. Standard methods for estimating allele-specific expression (ASE) consider gene-level expression, but the functional impact of a mutation is best assessed at isoform level. Hence our goal is to quantify the mutant–allele expression at isoform level. We have developed and implemented a method, named MAX, for quantifying mutant–allele expression given a list of mutations. For a gene of interest, a mutant reference is constructed by incorporating all possible mutant versions of the wild-type isoforms in the transcriptome annotation. The mutant reference is then used for the RNA-seq reads mapping, which in principle works similarly for any quantification tool. We apply an alternating EM algorithm to the read-count data from the mapping step. In a simulation study, MAX performs well against standard isoform-quantification methods. Also, MAX achieves higher accuracy than conventional gene-based ASE methods such as ASEP. An analysis of a real dataset of acute myeloid leukemia reveals a subgroup of NPM1-mutated patients responding well to a kinase inhibitor. Our findings indicate that quantification of mutant–allele expression at isoform level is feasible and has potential added values for assessing the functional impact of DNA mutations in cancers.
Collapse
Affiliation(s)
- Wenjiang Deng
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet , Stockholm, Sweden
| | - Tian Mou
- School of Biomedical Engineering, Shenzhen University , Shenzhen, China
| | - Yudi Pawitan
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet , Stockholm, Sweden
| | - Trung Nghia Vu
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet , Stockholm, Sweden
| |
Collapse
|
14
|
Sen A, Huo Y, Elster J, Zage PE, McVicker G. Allele-specific expression reveals genes with recurrent cis-regulatory alterations in high-risk neuroblastoma. Genome Biol 2022; 23:71. [PMID: 35246212 PMCID: PMC8896304 DOI: 10.1186/s13059-022-02640-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Accepted: 02/23/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Neuroblastoma is a pediatric malignancy with a high frequency of metastatic disease at initial diagnosis. Neuroblastoma tumors have few recurrent protein-coding mutations but contain extensive somatic copy number alterations (SCNAs) suggesting that mutations that alter gene dosage are important drivers of tumorigenesis. Here, we analyze allele-specific expression in 96 high-risk neuroblastoma tumors to discover genes impacted by cis-acting mutations that alter dosage. RESULTS We identify 1043 genes with recurrent, neuroblastoma-specific allele-specific expression. While most of these genes lie within common SCNA regions, many of them exhibit allele-specific expression in copy neutral samples and these samples are enriched for mutations that are predicted to cause nonsense-mediated decay. Thus, both SCNA and non-SCNA mutations frequently alter gene expression in neuroblastoma. We focus on genes with neuroblastoma-specific allele-specific expression in the absence of SCNAs and find 26 such genes that have reduced expression in stage 4 disease. At least two of these genes have evidence for tumor suppressor activity including the transcription factor TFAP2B and the protein tyrosine phosphatase PTPRH. CONCLUSIONS In summary, our allele-specific expression analysis discovers genes that are recurrently dysregulated by both large SCNAs and other cis-acting mutations in high-risk neuroblastoma.
Collapse
Affiliation(s)
- Arko Sen
- Integrative Biology Laboratory, Salk Institute for Biological Studies, La Jolla, California, USA
| | - Yuchen Huo
- Department of Pediatrics, Division of Hematology-Oncology, University of California San Diego, La Jolla, California, USA
| | - Jennifer Elster
- Department of Pediatrics, Division of Hematology-Oncology, University of California San Diego, La Jolla, California, USA.,Peckham Center for Cancer and Blood Disorders, Rady Children's Hospital-San Diego, San Diego, California, USA
| | - Peter E Zage
- Department of Pediatrics, Division of Hematology-Oncology, University of California San Diego, La Jolla, California, USA.,Peckham Center for Cancer and Blood Disorders, Rady Children's Hospital-San Diego, San Diego, California, USA
| | - Graham McVicker
- Integrative Biology Laboratory, Salk Institute for Biological Studies, La Jolla, California, USA.
| |
Collapse
|
15
|
cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior. Cell Rep 2022; 38:110360. [PMID: 35172153 DOI: 10.1016/j.celrep.2022.110360] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Revised: 11/24/2021] [Accepted: 01/19/2022] [Indexed: 12/13/2022] Open
Abstract
How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F1 hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.
Collapse
|
16
|
Valentini S, Marchioretti C, Bisio A, Rossi A, Zaccara S, Romanel A, Inga A. TranSNPs: A class of functional SNPs affecting mRNA translation potential revealed by fraction-based allelic imbalance. iScience 2021; 24:103531. [PMID: 34917903 PMCID: PMC8666669 DOI: 10.1016/j.isci.2021.103531] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 10/27/2021] [Accepted: 11/23/2021] [Indexed: 12/23/2022] Open
Abstract
Few studies have explored the association between SNPs and alterations in mRNA translation potential. We developed an approach to identify SNPs that can mark allele-specific protein expression levels and could represent sources of inter-individual variation in disease risk. Using MCF7 cells under different treatments, we performed polysomal profiling followed by RNA sequencing of total or polysome-associated mRNA fractions and designed a computational approach to identify SNPs showing a significant change in the allelic balance between total and polysomal mRNA fractions. We identified 147 SNPs, 39 of which located in UTRs. Allele-specific differences at the translation level were confirmed in transfected MCF7 cells by reporter assays. Exploiting breast cancer data from TCGA we identified UTR SNPs demonstrating distinct prognosis features and altering binding sites of RNA-binding proteins. Our approach produced a catalog of tranSNPs, a class of functional SNPs associated with allele-specific translation and potentially endowed with prognostic value for disease risk.
Collapse
Affiliation(s)
- Samuel Valentini
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
| | - Caterina Marchioretti
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
- Department of Biomedical Sciences (DBS), University of Padova, 35131 Padova, Italy
| | - Alessandra Bisio
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
| | - Annalisa Rossi
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
| | - Sara Zaccara
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
- Weill Medical College, Cornell University, New York 10065, NY, USA
| | - Alessandro Romanel
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
| | - Alberto Inga
- Department of Cellular, Computational and Integrative Biology (CIBIO), University of Trento, 38123 Trento, Italy
| |
Collapse
|
17
|
Sen A, Prager BC, Zhong C, Park D, Zhu Z, Gimple RC, Wu Q, Bernatchez JA, Beck S, Clark AE, Siqueira-Neto JL, Rich JN, McVicker G. Leveraging Allele-Specific Expression for Therapeutic Response Gene Discovery in Glioblastoma. Cancer Res 2021; 82:377-390. [PMID: 34903607 DOI: 10.1158/0008-5472.can-21-0810] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Revised: 09/13/2021] [Accepted: 12/02/2021] [Indexed: 11/16/2022]
Abstract
Glioblastoma is the most prevalent primary malignant brain tumor in adults and is characterized by poor prognosis and universal tumor recurrence. Effective glioblastoma treatments are lacking, in part due to somatic mutations and epigenetic reprogramming that alter gene expression and confer drug resistance. To investigate recurrently dysregulated genes in glioblastoma we interrogated allele-specific expression (ASE), the difference in expression between two alleles of a gene, in glioblastoma stem cells (GSC) derived from 43 patients. A total of 118 genes were found with recurrent ASE preferentially in GSCs compared to normal tissues. These genes were enriched for apoptotic regulators, including schlafen family member 11 (SLFN11). Loss of SLFN11 gene expression was associated with aberrant promoter methylation and conferred resistance to chemotherapy and PARP inhibition. Conversely, low SLFN11 expression rendered GSCs susceptible to the oncolytic flavivirus Zika. This discovery effort based upon ASE revealed novel points of vulnerability in GSCs, suggesting a potential alternative treatment strategy for chemotherapy resistant glioblastoma.
Collapse
Affiliation(s)
- Arko Sen
- Salk Institute for Biological Studies
| | - Briana C Prager
- Cleveland Clinic Lerner College of Medicine, Cleveland Clinic
| | | | | | - Zhe Zhu
- Medicine, University of California, San Diego
| | | | - Qiulian Wu
- Medicine, University of California - San Diego School of Medicine
| | - Jean A Bernatchez
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego
| | | | | | | | - Jeremy N Rich
- Department of Neurology, University of Pittsburgh Cancer Institute
| | | |
Collapse
|
18
|
Sherbina K, León-Novelo LG, Nuzhdin SV, McIntyre LM, Marroni F. Power calculator for detecting allelic imbalance using hierarchical Bayesian model. BMC Res Notes 2021; 14:436. [PMID: 34838135 PMCID: PMC8626927 DOI: 10.1186/s13104-021-05851-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Accepted: 11/15/2021] [Indexed: 11/10/2022] Open
Abstract
OBJECTIVE Allelic imbalance (AI) is the differential expression of the two alleles in a diploid. AI can vary between tissues, treatments, and environments. Methods for testing AI exist, but methods are needed to estimate type I error and power for detecting AI and difference of AI between conditions. As the costs of the technology plummet, what is more important: reads or replicates? RESULTS We find that a minimum of 2400, 480, and 240 allele specific reads divided equally among 12, 5, and 3 replicates is needed to detect a 10, 20, and 30%, respectively, deviation from allelic balance in a condition with power > 80%. A minimum of 960 and 240 allele specific reads divided equally among 8 replicates is needed to detect a 20 or 30% difference in AI between conditions with comparable power. Higher numbers of replicates increase power more than adding coverage without affecting type I error. We provide a Python package that enables simulation of AI scenarios and enables individuals to estimate type I error and power in detecting AI and differences in AI between conditions.
Collapse
Affiliation(s)
- Katrina Sherbina
- Quantitative and Computational Biology Section, University of Southern California, Los Angeles, CA, 90046, USA
| | - Luis G León-Novelo
- Department of Biostatistics and Data Science, The University of Texas Health Science Center at Houston-School of Public Health, Houston, TX, 77030, USA
| | - Sergey V Nuzhdin
- Molecular and Computational Biology Section, University of Southern California, Los Angeles, CA, 90046, USA
| | - Lauren M McIntyre
- Genetics Institute and Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL, 32603, USA
| | - Fabio Marroni
- Dipartimento di Scienze Agroalimentari, Ambientali e Animali, Università di Udine, 33100, Udine, Italy.
| |
Collapse
|
19
|
Alonso L, Piron A, Morán I, Guindo-Martínez M, Bonàs-Guarch S, Atla G, Miguel-Escalada I, Royo R, Puiggròs M, Garcia-Hurtado X, Suleiman M, Marselli L, Esguerra JLS, Turatsinze JV, Torres JM, Nylander V, Chen J, Eliasson L, Defrance M, Amela R, Mulder H, Gloyn AL, Groop L, Marchetti P, Eizirik DL, Ferrer J, Mercader JM, Cnop M, Torrents D. TIGER: The gene expression regulatory variation landscape of human pancreatic islets. Cell Rep 2021; 37:109807. [PMID: 34644572 PMCID: PMC8864863 DOI: 10.1016/j.celrep.2021.109807] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Revised: 07/23/2021] [Accepted: 09/16/2021] [Indexed: 12/30/2022] Open
Abstract
Genome-wide association studies (GWASs) identified hundreds of signals associated with type 2 diabetes (T2D). To gain insight into their underlying molecular mechanisms, we have created the translational human pancreatic islet genotype tissue-expression resource (TIGER), aggregating >500 human islet genomic datasets from five cohorts in the Horizon 2020 consortium T2DSystems. We impute genotypes using four reference panels and meta-analyze cohorts to improve the coverage of expression quantitative trait loci (eQTL) and develop a method to combine allele-specific expression across samples (cASE). We identify >1 million islet eQTLs, 53 of which colocalize with T2D signals. Among them, a low-frequency allele that reduces T2D risk by half increases CCND2 expression. We identify eight cASE colocalizations, among which we found a T2D-associated SLC30A8 variant. We make all data available through the TIGER portal (http://tiger.bsc.es), which represents a comprehensive human islet genomic data resource to elucidate how genetic variation affects islet function and translates into therapeutic insight and precision medicine for T2D.
Collapse
Affiliation(s)
- Lorena Alonso
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | - Anthony Piron
- ULB Center for Diabetes Research, Université Libre de Bruxelles, Brussels 1070, Belgium; Interuniversity Institute of Bioinformatics in Brussels (IB2), Brussels 1050, Belgium
| | - Ignasi Morán
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | - Marta Guindo-Martínez
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | - Sílvia Bonàs-Guarch
- Bioinformatics and Genomics Program, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona 08003, Spain; Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM) Barcelona 08013, Spain
| | - Goutham Atla
- Bioinformatics and Genomics Program, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona 08003, Spain; Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM) Barcelona 08013, Spain
| | - Irene Miguel-Escalada
- Bioinformatics and Genomics Program, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona 08003, Spain; Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM) Barcelona 08013, Spain
| | - Romina Royo
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | - Montserrat Puiggròs
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | - Xavier Garcia-Hurtado
- Bioinformatics and Genomics Program, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona 08003, Spain; Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM) Barcelona 08013, Spain
| | - Mara Suleiman
- Department of Clinical and Experimental Medicine and AOUP Cisanello University Hospital, University of Pisa, Pisa 56126, Italy
| | - Lorella Marselli
- Department of Clinical and Experimental Medicine and AOUP Cisanello University Hospital, University of Pisa, Pisa 56126, Italy
| | - Jonathan L S Esguerra
- Unit of Islet Cell Exocytosis, Lund University Diabetes Centre, Malmö 214 28, Sweden
| | | | - Jason M Torres
- Clinical Trial Service Unit and Epidemiological Studies Unit, Nuffield Department of Population Health, University of Oxford, Oxford OX3 7LF, UK; Wellcome Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford OX3 7LF, UK
| | - Vibe Nylander
- Oxford Centre for Diabetes, Endocrinology, and Metabolism, Radcliffe Department of Medicine, University of Oxford, Oxford OX3 7LE, UK
| | - Ji Chen
- Exeter Centre of Excellence for Diabetes Research (EXCEED), University of Exeter Medical School, Exeter EX4 4PY, UK
| | - Lena Eliasson
- Unit of Islet Cell Exocytosis, Lund University Diabetes Centre, Malmö 214 28, Sweden
| | - Matthieu Defrance
- ULB Center for Diabetes Research, Université Libre de Bruxelles, Brussels 1070, Belgium
| | - Ramon Amela
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | - Hindrik Mulder
- Unit of Molecular Metabolism, Lund University Diabetes Centre, Malmö 214 28, Sweden
| | - Anna L Gloyn
- Wellcome Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford OX3 7LF, UK; Oxford Centre for Diabetes, Endocrinology, and Metabolism, Radcliffe Department of Medicine, University of Oxford, Oxford OX3 7LE, UK; Division of Endocrinology, Department of Pediatrics, Stanford University School of Medicine, Stanford, CA 94304, USA; NIHR Oxford Biomedical Research Centre, Churchill Hospital, Oxford OX3 7DQ, UK; Stanford Diabetes Research Centre, Stanford University, Stanford, CA 94305, USA
| | - Leif Groop
- Unit of Islet Cell Exocytosis, Lund University Diabetes Centre, Malmö 214 28, Sweden; Unit of Molecular Metabolism, Lund University Diabetes Centre, Malmö 214 28, Sweden; Finnish Institute of Molecular Medicine Finland (FIMM), Helsinki University, Helsinki 00014, Finland
| | - Piero Marchetti
- Department of Clinical and Experimental Medicine and AOUP Cisanello University Hospital, University of Pisa, Pisa 56126, Italy
| | - Decio L Eizirik
- ULB Center for Diabetes Research, Université Libre de Bruxelles, Brussels 1070, Belgium; WELBIO, Université Libre de Bruxelles, Brussels 1050, Belgium
| | - Jorge Ferrer
- Bioinformatics and Genomics Program, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona 08003, Spain; Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM) Barcelona 08013, Spain; Section of Epigenomics and Disease, Department of Medicine, Imperial College London, London SW7 2AZ, UK
| | - Josep M Mercader
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain; Programs in Metabolism and Medical and Population Genetics, Broad Institute of Harvard and MIT, Cambridge, MA 02142, USA; Diabetes Unit and Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Medicine, Harvard Medical School, Boston, MA 02115, USA.
| | - Miriam Cnop
- ULB Center for Diabetes Research, Université Libre de Bruxelles, Brussels 1070, Belgium; Division of Endocrinology, Erasmus Hospital, Université Libre de Bruxelles, Brussels 1070, Belgium.
| | - David Torrents
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain; Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona 08010, Spain.
| |
Collapse
|
20
|
Abstract
Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.
Collapse
|
21
|
Berardi AE, Esfeld K, Jäggi L, Mandel T, Cannarozzi GM, Kuhlemeier C. Complex evolution of novel red floral color in Petunia. THE PLANT CELL 2021; 33:2273-2295. [PMID: 33871652 PMCID: PMC8364234 DOI: 10.1093/plcell/koab114] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2021] [Accepted: 04/12/2021] [Indexed: 05/20/2023]
Abstract
Red flower color has arisen multiple times and is generally associated with hummingbird pollination. The majority of evolutionary transitions to red color proceeded from purple lineages and tend to be genetically simple, almost always involving a few loss-of-function mutations of major phenotypic effect. Here we report on the complex evolution of a novel red floral color in the hummingbird-pollinated Petunia exserta (Solanaceae) from a colorless ancestor. The presence of a red color is remarkable because the genus cannot synthesize red anthocyanins and P. exserta retains a nonfunctional copy of the key MYB transcription factor AN2. We show that moderate upregulation and a shift in tissue specificity of an AN2 paralog, DEEP PURPLE, restores anthocyanin biosynthesis in P. exserta. An essential shift in anthocyanin hydroxylation occurred through rebalancing the expression of three hydroxylating genes. Furthermore, the downregulation of an acyltransferase promotes reddish hues in typically purple pigments by preventing acyl group decoration of anthocyanins. This study presents a rare case of a genetically complex evolutionary transition toward the gain of a novel red color.
Collapse
Affiliation(s)
- Andrea E. Berardi
- Institute of Plant Sciences, University of Bern, Bern 3013, Switzerland
| | - Korinna Esfeld
- Institute of Plant Sciences, University of Bern, Bern 3013, Switzerland
| | - Lea Jäggi
- Institute of Plant Sciences, University of Bern, Bern 3013, Switzerland
| | - Therese Mandel
- Institute of Plant Sciences, University of Bern, Bern 3013, Switzerland
| | | | - Cris Kuhlemeier
- Institute of Plant Sciences, University of Bern, Bern 3013, Switzerland
- Author for correspondence:
| |
Collapse
|
22
|
Thind AS, Monga I, Thakur PK, Kumari P, Dindhoria K, Krzak M, Ranson M, Ashford B. Demystifying emerging bulk RNA-Seq applications: the application and utility of bioinformatic methodology. Brief Bioinform 2021; 22:6330938. [PMID: 34329375 DOI: 10.1093/bib/bbab259] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2021] [Revised: 06/14/2021] [Accepted: 06/18/2021] [Indexed: 12/13/2022] Open
Abstract
Significant innovations in next-generation sequencing techniques and bioinformatics tools have impacted our appreciation and understanding of RNA. Practical RNA sequencing (RNA-Seq) applications have evolved in conjunction with sequence technology and bioinformatic tools advances. In most projects, bulk RNA-Seq data is used to measure gene expression patterns, isoform expression, alternative splicing and single-nucleotide polymorphisms. However, RNA-Seq holds far more hidden biological information including details of copy number alteration, microbial contamination, transposable elements, cell type (deconvolution) and the presence of neoantigens. Recent novel and advanced bioinformatic algorithms developed the capacity to retrieve this information from bulk RNA-Seq data, thus broadening its scope. The focus of this review is to comprehend the emerging bulk RNA-Seq-based analyses, emphasizing less familiar and underused applications. In doing so, we highlight the power of bulk RNA-Seq in providing biological insights.
Collapse
Affiliation(s)
- Amarinder Singh Thind
- University of Wollongong, Wollongong, Australia.,Illawarra Health and Medical Research Institute, Wollongong, Australia
| | - Isha Monga
- Columbia University, New York City, NY, USA
| | | | - Pallawi Kumari
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh, India
| | - Kiran Dindhoria
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh, India
| | | | - Marie Ranson
- University of Wollongong, Wollongong, Australia.,Illawarra Health and Medical Research Institute, Wollongong, Australia
| | - Bruce Ashford
- University of Wollongong, Wollongong, Australia.,Illawarra Health and Medical Research Institute, Wollongong, Australia
| |
Collapse
|
23
|
Abstract
Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.
Collapse
Affiliation(s)
- Siobhan Cleary
- School of Mathematics, Statistics and Applied Mathematics, National University of Ireland, Galway H91 H3CY, Ireland;
| | - Cathal Seoighe
- School of Mathematics, Statistics and Applied Mathematics, National University of Ireland, Galway H91 H3CY, Ireland;
| |
Collapse
|
24
|
Oh S, Shao J, Mitra J, Xiong F, D'Antonio M, Wang R, Garcia-Bassets I, Ma Q, Zhu X, Lee JH, Nair SJ, Yang F, Ohgi K, Frazer KA, Zhang ZD, Li W, Rosenfeld MG. Enhancer release and retargeting activates disease-susceptibility genes. Nature 2021; 595:735-740. [PMID: 34040254 PMCID: PMC11171441 DOI: 10.1038/s41586-021-03577-1] [Citation(s) in RCA: 67] [Impact Index Per Article: 22.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Accepted: 04/23/2021] [Indexed: 02/04/2023]
Abstract
The functional engagement between an enhancer and its target promoter ensures precise gene transcription1. Understanding the basis of promoter choice by enhancers has important implications for health and disease. Here we report that functional loss of a preferred promoter can release its partner enhancer to loop to and activate an alternative promoter (or alternative promoters) in the neighbourhood. We refer to this target-switching process as 'enhancer release and retargeting'. Genetic deletion, motif perturbation or mutation, and dCas9-mediated CTCF tethering reveal that promoter choice by an enhancer can be determined by the binding of CTCF at promoters, in a cohesin-dependent manner-consistent with a model of 'enhancer scanning' inside the contact domain. Promoter-associated CTCF shows a lower affinity than that at chromatin domain boundaries and often lacks a preferred motif orientation or a partnering CTCF at the cognate enhancer, suggesting properties distinct from boundary CTCF. Analyses of cancer mutations, data from the GTEx project and risk loci from genome-wide association studies, together with a focused CRISPR interference screen, reveal that enhancer release and retargeting represents an overlooked mechanism that underlies the activation of disease-susceptibility genes, as exemplified by a risk locus for Parkinson's disease (NUCKS1-RAB7L1) and three loci associated with cancer (CLPTM1L-TERT, ZCCHC7-PAX5 and PVT1-MYC).
Collapse
Affiliation(s)
- Soohwan Oh
- Howard Hughes Medical Institute, Department and School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Jiaofang Shao
- Department of Biochemistry and Molecular Biology, McGovern Medical School, University of Texas Health Science Center, Houston, TX, USA
| | - Joydeep Mitra
- Department of Genetics, Albert Einstein College of Medicine, New York, NY, USA
| | - Feng Xiong
- Department of Biochemistry and Molecular Biology, McGovern Medical School, University of Texas Health Science Center, Houston, TX, USA
| | - Matteo D'Antonio
- Department of Pediatrics, School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Ruoyu Wang
- Department of Biochemistry and Molecular Biology, McGovern Medical School, University of Texas Health Science Center, Houston, TX, USA
- Graduate School of Biomedical Sciences, University of Texas MD Anderson Cancer Center and UTHealth, Houston, TX, USA
| | - Ivan Garcia-Bassets
- Department of Medicine, School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Qi Ma
- Howard Hughes Medical Institute, Department and School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Xiaoyu Zhu
- Department of Biochemistry and Molecular Biology, McGovern Medical School, University of Texas Health Science Center, Houston, TX, USA
| | - Joo-Hyung Lee
- Department of Biochemistry and Molecular Biology, McGovern Medical School, University of Texas Health Science Center, Houston, TX, USA
| | - Sreejith J Nair
- Howard Hughes Medical Institute, Department and School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Feng Yang
- Howard Hughes Medical Institute, Department and School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Kenneth Ohgi
- Howard Hughes Medical Institute, Department and School of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Kelly A Frazer
- Department of Pediatrics, School of Medicine, University of California San Diego, La Jolla, CA, USA
- Institute for Genomic Medicine, University of California San Diego, La Jolla, CA, USA
| | - Zhengdong D Zhang
- Department of Genetics, Albert Einstein College of Medicine, New York, NY, USA
| | - Wenbo Li
- Department of Biochemistry and Molecular Biology, McGovern Medical School, University of Texas Health Science Center, Houston, TX, USA.
- Graduate School of Biomedical Sciences, University of Texas MD Anderson Cancer Center and UTHealth, Houston, TX, USA.
| | - Michael G Rosenfeld
- Howard Hughes Medical Institute, Department and School of Medicine, University of California San Diego, La Jolla, CA, USA.
| |
Collapse
|
25
|
Mendelevich A, Vinogradova S, Gupta S, Mironov AA, Sunyaev SR, Gimelbrant AA. Replicate sequencing libraries are important for quantification of allelic imbalance. Nat Commun 2021; 12:3370. [PMID: 34099647 PMCID: PMC8184992 DOI: 10.1038/s41467-021-23544-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Accepted: 04/30/2021] [Indexed: 12/13/2022] Open
Abstract
A sensitive approach to quantitative analysis of transcriptional regulation in diploid organisms is analysis of allelic imbalance (AI) in RNA sequencing (RNA-seq) data. A near-universal practice in such studies is to prepare and sequence only one library per RNA sample. We present theoretical and experimental evidence that data from a single RNA-seq library is insufficient for reliable quantification of the contribution of technical noise to the observed AI signal; consequently, reliance on one-replicate experimental design can lead to unaccounted-for variation in error rates in allele-specific analysis. We develop a computational approach, Qllelic, that accurately accounts for technical noise by making use of replicate RNA-seq libraries. Testing on new and existing datasets shows that application of Qllelic greatly decreases false positive rate in allele-specific analysis while conserving appropriate signal, and thus greatly improves reproducibility of AI estimates. We explore sources of technical overdispersion in observed AI signal and conclude by discussing design of RNA-seq studies addressing two biologically important questions: quantification of transcriptome-wide AI in one sample, and differential analysis of allele-specific expression between samples.
Collapse
Affiliation(s)
- Asia Mendelevich
- Skolkovo Institute of Science and Technology, Moscow, Russia.
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA.
| | - Svetlana Vinogradova
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA
| | - Saumya Gupta
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA
- Broad Institute of Harvard and MIT, Cambridge, USA
| | - Andrey A Mironov
- Lomonosov Moscow State University, Faculty of Bioengineering and Bioinformatics, Moscow, Russia
- Institute of Information Transmission Problems, Russian Academy of Sciences, Moscow, Russia
| | - Shamil R Sunyaev
- Department of Biomedical Informatics, Harvard Medical School, Boston, USA
- Division of Genetics, Brigham and Women's Hospital, Boston, USA
| | - Alexander A Gimelbrant
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA.
- Broad Institute of Harvard and MIT, Cambridge, USA.
| |
Collapse
|
26
|
A survey of cancer genome signatures identifies genes connected to distinct chromosomal instability phenotypes. THE PHARMACOGENOMICS JOURNAL 2021; 21:390-401. [PMID: 33731882 DOI: 10.1038/s41397-021-00217-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 01/16/2021] [Accepted: 01/27/2021] [Indexed: 01/31/2023]
Abstract
Certain breast and ovarian cancers are characterised by high levels of chromosomal instability. We established a suite of eleven SNP array-based signatures of various forms of chromosomal instability. To understand what biological mechanisms might underpin these signatures, we generated and assembled genetic-feature data including allele-specific expression, fusion genes, gene expression, methylation, somatic coding mutations and protein expression. For each signature, we extracted a compendium of significantly associated genetic features using machine learning. We established an association between subchromosomal allelic imbalance-based measures and DNA repair genes. Numerical chromosomal instability and chromothripsis were found to have distinct genetic associations but shared a relationship to mitotic processes, while whole-genome doubling was characterised by TP53 mutation, and high AURKA and GINS1 expression. Furthermore, we identified two genetically distinct forms of tandem duplicator phenotypes. These findings identify potentially novel genomic targets for validation and drug development for specific subsets of breast and ovarian cancer.
Collapse
|
27
|
Atak ZK, Taskiran II, Demeulemeester J, Flerin C, Mauduit D, Minnoye L, Hulselmans G, Christiaens V, Ghanem GE, Wouters J, Aerts S. Interpretation of allele-specific chromatin accessibility using cell state-aware deep learning. Genome Res 2021; 31:1082-1096. [PMID: 33832990 PMCID: PMC8168584 DOI: 10.1101/gr.260851.120] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Accepted: 04/05/2021] [Indexed: 12/26/2022]
Abstract
Genomic sequence variation within enhancers and promoters can have a significant impact on the cellular state and phenotype. However, sifting through the millions of candidate variants in a personal genome or a cancer genome, to identify those that impact cis-regulatory function, remains a major challenge. Interpretation of noncoding genome variation benefits from explainable artificial intelligence to predict and interpret the impact of a mutation on gene regulation. Here we generate phased whole genomes with matched chromatin accessibility, histone modifications, and gene expression for 10 melanoma cell lines. We find that training a specialized deep learning model, called DeepMEL2, on melanoma chromatin accessibility data can capture the various regulatory programs of the melanocytic and mesenchymal-like melanoma cell states. This model outperforms motif-based variant scoring, as well as more generic deep learning models. We detect hundreds to thousands of allele-specific chromatin accessibility variants (ASCAVs) in each melanoma genome, of which 15%-20% can be explained by gains or losses of transcription factor binding sites. A considerable fraction of ASCAVs are caused by changes in AP-1 binding, as confirmed by matched ChIP-seq data to identify allele-specific binding of JUN and FOSL1. Finally, by augmenting the DeepMEL2 model with ChIP-seq data for GABPA, the TERT promoter mutation, as well as additional ETS motif gains, can be identified with high confidence. In conclusion, we present a new integrative genomics approach and a deep learning model to identify and interpret functional enhancer mutations with allelic imbalance of chromatin accessibility and gene expression.
Collapse
Affiliation(s)
- Zeynep Kalender Atak
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Ibrahim Ihsan Taskiran
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Jonas Demeulemeester
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium.,Cancer Genomics Laboratory, The Francis Crick Institute, London NW1 1AT, United Kingdom
| | - Christopher Flerin
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - David Mauduit
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Liesbeth Minnoye
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Gert Hulselmans
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Valerie Christiaens
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Ghanem-Elias Ghanem
- Institut Jules Bordet, Université Libre de Bruxelles, 1000 Brussels, Belgium
| | - Jasper Wouters
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| | - Stein Aerts
- VIB-KU Leuven Center for Brain and Disease Research, 3000 Leuven, Belgium.,KU Leuven, Department of Human Genetics KU Leuven, 3000 Leuven, Belgium
| |
Collapse
|
28
|
Cheng Z, Vermeulen M, Rollins-Green M, DeVeale B, Babak T. Cis-regulatory mutations with driver hallmarks in major cancers. iScience 2021; 24:102144. [PMID: 33665563 PMCID: PMC7903341 DOI: 10.1016/j.isci.2021.102144] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2020] [Revised: 09/02/2020] [Accepted: 01/25/2021] [Indexed: 12/05/2022] Open
Abstract
Despite the recent availability of complete genome sequences of tumors from thousands of patients, isolating disease-causing (driver) non-coding mutations from the plethora of somatic variants remains challenging, and only a handful of validated examples exist. By integrating whole-genome sequencing, genetic data, and allele-specific gene expression from TCGA, we identified 320 somatic non-coding mutations that affect gene expression in cis (FDR<0.25). These mutations cluster into 47 cis-regulatory elements that modulate expression of their subject genes through diverse molecular mechanisms. We further show that these mutations have hallmark features of non-coding drivers; namely, that they preferentially disrupt transcription factor binding motifs, are associated with a selective advantage, increased oncogene expression and decreased tumor suppressor expression. Enrichment of functional non-coding somatic mutations predicts drivers Elevated variant allele frequencies are consistent with roles in tumorigenesis Putative non-coding drivers disrupt transcription factor binding motifs Predicted drivers associate with increased oncogene and decreased TSG expression
Collapse
Affiliation(s)
- Zhongshan Cheng
- Department of Biology, Queen's University, Kingston, ON K7L 3N6, Canada
| | - Michael Vermeulen
- Department of Biology, Queen's University, Kingston, ON K7L 3N6, Canada
| | | | - Brian DeVeale
- The Eli and Edythe Broad Center of Regeneration Medicine and Stem Cell Research, Center for Reproductive Sciences, University of California, San Francisco, San Francisco, CA 94143, USA
| | - Tomas Babak
- Department of Biology, Queen's University, Kingston, ON K7L 3N6, Canada
| |
Collapse
|
29
|
Fan J, Wang X, Xiao R, Li M. Detecting cell-type-specific allelic expression imbalance by integrative analysis of bulk and single-cell RNA sequencing data. PLoS Genet 2021; 17:e1009080. [PMID: 33661921 PMCID: PMC7963069 DOI: 10.1371/journal.pgen.1009080] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2020] [Revised: 03/16/2021] [Accepted: 02/09/2021] [Indexed: 12/27/2022] Open
Abstract
Allelic expression imbalance (AEI), quantified by the relative expression of two alleles of a gene in a diploid organism, can help explain phenotypic variations among individuals. Traditional methods detect AEI using bulk RNA sequencing (RNA-seq) data, a data type that averages out cell-to-cell heterogeneity in gene expression across cell types. Since the patterns of AEI may vary across different cell types, it is desirable to study AEI in a cell-type-specific manner. Although this can be achieved by single-cell RNA sequencing (scRNA-seq), it requires full-length transcript to be sequenced in single cells of a large number of individuals, which are still cost prohibitive to generate. To overcome this limitation and utilize the vast amount of existing disease relevant bulk tissue RNA-seq data, we developed BSCET, which enables the characterization of cell-type-specific AEI in bulk RNA-seq data by integrating cell type composition information inferred from a small set of scRNA-seq samples, possibly obtained from an external dataset. By modeling covariate effect, BSCET can also detect genes whose cell-type-specific AEI are associated with clinical factors. Through extensive benchmark evaluations, we show that BSCET correctly detected genes with cell-type-specific AEI and differential AEI between healthy and diseased samples using bulk RNA-seq data. BSCET also uncovered cell-type-specific AEIs that were missed in bulk data analysis when the directions of AEI are opposite in different cell types. We further applied BSCET to two pancreatic islet bulk RNA-seq datasets, and detected genes showing cell-type-specific AEI that are related to the progression of type 2 diabetes. Since bulk RNA-seq data are easily accessible, BSCET provides a convenient tool to integrate information from scRNA-seq data to gain insight on AEI with cell type resolution. Results from such analysis will advance our understanding of cell type contributions in human diseases. Detection of allelic expression imbalance (AEI), a phenomenon where the two alleles of a gene differ in their expression magnitude, is a key step towards the understanding of phenotypic variations among individuals. Existing methods detect AEI using bulk RNA sequencing (RNA-seq) data and ignore AEI variations among different cell types. Although single-cell RNA sequencing (scRNA-seq) has enabled the characterization of cell-to-cell heterogeneity in gene expression, the high costs have limited its application in AEI analysis. To overcome this limitation, we developed BSCET to characterize cell-type-specific AEI using the widely available bulk RNA-seq data by integrating cell-type composition information inferred from scRNA-seq samples. Since the degree of AEI may vary with disease phenotypes, we further extended BSCET to detect genes whose cell-type-specific AEIs are associated with clinical factors. Through extensive benchmark evaluations and analyses of two pancreatic islet bulk RNA-seq datasets, we demonstrated BSCET’s ability to refine bulk-level AEI to cell-type resolution, and to identify genes whose cell-type-specific AEIs are associated with the progression of type 2 diabetes. With the vast amount of easily accessible bulk RNA-seq data, we believe BSCET will be a valuable tool for elucidating cell type contributions in human diseases.
Collapse
Affiliation(s)
- Jiaxin Fan
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
| | - Xuran Wang
- Department of Statistics and Data Science, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of America
| | - Rui Xiao
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
- * E-mail: (RX); (ML)
| | - Mingyao Li
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
- * E-mail: (RX); (ML)
| |
Collapse
|
30
|
Tomlinson MJ, Polson SW, Qiu J, Lake JA, Lee W, Abasht B. Investigation of allele specific expression in various tissues of broiler chickens using the detection tool VADT. Sci Rep 2021; 11:3968. [PMID: 33597613 PMCID: PMC7889858 DOI: 10.1038/s41598-021-83459-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2020] [Accepted: 02/01/2021] [Indexed: 12/30/2022] Open
Abstract
Differential abundance of allelic transcripts in a diploid organism, commonly referred to as allele specific expression (ASE), is a biologically significant phenomenon and can be examined using single nucleotide polymorphisms (SNPs) from RNA-seq. Quantifying ASE aids in our ability to identify and understand cis-regulatory mechanisms that influence gene expression, and thereby assist in identifying causal mutations. This study examines ASE in breast muscle, abdominal fat, and liver of commercial broiler chickens using variants called from a large sub-set of the samples (n = 68). ASE analysis was performed using a custom software called VCF ASE Detection Tool (VADT), which detects ASE of biallelic SNPs using a binomial test. On average ~ 174,000 SNPs in each tissue passed our filtering criteria and were considered informative, of which ~ 24,000 (~ 14%) showed ASE. Of all ASE SNPs, only 3.7% exhibited ASE in all three tissues, with ~ 83% showing ASE specific to a single tissue. When ASE genes (genes containing ASE SNPs) were compared between tissues, the overlap among all three tissues increased to 20.1%. Our results indicate that ASE genes show tissue-specific enrichment patterns, but all three tissues showed enrichment for pathways involved in translation.
Collapse
Affiliation(s)
- M Joseph Tomlinson
- Department of Animal and Food Sciences, University of Delaware, 531 South College Ave, Newark, DE, 19716, USA.,Center for Bioinformatics and Computational Biology, University of Delaware, Newark, USA
| | - Shawn W Polson
- Department of Computer and Information Sciences, University of Delaware, Newark, USA.,Department of Biological Sciences, University of Delaware, Newark, USA.,Center for Bioinformatics and Computational Biology, University of Delaware, Newark, USA
| | - Jing Qiu
- Department of Applied Economics and Statistics, University of Delaware, Newark, USA.,Center for Bioinformatics and Computational Biology, University of Delaware, Newark, USA
| | - Juniper A Lake
- Department of Animal and Food Sciences, University of Delaware, 531 South College Ave, Newark, DE, 19716, USA.,Center for Bioinformatics and Computational Biology, University of Delaware, Newark, USA
| | - William Lee
- Maple Leaf Farms, Inc., Leesburg, IN, 46538, USA
| | - Behnam Abasht
- Department of Animal and Food Sciences, University of Delaware, 531 South College Ave, Newark, DE, 19716, USA. .,Center for Bioinformatics and Computational Biology, University of Delaware, Newark, USA.
| |
Collapse
|
31
|
McGirr JA, Martin CH. Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence. Mol Biol Evol 2021; 38:405-423. [PMID: 32877534 PMCID: PMC7826174 DOI: 10.1093/molbev/msaa218] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.
Collapse
Affiliation(s)
- Joseph A McGirr
- Environmental Toxicology Department, University of California, Davis, CA
| | - Christopher H Martin
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, CA
| |
Collapse
|
32
|
Robles-Espinoza CD, Mohammadi P, Bonilla X, Gutierrez-Arcelus M. Allele-specific expression: applications in cancer and technical considerations. Curr Opin Genet Dev 2021; 66:10-19. [PMID: 33383480 PMCID: PMC7985293 DOI: 10.1016/j.gde.2020.10.007] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Revised: 10/26/2020] [Accepted: 10/31/2020] [Indexed: 11/18/2022]
Abstract
Allele-specific gene expression can influence disease traits. Non-coding germline genetic variants that alter regulatory elements can cause allele-specific gene expression and contribute to cancer susceptibility. In tumors, both somatic copy number alterations and somatic single nucleotide variants have been shown to lead to allele-specific expression of genes, many of which are considered drivers of tumor growth. Here, we review recent studies revealing the pervasive presence of this phenomenon in cancer susceptibility and progression. Furthermore, we underscore the importance of careful experimental design and computational analysis for accurate allelic expression quantification and avoidance of false positives. Finally, we discuss additional methodological challenges encountered in cancer studies and in the burgeoning field of single-cell transcriptomics.
Collapse
Affiliation(s)
- Carla Daniela Robles-Espinoza
- Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Campus Juriquilla, Boulevard Juriquilla 3001, Santiago de Querétaro 76230, Mexico; Wellcome Sanger Institute, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Pejman Mohammadi
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, USA; Scripps Translational Science Institute, The Scripps Research Institute, La Jolla, CA, USA
| | - Ximena Bonilla
- Department of Computer Science, ETH Zurich, Universitätsstr. 6, 8092 Zürich, Switzerland; Swiss Institute of Bioinformatics, Quartier Sorge - Bâtiment Amphipôle, Lausanne 1015, Switzerland; University Hospital Zurich, Rämistrasse 100, 8091 Zürich, Switzerland
| | - Maria Gutierrez-Arcelus
- Center for Data Sciences, Brigham and Women's Hospital, Harvard Medical School, Boston, MA 02115, USA; Division of Genetics, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA 02115, USA; Division of Rheumatology, Inflammation and Immunity, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA 02115, USA; Program in Medical and Population Genetics, Broad Institute, Cambridge, MA 02142, USA; Division of Immunology, Department of Pediatrics, Boston Children's Hospital, Harvard Medical School, Boston, Massachusetts, USA.
| |
Collapse
|
33
|
Yarahmadov T, Robinson S, Hanemian M, Pulver V, Kuhlemeier C. Identification of transcription factors controlling floral morphology in wild Petunia species with contrasting pollination syndromes. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2020; 104:289-301. [PMID: 32780443 PMCID: PMC7693086 DOI: 10.1111/tpj.14962] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2019] [Accepted: 07/15/2020] [Indexed: 05/29/2023]
Abstract
Adaptation to different pollinators is an important driver of speciation in the angiosperms. Genetic approaches such as QTL mapping have been successfully used to identify the underlying speciation genes. However, these methods are often limited by widespread suppression of recombination due to divergence between species. While the mutations that caused the interspecific differences in floral color and scent have been elucidated in a variety of plant genera, the genes that are responsible for morphological differences remain mostly unknown. Differences in floral organ length determine the pollination efficiency of hawkmoths and hummingbirds, and therefore the genes that control these differences are potential speciation genes. Identifying such genes is challenging, especially in non-model species and when studying complex traits for which little prior genetic and biochemical knowledge is available. Here we combine transcriptomics with detailed growth analysis to identify candidate transcription factors underlying interspecific variation in the styles of Petunia flowers. Starting from a set of 2284 genes, stepwise filtering for expression in styles, differential expression between species, correlation with growth-related traits, allele-specific expression in interspecific hybrids, and/or high-impact polymorphisms resulted in a set of 43 candidate speciation genes. Validation by virus-induced gene silencing identified two MYB transcription factors, EOBI and EOBII, that were previously shown to regulate floral scent emission, a trait associated with pollination by hawkmoths.
Collapse
Affiliation(s)
- Tural Yarahmadov
- Institute of Plant SciencesUniversity of BernAltenbergrain 21BernCH‐3013Switzerland
- Department of BioMedical ResearchUniversity of BernBernCH‐3008Switzerland
| | - Sarah Robinson
- Institute of Plant SciencesUniversity of BernAltenbergrain 21BernCH‐3013Switzerland
- Sainsbury LaboratoryUniversity of CambridgeCambridgeCB2 1LRUK
| | - Mathieu Hanemian
- Institute of Plant SciencesUniversity of BernAltenbergrain 21BernCH‐3013Switzerland
- LIPMUniversité de ToulouseINRAECNRSCastanet‐TolosanFrance
| | - Valentin Pulver
- Institute of Plant SciencesUniversity of BernAltenbergrain 21BernCH‐3013Switzerland
| | - Cris Kuhlemeier
- Institute of Plant SciencesUniversity of BernAltenbergrain 21BernCH‐3013Switzerland
| |
Collapse
|
34
|
Dong L, Wang J, Wang G. BYASE: a Python library for estimating gene and isoform level allele-specific expression. Bioinformatics 2020; 36:4955-4956. [DOI: 10.1093/bioinformatics/btaa636] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2020] [Revised: 05/12/2020] [Accepted: 07/10/2020] [Indexed: 12/15/2022] Open
Abstract
Abstract
Summary
Allele-specific expression (ASE) is involved in many important biological mechanisms. We present a python package BYASE and its graphical user interface (GUI) tool BYASE-GUI for the identification of ASE from single-end and paired-end RNA-seq data based on Bayesian inference, which can simultaneously report differences in gene-level and isoform-level expression. BYASE uses both phased SNPs and non-phased SNPs, and supports polyploid organisms.
Availability and implementation
The source codes of BYASE and BYASE-GUI are freely available at https://github.com/ncjllld/byase and https://github.com/ncjllld/byase_gui.
Supplementary information
Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Lili Dong
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China
| | - Jianan Wang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China
| | - Guohua Wang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China
- State Key Laboratory of Tree Genetics and Breeding, Northeast Forestry University, Harbin, Heilongjiang 150001, China
| |
Collapse
|
35
|
Liu Y, Li C, Shen S, Chen X, Szlachta K, Edmonson MN, Shao Y, Ma X, Hyle J, Wright S, Ju B, Rusch MC, Liu Y, Li B, Macias M, Tian L, Easton J, Qian M, Yang JJ, Hu S, Look AT, Zhang J. Discovery of regulatory noncoding variants in individual cancer genomes by using cis-X. Nat Genet 2020; 52:811-818. [PMID: 32632335 PMCID: PMC7679232 DOI: 10.1038/s41588-020-0659-5] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Accepted: 06/05/2020] [Indexed: 12/30/2022]
Abstract
We developed cis-X, a computational method for discovering regulatory noncoding variants in cancer by integrating whole-genome and transcriptome sequencing data from a single cancer sample. cis-X first finds aberrantly cis-activated genes that exhibit allele-specific expression accompanied by an elevated outlier expression. It then searches for causal noncoding variants that may introduce aberrant transcription factor binding motifs or enhancer hijacking by structural variations. Analysis of 13 T-lineage acute lymphoblastic leukemias identified a recurrent intronic variant predicted to cis-activate the TAL1 oncogene, a finding validated in vivo by chromatin immunoprecipitation sequencing of a patient-derived xenograft. Candidate oncogenes include the prolactin receptor PRLR activated by a focal deletion that removes a CTCF-insulated neighborhood boundary. cis-X may be applied to pediatric and adult solid tumors that are aneuploid and heterogeneous. In contrast to existing approaches, which require large sample cohorts, cis-X enables the discovery of regulatory noncoding variants in individual cancer genomes.
Collapse
Affiliation(s)
- Yu Liu
- Pediatric Translational Medicine Institute, Shanghai Children's Medical Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China. .,Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA.
| | - Chunliang Li
- Department of Tumor Cell Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Shuhong Shen
- Pediatric Translational Medicine Institute, Shanghai Children's Medical Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China.,Key Laboratory of Pediatric Hematology & Oncology Ministry of Health, Department of Hematology & Oncology, Shanghai Children's Medical Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Xiaolong Chen
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Karol Szlachta
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Michael N Edmonson
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Ying Shao
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Xiaotu Ma
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Judith Hyle
- Department of Tumor Cell Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Shaela Wright
- Department of Tumor Cell Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Bensheng Ju
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Michael C Rusch
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Yanling Liu
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Benshang Li
- Pediatric Translational Medicine Institute, Shanghai Children's Medical Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China.,Key Laboratory of Pediatric Hematology & Oncology Ministry of Health, Department of Hematology & Oncology, Shanghai Children's Medical Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Michael Macias
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Liqing Tian
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - John Easton
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Maoxiang Qian
- Department of Pharmaceutical Sciences, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Jun J Yang
- Department of Pharmaceutical Sciences, St. Jude Children's Research Hospital, Memphis, TN, USA.,Hematological Malignancies Program, St. Jude Children's Research Hospital, Memphis, TN, USA.,Department of Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Shaoyan Hu
- Children's Hospital of Soochow University, Suzhou, China
| | - A Thomas Look
- Department of Pediatric Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA.,Division of Pediatric Hematology-Oncology, Boston Children's Hospital, Boston, MA, USA
| | - Jinghui Zhang
- Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, TN, USA.
| |
Collapse
|
36
|
McGirr JA, Martin CH. Ecological divergence in sympatry causes gene misexpression in hybrids. Mol Ecol 2020; 29:2707-2721. [PMID: 32557903 PMCID: PMC8209238 DOI: 10.1111/mec.15512] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 05/21/2020] [Accepted: 06/01/2020] [Indexed: 12/17/2022]
Abstract
Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.
Collapse
Affiliation(s)
- Joseph A. McGirr
- Department of Biology, University of North Carolina, Chapel
Hill, NC 27514
| | - Christopher H. Martin
- Department of Biology, University of North Carolina, Chapel
Hill, NC 27514
- Department of Integrative Biology and Museum of Vertebrate
Zoology, University of California, Berkeley, CA 94720
| |
Collapse
|
37
|
Fan J, Hu J, Xue C, Zhang H, Susztak K, Reilly MP, Xiao R, Li M. ASEP: Gene-based detection of allele-specific expression across individuals in a population by RNA sequencing. PLoS Genet 2020; 16:e1008786. [PMID: 32392242 PMCID: PMC7241832 DOI: 10.1371/journal.pgen.1008786] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Revised: 05/21/2020] [Accepted: 04/21/2020] [Indexed: 12/16/2022] Open
Abstract
Allele-specific expression (ASE) analysis, which quantifies the relative expression of two alleles in a diploid individual, is a powerful tool for identifying cis-regulated gene expression variations that underlie phenotypic differences among individuals. Existing methods for gene-level ASE detection analyze one individual at a time, therefore failing to account for shared information across individuals. Failure to accommodate such shared information not only reduces power, but also makes it difficult to interpret results across individuals. However, when only RNA sequencing (RNA-seq) data are available, ASE detection across individuals is challenging because the data often include individuals that are either heterozygous or homozygous for the unobserved cis-regulatory SNP, leading to sample heterogeneity as only those heterozygous individuals are informative for ASE, whereas those homozygous individuals have balanced expression. To simultaneously model multi-individual information and account for such heterogeneity, we developed ASEP, a mixture model with subject-specific random effect to account for multi-SNP correlations within the same gene. ASEP only requires RNA-seq data, and is able to detect gene-level ASE under one condition and differential ASE between two conditions (e.g., pre- versus post-treatment). Extensive simulations demonstrated the convincing performance of ASEP under a wide range of scenarios. We applied ASEP to a human kidney RNA-seq dataset, identified ASE genes and validated our results with two published eQTL studies. We further applied ASEP to a human macrophage RNA-seq dataset, identified genes showing evidence of differential ASE between M0 and M1 macrophages, and confirmed our findings by results from cardiometabolic trait-relevant genome-wide association studies. To the best of our knowledge, ASEP is the first method for gene-level ASE detection at the population level that only requires the use of RNA-seq data. With the growing adoption of RNA-seq, we believe ASEP will be well-suited for various ASE studies for human diseases.
Collapse
Affiliation(s)
- Jiaxin Fan
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
| | - Jian Hu
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
| | - Chenyi Xue
- Division of Cardiology, Department of Medicine, Columbia University Irving Medical Center, New York City, New York, United States of America
| | - Hanrui Zhang
- Division of Cardiology, Department of Medicine, Columbia University Irving Medical Center, New York City, New York, United States of America
| | - Katalin Susztak
- Departments of Medicine and Genetics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
| | - Muredach P. Reilly
- Division of Cardiology, Department of Medicine, Columbia University Irving Medical Center, New York City, New York, United States of America
- The Irving Institute for Clinical and Translational Research, Columbia University Irving Medical Center, New York City, New York, United States of America
| | - Rui Xiao
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
| | - Mingyao Li
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America
| |
Collapse
|
38
|
Przytycki PF, Singh M. Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations. Cell Syst 2020; 10:193-203.e4. [PMID: 32078798 PMCID: PMC7457951 DOI: 10.1016/j.cels.2020.01.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2019] [Revised: 12/04/2019] [Accepted: 01/22/2020] [Indexed: 01/23/2023]
Abstract
Identifying cancer-relevant mutations in noncoding regions is challenging due to the large numbers of such mutations, their low levels of recurrence, and difficulties in interpreting their functional impact. To uncover genes that are dysregulated due to somatic mutations in cis, we build upon the concept of differential allele-specific expression (ASE) and introduce methods to identify genes within an individual's cancer whose ASE differs from what is found in matched normal tissue. When applied to breast cancer tumor samples, our methods detect the known allele-specific effects of copy number variation and nonsense-mediated decay. Further, genes that are found to recurrently exhibit differential ASE across samples are cancer relevant. Genes with cis mutations are enriched for differential ASE, and we find 147 potentially functional noncoding mutations cis to genes that exhibit significant differential ASE. We conclude that differential ASE is a promising means for discovering gene dysregulation due to cis noncoding mutations.
Collapse
Affiliation(s)
- Pawel F Przytycki
- Department of Computer Science, Princeton University, Princeton, NJ 08544, USA; Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08544, USA
| | - Mona Singh
- Department of Computer Science, Princeton University, Princeton, NJ 08544, USA; Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08544, USA.
| |
Collapse
|
39
|
Clayton EA, Khalid S, Ban D, Wang L, Jordan IK, McDonald JF. Tumor suppressor genes and allele-specific expression: mechanisms and significance. Oncotarget 2020; 11:462-479. [PMID: 32064050 PMCID: PMC6996918 DOI: 10.18632/oncotarget.27468] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Accepted: 01/13/2020] [Indexed: 12/12/2022] Open
Abstract
Recent findings indicate that allele-specific expression (ASE) at specific cancer driver gene loci may be of importance in onset/progression of the disease. Of particular interest are loss-of-function (LOF) of tumor suppressor gene (TSGs) alleles. While LOF tumor suppressor mutations are typically considered to be recessive, if these mutant alleles can be significantly differentially expressed relative to wild-type alleles in heterozygotes, the clinical consequences could be significant. LOF TSG alleles are shown to be segregating at high frequencies in world-wide populations of normal/healthy individuals. Matched sets of normal and tumor tissues isolated from 233 cancer patients representing four diverse tumor types demonstrate functionally important changes in patterns of ASE in individuals heterozygous for LOF TSG alleles associated with cancer onset/progression. While a variety of molecular mechanisms were identified as potentially contributing to changes in ASE patterns in cancer, changes in DNA copy number and allele-specific alternative splicing possibly mediated by antisense RNA emerged as predominant factors. In conclusion, LOF TSGs are segregating in human populations at significant frequencies indicating that many otherwise healthy individuals are at elevated risk of developing cancer. Changes in ASE between normal and cancer tissues indicates that LOF TSG alleles may contribute to cancer onset/progression even when heterozygous with wild-type functional alleles.
Collapse
Affiliation(s)
- Evan A. Clayton
- Integrated Cancer Research Center, School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
| | - Shareef Khalid
- Integrated Cancer Research Center, School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
| | - Dongjo Ban
- Integrated Cancer Research Center, School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
| | - Lu Wang
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- PanAmerican Bioinformatics Institute, Cali, Colombia
| | - I. King Jordan
- Integrated Cancer Research Center, School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- PanAmerican Bioinformatics Institute, Cali, Colombia
- Applied Bioinformatics Laboratory, Atlanta, GA, USA
| | - John F. McDonald
- Integrated Cancer Research Center, School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
| |
Collapse
|
40
|
Banos G, Clark EL, Bush SJ, Dutta P, Bramis G, Arsenos G, Hume DA, Psifidi A. Genetic and genomic analyses underpin the feasibility of concomitant genetic improvement of milk yield and mastitis resistance in dairy sheep. PLoS One 2019; 14:e0214346. [PMID: 31765378 PMCID: PMC6876840 DOI: 10.1371/journal.pone.0214346] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Accepted: 10/31/2019] [Indexed: 11/19/2022] Open
Abstract
Milk yield is the most important dairy sheep trait and constitutes the key genetic improvement goal via selective breeding. Mastitis is one of the most prevalent diseases, significantly impacting on animal welfare, milk yield and quality, while incurring substantial costs. Our objectives were to determine the feasibility of a concomitant genetic improvement programme for enhanced milk production and resistance to mastitis. Individual records for milk yield, and four mastitis-related traits (milk somatic cell count, California Mastitis Test score, total viable bacterial count in milk and clinical mastitis presence) were collected monthly throughout lactation for 609 ewes of the Chios breed. All ewes were genotyped with a mastitis specific custom-made 960 single nucleotide polymorphism (SNP) array. We performed targeted genomic association studies, (co)variance component estimation and pathway enrichment analysis, and characterised gene expression levels and the extent of allelic expression imbalance. Presence of heritable variation for milk yield was confirmed. There was no significant genetic correlation between milk yield and mastitis traits. Environmental factors appeared to favour both milk production and udder health. There were no overlapping of SNPs associated with mastitis resistance and milk yield in Chios sheep. Furthermore, four distinct Quantitative Trait Loci (QTLs) affecting milk yield were detected on chromosomes 2, 12, 16 and 19, in locations other than those previously identified to affect mastitis resistance. Five genes (DNAJA1, GHR, LYPLA1, NUP35 and OXCT1) located within the QTL regions were highly expressed in both the mammary gland and milk transcriptome, suggesting involvement in milk synthesis and production. Furthermore, the expression of two of these genes (NUP35 and OXCT1) was enriched in immune tissues implying a potentially pleiotropic effect or likely role in milk production during udder infection, which needs to be further elucidated in future studies. In conclusion, the absence of genetic antagonism between milk yield and mastitis resistance suggests that simultaneous genetic improvement of both traits be achievable.
Collapse
Affiliation(s)
- Georgios Banos
- Scotland’s Rural College, Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
- School of Veterinary Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
| | - Emily L. Clark
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
| | - Stephen J. Bush
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
- Nuffield Department of Clinical Medicine, University of Oxford, John Radcliffe Hospital, Headington, Oxford, England, United Kingdom
| | - Prasun Dutta
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
| | - Georgios Bramis
- School of Veterinary Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
| | - Georgios Arsenos
- School of Veterinary Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
| | - David A. Hume
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
- Mater Research Institute-University of Queensland, Translational Research Institute, Woolloongabba, Australia
| | - Androniki Psifidi
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, Scotland, United Kingdom
- Royal Veterinary College, University of London, Hatfield, England, United Kingdom
- * E-mail: ,
| |
Collapse
|
41
|
D'Antonio M, Reyna J, Jakubosky D, Donovan MKR, Bonder MJ, Matsui H, Stegle O, Nariai N, D'Antonio-Chronowska A, Frazer KA. Systematic genetic analysis of the MHC region reveals mechanistic underpinnings of HLA type associations with disease. eLife 2019; 8:e48476. [PMID: 31746734 PMCID: PMC6904215 DOI: 10.7554/elife.48476] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Accepted: 11/19/2019] [Indexed: 02/06/2023] Open
Abstract
The MHC region is highly associated with autoimmune and infectious diseases. Here we conduct an in-depth interrogation of associations between genetic variation, gene expression and disease. We create a comprehensive map of regulatory variation in the MHC region using WGS from 419 individuals to call eight-digit HLA types and RNA-seq data from matched iPSCs. Building on this regulatory map, we explored GWAS signals for 4083 traits, detecting colocalization for 180 disease loci with eQTLs. We show that eQTL analyses taking HLA type haplotypes into account have substantially greater power compared with only using single variants. We examined the association between the 8.1 ancestral haplotype and delayed colonization in Cystic Fibrosis, postulating that downregulation of RNF5 expression is the likely causal mechanism. Our study provides insights into the genetic architecture of the MHC region and pinpoints disease associations that are due to differential expression of HLA genes and non-HLA genes.
Collapse
Affiliation(s)
- Matteo D'Antonio
- Institute for Genomic MedicineUniversity of California, San DiegoSan DiegoUnited States
- Department of PediatricsRady Children’s Hospital, University of California, San DiegoSan DiegoUnited States
| | - Joaquin Reyna
- Department of PediatricsRady Children’s Hospital, University of California, San DiegoSan DiegoUnited States
- Biomedical Sciences Graduate ProgramUniversity of California, San DiegoLa JollaUnited States
| | - David Jakubosky
- Biomedical Sciences Graduate ProgramUniversity of California, San DiegoLa JollaUnited States
- Bioinformatics and Systems Biology Graduate ProgramUniversity of California, San DiegoSan DiegoUnited States
| | - Margaret KR Donovan
- Bioinformatics and Systems Biology Graduate ProgramUniversity of California, San DiegoSan DiegoUnited States
- Department of Biomedical InformaticsUniversity of California, San DiegoSan DiegoUnited States
| | - Marc-Jan Bonder
- European Molecular Biology Laboratory, European Bioinformatics InstituteCambridgeUnited Kingdom
| | - Hiroko Matsui
- Institute for Genomic MedicineUniversity of California, San DiegoSan DiegoUnited States
| | - Oliver Stegle
- European Molecular Biology Laboratory, European Bioinformatics InstituteCambridgeUnited Kingdom
| | - Naoki Nariai
- Department of PediatricsRady Children’s Hospital, University of California, San DiegoSan DiegoUnited States
| | - Agnieszka D'Antonio-Chronowska
- Institute for Genomic MedicineUniversity of California, San DiegoSan DiegoUnited States
- Department of PediatricsRady Children’s Hospital, University of California, San DiegoSan DiegoUnited States
| | - Kelly A Frazer
- Institute for Genomic MedicineUniversity of California, San DiegoSan DiegoUnited States
- Department of PediatricsRady Children’s Hospital, University of California, San DiegoSan DiegoUnited States
| |
Collapse
|
42
|
Xie J, Ji T, Ferreira MAR, Li Y, Patel BN, Rivera RM. Modeling allele-specific expression at the gene and SNP levels simultaneously by a Bayesian logistic mixed regression model. BMC Bioinformatics 2019; 20:530. [PMID: 31660858 PMCID: PMC6819473 DOI: 10.1186/s12859-019-3141-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Accepted: 10/09/2019] [Indexed: 12/29/2022] Open
Abstract
BACKGROUND High-throughput sequencing experiments, which can determine allele origins, have been used to assess genome-wide allele-specific expression. Despite the amount of data generated from high-throughput experiments, statistical methods are often too simplistic to understand the complexity of gene expression. Specifically, existing methods do not test allele-specific expression (ASE) of a gene as a whole and variation in ASE within a gene across exons separately and simultaneously. RESULTS We propose a generalized linear mixed model to close these gaps, incorporating variations due to genes, single nucleotide polymorphisms (SNPs), and biological replicates. To improve reliability of statistical inferences, we assign priors on each effect in the model so that information is shared across genes in the entire genome. We utilize Bayesian model selection to test the hypothesis of ASE for each gene and variations across SNPs within a gene. We apply our method to four tissue types in a bovine study to de novo detect ASE genes in the bovine genome, and uncover intriguing predictions of regulatory ASEs across gene exons and across tissue types. We compared our method to competing approaches through simulation studies that mimicked the real datasets. The R package, BLMRM, that implements our proposed algorithm, is publicly available for download at https://github.com/JingXieMIZZOU/BLMRM . CONCLUSIONS We will show that the proposed method exhibits improved control of the false discovery rate and improved power over existing methods when SNP variation and biological variation are present. Besides, our method also maintains low computational requirements that allows for whole genome analysis.
Collapse
Affiliation(s)
- Jing Xie
- Department of Statistics, University of Missouri at Columbia, Columbia, 65211 MO USA
| | - Tieming Ji
- Department of Statistics, University of Missouri at Columbia, Columbia, 65211 MO USA
| | | | - Yahan Li
- Division of Animal Science, University of Missouri at Columbia, Columbia, 65211 MO USA
| | - Bhaumik N. Patel
- Division of Animal Science, University of Missouri at Columbia, Columbia, 65211 MO USA
| | - Rocio M. Rivera
- Division of Animal Science, University of Missouri at Columbia, Columbia, 65211 MO USA
| |
Collapse
|
43
|
Benaglio P, D'Antonio-Chronowska A, Ma W, Yang F, Young Greenwald WW, Donovan MKR, DeBoever C, Li H, Drees F, Singhal S, Matsui H, van Setten J, Sotoodehnia N, Gaulton KJ, Smith EN, D'Antonio M, Rosenfeld MG, Frazer KA. Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits. Nat Genet 2019; 51:1506-1517. [PMID: 31570892 PMCID: PMC6858543 DOI: 10.1038/s41588-019-0499-3] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2018] [Accepted: 08/15/2019] [Indexed: 12/15/2022]
Abstract
The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes.
Collapse
Affiliation(s)
- Paola Benaglio
- Department of Pediatrics, Rady Children's Hospital, Division of Genome Information Sciences, University of California, San Diego, La Jolla, CA, USA
| | | | - Wubin Ma
- Howard Hughes Medical Institute, Department of Medicine, University of California, San Diego, La Jolla, CA, USA
| | - Feng Yang
- Howard Hughes Medical Institute, Department of Medicine, University of California, San Diego, La Jolla, CA, USA
| | | | - Margaret K R Donovan
- Bioinformatics and Systems Biology, University of California, San Diego, La Jolla, CA, USA.,Department of Biomedical Informatics, University of California, San Diego, La Jolla, CA, USA
| | - Christopher DeBoever
- Bioinformatics and Systems Biology, University of California, San Diego, La Jolla, CA, USA
| | - He Li
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, USA
| | - Frauke Drees
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, USA
| | - Sanghamitra Singhal
- Department of Pediatrics, Rady Children's Hospital, Division of Genome Information Sciences, University of California, San Diego, La Jolla, CA, USA
| | - Hiroko Matsui
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, USA
| | - Jessica van Setten
- Department of Cardiology, University Medical Center Utrecht, University of Utrecht, Utrecht, the Netherlands
| | - Nona Sotoodehnia
- Department of Medicine, Cardiovascular Health Research Unit, Division of Cardiology, University of Washington, Seattle, WA, USA.,Department of Epidemiology, Cardiovascular Health Research Unit, Division of Cardiology, University of Washington, Seattle, WA, USA
| | - Kyle J Gaulton
- Department of Pediatrics, Rady Children's Hospital, Division of Genome Information Sciences, University of California, San Diego, La Jolla, CA, USA
| | - Erin N Smith
- Department of Pediatrics, Rady Children's Hospital, Division of Genome Information Sciences, University of California, San Diego, La Jolla, CA, USA
| | - Matteo D'Antonio
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, USA
| | - Michael G Rosenfeld
- Howard Hughes Medical Institute, Department of Medicine, University of California, San Diego, La Jolla, CA, USA.
| | - Kelly A Frazer
- Department of Pediatrics, Rady Children's Hospital, Division of Genome Information Sciences, University of California, San Diego, La Jolla, CA, USA. .,Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, USA.
| |
Collapse
|
44
|
Salavati M, Bush SJ, Palma-Vera S, McCulloch MEB, Hume DA, Clark EL. Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep. Front Genet 2019; 10:863. [PMID: 31608110 PMCID: PMC6761296 DOI: 10.3389/fgene.2019.00863] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 08/19/2019] [Indexed: 12/13/2022] Open
Abstract
Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.
Collapse
Affiliation(s)
- Mazdak Salavati
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| | - Stephen J. Bush
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| | - Sergio Palma-Vera
- Leibniz Institute for Farm Animal Biology (FBN), Institute for Reproductive Biology, Dummerstorf, Germany
| | - Mary E. B. McCulloch
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| | - David A. Hume
- Mater Research Institute-University of Queensland, Translational Research Institute, Woolloongabba, QLD, Australia
| | - Emily L. Clark
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| |
Collapse
|
45
|
McGirr JA, Martin CH. Hybrid gene misregulation in multiple developing tissues within a recent adaptive radiation of Cyprinodon pupfishes. PLoS One 2019; 14:e0218899. [PMID: 31291291 PMCID: PMC6619667 DOI: 10.1371/journal.pone.0218899] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 06/11/2019] [Indexed: 12/24/2022] Open
Abstract
Genetic incompatibilities constitute the final stages of reproductive isolation and speciation, but little is known about incompatibilities that occur within recent adaptive radiations among closely related diverging populations. Crossing divergent species to form hybrids can break up coadapted variation, resulting in genetic incompatibilities within developmental networks shaping divergent adaptive traits. We crossed two closely related sympatric Cyprinodon pupfish species–a dietary generalist and a specialized molluscivore–and measured expression levels in their F1 hybrids to identify regulatory variation underlying the novel craniofacial morphology found in this recent microendemic adaptive radiation. We extracted mRNA from eight day old whole-larvae tissue and from craniofacial tissues dissected from 17–20 day old larvae to compare gene expression between a total of seven F1 hybrids and 24 individuals from parental species populations. We found 3.9% of genes differentially expressed between generalists and molluscivores in whole-larvae tissues and 0.6% of genes differentially expressed in craniofacial tissue. We found that 2.1% of genes were misregulated in whole-larvae hybrids whereas 19.1% of genes were misregulated in hybrid craniofacial tissues, after correcting for sequencing biases. We also measured allele specific expression across 15,429 heterozygous sites to identify putative compensatory regulatory mechanisms underlying differential expression between generalists and molluscivores. Together, our results highlight the importance of considering misregulation as an early indicator of genetic incompatibilities in the context of rapidly diverging adaptive radiations and suggests that compensatory regulatory divergence drives hybrid gene misregulation in developing tissues that give rise to novel craniofacial traits.
Collapse
Affiliation(s)
- Joseph A. McGirr
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- * E-mail:
| | - Christopher H. Martin
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Department of Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America
| |
Collapse
|
46
|
Allelic Switching of DLX5, GRB10, and SVOPL during Colorectal Cancer Tumorigenesis. Int J Genomics 2019; 2019:1287671. [PMID: 31093489 PMCID: PMC6481143 DOI: 10.1155/2019/1287671] [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] [Received: 10/01/2018] [Revised: 01/28/2019] [Accepted: 02/07/2019] [Indexed: 11/26/2022] Open
Abstract
Allele-specific expression (ASE) is found in approximately 20-30% of human genes. During tumorigenesis, ASE changes due to somatic alterations that change the regulatory landscape. In colorectal cancer (CRC), many chromosomes show frequent gains or losses while homozygosity of chromosome 7 is rare. We hypothesized that genes essential to survival show allele-specific expression (ASE) on both alleles of chromosome 7. Using a panel of 21 recently established low-passage CRC cell lines, we performed ASE analysis by hybridizing DNA and cDNA to Infinium HumanExome-12 v1 BeadChips containing cSNPs in 392 chromosome 7 genes. The results of this initial analysis were extended and validated in a set of 89 paired normal mucosa and CRC samples. We found that 14% of genes showed ASE in one or more cell lines and identified allelic switching of the potential cell survival genes DLX5, GRB10, and SVOPL on chromosome 7, whereby the most abundantly expressed allele in the normal tissue is the lowest expressed allele in the tumor and vice versa. We established that this allelic switch does not result from loss of imprinting. The allelic switching of SVOPL may be a result of transcriptional downregulation, while the exact mechanisms resulting in the allelic switching of DLX5 and GRB10 remain to be elucidated. In conclusion, our results show that profound changes take place in allelic transcriptional regulation during the tumorigenesis of CRC.
Collapse
|
47
|
Greenwald WW, Li H, Benaglio P, Jakubosky D, Matsui H, Schmitt A, Selvaraj S, D'Antonio M, D'Antonio-Chronowska A, Smith EN, Frazer KA. Subtle changes in chromatin loop contact propensity are associated with differential gene regulation and expression. Nat Commun 2019; 10:1054. [PMID: 30837461 PMCID: PMC6401380 DOI: 10.1038/s41467-019-08940-5] [Citation(s) in RCA: 70] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Accepted: 02/04/2019] [Indexed: 12/13/2022] Open
Abstract
While genetic variation at chromatin loops is relevant for human disease, the relationships between contact propensity (the probability that loci at loops physically interact), genetics, and gene regulation are unclear. We quantitatively interrogate these relationships by comparing Hi-C and molecular phenotype data across cell types and haplotypes. While chromatin loops consistently form across different cell types, they have subtle quantitative differences in contact frequency that are associated with larger changes in gene expression and H3K27ac. For the vast majority of loci with quantitative differences in contact frequency across haplotypes, the changes in magnitude are smaller than those across cell types; however, the proportional relationships between contact propensity, gene expression, and H3K27ac are consistent. These findings suggest that subtle changes in contact propensity have a biologically meaningful role in gene regulation and could be a mechanism by which regulatory genetic variants in loop anchors mediate effects on expression.
Collapse
Affiliation(s)
- William W Greenwald
- Bioinformatics and Systems Biology Graduate Program, University of California, San Diego, La Jolla, CA, 92093, USA
| | - He Li
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, 92093, USA
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Paola Benaglio
- Department of Pediatrics and Rady Children's Hospital, University of California, San Diego, La Jolla, CA, 92093, USA
| | - David Jakubosky
- Biomedical Sciences Graduate Program, University of California, San Diego, La Jolla, CA, 92093, USA
- Department of Biomedical Sciences, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Hiroko Matsui
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, 92093, USA
| | | | | | - Matteo D'Antonio
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, 92093, USA
- Moores Cancer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | | | - Erin N Smith
- Department of Pediatrics and Rady Children's Hospital, University of California, San Diego, La Jolla, CA, 92093, USA.
| | - Kelly A Frazer
- Institute for Genomic Medicine, University of California, San Diego, La Jolla, CA, 92093, USA.
- Department of Pediatrics and Rady Children's Hospital, University of California, San Diego, La Jolla, CA, 92093, USA.
| |
Collapse
|
48
|
Abstract
Allele-specific expression arises when transcriptional activity at the different alleles of a gene differs considerably. Although extensive research has been carried out to detect and characterize this phenomenon, the landscape of allele-specific expression in cancer is still poorly understood. In this chapter, we describe a fast and reliable analysis pipeline to study allele-specific expression in cancer using next-generation sequencing data. The pipeline provides a gene-level analysis approach that exploits paired germline DNA and tumor RNA sequencing data and benefits from parallel computation resources when available.
Collapse
Affiliation(s)
- Alessandro Romanel
- Centre for Integrative Biology (CIBIO), University of Trento, Trento, Italy.
| |
Collapse
|
49
|
Esfeld K, Berardi AE, Moser M, Bossolini E, Freitas L, Kuhlemeier C. Pseudogenization and Resurrection of a Speciation Gene. Curr Biol 2018; 28:3776-3786.e7. [PMID: 30472000 DOI: 10.1016/j.cub.2018.10.019] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Revised: 09/05/2018] [Accepted: 10/05/2018] [Indexed: 12/13/2022]
Abstract
A persistent question in evolutionary biology is how complex phenotypes evolve and whether phenotypic transitions are reversible. Multiple losses of floral pigmentation have been documented in the angiosperms, but color re-gain has not yet been described, supporting that re-gain is unlikely. Pollinator-mediated selection in Petunia has resulted in several color shifts comprised of both losses and gains of color. The R2R3-MYB transcription factor AN2 has been identified as a major locus responsible for shifts in pollinator preference. Whereas the loss of visible color has previously been attributed to repeated pseudogenization of AN2, here, we describe the mechanism of an independent re-gain of floral color via AN2 evolution. In P. secreta, purple color is restored through the improbable resurrection of AN2 gene function from a non-functional AN2-ancestor by a single reading-frame-restoring mutation. Thus, floral color evolution in Petunia is mechanistically dependent on AN2 functionality, highlighting its role as a hotspot in color transitions and a speciation gene for the genus.
Collapse
Affiliation(s)
- Korinna Esfeld
- Institute of Plant Sciences, University of Bern, Altenbergrain 21, 3013 Bern, Switzerland
| | - Andrea E Berardi
- Institute of Plant Sciences, University of Bern, Altenbergrain 21, 3013 Bern, Switzerland
| | - Michel Moser
- Institute of Plant Sciences, University of Bern, Altenbergrain 21, 3013 Bern, Switzerland
| | - Eligio Bossolini
- Institute of Plant Sciences, University of Bern, Altenbergrain 21, 3013 Bern, Switzerland
| | - Loreta Freitas
- Department of Genetics, University Fed Rio Grande do Sul, POB 15053, Porto Alegre, 91501970 Rio Grande do Sul, Brazil
| | - Cris Kuhlemeier
- Institute of Plant Sciences, University of Bern, Altenbergrain 21, 3013 Bern, Switzerland.
| |
Collapse
|
50
|
Behavior-dependent cis regulation reveals genes and pathways associated with bower building in cichlid fishes. Proc Natl Acad Sci U S A 2018; 115:E11081-E11090. [PMID: 30397142 PMCID: PMC6255178 DOI: 10.1073/pnas.1810140115] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Many behaviors are associated with heritable genetic variation [Kendler and Greenspan (2006) Am J Psychiatry 163:1683-1694]. Genetic mapping has revealed genomic regions or, in a few cases, specific genes explaining part of this variation [Bendesky and Bargmann (2011) Nat Rev Gen 12:809-820]. However, the genetic basis of behavioral evolution remains unclear. Here we investigate the evolution of an innate extended phenotype, bower building, among cichlid fishes of Lake Malawi. Males build bowers of two types, pits or castles, to attract females for mating. We performed comparative genome-wide analyses of 20 bower-building species and found that these phenotypes have evolved multiple times with thousands of genetic variants strongly associated with this behavior, suggesting a polygenic architecture. Remarkably, F1 hybrids of a pit-digging and a castle-building species perform sequential construction of first a pit and then a castle bower. Analysis of brain gene expression in these hybrids showed that genes near behavior-associated variants display behavior-dependent allele-specific expression with preferential expression of the pit-digging species allele during pit digging and of the castle-building species allele during castle building. These genes are highly enriched for functions related to neurodevelopment and neural plasticity. Our results suggest that natural behaviors are associated with complex genetic architectures that alter behavior via cis-regulatory differences whose effects on gene expression are specific to the behavior itself.
Collapse
|