Gene duplication and biased functional retention of paralogs in bacterial genomes

Promoter propagation in prokaryotes

Transcriptional activation or ‘rewiring’ of silent genes is an important, yet poorly understood, phenomenon in prokaryotic genomes. Anecdotal evidence coming from experimental evolution studies in bacterial systems has shown the promptness of adaptation upon appropriate selective pressure. In many cases, a partial or complete promoter is mobilized to silent genes from elsewhere in the genome. We term hereafter such recruited regulatory sequences as Putative Mobile Promoters (PMPs) and we hypothesize they have a large impact on rapid adaptation of novel or cryptic functions. Querying all publicly available prokaryotic genomes (1362) uncovered >4000 families of highly conserved PMPs (50 to 100 long with ≥80% nt identity) in 1043 genomes from 424 different genera. The genomes with the largest number of PMP families are Anabaena variabilis (28 families), Geobacter uraniireducens (27 families) and Cyanothece PCC7424 (25 families). Family size varied from 2 to 93 homologous promoters (in Desulfurivibrio alkaliphilus). Some PMPs are present in particular species, but some are conserved across distant genera. The identified PMPs represent a conservative dataset of very recent or conserved events of mobilization of non-coding DNA and thus they constitute evidence of an extensive reservoir of recyclable regulatory sequences for rapid transcriptional rewiring.

Rhizome of life, catastrophes, sequence exchanges, gene creations, and giant viruses: how microbial genomics challenges Darwin

Darwin's theory about the evolution of species has been the object of considerable dispute. In this review, we have described seven key principles in Darwin's book The Origin of Species and tried to present how genomics challenge each of these concepts and improve our knowledge about evolution. Darwin believed that species evolution consists on a positive directional selection ensuring the “survival of the fittest.” The most developed state of the species is characterized by increasing complexity. Darwin proposed the theory of “descent with modification” according to which all species evolve from a single common ancestor through a gradual process of small modification of their vertical inheritance. Finally, the process of evolution can be depicted in the form of a tree. However, microbial genomics showed that evolution is better described as the “biological changes over time.” The mode of change is not unidirectional and does not necessarily favors advantageous mutations to increase fitness it is rather subject to random selection as a result of catastrophic stochastic processes. Complexity is not necessarily the completion of development: several complex organisms have gone extinct and many microbes including bacteria with intracellular lifestyle have streamlined highly effective genomes. Genomes evolve through large events of gene deletions, duplications, insertions, and genomes rearrangements rather than a gradual adaptative process. Genomes are dynamic and chimeric entities with gene repertoires that result from vertical and horizontal acquisitions as well as de novo gene creation. The chimeric character of microbial genomes excludes the possibility of finding a single common ancestor for all the genes recorded currently. Genomes are collections of genes with different evolutionary histories that cannot be represented by a single tree of life (TOL). A forest, a network or a rhizome of life may be more accurate to represent evolutionary relationships among species.

Meta-analysis of general bacterial subclades in whole-genome phylogenies using tree topology profiling

In the last two decades, a large number of whole-genome phylogenies have been inferred to reconstruct the Tree of Life (ToL). Underlying data models range from gene or functionality content in species to phylogenetic gene family trees and multiple sequence alignments of concatenated protein sequences. Diversity in data models together with the use of different tree reconstruction techniques, disruptive biological effects and the steadily increasing number of genomes have led to a huge diversity in published phylogenies. Comparison of those and, moreover, identification of the impact of inference properties (underlying data model, inference technique) on particular reconstructions is almost impossible. In this work, we introduce tree topology profiling as a method to compare already published whole-genome phylogenies. This method requires visual determination of the particular topology in a drawn whole-genome phylogeny for a set of particular bacterial clans. For each clan, neighborhoods to other bacteria are collected into a catalogue of generalized alternative topologies. Particular topology alternatives found for an ordered list of bacterial clans reveal a topology profile that represents the analyzed phylogeny. To simulate the inhomogeneity of published gene content phylogenies we generate a set of seven phylogenies using different inference techniques and the SYSTERS-PhyloMatrix data model. After tree topology profiling on in total 54 selected published and newly inferred phylogenies, we separate artefactual from biologically meaningful phylogenies and associate particular inference results (phylogenies) with inference background (inference techniques as well as data models). Topological relationships of particular bacterial species groups are presented. With this work we introduce tree topology profiling into the scientific field of comparative phylogenomics.

Protein-SIP enables time-resolved analysis of the carbon flux in a sulfate-reducing, benzene-degrading microbial consortium

Benzene is a major contaminant in various environments, but the mechanisms behind its biodegradation under strictly anoxic conditions are not yet entirely clear. Here we analyzed a benzene-degrading, sulfate-reducing enrichment culture originating from a benzene-contaminated aquifer by a metagenome-based functional metaproteomic approach, using protein-based stable isotope probing (protein-SIP). The time-resolved, quantitative analysis of carbon fluxes within the community supplied with either (13)C-labeled benzene or (13)C-labeled carbonate yielded different functional groups of organisms, with their peptides showing specific time dependencies of (13)C relative isotope abundance indicating different carbon utilization. Through a detailed analysis of the mass spectrometric (MS) data, it was possible to quantify the utilization of the initial carbon source and the metabolic intermediates. The functional groups were affiliated to Clostridiales, Deltaproteobacteria and Bacteroidetes/Chlorobi. The Clostridiales-related organisms were involved in benzene degradation, putatively by fermentation, and additionally used significant amounts of carbonate as a carbon source. The other groups of organisms were found to perform diverse functions, with Deltaproteobacteria degrading fermentation products and Bacteroidetes/Chlorobi being putative scavengers feeding on dead cells. A functional classification of identified proteins supported this allocation and gave further insights into the metabolic pathways and the interactions between the community members. This example shows how protein-SIP can be applied to obtain temporal and phylogenetic information about functional interdependencies within microbial communities.

Genome plasticity and systems evolution in Streptomyces

Background

Streptomycetes are filamentous soil-dwelling bacteria. They are best known as the producers of a great variety of natural products such as antibiotics, antifungals, antiparasitics, and anticancer agents and the decomposers of organic substances for carbon recycling. They are also model organisms for the studies of gene regulatory networks, morphological differentiation, and stress response. The availability of sets of genomes from closely related Streptomyces strains makes it possible to assess the mechanisms underlying genome plasticity and systems adaptation.

Results

We present the results of a comprehensive analysis of the genomes of five Streptomyces species with distinct phenotypes. These streptomycetes have a pan-genome comprised of 17,362 orthologous families which includes 3,096 components in the core genome, 5,066 components in the dispensable genome, and 9,200 components that are uniquely present in only one species. The core genome makes up about 33%-45% of each genome repertoire. It contains important genes for Streptomyces biology including those involved in gene regulation, secretion, secondary metabolism and morphological differentiation. Abundant duplicate genes have been identified, with 4%-11% of the whole genomes composed of lineage-specific expansions (LSEs), suggesting that frequent gene duplication or lateral gene transfer events play a role in shaping the genome diversification within this genus. Two patterns of expansion, single gene expansion and chromosome block expansion are observed, representing different scales of duplication.

Conclusions

Our results provide a catalog of genome components and their potential functional roles in gene regulatory networks and metabolic networks. The core genome components reveal the minimum requirement for streptomycetes to sustain a successful lifecycle in the soil environment, reflecting the effects of both genome evolution and environmental stress acting upon the expressed phenotypes. A better understanding of the LSE gene families will, on the other hand, bring a wealth of new insights into the mechanisms underlying strain-specific phenotypes, such as the production of novel antibiotics, pathogenesis, and adaptive response to environmental challenges.

