雄性不育基因ZmTGA10及其在创制玉米雄性不育系中的应用的制作方法

雄性不育基因zmtga10及其在创制玉米雄性不育系中的应用
技术领域
1.本发明属于植物生物技术育种领域，具体涉及雄性不育基因zmtga10及其在创制玉米雄性不育系中的应用。


背景技术：

2.玉米是我国第一大粮食作物，常年播种面积在5.5亿亩以上，玉米种业健康发展对保障国家粮食安全有着重大的战略意义。同时，玉米种业也是全球市值最大、商业化程度最高、科技含量最高的种业领域，是全球种业竞争的战略要地。我国玉米种业与国际领先水平相比，在科技创新和产业模式等方面，仍然存在巨大差距。一是受限于玉米雄性不育基础研究尚未取得根本性突破等因素，导致玉米自交系知识产权保护困难，造成近年来我国玉米种业长期存在跟随性、模仿性育种现象，重大新品种选育效率缓慢。二是玉米制种产业仍然处于依靠人工去雄为主的劳动密集型阶段，成本高，资源消耗巨大，种子质量难以保障。
3.玉米是杂种优势利用最成功的作物之一，雄性不育系是作物杂种优势利用和杂交制种的重要材料，根据雄蕊发育的特征，包括两种情况：一是不能产生正常的花粉，二是花粉育性正常但花药开裂异常不能及时散粉。根据控制的基因分类，雄性不育主要包括细胞质雄性不育（cms）和细胞核雄性不育（gms）。cms 是由线粒体基因和核基因共同控制的，虽然已应用于玉米育种和杂交种生产中，但存在资源利用率低、不育系胞质单一、易感病等问题。 gms 由核基因单独控制，可以克服 cms缺陷，但很难通过常规育种方法大量繁殖纯合不育系。近年来，随着生物技术的进步，通过基因工程和分子设计育种相结合创制的玉米多控不育技术以及植物通用型显性不育技术，可以有效地解决玉米隐性核雄性不育系的保持和繁殖问题。实现上述技术应用的重要前提是获得大量功能明确的控制玉米雄性发育的gms基因和对应的雄性不育材料。
4.与模式植物拟南芥和模式作物水稻相比，玉米中已克隆和鉴定的gms基因和创制的雄性不育材料均相对较少。crispr/cas9（clustered, regularly interspaced, short palindromicrepeats
‑
associated endonuclease 9）基因编辑技术由于具有成本低廉、操作简易和突变诱导率高等特点，被越来越广泛地应用于植物基因功能研究，作物遗传改良和育种等多个方面，应用前景十分广阔。利用crispr/cas9技术挖掘鉴定玉米雄性不育候选基因和创制雄性不育材料，可以快速丰富玉米gms基因和不育材料资源，从而促进玉米不育化育种和制种的推广和应用，最终可以有效解决我国玉米种业行业长期面临缺乏稳定的玉米不育系和突破性大品种的瓶颈问题。


技术实现要素：

