Skip to main content

High abundance of Ralstonia solanacearum changed tomato rhizosphere microbiome and metabolome



Rhizosphere microbiome is dynamic and influenced by environment factors surrounded including pathogen invasion. We studied the effects of Ralstonia solanacearum pathogen abundance on rhizosphere microbiome and metabolome by using high throughput sequencing and GC-MS technology.


There is significant difference between two rhizosphere bacterial communities of higher or lower pathogen abundance, and this difference of microbiomes was significant even ignoring the existence of pathogen. Higher pathogen abundance decreased the alpha diversity of rhizosphere bacterial community as well as connections in co-occurrence networks. Several bacterial groups such as Bacillus and Chitinophaga were negatively related to the pathogen abundance. The GC-MS analysis revealed significantly different metabolomes in two groups of rhizosphere soils, i.e., the rhizosphere soil of lower harbored more sugars such as fructose, sucrose and melibiose than that in high pathogen abundance.


The dissimilar metabolomes in two rhizosphere soils likely explained the difference of bacterial communities with Mantel test. Bacillus and Chitinophaga as well as sugar compounds negatively correlated with high abundance of pathogen indicated their potential biocontrol ability.


Microbiota associated with the host plant, especially roots, determine the infection of the soil-borne pathogen [29]. Assembling a self-serving rhizosphere microbiota is vital for both plant and pathogen. Plants recruit beneficial microbes to stimulate plant growth, elicit plant systemic defense, and antagonize pathogens [3, 32]. This kind of recruitment is realized by the release of specific compound, and by alteration of hormone expression levels in plant [4, 35]. For plant pathogens, recognition of host plant via root exudates is the first step of successful invasion [7, 12]; in addition, suppression of host defenses and acquisition of essential nutrients also master the invasion events [7, 15, 44]. The interactions between pathogens or between pathogen and other microbes may promote pathogenicity by substance exchange or effectors gene expression enhancement [16, 31]. Such studies indicated that there is a pathogen helper during the process of pathogenesis. Investigation of rhizosphere microbiome after pathogen invasion may give insights on interactions between the host plant and pathogen as well as relationship between pathogen and other microbes.

Soil microbial community is causally affected by abiotic factors such as physiochemical soil properties, pH, temperature, and soil moisture [13, 25]. Whereas, little attention is given to the effects of soil metabolome on soil microbiota, such as sugars, organic acids, and other low molecular weight compounds, which can serve as available nutrients for soil microbes and can be responsible for the interactions of plants with rhizosphere microbes [26]. Root exudates are the main source of rhizosphere soil low molecular weight molecules, and changes in the root exudates composition can affect the composition of rhizosphere soil microbial community [17, 26, 48]. The soil nutrient condition and health status of host plant can also shift the composition of rhizosphere microbiota by affecting root exudates [30].

Species invasion has a strong impact on the original ecosystem through competition of resources and niches with the indigenous microbes [43]. In rhizosphere soil, pathogen invasion can also affect the microbial community by indirectly influencing the host [38, 41], for example, Berendsen et al [4] showed that Arabidopsis thaliana specifically promotes three bacterial species in the rhizosphere upon foliar defense activation by the downy mildew pathogen Hyaloperonospora arabidopsidis. As a gram-negative phytopathogenic bacterium, Ralstonia solanacearum leads to the global bacterial wilt disease [51,52,53]. It showed the rapid capacity to spread to the aboveground after invasion of xylem vessels [54,55,56] and then excessively produce the extracellular polysaccharides to block the water flow [57, 58].

In this study, the rhizosphere soil samples were collected from a 7-year-monoculture greenhouse characterized by a high abundance of R. solanacearum with 100% disease incidence at the last of fruiting stage previously. The tomato plants together with the soils were sampled at the late developmental stage. Total pathogen density was detected by q-PCR, low molecular weight molecules of soil were determined by GC-MS technology, and high throughput sequencing was used to determine changes in the composition of the rhizosphere bacterial community. We set up this study to answer two questions: (1) does R. solanacearum invasion determine the composition of rhizosphere microbial community? (2) do low molecular weight molecules of the rhizosphere soil play a role in the changes of the composition of rhizosphere microbial community?


The abundance of Ralstonia solanacearum in tomato rhizosphere

We selected 12 soil samples and divided them into higher RS abundance group (HRS) and lower RS abundance group (LRS), thus having 6 replicates of each group. The abundance of RS was significantly different (p < 0.05, t-test) between two groups (Fig. 1). Then 12 samples were sequenced with 16S rRNA amplicon by Illumina Hiseq technology to unfold the rhizosphere bacterial microbial community.

Fig. 1

The abundance of RS in two rhizosphere soil samples measured by quantitative PCR

Bacterial community composition

