The Exometabolome of Xylella fastidiosa in Contact with Paraburkholderia phytofirmans Supernatant Reveals Changes in Nicotinamide, Amino Acids, Biotin, and Plant Hormones

Microbial competition within plant tissues affects invading pathogens’ fitness. Metabolomics is a great tool for studying their biochemical interactions by identifying accumulated metabolites. Xylella fastidiosa, a Gram-negative bacterium causing Pierce’s disease (PD) in grapevines, secretes various virulence factors including cell wall-degrading enzymes, adhesion proteins, and quorum-sensing molecules. These factors, along with outer membrane vesicles, contribute to its pathogenicity. Previous studies demonstrated that co-inoculating X. fastidiosa with the Paraburkholderia phytofirmans strain PsJN suppressed PD symptoms. Here, we further investigated the interaction between the phytopathogen and the endophyte by analyzing the exometabolome of wild-type X. fastidiosa and a diffusible signaling factor (DSF) mutant lacking quorum sensing, cultivated with 20% P. phytofirmans spent media. Liquid chromatography–mass spectrometry (LC-MS) and the Method for Metabolite Annotation and Gene Integration (MAGI) were used to detect and map metabolites to genomes, revealing a total of 121 metabolites, of which 25 were further investigated. These metabolites potentially relate to host adaptation, virulence, and pathogenicity. Notably, this study presents the first comprehensive profile of X. fastidiosa in the presence of a P. phytofirmans spent media. The results highlight that P. phytofirmans and the absence of functional quorum sensing affect the ratios of glutamine to glutamate (Gln:Glu) in X. fastidiosa. Additionally, two compounds with plant metabolism and growth properties, 2-aminoisobutyric acid and gibberellic acid, were downregulated when X. fastidiosa interacted with P. phytofirmans. These findings suggest that P. phytofirmans-mediated disease suppression involves modulation of the exometabolome of X. fastidiosa, impacting plant immunity.


