Evolutionary dynamics of full genome content in Escherichia coli

Characterizing the loss and acquisition of Tn2-like transposon in Klebsiella pneumoniae: Implications for carbapenemase gene dissemination

Over 1 year, two KPC-producing and two non-KPC-producing Klebsiella pneumoniae strains were isolated from a patient. Genome and DNA hybridization analyses revealed the first three strains as a clonal lineage, with carbapenem resistance changes due to a Tn2-like transposon on an IncR/IncFII plasmid. The fourth strain, carrying three plasmids, caused a lethal infection and represented a different lineage. All strains belonged to the ST11-SL47-OL101 type. This study highlights the Tn2-like transposon's role in carbapenemase gene spread and the importance of distinguishing between bacterial colonization and infection.

OUP accepted manuscript

Escherichia coli has a rich history as biology's ‘rock star’, driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance—including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Molecular insights into the genome dynamics and interactions between core and acquired genomes of Vibrio cholerae

Bacterial species are hosts to horizontally acquired mobile genetic elements (MGEs), which encode virulence, toxin, antimicrobial resistance, and other metabolic functions. The bipartite genome of Vibrio cholerae harbors sporadic and conserved MGEs that contribute in the disease development and survival of the pathogens. For a comprehensive understanding of dynamics of MGEs in the bacterial genome, we engineered the genome of V. cholerae and examined in vitro and in vivo stability of genomic islands (GIs), integrative conjugative elements (ICEs), and prophages. Recombinant vectors carrying the integration module of these GIs, ICE and CTXΦ, helped us to understand the efficiency of integrations of MGEs in the V. cholerae chromosome. We have deleted more than 250 acquired genes from 6 different loci in the V. cholerae chromosome and showed contribution of CTX prophage in the essentiality of SOS response master regulator LexA, which is otherwise not essential for viability in other bacteria, including Escherichia coli In addition, we observed that the core genome-encoded RecA helps CTXΦ to bypass V. cholerae immunity and allow it to replicate in the host bacterium in the presence of similar prophage in the chromosome. Finally, our proteomics analysis reveals the importance of MGEs in modulating the levels of cellular proteome. This study engineered the genome of V. cholerae to remove all of the GIs, ICEs, and prophages and revealed important interactions between core and acquired genomes.

Comparative Analysis of Tools and Approaches for Source Tracking Listeria monocytogenes in a Food Facility Using Whole-Genome Sequence Data

As WGS is increasingly used by food industry to characterize pathogen isolates, users are challenged by the variety of analysis approaches available, ranging from methods that require extensive bioinformatics expertise to commercial software packages. This study aimed to assess the impact of analysis pipelines (i.e., different hqSNP pipelines, a cg/wgMLST pipeline) and the reference genome selection on analysis results (i.e., hqSNP and allelic differences as well as tree topologies) and conclusion drawn. For these comparisons, whole genome sequences were obtained for 40 Listeria monocytogenes isolates collected over 18 years from a cold-smoked salmon facility and 2 other isolates obtained from different facilities as part of academic research activities; WGS data were analyzed with three hqSNP pipelines and two MLST pipelines. After initial clustering using a k-mer based approach, hqSNP pipelines were run using two types of reference genomes: (i) closely related closed genomes (“closed references”) and (ii) high-quality de novo assemblies of the dataset isolates (“draft references”). All hqSNP pipelines identified similar hqSNP difference ranges among isolates in a given cluster; use of different reference genomes showed minimal impacts on hqSNP differences identified between isolate pairs. Allelic differences obtained by wgMLST showed similar ranges as hqSNP differences among isolates in a given cluster; cgMLST consistently showed fewer differences than wgMLST. However, phylogenetic trees and dendrograms, obtained based on hqSNP and cg/wgMLST data, did show some incongruences, typically linked to clades supported by low bootstrap values in the trees. When a hqSNP cutoff was used to classify isolates as “related” or “unrelated,” use of different pipelines yielded a considerable number of discordances; this finding supports that cut-off values are valuable to provide a starting point for an investigation, but supporting and epidemiological evidence should be used to interpret WGS data. Overall, our data suggest that cgMLST-based data analyses provide for appropriate subtype differentiation and can be used without the need for preliminary data analyses (e.g., k-mer based clustering) or external closed reference genomes, simplifying data analyses needs. hqSNP or wgMLST analyses can be performed on the isolate clusters identified by cgMLST to increase the precision on determining the genomic similarity between isolates.

Metabolic Determinants of Enzyme Evolution in a Genome-Scale Bacterial Metabolic Network

Different genes and proteins evolve at very different rates. To identify the factors that explain these differences is an important aspect of research in molecular evolution. One such factor is the role a protein plays in a large molecular network. Here, we analyze the evolutionary rates of enzyme-coding genes in the genome-scale metabolic network of Escherichia coli to find the evolutionary constraints imposed by the structure and function of this complex metabolic system. Central and highly connected enzymes appear to evolve more slowly than less connected enzymes, but we find that they do so as a by-product of their high abundance, and not because of their position in the metabolic network. In contrast, enzymes catalyzing reactions with high metabolic flux-high substrate to product conversion rates-evolve slowly even after we account for their abundance. Moreover, enzymes catalyzing reactions that are difficult to by-pass through alternative pathways, such that they are essential in many different genetic backgrounds, also evolve more slowly. Our analyses show that an enzyme's role in the function of a metabolic network affects its evolution more than its place in the network's structure. They highlight the value of a system-level perspective for studies of molecular evolution.

Escherichia coli: an old friend with new tidings

Escherichia coli is one of the most-studied microorganisms worldwide but its characteristics are continually changing. Extraintestinal E. coli infections, such as urinary tract infections and neonatal sepsis, represent a huge public health problem. They are caused mainly by specialized extraintestinal pathogenic E. coli (ExPEC) strains that can innocuously colonize human hosts but can also cause disease upon entering a normally sterile body site. The virulence capability of such strains is determined by a combination of distinctive accessory traits, called virulence factors, in conjunction with their distinctive phylogenetic background. It is conceivable that by developing interventions against the most successful ExPEC lineages or their key virulence/colonization factors the associated burden of disease and health care costs could foreseeably be reduced in the future. On the other hand, one important problem worldwide is the increase of antimicrobial resistance shown by bacteria. As underscored in the last WHO global report, within a wide range of infectious agents including E. coli, antimicrobial resistance has reached an extremely worrisome situation that ‘threatens the achievements of modern medicine’. In the present review, an update of the knowledge about the pathogenicity, antimicrobial resistance and clinical aspects of this ‘old friend’ was presented.

Modelling Prokaryote Gene Content

