不同植物类型复合垂直流人工湿地根系微生物群落结构的研究

蒋旭瑶，吉喜燕，黄德英，张继彪

（1.复旦大学环境科学与工程系，上海 200433；2.复旦大学化学系，上海 200433）

近年来，人工湿地由于其造价低、便于维修管理、处理效率高、观赏性强等优点，被广泛应用于污水处理中[1-2]。然而，湿地在净化水质的同时会释放大量温室气体，越来越受到广泛关注[3]。微生物和植物是人工湿地的重要组成部分，其污水处理能力取决于湿地设计、微生物群落和植物种类等[4]。微生物在湿地生物地球化学循环中起着重要作用；不同的微生物群落会产生不同的生物化学反应，导致营养盐的消耗[5-6]。因此，湿地植物根系微生物群落结构的研究有助于改进人工湿地的设计以及管理，提高湿地水质净化能力，并对湿地温室气体排放的控制提供理论依据[7]。湿地中微生物的硝化、反硝化作用是去除氮的主要途径，硝化、反硝化细菌是反应的主要执行者[8]。因此，对湿地中硝化细菌、反硝化细菌展开分析有助于湿地污水脱氮机制的研究。

过去，人们针对湿地微生物种群的研究一直建立在传统的分离和培养方法上，这种方法费时费力，培养的种群数量也有限[9]。目前，关于湿地土壤微生物的研究大多运用传统分子生物学技术，如PCR-DG⁃GE和FISH[10]。PCR-DGGE技术是将PCR扩增和变性梯度凝胶电泳结合起来，可对同一长度的PCR扩增片断按序列的不同在变性梯度凝胶上分离，在分析微生物群落多样性方面具有明显的优势，但是它只适用于微生物比较丰富的样本[10]。FISH技术是根据已知微生物不同分类级别上种群特异的DNA序列，以荧光标记的特异寡聚核苷酸片段作为探针，与环境基因组中DNA分子杂交，该技术可以进行样品原位杂交，但受引物以及试样不合适的限制[10]。高通量DNA测序是一种测序深度较深、覆盖度较广的方法，该技术通过大量平行序列产生大量DNA数据，其中的可操作性分类单元（OTUs）可用于有效分析样本中微生物的多样性[10]。

本研究选择美人蕉、菖蒲和水葱三种典型的湿地植物，通过构建复合垂直流人工湿地研究不同季节、上行池和下行池土壤根际微生物群落结构的变化特征，分析根系微生物结构和功能的对应关系，旨在为人工湿地的改进和优化提供参考。

1 材料与方法

1.1 试验系统

复合垂直流人工湿地系统位于复旦大学邯郸校区研究基地，在自然光照下的露天大棚内运行，示意图如图1所示。每个系统由上、下行池组成，尺寸为150 cm×90 cm×90 cm（长×宽×高），下行池内基质自下而上分别为20 cm砾石（ϕ30～50 mm）、20 cm砾石（ϕ 10～20 mm）、20 cm细砂和10 cm土壤；上行池自下而上分别为20 cm砾石（ϕ30～50 mm）、20 cm砾石（ϕ10～20 mm）、10 cm细砂和10 cm土壤。每个系统种植一种植物，分别为菖蒲、美人蕉和水葱，每平方米分别种植20棵、20棵和20簇。模拟生活污水从湿地下行池上端流入，从上行池上端流出。模拟生活污水主要是往自来水中加入葡萄糖、淀粉、鱼粉蛋白胨、硝酸钾、磷酸氢二钾、硫酸镁等。进水水质总有机碳、铵态氮、总氮和总磷的浓度分别为68～84、20～28、50～59 mg·L-1和4～6 mg·L-1。进水采用间歇流的方式，一周配一次1000 L污水，污水放置在一个1000 L的塑料桶中，水力负荷为330 mm·d-1，停留时间为5 d。系统从2015年6月开始运行，到2017年5月份结束，运行稳定，植物长势良好。

图1 复合垂直流人工湿地系统示意图Figure 1 Picture of integrated vertical flow constructed wetland

1.2 样品采集

