Transduction, restriction and recombination patterns in Escherichia coli

Distribution of RecBCD and AddAB recombination-associated genes among bacteria in 33 phyla

Homologous recombination plays key roles in fundamental processes such as recovery from DNA damage and in bacterial horizontal gene transfer, yet there are still open questions about the distribution of recognized components of recombination machinery among bacteria and archaea. RecBCD helicase-nuclease plays a central role in recombination among Gammaproteobacteria like Escherichia coli; while bacteria in other phyla, like the Firmicute Bacillus subtilis, use the related AddAB complex. The activity of at least some of these complexes is controlled by short DNA sequences called crossover hotspot instigator (Chi) sites. When RecBCD or AddAB complexes encounter an autologous Chi site during unwinding, they introduce a nick such that ssDNA with a free end is available to invade another duplex. If homologous DNA is present, RecA-dependent homologous recombination is promoted; if not (or if no autologous Chi site is present) the RecBCD/AddAB complex eventually degrades the DNA. We examined the distribution of recBCD and addAB genes among bacteria, and sought ways to distinguish them unambiguously. We examined bacterial species among 33 phyla, finding some unexpected distribution patterns. RecBCD and addAB are less conserved than recA, with the orthologous recB and addA genes more conserved than the recC or addB genes. We were able to classify RecB vs. AddA and RecC vs. AddB in some bacteria where this had not previously been done. We used logo analysis to identify sequence segments that are conserved, but differ between the RecBC and AddAB proteins, to help future differentiation between members of these two families.

Pandemic fluoroquinolone resistant Escherichia coli clone ST1193 emerged via simultaneous homologous recombinations in 11 gene loci

Global growth in antibiotic resistance is a major social problem. A high level of resistance to fluoroquinolones requires the concurrent presence of at least 3 mutations in the target proteins-2 in DNA gyrase (GyrA) and 1 in topoisomerase IV (ParC), which occur in a stepwise manner. In the Escherichia coli chromosome, the gyrA and parC loci are positioned about 1 Mb away from each other. Here we show that the 3 fluoroquinolone resistance mutations are tightly associated genetically in naturally occurring strains. In the latest pandemic uropathogenic and multidrug-resistant E. coli clonal group ST1193, the mutant variants of gyrA and parC were acquired not by a typical gradual, stepwise evolution but all at once. This happened as part of 11 simultaneous homologous recombination events involving 2 phylogenetically distant strains of E. coli, from an uropathogenic clonal complex ST14 and fluoroquinolone-resistant ST10. The gene exchanges swapped regions between 0.5 and 139 Kb in length (183 Kb total) spread along 976 Kb of chromosomal DNA around and between gyrA and parC loci. As a result, all 3 fluoroquinolone resistance mutations in GyrA and ParC have simultaneously appeared in ST1193. Based on molecular clock estimates, this potentially happened as recently as <12 y ago. Thus, naturally occurring homologous recombination events between 2 strains can involve numerous chromosomal gene locations simultaneously, resulting in the transfer of distant but tightly associated genetic mutations and emergence of a both highly pathogenic and antibiotic-resistant strain with a rapid global spread capability.

The Patchy Distribution of Restriction⁻Modification System Genes and the Conservation of Orphan Methyltransferases in Halobacteria

Restriction⁻modification (RM) systems in bacteria are implicated in multiple biological roles ranging from defense against parasitic genetic elements, to selfish addiction cassettes, and barriers to gene transfer and lineage homogenization. In bacteria, DNA-methylation without cognate restriction also plays important roles in DNA replication, mismatch repair, protein expression, and in biasing DNA uptake. Little is known about archaeal RM systems and DNA methylation. To elucidate further understanding for the role of RM systems and DNA methylation in Archaea, we undertook a survey of the presence of RM system genes and related genes, including orphan DNA methylases, in the halophilic archaeal class Halobacteria. Our results reveal that some orphan DNA methyltransferase genes were highly conserved among lineages indicating an important functional constraint, whereas RM systems demonstrated patchy patterns of presence and absence. This irregular distribution is due to frequent horizontal gene transfer and gene loss, a finding suggesting that the evolution and life cycle of RM systems may be best described as that of a selfish genetic element. A putative target motif (CTAG) of one of the orphan methylases was underrepresented in all of the analyzed genomes, whereas another motif (GATC) was overrepresented in most of the haloarchaeal genomes, particularly in those that encoded the cognate orphan methylase.

Resistance and tolerance to foreign elements by prokaryotic immune systems - curating the genome

To engage in adaptive symbioses or genetic exchange, organisms must interact with foreign, non-self elements despite the risks of predation and parasitism. By surveying the interface between self and non-self, immune systems can help ensure the benevolence of these interactions without isolating their hosts altogether. In this Essay, we examine prokaryotic restriction-modification and CRISPR-Cas (clustered, regularly interspaced palindromic repeat-CRISPR-associated proteins) activities and discuss their analogy to mammalian immune pathways. We further explain how their capacities for resistance and tolerance are optimized to reduce parasitism and immunopathology during encounters with non-self.

Plasmid pEC156, a Naturally Occurring Escherichia coli Genetic Element That Carries Genes of the EcoVIII Restriction-Modification System, Is Mobilizable among Enterobacteria

Type II restriction-modification systems are ubiquitous in prokaryotes. Some of them are present in naturally occurring plasmids, which may facilitate the spread of these systems in bacterial populations by horizontal gene transfer. However, little is known about the routes of their dissemination. As a model to study this, we have chosen an Escherichia coli natural plasmid pEC156 that carries the EcoVIII restriction modification system. The presence of this system as well as the cis-acting cer site involved in resolution of plasmid multimers determines the stable maintenance of pEC156 not only in Escherichia coli but also in other enterobacteria. We have shown that due to the presence of oriT-type F and oriT-type R64 loci it is possible to mobilize pEC156 by conjugative plasmids (F and R64, respectively). The highest mobilization frequency was observed when pEC156-derivatives were transferred between Escherichia coli strains, Enterobacter cloacae and Citrobacter freundii representing coliform bacteria. We found that a pEC156-derivative with a functional EcoVIII restriction-modification system was mobilized in enterobacteria at a frequency lower than a plasmid lacking this system. In addition, we found that bacteria that possess the EcoVIII restriction-modification system can efficiently release plasmid content to the environment. We have shown that E. coli cells can be naturally transformed with pEC156-derivatives, however, with low efficiency. The transformation protocol employed neither involved chemical agents (e.g. CaCl₂) nor temperature shift which could induce plasmid DNA uptake.

