Next Article in Journal
Self-Amplifying RNA Vaccine Candidates: Alternative Platforms for mRNA Vaccine Development
Previous Article in Journal
Selene-Ethylenelacticamides and N-Aryl-Propanamides as Broad-Spectrum Leishmanicidal Agents
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 12
Issue 1
10.3390/pathogens12010137
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Evolution Compensated for the Plasmid Fitness Costs Brought by Specific Genetic Conflicts

by
Feifeng Li
1,†,
Jiong Wang
1,†,
Ying Jiang
2,
Yingyi Guo
1,
Ningjing Liu
1,
Shunian Xiao
1,
Likang Yao
1,
Jiahui Li
1,
Chuyue Zhuo
1,
Nanhao He
1,
Baomo Liu
3,* and
Chao Zhuo
1,*
1
State Key Laboratory of Respiratory Disease, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou 510030, China
2
Department of Respiratory Medicine, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou 510030, China
3
Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Sun Yat-sen Univesity, Guangzhou 510030, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Pathogens 2023, 12(1), 137; https://doi.org/10.3390/pathogens12010137
Submission received: 29 November 2022 / Revised: 30 December 2022 / Accepted: 9 January 2023 / Published: 13 January 2023
(This article belongs to the Section Bacterial Pathogens)
Browse Figures
Review Reports Versions Notes

Abstract

:
New Delhi metallo-β-lactamase (NDM)-carrying IncX3 plasmids is important in the transmission of carbapenem resistance in Escherichia coli. Fitness costs related to plasmid carriage are expected to limit gene exchange; however, the causes of these fitness costs are poorly understood. Compensatory mutations are believed to ameliorate plasmid fitness costs and enable the plasmid’s wide spread, suggesting that such costs are caused by specific plasmid–host genetic conflicts. By combining conjugation tests and experimental evolution with comparative genetic analysis, we showed here that the fitness costs related to ndm/IncX3 plasmids in E. coli C600 are caused by co-mutations of multiple host chromosomal genes related to sugar metabolism and cell membrane function. Adaptive evolution revealed that mutations in genes associated with oxidative stress, nucleotide and short-chain fatty acid metabolism, and cell membranes ameliorated the costs associated with plasmid carriage. Specific genetic conflicts associated with the ndm/IncX3 plasmid in E. coli C600 involve metabolism and cell-membrane-related genes, which could be ameliorated by compensatory mutations. Collectively, our findings could explain the wide spread of IncX3 plasmids in bacterial genomes, despite their potential cost.

1. Introduction

The emergence and dissemination of carbapenem resistance conferred by New Delhi metallo-β-lactamase (NDM) has become a major concern worldwide [1]. Of the plasmids reported to carry blandm, IncX3-type plasmids are the most common. The ndm/IncX3 plasmids have been detected in Asia, Europe, and North America, though strains from China carry the highest proportion of ndm/IncX3 (68.4%) [2,3].
Plasmids reduce the competitiveness of the host bacteria and so confer a fitness cost [4], which is considered a major barrier to plasmid maintenance in microbial populations [5]. This cost in fitness has been attributed to plasmid properties, such as the number of copies, size of the plasmid, and related metabolic burden [4,5,6,7]. The ndm/IncX3 plasmid was considered to confer a low fitness cost because of its low copy number, minor metabolic burden, and stable genetic backbone [8,9,10,11]; however, high fitness costs have been reported in some E. coli strains [12,13]. The mechanism of the diversity of fitness costs while carrying the ndm/IncX3 plasmid remains unclear and is potentially unrelated to the general properties of the plasmid.
Genetic conflicts, which emerge from plasmid–chromosome gene interactions as plasmid entry into the host, are proposed by some researchers as vital mechanisms for plasmid-related fitness costs [14,15]. Plasmids are independently replicating genetic elements, which differ from chromosomes. Genetic conflicts arise when plasmids and chromosomes attempt to independently maximize their own fitness [16]. Such genetic conflicts have manifested as mutations, such as missense variant, deletion, and stop-gained, occurring in chromosomes or plasmids. Frequent conflict between plasmid replication and the host chromosome comes at the cost of host survival, such as the fitness cost reported for Salmonella enterica serovar Typhimurium ATCC14028 carrying the MDR IncHI2 plasmid pJXP9 manifesting as chromosomal mutations [17]. These genetic conflicts would likely be strain-specific, depending on the genetic background of the strain. Hall et al. [14] proposed that specific genetic conflicts were responsible for the fitness cost reported for Pseudomonas fluorescens SBW25 carrying the megaplasmids pQBR57 and pQBR103; the specific genetic conflict occurred as the SOS-response-related gene PFLU4242 mutation. However, there have been no studies on the relationship between specific genetic conflicts and the fitness cost of carrying ndm/IncX3 plasmids. Our work has previously revealed that the ndm/IncX3 plasmids are highly conserved [9], and thus we suggest that specific genetic conflicts that may occur in the host chromosome confer a diverse fitness cost.
Following plasmid-related fitness cost emerging, this could have one of two outcomes: either that plasmids will be lost, which is clearly inconsistent with the current epidemic status of plasmids, or that plasmids and host strains can undergo compensatory evolution, which reduces the cost of plasmid carriage and allows for plasmid persistence. Similarly, for the plasmid to persist, the fitness cost imposed by genetic conflicts conferring targeted deletions and mutations should be compensated by adaptive evolution [14,17,18]. Therefore, one may expect coevolution between the host chromosomes and plasmids to correspond with a reduction in the deleterious effects of the plasmid. We suggest that host strains carrying the ndm/IncX3 plasmid initially have a fitness cost, which is compensated by adaptive evolution, and thus may lead to diverse levels of fitness cost, but the mechanism of the diverse fitness cost after coevolution remains unknown.
This study investigates the relationship between diverse fitness costs associated with the ndm/IncX3 plasmid and specific genetic conflicts, and how adaptive evolution compensates the fitness cost caused by this specific genetic conflict. We evaluated the mechanisms underlying the fitness costs specifically related to the ndm/IncX3 plasmid pHN330 in E. coli C600 and performed experimental evolution to explore the potential for compensatory evolution and its mechanism. We found that carrying the ndm/IncX3 plasmid generated an elevated fitness cost for E. coli C600, suggesting the occurrence of a specific genetic conflict. Adaptive evolution can compensate for the specific plasmid–host chromosome conflicts, ameliorating the cost of the C600 carrying the ndm/IncX3 plasmid.

2. Materials and Methods

2.1. Bacterial Strains and Plasmid

