Microbial pathway analysis tool or pipeline
2
0
Entering edit mode
14 months ago
ssko ▴ 20

Hello everyone!

I had 16S rRNA sequence data obtained from 21 samples which are separated into two groups. I performed functional community prediction with these 16S sequences using MicFunPred tool. Then, I obtained an ortholog table (KO table) like OTU table. Using this table, I did differential abundance analysis to see significantly abundant features between two groups. When I use fdr correction I got a few statistically significant KOs.

I was tasked to find pathways from those significant KOs. I searched many articles and they mostly used Ingenuity pathway analysis tool or other commercial tools. Does anyone know of free, open source tools for microbial pathway analysis? or I am open to any advice/idea on finding pathways with the KOs information I have obtained.

Thank you all!

pathways 16SrRNA microbiome • 1.1k views
ADD COMMENT
0
Entering edit mode

Not exactly my field of expertise, but maybe you can get some inspiration for tools from nf-core ampliseq pipeline. Possibly some of the tools from the funcscan pipeline will also run with the input that you already have?

ADD REPLY
1
Entering edit mode
14 months ago

Try paprica or PICRUSt2. Pathway prediction based on 16S rDNA metataxonomic reads is exactly what they do. I hope I understand your task correctly.

Also, there are some analogous tools listed and compared in “Inferring microbiota functions from taxonomic genes: a review”.

What concerns paprica, you can find some useful tutorials here: https://www.polarmicrobes.org/category/paprica/.

ADD COMMENT
0
Entering edit mode
3 months ago

I think you can simply load your KO's list up there in the KEGG website KEGG_mapper. Than you should see KEGG pathways drawn. I suggest you to read carefully instructions becouse i used it some times ago.

I have a question for you, i'm using MicFunPred tool too but i was not able to reproduce the last graph shown on the README file, where you?

ADD COMMENT
1
Entering edit mode

I couldn't reproduce the last graph either, so I created a similar one using the KO_taxon_contrib.tsv file for the relevant path using R.

ADD REPLY

Login before adding your answer.

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