四川部分地区禽源空肠弯曲菌MLST分型及遗传进化分析
姚学萍1,2, 张博1, 王寒1, 邬旭龙1, 周丽军1, 曹随忠1,2, 王印1,2, 杨泽晓1,2, 德西措姆3
1.四川农业大学动物医学院,成都 611130
2.动物疫病与人类健康四川省重点实验室,成都 611130
3.西藏昌都市左贡县畜牧站,左贡 854400

第一作者:姚学萍,Email:yaoxueping74@126.com; ORCID:0000-0002-6800-9912

摘要
目的 为了解四川部分地区分离到的48株禽源空肠弯曲菌间的分子特征,以对来源不同的空肠弯曲菌进行分子分型及遗传进化分析。方法 本研究参照MLST数据库,选取了空肠弯曲菌MLST分型7个管家基因 aspA, glnA, glyA, gltA, tkt, pgm uncA作为目的基因分别扩增及测序,并将所得序列经比对后进行遗传进化分析。结果 48株分离株中共含24种不同的序列型分属于9种不同的克隆系及11种独特型,其中22株分离株含有新序列型,占分离株的45.83%(22/48);9种不同的克隆系中以CC-21、CC-353、CC-464克隆系所含菌株相对较多,而克隆系CC-45、CC-48最少,均为1株分离株;11种独特型共存于21株分离株中,其中以ST-8675序列型最多,为8株。遗传进化树显示,属同一克隆系的不同型别间的分离株亲缘关系较近,如同属于CC-21克隆系的ST-298与ST-760型别分离株处于同一小分支、亲缘关系较近,但不同克隆系的分离株间遗传进化关系较远,鹌鹑源分离株(V13-20)单独处于一大分支与其余禽源分离株间的亲缘关系最远。总体而言,不同地区、宿主来源的分离株呈现遗传多样性。结论 本研究中分离株间具有基因多样性,MLST分型为了解该地区空肠弯曲菌的遗传特性和分布特点提供了参考依据。
关键词: 空肠弯曲菌; MLST分型; 管家基因; 遗传进化
中图分类号:S855.1 文献标志码:A 文章编号:1002-2694(2019)01-0005-06
MLST typing and genetic evolution of Campylobacter jejuni isolated from poultry in some areas of Sichuan,China
YAO Xue-ping1,2, ZHANG Bo1, WANG Han1, WU Xu-long1, ZHOU Li-jun1, CAO Sui-zhong1,2, WANG Yin1,2, YANG Ze-xiao1,2, DEXI Cuo-mu3
1.College of Veterinary Medicine, Sichuan Agricultural University, Chengdu 611130, China
2.Key Laboratory of Animal Disease and Human Health of Sichuan Province, Chengdu 611130, China
3. Animal Husbandry Station of Tibet Chang Zuogong County, Zuogong 854400, China

Anbout author:Yao Xue-ping and Zhang Bo contributed equally.

Abstract

To determine the molecular characteristics of 48 strains of Campylobacter jejuni isolated from some parts of Sichuan, the molecular typing and genetic evolution of Campylobacter jejuni were analyzed. According to MLST database the selection of Campylobacter jejuni MLST genotyping of 7 housekeeping genes aspA, glnA, glyA, gltA, tkt, pgm and uncA, the sequencing as target genes were amplified, and the sequences were compared with phylogenetic analysis. The result showed that there were 24 different sequences in 48 isolates, belonging to 9 different clones and 11 unique types, among which 22 of the isolates contained a new sequence, accounting for 45.83% (22/48). Among the 9 different clones, CC-21, CC-353, and CC-464 clones contained more strains relatively, while the clones CC-45 and CC-48 contained one isolated. The 11 unique types coexist in 21 isolates, of which 8 strains were the ST-8675 sequence types. Phylogenetic tree showed that the isolates belonging to different clones of the same clonal line have close genetic relationships, like the ST-298 and ST-760 isolates belonging to the CC-21 clonal line were on the same small branch and closely related. However, the genetic relationship between the isolates of different clonal lines was far away. The quail source isolates (V13-20) were in a large branch and has the most distant relationship with other poultry isolates. Overall, isolates from different regions and host sources showed genetic diversity. Gene polymorphism exists in the isolates in this study, and MLST typing provides a reference for understanding the genetic characteristics and distribution characteristics of Campylobacter jejuni in the area.

