Search

The ability of bacterial plasmids to adapt to novel hosts and thereby shift their host range is key to their long-term persistence in bacterial communities. Promiscuous plasmids of the incompatibility group P (IncP)-1 can colonize a wide range of hosts, but it is not known if and how they can contract, shift or further expand their host range. To understand the evolutionary mechanisms of host range shifts of IncP-1 plasmids, an IncP-1β mini-replicon was experimentally evolved in four hosts in which it was initially unstable. After 1000 generations in serial batch cultures under antibiotic selection for plasmid maintenance (kanamycin resistance), the stability of the mini-plasmid dramatically improved in all coevolved hosts. However, only plasmids evolved in Shewanella oneidensis showed improved stability in the ancestor, indicating that adaptive mutations had occurred in the plasmid itself. Complete genome sequence analysis of nine independently evolved plasmids showed seven unique plasmid genotypes that had various kinds of single mutations at one locus, namely, the N-terminal region of the replication initiation protein TrfA. Such parallel evolution indicates that this region was under strong selection. In five of the seven evolved plasmids, these trfA mutations resulted in a significantly higher plasmid copy number. Evolved plasmids were found to be stable in four other naive hosts, but could no longer replicate in Pseudomonas aeruginosa. This study shows that plasmids can specialize to a novel host through trade-offs between improved stability in the new host and the ability to replicate in a previously permissive host.

Plasmids are ubiquitous mobile elements that serve as a pool of many host beneficial traits such as antibiotic resistance in bacterial communities. To understand the importance of plasmids in horizontal gene transfer, we need to gain insight into the 'evolutionary history' of these plasmids, i.e. the range of hosts in which they have evolved. Since extensive data support the proposal that foreign DNA acquires the host's nucleotide composition during long-term residence, comparison of nucleotide composition of plasmids and chromosomes could shed light on a plasmid's evolutionary history. The average absolute dinucleotide relative abundance difference, termed delta-distance, has been commonly used to measure differences in dinucleotide composition, or 'genomic signature', between bacterial chromosomes and plasmids. Here, we introduce the Mahalanobis distance, which takes into account the variance-covariance structure of the chromosome signatures. We demonstrate that the Mahalanobis distance is better than the delta-distance at measuring genomic signature differences between plasmids and chromosomes of potential hosts. We illustrate the usefulness of this metric for proposing candidate long-term hosts for plasmids, focusing on the virulence plasmids pXO1 from Bacillus anthracis, and pO157 from Escherichia coli O157:H7, as well as the broad host range multi-drug resistance plasmid pB10 from an unknown host.

Little is known about the range of hosts in which broad-host-range (BHR) plasmids can persist in the absence of selection for plasmid-encoded traits, and whether this "long-term host range" can evolve over time. Previously, the BHR multidrug resistance plasmid pB10 was shown to be highly unstable in Stenotrophomonas maltophilia P21 and Pseudomonas putida H2. To investigate whether this plasmid can adapt to such unfavorable hosts, we performed evolution experiments wherein pB10 was maintained in strain P21, strain H2, and alternatingly in P21 and H2. Plasmids that evolved in P21 and in both hosts showed increased stability and decreased cost in ancestral host P21. However, the latter group showed higher variability in stability patterns, suggesting that regular switching between distinct hosts hampered adaptive plasmid evolution. The plasmids evolved in P21 were also equally or more stable in other hosts compared to pB10, which suggested true host-range expansion. The complete genome sequences of four evolved plasmids with improved stability showed only one or two genetic changes. The stability of plasmids evolved in H2 improved only in their coevolved hosts, not in the ancestral host. Thus a BHR plasmid can adapt to an unfavorable host and thereby expand its long-term host range.

