本发明涉及一种用于hpv分型及整合的高通量测序检测方法。
背景技术：
：人乳头瘤病毒（hpv）是一种小型dna病毒，是一种嗜黏膜和皮肤上皮性病毒，具有高度特异性。hpv病毒，主要由dna核心和蛋白衣壳组成。其基因组为双链环状dna，长约7.5～8.0kb，含8个开放阅读框（openreadingframes,orfs），依功能不同分3个区：（1）早期区（earlyregion,e区）：分别编码e1、e2、e4、e5、e6、e7等六个早期蛋白，参与病毒dna的复制、转录、翻译调控和转化等功能；（2）晚期区（lateregion,l区）：编码主要外壳蛋白l1和次要外壳蛋白l2；（3）非编码区（uncodingregion,ucr）又称长控区（lcr）或上游调控区(urr)：含有hpv基因组dna的复制起点和hpv表达所必需的调控元件。hpv的分类主要依据hpv病毒基因组中l1区的dna同源性。目前已经发现了170多种hpv型，有40多种可感生殖腔。hpv可以引起人全身各部位的皮肤或者粘膜感染，导致多种上皮组织疾病。例如皮肤寻常疣，扁平疣，外生殖器的尖锐湿疣等。宫颈上皮组织长期感染的某些型别的hpv可能引起基因突变，导致宫颈上皮内增生和宫颈癌。宫颈癌发病率居全球女性恶性肿瘤的第2位，发病率仅次于乳腺癌，每年约有50万宫颈癌新发病例，其中约1/3为死亡病例。据统计中国每年宫颈癌新发病例约13.5万，发病率为14.6/10万。根据其导致宫颈癌风险的高低将这些hpv分为高危型和低危型。高危型hpv可能引发妇女发生宫颈癌．主要包括hpv16、18、31、33、35、39、45、51、52、56、58、59、68、26、53、66、73、82等型别；低危型hpv主要引起尖锐湿疣，包括hpv6、11、40、42、43、44、54、61、81等型别。cp8304等亚型，与宫颈癌及宫颈上皮内高度病变（cinⅱ/ⅲ）的发生相关，尤其是hpv16和18型。hpv感染是宫颈癌发生的主要病因，持续高危型hpv的感染是癌前病变和浸润性宫颈癌发展的必要因素。由于不同亚型hpv的致病机制、致癌危险性及感染后的预后各不相同，因此对hpv各亚型进行分型以及是否整合到人类基因组上的检测有重要的临床意义。因此，临床上有必要建立一种简便、快速、灵敏、特异的分型及整合的检测方法。对hpv病毒感染进行准确的检测可提高宫颈癌前病变筛查的敏感性，同时改善妇女宫颈癌的防治。宫颈癌早期症状不明显，筛查预警对早期发现宫颈癌，降低死亡率极为重要。传统上，宫颈涂片检测是群体筛查的常规手段，属于细胞学检测方法。实验室的技术人员会于显微镜下观察不正常的细胞以诊断病人的病情，但由于这方法只以观察作为诊断基础，受取材、涂片制作质量、阅片技术影响，其准确率低，假阴性高，重复性差，得到的结果及分析漏诊率可达30%。在检查和筛选hpv的基因型方法中，实时荧光pcr法及杂交捕获法是最为常用的，但是对于不同种基因型hpv的诊断操作繁琐，检测通量低，难以满足临床上同时检测多个亚型的需要，不适合于对子宫颈癌的大规模筛查。商业化的杂交捕获法试剂盒无须pcr扩增即可检测出临床样品中的hpvdna，并且可区别高危型和低危型两类，但是此类试剂盒无法鉴定出感染hpv的具体基因型及是否整合，其假阴性和假阳性较多，灵敏性和特异性有待提高。但在分子诊断领域，最直接而明确的技术是基因测序，传统的基因测序方法是通过sanger终止法进行的，该方法一次只能读取单个基因的序列，而对于有众多基因亚型感染的疾病，很难获得理想的测序结果。现有的这些检测方法只能进行少量的hpv分型，而对于hpv是否整合到人类基因组上无法检测，而hpv是否整合到人类基因组上是发生hpv导致癌症的前提。第二代高通量测序测序技术实现了同时对几百万甚至数亿条dna序列进行测序，具有费用低、通量高、速度快、灵敏度高及准确度高等优点，更重要的是可以在检测分型的同时，检测hpv是否整合到人类基因组上，以及整合的位点是否会导致癌症的发生。本专利在第二代高通量测序测序技术的帮助下，可以对高危型hpv及低危型hpv进行精确分型，并对hpv是否整合到人类基因组上，以及整合的位点是否会导致癌症的发生进行精确地个体化评估，预防病变发生的风险，从而预防肿瘤的发生。技术实现要素：本发明提供一种准确度高、临床意义明确的hpv分型及整合检测方法。1.一种用于hpv分型及整合的检测方法，包括hpv分型及整合高通量测序测序检测方法。hpv分型及整合检测试剂盒的方法包括如下步骤：(1)根据hpv基因组，调取hpv基因组序列。常见hpv分型列表如下：hpv6hpv11hpv16hpv18hpv26hpv31hpv33hpv35hpv39hpv45hpv51hpv52hpv53hpv56hpv58hpv59hpv66hpv68hpv73hpv82hpv40hpv42hpv43hpv44hpv54hpv61hpv70hpv72hpv81调取包含列表不同亚型hpv基因组序列；(2)对每个区域中非重复区域设计120bp的探针序列，每个序列沿着基因位置挪动设计，探针之间的挪动大小60bp;(3)采用原位合成技术，在芯片上大量合成设计的探针，并利用多聚酶链式反应或转录的方法扩增出大量的带有生物术标记的探针，并制作hpv基因扫描试剂盒;所述试剂盒筛查方法包括如下步骤:从患者获取宫颈细胞样本并提取dna，然后利用hpv基因探针将目的基因序列捕获出来，再用测序仪(illuminanextseq500)进行高通量测序测序，进而分析，找出与hpv基因信息，从而得到患者hpv的基因型，以达到区分是否高危hpv及整合。根据权利要求1所述的一种用于hpv分型及整合的检测方法，采用测序仪进行高通量测序测序分析流程包括如下步骤:(1)测序仪(illuminanextseq500)获取原始短序列并去除测序数据中的接头;(2)用trimgalore进行qc并去除低质量数据;(3)把短序列用bwa-burrows-wheeleralignment(bwa)软件定位到相应的hpv基因组数据（例如alphapapillomavirus7,humanpapillomavirus18(hpv18),completegenomegi|60975|lcl|hpv18ref.1|hpv18ref）相应的位置上(所用到参数:bwaaln-l-l40-i10-k2-t7-e40-m3-f);(4)使用picardtools及自写脚本，针对不同样本不同类型的hpv基因进行归一化处理，统计相应的覆盖度及短序列数量，得到分型结果。所述基因整合分析流程包括如下步:(1)测序仪(illuminanextseq500)获取原始短序列;(2)用trimgalore进行qc并去除低质量数据;(3)把短序列用bwa-burrows-wheeleralignment(bwa)软件定位到hpv基因组数据相应的位置上(所用到参数:bwaaln-l-l40-i10-k2-t7-e40-m3-f);(4)用factera检测多个基因的全部或一部分的编码区首尾相连。包括染色体易位、中间缺失或染色体倒置所致的结果，得到基因整合结果。本发明一种hpv分型及整合的高通量测序检测方法的有益效果如下:本发明一种hpv分型及整合的高通量测序检测方法，是以dna文库的方式进行捕获测序，然后利用hpv基因探针将目的基因序列捕获出来，再用测序仪(illuminanextseq500)进行高通量测序测序，进而分析，找出与hpv基因信息，从而得到患者hpv的基因型，以达到区分是否高危hpv及整合。本方法采取高深度测序的方法（500x）检测hpv分型及整合。相比于临床上传统的病理学诊断而言，本方法有以下优势:1、相比于传统宫颈涂片检测等常规细胞学检测手段。采用二代高通量测序等直接测序鉴定hpv分型及整合，其准确率更高，且重复性更好；采取高深度测序，平均要达到500x，保证足够的检测深度，使假阴性更低。2、相比于传统的基因测序方法，一次只能读取单个基因的序列，本专利所用方法能高通量测序区分hpv分型及整合，一次检测即可区分下列hpv6，hpv11，hpv16，hpv18，hpv26，hpv31，hpv33，hpv35，hpv39，hpv45，hpv51，hpv52，hpv53，hpv56，hpv58，hpv59，hpv66，hpv68，hpv73，hpv82，hpv40，hpv42，hpv43，hpv44，hpv54，hpv61，hpv70，hpv72，hpv81。具体实施方式本发明的一种hpv分型及整合的高通量测序检测方法，包括如下步骤:1.样本文库制备:(1)、超声片段化：起始量为3ug，加无核酶的水到100ul稀释到30ng/ul。采用scientz08-ⅲ杯式超声波细胞粉碎机进行超声片段化，设定参数为：功率70%，打断3s，停止1s，循环30-60min。(2)、将磁珠提前半个小时拿出以恢复到室温，震荡混匀，取180ul磁珠加入超声打断后的pcr产物中，抽打混匀，温室孵育5min。(3)、将pcr单管放于磁力架上，静置5min，去上清，保持pcr管在磁力架上，加入80%乙醇（新配）200ul漂洗，静置30s后去上清，80%乙醇漂洗两遍，室温干燥10min（还有2min时用小抢头吸残液），直至无乙醇残留。（4）将pcr管从磁力架上取出，加入65ul超纯灭菌水，轻弹混匀，静置5min，短暂离心，将pcr单管置于磁力架上，静置5min后小心吸取上清60ul于新的pcr管中。（5）将40ulendprepmix加入（4）中60ul上清的pcr管中，吹打混匀，短暂离心，进行末端修复，反应条件为30℃、30min。(4)、将pcr管从磁力架中取出，加入20ul超纯灭菌水洗脱，轻弹混匀，静置5min，短暂离心，置于磁力架上，待溶液澄清后吸上清17.5ul于新的pcr管中，再在pcr管中加入12.5ulda-tailing（解冻后放4℃，离心后再使用），吹打混匀。进行末端加尾，反应条件为37℃30min,70℃5min,4℃5min。(5)、在上步加尾后的产物中加入ligationmix和adapter各2.5ul，吹打混匀，短暂离心，放入pcr仪中，反应条件为30℃10min，20℃20min。(6)、在上步连接产物中加入15ul的超纯灭菌水，再加入50ul磁珠，震荡混匀，室温孵育5min，短暂离心，放入磁力架上，静置5min，澄清后去上清，保持pcr管在磁力架上，用200ul的80%乙醇（新配）漂洗两次，室温干燥10min，直至无乙醇残留。(7)、将pcr管从磁力架上取出，加入25ul超纯灭菌水，轻弹混匀，静置5min，短暂离心，将pcr单管置于磁力架上，静置5min后小心吸取上清22ul于新的pcr管中。(8)、将pcrprimermix和amplificationmix2解冻后颠倒混匀，在上步连接产物中加入3ul的primermix和25ul的amplificationmix2，吹打混匀，短暂离心，进行pcr扩增，反应条件为第一步95℃3min第二步98℃20s第三步60℃15s第四步72℃30s第五步72℃5min第六步holdat4℃第二----四步扩增15cycle。(9)、用上步扩增产物进行琼脂糖凝胶电泳，观察条带亮度。(10)、将磁珠提前半个小时拿出以恢复到室温，震荡混匀，取50ul加入pcr扩增产物中，震荡混匀，室温孵育5min，短暂离心，置于磁力悬架上，静置5min，待溶液澄清后去上清，用200ul的80%乙醇（新配）漂洗两次，室温干燥10min，直至无乙醇残留。(11)、将pcr管从磁力架中取出，加入40ul超纯灭菌水洗脱，轻弹混匀，静置5min，短暂离心，置于磁力架上，待溶液澄清后吸上清38ul于新的ep管中，用qubit测溶液浓度并记录，-20℃保存。2.样本富集液相杂交(1)、用于杂交的建库后dna总量为500ng，根据检测的溶液浓度计算杂交时加入ep管的溶液体积。(2)、在上步的ep管中加入2.5ul的cot-1dna，2.5ul的ssdna，2ul的p5（500um），2ul的p7-x（500um），短暂离心。(3)、将上步混合好的溶液在60℃真空干燥30min，直至无液体残留。(4)、在上步干燥后的ep管中加入10ul超纯灭菌水，轻弹混匀，水化10min，短暂离心，取10ul加入pcr八连排中。(5)、将上步pcr八连排放入pcr仪中进行pcr反应，反应条件为第一步95℃5min第二步65℃5min第三步65℃无限在95℃结束后，将1ulrnasin和5ulrnabait混合（每个样品所取量），在65℃2.5min时将rnasin和rnabait混合液以及杂交液放入65℃pcr仪中，在65℃5min结束后先在八连排每个孔中加入6ul的混合液，加入时吹打混匀，2.5min后加入10ul（每个样品所取量）杂交液，加入后吹打混匀，65℃过夜。3.捕获(1)、洗磁珠，取dnabeads磁珠震荡混匀，根据样品数取磁珠（每个样品取30ul磁珠），用200ulbeads洗液洗3次。(2)、悬浮磁珠于bindingbuffer中，每个样品取165ul的bindingbuffer。(3)、吸取bindingbuffer和磁珠的结合液165ul于杂交样品中，震荡45min（每隔5min上下颠倒混匀），45min后置于磁力架上，待澄清后吸弃上清。(4)、加165ul的洗液1于上步磁珠中，洗15min，每5min颠倒混匀一次，15min后置于磁力架上，待澄清后吸弃上清。(5)、将洗液2提前放入70℃pcr仪预热，加165ul的洗液2于上步磁珠中，颠倒混匀，在70℃洗10min，10min后置于磁力架上，待澄清后吸弃上清，洗液2洗3次。(6)、在上步八连排的每个孔中加入22ul的超纯灭菌水，3ul的primer和25ul的amplificationmix2，颠倒混匀，进行pcr扩增，扩增程序为第一步95℃3min第二步98℃20s第三步60℃15s第四步72℃30s第五步72℃5min第六步holdat4℃第二----四步扩增15cycle。(7)、将上步八连排放在磁力架上，待溶液澄清后，吸上清进行琼脂糖凝胶电泳，观察条带亮度。(8)、将磁珠提前半个小时拿出以恢复到室温，震荡混匀，吸取上步八连排中的上清于新的pcr单管中，吸取50ul磁珠加入pcr单管中，震荡混匀，室温孵育5min，短暂离心，置于磁力悬架上，静置5min，待溶液澄清后去上清，用200ul的80%乙醇（新配）漂洗两次，室温干燥10min，直至无乙醇残留。(9)、将pcr管从磁力架中取出，加入40ul超纯灭菌水洗脱，轻弹混匀，静置5min，短暂离心，置于磁力架上，待溶液澄清后吸上清38ul于新的ep管中，用qubit测溶液浓度并记录，-20℃保存。以上得到的dna通过illuminanextseq500测序，得到测序的数据。4.hpv分型及整合的高通量测序检测方法数据分析流程(1)测序仪(illuminanextseq500)获取原始短序列并去除测序数据中的接头;(2)用trimgalore进行qc并去除低质量数据;(3)把短序列用bwa-burrows-wheeleralignment(bwa)软件定位到相应的hpv基因组数据（例如alphapapillomavirus7,humanpapillomavirus18(hpv18),completegenomegi|60975|lcl|hpv18ref.1|hpv18ref）相应的位置上(所用到参数:bwaaln-l-l40-i10-k2-t7-e40-m3-f);(4)使用picardtools及自写脚本，针对不同样本不同类型的hpv基因进行归一化处理，统计相应的覆盖度及及短序列数量，得到分型结果。所述基因整合分析流程包括如下步:(1)测序仪(illuminanextseq500)获取原始短序列;(2)用trimgalore进行qc并去除低质量数据;(3)把短序列用bwa-burrows-wheeleralignment(bwa)软件定位到hpv基因组数据相应的位置上(所用到参数:bwaaln-l-l40-i10-k2-t7-e40-m3-f);(4)用factera检测多个基因的全部或一部分的编码区首尾相连。包括染色体易位、中间缺失或染色体倒置所致的结果，得到基因整合结果。上面所述的实施例仅是对本发明的优选实施方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神前提下，本领域普通工程技术人员对本发明技术方案做出的各种变形和改进，均应落入本发明的权利要求书确定的保护范围内。当前第1页12