畜牧兽医学报  2023, Vol. 54 Issue (12): 5184-5197. DOI: 10.11843/j.issn.0366-6964.2023.12.028    PDF    
海南省奶源和养殖环境主要肠杆菌科细菌的分布及基因分型
郝若晨1,2, 唐敏嘉2, 刘光亮1,3, 张艳1, Muhammad Shoaib2, 尚若锋2, 曹宗喜1, 蒲万霞2     
1. 海南省农业科学院畜牧兽医研究所 海南省热带动物繁育与疫病研究重点实验室, 海口 571100;
2. 中国农业科学院兰州畜牧与兽药研究所/农业农村部兽用药物创制重点实验室/甘肃省新兽药工程重点实验室, 兰州 730050;
3. 中国农业科学院兰州兽医研究所, 兰州 730046
摘要:本研究从2021年海南省两个奶牛场采集的52份样本中分离肠杆菌科细菌, 并对分离菌株进行了毒力基因、质粒型和生物被膜形成能力检测及肠杆菌科细菌基因间重复序列分子分型(ERIC-PCR)研究, 以期明确牛奶及奶牛场环境中肠杆菌科细菌分布流行状况及其生物学特性。首先, 采用16S rDNA PCR和ERIC-PCR对分离的肠杆菌科细菌进行鉴定和去重, 结果显示52份样本中分离到肠杆菌科细菌49株, 其中大肠埃希菌(Escherichia coli)32株, 肺炎克雷伯菌(Klebsiella pneumoniae)13株, 阴沟肠杆菌(Enterobacter cloacae)3株, 产气肠杆菌(Enterobacter aerogenes)1株; 其次, 采用PCR方法进行了分离菌株的毒力基因检测, 结果显示, ompA的检出率最高; 之后, 采用CARATTOLI创建的基于PCR的质粒复制子分型, 结果显示K质粒的检出率最高; 采用刚果红琼脂法对分离菌株进行了生物被膜形成能力定性检测, 发现49株肠杆菌科细菌中有37株指示有生物被膜形成能力; 最后, 通过ERIC-PCR分型进行了肠杆菌科细菌的亲缘关系和遗传多样性分析, 发现49株肠杆菌科菌株可划分为17型(Ⅰ~ⅩⅦ), Ⅵ型为优势型(均为Escherichia coli), 共20株。整体研究表明肠杆菌科细菌在奶源环境中广泛存在、生物特性多样, 推测具有一定的致病力和耐药性, 引起人类和动物疾病风险较高, 提示养殖、运输中的卫生和规范操作至关重要。
关键词肠杆菌科细菌    分子分型    质粒型    生物被膜    
Distribution and Genotyping of Major Enterobacteriaceae Bacteria in Milk Sources and Dairy Farm Environment of Hainan Province
HAO Ruochen1,2, TANG Minjia2, LIU Guangliang1,3, ZHANG Yan1, MUHAMMAD Shoaib2, SHANG Ruofeng2, CAO Zongxi1, PU Wanxia2     
1. Hainan Provincial Key Laboratory of Tropical Animal Breeding and Disease Research, Institute of Animal Science and Veterinary Medicine, Hainan Academy of Agricultural Sciences, Haikou 571100, China;
2. Gansu Province Key Laboratory of New Animal Drug Project, Ministry of Agriculture and Rural Affairs Key Laboratory of Veterinary Drug Discovery/Lanzhou Institute of Husbandry and Pharmaceutical Sciences of CAAS, Lanzhou 730050, China;
3. Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou 730046, China
Abstract: In this study, major bacteria in Enterobacteriaceae were isolated from 52 samples collected from two dairy farms in Hainan Province in 2021. The virulence genes detection, plasmid type, biofilm-forming ability and enterobacterial repetitive intergenic consensus (ERIC) molecular typing of isolates were studied in order to clarify the distribution, prevalence and biological characteristics of the isolates in Enterobacteriaceae in milk and dairy farm environment. First, the isolates in Enterobacteriaceae were identified by 16S rDNA PCR. The results showed that 49 strains of Enterobacteriaceae were isolated from 52 samples, including 32 strains of Escherichia coli, 13 strains of Klebsiella pneumoniae, 3 strains of Enterobacter cloacae and 1 strain of Enterobacter aerogenes. The results of virulence gene testing showed that ompA had the highest detection rate. Plasid typing is performed by PCR-based replicon typing created by CARATTOLI. The results showed that the detection rate of K plasmid was the highest among 49 isolated bacteria. The biofilm-formating ability of isolates was qualitatively detected using the Congo Red Agar method. Thirty-seven isolates were positive. The genetic relationship and diversity of Enterobacteriaceae were determined by ERIC-PCR. The results showed that 49 strains of Enterobacteriaceae were divided into 17 types (I-XVI), and the dominant type was type VI (all Escherichia coli), with a total of 20 strains. The overall study showed that Enterobacteriaceae widely existed in the milk and environment, causing a high risk of human and animal diseases. It is suggested that hygiene and standardized operation in breeding and transportation are crucial.
Key words: Enterobacteriaceae    molecular typing    plasmid type    biofilm    

近十年来乳制品的生产和消费大幅增加,并存在强劲增长的势头,同时乳制品污染成为全世界非常关注的问题。乳制品含有丰富的营养成分,包括碳水化合物、蛋白质和矿物质,可以促进微生物的生长[1]。200多种已知的病原体可以通过食物传播,其中有几种肠杆菌科细菌属于对人类和动物会造成严重影响的重要食源性病原体,如产生毒素的沙门菌(Salmonella)、志贺菌(Shigella)和大肠埃希菌(Escherichia coliE. coli)[2]。食品中的致病性E. coli污染是一个公共卫生问题,这些细菌会造成严重的肠道疾病[3]

