新型植物碱基编辑器pAYBEs

新型植物碱基编辑器paybes
技术领域
1.本发明属于分子育种领域，具体涉及新型植物碱基编辑器paybes，特别涉及可实现植物基因组中单个碱基由腺嘌呤(a)颠换为胞嘧啶(c)或胸腺嘧啶(t)的植物a-to-y单碱基编辑器paybes及其应用。


背景技术：

2.在植物中许多重要的农艺性状是由基因组序列的一个或少数几个碱基的变化引起的。胞嘧啶碱基编辑器(cbe)和腺嘌呤碱基编辑器(abe)分别可以在编辑窗口内实现由胞嘧啶(c)到胸腺嘧啶(t)、腺嘌呤(a)到鸟嘌呤(g)的单个碱基转换。目前，胞嘧啶碱基编辑器和腺嘌呤碱基编辑器已经广泛应用于农作物重要基因功能验证和农作物改良。此外，使用尿嘧啶dna糖基化酶替换尿嘧啶糖基化酶抑制剂基因，还开发出了可以实现胞嘧啶(c)颠换为鸟嘌呤(g)的cgbe技术，并在植物中进行了应用。
3.近期，在腺嘌呤碱基编辑器基础之上，添加人源烷基腺嘌呤dna糖基化酶(hmpg)，开发了适宜动物细胞使用的aybe系统，在动物细胞系中可以实现靶标序列中腺嘌呤(a)到y(y＝胞嘧啶c或胸腺嘧啶t)的单碱基颠换。具体地，将人源烷基腺嘌呤dna糖基化酶(hmpg)与腺嘌呤编辑器融合，在靶点处诱导腺嘌呤脱氨形成肌苷，然后在糖基化酶作用下，移除肌苷产生ap位点(apurinic site)，经基因组错配修复后，可实现腺嘌呤(a)到胞嘧啶(c)或胸腺嘧啶(t)的颠换(图1)。研发植物aybe系统(paybe)将极大地拓展碱基编辑在农作物基因功能分析和遗传改良中的应用范围和潜力。然而，目前在植物中还没有关于植物aybe系统的报道和文献记载。


技术实现要素：

