High fitness costs and instability of gene duplications reduce rates of evolution of new genes by duplication-divergence mechanisms

Comprehensive Assessment of Initial Adaptation of Extended-Spectrum β-Lactamase-Positive ST131 Escherichia coli to Carbapenem Exposure

BACKGROUND

It remains unclear how high-risk Escherichia coli lineages, like sequence type (ST) 131, initially adapt to carbapenem exposure in their progression to carbapenem resistance.

METHODS

Carbapenem mutation frequency was measured in multiple subclades of extended-spectrum β-lactamase (ESBL)-positive ST131 clinical isolates using a fluctuation assay followed by whole genome sequencing (WGS) characterization. Genomic, transcriptomic, and porin analyses of the ST131 C2/H30Rx isolate MB1860, under prolonged, increasing carbapenem exposure was performed using 2 experimental evolutionary platforms to measure fast versus slow adaptation.

RESULTS

All 13 ESBL-positive ST131 strains selected from a diverse (n = 184) ST131 bacteremia cohort had detectable ertapenem (ETP) mutational frequencies, with a positive correlation between initial ESBL gene copy number and mutation frequency (r = 0.87, P < 1e-5). WGS analysis of mutants showed that initial response to ETP exposure resulted in significant increases in ESBL gene copy numbers or mutations in Omp genes in the absence of ESBL gene amplification with subclade-specific associations. In both experimental evolutionary platforms, MB1860 responded to initial ETP exposure by increasing blaCTX-M-15 copy numbers via modular, IS26-mediated pseudocompound transposons (PCTns). Increased transcript level of genes present within the PCTn was a conserved expression signal in both experimental evolutionary platforms. Stable mutations in Omp encoding genes were detected only after prolonged increasing carbapenem exposure, consistent with clinical observations.

CONCLUSIONS

ESBL gene amplification is a conserved response to initial carbapenem exposure, especially within the high-risk ST131 C2/H30Rx subclade. Targeting such amplification could assist with mitigating carbapenem resistance development.

The fitness connection of antibiotic resistance

More than three decades ago multidrug-resistant (MDR) clones of the pathogens: Staphylococcus aureus, Escherichia coli, Klebsiella pneumoniae, Clostridioides difficile, Enterococcus faecium, Pseudomonas aeruginosa and Acinetobacter baumannii have started to disseminate across wide geographical areas. A characteristic feature of all these MDR lineages is the carriage of some mutations in the quinolone resistance-determining regions (QRDRs) of DNA gyrase and topoisomerase IV which besides conferring resistance to fluoroquinolones are associated with a fitness benefit. Several lines of evidence strongly suggest that extra fitness conferred by these mutations facilitated the dissemination of the international MDR lineages. MDR pathogens require extra energy to cover the fitness cost conferred by the excess antibiotic resistance gene cargo. However, extra energy generated by upgraded metabolic activity was demonstrated to increase the uptake of antibiotics enhancing susceptibility. Accordingly, MDR bacteria need additional positive fitness schemes which, similarly to the QRDR advantage, will not compromise resistance. Some of these, not clone-specific effects are large genomes, the carriage of low-cost plasmids, the transfer of plasmid genes to the chromosome, the application of weak promoters in integrons and various techniques for the economic control of the activity of the integrase enzyme including a highly sophisticated system in A. baumannii. These impacts - among others - will confer a fitness advantage promoting the spread of MDR pathogens. However, even the potential of extra fitness generated by the combined effect of various schemes is not without limit and virulence-related genes or less relevant antibiotic resistance gene cargoes will often be sacrificed to permit the acquisition of high-priority resistance determinants. Accordingly major MDR clone strains are usually less virulent than susceptible isolates. In summary, a fitness approach to the research of antibiotic resistance is very useful since the fitness status of MDR bacteria seem to profoundly impact the capacity to disseminate in the healthcare setting.

A sequence context-based approach for classifying tumor structural variants without paired normal samples

Although several recent studies have characterized structural variants (SVs) in germline and cancer genomes independently, the genomic contexts of these SVs have not been comprehensively compared. We examined similarities and differences between 2 million germline and 115 thousand tumor SVs from a cohort of 963 patients from The Cancer Genome Atlas. We found significant differences in features related to their genomic sequences and localization that suggest differences between SV-generating processes and selective pressures. For example, our results show that features linked to transposon-mediated processes are associated with germline SVs, while somatic SVs more frequently show features characteristic of chromoanagenesis. These genomic differences enabled us to develop a classifier-the Germline and Tumor Structural Variant or "the great GaTSV" -that accurately distinguishes between germline and cancer SVs in tumor samples that lack a matched normal sample.

Rational Design Assisted by Evolutionary Engineering Allows (De)Construction and Optimization of Complex Phenotypes in Pseudomonas putida KT2440

Beyond the rational construction of genetic determinants to encode target functions, complex phenotype engineering requires the contextualisation of their expression within the metabolic and genetic background of the host strain. Furthermore, wherever metabolic complexity is involved, phenotype engineering demands standard, reliable, plug-and-play tools. We introduce GENIO (GENome Integration and fitness Optimization platform for Pseudomonas putida), a framework to optimise genetic circuit performance by means of (i) chromosome-location-based differential gene expression and (ii) subsequent fitness improvement through evolutionary engineering if needed. Using gene expression strength and cell-to-cell variation, we characterised 10 P. putida chromosomal loci (ppLPS) to show that genome context rather than distance to ORI is the main factor driving differential expression performance. We further contextualised ppLPS gene expression against well-known chromosomal integration sites and plasmids displaying different copy numbers. GENIO supports comprehensive exploration of the gene expression space across P. putida's genome while unlocking performance optimization of complex heterologous metabolic pathways through evolutionary engineering. To demonstrate the usability of GENIO, we restored P. putida's aromatic hydrocarbon metabolism by (de)constructing the toluene/m-xylene catabolic pathway coded in the pWW0 plasmid. We also showed that engineering complex phenotypes requires accurate contextualisation of the synthetic pathways involved, a process that benefits from biological robustness.

Machine learning detection of heteroresistance in Escherichia coli

BACKGROUND

Heteroresistance (HR) is a significant type of antibiotic resistance observed for several bacterial species and antibiotic classes where a susceptible main population contains small subpopulations of resistant cells. Mathematical models, animal experiments and clinical studies associate HR with treatment failure. Currently used susceptibility tests do not detect heteroresistance reliably, which can result in misclassification of heteroresistant isolates as susceptible which might lead to treatment failure. Here we examined if whole genome sequence (WGS) data and machine learning (ML) can be used to detect bacterial HR.

METHODS

We classified 467 Escherichia coli clinical isolates as HR or non-HR to the often used β-lactam/inhibitor combination piperacillin-tazobactam using pre-screening and Population Analysis Profiling tests. We sequenced the isolates, assembled the whole genomes and created a set of predictors based on current knowledge of HR mechanisms. Then we trained several machine learning models on 80% of this data set aiming to detect HR isolates. We compared performance of the best ML models on the remaining 20% of the data set with a baseline model based solely on the presence of β-lactamase genes. Furthermore, we sequenced the resistant sub-populations in order to analyse the genetic mechanisms underlying HR.

FINDINGS

The best ML model achieved 100% sensitivity and 84.6% specificity, outperforming the baseline model. The strongest predictors of HR were the total number of β-lactamase genes, β-lactamase gene variants and presence of IS elements flanking them. Genetic analysis of HR strains confirmed that HR is caused by an increased copy number of resistance genes via gene amplification or plasmid copy number increase. This aligns with the ML model's findings, reinforcing the hypothesis that this mechanism underlies HR in Gram-negative bacteria.

INTERPRETATION

