翁幸鐾，糜祖煌，王春新，朱健銘

?

多耐藥大腸埃希菌NB8株全基因組耐藥基因檢測(cè)

翁幸鐾1，糜祖煌2，王春新3，朱健銘4

目的 檢測(cè)1株多耐藥大腸埃希菌NB8株的遺傳學(xué)背景。方法 病原分離自2012年4月寧波市第一醫(yī)院住院患者尿液，作16S rDNA和gyrA基因測(cè)序、BLASTn比對(duì)確認(rèn)為大腸埃希菌，采用Illumina HiSeq與Ion Torrent PGM兩種大規(guī)模并行測(cè)序儀進(jìn)行全基因組分析,再行人工測(cè)序。最后NB8株和9株已完成全基因組測(cè)序的大腸埃希菌耐藥基因的比較基因組學(xué)研究。結(jié)果 多耐藥大腸埃希菌NB8株全基因組測(cè)序得到一條推定的染色體序列,長(zhǎng)4 550 369 bp (內(nèi)含14個(gè)缺口)，得到一條推定的質(zhì)粒序列,長(zhǎng)635 377 bp(內(nèi)含33個(gè)缺口)。從NB8株全基因組數(shù)據(jù)中挖掘出了β-內(nèi)酰胺類(lèi)、氨基糖苷類(lèi)、喹諾酮類(lèi)、大環(huán)內(nèi)酯類(lèi)、磺胺類(lèi)、四環(huán)素類(lèi)、多粘菌素的耐藥基因，并挖掘出接合性質(zhì)粒、轉(zhuǎn)座子、插入序列、整合子的遺傳標(biāo)記，以及5大類(lèi)外排泵基因。結(jié)論 多耐藥大腸埃希菌NB8株遺傳學(xué)背景明確，主動(dòng)外排機(jī)制增強(qiáng)也會(huì)導(dǎo)致或者增強(qiáng)NB8株的耐藥性。質(zhì)粒、轉(zhuǎn)座子和整合子等可移動(dòng)遺傳元件可以使細(xì)菌的耐藥性得以快速傳播，使受體菌表現(xiàn)為多重耐藥。

大腸埃希菌；全基因組；多耐藥；耐藥基因；可移動(dòng)遺傳元件；外排泵基因

每年超過(guò)800萬(wàn)美國(guó)女性受尿道致病性埃希菌(UPEC)感染而產(chǎn)生泌尿道感染(UTI)癥狀[1]，33%女性在24歲時(shí)已有UTI史并接受抗生素治療，超過(guò)半數(shù)的女性在一生中曾有過(guò)UTI史[2]，大于80%的UTI由UPEC引起[3]。國(guó)內(nèi)對(duì)UPEC的研究較少，我們?cè)鴪?bào)告28株多耐藥UPEC的耐藥基因和毒力基因等研究[4-7]。為深入了解UPEC的遺傳學(xué)背景，特以Illumina HiSeq與Ion Torrent PGM并行測(cè)序儀對(duì)一株多耐藥大腸埃希菌NB8株做全基因組分析,并從全基因組序列中挖掘耐藥基因、可移動(dòng)遺傳元件和外排泵基因，再和9株已完成全基因組測(cè)序的大腸埃希菌的耐藥基因作比較，現(xiàn)報(bào)告如下。

1 材料與方法

1.1 菌株 大腸埃希菌NB8株分離自2012年4月11日寧波市第一醫(yī)院一例83歲女性患者的尿液。菌株先以法國(guó)生物梅里埃公司VITEK 2-Compact全自動(dòng)微生物鑒定儀及配套革蘭陰性菌鑒定卡鑒定，再以16S rDNA和gyrA基因測(cè)序,經(jīng)GenBank比對(duì)確認(rèn)為大腸埃希菌。

