Hi there,
I wanted to perform DE of lncRNAs and mRNAs from the RNA-seq data using DEseq2. I am wondering if it is a good idea to have separate annotation files for each and then perform the differential expression (DE) or is it OK to perform DE analysis using the annotation files that contains both lncRNAs as well as protein-coding mRNAs and then filter based off the gene-biotype?
Your response or insights are highly appreciated. Thanks very much!
Regards, Bhumi
2 answers
The only reason to process them separately is if the lncRNAs are so lowly expressed that they keep getting excluded due to independent filtering. If that does happen to be the case, then perform everything up to and including computation of the size factors with the lncRNAs mixed with the other genes (this is (1) to ensure that the size factors remain constant and (2) to ensure robustness). Thereafter you can split and test separately.
This post will give you some idea,
A: Any One please provide protocol for Analysing long noncoding RNA illumina NGS da
I would recommend you to perform DE analysis separately (since protein coding expression might create biasness due to its abundance).
Log in to answer this question.
-
About bowtie for sRNA mapping
written by r00628112 •
-
A recent approach to identify lncRNAs from a large list of differentially expressed genes
written by diego1530 •
-
DESeq2 multiple group comparisons
written by Bhumi •
-
Small RNA sequencing adapter trimming and alignment query
written by Bhumi •
-
Filtering out lowly expressed genes
written by Bhumi •
-
Reference genome index for LncRNAs (from RNA-seq data)
written by SHN •
-
Classification of samples based on gene expression profiles
written by Bhumi •
-
How to perform lncRNA differential gene expression analysis (unstranded)?
written by mmitra •
-
Filtering in RNA-Seq to get Differentially Expressed lncRNA
written by conor93
-
DE analysis of lncRNA from C elegans
written by Assa Yeroslaviz
My preference is to keep them separate and to conduct the analysis separately. Protein coding mRNAs are typically expressed at higher levels than ncRNAs, with many ncRNAs being expressed at very low values (if at all, in your tissue of interest). For that reason, I feel that normalising them together would in fact bias the results, but I don't have data to back up this hypothesis.
I would be interested to hear other ideas.
It will be a good idea to process them separately starting from alignment to DE. The following research article suggests a similar approach.
http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0150628
Thanks all for your responses. I do agree that performing lncRNA and mRNA analysis separately is a better approach compared to processing them altogether.