肖 穎，武雅婷，趙婉妤,2，崔學(xué)文，黎 明，許 欣

沙門(mén)菌(Salmonella)是革蘭陰性人獸共患病原菌，與人類食源性疾病密切相關(guān)。近年來(lái)，由沙門(mén)菌引起的食源性疾病暴發(fā)事件與日俱增。2010-2018年，四川省食源性疾病暴發(fā)事件中，44.01%由沙門(mén)菌感染引起，其中成都市的發(fā)病人數(shù)居四川省首位[1]。了解成都市沙門(mén)菌的基因組特征對(duì)于監(jiān)測(cè)與防控沙門(mén)菌引發(fā)的食源性疾病具有重要意義，但目前未見(jiàn)成都市人源沙門(mén)菌全基因組分析及報(bào)道。

沙門(mén)菌的血清型種類復(fù)雜，截止目前，世界范圍內(nèi)已經(jīng)發(fā)現(xiàn)了2 659種不同的沙門(mén)菌血清型[2]。然而傳統(tǒng)的血清學(xué)分型方法工作量大、耗時(shí)長(zhǎng)，全基因組測(cè)序(whole genome sequencing，WGS)作為一種快速且分辨率更高的方法，可通過(guò)測(cè)序數(shù)據(jù)預(yù)測(cè)血清型、分析耐藥基因與毒力因子等，國(guó)外已有學(xué)者將其運(yùn)用在沙門(mén)菌的監(jiān)測(cè)與疫情調(diào)查中[3-5]。研究表明[6]，WGS與常規(guī)血清學(xué)分型方法的符合率為97.6%，具有廣闊的應(yīng)用前景。

因此，本研究對(duì)成都市35株人源沙門(mén)菌分離株進(jìn)行全基因組測(cè)序，預(yù)測(cè)其血清型、耐藥基因及可移動(dòng)遺傳元件(mobile genetic elements，MGE)，注釋并分析其毒力因子及分布，以探究成都市人源沙門(mén)菌耐藥與毒力分子特征，為深入研究沙門(mén)菌致病機(jī)理、監(jiān)測(cè)與預(yù)防沙門(mén)菌感染提供參考依據(jù)。

1 材料與方法

1.1 菌株來(lái)源 35株臨床沙門(mén)菌，分離自2020年5-6月的腹瀉患者糞便，菌株來(lái)自成都市青羊區(qū)(14株)、郫都區(qū)(7株)、成華區(qū)(4株)、金牛區(qū)(3株)、武侯區(qū)(2株)、天府新區(qū)(2株)、新都區(qū)(2株)、錦江區(qū)(1株)8個(gè)行政區(qū)，菌株分離與收集由成都市疾病預(yù)防控制中心完成。利用VITEK 2 COMPACT 全自動(dòng)微生物分析鑒定儀(法國(guó)BioMerienx公司)對(duì)菌株進(jìn)行沙門(mén)菌鑒定。

1.2 全基因組測(cè)序 收集到的沙門(mén)菌分離株送往中國(guó)科學(xué)院微生物研究所進(jìn)行全基因組測(cè)序，流程如下：提取菌株DNA，使用Nanodrop檢測(cè)DNA基因組濃度與質(zhì)量，進(jìn)行文庫(kù)構(gòu)建；使用IlluminaXten測(cè)序系統(tǒng)進(jìn)行雙端150 bp測(cè)序；利用Fastqc、Trimmomatic對(duì)下機(jī)原始數(shù)據(jù)進(jìn)行質(zhì)量檢測(cè)、統(tǒng)計(jì)；使用SPAdes的Unicycler短序列組裝軟件進(jìn)行組裝，并根據(jù)reads中paired-end、overlap關(guān)系對(duì)組裝結(jié)果進(jìn)行局部?jī)?yōu)化；使用GapClose 1.1.2軟件對(duì)組裝結(jié)果進(jìn)行修補(bǔ)，除去冗余序列獲得最終組裝結(jié)果。

