空肠弯曲菌(Campylobacter jejuni, C. jejuni)是一种全球范围内主要的食源性人兽共患病病原菌，可引起禽类的腹泻、肝炎等，且禽类是其主要的自然宿主之一。该菌也可引起人类腹泻，是格林巴利综合征的前趋因子^[1]。自二十世纪70年代起，空肠弯曲菌就是引起肠道感染的主要原因之一，由此引发的食品卫生安全问题日益突出，严重影响人畜健康，成为全球公共卫生的焦点^{[2, 3]}。Boes等研究表明，家禽是空肠弯曲菌最主要的传染源之一，消费禽肉与空肠弯曲菌的流行密切相关^{[4, 5]}。因此对家禽空肠弯曲菌流行情况的调查和分子特征研究是有必要的。

多位点序列分型是近年来发展很快的分子生物学分析方法,具有很高的分辨能力,既适于分子流行病学研究,也可用于分子进化学的研究。Dingle等在2001年最早提出使用MLST方法进行空肠弯曲菌的研究^[6]。2003年Manning等运用MLST对动物源和人源空肠弯曲菌进行了分子分型研究，从分子水平上阐述了两者之间的联系^[7]。国内相关研究相对较少，薛峰等^[8]和张毛军等^[9]分别用MLST方法对华东地区和华北地区空肠弯曲菌分离株进行分子分型，而华中地区暂无相关数据。

本研究使用MLST方法对湖北地区2013–2014年分离的47株禽源空肠弯曲菌进行分子分型研究，为分析我国空肠弯曲菌的流行和分子遗传特征提供科学数据。

2013–2014年湖北地区分离的47株禽源空肠弯曲菌(表1)，均经过PCR鉴定和生化鉴定。武汉的样品来源于武汉市某一大型活禽交易市场，活禽来源于武汉市周边多个养鸡场，而大冶样品来源于2个养鸡场和1个养鸭场、监利样品来源于1个养鸭场、仙桃的样品来源于3个养鸡场。

(1) 炭头孢哌酮去氧胆酸(mCCDA)培养基、抗生素添加剂及厌氧罐购自OXOID公司；微需氧产气袋购自日本MGC公司；布氏肉汤购自北京陆桥生物科技公司；(2) 细菌基因组提取试剂盒和DNA Marker购自宝生物工程(大连)有限公司，2×Taq Plus Master Mix (Dye Plus)购自南京诺唯赞生物科技有限公司。

–80 ℃保存的分离菌种复苏后接种单菌落于含10%小牛血清的布氏肉汤中微需氧条件(42 ℃，5% O₂，10% CO₂，85% N₂)培养48 h，使用细菌基因组提取试剂盒提取菌株的总DNA，–20 ℃保存备用。