Recombinant transfer in the basic genome of Escherichia coli

An approximation to the ∼4-Mbp basic genome shared by 32 strains of Escherichia coli representing six evolutionary groups has been derived and analyzed computationally. A multiple alignment of the 32 complete genome sequences was filtered to remove mobile elements and identify the most reliable ∼90% of the aligned length of each of the resulting 496 basic-genome pairs. Patterns of single base-pair mutations (SNPs) in aligned pairs distinguish clonally inherited regions from regions where either genome has acquired DNA fragments from diverged genomes by homologous recombination since their last common ancestor. Such recombinant transfer is pervasive across the basic genome, mostly between genomes in the same evolutionary group, and generates many unique mosaic patterns. The six least-diverged genome pairs have one or two recombinant transfers of length ∼40-115 kbp (and few if any other transfers), each containing one or more gene clusters known to confer strong selective advantage in some environments. Moderately diverged genome pairs (0.4-1% SNPs) show mosaic patterns of interspersed clonal and recombinant regions of varying lengths throughout the basic genome, whereas more highly diverged pairs within an evolutionary group or pairs between evolutionary groups having >1.3% SNPs have few clonal matches longer than a few kilobase pairs. Many recombinant transfers appear to incorporate fragments of the entering DNA produced by restriction systems of the recipient cell. A simple computational model can closely fit the data. Most recombinant transfers seem likely to be due to generalized transduction by coevolving populations of phages, which could efficiently distribute variability throughout bacterial genomes.

Naturally-occurring, dually-functional fusions between restriction endonucleases and regulatory proteins

Background

Restriction-modification (RM) systems appear to play key roles in modulating gene flow among bacteria and archaea. Because the restriction endonuclease (REase) is potentially lethal to unmethylated new host cells, regulation to ensure pre-expression of the protective DNA methyltransferase (MTase) is essential to the spread of RM genes. This is particularly true for Type IIP RM systems, in which the REase and MTase are separate, independently-active proteins. A substantial subset of Type IIP RM systems are controlled by an activator-repressor called C protein. In these systems, C controls the promoter for its own gene, and for the downstream REase gene that lacks its own promoter. Thus MTase is expressed immediately after the RM genes enter a new cell, while expression of REase is delayed until sufficient C protein accumulates. To study the variation in and evolution of this regulatory mechanism, we searched for RM systems closely related to the well-studied C protein-dependent PvuII RM system. Unexpectedly, among those found were several in which the C protein and REase genes were fused.

Results

The gene for CR.NsoJS138I fusion protein (nsoJS138ICR, from the bacterium Niabella soli) was cloned, and the fusion protein produced and partially purified. Western blots provided no evidence that, under the conditions tested, anything other than full-length fusion protein is produced. This protein had REase activity in vitro and, as expected from the sequence similarity, its specificity was indistinguishable from that for PvuII REase, though the optimal reaction conditions were different. Furthermore, the fusion was active as a C protein, as revealed by in vivo activation of a lacZ reporter fusion to the promoter region for the nsoJS138ICR gene.

Conclusions

Fusions between C proteins and REases have not previously been characterized, though other fusions have (such as between REases and MTases). These results reinforce the evidence for impressive modularity among RM system proteins, and raise important questions about the implications of the C-REase fusions on expression kinetics of these RM systems.

A bistable hysteretic switch in an activator-repressor regulated restriction-modification system

Restriction–modification (RM) systems are extremely widespread among bacteria and archaea, and are often specified by mobile genetic elements. In type II RM systems, where the restriction endonuclease (REase) and protective DNA methyltransferase (MTase) are separate proteins, a major regulatory challenge is delaying expression of the REase relative to the MTase after RM genes enter a new host cell. Basic understanding of this regulation is available for few RM systems, and detailed understanding for none. The PvuII RM system is one of a large subset in which the central regulatory role is played by an activator–repressor protein (called C, for controller). REase expression depends upon activation by C, whereas expression of the MTase does not. Thus delay of REase expression depends on the rate of C-protein accumulation. This is a nonlinear process, as C also activates transcription of its own gene. Mathematical modeling of the PvuII system led to the unexpected predictions of responsiveness to a factor not previously studied in RM system control—gene copy number—and of a hysteretic response. In this study, those predictions have been confirmed experimentally. The results may apply to many other C-regulated RM systems, and help explain their ability to spread so widely.

Diverse functions of restriction-modification systems in addition to cellular defense

Restriction-modification (R-M) systems are ubiquitous and are often considered primitive immune systems in bacteria. Their diversity and prevalence across the prokaryotic kingdom are an indication of their success as a defense mechanism against invading genomes. However, their cellular defense function does not adequately explain the basis for their immaculate specificity in sequence recognition and nonuniform distribution, ranging from none to too many, in diverse species. The present review deals with new developments which provide insights into the roles of these enzymes in other aspects of cellular function. In this review, emphasis is placed on novel hypotheses and various findings that have not yet been dealt with in a critical review. Emerging studies indicate their role in various cellular processes other than host defense, virulence, and even controlling the rate of evolution of the organism. We also discuss how R-M systems could have successfully evolved and be involved in additional cellular portfolios, thereby increasing the relative fitness of their hosts in the population.

Role of DNA methyltransferases in epigenetic regulation in bacteria

In prokaryotes, alteration in gene expression was observed with the modification of DNA, especially DNA methylation. Such changes are inherited from generation to generation with no alterations in the DNA sequence and represent the epigenetic signal in prokaryotes. DNA methyltransferases are enzymes involved in DNA modification and thus in epigenetic regulation of gene expression. DNA methylation not only affects the thermodynamic stability of DNA, but also changes its curvature. Methylation of specific residues on DNA can affect the protein-DNA interactions. DNA methylation in prokaryotes regulates a number of physiological processes in the bacterial cell including transcription, DNA mismatch repair and replication initiation. Significantly, many reports have suggested a role of DNA methylation in regulating the expression of a number of genes in virulence and pathogenesis thus, making DNA methlytransferases novel targets for the designing of therapeutics. Here, we summarize the current knowledge about the influence of DNA methylation on gene regulation in different bacteria, and on bacterial virulence.

Extensive genomic variation within clonal complexes of Neisseria meningitidis