The average read count of all samples was 15,143 (with the standard deviation (SD): 3416). The number of the OTUs could reflect the diversity of the bacterial communities and numbers of the OTUs range between 213 and 410 per sample with an average of 340 (SD: 59). The top four phylum with the highest relative abundance were Proteobacteria (76.72%), Bacteroidetes (13.63%), Cyanobacteria (4.36%), and Actinobacteria (1.98%) (Fig. 2a). There was a clear difference between two community composition evaluated with principal coordinates analysis (PCoA) and the difference of community composition were statistically significant through analysis of similarity (PERMANOVA) (p = 0.001, R2 = 0.44) (Fig. 2b). Upon closer inspection of these bacterial communities, one OTU belonging to Ralstonia was dominant among Proteobacteria (Fig. 2a). The sequencing data upon Ralstonia showed the similar result with the quantitative data that HRS group harbored a significant higher relative abundance than LRS groups (Fig. 2a). After the exclusion of Ralstonia, the differences among the two groups was still statistically significant (p = 0.001 in PERMANOVA, Supplementary Fig. 1).

Fig. 2

a The maximum-likelihood phylogenetic tree of all OTUs clustered. The outer ring shows the taxonomic group, the central ring shows the high abundance of RS group, and the inner ring shows the low abundance of RS group. The relative abundance of RS in two groups were present at the left bottom, and asterisk indicated the significant difference between two groups. b Principal coordinates analysis (PCoA) with Bray-Curtis dissimilarity of the rhizosphere bacterial community in low abundance of RS samples (LRS) and high abundance of RS samples (HRS)

The alpha diversity indicators were calculated considering all the observed OTUs, the high abundance of Ralstonia decreased the alpha diversity of rhizosphere bacterial communities of HRS compared to LRS group (Fig. 3).

Fig. 3

The alpha diversity indices of rhizosphere bacterial community in low abundance of RS (LRS) samples and high abundance of RS (HRS) samples. Shannon (a), Chao (b), and Faith’s PD (c). Asterisk indicated the significant difference between two groups

Comparing changes in the rhizosphere bacterial co-occurrence networks

To gain further insight into the effect of pathogen invasion on bacterial community, we constructed bacterial co-occurrence networks based on the correlation analysis of taxonomic profiles in high and low Ralstonia abundance rhizosphere of tomato (Fig. 4). The networks differed between two groups. In general, the higher abundance of Ralstonia led to a clear decrease in the number of connections. There were 197 and 310 nodes with 666 and 2153 links for HRS and LRS groups, respectively. The other global statistics of this network analysis are listed in Supplementary Table 1. Dominate hubs in two networks showed some differences, in detail, the major hub bacteria were Paenibacillus, Flavisolibacter, Chitinophagaceae, Sphingobacteriales, MND1, Solibacterales, Rhodospirillaceae and GOUTA19 in the HRS rhizosphere, and Xiphinematobacter, Xanthomonadaceae, iii1–15, H39, EW055, MND1 and DS-18 in the LRS rhizosphere. The Ralstonia was connected with only one OTU classified as Kaistobacter, which may be potentially important bacterial group related to Ralstonia invasion.

Fig. 4

Co-occurrence network analysis of the rhizosphere bacterial community in high abundance of RS (HRS) samples (a) and low abundance of RS (LRS) samples (b). Co-occurrence networks of tomato rhizobacteria were constructed based on correlation analysis of taxonomic profiles. Connections are drawn between nodes that were significantly (P < 0.01; Spearman’s rank correlation test) and highly (Spearman’s r > 0.96) correlated. Hub bacteria for each network are ranked according to the number of connections in the network. For (A) 1): O_Solibacterales; 2): G_Flavisolibacter; 3): G_Paenibacillus; 4): F_Chitinophagaceae; 5): O_Sphingobacteriales; 6): O_MND1; 7): O_Solibacterales; 8): O_Sphingobacteriales; 9): G_GOUTA19; 10): F_Rhodospirillaceae. (B) 1): O_DS-18; 2): O_iii1–15; 3): Unassigned; 4): G_Candidatus Xiphinematobacter; 5): O_H39; 6): O_EW055; 7): O_MND1; 8): F_Xanthomonadaceae

Potential bacterial taxonomies related to pathogen resistance

To investigate the potential bacterial groups involved in pathogen multiplication, we further screened the relative abundance at the genus level. Ten genera appeared with their relative abundance negatively related to the pathogen abundance, and their relative abundance were all significantly higher in the LRS group than in HRS group (Fig. 5). For example, Chitinophaga was significantly (p < 0.05, t-test) higher in LRS group than in HRS group with the relative abundance of 17.95 and 4.6%, respectively.

Fig. 5

The relationship between targeted bacterial groups and RS in tomato rhizosphere bacterial microbiome. x-axis means the relative abundance of targeted bacterial groups (%), and y-axis means the relative abundance of RS (%)

Characterization of rhizosphere soil metabolites profile and relationship with the composition of bacterial community

