ISSN1006-7167
CN31-1707/T
RESEARCHANDEXPLORATIONINLABORATORY
第40卷第1期 Vol.40No.1
2021年1月Jan.2021 
DOI:
10.19927/j.cnki.syyt.2021.01.006
高通量测序技术在水生态污染监测中的应用
景 明, 陈
 瑜, 高 昕, 李继影
(江苏省苏州环境监测中心,江苏苏州215000)
摘 要:针对张家港某河道水异常事件,通过高通量测序技术结合显微镜鉴定结果,判断引起河道生态污染的微生物种类,并分析引起该事件的原因。通过Illumina测序,对河道的原核生物及真核生物分别进行测序,真核生物量较多且落组成较为单一,原核生物量相对较少且落分布较为均匀,与显微镜检测结果一
致。对获得的OTU进行物种注释,得到引起河道生态污染的主要微生物为
Pedospumellaencystans。本研究为水生态监测提供了新的视角,也证明高通量测序技术在水生态监测中具有较好的应用潜力。关键词:高通量测序;真核生物;原核生物
中图分类号:X835   文献标志码:A   文章编号:
1006-7167(2021)01-0027-05ApplicationofHigh throughputSequencingTechnologyin
MonitoringofWaterEcologicalPollution
JINGMing, CHENYu, GAOXin,
 LIJiying
JiangsuSuzhouEnvironmentalMonitoringCenter,Suzhou215000,Jiangsu,China)Abstract:
ThisstudyfocusesontheabnormalwatercolorinariverinZhangjiagang.Thereasonforabnormalwateroftheriverisoutbreakofamicroorganismbyhigh throughputsequencingandmicroscope.Thepossiblecausesofthiseventwereanalyzedinthispaper.TheprokaryotesandeukaryotesintheriverweresequencedrespectivelybyIlluminasequencing.Theresultshowsthatthecommunitycompositionofprokaryotesisrelativelyrich
,andtheproportion
differenceofdominantmicroorganismsissmall,however,
thecompositionofeukaryoticcommunityisverysimpleandconsistentwiththemicroscopicresults.
OTUrepresentativesequencesweretaxonomicallyclassifiedandthe
microorganismthatcausethisecologicalpollutionwasPedospumellaencystans.Thisstudyprovidesanewperspectiveforwaterecologymonitoringandprovesthathigh throughputsequencingtechnologyhasgreatapplicationpotentialinwaterecologymonitoring.
Keywords:high throughputsequencing;eukaryote;
prokaryote
收稿日期:
2020 06 08基金项目:2018年度江苏省环境监测科研基金项目(1810
)作者简介:景 
明(1986-
),男,江苏苏州人,硕士,工程师,研究方向为环境生态监测。Tel.:0512 68338041;E mail:jing_m2012@sina.cn通信作者:李继影(1980-
),女,安徽阜南人,高级工程师,研究方向为湖泊水生态。Tel.:0512 68338012;E mail:
li_jiying@qq.com
0 
河水污染
 言  张家港位于江苏省苏州市,区域范围内水系发达,随着区域经济发展,人类活动一定程度上影响了河道
水质。本研究主要针对张家港某生态污染河道河水呈赤红,经显微镜初步鉴定发现该现象是由微生物爆发导致的生态污染事件。微生物是指示水生态环境变化的敏感指标[1
]。目前,微生物分类主要依靠光学根
据细胞的形态结构特征确定其种类[2
]。但显微镜检测适用于鉴定细胞体积较大,具有显著形态学的大型细胞,对于部分个体较小的细胞无法有效鉴定[3
]。
 第