Meningococcal disease is a widely distributed complex disease affecting all age categories. It can cause severe meningitis and septicemia, especially in unvaccinated infants and young children. The causative agent, Neisseria meningitidis (Nm), can be phenotypically and genetically differentiated into serogroups and sequence types (STs) and has a highly dynamic population structure. To obtain a deeper understanding of the epidemiology of Nm, we sequenced seven Nm genomes. Large-scale genomic analysis was conducted with these 7 Nm genomes, 27 additional Nm genomes from GenBank, and 4 other Neisseria genomes. We observed extensive homologous recombination in all gene functional categories among different Nm genomes. Homologous recombination is so frequent that it has resulted in numerous chimeric open reading frames, including genes in the capsule biosynthesis cluster and loci targeted by commercial vaccines. Our results reveal that, despite widespread use, evolutionary relationships inferred from the standard seven-gene multilocus sequence typing (MLST) method could not predict virulence gene content or strain phenotype. In fact, up to 28% of the virulence-associated genes could differ between strains of identical STs. Consistent with previous studies, we found that allelic recombination is also associated with alterations in antibiotic susceptibility. Overall, these findings emphasize the extensive genomic plasticity of Nm and the limitations of standard molecular methods to quantify this genotypic and phenotypic diversity.

Transformation of natural genetic variation into Haemophilus influenzae genomes

Many bacteria are able to efficiently bind and take up double-stranded DNA fragments, and the resulting natural transformation shapes bacterial genomes, transmits antibiotic resistance, and allows escape from immune surveillance. The genomes of many competent pathogens show evidence of extensive historical recombination between lineages, but the actual recombination events have not been well characterized. We used DNA from a clinical isolate of Haemophilus influenzae to transform competent cells of a laboratory strain. To identify which of the ∼40,000 polymorphic differences had recombined into the genomes of four transformed clones, their genomes and their donor and recipient parents were deep sequenced to high coverage. Each clone was found to contain ∼1000 donor polymorphisms in 3–6 contiguous runs (8.1±4.5 kb in length) that collectively comprised ∼1–3% of each transformed chromosome. Seven donor-specific insertions and deletions were also acquired as parts of larger donor segments, but the presence of other structural variation flanking 12 of 32 recombination breakpoints suggested that these often disrupt the progress of recombination events. This is the first genome-wide analysis of chromosomes directly transformed with DNA from a divergent genotype, connecting experimental studies of transformation with the high levels of natural genetic variation found in isolates of the same species.

The ability of bacteria to acquire genetic information from their relatives—called natural competence—poses a major health risk, since recombination between pathogenic bacterial lineages can help bacteria develop resistance to antibiotics and adapt to host defenses. In this study we transformed competent cells of the human pathogen Haemophilus influenzae with genomic DNA from a divergent clinical isolate and used deep sequencing to identify the recombination events in four transformed chromosomes. The results show that transformation of single competent cells is more extensive than expected, and suggests that transformation can be used as a tool to map traits that vary between clinical isolates.

Helicobacter pylori genome evolution during human infection

High genetic diversity is a hallmark of the gastric pathogen Helicobacter pylori. We used 454 sequencing technology to perform whole-genome comparisons for five sets of H. pylori strains that had been sequentially cultured from four chronically infected Colombians (isolation intervals=3-16 y) and one human volunteer experimentally infected with H. pylori as part of a vaccine trial. The four sets of genomes from Colombian H. pylori differed by 27-232 isolated SNPs and 16-441 imported clusters of polymorphisms resulting from recombination. Imports (mean length=394 bp) were distributed nonrandomly over the chromosome and frequently occurred in groups, suggesting that H. pylori first takes up long DNA fragments, which subsequently become partially integrated in multiple shorter pieces. Imports were present at significantly increased frequency in members of the hop family of outer membrane gene paralogues, some of which are involved in bacterial adhesion, suggesting diversifying selection. No evidence of recombination and few other differences were identified in the strain pair from an infected volunteer, indicating that the H. pylori genome is stable in the absence of mixed infection. Among these few differences was an OFF/ON switch in the phase-variable adhesin gene hopZ, suggesting strong in vivo selection for this putative adhesin during early colonization.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

The population genetics of commensal Escherichia coli

The primary habitat of Escherichia coli is the vertebrate gut, where it is the predominant aerobic organism, living in symbiosis with its host. Despite the occurrence of recombination events, the population structure is predominantly clonal, allowing the delineation of major phylogenetic groups. The genetic structure of commensal E. coli is shaped by multiple host and environmental factors, and the determinants involved in the virulence of the bacteria may in fact reflect adaptation to commensal habitats. A better characterization of the commensal niche is necessary to understand how a useful commensal can become a harmful pathogen. In this Review we describe the population structure of commensal E. coli, the factors involved in the spread of different strains, how the bacteria can adapt to different niches and how a commensal lifestyle can evolve into a pathogenic one.

M1.MboII and M2.MboII type IIS methyltransferases: different specificities, the same target

Methylation of a base in a specific DNA sequence protects the DNA from nucleolytic cleavage by restriction enzymes recognizing the same sequence. The MboII restriction-modification (R-M) system of Moraxella bovis ATCC 10900 consists of a restriction endonuclease gene and two methyltransferase genes. The enzymes encoded by this system recognize an asymmetrical sequence 5'-GAAGA-3'/3'-CTTCT-5'. M1.MboII modifies the last adenine in the recognition sequence 5'-GAAGA-3' to N(6)-methyladenine. A second methylase, M2.MboII, was cloned and purified to electrophoretic homogeneity using a four-step chromatographic procedure. It was demonstrated that M2.MboII modifies the internal cytosine in the recognition sequence 3'-CTTCT-5', yielding N(4)-methylcytosine, and moreover is able to methylate single-stranded DNA. The protein exists in solution as a monomer of molecular mass 30 000+/-1000 Da under denaturing conditions. Divalent cations (Ca(2+), Mg(2+), Mn(2+) and Zn(2+)) inhibit M2.MboII methylation activity. It was found that the isomethylomer M2.NcuI from Neisseria cuniculi ATCC 14688 behaves in the same manner. Functional analysis showed that the complete MboII R-M system, consisting of two methyltransferases genes and the mboIIR gene, is the most stable and the least harmful to bacterial cells.

Looking for Darwin's footprints in the microbial world

As we observe the 200th anniversary of Charles Darwin's birth, microbiologists interested in the application of Darwin's ideas to the microscopic world have a lot to celebrate: an emerging picture of the (mostly microbial) Tree of Life at ever-increasing resolution, an understanding of horizontal gene transfer as a driving force in the evolution of microbes, and thousands of complete genome sequences to help formulate and refine our theories. At the same time, quantitative models of the microevolutionary processes shaping microbial populations remain just out of reach, a point that is perhaps most dramatically illustrated by the lack of consensus on how (or even whether) to define bacterial species. Here, we summarize progress and prospects in bacterial population genetics, with an emphasis on detecting the footprint of positive Darwinian selection in microbial genomes.

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.

