20世纪70年代至21世纪初，竹类植物对缓解森林资源供需矛盾和协调森林经济效益做出了巨大贡献^[1]。第八次全国森林资源清查（2009-2013年）数据表明，我国竹林面积为6.01×10⁶ hm²，其中毛竹林占70%以上^[2]。毛竹（Phyllostachys edulis）是我国最主要的笋材两用竹种，其用途范围广泛，造林后可持续经营，已在南方诸省得到规模化推广。然而，现有毛竹林多数为中小径级低产林，且面临着林地衰退、经济下滑、劳动力和农资成本不断上升等问题^[2]，因此，对中小径级毛竹林进行大面积改造的意义不大。综合考虑其潜在的利用价值，改善其生态系统环境，探索新的“春笋冬出”地表覆盖模式已成为毛竹林可持续经营亟待解决的问题。

利用农林剩余物覆盖的技术兴起于澳大利亚和美国，目前在我国的应用亦十分广泛。我国农林剩余物多为有机物料，种类丰富，如砻糠和果壳、秸秆类（包括稻草、麦秆和玉米秆等）、树叶（包括竹叶）和锯末（包括木材和竹材）等^[3-4]。覆盖技术具有平衡土壤水分和肥力，改善土壤-大气的界面层气流和热量状况等多方面的生态效应^[5]，已成为许多国家和地区广泛采用的土壤管理调控技术之一^[6-7]。研究发现，地膜或杂草覆盖有蓄水保墒的作用，可促进土壤功能微生物群的活动^[3]。一方面，施肥和地膜覆盖使土壤微生物多样性受短期土壤养分环境改变的影响^[8]；另一方面，长期管理措施使土壤结构和组成成分发生明显的改变，从而影响微生物的群落结构、功能及代谢活性^[7]。针对中小径级毛竹林经济效益日趋衰退的现状，实施以冬季覆盖地表增温和重度施肥为技术核心的“春笋冬出”集约栽培可实现巨大的经济效益。然而，毛竹林地覆盖经营中大量引入有机覆盖物和化学肥料（氮肥为主），极易导致生态系统退化，爆发式萌发休眠芽，打破自身固有的生长发育与繁殖生活史节律，是一种被动适应^[2]，严重影响毛竹林的可持续经营。研究表明，土壤微生物是维持土壤系统稳定性、健康和质量的关键因素^[9]，对植物生长也具有重要作用^{[7, 10-11]}。

长期覆盖和重度施肥模式对环境和资源进行了高强度人为干扰，被大量引入毛竹林地的外源有机物料与毛竹林地自然凋落物的分解过程不同，而微生物活动是控制分解过程最为关键的因子^[7]。在这种强选择压下，人工毛竹林系统与近自然毛竹林生态系统存在怎样的差异、其土壤细菌多样性如何变化、轮休是否为覆盖作业的必要补充、以及如何提升土壤微生物多样性与群落结构等问题亟待解决。因此，本试验以中小径级的低产毛竹林为研究对象，基于覆盖和轮休的不同状态，分析覆盖毛竹林的土壤微生物群落变化特征，探讨不同经营状态对毛竹林生态系统土壤微生物结构的影响，为中小径级低产毛竹定向培育和衰退毛竹林生态修复技术的研究提供技术支撑和数据参考。

毛竹覆盖试验地位于浙江省湖州市长兴县小浦镇葡萄岕（119°45′51″ E，30°58′59″ N），属亚热带海洋性季风气候。年平均气温（15.6±0.5）℃，有效积温5 750 ℃，无霜期239 d，年平均降水量1 309 mm，年平均日照时数1 810.3 h，历年平均日照百分率为41%，光照分配较均匀。土壤以发育于酸性岩浆岩的红壤为主。试验地土壤层深厚，可达1.5 m。林地海拔77~81 m，坡度＜5°，坡向（东南坡，SE）161°~166°，土壤有机质平均质量分数为30.2 g/kg，全氮1.16~2.24 g/kg，有效磷98.00~ 101.76 mg/kg，速效钾54.75~92.61 mg/kg。

