RNA-Seq: a revolutionary tool for transcriptomics
TL;DR: The RNA-Seq approach to transcriptome profiling that uses deep-sequencing technologies provides a far more precise measurement of levels of transcripts and their isoforms than other methods.
Abstract: RNA-Seq is a recently developed approach to transcriptome profiling that uses deep-sequencing technologies. Studies using this method have already altered our view of the extent and complexity of eukaryotic transcriptomes. RNA-Seq also provides a far more precise measurement of levels of transcripts and their isoforms than other methods. This article describes the RNA-Seq approach, the challenges associated with its application, and the advances made so far in characterizing several eukaryote transcriptomes.
Citations
More filters
TL;DR: In this paper, a transcriptome analysis revealed that the slow symptomless growth was accompanied by minor metabolic responses and slightly suppressed defences in the host, whereas necrotrophic growth was associated with enhanced host responses involving energy metabolism, transport, signalling, defence and oxidative stress as well as a decrease in photosynthesis.
Abstract: The disease septoria leaf blotch of wheat, caused by fungal pathogen Septoria tritici, is of worldwide concern. The fungus exhibits a hemibiotrophic lifestyle, with a long symptomless, biotrophic phase followed by a sudden transition to necrotrophy associated with host necrosis. Little is known about the systematic interaction between fungal pathogenicity and host responses at specific growth stages and the factors triggering the transition. In order to gain some insights into global transcriptome alterations in both host and pathogen during the two phases of the compatible interaction, disease transition was monitored using pathogenesis-related gene markers and H2O2 signature prior to RNA-Seq. Transcriptome analysis revealed that the slow symptomless growth was accompanied by minor metabolic responses and slightly suppressed defences in the host, whereas necrotrophic growth was associated with enhanced host responses involving energy metabolism, transport, signalling, defence and oxidative stress as well as a decrease in photosynthesis. The fungus expresses distinct classes of stage-specific genes encoding potential effectors, probably first suppressing plant defence responses/facilitating the symptomless growth and later triggering life style transition and inducing host necrosis/facilitating the necrotrophic growth. Transport, signalling, anti-oxidative stress mechanisms and apoplastic nutrient acquisition play important roles in the entire infection process of S. tritici. Our findings uncover systematic S. tritici-induced expression profiles of wheat related to specific fungal infection strategies and provide a transcriptome resource for studying both hosts and pathogens in plant-Dothideomycete interactions.
93 citations
Cites background from "RNA-Seq: a revolutionary tool for t..."
...It can distinguish between paralogous genes, detect and quantify transcripts of low or high abundance and identify transcript sequence polymorphisms and novel trans-splicing and splice isoforms [17]....
[...]
TL;DR: This is the first study that characterized the integrated transcriptomic profiles during early development of penaeid shrimp, and these findings will serve as significant references for shrimp developmental biology and aquaculture research.
Abstract: Penaeid shrimp has a distinctive metamorphosis stage during early development. Although morphological and biochemical studies about this ontogeny have been developed for decades, researches on gene expression level are still scarce. In this study, we have investigated the transcriptomes of five continuous developmental stages in Pacific white shrimp (Litopenaeus vannamei) with high throughput Illumina sequencing technology. The reads were assembled and clustered into 66,815 unigenes, of which 32,398 have putative homologues in nr database, 14,981 have been classified into diverse functional categories by Gene Ontology (GO) annotation and 26,257 have been associated with 255 pathways by KEGG pathway mapping. Meanwhile, the differentially expressed genes (DEGs) between adjacent developmental stages were identified and gene expression patterns were clustered. By GO term enrichment analysis, KEGG pathway enrichment analysis and functional gene profiling, the physiological changes during shrimp metamorphosis could be better understood, especially histogenesis, diet transition, muscle development and exoskeleton reconstruction. In conclusion, this is the first study that characterized the integrated transcriptomic profiles during early development of penaeid shrimp, and these findings will serve as significant references for shrimp developmental biology and aquaculture research.
93 citations
Cites background from "RNA-Seq: a revolutionary tool for t..."
...The advent of RNA-Seq provides a far more high-throughput and precise measurement of levels of transcripts and their isoforms than other methods [14]....
[...]
...It has clear advantages over other approaches and is expected to revolutionize the manner in which eukaryotic transcriptomes are analyzed [14]....
[...]
TL;DR: The purpose of FusionHunter is to identify potential fusions with high sensitivity and specificity and to guide further functional validation in the laboratory, and it is shown that FusionHunter can accurately detect fusions that were previously confirmed by RT-PCR in a publicly available dataset.
Abstract: Motivation: Fusion transcripts can be created as a result of genome rearrangement in cancer. Some of them play important roles in carcinogenesis, and can serve as diagnostic and therapeutic targets. With more and more cancer genomes being sequenced by nextgeneration sequencing technologies, we believe an efficient tool for reliably identifying fusion transcripts will be desirable for many groups. Results: We designed and implemented an open-source software tool, called FusionHunter, which reliably identifies fusion transcripts from transcriptional analysis of paired-end RNA-seq. We show that FusionHunter can accurately detect fusions that were previously confirmed by RT-PCR in a publicly available dataset. The purpose of FusionHunter is to identify potential fusions with high sensitivity and specificity and to guide further functional validation in the laboratory. Availability: http://bioen-compbio.bioen.illinois.edu/FusionHunter/.
93 citations
TL;DR: The first target mRNAs were identified in archaea, including one sRNA that binds to the 5′-region of two mRN as in Methanosarcina mazei Gö1 and a few sRNAs that bind to 3′-UTRs in Sulfolobus solfataricus, three Pyrobaculum species, and Haloferax volcanii, indicating that archaeal sRN
Abstract: Small regulatory RNAs (sRNAs) are universally distributed in all three domains of life, Archaea, Bacteria, and Eukaryotes. In bacteria, sRNAs typically function by binding near the translation start site of their target mRNAs and thereby inhibit or activate translation. In eukaryotes, miRNAs and siRNAs typically bind to the 3′-untranslated region (3′-UTR) of their target mRNAs and influence translation efficiency and/or mRNA stability. In archaea, sRNAs have been identified in all species investigated using bioinformatic approaches, RNomics, and RNA-Seq. Their size can vary significantly between less than 50 to more than 500 nucleotides. Differential expression of sRNA genes has been studied using northern blot analysis, microarrays, and RNA-Seq. In addition, biological functions have been unraveled by genetic approaches, i.e., by characterization of designed mutants. As in bacteria, it was revealed that archaeal sRNAs are involved in many biological processes, including metabolic regulation, adaptation to extreme conditions, stress responses, and even in regulation of morphology and cellular behavior. Recently, the first target mRNAs were identified in archaea, including one sRNA that binds to the 5′-region of two mRNAs in Methanosarcina mazei Go1 and a few sRNAs that bind to 3′-UTRs in Sulfolobus solfataricus, three Pyrobaculum species, and Haloferax volcanii, indicating that archaeal sRNAs appear to be able to target both the 5′-UTR or the 3′-UTRs of their respective target mRNAs. In addition, archaea contain tRNA-derived fragments (tRFs), and one tRF has been identified as a major ribosome-binding sRNA in H. volcanii, which downregulates translation in response to stress. Besides regulatory sRNAs, archaea contain further classes of sRNAs, e.g., CRISPR RNAs (crRNAs) and snoRNAs.
93 citations
Cites methods from "RNA-Seq: a revolutionary tool for t..."
...Whole genome transcriptome analysis via high-throughput sequencing of cDNA libraries, called RNA-Seq, became available a few years ago and enabled the qualitative analysis of the RNA inventory of species as well as the quantitative analysis of differential transcript levels under various conditions.(39,40) The first RNA-Seq study with an archaeal species was performed for Methanosarcina mazei Gö1 under different nitrogen availabilities, leading to the identification of 242 intergenic and antisense sRNA, including six cis-antisense sRNAs overlapping with transposase genes and 40 sRNA candidates containing very short ORFs potentially encoding peptides smaller than 30 amino acids....
[...]
TL;DR: Several microRNAs and mRNAs involved in the drought response of tomato were identified using high-throughput sequencing, which will provide new insights into the complex regulatory network of plant adaption to drought stress.
Abstract: Abiotic stresses cause severe loss of crop production. Among them, drought is one of the most frequent environmental stresses, which limits crop growth, development and productivity. Plant drought tolerance is fine-tuned by a complex gene regulatory network. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success to improve plant yield and quality. Recent studies have demonstrated that microRNAs play critical roles in plant drought tolerance. However, little is known about the microRNA in drought response of the model plant tomato. Here, we described the profiling of drought-responsive microRNA and mRNA in tomato using high-throughput next-generation sequencing. Drought stress was applied on the seedlings of M82, a drought-sensitive cultivated tomato genotype, and IL9–1, a drought-tolerant introgression line derived from the stress-resistant wild species Solanum pennellii LA0716 and M82. Under drought, IL9–1 performed superior than M82 regarding survival rate, H2O2 elimination and leaf turgor maintenance. A total of four small RNA and eight mRNA libraries were constructed and sequenced using Illumina sequencing technology. 105 conserved and 179 novel microRNAs were identified, among them, 54 and 98 were differentially expressed upon drought stress, respectively. The majority of the differentially-expressed conserved microRNAs was up-regulated in IL9–1 whereas down-regulated in M82. Under drought stress, 2714 and 1161 genes were found to be differentially expressed in M82 and IL9–1, respectively, and many of their homologues are involved in plant stress, such as genes encoding transcription factor and protein kinase. Various pathways involved in abiotic stress were revealed by Gene Ontology and pathway analysis. The mRNA sequencing results indicated that most of the target genes were regulated by their corresponding microRNAs, which suggested that microRNAs may play essential roles in the drought tolerance of tomato. In this study, numerous microRNAs and mRNAs involved in the drought response of tomato were identified using high-throughput sequencing, which will provide new insights into the complex regulatory network of plant adaption to drought stress. This work will also help to exploit new players functioning in plant drought-stress tolerance.
92 citations
Cites background from "RNA-Seq: a revolutionary tool for t..."
...it also enable us to detect novel transcripts and alternative splicing events [58, 59]....
[...]
References
More filters
TL;DR: Although >90% of uniquely mapped reads fell within known exons, the remaining data suggest new and revised gene models, including changed or additional promoters, exons and 3′ untranscribed regions, as well as new candidate microRNA precursors.
Abstract: We have mapped and quantified mouse transcriptomes by deeply sequencing them and recording how frequently each gene is represented in the sequence sample (RNA-Seq). This provides a digital measure of the presence and prevalence of transcripts from known and previously unknown genes. We report reference measurements composed of 41–52 million mapped 25-base-pair reads for poly(A)-selected RNA from adult mouse brain, liver and skeletal muscle tissues. We used RNA standards to quantify transcript prevalence and to test the linear range of transcript detection, which spanned five orders of magnitude. Although >90% of uniquely mapped reads fell within known exons, the remaining data suggest new and revised gene models, including changed or additional promoters, exons and 3′ untranscribed regions, as well as new candidate microRNA precursors. RNA splice events, which are not readily measured by standard gene expression microarray or serial analysis of gene expression methods, were detected directly by mapping splice-crossing sequence reads. We observed 1.45 × 10 5 distinct splices, and alternative splices were prominent, with 3,500 different genes expressing one or more alternate internal splices. The mRNA population specifies a cell’s identity and helps to govern its present and future activities. This has made transcriptome analysis a general phenotyping method, with expression microarrays of many kinds in routine use. Here we explore the possibility that transcriptome analysis, transcript discovery and transcript refinement can be done effectively in large and complex mammalian genomes by ultra-high-throughput sequencing. Expression microarrays are currently the most widely used methodology for transcriptome analysis, although some limitations persist. These include hybridization and cross-hybridization artifacts 1–3 , dye-based detection issues and design constraints that preclude or seriously limit the detection of RNA splice patterns and previously unmapped genes. These issues have made it difficult for standard array designs to provide full sequence comprehensiveness (coverage of all possible genes, including unknown ones, in large genomes) or transcriptome comprehensiveness (reliable detection of all RNAs of all prevalence classes, including the least abundant ones that are physiologically relevant). Other
12,293 citations
TL;DR: Serial analysis of gene expression (SAGE) should provide a broadly applicable means for the quantitative cataloging and comparison of expressed genes in a variety of normal, developmental, and disease states.
Abstract: PROBLEM TO BE SOLVED: To provide a method for preparing a short nucleotide sequence (tag) which is useful to identify a cDNA oligonucleotide and is derived from a restricted position in a mRNA or a cDNA. SOLUTION: This is the method of preparing a tag for identifying the cDNA oligonucleotide. The above method comprises preparing the cDNA oligonucleotide bearing 5' and 3' terminals, collecting cDNA fragments by cutting the cDNA oligonucleotide with a restriction enzyme at the first restriction endonuclease site, separating a cDNA oligonucleotide bearing 5' or 3' terminal and connecting an oligonucleotide linker to the isolated cDNA fragment bearing the cDNA oligonucleotide 5' or 3' terminal. Here, the oligonucleotide linker contains the recognition site of the second restriction endonuclease enzyme and the isolated cDNA fragment is cut with the second restriction endonuclease enzyme which cuts the cDNA fragment in a section separated from the recognition site to obtain the tag for identifying the cDNA oligonucleotide.
4,437 citations
TL;DR: This work describes the software MAQ, software that can build assemblies by mapping shotgun short reads to a reference genome, using quality scores to derive genotype calls of the consensus sequence of a diploid genome, e.g., from a human sample.
Abstract: New sequencing technologies promise a new era in the use of DNA sequence. However, some of these technologies produce very short reads, typically of a few tens of base pairs, and to use these reads effectively requires new algorithms and software. In particular, there is a major issue in efficiently aligning short reads to a reference genome and handling ambiguity or lack of accuracy in this alignment. Here we introduce the concept of mapping quality, a measure of the confidence that a read actually comes from the position it is aligned to by the mapping algorithm. We describe the software MAQ that can build assemblies by mapping shotgun short reads to a reference genome, using quality scores to derive genotype calls of the consensus sequence of a diploid genome, e.g., from a human sample. MAQ makes full use of mate-pair information and estimates the error probability of each read alignment. Error probabilities are also derived for the final genotype calls, using a Bayesian statistical model that incorporates the mapping qualities, error probabilities from the raw sequence quality scores, sampling of the two haplotypes, and an empirical model for correlated errors at a site. Both read mapping and genotype calling are evaluated on simulated data and real data. MAQ is accurate, efficient, versatile, and user-friendly. It is freely available at http://maq.sourceforge.net.
2,927 citations
TL;DR: It is found that the Illumina sequencing data are highly replicable, with relatively little technical variation, and thus, for many purposes, it may suffice to sequence each mRNA sample only once (i.e., using one lane).
Abstract: Ultra-high-throughput sequencing is emerging as an attractive alternative to microarrays for genotyping, analysis of methylation patterns, and identification of transcription factor binding sites. Here, we describe an application of the Illumina sequencing (formerly Solexa sequencing) platform to study mRNA expression levels. Our goals were to estimate technical variance associated with Illumina sequencing in this context and to compare its ability to identify differentially expressed genes with existing array technologies. To do so, we estimated gene expression differences between liver and kidney RNA samples using multiple sequencing replicates, and compared the sequencing data to results obtained from Affymetrix arrays using the same RNA samples. We find that the Illumina sequencing data are highly replicable, with relatively little technical variation, and thus, for many purposes, it may suffice to sequence each mRNA sample only once (i.e., using one lane). The information in a single lane of Illumina sequencing data appears comparable to that in a single array in enabling identification of differentially expressed genes, while allowing for additional analyses such as detection of low-expressed genes, alternative splice variants, and novel transcripts. Based on our observations, we propose an empirical protocol and a statistical framework for the analysis of gene expression using ultra-high-throughput sequencing technology.
2,834 citations
TL;DR: The program SOAP is designed to handle the huge amounts of short reads generated by parallel sequencing using the new generation Illumina-Solexa sequencing technology, which supports multi-threaded parallel computing and has a batch module for multiple query sets.
Abstract: Summary: We have developed a program SOAP for efficient gapped and ungapped alignment of short oligonucleotides onto reference sequences. The program is designed to handle the huge amounts of short reads generated by parallel sequencing using the new generation Illumina-Solexa sequencing technology. SOAP is compatible with numerous applications, including single-read or pair-end resequencing, small RNA discovery and mRNA tag sequence mapping. SOAP is a command-driven program, which supports multi-threaded parallel computing, and has a batch module for multiple query sets. Availability: http://soap.genomics.org.cn Contact: soap@genomics.org.cn
2,729 citations
"RNA-Seq: a revolutionary tool for t..." refers methods in this paper
...There are several programs for mapping reads to the genome, including ELAND, SOA...
[...]