本试验于2016年7月、10月和2017年1月用五点采样法分别采集三种植物湿地系统的上行池、下行池10 cm处根系土壤。样品标记为：ACL表示菖蒲，CIL表示美人蕉，SVV表示水葱；7、10、1表示月份；D表示下行池，U表示上行池。样品名称分别为：ACL7D、ACL7U、CIL7D、CIL7U、SVV7D、SVV7U、ACL10D、 ACL10U、 CIL10D、 CIL10U、 SVV10D、SVV10U、ACL1D、ACL1U、CIL1D、CIL1U、SVV1D、SVV1U。

1.3 DNA提取

采用Omega E.Z.N.A.®Soil DNA Kit提取试剂盒提取根系土壤DNA，按说明书完成相关操作，然后将提取好的DNA保存在-20℃。

1.4 PCR扩增与测序

本实验采用引物520F（5′-barcode+GCACCTA⁃AYTGGGYDTAAAGNG-3′）和 802R（5′-TACNVGGGTATCTAATCC-3′）[11]，扩增长度约为 280 bp的 16S rRNA V4区用来测序。

PCR采用NEB Q5 DNA高保真聚合酶，各组成及用量如下：0.25 μLQ5 high-fidelity DNA polymerase，5 μL 5*Reaction Buffer，5 μL 5*High GC Buffer，0.5 μL dNTP（10 mmol·L-1），1 μL模板 DNA，1 μL 正向引物（10 μmol·L-1），1 μL反向引物（10 μmol·L-1）和11.25 μL水。PCR扩增程序如下：将PCR反应所需的成分配制完后，在PCR仪上于98℃预变性30 s，使模板DNA充分变性，然后进入扩增循环。在每一个循环中，先于98℃保持15 s使模板变性，然后将温度降到50℃，保持30 s，使引物与模板充分退火；在72℃保持30 s，使引物在模板上延伸，合成DNA，完成一个循环。重复这样的循环25～27次，使扩增的DNA片段大量累积。最后，在72℃保持5 min，使产物延伸完整，4℃保存。

利用Illumina公司的TruSeq Nano DNA LT Li⁃brary Prep Kit进行建库；在MiSeq机器上进行MiSeq Reagent Kit V3（600 cycles）2×300 bp的双端测序。

1.5 数据分析

对样品进行Alpha多样性分析时，在90%的最低测序深度水平，采取对OTU丰度矩阵中全体样本进行随机重复抽样，并使用Chao1指数来表示群落丰富度。样品进行Beta多样性分析时，对Unweighted Uni⁃Frac距离矩阵和Weighted UniFrac距离矩阵进行NMDS分析，用于描述群落样本的结构分布。

2 结果与讨论

图2 样品的稀疏曲线Figure 2 The rarefaction curves for samples

样品的稀疏曲线如图2所示，所有样品测序量在5260时，曲线已趋于平缓或者达到平台期，因此可以认为测序深度已经基本覆盖到样品中所有的物种。

三种植物湿地微生物多样性指数如图3所示。水葱湿地在实验期间，上行池与下行池的物种丰富度差异不大。除10月份美人蕉湿地下行池Chao1指数小于上行池Chao1指数外，菖蒲湿地与美人蕉湿地表现为其下行池的丰富度大于上行池。这是因为下行池为污水输入端，营养盐含量较高，适宜微生物代谢繁殖。流入上行池的污水经过下行池的降解处理，能被微生物利用的营养盐降低，导致微生物丰富度下降。

图3 三种植物湿地微生物多样性指数Figure 3 Diversity indices of microbes in three plant wetlands

本研究用NMDS分析（Nonmetric Multidimension⁃al Scaling），通过降维处理简化数据结构，在特定距离尺度下描述样本分布特征。图4a为Unweighted Uni⁃Frac NMDS分析的样本二维排序图。除ACL7U样品，样本组均匀分布在坐标上下两侧，表明在试验期间上行池和下行池样本组群落成员相似，但是两者群落成员组成之间有差异，存在各自独有的微生物。这可能是因为下行池消耗了大量氧气，导致上行池氧气较低，使得下行池中有氧细菌远多于上行池中有氧细菌数[12]。