试验地从2009年开始进行覆盖轮休（3年覆盖+2年轮休为1个周期）经营，面积共2 hm²。覆盖试验地母竹平均胸径9.4 cm。当年留养1~5年生母竹2 700~3 300株/hm²。5-6月进行垦复，垦复深度10~15 cm，并抛施行鞭肥[750 kg/hm²复合肥，硫酸钾型m（N）:m（P₂O）:m（K₂O）=15:15:15；以下所用复合肥都为此成分]；8-9月施孕笋肥（1 500 kg/hm²复合肥）；10-11月，满园灌水，竹林土壤湿度＞80%，然后施7 500 kg/hm²鸡粪有机肥（鸡粪干物质约含粗蛋白19.3%、脂肪2.4%、灰分13.1%、糖类10.3%、纤维7.2%，氮2.45%、磷2.66%、钾0.91%）+1 875 kg/hm²竹炭肥（C/N=14.7，pH 8.13~ 8.50，电导率5.47 mS/cm，含水率30.5%，N 11.9%，P₂O₅ 4.8%，K₂O 5.7%）^[2]，覆土；11-12月，对毛竹林地进行双层覆盖，下层稻草（15~20 cm）起发酵增温作用，同时施水，上层谷壳（20~25 cm）起保温作用。轮休区域1~5年生母竹2 250~3 000株/hm²，其中2、3、4年生竹比例达到80%，1年和5年生竹控制在20%；同时，5-6月初垦复并施行鞭肥（750 kg/hm²复合肥），8-9月施孕笋肥（750 kg/hm²复合肥），并平整地块。

样地划分为对照（从未实施覆盖和施肥）、第2年覆盖、第2年轮休和2009-2015年长期覆盖，样地状况详见表 1。2016年3月进行土壤取样。每块样地内设置5 m×5 m的样方5个，沿样方对角线每隔1~2 m确定取样点，去除表层凋落物和杂质后，用土钻取深度为0~20 cm土壤，样方内的土样混合均匀，去除根和石头，放置于干冰低温采样箱内。4~6 h内带回实验室，保存于-80 ℃超低温冰箱中，用于DNA的提取和后续试验。

表1(Table 1)

表1 样地经营状态 Table 1 Status of sample plots

表1 样地经营状态 Table 1 Status of sample plots 样地 Sample plots 状况 Status 面积 Area/hm2 pH 土壤含水量 Soil water content/% 土壤全氮 Soil total N/(g/kg) 土壤有机碳 Soil organic carbon/ (g/kg) 对照 CK 未施肥覆盖，自然状态 0.66 5.0 ± 0.2a 23.26± 0.71b 1.07 ± 0.07c 10.08 ± 2.21c 覆盖 Mulching 轮休2年，覆盖第2年 0.34 4.4 ± 0.1b 25.47±3.17ab 1.89 ± 0.07a 20.01 ± 0.66b 轮休 Alternate operation 覆盖3年，轮休第2年 0.34 4.5 ± 0.2b 20.40±1.27c 1.27 ± 0.10b 13.90±1.38c 长期覆盖 Long-term mulching 2009-2015年一直覆盖 0.66 4.1 ± 0.2c 25.73± 1.72a 2.09 ± 0.15a 25.38± 1.81a 同列数据后的不同小写字母表示在P＜0.05水平差异有统计学意义。 Values within a column followed by different lowercase letters represent statistically significant difference at the 0.05 probability level.

点击放大

DNA提取采用TIANamp Stool DNA试剂盒[购自天根生化科技（北京）有限公司]。对微生物基因组16S rDNA v3/4高变区进行扩增，引物序列如下。上游引物：CCTACGGGNGGCWGCAG；下游引物：GACTACHVGGGTATCTAATCC。