Keyword: Campylobacter jejuni; MLST typing; housekeeping genes; genetic evolution

空肠弯曲菌(Campylobacter jejuni)是广泛存在于家禽、家畜等动物体内的一种常见的人畜共患的食源性感染致病菌, 主要引起人的急性胃肠炎, 进入血液会引发败血症、脑炎、心内膜炎、心肌炎、反应性关节炎和泌尿系统感染等, 严重的还会引发神经系统损伤性疾病格林巴利综合征(Guillain-Barre syndrome, GBS)[1, 2]。禽类是空肠弯曲菌最主要的传染源, 消费禽肉及制品与该菌的流行紧密相关, 故而该菌早已成为全球公共卫生关注的重点[3], 因此对于空肠弯曲菌的分类研究, 以探究其致病相关性显得尤为重要。

多位点序列分型(Multilocus sequence typing, MLST)技术是一种基于核酸序列测定的新型细菌分型方法, 其与多位点酶电泳分析(Multilocus enzyme electrophoresis, MEE)原理相同, 不同的是MLST通过测定细菌7个管家基因的核苷酸序列后根据序列的不同, 经比对后分别命名不同序列号管家基因的等位基因值, 7个管家基因不同序列号的组合决定菌株特异的序列型(Sequence type, ST)[4, 5]。MLST数据库针对每个菌株个体都有相应的、唯一的ST型别且每个菌属的管家基因及引物均被标准化与统一, 以便全球菌株MLST结果的比较分析。

依此, 本研究以四川部分地区分离到的48株禽源空肠弯曲菌作为研究载体, 使用MLST分子分型技术对其进行分子分型研究及遗传进化分析, 以期进一步从分子水平研究该菌的遗传特性及分布特点。

1 材料与方法
1.1 材料

48株分离株源于对四川多地区采集的禽类盲肠及内容物样品的分离[6], 哥伦比亚培养基、Skirrow培养基及其对应抗生素添加剂、脑心浸液均购自于成都康迪生物科技有限公司, 细菌基因组提取试剂盒、胶回收试剂盒均购自于天根生化科技(北京)有限公司, 其他试剂均为国产分析纯产品。

1.2 分离株MLST分型方法、反应体系及参数

