CN105586417B - 一种玉米微生物耐药性基因的检测方法 - Google Patents
一种玉米微生物耐药性基因的检测方法 Download PDFInfo
- Publication number
- CN105586417B CN105586417B CN201610061033.6A CN201610061033A CN105586417B CN 105586417 B CN105586417 B CN 105586417B CN 201610061033 A CN201610061033 A CN 201610061033A CN 105586417 B CN105586417 B CN 105586417B
- Authority
- CN
- China
- Prior art keywords
- gene
- sample
- detected
- genes
- drug resistance
- 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.)
- Expired - Fee Related
Links
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 223
- 206010059866 Drug resistance Diseases 0.000 title claims abstract description 63
- 244000005700 microbiome Species 0.000 title claims abstract description 17
- 240000008042 Zea mays Species 0.000 title claims abstract description 16
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 title claims abstract description 16
- 235000002017 Zea mays subsp mays Nutrition 0.000 title claims abstract description 16
- 235000005822 corn Nutrition 0.000 title claims abstract description 16
- 238000001514 detection method Methods 0.000 title claims description 23
- 238000012163 sequencing technique Methods 0.000 claims abstract description 60
- 239000012634 fragment Substances 0.000 claims abstract description 55
- 230000003321 amplification Effects 0.000 claims abstract description 49
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 45
- 108020004707 nucleic acids Proteins 0.000 claims abstract description 30
- 102000039446 nucleic acids Human genes 0.000 claims abstract description 30
- 150000007523 nucleic acids Chemical class 0.000 claims abstract description 30
- 238000002474 experimental method Methods 0.000 claims abstract description 21
- 238000012165 high-throughput sequencing Methods 0.000 claims abstract description 17
- 239000003814 drug Substances 0.000 claims description 53
- 229940079593 drug Drugs 0.000 claims description 53
- 238000012360 testing method Methods 0.000 claims description 33
- 241001057636 Dracaena deremensis Species 0.000 claims description 6
- 206010048723 Multiple-drug resistance Diseases 0.000 abstract description 2
- 238000012258 culturing Methods 0.000 abstract 1
- 238000013461 design Methods 0.000 description 9
- 238000000605 extraction Methods 0.000 description 8
- 239000000203 mixture Substances 0.000 description 7
- 230000003115 biocidal effect Effects 0.000 description 6
- 239000003153 chemical reaction reagent Substances 0.000 description 6
- 108020004414 DNA Proteins 0.000 description 5
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 238000003752 polymerase chain reaction Methods 0.000 description 5
- 244000052616 bacterial pathogen Species 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 238000007403 mPCR Methods 0.000 description 4
- 239000003242 anti bacterial agent Substances 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 102000053602 DNA Human genes 0.000 description 2
- 102000012410 DNA Ligases Human genes 0.000 description 2
- 108010061982 DNA Ligases Proteins 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 229930193140 Neomycin Natural products 0.000 description 2
- 238000012408 PCR amplification Methods 0.000 description 2
- 229940088710 antibiotic agent Drugs 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001962 electrophoresis Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 229960000318 kanamycin Drugs 0.000 description 2
- 229930027917 kanamycin Natural products 0.000 description 2
- SBUJHOSQTJFQJX-NOAMYHISSA-N kanamycin Chemical compound O[C@@H]1[C@@H](O)[C@H](O)[C@@H](CN)O[C@@H]1O[C@H]1[C@H](O)[C@@H](O[C@@H]2[C@@H]([C@@H](N)[C@H](O)[C@@H](CO)O2)O)[C@H](N)C[C@@H]1N SBUJHOSQTJFQJX-NOAMYHISSA-N 0.000 description 2
- 229930182823 kanamycin A Natural products 0.000 description 2
- 239000011259 mixed solution Substances 0.000 description 2
- 229960004927 neomycin Drugs 0.000 description 2
- 108700022487 rRNA Genes Proteins 0.000 description 2
- 210000003705 ribosome Anatomy 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229920000936 Agarose Polymers 0.000 description 1
- 241000894006 Bacteria Species 0.000 description 1
- 108020004256 Beta-lactamase Proteins 0.000 description 1
- 102100038710 Capping protein-inhibiting regulator of actin dynamics Human genes 0.000 description 1
- 238000007400 DNA extraction Methods 0.000 description 1
- 108090000204 Dipeptidase 1 Proteins 0.000 description 1
- 101000957909 Homo sapiens Capping protein-inhibiting regulator of actin dynamics Proteins 0.000 description 1
- 229930182555 Penicillin Natural products 0.000 description 1
- JGSARLDLIJGVTE-MBNYWOFBSA-N Penicillin G Chemical compound N([C@H]1[C@H]2SC([C@@H](N2C1=O)C(O)=O)(C)C)C(=O)CC1=CC=CC=C1 JGSARLDLIJGVTE-MBNYWOFBSA-N 0.000 description 1
- 238000011529 RT qPCR Methods 0.000 description 1
- 238000003149 assay kit Methods 0.000 description 1
- 102000006635 beta-lactamase Human genes 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001647 drug administration Methods 0.000 description 1
- 239000000839 emulsion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012869 ethanol precipitation Methods 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 244000000010 microbial pathogen Species 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 229940049954 penicillin Drugs 0.000 description 1
- 230000026731 phosphorylation Effects 0.000 description 1
- 238000006366 phosphorylation reaction Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 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/6869—Methods for sequencing
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Analytical Chemistry (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
技术领域
本发明涉及生物技术领域，特别涉及一种玉米微生物耐药性基因的检测方法。
背景技术
抗生素类药物广泛应用于防治玉米病害，目前玉米病害在防治过程中，滥用抗生素类药物现象十分普遍，导致病原微生物进化出耐药性基因，产生耐药性，使得抗生素类药物在防治玉米病害时失效。因此，需要针对玉米中微生物耐药性基因针对性地用药，以避免病原菌产生耐药性，同时，减少农药用量，降低成本。
现有的检测耐药性基因的方法包括：预估待测样品上携带的需要检测的耐药性基因的病原菌，分离并纯化该病原菌，通过PCR(Polymease Chain Reaction，聚合酶链式反应)引物扩增该耐药性基因，利用电泳或实时定量PCR方法逐一判断待测样品中是否存在需要检测的耐药性基因。
在实现本发明的过程中，发明人发现现有技术至少存在以下问题中的一种：
需要培养、分离和纯化病原菌，使得检测的时间较长，延误治疗时机；无法检测不可培养的耐药菌中的耐药性基因；一次只能针对一种或少数几种病原菌进行检测；一次只能检测一种或少数几种耐药性基因，不能检测到待测样品中所有的耐药性基因，不能给予用药所需的全面信息；不能实现耐药性基因的定量检测。
发明内容
为了解决现有技术中检测耐药性基因的方法有待提高的问题，本发明实施例提供了一种玉米微生物耐药性基因的检测方法。所述技术方案如下：
本发明实施例提供了一种玉米微生物耐药性基因的检测方法，所述方法包括：
确定待测样品中需要检测的微生物的耐药性基因、所述待测样品中的内源标准基因和所述待测样品的外源标准基因，所述待测样品为整株玉米植株或部分所述玉米植株；
制备用于扩增所述耐药性基因、所述内源标准基因和所述外源标准基因的测试区域的多重扩增引物；
提取所述待测样品的基因组；
向所述待测样品的基因组中加入所述外源标准基因，获得混合核酸；
利用所述多重扩增引物对所述混合核酸进行扩增，得到扩增产物，利用所述扩增产物构建高通量测序文库；
对所述高通量测序文库进行高通量测序，得到测序片段组；
分析所述测序片段组，获得所述待测样品中所述耐药性基因的测序片段的数量、所述内源标准基因的测序片段的数量和所述外源标准基因的测序片段的数量；
根据所述外源标准基因的测序片段的数量和所述内源标准基因的测序片段的数量，判断实验是否成功；
若所述实验成功，则计算所述待测样品中所述耐药性基因的含量；
根据所述耐药性基因的含量判定所述待测样品是否含有耐药性基因。
具体地，所述耐药性基因至少为1个，所述内源标准基因至少为1个，所述外源标准基因至少为1个。
具体地，所述内源标准基因为所述待测样品的微生物中的基因。
具体地，所述外源标准基因不存在于已知基因组的生物中。
具体地，用于扩增所述耐药性基因的引物的设计区域或扩增区域与所述待测样品中的生物的基因组的其它区域的所有生物的基因组上除耐药性基因外的其它区域的同源性低于98％。
具体地，所述判断实验是否成功的方法为：当所述外源标准基因的测序片段的数量和所述内源标准基因的测序片段的数量均≥α1时，则实验成功；当所述外源标准基因的测序片段的数量或所述内源标准基因的测序片段的数量<α1时，则实验失败；其中，α1为判定阈值。
具体地，计算所述待测样品中所述耐药性基因的含量的方法为：第m种所述耐药性基因的含量的计算公式为
其中，i为第m种所述耐药性基因的第i个测试区域，n1为所述第m种所述耐药性基因的测试区域的个数，bi为所述第m种耐药性基因的第i个测试区域的测序片段的数量，k为第k种所述内源标准基因，n3为所述内源标准基因的个数，j为第k种所述内源标准基因的第j个测试区域，n2为第k种所述内源标准基因的测试区域的个数，aj为第k种所述内源标准基因的第j种测试区域的测序片段的数量；N为所有所述内源标准基因的测试区域的总数。
具体地，判定所述待测样品中是否含有所述耐药性基因的方法为：所述耐药性基因的含量≥α2时，判定所述待测样品含有所述耐药性基因；当所有所述耐药性基因的含量<α2时，判定所述待测样品不含有所述耐药性基因；其中，α2为判定阈值。
具体地，所述方法还包括：向所述待测样品中加入所述多重扩增引物不能扩增的外源核酸，将所述外源核酸与所述待测样品的基因组一同提取。
具体地，所述外源标准基因的质量与所述待测样品的基因组的总质量的比例大于1/100000。
本发明实施例提供的技术方案带来的有益效果是：本发明实施例提供的方法可以在不需要培养与纯化的前提下，一次性的、定量的、快速的、准确的检测任意数量和类型的微生物中的任意多种希望检测的耐药性基因，可以给予用药所需的全面信息。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将对本发明实施方式作进一步地详细描述。
实施例
玉米微生物的耐药性基因检测
待测样品可以为整株玉米植株或部分玉米植株，本实施例中待测样品为玉米叶片，本实施例中使用的玉米叶片取自武汉沌口开发区的玉米田，检测其中的微生物的耐药性基因是为了针对性地对玉米病害用药，避免耐药性的产生，同时，减少不必要的用药，节约成本。本实施例包括如下步骤：
步骤1、确定待测样品中需要检测的微生物的耐药性基因、待测样品中的内源标准基因和待测样品的外源标准基因，具体方法如下：
其中，耐药性基因至少为1个，内源标准基因至少为1个，外源标准基因至少为1个。内源标准基因为待测样品的微生物中的基因。外源标准基因不存在于已知基因组的生物中。
本实施例中的耐药性基因抗青霉素、β-内酰胺酶、头孢菌素酶、卡那霉素和新霉素，其中，耐药性基因APH(3')-Ia同时抗卡那霉素和新霉素。本实施例中内源标准基因为1种，具体地，内源标准基因为核糖体rRNA基因，该核糖体rRNA基因存在于绝大部分生物中且具有保守区域。在本实施例中，外源标准基因为1种，具体地，外源标准基因为ERCC-00004基因，其在NCBI(http://www.ncbi.nlm.nih.gov)上进行同源比对时，未发现其在已知生物参考基因组中存在同源序列，即外源标准基因不存在于已知的生物基因组中。表1为本实施例中被检测基因相关信息与检测结果，在表1中耐药性基因名称与序列编号与抗生素耐药性数据库(CRAD：The Comprehensive Antibiotic Resistance Datab ase，网址为http://arpcard.mcmaster.ca/)中的名称与序列编号一致。
表1 为本实施例中被检测基因相关信息与检测结果
表1中“/”表示无。
步骤2、制备用于扩增耐药性基因、内源标准基因和外源标准基因的测试区域的多重扩增引物，具体方法如下：
选择耐药性基因的保守区域设计多重扩增引物；用于扩增耐药性基因的引物的设计区域或扩增区域与待测样品中的生物的基因组的其它区域的同源性低于98％，这一要求确保了耐药性基因的检测不受待测样品中的生物的基因组的其它区域的干扰。具体地，耐药性基因的引物的设计区域在不同变体间是保守的，这样才能保证用相同的引物可以扩增出相同的耐药性基因的不同变体，若耐药性基因的引物的设计区域与待测样品中的生物的基因组的其它区域的同源性低于98％，则该引物不会扩增其它区域，因此，不会干扰耐药性基因的检测。否则，要求耐药性基因的引物的扩增区域与待测样品中的生物的基因组的其它区域的同源性低于98％，这样，获得的扩增产物可以与耐药性基因相区分，同样不会干扰耐药性基因的检测。同时，多重扩增引物中每对引物的扩增效率均在95％～105％之间。
具体地，在抗生素耐药性数据库中下载表1中的耐药性基因的序列，每一个下载的耐药性基因的序列的编号见表1。利用同源比对获得同一耐药性基因的不同序列间的保守区域，作为多重扩增引物的设计区域。若耐药性基因只包含一种序列，那么该序列的全部区域作为保守区域和多重扩增引物的设计区域。将以上获得的多重扩增引物的设计区域通过NCBI与所有生物的基因组上除耐药性基因外的其它区域进行同源比对，保留同源性低于98％的多重扩增引物的设计区域。核糖体rRNA基因的V5区以及外源标准基因ERCC-00004基因的全部序列作为保守区域和多重扩增引物的设计区域。
多重扩增引物获取过程如下：登录赛默飞世尔公司多重PCR引物在线设计网页https://ampliseq.com/，在“Application type”选项选择“DNA Hotspot designs(single-pool)”。将获得的多重扩增引物的设计区域用100个N连接起来，形成一个人工参考基因组，并在“Select the genome you wish to use”选项中选择“Custom”后，上传人工参考基因组。DNA Type选项选择“Standard DNA”。在多重扩增引物的设计区域内，随机挑选3个以上的不重叠的区域作为扩增区域，并填入Add Hotspot选项中，最后点击“Submittargets”按钮提交，获得用于扩增耐药性基因、内源标准基因和外源标准基因的多重扩增引物。由生工生物工程(上海)股份有限公司逐一合成多重扩增引物序列的每一重引物、以及由多重扩增引物的每一重引物所扩增的模板序列，模板序列是指上述多重扩增引物在填入Add Hotspot选项的扩增区域，按照美国赛默飞世尔公司的StepOne实时定量PCR仪的操作手册(Part Number 4376784Rev.E)检测每一多重扩增引物的扩增效率，仅保留扩增效率在95％～105％的多重扩增引物，在本实施例中，为节省成本，在满足条件的多重扩增引物中，仅选取最接近100％的一个引物作为最终入选的多重扩增引物，结果见表1。多重扩增引物由美国赛默飞世尔公司合成并以混合液体的形式提供使用。本实施例采用美国赛默飞世尔公司提供的多重PCR技术，其能够同时扩增多至12000个测试区域，因此，本发明有能力一次性检测现有的所有希望检测的耐药性基因。
步骤3、向待测样品中加入多重扩增引物不能扩增的外源核酸，将外源核酸与待测样品的基因组一同提取，得到混合样品，具体方法如下：
若待测样品中基因组含量较高，核酸提取正常，则不用添加外源核酸。若待测样品中基因组的含量较低，基因组提取困难，则需要向待测样品中加入多重扩增引物不能扩增的外源核酸后，可以方便提取待测样品中的基因组，一般情况下，需要加入1ug左右的外源核酸，即可保证待测样品中的基因组的正常提取。加入的外源核酸的量根据核酸的提取方法可以作适当调整，例如，大量提取时，加入量应有所增加，超微量提取时，可以适当地减少，在现有的基因组提取技术水平下，一般最少不得少于1ng。由于多重扩增引物不能扩增外源核酸，因此，所加入的外源核酸并不影响耐药性基因的检测。在本实施例中，外源核酸为ERCC-00024基因，外源核酸ERCC-00024基因的序列如序列表中SEQ ID NO:13所示，ERCC-00024基因由生工生物工程(上海)股份有限公司合成。本实施例中，在采用N96植物基因组DNA提取试剂盒提取基因组时，发现从待测样品中提取的核酸量正常，因此，不用向待测样品中加入的外源核酸。
步骤4、提取待测样品的基因组，具体方法如下：
按N96植物基因组DNA提取试剂盒(生产公司：天根生化科技(北京)有限公司，产品货号：DP338)的操作说明书提取混合样品的基因组，提取出的为混合样品的基因组核酸。利用分光光度计(美国Quawell公司生产，型号为Q5000)中的双链DNA程序，检测获得的待测样品的基因组的量。本实施例中，待测样品的基因组的总量为3520ng。
步骤5、向待测样品的基因组中加入外源标准基因，外源标准基因的质量与待测样品的基因组的总质量的比例大于1/100000，获得混合核酸，具体方法如下：
在本实施例中，外源标准基因的质量与待测样品的基因组的总质量的比例为1/1000，由于待测样品的基因组的总量为3520ng，因此，需要加入3.52ng外源标准基因，获得混合核酸。
步骤6、利用多重扩增引物对混合核酸进行扩增，得到扩增产物，利用扩增产物构建高通量测序文库，具体方法如下：
利用文库构建试剂盒2.0(由美国赛默飞世尔公司生产，货号为4475345)构建高通量测序文库。该文库构建试剂盒包括以下试剂：5×Ion AmpliSeqTM HiFi Mix、FuPa试剂、转换试剂、测序接头溶液和DNA连接酶。文库构建的方法按该文库构建试剂盒的操作手册《IonAmpliSeqTM Library Preparation》(出版号：MAN0006735，版本：A.0)进行。多重PCR的扩增体系如下：5×Ion AmpliSeqTM HiFi Mix 4μl、制备的多重扩增引物的混合液体4μl、混合核酸10ng和无酶水11μl。多重PCR的扩增程序如下：99℃，2分钟；(99℃，15秒；60℃，4分钟)×25个循环；10℃保温。利用FuPa试剂消化掉多重PCR扩增产物中多余的引物后，再进行磷酸化，具体方法为：向多重PCR的扩增产物中加入2μl FuPa试剂，混匀后，在PCR仪上按如下程序反应：50℃，10分钟；55℃，10分钟；60℃，10分钟；10℃保存，得到混合物a，混合物a为含有经过磷酸化的扩增产物溶液。将磷酸化的扩增产物连接上测序接头，具体方法为：向混合物a中加入转换试剂4μl、测序接头溶液2μl和DNA连接酶2μl，混匀后，在PCR仪上按如下程序反应：22℃，30分钟；72℃，10分钟；10℃保存，得到混合液b。利用标准的乙醇沉淀方法纯化混合液b后溶解于10μl无酶水中。利用美国Invitrigen公司生产的
dsDNA HS AssayKit(货号为Q32852)并按照其说明书进行检测，获得混合液b的质量浓度后，将纯化后的混合液b稀释至15ng/ml，得到浓度约100pM的高通量测序文库。
步骤7、对高通量测序文库进行高通量测序，得到测序片段组，具体方法如下：
利用获得的高通量测序文库和试剂盒Ion PI Template OT2 200 Kit v2(美国invirtrigen公司生产，货号为4485146)进行测序前的ePCR(Emulsion PCR，乳化聚合酶链反应)扩增，操作方法按该试剂盒的操作手册进行。利用ePCR产物和试剂盒Ion PISequencing 200 Kit v2(美国invirtrigen公司生产，货号为4485149)在Proton二代高通量测序仪上进行高通量测序，操作方法按该试剂盒的操作手册进行。在本实施例中，高通量测序量设置为1M测序片段(1M＝100万)，测序长度设置为500cycle(循环)，测序结束后，获得测序片段组。
步骤8、分析测序片段组，获得待测样品中耐药性基因的测序片段的数量、内源标准基因的测序片段的数量和外源标准基因的测序片段的数量，具体方法如下：
根据测序片段的引物，利用blastall(version 2.2.26)软件，按其默认的参数设置，将测序片段组比对到表1中的多重扩增引物对应的耐药性基因、内源标准基因和外源标准基因的检测区域上，保留与检测区域的同源性高于98％的测序片段(同源于低于98％的测序片段可能由非特异性扩增获得)，它们分别代表待测样品中耐药性基因的测序片段的数量、内源标准基因的测序片段的数量和外源标准基因的测序片段的数量，其结果见表1。
步骤9、根据外源标准基因的测序片段的数量和内源标准基因的测序片段的数量，判断实验是否成功，具体方法如下：
判断实验是否成功的方法为：当外源标准基因的测序片段的数量和内源标准基因的测序片段的数量均≥α1，则实验成功；当外源标准基因的测序片段的数量或内源标准基因的测序片段的数量<α1时，则实验失败；其中，α1为判定阈值。若外源标准基因的量过低，可能是高通量文库构建或高通量测序不成功，若内源标准基因量过低，可能是核酸提取不成功，若实验不成功则针对具体情况进行调整，直至实验成功为止。本实施例中，α1取10条测序片段，从表1可以看出，外源标准基因的测序片段的数量和内源标准基因的测序片段的数量均≥10条，因此，本实施例提供的实验是成功的。
步骤10、若实验成功，则计算待测样品中耐药性基因的含量，具体方法如下：
第m种耐药性基因的含量的计算公式为
其中，i为第m种耐药性基因的第i个测试区域，n1为第m种耐药性基因的测试区域的个数，bi为第m种耐药性基因的第i个测试区域的测序片段的数量，k为第k种内源标准基因，n3为内源标准基因的个数，j为第k种内源标准基因的第j个测试区域，n2为第k种内源标准基因的测试区域的个数，aj为第k种内源标准基因的第j种测试区域的测序片段的数量；N为所有内源标准基因的测试区域的总数。本发明采用平均数计算耐药性基因的含量，是为了减少扩增效率等因素对结果的影响，若检测了每一个多重扩增引物的扩增效率且保证扩增效率在95％-105％之间，那么扩增效率对结果的影响不大，因此，可以仅采用一个多重扩增引物扩增的检测区域来代表耐药性基因、内源标准基因和外源标准基因。
从表1可以看出，本实施例中，共检测了4种耐药性基因，所以m＝4，每一种耐药性基因检测了1个测试区域，所以，n1＝1，共检测了1个内源标准基因，所以n3＝1，每个内源标准基因共检测了1个测试区域，所以，n2＝1，N＝1，表1列出了每一种耐药性基因与每一个内源标准基因的测序片段的数量，将它们代入耐药性基因的含量的计算公式中可得到每一种耐药性基因的含量，其结果见表1。
步骤11、根据耐药性基因的含量判定待测样品是否含有耐药性基因，具体方法如下：
判定待测样品含有耐药性基因的方法为：耐药性基因的含量≥α2时，判定待测样品含有耐药性基因；当所有耐药性基因的含量<α2时，判定待测样品不含有耐药性基因；其中，α2为判定阈值。α2是根据要求的严格程度人为设定，当少量耐药性基因存在即可导致抗生素无效时，α2值应该偏低，相反，可偏高。当缺乏耐药性基因的量与抗生素效果关系时(大部分情况如此)，可将α2值设定为0.1％，以便于在统计误差与灵敏度间获得平衡。在本实施例中，将α2值设定为0.1％。从表1可以看出，本实施例子中，除耐药性基因CTX-M-14外，其它三种耐药性基因的含量均<α2＝0.1％，因此，判定待测样品为含有耐药性基因为CTX-M-14，其对应的抗生素类药物防治玉米病害时，由于耐药性的存在，可能效果不佳。本实施例通过4对多重扩增引物，一次性地定量检测了4种所耐药性基因的含量。在本实施例的基础上，只需要设置更多的多重扩增引物，即可实现一次定量检测任意多种希望检测的耐药性基因的含量。
验证本实施例的检测结果，具体方法如下：
人工合成表1中的多重扩增引物，逐一利用每一个多重扩增引物对待测样品的DNA进行PCR扩增，对扩增产物进行琼脂糖电泳检测，发现耐药性基因CTX-M-14存在扩增条带，判断为阳性，同时，内源标准基因和外源标准基因均有明显的扩增条带，判定为阳性，其它三种耐药性基因的引物的扩增产物无明显的条带，判定为阴性，该结果与本实施例的检测结果一致。
本发明实施例提供的方法可以检测不同的待测样品，通用性强。测序是核酸检测的终极标准，因此本发明实施例检测结果准确度高；测序可以区分单碱基差异，因此本发明实施例耐药性基因检测的分辨高。本发明实施例提供的方法可以在不需要培养与纯化的前提下进行检测，因此，检测是快速的，这对于病情的及时诊断和及时用药十分重要，同时，为不可培养病原物的检测提供了可能。本发明实施例可以检测任意多种耐药性基因，为寻找一种有效预防或治疗玉米病害的抗生素提供了全面信息。因此，本发明实施例可以一次性的、定量的、快速的、准确的检测任意数量和类型的微生物中的任意多种希望检测的耐药性基因，其效果是现有技术达不到的。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种玉米微生物耐药性基因的检测方法，其特征在于，所述方法包括：
确定待测样品中需要检测的微生物的耐药性基因、所述待测样品中的内源标准基因和所述待测样品的外源标准基因，所述待测样品为整株玉米植株或部分所述玉米植株；
制备用于扩增所述耐药性基因、所述内源标准基因和所述外源标准基因的测试区域的多重扩增引物；
向所述待测样品中加入所述多重扩增引物不能扩增的外源核酸，提取所述待测样品的基因组，将所述外源核酸与所述待测样品的基因组一同提取；
向所述待测样品的基因组中加入所述外源标准基因，获得混合核酸；
利用所述多重扩增引物对所述混合核酸进行扩增，得到扩增产物，利用所述扩增产物构建高通量测序文库；
对所述高通量测序文库进行高通量测序，得到测序片段组；
分析所述测序片段组，获得所述待测样品中所述耐药性基因的测序片段的数量、所述内源标准基因的测序片段的数量和所述外源标准基因的测序片段的数量，所述外源标准基因的质量与所述待测样品的基因组的总质量的比例大于1/100000；
根据所述外源标准基因的测序片段的数量和所述内源标准基因的测序片段的数量，判断实验是否成功，所述判断实验是否成功的方法为：当所述外源标准基因的测序片段的数量和所述内源标准基因的测序片段的数量均≥α1时，则实验成功；当所述外源标准基因的测序片段的数量或所述内源标准基因的测序片段的数量<α1时，则实验失败；其中，α1为判定阈值；
若所述实验成功，则计算所述待测样品中所述耐药性基因的含量，计算所述待测样品中所述耐药性基因的含量的方法为：第m种所述耐药性基因的含量的计算公式为
其中，i为第m种所述耐药性基因的第i个测试区域，n1为所述第m种所述耐药性基因的测试区域的个数，bi为所述第m种耐药性基因的第i个测试区域的测序片段的数量，k为第k种所述内源标准基因，n3为所述内源标准基因的个数，j为第k种所述内源标准基因的第j个测试区域，n2为第k种所述内源标准基因的测试区域的个数，aj为第k种所述内源标准基因的第j种测试区域的测序片段的数量；N为所有所述内源标准基因的测试区域的总数；
根据所述耐药性基因的含量判定所述待测样品是否含有耐药性基因，判定所述待测样品中是否含有所述耐药性基因的方法为：所述耐药性基因的含量≥α2时，判定所述待测样品含有所述耐药性基因；当所有所述耐药性基因的含量<α2时，判定所述待测样品不含有所述耐药性基因；其中，α2为判定阈值。
2.根据权利要求1所述的检测方法，其特征在于，所述耐药性基因至少为1个，所述内源标准基因至少为1个，所述外源标准基因至少为1个。
3.根据权利要求1所述的检测方法，其特征在于，所述内源标准基因为所述待测样品的微生物中的基因。
4.根据权利要求1所述的检测方法，其特征在于，所述外源标准基因不存在于已知基因组的生物中。
5.根据权利要求1所述的检测方法，其特征在于，用于扩增所述耐药性基因的引物的设计区域或扩增区域与所述待测样品中的生物的基因组的其它区域的同源性低于98％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610061033.6A CN105586417B (zh) | 2016-01-29 | 2016-01-29 | 一种玉米微生物耐药性基因的检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610061033.6A CN105586417B (zh) | 2016-01-29 | 2016-01-29 | 一种玉米微生物耐药性基因的检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105586417A CN105586417A (zh) | 2016-05-18 |
| CN105586417B true CN105586417B (zh) | 2020-01-31 |
Family
ID=55926316
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610061033.6A Expired - Fee Related CN105586417B (zh) | 2016-01-29 | 2016-01-29 | 一种玉米微生物耐药性基因的检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105586417B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102181575A (zh) * | 2011-03-18 | 2011-09-14 | 中国科学院武汉病毒研究所 | 用于检测乙型肝炎病毒耐药突变位点的引物及方法 |
| CN104561350A (zh) * | 2015-01-28 | 2015-04-29 | 深圳华大基因研究院 | 试剂盒及其用途 |
-
2016
- 2016-01-29 CN CN201610061033.6A patent/CN105586417B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102181575A (zh) * | 2011-03-18 | 2011-09-14 | 中国科学院武汉病毒研究所 | 用于检测乙型肝炎病毒耐药突变位点的引物及方法 |
| CN104561350A (zh) * | 2015-01-28 | 2015-04-29 | 深圳华大基因研究院 | 试剂盒及其用途 |
Non-Patent Citations (3)
| Title |
|---|
| Detection of cytomegalovirus drug resistance mutations by next-generation sequencing;Malaya K. Sahoo等;《Journal of Clinical Microbiology》;20131130;第51卷(第11期);3700–3710 * |
| 志贺菌1、2类整合子及ISCR1携带情况与耐药性的关系;董利娟等;《天津医药》;20150415;第63卷(第04期);400-403 * |
| 焦磷酸测序与PCR产物直接测序在乙型肝炎病毒耐药基因检测中的敏感性比较;孙树梅等;《中华医院感染学杂志》;20120720;第22卷(第14期);2991-2993 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105586417A (zh) | 2016-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Friedersdorf et al. | Advancing the functional utility of PAR-CLIP by quantifying background binding to mRNAs and lncRNAs | |
| CN110628890B (zh) | 测序质控标准品及其应用与产品 | |
| CN104120184B (zh) | 一种利用扩增dna片段长度多态性测定短链rna的方法 | |
| CN111808854B (zh) | 带有分子条码的平衡接头及快速构建转录组文库的方法 | |
| US20190390269A1 (en) | Method for detecting known nucleotide modifications in an rna | |
| DK3105324T3 (en) | NGS SYSTEM CONTROL AND PROCEDURES COMPREHENSIVE THESE | |
| CN104093854A (zh) | 表征组合物中的rna的方法和试剂盒 | |
| Hunter et al. | Broad variation in response of individual introns to splicing inhibitors in a humanized yeast strain | |
| US20210277458A1 (en) | Methods, systems, and aparatus for nucleic acid detection | |
| CN105603073B (zh) | 一种非诊断目的的痰液微生物耐药性基因的检测方法 | |
| CN110408612B (zh) | 一种低浓度dna标准物质的保护剂、保存方法及应用 | |
| CN105586417B (zh) | 一种玉米微生物耐药性基因的检测方法 | |
| CN116640833A (zh) | 一种利用毛细管电泳技术自动化批量检测mRNA加帽率的方法 | |
| CN105586410B (zh) | 一种土壤微生物耐药性基因的检测方法 | |
| CN105648060B (zh) | 一种非诊断目的的人体病原微生物耐药性基因的检测方法 | |
| CN105603080B (zh) | 一种非诊断目的的血液微生物耐药性基因的检测方法 | |
| CN105586415B (zh) | 一种非诊断目的的粪便微生物耐药性基因的检测方法 | |
| CN105586414B (zh) | 一种水稻微生物耐药性基因的检测方法 | |
| CN105624297B (zh) | 一种植物微生物耐药性基因的检测方法 | |
| US20240200154A1 (en) | Detection of sars-cov-2 variant | |
| Li et al. | Development of a novel reference plasmid for accurate quantification of genetically modified Kefeng6 rice DNA in food and feed samples | |
| CN113293202B (zh) | 一种定量检测机体中mRNA含量的实时荧光定量PCR试剂盒及检测方法和应用 | |
| CN111295713A (zh) | 稳定同位素标记的核酸作为内标物的核酸定量方法及其用途 | |
| CN104805186B (zh) | 一种测试玉米品种实质性派生关系的方法 | |
| CN116515976B (zh) | 一种转录组测序的校正方法及其试剂盒 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200131 |