聚合酶链式反应（polymerase chain reaction, PCR）体系：模板量约为10 ng；1×反应缓冲液（NEB Q5TM）、0.3 mmol/L dNTP、1 U HotStarTaq聚合酶共10 µL；10 µmol/L上游引物和10 µmol/L下游引物各1 µL（Invitrogen公司，美国）；然后加Ambion无核酸酶水（Thermo Fisher Scientific公司，美国）至20 µL。PCR反应条件：94 ℃预变性2 min；94 ℃变性20 s，52 ℃退火40 s，72 ℃延伸1 min，35个循环；72 ℃保持2 min；4 ℃保存。样本添加特异性标签序列，采用以下PCR反应体系：模板量1 µL（Amplicon PCR纯化后产物）；1×反应缓冲液（NEB Q5TM）、0.3 mmol/L dNTP、1 U Q5TM DNA聚合酶（NEB公司，美国）共10 µL；10 µmol/L上游引物和10 µmol/L下游引物各1 µL（Invitrogen公司，美国）；Ambion无核酸酶水（Thermo Fisher Scientific公司，美国）7 µL。PCR反应条件：98 ℃预变性30 s；94 ℃变性10 s，65 ℃退火30 s，72 ℃延伸30 s，11个循环；72 ℃保持5 min。每个样品的PCR产物割胶回收后用紫外-可见分光光度计（NanoDrop ND1000，美国）测量其浓度。

文库构建和Illumina MiSeq测序由杭州晶佰生物科技有限公司进行。

测序后的DNA序列采用PANDAseq进行拼接^[11]。拼接后的序列应用USEARCH除去嵌合体^[12]，然后用QIIME软件将所有样品序列进行整合^[13]。1个微生物可操作单元（operational taxonomic unit, OTU）为序列相似度大于97%的DNA序列。基于Greengenes v13_08数据库提取OTU并进行种属分类。去掉序列占比＜0.005%的OTU，然后进行后续分析^[14]。用QIIME软件计算多样性指数，包括Shannon指数、Chao1指数等。同时，比较不同样品原核生物群落的差异。利用SPSS 11.0软件的最小显著差数法（least-significant difference, LSD）进行多重比较分析；利用R语言软件的非加权UniFrac（unweighted UniFrac）进行土壤细菌群落结构的主坐标分析（principal co-ordinate analysis, PCoA）。

对照、覆盖、轮休和长期覆盖4类样地的样本覆盖度（最终测序结果占整个基因组的比例）分别为（98.99±0.12）%、（99.14±0.26）%、（98.92±0.20）%和（99.12±0.24）%，均大于95%，表明测序结果能代表样品中的细菌组成。经高通量测序发现，中小径级毛竹林4类经营状态下的样地细菌类群达到26门69纲118目197属。其中：对照样地的细菌类群为24门61纲99目132科156属，覆盖样地的细菌类群为24门62纲97目129科151属，轮休样地的细菌类群为24门60纲97目132科160属，长期覆盖样地的细菌类群为24门58纲88目123科150属。

由图 1可知：对照样地细菌类群在门的分类水平上的优势类群（平均所占比例大于5%）分别为酸杆菌门（Acidobacteria，44.22%）、AD3（14.21%）、变形菌门（Proteobacteria，8.23%）、绿弯菌门（Chloroflexi，7.91%）和疣微菌门（Verrucomicrobia，6.82%）；覆盖样地细菌类群门水平上的优势类群分别为酸杆菌门（42.38%）、变形菌门（17.79%）、TM7（6.11%）、绿弯菌门（5.77%）和放线菌门（Actinobacteria，5.07%）；轮休样地细菌类群门水平上的优势类群分别为酸杆菌门（30.85%）、AD3（15.18%）、变形菌门（14.72%）、放线菌门（8.39%）和绿弯菌门（7%）；长期覆盖样地细菌类群门水平上的优势类群分别为酸杆菌门（34.43%）、变形菌门（25.41%）、TM7（6.68%）和放线菌门（6.65%）。比较各菌门数量发现：上述菌门中仅酸杆菌门在4类样地中都大量存在；对照样地中变形菌门、放线菌门、TM7和拟杆菌门（Bacteroidetes）数量较其他类型样地低，AD3、疣微菌门和GAL15的数量较其他类型样地高；覆盖样地的AD3和疣微菌门数量较少；长期覆盖样地中绿弯菌门和AD3数量较少。

图1(Figure 1)

图1 4种经营类型毛竹林土壤细菌门分类水平分布 Fig. 1 Phylum distribution of soil bacteria in four types of sample plots

