Evolutionary trajectories of beta-lactamase CTX-M-1 cluster enzymes: predicting antibiotic resistance

Type IV-A3 CRISPR-Cas systems drive inter-plasmid conflicts by acquiring spacers in trans

Plasmid-encoded type IV-A CRISPR-Cas systems lack an acquisition module, feature a DinG helicase instead of a nuclease, and form ribonucleoprotein complexes of unknown biological functions. Type IV-A3 systems are carried by conjugative plasmids that often harbor antibiotic-resistance genes and their CRISPR array contents suggest a role in mediating inter-plasmid conflicts, but this function remains unexplored. Here, we demonstrate that a plasmid-encoded type IV-A3 system co-opts the type I-E adaptation machinery from its host, Klebsiella pneumoniae (K. pneumoniae), to update its CRISPR array. Furthermore, we reveal that robust interference of conjugative plasmids and phages is elicited through CRISPR RNA-dependent transcriptional repression. By silencing plasmid core functions, type IV-A3 impacts the horizontal transfer and stability of targeted plasmids, supporting its role in plasmid competition. Our findings shed light on the mechanisms and ecological function of type IV-A3 systems and demonstrate their practical efficacy for countering antibiotic resistance in clinically relevant strains.

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.

Slow Protein Dynamics Elicits New Enzymatic Functions by Means of Epistatic Interactions

Protein evolution depends on the adaptation of these molecules to different functional challenges. This occurs by tuning their biochemical, biophysical, and structural traits through the accumulation of mutations. While the role of protein dynamics in biochemistry is well recognized, there are limited examples providing experimental evidence of the optimization of protein dynamics during evolution. Here we report an NMR study of four variants of the CTX-M β-lactamases, in which the interplay of two mutations outside the active site enhances the activity against a cephalosporin substrate, ceftazidime. The crystal structures of these enzymes do not account for this activity enhancement. By using NMR, here we show that the combination of these two mutations increases the backbone dynamics in a slow timescale and the exposure to the solvent of an otherwise buried β-sheet. The two mutations located in this β-sheet trigger conformational changes in loops located at the opposite side of the active site. We postulate that the most active variant explores alternative conformations that enable binding of the more challenging substrate ceftazidime. The impact of the mutations in the dynamics is context-dependent, in line with the epistatic effect observed in the catalytic activity of the different variants. These results reveal the existence of a dynamic network in CTX-M β-lactamases that has been exploited in evolution to provide a net gain-of-function, highlighting the role of alternative conformations in protein evolution.

Prediction of Antibiotic Resistance Evolution by Growth Measurement of All Proximal Mutants of Beta-Lactamase

The antibiotic resistance crisis continues to threaten human health. Better predictions of the evolution of antibiotic resistance genes could contribute to the design of more sustainable treatment strategies. However, comprehensive prediction of antibiotic resistance gene evolution via laboratory approaches remains challenging. By combining site-specific integration and high-throughput sequencing, we quantified relative growth under the respective selection of cefotaxime or ceftazidime selection in ∼23,000 Escherichia coli MG1655 strains that each carried a unique, single-copy variant of the extended-spectrum β-lactamase gene blaCTX-M-14 at the chromosomal att HK022 site. Significant synergistic pleiotropy was observed within four subgenic regions, suggesting key regions for the evolution of resistance to both antibiotics. Moreover, we propose PEARP and PEARR, two deep-learning models with strong clinical correlations, for the prospective and retrospective prediction of blaCTX-M-14 evolution, respectively. Single to quintuple mutations of blaCTX-M-14 predicted to confer resistance by PEARP were significantly enriched among the clinical isolates harboring blaCTX-M-14 variants, and the PEARR scores matched the minimal inhibitory concentrations obtained for the 31 intermediates in all hypothetical trajectories. Altogether, we conclude that the measurement of local fitness landscape enables prediction of the evolutionary trajectories of antibiotic resistance genes, which could be useful for a broad range of clinical applications, from resistance prediction to designing novel treatment strategies.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Allogenous Selection of Mutational Collateral Resistance: Old Drugs Select for New Resistance Within Antibiotic Families

Allogeneous selection occurs when an antibiotic selects for resistance to more advanced members of the same family. The mechanisms of allogenous selection are (a) collateral expansion, when the antibiotic expands the gene and gene-containing bacterial populations favoring the emergence of other mutations, inactivating the more advanced antibiotics; (b) collateral selection, when the old antibiotic selects its own resistance but also resistance to more modern drugs; (c) collateral hyper-resistance, when resistance to the old antibiotic selects in higher degree for populations resistant to other antibiotics of the family than to itself; and (d) collateral evolution, when the simultaneous or sequential use of antibiotics of the same family selects for new mutational combinations with novel phenotypes in this family, generally with higher activity (higher inactivation of the antibiotic substrates) or broader spectrum (more antibiotics of the family are inactivated). Note that in some cases, collateral selection derives from collateral evolution. In this article, examples of allogenous selection are provided for the major families of antibiotics. Improvements in minimal inhibitory concentrations with the newest drugs do not necessarily exclude “old” antibiotics of the same family of retaining some selective power for resistance to the newest agents. If this were true, the use of older members of the same drug family would facilitate the emergence of mutational resistance to the younger drugs of the family, which is frequently based on previously established resistance traits. The extensive use of old drugs (particularly in low-income countries and in farming) might be significant for the emergence and selection of resistance to the novel members of the family, becoming a growing source of variation and selection of resistance to the whole family. In terms of future research, it could be advisable to focus antimicrobial drug discovery more on the identification of new targets and new (unique) classes of antimicrobial agents, than on the perpetual chemical exploitation of classic existing ones.

Ancient Antibiotics, Ancient Resistance

As the spread of antibiotic resistance threatens our ability to treat infections, avoiding the return of a preantibiotic era requires the discovery of new drugs. While therapeutic use of antibiotics followed by the inevitable selection of resistance is a modern phenomenon, these molecules and the genetic determinants of resistance were in use by environmental microbes long before humans discovered them. In this review, we discuss evidence that antibiotics and resistance were present in the environment before anthropogenic use, describing techniques including direct sampling of ancient DNA and phylogenetic analyses that are used to reconstruct the past. We also pay special attention to the ecological and evolutionary forces that have shaped the natural history of antibiotic biosynthesis, including a discussion of competitive versus signaling roles for antibiotics, proto-resistance, and substrate promiscuity of biosynthetic and resistance enzymes. Finally, by applying an evolutionary lens, we describe concepts governing the origins and evolution of biosynthetic gene clusters and cluster-associated resistance determinants. These insights into microbes' use of antibiotics in nature, a game they have been playing for millennia, can provide inspiration for discovery technologies and management strategies to combat the growing resistance crisis.

Primary and promiscuous functions coexist during evolutionary innovation through whole protein domain acquisitions

