Strain belonging to an emerging, virulent sublineage of ST131 Escherichia coli isolated in fresh spinach, suggesting that ST131 may be transmissible through agricultural products.

Food contamination with pathogenic Escherichia coli can cause severe disease. Here, we report the isolation of a multidrug resistant strain (A23EC) from fresh spinach. A23EC belongs to subclade C2 of ST131, a virulent clone of Extraintestinal Pathogenic E. coli (ExPEC). Most A23EC virulence factors are concentrated in three pathogenicity islands. These include PapGII, a fimbrial tip adhesin linked to increased virulence, and CsgA and CsgB, two adhesins known to facilitate spinach leaf colonization. A23EC also bears TnMB1860, a chromosomally-integrated transposon with the demonstrated potential to facilitate the evolution of carbapenem resistance among non-carbapenemase-producing enterobacterales. This transposon consists of two IS26-bound modular translocatable units (TUs). The first TU carries aac(6')-lb-cr, bla OXA-1, ΔcatB3, aac(3)-lle, and tmrB, and the second one harbors bla CXT-M-15. A23EC also bears a self-transmissible plasmid that can mediate conjugation at 20°C and that has a mosaic IncF [F(31,36):A(4,20):B1] and Col156 origin of replication. Comparing A23EC to 86 additional complete ST131 sequences, A23EC forms a monophyletic cluster with 17 other strains that share the following four genomic traits: (1) virotype E (papGII+); (2) presence of a PAI II536-like pathogenicity island with an additional cnf1 gene; (3) presence of chromosomal TnMB1860; and (4) frequent presence of an F(31,36):A(4,20):B1 plasmid. Sequences belonging to this cluster (which we named "C2b sublineage") are highly enriched in septicemia samples and their associated genetic markers align with recent reports of an emerging, virulent sublineage of the C2 subclade, suggesting significant pathogenic potential. This is the first report of a ST131 strain belonging to subclade C2 contaminating green leafy vegetables. The detection of this uropathogenic clone in fresh food is alarming. This work suggests that ST131 continues to evolve, gaining selective advantages and new routes of transmission. This highlights the pressing need for rigorous epidemiological surveillance of ExPEC in vegetables with One Health perspective.

Detection of Horizontal Gene Transfer Mediated by Natural Conjugative Plasmids in E. coli.

Conjugation represents one of the main mechanisms facilitating horizontal gene transfer in Gram-negative bacteria. This work describes methods for the study of the mobilization of naturally occurring conjugative plasmids, using two naturally-occurring plasmids as an example. These protocols rely on the differential presence of selectable markers in donor, recipient, and conjugative plasmid. Specifically, the methods described include 1) the identification of natural conjugative plasmids, 2) the quantification of conjugation rates in solid culture, and 3) the diagnostic detection of the antibiotic resistance genes and plasmid replicon types in transconjugant recipients by polymerase chain reaction (PCR). The protocols described here have been developed in the context of studying the evolutionary ecology of horizontal gene transfer, to screen for the presence of conjugative plasmids carrying antibiotic-resistance genes in bacteria found in the environment. The efficient transfer of conjugative plasmids observed in these experiments in culture highlights the biological relevance of conjugation as a mechanism promoting horizontal gene transfer in general and the spread of antibiotic resistance in particular.

Experimental and In Silico Analysis of TEM ß-Lactamase Adaptive Evolution.

Multiple mutations often have non-additive (epistatic) phenotypic effects. Epistasis is of fundamental biological relevance but is not well understood mechanistically. Adaptive evolution, i.e., the evolution of new biochemical activities, is rich in epistatic interactions. To better understand the principles underlying epistasis during genetic adaptation, we studied the evolution of TEM-1 ß-lactamase variants exhibiting cefotaxime resistance. We report the collection of a library of 487 observed evolutionary trajectories for TEM-1 and determine the epistasis status based on cefotaxime resistance phenotype for 206 combinations of 2-3 TEM-1 mutations involving 17 positions under adaptive selective pressure. Gain-of-function (GOF) mutations are gatekeepers for adaptation. To see if GOF phenotypes can be inferred based solely on sequence data, we calculated the enrichment of GOF mutations in the different categories of epistatic pairs. Our results suggest that this is possible because GOF mutations are particularly enriched in sign and reciprocal sign epistasis, which leave a major imprint on the sequence space accessible to evolution. We also used FoldX to explore the relationship between thermodynamic stability and epistasis. We found that mutations in observed evolutionary trajectories tend to destabilize the folded structure of the protein, albeit their cumulative effects are consistently below the protein's free energy of folding. The destabilizing effect is stronger for epistatic pairs, suggesting that modest or local alterations in folding stability can modulate catalysis. Finally, we report a significant relationship between epistasis and the degree to which two protein positions are structurally and dynamically coupled, even in the absence of ligand.