The patchy distribution of genes across the prokaryotes may be caused by multiple gene losses or lateral transfer. Probabilistic models of gene gain and loss are needed to distinguish between these possibilities. Existing models allow only single genes to be gained and lost, despite the empirical evidence for multi-gene events. We compare birth-death models (currently the only widely-used models, in which only one gene can be gained or lost at a time) to blocks models (allowing gain and loss of multiple genes within a family). We analyze two pairs of genomes: two E. coli strains, and the distantly-related Archaeoglobus fulgidus (archaea) and Bacillus subtilis (gram positive bacteria). Blocks models describe the data much better than birth-death models. Our models suggest that lateral transfers of multiple genes from the same family are rare (although transfers of single genes are probably common). For both pairs, the estimated median time that a gene will remain in the genome is not much greater than the time separating the common ancestors of the archaea and bacteria. Deep phylogenetic reconstruction from sequence data will therefore depend on choosing genes likely to remain in the genome for a long time. Phylogenies based on the blocks model are more biologically plausible than phylogenies based on the birth-death model.

Evolutionary and Functional Relationships of the dha Regulon by Genomic Context Analysis

3-hydroxypropionaldehyde (3-HPA) and 1,3-propanediol (1,3-PD) are subproducts of glycerol degradation and of economical interest as they are used for polymers synthesis, such as polyesters and polyurethanes. Some few characterized bacterial species (mostly from Firmicutes and Gamma-proteobacteria groups) are able to catabolize these monomers from glycerol using the gene products from the dha regulon. To expand our knowledge and direct further experimental studies on the regulon and related genes for the anaerobic glycerol metabolism, an extensive genomic screening was performed to identify the presence of the dha genes in fully sequenced prokaryotic genomes. Interestingly, this work shows that although only few bacteria species are known to produce 3-HPA or 1,3-PD, the incomplete regulon is found in more than 100 prokaryotic genomes. However, the complete pathway is found only in a few dozen species belonging to five different taxonomic groups, including one Archaea species, Halalkalicoccus jeotgali. Phylogenetic analysis and conservation of both gene synteny and primary sequence similarity reinforce the idea that these genes have a common origin and were possibly acquired by lateral gene transfer (LGT). Besides the evolutionary aspect, the identification of homologs from several different organisms may predict potential alternative targets for faster or more efficient biological synthesis of 3-HPA or 1,3-PD.

Inter-individual differences in the gene content of human gut bacterial species

BACKGROUND

Gene content differences in human gut microbes can lead to inter-individual phenotypic variations such as digestive capacity. It is unclear whether gene content variation is caused by differences in microbial species composition or by the presence of different strains of the same species; the extent of gene content variation in the latter is unknown. Unlike pan-genome studies of cultivable strains, the use of metagenomic data can provide an unbiased view of structural variation of gut bacterial strains by measuring them in their natural habitats, the gut of each individual in this case, representing native boundaries between gut bacterial populations. We analyzed publicly available metagenomic data from fecal samples to characterize inter-individual variation in gut bacterial species.

RESULTS

A comparison of 11 abundant gut bacterial species showed that the gene content of strains from the same species differed, on average, by 13% between individuals. This number is based on gene deletions only and represents a lower limit, yet the variation is already in a similar range as observed between completely sequenced strains of cultivable species. We show that accessory genes that differ considerably between individuals can encode important functions, such as polysaccharide utilization and capsular polysaccharide synthesis loci.

CONCLUSION

Metagenomics can yield insights into gene content variation of strains in complex communities, which cannot be predicted by phylogenetic marker genes alone. The large degree of inter-individual variability in gene content implies that strain resolution must be considered in order to fully assess the functional potential of an individual's human gut microbiome.

Two rep genes in small cryptic plasmid pKST21 of Escherichia coli

The complete nucleotide sequence of a small cryptic plasmid pKST21 from Escherichia coli was determined. This plasmid is 1,460 bp long with an overall GC content of 51 %. Based on sequence analysis, the presence of two segments with different average GC density was observed. The segment with higher GC content revealed 98-90 % similarity to several small plasmids of E. coli and to pCR1 from Gram-positive Corynebacterium renale. Plasmid pKST21 possesses two conversely oriented open reading frames encoding proteins with a high degree of amino acid identity to Rep proteins involved in replication. ORF1 encodes replication protein similar to RepA protein of Bartonella tribocorum or Bacillus cereus plasmids or to the putative plasmid Rep protein from ecologically close Selenomonas ruminantium. ORF2 similarly encodes a replication protein, which shares 97 % homology with Rep protein from C. renale. Genetic diversity observed in plasmid pKST21 indicates a mosaic structure of the plasmid with different segments acquired from different sources. Deletion analysis showed that both fragments carrying the repA and repB genes are necessary for the replication of pKST21 in E. coli. The presence of plasmid with the same gene composition was revealed in 14 % of tested E. coli isolates from the rumen of sheep. All these strains produced identical ERIC-PCR profiles indicating isogenic origin of the strain and lack of horizontal gene transfer of pKST21 plasmid.

Vaginal microbiome: rethinking health and disease

Vaginal microbiota form a mutually beneficial relationship with their host and have a major impact on health and disease. In recent years our understanding of vaginal bacterial community composition and structure has significantly broadened as a result of investigators using cultivation-independent methods based on the analysis of 16S ribosomal RNA (rRNA) gene sequences. In asymptomatic, otherwise healthy women, several kinds of vaginal microbiota exist, the majority often dominated by species of Lactobacillus, while others are composed of a diverse array of anaerobic microorganisms. Bacterial vaginosis is the most common vaginal condition and is vaguely characterized as the disruption of the equilibrium of the normal vaginal microbiota. A better understanding of normal and healthy vaginal ecosystems that is based on their true function and not simply on their composition would help better define health and further improve disease diagnostics as well as the development of more personalized regimens to promote health and treat diseases.

Evolution of variation in presence and absence of genes in bacterial pathways

Background

Bacterial genomes exhibit a remarkable degree of variation in the presence and absence of genes, which probably extends to the level of individual pathways. This variation may be a consequence of the significant evolutionary role played by horizontal gene transfer, but might also be explained by the loss of genes through mutation. A challenge is to understand why there would be variation in gene presence within pathways if they confer a benefit only when complete.

Results

Here, we develop a mathematical model to study how variation in pathway content is produced by horizontal transfer, gene loss and partial exposure of a population to a novel environment.

Conclusions

We discuss the possibility that variation in gene presence acts as cryptic genetic variation on which selection acts when the appropriate environment occurs. We find that a high level of variation in gene presence can be readily explained by decay of the pathway through mutation when there is no longer exposure to the selective environment, or when selection becomes too weak to maintain the genes. In the context of pathway variation the role of horizontal gene transfer is probably the initial introduction of a complete novel pathway rather than in building up the variation in a genome without the pathway.