40卷 随着基因测序技术不断地发展与普及,大数据分
析能力不断提高,基因组研究手段正越来越多地应用
于鉴定形态学难以辨别的物种[4
]。原核生物和真核
生物的基因组都含有保守区和可变区,其中可变区能够区分物种之间的差异,适用于进行测序和分类鉴定,是研究环境样品中微生物物种落的重要手段。目前高通量测序技术已经广泛应用于研究蓝藻水华的形成及其爆发原因的机理。本研究通过Illumina测序对该河道的原核生物及真核生物分别进行测序,分析其落结构及成因。
1 
材料与方法
1.1 
样品采集及前处理
本次采样地点为张家港污染河道(坐标为
120°38′15.43E,31°48′32.11N),采样时间2019 12
23,每个采样点使用垂直采样器在水面下0.5m处收集5L水样,储存在塑料瓶内,置于4℃冷
藏箱后运回实验室,在24h内完成水样中微生物的富集。采用标准方法[5
]测定以下常规水质指标:叶绿素a、氨氮、总磷、总氮、化学需氧量、度和pH。样品浓缩采用膜过滤法,将样品通过醋酸纤维素滤膜(孔径0.45μm、直径50mm),用PBS(0.01mol/L)对滤膜进行洗脱。然后将洗脱液进行高速离心浓缩(14000r/min,10min)
,最后去掉上层清液,用1mLPBS进行洗脱,然后进行DNA提取。1.2 DNA提取及高通量测序
DNA提取采用CTAB[6]法标准化操作流程,然后采用琼脂糖凝胶电泳进行DNA浓度和纯度检测。最后取适量的样本DNA于
离心管中,使用无菌水稀释样本至1ng/μL作为基因组DNA模
板。根据待测序区域的选择,使用带Barcode的特异引物(NewEngland
Biolabs公司的Phusion High FidelityPCRMasterMix
withGCBuffer)和高效高保真酶进行PCR,
确保扩增效率和准确性。本研究选择针对原核生物16SrRNA
基因中高度变化的V4区域[7
]以及真核生物18SrRNA
的V4区域[8
]的引物作为扩增子,详见表1。PCR的反应条件为:95℃5min,94℃60s,57℃45s,72℃60s,共34个循环,最后72℃终延伸10min,16℃5
min。PCR产物通过2%的琼脂凝胶电泳检测。根据
PCR产物的浓度对样品进行等量混合,然后使用1×TAE浓度2%的琼脂糖胶电泳对PCR产物进行纯化,将剪切得到的目标条带使用试剂盒回收产物,试剂盒为ThermoScientific公司GeneJET胶回收试剂盒。  使用建库试剂盒(
IonPlusFragmentLibraryKit48rxns,Thermofisher)
进行文库的构建,构建好的文库经过Qubit定
量和文库检测合格后,进行上机测序(IonS5TM
XL,Thermofisher)
1 PCR1.3 
使用Cutadapt[9
]先对原始数据的reads进行低质量部分剪切,去掉干扰数据,再根据Barcode从处理后的数据中拆分出各样品数据,然后进行去除嵌合体序列处理,即对截去Barcode和
引物序列初步质控得到原始数据
Reads序列[10
]通过与物种注释数据库进行比对检测嵌合体序列,并最终去除其中的嵌合体序列[11
],得到最终的有效数据(CleanReads)
。用
Uparse软件[12
]对样品进行聚类分析,以97%的一致性(Identity)将序列聚类成为OTUs(Operational
TaxonomicUnits),选取OTUs中出现频数最高的序列作为OTUs的代表序列。对OTUs序列进行物种注释,用RDPClassifier方法与Silva132数据库进行物种注释分析(设定阈值为0.6~1.0)。使用MUSCLE软件
进行快速多序列比对,得到所有OTUs序
列的系统发生关系。最后,以样品中数据量最少为标准对测序数据进行均一化处理,后续的Alpha多样性分析和Beta多样性分析都是基于均一化处理后的数据。
2 
结果与分析
2.1 
测序数据分析及多样性指数本次扩增的第1轮PCR扩增产物如图1所示,所扩增的基因片段条带单一,长度分别约为300bp及
350bp,
真核生物的生物量显著高于原核生物。(
a)16SrRNA 
b)18SrRNA图1 扩增产物电泳图
测序共获得16SrRNA基因V4可变区的88968
条原始序列,通过数据筛选后共得到80148个高质量序列;18SrRNA基因V4可变区的88510条原始序列,通过数据筛选后,共得到87841个
高质量序列。8
1期景 明,等:高通量测序技术在水生态污染监测中的应用 将序列进行随机抽样,并通过对所抽到的序列数和代表的
OTU数来构建稀释性曲线,如图2所示。由图可知,原核生物得到的注释为548种OTU,
真核生物得到的注释为33种,两个样品的稀释曲线趋于平坦,说明原核生物与真核生物的测序结果都较为完整,原核生物的物种分布较为均匀种组成较为丰富;真核生物的种组成比较单一
(a)原核生
(b)真核生物
2 OTUs
稀释曲线
对不同样本在97%一致性阈值下的Alpha
Diversity分析指数进行统计,结果见表2(均一化时选取的数据量:cutoff=65446(原核生物);cutoff=
79939(真核生物))。Shannon Wiener指数和Simpson指数表示个体分布的
均匀度,各物种之间个体分布越均匀,指数越高,如果每一个体都属于不同的物种,指数则达到极值,由表可知真核生物的物种组成比原核生物单一。
表2 AlphaIndices
统计表
2.2 原核生物的多样性和落结构分析对样品中的原核生物OTUs进行注释后发现,本次样品原核生物多样性丰富,涵盖了22门74种。其中
,变形菌门(Proteobacteria)21种,占54.5%;拟杆菌门
(Bacteroidetes)24种,占15.77%;放线菌门(Actinobacteria)
1种
,占12.37%:厚壁菌门(Firmicutes)15种,占8.19%,如图3所示。原核生物OTUs前6的属占总测序量的百分比见表3,6个属
OTUs总百分比仅为32.18%,表明原核生物种类分布较为均匀。
图3 原核生物门水平上的物种相对丰度饼状图
3 62.3 
真核生物的多样性和落结构分析
所得到真核生物的OTUs序列进行物种注释后发现,本次采集的样品中真核生物种类比较单一,共发现4个门12种,其中金藻门(Chrysophyta)占93.95%,
如图4所示,样品中核心OTU是Pedospumella
encystans占93.94%,3种微生物占真核生物的95%以上(见表4)
。图4 真核生物门水平上的物种相对丰度饼状图
 第
