SQANTI: extensive characterization of long-read transcript sequences for quality control in full-length transcriptome identification and quantification

Abstract

High-throughput sequencing of full-length transcripts using long reads has paved the way for the discovery of thousands of novel transcripts, even in well-annotated mammalian species. The advances in sequencing technology have created a need for studies and tools that can characterize these novel variants. Here, we present SQANTI, an automated pipeline for the classification of long-read transcripts that can assess the quality of data and the preprocessing pipeline using 47 unique descriptors. We apply SQANTI to a neuronal mouse transcriptome using Pacific Biosciences (PacBio) long reads and illustrate how the tool is effective in characterizing and describing the composition of the full-length transcriptome. We perform extensive evaluation of ToFU PacBio transcripts by PCR to reveal that an important number of the novel transcripts are technical artifacts of the sequencing approach and that SQANTI quality descriptors can be used to engineer a filtering strategy to remove them. Most novel transcripts in this curated transcriptome are novel combinations of existing splice sites, resulting more frequently in novel ORFs than novel UTRs, and are enriched in both general metabolic and neural-specific functions. We show that these new transcripts have a major impact in the correct quantification of transcript levels by state-of-the-art short-read-based quantification algorithms. By comparing our iso-transcriptome with public proteomics databases, we find that alternative isoforms are elusive to proteogenomics detection. SQANTI allows the user to maximize the analytical outcome of long-read technologies by providing the tools to deliver quality-evaluated and curated full-length transcriptomes.

References 82

1. 1.
   A survey of the sorghum transcriptome using single-molecule long reads · 2016

   No notice
2. 2.

   Amaral AI, Hadera MG, Tavares JM, Kotter MRN, Sonnewald U. 2016 Characterization of glucose-related metabolic pathways in differentiated rat oligodendrocyte lineage cells. Glia 64: 21–34.

   Not in our library10.1002/glia.22900

3. 3.
   Improving PacBio Long Read Accuracy by Short Read Alignment · 2012

   No notice
4. 4.

   Au KF, Sebastiano V, Afshar PT, Durruthy JD, Lee L, Williams BA, van Bakel H, Schadt EE, Reijo-Pera RA, Underwood JG, 2013 Characterization of the human ESC transcriptome by hybrid sequencing. Proc Natl Acad Sci 110: E4821–E4830.

   Not in our library10.1073/pnas.1320101110

5. 5.

   Black DL. 2003 Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem 72: 291–336.

   Not in our library10.1146/annurev.biochem.72.121801.161720

6. 6.

   Breiman L. 2001 Random forests. Mach Learn 45: 5–32.

   Not in our library

7. 7.

   Breiman L, Friedman JH, Olshen RA, Stone CJ. 1984 Classification and regression trees. Wadsworth & Brooks/Cole Advanced Books & Software, Monterey, CA.

   Not in our library

8. 8.

   Cai J, Cheng A, Luo Y, Lu C, Mattson MP, Rao MS, Furukawa K. 2004 Membrane properties of rat embryonic multipotent neural stem cells. J Neurochem 88: 212–226.

   Not in our library10.1046/j.1471-4159.2003.02184.x

9. 9.
   Pacific biosciences sequencing technology for genotyping and variation discovery in human data · 2012

   No notice
10. 10.

    Cartolano M, Huettel B, Hartwig B, Reinhardt R, Schneeberger K. 2016 cDNA library enrichment of full length transcripts for SMRT long read sequencing. PLoS One 11: e0157779.

    Not in our library10.1371/journal.pone.0157779

11. 11.

    Chen J, Weiss WA. 2014 Alternative splicing in cancer: implications for biology and therapy. Oncogene 34: 1–14.

    Not in our library10.1038/onc.2013.570

12. 12.

    Cliff TS, Dalton S. 2017 Metabolic switching and cell fate decisions: implications for pluripotency, reprogramming and development. Curr Opin Genet Dev 46: 44–49.

    Not in our library10.1016/j.gde.2017.06.008