After genomics, what proteomics tools could help us understand the antimicrobial resistance of Escherichia coli?

Proteomic approaches have been considerably improved during the past decade and have been used to investigate the differences in protein expression profiles of cells grown under a broad spectrum of growth conditions and with different stress factors including antibiotics. In Europe, the most significant disease threat remains the presence of microorganisms that have become resistant to antimicrobials and so it is important that different scientific tools are combined to achieve the largest amount of knowledge in this area of expertise. The emergence and spread of the antibiotic-resistant Gram-negative pathogens, such as Escherichia coli, can lead to serious problem public health in humans. E. coli, a very well described prokaryote, has served as a model organism for several biological and biotechnological studies increasingly so since the completion of the E. coli genome-sequencing project. The purpose of this review is to present an overview of the different proteomic approaches to antimicrobial-resistant E. coli that will be helpful to obtain a better knowledge of the antibiotic-resistant mechanism(s). This can also aid to understand the molecular determinants involved with pathogenesis, which is essential for the development of effective strategies to combat infection and to reveal new therapeutic targets. This article is part of a Special Issue entitled: Proteomics: The clinical link.

A quantitative account of genomic island acquisitions in prokaryotes

Background

Microbial genomes do not merely evolve through the slow accumulation of mutations, but also, and often more dramatically, by taking up new DNA in a process called horizontal gene transfer. These innovation leaps in the acquisition of new traits can take place via the introgression of single genes, but also through the acquisition of large gene clusters, which are termed Genomic Islands. Since only a small proportion of all the DNA diversity has been sequenced, it can be hard to find the appropriate donors for acquired genes via sequence alignments from databases. In contrast, relative oligonucleotide frequencies represent a remarkably stable genomic signature in prokaryotes, which facilitates compositional comparisons as an alignment-free alternative for phylogenetic relatedness.

In this project, we test whether Genomic Islands identified in individual bacterial genomes have a similar genomic signature, in terms of relative dinucleotide frequencies, and can therefore be expected to originate from a common donor species.

Results

When multiple Genomic Islands are present within a single genome, we find that up to 28% of these are compositionally very similar to each other, indicative of frequent recurring acquisitions from the same donor to the same acceptor.

Conclusions

This represents the first quantitative assessment of common directional transfer events in prokaryotic evolutionary history. We suggest that many of the resident Genomic Islands per prokaryotic genome originated from the same source, which may have implications with respect to their regulatory interactions, and for the elucidation of the common origins of these acquired gene clusters.

A close relationship between primary nucleotides sequence structure and the composition of functional genes in the genome of prokaryotes

Comparative genomics is an essential tool to unravel how genomes change over evolutionary time and to gain clues on the links between functional genomics and evolution. In prokaryotes, the large, good quality, genome sequences available in public databases and the recently developed large-scale computational methods, offer an unprecedent view on the ecology and evolution of microorganisms through comparative genomics. In this work, we examined the links among genome structure (i.e., the sequential distribution of nucleotides itself by detrended fluctuation analysis, DFA) and genomic diversity (i.e., gene functionality by Clusters of Orthologous Genes, COGs) in 828 full sequenced prokaryotic genomes from 548 different bacteria and archaea species. DFA scaling exponent α indicated persistent long-range correlations (fractality) in each genome analyzed. Higher resolution power was found when considering the sequential succession of purine (AG) vs. pyrimidine (CT) bases than either keto (GT) to amino (AC) forms or strongly (GC) vs. weakly (AT) bonded nucleotides. Interestingly, the phyla Aquificae, Fusobacteria, Dictyoglomi, Nitrospirae, and Thermotogae were closer to archaea than to their bacterial counterparts. A strong significant correlation was found between scaling exponent α and COGs distribution, and we consistently observed that the larger α the more heterogeneous was the gene distribution within each functional category, suggesting a close relationship between primary nucleotides sequence structure and functional genes composition.

Survival of Escherichia coli in the environment: fundamental and public health aspects

In this review, our current understanding of the species Escherichia coli and its persistence in the open environment is examined. E. coli consists of six different subgroups, which are separable by genomic analyses. Strains within each subgroup occupy various ecological niches, and can be broadly characterized by either commensalistic or different pathogenic behaviour. In relevant cases, genomic islands can be pinpointed that underpin the behaviour. Thus, genomic islands of, on the one hand, broad environmental significance, and, on the other hand, virulence, are highlighted in the context of E. coli survival in its niches. A focus is further placed on experimental studies on the survival of the different types of E. coli in soil, manure and water. Overall, the data suggest that E. coli can persist, for varying periods of time, in such terrestrial and aquatic habitats. In particular, the considerable persistence of the pathogenic E. coli O157:H7 is of importance, as its acid tolerance may be expected to confer a fitness asset in the more acidic environments. In this context, the extent to which E. coli interacts with its human/animal host and the organism's survivability in natural environments are compared. In addition, the effect of the diversity and community structure of the indigenous microbiota on the fate of invading E. coli populations in the open environment is discussed. Such a relationship is of importance to our knowledge of both public and environmental health.

Inferring bacterial genome flux while considering truncated genes

Bacterial gene content variation during the course of evolution has been widely acknowledged and its pattern has been actively modeled in recent years. Gene truncation or gene pseudogenization also plays an important role in shaping bacterial genome content. Truncated genes could also arise from small-scale lateral gene transfer events. Unfortunately, the information of truncated genes has not been considered in any existing mathematical models on gene content variation. In this study, we developed a model to incorporate truncated genes. Maximum-likelihood estimates (MLEs) of the new model reveal fast rates of gene insertions/deletions on recent branches, suggesting a fast turnover of many recently transferred genes. The estimates also suggest that many truncated genes are in the process of being eliminated from the genome. Furthermore, we demonstrate that the ignorance of truncated genes in the estimation does not lead to a systematic bias but rather has a more complicated effect. Analysis using the new model not only provides more accurate estimates on gene gains/losses (or insertions/deletions), but also reduces any concern of a systematic bias from applying simplified models to bacterial genome evolution. Although not a primary purpose, the model incorporating truncated genes could be potentially used for phylogeny reconstruction using gene family content.

Conditions for the evolution of gene clusters in bacterial genomes