参照MLST数据库(http://www.mlst.net/ databases/)推荐的方法选取空肠弯曲菌 MLST分型相关的7个管家基因(aspA, glnA, glyA, gltA, tkt, pgmuncA)作为目的基因, 分别进行扩增与测序。

反应体系及参数:反应体系为25 μ L:其中2× Taq PCR Mix 12.5 μ L, 引物(10 μ mol/L)各1 μ L, 模板2 μ L, ddH2O补足至25 μ L; 反应参数为:94 ℃ 5 min; 94 ℃ 30 s、50 ℃ 1 min、72 ℃ 1 min, 35个循环, 72 ℃ 10 min。反应结束后, 将产物经琼脂糖凝胶电泳后, 使用试剂盒回收条带单一且明亮的扩增片段送至华大基因测序。

1.3 序列类型(Sequence type, ST)的确定

对7个管家基因的双向测序图谱进行比对, 确保序列一致后将不同的基因序列在MLST数据库中提交, 进行检索比对。根据数据库提供的每种等位基因的序列号, 确定7个管家基因的组合, 即为该菌株的ST型, 并根据该菌株的ST型在数据库中的排列序号对其命名。经命名的ST型或出现新ST型, 其归因于两种情况, 存在新的等位基因序列或存在各等位基因位点的新组合[7]。相似的ST型菌株可组成一个ST序列群, 即一个克隆系(Clonal complex, CC)[8]

1.4 遗传进化分析

确定所有分离株的ST型别及克隆系后, 从MLST数据库网站(http://www.mlst.net/databases/)下载分离株及近5年国内外多个国家不同源分离株的7个等位基因的串联序列(按照aspA-glnA-gltA-glyA-pgm-tkt-uncA的顺序排列), 采用MEGA6.0软件中的Neighbor-jioning方法构建系统进化树, 分析本次分离株和其他地区分离株之间的遗传进化关系。

2 结 果
2.1 分离株管家基因扩增结果

使用空肠弯曲菌7个管家基因扩增引物对分离株的管家基因进行扩增, 产物经琼脂糖凝胶电泳后, 结果显示(图1, 其中1株分离株)各管家基因的片段大小与预期结果一致。

图1 空肠弯曲菌的MLST扩增结果Fig.1 MLST amplification results of Campylobacter jejuni
Note: M: DNA Marker; 1:glyA; 2:glnA; 3:gltA; 4:tkt; 5:uncA; 6:aspA; 7:pgm

2.2 MLST分型结果与分析

2.2.1 新序列型及新等位基因值 48株分离株中共含24种不同的序列型, 其中22株分离株含有新序列型属11个新序列型, 占分离株的45.83%。新序列型中眉山鹌鹑源的8株分离株和雅安鸡源的3株分离株全为新序列型; 温江鸡源的分离株中有7株属于新序列型, 占温江鸡源分离株的26.92%; 温江鸭源的分离株中有1株属于新序列型; 金堂鸡源的分离株有3株属于新序列型, 占金堂分离株的33.33%。

48株分离株中共有14株分离株含有新等位基因序列, 7个管家基因中的aspA, pgm, tkt, uncA基因存在新的等位基因序列, 而其余3个等位基因均未发现新的等位基因序列。14株分离株各存在1个新等位基因序列, 其中tkt位点共有10株分离株(V5、V6、V13-20)发现新等位基因序列, 分离株V26、V27的pgm位点为新等位基因序列, aspAuncA位点各有1株分离株为新等位基因序列; 分离株中的V26、V27所属的序列型有对应的克隆系均为ST-1332CC。除此之外, 另有8株分离株(C2、C4、V7、V8、V25、V33、V34、V47)存在等位基因位点的新组合, 属于新序列型且均无对应的克隆系, 为独特型(表1)。

表1 空肠弯曲菌分离株中的新等位基因及新序列型 Tab.1 New allele sequences and new series of alleles of Campylobacter jejuni isolates

2.2.2 48株分离株所属不同序列型和克隆系的数量 经序列比对后, 48株分离株共含有24种不同的序列型, 分属于9种克隆系和11种独特型。其中11种独特型共包含21株分离株, 占总分离株的43.75%。CC-21、CC-353和CC-1332均包含2种以上序列型, 分别占菌株总数的12.5%、8.33%、6.25%, 其余克隆系只包含一种序列型。所有的不同序列型中, 独特型ST-8675包含的菌株数最多, 共包含8株分离株, 占总菌株数的16.67%, ST-760、ST-464、ST-161和ST-354均包含3株以上的分离株, 其余序列型只包含1-2株分离株, 具体信息(表2)。

表2 48株分离株所属不同序列型和克隆系的数量 Tab.2 Different STs and clonal complex contain 48 isolates
2.3 48分离株MLST结果统计与分析

将分离株等位基因扩增结果统计后发现48株分离株中, 等位基因值的数量最多的是pgm, 为14个; 其次glnA, 为12个; gltAtktuncAaspA均为9个; 最少的glyA, 为8个。等位基因值数量的多少在某种程度上反应了该位点基因的多样性, 依此或可推知pgm基因最不稳定, 突变最多, 而glyA基因较稳定。依据本研究中出现的所有等位基因数值及其各数值所对应菌株的数量, 结果显示, glyA基因值为2的细菌数量最多, 为21株; uncA基因值为6的细菌数量次之, 为17株; 而出现的各管家基因值为1的菌株分布最为广泛。

2.4 不同克隆系、序列型包含的不同地区来源菌株的数量

将不同克隆系、序列型包含的不同地区来源的菌株的数量进行统计后, 结果显示, 源于(金堂养鸡场)的9株分离株率属于6种序列型; 源于(温江活禽市场2、3)地区的分离株均有8株率属于6种序列型; 源于(眉山鹌鹑养殖场)的8株分离株均归属于ST-8675序列型, 而源于(温江养鸡场2)的分离株仅有2株率属于2种序列型。针对克隆系而言, 以CC-21包含的分离株最多为6株, 以CC-45包含的分离株最少为1株。除此之外可发现同一个序列型包含2株及以上菌株的, 菌株多为同一地区分离株, 如ST-298包含2株(温江活禽市场1)的分离株, ST-760包含3株(温江养鸡场1)的分离株, 表明菌株的流行与地域特征有关。

2.5 分离株的遗传进化分析

使用MEGA6.0软件对48株分离株所含的24种序列型及国内2株腹泻病人空肠弯曲菌分离株、2株鸡源分离株、1株鸭源分离株, 以及英国、丹麦、波兰、新西兰、阿根廷、韩国、泰国、肯尼亚、印度鸡源分离株各1株, 共38株空肠弯曲菌的序列型进行遗传进化树的构建(图2)。结果显示, 泰国和丹麦的ST-8659和ST-8661型别分离株与其余36个型别分离株处于不同分支, 遗传距离较远; 眉山鹌鹑源的ST-8675型的8株分离株(V13-20)与其同属于一大分支的其余35个型别分离株遗传距离较远; 属于同一克隆系的2种型别菌株遗传距离近, 如同属于CC-21克隆系的ST-298与ST-760处于同一小分支、属于CC-1332克隆系的ST-8671与ST-696处于同一小分支; 来源于同一地区不同场的分离株遗传距离较近, 如来源于(温江活禽市场1)ST-5864序列型的V26分离株和(温江养鸡场2、活禽市场3)ST-8090序列型的V40、V43分离株处于同一小分支, 亲缘关系较近。除此之外, 本研究中序列型为ST-354的V24、V42、V46分离株与英国ST-354序列型分离株处于同一小分支并与印度ST-8657序列型菌株处于同一分支, 暗示着不同克隆系间的序列也存在密切相关性; 国内腹泻病人ST-8727序列型分离株与本研究中鸡源分离株C1亲缘关系最近, 或可推之鸡与人类感染有重要关系。

图2 本研究24种序列型和国内外14株分离株序列型的遗传进化树Fig.2 Phylogenetic tree of 24 kinds of sequence types in this study and 14 sequence types in home and abroad

3 讨 论

MLST作为一种常用且准确的分子分型方法, 具有分辨率高、可积累、重复性强等优势得到了研究者的广泛关注[9, 10]。本研究对四川部分地区分离的48株禽源空肠弯曲菌进行MLST分型, 共获得24种不同的序列型, 其中有22株分离株属于11种新序列型, 新序列型占比为45.83%, 发现的新序列型数量所占比率与薛峰[11]、董俊[12]等从禽类肠道分离的菌株、王新等[13]从禽肉分离的菌株以及Petra等[14]自肉鸡群分离的菌株新序列型所占比率均有差异, 或与不同地区分离株的遗传进化有关。在发现的新序列型中, 眉山鹌鹑源的8株分离株均为新序列型且变异度较小, 或与菌株间的基因水平转移、代谢类型相似等因素有关, 是否如此机制还有待于进一步的研究证实后以便为更好地揭示空肠弯曲菌致病性与分子分型之间的相关性。

空肠弯曲菌可以在多种不同宿主间进行传播, 但某一类宿主所携带的空肠弯曲菌则大多归为某一序列型或克隆系, 这可能与生态适应、地理隔离、宿主相关等因素均有关[15]。本研究中的不同源分离株均存在独特型序列, 不同源分离株中的独特型序列可反映出该源分离株基因间的多样性。对分离株所含的序列型进行遗传进化分析发现, 鹌鹑源分离株V13-20与鸭源、鸡源分离株的亲缘关系较远, 表明不同家禽中的空肠弯曲菌序列型存在差异具有多样性, 分离株中出现新序列型或与空肠弯曲菌MLST数据库中的数据量缺乏有关; 分离株所含的CC-21、CC-45、CC-353、CC-354克隆系与Ragimbeau等[16]自腹泻病人中分离到的空肠弯曲菌所含的克隆系有交叉重叠, 暗示着禽类对于人类致病有潜在相关性。综上所述, 本研究结果为MLST数据库的扩充及基于该方法探究不同地域、宿主来源的空肠弯曲菌致病性与分子分型的相关性提供了参考依据。

利益冲突:无

The authors have declared that no competing interests exist.

参考文献
[1] Yoo JH, Choi NY, Bae YM, et al. Development of a selective agar plate for the detection of Campylobacter spp. in fresh produce[J]. Internation J Food Microbiol, 2014, 189(7): 67-74. DOI:10.1016/j.ijfoodmicro [本文引用:1]
[2] Epps SV, Harvey RB, Hume ME, et al. Foodborne Campylobacter: infections, metabolism, pathogenesis and reservoirs[J]. Internation J Environment Research Public Health, 2013, 10(12): 6292-6304. DOI:10.3390/ijerph10126292 [本文引用:1]
[3] Gillespie IA, Ai E. A case-case comparison of Campylobacter coli and Campylobacter jejuni infection: a tool for generating hypotheses[J]. Emerg Infect Dis, 2002, 8(9): 937-942. DOI:10.3201/eid0809.010187 [本文引用:1]
[4] Marri PR, Harris LK, Houmiel K, et al. The effect of chromosome geometry on genetic diversity[J]. Genetics, 2008, 179(1): 511-516. DOI:10.1534/genetics.107.085621 [本文引用:1]
[5] Dingle KE, Colles FM, Wareing DR, et al. Multilocus sequence typing system for Campylobacter jejuni[J]. J Clinical Microbiol, 2001, 39(1): 14-23. DOI:10.1128/JCM.39.1.14-23.2001 [本文引用:1]
[6] 王寒, 张博, 莫亚霞, . 四川部分地区禽源空肠弯曲菌ERIC-PCR分型及毒力相关基因分析[J]. 浙江农业学报, 2017, 29(9): 1465-1473. DOI:10.3969/j.issn.1004-1524.2017.09.07 [本文引用:1]
[7] Choo KH, Lee HJ, Knight NJ, et al. Multilocus sequence typing (MLST) analysis of Cand ida albicans isolates colonizing acrylic dentures before and after denture replacement[J]. Medical Mycology, 2017, 55(6): 673-679. DOI:10.1093/mmy/myw128 [本文引用:1]
[8] Sails AD, Swaminathan B, Fields PI. Clonal complexes of Campylobacter jejuni identified by multilocus sequence typing correlate with strain associations identified by multilocus enzyme electrophoresis[J]. J Clinical Microbiol, 2003, 41(9): 4058-4067. DOI:10.1128/JCM.41.9.4058-4067.2003 [本文引用:1]
[9] Olkkola S, Nyka S, Raulo S, et al. Antimicrobial resistance and multilocus sequence types of finnish Campylobacter jejuni isolates from multiple sources[J]. Zoonoses Public Health, 2016, 63(1): 10-9. DOI:10.1111/zph.12198 [本文引用:1]
[10] Kovanen SM, Kivisto RI, Rossi M, et al. A combination of MLST and CRISPR typing reveals dominant Campylobacter jejuni types in organically farmed laying hens[J]. J Applied Microbiol, 2014, 117(1): 249-257. DOI:10.1111/jam.12503 [本文引用:1]
[11] 薛峰, 徐飞, 栾军, . 多位点序列分型分析空肠弯曲菌华东动物源分离株[J]. 微生物学报, 2010, 50(3): 298-303. DOI:10.13343/j.cnki.wsxb.2010.03.015 [本文引用:1]
[12] 董俊, 韩梅, 周康, . 47株空肠弯曲菌湖北禽源株的多位点序列分型[J]. 微生物学报, 2016, 56(1): 150-156. DOI:10.13343/j.cnki.wsxb.20150175 [本文引用:1]
[13] 王新, 周婷, 孟江洪. 禽肉空肠弯曲杆菌的多位点序列分型技术研究[J]. 中国食品学报, 2010, 10(2): 180-186. DOI:10.16429/j.1009-7848.2010.02.018 [本文引用:1]
[14] Griekspoor P, Engvall EO, Olsen B, et al. Multilocus sequence typing of Campylobacter jejuni from broilers[J]. Veterinary Microbiol, 2010, 140(1/2): 180-185. DOI:10.1016/j.vetmic.2009.07.022 [本文引用:1]
[15] Karenlampi R, Rautelin H, Schonbergnorio D, et al. Longitudinal study of finnish Campylobacter jejuni and C. coli isolates from humans, using multilocus sequence typing, including comparison with epidemiological data and isolates from poultry and cattle[J]. Applied Environment Microbiol, 2007, 73(1): 148-155. DOI:10.1128/AEM.01488-06 [本文引用:1]
[16] Ragimbeau C, Schneider F, Losch S, et al. Multilocus sequence typing, pulsed-field gel electrophoresis, and fla short variable region typing of clonal complexes of Campylobacter jejuni strains of human, bovine, and poultry origins in Luxembourg[J]. Applied Environment Microbiol, 2008, 74(24): 7715-7722. DOI:10.1128/AEM.00865-08 [本文引用:1]