This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to annotate draft genome

Hi all,

I have several bacterial draft genomes assembled by spades. After checking with checkM, I wanna annotate these bacteria which are with 'contaminatin' < 5%, in order to find out their species information. I aligned the draft genome to nt database, but I found one particular genome could be assigned to different bacteria. I am afraid that maybe this approach is not rational. Is there any approach competent in quick bacteria annotation. Thanks in advance.

contigs bacteria draft genome
• 2,541 views • link
updated by andres.firrincieli andres.firrincieli • written by zhangdengwei zhangdengwei

Can you clarify what

which are with 'contaminatin' < 5%,

means?

Are these mixed samples (e.g. metagenomes)? Are you expecting contamination? If you have single, clean draft genome assemblies, prokka is a tool of choice for annotation.

• 0 views • link
written by Joe Joe

Yes, my draft genomes were derived from metagenomes, and I try to split them into single bacteria. Afterward, I used checkM to determine whether it is clean. Each separated bacterium has hundreds of contigs. So now I want to know what the separated bacteria is.

• 0 views • link
written by zhangdengwei zhangdengwei

2 answers

Assuming that you have a computer with at least 128 Gb RAM (or a combo of RAM+swap > 128 Gb), the most consistent way of doing this is by using a GTDB toolkit. Below is an example of the final output for each bin that is more than 10% complete (it is truncated on the right side so as to not run too far out). As you can see, most metagenomic bins are classified down to the family or genus level, with couple of them having a species designation.

group_06        d__Archaea;p__Crenarchaeota;c__Thermoprotei;o__Desulfurococcales;f__Desulfurococcaceae;g__Thermosphaera;s__Thermosphaera aggregans
group_11        d__Archaea;p__Nanoarchaeota;c__Nanoarchaeia;o__Nanoarchaeales;f__Nanopusillaceae;g__Nanopusillus;s__
group_16        d__Archaea;p__Crenarchaeota;c__Thermoprotei;o__Desulfurococcales;f__Acidilobaceae;g__;s__
group_04        d__Archaea;p__Crenarchaeota;c__Thermoprotei;o__Thermoproteales;f__Thermoproteaceae;g__Pyrobaculum;s__
group_18        d__Archaea;p__Crenarchaeota;c__Thermoprotei;o__Thermofilales;f__Thermofilaceae;g__;s__ 
group_15        d__Bacteria;p__Deinococcota;c__Deinococci;o__Deinococcales;f__Thermaceae;g__Thermus;s__Thermus aquaticus
group_03        d__Bacteria;p__Aquificota;c__Aquificae;o__Aquificales;f__Aquificaceae;g__Thermocrinis;s__
group_17        d__Bacteria;p__Desulfobacterota;c__Thermodesulfobacteria;o__Thermodesulfobacteriales;f__Thermodesulfobacteriaceae;g__;s__
• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Many thanks. That's what I want.

• 0 views • link
written by zhangdengwei zhangdengwei

As Mensur said GTDB can be used to determine taxonomic affiliation of your bins. In case you do not have 128 Gb of RAM, GTDB has been implemented in KBase

• 0 views • link
written by andres.firrincieli andres.firrincieli

Log in to answer this question.

More posts like this