1.2 藥物敏感性試驗(yàn) 用上述的微生物鑒定儀及配套的藥敏卡做16種藥物敏感性檢測(cè)；紅霉素、四環(huán)素、美羅培南、頭孢哌酮/舒巴坦的藥敏試驗(yàn)采用紙片擴(kuò)散法，M-H瓊脂為法國(guó)生物梅里埃公司產(chǎn)品，藥敏紙片為英國(guó)Oxoid公司產(chǎn)品。再根據(jù)2012年美國(guó)臨床實(shí)驗(yàn)室標(biāo)準(zhǔn)化研究所(CLSI)的標(biāo)準(zhǔn)進(jìn)行抗菌藥物敏感性判斷[8]。NB8株對(duì)氨芐西林(≥32)、頭孢唑啉(≥64)、頭孢曲松(≥64)、頭孢他啶(≥64)、頭孢吡肟(≥64)、氨曲南(≥64)、頭孢替坦(≥64)、氨芐西林/舒巴坦(≥32)、慶大霉素(≥16)、妥布霉素(≥16)、丁胺卡那(≥64)、環(huán)丙沙星(≥4)、左氧氟沙星(≥8)、復(fù)方新諾明(≥320)、紅霉素(7mm)、四環(huán)素(7mm) 均耐藥；對(duì)哌拉西林/他唑巴坦(64)、頭孢哌酮/舒巴坦(17mm)中介；對(duì)亞胺培南(≤1)、美羅培南(27mm)敏感。

1.3 全基因組測(cè)序

1.3.1 測(cè)序策略 樣本DNA采取隨機(jī)散彈式打斷來(lái)構(gòu)建 DNA文庫(kù)，并使用Illumina HiSeq 與 Ion Torrent PGM 兩種平臺(tái)進(jìn)行測(cè)序。組裝之后再利用PCR等方法進(jìn)行補(bǔ)缺口。

1.3.2 測(cè)序流程 同參考文獻(xiàn)[9-10]。

1.3.3 組裝結(jié)果 使用軟件將 clean reads 先進(jìn)行de novo組裝后，得到的 contigs除了利用paired-end information得到 scaffolds之外，再與參考基因組：大腸埃希菌SE15 (Accession number: NC_013654) 進(jìn)行對(duì)照與排列順序，得到推定的一條染色體和一個(gè)質(zhì)粒。全基因組測(cè)序委托中國(guó)臺(tái)灣基龍米克斯生物科技股份有限公司進(jìn)行。

1.3.4 補(bǔ)缺口 針對(duì)缺口設(shè)計(jì)引物，得到PCR產(chǎn)物后再使用 Sanger 測(cè)序法，將大部分缺口補(bǔ)充完整。

1.4 生物信息分析

1.4.1 基因注釋 使用 Glimmer v3.02與 prokaryotic GeneMark.hmm v2.10f對(duì)推定的染色體與contigs 進(jìn)行 ORF 預(yù)測(cè)(Start codons 使用 ATG & GTT; Stop codons 使用 TAG & TGA & TAA)，并利用RBSfinder程序找出基因的轉(zhuǎn)錄起始點(diǎn)附近的核糖體連接位點(diǎn)。使用 RNAmmer (v1.2)與 tRNAscan-SE (v1.23)，進(jìn)行 rRNA & tRNA 基因的預(yù)測(cè)。

1.4.2 功能注釋?zhuān)簩㈩A(yù)測(cè)出的基因段落轉(zhuǎn)成蛋白質(zhì)序列分別與NCBI 的 nr & microbial RefSeq protein databases & COG databases 以及 KEGG的protein database進(jìn)行 BLASTP 比對(duì)，得到基因注釋。

2 結(jié) 果

2.1 全基因組組裝結(jié)果 多耐藥大腸埃希菌NB8株經(jīng)Illumina HiSeq與Ion Torrent PGM測(cè)序儀進(jìn)行平行測(cè)序,二種儀器測(cè)得序列合并組裝后用PCR法補(bǔ)缺口。最終得到一條推定的染色體序列,長(zhǎng)4 550 369 bp(內(nèi)含14個(gè)缺口)，和一條推定的質(zhì)粒序列,長(zhǎng)635 377 bp(內(nèi)含33個(gè)缺口)。美國(guó)GenBank登陸號(hào)：LBIS00000000。共挖掘出5 248個(gè)基因，5 023個(gè)CDS(蛋白質(zhì)編碼區(qū))，109個(gè)假基因，18個(gè)rRNAs (包含5S、16S、23S)，85個(gè)tRNAs，13個(gè)ncRNA。

