CN113005128A

CN113005128A - 雄性不育基因ZmMYB84及其在创制玉米雄性不育系中的应用

Info

Publication number: CN113005128A
Application number: CN202110272146.1A
Authority: CN
Inventors: 万向元; 颜廷玮; 安学丽; 江易林; 刘欣泽; 朱涛涛; 李紫文; 董振营; 吴锁伟; 李金萍
Original assignee: Beijing Shou Jia Li Hua Sci Tech Co ltd; Beijing Zhongzhi International Research Institute Of Biological Agriculture; University of Science and Technology Beijing USTB
Current assignee: Beijing Shou Jia Li Hua Sci Tech Co ltd; Beijing Zhongzhi International Research Institute Of Biological Agriculture; University of Science and Technology Beijing USTB
Priority date: 2021-03-12
Filing date: 2021-03-12
Publication date: 2021-06-22
Anticipated expiration: 2041-03-12
Also published as: CN113005128B

Abstract

本发明公开了雄性不育基因ZmMYB84及其在创制玉米雄性不育系中的应用，具有SEQ ID NO.1所示的核苷酸序列，该基因编码的蛋白质具有SEQ ID NO.2所示的氨基酸序列。本发明通过CRISPR/Cas9基因编辑技术在野生型玉米中定点突变该基因，可以导致完全的雄性不育，发现了ZmMYB84基因对玉米雄性生殖发育的调控功能。通过后代筛选，可以得到不含转基因成分的不育系，创制出稳定的玉米雄性不育系，对玉米雄性育性控制和杂交种子生产具有重要意义。本发明还针对获得的myb84雄性不育突变体设计了功能分子标记，在玉米雄性不育系培育、不育化杂交制种和分子标记辅助选择中具有重要的应用价值。

Description

雄性不育基因ZmMYB84及其在创制玉米雄性不育系中的应用

技术领域

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

背景技术

玉米是我国第一大粮食作物，常年播种面积在5.5亿亩以上，玉米种业健康发展对保障国家粮食安全有着重大的战略意义。同时，玉米种业也是全球市值最大、商业化程度最高、科技含量最高的种业领域，是全球种业竞争的战略要地。我国玉米种业与国际领先水平相比，在科技创新和产业模式等方面，仍然存在巨大差距。一是受限于玉米雄性不育基础研究尚未取得根本性突破等因素，导致玉米自交系知识产权保护困难，造成近年来我国玉米种业长期存在跟随性、模仿性育种现象，重大新品种选育效率缓慢。二是玉米制种产业仍然处于依靠人工去雄为主的劳动密集型阶段，成本高，资源消耗巨大，种子质量难以保障。

玉米是杂种优势利用最成功的作物之一，雄性不育系是作物杂种优势利用和杂交制种的重要材料，主要包括细胞质雄性不育（CMS）和细胞核雄性不育（GMS）。CMS 是由线粒体基因和核基因共同控制的，虽然已应用于玉米育种和杂交种生产中，但存在资源利用率低、不育系胞质单一、易感病等问题。 GMS 由核基因单独控制，可以克服 CMS缺陷，但很难通过常规育种方法大量繁殖纯合不育系。近年来，随着生物技术的进步，通过基因工程和分子设计育种相结合创制的玉米多控不育技术以及植物通用型显性不育技术，可以有效地解决玉米隐性核雄性不育系的保持和繁殖问题。实现上述技术应用的重要前提是获得大量功能明确的控制玉米雄性发育的GMS基因和对应的雄性不育材料。

与模式植物拟南芥和模式作物水稻相比，玉米中已克隆和鉴定的GMS基因和创制的雄性不育材料均相对较少。CRISPR/Cas9（Clustered, Regularly Interspaced, ShortPalindromicRepeats-associated Endonuclease 9）基因编辑技术由于具有成本低廉、操作简易和突变诱导率高等特点，被越来越广泛地应用于植物基因功能研究，作物遗传改良和育种等多个方面，应用前景十分广阔。利用CRISPR/Cas9技术挖掘鉴定玉米雄性不育候选基因和创制雄性不育材料，可以快速丰富玉米GMS基因和不育材料资源，从而促进玉米不育化育种和制种的推广和应用，最终可以有效解决我国玉米种业行业长期面临缺乏稳定的玉米不育系和突破性大品种的瓶颈问题。

发明内容

