Module 1: Introduction to Next-Generation Sequencing & RNA-Seq Fundamentals
Module 1: Introduction to Next-Generation Sequencing & RNA-Seq Fundamentals
Core concepts of high-throughput transcriptomics: Illumina Sequencing-by-Synthesis (SBS) chemistry, flowcell cluster generation, paired-end vs. single-end reads, and dynamic range of transcriptome quantification.
1. Introduction to Next-Generation Sequencing & RNA-Seq
RNA Sequencing (RNA-Seq) is the gold-standard high-throughput technology used to profile, quantify, and characterize the complete transcriptome of any biological sample at single-nucleotide resolution. Unlike legacy hybridization arrays (microarrays) that required pre-designed probe sets and suffered from narrow dynamic ranges, RNA-seq offers unbiased discovery of novel transcripts, alternative splicing isoforms, gene fusions, and single nucleotide polymorphisms (SNPs) within transcribed regions.
Key Advantage of RNA-Seq over Microarrays:
- Unbounded Dynamic Range: Quantifies low-abundance transcription factors alongside ultra-high abundance ribosomal or housekeeping mRNAs without sensor saturation.
- No Prior Genomic Sequence Required: Capable of de novo transcriptome assembly in non-model organisms or infectious pathogens.
- Structural Transcriptome Resolution: Directly identifies exon skipping, intron retention, alternative 5'/3' splice sites, and chimeric gene fusions.
How Illumina Sequencing-by-Synthesis (SBS) Works
The overwhelming majority of short-read transcriptomic datasets are generated using Illumina sequencing-by-synthesis chemistry. Double-stranded cDNA fragments flanked by platform adapters (P5/P7) hybridize to an automated glass flow cell coated with a dense lawn of complementary oligonucleotides:
- Bridge Amplification / Exclusion Amplification (ExAmp): Generates clonal clusters containing ∼1,000 identical copies of each cDNA molecule to boost optical signal intensity.
- Fluorescent Reversible Terminator Incorporation: In each cycle, all four dNTPs (each labeled with a distinct cleavable fluorophore and reversible 3'-OH block) compete for incorporation by DNA polymerase.
- High-Resolution Laser Imaging: The flow cell surface is scanned to detect fluorescent emission wavelengths for base identification (4-channel, 2-channel on NovaSeq/NextSeq, or 1-channel on iSeq).
- Cleavage & De-blocking: Chemical cleavage removes the fluorophore and regenerates the 3'-OH group, permitting the next sequential nucleotide addition.
Hands-On Summary Table: Short-Read vs. Long-Read Transcriptomics
| Feature | Illumina Short-Read RNA-Seq | PacBio Iso-Seq / Oxford Nanopore Direct RNA |
|---|---|---|
| Read Length | 50 – 150 bp (Paired-End 2x150) | Full-length transcripts (1,000 – 15,000+ bp) |
| Throughput & Cost | Ultra-high throughput (billions of reads), low cost/sample | Moderate throughput, higher cost per read |
| Primary Application | Accurate quantitative differential gene expression (DGE) | Full isoform discovery, novel splice junctions, direct RNA modifications |
| Computational Challenge | Spliced alignment across short exon junctions | Higher raw indel rates, transcript clustering |
Common Pitfalls & Troubleshooting Checklist
- Java Heap Exhaustion (GATK / Picard): If processes terminate with
java.lang.OutOfMemoryError, set your environment variable explicitly:export _JAVA_OPTIONS="-Xmx16g -Djava.io.tmpdir=/tmp". - Index Timestamp Mismatch: If tools complain that a BAM index is missing or out of date, regenerate with
samtools index -@ 4 sample.bam. - Conda Channel Conflicts: Avoid ABI incompatibilities between bioconda and conda-forge by setting strict priority:
conda config --set channel_priority strict. - Chromosome Naming Contamination: Always verify whether reference FASTA headers use
chr1or1; divergent naming causes silent alignment drops.