一种基因组DNA的PCR-free测序文库制备方法
技术领域tt
本发明涉及分子生物学中高通量测序技术的领域。更具体涉及一种基因组DNAttPCR-free测序文库制备方法，它适用于基因组DNA不经过PCR扩增就进行文库构建的样品。tt
背景技术tt
二代测序技术是现代分子生物学高通量测序研究中最常用的技术。由于传统的一tt代测序已经不能完全满足研究者的需要，对模式生物进行基因组重测序以及非模式生物基tt因组测序都需要费用更低、通量更高、速度更快的测序技术，因而催生了第二代测序技术tt(Next-generationsequencing)。现有第二代测序技术平台主要包括Roche/454FLX、ttIllumina/SolexaGenomeAnalyzer和AppliedBiosystemsSOLIDsystem三种，其核心原tt理是边合成边测序(SequencingbySynthesis)，即通过捕捉新合成的末端标记碱基信息tt来确定DNA的序列。这三种测序平台各有优点，其中454FLX的测序读长(Read)最长，Solexatt测序性价比最高，SOLID测序的准确度最高。这三种测序平台中Illumina公司的测序仪器占tt了全球市场七成的份额，其Solexa的升级版Hiseq系列是目前使用最多的测序仪。Hiseq测tt序平台单次反应的数据量由最初几十G已提升到1.5～1.8Tb，极大的降低了测序成本，使得tt一千美金测一个人类基因组已经成为了可能。目前第三代测序技术(单分子测序技术)仍旧tt存在准确度和通量不高等缺点，部分三代测序仪还处于研发和测试阶段，而二代测序技术tt自身也会在不断完善和发展，因此，在今后相当长的一段时间内二代测序还是占据主流。第tt二代测序技术对DNA测序的主要原理是先将的DNA进行片段化，片段化DNA末端修复，然后在tt两侧加上特异的接头，随后用不同的方法产生几百万个空间固定的PCR克隆阵列，利用引物tt杂交和酶延伸反应对每个延伸反应所掺入的荧光标记进行成像检测就可获取测序数据。tt
二代测序对DNA进行测序主要包括到文库制备和上机测序两个过程。在文库制备tt过程中，一般通过标准的PCR来扩增随机打断的基因组片段。但对于一些特殊模板，存在二tt级结构复杂或者热稳定性不好因素，从而影响模板PCR的扩增偏好，所以并非所有基因组序tt列都能在PCR扩增文库中同等体现。特别是对一些高GC或者高AT含量的模板，有时很难用ttPCR方法来扩增进行文库构建。PCR-free文库和常规PCR文库上机测序时并没区别，只是建tt库过程中不做PCR，PCR-free文库理论上可以改善数据读取分布并产生更均一的基因组覆tt盖。采用Illimina公司测序仪对文库上机测序时，PCR-free文库和常规PCR文库的簇扩增相tt似，不过也存在一些差异，其差异主要在于构建接头的方式。在Illimina测序仪的Flowttcell上簇扩增的接头包含一些额外的序列，能帮助与连接在流动槽表面的单链DNA进行分tt子杂交，流动槽本身会选择完全相连的模板分子且只扩增与接头序列完全连接的模板链，tt簇扩增便开展了类似普通PCR文库构建时的PCR富集扩增，这样可降低文库构建时PCR扩增tt引入的错误碱基和由此造成的偏向性。Illumina的TruSeqDNAPCR-free样品制备试剂盒tt采用了一种基于磁珠的流程，理论上可在一天内完成，并可与Illumina测序平台不断增加tt的读长兼容。但此试剂盒价格昂贵，单次建库成本高。另外，由于全程采用的是磁珠纯化的tt流程，其缺点是插入片段分布较宽，尤其在小片段区域分布较多，容易导致序列在两端测序tt时被测穿而产生数据浪费。综上，PCR-free文库是一种建库方法的补充，由于不进行PCR扩tt增而直接制备成上机文库，可提高一些高GC或者高AT区域的覆盖度，从而减少了PCR引入的tt错误碱基、数据偏向以及序列重复。tt
发明内容tt
本发明的目的在于针对现有技术的不足，提供一种新的快速高效的可应用于ttIllumina二代测序平台的基因组DNA的PCR-free测序文库制备方法。此方法不进行PCR扩增tt而直接制备成上机文库，可提高一些高GC或者高AT区域的覆盖度，从而减少了PCR引入的错tt误碱基、数据偏向以及序列重复。因此，此方法特别适用于高GC或者高AT含量的基因组DNAtt或者为了降低非模式生物DNADenovo测序中碱基错误率和序列重复的文库构建。tt
本发明提供以下技术方案：tt
本发明提供的基因组DNA的PCR-free测序文库制备方法具体流程为：具体包括以tt下步骤：tt
(1)对基因组DNA进行超声打断；tt
(2)对超声打断的片段进行纯化回收；tt
(3)利用NEB末端修复试剂对纯化产物进行末端修复；tt
(4)对末端修复产物进行纯化回收；tt
(5)利用NEB末端加A试剂对上述纯化产物进行DNA末端加A反应；tt
(6)对上述加A产物进行纯化回收；tt
(7)对上述纯化产物进行连接反应，本发明专利设计了多种兼容于Illumina二代tt测序仪的测序接头，接头序列含有标签(Index)序列，可同时对多个样品进行标记，接头引tt物合成采用HPLC，使用前需进行一定条件的处理；tt
接头序列为：tt
P5(5’-3’):AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATC*Ttt(SEQIDNo.1)tt
P7(5’-3’):/5’Pho/GATCGGAAGAGCACACGTCTGAACTCCAGTCACNNNNNNATCTCGTATGCCGttTCTTCTGCTTG(SEQIDNo.2)其中P7引物序列种的六碱基序列NNNNNN代表标签(Index)序tt列，标签(Index)序列用于标记不同的样品，可以将带不同标签(Index)序列的样品构建混tt合文库；tt
(8)对连接产物进行片段筛选，回收产物即为测序文库；tt
(9)对文库进行质检和上机测序，检测测序文库的浓度和片段大小，利用Illuminatt公司的二代测序仪对浓度和片段大小符合要求的文库进行高通量测序；tt
本发明利用NewEnglandBioLabs(NEB)公司的常规建库试剂，采用磁珠和切胶回tt收相结合的方法，对DNA建库流程进行了改进和优化。tt
本发明设计了P5端DNA接头序列如SEQIDNo.1所示。tt
本发明设计了长度为6bp独特的24个标签(Index)序列，标签(Index)序列分别插tt入DNA接头(SEQIDNo.2)中，形成24个P7接头，P7端DNA接头序列如SEQIDNo.3-SEQIDttNo.26所示。单个P7接头分别与P5组合，共24对DNAPCR-free标签接头。每一对标签接头可tt以单独用于构建一种测序文库。tt
本发明还提供了通过连接酶的催化作用，直接将标签(Index)序列与DNA接头连tt接，不需要经过PCR反应而成功构建DNA测序文库。tt
本发明提供的构建基因组DNA的PCR-free测序文库的起始量为1μg-6μg之间。tt
本发明提供的CovarisM220超声打断仪将起始DNA进行打断。根据需要打断的片tt段大小，调整相应的打断时间。打断范围为插入片段insert300bp和insert500bp两种，打tt断选用CovarisM220超声打断仪，打断参数依次是占空因数(Dutyfactor)设置为20％，峰tt值功率(Peakincidentpower)设置为50W，循环破碎系数(Cyclesperburst)设置为200，tt工作温度设置为20℃；两种不同的打断范围insert300bp打断持续时间70秒(sec)，inserttt500bp打断持续时间42sec，基因组DNA打断起始量为1μg-6μg之间，打断体系为130μL。tt
本发明提供的多个纯化步骤中都用到磁珠进行纯化，磁珠选用AMpureXPttBreads，纯化时磁珠用量均为1.6×(体积为所纯化的溶液体积的1.6倍)，其中，步骤(2)也tt可选用PCR产物试剂盒(QIAGEN(货号Cat.no.28704))进行纯化回收。tt
本发明提供的末端修复步骤采用采用NEB修复试剂，体系为修复试剂5μL，10×修tt复缓冲液试剂10μL，片段化DNA溶液85μL，总体系为100μL，修复条件为20℃，60分钟(min)；4tt℃，终止(hold)。tt
本发明提供的末端加A反应步骤采用NEB加A反应试剂，体系为10×加A反应缓冲液tt试剂5μL，克列诺片段(KlenowFragment，3′→5′exo-)3μL，修复DNA溶液42μL，总体系50μL，tt连接条件为37℃，60min；4℃，hold。tt
本发明提供的接头与DNA连接的反应体系，根据起始DNA的量调整连接体系中接头tt工作液的量:起始DNA量>5μg，加5μL接头液；起始DNA量2-5μg，加3μL接头工作液。同时，设计tt了不同含标签(Index)序列的特异接头引物(共24种)，接头引物合成采用HPLC，接头引物配tt制时，先将上述P5和P7接头引物用无核酸水(NFW)各自稀释成母液(浓度为100μM)，再配制tt10×退火缓冲液(溶液终浓度100mMTris-HCl，500mMNaCl，10mMEDTA，PH8.0)，分别取P5tt和P7引物母液各9μL，再加入上述2μL10×退火缓冲液中配制成接头工作液(20μL)，然后在ttPCR仪器上进行退火，退火反应条件为：97.5℃，2min；6℃/min(每分钟降低6℃)至16℃；16tt℃，hold，-20度长期保存，下次使用时需重新退火，连接反应采用NEB加A反应试剂，体系为tt10×T4DNA连接酶缓冲液5μL，加A连连接试剂23μL，接头工作液(浓度为45μM)5μL，连接增强tt试剂1.5μL，上述加A的DNA溶液54μL，总体系88.5μL，所述此处工作液加入量根据步骤(1)中ttDNA起始量进行调整，当步骤(1)中DNA起始量大于5ug，接头工作液可加5μL；当DNA起始量为tt2-5μg之间可加3μL；当起始量小于2μg可加1.5μL，连接反应条件为20℃，120min；4℃，hold。tt
本发明采用2％琼脂糖凝胶电泳分离不同片段大小的DNA，然后用胶回收的方法进tt行片段筛选。2％琼脂糖凝胶电泳条件为：120V，1.5小时(h)。根据需要的insert片段大小，tt切取相应范围的胶进行回收。所述的insert300bp切取300-500bp之间，insert500bp切取tt500-700bp之间。tt
本发明的文库质检用Qubit3.0对文库样品进行浓度检测、安捷伦2100核酸分析仪tt对片段范围进行检测和QPCR进行定量检测。tt
采用上述技术方案后，与现有技术相比，本发明具有以下有益效果：tt
1、适用性广；本发明适用于所有已知或者未知参考基因组的物种，可为一般普通tt分子生物学实验室采用。tt
2、操作简单；本发明利用常规分子操作，设计兼容Illumina二代测序仪的特异接tt头，并采用磁珠和胶回收的方法得到范围比较集中的测序文库。tt
3、性价高；本发明主要采用NEB公司的常规建库试剂，自建文库构建流程，使得单tt次建库成本极大降低。tt
4、重复序列和碱基错误率少；本发明由于不需要经过PCR反应而成功构建成DNA测tt序文库，因此重复序列和碱基错误少。tt
附图说明tt
下面结合附图和实施例对本发明进一步说明：tt
图1本发明的PCR-free文库建库流程示意图。左侧为实验流程图，右侧为文库序列tt信息，两侧信息对应。箭头代表流程实验方向。N代表碱基，“….”代表省略的碱基。P5和P7代tt表接头引物序列，圆圈代表磷酸化修饰位点，P7和P5引物经过处理形成Y型接头。tt
图2是实施示例中PCR-free文库安捷伦2100核酸分析仪片段范围检测结果。纵轴tt(FU)表示峰值大小，横轴(bp)表示片段大小范围，两次峰值为参照峰，中间峰为样本检测tt峰。旁边附图电泳图，上下两条带为参照带，中间为样本检测带。tt
图3是实施示例中混样文库qPCR浓度定量结果图。(A)为文库QPCR定量检测溶解曲tt线，图中对应的线条的为本实例中样本和标准品的溶解曲线。横轴(Tm)为退火温度，纵轴(-ttRn)为标准化报告Rn的倒数。(B)为qPCR产物2％琼脂糖电泳图，左边为文库扩增后产物，M为tt100bpDNAladder，从上到到下依次为1500bp、1000bp、900bp、800bp、700bp、600bp、500bp、tt400bp、300bp、200bp、100bp。中间最亮带为500bp。tt
图4是实施示例中Hiseq2500测序仪下机数据质量FastQC评估结果。纵轴轴代表Qtt值大小，横轴代表数据长度。Q20处为临界点，代表单个碱基测序质量正确率99％。tt
具体实施例tt
下面结合说明书附图和实施例，对本发明的具体实施例做进一步详细描述：tt
这些实施例仅用于说明本发明而不用于限制本发明要求保护的范围，任何对本流tt程采用相似的策略而只是进行部分修改的内容也在本专利的保护范围。下列实施例中未注tt明具体实验条件和方法，实验条件和方法参照相关试剂说明书。tt
实施示例一个茄科灯笼草(DC)DNAPCR-free文库的制备。tt
用上述茄科灯笼草(DC)DNA样品构建两种插入片段长度(即insert300bp和ttinsert500bp)的PCR-free文库包括以下步骤：tt
1.DNA样品超声打断。分别取两管DC样本的DNA，每管5μg起始量。用CovarisM220tt超声打断仪分别进行打断，打断体系为130μL。打断参数是：tt
2.对超声打断的片段进行纯化回收。纯化采用PCR产物试剂盒进行纯化回收，采用ttQIAGEN(货号Cat.no.28704)，用无核酸水(NFW)进行洗脱，洗脱体积为85μL。tt
3.利用NEB末端修复试剂对纯化产物进行末端修复。修复体系为(共100μL)：tt
NEB修复试剂(货号#E7371)5μLtt
10×修复缓冲液(货号#E7372)10μLtt
回收DNA溶液85μLtt
修复条件为：20℃，60min；4℃，终止(hold)。tt
4.对末端修复产物进行纯化回收。纯化时采用AMpureXPBreads磁珠，磁珠用量tt为1.6×，即160μL，回收体系为42μL。tt
5.利用NEB末端加A试剂(货号#E6053)对上述纯化产物进行DNA末端加A反应。反应tt体系为(共50μL)：tt
10×加A反应缓冲液试剂5μLtt
KlenowFragment(3′→5′exo-)3μLtt
上述修复DNA溶液42μLtt
6.对上述加A产物进行纯化回收。纯化时采用AMpureXPBreads磁珠，磁珠用量为tt1.6×，即90μL，回收体系为54μL。tt
7.对上述纯化产物进行连接反应。本步骤可分为以下三个小步骤：tt
(1)配制接头工作液：先将P5(序列为SEQIDNo.1)和P7接头(序列为SEQIDttNo.3)用无核酸水(NFW)各自稀释成母液(100μL)，再配制10×退火缓冲液(100mMTris-ttHCl，500mMNaCl，10mMEDTA，PH8.0)。接头工作液(共45μL)包括：tt
(2)退火：在PCR仪器上进行97.5度，2min；6℃/min(每分钟降低6℃)至16℃；16℃，tthold。-20度长期保存，下次使用时需重新退火一次。tt
(3)连接反应：采用NEB快速连接反应试剂(货号NEB#E7373)，体系(共88.5μL)为：tt
盖紧盖子后，在震荡仪瞬时震荡3s，小离心机上快速离心5s，在PCR上20℃，tt120min，4℃hold。
8.用2％琼脂糖凝胶电泳对连接产物进行片段筛选，电泳条件为：120V，1.5小时tt(h)。采用胶回收的方法进行片段筛选。insert300bp切取300-500bp之间的片段进行回收，ttinsert500bp切取500-700bp之间的片段，胶回收两种片段的文库洗脱均为15μL。回收产物tt即为测序文库。tt
9.质检和上机测序。文库质检用Qubit3.0对文库样品进行浓度检测、安捷伦2100tt核酸分析仪对片段范围进行检测(图2)和QPCR进行定量检测(图3)。将浓度和片段大小符合tt要求的文库利用Illumina公司的二代测序仪进行高通量测序。下机数据用FastQC软件进行tt数据质量和GC含量分析，重复序列(duplication)用Bowtie2.0软件进行分析，重复数据仅tt占极少部分。tt
申请人文库制备结果表明本方法发明能够制备种一种基因组DNAPCR-free测序tt文库的制备方法，可以进行后续的高通量测序，测序数据也可用后续分析。据此，本发明专tt利普遍适于常规分子生物学实验室进行高GC、高AT含量的基因组DNA或者为了减少测序中tt碱基错误率和序列重复的文库构建。tt
以上所述，仅是本发明的较佳实施例而已，并非对本发明的技术范围作出任何限tt制，故凡是依据本发明的技术实质对以上实施例所作的任何细微修改、等同变化与修饰，均tt仍属于本发明的技术方案范围内。tt