Genomic analysis of plasmid content in food isolates of E. coli strongly supports its role as a reservoir for the horizontal transfer of virulence and antibiotic resistance genes.

The link between E. coli strains contaminating foods and human disease is unclear, with some reports supporting a direct transmission of pathogenic strains via food and others highlighting their role as reservoirs for resistance and virulence genes. Here we take a genomics approach, analyzing a large set of fully-assembled genomic sequences from E. coli available in GenBank. Most of the strains isolated in food are more closely related to each other than to clinical strains, arguing against a frequent direct transmission of pathogenic strains from food to the clinic. We also provide strong evidence of genetic exchanges between food and clinical strains that are facilitated by plasmids. This is based on an overlapped representation of virulence and resistance genes in plasmids isolated from these two sources. We identify clusters of phylogenetically-related plasmids that are largely responsible for the observed overlap and see evidence of specialization, with some food plasmid clusters preferentially transferring virulence factors over resistance genes. Consistent with these observations, food plasmids have a high mobilization potential based on their plasmid taxonomic unit classification and on an analysis of mobilization gene content. We report antibiotic resistance genes of high clinical relevance and their specific incompatibility group associations. Finally, we also report a striking enrichment for adhesins in food plasmids and their association with specific IncF replicon subtypes. The identification of food plasmids with specific markers (Inc and PTU combinations) as mediators of horizontal transfer between food and clinical strains opens new research avenues and should assist with the design of surveillance strategies.

Mechanisms of Theta Plasmid Replication in Enterobacteria and Implications for Adaptation to Its Host.

Plasmids are autonomously replicating sequences that help cells adapt to diverse stresses. Theta plasmids are the most frequent plasmid class in enterobacteria. They co-opt two host replication mechanisms: replication at oriC, a DnaA-dependent pathway leading to replisome assembly (theta class A), and replication fork restart, a PriA-dependent pathway leading to primosome assembly through primer extension and D-loop formation (theta classes B, C, and D). To ensure autonomy from the host's replication and to facilitate copy number regulation, theta plasmids have unique mechanisms of replication initiation at the plasmid origin of replication (ori). Tight plasmid copy number regulation is essential because of the major and direct impact plasmid gene dosage has on gene expression. The timing of plasmid replication and segregation are also critical for optimizing plasmid gene expression. Therefore, we propose that plasmid replication needs to be understood in its biological context, where complex origins of replication (redundant origins, mosaic and cointegrated replicons), plasmid segregation, and toxin-antitoxin systems are often present. Highlighting their tight functional integration with ori function, we show that both partition and toxin-antitoxin systems tend to be encoded in close physical proximity to the ori in a large collection of Escherichia coli plasmids. We also propose that adaptation of plasmids to their host optimizes their contribution to the host's fitness while restricting access to broad genetic diversity, and we argue that this trade-off between adaptation to host and access to genetic diversity is likely a determinant factor shaping the distribution of replicons in populations of enterobacteria.

Growth rate assays reveal fitness consequences of ß-lactamases.

NON-TECHNICAL SUMMARY: Antibiotic resistance is a global human health problem. We partnered with Dignity Health Mercy Medical Center to study antibiotic resistance in clinical isolates. We tested whether growth rates, a sensitive assay used to measure the fitness of bacterial samples, correlate with a clinical test to measure antibiotic resistance. We found a strong correlation between these two methods suggesting that growth rates could be reliably applied to evolutionary studies of clinically relevant problems. Moreover, the sensitivity of the growth rates assay enabled us to identify fitness effects of specific antibiotic resistance genes.

Genetic control of ColE1 plasmid stability that is independent of plasmid copy number regulation.