从表 2中可以发现，在对照、覆盖、轮休和长期覆盖4块样地的20个土壤样品中，在属的分类水平上分别检测到（97.66 ± 0.67）%、（97.63 ± 0.33）%、（95.04±0.39）%和（98.58±0.51）%的明确分类信息。在属分类水平上（≥1.00%），Koribacteraceae科未定属（Koribacteraceae; g_）和Ellin6513目未定属（Ellin6513;f_; g_）在所有样地中都大量存在，且各样地间在统计学上差异显著（P＜0.05）；与对照样地相比较，覆盖、轮休与长期覆盖样地细菌比例明显上升（P＜0.05）的有酸杆菌科未定属（Acidobacteriaceae; g_）、Solibacteraceae科未定属（Solibacteraceae; g_）、酸微菌目未定属（Acidimicrobiales; f_;g_）、Gaiellaceae科未定属（Gaiellaceae; g_）、Chitinophagaceae科未定属（Chitinophagaceae; g_）、鞘脂杆菌科未定属（Sphingobacteriaceae; g_）、纤线杆菌纲未定属（Ktedonobacteria; o_;f_; g_）、WD2101目未定属（WD2101;f_; g_）、变形菌门所有8个属、TM7-1纲未定属（TM7-1;o_; f_;g_）、auto67-4W科未定属（auto67-4W;g_）和WPS-2门未定属（WPS-2;c_; o_; f_; g_）；而与其他3类样地比较（P＜0.05），对照样地细菌比例明显上升的有Koribacteraceae科未定属（Koribacteraceae; g_）、Ellin6513目未定属（Ellin6513;f_; g_）、AD3门2个未定属（ABS-6;o_; f_; g_和JG37-AG-4;o_; f_;g_）、B07-WMSP1目未定属（B07-WMSP1;f_; g_）、GAL15门未定属（GAL15;c_; o_; f_;g_）、芽单胞菌门Gemm-1目未定属（Gemm-1;o_; f_;g_）、Ellin515科未定属（Ellin515;g_）和DA101属（DA101）。

表2(Table 2)

表2 4类样地细菌属分类水平的百分比(≥1.00%) Table 2 Percentage of bacterial community composition at genus level in four types of sample plots (≥1.00%)