In order to establish a mechanistic explanation about the effects of pathogen abundance on the rhizosphere soil microbiome, rhizosphere soil metabolites of high and low RS abundance were collected and analyzed by gas chromatography-time of flight mass spectrometry (GC-TOF-MS). Principal component analysis (PCA) based on all of the 694 detected peaks revealed the clearly different patterns of rhizosphere soil metabolites between LRS and HRS treatment (AMOVA, p < 0.05, Fig. 6a). According to the structure of the molecules, 249 of the 694 detected peaks were divided into 8 categories, including sugar alcohols (8), sugars (38 compounds), sugar acid (4), low molecular weight organic acids (31), long-chain organic acids (37), esters (13), amides (40), amino acids, alcohols (22), or others (56). Differential analysis indicated that 16 compounds of them showed significant differences (t-test, p < 0.05) between the two groups (Fig. 6b). When evaluated them upon the above classification level, sugars and small molecule organic acids were higher in LRS treatment, but only sugars showed a significant difference (t-test, p < 0.05) (Fig. 6c). Overall, rhizosphere soil from LRS group harbored a significantly more abundance of sugars, especially fructose, sucrose and melibiose.

Fig. 6

a Principal component analysis of rhizosphere soil small molecular profiles. b Heatmap analysis of changes in abundance of representative compounds that were significantly different in abundance between two groups. c Abundance (cumulative peak area) of compound categories. Each bar represents the average of three replicates. Asterisk indicated the significant difference between two groups

Mantel tests showed a positive and significant correlation between rhizosphere soil metabolites and bacterial phylogeny (Person: r = 0.44, p = 0.006). Partial CCA-based variation partitioning analysis (VPA) demonstrated that rhizosphere soil metabolites explained 17.2% of variance.


Here, we report the variations of the tomato rhizosphere microbiome in association with high and low abundance of pathogen R. solanacearum. Results indicated that pathogen content in the rhizosphere significantly impacted the composition and diversity of bacterial community. The relative abundance of Bacteroidetes was significantly dominated by the enrichment of RS within Proteobacteria. Similar trend occurred when the tomato rhizosphere was infected by P. parasitica [21]. Generally, the low microbial diversity in the rhizosphere soil is beneficial for the pathogen invasion [24]. Indeed, in our study the success of pathogen invasion decreased the microbial diversity as shown by the significant decrease in Shannon, Chao and Faith’s PD indexes at high RS abundance samples, which indicated that pathogen invasion disrupted the host associated microbiome [43]. Principal coordinates analysis showed that the alteration of pathogen abundance altered the composition of bacterial community, and analysis after removing OTU of RS further supported this conclusion. The abundance of several beneficial genera (Pseudomonas, Burkholderia, and Bacillus) was negatively correlated with that of pathogen. The species of these beneficial genera are reported to play important roles in controlling soil-borne disease in disease-suppressive soils by producing antimicrobial secondary metabolites and volatile organic compounds [6, 26, 34, 45].

The correlation analysis of the abundances of beneficial microbes and pathogen has been used to recognize potential anti-pathogen microbes in some publications [9, 28]. However, we cannot consider them as the beneficial bacteria only by simple correlation analysis for almost all of the taxa were relatively decreased due to the higher abundance of the pathogen. Further investigations are required to design effective screening for isolation of anti- Ralstonia microbes.

Interactions among rhizosphere microbes are important not only for the stability of the community but also for the host development. Generally, more connections within the community indicate its stability with capacity of restraining pathogen invasion [39, 43]. In this study, the high RS abundance was found to destroy the bacterial connections of tomato rhizosphere. R. solanacearum with the highest relative abundance showed the least connections with other microbes. Only one OTU classified as Kaistobacter was connected to RS in both HRS and LRS samples. Kaistobacter has been reported to be more prevalent in suppressive soils [18, 20, 23], and it is a potential PGPR group because it can produce hormone [19]. Co-occurrence network analysis indicated that the RS invasion affected rhizosphere microbe networks, furthermore, the Flavisolibacter with negatively correlated to RS is one important hub bacteria in HRS networks. Both Kaistobacter and Flavisolibacter can biologically control the bacterial wilt disease caused by R. solanacearum, but this hypothesis needs further verification.

The rhizosphere microbiome which showed a direct effect on plant health and resistance to pathogens [27, 33] was mastered by rhizosphere soil metabolites [26, 47]. Root exudates depend on many abiotic and biotic factors, with the former including pathogen invasion, especially colonization. For example, Fusarium oxysporum colonization could stimulate the secretion of phenolic acids [46]. Foliar pathogen infection increased the concentration of some amino acids and fatty acids in root exudates [47]. While a little knowledge of rhizosphere soil metabolites impacted by pathogen especially RS was revealed till now. Our study showed that the higher pathogen abundance in tomato rhizosphere significantly altered the soil metabolite profile compared to the lower pathogen abundance samples. Sugars abundant in low RS samples indicated that high sugars were not beneficial for R. solanacearum growth. It was likely because sugars, such as fructose and sucrose, were not carbon source for R. solanacearum [44]. Indeed, fructose and sucrose were abundant in low RS rhizosphere samples. Amino acids were associated to R. solanacearum invasion [15], because they had a higher relative concentration in HRS than LRS rhizosphere soil samples. Several other compounds showed significantly positive or negative correlation with pathogen abundance but further studies should verify their effects on pathogen colonization.