ColE1-like plasmid vectors are widely used for expression of recombinant genes in E. coli. For these vectors, segregation of individual plasmids into daughter cells during cell division appears to be random, making them susceptible to loss over time when no mechanisms ensuring their maintenance are present. Here we use the plasmid pGFPuv in a recA relA strain as a sensitized model to study factors affecting plasmid stability in the context of recombinant gene expression. We find that in this model, plasmid stability can be restored by two types of genetic modifications to the plasmid origin of replication (ori) sequence: point mutations and a novel 269 nt duplication at the 5' end of the plasmid ori, which we named DAS (duplicated anti-sense) ori. Combinations of these modifications produce a range of copy numbers and of levels of recombinant expression. In direct contradiction with the classic random distribution model, we find no correlation between increased plasmid copy number and increased plasmid stability. Increased stability cannot be explained by reduced levels of recombinant gene expression either. Our observations would be more compatible with a hybrid clustered and free-distribution model, which has been recently proposed based on detection of individual plasmids in vivo using super-resolution fluorescence microscopy. This work suggests a role for the plasmid ori in the control of segregation of ColE1 plasmids that is distinct from replication initiation, opening the door for the genetic regulation of plasmid stability as a strategy aimed at enhancing large-scale recombinant gene expression or bioremediation.

Network Analysis of Protein Adaptation: Modeling the Functional Impact of Multiple Mutations.

The evolution of new biochemical activities frequently involves complex dependencies between mutations and rapid evolutionary radiation. Mutation co-occurrence and covariation have previously been used to identify compensating mutations that are the result of physical contacts and preserve protein function and fold. Here, we model pairwise functional dependencies and higher order interactions that enable evolution of new protein functions. We use a network model to find complex dependencies between mutations resulting from evolutionary trade-offs and pleiotropic effects. We present a method to construct these networks and to identify functionally interacting mutations in both extant and reconstructed ancestral sequences (Network Analysis of Protein Adaptation). The time ordering of mutations can be incorporated into the networks through phylogenetic reconstruction. We apply NAPA to three distantly homologous ß-lactamase protein clusters (TEM, CTX-M-3, and OXA-51), each of which has experienced recent evolutionary radiation under substantially different selective pressures. By analyzing the network properties of each protein cluster, we identify key adaptive mutations, positive pairwise interactions, different adaptive solutions to the same selective pressure, and complex evolutionary trajectories likely to increase protein fitness. We also present evidence that incorporating information from phylogenetic reconstruction and ancestral sequence inference can reduce the number of spurious links in the network, whereas preserving overall network community structure. The analysis does not require structural or biochemical data. In contrast to function-preserving mutation dependencies, which are frequently from structural contacts, gain-of-function mutation dependencies are most commonly between residues distal in protein structure.

Fluorescence-Based Reporters for Detection of Mutagenesis in E. coli.

Mutagenesis in model organisms following exposure to chemicals is used as an indicator of genotoxicity. Mutagenesis assays are also used to study mechanisms of DNA homeostasis. This chapter focuses on detection of mutagenesis in prokaryotes, which boils down to two approaches: reporter inactivation (forward mutation assay) and reversion of an inactivating mutation (reversion mutation assay). Both methods are labor intensive, involving visual screening, quantification of colonies on solid media, or determining a Poisson distribution in liquid culture. Here, we present two reversion reporters for in vivo mutagenesis that produce a quantitative output, and thus have the potential to greatly reduce the amount of test chemical and labor involved in these assays. This output is obtained by coupling a TEM ß lactamase-based reversion assay with GFP fluorescence, either by placing the two genes on the same plasmid or by fusing them translationally and interrupting the N-terminus of the chimeric ORF with a stop codon. We also describe a reporter aimed at facilitating the monitoring of continuous mutagenesis in mutator strains. This reporter couples two reversion markers, allowing the temporal separation of mutation events in time, thus providing information about the dynamics of mutagenesis in mutator strains. Here, we describe these reporter systems, provide protocols for use, and demonstrate their key functional features using error-prone Pol I mutagenesis as a source of mutations.

Mechanisms of Theta Plasmid Replication.

