Control of Streptomyces alfalfae XY25T Over Clubroot Disease and Its Effect on Rhizosphere Microbial Community in Chinese Cabbage Field Trials

Clubroot caused by Plasmodiophora brassicae is one of the most destructive diseases in cruciferous crops. Streptomyces alfalfae XY25^T, a biological control agent, exhibited great ability to relieve clubroot disease, regulate rhizosphere bacterial and fungal communities in Chinese cabbage, and promote its growth in greenhouse. Therefore, field experiments were carried out to investigate the effects of S. alfalfae XY25^T on clubroot and rhizosphere microbial community in Chinese cabbage. Results showed that the control efficiency of clubroot by S. alfalfae XY25^T was 69.4%. Applying the agent can alleviate soil acidification; increase the contents of soil organic matter, available nitrogen, available phosphorus, and available potassium; and enhance activities of invertase, urease, catalase, and alkaline phosphatase. During Chinese cabbage growth, bacterial diversity decreased first and then increased, and fungal diversity decreased gradually after inoculation with S. alfalfae XY25^T. High-throughput sequencing analysis showed that the main bacterial phyla were Proteobacteria, Bacteroidetes, Acidobacteria, and Planctomycetes, and the major fungal phyla were Ascomycota and Basidiomycota in rhizosphere soil. The dominant bacterial genera were Flavobacterium, Candidatus, Pseudomonas, Stenotrophomonas, Sphingomonas, Flavisolibacter, and Gemmatimonbacteria with no significant difference in abundance, and the major fungal genera were Monographella, Aspergillus, Hypocreales, Chytridiaceae, Fusarium, Pleosporales, Agaricales, Mortierella, and Pleosporales. The significant differences were observed among Pleosporales, Basidiomycota, Colletotrichum, two strains attributed to Agaricales, and another two unidentified fungi by using S. alfalfae XY25^T. Moreover, quantitative real-time PCR results indicated that P. brassicae content was significantly decreased after the agent inoculation. In conclusion, S. alfalfae XY25^T can affect rhizosphere microbial communities; therefore, applying the agent is an effective approach to reduce the damage caused by clubroot.

Introduction

Cruciferous plants are widely planted worldwide and occupy an important position in agricultural production (Botero et al., 2019). Clubroot disease is a highly destructive soil-borne disease caused by the pathogen Plasmodiophora brassicae (Dixon, 2014). Cruciferous crops are highly susceptible and are greatly affected by clubroot diseases. The prevention and treatment of cruciferous clubroot disease have always become the focus of research. Rotation, applying chemical agents such as cyclophosphamide, and breeding resistant varieties can inhibit the spore vitality of P. brassicae and reduce the occurrence of clubroot disease (Peng et al., 2014; Yang et al., 2020). However, these methods can only alleviate crop clubroot disease, but not completely eliminate it; moreover, the residue of chemical agents will endanger food safety and human health (Lee et al., 2008; Kowata-Dresch and May-De Mio, 2012; Botero et al., 2019; Mehraj et al., 2020). In addition, P. brassicae has strong pathogenicity and rapid mutation, and the resistance of crop-resistant varieties is weak (Hirai, 2006; Strelkov et al., 2018). Therefore, these control methods have certain restrictions.

Biological control is an effective control method by applying microbial fertilizers or agents to change the composition of the root soil microbial community and soil acidity and alkalinity, thus affecting the absorption of water, soil organic matter, and other nutrients by crops, and promote the growth of crops (O’Brien, 2017; Rui et al., 2017; Liu et al., 2018). Streptomyces is the largest genus of Actinomycetes, and the active ingredients produced by this strain play a key role in maintaining the balance of soil microorganisms; therefore, it was widely used as biological control agents (BCAs) for plant diseases (Minuto et al., 2005; Chen et al., 2016). S. alfalfae XY25^T is a BCA identified by She et al. (2016) and it was a Gram-positive, aerobic bacterium with extensive branched substrate mycelia and aerial. Our recent findings found that this strain could improve continuous cropping obstacles and promote crop growth, and it showed good application prospects in the prevention and treatment of clubroot diseases under greenhouse conditions (Liu et al., 2020).

Rhizosphere refers to the soil environment including the root system and its surrounding areas. It is a place where pathogens attack plants and establish parasitic relationships (Kumar et al., 2020). Like other habitats, the rhizosphere contains a large number of interacting microbial populations, and it plays an important role in maintaining plant health (Wei et al., 2019; Daval et al., 2020). Rhizosphere microorganisms prevent or resist the invasion of soil-borne diseases through antagonism, nutrient competition, parasitism, and group sensing (Mendes et al., 2013), and the imbalance of microbial community in rhizosphere soil is an important cause of soil-borne diseases. Improper fertilization, lack of soil nutrients, and accumulation of root exudates lead to changes in the soil microenvironment, resulting in selective adaptation of rhizosphere soil microorganisms, enriching certain specific microbial populations, breaking the soil microecological balance, and eventually causing plant disease (Sudini et al., 2011; Yang et al., 2020).

Cabbage, radish, and other cruciferous vegetables are mainly cultivated in the alpine region in western Hubei Province, which is the largest off-season vegetable production base in Hubei with a planting area of more than 145,400 ha in 2019 (Jiao et al., 2020). Vegetable cultivation is the pillar industry to increase local farmers’ income. However, in recent years, clubroot has spread all over the vegetable area, and the clubroot incidence in Chinese cabbage fields has reached up to more than 60%, which has become the main problem faced by the alpine vegetable industry development in western Hubei (Liu et al., 2020). Our previous pot experiments revealed some obvious inhibiting effects of S. alfalfae XY25^T on clubroots disease. However, there is still lack of evidence of S. alfalfae XY25^T to be used in fields. Therefore, Chinese cabbage field planting experiments in the clubroot high-risk area were conducted to evaluate the control effect of Streptomyces on clubroot disease, and quantitative real-time PCR and high-throughput sequencing analysis were performed to reveal the influence of this strain on rhizosphere microorganisms.