The plasmid pHN330 was isolated in our laboratory as reported in our previous study [19]. We performed conjugation experiments and successfully obtained E. coli J53 and C600 transformants containing pHN330; Table 1 summarizes the strain information. Susceptibility to antibiotics was assessed using a standard microdilution method according to the CLSI guidelines [20].

2.2. Growth Curve Assay

The growth curves of the strains were determined by measuring the optical density (OD) at 600 nm (BioTek, Vermont, USA). Briefly, cultures grown in Luria–Bertani (LB) broth overnight at 37 °C were diluted 1:1000 in fresh LB broth and grown at 37 °C with shaking (200 rpm) for 24 h. The OD600 of the culture was measured once every hour for a total of 24 h [21].

2.3. Competition Experiments

To assess the fitness impact of pHN330 carriage on the host, transformants and host strains were paired to perform competition experiments. The overnight bacterial suspension of paired strains was adjusted to 0.5 McFarland, mixed at a ~1:1 ratio, and 50 μL of the mixture was seeded into 5 mL of LB broth for 24 h of incubation (37 °C, 200 rpm). At 0 and 24 h, 100 μL of the culture was diluted and the diluted sample was transferred to LB agar or LB agar supplemented with meropenem (1 mg/L). After incubation at 37 °C for 18–24 h, colony forming units (CFUs) were counted and the competitive index (CI) was calculated as previously described [22]: CI = (transconjugant CFU/original CFU)/(inoculated transconjugant CFU/inoculated original CFU).

2.4. Experimental Evolution

Briefly, the C600T strain was inoculated into 5 mL of LB broth with or without imipenem (2 mg/L). These cultures were considered generation 0 for each lineage. We transferred 5 μL of the culture daily to 5 mL fresh LB broth with imipenem (2 mg/L) or without antibiotics over 50 d (10 generations per day, 500 generations total). PCR analysis confirmed that all the progeny strains were ndm-positive [19].

2.5. Whole-Genome Sequencing and Bioinformatic Analysis

Whole-genome sequencing was performed on all strains. Briefly, single colonies were grown for 16 h in LB broth at 37 °C and bacterial DNA was extracted using a bacterial DNA extraction kit (Omega Bio-Tek, USA) according to the manufacturer’s instructions. Whole-genome sequencing was performed by Novegene (Novogene BioTech, Beijing, China) on an Illumina Novaseq platform. Genomes were assembled using Shovill version 1.0.4 [23] and annotated using Prokka software [24]. Single nucleotide polymorphism (SNP) analysis was performed using Snippy [25]. The online tool DAVID [26] was used to annotate SNP genes according to biological function and cellular components via Gene Ontology analysis [27]. The genome of the E. coli strains were submitted to GenBank under accession number PRJNA907472.

2.6. Statistical Analysis

Data were analyzed using SPSS software version 25.0 (IBM SPSS, Chicago, IL, USA). CIs and growth rates were compared between groups using a paired t-test, with statistical significance set at p < 0.05.

3. Results

3.1. Fitness Cost of Carrying IncX3 Plasmid

The pHN330 plasmid harboring blandm and the studied strains acquired resistance to carbapenems after carrying the IncX3 plasmid. We performed antibiotic susceptibility tests (Table 2) and found that all strains carrying pHN330 were resistant to imipenem and meropenem, with MIC values 32- to 512-fold higher than the control strains.
Competition experiments and growth curves were used to assess the fitness cost of carrying pHN330 (Figure 1). We found that when J53 was the host, there were no significant differences in growth rates between the transformants and control bacteria (p > 0.05); when C600 was used as the host, the growth rate of the transformants was lower than that of control bacteria (2.47 vs. 2.16, p < 0.0001).
Further competition experiments (Figure 2) showed that the fitness cost of carrying the IncX3 plasmid was minor for the J53 strain (CI: 0.92 ± 0.03) but carrying the pHN330 had a clear fitness cost for the C600 strain (CI: 0.15 ± 0.04).

3.2. Comparative Analysis of E. coli C600 and E. coli C600/pHN330

To explore the mechanism conferring the fitness cost related to the IncX3 plasmid, we analyzed whether SNPs occurred in the plasmid or chromosomal genes by comparative genomics. We detected several SNPs in J53 (8) and C600 (29), but none were observed in the IncX3 plasmid (Table 3 and Table S1).
After J53 obtained the IncX3 plasmid, no SNPs were found in the coding sequence (CDS) regions, and all SNP positions were located in intergenic regions (Table 3). In contrast, after C600 obtained pHN330, 10 non-synonymous SNPs occurred in the CDS region, which were located within the ptrA, lpdA, bamA, melA, gatA, mdoB, yihV, mntH, trkA, and ydfK genes, which we termed SNP genes (Table 4).
Functional classification and clustering analyses with DAVID of all SNP genes in the C600 strain (Tables S3 and S4) indicated that the bamA, mdoB, mntH, and trkA genes were associated with the cell membrane, and lpdA, ptrA, gatA, yihV, ydfK, and melA were related to metabolism.

3.3. IncX3 Plasmid Maintenance in E. coli C600 after 500 Passages

Plasmid–host coevolution can gradually reduce the fitness cost of plasmids. To verify whether the same phenomenon exists in C600T, the experimental evolution of C600T was assessed over 50 d without antibiotics or with imipenem 2 mg/L. pHN330 was consistently present in all progeny strains, as confirmed by the PCR assay. The resistance profiles of C600T before and after laboratory evolution (Table 2) showed that the MIC of meropenem was unchanged after passage without antibiotics but decreased 4-fold after passage with imipenem (from 64 mg/L to 16 mg/L), while the MIC of imipenem decreased 2-fold after passage with and without antibiotics (from 64 mg/L to 32 mg/L).

3.4. Fitness Cost of Progeny Strains Decreased Compared to C600T Strain

After passages without antibiotics or with imipenem (2 mg/L), the growth rates of C600N500 and C600Y500 significantly increased compared to that of the C600T strain but did not recover to C600 levels. The CIs of C600N500 (0.22 vs. 0.14, t = −4.793, p = 0.003) and C600Y500 (0.28 vs. 0.14, t = −7.328, p = 0.001) were clearly higher than that of C600T (Figure 3).

3.5. SNPs Occurring in the CDS Region Reduced the Fitness Cost of C600T

