胡毅翔，張歡歡，余陳歡，金曉音，應(yīng)華忠

(浙江省醫(yī)學(xué)科學(xué)院 浙江省實(shí)驗(yàn)動(dòng)物與安全性研究重點(diǎn)實(shí)驗(yàn)室， 杭州 310013)

?

技術(shù)方法

基于高通量測(cè)序技術(shù)的實(shí)驗(yàn)動(dòng)物沙門氏菌檢測(cè)方法的建立與評(píng)價(jià)

胡毅翔，張歡歡，余陳歡，金曉音，應(yīng)華忠

(浙江省醫(yī)學(xué)科學(xué)院 浙江省實(shí)驗(yàn)動(dòng)物與安全性研究重點(diǎn)實(shí)驗(yàn)室， 杭州 310013)

目的 建立基于高通量測(cè)序技術(shù)的實(shí)驗(yàn)動(dòng)物沙門氏菌檢測(cè)方法，應(yīng)用于實(shí)驗(yàn)動(dòng)物沙門氏菌的檢測(cè)。方法 提取小鼠糞便DNA樣本，分別針對(duì)16S rDNA區(qū)域、23S rDNA區(qū)域、16S～23S rDNA IS、23S～5S rDNA IS、gyrB優(yōu)選區(qū)域設(shè)計(jì)通用引物，對(duì)各個(gè)引物進(jìn)行測(cè)試分析，優(yōu)化擴(kuò)增條件并建庫(kù)，利用Illumina高通量測(cè)序技術(shù)檢測(cè)區(qū)分42個(gè)樣本中的沙門氏菌，評(píng)價(jià)該方法的特異性和穩(wěn)定性。結(jié)果 篩選發(fā)現(xiàn)沙門氏菌的菌種分類優(yōu)選區(qū)域?yàn)間yrB基因，gyrB基因引物序列為FP5’-AACCACCGCAATCAGACCTT3’,FP5’-AGCCACGAAACCTTCACYA-3’。對(duì)引物進(jìn)行優(yōu)化，確定最佳擴(kuò)增條件及樣品上樣量并正式建庫(kù)，通過高通量測(cè)序和序列分析能檢測(cè)出42個(gè)樣品中極微量的沙門氏菌，檢測(cè)方法穩(wěn)定，靈敏度高，檢測(cè)限可達(dá)0～102的cfu。結(jié)論 本實(shí)驗(yàn)利用高通量測(cè)序技術(shù)，建立了一套完整檢測(cè)實(shí)驗(yàn)動(dòng)物沙門氏菌的生物體系，能檢測(cè)出實(shí)驗(yàn)動(dòng)物體內(nèi)極微量的沙門氏菌，檢測(cè)方法穩(wěn)定性好，靈敏度高，可沿用至其他種類病原微生物的檢測(cè)。

沙門氏菌；實(shí)驗(yàn)動(dòng)物；高通量測(cè)序

沙門氏菌(Salmonellae)是一類可以引起沙門氏菌病的典型人畜共患病原菌，屬腸桿菌科 (Enterobacteriaceae)，兼性厭氧型革蘭氏陰性菌[1，2]。自1885年Salmon等分離得到豬霍亂沙門氏菌以來，目前已檢出沙門氏菌血清型約2 500余種，其中我國(guó)發(fā)現(xiàn)的100多種[3]。沙門氏菌在分類學(xué)上被分為兩個(gè)種，分別為Salmonellaenterica和Salmonellabongori，可由多種途徑向人類進(jìn)行傳播，進(jìn)入人體后在腸道內(nèi)寄生并繁殖，所產(chǎn)生的沙門氏菌內(nèi)毒素能進(jìn)入血液循環(huán)，引發(fā)腸道黏膜腫脹、脫落、甚至潰瘍、壞死[4]。在實(shí)驗(yàn)室研究方面，攜帶或者潛在攜帶沙門氏菌的實(shí)驗(yàn)動(dòng)物不僅對(duì)實(shí)驗(yàn)人員造成潛在的傳染隱患，也會(huì)對(duì)后續(xù)科學(xué)實(shí)驗(yàn)產(chǎn)生較大的影響，因此對(duì)于實(shí)驗(yàn)動(dòng)物的沙門氏菌的檢測(cè)尤為重要[5]。

