布魯氏菌分離株MLVA分子分型鑒定

馬曉菁，葉 鋒，姚 剛，易新萍，劉 帥，谷文喜，吐爾洪·努爾，馬俊杰，鐘 旗

?

布魯氏菌分離株MLVA分子分型鑒定

馬曉菁1，葉 鋒1，姚 剛2，易新萍1，劉 帥2，谷文喜1，吐爾洪·努爾1，馬俊杰1，鐘 旗1

目的 對牛奶、羊奶以及羊流產胎兒中分離布魯氏菌進行分型鑒定。方法 采集新疆4個地區牛奶20份，羊奶20份，羊流產胎兒10份，分離培養后，運用VirB8-PCR和布魯氏菌分子分型PCR(AMOS-PCR)方法對分離株進行布魯氏菌屬、種的鑒定，同時采用MLVA方法對分離株進一步鑒定，并將測定結果與http://mlva.u-psud.fr/數據庫提供的布魯氏菌數據進行比較，結合軟件BioNumerics 6.6進行聚類分析。結果 6株分離株均為羊種3型布魯氏菌，與遼寧省分離羊種3型布魯氏菌在聚類分析中關系較近。結論 MLVA作為布魯氏菌基因分型鑒定的一種快速、可靠的方法，為今后布魯氏菌傳染源的追溯及防控措施的制定提供依據。

MLVA; 牛乳; 布魯氏菌; 鑒定

布魯氏菌病是由布魯氏菌引起的世界范圍內嚴重危害畜牧業和人類健康的人畜共患病。布魯氏菌屬分為9個種[1]，其中牛種和羊種布魯氏菌傳染性較強，也是我國主要流行的布魯氏菌菌種。新疆地區早在20世紀50～80年代初就是我國主要的布病疫區，90年代初新疆布病疫情得到了有效的控制，但此后布病疫情不斷回升[2]。人類可通過直接接觸患病動物感染布魯氏菌病, 患病動物的皮、毛、乳、肉以及內臟也是該病重要的傳染源[3-6]。新疆人間布病監測結果顯示2013年乳肉加工人員、畜產品銷售加工人員布病疫情上升較明顯[7]。畜間布病疫情回升不僅嚴重影響畜牧業的發展，同時危害人民健康。

布魯氏菌傳統的細菌學鑒定方法，雖可以對病畜做出確切的診斷，但鑒定時間長，且對實驗室操作人員具有較高的感染風險[8]。多位點可變數目串聯重復序列分析(Multiplelocus variable-number tandem-repeat analysis,MLVA)是近年發展起來的以PCR為基礎，根據被檢菌株散在基因組中不同位點的、變數量串聯重復序列(variable number randem repeat，VNTR)的重復單元拷貝數的差異來進行分子分型的技術。該方法具有操作簡單，重復性好，可比性強等優點[9]。本研究利用MLVA方法對牛奶、羊奶以及羊流產胎兒中分離布魯氏菌進行生物型的分析，為追溯傳染源，確定當地布魯氏菌流行趨勢制定相應的防控措施提供依據。

1 材料與方法

1.1 對象 分別從新疆不同地區采集牛奶20份、羊奶20份、羊流產胎兒10份，送實驗室-20 ℃保存。

1.2 主要試劑 布魯氏菌瓊脂培養基(BBLTMBrucella Agar)和布魯氏菌液體培養基均購自Becton and Dickinson company公司，布魯氏菌培養添加劑(OXO ID Ltd)、PCRmix、細菌基因組DNA提取試劑盒均購自北京莊盟生物有限公司，引物由上?；瞪锛夹g有限公司合成，PCR產物送北京鼎國生物有限公司測序。牛種布魯氏菌疫苗株A19，羊種布魯氏菌標準菌株16M均由新疆畜牧科學院獸醫研究所布病室保存。

1.3 方法