We demonstrate that a combination of WGS and ML can identify HR in bacteria with perfect sensitivity and high specificity. This improved detection would allow for better-informed treatment decisions and potentially reduce the occurrence of treatment failures associated with HR.

FUNDING

Funding provided to DIA from the Swedish Research Council (2021-02091) and NIH (1U19AI158080-01).

Heteroresistance in Enterobacter cloacae complex caused by variation in transient gene amplification events

Heteroresistance (HR) in bacteria describes a subpopulational phenomenon of antibiotic resistant cells of a generally susceptible population. Here, we investigated the molecular mechanisms and phenotypic characteristics underlying HR to ceftazidime (CAZ) in a clinical Enterobacter cloacae complex strain (ECC). We identified a plasmid-borne gene duplication-amplification (GDA) event of a region harbouring an ampC gene encoding a β-lactamase blaDHA-1 as the key determinant of HR. Individual colonies exhibited variations in the copy number of the genes resulting in resistance level variation which correlated with growth onset (lag times) and growth rates in the presence of CAZ. GDA copy number heterogeneity occurred within single resistant colonies, demonstrating heterogeneity of GDA on the single-cell level. The interdependence between GDA, lag time and antibiotic treatment and the strong plasticity underlying HR underlines the high risk for misdetection of antimicrobial HR and subsequent treatment failure.

Structural and biochemical insights of xylose MFS and SWEET transporters in microbial cell factories: challenges to lignocellulosic hydrolysates fermentation

The production of bioethanol from lignocellulosic biomass requires the efficient conversion of glucose and xylose to ethanol, a process that depends on the ability of microorganisms to internalize these sugars. Although glucose transporters exist in several species, xylose transporters are less common. Several types of transporters have been identified in diverse microorganisms, including members of the Major Facilitator Superfamily (MFS) and Sugars Will Eventually be Exported Transporter (SWEET) families. Considering that Saccharomyces cerevisiae lacks an effective xylose transport system, engineered yeast strains capable of efficiently consuming this sugar are critical for obtaining high ethanol yields. This article reviews the structure-function relationship of sugar transporters from the MFS and SWEET families. It provides information on several tools and approaches used to identify and characterize them to optimize xylose consumption and, consequently, second-generation ethanol production.

Comprehensive Assessment of Initial Adaptation of ESBL Positive ST131 Escherichia coli to Carbapenem Exposure

Background

It remains unclear how high-risk Escherichia coli lineages, like sequence type (ST) 131, initially adapt to carbapenem exposure in their progression to becoming carbapenem resistant.

Methods

Carbapenem mutation frequency was measured in multiple subclades of extended-spectrum β-lactamase (ESBL) positive ST131 clinical isolates using a fluctuation assay followed by whole genome sequencing (WGS) characterization. Genomic, transcriptomic, and porin analyses of ST131 C2/ H 30Rx isolate, MB1860, under prolonged, increasing carbapenem exposure was performed using two distinct experimental evolutionary platforms to measure fast vs. slow adaptation.

Results

All thirteen ESBL positive ST131 strains selected from a diverse (n=184) ST131 bacteremia cohort had detectable ertapenem (ETP) mutational frequencies with a statistically positive correlation between initial ESBL gene copy number and mutation frequency (r = 0.87, P -value <1e-5). WGS analysis of mutants showed initial response to ETP exposure resulted in significant increases in ESBL gene copy numbers or mutations in outer membrane porin (Omp) encoding genes in the absence of ESBL gene amplification with subclade specific associations. In both experimental evolutionary platforms, MB1860 responded to initial ETP exposure by increasing bla CTX-M-15 copy numbers via modular, insertion sequence 26 (IS 26 ) mediated pseudocompound transposons (PCTns). Transposase activity driven by PCTn upregulation was a conserved expression signal in both experimental evolutionary platforms. Stable mutations in Omp encoding genes were detected only after prolonged increasing carbapenem exposure consistent with clinical observations.

Conclusions

ESBL gene amplification is a conserved response to initial carbapenem exposure, especially within the high-risk ST131 C2/ H 30Rx subclade. Targeting such amplification could assist with mitigating carbapenem resistance development.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

Comparison of germline and somatic structural variants in cancers reveal systematic differences in variant generating and selection processes

Although several recent studies have characterized structural variants (SVs) in germline and cancer genomes, the features of SVs in these different contexts have not been directly compared. We examined similarities and differences between 2 million germline and 115 thousand tumor SVs from a cohort of 963 patients from The Cancer Genome Atlas (TCGA). We found significant differences in features related to their genomic sequences and localization that suggest differences between SV-generating processes and selective pressures. For example, we found that transposon-mediated processes shape germline much more than somatic SVs, while somatic SVs more frequently show features characteristic of chromoanagenesis. These differences were extensive enough to enable us to develop a classifier - "the great GaTSV" - that accurately distinguishes between germline and cancer SVs in tumor samples that lack a matched normal sample.

Generation and maintenance of the circularized multimeric IS26-associated translocatable unit encoding multidrug resistance

In gram-negative bacteria, IS26 often exists in multidrug resistance (MDR) regions, forming a pseudocompound transposon (PCTn) that can be tandemly amplified. It also generates a circular intermediate called the "translocatable unit (TU)", but the TU has been detected only by PCR. Here, we demonstrate that in a Klebsiella pneumoniae MDR clone, mono- and multimeric forms of the TU were generated from the PCTn in a preexisting MDR plasmid where the inserted form of the TU was also tandemly amplified. The two modes of amplification were reproduced by culturing the original clone under antimicrobial selection pressure, and the amplified state was maintained in the absence of antibiotics. Mono- and multimeric forms of the circularized TU were generated in a RecA-dependent manner from the tandemly amplified TU, which can be generated in RecA-dependent and independent manners. These findings provide novel insights into the dynamic processes of genome amplification in bacteria.

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Vancomycin heteroresistance caused by unstable tandem amplifications of the vanM gene cluster on linear conjugative plasmids in a clinical Enterococcus faecium

Vancomycin heteroresistance is prone to missed detection and poses a risk of clinical treatment failure. We encountered one clinical Enterococcus faecium strain, SRR12, that carried a complete vanM gene cluster but was determined as susceptible to vancomycin using the broth microdilution method. However, distinct subcolonies appeared within the clear zone of inhibition in the E-test assay, one of which, named SRR12-v1, showed high-level resistance to vancomycin. SRR12 was confirmed as heteroresistant to vancomycin using population analysis profiling and displayed "revive" growth curves with a lengthy lag phase of over 13 hours when exposed to 2-32 mg/L vancomycin. The resistant subcolony SRR12-v1 was found to carry an identical vanM gene cluster to that of SRR12 but a significantly increased vanM copy number in the genome. Long-read whole genome sequencing revealed that a one-copy vanM gene cluster was located on a pELF1-like linear plasmid in SRR12. In comparison, tandem amplification of the vanM gene cluster jointed with IS1216E was seated on a linear plasmid in the genome of SRR12-v1. These amplifications of the vanM gene cluster were demonstrated as unstable and would decrease accompanied by fitness reversion after serial passaging for 50 generations under increasing vancomycin pressure or without antibiotic pressure but were relatively stable under constant vancomycin pressure. Further, vanM resistance in resistant variants was verified to be carried by conjugative plasmids with variable sizes using conjugation assays and S1-pulsed field gel electrophoresis blotting, suggesting the instability/flexibility of vanM cluster amplification in the genome and an increased risk of vanM resistance dissemination.

Heteroresistance to piperacillin/tazobactam in Klebsiella pneumoniae is mediated by increased copy number of multiple β-lactamase genes

Background