Interactions between host plants and pathogen in vivo are research hotspots in the field of plant protection. While the interactions happened in rhizosphere become attractive in plant and soil studies. Pathogen invasion even colonization can alter the composition of rhizosphere soil metabolites as well as microbial community, while in turn, the altered metabolites and microbial community can affect the colonization of pathogen [40, 43]. This kind of dynamic cause-and-effect interaction results in a complex system. This study revealed the features of rhizosphere microbiome and metabolites, with the pathogen invasion reducing the diversity and changing relative abundance of some groups. The high abundance of pathogen invasion also reduced the concentration of many sugars in the rhizosphere soil. However, our study did not clarify how microbial groups or metabolites were changed, by pathogen invasion or by host plant? Further studies are needed to better understand the interactions in rhizosphere between pathogen and host plants.


Experiment description and soil sampling

The greenhouse experiment of tomato (“Shi-Ji-Hong-Guan” (production ID: 1880518) is a commercial variety in China purchased from Jiuquan Sanbao Seed CO., LTD, which has not been deposited in any publicly available herbarium) was carried out at the Nanjing Institute of Vegetable Science (31°43′ N, 118°46′ E), Nanjing, Jiangsu province, China, which was a 7-year-monoculture long-term field for studying tomato bacterial wilt. For this study, tomato plants were planted from 22nd, March 2017 to 20th, June 2017, and the rhizosphere samples were collected on 7th June 2017 before the harvest of tomato fruit. In detail, based on visual observations of tomato bacterial wilt symptoms, the roots of healthy and weak tomatoes (9 replicates for each group) from the experiment site with 100% disease incidence previously. Then, all root samples were aseptically transferred to storage bags and maintained on ice prior to transport to the laboratory immediately. The rhizosphere soil samples were obtained according to Fu et al [14]. Briefly, soil loosely adhered to the tomato roots was shaken off and discarded, then the root tissues with their associated rhizosphere soil were cut into 1 cm segments by using a sterile scalpel under aseptic conditions. Sterile water was used to rinse the soil which closely bound to the tomato root segments. For collection of low molecular weight molecules, half of the suspension was taken out for lyophilization. The rest of them were centrifuged at 15,000 g for 20 min and the precipitations were used to extracted DNA.

Soil genome DNA extraction and Ralstonia solanacearum quantification

The rhizospheric soil samples were subsequently extracted the total DNA using the UltraClean Soil DNA Isolation Kits (MoBio Laboratories Inc., Carlsbad, USA) according to the manufacturer’s protocol. Two technical replicates per sample were used to minimize the DNA extraction bias, negative controls were also performed because of uncontrolled DNA polluted [36, 42]. Samples were stored at − 20 °C, and performing polymerase chain reaction was performed with the pooled DNA samples of technical replicates. The NanoDrop ND-2000 spectrophotometer (NanoDrop, ND2000, Thermo Scientific, 111 Wilmington, DE) was used to assess all DNA sample quality based on the 260/280 nm and 260/230 nm absorbance ratios. The concentration of extracted DNA ranged between 40 ng/μl and 60 ng/μl. Ralstonia solanacearum abundance was determined by quantitative PCR using primers targeting the fliC gene that coding the flagella subunit (forward primer: 5′-GAA CGC CAA CGG TGC GAA CT-3′ and reverse primer: 5′-GGC GGC CTT CAG GGA GGT C-3′) on Applied Biosystems StepOne Plus (Applied Biosystems, CA, USA) [37]. The standard curves generation and qPCR assay were performed as described in our previous publication [49]. Samples were separated into two groups based on the abundance of RS (Ralstonia solanacearum): high abundance of RS (HRS) and low abundance of RS (LRS).

PCR amplification and sequencing

Primers 338F/806R were used to amplified the V3-V4 part of the 16S rRNA gene of the extracted DNA [11]. The PCR amplification with the following amplification cycles: 95 °C for 5 min -- 30 cycles of 94 °C for 30 s + 52 °C for 30 s + 72 °C for 30 s -- 72 °C for 10 min. For amplification, the 25 μL reaction mixtures with 1 μL each primer (10 μM), 5 μL 5 x Q5 GC high enhancer, 2 μL deoxynucleoside triphosphates (dNTPs), 9.75 μL of sterilized ultrapure water, 5 μL 5 x PCR buffer, 0.25 μL 5 U/μl of Q5 polymerase and [8] and 1 μL DNA (20 ng/μL). The amplicon process was performed based on the reported best practices [36, 42]. The bands were excised and purified from 1.2% agarose gels using the MinElute PCR Purification Kit (Qiagen, Germany), and the QiagenQIAquick Gel Extraction kit (Qiagen, Germany). The PCR products were mixed and paired-end sequencing on the Illumina MiSeq sequencing platform at Guangdong Magigene Biotechnology Co., Ltd. China.

Amplicon sequence processing and analysis

The 16S rRNA gene sequences were processed using QIIME1.9.1 [5], USEARCH 10.0 and in-house scripts. Paired-end Illumina reads were filtered by FastQC and joined by script. Based on the high-confidence 16S representative sequences, an OTU table was generated by USEARCH (−usearch_global and scripts, cutoff = 0.97). Another filter step was performed to remove non-bacterial 16S rRNA gene sequences by aligning representative sequences of all OTUs to the Greengenes_13.5 database using PyNAST ( script). The taxonomy of the OTUs was classified with the Greengene_13.5 classifier.

