Positive selection and compensatory adaptation interact to stabilize non-transmissible plasmids

Plasmid-bacteria associations in the clinical context

Antimicrobial resistance (AMR) is one of the most pressing global health problems, with plasmids playing a central role in its evolution and dissemination. Over the past decades, many studies have investigated the ecoevolutionary dynamics between plasmids and their bacterial hosts. However, what drives the epidemiological success of certain plasmid-bacterium associations remains unclear. In this opinion article, we review which factors influence these associations and underline that studying plasmid-host interactions of clinical relevance is critical for understanding the evolution and spread of AMR. We also highlight the increasing importance of integrating experimental research with bioinformatics and machine learning tools to study plasmid-bacteria dynamics. This combined approach will assist researchers to dissect the molecular mechanisms underlying successful plasmid-host associations and to design strategies to prevent and predict future high-risk associations.

Evolution of a Plasmid Regulatory Circuit Ameliorates Plasmid Fitness Cost

Plasmids promote adaptation of bacteria by facilitating horizontal transfer of diverse genes, notably those conferring antibiotic resistance. Some plasmids, like those of the incompatibility group IncP-1, are known to replicate and persist in a broad range of bacteria. We investigated a poorly understood exception, the IncP-1β plasmid pBP136 from a clinical Bordetella pertussis isolate, which quickly became extinct in laboratory Escherichia coli populations. Through experimental evolution, we found that the inactivation of a previously uncharacterized plasmid gene, upf31, drastically improved plasmid persistence in E. coli. The gene inactivation caused alterations in the plasmid regulatory system, including decreased transcription of the global plasmid regulators (korA, korB, and korC) and numerous genes in their regulons. This is consistent with our findings that Upf31 represses its own transcription. It also caused secondary transcriptional changes in many chromosomal genes. In silico analyses predicted that Upf31 interacts with the plasmid regulator KorB at its C-terminal dimerization domain (CTD). We showed experimentally that adding the CTD of upf31/pBP136 to the naturally truncated upf31 allele of the stable IncP-1β archetype R751 results in plasmid destabilization in E. coli. Moreover, mutagenesis showed that upf31 alleles encoded on nearly half of the sequenced IncP-1β plasmids also possess this destabilization phenotype. While Upf31 might be beneficial in many hosts, we show that in E. coli some alleles have harmful effects that can be rapidly alleviated with a single mutation. Thus, broad-host-range plasmid adaptation to new hosts can involve fine-tuning their transcriptional circuitry through evolutionary changes in a single gene.

Fitness and adaptive evolution of a Rhodococcus sp. harboring dioxin-catabolic plasmids

Catabolic plasmids are critical factors in the degradation of recalcitrant xenobiotics, such as dioxins. Understanding the persistence and evolution of native catabolic plasmids is pivotal for controlling their function in microbial remediation. Here, we track the fitness and evolution of Rhodococcus sp. strain p52 harboring dioxin-catabolic plasmids under nonselective conditions without contaminant. Growth curve analysis and competition experiments demonstrated that pDF01 imposed fitness costs, whereas pDF02 conferred fitness benefits. During stability tests, pDF01 tended to be lost from the population, while pDF02 maintained at least one copy in the cell until proliferation of the 400th generation. Genome-wide gene expression profiling combined with codon usage bias analysis revealed that the high expression of pDF01 genes involved in dibenzofuran catabolism and regulation caused metabolic burdens. In contrast, potential cooperation between the pDF02-encoded short-chain dehydrogenase/reductase family oxidoreductase and the redox cofactor mycofactocin, which synthetic genes are located on the chromosome, may explain the benefit of pDF02. The fitness cost imposed by pDF01 was alleviated during adaptive evolution and was associated with the transcriptional downregulation of the dibenzofuran degradation genes on pDF01, and the global regulation of genome-wide gene expression involving basic metabolism, transport, and signal transduction. This study broadens our understandings on the persistence and evolution of dioxin-catabolic mega-plasmids, thus paving the way for the bioremediation of recalcitrant xenobiotic pollution in the environment.

Protozoa-enhanced conjugation frequency alters the dissemination of soil antibiotic resistance

Protozoa, as primary predators of soil bacteria, represent an overlooked natural driver in the dissemination of antibiotic resistance genes (ARGs). However, the effects of protozoan predation on ARGs dissemination at the community level, along with the underlying mechanisms, remain unclear. Here we used fluorescence-activated cell sorting, qPCR, combined with metagenomics and reverse transcription quantitative PCR, to unveil how protozoa (Colpoda steinii and Acanthamoeba castellanii) influence the plasmid-mediated transfer of ARGs to soil microbial communities. Protozoan predation reduced the absolute abundance of plasmids but promoted the expression of conjugation-associated genes, leading to a 5-fold and 4.5-fold increase in conjugation frequency in the presence of C. steinii and A. castellanii, respectively. Excessive oxidative stress, increased membrane permeability, and the provoked SOS response closely associated with the increased conjugative transfer. Protozoan predation also altered the plasmid host range and selected for specific transconjugant taxa along with ARGs and virulence factors carried by transconjugant communities. This study underscores the role of protozoa in the plasmid-mediated conjugative transfer of ARGs, providing new insights into microbial mechanisms that drive the dissemination of environmental antibiotic resistance.

Abundance measurements reveal the balance between lysis and lysogeny in the human gut microbiome

The human gut contains diverse communities of bacteriophage, whose interactions with the broader microbiome and potential roles in human health are only beginning to be uncovered. Here, we combine multiple types of data to quantitatively estimate gut phage population dynamics and lifestyle characteristics in human subjects. Unifying results from previous studies, we show that an average human gut contains a low ratio of phage particles to bacterial cells (~1:100), but a much larger ratio of phage genomes to bacterial genomes (~4:1), implying that most gut phage are effectively temperate (e.g., integrated prophage, phage-plasmids, etc.). By integrating imaging and sequencing data with a generalized model of temperate phage dynamics, we estimate that phage induction and lysis occurs at a low average rate (~0.001-0.01 per bacterium per day), imposing only a modest fitness burden on their bacterial hosts. Consistent with these estimates, we find that the phage composition of a diverse synthetic community in gnotobiotic mice can be quantitatively predicted from bacterial abundances alone, while still exhibiting phage diversity comparable to native human microbiomes. These results provide a foundation for interpreting existing and future studies on links between the gut virome and human health.

Decoding the enigma: unveiling the transmission characteristics of waterfowl-associated bla NDM-5-positive Escherichia coli in select regions of China

Escherichia coli (E. coli) serves as a critical indicator microorganism for assessing the prevalence and dissemination of antibiotic resistance, notably harboring various antibiotic-resistant genes (ARGs). Among these, the emergence of the bla NDM gene represents a significant threat to public health, especially since carbapenem antibiotics are vital for treating severe infections caused by Gram-negative bacteria. This study aimed to characterize the antibiotic resistance features of bla NDM-5-positive E. coli strains isolated from waterfowl in several regions of China and elucidate the dissemination patterns of the bla NDM-5 gene. We successfully isolated 103 bla NDM-5-positive E. coli strains from 431 intestinal fecal samples obtained from waterfowl across five provincial-level units in China, with all strains exhibiting multidrug resistance (MDR). Notably, the bla NDM-5 gene was identified on plasmids, which facilitate efficient and stable horizontal gene transfer (HGT). Our adaptability assays indicated that while the bla NDM-5-positive plasmid imposed a fitness cost on the host bacteria, the NDM-5 protein was successfully induced and purified, exhibiting significant enzymatic activity. One strain, designated DY51, exhibited a minimum inhibitory concentration (MIC) for imipenem of 4 mg/L, which escalated to 512 mg/L following exposure to increasing imipenem doses. This altered strain demonstrated stable resistance to imipenem alongside improved adaptability, correlating with elevated relative expression levels of the bla NDM-5 and overexpression of efflux pumps. Collectively, this study highlights the horizontal dissemination of the bla NDM-5 plasmid among E. coli strains, confirms the associated fitness costs, and provides insights into the mechanisms underlying the stable increase in antibiotic resistance to imipenem. These findings offer a theoretical framework for understanding the dissemination dynamics of bla NDM-5 in E. coli, which is essential for developing effective strategies to combat carbapenem antibiotic resistance.

Porin deficiency or plasmid copy number increase mediated carbapenem-resistant Escherichia coli resistance evolution

