用于与BCR或ABL基因杂交的DNA探针库及采用其富集BCR-ABL基因片段的方法

摘要

本发明提供了一种用于与BCR或ABL基因杂交的DNA探针库，所述DNA探针库包括一个或多个能够与BCR或ABL基因杂交的DNA探针。本发明还提供了采用该DNA探针库富集BCR-ABL基因片段的方法。基于此，本发明还提供了一种检测是否存在BCR-ABL基因或其基因结构突变的方法。采用本发明的方法可以成数千倍地富集得到基因片段，并可以将其用于下一代测序技术以进行是否存在BCR-ABL基因或其基因结构突变的检测，包括单碱基突变、mRNA缺失或增多、mRNA结构颠换、mRNA剪接变化。

说明书

技术领域

本发明涉及基因检测领域，具体而言，本发明涉及一种BCR-ABL基因片段的富集和提取方法，所述方法可以精准地富集BCR-ABL基因片段，进而可以选择性地进行是否存在BCR-ABL基因以及其基因结构突变的检测。 

背景技术

DNA测序技术正处在天翻地覆的剧变中，其突出特点为同时对众多的位点进行观察分析（大规模平行），从而逐步实现测序通量的大幅增长，原始数据中每个碱基的测序成本的急剧下跌。基于此，以前高不可攀的奢侈性活动（如个人基因测序、宏基因组学研究），逐渐变得越来越切实可行。特别是随着科学的发展，由于传统的Sanger测序已经不能完全满足研究的需要，下一代测序技术（第二代测序技术，Next-generation sequencing）应运而生。 

下一代测序技术是同步化三磷酸核苷酸的洗脱方法和同步化的光学检测方法的结合。例如采用以鲁米那（Illumina）提供的仪器，以短的连续性的片段序列和测序阅读长度的形式，每周输出数亿以计的碱基的DNA序列。这是一种由DNA连接酶和聚合酶主导化学过程的第二代测序方法，DNA序列将会用生物信息学的方法拼接还原。在成本方面，下一代测序技术比第一代测序技术大体上降低了1000倍，并且还正以指数级的速度下降。目前下一代测序技术已经广泛应用于全基因组测序，但是在其它基因分析领域，例如特定基因结构突变的检测等方面尚未得到充分开发。 

基因检测是利用血液、其他体液或细胞对核酸进行检测的技术，其通过特定设备对被检测者核酸分子信息作检测，分析它所含有的各种基因情况，从而使人们能了解自己的基因信息，判断身体患有疾病的情况或风险，从而可适应性地选择疾病治疗方案，甚至预防疾病的发生。 

现有的基因突变检测方法还主要集中于以下3个技术： 

1）PCR突变检测 

基于PCR片段扩增和凝胶电泳分离的方法，从而分辨出野生型和突变性的差异。其缺点包括：其为单碱基突变类型检测，不能对mRNA进行定量，不能检测基因颠换；敏感性低；耗时长，单次只能检测一个基因突变；不能确定碱基的具体变化；无法进行高通量检测。 

2）Q-PCR（定量PCR）检测 

基于荧光和PCR片段扩增来对模板mRNA多少来进行定量。其缺点包括：不能检测单碱基类型突变；不能检测基因颠换；只能检测已知突变；耗时长，单次只能检测一个基因突变；不能确定碱基的具体变化。 

3）基因芯片技术 

用高精度技术将DNA片段打印在芯片上，然后通过DNA杂交结合特性来确定突变性质。其缺点包括：不能检测基因颠换；只能检测已知突变；成本高，通量较小；准确率低，通常需要重复2次以上才能确定结果。 

因此可知，目前在基因检测技术方面尚缺少更为全面、快速、简便、准确的检测方法，特别是在目前拥有下一代测序技术这种以高通量、低成本基因测序为特点的技术的情况下，如何基于下一代测序技术开发新型基因检测和筛查技术，是本领域研究人员广泛关注的问题。并且，如果应用下一代测序技术来检测基因的结构突变情况，如何获得能够满足下一代测序技术要求的基因检测样本，也是期待解决的一个技术问题。 

BCR-ABL基因常见于慢性粒细胞性白血病患者体内。慢性粒细胞性白血病是一种发生于造血干细胞的血液系统恶性克隆增生性疾病，部分患者的22号染色体和9号染色体的相互易位产生费城染色体，造成BCR与ABL基因融合，这种融合基因的产物蛋白产生酪氨酸蛋白激酶活动。研究表明，这种染色体易位现象与慢性粒细胞性白血病有高度的灵敏度和预测性，有95%的慢性粒细胞性白血病患者被检测有此染色体易位。此法也发现于大约20-30%的成人与2-10%的未成年患者急性白血病以及个别的急性骨髓性白血病案例中。 