13. 13.

    Cocquet J, Chong A, Zhang G, Veitia RA. 2006 Reverse transcriptase template switching and false alternative transcripts. Genomics 88: 127–131.

    Not in our library10.1016/j.ygeno.2005.12.013

14. 14.

    Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. 2005 Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676.

    Not in our library10.1093/bioinformatics/bti610

15. 15.

    Cortes C, Vapnik V. 1995 Support-vector networks. Mach Learn 20: 273–297.

    Not in our library

16. 16.

    Craig R, Beavis RC. 2004 TANDEM: matching proteins with tandem mass spectra. Bioinformatics 20: 1466–1467.

    Not in our library10.1093/bioinformatics/bth092

17. 17.

    Deutsch EW, Overall CM, Van Eyk JE, Baker MS, Paik YK, Weintraub ST, Lane L, Martens L, Vandenbrouck Y, Kusebauch U, 2016 Human Proteome Project Mass Spectrometry Data Interpretation Guidelines 2.1. J Proteome Res 15: 3961–3970.

    Not in our library10.1021/acs.jproteome.6b00392

18. 18.

    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. 2013 STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29: 15–21.

    Not in our library10.1093/bioinformatics/bts635

19. 19.

    D'Souza I, Poorkaj P, Hong M, Nochlin D, Lee VM, Bird TD, Schellenberg GD. 1999 Missense and silent τ gene mutations cause frontotemporal dementia with parkinsonism-chromosome 17 type, by affecting multiple alternative RNA splicing regulatory elements. Proc Natl Acad Sci 96: 5598–5603.

    Not in our library10.1073/pnas.96.10.5598

20. 20.
    The Human Pancreatic Islet Transcriptome: Expression of Candidate Genes for Type 1 Diabetes and the Impact of Pro-Inflammatory Cytokines · 2012

    No notice
21. 21.

    Eng JK, McCormack AL, Yates JR. 1994 An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J Am Soc Mass Spectrom 5: 976–989.

    Not in our library10.1016/1044-0305(94)80016-2

22. 22.

    Ezkurdia I, Rodriguez JM, Carrillo-De Santa Pau E, Vázquez J, Valencia A, Tress ML. 2015 Most highly expressed protein-coding genes have a single dominant isoform. J Proteome Res 14: 1880–1887.

    Not in our library10.1021/pr501286b

23. 23.

    Frankish A, Mudge JM, Thomas M, Harrow J. 2012 The importance of identifying alternative splicing in vertebrate genome annotation. Database 2012: bas014.

    Not in our library10.1093/database/bas014

24. 24.
    Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing · 2015

    No notice
25. 25.

    Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A. 2008 High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36: 3420–3435.

    Not in our library10.1093/nar/gkn176

26. 26.

    Hackl T, Hedrich R, Schultz J, Förster F. 2014 proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics 30: 1–8.

    Not in our library10.1093/bioinformatics/btu392

27. 27.

    Han D, Jin J, Woo J, Min H, Kim Y. 2014 Proteomic analysis of mouse astrocytes and their secretome by a combination of FASP and StageTip-based, high pH, reversed-phase fractionation. Proteomics 14: 1604–1609.

    Not in our library10.1002/pmic.201300495

28. 28.
    Apparent Non-Canonical Trans-Splicing Is Generated by Reverse Transcriptase In Vitro · 2010

    No notice
29. 29.

    Irimia M, Weatheritt RJ, Ellis JD, Parikshak NN, Gonatopoulos-Pournatzis T, Babor M, Quesnel-Vallières M, Tapial J, Raj B, O'Hanlon D, 2014 A highly conserved program of neuronal microexons is misregulated in autistic brains. Cell 159: 1511–1523.

    Not in our library10.1016/j.cell.2014.11.035

30. 30.

    Kanadia RN, Johnstone KA, Mankodi A, Lungu C, Thornton CA, Esson D, Timmers AM, Hauswirth WW, Swanson MS. 2003 A muscleblind knockout model for myotonic dystrophy. Science 302: 1978–1980.

    Not in our library10.1126/science.1088583

