View source for BioMicroCenter:RNAseq

{{BioMicroCenter}}
[[Image:RNA-Seq_pres.jpg|thumb|300px|right|Wang ET, et al. Nature 2008]]
The BioMicro Center supports several methods for RNAseq. The choice of method is highly dependent on the amount of input RNA and the quality of the input RNA. At this time, we only accept eukaryotic samples for RNAseq but we are willing to discuss exploring prokaryotic method development.

== RNAseq Methods ==

{| border=1
 ! Amount of RNA
 ! Quality of RNA
 ! Method Recommended
 |-
 | >100ng || RIN:9.0 || Illumina TruSeq
 |-
 | >1ug || RIN:5.0 || Epicenter RiboZero
 |- 
 | 10pg-100ng || RIN:9.0 || Clontech SMARTer Low Input
 |- 
 | 1ng-1ug || RIN:5.0 || NuGEN Ovation RNA-Seq System V2
|}

===[http://www.illumina.com/products/truseq_rna_sample_prep_kit_v2.ilmn Illumina TruSeq]===
{|
|- valign="top"
| 
Illumina's TruSeq chemistry is the primary RNAseq method used in the BioMicro Center. This chemistry uses polyT beads to isolate the mRNA from the rRNA and tRNA. The use of these beads requires that the RNA be of very high quality or only the 3' end of transcripts will be isolated. Purified mRNA is then fragmented with metal and random priming is used to convert the sample to cDNA. Once double-stranded cDNA is generated, the sample is transferred to the [[BioMicroCenter:SPRI-Works|SPRIworks]] for the remainder of sample preparation.
|
[[Image:BMC_RNAseqMethod.png|thumb|right|200px|TruSeq Chemistry]]
|
[[Image:BMC_IlluminaRNAseq.png|thumb|right|300px|Sample Data]]
|}
=== [http://www.epibio.com/item.asp?ID=589 EpiCenter RiboZero] ===
For samples with degraded RNA or samples where you are interested in looking at non-polyA RNAs, the BioMicro Center utilizes the Epicenter ribominus kit. This kit uses *check me* magnetic beads coupled to rRNA and tRNA sequences to remove these sequences from the solution. The remaining mRNA fragments can then be converted in to cDNA. Once double-stranded cDNA is generated, the sample is transferred to the [[BioMicroCenter:SPRI-Works|SPRIworks]] for the remainder of sample preparation.

=== [http://www.clontech.com/US/Products/cDNA_Synthesis_and_Library_Construction/cDNA_Synthesis_Kits/Ultra_Low_Input_RNA_cDNA_Synthesis Clontech SMARTer Low-Input] ===
{|
|- valign="top"
|
For samples with less then 100ng of input, the BioMicro Center utilizes the [http://www.clontech.com/US/Products/cDNA_Synthesis_and_Library_Construction/cDNA_Synthesis_Kits/Ultra_Low_Input_RNA_cDNA_Synthesis Clontech SMARTer system]. This system differs from the TruSeq chemisry in that it begins with cDNA generation using polyT priming and proprietary chemistry. The use of polyT priming requires the RNA to be of high quality. Full length double-stranded cDNAs are generated and amplified by PCR. These cDNAs are then fragmented and transferred to the [[BioMicroCenter:SPRI-Works|SPRIworks]] for the remainder of sample preparation. Data from this system is of similar quality to samples created with Illumina TruSeq chemistry.
|
[[Image:BMC_ClontechChemistry.png|thumb|right|200px|Clontech system.<BR>Image from Clonetech]]
|}
=== [http://www.nugeninc.com/nugen/index.cfm/products/cs/ngs/rna-seq-v2/ NuGEN Ovation RNAseq System V2] === 
{|
|- valign="top"
|
For samples with less then 100ng of input and restricted input amounts, our kit of choice is the NuGEN Ovation. This kit utilizes non-random nonamers designed to not amplify ribosomal RNA to create double stranded cDNA fragments. The fragmented cDNA is then transferred to the [[BioMicroCenter:SPRI-Works|SPRIworks]] for the remainder of sample preparation. The downside of this kit is that the non-random nonamers cannot amplify all areas of the genome and certain portions of genes are often lost. Still, for many samples, this kit is the only way to generate RNAseq libraries.
|
[[Image:BMC_SPIA.gif|thumb|right|200px|NuGEN SPIA Chemistry. <BR> Image from NuGen.]]
|}
== Additional Chemistries Available in the BioMicro Center ==
=== Strand Specific Sequencing ===
{|
 |-
 |For samples with high amounts of input, we can modify the chemistry of cDNA creation to allow detecting the strand of the RNA using the dUTP 2nd strand marking protocol that preformed best in J. Levin et al 2010. In this method, actinomycin D is added to the 1st strand synthesis to prevent reinitiation of the RT-PCR enzymes. Then, dUTP is substituted in for dTTP in the second strand, allowing a clear distinction between the forward and reverse strands. After library construction, but before amplification, the dUTP containing strand is degraded and only the reverse strand (template strand) remains. This method is highly efficient but requires significant amount of RNA as the actinomycin causes a significant reduction in 1st strand yield. 
 |-
 |[[Image:BMC_strandSpecific.gif|thumb|600px|left|Sample Data from Strand Specific RNAseq]]
|}
=== [[BioMicroCenter:PippinPrep|Size Selection]] ===
For some applications of RNAseq, such as splice choice determination, having a precise knowledge of the insert size is critical. While the [[BioMicroCenter:SPRI-Works|SPRIworks]] does provide some size selection (typically restricting fragments to between 150 and 350bp), this can be too wide for some methodologies. In these cases, after libraries are amplified, they can be run on the [[BioMicroCenter:PippinPrep|Sage BluePippin]] (either singly or pooled). Here the size distribution can be much tighter, with most of the DNA fragments being within a 50nt range.

== Comparison of the RNAseq methods ==

The BioMicro Center has done testing in head-to-head competitions of the TruSeq, NuGEN v2 and Clontech kits. These data were presented at AGBT 2012 as part of a poster. The authors were: Avanti Shrikumar, Zachary Banks, Manlin Luo, Ryan Sinapius, Paola Favaretto, Jessica Hurt, Chris Burge, and Stuart S. Levine. Selections of the poster are shown below:

{| cellpadding="10"
 |-
 |rowspan="2" width=350| '''Summary''' <BR>
Sequencing of the transcriptome (RNAseq) has become an increasingly important tool in the molecular biology toolkit and is rapidly replacing microarrays as the primary method for determining genome-wide expression levels. Several vendors have created pre-packaged kits for creating RNAseq libraries for Illumina sequencing. These kits differ significantly in methodology and in the amounts of input required. Here we provide a head to head test of five different RNAseq kits in a core setting. The kits were evaluated on two experimental samples with similar expression patterns, murine embryonic stem cells and the same cells with a single factor knocked down by RNAi, to determine the sensitivity of each method. Each kit was additionally evaluated across  three different concentrations of RNA input. We found that the different methodologies show different RPKM levels for each transcript and also vary in their technical reproducibility. The different methods resulted in small but largely distinct lists of differentially expressed genes that we compared to genes with known expression changes
 |[[Image:AGBT12_1.png|thumb|center|300px]] <BR>
 |align="center"|
 {|Border=1 style="text-align: center; width:150px;"
  !Kit Type !! Home Made (Burge Lab) !! Illumina TruSeq !! NuGen v1 !! NuGen v2 !! Clontech SMARTer
  |-
  |1ug || X || X ||  ||  ||  
  |-
  |100ng || || X ||  || X ||
  |-
  |10ng ||  || X* || X || X || X
  |-
  |1ng || || || X* || X || X
  |-
  |0.1ng || || || || X* || X
 |}
|-
|'''Experimental Design:''' ES cells were transfected with siRNA targeting a splicing factor or a control siRNA. The tested splicing factor normally blocks inclusion of specific exons that had been previously identified by RT-PCR. Reduction of the splicing factor’s levels should lead to an increase in the amounts of these specific transcripts. RNA was collected from the cells and analyzed by RNA-seq. Samples were sequenced to a depth of at least 7.5m reads of 40nt length on either a GAIIx or a HiSeq2000
|'''Testing Matrix:''' A single sample of control and splicing factor knockdown samples were serially diluted. Aliquots of the diluted samples were tested against the 5 methods to determine the ability of each kit to identify differentially expressed genes in a biologically challenging situation as well as to identify their sensitivity to different input amounts and quantify the amount of technical variation. * indicates the amount tested was below the minimum recommended input.
|-
| [[Image:AGBT12_2.png|thumb|center|300px]]
| [[Image:AGBT12_3.png|thumb|center|200px]]
| [[Image:AGBT12_4.png|thumb|center|300px]]
|-
|'''Fraction of reads aligned:''' Reads aligning once or multiple times to the mouse genome (mm9) are shown. Reads were aligned with Bowtie. NuGen samples show an increase in non-aligned reads at low input amounts.
|'''3’ Bias:''' 
Read densities were calculated along the exonic portions of each transcript. Transcripts were grouped by length and plotted as metagenes with both 5’ and 3’ ends locked. Clear 3’ bias can be observed in the samples processed using the Clontech kit.  
|'''Variation in RPKM score across the transcript.'''(right) The evenness of coverage within the transcript was measured by comparing the read density in largest two exons. Only exons greater than 200nt and genes with PRKMs over 10 were included in this analysis (n = 1,876 genes including 3,676 exons) (left) the average number of exons in the above  data set with very low coverage (<5 reads) is shown. 
|}

== Choosing a read length and read depth ==

== RNAseq References ==
A few notable references. Please feel free to add more:

<biblio>
#Paper1 pmid=18978772
#Paper2 pmid=19015660
#Paper3 pmid=18550803
#Paper4 pmid=20711195
<\biblio>