针对现有技术的不足，本发明的目的是提供雄性不育基因ZmMYB84及其在创制玉米雄性不育系中的应用，可以用来创制玉米雄性不育系，从而应用于玉米杂交育种和制种。

为实现上述目的，本发明提供了ZmMYB84基因在控制玉米雄性生殖发育中的应用，其特征在于，所述基因的氨基酸序列如SEQ ID NO.2所示。一般可以预见这些来自不同植物或不同玉米材料中的同源基因具有相同或者相似的功能，因此同样可以利用这些基因改良植物的农艺性状。进一步，即使不能预见这些基因的功能，本领域的一般技术人员可以根据本发明提供的方法和现有技术测定它们是否具有控制植物雄性育性的功能。

在另一方面，本发明还提供了根据权利要求1所述的应用，其特征在于，所述基因的核苷酸序列为SEQ ID NO.1所示。

在另一方面，本发明还提供了一种创制玉米雄性不育系的方法，其特征在于，抑制玉米中权利要求1或2所述基因的表达和/或活性，选择玉米雄性不育的植株。

在一些实施方案中，上述抑制蛋白表达和/或活性的方法包括基因编辑、RNA干扰、T-DNA插入中任一种。

在一些实施方案中，上述基因编辑采用CRISPR/Cas9方法。

在一些实施方案中，所述CRISPR/Cas9方法包括：在所述基因第一外显子处设计CRISPR/Cas9载体靶点，所述靶点的DNA序列如SEQ ID NO.3和SEQ ID NO.4所示。

在另一方面，本发明还提供一种获得myb84雄性不育系的方法，其特征在于，通过权利要求3-6之任一所述方法获得的 myb84雄性不育系与目标材料进行杂交和回交，从而使目标材料获得myb84雄性不育的性状和基因突变。

本发明还包括通过上述任一方法获得的myb84雄性不育系在杂交育种和制种中的应用。所述在杂交育种和制种中的应用是指将 myb84雄性不育系作为母本与其他父本进行杂交，或是将获得的myb84雄性不育系与其他目标材料进行杂交和回交，从而使目标材料获得 myb84雄性不育的性状和基因突变。

更进一步地，本发明还提供了一种玉米雄性不育系myb84的分子标记引物，引物ZmMYB84-F和ZmMYB84-R序列分别如SEQ ID NO.5和SEQ ID NO.6所示。

本发明的优点及有益效果如下：ZmMYB84（Zm00001d025664）基因及编码的蛋白调控玉米雄性生殖发育是之前没有报道过的。本发明通过利用CRISPR/Cas9方法突变玉米基因ZmMYB84（Zm00001d025664），发现了ZmMYB84（Zm00001d025664）基因对玉米雄穗发育的调控功能。利用CRISPR/Cas9基因编辑的方法和编辑后获得的雄性不育突变体，可以创制玉米雄性不育系，从而可以应用于玉米杂交育种和制种。针对三种myb84的雄性不育系开发共分离分子标记，可用于植株的育性等位基因鉴定、分子标记辅助育种中目标单株的筛选和种子纯度鉴定等。

附图说明

图1为ZmMYB84基因在玉米不同发育时期花药中的表达模式分析

S5，造孢细胞时期；S6, 小孢子母细胞时期；S7，减数分裂开始时期；S8a，减数分裂Ⅰ，二分体时期；S8b，减数分裂Ⅱ，四分体时期；S8b-9，四分体-单核小孢子时期；S9，单核小孢子时期；S9-10，单核小孢子-小孢子空泡化时期；S10，小孢子空泡化时期；S11，小孢子第一次不均等有丝分裂，二核小孢子时期；S12，小孢子第二次有丝分裂，三核小孢子时期。

图2为pCas9-ZmMYB84定点突变表达载体的物理图谱

pCas9-ZmMYB84：从T-DNA的左边界到右边界分别是除草剂抗性基因Bar的表达盒；核酸酶编码基因Cas9的表达盒；ZmMYB84基因靶标2（MT2）的表达盒；靶标1（MT1）的表达盒。

图3为野生型ZmMYB84与其不育突变体的基因结构和DNA序列分析

野生型ZmMYB84（WT- ZmMYB84）：基因全长1241 bp，包括3个外显子和2个内含子；myb84突变体ZmMYB84-Cas9-1：在第1个外显子40 bp-92 bp之间缺失53 bp；突变体ZmMYB84-Cas9-2：在第1个外显子40 bp-91 bp之间缺失52 bp；突变体ZmMYB84-Cas9-3：在第1个外显子40 bp-93 bp之间缺失54 bp。