Piperacillin/tazobactam is a β-lactam/β-lactamase inhibitor combination with a broad spectrum of activity that is often used as empirical and/or targeted therapy among hospitalized patients. Heteroresistance (HR) is a form of antibiotic resistance in which a minority population of resistant cells coexists with a majority susceptible population that has been found to be a cause of antibiotic treatment failure in murine models.

Objectives

To determine the prevalence of HR and mechanisms of HR to piperacillin/tazobactam among Klebsiella pneumoniae bloodstream infection (BSI) isolates.

Materials

From July 2018 to June 2021, K. pneumoniae piperacillin/tazobactam-susceptible BSI isolates were collected from two tertiary hospitals in Atlanta, GA, USA. Only first isolates from each patient per calendar year were included. Population analysis profiling (PAP) and WGS were performed to identify HR and its mechanisms.

Results

Among 423 K. pneumoniae BSI isolates collected during the study period, 6% (25/423) were found to be HR with a subpopulation surviving above the breakpoint. WGS of HR isolates grown in the presence of piperacillin/tazobactam at concentrations 8-fold that of the MIC revealed copy number changes of plasmid-located β-lactamase genes blaCTX-M-15, blaSHV33, blaOXA-1 and blaTEM-1 by tandem gene amplification or plasmid copy number increase.

Conclusions

Prevalence of HR to piperacillin/tazobactam among bloodstream isolates was substantial. The HR phenotype appears to be caused by tandem amplification of β-lactamase genes found on plasmids or plasmid copy number increase. This raises the possibility of dissemination of HR through horizontal gene transfer and requires further study.

Plasmid-mediated phenotypic noise leads to transient antibiotic resistance in bacteria

The rise of antibiotic resistance is a critical public health concern, requiring an understanding of mechanisms that enable bacteria to tolerate antimicrobial agents. Bacteria use diverse strategies, including the amplification of drug-resistance genes. In this paper, we showed that multicopy plasmids, often carrying antibiotic resistance genes in clinical bacteria, can rapidly amplify genes, leading to plasmid-mediated phenotypic noise and transient antibiotic resistance. By combining stochastic simulations of a computational model with high-throughput single-cell measurements of blaTEM-1 expression in Escherichia coli MG1655, we showed that plasmid copy number variability stably maintains populations composed of cells with both low and high plasmid copy numbers. This diversity in plasmid copy number enhances the probability of bacterial survival in the presence of antibiotics, while also rapidly reducing the burden of carrying multiple plasmids in drug-free environments. Our results further support the tenet that multicopy plasmids not only act as vehicles for the horizontal transfer of genetic information between cells but also as drivers of bacterial adaptation, enabling rapid modulation of gene copy numbers. Understanding the role of multicopy plasmids in antibiotic resistance is critical, and our study provides insights into how bacteria can transiently survive lethal concentrations of antibiotics.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Bacteria can compensate the fitness costs of amplified resistance genes via a bypass mechanism

Antibiotic heteroresistance is a phenotype in which a susceptible bacterial population includes a small subpopulation of cells that are more resistant than the main population. Such resistance can arise by tandem amplification of DNA regions containing resistance genes that in single copy are not sufficient to confer resistance. However, tandem amplifications often carry fitness costs, manifested as reduced growth rates. Here, we investigated if and how these fitness costs can be genetically ameliorated. We evolved four clinical isolates of three bacterial species that show heteroresistance to tobramycin, gentamicin and tetracyclines at increasing antibiotic concentrations above the minimal inhibitory concentration (MIC) of the main susceptible population. This led to a rapid enrichment of resistant cells with up to an 80-fold increase in the resistance gene copy number, an increased MIC, and severely reduced growth rates. When further evolved in the presence of antibiotic, these strains acquired compensatory resistance mutations and showed a reduction in copy number while maintaining high-level resistance. A deterministic model indicated that the loss of amplified units was driven mainly by their fitness costs and that the compensatory mutations did not affect the loss rate of the gene amplifications. Our findings suggest that heteroresistance mediated by copy number changes can facilitate and precede the evolution towards stable resistance.

Comparative genomics reveals probable adaptations for xylose use in Thermoanaerobacterium saccharolyticum

Second-generation ethanol, a promising biofuel for reducing greenhouse gas emissions, faces challenges due to the inefficient metabolism of xylose, a pentose sugar. Overcoming this hurdle requires exploration of genes, pathways, and organisms capable of fermenting xylose. Thermoanaerobacterium saccharolyticum is an organism capable of naturally fermenting compounds of industrial interest, such as xylose, and understanding evolutionary adaptations may help to bring novel genes and information that can be used for industrial yeast, increasing production of current bio-platforms. This study presents a deep evolutionary study of members of the firmicutes clade, focusing on adaptations in Thermoanaerobacterium saccharolyticum that may be related to overall fermentation metabolism, especially for xylose fermentation. One highlight is the finding of positive selection on a xylose-binding protein of the xylFGH operon, close to the annotated sugar binding site, with this protein already being found to be expressed in xylose fermenting conditions in a previous study. Results from this study can serve as basis for searching for candidate genes to use in industrial strains or to improve Thermoanaerobacterium saccharolyticum as a new microbial cell factory, which may help to solve current problems found in the biofuels' industry.

Role of riboflavin biosynthesis gene duplication and transporter in Aeromonas salmonicida virulence in marine teleost fish

Active flavins derived from riboflavin (vitamin B₂) are essential for life. Bacteria biosynthesize riboflavin or scavenge it through uptake systems, and both mechanisms may be present. Because of riboflavin's critical importance, the redundancy of riboflavin biosynthetic pathway (RBP) genes might be present. Aeromonas salmonicida, the aetiological agent of furunculosis, is a pathogen of freshwater and marine fish, and its riboflavin pathways have not been studied. This study characterized the A. salmonicida riboflavin provision pathways. Homology search and transcriptional orchestration analysis showed that A. salmonicida has a main riboflavin biosynthetic operon that includes ribD, ribE1, ribBA, and ribH genes. Outside the main operon, putative duplicated genes ribA, ribB and ribE, and a ribN riboflavin importer encoding gene, were found. Monocistronic mRNA ribA, ribB and ribE2 encode for their corresponding functional riboflavin biosynthetic enzyme. While the product of ribBA conserved the RibB function, it lacked the RibA function. Likewise, ribN encodes a functional riboflavin importer. Transcriptomics analysis indicated that external riboflavin affected the expression of a relatively small number of genes, including a few involved in iron metabolism. ribB was downregulated in response to external riboflavin, suggesting negative feedback. Deletion of ribA, ribB and ribE1 showed that these genes are required for A. salmonicida riboflavin biosynthesis and virulence in Atlantic lumpfish (Cyclopterus lumpus). A. salmonicida riboflavin auxotrophic attenuated mutants conferred low protection to lumpfish against virulent A. salmonicida. Overall, A. salmonicida has multiple riboflavin endowment forms, and duplicated riboflavin provision genes are critical for A. salmonicida infection.

Role of blaTEM and OmpC in the piperacillin-tazobactam resistance evolution by E. coli in patients with complicated intra-abdominal infection