Comparative genomics reveals adaptation by Alteromonas sp. SN2 to marine tidal-flat conditions: cold tolerance and aromatic hydrocarbon metabolism

Alteromonas species are globally distributed copiotrophic bacteria in marine habitats. Among these, sea-tidal flats are distinctive: undergoing seasonal temperature and oxygen-tension changes, plus periodic exposure to petroleum hydrocarbons. Strain SN2 of the genus Alteromonas was isolated from hydrocarbon-contaminated sea-tidal flat sediment and has been shown to metabolize aromatic hydrocarbons there. Strain SN2's genomic features were analyzed bioinformatically and compared to those of Alteromonas macleodii ecotypes: AltDE and ATCC 27126. Strain SN2's genome differs from that of the other two strains in: size, average nucleotide identity value, tRNA genes, noncoding RNAs, dioxygenase gene content, signal transduction genes, and the degree to which genes collected during the Global Ocean Sampling project are represented. Patterns in genetic characteristics (e.g., GC content, GC skew, Karlin signature, CRISPR gene homology) indicate that strain SN2's genome architecture has been altered via horizontal gene transfer (HGT). Experiments proved that strain SN2 was far more cold tolerant, especially at 5°C, than the other two strains. Consistent with the HGT hypothesis, a total of 15 genomic islands in strain SN2 likely confer ecological fitness traits (especially membrane transport, aromatic hydrocarbon metabolism, and fatty acid biosynthesis) specific to the adaptation of strain SN2 to its seasonally cold sea-tidal flat habitat.

A cost-effective and universal strategy for complete prokaryotic genomic sequencing proposed by computer simulation

BACKGROUND

Pyrosequencing techniques allow scientists to perform prokaryotic genome sequencing to achieve the draft genomic sequences within a few days. However, the assemblies with shotgun sequencing are usually composed of hundreds of contigs. A further multiplex PCR procedure is needed to fill all the gaps and link contigs into complete chromosomal sequence, which is the basis for prokaryotic comparative genomic studies. In this article, we study various pyrosequencing strategies by simulated assembling from 100 prokaryotic genomes.

FINDINGS

Simulation study shows that a single end 454 Jr. run combined with a paired end 454 Jr. run (8 kb library) can produce: 1) ~90% of 100 assemblies with < 10 scaffolds and ~95% of 100 assemblies with < 150 contigs; 2) average contig N50 size is over 331 kb; 3) average single base accuracy is > 99.99%; 4) average false gene duplication rate is < 0.7%; 5) average false gene loss rate is < 0.4%.

CONCLUSIONS

A single end 454 Jr. run combined with a paired end 454 Jr. run (8 kb library) is a cost-effective way for prokaryotic whole genome sequencing. This strategy provides solution to produce high quality draft assemblies for most of prokaryotic organisms within days. Due to the small number of assembled scaffolds, the following multiplex PCR procedure (for gap filling) would be easy. As a result, large scale prokaryotic whole genome sequencing projects may be finished within weeks.

Metabolically versatile large-genome prokaryotes

Although versatile microorganisms are critical in industrial applications where the ability to cope with change and carry out complex tasks is needed, very little is in fact known about the evolutionary and ecological meanings of versatility in prokaryotes. Testing the hypothesis that a large genome size is a prerequisite for versatility in prokaryotes, we found that putatively versatile prokaryotes are phylogenetically and ecologically diverse and indeed include many well known and commercially relevant versatile microorganisms. Despite individual differences in metabolic abilities, a common trait of large-genome prokaryotes appears that they have gained their large genomes as an evolutionary response to nutrient-scarce and/or variable environments. This insight seriously questions the ability of traditional microbiology methods to isolate versatile prokaryotes and casts doubt on the ecological relevance of knowledge based on the study of specialists.

Prokaryotic systematics in the genomics era

As an essential and basic biological discipline, prokaryotic systematics is entering the era of genomics. This paradigmatic shift is significant not only for understanding molecular phylogeny at the whole genome level but also in revealing the genetic or epigenetic basis that accounts for the phenotypic criteria used to classify and identify species. These developments provide an opportunity and a challenge for systematists to reanalyze the molecular mechanisms underlying the taxonomic characteristics of prokaryotes by drawing the knowledge from studies of genomics and/or functional genomics employing platform technologies and related bioinformatics tools. It is expected that taxonomic books, such as Bergey's Manual of Systematic Bacteriology may evolve into a systematics library indexed by phylogenomic information with an comprehensive understanding of prokaryotic speciation and associated increasing knowledge of biological phenomena.

Heterologous expression and characterization of two 1-hydroxy-2-naphthoic acid dioxygenases from Arthrobacter phenanthrenivorans

A protein fraction exhibiting 1-hydroxy-2-naphthoic acid (1-H2NA) dioxygenase activity was purified via ion exchange, hydrophobic interactions, and gel filtration chromatography from Arthrobacter phenanthrenivorans sp. nov. strain Sphe3 isolated from a Greek creosote-oil-polluted site. Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) and tandem MS (MS-MS) analysis revealed that the amino acid sequences of oligopeptides of the major 45-kDa protein species, as analyzed by SDS-PAGE and silver staining, comprising 29% of the whole sequence, exhibited strong homology with 1-H2NA dioxygenase of Nocardioides sp. strain KP7. A BLAST search of the recently sequenced Sphe3 genome revealed two putative open reading frames, named diox1 and diox2, showing 90% nucleotide identity to each other and 85% identity at the amino acid level with the Nocardia sp. homologue. diox1 was found on an indigenous Sphe3 plasmid, whereas diox2 was located on the chromosome. Both genes were induced by the presence of phenanthrene used as a sole carbon and energy source, and as expected, both were subject to carbon catabolite repression. The relative RNA transcription level of the chromosomal (diox2) gene was significantly higher than that of its plasmid (diox1) homologue. Both diox1 and diox2 putative genes were PCR amplified, cloned, and overexpressed in Escherichia coli. Recombinant E. coli cells expressed 1-H2NA dioxygenase activity. Recombinant enzymes exhibited Michaelis-Menten kinetics with an apparent K(m) of 35 μM for Diox1 and 29 μM for Diox2, whereas they showed similar kinetic turnover characteristics with K(cat)/K(m) values of 11 × 10(6) M(-1) s(-1) and 12 × 10(6) M(-1) s(-1), respectively. Occurrence of two diox1 and diox2 homologues in the Sphe3 genome implies that a replicative transposition event has contributed to the evolution of 1-H2NA dioxygenase in A. phenanthrenivorans.

Gentisate 1,2-dioxygenase, in the third naphthalene catabolic gene cluster of Polaromonas naphthalenivorans CJ2, has a role in naphthalene degradation