To investigate the mechanism underlying the decreased fitness cost of IncX3 plasmid-carrying strains after evolution, progeny strains and transformants were subjected to comparative genomic analysis (Table 5 and Table S2). SNPs were present in all progeny strains (C600N500:28; C600Y500:27). After passages without antibiotics or with imipenem (2 mg/L), the progeny strains had more SNPs in the CDS region (C600N500:17; C600Y500:11).
As shown in Table 6, C600N500 had 14 SNP genes, gudD, araC, dgt, bamA, rfaF, lptF, ycjX, bcsG, pepQ, mshA, sfmD, paaJ, mreC, and tfaE, and C600Y500 had 10 SNP genes, gloA, ydiF, fliY, hycF, caiA, bamA, rlmB, sbp, acnA, and yghQ. A reverse mutation in the bamA gene was observed in the absence of antibiotic passages and with imipenem passages (1882T > G p.Phe628Val, after passages, reverse).
Further gene functional classification and clustering analyses of the non-synonymous SNP genes were performed in the progeny strains. Among the SNP genes that occurred in C600N500, bamA, bcsG, lptF, mreC, sfmD, mshA, and rfaF were related to the cell membrane, and gudD, araC, pepQ, paaJ, and dgt were related to metabolism (Tables S3 and S5). Among the SNP genes found in C600Y500, fliY, bamA, yghQ, and sbp were related to the cell membrane, and caiA, gloA, acnA, rlmB, hycF, and ydiF were related to metabolism (Tables S3 and S6). The number of SNP genes and gene clusters in the progeny strains and transformants are shown in Figure 4.

4. Discussion

