Microbial strains, plasmids, primers, and media
Microbial strains and plasmids used in this study are listed in Table 1, and Extended Data Tables 1 and 4. E. amylovora strains, P. soli, P. savastanoi pv. savastanoi, P. syringae strains, R. solanacearum, Dickeya and Pectobacterium strains, C. michiganensis strains, X. campestris strains, and X. arboricola pv. juglandis were grown in Luria Bertani (LB) broth at 28°C. X. axonopodis pv. citri strains were grown in NA (nutrient broth) medium (beef extract, 3 g/L; yeast extract, 1 g/L; polypeptone, 5 g/L; and sucrose, 10 g/L) at 28°C. Fermentation of P. soli 0617-T307 was conducted in YM (yeast extract, 4 g/L and malt extract, 10 g/L) or YME medium (yeast extract, 4 g/L; glucose, 4 g/L; and malt extract, 10 g/L) at 16°C. Venturia inaequalis was grown on PDA (potato dextrose agar), while P. infestans was grown on RYE medium (dry rye berries, 60 g/L and sucrose, 20 g/L) at room temperature (22 ºC). Oligonucleotide primers used for cloning are listed in Supplementary Tables 1 and 2.
Genome sequencing, assembly, and annotation
High molecular weight genomic DNA of P. soli 0617-T307 was extracted, and the quality of the obtained DNA was checked by spectrophotometry at Next Generation Sequencing Core (UW-Madison, Madison, WI). Genome sequencing was conducted on the Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio) HiFi platforms. ONT facilitated the assembly of genomes with reads ranging from 50 to 120 kb in length. Consensus error corrections on the genomes and/or additional extrachromosomal elements were performed with PacBio reads at 8–14 kb size that were mapped against the assembly created from ONT reads. Genome sequencing resulted in a total of one single circular contig with a length in the 1 + MB range. For genome assembly and annotation, the polished contigs were compared against a BRC-curated subset of NCBI Prokaryotic RefSeq and GenBank accessions. This involved using a custom database of prokaryote sequences constructed by UW-Madison, sourced from NCBI on January 27, 2020. The five best matches are sorted by a BLAST + v2.8.0 bit score (blastn). A comparative analysis of the assembled contig and the highest scoring NCBI match is made using MUMmer4 56. Each contig × NCBI reference MUM comparison was filtered requiring an exact match length of at least 2 kb. The dotplot was generated with MUMmer4 by computing maximal exact matching, match clustering, and alignment extension between the contig and the single best-match NCBI sequence. The assembly features of the polished assembly were depicted by CIRCOS 57. Annotation of coding regions (genes on forward and reverse strands) including ORFs was determined by PROKKA 58. The boundary was defined by the gene dnaA, a protein that activates the initiation of DNA replication in nearly all bacteria (the genes dnaN and gyrB are usually associated with dnaA) 59. A low error-rate assembly was generated after three error-type corrections. The high-quality complete genome was deposited in Genbank with accession number CP151184 (BioProject accession: PRJNA1094439).
Species identification
For species identification, the genome sequences of representative Pseudomonas species were obtained from NCBI RefSeq database. The marker genes were parsed from the genome sequences and analyzed according to the guidelines established for Pseudomonas 32. The procedure for molecular phylogenetic analysis was based on that described previously 60. Briefly, the multiple sequence alignment was performed using MUSCLE v3.8.31 (10.1093/nar/gkh340) and the maximum likelihood phylogeny was inferred using PhyML v3.3.20180621 (10.1093/sysbio/syq010). The genome-wide average nucleotide identity was calculated using FastANI v1.1 33.
Construction of deletion and complementation strains
Deletion mutants Δras1 to Δras6 were generated using a double cross-over gene knock-out method as previously described 61. Sequences flanking ras1 at 714 bp upstream and 910 bp downstream were amplified by PCR using primers XbaI-ras1-UF/ras1UR and ras1-UF/EcoRI-ras1-DR (Supplementary Table 1), respectively. The two fragments were fused by overlapping PCR and cloned into a suicidal plasmid pEX18-Gm with restriction sites of XbaI and EcoRI. The construct was first transformed into E. coli S17-1 and conjugated with P. soli 0617-T307 on LB agar plate at 28°C. The cells on the plate were then rinsed off with 0.9% NaCl solution and spread on a selection plate with gentamicin (50 µg/ml) and carbenicillin (100 µg/ml) at 28°C for 2 days. P. soli 0617-T307 is naturally resistant to carbenicillin, while S17-1 does not. The selection pressure of gentamicin forced the integration of the plasmid into the genome through homologous recombination (first homologous recombination) at the ras1 upstream location. The positive clone was selected and incubated on YM medium with 12% sucrose for the second homologous recombination that forced the excision of the plasmid sequence. Due to the high G-C content in the P. soli 0617-T307 genome, the upstream and downstream primers were designed for efficient PCR to lower the annealing temperature by covering 15 bp upstream and 17 bp downstream sequences of ras1. Deletions of other ras genes were carried out in the same manner and the design for upstream and downstream fragments (bp) are listed in Supplementary Tables 1 and 2.
The Δras2Δras5 double mutant was made by deleting ras5 from the chromosome of mutant Δras2. The deletion construct pEX18-GmR-ras5 was delivered to E. coli S17-1 and transferred to mutant Δras2 by bi-parental mating. The deletion of targeted genes in all mutants was confirmed by PCR and DNA sequencing. Complementation of mutants was done by cloning the gene back to its original location through homologous recombination. The sequences upstream and downstream of the target gene used in the knock-out method, along with the target gene sequence, were amplified by PCR using XbaI-target gene-UF and EcoRI-target gene-DR (Supplementary Table 1). The amplified sequence was cloned into the suicide vector pEX18-GmR with the restriction sites XbaI and EcoRI. Subsequent steps of homologous recombination in deletion mutant strain were done as described in the knock-out method. The complementation of mutants was confirmed by PCR and DNA sequencing.
Greenhouse assays
Two-year-old cv. Snowdrift crabapple trees (Malus sp.) grafted onto cv. Dolgo rootstock and grown under greenhouse conditions (16-h light/8-h dark photoperiod at 28 ± 2°C and relative humidity of 60 ± 5%) were used in the experiments. At the 80% bloom stage, RAA solutions at 5 or 10 ppm with 0.12% (v/v) of surfactant (Regulaid; KALO, Overland Park, Kansas, U.S.) were sprayed directly onto the flowers in the evening. Flowers sprayed with water or 100 ppm of streptomycin (FireWall 50; AgroSource, Tequesta, FL) were used as negative or positive controls, respectively. A suspension of E. amylovora strain Ea110 (0.5 × 105 CFU/ml) was prepared from an overnight culture in LB broth at 28°C and sprayed onto the flowers the next morning, and was followed by a second application of RAA, water, or streptomycin in the evening. The incidence of fire blight was determined 5 days after inoculation. Flower clusters were counted as infected if more than one flower showed symptoms. Three trees were used per treatment. Statistical significance was determined using the least significant difference (LSD) method and one-way analysis of variance (ANOVA) analysis for comparisons between treatments (P < 0.05). The experiments were repeated three times independently.
Field trials
Over three years, field trials were conducted in California, Connecticut, Michigan, and New York utilizing local spray application tools on regional apple or pear cultivars (Supplementary Data 1). In 2022, Michigan field trials were conducted on 5-year-old ‘Buckeye Gala’ apple trees on M.9 rootstock at the Northwest Michigan Horticultural Research Center near Traverse City. Treatments were applied using 11.34-liter backpack sprayers (Model 473-P, Solo; Newport News, VA) on 16 May (70 to 80% bloom) and 18 May (full bloom). Inoculation with E. amylovora was done at 80% bloom on 17 May by spraying the outer perimeter of each tree with a backpack sprayer using an aqueous suspension of E. amylovora strain Ea110 (1.0 × 106 CFU/ml). Inoculation was conducted during the evening to ensure optimal conditions for bacterial survival. Blossom blight was assessed on 16 June. A total of 50 flower clusters from each of four replicate trees per treatment were evaluated for incidence of disease.
In 2023, a trial was conducted in New York at Cornell AgriTech in Geneva on 19-year-old ‘Idared’ apple trees on B.9 rootstock. Treatments were applied using a Solo 451 gas-powered mist blower (Solo Incorporated, Newport News, VA) calibrated to deliver 935.4 L ha− 1 (1.9 L/tree) at 80% bloom (5 May) and full bloom/early petal fall (9 May). Trees were inoculated at 80 to 90% bloom (8 May) with E. amylovora strain Ea273 (1 × 106 CFU/ml using a Solo 475-B backpack sprayer (Solo Incorporated). Disease was assessed on 2 June, and the incidence of fire blight was expressed as the number of blighted flowers out of five flowers in the cluster with 20 cluster assessments for six replicate trees per treatment for a total of 120 clusters per treatment.
In Connecticut, two trials were conducted at the Lockwood Farm of the Connecticut Agricultural Experiment Station in Hamden. In 2022, 40-year-old ‘Spartan’ apple trees were used. Treatments were sprayed using 18.92-liter backpack motorized sprayers (Solo 433, Newport News, VA) at 80–90% bloom on 6 May and at 100% bloom on 7 May, with approximately 1.9 liters per tree. Inoculation was conducted on 7 May by spraying E. amylovora strain MASHBO (1 × 106 CFU/ml) before the second application. Fire blight incidence was determined on 24 May by calculating the percentage of infected flower clusters of the total flower clusters. 100–300 flower clusters were evaluated on each tree. In 2023, experiments were conducted on 33-year-old ‘Early Macoun’ apple trees. Treatments were sprayed on 20 April (80–90% bloom) and 21 April (100% bloom), and inoculation was done on 21 April before the second spray. Disease was evaluated on 18 May.100–300 flower clusters were evaluated on each tree.
In California, the efficacy of RAA against fire blight was evaluated on approximately 25-year-old ‘Bartlett’ pear trees in Live Oak. Treatments were applied using a backpack air-blast sprayer (SR 430; Stihl Inc., Virginia Beach, VA, U.S.) at 935.4 L ha− 1 on 30 March (30% bloom) and 11 April (petal fall) 2023, and on 4 April (10% bloom) and 11 April (full bloom) 2024. Natural disease incidence was evaluated on 5 May 2023 or 18 April 2024, and the number of fire blight strikes on 100 flower clusters per replicate was counted.
Disease incidence data for all field trials were subjected to ANOVA for a randomized block design using Generalized Linear Mixed Model (GLIMMIX) procedures of SAS (version 9.4; SAS Institute Inc., Cary, NC). All percentage data were subjected to arcsine square root transformation prior to analysis. Multiple comparisons for significant fixed effects (P < 0.05) were determined using the LSMEANS procedure in SAS with an adjustment for Tukey’s HSD to control for family-wise error.
Fermentation and compound extraction
P. soli 0617-T307 was grown in 500 ml of YME media at 28°C for 24 h. Subsequently, the seed culture was inoculated into a 20-liter fermenter (BioFlo IV, New Brunswick Scientific Co., NJ) containing 12 liters of YME media. The fermentation proceeded at 16°C for 24 h. The agitation speed and the airflow rate were 200 rpm and 2 L/min, respectively.
Bacterial metabolites were extracted by ethyl acetate. The organic layer was separated and dried using sodium sulfate and rotary-evaporated at 35°C. Metabolites were then resuspended in 20 ml of methanol, and the methanol was evaporated in a fume hood. This resulted in 2.9 g of crude extract per 12-liter culture.
Compound separation and antimicrobial activity identification
The crude extract was dissolved in acetone and mixed with silica gel, which was loaded to a silica gel column (φ3.0 X 20 cm) on a flash chromatography system (Yamazen AI-580) equipped with a UV detector. The sample was eluted with 280 ml of each of the following solvents in order with increasing polarity: 100% hexane, 75% hexane/25% ethyl acetate, 50% hexane/50% ethyl acetate, 25% hexane/75% ethyl acetate, 100% ethyl acetate, 50% ethyl acetate/50% acetone, 100% acetone, and 100% methanol at a flow rate of 20 ml/min. The elute was monitored at UV 254 nm, fractions were collected by a time mode at 20 ml/tube, and 114 fractions/tubes (F1-113) were generated.
The fractions were used in antimicrobial plate assays. Aliquots of 1 ml of each fraction were first vacuum evaporated using a vacuum concentrator (Eppendorf, Enfield, CT) and re-dissolved in 50 µl DMSO. To test the antimicrobial activity, overnight cultures of E. amylovora were 1:100 diluted with water (~ 108 CFU/ml), spread onto LB agar plates, and 2 µl-droplets of each re-dissolved fraction were added equidistantly. DMSO alone was used as a negative control. The plates were then incubated at 28°C for 24 h, and the presence or absence of inhibition zones was observed.
The flash chromatographic fractions containing RAA (F38-40) and RAB (F50-54) were subjected to prep-HPLC purification on an Agilent C18 column (2.12 × 25 cm, 3.5 µm) with mobile phase A: Water with 0.1% formic acid, and mobile phase B: Methanol with 0.1% formic acid. The flow rate was 8.0 mL/min. The elute was monitored at 254 nm using a DAD detector. For RAA, the gradient program 40 to 100% B in 19 min was used. RAA was eluted at Rt 17.5 min. For RAB, the gradient program 20 to 60% B in 10 min was used. RAB was eluted at Rt 10.5 min.
RAA and RAB characterization
The structures of the two compounds were investigated by multiple analytical techniques, including high-resolution mass spectrometry (HR-MS), infrared spectroscopy (IR), ultraviolet spectroscopy (UV), 1D/2D nuclear magnetic resonance (NMR), and X-ray crystallography (Supplementary Data 2 and 3). RAA and RAB crystals were obtained through slow evaporation of their respective methanol solutions at room temperature. The X-ray diffraction analysis was performed at the Department of Chemistry, Marquette University, Milwaukee, WI, using an Oxford Diffraction SuperNova kappa-diffractometer equipped with dual microfocus Cu/Mo X-ray sources at 100K with Cu(Kα) radiation.
HPLC analytical methods
Analytical HPLC was done using an Agilent 1260 Infinity II system (Agilent, Santa Clara, CA). For the analysis of RAA (Method A), a PHENOMENEX 00B-4018-E0 3 µm, 50 x 4.6 mm column was used to achieve separation. Detection occurred at 406 nm with a retention time of 2.5 min. The mobile phase consisted of 10% acetonitrile (ACN) and 90% water + 0.1% formic acid. The flow rate was set to 0.6 ml/min, and the autosampler was configured to inject 10 µl aliquots of each sample. The standard curve of HPLC-purified RAA was used to determine RAA concentrations.
For HPLC analysis of RAB (Method B), a Phenomenex® Luna® Phenyl Hexyl HPLC Column 3 µm 150 X 4.6 mm 00f-4256-e0, was used. Detection occurred at 254 nm with a retention time of 9.8 min. The mobile phase consisted of 10% acetonitrile (ACN) and 90% water + 0.1% formic acid. The flow rate was set to 0.4 ml/min, and the autosampler was configured to inject 10-µl aliquots of each sample. The standard curve of HPLC-purified RAB was used to determine RAB concentration.
LC-MS analysis
LC-MS analysis was performed using a Shimadzu LCMS-2020 system (Shimadzu, Japan). Chromatographic separation was achieved on a Phenomenex Luna Phenyl Hexyl column (150 × 4.6 mm, 3 µm) maintained at 40°C. The mobile phases consisted of 0.1% formic acid in water (A) and acetonitrile (B). The flow rate was set at 0.4 mL/min with a gradient elution program: 10% B (0–1.0 min), 10–90% B (1.0–10.0 min), 90% B (10.0–14.9 min), 90 − 10% B (14.9–15.0 min), and 10% B (15.0–20.0 min). The injection volume was 5 µL, and the total run time was 20 min. Mass spectrometric detection was carried out using combined electrospray ionization (ESI) and atmospheric pressure chemical ionization (APCI) in positive mode. The desolvation line was heated to 250°C, and the heat block was at 400°C. Nebulizing gas flow was 1.3 L/min, and drying gas flow was 10.0 L/min. Data acquisition and analysis were performed using LabSolutions software (Shimadzu, Japan).
Antimicrobial Assay
To assess the growth inhibition of E. amylovora strains by P. soli, 100 µl of P. soli and 10 µl of E. amylovora at an optical density of OD5900 = 2.0 were co-cultivated in 10 ml of YM medium in 150-ml Erlenmeyer flasks at 28°C with shaking at 220 rpm. The population of E. amylovora was quantified after 18 h by colony forming unit assay.
The antimicrobial assay was conducted following the CLSI Antimicrobial Susceptibility Testing (AST) Standards, as outlined in CLSI document M07-A10 issued in January 2015. In brief, overnight bacterial cultures were diluted to an OD590 of 0.01 using the appropriate media and then aliquoted into 96-well plates with 200 µl per well. RAA or streptomycin was added to each well to final concentrations of 40, 20, 10, 5, 2.5, 1.25, 0.62, 0.31, 0.15, or 0.07 µg/ml, respectively. Control wells were supplemented with either water (for streptomycin) or DMSO (for RAA). The plates were then incubated at 28°C without agitation. MIC values, defined as the lowest concentration of the compound that resulted in no bacterial growth after 24 h, were determined.
V. inaequalis was cultivated on PDA agar in the dark at room temperature (22 ºC). A suspension containing conidia and mycelia in 0.01 M PBS was prepared, and 10-µl aliquots were applied to each RAA-amended PDA plate. Plates containing 0.4% DMSO served as a negative control, while those with 1000 µg/ml CuSO4 were used as a positive control. Plates were incubated at room temperature in the dark, and the diameter of each V. inaequalis colony was measured after 14 days. MIC values were determined after 7 days. P. infestans strains were cultured on RYE plates at room temperature. Mycelia were washed with 0.01 M PBS, centrifuged for 20 min at 3,500 rpm, and resuspended in sterile distilled water. Ten µl of this suspension was placed onto each of three equal-sized sections of each RAA-amended plate. RYE plates containing 0.4% DMSO served as the negative control. Plates were incubated at room temperature in the dark, diameters of P. infestans colonies were measured after 4 days, and MIC values were determined.