2.2 全基因組功能注釋結(jié)果 COG基因功能分類(lèi)提示NB8株共有3 922個(gè)功能基因，按功能分類(lèi)，可分為3類(lèi)。第1類(lèi)為代謝類(lèi)基因，1 643個(gè)(41.89%)，其中糖類(lèi)轉(zhuǎn)運(yùn)和代謝基因371個(gè)(9.46%)，氨基酸轉(zhuǎn)運(yùn)和代謝基因359個(gè)(9.15%)，產(chǎn)能和儲(chǔ)能基因275個(gè)(7.01%)，無(wú)機(jī)鐵轉(zhuǎn)運(yùn)和代謝基因245個(gè)(6.25%)，輔酶轉(zhuǎn)運(yùn)和代謝基因147個(gè)(3.75%)。第2類(lèi)為信息儲(chǔ)存和處理基因，629個(gè)(16.03%)，其中轉(zhuǎn)錄相關(guān)基因242個(gè)(6.17%)，復(fù)制、重組和修復(fù)基因208個(gè)(5.30%)，翻譯、染色體結(jié)構(gòu)和起源基因177個(gè)(4.51%)。第3類(lèi)為細(xì)胞過(guò)程和信號(hào)轉(zhuǎn)導(dǎo)基因，860個(gè)(21.92%)，其中細(xì)胞壁/膜/外殼起源基因239個(gè)(6.09%)，翻譯修飾，蛋白翻轉(zhuǎn)，分子伴侶基因基因142個(gè)(3.62%)，信號(hào)轉(zhuǎn)導(dǎo)基因137個(gè)(3.49%)，細(xì)胞內(nèi)運(yùn)輸，分泌和小泡運(yùn)輸基因129個(gè)(3.29%)，細(xì)胞運(yùn)動(dòng)基因122個(gè)(3.11%)。另有功能預(yù)測(cè)基因460個(gè)，功能未知基因330個(gè)。用KEGG 反應(yīng)代謝通路分析也得到類(lèi)似的結(jié)果。

同源蛋白來(lái)源物種表提示NB8株來(lái)源于大腸埃希菌的蛋白有5 003個(gè)，占97.5%，還有少量的蛋白可能來(lái)源于肺炎克雷伯菌、沙門(mén)菌和志賀菌等細(xì)菌，需要做進(jìn)一步的群體遺傳學(xué)分析。