目前慢性粒细胞白血病的治疗方法主要有酪氨酸激酶抑制剂，免疫疗法和造血干细胞移植。其中，抑制BCR-ABL激酶活性的格列卫（STI-571）是临床经常使用的标靶药物，在其之后还有新型抑制剂达沙替尼(Dasatinib)或尼罗替尼(Nilotinib)。伊马替尼是2-苯氨嘧啶衍生物，它可以选择性地阻断ATP与ABL激酶结合位点，使该酶失活，进而阻止了一系列的信号传导。 实验表明，伊马替尼不杀伤BCR-ABL-细胞，只杀伤BCR-ABL+细胞。伊马替尼通过抑制BCR-ABL活性，使得参与细胞周期、粘附骨架形成等生理过程的多种基因转录发生改变，引起含有BCR-ABL细胞分化的分化和凋亡。 

然而，格列卫并不能根治CML，因为它的作用原理是抑制BCR-ABL激酶活性，并不能抑制BCR-ABL的形成和表达。因此停药后易复发以及耐药性是两大存在问题。耐药的发生机制有：(1)ABL激酶区发生突变，改变了格列卫结合位点空间构象，导致格列卫不能结合；(2)BCR/ABL基因表达扩增或者融合蛋白过度表达，超过了格列卫的竞争结合能力。临床中80%的耐药是由于激酶区突变引起的，已报道的突变超过50种。部分患者在使用格列卫治疗前就携带有突变的克隆。约10%的耐药是由于BCR/ABL融合基因的过量表达导致的，可以用增加药物剂量来达到一定效果。 

目前发现的ABL激酶区突变主要集中于ATP结合区（第244～255氨基酸）、T315、M351及活化区四个区。部分突变患者可以对新型的激酶抑制剂达沙替尼(Dasatinib)或尼罗替尼(Nilotinib)有治疗反应。ATP结合区突变耐药性较强，如G250E、Y253H及E255K等。但是目前尤以T315I耐药程度最高，所有激酶抑制剂都没有效果。M351T、E355G、M244V等耐药程度不高的突变，可以提高剂量。因此对携有BCR-ABL融合基因的病人进行治疗前与治疗后的突变检测是十分重要的。 

因此，检测患者是否具有BCR-ABL基因以及检测BCR-ABL基因的突变，对判断疾病的发生发展、治疗药物的选择等方面均有重要意义。 

发明内容

针对上述技术问题，本发明人通过大量实验，开发了基于杂交选择而捕获BCR-ABL基因片段序列的方法，在受试者体内具有BCR-ABL基因的情况下，采用该方法可以获得成几千倍富集的BCR-ABL基因片段，该经富集的BCR-ABL基因片段样本可以选择性地应用于各种基因检测技术，特别是可以应用下一代测序技术进行基因突变、缺失、增加、和颠换等方面的检测。另外，通过富集后进行检测，还可鉴别受试者体内是否具有BCR-ABL基因。 

具体而言，本发明的技术方案如下： 

一方面，本发明提供了一种用于与BCR或ABL基因杂交的DNA探针库，所述DNA探针库包括一个或多个能够与BCR或ABL基因杂交的DNA探针，所述DNA探针包含以下序列： 

SEQ ID NO.1、SEQ ID NO.2、SEQ ID NO.3、SEQ ID NO.4、SEQ ID NO.5、SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8、SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11、SEQ ID NO.12、SEQ ID NO.13、SEQ ID NO.14、SEQ ID NO.15、SEQ ID NO.16、SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20、SEQ ID NO.21、SEQ ID NO.22、SEQ ID NO.23、SEQ ID NO.24、SEQ ID NO.25、SEQ ID NO.26、SEQ ID NO.27、SEQ ID NO.28、SEQ ID NO.29、SEQ ID NO.30、或SEQ ID NO.31； 

优选地，所述DNA探针的序列如以下序列所示： 

SEQ ID NO.1、SEQ ID NO.2、SEQ ID NO.3、SEQ ID NO.4、SEQ ID NO.5、SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8、SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11、SEQ ID NO.12、SEQ ID NO.13、SEQ ID NO.14、SEQ ID NO.15、SEQ ID NO.16、SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20、SEQ ID NO.21、SEQ ID NO.22、SEQ ID NO.23、SEQ ID NO.24、SEQ ID NO.25、SEQ ID NO.26、SEQ ID NO.27、SEQ ID NO.28、SEQ ID NO.29、SEQ ID NO.30、或SEQ ID NO.31。 

另一方面，本发明提供一种富集BCR-ABL基因片段的方法，所述方法包括以下步骤： 

1）获得受试者的DNA样本库； 

2）获得能够与BCR或ABL基因杂交的DNA探针库； 

3）使所述DNA探针库与所述DNA样本库进行杂交；和 

4）分离步骤3）的杂交产物，然后释放经杂交富集的基因片段。 

其中，所述步骤1）中的DNA样本库由双链DNA片段组成，并且，所述步骤1）包括： 

1-1）提取受试者的全基因组DNA，然后将其片段化；或者 

1-2）提取受试者的mRNA，将其片段化，然后以该经片段化的mRNA为模板合成双链cDNA； 

其中，所述受试者为哺乳动物，优选人，且从受试者的细胞、组织或体液样本中提取全基因组DNA或mRNA； 

优选地，所述DNA片段的长度为150-600bp； 

进一步优选地，所述DNA片段的长度为150-200bp。 