The ndm/IncX3 plasmid is thought to confer a low fitness cost that contributes to its worldwide occurrence [10,11]. However, some studies have also reported high fitness costs [12,13], though the variety of fitness costs and the mechanism involved remain unclear. Our group has been exploring IncX3 plasmid prevalence for several years [8,9,19]. In the present study, we found that the ndm/IncX3 plasmid led to an elevated fitness cost in the C600 strain, for which specific genetic conflicts were responsible and related to SNPs in genes of the cell membrane and sugar metabolism. Adaptive evolution can compensate for specific plasmid–host chromosome conflicts through SNPs in genes related to nucleotide metabolism, strain biofilm formation, LPS synthesis, antioxidant stress, cytoskeleton, and short-chain fatty acid metabolism, ameliorating the cost of the C600 carrying ndm/IncX3 plasmid. This indicates that diversity in fitness costs is likely common when evaluating plasmid carriage across hosts, and genetic conflicts may provide relevant insights.
In the present study, plasmid-carrying E. coli J53 and C600 showed different fitness costs, which is consist with previous studies that indicated that ndm/IncX3 plasmids cause diverse fitness costs [12,13]. Subsequent comparative genomic analysis revealed that the fitness cost of carrying ndm/IncX3 plasmids is associated with several genetic mutations. Plasmid-carrying E. coli C600 had more SNPs in CDS region genes and fewer in intergenic regions on the chromosome compared to E. coli J53, indicating that mutations in intergenic regions incur a low effect on fitness cost. For instance, Buckner et al. [28] showed that the SNPs of Klebsiella pneumoniae Ecl8 carrying the pKpQIL plasmid were mainly in the intergenic regions with low fitness cost. Wang et al. [29] found that mutations in gyrB resulted in considerable costs, while carrying the OqxAB-encoding plasmids pHNGC59 and pHNFD436. Thus, mutations in functional CDS regions play an important role in generating fitness costs.
We analyze what specific genetic conflicts are related to the function of the SNP genes mutated in the CDS region of E. coli C600. The functions of the non-synonymous SNP genes in C600 (ptrA, lpdA, bamA, melA, gatA, mdoB, yihV, mntH, trkA, and ydfK) were related to the cell membrane and sugar metabolism (Table S4). Overall, the SNPs occurring in the metabolism-related genes (lpdA, melA, gatA, and yihV) reduce the efficiency with which E. coli C600 metabolizes different sugars, resulting in slower growth. The lpdA gene is a component of the pyruvate dehydrogenase complex, which codes for lipoamide dehydrogenase and has been reported to promote the growth rate of E. coli by glucose metabolism [30]. MelA, gatA, and yihV, which are involved in the metabolism of other sugars (galactitol, melibiose, and sulfoquinovose) in addition to glucose, have been reported to restrict bacteria from utilizing essential nutrients and reduce growth rates [31,32,33]. Mutants occurring in the cell-membrane-related genes (bamA, trkA, and mdoB) affected the stability of E. coli C600 membranes, reduced resistance to osmotic pressure, and promoted potassium ion imbalance, resulting in the decrease in bacterial growth. The bamA gene is an outer membrane protein assembly factor, which encodes a component of the BAM complex that is essential for the construction and maintenance of the outer membrane [34,35]. Phosphoglycerol transferase I—encoded by the mdoB gene—is involved in the biosynthesis of membrane-derived oligosaccharides and has been reported to affect the regulation of osmolarity and the stability of the bacterial membrane [36]. TrkA, the NAD-binding component of the Trk potassium transporter, is responsible for maintaining potassium ion homeostasis under hypotonic conditions [37]. Consequently, the co-mutations in these metabolic and cell-membrane-related genes reduced the ability of pHN330 carrying E. coli C600 to utilize different sugars, affected the membrane stability, decreased the ability to resist various environmental stresses, and caused a potassium ion imbalance, resulting in fitness costs.
Plasmid fitness costs related to bacterial metabolism and cell membrane function have been previously reported. Alvaro et al. [38] found that carrying different plasmids altered the expression of metabolic genes in Pseudomonas aeruginosa PAO1, which altered the regulation of energy production, lipid transport, and carbohydrate transport. E. coli DH10 carrying the IncA/C2 plasmid pAR060302 has also been reported to have an elevated fitness cost [39], specifically through the upregulated expression of the plasma membrane and active transmembrane transport genes [40].
In nature, the cost of carrying plasmids can be compensated through adaptive evolution [41]. In the present study, we found that the cost could be compensated by co-evolution, but the underlying mechanisms of host–plasmid coevolution compensating for the fitness cost are still unknown. Comparative genomic analysis between progeny strains and transformants revealed that bamA showed a reverse mutation in the compensated strains C600N500 and C600Y500. BamA is an essential component of the Omp85 protein superfamily and is well conserved in all Gram-negative bacteria [34,42]. BamA is important for the activity of the conserved β-barrel assembly machinery, which is responsible for the assembly of the outer membrane β-barrel conformation in Gram-negative bacteria [43,44]. BamA, the component of β-barrel membrane proteins, is proven to be vital for nutrient import, signaling, motility, and survival of Gram-negative bacteria [45]. Wu et al. found that the BamA protein depleted in E coli JCM166 results in cell death [46]. Thus, our results are consistent with the function of the bamA gene; mutants that occurred in the bamA gene resulted in a decreased growth rate and, after compensation, the bamA gene reverse mutant increased the growth rate of the C600N500 and C600Y500 strains. Future studies are needed to explore the role of the bamA gene in determining the fitness cost of another IncX3 plasmid bearing E. coli.
In addition to the bamA gene, GO clustering analysis suggested that other non-synonymous SNP genes occurred in progeny strains C600N500 and C600Y500, which were also related to the cell membrane, nucleotide, and short-chain fatty acid metabolism (Tables S5 and S6). The metabolism-related non-synonymous SNP gene dgt occurred in C600N500 and is associated with nucleotide metabolism. The dgt gene encodes for dGTPase, responsible for the hydrolysis of dGTP to deoxyguanosine and triphosphate. Mutants in the dgt gene can cause mismatches during DNA replication [47]. The cell-membrane-related non-synonymous SNP genes bcsG, lptF, rfaF, and mreC occurred in C600N500 and are associated with biofilm formation, LPS synthesis, and cell cytoskeletal formation. The gene bcsG, which is involved in the production of exopolysaccharide cellulose, has been reported to affect biofilm formation in E. coli [48]. The lptF and rfaF genes are essential for lipopolysaccharide formation, which is critical for cell motility, intestinal colonization, and biofilm formation, and encode for major extracellular polymers [49,50,51]. The gene mreC encodes for cytoskeletal proteins that regulate cell morphology and survival and is highly conserved in bacteria [52].
Adaptive evolution can compensate for fitness costs through altered metabolism and cell membrane function. Elias et al. found that the fitness cost of K. pneumoniae after obtaining ceftazidime-avibactam resistance can be compensated through adaptive evolution through the downregulated expression of phospholipids and LPS biogenesis and modification [53]. Previous studies have found that P. fluorescens carries pQBR103 at an elevated cost and can be compensated by a mutant in the gacA/gacS two-component system [54], which regulates energy metabolism and central intermediary metabolism [55]. Ma et al. found that after evolution, mutant DNA replication in the dnaK gene compensated for the fitness cost of E. coli carrying the mcr/IncHI plasmid [56].
Oxidative stress resistance and the metabolism of short-chain fatty acids were related to the adaptive evolution of C600Y500 under antibiotics. The ydiF gene encodes acetate-CoA transferase, which is a short-chain acyl-CoA involved in the metabolism of short-chain fatty acids [57]. FliY is an E. coli L-cysteine/L-cystine importer that cooperates with the L-cysteine exporter EAMA in the L-cysteine/L-cystine shuttle system involved in oxidative stress resistance [58] and has been reported to affect the growth of E. coli in LB broth [59]. The gloA gene encodes glyoxalase I involved in the glutathione (GSH)-dependent glyoxalase I pathway [60], which influences the metabolism of methylglyoxal; a high concentration of methylglyoxal has been reported to lead to bacterial growth stasis and even death [61]. The fliY, gloA, and acnA gene mutations in C600Y500, which are related to different pathways involved in oxidative stress resistance, may be related to the selection pressure of imipenem on C600Y500. Wang et al. found that imipenem selection pressure in Acinetobacter baumannii stimulated the strain to produce hydroxyl radicals [62]. Another study found that ampicillin induced the formation of hydroxyl radicals through oxidative damage, which led to bacterial death [63].
Our results suggest that while IncX3 plasmids are already known to have a low fitness cost, occasional elevated fitness costs can occur. This suggests that IncX3 plasmids could have initially been high-fitness-cost plasmids; however, through host–plasmid coevolution, the IncX3 plasmid became a low-fitness-cost plasmid. Our findings indicate that the host–plasmid adaptive coevolution process is a trade-off between plasmid and host chromosome fitness, in which the IncX3 plasmid is an approximate perfect plasmid, and thus the plasmid and host reach a compromise manifested as host chromosome mutations alleviating the plasmid-borne fitness cost. S. Typhimurium ATCC 14028 and IncHI2 plasmid pJXP9 plasmid evolutionary analysis also confirmed the plasmid–host trade-off [17]. This study regarding the mechanism of the fitness cost while carrying the ndm/IncX3 plasmid provides insight about how bacteria maintenance plasmids. Our findings provide a reasonable representation of the evolution of the IncX3 plasmid and provide a better understanding of its prevalence.
The present study has some limitations. We found that bamA may be a vital factor compensating for the fitness cost related to specific genetic conflicts. The role of BamA mutations should be confirmed in complementation experiments. Fitness costs and compensatory mechanisms are related to diverse metabolic pathways; thus, in-depth metabolomic and transcriptomic analyses are needed to provide more specific insights across bacterial strains and species.
In conclusion, we report that the fitness cost of the ndm/IncX3 plasmid was related to specific genetic conflicts that were ameliorated by compensatory mutations following coevolution, specifically through the bamA gene reverse mutation. Our study provides a clearer understanding of the mechanisms related to the fitness cost of carrying plasmids and their adaptive evolution with host strains, which is essential for predicting the prevalence and evolution of plasmid-mediated horizontal gene transfer.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens12010137/s1, Table S1: SNPs found in chromosomes of E. coli strains carrying phn330, Table S2: SNPs found in chromosomes of the progeny strains of E. coli C600T, Table S3: Gene functional classification and clustering analyses with DAVID of non-synonymous SNP genes, Table S4: Detailed function of non-synonymous SNP genes in E. coli C600 carrying the pHN330 plasmid, Table S5: Detailed functions of non-synonymous SNP genes in E. coli C600N500, Table S6: Detailed functions of non-synonymous SNP genes in E. coli C600Y500.

Author Contributions