This study investigated resistance evolution mechanisms of conjugated plasmids and bacterial hosts under different concentrations of antibiotic pressure. Ancestral strain ECNX52 was constructed by introducing the blaNDM-5-carrying IncX3 plasmid into E. coli C600, and was subjected to laboratory evolution under different concentrations of meropenem pressure. Minimal inhibitory concentrations and conjugation frequency were determined. Fitness of these strains was assessed. Whole genome sequencing and transcriptional changes were performed. Ancestral host or plasmids were recombined with evolved hosts or plasmids to verify plasmid or host factors in resistance evolution. Role of the repA mutation on plasmid copy number was determined. Two out of the four clones (EM2N1 and EM2N3) exhibited four-fold increase in MIC when exposed to a continuous pressure of 2 μg/mL MEM (1/32 MIC), by down regulating expression of outer membrane protein ompF. Besides, all four clones displayed four-fold increase in MIC and higher conjugation frequency when subjected to a continuous pressure of 4 μg/mL MEM (1/16 MIC), attributing to increasing plasmid copy number generated by repA D140Y (GAT→TAT) mutation. Bacterial hosts and conjugative plasmids can undergo resistance evolution under certain concentrations of antimicrobial pressure by reducing the expression of outer membrane proteins or increasing plasmid copy numbers.

Evolution and maintenance of a large multidrug-resistant plasmid in a Salmonella enterica Typhimurium host under differing antibiotic selection pressures

The dissemination of antibiotic resistance genes (ARGs) through plasmids is a major mechanism for the development of bacterial antimicrobial resistance. The adaptation and evolution mechanisms of multidrug-resistant (MDR) plasmids with their hosts are not fully understood. Herein, we conducted experimental evolution of a 244 kb MDR plasmid (pJXP9) under various conditions including no antibiotics and mono- or combinational drug treatments of colistin (CS), cefotaxime (CTX), and ciprofloxacin (CIP). Our results showed that long-term with or without positive selections for pJXP9, spanning approximately 600 generations, led to modifications of the plasmid-encoded MDR and conjugative transfer regions. These modifications could mitigate the fitness cost of plasmid carriage and enhance plasmid maintenance. The extent of plasmid modifications and the evolution of plasmid-encoded antibiotic resistance depended on treatment type, particularly the drug class and duration of exposure. Interestingly, prolonged exposure to mono- and combinational drugs of CS and CIP resulted in a substantial loss of the plasmid-encoded MDR region and antibiotic resistance, comparable to the selection condition without antibiotic. By contrast, combinational treatment with CTX contributed to the maintenance of the MDR region over a long period of time. Furthermore, drug selection was able to maintain and even amplify the corresponding plasmid-encoded ARGs, with co-selection of ARGs in the adjacent regions. In addition, parallel mutations in chromosomal arcA were also found to be associated with pJXP9 plasmid carriage among endpoint-evolved clones from diverse treatments. Meanwhile, arcA deletion improved the persistence of pJXP9 plasmid without drugs. Overall, our findings indicated that plasmid-borne MDR region deletion and chromosomal arcA inactivation mutation jointly contributed to co-adaptation and co-evolution between MDR IncHI2 plasmid and Salmonella Typhimurium under different drug selection pressure.IMPORTANCEThe plasmid-mediated dissemination of antibiotic resistance genes has become a significant concern for human health, even though the carriage of multidrug-resistant (MDR) plasmids is frequently associated with fitness costs for the bacterial host. However, the mechanisms by which MDR plasmids and bacterial pairs evolve plasmid-mediated antibiotic resistance in the presence of antibiotic selections are not fully understood. Herein, we conducted an experimental evolution of a large multidrug-resistant plasmid in a Salmonella enterica Typhimurium host under single and combinatorial drug selection pressures. Our results show the adaptive evolution of plasmid-encoded antibiotic resistance through alterations of the MDR region in the plasmid, in particular substantial loss of the MDR region, in response to different positive selections, especially mono- and combinational drugs of colistin and ciprofloxacin. In addition, strong parallel mutations in chromosomal arcA were associated with pJXP9 carriage in Salmonella Typhimurium from diverse treatments. Our results thus highlight promoting the loss of the plasmid's MDR region could offer an alternative approach for combating plasmid-encoded antibiotic resistance.

The evolutionary landscape of prokaryotic chromosome/plasmid balance

The balance between chromosomal and plasmid DNAs determines the genomic plasticity of prokaryotes. Natural selections, acting on the level of organisms or plasmids, shape the abundances of plasmid DNAs in prokaryotic genomes. Despite the importance of plasmids in health and engineering, there have been rare systematic attempts to quantitatively model and predict the determinants underlying the strength of different selection forces. Here, we develop a metabolic flux model that describes the intracellular resource competition between chromosomal and plasmid-encoded reactions. By coarse graining, this model predicts a landscape of natural selections on chromosome/plasmid balance, which is featured by the tradeoff between phenotypic and non-phenotypic selection pressures. This landscape is further validated by the observed pattern of plasmid distributions in the vast collection of prokaryotic genomes retrieved from the NCBI database. Our results establish a universal paradigm to understand the prokaryotic chromosome/plasmid interplay and provide insights into the evolutionary origin of plasmid diversity.

Why do mobile genetic elements transfer DNA of their hosts?

The prokaryote world is replete with mobile genetic elements (MGEs) - self-replicating entities that can move within and between their hosts. Many MGEs not only transfer their own DNA to new hosts but also transfer host DNA located elsewhere on the chromosome in the process. This could potentially lead to indirect benefits to the host when the resulting increase in chromosomal variation results in more efficient natural selection. We review the diverse ways in which MGEs promote the transfer of host DNA and explore the benefits and costs to MGEs and hosts. In many cases, MGE-mediated transfer of host DNA might not be selected for because of a sex function, but evidence of MGE domestication suggests that there may be host benefits of MGE-mediated sex.

Insights into the ecology of the infant gut plasmidome

Plasmids are small DNA molecules that enable bacteria to share beneficial traits, influencing microbial communities. However, their role within the human gut microbiome remains largely unknown. In this study, we investigate the gut microbiomes of 34 mother-child cohorts, employing a plasmid analysis workflow to understand the impact of plasmids on the gut microbiome. We create a plasmid phylogenetic tree, devise a method for assigning plasmid hosts, and examine potential plasmid transfer networks. Our research discovers a wide variety of previously unidentified plasmid sequences, indicating that current databases do not fully represent the gut plasmidome. Interestingly, infants display greater plasmid diversity compared to mothers and other healthy adults. We find that Bacteroidota, a major bacterial phylum, serves as the primary host for gut plasmids and plays a dominant role in gut plasmid transfer events. Additionally, plasmids broaden the genetic capabilities of bacteria, with their influence on bacterial function becoming more apparent as children's gut microbiomes develop. This study sheds light on the role of plasmids in the infant gut microbiome, making a significant contribution to our understanding of plasmid biology.

Compensatory evolution of chromosomes and plasmids counteracts the plasmid fitness cost

Plasmids incur a fitness cost that has the potential to restrict the dissemination of resistance in bacterial pathogens. However, bacteria can overcome this disadvantage by compensatory evolution to maintain their resistance. Compensatory evolution can occur via both chromosomes and plasmids, but there are a few reviews regarding this topic, and most of them focus on plasmids. In this review, we provide a comprehensive overview of the currently reported mechanisms underlying compensatory evolution on chromosomes and plasmids, elucidate key targets regulating plasmid fitness cost, and discuss future challenges in this field. We found that compensatory evolution on chromosomes primarily arises from mutations in transcriptional regulatory factors, whereas compensatory evolution of plasmids predominantly involves three pathways: plasmid copy number regulation, conjugation transfer efficiency, and expression of antimicrobial resistance (AMR) genes. Furthermore, the importance of reasonable selection of research subjects and effective integration of diverse advanced research methods is also emphasized in our future study on compensatory mechanisms. Overall, this review establishes a theoretical framework that aims to provide innovative ideas for minimizing the emergence and spread of AMR genes.

Hijackers, hitchhikers, or co-drivers? The mysteries of mobilizable genetic elements

Mobile genetic elements shape microbial gene repertoires and populations. Recent results reveal that many, possibly most, microbial mobile genetic elements require helpers to transfer between genomes, which we refer to as Hitcher Genetic Elements (hitchers or HGEs). They may be a large fraction of pathogenicity and resistance genomic islands, whose mechanisms of transfer have remained enigmatic for decades. Together with their helper elements and their bacterial hosts, hitchers form tripartite networks of interactions that evolve rapidly within a parasitism-mutualism continuum. In this emerging view of microbial genomes as communities of mobile genetic elements many questions arise. Which elements are being moved, by whom, and how? How often are hitchers costly hyper-parasites or beneficial mutualists? What is the evolutionary origin of hitchers? Are there key advantages associated with hitchers' lifestyle that justify their unexpected abundance? And why are hitchers systematically smaller than their helpers? In this essay, we start answering these questions and point ways ahead for understanding the principles, origin, mechanisms, and impact of hitchers in bacterial ecology and evolution.