细菌的致病力与其携带的毒力基因种类与数量高度相关。毒力基因种类繁多,包括黏附素、铁摄取系统、脂多糖、多糖荚膜和侵袭素等,这些毒力基因通常位于致病岛(pathogenicity islands,PAIs)、质粒和其它移动遗传元件(mobile genetic elements,MGEs)上[4]。高致病性毒力岛(high-pathogenicity island,HPI)是一种重要的毒力基因,其核心功能区由fyuAirp2组成,因此fyuAirp2是检测HPI的标志性基因[5]ibeB是重要的致病因子[6],与免疫逃避相关的毒力基因主要包括荚膜、补体抗性蛋白和外膜蛋白酶等,ompAompTtraTcvaCiss分别编码外膜蛋白、外膜蛋白酶、补体抗体蛋白、定殖因子以及抗吞噬作用因子[7]。铁是许多细菌生长所必需的,摄铁系统是细菌重要的生存机制,iroN编码铁载体受体,iucD编码NADPH依赖性赖氨酸N(6′)-单加氧酶,该酶与铁载体生物合成蛋白相关[7]

质粒是能够在细胞之间进行自我传递的自主DNA分子,携带的基因对宿主的生长或生存是不必要的[8]。几乎所有的细菌都有质粒,它们可以垂直传播到后代,也可以水平传播。质粒通过转座子或插入序列等MGEs获得新基因,并能够在广泛的宿主中复制,使它们成为肠杆菌科细菌的耐药性(antimicrobial resistance,AMR)和致病性传播的完美载体[9]。六种主要的质粒家族已被证明可以介导肠道细菌物种之间的抗菌素耐药性传播,即Inc F、Inc A/C、Inc L/M、Inc N、Inc I和Inc HI2[10]。由于毒力因子可增加细菌存活率并可能影响抗生素耐药性表达[11],因此,对不同细菌宿主中毒力基因和质粒特征的鉴定是了解细菌致病力和耐药性的基础。

生物被膜是由细菌在生长过程中为适应生存环境而黏附于物体或活性组织表面并包被其自身而产生的细胞外多糖基质形成的,是细菌的特殊存在形式。它的形成一方面促进细菌逃避机体免疫系统的作用,阻碍细胞吞噬或减少吞噬作用与氧活性降低后的应激反应[12];另一方面可以阻止或延缓药物的渗透,生物被膜内细菌的生理学特点影响了细菌对药物的敏感性[13-14]。生物被膜结构可以破坏人类的防御系统,为微生物提供庇护所,导致免疫逃避和细菌耐药性[15]E. coli、肺炎克雷伯菌(Klebsiella pneumoniaeK. pneumoniae)、Salmonella、铜绿假单胞菌、金黄色葡萄球菌、炭疽杆菌等大多数细菌能产生生物被膜[16]。美国国立卫生研究院表示80%的慢性感染与生物被膜形成相关[17],因此检测细菌生物被膜形成能力是了解其致病力和耐药性的另一个主要的生物学特性。

了解特定地理区域内流行菌株的遗传关系对于疾病的预防是至关重要的。肠杆菌科细菌基因间重复一致序列-聚合酶链式反应分型(enterobacterial repetitive intergenic consensus-PCR, ERIC-PCR)技术是一种用于细菌流行病学分析和基因分型的分子方法,基于基因间重复共有序列进行PCR扩增,由条带数目及大小进行分型,可以评估菌株的亲缘关系和遗传多样性。我国要将海南省建设为面向全球的自由贸易试验区,食品安全问题就显得尤为重要,但海南省奶牛养殖业欠发达,目前很少有研究聚焦于牛奶生产环境中被细菌污染相关的风险因素。

因此本研究从海南省两个奶牛场采集鲜乳、工人手臂、挤奶设备、牛舍护栏、粪便和商品奶样本,共52份。对其进行了样本中肠杆菌科细菌的分离与鉴定、肠杆菌科细菌质粒型的检测、肠杆菌科细菌生物被膜表型的检测以及肠杆菌科细菌基因间重复序列分子分型,旨在探明海南省肠杆菌科细菌在生乳、商品奶及奶牛养殖环境的分布及分型情况,分析菌株的遗传关系。

1 材料与方法 1.1 试验材料

PCR和全自动紫外凝胶成像仪(美国Applied Biosystems公司);移液器、高速离心机和紫外分光光度计(德国Eppendorf公司);电泳槽和高压电泳仪(北京六一仪器厂);多功能酶标仪(美国Gene公司);琼脂糖购自西班牙Biowest公司;刚果红培养基(Congo Red Agar,CRA)购自上海中秦化学试剂有限公司;Premix TaqTM、10×Taq Buffer(Mg2+)和Taq酶(5 U·μL-1)均购自宝生物工程(大连)有限公司(大连TaKaRa公司)。

1.2 样品采集

样品于2021年采集自海南省两个奶牛养殖场,其中A牛场养殖规模为200头,采集健康奶牛鲜乳10份、工人手臂拭子4份、挤奶设备拭子4份、粪便5份和商品奶8份;B牛场养殖规模为100头,采集健康奶牛鲜乳9份、工人手臂拭子3份、牛舍护栏拭子7份和粪便2份。

奶样采集时,先将乳房擦拭消毒,丢弃前3把奶,分别采集每头牛4个乳区的奶3~4 mL,混合后为一份样品;采集工人手臂、牛舍护栏以及挤奶设备源样品时,先将灭菌棉签蘸取适量生理盐水,然后用棉签分别在工人手臂、牛舍护栏以及挤奶设备处涂抹,放入有培养基的采样管中;粪样采用五点采样法收集。

1.3 细菌分离鉴定和去重

