Current progress in the field of next-generation transcriptome sequencing have contributed significantly to the study of various malignancies including glioblastoma multiforme (GBM). Differential sequencing of transcriptomes of patients and non-tumor controls has a potential to reveal novel transcripts with significant role in GBM. One such candidate group of molecules are long non-coding RNAs (lncRNAs) which have been proved to be involved in processes such as carcinogenesis, epigenetic modifications and resistance to various therapeutic approaches. To maximize the value of transcriptome sequencing, a proper protocol for library preparation from tissue-derived RNA needs to be found which would produce high quality transcriptome sequencing data and increase the number of detected lncRNAs. It is important to mention that success of library preparation is determined by the quality of input RNA, which is in case of real-life tissue specimens very often altered in comparison to high quality RNA commonly used by manufacturers for development of library preparation chemistry. In the present study, we used GBM and non-tumor brain tissue specimens and compared three different commercial library preparation kits, namely NEXTflex Rapid Directional qRNA-Seq Kit (Bioo Scientific), SENSE Total RNA-Seq Library Prep Kit (Lexogen) and NEBNext Ultra II Directional RNA Library Prep Kit for Illumina (NEB). Libraries generated using SENSE kit were characterized by the most normal distribution of normalized average GC content, the least amount of over-represented sequences and the percentage of ribosomal RNA reads (0.3-1.5%) and highest numbers of uniquely mapped reads and reads aligning to coding regions. However, NEBNext kit performed better having relatively low duplication rates, even transcript coverage and the highest number of hits in Ensembl database for every biotype of our interest including lncRNAs. Our results indicate that out of three approaches the NEBNext library preparation kit was most suitable for the study of lncRNAs via transcriptome sequencing. This was further confirmed by highly consistent data reached in an independent validation on an expanded cohort.

1st Department of Pathological Anatomy St Anne's University Hospital and Faculty of Medicine Masaryk University Brno Czech Republic

Central European Institute of Technology Masaryk University Brno Czech Republic

Department of Comprehensive Cancer Care Masaryk Memorial Cancer Institute Faculty of Medicine Masaryk University Brno Czech Republic

Department of Neurosurgery St Anne's University Hospital and Faculty of Medicine Masaryk University Brno Czech Republic

Department of Neurosurgery University Hospital Brno Faculty of Medicine Masaryk University Brno Czech Republic

Department of Neurosurgery University Hospital Ostrava Ostrava Czech Republic

Department of Pathology University Hospital Brno Faculty of Medicine Masaryk University Brno Czech Republic

See more in PubMed

Prabhakar B, Zhong XB, Rasmussen TP. Exploiting Long Noncoding RNAs as Pharmacological Targets to Modulate Epigenetic Diseases. Yale J Biol Med. 2017;90(1):73–86. Epub 2017/03/31. PubMed PMC

Yu H, Xue Y, Wang P, Liu X, Ma J, Zheng J, et al. Knockdown of long non-coding RNA XIST increases blood-tumor barrier permeability and inhibits glioma angiogenesis by targeting miR-137. Oncogenesis. 2017;6(3):e303 Epub 2017/03/14. 10.1038/oncsis.2017.7 PubMed DOI PMC

Lu YF, Cai XL, Li ZZ, Lv J, Xiang YA, Chen JJ, et al. LncRNA SNHG16 Functions as an Oncogene by Sponging MiR-4518 and Up-Regulating PRMT5 Expression in Glioma. Cell Physiol Biochem. 2018;45(5):1975–85. Epub 2018/03/13. 10.1159/000487974 . PubMed DOI

Tyler AD, Christianson S, Knox NC, Mabon P, Wolfe J, Van Domselaar G, et al. Comparison of Sample Preparation Methods Used for the Next-Generation Sequencing of Mycobacterium tuberculosis. PLoS One. 2016;11(2):e0148676 Epub 2016/02/06. 10.1371/journal.pone.0148676 PubMed DOI PMC

Head SR, Komori HK, LaMere SA, Whisenant T, Van Nieuwerburgh F, Salomon DR, et al. Library construction for next-generation sequencing: overviews and challenges. Biotechniques. 2014;56(2):61–4, 6, 8, passim. Epub 2014/02/08. 10.2144/000114133 PubMed DOI PMC

Guil S, Esteller M. Cis-acting noncoding RNAs: friends and foes. Nat Struct Mol Biol. 2012;19(11):1068–75. Epub 2012/11/08. 10.1038/nsmb.2428 . PubMed DOI

Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464(7291):1071–6. Epub 2010/04/16. 10.1038/nature08975 PubMed DOI PMC

Gao YF, Wang ZB, Zhu T, Mao CX, Mao XY, Li L, et al. A critical overview of long non-coding RNA in glioma etiology 2016: an update. Tumour Biol. 2016;37(11):14403–13. Epub 2016/09/16. 10.1007/s13277-016-5307-4 . PubMed DOI

