qMYR2基因在调节或筛选稻米中油脂组分含量的应用的制作方法

本发明涉及水稻分子育种领域，更特别地，涉及qmyr2基因在调节和筛选稻米中油脂组分含量的应用。

背景技术：

植物种子中脂质主要以甘油三酯的形式储存。大米中的油脂对大米的储存、食用品质和消费者的健康都非常重要，其含油率一般为2-4％。米饭的食用品质和适口性受到大米中的脂肪酸组成与含量的影响。目前，本领域中对大米的脂质合成通路知之甚少，liu等证实了osfad3将亚油酸(c18:2)转化成亚麻酸(c18:2)。zaplin等证实osfad2-1将油酸转化成亚油酸。一些研究者还初步鉴定了几个数量性状位点(qtl)。

但是，目前鉴定到的基因和qtl都太少，仍然无法知晓大米油脂合成的遗传和分子基础，因此需要找到影响大米中的油脂组分及含量的主效基因，用做筛选标志或遗传操作的靶点。

技术实现要素：

在本发明的发明人通过研究发现，鉴定了23个与油含量和油脂组成显著相关的位点，并且通过连锁分析的验证得到了11个有可能参与到了水稻中与稻米油脂代谢有关的通路的候选基因。在进一步的研究后，我们发现这些位点中，4个基因对油脂组成的自然变异贡献巨大，并且在亚群中有分化。其中之一就是qmyr2基因，我们的研究显示，qmyr2是正向调控c14:0含量的主效基因。该基因的突变如果导致蛋白功能的降低或消失，使稻米中的c14:0含量降低。

基于以上研究，本发明提供了qmyr2基因在调节稻米中油脂组分含量的应用，所述qmyr2基因的序列如seqidno:1所示。

在一个优选实施方案中，所述油脂组分为c14:0。

本发明还提供了一种调节稻米中油脂组分含量的方法，包括导入qmyr2基因、提高qmyr2基因表达水平、降低qmyr2基因表达水平或者突变qmyr2基因使qmyr2基因表达的蛋白功能提高、降低或消失的步骤，所述qmyr2基因的序列如seqidno:1所示。

例如，如果目标稻种其他性状优良，但是稻米中c14:0含量太高，并且该稻种含有正常功能的qmyr2基因，可通过敲除或突变目标稻种中的qmyr2基因，使qmyr2基因不表达，或者表达的蛋白的肉豆蔻酰基-acp硫酯酶功能降低或没有肉豆蔻酰基-acp硫酯酶功能，即可降低稻米中的c14:0含量。

反之，如果目标稻种其他性状优良，但是稻米中c14:0含量太低，并且该稻种不含有正常功能的qmyr2基因，可通过导入qmyr2基因，或突变目标稻种中的不正常的qmyr2基因使其回复突变成具有正常肉豆蔻酰基-acp硫酯酶功能的qmyr2基因，即可提高稻米中的c14:0含量。

本发明还提供了qmyr2基因在水稻育种中作为稻米油脂组分含量的筛选标志的应用，所述qmyr2基因的序列如seqidno:1所示。

本发明还提供了一种水稻杂交育种方法，包括使用qmyr2基因或其突变体作为筛选标志来选择目标基因型的步骤。

在一个优选实施方案中，所述水稻杂交育种方法包括以下步骤：

s1：鉴定亲本中的qmyr2基因及其突变状态；

s2：选择qmyr2基因或其突变体作为筛选标志；

s3：将所述筛选标志与其他优势性状聚合到一起。

例如，已有稻种的其他性状优良，但稻米中c14:0含量太高，并且含有正常功能的qmyr2基因，要想通过杂交育种的方法降低稻米的c14:0含量，可先选择qmyr2基因发生了导致磷脂酰胆碱：甘油二脂转磷酸胆碱酶2功能降低或消失的突变的亲本与该稻种进行杂交，通过将qmyr2基因突变体作为分子标志来使该突变基因与其他优势性状聚合到一起。

反之，已有稻种的其他性状优良，但稻米中c14:0含量太低，并且不含正常功能的qmyr2基因，要想通过杂交育种的方法提高稻米的c14:0含量，可先选择具有正常肉豆蔻酰基-acp硫酯酶功能的qmyr2基因亲本与该稻种进行杂交，通过将该正常qmyr2基因作为分子标志来使该突变基因与其他优势性状聚合到一起。

附图说明

图1为通过gc-ms检测的533种稻米品种中的10种油脂组分的百分含量；

图2为533个稻米品种中澳洲稻(aus)、籼稻(ind)、粳稻(jap)的油脂相关性状的表型分布；

图3为籼稻亚群中c14:0含量关联分析的曼哈顿图(manhattanplot)，虚线表示显著性阈值(-log10(p)＝7.06)

图4为染色体2上的显著性峰周围的区域曼哈顿图和连锁不平衡(ld)热图，箭头表示候选基因的核苷酸变异位点，垂向虚线表示峰周围的候选区；

图5为qmyr2的基因结构以及该基因的多态性；