所述步骤2）中的DNA探针库为如上所述的DNA探针库。具体而言，所述DNA探针库包括一个或多个能够与BCR或ABL基因杂交的DNA探 针，所述DNA探针包含以下序列： 

SEQ ID NO.1、SEQ ID NO.2、SEQ ID NO.3、SEQ ID NO.4、SEQ ID NO.5、SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8、SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11、SEQ ID NO.12、SEQ ID NO.13、SEQ ID NO.14、SEQ ID NO.15、SEQ ID NO.16、SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20、SEQ ID NO.21、SEQ ID NO.22、SEQ ID NO.23、SEQ ID NO.24、SEQ ID NO.25、SEQ ID NO.26、SEQ ID NO.27、SEQ ID NO.28、SEQ ID NO.29、SEQ ID NO.30、或SEQ ID NO.31； 

优选地，所述DNA探针的序列如以下序列所示： 

SEQ ID NO.1、SEQ ID NO.2、SEQ ID NO.3、SEQ ID NO.4、SEQ ID NO.5、SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8、SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11、SEQ ID NO.12、SEQ ID NO.13、SEQ ID NO.14、SEQ ID NO.15、SEQ ID NO.16、SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20、SEQ ID NO.21、SEQ ID NO.22、SEQ ID NO.23、SEQ ID NO.24、SEQ ID NO.25、SEQ ID NO.26、SEQ ID NO.27、SEQ ID NO.28、SEQ ID NO.29、SEQ ID NO.30、或SEQ ID NO.31。 

此外，所述步骤3）包括： 

3-1）采用选择性标记标记DNA探针库中的DNA探针；和 

3-2）使所述DNA探针库与DNA样本库进行杂交； 

优选地，所述步骤3-1）中的选择性标记为生物素；进一步优选地，所述步骤3-2）包括在PCR扩增仪中，在65℃下将所述DNA探针库与DNA样本库孵育24小时。 

因此，所述方法的步骤4）中，优选利用DNA探针上的选择性标记分离杂交产物。进一步优选地，所述步骤3-1）中的选择性标记为生物素，所述步骤4）中利用链霉亲和素-生物素的亲和作用分离杂交产物。 

另一方面，本发明还提供一种检测BCR-ABL基因是否存在或检测其基因结构突变的方法，所述方法包括以下步骤： 

1）根据上述方法富集基因片段；和 

2）检测是否富集得到BCR-ABL基因或检测其基因结构突变。 

优选地，所述步骤2）中采用下一代测序技术，通过对富集得到的基因片段进行测序而检测是否富集得到BCR-ABL基因或其基因结构突变。 

又一方面，本发明提供一种用于富集BCR-ABL基因片段的试剂盒，所 述试剂盒包含上述的DNA探针库。具体而言，所述DNA探针库包括一个或多个能够与BCR或ABL基因杂交的DNA探针，所述DNA探针包含以下序列： 

SEQ ID NO.1、SEQ ID NO.2、SEQ ID NO.3、SEQ ID NO.4、SEQ ID NO.5、SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8、SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11、SEQ ID NO.12、SEQ ID NO.13、SEQ ID NO.14、SEQ ID NO.15、SEQ ID NO.16、SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20、SEQ ID NO.21、SEQ ID NO.22、SEQ ID NO.23、SEQ ID NO.24、SEQ ID NO.25、SEQ ID NO.26、SEQ ID NO.27、SEQ ID NO.28、SEQ ID NO.29、SEQ ID NO.30、或SEQ ID NO.31； 

优选地，所述DNA探针的序列如以下序列所示： 

SEQ ID NO.1、SEQ ID NO.2、SEQ ID NO.3、SEQ ID NO.4、SEQ ID NO.5、SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8、SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11、SEQ ID NO.12、SEQ ID NO.13、SEQ ID NO.14、SEQ ID NO.15、SEQ ID NO.16、SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20、SEQ ID NO.21、SEQ ID NO.22、SEQ ID NO.23、SEQ ID NO.24、SEQ ID NO.25、SEQ ID NO.26、SEQ ID NO.27、SEQ ID NO.28、SEQ ID NO.29、SEQ ID NO.30、或SEQ ID NO.31。 

以下是本发明的详细描述： 

本发明提供一种富集BCR-ABL基因片段的方法。具体而言，本发明的方法包括：从哺乳动物例如人的细胞、体液或组织样本中提取基因组DNA或mRNA，经处理或合成cDNA，从而获得片段化的双链DNA作为DNA样本库；此外，针对要富集的BCR-ABL基因，设计与BCR或ABL基因杂交的DNA探针，从中筛选出多个探针作为DNA探针库；然后，将该DNA样本库与DNA探针库进行杂交，从而从DNA样本库中富集得到基因片段。根据本发明的具体实施方式，可以先将DNA探针库中的各个探针进行生物素化，然后在杂交后用链霉亲和素磁珠吸附杂交产物，再从磁珠上释放出富集的基因片段。经适应性处理，可以采用下一代测序基因对基因片段进行检测，以鉴定是否富集得到BCR-ABL基因或其基因结构突变或者进行其基因结构突变的检测，包括单碱基突变、mRNA缺失或增多、mRNA结构颠换、mRNA剪接变化。 

