Differential exon usage analysis using DEXSeq for mRNA and lncRNA separately
Entering edit mode
Anubhav • 0
Last seen 3.3 years 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 • 1.0k views
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



Login before adding your answer.

Traffic: 698 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6