医院相关性ST5型耐甲氧西林金黄色葡萄球菌ZJ5499的全基因测序及序列分析

目的对1例医院获得性耐甲氧西林金黄色葡萄球菌(MRSA)肺炎患者血中分离的1株MRSA(ZJ5499)进行基因组序列信息解析。方法采用Illumina平台高通量测序和Sanger测序相结合对ZJ5499菌株进行全基因组测序,并使用相关软件对序列进行拼接、基因预测、功能注释、直系同源簇注释(COG)及毒力因子和耐药基因分析;并与国内常见流行序列型(ST型)菌株进行进化关系分析、毒力因子和耐药基因比较。结果 ZJ5499菌株基因组大小为2 888 783bp,GC含量32.84%,序列已提交至GenBank数据库,登录号为CP011685。该菌株基因组中含有大量与致病性相关的毒力因子,并含有spc、aadD、mecA、norA和erm(A)五个耐药相关基因,与同一ST型菌株基本一致。进化关系分析显示该菌株与同一ST型菌株关系较近。毒力因子与ST5型菌株没有较大差异,而较其他ST型明显增多。结论本研究报道了临床MRSA菌株ZJ5499的全基因组序列。序列分析发现该菌株的毒力和耐药性与ST5型的菌株相近,而ST5型菌株较其他国内流行ST型菌株携带较多毒力基因。     

创伤弧菌毒力相关因子的全基因组关联分析研究

【目的】创伤弧菌是致死率最高的弧菌物种,但目前尚无在全基因组层面挖掘毒力相关因子的研究。本研究以创伤弧菌分离来源(临床和环境)作为不同表型,通过与260株基因组序列进行关联分析,挖掘毒力相关因子,从而进一步了解创伤弧菌致病因素。【方法】对139株创伤弧菌分离株进行高通量测序,获取其全基因组序列;与公共数据库已公开发表的121株基因组整合,使用pyseer软件进行全基因组关联分析,对与不同分离来源显著相关的基因进行注释和解读。【结果】共发现11个基因与临床分离株显著相关,其中9个是本研究新发现的创伤弧菌潜在毒力相关因子。【结论】本研究使用群体基因组学和统计遗传学方法,在全基因组范围扫描挖掘了创伤弧菌毒力相关因子,为深入揭示该物种致病机制、设计新的疫苗和治疗靶点提供了重要依据。     

应用全基因组PCR扫描方法进行猪链球菌 基因组结构的比较分析

2005年, 在中国四川局部地区爆发人感染猪链球菌疫情, 因其感染人数多, 病死率高引起关注, 为了确定该致病菌是否发生变异, 通过应用全基因组PCR扫描方法(WGPScaning)、多位点序列分型技术(MLST)以及毒力相关基因的序列测定, 比较分别来自本次疫情中的病人和病猪、以前流行期间感染的病人分离的菌株以及网上公布的来自欧洲的猪链球菌基因组序列, 结果显示各菌株的基因组结构相似, 毒力相关基因没有差异, 所有菌株都属于ST1序列群, 说明本次引起四川疫情的菌株未发生明显的基因组结构的改变.     

2016-2017年宁夏地区牛源金黄色葡萄球菌的全基因组分析