下面以富集得到的基因片段用于基于下一代测序技术的基因结构突变 检测为例，示例性地说明本发明，其中总体工艺流程见图1。 

一、准备mRNA／DNA样本库 

1.准备基因组DNA样本（采用此种方式获得的DNA样本库称为“源自全基因组的DNA样本库”） 

1.1DNA提取 

DNA提取，包括新鲜组织，新鲜血液和细胞，固定和石蜡样本，商业化公司提取试剂盒。以上均按说明书指示方法操作。 

使用分光光度定量仪以及凝胶电泳系统检测DNA模板质量和浓度。dsDNA模板260nm吸光率大于0.05以上，吸光率A260/A280比值在1.8到2之间为合格。 

1.2DNA片段化 

将3微克高质量的基因组DNA用低TE缓冲液稀释至120微升。按照组织匀浆机使用说明书，将DNA片段化，片段长度为150-200碱基。 

DNA过柱纯化，商业化公司纯化试剂盒。 

1.3DNA样本库质量检测 

用生物分析仪进行DNA定性定量分析，确认DNA片段长度峰值合理。 

2.准备cDNA样本 

2.1mRNA提取（采用此种方式获得的DNA样本库称为“源自mRNA的DNA样本库”即cDNA样本库） 

mRNA提取，包括新鲜组织，新鲜血液和细胞，固定和石蜡样本，商业化公司提取试剂盒。以上均按说明书指示方法操作。 

使用分光光度定量仪以及凝胶电泳系统检测mRNA质量和浓度，吸光率A260/A280比值在1.8到2之间为合格。 

2.2mRNA片段化 

采用NEBNext RNA Fragmentation系统或者其他商业化公司mRNA片段化试剂盒。 

mRNA过柱纯化，商业化公司纯化试剂盒 

2.3用商业化公司cDNA合成试剂盒进行mRNA合成第一链以及第二链cDNA。 

cDNA过柱纯化，商业化公司纯化试剂盒。 

3.cDNA／DNA末端修补 

利用T4聚合酶及Klenow大肠杆菌聚合酶片断，对于cDNA/DNA5'突 出粘末端补平以及3'突出粘末端打平，产生平末端，用于后续的平端连接。反应在PCR扩增仪中进行，20摄氏度，30分钟。 

反应材料 体积 纯化后cDNA／DNA样本库 50微升 磷酸化反应缓冲液 10微升 脱氧碱基混合物dNTP（每种10mM） 4微升 T4DNA聚合酶 5微升 Klenow大肠杆菌聚合酶片段 1微升 T4多聚核苷酸激酶 5微升 无核酸酶水 总体积补至100微升

cDNA／DNA过柱纯化，商业化公司纯化试剂盒。 

4.在cDNA／DNA样本3'末端加上碱基A 

反应在PCR扩增仪中进行，37℃，30分钟。 

反应材料 体积 cDNA／DNA样本库 约30微升 10X Klenow大肠杆菌聚合酶缓冲液 5微升 脱氧碱基dATP（1mM） 10微升 Klenow大肠杆菌聚合酶片段 3微升 无核酸酶水 总体积补至50微升

cDNA／DNA过柱纯化，商业化公司纯化试剂盒。 

5.在cDNA／DNA两端加上接头 

反应材料 体积 cDNA／DNA样本库 约15微升 2X T4DNA连接酶缓冲液 5微升 DNA两端接头 6微升 T4DNA连接酶 3微升 无核酸酶水 总体积补至50微升

cDNA／DNA过柱纯化，商业化公司纯化试剂盒。 

如使用mRNA→cDNA，进行6和7； 

如果是用基因组DNA，直接跳到8。 

6.分离出合适长度的cDNA片段 

使用电泳凝胶，对照DNA梯度标准，在凝胶上剪切出150－250碱基cDNA片段。 

将含有cDNA样本库的凝胶样本过柱纯化，商业化公司纯化试剂盒。 

7.cDNA片段样本库质量检测 

使用生物分析仪，进行cDNA定性定量分析，并确认分离出的cDNA片断长度峰值合理。 

8.扩增DNA模板 

聚合酶链反应（PCR），在PCR扩增仪中进行。 

PCR条件：置于PCR扩增仪中，98℃预变性30秒，98℃变性30秒，65℃退火30秒，72℃延伸30秒，共循环15次（cDNA样本库）或者4－6次（DNA样本库）。最后在72℃延伸5分钟。 

PCR扩增产物过柱纯化，商业化公司纯化试剂盒。 

9.扩增后cDNA／DNA样本库质量检测 

使用生物分析仪，进行cDNA／DNA定性定量分析，并确认纯化后片段长度峰值合理，约200bp。 

对于得到的cDNA/DNA样本库，如果cDNA小于30纳克/微升，DNA浓度小于150纳克/微升，须将样品经过真空浓缩机低温干燥（低于45℃），再用无核酸酶水溶解至所需浓度。 

二、准备DNA探针库 

针对BCR和ABL1基因准备DNA探针库。 