Piperacillin-tazobactam resistance (P/T-R) is increasingly reported among Escherichia coli isolates. Although in vitro experiments have suggested that blaTEM gene plays a key role in the P/T-R acquisition, no clinical in vivo study has yet confirmed the role of blaTEM or other genes. Therefore, we aimed to identify the mechanisms underlying P/T-R by following up patients with E. coli complicated intra-abdominal infections (cIAI) who experienced P/T treatment failure. Four pairs of strains, clonally related from four patients, were isolated both before and after treatment with P/T dosed at 4 g/0.5 g intravenously. The P/T MIC was tested using broth microdilution, and β-lactamase activity was determined in these isolates. Whole-genome sequencing (WGS) was performed to decipher the role of blaTEM and other genes associated with P/T-R. Changes in the outer membrane protein (OMP) profile were analyzed using SDS-PAGE, and blaTEM and ompC transcription levels were measured by RT-qPCR. In addition, in vitro competition fitness was performed between each pairs of strains (P/T-susceptible vs. P/T-resistant). We found a higher copy number of blaTEM gene in P/T-R isolates, generated by three different genetic events: (1) IS26-mediated duplication of the blaTEM gene, (2) generation of a small multicopy plasmid (ColE-like) carrying blaTEM, and (3) adaptive evolution via reduction of plasmid size, leading to a higher plasmid copy number. Moreover, two P/T-R strains showed reduced expression of OmpC. This study describes the mechanisms involved in the acquisition of P/T-R by E. coli in patients with cIAI. The understanding of P/T-R evolution is crucial for effectively treating infected patients and preventing the spread of resistant microorganisms.

Cancer genomes tolerate deleterious coding mutations through somatic copy number amplifications of wild-type regions

Cancers evolve under the accumulation of thousands of somatic mutations and chromosomal aberrations. While most coding mutations are deleterious, almost all protein-coding genes lack detectable signals of negative selection. This raises the question of how tumors tolerate such large amounts of deleterious mutations. Using 8,690 tumor samples from The Cancer Genome Atlas, we demonstrate that copy number amplifications frequently cover haploinsufficient genes in mutation-prone regions. This could increase tolerance towards the deleterious impact of mutations by creating safe copies of wild-type regions and, hence, protecting the genes therein. Our findings demonstrate that these potential buffering events are highly influenced by gene functions, essentiality, and mutation impact and that they occur early during tumor evolution. We show how cancer type-specific mutation landscapes drive copy number alteration patterns across cancer types. Ultimately, our work paves the way for the detection of novel cancer vulnerabilities by revealing genes that fall within amplifications likely selected during evolution to mitigate the effect of mutations.

Gene loss during a transition to multicellularity

Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited novel interactions, suggesting that gene loss could play a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest gene loss, in addition to gene gain and co-option, may be important for the evolution developmental complexity.

How to Recruit a Promiscuous Enzyme to Serve a New Function

Promiscuous enzymes can be recruited to serve new functions when a genetic or environmental change makes catalysis of a novel reaction important for fitness or even survival. Subsequently, gene duplication and divergence can lead to evolution of an efficient and specialized new enzyme. Every organism likely has thousands of promiscuous enzyme activities that provide a vast reservoir of catalytic potential. However, much of this potential may not be accessible. We compiled kinetic parameters for promiscuous reactions catalyzed by 108 enzymes. The median value of k_cat/K_M is a very modest 31 M^-1 s^-1. Based upon the fluxes through metabolic pathways in E. coli, we estimate that many, if not most, promiscuous activities are too inefficient to impact fitness. However, mutations can elevate the level of an insufficient promiscuous activity by increasing enzyme expression, improving k_cat/K_M, or altering concentrations of the promiscuous and native substrates and allosteric regulators. Particularly in large bacterial populations, stochastic mutations may provide a viable pathway for recruitment of even inefficient promiscuous activities.

Class 1 integrons and multiple mobile genetic elements in clinical isolates of the Klebsiella pneumoniae complex from a tertiary hospital in eastern China

Background

The emergence of highly drug-resistant K. pneumoniae, has become a major public health challenge. In this work, we aim to investigate the diversity of species and sequence types (STs) of clinical Klebsiella isolates and to characterize the prevalence and structure of class 1 integrons.

Methods

Based on the whole genome sequencing, species identification was performed by 16S rRNA gene homology and average nucleotide identity (ANI) analysis. STs were determined in accordance with the international MLST schemes for K. pneumoniae and K. variicola. Integron characterization and comparative genomic analysis were performed using various bioinformatic tools.

Results

Species identification showed that the 167 isolates belonged to four species: K. pneumoniae, K. variicola subsp. variicola, K. quasipneumoniae and K. aerogenes. Thirty-six known and 5 novel STs were identified in K. pneumoniae, and 10 novel STs were identified in K. variicola subsp. variicola. Class 1 integrons were found in 57.49% (96/167) of the isolates, and a total of 169 resistance gene cassettes encoding 19 types of resistance genes, including carbapenem resistance gene (bla _IPM-4) and class D β-lactamases gene (bla _OXA-1 and bla _OXA-10), were identified. Among the 17 complete genomes, 29 class 1 integrons from 12 groups were found, only 1 group was encoded on chromosomes. Interestingly, one plasmid (pKP167-261) carrying two copies of approximately 19-kb IS26-Int1 complex resistance region that contains an integron and a multidrug resistance gene fragment.

Conclusion

The results of this work demonstrated that the species and STs of the clinical Klebsiella isolates were more complex by the whole genome sequence analysis than by the traditional laboratory methods. Finding of the new structure of MGEs related to the resistance genes indicates the great importance of deeply exploring the molecular mechanisms of bacterial multidrug resistance.

Draft genome of six Cuban Anolis lizards and insights into genetic changes during their diversification

Background

Detecting genomic variants and their accumulation processes during species diversification and adaptive radiation is important for understanding the molecular and genetic basis of evolution. Anolis lizards in the West Indies are good models for studying evolutionary mechanisms because of the repeated evolution of their morphology and the ecology. We performed de novo genome assembly of six Cuban Anolis lizards with different ecomorphs and thermal habitats (Anolis isolepis, Anolis allisoni, Anolis porcatus, Anolis allogus, Anolis homolechis, and Anolis sagrei). We carried out a comparative analysis of these genome assemblies to investigate the genetic changes that occurred during their diversification.

Results

We reconstructed novel draft genomes with relatively long scaffolds and high gene completeness, with the scaffold N50 ranging from 5.56 to 39.79 Mb and vertebrate Benchmarking Universal Single-Copy Orthologs completeness ranging from 77.5% to 86.9%. Comparing the repeat element compositions and landscapes revealed differences in the accumulation process between Cuban trunk-crown and trunk-ground species and separate expansions of several families of LINE in each Cuban trunk-ground species. Duplicated gene analysis suggested that the proportional differences in duplicated gene numbers among Cuban Anolis lizards may be associated with differences in their habitat ranges. Additionally, Pairwise Sequentially Markovian Coalescent analysis suggested that the effective population sizes of each species may have been affected by Cuba’s geohistory.

Conclusions

We provide draft genomes of six Cuban Anolis lizards and detected species and lineage-specific transposon accumulation and gene copy number changes that may be involved in adaptive evolution. The change processes in the past effective population size was also estimated, and the factors involved were inferred. These results provide new insights into the genetic basis of Anolis lizard diversification and are expected to serve as a stepping stone for the further elucidation of their diversification mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02086-7.

Adaptive evolvability through direct selection instead of indirect, second-order selection

Can evolvability itself be the product of adaptive evolution? To answer this question is challenging, because any DNA mutation that alters only evolvability is subject to indirect, "second order" selection on the future effects of this mutation. Such indirect selection is weaker than "first-order" selection on mutations that alter fitness, in the sense that it can operate only under restrictive conditions. Here I discuss a route to adaptive evolvability that overcomes this challenge. Specifically, a recent evolution experiment showed that some mutations can enhance both fitness and evolvability through a combination of direct and indirect selection. Unrelated evidence from gene duplication and the evolution of gene regulation suggests that mutations with such dual effects may not be rare. Through such mutations, evolvability may increase at least in part because it provides an adaptive advantage. These observations suggest a research program on the adaptive evolution of evolvability, which aims to identify such mutations and to disentangle their direct fitness effects from their indirect effects on evolvability. If evolvability is itself adaptive, Darwinian evolution may have created more than life's diversity. It may also have helped create the very conditions that made the success of Darwinian evolution possible.