1.3.1 細菌培養 將乳樣混勻取300 μL接種布魯氏菌瓊脂培養基，無菌取流產胎兒脾臟研磨后接種布魯氏菌瓊脂培養基，置于37 ℃，含5% CO2培養48 h～72 h,觀察細菌生長情況。

1.3.2 布魯氏菌屬、種PCR鑒定 取布魯氏菌瓊脂培養基生長疑似菌落，接種布魯氏菌肉湯培養基中，增菌培養。按照細菌基因組DNA提取試劑盒方法提取DNA作為模板，參照鐘旗[10]等建立的VirB8-PCR方法，擴增布魯氏菌VirB8基因，對分離株進行布魯氏菌屬的鑒定。根據參考文獻[11-12]采用AMOS-PCR方法對布魯氏菌分離株進行種型的鑒定。反應條件如下：94 ℃ 5 min；94 ℃ 60 s、60 ℃ 90 s、72 ℃ 60 s，40個循環；72 ℃延伸10 min。

1.3.3 MLVA鑒定 合成引物16對Panel 1(Bruce06、Bruce08、Bruce11、Bruce12、Bruce42、Bruce43、Bruce45、Bruce55),Panel 2(Bruce18、Bruce19、Bruce21、Bruce04、Bruce07、Bruce09、Bruce16、Bruce30 )見文獻[13-14]，反應體系：PCRmix 12.5 μL，上、下游引物(10 pmol/μL)各1 μL、分離菌DNA 1 μL雙蒸水補足25 μL，反應條件：96 ℃，3 min；96 ℃ 30 s、60 ℃ 30 s、72 ℃ 30 s，35個循環；72 ℃延伸5 min。Panel 1引物擴增PCR產物1%瓊脂糖凝膠電泳，Panel 2引物擴增產物毛細管電泳，并送公司測序。計算獲得各位點的重復單元數，整理成Excel表格，將結果輸入http://mlva.u-psud.fr/，與數據庫進行比較，以確定分型，通過軟件BioNumerics 6.6聚類分析，聚類方式采用平均連鎖聚類法(UPGMA)。

2 結 果

2.1 細菌分離培養 40份奶樣、10份羊流產胎兒脾臟接種布魯氏菌瓊脂培養基培養3～4 d，觀察菌落生長情況，其中1份牛奶樣，3份羊奶樣，2份流產胎兒脾臟接種的培養基上可見圓形、無色透明、表面凸起的菌落生長。

2.2 布魯氏菌屬VirB8-PCR鑒定 挑取6個培養基上的疑似菌落，增菌，提取基因組DNA作為模板，采用VirB8-PCR方法對疑似布魯氏菌PCR鑒定。分別以牛種布魯氏菌疫苗株A19和大腸桿菌DH5α作為陽性對照和陰性對照。PCR鑒定結果可見6株分離株均擴增出750 bp左右的條帶，與陽性對照一致，而陰性對照DH5α未擴增出相應條帶，見圖1。

M：DNA DL2000 Marker;1:E.coli DH5α；2：牛種布魯氏菌標準菌A19；3-8：奶樣、流產胎兒中布魯氏菌分離株M：DNA DL2000 Marker;1:E.coli DH5α；2：Standard strain of B.abortus;3-8:Brucella isolated in milk and aborted fetuses圖1 布魯氏菌分離株VirB8-PCR鑒定Fig.1 Identification of Brucella by VirB8-PCR

2.3 布魯氏菌種型AMOS-PCR鑒定 VirB8-PCR鑒定陽性分離株采用AMOS-PCR方法進一步種型的鑒定，牛種A19(1、2、3a、4型)、羊種16M和大腸桿菌DH5α作為陽性對照和陰性對照。鑒定結果顯示6株布魯氏菌分離株均擴增出731 bp條帶，與羊種布魯氏菌(1、2、3型)擴增結果一致，見圖2。表明6株分離株為羊種布魯氏菌。

