Detecting alternatively spliced transcript isoforms from single‐molecule long‐read sequences without a reference genome. (5th June 2017)

• Record Type: Journal Article
• Authors: Liu, Xiaoxian; Mei, Wenbin; Soltis, Pamela S.; Soltis, Douglas E.; Barbazuk, W. Brad
• Is Part Of: Molecular ecology resources. Volume 17:Number 6(2017:Nov.)
• Page Start: 1243
• Page End: 1256
• Publication Date: 2017-06-05
• Subjects: alternative splicing -- Amborella trichopoda -- pipeline -- single‐molecule long‐read sequencing -- transcriptomics; Molecular ecology -- Periodicals; 572.8
• DOI: 10.1111/1755-0998.12670
• Languages: English
• ISSNs: 1755-098X

Abstract

Abstract: Alternative splicing (AS) is a major source of transcript and proteome diversity, but examining AS in species without well‐annotated reference genomes remains difficult. Research on both human and mouse has demonstrated the advantages of using Iso‐Seq™ data for isoform‐level transcriptome analysis, including the study of AS and gene fusion. We applied Iso‐Seq™ to investigate AS in Amborella trichopoda, a phylogenetically pivotal species that is sister to all other living angiosperms. Our data show that, compared with RNA‐Seq data, the Iso‐Seq™ platform provides better recovery on large transcripts, new gene locus identification and gene model correction. Reference‐based AS detection with Iso‐Seq™ data identifies AS within a higher fraction of multi‐exonic genes than observed for published RNA‐Seq analysis (45.8% vs. 37.5%). These data demonstrate that the Iso‐Seq™ approach is useful for detecting AS events. Using the Iso‐Seq‐defined transcript collection in Amborella as a reference, we further describe a pipeline for detection of AS isoforms from PacBio Iso‐Seq™ without using a reference sequence ( de novo ). Results using this pipeline show a 66%–76% overall success rate in identifying AS events. This de novo AS detection pipeline provides a method to accurately characterize and identify bona fide alternatively spliced transcripts in any nonmodel system that lacks a reference genome sequence. Hence, our pipeline has huge potential applications and benefits to the broader biology community.