图4为野生型与myb84纯合突变体的雄穗、花药和花粉粒表型分析

上排为玉米野生型（WT）与 ZmMYB84-Cas9-1、ZmMYB84-Cas9-2、ZmMYB84-Cas9-3突变体雄穗的表型比较；第二排为WT及与ZmMYB84-Cas9-1、ZmMYB84-Cas9-2、ZmMYB84-Cas9-3突变体花药的表型比较；下排为WT与 ZmMYB84-Cas9-1、ZmMYB84-Cas9-2、ZmMYB84-Cas9-3突变体花粉粒的I₂-KI染色比较。

图5为野生型与myb84纯合突变体的花药扫描电镜（SEM）分析

从左至右依次为：野生型（WT）花药整体；myb84花药整体；剥开后的WT（上）和myb84（下）花药；WT的成熟花粉粒（上）和myb84无法扫描到花粉粒（下）；WT（上）和myb84（下）花药的外表皮角质层；WT（上）和myb84（下）花药的内表皮乌氏体。

图6为利用共分离标记对ZmMYB84-Cas9-1不育系的F₂代植株进行基因分型

共分离标记ZmMYB84-F/R对7株ZmMYB84-Cas9-1不育系F₂代植株的PCR和琼脂糖凝胶电泳鉴定结果：在纯合野生型（AA）植株中扩增出224 bp条带；在MYB84/ myb84杂合型（Aa）植株中扩增出224 bp和171 bp两条带；在myb84/ myb84纯合突变型（aa）植株中扩增出171 bp条带。

图7为利用共分离标记对ZmMYB84-Cas9-2不育系的F₂代植株进行基因分型

共分离标记ZmMYB84-F/R对6株ZmMYB84-Cas9-2不育系F₂代植株的PCR和琼脂糖凝胶电泳鉴定结果：在纯合野生型（AA）植株中扩增出224bp条带；MYB84/ myb84杂合型（Aa）植株中扩增出224 bp和170 bp两条带；在myb84/ myb84纯合突变型（aa）植株中扩增出170 bp条带。

图8为利用共分离标记对ZmMYB84-Cas9-3不育系的F₂代植株进行基因分型

共分离标记ZmMYB84-F/R对6株ZmMYB84-Cas9-3不育系F₂代植株的PCR和琼脂糖凝胶电泳鉴定结果：在纯合野生型（AA）植株中扩增出224 bp条带；在MYB84/ myb84杂合型（Aa）植株中扩增出224 bp和172 bp两条带；在myb84/ myb84纯合突变型（aa）植株中扩增出172 bp条带。

具体实施方式

下述实施例用来说明本发明，但不限制本发明的范围。在不背离本发明精神和实质的情况下，对本发明方法、步骤或条件所作的修改或替换，均属于本发明的范围。如无特殊说明，实施例中所用引物及基因的合成和测序均由生工生物工程（上海）股份有限公司完成。其他生化试剂非特别注明外均为常规市售试剂，实施例中所用的技术手段为本领域技术人员所熟知的常规手段。

实施例一玉米ZmMYB84（Zm00001d025664）基因序列和表达模式分析

在maizeGDB 库（https://www.maizegdb.org/）中，查询到玉米ZmMYB84（Zm00001d025664，GRMZM2G173633）基因，在B73中基因的核酸序列见SEQ ID NO.1所示，基因功能注释为MYB84转录因子（MYB-transcription factor 84- MYB84)，其编码蛋白包含356个氨基酸，序列如SEQ ID NO.2所示。

MYB类转录因子参与植物中众多生理过程的调节，由于ZmMYB84基因在玉米中的实际功能还没有公开的研究资料，为了研究该基因与玉米雄性生殖发育的关系，本发明首先利用qRT-PCR分析了该基因在玉米花药发育不同阶段的表达模式。具体步骤如下：

1、玉米花药的取样与发育时期鉴定

从玉米自交系B73处于不同发育阶段的雄穗中，按照花药的长度，收集不同长度的花药样品；每份样品采集20个长度相近的新鲜花药，其中3个固定在FAA溶液中（Coolaber，中国）通过树脂半薄切片实验确定具体的发育阶段，其余17个花药立即冷冻在液氮中，用于RNA的提取。

利用梯度乙醇（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）的细胞学特点，确定每份样品的具体发育时期。

、qRT-PCR分析