[目的] 了解宁夏地区奶牛乳腺炎金黄色葡萄球菌（Staphylococcus aureus,SA）代表菌株的基因组序列基本特征,进一步探究其耐药基因型、毒力及进化关系,为兽医临床防治提供理论依据。[方法] 采用纸片法对97株金黄色葡萄球菌临床分离株进行抗菌药物敏感性试验,同时进行葡萄球菌蛋白A（Staphylococcus aureus protein A,spa）分型、多位点序列分型（multilocus sequence typing,MLST）,根据分型结果选取16株代表菌株进行全基因组测序,并对获得的测序序列进行处理分析。[结果] 药敏试验结果显示97株分离株对18种抗菌药物存在不同程度的耐药,其中9株耐甲氧西林金黄色葡萄球菌（methicillin-resistantStaphylococcus aureus,MRSA）对青霉素、氨苄西林、苯唑西林、头孢噻呋、磺胺异噁唑、红霉素、庆大霉素和克林霉素等8种抗菌药物完全耐药,甲氧西林敏感金黄色葡萄球菌（methicillin-sensitiveStaphylococcus aureus,MSSA）菌株对青霉素、氨苄西林、磺胺异噁唑耐药率较高。耐药基因数据库（antibiotic resistance genes database,ARDB）注释分析显示16株代表菌株共携带21种耐药基因,其中norA、tet38、bacA、mepA的携带率较高,与药敏试验结果具有一定的相关性。毒力基因数据库（virulence factors of pathogenic bacteria,VFDB）注释分析显示所有菌株携带多种与粘附、宿主免疫逃逸、分泌、胞外酶编码、铁摄取等疾病相关的毒力基因,MRSA菌株均携带较多毒力因子,MSSA菌株携带毒力因子数目不等。基因岛预测结果显示16株代表菌株存在不同数量的基因岛且MRSA菌株携带基因岛数目及毒力基因岛较多,但耐药基因岛数目与MSSA差异不明显。SNP分析结果显示部分分离株同源性较高,同源性较高的两株MRSA的全基因组基本序列特征差异较小,携带的耐药、毒力基因情况相似。[结论] 宁夏地区牛源SA分离株耐药性情况严重且具有较高的毒力水平,本研究为家畜相关MRSA（livestock-associated MRSA,LA-MRSA）与MSSA基因组序列信息的比较分析及宁夏地区SA感染的临床防控提供参考依据。     

应用全基因组PCR扫描方法进行猪链球菌 基因组结构的比较分析

2005年, 在中国四川局部地区爆发人感染猪链球菌疫情, 因其感染人数多, 病死率高引起关注, 为了确定该致病菌是否发生变异, 通过应用全基因组PCR扫描方法(WGPScaning)、多位点序列分型技术(MLST)以及毒力相关基因的序列测定, 比较分别来自本次疫情中的病人和病猪、以前流行期间感染的病人分离的菌株以及网上公布的来自欧洲的猪链球菌基因组序列, 结果显示各菌株的基因组结构相似, 毒力相关基因没有差异, 所有菌株都属于ST1序列群, 说明本次引起四川疫情的菌株未发生明显的基因组结构的改变.     

一株羊源A型多杀性巴氏杆菌的全基因组测序及生物学特性分析

[背景] 多杀性巴氏杆菌（Pasteurella multocida,Pm）是一种革兰氏阴性菌,可引起动物和人类的呼吸道疾病和败血症等。本实验室前期分离鉴定一株A型Pm HN02菌株。[目的] 通过对HN02菌株的全基因组测序及生物信息学分析,扩充多杀性巴氏杆菌的基因组数据库信息;通过毒力基因鉴定和系统进化树分析,明确该菌株含有的毒力基因和遗传进化关系,为临床预防和诊断提供理论依据。[方法] 使用单分子实时测序（Single Molecule Real Time Sequencing,SMRT）技术对Pm HN02菌株进行全基因组测序,利用Illumina测序校正后进行基因功能注释和生物信息学分析。使用PCR鉴定菌株毒力基因,并构建进化树进行分析。[结果] Pm HN02菌株全基因组大小为2 333 292 bp,GC含量为40.15mol%,预测到的编码基因有2 389个,包含19个rRNA （6个23S rRNA、6个16S rRNA、7个5S rRNA）、62个tRNA基因、5个sRNA;含84个串联重复序列、66个小卫星DNA、2个微卫星DNA、9个基因岛、9个前噬菌体;分别有1 648、2 190和1 917个基因注释在GO、KEGG和COG数据库中,而且大部分富集于Pm的代谢过程;还有85个III型分泌系统效应蛋白、191个表型突变基因、165个毒力因子相关基因。根据分析结果绘制该菌株的全基因组圈图,并将基因组信息提交至NCBI后获得登录号cp037865。PCR鉴定发现该菌株含有fimA、toxA等14个毒力基因,缺失了tadD等毒力基因。系统进化树分析发现该菌株同北京的Pm3菌株（MH150895.1）进化关系最接近。[结论] 研究完成了A型Pm HN02株的全基因组测序和生物学特性鉴定,揭示了其同国内外Pm分离株的进化关系,为预防Pm疾病流行和探索Pm致病机制提供了参考。     

