Evolution of mosaic operons by horizontal gene transfer and gene displacement in situ

Phenotypic and genomic diversity of Lactobacillus plantarum strains isolated from various environmental niches

Lactobacillus plantarum is a ubiquitous microorganism that is able to colonize several ecological niches, including vegetables, meat, dairy substrates and the gastro-intestinal tract. An extensive phenotypic and genomic diversity analysis was conducted to elucidate the molecular basis of the high flexibility and versatility of this species. First, 185 isolates from diverse environments were phenotypically characterized by evaluating their fermentation and growth characteristics. Strains clustered largely together within their particular food niche, but human fecal isolates were scattered throughout the food clusters, suggesting that they originate from the food eaten by the individuals. Based on distinct phenotypic profiles, 24 strains were selected and, together with a further 18 strains from an earlier low-resolution study, their genomic diversity was evaluated by comparative genome hybridization against the reference genome of L. plantarum WCFS1. Over 2000 genes were identified that constitute the core genome of the L. plantarum species, including 121 unique L. plantarum-marker genes that have not been found in other lactic acid bacteria. Over 50 genes unique for the reference strain WCFS1 were identified that were absent in the other L. plantarum strains. Strains of the L. plantarum subspecies argentoratensis were found to lack a common set of 24 genes, organized in seven gene clusters/operons, supporting their classification as a separate subspecies. The results provide a detailed view on phenotypic and genomic diversity of L. plantarum and lead to a better comprehension of niche adaptation and functionality of the organism.

Operons

Operons (clusters of co-regulated genes with related functions) are common features of bacterial genomes. More recently, functional gene clustering has been reported in eukaryotes, from yeasts to filamentous fungi, plants, and animals. Gene clusters can consist of paralogous genes that have most likely arisen by gene duplication. However, there are now many examples of eukaryotic gene clusters that contain functionally related but non-homologous genes and that represent functional gene organizations with operon-like features (physical clustering and co-regulation). These include gene clusters for use of different carbon and nitrogen sources in yeasts, for production of antibiotics, toxins, and virulence determinants in filamentous fungi, for production of defense compounds in plants, and for innate and adaptive immunity in animals (the major histocompatibility locus). The aim of this article is to review features of functional gene clusters in prokaryotes and eukaryotes and the significance of clustering for effective function.

Horizontal gene transfer and diverse functional constrains within a common replication-partitioning system in Alphaproteobacteria: the repABC operon

Background

The repABC plasmid family, which is extensively present within Alphaproteobacteria, and some secondary chromosomes of the Rhizobiales have the particular feature that all the elements involved in replication and partitioning reside within one transcriptional unit, the repABC operon. Given the functional interactions among the elements of the repABC operon, and the fact that they all reside in the same operon, a common evolutionary history would be expected if the entire operon had been horizontally transferred. Here, we tested whether there is a common evolutionary history within the repABC operon. We further examined different incompatibility groups in terms of their differentiation and degree of adaptation to their host.

Results

We did not find a single evolutionary history within the repABC operon. Each protein had a particular phylogeny, horizontal gene transfer events of the individual genes within the operon were detected, and different functional constraints were found within and between the Rep proteins. When different repABC operons coexisted in the same genome, they were well differentiated from one another. Finally, we found different levels of adaptation to the host genome within and between repABC operons coexisting in the same species.

Conclusion

Horizontal gene transfer with conservation of the repABC operon structure provides a highly dynamic operon in which each member of this operon has its own evolutionary dynamics. In addition, it seems that different incompatibility groups present in the same species have different degrees of adaptation to their host genomes, in proportion to the amount of time the incompatibility group has coexisted with the host genome.

Lateral transfer of genes and gene fragments in prokaryotes

Lateral genetic transfer (LGT) involves the movement of genetic material from one lineage into another and its subsequent incorporation into the new host genome via genetic recombination. Studies in individual taxa have indicated lateral origins for stretches of DNA of greatly varying length, from a few nucleotides to chromosome size. Here we analyze 1,462 sets of single-copy, putatively orthologous genes from 144 fully sequenced prokaryote genomes, asking to what extent complete genes and fragments of genes have been transferred and recombined in LGT. Using a rigorous phylogenetic approach, we find evidence for LGT in at least 476 (32.6%) of these 1,462 gene sets: 286 (19.6%) clearly show one or more "observable recombination breakpoints" within the boundaries of the open reading frame, while a further 190 (13.0%) yield trees that are topologically incongruent with the reference tree but do not contain a recombination breakpoint within the open reading frame. We refer to these gene sets as observable recombination breakpoint positive (ORB(+)) and negative (ORB(-)) respectively. The latter are prima facie instances of lateral transfer of an entire gene or beyond. We observe little functional bias between ORB(+) and ORB(-) gene sets, but find that incorporation of entire genes is potentially more frequent in pathogens than in nonpathogens. As ORB(+) gene sets are about 50% more common than ORB(-) sets in our data, the transfer of gene fragments has been relatively frequent, and the frequency of LGT may have been systematically underestimated in phylogenetic studies.

Expression of phosphofructokinase in Neisseria meningitidis

Neisseria meningitidis serogroup B is a pathogen that can infect diverse sites within the human host. According to the N. meningitidis genomic information and experimental observations, glucose can be completely catabolized through the Entner-Doudoroff pathway and the pentose phosphate pathway. The Embden-Meyerhof-Parnas pathway is not functional, because the gene for phosphofructokinase (PFK) is not present. The phylogenetic distribution of PFK indicates that in most obligate aerobic organisms, PFK is lacking. We conclude that this is because of the limited contribution of PFK to the energy supply in aerobically grown organisms in comparison with the energy generated through oxidative phosphorylation. Under anaerobic or microaerobic conditions, the available energy is limiting, and PFK provides an advantage, which explains the presence of PFK in many (facultatively) anaerobic organisms. In accordance with this, in silico flux balance analysis predicted an increase of biomass yield as a result of PFK expression. However, analysis of a genetically engineered N. meningitidis strain that expressed a heterologous PFK showed that the yield of biomass on substrate decreased in comparison with a pfkA-deficient control strain, which was associated mainly with an increase in CO(2) production, whereas production of by-products was similar in the two strains. This might explain why the pfkA gene has not been obtained by horizontal gene transfer, since it is initially unfavourable for biomass yield. No large effects related to heterologous expression of pfkA were observed in the transcriptome. Although our results suggest that introduction of PFK does not contribute to a more efficient strain in terms of biomass yield, achievement of a robust, optimal metabolic network that enables a higher growth rate or a higher biomass yield might be possible after adaptive evolution of the strain, which remains to be investigated.

Evolutionary constraints permeate large metabolic networks

BACKGROUND

Metabolic networks show great evolutionary plasticity, because they can differ substantially even among closely related prokaryotes. Any one metabolic network can also effectively compensate for the blockage of individual reactions by rerouting metabolic flux through other pathways. These observations, together with the continual discovery of new microbial metabolic pathways and enzymes, raise the possibility that metabolic networks are only weakly constrained in changing their complement of enzymatic reactions.

RESULTS

To ask whether this is the case, I characterized pairwise and higher-order associations in the co-occurrence of genes encoding metabolic enzymes in more than 200 completely sequenced representatives of prokaryotic genera. The majority of reactions show constrained evolution. Specifically, genes encoding most reactions tend to co-occur with genes encoding other reaction(s). Constrained reaction pairs occur in small sets whose number is substantially greater than expected by chance alone. Most such sets are associated with single biochemical pathways. The respective genes are not always tightly linked, which renders horizontal co-transfer of constrained reaction sets an unlikely sole cause for these patterns of association.

