CN115323067B - 用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 - Google Patents
用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 Download PDFInfo
- Publication number
- CN115323067B CN115323067B CN202211133701.3A CN202211133701A CN115323067B CN 115323067 B CN115323067 B CN 115323067B CN 202211133701 A CN202211133701 A CN 202211133701A CN 115323067 B CN115323067 B CN 115323067B
- Authority
- CN
- China
- Prior art keywords
- drug
- enterobacter cloacae
- phenotype
- characteristic
- kpc
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 82
- 239000003814 drug Substances 0.000 title claims abstract description 79
- 229940079593 drug Drugs 0.000 title claims abstract description 71
- 241000588697 Enterobacter cloacae Species 0.000 title claims abstract description 57
- 230000035945 sensitivity Effects 0.000 title claims abstract description 32
- 238000012163 sequencing technique Methods 0.000 title abstract description 12
- 230000003115 biocidal effect Effects 0.000 title description 10
- WKDDRNSBRWANNC-UHFFFAOYSA-N Thienamycin Natural products C1C(SCCN)=C(C(O)=O)N2C(=O)C(C(O)C)C21 WKDDRNSBRWANNC-UHFFFAOYSA-N 0.000 claims abstract description 15
- 229960002182 imipenem Drugs 0.000 claims abstract description 15
- ZSKVGTPCRGIANV-ZXFLCMHBSA-N imipenem Chemical compound C1C(SCC\N=C\N)=C(C(O)=O)N2C(=O)[C@H]([C@H](O)C)[C@H]21 ZSKVGTPCRGIANV-ZXFLCMHBSA-N 0.000 claims abstract description 15
- 229960002260 meropenem Drugs 0.000 claims abstract description 13
- DMJNNHOOLUXYBV-PQTSNVLCSA-N meropenem Chemical compound C=1([C@H](C)[C@@H]2[C@H](C(N2C=1C(O)=O)=O)[C@H](O)C)S[C@@H]1CN[C@H](C(=O)N(C)C)C1 DMJNNHOOLUXYBV-PQTSNVLCSA-N 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims description 47
- 101100026178 Caenorhabditis elegans egl-3 gene Proteins 0.000 claims description 8
- 101000740455 Klebsiella pneumoniae Metallo-beta-lactamase type 2 Proteins 0.000 claims description 8
- 101100108294 Caenorhabditis elegans aex-5 gene Proteins 0.000 claims description 7
- 101100058532 Caenorhabditis elegans bli-4 gene Proteins 0.000 claims description 7
- 101100187206 Enterobacter cloacae nmcR gene Proteins 0.000 claims description 7
- 101000702393 Homo sapiens Signal peptide peptidase-like 2B Proteins 0.000 claims description 5
- 102100030404 Signal peptide peptidase-like 2B Human genes 0.000 claims description 5
- 108010009816 IMP-8 enzyme Proteins 0.000 claims description 4
- 239000003153 chemical reaction reagent Substances 0.000 claims description 4
- 230000035772 mutation Effects 0.000 abstract description 21
- MYSWGUAQZAJSOK-UHFFFAOYSA-N ciprofloxacin Chemical compound C12=CC(N3CCNCC3)=C(F)C=C2C(=O)C(C(=O)O)=CN1C1CC1 MYSWGUAQZAJSOK-UHFFFAOYSA-N 0.000 abstract description 16
- GSDSWSVVBLHKDQ-JTQLQIEISA-N Levofloxacin Chemical compound C([C@@H](N1C2=C(C(C(C(O)=O)=C1)=O)C=C1F)C)OC2=C1N1CCN(C)CC1 GSDSWSVVBLHKDQ-JTQLQIEISA-N 0.000 abstract description 12
- 239000003242 anti bacterial agent Substances 0.000 abstract description 12
- 229940088710 antibiotic agent Drugs 0.000 abstract description 12
- 229960003376 levofloxacin Drugs 0.000 abstract description 12
- 229960003405 ciprofloxacin Drugs 0.000 abstract description 8
- 208000015181 infectious disease Diseases 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 206010059866 Drug resistance Diseases 0.000 description 24
- 238000000034 method Methods 0.000 description 7
- 230000001404 mediated effect Effects 0.000 description 6
- LISFMEBWQUVKPJ-UHFFFAOYSA-N quinolin-2-ol Chemical compound C1=CC=C2NC(=O)C=CC2=C1 LISFMEBWQUVKPJ-UHFFFAOYSA-N 0.000 description 6
- -1 KPC-6 Proteins 0.000 description 5
- LVZWSLJZHVFIQJ-UHFFFAOYSA-N Cyclopropane Chemical compound C1CC1 LVZWSLJZHVFIQJ-UHFFFAOYSA-N 0.000 description 4
- 239000013612 plasmid Substances 0.000 description 4
- 101100111650 Enterobacter cloacae nmcA gene Proteins 0.000 description 2
- 229940124350 antibacterial drug Drugs 0.000 description 2
- 230000001580 bacterial effect Effects 0.000 description 2
- 244000052616 bacterial pathogen Species 0.000 description 2
- 229960000484 ceftazidime Drugs 0.000 description 2
- NMVPEQXCMGEDNH-TZVUEUGBSA-N ceftazidime pentahydrate Chemical compound O.O.O.O.O.S([C@@H]1[C@@H](C(N1C=1C([O-])=O)=O)NC(=O)\C(=N/OC(C)(C)C(O)=O)C=2N=C(N)SC=2)CC=1C[N+]1=CC=CC=C1 NMVPEQXCMGEDNH-TZVUEUGBSA-N 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012165 high-throughput sequencing Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 230000008261 resistance mechanism Effects 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 238000012070 whole genome sequencing analysis Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000894006 Bacteria Species 0.000 description 1
- 206010011409 Cross infection Diseases 0.000 description 1
- 238000007400 DNA extraction Methods 0.000 description 1
- 108090000204 Dipeptidase 1 Proteins 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- YCKRFDGAMUMZLT-UHFFFAOYSA-N Fluorine atom Chemical compound [F] YCKRFDGAMUMZLT-UHFFFAOYSA-N 0.000 description 1
- 206010064571 Gene mutation Diseases 0.000 description 1
- 206010048723 Multiple-drug resistance Diseases 0.000 description 1
- 206010029803 Nosocomial infection Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 108010013381 Porins Proteins 0.000 description 1
- 206010040047 Sepsis Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 102000006635 beta-lactamase Human genes 0.000 description 1
- 210000000013 bile duct Anatomy 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 108010068385 carbapenemase Proteins 0.000 description 1
- 229940041011 carbapenems Drugs 0.000 description 1
- 229960002100 cefepime Drugs 0.000 description 1
- HVFLCNVBZFFHBT-ZKDACBOMSA-N cefepime Chemical compound S([C@@H]1[C@@H](C(N1C=1C([O-])=O)=O)NC(=O)\C(=N/OC)C=2N=C(N)SC=2)CC=1C[N+]1(C)CCCC1 HVFLCNVBZFFHBT-ZKDACBOMSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 125000001559 cyclopropyl group Chemical group [H]C1([H])C([H])([H])C1([H])* 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 229960003444 immunosuppressant agent Drugs 0.000 description 1
- 239000003018 immunosuppressive agent Substances 0.000 description 1
- 230000000968 intestinal effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 244000052769 pathogen Species 0.000 description 1
- 230000001717 pathogenic effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 102000007739 porin activity proteins Human genes 0.000 description 1
- 238000004902 predicting drug resistance Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 208000013223 septicemia Diseases 0.000 description 1
- 238000011895 specific detection Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 210000001635 urinary tract Anatomy 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
- C12Q1/689—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for bacteria
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Analytical Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法，通过检测特征组合中的基因的存在与否或突变的存在与否判断阴沟肠杆菌对亚胺培南、美罗培南、环丙沙星和左氧氟沙星的药敏表型。本发明提供的特征组合预测药敏表型的准确率高，有利于临床开展对阴沟肠杆菌相关感染的精准治疗。
Description
技术领域
本发明涉及基因测序技术，具体说，是一种用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法。
背景技术
阴沟肠杆菌广泛分布于自然界中，是肠道正常菌群之一，属于条件致病菌。当机体免疫功能低下或长期应用抗菌药物造成机体菌群失调时，易导致阴沟肠杆菌感染，引起呼吸道、伤口、胆管、泌尿道感染及败血症等，近年来已成为医院内感染的重要致病菌之一。
阴沟肠杆菌的耐药机制十分复杂，可同时产Ampc酶、超广谱β-内酰胺酶(ESBLs)及碳青霉烯酶，从而导致多重耐药，引起难治性感染。近年来，随着广谱抗菌药物、免疫抑制剂等大量而不合理的使用，阴沟肠杆菌的感染率不断增加，给临床治疗带来极大困扰，且耐药株的不断出现也引起了医学界的广泛关注。目前阴沟肠杆菌耐碳青霉烯类的耐药基因以KPC-2、NDM-1、IMP-4为主。阴沟肠杆菌对喹诺酮类药物的耐药机制主要包括由染色体介导的药物作用靶位改变、外膜孔蛋白通透性下降、外排泵主动外排和由质粒介导的耐药机制。染色体介导的喹诺酮类耐药只能垂直传播,由亲代传给子代,不易在不同菌株间发生播散,传播速度较慢。而临床资料表明,阴沟肠杆菌对喹诺酮类耐药的上升速度很快,可能存在由可移动元件如质粒介导的耐药基因的水平转移。近年来的研究发现,质粒介导的喹诺酮类耐药(plasmid mediated quinolone resistance,PMQR)由三类基因qnr、qepA和aac(6’)-Ib-cr参与,是细菌对喹诺酮类药物耐药的新机制。然而，针对不同抗生素，文献中已有报道的耐药相关基因多达数十至数百种，其中有些基因对耐药表型并不相关或贡献度小，对药敏表型判断造成了误导。研究证实，通过现有权威耐药基因数据库Resfinder和CARD中收录的耐药基因来预测耐药表型，准确率仅有66％和52％(doi:10.1093/jac/dkaa257)。
目前针对阴沟肠杆菌的耐药性检测可分为表型检测和基因型检测方法。在表型检测方面，临床的主要方法是微生物培养+药敏实验，该方法存在培养时间长，培养阳性率低等局限性，而无法完全满足临床精准治疗的需求。在基因检测方面，主要是针对特定、有限的耐药基因进行检测，该检测内容较单一，并不能反应真实的药敏结果。因此亟需寻找一种快速，不依赖于培养，且准确性高的检测临床常用抗生素药敏表型的方法以指导临床精准治疗。
发明内容
本发明所要解决的技术问题是，提供一种用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法，可一次性分析亚胺培南、美罗培南、环丙沙星、左氧氟沙星4种抗生素耐药情况，具有较好的准确率、敏感度和特异性。
为了解决上述技术问题，本发明采用的技术方案是：一种用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合，抗生素为亚胺培南、美罗培南、环丙沙星、左氧氟沙星中的一种或几种组合，
预测阴沟肠杆菌对亚胺培南药敏表型的特征基因组合包括NDM-1、NmcR、KPC-2、KPC-4、KPC-6、KPC-3、IMP-4、IMP-8、IMP-70和VIM-1，基因同时进行检测，若检测结果均为阴性，推测为敏感；若任一基因检测结果为阳性，推测为耐药；和/或
预测阴沟肠杆菌对美罗培南药敏表型的特征基因组合包括NDM-1、NmcR、NmcA、TEM-54、KPC-3、KPC-2，基因同时进行检测，若检测结果均为阴性，推测为敏感；若任一基因检测结果为阳性，推测为耐药；和/或
预测阴沟肠杆菌对环丙沙星药敏表型的特征基因突变位点和特征基因包括Ecl_parC:250(G->A)、Ecl_parC:239(G->T)和QnrB17，基因同时进行检测，若检测结果均为阴性，推测为敏感；若任一基因检测结果为阳性，推测为耐药；和/或
预测阴沟肠杆菌对左氧氟沙星药敏表型的特征基因突变位点和特征基因包括Ecl_gyrA:248(C->A)、Ecl_gyrA:248(C->T)和QnrB17，若检测结果均为阴性，推测为敏感；若检测结果为阳性，推测为耐药。
优选，所述预测阴沟肠杆菌对亚胺培南药敏表型的特征基因组合还包括ACT-45、TEM-54、OXA-1基因。
优选，所述预测阴沟肠杆菌对美罗培南药敏表型的特征基因组合还包括KPC-4、KPC-6、ACT-45基因。
优选，所述预测阴沟肠杆菌对环丙沙星药敏表型的特征基因突变位点还包括Ecl_gyrA:260(A->C)。
优选，所述预测阴沟肠杆菌对左氧氟沙星药敏表型的特征基因突变位点还包括Ecl_gyrA:260(A->C)。
含有上述用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合检测试剂的试剂盒。
采用上述试剂盒对阴沟肠杆菌进行药敏表型预测的测序方法，采用全基因组测序方法或宏基因组测序方法。
本发明的有益效果是：本发明是基于核酸分子检测耐药的方法，可不依赖于细菌培养，直接针对临床标本或其他方式取得的少量阴沟肠杆菌基因组核酸进行检测，并根据本发明提供的特征基因组合预测其药敏表型，相比传统方法具有检测周期短，检测灵敏度高的特点。
附图说明
图1是用于亚胺培南药敏表型预测的特征基因在公共数据库下载的阴沟肠杆菌菌株中的检出结果(左半)与实际药敏表型结果和预测药敏表型结果之间的对应关系(右半)图。
图2是用于美罗培南药敏表型预测的特征基因在公共数据库下载的阴沟肠杆菌菌株中的检出结果(左半)与实际药敏表型结果和预测药敏表型结果之间的对应关系(右半)图。
图3是用于环丙沙星药敏表型预测的特征基因在公共数据库下载的阴沟肠杆菌菌株中的检出结果(左半)与实际药敏表型结果和预测药敏表型结果之间的对应关系(右半)图。
图4是用于左氧氟沙星药敏表型预测的特征基因突变位点在公共数据库下载的阴沟肠杆菌菌株中的检出结果(左半)与实际药敏表型结果和预测药敏表型结果之间的对应关系(右半)图。
具体实施方式
下面将结合本发明实施例，对本发明实施例中的技术方案进行清楚、完整地描述；显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例，基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合，抗生素为亚胺培南、头孢他啶、头孢吡肟、环丙沙星中的一种或几种组合，
预测阴沟肠杆菌对亚胺培南药敏表型的特征基因组合包括NDM-1、NmcR、KPC-2、KPC-4、KPC-6、KPC-3、IMP-4、IMP-8、IMP-70和VIM-1，基因同时进行检测，若检测结果均为阴性，推测为敏感；若任一基因检测结果为阳性，推测为耐药；和/或
预测阴沟肠杆菌对美罗培南药敏表型的特征基因组合包括NDM-1、NmcR、NmcA、TEM-54、KPC-3、KPC-2，基因同时进行检测，若检测结果均为阴性，推测为敏感；若任一基因检测结果为阳性，推测为耐药；和/或
预测阴沟肠杆菌对环丙沙星药敏表型的特征基因突变位点和特征基因包括Ecl_parC:250(G->A)、Ecl_parC:239(G->T)和QnrB17，若检测结果均为阴性，推测为敏感；若任一基因检测结果为阳性，推测为耐药；和/或
预测阴沟肠杆菌对左氧氟沙星药敏表型的特征基因突变位点和特征基因包括Ecl_gyrA:248(C->A)、Ecl_gyrA:248(C->T)和QnrB17，若检测结果均为阴性，推测为敏感；若检测结果为阳性，推测为耐药。
优选，所述预测阴沟肠杆菌对亚胺培南药敏表型的特征基因组合还包括ACT-45、TEM-54、OXA-1基因。
优选，所述预测阴沟肠杆菌对美罗培南药敏表型的特征基因组合还包括KPC-4、KPC-6、ACT-45基因。
优选，所述预测阴沟肠杆菌对环丙沙星药敏表型的特征基因突变位点还包括Ecl_gyrA:260(A->C)。
优选，所述预测阴沟肠杆菌对左氧氟沙星药敏表型的特征基因突变位点还包括Ecl_gyrA:260(A->C)。
含有上述用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合检测试剂的试剂盒。
采用上述试剂盒对阴沟肠杆菌进行药敏表型的测序方法，采用全基因组测序方法或宏基因组测序方法。
下面结合具体的检测方法和效果对比，对本发明进行详细说明：
实施例1用特征组合对公共数据库中阴沟肠杆菌的药敏表型进行预测
1.1数据收集：从公共数据库(NCBI NDARO数据库和PATRIC数据库)下载817株阴沟肠杆菌基因组信息及其对应抗生素药敏表型数据。其中亚胺培南(英文缩写imipenem)耐药株99株，敏感株59株；美罗培南(英文缩写meropenem)耐药株77株，敏感株70株；环丙沙星(英文缩写ciprofloxacin)耐药株68株，敏感株279株；左氧氟沙星(英文缩写levofloxacin)耐药株65株，敏感株100株。
1.2耐药基因及突变检测：用ncbi-blast(v2.9.0+)软件将组装好的基因组序列与耐药数据库进行比对(参数：-evalue 1e-5 -outfmt 0 -num_alignments 10000)，进行上述耐药基因和突变位点检测。对于基因有无特征，与该耐药基因的参考序列的比对一致率高于90％且覆盖率高于60％，则认为有该耐药基因检出。对于基因突变特征，若比对上该耐药基因的序列支持该突变位点，则认为有该突变检出。
1.3统计各阴沟肠杆菌的菌株中耐药基因和突变的检出情况。
1.4药敏结果预测：对于某一株菌的任一种抗生素的药敏预测，检测出特征组合中的任一种特征，则认为该菌株对该抗生素表型耐药；否则，则判断为敏感。如图1-图4所示，亚胺培南、美罗培南、环丙沙星、左氧氟沙星四种模型中，通过各菌株中的检出的特征预测得到的药敏结果，和实际药敏表型高度一致。预测性能如表1所示，特征组合预测亚胺培南、美罗培南、环丙沙星、左氧氟沙星的药敏表型准确率分别为0.835、0.823、0.942、0.891，特异性和敏感度也在较高的水平。各耐药特征的特征检出样本数、表型耐药样本数和阳性预测值(阳性预测值＝表型耐药样本数/特征检出样本数)见表2。上述结果说明，利用上述特征组合，对阴沟肠杆菌表型耐药和表型敏感有较好的区分效果。
表1.预测公共数据库来源菌株的药敏表型的性能
表2.特征基因和突变在公共数据库来源的菌株中的检出频数
实施例2用优选的重要特征组合对公共数据库中阴沟肠杆菌的药敏表型进行预测2.1综合实施例1中各特征的检出频率、阳性预测值和特征间共同检出关系，选择NDM-1、NmcR、KPC-2、KPC-4、KPC-6、KPC-3、IMP-4、IMP-8、IMP-70和VIM-1基因有无作为亚胺培南的重要特征组合，选择OXA-23、ADC-30、ADC-162、ADC-56基因有无作为头孢他啶的重要特征组合，选择Ecl_parC:250(G->A)、Ecl_parC:239(G->T)和QnrB17基因突变位点和基因有无作为环丙沙星的重要特征组合，选择Ecl_gyrA:248(C->A)、Ecl_gyrA:248(C->T)和QnrB17基因突变位点和基因有无作为左氧氟沙星的重要特征组合。
2.2针对实施例1中下载的阴沟肠杆菌菌株，根据1.2得出的特征检出结果，对3.1所述优选特征进行统计。
2.3检测出优选特征组合中的任一种特征，则认为该菌株对该抗生素表型耐药；否则，则判断为敏感。如表3所示，利用优选的耐药特征预测药敏表型有着较好的准确率、灵敏度和特异性。可见，利用优选的特征组合也可以较好的预测公共数据库中下载的阴沟肠杆菌的药敏表型。
表3.利用优选特征组合预测公共数据库来源菌株的药敏表型的性能
实施例3用优选的重要特征组合对临床标本中的阴沟肠杆菌进行药敏表型预测
3.1样本采集：从某医院采集了培养结果为阴沟肠杆菌阳性，且进行了药敏测试的临床标本1例，包括肺泡灌洗液。
3.2样本高通量测序测序：对样本进行DNA提取、通过Qubit检测，确认DNA质量可以满足后续测序要求，对提取的DNA进行文库建库及高通量测序(Illumina Nextseq 550SE75)。
3.3测序数据指控：使用fastp(v0.19.5)软件对得到的原始fastq序列数据进行过滤(参数设置：-q 15-u 40-l read_length*0.67),去除低质量序列和过短序列；同时使用komplexity(v0.3.6)软件计算序列信息复杂度(参数设置：-F-t0.4)，并过滤掉低复杂度的序列。
3.4人源序列去除：将质控过滤得到的clean序列，使用bowtie2(v2.3.4.3)软件进行与人参考基因组序列(human_38)进行比对(参数设置：--mm--very-sensitive-k1)，以过滤掉人源序列。
3.5物种注释：使用minimap2软件(v2.17)将Illumina reads序列与以上设定的目标病原参考基因组序列库(来源于NCBI genome数据库)进行比对(比对参数：-xsr-a--secondary＝no-L))，并采用LCA算法进行物种注释统计，最后统计检出阴沟肠杆菌的序列数和基因组覆盖度。
3.6耐药特征检测：使用blastn(版本2.9.0+)软件将reads序列与耐药基因参考序列进行比对。对于基因有无特征，参考序列的比对一致率高于90％的reads大于1条，则认为有该耐药基因检出。对于耐药基因突变特征，针对某一突变位点，如果支持该突变位点的reads比例大于0.2，则认为有该突变检出。
3.7药敏表型预测：用2.1所述的各抗生素的重要特征组合对临床标本中阴沟肠杆菌对各抗生素的药敏表型进行预测。针对某一个抗生素，有任一重要特征检出，则认为阴沟肠杆菌对该抗生素耐药；否则，阴沟肠杆菌对该药敏感。当检出的基因组覆盖度小于一定比例且没有耐药特征检出时，由于不能确定该特征是否存在(特征可能存在于未覆盖的区域)，因此无法给出预测。表4为部分临床标本统计结果：
表4.部分临床标本的检测信息
以上所述的实施例仅用于说明本发明的技术思想及特点，其目的在于使本领域内的技术人员能够理解本发明的内容并据以实施，不能仅以本实施例来限定本发明的专利范围，即凡本发明所揭示的精神所作的同等变化或修饰，仍落在本发明的专利范围内。
Claims (2)
1.检测特征基因组合的试剂在制备用于预测阴沟肠杆菌对亚胺培南耐药性和敏感性的试剂盒中的用途，所述特征基因组合包括NDM-1、NmcR、KPC-2、KPC-3、KPC-4、KPC-6、ACT-45、TEM-54、IMP-4、IMP-8、IMP-70、VIM-1和OXA-1；同时对所述组合的所有特征基因进行检测，若检测结果均为阴性，推测为敏感；若任一特征基因检测结果为阳性，推测为耐药。
2.检测特征基因组合的试剂在制备用于预测阴沟肠杆菌对美罗培南耐药性和敏感性的试剂盒中的用途，所述特征基因组合包括NDM-1、NmcR、KPC-2、KPC-3、KPC-4、KPC-6、ACT-45、TEM-54和NmcA，同时对所述组合的所有特征基因进行检测，若检测结果均为阴性，推测为敏感；若任一特征基因检测结果为阳性，推测为耐药。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211133701.3A CN115323067B (zh) | 2022-09-19 | 2022-09-19 | 用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211133701.3A CN115323067B (zh) | 2022-09-19 | 2022-09-19 | 用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115323067A CN115323067A (zh) | 2022-11-11 |
| CN115323067B true CN115323067B (zh) | 2023-11-17 |
Family
ID=83930968
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211133701.3A Active CN115323067B (zh) | 2022-09-19 | 2022-09-19 | 用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115323067B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104560981A (zh) * | 2015-01-23 | 2015-04-29 | 杭州迪安医学检验中心有限公司 | 一种检测细菌喹诺酮类耐药基因的引物及试剂盒 |
| CN114067912A (zh) * | 2021-11-23 | 2022-02-18 | 天津金匙医学科技有限公司 | 基于机器学习筛选细菌耐药表型相关重要特征基因的方法 |
| CN114333987A (zh) * | 2021-12-30 | 2022-04-12 | 天津金匙医学科技有限公司 | 一种基于宏基因组测序的预测耐药表型的数据分析方法 |
-
2022
- 2022-09-19 CN CN202211133701.3A patent/CN115323067B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104560981A (zh) * | 2015-01-23 | 2015-04-29 | 杭州迪安医学检验中心有限公司 | 一种检测细菌喹诺酮类耐药基因的引物及试剂盒 |
| CN114067912A (zh) * | 2021-11-23 | 2022-02-18 | 天津金匙医学科技有限公司 | 基于机器学习筛选细菌耐药表型相关重要特征基因的方法 |
| CN114333987A (zh) * | 2021-12-30 | 2022-04-12 | 天津金匙医学科技有限公司 | 一种基于宏基因组测序的预测耐药表型的数据分析方法 |
Non-Patent Citations (4)
| Title |
|---|
| Deborah S. Ashcraft 等.1292. Evaluation of Synergy with Piperacillin/Tazobactam plus Meropenem Against Carbapenemase-Producing Klebsiella pneumoniae and Enterobacter cloacae Using ETEST.Open Forum Infect Dis.2021,第8卷(第Suppl 1期),S735. * |
| Negin Bolourchi等.Comparative resistome and virulome analysis of clinical NDM-1–producing carbapenem-resistant Enterobacter cloacae complex.Journal of Global Antimicrobial Resistance.2022,第28卷254–263. * |
| 冯婷婷等.阴沟肠杆菌感染与耐药机制的研究进展.中国人兽共患病学报.2017,第33卷(第10期),933-937. * |
| 王华等.2株产IMP型碳青霉烯酶阴沟肠杆菌的耐药特点及全基因组测序分析.临床与病理杂志.2022,第42卷(第6期),1275-1283. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115323067A (zh) | 2022-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Olm et al. | inStrain profiles population microdiversity from metagenomic data and sensitively detects shared microbial strains | |
| Ellington et al. | The role of whole genome sequencing in antimicrobial susceptibility testing of bacteria: report from the EUCAST Subcommittee | |
| Hashem et al. | Biofilm formation in enterococci: genotype-phenotype correlations and inhibition by vancomycin | |
| Lieberman et al. | Genetic variation of a bacterial pathogen within individuals with cystic fibrosis provides a record of selective pressures | |
| Gibreel et al. | High metabolic potential may contribute to the success of ST131 uropathogenic Escherichia coli | |
| Nwaiwu et al. | An in silico analysis of acquired antimicrobial resistance genes in Aeromonas plasmids | |
| CN114898808B (zh) | 一种预测肺炎克雷伯菌对头孢吡肟敏感性的方法及系统 | |
| Salamzade et al. | Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance | |
| Boeck et al. | Mycobacterium abscessus pathogenesis identified by phenogenomic analyses | |
| Tchesnokova et al. | A novel 7-single nucleotide polymorphism-based clonotyping test allows rapid prediction of antimicrobial susceptibility of extraintestinal Escherichia coli directly from urine specimens | |
| Lepuschitz et al. | Analytical performance validation of next-generation sequencing based clinical microbiology assays using a K-mer analysis workflow | |
| WO2017012660A1 (en) | Genetic testing for predicting resistance of serratia species against antimicrobial agents | |
| Hayer et al. | Global Distribution of Extended Spectrum Cephalosporin and Carbapenem Resistance and Associated Resistance Markers in Escherichia coli of Swine Origin–A Systematic Review and Meta-Analysis | |
| Wang et al. | High prevalence of unstable antibiotic heteroresistance in cyanobacteria causes resistance underestimation | |
| Lou et al. | Using strain-resolved analysis to identify contamination in metagenomics data | |
| Rad et al. | Detection of NDM-1 producing Klebsiella pneumoniae ST15 and ST147 in Iran during 2019–2020 | |
| Zhao et al. | Whole-genome sequencing reveals high-risk clones of Pseudomonas aeruginosa in Guangdong, China | |
| Humphries et al. | Machine-learning model for prediction of cefepime susceptibility in Escherichia coli from whole-genome sequencing data | |
| CN115651990B (zh) | 用于预测大肠埃希菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 | |
| Liu et al. | Molecular epidemiology and clinical characteristics of the type VI secretion system in Klebsiella pneumoniae causing abscesses | |
| Liu et al. | Three novel sequence types carbapenem-resistant Klebsiella pneumoniae strains ST5365, ST5587, ST5647 isolated from two tertiary teaching general hospitals in Shanxi Province, in North China: molecular characteristics, resistance and virulence factors | |
| CN115305292B (zh) | 用于预测金黄色葡萄球菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 | |
| CN115323067B (zh) | 用于预测阴沟肠杆菌对抗生素药敏表型的特征基因组合、试剂盒及测序方法 | |
| Chilam et al. | Genomic surveillance of Pseudomonas aeruginosa in the Philippines, 2013–2014 | |
| Adjei et al. | Phylogenetic analysis of carbapenem-resistant Acinetobacter baumannii isolated from different sources using Multilocus Sequence Typing Scheme |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address |
Address after: 100071 Room 110, 1st Floor, Building 3, No. 2, Dongbinhe Road, You'anmenwai, Fengtai District, Beijing Patentee after: Beijing Jinshao Medical Laboratory Co.,Ltd. Country or region after: China Patentee after: Tianjin JinKe Medical Technology Co.,Ltd. Patentee after: Tianjin Huazhinuo Technology Co.,Ltd. Address before: 100071 Room 110, 1st Floor, Building 3, No. 2, Dongbinhe Road, You'anmenwai, Fengtai District, Beijing Patentee before: Beijing Jinshao Medical Laboratory Co.,Ltd. Country or region before: China Patentee before: Tianjin JinKe Medical Technology Co.,Ltd. Patentee before: Jinshi Zhizao (Tianjin) Medical Technology Co.,Ltd. |