鼠衣原体毒力基因的筛选和基因差异型菌株的建立

【目的】探讨鼠衣原体(Chlamydia muridarum,Cm)标准株与减毒株全基因组序列中存在的差异,筛选毒力基因并建立不同基因型的单克隆菌株,为后续致病机制的研究奠定基础。【方法】将Cm标准株G0和减毒株G28以高通量测序法进行全基因组测序,通过全基因组序列比对分析并筛选潜在的毒力相关基因;经空斑形成实验(Plaque assay)在混合菌G0和G28中大量挑取空斑,以毒力靶基因的PCR测序鉴定从G0和G28中筛选含有不同毒力基因型的单克隆菌株。【结果】全基因测序结果显示Cm G0与G28的TC0412、TC0237和TC0668基因明显不同。通过挑取空斑初步筛选了111个空斑样品,通过3个毒力基因的PCR测序鉴定最终获得了G0和G28来源的56个单克隆菌株,并根据TC0412蛋白型分成B3、D1和E1三组,每组中含有TC0237和TC0668基因差异性菌株。【结论】Cm的致病能力可能与TC0412、TC0237和TC0668基因有关,筛选获得的单克隆菌株通过分组匹配后可用于后续的毒力基因功能研究。     

大熊猫源致病大肠杆菌CCHTP全基因组测序及耐药和毒力基因分析

致病性大肠杆菌是引起动物泌尿系统感染的重要病原菌,本研究对泌尿生殖道感染出现潜血的大熊猫尿液中分离的一株致病性大肠杆菌(Escherichia coli CCHTP)进行全基因组测序,检测其中耐药基因和毒力因子的情况,同时对基因岛上耐药和毒力基因及其基因环境进行研究。研究发现,大肠杆菌CCHTP中存在多种类型的耐药基因,其中外排泵系统基因数量最多,包括mdf A、emr E和mdt N等介导多重耐药外排泵的基因。此外,该菌还携带166种毒力因子及563个相关毒力基因,其中属于黏附与侵袭类的毒力因子及相关基因数量最多。对19个基因岛分析发现,基因岛GIs011和GIs017中各有一段包含耐药和毒力基因的序列,两侧与可移动遗传元件(转座酶和插入序列)相连,这些结构可能介导耐药及毒力基因水平转移。本研究通过全基因组测序分析了大熊猫源致病性大肠杆菌中存在的耐药及毒力基因情况,对大熊猫相关疾病的科学治疗、合理用药有重要意义。     

一株兔源A型多杀性巴氏杆菌的分离鉴定和全基因组测序及分析