2.4 MLVA分型 分別對6株羊種布魯氏菌分離株16個位點重復單元PCR擴增,獲得相應片段，見圖3。

整理測序結果計算重復單元數，將結果輸入到http://mlva.u-psud.fr/，與數據庫進行比較，結果見表1。

M：DNA DL2000 Marker;1:E.coli DH5α；2：羊種標準布魯氏菌16M；3：牛種標準布魯氏菌A19；4-9：奶樣及流產胎兒中布魯氏菌分離株M：DNA DL2000 Marker;1:E.coli DH5α；2：Standard strain of B.melitensis;3:Standard strain of B.abortus;4-9:Brucella isolated in milk and aborted fetuses圖2 布魯氏菌分離株AMOS-PCR鑒定Fig.2 Identification of Brucella biotypes by AMOS-PCR

M：DNA DL2000 Marker;1:Bruce 11位點；2：Bruce 06位點；3：Bruce 12位點；4：Bruce 45位點；5：Bruce 08位點；6：Bruce 42位點；7：Bruce 55位點；8：Bruce 43位點M：DNA DL2000 Marker;1：Bruce 11;2:Bruce 06;3:Bruce 12;4:Bruce 45;5:Bruce 08;6:Bruce 42;7:Bruce 55;8:Bruce 43.圖3 布魯氏菌分離株MLVA位點PCR擴增Fig.3 Brucella isolates PCR amplification of MLVA

表1 6株布魯氏菌分離株MLVA分型結果Tab.1 Result of 6 Brucella strain by MLVA

樣品編號(Straincode)生物型(Biotype)B06B08B11B12B42B43B45B55B18B19B21B04B07B09B16B30NR1羊種3型(B.melitensisbiovar3)153132232420844365NR5羊種3型(B.melitensisbiovar3)1531322324208443610NR8羊種3型(B.melitensisbiovar3)1531322324208543611NR11羊種3型(B.melitensisbiovar3)153132232420844365NR24羊種3型(B.melitensisbiovar3)153132232420844367NR32羊種3型(B.melitensisbiovar3)153132232420844365

VirB8-PCR、AMOS-PCR及MLVA方法對分離株鑒定,結果表明6株分離株均為羊種3型布魯氏菌,見表2。

2.5 聚類分析 將擴增產物VNTR重復次數，經BioNumerics 6.6進行聚類分析，結果可見牛奶、羊奶以及羊流產胎兒中分離6株羊種布魯氏菌與遼寧省分離羊種3型布魯氏菌在聚類分析中關系較近,見圖4。

表2 6株布魯氏菌分離株鑒定結果
Tab.2 ientification result of 6 Brucella strain

樣品編號(Straincode)樣品來源(Samplesources)VirB8-PCR方法鑒定結果(ResultofVirB8-PCR)AMOS-PCR方法鑒定結果(ResultofAMOS-PCR)MLVA方法鑒定結果(ResultofMLVA)NR1牛奶(BovineMilk)布魯氏菌(Brucella)羊種布魯氏菌(B.melitensis)羊種3型布魯氏菌(B.melitensisbiovar3)NR5羊奶(SheepMilk)布魯氏菌羊種布魯氏菌羊種3型布魯氏菌NR8羊奶(SheepMilk)布魯氏菌羊種布魯氏菌羊種3型布魯氏菌NR11羊流產胎兒脾臟(sleepnfromabortedfetusesofsheep)布魯氏菌羊種布魯氏菌羊種3型布魯氏菌NR24羊流產胎兒脾臟(sleepnfromabortedfetusesofsheep)布魯氏菌羊種布魯氏菌羊種3型布魯氏菌NR32羊流產胎兒脾臟(sleepnfromabortedfetusesofsheep)布魯氏菌羊種布魯氏菌羊種3型布魯氏菌