Polaromonas naphthalenivorans strain CJ2 metabolizes naphthalene via the gentisate pathway and has recently been shown to carry a third copy of gentisate 1,2-dioxygenase (GDO), encoded by nagI3, within a previously uncharacterized naphthalene catabolic gene cluster. The role of this cluster (especially nagI3) in naphthalene metabolism of strain CJ2 was investigated by documenting patterns in regulation, transcription and enzyme activity. Transcriptional analysis of wild-type cells showed the third cluster to be polycistronic and that nagI3 was expressed at a relatively high level. Individual knockout mutants of all three nagI genes were constructed and their influence on both GDO activity and cell growth was evaluated. Of the three knockout strains, CJ2ΔnagI3 showed severely diminished GDO activity and grew slowest on aromatic substrates. These observations are consistent with the hypothesis that nagI3 may prevent toxic intracellular levels of gentisate from accumulating in CJ2 cells. All three nagI genes from strain CJ2 were cloned into Escherichia coli: the nagI2 and nagI3 genes were successfully overexpressed. The subunit mass of the GDOs were ~36–39 kDa, and their structures were deduced to be dimeric. The K m values of NagI2 and NagI3 were 31 and 10 µM, respectively, indicating that the higher affinity of NagI3 for gentisate may protect the wild-type cells from gentisate toxicity. These results provide clues for explaining why the third gene cluster, particularly the nagI3 gene, is important in strain CJ2. The organization of genes in the third gene cluster matched that of clusters in Polaromonas sp. JS666 and Leptothrix cholodnii SP-6. While horizontal gene transfer (HGT) is one hypothesis for explaining this genetic motif, gene duplication within the ancestral lineage is equally valid. The HGT hypothesis was discounted by noting that the nagI3 allele of strain CJ2 did not share high sequence identity with its homologues in Polaromonas sp. JS666 and L. cholodnii SP-6.

Biological consequences of ancient gene acquisition and duplication in the large genome of Candidatus Solibacter usitatus Ellin6076

Members of the bacterial phylum Acidobacteria are widespread in soils and sediments worldwide, and are abundant in many soils. Acidobacteria are challenging to culture in vitro, and many basic features of their biology and functional roles in the soil have not been determined. Candidatus Solibacter usitatus strain Ellin6076 has a 9.9 Mb genome that is approximately 2–5 times as large as the other sequenced Acidobacteria genomes. Bacterial genome sizes typically range from 0.5 to 10 Mb and are influenced by gene duplication, horizontal gene transfer, gene loss and other evolutionary processes. Our comparative genome analyses indicate that the Ellin6076 large genome has arisen by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, and widespread small-scale gene duplications, resulting in an increased number of paralogs. Low amino acid sequence identities among functional group members, and lack of conserved gene order and orientation in regions containing similar groups of paralogs, suggest that most of the paralogs are not the result of recent duplication events. The genome sizes of additional cultured Acidobacteria strains were estimated using pulsed-field gel electrophoresis to determine the prevalence of the large genome trait within the phylum. Members of subdivision 3 had larger genomes than those of subdivision 1, but none were as large as the Ellin6076 genome. The large genome of Ellin6076 may not be typical of the phylum, and encodes traits that could provide a selective metabolic, defensive and regulatory advantage in the soil environment.

A computational study of elongation factor G (EFG) duplicated genes: diverged nature underlying the innovation on the same structural template

Background

Elongation factor G (EFG) is a core translational protein that catalyzes the elongation and recycling phases of translation. A more complex picture of EFG's evolution and function than previously accepted is emerging from analyzes of heterogeneous EFG family members. Whereas the gene duplication is postulated to be a prominent factor creating functional novelty, the striking divergence between EFG paralogs can be interpreted in terms of innovation in gene function.

Methodology/Principal Findings

We present a computational study of the EFG protein family to cover the role of gene duplication in the evolution of protein function. Using phylogenetic methods, genome context conservation and insertion/deletion (indel) analysis we demonstrate that the EFG gene copies form four subfamilies: EFG I, spdEFG1, spdEFG2, and EFG II. These ancient gene families differ by their indispensability, degree of divergence and number of indels. We show the distribution of EFG subfamilies and describe evidences for lateral gene transfer and recent duplications. Extended studies of the EFG II subfamily concern its diverged nature. Remarkably, EFG II appears to be a widely distributed and a much-diversified subfamily whose subdivisions correlate with phylum or class borders. The EFG II subfamily specific characteristics are low conservation of the GTPase domain, domains II and III; absence of the trGTPase specific G2 consensus motif “RGITI”; and twelve conserved positions common to the whole subfamily. The EFG II specific functional changes could be related to changes in the properties of nucleotide binding and hydrolysis and strengthened ionic interactions between EFG II and the ribosome, particularly between parts of the decoding site and loop I of domain IV.

Conclusions/Significance

Our work, for the first time, comprehensively identifies and describes EFG subfamilies and improves our understanding of the function and evolution of EFG duplicated genes.

Dynamics of gene duplication in the genomes of chlorophyll d-producing cyanobacteria: implications for the ecological niche

Gene duplication may be an important mechanism for the evolution of new functions and for the adaptive modulation of gene expression via dosage effects. Here, we analyzed the fate of gene duplicates for two strains of a novel group of cyanobacteria (genus Acaryochloris) that produces the far-red light absorbing chlorophyll d as its main photosynthetic pigment. The genomes of both strains contain an unusually high number of gene duplicates for bacteria. As has been observed for eukaryotic genomes, we find that the demography of gene duplicates can be well modeled by a birth-death process. Most duplicated Acaryochloris genes are of comparatively recent origin, are strain-specific, and tend to be located on different genetic elements. Analyses of selection on duplicates of different divergence classes suggest that a minority of paralogs exhibit near neutral evolutionary dynamics immediately following duplication but that most duplicate pairs (including those which have been retained for long periods) are under strong purifying selection against amino acid change. The likelihood of duplicate retention varied among gene functional classes, and the pronounced differences between strains in the pool of retained recent duplicates likely reflects differences in the nutrient status and other characteristics of their respective environments. We conclude that most duplicates are quickly purged from Acaryochloris genomes and that those which are retained likely make important contributions to organism ecology by conferring fitness benefits via gene dosage effects. The mechanism of enhanced duplication may involve homologous recombination between genetic elements mediated by paralogous copies of recA.

Simultaneous identification of duplications and lateral gene transfers

The incongruency between a gene tree and a corresponding species tree can be attributed to evolutionary events such as gene duplication and gene loss. This paper describes a combinatorial model where so-called DTL-scenarios are used to explain the differences between a gene tree and a corresponding species tree taking into account gene duplications, gene losses, and lateral gene transfers (also known as horizontal gene transfers). The reasonable biological constraint that a lateral gene transfer may only occur between contemporary species leads to the notion of acyclic DTL-scenarios. Parsimony methods are introduced by defining appropriate optimization problems. We show that finding most parsimonious acyclic DTL-scenarios is NP-hard. However, by dropping the condition of acyclicity, the problem becomes tractable, and we provide a dynamic programming algorithm as well as a fixed-parameter tractable algorithm for finding most parsimonious DTL-scenarios.

Horizontal transfer, not duplication, drives the expansion of protein families in prokaryotes

Gene duplication followed by neo- or sub-functionalization deeply impacts the evolution of protein families and is regarded as the main source of adaptive functional novelty in eukaryotes. While there is ample evidence of adaptive gene duplication in prokaryotes, it is not clear whether duplication outweighs the contribution of horizontal gene transfer in the expansion of protein families. We analyzed closely related prokaryote strains or species with small genomes (Helicobacter, Neisseria, Streptococcus, Sulfolobus), average-sized genomes (Bacillus, Enterobacteriaceae), and large genomes (Pseudomonas, Bradyrhizobiaceae) to untangle the effects of duplication and horizontal transfer. After removing the effects of transposable elements and phages, we show that the vast majority of expansions of protein families are due to transfer, even among large genomes. Transferred genes--xenologs--persist longer in prokaryotic lineages possibly due to a higher/longer adaptive role. On the other hand, duplicated genes--paralogs--are expressed more, and, when persistent, they evolve slower. This suggests that gene transfer and gene duplication have very different roles in shaping the evolution of biological systems: transfer allows the acquisition of new functions and duplication leads to higher gene dosage. Accordingly, we show that paralogs share most protein-protein interactions and genetic regulators, whereas xenologs share very few of them. Prokaryotes invented most of life's biochemical diversity. Therefore, the study of the evolution of biology systems should explicitly account for the predominant role of horizontal gene transfer in the diversification of protein families.

