一种检测幽门螺杆菌基因及药物代谢型的Miseq测序方法与流程

本发明涉及分子生物学检测领域，具体涉及一种基于miseq高通量测序技术同时检测幽门螺杆菌特异性基因、幽门螺杆菌耐药基因和宿主药物代谢酶cyp2c19基因的方法。
背景技术：
：幽门螺杆菌(helicobacterpylori)是一种螺旋形弯曲、微需氧的革兰氏阴性菌，被广泛证实为多种上消化道疾病的致病因子，是慢性胃炎、消化性溃疡、胃黏膜相关淋巴组织(malt)淋巴瘤的最主要病因，也是胃癌发生的重要诱因之一。临床上根除幽门螺杆菌的常用方法为两种抗生素加上一种质子泵抑制剂(ppi)的三联疗法或者加上铋剂的四联疗法。但是，目前幽门螺杆菌的根除率正在不断下降，其中细菌的耐药性和ppi的使用剂量不当是根除失败的主要原因。幽门螺杆菌对抗生素耐药的分子机制已被广泛研究，其中较为重要的两种抗生素克拉霉素和喹诺酮类药物的耐药机制已被研究的较为清晰。研究表明超过90％的克拉霉素耐药的幽门螺杆菌都有23srrna的2142和2143位点发生突变，常见的突变有a2142g、a2143g、a2143c三种，也有一些其他的突变，如t2182c、c2245t、t2289c等，但都是集中在23srrna核酸可变区v区。国内的一项研究也显示，超过80％的喹诺酮类药物耐药的幽门螺杆菌都有gyra基因发生突变，其中常见的突变有四种asn87lys、ala88val、asp91gly/asn/ala/tyr、ala97val，asp86asn较为少见，87和91突变在所有的耐药菌株中发现，并且与高耐药性相关，但都是集中在喹诺酮类药物耐药决定区(qrdr)发生突变。在分子水平上可以很好地检测幽门螺杆菌对这两种抗生素的耐药情况，因此目前迫切需要开发一种幽门螺杆菌根除个体化的治疗方案，明确药物代谢酶类型和幽门螺杆菌耐药性，针对性指导幽门螺杆菌感染用药。高通量测序技术，又称大规模平行测序技术，将dna(cdna)随机片段化后、加接头、制备文库测序，通过对文库中数以万计的克隆继续可逆延伸反应，荧光信号采集和信号转化，最终获得dna序列。通过扩增子建库结合miseq测序也已被广泛运用科学和临床研究，最典型的运用为16srrna测序进行细菌种属鉴定。技术实现要素：本发明的主要目的是提供一种准确、快速、检测流程简单的方法，可以同时检测幽门螺 杆菌特异性、耐药基因和宿主cyp2c19基因多态性位点，并能够更加科学、安全、快速、准确指导幽门螺杆菌感染患者的用药。为实现以上目的，设计覆盖ureb、23srrna、gyra、cyp2c19*2、cyp2c19*3基因待测区域的特异性引物，扩增产物在360-460bp。进一步的，将各正向引物的5’端加上tcgtcggcagcgtcagatgtgtataagagacag-特异引物，各反向引物5’端加上gtctcgtgggctcggagatgtgtataagagacag-特异引物。进一步的，将5对引物按照特定的比例合并，形成primerpool。进一步的，裂解胃黏膜组织样本，同时提取宿主和幽门螺杆菌基因组。进一步的，利用primerpool扩增出待测样本耐药区域的扩增产物，并进行磁珠纯化。进一步的，使用index引物对纯化的扩增子进行第二次扩增，每个样本带上不同的index，此时样本建库完成。进一步的，通过对文库的定量和质控，合格后进行miseq双端测序；对测序结果进行质控和分析，滤去接头等序列，比对耐药区域和药物代谢酶基因多态性位点，根据突变reads/未突变reads，计算出耐药比例和药物代谢酶类型。本方法只需要简单的对胃黏膜样本进行裂解处理，即可获得幽门螺杆菌和宿主的全基因组。设计出幽门螺杆菌的特异性、耐药和宿主药物代谢酶cyp2c19基因，利用多重pcr技术将这些待测片段在一个pcr反应中实现扩增，结合index引物的扩增实现待测样本的文库构建，一次测序即可获得幽门螺杆菌感染情况、耐药情况及宿主药物代谢酶cyp2c19基因多态性类型。该方法设计合理、操作流程简单、检测结果快速、准确、通量极高，可以满足不同科研和临床的检测需求，有利于指导临床关于幽门螺杆菌的检测和根除，为幽门螺杆菌感染治疗提供分子基础。此外，本发明的方法配合使用illumina测序平台的miseq测序仪，即可完成文库测序。此款测序仪相比较于其他测序仪具有更加简单、快速、经济、灵活的优势，整个建库流程1天即可，测序也只需要在2天内即可完成，相对于其他测序平台更加快捷，而且miseq测序的准确度目前来说是最准确的。附图说明图1是文库构建示意图。①第一步pcr扩增，靶基因的特异扩增和接头添加。x：接头和测序序列；y：特异引物序列；f：正向特异引物；r：反向特异引物。②第二步pcr完成测序文库构建。fc：与测序芯片flowcell配对的序列；b：不同样本的index序列；f：正向index引物；r：反向index引物。图2是agilent2100dna1000kit文库检测图具体实施方式为了使相关领域的科研人员和临床工作人员能够更好的理解本发明方案，下面结合实施方式和附图予以详细说明。但本发明并不限于以下实施方式。实施例1设计覆盖ureb、23srrna、gyra、cyp2c19*2、cyp2c19*3基因待测区域的特异性引物，扩增产物在360-460bp，具体见表1。表1基因待测区域的特异性引物在上述各正向引物的5’端加上tcgtcggcagcgtcagatgtgtataagagacag接头序列，各反向引物5’端加上gtctcgtgggctcggagatgtgtataagagacag接头序列，具体见表2。表2含接头序列的引物3.将上述引物分别稀释至10μm后进行混合，primerpool中这些引物的混合比例为：a-23s∶a-cyp2c19*2∶a-cyp2c19*3∶a-gyra∶a-ureb＝2∶0.5∶0.5∶1∶1。4.合成illumina公司的二次扩增index引物，稀释至10um。目前常用的引物为20条，其中n70x系列12条，s50y系列8条，总共可以有96种引物组合，一般可以满足样本需求，具体引物序列见表3。表3二次扩增index引物引物名称index引物序列(5’-3’)n701caagcagaagacggcatacgagattcgccttagtctcgtgggctcggn702caagcagaagacggcatacgagatctagtacggtctcgtgggctcggn703caagcagaagacggcatacgagatttctgcctgtctcgtgggctcggn704caagcagaagacggcatacgagatgctcaggagtctcgtgggctcggn705caagcagaagacggcatacgagataggagtccgtctcgtgggctcggn706caagcagaagacggcatacgagatcatgcctagtctcgtgggctcggn707caagcagaagacggcatacgagatgtagagaggtctcgtgggctcggn708caagcagaagacggcatacgagatcctctctggtctcgtgggctcggn709caagcagaagacggcatacgagatagcgtagcgtctcgtgggctcggn710caagcagaagacggcatacgagatcagcctcggtctcgtgggctcggn711caagcagaagacggcatacgagattgcctcttgtctcgtgggctcggn712caagcagaagacggcatacgagattcctctacgtctcgtgggctcggs501aatgatacggcgaccaccgagatctacactagatcgctcgtcggcagcgtcs502aatgatacggcgaccaccgagatctacacctctctattcgtcggcagcgtcs503aatgatacggcgaccaccgagatctacactatcctcttcgtcggcagcgtcs504aatgatacggcgaccaccgagatctacacagagtagatcgtcggcagcgtcs505aatgatacggcgaccaccgagatctacacgtaaggagtcgtcggcagcgtcs506aatgatacggcgaccaccgagatctacacactgcatatcgtcggcagcgtcs507aatgatacggcgaccaccgagatctacacaaggagtatcgtcggcagcgtcs508aatgatacggcgaccaccgagatctacacctaagccttcgtcggcagcgtc实施例2提取胃黏膜中的基因组(1)在无菌的2ml圆底ep管中，加入100μl组织裂解液和一个直径3mm的无菌钢珠；(2)将胃黏膜组织从样本保存液中取出，转移至含有裂解液和钢珠的2ml圆底ep管中；(3)组织研磨仪研磨样本，参数为室温条件下50hz，60s，直至无组织块；(4)充分研磨后，短暂离心，将匀浆液体转移至新的无菌1.5mlep管中；(5)将含有匀浆组织的1.5mlep管置于95℃金属浴中孵育10分钟，中间可颠倒数次；(6)迅速将ep置于4℃冰箱中孵育30分钟；(7)冰箱中取出ep管，12000转离心，5分钟，收集上清到新的ep管中，即可得到宿主和细菌的全基因组。实施例3目的片段的扩增和纯化(1)用按照比例混合好的pcr引物配制如下反应体系：反应总体系25μl，其中2*phusionhfbufferpcrmastermix12.5μl，模板3.5μl，混合引物primerpool2μl，ddh2o7μl。反应条件为95℃预变性3分钟，循环内95℃变性30秒，60℃退火30秒，72℃延伸30秒，20个循环；循环结束后72℃延伸5分钟；10℃保存。(2)在25μl扩增产物中加入20μlagencourtampurexpbeads磁珠，吹打混匀，室温放置5min；放入磁架上，2-5min后待液体变清，吸弃上清，缓慢加入新鲜配制的200μl80％乙醇洗涤，30s后吸弃乙醇，并重复一次；静置磁架5-8min挥发乙醇；待磁珠中心开始干裂即加入32.5μl水进行洗脱；从磁架上移下来，反复吹吸均匀后室温静置5min；放入磁架，2-5min后待液体变清，吸取30μl上清至新的洁净离心管中，即得到纯化后的pcr产物。文库建库参考示意图参考图1实施例4测序文库的构建、纯化、定量、质控(1)利用pcr的方法完成文库构建，反应总体系25μl，包括2*phusionhfbufferpcrmastermix12.5μl，纯化后的扩增产物4μl，引物n70x1μl，引物s50y1μl，ddh2o6.5μl，反应条件为95℃预变性3分钟，循环内95℃变性30秒，62℃退火30秒，72℃延伸30秒，15个循环；循环结束后72℃延伸5分钟，10℃保存。(2)取以上反应中25μl扩增产物加入25μlagencourtampurexpbeads磁珠，吹打混匀，室温放置5min；放入磁架上，2-5min后待液体变清，吸弃上清，缓慢加入新鲜配制的200μl80％乙醇洗涤，30s后吸弃乙醇，并重复一次；静置磁架5-8min挥发乙醇；待磁珠中心开始干裂即加入52.5μl水进行洗脱；从磁架上移下来，反复吹吸均匀后室温静置5min；放入磁架，2-5min后待液体变清，吸取50μl上清至新的洁净离心管中，即可得到已纯化的文库。(3)文库定量和质控质控，配套使用agilent2100bioanalyzer分析仪和qubit2.0浓度测定仪，qubit2.0检测文库浓度，agilent2100检测文库大小以及浓度。本方法使用agilentdna1000kit进行浓度和大小测定，具体如下：1)将dna1000kit试剂盒中的染料(蓝色盖子)和胶(红色盖子)冰箱中取出，进行室温避光平衡30min，；2)将染料混合均匀后，取25μl染料加入500μl胶中，涡旋10s混匀；3)转移染料-胶混合物至试剂盒中的离心柱，室温离心2240g，15min；4)吸取9μl染料-胶混合物加入至dna1000芯片的最右侧第三个孔(有标记)，避免在移液过程中产生气泡；5)挤压注射器，从1ml位置至dna位置的卡槽，60s后松开注射器卡槽，注射器会自然上升，5s后松开底座的dna芯片；6)分别吸取9μl染料-胶混合物至dna1000芯片最右侧第一、二孔(有标记)，避免气泡产生；7)再剩余的13个孔中，分别加入5μldnamarker(绿色盖子)；8)吸取1μldnaladder(黄色盖子)至最右侧第四个孔；9)将待测样本加入到芯片的前3列，每个样本吸取1μl；10)在摇床中，固定好dna芯片，2240rpm，混合60s；11)根据芯片类型，设定好agilent2100bioanalyzer参数，即可进行浓度的测定。agilent2100bioanalyzerdna1000kit文库质控检测图参考图2。实施例5文库测序、数据分析(1)质控合格后的文库进行测序前的稀释，具体如下：1)根据agilent2100bioanalyzer和qubit2.0的检测结果对文库进行稀释，此时用水稀释即可，稀释到2nm，充分混合均匀；2)用变性剂0.2n的naoh按1∶1稀释2nm文库，此时文库浓度为1nm，稀释时避免剧烈混合，用移液器吹吸均匀即可；3)变性剂稀释后室温放置5min；4)用测序试剂的杂交bufferht1对变性后的文库进行终浓度稀释，一般在10-15pm均可，充分混合均匀后，放置4℃冰箱保存，使用时取出。(2)文库稀释到合适的浓度后，即可按照每个样本的数据量需求对不同文库进行混合，本方法要求每个样本的reads数不低于5000条，测序试剂选择miseqreagentkitv2，pe500-cycles，具体测序操作参考illumina平台的miseq使用说明书，简单来说：先对miseq测序仪的管路进行清洗，测序试剂室温化冻，测序仪参数设置，测序。(3)测序数据下机后，miseq测序仪对测序结果进行自动转化，生成fastq格式的文件，通过比对每个样本的测序结果，分析突变情况，得到不同胃黏膜样本的幽门螺杆菌的感染情况、耐药情况以及宿主药物代谢酶类型，深度发掘数据，计算突变菌株和未突变菌株的比例，可得到混合感染中耐药和敏感菌株的比例。当前第1页12