Materials and Methods

Strain Inoculation and Preparation

Streptomyces alfalfae XY25^T was obtained from State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan, China, and stored at the Chinese Type Culture Collection Center (CCTCCAA2015019T). The strain was inoculated in liquid Gauze’s synthetic medium at 28°C for 72 h on a rotary shaker (180 rpm). The cells were harvested immediately after centrifugation at 5,000 rpm for 10 min and diluted with sterile distilled water (SDW) for subsequent use.

Field Experiment and Sampling

The experiment was carried out from August to October 2018. The vegetable growing area of Miaoziling Village, Huoshaoping Township, Changyang County, Yichang City, Hubei Province, China (30.48N, 110.77E), was selected as the experimental site, where the altitude is 1,833 m with average annual temperature of 7.6°C and total rainfall of about 1,500 mm/year. In the selected field, severe clubroot disease had occurred, and it had lots of clubroot pathogens in the soil. The seedlings of Chinese cabbage at the four-leaf stage were transplanted into the field, which was divided into six plots. Three of the plots were randomly assigned to the treatment group (S), and the other three were assigned as the control group (C). An aliquot of 10 ml of S. alfalfae XY25^T (1 × 10⁶ CFU/ml) was inoculated into the root rhizosphere of each plant in the treatment group at 0, 14, and 28 days after transplantation. Meanwhile, the plants in the control group were treated with the same volume of SDW. Before each inoculation of BCA or SDW, the rhizosphere soil samples of plants were collected by root-shaking method (Bakker et al., 2015) at days 0, 14, 28, and 42, respectively. Two plants were randomly selected from each plot for investigation, and the soil samples in each group contained six biological replicates. Each sample was randomly divided into two parts. One part was immediately refrigerated for the determination of physicochemical index and enzyme activity, and the other part was immediately frozen and stored at −80°C for DNA extraction. At the end of the experiment (day 42), all the plants were pulled out, and their root disease incidence and disease index were evaluated according to the previous study (Liu et al., 2018).

Physicochemical Index and Enzyme Activity Determination

Soil pH was determined with a pH meter after shaking soil–water suspension (1:2.5) for 30 min. The soil samples were air dried, ground, and sieved (<2 mm). The contents of soil organic matter (SOM), available nitrogen (AN), available phosphorus (AP), and available potassium (AK) were measured by potassium dichromate titration, alkaline hydrolysis diffusion, molybdenum–antimony colorimetry, and flame photometer method, respectively (Yang et al., 2017). The 3,5-dinitrosalicylic acid colorimetry (Frankeberger and Johanson, 1983), indophenol blue colorimetry (Huang et al., 2012), disodium phosphate colorimetry, and potassium permanganate titration method (Johnson and Temple, 1964) were used to determine the activity of invertase, urease, alkaline phosphatase, and catalase, respectively.

DNA Extraction, High-Throughput Sequencing, and Data Analysis

Soil DNA was extracted from the soil samples with E.Z.N.A Soil DNA Kit (Omega, Norcross, Georgia, United States). The quality and quantity of the extracted DNA were checked by 1% agarose gel electrophoresis and a Nanodrop One spectrophotometer (Thermo, Waltham, MA, United States), respectively. The bacterial 16S rRNA was amplified by using the specific primers 341F (5′-CCTAYGGGRBGCASCAG-3′) and 806R (5′GGACTACNNGGGTATCTAAT-3′) targeting the variable V3–V4 regions. The specific primers ITS5-1737F (5′-GGAAGTAAAAGTCGTAACAAGG-3′) and ITS2-2043R (5′-GCTGCGTTCTTCATCGATGC-3′) were used for fungal library construction by targeting the ITS1 region. Each PCR mixture (50 μl) consisted of 3 μl of DNA template (approximately 60 ng), 25 μl of 2 × SYBR Premix Ex Taq (TaKaRa, Tokyo, Japan), 1 μl of each primer, and 20 μl nuclease-free water. PCR amplification was performed as follows: preincubation at 94°C for 5 min; 30 cycles of 94°C for 30 s, 52°C for 30 s, and 72°C for 30 s; and finally, an extension at 72°C for 10 min. The PCR-amplified products were analyzed by pair-end sequencing using the Illumna MiSeq platform.

The DNA library was constructed following the NEBNext^® Ultra^TM DNA Library Prep Kit for Illumina^® standard procedure (Schloss et al., 2009). The constructed amplicon library was sequenced by PE250 on the Illumina Hiseq2500 platform. The PE reads obtained by MiSeq sequencing were spliced according to the overlap between PE reads, and the sequences were optimized using the software FLASH and Trimmomatic.

Usearch software was used for operational taxonomic unit (OTU) cluster analysis (Edgar, 2010). Bioinformatics statistical analysis of OTUs with 99% sequence similarity was carried out; OTUs were normalized to the same number of reads (27,679 reads in 16S sequence and 59,470 in ITS sequence, the smallest read number of samples) in each sample. Each OTU represented a species (McDonald et al., 2012). Multiple diversity indexes of OTUs were analyzed using the software QIIME (Werner et al., 2012). Based on taxonomic information, a statistical analysis of the community composition was performed at each classification level, and the composition of each sample was counted at each classification level (Price et al., 2010). Bacterial 16S ribosomal DNA (rDNA) gene V3–V4 regions were compared against Greengenes database by using RDP Classifier. Fungal ITS1 regions were compared against Unite database by blast method. Species annotation and abundance analysis were performed using QIIME platform software and KRONA software (Catherine and Rob, 2005). The sequencing data of 16S rDNA have been submitted to the National Center for Biotechnology Information (NCBI) database with accession number PRJNA681733.

Quantitative Real-Time PCR