Plasmids are autonomously replicating pieces of DNA. This article discusses theta plasmid replication, which is a class of circular plasmid replication that includes ColE1-like origins of replication popular with expression vectors. All modalities of theta plasmid replication initiate synthesis with the leading strand at a predetermined site and complete replication through recruitment of the host's replisome, which extends the leading strand continuously while synthesizing the lagging strand discontinuously. There are clear differences between different modalities of theta plasmid replication in mechanisms of DNA duplex melting and in priming of leading- and lagging-strand synthesis. In some replicons duplex melting depends on transcription, while other replicons rely on plasmid-encoded trans-acting proteins (Reps); primers for leading-strand synthesis can be generated through processing of a transcript or in other replicons by the action of host- or plasmid-encoded primases. None of these processes require DNA breaks. The frequency of replication initiation is tightly regulated to facilitate establishment in permissive hosts and to achieve a steady state. The last section of the article reviews how plasmid copy number is sensed and how this feedback modulates the frequency of replication.

Mechanisms of Theta Plasmid Replication.

Plasmids are autonomously replicating pieces of DNA. This chapter discusses theta plasmid replication, which is class of circular plasmid replication that includes ColE1-like origins of replication popular with expression vectors. All modalities of theta plasmid replication initiate synthesis with the leading-strand at a pre-determined site and complete replication through recruitment of the host's replisome, which extends the leading-strand continuously while synthesizing the lagging-strand discontinuously. There are clear differences between different modalities of theta plasmid replication in mechanisms of DNA duplex melting and in priming of leading- and lagging-strand synthesis. In some replicons duplex melting depends on transcription, while other replicons rely on plasmid-encoded trans-acting proteins (Reps); primers for leading-strand synthesis can be generated through processing of a transcript or in other replicons by the action of host- or plasmid-encoded primases. None of these processes require DNA breaks. The frequency of replication initiation is tightly regulated to facilitate establishment in permissive hosts and to achieve a steady state. The last section of the chapter reviews how plasmid copy number is sensed and how this feedback modulates the frequency of replication.

Mechanisms of plasmid segregation: have multicopy plasmids been overlooked?

Plasmids are self-replicating pieces of DNA typically bearing non-essential genes. Given that plasmids represent a metabolic burden to the host, mechanisms ensuring plasmid transmission to daughter cells are critical for their stable maintenance in the population. Here we review these mechanisms, focusing on two active partition strategies common to low-copy plasmids: par systems type I and type II. Both involve three components: an adaptor protein, a motor protein, and a centromere, which is a sequence area in the plasmid that is recognized by the adaptor protein. The centromere-bound adaptor nucleates polymerization of the motor, leading to filament formation, which can pull plasmids apart (par I) or push them towards opposite poles of the cell (par II). No such active partition mechanisms are known to occur in high copy number plasmids. In this case, vertical transmission is generally considered stochastic, due to the random distribution of plasmids in the cytoplasm. We discuss conceptual and experimental lines of evidence questioning the random distribution model and posit the existence of a mechanism for segregation in high copy number plasmids that moves plasmids to cell poles to facilitate transmission to daughter cells. This mechanism would involve chromosomally-encoded proteins and the plasmid origin of replication. Modulation of this proposed mechanism of segregation could provide new ways to enhance plasmid stability in the context of recombinant gene expression, which is limiting for large-scale protein production and for bioremediation.

Random mutagenesis by error-prone pol plasmid replication in Escherichia coli.

Directed evolution is an approach that mimics natural evolution in the laboratory with the goal of modifying existing enzymatic activities or of generating new ones. The identification of mutants with desired properties involves the generation of genetic diversity coupled with a functional selection or screen. Genetic diversity can be generated using PCR or using in vivo methods such as chemical mutagenesis or error-prone replication of the desired sequence in a mutator strain. In vivo mutagenesis methods facilitate iterative selection because they do not require cloning, but generally produce a low mutation density with mutations not restricted to specific genes or areas within a gene. For this reason, this approach is typically used to generate new biochemical properties when large numbers of mutants can be screened or selected. Here we describe protocols for an advanced in vivo mutagenesis method that is based on error-prone replication of a ColE1 plasmid bearing the gene of interest. Compared to other in vivo mutagenesis methods, this plasmid-targeted approach allows increased mutation loads and facilitates iterative selection approaches. We also describe the mutation spectrum for this mutagenesis methodology in detail, and, using cycle 3 GFP as a target for mutagenesis, we illustrate the phenotypic diversity that can be generated using our method. In sum, error-prone Pol I replication is a mutagenesis method that is ideally suited for the evolution of new biochemical activities when a functional selection is available.