图6为通过连锁分析鉴定qmyr2作为88b2/hd9802s的重组自交(ril)群体中c14：0含量的主效位点的qtl扫描结果的统计图；

图7为不同qmyr2单倍型的c18:2含量统计图，其中a的测试对象为获取群体中的305个籼稻品种，b的测试对象为由88b2/hd9802s衍生的ril群体；

图8为野生型和敲除突变家系中的10种脂肪酸含量的统计图，*p<0.05,**p<0.01；***p<0.001。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。

1、样品来源

本发明使用的水稻材料包括533个栽培稻品种构成的自然群体和三个自己培育的重组自交群体。自然群体用于全基因组关联分析，三个重组自交群体用于连锁分析。

其中，533个栽培稻品种包括50个aus品种(aus)、305个籼稻品种(indica)和178个粳稻品种(japonica)。

三个重组近繁品种分别衍生自以下杂交系：珍汕97/明恢63(zs97/mh63)、88b-2/hd9802s和b805d/h6s。

2、样品处理

将稻谷脱壳成糙米，然后碾粉，过80目筛，于80℃下干燥24h，然后进行脂质提取。0.2g米粉加入10ml玻璃管中，加入4.5ml硫酸:甲醇溶液(1:19)和100μl溶于氯仿的十七烷酸标准溶液，混匀。十七烷酸标准溶液的配置为1.07g十七烷酸标样溶于10ml色谱纯氯仿。将混合物在85℃水浴中孵育2h。冷却至室温后，使用正己烷提取脂肪酸甲酯，并通过气相色谱-质谱联用仪(gc-ms)分析。所有鉴定的脂肪酸含量的总和作为含油量。每个样品进行两次测试，将平均值用于后续分析。

3、gc-ms分析结果

gc-ms分析结果如图1所示，鉴定出的10种脂肪酸中，软脂酸(c16:0)、油酸(c18:1)和亚油酸(c18:2)占总脂肪酸含量的90％以上。我们在研究中发现了大量油脂相关性状存在广泛的变异，从c18:2含量的1.5倍差异到肉豆蔻酸(c14:0)含量的7.4倍差异。总脂肪酸含量介于22.99-40.96mg/g的范围。

其中，澳洲稻(aus)中的c16:0含量最高，其次是籼稻(ind)和粳稻(jap)，c18:1和c18:2的含量则相反(图2)。

4、关联分析和连锁分析

全基因组关联研究(gwas)：分别对籼稻、粳稻和所有样品实施gwas。在每个分析组中，仅选择次要位点频率(ma)f＞5％、缺失率＜15％的snp进行关联分析。

连锁分析：对三个ril群体——zs97/mh63(96个株系)、88b-2/hd9802s(108个株系)和b805d/h6s(135个株系)——进行连锁分析。这些株系均通过二代测序，然后将50bp双末端读序比对到水稻参考基因组。其中比对质量≥40的序列和平均质量≥10的碱基用于鉴定snp。通过使用mpr算法基因分型来推断亲本株系的基因型，构建遗传bin图谱用于鉴定qtl。

根据以上工作，我们最终获得了520万个高质量snp，用于评估群体结构和亲缘关系，并在其中鉴定了23显著的关联位点。在23个位点中，我们发现了qmyr2是与c14:02含量相关的主效qtl(图3)。

如图4所示，qmyr2位于2号染色体的长臂上，我们基于点配对连锁不平衡(ld)相关系数(r²＜0.6)分析了25.86-26.02mb(162kb)区域内的序列，最终确定该qtl对应的候选基因为loc_os02g43090(序列为seqidno:1)，编码肉豆蔻酰基-acp硫酯酶其在拟南芥中的同源基因为fatb(酯酰基-acp硫酯酶b)基因。

loc_os02g43090基因1号外显子上第169位的g-a突变(ex1-g/a)存在较大的次等位基因频率，该突变导致氨基酸改变(图5)。

qmyr2的ex1-g/a位点的基因型在籼稻品种88b-2与hd9802s之间存在差异。并且，使用88b-2/hd9802s杂交系衍生的ril群体进行的qtl分析证实了qmyr2单倍型对c14:0含量的对比影响qmyr2(图6，表1)。

表1油脂相关性状的表型分析

根据突变位点的a和t碱基，我们将qmyr2划分为单倍型a(seqidno:1)和单倍型b(seqidno:2)。无论在305个籼稻品种群体中还是ril群体中，单倍型a株系的c14:0含量均高于单倍型b株系(图7)。

5、转基因分析

为了进一步确定qmyr2的功能，我们利用crispr-cas9系统敲除了zs97中的qmyr2。在制备crispr/cas9敲除构建体上，我们仔细地检查了基因编码区，设计了带有c末端ngg的23bp靶点。确定靶点后，我们将20bp靶点插入中间载体per8-cas9-u6或per8-cas9-u3，并克隆至pcxun-cas9，得到敲除载体，命名为ko-myr2。将敲除构建体导入大肠杆菌straintrans5α中进行测序，验证正确后，导入根瘤农杆菌agrobacteriumtumefaciensstraineha105，然后转入zs97。获得转基因植株后，通过测序验证敲除系的基因型。

