Magnetique: an interactive web application to explore transcriptome signatures of heart failure
Tóm tắt
Despite a recent increase in the number of RNA-seq datasets investigating heart failure (HF), accessibility and usability remain critical issues for medical researchers. We address the need for an intuitive and interactive web application to explore the transcriptional signatures of heart failure with this work. We reanalysed the Myocardial Applied Genomics Network RNA-seq dataset, one of the largest publicly available datasets of left ventricular RNA-seq samples from patients with dilated (DCM) or hypertrophic (HCM) cardiomyopathy, as well as unmatched non-failing hearts (NFD) from organ donors and patient characteristics that allowed us to model confounding factors. We analyse differential gene expression, associated pathway signatures and reconstruct signaling networks based on inferred transcription factor activities through integer linear programming. We additionally focus, for the first time, on differential RNA transcript isoform usage (DTU) changes and predict RNA-binding protein (RBP) to target transcript interactions using a Global test approach. We report results for all pairwise comparisons (DCM, HCM, NFD). Focusing on the DCM versus HCM contrast (DCMvsHCM), we identified 201 differentially expressed genes, some of which can be clearly associated with changes in ERK1 and ERK2 signaling. Interestingly, the signs of the predicted activity for these two kinases have been inferred to be opposite to each other: In the DCMvsHCM contrast, we predict ERK1 to be consistently less activated in DCM while ERK2 was more activated in DCM. In the DCMvsHCM contrast, we identified 149 differently used transcripts. One of the top candidates is the O-linked N-acetylglucosamine (GlcNAc) transferase (OGT), which catalyzes a common post-translational modification known for its role in heart arrhythmias and heart hypertrophy. Moreover, we reconstruct RBP – target interaction networks and showcase the examples of CPEB1, which is differentially expressed in the DCMvsHCM contrast. Magnetique (
https://shiny.dieterichlab.org/app/magnetique
) is the first online application to provide an interactive view of the HF transcriptome at the RNA isoform level and to include transcription factor signaling and RBP:RNA interaction networks. The source code for both the analyses (
https://github.com/dieterich-lab/magnetiqueCode2022
) and the web application (
https://github.com/AnnekathrinSilvia/magnetique
) is available to the public. We hope that our application will help users to uncover the molecular basis of heart failure.
Tài liệu tham khảo
Savarese G, Becher PM, Lund LH, Seferovic P, Rosano GMC, Coats A. Global burden of heart failure: A comprehensive and updated review of epidemiology. Cardiovasc Res. 2022 Feb 12.
Liu C-F, Ni Y, Moravec CS, Morley M, Ashley EA, Cappola TP, et al. Whole-Transcriptome Profiling of Human Heart Tissues Reveals the Potential Novel Players and Regulatory Networks in Different Cardiomyopathy Subtypes of Heart Failure. Circ Genom Precis Med. 2021 Feb 1;14(1):e003142.
Liu Y, Morley M, Brandimarto J, Hannenhalli S, Hu Y, Ashley EA, et al. RNA-Seq identifies novel myocardial gene expression signatures of heart failure. Genomics. 2015 Feb;105(2):83–9.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 2013 Jan;41(Database issue):D991-5.
Wilks C, Zheng SC, Chen FY, Charles R, Solomon B, Ling JP, et al. recount3: summaries and queries for large-scale RNA-seq expression and splicing. Genome Biol. 2021 Nov 29;22(1):323.
Roehr JT, Dieterich C, Reinert K. Flexbar 3.0 - SIMD and multicore parallelization. Bioinformatics. 2017 Sep 15;33(18):2941–2.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012 Mar;9(4)(4):357–9.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013 Jan 1;29(1):15–21.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009 Aug 15;25(16):2078–9.
Kovaka S, Zimin AV, Pertea GM, Razaghi R, Salzberg SL, Pertea M. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol. 2019 Dec;16(1):278. 20(.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
PVCA. https://www.bioconductor.org/packages/release/bioc/html/pvca.html. [Internet]. Accessed 7 Jun 2022.
Alexa A, Rahnenführer J, Lengauer T. Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics. 2006 Jul 1;22(13):1600–7.
Marini F, Binder H. pcaExplorer: an R/Bioconductor package for interacting with RNA-seq principal components. BMC Bioinformatics. 2019 Jun 13;20(1):331.
Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Statist. 2001 Aug;29(4):1165–88.
Marini F, Ludt A, Linke J, Strauch K. GeneTonic: an R/Bioconductor package for streamlining the interpretation of RNA-seq data. BMC Bioinformatics. 2021 Dec 23;22(1):610.
Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. [version 2; peer review: 2 approved]. F1000Res. 2015 Jan 1;4:1521.
Soneson C, Matthes KL, Nowicka M, Law CW, Robinson MD. Isoform prefiltering improves performance of count-based methods for analysis of differential transcript usage. Genome Biol. 2016 Jan 26;17:12.
Nowicka M, Robinson MD. DRIMSeq: a Dirichlet-multinomial framework for multivariate count outcomes in genomics. [version 2; peer review: 2 approved]. F1000Res. 2016 Jun 13;5:1356.
Goeman JJ, van de Geer SA, de Kort F, van Houwelingen HC. A global test for groups of genes: testing association with a clinical outcome. Bioinformatics. 2004 Jan;20(1)(1):93–9.
Benoit Bouvrette LP, Bovaird S, Blanchette M, Lécuyer E. oRNAment: a database of putative RNA binding protein target sites in the transcriptomes of model species. Nucleic Acids Res. 2020 Jan 8;48(D1):D166–73.
van Iterson M, Bervoets S, de Meijer EJ, Buermans HP, ’t Hoen PAC, Menezes RX, et al. Integrated analysis of microRNA and mRNA expression: adding biological significance to microRNA target predictions. Nucleic Acids Res. 2013 Aug;41(15):e146.
Liu A, Trairatphisan P, Gjerga E, Didangelos A, Barratt J, Saez-Rodriguez J. From expression footprints to causal pathways: contextualizing large signaling networks with CARNIVAL. NPJ Syst Biol Appl. 2019 Nov 11;5:40.
Türei D, Korcsmáros T, Saez-Rodriguez J. OmniPath: guidelines and gateway for literature-curated signaling pathway resources. Nat Methods. 2016 Nov 29;13(12):966–7.
Garcia-Alonso L, Holland CH, Ibrahim MM, Turei D, Saez-Rodriguez J. Benchmark and integration of resources for the estimation of human transcription factor activities. Genome Res. 2019 Aug;29(8):1363–75.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010 Jan 1;26(1):139–40.
Alvarez MJ, Shen Y, Giorgi FM, Lachmann A, Ding BB, Ye BH, et al. Functional characterization of somatic mutations in cancer using network-based inference of protein activity. Nat Genet. 2016 Aug;48(8):838–47.
Iorio F, Bernardo-Faura M, Gobbi A, Cokelaer T, Jurman G, Saez-Rodriguez J. Efficient randomization of biological networks while preserving functional characterization of individual nodes. BMC Bioinformatics. 2016 Dec 20;17(1):542.
ShinyProxy - Open. Source Enterprise Deployment for Shiny https://github.com/openanalytics/shinyproxy. Accessed 30 June 2022.
Morgan M, Obenchain V, Hester J, Pagès H. SummarizedExperiment: SummarizedExperiment container. R package version. Bioconductor; 2022.
Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS ONE. 2010 Nov 15;5(11):e13984.
Vitting-Seerup K, Sandelin A. The landscape of isoform switches in human cancers. Mol Cancer Res. 2017 Sep;15(9):1206–20.
Kong SW, Hu YW, Ho JWK, Ikeda S, Polster S, John R, et al. Heart failure-associated changes in RNA splicing of sarcomere genes. Circ Cardiovasc Genet. 2010 Apr;3(2):138–46.
Chen Q, Chen Y, Bian C, Fujiki R, Yu X. TET2 promotes histone O-GlcNAcylation during gene transcription. Nature. 2013 Jan 24;493(7433):561–4.
Erickson JR, Pereira L, Wang L, Han G, Ferguson A, Dao K, et al. Diabetic hyperglycaemia activates CaMKII and arrhythmias by O-linked glycosylation. Nature. 2013 Oct 17;502(7471):372–6.
Mu Y, Yu H, Wu T, Zhang J, Evans SM, Chen J. O-linked β-N-acetylglucosamine transferase plays an essential role in heart development through regulating angiopoietin-1. PLoS Genet. 2020 Apr 6;16(4):e1008730.
Park S-K, Zhou X, Pendleton KE, Hunter OV, Kohler JJ, O’Donnell KA, et al. A Conserved Splicing Silencer Dynamically Regulates O-GlcNAc Transferase Intron Retention and O-GlcNAc Homeostasis. Cell Rep. 2017 Aug 1;20(5):1088–99.
Riechert E, Kmietczyk V, Stein F, Schwarzl T, Sekaran T, Jürgensen L, et al. Identification of dynamic RNA-binding proteins uncovers a Cpeb4-controlled regulatory cascade during pathological cell growth of cardiomyocytes. Cell Rep. 2021 May 11;35(6):109100.
Ekure EN, Adeyemo A, Liu H, Sokunbi O, Kalu N, Martinez AF, et al. Exome Sequencing and Congenital Heart Disease in Sub-Saharan Africa. Circ Genom Precis Med. 2021 Feb;14(1):e003108.
Bava F-A, Eliscovich C, Ferreira PG, Miñana B, Ben-Dov C, Guigó R, et al. CPEB1 coordinates alternative 3’-UTR formation with translational regulation. Nature. 2013 Mar 7;495(7439):121–5.
Radke MH, Badillo-Lisakowski V, Britto-Borges T, Kubli DA, Jüttner R, Parakkat P, et al. Therapeutic inhibition of RBM20 improves diastolic function in a murine heart failure model and human engineered heart tissue. Sci Transl Med. 2021 Dec 1;13(622):eabe8952.
Kleinberg JM. Authoritative sources in a hyperlinked environment. J ACM (JACM). 1999 Sep;46(5)(1):604–32.
Gilbert CJ, Longenecker JZ, Accornero F. ERK1/2: an integrator of signals that alters cardiac homeostasis and growth. Biology (Basel). 2021 Apr 20;10(4).
Zhang Z, Li S, Cui M, Gao X, Sun D, Qin X, et al. Rosuvastatin enhances the therapeutic efficacy of adipose-derived mesenchymal stem cells for myocardial infarction via PI3K/Akt and MEK/ERK pathways. Basic Res Cardiol. 2013 Mar;108(2):333.
Backs J, Backs T, Bezprozvannaya S, McKinsey TA, Olson EN. Histone deacetylase 5 acquires calcium/calmodulin-dependent kinase II responsiveness by oligomerization with histone deacetylase 4. Mol Cell Biol. 2008 May;28(10):3437–45.
Müller R, Weirick T, John D, Militello G, Chen W, Dimmeler S, et al. ANGIOGENES: knowledge database for protein-coding and noncoding RNA genes in endothelial cells. Sci Rep. 2016 Sep;1:6:32475.
Ramirez Flores RO, Lanzer JD, Holland CH, Leuschner F, Most P, Schultz J-H, et al. Consensus Transcriptional Landscape of Human End-Stage Heart Failure. J Am Heart Assoc. 2021 Apr 6;10(7):e019667.
Gabisonia K, Prosdocimo G, Aquaro GD, Carlucci L, Zentilin L, Secco I, et al. MicroRNA therapy stimulates uncontrolled cardiac repair after myocardial infarction in pigs. Nature. 2019 May 8;569(7756):418–22.
Heinig M, Adriaens ME, Schafer S, van Deutekom HWM, Lodder EM, Ware JS, et al. Natural genetic variation of the cardiac transcriptome in non-diseased donors and patients with dilated cardiomyopathy. Genome Biol. 2017 Sep;14(1):170. 18(.