CONCLUSION

Even a limited number of available genomes suffices to show that metabolic network evolution is highly constrained by reaction combinations that are favored by natural selection. With increasing numbers of completely sequenced genomes, an evolutionary constraint-based approach may enable a detailed characterization of co-evolving metabolic modules.

Horizontal gene transfer of zinc and non-zinc forms of bacterial ribosomal protein S4

Background

The universal ribosomal protein S4 is essential for the initiation of small subunit ribosomal assembly and translational accuracy. Being part of the information processing machinery of the cell, the gene for S4 is generally thought of as being inherited vertically and has been used in concatenated gene phylogenies. Here we report the evolution of ribosomal protein S4 in relation to a broad sharing of zinc/non-zinc forms of the gene and study the scope of horizontal gene transfer (HGT) of S4 during bacterial evolution.

Results

In this study we present the complex evolutionary history of ribosomal protein S4 using 660 bacterial genomes from 16 major bacterial phyla. According to conserved characteristics in the sequences, S4 can be classified into C+ (zinc-binding) and C- (zinc-free) variants, with 26 genomes (mainly from the class Clostridia) containing genes for both. A maximum likelihood phylogenetic tree of the S4 sequences was incongruent with the standard bacterial phylogeny, indicating a departure from strict vertical inheritance. Further analysis using the genome content near the S4 genes, which are usually located in a conserved gene cluster, showed not only that HGT of the C- gene had occurred at various stages of bacterial evolution, but also that both the C- and C+ genes were present before the individual phyla diverged. To explain the latter, we theorize that a gene pool existed early in bacterial evolution from which bacteria could sample S4 gene variants, according to environmental conditions. The distribution of the C+/- variants for seven other zinc-binding ribosomal proteins in these 660 bacterial genomes is consistent with that seen for S4 and may shed light on the evolutionary pressures involved.

Conclusion

The complex history presented for "core" protein S4 suggests the existence of a gene pool before the emergence of bacterial lineages and reflects the pervasive nature of HGT in subsequent bacterial evolution. This has implications for both theoretical models of evolution and practical applications of phylogenetic reconstruction as well as the control of zinc economy in bacterial cells.

Functional replacement of a primary metabolic pathway via multiple independent eukaryote-to-eukaryote gene transfers and selective retention

Although lateral gene transfer (LGT) is now recognized as a major force in the evolution of prokaryotes, the contribution of LGT to the evolution and diversification of eukaryotes is less understood. Notably, transfers of complete pathways are believed to be less likely between eukaryotes, because the successful transfer of a pathway requires the physical clustering of functionally related genes. Here, we report that in one of the closest unicellular relatives of animals, the choanoflagellate, Monosiga, three genes whose products work together in the glutamate synthase cycle are of algal origin. The concerted retention of these three independently acquired genes is best explained as the consequence of a series of adaptive replacement events. More generally, this study argues that (i) eukaryote-to-eukaryote transfers of entire metabolic pathways are possible, (ii) adaptive functional replacements of primary pathways can occur, and (iii) functional replacements involving eukaryotic genes are likely to have also contributed to the evolution of eukaryotes. Lastly, these data underscore the potential contribution of algal genes to the evolution of nonphotosynthetic lineages.

Does gene translocation accelerate the evolution of laterally transferred genes?

Lateral gene transfer (LGT) and gene rearrangement are essential for shaping bacterial genomes during evolution. Separate attention has been focused on understanding the process of lateral gene transfer and the process of gene translocation. However, little is known about how gene translocation affects laterally transferred genes. Here we have examined gene translocations and lateral gene transfers in closely related genome pairs. The results reveal that translocated genes undergo elevated rates of evolution and gene translocation tends to take place preferentially in recently acquired genes. Translocated genes have a high probability to be truncated, suggesting that translocation followed by truncation/deletion might play an important role in the fast turnover of laterally transferred genes. Furthermore, more recently acquired genes have a higher proportion of genes on the leading strand, suggesting a strong strand bias of lateral gene transfer.

Insights into the evolution of sialic acid catabolism among bacteria

Background

Sialic acids comprise a family of nine-carbon amino sugars that are prevalent in mucus rich environments. Sialic acids from the human host are used by a number of pathogens as an energy source. Here we explore the evolution of the genes involved in the catabolism of sialic acid.

Results

The cluster of genes encoding the enzymes N-acetylneuraminate lyase (NanA), epimerase (NanE), and kinase (NanK), necessary for the catabolism of sialic acid (the Nan cluster), are confined 46 bacterial species, 42 of which colonize mammals, 33 as pathogens and 9 as gut commensals. We found a putative sialic acid transporter associated with the Nan cluster in most species. We reconstructed the phylogenetic history of the NanA, NanE, and NanK proteins from the 46 species and compared them to the species tree based on 16S rRNA. Within the NanA phylogeny, Gram-negative and Gram-positive bacteria do not form distinct clades. NanA from Yersinia and Vibrio species was most closely related to the NanA clade from eukaryotes. To examine this further, we reconstructed the phylogeny of all NanA homologues in the databases. In this analysis of 83 NanA sequences, Bacteroidetes, a human commensal group formed a distinct clade with Verrucomicrobia, and branched with the Eukaryotes and the Yersinia/Vibrio clades. We speculate that pathogens such as V. cholerae may have acquired NanA from a commensal aiding their colonization of the human gut. Both the NanE and NanK phylogenies more closely represented the species tree but numerous incidences of incongruence are noted. We confirmed the predicted function of the sialic acid catabolism cluster in members the major intestinal pathogens Salmonella enterica, Vibrio cholerae, V. vulnificus, Yersinia enterocolitica and Y. pestis.

Conclusion

The Nan cluster among bacteria is confined to human pathogens and commensals conferring them the ability to utilize a ubiquitous carbon source in mucus rich surfaces of the human body. The Nan region shows a mosaic evolution with NanA from Bacteroidetes, Vibrio and Yersinia branching closely together with NanA from eukaryotes.

Stochasticity in protein levels drives colinearity of gene order in metabolic operons of Escherichia coli

Gene order in some bacterial metabolic operons reflects ordering in the metabolic pathway. That this is true uniquely for operons expressed at low levels highlights the selective importance of fluctuations in protein levels.

In bacterial genomes, gene order is not random. This is most evident when looking at operons, these often encoding enzymes involved in the same metabolic pathway or proteins from the same complex. Is gene order within operons nonrandom, however, and if so why? We examine this issue using metabolic operons as a case study. Using the metabolic network of Escherichia coli, we define the temporal order of reactions. We find a pronounced trend for genes to appear in operons in the same order as they are needed in metabolism (colinearity). This is paradoxical as, at steady state, enzymes abundance should be independent of order within the operon. We consider three extensions of the steady-state model that could potentially account for colinearity: (1) increased productivity associated with higher expression levels of the most 5′ genes, (2) a faster metabolic processing immediately after up-regulation, and (3) metabolic stalling owing to stochastic protein loss. We establish the validity of these hypotheses by employing deterministic and stochastic models of enzyme kinetics. The stochastic stalling hypothesis correctly and uniquely predicts that colinearity is more pronounced both for lowly expressed operons and for genes that are not physically adjacent. The alternative models fail to find any support. These results support the view that stochasticity is a pervasive problem to a cell and that gene order evolution can be driven by the selective consequences of fluctuations in protein levels.