The abundance of different groups of microorganisms was determined by quantitative real-time PCR with the whole genomic DNA of microorganisms in each soil sample as a template. The primer sequences and sources are shown in Supplementary Table 1. Each PCR mixture (25 μl) consisted of 1 μl of DNA template (approximately 20 ng), 12.5 μl of 2 × SYBR Premix Ex Taq (Takara, Tokyo, Japan), 1 μl of each primer, and 9.5 μl dd H₂O. PCR amplification was performed using a Thermal Cycler PCR apparatus (Bio-Rad, Hercules, CA, United States) as follows: preincubation at 94°C for 5 min; 32 cycles of 94°C for 30 s, 54–66°C for 1 min, and 72°C for 1 min; and finally, an extension at 72°C for 10 min. Standard curves were, respectively, generated using serially diluted purified known genomic DNA standard by using the primers presented in Supplementary Table 1. The accuracy of the amplification results was tested by the melting curves.

Statistical Analysis

The data of soil physicochemical factors and qPCR were analyzed by SPSS (Version 26.0 and IBM Corp., Armonk, NY, United States). Tukey’s honestly significant difference (HSD) test was performed to reveal significant difference. Data were presented as mean ± standard deviation (SD). Origin Pro 2019b was used to draw the box plot of soil physicochemical factors.

Results

Control Efficiency

According to Table 1, at 42 days, the incidence of clubroot disease in the control group was 59.8% ± 4.3%, and the disease index was as high as 45.5% ± 2.8%, whereas in the treatment group, it was 21.8% ± 1.6% and 13.9% ± 1.78%, respectively. The control efficiency of S. alfalfae XY25^T on clubroot disease was 69.4%.

TABLE 1

Table 1. Statistics of disease rate, disease index, and control efficiency (day 42).

Physicochemical Factors

At the beginning of the growing season (day 0), there was no significant difference in physicochemical factors between the two groups (Figure 1). During the growth of Chinese cabbage, the soil pH of the control group continued to decrease, while that of the treatment group showed a slow upward trend. Statistical analysis results showed that the soil acidification trend was significantly alleviated by the application of S. alfalfae XY25^T (Figure 1A). The trends of other physicochemical factors in the treatment group and the control group were basically consistent (Figures 1B–I). The SOM, AN, AP, and AK contents of the treatment group were significantly higher than those of the control group (Figures 1B–E). Compared with the those in the control group, the contents of SOM, AN, AP, and AK in the treatment group at 42 days after transplantation increased by 30.9, 29.0, 11.3, and 50.5%, respectively (Figures 1B–E). The soil enzyme activities of the treatment group were also significantly higher than those of the control group. The soil activities of invertase, urease, alkaline phosphatase, and catalase significantly increased by 37.0, 36.9, 30.1, and 30.4%, respectively.

FIGURE 1

Figure 1. Effect of S. alfalfae XY25^T on physicochemical factors of soil. (A) pH value; (B) organic matter; (C) available nitrogen; (D) available phosphorus; (E) available potassium; (F) invertase; (G) urease; (H) alkaline phosphatase; (I) catalase. n = 6; different letters for the same item indicate significant differences at P < 0.05 level. The error bars represent the standard deviation (SD). The horizontal bars within boxes represent medians. The tops and bottoms of boxes represent the 75th and 25th percentiles, respectively. The upper and lower whiskers extend to data no more than 1.5 × the interquartile range from the upper edge and lower edge of the box, respectively.

Alpha Diversity

High-throughput sequencing was used to analyze the microbial alpha diversity of 24 soil samples (Table 2). The Chao1 index and Shannon index analysis of the bacterial community showed that there was no significant difference in alpha diversity between the two groups. The Shannon index of the control group increased first and then decreased. The change trend of the treatment group was complicated, and it exhibited a lower diversity than the control group. The Chao1 index and Shannon index analysis of the fungal community showed that there were no significant differences in these two indexes between the control groups, but significant differences were observed between the treatment groups. These results showed that the diversity of fungi in the rhizosphere soil was significantly reduced after the application of S. alfalfae XY25^T.

TABLE 2

Table 2. Analysis of alpha diversity estimators.

Microbial Community

During the whole growth of Chinese cabbage, the bacterial community was composed of 15 phyla, of which Proteobacteria was dominant, followed by Bacteroidetes, Acidobacteria, Planctomycetes, Thaumarchaeota, Chloroflexi, Gemmatimonadetes, Actinobacteria, Verrucomicrobia, Patescibacteria, etc. (Figure 2 and Table 3). A total of six phyla of fungal communities were detected at the phylum level. The most common one was Ascomycota, accounting for 45.7–72.9%, followed by Basidiomycota, Chytridiomycota, Zygomycota, Glomeromycota, Rozellomycota, and four unidentified phylum fungal communities (Figures 2, 4B). In the treatment group, the abundance of Basidiomycota decreased significantly with planting time (Table 4).

FIGURE 2

Figure 2. Effect of S. alfalfae XY25^T on microbial community. Bacteria at (A) phylum and (C) genus level and fungal at (B) phylum and (D) genus level. C0, C14, C28, and C42 represent the sampling days 0, 14, 28, and 42 in the control group, respectively; S0, S14, S28, and S42 represent the sampling days of the treatment group.

TABLE 3

Table 3. Main bacterial groups or the group with significant difference (%).

FIGURE 3

Figure 3. Principal components analysis (PCA) showed distribution of samples. (A) Bacterial community and (B) fungal community. The closer distance between points shows the higher similarity; n = 3.

FIGURE 4

Figure 4. Effect of S. alfalfae XY25^T on content of (A) total bacteria, (B) total fungi, (C) Plasmodiophora brassicae, and (D) Streptomyces. n = 3; the letter with * indicates significant differences at P < 0.05 level; the letters with ** indicate extremely significant differences at P < 0.01 level.

TABLE 4

Table 4. Main fungal groups or the group with significant difference (%).