【背景】多杀性巴氏杆菌可导致猪肺疫、牛出血性败血症和兔出血性败血症等多种疾病,严重威胁多种动物畜牧养殖业的健康发展。【目的】重庆某兔场送检一批病死兔,为研究其病原和治疗方法,对病原进行了微生物分离和全基因组测序分析。【方法】从2022年重庆某兔场送检兔病料中进行细菌分离纯化、生化试验、16S rRNA基因鉴定、荚膜血清型分型、药敏试验和毒力基因检测,同时通过全基因组测序结果进行毒力、耐药基因注释和遗传进化等分子生物学信息分析。【结果】该菌为兔源A:ST74多杀性巴氏杆菌,命名为LXSS001,基因组序列上传到NCBI数据库(登录号为CP119523.1),药敏试验显示该菌对四环素、杆菌肽、复方新诺明和磺胺异恶唑耐药,对头孢噻肟、头孢哌酮和丁胺卡那等药物敏感。全基因组长度为2 480 671 bp,并注释到了58个毒力基因和9类药物的靶向抗药基因。通过联合建树表明其与3480株一致性最高。【结论】本研究完成了一株A型多杀性巴氏杆菌的分离鉴定和全基因组测序,并揭示了其与国内外其他分离株的进化关系,为多杀性巴氏杆菌的后续研究提供了参考依据。     

ST59型社区获得性耐甲氧西林金黄色葡萄球菌毒力因子研究

了解ST59型社区获得性耐甲氧西林金黄色葡萄球菌(community-acquired methicillin-resistant Staphy-lococcus aureus,CA-MRSA)携带毒力因子的情况。用PCR方法扩增ST59型CA-MRSA PSMα、PVL、SEA、SEB、SEC、SED、SEE、TSST-1、ETA、ETB基因。5株CA-MRSA全部检出PSMα基因和PVL基因,均未检出SEA、SEB、SEC、SED、SEE、TSST-1、ETA、ETB基因。PSMα和PVL基因是ST59型社区获得性耐甲氧西林金黄色葡萄球菌常见的毒力因子。     

Complete Genome Sequence of a Pantón-Valentine Leukocidin-Negative Community-Associated Methicillin-Resistant Staphylococcus aureus Strain of Sequence type 72 from Korea

In the past decade, community-associated (CA-) infections with methicillin-resistant Staphylococcus aureus (MRSA) have emerged throughout the world. Different CA-MRSA strains dominate in different geographical locations. Many CA-MRSA lineages contain genes coding for the Pantón-Valentine leukocidin. However, the role of this leukotoxin in CA-MRSA pathogenesis is still controversial. The genome sequences of two key PVL-positive CA-MRSA strains (USA300, USA400) have been reported, but we lack information on the more recently found PVL-negative CA-MRSA strains. One such strain is the PVL-negative ST72, the main cause of CA-MRSA infections in Korea. Here, we report the entire genome sequence of CA-MRSA ST72 and analyze its gene content with a focus on virulence factors. Our results show that this strain does not have considerable differences in virulence factor content compared to other CA-MRSA strains (USA300, USA400), indicating that other toxins do not substitute for the lack of PVL in ST72. This finding is in accordance with the notion that differential expression of widespread virulence determinants, rather than the acquisition of additional virulence factors on mobile genetic elements, such as PVL, is responsible for the increased virulence of CA- compared to hospital-associated MRSA.     

The dominant Australian community-acquired methicillin-resistant Staphylococcus aureus clone ST93-IV [2B] is highly virulent and genetically distinct

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) USA300 has spread rapidly across North America, and CA-MRSA is also increasing in Australia. However, the dominant Australian CA-MRSA strain, ST93-IV [2B] appears distantly related to USA300 despite strikingly similar clinical and epidemiological profiles. Here, we compared the virulence of a recent Australian ST93 isolate (JKD6159) to other MRSA, including USA300, and found that JKD6159 was the most virulent in a mouse skin infection model. We fully sequenced the genome of JKD6159 and confirmed that JKD6159 is a distinct clone with 7616 single nucleotide polymorphisms (SNPs) distinguishing this strain from all other S. aureus genomes. Despite its high virulence there were surprisingly few virulence determinants. However, genes encoding α-hemolysin, Panton-Valentine leukocidin (PVL) and α-type phenol soluble modulins were present. Genome comparisons revealed 32 additional CDS in JKD6159 but none appeared to encode new virulence factors, suggesting that this clone's enhanced pathogenicity could lie within subtler genome changes, such as SNPs within regulatory genes. To investigate the role of accessory genome elements in CA-MRSA epidemiology, we next sequenced three additional Australian non-ST93 CA-MRSA strains and compared them with JKD6159, 19 completed S. aureus genomes and 59 additional S. aureus genomes for which unassembled genome sequence data was publicly available (82 genomes in total). These comparisons showed that despite its distinctive genotype, JKD6159 and other CA-MRSA clones (including USA300) share a conserved repertoire of three notable accessory elements (SSCmecIV, PVL prophage, and pMW2). This study demonstrates that the genetically distinct ST93 CA-MRSA from Australia is highly virulent. Our comparisons of geographically and genetically diverse CA-MRSA genomes suggest that apparent convergent evolution in CA-MRSA may be better explained by the rapid dissemination of a highly conserved accessory genome from a common source.     