To describe the composition of the rhizosphere community, 9519 sequences was extracted randomly for each sample to calculate alpha diversity index, included Shannon, Chao, and Faith’s PD indices, which estimated the script by Qiime. A phylogenetic tree was constructed using the script (with default settings – FastTree). Before calculation of beta diversity, we first used the CSS method to standardize the OTU profiles by script [59]. and Bray-Curtis similarity matrices were prepared using the script. The Principal coordinate analysis (PCoA) plots were generated from Bray-Curtis similarity matrices created using ggplot2 package in R version 3.4.3. For better evaluation the effect of the RS on rhizosphere bacterial community, we re-draw the PCoA plots excluding the 16S rRNA gene sequence of RS. Furthermore, the phylogenetic molecular ecological network analysis (MENA) [10, 50] was used to explore interaction networks in the two treatments, For MENA analysis, the OTU table was filtered by removing OTUs which relative abundance below 0.001 before uploading to the pipeline. The OTUs which at least in three of all samples were retain, some OTU which have missing value were kept blank and the Pearson correlation coefficient was selected for calculate similarity of OTUs. The correlations of the correlation matrix were filtered with correlation values ≥0.96, p-value < 0.05. The global properties of the networks were calculated in MENA pipeline. Gephi (v. 0.92) was used to visualize the networks finally [1].

Detection and analysis of low molecular weight molecules by GC-MS

Three rhizosphere soil samples from high RS abundance group or low RS abundance group were randomly picked (six in total) for the metabolite analysis. For the extraction of low molecular weight molecules, soil samples were divided into two parts (0.2 g for each) in 2 mL EP tubes and 24 μL of Adonitol (1 mg/mL stock in dH2O) were added as an internal standard. One part of the soil samples were homogenized in ball mill for 4 min at 45 Hz in 0.5 mL methanol solution (Vmethanol: VH2O = 3:1), then ultrasound treated for 5 min (incubated in ice water) for 5 times. The supernatant (0.4 mL) was transferred into a fresh centrifuge tube after centrifuged for 15 min at 10000×g at 4 °C. The residue was extracted with 0.5 mL ethyl acetate using the method mentioned above: then 0.4 mL of ethyl acetate extraction was transferred to the methanol extraction. The other part of the soil sample was extracted by ethyl acetate followed by methanol using the method reported above. Overall, 1.6 mL (0.4 × 4) solution was obtained from one soil sample, and 40 μL were taken and pooled as QC sample, and 1.2 mL were transferred into a fresh 2 mL GC/MS glass vial and dried with nitrogen gas. Then 20 μL methoxyamination hydrochloride (20 mg mL− 1 in pyridine) were added to the dried sample and incubated for 30 min at 80 °Cbefore being treated with 30 μL of the BSTFA (bis (trimethylsilyl) trifluoroacetamide) regent (1% TMCS (Trimethylchlorosilane),v/v); finally the mixture was incubated for 1.5 h at 70 °C. The GC-TOF-MS analysis and the raw peaks analysis were performed as the same as reported by Li et al. [22]

