一种TaqDNA聚合酶变体及其获得方法和在基因组编辑中的应用

一种taq dna聚合酶变体及其获得方法和在基因组编辑中的应用
1.本技术为申请号2021103206684、申请日2021年3月25日、发明名称“高特异性taq dna聚合酶变体及其在基因组编辑和基因突变检测中的应用”的分案申请。
技术领域
2.本发明属于生物技术领域，具体涉及一种taq dna聚合酶变体及其获得方法和在基因组编辑中的应用。


背景技术：

3.公开该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不必然被视为承认或以任何形式暗示该信息构成已经成为本领域一般技术人员所公知的现有技术。
4.crispr/cas9技术能够仅通过一小段引导rna在特定位点进行便捷的基因组编辑，现已广泛应用于功能基因组学研究，并在涉及遗传变异的疾病治疗中具有巨大潜力。目的基因组修饰主要有三种类型，包括由于双链断裂导致的易错非同源末端连接(nhej)修复，该修复会引起indels随机突变；使用dna模板进行同源介导的修复(hdr)或直接通过碱基编辑引起的精确碱基改变；以及通过招募转录因子或染色质修饰因子进行基因调控。对于基因组编辑应用，通常需要评估给定crispr靶标的编辑效率，并在某些情况下，对获得的单细胞克隆进行基因分型。目前已经开发了几种方法，包括gef-dpcr，getpcr和(act-pcr)，它们可在pcr扩增过程中将发生编辑修饰的dna与野生型序列区分开。但由于taq酶或taqman探针对dna突变的鉴别能力有限，实验需要仔细优化才能得到较为准确的结果。使用修饰的荧光探针或使用比野生型taq酶具有更好的错配选择能力的增强型dna聚合酶变体，都可以提高pcr检测的准确性。dna聚合酶变体能在不需要任何探针或引物修饰的情况下进行可靠的遗传变异检测，因此是提高基因变异检测准确度最经济有效的策略。
5.聚合酶与引物/模板双链dna在小沟处的相互作用对于复制起始复合物的组装至关重要，然而，这些相互作用力是高度冗余的，超过了有效dna复制起始的最低需求，取代这些氨基酸以破坏相应的相互作用可以提高错配延伸中dna聚合酶的选择性。基于此原理的dna聚合酶理性进化主要集中在基序c中少数几个极性氨基酸和碱性氨基酸的替换上，比如，在12个氨基酸位点处进行功能突变并通过对分子改组产生的组合文库中进行筛选，鉴定出了选择性提高的taq变体。但是，所有这些dna聚合酶突变体的理性设计都是以提高3
‘
末端单核苷酸错配延伸选择性为出发点的。然而，基因组编辑导致的插入缺失突变在很大程度上是复杂的以及不可预测的，这导致pcr检测引物和含有indel基因组dna之间的错配类型也是极其多样化的。因此，基因组编辑研究非常需要一种新的dna聚合酶变异体，它具有更好的识别基因组修饰引起的引物-模板错配的能力，该taq变体将使基因组编辑频率检测和单细胞克隆基因分型等实验更加准确和方便。


技术实现要素：