探针设计策略见图2。其中，侧翼序列长度为目标区域外的碱基数，大部分时候属于重复或者内含子区域。探针重叠区域为每个探针之间的重叠碱基数目。重叠区域越大，目标区域覆盖率越高；重叠区域越小，甚至探针间隔分布，测序成本则越低。探针长度越长，测序时对于单核苷酸多态性（SNP） 的容忍度越高。 

具体探针设计模式包括：探针长度与探针重叠/间隔区域长度固定，侧翼序列长度变化。这种模式能够保证一致覆盖，通过多个探针对单个基因的富集能够保证富集度高，相对脱靶率低。 

示例性地，最后将探针长度确定为120个碱基，以保证对SNP的容忍度和对基因颠换的敏感性。通过对primer软件的改进，对设计的探针进行分析，以准确的知道探针的退火温度、GC成分连续重复单基数量（如CCCCCCC)。 

首先使用了2倍重叠的方法，对BCR和ABL1基因（基于mRNA）进行全程覆盖。通过对探针质量的分析，选择了退火温度在90-100度之间，GC成分大致控制在40%-60%之间，并且连续单基数量较少的探针。同时通过对人类基因库的检索，保证所选探针的有较小的脱靶率。 

经过筛选，共有123个探针符合标准，通过IDT DNA Technologies，单独合成了每一个探针并用质谱分析保证质量，在5’端有生物素（Biotin）。 

使用每个探针分别对全基因组进行了富集和扩增，最后采用了49个富集效果明显，脱靶率低的探针。 

通过对外显子的分析，在保证基本每个外显子有一个相吻合的探针的情况下，对49个探针进行进一步筛选和优化，最后选定以下31个探针组成用于富集BCR-ABL基因片段的探针库： 

三、DNA捕获探针杂交 

1.将DNA样本库与生物素化的DNA探针库杂交 

将cDNA/DNA样本库与杂交缓冲液混合，反应条件为95℃5分钟，之后保持在65℃。反应在PCR扩增仪中进行。 

然后将该混合物与探针库混合，反应条件为65℃5分钟。将杂交反应置于PCR扩增仪中，65℃孵育24小时。 

四、得到经杂交富集的基因片段 

1.准备链霉亲和素（Streptavidin-Coated）磁珠 

使用Dynabeads链霉亲和素磁珠或者其它商业化公司链霉亲和素磁珠。将磁珠置于混匀仪上混匀，每个样本需要50微升磁珠。 

磁珠洗涤：混合50微升磁珠和200微升结合缓冲液，在混匀仪上混匀，使用Dynal磁选机或者其它商业化公司磁选机，将磁珠与缓冲液分离纯化，缓冲液弃掉不用。重复三次。 

加入200微升结合缓冲液。 

2.分离杂交产物 

混合1中的杂交反应混合物与2中的链霉亲和素磁珠，反复颠倒试管5次。在室温下振摇30分钟。使用Dynal磁选机或者其它商业化公司磁选机，将磁珠分离纯化。 

然后向磁珠中加入500微升洗涤缓冲液，在65℃孵育10分钟，每隔5分钟混匀一次。使用Dynal磁选机或者其它商业化公司磁选机，将磁珠分离纯化。 

以上步骤重复三次。 

3.cDNA/DNA富集样本释放 

将磁珠与50微升洗脱缓冲液混合，室温孵化10分钟，每隔5分钟混匀一次。使用Dynal磁选机或者其它商业化公司磁选机，将磁珠分离弃掉。此时上清液中即含有富集过的基因片段cDNA/DNA样本库。 

将样本库过柱纯化，商业化公司纯化试剂盒。 

经过RT-PCR对富集的BCR-ABL基因片段进行定量，以全基因组为对照，表明上述方法能够将BCR-ABL基因片段富集4598倍。 

五、PCR扩增与纯化 

将富集cDNA/DNA样本库进一步扩增，为测序仪器上样做准备。 

反应材料 体积 富集cDNA/DNA样本库 约30微升 10X高准确率超保真DNA聚合酶缓冲液 5微升 高准确率超保真DNA聚合酶 1微升 正引物 1微升 反引物 1微升 无核酸酶水 总体积补至50微升

PCR条件：置于PCR扩增仪中，98℃预变性30秒，98℃变性30秒，65℃退火30秒，72℃延伸30秒，共循环15次（cDNA样本库）或者4－6次（DNA样本库）。最后在72℃延伸5分钟。 

PCR扩增产物过柱纯化，商业化公司纯化试剂盒。 

六、采用下一代测序技术检测是否存在BCR-ABL基因或其基因结构突变 

使用下一代商业化的测序仪器进行测序，如Roche454、Illumina Hiseq 等。测序结果用已有的测序软件分析包进行分析。 

示例性地，使用TruSeq PE Cluster Kit v3-cBot-HS，使用桥式PCR对DNA样本库模板进行扩增：每个DNA样本片段将会在芯片上形成克隆簇，每条泳道上产生数百万这样的克隆簇。使用Illumina HiSeq2000下一代测序系统，PE-90bp其原理是边合成边测序。和传统Sanger方法相比，利用“可逆性末端终结反应”技术，四种dNTP碱基末端被保护基团封闭，并分别以不同颜色荧光标记。 