称取粪样100 g,加入400 mL灭菌BHI肉汤,混匀,吸取5 mL上清液接种于100 mL灭菌BHI液体培养基中,37 ℃ 180 r·min-1增菌15 h。吸取1 mL奶样,接种至9 mL BHI肉汤中进行增菌。无菌条件下蘸取增菌液,在MAC平皿上进行单划线接种,37 ℃静置培养18 h,挑取形态不一致的单菌落接种至MAC上培养,直至菌落形态单一。采用水煮法提取细菌基因组DNA,以纯化菌株的DNA为模板,进行ERIC-PCR扩增,ERIC-PCR反应体系为25 μL:10×Taq Buffer(Mg2+)2.5 μL、dNTP 2 μL、Taq酶1 μL、ERIC引物各0.5 μL(10 μmol·L-1),模板DNA 1.5 μL,ddH2O 17 μL。反应条件:94 ℃预变性5 min,94 ℃变性45 s,52 ℃退火1 min,72 ℃延伸5 min,35个循环;72 ℃最后延伸10 min。产物在2.0%琼脂糖凝胶上进行电泳并拍照,将分离自同一样品中条带数目、大小完全一致的菌株视为重复,仅保留一株。以纯化非重复菌株DNA为模板,扩增其16S rDNA序列,PCR反应体系(25 μL):Premix TaqTM12.5 μL,上下游引物各1 μL(10 μmol·L-1),模板1 μL,加ddH2O至25 μL。反应条件:94 ℃预变性5 min;94 ℃变性30 s,退火30 s,72 ℃下1 kb·min-1延伸,30个循环;72 ℃延伸5 min。扩增产物在1.0%琼脂糖凝胶上进行电泳并拍照记录结果。将符合预期大小的PCR产物送上海生工生物有限公司进行测序,在NCBI中将测序结果进行BLAST比对分析,按相似性≥99%判定,确定分离菌株的种属地位。经分离、纯化与鉴定的非重复性肠杆菌科细菌科菌株,于-70 ℃超低温冰箱保存菌种。

1.4 肠杆菌科细菌的毒力基因检测和质粒分型

根据参考文献合成了11种毒力基因的引物,引物由北京擎科生物有限公司合成,引物序列、扩增片段大小及退火温度见表 1。对49株肠杆菌科细菌菌株进行ibeBiucDiroNfyuAirp2、ompAompTtraTisscvaCtsh毒力基因的检测。与铁转运相关的fyuAirp2基因采用双重PCR检测,反应体系为20 μL:fyuAirp2基因的上下游引物各0.5 μL(10 μmol·L-1),Premix TaqTM 10 μL,模板1 μL,再补充去离子水至20 μL。其余毒力基因的检测均使用单反应PCR方法,反应体系与反应条件与“1.3 ”中16S rDNA扩增体系一致。以实验室保存的携带相应特异性基因的菌株作为阳性对照,灭菌超纯水为空白对照。在1%的琼脂糖凝胶上电泳,然后使用全自动凝胶成像系统拍照并记录结果。

表 1 毒力基因引物序列 Table 1 Primer sequences of virulence genes

采用Carattoli等[27]开发的基于肠杆菌科主要质粒不相容性组的复制子分型(PCR-based replicon typing,PBRT),分析所分离肠杆菌科菌株的质粒携带情况。根据参考文献合成了识别parA-parB、iterons和RNAI等目标位点的HI1、HI2、I1、X、L/M、N、FIA、FIB、W、Y、P、FIC、A/C、T、FII、F、K和B/O等18种常见的复制子引物[27],通过PCR方法对上述菌株携带的质粒型进行检测。该方法包括5组多重PCR和三个单重PCR反应,其引物序列、扩增片段大小及退火温度见表 2。以细菌总DNA为模板,使用可识别Inc复制子的18种引物进行基于PCR的复制子分型。

表 2 质粒分型引物信息 Table 2 Plasmid typing primers information

质粒分型基因检测的反应体系:上下游引物各0.4 μL(10 μmol·L-1),DNA模板1 μL,Premix TaqTM 10 μL,ddH2O 8.2 μL,总体积为20 μL。质粒分型PCR的set 1~5和7~8反应条件:预变性94 ℃ 5 min;94 ℃变性30 s,60 ℃退火30 s,72 ℃延伸1 min,35个循环;最后72 ℃延伸7 min。质粒分型PCR set 6的反应条件:预热变性94 ℃ 5 min;94 ℃变性30 s,52 ℃退火30 s,72 ℃延伸1 min,35个循环;最后72 ℃延伸7 min。

1.5 菌株生物被膜定性检测

复苏菌株,将复苏后的菌株划线接种于CRA上,37 ℃培养22~24 h,观察菌落在培养基上的颜色和菌落形态。如果菌落在CRA上呈干燥、黑色、光亮结晶和周边培养基褪色,指示为生物被膜阳性菌株;菌落在CRA上呈湿润、红色和周边培养基不褪色,指示为生物被膜阴性菌株。

1.6 肠杆菌科分离株的ERIC-PCR分子分型

以所分离菌株的基因组DNA做模板,利用ERIC引物对其进行PCR扩增,扩增产物在2%的琼脂糖凝胶上进行电泳,使用DL5000 DNA Marker,电泳后用紫外凝胶成像系统观察并拍照记录,并通过BioNumerics软件进行遗传进化关系分析。

2 结果 2.1 细菌分离鉴定

52份样本中有31份样本中能获得可培养的菌落,总分离率为59.6%(31/52),共分离77株菌。分别为A牛场分离率64.5%(20/31),分离得到57株菌;B牛场分离率52.4%(11/21),分离得到20株菌。粪样的分离率最高,为100.0%(7/7);其次为鲜乳样,分离率为78.9%(15/19),其中A牛场为90.0%(9/10),B牛场为66.7%(6/9);工人手臂源样本的分离率为57.1%(4/7),其中A牛场为75.0%(3/4),B牛场为33.3%(1/3);挤奶设备或牛舍护栏源样本的分离率为36.4%(4/11),其中A牛场为50.0%(2/4),B牛场为28.6%(2/7);商品奶样本的分离率为12.5%(1/8)(表 3)。整体而言,A牛场的分离率较B牛场高,需要提到的一点是,B牛场无商品奶。

