Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3′READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

[1]  Christopher W. J. Smith,et al.  Alternative splicing as a source of phenotypic diversity , 2022, Nature Reviews Genetics.

[2]  Shengqing Li,et al.  Mxi1-0 Promotes Hypoxic Pulmonary Hypertension Via ERK/c-Myc-dependent Proliferation of Arterial Smooth Muscle Cells , 2022, Frontiers in Genetics.

[3]  Xiang-Dong Fu,et al.  Alternative polyadenylation by sequential activation of distal and proximal PolyA sites , 2022, Nature Structural & Molecular Biology.

[4]  B. Shiotani,et al.  PRPF19 regulates p53-dependent cellular senescence by modulating alternative splicing of MDM4 mRNA , 2021, The Journal of biological chemistry.

[5]  B. Suess,et al.  SRSF3 and SRSF7 modulate 3′UTR length through suppression or activation of proximal polyadenylation sites and regulation of CFIm levels , 2021, Genome biology.

[6]  Shou-Jiang Gao,et al.  RNF167 activates mTORC1 and promotes tumorigenesis by targeting CASTOR1 for ubiquitination and degradation , 2021, Nature Communications.

[7]  P. Tan,et al.  Long-read transcriptome sequencing reveals abundant promoter diversity in distinct molecular subtypes of gastric cancer , 2021, Genome biology.

[8]  Byul A. Jee,et al.  Global spliceosome activity regulates entry into cellular senescence , 2020, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[9]  W. Sadee,et al.  Interpreting coronary artery disease GWAS results: A functional genomics approach assessing biological significance , 2020, bioRxiv.

[10]  S. Jagadish,et al.  Utilizing PacBio Iso-Seq for Novel Transcript and Gene Discovery of Abiotic Stress Responses in Oryza sativa L. , 2020, International journal of molecular sciences.

[11]  Yong Xia,et al.  Angiotensin II Decreases Endothelial Nitric Oxide Synthase Phosphorylation via AT1R Nox/ROS/PP2A Pathway , 2020, Frontiers in Physiology.

[12]  W. Hong,et al.  Histone methyltransferase Smyd3 is a new regulator for vascular senescence , 2020, Aging cell.

[13]  David P. Kreil,et al.  A comprehensive rat transcriptome built from large scale RNA-seq-based annotation , 2020, Nucleic acids research.

[14]  E. Hoffman,et al.  A long-read RNA-seq approach to identify novel transcripts of very large genes , 2020, Genome research.

[15]  W. Timens,et al.  Link between increased cellular senescence and extracellular matrix changes in COPD. , 2020, American journal of physiology. Lung cellular and molecular physiology.

[16]  Barbara J. Wold,et al.  A technology-agnostic long-read analysis pipeline for transcriptome discovery and quantification , 2019, bioRxiv.

[17]  S. Dutta,et al.  Alternative splicing in aging and longevity , 2019, Human Genetics.

[18]  E. Schuman,et al.  Full-length transcriptome reconstruction reveals a large diversity of RNA and protein isoforms in rat hippocampus , 2019, Nature Communications.

[19]  Gunnar Rätsch,et al.  A Pan-cancer Transcriptome Analysis Reveals Pervasive Regulation through Alternative Promoters , 2019, Cell.

[20]  Jianchun Cai,et al.  Tumor endothelial cell-derived cadherin-2 promotes angiogenesis and has prognostic significance for lung adenocarcinoma , 2019, Molecular Cancer.

[21]  Jinzhong Lin,et al.  Alternative polyadenylation dependent function of splicing factor SRSF3 contributes to cellular senescence , 2019, Aging.

[22]  D. Leckband,et al.  N-cadherin signaling via Trio assembles adherens junctions to restrict endothelial permeability , 2018, The Journal of cell biology.

[23]  M. Alique,et al.  Endothelial Cell Senescence in the Pathogenesis of Endothelial Dysfunction , 2018, Endothelial Dysfunction - Old Concepts and New Challenges.

[24]  G. Booz,et al.  Angiotensin II Signal Transduction: An Update on Mechanisms of Physiology and Pathophysiology. , 2018, Physiological reviews.

[25]  Xinhua Liu,et al.  HDAC4 regulates vascular inflammation via activation of autophagy , 2018, Cardiovascular research.