用Trizol试剂（Invitrogen，美国）提取上述鉴定的处于不同发育时期（S5-S12）的玉米花药总RNA；然后使用5X All-in-One RT Master Mix （ABM，加拿大）合成cDNA；采用TBGreen™PreMix Ex Taq™（TaKaRa，日本）在QuantStudio5 QuantStudio 5 Real-Time PCRSystem （ABI，美国）上进行定量逆转录聚合酶链反应检测，扩增引物为：qMYB84-F（5’-CGGAGAAAGGCAACGACAGC-3’）和qMYB84-R（5’- CGACTTCGCACGAACGGTATT-3’）；ZmActin1为参照基因，其扩增引物为：Actin1-F（5’- AAATGACGCAGATTATGTTTGA-3’）和Actin1-R（5’-GCTCGTAGTGAGGGAGTACC-3’）；每个发育时期包括三个生物学重复，每个样品有三个技术重复；数据用2^-ΔΔCt方法进行分析，并以均值±标准差（Means±SD）的形式给出定量结果。

ZmMYB84基因呈现花药发育时期特异表达的模式：在玉米花药发育的早期和晚期几乎没有表达，而在玉米花药发育的中后期（S8b）有一个非常高的表达峰值，随后表达迅速减弱（图1）。

实施例二玉米ZmMYB84（Zm00001d025664）基因的功能以及利用CRISPR/Cas9方法创制玉米雄性不育系

为了明确玉米ZmMYB84（Zm00001d025664）在玉米中的功能，本发明采用CRISPR/Cas9基因编辑方法突变 Zm00001d025664基因序列，敲除该基因在玉米中的功能。本发明选取玉米杂交种Hi II 作为基因编辑的受体材料。本发明分别选取基因保守区的SEQ IDNO.3和SEQ ID NO.4所示序列作为CRISPR/Cas9基因编辑的靶标区域。

、ZmMYB84的CRISPR/Cas9基因编辑载体的构建

本发明的基因编辑载体为pBUE411-MT1T2-Cas9，该载体的基础载体为pBUE411- Cas9，中间载体为pCBCmT1T2，提供gRNA。本发明通过在引物上设计靶点然后通过PCR获得MT-sgRNA进而通过酶切连接到基础载体中，具体的构建流程如下：

（1）靶标gRNA的设计。将ZmMYB84（Zm00001d025664）的基因序列输入 http://crispr.hzau.edu.cn/cgi-bin/CRISPR2/CRISPR 进行靶标设计。本发明选择的两个靶标区域的DNA序列如SEQ ID NO.3和SEQ ID NO.4所示。本发明的sgRNA骨架序列从中间载体pCBCmT1T2直接扩增获得。

（2）通过在引物上设计靶点然后PCR扩增获得MT-sgRNA。引物ZmMYB84-MT1-F和引物ZmMYB84-MT2-R扩增中间载体pCBCmT1T2，用于获得包含第一个和第二个靶标的sgRNA的片段，产物长度为891 bp。PCR 体系和条件如下：模板DNA （中间载体pCBCmT1T2≥30 ng/μL）1.2 μL；Primer F/R：各1.2 μL；灭菌ddH₂O：11.4 μL；2X MCLAB 酶（产品编号：I5HM-200）：15 μL。PCR的温度程序如下：① 98 ℃ 2分钟；② 98 ℃ 10秒；③ 58 ℃ 30秒；④ 72℃ 30秒；⑤ 从②－④循环34次；⑥ 72 ℃ 5分钟；⑦ 25 ℃ 10分钟。最后回收PCR产物。载体构建所需的引物序列如下：

ZmMYB84-MT1-F: 5’-ATATATGGTCTCTGGCGAGAAGGACAGCGTCAC

AAGCGCGTTTTAGAGCTAGAAATAGCAA-3’

ZmMYB84-MT2-R: 5’-ATTATTGGTCTCTAAACTACTGGGTGATGTAGG

AGATGCTTCTTGGTGCCGC-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，灭菌ddH₂O：6 μL。

图2所示为目的基因ZmMYB84（Zm00001d025664）的双靶标（对应第一个靶标和第二个靶标），标记基因Cas9和bar与骨架载体pBUE411-Cas9构建的表达载体pCas9-ZmMYB84。

、农杆菌介导的玉米遗传转化

将上述构建的pCas9-ZmMYB84通过热激法转入农杆菌EHA105中，PCR进行鉴定；然后将含有敲除载体的农杆菌加甘油混合，于-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个月后收获后代种子。