将疑似菌株的16S rRNA序列与NCBI数据库进行BLAST比对,鉴定结果如表 4所示。肠杆菌科为49株,其它菌株24株,NCBI数据库中未鉴定的有4株。在鉴定的所有物种中,E. coli的检出率最高,为41.6%(32/77);其次为K. pneumoniae,检出率为16.9%(13/77)。阴沟肠杆菌(Enterobacter cloacaeE. cloacae)的检出率为3.9%(3/77);产气肠杆菌(Enterobacter aerogenesE. aerogenes)检出率为1.3%(1/77)。

表 3 不同来源样本细菌分离率 Table 3 Bacterial isolation rate of samples from different sources
2.2 肠杆菌科细菌的毒力基因和质粒分型

毒力基因结果显示,ompA的检出率最高,为93.9%。A牛场鲜乳源菌株中ompA的检出率最高,为95.5%(21/22),未检出iroNiss基因。A牛场工人手臂源菌株中,ompA的检出率最高,为100.0%(7/7),未检出iucDfyuAirp2、ompTtraTisscvaC基因。A牛场挤奶设备或牛舍护栏源菌株中只检出了ibeBompA基因。A牛场粪源菌株中,ompA的检出率最高,为100%(15/15)。B牛场中工人手臂源菌株中只检出了ompA基因,粪源菌株中检出了ibeBompAompT基因。

质粒型检测结果显示(表 6),K质粒的阳性率最高,为44.9%(22/49);其次为W(40.8%,20/49)、FIB(36.7%,18/49)、Y(18.4%,9/49)、FrepB(14.3%,7/49)、FIA(6.1%,3/49)和HI2(2.0%,1/49)质粒。本研究中未检测到HI1、I1、L/M、N、X、P、FIC、T、A/C、FIIS和B/O质粒。2.0%的菌株携带6个质粒,4.1%的菌株携带5个质粒,12.2%的菌株携带4个质粒,20.4%的菌株携带3个质粒,6.1%的菌株携带2个质粒,8.2%的菌株携带1个质粒,46.9%的菌株未携带所检测的质粒。携带质粒的菌株中,96.2%(25/26)为E. coli,3.8%(1/26)为E. aerogenes

表 4 不同来源样本细菌分离鉴定情况 Table 4 Isolation and identification of bacteria from different sources
表 5 毒力基因在不同样品来源分离株中的分布 Table 5 Distribution of virulence genes in isolates from different sample sources
表 6 分离菌株携带质粒型情况 Table 6 plasmid type of isolated strains
2.3 菌株生物被膜定性

49株肠杆菌科细菌中有37株(75.5%,37/49)在CRA上有不同程度的黑色菌落出现,指示具有生物被膜形成能力,见图 1a。其中E. coli 20株(62.5%,20/32),K. pneumoniae 13株(100.0%,13/13),E. cloacae 3株(100.0%,3/3),E. aerogenes 1株(100.0%,1/1)。其余12株(24.5%,12/49)在CRA上为红色菌落,没有生物被膜形成能力,见图 1b,均为E. coli

a. 刚果红染色阳性;b..刚果红染色阴性 a. Congo red-positive; b. Congo red-negative 图 1 刚果红琼脂上培养的菌株外观 Fig. 1 The appearance of strains cultured on Congo red agar
2.4 肠杆菌科分离株的ERIC-PCR分子分型

以条件位置差异容许度为1.5%,优化值为1.5%对ERIC-PCR指纹图谱进行分析,依据UPGMA方式进行聚类,以70%为标准划分[28],最终49株肠杆菌科菌株被划分为Ⅰ-ⅩⅦ共17型,Ⅵ型为优势型,共20株,其次为ⅩⅤ型(6株)、Ⅷ(3株)、Ⅸ(3株)、Ⅳ(2株)、Ⅹ(2株)、ⅩⅠ(2株)和ⅩⅢ(2株),剩余的9株菌分布在其它9种型中,聚类关系见图 2。Ⅵ型全部为E. coli,其中19株分离自A牛场(包括鲜乳、工人手臂源和粪源),1株分离自B牛场粪源,说明A牛场菌株的聚类关系较近。ⅩⅤ型全部为K. pneumoniae,均分离自A牛场(包括鲜乳和工人手臂源)。

图 2 肠杆菌科细菌ERIC-PCR基因分型的聚类树图 Fig. 2 Cluster tree of ERIC-PCR genotype of Enterobacteriaceae
3 讨论

一些致病性肠杆菌科细菌不仅给畜牧业发展造成危害,还威胁人类健康,而牛是致病性肠杆菌科细菌引起人类感染(牛肉、乳制品、牛粪便污染)的重要来源[29]。因此,了解牛养殖环境肠杆菌科细菌流行情况对动物和人类预防疾病和维护健康至关重要,海南省的养殖业欠发达,对养殖环境和商品奶中肠杆菌科细菌的研究较少,故本研究以养殖环境和商品奶的肠杆菌科细菌为研究对象。本研究中鲜乳样细菌的分离率较高,为78.9%,其中A牛场更是高达90.0%。分离率较高的原因可能是样本来源牛场中奶牛饲养密度和活动量较大,接触污染源多,未及时清洗或未彻底清洗乳头;除此之外,海南省气候湿润,气温高,也容易滋生细菌。本研究中工人手臂源样本和挤奶设备或牛舍护栏源细菌的分离率分别为57.1%和36.4%,可能是因为挤奶工人或挤奶过程中消毒不当,这会导致牛奶或其衍生物被污染,必须引起我们足够的重视,因为摄入受病原体污染的牛奶或其衍生物可能会导致严重的人类感染,提示应对食品行业的工人提出安全和耐药性的建议和持续教育,以防止食品污染和耐药菌株的出现。对饲养环境勤消毒,勤通风与勤打扫,以及及时正确的清洁手臂和挤奶设备等。商品奶的分离率为12.5%,主要的原因是运输过程中样本意外暴露,但未及时处理导致样本被污染。

