Use of Kraken output for functional analysis
0
0
Entering edit mode
8 months ago
Srinka ▴ 20

Is there anyone available to assist me in utilizing Kraken output for functional analysis? I've been on the lookout for tools capable of processing Kraken report output for functional analysis.

metagenomics kraken shotgun • 1.1k views
ADD COMMENT
0
Entering edit mode

You're going to have to be more specific. "Functional analysis" is a catchall for a lot of different types of possible downstream analyses. Finding recent studies that have similar data is a good place to start in finding appropriate tools.

ADD REPLY
0
Entering edit mode

Hi, by functional analysis I meant the pathways the microbes are involved in.

ADD REPLY
0
Entering edit mode

Hello,

You can convert your classified Kraken report to the BIOM format, and then use the PICRUSt tool to predict the functions associated with KEGG and COG based on the OTU information present in the samples.

Thanks!

ADD REPLY
0
Entering edit mode

Hello, Thanks for your suggestion. I will surely try this.

ADD REPLY
0
Entering edit mode

Can you please suggest any tool to help convert the Kraken report to BIOM format?

ADD REPLY
0
Entering edit mode

Hello,

you can use kraken-biom [kraken_biom.py] tool to get your Biom file from input kreport

cheers,

ADD REPLY
0
Entering edit mode

Hi,

I've gone over the tutorial section of the Picrust GitHub and am confused about the study_seqs.fna file.

  • What should you use for the study_seqs.fna file, and are there any edits required to it?

Context: I ask this because Picrust requires the sequence IDs in the fna file to be linked to the first column of the OTU table. My OTU table is taxa x sample (row x col) and can be linked back to the kraken database I used. I am worried about the IDs in the database. After looking through the database fasta file, I found that multiple sequences linked to the same kraken id (taxa). My first assumption is to combine these sequences into one under the one ID, but I am worried this would negatively impact the alignment process in picrust. However, if I only choose one sequence then is there a way to decide which sequence is the best? If the fna file is supposed to be the ASVs linked to the specific sample, then I have no clue where to start to link those reads to the OTU table.

Many thanks.

ADD REPLY

Login before adding your answer.

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