医药化工废水厂中微生物群落结构分析Word文档格式.docx

医药化工废水厂中微生物群落结构分析Word文档格式.docx

《医药化工废水厂中微生物群落结构分析Word文档格式.docx》由会员分享，可在线阅读，更多相关《医药化工废水厂中微生物群落结构分析Word文档格式.docx（8页珍藏版）》请在冰豆网上搜索。

SpinKitforSoil（MPBiomedicals,CA,USA）试剂盒根据产品操作手册提取样品基因组DNA.采用1%琼脂凝胶电泳检测DNA样品的质量和完整性.用NanoDropTM2000分光光度计（ThermoFisherScientific,MA,USA）和Qubit®

2.0荧光光谱仪（ThermoFisherScientific,MA,USA）检测DNA的产量和纯度.基因组DNA样品用干冰保藏并立即送往上海美吉生物医药科技有限公司进行测序.

　　1.3ARGs生物信息学分析

　　采用IlluminaHiSeqXten平台进行PE150双端测序.首先对原始测序数据进行质量控制,去除碱基数大于40bp但是其质量分数低于38的序列,去掉歧义核苷酸超过10个的序列,去掉与adapter序列有15bp重叠的序列.采用SOAPdenovo软件组装cleanreads.当选择最高的scaffoldN50value来进行进一步分析时,选择不同的K-mer值（49、55和59）用于使用于不同的-d1、-M3、-u和-F参数.将得到的scaffold片段在N段打断,进而得到scaftig片段.仅仅大于500bp的scaftig片段被保留以作进一步分析.采用默认参数的MetaGeneMark软件从组装好的scaftig片段来预测开放阅读框（ORFs）,采用CD-HIT来获得非冗余基因.采用SOAPaligner软件将保留的cleanread与基因库进行比对,以计算匹配的数目.通过去除含有少于3个mappedreads映射片段的基因类别来获得unigene.用DIAMOND软件将unigene与NCBI-NR数据库进行比对,设置blastp参数e-value≤1e-5,并且采用MEGAN软件LCA算法进行注释.采用DIAMOND软件（blastp,evalue≤1e-5）将unigene对应的蛋白序列在KEGG数据库、eggNOG数据库和CAZy数据库中进行比对.利用BLAST软件（blastp,evalue≤1e-5）将unigene在CARD数据库进行比对,进行ARGs注释.采用Cytoscape软件进行网络分析,以研究ARGs与微生物分类的关系.具体联系污水宝或参见更多相关技术文档。

　　2结果与讨论

　　2.1医药化工污水厂中微生物群落结构

　　在该好氧活性污泥样品中,经过质量控制后,共得到reads数目109523076条,并由此组装得到55255个scaftig,并预测了105892个unigene.利用DIAMOND软件将unigene在NCBI-NR数据库中进行比对,后续用MEGAN软件进行注释,统计样本在门水平（表1）及属水平（图1）上的前10个丰度最大的群落组成.

表1废水厂中10个最丰富门的相对丰度

图1废水厂中10个最丰富属的相对丰度

　　根据分类结果,细菌域（Bacteria）是该医药化工园区废水厂中主要的微生物类群.并由表1可知,该废水厂丰度在前10的优势菌门是Proteobacteria、Chloroflexi、Bacteroidetes、Acidobacteria、Actinobacteria、Cyanobacteria、Deinococcus-Thermus、Firmicutes、Planctomycetes和CandidatusSaacharibacteria.其中最优势的菌门是Proteobacteria,这与其他关于污水厂报道相符合,丰度接近50%;

其次是Chloroflexi,丰度约占10%.Chloroflexi菌株可以在极端环境中生存,比如热泉和高盐环境,并且也经常从污水厂环境中分离得到,Chloroflexi类群中的丝状菌是污水处理厂活性污泥中的常见组成,其大量出现会产生污泥膨胀现象.细菌菌群丰度居第三的是Bacteroidetes,约占5%;

其他7个优势菌群的丰度在0.5%~3%之间.

　　在该混合污泥样品中,共检测到265个属,丰度在前10的优势属为Hyphomicrobium、Ardenticatena、Paracoccus、Thauera、Leptolyngbya、Filomicrobium、Methyloversatilis、Truepera、Ottowia和Aromatoleum,其中Hyphomicrobium是最主要的类群,属于革兰氏阴性菌,占据约12%;

Ardenticatena、Paracoccus和Thauera这3个优势属的丰度相似,均在4.5%左右,其中Ardenticatena丰度稍高,为4.7%;

其他6个优势属的丰度在0.2%~1.0%之间（图1）.在属水平上,Hyphomicrobium是最优势属,而在种水平上,Hyphomicrobiumzavarzinii是最优势种.根据报道,Hyphomicrobium属类菌株也经常出现在污水环境中,属于限制性的兼性甲基氧化菌,能够利用甲醇、甲基化胺或甲酸盐等一碳化合物作为唯一碳源生长.而Hyphomicrobiumzavarzinii能通过甲醇脱氢酶的过量表达,在污水甲基营养型代谢中起着关键作用.

