Search

Genes that undergo horizontal gene transfer (HGT) evolve in different genomic backgrounds as they move between hosts, in contrast to genes that evolve under strict vertical inheritance. Despite the ubiquity of HGT in microbial communities, the effects of host-switching on gene evolution have been understudied. Here, we present a novel framework to examine the consequences of host-switching on gene evolution by probing the existence and form of host-dependent mutational effects. We started exploring the effects of HGT on gene evolution by focusing on an antibiotic resistance gene (encoding a beta-lactamase) commonly found on conjugative plasmids in Enterobacteriaceae pathogens. By reconstructing the resistance landscape for a small set of mutationally connected alleles in three enteric species (Escherichia coli, Salmonella enterica, and Klebsiella pneumoniae), we uncovered that the landscape topographies were largely aligned with very low levels of host-dependent mutational effects. By simulating gene evolution with and without HGT using the species-specific empirical landscapes, we found that evolutionary outcomes were similar despite HGT. These findings suggest that the adaptive evolution of a mobile gene in one species can translate to adaptation in another species. In such a case, vehicles of cross-species HGT such as plasmids enable a distributed form of genetic evolution across a bacterial community, where species can ‘crowdsource’ adaptation from other community members. The role of evolutionary crowdsourcing on the evolution of bacteria merits further investigation.

Vibrio cholerae biotype El Tor is perpetuating the longest cholera pandemic in recorded history. The genomic islands VSP-1 and VSP-2 distinguish El Tor from previous pandemic V. cholerae strains. Using a co-occurrence analysis of VSP genes in >200,000 bacterial genomes we built gene networks to infer biological functions encoded in these islands. This revealed that dncV, a component of the cyclic-oligonucleotide-based anti-phage signalling system (CBASS) anti-phage defence system, co-occurs with an uncharacterized gene vc0175 that we rename avcD for anti-viral cytodine deaminase. We show that AvcD is a deoxycytidylate deaminase and that its activity is post-translationally inhibited by a non-coding RNA named AvcI. AvcID and bacterial homologues protect bacterial populations against phage invasion by depleting free deoxycytidine nucleotides during infection, thereby decreasing phage replication. Homologues of avcD exist in all three domains of life, and bacterial AvcID defends against phage infection by combining traits of two eukaryotic innate viral immunity proteins, APOBEC and SAMHD1.

Plasmids are mobile genetic elements that can act as mutualists or parasites to their bacterial hosts depending on their accessory genes and environment. Ecological network theory predicts that mutualists, such as plasmids with antimicrobial resistance (AMR) genes in the presence of antimicrobials, should act as generalists, while plasmids without beneficial genes are expected to be more specialised. Therefore, whether the relationship between plasmid and host is mutualistic or antagonistic is likely to have a strong impact on the formation of interaction network structures and the spread of AMR genes across microbial networks. Here we re-analyse Hi-C metagenome data from wastewater samples and identify plasmid signatures with machine learning to generate a natural host-plasmid network. We found that AMR-carrying plasmids indeed interacted with more hosts than non-AMR plasmids (on average 14 versus 3, respectively). The AMR plasmid-host subnetwork showed a much higher connectedness and nestedness than the subnetwork associated with non-AMR plasmids. The overall network was clustered around Proteobacteria and AMR-carrying plasmids giving them a crucial role in network connectivity. Therefore, by forming mutualistic networks with their hosts, beneficial AMR plasmids lead to more connected network structures that ultimately share a larger gene pool of AMR genes across the network.

Although animal manure is applied to agricultural fields for its nutrient value, it may also contain potential contaminants. To determine the variability in such contaminants as well as in valuable nutrients, nine uncomposted manure samples from Idaho dairies collected during 2.5 years were analyzed for macro- and micro-nutrients, hormones, phytoestrogens, antibiotics, veterinary drugs, antibiotic resistance genes, and genetic elements involved in the spread of antibiotic resistance. Total N ranged from 6.8 to 30.7 (C:N of 10 to 21), P from 2.4 to 9.0, and K from 10.2 to 47.7 g/kg manure. Zn (103 – 348 mg/kg) was more abundant than Cu (56 – 127 mg/kg) in all samples. Phytoestrogens were the most prevalent contaminants detected, with concentrations fluctuating over time, reflecting animal diets. This is the first study to document the presence of flunixin, a non-steroidal anti-inflammatory drug, in solid stacked manure from regular dairy operations. Monensin was the most frequently detected antibiotic. Progesterones and sulfonamides were regularly detected. We also investigated the relative abundance of several types of plasmids involved in the spread of antibiotic resistance in clinical settings. Plasmids belonging to the IncI, IncP, and IncQ1 incompatibility groups were found in almost all manure samples. IncQ1 plasmids, class 1 integrons, and sulfonamide resistance genes were the most widespread and abundant genetic element surveyed, emphasizing their potential role in the spread of antibiotic resistance. The benefits associated with amending agricultural soils with dairy manure must be carefully weighed against the potential negative consequences of any manure contaminants.

Promiscuous plasmids like IncP-1 plasmids play an important role in the bacterial adaptation to pollution by acquiring and distributing xenobiotic catabolic genes. However, most information comes from isolates and the role of plasmids in governing community-wide bacterial adaptation to xenobiotics and other adaptive forces is not fully understood. Current information on the contribution of IncP-1 plasmids in community adaptation is limited because methods are lacking that directly isolate and identify the plasmid borne adaptive functions in whole-community DNA. In this study, we optimized long-range PCR to directly access and identify the cargo carried by IncP-1 plasmids in environmental DNA. The DNA between the IncP-1 backbone genes trbP and traC, a main insertion site of adaptive trait determinants, is amplified and its content analyzed by high-throughput sequencing. The method was applied to DNA of an on-farm biopurification system (BPS), treating pesticide contaminated wastewater, to examine whether horizontal gene exchange of catabolic functions by IncP-1 plasmids is a main driver of community adaptation in BPS. The cargo recovered from BPS community DNA encoded catabolic but also resistance traits and various other (un)known functions. Unexpectedly, genes with catabolic traits composed only a minor fraction of the cargo, indicating that the IncP-1 region between trbP and traC is not a major contributor to catabolic adaptation of the BPS microbiome. Instead, it contains a functionally diverse set of genes which either may assist biodegradation functions, be remnants of random gene recruitment, or confer other crucial functions for proliferation in the BPS environment. IMPORTANCE This study presents a long-range PCR for direct and cultivation-independent access to the identity of the cargo of a major insertion hot spot of adaptive genes in IncP-1 plasmids and hence a new mobilome tool for understanding the role of IncP-1 plasmids in complex communities. The method was applied to DNA of an on-farm biopurification system (BPS) treating pesticide-contaminated wastewater, aiming at new insights on whether horizontal exchange of catabolic functions by IncP-1 plasmids is a main driver of community adaptation in BPS. Unexpectedly, catabolic functions represented a small fraction of the cargo genes while multiple other gene functions were recovered. These results show that the cargo of the target insertion hot spot in IncP-1 plasmids in a community, not necessarily relates to the main obvious selective trait imposed on that community. Instead, these functions might contribute to adaptation to unknown selective forces or represent remnants of random gene recruitment.