Conceptualization, C.Z. (Chao Zhuo) and B.L.; data curation, F.L. and J.W.; formal analysis, F.L., J.W., Y.J. and Y.G.; funding acquisition, C.Z. (Chao Zhuo); investigation, F.L. and J.W.; methodology, F.L., B.L. and N.L.; project administration, S.X. and J.W.; resources, C.Z. (Chao Zhuo) and C.Z. (Chuyue Zhuo); software, L.Y., J.L., C.Z. and N.H.; supervision, S.X. and B.L.; validation, F.L., L.Y. and J.L.; visualization, F.L., C.Z. (Chuyue Zhuo) and N.H.; writing—original draft, F.L.; writing—review and editing, F.L., Y.J., Y.G., B.L. and C.Z. (Chao Zhuo). All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Project No. 82172318) and the Special Foundation for National Science and Technology Basic Research Program of China (Project No. 2019FY101200).

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yong, D.; Toleman, M.A.; Giske, C.G.; Cho, H.S.; Sundman, K.; Lee, K.; Walsh, T.R. Characterization of a new metallo-beta-lactamase gene, bla(NDM-1), and a novel erythromycin esterase gene carried on a unique genetic structure in Klebsiella pneumoniae sequence type 14 from India. Antimicrob. Agents Chemother. 2009, 53, 5046–5054. [Google Scholar] [CrossRef] [Green Version]
  2. Zhao, Q.; Berglund, B.; Zou, H.; Zhou, Z.; Xia, H.; Zhao, L.; Nilsson, L.E.; Li, X. Dissemination of bla(NDM-5) via IncX3 plasmids in carbapenem-resistant Enterobacteriaceae among humans and in the environment in an intensive vegetable cultivation area in eastern China. Environ. Pollut. 2021, 273, 116370. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, W.; Feng, Y.; Tang, G.; Qiao, F.; McNally, A.; Zong, Z. NDM Metallo-β-Lactamases and Their Bacterial Producers in Health Care Settings. Clin. Microbiol. Rev. 2019, 32, e00115-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. San Millan, A.; MacLean, R.C. Fitness Costs of Plasmids: A Limit to Plasmid Transmission. Microbiol. Spectr. 2017, 5, 5.5.02. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Baltrus, D.A. Exploring the costs of horizontal gene transfer. Trends Ecol. Evol. 2013, 28, 489–495. [Google Scholar] [CrossRef]
  6. Zünd, P.; Lebek, G. Generation time-prolonging R plasmids: Correlation between increases in the generation time of Escherichia coli caused by R plasmids and their molecular size. Plasmid 1980, 3, 65–69. [Google Scholar] [CrossRef]
  7. Harrison, E.; Brockhurst, M.A. Plasmid-mediated horizontal gene transfer is a coevolutionary process. Trends Microbiol. 2012, 20, 262–267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Baomo, L.; Lili, S.; Moran, R.A.; van Schaik, W.; Chao, Z. Temperature-Regulated IncX3 Plasmid Characteristics and the Role of Plasmid-Encoded H-NS in Thermoregulation. Front. Microbiol. 2021, 12, 765492. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, B.; Guo, Y.; Liu, N.; Wang, J.; Li, F.; Yao, L.; Zhuo, C. In silico Evolution and Comparative Genomic Analysis of IncX3 Plasmids Isolated From China Over Ten Years. Front. Microbiol. 2021, 12, 725391. [Google Scholar] [CrossRef]
  10. Hammer-Dedet, F.; Aujoulat, F.; Jumas-Bilak, E.; Licznar-Fajardo, P. Persistence and Dissemination Capacities of a blaNDM-5-Harboring IncX-3 Plasmid in Escherichia coli Isolated from an Urban River in Montpellier, France. Antibiotics 2022, 11, 196. [Google Scholar] [CrossRef] [PubMed]
  11. Zhang, Y.; Liao, K.; Gao, H.; Wang, Q.; Wang, X.; Li, H.; Wang, R.; Wang, H. Decreased Fitness and Virulence in ST10 Escherichia coli Harboring blaNDM-5 and mcr-1 against a ST4981 Strain with blaNDM-5. Front. Cell Infect. Microbiol. 2017, 7, 242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. He, T.; Wei, R.; Zhang, L.; Sun, L.; Pang, M.; Wang, R.; Wang, Y. Characterization of NDM-5-positive extensively resistant Escherichia coli isolates from dairy cows. Vet. Microbiol. 2017, 207, 153–158. [Google Scholar] [CrossRef]
  13. Ma, T.; Fu, J.; Xie, N.; Ma, S.; Lei, L.; Zhai, W.; Shen, Y.; Sun, C.; Wang, S.; Shen, Z.; et al. Fitness Cost of blaNDM-5-Carrying p3R-IncX3 Plasmids in Wild-Type NDM-Free Enterobacteriaceae. Microorganisms 2020, 8, 377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hall, J.P.J.; Wright, R.C.T.; Harrison, E.; Muddiman, K.J.; Wood, A.J.; Paterson, S.; Brockhurst, M.A. Plasmid fitness costs are caused by specific genetic conflicts enabling resolution by compensatory mutation. PLoS Biol. 2021, 19, e3001225. [Google Scholar] [CrossRef]
  15. Haudiquet, M.; de Sousa, J.M.; Touchon, M.; Rocha, E.P.C. Selfish, promiscuous and sometimes useful: How mobile genetic elements drive horizontal gene transfer in microbial populations. Philos. Trans. R Soc. Lond. B Biol. Sci. 2022, 377, 20210234. [Google Scholar] [CrossRef]
  16. Mc Ginty, S.; Rankin, D.J. The evolution of conflict resolution between plasmids and their bacterial hosts. Evolution 2012, 66, 1662–1670. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, J.F.; Fang, L.X.; Chang, M.X.; Cheng, M.; Zhang, H.; Long, T.F.; Li, Q.; Lian, X.L.; Sun, J.; Liao, X.P.; et al. A Trade-Off for Maintenance of Multidrug-Resistant IncHI2 Plasmids in Salmonella enterica Serovar Typhimurium through Adaptive Evolution. mSystems 2022, 7, e0024822. [Google Scholar] [CrossRef] [PubMed]
  18. Werren, J.H. Selfish genetic elements, genetic conflict, and evolutionary innovation. Proc. Natl. Acad. Sci. USA 2011, 108 (Suppl. 2), 10863–10870. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, B.; Shui, L.; Zhou, K.; Jiang, Y.; Li, X.; Guan, J.; Li, Q.; Zhuo, C. Impact of Plasmid-Encoded H-NS-like Protein on blaNDM-1-Bearing IncX3 Plasmid in Escherichia coli. J. Infect. Dis. 2020, 221, S229–S236. [Google Scholar] [CrossRef] [PubMed]
  20. CLSI. Performance standards for antimicrobial susceptibility testing. CLSI Supplement M100; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2021. [Google Scholar]
  21. Yang, S.; Lopez, C.R.; Zechiedrich, E.L. Quorum sensing and multidrug transporters in Escherichia coli. Proc. Natl. Acad. Sci. USA 2006, 103, 2386–2391. [Google Scholar] [CrossRef]
  22. Liu, D.; Liu, Z.S.; Hu, P.; Cai, L.; Fu, B.Q.; Li, Y.S.; Lu, S.Y.; Liu, N.N.; Ma, X.L.; Chi, D.; et al. Characterization of surface antigen protein 1 (SurA1) from Acinetobacter baumannii and its role in virulence and fitness. Vet. Microbiol. 2016, 186, 126–138. [Google Scholar] [CrossRef]
  23. Trinetta, V.; Magossi, G.; Allard, M.W.; Tallent, S.M.; Brown, E.W.; Lomonaco, S. Characterization of Salmonella enterica Isolates from Selected U.S. Swine Feed Mills by Whole-Genome Sequencing. Foodborne Pathog. Dis. 2020, 17, 126–136. [Google Scholar] [CrossRef]
  24. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Olawoye, I.B.; Frost, S.D.W.; Happi, C.T. The Bacteria Genome Pipeline (BAGEP): An automated, scalable workflow for bacteria genomes with Snakemake. PeerJ 2020, 8, e10121. [Google Scholar] [CrossRef] [PubMed]
  26. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef]
  27. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Buckner, M.M.C.; Saw, H.T.H.; Osagie, R.N.; McNally, A.; Ricci, V.; Wand, M.E.; Woodford, N.; Ivens, A.; Webber, M.A.; Piddock, L.J.V. Clinically Relevant Plasmid-Host Interactions Indicate that Transcriptional and Not Genomic Modifications Ameliorate Fitness Costs of Klebsiella pneumoniae Carbapenemase-Carrying Plasmids. mBio 2018, 9, e02303-17. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, J.; Guo, Z.W.; Zhi, C.P.; Yang, T.; Zhao, J.J.; Chen, X.J.; Zeng, L.; Lv, L.C.; Zeng, Z.L.; Liu, J.H. Impact of plasmid-borne oqxAB on the development of fluoroquinolone resistance and bacterial fitness in Escherichia coli. J. Antimicrob. Chemother. 2017, 72, 1293–1302. [Google Scholar] [CrossRef] [Green Version]
  30. Murarka, A.; Clomburg, J.M.; Moran, S.; Shanks, J.V.; Gonzalez, R. Metabolic analysis of wild-type Escherichia coli and a pyruvate dehydrogenase complex (PDHC)-deficient derivative reveals the role of PDHC in the fermentative metabolism of glucose. J. Biol. Chem. 2010, 285, 31548–31558. [Google Scholar] [CrossRef] [Green Version]
  31. Schmitt, R. Analysis of melibiose mutants deficient in alpha-galactosidase and thiomethylgalactoside permease II in Escherichia coli K-12. J. Bacteriol. 1968, 96, 462–471. [Google Scholar] [CrossRef]
  32. Lengeler, J. Mutations affecting transport of the hexitols D-mannitol, D-glucitol, and galactitol in Escherichia coli K-12: Isolation and mapping. J. Bacteriol. 1975, 124, 26–38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Denger, K.; Weiss, M.; Felux, A.K.; Schneider, A.; Mayer, C.; Spiteller, D.; Huhn, T.; Cook, A.M.; Schleheck, D. Sulphoglycolysis in Escherichia coli K-12 closes a gap in the biogeochemical sulphur cycle. Nature 2014, 507, 114–117. [Google Scholar] [CrossRef] [Green Version]
  34. Voulhoux, R.; Bos, M.P.; Geurtsen, J.; Mols, M.; Tommassen, J. Role of a highly conserved bacterial protein in outer membrane protein assembly. Science 2003, 299, 262–265. [Google Scholar] [CrossRef]
  35. Rodríguez-Alonso, R.; Létoquart, J.; Nguyen, V.S.; Louis, G.; Calabrese, A.N.; Iorga, B.I.; Radford, S.E.; Cho, S.H.; Remaut, H.; Collet, J.F. Structural insight into the formation of lipoprotein-β-barrel complexes. Nat. Chem. Biol. 2020, 16, 1019–1025. [Google Scholar] [CrossRef]
  36. Fiedler, W.; Rotering, H. Properties of Escherichia coli mutants lacking membrane-derived oligosaccharides. J. Biol. Chem. 1988, 263, 14684–14689. [Google Scholar] [CrossRef]
  37. Bossemeyer, D.; Borchard, A.; Dosch, D.C.; Helmer, G.C.; Epstein, W.; Booth, I.R.; Bakker, E.P. K+-transport protein TrkA of Escherichia coli is a peripheral membrane protein that requires other trk gene products for attachment to the cytoplasmic membrane. J. Biol. Chem. 1989, 264, 16403–16410. [Google Scholar] [CrossRef] [PubMed]
  38. San Millan, A.; Toll-Riera, M.; Qi, Q.; Betts, A.; Hopkinson, R.J.; McCullagh, J.; MacLean, R.C. Integrative analysis of fitness and metabolic effects of plasmids in Pseudomonas aeruginosa PAO1. ISME J. 2018, 12, 3014–3024. [Google Scholar] [CrossRef] [PubMed]
  39. Subbiah, M.; Top, E.M.; Shah, D.H.; Call, D.R. Selection pressure required for long-term persistence of blaCMY-2-positive IncA/C plasmids. Appl. Environ. Microbiol. 2011, 77, 4486–4493. [Google Scholar] [CrossRef] [Green Version]
  40. Lang, K.S.; Johnson, T.J. Transcriptome modulations due to A/C2 plasmid acquisition. Plasmid 2015, 80, 83–89. [Google Scholar] [CrossRef] [Green Version]
  41. Loftie-Eaton, W.; Bashford, K.; Quinn, H.; Dong, K.; Millstein, J.; Hunter, S.; Thomason, M.K.; Merrikh, H.; Ponciano, J.M.; Top, E.M. Compensatory mutations improve general permissiveness to antibiotic resistance plasmids. Nat. Ecol. Evol. 2017, 1, 1354–1363. [Google Scholar] [CrossRef]
  42. Noinaj, N.; Gumbart, J.C.; Buchanan, S.K. The beta-barrel assembly machinery in motion. Nat. Rev. Microbiol. 2017, 15, 197–204. [Google Scholar] [CrossRef] [Green Version]
  43. Li, Y.; Zhu, X.; Zhang, J.; Lin, Y.; You, X.; Chen, M.; Wang, Y.; Zhu, N.; Si, S. Identification of a Compound That Inhibits the Growth of Gram-Negative Bacteria by Blocking BamA-BamD Interaction. Front. Microbiol. 2020, 11, 1252. [Google Scholar] [CrossRef]
  44. Storek, K.M.; Auerbach, M.R.; Shi, H.; Garcia, N.K.; Sun, D.; Nickerson, N.N.; Vij, R.; Lin, Z.; Chiang, N.; Schneider, K.; et al. Monoclonal antibody targeting the β-barrel assembly machine of Escherichia coli is bactericidal. Proc. Natl. Acad. Sci. USA 2018, 115, 3692–3697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Noinaj, N.; Kuszak, A.J.; Gumbart, J.C.; Lukacik, P.; Chang, H.; Easley, N.C.; Lithgow, T.; Buchanan, S.K. Structural insight into the biogenesis of β-barrel membrane proteins. Nature 2013, 501, 385–390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Wu, T.; Malinverni, J.; Ruiz, N.; Kim, S.; Silhavy, T.J.; Kahne, D. Identification of a multicomponent complex required for outer membrane biogenesis in Escherichia coli. Cell 2005, 121, 235–245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Gawel, D.; Hamilton, M.D.; Schaaper, R.M. A novel mutator of Escherichia coli carrying a defect in the dgt gene, encoding a dGTP triphosphohydrolase. J. Bacteriol. 2008, 190, 6931–6939. [Google Scholar] [CrossRef] [Green Version]
  48. Thongsomboon, W.; Serra, D.O.; Possling, A.; Hadjineophytou, C.; Hengge, R.; Cegelski, L. Phosphoethanolamine cellulose: A naturally produced chemically modified cellulose. Science 2018, 359, 334–338. [Google Scholar] [CrossRef] [Green Version]
  49. Sherman, D.J.; Xie, R.; Taylor, R.J.; George, A.H.; Okuda, S.; Foster, P.J.; Needleman, D.J.; Kahne, D. Lipopolysaccharide is transported to the cell surface by a membrane-to-membrane protein bridge. Science 2018, 359, 798–801. [Google Scholar] [CrossRef] [Green Version]
  50. Brabetz, W.; Müller-Loennies, S.; Holst, O.; Brade, H. Deletion of the heptosyltransferase genes rfaC and rfaF in Escherichia coli K-12 results in an Re-type lipopolysaccharide with a high degree of 2-aminoethanol phosphate substitution. Eur. J. Biochem. 1997, 247, 716–724. [Google Scholar] [CrossRef]
  51. Cote, J.M.; Taylor, E.A. The Glycosyltransferases of LPS Core: A Review of Four Heptosyltransferase Enzymes in Context. Int. J. Mol. Sci. 2017, 18, 2256. [Google Scholar] [CrossRef]
  52. Kruse, T.; Bork-Jensen, J.; Gerdes, K. The morphogenetic MreBCD proteins of Escherichia coli form an essential membrane-bound complex. Mol. Microbiol. 2005, 55, 78–89. [Google Scholar] [CrossRef] [PubMed]
  53. Eger, E.; Schwabe, M.; Schulig, L.; Hübner, N.O.; Bohnert, J.A.; Bornscheuer, U.T.; Heiden, S.E.; Müller, J.U.; Adnan, F.; Becker, K.; et al. Extensively Drug-Resistant Klebsiella pneumoniae Counteracts Fitness and Virulence Costs That Accompanied Ceftazidime-Avibactam Resistance Acquisition. Microbiol. Spectr. 2022, 10, e0014822. [Google Scholar] [CrossRef] [PubMed]
  54. Harrison, E.; Guymer, D.; Spiers, A.J.; Paterson, S.; Brockhurst, M.A. Parallel compensatory evolution stabilizes plasmids across the parasitism-mutualism continuum. Curr. Biol. 2015, 25, 2034–2039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Cheng, X.; de Bruijn, I.; van der Voort, M.; Loper, J.E.; Raaijmakers, J.M. The Gac regulon of Pseudomonas fluorescens SBW25. Environ. Microbiol. Rep. 2013, 5, 608–619. [Google Scholar] [CrossRef] [PubMed]
  56. Ma, K.; Feng, Y.; Zong, Z. Fitness cost of a mcr-1-carrying IncHI2 plasmid. PLoS One 2018, 13, e0209706. [Google Scholar] [CrossRef] [Green Version]
  57. Rangarajan, E.S.; Li, Y.; Ajamian, E.; Iannuzzi, P.; Kernaghan, S.D.; Fraser, M.E.; Cygler, M.; Matte, A. Crystallographic trapping of the glutamyl-CoA thioester intermediate of family I CoA transferases. J. Biol. Chem. 2005, 280, 42919–42928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Mironov, A.; Seregina, T.; Shatalin, K.; Nagornykh, M.; Shakulov, R.; Nudler, E. CydDC functions as a cytoplasmic cystine reductase to sensitize Escherichia coli to oxidative stress and aminoglycosides. Proc. Natl. Acad. Sci. USA 2020, 117, 23565–23570. [Google Scholar] [CrossRef]
  59. Ohtsu, I.; Wiriyathanawudhiwong, N.; Morigasaki, S.; Nakatani, T.; Kadokura, H.; Takagi, H. The L-cysteine/L-cystine shuttle system provides reducing equivalents to the periplasm in Escherichia coli. J. Biol. Chem. 2010, 285, 17479–17487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Clugston, S.L.; Barnard, J.F.; Kinach, R.; Miedema, D.; Ruman, R.; Daub, E.; Honek, J.F. Overproduction and characterization of a dimeric non-zinc glyoxalase I from Escherichia coli: Evidence for optimal activation by nickel ions. Biochemistry 1998, 37, 8754–8763. [Google Scholar] [CrossRef] [PubMed]
  61. MacLean, M.J.; Ness, L.S.; Ferguson, G.P.; Booth, I.R. The role of glyoxalase I in the detoxification of methylglyoxal and in the activation of the KefB K+ efflux system in Escherichia coli. Mol. Microbiol. 1998, 27, 563–571. [Google Scholar] [CrossRef] [PubMed]
  62. Wang, Y.; Bao, W.; Guo, N.; Chen, H.; Cheng, W.; Jin, K.; Shen, F.; Xu, J.; Zhang, Q.; Wang, C.; et al. Antimicrobial activity of the imipenem/rifampicin combination against clinical isolates of Acinetobacter baumannii grown in planktonic and biofilm cultures. World J. Microbiol. Biotechnol. 2014, 30, 3015–3025. [Google Scholar] [CrossRef] [PubMed]
  63. Kohanski, M.A.; Dwyer, D.J.; Hayete, B.; Lawrence, C.A.; Collins, J.J. A common mechanism of cellular death induced by bactericidal antibiotics. Cell 2007, 130, 797–810. [Google Scholar] [CrossRef] [PubMed]
