Main functions and taxonomic distribution of virulence genes in Brucella melitensis 16 M

Many virulence genes have been detected in attenuated mutants of Brucella melitensis 16 M; nevertheless, a complete report of these genes, including the main Cluster of Orthologous Groups (COG) represented as well as the taxonomical distribution among all complete bacterial and archaeal genomes, has not been analyzed. In this work a total of 160 virulence genes that have been reported in attenuated mutants in B. melitensis were included and analyzed. Additionally, we obtained 250 B. melitensis randomly selected genes as a reference group for the taxonomical comparisons. The COGs and the taxonomical distribution profile for 789 nonredundant bacterial and archaeal genomes were obtained and compared with the whole-genome COG distribution and with the 250 randomly selected genes, respectively. The main COGs associated with virulence genes corresponded to the following: intracellular trafficking, secretion and vesicular transport (U); cell motility (N); nucleotide transport and metabolism (F); transcription (K); and cell wall/membrane/envelope biogenesis (M). In addition, we found that virulence genes presented a higher proportion of orthologs in the Euryarchaeota and Proteobacteria phyla, with a significant decrease in Chlamydiae, Bacteroidetes, Tenericutes, Firmicutes and Thermotogae. In conclusion, we found that genes related to specific functions are more relevant to B. melitensis virulence, with the COG U the most significant. Additionally, the taxonomical distribution of virulence genes highlights the importance of these genes in the related Proteobacteria, being less relevant in distant groups of organisms with the exception of Euryarchaeota.