CN104264231A - 构建测序文库的方法及其应用 - Google Patents
构建测序文库的方法及其应用 Download PDFInfo
- Publication number
- CN104264231A CN104264231A CN201410519353.2A CN201410519353A CN104264231A CN 104264231 A CN104264231 A CN 104264231A CN 201410519353 A CN201410519353 A CN 201410519353A CN 104264231 A CN104264231 A CN 104264231A
- Authority
- CN
- China
- Prior art keywords
- sequencing data
- sequencing
- chain
- sequence
- optionally
- 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.)
- Granted
Links
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
公开了构建测序文库的方法及其应用，该方法包括：(a)在双链DNA片段的两端分别连接接头，以便获得连接产物；(b)将所述连接产物裂解为单链DNA片段；(c)利用探针对所述单链DNA片段进行筛选；(d)利用第一引物对所述单链DNA片段进行链延伸反应，以便获得链延伸产物；(e)对所述链延伸产物进行扩增，以便获得扩增产物，所述扩增产物构成所述测序文库。还公开了测序方法、确定核酸序列的方法、构建测序文库的装置、测序设备以及确定核酸序列的系统。
Description
技术领域
本发明涉及生物医学领域。具体而言，本发明涉及构建测序文库的方法、测序方法、确定核酸序列的方法、构建测序文库的装置、测序设备以及确定核酸序列的系统。
背景技术
高通量测序日益被关注，但是目前高通量测序用于低频率突变的检测仍有待改进。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此，根据本发明的实施例，本发明提出了用于构建测序文库的方法以及检测低频率突变的手段。
在本发明的第一方面，本发明提出了一种构建测序文库的方法。根据本发明的实施例，该方法包括：(a)在双链DNA片段的两端分别连接接头，以便获得连接产物，其中，所述接头包括第一链和第二链，所述第一链和第二链部分匹配并且所述第一链包含第一标签序列，以便所述接头上限定出双链区和两个单链尾部，所述两个单链尾部之一的序列中包含第一标签；(b)将所述连接产物裂解为单链DNA片段；(c)利用探针对所述单链DNA片段进行筛选，其中，所述探针特异性识别预定区域，其中，所述预定区域包括下列之一：(1)表1所示基因的至少之一；(2)(1)的CDS区域；以及(3)(2)的上下游至少10bp的区域；(d)利用第一引物对所述单链DNA片段进行链延伸反应，以便获得链延伸产物，其中，所述第一引物包括第二标签序列，并且所述第一引物适于与所述接头的第一链形成双链结构，只是所述第一标签序列与所述第二标签序列之间存在错配；(e)对所述链延伸产物进行扩增，以便获得扩增产物，所述扩增产物构成所述测序文库，其中，所述扩增采用适于同时扩增所述第一标签序列和所述第二标签序列的引物。。
由此，利用根据本发明实施例的构建测序文库的方法，能够有效地构建测序文库，同时，所构建的测序文库中，针对相同的双链DNA片段(在本文中也被称为“源序列”)的每条链，分别获得了具有第一标签序列和第二标签序列的扩增产物，由此，在后续测序结果的分析中，可以依据两种标签的测序结果进行互相校正，提高分析结果的可靠性。
根据本发明的实施例，所述双链DNA片段是通过下列步骤获得的：将核酸样本进行末端修复，以便获得经过修复的核酸样本；以及在所述核酸样本的5’末端添加碱基A，以便获得两端分别具有粘性末端碱基A的核酸样本，所述两端分别具有粘性末端碱基A的核酸样本构成所述双链DNA片段。由此，可以在后续操作中，方便地在所述双链DNA片段的两端添加接头。从而，提高了构建测序文库的效率。
根据本发明的实施例，所述核酸样本为人基因组DNA的至少一部分或游离核酸。根据本发明的实施例，所述人游离核酸是从患者的外周血提取的。根据本发明的实施例，所述患者患有妇科生殖道肿瘤(即妇科癌症)，所述妇科生殖道肿瘤为选自下列的至少之一：卵巢癌、子宫内膜癌和宫颈癌。由此，利用本发明实施例的方法，能够有效地对人类疾病患者的基因突变进行有效的分析，进而能够有效用于常见肿瘤的早诊、个体化用药、以及术后监控等。
根据本发明的实施例，所述人基因组DNA的至少一部分是通过对人基因组DNA进行随机打断而获得的。由此，可以在后续操作中，方便地在所述双链DNA片段的两端添加接头。从而，提高了构建测序文库的效率。
根据本发明的实施例，所述接头具有3’碱基T粘性末端。由此，可以在后续操作中，方便地在所述双链DNA片段的两端添加接头。从而，提高了构建测序文库的效率。
根据本发明的实施例，所述单链DNA片段是通过将所述连接产物进行变性处理获得的。由此，可以快速有效的获得单链DNA片段。根据本发明的一些实施例，所述变性处理可以为热变性处理或碱变性处理。
根据本发明的实施例，所述探针是以芯片的形式提供的。由此，可以提高探针筛选的效率。
根据本发明的实施例，在存在UDG酶/FPG酶时，进行所述链延伸反应。由此，可以有效地对存在损伤的DNA在链延伸过程中进行修复，减少假阳性的产生，提高构建测序文库的质量。
根据本发明的实施例，所述第一标签序列和所述第二标签序列分别独立地长度为4～10nt。根据本发明的实施例，所述第一标签序列和所述第二标签序列的长度均为8nt。根据本发明的实施例，所述第一标签序列和所述第二标签序列之间存在至少2nt的错配。发明人惊奇地发现，采用如此设置，能够有效地提高在后续分析中，利用第一标签序列和第二标签序列进行校正的效率。
根据本发明的实施例，所述接头的第一链具有SEQ ID NO：1所示的序列，所述接头的第二链具有SEQ ID NO：2所示的序列，所述第一标签具有SEQ ID NO：3-6中任一项所示的序列，所述第二标签具有SEQ ID NO：7-10中至少之一所示的序列，所述第一引物具有SEQ ID NO：11所示的序列，所述适于同时扩增所述第一标签序列和所述第二标签序列的引物具有SEQ ID NO：12和SEQ ID NO：13所示的序列。
其中，接头的第一链的序列中“XXXXXXXX”表示第一标签序列，第一引物中序列中的“XXXXXXXX”表示第二标签序列。
根据本发明的实施例，标签包括但不限于上述所述的4对，可以根据需要涉及多对标签以用于多样品的同时检测。
在本发明的第二方面，本发明提出了一种测序方法，该方法包括：根据前面所述的方法构建测序文库；对所述测序文库进行测序。
根据本发明的实施例，在Hiseq2000或Hiseq2500上进行所述测序。由此，可以有效地提高测序的效率。另外，前面关于构建测序文库的方法所描述的特征和优点，同样适用该测序方法，在此不再赘述。
在本发明的第三方面，本发明提出了一种确定核酸序列的方法，该方法包括：针对核酸样本，根据权利要求前面所述的方法进行测序，以便获得由多个测序数据构成的测序结果；基于所述测序结果，构建至少一个测序数据子集，其中，每个测序数据子集中的所有测序数据均对应核酸样本上相同的源序列；针对每一个测序数据子集，分别确定与所述第一标签序列对应的测序数据为正链测序数据，与所述第二标签序列对应的测序数据为负链测序数据；针对每一个测序数据子集，分别基于所述正链测序数据和所述负链测序数据，对测序数据进行校正，以便确定经过校正的测序数据；以及基于所述经过校正的测序数据，确定所述核酸样本的序列。由此，可以有效地基于正链测序数据和负链测序数据进行校正，提高分析结果的可靠性。
根据本发明的实施例，所述测序为双末端测序，所述测序结果由多对成对的测序数据构成。
根据本发明的实施例，基于所述测序结果，构建至少一个测序数据子集是通过下列步骤进行的：针对所述多对成对的测序数据的每一对，确定成对测序数据索引，所述成对测序数据索引由成对的测序数据的每一个的最初N个碱基构成，其中，N为10～20之间的整数；基于所述成对测序数据索引，构建至少一个初步测序数据子集，其中，所述初步测序数据子集中的每一个测序数据均具有相同的成对测序数据索引；以及基于所述初步测序数据子集中测序数据之间的汉明距离，对所述至少一个初步测序数据子集进行细分，以便获得多个所述测序数据子集。
根据本发明的实施例，N为12。
根据本发明的实施例，在所述多个测序数据子集的每一个中，任意两对成对测序数据的汉明距离不超过20。
根据本发明的实施例，在所述多个测序数据子集的每一个中，正链测序数据和负链测序数据分别为至少两个。
根据本发明的实施例，基于所述正链测序数据和所述负链测序数据，确定经过校正的测序数据是基于下列原则进行的：经过校正的测序数据中的每一个碱基同时获得至少50％正链测序数据和至少50％负链测序数据的支持。
根据本发明的实施例，经过校正的测序数据中的每一个碱基同时获得至少80％正链测序数据和至少80％负链测序数据的支持。
根据本发明的实施例，进一步包括：将所述经过校正的测序数据比对至参考序列上，并删除比对质量小于30的测序数据。
根据本发明的实施例，进一步包括：基于所述核酸样本的序列，进行SNV分析或Indel分析。
在本发明的第四方面，本发明提出了一种构建测序文库的装置。根据本发明的实施例，该装置包括：连接单元，用于在双链DNA片段的两端分别连接接头，以便获得连接产物，其中，所述接头包括第一链和第二链，所述第一链和第二链部分匹配并且所述第一链包含第一标签序列，以便所述接头上限定出双链区和两个单链尾部，所述两个单链尾部之一的序列中包含第一标签；裂解单元，用于将所述连接产物裂解为单链DNA片段；筛选单元，用于在进行所述链延伸之前，利用探针对所述单链DNA片段进行筛选，其中，所述探针特异性识别预定区域，其中，所述预定区域包括下列之一：(1)表1所示基因的至少之一；(2)(1)的CDS区域；以及(3)(2)的上下游至少10bp的区域；链延伸单元，用于利用第一引物对所述单链DNA片段进行链延伸反应，以便获得链延伸产物，其中，所述第一引物包括第二标签序列，并且所述第一引物适于与所述接头的第一链形成双链结构，只是所述第一标签序列与所述第二标签序列之间存在错配；扩增单元，用于对所述链延伸产物进行扩增，以便获得扩增产物，所述扩增产物构成所述测序文库，其中，所述扩增采用适于同时扩增所述第一标签序列和所述第二标签序列的引物。
根据本发明的实施例，上述装置能够有效地实施前面所描述的构建测序文库的方法，能够有效地构建测序文库，同时，所构建的测序文库中，针对相同的双链DNA片段(在本文中也被称为“源序列”)的每条链，分别获得了具有第一标签序列和第二标签序列的扩增产物，由此，在后续测序结果的分析中，可以依据两种标签的测序结果进行互相校正，提高分析结果的可靠性。
根据本发明的实施例，进一步包括：末端修复单元，用于将核酸样本进行末端修复，以便获得经过修复的核酸样本；以及末端修饰单元，用于在所述核酸样本的5’末端添加碱基A，以便获得两端分别具有粘性末端碱基A的核酸样本，所述两端分别具有粘性末端碱基A的核酸样本构成所述双链DNA片段。
根据本发明的实施例，所述探针是以芯片的形式提供的。
根据本发明的实施例，在存在UDG酶/FPG酶时，进行所述链延伸反应。由此，可以有效地对存在损伤的DNA在链延伸过程中进行修复，减少假阳性的产生，提高构建测序文库的质量。
根据本发明的实施例，所述第一标签序列和所述第二标签序列分别独立地长度为4～10nt。
根据本发明的实施例，所述第一标签序列和所述第二标签序列的长度均为8nt。
根据本发明的实施例，所述第一标签序列和所述第二标签序列之间存在至少2nt的错配。
根据本发明的实施例，所述接头的第一链具有SEQ ID NO：1所示的序列，所述接头的第二链具有SEQ ID NO：2所示的序列，所述第一标签具有SEQ ID NO：3-6中任一项所示的序列，所述第二标签具有SEQ ID NO：7-10中至少之一所示的序列，所述第一引物具有SEQ ID NO：11所示的序列，所述适于同时扩增所述第一标签序列和所述第二标签序列的引物具有SEQ ID NO：12和SEQ ID NO：13所示的序列。
根据本发明的实施例，标签包括但不限于上述所述的4对，可以根据需要涉及多对标签以用于多样品的同时检测。
本领域技术人员能够理解的是，前面针对构建测序文库的方法所描述的特征和优点，同样适用于该构建测序文库的装置，在此不再赘述。
在本发明的第五方面，本发明提出了一种测序设备。根据本发明的实施例，该测序设备包括：根据前面所述的构建测序文库的装置；测序装置，用于对所述测序文库进行测序。
由此，可以有效地提高测序的效率。另外，前面关于构建测序文库的方法和装置所描述的特征和优点，同样适用该测序设备，在此不再赘述。
根据本发明的实施例，所述测序装置为Hiseq2000或Hiseq2500。
在本发明的第六方面，本发明提出了一种确定核酸序列的系统。根据本发明的实施例，该系统包括：前面所述的测序设备，用于针对核酸样本进行测序，以便获得由多个测序数据构成的测序结果；测序数据子集构建设备，用于基于所述测序结果，构建至少一个测序数据子集，其中，每个测序数据子集中的所有测序数据均对应核酸样本上相同的源序列；测序数据分类设备，用于针对每一个测序数据子集，分别确定与所述第一标签序列对应的测序数据为正链测序数据，与所述第二标签序列对应的测序数据为负链测序数据；测序数据校正设备，用于针对每一个测序数据子集，分别基于所述正链测序数据和所述负链测序数据，对测序数据进行校正，以便确定经过校正的测序数据；以及序列确定设备，用于基于所述经过校正的测序数据，确定所述核酸样本的序列。由此，利用根据本发明实施例的确定核酸序列的系统，能够有效地实施前面确定核酸序列的方法。从而可以有效地基于正链测序数据和负链测序数据进行校正，提高分析结果的可靠性。
根据本发明的实施例，所述测序为双末端测序，所述测序结果由多对成对的测序数据构成。
根据本发明的实施例，测序数据子集构建设备包括：测序数据索引确定设备，用于针对所述多对成对的测序数据的每一对，确定成对测序数据索引，所述成对测序数据索引由成对的测序数据的每一个的最初N个碱基构成，其中，N为10～20之间的整数；初步筛选装置，用于基于所述成对测序数据索引，构建至少一个初步测序数据子集，其中，所述初步测序数据子集中的每一个测序数据均具有相同的成对测序数据索引；以及二次筛选装置，用于基于所述初步测序数据子集中测序数据之间的汉明距离，对所述至少一个初步测序数据子集进行细分，以便获得多个所述测序数据子集。
根据本发明的实施例，N为12。
根据本发明的实施例，在所述多个测序数据子集的每一个中，任意两对成对测序数据的汉明距离不超过20。
根据本发明的实施例，在所述多个测序数据子集的每一个中，正链测序数据和负链测序数据分别为至少两个。
根据本发明的实施例，基于所述正链测序数据和所述负链测序数据，确定经过校正的测序数据是基于下列原则进行的：经过校正的测序数据中的每一个碱基同时获得至少50％正链测序数据和至少50％负链测序数据的支持。
根据本发明的实施例，经过校正的测序数据中的每一个碱基同时获得至少80％正链测序数据和至少80％负链测序数据的支持。
根据本发明的实施例，进一步包括：将所述经过校正的测序数据比对至参考序列上，并删除比对质量小于30的测序数据。
根据本发明的实施例，进一步包括序列分析装置，所述序列分析装置用于基于所述核酸样本的序列，进行SNV分析或Indel分析。
本领域技术人员可以理解的是，前面关于确定核酸序列的方法所描述的优点和特征同样适用该确定核酸序列的系统，在此不再赘述。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1显示了根据本发明的实施例的构建测序文库的方法的流程图；
图2显示了根据本发明的一个实施例相同索引reads簇的分析结果；以及
图3显示了根据本发明的一个实施例的突变谱的分析结果。
具体实施方式
下面通过具体的实施例，对本发明进行说明，需要说明的是这些实施例仅仅是为了说明目的，而不能以任何方式解释成对本发明的限制。
一般方法
除非特别说明，在下面的实施例中，按照下面的一般方法进行：
一、设计探针
根据人类基因组HG19，调取相关基因的外显子序列，考虑到捕获区域的大小及成本，最终的芯片只涉及了上述基因的CDS区域，并对CDS区域前后延伸了20bp。芯片上覆盖有丰富的捕获探针，探针覆盖区域达98％，可以从复杂的基因组中富集目标DNA片段，在同一张芯片上以高特异性和高覆盖率捕获基因组区域。
二、构建测序文库和测序
参照图1，构建文库和测序的步骤如下：
1.抽取患者5ml外周血，离心分离血浆和白细胞，分别对血浆样本和白细胞样本提取DNA，白细胞提取的DNA之后将作为对照用于体细胞突变的检出。
2.血浆中提取出的游离循环DNA平均在170BP，之后直接按照常规建库方法进行3步酶促反应：末端修复，加“A”和连接特殊处理的测序接头(该接头上带有8BP的标签，被命名为index1，其不仅具有区别不同样品的功能，也将被用于之后正链的标记)。
3.获得的连接产物，进行WCNpan芯片杂交捕获，其洗脱的单链模板产物之后经过1轮1个循环的带有index2标记的引物扩增，使得反链被标记。同时在PCR过程中加入UDG/FPG酶进行孵育，以消除模板链中带有的DNA损伤，减少假阳性的产生。
4.正反链双index标记完成的产物，经过纯化后，进行第二轮PCR富集，完成文库的制备。
5.测序方法采用Hiseq 2000或Hiseq2500，根据测序量的不同和样本数，可以灵活选择合适的测序平台。
具体步骤包括：
1.cfDNA的提取
取5ml外周血分离出的血浆约2-3ml，按照QIAamp Circulating Nucleic Acid Kit提取试剂说明书，进行血浆cfDNA的提取。Qubit(Invitrogen，the Quant-iT TM dsDNA HS Assay Kit)定量所提取的DNA，总量约为5～50ng。
2.样品文库的制备：
血浆中提取的cfDNA，之后按照KAPA LTP Library Preparation Kit建库说明书，进行3步酶促反应。
1)末端修复
之后，加入Agencourt AMPure XP reagent 120μL，进行磁珠纯化，最后回溶42μLddH2O，带磁珠进行下一步反应。
2)加A
之后加入PEG/NaCl SPRI溶液90μL，充分混合，进行磁珠纯化，最后回溶(35-接头)μLddH2O，带磁珠进行下一步反应。
3)接头连接
之后分别加入PEG/NaCl SPRI溶液50μL 2次，进行2次磁珠纯化，最后回溶25μL ddH2O。
3芯片杂交捕获
本发明中采用发明人设计的针对妇科癌症的早筛芯片，参照芯片制造商提供的说明书进行杂交捕获。最后洗脱回溶21μL ddH2O带杂交洗脱磁珠。
4.双index正反链标记和富集：
总共进行2轮PCR，PCR 1进行反链标记和模板DNA损伤修复，PCR2进行扩增富集，完成文库制备。
1)PCR1
PCR1程序：
先除去杂交洗脱磁珠，然后加入Agencourt AMPure XP reagent 40μL，进行磁珠纯化，最后回溶20ul ddH2O，带磁珠进行下一步反应。
2)PCR2
PCR2程序：
先除去上一步磁珠，然后重新加入Agencourt AMPure XP reagent 50μL，进行磁珠纯化，最后回溶25μL ddH2O，进行QC及上机。
三、测序结果分析
1，将paired reads(成对测序数据)的reads1的前12bp碱基和reads2的前12bp碱基(即断点序列)连接成24bp的一条短序列，并且以这24bp作为paired reads的索引，并根据其index标记正链和反链。
2，对索引进行外部排序，以达到将同一个DNA模板的拷贝聚集到一起的目的。
3，对聚集起来的拥有相同索引的reads进行中心聚类，根据其序列之间的汉明距离，将每个有相同索引的大簇聚集成若干个小簇，每个小簇中任意两对paired reads的汉明距离不超过10，以达到区分开拥有相同索引却来自不同DNA模板的reads的目的。
4，对步骤3中获得的同一个DNA模板的拷贝簇进行筛选，若正链和反链的reads数都达到2对以上，则进行后续分析。
5，对满足4中条件的簇进行纠错，并产生一对无错的新reads，对于DNA模板的每一个测序碱基，若某种碱基型在正链的reads中的一致率达到80％，且在反链reads中的一致率也达到80％，则记新reads的这个碱基为此碱基型，否则记为N，这样便得到了代表原始DNA模板序列的新reads。
6，将新reads用bwa mem算法重新比对到基因组上，筛除比对质量小于30的reads。
7，SNV分析:
1)根据6中得到的reads进行统计，得到捕获区域内每个位点的碱基型分布，与主流碱基型(比例大于15％的碱基型)不一致的碱基型既为突变碱基型。统计目标区域覆盖大小、平均测序深度，正反链互配率，低频突变率等。
2)利用CCDS、人类基因组数据库(NCBI36.3)、dbSNP(v130)信息对SNP进行注释，确定突变位点发生的基因、坐标、mRNA位点、氨基酸改变、SNP功能(错义突变/无义突变/可变剪切位点)、SIFT预测SNP影响蛋白功能预测等；
3)根据患者样品与对照样品信息的比对，Call Somatic Mutation。同时在候选的SNV中去除掉在dbSNP、HAPMAP、1000人类基因组、其他外显子测序项目中出现的SNP，以作为最后疾病相关的候选SNV。
8，INDEL分析：
1)根据6中得到的reads中含有indel的reads进行统计，得到所有的indel并选择有2条及以上reads支持的indel作为可靠的突变indel，
2)利用CCDS、人类基因组数据库(NCBI36.3)、dbSNP(v130)信息对Indel进行注释，确定突变位点发生的基因、坐标、mRNA位点、编码区域序列的改变、对氨基酸的影响、InDel功能(氨基酸插入/氨基酸缺失/移码突变)；
3)根据患者样品与对照样品信息的比对，Call Somatic Mutation。同时在候选的Indel中去除掉在dbSNP以及其他外显子测序项目中出现的Indel，以作为最后疾病相关的候选Indel。
实施例1：妇科生殖道肿瘤
一、芯片设计
基于TCGA，ICGC，COSMIC等数据库和相关文献参考，发明人设计出针对妇科生殖道肿瘤早筛的基因芯片WCNpan。WCNpan芯片包括了：妇科生殖道肿瘤(宫颈癌、子宫内膜癌、卵巢癌)相关的Driver Gene(驱动基因)，高频突变基因，以及癌症12条信号通路中重要基因等，共计42个基因，300KB。
芯片具体设计过程：根据人类基因组HG19，调取上述42个基因的外显子序列，考虑到捕获区域的大小及成本，最终的芯片只涉及了上述基因的CDS区域，并对CDS区域前后延伸了20bp，芯片总计300kb。该芯片上覆盖有丰富的捕获探针，探针覆盖区域达98％，可以从复杂的基因组中富集目标DNA片段，在同一张芯片上以高特异性和高覆盖率捕获约300KB的基因组区域。
基因列表详情见表1。
表1
| AFF3 | BRCA2 | FBXW7 | MED12 | PDE4DIP | STK11 |
| AKAP9 | CDK12 | FGFR2 | MLL2 | PIK3CA | TP53 |
| AKT1 | CDKN2A | FGFR3 | MLL3 | PIK3R1 | |
| APC | CREBBP | FOXL2 | MSH6 | PPP2R1A | |
| ARID1A | CSMD3 | GNAS | NF1 | PTEN | |
| BCOR | CTNNB1 | HRAS | NFE2L2 | RB1 | |
| BRAF | EGFR | KIT | NRAS | RNF213 | |
| BRCA1 | FAT3 | KRAS | NSD1 | RNF43 |
二、测序分析
对1例宫颈不典型增生患者按照以上方法的步骤进行分析，结果如下：
表2：测序数据统计结果
注释：正反链互配率：基于3条reads以上正反链均有的簇/3条reads上总的簇的比值，以评估可用数据中正反链互配情况；有效数据利用率：基于至少满足2+/2-簇的reads纠错后的个数与总测序reads数的比值；平均测序深度：基于有效数据纠错后，对目标区域碱基的平均覆盖情况。
相同索引reads簇的分析结果见图2，其中，横坐标表示表簇的duplication(dup)个数，纵坐标表示满足某一dup个数的簇的总的reads个数。由图2的结果可知，dup簇绝大部分在8左右，较大部分簇内能满足2正+2反的条件，最终数据有效利用率为5.14％，平均测序深度为：1153.6X。
突变谱分析结果见图3，其中，互补的突变类型对于来源于双链的分子(DNA)，理论突变频率基本相同。横坐标代表碱基突变的类型；纵坐标代表突变的个数。图3的结果显示：突变碱基类型分布平衡，其突变频率(Mutations per nucleotide)为：1.7×10-6。
变异检测列表详情(基于exon区和非同义突变统计)见表1。
表1
| 基因 | cHGVS | pHGVS | 突变类型 | 突变频率 |
| PIK3CA | c.2119G>A | p.Glu707Lys | 错义突变 | 5.60％ |
| TP53 | c.817C>T | p.Arg273Cys | 错义突变 | 2.40％ |
| CDKN2A | c.217A>C | p.Ser73Arg | 错义突变 | 1.80％ |
| NRAS | c.182A>G | p.Gln61Arg | 错义突变 | 1.40％ |
结果分析：依据TCGA，COSMIC，ClinVar,HMGD等相关数据库以及文献资料，在患者血浆中检测到了PIK3CAp.Glu707Lys,TP53 p.Arg273Cys等Driver mutation(驱动突变)，预示着患者具有较高的癌症风险率，建议患者到相关医疗机构进行更全面的检测和采取相关干预措施。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。另外，需要说明的是，本领域技术人员能够理解，在本发明所提出的方案中所包含的步骤顺序，本领域技术人员可以进行调整，这也将包括在本发明的范围内。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
Claims (10)
1.一种构建测序文库的方法，其特征在于，包括：
(a)在双链DNA片段的两端分别连接接头，以便获得连接产物，其中，所述接头包括第一链和第二链，所述第一链和第二链部分匹配并且所述第一链包含第一标签序列，以便所述接头上限定出双链区和两个单链尾部，所述两个单链尾部之一的序列中包含第一标签；
(b)将所述连接产物裂解为单链DNA片段；
(c)利用探针对所述单链DNA片段进行筛选，其中，所述探针特异性识别预定区域，其中，所述预定区域包括下列之一：
(1)表1所示基因的至少之一；
(2)(1)的CDS区域；以及
(3)(2)的上下游至少10bp的区域；
(d)利用第一引物对所述单链DNA片段进行链延伸反应，以便获得链延伸产物，其中，所述第一引物包括第二标签序列，并且所述第一引物适于与所述接头的第一链形成双链结构，只是所述第一标签序列与所述第二标签序列之间存在错配；
(e)对所述链延伸产物进行扩增，以便获得扩增产物，所述扩增产物构成所述测序文库，其中，所述扩增采用适于同时扩增所述第一标签序列和所述第二标签序列的引物。
2.根据权利要求1所述的方法，其特征在于，所述双链DNA片段是通过下列步骤获得的：
将核酸样本进行末端修复，以便获得经过修复的核酸样本；以及
在所述核酸样本的5’末端添加碱基A，以便获得两端分别具有粘性末端碱基A的核酸样本，所述两端分别具有粘性末端碱基A的核酸样本构成所述双链DNA片段，
任选地，所述核酸样本为人基因组DNA的至少一部分或游离核酸，
任选地，所述人游离核酸是从患者的外周血提取的，
任选地，所述患者患有妇科癌症，所述妇科癌症为选自下列的至少之一：
卵巢癌、子宫内膜癌、以及宫颈癌，
任选地，所述人基因组DNA的至少一部分是通过对人基因组DNA进行随机打断而获得的，
任选地，所述接头具有3’碱基T粘性末端，
任选地，所述单链DNA片段是通过将所述连接产物进行变性处理获得的，
任选地，所述探针是以芯片的形式提供的，
任选地，在存在UDG酶/FPG酶时，进行所述链延伸反应，
任选地，所述第一标签序列和所述第二标签序列分别独立地长度为4～10nt，优选8nt，
任选地，所述第一标签序列和所述第二标签序列的长度均为8nt，
任选地，所述第一标签序列和所述第二标签序列之间存在至少2nt的错配，
任选地，所述接头的第一链具有SEQ ID NO：1所示的序列，所述接头的第二链具有SEQ ID NO：2所示的序列，所述第一标签具有SEQ ID NO：3-6中至少之一所示的序列，所述第二标签具有SEQ ID NO：7-10中至少之一所示的序列，所述第一引物具有SEQ IDNO：11所示的序列，所述第二引物具有SEQ ID NO：12所示的序列，所述第三引物具有SEQ ID NO：13所示的序列。
3.一种测序方法，其特征在于，包括：
根据权利要求1或2所述的方法构建测序文库；
对所述测序文库进行测序，
任选地，在Hiseq2000或Hiseq2500上进行所述测序。
4.一种确定核酸序列的方法，其特征在于，包括：
针对核酸样本，根据权利要求3所述的方法进行测序，以便获得由多个测序数据构成的测序结果；
基于所述测序结果，构建至少一个测序数据子集，其中，每个测序数据子集中的所有测序数据均对应核酸样本上相同的源序列；
针对每一个测序数据子集，分别确定与所述第一标签序列对应的测序数据为正链测序数据，与所述第二标签序列对应的测序数据为负链测序数据；
针对每一个测序数据子集，分别基于所述正链测序数据和所述负链测序数据，对测序数据进行校正，以便确定经过校正的测序数据；以及
基于所述经过校正的测序数据，确定所述核酸样本的序列。
5.根据权利要求4所述的方法，其特征在于，所述测序为双末端测序，所述测序结果由多对成对的测序数据构成，
任选地，基于所述测序结果，构建至少一个测序数据子集是通过下列步骤进行的：
针对所述多对成对的测序数据的每一对，确定成对测序数据索引，所述成对测序数据索引由成对的测序数据的每一个的最初N个碱基构成，其中，N为10～20之间的整数；
基于所述成对测序数据索引，构建至少一个初步测序数据子集，其中，所述初步测序数据子集中的每一个测序数据均具有相同的成对测序数据索引；以及
基于所述初步测序数据子集中测序数据之间的汉明距离，对所述至少一个初步测序数据子集进行细分，以便获得多个所述测序数据子集，
任选地，N为12，
任选地，在所述多个测序数据子集的每一个中，任意两对成对测序数据的汉明距离不超过20，
任选地，在所述多个测序数据子集的每一个中，正链测序数据和负链测序数据分别为至少两个，
任选地，基于所述正链测序数据和所述负链测序数据，确定经过校正的测序数据是基于下列原则进行的：
经过校正的测序数据中的每一个碱基同时获得至少50％正链测序数据和至少50％负链测序数据的支持，
任选地，经过校正的测序数据中的每一个碱基同时获得至少80％正链测序数据和至少80％负链测序数据的支持，
任选地，进一步包括：
将所述经过校正的测序数据比对至参考序列上，并删除比对质量小于30的测序数据，
任选地，基于所述核酸样本的序列，进行SNV分析或Indel分析。
6.一种构建测序文库的装置，其特征在于，包括：
连接单元，用于在双链DNA片段的两端分别连接接头，以便获得连接产物，其中，所述接头包括第一链和第二链，所述第一链和第二链部分匹配并且所述第一链包含第一标签序列，以便所述接头上限定出双链区和两个单链尾部，所述两个单链尾部之一的序列中包含第一标签；
裂解单元，用于将所述连接产物裂解为单链DNA片段；
筛选单元，用于在进行所述链延伸之前，利用探针对所述单链DNA片段进行筛选，其中，所述探针特异性识别预定区域，其中，所述预定区域包括下列之一：
(1)表1所示基因的至少之一；
(2)(1)的CDS区域；以及
(3)(2)的上下游至少10bp的区域；
链延伸单元，用于利用第一引物对所述单链DNA片段进行链延伸反应，以便获得链延伸产物，其中，所述第一引物包括第二标签序列，并且所述第一引物适于与所述接头的第一链形成双链结构，只是所述第一标签序列与所述第二标签序列之间存在错配；
扩增单元，用于对所述链延伸产物进行扩增，以便获得扩增产物，所述扩增产物构成所述测序文库，其中，所述扩增采用第二引物以及第三引物，所述第二引物识别所述接头的第二链，所述第三引物被设置为适于同时扩增所述第一标签序列和所述第二标签序列。
7.根据权利要求6所述的装置，其特征在于，进一步包括：
末端修复单元，用于将核酸样本进行末端修复，以便获得经过修复的核酸样本；以及
末端修饰单元，用于在所述核酸样本的5’末端添加碱基A，以便获得两端分别具有粘性末端碱基A的核酸样本，所述两端分别具有粘性末端碱基A的核酸样本构成所述双链DNA片段，
任选地，所述探针是以芯片的形式提供的，
任选地，在存在UDG酶/FPG酶时，进行所述链延伸反应，
任选地，所述第一标签序列和所述第二标签序列分别独立地长度为4～10nt，
任选地，所述第一标签序列和所述第二标签序列的长度均为8nt，
任选地，所述第一标签序列和所述第二标签序列之间存在至少2nt的错配，
任选地，所述接头的第一链具有SEQ ID NO：1所示的序列，所述接头的第二链具有SEQ ID NO：2所示的序列，所述第一标签具有SEQ ID NO：3-6中至少之一所示的序列，所述第二标签具有SEQ ID NO：7-10中至少之一所示的序列，所述第一引物具有SEQ IDNO：11所示的序列，所述第二引物具有SEQ ID NO：12所示的序列，所述第三引物具有SEQ ID NO：13所示的序列。
8.一种测序设备，其特征在于，包括：
根据权利要求6或7所述的构建测序文库的装置；
测序装置，用于对所述测序文库进行测序，
任选地，所述测序装置为Hiseq2000或Hiseq2500。
9.一种确定核酸序列的系统，其特征在于，包括：
权利要求8所述的测序设备，用于针对核酸样本进行测序，以便获得由多个测序数据构成的测序结果；
测序数据子集构建设备，用于基于所述测序结果，构建至少一个测序数据子集，其中，每个测序数据子集中的所有测序数据均对应核酸样本上相同的源序列；
测序数据分类设备，用于针对每一个测序数据子集，分别确定与所述第一标签序列对应的测序数据为正链测序数据，与所述第二标签序列对应的测序数据为负链测序数据；
测序数据校正设备，用于针对每一个测序数据子集，分别基于所述正链测序数据和所述负链测序数据，对测序数据进行校正，以便确定经过校正的测序数据；以及
序列确定设备，用于基于所述经过校正的测序数据，确定所述核酸样本的序列。
10.根据权利要求9所述的系统，其特征在于，所述测序为双末端测序，所述测序结果由多对成对的测序数据构成，
任选地，测序数据子集构建设备包括：
测序数据索引确定设备，用于针对所述多对成对的测序数据的每一对，确定成对测序数据索引，所述成对测序数据索引由成对的测序数据的每一个的最初N个碱基构成，其中，N为10～20之间的整数；
初步筛选装置，用于基于所述成对测序数据索引，构建至少一个初步测序数据子集，其中，所述初步测序数据子集中的每一个测序数据均具有相同的成对测序数据索引；以及
二次筛选装置，用于基于所述初步测序数据子集中测序数据之间的汉明距离，对所述至少一个初步测序数据子集进行细分，以便获得多个所述测序数据子集，
任选地，N为12，
任选地，在所述多个测序数据子集的每一个中，任意两对成对测序数据的汉明距离不超过20，
任选地，在所述多个测序数据子集的每一个中，正链测序数据和负链测序数据分别为至少两个，
任选地，基于所述正链测序数据和所述负链测序数据，确定经过校正的测序数据是基于下列原则进行的：
经过校正的测序数据中的每一个碱基同时获得至少50％正链测序数据和至少50％负链测序数据的支持，
任选地，经过校正的测序数据中的每一个碱基同时获得至少80％正链测序数据和至少80％负链测序数据的支持，
任选地，进一步包括：
将所述经过校正的测序数据比对至参考序列上，并删除比对质量小于30的测序数据，
任选地，进一步包括序列分析装置，所述序列分析装置用于基于所述核酸样本的序列，进行SNV分析或Indel分析。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410519353.2A CN104264231B (zh) | 2014-09-30 | 2014-09-30 | 构建测序文库的方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410519353.2A CN104264231B (zh) | 2014-09-30 | 2014-09-30 | 构建测序文库的方法及其应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104264231A true CN104264231A (zh) | 2015-01-07 |
| CN104264231B CN104264231B (zh) | 2017-04-19 |
Family
ID=52155799
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410519353.2A Active CN104264231B (zh) | 2014-09-30 | 2014-09-30 | 构建测序文库的方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104264231B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104293940A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN105087789A (zh) * | 2015-08-10 | 2015-11-25 | 北京吉因加科技有限公司 | 一种检测血浆cfDNA中BCR和TCR免疫组库的方法 |
| CN105525009A (zh) * | 2016-01-26 | 2016-04-27 | 绍兴华因生物科技有限公司 | 检测杂合性brca1/2基因缺失的方法及所用引物 |
| CN107523563A (zh) * | 2017-09-08 | 2017-12-29 | 杭州和壹基因科技有限公司 | 一种用于循环肿瘤dna分析的生物信息处理方法 |
| CN107577921A (zh) * | 2017-08-25 | 2018-01-12 | 云壹生物技术（大连）有限公司 | 一种肿瘤靶向基因测序数据解析方法 |
| CN110079592A (zh) * | 2018-01-26 | 2019-08-02 | 厦门艾德生物医药科技股份有限公司 | 用于检测基因突变和已知、未知基因融合类型的高通量测序靶向捕获目标区域的探针和方法 |
| WO2019191900A1 (en) * | 2018-04-03 | 2019-10-10 | Burning Rock Biotech | Compositions and methods for preparing nucleic acid libraries |
| CN113518830A (zh) * | 2019-01-31 | 2021-10-19 | 新加坡科技研究局 | 对具有非天然碱基对的核酸进行测序的方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101967476A (zh) * | 2010-09-21 | 2011-02-09 | 深圳华大基因科技有限公司 | 一种基于接头连接的DNA PCR-Free标签文库构建方法 |
| CN102127818A (zh) * | 2010-12-15 | 2011-07-20 | 张康 | 利用孕妇外周血建立胎儿dna文库的方法 |
| CN102409048A (zh) * | 2010-09-21 | 2012-04-11 | 深圳华大基因科技有限公司 | 一种基于高通量测序的dna标签文库构建方法 |
| CN102560688A (zh) * | 2010-12-15 | 2012-07-11 | 深圳华大基因科技有限公司 | 一种新的基于illumina测序平台的文库构建方法 |
| CN103103624A (zh) * | 2011-11-15 | 2013-05-15 | 深圳华大基因科技有限公司 | 高通量测序文库的构建方法及其应用 |
| CN103806111A (zh) * | 2012-11-15 | 2014-05-21 | 深圳华大基因科技有限公司 | 高通量测序文库的构建方法及其应用 |
-
2014
- 2014-09-30 CN CN201410519353.2A patent/CN104264231B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101967476A (zh) * | 2010-09-21 | 2011-02-09 | 深圳华大基因科技有限公司 | 一种基于接头连接的DNA PCR-Free标签文库构建方法 |
| CN102409048A (zh) * | 2010-09-21 | 2012-04-11 | 深圳华大基因科技有限公司 | 一种基于高通量测序的dna标签文库构建方法 |
| CN102127818A (zh) * | 2010-12-15 | 2011-07-20 | 张康 | 利用孕妇外周血建立胎儿dna文库的方法 |
| CN102560688A (zh) * | 2010-12-15 | 2012-07-11 | 深圳华大基因科技有限公司 | 一种新的基于illumina测序平台的文库构建方法 |
| CN103103624A (zh) * | 2011-11-15 | 2013-05-15 | 深圳华大基因科技有限公司 | 高通量测序文库的构建方法及其应用 |
| CN103806111A (zh) * | 2012-11-15 | 2014-05-21 | 深圳华大基因科技有限公司 | 高通量测序文库的构建方法及其应用 |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104293940A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN105087789A (zh) * | 2015-08-10 | 2015-11-25 | 北京吉因加科技有限公司 | 一种检测血浆cfDNA中BCR和TCR免疫组库的方法 |
| CN105087789B (zh) * | 2015-08-10 | 2018-08-03 | 北京吉因加科技有限公司 | 一种检测血浆cfDNA中BCR和TCR免疫组库的方法 |
| CN105525009B (zh) * | 2016-01-26 | 2018-10-16 | 绍兴华因生物科技有限公司 | 检测杂合性brca1/2基因缺失的方法及所用引物 |
| CN105525009A (zh) * | 2016-01-26 | 2016-04-27 | 绍兴华因生物科技有限公司 | 检测杂合性brca1/2基因缺失的方法及所用引物 |
| CN107577921A (zh) * | 2017-08-25 | 2018-01-12 | 云壹生物技术（大连）有限公司 | 一种肿瘤靶向基因测序数据解析方法 |
| CN107523563A (zh) * | 2017-09-08 | 2017-12-29 | 杭州和壹基因科技有限公司 | 一种用于循环肿瘤dna分析的生物信息处理方法 |
| CN110079592A (zh) * | 2018-01-26 | 2019-08-02 | 厦门艾德生物医药科技股份有限公司 | 用于检测基因突变和已知、未知基因融合类型的高通量测序靶向捕获目标区域的探针和方法 |
| CN110079592B (zh) * | 2018-01-26 | 2021-02-12 | 厦门艾德生物医药科技股份有限公司 | 用于检测基因突变和已知、未知基因融合类型的高通量测序靶向捕获目标区域的探针和方法 |
| WO2019191900A1 (en) * | 2018-04-03 | 2019-10-10 | Burning Rock Biotech | Compositions and methods for preparing nucleic acid libraries |
| WO2019192489A1 (en) * | 2018-04-03 | 2019-10-10 | Burning Rock Biotech | Compositions and methods for preparing nucleic acid libraries |
| CN113106145A (zh) * | 2018-04-03 | 2021-07-13 | 广州燃石医学检验所有限公司 | 用于制备核酸文库的组合物和方法 |
| CN113518830A (zh) * | 2019-01-31 | 2021-10-19 | 新加坡科技研究局 | 对具有非天然碱基对的核酸进行测序的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104264231B (zh) | 2017-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104264231A (zh) | 构建测序文库的方法及其应用 | |
| CN104293940A (zh) | 构建测序文库的方法及其应用 | |
| Chan et al. | Optical mapping reveals a higher level of genomic architecture of chained fusions in cancer | |
| AU2016293025A1 (en) | System and methodology for the analysis of genomic data obtained from a subject | |
| US12031186B2 (en) | Homologous recombination repair deficiency detection | |
| CN115679000B (zh) | 微小残留病灶的检测方法、装置、设备和存储介质 | |
| CN110093417B (zh) | 一种检测肿瘤单细胞体细胞突变的方法 | |
| CN104293941B (zh) | 构建测序文库的方法及其应用 | |
| JP2023523002A (ja) | 染色体近接実験における構造的変異検出 | |
| Kozarewa et al. | A modified method for whole exome resequencing from minimal amounts of starting DNA | |
| CN113748467A (zh) | 基于等位基因频率的功能丧失计算模型 | |
| JP2024056984A (ja) | エピジェネティック区画アッセイを較正するための方法、組成物およびシステム | |
| US20240141425A1 (en) | Correcting for deamination-induced sequence errors | |
| US20140364321A1 (en) | Method for analyzing DNA methylation based on MspJI cleavage | |
| CN102831331B (zh) | 基于酶切建库双末端测序的长度多态性标记的引物设计开发方法 | |
| Das et al. | Mapinsights: deep exploration of quality issues and error profiles in high-throughput sequence data | |
| JP2023524681A (ja) | 分配された核酸を使用した配列決定のための方法 | |
| CN114746560A (zh) | 改进甲基化多核苷酸结合的方法、组合物和系统 | |
| CN108304693B (zh) | 利用高通量测序数据分析基因融合的方法 | |
| Gorelick et al. | Respiratory complex and tissue lineage drive mutational patterns in the tumor mitochondrial genome | |
| US20130143746A1 (en) | Method for detecting gene region features based on inter-alu polymerase chain reaction | |
| WO2018219581A1 (en) | Method and system for nucleic acid sequencing | |
| CN117965734B (zh) | 一种用于检测硬纤维瘤的基因标志物、试剂盒、检测方法及应用 | |
| US20220223226A1 (en) | Methods for detecting and characterizing microsatellite instability with high throughput sequencing | |
| CN118006770A (zh) | 使用游离dna分子的正负链差异检测癌症 |
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 |