生化方法和免疫學(xué)方法是目前檢測(cè)沙門氏菌的主要方法，但這類方法檢驗(yàn)程序復(fù)雜繁瑣，耗時(shí)長(zhǎng)，部分檢驗(yàn)過程及結(jié)果判定存在實(shí)驗(yàn)人員的主觀判斷，因此新的檢測(cè)方法的引入勢(shì)在必行[6]。分子生物學(xué)方法具有快速、準(zhǔn)確的優(yōu)勢(shì)，現(xiàn)已開始被大量的應(yīng)用于檢測(cè)病原微生物，主要涉及到的分子生物學(xué)技術(shù)主要包括PCR技術(shù)、基因芯片和環(huán)介導(dǎo)等溫?cái)U(kuò)增技術(shù)(LAMP)等[7，8]。近年來，高通量測(cè)序技術(shù)開始引入病原微生物檢測(cè)領(lǐng)域，實(shí)現(xiàn)了更為高效的檢測(cè)手段。本次研究主要采用的是Illumina的Miseq高通量測(cè)序平臺(tái)，對(duì)于病原微生物檢測(cè)該平臺(tái)具有周期短、通量高、成本低、準(zhǔn)確率高等優(yōu)點(diǎn)[9，10]。本實(shí)驗(yàn)以分子生物學(xué)技術(shù)方法為基礎(chǔ)，采用Miseq高通量測(cè)序技術(shù)檢測(cè)實(shí)驗(yàn)室用小鼠所攜帶的沙門氏菌，通過本次研究，建立了一套完整檢測(cè)實(shí)驗(yàn)動(dòng)物沙門氏菌的技術(shù)體系。

1 材料和方法

1.1 動(dòng)物

SPF級(jí)ICR小鼠42只，體重18～20 g，由浙江省醫(yī)學(xué)科學(xué)院實(shí)驗(yàn)動(dòng)物中心提供 (實(shí)驗(yàn)動(dòng)物生產(chǎn)許可證號(hào)為 SCXK(浙) 2008-0034，使用許可證號(hào)為SYXK (浙) 2008-0114。

1.2 菌株

鼠傷寒沙門氏菌凍干株(編號(hào)： TA98)，甲型副傷寒沙門氏菌凍干株(編號(hào)： CMCC(B)50093)，腸炎沙門氏菌凍干株(編號(hào)： TDF918)，志賀菌凍干株(編號(hào)： TG212)為本實(shí)驗(yàn)室冷凍保存。

1.3 藥品和試劑

PremixExTaqTMHot Start Version(日本Takara公司，批號(hào)RR030A)，DL2000 DNA Marker(日本Takara公司，批號(hào)3427)，Axygen PCR Clean Up Kit(美國(guó)Axygen公司，批號(hào)AP-PCR-500)，QubitdsDNA HS Assay Kit(美國(guó)Invitrogen公司，批號(hào)Q32851)，PrimeScript? RT reagent Kit(日本Takara公司, 批號(hào)DRR037A)。高通量測(cè)序試劑盒(Illumina，美國(guó))，Trizol reagent(Invitrogen，美國(guó)，批號(hào)Q35438)。

1.4 主要儀器

Bioanalyzer 2100 生物分析儀(Agilent，美國(guó))；Hiseq 2500 高通量測(cè)序儀(Illumina，美國(guó))；熒光定量 PCR 儀(Bio-Rad，美國(guó))；Qubit 2.0核酸定量?jī)x(Life Technologies，美國(guó))；組織研磨儀(TissueLyser，Qiagen，德國(guó))。

1.5 實(shí)驗(yàn)方法

