This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Qiime file for Krona pie chart?

Hi,

I have analysis with Qiime of 16s reads, and I want to make a pie chart with Krona but do not know which file from Qiime to use for Krona?

qiime krona

Are you still trying to get it working? I know I had to make it work and figureD out how! Anyways without having to write / debug bash script! There is an R package that does Krona charts! Let me know if you wan tot know more! I'll report on my method if needed! :)

Dear Jeremieauger,

I also working using Krona but with MEGAN output. I face some problem because my command didn't work in Krona. I wonder what the R package that does Krona chart that you told?

Thank you,

Regards,

Lulu

3 answers

Hi, hope this answer is not too late. :)

convert biom to tsv format

biom convert -i single_sample.biom -o single_sample.tsv --to-tsv --table-type "OTU table" --header-key taxonomy

remove the first two rows with comment

# take a look at single_sample.tsv
$ head single_sample.tsv 
# Constructed from biom file
#OTU ID H2O taxonomy
346085  140.0   Bacteria; Proteobacteria; Alphaproteobacteria; Caulobacterales; Caulobacteraceae; Brevundimonas; Brevundimonas_bullata
10298   2.0 Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacteriales; Enterobacteriaceae; Photorhabdus; Photorhabdus_temperata
122823  3.0 Bacteria; Proteobacteria; Betaproteobacteria; Burkholderiales; Oxalobacteraceae; Massilia; EF516371_s
130468  2.0 Bacteria; Proteobacteria; Alphaproteobacteria; Sphingomonadales; Sphingomonadaceae; Sphingopyxis; Sphingopyxis_witflariensis
139977  38.0    Bacteria; Proteobacteria; Alphaproteobacteria; Sphingomonadales; Erythrobacteraceae; Erythrobacter; Erythrobacter_flavus
121751  2.0 Bacteria; Firmicutes; Clostridia; Clostridiales; Lachnospiraceae; Catonella; JX096343_s
96934   10.0    Bacteria; Proteobacteria; Gammaproteobacteria; Pseudomonadales; Pseudomonadaceae; Pseudomonas; Pseudomonas_guguanensis
95181   4.0 Bacteria; Proteobacteria; Gammaproteobacteria; Pseudomonadales; Moraxellaceae; Acinetobacter; Acinetobacter_radioresistens

# remove the comments
$ egrep -v "^#" single_sample.tsv > sample_no_comment.tsv

remove the first column

$ cut sample_no_comment.tsv -f 2- > sample.tsv
$ head -n2 sample.tsv 
140.0   Bacteria; Proteobacteria; Alphaproteobacteria; Caulobacterales; Caulobacteraceae; Brevundimonas; Brevundimonas_bullata
2.0 Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacteriales; Enterobacteriaceae; Photorhabdus; Photorhabdus_temperata

replace the "; " with "\t"

$ sed -i 's/;\s*/\t/g' sample.tsv
$ head -n2 sample.tsv 
140.0   Bacteria    Proteobacteria  Alphaproteobacteria Caulobacterales Caulobacteraceae    Brevundimonas   Brevundimonas_bullata
2.0 Bacteria    Proteobacteria  Gammaproteobacteria Enterobacteriales      Enterobacteriaceae   Photorhabdus    Photorhabdus_temperata

now we can make a pie chart with ktImportText(One of the Krona tools)

$ ktImportText -n enjoy_your_life -o sample.krona.html sample.tsv

open sample.krona.html with your web browser, you'll see this

Krona_snapshot_enjoy_your_life

all in one command(optional)

$ awk -F '\t|;' 'BEGIN{OFS="\t"} FNR > 2 {$1=""; print $0}' single_sample.tsv | cut -f 2- > sample.tsv

krona is awesome! I used it with kraken output glad to know qiime output can be used as well

Yeah, its really awesome and I'm glad my solution is useful to you:)

Hello badribio, Could you please explain how to create kraken output,what is the command to be used??

What have you done so far, and which step are you stuck?

I have finished qiime analysis,but I want to know what output file generated by qiime to provide as an input in kraken.

otu_table_L6.txt file generated by the script summarize_taxa.py will be the best file to be used for generating Krona pie chart for obvious reason (having information for all the taxonomic ranks right from phylum to genus).

Hi Vijay, I have installed Krona and Kraken what is the input file for kraken and krona??

You just have to import your qiime artefacts (table.qza, rooted-tree.qza, taxonomy.qza and sample-metadata.tsv) in the form of a Phyloseq. Then there is a neat package that generates Krona charts:

https://rdrr.io/github/cpauvert/psadd/man/plot_krona.html

I can further develop, but this post seems outdated!

Cheers, -JA

Log in to answer this question.