自然界中E. coli分布广泛,故本研究中,E. coli检出率最高,为41.6%。该菌可以通过粪便和废水等途径排放到环境中,由食物链进入人体,但大多数E. coli都是动物肠道中的正常寄居菌,只有很小一部分才会在一定条件下引起疾病。K. pneumoniae是一种革兰阴性菌,也是机会性和环境性病原体,广泛存在于动物黏膜、水和土壤等环境中,不仅导致奶牛乳腺炎,还会引起医院医疗相关感染[30]。在本研究中,A牛场的鲜乳、工人手臂源和粪样中均能分离得到K. pneumoniae。随着消费鲜乳的趋势,A牛场的鲜乳和工人手臂上分离到该菌必须引起我们足够的重视,因为携带该菌的鲜乳被人类食用后可能会引起严重的感染,而手臂上存在该菌也有可能造成携带者的感染或者进一步的传播。

为进一步确定分离株的致病潜力,又对其毒力基因进行了检测。有研究显示,只要同时含有papCiucDirp2、tshvataatAisscva 8个基因中4个以上即为强毒株[28],本研究中来自A牛场的粪源E. coli菌株HNAF202140同时携带了所检测的11个毒力基因(ibeBiucDiroNfyuAirp2、ompAompTtraTisscvaCtsh),说明该菌株存在着较大的致病风险。A和B牛场毒力基因检测率最高的是ompA。研究报道ompA是不同毒力基因中的流行率最高,最常见的基因,它在生物膜形成、维持膜完整性和抗生素转运中发挥着重要作用[31],这会加强细菌的抗生素耐药性,不仅如此,该基因中的L3位点可能与K. pneumoniae的致病性有关[32],所以ompA在本研究中的高检出率需要引起我们高度重视。目前,已有一些关于基因敲除技术的研究表明,iucDE. coli的致病力存在着某种联系[33]。贾青辉等[34]报道鸡源致病性E. coliiucD的检出率为95.2%,而在本研究中iucD的检出率为4.1%,明显低于上述报道,说明海南省奶源和养殖环境中致病性菌株较少,向环境、人或其它动物传播致病风险的可能性较低。徐政平[35]报道的APEC中的iss的检出率为69.8%;在Manita等[36]的研究中,APEC分离株中ompTiroN的检出率为100%;De Carli等[37]报道的APEC菌株中ompTiroN的检出率分别为100%和98.8%。本研究中issompTiroN的检出率分别为2.0%、8.2%和4.1%,远远低于上述报道。Asai等[38]发现这5种毒力基因(iutAhlyFissiroNompT)与APEC的高致病力相关,而本研究中的E. coli并未确定为APEC,这可能是issompTiroN的检出率较低的原因。总体来看本研究中毒力基因总体检出率相较于其它报道低,可能是因为本研究样品来源于养殖环境和鲜乳样。

病原菌对抗菌素的耐药情况可以从调查菌株中质粒携带情况入手。PBRT具有较高的特异性和灵敏性,但该方法也存在着一定的局限性,只能检测18种常见的质粒,对于一些新型的或变异的质粒可能无法识别。本研究中质粒类型多样,其中Inc K、Inc W和Inc FIB为优势质粒,这些优势质粒也经常在动物源E. coli中发现[39-40]。有研究显示携带AMR基因的质粒主要在典型的Inc组质粒上[30]。其中Inc K质粒主要与blaCMY-2blaCTX-M-14抗性基因在欧洲传播有关,尤其是西班牙和英国。而CMY型和CTX-M型β-内酰胺酶是E. coli中最常见的超广谱β-内酰胺酶(the extended-spectrum β-lactamases,ESBLs)家族,但是产生CMY和CTX-M的肠杆菌科的直接人畜共患食源性传播并不常见[41]。表明本研究中奶牛场的E. coli危害人类健康的可能不大。Inc W质粒被认为是最小的接合质粒,1980年代,在许多细菌中发现了Inc W质粒[42]。研究显示,Inc W质粒携带氯霉素、四环素、磺胺类、庆大霉素和甲氧苄啶等耐药基因(antimicrobial resistance genes,ARGs),也有研究显示携带blaKPC-2blaVIM-1抗性基因[43-44],这两个抗性基因都与碳青霉烯类药物耐药性有关[43]。Inc F组质粒是大小在45~200 kb的低拷贝接合质粒,Inc F质粒上最常见的耐药基因有编码ESBLs、碳青霉烯酶、氨基糖苷类修饰酶的基因[45]。这些质粒不仅会对细菌的生存做出贡献,其介导的耐药机制也是抗生素耐药性增加的原因。

CRA是生物被膜定性常用的方法之一。本研究生物被膜定性检测结果显示49株肠杆菌科细菌中有37株指示具有生物被膜形成的能力,在E. coliK. pneumoniaeE. cloacaeE. aerogenes中都能发现生物被膜形成的能力,而且E. coli中也有无生物被膜形成能力的菌株,表示生物被膜形成能力似乎不是菌种依赖性的。Staji等[46]发现,ERIC-PCR检测结果中同一来源分离菌株基因存在差异,不同来源也可有较近的亲缘关系。郑晓风等[47]用ERIC-PCR方法发现菌株间存在交叉传播。本研究ERIC-PCR检测结果显示49株肠杆菌科细菌被分为17型,参照我国华东地区的103株临床乳房炎E.coli被分为10型[48],可见本研究中肠杆菌科细菌广泛的DNA多样性。DNA多样性的原因可能是频繁使用抗生素导致的选择性压力引起的基因修饰。Ⅵ型菌株来源于鲜乳、粪便和工人手臂源,XⅤ型来源于鲜乳和工人手臂,说明粪便和工人手臂是鲜乳污染或者乳腺感染的主要来源,会进一步导致乳腺炎的发生。

