Getting OTUID rather than Taxonomy from Results function in DESeq2
Entering edit mode
ssingh59 • 0
Last seen 4.1 years ago

Hello All!

I am attempting to do differential abundance analysis using DESeq2 after doing the initial part of my analysis in QIIME 2. I was able to successfully export a .biom file, added sample metadata and taxonomy information, along with a tree and reference sequences using the import_biom function. I then used the code pasted with DESeq2 package version 1.18.1. However, the output always gave me the OTUID rather than the actually taxonomy... I even checked in phyloseq to see if the taxonomy was added correctly to the .biom file using tax_table(biom_file)[1000:1010, 1:7], and it returned taxonomic information... Any guidance is much appreciated!

subset <- subset_samples(biom_file, variable != "NA")
variabledds = phyloseq_to_deseq2(subset, ~ variable)

#a zero-tolerant variant of geometric mean:
gm_mean = function(x, na.rm=TRUE){
  exp(sum(log(x[x > 0]), na.rm=na.rm) / length(x))

geoMeans = apply(counts(variabledds), 1, gm_mean)
variabledds = estimateSizeFactors(variabledds, geoMeans = geoMeans)
variabledds = DESeq(variabledds, fitType="local")

alpha = 0.05

res0 = results(variabledds, contrast = c("variable", "Yes", "No"))
res0 = res0[order(res0$padj, na.last=NA), ]
sigtab0 = res0[(res0$padj < alpha), ]
deseq2 R microbiome phyloseq qiime • 754 views
Entering edit mode
Last seen 1 day ago
United States

This is a phyloseq question more than a DESeq2 question, as DESeq2 is not in control of the handoff of information. DESeq2 just receives rownames and count matrices from upstream software. And we don’t have any support for taxonomy or such structure on the rows of a DESeqDataSet, only rownames or what tabular data can go on the rowRanges.



Login before adding your answer.

Traffic: 225 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