经过QC筛选后，对测序结果使用了Bowtie对所得片段进行序列映射，百分之八十以上的测序片断能够被顺利映射。通过统计分析，百分之九十九以上的碱基被测序覆盖；百分之七十五以上碱基被覆盖六十次以上。 

利用Bioconductor软件，成功映射片段进行突变分析。 

综上所述，本发明人开发了基于杂交选择而捕获特定BCR-ABL基因片段序列的方法，采用该方法可以获得成几千倍富集的BCR-ABL基因片段，该经富集的BCR-ABL基因片段样本可以选择性地应用于各种基因检测技术，特别是可以应用下一代测序技术进行基因突变、缺失、增加、和颠换等方面的检测。这种基于杂交选择的靶向捕获，在很多领域极为有用。例如，对于捕获保存久远的DNA，对于临床样本中癌症基因的深度测序，对于罕见基因突变的检测等，由于其特定基因的富集作用，可以获得良好的目标基因检测样本，从而应用检测技术、例如下一代测序技术得以检测。 

并且，对于将通过本发明的方法富集得到的基因片段用于基于下一代测序技术的基因是否存在或基因结构突变检测的应用而言，还具有以下有益效果： 

使用本发明的基因富集方法及筛选得到的特定DNA探针库，能够成数千倍地富集基因片段，从而可以应用下一代测序技术、利用该基因片段的测序，而准确地测定是否存在BCR-ABL基因或其中的各种突变。并且，由于采用下一代测序技术，因此能够一次性检测多种类型基因突变；准确性高，传统技术例如基因芯片技术，通常需要重复两次以上才能确定检测结果，而本发明在一次反应中，对单个碱基进行反复测序，保证了数据的精准度，并且缩短了检测周期；敏感性高，和传统检测技术相比，本发明产生的数据能够达到碱基级的分辨率，使敏感度有了大幅度提高。 

附图说明

以下，结合附图来详细说明本发明的实施方案，其中： 

图1为本发明技术方案的示例性工艺流程图，其中富集得到目标基因，并用于基于下一代测序技术的基因结构突变检测。 

图2为本发明的探针设计策略的示意图。 

图3为实施例1的BCR：ABL融合基因检测结果，其中正对照为BCR-ABL融合基因质粒，负对照为正常人类B-cell全基因组DNA。 

具体实施方式

下面结合具体实施方式对本发明进行进一步的详细描述，给出的实施例仅为了阐明本发明，而不是为了限制本发明的范围。 

用以上描述的实验方法，针对BCR-ABL基因，对全基因组或基于mRNA合成的双链cDNA进行了杂交捕获。所用的细胞系为MDA-MB-231（乳腺癌），HT-29（直肠癌），K-562（血癌），HCT-116（直肠癌），NCI-H1975（非小细胞肺癌）等不同癌症的细胞系。 

下述实施例中的实验方法，如无特殊说明，均为常规方法。下述实施例中所用的试验材料，如无特殊说明，均为常规商店购买得到。 

实施例1：富集并检测K-562细胞系和MDA-MB-231细胞系的BCR-ABL融合基因 

一、分别准备待检测细胞系的DNA样本库 

1.提取K-562细胞系（人类慢性白血病细胞系，来自ATCC细胞库，货号：ATCC-CCL-243）和MDA-MB-231细胞系（人类乳腺癌细胞系，来自ATCC细胞库，货号：HTB-26）的全基因组DNA，然后将其片段化（采用此种方式获得的DNA样本库称为“源自全基因组的DNA样本库”） 

1.1DNA提取 

采用Qiagen Blood&Tissue DNeasy试剂盒（货号：69506），分别从2个细胞系提取全基因组DNA。按说明书指示方法操作。 

使用分光光度定量仪以及凝胶电泳系统检测DNA的质量和浓度。dsDNA的260nm吸光率大于0.05以上，吸光率A260/A280比值在1.8到2之间为合格。 

1.2DNA片段化 

将3微克高质量的基因组DNA用低TE缓冲液稀释至120微升。按照 组织匀浆机使用说明书，将DNA片段化，使片段长度为150-200碱基。 

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将DNA过柱纯化。 

1.3DNA样本库质量检测 

用生物分析仪进行DNA定性定量分析，确认DNA片段长度峰值合理。 

2.分别提取2个细胞系的mRNA，将其片段化，然后合成双链cDNA（采用此种方式获得的DNA样本库称为“源自mRNA的DNA样本库”即cDNA样本库） 

2.1mRNA提取 

采用Qiagen Rneasy试剂盒（货号：74106），从细胞系提取mRNA。按说明书指示方法操作。 

使用分光光度定量仪以及凝胶电泳系统检测mRNA质量和浓度，吸光率A260/A280比值在1.8到2之间为合格。 

2.2mRNA片段化 

采用NEBNext RNA Fragmentation系统或者其他商业化公司mRNA片段化试剂盒，将mRNA片段化。然后使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将mRNA过柱纯化。 

2.3合成双链cDNA 

以该经片段化的mRNA为模板，使用Life Technologies cDNA合成试剂盒（货号：AM1745）合成双链cDNA。 

