CN101921748B - 用于高通量检测人类乳头瘤病毒的dna分子标签 - Google Patents
用于高通量检测人类乳头瘤病毒的dna分子标签 Download PDFInfo
- Publication number
- CN101921748B CN101921748B CN2010102137369A CN201010213736A CN101921748B CN 101921748 B CN101921748 B CN 101921748B CN 2010102137369 A CN2010102137369 A CN 2010102137369A CN 201010213736 A CN201010213736 A CN 201010213736A CN 101921748 B CN101921748 B CN 101921748B
- Authority
- CN
- China
- Prior art keywords
- label
- seq
- sequence
- kinds
- pcr
- 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
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/70—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving virus or bacteriophage
- C12Q1/701—Specific hybridization probes
- C12Q1/708—Specific hybridization probes for papilloma
-
- 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)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Zoology (AREA)
- Engineering & Computer Science (AREA)
- Immunology (AREA)
- Wood Science & Technology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Virology (AREA)
- Biotechnology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Analytical Chemistry (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
本发明涉及人类乳头瘤病毒(Human Papilloma Virus，HPV)的检测方法，特别是基于Solexa测序法的HPV检测方法。
Description
技术领域
本发明涉及人类乳头瘤病毒(Human Papilloma Virus，HPV)的检测方法，特别是基于Solexa测序法的HPV检测方法。
背景技术
宫颈癌是世界上女性癌症的第二号杀手，仅次于乳腺癌。每年全世界约有50万新发病例，25万死亡病例，其中发展中国家占2/3。我国是宫颈癌的高发区，占全球宫颈癌总数的10％。研究表明，人乳头瘤病毒(HPV)与宫颈癌有着密切关系，是重要的致癌因子，同时也是引发宫颈癌的必要条件之一。已经表明，超过100种的HPV能够感染皮肤(皮肤类型)或呼吸道和肛门生殖道的粘膜(粘膜类型)，超过40种的HPV能够感染子宫颈。基于诱导的良性的，恶化前的或恶性的病变，将HPV分别分为低危型(例如HPV6，11，42，43和44)和高危型(例如HPV16，18，31，33和45)。因此，对HPV感染的及早发现和正确分型对宫颈癌防治显得至关重要。
目前HPV的检测方法主要有以下几种：①细胞学检查，其利用宫颈刮片细胞学检查或液基薄层细胞学检查，藉由细胞外观形态的转变来加以诊断。对于HPV感染，镜下可见挖空细胞、角化不良以及湿疣外底层细胞。其局限性在于，HPV感染的诊断的灵敏性和特异性较低。②免疫组化方法，其通过检测HPV的衣壳抗原来进一步明确HPV感染，其所得的阳性反应定位明确，判断可靠。但是，衣壳抗原仅在HPV-DNA复制成熟后才产生，因此，诊断为阴性的受试者并不能被确定为不被HPV感染，该方法灵敏性较低。③实时荧光定量PCR法(FQ-PCR)，其主要是应用荧光检测PCR仪，在PCR反应体系中加入荧光基团，利用荧光信号积累对PCR过程中每一个循环产生的扩增产物进行实时监测，从而对模板的初始浓度进行定量。此方法通量较低。④杂交捕获法(主要为HC-II系统)，其是目前唯一获得美国FDA认证的可用于临床的HPV DNA检测方法，并且其已获得欧洲CE及中国SDA认证。该方法使用专用的标本采集保存器；已获得专利保护的全长8000bp的RNA探针以及已获得专利保护的特异性第一抗体。其原理是将核酸探针杂交到受测检体的HPVDNA上，再通过化学荧光或酶反应借助增幅的信号进行检测。该方法所使用的核酸探针主要分为两种：针对低危险型HPV的核酸探针和针对高危险型HPV的核酸探针。该方法可用于HPV的初步筛检，但是不能确定HPV的特定类型，也不能确定多重感染的情况。
将上述检测方法联合使用，可以提高HPV检测的灵敏性并降低检测的假阴性率。然而，这些方法的组合成本相对较高，一般只适用于经济发达地区的HPV检测和宫颈癌筛查。对于经济较不发达地区，特别是在山区和广大农村地区，上述检测方法的联合使用存在着较大的局限性。因此，开发适合的、低成本的HPV检测方法，是亟待解决的问题。
另一方面，目前已知的HPV检测方法，例如上文描述的那些检测方法，通量都较低。当对大规模的样品进行HPV检测时，应用上述方法是耗时耗力的，并且成本高昂。因此，本领域迫切需要新的高通量的、低成本的HPV检测方法。
发明内容
本发明基于Solexa测序法以及PCR标签(PCR index)，开发了新型的HPV检测方法以及用于此的试剂盒。根据本发明的方法和试剂盒不仅可以实现HPV的高通量、低成本检测，而且可以实现HPV的精确分型。
定义
为了更好地理解本发明，下面提供相关术语的定义和解释。
如本文中所使用的，术语“PCR”是指聚合酶链式反应。
如本文中所使用的，术语“Solexa测序法”是指近几年开发的新一代DNA测序法，也称为第二代测序法。Solexa测序法与传统测序法(例如，Sanger测序法)的不同之处在于，其采用边合成边测序的原理进行DNA序列分析。Solexa测序法具有下述优点：1)成本低，仅为传统测序成本的1％；2)通量高，可以同时对多个样品进行测序，并且进行一次的Solexa测序法可以产生大约500亿(50G)个碱基的数据；3)精确性高(高于98.4％)，有效的解决了多聚重复序列的读取问题。另一方面，高测序通量在进行测序的序列的数目确定的情况下，又反过来提高了序列的测序深度(例如，针对每个序列，可以进行多次测序)，从而确保了测序结果的可靠性。如本文所使用的，术语“测序深度”是指一段DNA序列在测序数据中集中出现的次数。测序深度可以通过将测序量除以基因组长度来计算，例如测序深度为10，表示测了10次的整个基因组。
Solexa测序法的应用十分广泛。其可以用于基因组测序，基因分型，基因多态性研究等等。本发明的方法将Solexa测序法用于检测HPV：通过对待分析的样品进行针对HPV的测序，然后使用本领域已知的比对程序，例如BLAST和SOAP，将所得的测序结果与HPV数据库中的参考序列进行比对，从而实现对样品所感染的HPV的精确分型。本文中使用的HPV数据库包含本领域已知的HPV类型的序列，所述序列可见于例如公共数据库，例如NCBI数据库(http://www.ncbi.nlm.nih.gov/)。
如本文中可互换使用的，术语“PCR标签(PCR index)”、“标签(index)”或“引物标签(primer index)”是指添加在PCR引物5′末端的一小段碱基序列，其通过PCR扩增可以用于标记PCR产物，从而辨别不同模板来源的PCR产物的混合物中各个PCR产物的模板来源。通过在引物的5′末端添加标签，可以对PCR产物进行标记，从而可以将多个不同的PCR产物混合成一个文库，用于进一步的分析和处理。文库中各个不同的PCR产物各自具有独特的标签，从而根据各个PCR产物中独特的标签，可以将各个不同的PCR产物互相区分开来，并将其与PCR模板一一对应。
例如，当需要对多个样品进行测序时，可以在用于各个样品的引物的5′末端添加不同的标签，然后用添加了标签的引物分别对各个样品进行PCR反应，从而对各个样品(即，PCR产物)进行标记。PCR反应后，可以将来自各个样品的带有不同标签的PCR产物混合在一起组成一个文库，然后应用高通量的Solexa测序法同时对文库中的各个PCR产物进行测序。最终，在所得的测序数据中，通过独特的标签，可以将测序结果与各个PCR产物(从而样品模板)一一对应。
可以仅在用于PCR扩增的引物对的一条引物中引入标签，也可以在引物对的两条引物中都引入标签。当在引物对的两条引物中都引入标签时，每个PCR引物对与一对标签组合成一对标签引物，其中正向和反向PCR引物的5′端分别具有正向标签和反向标签，并且正反标签和正反引物序列是对应的，且正向标签和反向标签可以是相同的，或不同的。
设计标签时需要考虑多种因素，包括：1)标签序列中应当避免3个或3个以上的单碱基重复序列；2)所有标签的同一位点中碱基A和碱基C的总含量应在所有碱基含量的30％-70％之间，例如，当设计100条不同的标签序列时，每一条标签序列的第二个碱基(即，所谓的同一位点)中A和C占该100条序列第二个碱基总量的30％-70％；3)标签序列本身的GC含量应在40-60％之间；4)标签之间的序列差异应大于4个碱基；5)标签序列中应避免出现与用于测序的引物相似度高的序列；6)当标签序列添加到PCR扩增引物上后，应避免PCR扩增引物形成发卡结构和二聚体等二级结构。
如本文中所使用的，术语“标签引物(index primer)”是指带有标签的引物，其包含2个部分，标签部分和引物部分，其中标签部分用于在PCR扩增反应中标记PCR产物，而引物部分与模板碱基互补配对，用于扩增模板，并且其中标签部分，任选地通过连接序列，连接至引物部分的5′端。
如本文中所使用的，术语“接头(adapter)”或“文库接头(libraryadapter)”是指经设计的一段碱基序列，其可以连接至文库中的扩增后的PCR产物上，从而文库中的所有扩增后的PCR产物都可以借助该接头进行测序，例如，使用针对该接头设计的测序引物(而无需使用针对各PCR产物设计的特异性测序引物)进行测序。优选地，本发明的接头可以通过“PCR-FREE”方法连接至PCR产物上。
如本文中所使用的，术语“PCR-FREE”是指不进行PCR反应而直接将接头连接至PCR产物上，例如通过使用DNA连接酶将接头连接至PCR产物上。利用PCR-Free方法来构建测序文库是本领域技术人员已知的，参见例如，Nature Methods 6，291-295(2009)。“PCR-FREE”方法由于整个过程无需进行PCR而具有以下有利方面：1)减少了纯化步骤，降低了花费的时间和成本；2)减少了非特异性扩增；3)在含有许多序列高度相似的PCR产物的文库的构建过程中，避免由PCR引入错误，从而提高最后的测序结果的准确性。
如本文中所使用的，本发明的方法和试剂盒可以使用至少1种接头。不同的接头可以共有相同的一段序列(在本文中称为“测序序列”)，并且可进一步包含不同的特征序列，从而，不同的接头可以共用相同的引物(其针对相同的测序序列进行设计)进行测序，并且利用独特的特征序列可以辨别多个文库的混合物中各个PCR产物的文库来源，即，对不同文库来源的PCR产物进行进一步的标记。
将标签与具有不同特征序列的接头组合，可以大大提高标记的效率(参见图1)。例如，使用100种标签可以标记100个样品，而使用100种标签和200种具有不同特征序列的接头，可以标记100*200＝20000个样品。
因此，本发明的一个方面提供了一组标签，其包含至少10种，优选至少20种、至少30种、至少40种、至少50种、至少60种、至少70种、至少80种、至少90种或95种标签，并且所述标签具有选自SEQ ID NO：1-95的序列。在一个优选实施方案中，该组标签至少包括SEQ ID NO：1-10，或SEQ ID NO：11-20，或SEQ ID NO：21-20，或SEQ ID NO：31-40，或SEQ ID NO：41-50，或SEQ ID NO：51-60，或SEQ ID NO：61-70，或SEQ ID NO：71-80，或SEQ ID NO：81-90，或SEQ ID NO：91-95所示的标签，或者其任何两个或者多个的组合，例如SEQ ID NO：1-95所示的标签。
在本发明的一个优选实施方案中，本发明的标签用于标记SEQ IDNO：96-106所示的PCR引物，从而用于进行高通量HPV测序、检测或分型。
本发明的一个方面提供了一种标签引物组，其包含11种标签引物，所述标签引物的序列包含标签序列和PCR引物序列，并且所述标签序列，任选地通过连接序列，连接至所述PCR引物序列的5′端，其中
1)所述标签序列选自SEQ ID NO：1-95，且标签引物组中的11种标签引物的每一种的标签序列都相同，和
2)所述11种标签引物的PCR引物序列分别如SEQ ID NO：96-106所示。
本发明的标签引物组可扩增出至少16种大约170bp的产物，其对应于HPV基因组中最为保守的基因区(L1区)中的一段高度保守的DNA序列。因此，本发明的标签引物组可用于HPV的精确分型。
在一个优选实施方案中，本发明的标签引物组可用于HPV测序、检测或分型，从而其可用于医疗用途，例如诊断HPV的存在和确定HPV的类型等，以及非医疗用途，例如构建HPV数据库，鉴定HPV的新的类型和亚型，研究HPV类型分布的地域性特点，流行病学研究和疫苗研制等。在另一个优选实施方案中，本发明的标签引物组可用于制备试剂盒，所述试剂盒可用于HPV测序、检测或分型。
本发明的另一个方面提供了一种标签引物套组，其包含至少10个，优选至少20个、至少30个、至少40个、至少50个、至少60个、至少70个、至少80个、至少90个或95个上文描述的标签引物组。优选地，标签引物套组中，各个标签引物组所使用的标签序列互不相同。更优选，标签引物套组中所使用的标签序列至少包括SEQ IDNO：1-10，或SEQ ID NO：11-20，或SEQ ID NO：21-20，或SEQ ID NO：31-40，或SEQ ID NO：41-50，或SEQ ID NO：51-60，或SEQ ID NO：61-70，或SEQ ID NO：71-80，或SEQ ID NO：81-90，或SEQ ID NO：91-95所示的标签序列，或者它们任何两个或者多个的组合，例如SEQID NO：1-95所示的标签序列。
在一个优选实施方案中，本发明的标签引物套组可用于高通量HPV测序、检测或分型，从而其可用于医疗用途，例如大规模的HPV相关疾病的诊断和HPV的精确分型(为临床诊断和治疗方案选择提供依据)等，以及非医疗用途，例如构建HPV数据库，鉴定HPV的新的类型和亚型，研究HPV类型分布的地域性特点，流行病学研究和疫苗研制等。在另一个优选实施方案中，本发明的标签引物套组可用于制备试剂盒，所述试剂盒可用于高通量HPV测序、检测或分型。
本发明的另一个方面提供了一种试剂盒，其包含上文描述的标签引物组或标签引物套组。优选，本发明的试剂盒还包含至少1种，优选至少2种、至少10种、至少20种、至少30种、至少40种、至少50种、至少100种或至少200种接头。在一个优选实施方案中，所述接头适用于Solexa测序法，例如其可用于构建测序文库，例如所述接头可具有选自SEQ ID NO：121-132的序列。在一个优选实施方案中，通过PCR-FREE方法，例如DNA连接酶法，使用接头来构建测序文库。
在一个优选实施方案中，本发明的试剂盒可用于高通量HPV测序、检测或分型，并且如上所述，其可用于医疗用途和非医疗用途。
本发明的另一个方面提供了用于对一个或多个样品进行HPV测序、检测或分型的方法。所述方法包括使用上文描述的标签引物组或标签引物套组或试剂盒对各个样品的DNA进行扩增，然后进行测序以获得样品的序列的步骤。
本发明的另一个方面提供了用于对一个或多个样品进行HPV测序、检测或分型的方法，其包括下列步骤：
提供n个样品，n为大于等于1的整数，所述样品优选地来自哺乳动物，更优选是人，且优选是脱落细胞；可选地，将待分析的n个样品分成m个小组，m为整数且n≥m≥1；
1)对于每一个样品，使用一个标签引物组对样品的DNA进行扩增，其中
所述标签引物组包含11种标签引物，所述标签引物的序列包含标签序列和PCR引物序列，并且所述标签序列，任选地通过连接序列，连接至所述PCR引物序列的5′端，其中
i)所述标签序列选自SEQ ID NO：1-95，并且所述11种标签引物的每一种的标签序列都相同，和
ii)所述11种标签引物的PCR引物序列分别如SEQ ID NO：96-106所示，
其中，不同的样品使用的标签引物组可以相同或者不同，并且不同的标签引物组使用不同的标签序列；
2)将步骤1)中的使用不同标签引物组进行扩增所获得的扩增产物混合在一起，得到一个或多个PCR产物文库；
3)通过PCR-FREE方法，例如DNA连接酶法给步骤2)中获得的一个或多个PCR产物文库分别添加接头，从而构建一个或多个测序文库，其中不同的测序文库所使用的接头可以相同或者不同，且不同的接头共有相同的测序序列，但具有不同的特征序列；
4)任选地，将步骤3)中获得的使用不同接头的测序文库混合在一起，得到一个或多个文库混合物；
5)对步骤3)中获得的一个或多个测序文库或步骤4)中获得的一个或多个文库混合物分别进行测序，其中利用二代测序技术，优选Pair-End技术(例如Solexa、Illumina Hiseq 2000)进行测序；
6)根据标签引物组的标签引物序列，或者根据标签引物组的标签引物序列以及接头的特征序列，将测序结果与样品一一对应；
其中所述样品优选是脱落细胞，且优选来源于动物，例如人。
在优选的实施方案中，本发明的方法使用至少10个，优选至少20个、至少30个、至少40个、至少50个、至少60个、至少70个、至少80个、至少90个或95个上文描述的标签引物组。进一步优选地，所使用的标签序列至少包括SEQ ID NO：1-10，或SEQ ID NO：11-20，或SEQ ID NO：21-20，或SEQ ID NO：31-40，或SEQ ID NO：41-50，或SEQ ID NO：51-60，或SEQ ID NO：61-70，或SEQ ID NO：71-80，或SEQ ID NO：81-90，或SEQ ID NO：91-95所示的标签序列，或者它们任何两个或者多个的组合，例如SEQ ID NO：1-95所示的标签序列。
在本发明的方法的一个优选实施方案中，本发明的方法使用至少1种，优选至少2种、至少10种、至少20种、至少30种、至少40种、至少50种、至少100种或至少200种接头，例如所述接头可具有选自SEQ ID NO：121-132的序列。
在本发明的方法的一个优选实施方案中，在测序后，将获得的样品序列与上文描述的HPV数据库中的序列进行对比，从而对样品进行HPV精确分型。
在本发明的另一个方面，本发明基于Solexa测序法提供了用于对多个样品进行高通量HPV测序、检测或分型的方法，所述样品优选是脱落细胞，且优选来源于动物，例如人，所述方法包括下列步骤：
1)将待分析的样品分成m个小组，m为≥1的整数；
2)对每个样品小组进行如下步骤：
2a)从待分析的样品中提取DNA；
2b)根据用于扩增HPV DNA的引物组中的所有引物的序列，设计一套标签，标签的个数n等于该样品小组的样品数目；
2c)将2b)中设计的每一个标签分别添加至引物组的所有正向引物或所有反向引物或所有引物的序列的5′末端，从而提供n个标签引物组；
2d)使用2c)中提供的标签引物组对2a)中获得的样品DNA进行PCR扩增，从而提供PCR产物，其中针对每个样品DNA，使用不同的标签引物组；和
2e)将2d)中的所有PCR产物混合在一起，获得PCR产物文库；
3)给2)中得到的PCR产物文库加上接头，其中各个PCR产物文库使用不同的接头，从而构建m个测序文库，其中各个接头共有相同的测序序列，而具有互不相同的特征序列；
4)将m个测序文库混合在一起，利用二代测序技术，优选Pair-End技术(例如Solexa、Illumina Hiseq 2000)进行测序；得到所有样品的测序结果；
5)根据测序结果中的接头的特征序列、标签的序列以及引物的序列，将获得的测序结果与样品一一对应；和任选地，将每个样品的测序结果与HPV数据库进行比对，从而实现HPV测序、检测或分型。
在一个优选实施方案中，使用本领域技术人员熟知的方法进行DNA提取。例如，可以使用自动化的DNA提取仪和DNA提取试剂盒进行DNA提取，例如可商购获得的KingFisher自动提取仪，例如美国Thermo Scientific Kingfisher Flex全自动磁珠提取纯化系统。
在一个优选的实施方案中，2b)中的引物组包含11条引物，其序列分别如SEQ ID NO：96-106所示。这11条引物组成的引物组可扩增出至少16种大约170bp的产物，其对应于HPV基因组中最为保守的基因区(L1区)中的一段高度保守的DNA序列。因此，通过对该扩增产物进行精确测序，可以实现HPV的精确分型。
在又一个优选的实施方案中，2b)中设计的标签的个数为至少10个，优选至少20个，至少30个，至少40个，至少50个，至少60个，至少70个，至少80个，至少90个，或至少100个。优选，标签可具有选自SEQ ID NO：1-95的序列。在优选的实施方案中，不同样品小组之间的所使用的标签可以相同或者不同。在优选的实施方案中，引入正向引物的标签与引入反向引物的标签可以相同或者不同。在特别优选的实施方案中，2b)中设计的标签至少包括SEQ ID NO：1-10，或SEQ ID NO：11-20，或SEQ ID NO：21-20，或SEQ ID NO：31-40，或SEQ ID NO：41-50，或SEQ ID NO：51-60，或SEQ ID NO：61-70，或SEQ ID NO：71-80，或SEQ ID NO：81-90，或SEQ ID NO：91-95所示的标签，或者它们任何两个或者多个的组合。在特别优选的实施方案中，2b)中设计的标签如SEQ ID NO：1-95所示。
在一个优选的实施方案中，使用PCR-FREE方法，给PCR产物文库加上接头，例如使用DNA连接酶。特别地，在本发明的方法中，由于不同HPV类型之间的DNA序列高度相似，因此，本发明的测序文库的构建必须通过PCR-FREE方法来完成。相反，如果通过采用常规的pooling PCR将接头连接至PCR产物上来构建测序文库，那么所得到的文库中将含有大量的与原模板不一致的产物，导致不能对原模板进行准确测序。在一个优选的实施方案中，使用的接头的数目为至少1种，优选至少2种、至少10种、至少20种、至少30种、至少40种、至少50种、至少100种或至少200种，例如所述接头可具有选自SEQID NO：121-132的序列。
在本发明的方法的优选实施方案中，接头是可商购获得的接头，例如可购自Illumina公司的PCR-free Index Adapter Oligo Mix。在另外的实施方案中，本发明也可使用下列PCR-free接头(下划线部分为接头的特征序列)。
PCR-free接头1(SEQ ID NO：121)：ATCACG
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头2(SEQ ID NO：122)：CGATGT
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACCGATGTATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头3(SEQ ID NO：123)：TTAGGC
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACTTAGGCATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头4(SEQ ID NO：124)：TGACCA
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACTGACCAATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头5(SEQ ID NO：125)：ACAGTG
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACACAGTGATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头6(SEQ ID NO：126)：GCCAAT
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头7(SEQ ID NO：127)：CAGATC
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACCAGATCATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头8(SEQ ID NO：128)：ACTTGA
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACACTTGAATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头9(SEQ ID NO：129)：GATCAG
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头10(SEQ ID NO：130)：TAGCTT
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACTAGCTTATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头11(SEQ ID NO：131)：GGCTAC
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCTTCTGCTTG
PCR-free接头12(SEQ ID NO：132)：CTTGTA
5-Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCTTCTGCTTG
在一个优选的实施方案中，本发明的方法使用Solexa测序仪(例如，Illumina Genome Analyzer IIx测序仪)进行Solexa测序法。在另外的优选的实施方案中，HPV数据库包含本领域已知的HPV类型的序列，所述序列例如可见于例如公共数据库，例如NCBI数据库(http://www.ncbi.nlm.nih.gov/)。
在本发明的方法的一个优选实施方案中，样品可以是脱落细胞。在另一个优选实施方案中，样品可以来源于动物，优选哺乳动物，更优选人。
发明的有益效果
本发明的新型的HPV检测方法和用于此的试剂盒相对于现有技术具有以下有利方面。
1)高通量。通过使用标签与具有不同特征序列的接头，进行一次本发明的方法，可检测甚至10000份样品，从而本发明的方法可以广泛用于疾病普查工作，成为早期诊断疾病的有效手段。
2)低成本。本发明采用Solexa测序法进行测序，测序成本大大降低(仅为常规测序法的1％)，从而HPV检测成本大大降低。
3)可对HPV进行精确分型。通过使用多条引物(例如本发明的6条正向引物和5条反向引物)进行扩增，并将扩增产物的序列信息与HPV数据库进行比对，可以精确地确定HPV的类型，从而为临床诊断和治疗方案选择提供依据。
另外，使用本发明的方法还有利于发现新的HPV类型，包括已有的类型中的新的亚型和变异体，为科学研究提供更有效和方便的工具。
下面将结合附图和实施例对本发明的实施方案进行详细描述，但是本领域技术人员将理解，下列附图和实施例仅用于说明本发明，而不是对本发明的范围的限定。根据附图和优选实施方案的下列详细描述，本发明的各种目的和有利方面对于本领域技术人员来说将变得显然。
附图说明
图1为经标签和具有独特的特征序列的接头标记后的PCR产物的示意图。在本发明的示例性方法中，通过PCR在每个样品的PCR产物两端同时引入标签；把多个带有不同标签的PCR产物混合在一起，用于构建测序文库。在测序文库的构建过程中，当需要时，可以构建多个测序文库时，其中通过使用具有不同特征序列的接头，来标记各个测序文库。文库构建完毕后，经具有不同特征序列的接头标记的多个测序文库可以混合在一起，并用Solexa测序法进行同时测序(不同的测序文库之间所使用的标签可以相同或不同)。最后，根据测序结果中的接头的特征序列和标签的序列信息，可将测序结果与每个样品一一对应。
图2为部分PCR产物的琼脂糖凝胶电泳图。从电泳图上可观察到，PCR产物条带的大小为约170bp。其中，泳道M是50bp DNA分子梯，泳道1-14为随机挑选的HPV阳性样本的PCR产物。
具体实施方式
采用本发明的方法，对190份已知HC-II结果的样品进行HPV基因分型，结果发现：采用本发明的方法所得的结果不仅与已知的HC-II结果一致，而且实现了对HPV的精确分型。
实施例1：样品提取
根据厂商的说明书，使用KingFisher自动提取仪(美国ThermoScientific Kingfisher Flex全自动磁珠提取纯化系统)从190份已知HC-II结果的脱落细胞中提取DNA。使用程序“Bioeasy-200μl BloodDNA-KF.msz”，进行核酸提取。程序结束后，获得100μl左右的洗脱产物(提取的DNA)，用作下一步PCR扩增中的模板。
实施例2：PCR扩增
把实施例1中获得的190份DNA依次编号为1-190，并平均分为2组(HPV-1组：编号1-95；HPV-2组：编号96-190)。根据用于扩增HPV DNA的引物组(包括6条正向引物和5条反向引物)的各条引物的序列(表2，SEQ ID NO：96-107)，设计一套标签，共95个(表1，SEQ ID NO：1-95)。将设计的每一个标签分别添加至引物组的各条引物的序列的5′末端，从而获得95个标签引物组，其中每个标签引物组包括相应的6条正向标签引物和5条反向标签引物，并且不同的标签引物组使用不同的标签(即，95个标签引物组与95个标签一一对应)。
在96孔板中对所有样品进行PCR反应，共使用2块板(HPV-1组和HPV-2组各1块板)。使用实施例1中获得的DNA作为模板，并且在HPV-1组和HPV-2组(各95个样品)中，针对每个样品，使用不同的标签引物组进行PCR扩增(即，95个样品与95个标签引物组一一对应)。记录下每一个标签引物组(每一个标签)对应的样品的编号信息。每块板中还设置一个不添加模板的阴性对照。2块板中的阴性对照所用的引物分别与样品1和样品96所用的引物相同。
表1：标签与样品的相关信息。
| 标签编号 | 标签序列 | 对应96孔板位置 | 对应样品(组1) | 对应样品(组2) | SEQ ID NO： |
| PI-1 | GCTGCGACTC | A1 | 1 | 96 | 1 |
| PI-2 | GTGTAGATAC | A2 | 2 | 97 | 2 |
| PI-3 | CTGATATCTA | A3 | 3 | 98 | 3 |
| PI-4 | ACGATGCTAT | A4 | 4 | 99 | 4 |
| PI-5 | TAGACTAGAC | A5 | 5 | 100 | 5 |
| PI-6 | CTGTCTGTGT | A6 | 6 | 101 | 6 |
| PI-7 | GCATACTGAC | A7 | 7 | 102 | 7 |
| PI-8 | CTGCTCGCAT | A8 | 8 | 103 | 8 |
| PI-9 | CATGAGTAGA | A9 | 9 | 104 | 9 |
| PI-10 | TCTCACTATG | A10 | 10 | 105 | 10 |
| PI-11 | TGTACTACTA | A11 | 11 | 106 | 11 |
| PI-12 | GTAGACTAGT | A12 | 12 | 107 | 12 |
| PI-13 | ATATGCTACT | B1 | 13 | 108 | 13 |
| PI-14 | CACTCGCTGT | B2 | 14 | 109 | 14 |
| PI-15 | CATCACGCAC | B3 | 15 | 110 | 15 |
| PI-16 | AGCATGTGAT | B4 | 16 | 111 | 16 |
| PI-17 | AGCTAGTAGA | B5 | 17 | 112 | 17 |
| PI-18 | GCTATGTAGT | B6 | 18 | 113 | 18 |
| PI-19 | TACGATGATG | B7 | 19 | 114 | 19 |
| PI-20 | TACGCTGTAC | B8 | 20 | 115 | 20 |
| PI-21 | TATGTGTACT | B9 | 21 | 116 | 21 |
| PI-22 | TGACTCAGAC | B10 | 22 | 117 | 22 |
| PI-23 | TCGTAGCTCA | B11 | 23 | 118 | 23 |
| PI-24 | GAGACTCGTA | B12 | 24 | 119 | 24 |
| PI-25 | CTAGATGTCA | C1 | 25 | 120 | 25 |
| PI-26 | GATGACTCTC | C2 | 26 | 121 | 26 |
| PI-27 | TCAGTCGCAC | C3 | 27 | 122 | 27 |
| PI-28 | TGTAGTGAGT | C4 | 28 | 123 | 28 |
| PI-29 | TCATCGTAGA | C5 | 29 | 124 | 29 |
| PI-30 | TAGCATCTGT | C6 | 30 | 125 | 30 |
| PI-31 | TAGTAGTCGT | C7 | 31 | 126 | 31 |
| PI-32 | CTATACGTGC | C8 | 32 | 127 | 32 |
| PI-33 | CGACTGTAGA | C9 | 33 | 128 | 33 |
| PI-34 | GATGTCATGT | C10 | 34 | 129 | 34 |
| PI-35 | GTCTCGACTG | C11 | 35 | 130 | 35 |
| PI-36 | AGCTGACGAT | C12 | 36 | 131 | 36 |
| PI-37 | ATGATATAGT | D1 | 37 | 132 | 37 |
| PI-38 | ATGTGCTCTA | D2 | 38 | 133 | 38 |
| PI-39 | CTCACTCGAT | D3 | 39 | 134 | 39 |
| PI-40 | GCTGCGACTC | D4 | 40 | 135 | 40 |
| PI-41 | GAGTCATGTC | D5 | 41 | 136 | 41 |
| PI-42 | CATACGCTCA | D6 | 42 | 137 | 42 |
| PI-43 | CACTCTCGTC | D7 | 43 | 138 | 43 |
| PI-44 | GCACTAGATG | D8 | 44 | 139 | 44 |
| PI-45 | AGTACGCATG | D9 | 45 | 140 | 45 |
| PI-46 | TCTGTGACGT | D10 | 46 | 141 | 46 |
| PI-47 | TAGCTCATCT | D11 | 47 | 142 | 47 |
| PI-48 | AGCATACACT | D12 | 48 | 143 | 48 |
| PI-49 | GCTATAGTCA | E1 | 49 | 144 | 49 |
| PI-50 | CGTCTCATGC | E2 | 50 | 145 | 50 |
| PI-51 | GCTACTACGT | E3 | 51 | 146 | 51 |
| PI-52 | GAGTGTACTA | E4 | 52 | 147 | 52 |
| PI-53 | GTCATACGTG | E5 | 53 | 148 | 53 |
| PI-54 | TATGAGAGAT | E6 | 54 | 149 | 54 |
| PI-55 | ATCTGAGTAC | E7 | 55 | 150 | 55 |
| PI-56 | CGATAGCATC | E8 | 56 | 151 | 56 |
| PI-57 | ACTGATCTCA | E9 | 57 | 152 | 57 |
| PI-58 | CTCGATACTA | E10 | 58 | 153 | 58 |
| PI-59 | CATGTGACTG | E11 | 59 | 154 | 59 |
| PI-60 | CGCATCACTA | E12 | 60 | 155 | 60 |
| PI-61 | GCATATATCT | F1 | 61 | 156 | 61 |
| PI-62 | CTGATGCGAC | F2 | 62 | 157 | 62 |
| PI-63 | TCTCAGAGTC | F3 | 63 | 158 | 63 |
| PI-64 | CAGTGCGAGT | F4 | 64 | 159 | 64 |
| PI-65 | ATCTCTGATG | F5 | 65 | 160 | 65 |
| PI-66 | GCTAGTAGTC | F6 | 66 | 161 | 66 |
| PI-67 | ATGAGTCGTC | F7 | 67 | 162 | 67 |
| PI-68 | ATCACTCAGA | F8 | 68 | 163 | 68 |
| PI-69 | TCTCTCTGAT | F9 | 69 | 164 | 69 |
| PI-70 | CTCTAGTGCT | F10 | 70 | 165 | 70 |
| PI-71 | CGTCGTGCTA | F11 | 71 | 166 | 71 |
| PI-72 | CGACTACTAT | F12 | 72 | 167 | 72 |
| PI-73 | GCACGTCGAT | G1 | 73 | 168 | 73 |
| PI-74 | GTAGTGCTCT | G2 | 74 | 169 | 74 |
| PI-75 | CTGACGAGCT | G3 | 75 | 170 | 75 |
| PI-76 | CTATAGTCTA | G4 | 76 | 171 | 76 |
| PI-77 | ACACGCACTA | G5 | 77 | 172 | 77 |
| PI-78 | CTCGCACTAC | G6 | 78 | 173 | 78 |
| PI-79 | AGATCTCACT | G7 | 79 | 174 | 79 |
| PI-80 | ATACTAGTGT | G8 | 80 | 175 | 80 |
| PI-81 | ATATCTCGTA | G9 | 81 | 176 | 81 |
| PI-82 | TGACTGCGTA | G10 | 82 | 177 | 82 |
| PI-83 | TGTAGACGTA | G11 | 83 | 178 | 83 |
| PI-84 | AGAGACTATG | G12 | 84 | 179 | 84 |
| PI-85 | GTCGAGTCAC | H1 | 85 | 180 | 85 |
| PI-86 | TGACAGCTAC | H2 | 86 | 181 | 86 |
| PI-87 | CGCTAGACAT | H3 | 87 | 182 | 87 |
| PI-88 | CGTAGATATG | H4 | 88 | 183 | 88 |
| PI-89 | TGAGTCTGCT | H5 | 89 | 184 | 89 |
| PI-90 | TAGTCGTATG | H6 | 90 | 185 | 90 |
| PI-91 | CATACACGAC | H7 | 91 | 186 | 91 |
| PI-92 | CGCTCAGAGA | H8 | 92 | 187 | 92 |
| PI-93 | GTGAGTCTCA | H9 | 93 | 188 | 93 |
| PI-94 | GACAGATGAT | H10 | 94 | 189 | 94 |
| PI-95 | GCTGTGCGAC | H11 | 95 | 190 | 95 |
表2：未添加标签的用于扩增HPV DNA的引物组的各条引物的序列信息。
| 引物编号 | 引物序列 | SEQ ID NO： |
| F1 | TTTGTTACTGTGGTAGATACTAC | 96 |
| F2 | TTTGTTACTGTGGTGGATACTAC | 97 |
| F3 | TTTGTTACCGTTGTTGATACTAC | 98 |
| F4 | TTTGTTACTAAGGTAGATACCACTC | 99 |
| F5 | TTTGTTACTGTTGTGGATACAAC | 100 |
| F6 | TTTGTTACTATGGTAGATACCACAC | 101 |
| R1 | GAAAAATAAACTGTAAATCATATTCCT | 102 |
| R2 | GAAAAATAAATTGTAAATCATACTC | 103 |
| R3 | GAAATATAAATTGTAAATCAAATTC | 104 |
| R4 | GAAAAATAAACTGTAAATCATATTC | 105 |
| R5 | GAAAAATAAACTGCAAATCATATTC | 106 |
注：F表示正向引物，R表示反向引物。
使用下列PCR参数进行扩增：
95℃30秒→48℃30秒→72℃30秒(40个循环)
72℃10分钟→12℃∞
PCR反应体系为25μl，其组成是(所有试剂均购自Enzymatics公司)：
| 试剂 | 体积/反应 |
| 水(HPLC级) | 14.375μl |
| 10x Ex Taq缓冲液(Mg2+ plus) | 2.5μl |
| dNTP混合物(各2.5mM) | 2μl |
| 带有标签的F1/F2/F3/F4/F5/F6的混合物(各7.5pmol) | 0.5μl |
| 带有标签的R1/R2/R3/R4/R5的混合物(各7.5pmol) | 0.5μl |
| Ex Taq HS(5U/μl) | 0.125μl |
| 模板DNA | 5μl |
| 总体积 | 25μl |
PCR反应在Bio-Rad公司的PTC-200PCR仪上运行。PCR完成后，取3μl PCR产物在2.5％的琼脂糖凝胶上进行电泳检测(图2)。
实施例3：PCR产物的混合和纯化
把HPV-1组与HPV-2组中剩余的PCR产物各混合在一个3ml的EP管中(同样标记为HPV-1组和HPV-2组)，并震荡混匀。从2管混合物中各取出500μl DNA，并根据厂商的说明书，使用Qiagen DNAPurification试剂盒进行过柱纯化，得到200μl DNA。使用Nanodrop8000(Thermo Fisher Scientific公司)，测定纯化后的混合物的DNA浓度分别为98ng/μl(HPV-1组)和102ng/μl(HPV-2组)。
实施例4：Solexa测序文库的构建
4.1：末端修复反应
使用Thermomixer(Eppendorf公司)，对实施例3中获得的经纯化的两管DNA混合物分别进行DNA末端修复反应。修复反应的反应体系为100μl，其组成是(所有试剂均购自Enzymatics公司)：
| 试剂 | 体积/反应 |
| 上一步所得的DNA | 75μl |
| 20x多核苷酸激酶缓冲液(B904) | 10μL |
| dNTP混合物(各20mM) | 4μL |
| T4 DNA聚合酶 | 5μL |
| Klenow片段 | 1μL |
| T4多核苷酸激酶 | 5μL |
| 总体积 | 100μL |
反应条件为：20℃，30分钟。
根据厂商的说明书，使用QIAquick PCR Purification试剂盒纯化并回收DNA末端修复反应的产物。回收的产物溶于34μl EB(QIAGENElution Buffer)中。
4.2：3′末端加A反应
使用Thermomixer(Eppendorf公司)，对回收的DNA进行3′末端加A反应。反应体系为50ul，其组成是(所有试剂均购自Enzymatics公司)：
| 试剂 | 体积/反应 |
| 上一步所得的DNA | 32μl |
| dATP(1mM，GE公司) | 10μl |
| 10x Blue缓冲液 | 5μl |
| Klenow (3′-5′exo-) | 3μl |
| 总体积 | 50μl |
反应条件为：37℃，30分钟。
根据厂商的说明书，使用MiniElute PCR Purification Kit(QIAGEN公司)纯化并回收3′末端加A反应的产物。回收的产物溶于20μl的EB中。
4.3：添加Solexa接头
使用Thermomixer(Eppendorf公司)，对上一步获得的2种产物分别添加不同的接头以构建2个测序文库。记录下接头和文库的对应关系。
添加Solexa接头的反应体系为50ul，其组成是(所有试剂均购自Illumina公司)：
| 试剂 | 体积/反应 |
| 上一步所得的DNA | 11μL |
| 2x快速连接缓冲液 | 15μL |
| PCR-free Index Adapter Oligo Mix(25mM) | 1μL |
| T4DNA连接酶(快速，L603-HC-L) | 3μL |
| 总体积 | 30μL |
反应条件为：20℃，15分钟。
根据厂商的说明书，使用Ampure Beads(Beckman CoulterGenomics)纯化反应产物并将产物溶于17μl去离子水。使用AgilentBioanalyzer 2100(Agilent公司)和荧光定量PCR(QPCR)测定产物的DNA浓度，结果如下：
| 2100(nM) | qPCR(nM) | |
| HPV-1组 | 20.4 | 24.2 |
| HPV-2组 | 21.6 | 25.8 |
实施例5：Solexa测序
以Agilent Bioanalyzer 2100所测浓度为准，将上一步所得的2种产物等摩尔混合在一起(各取10pmol DNA)。根据厂商的说明书，使用Solexa测序仪(Illumina Genome Analyzer IIx测序仪)，用Solexa PE-75程序进行测序。
实施例6：结果分析
根据测序结果中的接头的特征序列以及标签引物(标签部分和引物部分)的序列的信息，将测序结果与每个样品一一对应。然后，使用本领域已知的比对程序，例如BLAST和SOAP，将每个样品的测序结果与HPV数据库进行比对，从而实现HPV检测，并对HPV进行精确分型。
所得到的检测结果与已知的结果完全一致(参见表3)，表明本发明方法可以用于精确检测样品中的HPV。
表3：190个样品的检测结果。
| 样品编号 | 原HC-II结果(RLU/CO值) | 本次检测的HPV结果 | 样品编号 | 原HC-II结果(RLU/CO值) | 本次检测的HPV结果 |
| 1 | 14.2 | HPV56 | 96 | 0.28 | 阴性 |
| 2 | 0.31 | 阴性 | 97 | 0.33 | 阴性 |
| 3 | 196.41 | HPV16 | 98 | 181.29 | HPV35，HPV6 |
| 4 | 5.76 | HPV18 | 99 | 77.32 | HPV16 |
| 5 | 0.35 | 阴性 | 100 | 91.22 | HPV39 |
| 6 | 99.86 | HPV18，HPV11，HPV16 | 1011 | 188.92 | HPV52 |
| 7 | 128.86 | HPV39 | 102 | 1352.83 | HPV35，HPV11，HPV39 |
| 8 | 35.12 | HPV18，HPV6 | 103 | 1.39 | HPV43 |
| 9 | 498.69 | HPV16，HPV56 | 104 | 119.5 | HPV45，HPV11 |
| 10 | 603.57 | HPV18，HPV31，HPV39 | 105 | 292.43 | HPV56，HPV31 |
| 11 | 0.27 | 阴性 | 106 | 2.91 | HPV68 |
| 12 | 3420.57 | HPV18 | 107 | 193.13 | HPV45 |
| 13 | 0.38 | 阴性 | 108 | 2.62 | HPV6 |
| 14 | 0.41 | 阴性 | 109 | 94.12 | HPV16 |
| 15 | 455.06 | HPV16 | 110 | 792.72 | HPV18，HPV31 |
| 16 | 8.93 | HPV18 | 111 | 31.76 | HPV11 |
| 17 | 0.6 | 阴性 | 112 | 0.25 | 阴性 |
| 18 | 0.41 | 阴性 | 113 | 0.23 | 阴性 |
| 19 | 0.29 | 阴性 | 114 | 750.82 | HPV56，HPV16 |
| 20 | 27.64 | HPV31 | 115 | 0.4 | 阴性 |
| 21 | 1985.41 | HPV56，HPV68 | 116 | 2.75 | HPV31 |
| 22 | 20.71 | HPV42 | 117 | 396.04 | HPV45 |
| 23 | 1795.83 | HPV11，HPV16，HPV52 | 118 | 354.76 | HPV18，HPV16 |
| 24 | 9.55 | HPV43 | 119 | 6.26 | HPV11 |
| 25 | 237.62 | HPV39 | 120 | 1719.67 | HPV16，HPV45 |
| 26 | 1.5 | HPV6 | 121 | 76.92 | HPV51 |
| 27 | 1478.98 | HPV68，HPV16 | 122 | 1318.02 | HPV56，HPV16，HPV42 |
| 28 | 115.31 | HPV44 | 123 | 0.28 | 阴性 |
| 29 | 419.31 | HPV16 | 124 | 0.33 | 阴性 |
| 30 | 1.81 | candHPV89 | 125 | 181.29 | HPV59 |
| 31 | 2013.61 | HPV52，HPV39 | 126 | 77.32 | HPV68 |
| 32 | 1379.09 | HPV54，HPV33 | 127 | 110.8 | HPV52 |
| 33 | 12.74 | HPV42 | 128 | 147.25 | HPV16 |
| 34 | 1695.31 | HPV16，candHPV89 | 129 | 0.24 | HPV26 |
| 35 | 1410.85 | HPV35 | 130 | 1.55 | HPV11，HPV53 |
| 36 | 1149.25 | HPV18 | 131 | 2.03 | HPV6，HPV66 |
| 37 | 0.24 | 阴性 | 132 | 8.45 | HPV43 |
| 38 | 1.55 | HPV11 | 133 | 0.2 | 阴性 |
| 39 | 2.03 | HPV11，HPV6 | 134 | 0.24 | 阴性 |
| 40 | 8.45 | HPV42 | 135 | 10.53 | HPV11 |
| 41 | 0.2 | 阴性 | 136 | 1410.85 | HPV16，HPV53，HPV70 |
| 42 | 0.22 | 阴性 | 137 | 1149.25 | HPV56，HPV81，HPV73 |
| 43 | 0.53 | 阴性 | 138 | 0.24 | 阴性 |
| 44 | 10.38 | HPV6 | 139 | 413.9 | HPV45 |
| 45 | 78.21 | HPV16 | 140 | 17.05 | HPV11 |
| 46 | 0.23 | 阴性 | 141 | 23.6 | HPV52 |
| 47 | 45.42 | HPV16，HPV18 | 142 | 3379.09 | HPV16，HPV35，HPV56 |
| 48 | 0.35 | 阴性 | 143 | 0.18 | 阴性 |
| 49 | 148.66 | HPV18，candHPV89 | 144 | 1.46 | HPV18 |
| 50 | 60.27 | HPV56 | 145 | 1.25 | HPV11，HPV26 |
| 51 | 0.28 | 阴性 | 146 | 2.13 | HPV6，HPV81 |
| 52 | 360.26 | HPV56，HPV68 | 147 | 872.52 | HPV16，HPV45，HPV52 |
| 53 | 50.31 | HPV18 | 148 | 1.5 | HPV18 |
| 54 | 0.18 | 阴性 | 149 | 4.33 | HPV16 |
| 55 | 0.31 | 阴性 | 150 | 0.82 | 阴性 |
| 56 | 196.41 | HPV16 | 151 | 60.35 | HPV59 |
| 57 | 5.76 | HPV51 | 152 | 0.24 | 阴性 |
| 58 | 0.23 | 阴性 | 153 | 0.23 | 阴性 |
| 59 | 0.88 | 阴性 | 154 | 0.18 | 阴性 |
| 60 | 0.16 | 阴性 | 155 | 1.46 | HPV51 |
| 61 | 870.63 | HPV52，HPV16 | 156 | 11.25 | HPV16 |
| 62 | 10.18 | HPV42 | 157 | 2.13 | HPV11 |
| 63 | 0.15 | 阴性 | 158 | 0.13 | 阴性 |
| 64 | 1.36 | HPV11 | 159 | 90.18 | HPV58 |
| 65 | 68.2 | HPV59 | 160 | 0.15 | 阴性 |
| 66 | 0.68 | 阴性 | 161 | 602.79 | HPV68，HPV16 |
| 67 | 130.41 | HPV45 | 162 | 132.68 | HPV56，HPV11 |
| 68 | 0.26 | 阴性 | 163 | 127.08 | HPV39，HPV54 |
| 69 | 5.25 | HPV6 | 164 | 602.79 | HPV33 |
| 70 | 0.46 | 阴性 | 165 | 276 | HPV18 |
| 71 | 8.23 | HPV40 | 166 | 243.6 | HPV45 |
| 72 | 0.28 | 阴性 | 167 | 229.44 | HPV51 |
| 73 | 100.16 | HPV43，HPV44 | 168 | 1384.92 | HPV16，HPV58，HPV72 |
| 74 | 450.13 | HPV41 | 169 | 172.64 | HPV58 |
| 75 | 127.08 | HPV39，HPV6 | 170 | 855.24 | HPV16，candHPV89 |
| 76 | 602.79 | HPV45 | 171 | 126.47 | HPV51 |
| 77 | 276 | HPV16 | 172 | 86.62 | HPV44，HPV11 |
| 78 | 243.6 | HPV6，HPV70，HPV39 | 173 | 879.37 | HPV18，HPV58 |
| 79 | 229.44 | HPV35 | 174 | 119.39 | HPV56 |
| 80 | 1384.92 | HPV52，HPV56，HPV11 | 175 | 0.61 | 阴性 |
| 81 | 172.64 | HPV26，HPV42 | 176 | 18.02 | HPV16 |
| 82 | 855.24 | HPV35，HPV6 | 177 | 16.06 | HPV18 |
| 83 | 620.69 | HPV52 | 178 | 60.69 | HPV56，HPV11 |
| 84 | 128.02 | HPV11 | 179 | 2.45 | HPV11 |
| 85 | 514.84 | HPV33 | 180 | 94.93 | HPV39 |
| 86 | 68.3 | HPV58 | 181 | 1635.3 | HPV16，HPV35，HPV51 |
| 87 | 402.15 | HPV59，HPV16 | 182 | 754.64 | HPV33，candHPV89 |
| 88 | 51.72 | HPV33 | 183 | 0.23 | HPV11 |
| 89 | 1.78 | HPV6 | 184 | 20.28 | HPV18 |
| 90 | 56.7 | HPV11，HPV31 | 185 | 0.16 | 阴性 |
| 91 | 186.06 | HPV16 | 186 | 0.13 | 阴性 |
| 92 | 0.02 | 阴性 | 187 | 60.18 | HPV59 |
| 93 | 386.06 | HPV18，HPV16 | 188 | 0.15 | 阴性 |
| 94 | 28.09 | HPV6，HPV44 | 189 | 1.36 | HPV43 |
| 95 | 186.06 | HPV68， | 190 | 0.28 | 阴性 |
此外，本发明的方法还可以对样品中的HPV进行精确分型。表4中提供了图2所示的泳道1-14所对应的样品的测序序列和分型结果。
参考文献
本文中用于举例说明本发明或提供关于本发明的实施的另外的详细内容的专利、出版物和其他材料通过引用合并入本文，并且为方便起见按下列文献目录提供。
[1].Pectasides D，Kanposioras K，Papaxoinis G et al.Chemotherapy for recurrent cervical cancer.Cancer TreatmentReviews，2008，34(7)：603-613.
[2].Brink，A.A.，P.J.Snijders，and C.J.Meijer.HPVdetection methods.Dis.Markers 2007，23：273-281.
[3].IARC.Handbooks of cancer prevention.Cervix cancerscreening[R].Lyon：IARC Press，2005.
[4].Doorbar，J.Molecular biology of human papillomavirusinfection and cervical cancer.Clin.Sci.2006，110：525-541.
[5].Cox T，Cuzick J.HPV DNA testing in cervical cancerscreening：From evidence to policies.Gyneeol Oncol，2006，103：8-11.
[6].Kulmala S，syIjhen.Human papillomavirus testing withthe hybrid capture assay and PCR as screening tools.ClinMicrobiol，2004，42(6)：2470-2475.
[7].Quail，M.et al.，A large genome center’s improvementsto the Illumina sequencing system.Nat.Methods，2008，5，1005-1010.
[8].Brown，C.G.et al.，Solexa/Illumina GAPipeline productand product documentation，Illumina Inc，2006.
[9].Lozano，R.Successfully integrating humanpapillomavirus testing into your practice.Arch.Pathol.LabMed，2003，127：991-994.
Claims (40)
1.一组标签，其包含至少10种标签，所述标签的序列选自SEQ IDNO：1-95。
2.权利要求1的一组标签，其包含至少20种标签。
3.权利要求1的一组标签，其包含至少30种标签。
4.权利要求1的一组标签，其包含至少40种标签。
5.权利要求1的一组标签，其包含至少50种标签。
6.权利要求1的一组标签，其包含至少60种标签。
7.权利要求1的一组标签，其包含至少70种标签。
8.权利要求1的一组标签，其包含至少80种标签。
9.权利要求1的一组标签，其包含至少90种标签。
10.权利要求1的一组标签，其包含95种标签。
11.权利要求1的一组标签，其至少包括SEQ ID NO：1-10，或SEQID NO：11-20，或SEQ ID NO：21-20，或SEQ ID NO：31-40，或SEQID NO：41-50，或SEQ ID NO：51-60，或SEQ ID NO：61-70，或SEQID NO：71-80，或SEQ ID NO：81-90，或SEQ ID NO：91-95所示的标签，或者其任何两个或者多个的组合。
12.权利要求1-11中任一项的一组标签的用途，其用于标记SEQID NO：96-106所示的PCR引物。
13.一种标签引物组，其包含11种标签引物，所述标签引物的序列包含标签序列和PCR引物序列，并且所述标签序列连接至所述PCR引物序列的5′端，其中
1)所述标签序列选自SEQ ID NO：1-95，并且所述11种标签引物的每一种的标签序列都相同，和
2)所述11种标签引物的PCR引物序列分别如SEQ ID NO：96-106所示。
14.权利要求13的标签引物组，其中所述标签序列通过连接序列连接至所述PCR引物序列的5′端。
15.权利要求13或14的标签引物组的用途，其用于制备试剂盒，所述试剂盒用于HPV测序、检测或分型。
16.一种标签引物套组，其包含至少10个权利要求13或14的标签引物组，其中，每个标签引物组所使用的标签序列互不相同。
17.权利要求16的标签引物套组，其包含至少20个权利要求13或14的标签引物组。
18.权利要求16的标签引物套组，其包含至少30个权利要求13或14的标签引物组。
19.权利要求16的标签引物套组，其包含至少40个权利要求13或14的标签引物组。
20.权利要求16的标签引物套组，其包含至少50个权利要求13或14的标签引物组。
21.权利要求16的标签引物套组，其包含至少60个权利要求13或14的标签引物组。
22.权利要求16的标签引物套组，其包含至少70个权利要求13或14的标签引物组。
23.权利要求16的标签引物套组，其包含至少80个权利要求13或14的标签引物组。
24.权利要求16的标签引物套组，其包含至少90个权利要求13或14的标签引物组。
25.权利要求16的标签引物套组，其包含95个权利要求13或14的标签引物组。
26.权利要求16的标签引物套组，所使用的标签序列至少包括如下的标签序列：如SEQ ID NO：1-10，或SEQ ID NO：11-20，或SEQID NO：21-20，或SEQ ID NO：31-40，或SEQ ID NO：41-50，或SEQID NO：51-60，或SEQ ID NO：61-70，或SEQ ID NO：71-80，或SEQID NO：81-90，或SEQ ID NO：91-95所示的标签序列，或者它们任何两个或者多个的组合。
27.权利要求16-26中任一项的标签引物套组的用途，其用于制备诊断试剂盒，所述试剂盒用于高通量HPV测序、检测或分型。
28.一种试剂盒，其包含权利要求13或14的标签引物组或权利要求16-26中任一项的标签引物套组。
29.权利要求28的试剂盒，其还包含至少1种接头。
30.权利要求28的试剂盒，其还包含至少2种接头。
31.权利要求28的试剂盒，其还包含至少10种接头。
32.权利要求28的试剂盒，其还包含至少20种接头。
33.权利要求28的试剂盒，其还包含至少30种接头。
34.权利要求28的试剂盒，其还包含至少40种接头。
35.权利要求28的试剂盒，其还包含至少50种接头。
36.权利要求28的试剂盒，其还包含至少100种接头。
37.权利要求28的试剂盒，其还包含至少200种接头。
38.权利要求29-37中任一项的试剂盒，其中所述接头适用于Solexa测序法。
39.权利要求29-37中任一项的试剂盒，其中所述接头适用于构建测序文库。
40.权利要求29的试剂盒，其中所述接头的序列选自SEQ ID NO：121-132。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102137369A CN101921748B (zh) | 2010-06-30 | 2010-06-30 | 用于高通量检测人类乳头瘤病毒的dna分子标签 |
| PCT/CN2010/001812 WO2012000149A1 (zh) | 2010-06-30 | 2010-11-12 | 用于高通量检测人类乳头瘤病毒的dna分子标签 |
| HK11103419.4A HK1149292A1 (en) | 2010-06-30 | 2011-04-05 | Dna molecular indexes for high-throughput detecting human papilloma virus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102137369A CN101921748B (zh) | 2010-06-30 | 2010-06-30 | 用于高通量检测人类乳头瘤病毒的dna分子标签 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101921748A CN101921748A (zh) | 2010-12-22 |
| CN101921748B true CN101921748B (zh) | 2012-11-14 |
Family
ID=43336965
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010102137369A Active CN101921748B (zh) | 2010-06-30 | 2010-06-30 | 用于高通量检测人类乳头瘤病毒的dna分子标签 |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN101921748B (zh) |
| HK (1) | HK1149292A1 (zh) |
| WO (1) | WO2012000149A1 (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101921874B (zh) | 2010-06-30 | 2013-09-11 | 深圳华大基因科技有限公司 | 基于Solexa测序法的检测人类乳头瘤病毒的方法 |
| CN102690809B (zh) * | 2011-03-24 | 2013-12-04 | 深圳华大基因科技服务有限公司 | Dna标签及其在构建和测序配对末端标签文库中的应用 |
| CN102978206A (zh) * | 2012-11-27 | 2013-03-20 | 北京诺禾致源生物信息科技有限公司 | 一种应用于混合建库的高通量测序接头及其建库方法 |
| WO2016050111A1 (en) * | 2014-09-30 | 2016-04-07 | Bgi Shenzhen | Biomarkers for rheumatoid arthritis and usage thereof |
| CN105524983B (zh) * | 2014-09-30 | 2019-06-21 | 大连晶泰生物技术有限公司 | 基于高通量测序的标记和捕获多个样本的一个或多个特定基因的方法和试剂盒 |
| WO2016049927A1 (en) * | 2014-09-30 | 2016-04-07 | Bgi Shenzhen Co., Limited | Biomarkers for obesity related diseases |
| WO2016049918A1 (en) * | 2014-09-30 | 2016-04-07 | Bgi Shenzhen Co., Limited | Biomarkers for coronary artery disease |
| WO2016049937A1 (en) * | 2014-09-30 | 2016-04-07 | Bgi Shenzhen Co., Limited | Biomarkers for rheumatoid arthritis and usage therof |
| WO2016049930A1 (en) * | 2014-09-30 | 2016-04-07 | Bgi Shenzhen Co., Limited | Biomarkers for rheumatoid arthritis and usage therof |
| CN109136217B (zh) * | 2017-06-27 | 2022-02-22 | 深圳华大基因股份有限公司 | 一种测序文库构建的方法、建库试剂及其应用 |
| CN108130323A (zh) * | 2017-12-20 | 2018-06-08 | 奥明（杭州）基因科技有限公司 | 一种适用于植物转座酶可接近性染色质分析的建库方法 |
| CN114410773B (zh) * | 2022-01-27 | 2024-05-03 | 宁波大学 | 用于预测或诊断抑郁症复发的标志物组合及其应用 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101624632A (zh) * | 2008-07-08 | 2010-01-13 | 天津生物芯片技术有限责任公司 | 人乳头瘤病毒基因芯片、制备方法、应用及检测用试剂盒与应用 |
| CN101921874A (zh) * | 2010-06-30 | 2010-12-22 | 深圳华大基因科技有限公司 | 基于Solexa测序法的检测人类乳头瘤病毒的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101250593B (zh) * | 2008-02-03 | 2011-01-26 | 广州益善生物技术有限公司 | 乳头瘤病毒检测与分型方法及其液相芯片 |
| CN101435002B (zh) * | 2008-12-12 | 2012-01-18 | 深圳华大基因科技有限公司 | 一种检测人类乳头瘤病毒基因型的方法 |
-
2010
- 2010-06-30 CN CN2010102137369A patent/CN101921748B/zh active Active
- 2010-11-12 WO PCT/CN2010/001812 patent/WO2012000149A1/zh active Application Filing
-
2011
- 2011-04-05 HK HK11103419.4A patent/HK1149292A1/xx unknown
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101624632A (zh) * | 2008-07-08 | 2010-01-13 | 天津生物芯片技术有限责任公司 | 人乳头瘤病毒基因芯片、制备方法、应用及检测用试剂盒与应用 |
| CN101921874A (zh) * | 2010-06-30 | 2010-12-22 | 深圳华大基因科技有限公司 | 基于Solexa测序法的检测人类乳头瘤病毒的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2012000149A1 (zh) | 2012-01-05 |
| CN101921748A (zh) | 2010-12-22 |
| HK1149292A1 (en) | 2011-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101921874B (zh) | 基于Solexa测序法的检测人类乳头瘤病毒的方法 | |
| CN101921748B (zh) | 用于高通量检测人类乳头瘤病毒的dna分子标签 | |
| CN101435002B (zh) | 一种检测人类乳头瘤病毒基因型的方法 | |
| US20090029346A1 (en) | Detection of human papilloma virus | |
| Dalstein et al. | Analytical evaluation of the PapilloCheck test, a new commercial DNA chip for detection and genotyping of human papillomavirus | |
| Schmitz et al. | Quantitative multiplex PCR assay for the detection of the seven clinically most relevant high-risk HPV types | |
| Militello et al. | Accurate human papillomavirus genotyping by 454 pyrosequencing | |
| Doorn et al. | Molecular detection and genotyping of human papillomavirus | |
| JP2014501535A (ja) | 高リスクヒトパピローマウイルスのジェノタイピング及び定量化のための材料及び方法 | |
| Wallace et al. | Facile, comprehensive, high-throughput genotyping of human genital papillomaviruses using spectrally addressable liquid bead microarrays | |
| CN104560962B (zh) | 用于检测人乳头瘤病毒的荧光pcr试剂盒及其专用引物组 | |
| Lin et al. | A simple method for the detection and genotyping of high-risk human papillomavirus using seminested polymerase chain reaction and reverse hybridization | |
| CN101942440B (zh) | 检测和定型食管hpv病毒的引物和方法 | |
| CN102367490A (zh) | 一种检测病毒的方法 | |
| CN101624632B (zh) | 人乳头瘤病毒基因芯片、制备方法、应用及检测用试剂盒与应用 | |
| Gharizadeh et al. | Sentinel-base DNA genotyping using multiple sequencing primers for high-risk human papillomaviruses | |
| CN106048081A (zh) | 一种hpv分型检测引物及其检测方法和应用 | |
| Li et al. | Detection of human papillomavirus genotypes with liquid bead microarray in cervical lesions of northern Chinese patients | |
| Jones et al. | Comparison of the PapilloCheck® DNA micro-array Human Papillomavirus detection assay with Hybrid Capture II and PCR-enzyme immunoassay using the GP5/6+ primer set | |
| CN103361442B (zh) | 一种用于检测hpv基因型别的方法及试剂盒 | |
| Cañadas et al. | Comparison of the f-HPV typing™ and Hybrid Capture II® assays for detection of high-risk HPV genotypes in cervical samples | |
| Larsson et al. | HPV Genotyping from the high risk mRNA Aptima assay-a direct approach using DNA from Aptima sample tubes | |
| KR20090129639A (ko) | 인유두종 바이러스 진단용 프로브 및 그 유전자형 확인방법 | |
| Travasso et al. | Human papillomavirus genotyping by multiplex pyrosequencing in cervical cancer patients from India | |
| CN112921033A (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 | ||
| REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1149292 Country of ref document: HK |
|
| ASS | Succession or assignment of patent right |
Owner name: BGI-SHENZHEN CO., LTD. Effective date: 20120327 Owner name: SHANGHAI BGI TECHNOLOGY CO., LTD. Free format text: FORMER OWNER: BGI-SHENZHEN CO., LTD. Effective date: 20120327 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| C53 | Correction of patent of invention or patent application | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yi Xin Inventor after: Liu Tao Inventor after: Xu Jiajia Inventor after: Yang Xiaonan Inventor before: Yi Xin Inventor before: Liu Tao Inventor before: Xu Jiajia |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: YI XIN LIU TAO XU JIAJIA TO: YI XIN LIU TAO XU JIAJIA YANG XIAONAN Free format text: CORRECT: ADDRESS; FROM: 518083 SHENZHEN, GUANGDONG PROVINCE TO: 201318 PUDONG NEW AREA, SHANGHAI |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20120327 Address after: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Applicant after: SHENZHEN BGI Co.,Ltd. Co-applicant after: BGI SHENZHEN Co.,Ltd. Address before: Beishan Industrial Zone Building in Yantian District of Shenzhen city of Guangdong Province in 518083 Applicant before: BGI SHENZHEN Co.,Ltd. |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1149292 Country of ref document: HK |
|
| ASS | Succession or assignment of patent right |
Owner name: SHENZHEN BEIJING GENOMICS INSTITUTE HEALTH TECHNOL Free format text: FORMER OWNER: BGI-SHENZHEN CO., LTD. Effective date: 20131115 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20131115 Address after: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Patentee after: SHENZHEN BGI Co.,Ltd. Patentee after: BGI HEALTH SERVICE Co.,Ltd. Address before: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Patentee before: SHENZHEN BGI Co.,Ltd. Patentee before: BGI SHENZHEN Co.,Ltd. |
|
| C56 | Change in the name or address of the patentee | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Patentee after: SHENZHEN BGI Co.,Ltd. Patentee after: BGI SHENZHEN Co.,Ltd. Address before: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Patentee before: SHENZHEN BGI Co.,Ltd. Patentee before: BGI DIAGNOSIS Co.,Ltd. Address after: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Patentee after: SHENZHEN BGI Co.,Ltd. Patentee after: BGI DIAGNOSIS Co.,Ltd. Address before: 201318, Pudong New Area, Shanghai new highway 3399 Lane 26 building first, 7, 8 Patentee before: SHENZHEN BGI Co.,Ltd. Patentee before: BGI HEALTH SERVICE Co.,Ltd. |