Introduction
Deciphering the molecular aspects of the interaction between Xylella fastidiosa and its plant hosts can provide important clues about disease development [1].One key aspect leading to disease is its ability to modulate colonization behavior in plants and in insect vector transmission through quorum sensing (QS) mechanisms [2,3].For this, X. fastidiosa uses a family of modified fatty acids known as diffusible signaling factors (DSFs) that control individual and collective behavior and expression of virulence factors [4].Many aspects contribute to the complexity of the disease at different molecular levels.These include the bacterium's capacity to become systemic, its ability to evade detection by the host immune system, and the subsequent overreaction of the plant immune system, causing, among other effects, plant water deprivation [1,[5][6][7].
Microorganisms, including symbiotic bacteria, have the capacity to influence their host behavior by regulating the synthesis of specific compounds, thus fulfilling the metabolic and protein requirements of their host organisms [8][9][10].Nevertheless, niche competition among microorganisms-a naturally environmental occurrence-such as a spatial dispute between a pathogen and an endophyte often alters the fitness dynamics within the host [11,12].Nutrient competition, nutrient provision, toxin secretion, and competitor predation-which are some relational possibilities among microorganisms and their hosts-can be unveiled by metabolite-level alteration assessed by metabolomics [10,[13][14][15][16].
X. fastidiosa causes diseases in several plant hosts of economic importance to world agriculture, including grapevines, almonds, citrus, and more recently, olive trees [1,6,17,18].This bacterium exclusively colonizes the lumen of xylem vessels of its plant hosts and is transmitted by insect vectors like sharpshooters throughout the Americas and the spittlebug Philaenus spumarius in Europe [19][20][21].X. fastidiosa produces biofilm and secretes several virulence factors such as cell wall-degrading enzymes (CWDEs) and lipases/esterases, among others [22][23][24].X. fastidiosa does not have the type III secretion system (T3SS), which in most bacterial pathogens is responsible for the delivery of effectors/virulence factors inside host cells [25][26][27], thus reflecting the lifestyle of this phytopathogen that colonizes xylem vessels, a tissue that consists mostly of dead lignified cells [28].
Paraburkholderia phytofirmans strain PsJN [29], known as a grapevine endophyte, multiplies within grapevines, similarly to X. fastidiosa [30,31].Co-inoculation with X. fastidiosa results in a significant reduction in the phytopathogenic population and disease symptoms as previously reported for other pathogens [32].This effect involves priming innate disease resistance pathways in grapevines, leading to reduced symptoms when infected with X. fastidiosa [33].Recently, in a subsequent study, P. phytofirmans PsJN, when topically applied with a surfactant, effectively controlled Pierce's disease (PD) in grapevines, demonstrating systemic protection against Xylella fastidiosa infection, although the extent of the protection was spatially limited, with potential implications for this endophyte's use as a biological control agent [34].
Exometabolomics (or footprinting), i.e., the analysis of secreted metabolites, is a very effective approach to tracking signatures of microorganisms, e.g., contamination signatures [35].Integrating metabolomics and genomics helps in tracking the origins of major trends in metabolite levels in different species, but also while interacting with each other [36,37].In this sense, the characterization of the exometabolome is a useful tool for further understanding the responses of X. fastidiosa to its environment and/or its interactions with other organisms.In addition, the analysis of exometabolomes has many analytical advantages over the analysis of intracellular compounds due to their lower turnover, higher stability, and, consequently, greater reproducibility of their metabolic footprint [38,39].
Here, we aimed to investigate the indirect interaction between a phytopathogen and endophyte by analyzing the exometabolome of X. fastidiosa, WT and ∆rpfF, both cultivated in a 20% P. phytofirmans spent medium.First, we showed that the P. phytofirmans spent medium was capable of disturbing X. fastidiosa mature biofilm formation in vitro.Additionally, we performed LC-MS analysis, revealing 121 metabolites in total, of which 25 were further investigated.Finally, using integration tools including MAGI, OmicsBox, and NetSeed, we mapped metabolites to genomes, assessed the metabolic networks, and annotated the functions of the mapped genes, respectively.This study provides further insights into the effects of QS and disease development.
For culture in liquid medium, after 7 days of culture in PD3 agar, bacteria were transferred to 100 mL of PD3 medium and maintained at 28 • C at 100 rpm for up to 7 days.Cultures were started with OD 600nm = 0.05.
P. phytofirmans PsJN was selected from a plate of King's B medium [44] containing rifampicin (KBR) at 100 µg/mL.Next, P. phytofirmans was transferred to 100 mL of PD3.X. fastidiosa and P. phytofirmans were grown at 28 • C and 100 rpm.Fresh culture was started from frozen stocks at −80 • C.
Alternatively, WT and ∆rpfF strains were grown in PD3 supplemented with 20% of a P. phytofirmans spent medium.P. phytofirmans PsJN grown for 1 day in PD3 medium was centrifugated at 4000× g for 30 min, 4 • C. The supernatant was collected, adjusted to pH 4, and subjected to vacuum filtration through a 0.22 µm membrane.The eluate was added to a regular PD3 medium to a final concentration of 20% and used for subsequent experiments.

X. fastidiosa Biofilm Measurement
Aliquots (2 µL, 20 µL, 200 µL, and 400 µL) of the P. phytofirmans PsJN spent medium eluate were added to a total 1.0 mL of PD3 liquid culture and then used to grow WT X. fastidiosa for 7 days at 28 • C and 200 rpm in a 24-well plate.After this period, the biofilm formed was quantified by staining with the crystal violet method, loosely, as previously described [45].Briefly, the biofilm deposited at the air-medium interface of the WT strain was washed with distilled water, stained for 20 min with 0.1% crystal violet reagent, washed again with distilled water, and then taken up in acetone:ethanol (6:4).Quantification was performed by measuring the absorbance at 595 nm (Abs 595nm ).WT X. fastidiosa culture without spent medium addition grown under the same conditions was used as the control.The experiment was repeated and this time an aliquot of 500 µL of the P. phytofirmans PsJN spent medium eluate was added to 2.5 mL of PD3 liquid culture in glass tubes and then used to grow either WT or ∆rpfF X. fastidiosa strains as described before.Cultures of WT or ∆rpfF X. fastidiosa strains without spent medium addition grown under the same conditions were used as the control.

Analysis of Metabolites of the Supernatant of X. fastidiosa Cultures by Liquid Chromatography Coupled to Mass Spectrometry (LC-MS)
Samples (1 mL) of spent medium from cultures and sterile control medium were centrifuged in 1.5 mL microcentrifuge tubes at 1700× g for 5 min.Supernatants, containing extracellular metabolites, were lyophilized and then resuspended in LC-MS-grade methanol (300 µL).Resuspended samples were centrifuged again to pellet salts.Supernatants were dried under vacuum (Savant SpeedVac Plus SC110A, Holbrook, NY, USA) and resuspended in 500 µL of LC-MS-grade methanol containing a mixture of internal standards (25 µM of 13C-15N-L-phenylalanine, 2-amino-3-bromo-5-methylbenzoic acid, 3,6-dihydroxy-4-methylpyridazine, d4 lysine, d5-benzoic acid, and 9-anthracene carboxylic acid).The samples were filtered through a 0.22 µm microcentrifuge filtration device and transferred to 1.5 mL borosilicate glass vials (Agilent, Santa Clara, CA, USA) for LC-MS analysis as previously described [46].Extraction blanks (blank microcentrifuge tubes taken through the entire extraction procedure to check for contaminants introduced during sample preparation) were included in the analysis.A quality control mixture (a defined mixture of common metabolites) and injection blanks (pure methanol) were also analyzed at the beginning and ends of each run to ensure no drift in retention times or signal abundances and no signs of column fouling or metabolite carryover.For polar metabolite analysis, an Agilent 1290 LC system equipped with a ZIC-pHILIC column (150 mm × 1 mm, 3.5 µm 100 Å, Merck SeQuant, Darmstadt, Germany) was used for metabolite separation with the following LC conditions: solvent A, 5 mM ammonium acetate; solvent B, 9:1 acetonitrile:H 2 O with 5 mM ammonium acetate; timetable, 0 min at 100% B, 1.5 min at 100% B, 21 min at 0% B, 27 min at 0% B, 33 min at 100% B, and 45 min at 100% B; 0.8 mL/min; column compartment temperature 40 • C. Mass spectrometry analyses were performed using a 6550 quadrupole time-of-flight mass spectrometer.Mass spectrometry data were collected at a mass range of 30-1200 m/z, drying gas rate of 11 L/min, and gas temperature of 290 • C. The nebulizer pressure was set at 30 psi and capillary voltage, at 3500 V. Data were acquired in both positive (+) and negative (−) polarities.The extracted samples were obtained in triplicate for both biological and technical replicates.Additionally, three replicate samples from the sterile control media and eleven samples from extraction blanks were processed under the same conditions.A total of 62 injections were performed, including 3-11 replicate injections.

LCMS Data Processing and Analysis
MassHunter qualitative analysis software (Agilent, Santa Clara, CA, USA) was used to inspect the raw data.Then, the raw data obtained were parsed using Python scripts within MetAtlas [47], which was used to extract out putative metabolite peaks using an in-house metabolite database containing m/z and retention time data [47,48].Next, the LC-MS data were analyzed using a custom code in Perl.Each of the detected peaks was evaluated to assign a level of confidence in the identification of the compound.Compounds that were identified had a corresponding retention time and m/z for a pure standard using the same methods as above.Negative and positive polarity mode lists of compounds were retrieved and merged.A quality control analysis of retention time and intensity stability was performed using internal and external standards; based on the quality control assessment, the negative mode was selected for peak integration due to better signal and retention time stability over positive mode.In order to confirm the identity of a compound, we checked the RT difference (i.e., abs [RT standard − RT µexperiment ]) and the ppm error mass (i.e., [( abs [m/z standard − m/z µexperiment ])/m/z standard ] × 10 6 ) (Supplementary Table S1).Compounds with an RT difference ≤ 0.05 and ppm error mass ≤ 6 were considered to have had their identities confirmed (Supplementary Table S1) [49].Compounds that did not pass one or both criteria (up to a ppm error mass ≤ 15) were considered putative.Next, compound abundance was normalized by strains OD in each biological replicate (Supplementary Table S2).
For further comparisons among the conditions, the compounds that satisfied the following criteria were considered expressed in the supernatants: [µi sample − µi blank ] − [µi medium − µi blank ] > 0, where µ stands for "average of peak area value" and i stands for any detected compound.Finally, an ANOVA test was conducted to acknowledge variability across all conditions (Supplementary Table S2).Compounds with a p-value ≤ 0.05 were considered statistically significant.Pairwise comparison among conditions was performed and a t-test with the Benjamini-Hochberg procedure at a false discovery rate (FDR) ≤ 0.05 was calculated for each compound.Compounds with statistical significance between two conditions were considered differentially expressed metabolites (DEMs) (Supplementary Table S3).Extracted ion chromatograms (EICs) were generated, one representative EIC of each extract or exometabolome was used for building a summary panel for each discussed compound, and all the remaining EICs from all biological and technical replicates were made available in the Supplementary Material (Supplementary Figures S1-S3).
A Python library for the Venn diagram was used to highlight unique and shared compounds across the conditions tested.These were referred to as Xf (X.fastidiosa WT Temecula1), ∆rpfF (X.fastidiosa ∆rpfF), Xf sm (X.fastidiosa Temecula1 WT cultivated in PD3 medium supplemented with 20% P. phytofirmans spent medium), ∆rpfF sm (X.fastidiosa ∆rpfF cultivated in PD3 medium supplemented with 20% P. phytofirmans spent medium), and P. phytofirmans PsJN (Pp).A matrix was built from three biological replicates of Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp.This matrix was given as an input to a multivariate data analysis in a principal coordinate analysis (PCoA) [50], taking as a parameter the Euclidean distance between the samples.To compute this dataset, we used the Python library skbio and the library emperor to build the PCoA plot.Permanova and permdisp statistical analysis were made with the PCoA input data using skbio modules permanova and permdisp, respectively, and the information was added later to the PCoA plot.
A heatmap based on clustering of the compounds across conditions was generated with the Python libraries matplotlib and seaborn.The expression values of all compounds under all conditions were normalized by a z-score.A color code and symbols were added to compound names to indicate their identification level and statistical significance.

Integration of Metabolites and Gene Set for Functional Annotation of Exometabolomes
Whole genome annotation (in FASTA format) from Xylella fastidiosa Temecula1 and Paraburkholderia phytofirmans PsJN together with a list of InChiKey identifiers from all the compounds assigned as expressed in previous steps were given as the input for a container (Docker) of MAGI v1.0 [51,52] for metabolite and gene integration.The output spreadsheets were filtered by metabolites expressed under each given condition.Next, all genes annotated by MAGI were retrieved and duplicates were filtered out.Each gene list representing an exometabolome was submitted to the software OmicsBox 3.1 [53] for functional annotation using gene ontology.Enrichment analysis was performed, annotations with an adjusted p-value ≤ 0.05 were collected, and the values were log 2 transformed and ranked.The top 5 under-and overenriched annotations from each condition were collected for comparison.The matplotlib and seaborn Python libraries were used to build a bar plot.

Bioinformatic Analysis
Expressed metabolites were compared in box plot and bar plot, coded by the matplotlib and seaborn Python libraries.Pairwise significance comparison was performed by the Python library statannotations and the t-test argument was used.When necessary, the precomputed Benjamini-Hochberg procedure with an FDR < 0.5% for pairwise comparison was used for annotations later (Supplementary Table S2).
Metabolic reactions of nicotinamide, nicotinate, and biotin were freely redrawn from annotated X. fastidiosa and P. phytofirmans Kegg pathways.Illustrations of gene operons were rebuilt based on the synteny presented in the browser IMG/DoE [54].Mummer [55] was used for the global alignment of genomic sequences.
For the metabolome and transcriptome integration, a bubble plot was created using the matplotlib and seaborn Python libraries.The metabolites' expression values were 2 z-scoretransformed in order to maximize expression differences and keep all values positive.

Cooperative Metabolic Interactions between X. fastidiosa and P. phytofirmans
In order to analyze the potential ecological microbial interactions [56] between X. fastidiosa and P. phytofirmans, a reverse ecology analysis was applied to predict the ecological structure of this symbiotic relationship [57].
For defining microbe-microbe cooperative and competitive potential in a pairwise manner, a local and customized NetSeed version [58] was used (to include the metabolic complementarity index and effective metabolic overlap index) [59,60].The metabolic complementarity index between two species measures the ratio (range 0-1) of acquired compounds in one species that are found in the metabolic network of the other species and are not used by it.The effective metabolic overlap index of two species is the fraction of compounds required by both species and is, therefore, a measure of potential nutritional competition [59], with values ranging from 0 (no competition) to 1 (high competition).

P. phytofirmans Interacts with X. fastidiosa through Its Exometabolome
In a prior study, Baccari and colleagues [33] demonstrated that direct contact between X. fastidiosa and P. phytofirmans was not necessary for reducing X. fastidiosa titer loads in xylem vessels.Here, we explore the indirect interaction between X. fastidiosa and P. phytofirmans, including the QS-insensitive mutant ∆rpfF.In our initial trials, we challenged X. fastidiosa WT Temecula1 with 0.2%, 2%, 20%, and 40% of P. phytofirmans spent media grown in 24-well plates (Figure 1A).Consequently, we established a specific condition for subsequent experiments.We used PD3 culture media supplemented with 20% P. phytofirmans PD3 spent media, which noticeably affected the X. fastidiosa phenotype.This impact was particularly evident in terms of biofilm intensity and deposition, as measured by the crystal violet assay (Figure 1).
were rebuilt based on the synteny presented in the browser IMG/DoE [54].Mummer [55] was used for the global alignment of genomic sequences.
For the metabolome and transcriptome integration, a bubble plot was created using the matplotlib and seaborn Python libraries.The metabolites' expression values were 2 z- score -transformed in order to maximize expression differences and keep all values positive.

Cooperative Metabolic Interactions between X. fastidiosa and P. phytofirmans
In order to analyze the potential ecological microbial interactions [56] between X. fastidiosa and P. phytofirmans, a reverse ecology analysis was applied to predict the ecological structure of this symbiotic relationship [57].
For defining microbe-microbe cooperative and competitive potential in a pairwise manner, a local and customized NetSeed version [58] was used (to include the metabolic complementarity index and effective metabolic overlap index) [59,60].The metabolic complementarity index between two species measures the ratio (range 0-1) of acquired compounds in one species that are found in the metabolic network of the other species and are not used by it.The effective metabolic overlap index of two species is the fraction of compounds required by both species and is, therefore, a measure of potential nutritional competition [59], with values ranging from 0 (no competition) to 1 (high competition).

P. phytofirmans Interacts with X. fastidiosa through Its Exometabolome
In a prior study, Baccari and colleagues [33] demonstrated that direct contact between X. fastidiosa and P. phytofirmans was not necessary for reducing X. fastidiosa titer loads in xylem vessels.Here, we explore the indirect interaction between X. fastidiosa and P. phytofirmans, including the QS-insensitive mutant ΔrpfF.In our initial trials, we challenged X. fastidiosa WT Temecula1 with 0.2%, 2%, 20%, and 40% of P. phytofirmans spent media grown in 24-well plates (Figure 1A).Consequently, we established a specific condition for subsequent experiments.We used PD3 culture media supplemented with 20% P. phytofirmans PD3 spent media, which noticeably affected the X. fastidiosa phenotype.This impact was particularly evident in terms of biofilm intensity and deposition, as measured by the crystal violet assay (Figure 1).WT strain was grown in a PD3 medium for 7 days with 0.2%, 2%, 20%, and 40% v/v of a 1-day PD3 P. phytofirmans (Pp) spent medium in a 24-well plate.WT strain was grown under same conditions without P. phytofirmans (Pp) spent medium and was set as a reference.(B) Temecula1 WT or ΔrpfF strains were grown in PD3 medium for 7 days with (Xf sm and ΔrpfF sm ) or without (Xf and ΔrpfF) 20% v/v of a 1-day PD3 P. phytofirmans (Pp) spent medium in a glass tube.Biofilm quantification was WT strain was grown in a PD3 medium for 7 days with 0.2%, 2%, 20%, and 40% v/v of a 1-day PD3 P. phytofirmans (Pp) spent medium in a 24-well plate.WT strain was grown under same conditions without P. phytofirmans (Pp) spent medium and was set as a reference.(B) Temecula1 WT or ∆rpfF strains were grown in PD3 medium for 7 days with (Xf sm and ∆rpfF sm ) or without (Xf and ∆rpfF) 20% v/v of a 1-day PD3 P. phytofirmans (Pp) spent medium in a glass tube.During further experiments, we cultivated both Xf and ∆rpfF in regular PD3 culture media or supplemented with 20% P. phytofirmans PsJN spent media (Xf sm or ∆rpfF sm ).The biofilm formation of Xf sm and ∆rpfF sm increased by 39% and 128%, respectively, in cultivations in glass tube (Figure 1B).Moreover, a larger amount of bacterial biomass was visually observed at the medium/air interface of Xf sm and ∆rpfF sm compared to Xf and ∆rpfF (Figure S1).As PD3 supplemented with 20% P. phytofirmans spent media led to quantifiable changes in the biofilm formation profile of both bacterial strains in in vitro cultures, we hypothesized that the secreted metabolomes (exometabolomes) of P. phytofirmans trigger a disturbance in the exometabolomes of the tested Xylella strains.

Exometabolome Variation among X. fastidiosa Strains in Response to P. phytofirmans
Subsequently, we investigated the exometabolomes of Xf and Pp to characterize the qualitative and quantitative signatures of their compounds.By employing an internal reference library, we were able to identify a total of 131 compounds across all strains and experimental conditions (Supplementary Tables S1-S3).Referring to the detailed criteria for detection outlined in the Section 2, we established a threshold and inferred the detection of 121 compounds among the entire spectrum of strains and conditions examined (Supplementary Table S3) [49].Specifically, within Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp, we detected 114, 119, 118, 115, and 48 compounds, respectively (Figure 2A).
There is a shared core of 43 compounds present in all five exometabolomes.Focusing exclusively on the X. fastidiosa strains and treatments (Xf, Xf sm , ∆rpfF, and ∆rpfF sm ), the Venn diagram highlights a core of 112 compounds.Notably, the Venn diagram also unveils a distinct subset of exclusive compounds (69) unique to X. fastidiosa (Xf, Xf sm , ∆rpfF, and ∆rpfF sm ).The preeminent class of compounds within the core exometabolome belongs to "nucleotides and derivatives" (comprising 32.6% of the total).Among the exclusively expressed compounds found in X. fastidiosa strains, amino acids and derivatives dominate (52.2%).It is noteworthy that nicotinamide stands as the sole compound exclusively identified in the Pp exometabolome, explored in further detail in the subsequent section.
Utilizing principal coordinate analysis (PCoA), we discerned a distinct clustering pattern within the exometabolomes, initially delineating separation based on bacterial species (Xf versus Pp), and subsequently segregating according to strains and treatments (Xf versus ∆rpfF, and Xf sm versus ∆rpfF sm ) (Figure 2B).This analysis yielded the formation of three primary and discernible clusters: (1) Pp, (2) Xf sm , ∆rpfF, and ∆rpfF sm , and (3) Xf.Notably, the implementation of permanova and permdisp analyses yielded statistically significant outcomes, denoted by p-values of <0.05 and <0.005, respectively, with respect to the distribution of samples.
Drawing from the outcomes of the PCoA, we postulate that the metabolites originating from Pp exert an influence on Xf exometabolomes, leading to a convergence in profile akin to that of the ∆rpfF metabolome.Interestingly, while ∆rpfF remains unaffected, as evidenced by the similarity to ∆rpfF sm samples, this transformative effect is conspicuously absent.The PCoA not only highlights sample correlation but also offers a platform for constructing a working hypothesis.In essence, the assimilation of Pp spent media by Xf potentially results in a convergence of exometabolomic profiles, aligning Xf more closely with the characteristics of ∆rpfF, at least at the exometabolome level.A noteworthy observation is the similarity apparent in the sample distribution between ∆rpfF and ∆rpfF sm .The distribution of exometabolome samples, as represented in Figure 2B, guides our hypothesis that the interaction of Xf with Pp spent media draws it near to the phenotype of ∆rpfF.The color code highlights the abundance of the metabolites under the conditions tested, ranging from low (blue) and medium (white) to high (red) expression.An additional column was added with the respective compound class as stated in PubChem.Two main clusters are highlighted, in green and red boxes, as they contain the compounds highly expressed in X. fastidiosa and P. phytofirmans, respectively.A z-score normalization was performed for each compound across the different conditions.Three biological replicates of each species and treatment were used for creating this heatmap.Compounds in black were confirmed by standard run comparison (regarding RT and m/z ppm), compounds in orange had a lower level of confirmation (either or both RT and m/z ppm), and compounds in red could not be confirmed by comparison RT.Low-confidence compounds were therefore considered as putative (possibly tautomers of the compounds assigned).(*) indicates ANOVA statistical significance across conditions.
Finally, we generated a heatmap illustrating the expression levels of compounds, normalized by 2 z-score .To enhance our comprehension of the functional attributes, variability, and expression levels of the metabolites, we organized the identified compounds into groups based on PubChem classes [61].These groups correspond to various chem-ical classes, including carbohydrates, amino acids (and derivatives), nucleotides (and derivatives), carboxylic acids, pyridines, amines, sulfurs, terpenes, polyols, parabens, and vitamins.Noteworthy are the compounds that serve as intermediates within metabolic pathways and the citric acid cycle (such as α-ketoglutarate, fumarate, lactate, and succinate), along with secondary metabolites (including suberic acid and shikimic acid), present in the exometabolomes of both X. fastidiosa and P. phytofirmans (Figure 2C and Supplementary Table S3).
Within Xf, 75 metabolites demonstrated significantly upregulated expression levels compared to Pp, ∆rpfF, Xf sm , and ∆rpfF sm .Of these, 44% were amino acids and derivatives.Similarly, Xf sm exhibited a majority of upregulated metabolites (5), accounting for 60% of amino acids and derivatives.Conversely, ∆rpfF and ∆rpfF sm displayed 8 and 12 upregulated compounds, respectively, with 50% and 41.7% being carboxylic acids.In contrast, the Pp exometabolome exhibited notably low metabolite detection, particularly when contrasted with those identified in Xf exometabolomes.Of the detected metabolites in Pp, 10 exhibited upregulation when compared to their counterparts in Xf exometabolomes, of which 50% were nucleotides or nucleotide derivatives.
The heatmap construction encompassed all compounds detected above the threshold of the blank (MeOH) and culture media (PD3 or PD3 sm ), totaling 121 compounds.Our annotation approach included two additional layers of information-statistical significance, reflecting variance among conditions, and the level of identification.The heatmap (Figure 2C) encompasses compounds that met the criteria of standard confirmation through RT and m/z ppm, and annotation also included assigning compounds demonstrating statistical significance (p < 0.05) based on ANOVA analysis (Supplementary Tables S2 and S3).A notable feature of the heatmap is the emergence of two distinct clusters: one primarily composed of compounds highly expressed in Pp and the other predominantly featuring compounds highly expressed in Xf strains.There exists a subset of compounds that are more highly expressed in Xf sm , ∆rpfF, and ∆rpfF sm in contrast to Xf. Upon a comprehensive review of the overall detected metabolites, we focused the investigation on a subset of 25 compounds (amino acids, vitamins, and plant hormones), as detailed in the subsequent sections.

X. fastidiosa Secretes High Amounts of Amino Acids and Vitamins
First, we chose to further investigate compounds that are well known for their role in the interaction between bacteria and some species of insect hosts [62].As such, we examined amino acid levels among the strains and conditions.
The nutritional statuses of amino acids within the exometabolomes of Xf, ∆rpfF, Xf sm , ∆rpfF sm , and Pp were assessed, primarily utilizing the ratio of glutamine (Gln) to glutamate (Glu).This ratio, commonly employed to gauge nitrogen (N) status in various eukaryotic cells, has been linked to N limitation when <0.2 and is indicative of N-replete cells when >0.5 [10,62].The computed median Gln:Glu ratios within the exometabolomes unveiled variations: Xf, ∆rpfF, Xf sm , and ∆rpfF sm exometabolomes all exhibited ratios predicting N-replete conditions, spanning from 0.52:1 in Xf to 8.33:1 in ∆rpfF (Figure 3A).Notably, Gln and Glu were not detected within the exometabolome of Pp under the assessed conditions.The 0.5 threshold for N-replete conditions was drawn in the plot, with all conditions (Xf, Xf sm , ∆rpfF, and ∆rpfF sm ) surpassing this threshold-except for Pp, which lies below, and Xf, which teeters at the threshold.Importantly, our analysis did not identify any statistically significant differences between Xf, Xf sm , ∆rpfF, and ∆rpfF sm .In order to provide a comprehensive understanding of compound levels within the experimental framework, we constructed a panel of EICs encompassing various controls, including the blank (MeOH) and culture media controls, PD3 and PD3 sm (Figure 3B,C; Supplementary Figure S1A-T).These control conditions were juxtaposed with the investigated states: Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp.As anticipated, the control conditions (culture media PD3 and PD3 sm ) exhibited detectable amino acid levels.Upon scrutinizing Xf, Xf sm , ∆rpfF, and ∆rpfF sm , a noticeable pattern emerged as the following: amino acid levels, whether pertaining to nEAAs or EAAs, consistently surpass those found in blank samples and even exceed the amino acid concentrations in culture media controls.This observation aligns with the knowledge that PD3, given its protein-rich composition, inherently contains a source of amino acids.Upon closer examination, this elevated amino acid content within Xf, Xf sm , ∆rpfF, and ∆rpfF sm challenges conventional assumptions, indicating a plausible role as an insect vector decoy within its presumed biological function.In stark contrast, the levels of EAAs and nEAAs in Pp are notably diminished, bordering on complete absence.
While the visual inspection of differences in the levels of essential amino acids (EAAs) and non-essential amino acids (nEAAs) does not consistently yield pronounced distinctions, we undertook individual assessments of these two amino acid categories.The nature of this evaluation displayed variability across species, strains, and treatments, with notable statistical significance emerging primarily between amino acid quantifications in Xf and Xf sm .Specifically, variations were evident in tyrosine, serine, proline, glutamine, aspartate, asparagine, alanine, and glycine within the nEAAs and valine, methionine, lysine, leucine, isoleucine, and arginine within the EAAs.Similarly, only histidine exhibited a discernible difference between ∆rpfF and ∆rpfF sm (Figure 3B,C).Notably, within the exometabolome of Pp, neither EAAs nor nEAAs were detected.
Given the disparities in ionization among amino acids, merging them into a unified value for comparison across diverse conditions was not viable.Nonetheless, a consistent observation emerged as the following: essential amino acids (EAAs) exhibit higher individual abundance within Xf, Xf sm , ∆rpfF, and ∆rpfF sm in comparison to non-essential amino acids (nEAAs).
We also investigated levels of complex B vitamins detected in the tested exometabolomes.These compounds are relevant to the relationship between an insect and host plant.Our findings illustrate the accumulation of nicotinamide in Pp, while it was notably absent in Xf, Xf sm , ∆rpfF, and ∆rpfF sm -this absence was evident across both normalized values (peak area intensity by OD) and EIC profiles, including in the comparison with the controls.Conversely, the presence of nicotinic acid follows an inverse pattern (Figure 4A-F, Supplementary Figure S2A-C).The product of pncA, namely nicotinic acid, exhibits high expression in Xf, along with downregulation and relatively consistent levels in Xf sm , ∆rpfF, and ∆rpfF sm .Biotin, on the other hand, was not detected in Xf, ∆rpfF, or Pp.In contrast, its presence was identified at relatively low levels in Xf sm and ∆rpfF sm .The pathways remain conserved in both X. fastidiosa and P. phytofirmans (Figure 4G,H).Main reactions in "nicotinate and nicotinamide metabolism" and in "biotin metabolism" according to KEGG pathways (www.kegg.jp/pathway/)(accessed on 1 July 2023) in (G) X. fastidiosa and (H) P. phytofirmans, respectively.The letter X in red indicates the absence of the enzyme that converts nicotinamide into nicotinic acid.

Two Plant Hormones Are Exclusively Secreted by X. fastidiosa
We conducted a comparative analysis of the levels of two compounds, 2-aminoisobutyric acid (AIB) and gibberellic acid (GA), detected within the exometabolomes of Xf, ΔrpfF, Xf sm , ΔrpfF sm , and Pp.These compounds possess implications for the interaction between X. fastidiosa and its plant hosts, potentially functioning as pathogenic metabolite effectors.
Across the conditions assessed in our metabolomic analysis, AIB and GA exhibit notably similar profiles.Both compounds form three distinct clusters: (1) Xf, (2) Xf sm , ΔrpfF, and EICs of representative assessments in the blank (MeOH), PD3, PD3 sm , Xf, ∆rpfF, Xf sm , ∆rpfF sm , and Pp extracts or exometabolomes.Main reactions in "nicotinate and nicotinamide metabolism" and in "biotin metabolism" according to KEGG pathways (www.kegg.jp/pathway/)(accessed on 1 July 2023) in (G) X. fastidiosa and (H) P. phytofirmans, respectively.The letter X in red indicates the absence of the enzyme that converts nicotinamide into nicotinic acid.

Two Plant Hormones Are Exclusively Secreted by X. fastidiosa
We conducted a comparative analysis of the levels of two compounds, 2-aminoisobutyric acid (AIB) and gibberellic acid (GA), detected within the exometabolomes of Xf, ∆rpfF, Xf sm , ∆rpfF sm , and Pp.These compounds possess implications for the interaction between X. fastidiosa and its plant hosts, potentially functioning as pathogenic metabolite effectors.
Across the conditions assessed in our metabolomic analysis, AIB and GA exhibit notably similar profiles.Both compounds form three distinct clusters: (1) Xf, (2) Xf sm , ∆rpfF, and ∆rpfF sm , and (3) Pp (Figure 5A,B, Supplementary Figure S3A,B).Comparatively, AIB and GA levels proved significantly higher in Xf than in any other strain or treatment.Xf sm displayed reduced AIB and GA expression, akin to the observed pattern in ∆rpfF.Conversely, ∆rpfF sm exhibited no significant alteration in AIB and GA expression compared to ∆rpfF.Notably, AIB remained undetected in the Pp exometabolome, while GA levels ranged from 6.63% to 26.89% of those detected in Xf, ∆rpfF, Xf sm , and ∆rpfF sm exometabolomes.
ΔrpfF sm , and (3) Pp (Figure 5A,B, Supplementary Figure S3A,B).Comparatively, AIB and GA levels proved significantly higher in Xf than in any other strain or treatment.Xf sm displayed reduced AIB and GA expression, akin to the observed pattern in ΔrpfF.Conversely, ΔrpfF sm exhibited no significant alteration in AIB and GA expression compared to ΔrpfF.Notably, AIB remained undetected in the Pp exometabolome, while GA levels ranged from 6.63% to 26.89% of those detected in Xf, ΔrpfF, Xf sm , and ΔrpfF sm exometabolomes.The chromatograms illustrate the contrast between the AIB levels in Pp, which appear lower than those in the controls, while the AIB levels in Xf consistently remain elevated (Figure 5C).GA was generally detected at low levels across all conditions, including Xf.However, both Xf and ΔrpfF exhibited levels 2 to 3 times higher than the controls.Xf sm and ΔrpfF sm exhibited increases of approximately 100% and 50%, respectively, compared to the culture medium controls.Notably, Pp's GA levels mirror those of the controls (Figure 5C,D).
We investigated the potential synthesis of AIB and GA by X. fastidiosa or P. phytofirmans through a search for gene clusters for their synthesis.We examined the presence of homologous genes or clusters akin to the reported aib operon from Rhodococcus wratislaviensis [63] and the Ga operon from Xanthomonas oryzae pv.oryzicola [9,64].While we identified orthologs of certain genes associated with AIB or GA synthesis (or metabolism) in X. fastidiosa or P. phytofirmans, we did not find them organized in operons as originally annotated in the genomes of Rhodococcus wratislaviensis and Xanthomonas oryzae pv.oryzicola, respectively (Supplementary Tables S14 and S15).The chromatograms illustrate the contrast between the AIB levels in Pp, which appear lower than those in the controls, while the AIB levels in Xf consistently remain elevated (Figure 5C).GA was generally detected at low levels across all conditions, including Xf.However, both Xf and ∆rpfF exhibited levels 2 to 3 times higher than the controls.Xf sm and ∆rpfF sm exhibited increases of approximately 100% and 50%, respectively, compared to the culture medium controls.Notably, Pp's GA levels mirror those of the controls (Figure 5C,D).

Exometabolome, Genome, and Transcriptome Integration for X. fastidiosa and P. phytofirmans
We investigated the potential synthesis of AIB and GA by X. fastidiosa or P. phytofirmans through a search for gene clusters for their synthesis.We examined the presence of homologous genes or clusters akin to the reported aib operon from Rhodococcus wratislaviensis [63] and the Ga operon from Xanthomonas oryzae pv.oryzicola [9,64].While we identified orthologs of certain genes associated with AIB or GA synthesis (or metabolism) in X. fastidiosa or P. phytofirmans, we did not find them organized in operons as originally annotated in the genomes of Rhodococcus wratislaviensis and Xanthomonas oryzae pv.oryzicola, respectively (Supplementary Tables S14 and S15).
3.5.Exometabolome, Genome, and Transcriptome Integration for X. fastidiosa and P. phytofirmans Firstly, we used NetSeed [58]-a tool for assessing metabolic network topology and determining a set of exogenously acquired compounds-to assess the metabolic networks of X. fastidiosa and P. phytofirmans, predicting them as profiles in cooperation or competition.
The analysis indicated a low degree of complementarity (<0.32) and a moderate-to-high level of competition (>0.49) between X. fastidiosa and P. phytofirmans.
In a second approach, using MAGI v.1.0[51], we searched for a connection between annotated genes in X. fastidiosa or P. phytofirmans genomes and the detected compounds, respectively, in Xf, Xf sm , ∆rpfF, ∆rpfF sm or Pp (Tables 2 and 3 and Supplementary Tables S4-S8).Out of the 121 identified compounds given as inputs to MAGI, 118 compounds returned results in gene-to-compound connection information from X. fastidiosa or P. phytofirmans.Specifically, all surveyed amino acids (both nEAAs and EAAs) exhibited connections with high reciprocity scores (≥2), i.e., compounds and annotated genes are both linked to the same reactions or metabolic pathways from the consulted databases RheA [65] or MetaCyc [66].CDS coding for amino acid metabolism was previously annotated in the genomes of X. fastidiosa [27] and P. phytofirmans [43].Results for AIB showed a high reciprocity score with 18 and 11 genes in X. fastidiosa and P. phytofirmans, respectively (Supplementary Tables S4-S7).Conversely, GA displayed a notably low reciprocity score with the only two associated genes in X. fastidiosa (PD0286, PD0716) and likewise with the three genes (Bphyt_2421, Bphyt_2594, and Bphyt_0368) in P. phytofirmans (Supplementary Table S8), suggesting a limited likelihood of synthesis by these genes and, respectively, reactions.In P. phytofirmans, Bphyt_5413 is linked to nicotinamide metabolism.In X. fastidiosa, PD0043, PD1071, and PD1494 are linked to biotin metabolism and PD0393 and PD1310, to nicotinic acid metabolism.Next, the resulting genes from the MAGI analysis as described above were used for functional annotation via gene ontology [67,68], through OmicsBox (Supplementary Tables S9-S13).From all the annotated functions, we focused on the top five over-enriched and top five underenriched functions (Figure 6A).Despite minor variations in ∆rpfF, Xf sm , and ∆rpfF smespecially when compared to Xf -over-represented functions generally appeared common across those conditions, such as "ligase activity, forming carbon-nitrogen bonds", "oxidoreductase activity, acting on NAD(P)H", and "ligase activity" (log 2 FC > 5).Functions like "vitamin B6 binding", "pyridoxal phosphate binding", and "tricarboxylic acid cycle" were enriched in Xf, ∆rpfF, Xf sm , and ∆rpfF sm , but underenriched in Pp.Conversely, "CoA-ligase activity" and "fatty acid ligase activity" were highly over-represented in Pp (log 2 FC > 5), whereas they were not enriched in Xf, ∆rpfF, Xf sm , or ∆rpfF sm (Figure 6A).
Finally, we lined up our metabolome data with a previously published RpfF regulon, i.e., transcriptomic data obtained with microarray technology [3].We did not directly and statistically integrate metabolome and transcriptome data, but, instead, we tracked genes directly linked to the detected compounds (obtained from MAGI) and crossreferenced them with the published differentially expressed genes (DEGs) in the RpfF regulon (Figure 6B).A consistent pattern exists where genes downregulated in the ∆rpfF strain were associated with lower levels of metabolites detected in ∆rpfF compared to Xf, including 2-aminoisobutyric acid, cAMP, cGMP, phenylalanine, and tyrosine.In contrast, the upregulated DEGs in the ∆rpfF strain showed no pattern at the metabolite level.Instead, metabolite expression within Xf sm , ∆rpfF, and ∆rpfF sm differed when compared to Xf, signifying a lack of regularity but rather an upregulation of specific metabolites in these conditions.Metabolites cross-checked as upregulated in the metabolome and as a consequence of upregulated genes in the RpfF Regulon include choline, cysteine, 2-deoxyadenosine-5 ′ -monophosphate, guanosine-5 ′ -monophosphate, citraconic acid, asparagine, glutamic acid, threonine, serine, homoserine, uridine-5 ′ -monophosphate, uridine, thymidine, 2-deoxyuridine, cellobiose, and galactose.

Discussion
In prior studies, it was demonstrated that the co-inoculation of P. phytofirmans along with X. fastidiosa leads to a reduction in the symptoms of leaf scorching induced by X. fastidiosa in grapevines.This ameliorative effect was attributed to the activation of plant defence genes, specifically involving pathways related to salicylic acid and ethylene [33].It was further reported that the topical application of P. phytofirmans with a surfactant was efficient and sufficient for PD control in grapevines, even though the protection was spatially limited [34].
In our investigation, we observed that the cultivation of X. fastidiosa in PD3 medium supplemented with 20% spent media from P. phytofirmans (Xf sm and ∆rpfF sm ) resulted in a perturbation of the X. fastidiosa biofilm structure.In particular, the biofilm exhibits an increased volume and extends across the air-liquid medium interface, deviating from the usual compacted ring pattern observed in the liquid culture of Xf.A higher amount of spent media from P. phytofirmans (40%) and different surfaces where X. fastidiosa was grown (plastic versus glass) also interfered drastically in the biofilm formation (Figure 1A,B).Additionally, Baccari and colleagues [33] previously highlighted a synergistic impact arising from the concurrent infection of X. fastidiosa and P. phytofirmans, leading to the activation of plant defence mechanisms.However, in that work, there was no exploration of a direct interaction between P. phytofirmans and X. fastidiosa at the metabolomic level, which could provide insights into their relevance for conventional virulence and plant colonization.
While comparing metabolite profiles, we found a considerable disparity in the number of identified metabolites in Pp in comparison to Xf, Xf sm , ∆rpfF, and ∆rpfF sm .The limited abundance of detected metabolites in Pp, coupled with reduced variability among replicates, contributed to its distinct distribution pattern observed in the principal coordinate analysis (PCoA).While replicates of the Xf exometabolome exhibit close clustering, there is noticeable variability within the Xf sm , ∆rpfF, and ∆rpfF sm exometabolomes, although statistical analyses such as permdisp and permanova confirmed significance across all conditions.The absence of a functional quorum-sensing (QS) system in the ∆rpfF strain was mirrored by its altered metabolite expression profile, akin to findings in other omic analyses such as transcriptomics and outer membrane vesicle (OMV) proteomics [3,69].Notably, the exometabolome profile of Xf sm shows greater similarity to that of ∆rpfF as opposed to Xf.Nevertheless, ∆rpfF sm exometabolomes maintain an overall likeness to ∆rpfF.Consequently, we report that Pp affects the Xf sm exometabolome, although this effect is notably diminished in the impaired DSF-producing strain (∆rpfF sm ).
X. fastidiosa is known to exhibit a hypersecretory trait, particularly regarding virulence proteins and extracellular vesicles (EVs), as established in prior reports [23,24,69,70].Expanding on this, our study targeted the characterization of the X. fastidiosa exometabolome, focusing on low-molecular-weight compounds.Our findings indicate that, in the overall comparison among the five different groups, Xf has a higher expression of compounds.Pairwise comparison of Xf versus Xf sm , ∆rpfF, ∆rpfF sm , or Pp showed a similar pattern, with at least 80% of compounds being upregulated in Xf.It remains an area for further investigation to determine whether the quorum-sensing system interferes with the ability of X. fastidiosa to secrete metabolites.The high secretion of certain metabolites, as will be discussed below, could benefit X. fastidiosa from plant defence mechanisms as the phytopathogen is still not established in a higher population under the influence of the QS system [2].
Another point to highlight is that the ∆rpfF sm exometabolome suggested a less responsive profile to P. phytofirmans spent media than Xf sm .Detected metabolite levels in ∆rpfF/∆rpfF sm were more similar to the levels in Xf /Xf sm (e.g., they had 114 compounds detected in common, from which only two were DEM downregulated in ∆rpfF/∆rpfF sm ).Taken together, these observations indicate that at least partially, the detected responses to P. phytofirmans occur by the disturbance of DSF signaling in Xf.Here, we were limited to metabolites being detected through the LCMS method in a HILIC column.Other coupled MS analyses could complement our findings, e.g., using a c18 column (focused on non-polar and large compounds) or through a surface-based MS technique such as NIMS, which would increase the detected metabolites as seen for other models [71][72][73][74].
From the exometabolomes of Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp, we further studied three sets of compounds: amino acids (EAAs and nEAAs), vitamin B complex, and hormones.The average Gln:Glu ratios were highly similar in Xf and ∆rpfF (0.52 and 0.57, respectively).These values (>0.5) are regarded in the literature as being indicative of N-replete cells [75].The average Gln:Glu ratios in Xf sm and ∆rpfF sm showed an increased trend from treatments without contact with P. phytofirmans spent media, although without statistical significance.Secondly, the overall EAA and nEAA levels were inspected using their chromatogram abundance profiles.EAAs are more abundant than nEAAs in Xf when compared to Xf sm .It is noteworthy that a high ratio of EAAs:nEAAs was reported in the bacteriomes of four xylem-feeding insects [10].We speculate that X. fastidiosa exometabolomes showing especially high EAA content might be an indication of attractiveness to insect vectors.Indeed, studies in the field that simultaneously evaluated the xylem of host plants and insects reported that adults chose plants with high amino acid concentrations [19].The levels of EAAs and nEAAs are significantly different between Xf and ∆rpfF.Another point based on a study conducted by Daugherty and colleagues [76] indicated that the insect prefers plants with lower symptoms than plants with disease symptoms.In fact, the ∆rpfF mutant reaches a systemic infection faster in grapevines compared to the WT strain, and ∆rpfF has a lower fitness for being less transmitted by the insect vector when compared to the WT strain.Moreover, ∆rpfF has a restricted ability to form cell aggregates, which is the successful state of X. fastidiosa when acquired by the insect vector feeding on xylem sap, as reported for the WT strain [41].In this case, it remains to be tested whether the lower level of amino acids secreted by ∆rpfF would add to its worse tendency to be acquired by the insect.P. phytofirmans had very low or totally absent detection of amino acids in its exometabolome under the conditions here tested.Nevertheless, P. phytofirmans has all the predicted pathways for amino acid production according to its genome annotation [43].
Similarly to amino acids, compounds within the vitamin B complex hold particular significance for insect-borne phytopathogens [77].Our metabolomic analysis unveiled the presence of nicotinamide, nicotinic acid, and biotin.Notably, the interplay between nicotinamide and nicotinic acid mirrors the pattern of coding sequences present or absent in X. fastidiosa and P. phytofirmans.The enzyme nicotinamidase is predicted to convert the compound nicotinamide into nicotinic acid (both being isoforms of vitamin B3).In fact, in the X. fastidiosa genome is predicted to lie the enzyme pncA (nicotinamidase), whereas in the P. phytofirmans genome, there is no predicted homolog of a nicotinamidase gene.This conversion was indeed corroborated by our metabolomic analysis (Figure 3, Supplementary Figure S4).Symbiont bacteria, e.g., Wolbachia and Baumannia, have been reported to convert precursors and provide B vitamins to their insect xylem sap-feeding hosts [77].Nevertheless, vitamin B3 is usually not among the B vitamins this symbiont bacteria are able to yield.It remains to be further investigated if the X. fastidiosa provision of this specific isoform of B vitamin increases a host plant's attractiveness towards an insect vector.We also detected biotin, another compound from vitamin B, in Xf sm and ∆rpfF sm .Nevertheless, there was a low level of expression in both exometabolomes and a low level of identification.
AIB, a non-proteinogenic amino acid, exists in two other isomeric forms: β-aminobutyric acid (BABA), known to induce plant disease resistance [78,79], and γ-aminobutyric acid (GABA), a neurotransmitter in animals that is also produced in plants, likely in a signaling role [80].GABA was detected as highly abundant in the xylem sap of grapevines infected with X. fastidiosa [16].AIB is reported as the immediate precursor of ethylene in higher plants.Specifically, AIB inhibits ethylene production because it acts through competitive inhibition of the conversion of 1-aminocyclopropane-1-carboxylic acid to ethylene [81].Very little is known about AIB metabolism in bacteria, although it is reported that the AIB catabolism in Rhodococcus wratislaviensis C31-06 leads to its conversion into α-methyl-D-serine, and after other downstream metabolic conversions, it results in pyruvic acid [63].In turn, gibberellic acid is a well-known plant growth promoter [82].In bacteria, it seems the role of gibberellins is linked to pathogenicity, where secreted gibberellins can act as virulence factors, most probably by suppressing jasmonic acid formation and, ultimately, impairing the host defence response [64], as reported in the rice pathogen Xanthomonas oryzae pv.oryzicola [83].Gibberellic acid (GA), on the other hand, functions as a phytohormone and is synthesized by both plants, fungi, and some bacterial species [84][85][86].
The identity of compounds AIB and GA was confirmed by comparing with standards by MS.Nevertheless, operons for AIB and GA synthesis are absent in X. fastidiosa and P. phytofirmans as seen in other bacteria, and it remains to be further investigated which genes are directly linked with their synthesis in X. fastidiosa or P. phytofirmans.Also, the fact that is well established that these hormones have a function in plant growth and development makes it even more interesting to address the roles of AIB and GA as virulence factors for X. fastidiosa during disease progression.AIB and GA do not seem to have a synergistic effect and therefore it remains unclear how their concomitant secretion would benefit X. fastidiosa.AIB has a higher expression than GA (3,444,408.4 intensity × OD −1 versus 506,103.7 intensity × OD −1 , respectively) and in case its function prevails, it would most probably induce an impairment of host plant growth.X. fastidiosa could take advantage of this by a slower response from the host plant in its vasculature system.Although a GA cluster in Burkholderia species [87], a closely related species, has been described, such a gene cluster was not found in the P. phytofirmans genome (Supplementary Figure S4).
The identification and characterization of changes in specialized metabolites caused by X. fastidiosa in plant hosts are still in the early stages.However, collectively, these findings present a promising outlook on the reprogramming of metabolism following the interaction of plants with the phytopathogen [88].For instance, in infected Leccino plants, higher amounts of salicylic acid were observed compared to Cellina di Nardò plants.The Leccino variety develops milder symptoms compared to those observed on the Cellina di Nardò variety [89].It remains to be investigated if X. fastidiosa interacts with its plant hosts through the specialized metabolites secreted, including AIB and GA.
Reverse ecology analysis revealed generally low complementarity (<0.32) and mediumto-high competition scores (>0.49) for X. fastidiosa towards P. phytofirmans and vice versa as stated in the results section.Therefore, X. fastidiosa and P. phytofirmans are predicted competitive bacterial species.This has implications for the search for a microorganism that can act as a biocontrol for the diseases caused by X. fastidiosa, as seen in the natural occurrence of other endophytes, e.g., M. mesophilicum in orange trees [90].
Overrepresented gene ontology (GO) annotations revealed distinctions not only between Xf and Pp but also for Xf sm , ∆rpfF, and ∆rpfF sm .Collectively, these findings suggest that P. phytofirmans induces a metabolic perturbation in Xf sm , rendering it more similar to ∆rpfF.This observation is further supported by the overall similarity between ∆rpfF sm and ∆rpfF.The specific metabolite or set of metabolites produced by P. phytofirmans, which is not the subject of investigation in this study and appears not to be synthesized by a ∆rpfF ortholog (as reported by Baccari and colleagues [33]), may still be exerting a disruptive effect on quorum sensing.Notably, quorum sensing disruption appears to be an effective strategy for disease suppression.For instance, in a study, both tobacco, a model host, and the orange tree, a natural host that produces DSF, were shown to have increased resistance to X. fastidiosa.The disruption in signaling confounds bacterial behavior and hinders disease development [91].

Conclusions
Overall, we demonstrated that metabolomics can be used for the prospection of X. fastidiosa footprinting.The secretion of EAAs and nEAAs, vitamin B complex, and the ratio between Gln:Glu are influenced by P. phytofirmans spent media and also, by the absence of functional QS signaling.Finally, metabolites with plant metabolism and growth properties, AIB and GA, were detected, which are downregulated in X. fastidiosa when in contact with

Figure 1 .
Figure 1.Biofilm increase in X. fastidiosa in response to exposure to P. phytofirmans secretome.(A) A screening of biofilm formation was performed by measurement of crystal violet staining (Abs595nm).WT strain was grown in a PD3 medium for 7 days with 0.2%, 2%, 20%, and 40% v/v of a 1-day PD3 P. phytofirmans (Pp) spent medium in a 24-well plate.WT strain was grown under same conditions without P. phytofirmans (Pp) spent medium and was set as a reference.(B) Temecula1 WT or ΔrpfF strains were grown in PD3 medium for 7 days with (Xf sm and ΔrpfF sm ) or without (Xf and ΔrpfF) 20% v/v of a 1-day PD3 P. phytofirmans (Pp) spent medium in a glass tube.Biofilm quantification was

Figure 2 .
Figure 2. X. fastidiosa and P. phytofirmans exometabolomes.(A) Venn diagram of Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp exometabolomes.Numbers represent shared or exclusive metabolites detected in the Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp exometabolomes.(B) Principal coordinate analysis (PCoA) of Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp exometabolomes.Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp were cultured in PD3 medium for 7 days without (dark yellow and dark green circles) or with (light yellow and light green circles) 20% supernatant of P. phytofirmans PsJN; Pp was cultured in PD3 medium for 1 day (red circles).Each circle represents three technical replicates with 116 putatively identified metabolites.Permanova and permdisp statistical analyses with 999 permutations were performed, and the statistical significances are shown in the figure.(C) Heatmap representing the abundance of metabolites in the Xf, Xf sm , ∆rpfF, ∆rpfF sm , and Pp exometabolomes.Each one of the 116 compounds detected in this study is indicated in the columns and each respective strain is indicated in the rows.The color code highlights the abundance of the metabolites under the conditions tested, ranging from low (blue) and medium (white) to high (red) expression.An additional column was added with the respective compound class as stated in PubChem.Two main clusters are highlighted, in green and red boxes, as they contain the compounds highly expressed in X. fastidiosa and P. phytofirmans, respectively.A z-score normalization was performed for each compound across the different conditions.Three biological replicates of each species and treatment were used for creating this heatmap.Compounds in black were confirmed by standard run comparison (regarding RT and m/z ppm), compounds in orange had a lower level of confirmation (either or both RT and m/z ppm), and compounds in red could not be confirmed by comparison RT.Low-confidence compounds were therefore considered as putative (possibly tautomers of the compounds assigned).(*) indicates ANOVA statistical significance across conditions.

Figure 6 .
Figure 6.Integration of metabolomic data with functional genomics and RpfF regulon (transcriptomic).(A) GO enrichment analysis performed with genes retrieved from MAGI analysis.GO annotations displayed are the top 5 over-and underrepresented under the 5 conditions tested, and all GO annotations shown passed the Benjamini-Hochberg correction with an FDR < 0.5% criterion.(B) Public data from the RpfF regulon were crossed with genes predicted to be directly involved with the synthesis of the compounds detected in the exometabolomes.The blue-to-red color code indicates down-and upregulated genes in the RpfF regulon.The shapes and sizes of circles indicate the expression levels in peak sizes of compounds detected among the conditions tested.Peak size numbers were normalized into squared z-scores in order to maximize expression level differences among conditions.Dark gray dashed circles around dots or "bubbles" indicate DEMs in relation to Xf.

Figure 6 .
Figure 6.Integration of metabolomic data with functional genomics and RpfF regulon (transcriptomic).(A) GO enrichment analysis performed with genes retrieved from MAGI analysis.GO annotations displayed are the top 5 over-and underrepresented under the 5 conditions tested, and all GO annotations shown passed the Benjamini-Hochberg correction with an FDR < 0.5% criterion.(B) Public data from the RpfF regulon were crossed with genes predicted to be directly involved with the synthesis of the compounds detected in the exometabolomes.The blue-to-red color code indicates down-and upregulated genes in the RpfF regulon.The shapes and sizes of circles indicate the expression levels in peak sizes of compounds detected among the conditions tested.Peak size numbers were normalized into squared z-scores in order to maximize expression level differences among conditions.Dark gray dashed circles around dots or "bubbles" indicate DEMs in relation to Xf.

Table 1 .
List of bacteria and strains used in this work.

Table 2 .
Top scoring genes-to-compounds in X. fastidiosa according to MAGI results.

Table 3 .
Top scoring genes-to-compounds in P. phytofirmans according to MAGI results.