1
|
Carpinteyro-Ponce J, Machado CA. The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes. Genome Biol Evol 2024; 16:evae047. [PMID: 38482945 PMCID: PMC10980976 DOI: 10.1093/gbe/evae047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/07/2024] [Indexed: 04/01/2024] Open
Abstract
Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.
Collapse
Affiliation(s)
| | - Carlos A Machado
- Department of Biology, University of Maryland, College Park, MD, USA
| |
Collapse
|
2
|
Zhu F, Lu J, Sun K, Deng C, Xu Y. Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing. BMC Genom Data 2024; 25:23. [PMID: 38408920 PMCID: PMC10895795 DOI: 10.1186/s12863-024-01208-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 02/14/2024] [Indexed: 02/28/2024] Open
Abstract
BACKGROUND Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates. METHODS In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species. RESULTS A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.
Collapse
Affiliation(s)
- Fei Zhu
- School of Life Sciences, Guizhou Normal University, 550025, Guiyang, Guizhou, China.
| | - Jing Lu
- School of Life Sciences, Guizhou Normal University, 550025, Guiyang, Guizhou, China
| | - Ke Sun
- School of Life Sciences, Guizhou Normal University, 550025, Guiyang, Guizhou, China
| | - Cao Deng
- Department of Bioinformatics, DNA Stories Bioinformatics Center, 610000, Chengdu, China
| | - Yu Xu
- School of Life Sciences, Guizhou Normal University, 550025, Guiyang, Guizhou, China
| |
Collapse
|
3
|
Ma J, Zhao X, Qi E, Han R, Yu T, Li G. Highly efficient clustering of long-read transcriptomic data with GeLuster. Bioinformatics 2024; 40:btae059. [PMID: 38310330 PMCID: PMC10881092 DOI: 10.1093/bioinformatics/btae059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 01/08/2024] [Accepted: 01/30/2024] [Indexed: 02/05/2024] Open
Abstract
MOTIVATION The advancement of long-read RNA sequencing technologies leads to a bright future for transcriptome analysis, in which clustering long reads according to their gene family of origin is of great importance. However, existing de novo clustering algorithms require plenty of computing resources. RESULTS We developed a new algorithm GeLuster for clustering long RNA-seq reads. Based on our tests on one simulated dataset and nine real datasets, GeLuster exhibited superior performance. On the tested Nanopore datasets it ran 2.9-17.5 times as fast as the second-fastest method with less than one-seventh of memory consumption, while achieving higher clustering accuracy. And on the PacBio data, GeLuster also had a similar performance. It sets the stage for large-scale transcriptome study in future. AVAILABILITY AND IMPLEMENTATION GeLuster is freely available at https://github.com/yutingsdu/GeLuster.
Collapse
Affiliation(s)
- Junchi Ma
- Research Center for Mathematics and Interdisciplinary Sciences (Frontiers Science Center for Nonlinear Expectations), Shandong University, Qingdao 266237, China
- School of Mathematics, Shandong University, Jinan, Shandong 250100, China
| | - Xiaoyu Zhao
- School of Mathematics, Shandong University, Jinan, Shandong 250100, China
| | - Enfeng Qi
- School of Mathematics and Statistics, Guangxi Normal University, Guilin 541000, China
| | - Renmin Han
- Research Center for Mathematics and Interdisciplinary Sciences (Frontiers Science Center for Nonlinear Expectations), Shandong University, Qingdao 266237, China
| | - Ting Yu
- Research Center for Mathematics and Interdisciplinary Sciences (Frontiers Science Center for Nonlinear Expectations), Shandong University, Qingdao 266237, China
| | - Guojun Li
- Research Center for Mathematics and Interdisciplinary Sciences (Frontiers Science Center for Nonlinear Expectations), Shandong University, Qingdao 266237, China
| |
Collapse
|
4
|
Wu Z, Miedzinska K, Krause JS, Pérez JH, Wingfield JC, Meddle SL, Smith J. A chromosome-level genome assembly of a free-living white-crowned sparrow (Zonotrichia leucophrys gambelii). Sci Data 2024; 11:86. [PMID: 38238322 PMCID: PMC10796373 DOI: 10.1038/s41597-024-02929-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Accepted: 01/03/2024] [Indexed: 01/22/2024] Open
Abstract
The white-crowned sparrow, Zonotrichia leucophrys, is a passerine bird with a wide distribution and it is extensively adapted to environmental changes. It has historically acted as a model species in studies on avian ecology, physiology and behaviour. Here, we present a high-quality chromosome-level genome of Zonotrichia leucophrys using PacBio and OmniC sequencing data. Gene models were constructed by combining RNA-seq and Iso-seq data from liver, hypothalamus, and ovary. In total a 1,123,996,003 bp genome was generated, including 31 chromosomes assembled in complete scaffolds along with other, unplaced scaffolds. This high-quality genome assembly offers an important genomic resource for the research community using the white-crowned sparrow as a model for understanding avian genome biology and development, and provides a genomic basis for future studies, both fundamental and applied.
Collapse
Affiliation(s)
- Zhou Wu
- The Roslin Institute and Royal (Dick) School of Veterinary Studies R(D)SVS, The University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, UK.
| | - Katarzyna Miedzinska
- The Roslin Institute and Royal (Dick) School of Veterinary Studies R(D)SVS, The University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, UK
| | - Jesse S Krause
- Department of Neurobiology, Physiology, and Behavior, University of California, Davis, CA, 95616, USA
- Department of Biology, University of Nevada Reno, Reno, NV, 89557, USA
| | - Jonathan H Pérez
- Department of Biology, University of South Alabama, Mobile, AL, 36688, USA
| | - John C Wingfield
- Department of Neurobiology, Physiology, and Behavior, University of California, Davis, CA, 95616, USA
| | - Simone L Meddle
- The Roslin Institute and Royal (Dick) School of Veterinary Studies R(D)SVS, The University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, UK
| | - Jacqueline Smith
- The Roslin Institute and Royal (Dick) School of Veterinary Studies R(D)SVS, The University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, UK.
| |
Collapse
|
5
|
Krusnauskas R, Stakaitis R, Steponaitis G, Almstrup K, Vaitkiene P. Identification and comparison of m6A modifications in glioblastoma non-coding RNAs with MeRIP-seq and Nanopore dRNA-seq. Epigenetics 2023; 18:2163365. [PMID: 36597408 PMCID: PMC9980576 DOI: 10.1080/15592294.2022.2163365] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
The most prominent RNA modification - N6-methyladenosine (m6A) - affects gene regulation and cancer progression. The extent and effect of m6A on long non-coding RNAs (lncRNAs) is, however, still not clear. The most established method for m6A detection is methylated RNA immunoprecipitation and sequencing (MeRIP-seq). However, Oxford Nanopore Technologies recently developed direct RNA-seq (dRNA-seq) method, allowing m6A identification at higher resolution and in its native form. We performed whole transcriptome sequencing of the glioblastoma cell line U87-MG with both MeRIP-seq and dRNA-seq. For MeRIP-seq, m6A peaks were identified using nf-core/chipseq, and for dRNA-seq - EpiNano pipeline. MeRIP-seq analysis revealed 5086 lncRNAs transcripts, while dRNA-seq identified 336 lncRNAs transcripts from which 556 and 198 were found to be m6A modified, respectively. While 24 lncRNAs with m6A overlapped between two methods. Gliovis database analysis revealed that the expression of the major part of identified overlapping lncRNAs was associated with glioma grade or patient survival prognosis. We found that the frequency of m6A occurrence in lncRNAs varied more than 9-fold throughout the provided list of 24 modified lncRNAs. The highest m6A frequency was detected in MIR1915HG, THAP9-AS1, MALAT1, NORAD1, and NEAT1 (49-88nt), while MIR99AHG, SNHG3, LOXL1-AS1, ILF3-DT showed the lowest m6A frequency (445-261nt). Taken together, (1) we provide a high accuracy list of 24 m6A modified lncRNAs of U87-MG cells; (2) we conclude that MeRIP-seq is more suitable for an initial m6A screening study, due to its higher lncRNA coverage, whereas dRNA-seq is most useful when more in-depth analysis of m6A quantity and precise location is of interest.Abbreviations: (dRNA-seq) direct RNA-seq, (GBM) glioblastoma, (LGG) low-grade glioma, (lncRNAs) long non-coding RNAs, (m6A) N6-methyladenosine, (MeRIP-seq) methylated RNA immunoprecipitation and sequencing, (ncRNA) non-coding RNA, (ONT) Oxford Nanopore Technologi; Lietuvos Mokslo Taryba.
Collapse
Affiliation(s)
- Raulas Krusnauskas
- Laboratory of Molecular Neurobiology, Neuroscience Institute, Lithuanian University of Health Sciences, Eiveniu str. 4, LT50161, Kaunas, Lithuania
| | - Rytis Stakaitis
- Laboratory of Molecular Neurooncology, Neuroscience Institute, Lithuanian University of Health Sciences, Eiveniu str. 4, LT50161, Kaunas, Lithuania
| | - Giedrius Steponaitis
- Laboratory of Molecular Neurooncology, Neuroscience Institute, Lithuanian University of Health Sciences, Eiveniu str. 4, LT50161, Kaunas, Lithuania
| | - Kristian Almstrup
- Department of Growth and Reproduction, Rigshospitalet, University of Copenhagen, GR-5064, Rigshospitalet, Blegdamsvej 9, DK-2100 Copenhagen, Denmark.,International Center for Research and Research Training in Endocrine Disruption of Male Reproduction and Child Health (Edmarc), Rigshospitalet, University of Copenhagen, GR-5064, Rigshospitalet, Blegdamsvej 9, DK-2100 Copenhagen, Denmark
| | - Paulina Vaitkiene
- Laboratory of Molecular Neurobiology, Neuroscience Institute, Lithuanian University of Health Sciences, Eiveniu str. 4, LT50161, Kaunas, Lithuania
| |
Collapse
|
6
|
Traenkner D, Shennib O, Johnson A, Weinbrom A, Taylor MR, Williams ME. Modular Splicing Is Linked to Evolution in the Synapse-Specificity Molecule Kirrel3. eNeuro 2023; 10:ENEURO.0253-23.2023. [PMID: 37977826 DOI: 10.1523/eneuro.0253-23.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 11/09/2023] [Accepted: 11/13/2023] [Indexed: 11/19/2023] Open
Abstract
Kirrel3 is a cell-adhesion molecule that instructs the formation of specific synapses during brain development in mouse and Kirrel3 variants may be risk factors for autism and intellectual disabilities in humans. Kirrel3 is predicted to undergo alternative splicing but brain isoforms have not been studied. Here, we present the first in-depth characterization of Kirrel3 isoform diversity in brain using targeted, long-read mRNA sequencing of mouse hippocampus. We identified 19 isoforms with predicted transmembrane and secreted forms and show that even rare isoforms generate detectable protein in the brain. We also analyzed publicly-available long-read mRNA databases from human brain tissue and found 11 Kirrel3 isoforms that, similar to mouse, encode transmembrane and secreted forms. In mice and humans, Kirrel3 diversity arises from alternative, independent use of protein-domain coding exons and alternative early translation-stop signals. Intriguingly, the alternatively spliced exons appear at branch points in the chordate phylogenetic tree, including one exon only found in humans and their closest living relatives, the great apes. Together, these results validate a simple pipeline for analyzing isoform diversity in genes with low expression and suggest that Kirrel3 function is fine-tuned by alternative splicing and may play a role in brain evolution.
Collapse
Affiliation(s)
- Dimitri Traenkner
- Department of Neurobiology, University of Utah, School of Medicine, Salt Lake City, UT 84112
| | - Omar Shennib
- Department of Neurobiology, University of Utah, School of Medicine, Salt Lake City, UT 84112
| | - Alyssa Johnson
- Department of Neurobiology, University of Utah, School of Medicine, Salt Lake City, UT 84112
| | - Adam Weinbrom
- Department of Neurobiology, University of Utah, School of Medicine, Salt Lake City, UT 84112
| | - Matthew R Taylor
- Department of Neurobiology, University of Utah, School of Medicine, Salt Lake City, UT 84112
| | - Megan E Williams
- Department of Neurobiology, University of Utah, School of Medicine, Salt Lake City, UT 84112
| |
Collapse
|
7
|
Torre D, Francoeur NJ, Kalma Y, Gross Carmel I, Melo BS, Deikus G, Allette K, Flohr R, Fridrikh M, Vlachos K, Madrid K, Shah H, Wang YC, Sridhar SH, Smith ML, Eliyahu E, Azem F, Amir H, Mayshar Y, Marazzi I, Guccione E, Schadt E, Ben-Yosef D, Sebra R. Isoform-resolved transcriptome of the human preimplantation embryo. Nat Commun 2023; 14:6902. [PMID: 37903791 PMCID: PMC10616205 DOI: 10.1038/s41467-023-42558-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 10/15/2023] [Indexed: 11/01/2023] Open
Abstract
Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.
Collapse
Affiliation(s)
- Denis Torre
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | | | - Yael Kalma
- Fertility and IVF Institute, Tel-Aviv Sourasky Medical Center, Affiliated to Tel Aviv University, Tel Aviv, 64239, Israel
| | - Ilana Gross Carmel
- Fertility and IVF Institute, Tel-Aviv Sourasky Medical Center, Affiliated to Tel Aviv University, Tel Aviv, 64239, Israel
| | - Betsaida S Melo
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Gintaras Deikus
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Kimaada Allette
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Ron Flohr
- Department of Cell and Developmental Biology, Sackler Faculty of Medicine, Sagol School of Neuroscience, Tel-Aviv University, Tel-Aviv, 69978, Israel
- CORAL - Center Of Regeneration and Longevity, Tel-Aviv Sourasky Medical Center, Tel Aviv, 64239, Israel
| | - Maya Fridrikh
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | | | - Kent Madrid
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Hardik Shah
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Ying-Chih Wang
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Shwetha H Sridhar
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Melissa L Smith
- Department of Biochemistry and Molecular Genetics, University of Louisville, Louisville, KY, 40202, USA
| | - Efrat Eliyahu
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Foad Azem
- Fertility and IVF Institute, Tel-Aviv Sourasky Medical Center, Affiliated to Tel Aviv University, Tel Aviv, 64239, Israel
| | - Hadar Amir
- Fertility and IVF Institute, Tel-Aviv Sourasky Medical Center, Affiliated to Tel Aviv University, Tel Aviv, 64239, Israel
| | - Yoav Mayshar
- Department of Molecular Cell Biology, Weizmann Institute of Science, 7610001, Rehovot, Israel
| | - Ivan Marazzi
- Department of Biological Chemistry, Center for Epigenetics and Metabolism, University of California, Irvine, CA, 92697, USA
| | - Ernesto Guccione
- Center for OncoGenomics and Innovative Therapeutics (COGIT); Department of Oncological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Black Family Stem Cell Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Eric Schadt
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Dalit Ben-Yosef
- Fertility and IVF Institute, Tel-Aviv Sourasky Medical Center, Affiliated to Tel Aviv University, Tel Aviv, 64239, Israel.
- Department of Cell and Developmental Biology, Sackler Faculty of Medicine, Sagol School of Neuroscience, Tel-Aviv University, Tel-Aviv, 69978, Israel.
- CORAL - Center Of Regeneration and Longevity, Tel-Aviv Sourasky Medical Center, Tel Aviv, 64239, Israel.
| | - Robert Sebra
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.
- Center for Advanced Genomics Technology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.
- Icahn Genomics Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.
- Black Family Stem Cell Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.
| |
Collapse
|
8
|
Triant DA, Walsh AT, Hartley GA, Petry B, Stegemiller MR, Nelson BM, McKendrick MM, Fuller EP, Cockett NE, Koltes JE, McKay SD, Green JA, Murdoch BM, Hagen DE, Elsik CG. AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes. Mamm Genome 2023; 34:418-436. [PMID: 37460664 PMCID: PMC10382368 DOI: 10.1007/s00335-023-10008-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2023] [Accepted: 06/29/2023] [Indexed: 07/30/2023]
Abstract
Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.
Collapse
Affiliation(s)
- Deborah A Triant
- Division of Animal Sciences, University of Missouri, Columbia, MO, 65211, USA
| | - Amy T Walsh
- Division of Animal Sciences, University of Missouri, Columbia, MO, 65211, USA
| | - Gabrielle A Hartley
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, 06269, USA
| | - Bruna Petry
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Morgan R Stegemiller
- Department of Animal, Veterinary and Food Sciences, University of Idaho, Moscow, ID, 83844, USA
| | - Benjamin M Nelson
- Division of Animal Sciences, University of Missouri, Columbia, MO, 65211, USA
| | - Makenna M McKendrick
- Department of Animal and Food Sciences, Oklahoma State University, Stillwater, OK, 74078, USA
| | - Emily P Fuller
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, 06269, USA
| | - Noelle E Cockett
- Department of Animal, Dairy, and Veterinary Sciences, Utah State University, Logan, UT, 84322, USA
| | - James E Koltes
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Stephanie D McKay
- Department of Animal and Veterinary Sciences, University of Vermont, Burlington, VT, 05405, USA
| | - Jonathan A Green
- Division of Animal Sciences, University of Missouri, Columbia, MO, 65211, USA
| | - Brenda M Murdoch
- Department of Animal, Veterinary and Food Sciences, University of Idaho, Moscow, ID, 83844, USA
| | - Darren E Hagen
- Department of Animal and Food Sciences, Oklahoma State University, Stillwater, OK, 74078, USA
| | - Christine G Elsik
- Division of Animal Sciences, University of Missouri, Columbia, MO, 65211, USA.
- Division of Plant Science & Technology, University of Missouri, Columbia, MO, 65211, USA.
- Institute for Data Science & Informatics, University of Missouri, Columbia, MO, 65211, USA.
| |
Collapse
|
9
|
Bao M, Wang X, Sun R, Wang Z, Li J, Jiang T, Lin A, Wang H, Feng J. Full-Length Transcriptome of the Great Himalayan Leaf-Nosed Bats ( Hipposideros armiger) Optimized Genome Annotation and Revealed the Expression of Novel Genes. Int J Mol Sci 2023; 24:ijms24054937. [PMID: 36902366 PMCID: PMC10003721 DOI: 10.3390/ijms24054937] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Revised: 02/26/2023] [Accepted: 02/28/2023] [Indexed: 03/08/2023] Open
Abstract
The Great Himalayan Leaf-nosed bat (Hipposideros armiger) is one of the most representative species of all echolocating bats and is an ideal model for studying the echolocation system of bats. An incomplete reference genome and limited availability of full-length cDNAs have hindered the identification of alternatively spliced transcripts, which slowed down related basic studies on bats' echolocation and evolution. In this study, we analyzed five organs from H. armiger for the first time using PacBio single-molecule real-time sequencing (SMRT). There were 120 GB of subreads generated, including 1,472,058 full-length non-chimeric (FLNC) sequences. A total of 34,611 alternative splicing (AS) events and 66,010 Alternative Polyadenylation (APA) sites were detected by transcriptome structural analysis. Moreover, a total of 110,611 isoforms were identified, consisting of 52% new isoforms of known genes and 5% of novel gene loci, as well as 2112 novel genes that have not been annotated before in the current reference genome of H. armiger. Furthermore, several key novel genes, including Pol, RAS, NFKB1, and CAMK4, were identified as being associated with nervous, signal transduction, and immune system processes, which may be involved in regulating the auditory nervous perception and immune system that helps bats to regulate in echolocation. In conclusion, the full-length transcriptome results optimized and replenished existing H. armiger genome annotation in multiple ways and offer advantages for newly discovered or previously unrecognized protein-coding genes and isoforms, which can be used as a reference resource.
Collapse
Affiliation(s)
- Mingyue Bao
- College of Life Science, Jilin Agricultural University, Changchun 130118, China
| | - Xue Wang
- College of Life Science, Jilin Agricultural University, Changchun 130118, China
| | - Ruyi Sun
- College of Life Science, Jilin Agricultural University, Changchun 130118, China
| | - Zhiqiang Wang
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun 130117, China
| | - Jiqian Li
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun 130117, China
| | - Tinglei Jiang
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun 130117, China
| | - Aiqing Lin
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun 130117, China
| | - Hui Wang
- College of Life Science, Jilin Agricultural University, Changchun 130118, China
- Correspondence: (H.W.); (J.F.)
| | - Jiang Feng
- College of Life Science, Jilin Agricultural University, Changchun 130118, China
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun 130117, China
- Correspondence: (H.W.); (J.F.)
| |
Collapse
|
10
|
Nelson TM, Ghosh S, Postler TS. L-RAPiT: A Cloud-Based Computing Pipeline for the Analysis of Long-Read RNA Sequencing Data. Int J Mol Sci 2022; 23:ijms232415851. [PMID: 36555493 PMCID: PMC9781625 DOI: 10.3390/ijms232415851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Revised: 12/07/2022] [Accepted: 12/11/2022] [Indexed: 12/15/2022] Open
Abstract
Long-read sequencing (LRS) has been adopted to meet a wide variety of research needs, ranging from the construction of novel transcriptome annotations to the rapid identification of emerging virus variants. Amongst other advantages, LRS preserves more information about RNA at the transcript level than conventional high-throughput sequencing, including far more accurate and quantitative records of splicing patterns. New studies with LRS datasets are being published at an exponential rate, generating a vast reservoir of information that can be leveraged to address a host of different research questions. However, mining such publicly available data in a tailored fashion is currently not easy, as the available software tools typically require familiarity with the command-line interface, which constitutes a significant obstacle to many researchers. Additionally, different research groups utilize different software packages to perform LRS analysis, which often prevents a direct comparison of published results across different studies. To address these challenges, we have developed the Long-Read Analysis Pipeline for Transcriptomics (L-RAPiT), a user-friendly, free pipeline requiring no dedicated computational resources or bioinformatics expertise. L-RAPiT can be implemented directly through Google Colaboratory, a system based on the open-source Jupyter notebook environment, and allows for the direct analysis of transcriptomic reads from Oxford Nanopore and PacBio LRS machines. This new pipeline enables the rapid, convenient, and standardized analysis of publicly available or newly generated LRS datasets.
Collapse
|
11
|
Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing. PLoS One 2022; 17:e0276956. [PMID: 36331983 PMCID: PMC9635732 DOI: 10.1371/journal.pone.0276956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Accepted: 10/18/2022] [Indexed: 11/06/2022] Open
Abstract
The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.
Collapse
|
12
|
Guan D, Halstead MM, Islas-Trejo AD, Goszczynski DE, Cheng HH, Ross PJ, Zhou H. Prediction of transcript isoforms in 19 chicken tissues by Oxford Nanopore long-read sequencing. Front Genet 2022; 13:997460. [PMID: 36246588 PMCID: PMC9561881 DOI: 10.3389/fgene.2022.997460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2022] [Accepted: 08/30/2022] [Indexed: 11/22/2022] Open
Abstract
To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).
Collapse
Affiliation(s)
- Dailu Guan
- Department of Animal Science, University of California Davis, Davis, CA, United States
| | - Michelle M. Halstead
- Department of Animal Science, University of California Davis, Davis, CA, United States
| | - Alma D. Islas-Trejo
- Department of Animal Science, University of California Davis, Davis, CA, United States
| | - Daniel E. Goszczynski
- Department of Animal Science, University of California Davis, Davis, CA, United States
| | - Hans H. Cheng
- USDA, ARS, USNPRC, Avian Disease and Oncology Laboratory, East Lansing, MI, United States
| | - Pablo J. Ross
- Department of Animal Science, University of California Davis, Davis, CA, United States
- *Correspondence: Pablo J. Ross, ; Huaijun Zhou,
| | - Huaijun Zhou
- Department of Animal Science, University of California Davis, Davis, CA, United States
- *Correspondence: Pablo J. Ross, ; Huaijun Zhou,
| |
Collapse
|
13
|
Zhang Y, Lou F, Chen J, Han Z, Yang T, Gao T, Song N. Single-molecule Real-time (SMRT) Sequencing Facilitates Transcriptome Research and Genome Annotation of the Fish Sillago sinica. MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2022; 24:1002-1013. [PMID: 36083383 DOI: 10.1007/s10126-022-10163-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 08/30/2022] [Indexed: 06/15/2023]
Abstract
As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species.
Collapse
Affiliation(s)
- Yuan Zhang
- Fishery College, Ocean University of China, Qingdao, 266003, China
| | - Fangrui Lou
- School of Ocean, Yantai University, Yantai, 264005, China
| | - Jianwei Chen
- BGI-Qingdao, BGI-Shenzhen, Qingdao, 266555, China
| | - Zhiqiang Han
- Fishery College, Zhejiang Ocean University, Zhoushan, 316022, China
| | - Tianyan Yang
- Fishery College, Zhejiang Ocean University, Zhoushan, 316022, China
| | - Tianxiang Gao
- Fishery College, Zhejiang Ocean University, Zhoushan, 316022, China
| | - Na Song
- Fishery College, Ocean University of China, Qingdao, 266003, China.
| |
Collapse
|
14
|
Sage SE, Nicholson P, Leeb T, Gerber V, Jagannathan V. Long-Read Transcriptome of Equine Bronchoalveolar Cells. Genes (Basel) 2022; 13:genes13101722. [PMID: 36292607 PMCID: PMC9602388 DOI: 10.3390/genes13101722] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Revised: 09/09/2022] [Accepted: 09/21/2022] [Indexed: 11/16/2022] Open
Abstract
We used Pacific Biosciences long-read isoform sequencing to generate full-length transcript sequences in equine bronchoalveolar lavage fluid (BALF) cells. Our dataset consisted of 313,563 HiFi reads comprising 805 Mb of polished sequence information. The resulting equine BALF transcriptome consisted of 14,234 full-length transcript isoforms originating from 7017 unique genes. These genes consisted of 6880 previously annotated genes and 137 novel genes. We identified 3428 novel transcripts in addition to 10,806 previously known transcripts. These included transcripts absent from existing genome annotations, transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. We provide transcript-level data for equine BALF cells as a resource to the scientific community.
Collapse
Affiliation(s)
- Sophie Elena Sage
- Swiss Institute of Equine Medicine, Department of Clinical Veterinary Medicine, Vetsuisse Faculty, University of Bern, 3001 Bern, Switzerland
| | - Pamela Nicholson
- Next Generation Sequencing Platform, University of Bern, 3001 Bern, Switzerland
| | - Tosso Leeb
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001 Bern, Switzerland
| | - Vinzenz Gerber
- Swiss Institute of Equine Medicine, Department of Clinical Veterinary Medicine, Vetsuisse Faculty, University of Bern, 3001 Bern, Switzerland
| | - Vidhya Jagannathan
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001 Bern, Switzerland
- Correspondence:
| |
Collapse
|
15
|
Zhang R, Kuo R, Coulter M, Calixto CPG, Entizne JC, Guo W, Marquez Y, Milne L, Riegler S, Matsui A, Tanaka M, Harvey S, Gao Y, Wießner-Kroh T, Paniagua A, Crespi M, Denby K, Hur AB, Huq E, Jantsch M, Jarmolowski A, Koester T, Laubinger S, Li QQ, Gu L, Seki M, Staiger D, Sunkar R, Szweykowska-Kulinska Z, Tu SL, Wachter A, Waugh R, Xiong L, Zhang XN, Conesa A, Reddy ASN, Barta A, Kalyna M, Brown JWS. A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis. Genome Biol 2022; 23:149. [PMID: 35799267 PMCID: PMC9264592 DOI: 10.1186/s13059-022-02711-0] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Accepted: 06/15/2022] [Indexed: 12/15/2022] Open
Abstract
BACKGROUND Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.
Collapse
Affiliation(s)
- Runxuan Zhang
- Information and Computational Sciences, James Hutton Institute, Dundee, DD2 5DA, Scotland, UK.
| | - Richard Kuo
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Midlothian, EH25 9RG, UK
| | - Max Coulter
- Plant Sciences Division, School of Life Sciences, University of Dundee at The James Hutton Institute, Invergowrie, Dundee, DD2 5DA, Scotland, UK
| | - Cristiane P G Calixto
- Plant Sciences Division, School of Life Sciences, University of Dundee at The James Hutton Institute, Invergowrie, Dundee, DD2 5DA, Scotland, UK
- Present address: Institute of Biosciences, University of São Paulo, São Paulo, 05508-090, Brazil
| | - Juan Carlos Entizne
- Plant Sciences Division, School of Life Sciences, University of Dundee at The James Hutton Institute, Invergowrie, Dundee, DD2 5DA, Scotland, UK
| | - Wenbin Guo
- Information and Computational Sciences, James Hutton Institute, Dundee, DD2 5DA, Scotland, UK
| | - Yamile Marquez
- Centre for Genomic Regulation, C/ Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Linda Milne
- Information and Computational Sciences, James Hutton Institute, Dundee, DD2 5DA, Scotland, UK
| | - Stefan Riegler
- Institute of Molecular Plant Biology, Department of Applied Genetics and Cell Biology, University of Natural Resources and Life Sciences (BOKU), Muthgasse 18, 1190, Vienna, Austria
- Present address: Institute of Science and Technology Austria, Am Campus 1, 3400, Klosterneuburg, Austria
| | - Akihiro Matsui
- Plant Genomic Network Research Team, RIKEN Center for Sustainable Resource Science, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan
| | - Maho Tanaka
- Plant Genomic Network Research Team, RIKEN Center for Sustainable Resource Science, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan
| | - Sarah Harvey
- Centre for Novel Agricultural Products (CNAP), Department of Biology, University of York Wentworth Way, York, YO10 5DD, UK
| | - Yubang Gao
- College of Forestry, Fujian Agriculture and Forestry University, Fuzhou, 350002, China
| | - Theresa Wießner-Kroh
- Center for Plant Molecular Biology (ZMBP), University of Tübingen, Auf der Morgenstelle 32, 72076, Tübingen, Germany
| | - Alejandro Paniagua
- Institute for Integrative Systems Biology (CSIC-UV), Spanish National Research Council, Paterna, Valencia, Spain
| | - Martin Crespi
- French National Centre for Scientific Research | CNRS INRAE-Universities of Paris Saclay and Paris, Institute of Plant Sciences Paris Saclay IPS2, Rue de Noetzlin, 91192, Gif sur Yvette, France
| | - Katherine Denby
- Centre for Novel Agricultural Products (CNAP), Department of Biology, University of York Wentworth Way, York, YO10 5DD, UK
| | - Asa Ben Hur
- Department of Computer Science, Colorado State University, 1873 Campus Delivery, Fort Collins, CO, 80523-1873, USA
| | - Enamul Huq
- Department of Molecular Biosciences, University of Texas at Austin, 100 East 24th St., Austin, TX, 78712-1095, USA
| | - Michael Jantsch
- Department of Cell and Developmental Biology, Center for Anatomy and Cell Biology, Medical University of Vienna, Schwarzspanierstrasse 17 A-1090, Vienna, Austria
| | - Artur Jarmolowski
- Department of Gene Expression, Adam Mickiewicz University, Poznań, Poland
| | - Tino Koester
- RNA Biology and Molecular Physiology, Faculty for Biology, Bielefeld University, Universitaetsstrasse 25, 33615, Bielefeld, Germany
| | - Sascha Laubinger
- Institut für Biologie und Umweltwissenschaften (IBU), Carl von Ossietzky Universität Oldenburg, Carl von Ossietzky-Str. 9-11, 26111, Oldenburg, Germany
- Institute of Biology, Department of Genetics, Martin Luther University Halle-Wittenberg, Halle (Saale), Germany
| | - Qingshun Quinn Li
- Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, CA, 91766, USA
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361102, Fujian, China
| | - Lianfeng Gu
- College of Forestry, Fujian Agriculture and Forestry University, Fuzhou, 350002, China
| | - Motoaki Seki
- Plant Genomic Network Research Team, RIKEN Center for Sustainable Resource Science, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan
| | - Dorothee Staiger
- RNA Biology and Molecular Physiology, Faculty for Biology, Bielefeld University, Universitaetsstrasse 25, 33615, Bielefeld, Germany
| | - Ramanjulu Sunkar
- Department of Biochemistry and Molecular Biology, Oklahoma State University, Stillwater, OK, 74078, USA
| | | | - Shih-Long Tu
- Institute of Plant and Microbial Biology, Academia Sinica, Taipei, Taiwan
| | - Andreas Wachter
- Center for Plant Molecular Biology (ZMBP), University of Tübingen, Auf der Morgenstelle 32, 72076, Tübingen, Germany
- Present address: Institute for Molecular Physiology, Johannes Gutenberg University Mainz, Hanns-Dieter-Hüsch-Weg 17, 55128, Mainz, Germany
| | - Robbie Waugh
- Cell and Molecular Sciences, James Hutton Institute, Dundee, DD2 5DA, Scotland, UK
| | - Liming Xiong
- Department of Biology, Hong Kong Baptist University, Hong Kong, China
| | - Xiao-Ning Zhang
- Biology Department, School of Arts and Sciences, St. Bonaventure University, 3261 West State Road, St. Bonaventure, NY, 14778, USA
| | - Ana Conesa
- Institute for Integrative Systems Biology (CSIC-UV), Spanish National Research Council, Paterna, Valencia, Spain
| | - Anireddy S N Reddy
- Department of Biology and Program in Cell and Molecular Biology, Colorado State University, Fort Collins, CO, 80523, USA
| | - Andrea Barta
- Max F. Perutz Laboratories, Medical University of Vienna, Center of Medical Biochemistry, Dr.-Bohr-Gasse 9/3, A-1030, Vienna, Austria
| | - Maria Kalyna
- Institute of Molecular Plant Biology, Department of Applied Genetics and Cell Biology, University of Natural Resources and Life Sciences (BOKU), Muthgasse 18, 1190, Vienna, Austria
| | - John W S Brown
- Plant Sciences Division, School of Life Sciences, University of Dundee at The James Hutton Institute, Invergowrie, Dundee, DD2 5DA, Scotland, UK
- Cell and Molecular Sciences, James Hutton Institute, Dundee, DD2 5DA, Scotland, UK
| |
Collapse
|
16
|
Schon MA, Lutzmayer S, Hofmann F, Nodine MD. Bookend: precise transcript reconstruction with end-guided assembly. Genome Biol 2022; 23:143. [PMID: 35768836 PMCID: PMC9245221 DOI: 10.1186/s13059-022-02700-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2021] [Accepted: 06/05/2022] [Indexed: 12/29/2022] Open
Abstract
We developed Bookend, a package for transcript assembly that incorporates data from different RNA-seq techniques, with a focus on identifying and utilizing RNA 5' and 3' ends. We demonstrate that correct identification of transcript start and end sites is essential for precise full-length transcript assembly. Utilization of end-labeled reads present in full-length single-cell RNA-seq datasets dramatically improves the precision of transcript assembly in single cells. Finally, we show that hybrid assembly across short-read, long-read, and end-capture RNA-seq datasets from Arabidopsis thaliana, as well as meta-assembly of RNA-seq from single mouse embryonic stem cells, can produce reference-quality end-to-end transcript annotations.
Collapse
Affiliation(s)
- Michael A Schon
- Cluster of Plant Developmental Biology, Laboratory of Molecular Biology, Wageningen University & Research, Wageningen, 6708, PB, The Netherlands.
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria.
| | - Stefan Lutzmayer
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria
| | - Falko Hofmann
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria
| | - Michael D Nodine
- Cluster of Plant Developmental Biology, Laboratory of Molecular Biology, Wageningen University & Research, Wageningen, 6708, PB, The Netherlands.
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria.
| |
Collapse
|
17
|
Zhang T, Wang T, Niu Q, Xu L, Chen Y, Gao X, Gao H, Zhang L, Liu GE, Li J, Xu L. Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle. BMC Biol 2022; 20:79. [PMID: 35351103 PMCID: PMC8966188 DOI: 10.1186/s12915-022-01269-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Accepted: 03/03/2022] [Indexed: 11/12/2022] Open
Abstract
BACKGROUND A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle. RESULTS Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits. CONCLUSIONS We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.
Collapse
Affiliation(s)
- Tianliu Zhang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Tianzhen Wang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Qunhao Niu
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Lei Xu
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Yan Chen
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Xue Gao
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Huijiang Gao
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Lupei Zhang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agricultural Research Service, Beltsville, Maryland 20705 USA
| | - Junya Li
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| | - Lingyang Xu
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193 People’s Republic of China
| |
Collapse
|
18
|
Healey HM, Bassham S, Cresko WA. Single-cell Iso-Sequencing enables rapid genome annotation for scRNAseq analysis. Genetics 2022; 220:iyac017. [PMID: 35143654 PMCID: PMC8893252 DOI: 10.1093/genetics/iyac017] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Accepted: 01/24/2022] [Indexed: 12/02/2022] Open
Abstract
Single-cell RNA sequencing is a powerful technique that continues to expand across various biological applications. However, incomplete 3'-UTR annotations can impede single-cell analysis resulting in genes that are partially or completely uncounted. Performing single-cell RNA sequencing with incomplete 3'-UTR annotations can hinder the identification of cell identities and gene expression patterns and lead to erroneous biological inferences. We demonstrate that performing single-cell isoform sequencing in tandem with single-cell RNA sequencing can rapidly improve 3'-UTR annotations. Using threespine stickleback fish (Gasterosteus aculeatus), we show that gene models resulting from a minimal embryonic single-cell isoform sequencing dataset retained 26.1% greater single-cell RNA sequencing reads than gene models from Ensembl alone. Furthermore, pooling our single-cell sequencing isoforms with a previously published adult bulk Iso-Seq dataset from stickleback, and merging the annotation with the Ensembl gene models, resulted in a marginal improvement (+0.8%) over the single-cell isoform sequencing only dataset. In addition, isoforms identified by single-cell isoform sequencing included thousands of new splicing variants. The improved gene models obtained using single-cell isoform sequencing led to successful identification of cell types and increased the reads identified of many genes in our single-cell RNA sequencing stickleback dataset. Our work illuminates single-cell isoform sequencing as a cost-effective and efficient mechanism to rapidly annotate genomes for single-cell RNA sequencing.
Collapse
Affiliation(s)
- Hope M Healey
- Institute of Ecology and Evolution, University of Oregon, Eugene, OR 97403, USA
| | - Susan Bassham
- Institute of Ecology and Evolution, University of Oregon, Eugene, OR 97403, USA
| | - William A Cresko
- Institute of Ecology and Evolution, University of Oregon, Eugene, OR 97403, USA
- Presidential Initiative in Data Science, University of Oregon, Eugene, OR 97403, USA
| |
Collapse
|
19
|
Karimi P, Bakhtiarizadeh MR, Salehi A, Izadnia HR. Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken. Sci Rep 2022; 12:2558. [PMID: 35169237 PMCID: PMC8847365 DOI: 10.1038/s41598-022-06528-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Accepted: 02/02/2022] [Indexed: 12/13/2022] Open
Abstract
Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.
Collapse
Affiliation(s)
- Parastoo Karimi
- Department of Animal and Poultry Science, College of Aburaihan, University of Tehran, Tehran, Iran
| | | | - Abdolreza Salehi
- Department of Animal and Poultry Science, College of Aburaihan, University of Tehran, Tehran, Iran
| | - Hamid Reza Izadnia
- Animal Science Improvement Research Department, Agricultural and Natural Resources Research and Education Center, Safiabad AREEO, Dezful, Iran
| |
Collapse
|
20
|
Hao W, Yang Z, Sun Y, Li J, Zhang D, Liu D, Yang X. Characterization of Alternative Splicing Events in Porcine Skeletal Muscles with Different Intramuscular Fat Contents. Biomolecules 2022; 12:biom12020154. [PMID: 35204660 PMCID: PMC8961525 DOI: 10.3390/biom12020154] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 01/11/2022] [Accepted: 01/14/2022] [Indexed: 02/06/2023] Open
Abstract
Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.
Collapse
Affiliation(s)
- Wanjun Hao
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Zewei Yang
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Yuanlu Sun
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Jiaxin Li
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Dongjie Zhang
- Institute of Animal Husbandry, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China;
| | - Di Liu
- Institute of Animal Husbandry, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China;
- Correspondence: (D.L.); (X.Y.); Tel.: +86-451-8667-7458 (D.L.); +86-451-5519-1738 (X.Y.)
| | - Xiuqin Yang
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
- Correspondence: (D.L.); (X.Y.); Tel.: +86-451-8667-7458 (D.L.); +86-451-5519-1738 (X.Y.)
| |
Collapse
|
21
|
Mueller RC, Ellström P, Howe K, Uliano-Silva M, Kuo RI, Miedzinska K, Warr A, Fedrigo O, Haase B, Mountcastle J, Chow W, Torrance J, Wood JMD, Järhult JD, Naguib MM, Olsen B, Jarvis ED, Smith J, Eöry L, Kraus RHS. A high-quality genome and comparison of short- versus long-read transcriptome of the palaearctic duck Aythya fuligula (tufted duck). Gigascience 2021; 10:giab081. [PMID: 34927191 PMCID: PMC8685854 DOI: 10.1093/gigascience/giab081] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 07/15/2021] [Accepted: 11/22/2021] [Indexed: 11/13/2022] Open
Abstract
BACKGROUND The tufted duck is a non-model organism that experiences high mortality in highly pathogenic avian influenza outbreaks. It belongs to the same bird family (Anatidae) as the mallard, one of the best-studied natural hosts of low-pathogenic avian influenza viruses. Studies in non-model bird species are crucial to disentangle the role of the host response in avian influenza virus infection in the natural reservoir. Such endeavour requires a high-quality genome assembly and transcriptome. FINDINGS This study presents the first high-quality, chromosome-level reference genome assembly of the tufted duck using the Vertebrate Genomes Project pipeline. We sequenced RNA (complementary DNA) from brain, ileum, lung, ovary, spleen, and testis using Illumina short-read and Pacific Biosciences long-read sequencing platforms, which were used for annotation. We found 34 autosomes plus Z and W sex chromosomes in the curated genome assembly, with 99.6% of the sequence assigned to chromosomes. Functional annotation revealed 14,099 protein-coding genes that generate 111,934 transcripts, which implies a mean of 7.9 isoforms per gene. We also identified 246 small RNA families. CONCLUSIONS This annotated genome contributes to continuing research into the host response in avian influenza virus infections in a natural reservoir. Our findings from a comparison between short-read and long-read reference transcriptomics contribute to a deeper understanding of these competing options. In this study, both technologies complemented each other. We expect this annotation to be a foundation for further comparative and evolutionary genomic studies, including many waterfowl relatives with differing susceptibilities to avian influenza viruses.
Collapse
Affiliation(s)
- Ralf C Mueller
- Department of Migration, Max Planck Institute of Animal Behavior, Radolfzell, 78315, Germany
- Department of Biology, University of Konstanz, Konstanz, 78457, Germany
| | - Patrik Ellström
- Department of Medical Sciences, Zoonosis Science Center, Uppsala University, Uppsala, SE-75185, Sweden
| | - Kerstin Howe
- Tree of Life, Wellcome Sanger Institute, Cambridge CB10 1SA, UK
| | | | - Richard I Kuo
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Katarzyna Miedzinska
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Amanda Warr
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Olivier Fedrigo
- Vertebrate Genome Laboratory, The Rockefeller University, New York, 10065, NY
| | - Bettina Haase
- Vertebrate Genome Laboratory, The Rockefeller University, New York, 10065, NY
| | | | - William Chow
- Tree of Life, Wellcome Sanger Institute, Cambridge CB10 1SA, UK
| | - James Torrance
- Tree of Life, Wellcome Sanger Institute, Cambridge CB10 1SA, UK
| | | | - Josef D Järhult
- Department of Medical Sciences, Zoonosis Science Center, Uppsala University, Uppsala, SE-75185, Sweden
| | - Mahmoud M Naguib
- Department of Medical Biochemistry and Microbiology, Zoonosis Science Center, Uppsala University, Uppsala, 75237, Sweden
| | - Björn Olsen
- Department of Medical Sciences, Zoonosis Science Center, Uppsala University, Uppsala, SE-75185, Sweden
| | - Erich D Jarvis
- Vertebrate Genome Laboratory and HHMI, The Rockefeller University, New York, 10065, NY
| | - Jacqueline Smith
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Lél Eöry
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Robert H S Kraus
- Department of Migration, Max Planck Institute of Animal Behavior, Radolfzell, 78315, Germany
- Department of Biology, University of Konstanz, Konstanz, 78457, Germany
| |
Collapse
|
22
|
Leung SK, Jeffries AR, Castanho I, Jordan BT, Moore K, Davies JP, Dempster EL, Bray NJ, O'Neill P, Tseng E, Ahmed Z, Collier DA, Jeffery ED, Prabhakar S, Schalkwyk L, Jops C, Gandal MJ, Sheynkman GM, Hannon E, Mill J. Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing. Cell Rep 2021; 37:110022. [PMID: 34788620 PMCID: PMC8609283 DOI: 10.1016/j.celrep.2021.110022] [Citation(s) in RCA: 59] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Revised: 07/30/2021] [Accepted: 10/28/2021] [Indexed: 12/05/2022] Open
Abstract
Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community. There is widespread transcript diversity in the cortex and many novel transcripts Some genes display big differences in isoform number between human and mouse cortex There is evidence of differential transcript usage between human fetal and adult cortex There are many novel isoforms of genes associated with human brain disease
Collapse
Key Words
- isoform, transcript, expression, brain, cortex, mouse, human, adult, fetal, long-read sequencing, alternative splicing
Collapse
Affiliation(s)
| | | | - Isabel Castanho
- University of Exeter, Exeter, UK; Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Ben T Jordan
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA, USA
| | | | | | | | | | | | | | | | | | - Erin D Jeffery
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA, USA
| | - Shyam Prabhakar
- Genome Institute of Singapore, Agency for Science, Technology and Research (A(∗)STAR), Singapore, Singapore
| | | | - Connor Jops
- Department of Psychiatry and Biobehavioral Sciences, Semel Institute for Neuroscience and Human Behavior, University of California Los Angeles, Los Angeles, CA, USA
| | - Michael J Gandal
- Department of Psychiatry and Biobehavioral Sciences, Semel Institute for Neuroscience and Human Behavior, University of California Los Angeles, Los Angeles, CA, USA
| | - Gloria M Sheynkman
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA, USA; Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA; UVA Cancer Center, University of Virginia, Charlottesville, VA, USA
| | | | | |
Collapse
|
23
|
Lagarrigue S, Lorthiois M, Degalez F, Gilot D, Derrien T. LncRNAs in domesticated animals: from dog to livestock species. Mamm Genome 2021; 33:248-270. [PMID: 34773482 PMCID: PMC9114084 DOI: 10.1007/s00335-021-09928-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Accepted: 10/19/2021] [Indexed: 11/29/2022]
Abstract
Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.
Collapse
Affiliation(s)
| | - Matthias Lorthiois
- Univ Rennes, CNRS, IGDR (Institut de Génétique et Développement de Rennes) - UMR 6290, 2 av Prof Leon Bernard, F-35000, Rennes, France
| | - Fabien Degalez
- INRAE, INSTITUT AGRO, PEGASE UMR 1348, 35590, Saint-Gilles, France
| | - David Gilot
- CLCC Eugène Marquis, INSERM, Université Rennes, UMR_S 1242, 35000, Rennes, France
| | - Thomas Derrien
- Univ Rennes, CNRS, IGDR (Institut de Génétique et Développement de Rennes) - UMR 6290, 2 av Prof Leon Bernard, F-35000, Rennes, France.
| |
Collapse
|
24
|
Naftaly AS, Pau S, White MA. Long-read RNA sequencing reveals widespread sex-specific alternative splicing in threespine stickleback fish. Genome Res 2021; 31:1486-1497. [PMID: 34131005 PMCID: PMC8327910 DOI: 10.1101/gr.274282.120] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Accepted: 06/15/2021] [Indexed: 01/07/2023]
Abstract
Alternate isoforms are important contributors to phenotypic diversity across eukaryotes. Although short-read RNA-sequencing has increased our understanding of isoform diversity, it is challenging to accurately detect full-length transcripts, preventing the identification of many alternate isoforms. Long-read sequencing technologies have made it possible to sequence full-length alternative transcripts, accurately characterizing alternative splicing events, alternate transcription start and end sites, and differences in UTR regions. Here, we use Pacific Biosciences (PacBio) long-read RNA-sequencing (Iso-Seq) to examine the transcriptomes of five organs in threespine stickleback fish (Gasterosteus aculeatus), a widely used genetic model species. The threespine stickleback fish has a refined genome assembly in which gene annotations are based on short-read RNA sequencing and predictions from coding sequence of other species. This suggests some of the existing annotations may be inaccurate or alternative transcripts may not be fully characterized. Using Iso-Seq we detected thousands of novel isoforms, indicating many isoforms are absent in the current Ensembl gene annotations. In addition, we refined many of the existing annotations within the genome. We noted many improperly positioned transcription start sites that were refined with long-read sequencing. The Iso-Seq-predicted transcription start sites were more accurate and verified through ATAC-seq. We also detected many alternative splicing events between sexes and across organs. We found a substantial number of genes in both somatic and gonadal samples that had sex-specific isoforms. Our study highlights the power of long-read sequencing to study the complexity of transcriptomes, greatly improving genomic resources for the threespine stickleback fish.
Collapse
Affiliation(s)
- Alice S Naftaly
- Department of Genetics, University of Georgia, Athens, Georgia 30602, USA
| | - Shana Pau
- Department of Genetics, University of Georgia, Athens, Georgia 30602, USA
- Department of Biology, University of Texas Arlington, Arlington, Texas 76019, USA
| | - Michael A White
- Department of Genetics, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
25
|
Perez CAG, Adachi S, Nong QD, Adhitama N, Matsuura T, Natsume T, Wada T, Kato Y, Watanabe H. Sense-overlapping lncRNA as a decoy of translational repressor protein for dimorphic gene expression. PLoS Genet 2021; 17:e1009683. [PMID: 34319983 PMCID: PMC8351930 DOI: 10.1371/journal.pgen.1009683] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Revised: 08/09/2021] [Accepted: 06/25/2021] [Indexed: 11/18/2022] Open
Abstract
Long noncoding RNAs (lncRNAs) are vastly transcribed and extensively studied but lncRNAs overlapping with the sense orientation of mRNA have been poorly studied. We analyzed the lncRNA DAPALR overlapping with the 5´ UTR of the Doublesex1 (Dsx1), the male determining gene in Daphnia magna. By affinity purification, we identified an RNA binding protein, Shep as a DAPALR binding protein. Shep also binds to Dsx1 5´ UTR by recognizing the overlapping sequence and suppresses translation of the mRNA. In vitro and in vivo analyses indicated that DAPALR increased Dsx1 translation efficiency by sequestration of Shep. This regulation was impaired when the Shep binding site in DAPALR was deleted. These results suggest that Shep suppresses the unintentional translation of Dsx1 by setting a threshold; and when the sense lncRNA DAPALR is expressed, DAPALR cancels the suppression caused by Shep. This mechanism may be important to show dimorphic gene expressions such as sex determination and it may account for the binary expression in various developmental processes. Long noncoding RNAs are vastly transcribed throughout the genome. Among them, RNAs overlapping the protein-coding RNA in sense orientation have been poorly studied because of the difficulty in differentiating their sequences from their overlapping coding RNAs although this class of RNAs has been reported to comprise the majority of the long noncoding RNAs. In the crustacean Daphnia magna, a long noncoding RNA, called DAPALR, is transcribed from the male determining gene, Doublesex1, and overlaps with the Doublesex1 5´ UTR. DAPALR activates Doublesex1 but this regulatory mechanism remains unknown. We found the RNA binding protein Shep bound to the Doublesex1 5´ UTR. In vitro and in vivo experiments indicated that Shep suppresses translation of the mRNA and DAPALR increases Doublesex1 translation efficiency by sequestration of Shep. Since male-specific expression of Doublesex1 is also regulated at the transcriptional level, we propose that Shep cancels the unexpected expression of Doublesex1 and maintains the feminized state for sexual dimorphism but DAPALR suppresses this repression by sequestration of Shep. We infer that this mechanism is not only for binary sex regulation but could function in the binary regulation of other genes in various biological processes.
Collapse
Affiliation(s)
| | - Shungo Adachi
- Cellular and Molecular Biotechnology Research Institute (CMB), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan
| | - Quang Dang Nong
- Department of Biotechnology, Graduate School of Engineering, Osaka University, Suita, Japan
| | - Nikko Adhitama
- Department of Biotechnology, Graduate School of Engineering, Osaka University, Suita, Japan
| | - Tomoaki Matsuura
- Department of Biotechnology, Graduate School of Engineering, Osaka University, Suita, Japan
| | - Toru Natsume
- Cellular and Molecular Biotechnology Research Institute (CMB), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan
| | - Tadashi Wada
- Department of Biotechnology, Graduate School of Engineering, Osaka University, Suita, Japan
| | - Yasuhiko Kato
- Department of Biotechnology, Graduate School of Engineering, Osaka University, Suita, Japan
- * E-mail: (HW); (YK)
| | - Hajime Watanabe
- Department of Biotechnology, Graduate School of Engineering, Osaka University, Suita, Japan
- * E-mail: (HW); (YK)
| |
Collapse
|
26
|
Overbey EG, Ng TT, Catini P, Griggs LM, Stewart P, Tkalcic S, Hawkins RD, Drechsler Y. Transcriptomes of an Array of Chicken Ovary, Intestinal, and Immune Cells and Tissues. Front Genet 2021; 12:664424. [PMID: 34276773 PMCID: PMC8278112 DOI: 10.3389/fgene.2021.664424] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 05/27/2021] [Indexed: 12/16/2022] Open
Abstract
While the chicken (Gallus gallus) is the most consumed agricultural animal worldwide, the chicken transcriptome remains understudied. We have characterized the transcriptome of 10 cell and tissue types from the chicken using RNA-seq, spanning intestinal tissues (ileum, jejunum, proximal cecum), immune cells (B cells, bursa, macrophages, monocytes, spleen T cells, thymus), and reproductive tissue (ovary). We detected 17,872 genes and 24,812 transcripts across all cell and tissue types, representing 73% and 63% of the current gene annotation, respectively. Further quantification of RNA transcript biotypes revealed protein-coding and lncRNAs specific to an individual cell/tissue type. Each cell/tissue type also has an average of around 1.2 isoforms per gene, however, they all have at least one gene with at least 11 isoforms. Differential expression analysis revealed a large number of differentially expressed genes between tissues of the same category (immune and intestinal). Many of these differentially expressed genes in immune cells were involved in cellular processes relating to differentiation and cell metabolism as well as basic functions of immune cells such as cell adhesion and signal transduction. The differential expressed genes of the different segments of the chicken intestine (jejunum, ileum, proximal cecum) correlated to the metabolic processes in nutrient digestion and absorption. These data should provide a valuable resource in understanding the chicken genome.
Collapse
Affiliation(s)
- Eliah G Overbey
- Department of Genome Sciences, Interdepartmental Astrobiology Program, University of Washington, Seattle, WA, United States
| | - Theros T Ng
- College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA, United States
| | - Pietro Catini
- College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA, United States
| | - Lisa M Griggs
- College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA, United States
| | - Paul Stewart
- College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA, United States
| | - Suzana Tkalcic
- College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA, United States
| | - R David Hawkins
- Department of Genome Sciences, Department of Medicine, University of Washington, Seattle, WA, United States
| | - Yvonne Drechsler
- College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA, United States
| |
Collapse
|
27
|
Eldefrawy F, Xu HS, Pusch E, Karkoura A, Alsafy M, Elgendy S, Williams SM, Navara K, Guo TL. Modulation of folliculogenesis in adult laying chickens by bisphenol A and bisphenol S: Perspectives on ovarian morphology and gene expression. Reprod Toxicol 2021; 103:181-190. [PMID: 34147626 DOI: 10.1016/j.reprotox.2021.06.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 05/25/2021] [Accepted: 06/15/2021] [Indexed: 10/21/2022]
Abstract
Both bisphenol A (BPA) and its analog bisphenol S (BPS) are industrial chemicals that have been used to make certain plastic products applied in chicken farms, including food and water containers. They are endocrine disrupting chemicals (EDCs) with xenoestrogenic activities and affect reproductive success in many ways. It was hypothesized that BPA and BPS could adversely affect the folliculogenesis in chickens due to their disruption of the estrogen responses, using either genomic or non-genomic mechanisms. This study investigated the deleterious effects of BPA and BPS on the ovaries when adult layer chickens were orally treated with these EDCs at 50 μg/kg body weight, the reference dose for chronic oral exposure of BPA established by the U.S. EPA. The chickens in both BPA and BPS-treated groups showed a decreased number of the preovulatory follicles. BPA-treated chickens showed a significant decrease in the diameter of F1. Additionally, both BPA and BPS treatments increased the infiltrations of lymphocytes and plasma cells in ovaries. Moreover, it was found that the ovaries of BPS-treated chickens weighed the most among the groups. RNA sequencing and subsequent pathway enrichment analysis of differentially expressed genes revealed that both BPA- and BPS-treatment groups showed significant changes in gene expression and pathways related to reproduction, immune function and carcinogenesis. Taken together, both BPA and BPS are potentially carcinogenic and have deleterious effects on the fertility of laying chickens by inducing inflammation, suggesting that BPS may not be a safe replacement for BPA.
Collapse
Affiliation(s)
- Fatma Eldefrawy
- Department of Anatomy and Embryology, College of Veterinary Medicine, Alexandria University, Egypt; Department of Veterinary Biomedical Sciences, University of Georgia, Athens, GA, United States
| | - Hannah Shibo Xu
- Department of Veterinary Biomedical Sciences, University of Georgia, Athens, GA, United States
| | - Elizabeth Pusch
- Department of Poultry Science, College of Agricultural and Environmental Sciences, University of Georgia, Athens, GA, United States
| | - Ashraf Karkoura
- Department of Anatomy and Embryology, College of Veterinary Medicine, Alexandria University, Egypt
| | - Mohamed Alsafy
- Department of Anatomy and Embryology, College of Veterinary Medicine, Alexandria University, Egypt
| | - Samir Elgendy
- Department of Anatomy and Embryology, College of Veterinary Medicine, Alexandria University, Egypt
| | - Susan M Williams
- Department of Population Health, Poultry Diagnostic and Research Center, University of Georgia, Athens, GA, United States
| | - Kristen Navara
- Department of Poultry Science, College of Agricultural and Environmental Sciences, University of Georgia, Athens, GA, United States
| | - Tai L Guo
- Department of Veterinary Biomedical Sciences, University of Georgia, Athens, GA, United States.
| |
Collapse
|
28
|
Tarazona S, Arzalluz-Luque A, Conesa A. Undisclosed, unmet and neglected challenges in multi-omics studies. NATURE COMPUTATIONAL SCIENCE 2021; 1:395-402. [PMID: 38217236 DOI: 10.1038/s43588-021-00086-z] [Citation(s) in RCA: 54] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 05/17/2021] [Indexed: 01/15/2024]
Abstract
Multi-omics approaches have become a reality in both large genomics projects and small laboratories. However, the multi-omics research community still faces a number of issues that have either not been sufficiently discussed or for which current solutions are still limited. In this Perspective, we elaborate on these limitations and suggest points of attention for future research. We finally discuss new opportunities and challenges brought to the field by the rapid development of single-cell high-throughput molecular technologies.
Collapse
Affiliation(s)
- Sonia Tarazona
- Department of Applied Statistics, Operations Research and Quality, Universitat Politècnica de València, Valencia, Spain
| | - Angeles Arzalluz-Luque
- Department of Applied Statistics, Operations Research and Quality, Universitat Politècnica de València, Valencia, Spain
| | - Ana Conesa
- Microbiology and Cell Science Department, Institute for Food and Agricultural Research, University of Florida, Gainesville, FL, USA.
- Genetics Institute, University of Florida, Gainesville, FL, USA.
- Institute for Integrative Systems Biology, Spanish National Research Council, Valencia, Spain.
| |
Collapse
|
29
|
Halstead MM, Islas-Trejo A, Goszczynski DE, Medrano JF, Zhou H, Ross PJ. Large-Scale Multiplexing Permits Full-Length Transcriptome Annotation of 32 Bovine Tissues From a Single Nanopore Flow Cell. Front Genet 2021; 12:664260. [PMID: 34093657 PMCID: PMC8173071 DOI: 10.3389/fgene.2021.664260] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Accepted: 04/06/2021] [Indexed: 12/18/2022] Open
Abstract
A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.
Collapse
Affiliation(s)
| | | | | | | | | | - Pablo J. Ross
- Department of Animal Science, University of California, Davis, Davis, CA, United States
| |
Collapse
|
30
|
Ren J, Li Q, Zhang Q, Clinton M, Sun C, Yang N. Systematic screening of long intergenic noncoding RNAs expressed during chicken embryogenesis. Poult Sci 2021; 100:101160. [PMID: 34058566 PMCID: PMC8170422 DOI: 10.1016/j.psj.2021.101160] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Revised: 03/05/2021] [Accepted: 03/22/2021] [Indexed: 11/29/2022] Open
Abstract
Long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes, including embryogenesis and development. To provide a systematic analysis of lncRNAs expressed during chicken embryogenesis, we used Iso-Seq and RNA-Seq to identify potential lncRNAs at embryonic stages from d 1 to d 8 of incubation: sequential stages covering gastrulation, somitogenesis, and organogenesis. The data characterized an expanded landscape of lncRNAs, yielding 45,410 distinct lncRNAs (31,282 genes). Amongst these, a set of 13,141 filtered intergenic lncRNAs (lincRNAs) transcribed from 9803 lincRNA gene loci, of which, 66.5% were novel, were further analyzed. These lincRNAs were found to share many characteristics with mammalian lincRNAs, including relatively short lengths, fewer exons, lower expression levels, and stage-specific expression patterns. Functional studies motivated by "guilt-by-association" associated individual lincRNAs with specific GO functions, providing an important resource for future studies of lincRNA function. Most importantly, a weighted gene co-expression network analysis suggested that genes of the brown module were specifically associated with the day 2 stage. LincRNAs within this module were co-expressed with proteins involved in hematopoiesis and lipid metabolism. This study presents the systematic identification of lincRNAs in developing chicken embryos and will serve as a powerful resource for the study of lincRNA functions.
Collapse
Affiliation(s)
- Junxiao Ren
- National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Quanlin Li
- National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Qinghe Zhang
- National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Michael Clinton
- Division of Developmental Biology, The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Roslin, Midlothian EH25 9RG, United Kingdom
| | - Congjiao Sun
- National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Ning Yang
- National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China.
| |
Collapse
|
31
|
Hayrabedyan S, Kostova P, Zlatkov V, Todorova K. Single-cell transcriptomics in the context of long-read nanopore sequencing. BIOTECHNOL BIOTEC EQ 2021. [DOI: 10.1080/13102818.2021.1988868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022] Open
Affiliation(s)
- Soren Hayrabedyan
- Laboratory of Reproductive OMICs Technologies, Institute of Biology and Immunology of Reproduction, Bulgarian Academy of Sciences, Sofia, Bulgaria
| | - Petya Kostova
- Gynecology Clinic, National Oncology Hospital, Sofia, Bulgaria
| | - Viktor Zlatkov
- Department of Obstetrics and Gynecology, Faculty of Medicine, Medical University of Sofia, Sofia, Bulgaria
| | - Krassimira Todorova
- Laboratory of Reproductive OMICs Technologies, Institute of Biology and Immunology of Reproduction, Bulgarian Academy of Sciences, Sofia, Bulgaria
| |
Collapse
|
32
|
Long-read RNA sequencing of human and animal filarial parasites improves gene models and discovers operons. PLoS Negl Trop Dis 2020; 14:e0008869. [PMID: 33196647 PMCID: PMC7704054 DOI: 10.1371/journal.pntd.0008869] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Revised: 11/30/2020] [Accepted: 10/09/2020] [Indexed: 01/01/2023] Open
Abstract
Filarial parasitic nematodes (Filarioidea) cause substantial disease burden to humans and animals around the world. Recently there has been a coordinated global effort to generate, annotate, and curate genomic data from nematode species of medical and veterinary importance. This has resulted in two chromosome-level assemblies (Brugia malayi and Onchocerca volvulus) and 11 additional draft genomes from Filarioidea. These reference assemblies facilitate comparative genomics to explore basic helminth biology and prioritize new drug and vaccine targets. While the continual improvement of genome contiguity and completeness advances these goals, experimental functional annotation of genes is often hindered by poor gene models. Short-read RNA sequencing data and expressed sequence tags, in cooperation with ab initio prediction algorithms, are employed for gene prediction, but these can result in missing clade-specific genes, fragmented models, imperfect mapping of gene ends, and lack of isoform resolution. Long-read RNA sequencing can overcome these drawbacks and greatly improve gene model quality. Here, we present Iso-Seq data for B. malayi and Dirofilaria immitis, etiological agents of lymphatic filariasis and canine heartworm disease, respectively. These data cover approximately half of the known coding genomes and substantially improve gene models by extending untranslated regions, cataloging novel splice junctions from novel isoforms, and correcting mispredicted junctions. Furthermore, we validated computationally predicted operons, manually curated new operons, and merged fragmented gene models. We carried out analyses of poly(A) tails in both species, leading to the identification of non-canonical poly(A) signals. Finally, we prioritized and assessed known and putative anthelmintic targets, correcting or validating gene models for molecular cloning and target-based anthelmintic screening efforts. Overall, these data significantly improve the catalog of gene models for two important parasites, and they demonstrate how long-read RNA sequencing should be prioritized for ongoing improvement of parasitic nematode genome assemblies. Filarial parasitic nematodes are vector-borne parasites that infect humans and animals. Brugia malayi and Dirofilaria immitis are transmitted by mosquitoes and cause human lymphatic filariasis and canine heartworm disease, respectively. Recent years have seen a dramatic increase in genomic and transcriptomic data sets and the concomitant increase in innovative strategies for drug target identification, validation, and screening. However, while the completeness of genome assemblies of filarial parasitic nematodes has seen steady improvements, the reliability of gene models has not kept pace, hindering cloning efforts. Long-read RNA sequencing technologies are uniquely able to improve gene models, but have not been widely used for the causative agents of neglected tropical diseases. Here, we report the improvement of gene models in both B. malayi and D. immitis by long-read RNA sequencing. We identified novel operons, deprecated false positive operons, identified dozens of novel genes, and described the parameters of polyadenylation. We also focused on putative anthelmintic targets, identifying novel isoforms and correcting gene models. These data substantially increase the trustworthiness of gene models in these two species and demonstrate how long-read sequencing approaches should be prioritized in the continued improvement of genome assemblies and their gene annotations.
Collapse
|
33
|
Kuo RI, Cheng Y, Zhang R, Brown JWS, Smith J, Archibald AL, Burt DW. Illuminating the dark side of the human transcriptome with long read transcript sequencing. BMC Genomics 2020; 21:751. [PMID: 33126848 PMCID: PMC7596999 DOI: 10.1186/s12864-020-07123-7] [Citation(s) in RCA: 68] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Accepted: 10/06/2020] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise. RESULTS We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA's genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models. CONCLUSIONS Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.
Collapse
Affiliation(s)
- Richard I Kuo
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Midlothian, EH25 9RG, UK.
| | - Yuanyuan Cheng
- The University of Queensland, St. Lucia, Brisbane, QLD, 4072, Australia
- School of Life and Environmental Sciences, University of Sydney, Sydney, New South Wales, Australia
| | - Runxuan Zhang
- Information and Computational Sciences, The James Hutton Institute, Invergowrie, Dundee, Scotland, UK
| | - John W S Brown
- Plant Sciences Division, School of Life Sciences, University of Dundee, Invergowrie, Dundee, Scotland, UK
- Cell and Molecular Sciences, The James Hutton Institute, Invergowrie, Dundee, Scotland, UK
| | - Jacqueline Smith
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Midlothian, EH25 9RG, UK
| | - Alan L Archibald
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Midlothian, EH25 9RG, UK
| | - David W Burt
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Midlothian, EH25 9RG, UK
- The University of Queensland, St. Lucia, Brisbane, QLD, 4072, Australia
| |
Collapse
|
34
|
Chen D, Du Y, Fan X, Zhu Z, Jiang H, Wang J, Fan Y, Chen H, Zhou D, Xiong C, Zheng Y, Xu X, Luo Q, Guo R. Reconstruction and functional annotation of Ascosphaera apis full-length transcriptome utilizing PacBio long reads combined with Illumina short reads. J Invertebr Pathol 2020; 176:107475. [PMID: 32976816 DOI: 10.1016/j.jip.2020.107475] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2019] [Revised: 08/26/2020] [Accepted: 09/16/2020] [Indexed: 01/19/2023]
Abstract
Ascosphaera apis is a widespread fungal pathogen of honeybee larvae that results in chalkbrood disease, leading to heavy losses for the beekeeping industry in China and many other countries. This work was aimed at generating a full-length transcriptome of A. apis using PacBio single-molecule real-time (SMRT) sequencing. Here, more than 23.97 Gb of clean reads was generated from long-read sequencing of A. apis mycelia, including 464,043 circular consensus sequences (CCS) and 394,142 full-length non-chimeric (FLNC) reads. In total, we identified 174,095 high-confidence transcripts covering 5141 known genes with an average length of 2728 bp. We also discovered 2405 genic loci and 11,623 isoforms that have not been annotated yet within the current reference genome. Additionally, 16,049, 10,682, 4520 and 7253 of the discovered transcripts have annotations in the Non-redundant protein (Nr), Clusters of Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Moreover, 1205 long non-coding RNAs (lncRNAs) were identified, which have less exons, shorter exon and intron lengths, shorter transcript lengths, lower GC percent, lower expression levels, and fewer alternative splicing (AS) evens, compared with protein-coding transcripts. A total of 253 members from 17 transcription factor (TF) families were identified from our transcript datasets. Finally, the expression of A. apis isoforms was validated using a molecular approach. Overall, this is the first report of a full-length transcriptome of entomogenous fungi including A. apis. Our data offer a comprehensive set of reference transcripts and hence contributes to improving the genome annotation and transcriptomic study of A. apis.
Collapse
Affiliation(s)
- Dafu Chen
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Yu Du
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Xiaoxue Fan
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Zhiwei Zhu
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Haibin Jiang
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Jie Wang
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Yuanchan Fan
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Huazhi Chen
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Dingding Zhou
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Cuiling Xiong
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Yanzhen Zheng
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Xijian Xu
- Jiangxi Province Institute of Apiculture, 330201 Nanchang, Jiangxi, China
| | - Qun Luo
- Jiangxi Province Institute of Apiculture, 330201 Nanchang, Jiangxi, China
| | - Rui Guo
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China; Engineering Research Center of Processing and Application of Bee Products of Ministry of Education, Fuzhou 350002, Fujian Province, China.
| |
Collapse
|
35
|
Smith J, Lipkin E, Soller M, Fulton JE, Burt DW. Mapping QTL Associated with Resistance to Avian Oncogenic Marek's Disease Virus (MDV) Reveals Major Candidate Genes and Variants. Genes (Basel) 2020; 11:genes11091019. [PMID: 32872585 PMCID: PMC7564597 DOI: 10.3390/genes11091019] [Citation(s) in RCA: 7] [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: 08/02/2020] [Revised: 08/20/2020] [Accepted: 08/26/2020] [Indexed: 01/13/2023] Open
Abstract
Marek’s disease (MD) represents a significant global economic and animal welfare issue. Marek’s disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify quantitative trait locus regions (QTLR) influencing resistance to MDV, including an F6 population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, microRNAs, long non-coding RNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.
Collapse
Affiliation(s)
- Jacqueline Smith
- The Roslin Institute and Royal (Dick) School of Veterinary Studies R(D)SVS, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Ehud Lipkin
- Department of Genetics, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Morris Soller
- Department of Genetics, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Janet E Fulton
- Hy-Line International, P.O. Box 310, 2583 240th St., Dallas Center, IA 50063, USA
| | - David W Burt
- The Roslin Institute and Royal (Dick) School of Veterinary Studies R(D)SVS, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| |
Collapse
|
36
|
Li Y, Liu Y, Yang H, Zhang T, Naruse K, Tu Q. Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis. Genome Res 2020; 30:924-937. [PMID: 32591361 PMCID: PMC7370878 DOI: 10.1101/gr.258871.119] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 06/17/2020] [Indexed: 12/13/2022]
Abstract
Medaka (Oryzias latipes) has become an important vertebrate model widely used in genetics, developmental biology, environmental sciences, and many other fields. A high-quality genome sequence and a variety of genetic tools are available for this model organism. However, existing genome annotation is still rudimentary, as it was mainly based on computational prediction and short-read RNA-seq data. Here we report a dynamic transcriptome landscape of medaka embryogenesis profiled by long-read RNA-seq, short-read RNA-seq, and ATAC-seq. By integrating these data sets, we constructed a much-improved gene model set including about 17,000 novel isoforms and identified 1600 transcription factors, 1100 long noncoding RNAs, and 150,000 potential cis-regulatory elements as well. Time-series data sets provided another dimension of information. With the expression dynamics of genes and accessibility dynamics of cis-regulatory elements, we investigated isoform switching, as well as regulatory logic between accessible elements and genes, during embryogenesis. We built a user-friendly medaka omics data portal to present these data sets. This resource provides the first comprehensive omics data sets of medaka embryogenesis. Ultimately, we term these three assays as the minimum ENCODE toolbox and propose the use of it as the initial and essential profiling genomic assays for model organisms that have limited data available. This work will be of great value for the research community using medaka as the model organism and many others as well.
Collapse
Affiliation(s)
- Yingshu Li
- State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Innovation Academy for Seed Design, Chinese Academy of Sciences, Beijing 100101, China.,Key Laboratory of Genetic Network Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yongjie Liu
- State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Innovation Academy for Seed Design, Chinese Academy of Sciences, Beijing 100101, China.,Key Laboratory of Genetic Network Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hang Yang
- State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Innovation Academy for Seed Design, Chinese Academy of Sciences, Beijing 100101, China.,Key Laboratory of Genetic Network Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Ting Zhang
- State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Innovation Academy for Seed Design, Chinese Academy of Sciences, Beijing 100101, China.,Key Laboratory of Genetic Network Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
| | - Kiyoshi Naruse
- Laboratory of Bioresources, National Institute for Basic Biology, Okazaki 444-8585, Aichi, Japan
| | - Qiang Tu
- State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Innovation Academy for Seed Design, Chinese Academy of Sciences, Beijing 100101, China.,Key Laboratory of Genetic Network Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
37
|
Chen T, Sun Q, Ma Y, Zeng W, Liu R, Qu D, Huang L, Xu H. A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing. Mol Genet Genomics 2020; 295:1227-1237. [PMID: 32524299 DOI: 10.1007/s00438-020-01691-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2019] [Accepted: 05/25/2020] [Indexed: 02/05/2023]
Abstract
The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171 bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.
Collapse
Affiliation(s)
- Tao Chen
- College of Biotechnology, Jiangsu University of Science and Technology, Zhenjiang, 212003, Jiangsu, China
- The Sericultural Research Institute, Chinese Academy of Agricultural Sciences, Nanjing, 212018, Jiangsu, China
| | - Qiwei Sun
- International Bioinformatics Center, BGI Genomics Co., Ltd, Shenzhen, 518083, Guangdong, China
| | - Yan Ma
- State Key Laboratory of Silkworm Genome Biology, College of Biotechnology, Southwest University, Chongqing, 400715, China
| | - Wenhui Zeng
- State Key Laboratory of Silkworm Genome Biology, College of Biotechnology, Southwest University, Chongqing, 400715, China
| | - Rongpeng Liu
- State Key Laboratory of Silkworm Genome Biology, College of Biotechnology, Southwest University, Chongqing, 400715, China
| | - Dawei Qu
- State Key Laboratory of Silkworm Genome Biology, College of Biotechnology, Southwest University, Chongqing, 400715, China
| | - Lihua Huang
- International Bioinformatics Center, BGI Genomics Co., Ltd, Shenzhen, 518083, Guangdong, China
| | - Hanfu Xu
- State Key Laboratory of Silkworm Genome Biology, College of Biotechnology, Southwest University, Chongqing, 400715, China.
| |
Collapse
|
38
|
de la Fuente L, Arzalluz-Luque Á, Tardáguila M, Del Risco H, Martí C, Tarazona S, Salguero P, Scott R, Lerma A, Alastrue-Agudo A, Bonilla P, Newman JRB, Kosugi S, McIntyre LM, Moreno-Manzano V, Conesa A. tappAS: a comprehensive computational framework for the analysis of the functional impact of differential splicing. Genome Biol 2020; 21:119. [PMID: 32423416 PMCID: PMC7236505 DOI: 10.1186/s13059-020-02028-w] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 04/23/2020] [Indexed: 12/26/2022] Open
Abstract
Recent advances in long-read sequencing solve inaccuracies in alternative transcript identification of full-length transcripts in short-read RNA-Seq data, which encourages the development of methods for isoform-centered functional analysis. Here, we present tappAS, the first framework to enable a comprehensive Functional Iso-Transcriptomics (FIT) analysis, which is effective at revealing the functional impact of context-specific post-transcriptional regulation. tappAS uses isoform-resolved annotation of coding and non-coding functional domains, motifs, and sites, in combination with novel analysis methods to interrogate different aspects of the functional readout of transcript variants and isoform regulation. tappAS software and documentation are available at https://app.tappas.org.
Collapse
Affiliation(s)
- Lorena de la Fuente
- Genomics of Gene Expression Laboratory, Prince Felipe Research Center, Valencia, Spain
- Present Address: Bioinformatics Unit, IIS Fundación Jiménez Díaz, Madrid, Spain
| | - Ángeles Arzalluz-Luque
- Department of Statistics and Operational Research, Polytechnical University of Valencia, Valencia, Spain
| | - Manuel Tardáguila
- Department of Microbiology and Cell Science, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL, USA
- Present Address: Human Genetics Department, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
| | - Héctor Del Risco
- Department of Microbiology and Cell Science, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL, USA
| | - Cristina Martí
- Genomics of Gene Expression Laboratory, Prince Felipe Research Center, Valencia, Spain
| | - Sonia Tarazona
- Department of Statistics and Operational Research, Polytechnical University of Valencia, Valencia, Spain
| | - Pedro Salguero
- Genomics of Gene Expression Laboratory, Prince Felipe Research Center, Valencia, Spain
| | - Raymond Scott
- Department of Microbiology and Cell Science, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL, USA
| | - Alberto Lerma
- Genomics of Gene Expression Laboratory, Prince Felipe Research Center, Valencia, Spain
| | - Ana Alastrue-Agudo
- Present Address: Human Genetics Department, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
| | - Pablo Bonilla
- Present Address: Human Genetics Department, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
| | - Jeremy R B Newman
- Genetics Institute, University of Florida, Gainesville, FL, USA
- Department of Pathology, University of Florida, Gainesville, FL, USA
| | - Shunichi Kosugi
- Genetics Institute, University of Florida, Gainesville, FL, USA
- Laboratory for Statistical and Translational Genetics, Center for Integrative Medical Sciences, RIKEN, Wako, Japan
| | - Lauren M McIntyre
- Genetics Institute, University of Florida, Gainesville, FL, USA
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL, USA
| | | | - Ana Conesa
- Department of Microbiology and Cell Science, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL, USA.
- Genetics Institute, University of Florida, Gainesville, FL, USA.
| |
Collapse
|
39
|
Abstract
Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.
Collapse
|
40
|
Kalenik BM, Góra-Sochacka A, Stachyra A, Olszewska-Tomczyk M, Fogtman A, Sawicka R, Śmietanka K, Sirko A. Response to a DNA vaccine against the H5N1 virus depending on the chicken line and number of doses. Virol J 2020; 17:66. [PMID: 32381003 PMCID: PMC7206725 DOI: 10.1186/s12985-020-01335-9] [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: 02/17/2020] [Accepted: 04/23/2020] [Indexed: 01/16/2023] Open
Abstract
Background Avian influenza virus infections cause significant economic losses on poultry farms and pose the threat of a possible pandemic outbreak. Routine vaccination of poultry against avian influenza is not recommended in Europe, however it has been ordered in some other countries, and more countries are considering use of the avian influenza vaccine as a component of their control strategy. Although a variety of such vaccines have been tested, most research has concentrated on specific antibodies and challenge experiments. Methods We monitored the transcriptomic response to a DNA vaccine encoding hemagglutinin from the highly pathogenic H5N1 avian influenza virus in the spleens of broiler and layer chickens. Moreover, in layer chickens the response to one and two doses of the vaccine was compared. Results All groups of birds immunized with two doses of the vaccine responded at the humoral level by producing specific anti-hemagglutinin antibodies. A response to the vaccine was also detected in the spleen transcriptomes. Differential expression of many genes encoding noncoding RNA and proteins functionally connected to the neuroendocrine-immune system was observed in different immunized groups. Conclusion Broiler chickens showed a higher number and wider range of fold-changes in the transcriptional response than laying hens.
Collapse
Affiliation(s)
- Barbara Małgorzata Kalenik
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5A, 02-106, Warsaw, Poland
| | - Anna Góra-Sochacka
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5A, 02-106, Warsaw, Poland
| | - Anna Stachyra
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5A, 02-106, Warsaw, Poland
| | - Monika Olszewska-Tomczyk
- Department of Poultry Diseases, National Veterinary Research Institute, Al. Partyzantow 57, 24-100, Puławy, Poland
| | - Anna Fogtman
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5A, 02-106, Warsaw, Poland
| | - Róża Sawicka
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5A, 02-106, Warsaw, Poland
| | - Krzysztof Śmietanka
- Department of Poultry Diseases, National Veterinary Research Institute, Al. Partyzantow 57, 24-100, Puławy, Poland
| | - Agnieszka Sirko
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5A, 02-106, Warsaw, Poland.
| |
Collapse
|
41
|
Sahlin K, Medvedev P. De Novo Clustering of Long-Read Transcriptome Data Using a Greedy, Quality Value-Based Algorithm. J Comput Biol 2020; 27:472-484. [PMID: 32181688 DOI: 10.1089/cmb.2019.0299] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Long-read sequencing of transcripts with Pacific Biosciences (PacBio) Iso-Seq and Oxford Nanopore Technologies has proven to be central to the study of complex isoform landscapes in many organisms. However, current de novo transcript reconstruction algorithms from long-read data are limited, leaving the potential of these technologies unfulfilled. A common bottleneck is the dearth of scalable and accurate algorithms for clustering long reads according to their gene family of origin. To address this challenge, we develop isONclust, a clustering algorithm that is greedy (to scale) and makes use of quality values (to handle variable error rates). We test isONclust on three simulated and five biological data sets, across a breadth of organisms, technologies, and read depths. Our results demonstrate that isONclust is a substantial improvement over previous approaches, both in terms of overall accuracy and/or scalability to large data sets.
Collapse
Affiliation(s)
- Kristoffer Sahlin
- Department of Computer Science and Engineering and Pennsylvania State University, University Park, Pennsylvania
| | - Paul Medvedev
- Department of Computer Science and Engineering and Pennsylvania State University, University Park, Pennsylvania.,Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania.,Center for Computational Biology and Bioinformatics, Pennsylvania State University, University Park, Pennsylvania
| |
Collapse
|
42
|
Ducrest A, Neuenschwander S, Schmid‐Siegert E, Pagni M, Train C, Dylus D, Nevers Y, Warwick Vesztrocy A, San‐Jose LM, Dupasquier M, Dessimoz C, Xenarios I, Roulin A, Goudet J. New genome assembly of the barn owl ( Tyto alba alba). Ecol Evol 2020; 10:2284-2298. [PMID: 32184981 PMCID: PMC7069322 DOI: 10.1002/ece3.5991] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Revised: 12/05/2019] [Accepted: 12/16/2019] [Indexed: 12/25/2022] Open
Abstract
New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single-Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.
Collapse
Affiliation(s)
- Anne‐Lyse Ducrest
- Department of Ecology and EvolutionUniversity of LausanneLausanneSwitzerland
| | | | | | - Marco Pagni
- Vital‐ITSwiss Institute of BioinformaticsLausanneSwitzerland
| | - Clément Train
- Department of Computational BiologyUniversity of LausanneLausanneSwitzerland
- Center for Integrative GenomicsUniversity of LausanneLausanneSwitzerland
- Swiss Institute of BioinformaticsLausanneSwitzerland
| | - David Dylus
- Department of Computational BiologyUniversity of LausanneLausanneSwitzerland
- Center for Integrative GenomicsUniversity of LausanneLausanneSwitzerland
- Swiss Institute of BioinformaticsLausanneSwitzerland
| | - Yannis Nevers
- Department of Computational BiologyUniversity of LausanneLausanneSwitzerland
- Center for Integrative GenomicsUniversity of LausanneLausanneSwitzerland
- Swiss Institute of BioinformaticsLausanneSwitzerland
| | - Alex Warwick Vesztrocy
- Center for Life's Origins and EvolutionDepartment of Genetics, Evolution and EnvironmentUniversity College LondonLondonUK
| | - Luis M. San‐Jose
- Laboratory Evolution and Biological DiversityUMR 5174CNRSUniversity of Toulouse III Paul SabatierToulouseFrance
| | | | - Christophe Dessimoz
- Department of Computational BiologyUniversity of LausanneLausanneSwitzerland
- Center for Integrative GenomicsUniversity of LausanneLausanneSwitzerland
- Swiss Institute of BioinformaticsLausanneSwitzerland
| | - Ioannis Xenarios
- Center for Integrative GenomicsUniversity of LausanneLausanneSwitzerland
| | - Alexandre Roulin
- Department of Ecology and EvolutionUniversity of LausanneLausanneSwitzerland
| | - Jérôme Goudet
- Department of Ecology and EvolutionUniversity of LausanneLausanneSwitzerland
- Swiss Institute of BioinformaticsLausanneSwitzerland
| |
Collapse
|
43
|
Luo Y, Liao X, Wu FX, Wang J. Computational Approaches for Transcriptome Assembly Based on Sequencing Technologies. Curr Bioinform 2020. [DOI: 10.2174/1574893614666190410155603] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Transcriptome assembly plays a critical role in studying biological properties and
examining the expression levels of genomes in specific cells. It is also the basis of many
downstream analyses. With the increase of speed and the decrease in cost, massive sequencing
data continues to accumulate. A large number of assembly strategies based on different
computational methods and experiments have been developed. How to efficiently perform
transcriptome assembly with high sensitivity and accuracy becomes a key issue. In this work, the
issues with transcriptome assembly are explored based on different sequencing technologies.
Specifically, transcriptome assemblies with next-generation sequencing reads are divided into
reference-based assemblies and de novo assemblies. The examples of different species are used to
illustrate that long reads produced by the third-generation sequencing technologies can cover fulllength
transcripts without assemblies. In addition, different transcriptome assemblies using the
Hybrid-seq methods and other tools are also summarized. Finally, we discuss the future directions
of transcriptome assemblies.
Collapse
Affiliation(s)
- Yuwen Luo
- School of Computer Science and Engineering, Central South University, Changsha, China
| | - Xingyu Liao
- School of Computer Science and Engineering, Central South University, Changsha, China
| | - Fang-Xiang Wu
- Division of Biomedical Engineering, University of Saskatchewan, Saskatchewan, Canada
| | - Jianxin Wang
- School of Computer Science and Engineering, Central South University, Changsha, China
| |
Collapse
|
44
|
Weirick T, Militello G, Hosen MR, John D, Moore JB, Uchida S. Investigation of RNA Editing Sites within Bound Regions of RNA-Binding Proteins. High Throughput 2019; 8:ht8040019. [PMID: 31795425 PMCID: PMC6970233 DOI: 10.3390/ht8040019] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Revised: 11/08/2019] [Accepted: 11/27/2019] [Indexed: 12/16/2022] Open
Abstract
Studies in epitranscriptomics indicate that RNA is modified by a variety of enzymes. Among these RNA modifications, adenosine to inosine (A-to-I) RNA editing occurs frequently in the mammalian transcriptome. These RNA editing sites can be detected directly from RNA sequencing (RNA-seq) data by examining nucleotide changes from adenosine (A) to guanine (G), which substitutes for inosine (I). However, a careful investigation of such nucleotide changes must be conducted to distinguish sequencing errors and genomic mutations from the genuine editing sites. Building upon our recent introduction of an easy-to-use bioinformatics tool, RNA Editor, to detect RNA editing events from RNA-seq data, we examined the extent by which RNA editing events affect the binding of RNA-binding proteins (RBP). Through employing bioinformatic techniques, we uncovered that RNA editing sites occur frequently in RBP-bound regions. Moreover, the presence of RNA editing sites are more frequent when RNA editing islands were examined, which are regions in which RNA editing sites are present in clusters. When the binding of one RBP, human antigen R [HuR; encoded by ELAV-like protein 1 (ELAV1)], was quantified experimentally, its binding was reduced upon silencing of the RNA editing enzyme adenosine deaminases acting on RNA (ADAR) compared to the control-suggesting that the presence of RNA editing islands influence HuR binding to its target regions. These data indicate RNA editing as an important mediator of RBP-RNA interactions-a mechanism which likely constitutes an additional mode of post-transcription gene regulation in biological systems.
Collapse
Affiliation(s)
- Tyler Weirick
- Cardiovascular Innovation Institute, University of Louisville, Louisville, KY 40202, USA
- RIKEN Center for Integrative Medical Sciences (IMS), 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Giuseppe Militello
- Cardiovascular Innovation Institute, University of Louisville, Louisville, KY 40202, USA
- Department of Molecular Cellular and Developmental Biology, Yale University, Yale Science Building-260 Whitney Avenue, New Haven, CT 06511, USA;
| | - Mohammed Rabiul Hosen
- Department of Internal Medicine-II, Molecular Cardiology, Biomedical Center (BMZ), University of Bonn, Sigmund-Freud-Str. 25, Bonn 53127, Germany;
| | - David John
- Institute of Cardiovascular Regeneration, Centre for Molecular Medicine, Goethe University Frankfurt, Theodor-Stern-Kai 7, Frankfurt am Main 60590, Germany;
| | - Joseph B. Moore
- The Christina Lee Brown Envirome Institute, Department of Medicine, University of Louisville, Louisville, KY 40202, USA;
- Diabetes and Obesity Center, University of Louisville, Louisville, KY 40202, USA
| | - Shizuka Uchida
- Cardiovascular Innovation Institute, University of Louisville, Louisville, KY 40202, USA
- The Christina Lee Brown Envirome Institute, Department of Medicine, University of Louisville, Louisville, KY 40202, USA;
- Diabetes and Obesity Center, University of Louisville, Louisville, KY 40202, USA
- Correspondence: ; Tel.: +1-502-854-0570
| |
Collapse
|
45
|
Sessegolo C, Cruaud C, Da Silva C, Cologne A, Dubarry M, Derrien T, Lacroix V, Aury JM. Transcriptome profiling of mouse samples using nanopore sequencing of cDNA and RNA molecules. Sci Rep 2019; 9:14908. [PMID: 31624302 PMCID: PMC6797730 DOI: 10.1038/s41598-019-51470-9] [Citation(s) in RCA: 56] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Accepted: 09/28/2019] [Indexed: 01/27/2023] Open
Abstract
Our vision of DNA transcription and splicing has changed dramatically with the introduction of short-read sequencing. These high-throughput sequencing technologies promised to unravel the complexity of any transcriptome. Generally gene expression levels are well-captured using these technologies, but there are still remaining caveats due to the limited read length and the fact that RNA molecules had to be reverse transcribed before sequencing. Oxford Nanopore Technologies has recently launched a portable sequencer which offers the possibility of sequencing long reads and most importantly RNA molecules. Here we generated a full mouse transcriptome from brain and liver using the Oxford Nanopore device. As a comparison, we sequenced RNA (RNA-Seq) and cDNA (cDNA-Seq) molecules using both long and short reads technologies and tested the TeloPrime preparation kit, dedicated to the enrichment of full-length transcripts. Using spike-in data, we confirmed that expression levels are efficiently captured by cDNA-Seq using short reads. More importantly, Oxford Nanopore RNA-Seq tends to be more efficient, while cDNA-Seq appears to be more biased. We further show that the cDNA library preparation of the Nanopore protocol induces read truncation for transcripts containing internal runs of T's. This bias is marked for runs of at least 15 T's, but is already detectable for runs of at least 9 T's and therefore concerns more than 20% of expressed transcripts in mouse brain and liver. Finally, we outline that bioinformatics challenges remain ahead for quantifying at the transcript level, especially when reads are not full-length. Accurate quantification of repeat-associated genes such as processed pseudogenes also remains difficult, and we show that current mapping protocols which map reads to the genome largely over-estimate their expression, at the expense of their parent gene.
Collapse
Affiliation(s)
- Camille Sessegolo
- Univ Lyon, Université Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Évolutive UMR5558, F-69622, Villeurbanne, France
- EPI ERABLE - Inria Grenoble, Rhône-Alpes, France
| | - Corinne Cruaud
- Genoscope, Institut de biologie François-Jacob, Commissariat a l'Energie Atomique (CEA), Université Paris-Saclay, F-91057, Evry, France
| | - Corinne Da Silva
- Genoscope, Institut de biologie François-Jacob, Commissariat a l'Energie Atomique (CEA), Université Paris-Saclay, F-91057, Evry, France
| | - Audric Cologne
- Univ Lyon, Université Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Évolutive UMR5558, F-69622, Villeurbanne, France
- EPI ERABLE - Inria Grenoble, Rhône-Alpes, France
| | - Marion Dubarry
- Genoscope, Institut de biologie François-Jacob, Commissariat a l'Energie Atomique (CEA), Université Paris-Saclay, F-91057, Evry, France
| | - Thomas Derrien
- Univ Rennes, CNRS, IGDR (Institut de génétique et développement de Rennes) - UMR 6290, F-35000, Rennes, France
| | - Vincent Lacroix
- Univ Lyon, Université Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Évolutive UMR5558, F-69622, Villeurbanne, France
- EPI ERABLE - Inria Grenoble, Rhône-Alpes, France
| | - Jean-Marc Aury
- Genoscope, Institut de biologie François-Jacob, Commissariat a l'Energie Atomique (CEA), Université Paris-Saclay, F-91057, Evry, France.
| |
Collapse
|
46
|
Beauclair L, Ramé C, Arensburger P, Piégu B, Guillou F, Dupont J, Bigot Y. Sequence properties of certain GC rich avian genes, their origins and absence from genome assemblies: case studies. BMC Genomics 2019; 20:734. [PMID: 31610792 PMCID: PMC6792250 DOI: 10.1186/s12864-019-6131-1] [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/11/2019] [Accepted: 09/23/2019] [Indexed: 12/14/2022] Open
Abstract
Background More and more eukaryotic genomes are sequenced and assembled, most of them presented as a complete model in which missing chromosomal regions are filled by Ns and where a few chromosomes may be lacking. Avian genomes often contain sequences with high GC content, which has been hypothesized to be at the origin of many missing sequences in these genomes. We investigated features of these missing sequences to discover why some may not have been integrated into genomic libraries and/or sequenced. Results The sequences of five red jungle fowl cDNA models with high GC content were used as queries to search publicly available datasets of Illumina and Pacbio sequencing reads. These were used to reconstruct the leptin, TNFα, MRPL52, PCP2 and PET100 genes, all of which are absent from the red jungle fowl genome model. These gene sequences displayed elevated GC contents, had intron sizes that were sometimes larger than non-avian orthologues, and had non-coding regions that contained numerous tandem and inverted repeat sequences with motifs able to assemble into stable G-quadruplexes and intrastrand dyadic structures. Our results suggest that Illumina technology was unable to sequence the non-coding regions of these genes. On the other hand, PacBio technology was able to sequence these regions, but with dramatically lower efficiency than would typically be expected. Conclusions High GC content was not the principal reason why numerous GC-rich regions of avian genomes are missing from genome assembly models. Instead, it is the presence of tandem repeats containing motifs capable of assembling into very stable secondary structures that is likely responsible.
Collapse
Affiliation(s)
- Linda Beauclair
- PRC, UMR INRA0085, CNRS 7247, Centre INRA Val de Loire, 37380, Nouzilly, France
| | - Christelle Ramé
- PRC, UMR INRA0085, CNRS 7247, Centre INRA Val de Loire, 37380, Nouzilly, France
| | - Peter Arensburger
- Biological Sciences Department, California State Polytechnic University, Pomona, CA, 91768, USA
| | - Benoît Piégu
- PRC, UMR INRA0085, CNRS 7247, Centre INRA Val de Loire, 37380, Nouzilly, France
| | - Florian Guillou
- PRC, UMR INRA0085, CNRS 7247, Centre INRA Val de Loire, 37380, Nouzilly, France
| | - Joëlle Dupont
- PRC, UMR INRA0085, CNRS 7247, Centre INRA Val de Loire, 37380, Nouzilly, France
| | - Yves Bigot
- PRC, UMR INRA0085, CNRS 7247, Centre INRA Val de Loire, 37380, Nouzilly, France.
| |
Collapse
|
47
|
Ren J, Sun C, Clinton M, Yang N. Dynamic Transcriptional Landscape of the Early Chick Embryo. Front Cell Dev Biol 2019; 7:196. [PMID: 31572727 PMCID: PMC6751280 DOI: 10.3389/fcell.2019.00196] [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/04/2019] [Accepted: 08/29/2019] [Indexed: 01/22/2023] Open
Abstract
Defining the dynamic transcriptome of the early embryo at high resolution would assist greatly in understanding vertebrate development. Here, we describe the dynamic transcription landscape of early chick embryo development using advanced single-molecule long-read isoform sequencing (Iso-Seq) and RNA-Seq technology. Our transcriptomic profiling reflected the time course of chicken embryonic development from day 1 to day 8 of incubation, a period encompassing gastrulation, somitogenesis, and organogenesis. This analysis identified transcriptional isoforms, alternative splicing (AS) events, fusion transcripts, alternative polyadenylation (APA) sites, and novel genes. Our results showed that intron retention (IR) represented the most abundant AS type and displayed distinct features and dynamic modulation during development. Moreover, we constructed a high-resolution expression profile across embryonic development. Our combined expression dataset correlates distinct gene clusters with specific morphological changes, and provides the first framework for the molecular basis of early chicken embryogenesis. Analysis of gene expression in the developing chicken embryo highlighted the dynamic nature and complexity of the chicken transcriptome and demonstrated that dramatically increased IR events are associated with distinct gene sets.
Collapse
Affiliation(s)
- Junxiao Ren
- Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China.,National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, China Agricultural University, Beijing, China
| | - Congjiao Sun
- Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China.,National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, China Agricultural University, Beijing, China
| | - Michael Clinton
- Division of Developmental Biology, The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Midlothian, United Kingdom
| | - Ning Yang
- Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China.,National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, China Agricultural University, Beijing, China
| |
Collapse
|
48
|
Wang R, Nambiar R, Zheng D, Tian B. PolyA_DB 3 catalogs cleavage and polyadenylation sites identified by deep sequencing in multiple genomes. Nucleic Acids Res 2019; 46:D315-D319. [PMID: 29069441 PMCID: PMC5753232 DOI: 10.1093/nar/gkx1000] [Citation(s) in RCA: 143] [Impact Index Per Article: 28.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2017] [Accepted: 10/12/2017] [Indexed: 12/11/2022] Open
Abstract
PolyA_DB is a database cataloging cleavage and polyadenylation sites (PASs) in several genomes. Previous versions were based mainly on expressed sequence tags (ESTs), which had a limited amount and could lead to inaccurate PAS identification due to the presence of internal A-rich sequences in transcripts. Here, we present an updated version of the database based solely on deep sequencing data. First, PASs are mapped by the 3′ region extraction and deep sequencing (3′READS) method, ensuring unequivocal PAS identification. Second, a large volume of data based on diverse biological samples increases PAS coverage by 3.5-fold over the EST-based version and provides PAS usage information. Third, strand-specific RNA-seq data are used to extend annotated 3′ ends of genes to obtain more thorough annotations of alternative polyadenylation (APA) sites. Fourth, conservation information of PAS across mammals sheds light on significance of APA sites. The database (URL: http://www.polya-db.org/v3) currently holds PASs in human, mouse, rat and chicken, and has links to the UCSC genome browser for further visualization and for integration with other genomic data.
Collapse
Affiliation(s)
- Ruijia Wang
- Department of Microbiology, Biochemistry and Molecular Genetics, Rutgers New Jersey Medical School and Rutgers Cancer Institute of New Jersey, Newark, NJ 07103, USA
| | - Ram Nambiar
- Department of Computer Science, New Jersey Institute of Technology, Newark, NJ 07102, USA
| | - Dinghai Zheng
- Department of Microbiology, Biochemistry and Molecular Genetics, Rutgers New Jersey Medical School and Rutgers Cancer Institute of New Jersey, Newark, NJ 07103, USA
| | - Bin Tian
- Department of Microbiology, Biochemistry and Molecular Genetics, Rutgers New Jersey Medical School and Rutgers Cancer Institute of New Jersey, Newark, NJ 07103, USA
| |
Collapse
|
49
|
A Premature Stop Codon in RAF1 Is the Priority Candidate Causative Mutation of the Inherited Chicken Wingless-2 Developmental Syndrome. Genes (Basel) 2019; 10:genes10050353. [PMID: 31075853 PMCID: PMC6562611 DOI: 10.3390/genes10050353] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Revised: 04/23/2019] [Accepted: 04/26/2019] [Indexed: 11/17/2022] Open
Abstract
The chicken wingless-2 (wg-2) mutation is inherited in an autosomal recessive fashion, and the resulting phenotype in mutant (wg-2/wg-2) individuals is a developmental syndrome characterized by absent wings, truncated legs, craniofacial as well as skin and feather defects, and kidney malformations. Mapping and genotyping established that the mutation resides within 227 kilobases (kb) of chromosome 12 in a wg-2 congenic inbred line. A capture array was designed to target and sequence the candidate region along with flanking DNA in 24 birds from the line. Many point mutations and insertions or deletions were identified, and analysis of the linked variants indicated a point mutation predicted to cause a premature stop codon in the RAF1 gene. Expression studies were conducted inclusive of all genes in the candidate region. Interestingly, RAF1 transcription was elevated, yet the protein was absent in the mutants relative to normal individuals. RAF1 encodes a protein integral to the Ras/Raf/MAPK signaling pathway controlling cellular proliferation, and notably, human RASopathies are developmental syndromes caused by germline mutations in genes of this pathway. Our work indicates RAF1 as the priority candidate causative gene for wg-2 and provides a new animal model to study an important signaling pathway implicated in limb development, as well as RASopathies.
Collapse
|
50
|
Beiki H, Liu H, Huang J, Manchanda N, Nonneman D, Smith TPL, Reecy JM, Tuggle CK. Improved annotation of the domestic pig genome through integration of Iso-Seq and RNA-seq data. BMC Genomics 2019; 20:344. [PMID: 31064321 PMCID: PMC6505119 DOI: 10.1186/s12864-019-5709-y] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 04/17/2019] [Indexed: 01/07/2023] Open
Abstract
BACKGROUND Our understanding of the pig transcriptome is limited. RNA transcript diversity among nine tissues was assessed using poly(A) selected single-molecule long-read isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) from a single White cross-bred pig. RESULTS Across tissues, a total of 67,746 unique transcripts were observed, including 60.5% predicted protein-coding, 36.2% long non-coding RNA and 3.3% nonsense-mediated decay transcripts. On average, 90% of the splice junctions were supported by RNA-seq within tissue. A large proportion (80%) represented novel transcripts, mostly produced by known protein-coding genes (70%), while 17% corresponded to novel genes. On average, four transcripts per known gene (tpg) were identified; an increase over current EBI (1.9 tpg) and NCBI (2.9 tpg) annotations and closer to the number reported in human genome (4.2 tpg). Our new pig genome annotation extended more than 6000 known gene borders (5' end extension, 3' end extension, or both) compared to EBI or NCBI annotations. We validated a large proportion of these extensions by independent pig poly(A) selected 3'-RNA-seq data, or human FANTOM5 Cap Analysis of Gene Expression data. Further, we detected 10,465 novel genes (81% non-coding) not reported in current pig genome annotations. More than 80% of these novel genes had transcripts detected in > 1 tissue. In addition, more than 80% of novel intergenic genes with at least one transcript detected in liver tissue had H3K4me3 or H3K36me3 peaks mapping to their promoter and gene body, respectively, in independent liver chromatin immunoprecipitation data. CONCLUSIONS These validated results show significant improvement over current pig genome annotations.
Collapse
Affiliation(s)
- H Beiki
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA
| | - H Liu
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA
| | - J Huang
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA.,College of Animal Science and Technology, Jiangxi Agricultural University, Nanchang, Jiangxi, People's Republic of China
| | - N Manchanda
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, 819 Wallace Road, Ames, IA, 50011, USA
| | - D Nonneman
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, 68933, USA
| | - T P L Smith
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, 68933, USA
| | - J M Reecy
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA
| | - C K Tuggle
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA.
| |
Collapse
|