4 结论

本研究样品的总分离率为59.6%(31/52),共分离出77株菌,不同来源的样品均有分布,粪便检出率最高,为100.0%;其次为鲜乳,检出率78.9%。经鉴定疑似菌株覆盖了5个科11个种,主要为肠杆菌科。所分离肠杆菌科细菌毒力基因的检测中ompA的检出率最高,为93.9%,其余毒力基因均有不同程度的检出;质粒类型多样,携带3个质粒型的菌株较多,携带质粒的菌株中,96.2%为E. coli,其中75.5%指示有生物被膜形成能力;分子分型中Ⅵ型为优势型,全部为E. coli;其次是ⅩⅤ型,全部为K. pneumoniae。肠杆菌科细菌在奶源环境中广泛存在、生物特性多样,推测具有一定的致病力和耐药性。

参考文献
[1]
REMENANT B, JAFFRÈS E, DOUSSET X, et al. Bacterial spoilers of food: behavior, fitness and functional properties[J]. Food Microbiol, 2015, 45: 45-53. DOI:10.1016/j.fm.2014.03.009
[2]
DEVLEESSCHAUWER B, PIRES S M, YOUNG I, et al. Associating sporadic, foodborne illness caused by Shiga toxin-producing Escherichia coli with specific foods: a systematic review and meta-analysis of case-control studies[J]. Epidemiol Infect, 2019, 147: e235. DOI:10.1017/S0950268819001183
[3]
CASTRO V S, CARVALHO R C T, CONTE-JUNIOR C A, et al. Shiga-toxin producing Escherichia coli: pathogenicity, supershedding, diagnostic methods, occurrence, and foodborne outbreaks[J]. Compr Rev Food Sci Food Saf, 2017, 16(6): 1269-1280. DOI:10.1111/1541-4337.12302
[4]
TIVENDALE K A, LOGUE C M, KARIYAWASAM S, et al. Avian-pathogenic Escherichia coli strains are similar to neonatal meningitis E. coli strains and are able to cause meningitis in the rat model of human disease[J]. Infect Immun, 2010, 78(8): 3412-3419.
[5]
SCHUBERT S, PICARD B, GOURIOU S, et al. Yersinia high-pathogenicity island contributes to virulence in Escherichia coli causing extraintestinal infections[J]. Infect Immun, 2002, 70(9): 5335-5337. DOI:10.1128/IAI.70.9.5335-5337.2002
[6]
HUANG S H, WAN Z S, CHEN Y H, et al. Further characterization of Escherichia coli brain microvascular endothelial cell invasion gene ibeA by deletion, complementation, and protein expression[J]. J Infect Dis, 2001, 183(7): 1071-1078. DOI:10.1086/319290
[7]
SAROWSKA J, FUTOMA-KOLOCH B, JAMA-KMIECIK A, et al. Virulence factors, prevalence and potential transmission of extraintestinal pathogenic Escherichia coli isolated from different sources: recent reports[J]. Gut Pathog, 2019, 11: 10. DOI:10.1186/s13099-019-0290-0
[8]
THOMAS C M, NIELSEN K M. Mechanisms of, and barriers to, horizontal gene transfer between bacteria[J]. Nat Rev Microbiol, 2005, 3(9): 711-721. DOI:10.1038/nrmicro1234
[9]
ARREDONDO-ALONSO S, ROGERS M R C, BRAAT J C, et al. Mlplasmids: a user-friendly tool to predict plasmid- and chromosome-derived sequences for single species[J]. Microb Genom, 2018, 4(11): e000224.
[10]
SCHLÜTER A, NORDMANN P, BONNIN R A, et al. IncH-type plasmid harboring blaCTX-M-15, blaDHA-1, and qnrB4 genes recovered from animal isolates[J]. Antimicrob Agents Chemother, 2014, 58(7): 3768-3773. DOI:10.1128/AAC.02695-14
[11]
NÚÑEZ C, PALAVECINO A, GONZÁLEZ I A, et al. Effective photodynamic therapy with Ir(Ⅲ) for virulent clinical isolates of extended-spectrum beta-lactamase Klebsiella pneumoniae[J]. Pharmaceutics, 2021, 13(5): 603. DOI:10.3390/pharmaceutics13050603
[12]
LEID J G, SHIRTLIFF M E, COSTERTON J W, et al. Human leukocytes adhere to, penetrate, and respond to Staphylococcus aureus biofilms[J]. Infect Immun, 2002, 70(11): 6339-6345. DOI:10.1128/IAI.70.11.6339-6345.2002
[13]
侯博, 王晨燕, 周伦江. 毒素-抗毒素系统在细菌生物被膜形成中的作用及调控机制[J]. 畜牧兽医学报, 2022, 53(10): 3326-3334.
HOU B, WANG C Y, ZHOU L J. The roles and regulatory mechanisms of toxin-antitoxin system in bacterial biofilm formation[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(10): 3326-3334. (in Chinese)
[14]
JESAITIS A J, FRANKLIN M J, BERGLUND D, et al. Compromised host defense on Pseudomonas aeruginosa biofilms: characterization of neutrophil and biofilm interactions[J]. J Immunol, 2003, 171(8): 4329-4339. DOI:10.4049/jimmunol.171.8.4329
[15]
MIRZAEI R, MOHAMMADZADEH R, ALIKHANI M Y, et al. The biofilm-associated bacterial infections unrelated to indwelling devices[J]. IUBMB Life, 2020, 72(7): 1271-1285. DOI:10.1002/iub.2266
[16]
喻华英, 李启峰. 大肠杆菌生物膜的筛选及生长曲线测定[J]. 中国农学通报, 2014, 30(5): 23-27.
YU H Y, LI Q F. Screening biofilm and determination growth biofilm of curve in Escherichia coli[J]. Chinese Agricultural Science Bulletin, 2014, 30(5): 23-27. (in Chinese)
[17]
DAVIES D. Understanding biofilm resistance to antibacterial agents[J]. Nat Rev Drug Discov, 2003, 2(2): 114-122. DOI:10.1038/nrd1008
[18]
张宇曦, 韩先干, 左佳坤, 等. 禽致病性大肠杆菌脂多糖核心型分布与毒力基因的相关性分析[J]. 微生物学通报, 2015, 42(8): 1619-1625.
ZHANG Y X, HAN X G, ZUO J K, et al. Distribution of lipopolysaccharide core type in avian pathogenic Escherichia coli and its correlation with virulence gene[J]. Microbiology China, 2015, 42(8): 1619-1625. (in Chinese)
[19]
王利勤. 鸡源致病性大肠埃希菌耐药基因及毒力基因检测研究[D]. 杨陵: 西北农林科技大学, 2013.
WANG L Q. Study on resistance genes and virulence genes of pathogenic Escherichia coli isolates from chickens[D]. Yangling: Northwest A&F University, 2013. (in Chinese)
[20]
孙畅, 赵丽丽, 徐丽晶, 等. 鸭源致病性大肠杆菌的分离鉴定及系统发育分析[J]. 中国兽医科学, 2017, 47(3): 286-294.
SUN C, ZHAO L L, XU L J, et al. Isolation, identification and phylogenetic analysis of avian pathogenic Escherichia coli from duck[J]. Chinese Veterinary Science, 2017, 47(3): 286-294. (in Chinese)
[21]
芮萍, 马增军, 杨彩然, 等. 狐貉貂源致病性大肠杆菌的分离与鉴定[J]. 中国兽医科学, 2015, 45(5): 517-521.
RUI P, MA Z J, YANG C R, et al. Isolation and identification of pathogenic Escherichia coli from fox, racoon dog and mink[J]. Chinese Veterinary Science, 2015, 45(5): 517-521. (in Chinese)
[22]
谭珊, 宋欣媛, 张明通, 等. 云南省昆明市动物园大肠杆菌分离株毒力基因分布特征[J]. 中国动物检疫, 2019, 36(4): 86-91.
TAN S, SONG X Y, ZHANG M T, et al. Distribution of virulence genes of Escherichia coli isolated from the zoo in Kunming City of Yunnan province[J]. China Animal Health Inspection, 2019, 36(4): 86-91. DOI:10.3969/j.issn.1005-944X.2019.04.020 (in Chinese)
[23]
侯晓. 某牦牛场感染性肺炎的病原学研究[D]. 兰州: 甘肃农业大学, 2020.
HOU X. Etiological study on pneumonia infectious in Yak farm[D]. Lanzhou: Gansu Agricultural University, 2020. (in Chinese)
[24]
周磊, 李泽伟, 孙起荣, 等. 54株猪源肠外致病性大肠杆菌血清型、系统进化群和基因型[J]. 微生物学通报, 2021, 48(4): 1182-1194.
ZHOU L, LI Z W, SUN Q R, et al. Serotypes, phylogenetic groups and genotypes of 54 extraintestinal pathogenic Escherichia coli from pig[J]. Microbiology China, 2021, 48(4): 1182-1194. (in Chinese)
[25]
范克伟, 余佳, 杨守深, 等. 华南虎源禽致病性大肠杆菌的分离鉴定及其生物学特性研究[J]. 西北农林科技大学学报: 自然科学版, 2018, 46(12): 1-10.
FAN K W, YU J, YANG S S, et al. Isolation, identification and biological characteristics of avian pathogenic Escherichia coli from South China tiger[J]. Journal of Northwest A&F University: Natural Science Edition, 2018, 46(12): 1-10. (in Chinese)
[26]
白灏, 冀辉, 韩先干, 等. 禽致病性大肠杆菌江苏、安徽分离株的生物学特性分析[J]. 微生物学通报, 2013, 40(7): 1315-1322.
BAI H, JI H, HAN X G, et al. Characterization of avian pathogenic Escherichia coli isolated from Jiangsu and Anhui provinces[J]. Microbiology China, 2013, 40(7): 1315-1322. (in Chinese)
[27]
CARATTOLI A, BERTINI A, VILLA L, et al. Identification of plasmids by PCR-based replicon typing[J]. J Microbiol Methods, 2005, 63(3): 219-228. DOI:10.1016/j.mimet.2005.03.018
[28]
IGWARAN A, OKOH A I. Campylobacteriosis agents in meat carcasses collected from two district municipalities in the Eastern Cape province, South Africa[J]. Foods, 2020, 9(2): 203. DOI:10.3390/foods9020203
[29]
FRIESEMA I H, VAN DE KASSTEELE J, DE JAGER C M, et al. Geographical association between livestock density and human Shiga toxin-producing Escherichia coli O157 infections[J]. Epidemiol Infect, 2011, 139(7): 1081-1087. DOI:10.1017/S0950268810002050
[30]
PODSCHUN R, ULLMANN U. Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors[J]. Clin Microbiol Rev, 1998, 11(4): 589-603. DOI:10.1128/CMR.11.4.589
[31]
YANG J, YUN S, PARK W. Blue light sensing BlsA-mediated modulation of meropenem resistance and biofilm formation in Acinetobacter baumannii[J]. mSystems, 2023, 8(1): e0089722. DOI:10.1128/msystems.00897-22
[32]
SAUREL O, IORDANOV I, NARS G, et al. Local and global dynamics in Klebsiella pneumoniae outer membrane protein a in lipid bilayers probed at atomic resolution[J]. J Am Chem Soc, 2017, 139(4): 1590-1597. DOI:10.1021/jacs.6b11565
[33]
陈艺丽. 禽病原性大肠杆菌E058株iucD、ftsk基因突变株的构建及其致病性评价[D]. 扬州: 扬州大学, 2009.
CHEN Y L. Construction of an iucD mutant and a ftsk mutant of avian pathogenic Escherichia coli and evaluation of their pathogenicity[D]. Yangzhou: Yangzhou University, 2009. (in Chinese)
[34]
贾青辉, 张召兴, 李蕴玉, 等. 河北地区鸡源致病性大肠杆菌毒力基因检测与耐药性分析[J]. 中国预防兽医学报, 2018, 40(5): 401-405.
JIA Q H, ZHANG Z X, LI Y Y, et al. Analysis of virulence genes and drug resistance of the pathogenic Escherichia coli isolates from chicken in Hebei province[J]. Chinese Journal of Preventive Veterinary Medicine, 2018, 40(5): 401-405. (in Chinese)
[35]
徐政平. 禽致病性大肠杆菌流行病学调查及耐药性分析[D]. 南京: 南京农业大学, 2017.
XU Z P. Epidemiological investigation and drug resistance of avian pathogenic Escherichia[D]. Nanjing: Nanjing Agricultural University, 2017. (in Chinese)
[36]
MANITA S, KUMAR B R, BHUMINAND D, et al. Antibiotic resistance pattern and virulence genes content in avian pathogenic Escherichia coli (APEC) from broiler chickens in Chitwan, Nepal[J]. BMC Vet Res, 2018, 14(1): 113. DOI:10.1186/s12917-018-1442-z
[37]
DE CARLI S, IKUTA N, LEHMANN F K M, et al. Virulence gene content in Escherichia coli isolates from poultry flocks with clinical signs of colibacillosis in Brazil[J]. Poult Sci, 2015, 94(11): 2635-2640. DOI:10.3382/ps/pev256
[38]
ASAI T, MASANI K, SATO C, et al. Phylogenetic groups and cephalosporin resistance genes of Escherichia coli from diseased food-producing animals in Japan[J]. Acta Vet Scand, 2011, 53(1): 52. DOI:10.1186/1751-0147-53-52
[39]
COTTELL J L, WEBBER M A, COLDHAM N G, et al. Complete sequence and molecular epidemiology of IncK epidemic plasmid encoding blaCTX-M-14[J]. Emerg Infect Dis, 2011, 17(4): 645-652. DOI:10.3201/eid1704.101009
[40]
MNIF B, KTARI S, RHIMI F M, et al. Extensive dissemination of CTX-M-1- and CMY-2-producing Escherichia coli in poultry farms in Tunisia[J]. Lett Appl Microbiol, 2012, 55(6): 407-413. DOI:10.1111/j.1472-765X.2012.03309.x
[41]
BALLASH G A, ALBERS A L, MOLLENKOPF D F, et al. Antimicrobial resistant bacteria recovered from retail ground meat products in the US include a Raoultella ornithinolytica co-harboring blaKPC-2 and blaNDM-5[J]. Sci Rep, 2021, 11(1): 14041. DOI:10.1038/s41598-021-93362-x
[42]
FERNÁNDEZ-LóPEZ R, GARCILLÁN-BARCIA M P, REVILLA C, et al. Dynamics of the IncW genetic backbone imply general trends in conjugative plasmid evolution[J]. FEMS Microbiol Rev, 2006, 30(6): 942-966. DOI:10.1111/j.1574-6976.2006.00042.x
[43]
ALMEIDA A C, DE SÁ CAVALCANTI F L, VILELA M A, et al. Escherichia coli ST502 and Klebsiella pneumoniae ST11 sharing an IncW plasmid harbouring the blaKPC-2 gene in an Intensive Care Unit patient[J]. Int J Antimicrob Agents, 2012, 40(4): 374-376. DOI:10.1016/j.ijantimicag.2012.05.022
[44]
MIRIAGOU V, DOUZINAS E E, PAPAGIANNITSIS C C, et al. Emergence of Serratia liquefaciens and Klebsiella oxytoca with metallo-β-lactamase-encoding IncW plasmids: further spread of the blaVIM-1-carrying integron In-e541[J]. Int J Antimicrob Agents, 2008, 32(6): 540-541. DOI:10.1016/j.ijantimicag.2008.06.022
[45]
ROZWANDOWICZ M, BROUWER M S M, FISCHER J, et al. Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae[J]. J Antimicrob Chemother, 2018, 73(5): 1121-1137. DOI:10.1093/jac/dkx488
[46]
STAJI H, BIRGANI S F, RAEISIAN B. Comparative clustering and genotyping of Campylobacter jejuni strains isolated from broiler and turkey feces by using RAPD-PCR and ERIC-PCR analysis[J]. Ann Microbiol, 2018, 68(11): 755-762. DOI:10.1007/s13213-018-1380-9
[47]
郑晓风, 张妍, 刘英玉, 等. 新疆部分地区牛羊源产志贺毒素大肠埃希菌菌株检测与分析[J]. 畜牧兽医学报, 2020, 51(10): 2518-2527.
ZHENG X F, ZHANG Y, LIU Y Y, et al. Detection and analysis of shiga toxin-producing Escherichia coli isolates from cattle and sheep sources in some regions of Xinjiang, China[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(10): 2518-2527. DOI:10.11843/j.issn.0366-6964.2020.10.020 (in Chinese)
[48]
JAVED M, JAM K, FAN H J. Serotypes, genotypes, virulence factors and antimicrobial resistance genes of Escherichia coli isolated in bovine clinical mastitis from Eastern China[J]. Pak Vet J, 2016, 36(4): 493-498.

(编辑   白永平)