Aim:
Assess the quality of raw datasets
Define quality trimming parameters prior running the complete RNAseq gene profiling pipeline
Run RNA-seq QC check
The pipeline requires preparing at least 2 files:
Metadata file (samplesheet.csv) thatspecifies the name of the samples, location of FASTQ files ('Read 1' and ‘Read 2’), and strandedness (forward, reverse, or auto. Note: auto is used when the strandedness of the data is unknown)
PBS Pro script (launch_nf-core_RNAseq_QC.pbs) with instructions to run the pipeline
Create the metadata file (samplesheet.csv):
Change to the data folder directory:
...
| Code Block |
|---|
cp /work/training/2024/rnaseq/scripts/create_samplesheet_nf-core_RNAseq_SEdata.sh $HOME/workshop/20242/session4_RNAseq/data |
Note: you could replace ‘$HOME/workshop/data’ with “.” A dot indicates ‘current directory’ and will copy the file to the directory where you are currently located
View the content of the script:
| Code Block |
|---|
cat create_samplesheet_nf-core_RNAseq_SEdata.sh |
Example for Single-End data (when only ‘Read 1’ is available):
...
NOTE: modify ‘read1_extension’ as appropriate for your data. For example: _1.fastq.gz or _R1_001.fastq.gz or _R1.fq.gz , etc
Let’s generate the metadata file by running the following command:
| Code Block |
|---|
sh create_samplesheet_nf-core_RNAseq_samplesheetSEdata.sh |
Check the newly created samplesheet.csv file:
| Code Block |
|---|
ls -l
cat samplesheet.cvs |
...
sample,fastq_1,fastq_2,strandedness SRR20622172,/work/training/rnaseq/data/SRR20622172.fastq.gz,,auto SRR20622173,/work/training/rnaseq/data/SRR20622173.fastq.gz,,auto SRR20622174,/work/training/rnaseq/data/SRR20622174.fastq.gz,,auto SRR20622175,/work/training/rnaseq/data/SRR20622175.fastq.gz,,auto SRR20622176,/work/training/rnaseq/data/SRR20622176.fastq.gz,,auto SRR20622177,/work/training/rnaseq/data/SRR20622177.fastq.gz,,auto SRR20622178,/work/training/rnaseq/data/SRR20622178.fastq.gz,,auto SRR20622179,/work/training/rnaseq/data/SRR20622179.fastq.gz,,auto SRR20622180,/work/training/rnaseq/data/SRR20622180.fastq.gz,,auto |
|---|
Copy the PBS Pro script for QC (launch_nf-core_RNAseq_QC.pbs)
Copy and paste the code below to the terminal:
| Code Block |
|---|
cp $HOME/workshop/2024/rnaseq/data/samplesheet.csv $HOME/workshop/2024/rnaseq/runs/run2_QC cp $HOME/workshop/scripts/launch_nf-core_RNAseq_QC.pbs $HOME/workshop/2024/rnaseq/runs/run2_QC cd $HOME/workshop/2024/rnaseq/runs/run2_QC |
Line 1: Copy the samplesheet.csv file to the working directory
Line 2: move to the working directory
Line 3: copy the launch_nf-core_RNAseq_QC.pbs submission script to the working directory
View the content of the launch_nf-core_RNAseq_QC.pbs script:
...
#!/bin/bash -l #PBS -N nfrnaseq_QC #PBS -l select=1:ncpus=2:mem=4gb #PBS -l walltime=24:00:00
#work on current directory cd $PBS_O_WORKDIR
#load java and set up memory settings to run nextflow module load java export NXF_OPTS='-Xms1g -Xmx4g'
nextflow run nf-core/rnaseq \ -profile singularity \ -r 3.14.0 \ --input samplesheet.csv \ --outdir results \ --genome GRCm38-local \ --skip_trimming \ --skip_alignment \ --skip_pseudo_alignment |
|---|
We recommend running the nextflow nf-core/rnaseq pipeline once and then assessing the fastqc results folder to assess if sequence biases are present in the 5'-end and 3'-end ends of the sequences.
Version 3.12.0 allows running the pipeline to do quality assessment only, without any alignment, read counting, or trimming.
Submitting the job
Once you have created the folder for the run, the samplesheet.csv file, and launch.pbs, you are ready to submit the job to the HPC scheduler:
...
Once the pipeline has finished running - Assess the QC report:
NOTE: To proceed, you need to be on QUT’s WiFi network or signed via VPN.
To browse the working folder in the HPC type in the file finder:
...
Assess QC reports (FastQC and MultiQC) to define how many nucleotides should be trimmed from the 5'-end and/or 3-end regions of the FASTQ reads (see Case 3 below).