Tuning the relative affinities for activating and repressing operators of a temporally regulated restriction-modification system

Most type II restriction-modification (R-M) systems produce separate endonuclease (REase) and methyltransferase (MTase) proteins. After R-M genes enter a new cell, MTase activity must appear before REase or the host chromosome will be cleaved. Temporal control of these genes thus has life-or-death consequences. PvuII and some other R-M systems delay endonuclease expression by cotranscribing the REase gene with the upstream gene for an autogenous activator/repressor (C protein). C.PvuII was previously shown to have low levels early, but positive feedback later boosts transcription of the C and REase genes. The MTase is expressed without delay, and protects the host DNA. C.PvuII binds to two sites upstream of its gene: O_L, associated with activation, and O_R, associated with repression. Even when symmetry elements of each operator are made identical, C.PvuII binds preferentially to O_L. In this study, the intra-operator spacers are shown to modulate relative C.PvuII affinity. In light of a recently reported C.Esp1396I-DNA co-crystal structure, in vitro and in vivo effects of altering O_L and O_R spacers were determined. The results suggest that the GACTnnnAGTC consensus is the primary determinant of C.PvuII binding affinity, with intra-operator spacers playing a fine-tuning role that affects mobility of this R-M system.

Mosaic DNA imports with interspersions of recipient sequence after natural transformation of Helicobacter pylori

Helicobacter pylori colonizes the gastric mucosa of half of the human population, causing gastritis, ulcers, and cancer. H. pylori is naturally competent for transformation by exogenous DNA, and recombination during mixed infections of one stomach with multiple H. pylori strains generates extensive allelic diversity. We developed an in vitro transformation protocol to study genomic imports after natural transformation of H. pylori. The mean length of imported fragments was dependent on the combination of donor and recipient strain and varied between 1294 bp and 3853 bp. In about 10% of recombinant clones, the imported fragments of donor DNA were interrupted by short interspersed sequences of the recipient (ISR) with a mean length of 82 bp. 18 candidate genes were inactivated in order to identify genes involved in the control of import length and generation of ISR. Inactivation of the antimutator glycosylase MutY increased the length of imports, but did not have a significant effect on ISR frequency. Overexpression of mutY strongly increased the frequency of ISR, indicating that MutY, while not indispensable for ISR formation, is part of at least one ISR-generating pathway. The formation of ISR in H. pylori increases allelic diversity, and contributes to the uniquely low linkage disequilibrium characteristic of this pathogen.

Population genomic analysis of strain variation in Leptospirillum group II bacteria involved in acid mine drainage formation

Deeply sampled community genomic (metagenomic) datasets enable comprehensive analysis of heterogeneity in natural microbial populations. In this study, we used sequence data obtained from the dominant member of a low-diversity natural chemoautotrophic microbial community to determine how coexisting closely related individuals differ from each other in terms of gene sequence and gene content, and to uncover evidence of evolutionary processes that occur over short timescales. DNA sequence obtained from an acid mine drainage biofilm was reconstructed, taking into account the effects of strain variation, to generate a nearly complete genome tiling path for a Leptospirillum group II species closely related to L. ferriphilum (sampling depth ∼20×). The population is dominated by one sequence type, yet we detected evidence for relatively abundant variants (>99.5% sequence identity to the dominant type) at multiple loci, and a few rare variants. Blocks of other Leptospirillum group II types (∼94% sequence identity) have recombined into one or more variants. Variant blocks of both types are more numerous near the origin of replication. Heterogeneity in genetic potential within the population arises from localized variation in gene content, typically focused in integrated plasmid/phage-like regions. Some laterally transferred gene blocks encode physiologically important genes, including quorum-sensing genes of the LuxIR system. Overall, results suggest inter- and intrapopulation genetic exchange involving distinct parental genome types and implicate gain and loss of phage and plasmid genes in recent evolution of this Leptospirillum group II population. Population genetic analyses of single nucleotide polymorphisms indicate variation between closely related strains is not maintained by positive selection, suggesting that these regions do not represent adaptive differences between strains. Thus, the most likely explanation for the observed patterns of polymorphism is divergence of ancestral strains due to geographic isolation, followed by mixing and subsequent recombination.

Communities of microbes in nature consist of a large number of distinct individuals. The variation in DNA sequence between these individuals contains a record of the evolutionary processes that have shaped each community. In most environments, however, the high number of distinct species makes obtaining information about the nature of this variation difficult or impossible. We obtained large amounts of sequence data for a natural community in an acid mine drainage system consisting of only a few species. This enabled us to reconstruct the genome of the dominant bacterium (Leptospirillum group II) and obtain detailed information about sequence variation between individuals, including differences in both gene content and gene sequence. Our analysis shows extensive recombination between closely related populations, as well as fewer instances of recombination between more distantly related individuals. Additionally, viruses and plasmids account for high variability in gene content between individuals. We conclude that sequence-level variation in this population is maintained through neutral processes (migration, recombination, and genetic drift) rather than natural selection. This suggests that closely related strains of the Leptospirillum group II population may not be ecologically distinct.

Deep sequencing of a low-complexity microbial community revealed extensive recombination as well as polymorphic and gene content variation between individuals of the dominant organism. We show that strains defined by linked polymorphisms are not maintained by positive selection; instead, they are predominantly maintained by the forces of migration and drift.

The Fidelity Index provides a systematic quantitation of star activity of DNA restriction endonucleases

Restriction endonucleases are the basic tools of molecular biology. Many restriction endonucleases show relaxed sequence recognition, called star activity, as an inherent property under various digestion conditions including the optimal ones. To quantify this property we propose the concept of the Fidelity Index (FI), which is defined as the ratio of the maximum enzyme amount showing no star activity to the minimum amount needed for complete digestion at the cognate recognition site for any particular restriction endonuclease. Fidelity indices for a large number of restriction endonucleases are reported here. The effects of reaction vessel, reaction volume, incubation mode, substrate differences, reaction time, reaction temperature and additional glycerol, DMSO, ethanol and Mn(2+) on the FI are also investigated. The FI provides a practical guideline for the use of restriction endonucleases and defines a fundamental property by which restriction endonucleases can be characterized.

Molecular characterization of the DNA methyltransferase M1.NcuI from Neisseria cuniculi ATCC 14688