表2 4类样地细菌属分类水平的百分比(≥1.00%) Table 2 Percentage of bacterial community composition at genus level in four types of sample plots (≥1.00%) 门 Phylum 属 Genus 对照 CK/% 覆盖 Mulching/% 轮休 Alternate operation/% 长期覆盖 Long-term mulching/% 酸杆菌门 Acidobacteria RB41;f_; g_ 2.18 ±0.09a 0.01±0.00b 0.01±0.00b 0.00±0.00b Acidobacteriaceae; g_ 0.67± 0.10c 11.13 ± 0.78a 6.12±1.13b 11.00 ± 1.61a Koribacteraceae; g_ 13.64±0.67a 10.91±1.79b 8.21±0.95b 6.69±0.42c Candidatus_Koribacter 1.84±0.22b 2.94±0.40a 1.56±0.07c 1.99±0.5 lb Ellin6513;f_; g_ 15.41±0.70a 8.98±0.59b 8.51±1.06b 6.03 ± 1.18c DS-18;f_; g_ 2.08±0.09a 0.05±0.00c 0.19± 0.02b 0.01 ± 0.00d Solibacterales; f_;g_ 2.73±0.40a 3.02±0.28a 1.95±0.27b 3.10±0.33a Solibacteraceae; g_ 0.09±0.01c 0.89±0.45ab 0.33±0.17b 1.59± 0.20a Candidatus_Solibacter 3.48±0.72a 3.62±0.84a 3.21±0.60a 2.96±0.37a 放线菌门 Actinobacteria Acidimicrobiales; f_;g_ 0.58±0.04c 0.97±0.1 lb 0.93±0.27b 1.37± 0.25a Actinomycetales; f_;g_ 2.02±0.57a 2.88±0.74a 3.81 ± 1.02a 2.22±0.60a Gaiellaceae; g_ 0.37± 0.05c 0.49 ± 0.08b 0.65±0.1 lb 1.91±0.05a AD3 ABS-6;o_; f_;g_ 10.25 ± 1.27a 1.92 ± 0.54b 9.55±1.02a 0.86±0.32c JG37-AG-4;o_; f_;g_ 3.96±0.77b 1.08±0.06c 5.62±0_53a 0.18±0.01d 拟杆菌门 Bacteroidetes Chitinophagaceae; g_ 0.42±0.12c 1.46±0.17b 1.06±0.3 lb 1.70±0.13a Sphingobacteriaceae; g_ 0.03±0.00c 1.71 ± 0.40a 0.28±0.03b 1.39±0.08a 绿弯菌门 Chloroflexi Ktedonobacteria; o_;f_; g_ 0.31±0.04c 1.75±0.21a 0.70±0.06b 0.68±0.11b Thermogemmatisporaceae; g_ 3.29 ± 0.60b 3.35 ± 0.47b 5.50±0.5 la 0.61±0.29c B07-WMSPl; f_;g_ 2.56±1.03a 0.15±0.01c 0.62±0.14b 0.18±0.00c GAL 15 c_; o_;f_; g_ 2.64±0.83a 0.15 ±0.00c 0.48±0.05b 0.09±0.00c 芽单胞菌门 Gemmatimonadetes Gemm-1;o_; f_;g_ 2.28±0.86a 0.44±0.05b 0.18±0.01c 0.60±0.12b 浮霉菌门 Planctomycetes WD2101;f_; g_ 0.70±0.09c 2.34±0.45b 1.62±0.27b 3.44±0.52a 变形菌门 Proteobacteria Ellin329;f_; g_ 0.82±0.25c 2.90±0.88b 1.71±0.35b 4.52±0.96a Rhodoplanes 0.92±0.04c 1.73±0.16b 1.58±0.21b 2.29 ± 0.24a Acetobacteraceae; g_ 0.09±0.01b 1.13±0.21a 0.71±0.24a 1.27± 0.23a Rhodospirillaceae; g_ 1.58±0.22c 3.57±0_53b 3.57±0.67b 4.81±0.39a Kaistobacter 0.27±0.05c 0.71±0.16b 0.93±0.12ab 1.34±0.32a Gammaproteobacteria; o_;f_; g_ 0.06±0.01c 1.10± 0.22a 0.69±0.08b 1.52±0.47a Sinobacteraceae; g_ 1.05±0.08c 1.72±0.14b 1.69±0.22b 3.62 ± 0.28a Xanthomonadaceae; g_ 0.30±0.05b 0.74±0.18a 0.98±0.21a 1.15±0.25a TM7 TM7-1;o_; f_;g_ 0.46 ± 0.03c 5.46±0.82ab 3.55 ± 1.37b 6.49 ± 0.48a 疣微菌门 Verrucomicrobia auto67-4W;g_ 0.07±0.02c 1.08±0.28b 0.71±0.15b 2.92±0.44a Ellin515;g_ 1.21±0.53a 0.26±0.07b 0.22±0.02b 0.17±0.00c DA101 4.70±0.61a 0.08±0.00c 1.03±0.26b 0.09±0.00c WPS-2 c_; o_;f_; g_ 2.60±0.24b 4.26±0.64a 4.18±1.07a 4.16± 0.70a 4类样地如有1类样地细菌属出现比例≥1.00%即进行多重比较（LSD）；同行数据后的不同小写字母表示在P＜0.05水平差异有统计学意义。 Multiple comparisons were performed if the quantitative proportion of one genus is equal or greater than 1.00% among the four types of sample plots. Values within a row followed by different lowercase letters represent statistically significant difference at the 0.05 probability level.

点击放大

α-多样性种类复杂度指数是对单个样本物种多样性的分析^[15]。本试验对原始序列条数、OTU数、Shannon指数和Chaol指数进行了评价，结果（图 2）显示：对照样地的原始序列条数最低，为（22 052± 370）条，与其他3类样在统计学上差异显著（P＜0.05），说明对照样地细菌群落在总体数量上处于弱势；而覆盖、轮休和长期覆盖的原始序列条数差异不显著，其中长期覆盖样地的原始序列条数平均值较低（24 437条）。OTU数表明，对照与长期覆盖样地差异不显著，而两者与覆盖和轮休样地都有显著差异（P＜0.05）。Shannon指数表明，对照样地的细菌物种多样性最低（8.33±0.33），且与其他3类样地有显著差异（P＜0.05），覆盖与轮休两者间差异不显著，说明覆盖使毛竹林土壤的细菌物种多样性上升。长期覆盖样地的Chaol指数最低（1 342±34），且与其他3类样地有显著差异（P＜0.05）；覆盖轮休样地的Chaol指数最高（1 662±36），而与覆盖样地差异不显著，与对照和长期覆盖样地均有显著差异（P＜0.05），说明长期覆盖样地使土壤细菌物种总数即丰度下降，而覆盖和轮休模式使土壤细菌丰度上升。