Multidrug resistance (MDR) of pathogens is an ongoing public health crisis exacerbated by the horizontal transfer of antibiotic resistance genes via conjugative plasmids. Factors that stabilize these plasmids in bacterial communities contribute to an even higher incidence of MDR, given the increased likelihood that a host will already contain a plasmid when it acquires another through conjugation. Here, we show one such stabilizing factor is host-plasmid coevolution under antibiotic selection, which facilitated the emergence of MDR via two distinct plasmids in communities consisting of Escherichia coli and Klebsiella pneumoniae once antibiotics were removed. In our system, evolution promoted greater stability of a plasmid in its coevolved host. Further, pleiotropic effects resulted in greater plasmid persistence in both novel host-plasmid combinations and, in some cases, multi-plasmid hosts. This evolved stability favoured the generation of MDR cells and thwarted their loss within communities with multiple plasmids. By selecting for plasmid persistence, the application of antibiotics may promote MDR well after their original period of use.

Bacterial plasmids substantially contribute to the rapid spread of antibiotic resistance, which is a crisis in healthcare today. Coevolution of plasmids and their hosts promotes this spread of resistance by ameliorating the cost of plasmid carriage. However, our knowledge of plasmid–bacteria coevolution is solely based on studies done in well-mixed liquid cultures, even though biofilms represent the main way of bacterial life on Earth and are responsible for most infections. The spatial structure and the heterogeneity provided by biofilms are known to lead to increased genetic diversity as compared with well-mixed liquids. Therefore, we expect that growth in this complex environment could affect the evolutionary trajectories of plasmid–host dyads. We experimentally evolved Shewanella oneidensis MR-1 with plasmid pBP136Gm in biofilms and chemostats and sequenced the genomes of clones and populations. Biofilm populations not only maintained a higher diversity of mutations than chemostat populations but contained a few clones with markedly more persistent plasmids that evolved via multiple distinct trajectories. These included the acquisition of a putative toxin–antitoxin transposon by the plasmid and chromosomal mutations. Some of these genetic changes resulted in loss of plasmid transferability or decrease in plasmid cost. Growth in chemostats led to a higher proportion of variants with decreased plasmid persistence, a phenomenon not detected in biofilms. We suggest that the presence of more stable plasmid–host dyads in biofilms reflects higher genetic diversity and possibly unknown selection pressures. Overall, this study underscores the importance of the mode of growth in the evolution of antibiotic-resistant bacteria.

The rapid spread of antibiotic resistance among bacterial pathogens is a serious human health threat. While a range of environments have been identified as reservoirs of antibiotic resistance genes (ARGs), we lack understanding of the origins of these ARGs and their spread from environment to clinic. This is partly due to our inability to identify the natural bacterial hosts of ARGs and the mobile genetic elements that mediate this spread, such as plasmids and integrons. Here we demonstrate that the in vivo proximity-ligation method Hi-C can reconstruct a known plasmid-host association from a wastewater community, and identify the in situ host range of ARGs, plasmids, and integrons by physically linking them to their host chromosomes. Hi-C detected both previously known and novel associations between ARGs, mobile genetic elements and host genomes, thus validating this method. We showed that IncQ plasmids and class 1 integrons had the broadest host range in this wastewater, and identified bacteria belonging to Moraxellaceae, Bacteroides, and Prevotella, and especially Aeromonadaceae as the most likely reservoirs of ARGs in this community. A better identification of the natural carriers of ARGs will aid the development of strategies to limit resistance spread to pathogens.

To increase our basic understanding of the ecology and evolution of conjugative plasmids, we need reliable estimates of their rate of transfer between bacterial cells. Current assays to measure transfer rate are based on deterministic modeling frameworks. However, some cell numbers in these assays can be very small, making estimates that rely on these numbers prone to noise. Here, we take a different approach to estimate plasmid transfer rate, which explicitly embraces this noise. Inspired by the classic fluctuation analysis of Luria and Delbrück, our method is grounded in a stochastic modeling framework. In addition to capturing the random nature of plasmid conjugation, our new methodology, the Luria–Delbrück method (“LDM”), can be used on a diverse set of bacterial systems, including cases for which current approaches are inaccurate. A notable example involves plasmid transfer between different strains or species where the rate that one type of cell donates the plasmid is not equal to the rate at which the other cell type donates. Asymmetry in these rates has the potential to bias or constrain current transfer estimates, thereby limiting our capabilities for estimating transfer in microbial communities. In contrast, the LDM overcomes obstacles of traditional methods by avoiding restrictive assumptions about growth and transfer rates for each population within the assay. Using stochastic simulations and experiments, we show that the LDM has high accuracy and precision for estimation of transfer rates compared to the most widely used methods, which can produce estimates that differ from the LDM estimate by orders of magnitude.

SUMMARYThe El Tor biotype of Vibrio cholerae is responsible for perpetuating the longest cholera pandemic in recorded history (1961-current). The genomic islands VSP-1 and -2 are two understudied genetic features that distinguish El Tor from previous pandemics. To understand their utility, we calculated the co-occurrence of VSP genes across bacterial genomes. This analysis predicted the previously uncharacterized vc0175, herein renamed deoxycytidylate deaminase Vibrio (dcdV), is in a gene network with dncV, a cyclic GMP-AMP synthase involved in phage defense. DcdV consists of two domains, a P-loop kinase and a deoxycytidylate deaminase, that are required for the deamination of dCTP and dCMP, inhibiting phage predation by corrupting cellular nucleotide concentrations. Additionally, DcdV is post-translationally inhibited by a unique noncoding RNA encoded 5’ of the dcdV locus. DcdV homologs are conserved in bacteria and eukaryotes and our results identify V. cholerae DcdV as the founding member of a previously undescribed bacterial phage defense system.