通過(guò)全基因組數(shù)據(jù)挖掘，NB8株對(duì)β-內(nèi)酰胺類(lèi)、氨基糖苷類(lèi)、喹諾酮類(lèi)、大環(huán)內(nèi)酯類(lèi)、磺胺類(lèi)、四環(huán)素類(lèi)的耐藥表型和基因型一一對(duì)應(yīng)。①NB8株對(duì)β-內(nèi)酰胺類(lèi)的耐藥表型：對(duì)氨芐西林、頭孢唑啉、頭孢曲松、頭孢他啶、頭孢吡肟、氨曲南、頭孢替坦、氨芐西林/舒巴坦耐藥，對(duì)哌拉西林/他唑巴坦、頭孢哌酮/舒巴坦中介，對(duì)亞胺培南、美羅培南敏感。基因型檢出產(chǎn)β-內(nèi)酰胺酶基因(3個(gè)A類(lèi)酶：blaTEM-1、2個(gè)blaCTX-M-3，1種B類(lèi)酶：假定的金屬β內(nèi)酰胺酶，1種C類(lèi)酶：ampC)；PBP靶位挖掘出pbp1B、pbp1C、pbp2、pbp4、pbp5/6、pbp7、pbpAmpH，膜孔蛋白挖掘出ompC、ompE、ompF，但無(wú)法在各大數(shù)據(jù)庫(kù)檢索到突變位點(diǎn)的資料，所以無(wú)法確定相關(guān)位點(diǎn)是否有突變。②NB8株對(duì)氨基糖苷類(lèi)的耐藥表型：對(duì)慶大霉素、妥布霉素、丁胺卡那均耐藥?；蛐蜋z出氨基糖苷類(lèi)修飾酶(包括1種核苷轉(zhuǎn)移酶基因：ant(3”)-Ⅰ、2種磷酸轉(zhuǎn)移酶基因：aph(3”)-Ⅰb、aph(6)-Ⅰd、1種乙酰轉(zhuǎn)移酶：aac(3)-Ⅱ)。③NB8株對(duì)喹諾酮類(lèi)的耐藥表型：對(duì)環(huán)丙沙星、左氧氟沙星耐藥。基因型檢出gyrA基因80位TCG(S)突變成TTG(L)，87位GAC(D)突變成AAC(N)；parC基因80位 AGC(S)突變成ATT(I)；挖掘出parE，但無(wú)法在各大數(shù)據(jù)庫(kù)檢索到突變位點(diǎn)的資料，所以無(wú)法確定相關(guān)位點(diǎn)是否有突變。④NB8株對(duì)大環(huán)內(nèi)酯類(lèi)的耐藥表型：對(duì)紅霉素耐藥?；蛐蜋z出2種2′-磷酸轉(zhuǎn)移酶：mph(2′)-Ⅰ和mph(2′)-Ⅰ調(diào)節(jié)子、1種大環(huán)內(nèi)酯類(lèi)轉(zhuǎn)運(yùn)ATP 結(jié)合/通透酶：macA。⑤NB8株對(duì)磺胺類(lèi)的耐藥表型：對(duì)復(fù)方新諾明耐藥?；蛐蜋z出2種二氫葉酸合成酶：sul1、sul2，1種二氫葉酸還原酶：dfrA17。6)NB8株對(duì)四環(huán)素類(lèi)的耐藥表型：對(duì)四環(huán)素耐藥?；蛐蜋z出3種耐四環(huán)素耐藥基因：tetA、2個(gè)tetC。7)NB8株對(duì)多粘菌素的耐藥表型；對(duì)多粘菌素耐藥。基因型檢出了arnA。

通過(guò)全基因組數(shù)據(jù)挖掘，NB8株共挖掘出20個(gè)接合性質(zhì)粒遺傳標(biāo)記(traA、traB、traC、traE、traF、traJ、traK、traL、traN、traP、traR、traU、traW、trbC、trbD、trbE、trbI和3個(gè)假定的接合轉(zhuǎn)移基因)、50個(gè)轉(zhuǎn)座酶基因遺傳標(biāo)記(YhgA-like轉(zhuǎn)座酶基因、insG、insI、insD、12個(gè)轉(zhuǎn)座酶基因、34個(gè)假定的轉(zhuǎn)座酶基因)、18個(gè)插入序列遺傳標(biāo)記(IS3、IS5、IS6、IS10、IS26、IS66、IS1294、IS6100、ISSd1、4個(gè)IS2、2個(gè)IS629、3個(gè)IS911)、25個(gè)整合子和噬菌體整合酶遺傳標(biāo)記(intI、整合酶基因xerC、整合酶基因xerD、前噬菌體 DLP12整合酶基因、噬菌體CP-933R 整合酶基因、前噬菌體CP-933X 整合酶基因、噬菌體整合酶基因、8個(gè)整合酶基因、2個(gè)假定的整合酶基因、5個(gè)假定的前噬菌體P4 整合酶基因、3個(gè)假定的噬菌體整合酶基因)、24個(gè)外排泵基因(1個(gè)ATP結(jié)合盒轉(zhuǎn)運(yùn)體家族(ABC)、1個(gè)多藥物與毒物外排家族(MATE)、1個(gè)小多藥外排家族(SMR)、12個(gè)主要易化因子超家族(MFS)、9個(gè)耐受-生節(jié)-分裂家族(RND))。

NB8 株染色體基因組片段序號(hào)5 048至5 052中發(fā)現(xiàn)1個(gè)Ⅰ類(lèi)整合子，包含dfrA17、ant(3”)-Ⅰ、SMR family、sul1 4個(gè)耐藥基因盒，見(jiàn)圖1。