圖4 6株布魯氏菌分離株MLVA聚類分析Fig.4 Dendrogram based on the MLVA genotyping assay of the 6 strains of Brucella

3 討 論

布病疫情的回升不僅與疫情程度、流行范圍有關，更重要的是對布病控制所采取有效的防治策略與措施相關[14]，針對不同地區布病流行情況因地制宜采取相應的措施能夠有效控制布病疫情。本研究通過分子生物學方法替代傳統方法完成分離株型的鑒定，與傳統的生物學鑒定方法相比操作過程簡單、快速，并降低操作人員感染布魯氏菌病的風險。多位點數目可重復序列分析(MLVA)通過分析VNTR分布進行基因分型研究。Le Flèche P和Al DahouK S[15-16]建立并完善的布魯氏菌MLVA-16位點，將16個位點分為3組，Panel 1包括8個位點，用于布魯氏菌種的比較。Panel 2A包括3個位點，Panel 2B包括5個位點，用于布魯氏菌進一步分型以及菌株間基因型的比較。楊杰[9]等初步建立了我國常見布魯氏菌MLVA分型參考標準操作方法，用于布魯氏菌鑒定、基因分型以及傳染源的追溯、布魯氏菌暴發和流行的分析。姜海[17]等通過MLVA方法結合多重PCR對24株犬種布魯氏菌遺傳關系研究，結果表明該方法可以用于布魯氏菌種(生物型)的鑒定，成為傳統表型鑒定方法的補充。毛玲玲[18]等采用MLVA方法對遼寧地區人間分離布魯氏菌33株分型分析，結果表明遼寧省布魯氏菌存在豐富的多態性，為布魯氏菌病的防治提供有力的分子生物學依據。本研究采用MLVA方法對分離布魯氏菌進行種型的鑒定，結果可見6株分離株均為羊種布魯氏菌生物3型，與遼寧省分離羊種3型布魯氏菌在聚類分析中關系較近。

新疆是我國五大牧區之一，也是布病的老疫區。我區布魯氏菌病的防治工作從1954年起開展，上世紀80～90年代畜間布病疫情大幅度下降，1993年對新疆分離布魯氏菌進行系統鑒定結果發現新疆流行的布魯氏菌有5個種8個生物型，分別為羊種布魯氏菌生物1、2、3型，其中以羊種1型最多，牛種布魯氏菌生物1、3型，其中牛種布魯氏菌生物1型最多，豬種布魯氏菌生物3型，犬種布魯氏菌、綿羊附睪種布魯氏菌[19]。近年來畜間布魯氏菌病疫情持續上升，2005—2007年家畜總布病陽性率為20世紀90年代的2倍，2010年家畜布病發病數是2006年的4.01倍[20]。本研究從牛奶、羊奶以及羊流產胎兒中分離出羊種布魯氏菌生物3型與近年來新疆地區及我國其它部分省市布魯氏菌流行菌型一致[21]，表明新疆地區布病疫情上升可能與家畜引進過程中隔離、復檢等監管措施不力有關，應加強家畜產地檢疫及流通環節的監督力度。

人間的布病來自畜間，乳和肉等畜產品是人間布魯氏菌病流行的重要傳染因子。因此，在制定合理畜間布病防治措施的同時應加強乳和肉等畜產品檢疫、監督管理。

[1] Cloeckaert A,Vizcaino N, Paquet JY, et al.Major outer membrane proteins ofBrucellaspp.past, present and future[J].Vet Microbiol,2002,90:229-247.

[2] Li L,Yan JH,Yan H,et al. Research of prevention technology and prevalent in livestock brucellosis in Xinjiang[J].Grass-Feeding Livestock,2012,(9):19-22.(in Chinese)

李玲，閆晶華，閆昊，等.新疆畜間布魯氏菌病流行情況及防制技術研究[J].草食家畜，2012(9)：19-22

[3] Adams L,Station T,NetLibrary I.Advances in brucellosis research[M].Texas: Texas A&M University,1990.