Genes encoding proteins in a common pathway are often found near each other along bacterial chromosomes. Several explanations have been proposed to account for the evolution of these structures. For instance, natural selection may directly favour gene clusters through a variety of mechanisms, such as increased efficiency of coregulation. An alternative and controversial hypothesis is the selfish operon model, which asserts that clustered arrangements of genes are more easily transferred to other species, thus improving the prospects for survival of the cluster. According to another hypothesis (the persistence model), genes that are in close proximity are less likely to be disrupted by deletions. Here we develop computational models to study the conditions under which gene clusters can evolve and persist. First, we examine the selfish operon model by re-implementing the simulation and running it under a wide range of conditions. Second, we introduce and study a Moran process in which there is natural selection for gene clustering and rearrangement occurs by genome inversion events. Finally, we develop and study a model that includes selection and inversion, which tracks the occurrence and fixation of rearrangements. Surprisingly, gene clusters fail to evolve under a wide range of conditions. Factors that promote the evolution of gene clusters include a low number of genes in the pathway, a high population size, and in the case of the selfish operon model, a high horizontal transfer rate. The computational analysis here has shown that the evolution of gene clusters can occur under both direct and indirect selection as long as certain conditions hold. Under these conditions the selfish operon model is still viable as an explanation for the evolution of gene clusters.

Genes involved in a common pathway or function are frequently found near each other on bacterial chromosomes. A number of hypotheses have been previously presented to explain this observation. A particularly influential theory is the selfish operon model, which posits that horizontal transfer could promote gene clustering by favouring transfer of arrangements of genes that are close together. Subsequent theoretical development and analysis of genomic data have contributed to the debate about the plausibility of this model. Here, by re-examining the evolutionary dynamics of gene clusters, we provide and discuss conditions under which gene clusters can evolve. We find that first, some form of bias for clustering is required for clusters to evolve. This bias can be in the form of bias in horizontal transfer towards genes that are close together, or direct natural selection for gene proximity. Our computational work does not present a theoretical obstacle to the selfish operon model as a possible explanation for the evolution of gene clusters.

Compensatory gene amplification restores fitness after inter-species gene replacements

Genes introduced by gene replacements and other types of horizontal gene transfer (HGT) represent a significant presence in many archaeal and eubacterial genomes. Most alien genes are likely to be neutral or deleterious upon arrival and their long-term persistence may require a mechanism that improves their selective contribution. To examine the fate of inter-species gene replacements, we exchanged three native S. typhimurium genes encoding ribosomal proteins with orthologues from various other microbes. The results show that replacement of each of these three genes reduces fitness to such an extent that it would provide an effective barrier against inter-species gene replacements in eubacterial populations. However, these fitness defects could be partially ameliorated by gene amplification that augmented the dosage of the heterologous proteins. This suggests that suboptimal expression is a common fitness constraint for inter-species gene replacements, with fitness costs conferred by either a lower expression level of the alien protein compared with the native protein or a requirement for an increased amount of the alien protein to maintain proper function. Our findings can explain the observation that duplicated genes are over-represented among horizontally transferred genes, and suggest a potential coupling between compensatory gene amplification after HGT and the evolution of new genes.

New approaches for functional genomic studies in staphylococci

Functional transcriptomics studies have resulted in interesting insights into Staphylococcus aureus diversity and pathogenicity. Here we review the principles, advantages and disadvantages of recent technical developments in the field of transcriptomics and their potential impact on S. aureus research.

Horizontal gene transfer of the secretome drives the evolution of bacterial cooperation and virulence

Background

Microbes engage in a remarkable array of cooperative behaviors, secreting shared proteins that are essential for foraging, shelter, microbial warfare, and virulence. These proteins are costly, rendering populations of cooperators vulnerable to exploitation by nonproducing cheaters arising by gene loss or migration. In such conditions, how can cooperation persist?

Results

Our model predicts that differential gene mobility drives intragenomic variation in investment in cooperative traits. More mobile loci generate stronger among-individual genetic correlations at these loci (higher relatedness) and thereby allow the maintenance of more cooperative traits via kin selection. By analyzing 21 Escherichia genomes, we confirm that genes coding for secreted proteins—the secretome—are very frequently lost and gained and are associated with mobile elements. We show that homologs of the secretome are overrepresented among human gut metagenomics samples, consistent with increased relatedness at secretome loci across multiple species. The biosynthetic cost of secreted proteins is shown to be under intense selective pressure, even more than for highly expressed proteins, consistent with a cost of cooperation driving social dilemmas. Finally, we demonstrate that mobile elements are in conflict with their chromosomal hosts over the chimeric ensemble's social strategy, with mobile elements enforcing cooperation on their otherwise selfish hosts via the cotransfer of secretome genes with “mafia strategy” addictive systems (toxin-antitoxin and restriction-modification).

Conclusion

Our analysis matches the predictions of our model suggesting that horizontal transfer promotes cooperation, as transmission increases local genetic relatedness at mobile loci and enforces cooperation on the resident genes. As a consequence, horizontal transfer promoted by agents such as plasmids, phages, or integrons drives microbial cooperation.

Evolutionary constraints permeate large metabolic networks

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Does gene translocation accelerate the evolution of laterally transferred genes?

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

Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the ∼18,000 families of orthologous genes, we found ∼2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genome's long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

Although abundant knowledge has been accumulated regarding the E. coli laboratory strain K-12, little is known about the evolutionary trajectories that have driven the high diversity observed among natural isolates of the species, which encompass both commensal and highly virulent intestinal and extraintestinal pathogenic strains. We have annotated or re-annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Although recombination rates are much higher than mutation rates, we were able to reconstruct a robust phylogeny based on the ∼2,000 genes common to all strains. Based on this phylogeny, we established the evolutionary scenario of gains and losses of thousands of specific genes, identifying functional classes under opposite selection pressures. This genome flux is confined to very few positions in the chromosome, which are the same for every genome. Notably, we identified few or no extraintestinal virulence-specific genes. We also defined a long-scale structure of recombination in the genome with lower recombination rates at the terminus of replication. These findings demonstrate that, despite a very high gene flow, genes can co-exist in an organised genome.

Uncovering rate variation of lateral gene transfer during bacterial genome evolution

BACKGROUND

Large scale genome arrangement, such as whole gene insertion/deletion, plays an important role in bacterial genome evolution. Various methods have been employed to study the dynamic process of gene insertions and deletions, such as parsimony methods and maximum likelihood methods. Previous maximum likelihood studies have assumed that the rate of gene insertions/deletions is constant over different genes. This assumption is unrealistic. For instance, it has been shown that informational genes are less likely to be laterally transferred than non-informational genes. However, how much of the variation in gene transfer rates is due to the difference between informational genes and non-informational genes is unclear. In this study, a Gamma-distribution was incorporated in the likelihood estimation by considering rate variation for gene insertions/deletions between genes. This makes it possible to address whether a difference between informational genes and non-informational genes is the main contributor to rate variation of lateral gene transfers.

RESULTS