31. 31.

    Kim S, Gupta N, Pevzner PA. 2008 Spectral probabilities and generating functions of tandem mass spectra: a strike against decoy databases. J Proteome Res 7: 3354–3363.

    Not in our library10.1021/pr8001244

32. 32.

    La Cognata V, Iemmolo R, D'Agata V, Scuderi S, Drago F, Zappia M, Cavallaro S. 2014 Increasing the coding potential of genomes through alternative splicing: the case of PARK2 gene. Curr Genomics 15: 203–216.

    Not in our library10.2174/1389202915666140426003342

33. 33.

    Ladd AN. 2013 CUG-BP, Elav-like family (CELF)-mediated alternative splicing regulation in the brain during health and disease. Mol Cell Neurosci 56: 456–464.

    Not in our library10.1016/j.mcn.2012.12.003

34. 34.
    Aberrant Alternative Splicing Is Another Hallmark of Cancer · 2013

    No notice
35. 35.

    Langmead B, Salzberg SL. 2012 Fast gapped-read alignment with Bowtie 2. Nat Methods 9: 357–359.

    Not in our library10.1038/nmeth.1923

36. 36.

    Lee JA, Damianov A, Lin CH, Fontes M, Parikshak NN, Anderson ES, Geschwind DH, Black DL, Martin KC. 2016 Cytoplasmic Rbfox1 regulates the expression of synaptic and autism-related genes. Neuron 89: 113–128.

    Not in our library10.1016/j.neuron.2015.11.025

37. 37.
    RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome · 2011

    No notice
38. 38.

    Li S, Yamada M, Han X, Ohler U, Benfey PN. 2016 High-resolution expression map of the Arabidopsis root reveals alternative splicing and lincRNA regulation. Dev Cell 39: 508–522.

    Not in our library10.1016/j.devcel.2016.10.012

39. 39.

    Liu S, Cheng C. 2013 Alternative RNA splicing and cancer. Wiley Interdiscip Rev RNA 4: 547–566.

    Not in our library10.1002/wrna.1178

40. 40.

    Loh W-Y, Shih Y-S. 1997 Split selection methods for classification trees. Stat Sin 7: 815–840.

    Not in our library

Cited by papers we hold 19

• SQANTI3: curation of long-read transcriptomes for accurate identification of known and novel isoforms · 2024
  No notice
• Systematic assessment of long-read RNA-seq methods for transcript identification and quantification · 2024
  No notice
• Functional annotation of the animal genomes: An integrated annotation resource for the horse. · 2023
  Not checked
• Single-cell isoform analysis in human immune cells · 2022
  Not checked
• Reference genome assemblies reveal the origin and evolution of allohexaploid oat · 2022
  No notice
• Error correction enables use of Oxford Nanopore technology for reference-free transcriptome analysis · 2021
  Corrected
• Detection of aberrant splicing events in RNA-seq data using FRASER · 2021
  Corrected
• isoCirc catalogs full-length circular RNA isoforms in human transcriptomes · 2021
  No notice
• Comprehensive characterization of single-cell full-length isoforms in human and mouse with long-read sequencing · 2021
  No notice
• Single-cell RNA-seq analysis of mouse preimplantation embryos by third-generation sequencing. · 2020
  Not checked
• Comprehensive identification of mRNA isoforms reveals the diversity of neural cell-surface molecules with roles in retinal development and disease · 2020
  Not checked
• Opportunities and challenges in long-read sequencing data analysis · 2020
  No notice
• High throughput error corrected Nanopore single cell transcriptome sequencing · 2020
  No notice
• Getting the Entire Message: Progress in Isoform Sequencing · 2019
  Not checked
• Realizing the potential of full-length transcriptome sequencing · 2019
  No notice
• Transcriptome assembly from long-read RNA-seq alignments with StringTie2 · 2019
  No notice
• A comprehensive examination of Nanopore native RNA sequencing for characterization of complex transcriptomes · 2019
  No notice
• The PRIDE database and related tools and resources in 2019: improving support for quantification data · 2018
  No notice
• Single-cell RNAseq for the study of isoforms—how is that possible? · 2018
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.