由图4b可知，除ACL7U样品，上行池样品组靠的较近，表明从群落成员丰度上来看，上行池样本相似性较下行池样本相似性高，这主要是因为下行池是污水净化的主要场所，下行池的营养盐含量和溶解氧含量等相关影响因素沿程变化梯度较大，造成微生物生存环境变化较大，三种植物湿地在不同环境影响下形成了差异略大的生态系统。而经过下行池的水流，其中的污染物浓度已经较低，其浓度、溶解氧等变化较为平缓，导致微生物所处环境较为相似[13]，因此上行池样品组相似性高于下行池。

图4 Unweighted UniFrac NMDS分析的样本二维排序图（a）和Weighted UniFrac NMDS分析的样本二维排序图（b）Figure 4 Unweighted UniFrac NMDS analysis of samples in two dimensions（a）and Weighted UniFrac NMDS analysis of samples in two dimensions（b）

门和属水平上微生物群落组成及丰度如图5所示，在菖蒲湿地中（图5左侧a），主要的微生物有变形菌门（Proteobacteria，38.4%）、酸杆菌门（Acidobacteria，16.9%）、绿弯菌门（Chloroflexi，13.1%）。变形菌门所占比例较大，这与已有研究一致[14-15]。酸杆菌门相对丰度在菖蒲湿地下行池夏秋冬季节、上行池夏秋冬季节占有较大比例，分别为15.6%、15.8%、17.0%、17.2%、16.5%、19.5%。绿弯菌门的相对丰度在菖蒲下行池夏秋冬季节、上行池夏秋冬季节分别占14.5%、11.7%、11.2%、11.0%、16.5%、13.8%。

美人蕉湿地的微生物群落结果如图5左侧b所示，下行池根系的变形菌门相对丰度在秋季最大，达到50.9%；上行池根系变形菌门相对丰度在夏季、秋季和冬季差异不大，维持在33.0%到34.2%之间。总体来说，下行池变形菌门相对丰度高于上行池。酸杆菌门在下行池、上行池的相对丰度夏季最高，分别为20.4%、18.7%，且下行池酸杆菌门相对丰度高于上行池。绿弯菌门下行池夏季较高，达到14.4%，上行池秋季较高，达到16.6%。

水葱湿地根系微生物主要类别同为变形菌门（34.6%）、酸杆菌门（17.5%）、绿弯菌门（10.8%）（图5左侧c）。随季节演替，变形菌门在下行池的相对丰度从46.3%降到31.1%。上行池相对丰度随季节变化不明显，基本维持在30.0%。下行池酸杆菌门的相对丰度夏季＞秋季＞冬季，分别为17.0%、16.1%、10.7%，上行池秋季＞夏季＞冬季，分别为21.7%、20.0%、19.6%。绿弯菌门下行池秋季＞夏季＞冬季，分别为13.0%、9.4%、6.8%。

三种湿地植物根系微生物组成在门水平上相似，主要有变形菌门、酸杆菌门、绿弯菌门等。变形菌门比例较高，包含许多微生物都与碳、氮、硫循环相关[16]，在去除碳、氮方面有重要的作用，同时对温室气体的生成有巨大影响。酸杆菌门是严格厌氧细菌，用来发酵芳香族混合物和醋酸盐[17]。绿弯菌门在营养盐的去除过程中发挥着重要作用[18]。同时湿地中还检测到了浮霉菌门（Planctomycetes）、芽单胞菌门（Gemmatimonadetes）、放线菌门（Actinobacteria）、拟杆菌门（Bacteroidetes）、泉古菌门（Crenarchaeota）、厚壁菌门（Firmicutes）等。浮霉菌门能将NO-2作为电子受体，将氨转化为氮气，对NH3-N的去除有重要影响[19]。芽单胞菌门中的一些菌种已被确认为革兰氏阴性、好氧、积累磷酸盐微生物[20]。放线菌门将氢气作为电子供体用来去除硝酸盐，对生成N2O有重要影响。疣微菌门相对丰度为0.8%，能氧化CH4，对其排放有抑制作用[21]。

