Bacterial Genome Evolution Research Articles

Two fundamentally different classes of microbial genes.

The evolution of bacterial and archaeal genomes is highly dynamic and involves extensive horizontal gene transfer and gene loss1-4. Furthermore, many microbial species appear to have open pangenomes, where each newly sequenced genome contains more than 10% ORFans, that is, genes without detectable homologues in other species5,6. Here, we report a quantitative analysis of microbial genome evolution by fitting the parameters of a simple, steady-state evolutionary model to the comparative genomic data on the gene content and gene order similarity between archaeal genomes. The results reveal two sharply distinct classes of microbial genes, one of which is characterized by effectively instantaneous gene replacement, and the other consists of genes with finite, distributed replacement rates. These findings imply a conservative estimate of the size of the prokaryotic genomic universe, which appears to consist of at least a billion distinct genes. Furthermore, the same distribution of constraints is shown to govern the evolution of gene complement and gene order, without the need to invoke long-range conservation or the selfish operon concept7.

Lateral gene transfer, bacterial genome evolution, and the Anthropocene.

Lateral gene transfer (LGT) has significantly influenced bacterial evolution since the origins of life. It helped bacteria generate flexible, mosaic genomes and enables individual cells to rapidly acquire adaptive phenotypes. In turn, this allowed bacteria to mount strong defenses against human attempts to control their growth. The widespread dissemination of genes conferring resistance to antimicrobial agents has precipitated a crisis for modern medicine. Our actions can promote increased rates of LGT and also provide selective forces to fix such events in bacterial populations. For instance, the use of selective agents induces the bacterial SOS response, which stimulates LGT. We create hotspots for lateral transfer, such as wastewater systems, hospitals, and animal production facilities. Conduits of gene transfer between humans and animals ensure rapid dissemination of recent transfer events, as does modern transport and globalization. As resistance to antibacterial compounds becomes universal, there is likely to be increasing selection pressure for phenotypes with adverse consequences for human welfare, such as enhanced virulence, pathogenicity, and transmission. Improved understanding of the ecology of LGT could help us devise strategies to control this fundamental evolutionary process.

A biophysical model of supercoiling dependent transcription predicts a structural aspect to gene regulation.

BackgroundTranscription in Escherichia coli generates positive supercoiling in the DNA, which is relieved by the enzymatic activity of gyrase. Recently published experimental evidence suggests that transcription initiation and elongation are inhibited by the buildup of positive supercoiling. It has therefore been proposed that intermittent binding of gyrase plays a role in transcriptional bursting. Considering that transcription is one of the most fundamental cellular processes, it is desirable to be able to account for the buildup and release of positive supercoiling in models of transcription.ResultsHere we present a detailed biophysical model of gene expression that incorporates the effects of supercoiling due to transcription. By directly linking the amount of positive supercoiling to the rate of transcription, the model predicts that highly transcribed genes’ mRNA distributions should substantially deviate from Poisson distributions, with enhanced density at low mRNA copy numbers. Additionally, the model predicts a high degree of correlation between expression levels of genes inside the same supercoiling domain.ConclusionsOur model, incorporating the supercoiling state of the gene, makes specific predictions that differ from previous models of gene expression. Genes in the same supercoiling domain influence the expression level of neighboring genes. Such structurally dependent regulation predicts correlations between genes in the same supercoiling domain. The topology of the chromosome therefore creates a higher level of gene regulation, which has broad implications for understanding the evolution and organization of bacterial genomes.Electronic supplementary materialThe online version of this article (doi:10.1186/s13628-016-0027-0) contains supplementary material, which is available to authorized users.

Limited dissemination of the wastewater treatment plant core resistome.

Horizontal gene transfer is a major contributor to the evolution of bacterial genomes and can facilitate the dissemination of antibiotic resistance genes between environmental reservoirs and potential pathogens. Wastewater treatment plants (WWTPs) are believed to play a central role in the dissemination of antibiotic resistance genes. However, the contribution of the dominant members of the WWTP resistome to resistance in human pathogens remains poorly understood. Here we use a combination of metagenomic functional selections and comprehensive metagenomic sequencing to uncover the dominant genes of the WWTP resistome. We find that this core resistome is unique to the WWTP environment, with <10% of the resistance genes found outside the WWTP environment. Our data highlight that, despite an abundance of functional resistance genes within WWTPs, only few genes are found in other environments, suggesting that the overall dissemination of the WWTP resistome is comparable to that of the soil resistome.

Multiple conversion between the genes encoding bacterial class-I release factors.