圖1 多耐藥大腸埃希菌NB8株染色體基因組發(fā)現(xiàn)一個(gè)Ⅰ類(lèi)整合子

Fig.1 Class 1 integron in chromosome of multi-drug resistantE.coliNB8 by whole genome sequencing

2.3 比較基因組學(xué)研究 NB8株和9株已完成全基因組測(cè)序的大腸埃希菌耐藥基因的比對(duì)，見(jiàn)表1。10株大腸埃希菌攜帶外排泵基因的陽(yáng)性率非常高，這和外排泵首要功能與細(xì)胞的正常物質(zhì)轉(zhuǎn)移和代謝有關(guān)，而外排抗生素等毒性物質(zhì)只不過(guò)是其次要和偶然的功能。但無(wú)法確定這些外排泵基因的實(shí)際表達(dá)水平。

表1 多耐藥大腸埃希菌NB8株和9株已完成全基因組測(cè)序的大腸埃希菌耐藥基因的比較基因組學(xué)研究
Tab.1 Comparative genomics of genes in multi-drug resistant to E. coli NB8 and other 9 strains of E. coli by whole genome sequencing

耐藥基因種類(lèi)Genetictype基因名稱(chēng)GeneNB8UTI89CFT073536E24377AHSK12O157∶H7EDL933O157∶H7str.SakaiAPECO1耐β-內(nèi)酰胺類(lèi)resistanttobeta-lactamsblaCTX-M+ampC++++++++耐氨基糖苷類(lèi)resistanttoaminoglycosidesaac(3)-Ⅱ+?+ant(3”)-Ⅰ+?+aph(3”)-Ⅰb++aph(6)-Id++耐大環(huán)內(nèi)酯類(lèi)resistanttomacrolidesmacA++++++++++耐磺胺類(lèi)resistanttosulfanilamidessul1++sul2+耐多粘菌素類(lèi)resistanttopolymyxinarnA++++++++++耐四環(huán)類(lèi)resistanttotetracyclicstetC++外排泵effluxpumpstolc++++++++++macb+++++++++mdte++++++++++mdtf++++++++++mdtg++++++++++mdth++++++++++mdtk++++++++++mdtl++++++++++mdtm+++++mdtn++++++++++mdto++++++++++mdtp++++++++++acra++++++++acrb++++++++++

3 討 論

與傳統(tǒng)的Sanger測(cè)序技術(shù)相比,新一代測(cè)序(next-generation sequencing,NGS)技術(shù)具有超高速、高通量、低成本和高效益的優(yōu)勢(shì)，理論上可測(cè)出細(xì)菌的全部染色體和質(zhì)粒序列,可讀出全部功能基因,包括臨床上深度關(guān)切的耐藥基因和毒力基因。故NGS技術(shù)在國(guó)內(nèi)外已逐步用在細(xì)菌耐藥機(jī)制研究和耐藥細(xì)菌傳播途徑分析中，如2011年5月份德國(guó)O104∶H4E.coli引起腹瀉爆發(fā)流行的兩株流行株TY2482株[11]和LB226692[12]株，再?gòu)娜蚪M序列中挖掘出耐藥基因和毒力基因，指導(dǎo)臨床治療。但NGS技術(shù)的缺憾為測(cè)序讀長(zhǎng)(reads)較短,通常會(huì)包含錯(cuò)誤。當(dāng)基因組中出現(xiàn)大量短重復(fù)序列時(shí)，無(wú)法準(zhǔn)確拼接，因此不能測(cè)通全長(zhǎng)而出現(xiàn)缺口。多耐藥大腸埃希菌NB8株經(jīng)Illumina HiSeq與Ion Torrent PGM兩種儀器測(cè)序、組裝后用PCR法補(bǔ)缺口，得到一條推定的染色體序列(內(nèi)含14個(gè)缺口)和一條推定的質(zhì)粒序列(內(nèi)含33個(gè)缺口)。