Interplay between base excision repair activity and toxicity of 3-methyladenine DNA glycosylases in an E. coli complementation system.

DNA glycosylases carry out the first step of base excision repair by removing damaged bases from DNA. The N3-methyladenine (3MeA) DNA glycosylases specialize in alkylation repair and are either constitutively expressed or induced by exposure to alkylating agents. To study the functional and evolutionary significance of constitutive versus inducible expression, we expressed two closely related yeast 3MeA DNA glycosylases - inducible Saccharomyces cerevisiae MAG and constitutive S. pombe Mag1 - in a glycosylase-deficient Escherichia coli strain. In both cases, constitutive expression conferred resistance to alkylating agent exposure. However, in the absence of exogenous alkylation, high levels of expression of both glycosylases were deleterious. We attribute this toxicity to excessive glycosylase activity, since suppressing spMag1 expression correlated with improved growth in liquid culture, and spMag1 mutants exhibiting decreased glycosylase activity showed improved growth and viability. Selection of a random spMag1 mutant library for increased survival in the presence of exogenous alkylation resulted in the selection of hypomorphic mutants, providing evidence for the presence of a genetic barrier to the evolution of enhanced glycosylase activity when constitutively expressed. We also show that low levels of 3MeA glycosylase expression improve fitness in our glycosylase-deficient host, implying that 3MeA glycosylase activity is likely necessary for repair of endogenous lesions. These findings suggest that 3MeA glycosylase activity is evolutionarily conserved for repair of endogenously produced alkyl lesions, and that inducible expression represents a common strategy to rectify deleterious effects of excessive 3MeA activity in the absence of exogenous alkylation challenge.

The mutagenic footprint of low-fidelity Pol I ColE1 plasmid replication in E. coli reveals an extensive interplay between Pol I and Pol III.

ColE1 plasmid replication is unidirectional and requires two DNA polymerases: DNA polymerase I (Pol I) and DNA polymerase III (Pol III). Pol I initiates leading-strand synthesis by extending an RNA primer, allowing the Pol III holoenzyme to assemble and finish replication of both strands. The goal of the present work is to study the interplay between Pol I and Pol III during ColE1 plasmid replication, to gain new insights into Pol I function in vivo. Our approach consists of using mutations generated by a low-fidelity mutant of Pol I (LF-Pol I) during replication of a ColE1 plasmid as a footprint for Pol I replication. This approach allowed mapping areas of Pol I replication on the plasmid with high resolution. In addition, we were able to approximate the strandedness of Pol I mutations throughout the plasmid, allowing us to estimate the spectrum of the LF-Pol I in vivo. Our study produced the following three mechanistic insights: (1) we identified the likely location of the polymerase switch at ~200 bp downstream of replication initiation; (2) we found evidence suggesting that Pol I can replicate both strands, supporting earlier studies indicating a functional redundancy between Pol I and Pol III (3) we found evidence pointing to a specific role of Pol I during termination of lagging-strand replication. In addition, we illustrate how our strand-specific footprinting approach can be used to dissect factors modulating Pol I fidelity in vivo.

Quantifying plasmid copy number to investigate plasmid dosage effects associated with directed protein evolution.

Our laboratory specializes in directed protein evolution, i.e., evolution of proteins under defined selective pressures in the laboratory. Our target genes are encoded in ColE1 plasmids to facilitate the generation of libraries in vivo. We have observed that when random mutations are not restricted to the coding sequence of the target genes, directed evolution results in a strong positive selection of plasmid origin of replication (ori) mutations. Surprisingly, this is true even during evolution of new biochemical activities, when the activity that is being selected was not originally present. The selected plasmid ori mutations are diverse and produce a range of plasmid copy numbers, suggesting a complex interplay between ori and coding mutations rather than a simple enhancement of level of expression of the target gene. Thus, plasmid dosage may contribute significantly to evolution by fine-tuning levels of activity. Here, we present examples illustrating these observations as well as our methods for efficient quantification of plasmid copy number.

Unraveling a connection between DNA demethylation repair and cancer.

In this issue of Molecular Cell, Dango and Mosammaparast discover that the human oxidative demethylase ALKBH3 functions in complex with a DNA helicase to eliminate N3-methylcytosine lesions from ssDNA and that specific cancer cell lines are dependent on this activity for proliferation (Dango et al., 2011).