、 T₀代植株CRISPR/Cas9突变结果检测

为确定T₀代植株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 μL75 %乙醇，洗涤沉淀两次，以10000 rpm速度离心5 min，弃掉上清；将样品放置在室温下自然干燥2-4小时，得到DNA沉淀，加入适量无菌水溶解，轻微震荡，充分溶解DNA。DNA样品存于-20 ℃保存。用Nanodrop检测DNA浓度，稀释至10 ng/L用作PCR模板。

然后根据ZmMYB84（Zm00001d025664）基因序列设计PCR引物。

检测靶标：MT1和MT2；产物大小：224 bp；引物序列如下：

ZmMYB84-T-F: 5’-TACCACTCCGAGACAGCAG-3’；

ZmMYB84-T-R: 5’-GAGGAACTACACGACACGC-3’。

提取基因组DNA，按照以下PCR参数扩增：

反应体系：15 μL MIX常规PCR体系，0.5 μL forward primer，0.5 μL reverseprimer，1 μL DNA，5.5 μL 灭菌ddH₂O，7.5 μL 2x taq mix（产品编号：10103ES ）。

反应程序：常规PCR：58 ℃退火、延伸30 s、32轮循环。

接着将PCR产物回收并连接T载体测序，通过测序多个T₀代独立阳性转化事件靶标区域的DNA序列，明确靶标区域是否发生基因编辑，最终发现3个T₀转化事件靶标区域序列发生了变化，且均为纯合突变，编辑前后的序列如图3所示，对应3个myb84纯合突变体：ZmMYB84-Cas9-1、ZmMYB84-Cas9-2、ZmMYB84-Cas9-3。与野生型序列比对显示ZmMYB84- Cas9-1、 ZmMYB84-Cas9-2和ZmMYB84-Cas9-3在靶标1和2处都发生了缺失突变。

对3个myb84纯合突变体中氨基酸序列进行比对分析发现，与未编辑的WT相比，突变后的品系ZmMYB84-Cas9-1和ZmMYB84-Cas9-2，其编码的的核苷酸在靶标1或2处的缺失使其氨基酸发生了移码突变，且随后的氨基酸发生提前终止。另外ZmMYB84-Cas9-3品系中编码的核苷酸在靶标1和2处的缺失使得发生了18个氨基酸的丢失。因此这些转化体中Zm00001d025664蛋白的功能均表现缺失。

、 F₁代植株的基因分型

由于在温室生长的玉米T₀代植株，经常雌穗和雄穗发育不协调，同时当编辑的基因与雄性发育相关时也会影响育性，因此为了繁殖T₀代植株，并使获得的基因编辑类型遗传下去，本发明使用玉米自交系郑58的野生型花粉为上述获得的ZmMYB84-Cas9-1、ZmMYB84-Cas9-2、ZmMYB84-Cas9-3的T₀代植株授粉，进而获得F₁代种子，生长的植株为F₁代植株。

F₁代植株包括2种分离类型，一种为Cas9-阳性植株（转基因植株），另一种为Cas9-阴性植株（非转基因植株），为了避免sgRNA和Cas9对杂交授粉引入的郑58野生型等位基因进行持续编辑，从而造成突变类型的复杂性，我们需要通过基因分型从F₁代植株中挑选出不含有Cas9基因，但含有T₀代突变类型的植株，这类植株自交后可以得到非转基因的F₂代。F₁代植株的基因分型步骤如下：

按照上述的CTAB法提取叶片DNA后，首先利用Cas9基因的特异引物Cas9-F（5’-CCCGGACAATAGCGATGT-3’）和Cas9-R（5’- GAGTGGGCCGACGTAGTA-3’）进行PCR扩增。PCR反应体系同上；反应程序：常规PCR：58℃退火、延伸1分钟、32轮循环。PCR产物进行琼脂糖凝胶电泳后，根据结果区分出Cas9-阳性植株和Cas9-阴性植株。

进一步针对Cas9-阴性植株，采用上述的检测MT1和MT2靶标的引物ZmMYB84-T-F和ZmMYB84-T-R进行PCR扩增；PCR产物纯化后，连接T载体，进行测序；根据测序结果分析确定T₀代突变类型的遗传情况。

实施例三myb84不育系的表型分析