Molecular examples of evolutionary innovation are scarce and generally involve point mutations. Innovation can occur through larger rearrangements, but here experimental data is extremely limited. Integron integrases innovated from double-strand- toward single-strand-DNA recombination through the acquisition of the I2 α-helix. To investigate how this transition was possible, we have evolved integrase IntI1 to what should correspond to an early innovation state by selecting for its ancestral activity. Using synonymous alleles to enlarge sequence space exploration, we have retrieved 13 mutations affecting both I2 and the multimerization domains of IntI1. We circumvented epistasis constraints among them using a combinatorial library that revealed their individual and collective fitness effects. We obtained up to 10⁴-fold increases in ancestral activity with various asymmetrical trade-offs in single-strand-DNA recombination. We show that high levels of primary and promiscuous functions could have initially coexisted following I2 acquisition, paving the way for a gradual evolution toward innovation.

Occurrence of Antibiotic-Resistant Bacteria and Genes in Two Drinking Water Treatment and Distribution Systems in the North-West Province of South Africa

With the increasing spread of antimicrobial resistance, there is growing attention to the contribution made by drinking water systems. The potential health impact of two drinking water treatment and distribution systems (A and B) in the North-West Province of South Africa was determined by investigating the water quality and occurrence of antimicrobial-resistant heterotrophic bacteria and genes in the raw and treated water over four seasons. Most of the physicochemical parameters except for electrical conductivity were within permissible limits. Coliform bacteria reduced from raw to potable water except for counts higher than the threshold recorded in Summer and Winter. A total of 203 heterotrophic bacterial isolates were recovered on chromogenic R2A medium and subjected to susceptibility testing to twelve antibiotics. Most of the isolates were resistant to β-lactam antibiotics and Trimethoprim, whereas they were susceptible to Ciprofloxacin, Erythromycin, and Neomycin. The proportions of Cephalothin and Kanamycin-resistant isolates were significantly higher (p < 0.05) after treatment for site A, compared to significantly lower β-lactam, Oxytetracycline, and Trimethoprim-resistant isolates for B. Over 50% of isolates were of high risk, indicating their origin from high antibiotic-use sources. Seventy-one (35%) isolates were multidrug-resistant, out of which the majority (53.5%, n = 38) possessed the strA gene, followed by strB 21 (29.6%), dfrB 13 (18.3%), aadA 11 (15.5%), bla_CTX-M 5 (7.0%), and tetA 3 (4.2%). The 16S rRNA gene sequences of the isolates revealed strains belonging to eight bacterial families, some of which are clinically important.

Rapid and robust evolution of collateral sensitivity in Pseudomonas aeruginosa antibiotic-resistant mutants

The analysis of trade-offs, as collateral sensitivity, associated with the acquisition of antibiotic resistance, is mainly based on the use of model strains. However, the possibility of exploiting these trade-offs for fighting already resistant isolates has not been addressed in depth, despite the fact that bacterial pathogens are frequently antibiotic-resistant, forming either homogeneous or heterogeneous populations. Using a set of Pseudomonas aeruginosa-resistant mutants, we found that ceftazidime selects pyomelanogenic tobramycin-hypersusceptible mutants presenting chromosomal deletions in the analyzed genetic backgrounds. Since pyomelanogenic resistant mutants frequently coexist with other morphotypes in patients with cystic fibrosis, we analyzed the exploitation of this trade-off to drive extinction of heterogeneous resistant populations by using tobramycin/ceftazidime alternation. Our work shows that this approach is feasible because phenotypic trade-offs associated with the use of ceftazidime are robust. The identification of conserved collateral sensitivity networks may guide the rational design of evolution-based antibiotic therapies in P. aeruginosa infections.

Improved bacterial recombineering by parallelized protein discovery

Exploiting bacteriophage-derived homologous recombination processes has enabled precise, multiplex editing of microbial genomes and the construction of billions of customized genetic variants in a single day. The techniques that enable this, multiplex automated genome engineering (MAGE) and directed evolution with random genomic mutations (DIvERGE), are however, currently limited to a handful of microorganisms for which single-stranded DNA-annealing proteins (SSAPs) that promote efficient recombineering have been identified. Thus, to enable genome-scale engineering in new hosts, efficient SSAPs must first be found. Here we introduce a high-throughput method for SSAP discovery that we call "serial enrichment for efficient recombineering" (SEER). By performing SEER in Escherichia coli to screen hundreds of putative SSAPs, we identify highly active variants PapRecT and CspRecT. CspRecT increases the efficiency of single-locus editing to as high as 50% and improves multiplex editing by 5- to 10-fold in E. coli, while PapRecT enables efficient recombineering in Pseudomonas aeruginosa, a concerning human pathogen. CspRecT and PapRecT are also active in other, clinically and biotechnologically relevant enterobacteria. We envision that the deployment of SEER in new species will pave the way toward pooled interrogation of genotype-to-phenotype relationships in previously intractable bacteria.

[General principles of antibiotic resistance evolution in bacteria (review of literature).]

Currently, the impact of antibiotic resistance on human health is a worldwide problem and its study is of great interest from a molecular genetic, environmental and clinical view-point. This review summarizes the latest data about antibiotic resistance, the classification of microorganisms as sensitive and resistant to the action of antibiotics, reveals the concept of minimum inhibitory concentration from modern positions. The resistance of microorganisms to antibacterial agents can be intrinsic and acquired, as well as being one of the examples of evolution that are currently available for study. Modern methods of whole-genome sequencing and complex databases of nucleotide-tagged libraries give an idea of the multifaceted nature of the mechanisms of intrinsic resistance to antibiotics and are able to provide information on genes encoding metabolic enzymes and proteins that regulate the basic processes of the physiology of bacteria. The article describes the main ways of spreading the resistance of microorganisms, reflects the concepts of "founder effect" and the fitness cost of bacteria, which underlie the emergence and evolution of antibiotic resistance. It is shown that the origin of antibiotic resistance genes that human pathogens currently possess can be traced by studying the surrounding not only clinical, but also non-clinical (ecological) habitats. As well as microorganisms of the surrounding ecosystems are the donors of resistance genes in horizontal gene transfer.

Antibiotic Resistance Evolution Is Contingent on the Quorum-Sensing Response in Pseudomonas aeruginosa

Different works have explored independently the evolution toward antibiotic resistance and the role of eco-adaptive mutations in the adaptation to a new habitat (as the infected host) of bacterial pathogens. However, knowledge about the connection between both processes is still limited. We address this issue by comparing the evolutionary trajectories toward antibiotic resistance of a Pseudomonas aeruginosa lasR defective mutant and its parental wild-type strain, when growing in presence of two ribosome-targeting antibiotics. Quorum-sensing lasR defective mutants are selected in P. aeruginosa populations causing chronic infections. Further, we observed they are also selected in vitro as a first adaptation for growing in culture medium. By using experimental evolution and whole-genome sequencing, we found that the evolutionary trajectories of P. aeruginosa in presence of these antibiotics are different in lasR defective and in wild-type backgrounds, both at the phenotypic and the genotypic levels. Recreation of a set of mutants in both genomic backgrounds (either wild type or lasR defective) allowed us to determine the existence of negative epistatic interactions between lasR and antibiotic resistance determinants. These epistatic interactions could lead to mutual contingency in the evolution of antibiotic resistance when P. aeruginosa colonizes a new habitat in presence of antibiotics. If lasR mutants are selected first, this would constraint antibiotic resistance evolution. Conversely, when resistance mutations (at least those studied in the present work) are selected, lasR mutants may not be selected in presence of antibiotics. These results underlie the importance of contingency and epistatic interactions in modulating antibiotic resistance evolution.

