- Research article
- Open Access
Organic mulching positively regulates the soil microbial communities and ecosystem functions in tea plantation
BMC Microbiology volume 20, Article number: 103 (2020)
Different mulches have variable effects on soil physicochemical characteristics, bacterial and fungal communities and ecosystem functions. However, the information about soil microbial diversity, community structure and ecosystem function in tea plantation under different mulching patterns was limited. In this study, we investigated bacterial and fungal communities of tea plantation soils under polyethylene film and peanut hull mulching using high-throughput 16S rRNA and ITS rDNA gene Illumina sequencing.
The results showed that the dominant bacterial phyla were Proteobacteria, Actinobacteria, Acidobacteria and Chloroflexi, and the dominant fungal phyla were Ascomycota, Mortierellomycota and Basidiomycota in all samples, but different mulching patterns affected the distribution of microbial communities. At the phylum level, the relative abundance of Nitrospirae in peanut hull mulching soils (3.24%) was significantly higher than that in polyethylene film mulching soils (1.21%) in bacterial communities, and the relative abundances of Mortierellomycota and Basidiomycota in peanut hull mulching soils (33.72, 21.93%) was significantly higher than that in polyethylene film mulching soils (14.88, 6.53%) in fungal communities. Peanut hull mulching increased the diversity of fungal communities in 0–20 cm soils and the diversity of bacterial communities in 20–40 cm soils. At the microbial functional level, there was an enrichment of bacterial functional features, including amino acid transport and metabolism and energy production and conversion, and there was an enrichment of fungal functional features, including undefined saprotrophs, plant pathogens and soils aprotrophs.
Unique distributions of bacterial and fungal communities were observed in soils under organic mulching. Thus, we believe that the organic mulching has a positive regulatory effect on the soil bacterial and fungal communities and ecosystem functions, and so, is more suitable for tea plantation.
Mulching is an important measure in agricultural production, which is mainly used for soil improvement and environmental protection. According to the type of materials, it could be categorized as inorganic and organic mulching. Inorganic mulching was widely used as a low-cost and water-saving measure in agriculture areas that were susceptible to drought. Although it could improve soil moisture, prevent soil nutrient loss and control crop pests and diseases [1, 2], it could change the soil biological characteristics and negatively impact on soil quality and sustainability, and even cause soil alkalization, resulting in injuries to plants . However, organic mulching is different, it is mainly derived from plant residues, which are proven to be better for soil health. The application of organic mulch on soils could not only inhibit weed germination, but also increase soil fertility and provide mineral elements for plants [4,5,6]. In addition, it also could increase the biodiversity of soil microecosystem .
Soil microorganism played important roles in nutrient cycling and structure maintain in the agro-ecosystem [8, 9]. Soil microbial diversity was involved in the circulation of matter and energy in the soil ecosystem . The changes of soil physicochemical properties during mulching process can drive the changes of the microbial community , which affect the microbial community structure and metabolic function [12, 13]. Soil microbial diversity and biomass were sensitive to the changes in the soil physicochemical properties and managements, and were useful for predicting the changes in the soil functions [14,15,16,17]. Therefore, investigating the diversity and composition of microbial communities in soils under mulching would be conducive to take effective measures to improve the fertility and productivity of soils and ensure the sustainable development of soil ecosystems.
Tea plant (Camellia sinensis L.), an evergreen leafy plant, prefers acidic soil according to the genetic characteristics, thus, the soil environment of tea plantation is different from that of other crops. Several studies reported that mulching could increase the richness and diversity of microbial community in other crop soils [18, 19]. However, to our knowledge, the information about soil microbial diversity and composition in tea plantation with inorganic or organic mulching was limited. In the present study, we performed one field trial to research the bacterial and fungal diversity, community composition and functional annotation of soils under two different materials (inorganic, polyethylene film; organic, peanut hull) on two depths of soils (0-20 cm depth; 20-40 cm depth). According to Linear discriminant analysis Effect Size (LEfSe) tool, soil bacterial and fungal communities with significant differences under different materials were identified. According to the spearman correlation, the relationships between soil bacterial, fungal community structure and soil property were discussed. This study not only compared the microbial effects of different materials, but also provided theoretical support for the rational application of different mulching materials in tea plantation. The results are of great significance not only for improving soil quality, but also for improving fertilizer efficiency.
The diversity of soil bacterial and fungal community
To ensure the quality of sequencing of all samples, we optimized the original sequences using the quality filter and chimera check. A total of 576,090 high-quality reads of bacteria and 886,572 high-quality reads of fungi were remained in the dataset with the average length of 418 bp and 234 bp, respectively. Rarefaction curve analysis showed that each curve was close to flat finally, indicating that the sample size of the sequencing was sufficient and the sequencing data were reasonable and acceptable (Additional file 1: Figure S1). Through clustering operations, the optimized sequences were classified into operational taxonomic units (OTUs) according to their similarity. With a 3% dissimilarity threshold, the sequences were classified into 6809 and 1684 OTUs in bacterial and fungal communities using the Ribosomal Database Project (RDP) classifier. In bacterial communities, the Venn diagram showed that the numbers of OTUs in P1 and P2 were more than that in F1 and F2 (Fig. 1a-b), but in fungal communities, the numbers of OTUs were not significant differences in three treatments (Fig. 1c-d).
To quantify the diversity and richness of microbial community of soils among three treatments, we calculated the indices of alpha diversity such as Chao1, ACE, Simpson and Shannon within a single microbial ecosystem (Table 1). The coverage indexes from 18 soil samples were greater than 0.96, indicating that the sequencing capacity was acceptable. Variance analysis showed that the different mulching materials and different soil depths affected the bacterial and fungal communities. In 0–20 cm soil, the richness of P1 was slightly higher than that in CK1 and F1 in bacterial communities, while there was no significant difference. However, in fungal communities, the richness of P1 was significantly lower than that in CK1 and F1. In 20–40 cm soil, the richness of P2 was significantly higher than that in CK2 and F2 in bacterial communities. Thus, peanut hull mulching greatly affected the diversity of fungal communities in 0-20 cm soil and bacterial communities in 20-40 cm soils.
To dig into the effects of different mulching materials on the distribution of soil microbial communities, we made the principal co-ordinates analysis (PCoA) based on the bray–curtis distance. The first principal coordination axis accounted for 37.45 and 25.04% in bacterial and fungal communities, respectively. The second principal coordination axis explained 13.53 and 15.74% in bacterial and fungal communities, respectively. In bacterial communities, there were significant differences between 0 and 20 cm soils and 20–40 cm soils: CK2, F2 and P2 were on the positive side of the x-axis, and CK1, F1and P1 were on the negative side. Furthermore, there were significant differences between the mulching plastic film and peanut shell on the topsoil: P1 and P2 were on the positive side of the y-axis, and CK1, CK2, F1 and F2 were on the negative side (Fig. 2a). Therefore, the distribution of bacterial communities was mainly affected by soil depths and mulching materials, especially organic mulches. In fungal communities, P1 and P2 clustered more closely together on the positive side of x-axis and CK1, CK2, F1 and F2 were on the negative side, indicating that the organic mulch could affect the distribution of fungal communities in soils (Fig. 2b).
The composition of soil bacterial and fungal community
To dissect the composition of soil microbial community, we aligned the top30 OTUs with the Silva 119 and Unite database at the phylum and genus levels. The dominant bacterial phyla were Proteobacteria, Actinobacteria, Acidobacteria and Chloroflexi (Fig. 3a), and the dominant fungal phyla were Ascomycota, Mortierellomycota and Basidiomycota (Fig. 3b). The dominant bacterial genus were norank_o__Gaiellales, Arthrobacter, norank_c__Subgroup_6, norank_c__KD4–96 and Sphingomonas (Fig. 3c), and the dominant fungal genus were Mortierella, Solicoccozyma, Fusarium and Clonostachys (Fig. 3d). In addition, there were significant differences in bacterial and fungal communities in soils under polyethylene film and peanut hull mulching (Fig. 4). At the phylum level, Nitrospirae presented significant differences in bacterial communities, and Mortierellomycota, Ascomycota and Basidiomycota presented significant differences in fungal communities. Therein, the relative abundance of Nitrospirae in P1 (3.24%) was significantly higher than that in F1 (1.21%) and CK1 (1.05%) of bacterial communities, and the relative abundances of Mortierellomycota (33.72%) and Basidiomycota (21.93%) in P1 was significantly higher than that in F1 and CK1of fungal communities. At the genus level, Sphingomonas, Nitrospira and Gemmayimonas presented significant differences in bacterial communities, Moretierella and solicoccozyma presented significant differences in fungal communities. Therein, the relative abundances of Nitrospira (3.24%) was significantly higher than that in F1 (1.21%) and CK1 (1.05%) of bacterial communities, and the relative abundances of Moretierella (33.66%) and solicoccozyma (10.60%) in P1 was significantly higher than that in F1 and CK1 of fungal communities.
To identify the bacterial and fungal communities with significant differences in soils among three treatments, we further analyzed microbial community in 0-20 cm soil from phylum to genus level using LEfSe tool (Additional file 2: Figure S2; Additional file 3: Figure S3). Eight groups of bacteria and two groups of fungi in P1 were significantly enriched, of which, bacteria included Gammaproteobacteria (class and the genus IS_44), Xanthobacteraceae (family), Parcubacteria (class), Amb_16S_1323 (the family and its genus), Candidatus_Magasanikbacteria (the order to genus), WS6_Dojkabacteria (the class to genus), UA_11 (the order to genus) and Paludibaculum (genus), and fungi included Piskurozymaceae (family and the genus Solicoccozyma) and Leucosporidiales (the order to genus). Three groups of bacteria and two groups of fungi in F1 were significantly enriched, of which, bacteria included Moraxellaceae (family), Solirubrobacteraceae (family) and Caulobacteraceae (family), and fungi included Acremonium (genus) and Emericellopsis (genus). The results showed that bacterial and fungal communities with significant differences in soils under different mulching were the noteworthy microbial communities in tea plantation.
Relationships between the preponderant phyla of microbial communities and the characteristics of soils
Firstly, to evaluate the physicochemical characteristics of soils under different mulching patterns, we analyzed the contents of nutrients in 0–20 cm and 20–40 cm depths soils (Table 2). In the topsoil, the content of TN in P1 was observably higher than that in CK1 and F1, and the soil moisture in P1 was significantly higher than that in CK1 and F1. While the contents of OM, AN and AP were slightly higher than those in CK1 and F1, but there was no significant difference. The pH value was F1 > P1 > CK1. In 20-40 cm soil, the contents of TN and TP in CK2 were observably higher than that in F2 and P2, and the soil moisture in P2 was significantly higher than that in CK2 and F2.
Secondly, to explore the influences of soil environmental factors on soil bacterial and fungal communities, we analyzed the relationships between microbial preponderant phyla and main soil properties using spearman correlation heatmap (Fig. 5). In bacterial communities, the soil pH was positively related with the abundance of Acidobacteria and Elusimicrobia. The soil moisture was positively related with the abundances of Nitrospirae and Rokubacteria, while it was negatively related with Cyanobacteria and Actinobacteria. The soil OM content was positively related with Bacteroidetes, while it was negatively related with the abundances of Firmicutes. The abundances of Firmicutes was positively interrelated with the contents of soil OM, TN, AN and AP. In fungal communities, the soil moisture was positively correlated to the abundances of Mortierellomycota, while it was negatively correlated to Chytridiomycota and Ascomycota. The soil OM was negatively correlated to Ascomycota. The AP and AK contents were positively correlated to Basidiomycota, while they were negatively correlated to the abundances of Ascomycota. The soil AN content was positively correlated to Mortierellomycota. The results showed that soil physicochemical characteristics had great influences on the structure of microbial communities.
Functional prediction of soil microbial communities
To better study the micro-ecological functions of the soil microorganisms in tea plantation under organic and inorganic mulching patterns, we analyzed the bacterial communities in the Cluster of Orthologous Groups (COG) database using PICRUSt (Fig. 6a). The results showed that the bacterial functional features were related to different treatments. It could be observed the metabolic functions were enriched in soil samples and the bacterial metabolism tended to be vigorous. These functional features included: General function prediction; Amino acid transport and metabolism; Energy production and conversion; Signal transduction mechanisms; Transcription; Cell wall/membrane/envelope biogenesis; Carbohydrate transport and metabolism; Inorganic ion transport and metabolism; Replication, recombination and repair. Therein, the relative abundances of amino acid transport and metabolism, Energy production and conversion and signal transduction mechanisms in soils under peanut hull mulching were higher than that in soils under polyethylene film mulching.
In addition, we assigned fungal OTUs to specific nutritional groups, and then subdivided them into specific ecological associations (Fig. 6b). FunGuild was a set of designations, which established fungal classification tool with rigorously defined and cited nutrition group assignments. Three ecological functions, including undefined saprotrophs, plant pathogens and soils aprotrophs, were significantly enriched. The relative abundances of undefined saprotroph and soil saprotroph in soils with peanut hull mulching were higher than other under ecological functions.
The mulching patterns on soils with different materials were widely adopted in agricultural cultivation for a long time. The general effect of mulching on soils has also been fully recognized. However, in terms of the influence of different mulching materials on the composition and function of microorganisms in soils, the depth of understanding is far from enough. In the study, in order to reveal the effects of mulching on microbial communities, we analyzed the bacterial and fungal diversity and community composition in tea plantation under polyethylene film or peanut hull using 16S rRNA and ITS rDNA, respectively. The results indicated that mulching could change soil properties and impact soil microbial diversity and community structure, and improve ecosystem functions of tea plantation.
The diversity of soil microbial community was critical to the integrity, stability and sustainability of soil ecosystems . It was affected by different agricultural management, such as different mulching patterns. The mulches of crop residues and straw on soils in the orchard had higher enzymatic activity, microbial biomass, showing that this management to soils under organic mulching was an effective and sustainable agricultural technology. As it produced the increase of the different fractions of organic carbon and microbial activity, it will translate into rapid improvement of soil quality . In an urban land cover experiment, the diversity of microbial community on bark-mulched soils was higher than that of gravel-mulched and unmanaged fields, showing that organic mulches could provide nutrient-rich substrates for fueling microbial richness in soils . However, little research had focused on the patterns of organic and inorganic mulching on soil in tea plantations could affect soil microbial diversity. In the present study, peanut hull mulching increased the diversity of fungal communities in 0–20 cm soil and the richness of bacterial communities in 20–40 cm soil (Table 1). The results indicated that the effect of mulching materials on the diversity of the bacteria was different from that of the fungi. High nutrient contents of peanut hull transported organic matter to soil system, which ultimately caused the quickly reproduction of soil microbes and increased the diversity of microbial communities. On the other hands, we analyzed the distribution of microbial communities in tea plantation soils under different mulching materials. The principal coordinates analysis of microbial communities further revealed the clear separation between organic and inorganic mulching, indicating that there were significant differences in the structures of bacterial and fungal communities in soils among organic and inorganic mulching patterns (Fig. 2). Thus, the organic and inorganic mulching on tea plantation soils constructed different relatively independent micro-ecological system. Compared with film mulching and no mulching, the organic mulching in tea plantation could obtain more abundant and diverse micro-ecological environment.
The composition structure of soil microbial communities were affected by different mulching patterns. In the present study, the results showed that the dominant bacterial phyla in tea plantation were Proteobacteria, Actinobacteria, Acidobacteria and Chloroflexi (Fig. 3a). Proteobacteria are the important phyla in the microbial community. They could be capable to exploit labile carbon sources and had higher relative abundances in nutrient-rich environments [23, 24]. Acidobacteria are the common eosinophil phyla in soils. They had the ability to degrade the cellulose and lignin of plant residues and played an important role in the carbon cycle [25, 26]. High abundances of Proteobacteria and Acidobacteria can enhance cycling of essential nutrients, which can improve soil fertility and sustainable utilization [27, 28]. Our results showed that the relative abundances of Proteobacteria and Acidobacteria in soil under peanut-hull mulching were higher than that with polyethylene-film mulching, indicating that peanut hull as crop residues could provide rich nutrients for Proteobacteria and the inside substance could be used by Acidobacteria (Fig. 4a). Actinobacteria are mostly saprophytes in the soil. They were functionally diverse and contributed to the decomposition of organic matter, which they had ability to break down complex substrates to give them a competitive advantage over other bacteria [29, 30]. Chloroflexi are the important phyla with green pigment. They exhibited oligotrophic lifestyles and grew at a lower rate under conditions of higher nutrient utilization . Our results showed that the relative abundance of Actinobacteria and Chloroflexi in soils with polyethylene film mulching was higher than that in soils with peanut hull mulching, indicating that the application of polyethylene film created a relatively nutrient-poor environment than peanut hull (Fig. 4a). In addition, we further digged into the significant accumulation of Sphingomonas and Nitrospira at on bacterial genus level (Fig. 4e). The genus Sphingomonas had extensive metabolic ability to degrade many aromatic compounds and the certain species of this genus could synthesize valuable extracellular biopolymers . Therefore, the accumulation of Sphingomonas in soils might be caused by mulching, which might affect the growth of tea plant by affecting the root secretions. The genus Nitrospira were involved in nitrification process in many natural environments , and they played important roles in biogeochemical nitrogen cycle in soils . The bacteria in the genus Nitrospira promoted the conversion of nitrite to nitrate by nitrification in soils [35, 36]. The soils covered by peanut hull had the potential advantage of exogenous nitrate intervention, which significantly increased the competitive advantage of Nitrospira. High relative abundance of Nitrospira in soils under peanut-hull mulching indicated that the organic mulching could enrich bacteria to involve in the nitrogen cycle and promote the biochemical cycle in tea plantations.
On the other hands, in the present study, the dominant fungal phyla were Ascomycota, Mortierellomycota and Basidiomycota (Fig. 3b), and the dominant fungal genus were Mortierella, Solicoccozyma, Fusarium and Clonostachys (Fig. 3d). Ascomycota were the main fungal decomposer in many soil ecosystems [37,38,39]. As biotrophs, they played the biggest role in recycling plant residues, and they might form symbioses with mycorrhizae or endophytes . Basidiomycota played an important role in degrading lignin under anaerobic conditions . Our results showed that the relative abundance of Basidiomycota in soils under peanut hull mulching was significantly higher than that in soils under polyethylene-film mulching, indicating that peanut hull could provide anaerobic and high lignin content for the growth of Basidiomycota, and it could make Basidiomycota attend the conversion of lignin under anaerobic conditions (Fig. 4e). At fungal genus level, Mortierella were known as saprobic and ubiquitous. Numbers of studies proved that they had ability to solubilize phosphorus, and were related to increase crop yields and establish symbiosis with plants [42, 43]. In our study, Mortierella showed significant accumulation in soils under peanut hull mulching, indicating that good quality soils in tea plantation were primarily related with high relative abundance of Mortierella. In addition, the genus Fusarium was important crop pathogens in the world. The disease caused by Fusarium could infect different host organs, including plant seeds and roots . They were capable of colonising root surfaces and produced some plant irritant, but they also produced some secondary metabolite toxin during growth and metabolism, which could cause negative plant-soil feedback [45, 46]. Our study showed Fusarium was the dominant fungal genus in tea plantation soil, but the relative abundance of Fusarium in soils under film mulching was higher than that in soils under peanut hull mulching. Therefore, covering with peanut shell might reduce the accumulation of soil borne pathogens and the occurrence of soil disease in tea plantation.
Agricultural practices in tea plantations could alter soil environmental factors . Soil environmental factors, in turn, further affected the composition of microbial communities [47, 48]. Previous studies in an apple orchard and maize system reported that the application of mulch enhanced water-holding capacity and soil nutrients [49, 50]. In the present study of tea plantation, soil moisture and soil TN contents with peanut-hull mulching were significantly higher than those in soils with plastic-film mulching on the topsoil (Table 2), indicating that the organic mulch had better adsorption effect on soil moisture and rainfall, and it could also provide nitrogen source for soil microorganisms. Spearman correlation heatmap analysis showed that the dominant bacterial and fungal phyla were significantly correlated with soil properties (Fig. 5). Soil moisture and soil total nitrogen content were the main physicochemical factors that affect the structure of bacterial community. Meanwhile, the soil moisture and available nitrogen contents were the main physicochemical factors that affect the structure of fungal community in tea plantation. Mulch was beneficial to adjust a microbial community with stronger interaction and could enhance the functional ability of the soil microbes. Therefore, different mulching methods altered the soil microbial communities of tea plantation by regulating soil properties.
Our study highlighted the influence of different mulches on soil microbial diversity and community structure in tea plantation. It evaluated the potential value of organic mulching for tea garden management from two aspects of microbial structure and ecological function. The results showed that organic mulching improved the physicochemical conditions of tea garden soils and induced great changes of soil bacterial and fungal communities, and also improved the soil micro ecological environment. Compared with film mulching and no mulching, the organic mulching in tea plantation could obtain more abundant and diverse micro-ecological environment. So, it is a better choice for tea plantation management. In addition, it is necessary to strengthen the research on the changes of rhizosphere microorganisms under different mulching treatments in the future. At the same time, because the metabolism of soil microorganism is related to the production performance of tea plant, it is also necessary to further research the effect of soil microorganism on the growth and development of new shoots of tea plant.
Field experiment and soil sampling
The field experiment was conducted from May 2019 at Rizhao in Shandong Province of China (35°04′N, 118°25′E). This region has atypical subtropical monsoon climate with an average annual precipitation of 1048 mm and an average temperature of 12.7 °C. The soil of the experimental field is brown soil, which is cleared and hoed manually before the experiment. The major soil parameters were as follows: pH 6.8, organic matter (OM, 11.92 g/kg), available nitrogen (AN, 91.18 mg/kg), available phosphorus (AP, 16.60 mg/kg) and available potassium (AK, 87.56 mg/kg).
The treatments were: (1) CK1: bare soil control, 0-20 cm soil; (2) CK2: bare soil control, 20-40 cm soil; (3) F1: mulching with polyethylene film, 0-20 cm soil; (4) F2: mulching with polyethylene film, 20-40 cm soil; (5) P1: mulching with peanut hull, 0-20 cm soil; and (6) P2: mulching with peanut hull, 20-40 cm soil. Each plot was 90 m2 under unified management and each treatment had three randomly replicate plots. The tea tree variety was ‘Zhongcha108’ with a row spacing of 1.5 m, and plant spacing of 0.33 m. The polyethylene film was black with a thickness of 0.004 mm, and was compacted with soil around the cover. The peanut hull was mulched in the soil surface only with a thickness of 10 cm. Prior to coring, mulch was removed from sampling area to prevent contamination of cores with surface organic matter. The soil sampling in every depth was homogenized by mixing together and passing through a 2-mm sieve for discarding mixed above-ground materials (plant residues and stones). Each soil sample with 3 replications was further divided into two parts: one part was frozen quickly in liquid nitrogen, and stored at − 80°Cuntil the microbial communities were analyzed; the other part was dried naturally in room and then analyzed for pH, available nitrogen (AN), available phosphorus (AP), available potassium (AK), organic matter (OM), total nitrogen (TN), total phosphorus (TP), total potassium (TK).
DNA extraction and PCR amplification
According to the manufacturer’s protocols, microbial DNA from 18 samples was extracted from 18 samples in soils and the E.Z.N.A.® Soil DNA Kit (Omega Bio-tek, Norcross, GA, U.S.) was used. The final DNA concentration and purified concentration and purification were determined by NanoDrop 2000 ultraviolet-visible spectrophotometer, and the DNA quality was checked by 1% agarose gel electrophoresis. The concentration and purification of final DNA were determined by NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Wilmington, USA), and the DNA quality was checked by 1% agarose gel electrophoresis. The DNA samples were individually amplified in V3-V4 hypervariable regions by PCR using primers 338F (5′-ACTCCTACGGGAGGCAGCAG − 3′) and 806R (5′-GGACTACHVGGGTWTCTAAT − 3′) for 16S rDNA in bacteria, and primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA − 3′) and ITS2R (5′-GCTGCGTTCTTCATCGATGC-3′) for ITS in fungi. The procedure used in the PCR reaction had been described in previous study . In short, the PCR reaction was repeated 3 times, and it should go through the procedures of denaturation, annealing, and extension. Then, the resulting PCR product was extracted from a 2% agarose gel and further purified and quantified using the AxyPrep DNA gel extraction kit (Axygen Biosciences, Union City, California, USA) and QuantiFluor™ -ST (according to the manufacturer’s protocol US Promega), respectively.
Illumina sequencing and processing of sequencing data
The purified amplicons were pooled on the Illumina MiSeq platform (Illumina, San Diego, USA) of equal molecular weight and paired-end sequencing (2 × 300) according to the standard protocol of MajorbioBio-Pharm Technology Co. Ltd. (Shanghai, China).
The original sequencing sequence was controlled using Trimmomatic software and merged by FLASH software. The specific criteria are consistent with previous study . UPARSE (version 7.1; http://drive5.com/uparse/) was used to cut the similarity of the operational classification units (OTUs) to 97%. http://drive5.com/uparse/) uses a novel “greedy” algorithm that could perform chimera filtering and OTU clustering at the same time. Using the confidence threshold of 70%, the classification of each 16S rRNA gene sequence was analyzed against the Silva database (Release132; http://www.arb-silva.de) through the RDP classifier algorithm (http://rdp.cme.msu.edu/). Using the confidence threshold of 70%, the classification of each ITS sequence was analyzed against the Unite database (version 7.2; http: //unite.ut.ee/index.php) through the RDP classifier algorithm.
The means and standard deviations of the data were calculated and statistically examined by ANOVA and Duncan’s multiple range tests by SPSS (SPSS, Inc., Chicago, IL, USA). Significant differences were considered as p < 0.05. Alpha diversity was applied to analyze microbial diversity, including the Chao1, Shannon, Ace, Simpson and rarefaction curve, which were calculated with Mothur software (version 1.30.1, http://www.mothur.org/). The principal co-ordinates analysis (PCoA) was applied to reduce the dimension of the original variables based on the Bray-Curtis distances using R2.1.3. The Circos graph was built using Circos-0.67-7 software (http://circos.ca/) . Linear discriminant analysis (LDA) coupled with effect size measurements (LEfSe) analysis was applied to search for statistically different biomarkers between different treatments using LEfSe software (http://huttenhower.sph.harvard.edu/galaxy/). The relationship between soil microbial community structures and soil environment factors was analyzed by Spearman correlation heatmap using pheatmap package.
Availability of data and materials
The raw sequencing data have been deposited in NCBI Sequence Read Archive (SRA) under accession number PRJNA614599 for bacteria, and PRJNA615029 for fungi.
Cluster of Orthologous Groups
Internal transcribed spacer
Linear discriminant analysis Effect Size
Operational Taxonomic Units
Principal co-ordinates analysis
Ribosomal Database Project
Wang YJ, Xie ZK, Malhi SS, Vera CL, Zhang YB, Wang JN. Effects of rainfall harvesting and mulching technologies on water use efficiency and crop yield in the semi-arid Loess Plateau, China. Agric Water Manag. 2009;96(3):374–82.
Chai Q, Gan YT, Turner NC, Zhang RZ, Yang C, Niu YN, Siddique KHM. Chapter two - water-saving innovations in Chinese agriculture. Adv Agron. 2014;126:149–201.
Ni X, Song WT, Zhang HC, Yang XL, Wang LG. Effects of mulching on soil properties and growth of tea olive (Osmanthus fragrans). PLoS One. 2016;11(8):e0158228.
Forge TA, Hogue E, Neilsen G, Neilsen D. Effects of organic mulches on soil microfauna in the root zone of apple: implications for nutrient fluxes and functional diversity of the soil food web. Appl Soil Ecol. 2003;22(1):39–54.
Yang YJ, Dungan R, Ibekwe AM, Valenzuela-Solano C, Crohn D, Crowley D. Effect of organic mulches on soil bacterial communities one year after application. Biol Fertil Soils. 2003;38(5):273–81.
Mulumba LN, Lal R. Mulching effects on selected soil physical properties. Soil Tillage Res. 2008;98(1):106–11.
Kołota E, Adamczewska-Sowińska K. Living mulches in vegetable crops production: Perspectives and limitations (a reviev). Acta Scientiarum Polonorum Hortorum Cultus. 2013;12(6):127–42.
Nacke H, Thürmer A, Wollherr A, Will C, Hodac L, Herold N, Schöning I, Schrumpf M, Daniel R. Pyrosequencing-based assessment of bacterial community structure along different management types in German forest and grassland soils. PLoS One. 2011;6(2):e17000.
Tian W, Wang L, Li Y, Zhuang KM, Li G, Zhang JB, Xiao XJ, Xi YG. Responses of microbial activity, abundance, and community in wheat soil after three years of heavy fertilization with manure-based compost and inorganic nitrogen. Agric Ecosyst Environ. 2015;213:219–27.
Mori AS, Isbell FI, Seidl R. β-diversity, community assembly, and ecosystem functioning. Trends Ecol Evol. 2018;33(7):549–64.
Qian HF, Zhang M, Liu GF, Lu T, Qu Q, Du BB, Pan XL. Effects of soil residual plastic film on soil microbial community structure and fertility. Water Air Soil Pollut. 2018;229(8):261.
Djemiel C, Grec S, Hawkins S. Characterization of bacterial and fungal community dynamics by high-throughput sequencing (HTS) metabarcoding during flax dew-retting. Front Microbiol. 2017;8:2052.
Dong WY, Si PF, Liu E, Yan CR, Zhang Z, Zhang YQ. Influence of film mulching on soil microbial community in a rainfed region of northeastern China. Sci Rep. 2017;7:8468.
Bending GD, Turner MK, Rayns F, Marx M-C, Wood M. Microbial and biochemical soil quality indicators and their potential for differentiating areas under contrasting agricultural management regimes. Soil Biol Biochem. 2004;36(11):1785–92.
Grandy S, Neff JC. Molecular C dynamics downstream: The biochemical decomposition sequence and its impact on soil organic matter structure and function. Sci Total Environ. 2008;404(2–3):297–307.
Lan GY, Li YW, Wu ZX, Xie GS. Soil bacterial diversity impacted by conversion of secondary forest to rubber or eucalyptus plantations: a case study of Hainan Island, South China. For Sci. 2017;63(1):87–93.
Romaniuk R, Giuffre L, Costantini AO, Nannipieri P. Assessment of soil microbial diversity measurements as indicators of soil functioning in organic and conventional horticulture systems. Ecol Indic. 2011;11(5):1345–53.
Wu JS, Lin HP, Meng CF, Jiang PK, Fu WJ. Effects of intercropping grasses on soil organic carbon and microbial community functional diversity under Chinese hickory (Carya cathayensis Sarg.) stands. Soil Res. 2014;52(6):575–83.
Cui H, Zhou Y, Gu ZH, Zhu HH, Fu SL, Yao Q. The combined effects of cover crops and symbiotic microbes on phosphatase gene and organic phosphorus hydrolysis in subtropical orchard soils. Soil Biol Biochem. 2015;82:119–26.
Zhao J, Ni T, Li Y, Xiong W, Ran W, Shen B, Shen QR, Zhang RF. Responses of bacterial communities in arable soils in a rice-wheat cropping system to different fertilizer regimes and sampling times. PLoS One. 2014;9(1):e85301.
García-Orenes F, Guerrero C, Roldán A, Mataix-Solera J, Cerdà A, Campoy M, Zornoza R, Zornoza GB, Caravaca F. Soil microbial biomass and activity under different agricultural management systems in a semiarid Mediterranean agroecosystem. Soil Tillage Res. 2010;109(2):110–5.
Martir-Torres MC, Bruns MA. Comparative diversity and abundance of ammonia monooxygenase genes in mulched and vegetated soils. Soil Biol Biochem. 2013;57:758–68.
Fierer N, Bradford MA, Jackson RB. Toward an ecological classification of soilbacteria. Ecology. 2007;88:1354–64.
Kuramae E, Yergeau E, Wong LC, Pijl A, Veen JAV, Kowalchuk G. Soil characteristics more strongly influence soil bacterial communities than land-use type. FEMS Microbiol Ecol. 2012;79:12–24.
Ward N, Challacombe J, Janssen PH, Henrissat B, Kuske CR, et al. Three genomes from the phylum acidobacteria provide insight into the lifestyles of these microorganisms in soils. Appl Environ Microbiol. 2009;75(7):2046–56.
Huang XQ, Liu LL, Wen T, Zhu R, Zhang JB, Cai ZC. Illumina MiSeq investigations on the changes of microbial community in the Fusarium oxysporum f.sp. cubense infected soil during and after reductive soil disinfestation. Microbiol Res. 2015;181:33–42.
Smit E, Leeflang P, Gommans S, Broek JVD, Mil SV, Wernars K. Diversity and seasonal fluctuations of the dominant members of the bacterial soil community in a wheat field as determined by cultivation and molecular methods. Appl Environ Microbiol. 2001;67:2284–91.
Torsvik V, Øvreås L. Microbial diversity and function in soil: from genes to ecosystems. Curr Opin Microbiol. 2002;5(3):240–5.
Kamau DM, Spiertz JHJ, Oenema O. Carbon and nutrient stocks of tea plantations differing in age, genotype and plant population density. Plant Soil. 2008;307:29–39.
Nemergut DR, Cleveland C, Wieder WR, Washenberger CL, Townsend A. Plot-scale manipulations of organic matter inputs to soils correlate with shifts in microbial community composition in a lowland tropical rain forest. Soil Biol Biochem. 2010;42(12):2153–60.
Tardy V, Chabbi A, Charrier X, Berranger CD, Reignier T, Dequiedt S, Faivre-Primot C, Terrat S, Ranjard L, Maron PA. Land use history shifts in situ fungal and bacterial successions following wheat straw input into the soil. PLoS One. 2015;10(6):e0130672.
Fredrickson JK, Balkwill D, Drake GR, Romine M, Ringelberg D, White DC. Aromatic-degrading sphingomonas isolates from the deep subsurface. Appl Environ Microbiol. 1995;61(5):1917–22.
Bartosch S, Hartwig C, Spieck E, Bock E. Immunological detection of nitrospira-like bacteria in various soils. Microb Ecol. 2002;43(1):26–33.
Hayatsu M, Tago K, Saito M. Various players in the nitrogen cycle: diversity and functions of the microorganisms involved in nitrification and denitrification. Soil Sci Plant Nutr. 2008;54:33–45.
Ahn JH, Yu R, Chandran K. Distinctive microbial ecology and biokinetics of autotrophic ammonia and nitrite oxidation in a partial nitrification bioreactor. Biotechnol Bioeng. 2008;100(6):1078–87.
Lücker S, Wagner M, Maixner F, Pelletier E, Koch H, Vacherie B, Rattei T, Damsté JSS, Spieck E, Paslier DL, Daims H. ANitrospirametagenome illuminates the physiology and evolution of globally important nitrite-oxidizing bacteria. Proc Natl Acad Sci U S A. 2010;107(30):13479–84.
Lei W, Yang XL, Rachel K, Wang Y, Tong DL, Ye M, Xin J. Combined use of alkaline slag and rapeseed cake to ameliorate soil acidity in an acid tea garden soil. Pedosphere. 2013;23(2):177–84.
Wang ZJ, Chen Q, Liu L, Wen XX, Liao YC. Responses of soil fungi to 5-year conservation tillage treatments in the drylands of northern China. Appl Soil Ecol. 2016;101:132–40.
Zhou J, Jiang X, Zhou BK, Zhao BS, Ma MC, Guan DW, Li J, Chen SF, Cao FM, Shen DL, Qin J. Thirty four years of nitrogen fertilization decreases fungal diversity and alters fungal community composition in black soil in northeast China. Soil Biol Biochem. 2016;95:135–43.
Ivancevic B, Karadelev M. Overview of fungi species in Prespa National Park (ALBANIA). Int J Ecosyst Ecol Sci. 2013;3(4):679.
Boer WD, Folman L, Summerbell R, Boddy L. Living in a fungal world: impact of fungi on soil bacterial niche development⋆. FEMS Microbiol Rev. 2005;29:795–811.
Fröhlich-Nowoisky J, Hill TCJ, Pummer BG, Franc GD, Pöschl U. Ice nucleation activity in the widespread soil fungus mortierella alpina. Biogeosciences. 2015;12(4):1057–71.
Grządziel J, Gałązka A. Fungal biodiversity of the most common types of polish soil in a long-term microplot experiment. Front Microbiol. 2019;10(6):1–12.
Chen T, Christensen M, Nan ZB, Hou FJ. Effects of grazing intensity on seed size, germination and fungal colonization of Lespedeza davurica in a semi-arid grassland of northwest China. J Arid Environ. 2017;144:91–7.
Mangan SA, Schnitzer S, Herre EA, Mack KML, Valencia MC, Sanchez EI, Bever JD. Negative plant-soil feedback predicts tree-species relative abundance in a tropical forest. Nature. 2010;466:752–5.
Yang M, Zhang XD, Xu YG, Mei XY, Jiang BB, Liao JJ, Yin ZB, Zheng JF, Zhao Z, Fan LM, He XH, Zhu YY, Zhu SS. Autotoxic ginsenosides in the rhizosphere contribute to the replant failure of panax notoginseng. PLoS One. 2015;10(2):e0118555.
Shen ZZ, Wang DS, Ruan YZ, Xue C, Zhang J, Li R, Shen QR. Deep 16S rRNA pyrosequencing reveals a bacterial community associated with banana fusarium wilt disease suppression induced by bio-organic fertilizer application. PLoS One. 2014;9(5):e98420.
Xun WB, Huang T, Zhao J, Ran W, Wang BR, Shen QR, Zhang RF. Environmental conditions rather than microbial inoculum composition determine the bacterial composition, microbial biomass and enzymatic activity of reconstructed soil microbial communities. Soil Biol Biochem. 2015;90:10–8.
Liu Y, Gao MS, Wu W, Tanveer SK, Wen XX, Liao YC. The effects of conservation tillage practices on the soil water-holding capacity of a non-irrigated apple orchard in the Loess Plateau, China. Soil Tillage Res. 2013;130:7–12.
Dong QG, Yang YC, Yu K, Feng H. Effects of straw mulching and plastic film mulching on improving soil organic carbon and nitrogen fractions, crop yield and water use efficiency in the Loess Plateau, China. Agric Water Manag. 2018;201:133–43.
Li P, Xue Y, Shi JL, Pan AH, Tang XM, Ming F. The response of dominant and rare taxa for fungal diversity within different root environments to the cultivation of Bt and conventional cotton varieties. Microbiome. 2018;6:184–98.
Yan MM, Chen SN, Huang TL, Li BQ, Li N, Liu KW, Zong RR, Miao YT, Huang X. Community compositions of phytoplankton and eukaryotes during the mixing periods of a drinking water reservoir: dynamics and interactions. Int J Environ Res Public Health. 2020;17:1128.
We would like to thank all authors for valuable discussions.
This research was subsidized by the Significant Application Projects of Agriculture Technology Innovation in Shandong Province (SD2019ZZ010), the Technology System of Modern Agricultural Industry in Shandong Province (SDAIT-19-01) and the Special Foundation for Distinguished Taishan Scholar of Shangdong Province (No.ts201712057), the Livelihood Project of Qingdao City (19–6–1-64-nsh), the Project of Agricultural Science and Technology Fund in Shandong Province (2019LY002, 2019YQ010), the Key Research and Development Project of Shandong province (2019LYXZ009).
Ethics approval and consent to participate
Not applicable. The tea plants used in this study were provided by Tea Research Institute of Qingdao Agricultural University. It does not require ethical approval.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, S., Wang, Y., Sun, L. et al. Organic mulching positively regulates the soil microbial communities and ecosystem functions in tea plantation. BMC Microbiol 20, 103 (2020). https://doi.org/10.1186/s12866-020-01794-8
- Camellia sinensis (L.) O. Kuntze
- Bacterial and fungal communities
- Tea plantation soil