A good source of information for this part is RNA-seqlopedia.

Considerations when obtaining your RNA.

The first step in a transcriptomic experiment is to obtain the RNA. After isolating total RNA from cells, one can directly sequence it. Nonetheless, the majority of the RNA in a cell is ribosomal RNA, which may need to be removed using specific kits. Moreover, total RNA also contains unprocessed immature transcripts and RNA targeted for degradation (at different stages of processing).

Therefore, unless one is interested in non-coding RNAs or other aspects related to transcription, it is usually better to apply protocols that extract the mature mRNAs (usually through the PolyA tails). Since most people are interested in coding-genes, it is more common to use mRNA-specific protocols.

Some protocols can also keep strand information. In this case, the reads have the same (or the reverse) strand as the transcribed RNA. This is particularly relevant when sequencing total RNA, noticeably to distinguish real transcripts from transcriptional activity resulting from stalled promoters or enhancers. It can also be useful to distinguish between overlapping genes.

Finally, we also need to consider the amount of material available. Are we dealing with samples with a lot of RNA (eg. cell cultures), or short amounts (eg. small tissue samples, single-cell) that are prone to amplification artifacts and presence of contaminant sequences?

The High Throughput Sequencing Workflow

Sanger sequencing brought about a technological revolution, as it allowed to directly read DNA molecules with relative ease and affordability. The Human Genome Project motivated further progress, leading to automated DNA-sequencing machines capable of sequencing up to 384 samples in a single batch using capillary electrophoresis.

At the moment, the high throughput sequencing technology most often used (by far) is Illumina. Similarly to the Sanger method, it is also based on the addition of nucleotides specifically modified to block DNA strand elongation, where each nucleotide is marked with a different color. Unlike the Sanger method, where a single DNA molecule is “read” at a time, modern illumina machines allow reading up to millions of DNA molecules simultaneously.

The following links are a good source of information regarding this sequencing technology:

• Illumina Sequencing by Synthesis.
• Elaine Mardis talk on NGS.

Commmon steps in most high throughput sequencing workflows:

• Extraction and purification of the DNA template (even RNA must usually be converted to cDNA)

• Fragmentation of the DNA template (into a size range that can be accommodated by the machine)

• Attachment of sequencing tags (to enable reading by the machine)

• Amplification of signal (usually trough PCR, often already in the machine)

• Reading of signal and conversion into nucleotide bases

Many sequencing machines can read both ends of a fragment. This is called paired-end sequencing.

When sending your samples to a sequencing facility, these are the most frequent parameters to consider:

• Single versus Paired-end sequencing

• Read Length

• Coverage (number of reads)

The following links are a good source of information regarding illumina sequencing technology:

• Illumina Sequencing by Synthesis.
• Elaine Mardis talk on NGS.

Designing your experiment for differential expression using RNAseq.

Longer read length, paired-end sequencing and strand-specific library preparation are particularly relevant to reveal gene structure. For example, on a non-model organism for which there is no genome sequenced, or the genes are poorly annotated. They are also relevant when alterative splicing is a factor to take into consideration. Discovering gene structure is a complex process and it would be the subject of an entire course on its own.

For this course, we will focus on the analysis of differential gene expression between conditions, on organisms for which gene annotation is available. Under these conditions, long reads, paired-end, and stranded library preparation methods are not as important. Therefore, for this type of experiments, we can safely go for the cheaper single-end sequencing and shorter read lengths (eg. 50bp or 76bp).

To infer genes differentially expressed between conditions, we need to obtain accurate measures of gene expression variance between the conditions. For this, we need replicates containing as much of the expected biological variance as possible. Chosing the number of replicates and depth of sequencing (number of reads) depends on the experiment. For highly controlled conditions (such as cell cultures), 2-3 replicates could be enough. In terms of coverage, 10-40 million reads should be enough to capture most “reasonably” expressed genes. Nonetheless, to be able to more accurately estimate how much is needed, one should always start from small pilot datasets.

At IGC we mainly use two library preparation methods (both unstranded): Smart-seq and QuantSeq. Smart-seq sequences full length cDNAs, while QuantSeq only sequences the ends of the transcripts (and thus requires less reads per sample). Smart-Seq can take samples with very low numbers of cells, including single-cell experiments which usually require less reads (5-10M), and many replicates. Specific analysis techniques are necessary for samples with very low cell numbers, which we will not be able to cover in this course.

2 - List steps in the analysis of RNAseq differential expression experiments

Steps in the analysis of RNA-Seq:

• QC of Raw Data; (Page 3)

• Preprocessing of Raw Data (if needed); (Page 4)

• Alignment of “clean” reads to reference genome (Page 5)

• QC of Aligments (Page 6)

• Generate table of counts of genes/transcripts (Page 7)

• Differential Analysis tests (Page 8)

• Post-analysis: Functional Enrichment (Page 9)

Back

Back to main page.