This is a test version of Biostars. For the public version, visit https://www.biostars.org.
taxonomic assignment of metagenomic bins and detect 16S rRNA gene in metagenomic bins

Hi, community,

I have Hiseq metagenomic data. I have used MEGAHIT assembly the reads, then I binned the contigs which results to 500 bins. I have used CheckM to assess the genome bins. All bins with <10% contamination and >70% completeness were saved and BLASTN to SILVA 132 SSU database. Most bins can detect 16S genes. However, some bins got hit to different taxonomy, such as one contig got hit with 'FUBE01000005.37491.40675 Bacteria;Proteobacteria;Gammaproteobacteria;Enterobacteriales;Enterobacteriaceae;Escherichia-Shigella;Shigella sonnei', one contig in the bin got hit to '> CP016073.2513631.2516362 Bacteria;Firmicutes;Bacilli;Bacillales;Stapyhlococcaceae;Staphylococcus;Staphylococcus pseudintermedius'. May I know what is the reason? Anything wrong with me?

next-gen assembly
• 2,827 views • link
updated by SilentGene • SilentGene • written by luyang1005 • luyang1005

1 answer

Hi there, apart from using blastn against nr databases to assign taxonomy of the metagenomic bins. I recommend you to try out the following 2 software:

  1. phylophlan. This software could help you integrate metagenomic bins in the 'tree of life' which has already involved over 3000 microbial genomes according to the universal conserved proteins in your bins. And you can probably know the taxonomy of the new bins according to the phylogenetic relationship shown by the supertree.

  2. JSpeciesWS. This is a really handy online tool that can search genomes against their entire genomes reference database GenomesDB according to the Tetra Correlation. The database in JSpeciesWS is really large, which allow you to find more potential of the taxonomy.

Cheers, H

• 1 views • link
written by SilentGene • SilentGene

Log in to answer this question.

More posts like this
  • Add taxonomy information of fileA to species list in fileB
    written by aenna_p • aenna_p
    Hello everyone, I have the following two .csv files (comma-separated): **A working_file:** genus,species,column3,column4,column5,column6,column7 Staphylococcus,aureus,40000,3.0,7.0,6.0,3.0 Neisseria,gonorrhoea,2300,40.0,1.0,3.0,4.4 Vibrio,cholerae,2961,0,47.7,0.0,3.1,0.8 Pseudomonas,aeruginosa,64404,0,66.6,0.0,2.8,8.0 ... **A taxonomy_file** domain,phylum,class,order,family,genus,species Bacteria,Firmicutes,Bacilli,Bacillales,Staphylococcaceae,Staphylococcus,aureus Bacteria,Firmicutes,Bacilli,Bacillales,Staphylococcaceae,Staphylococcus,capitis Bacteria,Firmicutes,Bacilli,Bacillales,Staphylococcaceae,Staphylococcus,saprophyticus Bacteria,Proteobacteria,Gammaproteobacteria,Pseudomonadales,Pseudomonadaceae,Pseudomonas,aeruginosa …
  • Taxonomy Of Blast Hits
    written by Darked89 Darked89
    Lets have 200k genomic contigs with some (unknown) bacterial contamination. I blasted (blastn vs nr) all of them, got tabulated output and passed the uniq …