To increase our basic understanding of the ecology and evolution of conjugative plasmids, we need a reliable estimate of their rate of transfer between bacterial cells. However, accurate estimates of plasmid transfer have remained elusive due to biological and experimental complexity. Current methods to measure transfer rate can be confounded by many factors. A notable example involves plasmid transfer between different strains or species where the rate that one type of cell donates the plasmid is not equal to the rate at which the other cell type donates. Asymmetry in these rates has the potential to bias or constrain current transfer estimates, thereby limiting our capabilities for estimating transfer in microbial communities. Inspired by the classic fluctuation analysis of Luria and Delbrück, we develop a novel approach, the Luria-Delbrück method (‘LDM’), for estimating plasmid transfer rate. Our new approach embraces the stochasticity of conjugation departing from the current deterministic population dynamic methods. In addition, the LDM overcomes obstacles of traditional methods by not being affected by different growth and transfer rates for each population within the assay. Using stochastic simulations and experiments, we show that the LDM has high accuracy and precision for estimation of transfer rates compared to the most widely used methods, which can produce estimates that differ from the LDM estimate by orders of magnitude.Significance StatementConjugative plasmids play significant roles in the ecology and evolution of microbial communities. Notably, antibiotic resistance genes are often encoded on conjugative plasmids. Thus, conjugation—the transfer of a plasmid copy from one cell to another—is a common way for antibiotic resistance to spread between important clinical pathogens. For both public health modeling and a basic understanding of microbial population biology, accurate estimates of this fundamental rate are of great consequence. We show that widely used methods can lead to biased estimates, deviating from true values by several orders of magnitude. Therefore, we developed a new approach, inspired by the classic fluctuation analysis of Luria and Delbrück, for accurately assessing the rate of plasmid conjugation under a variety of conditions.

By characterizing the trajectories of antibiotic resistance gene transfer in bacterial communities such as the gut microbiome, we will better understand the factors that influence this spread of resistance. Our aim was to investigate the host network of a multi-drug resistance broad-host-range plasmid in the culturable gut microbiome of zebrafish. This was done throughin vitroandin vivoconjugation experiments withEscherichia colias donor of the plasmid pB10::gfp. When this donor was mixed with the extracted gut microbiome, only transconjugants ofAeromonas veroniiwere detected. In separate matings between the same donor and four prominent isolates from the gut microbiome, the plasmid transferred to two of these four isolates,A. veronii and Plesiomonas shigelloides,but not toShewanella putrefaciensandVibrio mimicus. When theseA. veronii and P. shigelloidestransconjugants were the donors in matings with the same four isolates, the plasmid now also transferred fromA. veroniitoS. putrefaciens.P. shigelloideswas unable to donate the plasmid andV. mimicuswas unable to acquire it. Finally, when theE. colidonor was addedin vivoto zebrafish through their food, plasmid transfer was observed in the gut but only toAchromobacter sp.,a rare member of the gut microbiome. This work shows that the success of plasmid-mediated antibiotic resistance spread in a gut microbiome depends on the donor-recipient species combinations and therefore their spatial arrangement. It also suggests that rare gut microbiome members should not be ignored as potential reservoirs of multi-drug resistance plasmids from food.ImportanceTo understand how antibiotic resistance plasmids end up in human pathogens it is crucial to learn how, where and when they are transferred and maintained in members of bacterial communities such as the gut microbiome. To gain insight into the network of plasmid-mediated antibiotic resistance sharing in the gut microbiome, we investigated the transferability and maintenance of a multi-drug resistance plasmid among the culturable bacteria of the zebrafish gut. We show that the success of plasmid-mediated antibiotic resistance spread in a gut microbiome can depend on which species are involved, as some are important nodes in the plasmid-host network and others dead-ends. Our findings also suggest that rare gut microbiome members should not be ignored as potential reservoirs of multi-drug resistance plasmids from food.

Antibiotic resistance (AR) is a threat to modern medicine, and plasmids are driving the global spread of AR by horizontal gene transfer across microbiomes and environments. Determining the mobile resistome responsible for this spread of AR among environments is essential in our efforts to attenuate the current crisis. Biosolids are a wastewater treatment plant (WWTP) byproduct used globally as fertilizer in agriculture. Here, we investigated the mobile resistome of biosolids that are used as fertilizer. This was done by capturing resistance plasmids that can transfer to human pathogens and commensal bacteria. We used a higher-throughput version of the exogenous plasmid isolation approach by mixing several ESKAPE pathogens and a commensal Escherichia coli with biosolids and screening for newly acquired resistance to about 10 antibiotics in these strains. Six unique resistance plasmids transferred to Salmonella typhimurium, Klebsiella aerogenes, and E. coli. All the plasmids were self-transferable and carried 3–6 antibiotic resistance genes (ARG) conferring resistance to 2–4 antibiotic classes. These plasmids-borne resistance genes were further embedded in genetic elements promoting intracellular recombination (i.e., transposons or class 1 integrons). The plasmids belonged to the broad-host-range plasmid (BHR) groups IncP-1 or PromA. Several of them were persistent in their new hosts when grown in the absence of antibiotics, suggesting that the newly acquired drug resistance traits would be sustained over time. This study highlights the role of BHRs in the spread of ARG between environmental bacteria and human pathogens and commensals, where they may persist. The work further emphasizes biosolids as potential vehicles of highly mobile plasmid-borne antibiotic resistance genes.

