制药废水厂微生物群落和多种抗性基因相关性分析 袁立霞 罗晓 张文丽 蒋永丰Summary：为了研究多种生物处理工艺制药废水中微生物群落结构和抗性基因的分布特征、扩增情况及其相关性，采用Miseq高通量测序分析技术和荧光定量PCR技术对制药废水厂中活性污泥进行检测荧光定量结果表明：sul1，sul2，tetO，tetQ，tetW，OXA-1和可移动遗传元件int1在制药废水厂中各个阶段均能被广泛地检测到，总抗性基因浓度范围为3.09×108～2.26×109 copies/g（干重），基因总浓度上升了7.3倍Miseq测序结果表明：制药废水中主要优势菌门为Proteobacteria，Bacteroidetes，Firmicutes，Thermus和Gemmatimonadetes等门，其平均总相对丰度比例占到81.05%;冗余分析显示，Aeromicrobium与sul2呈较高程度的正相关性，可能是sul2在微生物群落中存在的可能的主要菌群;Rhodovulum和Rhodospirillaceae菌属与OXA-1和tetQ呈较高程度的正相关性，这些菌属是这些ARGs分布的主要菌属;Acinetobacter与tetO和tetW呈现较高程度的正相关性，可能是tetO和tetW在微生物群落中存在的可能的主要菌属。
因此，相关抗性基因的增值和分布与相关特定菌属有关，可以通过控制相关菌属的丰度来消减工业废水厂中的抗性基因Key：微生物生态学;抗性基因;荧光定量PCR;Miseq;微生物群落结构：X703文献标志码：AAbstract： In order to study the distribution， amplification and correlation of microbial community structure and resistance genes in pharmaceutical wastewater， a variety of biological treatment processes are used to treat pharmaceutical wastewater. Combining Miseq high throughput sequencing analysis technology and fluorescence quantitative PCR technology， activated sludge in pharmaceutical wastewater plants is detected. The results show that sul1， sul2， tetO， tetQ， tetW， OXA-1 and INT1 are widely detected at all stages of pharmaceutical wastewater treatment plant. The total resistance gene concentration ranges from 3.09×108 to 2.26×109 copies/g （dry weight）， and the total gene concentration increases by 7.3 times. Miseq sequencing results show that the dominant bacteria in pharmaceutical wastewater are Proteobacteria， Bacteroidetes， Firmicutes， Thermus and Gemmatimonadetes， with an average relative abundance ratio of 81.05%. Redundancy analysis shows that Aeromicrobium is positively correlated with sul2 to a high degree， which might be sul2 existing in microbial communities. Among the possible major microflora， Rhodovulum and Rhodospirillaceae are highly positively correlated with OXA-1 and tetQ， which are the main microflora of ARGs distribution. Acinetobacter has a high degree of positive correlation with tetO and tetW， which may be the main possible microflora of tetO and tetW in microbial community. The results show that appreciation and distribution of antibiotic resistance gens are connected with concerned specific bacterial genus. The controlled abundance of concerned bacterial genus could subduce resistance gene in industrial wastewater factory.Keywords：microbial ecology; antibiotic resistance genes;fuorescence quantitative PCR;Miseq;microbial community structure由于抗生素廣泛应用于人类和动物（畜禽/鱼）的疾病防治或动物的促进生长，大量残留抗生素进入环境，致使微生物在持续抗生素选择下产生了耐药性，形成了耐药菌[1-2]，其风险比抗生素药物本身的污染风险更高[3]，关于其行为特点和传播途径的研究日益增加。
抗生素抗性细菌对抗生素的耐药性机理包括：细菌外膜不渗透性障碍、细菌外排泵系统、抗生素作用的靶位变化和抗生素的钝化失活等[4]抗性基因是抗性菌具有抗性的主要原因之一，它可以通过多种形式的可移动遗传元件如质粒、整合子、转座子、插入序列等，突破细菌的种属关系广泛传播，加重抗性污染FORSBERG等[5]证明人类病源菌能传播抗性到土著细菌中，并且由土著细菌携带的抗性基因能转导到病原菌的宿主细胞抗性能潜在地影响环境中的微生物群落结构并对公共健康产生威胁[6]，且碳青霉烯类抗生素耐药肠杆菌科细菌患病率上升迅速[7]，这将对环境和人类健康构成严重的威胁[6，8]分析荧光定量PCR技术为研究环境中抗性基因和抗性微生物分布特征提供了有力的手段[8]已有研究表明制药废水处理系统的进水中往往含有高浓度的抗生素及生产原料和中间体[9]，对生物处理单元的微生物群落形成高强度选择性压力，为ARGs在污水中的持久存在和传播扩散提供了有利条件，并且制药废水处理系统出水中ARGs的浓度高于市政污水处理系统出水[10-11]22种类型的ARGs能对5种类型的抗生素产生抗性，包括sul 1，sul 2和sul 3能对磺胺类抗生素产生抗性;tetO，tetQ，tetW，tetM和tetT对四环素类抗生素产生抗性;OXA-1和OXA-2能对β-内酰胺类抗生素产生抗性。
可移动遗传元件（例如1类整合子）有利于形成多重耐药性并且加速地表水、土壤[12-13]和污水处理厂[14]中抗性基因的传播，并且可移动遗传元件（1类整合子）是细胞间转移抗性的指标[15]本文主要研究多种抗性基因在制药废水系统生物处理单元中增值和传播以及微生物群落的动态变化及两者的相关性选取河北省某抗生素生产废水为研究对象，通过现场采样分析，重点考察污泥中β-内酰胺类、四环素类和磺胺类抗性基因及微生物群落结构在废水处理过程中的分布和相关性特征，揭示废水厂活性污泥中土著抗性基因的分布规律以及与微生物群落结构之间的互相作用，以期为评估废水处理厂中抗性基因环境风险提供依据，并为制定有效的ARGs（antibiotic resistance genes）环境污染控制方法提供理论基础1材料与方法1.1样品采集河北省某制药废水处理厂，该厂设计处理水量300 m3/d，主要采用生物处理法处理生产废水，采样点和废水处理厂工艺流程如图1所示样品采用500 mL聚乙烯瓶，于2018年1月，采集不同污水处理单元活性污泥样品，分别标记为J（综合调节池）、IC（厌氧塔）、SBR（间歇式曝气池）、A（A池）、O（O池）、C（污泥池）、CH（出水）取样结束后，于冰盒中运送回实验室，离心图1采样点和废水处理厂工艺流程Fig.1Sampling point and process flow of wastewater treatment plant（5 min，11 000 r/min）后称取5 g冷冻于-80 ℃冰箱中，以备DNA提取。
1.2水样分析按照《水和废水监测分析方法》[16]国家标准方法分析常规化学指标，测定水中CODCr和氨氮，测定结果见表11.5抗性基因的定量检测采用SYBR-Green 实时定量PCR方法对各基因进行定量分析，检测仪器为StepOne型荧光定量PCR仪（ABI，美国）PCR产物经过克隆测序确认后，使用生工质粒提取试剂盒SK1131从阳性克隆子中提取质粒，用作标准曲线使用NanoDrop微量分光光度计（Thermo Scientific，美国）测定质粒浓度制作标准曲线时按照10倍梯度浓度稀释构建好的各质粒，于90 μL稀释液中加入10 μL质粒，做4～6个点，通过预实验选取合适标准品用于制备标准曲线标准质粒、环境样品、阴性对照均做3个平行，取平均值进行计算质粒拷贝数换算公式为Q=[C×10-9×6.02×1023]/[N×660]，式中：Q为质粒拷贝数，copies/μL;C为质粒浓度，ng/μL;N为克隆产物碱基数，（g/mol）;M1为载体碱基数，2 692 bp（pMD 18-T）;M2为载体碱基数，3 023 bp（pGEMX-T Easy）荧光定量PCR反应体系为12.5 μL 2×SYBR，1.0 μL DNA（10 ng/μL），0.5 μL引物F（10 mol/L），0.5 μL引物R（10 mol/L），10.5 μL水，总体积25.0 μL。
荧光定量PCR反应程序为1）50 ℃，2 min;2）95 ℃，5 min;3）95 ℃，20 s;4）退火，30 s;5）72 ℃，31 s;6）Plate read，重复1）～5），39次重复;7）Melt-curve分析60 ℃ to 95 ℃每隔0.2℃采集一次荧光以生成溶解曲线，根据溶解曲线的变化检测扩增结果的特异性退火温度和反应时间根据引物不同进行调整1.6统计分析抗性基因浓度结果和群落结构使用OriginPro 8.6软件（Origin Lab Corporation，USA）进行分析;使用MOTHUR软件计算各个样本Alpha多样性指标，以反映本次测序深度、物种均匀性等;使用R软件对样本绘制热图并分析;ARGs与群落结构相关性使用CANOCO 5.0（Microcomputer Power ，Ithaca，NY）软件分析，基因相关性使用SPSS 19.0软件进行相关性分析2结果讨论与分析2.1微生物菌群耐药基因定量检测利用荧光定量PCR的方法检测抗性基因的分布特征、变化情况和相对表达量，水处理阶段ARGs分布结果如图2所示2.2微生物群落多样性分析2。