Adaptation of cells to new environments

The evolutionary success of an organism is a testament to its inherent capacity to keep pace with environmental conditions that change over short and long periods. Mechanisms underlying adaptive processes are being investigated with renewed interest and excitement. This revival is partly fueled by powerful technologies that can probe molecular phenomena at a systems scale. Such studies provide spectacular insight into the mechanisms of adaptation, including rewiring of regulatory networks via natural selection of horizontal gene transfers, gene duplication, deletion, readjustment of kinetic parameters, and myriad other genetic reorganizational events. Here, we discuss advances in prokaryotic systems biology from the perspective of evolutionary principles that have shaped regulatory networks for dynamic adaptation to environmental change.

The prevalence of gene duplications and their ancient origin in Rhodobacter sphaeroides 2.4.1

BACKGROUND

Rhodobacter sphaeroides 2.4.1 is a metabolically versatile organism that belongs to α-3 subdivision of Proteobacteria. The present study was to identify the extent, history, and role of gene duplications in R. sphaeroides 2.4.1, an organism that possesses two chromosomes.

RESULTS

A protein similarity search (BLASTP) identified 1247 orfs (~29.4% of the total protein coding orfs) that are present in 2 or more copies, 37.5% (234 gene-pairs) of which exist in duplicate copies. The distribution of the duplicate gene-pairs in all Clusters of Orthologous Groups (COGs) differed significantly when compared to the COG distribution across the whole genome. Location plots revealed clusters of gene duplications that possessed the same COG classification. Phylogenetic analyses were performed to determine a tree topology predicting either a Type-A or Type-B phylogenetic relationship. A Type-A phylogenetic relationship shows that a copy of the protein-pair matches more with an ortholog from a species closely related to R. sphaeroides while a Type-B relationship predicts the highest match between both copies of the R. sphaeroides protein-pair. The results revealed that ~77% of the proteins exhibited a Type-A phylogenetic relationship demonstrating the ancient origin of these gene duplications. Additional analyses on three other strains of R. sphaeroides revealed varying levels of gene loss and retention in these strains. Also, analyses on common gene pairs among the four strains revealed that these genes experience similar functional constraints and undergo purifying selection.

CONCLUSIONS

Although the results suggest that the level of gene duplication in organisms with complex genome structuring (more than one chromosome) seems to be not markedly different from that in organisms with only a single chromosome, these duplications may have aided in genome reorganization in this group of eubacteria prior to the formation of R. sphaeroides as gene duplications involved in specialized functions might have contributed to complex genomic development.

Gene duplications in prokaryotes can be associated with environmental adaptation

Background

Gene duplication is a normal evolutionary process. If there is no selective advantage in keeping the duplicated gene, it is usually reduced to a pseudogene and disappears from the genome. However, some paralogs are retained. These gene products are likely to be beneficial to the organism, e.g. in adaptation to new environmental conditions. The aim of our analysis is to investigate the properties of paralog-forming genes in prokaryotes, and to analyse the role of these retained paralogs by relating gene properties to life style of the corresponding prokaryotes.

Results

Paralogs were identified in a number of prokaryotes, and these paralogs were compared to singletons of persistent orthologs based on functional classification. This showed that the paralogs were associated with for example energy production, cell motility, ion transport, and defence mechanisms. A statistical overrepresentation analysis of gene and protein annotations was based on paralogs of the 200 prokaryotes with the highest fraction of paralog-forming genes. Biclustering of overrepresented gene ontology terms versus species was used to identify clusters of properties associated with clusters of species. The clusters were classified using similarity scores on properties and species to identify interesting clusters, and a subset of clusters were analysed by comparison to literature data. This analysis showed that paralogs often are associated with properties that are important for survival and proliferation of the specific organisms. This includes processes like ion transport, locomotion, chemotaxis and photosynthesis. However, the analysis also showed that the gene ontology terms sometimes were too general, imprecise or even misleading for automatic analysis.

Conclusions

Properties described by gene ontology terms identified in the overrepresentation analysis are often consistent with individual prokaryote lifestyles and are likely to give a competitive advantage to the organism. Paralogs and singletons dominate different categories of functional classification, where paralogs in particular seem to be associated with processes involving interaction with the environment.

Impact of multiple beta-ketothiolase deletion mutations in Ralstonia eutropha H16 on the composition of 3-mercaptopropionic acid-containing copolymers

beta-Ketothiolases catalyze the first step of poly(3-hydroxybutyrate) [poly(3HB)] synthesis in bacteria by condensing two molecules of acetyl coenzyme A (acetyl-CoA) to acetoacetyl-CoA. Analyses of the genome sequence of Ralstonia eutropha H16 revealed 15 isoenzymes of PhaA in this bacterium. In this study, we generated knockout mutants of various phaA homologues to investigate their role in and contributions to poly(3HB) metabolism and to suppress biosynthesis of 3HB-CoA for obtaining enhanced molar 3-mercaptopriopionate (3MP) contents in poly(3HB-co-3MP) copolymers when cells were grown on gluconate plus 3-mercaptopropionate or 3,3'-dithiodipropionate. In silico sequence analysis of PhaA homologues, transcriptome data, and other aspects recommended the homologues phaA, bktB, H16_A1713/H16_B1771, H16_A1528, H16_B1369, H16_B0381, and H16_A0170 for further analysis. Single- and multiple-deletion mutants were generated to investigate the influence of these beta-ketothiolases on growth and polymer accumulation. The deletion of single genes resulted in no significant differences from the wild type regarding growth and polymer accumulation during cultivation on gluconate or gluconate plus 3MP. Deletion of phaA plus bktB (H16Delta2 mutant) resulted in approximately 30% less polymer accumulation than in the wild type. Deletion of H16_A1713/H16_B1771, H16_A1528, H16_B0381, and H16_B1369 in addition to phaA and bktB gave no differences in comparison to the H16Delta2 mutant. In contrast, deletion of H16_A0170 additionally to phaA and bktB yielded a mutant which accumulated about 30% poly(3HB) (wt/wt of the cell dry weight [CDW]). Although we were not able to suppress poly(3HB) biosynthesis completely, the copolymer compositions could be altered significantly with a lowered percentage ratio of 3HB constituents (from 85 to 52 mol%) and an increased percentage ratio of 3MP constituents (from 15 to 48 mol%), respectively. In this study, we demonstrated that PhaA, BktB, and H16_A0170 are majorly involved in poly(3HB) synthesis in R. eutropha H16. A fourth beta-ketothiolase or a combination of several of the other beta-ketothiolases contributed to a maximum of only 30% (wt/wt of CDW) of the remaining (co)polymer.

The complete genome sequence of Cupriavidus metallidurans strain CH34, a master survivalist in harsh and anthropogenic environments

