Question: NOOB education needed for microbiome functional analysis (specifically PICRUSt2).
0
gravatar for ariel
11 months ago by
ariel110
ariel110 wrote:

I give picrust2 a list of 16s ASVs that map to biota at various taxanomic levels. Picrust2 creates a table of inferred per-sample abundances of KEGG genome orthologs (KO) and one of KEGG enzyme classifications (EC). The main output is a list of per-sample inferred pathway abundances that seems to take both of these into account.

Could I have a basic description of how this works and why? I did RTFM, but I'm inexperienced with pathway analysis and totally green in the context of the microbiome. So, while I can recite the steps picrust2 uses to go from ASVs to pathways, I humbly admit that I don't really understand what I'm doing or what this information ultimately means.

ADD COMMENTlink written 11 months ago by ariel110
1

Which part of the process are you having problem with? Taxa -> KO or KO -> pathways?

ADD REPLYlink written 11 months ago by Asaf8.4k

Both, I think. Going from Taxa -> KO, am I ending up with all possible gene orthologs in that specific organism, or is it the orthologs represented across all the samples in the pool? Maybe that is the difference between "stratified" and "unstratified?" Then, is it just KO -> pathway, or is it KO + EC -> pathway. Otherwise, what is EC for?

ADD REPLYlink written 11 months ago by ariel110

KEGG pathways contain more than enzymes and are built using KOs. I guess EC numbers are there for you if you want to use them but not used to construct pathways. Although EC<->KO is pretty tight and should lead to same pathways.

ADD REPLYlink written 11 months ago by Asaf8.4k
1

After reading more docs, it turns out the picrust2 ONLY uses the EC numbers and gets the pathways from the MetaCyt database. It gives me the KO numbers, but does not match them with KEGG pathways.

ADD REPLYlink written 11 months ago by ariel110

Thanks for the update, good to know

ADD REPLYlink written 11 months ago by Asaf8.4k
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 2.3.0
Traffic: 1813 users visited in the last hour