游离DNA测序文库的构建方法
本发明涉及游离DNA测序文库的构建方法。
人体血液中，含有组织细胞释放的游离DNA(cell-free DNA，cfDNA)，其中包含肿瘤组织释放的肿瘤DNA(circulating tumor DNA，ctDNA)，癌症患者血浆游离DNA含有与肿瘤相关的多种突变基因。临床上，癌症患者的ctDNA主要来源于肿瘤细胞分泌、凋亡或坏死，其浓度平均可达180ng/mL，是正常人的20倍以上。研究表明，ctDNA的平均长度为140-170bp，其基因突变情况能直接反映肿瘤组织的基因突变情况[1,2]。随着第二代测序技术的发展，基于基因组变异检测外周血样本的测序技术已经应用于产前诊断等临床层面；而在血液学恶性肿瘤中，高度复发的体细胞结构变异已被证实是有效的分子诊断标记物，包括BCR-ABL原癌基因、免疫球蛋白基因、T细胞受体基因、维甲酸受体α基因等，基于实时定量PCR的分析方法可以定量检测外周血中游离癌细胞的目标基因重排，这些分子诊断标记物可用于实时监控肿瘤发展情况，检测手术后肿瘤残余量，以及长期的临床管理等。Nitzan Rosenfeld团队利用新技术-目标片段深度测序法在37例卵巢癌患者的62份血浆DNA样本中验证了该技术可以有效通过检测ctDNA基因突变监控肿瘤，其结果与digital PCR检测结果高度相符，在发现未知突变方面甚至更优于digital PCR，而该技术较digital PCR技术操作简单。随后，Nitzan Rosenfeld团队从52例转移性乳腺癌患者中选取30例具有PIK3CA或TP53基因点突变或基因结构变异的患者，对141份血浆DNA样本基因突变情况与另外两个检测biomarker，即CA15-3水平和CTC数量进行了比较，检测血浆DNA突变的灵敏度及和病情进展情况的吻合度都较另外两个biomarker高[3,4]。2013年，Murtaza等人跟踪检测了2例乳腺癌、3例卵巢癌和1例肺癌在血浆中游离DNA的变异情况，并阐述了在EGFR、RB1以及MED1等基因在不同时期药物治疗后的突变情况的改变，提示治疗后监控癌症患者外周血血浆中的基因变异能够指导医生及时发现复发与转移的情况并及时调整治疗方案[5]。
近年来肿瘤发病率和死亡率一直呈上升趋势，肿瘤治疗是世界医学重大挑战。许多肿瘤常规化疗效果差，预后不良是困扰临床的重要难题，而肿瘤耐药性则是肿瘤化疗失败的关键因素，不仅产生严重毒副作用，延误治疗，而且给患者带来经济负担。如，EGFR突变阳性肿瘤对厄洛替尼或吉非替尼有较好敏感性而对其他化疗药物敏感性较差，易瑞沙和特罗凯等化疗药物用药前需要对患者EGFR突变情况进行检测。因此，对肿瘤患者进行药物敏感性相关基因检测，针对不同个体选择最优的化疗药物和治疗方案，对提高肿瘤治疗效果、实现个体化医疗有重要意义。而临床上常采用组织穿刺活检等方法获取患者肿瘤组织再进行基因检测，不仅对患者造成严重创伤，而且操作繁琐，也不能达到药物敏感性实时监测的目的。
常规的ILLUMINA Hiseq2000建库需要大量gDNA(1-3微克)，不适用于血浆微tttttttt量片段化gDNA(约几纳克)建库。
参考文献
1.Diehl,F.and M.Li,et al.(2005)."Detection and quantification of mutations in the plasma of patients with colorectal tumors."PNAS 102:16368-16373.
2.Diehl,F.and K.Schmidt,et al.(2008)."Circulating mutant DNA to assess tumor dynamics."Nat Med 14(9):985-90.
3.Forshew,T.and M.Murtaza,et al.(2012)."Noninvasive identification and monitoring of cancer mutations by targeted deep sequencing of plasma DNA."SciTransl Med 4(136):136ra68.
4.Dawson,S.J.and D.W.Tsui,et al.(2013)."Analysis of circulating tumor DNA to monitor metastatic breast cancer."N Engl J Med 368(13):1199-209.
5.Murtaza,M.and S.J.Dawson,et al.(2013)."Non-invasive analysis of acquired resistance to cancer therapy by sequencing of plasma DNA."Nature 497(7447):108-12.
发明公开
本发明提供了游离DNA测序文库的构建方法。
本发明所提供的游离DNA测序文库的构建方法，包括：
1)对游离DNA进行末端修复，得到末端被修复的游离DNA；所述末端修复在含有DNA聚合酶和核苷酸激酶的末端修复体系中进行，所述DNA聚合酶在所述末端修复体系中的浓度可为1.5U-3U/50μL，具体可为1.5U/50μL，所述核苷酸激酶在所述末端修复体系中的浓度可为5U-10U/50μL，具体可为5U/50μL；
2)给所述末端被修复的游离DNA进行末端加“A”，得到末端加“A”的游离DNA；所述末端加“A”在含有Klenow片段(3’-5’exo-)(Klenow(3’-5’exo-))的末端加“A”体系中进行，所述Klenow片段(3’-5’exo-)在所述末端加“A”体系中的浓度可为2.5U-5U/50μL，具体可为5U/50μL；
3)给所述末端加“A”的游离DNA连接文库接头得到带接头的游离DNA；所述连接文库接头在含有文库接头和DNA连接酶的加接头体系中进行，所述文库接头在所述加接头体系中的浓度为20-40pmol/50μL，具体可为20pmol/50μL；所述DNA连接酶在所述加接头体系中的浓度可为1200U-2400U/50μL，具体可为1200U/50μL；
4)扩增所述带接头的游离DNA，得到游离DNA文库；所述扩增的循环数可为8-12，具体可为10。
上述游离DNA测序文库的构建方法中，所述末端修复体系的体积可为20-100μL，具体可为50μL。
上述游离DNA测序文库的构建方法中，所述DNA聚合酶为T4 DNA聚合酶,所述核苷酸激酶为T4多聚核苷酸激酶，所述DNA连接酶为T4 DNA连接酶。tttttttt
上述游离DNA测序文库的构建方法中，所述末端修复体系中含有Klenow片段(Klenow Fragment)，所述Klenow片段在所述末端修复体系的浓度可为2.5-5U/50μL，具体可为2.5U/50μL。
上述游离DNA测序文库的构建方法中，所述末端修复可在20℃进行30min。
上述游离DNA测序文库的构建方法中，所述末端加“A”可在37℃进行30min。
上述游离DNA测序文库的构建方法中，所述连接文库接头可在20℃进行15min。
上述游离DNA测序文库的构建方法中，4)所述扩增在含有Platinum Pfx DNA聚合酶(Platinum Pfx DNA polymerase)的扩增体系中进行，所述Platinum Pfx DNA聚合酶在所述扩增体系中的浓度为2U-4U/50μL，具体可为2U/50μL。
上述游离DNA测序文库的构建方法中，所述扩增的循环温度程序为先94℃20s，然后62℃40s，最后72℃50s。
上述游离DNA测序文库的构建方法中，所述构建方法可包括从体液中提取所述游离DNA的步骤。
上述游离DNA测序文库的构建方法在游离DNA的序列测定中的应用也属于本发明的保护范围。
上述游离DNA测序文库的构建方法在检测游离DNA变异情况中的应用也属于本发明的保护范围。
上文中，所述游离DNA可为体液中的游离DNA(离体的体液游离DNA)，如血液游离DNA(离体的血液游离DNA)、血浆游离DNA(离体的血浆游离DNA)。
实验证明，采用本发明的游离DNA测序文库的构建方法进行文库构建的测序分析方法，对外周血血浆游离DNA的目标区域捕获覆盖度平均达到98.7％，目标区域平均捕获效率达到了61.2％，与同类研究平均10-15％左右的捕获效率相比，该测序分析方法的捕获效率甚至达到细胞系/组织样本等有完整基因组DNA且质量较好样本的目标区域捕获建库的捕获效率。本发明的测序分析方法可以检出肺癌患者EGFR的点突变、短插入/缺失，MET和EGFR的拷贝数突变。本发明的测序分析方法可以检出肺癌患者外周血血浆中频率极低的突变的频率比例和药物耐药及敏感程度的突变频率比例。
图1为肺癌患者外周血血浆游离DNA文库2100检测峰图。
图2为肺癌患者外周血血浆游离DNA目标区域捕获覆盖度分布。
图3为肺癌患者外周血血浆游离DNA目标区域捕获效率分布。
图4为临床用药MET基因拷贝数在同一肺癌患者不同病程进展中MET基因拷贝数的改变情况：1代表病程进展Ⅰ期，2代表病程进展Ⅱ期，3代表病程进展Ⅲ期。
实施发明的最佳方式
以下的实施例便于更好地理解本发明，但并不限定本发明。下述实施例中的实验方法，如无特殊说明，均为常规方法。下述实施例中所用的试验材料，tttttttt如无特殊说明，均为自常规生化试剂商店购买得到的。以下实施例中的定量试验，均设置三次重复实验，结果取平均值。
下面提供的游离DNA测序文库的构建方法是对ILLUMINA Hiseq2000的DNA测序文库的构建方法中的如下步骤进行改进得到的适合构建游离DNA测序文库的方法：
1)对DNA进行末端修复的操作中，T4 DNA聚合酶在反应体系中的浓度分别为1.5U-3U/50μL，具体可为1.5U/50μL；T4多聚核苷酸激酶在反应体系中的浓度可为5U-10U/50μL，具体可为5U/50μL；
2)DNA进行末端加“A”的操作中，Klenow片段(3’-5’exo-)在反应体系中的浓度可为2.5U-5U/50μL，具体可为5U/50μL；
3)对DNA进行连接文库接头的操作中，文库接头在加接头体系中的浓度为20-40pmol/50μL，具体可为20pmol/50μL；T4 DNA连接酶在加接头体系中的浓度可为1200U-2400U/50μL，具体可为1200U/50μL；
4)扩增带接头DNA的操作中，扩增的循环数增至8-12个循环，具体可为10。
下面以肺癌患者为例阐明本申请的游离DNA测序文库的构建方法。
10×Polynucleotide Kinase(PNK)Buffer、dNTP Solution Set、T4 DNA Polymerase、Klenow Fragment、T4Polynucleotide Kinase、Klenow(3’-5’exo-)试剂盒、dATP、2×Rapid l igation buffer、T4 DNA Ligase(Rapid)均为Enzymatics公司的产品；PE Adapter oligo mix、Invitrogen Dynabeads M-280、
实施例1、肺癌患者外周血血浆游离DNA测序文库的构建方法及其应用
本实施例提供了外周血血浆游离DNA的测序分析方法，包括构建肺癌患者外周血血浆游离DNA测序文库和肺癌患者外周血血浆游离DNA测序文库的测序及分析这两个步骤。
文库构建方法为在Hiseq文库构建方法基础上改进得到的适合肺癌患者外周血血浆游离DNA测序文库的构建方法。
一、构建肺癌患者外周血血浆游离DNA测序文库
构建肺癌患者外周血血浆游离DNA测序文库包括下述步骤(一)、(二)、(三)、(四)和(五)。
(一)肺癌患者外周血血浆游离DNA的提取
采用SnoMagTM Circulating DNA试剂盒进行肺癌患者外周血血浆游离DNA的提取，具体如下:
1、取200μL肺癌患者外周血血浆于1.5mLEP管，加入600μL buffer LSBtttttttt(LSB缓冲液)。
2、加入20μLNanoMag Circulating Beads混匀，室温放置10min，每2-3min混匀一次。
3、将EP管置于磁力架上吸附1min，弃上清。
4、取下EP管加入150μL Buffer WA(WA缓冲液)，混匀。
5、将EP管置于磁力架上吸附1min，弃上清。
6、取下EP管加入150μL 75％乙醇，混匀。
7、将EP管置于磁力架上吸附1min，弃上清。
8、重复步骤6-7一次。
9、室温干燥EP管中的磁珠5min。
10、加入32μL elution buffer(洗脱液)混匀磁珠，室温静置5min。
11、将EP管置于磁力架上吸附1min，得到的上清液即为肺癌患者外周血血浆游离DNA溶液，转移该上清液至新的1.5mL EP管。测定16例肺癌患者外周血血浆游离DNA溶液中肺癌患者外周血血浆游离DNA的浓度。结果如表1所示：
表1、肺癌患者外周血血浆游离DNA的浓度
(二)对肺癌患者外周血血浆游离DNA进行末端修复，得到末端被修复的游离DNA
末端修复方法按照表2配制的50μL末端修复体系，该50μL末端修复体系中，T4 DNA聚合酶的浓度为1.5U/50μL，T4多聚核苷酸激酶的浓度为5U/50μL，Klenow片段的浓度为2.5U/50μL。
表2、末端修复方法的50μL末端修复体系
将50μL末端修复体系在Thermomixer恒温干浴器中20℃温育30min，得到末端被修复的游离DNA溶液。采用90μL AmpureXP beads纯化磁珠纯化获得的末端被修复的游离DNA，并用36μL Elution buffer(洗脱缓冲液)从AmpureXP beads洗脱该末端被修复的游离DNA得到纯化后的末端被修复的游离DNA溶液。
(三)对末端被修复的游离DNA进行末端加“A”，得到末端加A的游离DNA
末端加“A”方法按照表3配制的50μL末端加“A”体系，该50μL末端加“A”体系中，Klenow片段(3’-5’exo-)的浓度为5U/50μL。
表3、末端加“A”方法的50μL末端加“A”体系
将末端加“A”方法的50μL末端加“A”体系在Thermomixer恒温干浴器中37℃温育30min，得到末端加“A”的游离DNA溶液。采用90μL AmpureXP beads纯化获得的末端加“A”的游离DNA，并用24μL Elution buffer从AmpureXP beads洗脱该末端加“A”的游离DNA得到末端加“A”的游离DNA溶液。
(四)对末端加“A”的游离DNA的两端连接Illumina Hiseq2000文库接头，得到带接头的游离DNA
在反应体系中，文库接头(PE Adapter oligo)的浓度为20pmol/50μL，T4 DNA连接酶的浓度为1200U/50μL。
表4、两端连接文库接头方法的50μL加接头体系
将两端连接文库接头方法的50μL加接头体系在Thermomixer恒温干浴器中tttttttt20℃温育15min，得到带接头的游离DNA溶液。采用75μL AmpureXP beads纯化获得的带接头的游离DNA，并用34μL Elution buffer从AmpureXP beads洗脱该带接头的游离DNA得到纯化后的带接头的游离DNA溶液。
(五)对带接头的游离DNA进行文库PCR预扩增反应，构建肺癌患者外周血血浆游离DNA文库
文库PCR预扩增方法按照表5配制50μL文库PCR预扩增反应体系，该50μL文库PCR预扩增反应体系中Platinum Pfx DNA聚合酶(Platinum Pfx DNA polymerase)的浓度为2U/50μL。
表5、文库PCR预扩增方法的50μL文库PCR预扩增反应体系
将文库PCR预扩增方法的50μL文库PCR预扩增反应体系预先在PCR仪上于94℃预加热2min，使模板DNA充分变性然后进入扩增循环阶段；扩增的循环温度程序为94℃保持20s使模板变性，然后62℃保持40s使得引物与模板退火，最后72℃保持50s使引物在模板上延伸，合成DNA，完成一个循环，共进行10次这样的循环，使扩增的DNA片段大量积累；随后72℃保持10min，保证产物延伸完整，4℃保存。将经过PCR预扩增制备的肺癌患者外周血血浆游离DNA文库采用60μL AmpureXP beads纯化获得肺癌患者外周血血浆游离DNA文库，并用22μL Elution buffer(洗脱缓冲液)从AmpureXP beads(纯化磁珠)洗脱扩增的肺癌患者外周血血浆游离DNA，获得肺癌患者外周血血浆游离DNA文库(简称Pre-PCR library)。对肺癌患者外周血血浆游离DNA文库使用QubitdsDNA HS assay kit(Qubit双链DNA高敏检测试剂盒)定量确定该文库中DNA的浓度。
结果表明，肺癌患者外周血血浆游离DNA的Pre-PCR 16个文库的浓度如表6所示：
表6：血浆游离DNA的Pre-PCR文库浓度
二、肺癌患者外周血血浆游离DNA测序文库的测序及分析
(一)肺癌患者外周血血浆游离DNA文库与Nimblegen芯片杂交，进行目标区域DNA的捕获
1、将步骤一中(五)的Pre-PCR l ibrary，按照16例肺癌患者的Pre-PCR library混合在一起的方式进行混合得到一个混合池进行文库杂交，最终得到1个上机文库。一个混合池中DNA的含量是1000ng，每例肺癌患者的Pre-PCR library含量是62.5ng。
2、按照表7所示制备杂交文库体系，将该杂交文库体系使用真空浓缩仪于60℃浓缩15min，获得干燥的杂交文库产物。
表7、肺癌患者外周血血浆游离DNA与Nimblegen芯片杂交文库体系
3、按照表8配置杂交缓冲液反应体系。
表8、杂交缓冲液反应体系
4、向步骤2获得干燥的杂交文库产物中加入10.5μL配制好的步骤3的杂交缓冲液反应体系，于95℃处理10min，获得溶解于杂交缓冲液反应体系中的杂交文库体系溶液。
5、文库与探针杂交反应：取4.5μL Nimblegen芯片加入步骤4杂交文库体系溶液，于47℃杂交72h，获得杂交产物溶液。
(二)杂交后进行肺癌患者外周血血浆游离DNA目标区域DNA洗脱
1、将Invitrogen Dynabeads M-280磁珠剧烈震荡混匀重悬，备用。
2、添加100μL Invitrogen Dynabeads M-280磁珠至1.5mL EP管中，洗涤磁珠，步骤如下：
1)加入200μL Binding Buffer(Binding缓冲液)，剧烈震荡5s使得Invitrogen Dynabeads M-280重悬，置于磁力架上静置2min，待液体完全澄清后小心吸弃上清，获得Invitrogen Dynabeads M-280沉淀；
2)重复步骤1)1次，共洗涤2次；tttttttt
3)向装有Invitrogen Dynabeads M-280沉淀的1.5mL的EP管中加入100μL Beads Wash Buffer(磁珠洗涤缓冲液)重悬该磁珠。
3、把重悬后的磁珠转到含有步骤(一)的杂交产物溶液的样本管中，置于PCR仪中47℃孵育45min，每隔15min混匀一次，共三次。
4、将孵育结束后的磁珠悬浊液转移到新的1.5mLEP管中，置于磁力架上静置后弃上清。
5、向步骤4的磁珠沉淀中加入Wash Buffer I(洗涤缓冲液I)(47℃)100μL，混匀样本，置于磁力架上静置3-5min，待液体完全澄清，小心吸弃上清。
6、向步骤5的磁珠沉淀中加入Stringent buffer(Stringent缓冲液)(47℃)200μL，混匀，置于Thermomixer中47℃孵育5min，瞬时离心，置于磁力架上置于3-5min，待液体完全澄清，小心吸弃上清。
7、重复步骤6操作1次，共2次。
8、向步骤7的磁珠沉淀中加入Wash Buffer I(RT)200μL，置于Thermomixer中于20℃1400rpm离心2min，之后置于磁力架上静置3-5min，待液体完全澄清，小心吸弃上清。
9、重复步骤8操作2次，共3次。
10、向步骤9的磁珠沉淀中加入140μL Nuclease Free-water(无核酸水)重悬磁珠，得到含有捕获目标区域DNA的磁珠悬浊液。
(三)对捕获的目标区域DNA进行桥式PCR扩增
1、按照表9在1.5mL EP管中配置100μL桥式PCR扩增反应体系。
表9、桥式PCR扩增反应体系
将100μL桥式PCR扩增反应体系预先在PCR仪上于94℃预加热2min，使模板DNA充分变性然后进入扩增循环阶段；扩增的循环温度程序为94℃保持15s使模板变性，然后58℃保持30s使得引物与模板退火，最后72℃保持30s使引物在模板上延伸，合成DNA，完成一个循环，共进行15次这样的循环，使扩增的DNA片段大量积累；随后72℃保持5min，保证产物延伸完整，4℃保存，得到桥式PCR扩增产物。tttttttt
2、桥式PCR扩增产物纯化：将装有步骤1的桥式PCR扩增产物的PCR管置于磁力架上，静置3-5min至澄清，小心吸取上清转入新的1.5mLEP管中，加入120μL AmpureXP beads纯化上清液，最后用31μL Elution buffer(洗脱液)洗脱产物，得到最终上机文库。
3、上机文库质检：应用Agilent 2100Bioanalyzer和荧光定量PCR(QPCR)进行文库质量检测，2100检测峰图如图1所示，检测结果如表10所示：
表10：文库2100和QPCR检测结果
QPCR结果大于5nM,小于500nM，质量浓度大于0.5ng/μL,片段范围250-330bp即满足上机要求。该结果表明，文库的大部分片段集中在280-300bp，且浓度均符合上机测序要求。
(四)对捕获的目标区域DNA进行上机测序及数据分析
1、Illumina Hiseq2000上机测序：以QPCR检测浓度为准，按照Illumina Hiseq2000上机测序程序对步骤(三)中的纯化的桥式PCR产物进行上机测序。
2、下机数据分析：对上机测序数据进行数据过滤、数据比对以及变异检测分析。
2.1目标区域捕获覆盖度及捕获效率
测序分析方法的目标区域捕获覆盖度分布如图2所示，从图2中可以看出，目标区域捕获覆盖度平均达到98.7％。说明该方法得到的数据几乎覆盖到了所有目标区域，能够完全满足后续突变检测的需要。该测序分析方法的目标区域捕获效率分布如图3所示，平均捕获效率达到了61.2％，与同类研究平均10-15％左右的捕获效率相比，该测序分析方法的捕获效率甚至达到细胞系/组织样本等有完整基因组DNA且质量较好样本的目标区域捕获建库的捕获效率。
2.2突变检测结果
临床上，肺癌患者EGFR的点突变、短插入/缺失，MET和EGFR的拷贝数检测已经成为使用靶向药物易瑞沙和特罗凯治疗前必须进行的检查。采用该测序分析方法进行测序分析的肺癌患者外周血血浆中也检测到了这些基因的上述突变情况，同时将该结果与患者临床信息进行的比较，结果发现该测序分析方法的突变检测结果与临床结果一致。
临床上EGFR G719X是患者使用靶向药物敏感位点，该位点在患者群中发生突变的比例在3％左右。EGFR T790X突变是患者在使用靶向药物过程中可能产生的耐药继发突变，是患者产生抗药性的原因之一。该测序分析方法得到的突变检测结果表明，在由223例肺癌患者外周血血浆样本组成的群体中，EGFR T790X的突变频率是1.5％，EGFR G719X的突变频率是0.5％。说明该测序分析方法可以检出外周血血浆中频率极低的和药物耐药及敏感程度的突变频率比例。tttttttt
临床上随着靶向药物治疗的进行，MET等基因拷贝数的增加会引起患者对靶向药物的耐药，进而出现病程进展。采用该测序分析方法检测了一例肺癌患者在不同临床分期病程进展I期、病程进展II期和病程进展III期的外周血血浆游离DNA中的MET基因的拷贝数变化情况，结果表明从病程进展I期至病程进展III期，该例肺癌患者的外周血血浆游离DNA中的MET基因的拷贝数不断增加(图4)。临床诊断结果表明，该例肺癌患者从病程进展I期至病程进展III期出现了对靶向药物耐药的情况，可以从另一方面解释了该例肺癌患者的病程进展III期的血浆样本中检测cfDNA的MET的拷贝数异常增加可能是患者出现临床耐药症状的基因组证据之一。说明该测序分析方法可以很好地针对拷贝数的变化实时跟踪患者的病程进展情况。
肺癌患者ALK基因与其他基因的发生融合，也是肺癌患者的一种显著分子诊断标记物。采用该测序分析方法在一轮15个肺癌患者外周血血浆样本的测试中发现，ALK基因断裂点(breakpoint)的检出频率为93.3％(检测总样本数为15个，检出ALK基因断裂点样本数为14个)。
工业应用
鉴于ctDNA的基因突变情况能直接反映肿瘤组织的基因突变情况，可以采用本发明的游离DNA测序文库的构建方法进行游离DNA测序文库的构建，进而通过测序分析获得肿瘤患者的ctDNA突变情况及其含量的变化，以进行肿瘤的早期诊断检测、遗传杂合性评估、肿瘤动态变化监测、靶向治疗基因突变分析、肿瘤早期治疗反应评价、肿瘤微量残留检测和肿瘤耐药性实时检测等。tttttttt