与野生型zs97(单倍型a)相比，zs97敲除系ko-qmyr2的c18:1含量显著升高，c14:0含量显著降低(图8)。证明该位点可用于调节大米中的脂肪酸组成。

综合上面的实验可知，qmyr2是正调控c14:0含量的主效基因。该基因的突变如果导致蛋白功能的降低或消失，使稻米中的c14:0含量降低。

虽然以上的实例中仅展示了通过crispr-cas9系统得到敲除突变体和ex1-g/a突变体，但是，本领域技术人员在阅读本申请披露的内容后可知，只需要将该基因导入、敲除、敲降或突变成相关功能降低或无相关功能，即可达到调节稻米中c14:0含量的目的。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

序列表

<110>华中农业大学

湖北双水双绿生物科技有限公司

<120>qmyr2基因在调节或筛选稻米中油脂组分含量的应用

<160>2

<170>siposequencelisting1.0

<210>1

<211>897

<212>dna

<213>水稻(oryzasativa)

<400>1

atggcagggtctcttgccgcctcagcattcttcccaggtccaggctcatctcctgcagca60

tcagctagaagctccaagaatgctgctgttaccggcgaattgccggagaatttgagtgtc120

tgtggcattgtcgcaaagcctaacccacctcctgcagccatgcaagtaaaggcacaggct180

caaacccttcccaaggttaatggtacgaaggttaacctcaagacggtgaagcctgacatg240

gaggaaacggtgcctcacagtgctccaaagacgttctataaccaactgccggattggagc300

atgcttcttgcggctattacaaccatcttcctcgccgcagagaagcagtggacactgctt360

gattggaagccgaagaaacctgacatgcttgttgacacatttggctttggtaggatcatc420

caggacggtatggtgtttaggcagaacttcatgattcggtcctacgagattggcgctgat480

cgtacagcttctatagagacattgatgaatcatttacaggaaacggctcttaaccatgta540

aggactgctggtcttcttggagatggttttggggctacaccggagatgagcaaacggaac600

ttgatatgggttgtcagcaaaatccagcttcttgttgagcaataccccgcatggggagat660

atggttcaagttgacacatgggtcgctgctgctggcaaaaatggcatgcgtcgagactgg720

catgttcgtgactacaactctggccgaacaatcttgagagctacaagtgtttgggtgatg780

atgcacaagaaaactagaagactttcaaaaatgccagatgaagttagagctgaaataggc840

ccatatttcaatgaccgttcagctataacagaggagcagagtgaaaagttagcctag897

<210>2

<211>1281

<212>dna

<213>水稻(oryzasativa)

<400>2

atggcagggtctcttgctgcctcagcattcttcccaggtccaggctcatctcctgcagca60

tcagctagaagctccaagaatgctgctgttaccggcgaattgccggagaatttgagtgtc120

tgtggcattgtcgcaaagcctaacccacctcctgcagccatgcaagtaaaggcacaggct180

caaacccttcccaaggttaatggtacgaaggttaacctcaagacggtgaagcctgacatg240

gaggaaacggtgcctcacagtgctccaaagacgttctataaccaactgccggattggagc300

atgcttcttgcggctattacaaccatcttcctcgccgcagagaagcagtggacactgctt360

gattggaagccgaagaaacctgacatgcttgttgacacatttggctttggtaggattatc420

caggacggtatggtgtttaggcagaacttcatgatttggtcctacgagattggcgctgat480

cgtacagcttctatagagacattgatgaatcatttacaggaaacggctcttaaccatgta540

aggactgctggtcttcttggagatggttttggggctacaccggagatgagcaaacggaac600

ttgatatgggttgtcagcaaaatccagcttcttgttgagcaataccccgcatggggagat660

atggttcaagttgacacatgggtcgctgctgctggcaaaaatggcatgcgtcgagactgg720

catgttcgtgactacaactctggccgaacaatcttgagagctacaagtgtttgggtgatg780

atgcacaagaaaactagaagactttcaaaaatgccagatgaagttagagctgaaataggc840

ccatatttcaatgaccgttcagctataacagaggagcagagtgaaaagttagccaagaca900

ggaaataaagttggtgatgatgctacagagcaattcataagaaaggggctcactcctaga960

tggggtgacctcgatgtcaatcagcatgtgaacaatgttaaatatattgggtggatcctt1020

gaggtgagtgctccaatttcagtactggagaagcatgagcttgcaagcatgaccctggat1080

tacaggaaggagtgtggccgagacagcgtgctgcaatcacttaccaccgtgtcaggggaa1140

tgcaccagcattggcgccgacaagcaggcttctgccatccagtgcgaccatcttcttcag1200

cttgaatcaggagctgacattgtgaaggcacacacagagtggcgaccaaagcgatcgcat1260

gcagcagctgagaacgcgtaa1281