1. 引言
近年来,随着我国经济的高速发展及工业化进程的加快,土壤重金属污染的程度也不断加剧。据2014年的《全国土壤污染状况调查公报》报道,受Cd、Pb、Hg、Cr、Zn、Ni、As、Cu等重金属污染的土壤占19.4%,其中,受Cd、As、Pb等重金属污染的耕地面积达2 × 105 km2,约占耕地总面积的1/5 [1]。目前,土壤重金属污染治理主要有物理修复、化学修复和生物修复。物理修复成本较高,化学修复易对环境造成二次污染;生物修复操作简单,但占地面积大,修复时间较长,且目前发现的超累积植物(hyperaccumulator)数量不多,生物量小,因而该技术仍未得到大面积的推广应用。目前,对土壤重金属污染修复的研究大多集中在草本植物方面,而对木本植物少有涉及。
河池地处桂西北,矿产资源丰富,其中锌、锡、锑、铟产量在全国占有重要地位,所辖的南丹县是全国最大的锡和锌生产基地,还伴生有大量的镉、银、钨、铅、铜、锑、铟、镓、镉、硒等稀散金属及硫、砷等非金属矿产 [2]。本研究从乡土树种对土壤重金属的耐受性、吸收能力及根际土壤微生物多样性角度考虑,利用lltmina高通量测序技术,从细菌多样性、群落结构及其与土壤化学性质的相互关系等方面,比较分析矿区与非矿区乡土树种根际土壤细菌对重金属的响应程度,以期在今后的实际运用中,分别发挥植物、微生物修复的优点,构建利用乡土树种根际土壤微生物协同处理矿区或其他重金属异常区土壤重金属污染的技术方法。
2. 材料与方法
2.1. 样品采集与处理
采样点主要选择在河池市境内的选矿厂、冶炼厂、矿物运输的交通要道等区域。以南丹矿区为中心,尽量选择重金属污染比较严重的区域进行野外调研和样品采集,用GPS进行定位。
2019年10月中旬至2020年5月下旬,先后多次到河池有色金属矿区、冶炼厂及土壤重金属异常区周围进行调研采样。选择长势比较一致的植株,每种乡土树种均采集3株以上不同植株的根际土壤样品和叶片样品。为了更好地说明不同重金属胁迫对植物根际土壤微生物多样性的影响,同时采集非矿区相应乡土树种的根际土壤样品和叶片样品作为对照。
土壤样品采集时,按梅花形采样法采集相应植株树冠垂直投影范围内0~20 cm的根际土,用四分法取样,一部分土壤样品存于灭菌的塑料自封袋中,低温保存,尽快带回实验室,并放入−86℃超低温冰冻箱中保存,用于土壤微生物DNA序列测定及微生物多样性分析;另一部分土壤样品经自然风干、研磨过筛后,密封于塑料封口袋中,供土壤理化性质及重金属含量的测定。
采集叶片样品时,尽量采集不同部位的正常叶片,保持叶片完整,同时做好相应的采样记录。将叶片样品带回实验室后,及时用清水和去离子水洗涤,超声洗净后,置于105℃的电热鼓风恒温干燥箱中杀青30 min,再将温度调至60℃烘干,研磨粉碎,过100尼龙筛,装入塑料封口袋中,阴凉干燥处保存,备用。
2.2. 实验方法
2.2.1. 土壤理化性质的测定
水分含量测定用重量法,土壤全氮测定用半微量凯氏法 [3],土壤全磷测定用NaOH熔融—钼锑抗比色法 [4],土壤全钾的测定用NaOH熔融—火焰光度法 [4]。
2.2.2. 重金属含量的测定
土壤样品和叶片样品重金属含量的测定,用湿法消解,电感耦合等离子体质谱仪(ICP-MS)进行 [5] [6] [7] [8] [9]。
2.2.3. 土壤微生物DNA序列的测定
利用Fast kit DNA试剂盒提取根际土壤总DNA,储存于−20℃冰箱中保存备用。扩增16S rRNA基因的V3-V4区域,并送至MiSeq测序平台测序 [10]。其主要流程为:样品制备→核酸质检→PCR扩增→产物纯化制备→文库质检→MiSeq上机测序。
2.2.4. 数据处理
每个样品均做3个重复,结果取平均值,数据处理分析用Excel 2010和SPSS 22.0进行 [11]。高通量测序数据使用QIIME1 Pipeline V1.9.0处理和分析原始序列 [12]。ANCOM,ANOVA,Kruskal Wallis,LEfSe和DEseq2等方法被运用来鉴定分组和样本间丰度有差异的细菌 [13]。
Alpha多样性以及Beta多样性的分析主要用qiime2 diversity插件完成。Alpha多样性指数是对某个样品中物种多样性的分析,包含样品中的物种组成的丰富度和均匀度两个因素,通常用ObservedOTU,Shannon以及Faith’s Phylogenetic Diversity等指数来评估某个样本的物种多样性,指数越高,表明样本的多样性越复杂。
β多样性分析(β diversity)根据样本的OTU丰度信息计算Bray Curtis,Weighted Unifrac和Unweighted Unifrac距离来评估不同样品间的微生物群落结构差异。
3. 结果与分析
3.1. 不同乡土树种根际土壤的理化分析
不同乡土树种根际土壤基本理化特性的测定结果见表1。由此可知,不同乡土树种根际土壤pH为4.93~7.81,风干土含水量为1.61%~6.58%,土壤有机质含量为21.47 g∙kg−1~98.9 g∙kg−1,总氮含量为0.15 g∙kg−1~1.77 g∙kg−1,总磷含量为0.33 g∙kg−1~2.71 g∙kg−1,总钾含量为2.42 g∙kg−1~4.00 g∙kg−1;不同采样点之间,乡土树种根际土壤的理化特性有显著差异。
![](Images/Table_Tmp.jpg)
Table 1. Comparison of basic physical and chemical characteristics of rhizosphere soil of different native tree species
表1. 不同乡土树种根际土壤基本理化特性的比较
注:① KQ和FKQ分别为在矿区和非矿区采集的样品;② LFH、GBL、YL、SJ分别为枫香、银杏、紫玉兰和国槐,下同;③ 表中MC、SOM、TN、TP、TK分别指土壤含水量、有机质、全氮、全磷和全钾含量;④ 同一列中,不同大写字母表示相同数据的理化特性具有极显著差异,p < 0.01。
3.2. 不同乡土树种根际土壤及叶片样品的重金属含量
不同乡土树种根际土壤及叶片样品重金属含量的测定结果见表2。由此可知,生长在矿区的乡土树种根际土壤中的Mn、Zn、As、Pb含量比较高,其中,LFH-2中Mn含量达5160.00 mg/kg,LFH-7中Zn含量达8020 mg/kg,LFH-6中As含量达5910 mg/kg,Pb含量达2760 mg/kg;GBL-2中Mn含量达2920 mg/kg、Zn含量达4260 mg/kg,GBL-4中As含量达7120 mg/kg,GBL-3中Pb含量达2630 mg/kg;YL-5中Mn含量达4504 mg/kg;Zn含量达13821.23 mg/kg、As含量达9513.90 mg/kg、Pb含量达5066.17 mg/kg;SJ-6中Mn含量达达3284 mg/kg;Zn含量达16453.35 mg/kg、As含量达9184.52 mg/kg、Pb含量达4196.127 mg/kg。
与根际土相比,叶片样品的重金属含量相对较低;根际土重金属含量较高的乡土树种,其叶片中的重金属含量也相对较高;与矿区植物叶片样品相比,非矿区植物叶片样品的重金属含量相对较低。在不同的乡土植物中,LFH及GBL叶片中Mn含量和Zn含量高;YL及SJ叶片中Mn、Zn、As、Pb含量也比其他金属元素含量高。
![](Images/Table_Tmp.jpg)
Table 2. Comparison of heavy metal content in rhizosphere soil and leaf samples of different native tree species
表2. 不同乡土树种根际土壤及叶片样品重金属含量的比较
3.3. 土壤细菌测序数据分析
3.3.1. 乡土树种根际土壤细菌Alpha多样性指数统计分析
本项目研究的样本共24个,将其分为两组对细菌群落结构进行分析,其中对照组(FKQ)共7个样本,实验组(KQ)共17个样本。高通量焦磷酸测序共产生619,835条有效序列,其中实验组共17个样本,获得445,538条有效序列;对照组7个样本,共获得174,297条有效序列。利用独立样本t检验(双侧)对所有样本有效序列的物种丰富度指数(species richness)和物种多样性指数(species diversity)进行比较分析,不同组别样本对于shannon,chao1均不会表现出显著性(p > 0.05),但明显可看出细菌多样性指数FKQ > KQ。结果见图1和图2。
![](//html.hanspub.org/file/20-2181573x11_hanspub.png)
Figure 2. Box plot of Shannon index
图2. Shannon指数的箱型图
3.3.2. 乡土树种根际土壤细菌群落结构分析
β多样性分析(β diversity)是用来比较物种多样性方面存在的差异大小,本文采用Unweighted Unifrac指数衡量β多样性,结果见图3 (p < 0.05)。由此可推断KQ与FKQ组间菌群结构存在显著差异。利用该指数得到的分析图中两点之间的距离越远,表示两个样本之间的微生物群落结构相似度越低,差异越大。除个别样本存在偏差外,矿区乡土树种实验组(KQ)和非矿区乡土树种对照组(FKQ)的菌群差异基本呈现在相同区域。
![](//html.hanspub.org/file/20-2181573x13_hanspub.png)
Figure 3. Analysis of soil microbial β diversity
图3. 土壤微生物β多样性分析
3.3.3. 土壤细菌群落在门水平上的相对丰度
采用16SrRNA测序技术对矿区与非矿区乡土树种根际土壤样品细菌门水平多样性分析,结果见图4,由此可知,FKQ组中的优势菌门为变形菌门(26.95%)、放线菌门(29.21%)、酸杆菌门(18.17%)、绿弯菌门(9.93%)等,KQ组中的优势菌门为变形菌门(30.37%)、放线菌门(20.23%)、酸杆菌门(17.32%)、绿弯菌门(14.38%)、硝化螺旋菌门(4.51%)、芽单胞菌门(3.86%),其中,KQ组变形菌门、绿弯菌门、硝化螺旋菌门、芽单胞菌门的丰度明显高于FKQ组。
![](//html.hanspub.org/file/20-2181573x14_hanspub.png)
Figure 4. Column chart of the relative distribution of soil bacterial communities at the phylum level
图4. 土壤细菌群落在门水平上的相对分布情况柱形图
3.3.4. 土壤细菌群落在属水平的相对丰度
采用16SrRNA测序技术对矿区与非矿区乡土树种根际土壤样品细菌属水平多样性分析,结果见图5。从属的相对丰度来看,KQ组种的优势菌为红游动菌属(4.66%)、芽孢杆菌属(1.5%)、节细菌属(1.85%)、分枝杆菌属(1.01%);FKQ组种的优势菌为红游动菌属(3.99%)及一些未指明的菌属等。从属水平来看,KQ乡土树种根际特征微生物为根瘤菌属、硝化螺菌属,其相对丰度在对应分组中较高;FKQ乡土树种根际特征微生物为褐螺菌属、芽孢杆菌属、小月菌属,其相对丰度在对应分组中较高。具体情况见图6。
![](//html.hanspub.org/file/20-2181573x15_hanspub.png)
Figure 5. Column chart of the relative distribution of soil bacterial communities at the genus level
图5. 土壤细菌群落在属水平上的相对分布情况柱形图
![](//html.hanspub.org/file/20-2181573x16_hanspub.png)
Figure 6. LDA bar graph of LEfSe analysis
图6. LEfSe分析LDA柱形图
3.3.5. 微生物种群与环境因素相关性分析
基于RDA分析的结果表明,在所有被测化学性质中,pH和TP与细菌群落分布和种群分布间的相关性最大(P值分别为0.0005和0.0035);从重金属含量分析,As、Cu、Pb与细菌群落分布和种群分布间也有明显的相关性,其P值分别为0.0425、0.0035、0.0440。具体情况见图7。
注:在CCA/RDA物种排序图内,环境因子用箭头表示,箭头连线的长度代表某个环境因子与群落分布和种类分布间相关程度的大小(解释方差的大小),箭头越长,说明相关性越大,反之越小。
Figure 7. Correlation between different environmental factors and the genera of soil bacteria
图7. 不同环境因子与土壤细菌属水平的相关性
其中,生丝微菌属与Zn、Cu、Pb具有极显著的正相关,与Mn和As具有显著正相关性;硝化螺菌属和芽胞梭菌属与Zn、Pb具有极显著正相关性,具体情况见图8。
4. 讨论
4.1. 重金属对根际土壤细菌多样性的影响
随着矿物开采、金属冶炼加工、化工生产、化肥农药施用等人类活动强度的不断增加,重金属污染已成为当今人类社会面临的突出问题。土壤是一个复杂的体系,其容纳量有限,随着有毒有害物质的积累,最终必定会导致土壤环境的恶化。土壤微生物是土壤的重要组成部分,其活性和群落结构变化是衡量土壤质量的重要指标 [14]。相关研究结果表明,土壤微生物群落的变化,可以敏感地反映土壤重金属污染的状况 [15]。本研究结果表明,从细菌Alpha-多样性Chao指数和Shanno指数来看,FKQ均大于KQ。一些研究者发现,土壤受到重金属污染后,在一定时间内,微生物多样性呈下降趋势,但这种影响因土壤类型、重金属污染剂量和污染周期长短而异 [16] [17] [18],本研究所得的结果与之相似。
注:X轴上为环境因子,Y轴为物种。通过计算获得R值(秩相关)和P值。R值在图中以不同颜色展示,P值若小于0.05则用*号标出,右侧图例是不同R值的颜色区间,同时,左边的色条标注了物种所属门分类。*0.01 ≤ P < 0.05,**0.001 ≤ P < 0.01,***P < 0.001。
Figure 8. Heat map of the correlation between microbial species and phenotype at the genus level
图8. 属水平微生物物种与表型之间的相互关系热图
4.2. 重金属对根际土壤细菌群落结构的影响
实验结果表明,土壤重金属背景不同,即KQ与FKQ 对乡土树种根际细菌群落结构(Beta–多样性)有显著影响。基于PCoA结果发现,KQ与FKQ细菌群落沿Axis.1轴出现显著分离,说明KQ与FKQ根际土壤细菌群落结构的差异明显。
通过对实验结果的比较分析可知,在门水平上,在KQ组与FKQ组中,根际细菌群落变形菌门、放线菌门、绿弯菌门、硝化螺旋菌门、芽单胞菌门为优势菌门,这与赵立君 [19] 等人的相关研究结果一致;在属水平上,KQ组乡土树种根际特征微生物为根瘤菌属、硝化螺菌属,说明这两种特征菌属对重金属有较强的耐受性,由此推测,生活在矿区的乡土树种,之所以能够很好地适应矿区重金属污染的环境,并且生长良好,可能与其根际特征微生物密切相关。硝化螺菌属于硝化细菌,硝化细菌是完成氮循环的关键细菌 [20],它能将土壤中的氨或铵盐转变为硝酸盐,增加植物可利用的氮素;根瘤菌属是一类固氮细菌,可通过改变土壤中的氮环境来提高植物的氮循环,增加植物的抗逆性,从而对重金属矿区和工业废弃地的土壤污染起到修复作用,这与黄兴如 [21] 及张洁 [22] 等人的研究结果一致。
4.3. 根际细菌群落与土壤环境的相互关系
利用RDA分析细菌群落时发现,土壤理化性质中的pH和TP对细菌群落分布的有较大的影响;As、Cu、Zn、Pb等重金属元素对乡土树种根际土壤微生物具有较大影响,其中,硝化螺菌属和芽胞梭菌属与Zn、Pb具有极显著的正相关,这与王梦姣 [23] 等人的研究结果相似。
5. 结论
通过以上实验数据及分析讨论,可以初步得到以下结论:
1) 根际土壤微生物对植物吸收重金属具有重要的调控作用。
2) 乡土树种枫香、银杏、紫玉兰、国槐对重金属元素Mn、Zn、As、Pb有较强耐性,之所以能在重金属含量较高的土壤环境中正常生长,与其根际土壤细菌群落结构多样性密切相关。
3) 重金属对不同乡土树种根际土壤细菌群落具有较大的影响,其中,KQ组变形菌门、绿弯菌门、硝化螺旋菌门、芽单胞菌门的丰度明显高于FKQ组;在属水平上,矿区乡土树种(KQ组)根际土壤特征微生物为硝化螺菌属与根瘤菌属,可能通过改变土壤中的氮环境,增强植物的固氮作用和氮循环,提高植物代谢能力,从而强化植物对重金属的耐性。
基金项目
广西高校重点实验室项目(桂教科研〔2010〕6号),河池学院硕士专业学位建设基金课题(2016HJA007),河池学院高层次人才科研启动费项目(XJ2018GKQ016),河池学院重点实验室项目(校政发〔2016〕91号),广西高校大学生创新创业计划训练项目(202010605014)。
NOTES
*通讯作者。