40卷 
4 3 
5状,游动速度较快。且棒状藻体中含5~10个颗粒状素体。几乎所有的可见细胞均为该微生物,细胞颜为紫红,水体中大量出现,导致河道水面呈现为紫红,形态学鉴定无法准确判断该微生物种类。  根据高通量测序结果中真核生物和原核生物生物量以及物种丰度的
区别,导致本次河道生态污染的微
生物为Pedospumellaencystans。筛
选真核生物中特别关注的物种(相对丰度前
10的种)进行物种分类树统计,进行展示,其物种分类树如图6所示。
图5 显微镜观察结果
图6 样
本中特定物种分类树3 
 论3.1 
高通量测序技术的应用该河道微生物的优势种个体较小,通过显微镜难以准确判断其类别,而高通量测序技术具有能够检测
数量较少、体积微小且不宜培养的微生物的能力
[13
],用于检测环境样品中微生物落结构的生物多样性[4
]。目前江苏省浮游藻类监测主要建立在形态学鉴定的基础上,本研究采用高通量测序技术对河道原核生物及真核生物的落组成进行研究,同时结合形态学检测结果,发现引起本次河道生态污染的微生物为一种体积较小,难以通过传统显微镜鉴定种类的真核浮游藻类。高通量测序技术可以弥补显微镜观察中
难以判断的物种信息,从而提高了浮游藻类监测的综
合能力[14
]。3.2 水环境因子对生物落结构的影响本研究中,河道的水质参数如表5所示。由表5
可知,该河道水体属于劣5类水体,相比其他指标化学需氧量浓度较低(
5类水40mg/L),氨氮(5类水2.0mg/L)与总磷(5类水0.4mg/L)浓度较高,均超过5
类水要求5倍以上,可生化性较差(可生化处理的水质推荐的碳氮比为100∶5,该河道碳氮比为1∶1,推测有工业废水排入),对活性微生物繁殖造成影响,致使水体的自净能力得到影响,当有少数微生物能够适应该水质条件时,该微生物易成为优势种。
 第1期
