Forum "Annotathon: généralités, fonctionnement & bugs"

Thread subject: BLAST taxonomic report

[ Return to forums ]
BLAST taxonomic report
TheoRayane
2 Dec 2020 11:37
Non evaluated contribution

 

Bonjour,

 

J'ai un problème avec TaxReport. Voici une partie des résultats renvoyés:

 

Lineage report
.LUCA
. Bacteria
. .Proteobacteria
. . Gammaproteobacteria
. . .Gammaproteobacteria bacterium................ 1140 0.0    12 hits Bacteria:Proteobacteria:Gammaproteobacteria: SusC/RagA family TonB-linked outer membrane protein...            Gammaproteob... Gammaproteob...
. . .Alteromonas sp............................... 605  0.0    1 hit   Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Alteromonadaceae: SusC/RagA family TonB-linked outer membrane protein [Alteromon... Alteromonas sp. Alteromonas sp.
. . .Lacimicrobium alkaliphilum................... 295  8e-85  2 hits  Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Alteromonadaceae: SusC/RagA family TonB-linked outer membrane protein [Lutibacte... Lutibacter sp.  Lutibacter sp. 
. . .Marinobacter persicus........................ 288  5e-82  4 hits  Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Alteromonadaceae: TonB-dependent receptor [Flavobacterium sp. LC2016-01]            Flavobacteri... Flavobacteri...
. . .unclassified Colwellia....................... 331  7e-98  1 hit   Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Colwelliaceae: SusC/RagA family TonB-linked outer membrane protein...            Chitinophaga... Chitinophaga...
. . .Colwellia sp. MB02u-6........................ 327  4e-96  2 hits  Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Colwelliaceae: TonB-dependent receptor [Chryseobacterium soldanellicola]         Chryseobacte... Chryseobacte...
. . .Psychromonas sp. B3M02....................... 320  4e-94  2 hits  Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Psychromonadaceae: SusC/RagA family TonB-linked outer membrane protein...            Sediminibact... Sediminibact...
. . .Idiomarinaceae bacterium..................... 287  8e-82  1 hit   Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Idiomarinaceae: TonB-dependent receptor [Capnocytophaga cynodegmi]                Capnocytopha... Capnocytopha...
. . .Shewanella sp. 10N.286.52.C2................. 287  1e-81  2 hits  Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Shewanellaceae: TPA: SusC/RagA family TonB-linked outer membrane protein...       Chitinophaga... Chitinophaga...
. . .Pseudoalteromonas fuliginea.................. 275  3e-77  2 hits  Bacteria:Proteobacteria:Gammaproteobacteria:Alteromonadales:Pseudoalteromonadaceae: TonB-dependent receptor [Chitinophaga jiangningensis]             Chitinophaga... Chitinophaga...
 

On remarque rapidement dès les premières lignes qu'il y a un problème (ou alors je n'ai pas compris blast). Par exemple, nous avons un lacimicrobium alkaliphilum placé dans les proteobactéries alors qu'en réalité, la protéine renvoyée est un Lutibacter sp. (clade : bacteroidetes). Ce problème est récurrent sur de très nombreuses lignes (presque toutes). Ce n'est pas un problème de navigateur (j'ai essayé sur plusieurs navigateurs sur différents ordinateurs). Egalement, ce n'est pas un problème de report car j'ai essayé avec une autre protéine et les résultats étaient bons.

 

J'avais pensé à regardé le lineage taxonomy de Blastp pour contourner ce problème. Mais voici une partie des résultats renvoyés :

 

 

 

Je n'arrive pas à voir rapidement à quelle classe et à quel ordre appartienne chaque organisme. Il faut que j'utilise le taxonomy browser du NCBI pour le savoir pour chaque résultat...

 

Pourriez-vous me dire si j'ai manqué un élément (par manque de connaissances dans le domaine tel qu'un problème de ma part dans la lecture des résultats) ou si c'est un bug.

 

Je vous mets ci-dessous l'ORF étudié :

 MKNSLLKKPLLILMLLIGSISYGQSVKGVVSDADGSLPGVSITIQGTTTG
VQTDFDGNYTISAKQGDVLVYNYLGYKTEQRTVGSEAVINVTMTQDSTEL
DEIVVIGYGSTTVKDATGSVTSVKSEDFNKGVISSPEQLIQGKTAGVQIS
STSGAPGAGVNIRIRGSNSVRSNNNPLFVVDGIPLSGGATTAGANVIGLG
GNPAKNPLSFLNPNDIESISILKDASATAIYGSRGANGVVIVTTKQGRSK
KGEFNYQANVSSSKARQTYNLLNASQFISERQRLGLPVDAVSNVGSSTDW
QDVVFRTAQSHNQNLSYSKSHDKGSYIATLGFGNQTGIVENSNLERLTAR
LNLNQRFLDDKLTLSFQGSVSRVDDNLAPTGGTAGYKGDVLGAAFSANPT
WPNDPDYDGGTLVNPANLLAYTKGTTNTNRYLANVSLEYDLTSNLKAKVS
VGLDQSDSDNFSVSSSQALGLGDGVAGNGRGAYNVLEVKNELLEATISYN
KEYKNSKLDVLGGFSYQSFQNSGMYSGAWGFSTTDLDQMGSSLKNIVTGI
AGEISGQYQQFAYAPNFTGVAVRGWNGLDFTTTNISPSTSRASVRSLAVD
SYDNT

 

Bien cordialement,