IS1071, an insertion element that primarily flanks organic xenobiotic degradation genes in cultured isolates, is suggested to play a key role in the formation and distribution of bacterial catabolic pathway gene clusters. However, in environmental settings, the identity of the IS1071 genetic cargo and its correspondence to the local selective conditions remain unknown. To respond, we developed a long-range PCR approach amplifying accessory genes between two IS1071 copies from community DNA followed by amplicon sequencing. We applied this method to pesticide-exposed environments, i.e. linuron-treated agricultural soil and on-farm biopurification systems (BPS) treating complex agricultural wastewater, as to non-treated controls. Amplicons were mainly recovered from the pesticide-exposed environments and the BPS matrix showed a higher size diversity compared to the agricultural soil. Retrieved gene functions mirrored the main selection pressure as (i) a large fraction of the BPS amplicons contained a high variety of genes/gene clusters related to the degradation of organics including herbicides present in the wastewater and (ii) in the agricultural soil, recovered genes were associated with linuron degradation. Our metagenomic analysis extends observations from cultured isolates and provides evidence that IS1071 is a carrier of catabolic genes in xenobiotica stressed environments and contributes to community level adaptation towards pesticide biodegradation.

Horizontal gene transfer mediated by broad-host-range plasmids is an important mechanism of antibiotic resistance spread. While not all bacteria maintain plasmids equally well, plasmid persistence can improve over time, yet no general evolutionary mechanisms have emerged. Our goal was to identify these mechanisms and to assess if adaptation to one plasmid affects the permissiveness to others. We experimentally evolved Pseudomonas sp. H2 containing multidrug resistance plasmid RP4, determined plasmid persistence and cost using a joint experimental–modelling approach, resequenced evolved clones, and reconstructed key mutations. Plasmid persistence improved in fewer than 600 generations because the fitness cost turned into a benefit. Improved retention of naive plasmids indicated that the host evolved towards increased plasmid permissiveness. Key chromosomal mutations affected two accessory helicases and the RNA polymerase β-subunit. Our and other findings suggest that poor plasmid persistence can be caused by a high cost involving helicase–plasmid interactions that can be rapidly ameliorated. Plasmids facilitate the evolution of antibiotic resistance but little is known about bacteria–plasmid evolution. Here, the authors show that when bacteria adapt to one plasmid, they become generally permissive to plasmid carriage.

Produce is one of the most popular food commodities. Unfortunately, leafy greens can be a reservoir of transferable antibiotic resistance genes. We found that IncF and IncI plasmids were the most prevalent plasmid types in E. coli isolates from produce. This study highlights the importance of the rare microbiome associated with produce as a source of antibiotic resistance genes that might escape cultivation-independent detection, yet may be transferred to human pathogens or commensals., Produce is increasingly recognized as a reservoir of human pathogens and transferable antibiotic resistance genes. This study aimed to explore methods to characterize the transferable resistome of bacteria associated with produce. Mixed salad, arugula, and cilantro purchased from supermarkets in Germany were analyzed by means of cultivation- and DNA-based methods. Before and after a nonselective enrichment step, tetracycline (TET)-resistant Escherichia coli were isolated and plasmids conferring TET resistance were captured by exogenous plasmid isolation. TET-resistant E. coli isolates, transconjugants, and total community DNA (TC-DNA) from the microbial fraction detached from leaves or after enrichment were analyzed for the presence of resistance genes, class 1 integrons, and various plasmids by real-time PCR and PCR-Southern blot hybridization. Real-time PCR primers were developed for IncI and IncF plasmids. TET-resistant E. coli isolated from arugula and cilantro carried IncF, IncI1, IncN, IncHI1, IncU, and IncX1 plasmids. Three isolates from cilantro were positive for IncN plasmids and blaCTX-M-1. From mixed salad and cilantro, IncF, IncI1, and IncP-1β plasmids were captured exogenously. Importantly, whereas direct detection of IncI and IncF plasmids in TC-DNA failed, these plasmids became detectable in DNA extracted from enrichment cultures. This confirms that cultivation-independent DNA-based methods are not always sufficiently sensitive to detect the transferable resistome in the rare microbiome. In summary, this study showed that an impressive diversity of self-transmissible multiple resistance plasmids was detected in bacteria associated with produce that is consumed raw, and exogenous capturing into E. coli suggests that they could transfer to gut bacteria as well.

Self-transmissible and mobilizable plasmids contribute to the emergence and spread of multidrug-resistant bacteria by enabling the horizontal transfer of acquired antibiotic resistance. The objective of this study was to capture and characterize self-transmissible and mobilizable resistance plasmids from a coastal wetland impacted by urban stormwater runoff and human wastewater during the rainy season. Four plasmids were captured, two self-transmissible and two mobilizable, using both mating and enrichment approaches. Plasmid genomes, sequenced with either Illumina or PacBio platforms, revealed representatives of incompatibility groups IncP-6, IncR, IncN3, and IncF. The plasmids ranged in size from 36 to 144 kb and encoded known resistance genes for most of the major classes of antibiotics used to treat Gram-negative infections (tetracyclines, sulfonamides, β-lactams, fluoroquinolones, aminoglycosides, and amphenicols). The mobilizable IncP-6 plasmid pLNU-11 was discovered in a strain of Citrobacter freundii enriched from the wetland sediments with tetracycline and nalidixic acid, and encodes a novel AmpC-like β-lactamase (blaWDC-1), which shares less than 62% amino acid sequence identity with the PDC class of β-lactamases found in Pseudomonas aeruginosa. Although the IncR plasmid pTRE-1611 was captured by mating wetland bacteria with P. putida KT2440 as recipient, it was found to be mobilizable rather than self-transmissible. Two self-transmissible multidrug-resistance plasmids were also captured: the small (48 kb) IncN3 plasmid pTRE-131 was captured by mating wetland bacteria with Escherichia coli HY842 where it is seemed to be maintained at nearly 240 copies per cell, while the large (144 kb) IncF plasmid pTRE-2011, which was isolated from a cefotaxime-resistant environmental strain of E. coli ST744, exists at just a single copy per cell. Furthermore, pTRE-2011 bears the globally epidemic blaCTX-M-55 extended-spectrum β-lactamase downstream of ISEcp1. Our results indicate that urban coastal wetlands are reservoirs of diverse self-transmissible and mobilizable plasmids of relevance to human health.

Good annotation of plasmid genomes is essential to maximise the value of the rapidly increasing volume of plasmid sequences. This short review highlights some of the current issues and suggests some ways forward. Where a well-studied related plasmid system exists we recommend that new annotation adheres to the convention already established for that system, so long as it is based on sound principles and solid experimental evidence, even if some of the new genes are more similar to homologues in different systems. Where a well-established model does not exist we provide generic gene names that reflect likely biochemical activity rather than overall purpose particularly, for example, where genes clearly belong to a type IV secretion system but it is not known whether they function in conjugative transfer or virulence. We also recommend that annotators use a whole system naming approach to avoid ending up with an illogical mixture of names from other systems based on the highest scoring match from a BLAST search. In addition, where function has not been experimentally established we recommend using just the locus tag, rather than a function-related gene name, while recording possible functions as notes rather than in a provisional name.