At the genus level, the top 10 bacteria were Flavobacterium, Candidatus, Pseudomonas, Pyrinomonadaceae, Anaerolineaceae, Stenotrophomonas, Sphingomonas, Flavisolibacter, Nitrosomonadaceae, and Gemmatimonas (Figures 2, 4C). Statistical analysis showed that there was no significant difference in relative abundance among the top 10 genera. Some minor genera such as Acidovorax, Haliangium, Herbaspirillum, and Mizugakiibacter had significant difference in relative abundance. The abundance of Herbaspirillum in the treatment group increased significantly at 28 days. The top 20 genera of the fungi are shown in Figures 2, 4D. The top fungi included Monographella, Aspergillus, Hypocreales, Chytridiaceae, Fusarium, Pleosporales, Agaricales, Mortierella, Pleosporales, unidentified fungi, etc. (Figures 2, 4D). Analysis results showed that the most dominant genus, Monographella, had significant difference in relative abundance between the two groups. The abundance of Monographella in the treatment group was significantly higher than that in the control group at day 28. The significant differences were observed among Pleosporales, Basidiomycota, Colletotrichum, two strains attributed to Agaricales, and another two unidentified strains by using S. alfalfae XY25^T.

Principal Component Analysis

As shown in Figure 3A, for bacterial communities, the horizontal axis PCA1 accounted for 22.14%, and the vertical axis PCA2 occupied 14.82%. At the beginning (day 0), the bacterial community composition in the control group (C0) was similar to that in the treatment group (S0). With the growth of Chinese cabbage, the bacterial community composition in the rhizosphere soil changed, and the bacterial community compositions on days 14, 28, and 42 were far away from that at day 0. In the control group, the community compositions of the three replicates were similar. Compared with that in the control group, the bacterial community composition in the treatment group exhibited greater change, and the three replicates of treatment were farther apart (Figure 3A).

For the fungal community, the horizontal axis PCA1 accounted for 22.14%, and the vertical axis PCA2 occupied 16.93% (Figure 3B). The control group C0 and the treatment group S0 had similar fungal community composition. With the extension of growth time of Chinese cabbage at days 14, 28, and 42, there was a large difference between two groups. Compared with those in the control group, the difference among the three replicates in the treatment group was large (Figure 3B). These results indicated that the application of S. alfalfae XY25^T affected the microbial community in rhizosphere soil.

Correlations Between Microbial Community and Soil Edaphic Factors

As shown in Supplementary Figure 1A, the main bacterial communities that were significantly negatively correlated with soil pH were Ferruginibacter, Ellin6067, and UTCFX1. Flavobacterium was extremely positively correlated with AN, and this community was significantly positively correlated with urease, catalase, AP, and alkaline phosphatase. Pseudomonas and Chitinophaga were significantly positively correlated with SOM, AN, AP, AK, invertase, alkaline phosphatase, and urease. Stenotrophomonas was significantly positively correlated with SOM, AP, invertase, AN, and phosphatase. Flavisolibacter was significantly negatively correlated with urease and AN. Pedobacter was significantly positively correlated with urease, AN, and alkaline phosphatase. Other major communities such as Candidatus, Nitrosotalea, Sphingomonas, Gemmatimonas, Bryobacter, and Acidovorax exhibited no significant relationship with each factor.

As for fungi, soil pH was extremely significantly positively correlated with Junewangia, Rhizopus, and Phaeoacremonium; extremely negatively correlated with Leptosphaeria, Aphanoascus, and Robillarda; and significantly negatively correlated with Penicillium. Monographella was negatively correlated with SOM, AN, invertase, and alkaline phosphatase and negatively correlated with soil AP, AK, urease, and catalase. Aspergillus was extremely positively correlated with AK and positively correlated with alkaline phosphatase, AP, and SOM. Fusarium displayed a very significant positive correlation with alkaline phosphatase. Mortierella was significantly positively correlated with invertase and SOM and significantly positively correlated with AN and AK (Supplementary Figure 1B).

Microbial Population Determination Through Quantitative Real-Time PCR

During the growth of Chinese cabbage, the total bacterial content in the soil of the control group remained basically stable, whereas that of the treatment group increased first and then decreased, which was significantly higher on day 28 than on day 0 (Figure 4A). The total fungal content in the soil of the control group and the treatment group both increased first, and then, the control group reached the maximum value on day 28, followed by a decrease, whereas the treatment group reached the maximum value on day 42. However, no significant differences were observed between the control group and treatment group at each time point (Figure 4B). On day 14, the content of Plasmodiophora brassicae in the control group was significantly higher than that in the treatment group, indicating that there were signs of clubroot disease occurrence in Chinese cabbage on day 14 (Figure 4C). The content of Streptomyces in the control group remained stable, and its content in the treatment group gradually increased, which was significantly higher than that in the control group on day 42 after transplantation (Figure 4D).

Discussion

Chinese cabbage is often threatened by clubroot disease in its growth process. Generally, the clubroot disease occurs about at 30 days after the Chinese cabbage transplantation and eventually breaks out in an all-round way, resulting in crop yield reduction by more than 80% in severe cases (He et al., 2019). Biological methods have potential application prospects in the prevention and treatment of clubroot diseases. Lactobacillus plantarum PM411 and Lactobacillus plantarum TC92 can control the occurrence of plant diseases under relatively suitable humidity (Núria et al., 2018). One previous study has reported that Streptomyces platensis 3–10 inhibits the growth of spores of P. brassicae in vitro (Shakeel et al., 2016). The greenhouse experiment results have showed that the strains of endophytic actinomycetes can inhibit the occurrence of clubroot disease of Chinese cabbage (Lee et al., 2008; Liu et al., 2020). Based on the control efficiency of S. alfalfae XY25^T obtained in the greenhouse experiment, the field experiment was further carried out and showed that the clubroot disease incidence of Chinese cabbage was significantly reduced by applying S. alfalfae XY25^T, and the control effect was 69.4%. We hypothesized that the control efficiency of S. alfalfae XY25^T might be associated with changes in rhizosphere microbial community. Therefore, the rhizosphere samples of Chinese cabbage were collected to study the composition of microbial community by high-throughput sequencing and real-time PCR.

