- Research article
- Open Access
Genomes-based phylogeny of the genus Xanthomonas
BMC Microbiologyvolume 12, Article number: 43 (2012)
The genus Xanthomonas comprises several plant pathogenic bacteria affecting a wide range of hosts. Despite the economic, industrial and biological importance of Xanthomonas, the classification and phylogenetic relationships within the genus are still under active debate. Some of the relationships between pathovars and species have not been thoroughly clarified, with old pathovars becoming new species. A change in the genus name has been recently suggested for Xanthomonas albilineans, an early branching species currently located in this genus, but a thorough phylogenomic reconstruction would aid in solving these and other discrepancies in this genus.
Here we report the results of the genome-wide analysis of DNA sequences from 989 orthologous groups from 17 Xanthomonas spp. genomes available to date, representing all major lineages within the genus. The phylogenetic and computational analyses used in this study have been automated in a Perl package designated Unus, which provides a framework for phylogenomic analyses which can be applied to other datasets at the genomic level. Unus can also be easily incorporated into other phylogenomic pipelines.
Our phylogeny agrees with previous phylogenetic topologies on the genus, but revealed that the genomes of Xanthomonas citri and Xanthomonas fuscans belong to the same species, and that of Xanthomonas albilineans is basal to the joint clade of Xanthomonas and Xylella fastidiosa. Genome reduction was identified in the species Xanthomonas vasicola in addition to the previously identified reduction in Xanthomonas albilineans. Lateral gene transfer was also observed in two gene clusters.
Xanthomonas is a genus in the gamma division of Proteobacteria primarily constituted by pathogens to plants of considerable economic importance. These pathogens affect a wide variety of crops, including Citrus spp. (lime, orange, lemon and pomelo, among others), Oryza spp. (rice), crucifers (cabbage, broccoli, cauliflower, radish and Arabidopsis thaliana) and Manihot esculenta (cassava), with individual members showing a high degree of host specificity . Xanthomonas is among the few bacterial genera in which large DNA-DNA hybridization, RFLP and REP-PCR datasets are available [2–6] and have been employed for the taxonomical resolution of the group . In addition, the availability of more than ten genomes within the genus [8, 9] has allowed recent studies of comparative genomics and genome evolution [10, 11].
The genus Xanthomonas has been subject to numerous taxonomical and phylogenetic studies, starting with the description of Bacterium vesicatorium as the causal agent of bacterial spot on pepper and tomato  and its reclassification as Xanthomonas campestris [13, 14]. Xanthomonas was first described as a monotypic genus, and later divided in two groups, A and B [15, 16]. A subsequent study  classified 183 reported strains into 20 different species mainly based on DNA-DNA hybridization data. Since then, a general classification has been established based on polyphasic analysis [6, 17], while other analyses helped to clarify the classification in specific clades, mainly using Multi Locus Sequence Analysis (MLSA) and Amplified Fragment Length Polymorphism (AFLP) [18, 19]. This allowed the development of several typing and characterization tools (e.g., [11, 18–24]), which have revealed the diversity and complexity of the genus [23, 24], while showing the limitations of single locus analyses . However, during the last decade the taxonomy of this genus has still been subject to considerable debate. Genus-wide reclassifications have been proposed [26, 27], and frequent sub-specific reclassifications and proposals for new species have been published [19–21, 28–30].
A remarkable example of these conflicts is the classification of X. fuscans aurantifolii [26, 27], also known as X. axonopodis pv. "aurantifolii" [2, 6, 18, 31]. This taxon was originally identified as part of the DNA hybridization homology group "X. axonopodis" , but after its differentiation from other xanthomonads by DNA sequence-based molecular techniques, production of water-soluble brown pigment and host range, it was designated as X. fuscans . However, when these traits/methods were examined, none of them could individually differentiate X. fuscans from other pathovars within X. axonopodis [18, 31]. DNA-DNA reassociation assays, in turn, have differentiated X. fuscans from X. axonopodis, X. campestris and X. citri [2, 26, 27]. Additional host-range evidence has also been used to support the designation X. fuscans, separated from X. axonopodis and X. citri. Phaseolus vulgaris and Citrus spp. are infected by X. fuscans pvs. fuscans and aurantifolii, respectively, but are not infected by either X. axonopodis or X. campestris. Citrus spp., on the other hand, is also infected by X. citri . However, host range is usually a criterion to separate pathovars and not species. This example underscores the importance of a solid taxonomic classification with a phylogenetic basis.
Molecular phylogenetics has played an important role in the classification of the genus. Single locus analyses, including the use of 16S-23S rDNA spacers, the 16S rRNA gene and the DNA gyrase gyrB [32–35], generally agree with standing nomenclature but with low resolution below the species level. MLSA including sequences of protein-coding genes dnaK, fyuA and rpoD , has significantly extended previous results. In general, MLSA results suggest that X. citri and X. fuscans are closely related species and should be considered as a single species based on their 98.34% similarity in the proteins encoded by dnaK, fyuA, gyrB and rpoD . Recently, a phylogenomic approach was applied to resolve the phylogenetic relationships within the genus , although this work did not explore the phylogenetic distances between strains, and did not include sequences from X. axonopodis species. The general structure of the genus agreed with the standing nomenclature.
The use of genomic sequences as the basis for species delimitation has been explored as a new standard in bacteria in replacement of DNA-DNA hybridization [36, 37], particularly based on metrics such as the ANI (Average Nucleotide Identity) . The correspondence between DNA-DNA hybridization and sequence similarity has been exploited in Xanthomonas for the establishment of clades and species , but full genomic sequences have not been used so far for the resolution of the "X. axonopodis" clade (this is, including close relatives such as X. fuscans and X. euvesicatoria). Phylogenomic methods extend the analysis of primary sequence data from one or few loci (usually no more than twenty) to hundreds or thousands of loci at the same time, alleviating the problem of incongruence between characters [39, 40]. Here, we present a phylogeny of the genus based on seventeen complete and draft genomes, including five genomes from the "X. axonopodis" clade. We identified the orthologous genes and performed the phylogenetic inferences using a new library called Unus, which is briefly described here.
The automated selection of orthologous genes is consistent with manual selection
In order to compare a typical literature-based selection of genes for phylogenetic reconstruction in bacteria with the Unus automated method, using 989 genes in the genomes listed in Table 1, we evaluated the presence of the housekeeping genes used by AMPHORA . We found that several of these genes were absent in the draft genomes Xfa1, Xfa0 and Xvm0. In addition, in-paralogs (i.e., duplicated genes) were detected in the genome of XooK for several ribosomal proteins (large subunit; rplA, rplC, rplD, rplE, rplF, rplN) and were therefore discarded. This is possibly due to errors in the genome sequence, given that these genes are usually present as a single copy. Importantly, the absence of rpl genes in the XooK genome suggests that ribosomal proteins (from both the small and the large subunits) were located at mis-assembled regions of the genome sequence. Genes employed in the genus-wide analysis and used by AMPHORA include dnaG, nusA, pgk, pyrG, rplM, rplP, rplS, rplT, rpmA, rpoB, rpsB, rpsC, rpsE, rpsI, rpsK, rpsM and rpsS. Also, five out of the seven genes used by Pieretti et al.  (gyrB, recA, dnaK, atpD and glnA) were found in the constructed Orthology Groups (OG), while other two (groEL and efp) seemed to be absent in the draft genome of Xfa1. This underscores the importance of a flexible selection criterion of orthologous genes in a determined group of taxa, especially with unfinished genomes. A previous MLSA conducted by Young and collaborators  employed four protein-coding genes included in the previous lists plus the tonB-dependent receptor fyuA, also present in our selection. Another MLSA recently performed by Bui Thi Ngoc et al.  used the genes atpD, dnaK, efP and gyrB, all of which were present in our dataset. These data suggest that the automated selection using Bit Score Ratio (BSR) is in agreement with the classical selection of genes for phylogenetic studies. Therefore, some of the genes selected in this study can be used for future phylogenetic reconstructions.
The COG classification for the employed genes (Additional file 1) was compared among sets of genes obtained from automated selections at different taxonomical levels within the genus (Figure 1). COG categories related to central metabolism and ribosomal proteins presented a tendency to increase in representation (relative to other COG categories), as genomes from a wider taxonomical range were included (blue bars in Figure 1). Together, these categories covered 27% of the COG-classified genes and included genes that are frequently used for phylogenetic reconstruction. On the other hand, a reduction in the relative representation when including a wider taxonomical range of genomes was observed for categories related to peripheral metabolism and poorly characterized proteins (red bars in Figure 1). These categories covered 36.9% of the COG-classified genes and included clade-specific genes (without detectable orthologs in distant relatives) as well as genes absent in X. albilineans, which presents a notable genome size reduction . Pieretti and collaborators identified 131 ancestral genes potentially lost by pseudogenization or short deletions in X. albilineans and 480 potentially lost by both X. albilineans and Xylella fastidiosa . Most of the COG-classified genes putatively lost in X. albilineans or both X. albilineans and Xylella fastidiosa (56.2% and 56%, respectively) can be classified within these COG categories. The same tendency to increase in relative representation when increasing the number of taxa was displayed by genes without an assigned COG category (data not shown). The only category significantly impacted by discarding the in-paralogs was category L (replication, recombination and repair). This category covers 8.2% of the COG-classified genes, and 83.2% of those discarded by paralogy, suggesting frequent duplications of genes implicated in these processes. Putative transposases and inactive derivatives represent 76% of the discarded genes.
Phylogeny of the genus Xanthomonas
Our phylogenetic analysis was based on 989 OG (1,084,777 bp, Additional file 2), which included all markers used in previous Xanthomonas phylogenetic analyses. Both, the Maximum Likelihood tree and the Bayesian consensus tree reconstructed the same well-supported topology, with bootstrap supports of 100% for all the nodes (out of 1,001 replicates). The same relationships were also obtained with Maximum Parsimony (bootstrap support of 100% with 1,000 replicates).
A total of four clades were obtained in the phylogenomic reconstruction. The first clade includes X. oryzae, the second comprises X. vasicola, the third one groups together X. fuscans, X. euvesicatoria and X. axonopodis, and the fourth clade contains X. campestris (Figure 2a). These results agree with previous phylogenies of the genus [11, 17, 35, 42]. In order to further advance on the knowledge of the ancestral relationships of the genus Xanthomonas, and in particular the species Xylella fastidiosa, we performed a new analysis including three additional genomes in the Xanthomonadaceae family: Xylella fastidiosa str. 9a5c (GenBank entry AE003849.1), also a plant pathogen, but strictly transmitted by insect vectors; Pseudoxanthomonas suwonensis str. 11-1 (GenBank entry CP002446.1), a bacterium isolated from environmental samples but more commonly found in contaminated ones; and Stenotrophomonas maltophilia str. R551-3 (GenBank entry NC_011071.1), a common soil colonizer which has also been reported as a human opportunistic pathogen. These species are hereafter termed Xyf9, Pxs1 and StmR, respectively. This new analysis was based on a collection of 228 genes automatically compiled by the Unus library using Bit Score Ration (BSR). The resulting phylogeny revealed that the genus Xanthomonas is not monophyletic, with Xylella fastidiosa as its sister clade. X. albilineans should be placed in an independent genus in order for the taxonomy to match the phylogeny of the group (Figure 2b), as previously noted . This result differs from that presented by Pieretti and collaborators, based on seven housekeeping genes , where X. albilineans and X. fastidiosa form a single clade ancestral to all other Xanthomonas.
Genome evolution: gains and losses
The high number of pseudogenes and lost regions in X. albilineans suggests a reductive genome evolution in this species . This information, together with the position of the taxon in previous phylogenies [11, 42] and the reduced size of the close relative Xylella fastidiosa , could indicate either a reduced genome as the ancestral condition in the Xanthomonas genus or independent genome reductions in Xylella fastidiosa and X. albilineans. Pieretti and collaborators provide strong evidence supporting the latter hypothesis . However, the enrichment of phage-related regions in the Xylella genomes, as well as the presence of multiple Insertion Sequences (IS) in Xanthomonas reveal very active mobile elements in the Xanthomonadales order . To determine whether this reductive tendency extends to other genomes of the genus, we employed GenoPlast  for the detection of ancestral genomic gains and losses. The results (Figure 3 and Additional file 3) revealed that all the tip nodes in the X. oryzae species present net genomic losses compensated by genomic gains in ancestors of the species (i.e., internal nodes 20 and 24, as labeled in Additional file 3). Interestingly, the three genomes of the species X. vasicola presented large genomic gains (between 12.78% and 15.19% of the regions) after genomic losses exhibited by the most recent ancestral node of the species (11.47% of the regions). This level of genomic losses is almost twice as large as that exhibited by X. albilineans (5.92%), suggesting that the X. vasicola genomes are very dynamic, while maintaining a genome size comparable to other species in the genus.
Gene clusters and detection of putative gene transfer by orthology groups
In order to identify the distribution of OGs among taxa within Xanthomonas, a second set was constructed using OrthoMCL . Figure 4 depicts the general distribution, clustering by patterns of presence/absence among genomes, regardless of their relatedness. In general, the patterns presented by most of the OGs are monophyletic, as expected (blue columns in Figure 4). However, a few paraphyletic patterns were unexpectedly enriched. Further inspection revealed that most of the OGs in two of the most enriched paraphyletic patterns are clustered in the genomes and preserve synteny. We explored these patterns, and found two clusters of contiguous genes with paraphyletic distributions, suggesting horizontal transference of genetic material.
The first cluster (Figure 5a) is present in Xci3, Xeu8, Xcc8 and XccB, but absent in other genomes of X. campestris, in X. axonopodis and in X. fuscans. Similar genes were also found in Pseudomonas aeruginosa, Salmonella enterica and other species of the genera Pseudomonas, Salmonella and Acidovorax (Additional file 4). This cluster is mainly composed of putative secreted and membrane proteins, with few characterized orthologs. In Xanthomonas, only three of those genes have been characterized. The first two code for VirD4 and VirB4, which are proteins implicated in protein secretion by the Type IV secretion system in several bacteria, including Helicobacter, Agrobacterium and Bartonella [59, 60]. The third codes for RadC, a protein involved in DNA repair. The gene at the locus XCV2366_1 from Xeu8 presents homology with the oxidoreductase DbsA, an important protein for oxidative folding of disulphide-bonded proteins in Gram-negative bacteria . Only nine out of the nineteen genes in this cluster present a G+C content at least one standard deviation distant from the average for the coding regions within the Xeu8 genome (64.66 ± 3.91%). The values of Codon Adaptation Index (CAI) for the seventeen genes in the cluster were similar to the values obtained for other regions of the genome. The distribution of this cluster along the genus suggests flow of genetic material between different pathovars of Xanthomonas. However, G+C content and CAI analyses failed to relate this cluster to LGT. Furthermore, LGT regions predicted by AlienHunter  do not cover more than one gene in this region in any of the analysed genomes (data not shown). Interestingly, in all the genomes, predicted LGT regions surround the cluster at distances from one to eight Kbp.
The second cluster (Figure 5b) is present in XamC and Xfa0 but not in Xfa1, despite the high genome-wide similarity presented between Xfa1 and Xfa0 (Figure 2a). The classification of putative homologs of the genes in this cluster (see methods) revealed that it is mainly composed of sequences similar to proteins in Escherichia coli, Siphoviridae, Stenotrophomonas sp. SKA14, Salmonella enterica and Pseudomonas aeruginosa (Additional file 5). Moreover, members of the Siphoviridae viral family are known to be Pseudomonas and Xanthomonas phages, suggesting the presence of virus-mediated LGT. We cannot attribute the pattern to the mixture of chromosomal and plasmidic DNA in draft genomes (XamC and Xfa0), because none of the sequences presented similarity with genes in Xanthomonas plasmids. Note that the gene at the locus XAUC_17260_1 (Xfa0:1726 in Figure 5b) was originally annotated as yapH, but its product is a large protein of 1231 aa in Xfa0 and 1482 aa in XamC, putatively xenologous with a component of a phage tail (group COG4733 in the COG database). Two genes in the cluster (XamCg00977 and XamCg00978) presented a G+C content more than one standard deviation below the mean of the coding sequences in the XamC genome (i.e., 64.82 ± 3.31%), and a low CAI with respect to the whole predicted coding sequences (0.516 and 0.486, respectively). The other seven genes in the cluster presented average features, which would have precluded their identification as units potentially under LGT.
The results of the genome-based phylogenetic reconstruction suggest that certain changes should be considered in the nomenclature of the Xanthomonas genus. For instance, X. fuscans was recently proposed as a new species , but here we show that it should be considered as a later heterotypic synonym of X. citri, as previously suggested [18, 31]. Other clades in the standing bacterial nomenclature  within the Xanthonomonas genus were consistent with the phylogenetic reconstruction. Nevertheless, we observed a paralogy in the genus Xanthomonas when Xylella fastidiosa was included with X. albilineans outside the Xanthomonas group. Our results suggest that X. albilineans, probably along with other early-branching Xanthomonas, should be considered for a new genus designation. However, the relationships between X. albilineans, Xylella and the other Xanthomonas remain unclear. Another shared feature between Xylella fastidiosa and X. albilineans is the reduced genome. The reductions in these genomes were previously shown to be due to independent events . Here we show evidence suggesting that reductive genome evolution could also affect other clades in the genus such as X. vasicola.
The phylogenetic relationship between X. albilineans, Xylella fastidiosa and the rest of the taxa in the genus Xanthomonas is not clear. The genome of X. albilineans is part of the "early-branching species" , a group of species including X. albilineans and X. sacchari previously found to be basal in the phylogeny of the genus [7, 35]. The species is also a member of the "hyacinthii" group, a group of species with major differences in the 16S-23S rDNA Intergenic Spacer (ITS) with respect to the other members of the genus . Pieretti and collaborators  suggested that Xylella and X. albilineans form a monophyletic clade, which is basal to the rest of Xanthomonas. This is based on a Maximum Likelihood analysis with seven housekeeping genes. Our analyses with over two hundred genes suggest that X. albilineans is basal to Xylella and the rest of taxa in the genus Xanthomonas. Neither of the analyses obtains a good support value for these nodes. The most straightforward explanation for this is that certain regions of the genome support one topology and certain others support the second one. This could be due to a considerable number of LGT in these genomes. Alternatively, it could be due to the large amount of changes accumulated in Xylella fastidiosa, as revealed by the length of the corresponding branch (Figure 2b).
The phylogenetic tree presented in Figure 2a displays identical topology and similar relative branch lengths as inferred by different optimality criteria (Maximum Likelihood, Bayesian Inference, Maximum Parsimony). The tree supports monophyly in the species X. campestris, X. oryzae and X. vasicola. The clade "X. axonopodis" contains the species X. fuscans, X. citri, X. axonopodis and X. euvesicatoria. However, the lower coverage in terms of sequenced genomes of these species makes it difficult to support any further observation beyond the close relatedness within the clade with respect to other species.
Interestingly, the phylogeny displays a close relationship between the species X. fuscans and X. citri. In order to compare their similarity in the same framework of MLSA performed for other species of Xanthomonas (e.g., ), we constructed a matrix containing 989 loci employed for the phylogenetic inference (Table 2). According to the resulting matrix, a similarity threshold of 99% can differentiate bacteria recognized as belonging to the different pathovars (except in X. vasicola, for which pathovars vasculorum and musacearum display a similarity above 99%, possibly due to non-chromosomal sequences). All the species with currently accepted names  have similarities above 97%. This value (in accordance with previous MLSA calibrations ) also differentiate species outside the X. axonopodis clade, but fails to differentiate X. fuscans and X. citri, suggesting that the two pathovars conform a single species as previously suggested [18, 31]. This is also supported by the likelihood distances between these two taxa (Figure 2a, Table 2). Accordingly, we recommended that the species X. fuscans be regarded as a heterotypic synonym of X. citri.
Several robust methods for the identification of orthology, multiple sequence alignments and phylogenetic inferences have recently been developed (reviewed in ). However, a common flexible framework for their joint application in specialized phylogenetic studies and MLSA in general is still required. The BioPerl libraries, including the Bio::Phylo package [65, 66], provide valuable tools for the automation of analyses, but the connections between different steps are often not automated, making them time-consuming. Unus allows the execution of complete workflows in phylogenomics within a single interface, and its current functionalities and limitations underscore the need for a fully structured platform in the field, such as those available for other branches of genomics.
We compared the automatically selected OGs for the phylogenetic assessment with several lists of genes manually compiled. These comparisons indicated that, depending on the genome coverage and annotation of the drafts employed, our analyses broadly agree in the selection of OGs with those utilized previously for phylogenetic inference. Furthermore, the functional distribution of the automatically selected genes exhibits the expected behaviour at different taxonomical levels. Selections on broader taxonomical levels exhibit a larger representation of genes implicated in central-metabolism, while the proportion of clade-specific genes augments in narrower taxonomical levels.
The analysis of the distribution of COG categories shows that central metabolism and ribosomal proteins are favoured when comparing distant genomes, as they are in phylogenetic studies based on one or few loci. Genes in these categories are better suited than genes in other COG categories or unclassified genes because of two characteristics that are important for phylogenetic assessment. Firstly, genes implicated in central-metabolism and ribosomal genes are usually of single-copy. Genes with in-paralogs are normally avoided in phylogenetic inferences given the difficulty in identifying corresponding genes in sets of paralogy , despite some efforts to include them in phylogenetic analyses (e.g., ). Secondly, these genes are often present even in genomes from loosely related organisms. Although phylogenetic reconstructions based on gene content have proven successful (e.g., ), it is hard to achieve high resolution below species and it is not possible with incomplete draft genomes.
Additional genes suitable for phylogenetic analyses were detected through automated identification of orthologs, allowing a higher resolution among closely related taxa. These genes are usually not included in MLSA, although they can add important information about relationships within the group. For closely related bacteria (such as the X. oryzae pv. oryzae strains), the importance of such additional information resides on the low variability among genomes. Therefore, the option to select orthologs without a priori knowledge of the genes that will be included, allows for flexibility in terms of data availability, as well as the obtention of optimized phylogenetic resolution at any taxonomic level under study.
A previous study  suggested a reductive evolution in the genome of X. albilineans, revealed by the small genome (3.77 Mbp) and the high putative pseudogenization. We present evidence supporting the hypothesis that the reductive genome evolution occurs along the genus, and is not restricted to the species X. albilineans. In our analyses, the species X. albilineans effectively revealed large genomic reductions, but even larger reductions were presented by the species X. vasicola, with recent genomic gains only detected on tip nodes, suggesting a reductive evolution tendency followed by the acquisition of genomic regions. The genomic gains on tip nodes can be partly explained by the inclusion of non-chromosomal material in the draft genomes of X. vasicola, although this result was not found in other draft genomes in the study that have non-chromosomal material, such as XamC. An alternative explanation is that genomic gains have arisen by recent genetic exchange with other bacteria, as previously suggested for X. vasicola . However, the large ancestral losses cannot be explained by means of the incompleteness of the genomes, and may reflect an ancestral genomic reduction in the species. The size of the regions involved in such events, and whether they affect restricted functional categories of genes or random regions, is still to be determined.
We identified two clusters of genes with paraphyletic distribution, suggesting lateral gene transfer. One of the clusters, present in X. campestris and the "X. axonopodis" clade, exhibits interesting functional relationships with the Type IV Secretion System (T4SS), while most of the genes are annotated as coding for either putative secreted or membrane proteins. Identification of LGT events based only on intrinsic features such as the G+C content and the CAI would fail to identify both clusters, showcasing the usefulness the phylogenetic distribution of orthologs as a complement for the prediction of putative LGT events.
Currently, phylogenomic methods are finding a privileged place in phylogenetic inference and evolutionary studies, yet common frameworks for the flexible automation of workflows are not widely available. Here we used Unus, a package developed to facilitate the execution of phylogenetic workflows, to explore the phylogenetic structure of the genus Xanthomonas. We recovered a strongly supported phylogeny in accordance with previous results and high resolution in the closely related genomes of X. oryzae. The results also provide evidence for the reconsideration of the X. fuscans species, clarify relationships between X. citri, X. axonopodis and X. euvesicatoria, and show that the genus Xanthomonas is not a monophyletic clade. Our results allowed us to identify several interesting features in the evolution of Xanthomonas, including two large putative lateral gene transfer events, which would have been hard to detect by means of G+C content deviation or Codon Adaptation Index. We also detected evidence of an evolutionary tendency towards a reduction in genome size in at least two clades of the genus.
Seventeen Xanthomonas genomes were used in this study (Table 1). The names employed follow the list of prokaryotic names with standing nomenclature (LPSN) , although several additional names may exist in the scientific literature. Whenever possible, the strains have been tracked to the corresponding StrainInfo entry , in order to ease the resolution of strains deposited in different collections. Gene and gene product predictions were downloaded together with the genomes from NCBI (when available) and JCVI websites, except for the genome of X. axonopodis pv. manihotis str. CIO151 (unpublished), for which coding sequences (CDS) were predicted using Glimmer 3  trained with the X. euvesicatoria str. 85-10 CDS . All the genomes are referred to as stated in the abbreviation column in Table 1.
Generation of Unus, a new library for the execution of phylogenomic workflows
Unus is a Perl library that enables the easy execution of phylogenomic workflows including the detection of groups of orthologous genes, batch alignment of sequences, generation of files in a variety of formats and integration of accessory tests for recombination and models of evolution. The various possible workflows the user can go though in order to obtain a phylogenomic inference of the group of bacteria of interest are depicted in Figure 6. Fourteen Perl modules integrating the Unus package are available for download and code browsing at http://github.com/lmrodriguezr/Unus/. Figure 6 summarizes the different pipelines implemented with Unus and alternative programs that can be used.
Detection of orthologous groups
For the detection of Orthologous Groups (OG), we used the distribution of the Bits Score Ratio (BSR), a BLAST-based metric  essentially as previously described . Briefly, the BSR is defined as the proportion of the Bit Score of the alignment of the query sequence and the subject sequence, and the Bit Score of the alignment of the query sequence with itself (i.e., the maximum Bit Score for a given query). The histogram is usually bimodal (Additional file 6), and Unus detects the valley of the distribution as the threshold to accept a hit for each paired comparison. To avoid spurious results in distributions with shallow valleys or with no evident valley, the threshold for three distributions was set as the average threshold (as calculated for the other paired comparisons). This method accounts for the problems previously observed when considering the best hit only [73, 74], as in widely used methods such as the BLAST Reciprocal Best Match (RBM), also implemented for comparison (see Additional file 7 for the annotated pseudo-code).
Multiple sequence alignments were performed using MUSCLE  on each detected OG. Alignments were discarded when a strong signal of recombination was detected in the Phi test , i.e., p-value ≤ 0.01 under the null model of no recombination. Phylogenetic inference based on whole genomes used Maximum Likelihood (ML) optimality criterion, as implemented in RAxML v7.2.6 [77, 78] with the GTRCAT option, which takes the GTR model of nucleotide substitution, plus an approximation of the Gamma model of rate heterogeneity into account. Branch support was assessed using bootstrap sampling as previously reported . Analyses were performed with each gene in a separate partition to which an independent model of evolution was applied. The resulting ML phylogeny was compared with the consensus topology obtained from Bayesian Inference (BI) [79, 80], with exploration of parameters using the Metropolis-Coupled Monte Carlo Markov Chain (MC3) algorithm with one million generations, as implemented in MrBayes v3.1.2, sampling a tree every 1,000 generations. The log-likelihood scores of sampled points were plotted against generation time to determine when the chain became stationary. All sample points prior to this (300,000 trees) were discarded as burn-in samples. Data remaining after discarding burn-in samples were used to generate a majority rule consensus tree, where percentage of samples recovering any particular clade represented the posterior probability of that clade. Probabilities ≥ 95% were considered indicative of significant support. Branch lengths of the consensus tree were estimated by maximum likelihood . We performed additional phylogenetic reconstructions using Maximum Parsimony (MP) using the PAUP* package v4.0b10 . MP trees were obtained in an equal weighted heuristic search with tree-bisection-reconnection (TBR) branch swapping. The consensus tree was calculated using majority rule. Bootstrap (1,000 replicates, heuristic search TBR branch swapping) was used to assess support for each node. A similarity matrix of all the concatenated sequences was prepared using the DNADIST program of the PHYLIP package  using Kimura distance , in order to compare the distances within the "X. axonopodis" clade with previous MLSA.
Detection of genomic gains and losses
The genomic gains and losses were identified and quantified using GenoPlast  with 10,000 burn-in iterations followed by 100,000 additional iterations, 10 iterations between sampling and two independent runs with identical parameters. Analyses were performed assuming a single phylogenetic tree obtained by ML inference. The input multiple alignment was conducted with progressive Mauve , and post-processed with the tools for developers of Mauve  to first obtain a binary matrix of presence/absence by region, and afterwards a matrix of presence/absence patterns counts. GenoPlast processes this matrix for the calculation of probabilities of ancestral events of genomic gains and losses and implements a model-based method to infer the patterns of genome content evolution by Bayesian inference, assuming a Poisson distribution of genomic gains and losses. The phylogeny inferred here was used as scaffold.
Assignation of COG functional categories
Homology with entries in the Cluster of Orthologous Groups of proteins (COG) database  was determined by BLAST searches  against the COG sequences database. The BLAST search was performed using the default filtering algorithm and a minimum quality of alignments defined by a score of at least 250 bits, an identity of 50% of the aligned region or more, and an aligned region comprising 50% of the query gene or more. BLAST results were parsed and filtered using a custom Perl script with the above criteria. The Perl script also mapped the hits to the corresponding COG category, reporting the category or categories for each query sequence. Each set was analysed 1,000 times randomly sampling 75% of the query sequences to calculate the Standard Deviation (SD; Figure 1). For the characterization of OGs, each comprising one gene per genome, only genes present in the genome of X. euvesicatoria str. 85-10 were used as representative of the OG.
Taxonomical distribution of homologous sequences
BLAST searches against the non-redundant protein database of the NCBI (NR)  were performed in order to identify the homologs of one or more genes in other organisms, with default parameters and Expect value below 10-10. The BLAST result was subsequently parsed with a custom Perl script to extract the organisms, subsequently building a cumulative counts table and mapping these organisms to any fixed taxonomical level using the NCBI's Taxonomy database .
Hayward AC: The host of Xanthomonas. Xanthomonas. Edited by: Swings J-G, Civerolo EL. 1993, London: Chapman & Hall, 52-54.
Egel DS, Graham JH, Stall RE: Genomic relatedness of Xanthomonas campestris strains causing diseases of Citrus. Appl Environ Microbiol. 1991, 57: 2724-2730.
Louws FJ, Fulbright DW, Stephens CT, de Bruijn FJ: Specific genomic fingerprints of phytopathogenic Xanthomonas and Pseudomonas pathovars and strains generated with repetitive sequences and PCR. Appl Environ Microbiol. 1994, 60: 2286-2295.
Rademaker JLW, Hoste B, Louws FJ, et al: Comparison of AFLP and rep-PCR genomic fingerprinting with DNA-DNA homology studies: Xanthomonas as a model system. Int J Syst Evol Microbiol. 2000, 50: 665-677. 10.1099/00207713-50-2-665.
Simões THN, Gonçalves ER, Rosato YB, Mehta A: Differentiation of Xanthomonas species by PCR-RFLP of rpfB and atpD genes. FEMS Microbiol Lett. 2007, 271: 33-39. 10.1111/j.1574-6968.2007.00691.x.
Vauterin L, Hoste B, Kersters K, Swings J: Reclassification of Xanthomonas. Int J Syst Evol Microbiol. 1995, 45: 472-
Parkinson NM, Aritua V, Heeney J, et al: Phylogenetic analysis of Xanthomonas species by comparison of partial gyrase B gene sequences. Int J Syst Evol Microbiol. 2007, 57: 2881-2887. 10.1099/ijs.0.65220-0.
Koebnik R: The Xanthomonas Resource. [http://www.xanthomonas.org/]
Ryan RP, Vorhölter F-J, Potnis N, et al: Pathogenomics of Xanthomonas: understanding bacterium-plant interactions. Nature reviews. Microbiology. 2011, 9: 344-355.
Blom J, Albaum SP, Doppmeier D, et al: EDGAR: a software framework for the comparative analysis of prokaryotic genomes. BMC Bioinforma. 2009, 10: 154-10.1186/1471-2105-10-154.
Moreira LM, Almeida NF, Potnis N, et al: Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp. aurantifolii. BMC Genomics. 2010, 11: 238-10.1186/1471-2164-11-238.
Doidge EM: A tomato canker. Ann Appl Biol. 1921, 7: 407-430. 10.1111/j.1744-7348.1921.tb05528.x.
Dowson WJ: On the systematic position and generic names of the gram negative bacterial plant pathogens. Zentralblatt fur Bakteriologie, Parasitenkunde, Infektionskrankheiten und Hygiene. 1939, 177-193. Abteilung,
Dye DW: Genus IX. Xanthomonas. Dowson (1939). A Proposed Nomenclature and Classification for Plant Pathogenic Bacteria. Edited by: Young JM, Dye DW, Bradbury JF, Panagopoulos GC, Robbs CF. 1978, 153-177. N Z J Agric Res 21;
Stall RE, Beaulieu C, Egel DS, et al: Two genetically diverse groups of strains are included in Xanthomonas campestris pv. vesicatoria. Int J Syst Bacteriol. 1994, 44: 47-53. 10.1099/00207713-44-1-47.
Vauterin L, Swings J, Kersters K, et al: Towards an improved taxonomy of Xanthomonas. Int J Syst Bacteriol. 1990, 40: 312-316. 10.1099/00207713-40-3-312.
Rademaker JLW, Louws FJ, Schultz MH, et al: A comprehensive species to strain taxonomic framework for Xanthomonas. Phytopathology. 2005, 95: 1098-111. 10.1094/PHYTO-95-1098.
Ah-You N, Gagnevin L, Grimont PAD, et al: Polyphasic characterization of xanthomonads pathogenic to members of the Anacardiaceae and their relatedness to species of Xanthomonas. Int J Syst Evol Microbiol. 2009, 59: 306-318. 10.1099/ijs.0.65453-0.
Young JM, Wilkie JP, Park D-S, Watson DRW: New Zealand strains of plant pathogenic bacteria classified by multi-locus sequence analysis; proposal of Xanthomonas dyei sp. nov. Plant Pathol. 2010, 59: 270-281. 10.1111/j.1365-3059.2009.02210.x.
Aritua V, Parkinson NM, Thwaites R, et al: Characterization of the Xanthomonas sp. causing wilt of enset and banana and its proposed reclassification as a strain of X. vasicola. Plant Pathol. 2008, 57: 170-177.
Bui Thi Ngoc L, Vernière C, Jouen E, et al: Amplified fragment length polymorphism and multilocus sequence analysis-based genotypic relatedness among pathogenic variants of Xanthomonas citri pv. citri and Xanthomonas campestris pv. bilvae. Int J Syst Evol Microbiol. 2010, 60: 515-525. 10.1099/ijs.0.009514-0.
Rademaker JLW, Norman DJ, Forster RL, et al: Classification and identification of Xanthomonas translucens isolates, including those pathogenic to ornamental asparagus. Phytopathology. 2006, 96: 876-884. 10.1094/PHYTO-96-0876.
Valverde A, Hubert T, Stolov A, et al: Assessment of genetic diversity of Xanthomonas campestris pv. campestris isolates from Israel by various DNA fingerprinting techniques. Plant Pathol. 2007, 56: 17-25.
Vicente JG, Everett B, Roberts SJ: Identification of isolates that cause a leaf spot disease of brassicas as Xanthomonas campestris pv. raphani and pathogenic and genetic comparison with related pathovars. Phytopathology. 2006, 96: 735-745. 10.1094/PHYTO-96-0735.
Sawada H, Kunugi Y, Watauchi K, Kudo A, Sato T: Bacterial spot, a new disease of grapevine (Vitis vinifera) caused by Xanthomonas arboricola. Jpn J Phytopathol. 2011, 77: 7-22. 10.3186/jjphytopath.77.7.
Schaad NW, Postnikova E, Lacy GH, et al: Reclassification of Xanthomonas campestris pv. citri (ex Hasse 1915) Dye 1978 forms A, B/C/D, and E as X. smithii subsp. citri (ex Hasse) sp. nov. nom. rev. comb. nov., X. fuscans subsp. aurantifolii (ex Gabriel 1989) sp. nov. nom. rev. comb. nov., and X. alfalfae subsp. citrumelo (ex Riker and Jones) Gabriel et al., 1989 sp. nov. nom. rev. comb. nov.; X. campestris pv. malvacearum (ex Smith 1901) Dye 1978 as X. smithii subsp. smithii nov. comb. nov. nom. nov.; X. campestris pv. alfalfae (ex Riker and Jones, 1935) Dye 1978 as X. alfalfae subsp. alfalfae (ex Riker et al., 1935) sp. nov. nom. rev.; and "var. fuscans" of X. campestris pv. phaseoli (ex. Smith, 1987) Dye 1978 as X. fuscans subsp. fuscans sp. nov. Syst Appl Microbiol. 2005, 28: 494-518. 10.1016/j.syapm.2005.03.017.
Schaad NW, Postnikova E, Lacy GH, et al: Emended classification of xanthomonad pathogens on citrus. Syst Appl Microbiol. 2006, 29: 690-695. 10.1016/j.syapm.2006.08.001.
Ah-You N, Gagnevin L, Chiroleu F, et al: Pathological variations within Xanthomonas campestris pv. mangiferaeindicae support its separation into three distinct pathovars that can be distinguished by Amplified Fragment Length Polymorphism. Phytopathology. 2007, 97: 1568-1577. 10.1094/PHYTO-97-12-1568.
Fargier E, Manceau C: Pathogenicity assays restrict the species Xanthomonas campestris into three pathovars and reveal nine races within X. campestris pv. campestris. Plant Pathol. 2007, 56: 805-818. 10.1111/j.1365-3059.2007.01648.x.
Jones JB, Lacy GH, Bouzar H, Stall RE, Schaad NW: Reclassification of the xanthomonads associated with bacterial spot disease of tomato and pepper. Syst Appl Microbiol. 2004, 27: 755-762. 10.1078/0723202042369884.
Young JM, Park D-S, Shearman HM, Fargier E: A multilocus sequence analysis of the genus Xanthomonas. Syst Appl Microbiol. 2008, 31: 366-377. 10.1016/j.syapm.2008.06.004.
Gonçalves ER, Rosato YB: Phylogenetic analysis of Xanthomonas species based upon 16S-23S rDNA intergenic spacer sequences. Int J Syst Evol Microbiol. 2002, 52: 355-361.
Hauben L, Vauterin L, Swings J, Moore ER: Comparison of 16S ribosomal DNA sequences of all Xanthomonas species. Int J Syst Bacteriol. 1997, 47: 328-335. 10.1099/00207713-47-2-328.
Moore ER, Krüger AS, Hauben L, et al: 16S rRNA gene sequence analyses and inter- and intrageneric relationships of Xanthomonas species and Stenotrophomonas maltophilia. FEMS Microbiol Lett. 1997, 151: 145-153. 10.1111/j.1574-6968.1997.tb12563.x.
Parkinson NM, Cowie C, Heeney J, Stead DE: Phylogenetic structure of Xanthomonas determined by comparison of gyrB sequences. Int J Syst Evol Microbiol. 2009, 59: 264-274. 10.1099/ijs.0.65825-0.
Deloger M, El Karoui M, Petit M-A: A genomic distance based on MUM indicates discontinuity between most bacterial species and genera. J Bacteriol. 2009, 191: 91-99. 10.1128/JB.01202-08.
Richter M, Rosselló-Móra R: Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci USA. 2009, 106: 19126-19131. 10.1073/pnas.0906412106.
Konstantinidis KT, Tiedje JM: Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci USA. 2005, 102: 2567-2572. 10.1073/pnas.0409727102.
Rokas A, Williams BL, King N, Carroll SB: Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature. 2003, 425: 798-804. 10.1038/nature02053.
Philippe H, Delsuc F, Brinkmann H, Lartillot N: Phylogenomics. Annu Rev Ecol Evol Syst. 2005, 36: 541-562. 10.1146/annurev.ecolsys.35.112202.130205.
Wu M, Eisen JA: A simple, fast, and accurate method of phylogenomic inference. Genome Biol. 2008, 9: R151-10.1186/gb-2008-9-10-r151.
Pieretti I, Royer M, Barbe V, et al: The complete genome sequence of Xanthomonas albilineans provides new insights into the reductive genome evolution of the xylem-limited Xanthomonadaceae. BMC Genomics. 2009, 10: 616-10.1186/1471-2164-10-616.
Qian W, Jia Y, Ren S, et al: Comparative and functional genomic analyses of the pathogenicity of phytopathogen Xanthomonas campestris pv. campestris. Genome Res. 2005, 15: 757-767. 10.1101/gr.3378705.
da Silva A, Ferro J, Reinach F, et al: Comparison of the genomes of two Xanthomonas pathogens with differing host specificities. Nature. 2002, 417: 459-463. 10.1038/417459a.
Vorhölter F, Schneiker S, Goesmann A, et al: The genome of Xanthomonas campestris pv. campestris B100 and its use for the reconstruction of metabolic pathways involved in xanthan biosynthesis. J Biotechnol. 2008, 134: 33-45. 10.1016/j.jbiotec.2007.12.013.
Thieme F, Koebnik R, Bekel T, et al: Insights into genome plasticity and pathogenicity of the plant pathogenic bacterium Xanthomonas campestris pv. vesicatoria revealed by the complete genome sequence. J Bacteriol. 2005, 187: 7254-7266. 10.1128/JB.187.21.7254-7266.2005.
Studholme DJ, Kemen E, MacLean D, et al: Genome-wide sequencing data reveals virulence factors implicated in banana Xanthomonas wilt. FEMS Microbiol Lett. 2010, 310: 182-192. 10.1111/j.1574-6968.2010.02065.x.
Lee B, Park Y, Park D, et al: The genome sequence of Xanthomonas oryzae pathovar oryzae KACC10331, the bacterial blight pathogen of rice. Nucleic Acids Res. 2005, 33: 577-586. 10.1093/nar/gki206.
Ochiai H, Inoue Y, Takeya M, et al: Genome sequence of Xanthomonas oryzae pv. oryzae suggests contribution of large numbers of effector genes and Insertion Sequences to its race diversity. JARQ. 2005, 39: 275-287.
Salzberg S, Sommer D, Schatz M, et al: Genome sequence and rapid evolution of the rice pathogen Xanthomonas oryzae pv. oryzae PXO99A. BMC Genomics. 2008, 9: 204-10.1186/1471-2164-9-204.
Hötte B, Rath-Arnold I, Pühler A, Simon R: Cloning and analysis of a 35.3-kilobase DNA region involved in exopolysaccharide production by Xanthomonas campestris pv. campestris. J Bacteriol. 1990, 172: 2804-2807.
Kamoun S, Kado CI: Phenotypic switching affecting chemotaxis, xanthan production, and virulence in Xanthomonas campestris. Appl Environ Microbiol. 1990, 56: 3855-3860.
Restrepo S, Duque MC, Verdier V: Characterization of pathotypes among isolates of Xanthomonas axonopodis pv. manihotis in Colombia. Plant Pathol. 2000, 49: 680-687. 10.1046/j.1365-3059.2000.00513.x.
Mew TW, Cruz Vera CM, Medalla ES: Changes in race frequency of Xanthomonas oryzae pv. oryzae in response to rice cultivars planted in the Philippines. Plant Dis. 1992, 76: 1029-1032. 10.1094/PD-76-1029.
Simpson AJ, Reinach FC, Arruda P, et al: The genome sequence of the plant pathogen Xylella fastidiosa. The Xylella fastidiosa Consortium of the Organization for Nucleotide Sequencing and Analysis. Nature. 2000, 406: 151-159. 10.1038/35018003.
Monteiro-Vitorello CB, de Oliveira MC, Zerillo MM, et al: Xylella and Xanthomonas Mobil'omics. OMICS. 2005, 9: 146-159. 10.1089/omi.2005.9.146.
Didelot X, Darling ACE, Falush D: Inferring genomic flux in bacteria. Genome Res. 2009, 19: 306-317.
Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.
Atmakuri K, Cascales E, Christie PJ: Energetic components VirD4, VirB11 and VirB4 mediate early DNA transfer reactions required for bacterial type IV secretion. Mol Microbiol. 2004, 54: 1199-1211. 10.1111/j.1365-2958.2004.04345.x.
Kuldau GA, De Vos G, Owen J, McCaffrey G, Zambryski P: The virB operon of Agrobacterium tumefaciens pTiC58 encodes 11 open reading frames. Mol Gen Genet MGG. 1990, 221: 256-266.
Hu SH, Peek JA, Rattigan E, Taylor RK, Martin JL: Structure of TcpG, the DsbA protein folding catalyst from Vibrio cholerae. J Mol Biol. 1997, 268: 137-146. 10.1006/jmbi.1997.0940.
Langille MGI, Hsiao WWL, Brinkman FSL: Evaluation of genomic island predictors using a comparative genomics approach. BMC Bioinforma. 2008, 9: 329-10.1186/1471-2105-9-329.
Euzéby JPM: List of Prokaryotic names with Standing in Nomenclature. [http://www.bacterio.cict.fr/index.html]
Barton NH, Briggs DEG, Eisen JA, Goldstein DB, Patel NH: Phylogenetic Reconstruction. Evolution. 2007, New York: Cold Spring Harbo Laboratory Press
Stajich JE, Block D, Boulez K, et al: The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 12: 1611-1618. 10.1101/gr.361602.
Vos RA, Caravas J, Hartmann K, Jensen MA, Miller C: Bio::Phylo-phyloinformatic analysis using Perl. BMC Bioinforma. 2011, 12: 63-10.1186/1471-2105-12-63.
Fitch WM: Uses for evolutionary trees. Philos Trans R Soc Lond B Biol Sci. 1995, 349: 93-102. 10.1098/rstb.1995.0095.
Simmons MP, Donovan Bailey C, Nixon KC: Phylogeny reconstruction using duplicate genes. Mol Biol Evol. 2000, 17: 469-473.
Huson DH, Steel M: Phylogenetic trees based on gene content. Bioinformatics (Oxford, England). 2004, 20: 2044-2049. 10.1093/bioinformatics/bth198.
Dawyndt P, Vancanneyt M, De Meyer H, Swings J: Knowledge accumulation and resolution of data inconsistencies during the integration of microbial information sources. IEEE Trans Knowl Data Eng. 2005, 17: 1111-1126.
Delcher AL, Bratke KA, Powers EC, Salzberg SL: Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics (Oxford, England). 2007, 23: 673-679. 10.1093/bioinformatics/btm009.
Altschul SF, Madden TL, Schäffer AA, Zhang J: Gapped BLAST and PSI-BLAST: a new generation of protein database. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Koski LB, Golding GB: The closest BLAST hit is often not the nearest neighbor. J Mol Evol. 2001, 52: 540-542.
Moreno-Hagelsieb G, Latimer K: Choosing BLAST options for better detection of orthologs as reciprocal best hits. Bioinformatics (Oxford, England). 2008, 24: 319-324. 10.1093/bioinformatics/btm585.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Bruen TC, Philippe H, Bryant D: A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006, 172: 2665-2681.
Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.5c. 1993
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics (Oxford, England). 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.
Rannala B, Yang Z: Probability distribution of molecular evolutionary trees: a new method of phylogenetic inference. J Mol Evol. 1996, 43: 304-311. 10.1007/BF02338839.
Yang Z, Rannala B: Bayesian phylogenetic inference using DNA sequences: a Markov Chain Monte Carlo Method. Mol Biol Evol. 1997, 14: 717-724.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics (Oxford, England). 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Swofford DL: PAUP*. 2002
Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-1120. 10.1007/BF01731581.
Darling ACE, Mau B, Perna NT: progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010, 5: e11147-10.1371/journal.pone.0011147.
Darling ACE, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14: 1394-1403. 10.1101/gr.2289704.
Tatusov RL, Fedorova ND, Jackson JD, et al: The COG database: an updated version includes eukaryotes. BMC Bioinforma. 2003, 4: 41-10.1186/1471-2105-4-41.
Sayers EW, Barrett T, Benson DA, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2009, 37: D5-D15. 10.1093/nar/gkn741.
This project was funded by the Colombian administrative department of Science, Technology and Innovation (Colciencias) and the Vice-chancellor's Office of Research at the Universidad de Los Andes. We would like to thank Andrew Crawford, Ralf Koebnik and two anonymous reviewers for critical reading of the manuscript. We also thank Boris Szurek, Valérie Verdier, Kostantinos Konstantinidis, Catalina Arévalo and Camilo López for comments and discussion on the conception and development of this study.
LMR participated in the design and coordination of the study, acquired data, carried out the analysis and drafted the manuscript. AG participated in the design and coordination of the study, acquired data and critically revised the manuscript. MLA participated in the design and coordination of the analyses. CS participated in the design and coordination of the study and critically revised the manuscript, while SR participated in the design and coordination and critically revised the manuscript. AB conceived the study, participated in the design and coordination of the study, drafted and critically revised the manuscript. All authors read and approved the final manuscript.