全基因組測序技術在豬源沙門氏菌血清分型和耐藥性分析中的應用

張 璐，沈青春，徐錦濤，趙 琪，李 霆，程 敏，張純萍*

(1.中國獸醫(yī)藥品監(jiān)察所，北京100081；2.皇譽寵物食品有限公司，上海 200000)

沙門氏菌是自然界中最常見的人畜共患病原菌，能夠引起仔豬的副傷寒、敗血癥等疾病，而且豬感染后常成為隱性帶菌者，不定期排毒，從而感染其他畜禽。更為重要的是，豬所攜帶的沙門氏菌也是一種重要的食源性病原菌，在我國沙門氏菌食物中毒事件中，90%以上的是由肉類食品引起的，其中豬肉及其產(chǎn)品是主要污染源之一[1]。血清分型是實現(xiàn)病原體溯源分析的主要方法，但沙門氏菌的血清型眾多，常規(guī)血清分型方法較為繁瑣、耗時長，且對于出現(xiàn)交叉凝集或凝集現(xiàn)象不明顯的菌株，結(jié)果較難判斷[2]。沙門氏菌的預防和治療主要依賴于抗菌藥物，但抗菌藥物的頻繁使用，產(chǎn)生了大量的耐藥菌株。目前細菌耐藥性已經(jīng)成為全球共同關注的問題，尤其是在畜牧養(yǎng)殖行業(yè)，嚴重威脅著動物的健康，造成了巨大的經(jīng)濟損失[3]。隨著生物信息技術的發(fā)展，測序成本的大幅降低和多種公共數(shù)據(jù)庫的創(chuàng)建簡便了沙門氏菌血清分型方法和耐藥性的檢測[4]。國外研究表明，通過全基因組測序，可以簡單、快速和高效的檢測沙門氏菌的血清型和耐藥基因，且具有較好的敏感性和特異性[5-7]，但目前國內(nèi)相關文章較少。因此，本文選取了我國2011-2015年不同地區(qū)分離的60株豬沙門氏菌，利用全基因組測序技術分析沙門氏菌的血清型和耐藥性，并與常規(guī)檢驗結(jié)果進行了比較，從而分析全基因組測序技術(Whole genome sequencing，WGS)在沙門氏菌血清分型和耐藥性分析方面的應用能力。

1 材料和方法

1.1 菌株來源 豬源沙門氏菌共60株，分離時間為2011-2015年，其中上海14株、四川10株、廣東36株。質(zhì)控菌株為大腸桿菌ATCC 25922和ATCC 35218，均由中國獸醫(yī)藥品監(jiān)察所安全評價室保存。

1.2 主要試劑和儀器 沙門氏菌診斷血清購自丹麥國家血清研究院，革蘭氏陰性細菌鑒定板由本實驗室設計、上海星佰生物技術有限公司生產(chǎn)，細菌基因組DNA提取試劑盒購自天根生化科技(北京)有限公司，沙門氏菌全基因組測序分析由奧維森基因科技有限公司完成。

1.3 血清凝集試驗 參照沙門氏菌檢驗標準GBT4789.4-2016[8]，采用傳統(tǒng)玻片凝集法對沙門氏菌進行血清型鑒定。O抗原的鑒定需要先在潔凈的玻片上滴加1滴診斷血清，然后用潔凈的槍頭挑取營養(yǎng)瓊脂上培養(yǎng)的單菌落與O抗原多價血清充分混勻，晃動30～60 s，出現(xiàn)凝聚現(xiàn)象為陽性反應。如果是由于Vi抗原的存在而阻止了O凝集反應時，可挑取菌落在1 mL生理鹽水中做成濃菌液，在酒精燈火焰上煮沸后，若凝集則具有Vi抗原，同時設置無菌生理鹽水為對照，在生理鹽水中自凝者為粗糙形菌株，不能分型。H抗原的鑒定則需挑取swarm半固體瓊脂培養(yǎng)的周邊菌落與H抗原多價血清充分混勻，觀察凝集現(xiàn)象。最后根據(jù)O相與H相的凝集結(jié)果，查閱WHO的White-Kauffmann-Le Minor抗原表[9]確定血清型。