Pathogens 12 00137 g001 550
Figure 1. Growth curves of host strains and transformants. (A) Comparison of the growth curves for J53 and J53T showed no significant changes. (B) Comparison of growth curves for C600 and C600T; the growth rate was significantly decreased in C600T.
Figure 1. Growth curves of host strains and transformants. (A) Comparison of the growth curves for J53 and J53T showed no significant changes. (B) Comparison of growth curves for C600 and C600T; the growth rate was significantly decreased in C600T.
Pathogens 12 00137 g001
Pathogens 12 00137 g002 550
Figure 2. The competitive indices of the C600 and J53 strains. The transformants and host strains were paired to perform competition experiments and the competitive index (CI) was calculated. CI values < 1 indicate that the carried plasmid generates a fitness cost, with smaller CI values indicating higher fitness costs.
Figure 2. The competitive indices of the C600 and J53 strains. The transformants and host strains were paired to perform competition experiments and the competitive index (CI) was calculated. CI values < 1 indicate that the carried plasmid generates a fitness cost, with smaller CI values indicating higher fitness costs.
Pathogens 12 00137 g002
Pathogens 12 00137 g003 550
Figure 3. Experimental evolution of C600T without antibiotics or with imipenem (2 mg/L). (A,C) Comparison of the growth curves for C600N500, C600Y500, and C600T. (B,D) After evolution without or with antibiotics, the competitive index of C600N500 and C600Y500 increased significantly compared to that of C600T. * p = 0.003; ** p = 0.001.
Figure 3. Experimental evolution of C600T without antibiotics or with imipenem (2 mg/L). (A,C) Comparison of the growth curves for C600N500, C600Y500, and C600T. (B,D) After evolution without or with antibiotics, the competitive index of C600N500 and C600Y500 increased significantly compared to that of C600T. * p = 0.003; ** p = 0.001.
Pathogens 12 00137 g003
Pathogens 12 00137 g004 550
Figure 4. SNP gene function classification. SNP: single nucleotide polymorphism. GO clustering analysis was performed on non-synonymous SNP genes, which could be divided into metabolism- and cell-membrane-related genes. Numbers above the bars represent the number of non-synonymous SNP genes.
Figure 4. SNP gene function classification. SNP: single nucleotide polymorphism. GO clustering analysis was performed on non-synonymous SNP genes, which could be divided into metabolism- and cell-membrane-related genes. Numbers above the bars represent the number of non-synonymous SNP genes.
Pathogens 12 00137 g004
Table
Table 1. Strain descriptions.
Table 1. Strain descriptions.
StrainPlasmid RepliconDescription
C600Nonemodel E. coli strain
C600TIncX3C600 transformant
C600N500IncX3500th passage of C600T without antibiotics
C600Y500IncX3500th passage of C600T with 2 mg/L imipenem
C600D500None500th passage of C600
J53Nonemodel E. coli strain
J53TIncX3J53 transformant
Table
Table 2. Resistance profiles of the host strains and transformants.
Table 2. Resistance profiles of the host strains and transformants.
StrainMIC (mg/L)
MEMIPMFEPCAZATMAMKLEV
C6000.1250.250.12510.250.1250.125
C600T64
(512-fold) a		64
(256-fold) a		32
(256-fold) a		128
(128-fold) a		8
(32-fold) a		0.1250.125
C600N50032323212880.1250.125
C600Y50016323212880.1250.125
J530.1250.1250.1250.1250.540.125
J53T4
(32-fold) a		16
(128-fold) a		64
(512-fold) a		64
(512-fold) a		128
(512-fold) a		40.125
MEM—meropenem, IPM—imipenem, FEP—cefepime, CAZ—ceftazidime, ATM—aztreonam, AMK– amikacin, LEV—Levofloxacin. a Fold increase in MIC values of the transformants compared to the model strain.
Table
Table 3. Number of SNPs in the host strains carrying the pHN330 plasmid.
Table 3. Number of SNPs in the host strains carrying the pHN330 plasmid.
Strain SNP Counts
Total SNPsCDS CountIntergenic Sequence Count
J53T808
C600T291118
SNP—single nucleotide polymorphism, CDS—coding sequence. SNPs differed between J53 and C600, with J53 having no SNP in the CDS region and C600 having 11 SNPs in the CDS region. SNP gene—gene occurring in the CDS region.
Table
Table 4. Non-synonymous SNP genes in host strains carrying the pHN330 plasmid.
Table 4. Non-synonymous SNP genes in host strains carrying the pHN330 plasmid.
StrainNon-Synonymous SNP Gene
J53TNone
C600TptrA, lpdA, bamA, melA, gatA, mdoB, yihV, mntH, trkA, and ydfK
Table
Table 5. SNP counts in progeny strains.
Table 5. SNP counts in progeny strains.
Strain SNP Counts
Total SNPsCDS CountIntergenic Sequence Count
C600N500281711
C600Y500271116
SNP: single nucleotide polymorphism, CDS: coding sequence. After adaptive evolution, the SNP was located in the CDS region, and the intergenic sequence region occurred in the progeny strains.
Table
Table 6. Non-synonymous SNP genes in progeny strains.
Table 6. Non-synonymous SNP genes in progeny strains.
StrainNon-Synonymous SNP Gene
C600N500gudD, araC, dgt, bamA, rfaF, lptF, ycjX, bcsG, pepQ, mshA, sfmD, paaJ, mreC, and tfaE
C600Y500gloA, ydiF, fliY, hycF, caiA, bamA, rlmB, sbp, acnA, and yghQ
SNP gene—gene occurring in the CDS region.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, F.; Wang, J.; Jiang, Y.; Guo, Y.; Liu, N.; Xiao, S.; Yao, L.; Li, J.; Zhuo, C.; He, N.; et al. Adaptive Evolution Compensated for the Plasmid Fitness Costs Brought by Specific Genetic Conflicts. Pathogens 2023, 12, 137. https://doi.org/10.3390/pathogens12010137

AMA Style

Li F, Wang J, Jiang Y, Guo Y, Liu N, Xiao S, Yao L, Li J, Zhuo C, He N, et al. Adaptive Evolution Compensated for the Plasmid Fitness Costs Brought by Specific Genetic Conflicts. Pathogens. 2023; 12(1):137. https://doi.org/10.3390/pathogens12010137

Chicago/Turabian Style

Li, Feifeng, Jiong Wang, Ying Jiang, Yingyi Guo, Ningjing Liu, Shunian Xiao, Likang Yao, Jiahui Li, Chuyue Zhuo, Nanhao He, and et al. 2023. "Adaptive Evolution Compensated for the Plasmid Fitness Costs Brought by Specific Genetic Conflicts" Pathogens 12, no. 1: 137. https://doi.org/10.3390/pathogens12010137

APA Style

Li, F., Wang, J., Jiang, Y., Guo, Y., Liu, N., Xiao, S., Yao, L., Li, J., Zhuo, C., He, N., Liu, B., & Zhuo, C. (2023). Adaptive Evolution Compensated for the Plasmid Fitness Costs Brought by Specific Genetic Conflicts. Pathogens, 12(1), 137. https://doi.org/10.3390/pathogens12010137

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop