兩株布魯菌MLST新序列型(ST)的鑒定

，，， ，，，，，

布魯菌病(以下簡稱布病)是由布魯菌(Brucella)侵入機(jī)體引起的傳染——變態(tài)反應(yīng)性人獸共患傳染病，在我國屬于國家法定管理的乙類傳染病，嚴(yán)重危害人類健康和畜牧業(yè)的發(fā)展。布魯菌是一類兼性胞內(nèi)寄生、革蘭染色陰性的短小桿菌[1]，以牛、羊、豬、犬等家畜和野生動物為主要儲存宿主，也存在于海洋動物中，可經(jīng)呼吸道、消化道、生殖系統(tǒng)黏膜以及皮膚等多種途徑引起感染。人感染布魯菌可累及全身多個組織器官，臨床表現(xiàn)為發(fā)熱、多汗、乏力、關(guān)節(jié)肌肉疼痛等[2-3]，急性期可以治愈，但不易被確診，常遷延為慢性感染，慢性期布病患者可有肝、脾及睪丸腫大，關(guān)節(jié)和脊柱強(qiáng)直，肌腱攣縮變硬等，很難治愈，如引起腦炎和心肌炎者可致命，嚴(yán)重影響人類健康[4]。動物感染后主要引起流產(chǎn)，造成畜牧業(yè)減產(chǎn)，進(jìn)而導(dǎo)致極大的經(jīng)濟(jì)損失。布魯菌易借助氣溶膠的方式傳播，呈高度傳染性(僅吸入10-100個細(xì)菌即可導(dǎo)致人類患病[5])，因此布魯菌被認(rèn)為是一種潛在的生物戰(zhàn)劑，已被列入《國際禁止生物武器公約》的核查清單，也是美國反生物恐怖襲擊中的生物劑之一[1，6]。

近年來全球布病發(fā)病率上升，被WHO稱為“再度肆虐”的傳染病。我國自20世紀(jì)90年代中期以來，布病再度流行，疫情逐年上升，防控形勢嚴(yán)峻。2000年以后，北京市布病疫情與全國一樣，報告病例數(shù)呈逐年上升趨勢，疫區(qū)人群感染率不斷升高，并向普通人群擴(kuò)散。此外，畜間疫情頻發(fā)，并有局部暴發(fā)[7]，防控形勢非常嚴(yán)峻。了解布病病原特點和遺傳特征是布病防控的基礎(chǔ)，本項研究采用基于布魯菌屬外膜蛋白31基因的PCR(BCSP31-PCR)和基于插入序列IS711的牛、羊、綿羊附睪和豬種布魯菌特異的PCR(AMOS-PCR)技術(shù)對2013年和2014年北京地區(qū)分離自布病患者血液的2株布魯菌進(jìn)行屬和種(型)鑒定，采用多位點序列分型(Multiple locus sequence typing，MLST)技術(shù)對2株布魯菌菌株的19個管家基因、1個外膜蛋白基因及1個基因間區(qū)的序列進(jìn)行測定，將各個等位基因的序列與MLST數(shù)據(jù)庫中的等位基因序列進(jìn)行比對，確定菌株的等位基因譜及9位點序列型(sequence type，ST)和21位點序列型，分析其與布魯菌不同ST型的遺傳進(jìn)化關(guān)系，為北京地區(qū)布病的防控提供科學(xué)依據(jù)。

1 材料與方法

1.1實驗試劑 菌株培養(yǎng)基采用哥倫比亞血平板，購于英國OXOID公司，貨號PB0123A。PCR試劑Platinum PCR SuperMix，購于美國Invitrogen公司，貨號11306-016。核酸自動電泳用卡夾(QIAxcel DNA Screening Kit)與QX Alignment Marker 15bp/3Kb購自德國QIAGEN公司，貨號分別為929004和929522。DL2,000bp DNA Marker購自寶生物工程(大連)有限公司，貨號為D501A。引物由生工生物工程(上海)股份有限公司合成。