[4] Corbel MJ. Manual of diagnostic tests and vaccines for terrestrial animals[M].6thed. Paris: Office International des Epizootics,OIE.2008:36-37.

[6] Ren DK,Chang Q,Shi ML,et al. Epidemic situation and preventive measures ofBrucellain Xinjiang[J].Endem Dis Bull,2008,(5):36-37.(in Chinese)

任德坤，常青，師茂林，等.新疆布魯氏菌病疫情現狀與防治對策[J].地方病通報，2008,(5):36-37

[7] M Hasan, He HB,Tai XP,et al.Analysis of surveillance data and epidemic situation of human brucellosis in Xinjiang,2013[J]. Chin J Vector Biol Ctrl,2015,(2):86-88.(in Chinese)

木合塔爾·艾山，何海波，邰新平等. 新疆2013年人間布魯氏菌病監測結果及疫情分析[J].中國媒介生物學及控制雜志,2015(2):86-88

[8] Bricker BJ. PCR as a diagnostic tool for brucellosis[J].Vet Microbiol,2002,90:435-446.DOI: 10.1016/S0378-1135(02)00228-6

[9] Yang J,Zhao SL,Cui BY,et al. Establishment of standard operation procedures for multiple locus variable number tandem repeats analysis on typing ofBrucella[J].Dis Surveill, 2015,(2):86-88.(in Chinese)

楊杰，趙素蓮，崔步云，等. 多位點可變數目串聯重復序列分型方法檢測布魯氏菌分型標準化操作方法的建立和應用[J]. 疾病監測2012 ,27 (2 ).137-140.

[10] Zhong Q,Fan WX,Wu DL,et al.Establishment of method to detect the virulence factor ofBrucellathe VirB8-PCR assay[J]. Chin J Zoonoses,2008,24(1):50-54.(in Chinese)

鐘旗，范偉興，吳冬玲，等.布魯氏菌VirB8-PCR方法的建立[J]. 中國人獸共患病學報，2008，24(1)： 50-54.

[11] Zhong Q,Yi XP,Li B,et al.Study and application of polymerase chain reaction assay onBrucellaidentification[J].Chin J Zoonoses,2011,27(3): 241-245.(in Chinese)

鐘旗,易新萍,李博,等.布魯氏菌PCR鑒定方法的研究與應用[J]. 中國人獸共患病學報，2011，27(3): 241-245.

[12] Zhong Q,FanWX,He QN,et al.Use of the abbreviatedBrucellaAMOS-PCR assay for the species identification ofBrucella[J]. Chin J Zoonoses, 2007,23(7):683-686.(in Chinese)

鐘旗，范偉興，何倩倪，等.用AMOS—PCR對布魯氏菌種型鑒定的研究[J].中國人獸共患病學報，2007，23(7):683-686.

[13] Yang J.The establishment and application of theBrucellaMLVA and the HOOF standard operation procedures[D].Taiyuan:Shanxi Medical University,2012.(in Chinese)

楊杰.布魯氏菌MLVA和HOOF分型標準操作方法的建立及應用[D].山西醫科大學，2012.

[14] Shang DQ.Brucellaregarded as re-emerging disease and the reason[J].China J Ctrl Endem Dis,2001,16(1):29-34.(in Chinese)

尚德秋.布魯氏菌病再度肆虐及其原因[J].中國地方病防治雜志,2001,16(1):29-34

[15] Le Flèche P,Jacques I, Grayon M,et al. Evaluation and selection of tandem repeat loci for aBrucellaMLVA typing assay[J].BMC Microbiol,2006,6:9. DOI: 10.1186/1471-2180-6-9

[16] AL Dahouk S, Fleche PL,Nockler K,et al.Evaluation ofBrucellaMLVA typing for human brucellosis[J].Microbiological Methods,2007,69(1):137-145.DOI: 10.1016/j.mimet.2006.12.015