1.4 藥敏試驗 采用微量肉湯稀釋法對不同年代沙門氏菌進行耐藥性(磺胺異噁唑、氨芐西林、阿莫西林/克拉維酸、四環(huán)素、頭孢噻呋、阿奇霉素、復方新諾明、粘菌素、恩諾沙星、氟苯尼考、美羅培南)特征分析：挑取2～3個單菌落置5 mL滅菌生理鹽水中，用0.5麥氏比濁管進行比濁，調(diào)制菌液濃度為1.5×108cfu/mL左右；取菌液60 μL加入12 mL的藥敏接種培養(yǎng)液中，進行混勻稀釋，每孔100 μL，陰性對照孔單獨加入100 μL藥敏接種培養(yǎng)液，37 ℃恒溫培養(yǎng)箱培養(yǎng)18～20 h。在陰性對照孔清澈，陽性對照孔渾濁以及質(zhì)控菌株的最低抑菌濃度符合規(guī)定范圍的前提下，根據(jù)美國臨床實驗室標準化委員會標準[10]判斷被檢菌株是否耐藥。

1.5 WGS測序分析

1.5.1 測序 根據(jù)細菌基因組DNA提取劑盒說明書提取沙門氏菌DNA，DNA濃度高于100 ng/μL后，送公司經(jīng)Illumina測序平臺測序，然后用Trimmomatic(v0.36)[11]軟件進行數(shù)據(jù)質(zhì)量控制，使用SPAdes(v3.13.0)[12]軟件參考沙門氏菌P125109(ACCESSION:AM933172、基因組大小為4685848 bp)對測序數(shù)據(jù)進行組裝。

1.5.2 生物信息學分析 根據(jù)基因組流行病學中心網(wǎng)站中SeqSero2血清分型數(shù)據(jù)庫和ResFinder4抗性基因數(shù)據(jù)庫中的基因組信息建庫，采用 CentOS 7.8系統(tǒng)服務器創(chuàng)建相應的血清分型數(shù)據(jù)庫和耐藥基因數(shù)據(jù)庫，分析60株沙門氏菌的WGS 組裝數(shù)據(jù)，分別得到血清分型和耐藥基因攜帶結(jié)果。將WGS分型結(jié)果與常規(guī)血清分型結(jié)果進行比較，對結(jié)果有差異的進行再次常規(guī)血清分型驗證。對 WGS檢測的耐藥基因，使用 identity>90.0%，coverage>90.0%的閾值運行命令，進行基因型與耐藥表型進行對比分析[13]。

2 結(jié)果及分析

2.1 全基因組測序、組裝結(jié)果 經(jīng)Illumina測序平臺測序和Trimmomatic軟件質(zhì)控后，共得到約77.59 Gb的原始數(shù)據(jù)，組裝后的60株沙門氏菌的基因組大小在4648280～ 13443523 bp之間，GC含量為45.32%～52.24%，具體結(jié)果見表1。

表1 60株豬源沙門氏菌全基因組序列基本特征Table 1 The basic characteristics of whole genome sequence of 60 Salmonella strains from pigs

2.2 血清凝集分型與全基因組分型結(jié)果比較 60株沙門氏菌通過血清凝集試驗共鑒定出8種血清型，主要血清型為鼠傷寒和德爾卑；WGS分型共鑒定出9種血清型，優(yōu)勢血清型與常規(guī)血清分型結(jié)果一致；兩種方法共有54株分型結(jié)果一致，符合率為90.0%，具體分型結(jié)果見表2。

表2 常規(guī)血清分型與全基因組分型結(jié)果Tab 2 Results of routine serotyping and whole genome typing

2.3 藥敏試驗結(jié)果 沙門氏菌對11種抗菌藥物的耐藥性情況詳見表3。結(jié)果表明，豬源沙門氏菌對不同抗菌藥物的耐藥情況差異較大，其中對美羅培南、頭孢噻呋、阿奇霉素和粘菌素均為敏感，而對四環(huán)素、氨芐西林等7種抗菌藥物均產(chǎn)生了不同程度的耐藥性，尤其是對四環(huán)素(73.3%)的耐藥水平明顯高于其他抗菌藥物。

表3 60株沙門氏菌對11種抗菌藥物的耐藥情況Tab 3 Antimicrobial resistance of 60 Salmonella strains to 11 antimicrobial agents

2.4 全基因組測序預測沙門氏菌的耐藥基因型 基于60株沙門氏菌的WGS測序結(jié)果，使用抗性基因數(shù)據(jù)庫ResFinder4分析獲得的耐藥基因包括β-內(nèi)酰胺類(blaTEM-1、blaOXA-1)、四環(huán)素類(tetA和tetB)、磺胺類(sul1、sul2和sul3)、喹諾酮類(oqxA、oqxB和qnrS2)和酰胺醇類(floR)。