Bacteria require two class-I release factors, RF1 and RF2, that recognize stop codons and promote peptide release from the ribosome. RF1 and RF2 were most likely established through gene duplication followed by altering their stop codon specificities in the common ancestor of extant bacteria. This scenario expects that the two RF gene families have taken independent evolutionary trajectories after the ancestral gene duplication event. However, we here report two independent cases of conversion between RF1 and RF2 genes (RF1-RF2 gene conversion), which were severely examined by procedures incorporating the maximum-likelihood phylogenetic method. In both cases, RF1-RF2 gene conversion was predicted to occur in the region encoding nearly entire domain 3, of which functions are common between RF paralogues. Nevertheless, the ‘direction’ of gene conversion appeared to be opposite from one another—from RF2 gene to RF1 gene in one case, while from RF1 gene to RF2 gene in the other. The two cases of RF1-RF2 gene conversion prompt us to propose two novel aspects in the evolution of bacterial class-I release factors: (i) domain 3 is interchangeable between RF paralogues, and (ii) RF1-RF2 gene conversion have occurred frequently in bacterial genome evolution.

CRISPR-based screening of genomic island excision events in bacteria

Genomic analysis of Streptococcus thermophilus revealed that mobile genetic elements (MGEs) likely contributed to gene acquisition and loss during evolutionary adaptation to milk. Clustered regularly interspaced short palindromic repeats-CRISPR-associated genes (CRISPR-Cas), the adaptive immune system in bacteria, limits genetic diversity by targeting MGEs including bacteriophages, transposons, and plasmids. CRISPR-Cas systems are widespread in streptococci, suggesting that the interplay between CRISPR-Cas systems and MGEs is one of the driving forces governing genome homeostasis in this genus. To investigate the genetic outcomes resulting from CRISPR-Cas targeting of integrated MGEs, in silico prediction revealed four genomic islands without essential genes in lengths from 8 to 102 kbp, totaling 7% of the genome. In this study, the endogenous CRISPR3 type II system was programmed to target the four islands independently through plasmid-based expression of engineered CRISPR arrays. Targeting lacZ within the largest 102-kbp genomic island was lethal to wild-type cells and resulted in a reduction of up to 2.5-log in the surviving population. Genotyping of Lac(-) survivors revealed variable deletion events between the flanking insertion-sequence elements, all resulting in elimination of the Lac-encoding island. Chimeric insertion sequence footprints were observed at the deletion junctions after targeting all of the four genomic islands, suggesting a common mechanism of deletion via recombination between flanking insertion sequences. These results established that self-targeting CRISPR-Cas systems may direct significant evolution of bacterial genomes on a population level, influencing genome homeostasis and remodeling.

Signatures of host/symbiont genome coevolution in insect nutritional endosymbioses.

The role of symbiosis in bacterial symbiont genome evolution is well understood, yet the ways that symbiosis shapes host genomes or more particularly, host/symbiont genome coevolution in the holobiont is only now being revealed. Here, we identify three coevolutionary signatures that characterize holobiont genomes. The first signature, host/symbiont collaboration, arises when completion of essential pathways requires host/endosymbiont genome complementarity. Metabolic collaboration has evolved numerous times in the pathways of amino acid and vitamin biosynthesis. Here, we highlight collaboration in branched-chain amino acid and pantothenate (vitamin B5) biosynthesis. The second coevolutionary signature is acquisition, referring to the observation that holobiont genomes acquire novel genetic material through various means, including gene duplication, lateral gene transfer from bacteria that are not their current obligate symbionts, and full or partial endosymbiont replacement. The third signature, constraint, introduces the idea that holobiont genome evolution is constrained by the processes governing symbiont genome evolution. In addition, we propose that collaboration is constrained by the expression profile of the cell lineage from which endosymbiont-containing host cells, called bacteriocytes, are derived. In particular, we propose that such differences in bacteriocyte cell lineage may explain differences in patterns of host/endosymbiont metabolic collaboration between the sap-feeding suborders Sternorrhyncha and Auchenorrhynca. Finally, we review recent studies at the frontier of symbiosis research that are applying functional genomic approaches to characterization of the developmental and cellular mechanisms of host/endosymbiont integration, work that heralds a new era in symbiosis research.

Characterization of Insertion Sequence ISSau2 in the Human and Livestock-Associated Staphylococcus aureus.