A Newly Incompatibility F Replicon Allele (FIB81) in Extensively Drug-Resistant Escherichia coli Isolated from Diseased Broilers

Multiple drug resistance (MDR) has gained pronounced attention among Enterobacterales. The transfer of multiple antimicrobial resistance genes, frequently carried on conjugative incompatibility F (IncF) plasmids and facilitating interspecies resistance transmission, has been linked to Salmonella spp. and E. coli in broilers. In Egypt, the growing resistance is exacerbated by the limited clinical efficacy of many antimicrobials. In this study, IncF groups were screened and characterized in drug-resistant Salmonella spp. and E. coli isolated from broilers. The antimicrobial resistance profile, PCR-based replicon typing of bacterial isolates pre- and post-plasmid curing, and IncF replicon allele sequence typing were investigated. Five isolates of E. coli (5/31; 16.13%) and Salmonella spp. (5/36; 13.89%) were pan-susceptible to the examined antimicrobial agents, and 85.07% of tested isolates were MDR and extensively drug-resistant (XDR). Twelve MDR and XDR E. coli and Salmonella spp. isolates were examined for the existence of IncF replicons (FII, FIA, and FIB). They shared resistance to ampicillin, ampicillin/sulbactam, amoxicillin/clavulanate, doxycycline, cefotaxime, and colistin. All isolates carried from one to two IncF replicons. The FII-FIA-FIB+ and FII-FIA+FIB- were the predominant replicon patterns. FIB was the most frequently detected replicon after plasmid curing. Three XDR E. coli isolates that were resistant to 12-14 antimicrobials carried a newly FIB replicon allele with four nucleotide substitutions: C99→A, G112→T, C113→T, and G114→A. These findings suggest that broilers are a significant reservoir of IncF replicons with highly divergent IncF-FIB plasmid incompatibility groups circulating among XDR Enterobacterales. Supporting these data with additional comprehensive epidemiological studies involving replicons other than the IncF can provide insights for implementing efficient policies to prevent the spreading of new replicons to humans.

Fitness costs of Tn1546-type transposons harboring the vanA operon by plasmid type and structural diversity in Enterococcus faecium

BACKGROUND

This study analyzed the genetic traits and fitness costs of vancomycin-resistant Enterococcus faecium (VREfm) blood isolates carrying Tn1546-type transposons harboring the vanA operon.

METHODS

All E. faecium blood isolates were collected from eight general hospitals in South Korea during one-year study period. Antimicrobial susceptibility testing and vanA and vanB PCR were performed. Growth rates of E. faecium isolates were determined. The vanA-positive isolates were subjected to whole genome sequencing and conjugation experiments.

RESULTS

Among 308 E. faecium isolates, 132 (42.9%) were positive for vanA. All Tn1546-type transposons harboring the vanA operon located on the plasmids, but on the chromosome in seven isolates. The plasmids harboring the vanA operon were grouped into four types; two types of circular, nonconjugative plasmids (Type A, n = 50; Type B, n = 46), and two types of putative linear, conjugative plasmids (Type C, n = 16; Type D, n = 5). Growth rates of vanA-positive E. faecium isolates were significantly lower than those of vanA-negative isolates (P < 0.001), and reduction in growth rate under vancomycin pressure was significantly larger in isolates harboring putative linear plasmids than in those harboring circular plasmids (P = 0.020).

CONCLUSIONS

The possession of vanA operon was costly to bacterial hosts in antimicrobial-free environment, which provide evidence for the importance of reducing vancomycin pressure for prevention of VREfm dissemination. Fitness burden to bacterial hosts was varied by type and size of the vanA operon-harboring plasmid.

Tebuconazole promotes spread of a multidrug-resistant plasmid into soil bacteria to form new resistant bacterial strains

The development of antibiotic resistance threatens human and environmental health. Non-antibiotic stressors, including fungicides, may contribute to the spread of antibiotic resistance genes (ARGs). We determined the promoting effects of tebuconazole on ARG dissemination using a donor, Escherichia coli MG1655, containing a multidrug-resistant fluorescent plasmid (RP4) and a recipient (E. coli HB101). The donor was then incorporated into the soil to test whether tebuconazole could accelerate the spread of RP4 into indigenous bacteria. Tebuconazole promoted the transfer of the RP4 plasmid from the donor into the recipient via overproduction of reactive oxygen species (ROS), enhancement of cell membrane permeability and regulation of related genes. The dissemination of the RP4 plasmid from the donor to soil bacteria was significantly enhanced by tebuconazole. RP4 plasmid could be propagated into more genera of bacteria in tebuconazole-contaminated soil as the exposure time increased. These findings demonstrate that the fungicide tebuconazole promotes the spread of the RP4 plasmid into indigenous soil bacteria, revealing the potential risk of tebuconazole residues enhancing the dissemination of ARGs in soil environments.

Coexistence of a nonresistance-conferring IncI1 plasmid favors persistence of the blaCTX-M-bearing IncFII plasmid in Escherichia coli

The interaction between coexisting plasmids can affect plasmid-carried resistance gene persistence and spread. However, whether the persistence of the blaCTX-M gene in clinical Enterobacteriaceae is related to the interaction of coresident nonresistance-conferring plasmids has not been reported. This study was initiated to elucidate how a nonresistance-conferring IncI1 plasmid affected the blaCTX-M-bearing IncFII plasmid colocated on the same cell. Herein, we constructed three isogenic derivatives of E. coli C600, designated as C600FII, C600I1, and C600FII+I1, which harbored the blaCTX-M-IncFII plasmid and/or the nonresistance-IncI1 one. We discovered that strain C600FII+I1 conferred higher fitness advantages than strain C600FII; also, the stability of the blaCTX-M-IncFII plasmid was noticeably improved in an antibiotic-free environment when it coexisted with the IncI1 plasmid. To further explore why the IncI1 plasmid enhanced the persistence of the blaCTX-M-IncFII plasmid, we assessed the blaCTX-M-IncFII plasmid's copy numbers, conjugation frequencies, and rep gene expressions in strains C600FII and C600FII+I1. The results demonstrated that the rep expressions of the blaCTX-M-IncFII plasmid in strain C600FII+I1 was greatly decreased, along with the plasmid's copy numbers and mating efficiencies, compared to those in strain C600FII. Moreover, further study revealed that the intracellular ATP levels of strain C600FII+I1 were far lower than those of strain C600FII. Our findings confirmed that coexistence of the nonresistance-IncI1 plasmid can keep the blaCTX-M-IncFII plasmid more stable by increasing the fitness advantages of the host bacteria, which will pose a threat to preventing the long-term presence of the plasmid-carried blaCTX-M gene in clinical Enterobacteriaceae.

IMPORTANCE

So far, plasmid-carried blaCTX-M is still the most common extended-spectrum beta-lactamase (ESBL) genotype in clinical settings worldwide. Except for the widespread use of third-generation cephalosporins, the interaction between coexisting plasmids can also affect the long-term stable existence of the blaCTX-M gene; however, the study on that is still sparse. In the present study, we assess the interaction of coinhabitant plasmids blaCTX-M-IncFII and nonresistance-IncI1. Our results confirmed that the increased fitness advantages of strain C600FII+I1 were attributable to the cohabitant nonresistance-IncI1 plasmid, which largely reduced the intracellular ATP levels of host bacteria, thus decreasing the rep gene expression of the blaCTX-M-IncFII plasmid, its copy numbers, and mating efficiencies, while the higher fitness advantages of strain C600FII+I1 enhanced the persistence of the blaCTX-M-IncFII plasmid. The results indicate that the nonresistance-IncI1 plasmid contributes to the long-term existence of the blaCTX-M-IncFII plasmid, implying a potentially new strategy for controlling the spread of resistance plasmids in clinical settings by targeting nonresistance plasmids.

Costs of antibiotic resistance genes depend on host strain and environment and can influence community composition

Antibiotic resistance genes (ARGs) benefit host bacteria in environments containing corresponding antibiotics, but it is less clear how they are maintained in environments where antibiotic selection is weak or sporadic. In particular, few studies have measured if the direct effect of ARGs on host fitness is fixed or if it depends on the host strain, perhaps marking some ARG-host combinations as selective refuges that can maintain ARGs in the absence of antibiotic selection. We quantified the fitness effects of six ARGs in 11 diverse Escherichia spp. strains. Three ARGs (blaTEM-116, cat and dfrA5, encoding resistance to β-lactams, chloramphenicol, and trimethoprim, respectively) imposed an overall cost, but all ARGs had an effect in at least one host strain, reflecting a significant strain interaction effect. A simulation predicts these interactions can cause the success of ARGs to depend on available host strains, and, to a lesser extent, can cause host strain success to depend on the ARGs present in a community. These results indicate the importance of considering ARG effects across different host strains, and especially the potential of refuge strains to allow resistance to persist in the absence of direct selection, in efforts to understand resistance dynamics.