NB8株的遺傳學(xué)背景明確。NB8株對(duì)β-內(nèi)酰胺類(lèi)抗生素耐藥是因?yàn)楂@得了blaTEM-1、ampC、雙拷貝blaCTX-M-3等β-內(nèi)酰胺酶基因。對(duì)氨基糖苷類(lèi)抗生素耐藥是因?yàn)楂@得了ant(3”)-Ⅰ、aph(6)-Ⅰd、aph(3”)-Ⅰb、aac(3)-Ⅱ等氨基糖苷類(lèi)修飾酶基因。對(duì)喹諾酮類(lèi)抗生素耐藥是因?yàn)間yrA基因83位TCG(S)突變成TTG(L)，87位GAC(D)突變成AAC(N)，parC基因80位 AGC(S)突變成ATT(I)。對(duì)大環(huán)內(nèi)酯類(lèi)抗生素耐藥是因?yàn)楂@得了mph(2′)-Ⅰ調(diào)節(jié)子、mph(2′)-Ⅰ、macA。對(duì)磺胺類(lèi)抗生素耐藥是因?yàn)楂@得了二氫葉酸合成酶sul1、sul2和二氫葉酸還原酶dfrA17。對(duì)四環(huán)素類(lèi)抗生素耐藥是因?yàn)楂@得了tetA和雙拷貝的tetC。對(duì)多粘菌素耐藥是因?yàn)楂@得了arnA。并且還挖掘出大量介導(dǎo)基因水平轉(zhuǎn)移的接合性質(zhì)粒、轉(zhuǎn)座子、插入序列、整合子等可移動(dòng)遺傳元件和大量的外排泵基因。

[1] Bush LM, Chaparro-Rojas F, Okeh V, et al. Cumulative clinical experience from over a decade of use of levofloxacin in urinary tract infections: critical appraisal and role in therapy[J]. Etienne J Infect Drug Resist, 2011, 4: 177-189. DOI: 10.2147/IDR.S15610

[2] Fihn S. Acute uncomplicated urinary tract infection in women[J]. N Engl J Med, 2003, 349: 259-266. DOI: 10.1056/NEJMcp030027

[3] Foxman B. The epidemiology of urinary tract infection[J]. Nat Rev Urol, 2010, 7: 653. DOI: 10.1038/nrurol.2010.190

[4]WENG XB, MI ZH. Investigation of beta-lactamse genes, AmpC genes,ompC gene inEscherichiacoliisolated from urine[J]. Chin J Microbiol Immunol,2010,30(12):1083-1084.(in Chinese)

翁幸鐾,糜祖煌.大腸埃希菌尿液分離株β-內(nèi)酰胺酶基因、AmpC 酶基因和膜孔蛋白o(hù)mpC基因研究[J].中華微生物學(xué)和免疫學(xué)雜志,2010,30(12):1083-1084.

[5] WENG XB, MI ZH. New genotype of aac(6′)-Ⅰb ofEscherichiacoliisolated from urine[J]. Chin J Lab Med,2010,33 (2):149-151.(in Chinese)

翁幸鐾,糜祖煌.從尿液中分離出aac(6′)-Ⅰb 基因新亞型的大腸埃希菌[J].中華檢驗(yàn)醫(yī)學(xué)雜志,2010,33 (2):149-151.

[6] WENG XB, MI ZH. A new variant of gyrA gene inEscherichiacoliisolated from urine[J]. Chin J Microbiol Immunol, 2010,30 (1):11-16.(in Chinese)

翁幸鐾,糜祖煌.大腸埃希菌尿液分離株發(fā)現(xiàn)gyrA基因新變異株[J].中華微生物學(xué)和免疫學(xué)雜志,2010,30 (1):11-16.

[7] WENG XB, MI JR. New genotype of virulence geneompT in multidrug-resistantEscherichiacoliisolated from urine[J].Chin J Lab Med, 2012,35(2):176-179.(in Chinese)

翁幸鐾,糜家睿.多耐藥大腸埃希菌尿液分離株毒力基因檢出ompT基因新亞型[J].中華檢驗(yàn)醫(yī)學(xué)雜志,2012,35(2):176-179.

[8] CLSI. Performance standards for antimicrobial susceptibility testing[S]. 22ndInformational Supplement, 2012.