[26]  Xiaohui Wu,et al.  APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data , 2018, Bioinform..

[27]  Gloria M. Sheynkman,et al.  Full-length mRNA sequencing uncovers a widespread coupling between transcription initiation and mRNA processing , 2018, Genome Biology.

[28]  Q. Morris,et al.  QAPA: a new method for the systematic analysis of alternative polyadenylation from RNA-seq data , 2018, Genome Biology.

[29]  Miha Skalic,et al.  SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions , 2016, Genome Biology.

[30]  Tao Jiang,et al.  TAPAS: tool for alternative polyadenylation site analysis , 2018, Bioinform..

[31]  Wei Li,et al.  3′ UTR lengthening as a novel mechanism in regulating cellular senescence , 2018, Genome research.

[32]  L. McMillan,et al.  FMLRC: Hybrid long read error correction using an FM-index , 2018, BMC Bioinformatics.

[33]  Bin Tian,et al.  PolyA_DB 3 catalogs cleavage and polyadenylation sites identified by deep sequencing in multiple genomes , 2017, Nucleic Acids Res..

[34]  Heng Li,et al.  Minimap2: pairwise alignment for nucleotide sequences , 2017, Bioinform..

[35]  Wei Chen,et al.  Alternative Polyadenylation: Methods, Findings, and Impacts , 2017, Genom. Proteom. Bioinform..

[36]  J. Steitz,et al.  Comparative analysis reveals genomic features of stress-induced transcriptional readthrough , 2017, Proceedings of the National Academy of Sciences.

[37]  T. Tolvanen,et al.  Inhibition of SHIP2 in CD2AP-deficient podocytes ameliorates reactive oxygen species generation but aggravates apoptosis , 2017, Scientific Reports.

[38]  F. Baralle,et al.  Alternative splicing as a regulator of development and tissue identity , 2017, Nature Reviews Molecular Cell Biology.

[39]  T. Wieland,et al.  Identification of Cadherin 2 (CDH2) Mutations in Arrhythmogenic Right Ventricular Cardiomyopathy , 2017, Circulation. Cardiovascular genetics.

[40]  Kari Stefansson,et al.  Genetic variants associated with mosaic Y chromosome loss highlight cell cycle genes and overlap with cancer susceptibility , 2017, Nature Genetics.

[41]  Lennart Martens,et al.  1 SQANTI : extensive characterization of long read transcript sequences for quality control in 1 full-length transcriptome identification and quantification 2 3 , 2017 .

[42]  Geet Duggal,et al.  Salmon: fast and bias-aware quantification of transcript expression using dual-phase inference , 2017, Nature Methods.

[43]  B. Tian,et al.  Alternative polyadenylation of mRNA precursors , 2016, Nature Reviews Molecular Cell Biology.

[44]  Geng Wu,et al.  Structural mechanism for the arginine sensing and regulation of CASTOR1 in the mTORC1 signaling pathway , 2016, Cell Discovery.

[45]  Faye D. Schilkey,et al.  A survey of the sorghum transcriptome using single-molecule long reads , 2016, Nature Communications.

[46]  Lihua Huang,et al.  Mechanisms of resistance to EGFR tyrosine kinase inhibitors , 2015, Acta pharmaceutica Sinica. B.

[47]  Xiandong Meng,et al.  Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing , 2015, PloS one.

[48]  William Stafford Noble,et al.  The MEME Suite , 2015, Nucleic Acids Res..

[49]  Derek W Wright,et al.  Gateways to the FANTOM5 promoter level mammalian expression atlas , 2015, Genome Biology.

[50]  Wei Li,et al.  Dynamic analyses of alternative polyadenylation from RNA-seq reveal a 3′-UTR landscape across seven tumour types , 2014, Nature Communications.

[51]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[52]  Thomas Hackl,et al.  proovread: large-scale high-accuracy PacBio correction through iterative short read consensus , 2014, Bioinform..

[53]  Alexander G Williams,et al.  RNA‐seq Data: Challenges in and Recommendations for Experimental Design and Analysis , 2014, Current protocols in human genetics.

[54]  Alexandre Lomsadze,et al.  Identification of protein coding regions in RNA transcripts , 2014, BCB.