Evolution of a Plasmid Regulatory Circuit Ameliorates Plasmid Fitness Cost

Plasmids play a major role in rapid adaptation of bacteria by facilitating horizontal transfer of diverse genes, most notably those conferring antibiotic resistance. While most plasmids that replicate in a broad range of bacteria also persist well in diverse hosts, there are exceptions that are poorly understood. We investigated why a broad-host range plasmid, pBP136, originally found in clinical Bordetella pertussis isolates, quickly became extinct in laboratory Escherichia coli populations. Through experimental evolution we found that inactivation of a previously uncharacterized plasmid gene, upf31, drastically improved plasmid maintenance in E. coli. This gene inactivation resulted in decreased transcription of the global plasmid regulators (korA, korB, and korC) and numerous genes in their regulons. It also caused transcriptional changes in many chromosomal genes primarily related to metabolism. In silico analyses suggested that the change in plasmid transcriptome may be initiated by Upf31 interacting with the plasmid regulator KorB. Expression of upf31 in trans negatively affected persistence of pBP136Δupf31 as well as the closely related archetypal IncP-1β plasmid R751, which is stable in E. coli and natively encodes a truncated upf31 allele. Our results demonstrate that while the upf31 allele in pBP136 might advantageously modulate gene expression in its original host, B. pertussis, it has harmful effects in E. coli. Thus, evolution of a single plasmid gene can change the range of hosts in which that plasmid persists, due to effects on the regulation of plasmid gene transcription.

Prophage maintenance is determined by environment-dependent selective sweeps rather than mutational availability

Prophages, viral sequences integrated into bacterial genomes, can be beneficial and costly. Despite the risk of prophage activation and subsequent bacterial death, active prophages are present in most bacterial genomes. However, our understanding of the selective forces that maintain prophages in bacterial populations is limited. Combining experimental evolution with stochastic modeling, we show that prophage maintenance and loss are primarily determined by environmental conditions that alter the net fitness effect of a prophage on its bacterial host. When prophages are too costly, they are rapidly lost through environment-specific sequences of selective sweeps. Conflicting selection pressures that select against the prophage but for a prophage-encoded accessory gene can maintain prophages. The dynamics of prophage maintenance additionally depend on the sociality of this accessory gene. Prophage-encoded genes that exclusively benefit the lysogen maintain prophages at higher frequencies compared with genes that benefit the entire population. That is because the latter can protect phage-free "cheaters," reducing the benefit of maintaining the prophage. Our simulations suggest that environmental variation plays a larger role than mutation rates in determining prophage maintenance. These findings highlight the complexity of selection pressures that act on mobile genetic elements and challenge our understanding of the role of environmental factors relative to random chance events in shaping the evolutionary trajectory of bacterial populations. By shedding light on the key factors that shape microbial populations in the face of environmental changes, our study significantly advances our understanding of the complex dynamics of microbial evolution and diversification.

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.

Manipulating multi-level selection in a fungal entomopathogen reveals social conflicts and a method for improving biocontrol traits

Changes in parasite virulence are commonly expected to lead to trade-offs in other life history traits that can affect fitness. Understanding these trade-offs is particularly important if we want to manipulate the virulence of microbial biological control agents. Theoretically, selection across different spatial scales, i.e. between- and within-hosts, shapes these trade-offs. However, trade-offs are also dependent on parasite biology. Despite their applied importance the evolution of virulence in fungal parasites is poorly understood: virulence can be unstable in culture and commonly fails to increase in simple passage experiments. We hypothesized that manipulating selection intensity at different scales would reveal virulence trade-offs in a fungal pathogen of aphids, Akanthomyces muscarius. Starting with a genetically diverse stock we selected for speed of kill, parasite yield or infectivity by manipulating competition within and between hosts and between-populations of hosts over 7 rounds of infection. We characterized ancestral and evolved lineages by whole genome sequencing and by measuring virulence, growth rate, sporulation and fitness. While several lineages showed increases in virulence, we saw none of the trade-offs commonly found in obligately-killing parasites. Phenotypically similar lineages within treatments often shared multiple single-nucleotide variants, indicating strong convergent evolution. The most dramatic phenotypic changes were in timing of sporulation and spore production in vitro. We found that early sporulation led to reduced competitive fitness but could increase yield of spores on media, a trade-off characteristic of social conflict. Notably, the selection regime with strongest between-population competition and lowest genetic diversity produced the most consistent shift to early sporulation, as predicted by social evolution theory. Multi-level selection therefore revealed social interactions novel to fungi and showed that these biocontrol agents have the genomic flexibility to improve multiple traits-virulence and spore production-that are often in conflict in other parasites.

Horizontal gene transfer facilitates the molecular reverse-evolution of antibiotic sensitivity in experimental populations of H. pylori

Reversing the evolution of traits harmful to humans, such as antimicrobial resistance, is a key ambition of applied evolutionary biology. A major impediment to reverse evolution is the relatively low spontaneous mutation rates that revert evolved genotypes back to their ancestral state. However, the repeated re-introduction of ancestral alleles by horizontal gene transfer (HGT) could make reverse evolution likely. Here we evolve populations of an antibiotic-resistant strain of Helicobacter pylori in growth conditions without antibiotics while introducing an ancestral antibiotic-sensitive allele by HGT. We evaluate reverse evolution using DNA sequencing and find that HGT facilitates the molecular reverse evolution of the antibiotic resistance allele, and that selection for high rates of HGT drives the evolution of increased HGT rates in low-HGT treatment populations. Finally, we use a theoretical model and carry out simulations to infer how the fitness costs of antibiotic resistance, rates of HGT and effects of genetic drift interact to determine the probability and predictability of reverse evolution.

Adaptive evolution of plasmid and chromosome contributes to the fitness of a blaNDM-bearing cointegrate plasmid in Escherichia coli

Large cointegrate plasmids recruit genetic features of their parental plasmids and serve as important vectors in the spread of antibiotic resistance. They are now frequently found in clinical settings, raising the issue of how to limit their further transmission. Here, we conducted evolutionary research of a large blaNDM-positive cointegrate within Escherichia coli C600, and discovered that adaptive evolution of chromosome and plasmid jointly improved bacterial fitness, which was manifested as enhanced survival ability for in vivo and in vitro pairwise competition, biofilm formation, and gut colonization ability. From the plasmid aspect, large-scale DNA fragment loss is observed in an evolved clone. Although the evolved plasmid imposes a negligible fitness cost on host bacteria, its conjugation frequency is greatly reduced, and the deficiency of anti-SOS gene psiB is found responsible for the impaired horizontal transferability rather than the reduced fitness cost. These findings unveil an evolutionary strategy in which the plasmid horizontal transferability and fitness cost are balanced. From the chromosome perspective, all evolved clones exhibit parallel mutations in the transcriptional regulatory stringent starvation Protein A gene sspA. Through a sspA knockout mutant, transcriptome analysis, in vitro transcriptional activity assay, RT-qPCR, motility test, and scanning electron microscopy techniques, we demonstrated that the mutation in sspA reduces its transcriptional inhibitory capacity, thereby improving bacterial fitness, biofilm formation ability, and gut colonization ability by promoting bacterial flagella synthesis. These findings expand our knowledge of how cointegrate plasmids adapt to new bacterial hosts.

Structure and dispersion of the conjugative mobilome in surface ocean bacterioplankton

Mobile genetic elements (MGEs), collectively referred to as the "mobilome", can have a significant impact on the fitness of microbial communities and therefore on ecological processes. Marine MGEs have mainly been associated with wide geographical and phylogenetic dispersal of adaptative traits. However, whether the structure of this mobilome exhibits deterministic patterns in the natural community is still an open question. The aim of this study was to characterize the structure of the conjugative mobilome in the ocean surface bacterioplankton by searching the publicly available marine metagenomes from the TARA Oceans survey, together with molecular markers, such as relaxases and type IV coupling proteins of the type IV secretion system (T4SS). The T4SS machinery was retrieved in more abundance than relaxases in the surface marine bacterioplankton. Moreover, among the identified MGEs, mobilizable elements were the most abundant, outnumbering self-conjugative sequences. Detection of a high number of incomplete T4SSs provides insight into possible strategies related to trans-acting activity between MGEs, and accessory functions of the T4SS (e.g. protein secretion), allowing the host to maintain a lower metabolic burden in the highly dynamic marine system. Additionally, the results demonstrate a wide geographical dispersion of MGEs throughout oceanic regions, while the Southern Ocean appears segregated from other regions. The marine mobilome also showed a high similarity of functions present in known plasmid databases. Moreover, cargo genes were mostly related to DNA processing, but scarcely associated with antibiotic resistance. Finally, within the MGEs, integrative and conjugative elements showed wider marine geographic dispersion than plasmids.