The results show that models incorporating rate variation fit the data better than do constant rate models in many phylogenetic groups. Even though informational genes are less likely to be laterally transferred than non-informational genes, the degree of rate variation for insertions/deletions did not change dramatically and remained high even when informational genes were excluded from the study. This suggests that the variation in rate of insertions/deletions is not due mainly to the simple difference between informational genes and non-informational genes. Among genes that are not classified as informational and among the informational genes themselves, there are still large differences in the rates that these genes are inserted and deleted.

CONCLUSION

While the difference in informational gene rates contributes to rate variation, it is only a small fraction of the variation present; instead, a substantial amount of rate variation for insertions/deletions remains among both informational genes and among non-informational genes.

The emergence and fate of horizontally acquired genes in Escherichia coli

Bacterial species, and even strains within species, can vary greatly in their gene contents and metabolic capabilities. We examine the evolution of this diversity by assessing the distribution and ancestry of each gene in 13 sequenced isolates of Escherichia coli and Shigella. We focus on the emergence and demise of two specific classes of genes, ORFans (genes with no homologs in present databases) and HOPs (genes with distant homologs), since these genes, in contrast to most conserved ancestral sequences, are known to be a major source of the novel features in each strain. We find that the rates of gain and loss of these genes vary greatly among strains as well as through time, and that ORFans and HOPs show very different behavior with respect to their emergence and demise. Although HOPs, which mostly represent gene acquisitions from other bacteria, originate more frequently, ORFans are much more likely to persist. This difference suggests that many adaptive traits are conferred by completely novel genes that do not originate in other bacterial genomes. With respect to the demise of these acquired genes, we find that strains of Shigella lose genes, both by disruption events and by complete removal, at accelerated rates.

A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes

Background

Microarray based comparative genomic hybridisation (CGH) experiments have been used to study numerous biological problems including understanding genome plasticity in pathogenic bacteria. Typically such experiments produce large data sets that are difficult for biologists to handle. Although there are some programmes available for interpretation of bacterial transcriptomics data and CGH microarray data for looking at genetic stability in oncogenes, there are none specifically to understand the mosaic nature of bacterial genomes. Consequently a bottle neck still persists in accurate processing and mathematical analysis of these data. To address this shortfall we have produced a simple and robust CGH microarray data analysis process that may be automated in the future to understand bacterial genomic diversity.

Results

The process involves five steps: cleaning, normalisation, estimating gene presence and absence or divergence, validation, and analysis of data from test against three reference strains simultaneously. Each stage of the process is described and we have compared a number of methods available for characterising bacterial genomic diversity, for calculating the cut-off between gene presence and absence or divergence, and shown that a simple dynamic approach using a kernel density estimator performed better than both established, as well as a more sophisticated mixture modelling technique. We have also shown that current methods commonly used for CGH microarray analysis in tumour and cancer cell lines are not appropriate for analysing our data.

Conclusion

After carrying out the analysis and validation for three sequenced Escherichia coli strains, CGH microarray data from 19 E. coli O157 pathogenic test strains were used to demonstrate the benefits of applying this simple and robust process to CGH microarray studies using bacterial genomes.

Persistence drives gene clustering in bacterial genomes

BACKGROUND

Gene clustering plays an important role in the organization of the bacterial chromosome and several mechanisms have been proposed to explain its extent. However, the controversies raised about the validity of each of these mechanisms remind us that the cause of this gene organization remains an open question. Models proposed to explain clustering did not take into account the function of the gene products nor the likely presence or absence of a given gene in a genome. However, genomes harbor two very different categories of genes: those genes present in a majority of organisms - persistent genes - and those present in very few organisms - rare genes.

RESULTS

We show that two classes of genes are significantly clustered in bacterial genomes: the highly persistent and the rare genes. The clustering of rare genes is readily explained by the selfish operon theory. Yet, genes persistently present in bacterial genomes are also clustered and we try to understand why. We propose a model accounting specifically for such clustering, and show that indispensability in a genome with frequent gene deletion and insertion leads to the transient clustering of these genes. The model describes how clusters are created via the gene flux that continuously introduces new genes while deleting others. We then test if known selective processes, such as co-transcription, physical interaction or functional neighborhood, account for the stabilization of these clusters.

CONCLUSION

We show that the strong selective pressure acting on the function of persistent genes, in a permanent state of flux of genes in bacterial genomes, maintaining their size fairly constant, that drives persistent genes clustering. A further selective stabilization process might contribute to maintaining the clustering.

Extensive genomic diversity and selective conservation of virulence-determinants in enterohemorrhagic Escherichia coli strains of O157 and non-O157 serotypes

BACKGROUND

Enterohemorrhagic Escherichia coli (EHEC) O157 causes severe food-borne illness in humans. The chromosome of O157 consists of 4.1 Mb backbone sequences shared by benign E. coli K-12, and 1.4 Mb O157-specific sequences encoding many virulence determinants, such as Shiga toxin genes (stx genes) and the locus of enterocyte effacement (LEE). Non-O157 EHECs belonging to distinct clonal lineages from O157 also cause similar illness in humans. According to the 'parallel' evolution model, they have independently acquired the major virulence determinants, the stx genes and LEE. However, the genomic differences between O157 and non-O157 EHECs have not yet been systematically analyzed.

RESULTS

Using microarray and whole genome PCR scanning analyses, we performed a whole genome comparison of 20 EHEC strains of O26, O111, and O103 serotypes with O157. In non-O157 EHEC strains, although genome sizes were similar with or rather larger than O157 and the backbone regions were well conserved, O157-specific regions were very poorly conserved. Around only 20% of the O157-specific genes were fully conserved in each non-O157 serotype. However, the non-O157 EHECs contained a significant number of virulence genes that are found on prophages and plasmids in O157, and also multiple prophages similar to, but significantly divergent from, those in O157.

CONCLUSION

Although O157 and non-O157 EHECs have independently acquired a huge amount of serotype- or strain-specific genes by lateral gene transfer, they share an unexpectedly large number of virulence genes. Independent infections of similar but distinct bacteriophages carrying these virulence determinants are deeply involved in the evolution of O157 and non-O157 EHECs.

Silencing of xenogeneic DNA by H-NS--facilitation of lateral gene transfer in bacteria by a defense system that recognizes foreign DNA