In bacteria, different enzymes from the same metabolic pathway are often encoded within one transcriptional unit, an operon. There is also, we show, a tendency for the enzymes that are needed earlier in the pathway to feature earlier in the operon, so-called colinearity. Why might this be? We test three ideas, one old and two new. The prior suggestion supposes that proteins of genes early in operons will be at a higher dose. Although some operons are like this, in general, we see no relationship of protein dose with colinearity. We also find no evidence that operons that frequently need up-regulation are any more likely to be colinear. A third model is, however, supported. If an operon is rarely expressed, then all the proteins for this part of metabolism can be lost by chance. Rebooting such metabolism is fastest if the operon is colinear. This model predicts, correctly, that colinearity should be more frequent in operons that are expressed at a low level. This result is important for at least two reasons. First, it supports the view that chance events (such as protein loss) within cells are important on a day-to-day basis. Second, it challenges the supposition that natural selection will be weakest on lowly expressed genes. Where chance events are concerned, natural selection can be strong on genes expressed at a low level.

GO4genome: a prokaryotic phylogeny based on genome organization

Determining the phylogeny of closely related prokaryotes may fail in an analysis of rRNA or a small set of sequences. Whole-genome phylogeny utilizes the maximally available sample space. For a precise determination of genome similarity, two aspects have to be considered when developing an algorithm of whole-genome phylogeny: (1) gene order conservation is a more precise signal than gene content; and (2) when using sequence similarity, failures in identifying orthologues or the in situ replacement of genes via horizontal gene transfer may give misleading results. GO4genome is a new paradigm, which is based on a detailed analysis of gene function and the location of the respective genes. For characterization of genes, the algorithm uses gene ontology enabling a comparison of function independent of evolutionary relationship. After the identification of locally optimal series of gene functions, their length distribution is utilized to compute a phylogenetic distance. The outcome is a classification of genomes based on metabolic capabilities and their organization. Thus, the impact of effects on genome organization that are not covered by methods of molecular phylogeny can be studied. Genomes of strains belonging to Escherichia coli, Shigella, Streptococcus, Methanosarcina, and Yersinia were analyzed. Differences from the findings of classical methods are discussed.

Interdomain transfers of sugar transporters overcome barriers to gene expression

Horizontal gene transfer (HGT) is nature's mechanism for sharing evolved physiological traits among the members of microbial communities. The extent to which such transfers can be successful is best illustrated by the fact that Archaea-derived genes are found in many bacterial genomes, particularly those in the hyperthermophile Thermotoga maritima. The success of these intergenomic transfers depends upon the successful transcription of the newly acquired archaeal genes using a bacterial transcription machinery that does not recognize archaeal transcriptional signals. To examine how nature solves this problem, we looked to the T. maritima genome for examples of interdomain transfers. Here we lay the groundwork to examine this problem by more clearly delineating the phylogenetic history of Archaea-derived transporter genes in this genome. We find that five of these polysaccharide transporters were derived from the Archaea and one came from the Archaea after that lineage inherited it from the Bacteria. These data can be used for more detailed examinations of the recombinations that allowed these transporters to be expressed in a bacterial host. This work will guide examinations of the genome sequences from other members of the Thermotogales, which will become available.

Are protein domains modules of lateral genetic transfer?

Background

In prokaryotes and some eukaryotes, genetic material can be transferred laterally among unrelated lineages and recombined into new host genomes, providing metabolic and physiological novelty. Although the process is usually framed in terms of gene sharing (e.g. lateral gene transfer, LGT), there is little reason to imagine that the units of transfer and recombination correspond to entire, intact genes. Proteins often consist of one or more spatially compact structural regions (domains) which may fold autonomously and which, singly or in combination, confer the protein's specific functions. As LGT is frequent in strongly selective environments and natural selection is based on function, we hypothesized that domains might also serve as modules of genetic transfer, i.e. that regions of DNA that are transferred and recombined between lineages might encode intact structural domains of proteins.

Methodology/Principal Findings

We selected 1,462 orthologous gene sets representing 144 prokaryotic genomes, and applied a rigorous two-stage approach to identify recombination breakpoints within these sequences. Recombination breakpoints are very significantly over-represented in gene sets within which protein domain-encoding regions have been annotated. Within these gene sets, breakpoints significantly avoid the domain-encoding regions (domons), except where these regions constitute most of the sequence length. Recombination breakpoints that fall within longer domons are distributed uniformly at random, but those that fall within shorter domons may show a slight tendency to avoid the domon midpoint. As we find no evidence for differential selection against nucleotide substitutions following the recombination event, any bias against disruption of domains must be a consequence of the recombination event per se.

Conclusions/Significance

This is the first systematic study relating the units of LGT to structural features at the protein level. Many genes have been interrupted by recombination following inter-lineage genetic transfer, during which the regions within these genes that encode protein domains have not been preferentially preserved intact. Protein domains are units of function, but domons are not modules of transfer and recombination. Our results demonstrate that LGT can remodel even the most functionally conservative modules within genomes.

Recurring cluster and operon assembly for Phenylacetate degradation genes

Background

A large number of theories have been advanced to explain why genes involved in the same biochemical processes are often co-located in genomes. Most of these theories have been dismissed because empirical data do not match the expectations of the models. In this work we test the hypothesis that cluster formation is most likely due to a selective pressure to gradually co-localise protein products and that operon formation is not an inevitable conclusion of the process.

Results

We have selected an exemplar well-characterised biochemical pathway, the phenylacetate degradation pathway, and we show that its complex history is only compatible with a model where a selective advantage accrues from moving genes closer together. This selective pressure is likely to be reasonably weak and only twice in our dataset of 102 genomes do we see independent formation of a complete cluster containing all the catabolic genes in the pathway. Additionally, de novo clustering of genes clearly occurs repeatedly, even though recombination should result in the random dispersal of such genes in their respective genomes. Interspecies gene transfer has frequently replaced in situ copies of genes resulting in clusters that have similar content but very different evolutionary histories.

Conclusion

Our model for cluster formation in prokaryotes, therefore, consists of a two-stage selection process. The first stage is selection to move genes closer together, either because of macromolecular crowding, chromatin relaxation or transcriptional regulation pressure. This proximity opportunity sets up a separate selection for co-transcription.

The fate of new bacterial genes

Bacteria experience a continual influx of novel genetic material from a wide range of sources and yet their genomes remain relatively small. This aspect of bacterial evolution indicates that most newly arriving sequences are rapidly eliminated; however, numerous new genes persist, as evident from the presence of unique genes in almost all bacterial genomes. This review summarizes the methods for identifying new genes in bacterial genomes and examines the features that promote the retention and elimination of these evolutionary novelties.

NAD biosynthesis evolution in bacteria: lateral gene transfer of kynurenine pathway in Xanthomonadales and Flavobacteriales

The biosynthesis of quinolinate, the de novo precursor of nicotinamide adenine dinucleotide (NAD), may be performed by two distinct pathways, namely, the bacterial aspartate (aspartate-to-quinolinate) and the eukaryotic kynurenine (tryptophan-to-quinolinate). Even though the separation into eukaryotic and bacterial routes is long established, recent genomic surveys have challenged this view, because certain bacterial species also carry the genes for the kynurenine pathway. In this work, both quinolinate biosynthetic pathways were investigated in the Bacteria clade and with special attention to Xanthomonadales and Bacteroidetes, from an evolutionary viewpoint. Genomic screening has revealed that a small number of bacterial species possess some of the genes for the kynurenine pathway, which is complete in the genus Xanthomonas and in the order Flavobacteriales, where the aspartate pathway is absent. The opposite pattern (presence of the aspartate pathway and absence of the kynurenine pathway) in close relatives (Xylella ssp. and the order Bacteroidales, respectively) points to the idea of a recent acquisition of the kynurenine pathway through lateral gene transfer in these bacterial groups. In fact, sequence similarity comparison and phylogenetic reconstruction both suggest that at least part of the genes of the kynurenine pathway in Xanthomonas and Flavobacteriales is shared by eukaryotes. These results reinforce the idea of the role that lateral gene transfer plays in the configuration of bacterial genomes, thereby providing alternative metabolic pathways, even with the replacement of primary and essential cell functions, as exemplified by NAD biosynthesis.

