RNA Extraction from Yeast: Bead Mill Protocol | Transcriptomics Lab

cDNA synthesis is the process of converting RNA into complementary DNA (cDNA) using reverse transcriptase enzyme, which enables researchers to study gene expression through quantitative PCR (qPCR); the process involves measuring RNA concentration using spectrophotometry (A260/A280 ratio), calculating appropriate RNA input amounts, preparing a master mix with dNTPs, buffer, random primers, and reverse transcriptase, and performing thermal cycling to synthesize cDNA from the RNA template.

cDNA synthesis is the process of converting RNA into complementary DNA (cDNA) using reverse transcriptase enzyme, which is essential for quantitative PCR analysis since PCR requires DNA templates; the process involves measuring RNA concentration using spectrophotometry, calculating appropriate RNA amounts, preparing a master mix with nuclease-free water, RT buffer, random primers, and dNTPs, adding RNA samples, and performing thermal cycling to synthesize cDNA.

Quantitative RT-PCR (qRT-PCR) is the gold standard for validating gene expression levels from high-throughput data by measuring the inverse relationship between cycle threshold (CT) values and starting cDNA amount, requiring careful assay design spanning exon-exon junctions, validation of PCR efficiency through serial dilution experiments, selection of stable housekeeping genes for normalization, and systematic data analysis using ΔΔCT methodology to calculate fold changes between experimental conditions.

Reverse transcription converts RNA into complementary DNA (cDNA) using reverse transcriptase enzyme, naturally found in retroviruses like HIV and Hepatitis B. The enzyme has three activities: RNA-dependent DNA polymerase, RNase H, and DNA-dependent DNA polymerase. Three primer types are used: oligo(dT), random hexamers, and gene-specific primers. Reverse transcription requires specific buffer conditions and is performed at 37-54°C. After cDNA synthesis, the reaction is heated to 95°C to inactivate the enzyme. Two protocols exist: one-step (combining reverse transcription and PCR in a single tube) and two-step (performing each reaction separately). PCR in real-time (qPCR) is a quantitative technique that monitors DNA amplification during each cycle using two temperature cycles: denaturation at 94-95°C and annealing/extension at 60°C.

cDNA (complementary DNA) is synthesized from an RNA template through reverse transcription, a process that uses reverse transcriptase enzyme along with primers such as oligo d(T)-primers or random hexamers to convert RNA into double-stranded DNA; this technique is essential for gene expression analysis, including quantitative PCR and microarray studies, and allows researchers to study specific gene sequences by creating DNA copies from RNA templates.

RNA library preparation involves multiple steps: (1) RNA extraction and quality assessment (RIN > 8 recommended); (2) DNA removal from RNA samples; (3) mRNA enrichment (via poly-A selection using magnetic beads) or ribosomal RNA depletion; (4) RNA fragmentation; (5) Reverse transcription to cDNA; (6) Adapter ligation; (7) PCR amplification; (8) Library quality control. Each sample becomes a separate library, with barcodes allowing multiple samples to be sequenced together (multiplexing).

RNA sequencing enables analysis of gene expression by identifying active genes in a sample. The primary target is messenger RNAs (mRNAs), which constitute only a small fraction of total RNA, with ribosomal RNA dominating the cellular RNA pool. Effective library preparation requires three essential steps: (1) mRNA enrichment/depletion to isolate target transcripts, (2) reverse transcription converting RNA to complementary DNA (cDNA) since sequencers require DNA input, and (3) adapter addition for library completion. Two main enrichment strategies exist: poly-A capture using oligo-dT beads that bind mRNA tails, or rRNA depletion using biotinylated probes followed by magnetic bead removal. After enrichment, RNA fragmentation creates appropriately sized fragments for sequencing, followed by reverse transcription and second strand synthesis to generate double-stranded DNA libraries ready for amplification and sequencing.

The complete RNA-seq library preparation workflow includes RNA isolation, DNAse treatment to remove genomic DNA, cDNA synthesis, fragmentation, size selection, adapter addition, and amplification. Three main enrichment strategies exist: poly(A) selection captures polyadenylated mRNA, rRNA depletion removes abundant ribosomal RNAs, and exome capture targets known exonic sequences. Each strategy has distinct advantages and limitations depending on research objectives and sample characteristics.

