Tan Wang1,2†
Kexin Yang3†Qingyun Ma1,2Xu Jiang1Yiqing Zhou1Delong Kong1Zhiye Wang4Rebecca E. Parales5
Lin Li2*Xin Zhao3*
Zhiyong Ruan1,6,7*- 1CAAS-CIAT Joint Laboratory in Advanced Technologies for Sustainable Agriculture, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China
- 2State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan, China
- 3Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China
- 4Key Laboratory of Microbial Resources Exploitation and Application of Gansu Province, Institute of Biology, Gansu Academy of Sciences, Lanzhou, China
- 5Department of Microbiology and Molecular Genetics, College of Biological Sciences, University of California, Davis, Davis, CA, United States
- 6College of Resources and Environment, Tibet Agricultural and Animal Husbandry University, Linzhi, China
- 7College of Life Sciences, Yantai University, Yantai, China
As an ornamental flower crop, the long-term continuous monocropping of cut chrysanthemum causes frequent occurrence of diseases, seriously affecting the quality of cut chrysanthemum. The rhizosphere microbial community plays an important role in maintaining the healthy growth of plants, whereas the composition and dynamics of rhizosphere microbial community under continuous monocropping of cut chrysanthemum have not been fully revealed. In this study, the Illumina MiSeq high-throughput sequencing platform was used to monitor the dynamic changes of rhizosphere microbial communities in four varieties of cut chrysanthemum during 0–3 years of monocropping, and the soil physicochemical properties were also determined. Results showed that continuous monocropping significantly increased the fungal community richness and altered the profiles of the bacterial and fungal communities, leading to variation of community beta-diversity. With the increase of continuous cropping time, biocontrol bacteria decreased, while some plant pathogenic fungi were enriched in the rhizosphere of cut chrysanthemum. FAPROTAX-based functional prediction showed that the abundance of gene related to nitrogen and sulfur metabolism and chitin lysis was reduced in the rhizosphere of cut chrysanthemum. FUNGuild-based fungal function prediction showed that plant pathogenic fungal taxa were increasing in the rhizosphere of cut chrysanthemum, mainly Acremonium, Plectosphaerellaceae, Fusarium, and Cladosporium. Continuous cropping also reduced the content of ammonium nitrogen and increased soil salinity, resulting in deterioration of soil physical and chemical properties, which, together with the transformation of rhizosphere microbial community, became part of the reasons for the continuous cropping obstacle of cut chrysanthemum.
Introduction
Cut chrysanthemum is a popular flower crop because of its easy breeding, long bottle life, and storage transport. It is also known as a major cut flower worldwide, and its main production areas cover Netherlands, China, Japan, and Korea (Li et al., 2016). In China, the planting area of cut chrysanthemum exceeds 7,000 hectares, and the annual sales reach 2.65 billion sprays with a sales volume of 1.52 billion yuan (China Ministry of Agriculture, 2016). Cut chrysanthemum is generally grown in greenhouse sheds, which can be harvested up to five times a year, to obtain great economic benefits. However, similar to other crops, such as soybean, tomato, watermelon, and vanilla, long-term continuous monocropping will generate replanting problems. The replanting problem of cut chrysanthemum will induce adverse symptoms such as short stature, yellowing, and increased infection rate (Li et al., 2017). In commercial production, the replanting problem has led to a significant decline in the yield and quality of cut chrysanthemum, resulting in great economic losses for growers.
At present, a number of measures to address the replanting problem have been developed. Rotation can reduce weed and pest occurrence, facilitate soil nutrient utilization (Nevens and Reheul, 2001), and improve crop yield and quality (Zhao et al., 2020). Organic fertilizer contains a variety of organic acids (OA), peptides, and rich nutrient elements such as nitrogen, phosphorus, and potassium, which can not only provide comprehensive nutriment for crops but also increase and renew soil organic matter, promote microbial reproduction, and improve soil physical and chemical properties and biological activity (Chen et al., 2018). Long-term continuous monocropping will cause the accumulation of pathogenic microorganisms, and the application of soil sterilization (including high-temperature sterilization and reagent fumigation) can effectively reduce the number of pathogenic microorganisms (Li et al., 2019). However, soil sterilization may damage the soil microbial ecology, which is contrary to the objective of sustainable agriculture. Rhizosphere microbial communities are closely related to plant health (Raaijmakers et al., 2009). The root exudates from plants consist of phenolic acids (PA), OA, sugars, amino acids, fatty acids, alkaloids, terpenoids, enzymes, flavonoids, and proteins (De Weert et al., 2002; Bais et al., 2006). A variety of carbon and nitrogen sources obtain abundant microorganisms. Exploiting microorganisms isolated from the environment as antagonists or growth-promoting rhizobacteria can effectively reduce the occurrence of diseases and promote plant growth (Verbon and Liberman, 2016; Morales-Cedeño et al., 2021).
Nevertheless, the aforementioned method always transforms the soil microbiome to a healthier direction. As the second genome of plants, plant-associated microbiome plays an important role in plant growth. With regard to nutrient utilization, the symbiotic interactions between plants and arbuscular mycorrhizal fungi and Rhizobium bacteria lead to nutrient acquisition (Richardson and Simpson, 2011; Corné et al., 2014; Trivedi et al., 2016; Rachel et al., 2018). Non-symbiotic plant-growth-promoting bacteria can improve not only the bioavailability of insoluble minerals but also the root system architecture of host plants, thereby increasing the exploratory capacity of the root for water and minerals (Richardson and Simpson, 2011; Trivedi et al., 2016). A recent study suggested that the differences in nitrogen use efficiency in rice varieties are due to a high proportion of nitrogen-cycling-related bacterial supplementation, resulting in a more efficient transformation of nitrogen in the root environment in indica rice than in japonica rice varieties (Zhang et al., 2019). Rhizosphere microorganisms extend the bioavailability of nitrogen by nitrification, delay flowering time by converting tryptophan to the plant hormone indoleacetic acid, and stimulate plant growth by downregulating flowering genes (Lu et al., 2018). Plant-related microbial groups can mobilize nutrients that are not easily available to plants, such as inorganic phosphate and iron, through dissolution, mineralization, or excretion through iron-chelated siderophores. Moreover, plant-related microbiome helps plants fight diseases. The effect of plant defense based on natural microorganisms on plant health has been clearly demonstrated in suppressive soils. Plant root exudates stimulate, enrich, and support soil microorganisms as the first line of defense against soil-borne pathogens (Mendes et al., 2011). Community-based analyses of suppressive soils have demonstrated that no single phylum is uniquely associated with disease suppression (Mendes et al., 2011; Penton et al., 2014; Santhanam et al., 2015; Cha et al., 2016; Trivedi et al., 2017). A recent study has modeled general disease suppression of Fusarium oxysporum and identified bacterial abundance in Actinomycetes and Firmicutes as predictive markers for continent-scale soil disease suppression. Under abiotic or biological stress conditions, plants usually select a microbial community that promotes stress resistance (Wagner et al., 2014; Eida et al., 2018; Naylor and Coleman-Derr, 2018; Naylor et al., 2018; Timm et al., 2018; Fitzpatrick et al., 2019), which modifies plant evolutionary responses to environmental stress by altering the fitness of individual plant genotypes, the expression of plant traits related to fitness, and the strength or direction of natural selection occurring within populations that experience environmental stress through the effects of microorganisms on reproductive fitness (Lau and Lennon, 2012; Eida et al., 2018; Fitzpatrick et al., 2018, 2019). In the effect of fungal endophytes on plant growth performance under water stress, traits related to resource utilization and stress tolerance account for 26–53% (Kudjordjie et al., 2019). In addition, the complex microbial community can help plants deal with other adverse effects, such as heavy metals (Yang et al., 2020), herbicides (Fawzy et al., 2015), and pesticides (Li et al., 2020). The composition and structure of the microbiome must remain stable to continuously exert the functional potential of the plant-related microbiome, and the loss of its steady state may cause negative effects such as imbalance in soil nutrient conversion and the proliferation of pathogenic bacteria. Therefore, determining the dynamic changes of the rhizosphere microbial community under different continuous cropping periods is necessary to explore the mechanism of continuous cropping obstacles of cut chrysanthemum.
Based on previous studies, the rhizosphere microbial community of cut chrysanthemum is based on denaturing gradient gel electrophoresis (Duineveld et al., 2001; Song et al., 2013), which is not sufficient to reveal the overall dynamic changes of bacteria and fungi during continuous cropping of cut chrysanthemum. With the wide application of high-throughput sequencing technology, the microbial community landscape in the environment is gradually revealed. Here, we used Illumina MiSeq sequencing technology to detect the dynamic changes of microbial communities in the rhizosphere soil of four varieties of cut chrysanthemum during 0–3 years of monocropping period and measured the physical and chemical properties, nutrients, and other indicators of soil to evaluate the health status of the soil during continuous cropping. We hypothesized that (i) during continuous monocropping for 0–3 years, harmful microorganisms would accumulate in the rhizosphere. (ii) Differences in rhizosphere exudates of different varieties of cut chrysanthemum could affect the rhizosphere microbial community. (iii) With the increase of continuous monocropping time, soil physicochemical properties decline, and nutrient conversion slows down.
Materials and Methods
Study Site and Experimental Design
The experimental site was located in Nankou experimental base of the Chinese Academy of Agricultural Sciences, Beijing (116.57033, 40.237361), China. Cut chrysanthemums with different monocropping times were planted in different greenhouses at a distance of 8 m from each other. Each greenhouse is 80 m in length. In the greenhouse, each variety of cut chrysanthemum was randomly planted in three plots, separated by other cut chrysanthemum varieties. All plots adopted the same water and fertilizer management practices. During sampling, each plot randomly took soil samples of three chrysanthemum plants and mixed them into one replicate. The monoculture time selected in the experiment was 0 years (CK, before planting cut chrysanthemum), 0.5 years (H, planting for one crop), and 3 years (T, planting for six crops). The sampling period was the mature period. A total of four varieties were selected, namely, Glitter Orange (variety 1), Gem Yellow (variety 2), Cloud Red (variety 3), and Mona Lisa Pink (variety 4). Thus, the sample names were CK (0 years); H1, H2, H3, and H4 (0.5 years); and T1, T2, T3, and T4 (3 years). CK had four replicates; the others had three replicates. In obtaining rhizosphere soil, the fibrous roots and main roots of the plant were cut off with sterile scissors, and then the roots were placed into a sterilized 50 ml centrifuge tube. Subsequently, the root tissues were added with 0.2 × PBS buffer solution and swirled for 1 min. Afterward, the root tissues were removed and centrifuged at 8,000 rpm for 5 min. The supernate was discarded, and the residual part was the rhizosphere soil. The rhizosphere soil samples were placed at –80°C for subsequent DNA extraction.
Soil Physicochemical Properties
The planted soil of cut chrysanthemum (bulk soil) was collected to determine the physical and chemical properties. The sampling method of bulk soil referred to that of Edwards et al. (2014). The collected soil was air-dried on a kraft paper and screened by a 1 mm sieve, and the resulting soil was determined as dried soil. Soil pH was measured using a glass conductance meter in a 1:2.5 soil/water (W/V) suspension (Qiu et al., 2012). The content of soil total carbon, total nitrogen, and organic carbon was determined by using an elemental analyzer (Elementar Analysensysteme GmbH, Germany). Soil urease (Kandeler and Gerber, 1988), catalase (Johansson, 1988), sucrase (Gao et al., 2013), and alkaline phosphatase (Scandinavica et al., 2014) activities were determined using the soil enzyme activity assay kit from Solarbio Science & Technology Co. (Beijing, China). Soil available P was extracted using sodium bicarbonate and measured with the molybdenum blue method (Olsen, 1954). Soil available K was extracted using ammonium acetate and determined with flame photometry (nova 300, Analytic Jena, Germany) (Frank et al., 1998). Soil ammonium nitrogen and nitrate nitrogen were extracted with 1 M KCl and the indophenol-blue colorimetric and double wavelength (220 and 275 nm) methods, respectively, and their concentrations were measured using a spectrophotometer (UV-6000, China) (Zhao et al., 2017). Soil salinity was determined by an oven-drying method following the NY/T 1121.16-2006 standard water-soluble salt analysis of soil.
Soil DNA Extraction and High-Throughput Amplicon Sequencing
Total genome DNA from 0.5 g of fresh soil samples was extracted using a E.Z.N.A™ Mag-Bind Soil DNA Kit (OMEGA, Norcross, United States) following the manufacturer’s instructions. DNA integrity was detected on 1% agarose gel, and DNA concentration was quantitatively detected by Qubit. The V3–V4 region of the 16S rRNA gene was amplified using the PCR primers 341F (5′-CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTACHVGGGTATCTAATCC-3′) for bacterial community analysis, and the ITS1 region of the transcribed spacer was amplified using the PCR primers ITS1F (5′-CTTGGTCATTTA GAGGAAGTAA-3′) and ITS2R (5′-GCTGCGTTCTTCATCGATGC-3′). 16S rRNA genes were amplified using a specific primer with a barcode. The first round of PCRs was performed in 30 μl of reactions with 15 μl of 2 × Hieff Robust PCR Master Mix, 1 μl of Bar-PCR primer F and 1 μl of Primer R, and 10–20 ng of template DNA. PCR amplification cycles consisted of an initial denaturation at 94°C for 3 min, followed by 5 cycles of denaturation at 94°C for 30 s, annealing at 45°C for 20 s, and elongation at 65°C for 30 s; then 20 cycles of denaturation at 94°C for 20 s, annealing at 55°C for 20 s, and elongation at 72°C for 30 s; and finally elongation at 72°C for 5 min after cycles.
Illumina bridge PCR-compatible primers were introduced in the second round of amplification. The reaction system was identical to that of the first round; however, the primers were replaced by Primer F and Index-PCR Primer R. Library quality was assessed using a Qubit@ 3.0 fluorometer (Thermo Scientific, Waltham, MA, United States). Finally, the library was sequenced on an Illumina MiSeq platform, and 250 bp paired-end reads were generated.
Community Diversity and Function Analysis
USearch was used to cluster the sequences in the libraries. Operational taxonomic unit (OTU) clustering was performed on non-repeated sequences (excluding single sequences) based on 97% similarity, and the representative sequence of OTU was obtained by removing chimeras during clustering. Then, an RDP classifier was used to annotate taxonomic information for each representative sequence. In calculating alpha diversity, we refined the OTU table and calculated two indicators, namely, ACE index and Shannon even index, to estimate the species richness and species evenness of the community, respectively.
The 10 most abundant phyla were selected from all the samples included in the phylum-level species annotation and abundance information, and a histogram was drawn for each sample with regard to relative abundance. In addition, 30 genera with the highest abundance were selected to draw a heatmap. LEfSe (Segata et al., 2011) was used to compare data among samples and select biomarkers for each sample, and 2.0 was set as the threshold for the logarithmic LDA score for discriminating features. Principal co-ordinate analysis (PCoA) was performed to compare the community profiles of different samples in R1 by using the amplicon and phyloseq packages. Redundancy analysis (RDA) was performed to evaluate the relationship between microbial community profiles and soil physical and chemical properties by using the ggvegan package. Linear fitting analysis was used to evaluate the correlation between soil enzyme activity and soil nutrient element content and between soil enzyme activity and microbial community richness index. To further explore the impact of microbial community change, functional prediction with PICRUSt (Langille et al., 2013), FAPROTAX (Louca et al., 2016), and FUNGuild (Nguyen et al., 2015) annotation tools was performed. The predicted results were visualized by a principal component analysis (PCA) plot and boxplot. The figures were adjusted, combined, and modified by Adobe Illustrator 2021.
Results
Changes of Soil Physicochemical Properties of Cut Chrysanthemum During Continuous Monoculture
During continuous monocropping for 0–3 years, the contents of soil total carbon, total nitrogen, and organic carbon increased significantly in soil planted with “variety 1” and “variety 4,” whereas no evident change was observed in soil pH of all varieties (Table 1). Ammonium nitrogen content showed a decreasing trend in bulk soil of all varieties. In general, the content of nitrate nitrogen initially increased and then decreased during continuous monocropping, but it increased remarkably in soil planted with “variety 1” in the third year of monocropping. The content of available phosphorus and potassium also showed an increasing trend in replanted soil, and the increment in bulk soil planted with “variety 1” and “variety 4” was high. Likewise, the total salt content in soil increased in the third year of monocropping, and the maximum increment was observed in soil planted with “variety 1.”
Table 1. A summary of soil physicochemical properties in bulk soil of cut chrysanthemum.
During monocropping of cut chrysanthemum, the activities of four soil enzymes showed an upward trend in general. Urease activity of soil planted with “variety 1,” “variety 3,” and “variety 4” significantly increased in the third year, and a high increment was observed in bulk soil planted with “variety 1” (Figure 1). The activity of soil catalase also showed an increasing trend with the increase of monocropping time. The activity of sucrase was almost zero from 0 to 0.5 years, but it increased greatly in the third year (T1). Meanwhile, alkaline phosphatase activity increased from 0.5 to 3 years.
Figure 1. Changes of the activities of four pivotal soil enzymes with the continuous monoculture period (Tukey’s HSD test). (A) Urease. (B) Catalase. (C) Sucrase. (D) Alkaline phosphatase.
Diversity Changes of Rhizosphere Microbial Community With Continuous Monocropping Period
During the continuous monocropping period of 0–3 years, the richness and evenness of the rhizosphere bacterial community of the four varieties of cut chrysanthemum did not change significantly (Figures 2A,B). The richness of the fungal community increased remarkably in the third year, whereas the evenness remained constant (Figures 2C,D). During the monoculture period, the structure and composition of the rhizosphere microbial community changed dramatically. At the phylum level of bacteria (Figure 3A), Proteobacteria with the highest relative abundance maintained stable abundance with continuous monocropping. In another dominant phylum, Actinobacteria, its abundance decreased significantly in the third year of continuous monocropping, particularly in the rhizosphere of “variety 2.” On the contrary, the abundance of Bacteroidetes (Top4 phylum) showed a significant increase in the third year. Furthermore, the abundance of other dominant phyla, such as Verrucomicrobia, Gemmatimonadetes, and Chloroflexi, remained relatively stable in the rhizosphere. In the fungal community (Figure 3B), Ascomycota was the most abundant phylum. In year 0.5, the relative abundance of Ascomycota varied greatly in the rhizosphere of the four varieties of cut chrysanthemum, but it remained stable in the rhizosphere of the four varieties during the follow-up period. Basidiomycota (Top3 phylum) gradually enriched around the root with continuous monocropping, and it showed preference to “variety 1.” The relative abundance of Mortierellomycota (Top4 phylum) in CK soil was higher than that in the rhizosphere, and it occupied great proportion in the rhizosphere with continuous monocropping. At the genus level, Arthrobacter (Top1 genus) was largely enriched in the rhizosphere of 0.5 year samples, but it declined in the third year (Figure 3C). The relative abundance of Pseudomonas (Top4 genus) and Lysobacter (Top24 genus) also decreased in the third year of monoculture. Sphingobacterium (Top7 genus) was significantly enriched in the rhizosphere of all varieties. Some taxa, such as Gp6 (Top5 genus), Streptophyta (Top9 genus), and Ohtaekwangia (Top18 genus), only showed differences in some samples. In the fungal community (Figure 3D), Cladosporium (Top1 genus) was greatly enriched in the rhizosphere of cut chrysanthemum with continuous monoculture. Acremonium (Top2 genus) and Thyronectria (Top14 genus) showed the same trend. An unclassified genus (Top29), and Botryotrichum (Top30 genus), showed high abundance in the rhizosphere of cut chrysanthemum in 0.5 years. However, by the third year of monoculture, their abundance dropped dramatically.
Figure 2. Variations of the rhizosphere microbial community richness and evenness in different monoculture years. The ACE index indicates the richness of community, where the Shannon even index represents the evenness of the community. Box colors show the length of continuous monoculture time. (A,B) Prokaryotic community. (C,D) Fungal community.
Figure 3. The dynamics of the rhizosphere microbial community composition on phylum and genus levels during the continuous monoculture time series. (A) Bacterial community on phylum level. (B) Fungal community on phylum level. (C) Bacterial community on genus level. (D) Fungal community on genus level.
LEfSe identified the taxa that were significantly different among samples, which were defined as biomarkers (Figure 4). In the bacterial community, the abovementioned bacteria changed significantly in the third year of monoculture, which primarily belonged to Bacteroidetes and Acidobacteria (Figure 4A). At a low classification level, the families Sphingobacteriaceae (Top5 family), norank_Acidobacteria_Gp4 (Top7 family), Chitinophagaceae (Top4 family), Flavobacteriaceae (Top25 family), Planctomycetaceae (Top9 family), and Enterobacteriaceae (Top63 family) and the genera Sphingobacterium (Top6 genus), Gp6 (Top3 genus), Chryseobacterium (Top37 genus), Chryseolinea (Top29 genus), Gp7 (Top18 genus), Sphingopyxis (Top31 genus), Pelagibacterium (Top46 genus), Rhizobium (Top40 genus), Variovorax (Top33 genus), Flavobacterium (Top61 genus), Olivibacter (Top63 genus), Halomonas (Top51 genus), Stenotrophomonas (Top82 genus), Buttiauxella (Top106 genus), and Gp10 (Top35 genus) contributed to the occurrence of differences in the phylum level. For the fungal community, the difference in the rhizosphere community in the third year of monoculturing primarily occurred on Ascomycota and Rozellomycota (Figure 4B). At the species level, Cladosporium_sphaerospermum (Top1 species), unclassified_Acremonium (Top3 species), unclassified_Gibellulopsis (Top5 species), Acrostalagmus_luteoalbus (Top4 species), Thyronectria rhodochlora (Top7 species), and Podospora pyriformis (Top10 species) had high LDA scores, which partly resulted in the difference at the phylum level.
Figure 4. LEfSe of bacterial 16S rDNA (A) and fungal ITS rDNA (B) sequences with different abundances between continuous monoculture periods. The non-significantly different species are uniformly colored yellow, and the different species biomarkers are colored the same as the groups. The names of species represented by letters are shown in the legend on the bottom.
Principal co-ordinate analysis based on the Bray–Curtis distance showed that bacterial or fungal communities of different monocropping periods were separated on the first axis, and the differences accounted for 31% (Figure 5A) and 37% (Figure 5B), respectively. The PERMANOVA test showed significant differences among the three-period samples (P < 0.01). RDA of the rhizosphere microbial community and soil physicochemical properties during monocropping showed that (Figures 6A–C) the total contribution rates of soil total carbon (C), total nitrogen (N), organic carbon (OC), ammonium nitrogen (AN), nitrate nitrogen (NN), available phosphorus (AP), available potassium (AK), and soil salinity (SC) in shaping bacterial and fungal communities were 44.9 and 46.9%, respectively, of which the content of ammonium nitrogen had the highest contribution rate (13.0% for bacteria and 14.0% for fungi). On the first axis, the microbial community was primarily separated by AN, C, SC, and NN, whereas on the second axis, the microbial community was separated by OC, N, AP, and AK.
Figure 5. Difference of the rhizosphere microbial community profiles of various monoculture years. Points’ color represents different monoculture periods. The points are covered by a 95% confidence ellipse. Pairwise PERMANOVA is shown below the PCoA plots. (A) PCoA of prokaryotic communities. (B) PCoA of fungal communities.
Figure 6. Correlations between soil physicochemical properties and rhizosphere microbial species. (A,B) RDA of bacterial and fungal community profiles correlated to soil properties. (C) Explanation rate of soil physicochemical properties to the differences of community profiles (conditional effect).
Functional Prediction of Rhizosphere Microbial Community in Cut Chrysanthemum During Continuous Monocropping
Principal component analysis was performed with the results of COG and KEGG function annotations (Figure 7), and the results showed that the functional composition of rhizosphere microorganisms of cut chrysanthemum has changed significantly from 0.5 to 3 years (PERMANOVA, COG: P = 0.007, KO: P = 0.007). The five COG classifiers with the highest contribution rates that cause differences between samples were COG0583, COG1629, COG1595, COG2204, and COG0515. According to the distribution of sample points of different monocropping years on the PCA plot, it showed that COG1629 was the most influential COG classifier that caused differences of functional profiles. COG1629 was annotated as outer membrane receptor proteins, mostly referring to Fe transport. In the KEGG annotated results, the five most volatile KEGG classifications were K02026, K02004, K03088, K02027, and K02025. The K03088, belonging to RNA polymerase sigma-70 factor, ECF subfamily, had the greatest contribution to differences between groups (0.5 vs. 3 years).
Figure 7. Principal component analysis based on COG (A) and KEGG (B) functional annotations. The arrows represent the most volatile COG and KEGG classifiers in the functional prediction results. PERMANOVA was used to test the significance of the differences.
Based on the results of FAPROTAX and FUNGuild, significant changes in the abundance of some metabolic pathways of bacterial community (Figure 8A) and the plant pathogen of fungal community (Figure 8B) had occurred. Among these pathways, most of them involved nitrogen and sulfur metabolism, and one was associated with chitin lysis. Continuous monocropping obviously reduced the abundance of these metabolic genes. For the fungal community, all species were matched with different functional groups (Guild). Guilds which were part of plant pathogens had been observed in the rhizosphere, and some of them were massively enriched with continuous monocropping.
Figure 8. Changes in relative abundance of metabolic pathways in prokaryotic communities based on FAPROTAX (A) and variation of abundance of plant pathogen guilds based on FUNGuild (B). ANOVA was used to test the significance of differences.
Discussion
Continuous Monocropping Resulted in Soil Degradation of Cut Chrysanthemum
Here, we monitored the changes of soil-related physical and chemical properties of four varieties of cut chrysanthemum after continuous monocropping for 3 years, and the results showed that the content of soil carbon, nitrogen, available phosphorus, and available potassium showed different degrees of increment in the third year of monocropping (Table 1). The soil pH did not change dramatically, but a slight reduction occurred in the 3 year soil planted with “variety 1” and “variety 4.” Notably, the content of available nitrogen (including ammonium nitrogen and nitrate nitrogen) in the soil basically decreased with continuous cropping, which may cause nutrient imbalance and hinder the growth of cut chrysanthemum. In this study, the salt content in the soil showed an increasing trend with continuous monocropping. Soil salinity is also closely related to the growth of plants. Higher soil salinity produces higher osmotic pressure, which leads to the inability of crops to effectively absorb water and nutrients in the soil (Kluitenberg and Biggar, 1992; Tian et al., 2020), resulting in vegetable deficiency symptoms and decreasing crop quality and yield (Riaz et al., 2018). Soil enzyme, which is produced by plants, animals, and microorganisms, participates in the synthesis and oxygenolysis of various organic and inorganic substances, and it is an indicator for evaluating soil fertility (Dick and Kandeler, 2004; Wang et al., 2019). Soil urease (S-UE) is an enzyme that hydrolyzes amide organic nitrides into inorganic nitrides that can be directly utilized by plants (Singh and Ghoshal, 2013; Cai et al., 2015), while high urease activity rapidly hydrolyzes applied urea to ammonia which contributes to soil nitrogen losses and reduces N use efficiency of crop plants (Rana et al., 2021). Furthermore, ammonium accumulation in soil, particularly at high pH, may hinder the nitrification process at the midway and result in the accumulation of toxic levels of nitrites (Breuillin-sessoms et al., 2017). Soil catalase (S-CAT) is an important oxidoreductase, which is closely related to the activity of oxygen-consuming microbes. It can promote the decomposition of hydrogen peroxide—a metabolic intermediate—and alleviate its toxicity (Sun et al., 2013). The enzymatic activity of catalase can characterize the strength of soil humus and the conversion rate of organic matter (Lee et al., 2009). The soil sucrase activity (S-SC) reflects the maturation degree and fertility level of soil and is an important index for evaluating soil fertility (Stępniewska et al., 2009). Further, sucrase activity is directly related to the growth of crops and closely related to the transformation of organic matter and respiratory intensity (Gianfreda et al., 2005; Tripathi et al., 2007). Alkaline phosphatase enzymes (S-ALP) are believed to be mainly released by microorganisms (Nannipieri et al., 2011) and to be affected by plant species, microbial community composition, and rhizodeposition (Badri and Vivanco, 2009; Wasaki et al., 2018; Teixeira et al., 2019). Moreover, alkaline phosphatase enzymes are ubiquitous in soils, but their activity is somehow influenced by the quantity and quality of rhizodeposits (Ragot, 2016; Wasaki et al., 2018; Wu et al., 2018; Wei et al., 2019). In the current research, the activities of four soil enzymes all increased with continuous monocropping (Figure 1). Part of possible reasons were revealed by performing the linear fitting analysis of enzyme activity and soil physicochemical properties (Figure 9A) and ACE diversity index of microbial community (Figure 9B). Strong correlations were observed in S-UE and N (R2 = 0.36, P < 0.001), S-CAT and OC (R2 = 0.24, P = 0.008), S-SC and OC (R2 = 0.27, P = 0.005), and S-ALP and AP (R2 = 0.37, P < 0.001). The residual amount of nutrient elements differentiated due to the different utilization capacity of each cut chrysanthemum, and higher residual amounts always corresponded to higher enzyme activity. Still, the ACE index of fungal community increased with the increase of enzyme activity (Figures 9Be–h), which may indicate that cut chrysanthemum recruited more fungal taxa to increase enzyme production. The continuous increase of urease activity would be one of the reasons for the occurrence of replanting problems.
Figure 9. Linear fitting analysis of soil enzyme activity and soil physical and chemical properties (Aa–Ad) and microbial community richness index (Ba–Bh). R2 represented the magnitude of the correlation coefficient, and the p-value was generated by a t-test.
Continuous Monocropping Caused Changes in the Bacterial Community Composition and a Decrease in the Abundance of Nutrient Element Conversion Genes
The richness and evenness of bacterial community did not change significantly during continuous monocropping of the four varieties (Figures 2A,B). However, variations in the abundance of some taxa occurred in the rhizosphere of cut chrysanthemum (Figures 3A,C, 4A). Proteobacteria is the dominant taxon in different geographical regions and soil types (Janssen, 2006), and members of this phylum play important roles in phylogenetic, ecological, and pathological values and energy metabolism (Bryant and Frigaard, 2006; Mukhopadhya et al., 2012). Root exudates and plant species play a key role in shaping the rhizosphere microbial community, which leads to genotype-specific communities of plants within the same soil type (Grayston et al., 1998; Marschner et al., 2004; Chaparro et al., 2014). The differences in the abundance of Proteobacteria may be primarily attributed to the varieties of cut chrysanthemum (Figure 3A). Based on previous studies, Actinobacteria is associated with pathogen antagonism (Mendes et al., 2011; Xiong et al., 2017), which is enriched in large amounts in suppressive soil and significantly decreased with the increase of monocropping time (Xiong et al., 2015). In our study, the relative abundance of Actinobacteria reduced significantly in the third year of monoculture. Acidobacteria is also a member of a relatively rich and prevalent soil bacterial community (Janssen, 2006), and it plays an important role in organic matter degradation and nutrient cycling (Eichorst et al., 2018). Fierer et al. (2007) showed that Acidobacteria had the highest abundance in soil with low resource availability. In general, continuous monocropping reduces soil nutrients and soil pH, which may explain the high abundance of Acidobacteria in replanted soil of cut chrysanthemum. In addition, Bacteroidetes participates in carbon and nitrogen metabolism, such as organic matter degradation and nitrite oxidation (Clemens, 2017; Han et al., 2017; Eichorst et al., 2018; Fang et al., 2018), and it has higher abundance in environments with high carbon content. The proportion of Bacteroidetes in replanted soil increased probably because more debris were produced in the roots of cut chrysanthemum. Other dominant phyla, Acidobacteria and Chloroflexi, also play an important role in organic matter degradation and nutrient cycling (Eichorst et al., 2018; Fang et al., 2018). For low classification levels, some species of Pseudomonas are involved in plant pathogen inhibition in many soils, and it serves as PGPR to reduce plant diseases (Liu et al., 1995; Kinkel et al., 2012). In the rhizosphere of cut chrysanthemum, the abundance of Pseudomonas decreased remarkably in the third year of continuous monocropping. Many Bacillus species, such as Bacillus drentensis, Bacillus simplex, and Bacillus aryabhattai, secrete fengycin, bacillaene, difficidin, and iturins by nutritional competition to antagonize pathogens such as Rhizoctonia solani and Botrytis cinerea (Ongena and Jacques, 2008; Jeong et al., 2012). The relative abundance of Bacillus did not change significantly in our study. Rhizobium, as a nitrogen fixer, was largely enriched in the rhizosphere of cut chrysanthemum, which may be attributed to the consumption of ammonium nitrogen in bulk soil. However, it should not be ignored that the species of Pseudomonas are not all beneficial; some are even plant pathogenic bacteria. In order to further improve the recognition rate of species and strains, the non-clustering algorithm of DADA3 may be better than that of UPARSE.
Changes in the composition of the bacterial community led to variations in metabolic genes, which were predicted by FAPROTAX based on 16S rRNA sequences. Continuous monocropping significantly reduced the abundance of genes related to nitrogen and sulfur metabolism in the rhizosphere of cut chrysanthemum (Figure 8A). Nitrogen is an important component of plant cell proteins, nucleic acids, chlorophyll, enzymes, etc., which can improve photosynthetic efficiency and increase crop yield and quality. As an element equally important to plants as NPK, sulfur is essential for the function of plant proteins, enzymes, and other physiologically active substances. The lack of sulfur element may lead to smaller plants and yellow leaves, which are symptoms of nitrogen deficiency, and ultimately result in the loss of crop yield and quality. With the increase of continuous cropping years, the reduction of the related metabolic genes will reduce the conversion rate of the rhizosphere nutrient of cut chrysanthemum and gradually cause a state of “element deficiency.” This state may contribute to the occurrence of reduced buds, smaller flowers, and increased diseases in cut chrysanthemum.
Continuous Monocropping Significantly Increased the Number of Phytopathogenic Fungi in the Rhizosphere of Cut Chrysanthemum
During continuous monocropping, the richness of the rhizosphere fungal community significantly increased (Figure 2C). LEfSe showed significant changes in the abundance of multiple fungal taxa (Figure 4B). Ascomycota is the most dominant phylum in the rhizosphere of cut chrysanthemum. It contains many plant pathogens and is the main source of toxins that cause replanting diseases (Ecker et al., 2005; Li et al., 2014). It was significantly enriched in the rhizosphere of cut chrysanthemum with the increase of monocropping time. The third dominant phylum, Basidiomycota, has some species that are saprophytic, and many of its large saprophytic fungi cause wood rot. Nevertheless, other species can also promote the growth of plants by forming mycorrhizas with plants. In this study, Basidiomycota was significantly enriched in the rhizosphere of “variety 1.” At the genus level, Mortierella was involved in plant pathogen inhibition, which can inhibit the occurrence of the Chinese cabbage bulb disease and protect bananas against Fusarium wilt and root rot (Narisawa et al., 1998). It has higher abundance in vanilla suppressive soil, while it gradually accumulated in the rhizosphere of cut chrysanthemum with increasing planting time. Humicola also served as a biocontrol fungus (Sun et al., 2014), whose relative abundance increased significantly in the rhizosphere of cut chrysanthemum. Pathogenic fungi of some plant diseases, such as Alternaria, which can cause black spot in cabbage and Alternaria Panax in ginseng (Zhao et al., 1993; Ko et al., 2011), and Fusarium, which can cause wilt in chrysanthemums (Li et al., 2017), were enriched in the rhizosphere of cut chrysanthemum. Rhizoctonia, a pathogen that can cause blight in chrysanthemums, watermelons, etc. (Crater, 1992; Fehér, 1993), did not change significantly in our study.
To further reveal the changes of fungal community function, the FUNGuild method was used to divide fungal species into multiple guilds. A guild represents a more subdivided nutrition mode based on pathotroph, saprotroph, and symbiotroph, one of which is a plant pathogen. Guilds with significant changes were detected based on ANOVA. We found that some of the guilds belonging to plant pathogens were significantly enriched in the rhizosphere of cut chrysanthemum (Figure 8B). According to FUNGuild results, these guilds mostly belong to Acremonium (Hou et al., 2019), Plectosphaerellaceae (Giraldo and Crous, 2019), Fusarium (Ma et al., 2013), and Cladosporium (Carolina et al., 2021), which have been reported to be pathogens of certain crops. The pathogenic genus Acremonium (corresponding to Guild: Plant Pathogen D) has a relative abundance that even reached 50% in the rhizosphere of cut chrysanthemum (Top2 genus). Therefore, we speculated that the increased incidence of diseases during continuous monocropping of cut chrysanthemum was caused by not only Fusarium but also the joint action of multiple pathogens.
In the current research, four cultivars of cut chrysanthemum were selected to investigate the changes of the rhizosphere microbial community and soil physicochemical properties under different continuous cropping times. Continuous monoculture had less effect on the bacterial community, but it significantly reduced the number of bacteria that metabolized nitrogen and sulfur. Continuous cropping significantly increased the richness of rhizosphere fungi and decreased the abundance of chitin-lysis-related genes, prompting the rhizosphere of cut chrysanthemum to become a fungal community. The taxa largely enriched in rhizosphere belonged to phytopathogenic fungi, which may combine to disrupt the otherwise healthy microbial community and increase diseases in cut chrysanthemum. Continuous cropping did not always increase the number of pathogenic bacteria; some beneficial microorganisms were also enriched and promoted the increase of soil enzyme activity. But they seem unable to resist the trend of community degradation. In this study, ammonium nitrogen (AN) was one of the most important factors shaping microbial communities, and the variation of its content caused significant changes in bacterial and fungal communities. However, the total contribution rate of determined soil physicochemical properties to the community transformation was less than 50%, which suggested that there were other factors involved in shaping the rhizosphere microbial community. Allelochemicals are produced by leaves, roots, stems, etc., which have beneficial or adverse effects on other plants (Rice, 1984). Allelochemicals consist of OA, straight-chain alcohols, aldehydes or ketones, unsaturated lactones, fatty acids, etc. (Gross and Parthier, 1994; Einhellig, 1995; Rice, 1995; Seigler, 1996). If the donor and the recipient belong to the same species, this leads to intraspecific allelopathy, and the term used is autotoxicity (Singh et al., 2010). Autotoxicity is a major obstacle in continuous cropping under protected cultivation (Zhou et al., 2009), and it is also found in the cultivation of cut chrysanthemum. A variety of allelochemicals participate in the transformation of the composition and structure of the rhizosphere microbial community of cut chrysanthemum, and together with plant pathogens, they contribute to the replanting problem of cut chrysanthemum.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary files, further inquiries can be directed to the corresponding authors.
Author Contributions
ZR, XZ, and LL conceived and designed the experiments. TW, KY, and QM performed most of the experiments. TW and QM analyzed the data. TW wrote the manuscript. XJ and YZ helped to handle the plants. DK, ZW, and RP helped perform the analysis with constructive discussion. All authors read and approved the submitted version.
Funding
This research was partially supported by the National Natural Science Foundation of China (32070004 and 31670006). This study was also supported by the Fundamental Research Funds for Central Non-profit Scientific Institution (1610132020009), Central Public-interest Scientific Institution Basal Research Fund (Y2021GH18), Key Laboratory of Microbial Resources Exploitation and Application of Gansu Province (2019cx-04), and the Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture and Rural Affairs, China.
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to the Enpapers company for their editing services for the article and to our colleagues who helped us in the experiment.
Footnotes
References
Badri, D. V., and Vivanco, J. M. (2009). Regulation and function of root exudates. Plant Cell Environ. 32, 666–681. doi: 10.1111/j.1365-3040.2008.01926.x
Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., and Vivanco, J. M. (2006). The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol. 57, 233–266. doi: 10.1146/annurev.arplant.57.032905.105159
Breuillin-sessoms, F., Venterea, R. T., Sadowsky, M. J., Coulter, J. A., Clough, T. J., and Wang, P. (2017). Soil Biology & Biochemistry Nitri fi cation gene ratio and free ammonia explain nitrite and nitrous oxide production in urea-amended soils. Soil Biol. Biochem. 111, 143–153. doi: 10.1016/j.soilbio.2017.04.007
Bryant, D. A., and Frigaard, N. U. (2006). Prokaryotic photosynthesis and phototrophy illuminated. Trends Microbiol. 14, 488–496. doi: 10.1016/j.tim.2006.09.001
Cai, Z., Li, S., Zhang, W., and Ma, J. (2015). Effects of the novel pyrimidynyloxybenzoic herbicide ZJ0273 on enzyme activities, microorganisms and its degradation in Chinese soils. Environ. Sci. Pollut. Res. Int. 22, 4425–4433. doi: 10.1007/s11356-014-3674-1
Carolina, T., Javier, A., Adrián, C., and Graciela, P. (2021). Cladosporium species causing “‘ Cladosporium rot “’ on “‘ Bosc “’ pear fruit in Argentina. Rev. Argent. Microbiol. 53, 75–77. doi: 10.1016/j.ram.2019.11.006
Cha, J. Y., Han, S., Hong, H. J., Cho, H., Kim, D., Kwon, Y., et al. (2016). Microbial and biochemical basis of a Fusarium wilt-suppressive soil. ISME J. 10, 119–129. doi: 10.1038/ismej.2015.95
Chaparro, J. M., Badri, D. V., and Vivanco, J. M. (2014). Rhizosphere microbiome assemblage is affected by plant development. ISME J. 8, 790–803. doi: 10.1038/ismej.2013.196
Chen, W., Teng, Y., Li, Z., Liu, W., Ren, W., Luo, Y., et al. (2018). Mechanisms by which organic fertilizer and effective microbes mitigate peanut continuous cropping yield constraints in a red soil of south China. Appl. Soil Ecol. 128, 23–34. doi: 10.1016/j.apsoil.2018.03.018
China Ministry of Agriculture (2016). Available online at: http://zdscxx.moa.gov.cn:8080/nyb/pc/search.jsp
Clemens, S. (2017). The Environmental Detection of Comammox Nitrospira. Vienna: University of Vienna. doi: 10.25365/thesis.45220
Corné, M. J., Pieterse, Christos Zamioudis, Roeland, L., Berendsen, David, M., et al. (2014). Induced systemic resistance by beneficial microbes. Ann. Rev. Phytopathol. 52, 347–375.
Crater, G. D. (1992). “Potted Chrysanthemums,” in Introduction to Floriculture, (Berlin: Elsevier), 249–287. doi: 10.1016/b978-0-12-437651-9.50015-4
De Weert, S., Vermeiren, H., Mulders, I. H. M., Kuiper, I., Hendrickx, N., Bloemberg, G. V., et al. (2002). Flagella-driven chemotaxis towards exudate components is an important trait for tomato root colonization by Pseudomonas fluorescens. Mol. Plant-Microbe Interact. 15, 1173–1180. doi: 10.1094/MPMI.2002.15.11.1173
Dick, R. P., and Kandeler, E. (2004). “Enzymes in Soils,” in Encyclopedia of Soils in the Environment, (Netherland: Elsevier Inc.), 448–456. doi: 10.1016/B0-12-348530-4/00146-6
Duineveld, B. M., Kowalchuk, G. A., Keijzer, A., Van Elsas, J. D., and Van Veen, J. A. (2001). Analysis of bacterial communities in the rhizosphere of chrysanthemum via denaturing gradient gel electrophoresis of PCR-amplified 16S rRNA as well as DNA fragments coding for 16S rRNA. Appl. Environ. Microbiol. 67, 172–178. doi: 10.1128/AEM.67.1.172-178.2001
Ecker, D. J., Sampath, R., Willett, P., Wyatt, J. R., Samant, V., Massire, C., et al. (2005). The Microbial Rosetta Stone Database: A compilation of global and emerging infectious microorganisms and bioterrorist threat agents. BMC Microbiol. 5:19. doi: 10.1186/1471-2180-5-19
Edwards, J., Johnson, C., Santos-medellín, C., Lurie, E., and Kumar, N. (2014). Structure, variation, and assembly of the root-associated microbiomes of rice. PNAS 112, E911–E920. doi: 10.1073/pnas.1414592112
Eichorst, S. A., Trojan, D., Roux, S., Herbold, C., Rattei, T., and Woebken, D. (2018). Genomic insights into the Acidobacteria reveal strategies for their success in terrestrial environments. Environ. Microbiol. 20, 1041–1063. doi: 10.1111/1462-2920.14043
Eida, A. A., Ziegler, M., Lafi, F. F., Michell, C. T., Voolstra, C. R., Hirt, H., et al. (2018). Desert plant bacteria reveal host influence and beneficial plant growth properties. PLoS One. 13:e0208223. doi: 10.1371/journal.pone.0208223
Einhellig, F. A. (1995). “Allelopathy: current status and future goals,” in Allelopathy: Organisms, Processes and Applications, eds K. M. M. Inderjit Dakshini and F. A. Einhellig (Washington, DC: American Chemical Society), 1–24.
Fang, D., Zhao, G., Xu, X., Zhang, Q., Shen, Q., Fang, Z., et al. (2018). Microbial community structures and functions of wastewater treatment systems in plateau and cold regions. Bioresour. Technol. 249, 684–693. doi: 10.1016/j.biortech.2017.10.063
Fawzy, O., Maillard, E., Vuilleumier, S., Millet, M., and Imfeld, G. (2015). Science of the Total Environment Degradation of chloroacetanilide herbicides and bacterial community composition in lab-scale wetlands. Sci. Total Environ. 520, 222–231. doi: 10.1016/j.scitotenv.2015.03.061
Fehér, T. (1993). “Watermelon,” in Genetic Improvement of Vegetable Crops, (Netherland: Elsevier), 295–311. doi: 10.1016/b978-0-08-040826-2.50025-4
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Fitzpatrick, C. R., Copeland, J., Wang, P. W., Guttman, D. S., Kotanen, P. M., and Johnson, M. T. J. (2018). Assembly and ecological function of the root microbiome across angiosperm plant species. Proc. Natl. Acad. Sci. U.S.A. 115, E1157–E1165. doi: 10.1073/pnas.1717617115
Fitzpatrick, C. R., Mustafa, Z., and Viliunas, J. (2019). Soil microbes alter plant fitness under competition and drought. J. Evol. Biol. 32, 438–450. doi: 10.1111/jeb.13426
Frank, K., Beegle, D., and Denning, J. (1998). Recommended Chemical Soil test Procedures for the North Central Region. Columbia: University of Missour.
Gao, M., Song, W., Zhou, Q., Ma, X., and Chen, X. (2013). Interactive effect of oxytetracycline and lead on soil enzymatic activity and microbial biomass. Environ. Toxicol. Pharmacol. 36, 667–674. doi: 10.1016/j.etap.2013.07.003
Gianfreda, L., Antonietta, M., Piotrowska, A., Palumbo, G., and Colombo, C. (2005). Soil enzyme activities as affected by anthropogenic alterations?: intensive agricultural practices and organic pollution. Sci. Total Environ. 341, 265–279. doi: 10.1016/j.scitotenv.2004.10.005
Giraldo, A., and Crous, P. W. (2019). Studies in Mycology. Stud. Mycol. 92, 227–286. doi: 10.1016/j.simyco.2018.10.005
Grayston, S. J., Wang, S., Campbell, C. D., and Edwards, A. C. (1998). Selective influence of plant species on microbial diversity in the rhizosphere. Soil Biol. Biochem. 30, 369–378. doi: 10.1016/S0038-0717(97)00124-7
Gross, D., and Parthier, B. (1994). Novel natural substances acting in plant growth regulation. J. Plant Growth Regul. 13, 93–114. doi: 10.1007/bf00210953
Han, S., Luo, X., Liao, H., Nie, H., Chen, W., and Huang, Q. (2017). Nitrospira are more sensitive than Nitrobacter to land management in acid, fertilized soils of a rapeseed-rice rotation field trial. Sci. Total Environ. 599–600, 135–144. doi: 10.1016/j.scitotenv.2017.04.086
Hou, Y. M., Zhang, X., Zhang, N. N., Naklumpa, W., Zhao, W. Y., Liang, X. F., et al. (2019). Genera Acremonium and Sarocladium Cause Brown Spot on Bagged Apple Fruit in China. Plant Dis. 103, 1889–1901. doi: 10.1094/PDIS-10-18-1794-RE
Janssen, P. H. (2006). Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 72, 1719–1728. doi: 10.1128/AEM.72.3.1719
Jeong, H., Jeong, D. E., Kim, S. H., Song, G. C., Park, S. Y., Ryu, C. M., et al. (2012). Draft genome sequence of the plant growth-promoting bacterium Bacillus siamensis KCTC 13613T. J. Bacteriol. 194, 4148–4149. doi: 10.1128/JB.00805-12
Johansson, L. H. (1988). A spectrophotometric method for determination of catalase activity in small tissue samples. Anal. Biochem. 174, 331–336. doi: 10.1016/0003-2697(88)90554-4
Kandeler, E., and Gerber, H. (1988). Short-term assay of soil urease activity using colorimetric determination of ammonium. Biol. Fert Soils 6, 68–72.
Kinkel, L. L., Schlatter, D. C., Bakker, M. G., and Arenz, B. E. (2012). Streptomyces competition and co-evolution in relation to plant disease suppression. Res. Microbiol. 163, 490–499. doi: 10.1016/j.resmic.2012.07.005
Kluitenberg, G. J., and Biggar, J. W. (1992). Canopy temperature as a measure of salinity stress on sorghum. Irrig. Sci. 13, 115–121. doi: 10.1007/BF00191053
Ko, W. H., Yang, C. H., Lin, M. J., Chen, C. Y., and Tsou, Y. J. (2011). Humicola phialophoroides sp. nov. from soil with potential for biological control of plant diseases. Bot. Stud. 52, 197–202.
Kudjordjie, E. N., Sapkota, R., Steffensen, S. K., Fomsgaard, I. S., and Nicolaisen, M. (2019). Maize synthesized benzoxazinoids affect the host associated microbiome. Microbiome 7:59. doi: 10.1186/s40168-019-0677-7
Langille, M. G. I., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Lau, J. A., and Lennon, J. T. (2012). Rapid responses of soil microorganisms improve plant fitness in novel environments. Proc. Natl. Acad. Sci. U.S.A. 109, 14058–14062. doi: 10.1073/pnas.1202319109
Lee, S., Lee, J., Jeong, Y., and Kim, J. (2009). Chemosphere In situ stabilization of cadmium-, lead-, and zinc-contaminated soil using various amendments. Chemosphere 77, 1069–1075. doi: 10.1016/j.chemosphere.2009.08.056
Li, H., Qiu, Y., Yao, T., Ma, Y., Zhang, H., Yang, X., et al. (2020). Evaluation of seven chemical pesticides by mixed microbial culture (PCS-1): degradation ability, microbial community, and Medicago sativa phytotoxicity. J. Hazard. Mater. 389:121834. doi: 10.1016/j.jhazmat.2019.121834
Li, K., DiLegge, M. J., Minas, I. S., Hamm, A., Manter, D., and Vivanco, J. M. (2019). Soil sterilization leads to re-colonization of a healthier rhizosphere microbiome. Rhizosphere 12:100176. doi: 10.1016/j.rhisph.2019.100176
Li, P., Chen, J., Li, Y., Zhang, K., and Wang, H. (2017). Possible mechanisms of control of Fusarium wilt of cut chrysanthemum by Phanerochaete chrysosporium in continuous cropping fields: A case study. Sci. Rep. 7:15994. doi: 10.1038/s41598-017-16125-7
Li, P., Zhang, F., Chen, S., Jiang, J., Wang, H., Su, J., et al. (2016). Genetic diversity, population structure and association analysis in cut chrysanthemum (Chrysanthemum morifolium Ramat.). Mol. Genet. Genomics 291, 1117–1125. doi: 10.1007/s00438-016-1166-3
Li, X. G., Ding, C. F., Zhang, T. L., and Wang, X. X. (2014). Fungal pathogen accumulation at the expense of plant-beneficial fungi as a consequence of consecutive peanut monoculturing. Soil Biol. Biochem. 72, 11–18. doi: 10.1016/j.soilbio.2014.01.019
Liu, L., Kloepper, J. W., and Tuzun, S. (1995). Induction of systemic resistance in cucumber by plant growth-promoting rhizobacteria: duration of protection and effect of host resistance on protection and root colonization. Phytopathology 85, 1064–1068. doi: 10.1094/phyto-85-1064
Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277. doi: 10.1126/science.aaf4507
Lu, T., Ke, M., Lavoie, M., Jin, Y., Fan, X., Zhang, Z., et al. (2018). Rhizosphere microorganisms can influence the timing of plant flowering. Microbiome 6:231. doi: 10.1186/s40168-018-0615-0
Ma, L., Geiser, D. M., Proctor, R. H., Rooney, A. P., Donnell, K. O., Trail, F., et al. (2013). Fusarium Pathogenomics. Annu. Rev. Microbiol. 67, 399–416. doi: 10.1146/annurev-micro-092412-155650
Marschner, P., Crowley, D., and Yang, C. H. (2004). Development of specific rhizosphere bacterial communities in relation to plant species, nutrition and soil type. Plant Soil 261, 199–208. doi: 10.1023/B:PLSO.0000035569.80747.c5
Mendes, R., Kruijt, M., De Bruijn, I., Dekkers, E., Van Der Voort, M., Schneider, J. H. M., et al. (2011). Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science 332, 1097–1100. doi: 10.1126/science.1203980
Morales-Cedeño, L. R., Orozco-Mosqueda, M., del, C., Loeza-Lara, P. D., Parra-Cota, F. I, de los Santos-Villalobos, S., et al. (2021). Plant growth-promoting bacterial endophytes as biocontrol agents of pre- and post-harvest diseases: Fundamentals, methods of application and future perspectives. Microbiol. Res. 242:126612. doi: 10.1016/j.micres.2020.126612
Mukhopadhya, I., Hansen, R. El-Omar, E. M., and Hold, G. L. (2012). IBD-what role do Proteobacteria play? Nat. Rev. Gastroenterol. Hepatol. 9, 219–230. doi: 10.1038/nrgastro.2012.14
Nannipieri, P., Giagnoni, L., Landi, L., and Renella, G. (2011). “Role of phosphatase enzymes in soil,” in Phosphorus in Action, Vol. 26, eds E. Bünemann, A. Oberson, and E. Frossard (Berlin: Springer), 215–243.
Narisawa, K., Tokumasu, S., and Hashiba, T. (1998). Suppression of clubroot formation in Chinese cabbage by the root endophytic fungus, Heteroconium chaetospira. Plant Pathol. 47:206. doi: 10.1046/j.1365-3059.1998.00225.x
Naylor, D., and Coleman-Derr, D. (2018). Drought stress and root-associated bacterial communities. Front. Plant Sci. 8:223. doi: 10.3389/fpls.2017.02223
Naylor, D., Dong, Z., Simmons, T., Hixson, K. K., Kim, Y., Zink, E. M., et al. (2018). Erratum: Drought delays development of the sorghum root microbiome and enriches for monoderm bacteria (Proceedings of the National Academy of Sciences of the United States of America (2018) 115 (E4284-E4293) Proc. Natl. Acad. Sci. U.S.A. 115:E4952. doi: 10.1073/pnas.1807275115
Nevens, F., and Reheul, D. (2001). Crop rotation versus monoculture; yield, N yield and ear fraction of silage maize at different levels of mineral N fertilization. Netherlands J. Agric. Sci. 49, 405–425. doi: 10.1016/S1573-5214(01)80026-9
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2015). FUNGuild?: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Olsen, S. R. (1954). Estimation of available phosphorus in soils by extraction with sodium bicarbonate. Miscell. Paper Instit. Agricult. Res. Samaru. 418, 18–19.
Ongena, M., and Jacques, P. (2008). Bacillus lipopeptides: versatile weapons for plant disease biocontrol. Trends Microbiol. 16, 115–125. doi: 10.1016/j.tim.2007.12.009
Penton, C. R., Gupta, V. V. S. R., Tiedje, J. M., Neate, S. M., Ophel-Keller, K., Gillings, M., et al. (2014). Fungal community structure in disease suppressive soils assessed by 28S LSU gene sequencing. PLoS One 9:e93893. doi: 10.1371/journal.pone.0093893
Qiu, M., Zhang, R., Xue, C., Zhang, S., Li, S., Zhang, N., et al. (2012). Application of bio-organic fertilizer can control Fusarium wilt of cucumber plants by regulating microbial community of rhizosphere soil. Biol. Fertil. Soils 48, 807–816. doi: 10.1007/s00374-012-0675-4
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 321, 341–361. doi: 10.1007/s11104-008-9568-6
Rachel, Backer, J., Stefan, Rokem, Gayathri, et al. (2018). Plant growth-promoting rhizobacteria: context, mechanisms of action, and roadmap to commercialization of biostimulants for sustainable agriculture. Front. Plant Sci. 9:1473. doi: 10.3389/fpls.2018.01473
Ragot, S. (2016). Soil phoD and phoX alkaline phosphatase gene diversity responds to multiple environmental factors. FEMS Microbiol. Ecol. 93:fiw212. doi: 10.1093/femsec/fiw212
Rana, M. A., Mahmood, R., and Ali, S. (2021). Soil urease inhibition by various plant extracts. PLoS One 16:e0258568. doi: 10.1371/journal.pone.0258568
Riaz, M., Arif, M. S., Ashraf, M. A., Mahmood, R., Yasmeen, T., Shakoor, M. B., et al. (2018). “A comprehensive review on rice responses and tolerance to salt stress,” in Advances in Rice Research for Abiotic Stress Tolerance, (Berlin: Elsevier), 133–158. doi: 10.1016/B978-0-12-814332-2.00007-1
Rice, E. L. (1995). Biological Control of Weeds and Plant Diseases. Norman, USA: University of Oklahomka Press.
Richardson, A. E., and Simpson, R. J. (2011). Soil microorganisms mediating phosphorus availability update on microbial phosphorus. Plant Physiol. 156, 989–996. doi: 10.1104/pp.111.175448
Santhanam, R., Weinhold, A., Goldberg, J., Oh, Y., and Baldwin, I. T. (2015). Native root-associated bacteria rescue a plant from a sudden-wilt disease that emerged during continuous cropping. Proc. Natl. Acad. Sci. U.S.A. 112, E5013–E5020. doi: 10.1073/pnas.1505765112
Scandinavica, A. A., Soil, S. B., Tan, X., Kong, L., Yan, H., Wang, Z., et al. (2014). Influence of soil factors on the soil enzyme inhibition by Cd. Acta Agric. Scand. Sect. B 64, 666–674. doi: 10.1080/09064710.2014.953985
Segata, N., Izard, J., Walron, L., Gevers, D., Miropolsky, L., Garrett, W., et al. (2011). Metagenomic Biomarker Discovery and Explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Singh, A., and Ghoshal, N. (2013). European Journal of Soil Biology Impact of herbicide and various soil amendments on soil enzymes activities in a tropical rainfed agroecosystem. Eur. J. Soil Biol. 54, 56–62. doi: 10.1016/j.ejsobi.2012.10.003
Singh, H. P., Batish, D. R., Kohli, R. K., Singh, H. P., Batish, D. R., and Kohli, R. K. (2010). Autotoxicity: Concept, Organisms, and Ecological Significance. Crit. Rev. Plant Sci. Autot. 2689, 757–772. doi: 10.1080/07352689991309478
Song, A., Zhao, S., Chen, S., Jiang, J., Chen, S., Li, H., et al. (2013). The abundance and diversity of soil fungi in continuously monocropped chrysanthemum. Sci. World J. 2013:632920. doi: 10.1155/2013/632920
Stępniewska, Z., Wolińska, A., and Ziomek, J. (2009). Response of soil catalase activity to chromium contamination. J. Environ. Sci. 21, 1142–1147. doi: 10.1016/s1001-0742(08)62394-3
Sun, Y., Sun, G., Xu, Y., Wang, L., Liang, X., and Lin, D. (2013). Geoderma Assessment of sepiolite for immobilization of cadmium-contaminated soils. Geoderma 19, 149–155. doi: 10.1016/j.geoderma.2012.07.012
Sun, Z., Yang, L. M., Ma, X. J., and Zhang, Y. M. (2014). Screening and identification of antagonistic bacteria on alternaria panax Whetz. J. Jilin Agric. Univ. 36, 276–281.
Teixeira, P. L., Id, S. S., Law, T. F., and Jose, P. (2019). The effects of soil phosphorus content on plant microbiota are driven by the plant phosphate starvation response. PLoS Biol. 17:e3000534. doi: 10.1371/journal.pbio.3000534
Tian, F., Hou, M., Qiu, Y., Zhang, T., and Yuan, Y. (2020). Salinity stress effects on transpiration and plant growth under different salinity soil levels based on thermal infrared remote (TIR) technique. Geoderma 357:113961. doi: 10.1016/j.geoderma.2019.113961
Timm, C. M., Carter, K. R., Aa Carrell, Jun, S. R., and Weston, D. J. (2018). Abiotic stresses shift belowground populus -associated bacteria toward a core stress microbiome. mSystems 3, e70–e17. doi: 10.1128/mSystems.00070-17
Tripathi, S., Chakraborty, A., and Chakrabarti, K. (2007). Enzyme activities and microbial biomass in coastal soils of India. Soil Biol. Biochem. 39, 2840–2848. doi: 10.1016/j.soilbio.2007.05.027
Trivedi, P., Delgado-Baquerizo, M., Trivedi, C., Hamonts, K., Anderson, I. C., and Singh, B. K. (2017). Keystone microbial taxa regulate the invasion of a fungal pathogen in agro-ecosystems. Soil Biol. Biochem. 111, 10–14. doi: 10.1016/j.soilbio.2017.03.013
Trivedi, P., Trivedi, C., Grinyer, J., Anderson, I. C., and Singh, B. K. (2016). Harnessing host-vector microbiome for sustainable plant disease management of phloem-limited bacteria. Front. Plant Sci. 7:1423. doi: 10.3389/fpls.2016.01423
Verbon, E. H., and Liberman, L. M. (2016). Beneficial Microbes Affect Endogenous Mechanisms Controlling Root Development. Trends Plant Sci. 21, 218–229. doi: 10.1016/j.tplants.2016.01.013
Wagner, M. R., Lundberg, D. S., Coleman-Derr, D., Tringe, S. G., Dangl, J. L., and Mitchell-Olds, T. (2014). Natural soil microbes alter flowering phenology. Ecol. Lett. 17, 717–726. doi: 10.1111/ele.12276.Natural
Wang, X., Fan, J., Xing, Y., Xu, G., Wang, H., Deng, J., et al. (2019). “The Effects of Mulch and Nitrogen Fertilizer on the Soil Environment of Crop Plants,” in Advances in Agronomy, (New York, NY: Academic Press Inc), 121–173. doi: 10.1016/bs.agron.2018.08.003
Wasaki, J., Sakaguchi, J., Yamamura, T., Ito, S., Osaki, M., and Kandeler, E. (2018). Soil Science and Plant Nutrition P and N deficiency change the relative abundance and function of rhizosphere microorganisms during cluster root development of white lupin (Lupinus albus L.). Soil Sci. Plant Nutr. 00, 1–11. doi: 10.1080/00380768.2018.1536847
Wei, X., Hu, Y., Razavi, B. S., Zhou, J., Shen, J., Nannipieri, P., et al. (2019). Rare taxa of alkaline phosphomonoesterase-harboring microorganisms mediate soil phosphorus mineralization. Soil Biol. Biochem. 131, 62–70. doi: 10.1016/j.soilbio.2018.12.025
Wu, L., Kobayashi, Y., Wasaki, J., and Koyama, H. (2018). Soil Science and Plant Nutrition Organic acid excretion from roots?: a plant mechanism for enhancing phosphorus acquisition, enhancing aluminum tolerance, and recruiting beneficial rhizobacteria. Soil Sci. Plant Nutr. 00, 1–8. doi: 10.1080/00380768.2018.1537093
Xiong, W., Li, R., Ren, Y., Liu, C., Zhao, Q., Wu, H., et al. (2017). Distinct roles for soil fungal and bacterial communities associated with the suppression of vanilla Fusarium wilt disease. Soil Biol. Biochem. 107, 198–207. doi: 10.1016/j.soilbio.2017.01.010
Xiong, W., Zhao, Q., Zhao, J., Xun, W., Li, R., Zhang, R., et al. (2015). Different Continuous Cropping Spans Significantly Affect Microbial Community Membership and Structure in a Vanilla-Grown Soil as Revealed by Deep Pyrosequencing. Microb. Ecol. 70, 209–218. doi: 10.1007/s00248-014-0516-0
Yang, Y., Ding, J., Chi, Y., and Yuan, J. (2020). Characterization of bacterial communities associated with the exotic and heavy metal tolerant wetland plant Spartina alterniflora. Sci. Rep. 10:17985. doi: 10.1038/s41598-020-75041-5
Zhang, J., Liu, Y. X., Zhang, N., Hu, B., Jin, T., Xu, H., et al. (2019). Nrt1.1b is associated with root microbiota composition and nitrogen use in field-grown rice. Nat. Biotechnol. 37, 676–684. doi: 10.1038/s41587-019-0104-4
Zhao, J., Yang, Y., Zhang, K., Jeong, J., Zeng, Z., and Zang, H. (2020). Does crop rotation yield more in China? A meta-analysis. F. Crop. Res. 245:107659. doi: 10.1016/j.fcr.2019.107659
Zhao, Q., Jian, S., Nunan, N., Maestre, F. T., Tedersoo, L., He, J., et al. (2017). Altered precipitation seasonality impacts the dominant fungal but rare bacterial taxa in subtropical forest soils. Biol. Fertil. Soils 53, 231–245. doi: 10.1007/s00374-016-1171-z
Zhao, Y. F., Wu, L. J., Yang, Y. J., and Li, X. H. (1993). Study on the pathogon of ginseng phytophthoya blight. Acta Phytopathol. Sin. 3, 19–20.
Keywords: cut chrysanthemum, rhizosphere, community structure, soil physical and chemical property decline, continuous cropping barrier
Citation: Wang T, Yang K, Ma Q, Jiang X, Zhou Y, Kong D, Wang Z, Parales RE, Li L, Zhao X and Ruan Z (2022) Rhizosphere Microbial Community Diversity and Function Analysis of Cut Chrysanthemum During Continuous Monocropping. Front. Microbiol. 13:801546. doi: 10.3389/fmicb.2022.801546
Received: 25 October 2021; Accepted: 19 January 2022;
Published: 16 March 2022.
Edited by:
Khondoker M. G. Dastogeer, Bangladesh Agricultural University, BangladeshReviewed by:
Deep Chandra Suyal, Eternal University, IndiaMakoto Kanasugi, Nicolaus Copernicus University in Toruń, Poland
Copyright © 2022 Wang, Yang, Ma, Jiang, Zhou, Kong, Wang, Parales, Li, Zhao and Ruan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lin Li, lilin@mail.hzau.edu.cn; Xin Zhao, zhaoxin01@caas.cn; Zhiyong Ruan, ruanzhiyong@caas.cn
†These authors have contributed equally to this work