Wiki¶

We use the https://nf-co.re/rnaseq/3.13.2 to analyze the RNA seq data (fastq files). The pipeline removes the ribosomal RNA, check the quality of the reads, remove adapter and quality trim, removes genome contaminants, align the reads to the reference genome, sort and index the alignments, mark duplicates and perform quantification.

Data¶

The samples were downloaded from AWS to ilifu in the project folder

/cbio/projects/028/

The rawdata is in

/cbio/projects/028/rawdata/CleanData/

after git clonning the nextflow nfcore/rnaseq pipeline, we used rnaseq/bin/fastq_dir_to_samplesheet.py to create a samplesheet for the pipeline

 ./rnaseq/bin/fastq_dir_to_samplesheet.py /cbio/projects/028/rawdata/CleanData/ samplesheet.csv --strandedness auto --read1_extension "_1.fq.gz" --read2_extension "_2.fq.gz"

NB: We couldn't remove the ribosomal rRNA using the ensembl gtf file. As such we tried the gencode one (see, /cbio/projects/028/scripts/rnaseq24112023.sh).

I also encountered a problem with biomaRt, based on [[[https://github.com/grimbough/biomaRt/issues/89]]] this is due to the incompatibilities between dbplyr, BiocFileCache, and biomaRt as such I downgraded dbplyr to version 2.3.4. Otherwise one would upgrade their BiocManager to version 3.18 if they are using 3.17.