Lesson · 40 min · Free
RNA-seq Expression Analysis
RNA-seq Expression Analysis body { font-family: sans-serif; line-height: 1.6; color: #333; } h1, h2 { color: #0056b3; } pre { background-color: #f4f4f4; padding: 10px; border: 1px solid #ddd; overflow-x: auto; margin-bot
RNA-seq Expression Analysis
Welcome to this lesson on RNA-seq Expression Analysis, a fundamental technique in modern genomics for quantifying gene expression levels. RNA sequencing (RNA-seq) has revolutionized our understanding of transcriptomes, allowing us to identify and quantify RNAs present in a biological sample at a given time. This approach provides a comprehensive view of gene activity, enabling insights into disease mechanisms, drug responses, and developmental processes. At its core, RNA-seq involves sequencing fragmented RNA molecules, mapping these reads back to a reference genome or transcriptome, and then counting the number of reads that align to each gene. The more reads that map to a particular gene, the higher its expression level is presumed to be. This quantitative measure of gene expression is crucial for differential expression analysis, where we compare gene expression between different conditions (e.g., treated vs. control, healthy vs. diseased).
Key Steps in RNA-seq Expression Analysis
The workflow for RNA-seq expression analysis typically involves several critical steps, each with its own computational challenges and considerations: Read Quality Control: Raw sequencing reads often contain adapters, low-quality bases, and other artifacts that can negatively impact downstream analysis. Tools like FastQC are used to assess read quality, and trimming software (e.g., Trimmomatic, Cutadapt) is employed to remove unwanted sequences and low-quality regions. Read Alignment/Mapping: Cleaned reads are then aligned to a reference genome or transcriptome. This step identifies the genomic location from which each read originated. Popular aligners include STAR (Spliced Transcripts Alignment to a Reference) and HISAT2, which are designed to handle spliced reads characteristic of eukaryotic RNA. Read Quantification: After alignment, the next step is to count how many reads map to each gene or transcript. This is often done using tools like featureCounts or HTSeq-count, which assign aligned reads to genomic features (e.g., exons, genes) based on an annotation file (GTF/GFF). The output is typically a count matrix, where rows represent genes and columns represent samples, with values indicating the raw read counts. Normalization: Raw read counts cannot be directly compared across samples due to variations in sequencing depth, gene length, and RNA composition. Normalization methods (e.g., TMM, DESeq2's median-of-ratios) adjust these counts to make them comparable. This step is crucial for accurate differential expression analysis. Differential Expression Analysis: This is often the primary goal of RNA-seq. Statistical methods (e.g., DESeq2, edgeR, limma-voom) are used to identify genes that are significantly up- or down-regulated between different experimental conditions. These tools model the count data, often assuming a negative binomial distribution, and perform hypothesis testing to identify differentially expressed genes (DEGs). Downstream Functional Analysis: Once DEGs are identified, further analysis is performed to interpret their biological significance. This can include gene ontology (GO) enrichment analysis, pathway analysis (e.g., KEGG), and network analysis to understand the biological processes, molecular functions, and cellular components associated with the differentially expressed genes.
Example: Read Alignment with STAR
STAR is a highly efficient and accurate aligner for RNA-seq reads. Before aligning, you need to build a genome index. # Build STAR genome index (run once per reference genome) STAR --runThreadN 8 --runMode genomeGenerate \ --genomeDir /path/to/STAR_index \ --genomeFastaFiles /path/to/genome.fa \ --sjdbGTFfile /path/to/annotations.gtf \ --sjdbOverhang 100 # Align reads for a single sample STAR --runThreadN 8 --genomeDir /path/to/STAR_index \ --readFilesIn /path/to/sample_R1.fastq.gz /path/to/sample_R2.fastq.gz \ --outFileNamePrefix /path/to/output/sample_ \ --outSAMtype BAM SortedByCoordinate \ --outReadsUnmapped Fastx \ --outSAMstrandField intronMotif
Example: Gene Quantification with featureCounts
After alignment, featureCounts can be used to count reads mapping to genes based on a GTF annotation file. # Count reads for multiple BAM files featureCounts -p -t exon -g gene_id \ -a /path/to/annotations.gtf \ -o /path/to/output/counts.txt \ /path/to/sample1_Aligned.sortedByCoord.out.bam \ /path/to/sample2_Aligned.sortedByCoord.out.bam \ /path/to/sample3_Aligned.sortedByCoord.out.bam The -p flag indicates paired-end reads, -t exon specifies that features of type 'exon' should be counted, and -g gene_id specifies that counts should be summarized by 'gene_id'.
Key Takeaways
RNA-seq quantifies gene expression by sequencing and counting RNA molecules. The analysis workflow includes quality control, alignment, quantification, normalization, and differential expression analysis. STAR and HISAT2 are common aligners, while featureCounts and HTSeq-count are used for quantification. DESeq2 and edgeR are widely used for statistical differential expression analysis. Downstream functional analysis helps interpret biological meaning of differentially expressed genes.
Practice Exercise
Imagine you are analyzing RNA-seq data from two groups of patient samples: one treated with a novel anti-cancer drug and a control group receiving a placebo. After performing alignment and quantification, you have a count matrix. Describe the next two essential bioinformatics steps you would take to identify genes whose expression is significantly altered by the drug treatment, and briefly explain the purpose of each step. What type of biological insights would you hope to gain from the final list of genes?
Watch the full lesson — free
This topic is part of Bioinformatics & Computational Genomics, a complete AI-narrated video course. Press play once and watch the entire lecture like a movie.
Start the course free →