上述实施例二鉴定的不含有Cas9基因的F₁代植株自交后获得F₂代种子，三种myb84突变类型（ZmMYB84-Cas9-1、ZmMYB84-Cas9-2和ZmMYB84-Cas9-3）各取1个自交单穗进行穗行播种，在成熟期进行表型的调查。三种F₂株系中，可育株与不育株的比例，均符合3:1分离，进一步表明myb84不育系的不育性状由单个隐性基因控制，然后针对F₂代获得的稳定非转基因myb84不育系与野生型进行详细的表型比较。

、雄穗、花药和花粉活力的观察

在营养生长和雌穗的发育方面， myb84不育系（ZmMYB84-Cas9-1、ZmMYB84-Cas9-2和ZmMYB84-Cas9-3）的植株与野生型相比基本无差异；在雄穗发育方面，野生型能够正常抽雄，花药能够正常开裂、散粉，自交后可正常结实，而myb84不育系虽能正常抽雄，但是不能正常开花，花药颖壳不开裂，花药明显较小，并且发白干瘪不外露（图4）；进一步对野生型和突变体的花粉进行I₂-KI染色，发现野生型花粉发育正常，花粉粒染色后呈黑色，但突变体无花粉粒形成，并且突变体的花药在长度和体积上也小于野生型（图4）。这表明ZmMYB84（Zm00001d025664）基因控制玉米的雄性发育，通过基因编辑方法创制的myb84不育系为无花粉型不育系，呈现完全败育的特点。

、花药的扫描电镜（SEM）观察

为了深入解析myb84的细胞学特征，对野生型和突变体花药内外壁进行扫描电镜（SEM）分析。剥取处于成熟期（S13）的野生型和突变体花药，然后立即固定在FAA（Coolaber，中国）溶液中，固定液的体积不少于所取研究材料体积的20倍；对于突变体花药，可用解剖针在花药壁穿孔以提高固定液的渗透效果，或者反复抽真空至花药沉入固定液底部；室温固定2小时后，将材料置于4℃保存，或依次置于50%、60%、70%、80%、90%、100%的乙醇中进行脱水，每个梯度保持15分钟；材料可置于70%的乙醇中过夜或保存。脱水后的样品进行二氧化碳临界点干燥，然后镀金即可进行观察。发现myb84突变体的花药外表皮的角质层结构相比于野生型的外表皮更加致密；但是myb84突变体花药的内表皮表现的十分光滑，没有致密的颗粒状的乌氏小体形成(图5)。花药角质层是覆盖在花药表面的细胞外脂质层，保护花药免受外部非生物胁迫、内部组织失水和病原体的侵袭，位于花药内壁的乌氏体，被认为是孢粉素前体从绒毡层细胞到小孢子的转运载体。上述结果表明ZmMYB84（Zm00001d025664）基因突变后，可以阻断绒毡层内孢子花粉素前体物质的合成，并且该基因也可能参与调节花药角质层的形成。

实施例四myb84不育系鉴定的共分离分子标记开发和应用

1、共分离分子标记的开发

在本发明中，针对获得的三种myb84不育系的突变位点，利用Primer5.0软件进行引物设计，开发出一对共分离分子标记：ZmMYB84-F/R，结合PCR与琼脂糖凝胶电泳检测的方法，根据获得的条带及大小即可分离出突变体的基因型。

共分离分子标记ZmMYB84-F/R包括第一引物ZmMYB84-F和第二引物ZmMYB84-R；该标记能够特异性检测玉米ZmMYB84-Cas9-1、ZmMYB84-Cas9-2和ZmMYB84-Cas9-3突变体及由它们转育的玉米不育材料中的突变基因myb84，并能同时区分野生型MYB84基因和突变型myb84基因；针对突变基因myb84分别扩增出171 bp、170 bp和172 bp的条带，而对野生型MYB84基因扩增出224 bp的条带。引物序列如下：

ZmMYB84-F：5’-TACCACTCCGAGACAGCAG-3’

ZmMYB84-R：5’-GAGGAACTACACGACACGC-3’

2、共分离分子标记的应用

为了验证上述标记的有效性，以实施例三中获得的F₂株系为材料，进行myb84等位基因的检测。DNA提取方法、PCR扩增体系和条件同实施例二，PCR产物进行琼脂糖凝胶电泳分离。