[55]  Leena Salmela,et al.  LoRDEC: accurate and efficient long read error correction , 2014, Bioinform..

[56]  Hongzhe Li,et al.  A change-point model for identifying 3′UTR switching by next-generation RNA sequencing , 2014, Bioinform..

[57]  R. Touyz,et al.  Angiotensin II, NADPH oxidase, and redox signaling in the vasculature. , 2013, Antioxidants & redox signaling.

[58]  R. Elkon,et al.  Alternative cleavage and polyadenylation: extent, regulation and function , 2013, Nature Reviews Genetics.

[59]  Aaron A. Klammer,et al.  Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data , 2013, Nature Methods.

[60]  Larisa M Haupt,et al.  Review: Alternative Splicing (AS) of Genes As An Approach for Generating Protein Complexity , 2013, Current genomics.

[61]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[62]  C. Rembold,et al.  Free radical biology of the cardiovascular system. , 2012, Clinical science.

[63]  T. Babak,et al.  A quantitative atlas of polyadenylation in five mammals , 2012, Genome research.

[64]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[65]  P. Pandolfi,et al.  The functions and regulation of the PTEN tumour suppressor , 2012, Nature Reviews Molecular Cell Biology.

[66]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[67]  F. Blanco,et al.  The Splicing Factor SRSF1 as a Marker for Endothelial Senescence , 2012, Front. Physio..

[68]  Ji Min Kim,et al.  Mechanism of Ang II involvement in activation of NF-κB through phosphorylation of p65 during aging , 2012, AGE.

[69]  A. Hallberg,et al.  Angiotensin II, a Neuropeptide at the Frontier between Endocrinology and Neuroscience: Is There a Link between the Angiotensin II Type 2 Receptor and Alzheimer’s Disease? , 2011, Front. Endocrin..

[70]  Marcel Martin Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[71]  Michael Bader,et al.  Rat models of cardiovascular diseases. , 2010, Methods in molecular biology.

[72]  M. Gerstein,et al.  RNA-Seq: a revolutionary tool for transcriptomics , 2009, Nature Reviews Genetics.

[73]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[74]  H. Shan,et al.  Angiotensin II induces endothelial cell senescence via the activation of mitogen‐activated protein kinases , 2008, Cell biochemistry and function.

[75]  G. Tortora,et al.  EGFR antagonists in cancer treatment. , 2008, The New England journal of medicine.

[76]  Pearlly Yan,et al.  Genome-wide analysis of alternative promoters of human genes using a custom promoter tiling array , 2008, BMC Genomics.

[77]  H. Kajiyama,et al.  Functional Expression of the Angiotensin II Type1 Receptor in Human Ovarian Carcinoma Cells and Its Blockade Therapy Resulting in Suppression of Tumor Invasion, Angiogenesis, and Peritoneal Dissemination , 2005, Clinical Cancer Research.

[78]  A. Youkilis,et al.  Mxi1-0, an alternatively transcribed Mxi1 isoform, is overexpressed in glioblastomas. , 2004, Neoplasia.

[79]  Dixie L Mager,et al.  Complex controls: the role of alternative promoters in mammalian genomes. , 2003, Trends in genetics : TIG.

[80]  Chen Yan,et al.  Angiotensin II signaling pathways mediated by tyrosine kinases. , 2003, The international journal of biochemistry & cell biology.

[81]  Yuji Saito,et al.  Angiotensin II-mediated signal transduction pathways , 2002, Current hypertension reports.

[82]  A. G. Scicli,et al.  Angiotensin II stimulates migration of retinal microvascular pericytes: involvement of TGF-beta and PDGF-BB. , 2002, American journal of physiology. Heart and circulatory physiology.

[83]  B. Graveley Alternative splicing: increasing diversity in the proteomic world. , 2001, Trends in genetics : TIG.

[84]  S. Doggrell,et al.  Rat models of hypertension, cardiac hypertrophy and failure. , 1998, Cardiovascular research.

[85]  R. DePinho,et al.  Role of Mxi1 in ageing organ systems and the regulation of normal and neoplastic growth , 1998, Nature.

[86]  T. Ayoubi,et al.  Regulation of gene expression by alternative promoters , 1996, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[87]  J. Wilusz,et al.  Cleavage site determinants in the mammalian polyadenylation signal. , 1995, Nucleic acids research.