Systematic Analysis of Mobile Genetic Elements Mediating β-Lactamase Gene Amplification in Noncarbapenemase-Producing Carbapenem-Resistant Enterobacterales Bloodstream Infections

Noncarbapenemase-producing carbapenem-resistant Enterobacterales (non-CP-CRE) are increasingly recognized as important contributors to prevalent carbapenem-resistant Enterobacterales (CRE) infections. However, there is limited understanding of mechanisms underlying non-CP-CRE causing invasive disease. Long- and short-read whole-genome sequencing was used to elucidate carbapenem nonsusceptibility determinants in Enterobacterales bloodstream isolates at MD Anderson Cancer Center in Houston, Texas. We investigated carbapenem nonsusceptible Enterobacterales (CNSE) mechanisms (i.e., isolates with carbapenem intermediate resistance phenotypes or greater) through a combination of phylogenetic analysis, antimicrobial resistance gene detection/copy number quantification, porin assessment, and mobile genetic element (MGE) characterization. Most CNSE isolates sequenced were non-CP-CRE (41/79; 51.9%), whereas 25.3% (20/79) were Enterobacterales with intermediate susceptibility to carbapenems (CIE), and 22.8% (18/79) were carbapenemase-producing Enterobacterales (CPE). Statistically significant copy number variants (CNVs) of extended-spectrum β-lactamase (ESBL) genes (Wilcoxon Test; P-value < 0.001) were present in both non-CP-CR E. coli (median CNV = 2.6×; n = 17) and K. pneumoniae (median CNV = 3.2×, n = 17). All non-CP-CR E. coli and K. pneumoniae had predicted reduced expression of at least one outer membrane porin gene (i.e., ompC/ompF or ompK36/ompK35). Completely resolved CNSE genomes revealed that IS26 and ISEcp1 structures harboring blaCTX-M variants along with other antimicrobial resistance elements were associated with gene amplification, occurring in mostly IncFIB/IncFII plasmid contexts. MGE-mediated β-lactamase gene amplifications resulted in either tandem arrays, primarily mediated by IS26 translocatable units, or segmental duplication, typically due to ISEcp1 transposition units. Non-CP-CRE strains were the most common cause of CRE bacteremia with carbapenem nonsusceptibility driven by concurrent porin loss and MGE-mediated amplification of blaCTX-M genes. IMPORTANCE Carbapenem-resistant Enterobacterales (CRE) are considered urgent antimicrobial resistance (AMR) threats. The vast majority of CRE research has focused on carbapenemase-producing Enterobacterales (CPE) even though noncarbapenemase-producing CRE (non-CP-CRE) comprise 50% or more of isolates in some surveillance studies. Thus, carbapenem resistance mechanisms in non-CP-CRE remain poorly characterized. To address this problem, we applied a combination of short- and long-read sequencing technologies to a cohort of CRE bacteremia isolates and used these data to unravel complex mobile genetic element structures mediating β-lactamase gene amplification. By generating complete genomes of 65 carbapenem nonsusceptible Enterobacterales (CNSE) covering a genetically diverse array of isolates, our findings both generate novel insights into how non-CP-CRE overcome carbapenem treatments and provide researchers scaffolds for characterization of their own non-CP-CRE isolates. Improved recognition of mechanisms driving development of non-CP-CRE could assist with design and implementation of future strategies to mitigate the impact of these increasingly recognized AMR pathogens.

Genetic redundancy of 4-hydroxybenzoate 3-hydroxylase genes ensures the catabolic safety of Pigmentiphaga sp. H8 in 3-bromo-4-hydroxybenzoate-contaminated habitats

Genetic redundancy is prevalent in organisms and plays important roles in the evolution of biodiversity and adaptation to environmental perturbation. However, selective advantages of genetic redundancy in overcoming metabolic disturbance due to structural analogues have received little attention. Here, functional divergence of the three 4-hydroxybenzoate 3-hydroxylase (PHBH) genes (phbh1~3) was found in Pigmentiphaga sp. strain H8. The genes phbh1/phbh2 were responsible for 3-bromo-4-hydroxybenzoate (3-Br-4-HB, an anthropogenic pollutant) catabolism, whereas phbh3 was primarily responsible for 4-hydroxybenzoate (4-HB, a natural intermediate of lignin) catabolism. 3-Br-4-HB inhibited 4-HB catabolism by competitively binding PHBH3, and was toxic to strain H8 cells especially at high concentrations. The existence of phbh1/phbh2 not only enabled strain H8 to utilize 3-Br-4-HB, but also ensured the catabolic safety of 4-HB. Molecular docking and site-directed mutagenesis analyses revealed that Val199 and Phe384 of PHBH1/PHBH2 were required for the hydroxylation activity towards 3-Br-4-HB. Phylogenetic analysis indicated that phbh1 and phbh2 originated from a common ancestor and evolved specifically in strain H8 to adapt to 3-Br-4-HB-contaminated habitats, whereas phbh3 evolved independently. This study deepens our understanding of selective advantages of genetic redundancy in prokaryote's metabolic robustness and reveals the factors driving the divergent evolution of redundant genes in adaptation to environmental perturbation. This article is protected by copyright. All rights reserved.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Evolutionary Trajectories toward High-Level β-Lactam/β-Lactamase Inhibitor Resistance in the Presence of Multiple β-Lactamases

β-Lactam antibiotics are the first choice for the treatment of most bacterial infections. However, the increased prevalence of β-lactamases, in particular extended-spectrum β-lactamases, in pathogenic bacteria has severely limited the possibility of using β-lactam treatments. Combining β-lactam antibiotics with β-lactamase inhibitors can restore treatment efficacy by negating the effect of the β-lactamase and has become increasingly important against infections caused by β-lactamase-producing strains. Not surprisingly, bacteria with resistance to even these combinations have been found in patients. Studies on the development of bacterial resistance to β-lactam/β-lactamase inhibitor combinations have focused mainly on the effects of single, chromosomal or plasmid-borne, β-lactamases. However, clinical isolates often carry more than one β-lactamase in addition to multiple other resistance genes. Here, we investigate how the evolutionary trajectories of the development of resistance to three commonly used β-lactam/β-lactamase inhibitor combinations, ampicillin-sulbactam, piperacillin-tazobactam, and ceftazidime-avibactam, were affected by the presence of three common β-lactamases, TEM-1, CTX-M-15, and OXA-1. First-step resistance was due mainly to extensive gene amplifications of one or several of the β-lactamase genes where the amplification pattern directly depended on the respective drug combination. Amplifications also served as a stepping-stone for high-level resistance in combination with additional mutations that reduced drug influx or mutations in the β-lactamase gene bla_CTX-M-15. This illustrates that the evolutionary trajectories of resistance to β-lactam/β-lactamase inhibitor combinations are strongly influenced by the frequent and transient nature of gene amplifications and how the presence of multiple β-lactamases shapes the evolution to higher-level resistance.

Comparison of Two Distinct Subpopulations of Klebsiella pneumoniae ST16 Co-Occurring in a Single Patient