Many bacteria in the environment have adapted to the presence of toxic heavy metals. Over the last 30 years, this heavy metal tolerance was the subject of extensive research. The bacterium Cupriavidus metallidurans strain CH34, originally isolated by us in 1976 from a metal processing factory, is considered a major model organism in this field because it withstands milli-molar range concentrations of over 20 different heavy metal ions. This tolerance is mostly achieved by rapid ion efflux but also by metal-complexation and -reduction. We present here the full genome sequence of strain CH34 and the manual annotation of all its genes. The genome of C. metallidurans CH34 is composed of two large circular chromosomes CHR1 and CHR2 of, respectively, 3,928,089 bp and 2,580,084 bp, and two megaplasmids pMOL28 and pMOL30 of, respectively, 171,459 bp and 233,720 bp in size. At least 25 loci for heavy-metal resistance (HMR) are distributed over the four replicons. Approximately 67% of the 6,717 coding sequences (CDSs) present in the CH34 genome could be assigned a putative function, and 9.1% (611 genes) appear to be unique to this strain. One out of five proteins is associated with either transport or transcription while the relay of environmental stimuli is governed by more than 600 signal transduction systems. The CH34 genome is most similar to the genomes of other Cupriavidus strains by correspondence between the respective CHR1 replicons but also displays similarity to the genomes of more distantly related species as a result of gene transfer and through the presence of large genomic islands. The presence of at least 57 IS elements and 19 transposons and the ability to take in and express foreign genes indicates a very dynamic and complex genome shaped by evolutionary forces. The genome data show that C. metallidurans CH34 is particularly well equipped to live in extreme conditions and anthropogenic environments that are rich in metals.

Evolution by leaps: gene duplication in bacteria

BACKGROUND

Sequence related families of genes and proteins are common in bacterial genomes. In Escherichia coli they constitute over half of the genome. The presence of families and superfamilies of proteins suggest a history of gene duplication and divergence during evolution. Genome encoded protein families, their size and functional composition, reflect metabolic potentials of the organisms they are found in. Comparing protein families of different organisms give insight into functional differences and similarities.

RESULTS

Equivalent enzyme families with metabolic functions were selected from the genomes of four experimentally characterized bacteria belonging to separate genera. Both similarities and differences were detected in the protein family memberships, with more similarities being detected among the more closely related organisms. Protein family memberships reflected known metabolic characteristics of the organisms. Differences in divergence of functionally characterized enzyme family members accounted for characteristics of taxa known to differ in those biochemical properties and capabilities. While some members of the gene families will have been acquired by lateral exchange and other former family members will have been lost over time, duplication and divergence of genes and functions appear to have been a significant contributor to the functional diversity of today's microbes.

CONCLUSIONS

Protein families seem likely to have arisen during evolution by gene duplication and divergence where the gene copies that have been retained are the variants that have led to distinct bacterial physiologies and taxa. Thus divergence of the duplicate enzymes has been a major process in the generation of different kinds of bacteria.

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.

Genuine genetic redundancy in maleylacetate-reductase-encoding genes involved in degradation of haloaromatic compounds by Cupriavidus necator JMP134

Maleylacetate reductases (MAR) are required for biodegradation of several substituted aromatic compounds. To date, the functionality of two MAR-encoding genes (tfdF I and tfdF II) has been reported in Cupriavidus necator JMP134(pJP4), a known degrader of aromatic compounds. These two genes are located in tfd gene clusters involved in the turnover of 2,4-dichlorophenoxyacetate (2,4-D) and 3-chlorobenzoate (3-CB). The C. necator JMP134 genome comprises at least three other genes that putatively encode MAR (tcpD, hqoD and hxqD), but confirmation of their functionality and their role in the catabolism of haloaromatic compounds has not been assessed. RT-PCR expression analyses of C. necator JMP134 cells exposed to 2,4-D, 3-CB, 2,4,6-trichlorophenol (2,4,6-TCP) or 4-fluorobenzoate (4-FB) showed that tfdF I and tfdF II are induced by haloaromatics channelled to halocatechols as intermediates. In contrast, 2,4,6-TCP only induces tcpD, and any haloaromatic compounds tested did not induce hxqD and hqoD. However, the tcpD, hxqD and hqoD gene products showed MAR activity in cell extracts and provided the MAR function for 2,4-D catabolism when heterologously expressed in MAR-lacking strains. Growth tests for mutants of the five MAR-encoding genes in strain JMP134 showed that none of these genes is essential for degradation of the tested compounds. However, the role of tfdF I/tfdF II and tcpD genes in the expression of MAR activity during catabolism of 2,4-D and 2,4,6-TCP, respectively, was confirmed by enzyme activity tests in mutants. These results reveal a striking example of genetic redundancy in the degradation of aromatic compounds.

Haloarchaeal assimilatory nitrate-reducing communities from a saline alkaline soil

Assimilatory nitrate reduction (ANR) is a pathway wherein NO(3)(-) is reduced to NH(4)(+), an N species that can be incorporated into the biomass. There is little information about the ANR genes in Archaea and most of the known information has been obtained from cultivable species. In this study, the diversity of the haloarchaeal assimilatory nitrate-reducing community was studied in an extreme saline alkaline soil of the former lake Texcoco (Mexico). Genes coding for the assimilatory nitrate reductase (narB) and the assimilatory nitrite reductase (nirA) were used as functional markers. Primers to amplify and detect partial narB and nirA were designed. The analysis of these amplicons by cloning and sequencing showed that the deduced protein fragments shared >45% identity with other NarB and NirA proteins from Euryarchaeota and <38% identity with other nitrate reductases from Bacteria and Crenarchaeota. Furthermore, these clone sequences were clustered within the class Halobacteria with strong support values in both constructed dendrograms, confirming that desired PCR products were obtained. The metabolic capacity to assimilate nitrate by these haloarchaea seems to be important given that at pH 10 and higher, NH(4)(+) is mostly converted to toxic and volatile NH(3), and NO(3)(-) becomes the preferable N source.

On the need for widespread horizontal gene transfers under genome size constraint

Background

While eukaryotes primarily evolve by duplication-divergence expansion (and reduction) of their own gene repertoire with only rare horizontal gene transfers, prokaryotes appear to evolve under both gene duplications and widespread horizontal gene transfers over long evolutionary time scales. But, the evolutionary origin of this striking difference in the importance of horizontal gene transfers remains by and large a mystery.

Hypothesis

We propose that the abundance of horizontal gene transfers in free-living prokaryotes is a simple but necessary consequence of two opposite effects: i) their apparent genome size constraint compared to typical eukaryote genomes and ii) their underlying genome expansion dynamics through gene duplication-divergence evolution, as demonstrated by the presence of many tandem and block repeated genes. In principle, this combination of genome size constraint and underlying duplication expansion should lead to a coalescent-like process with extensive turnover of functional genes. This would, however, imply the unlikely, systematic reinvention of functions from discarded genes within independent phylogenetic lineages. Instead, we propose that the long-term evolutionary adaptation of free-living prokaryotes must have resulted in the emergence of efficient non-phylogenetic pathways to circumvent gene loss.

Implications

This need for widespread horizontal gene transfers due to genome size constraint implies, in particular, that prokaryotes must remain under strong selection pressure in order to maintain the long-term evolutionary adaptation of their "mutualized" gene pool, beyond the inevitable turnover of individual prokaryote species. By contrast, the absence of genome size constraint for typical eukaryotes has presumably relaxed their need for widespread horizontal gene transfers and strong selection pressure. Yet, the resulting loss of genetic functions, due to weak selection pressure and inefficient gene recovery mechanisms, must have ultimately favored the emergence of more complex life styles and ecological integration of many eukaryotes.

Reviewers

This article was reviewed by Pierre Pontarotti, Eugene V Koonin and Sergei Maslov.

Comparative genomic analyses of Streptococcus mutans provide insights into chromosomal shuffling and species-specific content

Background

Streptococcus mutans is the major pathogen of dental caries, and it occasionally causes infective endocarditis. While the pathogenicity of this species is distinct from other human pathogenic streptococci, the species-specific evolution of the genus Streptococcus and its genomic diversity are poorly understood.