1.5.1 小鼠沙門氏菌模型的建立及樣品收集： 小鼠適應(yīng)性飼養(yǎng)3 d后，隨機(jī)分為6組： 鼠傷寒沙門氏菌組(5只)、甲型副傷寒沙門氏菌組(5只)、志賀桿菌組(5只)、腸炎沙門氏菌和甲型副傷寒沙門氏菌(5只)、腸炎沙門氏菌組(18只)、陰性對(duì)照組(4只)。鼠傷寒沙門氏菌組每只小鼠灌服鼠傷寒沙門氏菌，甲型副傷寒沙門氏菌組每只小鼠灌服鼠傷寒沙門氏菌0.2 mL (5×103cfu/ mL)，志賀桿菌組每只小鼠灌服鼠傷寒沙門氏菌0.2 mL (5×103cfu/ mL)，腸炎沙門氏菌和甲型副傷寒沙門氏菌每只小鼠灌服腸炎沙門氏菌和甲型副傷寒沙門氏菌0.2 mL (5×103cfu/ mL)，腸炎沙門氏菌組每只小鼠灌服鼠傷寒沙門氏菌0.2 mL(0 cfu/mL、5×101cfu/mL、5×102cfu/mL、5×103cfu/mL、5×104cfu/mL、5×105cfu/mL)，每個(gè)濃度各3只，各組小鼠每隔2周灌胃1次，連續(xù)3次。取60 mg小鼠糞樣置于2 mL離心管中，采用PowerFecal DNA Isolation Kit 試劑盒(MoBio，美國(guó)) 抽提100 mg小鼠糞便樣本中微生物的總 DNA, 操作步驟按照說明書進(jìn)行，完成基因組DNA抽提后，放-20℃保存?zhèn)溆?，樣品信息見?。

表1 樣品信息

1.5.2 菌種分類優(yōu)選區(qū)域選擇： 通過文獻(xiàn)查找，針對(duì)沙門氏菌，尋找在研究范圍內(nèi)的幾個(gè)菌種相互之間差異變化大的區(qū)域。NCBI數(shù)據(jù)庫(kù)中查找沙門氏菌典型菌株，以鼠傷寒沙門氏菌(Salmonellaentericaserovar Typhimurium)典型菌株LT2、ST4/74等作為研究對(duì)象。通過序列比對(duì)，去除經(jīng)比對(duì)菌種間同源關(guān)系極高的區(qū)域，最終篩選得到5大類區(qū)域： ①16S rDNA區(qū)域，②23S rDNA區(qū)域，③16S～23S rDNA ISR，④23S～5S rDNA ISR，⑤gyrB優(yōu)選區(qū)域。

1.5.3 通用引物設(shè)計(jì)及優(yōu)化： 針對(duì)①16S rDNA區(qū)域，②23S rDNA區(qū)域③16S～23S rDNA ISR，④23S～5S rDNA ISR，⑤gyrB優(yōu)選區(qū)域設(shè)計(jì)通用引物，通用引物在多個(gè)菌種中都能PCR擴(kuò)增得到產(chǎn)物。對(duì)設(shè)計(jì)的引物進(jìn)行引物測(cè)試，淘汰不能擴(kuò)增的引物。在原有的通用引物基礎(chǔ)上加上建庫(kù)和測(cè)序需用到的接頭序列(Linker)，區(qū)分樣本使用的正反向Barcode序列(Index)，增加終文庫(kù)均一性和豐富度的連接堿基(heterogeneity spacer)，如此設(shè)計(jì)的引物長(zhǎng)度在90 bp～97 bp之間，該優(yōu)化引物可直接用于建庫(kù)，實(shí)際的建庫(kù)工作需要一步PCR反應(yīng)完成。優(yōu)化引物設(shè)計(jì)圖如圖1所示：

圖1 優(yōu)化引物結(jié)構(gòu)組成圖Fig.1 Diagram of optimization of the primer structure

1.5.4 沙門氏菌建庫(kù)PCR條件摸索： 建庫(kù)PCR采用的是90～97 bp長(zhǎng)度的長(zhǎng)鏈引物，這與普通16～25 bp的引物PCR擴(kuò)增條件會(huì)有很多不同，長(zhǎng)鏈引物擴(kuò)增條件比普通引物擴(kuò)增的條件要苛刻。建庫(kù)PCR條件摸索： 模板上樣量、退火溫度、循環(huán)、延伸時(shí)間，其中主要考查上樣量和退火溫度。(1)退火溫度： 考察的退火溫度為56℃、54℃、52℃。2) 模板上樣量： 針對(duì)模板上樣量設(shè)計(jì)11個(gè)模板濃度梯度實(shí)驗(yàn)，模板信息如表1所示。上樣量共分為12個(gè)梯度，依次為500、200、100、50、25、10、1 ng、100、10、1，100，10 fg。

1.5.5 高通量測(cè)序體系的特異性和敏感性檢測(cè)： 42個(gè)檢測(cè)樣本的建庫(kù)后采用 Illumina Miseq 高通量測(cè)序平臺(tái)進(jìn)行高通量測(cè)序。隨后對(duì)高質(zhì)量序列進(jìn)行提取, 提取方法: Miseq 測(cè)序得到的 PE reads 首先根據(jù) overla關(guān)系進(jìn)行拼接, 將成對(duì)的 reads 拼接成一條序列, 同時(shí)對(duì)reads 質(zhì)量和 merge 的效果進(jìn)行質(zhì)量質(zhì)控和過濾, 去除序列末端的后引物和接頭序列、多堿基 N、polyA/T尾巴及低質(zhì)量堿基; 去除所得序列的 barcode 標(biāo)簽序列、前引物序列；丟棄長(zhǎng)度短于 200 bp、模糊堿基數(shù)>0、序列平均質(zhì)量低于20的序列。高通量測(cè)序結(jié)果與樣本信息進(jìn)行對(duì)比，評(píng)價(jià)高通量測(cè)序體系檢測(cè)沙門氏菌的特異性和敏感性。

2 結(jié)果

2.1 通用引物設(shè)計(jì)及引物測(cè)試

針對(duì)①16S rDNA區(qū)域，②23S rDNA區(qū)域, ③16S～23S rDNA ISR，④23S～5S rDNA ISR，⑤gyrB優(yōu)選區(qū)域各設(shè)計(jì)3對(duì)通用引物。經(jīng)過測(cè)試，16S rDNA、23S rDNA、gyrB各有一對(duì)引物通過測(cè)試。但分析后發(fā)現(xiàn)16S rDNA、23S rDNA的引物實(shí)際建庫(kù)測(cè)序的結(jié)果與設(shè)計(jì)時(shí)的有偏差，它無法區(qū)分樣品中的腸炎沙門氏菌。最后確定沙門氏菌的菌種分類優(yōu)選區(qū)域?yàn)間yrB基因引物(表2)。gyrB基因通用引物引物測(cè)試電泳結(jié)果如圖2所示，由于gyrB基因不僅存在于沙門氏菌中，在志賀桿菌、大腸桿菌等腸科桿菌中都存在該基因，因此樣品中只要有以上任一菌種，都能擴(kuò)增出目的條帶。引物主要是優(yōu)勢(shì)擴(kuò)增沙門氏菌，但是該引物對(duì)于其他沙門氏菌和腸科桿菌的模板DNA還是有微弱的結(jié)合能力，因此即使是沙門氏菌陰性樣品，也能夠擴(kuò)增得到目的條帶。針對(duì)gyrB基因設(shè)計(jì)的1對(duì)引物(測(cè)試了鼠傷寒沙門氏菌和志賀桿菌其中的兩個(gè)樣本)，條帶單一，條帶大小約為456 bp，通過電泳圖可以看到此對(duì)引物對(duì)沙門菌的gyrB基因優(yōu)勢(shì)擴(kuò)增，條帶較亮，對(duì)屬于腸桿菌屬的志賀桿菌的gyrB基因擴(kuò)增能力有很大程度降低，條帶較淡。

表2 gyrB基因通用引物基本信息

M-DL2000 marker 1-鼠傷寒沙門氏菌樣本 2-志賀桿菌樣本 3-無模板陰性對(duì)照?qǐng)D2 gyrB基因通用引物引物測(cè)試電泳圖M-DL2000 marker. 1: Salmonella typhimurium sample; 2: Shigella sample; 3: Negative control.Fig.2 Electrophoregram of the universal primers for GyrB gene