The higher resistance rate to ceftazidime-avibactam (CZA) is mainly related to carbapenem resistance, especially New Delhi metallo-β-lactamase (NDM). The CZA-susceptible Klebsiella pneumoniae (K191663) and the later CZA-resistant isolates (K191724, K191725, K191773) co-producing NDM-4 and OXA-181 were obtained from the same hospitalized patient returning from Vietnam. Our study aims to elucidate the diversity of K. pneumoniae ST16 through comparative analysis of whole-genome sequencing (WGS) data and identify the potential evolution of plasmids by sequencing longitudinal clinical isolates during antibiotic treatment. Firstly, multilocus sequence typing analysis and phylogenic analysis suggested that these strains belong to the two lineages of K. pneumoniae ST16. Surprisingly, the CZA-resistant strains were closely related to K. pneumoniae ST16 described in South Korea, instead of the bla_NDM-4- or bla_OXA-181-carrying ST16 reported in Vietnam. Secondly, bla_NDM-4, bla_TEM-1B, and rmtB co-existed on a self-conjugative IncFII(Yp)-like plasmid, which played a significant role in CZA resistance. It could transfer into the recipient Escherichia coli J53 at high frequency, indicating the risk of mobile carbapenemases. In addition, the loss of 12-kbp fragment occurred in bla_NDM-4-positive isolate (K191773), which was likely caused by insertion sequence-mediated homologous recombination. Last but not least, as a repressor of acrAB operon system, acrR was truncated by a frameshift mutation in K191663. Thus, our study provided baseline information for monitoring the occurrence and development of bacterial resistance.

IMPORTANCE As a leading health care-acquired infection pathogen, Klebsiella pneumoniae is threatening a large number of inpatients due to its diverse antibiotic resistance and virulence factors. Heretofore, with a growing number of reports about the coexistence of several carbapenemases in carbapenem-resistant K. pneumoniae (CRKP), epidemiologic surveillance has been strengthened. Nevertheless, the nosocomial outbreaks by CRKP ST16 are gradually increasing worldwide. Our study provides a deeper insight into the diversification of clinical isolates of CRKP ST16 in China. In addition, the comparison analysis of resistant plasmids may reveal the transmission of carbapenemase-encoding genes. Furthermore, our study also highlights the importance of longitudinal specimen collection and continuous monitoring during the treatment, which play a crucial role in understanding the development of antibiotic resistance and the evolution of resistance plasmids.

The Role of Antibiotic Resistance Genes in the Fitness Cost of Multiresistance Plasmids

By providing the bacterial cell with protection against several antibiotics at once, multiresistance plasmids have an evolutionary advantage in situations where antibiotic treatments are common, such as in hospital environments. However, resistance plasmids can also impose fitness costs on the bacterium in the absence of antibiotics, something that may limit their evolutionary success. The underlying mechanisms and the possible contribution of resistance genes to such costs are still largely not understood. Here, we have specifically investigated the contribution of plasmid-borne resistance genes to the reduced fitness of the bacterial cell. The pUUH239.2 plasmid carries 13 genes linked to antibiotic resistance and reduces bacterial fitness by 2.9% per generation. This cost is fully ameliorated by the removal of the resistance cassette. While most of the plasmid-borne resistance genes individually were cost-free, even when overexpressed, two specific gene clusters were responsible for the entire cost of the plasmid: the extended-spectrum-β-lactamase gene bla_CTX-M-15 and the tetracycline resistance determinants tetAR. The bla_CTX-M-15 cost was linked to the signal peptide that exports the β-lactamase into the periplasm, and replacement with an alternative signal peptide abolished the cost. Both the tetracycline pump TetA and its repressor TetR conferred a cost on the host cell, and the reciprocal expression of these genes is likely fine-tuned to balance the respective costs. These findings highlight that the cost of clinical multiresistance plasmids can be largely due to particular resistance genes and their interaction with other cellular systems, while other resistance genes and the plasmid backbone can be cost-free. IMPORTANCE Multiresistance plasmids are one of the main drivers of antibiotic resistance development and spread. Their evolutionary success through the accumulation and mobilization of resistance genes is central to resistance evolution. In this study, we find that the cost of the introduction of a multiresistance plasmid was completely attributable to resistance genes, while the rest of the plasmid backbone is cost-free. The majority of resistance genes on the plasmid had no appreciable cost to the host cell even when overexpressed, indicating that plasmid-borne resistance can be cost-free. In contrast, the widespread genes bla_CTX-M-15 and tetAR were found to confer the whole cost of the plasmid by affecting specific cellular functions. These findings highlight how the evolution of resistance on plasmids is dependent on the amelioration of associated fitness costs and point at a conundrum regarding the high cost of some of the most widespread β-lactamase genes.

The evolution and role of eukaryotic-like domains in environmental intracellular bacteria: the battle with a eukaryotic cell

Intracellular pathogens that are able to thrive in different environments, such as Legionella spp. which preferentially live in protozoa in aquatic environments or environmental Chlamydiae which replicate either within protozoa or a range of animals, possess a plethora of cellular biology tools to influence their eukaryotic host. The host manipulation tools that evolved in the interaction with protozoa, confer these bacteria the capacity to also infect phylogenetically distinct eukaryotic cells, such as macrophages and thus they can also be human pathogens. To manipulate the host cell, bacteria use protein secretion systems and molecular effectors. Although these molecular effectors are encoded in bacteria, they are expressed and function in a eukaryotic context often mimicking or inhibiting eukaryotic proteins. Indeed, many of these effectors have eukaryotic-like domains. In this review we propose that the main pathways environmental intracellular bacteria need to subvert in order to establish the host eukaryotic cell as a replication niche are chromatin remodelling, ubiquitination signalling, and modulation of protein-protein interactions via tandem repeat domains. We then provide mechanistic insight into how these proteins might have evolved as molecular weapons. Finally, we highlight that in environmental intracellular bacteria the number of eukaryotic-like domains and proteins is considerably higher than in intracellular bacteria specialised to an isolated niche, such as obligate intracellular human pathogens. As mimics of eukaryotic proteins are critical components of host pathogen interactions, this distribution of eukaryotic-like domains suggests that the environment has selected them.

Ultra-deep long-read sequencing detects IS-mediated gene duplications as a potential trigger to generate arrays of resistance genes and a mechanism to induce novel gene variants such as blaCTX-M-243

BACKGROUND

Extended-spectrum β-lactamases (ESBLs) are enzymes that can render their hosts resistant to various β-lactam antibiotics. CTX-M-type enzymes are the most prevalent ESBLs and the main cause of resistance to third-generation cephalosporins in Enterobacteriaceae. The number of described CTX-M types is continuously rising, currently comprising over 240 variants. During routine screening we identified a novel blaCTX-M gene.

OBJECTIVES

To characterize a novel blaCTX-M variant harboured by a multidrug-resistant Escherichia coli isolate of sequence type ST354.

METHODS

Antibiotic susceptibilities were determined using broth microdilution. Genome and plasmid sequences were reconstructed using short- and long-read sequencing. The novel blaCTX-M locus was analysed using long-read and Sanger sequencing. Plasmid polymorphisms were determined in silico on a single plasmid molecule level.

RESULTS

The novel blaCTX-M-243 allele was discovered alongside a nearly identical blaCTX-M-104-containing gene array on a 219 kbp IncHI2A plasmid. CTX-M-243 differed from CTX-M-104 by only one amino acid substitution (N109S). Ultra-deep (2300-fold coverage) long-read sequencing revealed dynamic scaling of the blaCTX-M genetic contexts from one to five copies. Further antibiotic resistance genes such as blaTEM-1 also exhibited sequence heterogeneity but were stable in copy number.

CONCLUSIONS

We identified the novel ESBL gene blaCTX-M-243 and illustrate a dynamic system of varying blaCTX-M copy numbers. Our results highlight the constant emergence of new CTX-M family enzymes and demonstrate a potential evolutionary platform to generate novel ESBL variants and possibly other antibiotic resistance genes.

Evolutionary repair reveals an unexpected role of the tRNA modification m1G37 in aminoacylation