Results

We have sequenced the complete genome of S. mutans serotype c strain NN2025, and compared it with the genome of UA159. The NN2025 genome is composed of 2,013,587 bp, and the two strains show highly conserved core-genome. However, comparison of the two S. mutans strains showed a large genomic inversion across the replication axis producing an X-shaped symmetrical DNA dot plot. This phenomenon was also observed between other streptococcal species, indicating that streptococcal genetic rearrangements across the replication axis play an important role in Streptococcus genetic shuffling. We further confirmed the genomic diversity among 95 clinical isolates using long-PCR analysis. Genomic diversity in S. mutans appears to occur frequently between insertion sequence (IS) elements and transposons, and these diversity regions consist of restriction/modification systems, antimicrobial peptide synthesis systems, and transporters. S. mutans may preferentially reject the phage infection by clustered regularly interspaced short palindromic repeats (CRISPRs). In particular, the CRISPR-2 region, which is highly divergent between strains, in NN2025 has long repeated spacer sequences corresponding to the streptococcal phage genome.

Conclusion

These observations suggest that S. mutans strains evolve through chromosomal shuffling and that phage infection is not needed for gene acquisition. In contrast, S. pyogenes tolerates phage infection for acquisition of virulence determinants for niche adaptation.

Genomic and systems evolution in Vibrionaceae species

Background

The steadily increasing number of prokaryotic genomes has accelerated the study of genome evolution; in particular, the availability of sets of genomes from closely related bacteria has facilitated the exploration of the mechanisms underlying genome plasticity. The family Vibrionaceae is found in the Gammaproteobacteria and is abundant in aquatic environments. Taxa from the family Vibrionaceae are diversified in their life styles; some species are free living, others are symbiotic, and others are human pathogens. This diversity makes this family a useful set of model organisms for studying bacterial evolution. This evolution is driven by several forces, among them gene duplication and lateral gene transfer, which are believed to provide raw material for functional redundancy and novelty. The resultant gene copy increase in one genome is then detected as lineage-specific expansion (LSE).

Results

Here we present the results of a detailed comparison of the genomes of eleven Vibrionaceae strains that have distinct life styles and distinct phenotypes. The core genome shared by all eleven strains is composed of 1,882 genes, which make up about 31%–50% of the genome repertoire. We further investigated the distribution and features of genes that have been specifically expanded in one unique lineage of the eleven strains. Abundant duplicate genes have been identified in the eleven Vibrionaceae strains, with 1–11% of the whole genomes composed lineage specific radiations. These LSEs occurred in two distinct patterns: the first type yields one or more copies of a single gene; we call this a single gene expansion. The second pattern has a high evolutionary impact, as the expansion involves two or more gene copies in a block, with the duplicated block located next to the original block (a contiguous block expansion) or at some distance from the original block (a discontiguous block expansion). We showed that LSEs involve genes that are tied to defense and pathogenesis mechanisms as well as in the fundamental life cycle of Vibrionaceae species.

Conclusion

Our results provide evidence of genome plasticity and rapid evolution within the family Vibrionaceae. The comparisons point to sources of genomic variation and candidates for lineage-specific adaptations of each Vibrionaceae pathogen or nonpathogen strain. Such lineage specific expansions could reveal components in bacterial systems that, by their enhanced genetic variability, can be tied to responses to environmental challenges, interesting phenotypes, or adaptive pathogenic responses to host challenges.

Genesis, effects and fates of repeats in prokaryotic genomes

DNA repeats are causes and consequences of genome plasticity. Repeats are created by intrachromosomal recombination or horizontal transfer. They are targeted by recombination processes leading to amplifications, deletions and rearrangements of genetic material. The identification and analysis of repeats in nearly 700 genomes of bacteria and archaea is facilitated by the existence of sequence data and adequate bioinformatic tools. These have revealed the immense diversity of repeats in genomes, from those created by selfish elements to the ones used for protection against selfish elements, from those arising from transient gene amplifications to the ones leading to stable duplications. Experimental works have shown that some repeats do not carry any adaptive value, while others allow functional diversification and increased expression. All repeats carry some potential to disorganize and destabilize genomes. Because recombination and selection for repeats vary between genomes, the number and types of repeats are also quite diverse and in line with ecological variables, such as host-dependent associations or population sizes, and with genetic variables, such as the recombination machinery. From an evolutionary point of view, repeats represent both opportunities and problems. We describe how repeats are created and how they can be found in genomes. We then focus on the functional and genomic consequences of repeats that dictate their fate.

The life and death of gene families

One of the unique insights provided by the growing number of fully sequenced genomes is the pervasiveness of gene duplication and gene loss. Indeed, several metrics now suggest that rates of gene birth and death per gene are only 10-40% lower than nucleotide substitutions per site, and that per nucleotide, the consequent lineage-specific expansion and contraction of gene families may play at least as large a role in adaptation as changes in orthologous sequences. While gene family evolution is pervasive, it may be especially important in our own evolution since it appears that the "revolving door" of gene duplication and loss has undergone multiple accelerations in the lineage leading to humans. In this paper, we review current understanding of gene family evolution including: methods for inferring copy number change, evidence for adaptive expansion and adaptive contraction of gene families, the origins of new families and deaths of previously established ones, and finally we conclude with a perspective on challenges and promising directions for future research.

Origin and evolution of metabolic pathways

The emergence and evolution of metabolic pathways represented a crucial step in molecular and cellular evolution. In fact, the exhaustion of the prebiotic supply of amino acids and other compounds that were likely present in the ancestral environment, imposed an important selective pressure, favoring those primordial heterotrophic cells which became capable of synthesizing those molecules. Thus, the emergence of metabolic pathways allowed primitive organisms to become increasingly less-dependent on exogenous sources of organic compounds. Comparative analyses of genes and genomes from organisms belonging to Archaea, Bacteria and Eukarya revealed that, during evolution, different forces and molecular mechanisms might have driven the shaping of genomes and the arisal of new metabolic abilities. Among these gene elongations, gene and operon duplications undoubtedly played a major role since they can lead to the (immediate) appearance of new genetic material that, in turn, might undergo evolutionary divergence giving rise to new genes coding for new metabolic abilities. Gene duplication has been invoked in the different schemes proposed to explain why and how the extant metabolic pathways have arisen and shaped. Both the analysis of completely sequenced genomes and directed evolution experiments strongly support one of them, i.e. the patchwork hypothesis, according to which metabolic pathways have been assembled through the recruitment of primitive enzymes that could react with a wide range of chemically related substrates. However, the analysis of the structure and organization of genes belonging to ancient metabolic pathways, such as histidine biosynthesis and nitrogen fixation, suggested that other different hypothesis, i.e. the retrograde hypothesis or the semi-enzymatic theory, may account for the arisal of some metabolic routes.

Inferring functional relationships from conservation of gene order

The idea behind the gene neighbor method is that conservation of gene order in evolutionarily distant prokaryotes indicates functional association. The procedure presented here starts with the organization of all the genomes into pairs of adjacent genes. Then, pairs of genes in a genome of interest are mapped to their corresponding orthologs in other, informative, genomes. The final step is to determine whether the orthologs of each original pair of genes are also adjacent in the informative genome.

Multiple genome comparison within a bacterial species reveals a unit of evolution spanning two adjacent genes in a tandem paralog cluster