CTX-M-33, a CTX-M-15 derivative conferring reduced susceptibility to carbapenems

CTX-M-type extended-spectrum ß-lactamases (ESBL) are widespread among Enterobacterales worldwide. The most common variant is CTX-M-15 hydrolyzing ceftazidime at high rate, but sparing carbapenems. We identified here CTX-M-33, a point mutant derivative of CTX-M-15 (Asp to Ser substitution at Ambler position 109), exhibiting a low carbapenemase activity. ß-Lactamase CTX-M-33 was identified in a Klebsiella pneumoniae isolate belonging to ST405, lacking the outer membrane protein OmpK36, that was resistant to broad-spectrum cephalosporins and ß-lactam/ß-lactamase inhibitor combinations, and displayed a decreased susceptibility to carbapenems. Comparative hydrolytic activity assays showed that CTX-M-33 hydrolyzed ceftazidime at a lower level than CTX-M-15, but significantly hydrolyzed meropenem. In addition, CTX-M-33 showed higher Mutant Prevention Concentration values and wider mutant selection window in presence of meropenem, in accordance with its observed hydrolytic properties. We identified here the very first CTX-M enzyme possessing a weak carbapenemase activity, that may correspond to an emerging phenomenon when considering its possibility to evolve from the widespread ESBL CTX-M-15.

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.

Collateral sensitivity constrains resistance evolution of the CTX-M-15 β-lactamase

Antibiotic resistance is a major challenge to global public health. Discovery of new antibiotics is slow and to ensure proper treatment of bacterial infections new strategies are needed. One way to curb the development of antibiotic resistance is to design drug combinations where the development of resistance against one drug leads to collateral sensitivity to the other drug. Here we study collateral sensitivity patterns of the globally distributed extended-spectrum β-lactamase CTX-M-15, and find three non-synonymous mutations with increased resistance against mecillinam or piperacillin-tazobactam that simultaneously confer full susceptibility to several cephalosporin drugs. We show in vitro and in mice that a combination of mecillinam and cefotaxime eliminates both wild-type and resistant CTX-M-15. Our results indicate that mecillinam and cefotaxime in combination constrain resistance evolution of CTX-M-15, and illustrate how drug combinations can be rationally designed to limit the resistance evolution of horizontally transferred genes by exploiting collateral sensitivity patterns.

The intrinsic resistome of Klebsiella pneumoniae

Molecular epidemiology studies aiming at understanding the acquisition of antimicrobial resistance by clinical isolates of Klebsiella pneumoniae are regularly published; however, information on the genes that contribute to its characteristic phenotype of resistance to antibiotics (intrinsic resistome) is scarce. To fill this gap, a K. pneumoniae transposon mutant library was screened and 171 mutants presenting changes in their susceptibility to antibiotics were selected, in which the transposon insertion site was determined in 75. Twenty-seven mutants for which insertion points had been previously identified were included in the analysis. A total of 102 mutants were selected for further studies. In 70 mutants the transposon was inserted in a gene with a known function, whilst in 19 the insertion occurred in genes encoding proteins with unknown functions and 13 insertions occurred in intergenic regions. Moreover, 87 of the insertions were localised in the chromosome, with 15 insertions located in the two plasmids carried by this strain. Whereas some of the mutated genes are already known to be involved in antimicrobial resistance (ampG, acrB, tolC), several of them are involved in regular processes of bacterial physiology, including K. pneumoniae virulence. Together with results published for other organisms, these results support that determinants involved in basic processes of bacterial physiology may contribute to antimicrobial resistance. These findings also indicate that, besides acquired resistance genes, plasmids may harbour other genes belonging to their backbone that can also be involved in resistance.

Inferring the shape of global epistasis

Genotype-phenotype relationships are notoriously complicated. Idiosyncratic interactions between specific combinations of mutations occur and are difficult to predict. Yet it is increasingly clear that many interactions can be understood in terms of global epistasis. That is, mutations may act additively on some underlying, unobserved trait, and this trait is then transformed via a nonlinear function to the observed phenotype as a result of subsequent biophysical and cellular processes. Here we infer the shape of such global epistasis in three proteins, based on published high-throughput mutagenesis data. To do so, we develop a maximum-likelihood inference procedure using a flexible family of monotonic nonlinear functions spanned by an I-spline basis. Our analysis uncovers dramatic nonlinearities in all three proteins; in some proteins a model with global epistasis accounts for virtually all of the measured variation, whereas in others we find substantial local epistasis as well. This method allows us to test hypotheses about the form of global epistasis and to distinguish variance components attributable to global epistasis, local epistasis, and measurement error.

Defining Substrate Specificity in the CTX-M Family: the Role of Asp240 in Ceftazidime Hydrolysis

The natural diversification of CTX-M β-lactamases led to the emergence of Asp240Gly variants in the clinic that confer reduced susceptibility to ceftazidime (CAZ). In this study, we compared the impact of this substitution on CAZ and ceftazidime-avibactam (CZA) MICs against isogenic Escherichia coli strains with different porin deficiencies. Our results show a noticeable increase in CAZ resistance in clones expressing Asp240Gly-harboring CTX-M when combined with OmpF porin deficiency. Kinetic analysis revealed that the k_cat/K_m for CAZ was 5- to 15-fold higher for all Asp240Gly variants but remained 200- to 725-fold lower than that for cefotaxime (CTX). In vitro selection of CAZ-resistant clones yielded nonsusceptible CTX-M producers (MIC of >16 μg/ml) only after overnight incubation; the addition of avibactam (AVI) decreased MICs to a susceptible range against these variants. In contrast, the use of CZA as a selective agent did not yield resistant clones. AVI inactivated both CTX-M-12 and CTX-M-96, with an apparent inhibition constant comparable to that of SHV-2 and 1,000-fold greater than that of PER-2 and CMY-2, and k₂/K for CTX-M-12 was 24- and 35-fold higher than that for CTX-M-96 and CTX-M-15, respectively. Molecular modeling suggests that AVI interacts similarly with CTX-M-96 and CTX-M-15. We conclude that the impact of Asp240Gly in resistance may arise when other mechanisms are also present (i.e., OmpF deficiency). Additionally, CAZ selection could favor the emergence of CAZ-resistant subpopulations. These results define the role of Asp240 and the impact of the -Gly substitution and allow us to hypothesize that the use of CZA is an effective preventive strategy to delay the development of resistance in this family of extended-spectrum β-lactamases.

Ceftazidime Is the Key Diversification and Selection Driver of VIM-Type Carbapenemases

