一种长片段DNA文库构建方法
本发明涉及生物技术领域，尤其涉及一种长片段DNA文库构建方法。
长片段读取技术(Long Fragment Read,LFR)，是美国Complete Genomics公司提出的DNA文库构建测序技术(Methods and compositions for long fragment read sequencing,United States Patent 8,592,150)。其通过对基因组样品进行384孔板物理分隔的方法，将来自父本与母本的DNA长片段进行分离并分别加入不同的标签序列构建文库，在测序完成后可以完全定相基因组，确认突变位点是否在同一亲本染色体上。文库构建过程中对分离至孔板中的DNA长片段进行MDA扩增，在扩增过程中，掺入dUTP或其它dNTP类似物。后续通过内切酶的作用去除dNTP类似物并使用缺口平移的方法将扩增产物打断成较小的片段，然后使用带方向性的接头通过“连接接头1-延伸反应-连接接头2”的多步反应实现接头A的连接操作。完成以上步骤后进行环化并按照CG建库的后续流程连接测序接头B并完成文库构建。长片段读取技术可以有效改善测序的准确性并且降低起始DNA的投入量。
为满足文库构建实验的需求，分入384孔板中的DNA片段需要先进行扩增处理，CG公司使用多重链置换扩增(Multiple Displace Amplification,MDA)的方法对投入384孔板每个单孔中的DNA长片段进行扩增。而MDA方法容易形成非特异扩增，其在扩增过程中被置换的单链会与新的随机引物结合而形成高级复杂结构影响后续反应。并且，接头A连接的多步酶反应操作复杂繁琐。
发明公开
本发明的一个目的是提供一种长片段DNA文库构建方法。
本发明提供的方法，包括如下步骤:
1)对长片段DNA依次进行转座酶断裂、引入dUTP扩增和去除dUTP，得到断裂片段；
2)将带有不同标签的测序接头单链A和与其部分互补的带有不同标签的测序接头单链B以单链的形式分别加入含有所述断裂片段的体系中反tttttttt应，使所述断裂片段两端连接测序接头，通过测序接头单链A和测序接头单链B中标签序列的排列组合，使每份所述断裂片段对应的测序接头均相互区别，得到连接不同测序接头的产物；
所述带有不同标签的测序接头单链A和所述带有不同标签的测序接头单链B退火可成所述测序接头；
3)以所述连接测序接头后的产物为模板，以与所述测序接头匹配的引物，进行PCR扩增，得到的PCR扩增产物为连接不同测序接头的PCR扩增产物；
4)用所述连接不同测序接头的PCR扩增产物构建文库，即得到长片段DNA文库。
所述构建文库(具体见实施例)为将所述连接测序接头产物的PCR扩增产物依次进行消化dUTP、双链环化、EcoP15酶切、末端补平、去磷酸化、第二接头连接、第二接头连接产物扩增、分离得到单链，单链环化、即得到长片段DNA文库。
上述方法中，所述步骤1)的方法包括如下步骤：
(1)用转座酶将所述长片段DNA进行断裂、并在断裂后的片段两端均连上扩增接头，得到连接扩增接头的产物；
(2)以所述连接扩增接头的产物为模板，以与所述扩增接头匹配的DNA为引物，进行第一次PCR扩增，得到第一次PCR扩增产物；所述第一次PCR扩增过程中引入dUTP；
上述PCR扩增的反应体系(没有模板)：2x buffer 291.2uL，Primer B(20uM)8.48uL，Primer C(20uM)8.48uL，dNTP(各25mM)20.17uL，dUTP(4mM)5.04uL,Pfu turbo Cx聚合酶7.68uL,20％(体积百分比)Triton X-100水溶液20.8uL,10％(体积百分比)Tween20水溶液20.8uL，加入无核酶水至总体积560uL。
(3)去除所述第一次PCR扩增产物中的dUTP，得到缺口，进行缺口平移，并在3’末端加A，得到的产物为1)中所述断裂片段。
上述方法中，所述转座酶为包埋扩增接头的转座酶；
和/或，所述扩增接头为1种或2种，所述扩增接头由转座酶识别单链DNA分子和与其部分反向互补的单链DNA分子形成；tttttttt
和/或，相邻两个所述包埋扩增接头的转座酶的作用位置之间的片段大小为3-10kb。
和/或，(2)中与所述扩增接头匹配的DNA为引物为连接在所述断裂后的片段两端的扩增接头中与所述转座酶识别单链DNA分子反向互补的单链DNA分子中，除所述转座酶识别单链DNA分子反向互补序列外，剩余部分序列形成的引物对。
上述方法中，
所述扩增接头为接头1和接头2，
所述扩增接头为接头1和接头2，所述接头1为由转座酶识别单链DNA分子A和与其部分反向互补的单链DNA分子B组成，所述接头2为由所述转座酶识别单链DNA分子A和与其部分反向互补的单链DNA分子C组成；
所述与所述扩增接头匹配DNA分子作为引物由引物B和引物C组成，所述引物B为所述单链DNA分子B中除去与所述转座酶识别单链DNA分子A互补部分外的其余序列；所述引物C为所述单链DNA分子C中除去与所述转座酶识别单链DNA分子A互补部分外的其余序列；
和/或，去除所述第一次PCR扩增产物中的dUTP的方法包括如下：用尿嘧啶DNA糖基化酶和人源脱嘌呤脱嘧啶核酸内切酶对所述第一次PCR扩增产物进行酶切反应，得到酶切产物；再用聚合酶I、Taq聚合酶和dATP对所述酶切产物进行聚合反应，得到300-1200bp大小的DNA片段。
上述酶切反应体系(不含有模板)：UDG(2U/uL)14.56uL,APE 1(10U/uL)2.9uL,加无核酶水将总体积补充至560uL。
上述聚合反应体系(不含有模板)：聚合酶I(10U/uL)2.86uL,Taq聚合酶(5U/uL)5.7uL,dATP(100mM)50.4uL,加无核酶水将总体积补充至560uL。
上述方法中，所述标签序列是由n个碱基经过排列组合得到，所述碱基为A、G、C、T中至少一种，n大于等于8。
上述方法中，所述带有不同标签的测序接头单链A从5’至3’方向依次包括片段甲、片段乙、标签序列和片段丙；
所述带有不同标签的测序接头单链B从5’至3’方向依次包括与所述片段丙反向互补的片段、标签序列、片段丁和与所述片段甲反向互补的tttttttt片段；
所述步骤3)中，所述引物中的一条引物与所述测序接头单链A中的片段乙匹配(相同或反向互补)，所述引物中的另一条引物与单链B中的片段丁匹配(反向互补或相同)。
上述方法中，所述带有不同标签的测序接头单链A的条数为小于等于72，且各条单链标签序列不同；
所述带有不同标签的测序接头单链B的条数为小于等于72，且各条单链标签序列不同；
所述标签序列中n大于等于8且小于15。
上述方法中，步骤2)中，所述将带有标签的测序接头单链A和与其互补的带有不同标签的测序接头单链B以单链的形式分别加入含有所述断裂片段的体系中反应包括如下步骤：
2)-A，将72条所述带有不同标签序列的测序接头单链A分别加入含有所述断裂片段体系的芯片的纵向72列中；所述测序接头单链A反应缓冲液不含T4连接酶；
2)-B，将72条所述带有不同标签序列的测序接头单链B分别加入2)-A处理后的芯片的横向72行中，反应，得到连接测序接头产物。所述测序接头单链B反应缓冲液含有T4连接酶；
上述方法中，步骤3)中，所述引物中的一条引物与所述测序接头单链A中的片段乙相同或反向互补；
所述引物中的另一条引物与所述测序接头单链B中的片段丁的片段反向互补或相同。
所述PCR扩增是将5184孔板上各孔中所述连接不同测序接头的产物混合一起进行扩增；
上述方法中，所述方法的步骤1)的(2)、(3)和步骤2)均在5184孔芯片中进行。所述在所述断裂片段被分成若干份的情况下，是指将所述断裂片段分装在5184孔芯片的各孔中，分装可在步骤1)或步骤2)中进行。
上述方法中，所述长片段DNA为大于100Kb的片段。具体为大小为400kb的片段；tttttttt
所述转座酶为Tn5转座酶；本实施例中，包埋扩增接头的转座酶为Vazyme公司产品，命名为TruePrep mini DNA Sample Prep Kit(S302-01-B)转座酶试剂盒，最佳使用浓度为100倍稀释转座酶。包埋扩增接头的转座酶和长片段DNA的量见如下反应体系中，上述反应的反应体系包括所述单链DNA分子B、所述单链DNA分子C、反应缓冲液、dNTP、dUTP和聚合酶；2uL 100倍稀释后包埋扩增接头的转座酶、2uL 5x buffer(S302-01,Vazyme公司)、7ng人基因组DNA(长度为大于100kb)，加水补充至10uL。
上述方法中，所述方法的步骤1)的(2)、(3)和步骤2)中均采用wafergen MSND移液平台将各种物质加入芯片的孔中。
上述方法中，在所述用wafergen MSND移液平台将各种物质加入芯片的孔中后，还包括如下步骤：离心沉降各物质。
上述方法中，步骤1)的(1)和(2)中，还包括如下步骤：用变性试剂消化所述转座酶；所述变性试剂具体采用0.1-1％SDS溶液；所述消化的条件具体为25℃静置10分钟。
所述转座酶识别单链DNA分子A的核苷酸序列为序列表中序列5；
所述单链DNA分子B的核苷酸序列为序列表中序列6；
所述单链DNA分子C的核苷酸序列为序列表中序列7；
所述引物B的核苷酸序列为序列表中序列12；
所述引物C的核苷酸序列为序列表中序列13；
所述测序接头单链A的核苷酸序列为序列表中序列1或序列8；
所述测序接头单链B的核苷酸序列为序列表中序列2或序列9；
所述单链引物甲(一条引物)的核苷酸序列为序列表中序列3或序列10；
所述单链引物乙(另一条引物)的核苷酸序列为序列表中序列4或序列11。
上述方法中，步骤1)的(1)和(2)中，还包括如下步骤：用变性试剂消化所述转座酶；所述变性试剂具体采用0.1-1％SDS溶液；所述消化的条件具体为25℃静置10分钟。上述变性试剂为含有SDS溶液，具体由11.2uL 1％SDS(质量/体积百分比)和548.8uL 2x buffer(MP01137,ttttttttComplete Genomics)组成。
上述方法中，步骤1)的(1)中，所述转座酶断裂反应的条件为55℃5分钟；
步骤1)的(2)中，所述扩增的退火温度为68℃，退火时间为18min；
步骤1)的(3)中，所述酶切反应条件为37℃2小时，再65℃15分钟；所述聚合反应条件为23℃1小时，65℃30分钟；
步骤2)的2)-B中，所述反应条件为20℃2小时；
步骤3)中，所述PCR扩增的退火温度为68℃，退火时间为10min。
由上述方法制备得到的长片段DNA文库也是本发明保护的范围。
本发明的另一个目的是提供一种用于长片段DNA文库构建的连接测序接头产物制备方法。
本发明提供的方法，包括如下步骤：为上述方法中的步骤1)-2)。
上述的方法在构建长片段DNA文库中的应用也是本发明保护的范围。
本发明的变性试剂可以为NTbuffer、SDS或其他可以使转座酶变性并从DNA长片段中脱落的试剂。
图1为转座酶打断基因组片段示意图。
图2为wafergen MSND移液平台芯片示意图。
图3为转座酶在与基因组DNA结合示意图。
图4为转座酶脱落后加入与转座酶插入的接头1/2相匹配的引物序列进行扩增。
图5为dUTP在UDG酶及APE1酶的作用下被从PCR产物中切除且在产物上留下1bp的缺口。
图6为依靠dUTP的消化所产生的缺口平移。
图7为第一接头的各自带有标签序列的两条单链的加入方式示意图。
图8为第一接头的两条单链与插入片段连接方式示意图。
图9为第一接头连接完成后PCR扩增反应示意图。
图10为wafergen MSND两种喷液方式示意图。
图11为转座酶片段打断后PCR产物电泳胶图。
图12为芯片建库过程各步骤产物电泳图。tttttttt
图13为EcoP15内切酶切割产物。
图14为单链环化后电泳图。
图15为各孔中的数据量及对应频数的直方图。
图16测序深度分布曲线图。
实施发明的最佳方式
下述实施例中所使用的实验方法如无特殊说明，均为常规方法。
下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
本发明通过利用转座酶插入特异序列，借助该特异序列在带有5184孔的芯片中进行PCR，最终实现文库构建。
下述实施例中借助wafergen公司的MSND移液平台及带有5184孔的芯片，每次可以对样品或反应液进行不同方式的喷液，具体举例如下表1：
表1 为不同方式的喷液
实施例1、构建长片段DNA文库(适配Complete Genomics测序平台)
一、通过转座酶使长片段DNA断开为大小为3-10kb的目的片段
图1为转座酶打断基因组片段示意图，包埋有接头1和接头2的转座酶，在与基因组DNA混合后，随机与基因组中的位置结合，通过控制转座酶的使用量，控制相邻两个转座酶作用位置之间的片段大小约为3～10kb之间。
图2为wafergen MSND移液平台芯片示意图。该芯片横向共72行，纵向共72列，总微孔数量为5184个。
转座酶包埋1种或2种接头都可以。
本实施例采用的包埋扩增接头的转座酶包埋有两种接头，为将转座酶tttttttt识别单链DNA分子A和与A部分反向互补的单链DNA分子B退火得到的接头、转座酶识别单链DNA分子A和与A部分反向互补的单链DNA分子C退火得到的接头。这两种接头按各100uM的浓度等体积混合，然后与转座酶混合并进行包埋得到。本实施例采用的包埋扩增接头的转座酶为Vazyme公司产品Tn5转座酶，命名为TruePrep mini DNA Sample Prep Kit(S302-01-B)转座酶试剂盒。
上述序列信息如下：
转座酶识别单链DNA分子A:5’-CTGTCTCTTATACACATCT-3’(序列5)
单链DNA分子B:5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3’(序列6)
单链DNA分子C:5’-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3’(序列7)
在离心管中，控制包埋有接头序列的Tn5转座酶的使用量对基因组DNA进行作用，使两个相邻的转座酶作用位点在同一条DNA片段相距较长距离(约3-10K)。由于转座酶本身的特性，在作用完成后并不会立刻脱离DNA，而是仍然存在于作用位点，连接住被打断的长片段DNA，故此时整条DNA长片段仍然保持原有长度而非断开状态。
将上述转座酶作用产物稀释至一定浓度(使DNA长片段间可以实现充分分离)，然后利用wafergen的微量移液平台转移到芯片之中，实现DNA长片段的物理分隔。该芯片带有5184个微孔，即最多可将单个DNA样品分离成5184份。后续步骤加入的试剂均使用该移液平台进行。
1、摸索转座酶使用浓度
利用转座酶介导PCR反应，获得实验需求的片段长度3-10kb。
将包埋扩增接头的转座酶分别稀释50倍、100倍和150倍三个梯度。
在PCR反应管中加入如下反应体系：5x buffer 2uL(S302-01,Vazyme公司)，人基因组DNA(从离体血液细胞中提取得到基因组DNA；7ng，大于100kb)6uL，转座酶(不同稀释倍数)2uL。
将上述加有不同稀释倍数转座酶的反应体系混合后分别放入PCR仪中，55摄氏度5分钟反应，然后冷却至常温，反应混合物加入终浓度0.1％tttttttt的变性剂SDS混合，在室温下静置10分钟，得到大小为3-10kb的目的片段。然后将大小为3-10kb的目的片段稀释至终浓度约8pg/uL。
2)PCR反应
将上述1)得到的大小为3-10kb的目的片段加入如下PCR反应体系进行PCR扩增，得到目的片段扩增产物。
用于扩增的引物为连接在所述断裂后的片段两端的扩增接头中与所述转座酶识别单链DNA分子反向互补的单链DNA分子中，除与所述转座酶识别单链DNA分子反向互补序列外，剩余部分序列形成的引物对，具体如下：
引物B：5’-TCG
引物C：5’-GTCTCGTGGGCTCGG-3’(序列13)
上述100uL PCR反应体系包括：大小为3-10kb的目的片段(8pg/uL)5uL、2x buffer 50uL、引物B(20uM)0.5uL、引物C(20uM)0.5uL、DNA聚合酶0.5uL，补水至100uL。
使用表2所示的PCR程序进行扩增。
表2 为PCR程序
不同稀释倍数的转座酶结果如下：
100倍稀释转座酶条件下目的片段扩增产物的结果所示，可以看出，得到3-10kb的扩增产物，为所需目的片段。
50倍稀释转座酶条件下和150倍稀释转座酶条件下均未得到3-10kb的目的片段扩增产物。
可以看出，最优转座酶稀释倍数为100倍稀释。
2、最优转座酶稀释浓度片段化基因组DNA目的片段
采用上述1摸索的最优100倍稀释转座酶对基因组DNA进行处理，使tttttttt基因组DNA约3-10K左右距离间加入转座酶包埋带有的PCR接头序列，并以此接头序列为起点进行PCR扩增，具体如下：
1)、转座酶片段化DNA
将包埋扩增接头的转座酶和人基因组DNA按照如下反应体系进行反应，具体如下：
上述反应体系如下：2uL 100倍稀释后包埋扩增接头的转座酶、2uL 5x buffer(S302-01,Vazyme公司)、7ng人基因组DNA(长度为400kb)，加水补充至10uL。
反应体系可按照比例扩大或缩小。
上述反应为：反应体系放入PCR仪中，55摄氏度5分钟反应，然后冷却至常温，得到反应产物，将其稀释至终浓度约8pg/uL，得到稀释后的反应产物(每3-10KD插入转座酶的基因组DNA)。
使用wafergen MSND移液平台单样品分液流程将稀释后的片段化产物按照35nL/孔分入5184孔板各孔中(孔的溶剂为200nl)；再对芯片进行离心(4000rpm,5min)，使所喷液体沉降至芯片底部，得到含有反应产物的5184孔板。
图3为转座酶在与基因组DNA结合并插入接头1/2后，并不会立刻脱离而是附着在其作用位点中。在反应溶液中加入一定浓度的SDS，使得转座酶从DNA中脱落。
2)、消化转座酶
使用wafergen MSND移液平台单样品分液流程将变性试剂加入上述1)得到的含有反应产物的5184孔板各孔中，再对该芯片进行离心(时间4000rpm,5min)，使所喷液体沉降至芯片底部；离心完毕后室温下静置10分钟左右，用于消化转座酶，使其脱离基因组DNA并实现基因组DNA的片段化，得到3-10KD目的片段，含有其的芯片为含有3-10KD目的片段的5184孔芯片。
本操作可使用其它变性试剂及浓度，不限于SDS。
上述变性试剂由11.2uL 1％SDS(质量/体积百分比)和548.8uL 2x buffer组成；
3)、扩增目的片段得到目的片段扩增产物tttttttt
图4为在完成转座酶的脱落后，在PCR体系中加入引物B和引物C进行扩增。同时，在PCR体系中投入一定比例的dUTP(4％，dUTP:dATP+dTTP+dCTP+dGTP)，使扩增产物中掺入一定比例的dUTP。
使用设计好的引物，对上一步所得的产物进行PCR。在PCR体系之中掺入dUTP(dUTP：dNTP＝4：100)，使最终产物中部分的dTTP被dUTP代替，具体如下：
使用wafergen MSND移液平台单样品分液流程将下述PCR反应体系加入步骤2)得到的含有3-10KD目的片段的5184孔芯片各孔中，之后离心使所喷液体沉降，再将芯片板放入PCR仪中进行PCR扩增反应，得到含有dUTP的目的片段扩增产物的5184孔芯片。
在上述反应体系中加入一定比例的dUTP，使得在扩增中一部分的dTTP被dUTP替代，从而使扩增产物在后期可通过dUTP的去除而被再次打断。
上述PCR扩增反应体系：2x buffer 291.2uL，单链DNA分子B(20uM)8.48uL，单链DNA分子C(20uM)8.48uL，dNTP(各25mM)20.17uL，dUTP(4mM)5.04uL,Pfu turbo Cx聚合酶7.68uL,20％(体积百分比)Triton X-100水溶液20.8uL,10％(体积百分比)Tween20水溶液20.8uL，加入无核酶水至总体积560uL。
上述PCR扩增反应程序如表3所示：
表3 为PCR扩增反应程序
图11为转座酶片段打断后PCR产物电泳胶图。
反应完毕后，将芯片放置于37℃金属浴中静置15.5分钟蒸干多余液体。
二、3-10KD目的片段二次片段化为300-1200bp DNA短片段
1)去除dUTP形成单链缺口tttttttt
使用具有单独去除dUTP活性的酶(UDG、USER等)，对上步PCR产物作用，去除其中的dUTP及该位置的糖骨架，得到中间带有一些1个碱基长度的缺口的DNA双链产物(图5)，具体如下：
使用尿嘧啶DNA糖基化酶(UDG)及人源脱嘌呤脱嘧啶核酸内切酶(APE1)消化扩增中引入的dUTP。
具体如下：
按照wafergen MSND移液平台单样品分液流程将酶切反应体系加入上述一得到的含有dUTP的目的片段扩增产物的5184孔芯片的各孔中，之后离心使所喷液体沉降(4000rpm，5min)；再将芯片放入专用PCR仪中，进行酶切反应，得到含有酶切产物的5184孔芯片。
上述酶切反应体系：UDG(2U/uL)14.56uL,APE 1(10U/uL)2.9uL,加无核酶水将总体积补充至560uL。
上述酶切反应条件为，37℃2小时，65℃15分钟，之后恢复至常温放置。
2)聚合反应完成断裂
缺口平移打断DNA，并在末端加A。在反应体系中加入DNA聚合酶I和Taq酶。在23摄氏度下，利用DNA聚合酶I的作用活性(5’-3’外切酶活性和5’-3’聚合酶活性)，对上一步所得产物进行打断操作。DNA聚合酶I结合于去除dUTP后留下的缺口处，然后从5’-3’方向切除DNA碱基，并从5’-3’方向进行DNA碱基的聚合，实现缺口从5’-3’方向的平移。当DNA双链的两个缺口因为平移而相遇时，该DNA双链便会断开成更小的片段。同时在体系中加入的Taq聚合酶，在65摄氏度下对DNA双链进行补平，并在3’末端加入一个dATP碱基。而在该温度下，DNA聚合酶I因变性而失去活性。
图6为依靠dUTP的消化所产生的缺口，后来投入到反应体系中的DNA聚合酶I结合在该缺口位置发挥其5'-3'方向切割及5'-3'方向合成的作用，使得该缺口沿5'-3'方向平移。双链中位于正反两条链上的缺口同时进行平移并相遇，最终在相遇位置使整个双链断开。而同时加入到体系中的Taq酶，在65℃的温度下发挥其活性，使得产物的3’末端形成突出一个A碱基。tttttttt
利用聚合酶缺口平移的特性，依靠上述去除dUTP后的缺口进行缺口平移实现DNA片段打断，同时在反应体系中加入Taq聚合酶，在反应中使缺口平移后的DNA产物的3'末端加入dATP，具体如下：
按照wafergen MSND移液平台单样品分液流程将聚合反应体系喷液至上述1)含有酶切产物的5184孔芯片各孔中，之后离心使所喷液体沉降(4000rpm，5min)；再将芯片放入专用PCR仪中，进行聚合反应，得到含有300-1200bp大小的DNA片段5184孔芯片。
上述聚合反应体系：聚合酶I(10U/uL)2.86uL,Taq聚合酶(5U/uL)5.7uL,dATP(100mM)50.4uL,加无核酶水将总体积补充至560uL。
上述反应条件为，23℃1小时，65℃30分钟，之后恢复至室温。
三、连接测序接头
在上述二得到的300-1200bp大小的DNA片段两端连接带有部分测序引物的测序接头，测序接头由带有不同标签的测序接头单链A和带有不同标签的测序接头单链B退火形成，为了对芯片上每一个孔都可以在测序的过程中进行区分，设计在芯片中纵向的每一列(对应测序接头单链A)加入编号为1-72的带有不同标签测序接头单链A，在芯片中横向的每一行(对应测序接头单链B)加入编号为1-72的带有不同标签测序接头单链B，这样通过两个单独的标签序列便形成72x72的矩阵，使得每一个孔都获得特有的双标签序列组合。
测序接头单链A
其中，从5’至3’方向所有的粗体如下：
片段甲
片段乙
NNNNNNNNNN为10个碱基组成标签序列；
片段丙
*代表的是硫代磷酸，即在最后一个碱基合成时，使用的是硫代磷酸的核酸。tttttttt
测序接头单链B
/5Phos/
其中，从5’至3’方向所有的粗体如下：
图7为第一测序接头的各自带有标签序列的两条单链的加入方式示意图。为了使芯片中每一个微孔均带有特有的标签序列以进行相互区分，采用两个标签序列组成一个组合的形式。在实验中第一接头的两条单链分别加入。其中第一链的分别72种标签序列的单链(测序接头单链A)，以纵向形式加入，即每一纵列加入的第一链的标签序列相同。第二链的分别72种标签序列的单链(测序接头单链B)，以横向形式加入，即每一横行加入的第二链的标签序列相同。这样最终形成72x72中标签序列的矩阵，每一个孔获得唯一的两两标签序列组合。
为了使用两个单独的标签序列实现72x72种标签序列组合，在连接接头的步骤中，带有标签序列的测序接头单链A以纵向方式加入，带有标签序列的测序接头单链B以横向方式加入，具体如下：
1、测序接头单链A的添加
配制测序接头单链A测序接头单链A反应液如下：10x buffer(25％(体积比)PEG8000、500mM Tris-HCl、10mM ATP、100mM MaCl2)252uL，加无核酶水将总体积补充至957uL；
1)按照wafergen MSND移液平台“72样品35nL分液流程”要求将上述测序接头单链A反应液以每孔11.96uL加入到上述二得到的含有300-1200bp大小的DNA片段5184孔芯片的纵向72列孔中；
2)再将72种不同的测序接头单链A测序接头单链A(10uM)以每孔2.14uL的量按照wafergen MSND移液平台“72样品35nL分液流程”加入上述1)处理后的5184孔板的纵向72列孔中，离心沉降液体，备用，tttttttt得到添加测序接头单链A测序接头单链A的5184孔板。
2、测序接头单链B的添加
配制测序接头单链B测序接头单链B反应液：10xbuffer(25％(体积比)PEG8000、500mM Tris-HCl、10mM ATP、100mM MaCl2)296.4uL、T4连接酶(600U/uL，Enzymatics)125.97uL,加无核酶水将总体积补充至1186uL。
1)按照wafergen MSND移液平台“72样品50nL分液流程”要求将上述测序接头单链B测序接头单链B反应液以每孔13.95uL的体积加入到上述1得到的添加测序接头单链A的5184孔板的横向72行孔中；
2)再将72种不同的测序接头单链B(10uM)以每孔1.65uL的量按照wafergen MSND移液平台“72样品50nL分液流程”加入上述1)处理后的5184孔板的横向72行孔中，离心沉降液体，备用，得到添加测序接头单链B的5184孔板。
将上述添加测序接头单链B的5184孔板置于PCR仪中反应，反应条件为20℃2小时，得到含有连接测序接头产物的5184孔板。
图8为第一接头(测序接头)的两条单链与插入片段连接方式示意图。
图10为wafergen MSND两种喷液方式示意图。左图为“72样品35nL分液流程”配合测序接头单链A的加液方式，使得第一链的分别72种标签序列的单链(测序接头单链A)，实现纵向形式加入。右图为“72样品50nL分液流程”配合第一接头第二链的加液方式，使得第二链的分别72种标签序列的单链(测序接头单链B)，实现横向形式加入。
四、PCR扩增
图9为第一接头(测序接头)连接完成后PCR扩增反应示意图。
将上述三获得的含有连接测序接头产物的5184孔板所有孔的产物均通过离心的方式取出并混合于一离心管之中，然后进行PCR反应扩增，扩增过程会将接头不匹配的部分重新变成匹配的双链。由于产物同一条单链两端是带有不同标签序列的两条接头单链，所以即使PCR后接头的形状进行改变，但不会影响其标签序列组合。也因此，各个孔的产物混合在一起，也不会影响最后对每一个标签序列的组合的拆分，具体如下：
将上述三获得的含有连接测序接头产物的5184孔板所有孔的产物离tttttttt心收集于1.5mL离心管中，并使用1倍AmpureXP磁珠进行纯化，回收溶于100uL TE溶液中纯化产物，得到纯化后连接测序接头产物，对其进行PCR扩增反应，得到连接测序接头产物的PCR扩增产物。
上述PCR扩增采用的引物如下：
PCR引物第一链(上游引物)：GGUCGCCAGCCCUATGGC(序列3)
PCR引物第二链(下游引物)：AGGGCUGGCGACCUTGTCAG(序列4)
上述PCR扩增采用的反应体系如下：连接测序接头产物100uL，2x PfuCx buffer 275uL，PCR引物第一链(20uM)14uL，PCR引物第二链(20uM)14uL，PfuCx聚合酶11uL，加无核酶水补充至550uL。
上述PCR扩增的反应程序如下表4。
表4 为扩增反应程序
将连接测序接头产物的PCR扩增产物使用1倍AmpureXP磁珠进行纯化，纯化产物溶于60uL TE溶液中。
上述各步骤反应产物电泳检测如下图12所示，1、长片段PCR扩增产物(步骤一的2的3)得到的扩增产物)；2、长片段PCR阴性扩增产物(步骤一的2的3)得到的扩增产物，模板用水替代基因组DNA)；3、去除dUTP产物(步骤二的1))；4、去除dUTP反应阴性产物(步骤二的1)，含有dUTP的目的片段扩增产物替换为水)；5、缺口平移末端加A产物(步骤二的2)聚合反应产物)；6、缺口平移末端加A阴性产物(步骤二的2)聚合反应中的酶切产物替换为水)；7、接头连接产物1(步骤三的2的连接测序接头产物)；8接头连接产物2(步骤三的2的连接测序接头产物)；9、PCR产物1(步骤四的PCR扩增产物)；10、PCR产物2(步骤四的PCR扩增产物)；11、PCR阴性产物；M2、100bp分子量标记。从图12可以看tttttttt出各步骤得到的目标产物。
五、消化扩增产物的dUTP
上述连接测序接头产物的PCR扩增产物需要经过消化扩增引物中带有的dUTP，形成粘性末端，从而在后续实验中实现自身环化。dUTP消化反应如下，10x Taq buffer 11uL(RM00059,Complete Genomics),User enzyme(RM00017,Complete Genomics),纯化后PCR产物60uL，加无核酶水补充至110uL，37℃反应1小时，得到反应产物。
六、双链环化
上述五得到的反应产物加入10xTAbuffer(RM06601，Complete Genomics)180uL,加无核酶水至1810uL，混合物平均分成4管，于水浴60℃反应30分钟后，转至常温水浴20分钟。
预先配制反应混合物，无核酶水98uL，20x Circ mix(MP01134，Complete,Genomics)100uL,连接酶(L603-HC-1500,Enzymatics)2uL,混合后平均分成4份加入上述反应产物中混合，室温反应1小时以环化，得到环化产物。
对每一管反应产物，使用AmpureXP磁珠纯化回收，产物回收溶于70uL TE。
七、消化未环化DNA
环化反应中存在相当一部分未环化DNA，为了不影响后续实验结果，需要对未环化的DNA进行消化。消化反应如下，对于每管上述六得到的纯化环化产物，加入9x PS mix(MP01154,Complete Genomics)8.9uL,Plasmid-Safe(RM02046,Complete Genomics)10.4uL,加无核酶水补充至80uL，混合后于37℃反应1小时。反应产物使用AmpureXP磁珠纯化后溶于40uL TE溶液，得到消化后环化的产物。
八、使用EcoP15酶切并筛选目的片段
环化后的产物，最初加入的两端接头连接在一起，此处使用EcoP15酶切，使环化产物由接头两端的EcoP15限制酶切位点处向插入片段内部两端各切下约27bp，后续对酶切下来的产物进行片段筛选，获得中间带接头序列，两端带有酶切下来的插入片段的产物，具体如下：
酶切消化体系配制，上述七得到的消化后环化的产物37uL，5xEcoP 15ttttttttmix 3(MP01149,Complete Genomics)72uL,EcoP 15(RM00063,Complete Genomics)10.8uL,加无核酶水补充至360uL，于37℃反应16小时。
配制PEG32磁珠(按体积比混合，AmpureXP磁珠：0.5％Tween＝100:1，下同)
对于上述360uL酶反应产物，加入252uL PEG32磁珠结合，结合产物保留上清，并加入PEG32磁珠184uL，混合结合后回收磁珠并溶于52uL TE回收液(带有0.001％Tween20，体积比)，得到酶切产物。
图13为EcoP15内切酶切割，片段回收后产物。产物片段在140bp左右。
九、末端补平
对EcoP15酶切产物进行末端补平，以连接后续第二接头。
配制末端补平反应体系，取上述八得到的酶切产物44uL，10x NEB buffer 2(New England biolabs)5.4uL，25mM dNTP 0.8uL，10mg/mL BSA 0.4uL，T4 DNA polymerase(M0203-com,New England Biolabs),12℃反应20分钟，得到末端补平产物。
反应产物使用PEG32磁珠进行纯化并溶于48uL TE。
十、去磷酸化
对上述九得到的末端补平产物进行去磷酸化处理，以配合后续第二接头连接。配制反应体系如下，10x NEB buffer 2(New England Biolabs)5.75uL,Fast AP(EF0651,Fermentas)5.75uL,加入末端补平纯化后产物46uL，于37℃反应45分钟。反应产物使用75uL PEG32磁珠进行纯化，并溶于42uL TE溶液，得到去磷酸化产物。
十一、第二接头连接
第二接头进行定向连接，需要经过4步酶反应，先将上述十得到的去磷酸化产物进行3’末端接头序列的引入，再次进行末端补平并末端引入dATP后，连接引入5’末端接头序列，最后使用一段寡核苷酸序列替代其中一段序列并使用连接酶连接。
接头序列
3’末端接头
/phos/GTCTCCAGTCGAAGCCCGACG/3ddC/tttttttt
/3ddC/AGAGGTCAGCTTCG
5’末端接头
TTGGCCTCCGACT/3dT-Q/
/ddC/TGCTGGCGAACCGGAGGCTGA/5phos/
序列替换的两条序列
寡核苷酸序列1
/52bio/TCCTAAGACCGCTTGGCCTCCGACT
寡核苷酸序列2
/5phos/AGACAAGCTCGAGCTCGAGCGATCGGGCTTCGACTGGAGAC
PCR引物序列
PCR引物1
/52bio/TCCTAAGACCGCTTGGCCTCCGACT
PCR引物2
/5phos/AGACAAGCTCGAGCTCGAGCGATCGGGCTTCGACTGGAGAC
第二接头3’末端序列引入
配制反应混合液，3xHB(MP01139,Complete Genomics)24.7uL,无核酶水1.9uL，T4连接酶(600U/uL，Enzymatics)1.9uL,二接头3’端接头5.6uL(9uM)，加入去磷酸化的产物40uL。混合物于14℃反应2小时，反应完成使用63uL PEG32磁珠进行纯化，并溶于42uL TE溶液。
末端补平加A
配制反应混合液，5x klex NTA mix(MP01150,Complete Genomics)10.7uL,klenow(RM00066,Complete Genomics)1.1uL,无核酶水1.5uL，加入上步骤产物40uL。混合物于37℃反应1小时。反应完成使用69uL PEG32磁珠进行纯化，并溶于42uL TE溶液。
第二接头3’末端序列引入
配制反应混合液，3xHB 24.7uL,无核酶水1.9uL，T4连接酶1.9uL，二接头5’端接头5.6uL(9uM)，加入上步骤产物42uL。混合物于14℃tttttttt反应2小时，反应完成使用63uL PEG32磁珠进行纯化，并溶于42uL TE溶液。
序列替换
配制Oligo反应液，10x Taq buffer(RM00059,Complete Genomics)8uL,100mM ATP(MP01164,Complete Genomics)0.8uL,25mM dNTP 0.32uL,寡核苷酸序列1(20mM)1uL，寡核苷酸序列2(20mM)1uL，加入无核酶水至32uL。
配制酶反应液，10x Taq buffer 0.4uL,T4连接酶4.8uL,Taq聚合酶4.8uL。
实验中将32uL Oligo反应液加入到40uL上述纯化后DNA产物中，然后加入8uL酶反应液，混合后于37℃中孵育20分钟。反应完成使用80uL PEG32磁珠进行纯化，并溶于47uL TE溶液。纯化产物进行定量检测。
得到第二接头连接产物。
十二、第二接头连接产物扩增
往一1.5mL离心管中，加入2xPfuCx mix3 275uL，PfuCx聚合酶11uL，第二接头PCR引物1(20uM)13.75uL，第二接头PCR引物2(20uM)13.75uL，定量后的上述十一制备的第二接头连接产物60ng，加无核酶水补充至550uL。反应物混合后平均分至4个PCR管中进行PCR扩增反应，得到扩增产物。
扩增反应程序如下表5。
表5 为扩增反应程序
反应完成使用550uL PEG32磁珠进行纯化，并溶于85uL TE溶液。对纯化产物进行定量检测，取600ng进行后续单链环化步骤。tttttttt
十三、单链环化
配制1xBWB/tween混合液，往离心管中加入1xBWB(MP01111,Complete Genomics),0.5％Tween20 20uL,充分混合。
清洗磁珠。取stretavidin beads(MP01162，Complete Genomics)120uL，放置于在磁力架上使磁珠充分吸附后，去除上清液。使用1xBBB(MP01110，Complete Genomics)600uL清洗两次，然后加入120uL的1xBBB使磁珠重新悬浮，然后加入1％体积的0.5％Tween20。
使用清洗后的磁珠对产物进行单链分离。取上述十二得到的扩增产物600ng，加入TE补充至60uL，加入4xBBB(MP01145，Complete Genomics)20uL,30uL清洗后的stretavidin beads，充分混合后静置15分钟，然后放置于磁力架上使磁珠充分吸附，使用配制好的1xBWB/tween混合液清洗2次。清洗完成后加入75uL 0.1M NaOH溶液，混合放置2分钟后，放置于磁力架上使磁珠充分吸附，回收上清液。最后往回收的上清液中加入0.3M MOPS acid(MP01165，Complete Genomics)37.5uL，混合均匀。
上述步骤结束后获得112.5uL单链产物，产物进入单链环化步骤。预先配制以下反应液：
引物混合物，取Bridge Oligo(20uM)20uL，加入43uL无核酶水，混合均匀待用；
酶反应混合物，取无核酶水135.3uL，加入10xTA buffer(RM06601，Complete Genomics)35uL,100mM ATP 3.5uL,连接酶(600U/uL)1.2uL,混合均匀待用。
往上述112.5uL产物中，加入引物混合物63uL，酶反应混合物175uL，混合均匀后，置于37℃孵育1.5小时，得到单链环化产物。
图14为单链环化后电泳图。
十四、消化未环化单链产物
往上述十三得到的单链环化产物中，加入无核酶水1.5uL，10xTA buffer 3.7uL，Exo I(M0293L，New England Biolabs)11.1uL，Exo III(M0206L,New England Biolabs)3.7uL，充分混合后，置于37℃反应30分钟。反应完成后加入500mm EDTA 15.4uL,充分混合以终止反应。反应混合物加入500uL PEG32磁珠进行纯化，最终溶解于16uL TE溶液，得到长tttttttt片段文库。
文库制备过程结束。
实施例2、文库构建方法(适配illumina测序平台)
一、通过转座酶使长片段DNA断开为大小为3-10kb的目的片段
与实施例1的步骤相同。
二、目的片段3-10KD目的片段二次片段化为300-1200bp DNA短片段
与实施例1的步骤相同
三、连接测序接头
在上述二得到的300-1200bp大小的DNA片段上一步反应产物的两端连接带有部分测序引物的测序接头，该步骤的接头双链分别带有不同的标签序列，为了对芯片上每一个孔都可以在测序的过程中进行区分，设计在芯片中纵向的每一列(对应测序接头第一条链)加入编号为1-72的标签序列，在芯片中横向的每一行(对应接头第二条链)加入编号为1-72的标签序列。这样通过两个单独的标签序列便形成72x72的矩阵，使得每一个孔都获得特有的双标签序列组合。在上一步反应产物两端连接带有部分测序引物序列的接头。
测序接头单链A
从5’至3’方向依次由Flow cell接头P5(斜体)、10个碱基组成标签序列(N)和与Read 1测序接头互补的单链DNA分子(加粗，*代表的是硫代磷酸，即在最后一个碱基合成时，使用的是硫代磷酸的核酸)组成。
测序接头单链B
/5Phos/
从5’至3’方向依次由Read 2测序接头(加粗)、10个碱基组成标签序列(N)和与Flow cell接头P7互补的单链DNA分子(斜体，Phos即磷酸修饰，AmMo即氨基修饰)组成；72条3’端标签引物的10个碱基组成随机片段均不相同。
图7为第一接头的各自带有标签序列的两条单链的加入方式示意图。tttttttt为了使芯片中每一个微孔均带有特有的标签序列以进行相互区分，采用两个标签序列组成一个组合的形式。在实验中第一接头的两条单链分别加入。其中第一链的分别72种标签序列的单链，以纵向形式加入，即每一纵列加入的第一链的标签序列相同。第二链的分别72种标签序列的单链，以横向形式加入，即每一横行加入的第二链的标签序列相同。这样最终形成72x72中标签序列的矩阵，每一个孔获得唯一的两两标签序列组合。
为了使用两个单独的标签序列实现72x72种标签序列组合，在连接接头的步骤中，第一个带有标签序列的接头单链以纵向方式加入，第二个带有标签序列的接头单链以横向方式加入，具体如下：
1、测序接头单链A的添加
配制测序接头单链A反应液如下：10x buffer(25％(体积比)PEG8000、500mM Tris-HCl、10mM ATP、100mM MaCl2)252uL，加无核酶水将总体积补充至957uL；
1)按照wafergen MSND移液平台“72样品35nL分液流程”要求将上述测序接头单链A反应液以每孔11.96uL加入到上述二得到的含有300-1200bp大小的DNA片段5184孔芯片的纵向72列孔中；
2)再将72种不同的测序接头单链A(10uM)以每孔2.14uL的量按照wafergen MSND移液平台“72样品35nL分液流程”加入上述1)处理后的5184孔板的纵向72列孔中，离心沉降液体，备用，得到添加测序接头单链A的5184孔板。
2、测序接头单链B的添加
配制测序接头单链B反应液：10xbuffer(25％(体积比)PEG8000、500mM Tris-HCl、10mM ATP、100mM MaCl2)296.4uL、T4连接酶(600U/uL，Enzymatics)125.97uL,加无核酶水将总体积补充至1186uL。
1)按照wafergen MSND移液平台“72样品50nL分液流程”要求将上述测序接头单链B反应液以每孔13.95uL的体积加入到上述1得到的添加测序接头单链A的5184孔板的横向72行孔中；
2)再将72种不同的测序接头单链B(10uM)以每孔1.65uL的量按照wafergen MSND移液平台“72样品50nL分液流程”加入上述1)处理后的5184孔板的横向72行孔中，离心沉降液体，备用，得到添加测序接tttttttt头单链B的5184孔板。
将上述添加测序接头单链B的5184孔板置于PCR仪中反应，反应条件为20℃2小时，得到含有连接测序接头产物的5184孔板。
图8为第一接头的两条单链与插入片段连接方式示意图。
图10为wafergen MSND两种喷液方式示意图。左图为“72样品35nL分液流程”配合第一接头第一链的加液方式，使得第一链的分别72种标签序列的单链，实现纵向形式加入。右图为“72样品50nL分液流程”配合第一接头第二链的加液方式，使得第二链的分别72种标签序列的单链，实现横向形式加入。
四、PCR扩增
图9为第一接头连接完成后PCR扩增反应示意图。
将上述三获得的含有连接测序接头产物的5184孔板所有孔的产物均通过离心的方式取出并混合于一离心管之中，然后进行PCR反应扩增，扩增过程会将接头不匹配的部分重新变成匹配的双链。由于产物同一条单链两端是带有不同的两条接头单链/标签序列的，所以即使PCR后接头的形状进行改变，但不会影响其标签序列组合。也因此，各个孔的产物混合在一起，也不会影响最后对每一个标签序列的组合的拆分，具体如下：
将上述三获得的含有连接测序接头产物的5184孔板所有孔的产物离心收集于1.5mL离心管中，并使用1倍AmpureXP磁珠进行纯化，回收溶于100uL TE溶液中纯化产物，得到纯化后连接测序接头产物，对其进行PCR扩增反应，得到连接测序接头产物的PCR扩增产物。
上述PCR扩增采用的引物如下：
PCR引物第一链(上游引物)：AATGATACGGCGACCACCGAGATCT(序列10)
PCR引物第二链(下游引物)：CAAGCAGAAGACGGCATACGAGAT(序列11)
上述PCR扩增采用的反应体系如下：连接测序接头产物100uL，2x PfuCx buffer 275uL，PCR引物第一链(20uM)14uL，PCR引物第二链(20uM)14uL，PfuCx聚合酶11uL，加无核酶水补充至550uL。
上述PCR扩增的反应程序如下表6。
表6 为扩增反应程序tttttttt
将连接测序接头产物的PCR扩增产物使用1倍AmpureXP磁珠进行纯化，纯化产物溶于60uL TE溶液中。
上述各步骤反应产物电泳检测，得到目标产物。
至此文库制备结束。
实施例3、文库测序结果分析
一、测序
使用实施例2制备的DNA长片段文库在illumina平台上进行测序，测序深度40X。
二、比对
使用序列比对软件SOAP aligner 2.20(LiR,LiY,Kristiansen K,et al,SOAP:short oligonucleotide alignment program.Bioinformatics2008,24(5):713-714；LiR,YuC,LiY,et al,SOAP2:an improved ultrafast tool for short read alignment.Bioinformatics 2009,25(15):1966-1967；
http://soap.genomics.org.cn/soapaligner.html)把读取片段(reads)比对人类参考基因组Reference hg19(http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/)，比对存在多重结果时只选取一个比对结果(-r 1)。
三、按照标签组合信息，确定每个孔对应的reads。
四、统计各孔中的数据量及对应频数，绘制直方图(图15)。
数据统计如下：
单位：Mb
均值:36.9000标准差:13.55647tttttttt
最小值:0.3168最大值:123.2000中位数:36.7100
由上可知，本方法可以将所有5184个孔进行成功标记。
实施例4：不同文库构建方法的测序结果对比
一、文库构建和测序
使用CG公司的多重链置换扩增(Multiple Displace Amplification,MDA)的长片段DNA文库构建方法(Methods and compositions for long fragment read sequencing,United States Patent 8,592,150)制备DNA长片段文库，然后在Complete Genomincs平台上进行测序，测序深度80X。
使用实施例1制备DNA长片段文库，然后在Complete Genomincs平台上进行测序，测序深度60X。
二、比对
使用序列比对软件SOAP aligner 2.20(LiR,LiY,Kristiansen K,et al,SOAP:short oligonucleotide alignment program.Bioinformatics2008,24(5):713-714；LiR,YuC,LiY,et al,SOAP2:an improved ultrafast tool for short read alignment.Bioinformatics2009,25(15):1966-1967；
http://soap.genomics.org.cn/soapaligner.html)把测序后得到读取片段(reads)分别比对到人类参考基因组Reference hg19(http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/)，比对存在多重比对时只选取一个比对结果(-r 1)。
三、使用SOAP aligner 2.20中的soap.coverage计算Reads在Reference上的单点覆盖度(-cvg)。
四、绘制测序深度分布曲线图(图16)。
由图16可知，本方法的扩增较MDA扩增方法的均一性更好，更接近于常规测序方法的测序深度分布。
工业应用
本发明在LFR文库构建的原理基础上，取代原DNA扩增使用的MDA扩增方式以及接头A的连接方式，同时使用wafergen微量移液平台的芯片取代384孔板进行试验。本发明使用转座酶插入特定序列，并利用该特定tttttttt序列对每个单孔中的DNA长片段进行PCR扩增，扩增过程中为常规的“变性-退火-延伸”的方式，避免了MDA扩增形成复杂结构以及非特异扩增的问题。本发明中创造性地将两个独立的标签序列设计于双链接头的两个单链之中，然后以单链的形式分别加入，在连接的时候同时将两个标签序列连接到插入片段的两端。由于连接反应在常温下进行，本发明设计的引物在常温反应体系中会先进行退火，然后在连接酶的作用下，与插入片段的两个末端进行连接。在连接反应之前先将打断的插入片段3’末端延伸一个dATP，在连接反应中通过末端A-T配对的方式使接头的两个单链在插入片段的两端进行定向连接。
同时，本发明与wafergen微量移液平台相结合，将DNA长片段分离进入带有5184个细孔的芯片之中，单样品分离孔的数量为384孔板的10倍以上，分离效果更充分，由此获得更好的突变位点定位及组装效果。
本发明的实验证明，本发明通过与微孔芯片相结合的方式，提高了基因组中同源的长片段的分离概率，与Complete Genomics公司的384孔板分离方式相比，芯片中带有的5184孔分离效果更好。384孔板的分离，约有5％的同源长片段会被分入同一孔中而无法区分，而通过5184孔芯片的分离，该概率将进一步下降从而提升分析效率。
本发明使用长片段PCR的方式，对分离后每一个孔中的DNA片段进行扩增，取代了原方法中使用MDA扩增，因此避免了因MDA扩增导致的局部区域形成多级复杂结构的情况。在数据中，这些复杂结构所产生的数据因为无法被识别而去除，限制了所得的有效数据比例。使用长片段PCR扩增的方式，获得的数据比对效率更好。
接头连接过程中采用两接头单链分别带有标签序列并分别加入的方式，在连接过程中边进行退火边实现接头连接。相比较原方法中先进行接头连接-延伸-连接另一端接头的做法，本方法一步反应即可实现两端不同接头的连接并加入两端不同标签序列的组合，大大减少了操作的复杂度及操作时间。tttttttt