It has been assumed that an open reading frame (ORF) represents a unit of gene evolution as well as a unit of gene expression and function. In the present work, we report a case in which a unit comprising the 3' region of an ORF linked to a downstream intergenic region that is in turn linked to the 5' region of a downstream ORF has been conserved, and has served as the unit of gene evolution. The genes are tandem paralogous genes from the bacterium Staphylococcus aureus, for which more than ten entire genomes have been sequenced. We compared these multiple genome sequences at a locus for the lpl (lipoprotein-like) cluster (encoding lipoprotein homologs presumably related to their host interaction) in the genomic island termed nuSaalpha. A highly conserved nucleotide sequence found within every lpl ORF is likely to provide a site for homologous recombination. Comparison of phylogenies of the 5'-variable region and the 3'-variable region within the same ORF revealed significant incongruence. In contrast, pairs of the 3'-variable region of an ORF and the 5'-variable region of the next downstream ORF gave more congruent phylogenies, with distinct groups of conserved pairs. The intergenic region seemed to have coevolved with the flanking variable regions. Multiple recombination events at the central conserved region appear to have caused various types of rearrangements among strains, shuffling the two variable regions in one ORF, but maintaining a conserved unit comprising the 3'-variable region, the intergenic region, and the 5'-variable region spanning adjacent ORFs. This result has strong impact on our understanding of gene evolution because most gene lineages underwent tandem duplication and then diversified. This work also illustrates the use of multiple genome sequences for high-resolution evolutionary analysis within the same species.

The genome of the versatile nitrogen fixer Azorhizobium caulinodans ORS571

Background

Biological nitrogen fixation is a prokaryotic process that plays an essential role in the global nitrogen cycle. Azorhizobium caulinodans ORS571 has the dual capacity to fix nitrogen both as free-living organism and in a symbiotic interaction with Sesbania rostrata. The host is a fast-growing, submergence-tolerant tropical legume on which A. caulinodans can efficiently induce nodule formation on the root system and on adventitious rootlets located on the stem.

Results

The 5.37-Mb genome consists of a single circular chromosome with an overall average GC of 67% and numerous islands with varying GC contents. Most nodulation functions as well as a putative type-IV secretion system are found in a distinct symbiosis region. The genome contains a plethora of regulatory and transporter genes and many functions possibly involved in contacting a host. It potentially encodes 4717 proteins of which 96.3% have homologs and 3.7% are unique for A. caulinodans. Phylogenetic analyses show that the diazotroph Xanthobacter autotrophicus is the closest relative among the sequenced genomes, but the synteny between both genomes is very poor.

Conclusion

The genome analysis reveals that A. caulinodans is a diazotroph that acquired the capacity to nodulate most probably through horizontal gene transfer of a complex symbiosis island. The genome contains numerous genes that reflect a strong adaptive and metabolic potential. These combined features and the availability of the annotated genome make A. caulinodans an attractive organism to explore symbiotic biological nitrogen fixation beyond leguminous plants.

Adapting to environmental changes using specialized paralogs

When a bacterial species survives under changing environmental circumstances (e.g. salinity or temperature), its proteins might not function in all physicochemical conditions. We propose that prokaryotes cope with this problem by having two or more copies of the genes affected by environmental fluctuations, each one performing the same function under different conditions (i.e. ecoparalog). We identify potential examples in the bacterium Salinibacter ruber and in other species that experience wide environmental variations.

Molecular characterization of the Mycoplasma bovis p68 gene, encoding a basic membrane protein with homology to P48 of Mycoplasma agalactiae

Mycoplasmal lipoproteins are considered to be potential virulence determinants, which may carry out numerous important functions in pathogenesis including adhesion and immunomodulation. The prototype mycoplasmal immunomodulin is the macrophage-activating lipoprotein (MALP) of Mycoplasma fermentans. In this study, a homolog of the malp gene, designated p68, was identified and characterized in Mycoplasma bovis strain PG45 clonal variant #6. P68 belongs to the family of basic membrane proteins, which have been identified in diverse prokaryotes, including mycoplasmas. P68 revealed significant similarity and shared conserved selective lipoprotein-associated motifs with the highly immunogenic MALP-related lipoproteins P48 of M. bovis and P48 of Mycoplasma agalactiae. Determination of the genomic distribution of both M. bovis malp-homologs showed that p48 was present in all M. bovis strains tested, whereas the p68 gene was missing in some. Sequence comparison of the p68 genomic region in strains with and without this gene revealed that the region is very dynamic, with multiple genetic changes. Reverse-transcription PCR and primer extension analysis indicated that both p68 and p48 are transcribed in M. bovis under in vitro growth conditions. Mycoplasma bovis is the first mycoplasma species in which two malp-related genes have been identified.

Genome evolution and the emergence of fruiting body development in Myxococcus xanthus

Background

Lateral gene transfer (LGT) is thought to promote speciation in bacteria, though well-defined examples have not been put forward.

Methodology/Principle Findings

We examined the evolutionary history of the genes essential for a trait that defines a phylogenetic order, namely fruiting body development of the Myxococcales. Seventy-eight genes that are essential for Myxococcus xanthus development were examined for LGT. About 73% of the genes exhibit a phylogeny similar to that of the 16S rDNA gene and a codon bias consistent with other M. xanthus genes suggesting vertical transmission. About 22% have an altered codon bias and/or phylogeny suggestive of LGT. The remaining 5% are unique. Genes encoding signal production and sensory transduction were more likely to be transmitted vertically with clear examples of duplication and divergence into multigene families. Genes encoding metabolic enzymes were frequently acquired by LGT. Myxobacteria exhibit aerobic respiration unlike most of the δ Proteobacteria. M. xanthus contains a unique electron transport pathway shaped by LGT of genes for succinate dehydrogenase and three cytochrome oxidase complexes.

Conclusions/Significance

Fruiting body development depends on genes acquired by LGT, particularly those involved in polysaccharide production. We suggest that aerobic growth fostered innovation necessary for development by allowing myxobacteria access to a different gene pool from anaerobic members of the δ Proteobacteria. Habitat destruction and loss of species diversity could restrict the evolution of new bacterial groups by limiting the size of the prospective gene pool.

Double trouble: medical implications of genetic duplication and amplification in bacteria

Gene amplification allows organisms to adapt to changing environmental conditions. This type of increased gene dosage confers selectable benefits, typically by augmenting protein production. Gene amplification is a reversible process that does not require permanent genetic change. Although transient, altered gene dosage has significant medical impact. Recent examples of amplification in bacteria, described here, affect human disease by modifying antibiotic resistance, the virulence of pathogens, vaccine efficacy and antibiotic biosynthesis. Amplification is usually a two-step process whereby genetic duplication (step one) promotes further increases in copy number (step two). Both steps have important evolutionary significance for the emergence of innovative gene functions. Recent genome sequence analyses illustrate how genome plasticity can affect the evolution and immunogenic properties of bacterial pathogens.

A naturally occurring gene amplification leading to sulfonamide and trimethoprim resistance in Streptococcus agalactiae

