Abstract
Bacterial wilt caused by Ralstonia solanacearum is considered one of the most important diseases that cause economic losses to tomato. Currently, eco-friendly biocontrol agents have been increasingly considered as effective approaches to control tomato bacterial wilt. However, the specific mechanisms by which biocontrol bacteria with distinct functions exert their effects remain unclear. In this study, we employed a combination of amplicon sequencing, transcriptomics, and metabolomics analysis to investigate how Bacillus velezensis and Pseudomonas fluorescens affect the defense responses against R. solanacearum in tomato. We showed that the fermentation broth of these biocontrol agents inhibited the growth of R. solanacearum in vitro, and improves the ability of tomato plants against bacterial wilt. In general, different biocontrol agents protect plants from bacterial wilt in many ways, by recruiting specific microbial communities in rhizosphere soil and activating different synthetic/metabolic and signaling pathways. Collectively, our findings contribute to a more in-depth understanding in disease resistance mechanisms of biocontrol agents, and provide a theoretical foundation for the development of targeted strategies using beneficial microorganisms to suppress disease occurrence.
Introduction
Tomato (Solanum lycopersicum L.), the second-most consumed vegetable crop in the world, encounters many soil microbial infestations which reduce the quality of plants and lead to large losses in production (). Bacterial wilt caused by Ralstonia solanacearum (Rs) is a severely damaging vascular disease that spreads worldwide and infects a wide host range including tomato (). R. solanacearum enters mainly through root wounds, establishes colonies, reaches the xylem through the intercellular space, and eventually fills the entire xylem duct, resulting in necrosis and wilting of the infected plants ().
For many decades, chemical pesticides are commonly used to control crop disease and ensure agricultural harvest. However, the long-term continuous spraying of pesticides will cause pesticide residues of agricultural product, environmental pollution and risk on public health (). To solve this problem, environmentally friendly controls are recommended in order to avoid the adverse effects of chemical pesticides (; ). Among them, biological control by the application of natural antagonists of pathogenic bacteria can be an option for disease management. Many beneficial microbes such as Pseudomonas fluorescens, Bacillus velezensis, and Bacillus amyloliquefaciens, have been documented as effective biocontrol agents against R. solanacearum (; ; ; ).
These beneficial microbes endow to plants with antibacterial peptides (e.g., bacteriocins) and antibacterial metabolites (e.g., antibiotics) that protect plant from pathogenic attacks. Similarly, they also trigger a wealth of molecular processes in plants, including the production of osmolytes like soluble sugars, to increase the plant tolerance to biotic stresses (). Furthermore, many beneficial microbes stimulate host immune responses by priming jasmonic acid (JA)/ethylene (ET) and salicylic acid (SA) signaling pathways in plants, thereby playing an important role in controlling bacterial wilt disease (). For instance, the expression of plant defense marker genes such as PR1, PR2, PR3, PR5, PR12, ERF1, and LoxA in the SA and JA/ET signaling pathways, as well as genes like FLS2 and MKK9 in the MAPK signaling pathway, was significantly up-regulated in infected tomato plants treated with beneficial microbes (; ; ; ).
High-throughput sequencing techniques including genomics, metabolomics, transcriptomics, and proteomics were widespread used to unravel the intricacies involved in plant-pathogen interactions (). Among them, metabolome and transcriptome analysis provide a rapid and efficient way to identify key metabolites and regulatory pathways in plants responding to bacterial wilt (). Using these approaches, glutathione metabolism and phenylpropane pathways of tobacco were demonstrated as primary resistance pathways to R. solanacearum infection, and resistance-related genes responds to R. solanacearum were identified (; ). However, how different biocontrol agents regulate plant-R. solanacearum interaction and their specificity need to be further elucidated, which is crucial for developing effective control strategies.
In this study, we tested the antagonistic activity of biocontrol agents, including B. velezensis and P. fluorescens, against R. solanacearum in vitro and in the field. Multi-omics analysis including microbiome, transcriptome and metabolome was used to elucidated the mechanism underlying the ability of these biocontrol agents to promote the growth of tomato plants and confer resistance to bacterial wilt. The results showed that different biocontrol agents may recruit specific microbial communities and activate specific biological processes to improve the disease resistance of tomato plants. The identification of the pathways and compounds associated with defense responses will not only increase our understanding of the tomato-R. solanacearum interaction affected by biocontrol agents, but also facilitate the improvement of resistance breeding.
Materials and methods
Determining R. solanacearum pathogen inhibition by biocontrol bacteria
In order to study the negative effects of biocontrol bacterial on the pathogen, we compared the reduction in pathogen density when cultured in the fermentation broth of biocontrol bacteria versus when grown alone. We used R. solanacearum strain Rs01 (kindly provided by Prof. Xin Zhang, Institute of Plant Protection, Jiangsu Academy of Agricultural Sciences) as the tomato pathogenic bacterium, and biocontrol bacteria including B. velezensis strain LJBv226 and P. fluorescens strain LJPf01 in our experiments. These bacteria were routinely grown at 30°C in NB medium (glucose 10.0 g/L, peptone 5.0 g/L, yeast extract 0.5 g/L, beef extract 3.0 g/L) with shaking (200 rpm) before all the experiments.
The biocontrol bacteria were cultured in NB medium for 24 h as described above with the addition of centrifugation (10 min at 8000 rpm) and filtration (0.45 μm) steps to prepare the fermentation broth. The R. solanacearum (200 μL, 3.8×108 CFU/mL, OD600 = 2) was inoculated into these two kinds of fermentation broth (10 mL). After 0 h and 24 h of cultivation, pathogen densities were measured as optical density (OD600 (nm)) using a spectrophotometer. The R. solanacearum (200 μL, 3.8×108 CFU/mL, OD600 = 2) was also grown alone in NB medium (10 mL) as the control and all treatments were measured with 12 biological replicates.
Evaluation of the ability of biocontrol agents antagonistic tomato bacterial wilt in the field
The field experiment was performed from August to October 2024 at a tomato field in Huaiyin District of Jiangsu Province, China, where tomato bacterial wilt breaks out every year. The granules, produced by Chengdu Tepu Biotechnology Co., Ltd., contains 2×108 CFU/g B. velezensis or 5×108 CFU/g P. fluorescens. The tomato seedlings were transplanted into the field on August, 4. Before the onset of tomato bacterial wilt, The biocontrol agents were diluted to 1:150, 1:300 and 1:500, respectively, and sprayed to the roots (50 mL per tomato plant) for the first time (September, 12). and then twice again at the flowering stage (September, 15 and 18). Finally, the disease incidents, biocontrol efficiency, plant height, root length, stem diameter, and chlorophyll content (SPAD) were investigated at 15 d after the application of biocontrol agents (September, 27). Tomato plants without biocontrol agents were used as control. Disease incidence and biocontrol efficiency were calculated by formula as follows:
Disease incidence (%) = (Total number of infected plants/Total number of plants) ×100;
Biocontrol efficiency (%) = (Disease incidence of control - Disease incidence of treatment)/Disease incidence of control ×100
Rhizosphere sample collection and DNA extraction
Rhizosphere soil samples were collected 15 days after the application of biocontrol regents, according to the method described by Wei et al. with some modifications (). Specifically, tomato plants were carefully uprooted and the loose soil was shaken off, and then the soil attached to the root surface (within approximately 0–4 mm) was collected as rhizosphere soil. After removing the debris and residual fine roots, the soil samples collected in the field were sealed in sterile self-sealed bags, immediately placed in an ice box, and stored in a -80°C freezer until the determination of soil community structure.
Soil (0.25 g) from each sample was used for DNA extraction using the DNeasy PowerSoil Kit (QIAGEN, Inc., Netherlands), following the manufacturer’s instructions. The extracted DNA was checked on 1.2% (w/v) agarose gels, and quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).
Illumina Novaseq high-throughput sequencing and statistical analysis for microbial communities
PCR was performed to amplify the bacterial 16S rRNA genes V3–V4 region using the forward primer 338F (5’-ACTCCTACGGGAGGCAGCA-3’) and the reverse primer 806R (5’- GGACTACHVGGGTWTCTAAT-3’) () or amplify the fungal internal transcribed spacer (ITS) sequences using the forward primer ITS1 (5’-GGAAGTAAAAGTCGTAACAAGG-3’) and the reverse primer ITS2 (5’-GCTGCGTTCTTCATCGATGC-3’) (). Sample-specific 7-bp barcodes were incorporated into the primers for multiplex sequencing. The PCR components contained 5 μL of Q5 reaction buffer (5×), 5 μL of Q5 High-Fidelity GC buffer (5×), 0.25 μL of Q5 High-Fidelity DNA Polymerase (5U/μL), 2 μL (2.5 mM) of dNTPs, 1 μL (10 μM) of each Forward and Reverse primer, 2 μL of DNA Template, and 8.75 μL of ddH2O. Thermal cycling consisted of initial denaturation at 98°C for 2 min, followed by 25 cycles consisting of denaturation at 98°C for 15 s, annealing at 55 °C for 30 s, and extension at 72°C for 30 s, with a final extension of 5 min at 72°C. PCR products were measured through agarose gel electrophoresis, and purified with Vazyme VAHTSTM DNA Clean Beads. Library construction was performed according to the standard procedure of TruSeq Nano DNA LT Library Prep Kit for Illumina®. The constructed library was sequenced on an Illumina Nova 6000 platform at Irun Biological Technology Co., Ltd. (Shanxi, China), using Miseq.
Adaptors and primer sequences were removed using Cutadapt and FASTP (version 0.14.1) software to obtain clean reads. The raw sequences were then denoised by the use of DADA2. The Quantitative Insights Into Microbial Ecology (QIIME, version 2) pipeline was employed to process the sequencing data, as previously described (). Briefly, raw sequencing reads with exact matches to the barcodes were assigned to respective samples and identified as valid sequences. The low-quality sequences were filtered through following criteria (): sequences that had a length of <150 bp, sequences that had average Phred scores of <20, sequences that contained ambiguous bases, and sequences that contained mononucleotide repeats of >8 bp. Paired-end reads were assembled using FLASH (). After chimera detection, the remaining high-quality sequences were clustered into operational taxonomic units (OTUs) at 97% sequence identity by UCLUST (). A representative sequence was selected from each OTU using default parameters. OTU taxonomic classification was conducted by BLAST searching the representative sequences set against the Greengenes Database () using the best hit (). An OTU table was further generated to record the abundance of each OTU in each sample and the taxonomy of these OTUs. OTUs containing less than 0.001% of total sequences across all samples were discarded. To minimize the difference of sequencing depth across samples, an averaged, rounded rarefied OTU table was generated by averaging 100 evenly resampled OTU subsets under the 90% of the minimum sequencing depth for further analysis.
For bacterial 16S rRNA genes, the representative sequences were annotated against the SILVA database (http://www.arb-silva.de); for fungal ITS sequences, the representative sequences were annotated against the UNITE database (https://unite.ut.ee/). OTU-level alpha diversity indices, such as Chao1 richness estimator, Shannon diversity index, and Simpson index, were calculated using the OTU table in QIIME. Beta diversity analysis was performed to investigate the structural variation of microbial communities across samples using UniFrac distance metrics () and visualized via principal coordinate analysis (PCoA). A comparison of overall microbial distribution among samples was conducted on the relative abundances of phyla or genus using OTUs based on taxonomy. All the above analyses, including the community heatmap at the genus level, were performed using R software (v.3.5.3).
Transcriptome analysis
Root samples 15 d post the application of biocontrol agents at the concentration of 1:150 were collected, with three independent biological replicates per sample. Total RNA from each sample was extracted using the Total RNA Extraction Kit (AM1561, Ambion), following the manufacturer’s instructions. The quality of RNA was examined using Agilent Bioanalyzer 2100 system. A total of 1.5 μg of RNA per sample was used as the input materials to create sequencing libraries using NEBNext Ultra Directional RNA Library Prep kit (NEB, Beverly, MA, USA), and sequencing was performed at Biotree Technology Co. Ltd. (Shanghai, China) on an Illumina Novaseq 6000 platform. Low-quality reads and adapter sequences were removed from raw reads through fastp software before assembly. The sequencing clean reads were remapped to the reference sequence using Hisat2 v2.2.1 software.
DESeq2 v1.42.0 software was used to compute the differential expression post the application of biocontrol agents compared to tomato plants without treatment. Genes with a minimal 2-fold difference in expression (|log2 (fold change)| ≥ 1) and Padj ≤ 0.05 were considered as differential expressed genes (DEGs), and DEGs were used for the following KEGG(Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg/) pathway-enrichment analysis with clusterProfiler v4.8.1 software. Pathways with a corrected p value of <0.05 are considered to be significantly enriched.
Metabolome analysis
Root samples were harvested at 15 d post the application of biocontrol agents at the concentration of 1:150, and non-treated samples were used as control. Three biological replicates for each treatment were evaluated. The samples (20 mg ± 1 mg) were lyophilized, mixed with beads and 1000 μL of extraction solution (MeOH: ACN: H2O, 2:2:1 (v/v)). The extraction solution contains deuterated internal standards. The mixed solution was vortexed for 30 s. Then the mixed samples were homogenized (35 Hz, 4 min) and sonicated for 5 min in 4°C water bath, the step repeat for three times. The samples were incubated for 1 h at -40°C to precipitate proteins. Then the samples were centrifuged at 12000 rpm for 15 min at 4°C. Transfer 400 μL of liquid to the well of a protein precipitation plate. Place the plate on the manifold. Apply vacuum, 6 psi, 120 s. Take plate from the positive pressure device for analysis. The quality control (QC) of sample was prepared by mixing an equal aliquot of the supernatant of samples.
Metabolites were then analyzed through untargeted LC-MS/MS technology. Briefly, extraction of samples (2 μL) was injected into an Ultra Performance Liquid Chromatography (UPLC, Vanquish, Thermo Fisher Scientific) equipped with a tandem mass spectrometry (MS/MS, Orbitrap Exploris 120) bound with Waters ACQUITY UPLC BEH Amide column (2.1 mm × 100 mm, 1.7 μm). For polar metabolites, the mobile phase consisted of water (25 mmol/L ammonium acetate and 25 mmol/L ammonia water, A) and acetonitrile (B). For non-polar metabolites, the mobile phase consisted of water (0.01% acetic acid, A) and isopropyl alcohol: acetonitrile (1:1, v/v, B).
The detailed parameters of the mass spectrum are as follows: sheath gas flow rate, 50 Arb; aux gas flow rate, 15 Arb; capillary temperature, 320°C; full ms resolution, 60000; MS/MS resolution, 15000; collision energy, SNCE 20/30/40; spray voltage, 3.8 kV (positive) or -3.4 kV (negative). The raw data were processed using ProteoWizard software. The metabolites were identified by matching the secondary spectral information with the BiotreeDB (V3.0) database, and quantified by integrating peak areas. The multivariate analysis of identified metabolites was performed by Principal Component Analysis (PCA).
Transcriptome and metabolome correlation analysis
Differentially expressed genes (DEGs) mapped to the “Plant-pathogen interaction” pathway and differentially accumulated metabolites (DAMs) in pathways related to plant immunity were selected for correlation analysis. The correlation analysis was performed using Pearson’s correlation coefficient in the R language corr.test package to calculate the correlation coefficient and the p-value. The results were visualized by R language package.
Results
Bacillus velezensis or P. fluorescens had significant inhibitory effect on the growth of R. solanacearum
We tested the inhibitory effect of biocontrol agents on the growth of R. solanacearum in vitro. When R. solanacearum was cultured alone in NB medium for 24 h, its density reached to OD600 = 1.288. After cultured with fermentation broth of B. velezensis or P. fluorescens, the density of R. solanacearum significantly decreased to 0.755 or 0.707 (Figure 1), respectively, indicating the biocontrol potential of these beneficial bacteria against the R. solanacearum pathogen. There was no significant difference between these two kinds of fermentation broth in inhibiting the growth of R. solanacearum.
Figure 1
Bacillus velezensis or P. fluorescens effectively control tomato bacterial wilt and promote plant growth in the field
In order to explore the effect of controlling tomato bacterial wilt in the field, B. velezensis (2×108 CFU/g) or P. fluorescens (5×108 CFU/g) were diluted to 1:500, 1:300 and 1:150, and applied to tomato plants before the onset of tomato bacterial wilt and at the flowering stage, and then the control effect of disease was investigated at 15 d after the treatment. The results of the field test showed that all treatments significantly reduced disease severity of tomato bacterial wilt (Figure 2A). Compared with non-treated control plants, the biocontrol efficacy of B. velezensis or P. fluorescens of 1:500 reached 32.69% and 43.42%, respectively. With the increase of concentration, the biocontrol efficiency was significantly improved, with P. fluorescens of 1:150 reached 76.08% and B. velezensis of 1:150 reached 68.40% (Figure 2B).
Figure 2
Correspondingly, the application of these biocontrol agents significantly promoted the growth of tomato plants compared to control plants. The promotion effect was also gradually enhanced with the increase of concentration. When the dilution concentration of P. fluorescens reached 1:150, the average plant height, root length and stem thickness significantly increased from 96.47 cm, 13.14 cm and 9.06 mm to 105.49 cm, 15.48 cm and 11.75 mm, respectively. Meanwhile, B. velezensis of 1:150 showed similar effect on plant growth, which increased plant height, root length and stem diameter to 104.60 cm, 15.21 cm and 11.50 mm, respectively (Figures 2C–E).
Pathogen infection often leads to a decrease in photosynthetic assimilates production, which can be analyzed by monitoring chlorophyll content (SPAD value) (). Studies have shown that the chlorophyll content of plant leaves can be significantly increased after inoculation with biocontrol bacteria (; ). Compared with control plants, B. velezensis or P. fluorescens of 1:150 increased chlorophyll content from 2.41 SPAD to 2.98 SPAD or 3.00 SPAD, respectively, which is also concentration-dependent (Figure 2F), playing an important role in increasing the efficiency of photosynthesis in plants.
Therefore, these results indicate that biocontrol bacteria, such as B. velezensis and P. fluorescens, significantly improve the ability of tomato plants against bacterial wilt, thus promote their growth and development.
The composition of bacterial communities in rhizosphere soil changed significantly by biocontrol agents
Current studies mainly focus on the effect of bacterial wilt on the community composition in rhizosphere soil of multiple economically important crops, such as tobacco, pepper, and potato (; ). However, how this effect is influenced by biocontrol agents remains unclear. The diversity of bacterial communities in the rhizosphere soil of infected tomato plants treated with or without biocontrol agents was analyzed, and more than 60,000 high-quality 16S rRNA sequences were obtained for each sample. Principal coordinates analysis (PCoA) was conducted to confirm the differences among samples (Figure 3A). The degree of explanation of principal coordinate 1 (PC1) and PC2 for the soil sample differences was 16.9% and 14.3%, respectively. The results revealed that bacterial communities of different treatments were clustered, suggesting that B. velezensis or P. fluorescens had a great impact on the composition of bacterial communities.
Figure 3
Alpha diversity of bacterial communities was evaluated to more comprehensively understand the effects of these biocontrol agents on community structure. The Chao 1 results showed that the richness of bacterial species was not significantly changed among all treatments. By contrast, the Simpson index had a significant difference (p<0.05), and the Shannon diversity was also different although not significantly (p=0.081), indicating that B. velezensis or P. fluorescens affected the bacterial diversity in the rhizosphere of infected tomato (Figure 3B). Given that B. velezensis or P. fluorescens of 1:150 demonstrated a more pronounced improvement in plant growth and disease resistance, we further directly compared the bacterial community composition between the control and groups treated with biocontrol agents of high concentration. Similarly, no significant differences were observed in Chao 1 and Shannon indices, whereas the Simpson index was lower in groups treated with B. velezensis or P. fluorescens of 1:150, indicating a significantly reduced bacterial diversity (p=0.05, Supplementary Figure S1).
At the phylum level, Proteobacteria, Chloroflexi, Acidobacteriota, Actinobacteriota, and Gemmatimonadota were the dominant bacteria in all treatments, accounting for more than 75% of the total amount of bacterial. The community composition of these treatments was similar, but the relative abundance of each phylum did differ in the rhizosphere soil samples. Specifically, the relative abundance of the Proteobacteria, which contains the pathogenic species of bacterial wilt disease, was significantly decreased to 31.4%-34.1% in the groups treated with B. velezensis or P. fluorescens, compared to 36.7% in the group without biocontrol agent treatment (Figure 3C).
According to the relative abundance of each sample at the genus level, the top 50 taxon were selected to generate a bacterial community heatmap (Figure 3D). Compared with control, in which Ellin6067, Sphingomonas, MND1, TRA3-20, CCD24, mle1–7 and Rokubacteriales were more abundant, the effects of different concentrations and types of biocontrol agents on the composition of the bacterial community were not exactly the same. For example, the genera Gitt-GS-136, Iamia, Streptomyces, AKYG1722, and JG30-KF-CM45 were significantly enriched in groups treated with B. velezensis of 1:500, while the genera S0134 terrestrial group, Pseudoxanthomonas, KD4-96, IMCC26256, and Cupriavidus were significantly more abundant in groups treated with B. velezensis of 1:300. In contrast, Nocardioides, Luteimonas, Saccharimonadales, and Gemmatimonas were significantly enriched in groups treated with P. fluorescens of 1:300. These results suggest that the effects of different biocontrol agents on rhizosphere bacterial may not be identical, even though they all reduced the severity of tomato bacterial wilt and promoted plant growth.
Biocontrol agents had less effect on fungal community composition in rhizosphere soil of infected tomato plants
The involvement of soil fungus in antagonistic interaction with soil-borne pathogens have been comprehensively and precisely illustrated by previous studies (; ; ), thus we also examined the effect of biocontrol agents on fungal community structure in rhizosphere soil of infected tomato plants. After quality control, more than 80,000 high-quality fungal ITS sequences were obtained for each sample. Concerning fungal communities, the degree of explanation of PC1 and PC2 was only 6.7% and 6.8% (Figure 4A), indicating that sample discrimination was not obvious. Similarly, no obvious differences in species richness according to Chao 1 indices or species diversity according to Simpson and Shannon indices were observed among these samples (Figure 4B). Furthermore, no significant differences were observed in Chao1, Simpson, and Shannon indices between the control and groups treated with the highest concentration (Supplementary Figure S2), suggesting that the effect of these biocontrol agents on fungal community structure in rhizosphere soil was less than that on bacterial community.
Figure 4
However, the relative abundance of unique genera was different among samples. The application of biocontrol agents led to a remarkable enrichment of potentially beneficial genera (e.g., Lophotrichus, Botryotrichum), in stark contrast to the depletion of fungal pathogens (e.g., Ophiocordyceps, Plectosphaerella) (Figure 4C) (; ; ; ), emphasizing the potential significance of these biocontrol agents in preserving the health of tomato plants through specific fungal genera.
Moreover, different concentrations of biocontrol agents produce distinct dominant species, and fungal species affected by B. velezensis and P. fluorescens were also different (Figure 4D). For example, Ascobolus and Diversispora were enriched in the group treated with B. velezensis of 1:150, while Mendogia and Keratinophyton were dominant in the group treated with P. fluorescens of 1:150.
Biocontrol agents lead to transcriptional reprogramming in tomato infected with bacterial wilt
Considering that the high concentration (1:150) of biocontrol agents has a more obvious improvement in biocontrol efficiency and plant growth (Figure 2), transcriptome analysis was conducted using roots of infected tomato with or without biocontrol agents (1:150), to understand their underlying molecular mechanisms against tomato bacterial wilt. The number of clean reads for each sample ranged from 45.58 million to 52.50 million. Each library generated >6 Gb of clean bases. The percentage of Q20 was >97% and the percentage of Q30 was >93%.
Bacillus velezensis or P. fluorescens with high concentration induced intense transcriptome reprogramming in tomato plants (Figure 5A). Results showed that B. velezensis induced 3965 differentially expressed genes (DEGs), including 2273 up- and 1692 down-regulated DEGs (Figure 5B; Supplementary Table S1). In contrast, P. fluorescens induced 3709 DEGs, including 2025 up- and 1684 down-regulated genes (Figure 5B; Supplementary Table S2). Among them, 1886 DEGs are overlapped between different treatment, and others are biocontrol agents-specific (Figure 5C).
Figure 5
To investigate the biological processes and molecular functions that changes after the application of biocontrol agents, DEGs were mapped to the KEGG pathway database simultaneously under the criteria |log2 (fold change) | ≥ 1 and Padj ≤ 0.05. The top 20 terms were selected to draw scatter plots for display.
For the group treated with B. velezensis of 1:150, many biosynthesis and metabolism processes including “Phenylpropanoid biosynthesis”, “Arginine and proline metabolism”, “Alanine, aspartate and glutamate metabolism”, “alpha-Linolenic acid metabolism”, “Nitrogen metabolism” and “Biosynthesis of unsaturated fatty acids” were significantly enriched (Figure 6A), many of which play important roles in plant immunity (; ; ; ). In addition, the “Plant-pathogen interaction” pathway was also significantly upregulated in the sample (Figure 6A). These results suggest that high concentrations of this biocontrol agent activate a variety of metabolic and signaling pathways in tomato plants to help them increase the resistance against R. solanacearum, which is consistent with its significant improvement in disease resistance and plant growth (Figure 2).
Figure 6
In contrast, in addition to the same pathways activated by B. velezensis of 1:150, some other metabolism processes such as “Flavonoid biosynthesis” which are very important for plant defense against pathogens (), were also enriched in the group treated with P. fluorescens of high concentration (Figure 6B). Notably, “Plant hormone signal transduction” was activated by this treatment (Figure 6B). Previous studies revealed that salicylic acid (SA) and Jasmonate (JA) are major defense-related phytohormones, other phytohormones such as ethylene (ET), auxin, gibberellins (GAs), abscisic acid (ABA), brassinosteroids (BRs), and cytokinins (CKs) are also involved in defense responses (). We suppose that plant hormone signal transduction may be specifically activated by this beneficial microorganism, and thus regulate plant disease resistance coordinating with biosynthesis/metabolism and plant-pathogen interaction pathways, although it is not known which plant hormone plays a major role in this process.
Overall, the application of these biocontrol agents led to dramatic transcriptome reprogramming in tomatoes infected with bacterial wilt. Different types of biocontrol agents may improve plant resistance by regulating specific biosynthesis/metabolism processes or signaling pathways to maintain the balance between plant growth and biological stress.
Specific metabolites play an important role in the regulation of tomato bacterial wilt by biocontrol agents
To gain the overview of metabolites affected by these biocontrol agents in infected tomato plants, an untargeted metabolomics by LC-MS/MS was performed to evaluated metabolite profiles in detail. Similarly, the roots of infected tomatoes treated with high concentrations (1:150) of B. velezensis or P. fluorescens were analyzed. Metabolites with p value of <0.05 and the variable importance in project (VIP) > 1 were considered to be differentially accumulated metabolites (DAMs) (Supplementary Tables S3, S4), and the top 10 DAMs significantly up-regulated and down-regulated were selected for display (Figure 7).
Figure 7
Both of biocontrol agents resulted in a reduction of Ginkgolide J and the 2-Hydroxy-3-(3,4-dihydroxyphenyl)propenamide, whose superoxide scavenging effect contributes to the antioxidant properties (). Metabolite profiles revealed that specific metabolites were enriched or reduced after the application of these biocontrol agents. Particularly, 1-18:3-2-18:3-phosphatidylethanolamine, 8-amino-7-(carboxyamino)nonanoate, isopyridoxal and cholestan-3-one were significantly more abundant in the group treated with B. velezensis, while deacetylasperuloside acid methyl ester, daidzin, mycotoxin T-2 and typhaneoside et al. were enriched in the group treated with P. fluorescens (Figure 7). These results suggested that different biocontrol agents may induce the production of specific metabolites in infected tomato plants to resist bacterial wilt.
Metabolites with p value of <0.05 and the variable importance in project (VIP) > 1 were considered to be differentially accumulated metabolites (DAMs), and the top 10 DAMs significantly up-regulated and down-regulated were selected for display. Red and blue terms represent the up- and down-regulated metabolites, respectively.
Key pathways that are most relevant to the metabolite differences were analyzed, and results showed that “Pantothenate and CoA biosynthesis” has the highest correlation with different metabolites in the group treated with B. velezensis of 1:150, while “Phenylalanine metabolism” was the highest when treated with P. fluorescens of high concentration (Supplementary Figure S3), further clarifying the differences in the effects of these two biocontrol agents on the metabolism of infected tomato plants.
Correlation analysis of transcriptome and metabolome data
To further reveal the differences between each biocontrol bacteria as well as their impact on plant innate immunity in greater depth, we analyzed the correlation between DEGs mapped to the “Plant-pathogen interaction” pathway and DAMs in the metabolic pathways that may be involved in plant immunity. There are 51 and 53 DEGs mapped to the “Plant-pathogen interaction” pathway in the group treated with B. velezensis or P. fluorescens of 1:150, respectively (Supplementary Tables S5, S6). Among pathways significantly relevant to the metabolite differences, “Pantothenate and CoA biosynthesis”, “Glycine, serine and threonine metabolism”, and “Phenylalanine metabolism” have been previously reported to be involved in the defense-related responses (Supplementary Figure S3) (; ; ). Therefore, 14 and 12 DAMs mapped to these metabolic pathways in the group treated with B. velezensis or P. fluorescens, respectively, were selected for correlation analysis (Supplementary Tables S7, S8).
The results show that most of genes related to plant-pathogen interaction were positively correlated with pantothenic acid, 5, 10-methylene-THF, betaine, beta-alanine, 3-ureidopropionic acid, choline, alpha-ketoisovaleric acid and glycine in the group treated with B. velezensis. In contrast, except for 5, 10-methylene-THF and choline, most genes were only positively correlated with 2-phenylacetamide when treated with P. fluorescens (Figure 8). Although the “Plant-pathogen interaction” pathway was up-regulated in both groups, they might promote the production of different immune-related metabolites to promote tomato-R. solanacearum interaction.
Figure 8
Discussion
Biological control has become an increasingly important means to inhibit plant diseases due to its green and efficient advantages. Our understanding of the ways in which biocontrol agents protect plants from disease has developed considerably in recent years, which is crucial to the development of biocontrol strategies for plant diseases (). This study revealed that biocontrol agents, including B. velezensis and P. fluorescens, significantly inhibited the growth of R. solanacearum in vitro (Figure 1). At the same time, tomato bacterial wilt that is caused by R. solanacearum was efficiently controlled by the biocontrol agents, promoting the growth of tomato plants in the field (Figure 2), which makes these biocontrol agents an important strategy to resist tomato bacterial wilt.
We found that the composition of bacterial community in rhizosphere soil changed more significantly than that of fungal community after the application of biocontrol agents, with more abundant beneficial microorganisms and reduced harmful microorganisms (Figures 3, 4). In order to further investigate how these biocontrol agents improve the disease resistance of plants and whether their regulatory mechanisms are specific, we analyzed the transcriptome and metabolome changes of infected tomato roots. Results showed that B. velezensis activated the synthetic/metabolic pathway and plant-pathogen interaction pathway in infected tomato plants, while P. fluorescens also activated some other pathways such as plant hormone signal transduction (Figures 5, 6). Among them, secondary metabolites derived from biosynthetic pathways such as those for flavonoids and phenylpropanoid are essential in plant innate immunity and response to environmental stresses (; ). Further correlation analysis of transcriptome and metabolome data revealed the differences between each biocontrol bacteria as well as their impact on plant innate immunity (Figure 8).
Hormone-regulated genes play an essential role in regulating plant response to biotic and abiotic stress. In general, plant defense responses include SA-dependent defenses acting on biotrophs and JA/ET-dependent responses acting on necrotrophs (). Some beneficial microbes could trigger systemic resistance by regulating SA- and JA/ET-mediated defense genes, such as pathogenesis-related (PR) genes (e.g., PR1, PR2) and ethylene-responsive factors (ERFs), strengthening plant immunity against tomato bacterial wilt (; ). Our findings demonstrate that P. fluorescens activates plant hormone signal transduction, whereas B. velezensis does not. Further validation is required to determine which specific plant hormone signaling pathway and its associated genes are activated by P. fluorescens to enhance resistance against tomato bacterial wilt.
Previous studies have shown that B. velezensis could produce antibacterial substances such as lipopeptides and polyketides at ecological sites to exert antagonistic effects, inhibiting motility traits of R. solanacearum and damaging the cell wall (; ). In contrast, the 2,4-diacetylphloroglucinol (DAPG) from the crude metabolites of P. fluorescens has the antimicrobial potential against R. solanacearum (). Specific reactive oxygen species and hormone signaling cascades induced by P. fluorescens may also contribute to its role as beneficial microorganism (). We speculate that the differences in the regulatory mechanisms of these two kinds of biocontrol agents may be related to the substances they produce, although their biocontrol efficiency are similar. These biocontrol agents may recruit specific microorganisms in rhizosphere soil by producing different metabolites, and activated specific immune response through different pathways, thus protect tomato from bacterial wilt. Our study provides new insight into the specific mechanisms of different biocontrol agents against R. solanacearum, and establishes a theoretical basis for developing compound biocontrol agents.
Statements
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 authors.
Author contributions
JL: Writing – review & editing, Supervision, Funding acquisition, Resources. XD: Data curation, Investigation, Writing – original draft. TS: Writing – original draft, Investigation. WX: Data curation, Writing – original draft. TZ: Data curation, Writing – original draft. QW: Writing – original draft, Data curation. PG: Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by Modern agricultural industry technology research and development and promotion application project (Grant number: 23H010101238).
Conflict of interest
Author W-LX was employed by the company Chengdu Topu Biotechnology Co., LTD. Author TZ was employed by the company Chengdu Juzhuo Biotechnology Co., LTD. Author QW was employed by the company Sichuan Taikang BiotechnologyCo., LTD.
The remaining 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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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/fpls.2025.1620460/full#supplementary-material
References
1
AltschulS. F.MaddenT. L.SchäfferA. A.ZhangJ.ZhangZ.MillerW.et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res.25, 3389–3402. doi: 10.1093/nar/25.17.3389
2
BastíasD. A.BalestriniR.PollmannS.GundelP. E. (2022). Environmental interference of plant–microbe interactions. Plant Cell Environ.45, 3387–3398. doi: 10.1111/pce.14455
3
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. doi: 10.1038/nmeth.f.303
4
ChenH.JiangW. (2014). Application of high-throughput sequencing in understanding human oral microbiome related with health and disease. Front. Microbiol.5, 508. doi: 10.3389/fmicb.2014.00508
5
ClaessonM. J.O’SullivanO.WangQ.NikkiläJ.MarchesiJ. R.SmidtH.et al. (2009). Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PloS One4, e6669. doi: 10.1371/journal.pone.0006669
6
DermastiaM. (2019). Plant hormones in phytoplasma infected plants. Front. Plant Sci.10, 477. doi: 10.3389/fpls.2019.00477
7
DeSantisT. Z.HugenholtzP.LarsenN.RojasM.BrodieE. L.KellerK.et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microb.72, 5069–5072. doi: 10.1128/AEM.03006-05
8
DixonR. A. (2001). Natural products and plant disease resistance. Nature411, 843–847. doi: 10.1038/35081178
9
DongH.GaoR.DongY.YaoQ.ZhuH. (2023). Bacillus velezensis RC116 inhibits the pathogens of bacterial wilt and fusarium wilt in tomato with multiple biocontrol traits. Int. J. Mol. Sci.24, 8527. doi: 10.3390/ijms24108527
10
EdgarR. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics26, 2460–2461. doi: 10.1093/bioinformatics/btq461
11
HintonM. J.ParryD. W. (1993). Screening selected fungi for antagonism towards Pseudocercosporella herpotrichoides (fron) deighton, the cause of eyespot disease of cereals. Biocontrol Sci. Techn.3, 13–19. doi: 10.1080/09583159309355254
12
IraniS.ToddC. D.WeiY.Bonham-SmithP. C. (2019). Changes in phenylpropanoid pathway gene expression in roots and leaves of susceptible and resistant Brassica napus lines in response to Plasmodiophora brassicae inoculation. Physiol. Mol. Plant P.106, 196–203. doi: 10.1016/j.pmpp.2019.02.007
13
JamilA.MusheerN.KumarM. (2021). Evaluation of biocontrol agents for management of wilt disease of tomato incited by Fusarium oxysporum f. sp. lycopersici. Arch. Phytopathol. Plant Prot.54, 1722–1737. doi: 10.1080/03235408.2021.1938353
14
KimB.ParkA. R.SongC. W.SongH.KimJ. C. (2022). Biological control efficacy and action mechanism of Klebsiella pneumoniae JCK-2201 producing Meso-2,3-butanediol against tomato bacterial wilt. Front. Microbiol.13, 914589. doi: 10.3389/fmicb.2022.914589
15
KobmooN.MongkolsamritS.WutikhunT.TasanathaiK.KhonsanitA.ThanakitpipattanaD.et al. (2015). New species of Ophiocordyceps unilateralis, an ubiquitous pathogen of ants from Thailand. Fungal Biol.119, 44–52. doi: 10.1016/j.funbio.2014.10.008
16
KongH. G.SangM. K.AnJ. H.KimS.JinY. J.SongJ. (2022). Changes in the composition and microbial community of the pepper rhizosphere in field with bacterial wilt disease. Plant Pathol. J.38, 692–699. doi: 10.5423/PPJ.NT.09.2022.0136
17
LeK. D.KimJ.NguyenH. T.YuN. H.ParkA. R.LeeC. W.et al. (2021). Streptomyces sp. JCK-6131 protects plants against bacterial and fungal diseases via two mechanisms. Front. Plant Sci.12, 726266. doi: 10.3389/fpls.2021.726266
18
LiQ.LiuZ.JiangZ.JiaM.HouZ.DouD.et al. (2025). Phenylalanine metabolism-dependent lignification confers rhizobacterium-induced plant resistance. Plant Physiol.197, kiaf016. doi: 10.1093/plphys/kiaf016
19
LiY.WangL.SunG.LiX.ChenZ.FengJ.et al. (2021). Digital gene expression analysis of the response to Ralstonia solanacearum between resistant and susceptible tobacco varieties. Sci. Rep.11, 3887. doi: 10.1038/s41598-021-82576-8
20
LiZ.ZuC.WangC.YangJ.YuH.WuH. (2016). Different responses of rhizosphere and non-rhizosphere soil microbial communities to consecutive Piper nigrum L. monoculture. Sci. Rep.6, 35825. doi: 10.1038/srep35825
21
LozuponeC.KnightR. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microb.71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005
22
MaH.ShengC.QiaoL.ZhaoH.NiuD. (2020). A comparative proteomic approach to identify defence-related proteins between resistant and susceptible rice cultivars challenged with the fungal pathogen Rhizoctonia solani. Plant Growth Regul.90, 73–88. doi: 10.1007/s10725-019-00551-w
23
MagočT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics27, 2957–2963. doi: 10.1093/bioinformatics/btr507
24
NicolaS.TibaldiG.FontanaE. (2009). Tomato production systems and their application to the tropics. Acta Hortic.821, 27–33. doi: 10.17660/ActaHortic.2009.821.1
25
O’BrienP. A. (2017). Biological control of plant diseases. Australas. Plant Pathol.46, 293–304. doi: 10.1007/s13313-017-0481-4
26
PangZ.ChenJ.WangT.GaoC.LiZ.GuoL.et al. (2021). Linking plant secondary metabolites and plant microbiomes: A review. Front. Plant Sci.12, 621276. doi: 10.3389/fpls.2021.621276
27
RamlalA.AnitaR.AparnaN.CharuK.ReetaK.JitendraK.et al. (2023). Importance of omics approaches in plant-microbe interactions in plant disease control. Physiol. Mol. Plant P.128, 102153. doi: 10.1016/j.pmpp.2023.102153
28
RazaW.LingN.LiuD.WeiZ.HuangQ.ShenQ. (2016). Volatile organic compounds produced by Pseudomonas fluorescens WR-1 restrict the growth and virulence traits of Ralstonia solanacearum. Microbiol. Res.192, 103–113. doi: 10.1016/j.micres.2016.05.014
29
ScholtyssekH.DamerauW.WesselR.SchimkeI. (1997). Antioxidative activity of ginkgolides against superoxide in an aprotic environment. Chem-Biol. Interact.106, 183–190. doi: 10.1016/S0009-2797(97)00067-7
30
SharmaA.AbrahamianP.CarvalhoR.ChoudharyM.ParetM. L.ValladG. E.et al. (2022). Future of bacterial disease management in crop production. Annu. Rev. Phytopathol.60, 259–282. doi: 10.1146/annurev-phyto-021621-121806
31
ShiH.XuP.YuW.ChengY.DingA.WangW.et al. (2022). Metabolomic and transcriptomic analysis of roots of tobacco varieties resistant and susceptible to bacterial wilt. Genomics114, 110471. doi: 10.1016/j.ygeno.2022.110471
32
Siegel-HertzK.Edel-HermannV.ChapelleE.TerratS.RaaijmakersJ. M.SteinbergC. (2018). Comparative microbiome analysis of a fusarium wilt suppressive soil and a fusarium wilt conducive soil from the Châteaurenard region. Front. Microbiol.9, 568. doi: 10.3389/fmicb.2018.00568
33
SinghN.SiddiquiZ. A. (2015). Effects of Bacillus subtilis, Pseudomonas fluorescens and Aspergillus awamori on the wilt-leaf spot disease complex of tomato. Phytoparasitica43, 61–75. doi: 10.1007/s12600-014-0427-0
34
SoesantoL.MugiastutiE.KhoerurizaK. (2019). Granular formulation test of Pseudomonas fluorescens P60 for controling bacterial wilt (Ralstonia solanacearum) of tomato in planta. J. Agr. Sci.41, 513–523. doi: 10.17503/agrivita.v41i3.2318
35
SotoyamaK.AkutsuK.NakajimaM. (2017). Suppression of bacterial wilt of tomato by soil amendment with mushroom compost containing Bacillus amyloliquefaciens IUMC7. J. Gen. Plant Pathol.83, 51–55. doi: 10.1007/s10327-016-0690-7
36
SureshP.VarathrajuG.ShanmugaiahV.AlmaaryK. S.ElbadawiY. B.MubarakA. (2021). Partial purification and characterization of 2, 4-diacetylphloroglucinol producing Pseudomonas fluorescens VSMKU3054 against bacterial wilt disease of tomato. Saudi J. Biol. Sci.28, 2155–2167. doi: 10.1016/j.sjbs.2021.02.073
37
TangC.SunB.ZeeshanM.LiJ.ZhangX. (2023). Funneliformis mosseae–induced changes of rhizosphere microbial community structure enhance Capsicum annuum L. plant growth and fruit yield. Soil Sci. Soc Am. J.87, 843–855. doi: 10.1002/saj2.20528
38
TreutterD. (2005). Significance of flavonoids in plant resistance and enhancement of their biosynthesis. Plant Biol.7, 581–591. doi: 10.1055/s-2005-873009
39
VargasC.JiménezH. R.AlmarioC. G.AlmarioA. G. (2022). Effect of lactic acid bacteria on the control of Fusarium oxysporum and Ralstonia solanacearum on singly infected and co-infected tomato plants. J. Appl. Microbiol.134, 1–9. doi: 10.1093/jambio/lxac053
40
Villegas-EscobarV.González-JaramilloL. M.RamírezM.MoncadaR. N.Sierra-ZapataL.OrduzS.et al. (2018). Lipopeptides from Bacillus sp. EA-CB0959: Active metabolites responsible for in vitro and in vivo control of Ralstonia solanacearum. Biol. Control125, 20–28. doi: 10.1016/j.biocontrol.2018.06.005
41
WangJ.PengY.XieS.YuX.BianC.WuH.et al. (2023). Biocontrol and molecular characterization of Bacillus velezensis D against tobacco bacterial wilt. Phytopathol. Res.5, 50. doi: 10.1186/s42483-023-00204-x
42
WangW.Portal-GonzalezN.WangX.LiJ.LiH.PortielesR.et al. (2024). Metabolome-driven microbiome assembly determining the health of ginger crop (Zingiber officinale L. Roscoe) against rhizome rot. Microbiome12, 167. doi: 10.1186/s40168-024-01885-y
43
WeiC.LiangJ.WangR.ChiL.WangW.TanJ.et al. (2024a). Response of bacterial community metabolites to bacterial wilt caused by Ralstonia solanacearum: a multi-omics analysis. Front. Plant Sci.14, 1339478. doi: 10.3389/fpls.2023.1339478
44
WeiX.LongY.ChenY.MatakaS. N.JiangX.ZhouY.et al. (2024b). Cladophialophora guangxiense sp. nov., a new species of dark septate endophyte, mitigates tomato bacterial wilt and growth promotion activities. Agronomy14, 2771. doi: 10.3390/agronomy14122771
45
WuX.SunZ.QiF.LiuH.ZhaoM.WangJ.et al. (2023). Cytological and transcriptomic analysis to unveil the mechanism of web blotch resistance in Peanut. BMC Plant Biol.23, 518. doi: 10.1186/s12870-023-04545-9
46
XiaoS.MingY.HuQ.YeZ.SiH.LiuS.et al. (2023). GhWRKY41 forms a positive feedback regulation loop and increases cotton defence response against Verticillium dahliae by regulating phenylpropanoid metabolism. Plant Biotechnol. J.21, 961–978. doi: 10.1111/pbi.14008
47
XueL.GuM. Y.XuW. L.LuJ. J.XueQ. H. (2016). Antagonistic Streptomyces enhances defense-related responses in cotton for biocontrol of wilt caused by phytotoxin of Verticillium dahliae. Phytoparasitica44, 225–237. doi: 10.1007/s12600-016-0517-2
48
YangL.WeiZ.VallsM.DingW. (2022). Metabolic profiling of resistant and susceptible tobaccos response incited by Ralstonia pseudosolanacearum causing bacterial wilt. Front. Plant Sci.12, 780429. doi: 10.3389/fpls.2021.780429
49
YangB.ZhengM.DongW.XuP.ZhengY.YangW.et al. (2023). Plant disease resistance-related pathways recruit beneficial bacteria by remodeling root exudates upon Bacillus cereus AR156 treatment. Microbiol. Spectr.11, e0361122. doi: 10.1128/spectrum.03611-22
50
ZhangX.HuangT.WangQ.GuoY.ZhangP.XieH.et al. (2022). Mechanisms of resistance to spot blotch in Yunnan iron shell wheat based on metabolome and transcriptomics. Int. J. Mol. Sci.23, 5184. doi: 10.3390/ijms23095184
51
ZhengY.WangJ.ZhaoW.CaiX.XuY.ChenX.et al. (2022). Effect of bacterial wilt on fungal community composition in rhizosphere soil of tobaccos in tropical Yunnan. Plant Pathol. J.38, 203–211. doi: 10.5423/PPJ.OA.03.2022.0035
52
ZhouD.JingT.ChenY.WangF.QiD.FengR.et al. (2019). Deciphering microbial diversity associated with Fusarium wilt-diseased and disease-free banana rhizosphere soil. BMC Microbiol.19, 161. doi: 10.1186/s12866-019-1531-6
53
ZhouJ. Y.LiX.ZhaoD.Deng-WangM. Y.DaiC. C. (2016). Reactive oxygen species and hormone signaling cascades in endophytic bacterium induced essential oil accumulation in Atractylodes lancea. Planta244, 699–712. doi: 10.1007/s00425-016-2536-0
54
ZhuY.LiM.LiK.LiJ.LiuX.YangS.et al. (2024). Integrated metabolome and transcriptome analyses reveal the critical role of alpha-linolenic acid metabolism in Panax notoginseng root rot disease. Ind. Crop Prod.222, 120072. doi: 10.1016/j.indcrop.2024.120072
Summary
Keywords
bacterial wilt (Ralstonia solanacearum), biocontrol agent, transcriptomics, metabolomics, bacterial community
Citation
Du X-Q, Sun T-X, Xu W-L, Zhu T, Wang Q, Gu P-W and Lu J (2025) Multi-omics analysis reveals the specific role of biocontrol reagents against tomato bacterial wilt. Front. Plant Sci. 16:1620460. doi: 10.3389/fpls.2025.1620460
Received
29 April 2025
Accepted
23 June 2025
Published
14 July 2025
Volume
16 - 2025
Edited by
Charith Raj Adkar-Purushothama, Université de Sherbrooke, Canada
Reviewed by
Helena Gil Azinheira, University of Lisbon, Portugal
Dang Fengfeng, Yan’an University, China
Updates
Copyright
© 2025 Du, Sun, Xu, Zhu, Wang, Gu and Lu.
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: Pei-Wen Gu, gupeiwen2019@nxu.edu.cn; Jiang Lu, jiang.lu@sjtu.edu.cn
Disclaimer
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.