[9] Brozynska M,Furtado A,Henry RJ.Direct chloroplast sequencing: comparison of sequencing platforms and analysis tools for whole chloroplast barcoding[J]. PLoS One.2014, 9(10):e110387.DOI:10.1371/journal.pone. 0110387

[10] Jeong H,Lee DH,Ryu CM,et al.Toward complete bacterial genome sequencing through the combined use of multiple next-Generation sequencing platforms[J]. J Microbiol Biotechnol.2016, 26(1):207-212.DOI: 10.4014/ jmb.1507.07055

[11] Mellmann A, Harmsen D, Cummings CA, et al. Prospective genomic characterization of the German enterohemorrhagicEscherichiacoliO104:H4 outbreak by rapid next generation sequencing technology[J]. PLoS One, 2011, 6(7): e22751. DOI: 10.1371/journal.pone.0022751

[12] Rohde H, Qin J, Cui Y, et al. Open-source genomic analysis of Shiga-toxin-producingE.coliO104:H4[J]. N Engl J Med, 2011, 365(8): 718-724.DOI: 10.1056/NEJMoa1107643

Genes in multi-drug resistantEscherichiacoliNB8 by whole genome sequencing

WENG Xing-bei1, MI Zu-huang2, WANG Chun-xin3, ZHU Jian-ming4

(1.DepartmentofMedicalLaboratory,NingboFirstHospital,Ningbo315010China；2.DepartmentofBioinformation,WuxiClon-GenTechonologyInstitute,Wuxi214026;China；3.DepartmentofMedicalLaboratory,WuxiPeople’sTechnologyHospitalaffiliatedtoNanjingMedicalUniversity,Wuxi214023China；4.DepartmentofMedicalLaboratory,HangzhouYuhangHospitaloftraditionalChineseMedicine,Hangzhou311106China)

We investigated the genetic background of multi-drug resistantEscherichiacoli(E.coli)NB8. Multi-drug resistantE.coliNB8 was collected from urine sample of an inpatient in Ningbo First Hospital, in April, 2012. Sequencing 16S rDNA and gyrA, BLASTn algorithm were performed to identifyE.coli.Then, whole genome were sequenced by Illumina HiSeq and Ion Torrent PGM,and PCR was performed to fill gaps. Finally, comparative genomics of resistant genes in multi-drug resistantE.coliNB8 and other 9 strains ofE.coliwere performed. Results showed that by whole genome sequencing, a putative chromosome sequence (14 gaps inside) was 4 550 369 bp, a putative plasmid sequence (33 gaps inside) was 635 377 bp. Resistant genes to beta-lactams, aminoglycosides, quinolones, macrolides, sulfonamides, tetracyclines, polymyxin were positive. Furthermore, genetic markers of conjugative plasmids, transposons, insertion sequences, integrons and 5 classes of efflux genes were positive. In conclusion, genetic background of multi-drug resistantE.coliNB8 was clear.Enhancing efflux pumps contribute to resistance. Plasmids, transposons and integrons promote resistance among bacteria, and make recipient

Escherichiacoli; whole genome; multi-drug resistance; resistant gene; mobile genetic element; efflux pump gene

1.寧波市第一醫(yī)院檢驗(yàn)科，寧波 315010，wxb6006@hotmail.com；

2.無(wú)錫市克隆遺傳技術(shù)研究所生物信息學(xué)室，無(wú)錫 214026；

10.3969/j.issn.1002-2694.2016.08.011

R378

A

1002-2694(2016)08-0741-05

2015-09-01；

2016-02-15

寧波市自然科學(xué)基金(No.2012A610186)和浙江省中醫(yī)藥基金(No.2011ZB126)資助

3.南京醫(yī)科大學(xué)附屬無(wú)錫人民醫(yī)院醫(yī)學(xué)檢驗(yàn)科，無(wú)錫 214023；

4.杭州市余杭區(qū)中醫(yī)院檢驗(yàn)科，杭州 311106

Supported by the Ningbo Natural Science Foundation (No. 2012A610186) and the Zhejiang Traditional Chinese Medicine foundation (No. 2011ZB126).

bacterium to confer multidrug resistance.