IS1071, an insertion element that primarily flanks organic xenobiotic degradation genes in cultured isolates, is suggested to play a key role in the formation and distribution of bacterial catabolic pathway gene clusters. However, in environmental settings, the identity of the IS1071 genetic cargo and its correspondence to the local selective conditions remain unknown. To respond, we developed a long-range PCR approach amplifying accessory genes between two IS1071 copies from community DNA followed by amplicon sequencing. We applied this method to pesticide-exposed environments, i.e. linuron-treated agricultural soil and on-farm biopurification systems (BPS) treating complex agricultural wastewater, as to non-treated controls. Amplicons were mainly recovered from the pesticide-exposed environments and the BPS matrix showed a higher size diversity compared to the agricultural soil. Retrieved gene functions mirrored the main selection pressure as (i) a large fraction of the BPS amplicons contained a high variety of genes/gene clusters related to the degradation of organics including herbicides present in the wastewater and (ii) in the agricultural soil, recovered genes were associated with linuron degradation. Our metagenomic analysis extends observations from cultured isolates and provides evidence that IS1071 is a carrier of catabolic genes in xenobiotica stressed environments and contributes to community level adaptation towards pesticide biodegradation.

The emergence and spread of antibiotic resistance is a crisis in health care today. Antibiotic resistance is often horizontally transferred to susceptible bacteria by means of multi-drug resistance plasmids that may or may not persist in the absence of antibiotics. Because bacterial pathogens often grow as biofilms, there is a need to better understand the evolution of plasmid persistence in these environments. Here we compared the evolution of plasmid persistence in the pathogenAcinetobacter baumanniiwhen grown under antibiotic selection in biofilms versus well-mixed liquid cultures. After four weeks, clones in which the plasmid was more stably maintained in the absence of antibiotic selection were present in both populations. On average plasmid persistence increased more in liquid batch cultures, but variation in the degree of persistence was greater among biofilm-derived clones. The results of this study show for the first time that the persistence of MDR plasmids improves in biofilms.

Broad-host-range plasmids play a critical role in the spread of antibiotic resistance and other traits. In spite of increasing information about the genomic diversity of closely related plasmids, the relationship between sequence divergence and host range remains unclear. IncP-1 plasmids are currently classified into six subgroups based on the genetic distance of backbone genes. We investigated whether plasmids from two subgroups exhibit a different host range, using two IncP-1γ plasmids, an IncP-1β plasmid and their minireplicons. Efficiencies of plasmid establishment and maintenance were compared using five species that belong to the Alphaproteobacteria, Betaproteobacteria and Gammaproteobacteria. The IncP-1β plasmid replicated and persisted in all five hosts in the absence of selection. Of the two IncP-1γ plasmids, both were unable to replicate in alphaproteobacterial host Sphingobium japonicum, and one established itself in Agrobacterium tumefaciens but was very unstable. In contrast, both IncP-1γ minireplicons, which produced higher levels of replication initiation protein than the wild-type plasmids, replicated in all strains, suggesting that poor establishment of the native plasmids is in part due to suboptimal replication initiation gene regulation. The findings suggest that host ranges of distinct IncP-1 plasmids only partially overlap, which may limit plasmid recombination and thus result in further genome divergence.

In spite of the contribution of plasmids to the spread of antibiotic resistance in human pathogens, little is known about the transferability of various drug resistance plasmids in bacterial biofilms. The goal of this study was to compare the efficiency of transfer of 19 multidrug resistance plasmids into Escherichia coli recipient biofilms and determine the effects of biofilm age, biofilm-donor exposure time, and donor-to-biofilm attachment on this process. An E. coli recipient biofilm was exposed separately to 19 E. coli donors, each with a different plasmid, and transconjugants were determined by plate counting. With few exceptions, plasmids that transferred well in a liquid environment also showed the highest transferability in biofilms. The difference in transfer frequency between the most and least transferable plasmid was almost a million-fold. The 'invasibility' of the biofilm by plasmids, or the proportion of biofilm cells that acquired plasmids within a few hours, depended not only on the type of plasmid, but also on the time of biofilm exposure to the donor and on the ability of the plasmid donor to attach to the biofilm, yet not on biofilm age. The efficiency of donor strain attachment to the biofilm was not affected by the presence of plasmids. The most invasive plasmid was pHH2-227, which based on genome sequence analysis is a hybrid between IncU-like and IncW plasmids. The wide range in transferability in an E. coli biofilm among plasmids needs to be taken into account in our fight against the spread of drug resistance.

Summary Antibiotic selection drives adaptation of antibiotic resistance plasmids to new bacterial hosts, but the molecular mechanisms are still poorly understood. We previously showed that a broad-host-range plasmid was poorly maintained in Shewanella oneidensis, but rapidly adapted through mutations in the replication initiation gene trfA1. Here we examined if these mutations reduced the fitness cost of TrfA1, and whether this was due to changes in interaction with the host's DNA helicase DnaB. The strains expressing evolved TrfA1 variants showed a higher growth rate than those expressing ancestral TrfA1. The evolved TrfA1 variants showed a lower affinity to the helicase than ancestral TrfA1 and were no longer able to activate the helicase at the oriV without host DnaA. Moreover, persistence of the ancestral plasmid was increased upon overexpression of DnaB. Finally, the evolved TrfA1 variants generated higher plasmid copy numbers than ancestral TrfA1. The findings suggest that ancestral plasmid instability can at least partly be explained by titration of DnaB by TrfA1. Thus under antibiotic selection resistance plasmids can adapt to a novel bacterial host through partial loss of function mutations that simultaneously increase plasmid copy number and decrease unfavorably high affinity to one of the hosts' essential proteins.

Plasmids are extrachromosomal DNA molecules that can potentially be transferred horizontally as well as vertically in bacterial lineages. Plasmids have been shown to code for numerous bacterial phenotypes, the most concerning of which is resistance to commonly used antibiotics. Plasmid-mediated transfer of drug resistance genes has notably led to the rapid spread of resistant strains, which now presents a human health crisis. While the rapid evolution of bacteria due to horizontal plasmid transfer has historically dominated this research field, increasing attention is being paid to the (co)evolutionary interactions that shape the long-term population dynamics of bacteria and plasmids.