For rhizosphere soil metabolome analyses, principal component analysis (PCA) [60] was used to visualize the composition of two samples using the R package vegan. The pathway analysis was displayed by MetaboAnalyst (

Statistical analysis

Differences of alpha-diversity indexes between groups were tested using a nonparametric t-test in R revision 3.4.3. The differences of bacterial community composition and root exudates between the two treatments were tested using PERMANOVA (Permutational multivariate analyses of variance, Adonis, transformed data by Bray-Curtis, permutation = 999) and AMOVA (analysis of molecular variance). To determine the percent change in root exudates, we used t-tests for all compounds of rhizosphere soil metabolome with relative abundances to measure the significant difference in these abundances between the two samples. The P-values were corrected by the Benjamini-Hochberg FDR procedure for multiple comparisons [2]. All plots were created using R except the enrichment analysis of metabolic pathways plot.

Availability of data and materials

Raw sequences data were deposited in the NCBI Sequence Read Archive (SRA) database under accession number SRP246782 (



Operational taxonomic units


Gas Chromatography - Mass Spectrometry


Ralstonia solanacearum


lower RS abundance group


higher RS abundance group


Principal coordinates analysis


  1. 1.

    Bastian M, Heymann S, Jacomy M (2009) Gephi: an open source software for exploring and manipulating networks. Third international AAAI conference on weblogs and social media.

    Google Scholar 

  2. 2.

    Benjamini Y, Hochberg YJ, Jot Rsss B. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B. 1995;57:289–300.

  3. 3.

    Berendsen RL, Pieterse CM, Bakker PA. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86.

    CAS  Article  Google Scholar 

  4. 4.

    Berendsen RL, Vismans G, Yu K, Song Y, de Jonge R, Burgman WP, Burmolle M, Herschend J, Bakker P, Pieterse CMJ. Disease-induced assemblage of a plant-beneficial bacterial consortium. ISME J. 2018;12:1496–507.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. mBio. 2010;7:335–6.

  6. 6.

    Carrion VJ, Cordovez V, Tyc O, Etalo DW, de Bruijn I, de Jager VCL, Medema MH, Eberl L, Raaijmakers JM. Involvement of Burkholderiaceae and sulfurous volatiles in disease-suppressive soils. ISME J. 2018;12:2307–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Chuberre C, Plancot B, Driouich A, Moore JP, Bardor M, Gugi B, Vicre M. Plant immunity is compartmentalized and specialized in roots. Front Plant Sci. 2018;9:1692.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Claesson MJ, O'Sullivan O, Wang Q, Nikkilä J, Marchesi JR, Smidt H, de Vos WM, Ross RP, O'Toole PW. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS One. 2009;4:e6669.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Cross ML. Microbes versus microbes: immune signals generated by probiotic lactobacilli and their role in protection against microbial pathogens. FEMS Immunology & Medical Microbiology. 2002;34:245–53.

    CAS  Article  Google Scholar 

  10. 10.

    Deng Y, Jiang Y-H, Yang Y, He Z, Luo F, Zhou JJBB. Mol Ecol Netw Analyses. 2012;13:113.

    Article  Google Scholar 

  11. 11.

    Dennis KL, Wang Y, Blatner NR, Wang S, Saadalla A, Trudeau E, Roers A, Weaver CT, Lee JJ, Gilbert JA, Chang EB, Khazaie K. Adenomatous polyps are driven by microbe-instigated focal inflammation and are controlled by IL-10–producing T cells. Cancer Res. 2013.

  12. 12.

    Finkel OM, Castrillo G, Herrera Paredes S, Salas Gonzalez I, Dangl JL. Understanding and exploiting plant beneficial microbes. Curr Opin Plant Biol. 2017;38:155–63.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Flores-Rentería D, Rincón A, Valladares F, Curiel Yuste J. Agricultural matrix affects differently the alpha and beta structural and functional diversity of soil microbial communities in a fragmented Mediterranean holm oak forest. Soil Biol Biochem. 2016;92:79–90.

    CAS  Article  Google Scholar 

  14. 14.

    Fu L, Ruan Y, Tao C, Li R, Shen Q. Continous application of bioorganic fertilizer induced resilient culturable bacteria community associated with banana Fusarium wilt suppression. Sci Rep. 2016;6:27731.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Gu Y, Wei Z, Wang X, Friman V-P, Huang J, Wang X, Mei X, Xu Y, Shen Q, Jousset A. Pathogen invasion indirectly changes the composition of soil microbiome via shifts in root exudation profile. Biol Fertil Soils. 2016;52:997–1005.

    CAS  Article  Google Scholar 

  16. 16.

    Hanemian M, Barlet X, Sorin C, Yadeta KA, Keller H, Favery B, Simon R, Thomma BP, Hartmann C, Crespi M, Marco Y, Tremousaygue D, Deslandes L. Arabidopsis CLAVATA1 and CLAVATA2 receptors contribute to Ralstonia solanacearum pathogenicity through a miR169-dependent pathway. New Phytologist. 2016;211:502–15.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Huang X-F, Chaparro JM, Reardon KF, Zhang R, Shen Q, Vivanco JM. Rhizosphere interactions: root exudates, microbes, and microbial communities. Botany. 2014;92:267–75.

    Article  Google Scholar 

  18. 18.

    Jiang J, Yu M, Hou R, Li L, Ren X, Jiao C, Yang L, Xu H. Changes in the soil microbial community are associated with the occurrence of Panax quinquefolius L. root rot diseases. 2019;438(1):143-56.

  19. 19.

    Kolton M, Graber ER, Tsehansky L, Elad Y, Cytryn E. Biochar-stimulated plant performance is strongly linked to microbial diversity and metabolic potential in the rhizosphere. New Phytol. 2017;213:1393–404.

    CAS  Article  Google Scholar 

  20. 20.

    Kyselková M, Kopecký J, Frapolli M, Défago G, Ságová-Marečková M, Grundmann GL, Moënne-Loccoz Y. Comparison of rhizobacterial community composition in soil suppressive or conducive to tobacco black root rot disease. ISME J. 2009;3:1127.

    Article  Google Scholar 

  21. 21.

    Larousse M, Rancurel C, Syska C, Palero F, Etienne C, Industri B, Nesme X, Bardin M, Galiana E. Tomato root microbiota and Phytophthora parasitica-associated disease. Microbiome. 2017;5:56.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Li M, Song J, Ma Q, Kong D, Zhou Y, Jiang X, Parales R, Ruan Z, Zhang Q. Insight into the characteristics and new mechanism of nicosulfuron biodegradation by a Pseudomonas sp. LAM1902. J Agric Food Chem. 2020;68:826–37.

    CAS  Article  Google Scholar 

  23. 23.

    Liu X, Zhang S, Jiang Q, Bai Y, Shen G, Li S, Ding W. Using community analysis to explore bacterial indicators for disease suppression of tobacco bacterial wilt. Sci Rep. 2016;6:36773.

    CAS  Article  Google Scholar 

  24. 24.

    Locey KJ, Lennon JT. Scaling laws predict global microbial diversity. Proc Natl Acad Sci. 2016;113:5970–5.

    CAS  Article  Google Scholar 

  25. 25.

    Ma J, Ibekwe AM, Yang CH, Crowley DE. Bacterial diversity and composition in major fresh produce growing soils affected by physiochemical properties and geographic locations. Sci Total Environ. 2016;563-564:199–209.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Massalha H, Korenblum E, Tholl D, Aharoni A. Small molecules below-ground: the role of specialized metabolites in the rhizosphere. Plant J. 2017;90:788–807.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Mendes LW, Mendes R, Raaijmakers JM, Tsai SM. Breeding for soil-borne pathogen resistance impacts active rhizosphere microbiome of common bean. ISME J. 2018;12:3038–42.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Mendes R, Garbeva P, Raaijmakers JM. The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol Rev. 2013;37:634–63.

    CAS  Article  Google Scholar 

  29. 29.

    Mendes R, Raaijmakers JM. Cross-kingdom similarities in microbiome functions. ISME J. 2015;9:1905.

    CAS  Article  Google Scholar 

  30. 30.

    Mommer L, Kirkegaard J, van Ruijven J. Root-root interactions: towards a Rhizosphere framework. Trends Plant Sci. 2016;21:209–17.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Peyraud R, Cottret L, Marmiesse L, Genin S. Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen. Nat Commun. 2018;9:418.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Philippot L, Raaijmakers JM, Lemanceau P, Van Der Putten WH. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol. 2013;11:789.

    CAS  Article  Google Scholar 

  33. 33.

    Pieterse CM, Zamioudis C, Berendsen RL, Weller DM, Van Wees SC, Bakker PA. Induced systemic resistance by beneficial microbes. Annu Rev Phytopathol. 2014;52:347–75.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Raza W, Ling N, Yang L, Huang Q, Shen Q. Response of tomato wilt pathogen Ralstonia solanacearum to the volatile organic compounds produced by a biocontrol strain Bacillus amyloliquefaciens SQR-9. Sci Rep. 2016;6:24856.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Rudrappa T, Czymmek KJ, Pare PW, Bais HP. Root-secreted malic acid recruits beneficial soil bacteria. Plant Physiol. 2008;148:1547–56.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Scholer A, Jacquiod S, Vestergaard G, Schulz S, Schloter M. Analysis of soil microbial communities based on amplicon sequencing of marker genes. Biology and Fertility of Soils. 2017;53(5):485-89.

  37. 37.

    Schönfeld J, Heuer H, Van Elsas J, Smalla K. Specific and sensitive detection of Ralstonia solanacearum in soil on the basis of PCR amplification of fliC fragments. Appl Environ Microbiol. 2003;69:7248–56.

    Article  Google Scholar 

  38. 38.

    Scharf BE, Hynes MF, Alexandre GM. Chemotaxis signaling systems in model beneficial plant-bacteria associations. Plant Mol Biol. 2016;90:549–59.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Shi S, Nuccio EE, Shi ZJ, He Z, Zhou J, Firestone MK. The interconnected rhizosphere: high network complexity dominates rhizosphere assemblages. Ecol Lett. 2016;19:926–36.

    Article  PubMed  Google Scholar 

  40. 40.

    Thomashow LS, Kwak YS, Weller DM. Root-associated microbes in sustainable agriculture: models, Metabolites and Mechanisms. Pest Manag Sci. 2019.

  41. 41.

    Venturi V, Keel C. Signaling in the Rhizosphere. Trends Plant Sci. 2016;21:187–98.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Vestergaard G, Schulz S, Scholer A, Schloter M. Making big data smart—how to use metagenomics to understand soil quality. Biology and Fertility of Soils. 2017;53(5):479-84.

  43. 43.

    Wei Z, Hu J, Gu Y, Yin S, Xu Y, Jousset A, Shen Q, Friman VP. Ralstonia solanacearum pathogen disrupts bacterial rhizosphere microbiome during an invasion. Soil Biol Biochem. 2018a;118:8–17.

    CAS  Article  Google Scholar 

  44. 44.

    Wei Z, Yang T, Friman VP, Xu Y, Shen Q, Jousset A. Trophic network architecture of root-associated bacterial communities determines pathogen invasion and plant health. Nat Commun. 2015;6:8413.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Yuan J, Ruan Y, Wang B, Zhang J, Waseem R, Huang Q, Shen Q. Plant growth-promoting rhizobacteria strain Bacillus amyloliquefaciens NJN-6-enriched bio-organic fertilizer suppressed Fusarium wilt and promoted the growth of banana plants. J Agric Food Chem. 2013;61:3774–80.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Yuan J, Wu Y, Zhao M, Wen T, Huang Q, Shen Q. Effect of phenolic acids from banana root exudates on root colonization and pathogen suppressive properties of Bacillus amyloliquefaciens NJN-6. Biol Control. 2018a;125:131–7.

    CAS  Article  Google Scholar 

  47. 47.

    Yuan J, Zhao J, Wen T, Zhao M, Li R, Goossens P, Huang Q, Bai Y, Vivanco JM, Kowalchuk GA, Berendsen RL, Shen Q. Root exudates drive the soil-borne legacy of aboveground pathogen infection. Microbiome. 2018b;6:156.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Zhalnina K, Louie KB, Hao Z, Mansoori N, da Rocha UN, Shi S, Cho H, Zhao H, Karaoz U, Northen TR, Firestone MK. Dynamic root exudate chemistry and microbial substrate preferences drive patterns in rhizosphere microbial community assembly. Nat Microbiol. 2018;3:470.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Zhao M, Yuan J, Shen Z, Dong M, Liu H, Wen T, Li R, Shen Q. Predominance of soil vs root effect in rhizosphere microbiota reassembly. FEMS Microbiol Ecol. 2019;95.

  50. 50.

    Zhou J, Deng Y, Luo F, He Z, Tu Q, XJM Z. Functional molecular ecological networks. 2010;1:e00169–10.

  51. 51.

    Hayward AC. Biology and Epidemiology of Bacterial Wilt Caused by Pseudomonas Solanacearum. Ann Rev Phytopathol. 1991;29(1):65-87.

    CAS  Article  Google Scholar 

  52. 52.

    Jiang G, Wei Z, Xu J, Chen H, Zhang Y, She X, Macho AP, Ding W, Liao B. Bacterial wilt in China: history, current status, and future perspectives. Front Plant Sci. 2017;8.

  53. 53.

    Yabuuchi E, Kosako Y, Yano I, Hotta H, Nishiuchi Y. Transfer of two burkholderia and an alcaligenes species to Ralstonia gen. Nov.: proposal of Ralstonia pickettii (ralston, palleroni and doudoroff 1973) comb. Nov., Ralstonia solanacearum (smith 1896) comb. Nov. and Ralstonia eutropha (davis 1969) comb. Nov. Microbiol Immunol. 1995;39:897–904.

    CAS  Article  Google Scholar 

  54. 54.

    Dalsing BL, Truchon AN, Gonzalez-Orta ET, Milling AS, Allen C. Ralstonia solanacearum uses inorganic nitrogen metabolism for virulence, ATP production, and detoxification in the oxygen-limited host xylem environment. Mbio 2015;6.

  55. 55.

    Huang Q, Allen C. Polygalacturonases are required for rapid colonization and full virulence of Ralstonia solanacearum on tomato plants. Physiol Mol Plant Pathol. 2000;57:77–83.

    CAS  Article  Google Scholar 

  56. 56.

    Saile E, McGarvey JA, Schell MA, Denny TP. Role of extracellular polysaccharide and endoglucanase in root invasion and colonization of tomato plants by Ralstonia solanacearum. Phytopathology 1997;87:1264–1271.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Denny TP, Baek SR. Genetic-evidence that extracellular polysaccharide is a virulence factor of pseudomonas-solanacearum. Mol Plant Microbe Interact 1991;4:198–206.

    CAS  Article  Google Scholar 

  58. 58.

    Genin S, Denny TP. Pathogenomics of the Ralstonia solanacearum species complex. Ann Rev Phytopathol. 2012;50:67–89.

    CAS  Article  Google Scholar 

  59. 59.

    Paulson J, Stine O, Bravo H, Pop M, et al. Differential abundance analysis for microbial marker-gene surveys. Nat Methods 2013;10:1200–1202.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Worley B, Powers R. Multivariate Analysis in Metabolomics. Curr Metabolomics. 2013;1(1):92–107.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references




This study was financially supported by Natural Science Foundation of Jiangsu Province (BK20170724), the Fundamental Research Funds for the Central Universities (KJQN202017), Natural Science Foundation of China (31902107), National Postdoctoral Program for Innovative Talents (BX201600075), and Special Fund for Agro-scientific Research in the Public Interest: integrated management technology of crop wilt disease (201503110). The funding bodies had no role but provided the financial support to this research, including experimental implementation, sampling and data analysis.

Author information




JY, TW and TL: performed all experiments; JY and QS designed the study, and wrote the majority of the manuscript; MZ and QH: provided critical comments on the study, and helped to write the paper.

Corresponding author

Correspondence to Jun Yuan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Principal coordinates analysis (PCoA) with Bray-Curtis dissimilarity of the rhizosphere bacterial communities using the whole OTU table excluding the OTUs belonging to Ralstonia solanacearum.

Additional file 2: Table S1.

The characters of co-occurrence networks of two groups with different abundance of RS.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wen, T., Zhao, M., Liu, T. et al. High abundance of Ralstonia solanacearum changed tomato rhizosphere microbiome and metabolome. BMC Plant Biol 20, 166 (2020).

Download citation


  • Rhizosphere soil microbiome
  • Rhizosphere metabolome
  • Ralstonia solanacearum
  • Pathogen abundance
  • Co-occurrence network