Spizzo R, Almeida MI, Colombatti A, Calin GA. Long non-coding RNAs and cancer: a new frontier of translational research? Oncogene. 2012;31(43):4577–87. Epub 2012/01/24. 10.1038/onc.2011.621 PubMed DOI PMC

Rinn JL, Chang HY. Genome regulation by long noncoding RNAs. Annu Rev Biochem. 2012;81:145–66. Epub 2012/06/06. 10.1146/annurev-biochem-051410-092902 PubMed DOI PMC

Sultan M, Dokel S, Amstislavskiy V, Wuttig D, Sultmann H, Lehrach H, et al. A simple strand-specific RNA-Seq library preparation protocol combining the Illumina TruSeq RNA and the dUTP methods. Biochem Biophys Res Commun. 2012;422(4):643–6. Epub 2012/05/23. 10.1016/j.bbrc.2012.05.043 . PubMed DOI

Wang L, Si Y, Dedow LK, Shao Y, Liu P, Brutnell TP. A low-cost library construction protocol and data analysis pipeline for Illumina-based strand-specific multiplex RNA-seq. PLoS One. 2011;6(10):e26426 Epub 2011/11/01. 10.1371/journal.pone.0026426 PubMed DOI PMC

Zhou K, Zhang C, Yao H, Zhang X, Zhou Y, Che Y, et al. Knockdown of long non-coding RNA NEAT1 inhibits glioma cell migration and invasion via modulation of SOX2 targeted by miR-132. Mol Cancer. 2018;17(1):105 Epub 2018/07/29. 10.1186/s12943-018-0849-2 PubMed DOI PMC

Curwen V, Eyras E, Andrews TD, Clarke L, Mongin E, Searle SM, et al. The Ensembl automatic gene annotation system. Genome Res. 2004;14(5):942–50. Epub 2004/05/05. 10.1101/gr.1858004 PubMed DOI PMC

Sims D, Sudbery I, Ilott NE, Heger A, Ponting CP. Sequencing depth and coverage: key considerations in genomic analyses. Nat Rev Genet. 2014;15(2):121–32. Epub 2014/01/18. 10.1038/nrg3642 . PubMed DOI

Parekh S, Ziegenhain C, Vieth B, Enard W, Hellmann I. The impact of amplification on differential expression analyses by RNA-seq. Sci Rep. 2016;6:25533 Epub 2016/05/10. 10.1038/srep25533 PubMed DOI PMC

Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. Epub 2012/10/30. 10.1093/bioinformatics/bts635 PubMed DOI PMC

Li H, Yuan X, Yan D, Li D, Guan F, Dong Y, et al. Long Non-Coding RNA MALAT1 Decreases the Sensitivity of Resistant Glioblastoma Cell Lines to Temozolomide. Cell Physiol Biochem. 2017;42(3):1192–201. Epub 2017/07/03. 10.1159/000478917 . PubMed DOI

Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46(D1):D754–D61. Epub 2017/11/21. 10.1093/nar/gkx1098 PubMed DOI PMC

Sayols S, Scherzinger D, Klein H. dupRadar: a Bioconductor package for the assessment of PCR artifacts in RNA-Seq data. BMC Bioinformatics. 2016;17(1):428 Epub 2016/10/23. 10.1186/s12859-016-1276-2 PubMed DOI PMC

Andrews S. FastQC: a quality control tool for high throughput sequence data. Available from: http://www.bioinformatics.babraham.ac.uk/projects/fastqc.

Krishnaswami SR, Grindberg RV, Novotny M, Venepally P, Lacar B, Bhutani K, et al. Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons. Nat Protoc. 2016;11(3):499–524. Epub 2016/02/20. 10.1038/nprot.2016.015 PubMed DOI PMC

Qin N, Tong GF, Sun LW, Xu XL. Long Noncoding RNA MEG3 Suppresses Glioma Cell Proliferation, Migration, and Invasion by Acting as a Competing Endogenous RNA of miR-19a. Oncol Res. 2017;25(9):1471–8. Epub 2017/03/10. 10.3727/096504017X14886689179993 . PubMed DOI PMC

Gong X, Huang M. Long non-coding RNA MEG3 promotes the proliferation of glioma cells through targeting Wnt/beta-catenin signal pathway. Cancer Gene Ther. 2017;24(9):381–5. Epub 2017/10/14. 10.1038/cgt.2017.32 . PubMed DOI

Martin M. cutadapt Documentation, Release 1.17. Available from: https://media.readthedocs.org/pdf/cutadapt/stable/cutadapt.pdf.

Ewels P, Magnusson M, Lundin S, Kaller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8. Epub 2016/06/18. 10.1093/bioinformatics/btw354 PubMed DOI PMC