2.2 沙門氏菌建庫(kù)PCR條件優(yōu)化

2.2.1 退火溫度： 第1輪PCR，退火溫度56℃，無擴(kuò)增條帶，未成功；第2輪PCR，退火溫度54℃，無擴(kuò)增條帶，未成功；第3輪PCR，退火溫度52℃，有擴(kuò)增條帶，成功。

M： DL2000 marker模板量為1: 500 ng；2: 200 ng；3: 100 ng；4: 50 ng；5： 25 ng；6: 10 ng；7: 1 ng；8: 100 pg；9： 10 pg；10: 1 pg；11: 100 fg；12: 10 fg。圖3 沙門氏菌模板濃度梯度實(shí)驗(yàn)電泳結(jié)果圖M: DL2000 Marker. Loading quantity of samples. 1: 500 ng; 2: 200 ng; 3: 100 ng; 4: 50 ng; 5： 25 ng; 6: 10 ng; 7: 1 ng; 8: 100 pg; 9: 10 pg; 10: 1 pg; 11: 100 fg; 12: 10 fg.Fig.3 Results of concentration gradient experiment of the Salmonella templates

2.2.2 模板上樣量： 針對(duì)模板上樣量設(shè)計(jì)一個(gè)模板濃度梯度實(shí)驗(yàn)，共分為12個(gè)梯度，依次為500、200、100、50、25、10、1 ng、100、10、1、100、10 fg，實(shí)驗(yàn)結(jié)果如下電泳圖3。通過以上模板濃度梯度實(shí)驗(yàn)，可以發(fā)現(xiàn)模板在10 ng～200 ng范圍內(nèi)，都能擴(kuò)增得到目的片段。為了保證最后的產(chǎn)物有足夠的量，也就是要保證目的條帶足夠亮，選擇較為合適的模板濃度150 ng作為體系最終的模板濃度?？偨Y(jié)沙門氏菌gyrB基因建庫(kù)PCR體系為DNA 150 ng、PremixExTaq12.5 μL，上下游優(yōu)化引物 ( 10 pmol /μL) 各2.5 μL，補(bǔ)充ddH2O至20 μL；PCR最佳反應(yīng)條件為: 98℃ 預(yù)變性 30 s; 98℃ 變性 10 s，52℃退火30s， 72℃延伸 10 min，共 40個(gè)循環(huán); 72℃再延伸 5 min。

2.3 樣本正式建庫(kù)

采用以上摸索得到的建庫(kù)PCR條件將所有樣品進(jìn)行建庫(kù)，建庫(kù)PCR結(jié)果如圖4～6所示。結(jié)果顯示，42個(gè)樣本建庫(kù)成功，均可得到456 bp大小的目的條帶。

2.4 檢測(cè)結(jié)果的準(zhǔn)確性、特異性及重復(fù)性分析

42個(gè)樣本的高通量測(cè)序檢測(cè)結(jié)果如表3所示，各樣本檢測(cè)到的主體菌種與目的菌種基本上一致的。通過測(cè)序的方法在樣品中存在的極微量的菌種也被檢測(cè)到，因此在一些樣品中除了要檢測(cè)的主體菌種外，還存在極微量拷貝數(shù)的其他沙門氏菌菌種。在志賀桿菌中，特別設(shè)置了一個(gè)志賀桿菌的數(shù)據(jù)庫(kù)用于19～23號(hào)樣品的比對(duì)，發(fā)現(xiàn)無法檢測(cè)到該菌種，說明設(shè)計(jì)的引物可以特異擴(kuò)增沙門氏菌。在混合菌種中，該引物能將兩種不同的菌種進(jìn)行區(qū)分。針對(duì)腸炎沙門氏菌設(shè)置的cfu梯度，從檢測(cè)的結(jié)果可以看出，在0～102的cfu范圍內(nèi)，也能檢出，在其他cfu梯度樣品中都能穩(wěn)定檢出，準(zhǔn)確性及靈敏度高，各平行樣本之間的檢測(cè)結(jié)果一致，說明該方法具有較高重復(fù)性。通過高通量測(cè)序技術(shù)能夠?qū)z測(cè)出極微量的沙門氏菌，并能沙門氏菌很好的區(qū)分出來。