Plasmids of the incompatibility group IncP-1 can transfer and replicate in many genera of the Proteobacteria. They are composed of backbone genes that encode a variety of essential functions and accessory genes that have implications for human health and environmental remediation. Although it is well understood that the accessory genes are transferred horizontally between plasmids, recent studies have also provided examples of recombination in the backbone genes of IncP-1 plasmids. As a consequence, phylogeny estimation based on backbone genes is expected to produce conflicting gene tree topologies. The main goal of this study was therefore to infer the evolutionary history of IncP-1 plasmids in the presence of both vertical and horizontal gene transfer. This was achieved by quantifying the incongruence among gene trees and attributing it to known causes such as 1) phylogenetic uncertainty, 2) coalescent stochasticity, and 3) horizontal inheritance. Topologies of gene trees exhibited more incongruence than could be attributed to phylogenetic uncertainty alone. Species-tree estimation using a Bayesian framework that takes coalescent stochasticity into account was well supported, but it differed slightly from the maximum-likelihood tree estimated by concatenation of backbone genes. After removal of the gene that demonstrated a signal of intergroup recombination, the concatenated tree was congruent with the species-tree estimate, which itself was robust to inclusion/exclusion of the recombinant gene. Thus, in spite of horizontal gene exchange both within and among IncP-1 subgroups, the backbone genome of these IncP-1 plasmids retains a detectable vertical evolutionary history.

The CTX-M-type extended-spectrum β-lactamases (ESBLs) present a serious public health threat as they have become nearly ubiquitous among clinical gram-negative pathogens, particularly the enterobacteria. To aid in the understanding and eventual control of the spread of such resistance genes, we sought to determine the diversity of CTX-M ESBLs not among clinical isolates, but in the environment, where weaker and more diverse selective pressures may allow greater enzyme diversification. This was done by examining the CTX-M diversity in municipal wastewater and urban coastal wetlands in southern California, United States, by Sanger sequencing of polymerase chain reaction amplicons. Of the five known CTX-M phylogroups (1, 2, 8, 9, and 25), only genes from groups 1 and 2 were detected in both wastewater treatment plants (WWTPs), and group 1 genes were also detected in one of the two wetlands after a winter rain. The highest relative abundance of blaCTX-M group 1 genes was in the sludge of one WWTP (2.1 × 10(-4) blaCTX-M copies/16S rRNA gene copy). Gene libraries revealed surprisingly high nucleotide sequence diversity, with 157 new variants not found in GenBank, representing 99 novel amino acid sequences. Our results indicate that the resistomes of WWTPs and urban wetlands contain diverse blaCTX-M ESBLs, which may constitute a mobile reservoir of clinically relevant resistance genes.

The World Health Organization has declared the emergence of antibiotic resistance to be a global threat to human health. Broad-host-range plasmids have a key role in causing this health crisis because they transfer multiple resistance genes to a wide range of bacteria. To limit the spread of antibiotic resistance, we need to gain insight into the mechanisms by which the host range of plasmids evolves. Although initially unstable plasmids have been shown to improve their persistence through evolution of the plasmid, the host, or both, the means by which this occurs are poorly understood. Here, we sought to identify the underlying genetic basis of expanded plasmid host-range and increased persistence of an antibiotic resistance plasmid using a combined experimental-modeling approach that included whole-genome resequencing, molecular genetics and a plasmid population dynamics model. In nine of the ten previously evolved clones, changes in host and plasmid each slightly improved plasmid persistence, but their combination resulted in a much larger improvement, which indicated positive epistasis. The only genetic change in the plasmid was the acquisition of a transposable element from a plasmid native to the Pseudomonas host used in these studies. The analysis of genetic deletions showed that the critical genes on this transposon encode a putative toxin–antitoxin (TA) and a cointegrate resolution system. As evolved plasmids were able to persist longer in multiple naive hosts, acquisition of this transposon also expanded the plasmid’s host range, which has important implications for the spread of antibiotic resistance.

Although biofilms represent a common bacterial lifestyle in clinically and environmentally important habitats, there is scant information on the extent of gene transfer in these spatially structured populations. The objective of this study was to gain insight into factors that affect transfer of the promiscuous multidrug resistance plasmid pB10 in Escherichia coli biofilms. Biofilms were grown in different experimental settings, and plasmid transfer was monitored using laser scanning confocal microscopy and plate counting. In closed flow cells, plasmid transfer in surface-attached submerged biofilms was negligible. In contrast, a high plasmid transfer efficiency was observed in a biofilm floating at the air-liquid interface in an open flow cell with low flow rates. A vertical flow cell and a batch culture biofilm reactor were then used to detect plasmid transfer at different depths away from the air-liquid interface. Extensive plasmid transfer occurred only in a narrow zone near that interface. The much lower transfer frequencies in the lower zones coincided with rapidly decreasing oxygen concentrations. However, when an E. coli csrA mutant was used as the recipient, a thick biofilm was obtained at all depths, and plasmid transfer occurred at similar frequencies throughout. These results and data from separate aerobic and anaerobic matings suggest that oxygen can affect IncP-1 plasmid transfer efficiency, not only directly but also indirectly, through influencing population densities and therefore colocalization of donors and recipients. In conclusion, the air-liquid interface can be a hot spot for plasmid-mediated gene transfer due to high densities of juxtaposed donor and recipient cells.