Single-molecule sequencing reveals the molecular basis of multidrug-resistance in ST772 methicillin-resistant Staphylococcus aureus

Background

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of hospital-associated infection, but there is growing awareness of the emergence of multidrug-resistant lineages in community settings around the world. One such lineage is ST772-MRSA-V, which has disseminated globally and is increasingly prevalent in India. Here, we present the complete genome sequence of DAR4145, a strain of the ST772-MRSA-V lineage from India, and investigate its genomic characteristics in regards to antibiotic resistance and virulence factors.

Results

Sequencing using single-molecule real-time technology resulted in the assembly of a single continuous chromosomal sequence, which was error-corrected, annotated and compared to nine draft genome assemblies of ST772-MRSA-V from Australia, Malaysia and India. We discovered numerous and redundant resistance genes associated with mobile genetic elements (MGEs) and known core genome mutations that explain the highly antibiotic resistant phenotype of DAR4145. Staphylococcal toxins and superantigens, including the leukotoxin Panton-Valentinin Leukocidin, were predominantly associated with genomic islands and the phage φ-IND772PVL. Some of these mobile resistance and virulence factors were variably present in other strains of the ST772-MRSA-V lineage.

Conclusions

The genomic characteristics presented here emphasize the contribution of MGEs to the emergence of multidrug-resistant and highly virulent strains of community-associated MRSA. Antibiotic resistance was further augmented by chromosomal mutations and redundancy of resistance genes. The complete genome of DAR4145 provides a valuable resource for future investigations into the global dissemination and phylogeography of ST772-MRSA-V.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1599-9) contains supplementary material, which is available to authorized users.     

鲍曼不动杆菌携带前噬菌体的生物信息学分析

[目的] 预测并分析82株全基因组测序的鲍曼不动杆菌前噬菌体的携带情况，鉴定前噬菌体编码的抗生素耐药基因和毒力因子。[方法] 利用PHASTER （Phage Search Tool Enhanced Release）软件预测鲍曼不动杆菌携带的前噬菌体，采用CARD （The Comprehensive Antibiotic Research Database）和VFDB （Virulence Factors Database）在线分析软件预测前噬菌体编码的抗生素耐药基因和毒力因子。[结果] 预测到472条鲍曼不动杆菌前噬菌体，其中完整型前噬菌体201条，疑似型前噬菌体91条，缺陷型前噬菌体180条。平均每株鲍曼不动杆菌基因组中可携带至少2条完整型前噬菌体。每株鲍曼不动杆菌所携带的全部前噬菌体占其基因组比例约为4%-6%。29条前噬菌体携带77个耐药基因，耐药表型共有14种，分别来自15个不同的家族，涵盖6种抗生素耐药的作用机制。132条前噬菌体编码毒力基因，归类为38种毒力基因和34种毒力因子。不同类型的前噬菌体普遍携带1-2种毒力因子，少数前噬菌体携带3种及以上毒力因子。分析毒力因子可能的宿主来源构成比发现，除鲍曼不动杆菌外，脑膜炎奈瑟菌、痢疾志贺氏菌、嗜肺军团菌及其亚种等也有较高的结构比例，是可能的宿主来源。[结论] 鲍曼不动杆菌普遍携带前噬菌体，但前噬菌体基因在鲍曼不动杆菌基因组中所占比例不高。部分前噬菌体携带抗生素耐药基因，以氨基糖苷类、磺胺类及β-内酰胺类耐药为主。约30%的前噬菌体携带毒力基因。前噬菌体可能在鲍曼不动杆菌抗生素耐药性的获得、传播及致病性演变中发挥重要作用。     