1.2實驗菌株 實驗菌株BJ3與BJ91為分離自北京地區(qū)2013年和2014年布病患者血液的分離株。

1.3菌株DNA制備 將分離菌株接種于哥倫比亞血平板，置于37 ℃生化培養(yǎng)箱，培養(yǎng)72 h后，刮取1～2接種環(huán)菌苔到150 μL去離子水中，振蕩混勻制成菌懸液。再放入數(shù)控干浴器中，于100 ℃下煮沸10 min；12 000 r/min離心10 min，最后取上清液于-20 ℃冰箱中保存?zhèn)溆谩?/p>

1.4BCSP31-PCR法鑒定布魯菌屬 采用布魯菌屬特異性基因BCSP31作為定屬基因，按照參考文獻(xiàn)[8-9]提供的B4和B5引物序列和參數(shù)進(jìn)行布魯菌屬的鑒定。

1.5AMOS-PCR法鑒定布魯菌種/型 采用參考文獻(xiàn)[10-11]提供的以布魯菌屬IS711插入序列為基礎(chǔ)建立的AMOS-PCR引物及參數(shù)進(jìn)行布魯菌種/型的鑒定。

1.6MLST分型 采用參考文獻(xiàn)[12-13]提供的布魯菌19個管家基因(gap、aroA、glk、dnaK、gyrB、trpE、cobQ、int-hyp、omp25、prpE、caiA、csdB、soxA、leuA、mviM、fumC、fbaA、ddlA、putA、mutL和acnA)、1個外膜蛋白基因(omp25)和1個基因間區(qū)int-hyp作為MLST的靶標(biāo)基因，選取的靶標(biāo)基因、引物序列及擴(kuò)增長度見表1。PCR反應(yīng)條件為：95 ℃預(yù)變性5 min；94 ℃變性30 s，63 ℃退火30 s，72 ℃延伸1 min，30個循環(huán)；72 ℃延伸10 min。ddlA基因擴(kuò)增的退火溫度為55 ℃，其他條件不變。PCR產(chǎn)物采用核酸自動電泳儀電泳檢測，以DL 2 000 bp DNA Marker作為分子量標(biāo)準(zhǔn)。觀察電泳結(jié)果判斷PCR擴(kuò)增是否成功。若PCR成功，直接將電泳檢測后的樣品委托生工生物工程(上海)股份有限公司進(jìn)行純化及測序。運用MLST在線工具(網(wǎng)址https://pubmlst.org/bigsdb?db=pubmlst_brucella_seqdef&page=sequenceQuery)，分別將每株菌的21個特殊位點的測序結(jié)果與布魯菌MLST標(biāo)準(zhǔn)菌株等位基因序列進(jìn)行比對，明確序列是否與某個等位基因型的序列一致，如果一致，則定義為某個等位基因型，否則定義為一個新的等位基因型，從而獲得每株菌的等位基因譜(gap/aroA/glk/dnaK/gyrB/trpE/cobQ/int-hyp/omp25/prpE/caiA/csdB/soxA/leuA/mviM/fumC/fbaA/ddlA/putA/mutL/acnA)。將得到的等位基因譜與布魯菌MLST數(shù)據(jù)庫標(biāo)準(zhǔn)菌株的等位基因譜進(jìn)行比對，從而確定每株菌的ST型。

表1 用于擴(kuò)增和測序分析的21對引物序列
Tab.1 21 Primer sets for amplificstion and sequencing

位點推定功能引物序列5′?3′長度/bp位置gap3?磷酸甘油醛脫氫酶glyceraldehyde3?phosphatedehydrogenaseYGCCAAGCGCGTCATCGTGCGGYTGGAGAAGCCCCA589AE0172231685083?1685671aroA3?磷基酰亞胺1?羧乙烯轉(zhuǎn)移酶3?phosphoshikimate1?carboxyvinyltransferaseGACCATCGACGTGCCGGGYCATCAKGCCCATGAATTC565AE01722329974?30538glk葡萄糖激酶glucokinaseTATGGAAMAGATCGGCGGGGGCCTTGTCCTCGAAGG475#AE017224988660?989134dnaK伴侶蛋白質(zhì)chaperoneproteinCGTCTGGTCGAATATCTGGGCGTTTCAATGCCGAGCGA470AE0172232066742?2067211gyrBDNA回旋酶B亞基DNAgyraseBsubunitATGATTTCATCCGATCAGGTCTGTGCCGTTGCATTGTC469AE017223142378?141910trpE蒽酸合成酶anthranilatesynthaseGCGCGCMTGGTATGGCGCKCSCCGCCATAGGCTTC486AE0172231538194?1537709omp2525kD外膜蛋白25kDaouter?membraneproteinATGCGCACTCTTAAGTCTCGCCSAGGATGTTGTCCG490#AE017223710041?710530cobQ鈷啉胺酸合成酶cobyricacidsynthaseGCGGGTTTCAAATGCTTGGAGGCGTCAATCATGCCAGC422#AE0172231289341?1288920int?hyp假定蛋白的上游5′端upstreamandextreme5′ofhypotheticalproteinCAACTACTCTGTTGACCCGAGCAGCATCATAGCGACGGA430#AE0172231372708?1372279prpE丙酸?輔酶A連接酶propionate?CoAligaseGGTGCTGTTCACGCTGGAAAGGTTTTCGCAGGCGGCGAA468AE0172231687838?1688305caiA脂肪酰?輔酶A脫氫酶acyl?CoAdehydrogenaseTGTGTTCGGCAAGCCTTTGGGTCAAAAGACGTGCCACA449AE017224633492?633940csdB半胱氨酸脫巰基酶cysteinedesulfhydraseCGTCACTTCCTGGATCATTTCGCCACCGACGCTTATGAGAA487AE017223930109?929623soxA肌氨酸氧化酶α亞基sarcosineoxidasealphasubunitCCTCGTAAAGCGCCTTCCTGTTCGATGCCTCCACATTGG486AE017223245488?245003leuA2?異丙基蘋果酸合酶2?isopropylmalatesynthaseTCAACCGGATGAAGGAAGTCCCCTCGATAGTCTTGGTGACA482AE0172231534688?1534207表1(續(xù))位點推定功能引物序列5′?3′長度/bp位置mviM葡萄糖?果糖氧化還原酶前體glucose?fructoseoxidoreductaseprecursorATCGCCCGTTCGGTGACTGTTCGCCGTCCTTGTCC447#AE0172231990600?1990154fumC延胡索酸水合酶CfumaratehydrataseCCGACCATGTCAATATGAGCCGATATCGTTGGCGATCTTGAA452AE017224180504?180055fbaA果糖?二磷酸醛縮酶fructose?bisphosphatealdolaseCGTGAAATAACCTGATCTCACCATGCCGGTTTCAAGCGAAC458AE017224360206?360663ddlAD?丙氨酸?D?丙氨酸連接酶AD?alanine?D?alanineligaseATTTCAGTGCGCTCGAACAGGTTCTTCAATGATGAGATTAAA553AE0172231251848?1251296putA脯氨酸脫氫酶prolinedehydrogenaseGTGGGCGTGCAGCCTTTCGCCTGTGTGAGTACGAGCGG527AE017224512957?513483mutLDNA錯配修復(fù)蛋白DNAmismatchrepairproteinACATCCAAGCTGACCGACTCCCGTGCGATCACATCCGA549AE017224208009?207461acnA烏頭酸水合酶aconitatehydrataseGAAGGCCGCATCCCACTGGCGGCGAGGCAAGGTAAT490AE01722399997?99508

注：# 不同等位基因型有不同的片段長度

1.7聚類分析 應(yīng)用eBURST軟件將2株布魯菌的ST型與數(shù)據(jù)庫中的所有ST型進(jìn)行聚類分析，探討這兩個ST型與已知ST型的進(jìn)化關(guān)系。應(yīng)用MegAlign軟件對各ST型的核酸序列進(jìn)行聚類分析，探討各菌株間的進(jìn)化關(guān)系。

2 結(jié) 果

2.1BCSP31-PCR鑒定 BCSP31-PCR方法對分離株和對照菌株核酸進(jìn)行擴(kuò)增，產(chǎn)物采用核酸自動電泳儀電泳進(jìn)行檢測，結(jié)果顯示兩株分離菌株及陽性對照菌株核酸均擴(kuò)增出布魯菌屬特異的223 bp的條帶，而陰性對照無條帶出現(xiàn)。

2.2AMOS-PCR檢測結(jié)果 AMOS-PCR方法對分離株和對照菌株核酸進(jìn)行擴(kuò)增，產(chǎn)物采用核酸自動電泳儀電泳進(jìn)行檢測，結(jié)果顯示陽性羊種對照菌株核酸可擴(kuò)增出178 bp的條帶和羊種布魯菌特異的731 bp的條帶，2株分離菌株只有178 bp條帶，此外無其他條帶，而陰性對照無條帶出現(xiàn)。

2.3MLST分型結(jié)果 2株菌的19個管家基因、1個外膜蛋白基因(omp25)和1個基因間區(qū)int-hyp的擴(kuò)增片段測序結(jié)果與與布魯菌MLST數(shù)據(jù)庫等位基因信息比對，得到這2株分離株的21個位點(gap/aroA/glk/dnaK/gyrB/trpE/cobQ/int-hyp/omp25/prpE/caiA/csdB/soxA/leuA/mviM/fumC/fbaA/ddlA/putA/mutL/acnA)的等位基因序列號分別為2、1、2、2、1、3、1、4、1、13、2、2、2、2、1、1、3、7、6、2、3和2、1、2、2、1、3、1、4、29、13、2、2、2、2、1、1、3、7、6、2、3(表2，表3)。比對結(jié)果顯示，這2株布魯菌的9位點等位基因編號和21位點基因編號均為新組合，即新ST。

我們將這2株布魯菌的資料及其各等位基因序列編號通過互聯(lián)網(wǎng)錄入布魯菌MLST數(shù)據(jù)庫，已被接收并確認(rèn)為新ST，命名為ST74和ST75(9位點ST)與ST108和ST109(21位點ST)。

表2 2株分離株的等位基因譜(9位點ST)
Tab.2 Allele profile of 2 isolates (9 loci ST)

分離株STgaparoAglkdnaKgyrBtrpEcobQInt?hypomp25BJ3ST74212213141BJ91ST752122131429

表3 2株分離株的等位基因譜(21位點ST的其他12位點)
Tab.3 Allele profile of 2 isolates (addtional 12 loci of 21 loci ST)

分離株STprpEcaiAcsdBsoxAleuAmviMfumCfbaAddlAputAmutLacnABJ3ST1081322221137623BJ91ST1091322221137623

2.4聚類分析結(jié)果 聚類分析顯示，迄今已發(fā)現(xiàn)的73個9位點MLST分型的ST型所包含的布魯菌10個種呈現(xiàn)出相同種ST型菌株聚集較緊密，不同種之間遺傳距離及親緣關(guān)系較遠(yuǎn)的現(xiàn)象，這2株菌的9位點ST與其它幾個均屬于牛種布魯菌的ST型聚在一起，與ST2高度同源，遺傳關(guān)系較近，ST74與ST2只有1個等位基因int-hyp的差異，ST75與ST2有2個等位基因int-hyp、omp25的差異，提示ST74、ST75與ST2的親緣關(guān)系較近，ST2可暫定為1個可行度較高的始祖ST型(圖1)。運用MegAlign軟件(By ClustalW Method)對各ST型9個位點的核酸序列進(jìn)行聚類分析，結(jié)果顯示ST74、ST75型與屬于牛種布魯菌的其他18個ST型聚在一起(圖2)，對各ST型21個位點的核酸序列進(jìn)行聚類分析，結(jié)果顯示ST108、ST109型與屬于牛種布魯菌的其他30個ST型聚在一起(圖3)，提示這2株菌為牛種布魯菌，與AMOS-PCR結(jié)果吻合。

圖1 2種新ST型與已知布魯菌ST型的eBURST分析結(jié)果Fig.1 eBURST analysis of 2 novel STs and the known STs

圖2 2種新ST型與已知布魯菌ST型間的進(jìn)化關(guān)系(9位點ST)Fig.2 Phylogenetic tree of 2 novel STs and the known STs (9 loci ST)

圖3 2種新ST型與已知布魯菌ST型間的進(jìn)化關(guān)系(21位點ST)Fig.3 Phylogenetic tree of 2 novel STs and the known STs (21 loci ST)

3 討 論

1985年WHO根據(jù)在細(xì)菌培養(yǎng)過程中的生化以及噬菌體裂解和血清凝集特性，將布魯菌分為羊種布魯菌、牛種布魯菌、豬種布魯菌、綿羊附睪種布魯菌、沙林鼠種布魯菌和犬種布魯菌共6個種19個生物型，此為布魯菌的經(jīng)典分型。之后又陸續(xù)發(fā)現(xiàn)了4種新型布魯菌。目前認(rèn)為，布魯菌有10個種，它們分別為馬耳他布魯菌(羊種布魯菌B.melitensis，3個生物型)、流產(chǎn)布魯菌(牛種布魯菌B.abortus，8個生物型)、豬種布魯菌(B.suis，5個生物型)、綿羊附睪種布魯菌(B.ovis，1個生物型)、沙林鼠種布魯菌(B.eotomae，1個生物型)、犬種布魯菌(B.canis，1個生物型)，從海洋哺乳動物如海豹、鯨、海豚等中分離出的鯨型布魯菌(B.ceti)[14]和鰭型布魯菌(B.pinnipedialis)[15]，從野鼠中分離出的田鼠型布魯菌(B.microti)[16]，在人乳房植入體中發(fā)現(xiàn)的湖浪布魯菌(B.inopinata)[17]。其中4個種(羊種B.melitensis，牛種B.abortus，豬種B.suis和犬種B.canis)對人致病，致病作用大小順序一般是羊種大于牛種、豬種，大于犬種。布魯菌因其生物種、型和菌株的不同，致病力各異，對人致病力最強(qiáng)的為羊種菌，其次為豬種菌，牛種菌對人致病力弱，但亦時有病例發(fā)生。因此，了解布病病原學(xué)特征對布病的有效預(yù)防和控制至關(guān)重要。本研究將來源于患者的2株菌株采用屬特異性BCSP31-PCR進(jìn)行鑒定，結(jié)果顯示這兩株分離的可疑菌株均為布魯菌。AMOS-PCR結(jié)果只有178 bp條帶而無其他大小的目的條帶出現(xiàn)，提示分離菌株為為非羊種、非牛種1、2、4型、非豬種1型、非綿羊附睪種布魯菌。

傳統(tǒng)血清學(xué)方法鑒定只能將布魯菌鑒定到種(型)，無法對菌株進(jìn)行溯源分析。隨著分子生物學(xué)的發(fā)展，一些分子分型方法不斷出現(xiàn)，如PCR產(chǎn)物限制性片段長度多態(tài)性(PCR-RFIP)、隨機(jī)擴(kuò)增多態(tài)性DNA(RAP-PCR)、擴(kuò)增片斷長度多態(tài)性分析(AFLP)、脈沖場凝膠電泳(pulsed-field gel electrophoresis，PFGE)、多位點可變數(shù)目串聯(lián)重復(fù)序列分析(multi-locus variable-Humber tandem repeat analysis，MLVA-16)技術(shù)和多位點序列分型(multi-locus sequence typing，MLST)技術(shù)等，但是PCR-RFIP、RAP-PCR和AFLP分辨率較低，在分子流行病學(xué)中應(yīng)用價值十分有限。由于布魯菌屬種間DNA同源性高達(dá)90%，需要高分辨力的分型方法用來提供疾病傳播模式和分子流行病學(xué)信息。PFGE、MLVA和MLST在布病病的病原鑒定、傳染源追蹤與病原體的溯源等方面發(fā)揮了重要作用。但PFGE實驗過程中活菌操作繁多，實驗室感染風(fēng)險較大。相比MLVA-16，MLST在反映菌株遺傳特征方面更加敏感[18]，分辨率高，重復(fù)性好，數(shù)據(jù)可實現(xiàn)網(wǎng)絡(luò)共享，不同的實驗室之間具有良好的可比性[19]。2007年Whatmore等[12]根據(jù)7個管家基因、1個外膜蛋白基因及1個基因間區(qū)設(shè)計引物，對160株布魯菌進(jìn)行MLST分析，試驗將160株布氏菌分為27個序列型。為了提高分辨率，Whatmore[13]等又新增了12個靶標(biāo)基因，加上之前的9個特殊位點，形成了BruMLSA21分型方案。對收集自全球的500多株布魯菌進(jìn)行MLST分析，結(jié)果分了101個序列型。此后不斷有新的ST型被發(fā)現(xiàn)，數(shù)據(jù)庫不斷豐富，目前9位點序列型有73個ST型，21位點序列型有107個ST型。因此，本研究中我們采用布魯菌MLST分子分型方法對2013年和2014年分離自北京地區(qū)布病患者的2株菌株進(jìn)行分子分型分析并與已發(fā)現(xiàn)的ST型進(jìn)行比對，結(jié)果發(fā)現(xiàn)2株菌株的等位基因譜與目前已有的ST型的等位基因譜不完全一致，為新的ST型。

MLST數(shù)據(jù)庫(https://pubmlst.org/brucella/)是國際上公認(rèn)的登載布魯菌序列分型信息的重要數(shù)據(jù)庫，注冊登錄后將發(fā)現(xiàn)的2種新ST型菌株的等位基因譜信息提交，被接受并予以命名，現(xiàn)該庫關(guān)于Brucella的9位點ST更新為75種，21位點ST更新為109種。本研究中2株菌的ST分別被命名為ST74、ST75(9位點MLST)和ST108、ST109(21位點MLST)。ST74與ST2只有1個等位基因int-hyp的差異，ST75與ST2有2個等位基因int-hyp、omp25的差異，見表2。為了進(jìn)一步分析新發(fā)現(xiàn)的ST型，利用eBURST工具分析目前布魯菌數(shù)據(jù)庫中所有的73個ST型和新發(fā)現(xiàn)的ST型，并繪制出各ST型之間的關(guān)系圖。eBURST(Based Upon Related Sequence Type)是一種基于相關(guān)序列型的聚類分析方法。根據(jù)細(xì)菌的所有序列型的等位基因譜將細(xì)菌分為不同的組(group)，計算每個序列型的單位點變異型SLV(Single Locus Variants)個數(shù)、雙位點變異型DLV(Double Locus Variants)個數(shù)，鑒定可能的祖先序列型等。當(dāng)2個ST型之間只有1個等位基因數(shù)字不同時，這2個ST型互為單位點變異型SLV；如果有2個等位基因數(shù)字不同，這2個ST型則互為雙位點變異型DLV[20]。互為單位點變異型的序列型構(gòu)成一組，若該組有足夠多的ST型且其中的1個ST擁有最多數(shù)量的SLV，該ST型則可以暫定為1個可行度較高的始祖ST型。eBURST分析顯示(圖1)：迄今已發(fā)現(xiàn)的73個ST型所包含的布魯菌10個種呈現(xiàn)出相同種菌株的ST型聚集較緊密，而不同種之間遺傳距離及親緣關(guān)系較遠(yuǎn)的現(xiàn)象；以ST2、ST8、ST15、ST17、ST59各自為中心聚集在一起的ST型的親緣關(guān)系較近，這樣的菌株稱為克隆群(或稱克隆復(fù)合體clonal complexes)，分別為牛種、羊種、豬種、豬種和鯨型。ST74、ST75，ST108、ST109與其它幾個均屬于牛種布魯菌的ST型聚在一起，與ST2高度同源，遺傳關(guān)系較近，從遺傳學(xué)水平證實了這2株布魯菌為兩種新ST型牛種布魯菌。ST2可暫定為1個可行度較高的始祖ST型。運用MegAlign軟件(By Clustal W Method)對各ST型9個位點的核酸序列進(jìn)行聚類分析，探討不同ST型間的進(jìn)化關(guān)系。結(jié)果顯示(圖2)：羊種布魯菌幾乎所有的ST型(除ST70以外)均處于同一個進(jìn)化分支上，表現(xiàn)的高度集中，說明羊種布魯菌是一種高度保守的布魯菌，不同ST型間變異少，而國內(nèi)發(fā)現(xiàn)最多羊種ST型為ST8[21-22]。新發(fā)現(xiàn)的ST74、ST75型與ST2型遺傳關(guān)系很近，處于同一個進(jìn)化層次。ST74、75型與屬于牛種布魯菌的其他18個ST型聚在一起，21位點ST型的進(jìn)化樹分析顯示ST108、109型與屬于牛種布魯菌的其他30個ST型聚在一起(圖3)，提示這2株分離株為牛種布魯菌，與AMOS-PCR結(jié)果吻合。

近年來我國布病再度流行，防控形勢依然嚴(yán)峻。2000年以后，北京市布病疫情與全國一樣，報告病例數(shù)呈逐年上升，并由高危人群并向普通人群擴(kuò)散。此外，畜間疫情頻發(fā)，并有局部暴發(fā)[7]，防控形勢非常嚴(yán)峻。本研究結(jié)果顯示MLST分型具有較好的分辨率，能夠與其它方法互為補(bǔ)充，對布病的感染溯源具有重要作用，研究結(jié)果不僅充實了布魯菌MLST數(shù)據(jù)庫，也為北京地區(qū)布病的防控提供了科學(xué)依據(jù)，在保護(hù)人群健康和保障畜牧業(yè)發(fā)展方面具有重要意義。

參考文獻(xiàn)：

[1] Doganay GD, Doganay M.Brucellaas a potential agent of bioterrorism[J]. Recent Pat Antiinfect Drug Discov, 2013, 8(1): 27-33. DOI: 10.2174/1574891X11308010006

[2] Troy SB, Rickman LS, Davis CE. Brucellosis in San Diego: epidemiology and species-related differences in acute clinical presentations[J]. Medicine (Baltimore), 2005, 84(3): 174-187.

[3] Hasanjani Roushan MR, Ebrahimpour S. Human brucellosis: An overview[J]. Caspian J Intern Med, 2015, 6(1): 46-47. DOI: 10.1097/01.md.0000165659.20988.25

[4] Franco MP, Mulder M, Gilman RH, et al. Human brucellosis[J]. Lancet Infect Dis, 2007，7(12): 775-786. DOI: 10.1016/S1473-3099(07)70286-4

[5] Bellamy RJ, Freedman AR. Bioterrorism[J]. Qjm, 2001, 94(4): 227-234.

[6] Santis RD, Ciammaruconi A, Pomponi A, et al.Brucella: Molecular diagnostic techniques in response to bioterrorism threat [J]. J Bioterr Biodef, 2011, S2: 004. DOI: 10.4172/2157-2526.S2-004

[7] 李錫太, 何戰(zhàn)英, 陳艷偉,等. 運用集中度與圓形分布法分析北京市布魯氏菌病時間分布特征[J]. 公共衛(wèi)生與預(yù)防醫(yī)學(xué), 2016, 27(1):20-23.

[8] Baily GG, Krahn JB, Drasar BS et al. Detection ofBrucellamelitensisandBrucellaabortusby DNA amplification[J]. J Trop Med Hyg, 1992, 95(4): 271-275.

[9] Garshasbi M, Ramazani A, Sorouri R, et al. Molecular detection ofBrucellaspecies in patients suspicious of brucellosis from Zanjan, Iran[J]. Braz J Microbiol, 2014, 45(2): 533-538.

[10] Bricker BJ, Halling SM. Differentiation ofBrucellaabortusbv. 1, 2, and 4,Brucellamelitensis,Brucellaovis, andBrucellasuisbv. 1 by PCR[J]. J Clin Microbiol, 1994, 32(11): 2660-2666.

[11] Kang YX, Li XM, Piao DR, et al. Typing discrepancy between phenotypic and molecular characterization revealing an emerging biovar 9 variant of smooth phage-resistantB.abortusstrain 8416 in China[J]. Front Microbiol, 2015, 6(514): 1375. DOI: 10.3389/fmicb.2015.01375

[12] Whatmore AM, Perrett LL, MacMillan AP. Characterisation of the genetic diversity ofBrucellaby multilocus sequencing[J]. BMC Microbiol, 2007, 7(1): 34. DOI: 10.1186/1471-2180-7-34

[13] Whatmore AM, Koylass MS, Muchowski J, et al. Extended multilocus sequence analysis to describe the global population structure of the genusBrucella: phylogeography and relationship to biovars[J]. Front Microbiol, 2016, 7: 2049. DOI: 10.3389/fmicb.2016.02049

[14] Foster JT, Okinaka RT, Svensson R, et al. Real-time PCR assays of single-nucleotide polymorphisms defining the major Brucella clades[J]. J Clin Microbiol, 2008, 46(1): 296-301. DOI: 10.1128/JCM.01496-07

[15] Foster G, Osterman BS, Godfroid J, et al.Brucellaceti sp. nov. andBrucellapinnipedialissp. nov. forBrucellastrains with cetaceans and seals as their preferred hosts[J]. Int J Syst Evol Microbiol, 2007, 57(Pt 11): 2688-2693. DOI: 10.1099/ijs.0.65269-0

[16] Scholz HC, Hubalek Z, Sedlacek I, et al.Brucellamicrotisp. nov., isolated from the common vole Microtus arvalis[J]. Int J Syst Evol Microbiol, 2008, 58(Pt 2): 375-382. DOI: 10.1099/ijs.0.65356-0

[17] Scholz HC, Nockler K, Gollner C, et al.Brucellainopinatasp. nov., isolated from a breast implant infection[J]. Int J Syst Evol Microbiol, 2010, 60(Pt 4): 801-808. DOI: 10.1099/ijs.0.011148-0

[18] 李沛麗, 李世軍, 劉英,等. 貴州省兩株山羊源羊種布魯菌MLST分析[J]. 中國人獸共患病學(xué)報, 2016, 32(6):547-552.DOI：10.3969/j.issn.1002-2694.2016.06.008

[19] Vázquez JA, Berrón S. Multilocus sequence typing: the molecular marker of the Internet era[J]. Enfermedades Infecciosas Y Microbiología Clínica, 2004, 22(2): 113-120. DOI: 10.1016/S0213-005X(04)73045-1

[20] 趙宏, 諶志強(qiáng), 李君文.不同地理來源的食品中大腸埃希氏菌eBURST分析[J]. 解放軍預(yù)防醫(yī)學(xué)雜志, 2012, 30(5):325-328.

[21] Ke Y, Yuan X, Zhen Q, et al. Genome sequence ofBrucellamelitensisS66, an isolate of sequence type 8, prevalent in China[J]. J Bacteriol, 2012, 194(19): 5451. DOI: 10.1128/JB.01202-12

[22] Chen Y, Ke Y, Wang Y, et al. Changes of predominant species/biovars and sequence types ofBrucellaisolates, Inner Mongolia, China[J]. BMC Infect Dis, 2013, 13(1): 514. DOI: 10.1186/1471-2334-13-514