景 明,等:高通量测序技术在水生态污染监测中的应用 
5 
3.3 优势种Pedospumellaencystans
金藻是生态学和生态生理学的原生模式种,一部分种类是水体中的初级生产力,另一部分是水体中部分细菌的捕食者,通常最佳生长条件为温度较低,透明度较好且有机物含量较低的水体[15],在初春、晚秋或冬季多有发现,多数存活于水体的中下层,有研究表明在冰面下某些金藻种类仍能存活[16]。
近几年,系统发育树和大系统分类学显示了金藻门所包含种类的复杂性,许多种类为光合自养和异养生物之间的过渡[17],其中异养型金藻与自养型金藻相比,具有较小的基因组和细胞体积,相对表面积较大,能够高效利用环境中的营养盐且所需代谢时间短于其它藻类。Pedospumellaencystans属于金藻纲(Chrysophyceae),金藻目(Chromulinales),
Chromulinales科,Pedospumella属,属于一种异养型微型浮游藻类[15],细胞无细胞壁,仅含有一条鞭毛,
Pedospumellaencystans在中国发现较少,在已有研究中有报道,当水质的可生化性较差时,自养型微生物难以
生存的条件下该种类易成为优势种[18]。
4 结 论
本研究结果表明,运用高通量测序技术结合形态学鉴定能够揭示水环境中微生物的多样性及落组成,得到更客观、更准确的监测结果,在水生态监测和水质评价方面具有较好的应用潜力。
(1)对张家港生态污染河道进行形态学鉴定,落组成单一,可见的微生物基本为同一种类。(2)用高通量测序技术对18SrRNA基因V4区测序,结果显示数量较多且落组成单一,
Pedospumellaencystans占93.94%。用高通量测序技术对16SrRNA基因V4区测序,结果显示数量相对较
少且落分布较为均匀,由此得出主要的微生物为
Pedospumellaencystans。
(3)对该河道水质进行分析,推测该河道有工业废水混入,可生化性较差,同时冬季河道水质不利于常规藻类生长,使得Pedospumellaencystans成为优势种,导致了本次河道生态污染事件。
参考文献(References):
[1] 耿文华,陈继淼,冯剑丰,等.辽河浮游藻类落及生物多样性基准验证[J].中国环境科学,2014,34(1):239 245.[2] 杨萍萍.浮游藻类显微镜计数方法适用范围探究[J].资源节约与环保,2017(7):101 102.[3] MoralesEA,SiverPA,TrainorFR.Identificationofdiatoms(Bacillariophyceae)duringecologicalassessments:comparison
betweenlightmicroscopyandscanningelectronmicroscopy
techniques[J].ProceedingsoftheAcademyofNaturalSciencesof
Philadelphia,2001,151(1):95 103.
[4] 施军琼,张 明,杨燕君,等.基于高通量测序探讨大宁河不同水华期真核浮游生物落组成[J].西南大学学报(自然科学版),2020,42(2):1 7.
[5] 国家环境保护总局.水和废水监测分析方法[M].4版.北京:中国环境科学出版社,2002:88 284.
[6] 余智城,何雪娇,林秀香,等.柚叶片DNA提取方法比较[J].福建热作科技,2019,44(3):26 28.
[7] CaporasoJG,LauberCL,WaltersWA,etal.Globalpatternsof
16SrRNAdiversityatadepthofmillionsofsequencespersample
[J].ProceedingsoftheNationalAcademyofSciences,2011,108
(Sup.1):4516 4522.
[8] YoussefN,SheikCS,KrumholzLR,etal.Comparisonofspeciesrichnessestimatesobtainedusingnearlycompletefragmentsand
simulatedpyrosequencing generatedfragmentsin16SrRNAgene
basedenvironmentalsurveys.[J].Applied&Environmental
Microbiology,2009,75(16):5227.
[9] LiB,ZhangX,GuoF,etal.Characterizationoftetracyclineresistantbacterialcommunityinsalineactivatedsludgeusingbatch
stressincubationwithhigh throughputsequencinganalysis[J].
WaterResearch,2013,47(13):4207 4216.
[10] LundbergDS,YourstoneS,MieczkowskiP,etal.Practicalinnovationsforhigh throughputampliconsequencing[J].Nature
Methods,2013,10(10):999.
[11] JiaoS,LiuZ,LinY,etal.Bacterialcommunitiesinoilcontaminatedsoils:Biogeographyandco occurrencepatterns[J].
SoilBiologyandBiochemistry,2016,98:64 73.
[12] QinJ,LiY,CaiZ,etal.Ametagenome wideassociationstudyofgutmicrobiotaintype2diabetes[J].Nature,2012,490(7418):
55 60.
[13] 宋 伦,吴 景,李 楠,等.基于高通量测序的真核微藻多样性质控分析[J].水产科学,2019,38(6):804 812.
[14] 王靖淇,王书平,张 远,等.高通量测序技术研究辽河真核浮游藻类的落结构特征[J].环境科学,2017,38(4):
1403 1413.
[15] BeisserD,GraupnerN,BockC,etal.Comprehensivetranscriptomeanalysisprovidesnewinsightsintonutritionalstrategies
andphylogeneticrelationshipsofchrysophytes[J].PeerJournal,
2017(5):2832.
[16] 谢树莲,冯 佳.中国淡水金藻门植物研究进展[J].世界科技研究与发展,2007,29(6):1 6.
[17] OlefeldJL,MajdaS,AlbachDC,etal.Genomesizeofchrysophytesvarieswithcellsizeandnutritionalmode[J].
OrganismsDiversity&Evolution,2018,18(2):163 173.
[18] 梁采莹.环北部湾表层沉积物生源要素分布及萌发的浮游藻类落结构研究[D].广州:
暨南大学,2018.
13