M： DL2000 marker 1-17 :1-17號(hào)樣本圖4 1-17號(hào)樣本建庫(kù)PCR電泳圖M： DL2000 marker 1-17 : Sample 1-17Fig.4 PCR electrophoregram of the sequencing library of samples 1-17

M： DL2000 marker. 18-34:18-34號(hào)樣本圖5 18-34號(hào)樣本建庫(kù)PCR電泳圖M： DL2000 marker. 18-34: Samples 18-34Fig.5 PCR electrophoregram of the sequencing library of samples 18-34

表3 沙門氏菌樣本的數(shù)據(jù)分析結(jié)果

注： ST：Salmonellaentericaserovar Typhimurium；鼠傷寒沙門氏菌；SA：Salmonellaentericaserovar Paratyphi A；甲型副傷寒沙門氏菌；SE：Salmonellaentericaserovar Enteritidis；腸炎沙門氏菌；SH： Shigella；志賀菌

M： DL2000 marker. 35-42: 35-42號(hào)樣本圖6 35-42號(hào)樣本建庫(kù)PCR電泳圖M： DL2000 marker. 35-42: Sample 35-42Fig.6 PCR electrophoregram of the sequencing library of samples 35-42

3 討論

人畜共患病是指由同一種病原體引起，流行病學(xué)上相互關(guān)聯(lián)，在人類和動(dòng)物之間自然傳播的疫病。其病原包括病毒、細(xì)菌、支原體、螺旋體、立克次氏體、衣原體、真菌、寄生蟲等。這類疾病的特性在于很不容易消滅。因此這些病原微生物的存在會(huì)對(duì)人類健康、畜牧業(yè)安全生產(chǎn)、畜產(chǎn)品安全和公共衛(wèi)生造成重大危害[11]。沙門氏菌是一類可以引起沙門氏菌病的典型人畜共患病原菌, 并且在實(shí)驗(yàn)動(dòng)物中十分常見，可感染實(shí)驗(yàn)人員并對(duì)實(shí)驗(yàn)結(jié)果造成不可預(yù)估的影響，因此對(duì)于沙門氏菌的檢測(cè)就變得尤為重要。

如今，分子生物學(xué)技術(shù)已越來越多的被應(yīng)用于菌種鑒定，比如現(xiàn)在較為成熟的16S rDNA細(xì)菌菌種鑒定技術(shù)，ITS區(qū)真菌菌種鑒定技術(shù)等等，都體現(xiàn)出了簡(jiǎn)便、快速、高效的特點(diǎn)。菌種鑒定到屬的級(jí)別還是相對(duì)較容易的，現(xiàn)在也有很多技術(shù)可以達(dá)成，但是想要鑒定到種、亞種甚至是菌株的水平，鑒定的困難程度就要提高很多。16S rDNA能很好的將菌種鑒定到屬的層面, 但是由于16S rDNA的進(jìn)化速度非常慢，即太保守，雖對(duì)種系有一定的區(qū)分能力，但是難以區(qū)分某些種系非常接近的菌種及同一菌種的不同菌株。23S rDNA也可用于菌種鑒定，有些學(xué)者認(rèn)為其系統(tǒng)進(jìn)化分析能力甚至比16S rDNA更強(qiáng)，但是由于23S rDNA的分子大小較大，因此它的使用也受到了限制[12，13]。核糖體基因間隔區(qū)(intergenic spacer region, ISR)為rDNA區(qū)域之間高度可變的區(qū)域，它的進(jìn)化速度遠(yuǎn)遠(yuǎn)的高于rDNA區(qū)域，可以用于菌種間的鑒別。ISR區(qū)域主要有16S ～23S rDNA ISR區(qū)域和23S～5S rDNA ISR區(qū)域。ISR區(qū)域堿基變化程度與rDNA區(qū)域的程度相差不多，它優(yōu)于rDNA區(qū)域的原因主要是該區(qū)域存在片段插入或缺失，導(dǎo)致片段大小的變化從而產(chǎn)生了種間和種內(nèi)差異[14]。Morales等還發(fā)現(xiàn)23S～5S rDNA ISR區(qū)域比16S ～23S rDNA ISR區(qū)域具有更好的菌種區(qū)分能力[15]。然而，在某種程度上，ISR區(qū)域的種系鑒別能力也因其高度變異性和高度進(jìn)化率而受到限制。