参考MLST网站(http://www.mlst.net/)针对C. jejuni提供的7对PCR扩增引物序列和7对PCR测序引物序列^[6]，引物由生工生物工程(上海)股份有限公司合成。反应体系为2×Taq Plus Master Mix (Dye Plus) 12.5μL，模板DNA 2 μL (10 ng)，上下游引物各1 μL (10 μmol/L)，补水至25 μL。反应条件：(1) aspA、gltA、pgm和uncA基因反应条件为：94℃ 5 min；94 ℃ 30 s，50 ℃ 30s，72 ℃ 1 min，35个循环；72 ℃ 7 min；(2) glyA和glnA基因的反应条件为：94 ℃ 5min；94 ℃ 30s，52 ℃ 30s，72 ℃ 1 min，35个循环；72 ℃ 7 min；(3) tkt基因的反应条件为：94℃ 5 min；94 ℃ 30s，58 ℃ 30s，72 ℃ 1 min，35个循环；72 ℃ 7 min。

将PCR产物送往生工生物工程股份有限公司，使用7对测序引物进行测序。将测序结果上传到MLST数据库进行在线数据分析。获取每个菌株7个管家基因的等位基因编号(Allelic profile)，根据等位基因数得出序列型(Sequence type，ST)及克隆群(The clonal complex，CC)。与数据库比对后数据库中没有的Allelic profile和ST，则为新发现的Allelic profile和ST。根据获得的Allelic profile和ST使用START Version 2软件进行聚类分析构建系统进化树^{[10, 11, 12]}。

提取47株空肠弯曲菌分离株的基因组后，采用MLST网站提供的7对PCR扩增引物对C. jejuni基因组进行扩增，分别扩增出47株空肠弯曲菌分离株的7个管家基因，其中gltA为1012 bp，pgm为1150 bp，uncA为1120 bp，gly A为816 bp，tkt为1102 bp，aspA为899 bp，glnA为1262 bp。对PCR产物进行纯化后，使用MLST网站提供的7对测序引物对PCR产物进行双向测序。

将拼接好的序列上传MLST网站(http://www.mlst.net/)进行分析，获取相应等位基因编号和ST型及克隆群。如表2所示，本研究中共获得38个ST型，其中含25个新的ST型，2个新的等位基因编号(aspA384和glnA526)。其中32株属于以下克隆群：ST-1150 CC (4株)，ST-353 CC (8株)，ST-354 CC (3株)，ST-45 CC (2株)，ST-460 CC (1株)，ST-464 CC (9株)，ST-828 CC (2株)，ST-52 CC、ST-574 CC、ST-607 CC (各1株)，共有15株由于数据库数据有限未被分配(表2)。

47株C. jejuni多位点序列分型结果使用START Version 2软件聚类分析，制作遗传系统进化树(图1)。样本中出现的克隆群主要为ST-353 CC和ST-464 CC，大部分新发现的ST型未被分配克隆群，但与数据库已有的ST型1121、4258、2842、354和464亲缘关系上相近。本研究中包括了3株鸭源菌株(JLY1、JLY2、DYY1)和1株鸽源菌株(G2-14)，与2株鸡源菌株一起，均位于单独的一个进化分支上，与大部分鸡源菌株亲缘关系稍远。从地域来看，武汉大型活禽交易市场分离的30株空肠弯曲菌有27个不同的ST型，分布在不同的分支上，呈现较好的分子多样性。仙桃分离的8株空肠弯曲菌只有4个不同的ST型，7株处于同一分支上。而大冶分离的6株空肠弯曲菌为6个不同的ST型，无规律地分布在不同分支上。

图1 47株分离菌株的多位点序列系统进化树 Figure 1 Phylogenetic tree composed by means of the seven housekeeping genes of the MLST analysis. The dates show the name of isolates, sequence types and allelic profiles of the seven housekeeping genes (aspA, glnA, gltA, glyA, pgm, tkt, uncA), respectively.

图选项

多位点序列分型自1998年问世以来，以其高分辨率，可重复性好，可积累等优势而得到广泛应用^{[13, 14, 15]}。随着高通量测序技术以及分析软件的发展，数据库积累的加快，MLST在病原微生物分类方面的优势更加突出。传统分子分型方法有RAPD、RFLP、AFLP、PFGE等，这些分子分型方法一般采用电泳带型图谱作为分型的依据，图谱的识别需要专门的分析软件，而且研究结果在不同的实验室之间难以相互比对及共享，而MLST是一种可积累的方法，伴随着数据库的建立和使用，其结果可以在不同时间全球范围内进行数据共享^[16]。MLST数据库信息的准确性是建立在共享者正确使用的前提下，数据库要求实验者将新的型别上传，提供细菌进化相关信息，而已有型别的数据提交在分子流行病学统计上具有重要的意义^[8]。

据相关文献报道，我国华北地区分离的空肠弯曲菌克隆群有ST-353CC、ST-354CC、ST-45CC、ST-460CC、ST-52CC、ST-607CC、ST-49CC、ST-1034CC、ST-206CC、ST-257CC、ST-362CC、ST-21CC和ST-443CC，华东地区有ST-21CC、ST-49CC、ST-42CC和ST-45CC^{[8, 9]}。本研究共获得10个克隆群，与华北地区存在6个相同的克隆群，即ST-353 CC、ST-354 CC、ST-45 CC、ST-460 CC、ST-52 CC、ST-607 CC，与华东地区存在ST-45 CC相同，说明湖北地区空肠弯曲菌的流行与华东、华北地区空肠弯曲菌流行存在交叉性，这可能与武汉交通发达，人员及畜禽产品流动频繁密切相关；此外，本研究中出现的最多的克隆群ST-464CC在华东、华北地区都没出现过，而华东、华北均有的ST-21CC和ST-49CC在湖北地区没有出现过，说明我国空肠弯曲菌的流行存在地域特征性。

Manning等从人腹泻等疾病病人中分离到ST-21 CC、ST-45 CC、ST-283CC、ST-206 CC、ST-257 CC、ST-353 CC、ST-354 CC、ST-52 CC和ST-443 CC等克隆群，其中ST-21 CC和ST-45 CC数量较多，而家禽中报道较多的为ST-45 CC和ST-52 CC，人类与家禽重叠较多的为ST-45 CC^{[17, 18, 19, 20]}。本实验中分离的ST-45 CC、ST-52 CC、ST-353 CC和ST-354 CC菌株与以上研究报道的部分人源菌株具有相同的克隆群，推测可能对人存在潜在的致病性。我们也将进一步开展致病性研究，探讨空肠弯曲菌分子分型与致病性之间的关系。

用UPGMA方法对本实验得到的空肠弯曲菌ST型进行聚类分析，进化树显示新的ST型与数据库已有ST型在同一个分支。从物种上来看，鸭源JLY1、JLY2、DYY1菌株和鸽源G2-14菌株与其他大部分鸡源菌株亲缘关系稍远，显示不同家禽中空肠弯曲菌序列型存在一定的差异。从地域上来看，武汉大型活禽交易市场分离的空肠弯曲菌共占有多达27个不同的ST型，且分布在不同的分支上，呈现较好的分子多样性，可能与武汉菌株来自于不同的养殖场、不同的交易地有关。其中7株在同一分支上，类型相对比较单一，可能起源于共同的菌株。大冶分离株DYY1与武汉分离株G2-14在同一分支上，显示不同地域、不同宿主来源的菌株也存在亲缘关系很近的情况。

在本研究中出现了新的等位基因编号aspA384和glnA526以及25种新的ST型，可能是由于MLST数据库中空肠弯曲菌的数据库数据量不够大造成的。此前MLST数据库中的空肠弯曲菌样本来源主要是欧美和大洋州国家，亚洲国家的数据量相对较少。本实验的结果，不仅分析了湖北地区空肠弯曲菌的分子遗传特征，是对MLST数据库的进一步完善，也为研究不同区域不同来源的空肠弯曲菌的流行病学提供了科学数据。