4.本发明要解决的技术问题是如何在植物中实现靶标序列中腺嘌呤(a)到y(y＝胞嘧啶c或胸腺嘧啶t)的单碱基颠换。
5.为了解决上述技术问题，本发明首先提供了成套系统。
6.所述成套系统为如下m1)-m4)中任一种：
7.m1)所述成套系统包括腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和grna；
8.m2)所述成套系统包括腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和grna；
9.m3)所述成套系统包括反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和grna；
10.m4)所述成套系统包括反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和grna；
11.所述人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg为将所述人源烷基腺嘌呤dna糖基化酶hmpg氨基酸序列的第163位由甘氨酸突变为精氨酸，且将第169位由天冬酰胺突变
为丝氨酸，且将第198位由丝氨酸突变为丙氨酸，且将第202位由赖氨酸突变为丙氨酸，且将第203位由甘氨酸突变为丙氨酸，且将第206位由丝氨酸突变为丙氨酸，且将第210位由赖氨酸突变为丙氨酸后得到的蛋白质。
12.上述成套系统中，所述腺嘌呤脱氨酶tada8e为a1)或a2)：
13.a1)氨基酸序列是序列2所示的蛋白质；
14.a2)将序列2所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。
15.所述腺嘌呤脱氨酶tada8e的编码基因为a1)或a2)：
16.a1)序列1第2062-2559位所示的dna分子；
17.a2)与a1)限定的核苷酸序列具有75％或75％以上同一性，且编码所述腺嘌呤脱氨酶tada8e的dna分子。
18.所述cas9(d10a)缺刻酶为b1)或b2)：
19.b1)氨基酸序列是序列3所示的蛋白质；
20.b2)将序列3所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。
21.所述cas9(d10a)缺刻酶的编码基因为b1)或b2)：
22.b1)序列1第2656-6756位所示的dna分子；
23.b2)与b1)限定的核苷酸序列具有75％或75％以上同一性，且编码所述cas9(d10a)缺刻酶的dna分子。
24.所述人源烷基腺嘌呤dna糖基化酶hmpg为c1)或c2)：
25.c1)氨基酸序列是序列6所示的蛋白质；
26.c2)将序列6所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。
27.所述人源烷基腺嘌呤dna糖基化酶hmpg的编码基因为c1)或c2)：
28.c1)序列7所示的dna分子；
29.c2)与c1)限定的核苷酸序列具有75％或75％以上同一性，且编码所述人源烷基腺嘌呤dna糖基化酶hmpg的dna分子。
30.所述反式激活因子vp64为d1)或d2)：
31.d1)氨基酸序列是序列9所示的蛋白质；
32.d2)将序列9所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。
33.所述反式激活因子vp64的编码基因为d1)或d2)：
34.d1)序列8第1-150位所示的dna分子；
35.d2)与d1)限定的核苷酸序列具有75％或75％以上同一性，且编码所述反式激活因子vp64的dna分子。
36.所述grna靶向靶点序列；所述grna依次由所述靶点序列转录的rna和grna骨架组成，所述grna骨架为将序列1第7854-7929位中的t替换为u后得到的rna分子。
37.上述成套系统中，所述m1)的成套系统包括依次由腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶和人源烷基腺嘌呤dna糖基化酶hmpg融合而成的融合蛋白和grna；
38.所述m2)中的成套系统包括依次由腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶和人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg融合而成的融合蛋白和grna；
39.所述m3)中的成套系统包括依次由反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶和人源烷基腺嘌呤dna糖基化酶hmpg融合而成的融合蛋白和grna；
40.所述m4)中的成套系统包括依次由反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶和人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg融合而成的融合蛋白和grna。
41.进一步的，所述m1)的成套系统包括依次由核定位信号、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和核定位信号融合而成的融合蛋白和grna；
42.所述m2)中的成套系统包括依次由核定位信号、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和核定位信号融合而成的融合蛋白和grna；
43.所述m3)中的成套系统包括依次由核定位信号、反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和核定位信号融合而成的融合蛋白和grna；
44.所述m4)中的成套系统包括依次由核定位信号、反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和核定位信号融合而成的融合蛋白和grna。
45.更进一步的，所述成套系统还包括筛选剂抗性蛋白。所述筛选剂抗性蛋白可为本技术领域中公知的各种抗性蛋白，如卡那霉素抗性蛋白、除草剂抗性蛋白、磷酸甘露糖异构酶等。
46.所述m1)的成套系统包括依次由核定位信号、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和核定位信号融合而成的融合蛋白、grna和筛选剂抗性蛋白；
47.所述m2)中的成套系统包括依次由核定位信号、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和核定位信号融合而成的融合蛋白、grna和筛选剂抗性蛋白；
48.所述m3)中的成套系统包括依次由核定位信号、反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和核定位信号融合而成的融合蛋白、grna和筛选剂抗性蛋白；
49.所述m4)中的成套系统包括依次由核定位信号、反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和核定位信号融合而成的融合蛋白、grna和筛选剂抗性蛋白。
50.在本发明的一个具体实施例中，所述筛选剂抗性蛋白为潮霉素磷酸转移酶；所述潮霉素磷酸转移酶为e1)或e2)：
51.e1)氨基酸序列是序列4所示的蛋白质；
52.e2)将序列4所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。
53.所述潮霉素磷酸转移酶的编码基因为e1)或e2)：
54.e1)序列1第9018-10043位所示的dna分子；
55.e2)与e1)限定的核苷酸序列具有75％或75％以上同一性，且编码所述潮霉素磷酸转移酶的dna分子。
56.所述核定位信号为核定位信号nls。所述核定位信号nls为f1)或f2)：
57.f1)氨基酸序列是序列10所示的蛋白质；
58.f2)将序列10所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。
59.所述核定位信号nls的编码基因为f1)或f2)：
60.f1)序列1第2017-2037位所示的dna分子；
61.f2)与f1)限定的核苷酸序列具有75％或75％以上同一性，且编码所述核定位信号nls的dna分子。
62.上述任一所述蛋白质中，所述一个或几个氨基酸残基的取代和/或缺失和/或添加为不超过10个氨基酸残基的取代和/或缺失和/或添加或不超过9个氨基酸残基的取代和/或缺失和/或添加或不超过8个氨基酸残基的取代和/或缺失和/或添加或不超过7个氨基酸残基的取代和/或缺失和/或添加或不超过6个氨基酸残基的取代和/或缺失和/或添加或不超过5个氨基酸残基的取代和/或缺失和/或添加或不超过4个氨基酸残基的取代和/或缺失和/或添加或不超过3个氨基酸残基的取代和/或缺失和/或添加或不超过2个氨基酸残基的取代和/或缺失和/或添加或不超过1个氨基酸残基的取代和/或缺失和/或添加。
63.上述任一所述蛋白质可人工合成，也可先合成其编码基因，再进行生物表达得到。
64.上述任一所述编码基因中的同一性是指与天然核酸序列的序列相似性。所述同一性包括与本发明的编码序列2、序列3、序列4、序列6、序列9所示的氨基酸序列组成的蛋白质的核苷酸序列具有75％或更高，或80％或更高，或85％或更高，或90％或更高，或91％或更高，或92％或更高，或93％或更高，或94％或更高，或95％或更高，或96％或更高，或97％或更高，或98％或更高，或99％或更高同一性的核苷酸序列。同一性可以用肉眼或计算机软件进行评价。使用计算机软件，两个或多个序列之间的同一性可以用百分比(％)表示，其可以用来评价相关序列之间的同一性。
65.上述成套系统中，所述m1)的成套系统包括依次由核定位信号、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和核定位信号融合而成的融合蛋白表达盒、grna表达盒和筛选剂抗性蛋白表达盒；
66.所述m2)中的成套系统包括依次由核定位信号、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和核定位信号融合而成的融合蛋白表达盒、grna表达盒和筛选剂抗性蛋白表达盒；
67.所述m3)中的成套系统包括依次由核定位信号、反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和核定位信号融合而成的融合蛋白表达盒、grna表达盒和筛选剂抗性蛋白表达盒；
68.所述m4)中的成套系统包括依次由核定位信号、反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和核定位信号融合而成的融合蛋白表达盒、grna表达盒和筛选剂抗性蛋白表达盒。
69.进一步的，在所述成套系统的融合蛋白表达盒中，所述融合蛋白由ubi启动子驱动表达。所述ubi启动子的核苷酸序列如序列1第1-1992位所示。
70.在所述成套系统的grna表达盒中，所述grna由osu3启动子驱动表达。所述osu3启动子的核苷酸序列如序列1第7453-7833位所示。
71.在所述成套系统的筛选标记蛋白表达盒中，所述筛选标记蛋白由35s启动子驱动表达。所述35s启动子的核苷酸序列如序列1第8296-8973位所示。
72.更进一步的，所述成套系统由重组载体进行表达。所述成套系统包括的各个元件或各个表达盒可通过同一个载体进行表达，也可以通过多个载体进行表达。
73.在本发明的一个具体实施例中，所述m1)的成套系统通过下文中的重组载体paybev1-hmpg-osdep1、重组载体paybev1-hmpg-osepsps、重组载体paybev1-hmpg-osnrt1.1b、重组载体paybev1-hmpg-oswaxy-t1、重组载体paybev1-hmpg-oswaxy-t2或重组载体paybev1-hmpg-oswaxy-t3进行表达。
74.所述m2)的成套系统通过下文中的重组载体paybev2-mhmpg-osdep1、重组载体paybev2-mhmpg-osepsps、重组载体paybev2-mhmpg-osnrt1.1b、重组载体paybev2-mhmpg-oswaxy-t1、重组载体paybev2-mhmpg-oswaxy-t2或重组载体paybev2-mhmpg-oswaxy-t3进行表达。
75.所述m3)的成套系统通过下文中的重组载体paybev3-hmpg-vp64-osdep1、重组载体paybev3-hmpg-vp64-osepsps、重组载体paybev3-hmpg-vp64-osnrt1.1b、重组载体paybev3-hmpg-vp64-oswaxy-t1、重组载体paybev3-hmpg-vp64-oswaxy-t2或重组载体paybev3-hmpg-vp64-oswaxy-t3进行表达。
76.所述m4)的成套系统通过下文中的重组载体paybev4-mhmpg-vp64-osdep1、重组载体paybev4-mhmpg-vp64-osepsps、重组载体paybev4-mhmpg-vp64-osnrt1.1b、重组载体paybev4-mhmpg-vp64-oswaxy-t1、重组载体paybev4-mhmpg-vp64-oswaxy-t2或重组载体paybev4-mhmpg-vp64-oswaxy-t3进行表达。
77.为了解决上述技术问题，本发明又提供了上述成套系统的新用途。
78.本发明提供了上述成套系统在如下s1)-s9)任一种中的应用：
79.s1)植物基因组序列的编辑；
80.s2)制备植物基因组序列的编辑的产品；
81.s3)提高植物基因组序列的编辑效率；
82.s4)制备提高植物基因组序列的编辑效率的产品；
83.s5)拓展植物基因组序列的编辑窗口；
84.s6)制备拓展植物基因组序列的编辑窗口的产品；
85.s7)制备植物突变体；
86.s8)植物育种；
87.s9)制备植物育种的产品。
88.为了解决上述技术问题，本发明还提供了如下t1)-t4)任一所述的方法：
89.t1)植物基因组序列的编辑方法，包括如下步骤：使植物表达上述成套系统，以实现植物基因组序列的编辑；
90.t2)提高植物基因组序列的编辑效率的方法，包括如下步骤：使植物表达上述成套
系统，以实现提高植物基因组序列的编辑效率；
91.t3)拓展植物基因组序列的编辑窗口的方法，包括如下步骤：使植物表达上述成套系统，以实现拓展植物基因组序列的编辑窗口；
92.t4)植物突变体的制备方法，包括如下步骤：使植物表达上述成套系统，以获得植物突变体。
93.上述任一所述方法中，所述使植物表达上述成套系统的方法为将上述成套系统包括的各个元件或各个表达盒导入植物中。
94.进一步的，所述成套系统包括的各个元件或各个表达盒可通过同一个载体导入植物中，也可以通过多个载体导入植物中。
95.更进一步的，所述成套系统通过同一个重组载体导入植物中。
96.在本发明的具体实施例中，所述成套系统通过下文中的重组载体paybev1-hmpg-osdep1、重组载体paybev1-hmpg-osepsps、重组载体paybev1-hmpg-osnrt1.1b、重组载体paybev1-hmpg-oswaxy-t1、重组载体paybev1-hmpg-oswaxy-t2、重组载体paybev1-hmpg-oswaxy-t3、重组载体paybev2-mhmpg-osdep1、重组载体paybev2-mhmpg-osepsps、重组载体paybev2-mhmpg-osnrt1.1b、重组载体paybev2-mhmpg-oswaxy-t1、重组载体paybev2-mhmpg-oswaxy-t2、重组载体paybev2-mhmpg-oswaxy-t3、重组载体paybev3-hmpg-vp64-osdep1、重组载体paybev3-hmpg-vp64-osepsps、重组载体paybev3-hmpg-vp64-osnrt1.1b、重组载体paybev3-hmpg-vp64-oswaxy-t1、重组载体paybev3-hmpg-vp64-oswaxy-t2、重组载体paybev3-hmpg-vp64-oswaxy-t3、重组载体paybev4-mhmpg-vp64-osdep1、重组载体paybev4-mhmpg-vp64-osepsps、重组载体paybev4-mhmpg-vp64-osnrt1.1b、重组载体paybev4-mhmpg-vp64-oswaxy-t1、重组载体paybev4-mhmpg-vp64-oswaxy-t2或重组载体paybev4-mhmpg-vp64-oswaxy-t3导入植物中。
97.上述任一所述成套系统或应用或方法中，所述植物基因组序列的编辑包括植物基因组序列的碱基替换(如单碱基替换和小片段替换)、碱基插入(如单碱基插入和小片段插入)和碱基缺失(如单碱基缺失和小片段缺失)。所述单碱基替换包括单碱基转换(如碱基a转换为碱基g)和单碱基颠换(如碱基a颠换为碱基c或碱基t)。
98.上述任一所述成套系统或应用或方法中，所述植物为x1)或x2)或x3)：
99.x1)单子叶植物或双子叶植物；
100.x2)禾本科植物；
101.x3)水稻(如中花11)。
102.本发明设计并研发了一系列新型植物单碱基编辑器paybes。首先通过将腺嘌呤碱基编辑器tada8e-ncas9(d10a)与水稻密码子优化的hmpg融合，在植物中建立了第一代a-to-y单碱基编辑系统paybev1，该系统由腺嘌呤脱氨酶tada8e、ncas9(d10a)、人源烷基腺嘌呤dna糖基化酶(hmpg)以及基因组靶向引导rna(grna)组成，腺嘌呤脱氨酶将腺嘌呤(a)催化脱氨产生次黄嘌呤(i)，mpg蛋白可以移除次黄嘌呤形成无碱基位点ap site(apurinic site)并通过dna聚合酶错配修复dna双链，实现碱基a-to-y的碱基颠换。此外，本发明又将hmpg替换为mhmpg突变体建立了第二代a-to-y单碱基编辑系统paybev2，测试了mhmpg介导的a-to-y的单碱基编辑效率。在此基础上，本发明还将反式激活因子vp64融合在第一代paybev1系统的tada8e-ncas9(d10a)-hmpg的n端，建立了第三代a-to-y单碱基编辑系统
paybev3，提高了a-to-y的单碱基编辑效率。进一步，通过vp64耦合mhmpg突变体建立了第四代a-to-y单碱基编辑系统paybev4，进一步提高了a-to-y碱基编辑效率。
103.本发明还利用了水稻6个不同的内源靶点测试了4代不同版本的植物paybe编辑器的a-to-y编辑效率。原生质体测试结果表明，除了paybev1只在一个靶点发生了a-to-c外，paybev2、paybev3、paybev4在6个靶点处均实现了a-to-y碱基颠换。与paybev1相比，paybev2、paybev3、paybev4编辑效率分别得到显著提高，其中，paybev4载体介导的植物a-to-y碱基编辑效率最高。与paybev1、paybev2、paybev3相比，paybev4介导的植物a-to-c碱基编辑效率可分别提高9.78倍、3.06-5.15倍、1.27-5.67倍，最高a-to-c编辑效率可达4.40％。与paybev2和paybev3相比，a-to-t碱基编辑效率分别提高1.97-3.18倍和1.57-6.57倍，最高a-to-t编辑效率可达2.79％(表3-表5)。此外，paybev1、paybev2编辑窗口主要发生在a5-a6位点(图5)，paybev3、paybev4编辑窗口除了主要发生在a5-a6外，在a4和a8位也发生了不同程度的a-to-y编辑(图5)，表明耦合vp64后，可扩展paybe的编辑窗口。另外，paybev1、paybev2、paybev3、paybev4均可产生精准indel事件。此外，和对照pabe8e相比，四个不同版本的paybes均显著提高了的a-to-g编辑效率(图4、表3-表6)。以上结果表明，除了进行a-to-y单碱基编辑外，paybes尤其是paybev4-mhmpg-vp64，可用于农作物重要基因编码区或启动区的定向进化。
104.进一步稳定遗传转化结果表明，paybev1仅在oswaxy-t2靶点处实现了a-to-y碱基颠换，效率为2.02％(表9、表10)。paybev3在oswaxy-t2靶点处实现了a-to-t碱基颠换，效率为2.53％(表9、表10)。同时，在osnrt1.1b、oswaxy-t1、oswaxy-t3等靶点处实现了a-to-t碱基编辑，碱基编辑效率依次是1.32％、1.54％、9.46％。相比paybev1，paybev3在oswaxy-t2靶点的a-to-t效率提高了1.25倍(2.53％/2.02％)，同时在其他三个靶点，首次实现了a-to-t碱基编辑。与原生质体测试结果相似，paybev4在稳定转化中a-to-t编辑效率优于其他三种碱基编码器(图7、图8、图9、图10、表9、表10)。除osepsps靶点外，paybev4在osdep1、osnrt1.1b、oswaxy-t1、oswaxy-t2、oswaxy-t3等靶点处均实现了碱基a-to-t的编辑，碱基编辑效率依次是10.98％(9/82)、2.22％(2/90)、4.23％(3/71)、8.54％(7/82)、17.24％(10/58)。相比paybev3，paybev4在osdep1靶点首次实现a-to-t碱基编辑，在osnrt1.1b靶点的a-to-t效率提高了1.68倍(2.22％/1.32％)，在oswaxy-t1靶点的a-to-t效率提高了2.75倍(4.23％/1.54％)，在oswaxy-t2靶点的a-to-t效率提高了3.38倍(8.54％/2.53％)，在oswaxy-t3靶点的a-to-t效率提高了1.82倍(17.24％/9.46％)。根据幼苗基因型a-to-t的结果，确认a-to-t碱基编辑主要发生在a5和a6位。
105.由上述结果可知，本发明成功地设计了四个用于水稻a-to-y碱基编辑的paybes编辑器，并在水稻原生质体和稳定转化中对不同版本的paybes的编辑活性进行了测试，结果表明，带有hmpg的paybev1具有极低的a-to-y碱基编辑活性，与mhmpg结合的paybev2或与vp64偶联的paybev3均能显著提高a-to-y编辑效率；含有反式激活因子vp64的paybev3诱导a-to-y编辑的效果优于含有mhmpg的paybev2。此外，mhmpg和vp64在增强a-to-y编辑活性方面表现出协同效应。含有mhmpg和vp64的paybev4在a-to-y编辑方面优于其他版本的paybes，在水稻稳定系中实现了一些不可编辑的目标位点的a-to-y碱基颠换，a-to-y编辑效率最高达到17.24％。通过与abes、cbes和cgbes的结合，具有更高的a-to-y颠换效率的paybes将允许所有类型的碱基转换，从而极大地丰富植物单碱基编辑工具箱，为水稻和其
它农作物重要基因功能分析和遗传改良提供了重要技术支撑。
附图说明
106.图1为aybe系统a-to-y单碱基颠换原理示意图。aybe系统由腺嘌呤脱氨酶tada8e、crispr/ncas9(d10a)、烷基腺嘌呤dna糖基化酶(mpg)以及基因组靶向引导rna(grna)组成。腺嘌呤脱氨酶将腺嘌呤(a)催化脱氨产生次黄嘌呤(i)，mpg蛋白可以移除次黄嘌呤形成无碱基位点ap site(apurinic site)并通过dna聚合酶错配修复dna双链，实现碱基a-to-y的碱基颠换；因ncas9(d10a)切割靶标链，同时与非靶标链上ap site有机会形成双链断裂(double strand break，dsb)，在靶标位点形成indel编辑事件。tls(translesion synthesis)表示跨损伤修复与复制。
107.图2为人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和人源烷基腺嘌呤dna糖基化酶hmpg氨基酸序列的序列比对图。其中，红色方框内标注为mhmpg突变的7个氨基酸，分别是g163r、n169s、s198a、k202a、g203a、s206a、k210a。
108.图3为植物paybes载体的结构示意图。其中，tada8e为腺嘌呤脱氨酶；ncas9(d10a)为cas9(d10a)缺刻酶；nls为核定位信号；vp64为反式激活因子；hmpg为人源烷基腺嘌呤dna糖基化酶；mhmpg为人源烷基腺嘌呤dna糖基化酶七突突变体。
109.图4为水稻原生质体中测试paybe单碱基编辑器效率的统计图。图a为pabe8e在水稻原生质体中测试a-to-g编辑效率统计图。图b为paybev1-hmpg在水稻原生质体中测试a-to-c、a-to-t、a-to-g以及靶点处indel编辑效率的统计图。图c为paybev2-mhmpg在水稻原生质体中测试a-to-c、a-to-t、a-to-g以及靶点处indel编辑效率的统计图。图d为paybev3-hmpg-vp64在水稻原生质体中测试a-to-c、a-to-t、a-to-g以及靶点处indel编辑效率的统计图。图e为paybev4-mhmpg-vp64在水稻原生质体中测试a-to-c、a-to-t、a-to-g以及靶点处indel编辑效率的统计图。以上所有载体在测试时，每个靶点均在原生质体中进行了三次独立重复实验。
110.图5为paybe在水稻原生质体中测试paybes靶点编辑窗口的偏好性。其中，靶点序列以1-20位置排列，21-23位为pam序列；箱式图跨越25％-75％四分位数，其中水平线表示中位数，误差线延伸至最小值与最大值。图a为paybev1-hmpg靶点编辑窗口偏好性的统计图。图b为paybev2-mhmpg靶点编辑窗口偏好性的统计图。图c为paybev3-hmpg-vp64靶点编辑窗口偏好性的统计图。图d为paybev4-mhmpg-vp64靶点编辑窗口偏好性的统计图。
111.图6为水稻原生质体中paybev4-mhmpg-vp64产生的indel基因型。其中，绿色标注为pam位点，pam位点下方比值表示indel基因型在所有基因型事件的占比，黑色部分代表删除核酸片段，红色表示插入碱基片段。图a为paybev4-mhmpg-vp64在osnrt1.1b靶点编辑后深度测序所得的indel基因型。图b为paybev4-mhmpg-vp64在oswaxy-t1靶点编辑后深度测序所得的indel基因型。图c为paybev4-mhmpg-vp64在oswaxy-t2靶点编辑后深度测序所得的indel基因型。
112.图7为paybev1-hmpg介导稳定转化幼苗a-to-y编辑的sanger测序结果。其中，箭头表示碱基编辑发生位点。
113.图8为paybev3-hmpg-vp64介导稳定转化幼苗a-to-y编辑的sanger测序结果。其中，箭头表示碱基编辑发生位点。
114.图9为paybev4-mhmpg-vp64介导稳定转化幼苗a-to-y编辑的sanger测序结果。其中，箭头表示碱基编辑发生位点。
115.图10为paybev4-mhmpg-vp64介导indel编辑水稻植株的基因型分析。其中，绿色标注为pam位点，pam位点下方表示删除和插入基因型与所有转化事件的比例，黑色部分代表删除核酸片段，红色表示插入碱基片段，深灰色表示a-to-g编辑事件。图a为paybev4-mhmpg-vp64在osdep1靶点编辑后的indel基因型。图b为paybev4-mhmpg-vp64在osepsps靶点编辑后的indel基因型。图c为paybev4-mhmpg-vp64在osnrt1.1b靶点编辑后的indel基因型。图d为paybev4-mhmpg-vp64在oswaxy-t1靶点编辑后的indel基因型。图e为paybev4-mhmpg-vp64在oswaxy-t2靶点编辑后的indel基因型。图f为paybev4-mhmpg-vp64在oswaxy-t3靶点编辑后的indel基因型。
具体实施方式
116.下面结合具体实施方式对本发明进行进一步的详细描述，给出的实施例仅为了阐明本发明，而不是为了限制本发明的范围。以下提供的实施例可作为本技术领域普通技术人员进行进一步改进的指南，并不以任何方式构成对本发明的限制。
117.下述实施例中的实验方法，如无特殊说明，均为常规方法，按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
118.下述实施例中的1/2ms培养基(ph 5.8)是将ms&vitamins盐、蔗糖、植物凝胶和水混匀得到的培养基，其中各溶质在1/2ms培养基中的浓度分别为：2.15g/l ms&vitamins盐、15g/l蔗糖、2g/l植物凝胶。
119.下述实施例中的0.6m mannitol溶液(ph 5.8)是将mannitol和水混匀得到的溶液，其中mannitol浓度为：2.186g/20ml。
120.下述实施例中的酶解液(ph 5.8)是将cellulase rs、macerozyme r-10、mannitol、ph 5.7的mes、cacl2、bsa和水混匀得到的溶液，其中各溶质在酶解液中的浓度分别为：1.5％cellulase rs、0.75％ macerozyme r-10、0.6m mannitol、ph 5.7的mes、10mm cacl2、0.1％bsa。
121.下述实施例中的w5溶液(ph 5.8)是将nacl、cacl2、kcl、ph 5.7的mes和水混匀得到的溶液，其中各溶质在w5溶液中的浓度分别为：0.9％ nacl、125mm cacl2、5mm kcl、2mm mes。
122.下述实施例中的mmg溶液(ph 5.8)是将mannitol、mgcl2、mes和水混匀得到的溶液，其中各溶质在mmg溶液中的浓度分别为：0.6m mannitol、15mm mgcl2、4mm mes。
123.下述实施例中的peg-4000溶液(ph 5.8)是将peg4000、mannitol、cacl2和水混匀得到的溶液，其中各溶质在peg-4000溶液中的浓度分别为：40％(w/v)peg4000、0.6mmannitol、0.1m cacl2。
124.下述实施例中的lb液体培养基是将酵母粉、胰蛋白胨、nacl和水混匀得到的溶液，其中各溶质在lb液体培养基中的浓度分别为：0.5％(w/v)酵母粉、1％(w/v)胰蛋白胨、1％(w/v)nacl。
125.下述实施例中的aam培养基(ph 5.2)是将ms salts&vitamins盐、蔗糖、mes、葡萄
糖、酪蛋白氨基酸、乙酰丁香酮和100ml 10x aa amino acids混匀得到的培养基，其中各溶质在aam培养基中的浓度分别为4.3g/l ms salts&vitamins盐、68.5g/l蔗糖、0.5g/l mes、36g/l葡萄糖、500mg/l酪蛋白氨基酸、40mg/l乙酰丁香酮。上述10x aa amino acids溶液为将l-谷氨酰胺、l-天(门)冬氨酸、l-精氨酸、甘氨酸和水混匀得到的溶液，其中各溶质在10x aa amino acids溶液中的浓度为：8.76g/l l-谷氨酰胺、2.66g/l l-天(门)冬氨酸、1.74g/l l-精氨酸和75mg/l甘氨酸。
126.下述实施例中的r1培养基(ph 5.8)是将ms&vitamins盐、蔗糖、mes、酪蛋白氨基酸、l-脯氨酸、2，4-d、植物凝胶和水混匀得到的培养基，其中各溶质在r1培养基中的浓度分别为：4.3g/l ms&vitamins盐、30g/l蔗糖、0.5g/l mes、300mg/l酪蛋白氨基酸、2.8g/l l-脯氨酸、2mg/l 2，4-d、4g/l植物凝胶。
127.下述实施例中的r2培养基(ph 5.2)是将ms&vitamins盐、蔗糖、mes、酪蛋白氨基酸、2，4-d、植物凝胶、乙酰丁香酮和水混匀得到的培养基，其中各溶质在r2培养基中的浓度分别为：4.3g/l ms&vitamins盐、30g/l蔗糖、0.5g/l mes、300mg/l酪蛋白氨基酸、2mg/l 2，4-d、4g/l植物凝胶、20mg/ml乙酰丁香酮。
128.下述实施例中的r1筛选培养基(ph 5.8)是将ms&vitamins盐、蔗糖、mes、酪蛋白氨基酸、l-脯氨酸、2，4-d、植物凝胶和水混匀得到的培养基，其中各溶质在r1筛选培养基中的浓度为：4.3g/l ms&vitamins盐、30g/l蔗糖、0.5g/l mes、300mg/l酪蛋白氨基酸、2.8g/l l-脯氨酸、2mg/l 2，4-d、4g/l植物凝胶。
129.下述实施例中的r4分化培养基(ph 5.8)是将ms&vitamins盐、蔗糖、mes、酪蛋白氨基酸、山梨醇、激动素、naa、植物凝胶和水混匀得到的培养基，其中各溶质在r4分化培养基中的浓度分别为：4.3g/l ms&vitamins盐、30g/l蔗糖、0.5g/l mes、2g/l酪蛋白氨基酸、30g/l山梨醇、2mg/l激动素、1mg/l naa、4g/l植物凝胶。
130.下述实施例中的r5生根培养基(ph 5.8)是将ms&vitamins盐、蔗糖、mes、植物凝胶和水混匀得到的培养基，其中各溶质在r5培养基中的浓度分别为：2.15g/l ms&vitamins盐、15g/l蔗糖、0.5g/l mes、2g/l植物凝胶。
131.实施例1、不同植物a-to-y单碱基编辑器paybes的主要元件设计及其表达载体设计一、不同植物a-to-y单碱基编辑器paybes的主要元件设计
132.本发明中的植物a-to-y单碱基编辑器paybes有如下五种：pabe8e、paybev1-hmpg、paybev2-mhmpg、paybev3-hmpg-vp64和paybev4-mhmpg-vp64。各个植物a-to-y单碱基编辑器的主要元件分别如下：
133.单碱基编辑器pabe8e的主要元件如下：腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、grna和筛选标记蛋白(如hpt)。
134.单碱基编辑器paybev1-hmpg的主要元件如下：腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg、grna和筛选标记蛋白(如hpt)。
135.单碱基编辑器paybev2-mhmpg的主要元件如下：腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg、grna和筛选标记蛋白(如hpt)。
136.单碱基编辑器paybev3-hmpg-vp64的主要元件如下：反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg、grna和筛选标记蛋白(如hpt)。
137.单碱基编辑器paybev4-mhmpg-vp64的主要元件如下：反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg、grna和筛选标记蛋白(如hpt)。
138.上述人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg为将人源烷基腺嘌呤dna糖基化酶hmpg氨基酸序列的第163位由甘氨酸(g)突变为精氨酸(r)，且将第169位由天冬酰胺(n)突变为丝氨酸(s)，且将第198位由丝氨酸(s)突变为丙氨酸(a)，且将第202位由赖氨酸(k)突变为丙氨酸(a)，且将第203位由甘氨酸(g)突变为丙氨酸(a)，且将第206位由丝氨酸(s)突变为丙氨酸(a)，且将第210位由赖氨酸(k)突变为丙氨酸(a)后得到的蛋白质。
139.人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和人源烷基腺嘌呤dna糖基化酶hmpg的氨基酸序列比对图如图2所示。
140.二、不同植物a-to-y单碱基编辑器paybes的表达载体设计
141.上述各植物a-to-y单碱基编辑器paybes表达载体的结构示意图如图3所示。
142.单碱基编辑器pabe8e的表达载体包括融合蛋白tada8e&cas9(d10a)的表达盒、grna表达盒和筛选标记蛋白hpt表达盒。其中，融合蛋白tada8e&cas9(d10a)的表达盒依次包括ubi启动子、核定位信号nls、腺嘌呤脱氨酶tada8e的编码基因、cas9(d10a)缺刻酶的编码基因、核定位信号nls和e9t终止子。
143.单碱基编辑器paybev1-hmpg的表达载体包括融合蛋白tada8e&cas9(d10a)&hmpg的表达盒、grna表达盒和筛选标记蛋白hpt表达盒。其中，融合蛋白tada8e&cas9(d10a)&hmpg的表达盒依次包括ubi启动子、核定位信号nls、腺嘌呤脱氨酶tada8e的编码基因、cas9(d10a)缺刻酶的编码基因、人源烷基腺嘌呤dna糖基化酶hmpg的编码基因、核定位信号nls和e9t终止子。
144.单碱基编辑器paybev2-mhmpg的表达载体包括融合蛋白tada8e&cas9(d10a)&mhmpg的表达盒、grna表达盒和筛选标记蛋白hpt表达盒。其中，融合蛋白tada8e&cas9(d10a)&mhmpg的表达盒依次包括ubi启动子、核定位信号nls、腺嘌呤脱氨酶tada8e的编码基因、cas9(d10a)缺刻酶的编码基因、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg的编码基因、核定位信号nls和e9t终止子。
145.单碱基编辑器paybev3-hmpg-vp64的表达载体包括融合蛋白vp64&tada8e&cas9(d10a)&hmpg的表达盒、grna表达盒和筛选标记蛋白hpt表达盒。其中，融合蛋白vp64&tada8e&cas9(d10a)&hmpg的表达盒依次包括ubi启动子、核定位信号nls、反式激活因子vp64的编码基因、腺嘌呤脱氨酶tada8e的编码基因、cas9(d10a)缺刻酶的编码基因、人源烷基腺嘌呤dna糖基化酶hmpg的编码基因、核定位信号nls和e9t终止子。
146.单碱基编辑器paybev4-mhmpg-vp64的表达载体包括融合蛋白vp64&tada8e&cas9(d10a)&mhmpg的表达盒、grna表达盒和筛选标记蛋白hpt表达盒。其中，融合蛋白vp64&tada8e&cas9(d10a)&hmpg的表达盒依次包括ubi启动子、核定位信号nls、反式激活因子vp64的编码基因、腺嘌呤脱氨酶tada8e的编码基因、cas9(d10a)缺刻酶的编码基因、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg的编码基因、核定位信号nls和e9t终止子。
147.上述各单碱基编辑器中的grna表达盒均依次包括osu3启动子、grna的编码基因和sup4终止子。
148.上述各单碱基编辑器中的筛选标记蛋白hpt表达盒均依次包括35s启动子、筛选标
记蛋白hpt的编码基因和nos终止子。
149.实施例2、不同植物a-to-y单碱基编辑器paybes表达载体的构建及其在对水稻基因组进行碱基编辑中的应用
150.一、不同植物a-to-y单碱基编辑器paybes表达载体的构建
151.人工构建如下重组载体，各载体均为环状质粒：
152.单碱基编辑器pabe8e的表达载体共计6个，分别是pabe8e-osdep1、pabe8e-osepsps、pabe8e-osnrt1.1b、pabe8e-oswaxy-t1、pabe8e-oswaxy-t2、pabe8e-oswaxy-t3重组载体。
153.单碱基编辑器paybev1-hmpg的表达载体共计6个，分别是paybev1-hmpg-osdep1、paybev1-hmpg-osepsps、paybev1-hmpg-osnrt1.1b、paybev1-hmpg-oswaxy-t1、paybev1-hmpg-oswaxy-t2、paybev1-hmpg-oswaxy-t3重组载体。
154.单碱基编辑器paybev2-mhmpg的表达载体共计6个，分别是paybev2-mhmpg-osdep1、paybev2-mhmpg-osepsps、paybev2-mhmpg-osnrt1.1b、paybev2-mhmpg-oswaxy-t1、paybev2-mhmpg-oswaxy-t2、paybev2-mhmpg-oswaxy-t3重组载体。
155.单碱基编辑器paybev3-hmpg-vp64的表达载体共计6个，分别是paybev3-hmpg-vp64-osdep1、paybev3-hmpg-vp64-osepsps、paybev3-hmpg-vp64-osnrt1.1b、paybev3-hmpg-vp64-oswaxy-t1、paybev3-hmpg-vp64-oswaxy-t2、paybev3-hmpg-vp64-oswaxy-t3重组载体。
156.单碱基编辑器paybev4-mhmpg-vp64的表达载体共计6个，分别是paybev4-mhmpg-vp64-osdep1、paybev4-mhmpg-vp64-osepsps、paybev4-mhmpg-vp64-osnrt1.1b、paybev4-mhmpg-vp64-oswaxy-t1、paybev4-mhmpg-vp64-oswaxy-t2、paybev4-mhmpg-vp64-oswaxy-t3重组载体。
157.pabe8e-osdep1重组载体的核苷酸序列如序列表中序列1所示，其中，第1-1992位为ubi启动子序列，第2017-2037位为核定位信号nls的编码序列，编码序列10所示的核定位信号nls，第2062-2559位为腺嘌呤脱氨酶tada8e的编码基因序列，编码序列2所示的腺嘌呤脱氨酶tada8e，第2656-6756位为cas9(d10a)缺刻酶的编码基因序列，编码序列3所示的cas9(d10a)缺刻酶，第6757-6804位为核定位信号nls的编码序列，第6812-7446位为e9t终止子序列；第7453-7833位为osu3启动子序列、第7834-7853位为grna靶点序列，第7854-7929位为grna骨架的编码序列，第7930-7949位为sup4终止子序列；第8296-8973位为35s启动子序列，第9018-10043位为筛选标记蛋白hpt的编码基因序列，编码序列4所示的筛选标记蛋白hpt，第10090-10264位为nos终止子序列。
158.pabe8e-osepsps、pabe8e-osnrt1.1b、pabe8e-oswaxy-t1、pabe8e-oswaxy-t2、pabe8e-oswaxy-t3重组载体的核苷酸序列为将pabe8e-osdep1重组载体序列中的grna靶点序列分别替换为表1中osepsps、osnrt1.1b、oswaxy-t1、oswaxy-t2、oswaxy-t3对应的靶点序列，且保持其它序列不变后得到的序列。
159.paybev1-hmpg-osdep1、paybev1-hmpg-osepsps、paybev1-hmpg-osnrt1.1b、paybev1-hmpg-oswaxy-t1、paybev1-hmpg-oswaxy-t2、paybev1-hmpg-oswaxy-t3重组载体的核苷酸序列分别为将pabe8e-osdep1、pabe8e-osepsps、pabe8e-osnrt1.1b、pabe8e-oswaxy-t1、pabe8e-oswaxy-t2、pabe8e-oswaxy-t3重组载体序列中的第6756位和第6757位
之间插入序列5所示的dna分子，且保持其它序列不变后得到的序列。其中，序列5第40-930位为人源烷基腺嘌呤dna糖基化酶hmpg的编码序列，编码序列6所示的人源烷基腺嘌呤dna糖基化酶hmpg。
160.paybev2-mhmpg-osdep1、paybev2-mhmpg-osepsps、paybev2-mhmpg-osnrt1.1b、paybev2-mhmpg-oswaxy-t1、paybev2-mhmpg-oswaxy-t2、paybev2-mhmpg-oswaxy-t3重组载体的核苷酸序列分别为将paybev1-hmpg-osdep1、paybev1-hmpg-osepsps、paybev1-hmpg-osnrt1.1b、paybev1-hmpg-oswaxy-t1、paybev1-hmpg-oswaxy-t2、paybev1-hmpg-oswaxy-t3重组载体序列中的人源烷基腺嘌呤dna糖基化酶hmpg的编码序列替换为序列7所示的人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg的编码序列，且保持其它序列不变后得到的序列。
161.paybev3-hmpg-vp64-osdep1、paybev3-hmpg-vp64-osepsps、paybev3-hmpg-vp64-osnrt1.1b、paybev3-hmpg-vp64-oswaxy-t1、paybev3-hmpg-vp64-oswaxy-t2、paybev3-hmpg-vp64-oswaxy-t3重组载体分别为将paybev1-hmpg-osdep1、paybev1-hmpg-osepsps、paybev1-hmpg-osnrt1.1b、paybev1-hmpg-oswaxy-t1、paybev1-hmpg-oswaxy-t2、paybev1-hmpg-oswaxy-t3重组载体序列中的第2038位和第2061位之间的dna分子替换为序列8所示的dna分子，且保持其它序列不变后得到的序列。其中，序列8第1-150位为反式激活因子vp64的编码序列，编码序列9所示的反式激活因子vp64。
162.paybev4-mhmpg-vp64-osdep1、paybev4-mhmpg-vp64-osepsps、paybev4-mhmpg-vp64-osnrt1.1b、paybev4-mhmpg-vp64-oswaxy-t1、paybev4-mhmpg-vp64-oswaxy-t2、paybev4-mhmpg-vp64-oswaxy-t3重组载体的核苷酸序列分别为将paybev2-mhmpg-osdep1、paybev2-mhmpg-osepsps、paybev2-mhmpg-osnrt1.1b、paybev2-mhmpg-oswaxy-t1、paybev2-mhmpg-oswaxy-t2、paybev2-mhmpg-oswaxy-t3重组载体序列中的第2038位和第2061位之间的dna分子替换为序列8所示的dna分子，且保持其它序列不变后得到的序列。其中，序列8第1-150位为反式激活因子vp64的编码序列，编码序列9所示的反式激活因子vp64。
163.表1、靶点序列
164.基因靶点序列(5
′
to 3
′
)oswaxy-t1cctgacactggagttgattacaaoswaxy-t2cgccaagtacgacgcaaccacggoswaxy-t3ccatacttcaaaggaacttatggosnrt1.1bactagatatctaaaccattaaggosdep1agacaagcttggccctctttgggosepspsatgatatcctcctacatgtcagg
165.二、水稻原生质体的分离、转化以及精准编辑效率、编辑窗口的分析
166.将步骤一中构建的各个植物a-to-y单碱基编辑器paybes的重组表达载体(所有载体在测试时，每个靶点均在原生质体中进行了三次独立重复实验)分别按照如下步骤进行操作：
167.1、水稻原生质体的分离
168.将脱壳的水稻品种中花11种子先用75％的乙醇漂洗10分钟，再用20％的次氯酸钠
处理20分钟，无菌水洗涤5次以上，得到消毒后种子。然后将消毒后种子放在1/2ms培养基上培养2周左右，在26℃，12h光照(150μmol
·
m-2
·
s-1
)条件下培养，得到水稻幼苗。每个组培瓶可放30粒种子，60～90株幼苗可以做一次实验，分离出的原生质体量可以转化大约6个载体，每个质粒的转化设置3个重复。
169.选取水稻幼苗的茎干和叶鞘部分分离原生质体，具体步骤如下：用锋利的刀片切成大约0.5mm宽的细丝，可以20～30个放在一起切。将细丝放入0.6m mannitol溶液中，锡箔纸包裹避光处理10分钟。用75μm尼龙布过滤，将细丝放入50ml酶解液中，避光，真空泵抽真空30min，压强约50kpa，取出后放在室温摇床上消化5～6小时，转速10～20rpm。加等体积w5溶液稀释酶解产物，终止酶解反应后用尼龙布过滤酶解液于50ml离心管收集原生质体，28℃，150g水平离心，升降速调至3，离心5min，弃上清；用10ml w5溶液轻轻悬起，28℃，150g水平离心，升降速为3，离心5min，弃上清；加适量mmg溶液重悬，调整原生质体浓度为2
×
106/ml，血球细胞计数器计数后备用。
170.2、水稻原生质体的转化
171.加入约20μg无内毒素的高浓度aybe测试质粒(步骤一中构建的各个植物a-to-y单碱基编辑器paybes的重组表达载体)于2ml离心管，用剪刀剪去枪尖的枪头吸取200μl步骤1获得的原生质体沿管壁缓慢加入离心管中，轻轻吸打混匀，再加入约220μl peg-4000溶液，轻轻颠倒混匀，避光诱导转化10～20min；加入800μl w5溶液(室温)颠倒混匀，28℃，150g水平离心，升降速为3，离心5min，弃上清；加入1ml w5溶液，颠倒混匀，置28℃暗处培养。
172.转化后48h后，收集原生质体，弃上清。用枪吸打悬起原生质体于2ml离心管中，12000rpm离心1min后弃上清；提取原生质体基因组，具体方法同植物基因组提取方法，pcr扩增所选靶基因的目的片段进行hi-tom高通量测序，hi-tom测序使用引物如表2所示。
173.pcr扩增体系：10μl 2
×
taq酶，1μl正向引物，1μl反向引物，1μl基因组dna模板(60ng/μl)，7μl ddh2o。
174.pcr扩增程序：98℃，3min；98℃，15s；60℃，15s；72℃，15s；35个循环；72℃，5min；16℃保持。
175.表2、深度测序所用引物列表
[0176][0177][0178]
3、paybes在原生质体中编辑效率的统计分析结果
[0179]
本发明在水稻基因组中选定了6个内源靶点，以快速测试各个单碱基编辑器表达载体(下文简称载体)在水稻原生质体中的有效性。具体检测结果如图4、表3-表7所示。
[0180]
pabe8e载体的深度测序与预期结果一致，在6个靶点处均没有检测到有a-to-y碱基编辑事件，只检测到a-to-g碱基转换。其中，在osdep1靶点处，a-to-g碱基编辑效率最高为6.06％，最低为5.49％；在osepsps靶点处，a-to-g碱基编辑效率最高为3.29％，最低为1.46％；在osnrt1.1b靶点处，a-to-g碱基编辑最为有效，效率最高为6.56％，最低为5.71％；在oswaxy-t1靶点处，a-to-g碱基编辑效率最高为3.77％，最低为2.58％；在oswaxy-t2靶点处，a-to-g碱基编辑效率最高为1.34％，最低为1.28％；在oswaxy-t3靶点处，a-to-g碱基编辑效率最高为5.35％，最低为3.57％。
[0181]
paybev1-hmpg载体在osdep1靶点处没有检测到a-to-y碱基编辑事件的发生，a-to-g碱基编辑效率最高为4.05％，最低为3.64％；在osepsps靶点处没有检测到a-to-y碱基编辑事件的发生，a-to-g碱基编辑效率最高为3.28％，最低为1.92％；在osnrt1.1b靶点处深度测序首次在原生质体中检测到了a-to-c编辑事件，平均效率为0.44％，编辑效率最高为0.45％，最低为0.43％。a-to-g碱基编辑效率最高为6.65％，最低为5.80％；在oswaxy-t1靶点处没有检测到a-to-y碱基编辑事件的发生，a-to-g碱基编辑效率最高为2.98％，最低为2.31％；在oswaxy-t2靶点处没有检测到a-to-y碱基编辑事件的发生，a-to-g碱基编辑效率最高为1.67％，最低为1.09％；在oswaxy-t3靶点处没有检测到a-to-y碱基编辑事件的发生，a-to-g碱基编辑效率最高为4.80％，最低为2.38％。综上所述，植物a-to-y单碱基编辑
载体paybev1-hmpg在6个靶点处深度测序的结果显示，在osnrt1.1b靶点处检测到了a-to-c编辑事件，在其他靶点处没有检测到有a-to-y碱基编辑事件的出现。
[0182]
paybev2-mhmpg载体在6个靶点处均检测到了不同程度的a-to-y/a-to-g及靶点附近的indel。其中，在osdep1靶点处a-to-c碱基编辑效率最高为0.72％，最低为0.53％，平均效率为0.63％，靶点处a-to-t碱基编辑效率最高为0.90％，最低为0.49％，平均效率为0.63％，靶点处a-to-y碱基编辑效率最高为1.54％，最低为1.02％，平均效率为1.26％，靶点处indel事件效率最高为0.41％，最低为0.34％，靶点处a-to-g碱基编辑效率最高为25.50％，最低为14.31％；在osepsps检测到靶点处indel事件最高效率为1.95％，靶点处a-to-g碱基编辑效率最高为3.71％，最低为3.29％；在osnrt1.1b靶点处，a-to-c碱基编辑效率最高为1.44％，最低为0.62％，平均效率为0.94％，相比paybev1-hmpg载体，paybev2-mhmpg载体的a-to-c碱基编辑效率提高约3.20倍(1.44％/0.45％)，靶点处indel事件效率最高为0.91％，最低为0.78％，靶点处a-to-g碱基编辑效率最高为9.34％，最低为9.08％；在oswaxy-t1靶点处a-to-t碱基编辑效率最高为0.62％，最低为0.35％，平均效率为0.48％，靶点处a-to-g碱基编辑效率最高为4.49％，最低为4.15％；在oswaxy-t2仅一个重复检测到靶点处a-to-t碱基编辑效率为0.61％，靶点处indel事件效率最高为2.53％，最低为1.98％，靶点处a-to-g碱基编辑效率最高为4.76％，最低为3.79％；在oswaxy-t3靶点处a-to-t碱基编辑效率最高为1.31％，最低为0.62％，平均效率为0.89％，靶点处a-to-g碱基编辑效率最高为6.12％，最低为4.18％。综上所述，paybev2-mhmpg相比paybev1-hmpg，在6个靶点处，实现了a-to-y碱基颠换效率的提高，尤其是osnrt1.1b靶点处，a-to-c碱基编辑效率提高约3.2倍，并且在靶点处实现了小片段的碱基插入和缺失。
[0183]
相比paybev1-hmpg载体，paybev3-hmpg-vp64载体在6个靶点处均检测到了不同程度的a-to-y/a-to-g，与paybev2-mhmpg载体相比，a-to-y碱基编辑效率有不同程度的提升。其中，在osdep1靶点处a-to-c碱基编辑效率最高为2.92％，最低为2.27％，平均效率为2.49％，相比paybev2-mhmpg载体，a-to-c碱基编辑效率提高约4.06倍(2.92％/0.72％)，靶点处a-to-t碱基编辑效率最高为0.75％，最低为0.60％，平均效率为0.66％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率略微有所下降(0.75％/0.9％)，靶点处a-to-y碱基编辑效率最高为3.52％，最低为2.89％，平均效率为3.15％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约2.29倍(3.52％/1.54％)，靶点处a-to-g碱基编辑效率最高为10.59％，最低为6.40％；在osepsps靶点处a-to-c碱基编辑效率最高为0.54％，最低为0.37％，平均效率为0.44％，仅一个重复检测到靶点处a-to-t碱基编辑效率为0.31％，靶点处a-to-y碱基编辑效率最高为0.71％，最低为0.37％，平均效率为0.54％，靶点处a-to-g碱基编辑效率最高为1.01％，最低为0.67％；在osnrt1.1b靶点处a-to-c碱基编辑效率最高为2.06％，最低为1.93％，平均效率为1.99％，相比paybev1-hmpg载体，a-to-c碱基编辑效率提高约4.58倍(2.06％/0.45％)，相比paybev2-mhmpg载体，a-to-c碱基编辑效率提高约1.43倍(2.06％/1.44％)，靶点处a-to-t碱基编辑效率最高为0.48％，最低为0.45％，平均效率为0.47％，靶点处a-to-y碱基编辑效率最高为2.51％，最低为2.40％，平均效率为2.46％，相比paybev1-hmpg载体，a-to-y碱基编辑效率约提高约5.58倍(2.51％/0.45％)，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约1.74倍(2.51％/1.44％)，靶点处a-to-g碱基编辑效率最高为5.57％，最低为4.29％；在oswaxy-t1靶点处a-to-c碱基编辑效率最高为1.44％，最低为
0.34％，平均效率为0.78％，靶点处a-to-t碱基编辑效率最高为0.30％，最低为0.24％，平均效率为0.27％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率略微有所下降(0.30％/0.62％)，靶点处a-to-y碱基编辑效率最高为1.74％，最低为0.34％，平均效率为0.96％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约2.81倍(1.74％/0.62％)，靶点处a-to-g碱基编辑效率最高为3.18％，最低为1.81％；在oswaxy-t2靶点处a-to-c碱基编辑效率最高为1.53％，最低为0.95％，平均效率为1.30％，靶点处a-to-t碱基编辑效率最高为1.40％，最低为0.48％，平均效率为0.82％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率提高约2.30倍(1.40％/0.61％)，靶点处a-to-y碱基编辑效率最高为2.83％，最低为1.43％，平均效率为2.12％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约4.64倍(2.83％/0.61％)，靶点处a-to-g碱基编辑效率最高为5.03％，最低为3.68％；在oswaxy-t3靶点处a-to-c碱基编辑效率最高为0.75％，最低为0.65％，平均效率为0.69％，靶点处a-to-t碱基编辑效率最高为0.49％，最低为0.38％，平均效率为0.44％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率略微有所下降(0.49％/1.31％)，靶点处a-to-y碱基编辑效率最高为1.21％，最低为1.04％，平均效率为1.13％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率略微有所下降(1.24％/1.31％)，靶点处a-to-g碱基编辑效率最高为2.92％，最低为2.28％。综上所述，paybev3-hmpg-vp64载体相比paybev1-hmpg载体与paybev2-mhmpg载体，显著提高了a-to-y碱基颠换的效率，并且在靶点处出现了小片段的碱基插入和缺失。
[0184]
相比paybev1-hmpg载体、paybev2-mhmpg载体与paybev3-hmpg-vp64载体，paybev4-mhmpg-vp64载体在6个靶点处均显著提高了a-to-y的编辑效率。其中，在osdep1靶点处a-to-c碱基编辑效率最高为3.71％，最低为2.73％，平均效率为3.28％，相比paybev2-mhmpg载体，a-to-c碱基编辑效率提高约5.15倍(3.71％/0.72％)，相比paybev3-hmpg-vp64载体，a-to-c碱基编辑效率提高约1.27倍(3.71％/2.92％)，a-to-t碱基编辑效率最高为1.77％，最低为1.12％，平均效率为1.44％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率提高约1.97倍(1.77％/0.90％)，相比paybev3-hmpg-vp64载体，a-to-t碱基编辑效率提高约2.36倍(1.77％/0.75％)，靶点处a-to-y碱基编辑效率最高为5.13％，最低为4.50％，平均效率为4.72％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约3.33倍(5.13％/1.54％)，相比paybev3-hmpg-vp64载体，a-to-y碱基编辑效率提高约1.46倍(5.13％/3.52％)，a-to-g碱基编辑效率最高为9.61％，最低为6.08％；在osepsps靶点处a-to-c碱基编辑效率最高为3.06％，最低为1.58％，平均效率为2.26％，相比paybev3-hmpg-vp64载体，a-to-c碱基编辑效率提高约5.67倍(3.06％/0.54％)，a-to-t碱基编辑效率最高为1.37％，最低为0.64％，平均效率为1.09％，相比paybev3-hmpg-vp64载体，a-to-t碱基编辑效率提高约4.42倍(1.37％/0.31％)，靶点处a-to-y碱基编辑效率最高为3.70％，最低为2.84％，平均效率为3.35％，相比paybev3-hmpg-vp64载体，a-to-y碱基编辑效率提高约5.21倍(3.70％/0.71％)，a-to-g碱基编辑效率最高为7.00％，最低为5.01％；在osnrt1.1b靶点处a-to-c碱基编辑效率最高为4.40％，最低为3.09％，平均效率为3.66％，相比paybev1-hmpg载体，a-to-c碱基编辑效率提高约9.78倍(4.40％/0.45％)，相比paybev2-mhmpg载体，a-to-c碱基编辑效率提高约3.06倍(4.40％/1.44％)，相比paybev3-hmpg-vp64载体，a-to-c碱基编辑效率约提高约2.14倍(4.40％/2.06％)，a-to-t碱基编辑效率最高为2.42％，最低为0.82％，平均效率为1.47％，相比paybev3-hmpg-vp64载体，a-to-t碱基编辑效率提高约
5.04倍(2.42％/0.48％)，靶点处a-to-y碱基编辑效率最高为6.82％，最低为4.25％，平均效率为5.13％，相比paybev1-mhmpg载体，a-to-y碱基编辑效率提高约15.16倍(6.82％/0.45％)，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约4.74倍(6.82％/1.44％)，相比paybev3-hmpg-vp64载体，a-to-y碱基编辑效率提高约2.72倍(6.82％/2.51％)，靶点处indel事件效率最高为1.09％，最低为0.47％，a-to-g碱基编辑效率最高为15.27％，最低为13.08％；在oswaxy-t1靶点处a-to-c碱基编辑效率最高为3.22％，最低为2.43％，平均效率为2.84％，相比paybev3-hmpg-vp64载体，a-to-c碱基编辑效率约提高约2.24倍(3.22％/1.44％)，靶点处a-to-t碱基编辑效率最高为1.97％，最低为1.29％，平均效率为1.53％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率提高约3.18倍(1.97％/0.62％)，相比paybev3-hmpg-vp64载体，a-to-t碱基编辑效率提高约6.57倍(1.97％/0.30％)，靶点处a-to-y碱基编辑效率最高为4.55％，最低为4.17％，平均效率为4.37％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约7.34倍(4.55％/0.62％)，相比paybev3-hmpg-vp64载体，a-to-y碱基编辑效率提高约2.61倍(4.55％/1.74％)，靶点处indel事件效率最高为1.10％，最低为0.62％，靶点处a-to-g碱基编辑效率最高为8.18％，最低为6.15％；在oswaxy-t2靶点处a-to-c碱基编辑效率最高为3.69％，最低为2.30％，平均效率为3.00％，相比paybev3-hmpg-vp64载体，a-to-c碱基编辑效率约提高约2.41倍(3.69％/1.53％)，靶点处a-to-t碱基编辑效率最高为2.20％，最低为1.11％，平均效率为1.58％，相比paybev2-hmpg-vp64载体，a-to-t碱基编辑效率提高约3.61倍(2.20％/0.61％)，相比paybev3-hmpg-vp64载体，a-to-t碱基编辑效率提高约1.57倍(2.20％/1.40％)，靶点处a-to-y碱基编辑效率最高为5.89％，最低为3.41％，平均效率为4.58％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约9.66倍(5.89％/0.61％)；相比paybev3-hmpg-vp64载体，a-to-y碱基编辑效率提高约2.08倍(5.89％/2.83％)，靶点处indel事件效率最高为1.31％，最低为0.66％，靶点处a-to-g碱基编辑效率最高为12.34％，最低为7.74％；在oswaxy-t3靶点处a-to-c碱基编辑效率最高为1.41％，最低为1.21％，平均效率为1.34％，相比paybev3-hmpg-vp64载体，a-to-c碱基编辑效率提高约1.88倍(1.41％/0.75％)，靶点处a-to-t碱基编辑效率最高为2.79％，最低为1.84％，平均效率为2.25％，相比paybev2-mhmpg载体，a-to-t碱基编辑效率提高约2.13倍(2.79％/1.31％)，相比paybev3-hmpg-vp64载体，a-to-t碱基编辑效率提高约5.69倍(2.79％/0.49％)，靶点处a-to-y碱基编辑效率最高为4.20％，最低为3.05％，平均效率为3.59％，相比paybev2-mhmpg载体，a-to-y碱基编辑效率提高约3.21倍(4.20％/1.31％)，相比paybev3-hmpg-vp64载体，a-to-y碱基编辑效率提高约3.39倍(4.20％/1.24％)，靶点处a-to-g碱基编辑效率最高为5.09％，最低为3.90％。综上所述，与paybev2-mhmpg载体相比，paybev4-mhmpg-vp64载体的a-to-c编辑效率提高约3.06-5.15倍，a-to-t编辑效率提高约1.97-3.18倍，a-to-y编辑效率提高约3.21-7.43倍；相比paybev3-hmpg-vp64载体，paybev4-mhmpg-vp64载体的a-to-c编辑效率提高约1.27-5.67倍，a-to-t编辑效率提高约1.57-6.57倍，a-to-y编辑效率提高约1.46-5.21倍。
[0185]
表3、水稻原生质体中paybes在不同靶点a-to-c碱基编辑效率
[0186][0187]
表4、水稻原生质体中paybes在不同靶点a-to-t碱基编辑效率
[0188][0189]
表5、水稻原生质体中paybes在不同靶点a-to-y碱基编辑效率
[0190]
[0191][0192]
表6、水稻原生质体中paybes在不同靶点a-to-g碱基编辑效率
[0193][0194]
表7、水稻原生质体中paybes在不同靶点indel碱基编辑效率
[0195]
[0196][0197]
3、paybes在原生质体中编辑窗口的统计分析结果
[0198]
在原生质体测试结果中，paybes载体在osdep1、osepsps、osnrt1.1b、oswaxy-t1、oswaxy-t2、oswaxy-t3等靶点处均实现了碱基a-to-y的颠换。paybev2-mhmpg载体与paybev3-hmpg-vp64载体介导的碱基a-to-y的颠换主要发生在a5、a6位(图5)，其中，paybev3-hmpg-vp64载体在a5、a6位的a-to-c编辑效率大约是paybev2-mhmpg载体的2倍，a-to-t编辑效率没有明显差别。相比paybev2-mhmpg载体与paybev3-hmpg-vp64载体，paybev4-mhmpg-vp64载体在a5、a6位的a-to-y的颠换编辑效率显著提高，同时拓展了编辑窗口，表现为在a4、a8位的a-to-y编辑效率也有显著提高(图5)。
[0199]
4、paybes在原生质体中indel事件基因型的统计分析结果
[0200]
crispr/ncas9(d10a)在grna引导下只切割靶标链，从而导致dna在靶标链处造成缺口。植物paybe的腺嘌呤脱氨酶将腺嘌呤(a)催化脱氨产生次黄嘌呤(i)，mpg蛋白可以移除次黄嘌呤形成无碱基位点ap site(apurinic site)并通过dna聚合酶错配修复dna双链，实现碱基a-to-y的碱基颠换；同时，因ncas9(d10a)切割靶标链，同时与非靶标链上的ap site有机会形成交错的双链断裂(double strand break，dsb)，通过跨损伤修复与复制在靶标位点产生精准indel编辑事件(图1)。
[0201]
在原生质体测试结果中，paybes载体在osdep1、osepsps、osnrt1.1b、oswaxy-t1、oswaxy-t2等靶点均检测到小片段精准缺失事件。以paybev4-mhmpg-vp64载体为例(图6)，在osdep1靶点中，缺失片段从7bp到10bp不等，其中a8至a14位点之间片段缺失效率为
0.81％，a8位点至ncas9(d10a)切割位点之间缺失片段事件出现效率为0.47％。在osepsps靶点中，缺失片段从6bp到13bp不等，其中a5位点至ncas9(d10a)切割位点之间缺失片段事件出现效率为0.62％，a6至c14位点之间片段缺失效率为0.97％，a9至c14位点之间片段缺失效率为1.10％。在osnrt1.1b靶点中，缺失片段从4bp到12bp，其中g7位点至ncas9(d10a)切割位点之间缺失片段事件出现效率为0.87％，a9至a12位点之间片段缺失效率为0.92％，g7至c18位点之间片段缺失效率为0.66％，同时，在a16与a17位点之间，有着9bp的小片段插入，效率为0.44％。
[0202]
三、农杆菌介导的水稻稳定遗传转化及基因型与表型分析
[0203]
1、农杆菌介导的水稻稳定遗传转化
[0204]
选取饱满的水稻品种中花11种子，剥去种皮，灭菌洗涤浸泡后，均匀的点在灭菌r1固体培养基上，28℃黑暗处培养4～6周诱导愈伤形成；挑取状态较为紧实的愈伤颗粒至于新的r1继代培养基上，培养约1周后，用于农杆菌介导水稻愈伤遗传转化。
[0205]
已制备好的农杆菌(eha105感受态，购买于庄盟国际有限公司，货号为zc142)的单克隆接种于10ml含相应抗生素的lb液体培养基中，28℃振荡培养12小时；当培养至对数期时，收集农杆菌菌体，并用aam侵染液重悬农杆菌，将已培养好的水稻愈伤组织浸泡到此农杆菌重悬液中侵染30分钟；将侵染完的愈伤组织转入r2培养基上进行共培养3天；共培养结束后转至带有50mg/l潮霉素的r1筛选培养基上进行筛选培养，28℃培养4～6周；选取生长良好呈淡黄色的愈伤组织，用无菌镊子移至r4分化培养基，28℃持续光照培养；待分化出来的幼苗长至2-5厘米，转入r5生根培养基培养2～3周，获得抗性t0代植株(再生植株)，之后移入土中置于温室生长(温度28～30℃，16小时光照/8小时黑暗)。
[0206]
2、再生植株的基因型检测
[0207]
对步骤1中筛选获得的再生植株进行基因型检测及分析，使用ctab提取方法，提取t0代植株的基因组dna，以t0代植株基因组dna为模板，进行pcr扩增获得靶标基因片段，对筛选获得的再生植株进行基因型检测及分析。基因型检测所用引物如表8所示。
[0208]
表8、基因型检测所用引物列表
[0209]
[0210][0211]
pcr扩增体系：15μl 2
×
taq酶，1μl正向引物，1μl反向引物，1μl基因组dna模板(100ng/μl)，12μl ddh2o。
[0212]
pcr扩增程序：98℃，3min；98℃，15s；60℃，15s；72℃，30s；35个循环；72℃，5min；16℃保持。
[0213]
对pcr扩增产物进行sanger测序，与野生型基因组序列进行比对，鉴定编辑类型。
[0214]
3、paybes介导的水稻植株的基因型分析
[0215]
本发明在水稻植株中对不同paybes的编辑活性进行了测试。具体检测情况如图7、图8、图9、图10、表9、表10所示。
[0216]
paybev1-hmpg载体在osdep1靶点共得到86簇转基因t0代苗，其中有3簇幼苗(#12、#25、#55)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为3.49％，81簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为94.19％。在osepsps靶点共得到52簇转基因t0代苗，未检测到a-to-y编辑或indel事件，其中19簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为36.54％。在osnrt1.1b靶点共得到88簇转基因t0代苗，其中有3簇幼苗(#12、#15、#34)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为3.41％，剩下的79簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为89.77％。在oswaxy-t1靶点共得到53簇转基因t0代苗，未检测到a-to-y编辑或indel事件，其中34簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为64.15％。在oswaxy-t2靶点共得到99簇转基因t0代苗，其中有3簇幼苗(#31、#29、#48)在靶点处检测到了基因型为双等位状态的a-to-t的编辑事件，a-to-t编辑效率为2.02％，其中，#31有4株单株，#48和#29基因型相同，为同一愈伤来源，分别含有3株和1株单株。此外，有8簇幼苗(#23、#25、#27、#38、#48、#77、#78、#92)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为8.08％，其中9bp的插入小片段在靶点处是发生了复制，93簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为93.94％。在oswaxy-t3靶点共得到71簇转基因t0代苗，未检测到a-to-y编辑或indel事件，61簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为85.92％。同时由深度测序结果分析可知，精准删除与插入片段发生窗口集中在a4-a6至ncas9(d10a)切割位点，最多可以造成14bp的缺失与9bp的插入。这些结果表明，paybev1-hmpg载体可以实现a-to-y的碱基编辑，同时也能够造成小片段的精准插入和缺失。
[0217]
paybev3-hmpg-vp64载体在osdep1靶点共得到56簇转基因t0代苗，在靶点处未检测到a-to-y编辑事件的出现，其中有14簇幼苗(#49、#3、#12、#25、#32、#33、#34、#36、#38、#39、#52、#31、#37、#40)在靶点处检测到了indel事件的出现，indel事件发生效率为25.00％，有51簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为91.07％。paybev3-hmpg-vp64载体在osepsps靶点共得到76簇转基因t0代苗，在靶点处未检测到a-to-y编辑事件的出现，其中有1簇幼苗(#23)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为1.32％，有44簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为57.89％。paybev3-hmpg-vp64载体在osnrt1.1b靶点共得到76簇转基因t0代苗，其中有1簇幼苗(#5)基因型在靶点处检测到了a-to-t编辑事件，a-to-t编辑效率为1.32％，#5株系有3株单株，为双等位状态a6》t6编辑事件，有1簇幼苗(#3)在靶点处检测到了indel事件的出现，indel事件发生效率为1.32％，有64簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为84.21％。paybev3-hmpg-vp64载体在oswaxy-t1靶点共得到65簇转基因t0代苗，有1簇幼苗(#14)在靶点处检测到了基因型为a-to-t的编辑事件，a-to-t编辑效率为1.54％，#14株系含有单株2株，为a5》t5、a6》t6杂合状态编辑事件，此外，有13簇幼苗(#2、#25、#28、#30、#32、#55、#58、#49、#13、#20、#3、#7、#63)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为20.00％，有42簇幼苗均检测到了a-to-g编辑事件，a-to-g编
辑效率为64.62％。paybev3-hmpg-vp64载体在oswaxy-t2靶点共得到79簇转基因t0代苗，有3簇幼苗(#12、#3、#42)在靶点处检测到了基因型为a-to-t的编辑事件，a-to-t编辑效率为2.53％，相比paybev1-hmpg载体，在oswaxy-t2靶点a-to-t编辑效率提高1.25倍(2.53％/2.02％)。其中，#12株系为双等位状态a5》t5编辑事件，含有单株3株，#3与#42株系基因型相同视为同一个转化事件，为a5》t5的杂合编辑植株，分别含有单株1株、2株。此外，有15簇幼苗基因型在靶点处检测到了indel事件的出现，indel事件发生效率为18.99％，有66簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为83.54％。paybev3-hmpg-vp64载体在oswaxy-t3靶点共得到74簇转基因t0代苗，有7簇幼苗(#33、#40、#70、#43、#53、#54、#55)在靶点处检测到了基因型为a-to-t的编辑事件，a-to-t编辑效率为9.46％。其中，#33株系为含有t5》a5编辑事件的双等位编辑植株，含有单株2株，#40株系为含有t6》a6编辑事件的双等位编辑植株，含有单株2株，#70株系为含有t5》a5编辑事件的杂合编辑植株，含有单株2株，#43、#53、#54、#55等株系为嵌合状态的编辑植株，分别含有单株1株、2株、2株、1株单株；此外，有4簇幼苗基因型在靶点处检测到了indel事件的出现，indel事件发生效率为5.41％，有54簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为72.97％。
[0218]
paybev4-mhmpg-vp64载体在osdep1靶点共得到82簇转基因t0代苗，有9簇幼苗(#34-1、#35、#40-1、#42、#44-1、#44-2、#44-3、#65、#99)在靶点处检测到了a-to-t的编辑事件，a-to-t编辑效率为10.98％，相比paybev3-hmpg-vp64载体，在osdep1靶点首次实现稳定株系的编辑。其中，#42株系有1株单株，为含有a5》t5编辑事件的纯合编辑植株，#35株系有3株单株，与#44-1、#44-2、#44-3等株系为含有a5》t5编辑事件的杂合编辑植株，#65株系与#99株系分别含有3株与1株单株，#65与#34-1、#40-1等株系皆为含有a6》t6编辑事件的嵌合编辑植株。此外，有28簇幼苗(#1、#2、#3、#7、#8、#9、#11、#20、#22、#25、#31、#36、#38、#40、#41、#44、#45、#46、#47、#48、#55、#56、#57、#58、#63、#65、#71、#74)在靶点处检测到了indel事件的出现，indel事件发生效率为34.15％，其中多种类型的小片段插入是在靶点处发生了复制，删除片段多是集中在a5位点至ncas9切口处，有64簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为78.05％。paybev4-mhmpg-vp64载体在osepsps靶点共得到65簇转基因t0代苗，在靶点处未检测到a-to-y编辑事件的出现。其中有1簇幼苗(#23)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为1.54％，有33簇幼苗均检测到了a-to-g编辑事件的出现，a-to-g编辑效率为50.77％。paybev4-mhmpg-vp64载体在osnrt1.1b靶点共得到90簇转基因t0代苗，有2簇幼苗(#61、#89)在靶点处检测到了基因型为a-to-t的编辑事件，a-to-t编辑效率为2.22％，相比paybev3-hmpg-vp64载体，在osnrt1.1b靶点a-to-t编辑效率提高1.68倍(2.22％/1.32％)。其中，#61株系有4株单株，为含有a4》t4编辑事件的双等位编辑植株，#89株系有1株单株，为含有a4》t4编辑事件的嵌合编辑植株。此外，有5簇幼苗基因型在靶点处检测到了indel事件的出现，indel事件发生效率为5.56％，有74簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为82.22％。paybev4-mhmpg-vp64载体在oswaxy-t1靶点共得到71簇转基因t0代苗，有5簇幼苗(#64、#51、#2、#49、#50)在靶点处检测到了基因型为a-to-t的编辑事件，其中#2、#49、#50株系基因型一致，视为同一个转化事件，a-to-t编辑效率为4.23％，相比paybev3-hmpg-vp64载体，在oswaxy-t1靶点a-to-t编辑效率提高2.75倍(4.23％/1.54％)。#2、#49、#50株系分别含有单株3株、1株、2株，为含有t5》a5、t6》a6编辑事件的杂合编辑植株，#64株系有4株单株，为含有双t5》a5、t6》a6编辑事件的
双等位编辑植株，#51株系有2株单株，为含有t5》a5、t6》a6编辑事件的嵌合编辑植株。此外，有10簇幼苗(#2、#3、#7、#13、#20、#25、#28、#58、#59、#63)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为14.08％，有45簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为63.38％。paybev4-mhmpg-vp64载体在oswaxy-t2靶点共得到82簇转基因t0代苗，有7簇幼苗(#6、#20、#27、#32、#56、#39、#74)在靶点处检测到了基因型为a-to-t的编辑事件，a-to-t编辑效率为8.54％，相比paybev3-hmpg-vp64载体，在oswaxy-t2靶点a-to-t编辑效率提高3.38倍(8.54％/2.53％)。其中，#6、#20、#27、#32、#56为含有a5》t5编辑事件的双等位编辑植株，分别含有3株、1株、4株、3株、5株单株，#39、#74为含有a5》t5编辑事件的嵌合编辑植株，分别含有单株2株、1株。此外，有21簇幼苗(#2、#3、#4、#5、#7、#8、#10、#13、#15、#16、#20、#21、#25、#27、#28、#29、#81、#82)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为25.61％。有62簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为75.61％。paybev4-mhmpg-vp64载体在oswaxy-t3靶点共得到58簇转基因t0代苗，有10簇幼苗(#34、#6、#15、#37、#3、#4、#31、#36、#47、#51)在靶点处检测到了基因型为a-to-t的编辑事件，a-to-t编辑效率为17.24％，相比paybev3-hmpg-vp64载体，在oswaxy-t3靶点a-to-t编辑效率提高1.82倍(17.24％/9.46％)。其中，#34株系含有t6》a6编辑事件的双等位编辑植株，含有单株3株；#6、#15、#37等株系为含有t5》a5编辑事件的杂合编辑植株，分别含有2株、2株、3株单株；#3、#4、#31、#36、#47、#51株系含有t5》a5编辑事件的嵌合编辑植株，分别含有单株2株、1株、3株、2株、2株、1株。此外，有6簇幼苗(#15、#16、#25、#31、#47、#54)基因型在靶点处检测到了indel事件的出现，indel事件发生效率为10.34％，有37簇幼苗均检测到了a-to-g编辑事件，a-to-g编辑效率为63.79％。
[0219]
表9、水稻稳定转化株系paybes效率统计
[0220]
[0221][0222]
注：下划线处为pam序列。
[0223]
表10、a-to-y水稻稳定转化株系基因型列表
[0224]
[0225]
[0226]
[0227][0228]
注：下划线处为pam序列；括号内为每个株系含有单株数目；a-to-y碱基编辑标注为加粗；
[0229]
ho表示纯合；he表示杂合；bi表示双等位；chi表示嵌合。
[0230]
4、植物单碱基编辑a to y脱靶检测
[0231]
根据crispr-ge(http://skl.scau.edu.cn/)网站预测脱靶靶点，寻找序列并设计相应的引物对a-to-y编辑株系进行扩增鉴定，与osdep1、osepsps、osnrt1.1b、oswaxy-t1、oswaxy-t2、oswaxy-t3等靶点相似的脱靶位点，pcr扩增经sanger测序后，结果显示在预测的脱靶位点处并没有编辑事件的出现。详细脱靶位点序列信息见表11。
[0232]
表11、脱靶位点分析
[0233][0234]
注：下划线处为pam序列；加粗碱基为脱靶位点与靶点错配碱基。

技术特征：
1.成套系统，所述成套系统为如下m1)-m4)中任一种：m1)所述成套系统包括腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和grna；m2)所述成套系统包括腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和grna；m3)所述成套系统包括反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶hmpg和grna；m4)所述成套系统包括反式激活因子vp64、腺嘌呤脱氨酶tada8e、cas9(d10a)缺刻酶、人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg和grna；所述人源烷基腺嘌呤dna糖基化酶七突突变体mhmpg为将所述人源烷基腺嘌呤dna糖基化酶hmpg氨基酸序列的第163位由甘氨酸突变为精氨酸，且将第169位由天冬酰胺突变为丝氨酸，且将第198位由丝氨酸突变为丙氨酸，且将第202位由赖氨酸突变为丙氨酸，且将第203位由甘氨酸突变为丙氨酸，且将第206位由丝氨酸突变为丙氨酸，且将第210位由赖氨酸突变为丙氨酸后得到的蛋白质。2.根据权利要求1所述的成套系统，其特征在于：所述腺嘌呤脱氨酶tada8e为a1)或a2)：a1)氨基酸序列是序列2所示的蛋白质；a2)将序列2所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。3.根据权利要求1或2所述的成套系统，其特征在于：所述cas9(d10a)缺刻酶为b1)或b2)：b1)氨基酸序列是序列3所示的蛋白质；b2)将序列3所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。4.根据权利要求1-3任一所述的成套系统，其特征在于：所述人源烷基腺嘌呤dna糖基化酶hmpg为c1)或c2)：c1)氨基酸序列是序列6所示的蛋白质；c2)将序列6所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。5.根据权利要求1-4任一所述的成套系统，其特征在于：所述反式激活因子vp64为d1)或d2)：d1)氨基酸序列是序列9所示的蛋白质；d2)将序列9所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。6.根据权利要求1-5任一所述的成套系统，其特征在于：所述成套系统还包括筛选剂抗性蛋白；或，所述筛选剂抗性蛋白为潮霉素磷酸转移酶；或，所述潮霉素磷酸转移酶为e1)或e2)：e1)氨基酸序列是序列4所示的蛋白质；
e2)将序列4所示的氨基酸序列经过一个或几个氨基酸残基的取代和/或缺失和/或添加且具有相同功能的蛋白质。7.权利要求1-6任一所述的成套系统在如下s1)-s9)任一种中的应用：s1)植物基因组序列的编辑；s2)制备植物基因组序列的编辑的产品；s3)提高植物基因组序列的编辑效率；s4)制备提高植物基因组序列的编辑效率的产品；s5)拓展植物基因组序列的编辑窗口；s6)制备拓展植物基因组序列的编辑窗口的产品；s7)制备植物突变体；s8)植物育种；s9)制备植物育种的产品。8.如下t1)-t4)任一所述的方法：t1)植物基因组序列的编辑方法，包括如下步骤：使植物表达权利要求1所述的成套系统，以实现植物基因组序列的编辑；t2)提高植物基因组序列的编辑效率的方法，包括如下步骤：使植物表达权利要求1所述的成套系统，以实现提高植物基因组序列的编辑效率；t3)拓展植物基因组序列的编辑窗口的方法，包括如下步骤：使植物表达权利要求1所述的成套系统，以实现拓展植物基因组序列的编辑窗口；t4)植物突变体的制备方法，包括如下步骤：使植物表达权利要求1所述的成套系统，以获得植物突变体。9.根据权利要求1-6任一所述的成套系统或权利要求7所述的应用或权利要求8所述的方法，其特征在于：所述植物基因组序列的编辑包括将植物基因组序列中的碱基a颠换为碱基c或碱基t。10.根据权利要求1-6任一所述的成套系统或权利要求7所述的应用或权利要求8所述的方法，其特征在于：所述植物为x1)或x2)或x3)：x1)单子叶植物或双子叶植物；x2)禾本科植物；x3)水稻。

技术总结
本发明公开了新型植物碱基编辑器pAYBEs。本发明通过将腺嘌呤碱基编辑器与hMPG融合，在植物中建立了第一代A-to-Y单碱基编辑系统pAYBEv1，进一步通过将hMPG替换为mhMPG建立了第二代A-to-Y单碱基编辑系统pAYBEv2，在此基础上，又通过将反式激活因子Vp64分别融合在pAYBEv1和pAYBEv2中的TadA8e-nCas9(D10A)的N端，建立了第三代单碱基编辑系统pAYBEv3和第四代单碱基编辑系统pAYBEv4。本发明构建的植物AYBE单碱基编辑系统不仅丰富了植物单碱基编辑工具箱，也为水稻和其它农作物重要基因功能分析和遗传改良提供重要技术支撑。能分析和遗传改良提供重要技术支撑。


技术研发人员：夏兰琴 李玉才 李少雅 林勇 张阳军
受保护的技术使用者：中国农业科学院作物科学研究所
技术研发日：2023.05.22
技术公布日：2023/7/22