图5 门和属水平上微生物群落在样本中组成及丰度Figure 5 Taxonomic composition at phylum and genus level and abundance of microbes in each sample

属水平上的微生物群落组成及相对丰度如图5右侧图所示，三种湿地根系中主要的微生物属有Ni⁃trospira、Candidatus Nitrosophaera等。其中，Nitrospira是硝化细菌；Candidatus Nitrosophaera是亚硝化细菌。不同植物湿地优势菌属有类似，也有差异。硝化细菌最主要的微生物属为Nitrospira，下行池在秋季达到最高，相对丰度为2.44%。

菖蒲湿地在属水平上的硝化、反硝化细菌组成及其相对丰度如表1所示。Nitrospira和Candidatus Ni⁃trosophaera是主要的硝化细菌，下行池、上行池相对丰度均在秋季达到最高，分别为6.64%、4.33%。反硝化微生物有Dechloromonas、Flavobacterium、Pseudomo⁃nas属等，这与已有研究结果一致[22-25]。菖蒲下行池的根系反硝化微生物在夏季所占比例最多，约1.34%，秋季与冬季差异不明显，上行池反硝化细菌相对丰度随季节变化不明显，基本维持在0.62%～0.78%。下行池反硝化微生物相对丰度高于上行池。下行池硝化细菌、反硝化细菌相对丰度都比上行池高，这可能是因为下行池是进水端，含有丰富的氮源、碳源供微生物生长活动[13]。

在美人蕉湿地中，下行池植物根系的硝化细菌Nitrospira和Candidatus Nitrosophaera相对丰度在夏季达到最高3.60%，上行池秋季达到最高3.62%（表2）。反硝化细菌主要有Dechloromonas、Flavobacterium、Pseudomonas等。下行池和上行池根系反硝化细菌相对丰度分别在秋季和夏季达到最高，分别为2.30%和0.98%。

表1 菖蒲湿地硝化细菌、反硝化细菌的相对丰度表（%）Table 1 The relative abundance of nitrifying bacteria and denitrifying bacteria in ACL cell（%）

表2 美人蕉湿地硝化细菌、反硝化细菌的相对丰度表（%）Table 2 The relative abundance of nitrifying bacteria and denitrifying bacteria in CIL cell（%）

在水葱湿地中，硝化细菌为Nitrospira和Candida⁃tus Nitrosophaera（表3），上、下行池的相对丰度在秋季达到最高，分别占4.16%、4.78%。反硝化细菌有De⁃chloromonas、Flavobacterium、Pseudomonas等，下行池反硝化细菌相对丰度在秋季达到最高，为1.51%；上行池冬季达到最高，为1.55%。

总之，湿地植物根系硝化细菌属为Nitrospira、Candidatus Nitrosophaera，反硝化细菌属为Dechlo⁃romonas、Flavobacterium、Pseudomonas等。菖蒲和水葱湿地秋季硝化细菌丰度均高于夏季和冬季。

3 结论

（1）菖蒲湿地、美人蕉湿地和水葱湿地系统Chao1指数均表现为秋季＞夏季＞冬季，即微生物丰富度秋季＞夏季＞冬季。水葱湿地上行池与下行池物种丰富度差异不大；除秋季的美人蕉湿地，菖蒲湿地与美人蕉湿地Chao1指数表现为下行池＞上行池。

（2）在试验期间，各湿地上行池和下行池样本组群落组成均相似；但上、下行池群落成员组成之间有差异，上行池样本相似性较下行池样本高。

（3）门水平上，湿地根系微生物从大到小依次为变形菌门、酸杆菌门、绿弯菌门等。从属水平上，检测到了硝化细菌Nitrospira和Candidatus Nitrosophaera；反硝化细菌Dechloromonas、Flavobacterium和Pseudo⁃monas等。菖蒲湿地和水葱湿地根系硝化细菌秋季的相对丰度均高于夏季和冬季。

表3 水葱湿地硝化细菌、反硝化细菌的相对丰度表（%）Table 3 The relative abundance of nitrifying bacteria and denitrifying bacteria in SVV cell（%）