Microbial communities play an important role in soil structure and disease control of plants (Xue et al., 2018). In recent years, high-throughput sequencing technology has deepened the study of the diversity of soil microbial communities (Roesch et al., 2007; Blaalid et al., 2012; Liu et al., 2018). In this study, high-throughput sequencing and qPCR were used to study the microbial community during the occurrence of Chinese cabbage clubroot disease. The results showed that at the phylum level, Proteobacteria was the most dominant bacteria with an abundance of 23.9–40.6%, followed by Bacteroidetes, Acidobacteria, and Plantomytomytes. Acidobacteria has been reported to be one of the five largest taxa in the soil (Roesch et al., 2007; Liu et al., 2018), and previous studies have shown that Proteobacteria is the most dominant bacteria in the soil (Delgado-Baquerizo et al., 2018), followed by Basidiomycota (Chowdhury et al., 2013; Liu et al., 2018). Planctomycetes is also one of the main communities in this study, and this result is different from many other findings (Roesch et al., 2007). In this study, the largest bacterial group at the genus level was Flavobacterium, whose abundance in the treatment group was about 2% higher than that in the control group, but there was no significant difference between these two groups. Flavobacterium hercynium EPB-C313 has been found to have a significant inhibitory effect on P. brassicae (Hahm et al., 2012). A previous study has shown that Bacillus subtilis XF-1 can promote the growth of Flavobacterium and that Flavobacterium has a potential role in the prevention and treatment of clubroot disease (Liu et al., 2018). Our study indicated that the major microbial groups were “Candidatus Nitrosotalea,” Pseudomonas, Stenotrophomonas, Sphingomonas, Flavisolibacter, and Gemmatimonas in terms of abundance order. “Candidatus Nitrosotalea” is a representative of ammonia-oxidizing archaea, which is common in soil, and it plays an important role in the nitrogen cycle (Lehtovirta-Morley et al., 2016). Pseudomonas and Stenotrophomonas are common taxa in soil (Dai et al., 2020). Our data showed that after applying S. alfalfae XY25^T, the top 10 most abundant genera did not change significantly, but a few genera showed significant changes. For example, at 28 days after transplantation, the abundance of Herbaspirillum in the treatment group increased significantly. Herbaspirillum plays an important role in the circulation of nitrogen to promote the growth of sugarcane (Manoel da Silva et al., 2015).

In this study, the most dominant fungus at the phylum level was Ascomycota with its abundance ranging from 45.7 to 72.9%, which was consistent with a previous study finding that Ascomycota was the most abundant phylum in the soil, reaching 29.98–45.69% (Liu et al., 2018). Our data showed that the second dominant fungus was Basidiomycota. Interestingly, the abundance of Basidiomycota in the treatment group was 31.1% at day 0, but it was only 3.37% at day 42, indicating that the application of S. alfalfae XY25^T significantly decreased its Basidiomycota abundance as the crop grew. Our result was similar to a previous study report that Basidiomycota was the second-ranked phylum whose abundance was much lower than Ascomycota (Liu et al., 2018). However, some studies have reported different findings that the main phylum in rhizosphere soil was Basidiomycota (Buee et al., 2009; Blaalid et al., 2012). Our results indicated that at the genus level, the most abundant fungi were Monographella, Aspergillus, Hypocreales, Chytridiaceae, Fusarium, Pleosporales, Agaricales, Mortierella, and Pleosporales, which was not completely consistent with a previous study reporting that Fusarium, Cryptococcus, and Mortierella were the most abundant fungi (Liu et al., 2018). Our data indicated that the application effect of streptomyces on soil fungal communities is different from that on bacterial communities and that the abundance of many fungal communities changed significantly at different time points. For example, the abundance of Monographella was significantly higher in the treatment group than in the control group at 28 days. There were significant differences in abundance among other fungal communities belonging to Agaricales, Pleosporales, Basidiomycota, Colletotrichum, and two unknown fungi. There are many nonculturable microorganisms in the soil, of which these nonculturable fungi play an important role in the formation of clubroot disease. This study showed the application of S. alfalfae XY25^T had a great influence on the fungal community, which might be one of the reasons for preventing and controlling clubroot disease by this strain.

The imbalance of rhizosphere microbial community composition is an important cause of soil-borne diseases. Soil biodiversity and soil microbial composition determine ecosystem multifunctionality and soil bacterial diversity, which exert an inhibitory effect on the plant pathogen (van Elsas et al., 2012; Cameron et al., 2014). Generally, the richer the microbial community composition, the stronger the ability to inhibit pathogens (Raaijmakers et al., 2009). This study showed that during the growth of crops, soil bacterial diversity decreased first, and then increased, and fungal diversity significantly decreased after applying S. alfalfae XY25^T. PCA analysis results showed that the application of S. alfalfae XY25^T affected the microbial communities in the rhizosphere soil. Our qPCR results showed that the contents of total bacteria and total fungi remained unchanged or increased slightly. These results indicated that application of S. alfalfae XY25^T could cause some types of microorganisms to aggregate and grow. Moreover, qPCR results also showed that the content of P. brassicae in the treatment group was significantly higher than that in the control group on day 14 after transplantation. On day 28, clubroots were observed. The above results suggested that the application of S. alfalfae XY25^T reduced the incidence of clubroot disease of Chinese cabbage by inhibiting the growth of P. brassicae.

Soil microorganisms could effectively reflect soil fertility. Soil acidity, fertility, and enzyme activity are closely related to the activities of microorganisms, and these factors jointly affect the composition of microorganisms (Zhang et al., 2016; Rui et al., 2017; Xue et al., 2018). In general, the content of microorganisms is directly proportional to soil fertility. This study showed that the application of S. alfalfae XY25^T could significantly alleviate the trend of soil acidification and improve soil fertility and enzyme activity. Changes in soil enzyme activity have an important impact on the formation of soil organic matter and nutrient cycle (Liu et al., 2013). The application of microbial agents promotes the growth of microorganisms and enhances their activities so that more enzymes are produced, which is conducive to the production of organic matter.