Recent studies have indicated that the evolutionarily common catabolic gene clusters are loaded on structurally diverse toluene-catabolic (TOL) plasmids and their residing transposons. To elucidate the mechanisms supporting the diversification of catabolic plasmids and transposons, we determined here the complete 107,929 bp sequence of pWW53, a TOL plasmid from Pseudomonas putida MT53. pWW53 was found to belong to the IncP-7 incompatibility group that play important roles in the catabolism of several xenobiotics. pWW53 carried two distinct transposase-resolvase gene clusters (tnpAR modules), five short terminal inverted repeats (IRs), and three site-specific resolution (res) sites that are all typical of class II transposons. This organization of pWW53 suggested the four possible transposable regions, Tn4657 to Tn4660. The largest 86 kb region (Tn4657) spanned the three other regions, and Tn4657 and Tn4660 (62 kb) covered all of the 36 xyl genes for toluene catabolism. Our subsequent transposition experiments clarified that the three transposons, Tn4657 to Tn4659, indeed exhibit their transposability, and that pWW53 also generated another 37 kb toluene-catabolic transposon, Tn4656, which carried the two separated and inversely oriented segments of pWW53: the tnpRA-IR module of Tn4658 and a part of xyl gene clusters on Tn4657. The Tn4658 transposase was able to mediate the transposition of Tn4658, Tn4657, and Tn4656, while the Tn4659 transposase catalyzed only the transposition of Tn4659. Tn4656 was formed by the Tn4658 resolvase-mediated site-specific inversion between the two inversely oriented res sites on pWW53. These findings and comparison with other catabolic plasmids clearly indicate multiple copies of transposition-related genes and sites on one plasmid and their recombination activities contribute greatly to the diversification of plasmid structures as well as wide dissemination of the evolutionary common gene clusters in various plasmids.

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 Tn21Km into IncP-1beta plasmid pBP136 and found that most Tn21Km insertions (67%) were located around oriV. A similar experiment using the oriV region of IncP-1beta 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 Tn402/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.

Two different cultivation-independent approaches were applied to isolate genes for naphthalene dioxygenase (NDO) from oil-contaminated soil in Japan. One approach was the construction of a broad-host-range cosmid-based metagenomic DNA library, and the other was the so-called exogenous plasmid isolation technique. Our screening of NDO genes in both approaches was based on the functional complementation of Pseudomonas putida strains which contained Tn4655K, a transposon carrying the entire set of naphthalene-catabolic (nah) genes but lacking the NDO-encoding gene. We obtained in the former approach a cosmid clone (pSLX928-6) that carried an nah upper pathway operon for conversion of naphthalene to salicylate, and this operon showed a significantly high level of similarity to the corresponding operon on an IncP-9 naphthalene-catabolic plasmid, pDTG1. In the latter approach, the microbial fraction from the soil was mated with a plasmid-free P. putida strain containing a chromosomal copy of Tn4655K, and transconjugants were obtained that received either a 200- or 80-kb plasmid containing all the nah genes for the complete degradation of naphthalene. Subsequent analysis revealed that (1) both plasmids belong to the IncP-9 incompatibility group; (2) their nah upper pathway operons are significantly similar, but not completely identical, to those of pDTG1 and pSLX928-6; and (3) these plasmids carried genes for the salicylate metabolism by the meta-cleavage pathway.

The complete 41,268 bp nucleotide sequence of the IncP-1beta plasmid pBP136 from the human pathogen Bordetella pertussis, the primary aetiological agent of whooping cough, was determined and analysed. This plasmid carried a total of 46 ORFs: 44 ORFs corresponding to the genes in the conserved IncP-1beta backbone, and 2 ORFs similar to the XF1596 and XF1597 genes with unknown function of the plant pathogen Xylella fastidiosa. Interestingly, pBP136 had no accessory genes carrying genetic traits such as antibiotic or mercury resistance and/or xenobiotic degradation. Moreover, pBP136 had only two of the kle genes (kleAE) that have been reported to be important for the stability of IncP-1 plasmid in Pseudomonas aeruginosa. Phylogenetic analysis of the Kle proteins revealed that the KleA and KleE of pBP136 were phylogenetically distant from those of the present IncP-1 plasmids. In contrast, IncC1 and KorC, encoded upstream and downstream of the kle genes respectively, and the replication-initiation protein, TrfA, were closely related to those of the IncP-1beta 'R751 group'. These results suggest that (i) pBP136 without any apparent accessory genes diverged early from an ancestor of the present IncP-1beta plasmids, especially those of the R751 group, and (ii) the kle genes might be incorporated independently into the backbone region of the IncP-1 plasmids for their stable maintenance in various host cells.