In recent decades, carbapenems have been considered the last line of antibiotic therapy for Gram-negative bacterial infections. Unfortunately, strains carrying a high diversity of β-lactamases able to hydrolyze carbapenems have emerged in the clinical setting. Among them, VIM β-lactamases have diversified in a bloom of variants. The evolutionary reconstructions performed in this work revealed that, at the end of the 1980s, two independent events involving diversification from VIM-2 and VIM-4 produced at least 25 VIM variants. Later, a third event involving diversification from VIM-1 occurred in the mid-1990s. In a second approach to understanding the emergence of VIM carbapenemases, 44 mutants derived from VIM-2 and VIM-4 were obtained by site-directed mutagenesis based on those positions predicted to be under positive selection. These variants were expressed in an isogenic context. The more-evolved variants yielded increased levels of hydrolytic efficiency toward ceftazidime to a higher degree than toward carbapenems. In fact, an antagonist effect was frequently observed. These results led us to develop an experimental-evolution step. When Escherichia coli strains carrying VIM-2 or VIM-4 were submitted to serial passages at increasing concentrations of carbapenems or ceftazidime, more-efficient new variants (such as VIM-11 and VIM-1, with N165S [bearing a change from N to S at position 165] and R228S mutations, respectively) were only obtained when ceftazidime was present. Therefore, the observed effect of ceftazidime in the diversification and selection of VIM variants might help to explain the recent bloom of carbapenemase diversity, and it also represents another example of the potential universal effect exerted by ceftazidime in the selection of more-efficient β-lactamase variants, as in TEM, CTX-M, or KPC enzymes.

One of the objectives recently proposed by the World Health Organization (WHO) Assembly in the global plan on antimicrobial resistance was to improve the understanding and knowledge of antimicrobial resistance. In the present work, we paid attention to the drivers of diversification and selection of new carbapenemases in Gram-negative bacteria, which occupy one of the most critical places in the WHO priority list of antibiotic-resistant microorganisms. Based on evolutionary-reconstruction, site-directed-mutagenesis, and experimental-evolution approaches, we proposed a critical role of ceftazidime exposure in the selection of VIM carbapenemase variants. This surprising finding is also applicable to other β-lactamases, indicating that ceftazidime, and not other antibiotics, might have a universal effect in the diversification of β-lactamases. Our results might help to define future strategies to reconsider the extended use of ceftazidime.

Identification of Plasmid-Encoded sRNAs in a blaNDM-1-Harboring Multidrug-Resistance Plasmid pNDM-HK in Enterobacteriaceae

Small RNAs (sRNAs) play significant roles in regulating gene expression post-transcriptionally in response to environmental changes in bacteria. In this work, we identified and characterized six novel sRNAs from an emerging multidrug-resistance (MDR) plasmid pNDM-HK, a New Delhi metallo-β-lactamase 1 gene (bla_NDM−1)-carrying IncL/M plasmid that has caused worldwide threat in recent years. These sRNAs are located at different regions of pNDM-HK, such as replication, stability, and variable regions. Moreover, one of the plasmid-encoded sRNAs (NDM-sR3) functions in an Hfq-dependent manner and possibly plays roles in the fitness of pNDM-HK carrying bacteria. In addition, we attempted to construct the phylogenetic tree based on these novel sRNAs and surprisingly, the sRNA-phylogenetic tree provided significant information about the evolutionary pathway of pNDM-HK, including possible gene acquisition and insertion from relevant plasmids. Moreover, the sRNA-phylogenetic tree can specifically cluster the IncM2 type and distinguish it from other IncL/M subtypes. In summary, this is the first study to systematically identify and characterize sRNAs from clinically-isolated MDR plasmids. We believe that these newly found sRNAs could lead to further understanding and new directions to study the evolution and dissemination of the clinically MDR bacterial plasmids.

The number of genes encoding repeat domain-containing proteins positively correlates with genome size in amoebal giant viruses

Curiously, in viruses, the virion volume appears to be predominantly driven by genome length rather than the number of proteins it encodes or geometric constraints. With their large genome and giant particle size, amoebal viruses (AVs) are ideally suited to study the relationship between genome and virion size and explore the role of genome plasticity in their evolutionary success. Different genomic regions of AVs exhibit distinct genealogies. Although the vertically transferred core genes and their functions are universally conserved across the nucleocytoplasmic large DNA virus (NCLDV) families and are essential for their replication, the horizontally acquired genes are variable across families and are lineage-specific. When compared with other giant virus families, we observed a near–linear increase in the number of genes encoding repeat domain-containing proteins (RDCPs) with the increase in the genome size of AVs. From what is known about the functions of RDCPs in bacteria and eukaryotes and their prevalence in the AV genomes, we envisage important roles for RDCPs in the life cycle of AVs, their genome expansion, and plasticity. This observation also supports the evolution of AVs from a smaller viral ancestor by the acquisition of diverse gene families from the environment including RDCPs that might have helped in host adaption.

Functional Characterization of CTX-M-14 and CTX-M-15 β-Lactamases by In Vitro DNA Shuffling

This work investigated the molecular events driving the evolution of the CTX-M-type β-lactamases by the use of DNA shuffling of fragments of the bla_CTX-M-14 and bla_CTX-M-15 genes. Analysis of a total of 51 hybrid enzymes showed that enzymatic activity could be maintained in most cases, yet hybrids that were active possessed fewer amino acid substitutions than those that were inactive, suggesting that point mutations in the constructs rather than reshuffling of the fragments of the two target genes would more likely cause disruption of CTX-M activity. For example, the P⁶⁷L and L²⁶¹P changes in a CTX-M-14 fragment could completely abolish the activity of the enzyme on all antibiotics tested. Structural analysis showed that L²¹⁶ was located in the active-site β sheet and might interact with the adjacent hydrophobic residues to stabilize the active-site β sheet and maintain the integrity of the enzyme active site. Likewise, a single amino acid substitution, E⁶⁴K, was found to exhibit a significant suppressive effect on CTX-M-15 activity. Structural analysis showed that E⁶⁴ might form a salt bridge with R⁴⁴, disruption of which might affect CTX-M-15 activity. Further analysis of the structure-function relationship of a range of mutant enzymes confirmed that, as can be expected, unstable enzymes lose their activity and avoid selective events. These findings suggest that the distal pockets could also contribute to the activity of the enzymes and may be regarded as alternative targets for inhibitor development.

Genetic and biochemical characterisation of CTX-M-37 extended-spectrum β-lactamase from an Enterobacter cloacae clinical isolate from Mongolia

OBJECTIVES

The aims of this study were to determine the resistance level of a bla_CTX-M-37-carrying Enterobacter cloacae isolate from Mongolia, to analyse kinetic parameters of the purified enzyme and to compare the genetic environment of the gene.

METHODS

Minimum inhibitory concentrations (MICs) were determined using the Clinical and Laboratory Standards Institute (CLSI) agar dilution method. Purified CTX-M-37 enzyme was used to determined kinetic parameters. The genetic environment of the blaCTX-M-37 gene in E. cloacae was compared with a Kluyvera cryocrescens isolate.

RESULTS