6.针对现有技术存在的问题，本发明提供一种taq dna聚合酶变体及其获得方法和在基因组编辑中的应用。对野生型全长taq dna聚合酶进行了半理性的定向分子进化来提高其特异性。选取taq酶上与引物/模板复合物有直接相互作用的全部极性氨基酸进行逐个突变，获得40个taq变体，然后在这些变体及野生型序列的基础上进行广泛的随机诱变，生成taq突变体文库。在我们的qpcr筛选系统上，以基因组编辑indels质粒为模板，筛选出一系列具有高特异性的taq突变体，在crispr/cas9编辑效率评估和单细胞克隆基因分型中展现出了极大的优势，因此具有良好的实际应用之价值。
7.具体的，本发明涉及以下技术方案：
8.本发明的第一个方面，提供一种taq dna聚合酶变体，所述taq dna聚合酶变体选自下组的一个或多个位点发生突变：s577a、w645r、i707v、r405q、t569v、k354r、k531q、l441m、s543a、r630w、f692y、y719f、m4i、d371e、v518d、a798v、g32d、d238v、w398c、n485l、i503f、r771k、e284k、i614l、t588s、l789f、g59w、v155f、k508q、r229g、e255v、q489l、e90k、e132q、p369t、t513a、d151g、s515a、r741q、a294s、a675v、e688d、v740a、g173d、l500i、r37q、t140s、d365n、t140a、l538i、p10a、e303g、l484i、r492m、f272s、e794d、e170g、k508t、d578l、e818v、i799f、k206r、r229w、r249c、v390m、e404g、e267v、s577a、q680h、r328m、r469c、e159d、d181h、p387l、a61t、d91n、k100e、k131n、a777v、p194h、p369t、t514v、y719f、a118s、r435w、e708d、p6t、d177e、l252m、e465d、s699t、e135v、p316s、g422w、t385a、r137c、p685s、e818k、l828v、a414t、s515a、a600t、s36i、e171k、s576a、e57d、d222y、h28l、e112d、l245p、r630l、l351f、l657p、p816s，其中，氨基酸残基编号采用seq id no.1(野生型taq dna聚合酶的氨基酸序列)所示的编号。
9.所述taq dna聚合酶变体的氨基酸序列具有与seq id no.1相比至少80％的同源性；更优选地，具有至少90％的同源性；最优选地，具有至少95％的同源性；如具有至少96％、97％、98％、99％的同源性。
10.所述taq dna聚合酶变体中的突变位点的数量为1-6个，进一步优选为1-4个，如1个、2个、3个或4个。
11.所述taq dna聚合酶变体在seq id no.1所示的野生型taq dna聚合酶基础上进行突变，并且所述taq dna聚合酶变体选自下组中的突变体：
12.[0013][0014]
上述表格中的taq dna聚合酶变体是按照特异性从高到底排序，其中排在前十的变体为优秀变体，与野生型taq相比，它们检测indels错配的ct值至少多出7个循环，表示这些变体的选择性显著提高，其中突变体taq388拥有最好的选择性，提高了约23个循环。同时，taq388变异对源于indel和单核苷酸变异错配的pcr选择性均有极其显著的提高。在应用中，该taq变体显著提高了getpcr方法对单细胞克隆基因分型的准确性，同时也使as-qpcr snp基因分型成为更可行的方法。
[0015]
本发明的第二方面，提供一种多核苷酸分子，所述多核苷酸分子编码上述第一方面所述的taq dna聚合酶变体。
[0016]
本发明的第三方面，提供了一种重组表达载体，所述重组表达载体含有本发明第二方面所述的多核苷酸分子。
[0017]
具体的，所述重组表达载体通过上述多核苷酸分子有效地连接到表达载体上获得，所述表达载体为病毒载体、质粒、噬菌体、噬菌粒、黏粒、f黏粒、噬菌体或人工染色体中
的任意一种或多种；病毒载体可包括腺病毒载体、逆转录病毒载体或腺伴随病毒载体，人工染色体包括细菌人工染色体(bac)、噬菌体p1衍生的载体(pac)、酵母人工染色体(yac)或哺乳动物人工染色体(mac)。
[0018]
本发明的第四方面，提供了一种宿主细胞，所述宿主细胞含有本发明第三方面所述的载体或染色体整合有本发明第二方面所述的多核苷酸分子。
[0019]
所述宿主细胞可以为原核细胞或真核细胞。
[0020]
更具体的，所述宿主细胞是细菌细胞、真菌细胞或植物细胞中的任意一种或多种；
[0021]
其中所述细菌细胞为埃希氏菌属、农杆菌属、芽孢杆菌属、链霉菌属、假单胞菌属或葡萄球菌属内的任何种；
[0022]
更具体的，所述细菌细胞为大肠杆菌(如大肠杆菌dh5α)、根癌农杆菌(如gv3101)、发根农杆菌、乳酸乳球菌、枯草芽孢杆菌、蜡状芽孢杆菌或荧光假单胞菌。
[0023]
所述真菌细胞包括酵母菌。
[0024]
转基因植物包括拟南芥植株、玉米植株、高粱植株、马铃薯植株、番茄植株、小麦植株、油菜植株、油菜籽植株、大豆植株、稻植株、大麦植株或烟草植株。
[0025]
本发明的第五方面，提供了一种制备本发明第一方面所述taq dna聚合酶变体的方法，包括步骤：培养本发明第四方面所述的宿主细胞，从而表达出所述的taq dna聚合酶变体；和分离所述的taq dna聚合酶变体。
[0026]
本发明的第六方面，提供了一种试剂盒，所述试剂盒包含本发明第一方面所述的taq dna聚合酶变体。
[0027]
本发明的第七个方面，提供上述第一方面所述taq dna聚合酶变体、第二方面所述多核苷酸分子、第三方面所述重组表达载体、第四方面所述宿主细胞、第六方面所述试剂盒在如下任意一种或多种中的应用：
[0028]
1)基因组编辑检测(如基于crispr/cas9的基因组编辑)；
[0029]
2)基因突变检测(如单细胞克隆基因分型，和snp基因分型分析等)。
[0030]
以上一个或多个技术方案的有益技术效果：
[0031]
上述技术方案提供了一种高特异性taq酶变体及其在基因组编辑和基因突变检测中的应用。本发明对野生型全长taq dna聚合酶进行了半理性的定向分子进化来提高其特异性。选取taq酶上与引物/模板复合物有直接相互作用的全部极性氨基酸进行逐个突变，获得40个taq变体，然后在这些变体及野生型序列的基础上进行广泛的随机诱变，生成taq突变体文库。在我们的qpcr筛选系统上，以基因组编辑indels质粒为模板，筛选出一系列具有高特异性的taq突变体。其中，特异性最好的一个变体taq388在手掌区(s577a)和手指区(w645r和i707v)发生了三个氨基酸突变，在crispr/cas9编辑效率评估和单细胞克隆基因分型中展现出了极大的优势。此外，该变体在检测自然发生的遗传变异如snp也具有优良的性能，因此具有良好的实际应用之价值。
附图说明
[0032]
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
[0033]
图1为本发明高特异性taq定向进化策略的图解。
[0034]
(a)参与taq-引物/模板相互作用的40个极性氨基酸的示意图。极性氨基酸在序列上用箭头表示。(b)taq直接进化的原理和流程图。参与dna相互作用的40个氨基酸被单独突变，然后用易错pcr进行随机突变，并在筛选系统上对taq变体的活性和选择性进行了评估，筛选系统使用了26个在hoxb13基因sgrna靶1处含有indels的构建体，并给出了检测引物和退火区序列。与野生型taq相比，高选择性taq变异体具有更大的测试扩增ct值。
[0035]
图2为本发明高选择性taq变体的筛选
[0036]
(a)使用生长在含有iptg的lb琼脂平板中的菌落，评估40个taq变异体的酶活性以及在区分indel造成的错配中的选择性。ct值为45表示不再有聚合酶扩增活性。平均值
±
s.e.m，n＝3个技术重复。(b)在第一轮筛选中，对随机突变文库中的1316个转化子进行聚合酶活性和选择性评价。176个转化子保持了完整的聚合酶活性，并具有更高的特异性，并突出显示。(c)对176个转化子进行进一步的活性和选择性评估，选出了39个证实其选择性提高的转化子，并突出显示。(d)用纯化蛋白鉴定39个taq变异体。用箭头表示特异性最好的三个突变体。
[0037]
图3为本发明taq388对indel变异的选择性扩增能力分析。
[0038]
(a)基于taqman探针的qpcr系统中，taq388在qpcr反应种对hoxb13基因上模拟indels突变混合物造成的引物-模板错配的选择性评价。(b)taq388在sybr green qpcr体系中识别、选择上述indels能力评价。
[0039]
图4为本发明taq388识别单核苷酸错配的能力。
[0040]
(a)taq变异体对位于引物3
‘
末端最后一个核苷酸的引物-模板错配的敏感性评价，给出了引物和模板的序列。用匹配的模板计算相对pcr信号为100％。平均值
±
s.e.m，n＝3个独立的技术重复。(b)利用位于引物3
‘
末端的倒数第二个核苷酸的引物-模板不匹配来评价taq变异体的敏感性。平均值
±
s.e.m，n＝3个独立的技术重复。(c-d)在mcf7(c/c)(c)和t-47d(t/t)(d)基因组dna的等位基因特异性qpcr分析中，taq388区分乳腺癌风险snp rs4808611的不同等位基因的能力。
[0041]
图5为本发明taq388在getpcr检测基因组编辑中的应用。
[0042]
(a-b)比较taq388和野生型taq在qpcr扩增种对hoxb13基因上26个不同的indels的识别能力，taqman探针法(a)或sybr green法(b)检测携带每个indel的质粒。(c)比较taq388和野生型taq在hoxb13基因sgrna靶2进行基因组编辑的lenti-x 293t单细胞克隆的基因分型分析。所有20个克隆均含有先前确定的双等位基因插入缺失突变。(d)在dyrk1a基因sgrna靶1处进行基因组编辑的lenti-x 293t单细胞克隆的基因分型中,比较taq388和taq的特异性。所有编辑后的克隆都是双等位基因indel变异，sanger测序证实了这一点。检测引物中的观察碱基突出显示，pam序列“ngg”显示为浅色。ct值越大，酶的选择性越好。ct值为45表示无扩增信号。(平均值
±
s.e.m，n＝3个独立的技术重复)。
[0043]
图6为本发明taq变体在snp基因分型中的应用。
[0044]
(a-e)使用taq388通过qpcr对30个基因组dna样本上的5个snp位点rs2236007(a)、rs4808611(b)、rs11055880(c)、rs2290203(d)和rs2046210(e)进行的基因分型，并与野生型taq进行比较。使用公式：等位基因1％＝2-ct(allele1)/
(2-ct(allele1)
2-ct(allele2)
)计算每个等位基因的百分比含量。轴上的点是纯合基因型，轴之间的点是杂合基因型。taq388能成功区分每个基因型，但野生taq因其特异性差而无法判定样品的基因型。(f-j)taq388和野生型
taq对5个snp的等位基因特异性qpcr分析的终点荧光散点图。原点附近的灰点是用于对照的无模板扩增样品。
[0045]
图7为本发明高特异度taq的进化。
[0046]
(a)sanger测序确定的39个taq变异体的氨基酸突变，阴影标示的克隆为10个选择性最好的变体。(b)对从大肠杆菌中表达和纯化的39个taq突变体进行sds-page分析。(c)pcr扩增过程中野生型taq和taq388的突变频率，通过sanger测序分析确定。将taq388变异体扩增出的taq编码序列克隆到质粒中，对每个taq突变体的20个单细胞克隆进行测序，以鉴定突变。(d)利用taq388和野生型taq进行pcr扩增时产生的突变类型。
[0047]
图8为本发明taq变体对错配的敏感性。
[0048]
(a-c)在t-47d细胞(g/g)和vcap细胞(a/a)基因组dna的等位基因特异性qpcr分析中，taq388区分乳腺癌风险snp rs2236007的不同等位基因的能力。以及两种肿瘤细胞系中rs2236007位点基因型的sanger测序分析。
[0049]
(d)taq388与图中标示的五种商用qpcr检测预混产品比较区分indels的能力；taq388与图中标记的五种商品化qpcr主体混合物比较区分rs2236007的snp等位基因的能力。
[0050]
图9为本发明taq388与其他提高pcr选择性策略在snp检测中的比较。
[0051]
(a)用as-qpcr检测sw620基因组dna中tp53-g818a的遗传变异。taq388与3
‘
末端带有ddc的封闭型引物进行比较。(b)通过as-qpcr检测mda-mb-231基因组dna中tp53-g839a的变异。taq388与3
‘
末端带有ddc的封闭型引物进行比较。(c)用as-qpcr检测sw620基因组dna中tp53-g818a变异。将taq388与3
‘
末端含lna的引物进行比较。(d)用as-qpcr检测mda-mb-231基因组dna中的tp53-g839a。taq388与lna引物进行比较。(e)通过qpcr从mda-mb-231细胞中扩增出tp53-g839a。将taq388与3
‘
末端磷酸化的封闭型引物进行比较。
[0052]
图10为本发明在端点snp基因分型中野生taq的评估。
[0053]
(a-e)七种dna样品的sanger测序色谱，在对这五种样品进行qpcr snp基因分型时表现出广泛差异化的不同等位基因含量。sanger测序结果与qpcr结果高度一致。
具体实施方式
[0054]
应该指出，以下详细说明都是例示性的，旨在对本技术提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本技术所属技术领域的普通技术人员通常理解的相同含义。
[0055]
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本技术的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。下列具体实施方式中如果未注明具体条件的实验方法，通常按照本领域技术内的分子生物学的常规方法和条件，这种技术和条件在文献中有完整解释。参见例如sambrook等人，《分子克隆：实验手册》中所述的技术和条件，或按照制造厂商所建议的条件。
[0056]
以下通过实施例对本发明做进一步解释说明，但不构成对本发明的限制。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围。
[0057]
实施例
[0058]
1.实验材料和方法
[0059]
1.1 taq聚合酶的定点和随机诱变
[0060]
细菌表达taq聚合酶所用的质粒paktaq(addgene#25712)购自addgene网站。通过在paktaq的基础上进行定点突变pcr，逐个对参与taq酶-dna相互作用的40个极性氨基酸进行氨基酸替换(图1a)。在20μl定点突变pcr反应中含有4pmol定点突变引物和10μl 2x prime star max premix(takara)，pcr程序为98℃预变性15秒，然后98℃变性10秒，72℃延伸2分钟，循环25次，最后在72℃延伸5分钟。pcr产物中加入fastdigest dpni(thermo fisher scientific)在37℃切割2小时后直接用于转化dh5α感受态细胞，并将其涂布于含有氨苄青霉素的lb琼脂平板，在37℃培养箱中倒置培养过夜。第二天挑取单克隆接种到lb培养基中，37℃ 250rpm震荡培养过夜，从中提取质粒后用于sanger测序。
[0061]
将sanger测序确认的这40个突变体等比例混合，并将其以1：1的比例与paktaq混合，作为模板使用genemorph ii random mutagenesis kit(agilent technologies)通过易错pcr方法进行随机诱变。在25μl易错pcr反应体系中含有2.5μl 10xmutazyme ii reaction buffer,0.5μl 40mm dntp mix，1pmol上游和下游引物，0.5μlmutazme ii dna polymerase(2.5u/μl)和15ng模板质粒。pcr程序为95℃预变性2分钟，然后95℃变性30秒，60℃退火30秒，72℃延伸3分钟，循环10次，最后在72℃延伸10分钟。pcr产物通过ecori/sali双酶切克隆到原表达载体中。通过单克隆sanger测序确定转化子发生突变的频率，我们依据产品说明书对易错pcr的模板量和循环次数进行了调整，直到获得达到我们的要求的突变频率。
[0062]
1.2菌落qpcr筛选高特异性taq变体
[0063]
用随机突变文库质粒转化大肠杆菌dh5α感受态细胞，诱导taq突变体在含有氨苄青霉素和iptg的lb固体培养基中表达蛋白。为了确定不同taq变异体的活性和特异性，我们以26个基于pcdna3.1载体的带有模拟crispr/cas9基因编辑indel的hoxb13基因质粒作为pcr模板，采用菌落实时定量pcr方法进行筛选。在单管qpcr反应中包含两个扩增子，即检测扩增子和对照扩增子。检测扩增子的上游引物横跨模拟的基因组编辑位点，用来考察taq酶对于indel造成的引物-模板错配的选择性，检测扩增子使用一条fam标记的taqman探针。对照扩增则匹配临近未发生突变的序列，用来衡量taq酶变体的聚合酶活性是否受到了影响，对应一条vic标记的taqman探针，这里使用的引物都是根据getpcr策进行设计的，值得注意的是，质粒用fast digest noti(thermo science
tm
，cat#fd0593)线性化，以避免两个探针之间的荧光信号干扰。挑取生长在含有iptg的lb琼脂平板上的表达taq变体的单克隆菌落，加入10μl的1xtaq enzyme screening buffer(50mmtris-hcl[ph8.8]，16mm[nh4]2so4，0.1％[v/v]20，2.5μm mgcl2，每个dntp 0.25mm)混合均匀后，加入7μl到20μl的qpcr体系中。每个引物和探针的工作浓度分别为0.2μm和0.1μm。定量pcr程序为：95℃预变性5分钟，然后95℃变性30秒，68℃退火30秒，72℃延伸10秒，循环45次。检测扩增子ct值增加而对照扩增子ct值不变的taq变体时希望得到的特异性增加的变体。
[0064]
1.3 taq变体的纯化
[0065]
经过两轮菌落qpcr筛选，最终获得了39个改进的变体，通过sanger测序分析确定每个变体的突变氨基酸，并在大肠杆菌中进行表达和纯化。对于每个克隆，将其对应的100μ
l的过夜培养物，转接到含氨苄青霉素抗性的4ml lb液体培养基中，并在37℃和转速为250rpm活化约4h，当od600 nm达到0.8时，加入终浓度为1mm的iptg诱导蛋白表达，在37℃和转速为250rpm孵育12h。5000rpm离心3min收集菌体，用400μl缓冲液(50mmtris-hcl[ph7.9]，50mm蔗糖，1mmedta[ph8.0])重悬菌体沉淀，室温下5000rpm离心3min收集菌体。用200μl的预裂解液(50mmtris-hcl[ph7.9]，50mm蔗糖，1mmedta[ph8.0]，4mg/ml溶菌酶[amresco])，室温孵育15min。然后，将菌体悬液放入-80℃的冰箱中冷冻30min，然后将其置于室温下至完全解冻。在重复之前的冻融操作一次后，立即将此溶液放在37℃水浴中孵育15min。然后加入1μl的5mg/ml dnasei、1μl的1mcacl2和2μl的1mmncl2，混合均匀。继续在37℃孵育30min后，加入200μl的裂解缓冲液(10mmtris-hcl[ph7.9]，50mmkcl，1mmedta[ph8.0]，0.5％[v/v]20，0.5％[v/v]np40)并混合均匀，然后将该裂解液在75℃下孵育1h，随后在4℃下15000rpm离心10min，收集上清溶液。在其中加入0.12g固体(nh4)2so4，于4℃下旋转孵育30min。然后该溶液在4℃下15000rpm离心20min收集沉淀，将其重悬在300μl储存缓冲液中(50mmtris-hci[ph7.9]，50mmkcl，0.1mmedta[ph8.0]，1xpi，0.1％[v/v]20，50％[v/v]glycerol)，并在-20℃下储存。
[0066]
最后用sds-page电泳检测蛋白样品中taq突变体的含量，即将蛋白样品加入到由12％分离胶和5％浓缩胶组成的凝胶中，跑完电泳并用estaintml1蛋白染色剂(genscript)染色，并用quantum-st5(vilber lourmat,france)进行凝胶成像分析。
[0067]
1.4 taq388突变体的扩增保真度分析
[0068]
为了比较taq388和野生型taq的保真度，我们以质粒paktaq中的taq聚合酶编码序列为模板使用10x taq enzyme screening buffer进行pcr扩增。pcr产物用fastdigest ecori(thermo)和fastdigest sali(thermo)双酶切，然后插入同样双酶切的载体paktaq中。将连接产物转化大肠杆菌dh5α感受态细胞，挑选20个单细胞克隆进行sanger测序，计算每个克隆中扩增子序列的突变碱基数以得到突变频率。
[0069]
1.5 getpcr分析条件
[0070]
在基于sybr green的getpcr方法中，15μl的反应体系里含有7.5μl的2x taqbuffer，3pmol的每个引物，0.005ng质粒dna或3ng基因组作为模板，1μl的taq聚合酶。在qpcr仪rotor-gene q 2plex，qiagen上进行的分析，程序为：95℃初始变性5min，95℃变性30s，引物退火64～70℃，30s，72℃延伸10s，在96热循环仪(roche applied science,germany)上进行的分析则使用以下条件：95℃初始变性5min。
[0071]
在使用taqman探针的getpcr方法中，反应体系为20μl，包括2μl的10x taq enzyme screening buffer，0.1ng的质粒dna或10ng的基因组作为模板，4pmol的引物和2pmol的探针，1μl的taq聚合酶。在qpcr仪(rotor-gene q 2plex，qiagen)中采用以下程序进行实时pcr：在95℃下初始变性5min，然后在95℃变性30s，引物退火在64～70℃下退火30s，在72℃延伸10s，当使用96热循环仪(roche applied science,germany)时，使用了以下条件：初始变性循环(95℃，5min)，然后是45个pcr循环(95℃，15s，64-70℃，15s，72℃，15s)。
[0072]
1.6 taq388在indel检测中的选择性分析
[0073]
在sybr green和taqman探针法qpcr系统中检测taq388对于indel导致的引物-模板错配的选择性。这里使用的pcr模板是在taq变体筛选系统中用到的26个模拟indel的质
粒。这26个质粒混合在一起时模拟基因组编辑产生的indels混合物，而每个质粒单独作为模板则代表在基因组编辑实验中分离的带有纯合indel的单细胞克隆。对于taqman探针法qpcr检测，在一个20μl的反应体系中使用1对检测引物及1个对应的taqman检测探针，1对对照引物及1个对照taqman探针。sybr green方法有所不同，它不使用taqman探针，需要在两个反应管中分别进行检测扩增和对照扩增。
[0074]
在基因组编辑的实际应用场景中检测taq388的选择性时，使用了31个经过crispr/cas9基因组编辑的lenti-x293t单克隆细胞基因组dna，其中，20个单克隆细胞为hoxb13基因发生了双等位基因编辑，11个为dyrk1a基因发生双等位基因。未发生编辑的lenti-x 293t细胞系基因组用作两个系列的内部参照物，结合sybr green或taqman探针的qpcr用96仪器(roche)进行检测(图5c，d)。这里的pcr条件和程序参见getpcr分析条件部分。
[0075]
1.7 taq388在snp基因分型中的应用
[0076]
检测中使用了30份基因组dna样本，其中10个来自乳腺癌细胞系(mcf7，t47d，mda-mb-231，bt-474，bt-20，bt-549，sk-br-3，zr-75-1，mda-mb-468，mda-mb-453)，5个来自前列腺癌细胞系(lncap，du 145，pc3，22rv1，vcap)和4个来自其他类型的细胞系(hek293t,jurkat,hl-60,k562)，11个是来自研究者自身经过隐去个人信息处理的基因组dna。pcr反应中使用了针对5个snp位点(rs2046210[c/t]、rs2290203[c/t]、rs11055880[c/t]、rs4808611[c/t]和rs2236007[ga/ct])设计的等位基因特异性引物。在qpcr进行snp基因分型分析时，一方面我们依据qpcr得到的等位基因特异性的ct值，来计算在样品中该位点处每个等位基因的百分比含量，据此确定其基因型，以rs4808611为例，从qpcr反应中获得c等位基因特异性引物以及t等位基因特异性引物的ct值，然后使用公式分别计算两个等位基因的比例，c等位基因[c％＝2^-ct(c)/(2^-ct(c) 2^-ct(t)]和t等位基因[t％＝2^-ct(t)/(2^-ct(c) 2^-ct(t)]的比值；另一方面，我们可以将被测等位基因的荧光值直接绘制成散点图，直观地显示这些细胞系的基因型。这里的pcr条件和程序参见getpcr分析条件部分。作为对比，五种商品化产品也在rs2236007位点的基因分型中进行了应用，它们是2x ultra sybr mix、thunderbird sybr qpcr mix、select master mix、life power和2x t5fast qpcr，每种商品的扩增条件参照各自的产品说明书进行。
[0077]
1.8封闭型引物或lna引物的pcr
[0078]3‘
端含有ddc或磷酸基团的封闭引物和lna引物可用来提高等位基因扩增的选择性，我们针对sw620细胞基因组含有的纯合tp53-g818a位点和mda-mb-231细胞基因组上含有的tp53-g839a位点，设计等位基因特异性引物、对照扩增引物以及阻断引物，对它们提高pcr选择性进行了评估。在15μlqpcr反应体系中含有1xtaqbuffer，3pmol上、下游引物，以及0.005ng带有变异位点的pcr产物作为模板，pcr扩增程序为95℃预变性5分钟，然后45个循环的95℃15s,68℃15s,72℃15s，最后跟着一个标准的融解曲线程序。
[0079]
2.结果
[0080]
2.1高特异性taq定向进化的合理设计
[0081]
尽管5’核酸外切酶缺失的大片段(klentaq)能提高保真度和热稳定性，但是为了使最终的dna聚合酶变体同时适用于基于sybr green和基于taqman探针的qpcr分析，我们选择了全长的thermus aquaticus(taq)dna聚合酶(seq id no.1)而不是klentaq作为分子
进化的起始分子。科研人员认识到，替换掉与引物/模板复合物直接相互作用的氨基酸或影响结合口袋几何形状的氨基酸可以改变聚合酶的选择性。在以前的研究中，研究人员仅选择了一部分接触引物/模板的氨基酸进行突变。在本研究中，为了选出用于理性设计的候选氨基酸，我们调查了dna聚合酶开放形式和闭合形式的晶体结构，并选择了所有与引物/模板双链直接接触的40个极性氨基酸作为突变的靶标(图1a)。其中17个残基与引物链接触，24个残基与模板链接触，1个残基arg573同时与两者接触。针对这些选出的氨基酸，我们首先进行了定点突变，用含有非极性侧链的亮氨酸、丙氨酸或缬氨酸取代了40个极性氨基酸残基，同时尽量保持它们的空间几何形状不变。具体地说，将氨基酸n、r、q、e、k、y、d、m和h替换为l，将s和t分别替换为a和v(见下表)。由于通常氨基酸的极性侧链是直接参与接触的基团，非极性氨基酸残基的取代将有效破坏相应的相互作用，从而使taq聚合酶对引物/模板错配更加敏感，因此有望提高聚合酶在错配延伸方面的选择性。
[0082]
我们直接将生长在含有iptg的lb琼脂平板上的转化子用于高通量筛选，无需复杂的蛋白质纯化操作。首先在基于taqman探针的菌落qpcr系统上评估了40个taq变体的活性和选择性，该筛选系统使用26种在hoxb13基因上模拟indel的质粒作为模板。在该系统中，我们在一个反应管中设计了两个扩增子，其一是用来评估聚合酶选择性的检测扩增子，这里检测引物可退火到野生型dna序列上，这是发生基因组编辑产生indels的区域；另一个是用来评价聚合酶活性的对照扩增子，扩增引物退火到邻近区域(图1b)。26个indels会导致与检测引物的各种不匹配，与野生型taq相比，检测扩增子ct值的增加可以表明突变体选择性的增强。同时，如果对照扩增子ct值保持不变，则表明被测试的taq突变体活性不受突变的影响。
[0083]
我们发现其中9个变体严重丧失了聚合酶活性，包括r536l，y545l，r573l，n580l，n583l，y671l，n750l，q754l和h784l。与野生型taq相比，19个变体显示出更好的选择性，具有统计学意义，其中8个变体比野生型taq多出5个循环，表明这几个变体具有更好的选择性(图2a)。但是，即使保留完整活性并具有最高选择性的变体t206v也只能提高13.9个循环，还存在很大的局限性。
[0084]
[0085][0086]
2.2高选择性taq酶的广泛诱变分子进化
[0087]
进一步，我们在这40个变异体以及野生型taq的基础上做了广泛的随机突变，以筛选特异性更好的taq变体。用genemorph ii随机突变试剂盒将野生型taq表达载体与40个突变体混合后进行易错pcr，该试剂盒能以最小的突变偏倚引入合理水平的突变率。对于通过随机突变进行的定向蛋白质进化，通常每个构建体有2-7个核苷酸突变，对应1-3个氨基酸突变。通过调整输入模板量和循环次数，我们得到了taq基因编码区上含平均5.3个突变的taq突变体文库。然后将易错pcr产物克隆到原核表达质粒paktaq中，并直接将生长在含有iptg的lb琼脂平板上的单细胞菌落应用到qpcr筛选系统进行筛选。
[0088]
我们总共筛选了1316个克隆(图2b)，其中1001个克隆(76.1％)的扩增曲线在x轴上向右偏移并超过5个循环表示它们失去了大部分或全部的聚合酶活性，有101个克隆(7.7％)不仅保持了完整的活性，还表现出极高的选择性，甚至对于检测indel错配的扩增反应完全没有扩增信号。为了进一步证实这些高选择性taq变体的特异性，除了101个克隆外，我们扩大了范围，额外选择了75个克隆，这些克隆满足ct(ctrl)《14.5和ct(test)》30的标准(图2c中的色点)。这次，我们在含有lb琼脂平板的iptg上划线，收集直径大于2mm的菌落，并在qpcr筛选系统中进行了评估。我们发现只有62个菌落(35.2％)仍然符合ct(ctrl)《14.5和ct(test)》30的高特异性标准，这可能反映之前菌落qpcr体系稳定性较差。此时，我们选择了39个符合较高标准(ct(ctrl)《14.5和ct(test)》40)的克隆进行sanger测序，并对这些taq酶变体(见下表)在大肠杆菌中进行了蛋白表达和纯化，进一步用纯化的taq聚合酶(图2c中的圆点)进一步验证。有趣的是，我们发现这39个变体中，只有13个变体的氨基酸替换涉及了taq聚合酶和引物/模板复合物之间的直接接触(图7a)。
[0089]
[0090][0091]
2.3 taq变体的纯化及其选择性的验证
[0092]
如上所述，我们在大肠杆菌中表达并纯化了这39个特异性提高的taq变体。它们在
sds-page分析中表现出相似的纯度，表观分子量均为94kda(图7b)。我们在qpcr筛选系统中评估了这些变体在indels检测系统中的聚合酶活性和选择性，最后鉴定出10个优秀变体，与野生型taq相比，它们检测indels错配的ct值至少多出7个循环，表示这些变体的选择性显著提高(p《0.05)(图2d中的色点)，其中突变体taq388拥有最好的选择性，提高了约23个循环，在随后的实验中我们选择使用了该变体进行了系统的评估和应用。
[0093]
随后，我们通过sanger测序来评价taq388变异体在pcr扩增中的保真度。用taq388扩增taq编码序列，并将其克隆到原载体中，转化进大肠杆菌后，挑取单克隆进行sanger测序分析由于pcr扩增产生的dna突变。我们发现taq388的保真度提高了4.7倍(图7c)。值得注意的是，野生型taq发生了3种类型的突变，包括56.5％的转换、39.1％的颠换和4.4％的缺失，而taq388只产生了转换类型的突变(图7d)。简而言之，我们获得了多个增强型taq酶变异体，它在扩增indel引起的引物/模板错配时具有显著增强的选择性，并且在pcr扩增中保真度也提高了4.7倍。
[0094]
2.4增强型taq鉴别错配的能力
[0095]
然后，我们系统地评估了taq388变体对各种类型的引物/模板错配的区分能力。首先，在基于taqman探针的qpcr筛查系统上测试了其区分indels错配的能力。结果表明，taq388比野生型taq聚合酶的选择性提高了23个循环，这在筛选过程中已经得到了体现(图3a)。当在基于sybr green的qpcr系统中使用相同的引物和模板进行测试时，该变体区分indels错配的能力也大大提高，但程度低于基于taqman探针的系统(图3b)。进一步，我们系统地研究了该变体在引物3'末端最后一个或倒数第二个位置上识别单核苷酸错配的能力。为了产生单核苷酸错配，我们构建了在hoxb13 c.251g位置含有三种类型的单核苷酸变异的质粒作为qpcr模板，这三种类型包括c.251g》a、c.251g》t、c.251g》c(图4a、b)。我们使用4种仅3
′
末端核苷酸有差异的引物进行基于sybr green的qpcr分析，发现与野生型taq相比，taq388聚合酶变体在所有12种错配类型中都大大降低了来自错配模板的扩增信号(图4a)。类似地，使用具有不同的3
‘
末端倒数第二位核苷酸的引物进行qpcr分析表明，taq388变体也比野生型taq在引物3’末端的倒数第二位错配时具有更高的选择性(图4b)
[0096]
接下来，我们在基因组dna的实际应用场景中评估了taq变体对单核苷酸错配的扩增选择性。我们使用3
‘
末端靶向rs4808611位点的等位基因特异性引物，对snp位点基因型分别为c/c和t/t的mcf7细胞(图4c)和t-47d细胞(图4d)的基因组dna进行了qpcr分析。我们发现，对于这两个等位基因特异性引物，taq388变异体均比野生型taq具有更高的选择性。具体地说，对于t等位基因引物，来自c/c基因型的mcf7基因组dna的taq388变体的错配脱靶扩增强度减少了大约10个循环(图4c)，而对于c等位基因引物，来自t/t基因型t-47d基因组dna的扩增水平比taq(图4d)减少了10个循环以上。此外，我们在另一个snp位点rs2236007上观察到了类似的结果。具体地说，对于a等位基因特异性引物，用taq388变体对g/g基因型t-47d基因组dna的扩增水平减少了10.5个循环(图8a)，而对于g等位基因引物，来自a/a基因型vcap基因组dna的扩增水平比taq减少了多达7个循环(图8b)。
[0097]
此外，我们还将taq388变体与5种商用的基于sybr green的qpcr预混液产品进行了比较。值得注意的是，taq388聚合酶对indel造成的引物/模板错配表现出了比所列的所有商业产品更高的选择性(图8c)。此外，该变体在使用g/g和a/a基因型的基因组dna样本进行rs2236007位点的等位基因特异性pcr扩增中表现出比商品化产品更好的选择性(图8d)。
[0098]
2.5 taq388在基因组编辑单细胞克隆基因分型中的应用
[0099]
在功能基因组学研究中，我们通常需要在基因组编辑实验之后筛选大量子代个体或单细胞克隆，才能获得含有目标基因修饰的实验材料，而具有更高选择性的增强型taq聚合酶能够大大提高基因分型的准确性。因此，我们将taq388应用于单克隆的基因分型分析，模板为在筛选系统中作为模板使用的26个质粒。在基于taqman探针的qpcr分析中，使用野生型序列特异性测试引物，与野生型taq聚合酶相比，taq388区分插入/缺失的能力大大提高了，26个indel模板dna平均提高了16.9个循环(图5a)，其中23个indels模板甚至完全没有扩增信号。这表明taq388在识别和区分indels导致的引物/模板错配方面拥有极其出色的能力。当在基于sybr green的qpcr分析中，taq388在区分这26个indels与野生型的能力平均提高了10.7个循环，同样也显示出比野生taq更强的扩增特异性(图5b)。尽管不如在基于taqman探针的qpcr分析中那样出色，但基于sybr green的qpcr分析中野生型构建体和插入缺失构建体之间的最小ct值差异仍超过9个循环，这足以用来准确鉴定插入缺失序列的单细胞克隆。
[0100]
接下来，我们在实际应用场景中用基因组dna作为模板评估了taq388在31个单细胞克隆的基因分型分析中的表现，这些克隆是在lenti-x 293t上针对hoxb13基因和dyrk1a基因进行了crispr/cas9介导的基因组编辑7。sanger测序表明，其中二十个克隆在hoxb13基因上产生了双等位基因的indel突变，十一个单细胞克隆在dyrk1a基因上产生了双等位基因indel突变。qpcr基因分型分析结果表明，不管是在hoxb13基因上发生的基因编辑，还是在dyrk1a基因上发生的基因编辑，taq388都表现出了比taq聚合酶更好的区分indel序列与野生型序列的能力(图5c，d)。对于在hoxb13 sgrna靶标2上进行的基因组编辑，taq388和taq聚合酶区分indels与野生序列的能力的平均δct值分别为14.2和10.1个循环(图5c)。具体来说，在检测ht2-04克隆时，taq聚合酶仅给出4个循环的δct值，但是taq388在全部45个pcr循环运行结束都没有检测到有效扩增信号。关于在dyrk1a sgrna靶标1上的基因组编辑，由taq388和taq聚合酶测定的indels突变引起的δct值分别为9.5和2.6个循环(图5d)。这表明，taq388的应用可使基因组编辑检测更加准确可靠。
[0101]
2.6 taq388在snp基因分型中的应用
[0102]
作为第三代分子标记，snp位点具有许多优势，包括分布广泛和遗传稳定性高。它已广泛用于分子生物学，疾病预测和治疗等领域。但是，snp检测很大程度上也受限于dna聚合酶的特异性。因此，我们接下来测试了taq388在snp基因分型分析中的应用潜力，测试中使用了30个基因组dna样品，其中有19个来自从atcc购买的细胞系，11个样品来自发明人，经过随机打乱后排列和编号以隐藏个人信息。我们使用taq388进行等位基因特异性sybr green qpcr扩增，针对五个snp位点rs2236007，rs4808611，rs11055880，rs2290203和rs2046210进行基因分型分析，并通过sanger测序确定了这30个样品的snp基因型。
[0103]
我们采用了两种方法来确定样品的基因型。首先，我们通过图6附图说明中描述的方法，使用等位基因特异性ct值来计算该等位基因的比例，并据此确定snp基因型。理论上，对于等位基因1纯合的样本，计算得到的等位基因1和等位基因2的含量应该分别为100％和0％，杂合样本的两个等位基因的百分含量应介于这两个数值之间。对于snp位点rs2236007，使用taq388进行的qpcr分析表明，所有样本的snp基因型都能准确鉴定。其中a/a样本和g/g样本位于相应的坐标轴上，g/a样本位于两者之间(图6a)。出乎意料的是，这10
个g/a样品分布在一个相当分散而不是聚焦在50％附近的区域。我们检查了相应样品的sanger测序色谱图，发现这些样品的等位基因比例与sanger测序峰图中的相对峰高高度相关(图10a)。例如，sk-br-3细胞株具有最高的a等位基因比例，同时在sanger测序中也表现为a峰远远高于g峰，这表明，我们用taq388 qpcr基因分型计算出来的等位基因比例真实地反映了该样品的基因型。相反，在用野生taq聚合酶进行的qpcr分析中，所有样本点都堆积在第一象限，无法确定每个样本的基因型(图6a)。使用taq388聚合酶对剩下的四个snp位点rs4808611(图6b)、rs11055880(图6c)、rs2290203(图6d)和rs2046210(图6e)进行的基因分型，成功对每一个样品的snp基因型进行了确定。并且，杂合基因型样本的分散布局特性也与sanger测序中相应的峰高有很好的相关性(图10b-e)。
[0104]
常用的终点法snp基因分型技术使用taqman探针或等位基因特异性引物来区分不同的等位基因，现有情况下，为了准确进行snp基因分型，pcr对等位基因之间的选择性仍急需进一步的提高。因此，接下来我们评估了taq388在终点法基因分型方法中的应用，即在等位基因特异性pcr循环步骤结束后读取sybr green荧光，据此来确定样本的基因型。对rs2236007位点的分析结果表明，与野生型taq聚合酶相比，taq388的qpcr扩增能完全区别基因型为g/g，g/a和a/a的三组样品(图6f)，而野生型taq qpcr扩增后三种基因型的样品完全堆积在一起，无法区分。同样，我们还使用taq388聚合酶成功地对其他四个snp位点rs4808611(图6g)、rs11055880(图6h)、rs2290203(图6i)和rs2046210(图6j)进行了基因分型。
[0105]
本发明中，对全长taq聚合酶进行了半理性的定向进化，以提高其在pcr扩增中区分基因组编辑突变序列造成的引物-模板错配的能力。首先，我们对taq聚合酶上与引物/模板双链结构有直接相互作用的40个极性氨基酸进行了逐个定点诱变。然后，在这些变体以及野生型taq序列的基础上进行广泛的随机突变，生成全面的taq突变体文库。以带有indel的hoxb13基因质粒为pcr扩增模板，在qpcr平台上经过三轮筛选和验证，筛选出了多个特异性显著提高的taq变异体，其中带有s577a、w645r和i707v替换的taq388变异体表现最好。taq388变异对源于indel和单核苷酸变异错配的pcr选择性均有极其显著的提高。在应用中，该taq变体显著提高了getpcr方法对单细胞克隆基因分型的准确性，同时也使as-qpcr snp基因分型成为更可行的方法。
[0106]
之前所有提高dna聚合酶特异性的尝试都集中在辨别单核苷酸错配的能力上。本发明第一次针对基因组编辑indel导致的引物/模板错配，通过广泛的定向进化来获得性能更好的taq聚合酶变体。此外，作为起始分子，我们使用了全长taq聚合酶而不是其他研究中常用的klenow片段，这使得taq388变体不仅适用于基于sybr green的qpcr也适用于基于taqman探针的qpcr应用。
[0107]
并且，以前的研究大多是局限性的理性设计，侧重并局限于与引物/模板复合物相互作用的一部分极性氨基酸残基，以及它们之间进一步的简单组合应用。在这里，我们不仅包括了所有40个与引物/模板双链直接接触的极性氨基酸残基，并且在此基础上进一步进行了广泛的随机突变，以建立一个更全面的taq突变体文库。值得注意的是，在最终的39个变异体中，只有13个变体的氨基酸替换涉及了引物/模板接触的残基，并且，所有这些筛选出来的改良的变体都包含了未参与这种接触的氨基酸突变。此外，在我们最终获得的最好的10个变体中，多达5个taq变体的氨基酸突变完全不涉及那些参与酶/引物/模板相互作用
的氨基酸。这表明，这些引物/模板非接触性氨基酸的替代也有助于提高dna聚合酶的选择性，为dna聚合酶的进化提供了新的方向。
[0108]
当应用于检测基因组编辑突变时，taq388变体显示出极强的区分基因编辑序列和野生型序列的能力。这将使getpcr在基因组编辑实验中检测基因组编辑效率和单细胞克隆的基因分型变得更加准确和方便。当应用于检测那些自然发生的遗传变异时，taq388在as-qpcr分析中也显示出出色的snp等位基因识别能力。我们得益于taq388在pcr反应中的优秀的等位基因选择能力，实现了两种简单而有效的snp基因分型方法，即利用等位基因特异性ct值计算等位基因比例或绘制等位基因特异性pcr扩增的终点荧光散点图。对于这两种方法，三种基因型的样品都可以实现轻松、准确的鉴定。
[0109]
综上所述，通过半理性的定向进化，我们开发了多个taq聚合酶变体，它们对来自基因组编辑indels的引物/模板错配具有显著提高的选择性，其中最好的突变体taq388在基因组编辑测试和遗传变异检测方面显示出极大的潜力，该策略的成功为dna聚合酶的进化提供了新的思路。
[0110]
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。