CN113249507B - 一种病原耐药基因的存在及表达情况的共检测方法 - Google Patents

一种病原耐药基因的存在及表达情况的共检测方法 Download PDF

Info

Publication number
CN113249507B
CN113249507B CN202110754521.6A CN202110754521A CN113249507B CN 113249507 B CN113249507 B CN 113249507B CN 202110754521 A CN202110754521 A CN 202110754521A CN 113249507 B CN113249507 B CN 113249507B
Authority
CN
China
Prior art keywords
sequencing
dna
detection method
drug
library
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
Application number
CN202110754521.6A
Other languages
English (en)
Other versions
CN113249507A (zh
Inventor
林东旭
袁光孝
朱方何
向亮
杨懿婧
陈杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Saizhe Zhizao Technology Co.,Ltd.
Original Assignee
Hunan Saizhe Zhizao Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan Saizhe Zhizao Technology Co ltd filed Critical Hunan Saizhe Zhizao Technology Co ltd
Priority to CN202110754521.6A priority Critical patent/CN113249507B/zh
Publication of CN113249507A publication Critical patent/CN113249507A/zh
Application granted granted Critical
Publication of CN113249507B publication Critical patent/CN113249507B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/6844—Nucleic acid amplification reactions
    • C12Q1/686—Polymerase chain reaction [PCR]
    • 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

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (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)
  • Chemical Kinetics & Catalysis (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明提供了一种病原耐药基因的存在及表达情况的共检测方法。所述共检测方法使用的反转录随机引物5'端包含固定碱基和随机碱基，用于区分DNA和RNA；所述共检测方法包括：样本核酸的提取、反转录、片段化处理和连接测序接头等文库构建步骤，再对所得文库进行宏基因组测序或捕获测序，进而得到病原耐药基因的存在情况和表达情况。所述共检测方法无论是通过宏基因组测序还是捕获测序，都可以通过单次实验测序得到耐药相关基因在DNA和RNA两个层面的测序数据，对指导临床用药具有更加重要的意义。

Description

一种病原耐药基因的存在及表达情况的共检测方法
技术领域
本发明属于微生物检测技术领域，具体涉及一种病原耐药基因的存在及表达情况的共检测方法，尤其涉及一种基于二代测序平台共检测耐药基因的存在及表达情况的方法。
背景技术
抗菌药物（抗生素）通过杀灭细菌发挥治疗感染的作用，包括β-内酰胺类、氨基糖苷类、四环素类、氟喹诺酮类、叶酸途径抑制剂类、多肽类和大环内酯类等，在治疗时需根据不同的感染性疾病进行合理选择。细菌作为一类广泛存在的生物体，可以通过多种形式获得对抗菌药物的抵抗作用，逃避被杀灭的危险，这种抵抗作用被称为“细菌耐药”（DrugResistance），获得耐药能力的细菌就被称为“耐药菌”。耐药菌对人类健康构成重大威胁，其常出现在医疗机构内，医院科室中内、外科均以大肠埃希菌（Escherichia coli）检出率最高，而ICU最常见病原菌是鲍曼不动杆菌（Acinetobacter baumannii），且ICU耐药率普遍高于非ICU。同时，越来越多的同一菌株携带多种耐药基因，这也进一步警示人们要合理地使用抗菌药物。
目前，临床常用的针对耐药菌的耐药基因检测方法主要是两种：1、药敏鉴定，该方法虽然是耐药基因检测的金标准，但是存在需要培养，鉴定周期长，灵敏度不高等缺点；2、荧光PCR，相较于药敏鉴定，虽然不依赖培养，提高了灵敏度，但也存在通量不高限制，一次只能检测几个耐药位点等问题。
宏基因组测序技术（mNGS）作为一种无差别测序方法，也可以应用于临床病原抗生素耐药基因检测和分析。但是，目前mNGS在临床方面的反馈效果并不理想，主要原因包括基因碎片化过于严重导致漏检或者检测到的片段临床意义不明确等问题。
靶向测序（Target sequencing）是一种对感兴趣的基因或区域进行特异富集并测序的方法，包括多重扩增测序法和探针捕获测序法两种方式。使用靶向探针捕获抗生素耐药基因或位点的高通量测序法来检测耐药基因已成为一个有潜力的方法。该方法和多重扩增测序法相比有以下几个优势：1、捕获探针的设计难度低于多重扩增；2、探针捕获的覆盖区域大，可以一次性捕获检测更多的耐药相关基因；3、部分耐药基因存在序列相似度高的特点，特别是归属于同一耐药基因家族的不同成员，如果使用多重扩增测序法难度较大。
目前为止，尚未开发出一种既能检测耐药基因片段是否存在又能够同时检测耐药基因转录表达情况的测序方法。
发明内容
针对现有技术存在的不足，本发明的目的在于提供一种病原耐药基因的存在及表达情况的共检测方法。所述方法通过单次实验测序即可得到耐药相关基因在DNA和RNA两个层面的测序数据，既能够检测耐药基因片段是否存在、又能够同时检测耐药基因转录表达情况。
为达此目的，本发明采用以下技术方案：
第一方面，本发明提供一种用于检测病原耐药基因的存在及表达情况的反转录随机引物，所述反转录随机引物包括5~10个随机碱基和2~5个固定碱基，所述固定碱基包括ACT。
由于RNA需要进行反转录步骤操作生成cDNA，所以，在测序序列含有固定碱基标记的序列，其原始序列均为RNA。本发明中，使用固定碱基ACT进行标记，可以防止其与其他扩增得到的测序序列混淆，区分度较高，反转录结束后就可以与核酸中的DNA合并打断建库和测序。
本发明中，所述反转录随机引物中随机碱基的个数可以为5个、6个、7个、8个、9个或者10个，固定碱基的个数可以为2个、3个、4个或者5个。
同时，本发明中，固定碱基的选择不考虑三者完全相同的碱基，如AAA，CCC等，其会影响引物整体的GC含量，也会影响引物的碱基复杂度；且本发明中，选择的固定碱基为ACT，其反转录后共检测的效果最好。
本发明中，所述固定碱基可以直接位于5'末端，或者与5'末端间隔若干个随机碱基；其中，优选为间隔1个随机碱基。
优选地，所述反转录随机引物自5'端起依次包括1~2个随机碱基、3个固定碱基和5~8个随机碱基。
优选地，所述反转录随机引物的序列如SEQ ID NO.1所示，具体为：
5'-NACTNNNNNN-3'。
随机引物中，N表示随机碱基A、T、G或C中的任意一个。
第二方面，本发明还包括如第一方面所述的反转录随机引物在实时荧光定量PCR检测、宏基因组测序或捕获测序中的应用。
第三方面，本发明提供一种病原耐药基因的存在及表达情况的共检测方法，所述共检测方法包括如下步骤：
（1）提取样本中的核酸；
（2）对步骤（1）所得核酸进行片段化处理；
（3）使用第一方面所述的反转录随机引物和逆转录酶对步骤（2）中所得片段化核酸进行反转录，再对片段化的核酸进行末端修复和加A尾，而后连接测序接头，纯化；
（4）对步骤（3）所得纯化后的产物进行PCR文库扩增并加index，纯化得到DNA文库；
（5）对所得DNA文库进行宏基因组测序，得到病原耐药基因的存在情况和表达情况；
或者，将所得DNA文库和捕获探针混合进行文库杂交反应，纯化捕获片段，再对所述纯化后的捕获片段进行二代高通量测序，获得病原耐药基因的存在情况和表达情况。
需要注意的是，本发明中所述“样本”或“样品”是指需要进行检测分析的物质，其可以来自于个体（如人或动物，更具体地如全血、血浆、唾液、脑髓液、关节液、痰液、肺泡灌洗液、尿液、粪便、腹胸腹水或组织等），也可以是其它来源的，例如一些经处理或未经处理的实验室材料，或人工合成的测试品。应理解，对“样本”或“样品”的检测并非仅涉及诊断目的，还可涉及其它非诊断目的。
作为本发明优选的技术方案，步骤（5）所述捕获探针检测的基因包括vanA、KPC、NDM-1、IMP、IMI、VIM、IND、SME-1、embB、ethA、gyrA、katG、ahpC、panD、pncA、murM、rpoB、rpsA、CTX-M、rpsL、rrs、tetM、mecA、mecB、mecC、CAU-1、FEZ-1、CphA、Sfh-1、OXA-2、mcr-1、ACC或MOX中的任意一种或至少两种的组合。
本发明中上述检测基因对应耐受的药物包括碳青霉烯类、头孢菌素类、甲氧西林类、β-内酰胺类、氟醛诺酮类、多肽类、氨基糖苷类、大环内酯类、四环素类、多粘菌素类或结核相关药物等。
作为本发明优选的技术方案，步骤（5）所述捕获探针还包含对照探针，所述对照探针用于检测人内源基因，包括GAPDH、ATM、RAD51或β-actin中的任意一种或至少两种的组合。
优选地，步骤（5）所述捕获探针带有生物素标记，长度为80~120 bp，例如80 bp、85bp、90 bp、95 bp、100 bp、105 bp、110 bp、115 bp或120 bp等。
本发明中，所述捕获在目标基因的外显子上均匀分布；同时，本发明还对设计的所有探针序列进行比对去除冗余探针，形成bed文件。
此外，由于捕获测序在方法学上的优势，本发明可以在上述基因的基础上进一步增加捕获基因或者扩大区域范围，进而得到更多的测序信息。
作为本发明优选的技术方案，步骤（5）所述二代高通量测序的数据分析方法包括：
使用Fastp对低质量的数据进行过滤，使用BWA比对耐药基因数据库，使用samtools进行去重并实现DNA和RNA数据的分类。
本发明中，所述数据分析流程可以按照如下流程进行：
1）下机数据质控：使用Fastp( https://github.com/OpenGene/fastp)对低质量的数据进行过滤；
2）序列比对：使用BWA(http://bio-bwa.sourceforge.net/)比对耐药基因数据库；
3）数据处理：数据使用samtools (http://samtools.sourceforge.net/)和gencore (https://github.com/OpenGene/gencore)进行DNA和RNA数据sort和去除重复（duplicate）。
本发明中，所述共检测方法可以按照如下步骤进行：
（1）样本提取：对样本中的DNA和/或RNA进行提取，所述样本包括全血、血浆、唾液、脑髓液、关节液、痰液、肺泡灌洗液、尿液、粪便、腹胸腹水和组织等，也可包括经处理或未经处理的实验室材料或人工合成的测试品等；
（2）片段化方式可以是超声法或酶切法等；例如，本发明中可使用Bioruptor pico打断仪进行超声打断；
打断条件：4℃条件下，设置超声30 s，间隔30 s，4~6个循环，取出样品，振荡混匀，短暂离心，重复上述打断步骤8~10次；
（3）cDNA第一条链合成：使用反转录随机引物和逆转录酶对片段化的核酸进行反转录；
再对片段化的核酸进行末端修复和加碱基A，然后两端加上测序接头，并用磁珠纯化以除去未连接的接头序列；
（4）PCR文库扩增加index区分不同样本，纯化后并进行文库质控；构建好的文库通过Qubit3.0和Agilent 2100进行浓度检测和长度分布检测来进行质量评价；
随后，可以使用mNGS直接进行测序，例如利用illumina Nextseq550平台进行PE75测序，即可得到病原耐药基因的存在情况和表达情况。
若需要进行序列捕获则进行后续步骤；
（5）将DNA文库和探针进行文库杂交反应，杂交时间为0.5~5 h（例如0.5 h、1 h、1.5 h、2 h、2.5 h、3 h、3.5 h、4 h、4.5 h或5 h等），杂交后再PCR然后纯化捕获片段；
（6）将捕获的目标DNA文库进行二代高通量测序，对低质量的数据进行过滤，比对耐药基因数据库去除重复并结合人参照基因的测序深度，实现DNA和RNA数据的分类，进而得到病原耐药基因的存在情况和表达情况。
第四方面，本发明还提供一种如第三方面所述的共检测方法在非疾病的诊断和治疗目的的、病原微生物的鉴定和/或检测中的应用。
本发明所述的数值范围不仅包括上述列举的点值，还包括没有列举出的上述数值范围之间的任意的点值，限于篇幅及出于简明的考虑，本发明不再穷尽列举所述范围包括的具体点值。
与现有技术相比，本发明的有益效果为：
（1）本发明提供一种反转录随机引物，其在5'端增加固定碱基序列标记，能够并不影响反转录的进行，同时还能标记反转录后的RNA，达到区分耐药基因的DNA和表达的RNA的目的；该结构为RT-qPCR和宏基因组测序中DNA和RNA的区分提供了实验基础；
（2）本发明提供的共检测方法，配合反转录随机引物，无论是通过mNGS还是捕获测序，都可以通过单次实验测序得到耐药相关基因在DNA和RNA两个层面的测序数据，并根据所得数据进行分析；且结合本发明中实施例所得的数据可知，在多数情况下耐药基因DNA存在与否和药敏相关性较低，所以判断耐药基因是否转录和表达对指导临床用药更有意义；
（3）由于部分耐药基因存在序列相似度高的特点，特别是归属于同一耐药基因家族的不同成员，如果使用多重扩增测序法难度较大，但设计探针捕获可以有效避开；本发明还结合快速杂交，在保证质量的前提下，尽可能缩短监测时间，满足临床需求；同时，鉴于捕获测序的方法学优势，本发明还可以在所述的耐药基因基础上进一步增加捕获基因或扩大区域范围，得到更多测序信息，因此，本发明提供的共检测方法在病原微生物的鉴定或检测中具有广阔的应用前景。
附图说明
图1为实施例3中捕获的目的文库的质控结果图。
图2为实施例3中下机测序数据的分析处理流程图。
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案，但下述的实例仅仅是本发明的简易例子，并不代表或限制本发明的权利保护范围，本发明的保护范围以权利要求书为准。
以下实施例中，若无特殊说明，所用试剂及耗材均购自本领域常规试剂厂商；若无特殊说明，所用实验方法和技术手段均为本领域常规的方法和手段。
实施例1
本实施例中使用反转录随机引物进行RT-qPCR检测，以验证反转录随机引物能够正常使用完成反转录步骤。
（1）以三份正常人新鲜血液作为待测样本；
（2）RNA提取：使用Qiagen QIAamp RNA Blood Mini Kit提取三份样本的全血RNA，其操作步骤参考试剂盒的使用说明书；
（3）cDNA第一条链合成：反转录中使用的反转录酶为Invitrogen公司的SuperScript IV逆转录酶，具体反应体系如下表1所示：
Figure 913250DEST_PATH_IMAGE001
其中，随机引物中N表示随机碱基A、T、G或C中的任意一个，随机引物的长度10 bp。
在冰上配制体系，混合均匀后短暂离心，然后在23℃，10 min；55℃，10 min；80℃，10 min进行反转录；
（4）qPCR扩增：针对人同源重组蛋白RAD51基因设计引物和探针，引物探针列表如下表2所示，每个样本三个重复，一共9个复孔。
Figure 361549DEST_PATH_IMAGE002
以上一步的反转录产物为模板在ABI7900上进行定量荧光PCR，反应体系和条件如下表3和表4所示：
Figure 173909DEST_PATH_IMAGE003
Figure 834698DEST_PATH_IMAGE004
（5）结果分析：
样本1的平均Ct值（Mean Ct）为27.8，样本2的平均Ct值（Mean Ct）为28.1，样本3的平均Ct值（Mean Ct）为27.9。
由实验结果可知，检测得到的Ct值正常且不同样本的结果接近，证明本发明设计的随机引物可以正常使用完成反转录实验。
实施例2
本实施例中以实施例1中提供的反转录随机引物配合mNGS检测RNA病毒及其耐药基因。
（1）样本选择：以5例明确RNA病毒阳性的肺泡灌洗液为样本；
（2）核酸提取和片段化：本实施例使用广州赛哲生物科技股份有限公司的DNA和RNA共提取试剂盒（粤穗械备20200771）进行核酸提取，
核酸片段化：取上一步处理的核酸，置于打断管中，使用Bioruptor pico打断仪进行打断；
打断条件：4℃条件下，设置超声30 s，间隔30 s，5个循环，取出样品，振荡混匀，短暂离心，重复上述打断步骤10次。
（3）反转录：cDNA第一条链合成参考实施例1中所述的步骤；
（4）末端修复：采用KAPA Hyper文库构建试剂盒，其中反应体系如下表5所示：
Figure 180229DEST_PATH_IMAGE005
在冰上配制上述反应体系，混合均匀后短暂离心，然后在20℃反应30 min；65℃反应30 min；
（5）接头连接：按照下表6配制体系；
Figure 268270DEST_PATH_IMAGE006
混合均匀，短暂离心后20℃反应15 min，4℃保存；
（6）连接产物纯化：提前半小时将Agencourt AMPure XP磁珠放在室温下平衡，使用0.6×磁珠对连接产物进行纯化，用22 μL无核酶水溶解，取20 μL到下一步反应中；
（7）文库PCR扩增：在冰上配制如下表7所示的反应体系；
Figure 800883DEST_PATH_IMAGE007
将配制好的体系混合均匀，短暂离心，然后运行如下表8所示的反应程序；
Figure 530941DEST_PATH_IMAGE008
（8）文库纯化：使用1.0×AMPure XP磁珠对PCR产物进行纯化，用22 μL无核酶水溶解，吸取20 μL纯化后的文库到干净的1.5 mL离心管中。
（9）质控和测序：构建好的文库通过Qubit3.0和Agilent 2100进行浓度检测和长度分布检测来进行质量评价；纯化的文库利用illumina Nextseq550平台进行PE75测序。
（10）测序结果分析：所得结果如下表9所示，由表可知，mNGS宏基因组检测和临床结果一致，RNA病毒确诊病例中有耐药基因检出但多数不表达。
Figure 465399DEST_PATH_IMAGE009
实施例3
本实施例中以实施例1中提供的反转录随机引物配合捕获测序检测肺泡灌洗液样本中的耐药基因及其表达情况。
（1）样本选择：以4例抗生素耐受的肺泡灌洗液作为待测样本；
（2）杂交捕获前步骤同实施例2；
（3）杂交捕获测序：在1.5 mL LoBind离心管中将4个文库（500 ng/文库）混合成一个pool，充分拨弹混匀后离心将混合物收集到管底，然后用封口膜封住离心管口，并在封口膜上刺几个小孔，最后用抽真空浓缩仪抽干；
按照表10加入各组分到离心管底，吹吸、混匀后全部转移到200 μL PCR管中，95℃变性10 min；
Figure 989922DEST_PATH_IMAGE010
加入2 μL探针Mix后涡旋大约3 s混合均匀（最后体积21 μL），65℃经过0.5 h快速杂交捕获；同时将Dynabeads M270链酶亲和素T1磁珠从4℃冰箱取出，室温平衡30 min；
所述探针序列的序列由表11中SEQ ID NO.5~37所示；
Figure 9830DEST_PATH_IMAGE011
Figure 278001DEST_PATH_IMAGE012
Figure 66965DEST_PATH_IMAGE013
完成杂交反应后用T1磁珠吸附捕获文库并纯化。
将捕获的目的文库用通用引物进行8个循环扩增，用1.0×AMPure XP磁珠纯化，质控结果如图1所示。
如图1所示，所得文库的片段长度集中在300~500 bp，所得文库质量较为稳定。
（4）测序数据分析处理：
具体流程如图2所示：
S1）测序下机reads：纯化的文库利用illumina Miniseq平台进行PE75测序；
S2）低质量序列过滤：将测序下机reads使用Fastp ( https://github.com/OpenGene/fastp)对低质量的数据进行过滤；
S3）捕获基因序列比对，该步骤包括：
S4）使用BWA(http://bio-bwa.sourceforge.net/)比对耐药基因数据库，去除重复的耐药基因序列；和S5）结合人参照内源基因（即人源基因GAPDH、ATM、RAD51和β-actin）的测序深度；
S6）区分DNA和RNA reads，即进行DNA和RNA数据分类。
所得结果如下表12所示，其中包括总映射率（Total mapping rate，%）、目标上限率（Target cap rate，%）、20×覆盖率（20×coverage，%）、存在的DNA片段、表达的RNA和临床表现的耐药种类。
Figure 526503DEST_PATH_IMAGE014
由表可知，待测样本中测得的DNA片段明显多于RNA片段，说明在耐药菌中可能含有多种耐药基因，但是多数的耐药基因并未表达；而且，通过上述方法测得的表达的RNA，其对应的耐药基因基本与药敏试验结果一致。
综上所述，本发明提供的共检测方法，可以通过单次实验测序得到耐药相关DNA和RNA两个层面的测序数据进行分析，且从结果可知，多数情况下耐药基因DNA的存在与否和药敏相关性较低，需要依据其表达情况即RNA来判断其药敏类型，因此，本发明所述的共检测方法对指导临床用药具有重要的意义。
申请人声明，以上所述仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，所属技术领域的技术人员应该明了，任何属于本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，均落在本发明的保护范围和公开范围之内。
序列表
<110> 广州赛哲生物科技股份有限公司
<120> 一种病原耐药基因的存在及表达情况的共检测方法
<130> 20210618
<160> 37
<170> SIPOSequenceListing 1.0
<210> 1
<211> 10
<212> DNA
<213> 人工合成()
<220>
<221> misc_feature
<222> (1)..(1)
<223> n is a, c, g, or t
<220>
<221> misc_feature
<222> (5)..(10)
<223> n is a, c, g, or t
<400> 1
nactnnnnnn 10
<210> 2
<211> 18
<212> DNA
<213> 人工合成()
<400> 2
gaattagtcc ttacgcaa 18
<210> 3
<211> 21
<212> DNA
<213> 人工合成()
<400> 3
ttggcgggaa ttctgaaagc c 21
<210> 4
<211> 17
<212> DNA
<213> 人工合成()
<400> 4
cgtgaccctg ggcgaga 17
<210> 5
<211> 94
<212> DNA
<213> 人工合成()
<400> 5
atacgaaggt gaacatcata gacacgccag gacatatgga tttctttata tcgttcatta 60
tcagttttag atggggcaat tctactgatt tctg 94
<210> 6
<211> 91
<212> DNA
<213> 人工合成()
<400> 6
ttcttcacca cccagggacg ggtttatcgg gccaaggcct acgacttgcc cgaggcctcc 60
cggacggcgc gcgggcagca cgtggccaac c 91
<210> 7
<211> 100
<212> DNA
<213> 人工合成()
<400> 7
tctcttgcca ctttcgcgca tgcgcaagaa ggcacgctag aacgttctga ctggaggaag 60
tttttcagcg aatttcaagc caaaggcacg atagttgtgg 100
<210> 8
<211> 97
<212> DNA
<213> 人工合成()
<400> 8
gtttagcttt gtcggcattt aatgctcatg ctaacaaaag tgatgctgcg gcaaaacaaa 60
taaaaaaatt agaggaagac tttgatggga ggattgg 97
<210> 9
<211> 102
<212> DNA
<213> 人工合成()
<400> 9
ctatggtaat gctcttgcac gagcaaataa agaatatgtc cctgcatcaa catttaagat 60
gctaaatgct ttaatcgggc tagaaaatca taaagcaaca ac 102
<210> 10
<211> 106
<212> DNA
<213> 人工合成()
<400> 10
tgtactctcg ccggataaaa aaatgcacgg attacttgtt aaaaagaacc atgaatatga 60
aatcaaccat gttgatgtag cattttcagc tttgcatggc aagtca 106
<210> 11
<211> 100
<212> DNA
<213> 人工合成()
<400> 11
gtggaaagcg ctttagcggt gtgggatgag cccgcggcct atcagcttgt tggtggggtg 60
acggcctacc aaggcgacga cgggtagccg gcctgagagg 100
<210> 12
<211> 103
<212> DNA
<213> 人工合成()
<400> 12
gattaattcg atgcaacgcg aagaacctta cctgggtttg acatgcacag gacgcgtcta 60
gagataggcg ttcccttgtg gcctgtgtgc aggtggtgca tgg 103
<210> 13
<211> 106
<212> DNA
<213> 人工合成()
<400> 13
aatagtgatc ttgaggctaa tgttcccatg attaaagcaa gtataaaaaa actaggtttt 60
aaattcagtg atactaaaat tttgctgatt agccatgctc attttg 106
<210> 14
<211> 102
<212> DNA
<213> 人工合成()
<400> 14
ccatactctc cccggatagg aaaacgcatg ggctgcttgt catgaaagaa agcgaatacg 60
aaacacggcg tattgatgtg gctttcccgg ttttgcatgg ca 102
<210> 15
<211> 101
<212> DNA
<213> 人工合成()
<400> 15
tcgggctgtg aggtcgggtg tgcggtcatg gggaacgagg atgatttgat tgtcggcgaa 60
gtggatcaaa tccggctgag ccacggtatc ttccgcatcc a 101
<210> 16
<211> 100
<212> DNA
<213> 人工合成()
<400> 16
accgcatcgt catggatcgg gaggtcgtca agctgggcgg catcgcgttc accgcccatt 60
tcatgcccgg gcacacgccg gggagcaccg cctggacctg 100
<210> 17
<211> 97
<212> DNA
<213> 人工合成()
<400> 17
cgcgcacacc gatcacgccg gccctgtggc cgagatcaag cgccgcacgg gggcgcaggt 60
ggtcgtgaac gccgaaaccg cagtgctgct ggcccga 97
<210> 18
<211> 101
<212> DNA
<213> 人工合成()
<400> 18
gatgaacatc ggccagattt tggagaccca cctgggttgg tgtgcccaca gcggctggaa 60
ggtcgacgcc gccaaggggg ttccggactg ggccgccagg c 101
<210> 19
<211> 104
<212> DNA
<213> 人工合成()
<400> 19
ggacatgacg tacgcggctc cactgttcgt caccgccgag ttcatcaaca acaacaccgg 60
tgagatcaag agtcagacgg tgttcatggg tgacttcccg atga 104
<210> 20
<211> 102
<212> DNA
<213> 人工合成()
<400> 20
cacagaaagc gcctgagtaa gcgtttgctg ctgcccatgt tgctgtgttt attggctcaa 60
caaacgcagg ctgtggcagc tgagcagacc aaggtcagtg ac 102
<210> 21
<211> 104
<212> DNA
<213> 人工合成()
<400> 21
gatgaaacca gtcgttttag ctgggacgga aaggttcgcg aaattgccgt ctggaacagg 60
gaccagagtt ttcgcaccgc aatgaagtac tctgtggtgc ctgt 104
<210> 22
<211> 104
<212> DNA
<213> 人工合成()
<400> 22
ataagtattg gggcttgtgc tgagcatagt atgagtcgag caaaaacaag tacaattcca 60
caagtgaata actcaatcat cgatcagaat gttcaagcgc tttt 104
<210> 23
<211> 96
<212> DNA
<213> 人工合成()
<400> 23
ccatgcaagc atctacagtg cctgtatatc aagaattggc acgtcgtatt ggtccaagct 60
taatgcaaag tgaattgcaa cgtattggtt atggca 96
<210> 24
<211> 102
<212> DNA
<213> 人工合成()
<400> 24
tcaagatact gaacacgctg atctcgctcg aagaaaaagc catctctgga aaagacgatg 60
tgctgaaatg ggacgggcat atttacgatt ttccagattg ga 102
<210> 25
<211> 104
<212> DNA
<213> 人工合成()
<400> 25
caagatactg aacacgctga tctcgctcga agaaaaagcc atctctggaa aagacgatgt 60
gctgaaatgg gacgggcata tttacgattt tccagattgg aatc 104
<210> 26
<211> 104
<212> DNA
<213> 人工合成()
<400> 26
tgcctcgacg ctactgatgg tatcttttag ttcattcgaa acctcagcgc aatccccact 60
gttaaaagag caaattgaat ccatagtcat tggaaaaaaa gcca 104
<210> 27
<211> 101
<212> DNA
<213> 人工合成()
<400> 27
tagccgtcaa cgacataggc tctagcgagg actttctcgc cgcaatagac aaaacgatca 60
agtacttcaa cgatggcgac atcgtcgaag gcaccatcgt c 101
<210> 28
<211> 100
<212> DNA
<213> 人工合成()
<400> 28
agggtgtcgt gtcctcgatc gtcaacttcg gcgcgttcgt cgatctcggc ggtgtggacg 60
gtctggtgca tgtctccgag ctatcgtgga agcacatcga 100
<210> 29
<211> 105
<212> DNA
<213> 人工合成()
<400> 29
agcagagtct ttgccagatt taaaaattga aaagcttgat gaaggcgttt atgttcatac 60
ttcgtttgaa gaagttaacg ggtggggcgt tgttcctaaa catgg 105
<210> 30
<211> 106
<212> DNA
<213> 人工合成()
<400> 30
ttggcgaagt aagaatgcag aaatacctta aaaaattttc ctatggcaac cagaatatca 60
gtggtggcat tgacaaattc tggttggaag gccagcttag aatttc 106
<210> 31
<211> 105
<212> DNA
<213> 人工合成()
<400> 31
gaaaaccagc taatagtaaa agaggctttg gtaacggagg cggcacctga atatctagtg 60
cattcaaaaa ctggtttttc tggtgtggga actgagtcaa atcct 105
<210> 32
<211> 107
<212> DNA
<213> 人工合成()
<400> 32
ctctttgatg ttccatggga aaaaatacag taccaaagcc tcatggatac cattaaaaaa 60
cgtcataatt taccggttgt agcggtattt gccacacact cccatga 107
<210> 33
<211> 105
<212> DNA
<213> 人工合成()
<400> 33
gaatttcggc cacggagcaa atcagctttt taagaaagct gtatcacaat aagttacacg 60
tatcggagcg cagccagcgt attgtcaaac aagccatgct gaccg 105
<210> 34
<211> 105
<212> DNA
<213> 人工合成()
<400> 34
cgatgaaata ttcagttgtg cctgtttatc aagaatttgc ccgccaaatt ggcgaggcac 60
gtatgagcaa gatgctacat gctttcgatt atggtaatga ggaca 105
<210> 35
<211> 106
<212> DNA
<213> 人工合成()
<400> 35
tcgatggtgg tattcgcatt tcggctaccc agcaaatcgc ttttttacgc aagctgtatc 60
acaacaagct gcacgtttct gagcgtagtc agcgcatcgt gaaaca 106
<210> 36
<211> 106
<212> DNA
<213> 人工合成()
<400> 36
ccgcgatgaa gtactcagtt gtgcctgttt atcaagaatt tgcccgccaa attggtgagg 60
cacgtatgag taaaatgctg cacgccttcg attatggcaa tgagga 106
<210> 37
<211> 105
<212> DNA
<213> 人工合成()
<400> 37
gcatctttgt ctacttcccc gaccagcagg tgctctatgg caactgcatc ctcaaggaga 60
agctgggcaa cctgagcttt gccgatgtga aggcctatcc gcgca 105

Claims (6)

1.一种非诊断和治疗目的的病原耐药基因的存在及表达情况的共检测方法，其特征在于，所述共检测方法包括如下步骤：
（1）提取样本中的核酸，包括DNA和RNA；
（2）对步骤（1）所得核酸进行片段化处理；
（3）使用反转录随机引物和逆转录酶对步骤（2）中所得片段化核酸进行反转录，再对片段化的核酸进行末端修复和加A尾，而后连接测序接头，纯化；
所述反转录随机引物自5'端起依次包括1个随机碱基、3个固定碱基和6个随机碱基，所述反转录随机引物的序列为：5'-NACTNNNNNN-3'；
（4）对步骤（3）所得纯化后的产物进行PCR文库扩增并加index，纯化得到核酸文库；
（5）对所得核酸文库进行宏基因组测序，得到病原耐药基因的存在情况和表达情况；
或者，将所得核酸文库和捕获探针混合进行文库杂交反应，纯化捕获片段，再对所述纯化后的捕获片段进行二代高通量测序，获得病原耐药基因的存在情况和表达情况；
步骤（5）所述捕获探针检测的基因包括vanA、KPC、NDM-1、IMP、IMI、VIM、IND、SME-1、embB、ethA、gyrA、katG、ahpC、panD、pncA、murM、rpoB、rpsA、CTX-M、rpsL、rrs、tetM、mecA、mecB、mecC、CAU-1、FEZ-1、CphA、Sfh-1、OXA-2、mcr-1、ACC或MOX中的任意一种或至少两种的组合。
2.根据权利要求1所述的共检测方法，其特征在于，步骤（5）所述捕获探针还包含对照探针，所述对照探针用于检测人内源基因，包括GAPDH、ATM、RAD51或β-actin中的任意一种或至少两种的组合。
3.根据权利要求1所述的共检测方法，其特征在于，步骤（5）所述捕获探针带有生物素标记，长度为80~120 bp。
4.根据权利要求1所述的共检测方法，其特征在于，步骤（5）所述二代高通量测序的数据分析方法包括：
使用Fastp对低质量的数据进行过滤，使用BWA比对耐药基因数据库，使用samtools和gencore进行去重并实现DNA和RNA数据的分类。
5.根据权利要求1所述的共检测方法，其特征在于，步骤（5）所述文库杂交反应中杂交时间为0.5~5 h。
6.如权利要求1~5任一项所述的共检测方法在非疾病的诊断和治疗目的方面的病原微生物的鉴定或检测中的应用。
CN202110754521.6A 2021-07-05 2021-07-05 一种病原耐药基因的存在及表达情况的共检测方法 Active CN113249507B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110754521.6A CN113249507B (zh) 2021-07-05 2021-07-05 一种病原耐药基因的存在及表达情况的共检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110754521.6A CN113249507B (zh) 2021-07-05 2021-07-05 一种病原耐药基因的存在及表达情况的共检测方法

Publications (2)

Publication Number Publication Date
CN113249507A CN113249507A (zh) 2021-08-13
CN113249507B true CN113249507B (zh) 2021-12-10

Family

ID=77190675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110754521.6A Active CN113249507B (zh) 2021-07-05 2021-07-05 一种病原耐药基因的存在及表达情况的共检测方法

Country Status (1)

Country Link
CN (1) CN113249507B (zh)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102732517A (zh) * 2012-07-19 2012-10-17 光明乳业股份有限公司 一种pcr步移技术及其所使用的试剂盒
CN102757960A (zh) * 2012-07-30 2012-10-31 光明乳业股份有限公司 一种扩增未知序列的半随机pcr技术及其引物和试剂盒
CN105132409A (zh) * 2015-09-02 2015-12-09 中国科学院北京基因组研究所 RNA 5mC重亚硫酸盐测序的文库构建方法及其应用
CN105143525A (zh) * 2013-03-15 2015-12-09 纽亘技术公司 用于生成成链rna或dna文库的方法、组合物和试剂盒
CN105420371A (zh) * 2015-12-21 2016-03-23 张明 一种多病原及耐药基因检测方法
CN108334750A (zh) * 2018-04-19 2018-07-27 江苏先声医学诊断有限公司 一种宏基因组数据分析方法及系统
WO2020106630A1 (en) * 2018-11-19 2020-05-28 The Regents Of The University Of California Methods for detecting and sequencing a target nucleic acid
CN111471676A (zh) * 2020-03-13 2020-07-31 广州市达瑞生物技术股份有限公司 一种宏基因组二代测序的建库样本的制备方法
CN112226488A (zh) * 2020-10-13 2021-01-15 上海宝藤生物医药科技股份有限公司 一种基于总核酸和宏基因组学构建病原微生物测序文库的方法及试剂盒
CN112522792A (zh) * 2020-01-22 2021-03-19 微岩医学科技(北京)有限公司 一种rna测序文库的构建方法
CN112530519A (zh) * 2020-12-14 2021-03-19 广东美格基因科技有限公司 一种检测样本中微生物和耐药基因的方法和系统
CN112795654A (zh) * 2021-02-09 2021-05-14 上海睿璟生物科技有限公司 用于生物体融合基因检测与融合丰度定量的方法及试剂盒
CN112813196A (zh) * 2020-12-31 2021-05-18 广州微远基因科技有限公司 用于检测病原微生物的捕获探针组、方法、试剂盒及应用
CN112967753A (zh) * 2021-02-25 2021-06-15 美格医学检验所(广州)有限公司 一种基于纳米孔测序的病原微生物检测系统和方法
CN113046415A (zh) * 2021-03-18 2021-06-29 微岩医学科技(北京)有限公司 一种rna测序文库的构建方法及其应用

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106884037B (zh) * 2015-12-16 2020-11-03 博奥生物集团有限公司 一种检测细菌耐药基因的基因芯片试剂盒

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102732517A (zh) * 2012-07-19 2012-10-17 光明乳业股份有限公司 一种pcr步移技术及其所使用的试剂盒
CN102757960A (zh) * 2012-07-30 2012-10-31 光明乳业股份有限公司 一种扩增未知序列的半随机pcr技术及其引物和试剂盒
CN105143525A (zh) * 2013-03-15 2015-12-09 纽亘技术公司 用于生成成链rna或dna文库的方法、组合物和试剂盒
CN105132409A (zh) * 2015-09-02 2015-12-09 中国科学院北京基因组研究所 RNA 5mC重亚硫酸盐测序的文库构建方法及其应用
CN105420371A (zh) * 2015-12-21 2016-03-23 张明 一种多病原及耐药基因检测方法
CN108334750A (zh) * 2018-04-19 2018-07-27 江苏先声医学诊断有限公司 一种宏基因组数据分析方法及系统
WO2020106630A1 (en) * 2018-11-19 2020-05-28 The Regents Of The University Of California Methods for detecting and sequencing a target nucleic acid
CN112522792A (zh) * 2020-01-22 2021-03-19 微岩医学科技(北京)有限公司 一种rna测序文库的构建方法
CN111471676A (zh) * 2020-03-13 2020-07-31 广州市达瑞生物技术股份有限公司 一种宏基因组二代测序的建库样本的制备方法
CN112226488A (zh) * 2020-10-13 2021-01-15 上海宝藤生物医药科技股份有限公司 一种基于总核酸和宏基因组学构建病原微生物测序文库的方法及试剂盒
CN112530519A (zh) * 2020-12-14 2021-03-19 广东美格基因科技有限公司 一种检测样本中微生物和耐药基因的方法和系统
CN112813196A (zh) * 2020-12-31 2021-05-18 广州微远基因科技有限公司 用于检测病原微生物的捕获探针组、方法、试剂盒及应用
CN112795654A (zh) * 2021-02-09 2021-05-14 上海睿璟生物科技有限公司 用于生物体融合基因检测与融合丰度定量的方法及试剂盒
CN112967753A (zh) * 2021-02-25 2021-06-15 美格医学检验所(广州)有限公司 一种基于纳米孔测序的病原微生物检测系统和方法
CN113046415A (zh) * 2021-03-18 2021-06-29 微岩医学科技(北京)有限公司 一种rna测序文库的构建方法及其应用

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Metagenomic next-generation sequencing of rectal swabs for the surveillance of antimicrobial-resistant organisms on the Illumina Miseq and Oxford MinION platforms;Rebecca Yee等;《European Journal of Clinical Microbiology & Infectious Diseases》;20200811;第40卷(第1期);第95-102页 *
宏基因组技术在感染性疾病中的应用;孙波等;《中国临床新医学》;20210130;第14卷(第1期);第19-22页 *

Also Published As

Publication number Publication date
CN113249507A (zh) 2021-08-13

Similar Documents

Publication Publication Date Title
KR102535489B1 (ko) 미생물 검출 방법 및 특성 규명 방법
CN106906211B (zh) 一种分子接头及其应用
CN109762915B (zh) 一种细菌耐药基因的检测方法及其专用试剂盒
CN110964840A (zh) 用于检测5种血流感染病原体的引物组、试剂盒及建库方法
CN111621466A (zh) 一种肺动脉组织单细胞悬液的制备方法及应用
JP2022137019A (ja) Ｍｉｔｒａチップ抽出を使用して嚢胞性線維症の変異を検出するための方法
WO2021179469A1 (zh) 检测病原体的组合物、试剂盒及方法
CN115927565A (zh) 基于mNGS检测CAR-T细胞制品中病原微生物的方法建立及应用
CN111662959A (zh) 一种真菌高通量快速检测方法
CN113249440B (zh) 基于纳米孔测序平台的组织样本宏基因组建库和检测方法
CN110777154A (zh) 用于结核分枝杆菌耐药检测的突变基因、及其检测方法和试剂盒
WO2021250617A1 (en) A rapid multiplex rpa based nanopore sequencing method for real-time detection and sequencing of multiple viral pathogens
CN113265452A (zh) 一种基于Nanopore宏基因组RNA-seq的生物信息学检测病原体的方法
CN113249507B (zh) 一种病原耐药基因的存在及表达情况的共检测方法
CN113249375A (zh) 一种快速高效富集粪便病毒的高通量检测方法
CN113201599A (zh) 一种基于PCR和nanopore测序检测脑脊液感染哪些病原体的方法
CN113817837A (zh) 一段来源于细粒棘球绦虫的游离dna序列及其应用
CN114107325B (zh) 宏基因组内参及其制备方法和应用以及宏基因组血流病原体检测方法
CN113549709A (zh) 利用巢式RPA技术检测SARS-CoV-2的引物对、探针、试剂盒及其应用
CN111254226A (zh) 用于快速检测犬瘟热病毒的引物组、荧光可视化rt-lamp试剂盒及方法
WO2024048236A1 (ja) ポリヌクレオチド、キット、及び、診断方法
TWI797593B (zh) 基於總體基因體定序之定序文庫的建構方法、病原檢測方法及其病原檢測系統(一)
CN112410343B (zh) 基于crispr的试剂盒及其应用
CN112941214B (zh) 一种用于革兰氏阴性菌耐药基因高通量扩增子测序的引物组及应用
US20230079748A1 (en) Preparation method, product, and application of circulating tumor dna reference samples

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20211027

Address after: 410000 1st, 6th and 7th floors, building A9, Hunan Inspection and testing characteristic industrial park, No. 336, Shixue Road, Shixue street, Yuelu District, Changsha City, Hunan Province

Applicant after: Hunan Saizhe Zhizao Technology Co.,Ltd.

Address before: 510000 units 304 and 305 on the third floor of R & D zone a and units 303, 304, 305 and 307 on the third floor of R & D zone B, No. 1, spiral Fourth Road, International Biological Island, Guangzhou, Guangdong

Applicant before: GUANGZHOU SAGENE BIOTECH Co.,Ltd.

GR01 Patent grant
GR01 Patent grant