Plasmids, a molecular cornerstone of antimicrobial resistance in the One Health era

Antimicrobial resistance (AMR) poses a substantial threat to human health. The widespread prevalence of AMR is, in part, due to the horizontal transfer of antibiotic resistance genes (ARGs), typically mediated by plasmids. Many of the plasmid-mediated resistance genes in pathogens originate from environmental, animal or human habitats. Despite evidence that plasmids mobilize ARGs between these habitats, we have a limited understanding of the ecological and evolutionary trajectories that facilitate the emergence of multidrug resistance (MDR) plasmids in clinical pathogens. One Health, a holistic framework, enables exploration of these knowledge gaps. In this Review, we provide an overview of how plasmids drive local and global AMR spread and link different habitats. We explore some of the emerging studies integrating an eco-evolutionary perspective, opening up a discussion about the factors that affect the ecology and evolution of plasmids in complex microbial communities. Specifically, we discuss how the emergence and persistence of MDR plasmids can be affected by varying selective conditions, spatial structure, environmental heterogeneity, temporal variation and coexistence with other members of the microbiome. These factors, along with others yet to be investigated, collectively determine the emergence and transfer of plasmid-mediated AMR within and between habitats at the local and global scale.

Antimicrobial resistance level and conjugation permissiveness shape plasmid distribution in clinical enterobacteria

Conjugative plasmids play a key role in the dissemination of antimicrobial resistance (AMR) genes across bacterial pathogens. AMR plasmids are widespread in clinical settings, but their distribution is not random, and certain associations between plasmids and bacterial clones are particularly successful. For example, the globally spread carbapenem resistance plasmid pOXA-48 can use a wide range of enterobacterial species as hosts, but it is usually associated with a small number of specific Klebsiella pneumoniae clones. These successful associations represent an important threat for hospitalized patients. However, knowledge remains limited about the factors determining AMR plasmid distribution in clinically relevant bacteria. Here, we combined in vitro and in vivo experimental approaches to analyze pOXA-48-associated AMR levels and conjugation dynamics in a collection of wild-type enterobacterial strains isolated from hospitalized patients. Our results revealed significant variability in these traits across different bacterial hosts, with Klebsiella spp. strains showing higher pOXA-48-mediated AMR and conjugation frequencies than Escherichia coli strains. Using experimentally determined parameters, we developed a simple mathematical model to interrogate the contribution of AMR levels and conjugation permissiveness to plasmid distribution in bacterial communities. The simulations revealed that a small subset of clones, combining high AMR levels and conjugation permissiveness, play a critical role in stabilizing the plasmid in different polyclonal microbial communities. These results help to explain the preferential association of plasmid pOXA-48 with K. pneumoniae clones in clinical settings. More generally, our study reveals that species- and strain-specific variability in plasmid-associated phenotypes shape AMR evolution in clinically relevant bacterial communities.

Decoding the enigma: unveiling the molecular transmission of avian-associated tet(X4)-positive E. coli in Sichuan Province, China

Tigecycline is considered one of the "last resort antibiotics" for treating complex infections caused by multidrug-resistant (MDR) bacteria, especially for combating clinical resistant strains that produce carbapenemases. However, the tet(X4) gene, which carried by different plasmids can mediate high levels of bacterial resistance to tigecycline, was first reported in 2019. Here, we report the emergence of the plasmid-mediated tet(X4) in avian environment of Sichuan Province. A total of 21 tet(X4)-positive Escherichia coli (E. coli) strains were isolated and identified from avian samples in selected regions, with an isolation rate of 1.6% (21/1,286), and all of them were MDR strains. Multilocus Sequence Typing (MLST) method was used to classify the 21 tet(X4)-positive E. coli into the ST206, ST761, ST155, ST1638, ST542, and ST767 types, which also belong to the 3 phylogenetic subgroups A, B1, and C. Tet(X4) is located on mobile plasmids that can be efficiently and stably propagated. The results of fitness cost experiments showed that tet(X4)-positive plasmids may incur some fitness cost to host bacteria, but different tet(X4)-positive plasmids bring about differential fitness costs. Whole-genome sequencing further confirmed the tet(X4) gene can be located on IncX1-type plasmids and the core genetic structures are ISVsa3-rdmc-tet(X4) or rdmc-tet(X4)-ISVsa3, the former is a 7 copies tandem repeat structure. In this study, we isolated and identified tet(X4)-positive E. coli from the avian origin in Sichuan, analyzed the mobility of the tet(X4) by conjugational transfer and S1-PFGE, and evaluated the biological characteristics of the tet(X4)-positive plasmid using the results of conjugational frequency, plasmid stability, and fitness costs. Finally, combined with the third-generation whole-genome sequencing analysis, the molecular transmission characteristics of the tet(X4) were preliminarily clarified, providing a scientific basis for guiding veterinary clinical use in this area, as well as risk assessment and prevention of the transfer and spread of tigecycline resistant strains or genes.

The honeybee gut resistome and its role in antibiotic resistance dissemination

There is now general concern about widespread antibiotic resistance, and growing evidence indicates that gut microbiota is critical in providing antibiotic resistance. Honeybee is an important pollinator; the incidence of antibiotic resistance genes in honeybee gut causes potential risks to not only its own health but also to public and animal health, for its potential disseminator role, thus receiving more attention from the public. Recent analysis results reveal that the gut of honeybee serves as a reservoir of antibiotic resistance genes, probably due to antibiotics application history in beekeeping and horizontal gene transfer from the highly polluted environment. These antibiotic resistance genes accumulate in the honeybee gut and could be transferred to the pathogen, even having the potential to spread during pollination, tending, social interactions, etc. Newly acquired resistance traits may cause fitness reduction in bacteria whereas facilitating adaptive evolution as well. This review outlines the current knowledge about the resistome in honeybee gut and emphasizes its role in antibiotic resistance dissemination.

Regulatory fine-tuning of mcr-1 increases bacterial fitness and stabilises antibiotic resistance in agricultural settings

Antibiotic resistance tends to carry fitness costs, making it difficult to understand how resistance can be maintained in the absence of continual antibiotic exposure. Here we investigate this problem in the context of mcr-1, a globally disseminated gene that confers resistance to colistin, an agricultural antibiotic that is used as a last resort for the treatment of multi-drug resistant infections. Here we show that regulatory evolution has fine-tuned the expression of mcr-1, allowing E. coli to reduce the fitness cost of mcr-1 while simultaneously increasing colistin resistance. Conjugative plasmids have transferred low-cost/high-resistance mcr-1 alleles across an incredible diversity of E. coli strains, further stabilising mcr-1 at the species level. Regulatory mutations were associated with increased mcr-1 stability in pig farms following a ban on the use of colistin as a growth promoter that decreased colistin consumption by 90%. Our study shows how regulatory evolution and plasmid transfer can combine to stabilise resistance and limit the impact of reducing antibiotic consumption.

Linking within- and between-host scales for understanding the evolutionary dynamics of quantitative antimicrobial resistance

Understanding both the epidemiological and evolutionary dynamics of antimicrobial resistance is a major public health concern. In this paper, we propose a nested model, explicitly linking the within- and between-host scales, in which the level of resistance of the bacterial population is viewed as a continuous quantitative trait. The within-host dynamics is based on integro-differential equations structured by the resistance level, while the between-host scale is additionally structured by the time since infection. This model simultaneously captures the dynamics of the bacteria population, the evolutionary transient dynamics which lead to the emergence of resistance, and the epidemic dynamics of the host population. Moreover, we precisely analyze the model proposed by particularly performing the uniform persistence and global asymptotic results. Finally, we discuss the impact of the treatment rate of the host population in controlling both the epidemic outbreak and the average level of resistance, either if the within-host scale therapy is a success or failure. We also explore how transitions between infected populations (treated and untreated) can impact the average level of resistance, particularly in a scenario where the treatment is successful at the within-host scale.

Promotion of plasmid maintenance by heterogeneous partitioning of microbial communities