Multidrug resistance bla CMY-2 plasmids that confer resistance to expanded-spectrum cephalosporins have been found in multiple bacterial species collected from different hosts worldwide. The widespread distribution of bla CMY-2 plasmids may be driven by antibiotic use that selects for the dissemination and persistence of these plasmids. Alternatively, these plasmids may persist and spread in bacterial populations in the absence of selection pressure if a balance exists among conjugative transfer, segregation loss during cell division, and fitness cost to the host. We conducted a series of experiments (both in vivo and in vitro ) to study these mechanisms for three bla CMY-2 plasmids, peH4H, pAR060302, and pAM04528. Results of filter mating experiments showed that the conjugation efficiency of bla CMY-2 plasmids is variable, from −7 for pAM04528 and peH4H to ∼10 −3 for pAR060302. Neither peH4H nor pAM04528 was transferred from Escherichia coli strain DH10B, but peH4H was apparently mobilized by the coresident trimethoprim resistance-encoding plasmid pTmpR. Competition studies showed that carriage of bla CMY-2 plasmids imposed a measurable fitness cost on the host bacteria both in vitro (0.095 to 0.25) and in vivo (dairy calf model). Long-term passage experiments in the absence of antibiotics demonstrated that plasmids with limited antibiotic resistance phenotypes arose, but eventually drug-sensitive, plasmid-free clones dominated the populations. Given that plasmid decay or loss is inevitable, we infer that some level of selection is required for the long-term persistence of bla CMY-2 plasmids in bacterial populations.

Bioactive amounts of antibiotics as well as resistant bacteria reach the soil through manure fertilization. We investigated plasmids that may stimulate the environmental spread and interspecies transfer of antibiotic resistance. After treatment of two soils with manure, either with or without the sulfonamide antibiotic sulfadiazine, a significant increase in copies of the sulfonamide resistance gene sul2 was detected by qPCR. All sul2 carrying plasmids, captured in Escherichia coli from soil, belonged to a novel class of self-transferable replicons. Manuring and sulfadiazine significantly increased the abundance of this replicon type in a chemically fertilized but not in an annually manured soil, as determined by qPCR targeting a transfer gene. Restriction patterns and antibiograms showed a considerable diversity within this novel plasmid group. Analysis of three complete plasmid sequences revealed a conserved 30 kbp backbone with only 36% G+C content, comprised of transfer and maintenance genes with moderate homology to plasmid pIPO2 and a replication module (rep and oriV) of other descent. The plasmids differed in composition of the 27.0-28.3 kbp accessory region, each of which carried ISCR2 and several resistance genes. Acinetobacter spp. was identified as a potential host of such LowGC-type plasmids in manure and soil.

In spite of the importance of plasmids in bacterial adaptation, we have a poor understanding of their dynamics. It is not known if or how plasmids persist in and spread through (invade) a bacterial population when there is no selection for plasmid-encoded traits. Moreover, the differences in dynamics between spatially structured and mixed populations are poorly understood. Through a joint experimental/theoretical approach, we tested the hypothesis that self-transmissible IncP-1 plasmids can invade a bacterial population in the absence of selection when initially very rare, but only in spatially structured habitats and when nutrients are regularly replenished. Using protocols that differed in the degree of spatial structure and nutrient levels, the invasiveness of plasmid pB10 in Escherichia coli was monitored during at least 15 days, with an initial fraction of plasmid-bearing (p(+)) cells as low as 10(-7). To further explore the mechanisms underlying plasmid dynamics, we developed a spatially explicit mathematical model. When cells were grown on filters and transferred to fresh medium daily, the p(+) fraction increased to 13%, whereas almost complete invasion occurred when the population structure was disturbed daily. The plasmid was unable to invade in liquid. When carbon source levels were lower or not replenished, plasmid invasion was hampered. Simulations of the mathematical model closely matched the experimental results and produced estimates of the effects of alternative experimental parameters. This allowed us to isolate the likely mechanisms most responsible for the observations. In conclusion, spatial structure and nutrient availability can be key determinants in the invasiveness of plasmids.

Conjugative plasmids play a very important role in bacterial adaptation through the dissemination of useful traits. Incompatibility group P-1 (IncP-1) plasmids exhibit an extreme broad-host-range among Gram-negative bacteria and known to be one of the major agents to disseminate various phenotypic traits such as antibiotic resistance and xenobiotic degradation. Although the plasmids are believed to be very stable in most Gram-negative bacteria, little is known about the factors that affect their stability in various hosts, allowing their persistence in bacterial population. Here we show that the stability of the cryptic IncP-1β plasmid pBP136 differed greatly in four different Escherichia coli K12 host backgrounds (MG1655, DH5α, EC100, and JM109), whereas the closely related plasmid pB10 was stable in all four strains. The supply of the kleF gene, which is involved in the stability of IncP-1 plasmids but absent in pBP136, did not improve the stability of the plasmid. Our findings suggest that persistence of IncP-1 plasmids in the absence of selection is affected by strain-specific factors.

The overall architecture of IncP-1 plasmids is very conserved in that the accessory genes are typically located in one or two specific regions: between oriV and trfA and between the tra and trb operons. Various hypotheses have been formulated to explain this, but none have been tested experimentally. We investigated whether this structural similarity is due to region-specific transposition alone or also is reliant on selection for plasmids with insertions limited to these two regions. We first examined the transposition of Tn 21 Km into IncP-1β plasmid pBP136 and found that most Tn 21 Km insertions (67%) were located around oriV . A similar experiment using the oriV region of IncP-1β plasmid pUO1 confirmed these results. We then tested the transferability, stability, and fitness cost of different pBP136 derivatives to determine if impairment of these key plasmid characters explained the conserved plasmid architecture. Most of the pBP136 derivatives with insertions in transfer genes were no longer transferable. The plasmids with insertions in the oriV - trfA and tra-trb regions were more stable than other plasmid variants, and one of these also showed a significantly lower fitness cost. In addition, our detailed sequence analysis of IncP-1 plasmids showed that Tn 402 / 5053 -like transposons are situated predominantly between the tra and trb operons and close to the putative resolution site for the ParA resolvase, a potential hot spot for those transposons. Our study presents the first empirical evidence that region-specific insertion of transposons in combination with selection for transferable and stable plasmids explains the structural similarity of IncP-1 plasmids.

IncP-1 plasmids are known to be promiscuous, but it is not understood if they are equally well adapted to various species within their host range. Moreover, little is known about their fate in bacterial communities. We determined if the IncP-1β plasmid pB10 was unstable in some Proteobacteria, and whether plasmid stability was enhanced after long-term carriage in a single host and when regularly switched between isogenic hosts. Plasmid pB10 was found to be very unstable in Pseudomonas putida H2, and conferred a high cost (c. 20% decrease in fitness relative to the plasmid-free host). H2(pB10) was then evolved under conditions that selected for plasmid maintenance, with or without regular plasmid transfer (host-switching). When tested in the ancestral host, the evolved plasmids were more stable and their cost was significantly reduced (9% and 16% for plasmids from host-switched and nonswitched lineages, respectively). Our findings suggest that IncP-1 plasmids can rapidly adapt to an unfavorable host by improving their overall stability, and that regular conjugative transfer accelerates this process.