Metagenomic and stable isotopic analyses of modern freshwater microbialites in Cuatro Ciénegas, Mexico

Ancient biologically mediated sedimentary carbonate deposits, including stromatolites and other microbialites, provide insight into environmental conditions on early Earth. The primary limitation to interpreting these records is our lack of understanding regarding microbial processes and the preservation of geochemical signatures in contemporary microbialite systems. Using a combination of metagenomic sequencing and isotopic analyses, this study describes the identity, metabolic potential and chemical processes of microbial communities from living microbialites from Cuatro Ciénegas, Mexico. Metagenomic sequencing revealed a diverse, redox-dependent microbial community associated with the microbialites. The microbialite community is distinct from other marine and freshwater microbial communities, and demonstrates extensive environmental adaptation. The microbialite metagenomes contain a large number of genes involved in the production of exopolymeric substances and the formation of biofilms, creating a complex, spatially structured environment. In addition to the spatial complexity of the biofilm, microbial activity is tightly controlled by sensory and regulatory systems, which allow for coordination of autotrophic and heterotrophic processes. Isotopic measurements of the intracrystalline organic matter demonstrate the importance of heterotrophic respiration of photoautotrophic biomass in the precipitation of calcium carbonate. The genomic and stable isotopic data presented here significantly enhance our evolving knowledge of contemporary biomineralization processes, and are directly applicable to studies of ancient microbialites.

Phylogeny vs genome reshuffling: horizontal gene transfer

The evolutionary events in organisms can be tracked to the transfer of genetic material. The inheritance of genetic material among closely related organisms is a slow evolutionary process. On the other hand, the movement of genes among distantly related species can account for rapid evolution. The later process has been quite evident in the appearance of antibiotic resistance genes among human and animal pathogens. Phylogenetic trees based on such genes and those involved in metabolic activities reflect the incongruencies in comparison to the 16S rDNA gene, generally used for taxonomic relationships. Such discrepancies in gene inheritance have been termed as horizontal gene transfer (HGT) events. In the post-genomic era, the explosion of known sequences through large-scale sequencing projects has unraveled the weakness of traditional 16S rDNA gene tree based evolutionary model. Various methods to scrutinize HGT events include atypical composition, abnormal sequence similarity, anomalous phylogenetic distribution, unusual phyletic patterns, etc. Since HGT generates greater genetic diversity, it is likely to increase resource use and ecosystem resilience.

Comparative genomic analysis of T-box regulatory systems in bacteria

T-box antitermination is one of the main mechanisms of regulation of genes involved in amino acid metabolism in Gram-positive bacteria. T-box regulatory sites consist of conserved sequence and RNA secondary structure elements. Using a set of known T-box sites, we constructed the common pattern and used it to scan available bacterial genomes. New T-boxes were found in various Gram-positive bacteria, some Gram-negative bacteria (delta-proteobacteria), and some other bacterial groups (Deinococcales/Thermales, Chloroflexi, Dictyoglomi). The majority of T-box-regulated genes encode aminoacyl-tRNA synthetases. Two other groups of T-box-regulated genes are amino acid biosynthetic genes and transporters, as well as genes with unknown function. Analysis of candidate T-box sites resulted in new functional annotations. We assigned the amino acid specificity to a large number of candidate amino acid transporters and a possible function to amino acid biosynthesis genes. We then studied the evolution of the T-boxes. Analysis of the constructed phylogenetic trees demonstrated that in addition to the normal evolution consistent with the evolution of regulated genes, T-boxes may be duplicated, transferred to other genes, and change specificity. We observed several cases of recent T-box regulon expansion following the loss of a previously existing regulatory system, in particular, arginine regulon in Clostridium difficile and methionine regulon in Lactobacillaceae. Finally, we described a new structural class of T-boxes containing duplicated terminator-antiterminator elements and unusual reduced T-boxes regulating initiation of translation in the Actinobacteria.

The impact of horizontal gene transfer in shaping operons and protein interaction networks--direct evidence of preferential attachment

BACKGROUND

Despite the prevalence of horizontal gene transfer (HGT) in bacteria, to this date there were few studies on HGT in the context of gene expression, operons and protein-protein interactions. Using the recently available data set on the E. coli protein-protein interaction network, we sought to explore the impact of HGT on genome structure and protein networks.

RESULTS

We classified the E. coli genes into three categories based on their evolutionary conservation: a set of 2158 Core genes that are shared by all E. coli strains, a set of 1044 Non-core genes that are strain-specific, and a set of 1053 genes that were putatively acquired by horizontal transfer. We observed a clear correlation between gene expressivity (measured by Codon Adaptation Index), evolutionary rates, and node connectivity between these categories of genes. Specifically, we found the Core genes are the most highly expressed and the most slowly evolving, while the HGT genes are expressed at the lowest level and evolve at the highest rate. Core genes are the most likely and HGT genes are the least likely to be member of the operons. In addition, we found the Core genes on average are more highly connected than Non-core and HGT genes in the protein interaction network, however the HGT genes displayed a significantly higher mean node degree than the Core and Non-core genes in the defence COG functional category. Interestingly, HGT genes are more likely to be connected to Core genes than expected by chance, which suggest a model of differential attachment in the expansion of cellular networks.

CONCLUSION

Results from our analysis shed light on the mode and mechanism of the integration of horizontally transferred genes into operons and protein interaction networks.

Double illegitimate recombination events integrate DNA segments through two different mechanisms during natural transformation of Acinetobacter baylyi

Acquisition of foreign DNA by horizontal gene transfer is seen as a major source of genetic diversity in prokaryotes. However, strongly divergent DNA is not genomically integrated by homologous recombination and would depend on illegitimate recombination (IR) events which are rare. We show that, by two mechanisms, during natural transformation of Acinetobacter baylyi two IR events can integrate DNA segments. One mechanism is double illegitimate recombination (DIR) acting in the absence of any homology (frequency: 7 x 10(-13) per cell). It occurs about 10(10)-fold less frequent than homologous transformation. The other mechanism is homology-facilitated double illegitimate recombination (HFDIR) being about 440-fold more frequent (3 x 10(-10) per cell) than DIR. HFDIR depends on a homologous sequence located between the IR sites and on recA(+). In HFDIR two IR events act on the same donor DNA molecule as shown by the joint inheritance of molecular DNA tags. While the IR events in HFDIR occurred at microhomologies, in DIR microhomologies were not used. The HFDIR phenomenon indicates that a temporal recA-dependent association of donor DNA at a homology in recipient DNA may facilitate two IR events on the 5' and 3' heterologous parts of the transforming DNA molecule.

Being pathogenic, plastic, and sexual while living with a nearly minimal bacterial genome

Mycoplasmas are commonly described as the simplest self-replicating organisms, whose evolution was mainly characterized by genome downsizing with a proposed evolutionary scenario similar to that of obligate intracellular bacteria such as insect endosymbionts. Thus far, analysis of mycoplasma genomes indicates a low level of horizontal gene transfer (HGT) implying that DNA acquisition is strongly limited in these minimal bacteria. In this study, the genome of the ruminant pathogen Mycoplasma agalactiae was sequenced. Comparative genomic data and phylogenetic tree reconstruction revealed that ∼18% of its small genome (877,438 bp) has undergone HGT with the phylogenetically distinct mycoides cluster, which is composed of significant ruminant pathogens. HGT involves genes often found as clusters, several of which encode lipoproteins that usually play an important role in mycoplasma–host interaction. A decayed form of a conjugative element also described in a member of the mycoides cluster was found in the M. agalactiae genome, suggesting that HGT may have occurred by mobilizing a related genetic element. The possibility of HGT events among other mycoplasmas was evaluated with the available sequenced genomes. Our data indicate marginal levels of HGT among Mycoplasma species except for those described above and, to a lesser extent, for those observed in between the two bird pathogens, M. gallisepticum and M. synoviae. This first description of large-scale HGT among mycoplasmas sharing the same ecological niche challenges the generally accepted evolutionary scenario in which gene loss is the main driving force of mycoplasma evolution. The latter clearly differs from that of other bacteria with small genomes, particularly obligate intracellular bacteria that are isolated within host cells. Consequently, mycoplasmas are not only able to subvert complex hosts but presumably have retained sexual competence, a trait that may prevent them from genome stasis and contribute to adaptation to new hosts.

Mycoplasmas are cell wall–lacking prokaryotes that evolved from ancestors common to Gram-positive bacteria by way of massive losses of genetic material. With their minimal genome, mycoplasmas are considered to be the simplest free-living organisms, yet several species are successful pathogens of man and animal. In this study, we challenged the commonly accepted view in which mycoplasma evolution is driven only by genome down-sizing. Indeed, we showed that a significant amount of genes underwent horizontal transfer among different mycoplasma species that share the same ruminant hosts. In these species, the occurrence of a genetic element that can promote DNA transfer via cell-to-cell contact suggests that some mycoplasmas may have retained or acquired sexual competence. Transferred genes were found to encode proteins that are likely to be associated with mycoplasma–host interactions. Sharing genetic resources via horizontal gene transfer may provide mycoplasmas with a means for adapting to new niches or to new hosts and for avoiding irreversible genome erosion.

Production of the non-ribosomal peptide plipastatin in Bacillus subtilis regulated by three relevant gene blocks assembled in a single movable DNA segment

Methods that allow the assembly of genes in one single DNA segment are of great use in bioengineering and synthetic biology. The biosynthesis of plipastatin, a lipopeptide antibiotic synthesized non-ribosomally by Bacillus subtilis 168, requires three gene blocks at different genome loci, i.e. the peptide synthetase operon ppsABCDE (38-kb), degQ (0.6kb), and sfp (1.0kb). We applied a DNA assembly protocol in B. subtilis, named ordered gene assembly in B. subtilis (OGAB) method, to incorporate those three gene blocks into a one-unit plasmid via one ligation-reaction. High yields of correct assembly, above 87%, allowed us to screen for the plasmid that produced plipastatin at a level approximately 10-fold higher than in the wild-type. In contrast to that recombinogenic technologies used in E. coli require repetitive assembly steps and/or several selection markers, our method features high fidelity and efficiency, is completed in one ligation using only one selection marker associating with plasmid vector, and is applicable to DNA fragments larger than 40kb.

Proteorhodopsin photosystem gene clusters exhibit co-evolutionary trends and shared ancestry among diverse marine microbial phyla

Since the recent discovery of retinylidene proteins in marine bacteria (proteorhodopsins), the estimated abundance and diversity of this gene family has expanded rapidly. To explore proteorhodopsin photosystem evolutionary and distributional trends, we identified and compared 16 different proteorhodopsin-containing genome fragments recovered from naturally occurring bacterioplankton populations. In addition to finding several deep-branching proteorhodopsin sequences, proteorhodopsins were found in novel taxonomic contexts, including a betaproteobacterium and a planctomycete. Approximately one-third of the proteorhodopsin-containing genome fragments analysed, as well as a number of recently reported marine bacterial whole genome sequences, contained a linked set of genes required for biosynthesis of the rhodopsin chromophore, retinal. Phylogenetic analyses of the retinal biosynthetic genes suggested their co-evolution and probable coordinated lateral gene transfer into disparate lineages, including Euryarchaeota, Planctomycetales, and three different proteobacterial lineages. The lateral transfer and retention of genes required to assemble a functional proteorhodopsin photosystem appears to be a coordinated and relatively frequent evolutionary event. Strong selection pressure apparently acts to preserve these light-dependent photosystems in diverse marine microbial lineages.

Modelling the evolution of the archeal tryptophan synthase

Background

Microorganisms and plants are able to produce tryptophan. Enzymes catalysing the last seven steps of tryptophan biosynthesis are encoded in the canonical trp operon. Among the trp genes are most frequently trpA and trpB, which code for the alpha and beta subunit of tryptophan synthase. In several prokaryotic genomes, two variants of trpB (named trpB1 or trpB2) occur in different combinations. The evolutionary history of these trpB genes is under debate.

Results

In order to study the evolution of trp genes, completely sequenced archeal and bacterial genomes containing trpB were analysed. Phylogenetic trees indicated that TrpB sequences constitute four distinct groups; their composition is in agreement with the location of respective genes. The first group consisted exclusively of trpB1 genes most of which belonged to trp operons. Groups two to four contained trpB2 genes. The largest group (trpB2_o) contained trpB2 genes all located outside of operons. Most of these genes originated from species possessing an operon-based trpB1 in addition. Groups three and four pertain to trpB2 genes of those genomes containing exclusively one or two trpB2 genes, but no trpB1. One group (trpB2_i) consisted of trpB2 genes located inside, the other (trpB2_a) of trpB2 genes located outside the trp operon. TrpA and TrpB form a heterodimer and cooperate biochemically. In order to characterise trpB variants and stages of TrpA/TrpB cooperation in silico, several approaches were combined. Phylogenetic trees were constructed for all trp genes; their structure was assessed via bootstrapping. Alternative models of trpB evolution were evaluated with parsimony arguments. The four groups of trpB variants were correlated with archeal speciation. Several stages of TrpA/TrpB cooperation were identified and trpB variants were characterised. Most plausibly, trpB2 represents the predecessor of the modern trpB gene, and trpB1 evolved in an ancestral bacterium.

Conclusion

In archeal genomes, several stages of trpB evolution, TrpA/TrpB cooperation, and operon formation can be observed. Thus, archeal trp genes may serve as a model system for studying the evolution of protein-protein interactions and operon formation.

The power of phylogenetic approaches to detect horizontally transferred genes

BACKGROUND

Horizontal gene transfer plays an important role in evolution because it sometimes allows recipient lineages to adapt to new ecological niches. High genes transfer frequencies were inferred for prokaryotic and early eukaryotic evolution. Does horizontal gene transfer also impact phylogenetic reconstruction of the evolutionary history of genomes and organisms? The answer to this question depends at least in part on the actual gene transfer frequencies and on the ability to weed out transferred genes from further analyses. Are the detected transfers mainly false positives, or are they the tip of an iceberg of many transfer events most of which go undetected by current methods?

RESULTS

Phylogenetic detection methods appear to be the method of choice to infer gene transfers, especially for ancient transfers and those followed by orthologous replacement. Here we explore how well some of these methods perform using in silico transfers between the terminal branches of a gamma proteobacterial, genome based phylogeny. For the experiments performed here on average the AU test at a 5% significance level detects 90.3% of the transfers and 91% of the exchanges as significant. Using the Robinson-Foulds distance only 57.7% of the exchanges and 60% of the donations were identified as significant. Analyses using bipartition spectra appeared most successful in our test case. The power of detection was on average 97% using a 70% cut-off and 94.2% with 90% cut-off for identifying conflicting bipartitions, while the rate of false positives was below 4.2% and 2.1% for the two cut-offs, respectively. For all methods the detection rates improved when more intervening branches separated donor and recipient.

CONCLUSION

Rates of detected transfers should not be mistaken for the actual transfer rates; most analyses of gene transfers remain anecdotal. The method and significance level to identify potential gene transfer events represent a trade-off between the frequency of erroneous identification (false positives) and the power to detect actual transfer events.

Positional homology in bacterial genomes

In comparative genomic studies, syntenic groups of homologous sequence in the same order have been used as supplementary information that can be used in helping to determine the orthology of the compared sequences. The assumption is that orthologous gene copies are more likely to share the same genome positions and share the same gene neighbors. In this study we have defined positional homologs as those that also have homologous neighboring genes and we investigated the usefulness of this distinction for bacterial comparative genomics. We considered the identification of positionaly homologous gene pairs in bacterial genomes using protein and DNA sequence level alignments and found that the positional homologs had on average relatively lower rates of substitution at the DNA level (synonymous substitutions) than duplicate homologs in different genomic locations, regardless of the level of protein sequence divergence (measured with non-synonymous substitution rate). Since gene order conservation can indicate accuracy of orthology assignments, we also considered the effect of imposing certain alignment quality requirements on the sensitivity and specificity of identification of protein pairs by BLAST and FASTA when neighboring information is not available and in comparisons where gene order is not conserved. We found that the addition of a stringency filter based on the second best hits was an efficient way to remove dubious ortholog identifications in BLAST and FASTA analyses. Gene order conservation and DNA sequence homology are useful to consider in comparative genomic studies as they may indicate different orthology assignments than protein sequence homology alone.

Gene Cluster Analysis Method Identifies Horizontally Transferred Genes with High Reliability and Indicates that They Provide the Main Mechanism of Operon Gain in 8 Species of γ-Proteobacteria

The formation mechanism of operons remains unresolved: operons may form by rearrangements within a genome or by acquisition of genes from other species, that is, horizontal gene transfer (HGT). One hindrance to its elucidation is the unavailability of a method to accurately identify HGT, although it is generally considered to occur. It is critically important first to select horizontally transferred (HT) genes reliably and then to determine the extent to which HGT is involved in operon formation. For this purpose, we considered indels in terms of gene clusters instead of individual genes and chose candidates of HT genes in 8 species of Escherichia, Shigella, and Salmonella based on the minimization of indels. To select a benchmark set of positively HT genes against which we can evaluate the candidate set, we devised another procedure using intergenetic alignments. Comparison with the benchmark set demonstrated the absence of a significant number of false positives in the candidate set, showing the high reliability of the method. Analyses of Escherichia coli K-12 operons revealed that although approximately 20 operons were probably gained from the last common ancestor of the 8 gamma-proteobacteria, deletion of intervening genes accounts for the formation of no operons, whereas horizontal transfer expanded 2 operons and introduced 4 entire operons. Based on these observations and reasoning, we suggest that the main mechanism of operon gain is HGT rather than intragenomic rearrangements. We propose that genes with related essential functions tend to reside in conserved operons, whereas genes in nonconserved operons mostly confer slight advantage to the organisms and frequently undergo horizontal transfer and decay. HT genes constitute at least 5.5% of the genes in the 8 species and approximately 45% of which originate from other gamma-proteobacteria. Genes involved in viral functions and mobile and extrachromosomal element functions are HT more often than expected. This finding indicates frequent mediation of HGT by bacteriophages. On the other hand, not only informational genes (those involved in transcription, translation, and related processes) but also operational genes (those involved in housekeeping) are HT less frequently than expected.

Evolution of the syntrophic interaction between Desulfovibrio vulgaris and Methanosarcina barkeri: Involvement of an ancient horizontal gene transfer

The sulfate reducing bacteria Desulfovibrio vulgaris and the methanogenic archaea Methanosarcina barkeri can grow syntrophically on lactate. In this study, a set of three closely located genes, DVU2103, DVU2104, and DVU2108 of D. vulgaris, was found to be up-regulated 2- to 4-fold following the lifestyle shift from syntroph to sulfate reducer; moreover, none of the genes in this gene set were differentially regulated when comparing gene expression from various D. vulgaris pure culture experiments. Although exact function of this gene set is unknown, the results suggest that it may play roles related to the lifestyle change of D. vulgaris from syntroph to sulfate reducer. This hypothesis is further supported by phylogenomic analyses showing that homologies of this gene set were only narrowly present in several groups of bacteria, most of which are restricted to a syntrophic lifestyle, such as Pelobacter carbinolicus, Syntrophobacter fumaroxidans, Syntrophomonas wolfei, and Syntrophus aciditrophicus. Phylogenetic analysis showed that all three individual genes in the gene set tended to be clustered with their homologies from archaeal genera, and they were rooted on archaeal species in the phylogenetic trees, suggesting that they were horizontally transferred from archaeal methanogens. In addition, no significant bias in codon and amino acid usages was detected between these genes and the rest of the D. vulgaris genome, suggesting the gene transfer may have occurred early in the evolutionary history so that sufficient time has elapsed to allow an adaptation to the codon and amino acid usages of D. vulgaris. This report provides novel insights into the origin and evolution of bacterial genes linked to the lifestyle change of D. vulgaris from a syntrophic to a sulfate-reducing lifestyle.

The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Determining the influence of horizontal gene transfer (HGT) on phylogenomic analyses and the retrieval of a tree of life is relevant for our understanding of microbial genome evolution. It is particularly difficult to differentiate between phylogenetic incongruence due to noise and that resulting from HGT. We have performed a large-scale, detailed evolutionary analysis of the different phylogenetic signals present in the genomes of Xanthomonadales, a group of Proteobacteria. We show that the presence of phylogenetic noise is not an obstacle to infer past and present HGTs during their evolution. The scenario derived from this analysis and other recently published reports reflect the confounding effects on bacterial phylogenomics of past and present HGT. Although transfers between closely related species are difficult to detect in genome-scale phylogenetic analyses, past transfers to the ancestor of extant groups appear as conflicting signals that occasionally might make impossible to determine the evolutionary origin of the whole genome.

Optimal gene partition into operons correlates with gene functional order

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.

Carotenoid biosynthesis in cyanobacteria: structural and evolutionary scenarios based on comparative genomics

Carotenoids are widely distributed pigments in nature and their biosynthetic pathway has been extensively studied in various organisms. The recent access to the overwhelming amount genomic data of cyanobacteria has given birth to a novel approach called comparative genomics. The putative enzymes involved in the carotenoid biosynthesis among the cyanobacteria were determined by similarity-based tools. The reconstruction of biosynthetic pathway was based on the related enzymes. It is interesting to find that nearly all the cyanobacteria share quite similar pathway to synthesize β-carotene except for Gloeobacter violaceus PCC 7421. The enzymes, crtE-B-P-Qb-L, involved in the upstream pathway are more conserved than the subsequent ones (crtW-R). In addition, many carotenoid synthesis enzymes exhibit diversity in structure and function. Such examples in the families of ζ –carotene desaturase, lycopene cylases and carotene ketolases were described in this article. When we mapped these crt genes to the cyanobacterial genomes, the crt genes showed great structural variation among species. All of them are dispersed on the whole chromosome in contrast to the linear adjacent distribution of the crt gene cluster in other eubacteria. Moreover, in unicellular cyanobacteria, each step of the carotenogenic pathway is usually catalyzed by one gene product, whereas multiple ketolase genes are found in filamentous cyanobacteria. Such increased numbers of crt genes and their correlation to the ecological adaptation were carefully discussed.