理论上，ZmMYB84-F/R在MYB84/ MYB84纯合野生型（AA）DNA中能扩增224 bp的条带，在myb84/myb84纯合突变型材料（aa）DNA中分别扩增出171 bp、170 bp和172 bp的条带，而在MYB84/ myb84杂合型（Aa）材料中，则能分别同时扩增出对应的两条条带。针对三种突变类型的F₂植株，ZmMYB84-F/R分子标记的验证结果如图6、图7和图8所示，结果显示设计的功能性分子标记对F₂植株的检测结果完全符合预期，在MYB84/ MYB84纯合野生型（AA）、MYB84/ myb84杂合型（Aa）和myb84/ myb84纯合突变型材料（aa）中分别扩增出对应大小的条带，可以作为MYB84，myb84等位基因检测的理想标记。

这些分子标记有助于在开花和授粉前确定突变基因型，以便在不同遗传背景下进行杂交和回交选育雄性不育系，具有重要的应用价值。

序列表

<110> 北京科技大学

北京中智生物农业国际研究院

北京首佳利华科技有限公司

<120> 雄性不育基因ZmMYB84及其在创制玉米雄性不育系中的应用

<160> 6

<170> SIPOSequenceListing 1.0

<210> 1

<211> 1241

<212> DNA

<213> 玉米(Zea mays)

<400> 1

atggggcgga tcccgtgctg cgagaaggac agcgtcaagc gcgggcagtg gacgcccgag 60

gaggacaaca agctgctctc ctacatcacc cagtacggca cgcgcaactg gcgcctcatc 120

cccaagaatg ccggtacgtc ctagcgaccg tccaatgtcc gcgtgtcgtg tagttcctct 180

caagtgtgac cgagtccgcg tgatgggatt gcaggactgc agcgatgcgg gaagagctgc 240

cggctccggt ggaccaacta cctgcgtccc gacctcaagc acggtgagtt cacggacacc 300

gaggagcaga ccatcatcaa gctgcactcc gtcgttggca acaggtggcc cgcccgtcct 360

cccctctcgt cgtgcctctt ctgatgggtg cctgttctga ctacaagcaa tgtgcaacaa 420

acgcgcgcgc gtgcaggtgg tcggtgatcg cggcgcagct gccgggtcgg acggacaacg 480

acgtcaagaa ccactggaac accaagctga agaagaagct gtccgggatg ggcatcgacc 540

ccatcacgca caagtccttc tcgcacctca tggccgagat cgccaccacg ctggcgccgc 600

cgcaggtggc ccacctcgcc gaggccgcgc tggggtgctt caaggacgag atgctccacc 660

tcctcaccaa gaagcgcccc accgacttcc cgtcgcccgc ggtgcccgac atgtcggcga 720

tcgcgggcgg ctccggcgtc gcggcgccct gcggcttccc ggcgccgccc cagaccgacg 780

acaccatcga gcgcatcaag ctgggcctgt cccgcgccat catgagcgag cccgccgcgc 840

cccccggcaa gcaggagcag ccctgggcgc cggccgactt gccggagggg ctgccgggga 900

tgtacgccac gtacaatccc gcctcgcacg gacacgaaga gttccgctac gacaacggga 960

cagtgccgga gtacgtcctc ggcggcggcg gcggcgcgga ccagggcacg tcgatgtgga 1020

gccaccagag catgtacagc gggagttcgg ccacggaggc cgcgcccagg ccggcggagg 1080

tgttgccgga gaaaggcaac gacagcgtcg ggagcagcgg cggcggcgag gaggcggacg 1140

acgtcaagga cggcgggaaa ggcggctccg atatgtccgg cctgtttgga tccgactgcg 1200

tactttggga cttgcccgac gagctgacca atcacatggt g 1241

<210> 2

<211> 356

<212> PRT

<213> 玉米(Zea mays)

<400> 2

Met Gly Arg Ile Pro Cys Cys Glu Lys Asp Ser Val Lys Arg Gly Gln

1 5 10 15

Trp Thr Pro Glu Glu Asp Asn Lys Leu Leu Ser Tyr Ile Thr Gln Tyr

20 25 30

Gly Thr Arg Asn Trp Arg Leu Ile Pro Lys Asn Ala Gly Leu Gln Arg

35 40 45

Cys Gly Lys Ser Cys Arg Leu Arg Trp Thr Asn Tyr Leu Arg Pro Asp

50 55 60

Leu Lys His Gly Glu Phe Thr Asp Thr Glu Glu Gln Thr Ile Ile Lys

65 70 75 80

Leu His Ser Val Val Gly Asn Arg Trp Ser Val Ile Ala Ala Gln Leu

85 90 95

