中国农业科学 ›› 2022, Vol. 55 ›› Issue (7): 1469-1478.doi: 10.3864/j.issn.0578-1752.2022.07.017

• 畜牧·兽医·资源昆虫 • 上一篇    

多重耐药大肠杆菌中前噬菌体的分布特征及诱导分离

刘教(),刘畅,陈进,王勉之,熊文广,曾振灵()   

  1. 华南农业大学/广东省兽药研制与安全评价重点实验室/国家兽医微生物耐药性风险评估实验室,广州 510642
  • 收稿日期:2021-02-06 接受日期:2021-09-30 出版日期:2022-04-01 发布日期:2022-04-18
  • 通讯作者: 曾振灵
  • 作者简介:刘教,E-mail: 1766612575@qq.com
  • 基金资助:
    国家自然科学基金面上项目(31872524);广东省普通高校省级重大科研项目(2017KZDXM006)

Distribution Characteristics of Prophage in Multidrug Resistant Escherichia coli as well as Its Induction and Isolation

LIU Jiao(),LIU Chang,CHEN Jin,WANG MianZhi,XIONG WenGuang,ZENG ZhenLing()   

  1. South China Agriculture University/Guangdong Provincial Key Laboratory of Veterinary Pharmaceutics Development and Safety Evaluation/National Risk Assessment Laboratory for Antimicrobial Resistance of Animal Original Bacteria, Guangzhou 510642
  • Received:2021-02-06 Accepted:2021-09-30 Online:2022-04-01 Published:2022-04-18
  • Contact: ZhenLing ZENG

摘要:

【目的】通过调查前噬菌体在多重耐药大肠杆菌中的分布特征、诱导分离以及前噬菌体中耐药基因与毒力基因的流行状况,为研究前噬菌体介导耐药基因在细菌的传播提供科学依据。【方法】挑选前期保存的2018—2019年广东省分离的131株禽源多重耐药大肠杆菌进行核酸提取及全基因组测序,将二代测序的结果组装拼接成全基因组序列,上传至噬菌体PHASTER网络数据库与数据库中已有的噬菌体基因组序列进行比对分析。利用CGE数据库比对耐药基因与毒力基因,从而获得在前噬菌体上耐药基因与毒力基因的分布情况。温和性噬菌体由丝裂霉素C诱导并使用双层平板法分离纯化。【结果】131株大肠杆菌药物敏感性试验的结果显示,氨苄西林、四环素、氟苯尼考、复方新诺明的耐药率均高达90%以上,其次是头孢类抗生素以及庆大霉素、环丙沙星、美罗培南和黏菌素均在50%左右,替加环素的耐药率达到了0.2%,所有菌株都呈现出多重耐药的现象,均为多重耐药大肠杆菌。131株多重耐药大肠杆菌中共检出736个前噬菌体片段,其中包含329个完整型前噬菌体,其与40个已知数据库噬菌体物种以不同百分比匹配上;可疑型噬菌体有66个,其与20个已知数据库噬菌体物种以不同百分比匹配上;不完整型噬菌体有341个,其与52个已知数据库噬菌体物种以不同百分比匹配上,完整型前噬菌体的基因序列显示出与已知的噬菌体物种的序列相似性最高,平均为58.53%;131株大肠杆菌中平均前噬菌体数量为5.6个,平均总含量为152.4 kb。前噬菌体基因组占其宿主基因组的比例分布在0.58%—5.87%,以3.0%为主。前噬菌体基因组长度范围在2.8—107.9 kb,其中13.0 kb的前噬菌体出现的频次最高,占所有前噬菌体的9.1%。CGE比对结果表明,131株多重耐药大肠杆菌的基因组共在18株前噬菌体序列检测到耐药基因mdf(A)、lnu(G)和mcr-1,其中mdf(A)、lnu(G)和mcr-1检出数分别为16、1和1。71株多重耐药大肠杆菌前噬菌体中携带有6种不同的毒力基因,其中存在部分菌株携带2种或者3种毒力基因,有62株前噬菌体携带端粒酶RNA基因terC,16株前噬菌体携带血清存活率增加基因iss,外膜蛋白酶ompT、黏附素基因ihacvaC和ABC转运蛋白基因mchF分别在2、2、1和1株前噬菌体中检出。mdf(A)和terC分别是前噬菌体中最常见的耐药基因和毒力基因。温和性噬菌体诱导试验结果显示,前噬菌体的诱导成功率为84.0%,但出现噬菌斑的概率仍比较低。【结论】前噬菌体在多重耐药大肠杆菌中分布广泛且携带有多种耐药基因和毒力基因,温和性噬菌体诱导成功率高,具有携带耐药基因及毒力基因水平传播的风险,需要加强和持续监测。