Transferable plasmids play a critical role in shaping the functions of microbial communities. Previous studies suggested multiple mechanisms underlying plasmid persistence and abundance. Here, we focus on the interplay between heterogeneous community partitioning and plasmid fates. Natural microbiomes often experience partitioning that creates heterogeneous local communities with reduced population sizes and biodiversity. Little is known about how population partitioning affects the plasmid fate through the modulation of community structure. By modeling and experiments, we show that heterogeneous community partitioning can paradoxically promote the persistence of a plasmid that would otherwise not persist in a global community. Among the local communities created by partitioning, a minority will primarily consist of members able to transfer the plasmid fast enough to support its maintenance by serving as a local plasmid haven. Our results provide insights into plasmid maintenance and suggest a generalizable approach to modulate plasmid persistence for engineering and medical applications.

Comparative genomics unveils extensive genomic variation between populations of Listeria species in natural and food-associated environments

Comprehending bacterial genomic variation linked to distinct environments can yield novel insights into mechanisms underlying differential adaptation and transmission of microbes across environments. Gaining such insights is particularly crucial for pathogens as it benefits public health surveillance. However, the understanding of bacterial genomic variation is limited by a scarcity of investigations in genomic variation coupled with different ecological contexts. To address this limitation, we focused on Listeria, an important bacterial genus for food safety that includes the human pathogen L. monocytogenes, and analyzed a large-scale genomic dataset collected by us from natural and food-associated environments across the United States. Through comparative genomics analyses on 449 isolates from the soil and 390 isolates from agricultural water and produce processing facilities representing L. monocytogenes, L. seeligeri, L. innocua, and L. welshimeri, we find that the genomic profiles strongly differ by environments within each species. This is supported by the environment-associated subclades and differential presence of plasmids, stress islands, and accessory genes involved in cell envelope biogenesis and carbohydrate transport and metabolism. Core genomes of Listeria species are also strongly associated with environments and can accurately predict isolation sources at the lineage level in L. monocytogenes using machine learning. We find that the large environment-associated genomic variation in Listeria appears to be jointly driven by soil property, climate, land use, and accompanying bacterial species, chiefly representing Actinobacteria and Proteobacteria. Collectively, our data suggest that populations of Listeria species have genetically adapted to different environments, which may limit their transmission from natural to food-associated environments.

Mixed strain pathogen populations accelerate the evolution of antibiotic resistance in patients

Antibiotic resistance poses a global health threat, but the within-host drivers of resistance remain poorly understood. Pathogen populations are often assumed to be clonal within hosts, and resistance is thought to emerge due to selection for de novo variants. Here we show that mixed strain populations are common in the opportunistic pathogen P. aeruginosa. Crucially, resistance evolves rapidly in patients colonized by multiple strains through selection for pre-existing resistant strains. In contrast, resistance evolves sporadically in patients colonized by single strains due to selection for novel resistance mutations. However, strong trade-offs between resistance and growth rate occur in mixed strain populations, suggesting that within-host diversity can also drive the loss of resistance in the absence of antibiotic treatment. In summary, we show that the within-host diversity of pathogen populations plays a key role in shaping the emergence of resistance in response to treatment.

The distribution of fitness effects of plasmid pOXA-48 in clinical enterobacteria

Antimicrobial resistance (AMR) in bacteria is a major public health problem. The main route for AMR acquisition in clinically important bacteria is the horizontal transfer of plasmids carrying resistance genes. AMR plasmids allow bacteria to survive antibiotics, but they also entail physiological alterations in the host cell. Multiple studies over the last few years have indicated that these alterations can translate into a fitness cost when antibiotics are absent. However, due to technical limitations, most of these studies are based on analysing new associations between plasmids and bacteria generated in vitro, and we know very little about the effects of plasmids in their native bacterial hosts. In this study, we used a CRISPR-Cas9-tool to selectively cure plasmids from clinical enterobacteria to overcome this limitation. Using this approach, we were able to study the fitness effects of the carbapenem resistance plasmid pOXA-48 in 35 pOXA-48-carrying isolates recovered from hospitalized patients. Our results revealed that pOXA-48 produces variable effects across the collection of wild-type enterobacterial strains naturally carrying the plasmid, ranging from fitness costs to fitness benefits. Importantly, the plasmid was only associated with a significant fitness reduction in four out of 35 clones, and produced no significant changes in fitness in the great majority of isolates. Our results suggest that plasmids produce neutral fitness effects in most native bacterial hosts, helping to explain the great prevalence of plasmids in natural microbial communities.

A mathematician's guide to plasmids: an introduction to plasmid biology for modellers

Plasmids, extrachromosomal DNA molecules commonly found in bacterial and archaeal cells, play an important role in bacterial genetics and evolution. Our understanding of plasmid biology has been furthered greatly by the development of mathematical models, and there are many questions about plasmids that models would be useful in answering. In this review, we present an introductory, yet comprehensive, overview of the biology of plasmids suitable for modellers unfamiliar with plasmids who want to get up to speed and to begin working on plasmid-related models. In addition to reviewing the diversity of plasmids and the genes they carry, their key physiological functions, and interactions between plasmid and host, we also highlight selected plasmid topics that may be of particular interest to modellers and areas where there is a particular need for theoretical development. The world of plasmids holds a great variety of subjects that will interest mathematical biologists, and introducing new modellers to the subject will help to expand the existing body of plasmid theory.

There are three major Neisseria gonorrhoeae β-lactamase plasmid variants which are associated with specific lineages and carry distinct TEM alleles

Neisseria gonorrhoeae is a significant threat to global health with an estimated incidence of over 80 million cases each year and high levels of antimicrobial resistance. The gonococcal β-lactamase plasmid, pbla, carries the TEM β-lactamase, which requires only one or two amino acid changes to become an extended-spectrum β-lactamase (ESBL); this would render last resort treatments for gonorrhoea ineffective. Although pbla is not mobile, it can be transferred by the conjugative plasmid, pConj, found in N. gonorrhoeae. Seven variants of pbla have been described previously, but little is known about their frequency or distribution in the gonococcal population. We characterised sequences of pbla variants and devised a typing scheme, Ng_pblaST that allows their identification from whole genome short-read sequences. We implemented Ng_pblaST to assess the distribution of pbla variants in 15 532 gonococcal isolates. This demonstrated that only three pbla variants commonly circulate in gonococci, which together account for >99 % of sequences. The pbla variants carry different TEM alleles and are prevalent in distinct gonococcal lineages. Analysis of 2758 pbla-containing isolates revealed the co-occurrence of pbla with certain pConj types, indicating co-operativity between pbla and pConj variants in the spread of plasmid-mediated AMR in N. gonorrhoeae. Understanding the variation and distribution of pbla is essential for monitoring and predicting the spread of plasmid-mediated β-lactam resistance in N. gonorrhoeae.

Insertion Sequences Determine Plasmid Adaptation to New Bacterial Hosts

Plasmids facilitate the vertical and horizontal spread of antimicrobial resistance genes between bacteria. The host range and adaptation of plasmids to new hosts determine their impact on the spread of resistance. In this work, we explore the mechanisms driving plasmid adaptation to novel hosts in experimental evolution. Using the small multicopy plasmid pB1000, usually found in Pasteurellaceae, we studied its adaptation to a host from a different bacterial family, Escherichia coli. We observed two different mechanisms of adaptation. One mechanism is single nucleotide polymorphisms (SNPs) in the origin of replication (oriV) of the plasmid, which increase the copy number in E. coli cells, elevating the stability, and resistance profile. The second mechanism consists of two insertion sequences (ISs), IS1 and IS10, which decrease the fitness cost of the plasmid by disrupting an uncharacterized gene on pB1000 that is harmful to E. coli. Both mechanisms increase the stability of pB1000 independently, but only their combination allows long-term maintenance. Crucially, we show that the mechanisms have a different impact on the host range of the plasmid. SNPs in oriV prevent the replication in the original host, resulting in a shift of the host range. In contrast, the introduction of ISs either shifts or expands the host range, depending on the IS. While IS1 leads to expansion, IS10 cannot be reintroduced into the original host. This study gives new insights into the relevance of ISs in plasmid-host adaptation to understand the success in spreading resistance. IMPORTANCE ColE1-like plasmids are small, mobilizable plasmids that can be found across at least four orders of Gammaproteobacteria and are strongly associated with antimicrobial resistance genes. Plasmid pB1000 carries the gene bla_ROB-1, conferring high-level resistance to penicillins and cefaclor. pB1000 has been described in various species of the family Pasteurellaceae, for example, in Haemophilus influenzae, which can cause diseases such as otitis media, meningitis, and pneumonia. To understand the resistance spread through horizontal transfer, it is essential to study the mechanisms of plasmid adaptation to novel hosts. In this work we identify that a gene from pB1000, which encodes a peptide that is toxic for E. coli, and the low plasmid copy number (PCN) of pB1000 in E. coli cells are essential targets in the described plasmid-host adaptation and therefore limit the spread of pB1000-encoded bla_ROB-1. Furthermore, we show how the interplay of two adaptation mechanisms leads to successful plasmid maintenance in a different bacterial family.