然后使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将cDNA过柱纯化。 

3.cDNA／DNA末端修补 

利用T4聚合酶及Klenow大肠杆菌聚合酶片断，对于cDNA/DNA5'突出粘末端补平以及3'突出粘末端打平，产生平末端，用于后续的平端连接。反应在PCR扩增仪中进行，20摄氏度，30分钟。 

反应材料 体积 纯化后cDNA／DNA样本库 50微升 磷酸化反应缓冲液 10微升 脱氧碱基混合物dNTP（每种10mM） 4微升 T4DNA聚合酶 5微升 Klenow大肠杆菌聚合酶片段 1微升

[0174] 

T4多聚核苷酸激酶 5微升 无核酸酶水 总体积补至100微升

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将cDNA／DNA过柱纯化。 

4.在cDNA／DNA样本3'末端加上碱基A 

反应在PCR扩增仪中进行，37℃，30分钟。 

反应材料 体积 cDNA／DNA样本库 约30微升 10X Klenow大肠杆菌聚合酶缓冲液 5微升 脱氧碱基dATP（1mM） 10微升 Klenow大肠杆菌聚合酶片段 3微升 无核酸酶水 总体积补至50微升

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将cDNA／DNA过柱纯化。 

5.在cDNA／DNA两端加上接头 

反应材料 体积 cDNA／DNA样本库 约15微升 2X T4DNA连接酶缓冲液 5微升 DNA两端接头 6微升 T4DNA连接酶 3微升 无核酸酶水 总体积补至50微升

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将cDNA／DNA过柱纯化。 

6.从获得的cDNA库中分离出合适长度的cDNA片段 

使用电泳凝胶，对照DNA梯度标准，在凝胶上剪切出150－250碱基cDNA片段。 

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将含有cDNA样本库的凝胶样本过柱纯化。 

7.cDNA片段样本库质量检测 

使用生物分析仪，进行cDNA定性定量分析，并确认分离出的cDNA片断长度峰值合理。 

8.扩增步骤5获得的DNA片段样本库或步骤7获得的cDNA片段样本库 

聚合酶链反应（PCR），在PCR扩增仪中进行。 

PCR条件：置于PCR扩增仪中，98℃预变性30秒，98℃变性30秒，65℃退火30秒，72℃延伸30秒，共循环15次（cDNA样本库）或者4－6次（DNA样本库）。最后在72℃延伸5分钟。 

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将PCR扩增产物过柱纯化。 

9.扩增后cDNA／DNA样本库的质量检测 

使用生物分析仪，进行cDNA／DNA定性定量分析，并确认纯化后片段长度峰值合理，约200bp。因此，分别获得了源自全基因组和mRNA的2个细胞系的DNA样本库。 

对于得到的cDNA/DNA样本库，如果cDNA小于30纳克/微升，DNA浓度小于150纳克/微升，须将样品经过真空浓缩机低温干燥（低于45℃），再用无核酸酶水溶解至所需浓度。本实施例的下文将采用获得的源自全基因组的DNA样本库进行富集和检测。 

二、针对BCR和ABL1基因准备DNA探针库 

探针设计策略见图2。侧翼序列长度为目标区域外的碱基数，大部分时候属于重复或者内含子区域。探针重叠区域为每个探针之间的重叠碱基数目。重叠区域越大，目标区域覆盖率越高；重叠区域越小，甚至探针间隔分布，测序成本则越低。探针长度越长，测序时对于单核苷酸多态性（SNP）的容忍度越高。 

探针设计模式：探针长度与探针重叠/间隔区域长度固定，侧翼序列长度变化。这种模式能够保证一致覆盖，通过多个探针对单个基因的富集能够 保证富集度高，相对脱靶率低。 

探针长度最后确定为120个碱基，以保证对SNP的容忍度和对基因颠换的敏感性。通过对primer5软件的改进，我们能够对设计的探针进行分析，以准确的知道探针的退火温度，GC成分，连续重复单基数量（如CCCCCC)。 

使用了2倍重叠的方法，对BCR和ABL1基因（基于mRNA；BCR NCBI登录号：NM_004327.3，ABL1NCBI登录号：NM_005157.4）进行全程覆盖。通过对探针质量的分析，选择了退火温度在90-100度之间，GC成分大致控制在40%-60%之间，并且连续单基数量较少的探针。同时通过对人类基因库的检索，保证所选探针的有较小的脱靶率。 

经过筛选，共有123个探针符合标准，通过IDT DNA Technologies，单独合成了每一个探针并用质谱分析保证质量，在5’端有生物素（Biotin）。 

使用每个探针分别对全基因组进行了富集和扩增，最后采用了49个富集效果明显，脱靶率低的探针。 

之后，通过对外显子的分析，在保证基本每个外显子有一个相吻合的探针的情况下，对49个探针进行了筛选和优化，最终共筛选得到31个DNA探针，分别为SEQ ID NO.1至SEQ ID NO.31。 

商业合成这些探针。 

三、将DNA样本库与生物素化的DNA探针库杂交 