关键词: 前噬菌体, 大肠杆菌, 诱导, 分布特征, 耐药基因

Abstract:

【Objective】 This study investigated the distribution characteristics of prophage in multi-drug resistant Escherichia coli, induction and isolation, as well as the prevalence of drug resistance and virulence genes in prophage, so as to provide a scientific basis for the study of prophage-mediated resistance genes in the spread of bacteria. 【Method】 131 multi-drug resistant E. coli isolating from poultry origin in Guangdong Province from 2018 to 2019 were selected in the laboratory for nucleic acid extraction and whole-genome sequencing. The results of second-generation sequencing were assembled and spliced into a whole-genome sequence and uploaded to the phage. The PHASTER network database was compared and analyzed with the existing phage genome sequences in the database. Drug resistance genes and virulence genes were compared on the CGE database, and then the distribution of drug resistance genes and virulence genes on the prophage were obtained. The mild phage was induced by mitomycin C, separated and purified by using the double-layer plate method. 【Result】 The results of the drug sensitivity test of 131 strains of Escherichia coli showed that the drug resistance rates of ampicillin, tetracycline, florfenicol and compound trimethoprim were all more than 90%, followed by cephalosporin antibiotics, gentamicin, ciprofloxacin, meropenem and colistin with all around 50%, and the resistance rate of tigecycline reached 0.2%. All strains showed multi-drug resistance, and they were all multi-drug resistant Escherichia coli. A total of 736 prophage fragments were detected in 131 strains of multi-drug resistant E. coli, including 329 complete prophage, 66 suspicious phages and 341 incomplete phage, which matched with 40, 20 and 52 known database phage species in different percentages, respectively. The gene sequence of the complete prophage showed that it matched the known phage species better, and the sequence similarity was the highest, with an average of 58.53%. The average number of prophages in 131 strains of E. coli was 5.6, and the average total content was 152.4 kb. Prophage genome accounted for 0.58% to 5.87% of its host genome, with 3.0% being the dominant. The length of the prophage genome ranged from 2.8 to 107.9 kb, and the 13.0 kb prophage had the highest frequency, accounting for 9.1% of all prophages. CGE comparison results showed that the genomes of 131 strains of multi-drug-resistant E. coli detected resistance genes mdf (A), lnu (G) and mcr-1 in 18 prophage sequences. The detected numbers of mdf (A), lnu (G) and mcr-1 were 16, 1, and 1, respectively. 71 strains of multi-drug resistant E. coli prophage carried 6 different virulence genes, and some strains carried 2 or 3 virulence genes. There were 62 prophages carrying the telomerase RNA gene terC, 16 prophages carrying the serum survival increasing gene iss, and the outer membrane protease ompT, among which the adhesin gene iha, the cvaC gene and the ABC transporter gene mchF were at 2, 2, 1, and 1, respectively. Mcr-a gene were detected in prophage of 1 strain multi-drug resistant Escherichia coli. The mdf (A) gene and terC gene were the most common resistance genes and virulence genes in prophage, respectively. The results of mild phage induction experiments showed that the success rate of prophage induction was 84.0%, but the probability of plaque appearance was still relatively low. 【Conclusion】Prophages were widely distributed in multi-drug resistant E. coli and carried a variety of resistance genes and virulence genes. Mild phages had a high induction rate, and have the risk of horizontal transmission of resistance genes and virulence genes, and need to be strengthened and sustained monitor.

Key words: prophage, Escherichia coli, induction, distribution characteristics, resistance genes