Ribosomal protein gene cluster analysis in eubacterium genomics: homology between Sinorhizobium meliloti strain 1021 and Bacillus subtilis

Was the Last Bacterial Common Ancestor a Monoderm after All?

The very nature of the last bacterial common ancestor (LBCA), in particular the characteristics of its cell wall, is a critical issue to understand the evolution of life on earth. Although knowledge of the relationships between bacterial phyla has made progress with the advent of phylogenomics, many questions remain, including on the appearance or disappearance of the outer membrane of diderm bacteria (also called Gram-negative bacteria). The phylogenetic transition between monoderm (Gram-positive bacteria) and diderm bacteria, and the associated peptidoglycan expansion or reduction, requires clarification. Herein, using a phylogenomic tree of cultivated and characterized bacteria as an evolutionary framework and a literature review of their cell-wall characteristics, we used Bayesian ancestral state reconstruction to infer the cell-wall architecture of the LBCA. With the same phylogenomic tree, we further revisited the evolution of the division and cell-wall synthesis (dcw) gene cluster using homology- and model-based methods. Finally, extensive similarity searches were carried out to determine the phylogenetic distribution of the genes involved with the biosynthesis of the outer membrane in diderm bacteria. Quite unexpectedly, our analyses suggest that all cultivated and characterized bacteria might have evolved from a common ancestor with a monoderm cell-wall architecture. If true, this would indicate that the appearance of the outer membrane was not a unique event and that selective forces have led to the repeated adoption of such an architecture. Due to the lack of phenotypic information, our methodology cannot be applied to all extant bacteria. Consequently, our conclusion might change once enough information is made available to allow the use of an even more diverse organism selection.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Operon Concatenation Is an Ancient Feature That Restricts the Potential to Rearrange Bacterial Chromosomes

The last common ancestor of the Gammaproteobacteria carried an important 40-kb chromosome section encoding 51 proteins of the transcriptional and translational machinery. These genes were organized into eight contiguous operons (rrnB-tufB-secE-rpoBC-str-S10-spc-alpha). Over 2 Gy of evolution, in different lineages, some of the operons became separated by multigene insertions. Surprisingly, in many Enterobacteriaceae, much of the ancient organization is conserved, indicating a strong selective force on the operons to remain colinear. Here, we show for one operon pair, tufB-secE in Salmonella, that an interruption of contiguity significantly reduces growth rate. Our data show that the tufB-secE operons are concatenated by an interoperon terminator–promoter overlap that plays a significant role regulating gene expression. Interrupting operon contiguity interferes with this regulation, reducing cellular fitness. Six operons of the ancestral chromosome section remain contiguous in Salmonella (tufB-secE-rpoBC and S10-spc-alpha) and, strikingly, each of these operon pairs is also connected by an interoperon terminator–promoter overlap. Accordingly, we propose that operon concatenation is an ancient feature that restricts the potential to rearrange bacterial chromosomes and can select for the maintenance of a colinear operon organization over billions of years.

Taxonomic assignment of uncultured prokaryotes with long range PCR targeting the spectinomycin operon

The taxonomic assignment of uncultured prokaryotes to known taxa is a major challenge in microbial systematics. This relies usually on the phylogenetic analysis of the ribosomal small subunit RNA or a few housekeeping genes. Recent works have disclosed ribosomal proteins as valuable markers for systematics and, due to the boom in complete genome sequencing, their use has become widespread. Yet, in the case of uncultured strains, for which complete genome sequences cannot be easily obtained, sequencing many markers is complicated and time consuming. Taking the advantage of the organization of ribosomal protein coding genes in large gene clusters, we amplified a 32 kb conserved region encompassing the spectinomycin (spc) operon using long range PCR from isolated and from uncultured nodular endophytic Frankia strains. The phylogenetic analysis of the 27 ribosomal protein genes contained in this region provided a robust phylogenetic tree consistent with phylogenies based on larger set of markers, indicating that this subset of ribosomal proteins contains enough phylogenetic signal to address systematic issues. This work shows that using long range PCR could break down the barrier preventing the use of ribosomal proteins as phylogenetic markers when complete genome sequences cannot be easily obtained.

Metabolic networks for nitrogen utilization in Prevotella ruminicola 23

Nitrogen metabolism in gut systems remains poorly studied in spite of its importance for microbial growth and its implications for the metabolism of the host. Prevotella spp. are the most predominant bacteria detected in the rumen, but their presence has also been related to health and disease states in the human gut and oral cavity. To explore the metabolic networks for nitrogen assimilation in this bacterium, changes in gene expression profiles in response to variations in the available nitrogen source and to different concentrations of ammonium were analyzed by microarray and reverse transcription quantitative PCR, and linked with function by further proteomic analysis. The observed patterns of transcript abundances for genes involved in ammonium assimilation differed from the classical “enteric paradigm” for nitrogen utilization. Expression of genes encoding high substrate affinity nitrogen assimilation enzymes (GS-GOGAT system) was similar in growth-limiting and non-limiting nitrogen concentrations in P. ruminicola 23, whereas E. coli and Salmonella spp. responses to excess nitrogen involve only low substrate affinity enzymes. This versatile behavior might be a key feature for ecological success in habitats such as the rumen and human colon where nitrogen is rarely limiting for growth, and might be linked to previously reported Prevotella spp. population imbalances relative to other bacterial species in gut systems.