一株牛源都柏林沙门氏菌的全基因组测序及毒力与耐药性分析

【背景】沙门氏菌(Salmonella spp.)是重要的人畜共患病原菌,其毒力和耐药性的不断增强引起广泛关注。【目的】了解从通辽市一犊牛死亡病例中所分离牛源都柏林沙门氏菌的毒力及耐药性情况。【方法】以病死犊牛肺脏为材料,经细菌分离纯化及16S rRNA基因测序,鉴定病原为沙门氏菌。采用动物试验、药敏试验和PCR方法对分离菌进行毒力、耐药性,以及毒力基因和耐药基因检测,并对其进行全基因组测序分析。【结果】分离菌具有较强毒力,对小鼠半数致死量为2.8×10⁶ CFU/mL。分离菌为多重耐药菌,仅对多粘菌素B和噻孢霉素敏感,对强力霉素和恩诺沙星中度敏感。检测13种沙门氏菌常见毒力基因,检出率为92.3%。对分离菌进行全基因组测序分析,该菌株为都柏林沙门氏菌,基因组大小为4 965 370 bp,GC含量为52.12%,同时携带2个质粒,大小分别为79 524 bp (pTLS-1)和45 301 bp (pTLS-2)。分离菌中共携带996个毒力基因和24个毒力岛;共携带42个耐药基因,其中4个为可水平转移基因,基因组中存在9个可移动遗传元件,包括插入序列和转座子等。【结论】分离牛源都柏林沙门氏菌菌株具有较强毒力且为多重耐药株,携带大量毒力基因及耐药基因。     

Molecular nature of methicillin-resistant Staphylococcus aureus derived from explosive nosocomial outbreaks of the 1980s in Japan

Community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA) with Panton-Valentine leukocidin (PVL) genes is increasing worldwide. Nosocomial outbreak-derived (hospital-acquired) MRSA (HA-MRSA) in Japan in the 1980s was also largely PVL(+). PVL(+) HA-MRSA and CA-MRSA shared the same multi-locus sequence type (ST30) and methicillin resistance cassette (SCCmecIV), but were divergent in oxacillin resistance, spa typing, PFGE analysis or clfA gene analysis. PVL(+) HA-MRSA, which probably originated in PVL(+)S. aureus ST30, was highly adhesive (carrying cna and bbp genes), highly-toxic (carrying luk(PV) and sea genes) and highly drug-resistant. PVL(+) HA-MRSA was once replaced by other PVL(-) HA-MRSA (e.g., ST5), and is re-emerging as CA-MRSA.     

Analysis of ultra low genome conservation in Clostridium difficile

Microarray-based comparative genome hybridisations (CGH) and genome sequencing of Clostridium difficile isolates have shown that the genomes of this species are highly variable. To further characterize their genome variation, we employed integration of data from CGH, genome sequencing and putative cellular pathways. Transcontinental strain comparison using CGH data confirmed the emergence of a human-specific hypervirulent cluster. However, there was no correlation between total toxin production and hypervirulent phenotype, indicating the possibility of involvement of additional factors towards hypervirulence. Calculation of C. difficile core and pan genome size using CGH and sequence data estimated that the core genome is composed of 947 to 1,033 genes and a pan genome comprised of 9,640 genes. The reconstruction, annotation and analysis of cellular pathways revealed highly conserved pathways despite large genome variation. However, few pathways such as tetrahydrofolate biosynthesis were found to be variable and could be contributing to adaptation towards virulence such as antibiotic resistance.     