The tRNA modification m¹G37, introduced by the tRNA methyltransferase TrmD, is thought to be essential for growth in bacteria because it suppresses translational frameshift errors at proline codons. However, because bacteria can tolerate high levels of mistranslation, it is unclear why loss of m¹G37 is not tolerated. Here, we addressed this question through experimental evolution of trmD mutant strains of Escherichia coli. Surprisingly, trmD mutant strains were viable even if the m¹G37 modification was completely abolished, and showed rapid recovery of growth rate, mainly via duplication or mutation of the proline-tRNA ligase gene proS. Growth assays and in vitro aminoacylation assays showed that G37-unmodified tRNA^Pro is aminoacylated less efficiently than m¹G37-modified tRNA^Pro, and that growth of trmD mutant strains can be largely restored by single mutations in proS that restore aminoacylation of G37-unmodified tRNA^Pro. These results show that inefficient aminoacylation of tRNA^Pro is the main reason for growth defects observed in trmD mutant strains and that proS may act as a gatekeeper of translational accuracy, preventing the use of error-prone unmodified tRNA^Pro in translation. Our work shows the utility of experimental evolution for uncovering the hidden functions of essential genes and has implications for the development of antibiotics targeting TrmD.

High prevalence of unstable antibiotic heteroresistance in cyanobacteria causes resistance underestimation

Both cyanobacterial bloom and antibiotic resistance have aggravated worldwide and posed a global threat to public health in recent years. Cyanobacteria can exhibit discrepancy between their resistance genotype and susceptible phenotype due to antibiotic heteroresistance, which leads to difficulties in unambiguously classifying cyanobacterial strains as susceptible or resistant. Here we profiled the prevalence and mechanisms of antibiotic heteroresistance in cyanobacterial strains isolated from 50 sites across four eutrophicated lakes in China. Among 300 cyanobacterial isolates tested against 19 different antibiotics, over 90% of cyanobacterial isolates exhibited HR to multiple antibiotics and 19.5% of isolate/antibiotic interactions classified as susceptible by traditional minimum inhibitory concentration (MIC) estimates were designated heteroresistant. Over 97% of these monoclonal HR cases were unstable, with an increased resistance of subpopulations due to amplification of known resistance genes with high fitness cost. Wide-type cyanobacterial isolates of Synechococcus, Synechocystis, Anabaena and Microcystis aeruginosa exposed to sub-MIC level of four antibiotics evolved high-level resistance with little fitness cost, resulting in stable polyclonal HR. Both stable polyclonal HR and unstable monoclonal HR observed in different cyanobacterial strains can be promoted under environmental levels of antibiotic pressure. The highly prevalent and unstable monoclonal HR with the potential for susceptibility misclassification highlighted underestimation of cyanobacteria-derived antibiotic resistance. Cost-effective strategies should be developed to identify heteroresistance in cyanobacteria and to avoid false positive or negative results in traditional susceptibility testing.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

The highly dynamic nature of bacterial heteroresistance impairs its clinical detection

Many bacterial species and antibiotic classes exhibit heteroresistance, a phenomenon in which a susceptible bacterial isolate harbors a resistant subpopulation that can grow in the presence of an antibiotic and cause treatment failure. The resistant phenotype is often unstable and without antibiotic selection it reverts back to susceptibility. Here we studied the dynamics by which these resistant subpopulations are enriched in the presence of antibiotic and recede back to their baseline frequency in the absence of selection. An increasing understanding of this instability will allow more effective diagnostics and treatment of infections caused by heteroresistant bacteria. We show for clinical isolates of Escherichia coli and Salmonella enterica that different antibiotics at levels below the MIC of the susceptible main population can cause rapid enrichment of resistant subpopulations with increased copy number of genes that cause resistance. Modelling and growth rate measurements of bacteria with increased gene copy number in cultures and by microscopy of single-cells in a microfluidic chip show that the fitness cost of gene amplifications and their intrinsic instability drives their rapid loss in the absence of selection. Using a common antibiotic susceptibility test, we demonstrate that this test strongly underestimates the occurrence of heteroresistance in clinical isolates.

Bacterial populations can show heteroresistance, where an antibiotic resistant subpopulation is part of a susceptible one. Pereira et al. show in Escherichia coli and Salmonella enterica that disk diffusion, a common antibiotic susceptibility test, underestimates the occurrence of heteroresistance in clinical isolates.

Discerning the role of polymicrobial biofilms in the ascent, prevalence, and extent of heteroresistance in clinical practice

Antimicrobial therapy is facing a worrisome and underappreciated challenge, the phenomenon of heteroresistance (HR). HR has been gradually documented in clinically relevant pathogens (e.g. Pseudomonas aeruginosa, Staphylococcus aureus, Burkholderia spp., Acinetobacter baumannii, Klebsiella pneumoniae, Candida spp.) towards several drugs and is believed to complicate the clinical picture of chronic infections. This type of infections are typically mediated by polymicrobial biofilms, wherein microorganisms inherently display a wide range of physiological states, distinct metabolic pathways, diverging refractory levels of stress responses, and a complex network of chemical signals exchange. This review aims to provide an overview on the relevance, prevalence, and implications of HR in clinical settings. Firstly, related terminologies (e.g. resistance, tolerance, persistence), sometimes misunderstood and overlapped, were clarified. Factors generating misleading HR definitions were also uncovered. Secondly, the recent HR incidences reported in clinically relevant pathogens towards different antimicrobials were annotated. The potential mechanisms underlying such occurrences were further elucidated. Finally, the link between HR and biofilms was discussed. The focus was to recognize the presence of heterogeneous levels of resistance within most biofilms, as well as the relevance of polymicrobial biofilms in chronic infectious diseases and their role in resistance spreading. These topics were subject of a critical appraisal, gaining insights into the ascending clinical implications of HR in antimicrobial resistance spreading, which could ultimately help designing effective therapeutic options.

IS26-mediated amplification of blaOXA-1 and blaCTX-M-15 with concurrent outer membrane porin disruption associated with de novo carbapenem resistance in a recurrent bacteraemia cohort

BACKGROUND

Approximately half of clinical carbapenem-resistant Enterobacterales (CRE) isolates lack carbapenem-hydrolysing enzymes and develop carbapenem resistance through alternative mechanisms.

OBJECTIVES

To elucidate development of carbapenem resistance mechanisms from clonal, recurrent ESBL-positive Enterobacterales (ESBL-E) bacteraemia isolates in a vulnerable patient population.

METHODS

This study investigated a cohort of ESBL-E bacteraemia cases in Houston, TX, USA. Oxford Nanopore Technologies long-read and Illumina short-read sequencing data were used for comparative genomic analysis. Serial passaging experiments were performed on a set of clinical ST131 Escherichia coli isolates to recapitulate in vivo observations. Quantitative PCR (qPCR) and qRT-PCR were used to determine copy number and transcript levels of β-lactamase genes, respectively.

RESULTS

Non-carbapenemase-producing CRE (non-CP-CRE) clinical isolates emerged from an ESBL-E background through a concurrence of primarily IS26-mediated amplifications of blaOXA-1 and blaCTX-M-1 group genes coupled with porin inactivation. The discrete, modular translocatable units (TUs) that carried and amplified β-lactamase genes mobilized intracellularly from a chromosomal, IS26-bound transposon and inserted within porin genes, thereby increasing β-lactamase gene copy number and inactivating porins concurrently. The carbapenem resistance phenotype and TU-mediated β-lactamase gene amplification were recapitulated by passaging a clinical ESBL-E isolate in the presence of ertapenem. Clinical non-CP-CRE isolates had stable carbapenem resistance phenotypes in the absence of ertapenem exposure.

CONCLUSIONS

These data demonstrate IS26-mediated mechanisms underlying β-lactamase gene amplification with concurrent outer membrane porin disruption driving emergence of clinical non-CP-CRE. Furthermore, these amplifications were stable in the absence of antimicrobial pressure. Long-read sequencing can be utilized to identify unique mobile genetic element mechanisms that drive antimicrobial resistance.

