- Research article
- Open Access
Evolutionary, structural and functional relationships revealed by comparative analysis of syntenic genes in Rhizobiales
BMC Evolutionary Biology volume 5, Article number: 55 (2005)
Comparative genomics has provided valuable insights into the nature of gene sequence variation and chromosomal organization of closely related bacterial species. However, questions about the biological significance of gene order conservation, or synteny, remain open. Moreover, few comprehensive studies have been reported for rhizobial genomes.
We analyzed the genomic sequences of four fast growing Rhizobiales (Sinorhizobium meliloti, Agrobacterium tumefaciens, Mesorhizobium loti and Brucella melitensis). We made a comprehensive gene classification to define chromosomal orthologs, genes with homologs in other replicons such as plasmids, and those which were species-specific. About two thousand genes were predicted to be orthologs in each chromosome and about 80% of these were syntenic. A striking gene colinearity was found in pairs of organisms and a large fraction of the microsyntenic regions and operons were similar. Syntenic products showed higher identity levels than non-syntenic ones, suggesting a resistance to sequence variation due to functional constraints; also, an unusually high fraction of syntenic products contained membranal segments. Syntenic genes encode a high proportion of essential cell functions, presented a high level of functional relationships and a very low horizontal gene transfer rate. The sequence variability of the proteins can be considered the species signature in response to specific niche adaptation. Comparatively, an analysis with genomes of Enterobacteriales showed a different gene organization but gave similar results in the synteny conservation, essential role of syntenic genes and higher functional linkage among the genes of the microsyntenic regions.
Syntenic bacterial genes represent a commonly evolved group. They not only reveal the core chromosomal segments present in the last common ancestor and determine the metabolic characteristics shared by these microorganisms, but also show resistance to sequence variation and rearrangement, possibly due to their essential character. In Rhizobiales and Enterobacteriales, syntenic genes encode a high proportion of essential cell functions and presented a high level of functional relationships.
A huge amount of information has been obtained from sequencing projects. More than two hundred complete microbial genomes are available to date in public databases and sequencing of a similar number is in progress . Many questions remain unsolved. For example, what is the biological meaning, if any, of gene arrangement in the bacterial chromosome?
Changes in gene sequence and chromosomal rearrangements constitute the main sources of genomic variability. Nonsynonymous substitutions in the first or second nucleotides of the codon change the encoded residue and are thus a driving force of natural selection. Genomic studies in bacteria regarding synonymous and nonsynonymous substitution rates have been published elsewhere [2, 3]. Chromosomes show constraints on rearrangement and works dealing with that aspect were recently reviewed by Rocha ; he suggested that there is a balance between conservation and change in the organization of the chromosome.
The operon represents the first level of the gene organization. Neighboring genes, especially those in co-directional and in divergent orientation, represent a second organization level because they show a certain functional association revealed by genomic context analysis . Regarding comparisons among closely related species, the gene order conservation, or synteny, represents a third level of organization. Synteny depends on shared ancestry and inter- and intrachromosomal exchanges, and represents a higher relationship between taxa. It was suggested that physiologically important gene clusters could be positively selected, and synteny perhaps reveals the functional constraints of these genes . For the detection of synteny it is necessary first to determine the set of orthologous genes in pairs of organisms and recently an inverse method has proven useful for this .
Recombination/transposition events can easily disrupt synteny. Species of Buchnera and Corynebacterium have low levels of chromosomal rearrangements and lack recA and recBCD orthologs, respectively [8, 9], thus suggesting that recombination is an important factor for loss of synteny. Synteny studies have focused in short gene clusters in eukaryotes [10, 11] while whole chromosome comparative analysis has been done in bacteria and archaea [12–18]. For example, there is a striking conservation between the chromosomes of Escherichia coli and Salmonella typhimurium , and also between those of Brucella melitensis, B. suis and Mesorhizobium loti . However, the synteny analyses for alpha proteobacteria such as those reported in the genome sequence determinations of Brucella melitensis, Sinorhizobium meliloti and Agrobacterium tumefaciens are highly schematic [20–22]. A thorough analysis of Rhizobiales genomes would determine if the key set of genes covering the most important metabolic functions in these organisms are syntenic.
Another factor affecting chromosomal rearrangement is horizontal gene transfer (HGT) which occurs in bacteria [23, 24], but estimating its impact on genome organization has proved a daunting task for two reasons. First, the reliability of compositional methods to detect HGT events has been questioned [25, 26]. Second, phylogenetic methods, albeit more reliable, are not always applicable and can easily be misleading without proper care. The results of a recently published analysis suggest that codon usage compatibility between alien genes and recipient genomes  is a prerequisite for successful HGT events. This premise has been supported by other evidence [28, 29].
Among the most accepted methods to deduce functional relationships of proteins are phylogenetic profile , gene neighboring , and the Rossetta stone method . These methods can give additive information about metabolic networks existing in organisms . ProLinks is a program based on these methods with an extensive library of predicted functional interactions from 83 genomes . Von Mering et al. , also applying these approaches with the STRING program, found global modularities in functional protein networks. A question remains about whether functional linkage differences exist in syntenic and non-syntenic gene clusters.
Rhizobiales is a prokaryotic order belonging to the alpha proteobacteria subdivision; some rhizobial species are intensively studied for their nitrogen-fixing ability when in symbiosis with leguminous plants. The order comprises both plant symbionts and plant and animal pathogens such as Rhizobium, Agrobacterium and Brucella, respectively. In rhizobia, genes responsible for the symbiotic interaction are commonly found on large plasmids or incorporated in a particular stretch of the chromosome called the symbiotic island [36–38]. The physiological potential of the rhizobial chromosome allows cell survival under different conditions. For example, an A. tumefaciens strain containing the symbiotic plasmid from Rhizobium etli induced nodules on legume plants [39, 40], and conversely, an S. meliloti derivative strain with Ri, the rhizogenic induction plasmid, formed root mats on alfalfa plants . Additionally, there are similarities in the parasitic/symbiotic strategies employed by species of the Rhizobiales . Also, it is possible to find diverse life-styles among the members of Enterobacteriales order (gamma proteobacteria): for example, Buchnera is an obligate aphid symbiont; E. coli and S. typhimurium are common gut inhabitants in mammals; and Shigella flexneri, Yersinia pestis and Erwinia carotovora are pathogens, either for animals or plants .
The complete genome sequences of seven species of Rhizobiales were available by 2004, namely S. meliloti , Mesorhizobium loti , Bradyrhizobium japonicum , A. tumefaciens [22, 45], B. melitensis , B. suis  and Rhodopseudomonas palustris . Although their genomes show a certain degree of conservation, variability corresponding to their evolutionary divergence points, microbial life styles and ecological niches was also found. Comparative genomics has captured the attention of researchers as a way of achieving a better understanding of the molecular basis underlying phenomena such as symbiosis and pathogenesis.
We classified the genes of several rhizobial species in order to gain a comprehensive insight into chromosomal conservation and genome rearrangement. Conserved genes among these species can reveal phylogenetic relationships, but also show metabolic strategies useful in understanding the niche diversity in which these organisms usually grow. In particular, syntenic/non-syntenic genes among these species were analyzed in terms of their sequence identity/similarity, physical characteristics of the encoded products and functional relationships among them. Additionally, in order to find more general trends, we compared these results with an analysis performed on genomes belonging to the Enterobacteriales.
Approach, strategy and outline
Our main objective was to enhance our understanding of the functional meaning of the gene arrangement on the bacterial chromosome, taking as examples some genomes from the Rhizobiales and Enterobacteriales. We consider that gene neighboring is not a random trait and gives an adaptive advantage to the cell because the proteins produced are likely to perform related functions. Our belief is that the coordinated expression of genes, organized on the chromosome either as operons or clusters, permits the correct integration of metabolic functions.
Our approaches were: i) to obtain a comprehensive gene classification, applicable to each of the species analyzed and suitable to make comparisons among them, and ii) to detect specific gene characteristics (if any) of each of the classes. In the first approach we identified orthologs among chromosomes, defined those that were syntenic, those in a different replicon, and those that were species-specific. For the second approach, we analyzed gene/protein sequences for identity, calculated the horizontal gene transfer rate for each class and the predicted molecular weight and isoelectric point of the peptides, and inferred the functional relationship in syntenic or non-conserved chromosomal regions. The results are presented in the following order: 1) Synteny in Rhizobiales (gene classification of Rhizobiales; gene organization and microsyntenic region formation; synteny and insertion sequences, horizontal gene transfer and codon usage), 2) Synteny in Enterobacteriales, 3) Sequence analysis of the chromosomal predicted orthologs, 4) Physical characteristics of the translated products of syntenic genes, and 5) Functional roles and linkage of chromosomal predicted orthologs. S. meliloti was taken as reference organism for the comparisons with each A. tumefaciens, M. loti, B. melitensis and E. coli. E. coli was used as base to compare with S. typhimurium, E. carotovora and S. meliloti.
1) Synteny in Rhizobiales
Gene classification of Rhizobiales
The selection criteria mentioned in Methods were applied to the chromosomes of S. meliloti (Sm), A. tumefaciens (At), M. loti (Ml), and B. melitensis (Bm). As compared to the chromosome of S. meliloti, we found that more than 60% of genes were chromosomal predicted orthologs in the At circular chromosome (At-C) and Bm chromosome I (Bm-I), one third in Ml chromosome and Bm chromosome II (Bm-II) and one quarter in the At linear chromosome (At-L); however, the number of chromosomal predicted orthologs in each organism was similar, about two thousand (Table 1). The Sm chromosome presents 3341 genes in 3.65 Mb.
We assessed the chromosomal genes with conserved order or synteny (see Methods). Syntenic genes represented about 70–80% of chromosomal predicted orthologs (see Table 1). That is, the conserved chromosomal order of these genes seems favored. The remarkable synteny level is highlighted by a group of 1038 common syntenic genes in all these species. Non-syntenic genes represented from 17% to 35% of the chromosomal predicted orthologs in these organisms. Only 98 non-syntenic genes were common in the four Rhizobiales.
The remaining categories obtained in this analysis were: homologs present in plasmids, homologs with the rest of rhizobial chromosomes, and those with no orthologs in the public databases (species-specific genes, also known as orphan). Homologs in plasmids were more abundant in the At-L,Ml and Bm-II chromosomes (Table 1). These replicons have special features as commented below and in Discussion. Aside from the two organisms being compared, some genes also matched with other rhizobial chromosomes (including B. japonicum USDA110) and comprised about a quarter of the chromosomal genes in these species (Table 1). Some of them only matched with unidirectional best hits (homologs). In this context, in an additional analysis between S. meliloti and A. tumefaciens, we found that 129 genes matched with unidirectional best hit (white fraction, Figure 1), but the rest were predicted orthologs with bidirectional best hits, either syntenic or non-syntenic, and are represented in the red and blue striped bars (denoted as rest, Figure 1), respectively. Species-specific genes were especially abundant in the M. loti chromosome, while in the other replicons this class covered at most 20%. One half of these species-specific genes were not present in the COG database  and the rest were denoted as hypothetical (data not shown).
In this way, all chromosomal genes were assigned to the categories mentioned, as shown in Figure 1 for the comparison of S. meliloti with A. tumefaciens. This approach gives a panoramic view about shared and unshared genes with other rhizobial species. Additional file 1 shows the categories for the other comparisons.
A schematic representation of rhizobial chromosomes was obtained relative to the classification of the genes. A high proportion of chromosomal predicted orthologs of the analyzed species were syntenic; however, in the A. tumefaciens linear, M. loti and B. melitensis I chromosomes, plasmidic homologs and species-specific genes were particularly abundant.
Gene organization in operons and microsyntenic regions in Rhizobiales
A relationship between predicted orthologs and operons was also explored. For S. meliloti (in comparison with A. tumefaciens), one half of the syntenic genes was found organized in 303 syntenic operons, a half of the total predicted operons (606); similar proportions were found for the A. tumefaciens circular and B. melitensis I chromosomes. In the A. tumefaciens linear, B. melitensis II and M. loti chromosomes, the proportion was about 21% (Table 1). The first two can be considered as accessory chromosomes and the last is the largest chromosome of the analyzed species. Non-syntenic genes were found organized in operons in a very small proportion in all analyzed Rhizobiales (Table 1).
A relevant level of operon conservation was found. In the comparison of S. meliloti with A. tumefaciens, 50% of the syntenic genes organized in operons were in identical operons, and taken together with those in similar operons (differing by one or two genes), 82% of total syntenic genes were in conserved operons. Similar proportions were found for the other comparisons (data not shown).
The operons formed with plasmidic homologs constituted a small fraction of the predicted operons (Table 1). Also, a reduced proportion was found for species-specific operons, except for Ml. Finally, mixed operons were present in a higher amount and ranged from 37 to 64% of the predicted operons (Table 1). This category revealed the highest rate of chromosomal rearrangements in these organisms. The mixed operons contained 17% of the syntenic genes, on average.
Syntenic genes were found in clusters and were assigned to microsyntenic regions (see Methods). When comparing S. meliloti with both A. tumefaciens chromosomes, 205 regions were common. In particular, At-C regions in common with Sm chromosome were found along all the chromosome, except in the third quarter (Figure 2). The third quarter had colinearity with At-L (Figure 2, see also Figure 3 panel a). In the other comparisons, a similar amount of common regions were observed (Table 1, see also Additional files 2 and 3). 146 regions were shared in Sm, At, and Ml and 94 regions were common to the four Rhizobiales. About 90% of syntenic genes were located in the microsyntenic regions (Table 1).
When At-C syntenic genes were compared with Sm chromosome, a high level of colinearity was found (see Additional file 4 panel a). In the case of Ml, the synteny was disrupted possibly due to a conflicting annotation; for Bm-I, an inverse colinearity was obtained, possibly by oriC inversion relative to Sm (see Additional file 4 panels b and c, respectively). In regard to the rearrangement of microsyntenic regions, an extensive chromosomal colinearity was observed in long tracts. For example, in the comparison of Sm and At circular chromosomes, 93 microsyntenic regions were colinear, 24 almost colinear, 18 with drastic changes, and 13 inverted. The schematic representation is shown in Figure 3 panel a, lower part. In the comparison with the At linear and Bm-II chromosomes, highly rearranged structures were found (Figure 3 panel a, upper part and panel c, upper part). Rearrangement of microsyntenic regions on the chromosomes of M. loti (Figure 3 panel b) and B. melitensis I (Figure 3 panel c, lower part), showed a high level of colinearity when the conflicting annotation and the oriC inversion, respectively, were modified (see Methods).
Syntenic genes were organized mainly at two levels: operons and microsyntenic regions. Syntenic operons were as abundant as mixed operons. Extensive blocks of chromosomal colinearity were found, despite rearrangements.
Synteny and insertion sequences, horizontal gene transfer and codon usage
The mobile elements play an important role in chromosomal rearrangement. To determine how these elements were dispersed among microsyntenic regions, insertion sequence (IS) and transposase locations were analyzed. The S. meliloti chromosome contains 51 IS and 68 transposases belonging to diverse families . Of the total transposases, 40 were found with homologs in plasmids, 20 were common with other rhizobial chromosomes, 6 were denoted as species-specific and 2 were common with the A. tumefaciens linear chromosome. Only 14 pairs of IS/transposases (27% of the total) were found inside microsyntenic regions.
We assessed the influence of horizontal gene transfer (HGT) on the genomic structure of rhizobial genomes (see Methods). Table 2 shows calculated HGT events for each of the classes syntenic, non-syntenic and plasmid homologs. Even though syntenic genes were in the largest class analyzed, they displayed the lowest number of HGT events (Table 2). Predicted HGT rates of non-syntenic genes were 2 to 8 times higher than those of syntenic genes. On the other hand, unlike any other gene class, species-specific genes had the strongest bias toward a low codon richness index in the four Rhizobiales (denoted with crosses, Figure 4 panels a to d). No significant differences were found for syntenic or non-syntenic genes in each organism.
2) Synteny in Enterobacteriales (gamma proteobacteria)
To assess the adequacy of applying our synteny analysis to another bacterial clade, we chose two members from the Enterobacteriales (gamma proteobacteria), the closely related E. coli and S. typhimurium genomes and defined their orthologous genes. These organisms contained 3092 predicted orthologs and 95% of them (2943) were syntenic. An extensive chromosomal colinearity with few rearrangements was found (data not shown).
Also, we chose a more phylogenetically distant species, Erwinia carotovora subsp. atroseptica to compare with E. coli. Their genomes comprise 4254 and 4477 genes, respectively. They shared 2477 orthologous genes and these represented about half the total genes in each chromosome. In the detection of synteny between these organisms, 1993 genes (80.4% of the orhologs) fulfilled our requirement and the rest, 484, were classified as non-syntenic. When the genes were assigned to microsyntenic regions, 230 regions were found and contained 92.8% (1849) of total syntenic genes and the rest, 7.2% (144 genes), were detected in the non-conserved tracts (see Additional file 5). 172 non-orthologous genes also formed part of microsyntenic regions. The 230 non-conserved regions contained 2233 genes.
To define the conservation of orthology and synteny of Enterobacteriales genomes in comparison with the Rhizobiales, we compared the S. meliloti and E. coli chromosomes. We found 777 predicted orthologs between them, a proportion that represents only a third of the genes shared in the Rhizobiales. By visual examination no synteny was apparent between the Sm and E. coli chromosomes (see Additional file 4 panel d). With an algorithm, only 198 syntenic genes were assigned to 65 microsyntenic regions.
3) Sequence analysis of the chromosomal predicted orthologs
To determine the level of sequence identity among chromosomal predicted orthologs of the Rhizobiales, distribution curves of their translation products were obtained. Two different types of curves were found. By comparing S. meliloti with both A. tumefaciens chromosomes, the syntenic products presented a Gaussian distribution, with a tendency to high identity levels (asymmetry coefficient g1 = -0.48, significant at P < 0.001) (Figure 5 panel a) and a mean value of 71.3%. The non-syntenic products showed a non- asymmetric distribution (g1 = -0.09) with a lower mean value (61.9%) (Figure 5 panel a). Similar curves were found for the chromosomal predicted orthologs of the other comparisons (Additional file 6 panels a and b, for Sm-Ml and Sm-Bm comparisons, respectively). When At replicons were separately compared with Sm, syntenic products of both At-C or At-L showed a similar bias toward high identity levels; in contrast, non-syntenic products of At-L showed a strongly deviated distribution to low identity values (data not shown). The tendency of syntenic products to higher identity levels reflects not only restriction to change but also functional constraints, possibly due to an essential character. Conversely, the lower identity levels of non-syntenic products represent lower restrictions to change and higher functional versatility.
A comparison of E. coli and S. typhimurium genomes (belonging to the Enterobacteriales) was performed. A very asymmetric distribution curve with a tendency to high identity levels for syntenic products was obtained (Figure 5 panel b); it remarkably resembled that obtained in the comparisons of syntenic products of the Rhizobiales. Interestingly, in the comparison of S. meliloti with E. coli both chromosomal non-syntenic and syntenic products had asymmetric curves with strong tendency to low identity levels (Figure 5, panel c).
To obtain a complete view of sequence variation of chromosomal predicted orthologs in the four Rhizobiales, either syntenic or non-syntenic genes were graphed in relation to their identity levels in comparison with S. meliloti. Figure 6 panel a shows the sequence identity of 1038 common syntenic genes in the four Rhizobiales, for Sm-At, Sm-Ml and Sm-Bm comparisons. We referred the comparison to the identity of Sm-Ml genes ordered progressively. We found that each comparison showed a particular clustering, probably related to the phylogenetic distances between these organisms. The Pearson correlation coefficients were r = 0.81 and r = 0.88 for Sm-Bm and Sm-At comparisons, respectively. In the clustering of 98 non-syntenic genes common to the four Rhizobiales (data not shown), lower correlation coefficients were obtained (r = 0.71 and r = 0.79, for Sm-Bm and Sm-At comparisons, respectively). In Figure 6 panel b the sequence identity of 140 common syntenic genes in the Rhizobiales and Enterobacteriales are shown. While the comparison for Rhizobiales follows a similar tendency to the previously observed, for E. coli and Salmonella there is a higher relatedness level. However, when comparing S. meliloti and E. coli a different tendency with very low identity level is observed. A high level of clustering among closely related species belonging to the same subdivision is visible and also, a tendency to reduced relatedness in members of different subdivisions; that is, synteny is a common trait for members of each subdivision.
To determine the meaning of sequence differences we analyzed the translated syntenic product ArgC, which participates in the arginine biosynthetic pathway. The alignment presented in Additional file 7 panel a shows 111 positions with identical residues and a range from 33 to 113 different residues, particular for each of the species compared. However, sequences from Brucella melitensis and Brucella suis showed only one difference between them (not shown). Synteny is almost absolute in these organisms (Additional file 4 panel e). Changing residues (possibly species signatures) were dispersed along the sequences and varied according to the identity level.
To determine more comprehensively sequence differences and similarities in Rhizobiales (alpha proteobacteria) and Enterobacteriales (gamma proteobacteria), we selected five organisms from each: S. meliloti, A. tumefaciens, M. loti, B. melitensis and R. palustris for the first and E. coli, S. typhimurium, S. flexneri, Buchnera sp. and E. carotovora for the last. We chose some syntenic products from the arginine biosynthetic pathway, namely ArgB, ArgC, ArgD, ArgF, ArgG, and ArgH. The alignment belonging to ArgC is shown in Additional file 7 panel b. As can be observed, there are sequence identities and differences among species from the same order ("species signatures"), but also similarities between species of different orders, albeit at smaller level. Additionally, we found an interesting pattern: proteins from Enterobacteriales showed an almost uniform level of sequence identity and similarity (about 50 and 70%, on average, respectively), however, sequences from Rhizobiales showed a clear increasing tendency, from 25 to 62% in identity, and from 44 to 81% in similarity (see Additional file 7 panel c). These different profiles possibly are related with the particular conditions of the niches occupied by these organisms.
4) Physical characteristics of the translated syntenic genes
Molecular weight (MW) and isoelectric point (pI) are the main traits for proteomic comparisons. To determine whether syntenic and non-syntenic genes could present differential protein characteristics, we graphed pairs of translated predicted orthologs for both MW and pI parameters. In all comparisons, MW graphs showed lower dispersion than pI ones. Figure 7 shows the pI graph for the S. meliloti-A. tumefaciens comparison (Sm-M. loti and Sm-B. melitensis comparisons are in Additional file 8 panels a and b, respectively). A large group of proteins (82% in each comparison) was located on the diagonal. The correlation coefficient for this group was r = 0.96 in all comparisons. However, the rest of the predicted proteins showed differential pI's and were assigned to sectors (Figure 7). Sector I had acidic proteins in S. meliloti and basic proteins in the organism compared with it. Sector II had basic proteins in S. meliloti and acidic in organism compared to it. In sector III were neutral proteins in S. meliloti and covered all pI range in each of the organisms compared. Proteins in sector IV were neutral in the comparing organism and covered all pI ranges in S. meliloti. In Additional file 9 there is a summary of pI variability of common syntenic products from comparisons with the chromosomes of S. meliloti, A. tumefaciens and M. loti. About 75% of products showed similar pI and 14% and 8% presented high and low variation, respectively. High level is defined as pI variation from acid in one or two organisms and basic in the other, and viceversa. Low level is defined as pI variation from neutral to acid or basic (Additional file 9). Proteins of the sectors mainly corresponded to the functional categories of energy generation, post-translational modification, and transport. In the case of non-syntenic products, a pattern similar to that described above was found (see Additional file 10). Proteins in the diagonal were the most conserved group with subtle pI changes possibly responding to species adaptation, whereas those with deviated pI's may represent a group with higher functional versatility.
Since the functional categories mentioned above for sectors are known to often interact with the cell membrane, a membrane prediction for all syntenic products from the S. meliloti-A. tumefaciens circular chromosomes comparison was assessed. Strikingly, 790 syntenic products were predicted to contain membranal segments. This amount represents almost all membranal proteins coded in the S. meliloti chromosome, considering that bacterial genomes have 18–28% membranal proteins [49–51]. About 70% of predicted membranal proteins with assigned function belonged to transport, energy generation, post-translational modification, cofactor synthesis, amino acid metabolism and central intermediary metabolism categories (Additional file 11). There are reports about membrane-interacting proteins with functions such as amino acid and cofactor biosynthesis and central intermediary metabolism [52–56].
To determine whether the charged amino acid residues were clustered in proteins from the sectors described above, we selected several proteins from each. The residues determining radical changes in pI were observed scattered along the sequences (data not shown).
5) Functional roles and linkage of the chromosomal predicted orthologs
To define the metabolic participation of the chromosomal predicted orthologs, a functional classification was made with chromosomal syntenic and non-syntenic genes. As shown in Figure 8, syntenic genes of S. meliloti-A. tumefaciens chromosomes contained a high proportion of the most important house-keeping functions. The relative proportion of non-syntenic genes grew with decreasing functional essentiality, for example transport and binding proteins, cellular processes and regulatory functions (Figure 8). Furthermore, the syntenic products covered 85% of the main metabolic pathways as defined in MetaCyc for S. meliloti (see Additional file 12). Similar results in functional coverage were obtained with the other comparisons (see Additional file 13 panels a and b). Common syntenic genes in the four Rhizobiales also included a large fraction of the house-keeping functions (lower segment of the red bars, Figure 8). In the case of the comparison between E. coli and E. carotovora, syntenic genes also covered a high proportion of essential functions, however, the first two positions were occupied by cofactor and nucleotide synthesis (see Additional file 13 panel c). In this regard, it is important to note that E. coli and E. carotovora possess no large plasmids.
When functional classes of genes in the microsyntenic regions were divided into Informational, Operational and Cellular processes superclasses and graphed for the S. meliloti chromosome with a 100 Kb window, an interesting pattern was observed (Figure 9), with the majority of the peaks belonging to a superclass matching with valleys of the other(s). This could represent functionally specialized blocks of chromosomal tracts, which were part of the ancestral rhizobial chromosome. For instance, the existence of genomic domains is accepted in eukaryotes .
The ProLinks program was used to determine how the chromosomal predicted orthologs are functionally related. Functional links were calculated for all genes in the S. meliloti and E. coli chromosomes and then correlated to their neighbors. For the S. meliloti-A.tumefaciens comparison, the microsyntenic regions presented, on average, 3.68 connections per node (a protein in a network), almost twice the value obtained in the non-conserved regions (2.06). From 205 microsyntenic regions, 104 had functional networks; in the case of non-conserved regions, only 35 presented networks. Networks with less than 6 connections were omitted. The networks of microsyntenic regions presented 1057 syntenic genes and from these 795 (75.2%) were organized in operons. 99 non-syntenic genes were in the networks of the non-conserved tracts, and only 56 (56.5%) were in operons. In the case of the synteny comparison between E. carotovora-E. coli (Enterobacteriales), 230 microsyntenic regions were obtained and from these, 161 presented functional networks, with a connectivity average of 5.33 connections/node. The non conserved regions with networks were 71 with a connectivity average of 2.04 connections/node. From 1497 syntenic genes in the networks, 1106 (73.8%) formed part of operons. Network connectivity obtained in S. meliloti and E. coli is shown in the graph of Additional file 14 (panels a and b, respectively). There is a striking difference in connectivity level in the networks from syntenic (gray bars) or non-conserved regions (black bars). The connectivity levels in syntenic vs non-conserved regions in both organisms were similar using the STRING program (data not shown).
The comparative genomic analysis reported here was useful in finding interesting gene properties. Orthologs with conserved replicon and neighborhood were the principal component of the chromosomes. Compared with non-syntenic genes, syntenic ones had higher identity levels, lower horizontal gene transfer (HGT) rates, showed strongly organized structures as operons and microsyntenic regions and a relative absence of mobile elements. Thus, the syntenic genes can be considered as the chromosomal backbone of the order. Plasmidic homologs were scattered on the chromosomes, and higher HGT rate and linkage to transposases support their extrachromosomal origin. Species-specific genes had the lowest codon richness index, and possibly were acquired in the evolutionary history of each of the species.
In this way, a rhizobial chromosomal origin can be envisioned. The chromosomal orthologs were the gene set derived from the common cenancestor. From these, syntenic genes conserved a relative chromosomal order (and operonic organization) and encode the essential functions of the cell; non-syntenic genes lost the clustering and possibly some came from HGT events. The plasmidic homologs were obtained possibly by mobilization throughout replicons, a nonrare process in the rhizobial phylogenetic branch. The species-specific genes represent the particular gene set of the species and are the most intriguing group due to their unknown functional roles and origin. Work with members of the last group will help to define traits not shared with other species.
The Rhizobiales species analyzed showed a striking proportion of orthologous genes, mainly chromosomal syntenic; non-syntenic genes were found in lower proportion. A large fraction of the first class was common in the four organisms. Syntenic genes had a strong tendency to form operons and almost all were clustered in microsyntenic regions. Additionally, these operons were conserved in pairs of organisms. Therefore, a strong restriction for chromosomal rearrangement is visible. Given that these organisms cover a wide spectrum of environmental distributions, from plant rhizosphere to animal host, the conserved chromosomal tracts may be important to determine the metabolic properties common to the order. Similar results were observed in the Enterobacteriales comparison. In a recent report, using computational inference, Boussau et al.  proposed a common ancestral set of about 3000 genes for proteobacterial genomes.
Although rhizobial genomes shared common traits, important differences were also observed. For example, abundance of plasmidic homologs and species-specific genes in A. tumefaciens (linear), B. melitensis II and M. loti chromosomes confirmed their complex evolutionary histories [20, 22, 59, 60]. The M. loti chromosome presents intensive incorporation of foreign genes by horizontal transfer, such as those belonging to the symbiotic island . In regard to species-specific genes, a large number presented low codon richness index. This category will be reduced with incorporation of other rhizobial genomes into the databases. For example, R. leguminosarum biovar viciae 3841, R. tropici PRF81, R. sp. ANU265 and R. etli CFN42 genomes soon will be available . Mobilization elements participate in chromosomal rearrangement and are abundant in Rhizobiales. These elements can decompose the microsyntenic regions where they are. In this way, all synteny approaches consist of snapshots in chromosomal evolution.
From the comparison of sequence identities among chromosomal predicted orthologs in the rhizobial species analyzed, an interesting characteristic was the differential distribution curves obtained. The asymmetric curves of syntenic products, deviated to high identity levels (with peaks at 70–75%, Figure 5 panel a and Additional file 6 panels a and b), possibly reveal their essentiality and can be compared with those from gamma proteobacteria with high identity (with peak at 95%, Figure 5 panel b). Conversely, non-syntenic genes had curves with lower sequence identity levels reflecting a higher functional versatility. In the case of S. meliloti-E. coli comparison, the identical curves for syntenic and non-syntenic genes with the majority at very low identity values (35%, Figure 5 panel c) reflect their greater phylogenetic distances.
The high identity of syntenic genes indirectly reveals their essential character; for the non-syntenic genes the low identity could represent adaptability to the ecological niche of the species. The identity relatedness in the syntenic genes among rhizobial species (Figure 6 panel a) revealed a cohesively evolved group; additionally the sequence differences were reflected in the theoretical pI plots of the proteins encoded by these genes, with a majority in the diagonal and the rest in sectors with strong pI deviation. Species signatures of the sequences (see Additional file 7) showed a differential level of changed residues and these could represent functional adaptation to a niche; this proposal is supported with the almost identical sequences of B. melitensis and B. suis. On the other hand, invariant peptides perhaps contribute to structural conformation . Experiments in progress in our lab will determine the validity of our proposal. Recently, a conservative change which altered the function of the transcriptional regulator BosR , and pathogenic differences in enteric bacteria due to the expression of PmrD regulators with divergent sequences  were reported.
A typical trait of Enterobacteriales is its pathogenic character, and this means a very intimate, frequent contact with their hosts and almost constant, homogeneous conditions: the extreme case is that of Buchnera sp., an aphid-obligate symbiont. In contrast, Rhizobiales are commonly found in the soil in saprophytic living style, and occasionally they associate with hosts, and therefore face more heterogeneous, variable environments. These features were reflected in the sequence comparisons of proteins from the arginine biosynthetic pathway (see Additional file 7). The syntenic gene organization was different in the Rhizobiales compared with Enterobacteriales, however it is important to note the high orthology and synteny degrees between members of each clade, the essentiality of functions covered by the syntenic genes and the high functional linkage in the microsyntenic regions in each chromosome (see below). From the previous observations we can obtain a general trend: bacterial clades present a particular chromosomal gene arrangement and such plasticity possibly was selected in relation to the niches occupied by these organisms.
As have others, some time ago we found a proteomic bimodal distribution when MW and pI were plotted for rhizobial gene products. The main fraction was located at acid and basic pI's, in a shape resembling butterfly's wings and the body, at neutral pI, presenting a low number of proteins. Recently Knight et al. , in a vast genomic study, reported similar results and additionally, related proteome similarities to shared metabolic features. In this respect, we plotted pI of syntenic products from pair comparisons of studied species in order to detect proteomic similarities in these organisms. It was possible to analyze the pI variability in each of the common syntenic products; however, it is necessary to carry out experimentation to determine the biological role of pI variation. The overrepresented membranal fraction in the syntenic products must be further explored to determine whether they respond to different extracellular/intracellular signals. A rapid evolution for the membranal proteomic fraction was suggested in the same study .
The high proportion of syntenic operons in the microsyntenic regions and the duplicated connections per gene in the network-forming regions (in regard to non-conserved regions), in Rhizobiales but also in Enterobacteriales, supported the functional linkage and interaction of these genes in the conserved tracts (see Additional file 14), and this is a factor which could help to define the role of selective pressure in maintaining the gene order.
Functional characterization of the predicted orthologs deepened our understanding of their cellular roles. In both Rhizobiales and Enterobacteriales, a clear essentiality was observed for chromosomal syntenic genes, agreeing with their sequence restrictions for change. Non-syntenic genes, on the other hand, appeared abundant in functions granting metabolic versatility to the cell. By calculating nonsynonymous/synonymous substitution rates, other authors have shown that in bacteria most conserved genes cover the essential functions of the cell .
Our synteny analysis defined a multi-level gene organization in the bacterial chromosome. Restriction of sequence variation in these genes, with clear essential functional roles, appeared extended to the conservation of chromosomal arrangement. In this way, synteny possibly has an important biological significance in these organisms.
Identification of orthologs
Available genome sequences of the fast growing Rhizobiales, S. meliloti 1021 [accession number GenBank: NC_003047], A. tumefaciens C58 (Cereon) [GenBank: NC_003062 and NC_003063], M. loti MAFF303099 [GenBank: NC_002678] and B. melitensis 16 M [GenBank: NC_003317 and NC_003318] were obtained from the Genome division of the NCBI Entrez system . Genomes of E. coli K12, S. typhimurium LT2 and Erwinia carotovora subsp. atroseptica SCRI1043 [accession numbers GenBank: NC_00913NC_003197, and NC_004547], were used for synteny and functional analysis. The genome of B. suis 1330 [accession numbers GenBank: NC_004310 and NC_004311 for chromosomes I and II, respectively] was used for synteny and sequence identity analysis. Genomes of Bradyrhizobium japonicum USDA110 [accession number GenBank: NC_004463] and Rhodopseudomonas palustris CGA009, also belonging to the alpha proteobacteria, were not considered in the main analysis because of their more distant phylogenetic relationship with the fast-growing rhizobia. To obtain a comprehensive view of shared genes among Rhizobiales, we differentiated genes by orthology and their presence in the same or different replicons in the analyzed species. Chromosomal orthologs were assigned by the best bidirectional hit between pairs of organisms, using the Fasta34 program . Unidirectional best hits (homologs) were considered to cover complete chromosomal gene number (see Results) and for detection of chromosomal genes with plasmidic homologs. Parameters were: an identity of at least 50%, overlapping by at least in 150 nt and an expectance (E) score of <10-3. The base organism for comparisons was S. meliloti 1021. GenBank accession numbers of proteins used for alignments in the order ArgB, ArgC, ArgD, ArgF, ArgG and ArgH, were as follows. Rhizobiales. R. palustris: NP_945982.1, NP_947833.1, NP_950107.1, NP_950106.1, NP_945745.1, NP_950077.1; B. melitensis: NP_541250.1, NP_540088.1, NP_540538.1, NP_540537.1, NP_540787.1, NP_539004.1; M. loti: NP_105609.1, NP_108547.1, NP_106269.1, NP_106270.1, NP_105253.1, NP_104594.1; A. tumefaciens: NP_353412.1, NP_354256.1, NP_353456.1, NP_353457.1, NP_355604.1, NP_357013.1; S. meliloti: NP_384545.1, NP_385346.1, NP_384623.1, NP_384624.1, NP_387315.1, NP_386753.1. Enterobacteriales. Buchnera sp.: NP_239886.1, NP_239885.1, NP_240341.1, NP_240186.1, NP_239887.1, NP_239888.1; E. coli: NP_418394.1, NP_418393.1, NP_417818.1, NP_414807.1, NP_417640.1, NP_418395.1; E. carotovora: YP_048320.1, YP_048319.1, YP_052152.1, YP_048510.1, YP_048232.1, YP_048321.1S. typhimurium: NP_457937.1, NP_457938.1, NP_458434.1, NP_458882.1, NP_457671.1, NP_457936.1; S. flexneri: NP_838925.1, NP_838926.1, NP_839526.1, NP_839629.1, NP_838682.1, NP_838924.1. All data sets are available on request.
Procedures for detection of syntenic genes, microsyntenic region formation and operon similarity
To consider chromosomal orthologs as syntenic among pairs of organisms, at least two genes must remain contiguous in both chromosomes. Microsyntenic region formation and extension fulfilled the following criterion: a pair of predicted orthologs separated from at least one other by no more than three genes (from the rest of categories). The minimal region was formed by a stretch containing three syntenic genes. Operon prediction was performed as reported by Moreno-Hagelsieb and Collado-Vides . Rearrangements were graphed using initial positions of the microsyntenic regions in each chromosome. The syntenic regions of M. loti and B. melitensis I chromosomes were graphed so as to increase the colinearity in these replicons. The M. loti chromosome was segmented into two halves at 3.5 Mb position and the fragment covering from 3.5 to 7.0 Mb was located in the first position and then both halves were aligned with microsyntenic regions of S. meliloti. In the case of Brucella chromosome I, the origin was inverted. Graphs were obtained using the GenVision program (DNAStar Inc., Madison, WI). For operon similarity calculation, a limit of three different genes in each operon was allowed. All data sets are available on request.
Detection of horizontal gene transfer
Prediction of horizontally transferred genes in the Rhizobiales genomes was performed using the method described by Medrano-Soto et al. . Briefly, it on based in similar gene length, maximum global protein identity, conflicting phylogenies and codon usage of xenologous genes. In the case of A. tumefaciens, this rate was calculated using sequences and annotation obtained from the U. of Washington Sequencing Project . Species-specific (or orphan) genes were not considered by two reasons: (1) the lack of orthologs in other genomes precluded phylogenetic analysis, and (2) impossibility of correlating these genes to a synteny.
Sequence comparison of chromosomal orthologs in Rhizobiales
The identity of peptidic sequences of chromosomal predicted orthologs were used to graph the distribution curves. The asymmetry of distribution curves, or skewness, was calculated by the asymmetry coefficient of Pearson (g1) as described elsewhere . To correlate nucleotide sequence identities of the chromosomal predicted orthologs in the four rhizobial genomes, the gene identity of the predicted orthologs from the S. meliloti-M. loti comparison was graphed in progressive order. Then, corresponding predicted orthologs of the other comparisons were located at their corresponding identity percentages. Correlation coefficient values were calculated by the Pearson method. Plasmidic homologs and species-specific genes were not graphed because they have no counterparts in the pairs of analyzed genomes.
Theoretical proteome and transmembranal protein prediction
Theoretical proteomes were obtained by calculating molecular weight and isoelectric point for each translated chromosomal predicted ortholog. Both parameters were estimated with the pI/MW prediction tool of the Laquip Proteomic Team page . To determine the set of orthologs coding proteins predicted to interact with the cell membrane, we used the TMAP program, version 46 , available in the EMBOSS package (European Molecular Biology Open Software Suite, ). Correlation coefficient values were obtained with the Pearson method. Alignments were performed with ClustalW .
Functional classification of chromosomal predicted orthologs
Chromosomal predicted orthologs were assigned to the functional classes used for Agrobacterium tumefaciens C58 in the U. of Washington genome report . For the distribution of functions coded in the microsyntenic regions along the S. meliloti chromosome, classes were assigned into Operational (including amino acid, fatty acid, carbohydrate and nucleotide metabolism, energy generation, central intermediary metabolism, transport and cofactor synthesis), Informational (DNA metabolism, transcription, translation and regulatory functions), and Cellular processes (cell envelope, cell division, secretion and chemotaxis) superclasses. This grouping, except for Cellular processes, is similar to that of Rivera et al. . For functional relationship inference, the ProLinks  and STRING databases  were used with permission. The confidence levels used were 0.6 and 0.9, respectively. Resulting networks with ProLinks with less than 6 links were omitted from the count. Networks were constructed with the Pajek Program (written by A. Vlado), version 1.02, available in the web . Assignment of metabolic pathways was performed using the MetaCyc database , with permission.
Genomes Online Database. [http://www.genomesonline.org]
Itoh T, Martin W, Nei M: Acceleration of genomic evolution caused by enhanced mutation rate in endocellular symbionts. Proc Natl Acad Sci U S A. 2002, 99 (20): 12944-12948. 10.1073/pnas.192449699.
Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Microevolutionary Genomics of Bacteria. Theoretical Population Biology. 2002, 61 (4): 435-10.1006/tpbi.2002.1588.
Rocha EP: Order and disorder in bacterial genomes. Curr Opin Microbiol. 2004, 7 (5): 519-527. 10.1016/j.mib.2004.08.006.
Korbel JO, Jensen LJ, von Mering C, Bork P: Analysis of genomic context: prediction of functional associations from conserved bidirectionally transcribed gene pairs. Nat Biotechnol. 2004, 22 (7): 911-917. 10.1038/nbt988.
Tamames J: Evolution of gene order conservation in prokaryotes. Genome Biol. 2001, 2 (6): R0020-10.1186/gb-2001-2-6-research0020.
Zheng XH, Lu F, Wang ZY, Zhong F, Hoover J, Mural R: Using shared genomic synteny and shared protein functions to enhance the identification of orthologous gene pairs. Bioinformatics. 2005, 21 (6): 703-710. 10.1093/bioinformatics/bti045.
Nakamura Y, Nishio Y, Ikeo K, Gojobori T: The genome stability in Corynebacterium species due to lack of the recombinational repair system. Gene. 2003, 317 (1-2): 149-155. 10.1016/S0378-1119(03)00653-X.
Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H: Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature. 2000, 407 (6800): 81-86. 10.1038/35024074.
Cannon SB, McCombie WR, Sato S, Tabata S, Denny R, Palmer L, Katari M, Young ND, Stacey G: Evolution and microsynteny of the apyrase gene family in three legume genomes. Mol Genet Genomics. 2003, 270 (4): 347-361. 10.1007/s00438-003-0928-x.
Gualtieri G, Bisseling T: Microsynteny between the Medicago truncatula SYM2-orthologous genomic region and another region located on the same chromosome arm. Theor Appl Genet. 2002, 105 (5): 771-779. 10.1007/s00122-002-0949-6.
Bankier AT, Spriggs HF, Fartmann B, Konfortov BA, Madera M, Vogel C, Teichmann SA, Ivens A, Dear PH: Integrated mapping, chromosomal sequencing and sequence analysis of Cryptosporidium parvum. Genome Res. 2003, 13 (8): 1787-1799.
Bannantine JP, Zhang Q, Li LL, Kapur V: Genomic homogeneity between Mycobacterium avium subsp. avium and Mycobacterium avium subsp. paratuberculosis belies their divergent growth rates. BMC Microbiol. 2003, 3 (1): 10-10.1186/1471-2180-3-10.
Buchrieser C, Rusniok C, Kunst F, Cossart P, Glaser P: Comparison of the genome sequences of Listeria monocytogenes and Listeria innocua: clues for evolution and pathogenicity. FEMS Immunol Med Microbiol. 2003, 35 (3): 207-213. 10.1016/S0928-8244(02)00448-0.
Eppinger M, Baar C, Raddatz G, Huson DH, Schuster SC: Comparative analysis of four Campylobacterales. Nat Rev Microbiol. 2004, 2 (11): 872-885. 10.1038/nrmicro1024.
Horimoto K, Fukuchi S, Mori K: Comprehensive comparison between locations of orthologous genes on archaeal and bacterial genomes. Bioinformatics. 2001, 17 (9): 791-802. 10.1093/bioinformatics/17.9.791.
Moreira LM, de Souza RF, Almeida NFJ, Setubal JC, Oliveira JC, Furlan LR, Ferro JA, da Silva AC: Comparative genomics analyses of citrus-associated bacteria. Annu Rev Phytopathol. 2004, 42: 163-184. 10.1146/annurev.phyto.42.040803.140310.
Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, White O, Hickey EK, Peterson J, Utterback T, Berry K, Bass S, Linher K, Weidman J, Khouri H, Craven B, Bowman C, Dodson R, Gwinn M, Nelson W, DeBoy R, Kolonay J, McClarty G, Salzberg SL, Eisen J, Fraser CM: Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 2000, 28 (6): 1397-1406. 10.1093/nar/28.6.1397.
Krawiec S, Riley M: Organization of the bacterial chromosome. Microbiol Rev. 1990, 54 (4): 502-539.
Paulsen IT, Seshadri R, Nelson KE, Eisen JA, Heidelberg JF, Read TD, Dodson RJ, Umayam L, Brinkac LM, Beanan MJ, Daugherty SC, Deboy RT, Durkin AS, Kolonay JF, Madupu R, Nelson WC, Ayodeji B, Kraul M, Shetty J, Malek J, Van Aken SE, Riedmuller S, Tettelin H, Gill SR, White O, Salzberg SL, Hoover DL, Lindler LE, Halling SM, Boyle SM, Fraser CM: The Brucella suis genome reveals fundamental similarities between animal and plant pathogens and symbionts. Proc Natl Acad Sci U S A. 2002, 99 (20): 13148-13153. 10.1073/pnas.192319099.
Galibert F, Finan TM, Long SR, Puhler A, Abola P, Ampe F, Barloy-Hubler F, Barnett MJ, Becker A, Boistard P, Bothe G, Boutry M, Bowser L, Buhrmester J, Cadieu E, Capela D, Chain P, Cowie A, Davis RW, Dreano S, Federspiel NA, Fisher RF, Gloux S, Godrie T, Goffeau A, Golding B, Gouzy J, Gurjal M, Hernandez-Lucas I, Hong A, Huizar L, Hyman RW, Jones T, Kahn D, Kahn ML, Kalman S, Keating DH, Kiss E, Komp C, Lelaure V, Masuy D, Palm C, Peck MC, Pohl TM, Portetelle D, Purnelle B, Ramsperger U, Surzycki R, Thebault P, Vandenbol M, Vorholter FJ, Weidner S, Wells DH, Wong K, Yeh KC, Batut J: The composite genome of the legume symbiont Sinorhizobium meliloti. Science. 2001, 293 (5530): 668-672.
Goodner B, Hinkle G, Gattung S, Miller N, Blanchard M, Qurollo B, Goldman BS, Cao Y, Askenazi M, Halling C, Mullin L, Houmiel K, Gordon J, Vaudin M, Iartchouk O, Epp A, Liu F, Wollam C, Allinger M, Doughty D, Scott C, Lappas C, Markelz B, Flanagan C, Crowell C, Gurson J, Lomo C, Sear C, Strub G, Cielo C, Slater S: Genome Sequence of the Plant Pathogen and Biotechnology Agent Agrobacterium tumefaciens C58. Science. 2001, 294 (5550): 2323-2328. 10.1126/science.1066803.
Doolittle WF: Lateral genomics. Trends Cell Biol. 1999, 9 (12): M5-8. 10.1016/S0962-8924(99)01664-5.
Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405 (6784): 299-304. 10.1038/35012500.
Koski LB, Morton RA, Golding GB: Codon bias and base composition are poor indicators of horizontally transferred genes. Mol Biol Evol. 2001, 18 (3): 404-412.
Wang B: Limitations of compositional approach to identifying horizontally transferred genes. J Mol Evol. 2001, 53 (3): 244-250. 10.1007/s002390010214.
Medrano-Soto A, Moreno-Hagelsieb G, Vinuesa P, Christen JA, Collado-Vides J: Successful lateral transfer requires codon usage compatibility between foreign genes and recipient genomes. Mol Biol Evol. 2004, 21 (10): 1884-1894. 10.1093/molbev/msh202.
Hao W, Golding GB: Patterns of bacterial gene movement. Mol Biol Evol. 2004, 21 (7): 1294-1307. 10.1093/molbev/msh129.
Jain R, Rivera MC, Moore JE, Lake JA: Horizontal gene transfer accelerates genome innovation and evolution. Mol Biol Evol. 2003, 20 (10): 1598-1602. 10.1093/molbev/msg154.
Pellegrini M, Marcotte EM, Thompson MJ, Eisenberg D, Yeates TO: Assigning protein functions by comparative genome analysis: protein phylogenetic profiles. Proc Natl Acad Sci U S A. 1999, 96 (8): 4285-4288. 10.1073/pnas.96.8.4285.
Overbeek R, Fonstein M, D'Souza M, Pusch GD, Maltsev N: The use of gene clusters to infer functional coupling. Proc Natl Acad Sci U S A. 1999, 96 (6): 2896-2901. 10.1073/pnas.96.6.2896.
Marcotte EM, Pellegrini M, Ho-Leung N, Rice D, Yeates T, Elsenberg D: Detecting protein function and protein-protein interactions from genome sequences. Science. 1999, 285 (5428): 751-753. 10.1126/science.285.5428.751.
Yanai I, DeLisi C: The society of genes: networks of functional links between genes from comparative genomics. Genome Biol. 2002, 3 (11): R0064-10.1186/gb-2002-3-11-research0064.
Bowers PM, Pellegrini M, Thompson MJ, Fierro J, Yeates TO, Eisenberg D: Prolinks: a database of protein functional linkages derived from coevolution. Genome Biol. 2004, 5 (5): R35-10.1186/gb-2004-5-5-r35.
von Mering C, Zdobnov EM, Tsoka S, Ciccarelli FD, Pereira-Leal JB, Ouzounis CA, Bork P: Genome evolution reveals biochemical networks and functional modules. Proc Natl Acad Sci U S A. 2003, 100 (26): 15428-15433. 10.1073/pnas.2136809100.
Freiberg C, Fellay R, Bairoch A, Broughton WJ, Rosenthal A, Perret X: Molecular basis of symbiosis between Rhizobium and legumes. Nature. 1997, 387 (6631): 394-401. 10.1038/387394a0.
Gonzalez V, Bustos P, Ramirez-Romero MA, Medrano-Soto A, Salgado H, Hernandez-Gonzalez I, Hernandez-Celis JC, Quintero V, Moreno-Hagelsieb G, Girard L, Rodriguez O, Flores M, Cevallos MA, Collado-Vides J, Romero D, Davila G: The mosaic structure of the symbiotic plasmid of Rhizobium etli CFN42 and its relation to other symbiotic genome compartments. Genome Biol. 2003, 4 (6): R36-10.1186/gb-2003-4-6-r36.
Gottfert M, Rothlisberger S, Kundig C, Beck C, Marty R, Hennecke H: Potential symbiosis-specific genes uncovered by sequencing a 410-kilobase DNA region of the Bradyrhizobium japonicum chromosome. J Bacteriol. 2001, 183 (4): 1405-1412. 10.1128/JB.183.4.1405-1412.2001.
Hirsch AM, Drake D, Jacobs TW, Long SR: Nodules are induced on alfalfa roots by Agrobacterium tumefaciens and Rhizobium trifolii containing small segments of the Rhizobium meliloti nodulation region. J Bacteriol. 1985, 161 (1): 223-230.
Martinez E, Palacios R, Sanchez F: Nitrogen-fixing nodules induced by Agrobacterium tumefaciens harboring Rhizobium phaseoli plasmids. J Bacteriol. 1987, 169 (6): 2828-2834.
Weller SA, Stead DE, Young JP: Acquisition of an Agrobacterium Ri plasmid and pathogenicity by other alpha-Proteobacteria in cucumber and tomato crops affected by root mat. Appl Environ Microbiol. 2004, 70 (5): 2779-2785. 10.1128/AEM.70.5.2779-2785.2004.
Moran NA, Wernegreen JJ: Lifestyle evolution in symbiotic bacteria: insights from genomics. Trends In Ecology And Evolution. 2000, 15 (8): 321-326. 10.1016/S0169-5347(00)01902-9.
Kaneko T, Nakamura Y, Sato S, Asamizu E, Kato T, Sasamoto S, Watanabe A, Idesawa K, Ishikawa A, Kawashima K, Kimura T, Kishida Y, Kiyokawa C, Kohara M, Matsumoto M, Matsuno A, Mochizuki Y, Nakayama S, Nakazaki N, Shimpo S, Sugimoto M, Takeuchi C, Yamada M, Tabata S: Complete genome structure of the nitrogen-fixing symbiotic bacterium Mesorhizobium loti. DNA Res. 2000, 7 (6): 331-338. 10.1093/dnares/7.6.331.
Kaneko T, Nakamura Y, Sato S, Minamisawa K, Uchiumi T, Sasamoto S, Watanabe A, Idesawa K, Iriguchi M, Kawashima K, Kohara M, Matsumoto M, Shimpo S, Tsuruoka H, Wada T, Yamada M, Tabata S: Complete genomic sequence of nitrogen-fixing symbiotic bacterium Bradyrhizobium japonicum USDA110. DNA Res. 2002, 9 (6): 189-197. 10.1093/dnares/9.6.189.
Wood DW, Setubal JC, Kaul R, Monks DE, Kitajima JP, Okura VK, Zhou Y, Chen L, Wood GE, Almeida NFJ, Woo L, Chen Y, Paulsen IT, Eisen JA, Karp PD, Bovee DS, Chapman P, Clendenning J, Deatherage G, Gillet W, Grant C, Kutyavin T, Levy R, Li MJ, McClelland E, Palmieri A, Raymond C, Rouse G, Saenphimmachak C, Wu Z, Romero P, Gordon D, Zhang S, Yoo H, Tao Y, Biddle P, Jung M, Krespan W, Perry M, Gordon-Kamm B, Liao L, Kim S, Hendrick C, Zhao ZY, Dolan M, Chumley F, Tingey SV, Tomb JF, Gordon MP, Olson MV, Nester EW: The genome of the natural genetic engineer Agrobacterium tumefaciens C58. Science. 2001, 294 (5550): 2317-2323. 10.1126/science.1066804.
DelVecchio VG, Kapatral V, Redkar RJ, Patra G, Mujer C, Los T, Ivanova N, Anderson I, Bhattacharyya A, Lykidis A, Reznik G, Jablonski L, Larsen N, D'Souza M, Bernal A, Mazur M, Goltsman E, Selkov E, Elzer PH, Hagius S, O'Callaghan D, Letesson JJ, Haselkorn R, Kyrpides N, Overbeek R: The genome sequence of the facultative intracellular pathogen Brucella melitensis. Proc Natl Acad Sci U S A. 2002, 99 (1): 443-448. 10.1073/pnas.221575398.
Larimer FW, Chain P, Hauser L, Lamerdin J, Malfatti S, Do L, Land ML, Pelletier DA, Beatty JT, Lang AS, Tabita FR, Gibson JL, Hanson TE, Bobst C, Torres JL, Peres C, Harrison FH, Gibson J, Harwood CS: Complete genome sequence of the metabolically versatile photosynthetic bacterium Rhodopseudomonas palustris. Nat Biotechnol. 2004, 22 (1): 55-61. 10.1038/nbt923.
COGs: Phylogenetic classification of protein encoded in complete genomes database. [http://www.ncbi.nlm.nih.gov/COG]
Gerstein M: Patterns of protein-fold usage in eight microbial genomes: a comprehensive structural census. Proteins. 1998, 33 (4): 518-534. 10.1002/(SICI)1097-0134(19981201)33:4<518::AID-PROT5>3.0.CO;2-J.
Gerstein M, Levitt M: A structural census of the current population of protein sequences. Proc Natl Acad Sci U S A. 1997, 94 (22): 11911-11916. 10.1073/pnas.94.22.11911.
Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 305 (3): 567-580. 10.1006/jmbi.2000.4315.
Fernandez-Herrero LA, Badet-Denisot MA, Badet B, Berenguer J: glmS of Thermus thermophilus HB8: an essential gene for cell-wall synthesis identified immediately upstream of the S-layer gene. Mol Microbiol. 1995, 17 (1): 1-12.
Mac Siomoin RA, Nakata N, Murai T, Yoshikawa M, Tsuji H, Sasakawa C: Identification and characterization of ispA, a Shigella flexneri chromosomal gene essential for normal in vivo cell division and intracellular spreading. Mol Microbiol. 1996, 19 (3): 599-609. 10.1046/j.1365-2958.1996.405941.x.
Maggio-Hall LA, Claas KR, Escalante-Semerena JC: The last step in coenzyme B(12) synthesis is localized to the cell membrane in bacteria and archaea. Microbiology. 2004, 150 (Pt 5): 1385-1395. 10.1099/mic.0.26952-0.
Wu G, Williams HD, Gibson F, Poole RK: Mutants of Escherichia coli affected in respiration: the cloning and nucleotide sequence of ubiA, encoding the membrane-bound p-hydroxybenzoate:octaprenyltransferase. J Gen Microbiol. 1993, 139 (8): 1795-1805.
Xia T, Song J, Zhao G, Aldrich H, Jensen RA: The aroQ-encoded monofunctional chorismate mutase (CM-F) protein is a periplasmic enzyme in Erwinia herbicola. J Bacteriol. 1993, 175 (15): 4729-4737.
Kosak ST, Groudine M: Gene order and dynamic domains. Science. 2004, 306 (5696): 644-647. 10.1126/science.1103864.
Boussau B, Karlberg EO, Frank AC, Legault BA, Andersson SG: Computational inference of scenarios for alpha-proteobacterial genome evolution. Proc Natl Acad Sci U S A. 2004, 101 (26): 9722-9727. 10.1073/pnas.0400975101.
Sullivan JT, Ronson CW: Evolution of rhizobia by acquisition of a 500-kb symbiosis island that integrates into a phe-tRNA gene. Proc Natl Acad Sci U S A. 1998, 95 (9): 5145-5149. 10.1073/pnas.95.9.5145.
Wong K, Golding GB: A phylogenetic analysis of the pSymB replicon from the Sinorhizobium meliloti genome reveals a complex evolutionary history. Can J Microbiol. 2003, 49 (4): 269-280. 10.1139/w03-037.
Prakash T, Ramakrishnan C, Dash D, Brahmachari SK: Conformational analysis of invariant peptide sequences in bacterial genomes. J Mol Biol. 2005, 345 (5): 937-955. 10.1016/j.jmb.2004.11.008.
Seshu J, Boylan JA, Hyde JA, Swingle KL, Gherardini FC, Skare JT: A conservative amino acid change alters the function of BosR, the redox regulator of Borrelia burgdorferi. Mol Microbiol. 2004, 54 (5): 1352-1363. 10.1111/j.1365-2958.2004.04352.x.
Winfield MD, Groisman EA: Phenotypic differences between Salmonella and Escherichia coli resulting from the disparate regulation of homologous genes. Proc Natl Acad Sci U S A. 2004, 101 (49): 17162-17167. 10.1073/pnas.0406038101.
Knight CG, Kassen R, Hebestreit H, Rainey PB: Global analysis of predicted proteomes: functional adaptation of physical properties. Proc Natl Acad Sci U S A. 2004, 101 (22): 8390-8395. 10.1073/pnas.0307270101.
Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Essential genes are more evolutionarily conserved than are nonessential genes in bacteria. Genome Res. 2002, 12 (6): 962-968. 10.1101/gr.87702. Article published online before print in May 2002.
Completed Microbial Genomes. [http://www.ncbi.nlm.nih.gov/genomes/MICROBES/Complete.html]
Persson B, Argos P: Prediction of transmembrane segments in proteins utilising multiple sequence alignments. J Mol Biol. 1994, 237 (2): 182-192. 10.1006/jmbi.1994.1220.
Moreno-Hagelsieb G, Collado-Vides J: A powerful non-homology method for the prediction of operons in prokaryotes. Bioinformatics. 2002, 18 Suppl 1: S329-36.
Sokal RR, Rohlf FJ: Biometry: The pinciples and practice of statistics in biological research. Edited by: 3rd . 2003, New York , WH Freeman and Company, 109-117.
Laquip Proteomic Team. [http://proteome.ibi.unicamp.br/tools]
Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Natl Acad Sci U S A. 1988, 85 (8): 2444-2448.
Software EMBO: EMBOSS. [http://www.emboss.org]
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680.
Rivera MC, Jain R, Moore JE, Lake JA: Genomic evidence for two functionally distinct gene classes. Proc Natl Acad Sci U S A. 1998, 95 (11): 6239-6244. 10.1073/pnas.95.11.6239.
Krieger CJ, Zhang P, Mueller LA, Wang A, Paley S, Arnaud M, Pick J, Rhee SY, Karp PD: MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res. 2004, 32 (Database issue): D438-42. 10.1093/nar/gkh100.
This work was partially supported by the grant 028 from the Consejo Nacional de Ciencia y Tecnología-México.
We wish thank to Gabriel Moreno-Hagelsieb for providing the operon prediction, César Bonavides for help with metabolic pathway assignment, Rafael Palacios, Michael Dunn and Yolanda Mora for manuscript review, and an anonymous referee for valuable suggestions.
GG and AA performed the computer predictions. HP made the functional analysis, interpreted the data and wrote the paper. RD and MAV participated in the work design. AM-S made the HGT and CRI analysis. JM conceived and directed the project. All authors read and approved the final manuscript.
Gabriela Guerrero, Humberto Peralta contributed equally to this work.
Electronic supplementary material
Additional File 1: Schematic representation of the Rhizobiales chromosomes in comparison with S. meliloti, according to the gene classification of predicted orthologs and homologs. Panels: (a), S. meliloti-A. tumefaciens comparison. (b), S. meliloti-M. loti comparison. (c), S. meliloti-B. melitensis comparison. Red striped bars, syntenic genes with the organism in comparison. Blue striped bars, non-syntenic genes with the organism in comparison. White bars, homologs with other Rhizobiales chromosomes (for S. meliloti, compare with Fig. 1, white fraction). Green bars, homologs in plasmids. Gray bars, species-specific genes. In panels a and c, the S. meliloti chromosome shows syntenic and non-syntenic genes with both replicons of the organisms under comparison. Red striped bars, syntenic genes, lower fraction: with (a) At-C and (c) Bm-I chromosomes; upper fraction: with (a) At-L and (c) Bm-II chromosomes. Blue striped bars, non-syntenic genes, lower fraction: with (a) At-C and (c) Bm-I chromosomes; upper fraction: with (a) At-L and (c) Bm-II chromosomes. (JPEG 160 KB)
Additional File 2: Synteny histogram of S. meliloti in comparison with M. loti chromosome. Red bars, syntenic genes. Surrounded with yellow boxes, microsyntenic regions. Microsyntenic regions are denoted by letters (and numbers) in progressive order. Blue bars, non-syntenic genes. Green bars, homologs in plasmids. Gray bars, species-specific genes. White bars, homologs with other Rhizobiales chromosomes. Direction of transcription is denoted by upper (plus) or lower (minus) positions in respect to central line. Predicted operons are denoted by red arrows. Scale in bp. (EPS 648 KB)
Additional File 3: Synteny histogram of S. meliloti in comparison with B. melitensis chromosomes. Red bars, syntenic genes. Surrounded with yellow boxes, microsyntenic regions with B. melitensis chromosome I (Bm-I). Surrounded with light green boxes, microsyntenic regions with B. melitensis chromosome II (Bm-II). Microsyntenic regions are denoted by letters (and numbers) in progressive order. Dark blue bars, non-syntenic genes with Bm-I. Light blue bars, nonsyntenic genes with Bm-II. Green bars, homologs in plasmids. Gray bars, species-specific genes. White bars, homologs with other Rhizobiales chromosomes. Direction of transcription is denoted by upper (plus) or lower (minus) positions in respect to central line. Predicted operons are denoted by red arrows. Scale in bp. (EPS 632 KB)
Additional File 4: Synteny of Rhizobiales and Enterobacteriales. Panels: (a), S. meliloti-A. tumefaciens circular chromosomes comparison. (b), S. meliloti-M. loti comparison. (c), S. meliloti-B. melitensis chromosome I comparison. (d), S. meliloti-E. coli comparison. (e), B. suis-B. melitensis chromosomes I comparison. Red dots, syntenic genes. Blue dots, non-syntenic genes. Scales in bp. (PPT 978 KB)
Additional File 5: Synteny histogram of E. coli in comparison with E. carotovora chromosome. Red bars, syntenic genes. Surrounded with yellow boxes, microsyntenic regions with E. carotovora chromosome Microsyntenic regions are denoted by letters (and numbers) in progressive order. Blue bars, non-syntenic genes. Gray bars, non-orthologous genes. Scale in bp. (EPS 728 KB)
Additional File 6: Sequence identity distribution of chromosomal translated orthologs. Panels: (a), syntenic and non-syntenic products from the S. meliloti-M. loti comparison. (b), syntenic and non-syntenic products from the S. meliloti-B. melitensis (chromosomes I and II) comparison. Y-axis, relative proportions. Red bars, syntenic genes. Blue bars, non-syntenic genes. Y-axis, relative proportions. (PPT 53 KB)
Additional File 7: Sequence alignments and data from the alignments of proteins from the arginine biosynthetic pathway in Rhizobiales and Enterobacteriales. Panels: (a), ArgC in Rhizobiales. Identical residues for aech position are marked with yellow. Least abundant residues for a given position are denoted with an specific color for each of the species: dark blue, differences in R. palustris; green, differences in B. melitensis; red, differences in M. loti; gray, differences in A. tumefaciens; violet, S. meliloti. (b), ArgC in Rhizobiales and Enterobacteriales. Identical residues for aech position are marked with yellow. Least abundant residues for a given position are denoted with an specific color for each of the species: Rhizobiales, same code of panel (a). Enterobacteriales: brown, Buchnera; pink, E. carotovora; blue, S. typhimurium. E. coli and S. flexneri, none. (c), data of the identity (*) and similarity (:*) in residues and in percentage (bold) of the alignments of the proteins in Rhizobiales and Enterobacteriales. (PDF 157 KB)
Additional File 8: Theoretical isoelectric points (pI) of syntenic products. Panels: (a), S. meliloti-M. loti comparison. (b), S. meliloti-B. melitensis (chromosomes I and II) comparison. Dots represent translated products. Scales in pH units. (PPT 607 KB)
Additional File 9: Summary of proteins with differential pI's from comparisons with chromosomes of S. meliloti, A. tumefaciens (circular) and M. loti. (DOC 37 KB)
Additional File 10: Theoretical isoelectric points (pI) of nonsyntenic products. Panels: (a), S. meliloti-A. tumefaciens (both chromosomes) comparison. (b), S. meliloti-M. loti comparison. (c), S. meliloti-B. melitensis (both chromosomes) comparison. Dots represent translated products. Scales in pH units. (PPT 297 KB)
Additional File 11: Functional categories of syntenic products of the membranal prediction of S. meliloti- A. tumefaciens (circular chromosome) comparison. (DOC 28 KB)
Additional File 12: Metabolic pathways covered by syntenic products of the S. meliloti-A. tumefaciens comparison. Scheme belongs to the MetaCyc pathways for S. meliloti chromosome, used with permission. Highlighted with green, reactions covered with syntenic products. (EPS 388 KB)
Additional File 13: Coverage of functional classes with syntenic and non-syntenic genes. Panels: (a), S. meliloti-M. loti comparison. Classes: 1) Translation, 2) Transcription, 3) Purine, pyrimidine, nucleoside and nucleotide metabolism, 4) Cellular processes, 5) Energy metabolism, 6) Cell envelope, 7) Fatty acid and phospholipid metabolism, 8) Biosynthesis of cofactors, prosthetic groups and carriers, 9) Transport and ATP binding proteins, 10) Amino acid metabolism, 11) DNA metabolism, 12) Regulatory functions, 13) Central intermediary metabolism. (b), S. meliloti-B. melitensis comparison. Classes: 1) Transcription, 2) Translation, 3) Cellular processes, 4) Biosynthesis of cofactors, prosthetic groups and carriers, 5) Cell envelope, 6) Energy metabolism, 7) Fatty acid and phospholipid metabolism, 8) Purine, pyrimidine, nucleoside and nucleotide metabolism, 9) Amino acid metabolism, 10) Transport and ATP binding proteins, 11) DNA metabolism, 12) Regulatory functions, 13) Central intermediary metabolism. (c), E. coli-E. carotovora comparison. Classes: 1) Biosynthesis of cofactors, prosthetic groups and carriers, 2) Purine, pyrimidine, nucleoside and nucleotide metabolism, 3) Translation, 4) Fatty acid and phospholipid metabolism, 5) Transcription, 6) Cellular processes, 7) DNA metabolism, 8) Energy metabolism, 9) Amino acid metabolism. 10) Cell envelope, 11) Regulatory functions, 12) Transport and ATP binding proteins, 13) Central intermediary metabolism. Note that order of classes is different to that in Fig. 8. Red bars, syntenic genes. Blue bars, non-syntenic genes. (PPT 87 KB)
Additional File 14: Connectivity values from networks formed by microsyntenic and non-conserved regions in (a) S. meliloti (in comparison with A. tumefaciens) and (b) E. coli (in comparison with E. carotovora). Y-axis, connections per network. First syntenic networks, with 1060 (S. meliloti) and 810 (E. coli) connections, were omitted for clarity. Arranged in decrecent connectivity order. Gray bars, microsyntenic regions. Black bars, non-conserved regions. Successive networks, with connectivity values lower than 6, were omitted. (PPT 96 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Guerrero, G., Peralta, H., Aguilar, A. et al. Evolutionary, structural and functional relationships revealed by comparative analysis of syntenic genes in Rhizobiales. BMC Evol Biol 5, 55 (2005). https://doi.org/10.1186/1471-2148-5-55
- Horizontal Gene Transfer
- Horizontal Gene Transfer Event
- Brucella Melitensis
- Syntenic Gene
- Gamma Proteobacteria