A platform independent RNA-Seq protocol for the detection of transcriptome complexity

BackgroundRecent studies have demonstrated an unexpected complexity of transcription in eukaryotes. The majority of the genome is transcribed and only a little fraction of these transcripts is annotated as protein coding genes and their splice variants. Indeed, most transcripts are the result of antisense, overlapping and non-coding RNA expression. In this frame, one of the key aims of high throughput transcriptome sequencing is the detection of all RNA species present in the cell and the first crucial step for RNA-seq users is represented by the choice of the strategy for cDNA library construction. The protocols developed so far provide the utilization of the entire library for a single sequencing run with a specific platform.ResultsWe set up a unique protocol to generate and amplify a strand-specific cDNA library representative of all RNA species that may be implemented with all major platforms currently available on the market (Roche 454, Illumina, ABI/SOLiD). Our method is reproducible, fast, easy-to-perform and even allows to start from low input total RNA. Furthermore, we provide a suitable bioinformatics tool for the analysis of the sequences produced following this protocol.ConclusionWe tested the efficiency of our strategy, showing that our method is platform-independent, thus allowing the simultaneous analysis of the same sample with different NGS technologies, and providing an accurate quantitative and qualitative portrait of complex whole transcriptomes.

[1]  B. Milash,et al.  Deep Sequencing for the Detection of Virus-Like Sequences in the Brains of Patients with Multiple Sclerosis: Detection of GBV-C in Human Brain , 2012, PloS one.

[2]  S Rozen,et al.  Primer3 on the WWW for general users and for biologist programmers. , 2000, Methods in molecular biology.

[3]  Piero Carninci,et al.  Whole genome transcriptome analysis , 2009, RNA biology.

[4]  G. Church,et al.  Sequencing genomes from single cells by polymerase cloning , 2006, Nature Biotechnology.

[5]  Giovanni Romeo,et al.  A mutation threshold distinguishes the antitumorigenic effects of the mitochondrial gene MTND1, an oncojanus function. , 2011, Cancer research.

[6]  Daniel H. Huson,et al.  48. MetaSim: A Sequencing Simulator for Genomics and Metagenomics , 2011 .

[7]  Thomas D. Wu,et al.  GMAP: a genomic mapping and alignment program for mRNA and EST sequence , 2005, Bioinform..

[8]  E. Wagner,et al.  Metabolism and regulation of canonical histone mRNAs: life without a poly(A) tail , 2008, Nature Reviews Genetics.

[9]  Z. Ning,et al.  Amplification-free Illumina sequencing-library preparation facilitates improved mapping and assembly of GC-biased genomes , 2009, Nature Methods.

[10]  Mehrdad Hajibabaei,et al.  Next‐generation sequencing technologies for environmental DNA research , 2012, Molecular ecology.

[11]  Lin Liu,et al.  Comparison of Next-Generation Sequencing Systems , 2012, Journal of biomedicine & biotechnology.

[12]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[13]  Rameen Beroukhim,et al.  Genome coverage and sequence fidelity of phi29 polymerase-based multiple strand displacement whole genome amplification. , 2004, Nucleic acids research.

[14]  F. Dean,et al.  Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle amplification. , 2001, Genome research.

[15]  Rickard Sandberg,et al.  How to analyze gene expression using RNA-sequencing data. , 2012, Methods in molecular biology.

[16]  Graziano Pesole,et al.  NGS-Trex: Next Generation Sequencing Transcriptome profile explorer , 2012, BMC Bioinformatics.

[17]  Erez Lieberman Aiden,et al.  The expanding scope of DNA sequencing , 2012, Nature Biotechnology.

[18]  Timothy B. Stockwell,et al.  Mechanism of chimera formation during the Multiple Displacement Amplification reaction , 2007, BMC biotechnology.

[19]  Informazioni Generali.  Genetica Medica , 1955 .

[20]  G. Pesole,et al.  Identification of tumor-associated cassette exons in human cancer through EST-based computational prediction and experimental validation , 2010, Molecular Cancer.

[21]  E. Bonora,et al.  Defective oxidative phosphorylation in thyroid oncocytic carcinoma is associated with pathogenic mitochondrial DNA mutations affecting complexes I and III. , 2006, Cancer research.

[22]  Daniel H. Huson,et al.  MetaSim—A Sequencing Simulator for Genomics and Metagenomics , 2008, PloS one.

[23]  R. Jensen,et al.  Differentially expressed alternatively spliced genes in Malignant Pleural Mesothelioma identified using massively parallel transcriptome sequencing , 2009, BMC Medical Genetics.

[24]  Data production leads,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[25]  Raymond K. Auerbach,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[26]  M. Deutscher,et al.  An important role for RNase R in mRNA decay. , 2005, Molecular cell.

[27]  James B. Brown,et al.  Global patterns of tissue-specific alternative polyadenylation in Drosophila. , 2012, Cell reports.

[28]  Michael Q. Zhang,et al.  Poly A- Transcripts Expressed in HeLa Cells , 2008, PloS one.

[29]  R. Daniel,et al.  Achievements and new knowledge unraveled by metagenomic approaches , 2009, Applied Microbiology and Biotechnology.

[30]  Ernesto Picardi,et al.  Bioinformatics approaches for genomics and post genomics applications of next-generation sequencing , 2010, Briefings Bioinform..

[31]  Amrutlal K. Patel,et al.  Identification of novel transcripts deregulated in buccal cancer by RNA-seq. , 2012, Gene.

[32]  B. Williams,et al.  Mapping and quantifying mammalian transcriptomes by RNA-Seq , 2008, Nature Methods.

[33]  O. Emanuelsson,et al.  Analysis of transcript and protein overlap in a human osteosarcoma cell line , 2010, BMC Genomics.

[34]  T. Steitz,et al.  A specific subdomain in phi29 DNA polymerase confers both processivity and strand-displacement capacity. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[35]  M. Attimonelli,et al.  The genetic and metabolic signature of oncocytic transformation implicates HIF1alpha destabilization. , 2010, Human molecular genetics.

[36]  Michael Zuker,et al.  Mfold web server for nucleic acid folding and hybridization prediction , 2003, Nucleic Acids Res..

[37]  Xi-rong Guo,et al.  microRNA expression profiling in fetal single ventricle malformation identified by deep sequencing. , 2011, International journal of molecular medicine.

[38]  M. Lovett,et al.  Combined deep microRNA and mRNA sequencing identifies protective transcriptomal signature of enhanced PI3Kα signaling in cardiac hypertrophy. , 2012, Journal of molecular and cellular cardiology.