Gene amplifications have been detected as a transitory phenomenon in bacterial cultures. They are predicted to contribute to rapid adaptation by simultaneously increasing the expression of genes clustered on the chromosome. However, genome amplifications have rarely been described in natural isolates. Through DNA array analysis, we have identified two Streptococcus agalactiae strains carrying tandem genome amplifications: a fourfold amplification of 13.5 kb and a duplication of 92 kb. Both amplifications were located close to the terminus of replication and originated independently from any long repeated sequence. They probably arose in the human host and showed different stabilities, the 13.5-kb amplification being lost at a frequency of 0.003 per generation and the 92-kb tandem duplication at a frequency of 0.035 per generation. The 13.5-kb tandem amplification carried the five genes required for dihydrofolate biosynthesis and led to both trimethoprim (TMP) and sulfonamide (SU) resistance. Resistance to SU probably resulted from the increased synthesis of dihydropteroate synthase, the target of this antibiotic, whereas the amplification of the whole pathway was responsible for TMP resistance. This revealed a new mechanism of resistance to TMP involving an increased dihydrofolate biosynthesis. This is, to our knowledge, the first reported case of naturally occurring antibiotic resistance resulting from genome amplification in bacteria. The low stability of DNA segment amplifications suggests that their role in antibiotic resistance might have been underestimated.

Secrets of soil survival revealed by the genome sequence of Arthrobacter aurescens TC1

Arthrobacter sp. strains are among the most frequently isolated, indigenous, aerobic bacterial genera found in soils. Member of the genus are metabolically and ecologically diverse and have the ability to survive in environmentally harsh conditions for extended periods of time. The genome of Arthrobacter aurescens strain TC1, which was originally isolated from soil at an atrazine spill site, is composed of a single 4,597,686 basepair (bp) circular chromosome and two circular plasmids, pTC1 and pTC2, which are 408,237 bp and 300,725 bp, respectively. Over 66% of the 4,702 open reading frames (ORFs) present in the TC1 genome could be assigned a putative function, and 13.2% (623 genes) appear to be unique to this bacterium, suggesting niche specialization. The genome of TC1 is most similar to that of Tropheryma, Leifsonia, Streptomyces, and Corynebacterium glutamicum, and analyses suggest that A. aurescens TC1 has expanded its metabolic abilities by relying on the duplication of catabolic genes and by funneling metabolic intermediates generated by plasmid-borne genes to chromosomally encoded pathways. The data presented here suggest that Arthrobacter's environmental prevalence may be due to its ability to survive under stressful conditions induced by starvation, ionizing radiation, oxygen radicals, and toxic chemicals.

Soil systems contain the greatest diversity of microorganisms on earth, with 5,000–10,000 species of microorganism per gram of soil. Arthrobacter sp. strains have a primitive life cycle and are among the most frequently isolated, indigenous soil bacteria, found in common and deep subsurface soils, arctic ice, and environments contaminated with industrial chemicals and radioactive materials. To better understand how these bacteria survive in environmentally harsh conditions, the authors used a structural genomics approach to identify genes involved in soil survival of Arthrobacter aurescens strain TC1, a bacterium originally isolated for its ability to degrade the herbicide atrazine. They found that the genome of this bacterium comprises a single circular chromosome and two plasmids that encode for a large number proteins involved in stress responses due to starvation, desiccation, oxygen radicals, and toxic chemicals. A. aurescens' metabolic versatility is in part due to the presence of duplicated catabolic genes and its ability to funnel plasmid-derived intermediates into chromosomally encoded pathways. Arthrobacter's array of genes that allow for survival in stressful conditions and its ability to produce a temperature-tolerant “cyst”-like resting cell render this soil microorganism able to survive and prosper in a variety of environmental conditions.

Natural merodiploidy of the lux-rib operon of Photobacterium leiognathi from coastal waters of Honshu, Japan

Sequence analysis of the bacterial luminescence (lux) genes has proven effective in helping resolve evolutionary relationships among luminous bacteria. Phylogenetic analysis using lux genes, however, is based on the assumptions that the lux genes are present as single copies on the bacterial chromosome and are vertically inherited. We report here that certain strains of Photobacterium leiognathi carry multiple phylogenetically distinct copies of the entire operon that codes for luminescence and riboflavin synthesis genes, luxCDABEG-ribEBHA. Merodiploid lux-rib strains of P. leiognathi were detected during sequence analysis of luxA. To define the gene content, organization, and sequence of each lux-rib operon, we constructed a fosmid library of genomic DNA from a representative merodiploid strain, lnuch.13.1. Sequence analysis of fosmid clones and genomic analysis of lnuch.13.1 defined two complete, physically separate, and apparently functional operons, designated lux-rib1 and lux-rib2. P. leiognathi strains lelon.2.1 and lnuch.21.1 were also found to carry lux-rib1 and lux-rib2, whereas ATCC 25521T apparently carries only lux-rib1. In lnuch.13.1, lelon.2.1, lnuch.21.1, and ATCC 25521T, lux-rib1 is flanked upstream by lumQ and putA and downstream by a gene for a hypothetical multidrug efflux pump. In contrast, transposase genes flank lux-rib2 of lnuch.13.1, and the chromosomal location of lux-rib2 apparently differs in lnuch.13.1, lelon.2.1, and lnuch.21.1. Phylogenetic analysis demonstrated that lux-rib1 and lux-rib2 are more closely related to each other than either one is to the lux and rib genes of other bacterial species, which rules out interspecies lateral gene transfer as the origin of lux-rib2 in P. leiognathi; lux-rib2 apparently arose within a previously unsampled or extinct P. leiognathi lineage. Analysis of 170 additional strains of P. leiognathi, for a total of 174 strains examined from coastal waters of Japan, Taiwan, the Philippine Islands, and Thailand, identified 106 strains that carry only a single lux-rib operon and 68 that carry multiple lux-rib operons. Strains bearing a single lux-rib operon were obtained throughout the geographic sampling range, whereas lux-rib merodiploid strains were found only in coastal waters of central Honshu. This is the first report of merodiploidy of lux or rib genes in a luminous bacterium and the first indication that a natural merodiploid state in bacteria can correlate with geography.

Evolution of ribonuclease H genes in prokaryotes to avoid inheritance of redundant genes

BACKGROUND

A theoretical model of genetic redundancy has proposed that the fates of redundant genes depend on the degree of functional redundancy, and that functionally redundant genes will not be inherited together. However, no example of actual gene evolution has been reported that can be used to test this model. Here, we analyzed the molecular evolution of the ribonuclease H (RNase H) family in prokaryotes and used the results to examine the implications of functional redundancy for gene evolution.

RESULTS

In prokaryotes, RNase H has been classified into RNase HI, HII, and HIII on the basis of amino acid sequences. Using 353 prokaryotic genomes, we identified the genes encoding the RNase H group and examined combinations of these genes in individual genomes. We found that the RNase H group may have evolved in such a way that the RNase HI and HIII genes will not coexist within a single genome--in other words, these genes are inherited in a mutually exclusive manner. Avoiding the simultaneous inheritance of the RNase HI and HIII genes is remarkable when RNase HI contains an additional non-RNase H domain, double-stranded RNA, and an RNA-DNA hybrid-binding domain, which is often observed in eukaryotic RNase H1. This evolutionary process may have resulted from functional redundancy of these genes, because the substrate preferences of RNase HI and RNase HIII are similar.

CONCLUSION

We provide two possible evolutionary models for RNase H genes in which functional redundancy contributes to the exclusion of redundant genes from the genome of a species. This is the first empirical study to show the effect of functional redundancy on changes in gene constitution during the course of evolution.

Gene conversion is a convergent strategy for pathogen antigenic variation

Recent studies on three unrelated vector-borne pathogens, Anaplasma marginale, Borrelia hermsii and Trypanosoma brucei, illustrate the central importance of gene conversion as a mechanism for antigenic variation, which results in subsequent evasion of the immune response and persistence in the reservoir host. The combination of genome sequence data and in vivo studies tracking variant emergence not only provides insight into the genetic mechanisms for variant generation and hierarchy in variant expression but also highlights gaps in our knowledge regarding variant capacity and usage in vivo.