細(xì)菌中的一些基因也可以作為系統(tǒng)發(fā)育標(biāo)記基因，比如gyrB基因，該基因編碼DNA解旋酶(gyrase)的B亞單位，具有DNA依賴的ATP酶活性，催化ATP的水解。Fukushima應(yīng)用gyrB基因?qū)Υ竽c埃希菌、沙門氏菌和志賀桿菌進(jìn)行系統(tǒng)發(fā)育分析，證實(shí)了gyrB基因優(yōu)于16S rDNA基因[14]。由于gyrB基因的進(jìn)化速度比16S rDNA等其他基因快，因此更適用于親緣關(guān)系較近的種內(nèi)及種間的鑒定和分類。

在本研究中，針對(duì)16S rDNA、3S rDN、16S～23S rDNA IS、23S～5S rDNA IS、gyrB優(yōu)選區(qū)域設(shè)計(jì)通用引物，對(duì)各個(gè)引物進(jìn)行測(cè)試分析, 最終篩選發(fā)現(xiàn)沙門氏菌的菌種分類優(yōu)選區(qū)域?yàn)間yrB基因，gyrB基因引物序列FP 5’-AACCACCGCAATCAGACCTT-3’，RP 5’-AGCCACGAAACCTTCACYA-3’。優(yōu)化擴(kuò)增條件并建庫(kù)，利用Illumina高通量測(cè)序技術(shù)檢測(cè)區(qū)分42個(gè)樣本進(jìn)行通過高通量測(cè)序和序列分析，能夠檢測(cè)出樣本中極微量的沙門氏菌，檢測(cè)方法穩(wěn)定，檢測(cè)限可達(dá)0～102的cfu，檢測(cè)方法穩(wěn)定性好，靈敏度高，可沿用至其他科學(xué)研究用的動(dòng)物(如大鼠、豚鼠、家兔等)病原微生物的檢測(cè)，同時(shí)還可將該實(shí)驗(yàn)體系應(yīng)用于其他種類病原微生物的檢測(cè)。

[1] 王軍, 鄭增忍, 王晶鈺. 動(dòng)物源性食品中沙門氏茵的風(fēng)險(xiǎn)評(píng)估 [J]. 中國(guó)動(dòng)物檢疫, 2007, 24(4): 23-25.

[2] 張嘉寧, 藏富妍, 顧為望. 應(yīng)用PCR法檢測(cè)小鼠沙門氏菌 [J]. 中國(guó)比較醫(yī)學(xué)雜志, 1999, 9(1): 45-47.

[3] McGuinness S, McCabe E, O’Regan E, et al. Development and validation of a rapid real-time PCR based method for the specific detection of Salmonella on fresh meat [J]. Meat Sci, 2009, 83(3): 555-562.

[4] Liu FP, Tong HM.GyrAandParCgene mutation of clinically isolated fluoroquinolone-resistant strain of salmonella [J]. J Northeast Agr Univ, 2006,13(1): 47-50.

[5] Yan SS, Pendrak ML, Abela-Ridder B, et al. An overview of Salmonella typing — public health perspectives [J]. Clin Appl Immunol Rev, 2003, 4(3): 189-204.

[6] 范媛媛, 鄭冰, 王樹祥, 等. 一種新型培養(yǎng)基在沙門氏菌檢測(cè)中的效果評(píng)價(jià) [J]. 食品工業(yè)科技, 2008, 8: 415-417.

[7] 楊保偉, 盛敏, 席美麗, 等. 食源性沙門氏菌耐藥性檢測(cè)及相關(guān)質(zhì)粒 [J]. 微生物學(xué)報(bào), 2008, 48(8): 1006-1012.

