Question: Distinguishing protein transcript encoded within lincRNA in RNA-seq DE analysis
0
gravatar for sysboolean
8 months ago by
sysboolean40
Denmark
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 https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5520851/ 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 http://www.ensembl.org/Homo_sapiens/Gene/Summary?db=core;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
1
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.

Kevin

ADD COMMENTlink written 8 months ago by Kevin Blighe69k
1

Thanks Kevin, I will look into your suggestions !

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