Mobile genetic elements play important roles in evolution and diversification of bacterial genomes. ISSau2 is 1660bp in length with terminal 5’-TG and CA-3’ dinucleotides and has two overlapping reading frames orfA and orfB. It has been found in a wide range of S. aureus, such as HA-MRSA252, LGA251, MRSA S0385 and ED133. To determine distribution of ISSau2, 164 S. aureus isolates from milk samples of mastitic cows from our laboratory and all the S. aureus strains from the National Center for Biotechnology Information (NCBI) database were screened for the presence of ISSau2. Next, in order to explore a potential relationship among S. aureus ISSau2-containing strains and isolates, a relationship among 10 ISSau2-positive S. aureus isolates and 27 ISSau2-positive S. aureus strains was investigated by a phylogenetic analysis. These ISSau2 isolates and strains could be classified into four groups (A, B, C and D). The strains or isolates in Group D were all isolated from mammary glands, suggesting tissue specificity. All strains in Group B had an identical ISSau2 derivative, termed ISSau2 1628, with 32bp deletion at the 3’ terminus. ISSau2 1628 in strain ST398 from Group B was closely related to ISSau2 in strain LGA251 from Group D.

Bacterial genome evolution in superspecies systems: An approach to the reconstruction of symbiogenesis processes

Bacteria form a broad spectrum of symbioses with eukaryotes. This permits reconstruction of the symbiogenesis processes providing the transformation of free-living microorganisms into cellular organelles. In ecologically (conditionally) obligate symbioses, an increase in the size and complexity of the bacterial genome structure was observed. This was associated with segregation of the regions controlling symbiosis into gene clusters, islands, and plasmids. In genetically (strictly) obligate symbioses, a reduction of "nonsymbiotic" regions of microbial genome occurs, which could begin from genes encoding metabolic and regulatory functions. It is extended towards genes encoding template processes. Conditionally obligate symbioses are characterised by the activation of horizontal gene transfer between various forms of microsymbionts, while for strictly obligate intracellular symbioses an activation of endo-symbiotic gene transfer between microsymbionts and their hosts was detected. The latter is responsible for bacterial transition from the functional (based on gene cross-regulation) to structural (based on recombination) genetic integration with hosts, which later could be followed by the complete assimilation of microbial genomes. In α-proteobacteria this evolutionary pathway could result in the formation of cellular organelles that are deficient in their own genomes but capable of preserving proteomic and cytological traits as a result of the gene-product import synthesized in cytosol (hydrogenosomes and mitosomes). The symbiogenic evolution of cyanobacteria could result in the loss of the plasmids generated from them, while the host maintains a significant part of their genome in nuclear chromosomes.

The CRISPR-Cas Immune System and Genetic Transfers: Reaching an Equilibrium.

Horizontal gene transfer drives the evolution of bacterial genomes, including the adaptation to changing environmental conditions. Exogenous DNA can enter a bacterial cell through transformation (free DNA or plasmids) or through the transfer of mobile genetic elements by conjugation (plasmids) and transduction (bacteriophages). Favorable genes can be acquired, but undesirable traits can also be inadvertently acquired through these processes. Bacteria have systems, such as clustered regularly interspaced short palindromic repeat CRISPR-associated genes (CRISPR-Cas), that can cleave foreign nucleic acid molecules. In this review, we discuss recent advances in understanding CRISPR-Cas system activity against mobile genetic element transfer through transformation and conjugation. We also highlight how CRISPR-Cas systems influence bacterial evolution and how CRISPR-Cas components affect plasmid replication.

GC-Content Evolution in Bacterial Genomes: The Biased Gene Conversion Hypothesis Expands

The characterization of functional elements in genomes relies on the identification of the footprints of natural selection. In this quest, taking into account neutral evolutionary processes such as mutation and genetic drift is crucial because these forces can generate patterns that may obscure or mimic signatures of selection. In mammals, and probably in many eukaryotes, another such confounding factor called GC-Biased Gene Conversion (gBGC) has been documented. This mechanism generates patterns identical to what is expected under selection for higher GC-content, specifically in highly recombining genomic regions. Recent results have suggested that a mysterious selective force favouring higher GC-content exists in Bacteria but the possibility that it could be gBGC has been excluded. Here, we show that gBGC is probably at work in most if not all bacterial species. First we find a consistent positive relationship between the GC-content of a gene and evidence of intra-genic recombination throughout a broad spectrum of bacterial clades. Second, we show that the evolutionary force responsible for this pattern is acting independently from selection on codon usage, and could potentially interfere with selection in favor of optimal AU-ending codons. A comparison with data from human populations shows that the intensity of gBGC in Bacteria is comparable to what has been reported in mammals. We propose that gBGC is not restricted to sexual Eukaryotes but also widespread among Bacteria and could therefore be an ancestral feature of cellular organisms. We argue that if gBGC occurs in bacteria, it can account for previously unexplained observations, such as the apparent non-equilibrium of base substitution patterns and the heterogeneity of gene composition within bacterial genomes. Because gBGC produces patterns similar to positive selection, it is essential to take this process into account when studying the evolutionary forces at work in bacterial genomes.

