Differential exon usage analysis using DEXSeq for mRNA and lncRNA separately
0
0
Entering edit mode
Anubhav • 0
@Anubhav-24483
Last seen 22 months ago

I want to do a Differential Exon Usage (DEU) analysis of both mRNA and lncRNA from the same RNAseq samples using DEXSeq. The raw reads were aligned using STAR and GENCODE human v33 as the annotation file. For the DEU, I would like to use a slightly older version (GENCODE v30) for counting reads, since it contains lncRNAs separated into various classes (lincRNA, sense_overlapping, antisense, etc.). Now, I can think of two strategies for the analysis :

  1. Separate the original comprehensive gtf file into subfiles for each category (mRNA and lncRNA) using the "gene type" argument from the 9th column and "exon" from 3rd column of the GTF file, followed by counting of each GTF file using featureCounts (or HTSeq) and subsequent input of the count matrices into R for DEXSeq analysis.

  2. Run featureCounts (or HTSeq) on the full GTF file, read it into R, do DEXSeq analysis, and then separate the results into mRNA or lncRNA by mapping the ENSEMBLID to the "gene type" from the GTf file.

Which one of these approaches would be the correct way to go about it?

DEXSeq • 499 views
ADD COMMENT
1
Entering edit mode

Hi Anubhav,

Thanks for your interest. Both approaches are correct, I'd just try both and have a look at the resuts. I expect the results to be very similar

Alejandro

ADD REPLY

Login before adding your answer.

Traffic: 426 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6