RNA-seq requires systematic experimental design matching research objectives: poly(A) selection or rRNA depletion for mRNA enrichment, with the latter mandatory for prokaryotes lacking poly(A) tails. The workflow progresses through library preparation (RNA quality assessment, cDNA conversion, fragmentation, adapter ligation), bridge PCR for cluster generation, and sequencing by synthesis using fluorescently labeled nucleotides. Quality control employs Qubit, Agilent Bioanalyzer, and qPCR to verify library integrity. Data analysis transforms raw FASTQ files through alignment to reference genomes, followed by critical normalization using FPKM to account for unequal sequencing depth and varying gene lengths. Interpretive techniques include heat maps for visualizing expression changes, PCA for identifying correlated gene expression patterns, and functional annotation to reveal enriched biological pathways among differentially expressed genes.

RNA-seq library preparation involves mRNA enrichment (poly-A selection captures ~70% of transcriptome), fragmentation, cDNA synthesis using reverse transcriptase and random hexamers, end repair, A-tailing, adapter ligation, and PCR amplification. Unique Molecular Identifiers (UMIs) correct for PCR amplification bias by collapsing reads from the same original molecule. Experimental design must account for batch effects from different sequencing runs, kits, or personnel. Samples should be randomized across conditions and processed in the same batch when possible. Technical and biological replicates are essential for capturing natural variability and identifying technical artifacts.

The differential expression pipeline configures DESeq2 with factor names and levels representing experimental conditions, then runs StringTie again with merged transcripts as reference to prepare gene counts. DESeq2 performs statistical testing to identify significantly changing genes, generating diagnostic plots including PCA, dispersion estimates, and MA plots showing log2 fold changes versus expression levels. Results tables contain gene IDs, base means, log2 fold changes (positive = upregulated, negative = downregulated), p-values, and adjusted p-values. AnnotateDESeq2 converts gene IDs to human-readable names using merged transcript files. Filter results by removing genes with 'NA' gene names and selecting significant adjusted p-values (≤0.05), then separate upregulated and downregulated genes for focused analysis.

This section presents the complete DESeq2 statistical workflow for differential expression analysis. The pipeline includes: (1) Setting reference levels for contrasts using factor() with explicit level ordering, (2) Filtering low-count genes using rowSums(counts) >= 5 to remove unreliable estimates, (3) Running DESeq() to perform the core statistical analysis fitting a generalized linear model to count data, and (4) Extracting results using results() to obtain p-values, adjusted p-values, and log2 fold change estimates. The results object contains essential columns including baseMean, log2FoldChange, lfcSE, pvalue, and padj for downstream analysis and interpretation.

Differential gene expression analysis using RNA-seq data requires specialized normalization methods (CPM, RPKM/FPKM, TPM) and statistical models (negative binomial distribution) to account for sequencing depth, gene length, and biological variation. DESeq2 and EdgeR are the two primary R-based tools that implement these approaches, using housekeeping genes for normalization and controlling false discovery rate (FDR) through multiple testing correction. The analysis pipeline includes quality control, alignment with splice-aware tools (STAR, HISAT2), counting, normalization, statistical testing, and visualization through volcano plots and heatmaps.

DESeq2 performs differential gene expression analysis by comparing gene quantification files between experimental conditions. Users specify factor names and levels, select appropriate quantification files, and provide transcript-to-gene annotation. The analysis generates PCA plots showing sample clustering and heatmaps displaying expression patterns. Results are filtered using p-value thresholds (<0.05) and fold change cutoffs (>2-fold) to identify significantly regulated genes for downstream functional analysis.

Dispersion measures variability between biological replicates and is essential for accurate statistical testing. DESeq2 estimates dispersion through three steps: (1) Maximum likelihood estimation for each gene's dispersion, (2) Fitting a curve to gene-wise dispersions against mean expression, (3) Shrinking individual estimates toward the fitted curve. Genes with extremely high dispersion are not shrunk, indicating they have additional unexplained variability beyond biological/technical factors.