Discrimination of Escherichia coli O157, O26 and O111 from other serovars by MALDI-TOF MS based on the S10-GERMS method

Enterohemorrhagic Escherichia coli (EHEC), causes a potentially life-threatening infection in humans worldwide. Serovar O157:H7, and to a lesser extent serovars O26 and O111, are the most commonly reported EHEC serovars responsible for a large number of outbreaks. We have established a rapid discrimination method for E. coli serovars O157, O26 and O111 from other E. coli serovars, based on the pattern matching of mass spectrometry (MS) differences and the presence/absence of biomarker proteins detected in matrix-assisted laser desorption/ionization time-of-flight MS (MALDI-TOF MS). Three biomarkers, ribosomal proteins S15 and L25, and acid stress chaperone HdeB, with MS m/z peaks at 10138.6/10166.6, 10676.4/10694.4 and 9066.2, respectively, were identified as effective biomarkers for O157 discrimination. To distinguish serovars O26 and O111 from the others, DNA-binding protein H-NS, with an MS peak at m/z 15409.4/15425.4 was identified. Sequence analysis of the O157 biomarkers revealed that amino acid changes: Q80R in S15, M50I in L25 and one mutation within the start codon ATG to ATA in the encoded HdeB protein, contributed to the specific peak pattern in O157. We demonstrated semi-automated pattern matching using these biomarkers and successfully discriminated total 57 O157 strains, 20 O26 strains and 6 O111 strains with 100% reliability by conventional MALDI-TOF MS analysis, regardless of the sample conditions. Our simple strategy, based on the S10-spc-alpha operon gene-encoded ribosomal protein mass spectrum (S10-GERMS) method, therefore allows for the rapid and reliable detection of this pathogen and may prove to be an invaluable tool both clinically and in the food industry.

Mobile elements in a single-filament orange Guaymas Basin Beggiatoa ("Candidatus Maribeggiatoa") sp. draft genome: evidence for genetic exchange with cyanobacteria

The draft genome sequence of a single orange Beggiatoa ("Candidatus Maribeggiatoa") filament collected from a microbial mat at a hydrothermal site in Guaymas Basin (Gulf of California, Mexico) shows evidence of extensive genetic exchange with cyanobacteria, in particular for sensory and signal transduction genes. A putative homing endonuclease gene and group I intron within the 23S rRNA gene; several group II catalytic introns; GyrB and DnaE inteins, also encoding homing endonucleases; multiple copies of sequences similar to the fdxN excision elements XisH and XisI (required for heterocyst differentiation in some cyanobacteria); and multiple sequences related to an open reading frame (ORF) (00024_0693) of unknown function all have close non-Beggiatoaceae matches with cyanobacterial sequences. Sequences similar to the uncharacterized ORF and Xis elements are found in other Beggiatoaceae genomes, a variety of cyanobacteria, and a few phylogenetically dispersed pleiomorphic or filamentous bacteria. We speculate that elements shared among filamentous bacterial species may have been exchanged in microbial mats and that some of them may be involved in cell differentiation.

Organisms can essentially be classified according to two codon patterns

We recently classified 23 bacteria into two types based on their complete genomes; "S-type" as represented by Staphylococcus aureus and "E-type" as represented by Escherichia coli. Classification was characterized by concentrations of Arg, Ala or Lys in the amino acid composition calculated from the complete genome. Based on these previous classifications, not only prokaryotic but also eukaryotic genome structures were investigated by amino acid compositions and nucleotide contents. Organisms consisting of 112 bacteria, 15 archaea and 18 eukaryotes were classified into two major groups by cluster analysis using GC contents at the three codon positions calculated from complete genomes. The 145 organisms were classified into "AT-type" and "GC-type" represented by high A or T (low G or C) and high G or C (low A or T) contents, respectively, at every third codon position. Reciprocal changes between G or C and A or T contents at the third codon position occurred almost synchronously in every codon among the organisms. Correlations between amino acid concentrations (Ala, Ile and Lys) and the nucleotide contents at the codon position were obtained in both "AT-type" and "GC-type" organisms, but with different regression coefficients. In certain correlations of amino acid concentrations with GC contents, eukaryotes, archaea and bacteria showed different behaviors; thus these kingdoms evolved differently. All organisms are basically classifiable into two groups having characteristic codon patterns; organisms with low GC and high AT contents at the third codon position and their derivatives, and organisms with an inverse relationship.

Gene teams: a new formalization of gene clusters for comparative genomics

This paper describes an efficient algorithm based on a new concept called gene team for detecting conserved gene clusters among an arbitrary number of chromosomes. Within the clusters, neither the order of the genes nor their orientation need be conserved. In addition, insertion of foreign genes within the clusters are permitted to a user-defined extent. This algorithm has been implemented in a publicly available TEAM software that proves to be an efficient tool for systematic searches of conserved gene clusters. Examples of actual biological results are provided. The software is downloadable from http://www-igm.univ-mlv.fr/ approximately raffinot/geneteam.html.