1.3 血清學(xué)預(yù)測(cè)與MLST分型 基于ZOOM TECH 微生物數(shù)據(jù)分析平臺(tái)(http://39.107.235.93:8080/#/login)，使用SeqSero根據(jù)特定抗原的多態(tài)性完成血清型預(yù)測(cè)，使用PubMLST對(duì)分離株進(jìn)行MLST分型。

1.4 基因預(yù)測(cè)、注釋與分析

1.4.1 基因組組分分析 使用 GeneMarkS 分析工具(http://topaz.gatech.edu/GeneMark/genemarks.cgi)預(yù)測(cè)沙門(mén)菌分離株編碼基因。

1.4.2 耐藥基因預(yù)測(cè)、注釋與MGE分析 基于CARD(Comprehensive Antibiotic Research Database)數(shù)據(jù)庫(kù)對(duì)沙門(mén)菌分離株耐藥基因進(jìn)行預(yù)測(cè)及注釋。利用Center for Genomic Epidemiology平臺(tái)的Mobile Element Finder分析工具(https://cge.cbs.dtu.dk/services/MobileElementFinder)預(yù)測(cè)菌株的可移動(dòng)遺傳元件，包括質(zhì)粒與插入片段。

1.4.3 毒力因子預(yù)測(cè)、注釋與分析 基于VFDB(Virulence Factor Database)數(shù)據(jù)庫(kù)對(duì)沙門(mén)菌分離株毒力因子進(jìn)行預(yù)測(cè)及注釋。使用R 4.1.0軟件構(gòu)建毒力因子數(shù)目分析熱圖，對(duì)每個(gè)基因組的毒力因子數(shù)目進(jìn)行Scale參數(shù)標(biāo)準(zhǔn)化。

2 結(jié) 果

2.1 血清學(xué)預(yù)測(cè)與MLST分型 35株人源沙門(mén)菌分離株共預(yù)測(cè)出5種血清型(表1)，其中I 4,[5],12:i:-血清型沙門(mén)菌最多(15株，42.9%)，其次為鼠傷寒血清型(Typhimurium)沙門(mén)菌(12株，34.3%)。

表1 35株沙門(mén)菌分離株血清型預(yù)測(cè)結(jié)果Tab.1 Serotype prediction results for 35 Salmonella isolates

MLST分析共鑒定出5種ST型。如圖1所示，血清型預(yù)測(cè)與ST型結(jié)果對(duì)應(yīng)良好，I 4,[5],12:i:-血清型預(yù)測(cè)菌株除2株未在數(shù)據(jù)庫(kù)中查得所對(duì)應(yīng)的ST型，其余13株均為ST34型，鼠傷寒血清型(Typhimurium)預(yù)測(cè)菌株均為ST19型，腸炎血清型(Enteritidis)預(yù)測(cè)菌株為ST11型，德比血清型(Derby)預(yù)測(cè)菌株為ST40型，火雞血清型(Meleagridis)預(yù)測(cè)菌株為ST463型。

“-”表示未在數(shù)據(jù)庫(kù)中查得所對(duì)應(yīng)的ST型。圖1 35株沙門(mén)菌分離株血清型預(yù)測(cè)與MLST分型結(jié)果對(duì)比Fig.1 Comparison between serotype prediction and MLST typing of 35 Salmonella isolates

2.2 基因組組分分析 35株沙門(mén)菌分離株基于測(cè)序數(shù)據(jù)的基因組組分結(jié)果見(jiàn)表2，基因組范圍為4.7～5.2 Mb；GC含量為52%左右；編碼基因數(shù)量4 600～5 200 個(gè)，平均長(zhǎng)度約 900 bp，總長(zhǎng)度約占基因組的 88%。I 4,[5],12:i:-血清型沙門(mén)菌的基因組大小和編碼基因數(shù)目高于其他血清型菌株。

表2 35株沙門(mén)菌分離株基因組組分預(yù)測(cè)Tab.2 Genome component prediction for 35 Salmonella isolates

2.3 耐藥基因與MGE分析 35株沙門(mén)菌分離株基于CARD數(shù)據(jù)庫(kù)共預(yù)測(cè)到10類42種不同的耐藥基因(圖2)，以氨基糖苷類耐藥基因(42.68%，102/239)最多，其中aac(6′)-Iaa攜帶率為100%(35/35)。不同血清型沙門(mén)菌的耐藥基因攜帶情況存在差異。I 4,[5],12:i:-沙門(mén)菌攜帶耐藥基因數(shù)目更多，且硫胺類sul2、四環(huán)素類tet(B)攜帶率為100%(15/15)。

35株沙門(mén)菌分離株質(zhì)粒攜帶率為97.14%(34/35)，單個(gè)菌株最多同時(shí)攜帶5個(gè)質(zhì)粒。共預(yù)測(cè)出17種不同的質(zhì)粒，包括IncQ1(40.00%，14/35)、IncFII(S)與IncFIB(S)(37.14%，13/35)、Col440I(22.86%，8/35)等。其中IncQ1質(zhì)粒只在I 4,[5],12:i:-沙門(mén)菌中攜帶，腸炎沙門(mén)菌和部分鼠傷寒(8株)均同時(shí)攜帶IncFII(S)和IncFIB(S)質(zhì)粒。MGE分析結(jié)果如表3，部分質(zhì)粒與耐藥基因關(guān)系密切，且存在插入序列。I 4,[5],12:i:-沙門(mén)菌攜帶更多的質(zhì)粒與插入序列，其攜帶的IncQ1質(zhì)粒(93.33%，14/15)與多種耐藥基因相關(guān)。

表3 35株沙門(mén)菌分離株耐藥基因及相關(guān)可移動(dòng)遺傳元件的攜帶情況Tab.3 Drug resistance genes and mobile genetic elements carried by 35 Salmonella isolates

2.4 毒力因子分析 基于VFDB數(shù)據(jù)庫(kù)對(duì)35株沙門(mén)菌分離株的毒力因子注釋(Identity≥80)，共檢出120種3 998個(gè)毒力基因，菌株毒力因子數(shù)目為103～123個(gè)(表4)。分泌系統(tǒng)相關(guān)的毒力因子數(shù)目最多，為76～79個(gè)。不同血清型的菌株之間毒力因子數(shù)目不同，I 4,[5],12:i:-、鼠傷寒和腸炎沙門(mén)菌的毒力因子數(shù)目大于其他血清型沙門(mén)菌。Spv、Pef、Rck毒力因子均只在13株沙門(mén)菌中存在，即全部的腸炎沙門(mén)菌和8株鼠傷寒沙門(mén)菌，且僅該8株鼠傷寒沙門(mén)菌保留GrvA毒力因子。僅德比和火雞沙門(mén)菌缺失SodCI毒力因子。Fim、Agf、Csg、Mgt、Inv、misL、avrA、Org、Ssa、ompA、mig-14、Ste、Sop等毒力因子在所有菌株中數(shù)目相等。

表4 35株沙門(mén)菌分離株基因組VFDB毒力因子注釋個(gè)數(shù)Tab.4 Genome annotation numbers of virulence factors in VFDB for 35 Salmonella isolates

毒力因子分布熱圖(圖3)所示，35株沙門(mén)菌基于毒力因子分布分成3簇，其中德比沙門(mén)菌單獨(dú)成簇，腸炎沙門(mén)菌和部分鼠傷寒沙門(mén)菌(8株)共同成簇，剩余鼠傷寒沙門(mén)菌(4株)與火雞、I 4,[5],12:i:-沙門(mén)菌共同成簇。鼠傷寒沙門(mén)菌均具有更高的鞭毛毒力因子，且呈現(xiàn)出兩種毒力因子分布，其中，部分鼠傷寒沙門(mén)菌(8株)呈現(xiàn)出更高的毒力質(zhì)粒、質(zhì)粒編碼菌毛、血清抗性3種毒力因子分布。腸炎沙門(mén)菌的毒力質(zhì)粒、質(zhì)粒編碼菌毛、血清抗性毒力因子分布與部分鼠傷寒沙門(mén)菌(8株)類似，但鞭毛毒力因子分布更少。德比沙門(mén)菌纖維粘附素相關(guān)毒力因子Lpf(Long polar fimbriae)分布低于其他血清型。火雞、I 4,[5],12:i:-沙門(mén)菌毒力因子分布較為相似。

圖3 35株沙門(mén)菌分離株毒力因子分布熱圖Fig.3 Virulence factor annotation heatmap for 35 Salmonella isolates

3 討 論

I 4,[5],12:i:-沙門(mén)菌是過(guò)去20年里世界范圍內(nèi)迅速出現(xiàn)的新亞型，被認(rèn)為是鼠傷寒沙門(mén)菌第Ⅱ相鞭毛蛋白不表達(dá)的變種[7]。國(guó)內(nèi)已有多地報(bào)道[8-10]I 4,[5],12:i:-為臨床沙門(mén)菌分離株的優(yōu)勢(shì)血清型，且在多種動(dòng)物中發(fā)現(xiàn)該血清型沙門(mén)菌的流行率逐年增高[11-12]，提示I 4,[5],12:i:-已逐漸成為人與動(dòng)物沙門(mén)菌感染的主要血清型之一。本研究I 4,[5],12:i:-沙門(mén)菌主要為ST34型，這與國(guó)內(nèi)外學(xué)者的研究結(jié)果一致[3,13]。

本研究I 4,[5],12:i:-沙門(mén)菌耐藥基因以氨基糖苷類aac(6′)-Iaa、硫胺類sul2、四環(huán)素類tet(B)為主，這與目前國(guó)內(nèi)研究結(jié)果一致[13]。本研究幾乎全部I 4,[5],12:i:-沙門(mén)菌攜帶lncQ1質(zhì)粒，已有研究發(fā)現(xiàn)[14]lncQ1質(zhì)粒是sul2、aph(6)-Id耐藥基因的載體，可在不同來(lái)源的多種細(xì)菌中穩(wěn)定維持和轉(zhuǎn)移，這提示I 4,[5],12:i:-血清型可能存在該質(zhì)?；蚪閷?dǎo)的獲得性耐藥。此外，本研究MGE分析發(fā)現(xiàn)aadA1、cmlA1等耐藥基因與lncQ1質(zhì)粒相關(guān)，且存在插入片段，其具體關(guān)系目前尚未得到充分研究與闡述，有待于進(jìn)一步探明。

Spv毒力質(zhì)粒(SalmonellaPlasmid Virulence)是沙門(mén)菌重要的毒力因子，其基因包括spvR、spvA、spvB、spvC和spvD[15]，其中，spvB、spvC和spvR是沙門(mén)菌的毒力表型所必需的[16]。spvB具有細(xì)胞毒性作用，與沙門(mén)菌感染期間的宿主細(xì)胞晚期凋亡有關(guān)[17]；spvC具有蘇氨酸裂解酶活性，可以抑制宿主細(xì)胞的MAP激酶信號(hào)[18-19]；spvR對(duì)促進(jìn)下游spvA、spvB、spvC和spvD基因的轉(zhuǎn)錄具有核心作用[16,20]。本研究中，spvB、spvC、spvR均只在腸炎血清型菌株和部分鼠傷寒血清型菌株中發(fā)現(xiàn)(37.1%，13/35)。這與目前研究結(jié)論一致[19,21]，即只有腸炎和鼠傷寒沙門(mén)菌攜帶Spv毒力質(zhì)粒。研究表明Spv操縱子通常位于IncFII質(zhì)?；騃ncFII-InFIB雜交質(zhì)粒上[22]，本研究同樣發(fā)現(xiàn)僅預(yù)測(cè)出Spv毒力基因的沙門(mén)菌攜帶IncFII(S)和IncFIB(S)質(zhì)粒。此外，本研究發(fā)現(xiàn)，pefA、pefB、pefC、pefD、rck毒力因子也只在上述腸炎和鼠傷寒沙門(mén)菌中保留。Pef為質(zhì)粒編碼菌毛(Plasmid-encoded fimbriae)相關(guān)基因，可調(diào)節(jié)沙門(mén)菌對(duì)腸上皮細(xì)胞的粘附[23]；Rck為血清抗性相關(guān)基因，與沙門(mén)菌抵抗補(bǔ)體殺傷作用有關(guān)[24]。一項(xiàng)體外實(shí)驗(yàn)顯示沙門(mén)菌毒力與Spv、Pef、Rck毒力基因的攜帶呈正相關(guān)[21]。Spv與Pef、Rck毒力因子可能存在協(xié)同致病機(jī)制，值得在轉(zhuǎn)錄組和蛋白組進(jìn)一步研究。

SodCI毒力因子是Cu-Zn 超氧化物歧化酶，被認(rèn)為有助于沙門(mén)菌在巨噬細(xì)胞中的生存[25]，且是鼠傷寒沙門(mén)菌發(fā)揮毒力所必需的[26]。本研究中，SodCI毒力因子在鼠傷寒、腸炎、I 4,[5],12:i:-沙門(mén)菌中保留，德比和火雞沙門(mén)菌缺失，這與Kuijpers等[21]的結(jié)果一致，SodCI毒力因子可能與這些沙門(mén)菌的較強(qiáng)致病力有關(guān)。

本研究系首次采用WGS技術(shù)對(duì)成都市人源沙門(mén)菌進(jìn)行分型、耐藥基因與毒力因子注釋及分析，描述其基因組特征，豐富了該地沙門(mén)菌監(jiān)測(cè)數(shù)據(jù)庫(kù)，在基因水平為研究沙門(mén)菌耐藥、毒力機(jī)制提供了資料，其致病機(jī)理還需在基因組、轉(zhuǎn)錄組與蛋白組不斷探究及驗(yàn)證。

利益沖突：無(wú)

引用本文格式：肖穎，武雅婷，趙婉妤，等. 成都市人源沙門(mén)菌基因組特征分析[J].中國(guó)人獸共患病學(xué)報(bào)，2022，38(5)：433-440，446. DOI:10.3969/j.issn.1002-2694.2022.00.059