Conclusion

In summary, S. alfalfae XY25^T could be used to reduce the clubroot incidence on Chinese cabbage in the outbreak area. After applying S. alfalfae XY25^T, the community of rhizosphere microorganisms in Chinese cabbage changed, and the fungal diversity decreased significantly. Furthermore, the application of S. alfalfae XY25^T also improved the soil fertility, reduced soil acidification, and inhibited the growth of P. brassicae. This study provides a better understanding of the importance of the microbial community in soil clubroot disease, and it also offers an effective strategy and theoretical basis for biological control of clubroot disease.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA681733.

Author Contributions

YL, SZ, and ZqZ designed the study. LQ and KL performed the experiments. YlH and YmH drafted the manuscript. ZjZ, XX, and SX analyzed the data. YlH and ZjZ contributed to the manuscript revision. ZqZ, SX, and SZ were the responsible for the visualization, supervision, project administration, and funding acquisition. YmH and YL contributed to the overall support of this study. All authors read and approved the final manuscript.

Funding

This work was supported by grant 41601257 from the NSFC, by grant 2020CFB527 of the NSF of Hubei Province, by grant 2018ABA092 and 2020BAB095 from the Hubei Major Project of Technique Innovation, by grant 2018skjcx01 from the Agricultural Science and Technology Innovation Project in Hubei Province, and by grant EWPL201709 from the Hubei Key Laboratory of Edible Wild Plants Conservation and Utilization.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.641556/full#supplementary-material

Supplementary Figure 1 | Correlations between soil physicochemical factors and genus of microbial community.

Supplementary Table 1 | Primers used for quantitative real-time PCR.

References

Bakker, M. G., Chaparro, J. M., Manter, D. K., and Vivanco, J. M. (2015). Impacts of bulk soil microbial community structure on rhizosphere microbiomes of Zea mays. Plant Soil 392, 115–126. doi: 10.1007/s11104-015-2446-0

CrossRef Full Text | Google Scholar

Blaalid, R., Carlsen, T., Kumar, S., Halvorsen, R., Ugland, K. I., Fontana, G., et al. (2012). Changes in the root-associated fungal communities along a primary succession gradient analysed by 454 pyrosequencing. Mol. Ecol. 21, 1897–1908. doi: 10.1111/j.1365-294X.2011.05214.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Botero, A., Garcia, C., Gossen, B. D., Strelkov, S. E., Todd, C. D., Bonham-Smith, P. C., et al. (2019). Clubroot disease in Latin America: distribution and management strategies. Plant Pathol. 68, 827–833. doi: 10.1111/ppa.13013

CrossRef Full Text | Google Scholar

Buee, M., Reich, M., Murat, C., Morin, E., Nilsson, R., Uroz, S., et al. (2009). 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol. 184, 449–456. doi: 10.1111/j.1469-8137.2009.03003.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameron, W., Franz, B. S., Franco, W., and Van der Heijden, M. G. A. (2014). Soil biodiversity and soil community composition determine ecosystem multifunctionality. Proc. Natl. Acad. Sci. U.S.A. 111, 5266–5270. doi: 10.1073/pnas.1320054111

PubMed Abstract | CrossRef Full Text | Google Scholar

Catherine, L., and Rob, K. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235. doi: 10.1128/aem.71.12.8228-8235.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. Y., Chen, P. C., and Tsay, T. T. (2016). The biocontrol efficacy and antibiotic activity of Streptomyces plicatus on the oomycete Phytophthora capsici. Biol. Control 98, 34–42. doi: 10.1016/j.biocontrol.2016.02.011

CrossRef Full Text | Google Scholar

Chowdhury, S. P., Dietel, K., Rändler, M., Schmid, M., Junge, H., Borriss, R., et al. (2013). Effects of Bacillus amyloliquefaciens FZB42 on lettuce growth and health under pathogen pressure and its impact on the rhizosphere bacterial community. PLoS One 8:e68818. doi: 10.1371/journal.pone.0068818

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Yang, F., Zhang, L., Xu, Z., Fan, X., Tian, Y., et al. (2020). Wheat−associated microbiota and their correlation with stripe rust reaction. J. Appl. Microbiol. 128, 544–555. doi: 10.1111/jam.14486

PubMed Abstract | CrossRef Full Text | Google Scholar

Daval, S., Gazengel, K., Belcour, A., Linglin, J., Guillerm-Erckelboudt, A.-Y., Sarniguet, A., et al. (2020). Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and brassica napus transcriptomes. Microb. Biotechnol. 13, 1648–1672. doi: 10.1111/1751-7915.13634

PubMed Abstract | CrossRef Full Text | Google Scholar

Delgado-Baquerizo, M., Oliverio, A. M., Brewer, T. E., Benavent-González, A., Eldridge, D. J., Bardgett, R. D., et al. (2018). A global atlas of the dominant bacteria found in soil. Science 359, 320–325. doi: 10.1126/science.aap9516

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon, G. R. (2014). Clubroot (Plasmodiophora brassicae Woronin)–an agricultural and biological challenge worldwide. Can. J. Plant. Pathol. 36, 5–18. doi: 10.1080/07060661.2013.875487

CrossRef Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

Frankeberger, W. T., and Johanson, J. B. (1983). Method of measuring invertase activity in soils. Plant Soil 74, 301–311. doi: 10.1007/BF02181348

CrossRef Full Text | Google Scholar

Hahm, S. S., Kim, J. T., Han, K. S., Kim, B. R., Kim, H. K., Nam, Y. K., et al. (2012). Biocontrol efficacy of endophytic bacteria Flavobacterium hercynim EPB-C313 for control of Chinese cabbage clubroot. Res. Plant Dis. 18, 210–216. doi: 10.5423/RPD.2012.18.3.210

