式中，N代表阳性克隆子总数，n1代表仅含单个克隆子的数量.

　　所有测序结果在GenBank数据库进行Blast比对后下载同源性序列，系统发育树的绘制借助MEGA5.0软件完成.

　　3 结果与讨论

　　3.1 可培养耐红霉素细菌的计数与分离鉴定

　　根据抗菌药物敏感性实验执行标准，红霉素的最低抑菌浓度(MIC)为8 μg · mL-1，本研究选择略高于其MIC值的抗生素浓度，即红霉素的终浓度为16 μg · mL-1.利用平板菌落计数法对抗性细菌的计数结果为9.86×108 CFU · g-1.将平板上形态各异的耐药菌株进行反复纯化分离后，得到了3株独立的菌株，其各自的形态如表 1中所示.

表1 可培养耐红霉素菌菌落形态

　　分别将3株菌株的16S rDNA序列在NCBI核酸数据库中进行BLAST同源性比对，确定与其最接近的种系群体.结合各菌株形态学特征及16S rDNA 序列分析，初步鉴定为赖氨酸芽孢杆菌(Lysinibacillus sp.)、土壤芽孢杆菌(Solibacillus silvestris)以及蜡样芽孢杆菌(Bacillus cereus)，各自在GenBank中的登录号分为KP694229、KP694230、KP694231，且上述3种可培养菌均属于芽孢杆菌纲.需要注意的是，蜡样芽孢杆菌在自然界分布广泛，是一种食源性致病菌，容易引起食物中毒，对红霉素极度敏感，对青霉素耐药.然而，本文的耐药实验发现，从底泥中分离出的蜡样芽孢杆菌GD2C对红霉素已表现出一定的耐药性，这可能是菌株分离自长期受低浓度红霉素污染的环境，这意味着低浓度红霉素的诱导作用可使蜡样芽孢杆菌具备一定的耐药性.

　　3.2 16S rDNA克隆文库的分析

　　随机挑取83个阳性克隆子进行测序(序列长均为1500 bp左右)，将测序结果在NCBI的核酸数据库中进行Blast比对，将序列比对结果相似度大于等于97%的定义为一个操作单元(Operational Taxonomic Unit，OTU)，共得到6个OTU，其中有2个OTU类型代表单克隆，最多的OTU类型包括36个克隆子.随后以文库覆盖率C及稀有性曲线来评价所构建的文库，发现该文库覆盖率为97.6%，且稀有性曲线趋于平缓，这意味着所构建的克隆文库库容较大，所挑取的克隆子数量能够较为完整的反映样品中微生物的多样性.上述两项分析表明该克隆文库所含83个克隆子能够较全面地反映底泥中耐红霉素细菌种群结构的多样性.

 图1 底泥中耐红霉素细菌16S rDNA克隆文库的稀缺性曲线 

　　样品中16S rDNA克隆文库结果如表 2所示，其中4个OTU的克隆子序列与NCBI核酸数据库中已知细菌的16S rDNA序列相似度在97%以上，2个OTU的相似度低于97%.由比对结果可知，长期受低浓度红霉素污染的底泥中，耐红霉素的细菌类群可分为三大类群，即：未获培养细菌(Uncultured bacterium)、芽孢杆菌纲细菌(Bacilli)和梭菌纲细菌(Clostridia)，各涉及到2个OTU(54个克隆子)、3个OTU(28个克隆子)和1个OTU(1个克隆子).底泥样品中耐红霉素微生物中的优势菌群为某种未获培养的细菌，在文库中占43.37%，其次是苏云金芽孢杆菌，占文库的24.10%.

表2 底泥中耐红霉素细菌16S rDNA克隆文库结果 

　　苏云金芽孢杆菌是一种有益杆菌，多用做杀虫剂，目前对于其耐药性的报道较少，但从文库结果来看，苏云金芽孢杆菌已对红霉素产生了一定的耐药性;另外，红霉素原本对梭状芽孢杆菌有较强的抑制作用，在本研究中也表现出了一定的耐药性，这可能也是因为环境中残留的低浓度红霉素对其诱导，使其对红霉素产生了一定的抗性，这也进一步说明了环境中残留抗生素存在潜在影响.由克隆文库结果可以看出，在已知菌属中，芽孢杆菌属所占比例最大，这与上述所筛选分离出的可培养的耐药菌均为芽孢杆菌属细菌的结果是互相吻合的.

　　另外，构建的16S rDNA文库中有25.3%(21个克隆子)的细菌与GenBank中已知细菌的16S rDNA序列同源性低于97%，而在序列同源性研究中，一般认为细菌16S rDNA序列同源性低于97%即属于不同的种，这说明底泥中蕴藏着许多未知的耐药细菌，有待于进一步深入研究.