The life-cycle of operons

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution may be driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysA and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although genes within operons are usually closely spaced because of a neutral bias toward deletion and because of selection against large overlaps, genes in highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution may be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.

Genome-wide detection and analysis of homologous recombination among sequenced strains of Escherichia coli

BACKGROUND

Comparisons of complete bacterial genomes reveal evidence of lateral transfer of DNA across otherwise clonally diverging lineages. Some lateral transfer events result in acquisition of novel genomic segments and are easily detected through genome comparison. Other more subtle lateral transfers involve homologous recombination events that result in substitution of alleles within conserved genomic regions. This type of event is observed infrequently among distantly related organisms. It is reported to be more common within species, but the frequency has been difficult to quantify since the sequences under comparison tend to have relatively few polymorphic sites.

RESULTS

Here we report a genome-wide assessment of homologous recombination among a collection of six complete Escherichia coli and Shigella flexneri genome sequences. We construct a whole-genome multiple alignment and identify clusters of polymorphic sites that exhibit atypical patterns of nucleotide substitution using a random walk-based method. The analysis reveals one large segment (approximately 100 kb) and 186 smaller clusters of single base pair differences that suggest lateral exchange between lineages. These clusters include portions of 10% of the 3,100 genes conserved in six genomes. Statistical analysis of the functional roles of these genes reveals that several classes of genes are over-represented, including those involved in recombination, transport and motility.

CONCLUSION

We demonstrate that intraspecific recombination in E. coli is much more common than previously appreciated and may show a bias for certain types of genes. The described method provides high-specificity, conservative inference of past recombination events.

Gene fusion/fission is a major contributor to evolution of multi-domain bacterial proteins

Most proteins comprise one or several domains. New domain architectures can be created by combining previously existing domains. The elementary events that create new domain architectures may be categorized into three classes, namely domain(s) insertion or deletion (indel), exchange and repetition. Using 'DomainTeam', a tool dedicated to the search for microsyntenies of domains, we quantified the relative contribution of these events. This tool allowed us to collect homologous bacterial genes encoding proteins that have obviously evolved by modular assembly of domains. We show that indels are the most frequent elementary events and that they occur in most cases at either the N- or C-terminus of the proteins. As revealed by the genomic neighbourhood/context of the corresponding genes, we show that a substantial number of these terminal indels are the consequence of gene fusions/fissions. We provide evidence showing that the contribution of gene fusion/fission to the evolution of multi-domain bacterial proteins is lower-bounded by 27% and upper-bounded by 64%. We conclude that gene fusion/fission is a major contributor to the evolution of multi-domain bacterial proteins.

Post-genomics of the model haloarchaeon Halobacterium sp. NRC-1

Halobacteriumsp. NRC-1 is an extremely halophilic archaeon that is easily cultured and genetically tractable. Since its genome sequence was completed in 2000, a combination of genetic, transcriptomic, proteomic, and bioinformatic approaches have provided insights into both its extremophilic lifestyle as well as fundamental cellular processes common to all life forms. Here, we review post-genomic research on this archaeon, including investigations of DNA replication and repair systems, phototrophic, anaerobic, and other physiological capabilities, acidity of the proteome for function at high salinity, and role of lateral gene transfer in its evolution.

Reconstruction and regulation of the central catabolic pathway in the thermophilic propionate-oxidizing syntroph Pelotomaculum thermopropionicum

Obligate anaerobic bacteria fermenting volatile fatty acids in syntrophic association with methanogenic archaea share the intermediate bottleneck step in organic-matter decomposition. These organisms (called syntrophs) are biologically significant in terms of their growth at the thermodynamic limit and are considered to be the ideal model to address bioenergetic concepts. We conducted genomic and proteomic analyses of the thermophilic propionate-oxidizing syntroph Pelotomaculum thermopropionicum to obtain the genetic basis for its central catabolic pathway. Draft sequencing and subsequent targeted gap closing identified all genes necessary for reconstructing its propionate-oxidizing pathway (i.e., methylmalonyl coenzyme A pathway). Characteristics of this pathway include the following. (i) The initial two steps are linked to later steps via transferases. (ii) Each of the last three steps can be catalyzed by two different types of enzymes. It was also revealed that many genes for the propionate-oxidizing pathway, except for those for propionate coenzyme A transferase and succinate dehydrogenase, were present in an operon-like cluster and accompanied by multiple promoter sequences and a putative gene for a transcriptional regulator. Proteomic analysis showed that enzymes in this pathway were up-regulated when grown on propionate; of these enzymes, regulation of fumarase was the most stringent. We discuss this tendency of expression regulation based on the genetic organization of the open reading frame cluster. Results suggest that fumarase is the central metabolic switch controlling the metabolic flow and energy conservation in this syntroph.

Orthologs, paralogs, and evolutionary genomics

Orthologs and paralogs are two fundamentally different types of homologous genes that evolved, respectively, by vertical descent from a single ancestral gene and by duplication. Orthology and paralogy are key concepts of evolutionary genomics. A clear distinction between orthologs and paralogs is critical for the construction of a robust evolutionary classification of genes and reliable functional annotation of newly sequenced genomes. Genome comparisons show that orthologous relationships with genes from taxonomically distant species can be established for the majority of the genes from each sequenced genome. This review examines in depth the definitions and subtypes of orthologs and paralogs, outlines the principal methodological approaches employed for identification of orthology and paralogy, and considers evolutionary and functional implications of these concepts.

Clustered genes related to sulfate respiration in uncultured prokaryotes support the theory of their concomitant horizontal transfer

The dissimilatory reduction of sulfate is an ancient metabolic process central to today's biogeochemical cycling of sulfur and carbon in marine sediments. Until now its polyphyletic distribution was most parsimoniously explained by multiple horizontal transfers of single genes rather than by a not-yet-identified "metabolic island." Here we provide evidence that the horizontal transfer of a gene cluster may indeed be responsible for the patchy distribution of sulfate-reducing prokaryotes (SRP) in the phylogenetic tree. We isolated three DNA fragments (32 to 41 kb) from uncultured, closely related SRP from DNA directly extracted from two distinct marine sediments. Fosmid ws39f7, and partially also fosmids ws7f8 and hr42c9, harbored a core set of essential genes for the dissimilatory reduction of sulfate, including enzymes for the reduction of sulfur intermediates and synthesis of the prosthetic group of the dissimilatory sulfite reductase. Genome comparisons suggest that encoded membrane proteins universally present among SRP are critical for electron transfer to cytoplasmic enzymes. In addition, novel, conserved hypothetical proteins that are likely involved in dissimilatory sulfate reduction were identified. Based on comparative genomics and previously published experimental evidence, a more comprehensive model of dissimilatory sulfate reduction is presented. The observed clustering of genes involved in dissimilatory sulfate reduction has not been previously found. These findings strongly support the hypothesis that genes responsible for dissimilatory sulfate reduction were concomitantly transferred in a single event among prokaryotes. The acquisition of an optimized gene set would enormously facilitate a successful implementation of a novel pathway.

The Evolution of Microbial Phosphonate Degradative Pathways