The methyltransferase M1.NcuI is a member of the restriction-modification system in Neisseria cuniculi ATCC14688 and recognizes the asymmetric pentanucleotide sequence 5'-GAAGA-3'/3'-CTTCT-5'. We purified M1.NcuI to electrophoretic homogeneity using a four-step chromatographic procedure. M1.NcuI is a protein with M(r)=32,000+/-1000 under denaturing conditions. It modifies the recognition sequence by transferring the methyl group from S-adenosyl-l-methionine to the 3' adenine of the pentanucleotide sequence 5'-GAAGA-3'. M1.NcuI, like many other methyltransferases, occurs as a monomer in solution, as determined by gel filtration. Divalent cations inhibit the methylation activity of M1.NcuI. Optimal enzyme activity was observed at a pH of 8.0. M1.NcuI cross-reacted with anti-M1.MboII serum which reflects the similarity of M1.NcuI with M1.MboII at the amino acid level. The gene coding for the enzyme, designated ncuIM1, was cloned, sequenced and overexpressed in Escherichia coli. The structural gene is 780 nucleotides in length coding for a protein of 259 amino acids (M(r) 30,098). The presence and distribution of nine highly conserved amino acid sequence motifs and a putative target recognition domain in the enzyme structure suggest that M1.NcuI, similar to M1.MboII and M1.HpyAII, belongs to N(6)-adenine beta-class DNA methyltransferases.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

N6-methyl-adenine: an epigenetic signal for DNA-protein interactions

N(6)-methyl-adenine is found in the genomes of bacteria, archaea, protists and fungi. Most bacterial DNA adenine methyltransferases are part of restriction-modification systems. Certain groups of Proteobacteria also harbour solitary DNA adenine methyltransferases that provide signals for DNA-protein interactions. In gamma-proteobacteria, Dam methylation regulates chromosome replication, nucleoid segregation, DNA repair, transposition of insertion elements and transcription of specific genes. In Salmonella, Haemophilus, Yersinia and Vibrio species and in pathogenic Escherichia coli, Dam methylation is required for virulence. In alpha-proteobacteria, CcrM methylation regulates the cell cycle in Caulobacter, Rhizobium and Agrobacterium, and has a role in Brucella abortus infection.

Nature of the promoter activated by C.PvuII, an unusual regulatory protein conserved among restriction-modification systems

A widely distributed family of small regulators, called C proteins, controls a subset of restriction-modification systems. The C proteins studied to date activate transcription of their own genes and that of downstream endonuclease genes; this arrangement appears to delay endonuclease expression relative to that of the protective methyltransferase when the genes enter a new cell. C proteins bind to conserved sequences called C boxes. In the PvuII system, the C boxes have been reported to extend from -23 to +3 relative to the transcription start for the gene for the C protein, an unexpected starting position relative to a bound activator. This study suggests that transcript initiation within the C boxes represents initial, C-independent transcription of pvuIICR. The major C protein-dependent transcript appears to be a leaderless mRNA starting farther downstream, at the initiation codon for the pvuIIC gene. This conclusion is based on nuclease S1 transcript mapping and the effects of a series of nested deletions in the promoter region. Furthermore, replacing the region upstream of the pvuIIC initiation codon with a library of random oligonucleotides, followed by selection for C-dependent transcription, yielded clones having sequences that resemble -10 promoter hexamers. The -35 hexamer of this promoter would lie within the C boxes. However, the spacing between C boxes/-35 and the apparent -10 hexamer can be varied by +/-4 bp with little effect. This suggests that, like some other activator-dependent promoters, PpvuIICR may not require a -35 hexamer. Features of this transcription activation system suggest explanations for its broad host range.

Type I restriction enzyme with RecA protein promotes illegitimate recombination

Illegitimate (non-homologous) recombination requires little or no sequence homology between recombining DNAs and has been regarded as being a process distinct from homologous recombination, which requires a long stretch of homology between recombining DNAs. However, we have found a type of illegitimate recombination that requires an interaction between long homologous DNA sequences. It was detected when a plasmid that carried 2-kb-long inverted repeats was subjected to type I (EcoKI) restriction in vivo within a special mutant strain of Escherichia coli. In the present work, we analyzed genetic requirements for this type of illegitimate recombination in well-defined genetic backgrounds. Our analysis demonstrated dependence on RecA function and on the presence of two EcoKI sites on the substrate DNA. These results are in harmony with a model in which EcoKI restriction enzyme attacks an intermediate of homologous recombination to divert it to illegitimate recombination.

Characterization of the LlaCI methyltransferase from Lactococcus lactis subsp. cremoris W15 provides new insights into the biology of type II restriction-modification systems

The gene encoding the LlaCI methyltransferase (M.LlaCI) from Lactococcus lactis subsp. cremoris W15 was overexpressed in Escherichia coli. The enzyme was purified to apparent homogeneity using three consecutive steps of chromatography on phosphocellulose, blue-agarose and Superose 12HR, yielding a protein of M(r) 31 300+/-1000 under denaturing conditions. The exact position of the start codon AUG was determined by protein microsequencing. This enzyme recognizes the specific palindromic sequence 5'-AAGCTT-3'. Purified M.LlaCI was characterized. Unlike many other methyltransferases, M.LlaCI exists in solution predominantly as a dimer. It modifies the first adenine residue at the 5' end of the specific sequence to N(6)-methyladenine and thus is functionally identical to the corresponding methyltransferases of the HindIII (Haemophilus influenzae Rd) and EcoVIII (Escherichia coli E1585-68) restriction-modification systems. This is reflected in the identity of M.LlaCI with M.HindIII and M.EcoVIII noted at the amino acid sequence level (50 % and 62 %, respectively) and in the presence of nine sequence motifs conserved among N(6)-adenine beta-class methyltransferases. However, polyclonal antibodies raised against M.EcoVIII cross-reacted with M.LlaCI but not with M.HindIII. Restriction endonucleases require Mg(2+) for phosphodiester bond cleavage. Mg(2+) was shown to be a strong inhibitor of the M.LlaCI enzyme and its isospecific homologues. This observation suggests that sensitivity of the M.LlaCI to Mg(2+) may strengthen the restriction activity of the cognate endonuclease in the bacterial cell. Other biological implications of this finding are also discussed.

Maintenance of species identity and controlling speciation of bacteria: a new function for restriction/modification systems?

Bacteria frequently exchange DNA among each other by horizontal gene transfer. However, maintenance of species identity and in particular speciation requires a certain barrier against an unregulated uptake of foreign DNA. Here it is suggested that formation of such a barrier is one important biological function of restriction/modification systems, in addition to the classical function of protection of bacteria against bacteriophage infection. This model explains the extreme variability and wide distribution of restriction/modification systems among prokaryotes, the prevalence of RM-systems in pathogenic bacteria and the existence of several RM-systems in single bacterial strains.

