Tool for metagenomic functional analysis with 4-factor study design
1
0
Entering edit mode
Aditya • 0
@d153652d
Last seen 11 weeks ago
United States

Hi all,

We've got some raw KEGG ortholog counts from shotgun metagenome sequencing for a somewhat unique (as least as far as I've seen) dataset. There are multiple different study arms, each with a respective treatment group and a control group. These are diet related studies, so the control group in each arm received a certain diet, while the treatment group received the same diet with an additional food of an interest. The control diets were not the same across the different study arms. Sequencing was performed for each sample prior to receiving the intervention (treatment or control), and afterwards. Note that the subjects undergoing the treatment/control interventions are not paired, so a different set of subjects served as the control.

We are primarily interested in finding genes of interest that are differentially expressed within each study arm's treatment group. We want to use the pre-intervention timepoints to control for the initial state of each subject's gene expression, and the control group samples to control for the effect of background diet that was provided. So far we've looked at DESeq2/edgeR for performing this analysis, but have seen both across this forum and the literature that these tools might not be suited for our use case (the same applies to LEFse, ANCOM-BC, etc).

Does anyone have any thoughts on specific pre-existing tools that might be fit for this use case, or would a bespoke analysis be the best route to go down? We're having trouble finding examples of people using this sort of study design, and specifically in the microbiome metagenome space, so I figured I would ask here.

Any input is appreciated - thanks!

Metagenomics DESeq2 metagenomeSeq edgeR KEGG • 122 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

I have no experience with microbiome or metagenomics, but I have reason to think DESeq2 is not an appropriate model for this type of data.

ADD COMMENT

Login before adding your answer.

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