DESeq2 for normalizing metagenomic counts
1
0
Entering edit mode
@ajgasparrini-23549
Last seen 2.0 years ago
United States

I'm interested in normalizing data from a metagenomic experiment. It's similar to a transcriptomic experiment in that I have counts for several thousand genes across several hundred samples. Two key differences that I'm wary of are (a) this data is much sparser than you would encounter in a transcriptomic experiment and (b) much greater inter sample variability than you would expect in a transcriptomic experiment

Would the mean ratio method (i.e. as called by estimateSizeFactors()) be appropriate for normalization in this case? Are there known core housekeeping genes in a metagenome which can be used as controls for size factor estimation?

deseq2 • 518 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 14 hours ago
United States

You can search here with google for some threads where i talk about DESeq2 for metagenomics / microbiome.

I’m not convinced it’s an optimal method for various reasons but check the threads. I know there was one where i linked to a recent study by Levi Waldron et al.

For size factors, poscounts should be better than default, see the man pages. But it’s also not perfect.

ADD COMMENT

Login before adding your answer.

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