The naphthalene-catabolic (nah) genes on the incompatibility group P-9 (IncP-9) self-transmissible plasmid NAH7 from Pseudomonas putida G7 are some of the most extensively characterized genetic determinants for bacterial aerobic catabolism of aromatic hydrocarbons. In contrast to the detailed studies of its catabolic cascade and enzymatic functions, the biological characteristics of plasmid NAH7 have remained unclear. Our sequence determination in this study together with the previously deposited sequences revealed the entire structure of NAH7 (82,232 bp). Comparison of NAH7 with two other completely sequenced IncP-9 catabolic plasmids, pDTG1 and pWW0, revealed that the three plasmids share very high nucleotide similarities in a 39-kb region encoding the basic plasmid functions (the IncP-9 backbone). The backbone of NAH7 is phylogenetically more related to that of pDTG1 than that of pWW0. These three plasmids carry their catabolic gene clusters at different positions on the IncP-9 backbone. All of the NAH7-specified nah genes are located on a class II transposon, Tn4655. Our analysis of the Tn4655-encoded site-specific recombination system revealed that (i) a novel tyrosine recombinase, TnpI, catalyzed both the intra- and intermolecular recombination between two copies of the attI site, (ii) the functional attI site was located within a 119-bp segment, and (iii) the site-specific strand exchange occurred within a 30-bp segment in the 41-bp CORE site. Our results and the sequence data of other naphthalene-catabolic plasmids, pDTG1 and pND6-1, suggest a potential role of the TnpI-attI recombination system in the establishment of these catabolic plasmids.

The self-transmissible plasmid pUO1 from Delftia acidovorans strain B carries two haloacetate-catabolic transposons, TnHad1 and TnHad2, and the mer genes for resistance to mercury. The complete 67,066-bp sequence of pUO1 revealed that the mer genes were also carried by two Tn402/Tn5053-like transposons, Tn4671 and Tn4672, and that the pUO1 backbone regions shared 99% identity to those of the archetype IncP-1beta plasmid R751. Comparison of pUO1 with three other IncP-1beta plasmids illustrated the importance of transposon insertion in the diversity and evolution of this group of plasmids. Mutational analysis of the four outermost residues in the inverted repeats (IRs) of TnHad2, a Tn21-related transposon, revealed a crucial role of the second residue of its IRs in transposition.

Various xenobiotic-degrading genes on many catabolic plasmids are often flanked by two copies of an insertion sequence, IS1071. This 3.2-kb IS element has long (110-bp) terminal inverted repeats (IRs) and a transposase gene that are phylogenetically related to those of the class II transposons. However, the transposition mechanism of IS1071 has remained unclear. Our study revealed that IS1071 was only able to transpose at high frequencies in two environmental β-proteobacterial strains, Comamonas testosteroni and Delftia acidovorans, and not in any of the bacteria examined which belong to the α- and γ-proteobacteria. IS1071 was found to have the functional features of the class II transposons in that (i) the final product of the IS1071 transposition was a cointegrate of its donor and target DNA molecules connected by two directly repeated copies of IS1071, one at each junction; (ii) a 5-bp duplication of the target sequence was observed at the insertion site; and (iii) a tnpA mutation of IS1071 was efficiently complemented by supplying the wild-type tnpA gene in trans. Deletion analysis of the IS1071 IR sequences indicated that nearly the entire region of the IRs was required for its transposition, suggesting that the interaction between the transposase and IRs of IS1071 might be different from that of the other well-characterized class II transposons. [ABSTRACT FROM AUTHOR]