CrossRef Full Text | Google Scholar

He, P., Cui, W., Munir, S., He, P., Li, X., Wu, Y., et al. (2019). Plasmodiophora brassicae root hair interaction and control by Bacillus subtilis XF-1 in Chinese cabbage. Biol. Control 128, 56–63. doi: 10.1016/j.biocontrol.2018.09.020

CrossRef Full Text | Google Scholar

Hirai, M. (2006). Genetic analysis of clubroot resistance in Brassica crops. Breed. Sci. 56, 223–229. doi: 10.1270/jsbbs.56.223

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Zhen, L. I., and Jian, Z. (2012). Improvement of indophenol blue colorimetric method on activity of urease in soil. J. Civ. Archit. Environ. Eng. 1, 102–107. doi: 10.11835/j.issn.1674-4764.2012.01.020

CrossRef Full Text | Google Scholar

Jiao, Z., Jiao, Z., Wu, J., Zhu, F., Guo, F., and Zhou, J., et al. (2020). Analysis of pest species and population dynamics of highland solanaceous vegetables in Hubei Province. China Veg. 375, 70–76.

Google Scholar

Johnson, J. L., and Temple, K. L. (1964). Some variables affecting the measurement of “Catalase activity” in soil. Soil Sci. Soc. Am. J. 28, 207–209. doi: 10.2136/sssaj1964.03615995002800020024x

CrossRef Full Text | Google Scholar

Kowata-Dresch, L., and May-De Mio, L. (2012). Clubroot management of highly infested soils. Crop Prot. 35, 47–52. doi: 10.1016/j.cropro.2011.12.012

CrossRef Full Text | Google Scholar

Kumar, P., Singh, J., Saini, A., and Singh, B. (2020). Antifungal activity exhibited bacterial strains isolated from wheat-rice rhizosphere soils and their antagonistic nature toward various fungal isolates. J. Pharmacogn. Phytochem. 9, 1521–1525.

Google Scholar

Lee, S. O., Choi, G. J., Choi, Y. H., Jang, K. S., Park, D. J., Kim, C. J., et al. (2008). Isolation and characterization of endophytic actinomycetes from Chinese cabbage roots as antagonists to Plasmodiophora brassicae. J. Microbiol. Biotechn. 18, 1741–1746. doi: 10.4014/jmb.0800.108

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehtovirta-Morley, L. E., Sayavedra-Soto, L. A., Gallois, N., Schouten, S., Stein, L. Y., Prosser, J. I., et al. (2016). Identifying potential mechanisms enabling acidophily in the ammonia-oxidizing archaeon “Candidatus Nitrosotalea devanaterra”. Appl. Environ. Microbiol. 82, 2608–2619. doi: 10.1128/aem.04031-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Yang, Z., He, P., Munir, S., Wu, Y., Ho, H., et al. (2018). Deciphering the bacterial and fungal communities in clubroot-affected cabbage rhizosphere treated with Bacillus subtilis XF-1. Agric. Ecosyst. Environ. 256, 12–22. doi: 10.1016/j.agee.2018.01.001

CrossRef Full Text | Google Scholar

Liu, K., Ji, Q., Jing, C., Liang, Y., Hu, Y., and Hu, Y. (2020). Effect of Streptomyces alfalfae XY25T application on physicochemical properties and microflora in clubroot-diseased soil. Microbiol. China 47, 97–108. doi: 10.13344/j.microbiol.china.190244

CrossRef Full Text | Google Scholar

Liu, L., Huang, B., Sun, J., Guo, S., Li, L., and Guo, H. (2013). Relationship between soil microbial quantity, enzyme activity and soil fertility in hot pepper greenhouse soils of different continuous cropping years. Soil Fert. Sci. China 2, 5–10.

Google Scholar

Manoel da Silva, J., Carvalho dos Santos, T. M., Santos de Albuquerque, L., Coentro Montaldo, Y., Ubaldo Lima de Oliveira, J., Mesquita da Silva, S. G., et al. (2015). Potential of the endophytic bacteria (Herbaspirillum spp. And Bacillus spp.) to promote sugarcane growth. Aust. J. Crop Sci. 9, 754–760.

Google Scholar

McDonald, R., Schreier, H. J., and Watts, J. E. (2012). Phylogenetic analysis of microbial communities in different regions of the gastrointestinal tract in Panaque nigrolineatus, a wood-eating fish. PLoS One 7:e48018. doi: 10.1371/journal.pone.0048018

PubMed Abstract | CrossRef Full Text | Google Scholar

Mehraj, H., Akter, A., Miyaji, N., Miyazaki, J., Shea, D. J., Fujimoto, R., et al. (2020). Genetics of clubroot and Fusarium wilt disease resistance in brassica vegetables: the application of marker assisted breeding for disease resistance. Plants 9:726. doi: 10.3390/plants9060726

PubMed Abstract | CrossRef Full Text | Google Scholar

Mendes, R., Garbeva, P., and Raaijmakers, J. M. (2013). The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663. doi: 10.1111/1574-6976.12028

PubMed Abstract | CrossRef Full Text | Google Scholar

Minuto, A., Spadaro, D., Garibaldi, A., and Gullino, M. L. (2005). Control of soilborne pathogens of tomato using a commercial formulation of Streptomyces griseoviridis and solarization. Crop Prot. 25, 601–613. doi: 10.1016/j.cropro.2004.11.003

CrossRef Full Text | Google Scholar

Núria, D., Esther, B., Jesús, F., Emilio, M., and Anna, B. (2018). Enhancing water stress tolerance improves fitness in biological control strains of Lactobacillus plantarum in plant environments. PLoS One 13:e0190931. doi: 10.1371/journal.pone.0190931

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Brien, P. A. (2017). Biological control of plant diseases. Australas. Plant Pathol. 46, 293–304. doi: 10.1007/s13313-017-0481-4

