Question: How to perform lncRNA differential gene expression analysis (unstranded)?
gravatar for mmitra
21 months ago by
Los Angeles, United States
mmitra30 wrote:

Hi all, I am trying to perform differential gene expression analysis between two conditions. Primarily, I am interested in the lncRNAs that are differentially expressed. My workflow is: Salmon--> Txtimport--> DEseq2. I have done the analysis using two different annotations (both Gencode v22): comprehensive gene annotation (all the genes including lncRNAs) and long non coding RNA gene annotation (lncRNAs only). The RNA-seq library preparation method is of unstranded type, so it would be a problem to unequivocally determine the expression of lncRNAs that are antisense to protein-coding genes. I have two questions:

  1. If I am primarily interested in lncRNAs that are differentially expressed, would it matter if I do the analysis using comprehensive gene annotation (all the genes) and long non coding RNA gene annotation (lncRNA genes only)?

  2. As the analysis is for unstranded, shall I just not consider the antisense lncRNA genes after the differential expression analysis? Would it be better to start the analysis with an annotation that contains only long intergenic RNAs (lincRNAs) that do not overlap with any protein coding genes? I am not sure how to modify the gencode annotation to do the latter.

I would appreciate any help. Thanks in advance!

ADD COMMENTlink modified 21 months ago by JC11k • written 21 months ago by mmitra30
gravatar for JC
21 months ago by
JC11k wrote:
  1. I consider that you need to use all genes, otherwise a mapping bias could exist.
  2. That is a tricky question, definitively you need to run some run tests and check the outputs to see if makes sense and if possible, to validate in the lab.
ADD COMMENTlink written 21 months ago by JC11k

Thanks for your comments. Very helpful!

ADD REPLYlink written 21 months ago by mmitra30
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1301 users visited in the last hour