Pro Gly Arg Thr Asp Asn Asp Val Lys Asn His Trp Asn Thr Lys Leu

100 105 110

Lys Lys Lys Leu Ser Gly Met Gly Ile Asp Pro Ile Thr His Lys Ser

115 120 125

Phe Ser His Leu Met Ala Glu Ile Ala Thr Thr Leu Ala Pro Pro Gln

130 135 140

Val Ala His Leu Ala Glu Ala Ala Leu Gly Cys Phe Lys Asp Glu Met

145 150 155 160

Leu His Leu Leu Thr Lys Lys Arg Pro Thr Asp Phe Pro Ser Pro Ala

165 170 175

Val Pro Asp Met Ser Ala Ile Ala Gly Gly Ser Gly Val Ala Ala Pro

180 185 190

Cys Gly Phe Pro Ala Pro Pro Gln Thr Asp Asp Thr Ile Glu Arg Ile

195 200 205

Lys Leu Gly Leu Ser Arg Ala Ile Met Ser Glu Pro Ala Ala Pro Pro

210 215 220

Gly Lys Gln Glu Gln Pro Trp Ala Pro Ala Asp Leu Pro Glu Gly Leu

225 230 235 240

Pro Gly Met Tyr Ala Thr Tyr Asn Pro Ala Ser His Gly His Glu Glu

245 250 255

Phe Arg Tyr Asp Asn Gly Thr Val Pro Glu Tyr Val Leu Gly Gly Gly

260 265 270

Gly Gly Ala Asp Gln Gly Thr Ser Met Trp Ser His Gln Ser Met Tyr

275 280 285

Ser Gly Ser Ser Ala Thr Glu Ala Ala Pro Arg Pro Ala Glu Val Leu

290 295 300

Pro Glu Lys Gly Asn Asp Ser Val Gly Ser Ser Gly Gly Gly Glu Glu

305 310 315 320

Ala Asp Asp Val Lys Asp Gly Gly Lys Gly Gly Ser Asp Met Ser Gly

325 330 335

Leu Phe Gly Ser Asp Cys Val Leu Trp Asp Leu Pro Asp Glu Leu Thr

340 345 350

Asn His Met Val

355

<210> 3

<211> 19

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 3

gaaggacagc gtcaagcgc 19

<210> 4

<211> 19

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 4

tctcctacat cacccagta 19

<210> 5

<211> 19

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 5

taccactccg agacagcag 19

<210> 6

<211> 19

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 6

gaggaactac acgacacgc 19

Claims

1.ZmMYB84基因在控制玉米雄性生殖发育中的应用，其特征在于，所述基因的氨基酸序列如SEQ ID NO.2所示。

2.根据权利要求1所述的应用，其特征在于，所述基因的核苷酸序列为SEQ ID NO.1所示。

3.一种创制玉米雄性不育系的方法，其特征在于，抑制玉米中权利要求1或2所述基因的表达和/或活性，选择玉米雄性不育的植株。

4.根据权利要求3所述创制玉米雄性不育系的方法，其特征在于，所述抑制基因表达和/或活性的方法包括基因编辑、RNA干扰、T-DNA插入中任一种。

5.根据权利要求4所述创制玉米雄性不育系的方法，其特征在于，所述基因编辑采用CRISPR/Cas9方法。

6.根据权利要求5所述创制玉米雄性不育系的方法，其特征在于，所述CRISPR/Cas9方法包括：在所述基因第一外显子处设计CRISPR/Cas9载体靶点，所述靶点的DNA序列如SEQID NO.3和SEQ ID NO.4所示。

7.一种获得myb84雄性不育系的方法，其特征在于，通过权利要求3-6之任一所述方法获得的 myb84 雄性不育系与目标材料进行杂交和回交，从而使目标材料获得myb84雄性不育的性状和基因突变。

8.通过权利要求3-6之任一所述的方法获得的 myb84不育系在杂交育种和制种中的应用。

9.根据权利要求8所述的应用，其中所述的杂交育种和制种是指将 myb84不育系作为母本与其他父本进行杂交。

10.根据权利要求8所述的应用，包括将获得的myb84雄性不育系与其他目标材料进行杂交和回交，从而使目标材料获得 myb84雄性不育的性状和基因突变。

11.一种玉米雄性不育系myb84的分子标记引物ZmMYB84-F 和ZmMYB84-R，其特征在于，引物ZmMYB84-F和ZmMYB84-R序列分别如SEQ ID NO.5和SEQ ID NO.6所示。