- 1Institute for Multidisciplinary Research in Applied Biology, Universidad Pública de Navarra (UPNA), Edificio de Agrobiotecnología, Mutilva Baja, Spain
- 2Área de Genética, Facultad de Ciencias, Universidad de Málaga, Málaga, Spain
- 3Instituto de Hortofruticultura Subtropical y Mediterránea “La Mayora”, Consejo Superior de Investigaciones Científicas (IHSM-UMA-CSIC), Málaga, Spain
Pseudomonas syringae pv. savastanoi NCPPB 3335 is the causal agent of olive knot disease and contains three virulence plasmids: pPsv48A (pA), 80 kb; pPsv48B (pB), 45 kb, and pPsv48C (pC), 42 kb. Here we show that pB contains a complete MPFT (previously type IVA secretion system) and a functional origin of conjugational transfer adjacent to a relaxase of the MOBP family; pC also contains a functional oriT-MOBP array, whereas pA contains an incomplete MPFI (previously type IVB secretion system), but not a recognizable oriT. Plasmid transfer occurred on solid and in liquid media, and on leaf surfaces of a non-host plant (Phaseolus vulgaris) with high (pB) or moderate frequency (pC); pA was transferred only occasionally after cointegration with pB. We found three plasmid-borne and three chromosomal relaxase genes, although the chromosomal relaxases did not contribute to plasmid dissemination. The MOBP relaxase genes of pB and pC were functionally interchangeable, although with differing efficiencies. We also identified a functional MOBQ mobilization region in pC, which could only mobilize this plasmid. Plasmid pB could be efficiently transferred to strains of six phylogroups of P. syringae sensu lato, whereas pC could only be mobilized to two strains of phylogroup 3 (genomospecies 2). In two of the recipient strains, pB was stably maintained after 21 subcultures in liquid medium. The carriage of several relaxases by the native plasmids of P. syringae impacts their transfer frequency and, by providing functional diversity and redundancy, adds robustness to the conjugation system.
Introduction
Approximately three quarters of the protein families in proteobacterial genomes have been acquired by horizontal gene transfer (HGT) during evolution (Popa and Dagan, 2011). In more recent evolutionary times, HGT also has a substantial impact on human life, having contributed to the emergence and dissemination of clinical bacteria with multiple resistance to antibiotics, as well as devastating plant pathogens (O'Brien et al., 2012; Martínez, 2018). Comparable accelerated evolution by HGT has also allowed the emergence of new bacterial plant diseases and the distribution of resistance to phytosanitary compounds in agricultural settings (Vivian et al., 2001; Baltrus et al., 2011; Firrao et al., 2018).
Plasmids are the main mobile genetic elements (MGE) driving HGT (Halary et al., 2010), and they can be transferred among distantly-related organisms by conjugation, a process mediating the transfer of DNA and proteins from a donor to a recipient through cell-to-cell contact (de la Cruz et al., 2010; Guzmán-Herrador and Llosa, 2019). Generally, conjugative transfer requires the expression of genes for three separate and specialized systems: (1) a system involved in mating-pair formation (MPF), mediating the physical transfer of DNA from the donor cell to the recipient cell, (2) a DNA transfer and replication (Dtr) delivery system, allowing the specific identification of mobilizable DNA and its delivery to the recipient cell through the MPF system, and (3) a type IV coupling protein, which is an ATPase that couples the two previous systems and is thought to energize the DNA transfer process.
The MPF system involves specialized type IV secretion systems (T4SS), which secrete proteins and nucleoprotein complexes across membranes (Backert and Grohmann, 2017; Rapisarda and Fronzes, 2018). There are eight main classes of MPF systems, class MPFT (previously named type IVA secretion system, or T4ASS) being the most abundant (Guglielmini et al., 2014). The members of these classes are functionally and structurally similar and their constituent proteins can often be exchanged among them. These systems consist of a highly variable number of genes, coding for proteins that assemble into a barrel-like complex spanning the cell envelope and finishing in a conjugation pilus. This pilus makes contact with the recipient cell and serves as a conduit to deliver the DNA as a nucleoprotein adduct.
The main components of the Dtr system are the origin of transfer (oriT) and the cognate relaxase, which, together with diverse auxiliary proteins, form the relaxosome. The oriT spans no more than one hundred to a few hundred base pairs and contains a few conserved motifs (Zrimec and Lapanje, 2018), specifically identifying the DNA that can be transferred by conjugation. The relaxase is an endonuclease that cleaves the oriT in a site- and strand-specific manner forming a covalent bond with its 5′ end, delivering it to the MPF system and piloting the transfer of this single-stranded DNA to the recipient cell. The DNA mobilization region usually contains the oriT located in close association with the relaxase gene and a gene for a relaxase accessory protein (RAP). Unlike oriT and the relaxase, RAPs are usually not essential for conjugation but increase the mobilization frequency (Backert and Grohmann, 2017; Garcillán-Barcia et al., 2019; Guzmán-Herrador and Llosa, 2019). Relaxases are organized into nine classes, with classes MOBP and MOBQ being prominent in Gammaproteobacteria (Garcillán-Barcia et al., 2020). In general, relaxases are highly specific for their cognate oriT sequences, although they can recognize heterologous sequences with a lower efficiency and, with few exceptions, can act both in cis and in trans (Garcillán-Barcia et al., 2019; Guzmán-Herrador and Llosa, 2019). Therefore, mobilizable plasmids are those containing an oriT, and generally the cognate RAP and relaxase genes, whereas conjugative plasmids contain both the MPF and the Dtr systems (Smillie et al., 2010).
The plant pathogenic gammaproteobacterium Pseudomonas syringae sensu lato (hereafter P. syringae) is a taxonomic complex of utmost economic and research importance, causing relevant diseases in nearly all cultivated herbaceous and woody plants (Mansfield et al., 2012). Individual isolates, however, usually have a narrow and distinctive plant host-range that determines their separation into more than 60 pathovars. The rather large pangenome of P. syringae (Gomila et al., 2017; Dillon et al., 2019) suggests a very high rate of gene acquisition by horizontal transfer (Brockhurst et al., 2019), which could be mediated by plasmids and/or other MGEs (Sundin and Murillo, 2009; Jackson et al., 2011). Most strains of P. syringae usually contain one or more native plasmids that increase bacterial fitness by carrying genes for virulence, conferring resistance to chemicals used in disease control (such as streptomycin or copper), or resistance to UV light, among many others traits. Field and laboratory experiments demonstrated the conjugative transfer of plasmids carrying copper and streptomycin resistance genes among P. syringae strains. Additionally, comparative phylogenetic analyses demonstrate an exchange of genes and plasmid fragments within populations and species of P. syringae (Vivian et al., 2001; Ma et al., 2007; Sundin and Murillo, 2009). In fact, complete MPFT and MPFI (previously named type IVB secretion system, or T4BSS) systems, as well as genes from both classes, have been found in native plasmids from P. syringae (Stavrinides and Guttman, 2004; Zhao et al., 2005; Pérez-Martínez et al., 2008; Gutiérrez-Barranquero et al., 2017). However, neither the functionality of the MPF and Dtr systems and their role in conjugation and mobilization of native plasmids from P. syringae, nor the host range of P. syringae plasmids have been explored.
P. syringae pv. savastanoi NCPPB 3335 causes aerial tumors in olive (Olea europaea) plants, and contains three virulence plasmids: pPsv48A (pA, 80 kb), pPsv48B (pB, 45 kb), and pPsv48C (pC, 42 kb). The three plasmids are individually required to induce full symptoms and to reach high populations in planta (Bardaji et al., 2011b; Castañeda-Ojeda et al., 2017a,b; Añorga et al., 2020). Sequencing of the plasmid complement of NCPPB 3335 identified 15 coding sequences (CDSs) that might constitute a complete MPFT in pB, plus a well-conserved putative origin of transfer (oriT) in pB and in pC, and an incomplete MPFI system, but not a recognizable oriT, in pA (Bardaji et al., 2011b).
Here we demonstrated that the MPFT system of pB participates in the conjugal transfer of pB and pC in solid and liquid media and on plant leaves. The MOBP relaxases of pB and pC were interchangeable, although with differing efficiencies, and we also identified a functional MOBQ mobilization region in pC. pB could be mobilized to diverse phylogroups of P. syringae, where they were stably maintained after 21 subcultures.
Materials and methods
Bacterial strains, plasmids and growth conditions
Relevant bacterial strains, plasmids and constructions used in this study are listed in Supplementary Table S1. Pseudomonas (at 28°C) and E. coli (at 37°C) strains were propagated using LB medium (Sambrook et al., 1989), due to the generation of small and round, well-defined colonies on this medium favoring colony counting. Liquid cultures were always incubated with orbital shaking (200 rpm). Counterselection of cells carrying the sacB gene, which confers lethality in the presence of sucrose, was done in nutrient agar (NA) medium (Oxoid, Basingstoke, UK) supplemented with 5% sucrose (medium SNA). Conjugations were carried out onto plates of LB, medium B (KMB) (King et al., 1954) or Minimal A Medium (MAM) (Miller, 1992); MAM was routinely supplemented with a final concentration of 0.04 g L−1 of casamino acids and 2 g L−1 of glucose (MacDonald et al., 1986). When necessary, and unless otherwise indicated, media were supplemented with (final concentrations, in μg mL−1): ampicillin (Amp), 100; gentamicin (Gm), 12.5; kanamycin (Km), 7 for P. syringae and 25 for E. coli; tetracycline (Tc), 3 for P. syringae and 12.5 for E. coli; spectinomycin (Sp), 25; streptomycin (Sm), 50.
General molecular procedures and bioinformatics
DNA was amplified using a high fidelity polimerase (PrimeStar HS, Takara Bio Inc., Japan or Q5 High-Fidelity DNA polymerase, New England Biolabs), or a standard enzyme (BIOTaq, Bioline, UK), using primers detailed in Supplementary Table S2, and cloned using the CloneJET PCR Cloning Kit (Thermo Scientific) or the pGEM-T Easy Vector System (Promega). Purification of plasmids from E. coli was done following a boiling method (Holmes and Quigley, 1981) or using a commercial kit (Illustra PlasmidPrep Mini Spin Kit, GE Healthcare). For plasmid profile gels, DNA was purified by alkaline lysis and separated by electrophoresis in 0.8% agarose gels with 1xTAE as described (Sesma et al., 2000). Plasmids were transferred to P. syringae by electroporation (Choi et al., 2006).
The NCBI BLAST programs (Johnson et al., 2008) and EMBL-EBI (https://www.ebi.ac.uk/services/) server tools (Hubbard et al., 2008) were used to compare and search for homology among sequences and to produce sequence alignments. Alignments were shaded using the Sequence Manipulation suite (Stothard, 2000). BProm (Solovyev and Salamov, 2011) and Mfold (Zuker, 2003) were used to predict binding sites and inverted repeats in putative oriT sequences. Search for protein motifs was done using the InterPro interface (Mitchell et al., 2015) (http://www.ebi.ac.uk/interpro/). Genome and nucleotide sequence visualization and manipulation was done using the Artemis genome browser and ACT (Carver et al., 2008). Primers were designed using the Primer3plus software (Untergasser et al., 2012).
Search for T4SSs was done using the web servers for SecReT4 2.0 (https://bioinfo-mml.sjtu.edu.cn/SecReT4/) (Bi et al., 2013) and oriTfinder (https://bioinfo-mml.sjtu.edu.cn/oriTfinder/) (Li et al., 2018), using default parameters. Hits were considered significant with an E-value lower than 10−5. We uploaded nucleotide fasta files for all searches because we noticed that, at least for SecReT4, detection was more sensitive in this way than using annotated nucleotide files. To identify relaxes MOB families, the complete set of annotated CDSs of plasmids pA, pB and pC were analyzed with the web server for MOBscan (https://castillo.dicom.unican.es/mobscan/) (Garcillán-Barcia et al., 2020). Search for protein motifs and fold recognition was done using the Pfam (http://pfam.xfam.org/; database version 35.0) (Mistry et al., 2020) and the Phyre2 v 2.0 (Kelley et al., 2015) web servers.
Construction of mutants
For the construction of mutants by marker-exchange, DNA regions of 0.7–1 kb flanking the target sequence for mutagenesis were selected to insert an antibiotic resistance gene between them and either introduce a deletion or interrupt the gene of interest. The flanking DNA fragments were amplified by PCR using a high-fidelity polimerase and primers containing appropriate restriction sites on both ends (Supplementary Table S2). Amplicons were separately cloned using the CloneJET PCR Cloning Kit (Thermo Scientific) and then cloned in their native orientation into pK18mobsacB (KmR) or pJQ200SK (GmR) as a single fragment, with an appropriate antibiotic resistance gene separating the two fragments. For this, we used the KmR cassette from pK18mobsacB (Schäfer et al., 1994), the GmR cassette from pJQ200SK (Quandt and Hynes, 1993), the SmR/SpR cassette from pHP45Ω (Prentki and Krisch, 1984), or the TcR cassette from pHP45Ω-Tc (Fellay et al., 1987), the last two of which introduce transcriptional stops and/or translational stops in all reading frames. The constructs were transferred to P. syringae by either conjugation or electroporation, and marker exchange clones were selected as being resistant to sucrose and sensitive to the antibiotic resistance conferred by the vector. All mutations were confirmed by PCR and sequencing. The recipient strain ΔABCTc contains a deletion of 223 nt (positions 1216031..1216253 in accession no. NZ_CP008742) into the pseudogene PSA3335_RS06050 that was replaced by the TcR cassette from pHP45Ω-Tc. Other mutations are described in Supplementary Table S1.
Transposon Tn5-GDYN1 confers diverse antibiotic resistances and sensitivity to sucrose (Supplementary Table S1), mediated by gene sacB, allowing for the selection of derivatives cured of tagged plasmids by selection on media with sucrose (Flores et al., 1993). Mutagenesis with this transposon, selection of tagged plasmids and identification of the insertion point were done as described previously (Bardaji et al., 2019).
Conjugation conditions
Conditions for conjugation are as detailed below. In all cases, except when indicated, experiments were repeated at least three times each with four technical replicates, and means were compared by an analysis of variance (one-way ANOVA p > 0.05), followed by Duncan's multiple range test (p < 0.05). For the estimation of transfer frequencies, the detection limit was around 10−7 transconjugants per recipient. Routinely, we confirmed the transfer of plasmids by examining the plasmid profile by gel electrophoresis of at least 10 randomly chosen transconjugants.
Transfer on solid media
For a simple conjugation protocol, cells from liquid cultures in LB grown overnight with selecting antibiotics were collected, washed with one volume of ¼ strength Ringer's solution (¼ Ringer's; Oxoid, Basingstoke, UK) and adjusted to an optical density at 600 nm (OD600) of 1 in this buffer. Frequencies of transfer using this protocol could be very high, but were highly variable among experiments. For a standardized conjugation protocol, producing repetitive conjugation transfer frequencies, cells grown overnight in liquid LB with antibiotic selection were diluted in the same medium to an OD600 of 0.1, and incubated until they reached an OD600 of 0.25. Cultures were then collected, washed and resuspended in the same volume of ¼ Ringer's. After this, we proceed equally for both protocols. Equal volumes of suspensions of donor and recipient strains were mixed and 10 μL were deposited onto plates of LB, KMB or MAM, dried and incubated at 28°C. Each 10 μL spot was individually collected after 24 h (LB and KMB) or 48 h (MAM) and resuspended in ¼ Ringer's; then, 25 μL from 10-fold serial dilutions were spotted onto LB plates with appropriate antibiotics to estimate the transfer frequency. All liquid cultures were incubated at 28°C with shaking (200 rpm).
To evaluate the transfer of native plasmids to other pseudomonads, we used as donor strain FAM-098 (Supplementary Table S1), an arginine auxotroph derivative of strain NCPPB 3335, transformed with either pB::Tn5-GDYN1 or pC::Tn5-GDYN1. A loopful each of donor and recipient bacteria grown on LB plates overnight were mixed in 250 μL of ¼ Ringer's, and 10 μL of the suspension spotted onto MAM medium. After 48 h at 28°C, the conjugation mix was stripped onto the minimal medium SSM (Meyer and Abdallah, 1978) plus kanamycin to select for transconjugants. The transfer experiments were repeated twice, each with four technical replicates. Transfer of pB and pC was confirmed by PCR and the examination of the plasmid profile of at least two transconjugants per recipient.
Transfer on liquid media
Cells from liquid cultures in LB grown overnight with selecting antibiotics were collected, washed with one volume of ¼ Ringer's (Oxoid, Basingstoke, UK) and adjusted to an OD600 of 1 in this buffer. Suspensions containing equal volumes of donor and recipient strains were diluted 1:10 in liquid medium LB, KMB or MAM and grown in the dark at 28°C with shaking (200 rpm) until reaching an OD600 of 1–1.5. Then, transfer frequency was estimated from population counts by spotting 25 μL of serial 10-fold dilutions in ¼ Ringer's onto LB plus appropriate antibiotic selection.
Transfer on plant surfaces
Bean plants (Phaseolus vulgaris cv. Canadian Wonder) were grown in peat/perlite (2:1) until the primary leaves were fully expanded (around 10 d), in a growth chamber with a 16/8 h (light/dark) photoperiod at 22/18°C. Strains NCPPB 3335 and ΔABCTc were scraped off LB plates incubated overnight at 28°C, washed and resuspended in 10 mM MgCl2 to an OD600 of 0.1 (ca. 107 cfu mL−1). A 1:1 mix of both strains was sprayed onto both sides of primary leaves using an airbrush (VL Double Action; Paasche Airbrush Company, Kenosha, WI, USA) with an air pressure of 1 bar, until runoff; an experimental replica consisted of two inoculated plants and the experiment was repeated three times. Five days after inoculation, 2 leaf discs were excised per plant with a No. 4 cork borer (0.6 cm diameter), homogenized in 1 mL of ¼ Ringer's and serial dilutions plated onto LB with or without antibiotic selection for population counts.
Plasmid stability assays
The stability of pB in relation with its conjugative ability was estimated using strain ΔABC containing either pB::Tn5-GDYN1 or the non-conjugative variant pB (virB4::Km-sacB), containing a KmR-sacB cassette inactivating gene virB4. Overnight cultures of these strains on MAM plates were used to prepare suspensions adjusted to an OD600 of 0.1 (ca. 107 cfu mL−1), and 20 μL were deposited onto MAM plates. After 24 h incubation at 28°C, the bacterial growth was scraped off the plates and resuspended in 1 mL of liquid MAM, and 20 μL of these suspensions was deposited onto MAM plates. The strains were serially transferred in this way for a total of four transfers, and plasmid loss was examined from the starting cultures and after the first and fourth transfer. For this, serial dilutions of the bacterial suspensions were plated onto NA and 100 well-isolated colonies for each strain and time point were replica-plated onto NA, NA+Km, and SNA. The proportion of clones that had lost plasmid pB was estimated from the fraction of clones that were resistant to sucrose and sensitive to kanamycin. Loss of pB was confirmed in random samples of clones by plasmid profile analysis. The experiment was repeated three times with three replicates.
We tested long-term plasmid stability of plasmid pB on different bacteria as described previously (De Gelder et al., 2007; Bardaji et al., 2019). Newly obtained transconjugants of P. syringae pv. phaseolicola 1448A and P. viridiflava CECT458 containing pB::Tn5-GDYN1, which were maintained at −80°C, were grown overnight on LB plus kanamycin, and a single colony was used to start a culture in 3 mL of LB. After overnight growth, 10 μL of this culture were transferred to 3 mL of LB and incubated in the same conditions for 24 h. We estimated populations of ca. 107 cfu in the initial inoculum and of ca. 109 cfu after 24 h of incubation; therefore, we estimated that ~7 generations were obtained per growth cycle. After 21 serial transfers in this way, with a conservative estimate of at least 150 generations, cultures were plated on LB and 100 well-isolated colonies were replica-plated on LB plates with and without kanamycin, to estimate the proportion of plasmid-free cells. The experiment was repeated three independent times, each with four technical replicates. In each assay, we examined the plasmid profile of at least 10 random KmR colonies to ensure that the plasmids had not undergone any structural changes or become integrated in the chromosome.
Accession numbers
The accession numbers of the relaxases investigated here are: for pB, UniProtKB H1ZX26 (MobB, MOBP family); and for pC, UniProtKB H1ZX53 (MobB, MOBP family) and UniProtKB H1ZX69 (TraA, MOBQ family).
Results
Bioinformatics identification of conjugation genes
The sequences of plasmids A, B and C were analyzed with the software SecReT4 and oriTfinder, as well as by manual BLAST, allowing us to refine the previous annotation of genes involved in conjugation (Bardaji et al., 2011b), eliminating misidentifications and identifying overlooked gene fragments (Figure 1, Supplementary Table S3). As previously described (Bardaji et al., 2011b), we confirmed that pA contains a partial set of MPFI genes whereas pB contains a complete MPFT system (Supplementary Table S3). Furthermore, pB contains a gene coding for a putative type IA topoisomerase homologous to TraE from plasmid RP4 (Supplementary Table S3), which could participate in resolution of cointegrates during vegetative or conjugative replication (Li et al., 1997).
Figure 1. Relevant genes in the three virulence plasmids of Pseudomonas syringae pv. savastanoi NCPPB3335. Location and extent of genes involved in replication, virulence, synthesis of type IV secretion systems, and relaxases and associated RAPs. The location and extent of the origin of transfer associated to the traA genes is speculative and based on other members of the MOBQ family relaxases.
For the identification and classification of relaxases in the three virulence plasmids, we used the two previous programs and the MOBscan software (Figure 1, Supplementary Table S3), which offered partial and conflicting results. Alleles of gene mobB are present in both pB and pC (Figure 1, Supplementary Figure S1), and MOBscan identified their deduced products as relaxases of the MOBP family. Remarkably, SecReT4 and oriTfinder identified the mobB gene from pB as a virB8 conjugation gene but did not report the mobB gene from pC, despite the high similarity of the MobB products (86.6% aa identity; Supplementary Figure S1). A Pfam search of the MobB products found a Relaxase/Mobilization nuclease domain (PF03432), which Pfam reported significant for pC (E value 1.0e-6) but not for pB (E value 7.8e-3), and a significant (E value ≤ 1.7e-14) Large polyvalent protein-associated domain (PF18821), which is widely present in proteins transmitted by conjugation (Iyer et al., 2017). The structure of the N-terminal residues of both MobB products could be modeled by Phyre2 with more than 90% confidence, being similar to diverse relaxase domains and a replication initiator protein. This is not surprising because it is described that relaxases share motifs with replication proteins (Guglielmini et al., 2014; Guzmán-Herrador and Llosa, 2019). We therefore considered that the mobB genes from pB and pC coded for bona fide relaxases of the MOBP family.
The three programs also predicted a second relaxase gene, gene traA, in both plasmids, coding for a MOBQ family relaxase, except that the gene in pB (PSA3335_RS00025) is truncated and likely not functional (Supplementary Figure S1, Supplementary Table S3).
oriTfinder identified the oriT sequences accompanying the MOBP relaxase genes in pB and pC but not a putative oriT associated to the MOBQ relaxase genes (Supplementary Figure S2). None of the programs identified any of the mobC RAP genes accompanying the relaxase genes; at least for the MOBP RAPs this could be due to their low sequence conservation (Supplementary Figure S1B). No relaxase genes nor oriT were found in plasmid pA.
Previous research demonstrated the chromosomal integration of native plasmids in P. syringae strains (Szabo and Mills, 1984; Neale et al., 2013). We therefore examined the possible presence of relaxases in the chromosome of strain NCPPB 3335, where we have previously found remnants of plasmid genomes (J. Murillo and C. Ramos, unpublished data). A blastp search using the pB and pC relaxases as query against the genome of strain NCPPB 3335 (accession no. NZ_CP008742), as well as a simple text search of its annotation, identified four CDSs annotated as relaxases or TraI domain-containing protein, which were examined with MOBscan. Two of them, PSA3335_RS05310 and PSA3335_RS09015, belong to the MOBP family and are preceded by a RAP gene. In blastp comparisons, PSA3335_RS05310 and PSA3335_RS09015 showed 30% identity with around 80% coverage, whereas only PSA3335_RS09015 showed detectable sequence homology with the MOBP relaxases from plasmids B and C, returning alignments with around 29% identity over 157 amino acids or less. CDS PSA3335_RS07700 is a pseudogene and PSA3335_RS23785 belongs to the MOBH family and codes for a putative relaxase from an integrating conjugative element.
Two virulence plasmids are mobilizable
We evaluated the conjugational transfer of the three virulence plasmids from strain NCPPB 3335 using diverse strains as donors. These strains contained combinations of one to three plasmids with only one of them tagged with transposon Tn5-GDYN1 (KmR; Supplementary Table S1). As the recipient, we routinely used strain ΔABCTc, a derivative of strain NCPPB 3335 cured of the three plasmids and tagged with a tetracycline-resistance (TcR) cassette (Supplementary Table S1). Using a simple conjugation protocol, we initially observed conjugation frequencies of pB as high as 3.7 ± 0.8 × 10−1 transconjugants per recipient, but with experiment-to-experiment variations of two orders of magnitude or more. We therefore used a standardized conjugation protocol (see materials and methods), which resulted in lower conjugation frequencies but provided consistent results in different experiments.
In standard mating experiments on solid media, we did not detect the transfer of pA (Table 1). However, using donors with a tagged pA in matings with the simple conjugation protocol, which allows for higher transfer frequencies, we occasionally obtained one or two TcR-KmR transconjugants per mating (transfer frequency ≤ 0.8 ± 0.2 × 10−6 transconjugants per recipient). We examined the plasmid profiles of 20 randomly chosen transconjugants to confirm plasmid transfer (Figure 2). All clones suspected to have received pA contained a plasmid of similar or larger size than pA but that hybridized to a pB-specific probe (Figure 2). One of the clones was apparently plasmidless and did not hybridize with the pB probe; we therefore assumed that during the mating process this clone had received a copy of Tn5-GDYN1, which is functional for transposition. By PCR with specific primers, we determined that four of these 20 clones contained both virulence genes ptz (from pA) and hopAO1 (from pB), whereas some other clones contained assorted combinations of maintenance genes from pA and pB (data not shown). This suggests that chimeric molecules between plasmids pA and pB might be formed and subsequently mobilized with very low frequency, and that pA is not mobilizable in its native form. Therefore, we did not examine the transfer of pA in any further experiments.
Table 1. Frequency of transfer of virulence plasmids from P. syringae pv. savastanoi NCPPB 3335 in solid and liquid media.
Figure 2. Plasmid profiles of transconjugants and their hybridization to a probe specific for pB. Conjugation was carried out on solid minimal medium A using the plasmidless strain ΔABCTc as recipient and as donors the following strains: (A) strain NCPPB 3335 containing pPsv48A::Tn5-GDYN1 (pA*); (B) strain ΔABC containing pPsv48B::Tn5-GDYN1 (pB*); (C) strain NCPPB 3335 containing pPsv48C::Tn5-GDYN1 (pC*). The bottom picture in (A,C) correspond to the gels shown above after being subjected to Southern hybridization using as probe a region from pB (positions 14062..15310, accession no. FR820586) not repeated elsewhere in the genome of strain NCPPB 3335. WT, wild type strain NCPPB 3335; D, donor strain; R, receptor strain ΔABCTc; C+, hybridization control (strain ΔABC with pB*); clp, chromosome and linear plasmids.
Conversely, we routinely obtained TcR-KmR transconjugants with tagged variants of pB and pC (Table 1), indicating that these native plasmids can be mobilized. The frequency of conjugation was high for pB, reaching over 10−2 transconjugants per recipient, whereas the transfer frequency of pC was always at least one order of magnitude lower. Transfer of pC only occurred when the donor also contained pB (Table 1 and data not shown). This indicates that pB contains all the machinery for its conjugative transfer, whereas pC is mobilizable.
Plasmid transfer occurred both on solid and in liquid media, although frequency of transfer was significantly higher on solid media (Table 1). Additionally, transfer was highest on minimal medium (MAM), whereas on rich media transfer was consistently higher on KMB than on LB (Table 1 and data not shown). Therefore, all matings were hereafter carried out on MAM plates.
We routinely examined the plasmid profiles of selected transconjugants to confirm plasmid transfer. All examined transconjugants from matings with the tagged plasmids B and C contained the plasmid band of the corresponding size, indicating that these plasmids can be efficiently mobilized (Figure 2). However, a sizable portion of the transconjugants selected for the acquisition of pC::Tn5-GDYN1 also contained pB. To estimate the frequency of co-transfer, we used as donor strain NCPPB 3335 (pC::Tn5-GDYN1) (KmR) and the plasmidless strain ΔABCTc (TcR) as recipient in mating experiments. The plasmid profiles of 300 KmR-TcR transconjugants from three independent experiments were analyzed visually and by Southern hybridization with a probe specific for plasmid pB (Figure 2 and data not shown). This showed that an average of 80% of the transconjugants had also acquired pB, indicating a frequent co-transfer of plasmid pB when pC is mobilized. Taking into account the data in Table 1, this may also imply that when pB is transferred by conjugation, pC is also transferred between 2 and 6% of the time.
The predicted oriT from pB and pC associated to MOBP relaxases are functional
Our previous annotation (Bardaji et al., 2011b) and the software oriTfinder predicted the presence of an oriT associated to the MOBP-family relaxase gene in both pB (107 nt) and pC (146 nt) (Supplementary Figure S2). To evaluate their functionality, fragments of 221 (pB) or 220 nt (pC) containing the predicted oriT (Supplementary Figure S2) were cloned into the non-mobilizable plasmid pKMAGJ and electroporated into strain NCPPB 3335. Both cloned fragments were transferred to strain ΔABCTc with average frequencies of 2.2 ± 1.5 × 10−5, for oriT-pB, and 3.3 ± 0.4 × 10−4 transconjugants per recipient for oriT-pC. These results indicate that the cloned sequences contain functional oriT sequences. Their lower transfer frequency compared to those of the whole plasmids (Table 1) suggest that (1) they might be missing sequences important for the transfer process, (2) that they compete inefficiently for the relaxase with the native copy of the oriT, and/or (3) that the MobB relaxases might act more efficiently in cis for their cognate oriT.
Despite repeated efforts, we were not able to clone the oriT associated to traA in pC for functional testing.
The MPFT from pB is essential for the mobilization of the virulence plasmids
Gene virB4 is part of the MPFT from pB and codes for a highly conserved protein widely present in different MPF clusters (a cytoplasmic ATPase providing energy for pilus assembly), thus being essential for conjugation (Peña et al., 2012). We constructed a derivative of pB containing a Gm-resistance cassette replacing the last 1092 nt of the virB4 CDS, causing the deletion of 42.8% of the C-terminal end of the VirB4 protein, resulting in plasmid pB(virB4::Gm) (Supplementary Tables S1, S2). Nevertheless, and since virB4 is part of a large operon, it is feasible that this mutation could exert polar effects and alter the expression of other downstream genes involved in conjugation.
In mating experiments with strain ΔABCTc as recipient, a strain containing pB (virB4::Gm) did not transfer either this plasmid or pC::sacB at detectable levels. These results then indicate that the mobilization of the virulence plasmids of strains NCPPB 3335 depends entirely on the functionality of the MPFT from pB.
Mobilization of pB does not appear to contribute to plasmid stability
Low-copy number plasmids generally contain diverse maintenance determinants ensuring their stable vertical inheritance (plasmid stability) (Nordström and Austin, 1989). However, our previous analyses did not find any relevant maintenance determinants in pB(Bardaji et al., 2019). Therefore, we evaluated the possibility that conjugative transfer might contribute to the stable inheritance of pB, as was described for other plasmids (De Gelder et al., 2007). Since the sacB gene causes lethality in the presence of sucrose, we could estimate the loss rate of pB::Tn5-GDYN1 or of the non-conjugative variant pB(virB4::Km-sacB), both containing sacB, from the proportion of sucrose-resistant colonies growing on sucrose-containing media, as previously described (Flores et al., 1993; Bardaji et al., 2019). To avoid false negatives, from counting those cells containing a plasmid with a spontaneous mutation inactivating gene sacB, we only considered that the plasmid was lost when the cells were resistant to sucrose (SucR) and sensitive to kanamycin (KmS). To also eliminate any possible interference, we evaluated the loss of pB in strains containing only this plasmid.
After one day of growth on solid MAM, the fraction of SucR-KmS cells (around 10−6 cfu mL−1; Figure 3) was not significantly different between the cultures containing pB::Tn5-GDYN1 or of the non-conjugative variant pB(virB4::Km-sacB). We obtained similar results after 4 days of daily transfers of each culture on solid fresh medium (Figure 3). In our particular experimental conditions, therefore, the conjugative ability of pB does not significantly contribute to the maintenance of this plasmid in the bacterial population.
Figure 3. Effect of the conjugative ability of pB (pPsv48B) on plasmid loss. Colored bars indicate the frequency of plasmid loss (mean ± sd), measured as the number of sucrose-resistant and kanamycin-sensitive colonies per total colony forming units in unsupplemented medium. The plasmidless strain P. syringae pv. savastanoi ΔABC containing tagged derivatives of pB that were either conjugative (pB::Tn5-GDYN1, dark orange) or non-conjugative (pB[virB4::Km-sacB], light orange) were subcultured every 24 h on solid medium MAM. Values are the average of at least three independent experiments, each with three replicates. Means with different letters are significantly different (one-way ANOVA and Duncan's multiple range test; p < 0.05).
Plasmids pB and pC can be mobilized by two to three different relaxases
Plasmids pB and pC each contain a MOBP relaxase gene, whereas pC contains an additional full-length MOBQ relaxase gene. To explore the role of these three relaxases in conjugative transfer, we constructed diverse strains containing insertions of an antibiotic resistance cassette in one or more of the three relaxase genes. As donors in mating experiments we used strains containing plasmids pB and pC, with appropriate relaxase mutations, and the plasmidless strain ΔABCTc as a recipient (Figure 4).
Figure 4. Role of the three relaxase genes in the conjugal mobilization of pB and pC evidenced by site-directed mutagenesis. Bars indicate the mobilization frequency (mean ± sd) of plasmids pB (orange) and pC (green) in strains containing functional (+) or mutated (–) versions of the three relaxase genes, as indicated in the panel below the graph. x, no transconjugants detected. Values are the average of at least three independent experiments, each with four replicates. Means with different letters are significantly different (one-way ANOVA and Duncan's multiple range test; p < 0.05).
The two MOBP relaxases were interchangeable and could mobilize both plasmids pB and pC (Figure 4). Plasmid pC was mobilized with the same low frequency by both relaxases; however, they were not equally efficient for the mobilization of pB, with a drop of at least two orders of magnitude in the frequency of transfer with the relaxase from pC (Figure 4).
The MOBQ relaxase mobilized only plasmid pC, and with a frequency of transfer that was not significantly different to those of the MOBP relaxases. The lack of mobilization of pB could be due to the fact that pB contains only a partial putative oriT associated to its truncated MOBQ relaxase gene (Supplementary Figure S2). Additionally, it is likely that the MOBQ relaxase show a cis-acting preference, as is typical of this protein family (Garcillán-Barcia et al., 2019), which could reduce the potential frequency of transfer of pB to undetectable levels.
Remarkably, the coexistence of different relaxases had a subtle but significant effect, with different combinations of functional relaxases increasing or reducing the transfer frequency of the individual plasmids (Figure 4). For instance, the frequency of transfer of pC was maximized with only two functional relaxase genes when the MOBP gene in pB was mutated.
Finally, we could not detect the mobilization of plasmids B or C in strains containing mutated versions of the three relaxase genes (Figure 4), even if the donor strain also contained an intact copy of plasmid pA (data not shown). This data indicates that the conjugative transfer of plasmids B and C is exclusively mediated by their own relaxase genes.
Plasmid pB also transfers among epiphytic populations
We wanted to evaluate the possibility that the virulence plasmids of NCPPB 3335 could be transferred among bacteria on plant surfaces. We decided to use bean plants for these experiments because of the advantages of using a herbaceous plant and because bean is a non-host plant for strain NCPPB 3335; additionally, strain NCPPB 3335 is phylogenetically close to different pathovars infecting bean.
To test this, we spray inoculated bacterial suspensions containing a 1:1 mixture (~2.5 × 107 cfu mL−1 each) of strains ΔABC(pB::Tn5-GDYN1) and ΔABCTc onto bean leaves. The strains survived only poorly on bean leaves, and at 5 dpi we could find bacterial populations only in two out of the six inoculated plants, ranging from 0.05 × 10−2 to 3 × 10−2 cfu leaf disc−1. Nevertheless, in these plants we observed the transfer of pB at a frequency of around 2 × 10−2 transconjugants per recipient, demonstrating that plasmid pB can be transferred in natural epiphytic populations. Since bean is a non-host plant, and because P. syringae strains can persist on many different plants and other reservoirs (Morris et al., 2019), this suggests that strain NCPPB 3335 could potentially transfer its virulence plasmids to a wide bacterial population in different environments.
pB mobilizes to different bacteria
Plasmids, like other MGEs, are usually distributed preferentially to phylogenetically related bacteria (Smillie et al., 2010; Jackson et al., 2011; Bardaji et al., 2017b). To test the transfer of pB to other bacteria, we transformed pB::Tn5-GDYN1 in an NCPPB 3335 derivative auxotrophic for arginine, strain FAM-098 (Matas et al., 2012), to be used as donor. As recipients, we tested 46 strains from eight phylogroups of P. syringae (Dillon et al., 2019) and other Pseudomonas species, and transconjugants were selected on SSM plus kanamycin.
Plasmid pB was efficiently transferred to 22 strains belonging to six phylogroups of P. syringae (Table 2 and Supplementary Table S4). However, we did not detect its transfer to any of the remaining strains tested. We evaluated the stability of pB in two recipient strains chosen as examples, P. syringae pv. phaseolicola 1448A and P. viridiflava CECT458. After 21 serial transfers on liquid LB medium, where pB is not transferred by conjugation (Table 1), we did not observe the generation of plasmid-free segregants. These results indicate that pB can potentially be distributed widely at least within P. syringae, and that it would be stably maintained once acquired.
Table 2. Conjugative transfer of plasmids pB and pC from P. syringae pv. savastanoi NCPPB 3335 to different strains and species of Pseudomonasa.
Conversely, we observed the transfer of pC to only two strains out of the 20 tested, and both from the same phylogroup that includes the plasmid parent strain P. syringae pv. savastanoi NCPPB 3335.
Discussion
P. syringae pv. savastanoi NCPPB 3335 contains three virulence plasmids: plasmid pA contains the cytokinin biosynthesis gene ptz and the type III effector hopAF1, pB contains the type III effector gene hopAO1, and pC carries the cytokinin biosynthesis gene idi (Bardaji et al., 2011b; Castañeda-Ojeda et al., 2017b; Añorga et al., 2020). Here we showed that the two small virulence plasmids, pB and pC, are transferable horizontally by conjugation in free living conditions, both in solid and liquid media, and on epiphytic populations of a non-host plant. Additionally, we observed that fragments of the large virulence plasmid pA can also be transferred at very low frequency after being cointegrated with plasmid pB. Therefore, the plasmid-borne virulence genes of P. syringae pv. savastanoi can be efficiently transferred horizontally by conjugation.
As occurs with other P. syringae strains (Zhao et al., 2005), native plasmids of strain NCPPB 3335 contained a variety of MPF systems, with an incomplete version of an MPFI (in pA) and a complete MPFT (in pB). Not surprisingly, the incomplete system of plasmid pA appears to be non-functional and mobilization of all plasmids of strain NCPPB 3335 was dependent on the complete MPFT system of pB. The evaluation of specific knockout mutants (Figure 4) also showed that the MPFT system and the T4C protein VirD4 from pB are promiscuous, being able to efficiently recognize both the MOBP and MOBQ relaxosomes. This is important, because not all MPFT systems can mobilize plasmids with MOBQ relaxases (Garcillán-Barcia et al., 2019), and suggests that these plasmids have undergone selection for the promiscuous use of conjugation systems to maximize their chances for horizontal distribution.
The native plasmids of strain NCPPB 3335 contained three functional plasmid-borne relaxase genes (Figures 1, 4), one located in pB and two in pC. The MOBP relaxases found in pB and pC were functionally interchangeable, although with differing efficiency that was determined both by the relaxase and the oriT (Figure 4). Thus, pB was transferred with high frequency by its own MOBP relaxase but with a frequency at least two orders of magnitude lower by the MOBP relaxase from pC, suggesting a high specificity with its cognate oriT. Conversely, pC was always transferred with low frequency by both MOBP relaxases. The third relaxase, from the MOBQ family, only allowed the transfer of pC, which was expected because pB did not appear to contain a cognate functional origin of transfer. Remarkably, our data indicate that different combinations of functional relaxases can impact the transfer ability of the native plasmids both positively and negatively (Figure 4).
We also found three other relaxase genes in the chromosome of strain NCPPB 3335, two MOBP and one MOBH. However, they do not seem to contribute to the mobilization of the native plasmids, because derivative strains containing mutations in all three plasmidic relaxases were non-conjugative (Figure 4). This is unexpected because MOBP relaxases can act in trans (Guzmán-Herrador and Llosa, 2019), and therefore it might indicate that the chromosomal relaxases preferentially act in cis or also in trans but with a very low frequency, or that they are non-functional.
The location of two relaxase genes in plasmid pC is mirrored by its possession of two functional replication initiation protein (Rep) genes (Bardaji et al., 2017a). Therefore, and because Rep genes are essential for plasmid replication, this suggest that plasmid pC originated from the cointegration of two previously independent plasmids. The presence of a truncated traA gene in pB, coding for the MOBQ relaxase, indicate that this plasmid might have also originated from the fusion of two plasmids. The reason for the occurrence of multiple relaxase genes in these plasmids is intriguing, because plasmids with more than one replicon are common but only <2.6% of the 11386 plasmids analyzed in a previous work contain more than one relaxase gene (Coluzzi et al., 2022). It is possible then that the occurrence of these relaxases does not impart any selective advantage to strain NCPPB 3335 or its plasmids, but that it is only a neutral side effect of the fortuitous cointegration of plasmids to generate a surviving chimeric molecule. Nevertheless, we have previously shown that the presence of an additional replicon in pC contributes to a very high plasmid stability (Bardaji et al., 2019). Additionally, the data presented here support the hypothesis that the occurrence of multiple relaxase genes likely increases the chances for the horizontal transfer of virulence genes.
One possible advantage of carrying a duplicate set of determinants for replication and horizontal transfer is to maximize plasmid carriage and distribution, because of the availability of spare functional copies in the event of gene inactivation. Plasmids of P. syringae are subjected to repeated injury from mobile genetic elements, which could accumulate to close to a third of the plasmid sequences and can lead to gene inactivation, sequence reorganizations and internal deletions (Bardaji et al., 2011a,b). In particular, the ubiquitous insertion sequence IS801 and diverse miniature transposable elements are the cause of frequent deletions and reorganizations of native plasmids in the P. syringae complex (Bardaji et al., 2011a, 2019). Additionally, incoming conjugative plasmids might differentially recognize the oriT-relaxase arrays of the native plasmids of strain NCPPB 3335 and facilitate their mobilization to other bacterial hosts, possibly even with a higher conjugation frequency than that achieved by the wild type strain. It is also possible that the occurrence of multiple relaxases is favored because they might impact the bacterial life cycle in other ways, instead of or in addition to their ability to transfer horizontally. In fact, relaxases can play several other functions besides their participation in conjugative transfer of their cognate oriT (Guzmán-Herrador and Llosa, 2019); thus, MOBP and MOBQ relaxases can have in trans activity on heterologous oriT and act as plasmid recombination enzymes. In particular, plasmid recombination may take place between related oriT sequences and lead to the formation or separation of cointegrates between two independent replicons (Wawrzyniak et al., 2017).
Conjugative transfer was shown to increase the stability of diverse plasmids in their native bacterial host through mobilization from plasmid-bearing to plasmid-free cells (De Gelder et al., 2007). In our conditions, however, the high rates of conjugative transfer of pB do not appear to contribute to plasmid stability. Therefore, it is likely that the main purpose of the mobilization system of strain NCPPB 3335 is the horizontal dissemination of its native plasmids. In fact, our in vitro experiments showed that pB can be transferred to different species of pseudomonads, all of which belong to P. syringae sensu lato. Additionally, pC was transferred to only two of the 20 bacteria assayed, and both belonging to the same genomospecies as strain NCPPB 3335. This suggests that the horizontal transfer of the native plasmids of strains NCPPB 3335 might be restricted or favored to closely related bacteria, as occurs with other mobile elements within the P. syringae complex or native plasmids from other bacteria (Smillie et al., 2010; Jackson et al., 2011; Bardaji et al., 2017b).
The acquisition of DNA by horizontal transfer has been revealed as a main factor leading to the diversification of P. syringae lineages (Guttman et al., 2008; Nowell et al., 2014; Dillon et al., 2019). In fact, previous works revealed incongruences between the phylogeny of diverse native plasmids and their P. syringae hosts (Ma et al., 2007; Bardaji et al., 2011b, 2017a; Gutiérrez-Barranquero et al., 2017), illustrating the horizontal exchange of plasmids among different pathovars. Although it is likely that this exchange was mediated by conjugation, the frequency and extent of plasmid transfer among P. syringae populations is unknown. Surprisingly, the conjugative transfer of native plasmids has been demonstrated in only a few strains of P. syringae pathovars actinidiae (Zhao et al., 2018), glycinea (Watanabe et al., 1998), papulans (Burr et al., 1988), syringae (Sundin and Bender, 1993; Cazorla et al., 2002), and tomato (Bender and Cooksey, 1986), and nearly always related to resistance to antibacterial compounds. Nevertheless, many native plasmids of P. syringae carry genes coding for an MPF. A hybridization analysis of 31 native plasmids of P. syringae showed that most of them contained genes for at least one MPF, the majority of the class MPFT (previously TIVASS), although in many cases the gene clusters were incomplete (Zhao et al., 2005). The analysis of completely sequenced plasmids also supported this trend (Buell et al., 2003; Stavrinides and Guttman, 2004; Joardar et al., 2005; Bardaji et al., 2011b; Gutiérrez-Barranquero et al., 2017). These results make it likely that most native plasmids of P. syringae are equipped for their horizontal exchange. However, we should keep in mind that these MPFs could participate in other cellular functions instead of conjugation, such as the secretion of DNA, effectors or toxins to eukaryotic cells (Wallden et al., 2010; Souza et al., 2015). More importantly, mobilization of plasmid DNA is only possible for molecules containing an oriT and in the presence of a compatible relaxase, and there is no information on the abundance and distribution of origins of transfer and relaxases among P. syringae plasmids. To estimate the potential for horizontal gene transfer, future works should strive to evaluate the distribution of oriT and relaxases among native plasmids of strains of the P. syringae complex.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
JM and CR conceived the study and designed the experiments. MA and MU performed the experiments. MA, MU, CR, and JM analyzed the data and interpreted the results. JM drafted the manuscript with contributions from MA, MU, and CR. All authors read and approved the final manuscript.
Funding
This work was funded by the Spanish Plan Nacional I+D+i grants PID2020-115177RB-C21 and PID2020-115177RB-C22 financed by the Ministerio de Ciencia e Innovación and the Agencia Estatal de Investigación, MCIN/AEI/10.13039/501100011033/, Spain. The funders had no role in study design, data collection, and interpretation or the decision to submit the work for publication.
Acknowledgments
We are grateful to Theresa H. Osinga for revising the manuscript and help with English usage. We are also indebted to F. Cazorla, M. M. López, E. López-Solanilla, and S. Tegli for their kind donation of bacterial strains.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1076710/full#supplementary-material
References
Añorga, M., Pintado, A., Ramos, C., De Diego, N., Ugena, L., Novák, O., et al. (2020). Genes ptz and idi, coding for cytokinin biosynthesis enzymes, are essential for tumorigenesis and in planta growth by P. syringae pv. savastanoi NCPPB 3335. Front. Plant Sci. 11, 1294. doi: 10.3389/fpls.2020.01294
Backert, S., and Grohmann, E., (eds.). (2017). Type IV Secretion in Gram-Negative and Gram-Positive Bacteria. Cham: Springer International Publishing.
Baltrus, D. A., Nishimura, M. T., Romanchuk, A., Chang, J. H., Mukhtar, M. S., Cherkis, K., et al. (2011). Dynamic evolution of pathogenicity revealed by sequencing and comparative genomics of 19 Pseudomonas syringae isolates. PLoS Pathog. 7, e1002132. doi: 10.1371/journal.ppat.1002132
Bardaji, L., Añorga, M., Echeverría, M., Ramos, C., and Murillo, J. (2019). The toxic guardians - multiple toxin-antitoxin systems provide stability, avoid deletions and maintain virulence genes of Pseudomonas syringae virulence plasmids. Mobile DNA 10, 7. doi: 10.1186/s13100-019-0149-4
Bardaji, L., Añorga, M., Jackson, R. W., Martínez-Bilbao, A., Yanguas, N., and Murillo, J. (2011a). Miniature transposable sequences are frequently mobilized in the bacterial plant pathogen Pseudomonas syringae pv. phaseolicola. PLoS ONE. 6, e25773. doi: 10.1371/journal.pone.0025773
Bardaji, L., Añorga, M., Ruiz-Masó, J. A., del Solar, G., and Murillo, J. (2017a). Plasmid replicons from Pseudomonas are natural chimeras of functional, exchangeable modules. Front. Microbiol. 8, 190. doi: 10.3389/fmicb.2017.00190
Bardaji, L., Echeverría, M., Rodríguez-Palenzuela, P., Martínez-García, P. M., and Murillo, J. (2017b). Four genes essential for recombination define GInts, a new type of mobile genomic island widespread in bacteria. Sci. Rep. 7, 46254. doi: 10.1038/srep46254
Bardaji, L., Pérez-Martínez, I., Rodríguez-Moreno, L., Rodríguez-Palenzuela, P., Sundin, G. W., Ramos, C., et al. (2011b). Sequence and role in virulence of the three plasmid complement of the model tumor-inducing bacterium Pseudomonas savastanoi pv. savastanoi NCPPB 3335. PLoS ONE 6, e25705. doi: 10.1371/journal.pone.0025705
Bender, C., and Cooksey, D. (1986). Indigenous plasmids in Pseudomonas syringae pv. tomato: conjugative transfer and role in copper resistance. J. Bacteriol. 165, 534–541. doi: 10.1128/jb.165.2.534-541.1986
Bi, D., Liu, L., Tai, C., Deng, Z., Rajakumar, K., and Ou, H.-Y. (2013). SecReT4: a web-based bacterial type IV secretion system resource. Nucleic Acids Res. 41, D660–D665. doi: 10.1093/nar/gks1248
Brockhurst, M. A., Harrison, E., Hall, J. P. J., Richards, T., McNally, A., and MacLean, C. (2019). The ecology and evolution of pangenomes. Curr. Biol. 29, R1094–R1103. doi: 10.1016/j.cub.2019.08.012
Buell, C. R., Joardar, V., Lindeberg, M., Selengut, J., Paulsen, I. T., Gwinn, M. L., et al. (2003). The complete genome sequence of the Arabidopsis and tomato pathogen Pseudomonas syringae pv. tomato DC3000. Proc. Natl. Acad. Sci. U.S.A. 100, 10181–10186. doi: 10.1073/pnas.1731982100
Bull, C. T., and Koike, S. T. (2015). Practical benefits of knowing the enemy: modern molecular tools for diagnosing the etiology of bacterial diseases and understanding the taxonomy and diversity of plant-pathogenic bacteria. Annu. Rev. Phytopathol. 53, 157–180. doi: 10.1146/annurev-phyto-080614-120122
Burr, T. J., Norelli, J. L., Katz, B., Wilcox, W. F., and Hoying, S. A. (1988). Streptomycin resistance of Pseudomonas syringae pv. papulans in apple orchards and its association with a conjugative plasmid. Phytopathology 78, 410–413. doi: 10.1094/Phyto-78-410
Carver, T., Berriman, M., Tivey, A., Patel, C., Bohme, U., Barrell, B. G., et al. (2008). Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics 24, 2672–2676. doi: 10.1093/bioinformatics/btn529
Castañeda-Ojeda, M. P., López-Solanilla, E., and Ramos, C. (2017a). Differential modulation of plant immune responses by diverse members of the Pseudomonas savastanoi pv. savastanoi HopAF type III effector family. Mol. Plant Pathol. 18, 625–634. doi: 10.1111/mpp.12420
Castañeda-Ojeda, M. P., Moreno-Pérez, A., Ramos, C., and López-Solanilla, E. (2017b). Suppression of plant immune responses by the Pseudomonas savastanoi pv. savastanoi NCPPB 3335 type III effector tyrosine phosphatases HopAO1 and HopAO2. Front. Plant Sci. 8, 680. doi: 10.3389/fpls.2017.00680
Cazorla, F. M., Arrebola, E., Sesma, A., Pérez-García, A., Codina, J. C., Murillo, J., et al. (2002). Copper resistance in Pseudomonas syringae strains isolated from mango is encoded mainly by plasmids. Phytopathology 92, 909–916. doi: 10.1094/PHYTO.2002.92.8.909
Choi, K. H., Kumar, A., and Schweizer, H. P. (2006). A 10-min method for preparation of highly electrocompetent Pseudomonas aeruginosa cells: application for DNA fragment transfer between chromosomes and plasmid transformation. J. Microbiol. Meth. 64, 391–397. doi: 10.1016/j.mimet.2005.06.001
Coluzzi, C., Garcillán-Barcia, M. P., de la Cruz, F., and Rocha, E. P. C. (2022). Evolution of plasmid mobility: origin and fate of conjugative and nonconjugative plasmids. Mol. Biol. Evol. 39, msac115. doi: 10.1093/molbev/msac115
De Gelder, L., Ponciano, J. M., Joyce, P., and Top, E. M. (2007). Stability of a promiscuous plasmid in different hosts: no guarantee for a long-term relationship. Microbiology 153, 452–463. doi: 10.1099/mic.0.2006/001784-0
de la Cruz, F., Frost, L. S., Meyer, R. J., and Zechner, E. L. (2010). Conjugative DNA metabolism in Gram-negative bacteria. FEMS Microbiol. Rev. 34, 18–40. doi: 10.1111/j.1574-6976.2009.00195.x
Dillon, M. M., Thakur, S., Almeida, R. N. D., Wang, P. W., Weir, B. S., and Guttman, D. S. (2019). Recombination of ecologically and evolutionarily significant loci maintains genetic cohesion in the Pseudomonas syringae species complex. Genome Biol. 20, 3. doi: 10.1186/s13059-018-1606-y
Fellay, R., Frey, J., and Krisch, H. (1987). Interposon mutagenesis of soil and water bacteria: a family of DNA fragments designed for in vitro insertional mutagenesis of Gram-negative bacteria. Gene 52, 147–154. doi: 10.1016/0378-1119(87)90041-2
Firrao, G., Torelli, E., Polano, C., Ferrante, P., Ferrini, F., Martini, M., et al. (2018). Genomic structural variations affecting virulence during clonal expansion of Pseudomonas syringae pv. actinidiae biovar 3 in Europe. Front. Microbiol. 9, 656. doi: 10.3389/fmicb.2018.00656
Flores, M., Brom, S., Stepkowski, T., Girard, M. L., Dávila, G., Romero, D., et al. (1993). Gene amplification in Rhizobium: identification and in vivo cloning of discrete amplifiable DNA regions (amplicons) from Rhizobium leguminosarum biovar phaseoli. Proc. Natl. Acad. Sci. U.S.A. 90, 4932–4936. doi: 10.1073/pnas.90.11.4932
Garcillán-Barcia, M. P., Cuartas-Lanza, R., Cuevas, A., and de la Cruz, F. (2019). cis-acting relaxases guarantee independent mobilization of MOBQ4 plasmids. Front. Microbiol. 10, 2557. doi: 10.3389/fmicb.2019.02557
Garcillán-Barcia, M. P., Redondo-Salvo, S., Vielva, L., and de la Cruz, F. (2020). “MOBscan: automated annotation of MOB relaxases,” in Horizontal Gene Transfer. Methods in Molecular Biology, Vol. 2075, ed F. de la Cruz (New York, NY: Humana), p. 295–308. doi: 10.1007/978-1-4939-9877-7_21
Gardan, L., Shafik, H., Belouin, S., Broch, R., Grimont, F., and Grimont, P. A. D. (1999). DNA relatedness among the pathovars of Pseudomonas syringae and description of Pseudomonas tremae sp. nov. and Pseudomonas cannabina sp. nov. (ex Sutic and Dowson 1959). Int. J. Syst. Bacteriol. 49, 469–478. doi: 10.1099/00207713-49-2-469
Gomila, M., Busquets, A., Mulet, M., García-Valdés, E., and Lalucat, J. (2017). Clarification of taxonomic status within the Pseudomonas syringae species group based on a phylogenomic analysis. Front. Microbiol. 8, 2422. doi: 10.3389/fmicb.2017.02422
Guglielmini, J., Néron, B., Abby, S. S., Garcillán-Barcia, M. P., la Cruz, F. d., et al. (2014). Key components of the eight classes of type IV secretion systems involved in bacterial conjugation or protein secretion. Nucleic Acids Res. 42, 5715–5727. doi: 10.1093/nar/gku194
Gutiérrez-Barranquero, J. A., Cazorla, F. M., de Vicente, A., and Sundin, G. W. (2017). Complete sequence and comparative genomic analysis of eight native Pseudomonas syringae plasmids belonging to the pPT23A family. BMC Genom. 18, 365. doi: 10.1186/s12864-017-3763-x
Guttman, D. S., Morgan, R. L., and Wang, P. W. (2008). “The evolution of the pseudomonads,” in Pseudomonas syringae Pathovars and Related Pathogens - Identification, Epidemiology and Genomics, eds M. B. Fatmi, A. Collmer, N. S. Iacobellis, J. W. Mansfield, J. Murillo, N. W. Schaad, et al. (Dordrecht: Springer Netherlands), 307–319.
Guzmán-Herrador, D. L., and Llosa, M. (2019). The secret life of conjugative relaxases. Plasmid 104, 102415. doi: 10.1016/j.plasmid.2019.102415
Halary, S., Leigh, J. W., Cheaib, B., Lopez, P., and Bapteste, E. (2010). Network analyses structure genetic diversity in independent genetic worlds. Proc. Natl. Acad. Sci. U. S. A. 107, 127–132. doi: 10.1073/pnas.0908978107
Holmes, D. S., and Quigley, M. (1981). A rapid boiling method for the preparation of bacterial plasmids. Anal. Biochem. 114, 193–197. doi: 10.1016/0003-2697(81)90473-5
Hubbard, T. J., Aken, B. L., Ayling, S., Ballester, B., Beal, K., Bragin, E., et al. (2008). Ensembl 2009. Nucleic Acids Res. 37, D690–D697. doi: 10.1093/nar/gkn828
Iyer, L. M., Burroughs, A. M., Anand, S., de Souza, R. F., and Aravind, L. (2017). Polyvalent proteins, a pervasive theme in the intergenomic biological conflicts of bacteriophages and conjugative elements. J. Bacteriol. 199, e00245-17. doi: 10.1128/JB.00245-17
Jackson, R. W., Vinatzer, B., Arnold, D. L., Dorus, S., and Murillo, J. (2011). The influence of the accessory genome on bacterial pathogen evolution. Mob. Genet. Elements 1, 55–65. doi: 10.4161/mge.1.1.16432
Joardar, V., Lindeberg, M., Jackson, R. W., Selengut, J., Dodson, R., Brinkac, L. M., et al. (2005). Whole-genome sequence analysis of Pseudomonas syringae pv. phaseolicola 1448A reveals divergence among pathovars in genes involved in virulence and transposition. J. Bacteriol. 187, 6488–6498. doi: 10.1128/JB.187.18.6488-6498.2005
Johnson, M., Zaretskaya, I., Raytselis, Y., Merezhuk, Y., McGinnis, S., and Madden, T. L. (2008). NCBI BLAST: a better web interface. Nucleic Acids Res. 36, W5–W9. doi: 10.1093/nar/gkn201
Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N., and Sternberg, M. J. E. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10, 845–858. doi: 10.1038/nprot.2015.053
King, E. O., Ward, N. K., and Raney, D. E. (1954). Two simple media for the demonstration of pyocyanin and fluorescin. J. Lab. Clin. Med. 44, 301–307.
Li, X., Xie, Y., Liu, M., Tai, C., Sun, J., Deng, Z., et al. (2018). oriTfinder: a web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements. Nucleic Acids Res. 46, W229–W234. doi: 10.1093/nar/gky352
Li, Z., Hiasa, H., Kumar, U., and DiGate, R. J. (1997). The traE gene of plasmid RP4 encodes a homologue of Escherichia coli DNA topoisomerase III. J. Biol. Chem. 272, 19582–19587. doi: 10.1074/jbc.272.31.19582
Ma, Z., Smith, J. J., Zhao, Y., Jackson, R. W., Arnold, D. L., Murillo, J., et al. (2007). Phylogenetic analysis of the pPT23A plasmid family of Pseudomonas syringae. Appl. Environ. Microbiol. 73, 1287–1295. doi: 10.1128/AEM.01923-06
MacDonald, E. M. S., Powell, G. K., Regier, D. A., Glass, N. L., Roberto, F., Kosuge, T., et al. (1986). Secretion of zeatin, ribosylzeatin, and ribosyl-1”-methylzeatin by Pseudomonas savastanoi: plasmid-coded cytokinin biosynthesis. Plant Physiol. 82, 742–747. doi: 10.1104/pp.82.3.742
Mansfield, J., Genin, S., Magori, S., Citovsky, V., Sriariyanum, M., Ronald, P., et al. (2012). Top 10 plant pathogenic bacteria in molecular plant pathology. Mol. Plant Pathol. 13, 614–629. doi: 10.1111/j.1364-3703.2012.00804.x
Martínez, J. L. (2018). Ecology and evolution of chromosomal gene transfer between environmental microorganisms and pathogens. Microbiol. Spectr. 6, MTBP-0006-2016. doi: 10.1128/microbiolspec.MTBP-0006-2016
Matas, I. M., Lambertsen, L., Rodríguez-Moreno, L., and Ramos, C. (2012). Identification of novel virulence genes and metabolic pathways required for full fitness of Pseudomonas savastanoi pv. savastanoi in olive (Olea europaea) knots. New Phytol. 196, 1182–1196. doi: 10.1111/j.1469-8137.2012.04357.x
Meyer, J. M., and Abdallah, M. A. (1978). The fluorescent pigment of Pseudomonas fluorescens: biosynthesis, purification and physicochemical properties. J. Gen. Microbiol. 107, 319–328. doi: 10.1099/00221287-107-2-319
Miller, J. H. (1992). A Short Course in Bacterial Genetics - A Laboratory Manual and Handbook for Escherichia coli and Related Bacteria. Plainview, NY: Cold Spring Harbor Laboratory Press.
Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Salazar, G. A., Sonnhammer, E. L. L., et al. (2020). Pfam: The protein families database in 2021. Nucleic Acids Res. 49, D412–D419. doi: 10.1093/nar/gkaa913
Mitchell, A., Chang, H.-Y., Daugherty, L., Fraser, M., Hunter, S., Lopez, R., et al. (2015). The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Res. 43, D213–D221. doi: 10.1093/nar/gku1243
Morris, C. E., Lamichhane, J. R., Nikolić, I., Stanković, S., and Moury, B. (2019). The overlapping continuum of host range among strains in the Pseudomonas syringae complex. Phytopathol. Res. 1, 4. doi: 10.1186/s42483-018-0010-6
Neale, H. C., Slater, R. T., Mayne, L.-M., Manoharan, B., and Arnold, D. L. (2013). In planta induced changes in the native plasmid profile of Pseudomonas syringae pathovar phaseolicola strain 1302A. Plasmid 70, 420–424. doi: 10.1016/j.plasmid.2013.07.002
Nordström, K., and Austin, S. J. (1989). Mechanisms that contribute to the stable segregation of plasmids. Annu. Rev. Genet. 23, 37–69. doi: 10.1146/annurev.ge.23.120189.000345
Nowell, R. W., Green, S., Laue, B. E., and Sharp, P. M. (2014). The extent of genome flux and its role in the differentiation of bacterial lineages. Genome Biol. Evol. 6, 1514–1529. doi: 10.1093/gbe/evu123
O'Brien, H. E., Thakur, S., Gong, Y., Fung, P., Zhang, J., Yuan, L., et al. (2012). Extensive remodeling of the Pseudomonas syringae pv. avellanae type III secretome associated with two independent host shifts onto hazelnut. BMC Microbiol. 12, 141. doi: 10.1186/1471-2180-12-141
Peña, A., Matilla, I., Martín-Benito, J., Valpuesta, J. M., Carrascosa, J. L., de la Cruz, F., et al. (2012). The hexameric structure of a conjugative VirB4 protein ATPase provides new insights for a functional and phylogenetic relationship with DNA translocases. J. Biol. Chem. 287, 39925–39932. doi: 10.1074/jbc.M112.413849
Pérez-Martínez, I., Zhao, Y., Murillo, J., Sundin, G. W., and Ramos, C. (2008). Global genomic analysis of Pseudomonas savastanoi pv. savastanoi plasmids. J. Bacteriol. 190, 625–635. doi: 10.1128/JB.01067-07
Popa, O., and Dagan, T. (2011). Trends and barriers to lateral gene transfer in prokaryotes. Curr. Opin. Microbiol. 14, 615–623. doi: 10.1016/j.mib.2011.07.027
Prentki, P., and Krisch, H. M. (1984). In vitro insertional mutagenesis with a selectable DNA fragment. Gene 29, 303–313. doi: 10.1016/0378-1119(84)90059-3
Quandt, J., and Hynes, M. F. (1993). Versatile suicide vectors which allow direct selection for gene replacement in Gram-negative bacteria. Gene 127, 15–21. doi: 10.1016/0378-1119(93)90611-6
Rapisarda, C., and Fronzes, R. (2018). Secretion systems used by bacteria to subvert host functions. Curr. Issues Mol. Biol. 25, 1–42. doi: 10.21775/cimb.025.001
Sambrook, J., Fritsch, E. F., and Maniatis, T. (1989). Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory.
Schäfer, A., Tauch, A., Jager, W., Kalinowski, J., Thierbach, G., and Pühler, A. (1994). Small mobilizable multipurpose cloning vectors derived from the Escherichia coli plasmids pK18 and pK19: selection of defined deletions in the chromosome of Corynebacterium glutamicum. Gene 145, 69–73. doi: 10.1016/0378-1119(94)90324-7
Sesma, A., Sundin, G. W., and Murillo, J. (2000). Phylogeny of the replication regions of pPT23A-like plasmids from Pseudomonas syringae. Microbiology 146, 2375–2384. doi: 10.1099/00221287-146-10-2375
Smillie, C., Garcillán-Barcia, M. P., Francia, M. V., Rocha, E. P. C., and de la Cruz, F. (2010). Mobility of plasmids. Microbiol. Mol. Biol. Rev. 74, 434–452. doi: 10.1128/MMBR.00020-10
Solovyev, V., and Salamov, A. (2011). “Automatic annotation of microbial genomes and metagenomic sequences,” in Metagenomics and Its Applications in Agriculture, Biomedicine and Environmental Studies, ed R. W. Li (Hauppauge, NY: Nova Science Publishers), 61–78.
Souza, D. P., Oka, G. U., Alvarez-Martinez, C. E., Bisson-Filho, A. W., Dunger, G., Hobeika, L., et al. (2015). Bacterial killing via a type IV secretion system. Nat. Commun. 6, 6453. doi: 10.1038/ncomms7453
Stavrinides, J., and Guttman, D. (2004). Nucleotide sequence and evolution of the five-plasmid complement of the phytopathogen Pseudomonas syringae pv. maculicola ES4326. J. Bacteriol. 186, 5101–5115. doi: 10.1128/JB.186.15.5101-5115.2004
Stothard, P. (2000). The sequence manipulation suite: JavaScript programs for analyzing and formatting protein and DNA sequences. Biotechniques 28, 1102–1104. doi: 10.2144/00286ir01
Sundin, G. W., and Bender, C. L. (1993). Ecological and genetic analysis of copper and streptomycin resistance in Pseudomonas syringae pv. syringae. Appl. Environ. Microbiol. 59, 1018–1024. doi: 10.1128/aem.59.4.1018-1024.1993
Sundin, G. W., and Murillo, J. (2009). “Gene traders: characteristics of native plasmids from plant pathogenic bacteria,” in Plant Pathogenic Bacteria: Genomics and Molecular Biology, ed R. W. Jackson (Cambs: Caister Academic Press), 295–310.
Szabo, L. J., and Mills, D. (1984). Integration and excision of pMC7105 in Pseudomonas syringae pv. phaseolicola: involvement of repetitive sequences. J. Bacteriol. 157, 821–827. doi: 10.1128/jb.157.3.821-827.1984
Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M., et al. (2012). Primer3-new capabilities and interfaces. Nucleic Acids Res. 40, e115. doi: 10.1093/nar/gks596
Vivian, A., Murillo, J., and Jackson, R. W. (2001). The role of plasmids in phytopathogenic bacteria: mobile arsenals? Microbiology 147, 763–780. doi: 10.1099/00221287-147-4-763
Wallden, K., Rivera-Calzada, A., and Waksman, G. (2010). Microreview: type IV secretion systems: versatility and diversity in function. Cell. Microbiol. 12, 1203–1212. doi: 10.1111/j.1462-5822.2010.01499.x
Watanabe, K., Nagahama, K., and Sato, M. (1998). A conjugative plasmid carrying the efe gene for the ethylene-forming enzyme isolated from Pseudomonas syringae pv. glycinea. Phytopathology 88, 1205–1209. doi: 10.1094/PHYTO.1998.88.11.1205
Wawrzyniak, P., Płucienniczak, G., and Bartosik, D. (2017). The different faces of rolling-circle replication and its multifunctional initiator proteins. Front. Microbiol. 8, 2353. doi: 10.3389/fmicb.2017.02353
Zhao, M., Butler, M., Taiaroa, G., and Poulter, R. (2018). Conjugation in Pseudomonas syringae pv. actinidiae (Psa). Acta Hortic. 1218, 327–332. doi: 10.17660/ActaHortic.2018.1218.45
Zhao, Y. F., Ma, Z. H., and Sundin, G. W. (2005). Comparative genomic analysis of the pPT23A plasmid family of Pseudomonas syringae. J. Bacteriol. 187, 2113–2126. doi: 10.1128/JB.187.6.2113-2126.2005
Zrimec, J., and Lapanje, A. (2018). DNA structure at the plasmid origin-of-transfer indicates its potential transfer range. Sci. Rep. 8, 1820. doi: 10.1038/s41598-018-20157-y
Keywords: cytokinin biosynthesis, type III effectors, olive knot disease, type IV secretion system, mating pair formation, Pseudomonas amygdali, Pseudomonas savastanoi, conjugation (mating)
Citation: Añorga M, Urriza M, Ramos C and Murillo J (2022) Multiple relaxases contribute to the horizontal transfer of the virulence plasmids from the tumorigenic bacterium Pseudomonas syringae pv. savastanoi NCPPB 3335. Front. Microbiol. 13:1076710. doi: 10.3389/fmicb.2022.1076710
Received: 21 October 2022; Accepted: 21 November 2022;
Published: 12 December 2022.
Edited by:
Marco Scortichini, Council for Agricultural and Economics Research (CREA), ItalyReviewed by:
Marcus Dillon, University of Toronto Mississauga, CanadaYong-Qiang He, Guangxi University, China
Copyright © 2022 Añorga, Urriza, Ramos and Murillo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jesús Murillo, jesus.murillo@unavarra.es