51
|
Agarwal V, Lopez-Darwin S, Kelley DR, Shendure J. The landscape of alternative polyadenylation in single cells of the developing mouse embryo. Nat Commun 2021; 12:5101. [PMID: 34429411 PMCID: PMC8385098 DOI: 10.1038/s41467-021-25388-8] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Accepted: 08/06/2021] [Indexed: 02/07/2023] Open
Abstract
3′ untranslated regions (3′ UTRs) post-transcriptionally regulate mRNA stability, localization, and translation rate. While 3′-UTR isoforms have been globally quantified in limited cell types using bulk measurements, their differential usage among cell types during mammalian development remains poorly characterized. In this study, we examine a dataset comprising ~2 million nuclei spanning E9.5–E13.5 of mouse embryonic development to quantify transcriptome-wide changes in alternative polyadenylation (APA). We observe a global lengthening of 3′ UTRs across embryonic stages in all cell types, although we detect shorter 3′ UTRs in hematopoietic lineages and longer 3′ UTRs in neuronal cell types within each stage. An analysis of RNA-binding protein (RBP) dynamics identifies ELAV-like family members, which are concomitantly induced in neuronal lineages and developmental stages experiencing 3′-UTR lengthening, as putative regulators of APA. By measuring 3′-UTR isoforms in an expansive single cell dataset, our work provides a transcriptome-wide and organism-wide map of the dynamic landscape of alternative polyadenylation during mammalian organogenesis. Alternative polyadenylation regulates localization, half-life and translation of mRNA isoforms. Here the authors investigate alternative polyadenylation using single cell RNA sequencing data from mouse embryos and identify 3’-UTR isoforms that are regulated across cell types and developmental time.
Collapse
Affiliation(s)
| | | | | | - Jay Shendure
- Department of Genome Sciences, University of Washington, Seattle, WA, USA. .,Howard Hughes Medical Institute, Seattle, WA, USA. .,Brotman Baty Institute for Precision Medicine, University of Washington, Seattle, WA, USA. .,Allen Discovery Center for Cell Lineage Tracing, Seattle, WA, USA.
| |
Collapse
|
52
|
He Y, Yu H, Zhao H, Zhu H, Zhang Q, Wang A, Shen Y, Xu X, Li J. Transcriptomic analysis to elucidate the effects of high stocking density on grass carp (Ctenopharyngodon idella). BMC Genomics 2021; 22:620. [PMID: 34399686 PMCID: PMC8369720 DOI: 10.1186/s12864-021-07924-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Accepted: 08/06/2021] [Indexed: 01/23/2023] Open
Abstract
Background Grass carp (Ctenopharyngodon idella) is one of the most widely cultivated fishes in China. High stocking density can reportedly affect fish growth and immunity. Herein we performed PacBio long-read single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing to evaluate the effects of high stocking density on grass carp transcriptome. Results SMRT sequencing led to the identification of 33,773 genes (14,946 known and 18,827 new genes). From the structure analysis, 8,009 genes were detected with alternative splicing events, 10,219 genes showed alternative polyadenylation sites and 15,521 long noncoding RNAs. Further, 1,235, 962, and 213 differentially expressed genes (DEGs) were identified in the intestine, muscle, and brain tissues, respectively. We performed functional enrichment analyses of DEGs, and they were identified to be significantly enriched in nutrient metabolism and immune function. The expression levels of several genes encoding apolipoproteins and activities of enzymes involved in carbohydrate enzymolysis were found to be upregulated in the high stocking density group, indicating that lipid metabolism and carbohydrate decomposition were accelerated. Besides, four isoforms of grass carp major histocompatibility complex class II antigen alpha and beta chains in the aforementioned three tissue was showed at least a 4-fold decrease. Conclusions The results suggesting that fish farmed at high stocking densities face issues associated with the metabolism and immune system. To conclude, our results emphasize the importance of maintaining reasonable density in grass carp aquaculture. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07924-4.
Collapse
Affiliation(s)
- Yan He
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China.,National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China.,Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China
| | - Hongyan Yu
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China.,National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China
| | - Honggang Zhao
- Department of Natural Resources, Cornell University, 14853, Ithaca, New York, USA
| | - Hua Zhu
- Beijing Key Laboratory of Fishery Biotechnology, Beijing Fisheries Research Institute, 100068, Beijing, China
| | - Qingjing Zhang
- Beijing Key Laboratory of Fishery Biotechnology, Beijing Fisheries Research Institute, 100068, Beijing, China
| | - Anqi Wang
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China
| | - Yubang Shen
- National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China
| | - Xiaoyan Xu
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China. .,National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China. .,Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China.
| | - Jiale Li
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China. .,National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China. .,Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China.
| |
Collapse
|
53
|
Van Essen D, Alcaraz AJG, Miller JGP, Jones PD, Doering JA, Wiseman S. The brominated flame retardant, TBCO, impairs oocyte maturation in zebrafish (Danio rerio). AQUATIC TOXICOLOGY (AMSTERDAM, NETHERLANDS) 2021; 238:105929. [PMID: 34375885 DOI: 10.1016/j.aquatox.2021.105929] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Revised: 07/15/2021] [Accepted: 07/27/2021] [Indexed: 06/13/2023]
Abstract
The brominated flame retardant, 1,2,5,6-tetrabromocyclooctane (TBCO), has been shown to decrease fecundity in Japanese medaka (Oryzias latipes) and there is indirect evidence from analysis of the transcriptome and proteome that this effect might be due to impaired oogenesis. An assay for disruption of oocyte maturation by chemical stressors has not been developed in Japanese medaka. Thus, using zebrafish (Danio rerio) as a model, objectives of the present study were to determine whether exposure to TBCO has effects on maturation of oocytes and to investigate potential mechanisms. Sexually mature female zebrafish were given a diet of 35.3 or 628.8 μg TBCO / g food for 14 days after which, stage IV oocytes were isolated to assess maturation in response to maturation inducing hormone. To explore potential molecular mechanisms, abundances of mRNAs of a suite of genes that regulate oocyte maturation were quantified by use of quantitative real-time PCR, and abundances of microRNAs were determined by use of miRNAseq. Ex vivo maturation of oocytes from fish exposed to TBCO was significantly less than maturation of oocytes from control fish. The percentage of oocytes which matured from control fish and those exposed to low and high TBCO were 89, 71, and 67%, respectively. Among the suite of genes known to regulate oocyte maturation, mRNA abundance of insulin like growth factor-3 was decreased by 1.64- and 3.44-fold in stage IV oocytes from females given the low and high concentrations of TBCO, respectively, compared to the control group. Abundances of microRNAs regulating the expression of proteins that regulate oocyte maturation, including processes related to insulin-like growth factor, were significantly different in stage IV oocytes from fish exposed to TBCO. Overall, results of this study indicated that impaired oocyte maturation might be a mechanism of reduced reproductive performance in TBCO-exposed fish. Results also suggested that effects of TBCO on oocyte maturation might be due to molecular perturbations on insulin-like growth factor signaling and expression of microRNAs.
Collapse
Affiliation(s)
- Darren Van Essen
- Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, T1K 3M4, Canada
| | | | - Justin G P Miller
- Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, T1K 3M4, Canada
| | - Paul D Jones
- Toxicology Centre, University of Saskatchewan, Saskatoon, SK, S7N 5B3, Canada; School of Environment and Sustainability, University of Saskatchewan, Saskatoon, SK, S7N 5B3, Canada
| | - Jon A Doering
- Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, T1K 3M4, Canada; Intersectoral Centre for Endocrine Disruptor Analysis (ICEDA), Institut National de la Recherche Scientifique (INRS), Centre Eau Terre Environnement, Québec City, QC, G1K 9A9, Canada
| | - Steve Wiseman
- Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, T1K 3M4, Canada; Intersectoral Centre for Endocrine Disruptor Analysis (ICEDA), Institut National de la Recherche Scientifique (INRS), Centre Eau Terre Environnement, Québec City, QC, G1K 9A9, Canada; Water Institute for Sustainable Environments, University of Lethbridge, Lethbridge, AB, T1K 3M4, Canada.
| |
Collapse
|
54
|
Meng X, Kuang K, Zhang Y, Guan K, Liu B, Zhou X. Alternative polyadenylation events differ dramatically between Tongcheng and Large White pigs in response to PRRSV infection. Anim Genet 2021; 52:744-748. [PMID: 34309053 DOI: 10.1111/age.13125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/05/2021] [Indexed: 12/01/2022]
Abstract
Alternative polyadenylation (APA) is a widespread post-transcriptional regulation mechanism that increases the biological complexity of transcriptome and proteome. However, it is unclear whether APA regulation plays a role in genetic resistance to porcine reproductive and respiratory syndrome virus (PRRSV). Here, we reported genome-wide APA regulation of porcine alveolar macrophages in PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs upon PRRSV infection. Using 3' mRNA sequencing strategy, we detected 75 981 high-quality APA sites in porcine alveolar macrophages of TC and LW pigs. Furthermore, 1202 and 1089 differentially expressed APA sites, as well as 79 and 117 untranslated region-APA switching genes were identified in TC pigs and LW pigs upon PRRSV infection respectively. The APA events in TC pigs and LW pigs were involved in different biological pathways, while APA events in TC pigs are directly associated with the immune response to PRRSV infection. In addition, we identified genetic variations affecting polyadenylation signal between TC pigs and LW pigs. These findings would provide helpful information on APA regulation for further understanding of genetic resistance to PRRSV.
Collapse
Affiliation(s)
- X Meng
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - K Kuang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Y Zhang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - K Guan
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - B Liu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China.,The Engineering Technology Research Center of Hubei Province Local Pig Breed Improvement, Wuhan, 430070, China
| | - X Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China.,The Engineering Technology Research Center of Hubei Province Local Pig Breed Improvement, Wuhan, 430070, China
| |
Collapse
|
55
|
Binzel DW, Li X, Burns N, Khan E, Lee WJ, Chen LC, Ellipilli S, Miles W, Ho YS, Guo P. Thermostability, Tunability, and Tenacity of RNA as Rubbery Anionic Polymeric Materials in Nanotechnology and Nanomedicine-Specific Cancer Targeting with Undetectable Toxicity. Chem Rev 2021; 121:7398-7467. [PMID: 34038115 PMCID: PMC8312718 DOI: 10.1021/acs.chemrev.1c00009] [Citation(s) in RCA: 56] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
RNA nanotechnology is the bottom-up self-assembly of nanometer-scale architectures, resembling LEGOs, composed mainly of RNA. The ideal building material should be (1) versatile and controllable in shape and stoichiometry, (2) spontaneously self-assemble, and (3) thermodynamically, chemically, and enzymatically stable with a long shelf life. RNA building blocks exhibit each of the above. RNA is a polynucleic acid, making it a polymer, and its negative-charge prevents nonspecific binding to negatively charged cell membranes. The thermostability makes it suitable for logic gates, resistive memory, sensor set-ups, and NEM devices. RNA can be designed and manipulated with a level of simplicity of DNA while displaying versatile structure and enzyme activity of proteins. RNA can fold into single-stranded loops or bulges to serve as mounting dovetails for intermolecular or domain interactions without external linking dowels. RNA nanoparticles display rubber- and amoeba-like properties and are stretchable and shrinkable through multiple repeats, leading to enhanced tumor targeting and fast renal excretion to reduce toxicities. It was predicted in 2014 that RNA would be the third milestone in pharmaceutical drug development. The recent approval of several RNA drugs and COVID-19 mRNA vaccines by FDA suggests that this milestone is being realized. Here, we review the unique properties of RNA nanotechnology, summarize its recent advancements, describe its distinct attributes inside or outside the body and discuss potential applications in nanotechnology, medicine, and material science.
Collapse
Affiliation(s)
- Daniel W Binzel
- Center for RNA Nanobiotechnology and Nanomedicine, College of Pharmacy, Dorothy M. Davis Heart and Lung Research Institute, James Comprehensive Cancer Center, College of Medicine, The Ohio State University, Columbus, Ohio 43210, United States
| | - Xin Li
- Center for RNA Nanobiotechnology and Nanomedicine, College of Pharmacy, Dorothy M. Davis Heart and Lung Research Institute, James Comprehensive Cancer Center, College of Medicine, The Ohio State University, Columbus, Ohio 43210, United States
| | - Nicolas Burns
- Center for RNA Nanobiotechnology and Nanomedicine, College of Pharmacy, Dorothy M. Davis Heart and Lung Research Institute, James Comprehensive Cancer Center, College of Medicine, The Ohio State University, Columbus, Ohio 43210, United States
| | - Eshan Khan
- Department of Cancer Biology and Genetics, The Ohio State University Comprehensive Cancer Center, College of Medicine, Center for RNA Biology, The Ohio State University, Columbus, Ohio 43210, United States
| | - Wen-Jui Lee
- TMU Research Center of Cancer Translational Medicine, School of Medical Laboratory Science and Biotechnology, College of Medical Science and Technology, Graduate Institute of Medical Sciences, College of Medicine, Taipei Medical University, Department of Laboratory Medicine, Taipei Medical University Hospital, Taipei 110, Taiwan
| | - Li-Ching Chen
- TMU Research Center of Cancer Translational Medicine, School of Medical Laboratory Science and Biotechnology, College of Medical Science and Technology, Graduate Institute of Medical Sciences, College of Medicine, Taipei Medical University, Department of Laboratory Medicine, Taipei Medical University Hospital, Taipei 110, Taiwan
| | - Satheesh Ellipilli
- Center for RNA Nanobiotechnology and Nanomedicine, College of Pharmacy, Dorothy M. Davis Heart and Lung Research Institute, James Comprehensive Cancer Center, College of Medicine, The Ohio State University, Columbus, Ohio 43210, United States
| | - Wayne Miles
- Department of Cancer Biology and Genetics, The Ohio State University Comprehensive Cancer Center, College of Medicine, Center for RNA Biology, The Ohio State University, Columbus, Ohio 43210, United States
| | - Yuan Soon Ho
- TMU Research Center of Cancer Translational Medicine, School of Medical Laboratory Science and Biotechnology, College of Medical Science and Technology, Graduate Institute of Medical Sciences, College of Medicine, Taipei Medical University, Department of Laboratory Medicine, Taipei Medical University Hospital, Taipei 110, Taiwan
| | - Peixuan Guo
- Center for RNA Nanobiotechnology and Nanomedicine, College of Pharmacy, Dorothy M. Davis Heart and Lung Research Institute, James Comprehensive Cancer Center, College of Medicine, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
56
|
Goering R, Engel KL, Gillen AE, Fong N, Bentley DL, Taliaferro JM. LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival. BMC Genomics 2021; 22:476. [PMID: 34174817 PMCID: PMC8234626 DOI: 10.1186/s12864-021-07781-1] [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: 01/14/2021] [Accepted: 06/07/2021] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND The sequence content of the 3' UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging. RESULTS To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated. CONCLUSIONS We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.
Collapse
Affiliation(s)
- Raeann Goering
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
- RNA Bioscience Initiative, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - Krysta L Engel
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - Austin E Gillen
- RNA Bioscience Initiative, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
- Division of Hematology, University of Colorado School of Medicine, Aurora, CO, USA
| | - Nova Fong
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - David L Bentley
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
- RNA Bioscience Initiative, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - J Matthew Taliaferro
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO, USA.
- RNA Bioscience Initiative, University of Colorado Anschutz Medical Campus, Aurora, CO, USA.
| |
Collapse
|
57
|
A CRISPR RNA-binding protein screen reveals regulators of RUNX1 isoform generation. Blood Adv 2021; 5:1310-1323. [PMID: 33656539 DOI: 10.1182/bloodadvances.2020002090] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2020] [Accepted: 01/19/2021] [Indexed: 11/20/2022] Open
Abstract
The proper balance of hematopoietic stem cell (HSC) self-renewal and differentiation is critical for normal hematopoiesis and is disrupted in hematologic malignancy. Among regulators of HSC fate, transcription factors have a well-defined central role, and mutations promote malignant transformation. More recently, studies have illuminated the importance of posttranscriptional regulation by RNA-binding proteins (RBPs) in hematopoiesis and leukemia development. However, the RBPs involved and the breadth of regulation are only beginning to be elucidated. Furthermore, the intersection between posttranscriptional regulation and hematopoietic transcription factor function is poorly understood. Here, we studied the posttranscriptional regulation of RUNX1, a key hematopoietic transcription factor. Alternative polyadenylation (APA) of RUNX1 produces functionally antagonistic protein isoforms (RUNX1a vs RUNX1b/c) that mediate HSC self-renewal vs differentiation, an RNA-processing event that is dysregulated in malignancy. Consequently, RBPs that regulate this event directly contribute to healthy and aberrant hematopoiesis. We modeled RUNX1 APA using a split GFP minigene reporter and confirmed the sensitivity of our model to detect changes in RNA processing. We used this reporter in a clustered regularly interspaced short palindromic repeats (CRISPR) screen consisting of single guide RNAs exclusively targeting RBPs and uncovered HNRNPA1 and KHDRBS1 as antagonistic regulators of RUNX1a isoform generation. Overall, our study provides mechanistic insight into the posttranscriptional regulation of a key hematopoietic transcription factor and identifies RBPs that may have widespread and important functions in hematopoiesis.
Collapse
|
58
|
Mitschka S, Fansler MM, Mayr C. Generation of 3'UTR knockout cell lines by CRISPR/Cas9-mediated genome editing. Methods Enzymol 2021; 655:427-457. [PMID: 34183132 DOI: 10.1016/bs.mie.2021.03.014] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In addition to the protein code, messenger RNAs (mRNAs) also contain untranslated regions (UTRs). 3'UTRs span the region between the translational stop codon and the poly(A) tail. Sequence elements located in 3'UTRs are essential for pre-mRNA processing. 3'UTRs also contain elements that can regulate protein abundance, localization, and function. At least half of all human genes use alternative cleavage and polyadenylation (APA) to further diversify the regulatory potential of protein functions. Traditional gene editing approaches are designed to disrupt the production of functional proteins. Here, we describe a method that allows investigators to manipulate 3'UTR sequences of endogenous genes for both single- 3'UTR and multi-3'UTR genes. As 3'UTRs can regulate individual functions of proteins, techniques to manipulate 3'UTRs at endogenous gene loci will help to disentangle multi-functionality of proteins. Furthermore, the ability to directly examine the impact of gene regulatory elements in 3'UTRs will provide further insights into their functional significance.
Collapse
Affiliation(s)
- Sibylle Mitschka
- Cancer Biology and Genetics Program, Memorial Sloan Kettering Cancer Center, New York, NY, United States
| | - Mervin M Fansler
- Cancer Biology and Genetics Program, Memorial Sloan Kettering Cancer Center, New York, NY, United States; Tri-Institutional Training Program in Computational Biology and Medicine, Weill-Cornell Graduate College, New York, NY, United States
| | - Christine Mayr
- Cancer Biology and Genetics Program, Memorial Sloan Kettering Cancer Center, New York, NY, United States; Tri-Institutional Training Program in Computational Biology and Medicine, Weill-Cornell Graduate College, New York, NY, United States.
| |
Collapse
|
59
|
Buchumenski I, Holler K, Appelbaum L, Eisenberg E, Junker JP, Levanon EY. Systematic identification of A-to-I RNA editing in zebrafish development and adult organs. Nucleic Acids Res 2021; 49:4325-4337. [PMID: 33872356 PMCID: PMC8096273 DOI: 10.1093/nar/gkab247] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2020] [Revised: 03/05/2021] [Accepted: 04/15/2021] [Indexed: 11/18/2022] Open
Abstract
A-to-I RNA editing is a common post transcriptional mechanism, mediated by the Adenosine deaminase that acts on RNA (ADAR) enzymes, that increases transcript and protein diversity. The study of RNA editing is limited by the absence of editing maps for most model organisms, hindering the understanding of its impact on various physiological conditions. Here, we mapped the vertebrate developmental landscape of A-to-I RNA editing, and generated the first comprehensive atlas of editing sites in zebrafish. Tens of thousands unique editing events and 149 coding sites were identified with high-accuracy. Some of these edited sites are conserved between zebrafish and humans. Sequence analysis of RNA over seven developmental stages revealed high levels of editing activity in early stages of embryogenesis, when embryos rely on maternal mRNAs and proteins. In contrast to the other organisms studied so far, the highest levels of editing were detected in the zebrafish ovary and testes. This resource can serve as the basis for understanding of the role of editing during zebrafish development and maturity.
Collapse
Affiliation(s)
- Ilana Buchumenski
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel
| | - Karoline Holler
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Lior Appelbaum
- The Faculty of Life Sciences and the Multidisciplinary Brain Research Center, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Eli Eisenberg
- Raymond and Beverly Sackler School of Physics and Astronomy and Sagol School of Neuroscience, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Jan Philipp Junker
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Erez Y Levanon
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel
| |
Collapse
|
60
|
Gourain V, Armant O, Lübke L, Diotel N, Rastegar S, Strähle U. Multi-Dimensional Transcriptome Analysis Reveals Modulation of Cholesterol Metabolism as Highly Integrated Response to Brain Injury. Front Neurosci 2021; 15:671249. [PMID: 34054419 PMCID: PMC8162057 DOI: 10.3389/fnins.2021.671249] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2021] [Accepted: 04/16/2021] [Indexed: 12/14/2022] Open
Abstract
Zebrafish is an attractive model to investigate regeneration of the nervous system. Despite major progress in our understanding of the underlying processes, the transcriptomic changes are largely unknown. We carried out a computational analysis of the transcriptome of the regenerating telencephalon integrating changes in the expression of mRNAs, their splice variants and investigated the putative role of regulatory RNAs in the modulation of these transcriptional changes. Profound changes in the expression of genes and their splice variants engaged in many distinct processes were observed. Differential transcription and splicing are important processes in response to injury of the telencephalon. As exemplified by the coordinated regulation of the cholesterol synthesizing enzymes and transporters, the genome responded to injury of the telencephalon in a multi-tiered manner with distinct and interwoven changes in expression of enzymes, transporters and their regulatory molecules. This coordinated genomic response involved a decrease of the mRNA of the key transcription factor SREBF2, induction of microRNAs (miR-182, miR-155, miR-146, miR-31) targeting cholesterol genes, shifts in abundance of splice variants as well as regulation of long non-coding RNAs. Cholesterol metabolism appears to be switched from synthesis to relocation of cholesterol. Based on our in silico analyses, this switch involves complementary and synergistic inputs by different regulatory principles. Our studies suggest that adaptation of cholesterol metabolism is a key process involved in regeneration of the injured zebrafish brain.
Collapse
Affiliation(s)
- Victor Gourain
- Institute of Biological and Chemical Systems-Biological Information Processing (IBCS-BIP), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.,UMR 1064 Centre de Recherche en Transplantation en Immunologie, Nantes, France
| | - Olivier Armant
- Institute of Biological and Chemical Systems-Biological Information Processing (IBCS-BIP), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.,PSE-ENV/SRTE/LECO, Institut de Radioprotection et de Sûreté Nucléaire (IRSN), Cadarache, Saint-Paul-Lez-Durance, France
| | - Luisa Lübke
- Institute of Biological and Chemical Systems-Biological Information Processing (IBCS-BIP), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Nicolas Diotel
- Institute of Biological and Chemical Systems-Biological Information Processing (IBCS-BIP), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.,UMR 1188, Diabète Athérothrombose Thérapies Réunion Océan Indien CYROI, Saint-Denis, France
| | - Sepand Rastegar
- Institute of Biological and Chemical Systems-Biological Information Processing (IBCS-BIP), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Uwe Strähle
- Institute of Biological and Chemical Systems-Biological Information Processing (IBCS-BIP), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.,COS, University Heidelberg, Heidelberg, Germany
| |
Collapse
|
61
|
Nascimento-Oliveira JI, Fantinatti BEA, Wolf IR, Cardoso AL, Ramos E, Rieder N, de Oliveira R, Martins C. Differential expression of miRNAs in the presence of B chromosome in the cichlid fish Astatotilapia latifasciata. BMC Genomics 2021; 22:344. [PMID: 33980143 PMCID: PMC8117508 DOI: 10.1186/s12864-021-07651-w] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Accepted: 04/23/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND B chromosomes (Bs) are extra elements observed in diverse eukaryotes, including animals, plants and fungi. Although Bs were first identified a century ago and have been studied in hundreds of species, their biology is still enigmatic. Recent advances in omics and big data technologies are revolutionizing the B biology field. These advances allow analyses of DNA, RNA, proteins and the construction of interactive networks for understanding the B composition and behavior in the cell. Several genes have been detected on the B chromosomes, although the interaction of B sequences and the normal genome remains poorly understood. RESULTS We identified 727 miRNA precursors in the A. latifasciata genome, 66% which were novel predicted sequences that had not been identified before. We were able to report the A. latifasciata-specific miRNAs and common miRNAs identified in other fish species. For the samples carrying the B chromosome (B+), we identified 104 differentially expressed (DE) miRNAs that are down or upregulated compared to samples without B chromosome (B-) (p < 0.05). These miRNAs share common targets in the brain, muscle and gonads. These targets were used to construct a protein-protein-miRNA network showing the high interaction between the targets of differentially expressed miRNAs in the B+ chromosome samples. Among the DE-miRNA targets there are protein-coding genes reported for the B chromosome that are present in the protein-protein-miRNA network. Additionally, Gene Ontology (GO) terms related to nuclear matrix organization and response to stimulus are exclusive to DE miRNA targets of B+ samples. CONCLUSIONS This study is the first to report the connection of B chromosomes and miRNAs in a vertebrate species. We observed that the B chromosome impacts the miRNAs expression in several tissues and these miRNAs target several mRNAs involved with important biological processes.
Collapse
Affiliation(s)
- Jordana Inácio Nascimento-Oliveira
- Department of Structural and Functional Biology, Institute of Bioscience at Botucatu, São Paulo State University (UNESP), Botucatu, SP, 18618-689, Brazil
| | | | - Ivan Rodrigo Wolf
- Department of Structural and Functional Biology, Institute of Bioscience at Botucatu, São Paulo State University (UNESP), Botucatu, SP, 18618-689, Brazil
| | - Adauto Lima Cardoso
- Department of Structural and Functional Biology, Institute of Bioscience at Botucatu, São Paulo State University (UNESP), Botucatu, SP, 18618-689, Brazil
| | - Erica Ramos
- Department of Structural and Functional Biology, Institute of Bioscience at Botucatu, São Paulo State University (UNESP), Botucatu, SP, 18618-689, Brazil
| | - Nathalie Rieder
- Faculty of Mathematics and Natural Sciences, University of Bonn, Bonn, Germany
| | - Rogerio de Oliveira
- Department of Biostatistics, Plant Biology, Parasitology and Zoology, Institute of Bioscience at Botucatu, São Paulo State University (UNESP), Botucatu, SP, Brazil
| | - Cesar Martins
- Department of Structural and Functional Biology, Institute of Bioscience at Botucatu, São Paulo State University (UNESP), Botucatu, SP, 18618-689, Brazil.
| |
Collapse
|
62
|
Bhat P, Burkard TR, Herzog VA, Pauli A, Ameres SL. Systematic refinement of gene annotations by parsing mRNA 3' end sequencing datasets. Methods Enzymol 2021; 655:205-223. [PMID: 34183122 DOI: 10.1016/bs.mie.2021.03.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
Alternative cleavage and polyadenylation generates mRNA 3' isoforms in a cell type-specific manner. Due to finite available RNA sequencing data of organisms with vast cell type complexity, currently available gene annotation resources are incomplete, which poses significant challenges to the comprehensive interpretation and quantification of transcriptomes. In this chapter, we introduce 3'GAmES, a stand-alone computational pipeline for the identification and quantification of novel mRNA 3'end isoforms from 3'mRNA sequencing data. 3'GAmES expands available repositories and improves comprehensive gene-tag counting by cost-effective 3' mRNA sequencing, faithfully mirroring whole-transcriptome RNAseq measurements. By employing R and bash shell scripts (assembled in a Singularity container) 3'GAmES systematically augments cell type-specific 3' ends of RNA polymerase II transcripts and increases the sensitivity of quantitative gene expression profiling by 3' mRNA sequencing. Public access: https://github.com/AmeresLab/3-GAmES.git.
Collapse
Affiliation(s)
- Pooja Bhat
- Institute of Molecular Biotechnology (IMBA), Vienna BioCenter (VBC), Vienna, Austria; Vienna BioCenter PhD Program, Doctoral School of the University at Vienna and Medical University of Vienna, Vienna, Austria
| | - Thomas R Burkard
- Institute of Molecular Biotechnology (IMBA), Vienna BioCenter (VBC), Vienna, Austria
| | - Veronika A Herzog
- Institute of Molecular Biotechnology (IMBA), Vienna BioCenter (VBC), Vienna, Austria
| | - Andrea Pauli
- Research Institute of Molecular Pathology (IMP), Vienna BioCenter (VBC), Vienna, Austria
| | - Stefan L Ameres
- Institute of Molecular Biotechnology (IMBA), Vienna BioCenter (VBC), Vienna, Austria; Max Perutz Labs, University of Vienna, Vienna BioCenter (VBC), Vienna, Austria.
| |
Collapse
|
63
|
Computational analysis of alternative polyadenylation from standard RNA-seq and single-cell RNA-seq data. Methods Enzymol 2021; 655:225-243. [PMID: 34183123 DOI: 10.1016/bs.mie.2021.03.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Alternative polyadenylation (APA) is a major mechanism of post-transcriptional regulation in various cellular processes including cell proliferation and differentiation. Since conventional APA profiling methods have not been widely adopted, global APA studies are very limited. In this chapter, we summarize current computational methods for analyzing APA in standard RNA-seq and scRNA-seq data and describe two state-of-the-art bioinformatic algorithms DaPars and scDaPars in detail. The bioinformatic pipelines for both DaPars and scDaPars are presented and the application of both algorithms are highlighted.
Collapse
|
64
|
Wright K, de Silva K, Plain KM, Purdie AC, Blair TA, Duggin IG, Britton WJ, Oehlers SH. Mycobacterial infection-induced miR-206 inhibits protective neutrophil recruitment via the CXCL12/CXCR4 signalling axis. PLoS Pathog 2021; 17:e1009186. [PMID: 33826679 PMCID: PMC8055004 DOI: 10.1371/journal.ppat.1009186] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Revised: 04/19/2021] [Accepted: 03/29/2021] [Indexed: 12/22/2022] Open
Abstract
Pathogenic mycobacteria actively dysregulate protective host immune signalling pathways during infection to drive the formation of permissive granuloma microenvironments. Dynamic regulation of host microRNA (miRNA) expression is a conserved feature of mycobacterial infections across host-pathogen pairings. Here we examine the role of miR-206 in the zebrafish model of Mycobacterium marinum infection, which allows investigation of the early stages of granuloma formation. We find miR-206 is upregulated following infection by pathogenic M. marinum and that antagomir-mediated knockdown of miR-206 is protective against infection. We observed striking upregulation of cxcl12a and cxcr4b in infected miR-206 knockdown zebrafish embryos and live imaging revealed enhanced recruitment of neutrophils to sites of infection. We used CRISPR/Cas9-mediated knockdown of cxcl12a and cxcr4b expression and AMD3100 inhibition of Cxcr4 to show that the enhanced neutrophil response and reduced bacterial burden caused by miR-206 knockdown was dependent on the Cxcl12/Cxcr4 signalling axis. Together, our data illustrate a pathway through which pathogenic mycobacteria induce host miR-206 expression to suppress Cxcl12/Cxcr4 signalling and prevent protective neutrophil recruitment to granulomas.
Collapse
Affiliation(s)
- Kathryn Wright
- Tuberculosis Research Program at the Centenary Institute, The University of Sydney, Camperdown, New South Wales, Australia
- The University of Sydney, Faculty of Science, Sydney School of Veterinary Science, Sydney, New South Wales, Australia
| | - Kumudika de Silva
- The University of Sydney, Faculty of Science, Sydney School of Veterinary Science, Sydney, New South Wales, Australia
| | - Karren M. Plain
- The University of Sydney, Faculty of Science, Sydney School of Veterinary Science, Sydney, New South Wales, Australia
| | - Auriol C. Purdie
- The University of Sydney, Faculty of Science, Sydney School of Veterinary Science, Sydney, New South Wales, Australia
| | - Tamika A. Blair
- ithree Institute, University of Technology Sydney, Ultimo, New South Wales, Australia
| | - Iain G. Duggin
- ithree Institute, University of Technology Sydney, Ultimo, New South Wales, Australia
| | - Warwick J. Britton
- Tuberculosis Research Program at the Centenary Institute, The University of Sydney, Camperdown, New South Wales, Australia
- Department of Clinical Immunology, Royal Prince Alfred Hospital, Camperdown, New South Wales, Australia
| | - Stefan H. Oehlers
- Tuberculosis Research Program at the Centenary Institute, The University of Sydney, Camperdown, New South Wales, Australia
- The University of Sydney, Faculty of Medicine and Health & Marie Bashir Institute, Camperdown, New South Wales, Australia
| |
Collapse
|
65
|
Lee S, Wei L, Zhang B, Goering R, Majumdar S, Wen J, Taliaferro JM, Lai EC. ELAV/Hu RNA binding proteins determine multiple programs of neural alternative splicing. PLoS Genet 2021; 17:e1009439. [PMID: 33826609 PMCID: PMC8055025 DOI: 10.1371/journal.pgen.1009439] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Revised: 04/19/2021] [Accepted: 02/22/2021] [Indexed: 12/15/2022] Open
Abstract
ELAV/Hu factors are conserved RNA binding proteins (RBPs) that play diverse roles in mRNA processing and regulation. The founding member, Drosophila Elav, was recognized as a vital neural factor 35 years ago. Nevertheless, little was known about its impacts on the transcriptome, and potential functional overlap with its paralogs. Building on our recent findings that neural-specific lengthened 3' UTR isoforms are co-determined by ELAV/Hu factors, we address their impacts on splicing. While only a few splicing targets of Drosophila are known, ectopic expression of each of the three family members (Elav, Fne and Rbp9) alters hundreds of cassette exon and alternative last exon (ALE) splicing choices. Reciprocally, double mutants of elav/fne, but not elav alone, exhibit opposite effects on both classes of regulated mRNA processing events in larval CNS. While manipulation of Drosophila ELAV/Hu RBPs induces both exon skipping and inclusion, characteristic ELAV/Hu motifs are enriched only within introns flanking exons that are suppressed by ELAV/Hu factors. Moreover, the roles of ELAV/Hu factors in global promotion of distal ALE splicing are mechanistically linked to terminal 3' UTR extensions in neurons, since both processes involve bypass of proximal polyadenylation signals linked to ELAV/Hu motifs downstream of cleavage sites. We corroborate the direct action of Elav in diverse modes of mRNA processing using RRM-dependent Elav-CLIP data from S2 cells. Finally, we provide evidence for conservation in mammalian neurons, which undergo broad programs of distal ALE and APA lengthening, linked to ELAV/Hu motifs downstream of regulated polyadenylation sites. Overall, ELAV/Hu RBPs orchestrate multiple broad programs of neuronal mRNA processing and isoform diversification in Drosophila and mammalian neurons.
Collapse
Affiliation(s)
- Seungjae Lee
- Developmental Biology Program, Sloan Kettering Institute, New York City, New York, United States of America
| | - Lu Wei
- Developmental Biology Program, Sloan Kettering Institute, New York City, New York, United States of America
| | - Binglong Zhang
- Developmental Biology Program, Sloan Kettering Institute, New York City, New York, United States of America
| | - Raeann Goering
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, Colorado, United States of America
- RNA Bioscience Initiative University of Colorado Anschutz Medical Campus, Aurora, Colorado, United States of America
| | - Sonali Majumdar
- Developmental Biology Program, Sloan Kettering Institute, New York City, New York, United States of America
| | - Jiayu Wen
- Department of Genome Sciences, The John Curtin School of Medical Research, The Australian National University, Canberra, Australia
| | - J. Matthew Taliaferro
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, Colorado, United States of America
- RNA Bioscience Initiative University of Colorado Anschutz Medical Campus, Aurora, Colorado, United States of America
| | - Eric C. Lai
- Developmental Biology Program, Sloan Kettering Institute, New York City, New York, United States of America
| |
Collapse
|
66
|
Function analysis and molecular characterization of cyclin A in ovary development of oriental river prawn, Macrobrachium nipponense. Gene 2021; 788:145583. [PMID: 33753150 DOI: 10.1016/j.gene.2021.145583] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Revised: 11/26/2020] [Accepted: 03/12/2021] [Indexed: 11/24/2022]
Abstract
Macrobrachium nipponense has the characteristics of fast ovarian development cycle, which leads to the coexistence of multiple generations, the reduction of commodity specifications and the low economic benefit. Therefore, the study on the mechanism of ovarian development is of great significance to the development of industry. Cyclin A (CycA)is a key gene regulating ovarian development in vertebrates, but little information was available for its function in crustaceans. In this study, the full-length cDNA of Mn-CycA was obtained from the ovary. The full-length cDNA (2033 bp) with an open reading frame of 1368 bp, encoded a 456-amino acid protein. qRT-PCR revealed tissue-specific expression pattern of Mn-CycA, with abundant expression in the ovary. Results in different developmental stages of ovary indicated that Mn-CycA expression is positively correlated with ovarian maturation. qRT-PCR In different developmental stages, the expression of Mn-CycA mRNA gradually increased during the embryonic stage and decreased significantly on the first day of the hatching stage. At the 25th day of the metamorphosis stage, the expression level of Mn-CycAmRNA in female shrimp was 3.5 times higher than that in male shrimp, which may be related to the proliferation of oogonia and the formation of oocytes. In situ hybridization (ISH) of ovary showed Mn-CycA was examined in all stages and was mainly located in oogonia and oocytes. Compared with the control group, the obvious change of gonad somatic index (GSI) proved that injection of Mn-CycA dsRNA could delay the ovarian development cycle, which provided strong evidence for the involvement of Mn-CycA in ovarian maturation and oogenesis, and expanded a new perspective for studying the fast ovarian development cycle in M. nipponense.
Collapse
|
67
|
Luu I, Ikert H, Craig PM. Chronic exposure to anthropogenic and climate related stressors alters transcriptional responses in the liver of zebrafish (Danio rerio) across multiple generations. Comp Biochem Physiol C Toxicol Pharmacol 2021; 240:108918. [PMID: 33141083 DOI: 10.1016/j.cbpc.2020.108918] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/16/2020] [Revised: 10/09/2020] [Accepted: 10/18/2020] [Indexed: 01/04/2023]
Abstract
The antidepressant, venlafaxine (VFX), and climate change stressors, such as increased water temperature and decreased dissolved oxygen, are current threats to aquatic environments. This study aimed to determine how microRNAs (miRNAs) and predicted targeted transcripts were altered in livers of zebrafish exposed to these stressors, and livers of their un-exposed F1 and F2 offspring. Following a 21 day exposure to multiple stressors (1 μg/L VFX, +5 °C ambient, 50% O2), then a subsequent 21 day recovery, relative abundances of cyp3a65, hsp70, hsp90, and ppargc1a and miRNAs predicted to target them (miR-142a, miR-16c, miR-181c, and miR-129, respectively) were measured in the liver via quantitative PCR (RT-qPCR). There were significant decreases in miR-142a in the exposed F0 generation and the exposed F1 generation. While there were no changes detected in cyp3a65 relative abundance, there was a significant inverse relationship between cyp3a65 and miR-142a. Hsp70 expression significantly increased in the F1 generation, which persisted to the F2 generation and the relative abundance of hsp90 significantly increased in all generations. There was a significant reduction in miR-181c in the F1 generation, but there was no significant relationship between miR-181c and hsp90. Finally, there was a significant decrease in ppargc1a relative abundance in the F1 generation which was associated with an increase in miR-129. Combined, these results suggest that parental exposure to multiple, environmentally relevant stressors can confer transcriptional and epigenetic responses in the F1 and F2 generations, although identifying which stressor is a driving force becomes unclear.
Collapse
Affiliation(s)
- Ivy Luu
- Department of Biology, University of Waterloo, 200 University Ave. W., Waterloo N2L 3G1, Ontario, Canada
| | - Heather Ikert
- Department of Biology, University of Waterloo, 200 University Ave. W., Waterloo N2L 3G1, Ontario, Canada.
| | - Paul M Craig
- Department of Biology, University of Waterloo, 200 University Ave. W., Waterloo N2L 3G1, Ontario, Canada
| |
Collapse
|
68
|
Jin W, Zhu Q, Yang Y, Yang W, Wang D, Yang J, Niu X, Yu D, Gong J. Animal-APAdb: a comprehensive animal alternative polyadenylation database. Nucleic Acids Res 2021; 49:D47-D54. [PMID: 32986825 PMCID: PMC7779049 DOI: 10.1093/nar/gkaa778] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Revised: 08/27/2020] [Accepted: 09/08/2020] [Indexed: 12/31/2022] Open
Abstract
Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that recognizes different polyadenylation signals on transcripts, resulting in transcripts with different lengths of 3′ untranslated regions and thereby influencing a series of biological processes. Recent studies have highlighted the important roles of APA in human. However, APA profiles in other animals have not been fully recognized, and there is no database that provides comprehensive APA information for other animals except human. Here, by using the RNA sequencing data collected from public databases, we systematically characterized the APA profiles in 9244 samples of 18 species. In total, we identified 342 952 APA events with a median of 17 020 per species using the DaPars2 algorithm, and 315 691 APA events with a median of 17 953 per species using the QAPA algorithm in these 18 species, respectively. In addition, we predicted the polyadenylation sites (PAS) and motifs near PAS of these species. We further developed Animal-APAdb, a user-friendly database (http://gong_lab.hzau.edu.cn/Animal-APAdb/) for data searching, browsing and downloading. With comprehensive information of APA events in different tissues of different species, Animal-APAdb may greatly facilitate the exploration of animal APA patterns and novel mechanisms, gene expression regulation and APA evolution across tissues and species.
Collapse
Affiliation(s)
- Weiwei Jin
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China
| | - Qizhao Zhu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, P.R. China
| | - Yanbo Yang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China
| | - Wenqian Yang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China
| | - Dongyang Wang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China
| | - Jiajun Yang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China
| | - Xiaohui Niu
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China
| | - Debing Yu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, P.R. China
| | - Jing Gong
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, P.R. China.,College of Biomedicine and Health, Huazhong Agricultural University, Wuhan 430070, P.R. China
| |
Collapse
|
69
|
Zygotic Genome Activation: Critical Prelude to the Most Important Time of Your Life. Methods Mol Biol 2021; 2218:319-329. [PMID: 33606242 DOI: 10.1007/978-1-0716-0970-5_25] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Activation of the embryonic genome during development represents a major developmental transition in all species. The history of its exploration began in the 1950s-1960s, when this idea was put forward and proven experimentally by Alexander Neyfakh. He observed the aberrant development of fish embryos upon X-ray irradiation and noted the different developmental outcomes depending on the stage when fertilized eggs were subjected to irradiation. Neyfakh also discriminated a regional difference of X-irradiation between the nucleus and the cytoplasm. By selecting the X-ray dose causing nuclear damage, he determined the beginning of zygotic transcription, which at that time became known as the morphogenetic function of nuclei. His team defined the link of zygotic transcription with the asynchronization of cell division and cell migration, the two other hallmarks, which along with the morphogenetic function (or the zygotic genome activation), are at the core of the mid-blastula transition during development. Within this framework, current studies using maternal mutants and application of modern methods of whole-embryo and single-cell transcriptomics begin to decipher the molecular mechanisms of the mid-blastula transition (or the maternal-zygotic transition).
Collapse
|
70
|
Kontur C, Jeong M, Cifuentes D, Giraldez AJ. Ythdf m 6A Readers Function Redundantly during Zebrafish Development. Cell Rep 2020; 33:108598. [PMID: 33378672 PMCID: PMC11407899 DOI: 10.1016/j.celrep.2020.108598] [Citation(s) in RCA: 66] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Revised: 11/09/2020] [Accepted: 12/14/2020] [Indexed: 12/19/2022] Open
Abstract
During the maternal-to-zygotic transition (MZT), multiple mechanisms precisely control massive decay of maternal mRNAs. N6-methyladenosine (m6A) is known to regulate mRNA decay, yet how this modification promotes maternal transcript degradation remains unclear. Here, we find that m6A promotes maternal mRNA deadenylation. Yet, genetic loss of m6A readers Ythdf2 and Ythdf3 did not impact global maternal mRNA clearance, zygotic genome activation, or the onset of gastrulation, challenging the view that Ythdf2 alone is critical to developmental timing. We reveal that Ythdf proteins function redundantly during zebrafish oogenesis and development, as double Ythdf2 and Ythdf3 deletion prevented female gonad formation and triple Ythdf mutants were lethal. Finally, we show that the microRNA miR-430 functions additively with methylation to promote degradation of common transcript targets. Together these findings reveal that m6A facilitates maternal mRNA deadenylation and that multiple pathways and readers act in concert to mediate these effects of methylation on RNA stability.
Collapse
Affiliation(s)
- Cassandra Kontur
- Department of Genetics, Yale University School of Medicine, New Haven, CT 06510, USA.
| | - Minsun Jeong
- Chey Institute for Advanced Studies, Seoul 06141, Republic of Korea
| | - Daniel Cifuentes
- Department of Biochemistry, Boston University School of Medicine, Boston, MA 02118, USA
| | - Antonio J Giraldez
- Department of Genetics, Yale University School of Medicine, New Haven, CT 06510, USA; Yale Stem Cell Center, Yale University School of Medicine, New Haven, CT 06510, USA; Yale Cancer Center, Yale University School of Medicine, New Haven, CT 06510, USA.
| |
Collapse
|
71
|
Dicks S, Jürgensen L, Leuschner F, Hassel D, Andrieux G, Boerries M. Cardiac Regeneration and Tumor Growth-What Do They Have in Common? Front Genet 2020; 11:586658. [PMID: 33362851 PMCID: PMC7756072 DOI: 10.3389/fgene.2020.586658] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 11/04/2020] [Indexed: 11/13/2022] Open
Abstract
Acute myocardial infarction is a leading cause of death. Unlike most adult mammals, zebrafish have the capability to almost fully regenerate their hearts after injury. In contrast, ischemic damage in adult human and mouse hearts usually results in scar tissue. mRNA-Sequencing (Seq) and miRNA-Seq analyses of heart regeneration in zebrafish over time showed that the process can be divided into three phases: the first phase represents dedifferentiation and proliferation of cells, the second phase is characterized by migration, and in the third phase cell signals indicate heart development and differentiation. The first two phases seem to share major similarities with tumor development and growth. To gain more insight into these similarities between cardiac regeneration and tumor development and growth, we used patient matched tumor normal (“healthy”) RNA-Seq data for several tumor entities from The Cancer Genome Atlas (TCGA). Subsequently, RNA data were processed using the same pipeline for both the zebrafish samples and tumor datasets. Functional analysis showed that multiple Gene Ontology terms (GO terms) are involved in both early stage cardiac regeneration and tumor development/growth across multiple tumor entities. These GO terms are mostly associated with cell cycle processes. Further analysis showed that orthologous genes are the same key players that regulated these changes in both diseases. We also observed that GO terms associated with heart development in the third late phase of cardiac regeneration are downregulated in the tumor entities. Taken together, our analysis illustrates similarities between cardiac remodeling and tumor progression.
Collapse
Affiliation(s)
- Severin Dicks
- Faculty of Medicine, Institute of Medical Bioinformatics and Systems Medicine, Medical Center-University of Freiburg, Albert-Ludwigs-Universität, Freiburg, Germany.,Faculty of Biology, Albert-Ludwigs-Universität, Freiburg, Germany
| | - Lonny Jürgensen
- Department of Cardiology, Angiology and Pneumology, University of Heidelberg, Heidelberg, Germany.,DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Florian Leuschner
- Department of Cardiology, Angiology and Pneumology, University of Heidelberg, Heidelberg, Germany.,DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - David Hassel
- Department of Cardiology, Angiology and Pneumology, University of Heidelberg, Heidelberg, Germany.,DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Geoffroy Andrieux
- Faculty of Medicine, Institute of Medical Bioinformatics and Systems Medicine, Medical Center-University of Freiburg, Albert-Ludwigs-Universität, Freiburg, Germany.,German Cancer Consortium (DKTK), German Cancer Research Center (DKFZ), Freiburg, Germany
| | - Melanie Boerries
- Faculty of Medicine, Institute of Medical Bioinformatics and Systems Medicine, Medical Center-University of Freiburg, Albert-Ludwigs-Universität, Freiburg, Germany.,German Cancer Consortium (DKTK), German Cancer Research Center (DKFZ), Freiburg, Germany.,Comprehensive Cancer Center Freiburg (CCCF), Medical Center-University of Freiburg, Freiburg, Germany
| |
Collapse
|
72
|
Alternative Polyadenylation: a new frontier in post transcriptional regulation. Biomark Res 2020; 8:67. [PMID: 33292571 PMCID: PMC7690165 DOI: 10.1186/s40364-020-00249-6] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 11/16/2020] [Indexed: 12/13/2022] Open
Abstract
Polyadenylation of pre-messenger RNA (pre-mRNA) specific sites and termination of their downstream transcriptions are signaled by unique sequence motif structures such as AAUAAA and its auxiliary elements. Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that processes RNA products depending on its 3'-untranslated region (3'-UTR) specific sequence signal. APA processing can generate several mRNA isoforms from a single gene, which may have different biological functions on their target gene. As a result, cellular genomic stability, proliferation capability, and transformation feasibility could all be affected. Furthermore, APA modulation regulates disease initiation and progression. APA status could potentially act as a biomarker for disease diagnosis, severity stratification, and prognosis forecast. While the advance of modern throughout technologies, such as next generation-sequencing (NGS) and single-cell sequencing techniques, have enriched our knowledge about APA, much of APA biological process is unknown and pending for further investigation. Herein, we review the current knowledge on APA and how its regulatory complex factors (CFI/IIm, CPSF, CSTF, and RBPs) work together to determine RNA splicing location, cell cycle velocity, microRNA processing, and oncogenesis regulation. We also discuss various APA experiment strategies and the future direction of APA research.
Collapse
|
73
|
Liyanage TD, Nikapitiya C, Lee J, De Zoysa M. Molecular insight into regulation of miRNAs in the spleen of zebrafish (Danio rerio) upon pathogenic Streptococcus parauberis infection. FISH & SHELLFISH IMMUNOLOGY 2020; 106:898-909. [PMID: 32889099 DOI: 10.1016/j.fsi.2020.08.045] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2020] [Revised: 07/10/2020] [Accepted: 08/27/2020] [Indexed: 06/11/2023]
Abstract
MicroRNAs (miRNAs) constitute a group of small non-coding RNAs (~22 nucleotides) and one of their main functions is to regulate the immune responses. Gram-positive bacterium, Streptococcus parauberis is the main causative agent of "Streptococcosis" in wide range of fish species. In this study, we performed high throughput sequencing analysis to identify the miRNA profile against S. parauberis infection in the spleen of zebrafish (Danio rerio). Overall, 349 known and 151 novel miRNAs were discovered. Among them, 12 known miRNAs (dre-miR-34b, dre-miR-135a, dre-miR-200b-5p, dre-miR-146b, dre-miR-31, dre-miR-17a-3p, dre-miR-222a-3p, dre-miR-731, dre-miR-301b-3p and dre-miR-30a-3p) and 9 novel miRNAs were differentially expressed (DE) in the spleen of S. parauberis challenged zebrafish. The identified 12 DE miRNAs were predicted to regulate 721 target genes. We confirmed the miRNA expression results by validating selected known and novel DE miRNAs using qRT-PCR. Gene Ontology (GO), Kyoto Encyclopedia of Genes (KEGG) pathway analysis and miRNA-mRNA interactions implies that specific target genes of DE miRNAs are associated with immune responses. The enriched pathways included Toll-like receptor (TLR), C-type lectin, NOD-like receptor, and RIG-I-like receptor signaling pathways, etc. Especially, dre-miR-200b-5p, dre-miR-146b, dre-miR-731, dre-miR-222a-3p, and dre-miR-34b were able to target potential immune-related genes such as il10, irak1, traf6, hspa8 and ikbke upon S. parauberis challenge. Thus, overall results could lay a foundation to understand the underlying immune regulatory role of miRNAs in response to pathogenic S. parauberis infection in teleosts.
Collapse
Affiliation(s)
- T D Liyanage
- College of Veterinary Medicine, Chungnam National University, Yuseong-gu Daejeon, 34134, Republic of Korea
| | - Chamilani Nikapitiya
- College of Veterinary Medicine, Chungnam National University, Yuseong-gu Daejeon, 34134, Republic of Korea
| | - Jehee Lee
- Department of Marine Life Sciences, Jeju National University, Jeju Self-Governing Province 63243, Republic of Korea
| | - Mahanama De Zoysa
- College of Veterinary Medicine, Chungnam National University, Yuseong-gu Daejeon, 34134, Republic of Korea.
| |
Collapse
|
74
|
Wei L, Lee S, Majumdar S, Zhang B, Sanfilippo P, Joseph B, Miura P, Soller M, Lai EC. Overlapping Activities of ELAV/Hu Family RNA Binding Proteins Specify the Extended Neuronal 3' UTR Landscape in Drosophila. Mol Cell 2020; 80:140-155.e6. [PMID: 33007254 DOI: 10.1016/j.molcel.2020.09.007] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 07/31/2020] [Accepted: 09/04/2020] [Indexed: 12/21/2022]
Abstract
The tissue-specific deployment of highly extended neural 3' UTR isoforms, generated by alternative polyadenylation (APA), is a broad and conserved feature of metazoan genomes. However, the factors and mechanisms that control neural APA isoforms are not well understood. Here, we show that three ELAV/Hu RNA binding proteins (Elav, Rbp9, and Fne) have similar capacities to induce a lengthened 3' UTR landscape in an ectopic setting. These factors promote accumulation of chromatin-associated, 3' UTR-extended, nascent transcripts, through inhibition of proximal polyadenylation site (PAS) usage. Notably, Elav represses an unannotated splice isoform of fne, switching the normally cytoplasmic Fne toward the nucleus in elav mutants. We use genomic profiling to reveal strong and broad loss of neural APA in elav/fne double mutant CNS, the first genetic background to largely abrogate this distinct APA signature. Overall, we demonstrate how regulatory interplay and functionally overlapping activities of neural ELAV/Hu RBPs drives the neural APA landscape.
Collapse
Affiliation(s)
- Lu Wei
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA
| | - Seungjae Lee
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA
| | - Sonali Majumdar
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA
| | - Binglong Zhang
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA
| | - Piero Sanfilippo
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA; Louis V. Gerstner, Jr. Graduate School of Biomedical Sciences, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA
| | - Brian Joseph
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA; Louis V. Gerstner, Jr. Graduate School of Biomedical Sciences, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA
| | - Pedro Miura
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA; Department of Biology, University of Nevada, Reno, Reno, NV 89557, USA
| | - Matthias Soller
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, UK
| | - Eric C Lai
- Department of Developmental Biology, Sloan Kettering Institute, New York, NY 10065, USA.
| |
Collapse
|
75
|
Carrasco J, Rauer M, Hummel B, Grzejda D, Alfonso-Gonzalez C, Lee Y, Wang Q, Puchalska M, Mittler G, Hilgers V. ELAV and FNE Determine Neuronal Transcript Signatures through EXon-Activated Rescue. Mol Cell 2020; 80:156-163.e6. [PMID: 33007255 DOI: 10.1016/j.molcel.2020.09.011] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 07/03/2020] [Accepted: 08/12/2020] [Indexed: 12/22/2022]
Abstract
The production of alternative RNA variants contributes to the tissue-specific regulation of gene expression. In the animal nervous system, a systematic shift toward distal sites of transcription termination produces transcript signatures that are crucial for neuron development and function. Here, we report that, in Drosophila, the highly conserved protein ELAV globally regulates all sites of neuronal 3' end processing and directly binds to proximal polyadenylation sites of target mRNAs in vivo. We uncover an endogenous strategy of functional gene rescue that safeguards neuronal RNA signatures in an ELAV loss-of-function context. When not directly repressed by ELAV, the transcript encoding the ELAV paralog FNE acquires a mini-exon, generating a new protein able to translocate to the nucleus and rescue ELAV-mediated alternative polyadenylation and alternative splicing. We propose that exon-activated functional rescue is a more widespread mechanism that ensures robustness of processes regulated by a hierarchy, rather than redundancy, of effectors.
Collapse
Affiliation(s)
- Judit Carrasco
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany; Faculty of Biology, Albert Ludwig University, 79104 Freiburg, Germany; International Max Planck Research School for Molecular and Cellular Biology (IMPRS-MCB), 79108 Freiburg, Germany
| | - Michael Rauer
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany
| | - Barbara Hummel
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany
| | - Dominika Grzejda
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany; Faculty of Biology, Albert Ludwig University, 79104 Freiburg, Germany; International Max Planck Research School for Molecular and Cellular Biology (IMPRS-MCB), 79108 Freiburg, Germany
| | - Carlos Alfonso-Gonzalez
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany; Faculty of Biology, Albert Ludwig University, 79104 Freiburg, Germany; International Max Planck Research School for Immunology, Epigenetics and Metabolism (IMPRS-IEM), 79108 Freiburg, Germany
| | - Yeon Lee
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Qingqing Wang
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Monika Puchalska
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany
| | - Gerhard Mittler
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany
| | - Valérie Hilgers
- Max-Planck-Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany.
| |
Collapse
|
76
|
Bae B, Gruner HN, Lynch M, Feng T, So K, Oliver D, Mastick GS, Yan W, Pieraut S, Miura P. Elimination of Calm1 long 3'-UTR mRNA isoform by CRISPR-Cas9 gene editing impairs dorsal root ganglion development and hippocampal neuron activation in mice. RNA (NEW YORK, N.Y.) 2020; 26:1414-1430. [PMID: 32522888 PMCID: PMC7491327 DOI: 10.1261/rna.076430.120] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Accepted: 06/06/2020] [Indexed: 05/04/2023]
Abstract
The majority of mouse and human genes are subject to alternative cleavage and polyadenylation (APA), which most often leads to the expression of two or more alternative length 3' untranslated region (3'-UTR) mRNA isoforms. In neural tissues, there is enhanced expression of APA isoforms with longer 3'-UTRs on a global scale, but the physiological relevance of these alternative 3'-UTR isoforms is poorly understood. Calmodulin 1 (Calm1) is a key integrator of calcium signaling that generates short (Calm1-S) and long (Calm1-L) 3'-UTR mRNA isoforms via APA. We found Calm1-L expression to be largely restricted to neural tissues in mice including the dorsal root ganglion (DRG) and hippocampus, whereas Calm1-S was more broadly expressed. smFISH revealed that both Calm1-S and Calm1-L were subcellularly localized to neural processes of primary hippocampal neurons. In contrast, cultured DRG showed restriction of Calm1-L to soma. To investigate the in vivo functions of Calm1-L, we implemented a CRISPR-Cas9 gene editing strategy to delete a small region encompassing the Calm1 distal poly(A) site. This eliminated Calm1-L expression while maintaining expression of Calm1-S Mice lacking Calm1-L (Calm1ΔL/ΔL ) exhibited disorganized DRG migration in embryos, and reduced experience-induced neuronal activation in the adult hippocampus. These data indicate that Calm1-L plays functional roles in the central and peripheral nervous systems.
Collapse
Affiliation(s)
- Bongmin Bae
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Hannah N Gruner
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Maebh Lynch
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Ting Feng
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Kevin So
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Daniel Oliver
- Department of Physiology and Cell Biology, University of Nevada, Reno School of Medicine, Reno, Nevada 89557, USA
| | - Grant S Mastick
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Wei Yan
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
- Department of Physiology and Cell Biology, University of Nevada, Reno School of Medicine, Reno, Nevada 89557, USA
| | - Simon Pieraut
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| | - Pedro Miura
- Department of Biology, University of Nevada, Reno, Nevada 89557, USA
| |
Collapse
|
77
|
Levin M, Zalts H, Mostov N, Hashimshony T, Yanai I. Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms. Nucleic Acids Res 2020; 48:5926-5938. [PMID: 32421815 PMCID: PMC7293032 DOI: 10.1093/nar/gkaa359] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2020] [Revised: 04/11/2020] [Accepted: 04/27/2020] [Indexed: 01/08/2023] Open
Abstract
Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.
Collapse
Affiliation(s)
- Michal Levin
- Quantitative Proteomics, Institute of Molecular Biology, Mainz 55128, Germany
| | - Harel Zalts
- Faculty of Biology, Technion - Israel Institute of Technology, Haifa 3200003, Israel
| | - Natalia Mostov
- Faculty of Biology, Technion - Israel Institute of Technology, Haifa 3200003, Israel
| | - Tamar Hashimshony
- Faculty of Biology, Technion - Israel Institute of Technology, Haifa 3200003, Israel
| | - Itai Yanai
- Institute for Computational Medicine, NYU Grossman School of Medicine, New York 10016, USA
| |
Collapse
|
78
|
Lawson ND, Li R, Shin M, Grosse A, Yukselen O, Stone OA, Kucukural A, Zhu L. An improved zebrafish transcriptome annotation for sensitive and comprehensive detection of cell type-specific genes. eLife 2020; 9:55792. [PMID: 32831172 PMCID: PMC7486121 DOI: 10.7554/elife.55792] [Citation(s) in RCA: 75] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 08/21/2020] [Indexed: 02/07/2023] Open
Abstract
The zebrafish is ideal for studying embryogenesis and is increasingly applied to model human disease. In these contexts, RNA-sequencing (RNA-seq) provides mechanistic insights by identifying transcriptome changes between experimental conditions. Application of RNA-seq relies on accurate transcript annotation for a genome of interest. Here, we find discrepancies in analysis from RNA-seq datasets quantified using Ensembl and RefSeq zebrafish annotations. These issues were due, in part, to variably annotated 3' untranslated regions and thousands of gene models missing from each annotation. Since these discrepancies could compromise downstream analyses and biological reproducibility, we built a more comprehensive zebrafish transcriptome annotation that addresses these deficiencies. Our annotation improves detection of cell type-specific genes in both bulk and single cell RNA-seq datasets, where it also improves resolution of cell clustering. Thus, we demonstrate that our new transcriptome annotation can outperform existing annotations, providing an important resource for zebrafish researchers.
Collapse
Affiliation(s)
- Nathan D Lawson
- Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
| | - Rui Li
- Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
| | - Masahiro Shin
- Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
| | - Ann Grosse
- Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, United States
| | - Onur Yukselen
- Bioinformatics Core, University of Massachusetts Medical School, Worcester, United States
| | - Oliver A Stone
- Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford, United Kingdom
| | - Alper Kucukural
- Program in Molecular Medicine, University of Massachusetts Medical School, Worcester, United States.,Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, United States
| | - Lihua Zhu
- Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, United States.,Program in Molecular Medicine, University of Massachusetts Medical School, Worcester, United States.,Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, United States
| |
Collapse
|
79
|
Liu Y, Zhu Z, Ho IHT, Shi Y, Li J, Wang X, Chan MTV, Cheng CHK. Genetic Deletion of miR-430 Disrupts Maternal-Zygotic Transition and Embryonic Body Plan. Front Genet 2020; 11:853. [PMID: 32849832 PMCID: PMC7417628 DOI: 10.3389/fgene.2020.00853] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Accepted: 07/13/2020] [Indexed: 11/13/2022] Open
Abstract
MiR-430 is considered an important regulator during embryonic development, but genetic loss-of-function study is still lacking. Here we demonstrated that genetic deletion of the miR-430 cluster resulted in developmental defects in cell movement, germ layer specification, axis patterning and organ progenitor formation in zebrafish. Transcriptome analysis indicated that the maternally provided transcripts were not properly degraded whereas the zygotic genome expressed genes were not fully activated in the miR-430 mutants. We further found that a reciprocal regulatory loop exists between miR-430 and maternally provided transcripts: the maternally provided transcripts (Nanog, Dicer1, Dgcr8, and AGOs) are required for miR-430 biogenesis and function, whereas miR-430 is required for the clearance of these maternally provided transcripts. These data provide the first genetic evidence that miR-430 is required for maternal-zygotic transition and subsequent establishment of embryonic body plan.
Collapse
Affiliation(s)
- Yun Liu
- School of Biomedical Sciences, The Chinese University of Hong Kong, Hong Kong, China.,Shenzhen Zelong Biological Technology Limited Cooperation, Shenzhen, China
| | - Zeyao Zhu
- Department of Anatomy, Histology and Developmental Biology, School of Basic Medical Sciences, Shenzhen University, Shenzhen, China
| | - Idy H T Ho
- School of Biomedical Sciences, The Chinese University of Hong Kong, Hong Kong, China.,Department of Anesthesia and Intensive Care, Prince of Wales Hospital, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, China
| | - Yujian Shi
- School of Biomedical Sciences, The Chinese University of Hong Kong, Hong Kong, China
| | - Jianzhen Li
- College of Life Sciences, Northwest Normal University, Lanzhou, China
| | - Xia Wang
- Department of Anatomy, Histology and Developmental Biology, School of Basic Medical Sciences, Shenzhen University, Shenzhen, China
| | - Matthew T V Chan
- Department of Anesthesia and Intensive Care, Prince of Wales Hospital, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, China
| | - Christopher H K Cheng
- School of Biomedical Sciences, The Chinese University of Hong Kong, Hong Kong, China
| |
Collapse
|
80
|
Genome-Wide Analysis Reveals Changes in Long Noncoding RNAs in the Differentiation of Canine BMSCs into Insulin-Producing Cells. Int J Mol Sci 2020; 21:ijms21155549. [PMID: 32756402 PMCID: PMC7432238 DOI: 10.3390/ijms21155549] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 07/24/2020] [Accepted: 07/31/2020] [Indexed: 12/27/2022] Open
Abstract
Long noncoding RNAs (lncRNAs) have been extensively explored over the past decade, including mice and humans. However, their impact on the transdifferentiation of canine bone marrow mesenchymal stem cells (cBMSCs) into insulin-producing cells (IPCs) is largely unknown. In this study, we used a three-step induction procedure to induce cBMSCs into IPCs, and samples (two biological replicates each) were obtained after each step; the samples consisted of “BMSCs” (B), “stage 1” (S1), “stage 2” (S2), “stage 3” (S3), and “islets” (I). After sequencing, 15,091 lncRNAs were identified, and we screened 110, 41, 23, and 686 differentially expressed lncRNAs (padjusted < 0.05) in B vs. S1, S1 vs. S2, S2 vs. S3, and I vs. S3 pairwise comparisons, respectively. In lncRNA target prediction, there were 166,623 colocalized targets and 2,976,362 correlated targets. Gene Ontology (GO) analysis showed that binding represented the main molecular functions of both the cis- and trans-modes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the insulin signaling pathway, Rap1 signaling pathway, tight junctions, MAPK signaling pathway, and cell cycle were enriched for these relative genes. The expression of lncRNAs was verified using qRT-PCR. This study provides a lncRNA catalog for future research concerning the mechanism of the transdifferentiation of cBMSCs into IPCs.
Collapse
|
81
|
Zhang Z, So K, Peterson R, Bauer M, Ng H, Zhang Y, Kim JH, Kidd T, Miura P. Elav-Mediated Exon Skipping and Alternative Polyadenylation of the Dscam1 Gene Are Required for Axon Outgrowth. Cell Rep 2020; 27:3808-3817.e7. [PMID: 31242415 PMCID: PMC7092717 DOI: 10.1016/j.celrep.2019.05.083] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2018] [Revised: 04/18/2019] [Accepted: 05/22/2019] [Indexed: 12/19/2022] Open
Abstract
Many metazoan genes express alternative long 3′ UTR isoforms in the nervous system, but their functions remain largely unclear. In Drosophila melanogaster, the Dscam1 gene generates short and long (Dscam1-L) 3′ UTR isoforms because of alternative polyadenylation (APA). Here, we found that the RNA-binding protein Embryonic Lethal Abnormal Visual System (Elav) impacts Dscam1 biogenesis at two levels, including regulation of long 3′ UTR biogenesis and skipping of an upstream exon (exon 19). MinION long-read sequencing confirmed the connectivity of this alternative splicing event to the long 3′ UTR. Knockdown or CRISPR deletion of Dscam1-L impaired axon outgrowth in Drosophila. The Dscam1 long 3′ UTR was found to be required for correct Elav-mediated skipping of exon 19. Elav thus co-regulates APA and alternative splicing to generate specific Dscam1 transcripts that are essential for neural development. This coupling of APA to alternative splicing might represent a new class of regulated RNA processing. Like most metazoan genes, Dscam1 expresses alternative short and long 3′ UTR mRNAs. Zhang et al. find that loss of Dscam1 long 3′ UTR transcripts impairs axon outgrowth in Drosophila. Long-read sequencing reveals that these long 3′ UTR mRNAs preferentially skip an upstream exon, altering Dscam1 amino acid sequence.
Collapse
Affiliation(s)
- Zhiping Zhang
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Kevin So
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Ryan Peterson
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Matthew Bauer
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Henry Ng
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Yong Zhang
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Jung Hwan Kim
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Thomas Kidd
- Department of Biology, University of Nevada, Reno, Reno, NV, USA
| | - Pedro Miura
- Department of Biology, University of Nevada, Reno, Reno, NV, USA.
| |
Collapse
|
82
|
Ren C, Deng K, Wang Z, Deng M, Fan Y, Zhang Y, Ma J, Wang S, Liu Z, Wang F. Reinterpreting sheep muscle strand-specific RNA sequencing data showing extensive 3'UTR extensions. Anim Genet 2020; 51:788-798. [PMID: 32696483 DOI: 10.1111/age.12987] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 05/16/2020] [Accepted: 06/23/2020] [Indexed: 12/01/2022]
Abstract
The more complex 3' UTR in higher organisms may have the function of increasing post-transcriptional gene regulation. Recent RNA sequencing technologies have provided us with the possibility to capture the complete 3' UTR landscape of different species and cells. However, no systematic analysis of sheep-related 3' UTR has been performed. Here, we conducted a detailed analysis of the 3' UTR with the primary goal of identifying intact 3' UTR landscapes in the sheep muscles of the three developmental stages. Based on strand-specific RNA sequencing (ssRNA-seq) data, we found that thousands of genes in sheep muscle are continuously transcribed after the UTR of the reference genome (Oar_v4.0). More than 66% of the 3' UTR extensions exhibit similar expression trends to their upstream gene exons. These 3' UTR extensions strongly enrich thousands of conserved microRNA binding sites. The 3' UTR extension-associated RNA of PFKM (PuaRNA) was predicted to be derived from the 3' UTR of PFKM. In sheep myocytes, myotubes and various tissues, the expression pattern of PuaRNA is positively correlated with that of PFKM. Taken together, these new 3' UTR annotations greatly extend the range of mammalian post-transcriptional regulatory networks, which have a particular impact on the regulation of sheep muscle development.
Collapse
Affiliation(s)
- Caifang Ren
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China.,School of Medicine, Jiangsu University, Zhengjiang, Jiangsu, 212013, China
| | - Kaiping Deng
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Zhibo Wang
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Mingtian Deng
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Yixuan Fan
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Yanli Zhang
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Jianyu Ma
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Shuting Wang
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Zifei Liu
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| | - Feng Wang
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Weigang Street, Xuanwu District, Nanjing, Jiangsu, 210095, China
| |
Collapse
|
83
|
miRTil: An Extensive Repository for Nile Tilapia microRNA Next Generation Sequencing Data. Cells 2020; 9:cells9081752. [PMID: 32707870 PMCID: PMC7465656 DOI: 10.3390/cells9081752] [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: 06/25/2020] [Revised: 07/16/2020] [Accepted: 07/20/2020] [Indexed: 12/04/2022] Open
Abstract
Nile tilapia is the third most cultivated fish worldwide and a novel model species for evolutionary studies. Aiming to improve productivity and contribute to the selection of traits of economic impact, biotechnological approaches have been intensively applied to species enhancement. In this sense, recent studies have focused on the multiple roles played by microRNAs (miRNAs) in the post-transcriptional regulation of protein-coding genes involved in the emergence of phenotypes with relevance for aquaculture. However, there is still a growing demand for a reference resource dedicated to integrating Nile Tilapia miRNA information, obtained from both experimental and in silico approaches, and facilitating the analysis and interpretation of RNA sequencing data. Here, we present an open repository dedicated to Nile Tilapia miRNAs: the “miRTil database”. The database stores data on 734 mature miRNAs identified in 11 distinct tissues and five key developmental stages. The database provides detailed information about miRNA structure, genomic context, predicted targets, expression profiles, and relative 5p/3p arm usage. Additionally, miRTil also includes a comprehensive pre-computed miRNA-target interaction network containing 4936 targets and 19,580 interactions.
Collapse
|
84
|
Nourse J, Spada S, Danckwardt S. Emerging Roles of RNA 3'-end Cleavage and Polyadenylation in Pathogenesis, Diagnosis and Therapy of Human Disorders. Biomolecules 2020; 10:biom10060915. [PMID: 32560344 PMCID: PMC7356254 DOI: 10.3390/biom10060915] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2020] [Revised: 06/10/2020] [Accepted: 06/13/2020] [Indexed: 12/11/2022] Open
Abstract
A crucial feature of gene expression involves RNA processing to produce 3′ ends through a process termed 3′ end cleavage and polyadenylation (CPA). This ensures the nascent RNA molecule can exit the nucleus and be translated to ultimately give rise to a protein which can execute a function. Further, alternative polyadenylation (APA) can produce distinct transcript isoforms, profoundly expanding the complexity of the transcriptome. CPA is carried out by multi-component protein complexes interacting with multiple RNA motifs and is tightly coupled to transcription, other steps of RNA processing, and even epigenetic modifications. CPA and APA contribute to the maintenance of a multitude of diverse physiological processes. It is therefore not surprising that disruptions of CPA and APA can lead to devastating disorders. Here, we review potential CPA and APA mechanisms involving both loss and gain of function that can have tremendous impacts on health and disease. Ultimately we highlight the emerging diagnostic and therapeutic potential CPA and APA offer.
Collapse
Affiliation(s)
- Jamie Nourse
- Institute for Clinical Chemistry and Laboratory Medicine, University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany; (J.N.); (S.S.)
- Center for Thrombosis and Hemostasis (CTH), University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany
| | - Stefano Spada
- Institute for Clinical Chemistry and Laboratory Medicine, University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany; (J.N.); (S.S.)
- Center for Thrombosis and Hemostasis (CTH), University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany
| | - Sven Danckwardt
- Institute for Clinical Chemistry and Laboratory Medicine, University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany; (J.N.); (S.S.)
- Center for Thrombosis and Hemostasis (CTH), University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany
- German Center for Cardiovascular Research (DZHK), Rhine-Main, Germany
- Correspondence:
| |
Collapse
|
85
|
Yu F, Zhang Y, Cheng C, Wang W, Zhou Z, Rang W, Yu H, Wei Y, Wu Q, Zhang Y. Poly(A)-seq: A method for direct sequencing and analysis of the transcriptomic poly(A)-tails. PLoS One 2020; 15:e0234696. [PMID: 32544193 PMCID: PMC7297374 DOI: 10.1371/journal.pone.0234696] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Accepted: 06/02/2020] [Indexed: 12/19/2022] Open
Abstract
Poly(A) tails at the 3' end of eukaryotic messenger RNAs control mRNA stability and translation efficiency. Facilitated by various NGS methods, alternative polyadenylation sites determining the 3'-UTR length of gene transcripts have been extensively studied. However, poly(A) lengths demonstrating dynamic and developmental regulation remain largely unexplored. The recently developed NGS-based methods for genome-wide poly(A) profiling have promoted the study of genom-wide poly(A) dynamics. Here we present a straight forward NGS-method for poly(A) profiling, which applies a direct 3'-end adaptor ligation and the template switching for 5'-end adaptor ligation for cDNA library construction. Poly(A) lengths are directly calculated from base call data using a self-developed pipeline pA-finder. The libraries were directly sequenced from the 3'-UTR regions into the followed poly(A) tails, firstly on NextSeq 500 to produce single-end 300-nt reads, demonstrating the method feasibility and that optimization of the fragmented RNA size for cDNA library construction could detecting longer poly (A) tails. We next applied Poly(A)-seq cDNA libraries containing 40-nt and 120-nt poly(A) tail spike-in RNAs on HiSeq X-ten and NovaSeq 6000 to obtain 150-nt and 250-nt pair-end reads. The sequencing profiles of the spike-in RNAs demonstrated both high accuracy and high quality score in reading poly(A) tails. The poly(A) signal bleeding into the 3' adaptor sequence and a sharp decreased quality score at the junction were observed, allowing the modification of pA-finder to remove homopolymeric signal bleeding. We hope that wide applications of Poly(A)-seq help facilitate the study of the development- and disease-related poly(A) dynamics and regulation, and of the recent emerging mixed tailing regulation.
Collapse
Affiliation(s)
- Fengyun Yu
- Laboratory for Genomics Regulation and Human Health, ABLife Inc., Wuhan, PR China
- ABLife BioBigData Instibute, Wuhan, PR China
| | - Yu Zhang
- Center for Genomics Analysis, ABLife Inc., Wuhan, PR China
| | - Chao Cheng
- ABLife BioBigData Instibute, Wuhan, PR China
- Center for Genomics Analysis, ABLife Inc., Wuhan, PR China
| | - Wenqing Wang
- Center for Genomics Analysis, ABLife Inc., Wuhan, PR China
| | - Zisong Zhou
- Center for Genomics Analysis, ABLife Inc., Wuhan, PR China
| | - Wenliang Rang
- Laboratory for Genomics Regulation and Human Health, ABLife Inc., Wuhan, PR China
| | - Han Yu
- Laboratory for Genomics Regulation and Human Health, ABLife Inc., Wuhan, PR China
| | - Yaxun Wei
- Center for Genomics Analysis, ABLife Inc., Wuhan, PR China
| | - Qijia Wu
- Laboratory for Genomics Regulation and Human Health, ABLife Inc., Wuhan, PR China
| | - Yi Zhang
- Laboratory for Genomics Regulation and Human Health, ABLife Inc., Wuhan, PR China
- ABLife BioBigData Instibute, Wuhan, PR China
- Center for Genomics Analysis, ABLife Inc., Wuhan, PR China
| |
Collapse
|
86
|
Bradshaw WJ, Valenzano DR. Extreme genomic volatility characterizes the evolution of the immunoglobulin heavy chain locus in cyprinodontiform fishes. Proc Biol Sci 2020; 287:20200489. [PMID: 32396805 PMCID: PMC7287348 DOI: 10.1098/rspb.2020.0489] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Accepted: 04/14/2020] [Indexed: 12/30/2022] Open
Abstract
The evolution of the adaptive immune system has provided vertebrates with a uniquely sophisticated immune toolkit, enabling them to mount precise immune responses against a staggeringly diverse range of antigens. Like other vertebrates, teleost fishes possess a complex and functional adaptive immune system; however, our knowledge of the complex antigen-receptor genes underlying its functionality has been restricted to a small number of experimental and agricultural species, preventing systematic investigation into how these crucial gene loci evolve. Here, we analyse the genomic structure of the immunoglobulin heavy chain (IGH) gene loci in the cyprinodontiforms, a diverse and important group of teleosts present in many different habitats across the world. We reconstruct the complete IGH loci of the turquoise killifish (Nothobranchius furzeri) and the southern platyfish (Xiphophorus maculatus) and analyse their in vivo gene expression, revealing the presence of species-specific splice isoforms of transmembrane IGHM. We further characterize the IGH constant regions of 10 additional cyprinodontiform species, including guppy, Amazon molly, mummichog and mangrove killifish. Phylogenetic analysis of these constant regions suggests multiple independent rounds of duplication and deletion of the teleost-specific antibody class IGHZ in the cyprinodontiform lineage, demonstrating the extreme volatility of IGH evolution. Focusing on the cyprinodontiforms as a model taxon for comparative evolutionary immunology, this work provides novel genomic resources for studying adaptive immunity and sheds light on the evolutionary history of the adaptive immune system.
Collapse
Affiliation(s)
- William J. Bradshaw
- Max Planck Institute for Biology of Ageing, Joseph-Stelzmann-Str. 296, 50937 Cologne, Germany
- CECAD Research Center, University of Cologne, Joseph-Stelzmann-Str. 26, 50937 Cologne, Germany
| | - Dario Riccardo Valenzano
- Max Planck Institute for Biology of Ageing, Joseph-Stelzmann-Str. 296, 50937 Cologne, Germany
- CECAD Research Center, University of Cologne, Joseph-Stelzmann-Str. 26, 50937 Cologne, Germany
| |
Collapse
|
87
|
Characterization and Functional Analysis of Polyadenylation Sites in Fast and Slow Muscles. BIOMED RESEARCH INTERNATIONAL 2020; 2020:2626584. [PMID: 32258109 PMCID: PMC7102456 DOI: 10.1155/2020/2626584] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 12/01/2019] [Accepted: 01/16/2020] [Indexed: 12/05/2022]
Abstract
Many increasing documents have proved that alternative polyadenylation (APA) events with different polyadenylation sites (PAS) contribute to posttranscriptional regulation. However, little is known about the detailed molecular features of PASs and its role in porcine fast and slow skeletal muscles through microRNAs (miRNAs) and RNA binding proteins (RBPs). In this study, we combined single-molecule real-time sequencing and Illumina RNA-seq datasets to comprehensively analyze polyadenylation in pigs. We identified a total of 10,334 PASs, of which 8734 were characterized by reference genome annotation. 32.86% of PAS-associated genes were determined to have more than one PAS. Further analysis demonstrated that tissue-specific PASs between fast and slow muscles were enriched in skeletal muscle development pathways. In addition, we obtained 1407 target genes regulated by APA events through potential binding 69 miRNAs and 28 RBPs in variable 3′ UTR regions and some are involved in myofiber transformation. Furthermore, the de novo motif search confirmed that the most common usage of canonical motif AAUAAA and three types of PASs may be related to the strength of motifs. In summary, our results provide a useful annotation of PASs for pig transcriptome and suggest that APA may serve as a role in fast and slow muscle development under the regulation of miRNAs and RBPs.
Collapse
|
88
|
Yang F, Wang W, Cetinbas M, Sadreyev RI, Blower MD. Genome-wide analysis identifies cis-acting elements regulating mRNA polyadenylation and translation during vertebrate oocyte maturation. RNA (NEW YORK, N.Y.) 2020; 26:324-344. [PMID: 31896558 PMCID: PMC7025505 DOI: 10.1261/rna.073247.119] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2019] [Accepted: 12/30/2019] [Indexed: 05/10/2023]
Abstract
Most cells change patterns of gene expression through transcriptional regulation. In contrast, oocytes are transcriptionally silent and regulate mRNA poly(A) tail length to control protein production. However, the genome-wide relationship of poly(A) tail changes to mRNA translation during vertebrate oocyte maturation is not known. We used Tail-seq and polyribosome analysis to measure poly(A) tail and translational changes during oocyte maturation in Xenopus laevis We identified large-scale poly(A) and translational changes during oocyte maturation, with poly(A) tail length changes preceding translational changes. Proteins important for completion of the meiotic divisions and early development exhibited increased polyadenylation and translation during oocyte maturation. A family of U-rich sequence elements was enriched near the polyadenylation signal of polyadenylated and translationally activated mRNAs. We propose that changes in mRNA polyadenylation are a conserved mechanism regulating protein expression during vertebrate oocyte maturation and that these changes are controlled by a spatial code of cis-acting sequence elements.
Collapse
Affiliation(s)
- Fei Yang
- Department of Molecular Biology, Massachusetts General Hospital, Boston, Massachusetts 02114, USA
- Department of Genetics, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Wei Wang
- Department of Molecular Biology, Massachusetts General Hospital, Boston, Massachusetts 02114, USA
- Department of Genetics, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Murat Cetinbas
- Department of Molecular Biology, Massachusetts General Hospital, Boston, Massachusetts 02114, USA
- Department of Genetics, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Ruslan I Sadreyev
- Department of Molecular Biology, Massachusetts General Hospital, Boston, Massachusetts 02114, USA
- Department of Pathology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Michael D Blower
- Department of Molecular Biology, Massachusetts General Hospital, Boston, Massachusetts 02114, USA
- Department of Genetics, Harvard Medical School, Boston, Massachusetts 02115, USA
| |
Collapse
|
89
|
Zhu S, Ye W, Ye L, Fu H, Ye C, Xiao X, Ji Y, Lin W, Ji G, Wu X. PlantAPAdb: A Comprehensive Database for Alternative Polyadenylation Sites in Plants. PLANT PHYSIOLOGY 2020; 182:228-242. [PMID: 31767692 PMCID: PMC6945835 DOI: 10.1104/pp.19.00943] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Accepted: 11/18/2019] [Indexed: 05/23/2023]
Abstract
Alternative cleavage and polyadenylation (APA) is increasingly recognized as an important regulatory mechanism in eukaryotic gene expression and is dynamically modulated in a developmental, tissue-specific, or environmentally responsive manner. Given the functional importance of APA and the rapid accumulation of APA sites in plants, a comprehensive and easily accessible APA site database is necessary for improved understanding of APA-mediated gene expression regulation. We present a database called PlantAPAdb that catalogs the most comprehensive APA site data derived from sequences from diverse 3' sequencing protocols and biological samples in plants. Currently, PlantAPAdb contains APA sites in six species, Oryza sativa (japonica and indica), Arabidopsis (Arabidopsis thaliana), Medicago truncatula, Trifolium pratense, Phyllostachys edulis, and Chlamydomonas reinhardtii APA sites in PlantAPAdb are available for bulk download and can be queried in a Google-like manner. PlantAPAdb provides rich information of the whole-genome APA sites, including genomic locations, heterogeneous cleavage sites, expression levels, and sample information. It also provides comprehensive poly(A) signals for APA sites in different genomic regions according to distinct profiles of cis-elements in plants. In addition, PlantAPAdb contains events of 3' untranslated region shortening/lengthening resulting from APA, which helps to understand the mechanisms underlying systematic changes in 3' untranslated region lengths. Additional information about conservation of APA sites in plants is also available, providing insights into the evolutionary polyadenylation configuration across species. As a user-friendly database, PlantAPAdb is a large and extendable resource for elucidating APA mechanisms, APA conservation, and gene expression regulation.
Collapse
Affiliation(s)
- Sheng Zhu
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
- National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, Fujian 3611002, China
| | - Wenbin Ye
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
- National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, Fujian 3611002, China
| | - Lishan Ye
- National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, Fujian 3611002, China
- Xiamen Health and Medical Big Data Center, Xiamen, Fujian 361008, China
| | - Hongjuan Fu
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
- National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, Fujian 3611002, China
| | - Congting Ye
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian 361102, China
| | - Xuesong Xiao
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
| | - Yuanhaowei Ji
- School of Mathematics, Northwest University, Xi'an, Shanxi 710127, China
| | - Weixu Lin
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
| | - Guoli Ji
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
- National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, Fujian 3611002, China
| | - Xiaohui Wu
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
- National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, Fujian 3611002, China
| |
Collapse
|
90
|
Ikert H, Craig PM. Chronic exposure to venlafaxine and increased water temperature reversibly alters microRNA in zebrafish gonads (Danio rerio). COMPARATIVE BIOCHEMISTRY AND PHYSIOLOGY D-GENOMICS & PROTEOMICS 2019; 33:100634. [PMID: 31715506 DOI: 10.1016/j.cbd.2019.100634] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Revised: 10/09/2019] [Accepted: 10/12/2019] [Indexed: 01/28/2023]
Abstract
MicroRNA (miRNA) are short, non-coding RNA that act by downregulating targeted mRNA transcripts. Only recently have they been used as endpoints in studies of aquatic toxicology. The purpose of this study was to determine the effect of an antidepressant contaminant, venlafaxine (VFX), and increased temperature on specific microRNA levels in zebrafish (Danio rerio) reproductive tissue. Adult zebrafish were exposed to one of four conditions; control, 1 μg/L VFX (VFX), 32 °C (Temp), or 1 μg/L VFX + 32 °C (VFX & Temp) for 21 days. Half of the fish were returned to control conditions for a 21-day recovery period. RT-qPCR was performed to measure relative abundances of several miRNAs known to respond to antidepressant exposure: dre-miR-22b-3p, dre-miR-301a, dre-miR-140-5p, dre-let-7d-5p, dre-miR-210-5p, and dre-miR-457b-5p. After the exposure period, dre-miR-22b-3p and dre-miR-301a showed a significant downregulation in response to all treatments. In contrast, after the recovery period, there were no significant differences in microRNA abundance. These altered microRNA are predicted to target several genes, including phosphofructokinase, and are associated with ovarian pathologies. Combined, we have shown that VFX and increased water temperature alter miRNA abundances in zebrafish reproductive tissue, an effect correlated with a functional stress response and cell cycle dysregulation.
Collapse
Affiliation(s)
- Heather Ikert
- Department of Biology, University of Waterloo, 200 University Ave. W., Waterloo, Ontario N2L 3G1, Canada.
| | - Paul M Craig
- Department of Biology, University of Waterloo, 200 University Ave. W., Waterloo, Ontario N2L 3G1, Canada.
| |
Collapse
|
91
|
Ciolli Mattioli C, Rom A, Franke V, Imami K, Arrey G, Terne M, Woehler A, Akalin A, Ulitsky I, Chekulaeva M. Alternative 3' UTRs direct localization of functionally diverse protein isoforms in neuronal compartments. Nucleic Acids Res 2019; 47:2560-2573. [PMID: 30590745 PMCID: PMC6411841 DOI: 10.1093/nar/gky1270] [Citation(s) in RCA: 83] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Revised: 11/27/2018] [Accepted: 12/13/2018] [Indexed: 01/01/2023] Open
Abstract
The proper subcellular localization of RNAs and local translational regulation is crucial in highly compartmentalized cells, such as neurons. RNA localization is mediated by specific cis-regulatory elements usually found in mRNA 3′UTRs. Therefore, processes that generate alternative 3′UTRs—alternative splicing and polyadenylation—have the potential to diversify mRNA localization patterns in neurons. Here, we performed mapping of alternative 3′UTRs in neurites and soma isolated from mESC-derived neurons. Our analysis identified 593 genes with differentially localized 3′UTR isoforms. In particular, we have shown that two isoforms of Cdc42 gene with distinct functions in neuronal polarity are differentially localized between neurites and soma of mESC-derived and mouse primary cortical neurons, at both mRNA and protein level. Using reporter assays and 3′UTR swapping experiments, we have identified the role of alternative 3′UTRs and mRNA transport in differential localization of alternative CDC42 protein isoforms. Moreover, we used SILAC to identify isoform-specific Cdc42 3′UTR-bound proteome with potential role in Cdc42 localization and translation. Our analysis points to usage of alternative 3′UTR isoforms as a novel mechanism to provide for differential localization of functionally diverse alternative protein isoforms.
Collapse
Affiliation(s)
- Camilla Ciolli Mattioli
- Non-coding RNAs and mechanisms of cytoplasmic gene regulation, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Aviv Rom
- Weizmann Institute of Science, Rehovot, Israel
| | - Vedran Franke
- BIMSB Bioinformatics platform, Max Delbrück Center for Molecular Medicine, Berlin 13125, Germany
| | - Koshi Imami
- Proteome Dynamics, Max Delbrück Center for Molecular Medicine, Berlin 13125, Germany
| | - Gerard Arrey
- Non-coding RNAs and mechanisms of cytoplasmic gene regulation, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Mandy Terne
- Developmental Biology / Signal Transduction, Max Delbrück Center for Molecular Medicine, Berlin 13125, Germany
| | - Andrew Woehler
- BIMSB Light Microscopy platform, Max Delbrück Center for Molecular Medicine, Berlin 13125, Germany
| | - Altuna Akalin
- BIMSB Bioinformatics platform, Max Delbrück Center for Molecular Medicine, Berlin 13125, Germany
| | | | - Marina Chekulaeva
- Non-coding RNAs and mechanisms of cytoplasmic gene regulation, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| |
Collapse
|
92
|
Yu Z, Lin J, Li QQ. Transcriptome Analyses of FY Mutants Reveal Its Role in mRNA Alternative Polyadenylation. THE PLANT CELL 2019; 31:2332-2352. [PMID: 31427469 PMCID: PMC6790095 DOI: 10.1105/tpc.18.00545] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 07/15/2019] [Accepted: 08/19/2019] [Indexed: 05/10/2023]
Abstract
A crucial step for mRNA polyadenylation is poly(A) signal recognition by trans-acting factors. The mammalian cleavage and polyadenylation specificity factor (CPSF) complex components CPSF30 and WD repeat-containing protein33 (WDR33) recognize the canonical AAUAAA for polyadenylation. In Arabidopsis (Arabidopsis thaliana), the flowering time regulator FY is the homolog of WDR33. However, its role in mRNA polyadenylation is poorly understood. Using poly(A) tag sequencing, we found that >50% of alternative polyadenylation (APA) events are altered in fy single mutants or double mutants with oxt6 (a null mutant of AtCPSF30), but mutation of the FY WD40-repeat has a stronger effect than deletion of the plant-unique Pro-Pro-Leu-Pro-Pro (PPLPP) domain. fy mutations disrupt AAUAAA or AAUAAA-like poly(A) signal recognition. Notably, A-rich signal usage is suppressed in the WD40-repeat mutation but promoted in PPLPP-domain deficiency. However, fy mutations do not aggravate the altered signal usage in oxt6 Furthermore, the WD40-repeat mutation shows a preference for 3' untranslated region shortening, but the PPLPP-domain deficiency shows a preference for lengthening. Interestingly, the WD40-repeat mutant exhibits shortened primary roots and late flowering with alteration of APA of related genes. Importantly, the long transcripts of two APA genes affected in fy are related to abiotic stress responses. These results reveal a conserved and specific role of FY in mRNA polyadenylation.
Collapse
Affiliation(s)
- Zhibo Yu
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian, China 361102
| | - Juncheng Lin
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian, China 361102
- Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, California 91766
| | - Qingshun Quinn Li
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian, China 361102
- Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, California 91766
| |
Collapse
|
93
|
Shui Y, Liu GF, Xu ZH, Zhu GY. Exploring potential proteins associated with cyclin B 3'UTR in Procambarus clarkii oocytes. Biochem Biophys Res Commun 2019; 517:458-462. [PMID: 31376940 DOI: 10.1016/j.bbrc.2019.07.104] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 07/26/2019] [Indexed: 10/26/2022]
Abstract
Cyclin B is a ubiquitous regulatory molecule and has been implicated in mitosis and meiosis in oocytes. Phenomenon that differ in the length of cyclin B 3'UTR in crustacean has attracted much attention, although molecular details are poorly understood. The study of 3'UTR-interacting proteins could yield much information in translational regulation and the mRNA localization process. Previous studies on crayfish suggested that the 3'UTR (1300 bp) probably contains the potential regulatory sequence/motifs such as CPEs and K-box et al. In present study, using pull-down assay coupled with mass spectrometry approach allowing us to explore the potential proteins associated with the 3'UTR. We finally identified four candidate proteins including Hspg 2, Vtg, eef1a and Tuba1a, which annotated as significant roles involved in cell differentiation, lipid transporter activity, and meiotic cell cycle process. The preliminary results will contribute to the advance in understanding the translational activation of cyclin B in oocyte maturation regulation in crustacean.
Collapse
Affiliation(s)
- Yan Shui
- Key Laboratory of Freshwater Fisheries and Germplasm Resources Utilization, Ministry of Agriculture, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi, 214081, PR China; Wuxi Fishery College, Nanjing Agricultural University, Wuxi, 214081, PR China.
| | - Guo-Feng Liu
- Key Laboratory of Freshwater Fisheries and Germplasm Resources Utilization, Ministry of Agriculture, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi, 214081, PR China; Wuxi Fishery College, Nanjing Agricultural University, Wuxi, 214081, PR China
| | - Zeng-Hong Xu
- Key Laboratory of Freshwater Fisheries and Germplasm Resources Utilization, Ministry of Agriculture, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi, 214081, PR China; Wuxi Fishery College, Nanjing Agricultural University, Wuxi, 214081, PR China
| | - Guang-Yan Zhu
- Key Laboratory of Freshwater Fisheries and Germplasm Resources Utilization, Ministry of Agriculture, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi, 214081, PR China
| |
Collapse
|
94
|
Arefeen A, Liu J, Xiao X, Jiang T. TAPAS: tool for alternative polyadenylation site analysis. Bioinformatics 2019; 34:2521-2529. [PMID: 30052912 DOI: 10.1093/bioinformatics/bty110] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2017] [Accepted: 02/22/2018] [Indexed: 01/08/2023] Open
Abstract
Motivation The length of the 3' untranslated region (3' UTR) of an mRNA is essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, correlation between diseases and the shortening (or lengthening) of 3' UTRs has been reported in the literature. This length is largely determined by the polyadenylation cleavage site in the mRNA. As alternative polyadenylation (APA) sites are common in mammalian genes, several tools have been published recently for detecting APA sites from RNA-Seq data or performing shortening/lengthening analysis. These tools consider either up to only two APA sites in a gene or only APA sites that occur in the last exon of a gene, although a gene may generally have more than two APA sites and an APA site may sometimes occur before the last exon. Furthermore, the tools are unable to integrate the analysis of shortening/lengthening events with APA site detection. Results We propose a new tool, called TAPAS, for detecting novel APA sites from RNA-Seq data. It can deal with more than two APA sites in a gene as well as APA sites that occur before the last exon. The tool is based on an existing method for finding change points in time series data, but some filtration techniques are also adopted to remove change points that are likely false APA sites. It is then extended to identify APA sites that are expressed differently between two biological samples and genes that contain 3' UTRs with shortening/lengthening events. Our extensive experiments on simulated and real RNA-Seq data demonstrate that TAPAS outperforms the existing tools for APA site detection or shortening/lengthening analysis significantly. Availability and implementation https://github.com/arefeen/TAPAS. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Ashraful Arefeen
- Department of Computer Science and Engineering, University of California, Riverside, CA, USA
| | - Juntao Liu
- School of Mathematics, Shandong University, Jinan, Shandong, China
| | - Xinshu Xiao
- Department of Integrative Biology and Physiology, University of California, Los Angeles, CA, USA
| | - Tao Jiang
- Department of Computer Science and Engineering, University of California, Riverside, CA, USA.,Institute of Integrative Genome Biology, University of California, Riverside, CA, USA.,MOE Key Lab of Bioinformatics and Bioinformatics Division, TNLIST/Department of Computer Science and Technology, Tsinghua University, Beijing, China
| |
Collapse
|
95
|
Genome-Wide Profiling of Polyadenylation Events in Maize Using High-Throughput Transcriptomic Sequences. G3-GENES GENOMES GENETICS 2019; 9:2749-2760. [PMID: 31239292 PMCID: PMC6686930 DOI: 10.1534/g3.119.400196] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Polyadenylation is an essential post-transcriptional modification of eukaryotic transcripts that plays critical role in transcript stability, localization, transport, and translational efficiency. About 70% genes in plants contain alternative polyadenylation (APA) sites. Despite availability of vast amount of sequencing data, to date, a comprehensive map of the polyadenylation events in maize is not available. Here, 9.48 billion RNA-Seq reads were analyzed to characterize 95,345 Poly(A) Clusters (PAC) in 23,705 (51%) maize genes. Of these, 76% were APA genes. However, most APA genes (55%) expressed a dominant PAC rather than favoring multiple PACs equally. The lincRNA genes with PACs were significantly longer in length than the genes without any PAC and about 48% genes had APA sites. Heterogeneity was observed in 52% of the PACs supporting the imprecise nature of the polyadenylation process. Genomic distribution revealed that the majority of the PACs (78%) were located in the genic regions. Unlike previous studies, large number of PACs were observed in the intergenic (n = 21,264), 5′-UTR (735), CDS (2,542), and the intronic regions (12,841). The CDS and introns with PACs were longer in length than without PACs, whereas intergenic PACs were more often associated with transcripts that lacked annotated 3′-UTRs. Nucleotide composition around PACs demonstrated AT-richness and the common upstream motif was AAUAAA, which is consistent with other plants. According to this study, only 2,830 genes still maintained the use of AAUAAA motif. This large-scale data provides useful insights about the gene expression regulation and could be utilized as evidence to validate the annotation of transcript ends.
Collapse
|
96
|
Sarropoulou E, Kaitetzidou E, Papandroulakis N, Tsalafouta A, Pavlidis M. Inventory of European Sea Bass ( Dicentrarchus labrax) sncRNAs Vital During Early Teleost Development. Front Genet 2019; 10:657. [PMID: 31404269 PMCID: PMC6670005 DOI: 10.3389/fgene.2019.00657] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 06/21/2019] [Indexed: 01/25/2023] Open
Abstract
During early animal ontogenesis, a plethora of small non-coding RNAs (sncRNAs) are greatly expressed and have been shown to be involved in several regulatory pathways vital to proper development. The rapid advancements in sequencing and computing methodologies in the last decade have paved the way for the production of sequencing data in a broad range of organisms, including teleost species. Consequently, this has led to the discovery of sncRNAs as well as the potentially novel roles of sncRNA in gene regulation. Among the several classes of sncRNAs, microRNAs (miRNAs) have, in particular, been shown to play a key role in development. The present work aims to identify the miRNAs that play important roles during early European sea bass (Dicentrarchus labrax) development. The European sea bass is a species of high commercial impact in European and especially Mediterranean aquaculture. This study reports, for the first time, the identification and characterization of small RNAs that play a part in the 10 developmental stages (from morula to all fins) of the European sea bass. From 10 developmental stages, more than 135 million reads, generated by next-generation sequencing, were retrieved from publicly available databases as well as newly generated. The analysis resulted in about 2,000 sample grouped reads, and their subsequently annotation revealed that the majority of transcripts belonged to the class of miRNAs followed by small nuclear RNAs and small nucleolar RNAs. The analysis of small RNA expression among the developmental stages under study revealed that miRNAs are active throughout development, with the main activity occurring after the earlier stages (morula and 50% epiboly) and at the later stages (first feeding, flexion, and all fins). Furthermore, investigating miRNAs exclusively expressed in one of the stages unraveled five miRNAs with a higher abundance only in the morula stage (miR-155, miR-430a, d1, d2, and miR-458), indicating possible important key roles of those miRNAs in further embryonic development. An additional target search showed putative miRNA-mRNA interactions with possible direct and indirect regulatory functions of the identified miRNAs.
Collapse
Affiliation(s)
- Elena Sarropoulou
- Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Center for Marine Research, Heraklion, Greece
| | - Elizabet Kaitetzidou
- Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Center for Marine Research, Heraklion, Greece
| | - Nikos Papandroulakis
- Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Center for Marine Research, Heraklion, Greece
| | | | | |
Collapse
|
97
|
Chen M, Ji G, Fu H, Lin Q, Ye C, Ye W, Su Y, Wu X. A survey on identification and quantification of alternative polyadenylation sites from RNA-seq data. Brief Bioinform 2019; 21:1261-1276. [PMID: 31267126 DOI: 10.1093/bib/bbz068] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2019] [Revised: 05/03/2019] [Accepted: 05/14/2019] [Indexed: 12/13/2022] Open
Abstract
Alternative polyadenylation (APA) has been implicated to play an important role in post-transcriptional regulation by regulating mRNA abundance, stability, localization and translation, which contributes considerably to transcriptome diversity and gene expression regulation. RNA-seq has become a routine approach for transcriptome profiling, generating unprecedented data that could be used to identify and quantify APA site usage. A number of computational approaches for identifying APA sites and/or dynamic APA events from RNA-seq data have emerged in the literature, which provide valuable yet preliminary results that should be refined to yield credible guidelines for the scientific community. In this review, we provided a comprehensive overview of the status of currently available computational approaches. We also conducted objective benchmarking analysis using RNA-seq data sets from different species (human, mouse and Arabidopsis) and simulated data sets to present a systematic evaluation of 11 representative methods. Our benchmarking study showed that the overall performance of all tools investigated is moderate, reflecting that there is still lot of scope to improve the prediction of APA site or dynamic APA events from RNA-seq data. Particularly, prediction results from individual tools differ considerably, and only a limited number of predicted APA sites or genes are common among different tools. Accordingly, we attempted to give some advice on how to assess the reliability of the obtained results. We also proposed practical recommendations on the appropriate method applicable to diverse scenarios and discussed implications and future directions relevant to profiling APA from RNA-seq data.
Collapse
Affiliation(s)
- Moliang Chen
- Department of Automation, Xiamen University, Xiamen 361005, China.,Xiamen Research Institute of National Center of Healthcare Big Data, Xiamen 361005, China
| | - Guoli Ji
- Department of Automation, Xiamen University, Xiamen 361005, China.,Xiamen Research Institute of National Center of Healthcare Big Data, Xiamen 361005, China
| | - Hongjuan Fu
- Department of Automation, Xiamen University, Xiamen 361005, China.,Xiamen Research Institute of National Center of Healthcare Big Data, Xiamen 361005, China
| | - Qianmin Lin
- Xiang' an hospital of Xiamen university, Xiamen 361005, China
| | - Congting Ye
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian 361102, China
| | - Wenbin Ye
- Department of Automation, Xiamen University, Xiamen 361005, China.,Xiamen Research Institute of National Center of Healthcare Big Data, Xiamen 361005, China
| | - Yaru Su
- College of Mathematics and Computer Science, Fuzhou University, Fuzhou 350116, China
| | - Xiaohui Wu
- Department of Automation, Xiamen University, Xiamen 361005, China.,Xiamen Research Institute of National Center of Healthcare Big Data, Xiamen 361005, China
| |
Collapse
|
98
|
Ye C, Long Y, Ji G, Li QQ, Wu X. APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data. Bioinformatics 2019; 34:1841-1849. [PMID: 29360928 DOI: 10.1093/bioinformatics/bty029] [Citation(s) in RCA: 74] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2017] [Accepted: 01/17/2018] [Indexed: 12/28/2022] Open
Abstract
Motivation Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites. Results We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome. Availability and implementation Freely available for download at https://apatrap.sourceforge.io. Contact liqq@xmu.edu.cn or xhuister@xmu.edu.cn. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Congting Ye
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian 361102, China
| | - Yuqi Long
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
| | - Guoli Ji
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
| | - Qingshun Quinn Li
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, Fujian 361102, China.,Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, CA 91766, USA
| | - Xiaohui Wu
- Department of Automation, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
99
|
Zhou X, Zhang Y, Michal JJ, Qu L, Zhang S, Wildung MR, Du W, Pouchnik DJ, Zhao H, Xia Y, Shi H, Ji G, Davis JF, Smith GD, Griswold MD, Harland RM, Jiang Z. Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis. Cell Mol Life Sci 2019; 76:2185-2198. [PMID: 30729254 PMCID: PMC6597005 DOI: 10.1007/s00018-019-03036-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Revised: 01/09/2019] [Accepted: 01/30/2019] [Indexed: 12/27/2022]
Abstract
RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.
Collapse
Affiliation(s)
- Xiang Zhou
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, 99164-7620, USA
- College of Animal Sciences and Veterinary Medicine, Huazhong Agricultural University, Wuhan, China
| | - Yangzi Zhang
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, 99164-7620, USA
| | - Jennifer J Michal
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, 99164-7620, USA
| | - Lujiang Qu
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, 99164-7620, USA
- College of Animal Sciences and Technology, China Agricultural University, Beijing, China
| | - Shuwen Zhang
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, 99164-7620, USA
| | - Mark R Wildung
- Laboratory for Biotechnology and Bioanalysis, Center for Reproductive Biology, Washington State University, Pullman, WA, USA
| | - Weiwei Du
- Laboratory for Biotechnology and Bioanalysis, Center for Reproductive Biology, Washington State University, Pullman, WA, USA
| | - Derek J Pouchnik
- Laboratory for Biotechnology and Bioanalysis, Center for Reproductive Biology, Washington State University, Pullman, WA, USA
| | - Hui Zhao
- School of Biomedical Sciences, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, China
| | - Yin Xia
- School of Biomedical Sciences, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, China
| | - Honghua Shi
- State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai, China
| | - Guoli Ji
- Department of Automation, Xiamen University, Xiamen, China
| | - Jon F Davis
- Department of Integrative Physiology and Neuroscience, Washington State University, Pullman, WA, USA
| | - Gary D Smith
- Departments of OB/GYN, Physiology, and Urology, University of Michigan, Ann Arbor, MI, USA
| | - Michael D Griswold
- Laboratory for Biotechnology and Bioanalysis, Center for Reproductive Biology, Washington State University, Pullman, WA, USA
| | - Richard M Harland
- Department of Molecular and Cell Biology, University of California Berkeley, Berkeley, CA, USA
| | - Zhihua Jiang
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, 99164-7620, USA.
| |
Collapse
|
100
|
Cao J, Ye C, Hao G, Dabney-Smith C, Hunt AG, Li QQ. Root Hair Single Cell Type Specific Profiles of Gene Expression and Alternative Polyadenylation Under Cadmium Stress. FRONTIERS IN PLANT SCIENCE 2019; 10:589. [PMID: 31134121 PMCID: PMC6523994 DOI: 10.3389/fpls.2019.00589] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2019] [Accepted: 04/23/2019] [Indexed: 05/15/2023]
Abstract
Transcriptional networks are tightly controlled in plant development and stress responses. Alternative polyadenylation (APA) has been found to regulate gene expression under abiotic stress by increasing the heterogeneity at mRNA 3'-ends. Heavy metals like cadmium pollute water and soil due to mining and industry applications. Understanding how plants cope with heavy metal stress remains an interesting question. The Arabidopsis root hair was chosen as a single cell model to investigate the functional role of APA in cadmium stress response. Primary root growth inhibition and defective root hair morphotypes were observed. Poly(A) tag (PAT) libraries from single cell types, i.e., root hair cells, non-hair epidermal cells, and whole root tip under cadmium stress were prepared and sequenced. Interestingly, a root hair cell type-specific gene expression under short term cadmium exposure, but not related to the prolonged treatment, was detected. Differentially expressed poly(A) sites were identified, which largely contributed to altered gene expression, and enriched in pentose and glucuronate interconversion pathways as well as phenylpropanoid biosynthesis pathways. Numerous genes with poly(A) site switching were found, particularly for functions in cell wall modification, root epidermal differentiation, and root hair tip growth. Our findings suggest that APA plays a functional role as a potential stress modulator in root hair cells under cadmium treatment.
Collapse
Affiliation(s)
- Jingyi Cao
- Department of Biology, Miami University, Oxford, OH, United States
- Cell, Molecular, Structural Biology Graduate Program, Miami University, Oxford, OH, United States
| | - Congting Ye
- Key Laboratory of Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, China
| | - Guijie Hao
- Department of Plant and Soil Sciences, University of Kentucky, Lexington, KY, United States
| | - Carole Dabney-Smith
- Cell, Molecular, Structural Biology Graduate Program, Miami University, Oxford, OH, United States
- Department of Chemistry and Biochemistry, Miami University, Oxford, OH, United States
| | - Arthur G. Hunt
- Department of Plant and Soil Sciences, University of Kentucky, Lexington, KY, United States
| | - Qingshun Q. Li
- Department of Biology, Miami University, Oxford, OH, United States
- Cell, Molecular, Structural Biology Graduate Program, Miami University, Oxford, OH, United States
- Key Laboratory of Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, China
- Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, CA, United States
| |
Collapse
|