No history yet

RNA-seq Data Preparation

Refining Raw SCD Transcriptomes

The first step in analysing vaso-occlusive crisis (VOC) transcriptomics is rigorous quality control. Raw data from sequencers comes in FASTQ format, which includes the nucleotide sequence and an associated quality score for each base. We use tools like FastQC to get a summary of this quality. This initial check looks for biases, adapter contamination, and overall read quality across the length of the sequences. For Sickle Cell Disease (SCD) datasets, this step is particularly important because blood samples can be challenging to work with, leading to variable RNA quality.

After assessing quality, we perform trimming. This process removes adapter sequences left over from the library preparation stage and cuts off low-quality bases, typically from the end of reads. Tools like Trimmomatic or Cutadapt are standard for this task. The goal is to clean the data, ensuring that only high-quality biological sequences are passed to the next stage. In SCD studies, a major challenge is the overwhelming presence of RNA, which can mask the signals from other cell types involved in a crisis. Careful QC helps confirm the data's integrity before we tackle this biological complexity.

Alignment and Quantification

With clean reads, the next step is alignment. We map the short RNA sequences to a reference genome to determine their origin. For RNA-seq data, this requires a splice-aware aligner, as reads may span introns that have been spliced out of mature mRNA. The most common tools for this are STAR and HISAT2. They are designed to handle spliced alignments efficiently and accurately. For SCD, which involves complex genetic variants, using a comprehensive genome reference with appropriate annotations is key to correctly mapping reads, especially those from erythroid precursor cells.

Once reads are aligned, we perform quantification to create a gene count matrix. This matrix tabulates the number of reads mapped to each gene for every sample. Tools like featureCounts or HTSeq-count are typically used. The output is a raw count table, which is the foundation for all subsequent differential expression analysis. This raw data isn't directly comparable between samples due to variations in sequencing depth and library composition, which is why the next step is critical.

The final step before analysis is normalisation. This process adjusts the raw counts to account for differences in library size (the total number of reads per sample) and other technical variations. Without it, you might incorrectly conclude a gene is highly expressed in one sample when it only appears to be because that sample was sequenced more deeply.

In longitudinal SCD studies, where samples are collected from the same patient during and between crises, are a significant concern. These are technical variations arising from processing samples at different times or with different reagent batches. They can introduce systematic noise that obscures true biological differences. Normalisation techniques, such as TMM (Trimmed Mean of M-values) used by edgeR or the median of ratios method in DESeq2, are designed to handle this. For high-variance transcriptomes seen in SCD, these methods are more robust than simpler ones like FPKM or TPM.

# --- A Simplified SCD RNA-seq Pipeline ---

# Input: sample1.R1.fastq.gz, sample1.R2.fastq.gz
# Reference: genome_index/, genome.gtf

# 1. Quality Control
fastqc sample1.R1.fastq.gz sample1.R2.fastq.gz

# 2. Trimming with Trimmomatic
trimmomatic PE -threads 4 \
  sample1.R1.fastq.gz sample1.R2.fastq.gz \
  s1.R1.paired.fq.gz s1.R1.unpaired.fq.gz \
  s1.R2.paired.fq.gz s1.R2.unpaired.fq.gz \
  ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36

# 3. Alignment with STAR
STAR --genomeDir genome_index/ --readFilesIn s1.R1.paired.fq.gz s1.R2.paired.fq.gz --readFilesCommand zcat --runThreadN 8 --outFileNamePrefix sample1_

# 4. Quantification with featureCounts
featureCounts -T 4 -p -t exon -g gene_id -a genome.gtf -o counts.txt sample1_Aligned.out.sam

# Output: counts.txt (Gene Count Matrix)

Let's review the key terms from this data preparation pipeline.

Now, test your understanding of the process.

Quiz Questions 1/5

What is the primary purpose of using a tool like FastQC in the initial step of transcriptomics analysis?

Quiz Questions 2/5

Why is a 'splice-aware' aligner, such as STAR, necessary for analysing RNA-seq data?

With a normalised gene count matrix, you are now ready to proceed to differential gene expression analysis, which will help identify the specific genes that change during a vaso-occlusive crisis.