A molecular phylogeny of enteric bacteria and implications for a bacterial species concept

A molecular phylogeny for seven taxa of enteric bacteria (Citrobacter freundii, Enterobacter cloacae, Escherichia coli, Hafnia alvei, Klebsiella oxytoca, Klebsiella pneumoniae, and Serratia plymuthica) was made from multiple isolates per taxa taken from a collection of environmental enteric bacteria. Sequences from five housekeeping genes (gapA, groEL, gyrA, ompA, and pgi) and the 16S rRNA gene were used to infer individual gene trees and were concatenated to infer a composite molecular phylogeny for the species. The isolates from each taxa formed tight species clusters in the individual gene trees, suggesting the existence of 'genotypic' clusters that correspond to traditional species designations. These sequence data and the resulting gene trees and consensus tree provide the first data set with which to assess the utility of the recently proposed core genome hypothesis (CGH). The CGH provides a genetically based approach to applying the biological species concept to bacteria.

DNA restriction is a barrier to natural transformation in Pseudomonas stutzeri JM300

Natural transformation is a mechanism for intra- and interspecific transfer of chromosomal DNA in Pseudomonas stutzeri. During this process a single strand derived from duplex DNA is transported into the cytoplasm and recombined with resident DNA. By electroporation, which introduces duplex DNA into cells, 100-fold lower transformation frequencies of P. stutzeri JM300 were observed with shuttle vector or broad-host-range plasmid DNA when the plasmids had replicated in Escherichia coli and not in P. stutzeri JM300. Moreover, the natural transformation with cloned chromosomal P. stutzeri JM300 DNA was reduced about 40-fold when the DNA had not been propagated in P. stutzeri JM300 but in E. coli. Restriction was also active during natural transformation by single-stranded DNA. Restriction during natural transformation and electroporation was abolished in mutants isolated from mutagenized JM300 cells after applying a multiple plasmid electroporation strategy for the enrichment of restriction-defective strains. The mutants had retained the ability for DNA modification. The P. stutzeri strain ATCC 17587 was found to have no restriction-modification system as seen in JM300. It is discussed whether restriction during natural transformation acts at presynaptic or postsynaptic stages of transforming DNA. Restriction as a barrier to transformation apparently contributes to sexual isolation and therefore may promote speciation in the highly diverse species P. stutzeri.

Molecular evolution of the Escherichia coli chromosome. VI. Two regions of high effective recombination

Two 6- to 8-min regions, centered respectively near 45 min (O-antigen region) and 99 min (restriction-modification region) on the Escherichia coli chromosome, display unusually high variability among 11 otherwise very similar strains. This variation, revealed by restriction fragment length polymorphism (RFLP) and nucleotide sequence comparisons, appears to be due to a great local increase in the retention frequency of recombinant replacements. We infer a two-step mechanism. The first step is the acquisition of a small stretch of DNA from a phylogenetically distant source. The second is the successful retransmission of the imported DNA, together with flanking native DNA, to other strains of E. coli. Each cell containing the newly transferred DNA has a very high selective advantage until it reaches a high frequency and (in the O-antigen case) is recognized by the new host's immune system. A high selective advantage increases the probability of retention greatly; the effective recombination rate is the product of the basic recombination rate and the probability of retention. Nearby nucleotide sequences clockwise from the O-antigen (rfb) region are correlated with specific O antigens, confirming local hitchhiking. Comparable selection involving imported restriction endonuclease genes is proposed for the region near 99 min.

What are bacterial species?

Bacterial systematics has not yet reached a consensus for defining the fundamental unit of biological diversity, the species. The past half-century of bacterial systematics has been characterized by improvements in methods for demarcating species as phenotypic and genetic clusters, but species demarcation has not been guided by a theory-based concept of species. Eukaryote systematists have developed a universal concept of species: A species is a group of organisms whose divergence is capped by a force of cohesion; divergence between different species is irreversible; and different species are ecologically distinct. In the case of bacteria, these universal properties are held not by the named species of systematics but by ecotypes. These are populations of organisms occupying the same ecological niche, whose divergence is purged recurrently by natural selection. These ecotypes can be discovered by several universal sequence-based approaches. These molecular methods suggest that a typical named species contains many ecotypes, each with the universal attributes of species. A named bacterial species is thus more like a genus than a species.

Prokaryotic evolution in light of gene transfer

Accumulating prokaryotic gene and genome sequences reveal that the exchange of genetic information through both homology-dependent recombination and horizontal (lateral) gene transfer (HGT) is far more important, in quantity and quality, than hitherto imagined. The traditional view, that prokaryotic evolution can be understood primarily in terms of clonal divergence and periodic selection, must be augmented to embrace gene exchange as a creative force, itself responsible for much of the pattern of similarities and differences we see between prokaryotic microbes. Rather than replacing periodic selection on genetic diversity, gene loss, and other chromosomal alterations as important players in adaptive evolution, gene exchange acts in concert with these processes to provide a rich explanatory paradigm-some of whose implications we explore here. In particular, we discuss (1) the role of recombination and HGT in giving phenotypic "coherence" to prokaryotic taxa at all levels of inclusiveness, (2) the implications of these processes for the reconstruction and meaning of "phylogeny," and (3) new views of prokaryotic adaptation and diversification based on gene acquisition and exchange.

Gene transfer in bacteria: speciation without species?

Although Bacteria and Archaea reproduce by binary fission, exchange of genes among lineages has shaped the diversity of their populations and the diversification of their lineages. Gene exchange can occur by two distinct routes, each differentially impacting the recipient genome. First, homologous recombination mediates the exchange of DNA between closely related individuals (those whose sequences are sufficient similarly to allow efficient integration). As a result, homologous recombination mediates the dispersal of advantageous alleles that may rise to high frequency among genetically related individuals via periodic selection events. Second, lateral gene transfer can introduce novel DNA into a genome from completely unrelated lineages via illegitimate recombination. Gene exchange by this route serves to distribute genes throughout distantly related clades and therefore may confer complex abilities--not otherwise found among closely related lineages--onto the recipient organisms. These two mechanisms of gene exchange play complementary roles in the diversification of microbial populations into independent, ecologically distinct lineages. Although the delineation of microbial "species" then becomes difficult--if not impossible--to achieve, a cogent process of speciation can be predicted.

Recombination and mutation during long-term gastric colonization by Helicobacter pylori: estimates of clock rates, recombination size, and minimal age

The bacterium Helicobacter pylori colonizes the gastric mucosa of half of the human population, resulting in chronic gastritis, ulcers, and cancer. We sequenced ten gene fragments from pairs of strains isolated sequentially at a mean interval of 1.8 years from 26 individuals. Several isolates had acquired small mosaic segments from other H. pylori or point mutations. The maximal mutation rate, the import size, and the frequency of recombination were calculated by using a Bayesian model. The calculations indicate that the last common ancestor of H. pylori existed at least 2,500-11,000 years ago. Imported mosaics have a median size of 417 bp, much smaller than for other bacteria, and recombination occurs frequently (60 imports spanning 25,000 bp per genome per year). Thus, the panmictic population structure of H. pylori results from very frequent recombination during mixed colonization by unrelated strains.

Structure and function of type II restriction endonucleases

More than 3000 type II restriction endonucleases have been discovered. They recognize short, usually palindromic, sequences of 4-8 bp and, in the presence of Mg(2+), cleave the DNA within or in close proximity to the recognition sequence. The orthodox type II enzymes are homodimers which recognize palindromic sites. Depending on particular features subtypes are classified. All structures of restriction enzymes show a common structural core comprising four beta-strands and one alpha-helix. Furthermore, two families of enzymes can be distinguished which are structurally very similar (EcoRI-like enzymes and EcoRV-like enzymes). Like other DNA binding proteins, restriction enzymes are capable of non-specific DNA binding, which is the prerequisite for efficient target site location by facilitated diffusion. Non-specific binding usually does not involve interactions with the bases but only with the DNA backbone. In contrast, specific binding is characterized by an intimate interplay between direct (interaction with the bases) and indirect (interaction with the backbone) readout. Typically approximately 15-20 hydrogen bonds are formed between a dimeric restriction enzyme and the bases of the recognition sequence, in addition to numerous van der Waals contacts to the bases and hydrogen bonds to the backbone, which may also be water mediated. The recognition process triggers large conformational changes of the enzyme and the DNA, which lead to the activation of the catalytic centers. In many restriction enzymes the catalytic centers, one in each subunit, are represented by the PD. D/EXK motif, in which the two carboxylates are responsible for Mg(2+) binding, the essential cofactor for the great majority of enzymes. The precise mechanism of cleavage has not yet been established for any enzyme, the main uncertainty concerns the number of Mg(2+) ions directly involved in cleavage. Cleavage in the two strands usually occurs in a concerted fashion and leads to inversion of configuration at the phosphorus. The products of the reaction are DNA fragments with a 3'-OH and a 5'-phosphate.

Antibiotic resistance in the ECOR collection: integrons and identification of a novel aad gene

The 72 Escherichia coli strains of the ECOR collection were examined for resistance to 10 different antimicrobial agents including ampicillin, tetracycline, mercury, trimethoprim, and sulfonamides. Eighteen strains were resistant to at least one of the antibiotics tested, and nearly 20% (14 of 72) were resistant to two or more. Several of the resistance determinants were shown to be carried on conjugative elements. The collection was screened for the presence of the three classes of integrons and for the sul1 gene, which is generally associated with class 1 integrons. The four strains found to carry a class 1 integron also had Tn21-encoded mercury resistance. One of the integrons encoded a novel streptomycin resistance gene, aadA7, with an attC site (or 59-base element) nearly identical to the attC site associated with the qacF gene cassette found in In40 (M.-C. Ploy, P. Courvalin, and T. Lambert, Antimicrob. Agents Chemother. 42:2557-2563, 1998). The conservation of associated attC sites among unrelated resistance cassettes is similar to arrangements found in the Vibrio cholerae superintegrons (D. Mazel, B. Dychinco, V. A. Webb, and J. Davies, Science 280:605-608, 1998) and supports the hypothesis that resistance cassettes are picked up from superintegron pools and independently assembled from unrelated genes and related attC sites.

Type I restriction systems: sophisticated molecular machines (a legacy of Bertani and Weigle)

Restriction enzymes are well known as reagents widely used by molecular biologists for genetic manipulation and analysis, but these reagents represent only one class (type II) of a wider range of enzymes that recognize specific nucleotide sequences in DNA molecules and detect the provenance of the DNA on the basis of specific modifications to their target sequence. Type I restriction and modification (R-M) systems are complex; a single multifunctional enzyme can respond to the modification state of its target sequence with the alternative activities of modification or restriction. In the absence of DNA modification, a type I R-M enzyme behaves like a molecular motor, translocating vast stretches of DNA towards itself before eventually breaking the DNA molecule. These sophisticated enzymes are the focus of this review, which will emphasize those aspects that give insights into more general problems of molecular and microbial biology. Current molecular experiments explore target recognition, intramolecular communication, and enzyme activities, including DNA translocation. Type I R-M systems are notable for their ability to evolve new specificities, even in laboratory cultures. This observation raises the important question of how bacteria protect their chromosomes from destruction by newly acquired restriction specifities. Recent experiments demonstrate proteolytic mechanisms by which cells avoid DNA breakage by a type I R-M system whenever their chromosomal DNA acquires unmodified target sequences. Finally, the review will reflect the present impact of genomic sequences on a field that has previously derived information almost exclusively from the analysis of bacteria commonly studied in the laboratory.

Cellular responses to postsegregational killing by restriction-modification genes

Plasmids that carry one of several type II restriction modification gene complexes are known to show increased stability. The underlying mechanism was proposed to be the lethal attack by restriction enzyme at chromosomal recognition sites in cells that had lost the restriction modification gene complex. In order to examine bacterial responses to this postsegregational cell killing, we analyzed the cellular processes following loss of the EcoRI restriction modification gene complex carried by a temperature-sensitive plasmid in an Escherichia coli strain that is wild type with respect to DNA repair. A shift to the nonpermissive temperature blocked plasmid replication, reduced the increase in viable cell counts and resulted in loss of cell viability. Many cells formed long filaments, some of which were multinucleated and others anucleated. In a mutant defective in RecBCD exonuclease/recombinase, these cell death symptoms were more severe and cleaved chromosomes accumulated. Growth inhibition was also more severe in recA, ruvAB, ruvC, recG, and recN mutants. The cells induced the SOS response in a RecBC-dependent manner. These observations strongly suggest that bacterial cells die as a result of chromosome cleavage after loss of a restriction modification gene complex and that the bacterial RecBCD/RecA machinery helps the cells to survive, at least to some extent, by repairing the cleaved chromosomes. These and previous results have led us to hypothesize that the RecBCD/Chi/RecA system serves to destroy restricted "nonself" DNA and repair restricted "self" DNA.

