Question: How to perform lncRNA differential gene expression analysis (unstranded)?
0
gravatar for mmitra
6 months ago by
mmitra30
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 6 months ago by JC8.2k • written 6 months ago by mmitra30
1
gravatar for JC
6 months ago by
JC8.2k
Mexico
JC8.2k 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 6 months ago by JC8.2k

Thanks for your comments. Very helpful!

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

Help
Access

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