Lateral gene transfer has played a prominent role in bacterial evolution, but the mechanisms allowing bacteria to tolerate the acquisition of foreign DNA have been incompletely defined. Recent studies show that H-NS, an abundant nucleoid-associated protein in enteric bacteria and related species, can recognize and selectively silence the expression of foreign DNA with higher adenine and thymine content relative to the resident genome, a property that has made this molecule an almost universal regulator of virulence determinants in enteric bacteria. These and other recent findings challenge the ideas that curvature is the primary determinant recognized by H-NS and that activation of H-NS-silenced genes in response to environmental conditions occurs through a change in the structure of H-NS itself. Derepression of H-NS-silenced genes can occur at specific promoters by several mechanisms including competition with sequence-specific DNA-binding proteins, thereby enabling the regulated expression of foreign genes. The possibility that microorganisms maintain and exploit their characteristic genomic GC ratios for the purpose of self/non-self-discrimination is discussed.

Genome content determination in methicillin-resistant Staphylococcus aureus

Staphylococcus aureus is a major pathogen responsible for both nosocomial and community-acquired infections. While the first S. aureus isolates displaying resistance to methicillin were reported in the early 1960s, endemic strains of methicillin-resistant S. aureus (MRSA) carrying multiple resistance determinants only became a worldwide nosocomial problem in the early 1980s, carrying a threefold attributable cost and a threefold excess length of hospital stay when compared with methicillin-susceptible S. aureus bacteremia. Recent efforts in the field of high-throughput sequencing resulted in the release of several MRSA genome sequences enabling the development of massively parallel tools to study clinical isolates of MRSA at the organism scale. Microarrays covering whole genomes and high-throughput sequencing devices are the two main techniques currently utilizable for whole-genome characterization. These tools not only provide information for the development of genotyping assays but also allow evaluation of potential virulence of the strains, by enumerating genetic-encoded resistance markers and toxin content. This appears particularly attractive for understanding the epidemiology of MRSA and the relationship between genome content on one side and virulence potential or epidemicity on the other side. In addition, sequence information is mandatory for the development of molecular tests allowing the rapid identification, genotyping and characterization of clinical isolates.

Modelling prokaryote gene content

The patchy distribution of genes across the prokaryotes may be caused by multiple gene losses or lateral transfer. Probabilistic models of gene gain and loss are needed to distinguish between these possibilities. Existing models allow only single genes to be gained and lost, despite the empirical evidence for multi-gene events. We compare birth-death models (currently the only widely-used models, in which only one gene can be gained or lost at a time) to blocks models (allowing gain and loss of multiple genes within a family). We analyze two pairs of genomes: two E. coli strains, and the distantly-related Archaeoglobus fulgidus (archaea) and Bacillus subtilis (gram positive bacteria). Blocks models describe the data much better than birth-death models. Our models suggest that lateral transfers of multiple genes from the same family are rare (although transfers of single genes are probably common). For both pairs, the estimated median time that a gene will remain in the genome is not much greater than the time separating the common ancestors of the archaea and bacteria. Deep phylogenetic reconstruction from sequence data will therefore depend on choosing genes likely to remain in the genome for a long time. Phylogenies based on the blocks model are more biologically plausible than phylogenies based on the birth-death model.

Markov models of genome segmentation

We introduce Markov models for segmentation of symbolic sequences, extending a segmentation procedure based on the Jensen-Shannon divergence that has been introduced earlier. Higher-order Markov models are more sensitive to the details of local patterns and in application to genome analysis, this makes it possible to segment a sequence at positions that are biologically meaningful. We show the advantage of higher-order Markov-model-based segmentation procedures in detecting compositional inhomogeneity in chimeric DNA sequences constructed from genomes of diverse species, and in application to the E. coli K12 genome, boundaries of genomic islands, cryptic prophages, and horizontally acquired regions are accurately identified.

The next generation of microarray research: applications in evolutionary and ecological genomics

Microarray technology is one of the key developments in recent years that has propelled biological research into the post-genomic era. With the ability to assay thousands to millions of features at the same time, microarray technology has fundamentally changed how biological questions are addressed, from examining one or a few genes to a collection of genes or the whole genome. This technology has much to offer in the study of genome evolution. After a brief introduction on the technology itself, we then focus on the use of microarrays to examine genome dynamics, to uncover novel functional elements in genomes, to unravel the evolution of regulatory networks, to identify genes important for behavioral and phenotypic plasticity, and to determine microbial community diversity in environmental samples. Although there are still practical issues in using microarrays, they will be alleviated by rapid advances in array technology and analysis methods, the availability of many genome sequences of closely related species and flexibility in array design. It is anticipated that the application of microarray technology will continue to better our understanding of evolution and ecology through the examination of individuals, populations, closely related species or whole microbial communities.

Gene expression analysis indicates extensive genotype-specific crosstalk between the conjugative F-plasmid and the E. coli chromosome

Background

Plasmids are an important component of the bacterial genome, but the crosstalk between genes encoded on the chromosome and on the plasmid is still poorly understood.

Results

We performed a large-scale survey for genes on the E. coli chromosome that are affected by the presence of the conjugative F-plasmid (crosstalk). The expression pattern of about 4% (107 genes) of the genes encoded by the chromosome was affected by the presence of the F-plasmid. Comparing two different Escherichia coli strains, MG1655 and DH5α, we found a strong host genotype-specific crosstalk of the host chromosome with the F-plasmid. About 88% of the genes affected by the presence of the F-plasmid showed a significant plasmid by host genotype interaction, i.e. the presence of the F-plasmid resulted in a different gene expression in the two host genotypes. Less than 12% of the genes showed an additive effect of gene expression, i.e. host genotype independent crosstalk between plasmid and host chromosome.

Conclusion

We propose that epistatic effects also contribute to the maintenance of F-plasmids in natural populations.

Design of a seven-genome Escherichia coli microarray for comparative genomic profiling

We describe the design and evaluate the use of a high-density oligonucleotide microarray covering seven sequenced Escherichia coli genomes in addition to several sequenced E. coli plasmids, bacteriophages, pathogenicity islands, and virulence genes. Its utility is demonstrated for comparative genomic profiling of two unsequenced strains, O175:H16 D1 and O157:H7 3538 (Deltastx(2)::cat) as well as two well-known control strains, K-12 W3110 and O157:H7 EDL933. By using fluorescently labeled genomic DNA to query the microarrays and subsequently analyze common virulence genes and phage elements and perform whole-genome comparisons, we observed that O175:H16 D1 is a K-12-like strain and confirmed that its phi3538 (Deltastx(2)::cat) phage element originated from the E. coli 3538 (Deltastx(2)::cat) strain, with which it shares a substantial proportion of phage elements. Moreover, a number of genes involved in DNA transfer and recombination was identified in both new strains, providing a likely explanation for their capability to transfer phi3538 (Deltastx(2)::cat) between them. Analyses of control samples demonstrated that results using our custom-designed microarray were representative of the true biology, e.g., by confirming the presence of all known chromosomal phage elements as well as 98.8 and 97.7% of queried chromosomal genes for the two control strains. Finally, we demonstrate that use of spatial information, in terms of the physical chromosomal locations of probes, improves the analysis.