CrossRef Full Text | Google Scholar

Peng, G., Lahlali, R., Hwang, S. F., Pageau, D., Hynes, R. K., McDonald, M. R., et al. (2014). Crop rotation, cultivar resistance, and fungicides/biofungicides for managing clubroot (Plasmodiophora brassicae) on canola. Can. J. Plant. Pathol. 36, 99–112. doi: 10.1080/07060661.2013.860398

CrossRef Full Text | Google Scholar

Price, L. B., Liu, C. M., Johnson, K. E., Aziz, M., Lau, M. K., Bowers, J., et al. (2010). The effects of circumcision on the penis microbiome. PLoS One 5:e8422. doi: 10.1371/journal.pone.0008422

PubMed Abstract | CrossRef Full Text | Google Scholar

Raaijmakers, J. M., Paulitz, T. C., Steinberg, C., Alabouvette, C., and Moënne-Loccoz, Y. (2009). The rhizosphere: a playground and battlefield for soilborne pathogens and beneficial microorganisms. Plant. Soil. Fertil. Sci. China 321, 341–361. doi: 10.1007/s11104-008-9568-6

CrossRef Full Text | Google Scholar

Roesch, L. F., Fulthorpe, R. R., Riva, A., Casella, G., Hadwin, A. K., Kent, A. D., et al. (2007). Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283–290. doi: 10.1038/ismej.2007.53

PubMed Abstract | CrossRef Full Text | Google Scholar

Rui, W., Hongchun, Z., Liguang, S., Gaofu, Q., Shu, C., and Xiuyun, Z. (2017). Microbial community composition is related to soil biological and chemical properties and bacterial wilt outbreak. Sci. Rep. 7:343. doi: 10.1038/s41598-017-00472-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Shakeel, Q., Lyu, A., Zhang, J., Wu, M., Chen, S., Chen, W., et al. (2016). Optimization of the cultural medium and conditions for production of antifungal substances by Streptomyces platensis 3-10 and evaluation of its efficacy in suppression of clubroot disease (Plasmodiophora brassicae) of oilseed rape. Biol. Control 101, 59–68. doi: 10.1016/j.biocontrol.2016.06.007

CrossRef Full Text | Google Scholar

She, W., Sun, Z., Yi, L., Zhao, S., and Liang, Y. (2016). Streptomyces alfalfae sp. nov. and comparisons with its closest taxa Streptomyces silaceus, Streptomyces flavofungini and Streptomyces intermedius. Int. J. Syst. Evol. Microbiol. 66, 44–49. doi: 10.1099/ijsem.0.000671

PubMed Abstract | CrossRef Full Text | Google Scholar

Strelkov, S. E., Hwang, S.-F., Manolii, V. P., Cao, T., Fredua-Agyeman, R., Harding, M. W., et al. (2018). Virulence and pathotype classification of Plasmodiophora brassicae populations collected from clubroot resistant canola (Brassica napus) in Canada. Can. J. Plant. Pathol. 40, 284–298. doi: 10.1080/07060661.2018.1459851

CrossRef Full Text | Google Scholar

Sudini, H., Liles, M. R., Arias, C. R., Bowen, K. L., and Huettel, R. N. (2011). Exploring soil bacterial communities in different peanut-cropping sequences using multiple molecular approaches. Phytopathology 101, 819–827. doi: 10.1094/phyto-11-10-0310

PubMed Abstract | CrossRef Full Text | Google Scholar

van Elsas, J. D., Chiurazzi, M., Mallon, C. A., Elhottovā, D., Krištůfek, V., and Salles, J. F. (2012). Microbial diversity determines the invasion of soil by a bacterial pathogen. Proc. Nat. Acad. Sci. U.S.A. 109, 1159–1164. doi: 10.1073/pnas.1109326109

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Z., Gu, Y., Friman, V.-P., Kowalchuk, G. A., Xu, Y., Shen, Q., et al. (2019). Initial soil microbiome composition and functioning predetermine future plant health. Sci. Adv. 5:eaaw0759. doi: 10.1126/sciadv.aaw0759

PubMed Abstract | CrossRef Full Text | Google Scholar

Werner, J. J., Koren, O., Hugenholtz, P., DeSantis, T. Z., Walters, W. A., Caporaso, J. G., et al. (2012). Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys. ISME J. 6, 94–103. doi: 10.1038/ismej.2011.82

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, C., Penton, C. R., Zhu, C., Chen, H., Duan, Y., Peng, C., et al. (2018). Alterations in soil fungal community composition and network assemblage structure by different long-term fertilization regimes are correlated to the soil ionome. Biol. Fertil. Soils 54, 95–106. doi: 10.1007/s00374-017-1241-x

CrossRef Full Text | Google Scholar

Yang, W., Guo, Y., Wang, X., Chen, C., Hu, Y., Cheng, L., et al. (2017). Temporal variations of soil microbial community under compost addition in black soil of Northeast China. Appl. Soil Ecol. 121, 214–222. doi: 10.1016/j.apsoil.2017.10.005

CrossRef Full Text | Google Scholar

Yang, X.-x., Huang, X.-q., Wu, W.-x., Xiang, Y.-j., Du, L., Zhang, L., et al. (2020). Effects of different rotation patterns on the occurrence of clubroot disease and diversity of rhizosphere microbes. J. Integr. Agric. 19, 2265–2273. doi: 10.1016/s2095-3119(20)63186-0

CrossRef Full Text | Google Scholar

Zhang, T., Wang, N. F., Liu, H. Y., Zhang, Y. Q., and Yu, L. Y. (2016). Soil pH is a key determinant of soil fungal community composition in the Ny-Ålesund Region, Svalbard (High Arctic). Front. Microbiol. 7:227. doi: 10.3389/fmicb.2016.00227

PubMed Abstract | CrossRef Full Text | Google Scholar