图2(Figure 2)

短栅上不同小写字母表示在P＜0.05水平差异有统计学意义。 Different lowercase letters above bars indicate significant differences at the 0.05 probability level. 图2 4类样地物种种类复杂度 Fig. 2 Species complexity in four types of sample plots

PCoA是一种研究数据相似性或差异性的可视化方法，由物种组成的距离矩阵计算得到^[16]。从图 3中可以看出：PC₁、PC₂和PC₃维度分别解释了70.69%、8.47%和7.71%的信息；供试的20个样地被分成2大部分，其中长期覆盖、覆盖和轮休样地从PC₂由低向高分布，距离较近，而且各类型样地在图中的样点明显聚集，说明三者的土壤细菌群落结构组成有相似之处，且每一类型样地内部的细菌群落结构十分接近；而对照样地与前三者的距离较远，说明其与前三者的土壤细菌结构及物种组成有较大差异。

图3(Figure 3)

图3 4种经营类型土壤细菌群落结构的主坐标分析 Fig. 3 Principal co-ordinate analysis of soil bacterial community diversity in four types of sample plots

森林土壤微生物数量和群落结构受林地管理制度、土壤类型、土壤理化性质、植被类型等多种因素的影响^{[9, 17]}，特别是不同林地管理方式对土壤的营养平衡和物质组成等产生作用，直接影响了土壤中微生物的多样性、群落组成和数量^[18-20]，进而对土壤中物质、能量的转化和营养物质的有效性产生影响^[18-19]。本试验结果表明，在中小径级毛竹林的对照、覆盖、轮休和长期覆盖4类不同经营措施样地中，土壤细菌种群多样性覆盖了26门69纲118目197属，其中酸杆菌在4类样地中占优势（30.85%~ 44.22%），其次是变形菌（8.23%~25.41%）、绿弯菌（2.48%~7.91%）、疣微菌（2.39%~6.82%）、AD3（1.05%~15.18%）、放线菌（3.99%~8.39%）和TM7（0.61%~6.68%）。这与KAMLESH等^[17]对不同类型耕作土壤细菌种群结构研究时发现的变形细菌居于主要地位而酸杆菌其次不同。本研究中酸杆菌占绝对优势而变形菌其次，而HANSEL等^[20]在对不同地带性土壤的研究中也发现变形菌优势比酸杆菌大，造成这种差异的原因可能与毛竹林土壤原始的细菌类群组成结构有关，亦可能与森林土壤经营方式和土地利用方式不同有关^[16]。本研究还发现，除变形菌和酸杆菌外，绿弯菌、疣微菌、AD3、放线菌和TM7在毛竹林土壤中占有一定优势。同时，在属的分类水平上，同一门的不同属对毛竹林覆盖技术的响应不同，如酸杆菌，这可能是由于酸杆菌门总丰度与pH值呈反比，但在不同pH值土壤中存在的酸杆菌门群组不同^[21]。

土壤微生物数量巨大、种类繁多，传统分离培养的方法遗漏了土壤中绝大多数微生物，难以反映土壤微生物群落的丰度及其作用机制^[12-13]。本试验对中小径级覆盖毛竹林土壤总DNA的16S rDNA v3/4高变区基因进行PCR扩增，利用Illumina Miseq高通量测序技术和生物信息学分析不同经营类型土壤细菌群落组成差异。研究发现，与其他3类样地相比，长期覆盖样地酸杆菌门3个未定属（Koribacteraceae; g_，Ellin6513;f_; g_和DS-18;f_; g_）、AD3门2个未定属（ABS-6;o_; f_;g_和JG37-AG-4;o_; f_;g_）、绿弯菌门的Thermogemmatisporaceae; g_和疣微菌门的Ellin515; g_数量显著下降，表明土壤微生物的分布与土壤环境密切相关，用土壤微生物相关指标来反映土壤肥力具有较大的应用潜力。

覆盖技术（包括重度施肥）对土壤的主要物理化学性质和微生物群落结构会产生显著影响，如造成土壤养分失衡、竹龄结构不合理、酶活性异常、重金属含量增加等一系列负面效应^[22]，严重影响竹林的可持续经营。本试验结果显示：与对照样地相比，其他3类样地土壤全氮和有机碳含量均提高。覆盖、轮休和长期覆盖3类样地的pH值显著下降，而且随着覆盖时间延长，土壤酸化越明显，这与多数竹林覆盖研究相同^{[5, 22]}。对照样地近自然经营状态，没有施肥，仅有采笋作业，故其土壤全氮和有机碳含量最低，覆盖和长期覆盖样地在大量使用有机覆盖物、有机肥和复合肥后，其土壤全氮和土壤有机碳含量明显增高，但长期覆盖的细菌数量、种群结构多样性却下降，这可能与土壤微生物对有机碳的利用效率较高有关^{[19, 23]}；轮休样地在去除毛竹林内剩余有机质和浅层翻耕后，土壤全氮和有机碳含量较覆盖和长期覆盖都有所下降，对重度施肥后维持毛竹林土壤碳氮平衡颇有裨益。不同土壤施肥管理制度与外界环境条件的综合作用，造成了土壤理化性质的显著差异，进而影响土壤微生物群落结构和多样性^{[4-5, 7, 17]}。本试验覆盖（主要为有机覆盖物、有机肥和复合肥）后的土壤pH值下降，这为酸杆菌提供了适宜的生长环境，使其多样性和相对丰度增加；同时，土壤有机质和养分含量提高，变形菌大部分类群为异养或者兼性营养菌，有机底物和养分为这部分微生物类群的生长提供了可溶性营养物质和能量^[19]。此外，由于土壤有机质及矿物质颗粒共同参与土壤团聚体形成，因此，在一定经营时间内，覆盖技术改良了毛竹林土壤，为土壤微生物提供了良好的栖居环境，调控了土壤微生物群落的代谢功能和结构^{[3, 24]}。

土壤养分是土壤微生物的能量源，同时依赖土壤微生物的转化而得到补充，因此微生物和土壤质量之间存在复杂的相互作用^{[16, 18]}。张瑞等^[25]认为，添加有机肥有利于土壤微生物代谢活性和微生物群落功能的多样性，单施氮磷钾复合肥和氮肥促进了少数几类微生物生长，而其他类群微生物相对丰度降低。WESSEN等^[26]发现，在裸地使用化肥可提高土壤微生物群落的多样性和丰度，但长期施化肥则使细菌和古菌多样性下降。袁红朝等^[19]认为，季节性林地覆盖对土壤固氮菌数量的影响不显著，对土壤细菌、真菌、放线菌及氨化菌数量的影响极显著。本试验中的α-多样性分析和PCoA分析表明，覆盖经营对土壤细菌的种群多样性和群落结构产生了较大影响。覆盖提高了细菌数量、种类和丰度，产生这种差异的原因在于大量C/N比高的有机质为细菌提供了基本的碳源和能源，有利于其数量和丰度的增加；而土壤表面浅层翻耕处理使其不受机械干扰，短时间内形成0~20 cm厚的半腐熟有机物质层，同样为细菌繁殖提供能源，这与高云超等^[27]在半干旱半湿润地区石灰性土壤试验中发现的浅层翻耕土壤有利于细菌生长的结果类似。同时，有机物料和有机肥（腐熟鸡粪）自身携带的微生物也增加了覆盖后竹林土壤细菌的多样性，且土壤里残余的大量纤维素在分解过程中产生的中间产物（有机酸、糖类等）能促进细菌的增殖，说明在一定时间范围内，重度配比施化肥、氮磷钾复合肥和有机肥可以促进土壤微生物生长。生态系统生物多样性对于系统稳定性有所影响，就群落水平而言，种群多样性高的群落，系统稳定性强，抵抗干扰的能力强^{[16, 20, 24]}。本试验发现，轮休样地的细菌数量、种类和丰度均为最高，说明轮休恢复土壤细菌群落多样性的效果明显，长期覆盖在细菌数量、种类和丰度上均有所下降，其Chao1指数甚至显著低于对照样地，说明长期覆盖样地对土壤细菌群落已产生抑制作用，干扰强度过大已超出微生物的恢复能力^[23]，不利于长期可持续经营。