This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Clusters of orthologous research tools

Hey, I am annotating 10 bacterial genomes belonging to the same species and I am looking for a way to analyze the clusters of orthologous genes. Can you tell me which path to follow to identify them and perhaps gather them in a table? Is there a tool that allows you to compare the protein sequences of the various genomes and identify the shared COGs?

Thanks for any help,

A.

genome annotation
• 1,347 views • link
updated by dthorbur dthorbur • written by anna • anna

2 answers

Hi Anna, there are many options for this kind of gene-based pangenome analysis.

I will start with our own https://github.com/eead-csic-compbio/get_homologues, which you can use following the bacteria tutorial and manual . Compared to other options out there I guess it has two advantages:

  • It supports different clustering algorithms (BDBH, OMCL, COGS) and you can get consensus core clusters.
  • You can pipe resulting clusters to https://github.com/vinuesa/get_phylomarkers to compute robust phyologenies out of sequence clusters.

GET_HOMOLOGUES can be installed with conda and can use a HPC cluster to parallelize tasks. A Docker container shipping with both GET_HOMOLOGUES and GET_PHYLOMARKERS is available at https://hub.docker.com/r/csicunam/get_homologues .

Another obvious software choice would be https://sanger-pathogens.github.io/Roary for instance.

• 0 views • link
written by b.contreras.moreira • b.contreras.moreira

I'm a big fan of OrthoFinder. Whilst it's designed for interspecific comparisons, I have used it for different assemblies of the same species before with success.

It even provides gene trees for each orthogroup, and an overall species tree based on all single-copy orthologs trees.

• 0 views • link
written by dthorbur dthorbur

Log in to answer this question.

More posts like this