Phosphonate utilization by microbes provides a potential source of phosphorus for their growth. Homologous genes for both C-P lyase and phosphonatase degradative pathways are distributed in distantly related bacterial species. The phn gene clusters for the C-P lyase pathway show great structural and compositional variation among organisms, but all contain phnG-phnM genes that are essential for C-P bond cleavage. In the gamma-proteobacterium Erwinia carotovora, genes common to phosphonate biosyntheses were found in neighboring positions of those for the C-P lyase degradative pathway and in the same transcriptional direction. A gene encoding a hypothetical protein DUF1045 was found predominantly associated with the phn gene cluster and was predicted functionally related to C-P bond cleavage. Genes for phosphonate degradation are frequently located in close proximity of genes encoding transposases or other mobile elements. Phylogenetic analyses suggest that both degradative pathways have been subject to extensive lateral gene transfers during their evolution. The implications of plasmids and transposition in the evolution of phosphonate degradation are also discussed.

Comparative genomics of Thermus thermophilus and Deinococcus radiodurans: divergent routes of adaptation to thermophily and radiation resistance

Background

Thermus thermophilus and Deinococcus radiodurans belong to a distinct bacterial clade but have remarkably different phenotypes. T. thermophilus is a thermophile, which is relatively sensitive to ionizing radiation and desiccation, whereas D. radiodurans is a mesophile, which is highly radiation- and desiccation-resistant. Here we present an in-depth comparison of the genomes of these two related but differently adapted bacteria.

Results

By reconstructing the evolution of Thermus and Deinococcus after the divergence from their common ancestor, we demonstrate a high level of post-divergence gene flux in both lineages. Various aspects of the adaptation to high temperature in Thermus can be attributed to horizontal gene transfer from archaea and thermophilic bacteria; many of the horizontally transferred genes are located on the single megaplasmid of Thermus. In addition, the Thermus lineage has lost a set of genes that are still present in Deinococcus and many other mesophilic bacteria but are not common among thermophiles. By contrast, Deinococcus seems to have acquired numerous genes related to stress response systems from various bacteria. A comparison of the distribution of orthologous genes among the four partitions of the Deinococcus genome and the two partitions of the Thermus genome reveals homology between the Thermus megaplasmid (pTT27) and Deinococcus megaplasmid (DR177).

Conclusion

After the radiation from their common ancestor, the Thermus and Deinococcus lineages have taken divergent paths toward their distinct lifestyles. In addition to extensive gene loss, Thermus seems to have acquired numerous genes from thermophiles, which likely was the decisive contribution to its thermophilic adaptation. By contrast, Deinococcus lost few genes but seems to have acquired many bacterial genes that apparently enhanced its ability to survive different kinds of environmental stresses. Notwithstanding the accumulation of horizontally transferred genes, we also show that the single megaplasmid of Thermus and the DR177 megaplasmid of Deinococcus are homologous and probably were inherited from the common ancestor of these bacteria.

Assessing Protein Co-evolution in the Context of the Tree of Life Assists in the Prediction of the Interactome

The identification of the whole set of protein interactions taking place in an organism is one of the main tasks in genomics, proteomics and systems biology. One of the computational techniques used by many investigators for studying and predicting protein interactions is the comparison of evolutionary histories (phylogenetic trees), under the hypothesis that interacting proteins would be subject to a similar evolutionary pressure resulting in a similar topology of the corresponding trees. Here, we present a new approach to predict protein interactions from phylogenetic trees, which incorporates information on the overall evolutionary histories of the species (i.e. the canonical "tree of life") in order to correct by the expected background similarity due to the underlying speciation events. We test the new approach in the largest set of annotated interacting proteins for Escherichia coli. This assessment of co-evolution in the context of the tree of life leads to a highly significant improvement (P(N) by sign test approximately 10E-6) in predicting interaction partners with respect to the previous technique, which does not incorporate information on the overall speciation tree. For half of the proteins we found a real interactor among the 6.4% top scores, compared with the 16.5% by the previous method. We applied the new method to the whole E.coli proteome and propose functions for some hypothetical proteins based on their predicted interactors. The new approach allows us also to detect non-canonical evolutionary events, in particular horizontal gene transfers. We also show that taking into account these non-canonical evolutionary events when assessing the similarity between evolutionary trees improves the performance of the method predicting interactions.

Horizontal transfer of two operons coding for hydrogenases between bacteria and archaea

Using a phylogenetic approach, we discovered three putative horizontal transfers between bacterial and archaeal species involving large clusters of genes. One transfer involves an operon of 13 genes, called mbx, which probably was transferred into the genome of Thermotoga maritima from a species belonging or close to the Pyrococcus genus. The two others implied an operon of six genes, called ech, transferred independently to the genomes of Thermoanaerobacter tengcongensis and Desulfovibrio gigas, from a species belonging or close to the Methanosarcina genus. All these transfers affected operons coding for multisubunit membrane-bound (NiFe) hydrogenases involved in the energy metabolism of the donor genomes. The functionality of the transferred operons has not been experimentally demonstrated for T. maritima, whereas in D. gigas and T. tengcongensis the encoded multisubunit hydrogenase could have a role in energy conservation. This report adds several cases of horizontal gene transfers among hydrogenases already described.

Operon formation is driven by co-regulation and not by horizontal gene transfer

The organization of bacterial genes into operons was originally ascribed to the benefits of co-regulation. More recently, the "selfish operon" model, in which operons are formed by repeated gain and loss of genes, was proposed. Indeed, operons are often subject to horizontal gene transfer (HGT). On the other hand, non-HGT genes are particularly likely to be in operons. To clarify whether HGT is involved in operon formation, we identified recently formed operons in Escherichia coli K12. We show that genes that have homologs in distantly related bacteria but not in close relatives of E. coli--indicating HGT--form new operons at about the same rates as native genes. Furthermore, genes in new operons are no more likely than other genes to have phylogenetic trees that are inconsistent with the species tree. In contrast, essential genes and ubiquitous genes without paralogs--genes believed to undergo HGT rarely--often form new operons. We conclude that HGT is not a cause of operon formation but instead promotes the prevalence of pre-existing operons. To explain operon formation, we propose that new operons reduce the amount of regulatory information required to specify optimal expression patterns and infer that operons should be more likely to evolve than independent promoters when regulation is complex. Consistent with this hypothesis, operons have greater amounts of conserved regulatory sequences than do individually transcribed genes.

Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea

We report the complete sequence of the 4,274,642-bp genome of Haloarcula marismortui, a halophilic archaeal isolate from the Dead Sea. The genome is organized into nine circular replicons of varying G+C compositions ranging from 54% to 62%. Comparison of the genome architectures of Halobacterium sp. NRC-1 and H. marismortui suggests a common ancestor for the two organisms and a genome of significantly reduced size in the former. Both of these halophilic archaea use the same strategy of high surface negative charge of folded proteins as means to circumvent the salting-out phenomenon in a hypersaline cytoplasm. A multitiered annotation approach, including primary sequence similarities, protein family signatures, structure prediction, and a protein function association network, has assigned putative functions for at least 58% of the 4242 predicted proteins, a far larger number than is usually achieved in most newly sequenced microorganisms. Among these assigned functions were genes encoding six opsins, 19 MCP and/or HAMP domain signal transducers, and an unusually large number of environmental response regulators-nearly five times as many as those encoded in Halobacterium sp. NRC-1--suggesting H. marismortui is significantly more physiologically capable of exploiting diverse environments. In comparing the physiologies of the two halophilic archaea, in addition to the expected extensive similarity, we discovered several differences in their metabolic strategies and physiological responses such as distinct pathways for arginine breakdown in each halophile. Finally, as expected from the larger genome, H. marismortui encodes many more functions and seems to have fewer nutritional requirements for survival than does Halobacterium sp. NRC-1.