Genetic Variation in Spatio-Temporal Confined USA300 Community-Associated MRSA Isolates: A Shift from Clonal Dispersion to Genetic Evolution?

Introduction

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) are increasingly isolated, with USA300-0114 being the predominant clone in the USA. Comparative whole genome sequencing of USA300 isolates collected in 2002, 2003 and 2005 showed a limited number of single nucleotide polymorphisms and regions of difference. This suggests that USA300 has undergone rapid clonal expansion without great genomic diversification. However, whole genome comparison of CA-MRSA has been limited to isolates belonging to USA300. The aim of this study was to compare the genetic repertoire of different CA-MRSA clones with that of HA-MRSA from the USA and Europe through comparative genomic hybridization (CGH) to identify genetic clues that may explain the successful and rapid emergence of CA-MRSA.

Materials and Methods

Hierarchical clustering based on CGH of 48 MRSA isolates from the community and nosocomial infections from Europe and the USA revealed dispersed clustering of the 19 CA-MRSA isolates. This means that these 19 CA-MRSA isolates do not share a unique genetic make-up. Only the PVL genes were commonly present in all CA-MRSA isolates. However, 10 genes were variably present among 14 USA300 isolates. Most of these genes were present on mobile elements.

Conclusion

The genetic variation present among the 14 USA300 isolates is remarkable considering the fact that the isolates were recovered within one month and originated from a confined geographic area, suggesting continuous evolution of this clone.     

Multidrug resistance and high virulence genotype in uropathogenic Escherichia coli due to diffusion of ST131 clonal group producing CTX-M-15: an emerging problem in a Tunisian hospital

A collection of 201 Escherichia coli strains isolated from urine of patients in a Tunisian hospital between January 2006 and July 2008 was studied. Microbial identification was done by conventional methods, and antibiotic susceptibility with disk diffusion method was performed according to the Clinical Laboratory and Standards Institute guidelines. Detection of extended-spectrum beta-lactamase (ESBL) was performed by double-disk synergy test (DDST) and identification was done by PCR and sequencing. ESBL-producing isolates were subjected to molecular typing by random amplified polymorphic DNA (RAPD) and ST131 detection by PCR. Four phylogenetic groups (A, B1, B2 and D), 18 virulence genes and CTX-M group were individualized using PCR. Statistical analysis was done by Pearson χ2 test and Mann–Whitney U test. The strains were recovered primarily from urology (28 %), maternity (19 %) and medicine (16 %) wards. Antibiotic resistance rates were ampicilin (72.1 %), nalidixic acid (41.8 %), ciprofloxacin (38.8 %), gentamicin (23.9 %) and cefotaxime (17.4 %). Thirty-one of cefotaxime-resistant isolates (n?=?35) had a positive DDST and harboured bla _CTX-M-15 gene. Twenty of them (64.5 %) belonged to the ST131 clone and showed the same RAPD DNA profile. Ciprofloxacin- and cotrimoxazole-susceptible isolates were significantly associated with phylogenetic group B2, whereas isolates that were resistant to these molecules were associated with B1 and D phylogenetic groups, respectively. Virulence genes were significantly more frequent among ciprofloxacin- and cotrimoxazole-susceptible strains than those resistant to these antibiotics. However, CXT-M-15-producing isolates were associated with many virulence genes. Isolates concomitantly susceptible to the three antimicrobials agents (ciprofloxacin, cefotaxime and cotrimoxazole) were significantly associated with group B2 and high virulence score, whereas isolates with resistance patterns especially those including resistance to ciprofloxacin belonged predominantly to B1 phylogroup and haboured few virulence genes. The emergence of virulent and multidrug-resistant E. coli is a concerning development that deserves close attention in our institution.