通過比較，WGS預測的沙門氏菌耐藥性與藥敏試驗結(jié)果的符合率較高，其中對恩諾沙星、美羅培南、頭孢噻呋和阿奇霉素預測的結(jié)果完全一致；對氨芐西林、磺胺異噁唑、四環(huán)素預測的結(jié)果符合率均高于90.0%，對復方新諾明、氟苯尼考、阿莫西林/克拉維酸預測符合率也大于75%，具體結(jié)果詳見表4。

表4 沙門氏菌藥敏試驗結(jié)果與全基因組測序結(jié)果的比較Tab 4 Comparison of antimicrobial sensitivity test results and whole genome sequencing results of Salmonella

3 討論與結(jié)論

隨著全基因組測序技術和生物信息分析方法的不斷進步，測序時間與成本的不斷降低，使得公共數(shù)據(jù)庫中存儲了大量的病原菌基因組序列，極大提升了細菌基因組測序數(shù)據(jù)的準確性與可用性。與此同時，簡便易用的在線數(shù)據(jù)分析與可視化工具也在不斷問世，如血清分型數(shù)據(jù)庫SeqSero2和耐藥基因數(shù)據(jù)庫ResFinder4，二者均可通過CGE網(wǎng)站直接上傳WGS組裝序列，即可預測相應的血清型和耐藥基因，簡便了操作流程[14-15]。

本研究利用SeqSero2分析獲得的血清型結(jié)果與常規(guī)分型結(jié)果的符合率為90.0%，這與他人的研究結(jié)果基本一致。如BANERJI等[16]的研究中，WGS預測520株沙門氏菌血清型結(jié)果與原始分型結(jié)果的符合率為98.0%。DIEP等[17]研究的98株沙門氏菌的血清型與常規(guī)分型結(jié)果的符合率為98.0%。在兩種分型結(jié)果不一致的血清型中，大部分是需要特殊培養(yǎng)基培養(yǎng)和多種血清鑒定的罕見血清型，如奧蘭、奧茲馬斯等。這可能是因為常規(guī)血清分型結(jié)果的判讀會因個人判讀標準不同而產(chǎn)生差異，亦或者是由于H相抗原基因中的等位基因在相變機制中發(fā)生突變而不表達而產(chǎn)生不一致的結(jié)果[18]。

為評估WGS在沙門氏菌抗菌藥物耐藥性中的鑒定能力，應用ResFinder4抗性基因數(shù)據(jù)庫分析沙門氏菌所攜帶的耐藥基因，從耐藥基因的攜帶推斷抗菌藥物的耐藥性來看，應用WGS預測的恩諾沙星、美羅培南、頭孢噻呋和阿奇霉素的耐藥性與藥敏試驗結(jié)果完全一致，對氨芐西林、磺胺異噁唑、四環(huán)素預測的結(jié)果符合率也均大于90.0%。ZANKARI等[19]通過WGS預測的49株鼠傷寒沙門氏菌基因型與表型藥敏試驗的結(jié)果完全一致。ZHU等[20]對189株沙門氏菌的耐藥性檢測中，WGS預測的磺胺異噁唑、氨芐西林和四環(huán)素的耐藥性與其耐藥表型的符合率分別為97.8%、94.6%和85.7%。NEUERT等[20]研究的3415株沙門氏菌對15種抗菌藥物的表型耐藥性與基因型耐藥性結(jié)果的符合率為97.8%。目前WGS 只能檢測到質(zhì)粒介導的耐藥基因，而對于染色體突變產(chǎn)生的耐藥基因無法檢出。在對阿莫西林/克拉維酸耐藥基因型檢測中，可能由于基因組測序過程中部分區(qū)域的覆蓋率較低，阻止了突變位點的檢測，亦或者出現(xiàn)了新的耐藥基因突變，從而出現(xiàn)了表型與基因型的符合率相對較低的結(jié)果[21]。

基于全基因組數(shù)據(jù)分析的沙門氏菌血清分型和耐藥性方法具有操作簡便、效率高的優(yōu)勢，可以在常規(guī)血清凝集試驗難以開展或沙門氏菌鞭毛基因不表達的情況下，替代血清凝集試驗，也可快速識別菌株攜帶的所有耐藥基因。為解決沙門氏菌耐藥性問題提供重要參考，在沙門氏菌流行病學的研究中具有較高的應用價值。