Greater Invasion and Persistence of mcr-1-Bearing Plasmids in Escherichia coli than in Klebsiella pneumoniae

The emergence of the plasmid-borne polymyxin resistance gene mcr-1 threatens the clinical utility of last-line polymyxins. Although mcr-1 has disseminated to various Enterobacterales species, the prevalence of mcr-1 is the highest among Escherichia coli isolates while remaining low in Klebsiella pneumoniae. The reason for such a difference in prevalence has not been investigated. In this study, we examined and compared the biological characteristics of various mcr-1 plasmids in these two bacterial species. Although mcr-1-bearing plasmids were stably maintained in both E. coli and K. pneumoniae, the former presented itself to be superior by demonstrating a fitness advantage while carrying the plasmid. The inter- and intraspecies transferability efficiencies were evaluated for common mcr-1-harboring plasmids (IncX4, IncI2, IncHI2, IncP, and IncF types) with native E. coli and K. pneumoniae strains as donors. Here, we found that the conjugation frequencies of mcr-1 plasmids were significantly higher in E. coli than in K. pneumoniae, regardless of the donor species and Inc types of the mcr-1 plasmids. Plasmid invasion experiments revealed that mcr-1 plasmids displayed greater invasiveness and stability in E. coli than in K. pneumoniae. Moreover, K. pneumoniae carrying mcr-1 plasmids showed a competitive disadvantage when cocultured with E. coli. These findings indicate that mcr-1 plasmids could spread more easily among E. coli than among K. pneumoniae isolates and that mcr-1 plasmid-carrying E. coli has a competitive advantage over K. pneumoniae, leading to E. coli being the main mcr-1 reservoir. IMPORTANCE As infections caused by multidrug-resistant "superbugs" are increasing globally, polymyxins are often the only viable therapeutic option. Alarmingly, the wide spread of the plasmid-mediated polymyxin resistance gene mcr-1 is restricting the clinical utility of this last-line treatment option. With this, there is an urgent need to investigate the factors contributing to the spread and persistence of mcr-1-bearing plasmids in the bacterial community. Our research highlights that the higher prevalence of mcr-1 in E. coli than in K. pneumoniae is attributed to the greater transferability and persistence of mcr-1-bearing plasmid in the former species. By gaining these important insights into the persistence of mcr-1 in different bacterial species, we will be able to formulate effective strategies to curb the spread of mcr-1 and prolong the clinical life span of polymyxins.

Factors associated with plasmid antibiotic resistance gene carriage revealed using large-scale multivariable analysis

Plasmids are major vectors of bacterial antibiotic resistance, but understanding of factors associated with plasmid antibiotic resistance gene (ARG) carriage is limited. We curated > 14,000 publicly available plasmid genomes and associated metadata. Duplicate and replicate plasmids were excluded; where possible, sample metadata was validated externally (BacDive database). Using Generalised Additive Models (GAMs) we assessed the influence of 12 biotic/abiotic factors (e.g. plasmid genetic factors, isolation source, collection date) on ARG carriage, modelled as a binary outcome. Separate GAMs were built for 10 major ARG types. Multivariable analysis indicated that plasmid ARG carriage patterns across time (collection years), isolation sources (human/livestock) and host bacterial taxa were consistent with antibiotic selection pressure as a driver of plasmid-mediated antibiotic resistance. Only 0.42% livestock plasmids carried carbapenem resistance (compared with 12% human plasmids); conversely, tetracycline resistance was enriched in livestock vs human plasmids, reflecting known prescribing practices. Interpreting results using a timeline of ARG type acquisition (determined by literature review) yielded additional novel insights. More recently acquired ARG types (e.g. colistin and carbapenem) showed increases in plasmid carriage during the date range analysed (1994-2019), potentially reflecting recent onset of selection pressure; they also co-occurred less commonly with ARGs of other types, and virulence genes. Overall, this suggests that following acquisition, plasmid ARGs tend to accumulate under antibiotic selection pressure and co-associate with other adaptive genes (other ARG types, virulence genes), potentially re-enforcing plasmid ARG carriage through co-selection.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Biofilms preserve the transmissibility of a multi-drug resistance plasmid

Self-transmissible multidrug resistance (MDR) plasmids are a major health concern because they can spread antibiotic resistance to pathogens. Even though most pathogens form biofilms, little is known about how MDR plasmids persist and evolve in biofilms. We hypothesize that (i) biofilms act as refugia of MDR plasmids by retaining them in the absence of antibiotics longer than well-mixed planktonic populations and that (ii) the evolutionary trajectories that account for the improvement of plasmid persistence over time differ between biofilms and planktonic populations. In this study, we evolved Acinetobacter baumannii with an MDR plasmid in biofilm and planktonic populations with and without antibiotic selection. In the absence of selection, biofilm populations were better able to maintain the MDR plasmid than planktonic populations. In planktonic populations, plasmid persistence improved rapidly but was accompanied by a loss of genes required for the horizontal transfer of plasmids. In contrast, in biofilms, most plasmids retained their transfer genes, but on average, plasmid, persistence improved less over time. Our results showed that biofilms can act as refugia of MDR plasmids and favor the horizontal mode of plasmid transfer, which has important implications for the spread of MDR.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

A Trade-Off for Maintenance of Multidrug-Resistant IncHI2 Plasmids in Salmonella enterica Serovar Typhimurium through Adaptive Evolution

Understanding the fitness costs associated with plasmid carriage is a key to better understanding the mechanisms of plasmid maintenance in bacteria. In the current work, we performed multiple serial passages (63 days, 627.8 generations) to identify the compensatory mechanisms that Salmonella enterica serovar Typhimurium ATCC 14028 utilized to maintain the multidrug-resistant (MDR) IncHI2 plasmid pJXP9 in the presence and absence of antibiotic selection. The plasmid pJXP9 was maintained for hundreds of generations even without drug exposure. Endpoint evolved (the endpoint of evolution) S. Typhimurium bearing evolved plasmids displayed decreased growth lag times and a competitive advantage over ancestral pJXP9 plasmid-carrying ATCC 14028 strains. Genomic and transcriptomic analyses revealed that the fitness costs of carrying pJXP9 were derived from both specific plasmid genes and particularly the MDR regions and conjugation transfer region I and conflicts resulting from chromosome-plasmid gene interactions. Correspondingly, plasmid deletions of these regions could compensate for the fitness cost that was due to the plasmid carriage. The deletion extent and range of large fragments on the evolved plasmids, as well as the trajectory of deletion mutation, were related to the antibiotic treatment conditions. Furthermore, it is also adaptive evolution that chromosomal gene mutations and altered mRNA expression correlated with changed physiological functions of the bacterium, such as decreased flagellar motility, increased oxidative stress, and fumaric acid synthesis but increased Cu resistance in a given niche. Our findings indicated that plasmid maintenance evolves via a plasmid-bacterium adaptative evolutionary process that is a trade-off between vertical and horizontal transmission costs along with associated alterations in host bacterial physiology. IMPORTANCE The current idea that compensatory evolution processes can account for the "plasmid paradox" phenomenon associated with the maintenance of large costly plasmids in host bacteria has attracted much attention. Although many compensatory mutations have been discovered through various plasmid-host bacterial evolution experiments, the basis of the compensatory mechanisms and the nature of the bacteria themselves to address the fitness costs remain unclear. In addition, the genetic backgrounds of plasmids and strains involved in previous research were limited and clinical drug resistance such as the poorly understood compensatory evolution among clinically dominant multidrug-resistant plasmids or clones was rarely considered. The IncHI2 plasmid is widely distributed in Salmonella Typhimurium and plays an important role in the emergence and rapid spread of its multidrug resistance. In this study, the predominant multidrug-resistant IncHI2 plasmid pJXP9 and the standard Salmonella Typhimurium ATCC 14028 bacteria were used for evolution experiments under laboratory conditions. Our findings indicated that plasmid maintenance through experimental evolution of plasmid-host bacteria is a trade-off between increasing plasmid vertical transmission and impairing its horizontal transmission and bacterial physiological phenotypes, in which compensatory mutations and altered chromosomal expression profiles collectively contribute to alleviating plasmid-borne fitness cost. These results provided potential insights into understanding the relationship of coexistence between plasmids encoding antibiotic resistance and their bacterial hosts and provided a clue to the adaptive forces that shaped the evolution of these plasmids within bacteria and to predicting the evolution trajectory of antibiotic resistance.