　　相比之下,城市污水处理厂（wastewatertreatmentplants,WWTPs）中处理的废水大部分为生活污水,可生化性好,废水中的抗生素一般来源于医院排放的医疗废水,含量较少,尤其经大量其他废水稀释之后浓度一般为痕量.从WWTPs中分析到的主要细菌种类属于粪便污染的常见指标：

大肠杆菌、总大肠菌群和肠球菌.监测到的主要耐药菌包括耐甲氧的Staphylococcusaureus、耐万古霉素的Enterococcusspp.及革兰氏阴性菌（如Enterobacteria、Pseudomonads和Acinetobacter）,这些菌对氟喹诺酮类、碳青霉烯类和β-内酰胺类等抗生素具有抗性.医药化工废水处理厂的抗生素浓度高于WWTPs的10~100倍,因此其抗性菌种的种类和耐药性也一般强于WWTPs.

　　2.2医药化工废水厂中微生物群落功能注释

　　用DIAMOND软件将unigene在KEGG数据库进行比对,共注释了54564个unigene,结果如图2所示.从中可知,该微生物菌群的功能通路相关的基因数量依次为：

代谢>

遗传信息处理>

环境信息处理>

人类疾病>

细胞过程>

生物系统,其中最主要的功能通路是代谢相关,其次是遗传信息处理,第三是环境信息处理.属于代谢类别的功能分类中丰度依次为：

碳代谢>

氨基酸代谢>

能量代谢>

核苷酸代谢>

辅因子和维生素的代谢>

异生素生物降解与新陈代谢>

脂质代谢>

萜类和聚酮化合物代谢>

其他氨基酸代谢>

其他次生代谢物生物合成>

多糖生物合成与代谢,其中碳代谢、氨基酸代谢和能量代谢的基因数量尤为丰富,分别为7788、6155和5067个.属于遗传信息处理类别的功能分类丰度依次为：

复制与修复>

折叠、分类与降解>

翻译>

转录,其中复制与修复是最主要部分,基因数量为2535个.该微生物菌群功能通路中还有一部分的环境信息处理类别的功能,其分类类别丰度依次为：

膜转运>

信号转导>

信号分子与相互作用,其中最主要的是膜转运,基因数据为2678个,其他两个类别的基因合计100个,数目甚少.属于人类疾病类别丰度依次为：

传染性疾病>

神经退行性疾病>

内分泌和代谢疾病>

癌症>

物质依赖>

免疫疾病>

心血管疾病,其中细菌传染性疾病相关的基因数目尤其突出,为517个,是其他人类疾病基因总数的约1.6倍.而在传染病疾病中,细菌感染性疾病相关的基因占据了大部分,约80%.

图2基于KEGG数据库注释的废水厂微生物群落功能基因分布 

　　在该废水厂中细菌传染疾病涉及到9类（图3）,其丰度依次为军团杆菌病>

结核病>

沙门氏菌感染>

幽门螺杆菌感染的上皮细胞信号传导>

金黄色葡萄球菌感染>

上皮细胞的细菌入侵>

百日咳>

致病性大肠杆菌感染>

霍乱弧菌感染,其中最主要的传染病为军团杆菌病,其基因数目为136个;

其次是结核病,基因数目为108个;

沙门氏菌感染和幽门螺杆菌感染的上皮细胞信号传导,基因数目均在60个左右;

其他疾病的基因数目在15个以下.

图3废水厂中细菌感染疾病相关的基因数量

　　根据功能注释结果,该废水厂菌群中最主要的功能通路是代谢相关,并且还存在许多与人类疾病相关的基因,其中主要是细菌感染性疾病.到目前为止,全球新出现的传染病事件的发生率有所增加,废水厂等特殊环境可以作为许多病原体的天然宿主对人类健康存在潜在威胁.

　　2.3医药化工废水厂中ARGs及耐药机制

　　用BLAST软件将unigene在抗性基因数据库CARD中进行比对,结果如图4所示.经过氨基酸大于50%标准控制后,共注释到74类抗生素功能基因,涉及包括β-内酰胺类、头孢菌素类、氨基糖苷类、大环内酯类、多肽类、喹诺酮类及磺胺类等多种抗生素抗性基因.

图4污水厂中不同ARGs的相对丰度