[17] Jiang H,Cui BY,Li LY,et al. TypingBrucellacanis strains gene with multiple-locus variable-number tandem-repeat analysis and multiplex PCR assay[J].Lett Biotechnol, 2009,20(3):336-338.(in Chinese)

姜海，崔步云，李蘭玉，等. 用MLVA技術和多重PCR對犬種布氏菌基因分型[J].生物技術通訊，2009,20(3):336-338

[18] Mao LL,Jiang H,Lei L,et al.Brucellosis epidemic situation and typing with MLVA in Liaoning province[J].Chin J Public Health,2012,28(12):1650-1652.(in chinese)

毛玲玲，姜海，雷露，等. 遼寧省2010 年布魯氏菌病疫情及MLVA分型[J].中國公共衛生.2012,28(12):1650-1652.

[19] Yan SD,Liu ZW,Jing GY,et al.The control of epizootic brucellosis in Xinjiang[J].Eendm Dis Bull,1993,8(4):111-113.(in Chinese)

閻守敦，劉志文，金根源. 新疆動物布魯氏菌病的防制[J].地方病通報，1993,8(4):111-113

[20] Zheng Q,Guan XD,CaoY,et al.Time and spatial distribution of human brucellosis in Xinjiang in the last 10 years[J].Bull Dis Ctrl Prev,2015, 30(3):8-12.(in Chinese)

鄭強，關曉冬，曹巖，等. 近十年新疆人間布魯氏菌病時空分布分析[J]. 疾病預防控制通報，2015，30(3)：8-12

[21] Yi XP, Li JP,Ye F,et al.Species and biotype analysis of isolated strains ofBrucellain Xinjiang[J].Xingjiang Agr Sci,2015,52(2):339-343.(in Chinese)

易新萍，李金平，葉鋒，等.新疆牛羊布魯氏菌流行株種及生物型鑒[J].2015,52(2):339-343

Zhong Qi,Email: yyyzqok@sina.com

Identification and classification ofBrucellastrain with MLVA

MA Xiao-jing1,YE Feng1,YAO Gang2,YI Xin-ping1,LIU Shuai2,GU Wen-xi1,NUER ·Tuerhong1,MA Jun-jie1, ZHONG Qi1

(1.InstituteofVeterinaryMedicine,XinjiangAcademyofAnimalScience,Urumqi830011,China;2.CollegeofAnimalMedicine,XinjiangAgriculturalUniversity,Urumqi, 830052,China)

We identified the type ofBrucellastrain from milk and abortus.The 20 portions of milk from dairy cows,20 portions of milk and 10 aborted fetuses from sheep were collected from four different area in Xinjiang.The suspected strains were identified by VirB8-PCR and AMOS-PCR.The further comparative analysis was conducted between our results andBrucellafrom database.(http://mlva.u-psud.fr/),and further phylogenetic tree was constructed by BioNumerics 6.6.Results showed that sixBrucellastrain were identified asBrucellamelitensis 3 by MLVA, and were close to theBrucellastrain belonging to Liaoning Province. Multiple Loci VNTR Analysis (MLVA )assays could be applied fast and reliably to detectBrucella.The result of MLVA provide the basis for tracing the source of infection and the scientific basis for the establishment of comprehensive measures for prevention and control of brucellosis.

MLVA;milk;Brucella;identification

鐘 旗，Email：yyyzqok@sina.com

1.新疆畜牧科學院獸醫研究所，烏魯木齊 830011；

2.新疆農業大學動物科學院，烏魯木齊 830052

10.3969/j.issn.1002-2694.2016.08.008

378

A

1002-2694(2016)08-0722-06

2016-04-20；

2016-05-29

新疆維吾爾自治區公益性科研院所基本科研業務經費項目(No.KY2014008)資助

Supported by the Basic Research Funds in Public Welfare Institution of Xinjiang(No. KY2014008)