Persistence of plasmid and tet(X4) in an Escherichia coli isolate coharboring blaNDM-5 and mcr-1 after acquiring an IncFII tet(X4)-positive plasmid

The prevalence of plasmid-mediated tigecycline resistance gene tet(X4) is presenting an increasing trend. Once tet(X4)-bearing plasmids are captured by multidrug-resistant bacteria, such as bla_NDM and mcr-coharboring bacteria, it will promote bacteria to develop an ultra-broad resistance spectrum, limiting clinical treatment options. However, little is known about the destiny of such bacteria or how they will evolve in the future. Herein, we constructed a multidrug-resistant bacteria coharboring tet(X4), bla_NDM-5, and mcr-1 by introducing a tet(X4)-bearing plasmid into a bla_NDM-5 and mcr-1 positive E. coli strain. Subsequently, the stability of tet(X4) and the plasmid was measured after being evolved under tigecycline or antibiotic-free circumstance. Interestingly, we observed both tet(X4)-bearing plasmids in tigecycline treated strains and non-tigecycline treated strains were stable, which might be jointly affected by the increased conjugation frequency and the structural alterations of the tet(X4)-positive plasmid. However, the stability of tet(X4) gene showed different scenarios in the two types of evolved strains. The tet(X4) gene in non-tigecycline treated strains was stable whereas the tet(X4) gene was discarded rapidly in tigecycline treated strains. Accordingly, we found the expression levels of tet(X4) gene in tigecycline-treated strains were several times higher than in non-tigecycline treated strains and ancestral strains, which might in turn impose a stronger burden on the host bacteria. SNPs analysis revealed that a myriad of mutations occurred in genes involving in conjugation transfer, and the missense mutation of marR gene in chromosome of tigecycline treated strains might account for the completely different stability of tet(X4)-bearing plasmid and tet(X4) gene. Collectively, these findings shed a light on the possibility of the emergence of multidrug resistant bacteria due to the transmission of tet(X4)-bearing plasmid, and highlighted that the antibiotic residues may be critical to the development of such bacteria.

IS26-mediated plasmid reshuffling results in convergence of toxin-antitoxin systems but loss of resistance genes in XDR Klebsiella pneumoniae from a chronic infection

Carbapenem-resistant Enterobacterales pose an urgent threat to human health worldwide. Klebsiella pneumoniae sequence type (ST) 14, initially identified in the Middle East and South-Asia and co-harbouring the carbapenemase genes bla _OXA-232 and bla _NDM-1, is now emerging globally. One such strain was detected in the USA in 2013 from a patient initially treated in India that also carried armA, a 16S rRNA methyltransferase that confers resistance to all clinically relevant aminoglycosides. Genetic and phenotypic changes were observed in 14 serial isolates collected from this chronically infected patient. The index isolate carried five plasmids, including an IncFIB-IncHI1B (harbouring armA and bla _NDM-1), an IncFIA (bla _CTX-M-15) and a ColE-like (bla _OXA-232), and was extensively resistant to antibiotics. Four years later, a subsequent isolate had accumulated 34 variants, including a loss-of-function mutation in romA, resulting in tigecycline non-susceptibility. Importantly, this isolate now only carried two plasmids, including a large mosaic molecule made of fragments, all harbouring distinct toxin-antitoxin systems, from three of the canonical plasmids. Of the original acquired antibiotic resistance genes, this isolate only retained bla _CTX-M-15, and as a result susceptibility to the carbapenems and amikacin was restored. Long-read sequencing of a subset of five representative isolates, collected between 2013 and 2017, allowed for the elucidation of the complex plasmid patterns and revealed the role of IS26-mediated plasmid reshuffling in the evolution of this clone. Such investigations of the mechanisms underlying plasmid stability, together with global and local surveillance programmes, are key to a better understanding of plasmid host range and dissemination.

A limit on the evolutionary rescue of an Antarctic bacterium from rising temperatures

Climate change is gradual, but it can also cause brief extreme heat waves that can exceed the upper thermal limit of any one organism. To study the evolutionary potential of upper thermal tolerance, we evolved the cold-adapted Antarctic bacterium Pseudoalteromonas haloplanktis to survive at 30°C, beyond its ancestral thermal limit. This high-temperature adaptation occurred rapidly and in multiple populations. It involved genomic changes that occurred in a highly parallel fashion and mitigated the effects of protein misfolding. However, it also confronted a physiological limit, because populations failed to grow beyond 30°C. Our experiments aimed to facilitate evolutionary rescue by using a small organism with large populations living at temperatures several degrees below their upper thermal limit. Larger organisms with smaller populations and living at temperatures closer to their upper thermal tolerances are even more likely to go extinct during extreme heat waves.

Virulence plasmid pINV as a genetic signature for Shigella flexneri phylogeny

Shigella flexneri is a major health burden in low- and middle-income countries, where it is a leading cause of mortality associated with diarrhoea in children, and shows an increasing incidence among travellers and men having sex with men. Like all Shigella spp., S. flexneri has evolved from commensal Escherichia coli following the acquisition of a large plasmid pINV, which contains genes essential for virulence. Current sequence typing schemes of Shigella are based on combinations of chromosomal genetic loci, since pINV-encoded virulence genes are often lost during growth in the laboratory, making these elements inappropriate for sequence typing. By performing comparative analysis of pINVs from S. flexneri strains isolated from different geographical regions and belonging to different serotypes, we found that in contrast to plasmid-encoded virulence genes, plasmid maintenance genes are highly stable pINV-encoded elements. For the first time, to our knowledge, we have developed a S. flexneri plasmid multilocus sequence typing (pMLST) method based on different combinations of alleles of the vapBC and yacAB toxin–antitoxin (TA) systems, and the parAB partitioning system. This enables typing of S. flexneri pINV plasmids into distinct ‘virulence sequence types’ (vSTs). Furthermore, the phylogenies of vST alleles and bacterial host core genomes suggests an intimate co-evolution of pINV with the chromosome of its bacterial host, consistent with previous findings. This work demonstrates the potential of plasmid maintenance loci as genetic characteristics to study as well as to trace the molecular phylogenesis of S. flexneri pINV and the phylogenetic relationship of this plasmid with its bacterial host.

High plasmidome diversity of extended-spectrum beta-lactam-resistant Escherichia coli isolates collected during one year in one community hospital

Plasmid-encoded antibiotic resistance encompasses many classes of currently used antibiotics. In globally distributed Escherichia coli lineages plasmids, which spread via horizontal gene transfer, are responsible for the dissemination of genes encoding extended-spectrum β-lactamases (ESBL). In this study, we combined 2nd and 3rd generation sequencing techniques to reconstruct the plasmidome of overall 97 clinical ESBL-E. coli isolates. Our results highlight the enormous plasmid diversity in respect to size, replicon-type and genetic content. Furthermore, we emphasize the diverse plasmid distribution patterns among the clinical isolates and the high intra- and extracellular mobility potential of resistance conferring genes. While the majority of resistance conferring genes were located on large plasmids of known replicon type, small cryptic plasmids seem to be underestimated resistance gene vectors. Our results contribute to a better understanding of the dissemination of resistance-conferring genes through horizontal gene transfer as well as clonal spread.

Recombination resolves the cost of horizontal gene transfer in experimental populations of Helicobacter pylori

Horizontal gene transfer (HGT)—the transfer of DNA between lineages—is responsible for a large proportion of the genetic variation that contributes to evolution in microbial populations. While HGT can bring beneficial genetic innovation, the transfer of DNA from other species or strains can also have deleterious effects. In this study, we evolve populations of the bacteria Helicobacter pylori and use DNA sequencing to identify over 40,000 genetic variants transferred by HGT. We measure the cost of many of these and find that both strongly beneficial mutations and deleterious mutations are genetic variants transferred by natural transformation. Importantly, we also show how recombination that separates linked beneficial and deleterious mutations resolves the cost of HGT.

Horizontal gene transfer (HGT) is important for microbial evolution, yet we know little about the fitness effects and dynamics of horizontally transferred genetic variants. In this study, we evolve laboratory populations of Helicobacter pylori, which take up DNA from their environment by natural transformation, and measure the fitness effects of thousands of transferred genetic variants. We find that natural transformation increases the rate of adaptation but comes at the cost of significant genetic load. We show that this cost is circumvented by recombination, which increases the efficiency of selection by decoupling deleterious and beneficial genetic variants. Our results show that adaptation with HGT, pervasive in natural microbial populations, is shaped by a combination of selection, recombination, and genetic drift not accounted for in existing models of evolution.