　　在这74类ARGs中sav1866、dfrE和mfd是最主要ARGs.sav1866基因属于多重耐药基因,其编码的蛋白是ABC膜转运蛋白,其耐药机制属于ABC外排泵.在细菌中,许多ABC家族转运蛋白都与抗生素耐受性有关.sav1866被检测为丰度最大的ARGs,推测这与该废水厂处理多种抗生素中间体废水及多种农药中间体废水密切相关.dfrE基因首次被发现在粪肠球菌中,其编码的蛋白是一种耐甲氧苄啶二氢叶酸还原酶（dihydrofolatereductases,DHFRs）,该种DHFRs具有高的甲氧苄啶抗性,其耐药机制是抗生素靶位点改变.mfd基因的产物Mfd是转录-修复偶联因子蛋白,属于DNA移位酶,也是DNA修复的必须因子,其耐药机制是抗生素靶位点保护.Mfd主要功能是促进DNA修复,但在某些情况下,如在复制-转录冲突区域和在稳定阶段的突变中,Mfd反而增加了突变,如Mfd促使空肠弯曲杆菌产生了对喹诺酮药物的抗性.在新的研究结果中,Mfd被认为是一种促进突变的进化因子,加速了抗生素耐药性的发展,是不同类群的细菌快速获得各种抗生素耐药性所必须的[35].由此推测,mfd基因的大量存在及积累可能会促进废水厂中微生物产生更多抗生素抗性.

　　统计这些ARGs的耐药机制,以调查该废水厂的耐药机制情况,如图5中（a）所示.在该医药化工废水好氧活性污泥处理系统中,微生物存在耐药的机制有5种,分别是抗生素失活（8%）、抗生素靶位点改变（4%）、抗生素靶位点保护（7%）、抗生素靶位点替换（11%）和抗生素特异的外排泵（70%）,其中抗生素特异的外排泵是主要的耐药机制,占比50%以上.在所研究的菌群中的外排泵组成如图5（b）所以,主要由ABC外排泵（18%）、MFS外排泵（10%）、RND外排泵（34%）和SMR外排泵（3%）组成,其中RND外排泵是该微生物菌群抗生素外排泵耐药机制的最主要类型,其次是ABC外排泵.据报道,细菌细胞膜上的主动外排泵可将抗菌药物排至细菌外,从而阻碍抗菌药物与细菌内的靶点结合,引起耐药.并且,外排泵在细菌中广泛存在,在介导细菌多重耐药中起着重要作用.本研究表明,氯霉素、四环素、磺胺类、氟喹诺酮类和大环内酯类抗生素均能通过细菌的主动外排机制而被排出细胞外.而该研究的医药化工污水系统所具有的外排泵中,RND外排泵是主要机制.RND外排泵的外排底物比较广泛,包括多种抗菌药物、消毒剂和染料等,由此可推测该菌群对多种抗生素具有抗性.结合图4,属于RND外排泵的ARGs有smeD、mdtB、mexJ、smeR、cmeB和mexQ等,涉及氨基香豆素抗生素、四环素抗生素、大环内酯类抗菌素、氨基糖苷抗生素、三氯沙抗生素和氟喹诺酮类抗生素、头孢菌素等ARGs;

属于ABC外排泵的ARGs有lmrD、otrC、novA、macB和sav1866,涉及林可酰胺类抗生素、四环素类抗生素和大环内酯类抗生素的耐药.

（a）微生物的抗生素耐药机制;

（b）抗生素特异的外排泵组成

图5污水厂中抗生素耐药机制类型占比

　　2.4废水厂中ARGs与微生物类群的关联

　　对废水厂中ARGs与微生物类群的共存模式进行网络分析（图6）,结果显示有44类ARGs与污水厂中属级的微生物分类群高度相关,并且主要类群的属都是多重耐药菌株,同时携带多种ARGs,涉及到4大类ARGs耐药机制.Hyphomicrobium属和Thauera属均可能携带有近20个ARGs,Paracoccus属和Methyloversatilis属可能携带有近10个ARGs.如Hyphomicrobium属中可能携带的抗性基因有dfrA26、mdsB、qepA、tetT、mexQ、mexF、mexL、rosB、sul2、dfrE、mfd和bacA等;

Thauera属则可能携带有mdtD、bacA、mexJ、tet36、oprM、acrB和novA等.

图6污水厂中ARGs与微生物类群的相关性网络分析

　　3结论

（1）该医药化工废水处理厂的菌群组成结构中,细菌域是主要微生物类类群,并且最优势的菌门是Proteobacteria,其次是Chloroflexi.最优势属是Hyphomicrobium,最优势种为Hyphomicrobiumzavarzinii.

（2）该医药化工废水厂中微生物群落最主要的功能通路是代谢相关,并还存在许多与人类疾病相关的基因,其中主要是细菌感染性疾病.

　　（3）该医药化工废水厂中ARGs共检测到74类ARGs,涉及包括β-内酰胺类、头孢菌素类、氨基糖苷类、大环内酯类、多肽类、喹诺酮类、磺胺类等多种抗生素抗性基因,其中sav1866、dfrE和mfd是最主要ARGs,而最主要的Hyphomicrobium属也同时携带这3种主要的ARGs.该废水厂菌群的耐药机制主要是抗生素特异的外排泵,其中RND外排泵是主要机制.

　　（4）该废水厂活性污泥中微生物属级分类与ARGs高度相关,并且主要类群的属都是多重耐药菌株.