Question: Distinguishing protein transcript encoded within lincRNA in RNA-seq DE analysis
gravatar for sysboolean
8 months ago by
sysboolean40 wrote:

Hi all,

I would appreciate some advice on a problem I have run into. My RNA-seq pipeline is fastp (QC) -> Salmon -> tximeta -> DESeq2 to do gene level DEA. My problem is related to the theme of this article i.e. protein coding genes within lincRNA. I am interested in differentially expressed lincRNA in my study system, using RNA-seq to discover candidate lincRNA that we then individually follow up with functional studies. If a lincRNA gene has both coding and non-coding splice variants (see for example;g=ENSG00000235387;r=9:35909490-35937153), I only see a single gene entry in my Salmon processed DESeq2 data. I don't know if the protein transcript or the lincRNA transcript is being differentially expressed.

Any ideas on how I can distinguish which splice variants are differentially expressed (padj <0.05, |log2FC| > 2) in my RNA-seq analysis ? Would using a splice-aware aligner like STAR followed by DESeq help ? Thanks for any input and sorry if this has already been asked before !

ADD COMMENTlink modified 8 months ago by Kevin Blighe69k • written 8 months ago by sysboolean40
gravatar for Kevin Blighe
8 months ago by
Kevin Blighe69k
Republic of Ireland
Kevin Blighe69k wrote:

I think that a splice-aware aligner followed by DEX-seq or IsoformSwitchAnalyzeR is what you need.


ADD COMMENTlink written 8 months ago by Kevin Blighe69k

Thanks Kevin, I will look into your suggestions !

ADD REPLYlink written 8 months ago by sysboolean40
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: 2574 users visited in the last hour