Question: DEXSeq on FeatureCounts output
0
gravatar for Vivek.b
4.3 years ago by
Vivek.b100
Germany
Vivek.b100 wrote:

Dear DEXSeq and FeatureCounts team

I was wondering if it is possible to better integrate DEXSeq with featureCounts (either the RSubread or the command line version) output, so that users don't have to modify the FeatureCounts output to convert into DEXSeq acceptable format. I have been doing this by adding Exon IDs (E1,E2..) to the gene names followed by  inputting it to DEXSeqDataSet command. However still I found that some of the plots didn't work properly (probably that's not the best way to do it). It would be great to have something like DEXSeqDataSetFromFeatureCounts function, or probably getting RSubread output in a SummarizedExperiment format. However, I would appreciate any advice on this issue.

Thanks

dexseq featurecounts • 1.7k views
ADD COMMENTlink modified 4.3 years ago • written 4.3 years ago by Vivek.b100

Hi Vivek.b, 

Would be happy to implement this integration in DEXSeq. Do you have a prototype code from the featureCounts part? Does it follow the same "counting rules" described in the DEXSeq manuscript/vignette?

Alejandro

ADD REPLYlink written 4.2 years ago by Alejandro Reyes1.7k

Hi Alejandro

Thanks for your reply (and sorry for my late reply). I had to make a couple of changes to make it work:

  • DEXSeq-flattened GFF converted to GTF for featurecounts.
  • Adding exon IDs to featurecounts output.
  • Made a DEXSeqDataSetFromFeatureCounts function to read the converted output into dexSeq.

After running DEXSeq, the output from featureCounts (if we also count reads overlapping more than one feature), is very similar to that from DEXSeq_count. Just that featureCounts is much faster. 

Now I can also visualize my output properly using plotDEXSeq and DEXSeqHTML functions. So things are working for me.

I have put my scripts on Github here.. I hope this can be integrated into DEXSeq.

ADD REPLYlink modified 4.1 years ago • written 4.1 years ago by Vivek.b100
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 16.09
Traffic: 168 users visited in the last hour