Jing Xiong1†
Shuguang Peng2†Yongjun Liu2
Huaqun Yin1Lei Zhou3
Zhicheng Zhou2Ge Tan1Yabing Gu1Hetian Zhang1Jingyi Huang1
Delong Meng1*- 1School of Minerals Processing and Bioengineering, Central South University, Changsha, China
- 2Tobacco Research Institute of Hunan Province, Changsha, China
- 3Beijing Research Institute of Chemical Engineering and Metallurgy, Beijing, China
Continuous cropping barriers lead to huge agriculture production losses, and fumigation and biological agents are developed to alleviate the barriers. However, there is a lack of literature on the differences between strong chemical fumigant treatment and moderate biological agent treatment. In this study, we investigated those differences and attempted to establish the links between soil properties, rhizosphere microbial community, and plant performance in both fumigation- and bioagent-treated fields. The results showed that the fumigation had a stronger effect on both soil functional microbes, i.e., ammonia oxidizers and soil-borne bacterial pathogens, and therefore, led to a significant change in soil properties, higher fertilizer efficiency, lower disease infections, and improved plant growth, compared with untreated control fields. Biological treatment caused less changes to soil properties, rhizosphere bacterial community, and plant physiology. Correlation and modeling analyses revealed that the bioagent effect was mainly direct, whereas fumigation resulted in indirect effects on alleviating cropping barriers. A possible explanation would be the reconstruction of the soil microbial community by the fumigation process, which would subsequently lead to changes in soil characteristics and plant performance, resulting in the effective alleviation of continuous cropping barriers.
Introduction
Long-term monoculture of the same cash crop is a common land use pattern in large-scale agricultural and horticultural intensive production (Tilman et al., 2011). The problems caused by the pattern are referred to as continuous cropping barriers (or obstacles) (Li et al., 2019a). The continuous cropping barriers include many aspects, such as land degradation and plant soil-borne disease infection, which results in crop yield reduction (Zeng et al., 2020). Continuous cropping alters soil physiochemical properties and microbial community, causes soil salinization and acidification, accumulates harmful microbes, reduces fertilizer efficiency, and leads to severe soil-borne diseases (Trabelsi et al., 2012), resulting in yield reduction and huge economic losses in agriculture production (Gentry et al., 2013). Soil microbial communities are closely related to plant nutrient uptake, disease infection (Yang et al., 2017), and immune system and hence, play a crucial role in the formation and alleviation of continuous cropping obstacles (Tan et al., 2021). For example, native root-associated bacteria could inhibit the wilt disease during continuous cropping (Santhanam et al., 2015). Plants also recruit microbes to help them with nutrient acquisition, disease resistance, and stress tolerance (Philippot et al., 2013). However, soil microbiome and their ecological functions, such as biogeochemical process, are affected by the long-term continuous cropping. In a monoculture system, the diversity of the functional group, i.e., ammonia, oxidizers decreases with time, while ammonia oxidation genes (amo) are enriched, leading to the rapid oxidation of ammonia to nitrate (Sheng et al., 2013), which resulted in the nitrogen loss and low efficiency of nitrogen nutrients. Previous studies have also revealed that continuous cropping decreased soil bacterial diversity and relative abundance of plant-beneficial microorganisms and increased the abundance of pathogenic Ralstonia, which lead to bacterial wilt disease outbreak (She et al., 2017; Xiao et al., 2018). Considering the importance of plant and soil microbiomes in regulating plant productivity (Brooker et al., 2015), it is necessary to study the plant-associated microbial community, such as the rhizosphere microbial community, in continuous cropping field. The importance of soil and rhizosphere microbial community in alleviating continuous cropping barriers has been studied by Zhao et al. (2014). The changes in soil microbial community composition over time were related to changes in the soil nutrient availability during anaerobic soil disinfestation (ASD) treatment to alleviate continuous cropping barriers (Guo et al., 2018). Diverse soil microbial community tended to suppress the soil-borne disease (Yang et al., 2017), while plant growth-promoting rhizobacteria (PGPR) enhanced plant immunity to soil-borne pathogens and increased crop yield in continuous cropping land (Lugtenberg and Kamilova, 2009). Although numerous studies have focused on the role of microbial community in mediating the plant growth and plant health in continuous cropping field, the ecological roles of the rhizosphere microbial community during the alleviation of continuous cropping barriers using biological or chemical approaches have been rarely studied.
Many technologies have been developed to improve soil quality and alleviate the barriers caused by continuous cropping. These include physical, chemical, and biological technologies, such as intercropping (Miao et al., 2011), rotation (Luca et al., 2015; Yang et al., 2016), stereoscopic cultivation (Liao et al., 2018), application of plant beneficial microorganisms, e.g., PGPR (Singh et al., 2017; Goswami and Deka, 2020), and fumigation treatment (Buena et al., 2007; Li et al., 2016, 2019a). Climate and geography, as well as the type of cash crop, however, have limited the agronomic practices of intercropping, rotation, and stereoscopic cultivation (Tilman, 2020). Additionally, the economic returns of the associated plant have also limited the use of these strategies. Therefore, novel methods of efficient and high economic returns are necessary. As an alternative, soil fumigation has been globally employed wherever high-value agricultural crops are subjected to long-term continuous cropping cultures. The fumigants including both biological and chemical have been applied to combat a wide range of pests, including nematodes, fungi, and bacteria, in the soil that have damaging effects on plant growth and crop production (Yagi et al., 1993; Wang et al., 2014). However, fumigation may have adverse health effects on people if not properly administered. Although the fumigation treatment has been effective in alleviating continuous crop obstacle, its application has not been widely adopted (Li et al., 2016). Common and reliable fumigation treatments include chloropicrin (FM1) and dazomet (FM2) (Wang et al., 2013). In addition, soil fumigants, such as chloropicrin, have been shown to affect soil microbial community diversity and composition, change soil nutrient transformation, and affect soil fertility and agricultural productivity (Zhang D. et al., 2019). In this study, soil fumigants were used because of their effectiveness in eliminating continuous cropping barriers. Since microbial communities are among the major obstacles to continuous crop production, several efforts are being made through the use of (micro-)biological techniques (e.g., high-throughput sequencing and high-throughput isolation) to develop beneficial microorganisms to solve the problem of continuous cropping obstacles.
The effects of fumigation or biological agent treatment on continuous cropping problems have been extensively studied. It has been demonstrated that fumigation decreased soil microbial diversity and improved plant growth (Zhang et al., 2020). However, the inherent, complex link between different technologies, underground micro ecosystem, and above-ground plants during alleviating of continuous cropping obstacles is still unclear.
In this study, we employed fumigation and biological agent treatments to alleviate continuous cropping barriers. Rhizosphere bacterial community was investigated through high-throughput sequencing and was linked to plant physiology including phytochrome and growth conditions in different treatments. We compared the effects of biological agent treatments and that of the fumigation treatment, to construct a strong link between microbial community and plant performance. The result of this study will offer theoretical support for developing new approaches for solving the continuous cropping problems.
Materials and methods
Experimental design and soil property determination
All field experiments were carried out at the Huayuan Agricultural Science Park, located in Xiangxi, China (109°27′5″ E and 28°24’57″ N). Detailed field experimental design is shown in the Supplementary Material. A field suffering from serious continuous cropping barriers was employed to carry out fumigation and biological agent treatment. The field was divided into 36 equal plots (6 treatments × 6 replicates), and 168 seedlings were planted in each plot. Three fumigation treatments, namely, chloropicrin (FM1), dazomet (FM2), and untreated control (CK_FM), and three biological treatments, namely, two biological agents offered by Prof. Jian Ye from the Institute of Microbiology, CAS, and untreated control (CK_AG). Fumigation and biological treatments were performed before transplanting tobacco seedlings in April 2019. Soil-borne disease investigation, rhizosphere soil sampling, plant sampling, and plant physiological investigation were carried out in August 2019. One sample was collected from each plot by combining rhizosphere soils from 5 plants, resulting in a total of 36 samples. Detailed sampling methods and experimental design are shown in the Supplementary Material (“Materials and Methods” section). Soil pH was determined in deionized water extract (soil/water ratio of 5 g/25 ml) using a pH meter (portable ORP meter, BPH-220, Bell, Instrument Equipment Co., Ltd., China) Organic nitrogen was determined by Kjeldahl digestion; organic carbon was determined by potassium dichromate oxidation-ferrous sulfate titration. Ammonia N and nitrate N were extracted using 2 M KCl and determined with the FIAstar 5000 Analyzer. Available P was extracted with 0.5 M NaHCO3 and determined with a spectrophotometer. Available K was extracted using 1 M NH4OAc and determined with flame spectrophotometry.
Phytochrome analysis
The plant pigments were determined using the high-performance liquid chromatography (HPLC)-based method (Schwartz and Von Elbe, 1982). Briefly, a volume of 25 ml of 90% acetone solution was used to extract phytochrome from 0.2 g of liquid nitrogen ground plant leaf material. For extracting the phytochrome, leaf material in acetone solution was ultrasonically crushed for 20 min in an ice bath and then filtered using a 0.45 μm filter membrane. The filtrate was used to measure phytochrome using HPLC. The column used for HPLC was Waters Nova-Pak-C18 (3.9 mm × 150 mm, 4 μm). Analysis condition was set as follows: 30°C, 0.5 ml/min, 5 min for equilibrium, acetone as mobile phase A, and 80% acetonitrile water (V/V) solution as mobile phase A.
DNA extraction, PCR amplification, and sequencing
Soil total DNA was extracted using Soil FastDNA® Spin Kit (MP Biomedicals, Santa Ana, CA, United States) with 0.5 g of liquid nitrogen ground soil sample, following the manufacturer’s instructions. The V3–V4 region of 16S rDNA was amplified using primer pair 341F (5′-CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTACHVGGGTATCTAATCC-3′) with 12 bp barcode sequences. PCR amplification was performed on Bio-Rad S1000 (Bio-Rad Laboratory, CA, United States). The PCR reaction included 25 μl 2 × Premix Taq (Takara Biotechnology, Dalian Co., Ltd., China), 1 μl each primer (10 μM), and 3 μl DNA (20 ng/μl) template in a volume of 50 μl. The PCR condition was set as follows: 5 min at 94°C for initialization; 30 cycles of 30 s denaturation at 94°C, 30 s annealing at 52°C, and 30 s extension at 72°C; followed by 10 min final elongation at 72°C. After PCR amplification, the PCR products were extracted using 1% agarose gel electrophoresis, and the products were recovered and mixed in equidensity ratios according to the GeneTools Analysis Software (Version 4.03.05.0, SynGene). The mixed PCR product was purified using the E.Z.N.A. Gel Extraction Kit (Omega, United States). Sequencing libraries were generated using the NEBNext® Ultra™ II DNA Library Prep Kit for Illumina® (New England Biolabs, Ipswich, MA, United States) following the manufacturer’s instructions. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, United States). Finally, the library was sequenced on an Illumina Nova6000 platform, and 250 bp paired-end reads were generated (Guangdong Magigene Biotechnology Co., Ltd., Guangzhou, China).
Data analysis and statistical analysis
FASTQ format and library spliced raw data were obtained from the Guangdong Magigene Biotechnology Co., Ltd., for further analysis. Raw data were processed on the Galaxy pipeline developed by the Institute for Environmental Genomics, The University of Oklahoma1 to generate OTU table and represent sequences, as described in our previous studies (Li X. et al., 2017; Gu et al., 2019; Meng et al., 2019). Briefly, the forward and reverse reads were assembled with 10–200 bp overlapping using Flash (Version 1.0) (Magoc and Salzberg, 2011). Lower quality sequences (QC score < 20 and length < 250 bp) were trimmed using BTRIM (Version 1.0) (Kong, 2014), and sequences containing “N” were also removed, and only sequences with 400–440 in length were used for further bioinformatics analysis. Finally, chimeras were removed, and sequences with 97% identity were assigned to the same operational taxonomic unit (OTU) using UPARSE (version usearch v7.01001_i86linux64) (Edgar, 2013), and singletons that have no similar sequences were removed. Sequence number for each sample ranged from 33,953 to 45,904; therefore, we rarefied all samples to 33,953 by randomly choosing sequences. All downstream analyses were carried out using the rarefied OTU table. The raw data were submitted to NCBI SRA database, with the BioProject accession number PRJNA687637. OTU taxonomic assignment was performed by blasting representative sequences to 16S rRNA training set 18 with RDP classifier2. Bacterial diversity indexes calculation and beta diversity analyses were performed on the R statistic platform (Version 4.0.3) vegan package (Version 2.5–7). Linear discriminant analysis (LDA) effect size (LEfSe) analysis was carried out to determine the significant difference in the taxa between treatments. Analysis of variance (ANOVA) with the Tukey’s test was performed to identify the significant difference among group means, using aov and TukeyHSD functions on the R software. A p-value of less than 0.05 was considered significant. Correlation analysis was performed using the “corrplot” R package with “Pearson” method, and partial least squares path modeling (PLSPM) models were constructed with “plspm” R package.
Results and discussion
Fumigants and bioagents treatment effect on plant performance and soil properties
Continuous cropping was a leading cause of soil acidification and fertility loss, which resulted in loss in crop yield (Chen et al., 2013), plant physiology, and disease outbreaks (Li S. et al., 2017). Soil properties including pH, AP, AN, NN, and AK were significantly changed by fumigation treatment, particularly, the availability of potassium and phosphorate were significantly increased by both fumigation treatments (Table 1). The increase in plant available nutrients, i.e., potassium and phosphorate, therefore, promoted the growth of plants (Table 2). The plant height, length, surface area of the longest leaf, and number of leaves were significantly increased in fumigation treated fields. The application of bioagents did not cause a significant change in soil properties except for the soil pH, which was increased from 5.90 to 6.16 and 6.29 for AG1 and AG2, respectively. This result is consistent with previous studies, which reported on the increase of available nutrients such as N, P, and K by fumigation (Wang et al., 2014; Li et al., 2019a). Soil nutrient change, particularly changes in available elements (i.e., N, P, and K), indicated that the biogeochemical cycling of NPK was changed by fumigation treatments. For example, the ammonia oxidation tended to be inhibited in fumigation-treated soils because compared with the untreated control soil, there was higher ammonia N content, but lower nitrate N content in fumigation-treated soils. The nitrification ability was demonstrated to be enhanced in long-term continuous cropping and chemical fertilized soils, which is attributed to the enrichment of ammonia oxidizers (Lage et al., 2010; Sheng et al., 2013). Fumigation chemicals have been reported to be effective inhibitors of nitrification (Yan et al., 2017), as many chemicals can act as metal chelators and bind to the active site of the nitrification enzymes such as AMOB subunit (Beeckman et al., 2018). Therefore, one could hypothesize that fumigation may have a strong inhibition effect on soil nitrifiers. Although the biological agent treatment did not cause any change in soil nutrients (N, P, and K contents) or plant performance, it inhibited soil-borne pathogen infection (Figure 1). This may be ascribed to the increased pH as well as the significant change in microbial groups (e.g., beneficial and pathogenic microbial groups) with specific functions in biologically treated soils. However, despite no changes in the microbial community and diversity by the biological agents, the Aridibacter, Arthrobacter, Ralstonia, and Ramlibacter were enriched in AG1-treated field, while Solitalea was enriched in AG2-treated field. Additionally, the phytochrome, such as the chlorophyll a, was increased by the biological agent treatment (Figure 1). Phytochrome content is essential for crop yield and physiology since the rate of photosynthesis is determined by the phytochrome (Boccalandro et al., 2003). The phytochrome could also regulate the immunity of crops (Moreno and Ballare, 2014). Other than the chlorophyll a, the violaxanthin was also significantly increased by the biological agent treatment.
Table 1. Soil properties.
Table 2. Plant physiology properties.
Figure 1. Plant soil-borne disease infection, (A) infected plant numbers and (B) disease infection rate, and plant phytochrome content in leaves (C). Results are means and SD of six replicated plots; different letters above the columns indicate the differences are significant at p < 0.05 level.
Compared with the biological agent treatment, fumigation treatment had stronger disease inhibition effects. Furthermore, not only did the fumigation treatment increase pH, but it also increased nutrient availability in soil and stimulated crop growth. Previous studies have reported that nutrient availability and plant physiology enhances plant immunity to pathogen infections (Moreno and Ballare, 2014). However, the fumigation only caused a significant increase in chlorophyll a but not in violaxanthin. This result indicated that a direct way of alleviating continuous barriers exists through fumigation, i.e., through the direct change in soil properties. However, the use of biological agents, which alter some specific microbial groups and phytochrome, could be an indirect way of alleviating continuous barriers. Similar results were obtained by Singh et al. (2018), who reported an increase in the leaf chlorophyll content in combined fumigation and PGPR-treated fields (Singh et al., 2018). Increase in phytochrome would lead to strong photosynthetic rate and therefore, was responsible for the plant growth promotion (Boccalandro et al., 2003).
Rhizosphere bacterial community effect on alleviating continuous cropping
Considering that the agriculture management patterns (e.g., fertilization and irrigation) were the same for all plots, the changes in soil nutrients, plant performance, and disease infection caused by fumigation treatments could be attributed to the influence of soil microbial community (Li et al., 2019b; Akimbekov et al., 2020). The bacterial community was investigated by high-throughput sequencing of 16S rDNA amplicons. The bacterial diversity was significantly reduced by the application of fumigants (Figure 2), compared with the biological agents, which caused a significant increase in bacterial diversity. This could be explained by the fact that fumigation chemicals are often toxicants to organisms; therefore, microbes that were sensitive to the fumigants were stressed. A significant decrease in bacterial richness (including observed OTU number and Chao1 index) indicated a decline in some bacterial taxa from the soil after fumigation treatment. Ammonia oxidizers might be such a group that was sensitive to fumigation chemicals. This explains the increase in ammonia nitrogen and decrease in nitrate nitrogen following fumigation treatment. Fumigants are well-known poisonous substances used in killing insects, nematodes, and other animals or plants that cause damage to foods, seeds, or human dwelling. Soil fumigants, such as chloropicrin and dazomet, used in this research may kill soil organisms by predation or infection (Yan et al., 2017). It is widely accepted that the diverse soil microbial community would benefit above-ground crops against disease infection and promote plant growth (Brooker et al., 2015; Tao et al., 2017; Yang et al., 2017; Xiao et al., 2018; Mathes et al., 2020). However, this was in contrast to the results obtained in this study for fumigation treatment. This can be ascribed to the fumigation treatment disrupting the entire soil microbial ecosystem (Griffiths et al., 2000), and the crop benefiting from the reassembled rhizosphere microbiome (Nicola et al., 2017; Zhang et al., 2020).
Figure 2. Bacterial diversity index in soil with different treatments. AG1, microbial agent 1 treatment; AG2, microbial agent 2 treatment; CK_AG, untreated control in biological treatment experiment. FM1, CoCl2, cobalt dichloride fumigation treatment; FM2, dazomet fumigation treatment; CK_FM, untreated control in fumigation experiment. Different letters above the box indicated the differences between treatments were significant at p < 0.05 level.
Application of microbial agents in agriculture system has received widespread attention. Previous studies have shown that microbial flora was effective in improving crop yield (Gu et al., 2019; Wang et al., 2020) and also suppressing disease outbreak (Raaijmakers and Mazzola, 2012; Verbruggen et al., 2013; Qin et al., 2019; Wang et al., 2020). Although the biological agent treatments showed some inhibition effects on plant disease, they had less effect on soil bacterial diversity, indicating the disease inhibition effect by biological agents might be a result of other factors such as soil pH, bacterial interactions, and/or plant immunity system.
Similar to the changes in bacterial diversity, fumigation had a significant effect on the changes in rhizosphere microbial community composition as compared with the biological agent (Figures 3A–C). Taking together, the Venn diagram showed there were 2,260 core OTUs in all samples (Figure 3D). Compared with the untreated control rhizosphere soil (4,051 for CK_AG and 4,102 for CK_FM), unique OTUs were slightly increased in microbial agent-treated rhizosphere soil (4,239 for AG1 and 4,156 for AG2), whereas the fumigation treatments decreased in unique OTUs (3,655 for FM1 and 3,461 for FM2). The results further confirmed that the fumigation chemicals were toxic to microorganisms. It is worthy to note that, although individual biological agent (Figure 3A) did not cause a significant shift in bacterial community composition, when two agents were combined, they changed the bacterial community composition significantly compared with the untreated control (ADNOIS F = 1.83 and p = 0.027, based on Bray-Curtis distance). A significant change in microbial community composition and a slight increase in the number of unique OTUs after microbial agent treatment would be responsible for the decrease in soil-borne plant disease infection, as literature reports that the rhizosphere microbial community is closely related to plant health and diverse rhizosphere microbial community would suppress soil-borne diseases (Sanguin et al., 2009; Bakker et al., 2010). Proteobacteria was the most abundant and diverse phylum (averagely 38.63% of the bacterial community) among all 29 phyla in rhizosphere soil, as shown in the bubble plot (Figure 4). Gemmatimonas was the most abundant genus, which averagely accounted for 6.35% of the bacterial community, followed by Gp3 (5.08%), Gp1 (3.30%), Sphingobium (3.22%), Saccharibacteria_genera_incertae_sedis (3.03%), and Sphingomonas (2.96%). It was also discovered that typical nitrifier, Nitrospira, was an abundant genus in rhizosphere soil, and the bubble size indicated the abundance of Nitrospira decreased in FM1 and FM2 soil (see also in Supplementary Table 1). We carried out the LEfSe analyses to reveal the specific changes in bacterial taxa after microbial agent and fumigation treatments (Figure 5). Results showed a similar pattern as the change in bacteria diversity and composition. More taxa with significant change (LDA > 3.5) were detected in fumigation treatments than those detected in microbial agent treatments, indicating that fumigation chemicals had stronger effects on soil microbial community than the biological agents. Biological agent AG1 enriched four genera, namely, Aridibacter, Arthrobacter, Ralstonia, and Ramlibacter. The AG2 only enriched the genus Solitalea. Fumigation treatments had an intense effect on rhizosphere bacteria. FM1 treatment enriched Arthrobacter, Streptomyces, Chryseobacterium, Flavobacterium, and Taibaiella, while FM2 treatment enriched Micromonospora, Nocardioides, Microbispora, Nonomuraea, Actinomadura, Ohtaekwangia, and Mucilaginibacter. Ordinary rhizosphere microbes such as Gp1, Gp16, Gp2, Gp3, Gp4, Gp6, and Gp7 in the Phylum were found to be more abundant in the untreated control. Acidobacteria were more abundant in control soil (Lee et al., 2008). Due to the toxicity of fumigation chemicals to organisms, microbes that lack special strategies to live in extreme conditions die or are inhibited. Previous studies reported that members in the Phylum of Acidobacteria are often sensitive to environmental changes, hence, their inhibition by fumigation treatments in this study (Barns et al., 1999). However, it is interesting to note that microbes enriched in fumigation treatment naturally have strong stress-resistant ability, e.g., Streptomyces (Thompson et al., 1980), or special structures, such as endospore (e.g., Microbispora and Micromonospora), that aids them to withstand adverse conditions (Nicholson et al., 2000; Henkin, 2016). Fumigants have been reported to have a strong inhibition effect on soil nitrifiers (Yan et al., 2017); therefore, we selected typical nitrifiers from the OTU table and compared their abundance among fumigation treatments (Supplementary Table 1). The results showed that typical nitrifiers, such as Nitrospira and Nitrospirillum, decreased significantly in fumigation-treated soil (Supplementary Table 1), and this is consistent with previous reports, which showed the relative abundance of nitrifiers decreased significantly following fumigation treatment (Yan et al., 2017). This could be the reason for the less nitrate nitrogen but more ammonia nitrogen in the fumigation-treated soils.
Figure 3. Bacterial community composition in rhizosphere soil after treated with biological agents (A) or fumigation chemicals (B). (C) An overall look; (D) Venn diagram. AG1, microbial agent 1 treatment; AG2, microbial agent 2 treatment; CK_AG, untreated control in biological treatment experiment. FM1, CoCl2, cobalt dichloride fumigation treatment; FM2, dazomet fumigation treatment; CK_FM, untreated control in fumigation experiment.
Figure 4. Bacterial community structure at genus and phylum level as indicated by bubble plot. Bubble size represents the abundance (0–1) of each genus, and the color of the burble indicates the phylum. AG1, microbial agent 1 treatment; AG2, microbial agent 2 treatment; CK_AG, untreated control in biological treatment experiment. FM1, CoCl2, cobalt dichloride fumigation treatment; FM2, dazomet fumigation treatment; CK_FM, untreated control in fumigation experiment.
Figure 5. LEfSe analysis showing different taxa among bacterial agent treatments (A) or fumigation treatments (B). A taxon with the LDA value of larger than 3.2 was selected as different taxa.
The effects of bioagent and fumigation on continuous cropping barriers
Fumigation treatment has stronger effect, on soil and plant than biological agents. The correlation analysis (Figures 6A,B) showed a similar trend. Fumigation treatment had a significant effect on either soil properties or rhizosphere microbial community, which led to subsequent changes in other consortia. Therefore, soil properties, microbial community, and plant growth showed significant positive or negative correlations in this study. On the contrary, the consortia in biological treatments showed no correlations.
Figure 6. Correlation relationship and PLSPM model of soil physiochemical properties, microbial diversity, phytochrome contents, disease infection and plant growth. (A,C) Biological treatment and (B,D) fumigation treatment. ENV, soil physiochemical properties. *0.01 ≤ p < 0.05; **0.001 ≤ p < 0.01; ***p < 0.001.
Soil acidification has been identified as one of the main factors causing soil-borne disease outbreaks; therefore, it was not surprising that the soil pH showed a negative correlation with the plant soil-borne disease infection in both fumigation and bioagent-treated fields. Leaf chlorophyll a was negatively correlated with plant disease, and this could be attributed to either disease infection inhibiting the synthesis of phytochrome (Berger et al., 2007) or the synthesis of phytochrome enhancing plant immunity to pathogen infections (Kariola et al., 2005). However, the correlation between bacterial diversity (i.e., Richness and Shannon diversity) and plant growth differed between bioagent-treated fields and fumigation-treated soils. The negative correlation between bacterial diversity and plant growth is in sharp contrast with other studies, which found that plants benefited from the diverse rhizosphere microbial community (Trivedi et al., 2020). The discrepancy between this study’s findings and those of other studies might be explained by the reconstruction of the bacterial community, decreased diversity, change in functional group, and compositional biases following fumigation treatment. This further explains the unusual positive correlation between microbial diversity and disease infection in fumigation treatment. Consequently, to explore the effect of the most important factor on continuous cropping barriers, the PLSPM analyses was constructed to evaluate the direct and indirect effects between indicators and latent constructs (Figures 6C,D). Generally, it was found that bioagent-treated fields had more positive effects than fumigation-treated fields. In bioagent-treated fields, the microbial diversity had a strong direct effect on plant growth (0.6166). On the contrary, microbial diversity had weak direct effects on plant growth (−0.0061) in fumigation-treated fields. However, the effect of microbial diversity on soil properties was weak (0.1624) in bioagent-treated fields but strong (−0.4783) in fumigation-treated fields. The results indicated that the PGPR in bioagent may directly affect plant performance and, therefore, promote plant growth, while fumigation may shift the microbial community that promotes nutrient transformation and plant health and regulates plant growth (Zhang L. et al., 2019).
Agriculture development requires green and sustainable technologies. The need to develop eco-friendly and cost-effective technologies is, therefore, a pressing necessity. This study offers theoretical support for chemical and biological approaches to alleviating the continuous barriers. However, further studies are required to establish the link between agricultural management and soil health, as well as the functional profiles of the soil microbial community. Additionally, the succession pattern of soil microbial community and absolute quantification of soil microbial functional groups needs to be elucidated.
Conclusion
The effects of bioagent and fumigation on continuous cropping barriers differed in the following three aspects. (i) The effects of bioagent on soil and crop were moderate, as compared with fumigation treatment, which had significantly strong effects on both soil and crop. (ii) Fumigation enhanced soil nutrient availability, decreased rhizosphere microbial diversity, balanced microbial community compositional biases, suppressed soil-borne disease, promoted plant growth, and thus, alleviated continuous cropping obstacles. The bioagents had some effects on soil acidification and, therefore, could effectively suppress soil-borne diseases. (iii) By contrast, the bioagent had strong direct effects, whereas fumigation had an indirect effect on alleviating continuous cropping barriers, accounting for the fact that fumigation first reconstructed soil microbial community, followed by a shift in soil properties and plant performance, resulting in an effective alleviation of continuous cropping obstacles.
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 in the article/Supplementary Material.
Author contributions
JX, SP, ZZ, LZ, and HY conceived and designed the research. JX, YL, GT, HZ, DM, and YG conducted the experiments. SP, YL, and ZZ contributed to new reagents or analytical tools. YL, JH, DM, and YG analyzed the data. DM wrote the manuscript. All authors read and approved the manuscript.
Funding
This research was supported by the key project of Science and Technology of Hunan Branch of China National Tobacco Corporation (XX2022-2024Aa01, HN2021KJ05, 20-22A02, HN2020KJ02, and 202104).
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 thank the National Natural Science Foundation of China (NSFC) (No. 41807332) and Hunan International Scientific and Technological Cooperation Base of Environmental Microbiome and Application (No. 2018WK4019) for funding the work. We also like to thank the Jian Ye from Institute of Microbiology, Chinese Academy of Sciences, for offering the biological agents isolated from the soils suffering continuous cropping barriers, Zhengzhou Tobacco Research Institute for measuring the phytochrome contents, and Emmanuel Konadu Sarkodie for improving the quality of English.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.923405/full#supplementary-material
Footnotes
References
Akimbekov, N., Qiao, X. H., Digel, Y., Abdieva, G. Z. M., Ualieva, E. R. Z., and Zhubanova, A. Z. R. (2020). The effect of leonardite-derived amendments on soil microbiome structure and potato yield. Agric. Basel 10:147. doi: 10.3390/agriculture10050147
Bakker, M. G., Glover, J. D., Mai, J. G., and Kinkel, L. L. (2010). Plant community effects on the diversity and pathogen suppressive activity of soil streptomycetes. Appl. Soil Ecol. 46, 35–42. doi: 10.1016/j.apsoil.2010.06.003
Barns, S. M., Takala, S. L., and Kuske, C. R. (1999). Wide distribution and diversity of members of the bacterial kingdom acidobacterium in the environment. Appl. Environ. Microbiol. 65, 1731–1737. doi: 10.1128/aem.65.4.1731-1737.1999
Beeckman, F., Motte, H., and Beeckman, T. (2018). Nitrification in agricultural soils: impact, actors and mitigation. Curr. Opin. Biotechnol. 50, 166–173. doi: 10.1016/j.copbio.2018.01.014
Berger, S., Sinha, A. K., and Roitsch, T. (2007). Plant physiology meets phytopathology: plant primary metabolism and plant–pathogen interactions. J. Exp. Bot. 58, 4019–4026. doi: 10.1093/jxb/erm298
Boccalandro, H. E., Ploschuk, E. L., Yanovsky, M. J., Sanchez, R. A., Gatz, C., and Casal, J. J. (2003). Increased phytochrome B alleviates density effects on tuber yield of field potato crops. Plant Physiol. 133, 1539–1546. doi: 10.1104/pp.103.029579
Brooker, R. W., Bennett, A. E., Cong, W. F., Daniell, T. J., George, T. S., Hallett, P. D., et al. (2015). Improving intercropping: a synthesis of research in agronomy, plant physiology and ecology. New Phytol. 206, 107–117. doi: 10.1111/nph.13132
Buena, A. P., Garcia-Alvarez, A., Diez-Rojo, M. A., Ros, C., Fernandez, P., Lacasa, A., et al. (2007). Use of pepper crop residues for the control of root-knot nematodes. Bioresour. Technol. 98, 2846–2851. doi: 10.1016/j.biortech.2006.09.042
Chen, D., Lan, Z., Bai, X., Grace, J. B., and Bai, Y. (2013). Evidence that acidification-induced declines in plant diversity and productivity are mediated by changes in below-ground communities and soil properties in a semi-arid steppe. J. Ecol. 101, 1322–1334. doi: 10.1111/1365-2745.12119
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10:996. doi: 10.1038/NMETH.2604
Gentry, L. F., Ruffo, M. L., and Below, F. E. (2013). Identifying factors controlling the continuous corn yield penalty. Agron. J. 105, 295–303. doi: 10.2134/agronj2012.0246
Goswami, M., and Deka, S. (2020). Plant growth-promoting rhizobacteria-alleviators of abiotic stresses in soil: a review. Pedosphere 30, 40–61. doi: 10.1016/s1002-0160(19)60839-8
Griffiths, B. S., Ritz, K., Bardgett, R. D., Cook, R., Christensen, S., Ekelund, F., et al. (2000). Ecosystem response of pasture soil communities to fumigation-induced microbial diversity reductions: an examination of the biodiversity–ecosystem function relationship. Oikos 90, 279–294. doi: 10.1034/j.1600-0706.2000.900208.x
Gu, Y., Meng, D., Yang, S., Xiao, N., Li, Z., Liu, Z., et al. (2019). Invader-resident community similarity contribute to the invasion process and regulate biofertilizer effectiveness. J. Clean Prod. 241:118278. doi: 10.1016/j.jclepro.2019.118278
Guo, H., Zhao, X., Rosskopf, E. N., Di Gioia, F., Hong, J. C., and McNear, D. H. (2018). Impacts of anaerobic soil disinfestation and chemical fumigation on soil microbial communities in field tomato production system. Appl. Soil Ecol. 126, 165–173. doi: 10.1016/j.apsoil.2017.12.018
Henkin, T. M. (2016). Classic spotlight: bacterial endospore resistance, structure, and genetics. J. Bacteriol. 198, 1904–1904. doi: 10.1128/jb.00312-16
Kariola, T., Brader, G., Li, J., and Palva, E. T. (2005). Chlorophyllase 1, a damage control enzyme, affects the balance between defense pathways in plants. Plant Cell 17, 282–294. doi: 10.1105/tpc.104.025817
Kong, Y. (2014). Btrim: a fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics 98, 152–153. doi: 10.1016/j.ygeno.2011.05.009
Lage, M. D., Reed, H. E., Weihe, C., Crain, C. M., and Martiny, J. B. H. (2010). Nitrogen and phosphorus enrichment alter the composition of ammonia-oxidizing bacteria in salt marsh sediments. ISME J. 4, 933–944. doi: 10.1038/ismej.2010.10
Lee, S.-H., Ka, J.-O., and Cho, J.-C. (2008). Members of the phylum acidobacteria are dominant and metabolically active in rhizosphere soil. FEMS Microbiol. Lett. 285, 263–269. doi: 10.1111/j.1574-6968.2008.01232.x
Li, H., Yuan, G., Zhu, C., Zhao, T., Zhang, R., Wang, X., et al. (2019a). Soil fumigation with ammonium bicarbonate or metam sodium under high temperature alleviates continuous cropping-induced Fusarium wilt in watermelon. Sci. Horticult. 246, 979–986. doi: 10.1016/j.scienta.2018.11.090
Li, H., Zhang, Y. Y., Yang, S., Wang, Z. R., Feng, X., Liu, H. Y., et al. (2019b). Variations in soil bacterial taxonomic profiles and putative functions in response to straw incorporation combined with N fertilization during the maize growing season. Agric. Ecosyst. Environ. 283:106578. doi: 10.1016/j.agee.2019.106578
Li, R., Shen, Z., Sun, L., Zhang, R., Fu, L., Deng, X., et al. (2016). Novel soil fumigation method for suppressing cucumber Fusarium wilt disease associated with soil microflora alterations. Appl. Soil Ecol. 101, 28–36. doi: 10.1016/j.apsoil.2016.01.004
Li, S., Liu, Y., Wang, J., Yang, L., Zhang, S. T., Xu, C., et al. (2017). Soil acidification aggravates bacterial wilt in South China. Front. Microbiol. 8:703. doi: 10.3389/fmicb.2017.00703
Li, X., Meng, D., Li, J., Yin, H., Liu, H., Liu, X., et al. (2017). Response of soil microbial communities and microbial interactions to long-term heavy metal contamination. Environ. Pollut. 231, 908–917. doi: 10.1016/j.envpol.2017.08.057
Liao, P. R., Liu, P. F., Wang, Y. L., Huang, C. M., Lan, L., Yang, Y., et al. (2018). Stereoscopic cultivation of Panax nowginseng: a new approach to overcome the continuous cropping obstacle. Ind. Crops Prod. 126, 38–47. doi: 10.1016/j.indcrop.2018.09.042
Luca, C., Pilu, R., Tambone, F., Scaglia, B., and Adani, F. (2015). New energy crop giant cane (Arundo donax L.) can substitute traditional energy crops increasing biogas yield and reducing costs. Bioresour. Technol. 191, 197–204. doi: 10.1016/j.biortech.2015.05.015
Lugtenberg, B., and Kamilova, F. (2009). Plant-growth-promoting rhizobacteria. Annu. Rev. Microbiol. 63, 541–556. doi: 10.1146/annurev.micro.62.081307.162918
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Mathes, F., Murugaraj, P., Bougoure, J., Pham, V. T. H., Truong, V. K., Seufert, M., et al. (2020). Engineering rhizobacterial community resilience with mannose nanofibril hydrogels towards maintaining grain production under drying climate stress. Soil Biol. Biochem. 142:107715. doi: 10.1016/j.soilbio.2020.107715
Meng, D., Li, J., Liu, T., Liu, Y., Yan, M., Hu, J., et al. (2019). Effects of redox potential on soil cadmium solubility: insight into microbial community. J. Environ. Sci. China 75, 224–232. doi: 10.1016/j.jes.2018.03.032
Miao, Y. X., Stewart, B. A., and Zhang, F. S. (2011). Long-term experiments for sustainable nutrient management in China. A review. Agron. Sustain. Dev. 31, 397–414. doi: 10.1051/agro/2010034
Moreno, J. E., and Ballare, C. L. (2014). Phytochrome regulation of plant immunity in vegetation canopies. J. Chem. Ecol. 40, 848–857. doi: 10.1007/s10886-014-0471-8
Nicholson, W. L., Munakata, N., Horneck, G., Melosh, H. J., and Setlow, P. (2000). Resistance of Bacillus endospores to extreme terrestrial and extraterrestrial environments. Microbiol. Mol. Biol. Rev. 64, 548–572. doi: 10.1128/mmbr.64.3.548-572.2000
Nicola, L., Turco, E., Albanese, D., Donati, C., Thalheimer, M., Pindo, M., et al. (2017). Fumigation with dazomet modifies soil microbiota in apple orchards affected by replant disease. Appl. Soil Ecol. 113, 71–79. doi: 10.1016/j.apsoil.2017.02.002
Philippot, L., Raaijmakers, J. M., Lemanceau, P., and van der Putten, W. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109
Qin, C., Tao, J., Liu, T., Liu, Y., Xiao, N., Li, T., et al. (2019). Responses of phyllosphere microbiota and plant health to application of two different biocontrol agents. Amb. Express 9:42. doi: 10.1186/s13568-019-0765-x
Raaijmakers, J. M., and Mazzola, M. (2012). “Diversity and natural functions of antibiotics produced by beneficial and plant pathogenic bacteria,” in Annual Review of Phytopathology, Vol. 50, eds N. K. VanAlfen, J. E. Leach, and S. Lindow (Palo Alto, CA: Annual Reviews Inc), 403–424.
Sanguin, H., Sarniguet, A., Gazengel, K., Moenne-Loccoz, Y., and Grundmann, G. L. (2009). Rhizosphere bacterial communities associated with disease suppressiveness stages of take-all decline in wheat monoculture. New Phytol. 184, 694–707. doi: 10.1111/j.1469-8137.2009.03010.x
Santhanam, R., Luu, V. T., 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
Schwartz, S. J., and Von Elbe, J. H. (1982). High performance liquid chromatography of plant pigments – a review. J. Liq. Chromatogr. 5(Suppl. 1) 43–73. doi: 10.1080/01483918208082215
She, S., Niu, J., Zhang, C., Xiao, Y., Chen, W., Dai, L., et al. (2017). Significant relationship between soil bacterial community structure and incidence of bacterial wilt disease under continuous cropping system. Arch. Microbiol. 199, 267–275. doi: 10.1007/s00203-016-1301-x
Sheng, R., Meng, D., Wu, M., Di, H., Qin, H., and Wei, W. (2013). Effect of agricultural land use change on community composition of bacteria and ammonia oxidizers. J. Soils Sediments 13, 1246–1256. doi: 10.1007/s11368-013-0713-3
Singh, N., Sharma, D. P., and Thakur, K. K. (2018). Effect of rootstocks and soil management on growth and physiological parameters in new plantations of apple under replant conditions. Indian J. Horticult. 75, 392–398. doi: 10.5958/0974-0112.2018.00067.1
Singh, V. K., Singh, A. K., and Kumar, A. (2017). Disease management of tomato through PGPB: current trends and future perspective. 3 Biotech. 7:255. doi: 10.1007/s13205-017-0896-1
Tan, G., Liu, Y., Peng, S., Yin, H., Meng, D., Tao, J., et al. (2021). Soil potentials to resist continuous cropping obstacle: three field cases. Environ. Res. 200:111319. doi: 10.1016/j.envres.2021.111319
Tao, J., Liu, X., Liang, Y., Niu, J., Xiao, Y., Gu, Y., et al. (2017). Maize growth responses to soil microbes and soil properties after fertilization with different green manures. Appl. Microbiol. Biotechnol. 101, 1289–1299. doi: 10.1007/s00253-016-7938-1
Thompson, C. J., Ward, J. M., and Hopwood, D. A. (1980). DNA cloning in Streptomyces: resistance genes from antibiotic-producing species. Nature 286, 525–527. doi: 10.1038/286525a0
Tilman, D. (2020). Benefits of intensive agricultural intercropping. Nat. Plants 6, 604–605. doi: 10.1038/s41477-020-0677-4
Tilman, D., Balzer, C., Hill, J., and Befort, B. L. (2011). Global food demand and the sustainable intensification of agriculture. Proc. Natl. Acad. Sci. U.S.A. 108, 20260–20264. doi: 10.1073/pnas.1116437108
Trabelsi, D., Ben Ammar, H., Mengoni, A., and Mhamdi, R. (2012). Appraisal of the crop-rotation effect of rhizobial inoculation on potato cropping systems in relation to soil bacterial communities. Soil Biol. Biochem. 54, 1–6. doi: 10.1016/j.soilbio.2012.05.013
Trivedi, P., Leach, J. E., Tringe, S. G., Sa, T., and Singh, B. K. (2020). Plant-microbiome interactions: from community assembly to plant health. Nat. Rev. Microbiol. 18, 607–621. doi: 10.1038/s41579-020-0412-1
Verbruggen, E., van der Heijden, M. G. A., Rillig, M. C., and Kiers, E. T. (2013). Mycorrhizal fungal establishment in agricultural soils: factors determining inoculation success. New Phytol. 197, 1104–1109. doi: 10.1111/j.1469-8137.2012.04348.x
Wang, C. Q., Zhao, D. Y., Qi, G. Z., Mao, Z. Q., Hu, X. N., Du, B. H., et al. (2020). Effects of Bacillus velezensis FKM10 for promoting the growth of Malus hupehensis Rehd. and inhibiting Fusarium verticillioides. Front. Microbiol. 10:2889. doi: 10.3389/fmicb.2019.02889
Wang, Q.-X., Yan, D.-D., Mao, L.-G., Ma, T.-T., Liu, P.-F., Wu, Z.-F., et al. (2013). Efficacy of 1,3-dichloropropene plus chloropicrin gelatin capsule formulation for the control of soilborne pests. Crop Prot. 48, 24–28. doi: 10.1016/j.cropro.2013.02.002
Wang, Q., Ma, Y., Yang, H., and Chang, Z. (2014). Effect of biofumigation and chemical fumigation on soil microbial community structure and control of pepper phytophthora blight. World J. Microbiol. Biotechnol. 30, 507–518. doi: 10.1007/s11274-013-1462-6
Xiao, Y., Liu, X., Meng, D., Tao, J., Gu, Y., Yin, H., et al. (2018). The role of soil bacterial community during winter fallow period in the incidence of tobacco bacterial wilt disease. Appl. Microbiol. Biotechnol. 102, 2399–2412. doi: 10.1007/s00253-018-8757-3
Yagi, K., Williams, J., Wang, N. Y., and Cicerone, R. J. (1993). Agricultural soil fumigation as a source of atmospheric methyl bromide. Proc. Natl. Acad. Sci. U.S.A. 90, 8420–8423. doi: 10.1073/pnas.90.18.8420
Yan, D., Wang, Q., Li, Y., Ouyang, C., Guo, M., and Cao, A. (2017). Analysis of the inhibitory effects of chloropicrin fumigation on nitrification in various soil types. Chemosphere 175, 459–464. doi: 10.1016/j.chemosphere.2017.02.075
Yang, H., Li, J., Xiao, Y., Gu, Y., Liu, H., Liang, Y., et al. (2017). An integrated insight into the relationship between soil microbial community and tobacco bacterial wilt disease. Front. Microbiol. 8:2179. doi: 10.3389/fmicb.2017.02179
Yang, R., Mo, Y., Liu, C., Wang, Y., Ma, J., Zhang, Y., et al. (2016). The effects of cattle manure and garlic rotation on soil under continuous cropping of watermelon (Citrullus lanatus L.). PLoS One 11:e0156515. doi: 10.1371/journal.pone.0156515
Zeng, J., Liu, J., Lu, C., Ou, X., Luo, K., Li, C., et al. (2020). Intercropping with turmeric or ginger reduce the continuous cropping obstacles that affect Pogostemon cablin (Patchouli). Front. Microbiol. 11:579719. doi: 10.3389/fmicb.2020.579719
Zhang, D. Q., Yan, D. D., Cheng, H. Y., Fang, W. S., Huang, B., Wang, X. L., et al. (2020). Effects of multi-year biofumigation on soil bacterial and fungal communities and strawberry yield. Environ. Pollut. 256:113415. doi: 10.1016/j.envpol.2019.113415
Zhang, D., Yan, D. D., Fang, W. S., Huang, B., Wang, X. L., Wang, X. N., et al. (2019). Chloropicrin alternated with biofumigation increases crop yield and modifies soil bacterial and fungal communities in strawberry production. Sci. Total Environ. 675, 615–622. doi: 10.1016/j.scitotenv.2019.04.222
Zhang, L., Dong, H., Zhu, Y., Zhang, J., Zeng, G., Yuan, Y., et al. (2019). Evolutions of different microbial populations and the relationships with matrix properties during agricultural waste composting with amendment of iron (hydr)oxide nanoparticles. Bioresour. Technol. 289:121697. doi: 10.1016/j.biortech.2019.121697
Keywords: continuous cropping barriers, rhizosphere bacterial community, fumigation, bioagents treatment, plant performance
Citation: Xiong J, Peng S, Liu Y, Yin H, Zhou L, Zhou Z, Tan G, Gu Y, Zhang H, Huang J and Meng D (2022) Soil properties, rhizosphere bacterial community, and plant performance respond differently to fumigation and bioagent treatment in continuous cropping fields. Front. Microbiol. 13:923405. doi: 10.3389/fmicb.2022.923405
Received: 19 April 2022; Accepted: 27 June 2022;
Published: 22 July 2022.
Edited by:
Yurong Liu, Huazhong Agricultural University, ChinaReviewed by:
Akansha Jain, Banaras Hindu University, IndiaDiptaraj Sangramsing Chaudhari, University of South Florida, United States
Copyright © 2022 Xiong, Peng, Liu, Yin, Zhou, Zhou, Tan, Gu, Zhang, Huang and Meng. 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: Delong Meng, delong.meng@csu.edu.cn
†These authors have contributed equally to this work