Network models of TEM ß-lactamase mutations coevolving under antibiotic selection show modular structure and anticipate evolutionary trajectories.

Understanding how novel functions evolve (genetic adaptation) is a critical goal of evolutionary biology. Among asexual organisms, genetic adaptation involves multiple mutations that frequently interact in a non-linear fashion (epistasis). Non-linear interactions pose a formidable challenge for the computational prediction of mutation effects. Here we use the recent evolution of ß-lactamase under antibiotic selection as a model for genetic adaptation. We build a network of coevolving residues (possible functional interactions), in which nodes are mutant residue positions and links represent two positions found mutated together in the same sequence. Most often these pairs occur in the setting of more complex mutants. Focusing on extended-spectrum resistant sequences, we use network-theoretical tools to identify triple mutant trajectories of likely special significance for adaptation. We extrapolate evolutionary paths (n = 3) that increase resistance and that are longer than the units used to build the network (n = 2). These paths consist of a limited number of residue positions and are enriched for known triple mutant combinations that increase cefotaxime resistance. We find that the pairs of residues used to build the network frequently decrease resistance compared to their corresponding singlets. This is a surprising result, given that their coevolution suggests a selective advantage. Thus, ß-lactamase adaptation is highly epistatic. Our method can identify triplets that increase resistance despite the underlying rugged fitness landscape and has the unique ability to make predictions by placing each mutant residue position in its functional context. Our approach requires only sequence information, sufficient genetic diversity, and discrete selective pressures. Thus, it can be used to analyze recent evolutionary events, where coevolution analysis methods that use phylogeny or statistical coupling are not possible. Improving our ability to assess evolutionary trajectories will help predict the evolution of clinically relevant genes and aid in protein design.

Roles of DNA polymerase I in leading and lagging-strand replication defined by a high-resolution mutation footprint of ColE1 plasmid replication.

DNA polymerase I (pol I) processes RNA primers during lagging-strand synthesis and fills small gaps during DNA repair reactions. However, it is unclear how pol I and pol III work together during replication and repair or how extensive pol I processing of Okazaki fragments is in vivo. Here, we address these questions by analyzing pol I mutations generated through error-prone replication of ColE1 plasmids. The data were obtained by direct sequencing, allowing an accurate determination of the mutation spectrum and distribution. Pol I's mutational footprint suggests: (i) during leading-strand replication pol I is gradually replaced by pol III over at least 1.3 kb; (ii) pol I processing of Okazaki fragments is limited to ∼20 nt and (iii) the size of Okazaki fragments is short (∼250 nt). While based on ColE1 plasmid replication, our findings are likely relevant to other pol I replicative processes such as chromosomal replication and DNA repair, which differ from ColE1 replication mostly at the recruitment steps. This mutation footprinting approach should help establish the role of other prokaryotic or eukaryotic polymerases in vivo, and provides a tool to investigate how sequence topology, DNA damage, or interactions with protein partners may affect the function of individual DNA polymerases.

Mutagenesis and functional selection protocols for directed evolution of proteins in E. coli.

The efficient generation of genetic diversity represents an invaluable molecular tool that can be used to label DNA synthesis, to create unique molecular signatures, or to evolve proteins in the laboratory. Here, we present a protocol that allows the generation of large (>10(11)) mutant libraries for a given target sequence. This method is based on replication of a ColE1 plasmid encoding the desired sequence by a low-fidelity variant of DNA polymerase I (LF-Pol I). The target plasmid is transformed into a mutator strain of E. coli and plated on solid media, yielding between 0.2 and 1 mutations/kb, depending on the location of the target gene. Higher mutation frequencies are achieved by iterating this process of mutagenesis. Compared to alternative methods of mutagenesis, our protocol stands out for its simplicity, as no cloning or PCR are involved. Thus, our method is ideal for mutational labeling of plasmids or other Pol I templates or to explore large sections of sequence space for the evolution of activities not present in the original target. The tight spatial control that PCR or randomized oligonucleotide-based methods offer can also be achieved through subsequent cloning of specific sections of the library. Here we provide protocols showing how to create a random mutant library and how to establish drug-based selections in E. coli to identify mutants exhibiting new biochemical activities.