一种构建大片段文库的方法及其应用
技术领域
本发明属于分子生物学领域,具体涉及一种大片段DNA文库的构建方法,由该方法构建的大片段DNA文库以及该大片段DNA文库在测序中的应用。
背景技术
对基因组序列未知或没有近源物种基因组信息的某个物种,对其不同长度基因组DNA片段及其文库进行序列测定,然后用生物信息学方法进行拼接、组装和注释,从而获得该物种完整的基因组序列图谱,称为基因组从头测序,也叫de novo测序。在基因组学飞速发展的今天,de novo测序与比较基因组方法联合可以用来探索该物种的起源和进化,研究其生长发育、形状生产以及环境适应的分子机制,是快速了解一个物种的途径之一。一个物种基因组序列图谱的完成,也将带动这个物种下游一系列研究的开展,可以构建该物种的基因组数据库,为该物种的后基因组学研究搭建一个高效的平台;为后续的基因挖掘、功能验证提供DNA序列信息。
高通量测序也叫二代测序,它可以对上千万条DNA片段同时进行测序,通量、时效和单碱基成本上较Sanger测序有着突破性的提高,这些优点使得目前de novo测序主要由二代测序来完成。二代测序的缺点是读长短,二代测序仪里读长最长的罗氏454也仅能进行最长400bp的DNA测序,这使得二代测序在进行de novo项目时要先构建小片段文库,根据不同小片段上的重合序列将其测序结构拼接成大小不一的contig(片段重叠群)。这些contig进一步拼接成基因组,则需要一种称为大片段文库的DNA文库,来确定不能重叠的contig在基因组中的位置和距离。
大片段文库是一种从一定片段大小的随机打断基因组DNA构建过来的文库,这些片段大小通常从2k到200k不等。由于这些片段大小已经远超二代测序读长的极限,因此在测序时不能采用全部测通的方法,而是采用双端配对测序,仅对这些片段的两端进行成对测序,以此判断测出的成对序列在基因组上的距离。为达到这个目的,构建大片段文库时,通常的程序可以如非专利文献1中报道的那样,先回收固定大小的DNA片段,然后将片段末端加上带生物素标记的接头,处理后的DNA片段被环化,再进行随机打断,使用链霉亲和素磁珠将生物素标记的打断片段钓取出来,建成测序文库进行双端配对测序。大片段文库的插入片段来源于最初构建文库的固定大小的DNA片段的两端,所以测出来的两端序列在基因组的上距离为构建文库前的DNA片段长度减去大片段文库的插入片段长度,以此来通过大片段文库辅助基因组装。
在大片段文库的构建过程中,希望得到的待测序片段是含有大片段双末端序列信息的重组片段,这样的重组片段一般由原始DNA片段环化产生并带有生物素标记。因此通常用生物素钓取的方式富集。由于原始DNA片段并不能完全环化,部分未环化的片段末端也含有生物素标记,如果不能很好地去除含有生物素的未环化的DNA片段,会造成钓取得到的含有生物素标记的待测序片段混有不含有大片段双末端序列信息的那些片段,导致双端测序文库测序数据量的浪费甚至影响最终结果。因此,如何在不影响环化DNA片段产率的前提下高效地消化掉被生物素标记的未环化DNA片段,成为大片段文库构建成功的一个关键步骤。目前的线性消化方法大多采用plasmid safe ATP-Dependent DNase来对线性DNA进行消化。plasmid safe ATP-Dependent DNase可以将双链DNA消化成脱氧核糖核苷酸,而对缺刻和环状双链DNA无活性。该消化方法不能特异地消化含有生物素标记的未环化DNA片段,是一种无特异性的线性消化方法。而且,采用plasmid safe ATP-Dependent DNase的消化方法需要添加ATP作为反应动能,以及核酸外切酶I作为辅助才能实施,会带来较高的人力、试剂和试剂储存成本。
非专利文献
1.Illumina.(2009)Mate Pair Library v2 Sample Preparation Guide For 2-5kb Libraries.
2.Filip Van Nieuwerburgh,Ryan C.Thompson,Jessica Ledesma,etal.Illumina mate-paired DNA sequencing-library preparation using Cre-Loxrecombination.Nucleic Acids Research.2012,40(3):e24.
发明内容
本发明的发明人为解决上述技术问题进行了深入研究,结果发现:通过先将DNA大片段加上5'端磷酸化接头,再对加有5'端磷酸化接头的DNA大片段进行环化,然后采用具有5'端磷酸化DNA特异性消化功能的核酸外切酶作为消化酶消化未环化的DNA大片段。由于消化酶具有5'端磷酸化DNA特异性消化功能,因此可以对加有5'端磷酸化接头的未环化的DNA大片段进行特异的消化,同时不需要添加ATP作为反应动能,也不需要核酸外切酶I作为辅助即可实施,从而完成了本发明。
即,基于上述认识的本发明包括:
1.一种大片段DNA文库的构建方法,包括:
步骤B:将DNA大片段进行末端修复,得到平末端DNA大片段;
步骤C:将步骤B得到的平末端DNA大片段加5'端磷酸化接头,得到加接头DNA大片段,所述5'端磷酸化接头带有生物素标记;
步骤D-1:将步骤C得到的加接头DNA大片段进行环化,得到环化DNA大片段与未环化DNA大片段的混合物;
步骤D-2:将步骤D-1得到的混合物进行线状DNA消化,得到环化DNA大片段,所述线状DNA消化采用具有5'端磷酸化DNA特异性消化功能的核酸外切酶。
2.根据项1所述的大片段DNA文库的构建方法,还包括
在所述步骤B之前进行的步骤A:将希望进行测序的DNA片段进行打断,得到DNA大片段;
在所述将步骤D-2之后进行的步骤E:将所述环化DNA大片段进行打断,得到测序用DNA片段;
步骤F:将测序用DNA片段进行片段捕获,采用链霉亲和素化固体载体钓取带有生物素标记的测序用DNA片段;
步骤G:将带有生物素标记的测序用DNA片段进行末端修复,得到平末端DNA片段;
步骤H:将步骤G中得到的平末端DNA片段进行加接头,得到加接头DNA片段;和
步骤I:将加接头DNA片段进行PCR扩增,得到扩增产物,从而构建大片段DNA文库。
3.根据项1或2所述的大片段DNA文库的构建方法,所述具有5'端磷酸化DNA特异性消化功能的外切酶为λ核酸外切酶和/或T7核酸外切酶。
4.根据项1-3任一项所述的大片段DNA文库的构建方法,所述DNA大片段的片段大小为1k~200kbp,优选1.5k~30kbp,更优选2k~20kbp,更优选5k~10kbp。
5.根据项1-4任一项所述的大片段DNA文库的构建方法,所述测序用DNA片段的片段大小为400~600bp。
6.根据项1-5任一项所述的大片段DNA文库的构建方法,所述5'端磷酸化接头为包含LoxP位点的接头。
7.根据项1-6所述的大片段DNA文库的构建方法,所述固体载体为磁珠。
8.一种大片段DNA文库,其是通过项1-7中任一项所述的大片段DNA文库的构建方法构建得到的。
9.一种大片段DNA文库的测序方法,其以项8所述的大片段DNA文库作为对象进行测序。
10.根据项9所述的测序方法,其中,所述测序为双端测序。
11.根据项9所述的测序方法,其中,所述测序利用Illumina平台进行。
12.一种用于构建大片段DNA文库的试剂盒,包括带有生物素标记的5'端磷酸化的接头,以及具有5'端磷酸化DNA特异性消化功能的核酸外切酶。
13.根据项12所述的试剂盒,还包括选自下组试剂中的至少一种或两种以上:用于末端修复的试剂,用于加接头的试剂,用于DNA片段环化的试剂,用于生物素钓取的试剂,用于DNA片段扩增的试剂。
与现有技术相比,本发明的有益效果是:采用具有5'端磷酸化DNA特异性消化功能的核酸外切酶,如Lambda Exonuclease(λ核酸外切酶)和/或T7核酸外切酶,选择性地对5'端磷酸化的双链DNA片段进行消化,可以特异性地除去经生物素标记的未环化的DNA片段,避免钓取生物素标记的片段中混入的那些经生物素标记但未环化的片段,即不含有末端配对信息的片段。此外,采用具有5'端磷酸化DNA特异性消化功能的核酸外切酶进行消化的步骤无需使用ATP或者引入其他酶,因此本发明提出的线性DNA消化系统从操作简易度和成本上都要优于传统的plasmid safe ATP-Dependent DNase线性DNA消化体系。
附图说明
图1是实施例1得到的文库的插入片段大小分布图。
具体实施方式
本说明书中提及的科技术语具有与本领域技术人员通常理解的含义相同的含义,如有冲突以本说明书中的定义为准。
首先,在一个方面,本发明提供一种大片段DNA文库的构建方法(本发明的文库构建方法),包括:
步骤B:将DNA大片段进行末端修复,得到平末端DNA大片段;
步骤C:将步骤B得到的平末端DNA大片段加5'端磷酸化接头,得到加接头DNA大片段,所述5'端磷酸化接头带有生物素标记;
步骤D-1:将步骤C得到的加接头DNA大片段进行环化,得到环化DNA大片段与未环化DNA大片段的混合物;
步骤D-2:将步骤D-1得到的混合物进行线状DNA消化,得到环化DNA大片段,所述线状DNA消化采用具有5'端磷酸化DNA特异性消化功能的核酸外切酶。
在传统的大片段DNA文库的构建方法中,是将未环化的DNA片段,无论加生物素标记与否,都进行线性消化。而本发明的文库构建方法与传统方法不同在于:将带有生物素标记的未环化的DNA片段进行线性消化。通过上述步骤,可以使经过生物素标记的未环化的DNA片段得到特异性的消化。
上述DNA大片段的末端修复可以通过任何本领域技术人员已知的方法来实现,例如,可以通过使用具有上述功能的DNA聚合酶(例如T4DNA聚合酶)来进行。
需要说明的是,在本说明书中对于步骤B中的“DNA大片段”的大小没有特殊限制,从构建De Novo测序文库的需求或者其他任何需要进行大片段文库的要求来看,可以是1k~100kbp左右的片段,优选是1.5k~30kbp左右的片段,更优选是2k~20kbp左右的片段,更优选是2k~10kbp左右的片段。
上述平末端DNA大片段的加5'端磷酸化接头可以通过任何本领域技术人员已知的方法来实现,例如,可以通过例如使用具有平末端连接功能的DNA连接酶(例如T4DNA连接酶、T3DNA连接酶)来进行。5'端磷酸化的接头本领域技术人员可以采用任何已知的5'端经过磷酸化修饰的接头(例如,非专利文献2中披露的LoxP接头),也可以通过本领域技术人员已知的方法从商业途径可获得的产品制备,还可以原样使用商业途径可获得的产品。上述加接头DNA大片段的环化可以通过任何本领域技术人员已知的方法来实现,例如,可以通过使用具有环化功能的DNA环化酶(例如Cre重组酶)来进行。
上述将混合物中未环化DNA大片段的消化,可以通过具有5'端磷酸化DNA消化活性的外切酶或消化5'端磷酸化的DNA双链片段的方法来实现。具有5'端磷酸化DNA消化活性的外切酶,优选地可以采用例如,λ核酸外切酶,T7核酸外切酶。
优选地,本发明的文库构建方法还包括:
在所述步骤B之前进行的步骤A:将希望进行测序的DNA片段进行打断,得到DNA大片段;
在所述将步骤D-2之后进行的步骤E:将所述环化DNA大片段进行打断,得到测序用DNA片段;
步骤F:将测序用DNA片段进行片段捕获,采用链霉亲和素化固相载体钓取带有生物素标记的测序用DNA片段;
步骤G:将带有生物素标记的测序用DNA片段进行末端修复,得到平末端DNA片段;
步骤H:将步骤G中得到的平末端DNA片段进行加接头,得到加接头DNA片段;和
步骤I:将加接头DNA片段进行PCR扩增,得到扩增产物,从而构建大片段DNA文库。
上述步骤A、E、F、G、H和I均可采用本技术领域的常规方法来进行。
需要说明的是,在本发明中对于步骤E中的“测序用DNA片段”没有特殊限制,从测序仪可接受的角度来看,优选是10~1000bp左右,更优选是20~800bp左右,更优选是30~750bp左右,更优选是40~700bp左右,更优选是50~650bp左右,更优选是100~600bp左右,更优选是150~550bp左右,更优选是300~400bp左右的DNA片段。
优选地,在本发明的文库构建方法的各个步骤之间例如步骤A与步骤B之间、步骤B与步骤C之间、步骤C与步骤D-1之间、步骤D-1与步骤D-2之间、步骤D-2与步骤E之间、步骤E与步骤F之间、步骤G与步骤H之间、步骤H与步骤I之间,和/或在步骤I之后可以加入对DNA片段进行纯化的步骤。
该纯化步骤可以采用本技术领域的常规方法来进行,例如可以通过使用纯化磁珠来进行。
在一个方面,本发明提供一种大片段DNA文库(本发明的文库),其可以采用例如本发明的文库构建方法来构建。
此外,在一个方面中,本发明提供一种大片段DNA文库的测序方法(本发明的测序方法),其中,以本发明的大片段DNA文库作为对象进行测序。
除了以本发明的大片段DNA文库作为对象之外,本发明的测序方法可以采用本技术领域的常规方法来进行。优选地,本发明的测序方法可以采用双端测序,例如利用Illumina平台(例如HiSeq2500或NextSeq500)进行的双端测序。但在一次测序可以将待测序DNA片段测通的情况下,也可以是单端测序。
本发明的建库方法可以例如使用试剂盒来实施,因此,在另一个方面中,本发明提供一种用于构建大片段DNA文库的试剂盒(本发明的建库试剂盒),其可用于实施本发明的文库构建方法,其包括带有生物素标记的5'端磷酸化的接头,以及具有5'端磷酸化DNA特异性消化功能的核酸外切酶。
优选地,带有生物素标记的5'端磷酸化的接头选自下组中的至少一种或两种以上:包含loxp位点(序列)的接头、包含限制性内切酶位点的接头、平末端接头、粘末端接头、双链接头、单链退火接头、单链发卡回环接头。Loxp位点即LoxP序列,由两部分组成,有两个13bp反向重复序列和中间间隔的8bp序列组成,8bp的间隔序列同时也确定了LoxP的方向。具有5'端磷酸化DNA特异性消化功能的核酸外切酶选自下组中的至少一种或两种以上:λ核酸外切酶、T7核酸外切酶、核酸外切酶Ⅰ、核酸外切酶Ⅱ、核酸外切酶Ⅱ、核酸外切酶Ⅳ、核酸外切酶Ⅴ、核酸外切酶Ⅵ、核酸外切酶Ⅶ、核酸外切酶Ⅷ、T5核酸外切酶、S1核酸酶、绿豆核酸酶、微球菌核酸酶、BAL-31核酸酶、RecJf核酸外切酶、脱氧核糖核酸酶Ⅰ、Ⅰ型限制性内切酶、Ⅱ型限制性内切酶、Ⅲ型限制性内切酶。
优选地,本发明的建库试剂盒还包括选自下组试剂中的至少一种或两种以上:用于末端修复的试剂,用于加接头的试剂,用于DNA片段环化的试剂,用于生物素钓取的试剂,用于DNA片段扩增的试剂。上述试剂可以采用任何本领域技术人员已知的试剂,例如,T4DNA聚合酶、Klenow片段、Klenow缓冲液、DNA连接酶缓冲液、DNA连接酶、Taq酶、dNTP、T4多聚核苷酸激酶、以及T4多聚核苷酸激酶缓冲液。
在本发明的建库试剂盒中,各试剂或装置优选单独包装,但在不影响本发明的实施的前提下,也可以混合包装。
实施例
以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例1大片段DNA文库构建
1.DNA样本的大片段处理
1.1 取2mL健康人全血按血液基因组DNA提取系统(0.1-20mL)(天根生化科技<北京>有限公司)操作说明书进行提取,得到DNA样本。取10μgDNA样本,使用HydroShearPlusDNA片段化仪按照说明书进行DNA样本的大片段打断处理,目标大片段长度为10kb。
1.2 配制0.8%的琼脂糖凝胶,使用TIANGEN公司1Kb DNA Ladder为分子量标准,100V电泳2小时。电泳结束后取出凝胶,放入含EB染料的TAE中染色20分钟。紫外照射下对9-12kb左右的片段进行切胶。
1.3 将切下的胶块放入已称重的干净的2.0mL离心管中,用琼脂糖凝胶回收试剂盒进行胶纯化回收DNA。
1.4 回收得到大片段DNA样本,用于下步反应或置于-80℃储存。
2. 末端修复
2.1 在1.5mL的离心管中配制末端修复反应体系:
dNTP Mix(dNTP混合物)(10mM each)是含有dATP,dCTP,dGTP和dTTP的钠盐的预混溶液,各自的浓度分别为10mM,总浓度为40mM(pH7.5)。
2.2 在Thermomixer C恒温混匀仪(EPPENDORF公司,以下简称恒温仪)中20℃温浴30分钟。
2.3 用AgencourtAMPure XP beads(BECKMAN COULTER公司,以下简称纯化磁珠)100μL纯化回收DNA,用35μL的EB缓冲液洗脱得到末端修复DNA样品。
3 接头连接
本步的接头采用非专利文献1中的5'端磷酸化的双链LoxP接头(即非专利文献中1的LoxP adapter oligos)。
3.1 在1.5ml的离心管中配制加接头反应体系:
3.2 在恒温仪中25℃温浴15分钟。
3.3 用1.0纯化磁珠回收纯化反加接头应体系中的DNA,用40μL的EB缓冲液洗脱得到加接头DNA样品。
4. 缺刻平移反应(Fill-In Reaction)
4.1 在1.5ml的离心管中配制缺刻平移反应体系:
4.2 在恒温仪中50℃温浴15分钟。
4.3 采用2.0 Flurometer(Life Technologies,CA,USA,以下简称Qubit2.0)检测步骤4.2的结果产物浓度,对浓度进行质量检测,确定回收DNA量大于400ng后进行下步反应。
5. DNA环化
5.1 在1.5mL的离心管中配制环化反应体系:
5.2 在恒温仪中37℃温浴50分钟,后70℃温育10分钟。
6. 线性DNA的消化
6.1 在反应完的步骤5.1的反应体系中依次加入以下试剂:
10×Lambda Exonuclease Buffer 6μL
Lambda Exonuclease 4μL
6.2 在恒温仪中37℃温浴30分钟消化线性DNA。
6.3 在恒温仪中75℃温浴10分钟酶失活后置冰上。
6.4 加入2μL EDTA(0.5M)到上述反应液中,充分终止消化作用,得到环化DNA。
7. 环化DNA片段化
7.1 使用Bioruptor DNA打断仪按下表打断程序打断环化DNA,得到打断DNA。
打断参数:
目标片段大小 |
ON/OFF时间(秒) |
循环数 |
300-400bp |
15/30 |
10 |
具体操作方法见《Bioruptor标准操作流程》
7.2 用1×纯化磁珠回收纯化步骤7.1中的打断DNA,溶于约50μL的EB缓冲液中用于下步反应或置于-80℃储存。
8. 纯化生物素标记的DNA
8.1 振荡重悬M-280 Streptavidin magnetic beads(链霉亲和素磁珠)。
8.2 吸取20μL重悬的磁珠置于1.5mL离心管,将离心管放置在磁分离架上等待1分钟,小心吸取弃上清。
8.3 用50μLBead Binding Buffer(磁珠结合缓冲液)洗涤磁珠。小心的重悬沉淀,将离心管放置在磁分离架上,等待1分钟,弃上清。重复此步骤一次。
8.4 用50μL磁珠结合缓冲液重悬磁珠。
8.5 加入来自步骤7.2的结果产物50μL,在恒温仪中20℃温浴15分钟(每2分钟震荡15秒,600rpm)。
8.6 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的Bead WashBufferⅠ(磁珠洗涤缓冲液)洗涤磁珠三次。
8.7 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的EB缓冲液洗涤磁珠两次。
8.8 移去最后一次洗涤的EB缓冲液,使用75μL的EB缓冲液重悬磁珠。
9. 末端修复
9.1 按照下表配制末端修复反应体系:
9.2 在恒温仪中20℃温浴30分钟(每2分钟震荡15秒,600rpm)。
9.3 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的磁珠洗涤缓冲液洗涤磁珠三次。
9.4 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的EB缓冲液洗涤磁珠两次。
9.5 移去最后一次洗涤的EB缓冲液,使用32μL的EB缓冲液重悬磁珠。
10.末端加“A”
10.1 按照下表配制加“A”反应体系:
10.2 置于恒温仪中37℃温浴30分钟(每2分钟震荡15秒,600rpm)。
10.3 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的磁珠洗涤缓冲液洗涤磁珠三次。
10.4 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的EB缓冲液洗涤磁珠两次。
10.5 移去最后一次洗涤的EB缓冲液,使用18μL的EB缓冲液重悬磁珠。
11.Illumina测序接头连接:
11.1 按照下表配制接头连接反应体系:
PE Adapters的碱基序列如下所示:
11.2 置于恒温仪中20℃温浴15分钟(每2分钟震荡15秒,600rpm)。
11.3 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的磁珠洗涤缓冲液洗涤磁珠三次。
11.4 将离心管放置在磁分离架上,等待1分钟,舍弃上清,用200μL的EB缓冲液洗涤磁珠两次。
11.5 移去最后一次洗涤的EB缓冲液,使用21μL的EB缓冲液重悬磁珠。
12.文库扩增
12.1 按照下表配制文库扩增反应体系:
引物1及引物2的碱基序列如下所示:
12.2 PCR反应的程序设定如下:
12.3 扩增产物使用琼脂糖电泳切胶回收400-600bp范围内的片段,获得大片段测序文库。
12.4 文库质检,片段大小与切胶回收一致,浓度测试符合上机测序要求。
13 上机测序及测序结果分析
13.1 将检测合格的步骤12.3得到的大片段测序文库,在HiSeq 2500测序平台上运行双端测序程序(PE150),得到下机数据如表1所示。
表1
通过表1可知:Clean reads占Original reads比例较高,Q30较高,整体测序结果较好。
13.2 截取100bp(PE100)的序列reads进行分析,通过去loxp接头的软件DeLoxer(非专利文献2),将Clean reads数据进行分类,分类结果如表2所示。
表2:
Mate-paired数据为通过DeLoxer软件处理后得到的大片段文库中真实的末端配对片段的数值和比率,其是在没有参考基因组情况下可以用于构建基因组骨架的有效数据。本实施例中,有效数据率为可比对数据(Mate-paired pair number)占原始数据(Cleanreads number的一半)的百分比。如表2所示的步骤12.3所得文库的有效数据率为37%,而这个比率在参考文献2中仅有28.62%。
13.3 将步骤13.2所得Mate-paired数据比对到人类参考基因组HG19,每对Mate-paired数据之间的距离对应一个大片段文库中原始大片段的大小(即插入片段大小,Insert size),Insert size出现的频次强度如图1所示。
由图1可知,本实施例构建的文库的插入片段主峰为9035bp,在期望建库大小10kb±10%以内,与目标大片段长度10kb基本一致。插入片段主峰±20%的范围内插入片段比例为83.6%,即所得文库的插入片段大小主要分布在目标大片段长度10kb附近。
上述说明示出并描述了本发明的优选实施例,如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
序列表
<110> 安诺优达基因科技(北京)有限公司
<120> 一种构建大片段文库的方法及其应用
<130> 1602SGCN
<160> 4
<170> PatentIn version 3.3
<210> 1
<211> 32
<212> DNA
<213> 人工序列
<400>
GATCGGAAGAGCGGTTCAGCAGGAATGCCGAG 32
<210> 2
<211> 32
<212> DNA
<213> 人工序列
<400>
ACACTCTTTCCCTACACGACGCTCTTCCGATC 32
<210> 3
<211> 38
<212> DNA
<213> 人工序列
<400> 引物1
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCT 38
<210> 4
<211> 38
<212> DNA
<213> 人工序列
<400> 引物2
CAAGCAGAAGACGGCATACGAGATCGGTCTCGGCATTC 38