　　3.3 底泥样品16S rDNA系统发育分析

　　分子系统发育分析，特别是对细菌16S rDNA基因的同源性进行分析有助于了解微生物群落结构的多样性.将表 1和表 2中的3株单菌和6个OTU的碱基序列及其相应的同源性序列进行聚类分析后，以嗜热链球菌Streptococcus thermophilus(NR042778)为外类群，绘制出底泥中耐红霉素细菌系统发育树，见图 2.

　　由图 2可见，土壤芽孢杆菌GD2B的碱基序列与OTU4的碱基序列之间相似度较高，且与OTU1存在较好的亲缘关系.因而将OTU1和OTU4的序列进一步在RDP数据库中进行比对，发现其最大相似菌株Uncultured bacterium clone(DQ337010)和Uncultured bacterium clone(AB506332)可以进一步划分到某未分类的动性球菌属(Unclassified Planococcaceae)，而土壤芽孢杆菌和动性球菌属均属于动球菌科(Planococcaceae)，所以OTU1、GD2B以及OTU4在系统发育树中聚集到一起.蜡样芽孢杆菌GD2C与OTU5的最大相似菌株苏云金芽孢杆菌存在较好的亲缘关系，这是由于苏云金芽孢杆菌和蜡样芽孢杆菌均属于蜡样杆菌属，两者之间具有高度的同源性.

 图2 底泥中耐红霉素细菌系统发育树 

　　红霉素作为第一代大环内酯类抗生素，自1952年应用于临床以来，一直被人们长时间大剂量使用，使得多种原本敏感的微生物对其产生了抗性(李喆宇等，2013)，比如肺炎链球菌、金黄色葡萄球菌、粪肠球菌等.而本研究发现的芽孢杆菌，相关报道尚不多见.于雅琼等(2007)发现，分离自土壤、饲料和肥料的芽孢杆菌对大环内酯类抗生素表现出一定的敏感性.本研究显示，经过为期1年的低浓度红霉素胁迫，芽孢杆菌对红霉素已表现出一定的耐药性，3株可培养的红霉素耐药菌及红霉素耐药菌克隆文库中有33.73%均属于芽孢杆菌.这可能是由于：①不同剂量使用大环内酯类药物，在机体内存在诱导出耐药株的可能性(刘禧杰，2008);②微生物细胞可以通过基因的水平转移或者可移动遗传因子获得抵抗力(李晶，2013)，通过抗性基因的水平转移获得耐受红霉素基因，从而获得对红霉素的抵抗力.当细菌具有抗生素抗性后，可能会对细菌代谢产生特定的改变(如抗重金属污染等特性)，甚至可使细菌在某些环境中更具有生存优势;而存在于基因转移单位的抗性基因，则可以自我复制，从而一直存在于微生物群落里(杨颖，2010)，这必将引发巨大的潜在危害，并威胁到环境及人类健康.

　　4 结论

　　1)用传统方法从底泥中筛选出的3株芽孢杆菌纲耐红霉素菌株，结合各菌株形态学特征及16S rDNA序列分析，初步鉴定为赖氨酸芽孢杆菌(Lysinibacillus sp.)、土壤芽孢杆菌(Solibacillus silvestris)以及蜡样芽孢杆菌(Bacillus cereus).其中原本对红霉素极度敏感的蜡样芽孢杆菌，由于低浓度红霉素的驯化，对红霉素产生了一定的耐药性.

　　2)长期受低浓度红霉素污染的底泥中，耐红霉素的细菌类群分三大类群，即：未获培养的细菌、芽孢杆菌纲和梭状芽孢杆菌纲.其中未获培养的细菌在整个文库中比例最大，可培养的耐药菌的优势类群为芽孢杆菌纲.结合系统发育树的分析，发现所占比例最大的某未获培养的细菌为某种未分类的动性球菌(43.37%)，可初步划分到芽孢杆菌纲(Bacilli)，尚需更多地分类鉴定.在16S rDNA克隆文库中，有25.3%的细菌与GenBank中已知细菌的序列同源性低于97%，这说明底泥中蕴藏着许多未知的耐药细菌，有待于进一步深入研究.

　　3)通过对16S rDNA文库研究发现，对红霉素具有一定耐药性的细菌既有致病菌(例如蜡样芽孢杆菌)，也有有益菌(例如苏云金杆菌).随着各种抗生素使用量的增加，一方面，使得原本对抗生素有一定敏感度的致病菌表现出了一定的耐药性，对疾病的防治加大了难度;另一方面，也使得环境中某些有益菌获得了一定的耐药性，使其可以更好的生存.在抗生素污染问题越来越严重的今天，有益的芽孢杆菌的应用研究，可能是解决抗生素问题的一个有效方案.