The E. cloacae isolate showed relatively low-level resistance to cefotaxime (MIC=16mg/L) compared with a CTX-M-3-producing strain (MIC=256mg/L), and CTX-M-37 had a lower k_cat/K_m value for cefotaxime (2.0μM^-1s^-1) compared with CTX-M-3 (3.5μM^-1s^-1), possibly due to Asn114Asp substitution. The bla_CTX-M-37 gene in the E. cloacae isolate was carried on a conjugative plasmid and was associated with an ISEcp1 element containing the -35 and -10 putative promoter sequences TTGAAA and TACAAT, respectively, unlike in the K. cryocrescens isolate.

CONCLUSIONS

The CTX-M-37-producing E. cloacae isolate showed relatively low-level resistance to cefotaxime and the purified enzyme had lower kinetic parameters as the result of Asn114Asp substitution. Presence of an ISEcp1 element and putative promoters upstream of the bla_CTX-M-37 gene in E. cloacae, but not in the K. cryocrescens isolate, indicated their roles in mobilisation and expression of the gene.

Beyond the Hypercube: Evolutionary Accessibility of Fitness Landscapes with Realistic Mutational Networks

Evolutionary pathways describe trajectories of biological evolution in the space of different variants of organisms (genotypes). The probability of existence and the number of evolutionary pathways that lead from a given genotype to a better-adapted genotype are important measures of accessibility of local fitness optima and the reproducibility of evolution. Both quantities have been studied in simple mathematical models where genotypes are represented as binary sequences of two types of basic units, and the network of permitted mutations between the genotypes is a hypercube graph. However, it is unclear how these results translate to the biologically relevant case in which genotypes are represented by sequences of more than two units, for example four nucleotides (DNA) or 20 amino acids (proteins), and the mutational graph is not the hypercube. Here we investigate accessibility of the best-adapted genotype in the general case of K > 2 units. Using computer generated and experimental fitness landscapes we show that accessibility of the global fitness maximum increases with K and can be much higher than for binary sequences. The increase in accessibility comes from the increase in the number of indirect trajectories exploited by evolution for higher K. As one of the consequences, the fraction of genotypes that are accessible increases by three orders of magnitude when the number of units K increases from 2 to 16 for landscapes of size N ∼ 10⁶ genotypes. This suggests that evolution can follow many different trajectories on such landscapes and the reconstruction of evolutionary pathways from experimental data might be an extremely difficult task.

Biological evolution is driven by heritable, genetic alterations that affect the fitness of organisms. However, the pool of “fitter” variants (genotypes) is often restricted and it is not at all obvious how evolution finds its way from low-fitness to high-fitness genotypes in a complex, multidimensional “fitness landscapes” with many peaks (fit organisms) and valleys (unfit ones). To address this question we investigate how likely it is for biological evolution to find a way “uphill” from a lower-fitness organism to the best adapted organism. We discover that the accessibility of the fittest organism depends on the number of types of basic “units” used to encode genotypes. These units can be, for example, the four DNA nucleotides A,T,C,G, or the ∼20 amino acids used for synthesizing proteins, and the choice of the most appropriate unit is dictated by how the genotypes and the fitnesses are related—a relationship that researchers have begun to unveil only recently. We find that increasing the number of units strongly increases the probability that there will be at least one uphill path to the best-adapted genotype, and the number of evolutionary pathways leading to it. Our findings suggest that biological evolution can follow many more pathways than previously thought.

Functional Trade-Offs in Promiscuous Enzymes Cannot Be Explained by Intrinsic Mutational Robustness of the Native Activity

The extent to which an emerging new function trades off with the original function is a key characteristic of the dynamics of enzyme evolution. Various cases of laboratory evolution have unveiled a characteristic trend; a large increase in a new, promiscuous activity is often accompanied by only a mild reduction of the native, original activity. A model that associates weak trade-offs with “evolvability” was put forward, which proposed that enzymes possess mutational robustness in the native activity and plasticity in promiscuous activities. This would enable the acquisition of a new function without compromising the original one, reducing the benefit of early gene duplication and therefore the selection pressure thereon. Yet, to date, no experimental study has examined this hypothesis directly. Here, we investigate the causes of weak trade-offs by systematically characterizing adaptive mutations that occurred in two cases of evolutionary transitions in enzyme function: (1) from phosphotriesterase to arylesterase, and (2) from atrazine chlorohydrolase to melamine deaminase. Mutational analyses in various genetic backgrounds revealed that, in contrast to the prevailing model, the native activity is less robust to mutations than the promiscuous activity. For example, in phosphotriesterase, the deleterious effect of individual mutations on the native phosphotriesterase activity is much larger than their positive effect on the promiscuous arylesterase activity. Our observations suggest a revision of the established model: weak trade-offs are not caused by an intrinsic robustness of the native activity and plasticity of the promiscuous activity. We propose that upon strong adaptive pressure for the new activity without selection against the original one, selected mutations will lead to the largest possible increases in the new function, but whether and to what extent they decrease the old function is irrelevant, creating a bias towards initially weak trade-offs and the emergence of generalist enzymes.

Understanding how enzymes evolve is a fundamental question that can help us decipher not only the mechanisms of evolution on a higher level, i.e., whole organisms, but also advances our knowledge of sequence-structure-function relationships as a guide to artificial evolution in the test tube. An important yet unexplained phenomenon occurs during the evolution of a new enzymatic function; it has been observed that new and ancestral functions often trade-off only weakly, meaning the original native activity is initially maintained at a high level despite drastic improvement of the new promiscuous activity. It has previously been proposed that weak trade-offs occur because the native activity is robust to mutations while the promiscuous activity is not. However, the present work contradicts this hypothesis, based on the detailed characterization of mutational effects on both activities in two examples of enzyme evolution. We propose an alternative explanation: the weak activity trade-off is consistent with being a by-product of strong selection for the new activity rather than an intrinsic property of the native activity.

Co-diversification of Enterococcus faecium Core Genomes and PBP5: Evidences of pbp5 Horizontal Transfer

Ampicillin resistance has greatly contributed to the recent dramatic increase of a cluster of human adapted Enterococcus faecium lineages (ST17, ST18, and ST78) in hospital-based infections. Changes in the chromosomal pbp5 gene have been associated with different levels of ampicillin susceptibility, leading to protein variants (designated as PBP5 C-types to keep the nomenclature used in previous works) with diverse degrees of reduction in penicillin affinity. Our goal was to use a comparative genomics approach to evaluate the relationship between the diversity of PBP5 among E. faecium isolates of different phylogenomic groups as well as to assess the pbp5 transferability among isolates of disparate clonal lineages. The analyses of 78 selected E. faecium strains as well as published E. faecium genomes, suggested that the diversity of pbp5 mirrors the phylogenomic diversification of E. faecium. The presence of identical PBP5 C-types as well as similar pbp5 genetic environments in different E. faecium lineages and clones from quite different geographical and environmental origin was also documented and would indicate their horizontal gene transfer among E. faecium populations. This was supported by experimental assays showing transfer of large (≈180–280 kb) chromosomal genetic platforms containing pbp5 alleles, ponA (transglycosilase) and other metabolic and adaptive features, from E. faecium donor isolates to suitable E. faecium recipient strains. Mutation profile analysis of PBP5 from available genomes and strains from this study suggests that the spread of PBP5 C-types might have occurred even in the absence of a significant ampicillin resistance phenotype. In summary, genetic platforms containing pbp5 sequences were stably maintained in particular E. faecium lineages, but were also able to be transferred among E. faecium clones of different origins, emphasizing the growing risk of further spread of ampicillin resistance in this nosocomial pathogen.