Plasmids are vectors for redundant chromosomal genes in the Bacillus cereus group.

BackgroundProkaryotic plasmids have played significant roles in the evolution of bacterial genomes and have a great impact on the metabolic functions of the host cell. Many bacterial strains contain multiple plasmids, but the relationships between bacterial plasmids and chromosomes are unclear. We focused on plasmids from the Bacillus cereus group because most strains contain several plasmids.ResultsWe collected the genome sequences of 104 plasmids and 20 chromosomes from B. cereus group strains, and we studied the relationships between plasmids and chromosomes by focusing on the pan-genomes of these plasmids and chromosomes. In terms of basic features (base composition and codon usage), the genes on plasmids were more similar to the chromosomal variable genes (distributed genes and unique genes) than to the chromosomal core genes. Although all the functional categories of the chromosomal genes were exhibited by the plasmid genes, the proportions of each category differed between these two gene sets. The 598 gene families shared between chromosomes and plasmids displayed a uniform distribution between the two groups. A phylogenetic analysis of the shared genes, including the chromosomal core gene set, indicated that gene exchange events between plasmids and chromosomes occurred frequently during the evolutionary histories of the strains and species in this group. Moreover, the shared genes between plasmids and chromosomes usually had different promoter and terminator sequences, suggesting that they are regulated by different elements at the transcriptional level.ConclusionsWe speculate that for the entire B. cereus group, adaptive genes are preserved on both plasmids and chromosomes; however, in a single cell, homologous genes on plasmids and the chromosome are controlled by different regulators to reduce the burden of maintaining redundant genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-014-1206-5) contains supplementary material, which is available to authorized users.

The conjugative DNA-transfer apparatus of Streptomyces

Conjugation is a major route of horizontal gene transfer, an important driving force in the evolution of bacterial genomes. Since antibiotic producing streptomycetes represent a natural reservoir of antibiotic resistance genes, the Streptomyces conjugation system might have a particular role in the dissemination of the resistance genes.Streptomycetes transfer DNA in a unique process, clearly distinguished from the well-known DNA-transfer by type IV secretion systems. A single plasmid-encoded DNA-translocase, TraB, transfers a double-stranded DNA-molecule to the recipient. Elucidation of the structure, pore forming ability and DNA binding characteristics of TraB indicated that the TraB conjugation system is derived from an FtsK-like ancestor protein suggesting that Streptomyces adapted the FtsK/SpoIIIE chromosome segregation system to transfer DNA between two distinct Streptomyces cells. Following the primary transfer, a multi-protein DNA-translocation apparatus consisting of TraB and several Spd-proteins spreads the newly transferred DNA to the neighbouring mycelial compartments resulting in the rapid colonization of the recipient mycelium by the donor DNA.

The stability region of the Streptomyces lividans plasmid pIJ101 encodes a DNA-binding protein recognizing a highly conserved short palindromic sequence motif.

Conjugation is a driving force in the evolution and shaping of bacterial genomes. In antibiotic producing streptomycetes even small plasmids replicating via the rolling-circle mechanism are conjugative. Although they encode only genes involved in replication and transfer, the molecular function of most plasmid encoded proteins is unknown. In this work we show that the conjugative plasmid pIJ101 encodes an overlooked protein, SpdA2. We show that SpdA2 is a DNA binding protein which specifically recognizes a palindromic DNA sequence (sps). sps is localized within the spdA2 coding region and highly conserved in many Streptomyces plasmids. Elimination of the palindrome or deletion of spdA2 in plasmid pIJ303 did not interfere with conjugative plasmid transfer or pock formation, but affected segregational stability.

Rapid editing and evolution of bacterial genomes using libraries of synthetic DNA.

Multiplex automated genome engineering (MAGE) is a powerful technology for in vivo genome editing that uses synthetic single-stranded DNA (ssDNA) to introduce targeted modifications directly into the Escherichia coli chromosome. MAGE is a cyclical process that involves transformation of ssDNA (by electroporation) followed by outgrowth, during which bacteriophage homologous recombination proteins mediate annealing of ssDNAs to their genomic targets. By iteratively introducing libraries of mutagenic ssDNAs targeting multiple sites, MAGE can generate combinatorial genetic diversity in a cell population. Alternatively, MAGE can introduce precise mutant alleles at many loci for genome-wide editing or for recoding projects that are not possible with other methods. In recent technological advances, MAGE has been improved by strain modifications and selection techniques that enhance allelic replacement. This protocol describes the manual execution of MAGE wherein each cycle takes ≈ 2.5 h, which, if carried out by two people, allows ≈ 10 continuous cycles of MAGE-based mutagenesis per day.