Evolution of new enzymes by gene duplication and divergence

Thousands of new metabolic and regulatory enzymes have evolved by gene duplication and divergence since the dawn of life. New enzyme activities often originate from promiscuous secondary activities that have become important for fitness due to a change in the environment or a mutation. Mutations that make a promiscuous activity physiologically relevant can occur in the gene encoding the promiscuous enzyme itself, but can also occur elsewhere, resulting in increased expression of the enzyme or decreased competition between the native and novel substrates for the active site. If a newly useful activity is inefficient, gene duplication/amplification will set the stage for divergence of a new enzyme. Even a few mutations can increase the efficiency of a new activity by orders of magnitude. As efficiency increases, amplified gene arrays will shrink to provide two alleles, one encoding the original enzyme and one encoding the new enzyme. Ultimately, genomic rearrangements eliminate co-amplified genes and move newly evolved paralogs to a distant region of the genome.

Piperacillin/tazobactam resistance in a clinical isolate of Escherichia coli due to IS26-mediated amplification of blaTEM-1B

A phenotype of Escherichia coli and Klebsiella pneumoniae, resistant to piperacillin/tazobactam (TZP) but susceptible to carbapenems and 3rd generation cephalosporins, has emerged. The resistance mechanism associated with this phenotype has been identified as hyperproduction of the β-lactamase TEM. However, the mechanism of hyperproduction due to gene amplification is not well understood. Here, we report a mechanism of gene amplification due to a translocatable unit (TU) excising from an IS26-flanked pseudo-compound transposon, PTn6762, which harbours blaTEM-1B. The TU re-inserts into the chromosome adjacent to IS26 and forms a tandem array of TUs, which increases the copy number of blaTEM-1B, leading to TEM-1B hyperproduction and TZP resistance. Despite a significant increase in blaTEM-1B copy number, the TZP-resistant isolate does not incur a fitness cost compared to the TZP-susceptible ancestor. This mechanism of amplification of blaTEM-1B is an important consideration when using genomic data to predict susceptibility to TZP.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Biological amplification of low frequency mutations unravels laboratory culture history of the bio-threat agent Francisella tularensis

Challenges of investigating a suspected bio attack include establishing if microorganisms have been cultured to produce attack material and to identify their source. Addressing both issues, we have investigated genetic variations that emerge during laboratory culturing of the bacterial pathogen Francisella tularensis. Key aims were to identify genetic variations that are characteristic of laboratory culturing and explore the possibility of using biological amplification to identify genetic variation present at exceedingly low frequencies in a source sample. We used parallel serial passage experiments and high-throughput sequencing of F. tularensis to explore the genetic variation. We found that during early laboratory culture passages of F. tularensis, gene duplications emerged in the pathogen genome followed by single-nucleotide polymorphisms in genes for bacterial capsule synthesis. Based on a biological enrichment scheme and the use of high-throughput sequencing, we identified genetic variation that likely pre-existed in a source sample. The results support that capsule synthesis gene mutations are common during laboratory culture, and that a biological amplification strategy is useful for linking a F. tularensis sample to a specific laboratory variant among many highly similar variants.

How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation

Promiscuity is the coincidental ability of an enzyme to catalyze its native reaction and additional reactions that are not biological functions in the same active site. Promiscuity plays a central role in enzyme evolution and is thus a useful property for protein and metabolic engineering. This review examines enzyme evolution holistically, beginning with evaluating biochemical support for four enzyme evolution models. As expected, there is strong biochemical support for the subfunctionalization and innovation-amplification-divergence models, in which promiscuity is a central feature. In many cases, however, enzyme evolution is more complex than the models indicate, suggesting much is yet to be learned about selective pressures on enzyme function. A complete understanding of enzyme evolution must also explain the ability of metabolic networks to integrate new enzyme activities. Hidden within metabolic networks are underground metabolic pathways constructed from promiscuous activities. We discuss efforts to determine the diversity and pervasiveness of underground metabolism. Remarkably, several studies have discovered that some metabolic defects can be repaired via multiple underground routes. In prokaryotes, metabolic innovation is driven by connecting enzymes acquired by horizontal gene transfer (HGT) into the metabolic network. Thus, we end the review by discussing how the combination of promiscuity and HGT contribute to evolution of metabolism in prokaryotes. Future studies investigating the contribution of promiscuity to enzyme and metabolic evolution will need to integrate deeper probes into the influence of evolution on protein biophysics, enzymology, and metabolism with more complex and realistic evolutionary models. ENZYMES: lactate dehydrogenase (EC 1.1.1.27), malate dehydrogenase (EC 1.1.1.37), OSBS (EC 4.2.1.113), HisA (EC 5.3.1.16), TrpF, PriA (EC 5.3.1.24), R-mandelonitrile lyase (EC 4.1.2.10), Maleylacetate reductase (EC 1.3.1.32).

Mutations that improve efficiency of a weak-link enzyme are rare compared to adaptive mutations elsewhere in the genome

New enzymes often evolve by gene amplification and divergence. Previous experimental studies have followed the evolutionary trajectory of an amplified gene, but have not considered mutations elsewhere in the genome when fitness is limited by an evolving gene. We have evolved a strain of Escherichia coli in which a secondary promiscuous activity has been recruited to serve an essential function. The gene encoding the ‘weak-link’ enzyme amplified in all eight populations, but mutations improving the newly needed activity occurred in only one. Most adaptive mutations occurred elsewhere in the genome. Some mutations increase expression of the enzyme upstream of the weak-link enzyme, pushing material through the dysfunctional metabolic pathway. Others enhance production of a co-substrate for a downstream enzyme, thereby pulling material through the pathway. Most of these latter mutations are detrimental in wild-type E. coli, and thus would require reversion or compensation once a sufficient new activity has evolved.

Resistance to piperacillin/tazobactam in Escherichia coli resulting from extensive IS26-associated gene amplification of blaTEM-1

Background

bla TEM-1 encodes a narrow-spectrum β-lactamase that is inhibited by β-lactamase inhibitors and commonly present in Escherichia coli. Hyperproduction of blaTEM-1 may cause resistance to penicillin/β-lactamase inhibitor (P/BLI) combinations.

Objectives

To characterize EC78, an E. coli bloodstream isolate, resistant to P/BLI combinations, which contains extensive amplification of blaTEM-1 within the chromosome.

Methods

EC78 was sequenced using Illumina and Oxford Nanopore Technology (ONT) methodology. Configuration of blaTEM-1 amplification was probed using PCR. Expression of blaTEM-1 mRNA was determined using quantitative PCR and β-lactamase activity was determined spectrophotometrically in a nitrocefin conversion assay. Growth rate was assessed to determine fitness and stability of the gene amplification was assessed by passage in the absence of antibiotics.

Results

Illumina sequencing of EC78 identified blaTEM-1B as the only acquired β-lactamase preceded by the WT P3 promoter and present at a copy number of 182.6 with blaTEM-1B bracketed by IS26 elements. The chromosomal location of the IS26-blaTEM-1B amplification was confirmed by ONT sequencing. Hyperproduction of blaTEM-1 was confirmed by increased transcription of blaTEM-1 and β-lactamase activity and associated with a significant fitness cost; however, the array was maintained at a relatively high copy number for 150 generations. PCR screening for blaTEM amplification of isolates resistant to P/BLI combinations identified an additional strain containing an IS26-associated amplification of a blaTEM gene.

Conclusions

IS26-associated amplification of blaTEM can cause resistance to P/BLI combinations. This adaptive mechanism of resistance may be overlooked if simple methods of genotypic prediction (e.g. gene presence/absence) are used to predict antimicrobial susceptibility from sequencing data.