Broad-host-range (BHR) IncP-1 plasmids have the ability to transfer between and replicate in nearly all species of the Alpha-, Beta- and Gammaproteobacteria, but surprisingly few data are available on the stability of these plasmids in strains within their host range. Moreover, even though molecular interactions between the bacterial host and its plasmid(s) exist, no systematic study to date has compared the stability of the same plasmid among different hosts. The goal of this study was to examine whether the stability characteristics of an IncP-1 plasmid can be variable between strains within the host range of the plasmid. Therefore, 19 strains within the Alpha-, Beta- or Gammaproteobacteria carrying the IncP-1beta plasmid pB10 were serially propagated in non-selective medium and the fraction of segregants was monitored through replica-picking. Remarkably, a large variation in the stability of pB10 in different strains was found, even between strains within the same genus or species. Ten strains showed no detectable plasmid loss over about 200 generations, and in two strains plasmid-free clones were only sporadically observed. In contrast, three strains, Pseudomonas koreensis R28, Pseudomonas putida H2 and Stenotrophomonas maltophilia P21, exhibited rapid plasmid loss within 80 generations. Parameter estimation after mathematical modelling of these stability data suggested high frequencies of segregation (about 0.04 per generation) or high plasmid cost (i.e. a relative fitness decrease in plasmid-bearing cells of about 15 and 40 %), which was confirmed experimentally. The models also suggested that plasmid reuptake by conjugation only played a significant role in plasmid stability in one of the three strains. Four of the 19 strains lost the plasmid very slowly over about 600 generations. The erratic decrease of the plasmid-containing fraction and simulation of the data with a new mathematical model suggested that plasmid cost was variable over time due to compensatory mutations. The findings of this study demonstrate that the ability of a so-called 'BHR' plasmid to persist in a bacterial population is influenced by strain-specific traits, and therefore observations made for one strain should not be generalized for the entire species or genus.

It has been predicted that precipitation and atmospheric nitrogen (N) deposition will increase in northern China; yet, ecosystem responses to the interactive effects of water and N remain largely unknown. In particular, responses of belowground microbial community to projected global change and their potential linkages to aboveground macro-organisms are rarely studied. In this study, we examined the responses of soil bacterial diversity and community composition to increased precipitation and multi-level N deposition in a temperate steppe in Inner Mongolia, China, and explored the diversity linkages between aboveground and belowground communities. It was observed that N addition caused the significant decrease in bacterial alpha-diversity and dramatic changes in community composition. In addition, we documented strong correlations of alpha- and beta-diversity between plant and bacterial communities in response to N addition. It was found that N enriched the so-called copiotrophic bacteria, but reduced the oligotrophic groups, primarily by increasing the soil inorganic N content and carbon availability and decreasing soil pH. We still highlighted that increased precipitation tended to alleviate the effects of N on bacterial diversity and dampen the plant-microbe connections induced by N. The counteractive effects of N addition and increased precipitation imply that even though the ecosystem diversity and function are predicted to be negatively affected by N deposition in the coming decades; the combination with increased precipitation may partially offset this detrimental effect.

Plasmids are important vehicles for rapid adaptation of bacterial populations to changing environmental conditions. It is thought that to reduce the cost of plasmid carriage, only a fraction of a local population carries plasmids or is permissive to plasmid uptake. Plasmids provide various accessory traits which might be beneficial under particular conditions. The genetic variation generated by plasmid carriage within populations ensures the robustness toward environmental changes. Plasmid-mediated gene transfer plays an important role not only in the mobilization and dissemination of antibiotic resistance genes but also in the spread of degradative pathways and pathogenicity determinants of pathogens. Here we summarize the state-of-the-art methods to study the occurrence, abundance, and diversity of plasmids in environmental bacteria. Increasingly, cultivation-independent total-community DNA-based methods are being used to characterize and quantify the diversity and abundance of plasmids in relation to various biotic and abiotic factors. An improved understanding of the ecology of plasmids and their hosts is crucial in the development of intervention strategies for antibiotic-resistance-gene spread. We discuss the potentials and limitations of methods used to determine the host range of plasmids, as the ecology of plasmids is tightly linked to their hosts. The recent advances in sequencing technologies provide an enormous potential for plasmid classification, diversity, and evolution studies, but numerous challenges still exist.

We report the draft genome sequence of Pseudomonas sp. nov. H2, isolated from creek sediment in Moscow, ID, USA. The strain is most closely related to Pseudomonas putida . However, it has a slightly smaller genome that appears to have been impacted by horizontal gene transfer and poorly maintains IncP-1 plasmids.

A self-transmissible broad-host-range (BHR) plasmid pSFA231 was isolated from petroleum-contaminated sediment in Shen-fu wastewater irrigation zone, China, using the triparental mating exogenous plasmid capture method. Based on its complete sequence the plasmid has a size of 41.5 kb and codes for 50 putative open reading frames (orfs), 28 of which represent genes involved in replication, partitioning and transfer functions of the plasmid. Phylogenetic analysis grouped pSFA231 into the newly defined PromA plasmid family, which currently includes five members. Further comparative genomic analysis shows that pSFA231 shares the common backbone regions with the other PromA plasmids, i.e., genes involved in replication, maintenance and control, and conjugative transfer. Nevertheless, phylogenetic divergence was found in specific gene products. We propose to divide the PromA group into two subgroups, PromA-α (pMRAD02, pSB102) and PromA-β (pMOL98, pIPO2T, pSFA231, pTer331), based on the splits network analysis of the RepA protein. Interestingly, a cluster of hypothetical orfs located between parA and traA of pSFA231 shows high similarity with the corresponding regions on pMOL98, pIPO2T and pTer331, suggesting these hypothetical orfs may represent ‘‘essential’’ plasmid backbone genes for the PromA-β subgroup. Alternatively, they may also be accessory genes that were first acquired and then stayed as the plasmid diverged. Our study increases the available collection of complete genome sequences of BHR plasmids, and since pSFA231 is the only characterized PromA plasmid from China, our findings also enhance our understanding of the genetic diversity of this plasmid group in different parts of the world.