Antimicrobial Susceptibility and Molecular Epidemiology of Multidrug-Resistant Klebsiella pneumoniae in Central China

Extended-spectrum-beta-lactamase-producing and carbapenemase-producing Klebsiella pneumoniae strains have rapidly spread through clinical units worldwide. This study investigated the epidemiology and resistance profiles of K. pneumoniae strains isolated in central China between 2009 and 2014. Antimicrobial susceptibility testing and polymerase chain reaction were used to investigate the prevalence of extended-spectrum beta-lactamases (ESBL) and carbapenemase production by these K. pneumoniae strains, and the prevalence of K. pneumoniae carbapenemase (KPC)-producing K. pneumoniae was investigated by multilocus sequence typing. Carbapenem resistance has emerged as a major concern in K. pneumoniae infections, as phenotype testing has detected carbapenemases in nearly 20% of isolates. KPC-producing isolates in a local epidemic were clonally related, with ST11 being the reservoir for the bla_KPC-2 gene and ESBL genes. During the 6-year collection period, the prevalence of ESBLs was dynamic, and suggested that bla_CTX-M-55 might become prevalent in the future. Our findings demonstrate the high prevalence of carbapenemase- and ESBL-producing K. pneumoniae in central China and predict a future local epidemic of KPC-2 and CTX-M-55.

Metabolic modelling in a dynamic evolutionary framework predicts adaptive diversification of bacteria in a long-term evolution experiment

BACKGROUND

Predicting adaptive trajectories is a major goal of evolutionary biology and useful for practical applications. Systems biology has enabled the development of genome-scale metabolic models. However, analysing these models via flux balance analysis (FBA) cannot predict many evolutionary outcomes including adaptive diversification, whereby an ancestral lineage diverges to fill multiple niches. Here we combine in silico evolution with FBA and apply this modelling framework, evoFBA, to a long-term evolution experiment with Escherichia coli.

RESULTS

Simulations predicted the adaptive diversification that occurred in one experimental population and generated hypotheses about the mechanisms that promoted coexistence of the diverged lineages. We experimentally tested and, on balance, verified these mechanisms, showing that diversification involved niche construction and character displacement through differential nutrient uptake and altered metabolic regulation.

CONCLUSION

The evoFBA framework represents a promising new way to model biochemical evolution, one that can generate testable predictions about evolutionary and ecosystem-level outcomes.

Predicting and preventing antimicrobial resistance utilizing pharmacodynamics: Part I gram positive bacteria

INTRODUCTION

Antimicrobial resistance is a potentially inevitable consequence of widespread use of antibiotics in the healthcare system. An enhanced understanding of pharmacodynamic (PD) targets that prevent antimicrobial resistance development will improve currently availably therapies and help to guide future drug development strategies. Current in vitro methods to predict bacterial resistance to antimicrobials consist of serial dilution experiments, determination of the mutant prevention concentration (MPC), mutant selection window (MSW), and human simulated pharmacodynamics studies. Clinical trial data and real -world surveillance studies can help validate or disprove in vitro modeling.

AREAS COVERED

This review will discuss methods of predicting development of resistance and how the use of pharmacodynamics can reduce or eliminate the emergence of resistance among Staphylococcus aureus, Streptococcus pneumoniae, and Enterococcus species.

EXPERT OPINION

Pharmacodynamic targets can be used successfully to guide antimicrobial therapy to prevent resistance development. Currently, PD targets do not take into consideration horizontal resistance gene transfer and various factors may lead to different PD targets based on sites of infection. Further research is necessary to guide future drug development strategies and optimize new drug therapies.

Structural and Kinetic Insights into the "Ceftazidimase" Behavior of the Extended-Spectrum β-Lactamase CTX-M-96

Diversification of the CTX-M β-lactamases led to the emergence of variants responsible for decreased susceptibility to ceftazidime, like the Asp240Gly-harboring "ceftazidimases". We solved the crystallographic structure of the Asp240Gly variant CTX-M-96 at 1.2 Å and evaluated the role of Asp240 in the activity toward oxyimino-cephalosporins through simulated models and kinetics. There seem to be subtle changes in the conformation of the active site cavity of CTX-M-96, compared to enzyme variants harboring the Asp240, and these small rearrangements could be due to localized shifts in the environment of the β3 strand. According to the crystallographic evidence, CTX-M-96 presents a "compact" active site, which in spite of its reduced cavity seems to allow the proper interaction with oxyimino-cephalosporins, as suggested by simulated models. The term "ceftazidimases" that is currently applied for the Asp240Gly-harboring CTX-M variants should be used carefully. Structural differences between CTX-M harboring the Asp240Gly mutation (and also probably others like those at Pro167) do not seem to be conclusive to determine the "ceftazidimase" behavior observed in vivo, which is in turn partially supported by the mild improvement in the catalytic efficiency toward ceftazidime by CTX-M-96 and similar enzymes, compared to "parental" Asp240-harboring variants. In addition, it is observed that alterations in OmpF expression could act synergistically with CTX-M-96 for yielding clinical resistance toward ceftazidime. We therefore propose that the observed resistance in vivo is due to the sum of synergic mechanisms, and the term "cefotaximases associated with ceftazidime resistance" could be conveniently used to describe CTX-M harboring the Asp240Gly substitution.

Delayed commitment to evolutionary fate in antibiotic resistance fitness landscapes

Predicting evolutionary paths to antibiotic resistance is key for understanding and controlling drug resistance. When considering a single final resistant genotype, epistatic contingencies among mutations restrict evolution to a small number of adaptive paths. Less attention has been given to multi-peak landscapes, and while specific peaks can be favoured, it is unknown whether and how early a commitment to final fate is made. Here we characterize a multi-peaked adaptive landscape for trimethoprim resistance by constructing all combinatorial alleles of seven resistance-conferring mutations in dihydrofolate reductase. We observe that epistatic interactions increase rather than decrease the accessibility of each peak; while they restrict the number of direct paths, they generate more indirect paths, where mutations are adaptively gained and later adaptively lost or changed. This enhanced accessibility allows evolution to proceed through many adaptive steps while delaying commitment to genotypic fate, hindering our ability to predict or control evolutionary outcomes.

blaCTX-M-1/9/1 Hybrid Genes May Have Been Generated from blaCTX-M-15 on an IncI2 Plasmid

Three hybrid CTX-M β-lactamases, CTX-M-64, CTX-M-123, and CTX-M-132, with N and C termini matching CTX-M-1 group enzymes and centers matching CTX-M-9 group enzymes, have been identified. The hybrid gene sequences suggested recombination between blaCTX-M-15 and blaCTX-M-14, the two most common blaCTX-M variants worldwide. However, blaCTX-M-64 and blaCTX-M-123 are found in an ISEcp1-blaCTX-M transposition unit with a 45-bp "spacer," rather than the 48 bp usually associated with blaCTX-M-15, and 112 bp of IncA/C plasmid backbone. This is closer to the context of blaCTX-M-55, which has one nucleotide difference from blaCTX-M-15, on IncI2 plasmid pHN1122-1. Here, we characterized an IncI2 plasmid carrying blaCTX-M-15 with a 45-bp spacer (pHNY2-1) by complete sequencing and also sequenced IncI2 plasmids carrying blaCTX-M-64 (pHNAH46-1) or blaCTX-M-132 (pHNLDH19) and an IncI1 plasmid carrying blaCTX-M-123 (pHNAH4-1). pHNY2-1 has the same ISEcp1-blaCTX-M-IncA/C insertion as pHN1122-1, pHNAH46-1, and pHNLDH19, and all four plasmid backbones are almost identical. pHNAH4-1 (IncI1 sequence type 108 [ST108]) carries a transposition unit that includes a 2,720-bp fragment of the IncI2 backbone, suggesting ISEcp1-mediated transfer of blaCTX-M-IncA/C-IncI2 to an IncI1 plasmid. All three hybrid blaCTX-M genes may have resulted from recombination between blaCTX-M-14 and blaCTX-M-15 with a 45-bp spacer on an IncI2 plasmid. Five additional Escherichia coli isolates of different sequence types from different provinces, farms, and/or animals had blaCTX-M-64 on a pHNAH46-1-like IncI2 plasmid and 9 had blaCTX-M-123 on a pHNAH4-1-like IncI1 ST108 plasmid. Thus, epidemic IncI plasmids may be responsible for the spread of blaCTX-M-64 and blaCTX-M-123 between different animals and different locations in China.

How structural and physicochemical determinants shape sequence constraints in a functional enzyme

The need for interfacing structural biology and biophysics to molecular evolution is being increasingly recognized. One part of the big problem is to understand how physics and chemistry shape the sequence space available to functional proteins, while satisfying the needs of biology. Here we present a quantitative, structure-based analysis of a high-resolution map describing the tolerance to all substitutions in all positions of a functional enzyme, namely a TEM lactamase previously studied through deep sequencing of mutants growing in competition experiments with selection against ampicillin. Substitutions are rarely observed within 7 Å of the active site, a stringency that is relaxed slowly and extends up to 15–20 Å, with buried residues being especially sensitive. Substitution patterns in over one third of the residues can be quantitatively modeled by monotonic dependencies on amino acid descriptors and predictions of changes in folding stability. Amino acid volume and steric hindrance shape constraints on the protein core; hydrophobicity and solubility shape constraints on hydrophobic clusters underneath the surface, and on salt bridges and polar networks at the protein surface together with charge and hydrogen bonding capacity. Amino acid solubility, flexibility and conformational descriptors also provide additional constraints at many locations. These findings provide fundamental insights into the chemistry underlying protein evolution and design, by quantitating links between sequence and different protein traits, illuminating subtle and unexpected sequence-trait relationships and pinpointing what traits are sacrificed upon gain-of-function mutation.

A method for the detection of antibiotic resistance markers in clinical strains of Escherichia coli using MALDI mass spectrometry

Matrix-assisted laser-desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF MS) is one of the most widely used mass spectrometry based approaches for bacterial identification and classification. The relatively simple sample preparation requirements and the speed of analysis which can usually be completed within a few minutes have resulted in the adoption and assimilation of MALDI-TOF MS into the routine diagnostic workflow of Clinical microbiology laboratories worldwide. This study describes the facilitation of bacterial discrimination based on antibiotic resistance markers through the implementation of MALDI-TOF MS. The periplasmic compartment of whole bacterial cells contains several proteins which confer antibiotic resistance in the Enterobacteriaceae. In order to reduce the complexity of the sample to be analysed via MALDI-TOF MS, the periplasm was extracted and subjected to in solution tryptic digestion followed by nano-LC separation. This method, established that peptide sequence biomarkers from several classes of antibiotic resistance proteins could be predicted using protein/peptide database tools such as Mascot. Biomarkers for a CTX-M-1 group extended spectrum β-lactamase, CMY-2 an Amp-C β-lactamase, VIM a metallo-β-lactamase, TEM a β-lactamase and KanR an aminoglycoside modifying enzyme were detected. This allowed for discrimination at a species level and at an almost identical strain level where the only difference between strains was the carriage of a modified antibiotic resistance carrying plasmid. This method also was able to detect some of these biomarkers in clinical strains where multiple resistance mechanisms were present.

In vitro prediction of the evolution of GES-1 β-lactamase hydrolytic activity

Resistance to β-lactams is constantly increasing due to the emergence of totally new enzymes but also to the evolution of preexisting β-lactamases. GES-1 is a clinically relevant extended-spectrum β-lactamase (ESBL) that hydrolyzes penicillins and broad-spectrum cephalosporins but spares monobactams and carbapenems. However, several GES-1 variants (i.e., GES-2 and GES-5) previously identified among clinical isolates display an extended spectrum of activity toward carbapenems. To study the evolution potential of the GES-1 β-lactamase, this enzyme was submitted to in vitro-directed evolution, with selection on increasing concentrations of the cephalosporin cefotaxime, the monobactam aztreonam, or the carbapenem imipenem. The highest resistance levels were conferred by a combination of up to four substitutions. The A6T-E104K-G243A variant selected on cefotaxime and the A6T-E104K-T237A-G243A variant selected on aztreonam conferred high resistance to cefotaxime, ceftazidime, and aztreonam. Conversely, the A6T-G170S variant selected on imipenem conferred high resistance to imipenem and cefoxitin. Of note, the A6T substitution involved in higher MICs for all β-lactams is located in the leader peptide of the GES enzyme and therefore is not present in the mature protein. Acquired cross-resistance was not observed, since selection with cefotaxime or aztreonam did not select for resistance to imipenem, and vice versa. Here, we demonstrate that the β-lactamase GES-1 exhibits peculiar properties, with a significant potential to gain activity against broad-spectrum cephalosporins, monobactams, and carbapenems.

General principles of antibiotic resistance in bacteria

Given the impact of antibiotic resistance on human health, its study is of great interest from a clinical view- point. In addition, antibiotic resistance is one of the few examples of evolution that can be studied in real time. Knowing the general principles involved in the acquisition of antibiotic resistance is therefore of interest to clinicians, evolutionary biologists and ecologists. The origin of antibiotic resistance genes now possessed by human pathogens can be traced back to environmental microorganisms. Consequently, a full understanding of the evolution of antibiotic resistance requires the study of natural environments as well as clinical ecosystems. Updated information on the evolutionary mechanisms behind resistance, indicates that ecological connectivity, founder effect and fitness costs are important bottle- necks that modulate the transfer of resistance from environmental microorganisms to pathogens.

Public health evolutionary biology of antimicrobial resistance: priorities for intervention

The three main processes shaping the evolutionary ecology of antibiotic resistance (AbR) involve the emergence, invasion and occupation by antibiotic-resistant genes of significant environments for human health. The process of emergence in complex bacterial populations is a high-frequency, continuous swarming of ephemeral combinatory genetic and epigenetic explorations inside cells and among cells, populations and communities, expanding in different environments (migration), creating the stochastic variation required for evolutionary progress. Invasion refers to the process by which AbR significantly increases in frequency in a given (invaded) environment, led by external invaders local multiplication and spread, or by endogenous conversion. Conversion occurs because of the spread of AbR genes from an exogenous resistant clone into an established (endogenous) bacterial clone(s) colonizing the environment; and/or because of dissemination of particular resistant genetic variants that emerged within an endogenous clonal population. Occupation of a given environment by a resistant variant means a permanent establishment of this organism in this environment, even in the absence of antibiotic selection. Specific interventions on emergence influence invasion, those acting on invasion also influence occupation and interventions on occupation determine emergence. Such interventions should be simultaneously applied, as they are not simple solutions to the complex problem of AbR.

Low rates of antimicrobial-resistant Enterobacteriaceae in wildlife in Taï National Park, Côte d'Ivoire, surrounded by villages with high prevalence of multiresistant ESBL-producing Escherichia coli in people and domestic animals

Antimicrobial resistance genes can be found in all ecosystems, including those where antibiotic selective pressure has never been exerted. We investigated resistance genes in a collection of faecal samples of wildlife (non-human primates, mice), people and domestic animals (dogs, cats) in Côte d’Ivoire; in the chimpanzee research area of Taï National Park (TNP) and adjacent villages. Single bacteria isolates were collected from antibiotic-containing agar plates and subjected to molecular analysis to detect Enterobacteriaceae isolates with plasmid-mediated genes of extended-spectrum beta-lactamases (ESBLs) and plasmid-mediated quinolone resistance (PMQR). While the prevalence of ESBL-producing E. coli in the villages was 27% in people (n = 77) and 32% in dogs (n = 38), no ESBL-producer was found in wildlife of TNP (n = 75). PMQR genes, mainly represented by qnrS1, were also present in human- and dog-originating isolates from the villages (36% and 42% in people and dogs, respectively), but no qnrS has been found in the park. In TNP, different variants of qnrB were detected in Citrobacter freundii isolates originating non-human primates and mice. In conclusion, ESBL and PMQR genes frequently found in humans and domestic animals in the villages were rather exceptional in wildlife living in the protected area. Although people enter the park, the strict biosecurity levels they are obliged to follow probably impede transmission of bacteria between them and wildlife.

What can we learn from fitness landscapes?

A combinatorially complete data set consists of studies of all possible combinations of a set of mutant sites in a gene or mutant alleles in a genome. Among the most robust conclusions from these studies is that epistasis between beneficial mutations often shows a pattern of diminishing returns, in which favorable mutations are less fit when combined than would be expected. Another robust inference is that the number of adaptive evolutionary paths is often limited to a relatively small fraction of the theoretical possibilities, owing largely to sign epistasis requiring evolutionary steps that would entail a decrease in fitness. Here we summarize these and other results while also examining issues that remain unresolved and future directions that seem promising.

Empirical fitness landscapes and the predictability of evolution

The genotype-fitness map (that is, the fitness landscape) is a key determinant of evolution, yet it has mostly been used as a superficial metaphor because we know little about its structure. This is now changing, as real fitness landscapes are being analysed by constructing genotypes with all possible combinations of small sets of mutations observed in phylogenies or in evolution experiments. In turn, these first glimpses of empirical fitness landscapes inspire theoretical analyses of the predictability of evolution. Here, we review these recent empirical and theoretical developments, identify methodological issues and organizing principles, and discuss possibilities to develop more realistic fitness landscape models.

Inferring fitness landscapes by regression produces biased estimates of epistasis

The genotype-fitness map plays a fundamental role in shaping the dynamics of evolution. However, it is difficult to directly measure a fitness landscape in practice, because the number of possible genotypes is astronomical. One approach is to sample as many genotypes as possible, measure their fitnesses, and fit a statistical model of the landscape that includes additive and pairwise interactive effects between loci. Here, we elucidate the pitfalls of using such regressions by studying artificial but mathematically convenient fitness landscapes. We identify two sources of bias inherent in these regression procedures, each of which tends to underestimate high fitnesses and overestimate low fitnesses. We characterize these biases for random sampling of genotypes as well as samples drawn from a population under selection in the Wright-Fisher model of evolutionary dynamics. We show that common measures of epistasis, such as the number of monotonically increasing paths between ancestral and derived genotypes, the prevalence of sign epistasis, and the number of local fitness maxima, are distorted in the inferred landscape. As a result, the inferred landscape will provide systematically biased predictions for the dynamics of adaptation. We identify the same biases in a computational RNA-folding landscape as well as regulatory sequence binding data treated with the same fitting procedure. Finally, we present a method to ameliorate these biases in some cases.

Impact of treatment strategies on cephalosporin and tetracycline resistance gene quantities in the bovine fecal metagenome

The study objective was to determine the effects of two treatment regimens on quantities of ceftiofur and tetracycline resistance genes in feedlot cattle. The two regimens were ceftiofur crystalline-free acid (CCFA) administered to either one or all steers within a pen and subsequent feeding/not feeding of therapeutic doses of chlortetracycline. A 26-day randomized controlled field trial was conducted on 176 steers. Real-time PCR was used to quantify bla_CMY-2, bla_CTX-M, tet(A), tet(B), and 16S rRNA gene copies/gram of feces from community DNA. A significant increase in ceftiofur resistance and a decrease in tetracycline resistance elements were observed among the treatment groups in which all steers received CCFA treatment, expressed as gene copies/gram of feces. Subsequent chlortetracycline administration led to rapid expansion of both ceftiofur and tetracycline resistance gene copies/gram of feces. Our data suggest that chlortetracycline is contraindicated when attempting to avoid expansion of resistance to critically important third-generation cephalosporins.

Live genomics for pathogen monitoring in public health

Whole genome analysis based on next generation sequencing (NGS) now represents an affordable framework in public health systems. Robust analytical pipelines of genomic data provides in short laps of time (hours) information about taxonomy, comparative genomics (pan-genome) and single polymorphisms profiles. Pathogenic organisms of interest can be tracked at the genomic level, allowing monitoring at one-time several variables including: epidemiology, pathogenicity, resistance to antibiotics, virulence, persistence factors, mobile elements and adaptation features. Such information can be obtained not only at large spectra, but also at the "local" level, such as in the event of a recurrent or emergency outbreak. This paper reviews the state of the art in infection diagnostics in the context of modern NGS methodologies. We describe how actuation protocols in a public health environment will benefit from a "streaming approach" (pipeline). Such pipeline would NGS data quality assessment, data mining for comparative analysis, searching differential genetic features, such as virulence, resistance persistence factors and mutation profiles (SNPs and InDels) and formatted "comprehensible" results. Such analytical protocols will enable a quick response to the needs of locally circumscribed outbreaks, providing information on the causes of resistance and genetic tracking elements for rapid detection, and monitoring actuations for present and future occurrences.