Molecular evolution of the Escherichia coli chromosome. V. Recombination patterns among strains of diverse origin

Incorporation patterns of donor DNA into recipient chromosomes following transduction or conjugation have been studied in the progeny of a variety of Escherichia coli crosses in which donor and recipient nucleotide sequences differ by 1-3%. Series of contiguous or variously spaced PCR fragments have been amplified from each recombinant chromosome and digested with a commercial restriction endonuclease previously shown to distinguish the respective parents in a given fragment. We conclude that entering donor DNA fragments are frequently abridged (cut and shortened) before incorporation, the cutting being due to restriction systems, and the shortening presumably due to exonuclease activity. Analysis of several backcrosses confirms, and extends to conjugation, the importance of restriction in E. coli recombination in nature. The transmission patterns in conjugation are similar to those of transduction, but (as expected) on a much larger scale. Asymmetric results of reciprocal crosses imply that mismatch frequency is not a major factor. Marked differences among the results of simple crosses according to parental strain combinations are consistent with observations that E. coli strains in nature vary dramatically in their restriction-modification systems.

PCR-based subtractive hybridization and differences in gene content among strains of Helicobacter pylori

Genes that are characteristic of only certain strains of a bacterial species can be of great biologic interest. Here we describe a PCR-based subtractive hybridization method for efficiently detecting such DNAs and apply it to the gastric pathogen Helicobacter pylori. Eighteen DNAs specific to a monkey-colonizing strain (J166) were obtained by subtractive hybridization against an unrelated strain whose genome has been fully sequenced (26695). Seven J166-specific clones had no DNA sequence match to the 26695 genome, and 11 other clones were mixed, with adjacent patches that did and did not match any sequences in 26695. At the protein level, seven clones had homology to putative DNA restriction-modification enzymes, and two had homology to putative metabolic enzymes. Nine others had no database match with proteins of assigned function. PCR tests of 13 unrelated H. pylori strains by using primers specific for 12 subtracted clones and complementary Southern blot hybridizations indicated that these DNAs are highly polymorphic in the H. pylori population, with each strain yielding a different pattern of gene-specific PCR amplification. The search for polymorphic DNAs, as described here, should help identify previously unknown virulence genes in pathogens and provide new insights into microbial genetic diversity and evolution.

The restriction-modification genes of Escherichia coli K-12 may not be selfish: they do not resist loss and are readily replaced by alleles conferring different specificities

Type II restriction and modification (R-M) genes have been described as selfish because they have been shown to impose selection for the maintenance of the plasmid that encodes them. In our experiments, the type I R-M system EcoKI does not behave in the same way. The genes specifying EcoKI are, however, normally residents of the chromosome and therefore our analyses were extended to monitor the deletion of chromosomal genes rather than loss of plasmid vector. If EcoKI were to behave in the same way as the plasmid-encoded type II R-M systems, the loss of the relevant chromosomal genes by mutation or recombination should lead to cell death because the cell would become deficient in modification enzyme and the bacterial chromosome would be vulnerable to the restriction endonuclease. Our data contradict this prediction; they reveal that functional type I R-M genes in the chromosome are readily replaced by mutant alleles and by alleles encoding a type I R-M system of different specificity. The acquisition of allelic genes conferring a new sequence specificity, but not the loss of the resident genes, is dependent on the product of an unlinked gene, one predicted [Prakash-Cheng, A., Chung, S. S. & Ryu, J. (1993) Mol. Gen. Genet. 241, 491-496] to be relevant to control of expression of the genes that encode EcoKI. Our evidence suggests that not all R-M systems are evolving as "selfish" units; rather, the diversity and distribution of the family of type I enzymes we have investigated require an alternative selective pressure.

E. coli translation initiation factor IF2--an extremely conserved protein. Comparative sequence analysis of the infB gene in clinical isolates of E. coli

The functionally uncharacterised N-terminal of translation initiation factor IF2 has been found to be extremely variable when comparing different bacterial species. In order to study the intraspecies variability of IF2 the 2670 basepairs nucleotide sequence of the infB gene (encoding IF2) was determined in 10 clinical isolates of E. coli. The N-terminal domains (I, II and III) were completely conserved indicating a specific function of this region of IF2. Only one polymorphic position was found in the deduced 890 amino acid sequence. This Gln/Gly490 is located within the central GTP/GDP-binding domain IV of IF2. The results are further evidence that IF2 from E. coli has reached a highly defined level of structural and functional development.

A new type of illegitimate recombination is dependent on restriction and homologous interaction

Illegitimate (nonhomologous) recombination requires little or no sequence homology between recombining DNAs and has been regarded as being a process distinct from homologous recombination, which requires a long stretch of homology between recombining DNAs. Under special conditions in Escherichia coli, we have found a new type of illegitimate recombination that requires an interaction between homologous DNA sequences. It was detected when a plasmid that carried 2-kb-long inverted repeats was subjected to type II restriction in vitro and type I (EcoKI) restriction in vivo within a delta rac recBC recG ruvC strain. Removal of one of the repeats or its replacement with heterologous DNA resulted in a reduction in the level of recombination. The recombining sites themselves shared, at most, a few base pairs of homology. Many of the recombination events joined a site in one of the repeats with a site in another repeat. In two of the products, one of the recombining sites was at the end of one of the repeats. Removal of one of the EcoKI sites resulted in decreased recombination. We discuss the possibility that some structure made by homologous interaction between the long repeats is used by the EcoKI restriction enzyme to promote illegitimate recombination. The possible roles and consequences of this type of homologous interaction are discussed.

A coalescent estimator of the population recombination rate

Population genetic models often use a population recombination parameter 4Nc, where N is the effective population size and c is the recombination rate per generation. In many ways 4Nc is comparable to 4Nu, the population mutation rate. Both combine genome level and population level processes, and together they describe the rate of production of genetic variation in a population. However, 4Nc is more difficult to estimate. For a population sample of DNA sequences, historical recombination can only be detected if polymorphisms exist, and even then most recombination events are not detectable. This paper describes an estimator of 4Nc, hereafter designated gamma (gamma), that was developed using a coalescent model for a sample of four DNA sequences with recombination. The reliability of gamma was assessed using multiple coalescent simulations. In general gamma has low to moderate bias, and the reliability of gamma is comparable, though less, than that for a widely used estimator of 4Nu. If there exists an independent estimate of the recombination rate (per generation, per base pair), gamma can be used to estimate the effective population size or the neutral mutation rate.