This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DEseq2 input from nf-core/rnaseq

I used the nf-core/rnaseq pipeline in the cluster human genome, and I got some files as the output of Star salmon.

My question is which file I use as input to DESeq2? and where is the full file that contains 58,000 genes?

salmon.merged.gene_counts_length.tsv (29744 genes)
salmon.merged.gene_counts_length_scaled.tsv (29744 genes)
salmon.merged.gene_counts_scaled.tsv (29744 genes)
salmon.merged.gene_counts.tsv (29744 genes)
salmon.merged.gene_temp.tsv (29744 genes)
deseq2

Looks like the salmon reference you used contained 29744 transcripts. There are currently 19396 protein coding genes in human genome as of this month.

Thanks GenoMax. So, what is the correct salmon reference that I should use?

How did you run the pipeline? Which reference did you use?

I ran the pipeline using the nextflow and the reference I used for GRCh38.

From where I can download the latest version of the human reference genome (I need to use wget in linux but I need the link). I need to use it for RNA-seq analysis

1 answer

You should use salmon.merged.gene_counts_length_scaled.tsv and here is the reference.

While feeding to DESeqDataSetFromMatrix() use round() for inputing counts.

Thank you bioinfo8. Also, I am not sure if I use the correct reference genome. So, from where I can download the latest version of the human reference genome (I need to use wget in linux but I need the link). I need to use it for RNA-seq analysis

Log in to answer this question.