The fate of laterally transferred genes: life in the fast lane to adaptation or death

Large-scale genome arrangement plays an important role in bacterial genome evolution. A substantial number of genes can be inserted into, deleted from, or rearranged within genomes during evolution. Detecting or inferring gene insertions/deletions is of interest because such information provides insights into bacterial genome evolution and speciation. However, efficient inference of genome events is difficult because genome comparisons alone do not generally supply enough information to distinguish insertions, deletions, and other rearrangements. In this study, homologous genes from the complete genomes of 13 closely related bacteria were examined. The presence or absence of genes from each genome was cataloged, and a maximum likelihood method was used to infer insertion/deletion rates according to the phylogenetic history of the taxa. It was found that whole gene insertions/deletions in genomes occur at rates comparable to or greater than the rate of nucleotide substitution and that higher insertion/deletion rates are often inferred to be present at the tips of the phylogeny with lower rates on more ancient interior branches. Recently transferred genes are under faster and relaxed evolution compared with more ancient genes. Together, this implies that many of the lineage-specific insertions are lost quickly during evolution and that perhaps a few of the genes inserted by lateral transfer are niche specific.

Over-representation of Chi sequences caused by di-codon increase in Escherichia coli K-12

Chi sequences (5'-GCTGGTGG-3') are cis-acting 8 bp sequence elements that enhance homologous recombination promoted by the RecBCD pathway in Escherichia coli. The genome of E. coli K-12 MG1655 contains 1009 Chi sequences and this frequency far exceeds the expected value for occurrence of an 8 bp sequence in a genome of this size. It is generally thought that the over-representation of Chi sequences indicates that they have been selected for during evolution because of their function in recombination. The genes from three E. coli strains (K-12, O157 and CFT) were classified into three categories (island, match to other E. coli, and backbone). Island genes have a different base composition and codon usage in comparison with those in the backbone genes, therefore they were relatively new and not yet adapted to the base composition patterns and codon usage typical of the recipient genome. The over-representation of Chi sequences was examined by comparing Chi frequencies and codon frequencies between island and backbone genes. The difference in the CTGGTG di-codon frequency between the backbone and island genes was correlated with the frequency of Chi sequences which were translated in the Leu-Val (-G/CTG/GTG/G-) reading frame in the K-12 strain. These results suggest that the main reading frame of Chi sequences increased as a result of the di-codon CTG-GTG increasing under a genome-wide pressure for adapting to the codon usage and base composition of the E. coli K-12 strain, and that the RecBCD recombinase might adjust its recognition sequence to a frequently occurring oligomer such as G-CTG-GTG-G.

A genome rearrangement has orphaned theEscherichia coliK-12 AcpT phosphopantetheinyl transferase from its cognateEscherichia coliO157:H7 substrates

Phosphopantetheinyl transferases (PPTases) are enzymes that catalyse the transfer of a 4'-phosphopantetheine moiety from CoA to a conserved serine residue of a carrier protein. These carrier proteins use the 4'-phosphopantetheine thiol to shuttle intermediates between the active sites of biosynthetic enzymes involved in fatty acid, non-ribosomal peptide and polyketide synthesis. Three PPTases have been previously been identified in Escherichia coli K-12 and other E. coli strains by homology searches and are encoded by the genes acpS, entD and acpT. Both AcpS and EntD have been well studied whereas the function of AcpT has been an enigma because no carrier protein substrate could be found. We report genetic and biochemical evidence that AcpT modifies two carrier proteins encoded in O-island 138, a cluster of fatty acid biosynthesis-like genes located adjacent to acpT in the genome of the pathogenic E. coli strain O157:H7 (E. coli K-12 and several other sequenced E. coli and Shigella strains lack O-island 138). The two carrier proteins of O-island 138 of strain O157:H7 are not modified (or only very poorly modified) by AcpS, the PPTase responsible for 4'-phosphopantetheine attachment to the acyl carrier protein (AcpP) of fatty acid synthesis. We demonstrate that AcpT cannot functionally replace AcpS in E. coli K-12 either in its native chromosomal location or upon insertion of acpT into the acpS chromosomal location. However, in the absence of AcpS activity AcpT does allow very slow growth thus providing a rationale for its retention in the absence of its cognate substrates. These results together with phylogenetic analyses and comparisons of the E. coli and Shigella strains of known genome sequence strongly argue that AcpT has been orphaned from its cognate substrates by a deletion event that occurred in a common ancestor of these organisms. This seems one of the few cases where a chromosomal rearrangement has been functionally demonstrated to be a deletion event rather than an insertion event in the reference organism. We also show that the previously reported suppression of an acpS mutation by the deletion of Lon protease is an artifact of the increased capsular polysaccharide production of lon strains.

Genomic patterns of pathogen evolution revealed by comparison of Burkholderia pseudomallei, the causative agent of melioidosis, to avirulent Burkholderia thailandensis

Background

The Gram-negative bacterium Burkholderia pseudomallei (Bp) is the causative agent of the human disease melioidosis. To understand the evolutionary mechanisms contributing to Bp virulence, we performed a comparative genomic analysis of Bp K96243 and B. thailandensis (Bt) E264, a closely related but avirulent relative.

Results

We found the Bp and Bt genomes to be broadly similar, comprising two highly syntenic chromosomes with comparable numbers of coding regions (CDs), protein family distributions, and horizontally acquired genomic islands, which we experimentally validated to be differentially present in multiple Bt isolates. By examining species-specific genomic regions, we derived molecular explanations for previously-known metabolic differences, discovered potentially new ones, and found that the acquisition of a capsular polysaccharide gene cluster in Bp, a key virulence component, is likely to have occurred non-randomly via replacement of an ancestral polysaccharide cluster. Virulence related genes, in particular members of the Type III secretion needle complex, were collectively more divergent between Bp and Bt compared to the rest of the genome, possibly contributing towards the ability of Bp to infect mammalian hosts. An analysis of pseudogenes between the two species revealed that protein inactivation events were significantly biased towards membrane-associated proteins in Bt and transcription factors in Bp.

Conclusion

Our results suggest that a limited number of horizontal-acquisition events, coupled with the fine-scale functional modulation of existing proteins, are likely to be the major drivers underlying Bp virulence. The extensive genomic similarity between Bp and Bt suggests that, in some cases, Bt could be used as a possible model system for studying certain aspects of Bp behavior.

Co-evolution of the outer surface protein C gene (ospC) and intraspecific lineages of Borrelia burgdorferi sensu stricto in the northeastern United States