Symbiosis plasmids bring their own mutagen to the wedding party.

Symbiosis plasmids bring their own mutagen to the wedding party.

Dynamic recruitment of amino acid transporters to the insect/symbiont interface

Symbiosis is well known to influence bacterial symbiont genome evolution and has recently been shown to shape eukaryotic host genomes. Intriguing patterns of host genome evolution, including remarkable numbers of gene duplications, have been observed in the pea aphid, a sap-feeding insect that relies on a bacterial endosymbiont for amino acid provisioning. Previously, we proposed that gene duplication has been important for the evolution of symbiosis based on aphid-specific gene duplication in amino acid transporters (AATs), with some paralogs highly expressed in the cells housing symbionts (bacteriocytes). Here, we use a comparative approach to test the role of gene duplication in enabling recruitment of AATs to bacteriocytes. Using genomic and transcriptomic data, we annotate AATs from sap-feeding and non sap-feeding insects and find that, like aphids, AAT gene families have undergone independent large-scale gene duplications in three of four additional sap-feeding insects. RNA-seq differential expression data indicate that, like aphids, the sap-feeding citrus mealybug possesses several lineage-specific bacteriocyte-enriched paralogs. Further, differential expression data combined with quantitative PCR support independent evolution of bacteriocyte enrichment in sap-feeding insect AATs. Although these data indicate that gene duplication is not necessary to initiate host/symbiont amino acid exchange, they support a role for gene duplication in enabling AATs to mediate novel host/symbiont interactions broadly in the sap-feeding suborder Sternorrhyncha. In combination with recent studies on other symbiotic systems, gene duplication is emerging as a general pattern in host genome evolution.

PanGP: A tool for quickly analyzing bacterial pan-genome profile

Summary: Pan-genome analyses have shed light on the dynamics and evolution of bacterial genome from the point of population. The explosive growth of bacterial genome sequence also brought an extremely big challenge to pan-genome profile analysis. We developed a tool, named PanGP, to complete pan-genome profile analysis for large-scale strains efficiently. PanGP has integrated two sampling algorithms, totally random (TR) and distance guide (DG). The DG algorithm drew sample strain combinations on the basis of genome diversity of bacterial population. The performance of these two algorithms have been evaluated on four bacteria populations with strain numbers varying from 30 to 200, and the DG algorithm exhibited overwhelming advantage on accuracy and stability than the TR algorithm.Availability: PanGP was developed with a user-friendly graphic interface and it was available at http://PanGP.big.ac.cn.Contact: xiaojingfa@big.ac.cn or wujy@big.ac.cnSupplementary information: Supplementary data are available at Bioinformatics online.

Predicting genetic and ecological characteristics of bacterial species by comparing the trajectories of dN/dS and dI/dS in bacterial genomes

Indel (insertion/deletion) causes gene disruption and is considered to be deleterious like non-synonymous mutation during the evolution of bacterial genomes. The trajectory of dN/dS (the ratio of non-synonymous to synonymous mutation) has been found to decrease exponentially over time, but the trajectory of dI/dS (the ratio of indel to synonymous mutation) has not been thoroughly explored yet. Here we compared the patterns of dN/dS and dI/dS for several bacterial species. The majority of them showed a much steeper dI/dS trajectory than the dN/dS trajectory, suggesting that indel was more deleterious than non-synonymous mutation and therefore was more difficult to fix in genomes. However, the naturally competent bacteria, or those with a much lower genetic barrier for DNA exchange, presented a reverse relationship between dI/dS and dN/dS, indicative of an exceptional ability to tolerate horizontal genetic transfer. The result suggests that plotting of dN/dS and dI/dS trajectories can help to predict the genetic and ecological characteristics of bacterial species.

An algebraic view of bacterial genome evolution

Rearrangements of bacterial chromosomes can be studied mathematically at several levels, most prominently at a local, or sequence level, as well as at a topological level. The biological changes involved locally are inversions, deletions, and transpositions, while topologically they are knotting and catenation. These two modelling approaches share some surprising algebraic features related to braid groups and Coxeter groups. The structural approach that is at the core of algebra has long found applications in sciences such as physics and analytical chemistry, but only in a small number of ways so far in biology. And yet there are examples where an algebraic viewpoint may capture a deeper structure behind biological phenomena. This article discusses a family of biological problems in bacterial genome evolution for which this may be the case, and raises the prospect that the tools developed by algebraists over the last century might provide insight to this area of evolutionary biology.