将DNA样本库与杂交缓冲液（10mM Tris-HCl,2%牛血清白蛋白,pH8.0）混合（混合后，DNA样本库浓度至多不超过50ng/ul），反应条件为95℃5分钟，之后保持在65℃。反应在PCR扩增仪中进行。 

然后以DNA样本库：探针库为1:100的摩尔比，将探针库加入上述混合物，反应条件为65℃5分钟。将杂交反应置于PCR扩增仪中，65℃孵育24小时。 

四、得到经杂交富集的基因片段 

1.准备链霉亲和素磁珠 

使用Dynabeads（Life technologies,货号：11206D）链霉亲和素磁珠或者其它商业化公司链霉亲和素磁珠。将磁珠置于混匀仪上混匀。 

磁珠洗涤：混合50微升磁珠和200微升结合缓冲液（10mM Tris-HCl,2%牛血清白蛋白,pH8.0），在混匀仪上混匀，使用Dynal磁选机或者其它商业化公司磁选机，将磁珠与缓冲液分离纯化，缓冲液弃掉不用。重复三次，每次加入200微升结合缓冲液。 

2.分离杂交产物 

混合步骤三中得到的杂交反应混合物与步骤四的1中得到的链霉亲和素磁珠，反复颠倒试管5次。在室温下振摇30分钟。使用Dynal磁选机或者其它商业化公司磁选机，将磁珠分离纯化。 

然后向磁珠中加入500微升洗涤缓冲液（磷酸缓冲液，0.1%Tween-20,0.1%SDS,pH7.4），在65℃孵育10分钟，每隔5分钟混匀一次。使用Dynal磁选机或者其它商业化公司磁选机，将磁珠分离纯化。以上步骤重复三次。 

3.DNA富集样本释放 

将磁珠与50微升洗脱缓冲液（10mM氢氧化钠溶液）混合，室温孵化10分钟，每隔5分钟混匀一次。使用Dynal磁选机或者其它商业化公司磁选机，将磁珠分离弃掉。此时上清液中即含有富集过的基因片段DNA样本库。 

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将样本库过柱纯化。 

采用RT-PCR荧光定量分析对经探针富集的样本进行基因片段富集效果的检验。其中，以未富集的全基因组DNA样本库作为对照，内参为b-actin，RT-PCR引物序列为： 

SEQ ID NO.32（正向）：TGAGGATACTGTGCTATGAAA 

SEQ ID NO.33（反向）：CCTTGCCCGTCGTGTGCTGAAG 

经RT-PCR检测，本发明的探针能够将BCR-ABL基因片段富集4598倍。因此，能够从全基因组中有选择性的检测BCR-ABL基因及其相关结构突变。 

五、PCR扩增与纯化 

将富集cDNA/DNA样本库进一步扩增，为测序仪器上样做准备。 

反应材料 体积 富集cDNA/DNA样本库 约30微升 10X高准确率超保真DNA聚合酶缓冲液 5微升 高准确率超保真DNA聚合酶 1微升 正引物 1微升 反引物 1微升 无核酸酶水 总体积补至50微升

PCR条件：置于PCR扩增仪中，98℃预变性30秒，98℃变性30秒，65℃退火30秒，72℃延伸30秒，共循环15次（cDNA样本库）或者4－6 次（DNA样本库）。最后在72℃延伸5分钟。 

使用Beckman Coulter Ampure Beads试剂盒（货号：A63880）将PCR扩增产物过柱纯化。 

六、采用下一代测序技术检测是否存在BCR-ABL基因或其基因结构突变 

使用TruSeq PE Cluster Kit v3-cBot-HS，使用桥式PCR对DNA样本库模板进行扩增：每个DNA样本片段将会在芯片上形成克隆簇，每条泳道上产生数百万这样的克隆簇。使用Illumina HiSeq2000下一代测序系统，PE-90bp其原理是边合成边测序。和传统Sanger方法相比，利用“可逆性末端终结反应”技术，四种dNTP碱基末端被保护基团封闭，并分别以不同颜色荧光标记。 

经过QC筛选后，对测序结果使用了Bowtie对所得片段进行序列映射，百分之八十以上的测序片断能够被顺利映射。通过统计分析，百分之九十九以上的碱基被测序覆盖；百分之七十五以上碱基被覆盖六十次以上。 

利用Bioconductor软件，成功映射片段进行突变分析。 

测序结果表明，在K-562细胞系中发现基因融合BCR：ABL，而在MDA-MB-231细胞系中没有检测到突变。通过本发明的方法，K-562的融合点可以准确定位在Chr22:23632743、Chr9:133614146。结果见图3。 

经过DNA提取、PCR扩增和Sanger测序法，以上突变已被验证。 

实施例2：富集并检测KU812细胞系的BCR-ABL融合基因 

利用实施例1中得到的由SEQ ID NO.1至SEQ ID NO.31组成的探针库，采用和实施例1中相同的方法检测KU812细胞系（人类慢性白血病细胞系，来自ATCC细胞库，货号：ATCC-CRL-210）的BCR-ABL基因，在在KU812细胞系中发现BCR:ABL基因融合，融合点为Chr22:23632851、Chr9:133650197。 

经过DNA提取，PCR扩增和Sanger测序法，以上突变已被验证。