Clinical and tick isolates of Borrelia burgdorferi sensu stricto, the bacterial agent of Lyme disease, from the northeastern United States were sequenced at 12 loci located on the main chromosome and 7 plasmids (lp54, cp26, cp9, lp17, lp25, lp28-2, and lp38). The outer surface protein C gene (ospC) showed the highest number (12) of major alleles (defined as alleles differing by 5% or more in nucleotide sequence) while other ORFs had only two to four major alleles. A non-recombining chromosomal marker, the rrs-rrlA ribosomal RNA spacer, was used to infer the intraspecific phylogeny among these B. burgdorferi isolates. We were thus able to analyze the multilocus genotypes in the context of a B. burgdorferi intraspecific phylogeny. Except for ospC, sequence variations at plasmid-borne loci showed broad inconsistency with the intraspecific phylogeny, supporting DNA exchanges mediated by plasmid transfers. The multilocus linkage frequently observed in B. burgdorferi populations is due more likely to a "founder effect" than to a lack of recombination. The exceptional phylogenetic consistency of ospC, in conjunction with its selectively maintained high intraspecific diversity, suggested a dominant role ospC plays in the initiation and maintenance of adaptive differentiation in B. burgdorferi.

On-chip non-equilibrium dissociation curves and dissociation rate constants as methods to assess specificity of oligonucleotide probes

Nucleic acid hybridization serves as backbone for many high-throughput systems for detection, expression analysis, comparative genomics and re-sequencing. Specificity of hybridization between probes and intended targets is always critical. Approaches to ensure and evaluate specificity include use of mismatch probes, obtaining dissociation curves rather than single temperature hybridizations, and comparative hybridizations. In this study, we quantify effects of mismatch type and position on intensity of hybridization signals and provide a new approach based on dissociation rate constants to evaluate specificity of hybridized signals in complex target mixtures. Using an extensive set of 18mer oligonucleotide probes on an in situ synthesized biochip platform, we demonstrate that mismatches in the center of the probe are more discriminating than mismatches toward the extremities of the probe and mismatches toward the attached end are less discriminating than those toward the loose end. The observed destabilizing effect of a mismatch type agreed in general with predictions using the nearest neighbor model. Use of a new parameter, specific dissociation temperature (T_d-w, temperature of maximum specific dissociation rate constant), obtained from probe–target duplex dissociation profiles considerably improved the evaluation of specificity. These results have broad implications for hybridization data obtained from complex mixtures of nucleic acids.

New multiple-deletion method for the Corynebacterium glutamicum genome, using a mutant lox sequence

Due to the difficulty of multiple deletions using the Cre/loxP system, a simple, markerless multiple-deletion method based on a Cre/mutant lox system combining a right-element (RE) mutant lox site with a left-element (LE) mutant lox site was employed for large-scale genome rearrangements in Corynebacterium glutamicum. Eight distinct genomic regions that had been identified previously by comparative analysis of C. glutamicum R and C. glutamicum 13032 genomes were targeted for deletion. By homologous recombination, LE and RE mutant lox sites were integrated at each end of a target region. Highly efficient and accurate deletions between the two chromosomal mutant lox sites in the presence of Cre recombinase were realized. A deletion mutant lacking 190 kb of chromosomal regions, encoding a total of 188 open reading frames (ORFs), was obtained. These deletions represent the largest genomic excisions in C. glutamicum reported to date. Despite the loss of numerous predicted ORFs, the mutant exhibited normal growth under standard laboratory conditions. The Cre/loxP system using a pair of mutant lox sites provides a new, efficient genome rearrangement technique for C. glutamicum. It should facilitate the understanding of genome functions of microorganisms.

(Patho-)Genomics of Escherichia coli

Escherichia coli represents a versatile and diverse enterobacterial species which can be subdivided into (i) nonpathogenic, commensal, (ii) intestinal pathogenic and (iii) extraintestinal pathogenic strains. This classification is mainly based on the presence or absence of DNA regions which are frequently associated with certain pathotypes. In most cases, this genetic information has been horizontally acquired and belongs to the flexible E. coli genome, such as plasmids, bacteriophages and genomic islands. These genomic regions contribute to the rapid evolution of E. coli variants as they are frequently subject to rearrangements, excision and transfer as well as further acquisition of additional DNA thus contributing to the creation of new (pathogenic) variants. Genetic diversity and genome plasticity of E. coli has been underestimated. The accumulating amount of sequence information generated in the era of "genomics" helps to increase our understanding of factors and mechanisms that are involved in diversification of this bacterial species as well as in those that may direct host specificity.

Diversification of DNA sequences in the symbiotic genome of Rhizobium etli

Bacteria of the genus Rhizobium and related genera establish nitrogen-fixing symbioses with the roots of leguminous plants. The genetic elements that participate in the symbiotic process are usually compartmentalized in the genome, either as independent replicons (symbiotic plasmids) or as symbiotic regions or islands in the chromosome. The complete nucleotide sequence of the symbiotic plasmid of Rhizobium etli model strain CFN42, symbiont of the common bean plant, has been reported. To better understand the basis of DNA sequence diversification of this symbiotic compartment, we analyzed the distribution of single-nucleotide polymorphisms in homologous regions from different Rhizobium etli strains. The distribution of polymorphisms is highly asymmetric in each of the different strains, alternating regions containing very few changes with regions harboring an elevated number of substitutions. The regions showing high polymorphism do not correspond with discrete genetic elements and are not the same in the different strains, indicating that they are not hypervariable regions of functional genes. Most interesting, some highly polymorphic regions share exactly the same nucleotide substitutions in more than one strain. Furthermore, in different regions of the symbiotic compartment, different sets of strains share the same substitutions. The data indicate that the majority of nucleotide substitutions are spread in the population by recombination and that the contribution of new mutations to polymorphism is relatively low. We propose that the horizontal transfer of homologous DNA segments among closely related organisms is a major source of genomic diversification.

Multiple large segment deletion method for Corynebacterium glutamicum

A precise and scarless genome excision method, employing the Cre/loxP system in concert with double-strand break (DSB)-stimulated intramolecular recombination was developed. The DSBs were mediated by the restriction endonuclease, I-SceI. It permitted multiple deletions of independent 14-, 43-, and 10-kb-long genomic regions on the Corynebacterium glutamicum genome. Accuracy of deletion was confirmed by the loss of marker genes, PCR, and sequencing of new genome joints. Eleven, 58, and 4 genes were predicted on the 14-, 43-, and 10-kb deleted regions, respectively. Although the resultant mutant lost a total of 67 kb encoding 73 genes, it still exhibited normal growth under standard laboratory conditions. Such a large segment deletion method in which multiple, successive deletions are possible is useful for genome engineering.