[8] He LL, Sok D, Azadnia P, et al. Toward a more accurate view of human B-cell repertoire by next-generation sequencing, unbiased repertoire capture and single-molecule barcoding [J].Sci Rep, 2014, 4: 67-78.

[9] Weinstein JA, Jiang N, White III RA, et al. High-throughput sequencing of the Zebrafish antibody repertoire [J]. Science, 2009, 324(5928): 807-810.

[10] DeKosky BJ, Ippolito GC, Deschner RP, et al. High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire [J]. Nat Biotech, 2013, 31(2):166-169.

[11] 王章云, 滕煥昭, 李柏桂, 等. 腸炎沙門氏菌引起的食物中毒細(xì)菌學(xué)調(diào)查 [J]. 中國(guó)人獸共患病雜志, 1999, 15(3): 115.

[12] Henrik C, Steen N, John EO. Phylogenetic relationships of Salmonella based on rRNA sequences [J]. Int J Systematic Bacteriol, 1998, 48: 605-610.

[13] Sara PL, Francisco R, Ruiting L. Variation of the ribosomal operon 16S-23S gene spacer region in representatives ofSalmonellaentericasubspecies [J]. J Bacteriol, 1998, 180(8) :2144-2151.

[14] Cesar AM, Richard G, Jean G. Linkage of avian and reproductive tract tropism with sequence divergence adjacent to the 5S ribosomal subunitrrfHofSalmonellaenterica[J]. FEMS Microbiol Lett, 2006, 264: 48-58.

[15] Fukushima M, Kakinuma R. Phylogenetic analysis of Salmonella, Shigella, and Escherichia coli strains on the basis of thegyrBgene sequence [J]. J Clin Microbiol, 2002, 40(8): 2779-2785.

Establishment and evaluation of a high throughput sequencing technology for detection ofSalmonellain laboratory animals

HU Yi-xiang, ZHANG Huan-huan, YU Chen-huan, JIN Xiao-yin, YING Hua-zhong

(Zhejiang Key Laboratory of Experimental Animal and Safety Evaluation, Zhejiang Academy of Medical Sciences, Hangzhou 310013, China)

Objective To establish a detection method ofSalmonellain laboratory animals based on a high-throughput sequencing technology, and to apply it in detection ofSalmonellain laboratory animals. Methods DNA samples were extracted from mouse feces. Universal primers for 16S rDNA, 23S rDNA, 16S-23S rDNA, 23S-5S rDNA region, gyrB preferred area were designed, respectively. Each primer was tested and analyzed to determine the best amplification conditions and build a database. Forty-two samples ofSalmonellawere assayed by Illumina high-throughput sequencing technology and evaluated the specificity and stability of this method. Results The species preferred region ofSalmonellawas gyrB gene region. The primers for gyrB gene were FP5’-AACCACCGCAATCAGACCTT3‘ and FP5’-AGCCACGAAACCTTCACYA-3’. The primers were optimized and determined, through a high-throughput sequencing, and the sequence analysis detected very small amount ofSalmonellain the 42 samples, indicating that this detection method is stable, highly sensitive, and the limit of detection reached to 0-102CFU. Conclusions We have established a complete detection system for detection ofSalmonellain laboratory animals based on a high-throughput sequencing technology, This system can detect trace amounts ofSalmonellain laboratory animals, and this detection method is stable and highly sensitive, which can be also used in detection of other kinds of pathogenic microorganism in laboratory animals.

Salmonella; Laboratory animal; High-throughput sequencing

浙江省衛(wèi)生高層次人才(2012F30026)，浙江省科技廳院所專項(xiàng)(2016F30001)。

胡毅翔(1989-)，男，碩士研究生，研究方向?yàn)橹兴幩幚韺W(xué)。Email: aahyxaa@163.com。

應(yīng)華忠(1968-)，男，研究員，從事動(dòng)物模型與比較醫(yī)學(xué)研究。 Email: YHZ0101@126.com。

R-33

A

1671-7856(2016)10-0072-07

10.3969.j.issn.1671-7856.2016.10.014

2016-03-08