5.针对现有技术的不足，本发明的目的是提供雄性不育基因zmtga10及其在创制玉米雄性不育系中的应用，可以用来创制玉米雄性不育系，从而应用于玉米杂交育种和制种。
6.为实现上述目的，本发明提供了zmtga10基因在控制玉米雄性生殖发育中的应用，其特征在于，所述基因的氨基酸序列如seq id no.2所示。一般可以预见这些来自不同植物
或不同玉米材料中的同源基因具有相同或者相似的功能，因此同样可以利用这些基因改良植物的农艺性状。进一步，即使不能预见这些基因的功能，本领域的一般技术人员可以根据本发明提供的方法和现有技术测定它们是否具有控制植物雄性育性的功能。
7.在另一方面，本发明还提供了zmtga10基因在控制玉米雄性生殖发育中的应用，其特征在于，所述基因的核苷酸序列为seq id no.1所示。
8.在另一方面，本发明还提供了一种创制玉米雄性不育系的方法，其特征在于，抑制玉米中zmtga10基因的表达和/或活性，选择玉米雄性不育的植株。
9.在一些实施方案中，上述抑制基因表达和/或活性的方法包括基因编辑、rna干扰、t
‑
dna插入中任一种。
10.在一些实施方案中，上述基因编辑采用crispr/cas9方法。
11.在一些实施方案中，所述crispr/cas9方法包括：在所述基因第四外显子、第8内含子和第九外显子处设计crispr/cas9载体靶点，所述靶点的dna序列如seq id no.3，seq id no.4，seq id no.5或seq id no.6所示。
12.在另一方面，本发明还提供一种获得tga10雄性不育系的方法，将通过上述方法获得的tga10雄性不育系与目标材料进行杂交和回交，从而使目标材料获得tga10雄性不育的性状和基因突变。
13.本发明还包括通过上述任一方法获得的tga10雄性不育系在杂交育种和制种中的应用。所述在杂交育种和制种中的应用是指将tga10雄性不育系作为母本与其他父本进行杂交，或是将获得的tga10雄性不育系与其他目标材料进行杂交和回交，从而使目标材料获得tga10雄性不育的性状和基因突变。
14.更进一步地，本发明还提供了两种玉米雄性不育系tga10的分子标记引物，引物zmtga10
‑
f1和zmtga10
‑
r1序列分别如seq id no.7和seq id no.8所示；引物zmtga10
‑
f2和zmtga10
‑
r2序列分别如seq id no.9和seq id no.10所示。
15.本发明的优点及有益效果如下：zmtga10（zm00001d020938）基因及编码的蛋白调控玉米雄性生殖发育是之前没有报道过的。本发明利用crispr/cas9方法突变玉米基因tga10（zm00001d020938），发现了zmtga10（zm00001d020938）基因对玉米花药发育的调控功能。利用crispr/cas9基因编辑的方法和编辑后获得的雄性不育突变体，可以创制玉米雄性不育系，从而可以应用于玉米杂交育种和制种。针对三种tga10的雄性不育系开发共分离分子标记，可用于植株的育性等位基因鉴定、分子标记辅助育种中目标单株的筛选和种子纯度鉴定等。
附图说明
16.图1为zmtga10基因在玉米不同发育时期花药中的表达模式分析s5，造孢细胞时期；s6, 小孢子母细胞时期；s7，减数分裂开始时期；s8a，减数分裂ⅰ，二分体时期；s8b，减数分裂ⅱ，四分体时期；s8b
‑
9，四分体
‑
单核小孢子时期；s9，单核小孢子时期；s9
‑
10，单核小孢子
‑
小孢子空泡化时期；s10，小孢子空泡化时期；s11，小孢子第一次不均等有丝分裂，二核小孢子时期；s12，小孢子第二次有丝分裂，三核小孢子时期。
17.图2为pcas9
‑
zmtga10定点突变表达载体的物理图谱pcas9
‑
zmtga10
‑
a：从t
‑
dna的左边界到右边界分别是除草剂抗性基因bar的表达
盒；核酸酶编码基因cas9的表达盒；zmtga10基因靶标2（mt2）的表达盒；靶标1（mt1）的表达盒。pcas9
‑
zmtga10
‑
b：从t
‑
dna的左边界到右边界分别是除草剂抗性基因bar的表达盒；核酸酶编码基因cas9的表达盒；zmtga10基因靶标4（mt4）的表达盒；靶标3（mt3）的表达盒。
18.图3为野生型zmtga10与其不育突变体的基因结构和dna序列分析野生型zmtga10（wt
‑ꢀ
zmtga10）：基因全长9447 bp，包括10个外显子和9个内含子；tga10突变体zmtga10
‑
cas9
‑
1：在第8个内含子和第9个外显子4368 bp
‑
4534 bp处缺失167个碱基；突变体zmtga10
‑
cas9
‑
2：在第8个内含子和第9个外显子4373 bp
‑
4534 bp之间缺失162个碱基；突变体zmtga10
‑
cas9
‑
3：在第4个外显子3114 bp处插入1个碱基（c）和在3246 bp处插入1个碱基（c）。
19.图4为野生型与tga10纯合突变体的雄穗、花药和花粉粒表型分析上排为玉米野生型（wt）与zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2、zmtga10
‑
cas9
‑
3突变体雄穗的表型比较；第二排为wt与zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2、zmtga10
‑
cas9
‑
3突变体花药的表型比较；下排为wt与zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2、zmtga10
‑
cas9
‑
3突变体花粉粒的i2‑
ki染色比较。
20.图5为野生型与tga10纯合突变体的花药扫描电镜（sem）分析从左至右依次为：野生型（wt）花药整体；tga10花药整体；剥开后的wt（上）和tga10（下）花药；wt的成熟花粉粒（上）和tga10无法扫描到花粉粒（下）；wt（上）和tga10（下）花药的外表皮角质层；wt（上）和tga10（下）花药的内表皮乌氏体。
21.图6为利用共分离标记对zmtga10
‑
cas9
‑
1不育系的f2代植株进行基因分型共分离标记zmtga10
‑
f1/r1对8株zmtga10
‑
cas9
‑
1不育系f2代植株的pcr和琼脂糖凝胶电泳鉴定结果：在纯合野生型（aa）植株中扩增出445 bp条带；在tga10/ tga10杂合型（aa）植株中扩增出445 bp和278 bp两条带；在tga10/ tga10纯合突变型（aa）植株中扩增出278 bp条带。
22.图7为利用共分离标记对zmtga10
‑
cas9
‑
2不育系的f2代植株进行基因分型共分离标记zmtga10
‑
f1/r1对6株zmtga10
‑
cas9
‑
2不育系f2代植株的pcr和琼脂糖凝胶电泳鉴定结果：在纯合野生型（aa）植株中扩增出445 bp条带；在tga10/ tga10杂合型（aa）植株中扩增出445 bp和283 bp两条带；在tga10/ tga10纯合突变型（aa）植株中扩增出283 bp条带。
23.图8为利用共分离标记对zmtga10
‑
cas9
‑
3不育系的f2代植株进行基因分型共分离标记zmtga10
‑
f2/r2对5株zmtga10
‑
cas9
‑
3不育系f2代植株的pcr和聚丙烯酰胺凝胶电泳（page）电泳鉴定结果：在纯合野生型（aa）植株中扩增出74 bp条带；在tga10/ tga10杂合型（aa）植株中扩增出74 bp和75 bp两条带；在tga10/ tga10纯合突变型（aa）植株中扩增出75 bp条带。
具体实施方式
24.下述实施例用来说明本发明，但不限制本发明的范围。在不背离本发明精神和实质的情况下，对本发明方法、步骤或条件所作的修改或替换，均属于本发明的范围。如无特殊说明，实施例中所用引物及基因的合成和测序均由生工生物工程（上海）股份有限公司完成。其他生化试剂非特别注明外均为常规市售试剂，实施例中所用的技术手段为本领域技
术人员所熟知的常规手段。
25.实施例一玉米tga10（zm00001d020938）基因序列和表达模式分析在maizegdb 库（https://www.maizegdb.org/）中，查询到玉米tga10（zm00001d020938，grmzm2g006578）基因，在b73中的核酸序列见seq id no.1所示，基因功能注释为bzip7转录因子（bzip
‑
transcription factor 7，bzip7)，其编码蛋白包含486个氨基酸，序列如seq id no.2所示。
26.bzip转录因子参与植物中众多生理过程的调节，由于zmtga10基因在玉米中的实际功能还没有公开的研究资料，为了研究该基因与玉米雄性生殖发育的关系，本发明首先利用qrt
‑
pcr分析了该基因在玉米花药发育不同阶段的表达模式。具体步骤如下：1、玉米花药的取样与发育时期鉴定从玉米自交系b73处于不同发育阶段的雄穗中，按照花药的长度，收集不同长度的花药样品；每份样品采集20个长度相近的新鲜花药，其中3个固定在faa溶液中（coolaber，中国）通过树脂半薄切片实验确定具体的发育阶段，其余17个花药立即冷冻在液氮中，用于rna的提取。
27.利用梯度乙醇（50%、70%、90%、100%）对用于树脂切片的固定花药进行脱水，每步15
‑
30分钟。脱水过程中花药可置于70%乙醇中长期保存；为便于后期包埋，可在90%乙醇中加入0.1%的伊红对材料进行染色；为保证脱水彻底，材料须在无水乙醇中脱水2
‑
3次。随后进行树脂替换，将花药依次置于乙醇与spurr树脂体积比为3:1、1:1、1:3液体中2
‑
4小时，最后置于纯树脂中过夜。树脂置换完成后，将花药置于模具中，加入200
ꢀµ
l spurr树脂，置于烘箱中，70℃聚合过夜。随后进行修块，然后可利用德国莱卡切片机进行切片，切片厚度为2
ꢀµ
m；用镊子夹取切好的片子，置于载玻片中央的无菌水中，42℃条件下展片过夜。将固定有样品的载玻片浸入0.1%的甲苯胺蓝染液中，染色1分钟，然后用去离子水冲洗，再置于展片台上，烘干后即可用于显微观察；也可以封片后长期保存。分析树脂切片的结果，根据玉米14个不同发育时期（stage1
‑
stage14：s1
‑
s14）的细胞学特点，确定每份样品的具体发育时期。
28.2、qrt
‑
pcr分析用trizol试剂（invitrogen，美国）提取上述鉴定的处于不同发育时期（s5
‑
s12）的玉米花药总rna；然后使用5x all
‑
in
‑
one rt master mix （abm，加拿大）合成cdna；采用tb green
™
premix ex taq
™
（takara，日本）在quantstudio5quantstudio 5 real
‑
time pcr system （abi，美国）上进行定量逆转录聚合酶链反应检测，扩增引物为：qtga10
‑
f（5
’‑ꢀ
tccgcctcttgttccttgtc
‑3’
）和qtga10
‑
r（5
’‑ꢀ
gacctaggctgcttggatcc
‑3’
）；zmactin1为参照基因，其扩增引物为：actin1
‑
f（5
’‑ꢀ
aaatgacgcagattatgtttga
‑3’
）和actin1
‑
r（5
’‑ꢀ
gctcgtagtgagggagtacc
‑3’
）；每个发育时期包括三个生物学重复，每个样品有三个技术重复；数据用2
‑
δδct
方法进行分析，并以均值
±
标准差（means
±
sd）的形式给出定量结果。
29.zmtga10基因呈现花药发育时期特异表达的模式：在玉米花药发育的早期比如s5时期有较高表达，随后开始减弱，然后到花药发育的中后期（s9）几乎没有表达（图1）。
30.实施例二玉米tga10（zm00001d020938）基因的功能以及利用crispr/cas9方法创制玉米雄性不育系为了明确玉米tga10（zm00001d020938）在玉米中的功能，本发明采用crispr/cas9
基因编辑方法突变zm00001d020938基因序列，敲除该基因在玉米中的功能。本发明选取玉米杂交种hi ii 作为基因编辑的受体材料。本发明分别选取基因保守区的seq id no.3、seq id no.4、seq id no.5和seq id no.6所示序列作为crispr/cas9基因编辑的靶标区域。
31.1、zmtga10的crispr/cas9基因编辑载体的构建本发明的基因编辑载体为pbue411
‑
mt1t2
‑
cas9，该载体的基础载体为pbue411
‑
cas9，中间载体为pcbcmt1t2，提供grna。本发明通过在引物上设计靶点然后通过pcr获得mt
‑
sgrna进而通过酶切连接到基础载体中，具体的构建流程如下：（1）靶标grna的设计。将zmtga10（zm00001d020938）的基因序列输入http://crispr.hzau.edu.cn/cgi
‑
bin/crispr2/crispr进行靶标设计。本发明选择的四个靶标区域的dna序列如seq id no.3、seq id no.4、seq id no.5和seq id no.6所示。本发明的sgrna骨架序列从中间载体pcbcmt1t2直接扩增获得。
32.（2）通过在引物上设计靶点然后pcr扩增获得mt
‑
sgrna。引物zmtga10
‑
mt1
‑
f和引物zmtga10
‑
mt2
‑
r扩增中间载体pcbcmt1t2，用于获得包含第一个和第二个靶标的sgrna的片段，同样的步骤使用引物zmtga10
‑
mt3
‑
f和引物zmtga10
‑
mt4
‑
r用于扩增获得包含第三个和第四个靶标的sgrna的片段，产物长度均为891 bp。pcr 体系和条件如下：模板dna （中间载体pcbcmt1t2 ≥30 ng/μl）1.2 μl；primer f/r：各1.2 μl；灭菌ddh2o：11.4 μl；2x mclab 酶（产品编号：i5hm
‑
200）：15 μl。pcr的温度程序如下：
①
98 ℃ 2分钟；
②
98 ℃ 10秒；
③
58 ℃ 30秒；
④
72 ℃ 30秒；
⑤ꢀ
从
②
－
④
循环34次；
⑥
72 ℃ 5分钟；
⑦
25 ℃ 10分钟。最后回收pcr产物。载体构建所需的引物序列如下：zmtga10
‑
mt1
‑
f:5
’‑
atatatggtctctggcgagtcctgtccccggcgccgggttttagagctagaaatagcaa
‑3’
zmtga10
‑
mt2
‑
r:5
’‑
attattggtctctaaacctgctgagccagcagcacctgcttcttggtgccgc
‑3’
zmtga10
‑
mt3
‑
f:5
’‑
atatatggtctctggcgatacgtcggccatatatacggttttagagctagaaatagcaa
‑3’
zmtga10
‑
mt4
‑
r:5
’‑
attattggtctctaaacccatctcgtcgtggtgagctgcttcttggtgccgc
‑3’
（3）通过酶切连接构建到骨架载体。将pbue411
‑
cas9载体和回收的带有靶标的sgrna片段用bsai消化，同时加入t4连接酶将载体和sgrna片段连接。15 μl的酶切连接体系如下，sgrna片段：2 μl，pbue411
‑
cas9载体（≥60 ng/μl) ：2 μl，10 x neb buffer：1.5 μl，bsai内切酶（产品编号：#r3733s ）：1 μl，t4 连接酶（产品编号：#m0202m ）：1 μl，灭菌ddh2o：6 μl。
33.图2所示为目的基因zmtga10（zm00001d020938）的双靶标（对应第一个靶标和第二个靶标），标记基因cas9和bar与骨架载体pbue411
‑
cas9构建的表达载体pcas9
‑
zmtga10
‑
a；第三个和第四个靶标对应的最终表达载体pcas9
‑
zmtga10
‑
b除靶标序列有差异外其余与pcas9
‑
zmtga10
‑
a相同。
34.2、农杆菌介导的玉米遗传转化将上述构建的pcas9
‑
zmtga10
‑
a和pcas9
‑
zmtga10
‑
b载体分别通过热激法转入农
杆菌eha105中，pcr进行鉴定；然后以1:1的浓度混合分别含有两种敲除载体的农杆菌，混合后加甘油于
‑
80 ℃保存菌液。以新鲜剥离的1.5 mm左右的玉米杂交种hi ii的幼胚作为受体材料，将剥离的玉米胚放入含有1.8 ml悬浮液的2 ml塑料离心管中，放置时间不超过1小时，每个离心管中大约放入100个幼胚；吸去悬浮液，并用新的悬浮液将幼胚清洗2遍，管底保留少量可没过幼胚的悬浮液，然后43 ℃热激2分钟，之后再冰浴1分钟，用移液枪吸尽管底残留的洗液，并加入1.0 ml的农杆菌侵染液，轻摇30秒，然后黑暗静置8分钟。接下来将离心管中的幼胚和侵染液倒入共培养基上，晃匀后用移液枪吸出多余的侵染液，所有幼胚的盾片朝上，于23 ℃黑暗共培养3天。共培养结束后，用无菌的镊子将幼胚转移至恢复培养基，于28 ℃培养7
‑
14天，中间过程需留意及时去掉幼胚上长出的幼芽。恢复培养结束后将幼胚放至含1.5 mg/l bialaphos筛选培养基上筛选培养3轮，每轮筛选2周，然后转到2 mg/l bialaphos筛选培养基上筛选培养2轮，每轮筛选2周。将抗性愈伤组织转至扩繁培养基，28 ℃，暗培养2周。随后将扩繁好的抗性愈伤组织转移至诱导培养基，于28 ℃黑暗环境下培养2周。然后转到分化培养基中，25 ℃，5000 lx，光照培养2周。培养结束后将分化出的苗簇分离出单个幼苗放置于生根培养基中，25 ℃，5000 lx，光照培养直到生根；将小苗转入小营养钵中生长，生长成活后移栽于温室中，3
‑
4个月后收获后代种子。
35.3、t0代植株crispr/cas9突变结果检测为确定t0代植株crispr/cas9突变结果，采取如下步骤进行：本发明首先采用ctab法提取玉米叶片dna，具体方法如下：剪取2厘米左右长度的幼苗叶片，放入装有钢珠的2 ml离心管；将装有叶片的离心管放入液氮浸没5分钟，然后利用研磨仪打碎叶片样品；向离心管中加入700 μl ctab提取缓冲液（其中含有1%的β
‑
巯基乙醇），并用力震荡混匀， 65 ℃恒温水浴中预热20
‑
30 min，（期间取出颠倒1
‑
2次，并注意实验样品编号的对应）；离心管冷却至室温后加入700 μl氯仿：异戊醇（24:1）萃取液，用力摇晃30s之后在室温条件下静置片刻；在4 ℃条件下，以12000 rpm速度离心5 min，取离心完成后的500 μl上清液于新的1.5 ml离心管中；将等体积的异丙醇加入到盛有上清液的离心管中轻轻震荡混匀，室温条件下静置10 min左右；随后将盛有样品的离心管放入4 ℃离心机中，以12000 rpm速度离心10 min之后，轻轻吸取上清液，弃上清，保留沉淀；加入800 μl 75 %乙醇，洗涤沉淀两次，以10000 rpm速度离心5 min，弃掉上清；将样品放置在室温下自然干燥2
‑
4小时，得到dna沉淀，加入适量无菌水溶解，轻微震荡，充分溶解dna。dna样品存于
‑
20 ℃保存。用nanodrop检测dna浓度，稀释至10 ng/l用作pcr模板。
36.然后根据zmtga10（zm00001d020938）基因序列设计pcr引物。
37.（1）检测靶标：mt1和mt2；产物大小：435 bp；引物序列如下：zmtga10
‑
t
‑
f1:5
’‑
cagatttgttcaagcactcgatgtc
‑3’
；zmtga10
‑
t
‑
r1:5
’‑
cagcaatggtcaaatcgagagag
‑3’
。
38.（2）检测靶标：mt3和mt4；产物大小：445bp；引物序列如下：zmtga10
‑
t
‑
f2:5
’‑
cggcctcagctcaggtacaa
‑3’
；zmtga10
‑
t
‑
r2:5
’‑
gtaagggagggagggacacgt
‑3’
。
39.提取基因组dna，按照以下pcr参数扩增：反应体系：15 μl mix常规pcr体系，0.5 μl forward primer，0.5 μl reverse primer，1 μl dna，5.5 μl 灭菌ddh2o，7.5 μl 2x taq mix（产品编号：10103es ）。
40.反应程序：常规pcr：58 ℃退火、延伸1分钟、32轮循环。
41.接着将pcr产物回收并连接t载体测序，通过测序多个t0代独立阳性转化事件靶标区域的dna序列，明确靶标区域是否发生基因编辑，最终发现3个t0转化事件靶标区域序列发生了变化，且均为同源突变，编辑前后的序列如图3所示，对应3个tga10突变体：zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2、zmtga10
‑
cas9
‑
3。与野生型序列比对显示zmtga10
‑
cas9
‑
1在靶标3和4处发生了缺失突变，zmtga10
‑
cas9
‑
2在靶标3和4处同样发生了缺失突变，而zmtga10
‑
cas9
‑
3在靶标1和2处发生了插入突变。
42.对3个tga10突变体中氨基酸序列进行比对分析发现，与未编辑的wt相比，突变后的品系zmtga10
‑
cas9
‑
1和zmtga10
‑
cas9
‑
2，其编码的的核苷酸在靶标3或4处的缺失一部分发生在第8个内含子上，在第9个外显子上共同有147 bp 的碱基缺失，使其丢失了47个氨基酸。另外zmtga10
‑
cas9
‑
3品系中编码的核苷酸在靶标1和2处的插入和缺失也使得氨基酸发生移码并导致氨基酸翻译提前终止。因此这些转化体中zm00001d020938蛋白的功能均表现缺失。
43.4、f1代植株的基因分型由于在温室生长的玉米t0代植株，经常雌穗和雄穗发育不协调，同时当编辑的基因与雄性发育相关时也会影响育性，因此为了繁殖t0代植株，并使获得的基因编辑类型遗传下去，本发明使用玉米自交系郑58的野生型花粉为上述获得的zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2、zmtga10
‑
cas9
‑
3的t0代植株授粉，进而获得f1代种子，生长的植株为f1代植株。
44.f1代植株包括2种分离类型，一种为cas9
‑
阳性植株（转基因植株），另一种为cas9
‑
阴性植株（非转基因植株），为了避免sgrna和cas9对杂交授粉引入的郑58野生型等位基因进行持续编辑，从而造成突变类型的复杂性，我们需要通过基因分型从f1代植株中挑选出不含有cas9基因，但含有t0代突变类型的植株，这类植株自交后可以得到非转基因的f2代。f1代植株的基因分型步骤如下：按照上述的ctab法提取叶片dna后，首先利用cas9基因的特异引物cas9
‑
f（5
’‑
cccggacaatagcgatgt
‑3’
）和cas9
‑
r（5
’‑ꢀ
gagtgggccgacgtagta
‑3’
）进行pcr扩增。pcr反应体系同上；反应程序：常规pcr：58℃退火、延伸1分钟、32轮循环。pcr产物进行琼脂糖凝胶电泳后，根据结果区分出cas9
‑
阳性植株和cas9
‑
阴性植株。
45.进一步针对cas9
‑
阴性植株，采用上述的检测mt1和mt2靶标的引物zmtga10
‑
t
‑
f1 和zmtga10
‑
t
‑
r1，以及检测mt3和mt4靶标的引物zmtga10
‑
t
‑
f2 和zmtga10
‑
t
‑
r2进行pcr扩增；pcr产物纯化后，连接t载体，进行测序；根据测序结果分析确定t0代突变类型的遗传情况。
46.实施例三tga10不育系的表型分析上述实施例二鉴定的不含有cas9基因的f1代植株自交后获得f2代种子，三种tga10突变类型（zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2和zmtga10
‑
cas9
‑
3）各取1个自交单穗进行穗行播种，在成熟期进行表型的调查。三种f2株系中，可育株与不育株的比例，均符合3:1分离，进一步表明tga10不育系的不育性状由单个隐性基因控制，然后针对f2代获得的稳定非转基因tga10不育系与野生型进行详细的表型比较。
47.1、雄穗、花药和花粉活力的观察
在营养生长和雌穗的发育方面，tga10不育系（zmtga10
‑
cas9
‑
1、zmtga10
‑
cas9
‑
2和zmtga10
‑
cas9
‑
3）的植株与野生型相比基本无差异；在雄穗发育方面，野生型能够正常抽雄，花药能够正常开裂、散粉，自交后可正常结实，而tga10不育系虽能正常抽雄，也能正常开花，但是花药完全不开裂，不能正常散粉和自交（图4）；进一步对野生型和突变体的花粉进行i2‑
ki染色，发现突变体和野生型花粉发育正常，花粉粒染色后呈黑色（图4）。这表明zmtga10（zm00001d020938）基因通过影响花药的开裂来控制玉米的雄性发育，通过基因编辑方法创制的tga10不育系为花药不开裂型的不育系。
48.2、花药的扫描电镜（sem）观察为了深入解析tga10的细胞学特征，对野生型和突变体花药内外壁进行扫描电镜（sem）分析。剥取处于成熟期（s13）的野生型和突变体花药，然后立即固定在faa（coolaber，中国）溶液中，固定液的体积不少于所取研究材料体积的20倍；对于突变体花药，可用解剖针在花药壁穿孔以提高固定液的渗透效果，或者反复抽真空至花药沉入固定液底部；室温固定2小时后，将材料置于4℃保存，或依次置于50%、60%、70%、80%、90%、100%的乙醇中进行脱水，每个梯度保持15分钟；材料可置于70%的乙醇中过夜或保存。脱水后的样品进行二氧化碳临界点干燥，然后镀金即可进行观察。发现tga10突变体花药内的花粉形态完整和饱满，但表面附着有部分残渣。突变体的花药外表皮和野生型一样能形成致密的网状角质层结构； tga10突变体花药的内表皮同样附着有少量残渣，但颗粒状的乌氏小体能正常形成。花药角质层是覆盖在花药表面的细胞外脂质层，保护花药免受外部非生物胁迫、内部组织失水和病原体的侵袭，位于花药内壁的乌氏体，被认为是孢粉素前体从绒毡层细胞到小孢子的转运载体。上述结果表明zmtga10（zm00001d020938）基因突变后，不影响花药角质层的形成和以及绒毡层内孢子花粉素前体物质的合成，只是影响花药的开裂。
49.实施例四tga10不育系鉴定的共分离分子标记开发和应用1、共分离分子标记的开发在本发明中，针对获得的三种tga10不育系的突变位点，利用primer5.0软件进行引物设计，开发出两对共分离分子标记：zmtga10
‑
f1/r1和zmtga10
‑
f2/r2，结合pcr与聚丙烯酰胺凝胶电泳（page）或琼脂糖凝胶电泳检测的方法，根据获得的条带及大小即可分离出突变体的基因型。
50.共分离分子标记zmtga10
‑
f1/r1包括第一引物zmtga10
‑
f1和第二引物zmtga10
‑
r1；该标记能够特异性检测玉米zmtga10
‑
cas9
‑
1突变体及由其转育的玉米不育材料中的突变基因tga10，并能同时区分野生型tga10基因和突变型tga10基因；针对突变基因tga10中扩增出278 bp的条带，而对野生型tga10基因扩增出445 bp的条带；该标记也能够特异性检测玉米zmtga10
‑
cas9
‑
2突变体及由其转育的玉米不育材料中的突变基因，并能同时区分野生型tga10基因和突变型tga10基因；针对突变基因tga10中扩增出283 bp的条带，而对野生型tga10基因扩增出445 bp的条带引物序列如下：zmtga10
‑
f1：5
’‑
cggcctcagctcaggtacaa
‑3’
zmtga10
‑
r1：5
’‑
gtaagggagggagggacacgt
‑3’
共分离分子标记zmtga10
‑
f2/r2包括第一引物zmtga10
‑
f2和第二引物zmtga10
‑
r2，该标记能够特异性检测玉米zmtga10
‑
cas9
‑
3突变体及由其转育的玉米不育材料中的突变基因tga10，并能同时区分野生型tga10基因和突变型tga10基因；针对突变基因tga10中
扩增出75 bp的条带，而对野生型tga10基因扩增出74bp的条带。引物序列如下：zmtga10
‑
f2：5
’‑
tggtagggggagattagtttagt
‑3’
zmtga10
‑
r2：5
’‑
tgatggtccggtgacctag
‑3’
2、共分离分子标记的应用为了验证上述标记的有效性，以实施例三中获得的f2株系为材料，进行tga10等位基因的检测。dna提取方法、pcr扩增体系和条件同实施例二，pcr产物进行page或琼脂糖凝胶电泳分离。
51.理论上，zmtga10
‑
f1/r1和zmtga10
‑
f2/r2在tga10/ tga10纯合野生型（aa）dna中能分别扩增出445 bp 、445 bp和74 bp的条带，在tga10/ tga10纯合突变型材料（aa）dna中分别扩增出278 bp、283 bp和75 bp的条带，而在tga10/ tga10杂合型（aa）材料中，则能分别同时扩增出对应的两条条带。zmtga10
‑
f1/r1和zmtga10
‑
f2/r2分子标记的验证结果如图6、图7和图8所示，结果显示设计的2对功能性分子标记对f2植株的检测结果完全符合预期，在tga10/ tga10纯合野生型（aa）、tga10/ tga10杂合型（aa）和tga10/ tga10纯合突变型材料（aa）中分别扩增出对应大小的条带，可以作为tga10等位基因检测的理想标记。
52.这些分子标记有助于在开花和授粉前确定突变基因型，以便在不同遗传背景下进行杂交和回交选育雄性不育系，具有重要的应用价值。