Bonjour,
tous les choix de groupes d'étude et externe que vous mentionnez ci-dessus sont tous cohérents, car vous prenez bien un groupe de même rang taxonomique que le groupe d'étude pour faire votre groupe externe.
- Classe des gamma- vs classe de alpha-, au sein du phylum des proteo-
- Classe des flavo- vs cytophagia au sein du phylum des bacteroidetes
- 2 classes différentes au sein du phylum des Cyanobacytéries
Et qu'il fallait plutôt prendre 2 groupes étude et externe avec une seule difference( le dernier ordre). Est ce exact?
=> Le groupe externe se choisit en fonction du groupe d'étude, il doit être de même rang taxonomique.
=> Si vous n'avez pas suffisamment de séq au même rang pour faire un groupe externe, et uniquement dans ce cas là, alors vous pouvez remonter d'un cran dans la taxo de réf. pour faire un groupe externe.
Malgré tout ces choix nous avions toujours le même problème au niveau des arbres: un mélange entre groupe externe et groupe étude.
=> Je pense que l'orf n'était pas non plus intégré dans les groupes d'études choisis ?
On a tout de même essayer avec deux groupes ayant 2 ordre différent : LUCA; Bacteria;Spirochaetes:Spirochaetia et LUCA:Bacteria:Cyanobacteria:Cyanophyceae
=> Non, là vous considérer les groupes d'étude et externe au niveau du phylum bactérien : Spirochaetes vs Cyanobacteria.
Vous nous aviez demandé de le refaire avec cette fois ci 5000 séquences
=> Oui, car l'e-value la plus faible, e-51, et donc la séq la plus proche de votre séq étudiée, est dans le groupe des spirochaetes, qui pouvait donc être le groupe d'appartenance. Cependant, vous n'aviez que 17 séq avec le blast avec 1000 seq cibles. L'idée était donc de voir si vous pouviez obtenir plus de seq dans ce groupe taxonomique en augmentant le nombre de cibles demandées.
Ce qui semble être le cas, puisque dans le dernier arbre que vous me présentez, le groupe d'étude est le groupe des spirochaetes, il y a plus que les 17 seq de départ, et votre ORF émerge bien dans ce groupe : il s'agit donc bien du bon groupe d'appartenance.
Si votre GE est le phylum des spirochaetes, alors le groupe externe, est un autre phylum bactérien, cyanobactérie par exemple, oui : le résultat présenté est juste.
Par contre, avec le résultat de blast avec 5000 séq cibles, ne pouviez-vous pas être plus précis dans le choix du groupe d'étude et groupe externe ? (Le but est de trouver le groupe d'appartenance le plus précis possible).
Avez-vous peut-être plusieurs classes, ou plusieurs ordres, ou plusieurs familles différent(e)s dans le phylum des spirochaetes ? Ce qui vous permettraient de choisir un GE plus précis ?
Là, pour en juger, il me faut la "Taxonomy Summury Table" du blast vs 5000 seq.
Celle que vous présentez ici, est celle du blast vs 1000 seq.
Le problème maintenant c'est l'alignement au format ClustalW ( disponible sur notre partie alignement; résultats bruts),
=> Il faut d'abord voir si vous pouvez être plus précis dans votre choix du groupe d'étude et du groupe externe. cf "Taxonomy Summury Table" du blast vs 5000 seq.
Dans ce cas, vous auriez effectivement des séq plus proches incluses dans l'alignement multiple.
Bon travail,
BW