• ISSN 1008-505X
  • CN 11-3996/S

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

正常供磷及低磷胁迫下油菜种子植酸浓度与含量的全基因组关联分析

刘海疆 李晓娟 谢毅文 吴海亚 鲁明星 丁广大 徐芳森 石磊

引用本文:
Citation:

正常供磷及低磷胁迫下油菜种子植酸浓度与含量的全基因组关联分析

    作者简介: 刘海疆 E-mail: liuhaijiang@webmail.hzau.edu.cn;
    通讯作者: 石磊, E-mail:leish@mail.hzau.edu.cn
  • 基金项目: 国家自然科学基金项目(32172662,31972498);中央高校基本科研业务费专项基金资助项目(2662019PY013)。

Genome wide association study of seed phytate concentration and content in Brassica napus under contrasting phosphate conditions

    Corresponding author: SHI Lei, E-mail:leish@mail.hzau.edu.cn
  • 摘要:   【目的】  鉴定不同磷水平油菜种子植酸浓度 (PA_Conc)和含量 (PA_Cont) 显著关联的 SNP 位点,挖掘候选基因及其优异单倍型,为油菜低植酸品种的培育提供遗传基础和优异基因资源。  【方法】  采用田间试验,测定正常供磷(施磷,P2O5 90 kg/hm2)及低磷(不施磷,P2O5 0 kg/hm2)条件下403个甘蓝型油菜品种种子中的PA_Conc和PA_Cont,利用基于全基因组重测序技术获得的530多万个高质量SNP (single nucleotide polymorphism)标记,用 TASSEL 5.0 软件中的一般线性模型(general linear model,GLM)和混合线性模型(mixed linear model,MLM)进行全基因组关联分析(genome-wide association study, GWAS),鉴定种子植酸浓度和含量显著关联的SNP位点及其候选基因。采用Bcftools提取BnaA07.MRP13基因的SNP,利用Haploview 4.2软件进行单倍型分析。  【结果】  关联群体正常供磷处理PA_Conc变异范围为17.52~81.32 mg/g,PA_Cont变异范围为0.30~0.99 mg/5粒种子;低磷处理PA_Conc变异范围为17.47~49.05 mg/g,PA_Cont变异范围为0.33~0.85 mg/5粒种子。与正常供磷相比,关联群体低磷胁迫下的PA_Conc和PA_Cont均值分别降低13.5%和10.8%。GLM模型检测到56个与油菜种子植酸浓度和含量显著关联的SNP位点,分布在除A08和C03外的17条染色体上,它们对表型变异的解释率为7.33%~14.28%。其中A03染色体上含有最多的SNP (9个)。MLM模型共检测到23个与油菜种子植酸浓度和含量显著关联的SNP位点,分布在A03、A05、A06、A07、A09、C01、C02、C04、C05、C06和C07等11条染色体上,它们对表型变异的解释率为9.6%~14.28%。其中A03染色体上含有8个SNP。对显著SNP位点上下游300 kb范围内的基因进行分析,挖掘出可能与种子中植酸合成相关的14个候选基因,其中包括拟南芥植酸合成通路中的MRPPLCPIP5KSULTR4等同源基因。对BnaA07.MRP13进行候选基因关联及单倍型分析,鉴定出高植酸单倍型“BnaA07.MRP13ContHap1”(TTTA)和低植酸单倍型“BnaA07.MRP13ContHap2”(CCCT)。自然群体中18%的品种(72/403)为低植酸单倍型,种子平均植酸含量为 0.56 mg/5粒种子,极显著低于高植酸单倍型 (平均植酸含量为0.62 mg/5粒种子)。  【结论】  正常供磷和低磷胁迫下,油菜种子植酸浓度和含量均具有广泛的遗传变异。GLM模型和MLM模型分别检测到56和23个与油菜种子植酸浓度和含量显著关联的SNP位点,预测的与油菜种子植酸含量相关的候选基因有14个,发现的BnaA07.MRP13的高和低植酸单倍型为后续油菜种子植酸含量的遗传改良奠定了基础。
  • 图 1  正常磷处理和低磷处理甘蓝型油菜关联分析群体种子植酸浓度和含量的频率分布直方图

    Figure 1.  Frequency distribution of phytate concentration and content of an association panel of Brassica napus under sufficient and deficient phosphorus supplies

    图 2  正常磷和低磷处理油菜种子植酸浓度和含量全基因组关联分析的曼哈顿图和QQ图

    Figure 2.  Manhattan and QQ plots of GWAS on seed phytate concentration and content of Brassica napus under sufficient and deficient P supply

    图 3  BnaA07.MRP13的候选基因关联分析(a)和单倍型分析(b)

    Figure 3.  Candidate genes association analysis (a) and haplotype analysis (b) of BnaA07.MRP13

    表 1  正常磷处理及低磷处理甘蓝型油菜关联分析群体种子植酸浓度和含量

    Table 1.  Seed phytate concentration and content in an association panel of Brassica napus under sufficient and deficient phosphorus supplies

    性状 Trait最大值 Max最小值 Min均值 Mean峰度 Kurtosis偏度 Skewness变异系数 CV (%)
    正常磷种子植酸浓度 PA_Conc_SP (mg/g)81.3217.5236.871.230.6721.37
    低磷种子植酸浓度 PA_Conc_DP (mg/g)49.0517.4731.91–0.440.3824.57
    正常磷种子植酸含量 PA_Cont_SP (mg/5 seeds)0.990.300.65–0.290.1818.85
    低磷种子植酸含量 PA_Cont_DP (mg/5 seeds)0.850.330.58–0.03–0.4018.91
    注(Note):PA—植酸 Phytate; Conc—浓度 Concentration; Cont—含量 Content; SP—充足磷 Sufficient P supply; DP—低磷 Deficient P supply.
    下载: 导出CSV

    表 2  正常磷和低磷处理甘蓝型油菜种子植酸浓度和植酸含量显著关联的SNP位点

    Table 2.  Significant SNPs associated with phytate concentration and content in Brassica napus seed under sufficient and deficient phosphorus supplies

    性状 Trait显著关联 SNP Significant SNPPP-value表型解释方差 PVE (%)模型 Model
    低磷植酸浓度
    PA_Conc_DP
    chrA03__273392794.34×10–710.72GLM
    chrA03__273394015.56×10–914.28GLM、MLM
    chrA10__127752942.38×10–712.36GLM
    chrAnn_random__242595924.24×10–711.60GLM
    chrAnn_random__242596053.63×10–711.52GLM
    chrCnn_random__780804851.19×10–79.36GLM
    低磷植酸含量
    PA_Cont_DP
    chrA02__214547667.09×10–79.06GLM
    chrA02__225013374.46×10–711.72GLM
    chrA02__242071274.62×10–712.10GLM
    chrA07__78092971.20×10–711.22GLM、MLM
    chrA07__78093001.20×10–711.18GLM、MLM
    chrA07__78093162.65×10–710.88GLM
    chrUnn_random__62233776.45×10–78.37GLM
    正常磷植酸浓度
    PA_Conc_SP
    chrA01__118142605.36×10–79.50GLM
    chrA02__241708813.53×10–79.60GLM
    chrA03__170574452.73×10–912.70GLM、MLM
    chrA03__191645673.12×10–913.05GLM、MLM
    chrA03__191646571.14×10–913.86GLM、MLM
    chrA03__255851664.12×10–710.07GLM、MLM
    chrA03_random__12899641.29×10–913.89GLM、MLM
    chrA03_random__12911852.46×10–811.56GLM、MLM
    chrA03_random__12912071.71×10–812.04GLM、MLM
    chrA04__44209445.24×10–710.01GLM
    chrA05__48486977.38×10–811.02GLM、MLM
    chrA06__163346619.56×10–912.35GLM、MLM
    chrA06__163346677.38×10–912.56GLM、MLM
    chrA09__62094155.27×10–810.01GLM、MLM
    chrA10__40489105.80×10–79.43GLM
    chrA10__40489196.01×10–79.40GLM
    chrA10__40489304.83×10–79.54GLM
    chrC01__74735651.35×10–811.66GLM、MLM
    chrC02__307762427.91×10–811.06GLM、MLM
    chrC02__316941602.42×10–79.60GLM、MLM
    chrC04__371174574.70×10–78.92GLM
    chrC04__434268206.38×10–710.23GLM
    chrC04__442705494.62×10–810.82GLM、MLM
    chrC05__4340777.13×10–811.07GLM、MLM
    chrC06__290649443.25×10–810.67GLM、MLM
    chrC06__306154366.57×10–79.03GLM
    chrC06__322088457.19×10–79.30GLM
    下载: 导出CSV
    续表 2 Table 2 continued
    性状 Trait显著关联 SNP Significant SNPPP-value表型解释方差 PVE (%)模型 Model
    chrC06__354253063.63×10–811.00GLM、MLM
    chrC06__354253126.56×10–810.71GLM、MLM
    chrC06_random__30871893.65×10–78.97GLM
    chrC07__348576545.82×10–79.73GLM
    chrC07__349563611.74×10–710.32GLM、MLM
    chrC07__371987295.14×10–79.14GLM
    chrC08__314146545.27×10–79.53GLM
    chrC08_random__44363911.67×10–812.11GLM
    chrC09__30881134.45×10–710.11GLM
    chrCnn_random__105359003.69×10–78.90GLM
    chrCnn_random__167746242.10×10–811.47GLM
    chrCnn_random__297580175.64×10–78.92GLM
    chrCnn_random__322494451.58×10–79.16GLM
    chrCnn_random__331278745.16×10–79.94GLM
    chrCnn_random__731820295.79×10–79.12GLM
    正常磷植酸含量
    PA_Cont_SP
    chrC09__373436224.71×10–77.33GLM
    注(Note):PA—植酸 Phytate; Conc—浓度 Concentration; Cont—含量Content; SP—充足磷 Sufficient P supply; DP—低磷 Deficient P supply; PVE—表型解释方差 Phenotypic variance explained.
    下载: 导出CSV

    表 3  正常磷和低磷处理甘蓝型油菜种子植酸浓度和含量显著关联 SNP 的候选基因

    Table 3.  Candidate genes underlying the significant SNPs for seed phytate concentration and content of Brassica napus under sufficient and deficient phosphorus supply conditions

    性状Trait显著 SNP
    Significant SNP
    油菜基因
    Genes in Brassica napus
    模型
    Model
    基因功能注释
    Gene function annotation
    正常磷植酸含量
    PA_Cont_SP
    chrC06__35245326BnaC06g37190DGLM、MLM磷酸盐转运蛋白1;9
    Phosphate transporter 1;9 (PHT1;9)
    chrC09__37343622BnaC09g33980DGLM、MLM磷脂酰肌醇特异性磷脂酶C1
    Phospholipase C1 (PLC1)
    chrC09__37343622BnaC09g34000DGLM、MLM磷脂酰肌醇特异性磷脂酶C5
    Phospholipase C5 (PLC5)
    chrC09__37343622BnaC09g34010DGLM、MLM磷脂酰肌醇特异性磷脂酶C4
    Phospholipase C4 (PLC4)
    低磷植酸含量
    PA_Cont_DP
    chrC04__9656140BnaC04g24700DGLM磷脂酰肌醇单磷酸5-激酶-4
    Phosphatidylinositol mono phosphate5-kinase-4 (PIP5K4)
    chrC07__34857654BnaC07g30200DGLM磷脂酰肌醇-3, 4-激酶
    Phosphatidylinositol -3- and 4-kinase
    chrC07__34956361BnaC07g30170DGLM液泡膜铁转运家族蛋白
    Vacuolar iron transporter family protein (VIT)
    chrA03__19164557BnaA03g04410DMLM硫酸盐转运蛋白4;1
    Sulfate transporter 4;1 (SULTR4;1)
    chrA07__7809297BnaA07g07310DGLM、MLM多药耐药相关蛋白13
    Multidrug resistance-associated protein 13 (MRP13)
    chrA07__7809297BnaA07g07320DGLM、MLM多药耐药相关蛋白13
    Multidrug resistance-associated protein 13 (MRP13)
    chrA07__7809297BnaA07g07330DGLM、MLM多药耐药相关蛋白12
    Multidrug resistance-associated protein 12 (MRP12)
    chrA02__24207127BnaA02g34000DGLM、MLM3'(2'),5'-二磷酸核苷酸酶活性
    3'(2'),5'-bisphosphate nucleotidase activity
    chrA02__24207127BnaUnng04100DGLM、MLM磷酸甘油酸变位酶家族蛋白,参与InsP6到InsP7的转换和InsP7到InsP8的转换
    Phosphoglycerate mutase-like family protein
    低磷植酸浓度
    PA_Conc_DP
    chrA10__12775294BnaA10g17710DGLM、MLM肌醇1,3,4磷酸-5/6-激酶蛋白
    Inositol 1,3,4-trisphosphate -5/6-kinase family protein
    注(Note):PA—植酸 Phytate; Conc—浓度 Concentration; Cont—含量 Content; SP—充足磷 Sufficient P supply; DP—低磷 Deficient P supply.
    下载: 导出CSV
  • [1] 中华人民共和国国家统计局. 中国统计年鉴[M]. 北京: 中国统计出版社, 2020.

    National Bureau of Statistics of China. China statistical yearbook[M]. Beijing: China Statistics Press, 2020.
    [2] 王汉中. 以新需求为导向的油菜产业发展战略[J]. 中国油料作物学报, 2018, 40(5): 613–617. Wang H Z. New-demand oriented oilseed rape industry developing strategy[J]. Chinese Journal of Oil Crop Sciences, 2018, 40(5): 613–617.
    [3] Kochian L V. Plant nutrition: Rooting for more P[J]. Nature, 2012, 488(7412): 466–467.
    [4] Lickfett T, Mattha B, Velasco L, et al. Seed yield, oil and phytate concentration in the seeds of two oilseed rape cultivars as affected by different phosphorus supply[J]. European Journal of Agronomy, 1999, 11(3–4): 293–299.
    [5] 张倩雯, 丁广大, 王效华, 等. 植物种子植酸研究进展[J]. 植物科学学报, 2016, 34(5): 814–820. Zhang Q W, Ding G D, Wang X H, et al. Research progress on plant seed phytate[J]. Plant Science Journal, 2016, 34(5): 814–820.
    [6] 李晓娟. 甘蓝型油菜高和低植酸品种种子植酸积累的差异及植酸含量相关候选基因的挖掘[D]. 武汉: 华中农业大学硕士学位论文, 2020.

    Li X J. Differences in phytate accumulation in seeds between high and low phytate Brassica napus cultivars and identification of candidate genes for phytate synthesis[D]. Wuhan: MS Dissertation of Huazhong Agricultural University, 2020.
    [7] Raboy V. Seeds for a better future: ‘low phytate’ grains help to overcome malnutrition and reduce pollution[J]. Trends in Plant Science, 2001, 6(10): 458–462.
    [8] Raboy V. Approaches and challenges to engineering seed phytate and total phosphorus[J]. Plant Science, 2009, 177(4): 281–296.
    [9] Otegui M S, Capp R, Staehelin L A. Developing seeds of Arabidopsis store different minerals in two types of vacuoles and in the endoplasmic reticulum[J]. Plant Cell, 2002, 14(6): 1311.
    [10] Xu X H, Zhao H J, Liu, Q L, et al. Mutations of the multi–drug resistance-associated protein ABC transporter gene 5 result in reduction of phytic acid in rice seeds[J]. Theoretical & Applied Genetics, 2009, 119(1): 75–83.
    [11] Laha D, Parvin N, Hofer A, et al. Arabidopsis ITPK1 and ITPK2 have an evolutionarily conserved phytic acid kinase activity[J]. ACS Chemical Biology, 2019, 14(10): 2127–2133.
    [12] Sashidhar N, Harloff H J , Jung C. Identification of phytic acid mutants in oilseed rape (Brassica napus) by large-scale screening of mutant populations through amplicon sequencing[J]. New Phytologist, 2019, 225(5): 2022–2034.
    [13] Sashidhar N, Harloff H J, Potgieter L, et al. Gene editing of three BnITPK genes in tetraploid oilseed rape leads to significant reduction of phytic acid in seeds[J]. Plant Biotechnology Journal, 2020, 18(11): 2241–2250.
    [14] Yamaji N, Takemoto Y, Miyaji T, et al. Reducing phosphorus accumulation in rice grains with an impaired transporter in the node[J]. Nature, 2017, 541(7635): 92.
    [15] Stangoulis J C R, Huynh B L, Welch R M, et al. Quantitative trait loci for phytate in rice grain and their relationship with grain micronutrient content[J]. Euphytica, 2007, 154(3): 289–294.
    [16] Bentsink L, Yuan K, Koornneef M, et al. The genetics of phytate and phosphate accumulation in seeds and leaves of Arabidopsis thaliana, using natural variation[J]. Theoretical & Applied Genetics, 2003, 106(7): 1234–1243.
    [17] Xiao Y, Liu H, Wu L, et al. Genome-wide association studies in maize: Praise and stargaze[J]. Molecular Plant, 2017, 10(3): 359–374.
    [18] Archana, Joshi-Saha, Kandali, et al. Repeat length variation in the 5'UTR of myo-inositol monophosphatase gene is related to phytic acid content and contributes to drought tolerance in chickpea (Cicer arietinum L.)[J]. Journal of Experimental Botany, 2015, 66(19): 5683–5690.
    [19] Perera I, Fukushima A, Akabane T, et al. Expression regulation of myo-inositol 3-phosphate synthase 1 (INO1) in determination of phytic acid accumulation in rice grain[J]. Scientific Reports, 2019, 9(1): 14866.
    [20] Haung W, Lantzsch H J. Sensitive method for the rapid determination of phytate in cereals and cereal products[J]. Journal of the Science of Food & Agriculture, 1983, 34: 1423–1426.
    [21] Tang S, Zhao H, Lu S, et al. Genome- and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus[J]. Molecular Plant, 2021, 14(20): 470–487.
    [22] Alexander D H, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals[J]. Genome Research, 2009, 19(9): 1655–1664.
    [23] Bradbury P J, Zhang Z, Kroon D E, et al. TASSEL: Software for association mapping of complex traits in diverse samples[J]. Bioinformatics, 2007, 23(19): 2633–2635.
    [24] Flint-Garcia, Sherry A, Thornsherry, et al. Structure of linkage disequilibrium in plants[J]. Annual Review of Plant Biology, 2003, 54: 357–374.
    [25] Wang B, Wu Z K, Li Z H, et al. Dissection of the genetic architecture of three seed-quality traits and consequences for breeding in Brassica napus[J]. Plant Biotechnology Journal, 2018, 16(7): 1336–1348.
    [26] Li X K, Guo Z L, Lü Y, et al. Genetic control of the root system in rice under normal and drought stress conditions by genome-wide association study[J]. PLoS Genetics, 2017, 13(7): e1006889.
    [27] Raboy V, Dickinson D B. Effect of phosphorus and zinc nutrition on soybean seed phytic acid and zinc[J]. Plant Physiology, 1984, 75(4): 1094–1098.
    [28] Jackson C A, Windes J M, Obert D, et al. Phosphorus fertility effects on the expression of the low phytic acid barley phenotype[J]. Crop Science, 2009, 49(5): 1800–1806.
    [29] Goldberg S, Sposito G. Achemical model of phosphate adsorption by soils: II. Noncalcareous soils[J]. Soil Science Society of America Journal, 1984, 48(4): 779–783.
    [30] Nagy R, Grob H, Weder B, et al. The Arabidopsis ATP-binding cassette protein AtMRP5/AtABCC5 is a high affinity inositol hexakisphosphate transporter involved in guard cell signaling and phytate storage[J]. Journal of Biological Chemistry, 2009, 284(48): 33614–33622.
    [31] Li W X, Zhao H J, Pang W Q, et al. Seed-specific silencing of OsMRP5 reduces seed phytic acid and weight in rice[J]. Transgenic Research, 2014, 23(4): 585.
    [32] Raboy V. Low phytic acid crops: Observations based on four decades of research[J]. Plants, 2020, 9(2): 140.
    [33] Matthus B, Lsingand R, Fiebig H J, et al. Determination of inositol phosphates IP3-IP6 in rapeseed and rapeseed meal by an HPLC method, Part 2: Investigations of rapeseed and rapeseed meal and comparison with other methods[J]. European Journal of Lipid Science & Technology, 2010, 97(10): 372–374.
  • [1] 张鹏霞周秀文梁雪郭营赵岩李斯深孔凡美 . 小麦产量及氮效率相关性状的全基因组关联分析. 植物营养与肥料学报, 2021, 27(6): 991-1003. doi: 10.11674/zwyf.20549
    [2] 梁桂红华营鹏宋海星张振华 . CACTFTPPCA1 (YACT)、Dof (AAAG)、MYB可能参与甘蓝型油菜对氮胁迫的响应. 植物营养与肥料学报, 2020, 26(2): 338-353. doi: 10.11674/zwyf.19112
    [3] 罗肖艳毕惠惠张旭睿毛培钧王航辉闰永行程西永许海霞 . 小麦耐低钾性状的全基因组关联分析. 植物营养与肥料学报, 2020, 26(3): 401-412. doi: 10.11674/zwyf.19176
    [4] 刘涛蔡秋燕张锡洲李廷轩余海英郭静怡陈光登 . 磷高效型野生大麦根系形态和根系分泌物对低水平植酸态有机磷的响应特征. 植物营养与肥料学报, 2016, 22(6): 1538-1547. doi: 10.11674/zwyf.16006
    [5] 郭丽郭程瑾路文静李小娟肖凯 . 磷转运蛋白基因TaPht1;4 的染色体定位及其与在低磷下与小麦吸磷能力的关系. 植物营养与肥料学报, 2014, 20(4): 877-884. doi: 10.11674/zwyf.2014.0409
    [6] 闫贵欣陈碧云许鲲高桂珍吕培军伍晓明李锋李俊 . 甘蓝型油菜ACCase、 DGAT2和PEPC基因对氮素用量的应答. 植物营养与肥料学报, 2012, 18(6): 1370-1378. doi: 10.11674/zwyf.2012.12094
    [7] . 同磷水平下甘蓝型油菜光合特性的基因型差异研究. 植物营养与肥料学报, 2010, 16(5): 1196-1202. doi: 10.11674/zwyf.2010.0521
    [8] 左青松唐瑶石剑飞杨光惠飞虎冷锁虎* . 甘蓝型油菜不同氮素子粒生产效率品种的氮素分配特性研究. 植物营养与肥料学报, 2009, 15(6): 1395-1400. doi: 10.11674/zwyf.2009.0621
    [9] 石磊梁宏玲徐芳森王运华 . 甘蓝型油菜幼苗体内磷组分差异与磷高效关系的研究 . 植物营养与肥料学报, 2008, 14(2): 351-356. doi: 10.11674/zwyf.2008.0223
    [10] 吴照辉贺立源左雪冬杨建峰门玉英 . 低磷胁迫对不同基因型水稻阶段生物学特征的影响 . 植物营养与肥料学报, 2008, 14(2): 227-234. doi: 10.11674/zwyf.2008.0205
    [11] 年夫照彭慧元徐芳森石磊王运华 . 不同硼效率甘蓝型油菜苗期对硼镁营养的反应. 植物营养与肥料学报, 2004, 10(5): 511-515. doi: 10.11674/zwyf.2004.0512
    [12] 石磊徐芳森刘峰孟金陵赵建伟王运华 . 甘蓝型油菜硼高效基因等位性检测. 植物营养与肥料学报, 2003, 9(2): 229-232. doi: 10.11674/zwyf.2003.0217
    [13] 杜昌文王运华徐芳森王火焰 . 不同硼效率甘蓝型油菜品种中硼的形态及其相互关系. 植物营养与肥料学报, 2002, 8(1): 105-109. doi: 10.11674/zwyf.2002.0120
    [14] 杨玉华王运华杜昌文吴礼树 . 硼对不同硼效率甘蓝型油菜品种细胞壁组成的影响. 植物营养与肥料学报, 2002, 8(3): 340-343. doi: 10.11674/zwyf.2002.0315
    [15] 段海燕徐芳森王运华 . 甘蓝型油菜不同磷效率品种苗期根系生长及磷营养的差异. 植物营养与肥料学报, 2002, 8(1): 65-69. doi: 10.11674/zwyf.2002.0112
    [16] 徐芳森王运华尹文华 . 甘蓝型油菜硼营养高效的遗传分析. 植物营养与肥料学报, 2001, 7(4): 429-434. doi: 10.11674/zwyf.2001.0412
    [17] 杨玉华吴礼树王运华曹享云 . 硼对不同甘蓝型油菜品种细胞壁酶活性的影响. 植物营养与肥料学报, 1999, 5(4): 341-346. doi: 10.11674/zwyf.1999.0408
    [18] 徐芳森王运华李建春 . 甘蓝型油菜硼营养高效在F1代的遗传研究. 植物营养与肥料学报, 1998, 4(3): 305-310. doi: 10.11674/zwyf.1998.0316
    [19] 刘国栋李继云李振声 . 低磷胁迫下小麦根系反应的基因型差异. 植物营养与肥料学报, 1996, 2(3): 212-218. doi: 10.11674/zwyf.1996.0304
    [20] 王震宇张福锁 . 硼对甘蓝型油菜蛋白质的影响. 植物营养与肥料学报, 1995, 1(1): 52-56. doi: 10.11674/zwyf.1995.0108
  • 加载中
图(3)表(4)
计量
  • 文章访问数:  135
  • HTML全文浏览量:  375
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-17
  • 网络出版日期:  2021-11-25
  • 刊出日期:  2021-10-25

正常供磷及低磷胁迫下油菜种子植酸浓度与含量的全基因组关联分析

    作者简介:刘海疆 E-mail: liuhaijiang@webmail.hzau.edu.cn
    通讯作者: 石磊, leish@mail.hzau.edu.cn
  • 1. 华中农业大学作物遗传改良国家重点实验室, 湖北武汉 430070
  • 2. 华中农业大学微量元素研究中心, 湖北武汉430070
  • 3. 武穴市农业农村局, 湖北武穴435400
  • 4. 湖北省农业农村厅油菜办公室, 湖北武汉430070
  • 基金项目: 国家自然科学基金项目(32172662,31972498);中央高校基本科研业务费专项基金资助项目(2662019PY013)。
  • 摘要:   【目的】  鉴定不同磷水平油菜种子植酸浓度 (PA_Conc)和含量 (PA_Cont) 显著关联的 SNP 位点,挖掘候选基因及其优异单倍型,为油菜低植酸品种的培育提供遗传基础和优异基因资源。  【方法】  采用田间试验,测定正常供磷(施磷,P2O5 90 kg/hm2)及低磷(不施磷,P2O5 0 kg/hm2)条件下403个甘蓝型油菜品种种子中的PA_Conc和PA_Cont,利用基于全基因组重测序技术获得的530多万个高质量SNP (single nucleotide polymorphism)标记,用 TASSEL 5.0 软件中的一般线性模型(general linear model,GLM)和混合线性模型(mixed linear model,MLM)进行全基因组关联分析(genome-wide association study, GWAS),鉴定种子植酸浓度和含量显著关联的SNP位点及其候选基因。采用Bcftools提取BnaA07.MRP13基因的SNP,利用Haploview 4.2软件进行单倍型分析。  【结果】  关联群体正常供磷处理PA_Conc变异范围为17.52~81.32 mg/g,PA_Cont变异范围为0.30~0.99 mg/5粒种子;低磷处理PA_Conc变异范围为17.47~49.05 mg/g,PA_Cont变异范围为0.33~0.85 mg/5粒种子。与正常供磷相比,关联群体低磷胁迫下的PA_Conc和PA_Cont均值分别降低13.5%和10.8%。GLM模型检测到56个与油菜种子植酸浓度和含量显著关联的SNP位点,分布在除A08和C03外的17条染色体上,它们对表型变异的解释率为7.33%~14.28%。其中A03染色体上含有最多的SNP (9个)。MLM模型共检测到23个与油菜种子植酸浓度和含量显著关联的SNP位点,分布在A03、A05、A06、A07、A09、C01、C02、C04、C05、C06和C07等11条染色体上,它们对表型变异的解释率为9.6%~14.28%。其中A03染色体上含有8个SNP。对显著SNP位点上下游300 kb范围内的基因进行分析,挖掘出可能与种子中植酸合成相关的14个候选基因,其中包括拟南芥植酸合成通路中的MRPPLCPIP5KSULTR4等同源基因。对BnaA07.MRP13进行候选基因关联及单倍型分析,鉴定出高植酸单倍型“BnaA07.MRP13ContHap1”(TTTA)和低植酸单倍型“BnaA07.MRP13ContHap2”(CCCT)。自然群体中18%的品种(72/403)为低植酸单倍型,种子平均植酸含量为 0.56 mg/5粒种子,极显著低于高植酸单倍型 (平均植酸含量为0.62 mg/5粒种子)。  【结论】  正常供磷和低磷胁迫下,油菜种子植酸浓度和含量均具有广泛的遗传变异。GLM模型和MLM模型分别检测到56和23个与油菜种子植酸浓度和含量显著关联的SNP位点,预测的与油菜种子植酸含量相关的候选基因有14个,发现的BnaA07.MRP13的高和低植酸单倍型为后续油菜种子植酸含量的遗传改良奠定了基础。

    English Abstract

    • 油菜是我国也是世界上重要的油料作物,每年种植面积约533万hm2,总产1348多万t[1-2]。菜籽油是国产食用油的第一大来源,在我国食用油市场中具有举足轻重的地位,因此,推动油菜产业发展是我国植物油供给安全的重要保障[2]。磷是植物生长发育过程中所必需的大量营养元素之一[3]。植酸是种子中磷的主要贮存形式。植酸通常以植酸钾、植酸钙和植酸镁等植酸盐的形式积累在种子中。在种子萌发时,植酸水解成磷酸盐和肌醇,并释放出金属离子,为种子萌发提供养分。油菜种子植酸磷占总磷80%左右,缺磷时油菜种子植酸磷和总磷均显著降低[4-6]。菜籽饼是油菜榨油后的副产品,磷含量较高,其中65%是植酸磷[4]。但是,种子和饼粕中植酸并不能被人和其他非反刍类动物消化吸收,造成了磷和其他矿质营养元素的浪费[7-8]。此外,植酸进入水体还可能导致水体富营养化。因此,培育低植酸油菜品种成为油菜育种的一个重要目标,解析其遗传基础对培育低植酸油菜品种具有重要意义。

      植酸的合成是由葡萄糖-6-磷酸通过肌醇-3-磷酸合成酶(MIPS)氧化环化生成肌醇-3-磷酸开始的。随后,肌醇单磷酸酶(IMP)通过去磷酸化作用将肌醇-3-磷酸转化为游离肌醇。同时,在肌醇激酶(MIK)的作用下,游离肌醇会再次被磷酸化生成肌醇-3-磷酸。之后,肌醇-3-磷酸通过不依赖脂类途径合成肌醇-1,3,4,5,6-五磷酸,而肌醇通过依赖脂类途径也形成肌醇-1,3,4,5,6-五磷酸。在肌醇多磷酸-2-激酶(IPK1)的催化作用下,肌醇-1,3,4,5,6-五磷酸生成植酸。最终,植酸通过多药耐药相关蛋白(MRP5)运输到蛋白贮存液泡(PSV)中,与金属阳离子螯合,以植酸钙、植酸镁和植酸钾等植酸盐形式储存起来[9]。目前,在拟南芥、水稻、玉米以及油菜中已鉴定出几种参与植酸生物合成的酶,如AtITPK1AtITPK2OsMRP5ZmMRP4BnPGK2BnITPK1BnITPK2[10-13]。除植酸合成和运输相关基因外,一些影响种子植酸浓度的磷转运蛋白基因突变体也被鉴定出来。有研究表明,水稻OsSPDT基因影响磷向籽粒的运输,从而影响种子中植酸的积累[14]OsSPDT主要在茎节中表达,负责控制磷向籽粒的分配。与野生型相比,OsSPDT突变体种子中植酸磷降低25%~32%,但种子萌发与成熟期产量并未受到影响[14]

      前人利用遗传群体,分别在水稻、拟南芥、绿豆等植物中鉴定到多个与植酸相关的遗传位点(QTL)[15-16]。在水稻中,Stangoulis等[15]利用粳稻Azucena和籼稻 IR64杂交培育的 DH遗传群体鉴定到2个控制种子植酸含量的QTL,分别位于第5和第12条染色体上,表型解释率分别为24%和15%,其中在第5条染色体上的QTL除与种子植酸含量相关外还和种子中总磷含量密切相关,对总磷含量的表型解释率为20%。在拟南芥中,Bentsink等[16]利用拟南芥Cvi/Ler RIL 群体检测到一个位于第3条染色体上的调控种子植酸含量和叶片植酸含量的多效QTL,该QTL可以分别解释拟南芥种子和叶片植酸含量的61.8%和22.5%的表型变异。相比于连锁分析,关联分析具有定位精度高,无需构建专门的遗传群体,节省时间等优点[17]。Archana等[18]以鹰嘴豆自然群体为材料,通过全基因组关联分析的方法,鉴定到1个与植酸含量相关的位于肌醇单磷酸酶基因的(CaIMP) 5'UTR 区域的简单重复序列标记NCPGR90,该标记简单重复序列长度的变化与鹰嘴豆种子植酸含量的高低相关。Perera等[19]通过GWAS分析鉴定到候选基因肌醇-3-磷酸合酶1 (OsINO1),在植酸积累过程中,与低植酸水稻品种WRC5相比,高植酸水稻品种WRC6的 OsINO1的表达量上调了1.8倍,表明INO1在水稻种子植酸的积累过程中起着重要作用。

      尽管植酸相关QTL在上述作物中已有报道,但在油菜中鲜有研究。本研究以来源于世界各地的403份甘蓝型油菜为关联分析群体,结合全基因组重测序获得的530余万高质量SNP (single nucleotide polymorphism)标记,对田间试验正常磷和低磷处理油菜种子中的PA_Conc和PA_Cont进行全基因组关联分析,鉴定到与植酸合成相关的SNP位点及其候选基因,开展候选基因的关联分析和单倍型分析,鉴定优异单倍型。研究结果将有助于解析甘蓝型油菜种子植酸含量的遗传机制,为低植酸油菜育种提供理论依据。

      • 关联分析群体包括不同来源的油菜品种403份,其中350个半冬性、44个春性、8个冬性和1个未知种性品种。这些品种中,361个来源于中国,21个来自欧洲,8个来自日本,5个来自加拿大,4个来自澳大利亚,3个来自韩国,1个未知。

        供试土壤为水稻土,其基本理化性质:pH为5.61、有机质为7.37 g/kg、全氮为0.91 g/kg、速效磷为11.56 mg/kg、土壤全磷为0.57 g/kg、土壤全钾为13.2 g/kg。

      • 上述油菜自然群体于2018年10月至2019年5月在湖北省武穴市梅川镇(29.85°N,115.55°E)种植,试验设置正常磷(P2O5 90 kg/hm2)和低磷(不施磷肥,P2O5 0 kg/hm2)两个处理,每个处理设置3次重复,采用随机区组排列。每个重复每个品种种植4行,每行8个单株,共32株,行距20 cm,株距20 cm,每17个品种为一个小区,共144个小区。播种前,正常磷处理田块的氮、磷、钾和硼肥按108、90、120 和 1.6 kg/hm2作为基肥施入。施用肥料品种分别为尿素(N 46.6%)、过磷酸钙(P2O5 12%)、氯化钾(KCl 60%)和硼砂(B 15%)。在播种60天后,按72 kg/hm2追施氮素。不施磷处理除不施磷肥外,其余肥料施用与正常磷处理相同,试验田的田间管理措施与当地油菜田常规管理措施一致。待成熟时,选取3株长势一致的植株进行取样,测定种子中植酸浓度并计算植酸含量。

        本试验采用铁沉淀法测定油菜种子中植酸浓度[20],具体步骤如下:用万分之一天平称量5粒种子,记录种子重量,放在1.5 mL离心管中,添加1 mL 0.4 mol/L HCl和1粒钨球,利用TissueLyser ball mill 打样机(Qiagen, Germany)碾磨3 min,以15000 r/min的转速离心5 min。离心后移取500 µL的上清液到1.5 mL离心管。再以15000 r/min的转速离心5 min。最后吸取200 µL的离心后上清液于1.5 mL离心管备用,称为原始样品提取液。同时制备植酸标准溶液,将50% (w/w)植酸溶液稀释得到植酸溶液A (浓度2.096 mg/mL),将植酸溶液A、0.4 mol/L HCl和0.2 mol/L HCl按照1∶1∶14的比例配制为0.131 mg/mL植酸标准溶液,利用0.2 mol/L盐酸按照梯度稀释法,稀释得到4个浓度水平的植酸标准溶液:0.131、0.0655、0.03275、0.016375 mg/mL。向96孔PCR板中,每个孔中添加15 µL水和30 µL 0.2 mol/L HCl,同时添加15 µL的样品提取液。在不同的孔中,添加60 µL植酸标准溶液,每式3份。将96孔PCR板以300 r/min的转速离心20 s,向每个孔中添加120 µL 0.02% (w/v)硫酸铁铵-0.2 mol/L HCl。96孔PCR板盖帽后,在99℃加热30 min,加热后冰浴冷却10 min。随后样品以3000 g相对离心力离心30 min。离心后将96孔PCR板中各样品和植酸标准溶液转移80 µL到96孔酶标板(平底)中。每个孔中添加120 µL显色剂1% (w/v) 2,2'-联吡啶-1% (v/v)巯基乙酸,置于摇床混合10 min。利用KC4酶标仪(Bio-Tek Instruments, USA)读取在519 nm的吸光度,每个孔读取3次数值,以3次重复的平均值作为吸光度数值。计算后得到种子植酸浓度,以mg/g表示,含量以mg/5粒种子为单位。利用 Microsoft Excel 2007软件初步整理表型数据,计算每个表型性状的平均值、标准差、变异系数,利用R中的psych软件包(https://cran.r-project.org/web/packages/psych/)计算每个性状的峰度与偏度,进行正态分布检验。

      • 甘蓝型油菜关联分析群体全基因组重测序共获得一千余万的高质量SNP标记[21]。本研究采用MAF > 5%和缺失率<20%为标准对该群体的所有SNP进行过滤,最终得到530万SNP位点用于后续关联分析。对SNP的命名采用以下方法:chr (染色体chromosome的缩写)开头,加上该SNP所在的染色体,再接下划线“__”,再加上该SNP在染色体对应的物理距离。如,“chrA07__7809297”表示位于A07染色体上7809297bp处的一个SNP位点。为了降低群体结构和亲缘关系对关联分析结果的影响,使用Admixture软件(Alexander et al 2009)进行群体结构分析[22]。使用Tassel 5.0软件对该群体进行亲缘关系分析,计算用于矫正关联分析的K矩阵[23]

      • 在关联分析中,群体结构(Q)及品种间的亲缘关系(K)会导致关联分析结果存在一定的假阳性[24]。本研究采用TASSEL 5.0软件中的只矫正群体结构(Q)的一般线性模型和同时矫正群体结构(Q)和亲缘关系(K)的混合线性模型来进行关联分析。用于筛选显著关联SNP位点的阈值设为1/nn为SNP标记数。根据关联分析所计算的所有–lg(P)观察值和期望值,利用R包GGplot2 (https://cran.r-project.org/web/packages/ggplot2/)和CMplot (https://cran.r-project.org/web/packages/CMplot/)分别绘制Manhattan图和Quantile-Quantile散点图(QQ plot)。

      • 根据 Wang等[25]的研究方法,将目标SNP上下游300 kb距离范围内的基因作为初步候选基因,结合拟南芥同源基因的功能注释,筛选确定候选基因。

      • 根据 Li等[26]的研究方法,利用软件bcftools对目标候选基因BnaA07.MRP13的CDS区域及其前面2 kb范围的SNP进行提取(http://samtools.github.io/bcftools/bcftools.html),将所获得的SNP按照最小基因型频率(minor allele frequency,MAF) > 5%和缺失率(missing rate) < 20%为标准进行过滤,获得高质量SNP标记用于后续的关联分析。使用GLM模型对候选基因内SNP和低磷胁迫下的PA_Cont进行候选基因关联分析,挖掘BnaA07.MRP13内与目标性状相关联的显著SNP,随后利用软件Haploview.4.2对候选基因内的显著SNP进行单倍型分析。

      • 低磷和正常磷处理甘蓝型油菜403个品种种子中PA_Conc和PA_Cont均呈连续性的正态分布,表现出广泛的遗传变异(表1图1)。其中,正常磷处理PA_Conc变异范围为17.52~81.32 mg/g,均值为36.87 mg/g,变异系数为21.37%;低磷处理PA_Conc变异范围为17.47 ~49.05 mg/g,均值为31.91 mg/g,变异系数为24.57%,与正常供磷相比,PA_Conc的均值下降13.5%。正常磷处理每5粒种子PA_Cont变异范围为0.30~0.99 mg,均值为0.65 mg,变异系数为18.85% (表1)。低磷处理每5粒种子PA_Cont的变异范围为0.33~0.85 mg,均值为0.58 mg,变异系数为18.91% (表1)与正常供磷相比,下降10.8%,PA_Cont的均值下降10.8%。与正常磷处理比较,低磷胁迫PA_Conc和PA_Cont分别下降4.96 mg/g和0.07 mg/5粒种子,但PA_Conc变异系数明显增加(表1图1)。

        图  1  正常磷处理和低磷处理甘蓝型油菜关联分析群体种子植酸浓度和含量的频率分布直方图

        Figure 1.  Frequency distribution of phytate concentration and content of an association panel of Brassica napus under sufficient and deficient phosphorus supplies

        表 1  正常磷处理及低磷处理甘蓝型油菜关联分析群体种子植酸浓度和含量

        Table 1.  Seed phytate concentration and content in an association panel of Brassica napus under sufficient and deficient phosphorus supplies

        性状 Trait最大值 Max最小值 Min均值 Mean峰度 Kurtosis偏度 Skewness变异系数 CV (%)
        正常磷种子植酸浓度 PA_Conc_SP (mg/g)81.3217.5236.871.230.6721.37
        低磷种子植酸浓度 PA_Conc_DP (mg/g)49.0517.4731.91–0.440.3824.57
        正常磷种子植酸含量 PA_Cont_SP (mg/5 seeds)0.990.300.65–0.290.1818.85
        低磷种子植酸含量 PA_Cont_DP (mg/5 seeds)0.850.330.58–0.03–0.4018.91
        注(Note):PA—植酸 Phytate; Conc—浓度 Concentration; Cont—含量 Content; SP—充足磷 Sufficient P supply; DP—低磷 Deficient P supply.
      • 采用GLM和MLM模型对正常磷和低磷处理油菜种子的PA_Conc和PA_Cont进行全基因组关联分析(图2)。其中,利用GLM模型鉴定到正常磷处理与PA_Conc显著关联的SNP位点42个,分别分布在A01、A02、A03、A04、A05、A06、A09、A10、C01、C02、C04、C05、C06、C07、C08、C09等16条染色体上,解释表型变异范围为8.90%~13.89% (表2);与PA_Cont显著关联的SNP位点1个,解释表型变异为7.33% (表2)。低磷处理鉴定到与PA_Conc显著关联的SNP位点6个,其中3个SNP位点分别分布在A03和A10两条染色体上,解释表型变异范围为9.36%~14.28% (表2);与PA_Cont显著关联SNP位点7个,其中6个SNP 位点分别分布在A02和A07两条染色体上,解释表型变异范围为8.37%~12.10% (表2)。

        图  2  正常磷和低磷处理油菜种子植酸浓度和含量全基因组关联分析的曼哈顿图和QQ图

        Figure 2.  Manhattan and QQ plots of GWAS on seed phytate concentration and content of Brassica napus under sufficient and deficient P supply

        表 2  正常磷和低磷处理甘蓝型油菜种子植酸浓度和植酸含量显著关联的SNP位点

        Table 2.  Significant SNPs associated with phytate concentration and content in Brassica napus seed under sufficient and deficient phosphorus supplies

        性状 Trait显著关联 SNP Significant SNPPP-value表型解释方差 PVE (%)模型 Model
        低磷植酸浓度
        PA_Conc_DP
        chrA03__273392794.34×10–710.72GLM
        chrA03__273394015.56×10–914.28GLM、MLM
        chrA10__127752942.38×10–712.36GLM
        chrAnn_random__242595924.24×10–711.60GLM
        chrAnn_random__242596053.63×10–711.52GLM
        chrCnn_random__780804851.19×10–79.36GLM
        低磷植酸含量
        PA_Cont_DP
        chrA02__214547667.09×10–79.06GLM
        chrA02__225013374.46×10–711.72GLM
        chrA02__242071274.62×10–712.10GLM
        chrA07__78092971.20×10–711.22GLM、MLM
        chrA07__78093001.20×10–711.18GLM、MLM
        chrA07__78093162.65×10–710.88GLM
        chrUnn_random__62233776.45×10–78.37GLM
        正常磷植酸浓度
        PA_Conc_SP
        chrA01__118142605.36×10–79.50GLM
        chrA02__241708813.53×10–79.60GLM
        chrA03__170574452.73×10–912.70GLM、MLM
        chrA03__191645673.12×10–913.05GLM、MLM
        chrA03__191646571.14×10–913.86GLM、MLM
        chrA03__255851664.12×10–710.07GLM、MLM
        chrA03_random__12899641.29×10–913.89GLM、MLM
        chrA03_random__12911852.46×10–811.56GLM、MLM
        chrA03_random__12912071.71×10–812.04GLM、MLM
        chrA04__44209445.24×10–710.01GLM
        chrA05__48486977.38×10–811.02GLM、MLM
        chrA06__163346619.56×10–912.35GLM、MLM
        chrA06__163346677.38×10–912.56GLM、MLM
        chrA09__62094155.27×10–810.01GLM、MLM
        chrA10__40489105.80×10–79.43GLM
        chrA10__40489196.01×10–79.40GLM
        chrA10__40489304.83×10–79.54GLM
        chrC01__74735651.35×10–811.66GLM、MLM
        chrC02__307762427.91×10–811.06GLM、MLM
        chrC02__316941602.42×10–79.60GLM、MLM
        chrC04__371174574.70×10–78.92GLM
        chrC04__434268206.38×10–710.23GLM
        chrC04__442705494.62×10–810.82GLM、MLM
        chrC05__4340777.13×10–811.07GLM、MLM
        chrC06__290649443.25×10–810.67GLM、MLM
        chrC06__306154366.57×10–79.03GLM
        chrC06__322088457.19×10–79.30GLM
        续表 2 Table 2 continued
        性状 Trait显著关联 SNP Significant SNPPP-value表型解释方差 PVE (%)模型 Model
        chrC06__354253063.63×10–811.00GLM、MLM
        chrC06__354253126.56×10–810.71GLM、MLM
        chrC06_random__30871893.65×10–78.97GLM
        chrC07__348576545.82×10–79.73GLM
        chrC07__349563611.74×10–710.32GLM、MLM
        chrC07__371987295.14×10–79.14GLM
        chrC08__314146545.27×10–79.53GLM
        chrC08_random__44363911.67×10–812.11GLM
        chrC09__30881134.45×10–710.11GLM
        chrCnn_random__105359003.69×10–78.90GLM
        chrCnn_random__167746242.10×10–811.47GLM
        chrCnn_random__297580175.64×10–78.92GLM
        chrCnn_random__322494451.58×10–79.16GLM
        chrCnn_random__331278745.16×10–79.94GLM
        chrCnn_random__731820295.79×10–79.12GLM
        正常磷植酸含量
        PA_Cont_SP
        chrC09__373436224.71×10–77.33GLM
        注(Note):PA—植酸 Phytate; Conc—浓度 Concentration; Cont—含量Content; SP—充足磷 Sufficient P supply; DP—低磷 Deficient P supply; PVE—表型解释方差 Phenotypic variance explained.

        相比于GLM模型,MLM模型对群体结构和亲缘关系均进行了矫正,鉴定到的显著关联的位点较少,但MLM模型鉴定到的23个显著SNP位点均与GLM共定位(表2)。其中,正常磷处理与PA_Conc显著关联的SNP位点20个,分别分布在A03、A05、A06、A09、C01、C02、C04、C05、C06、C07等10条染色体上,解释表型变异范围为9.60%~13.89% (表2)。MLM模型没有检测到与正常磷PA_Cont显著关联的SNP位点,低磷处理鉴定到1个位于A03染色体上与PA_Conc显著关联的SNP位点,该位点可解释的表型变异为14.28% (表2);2个与PA_Cont显著关联的SNP位点,解释的表型变异分别为11.18%和11.22% (表2)。

      • 基于甘蓝型油菜Darmor-bzh参考基因组序列,对显著关联SNP位点前后300 kb范围内的基因进行分析,同时结合拟南芥同源基因的基因功能注释,最终确定14个参与植酸合成和转运途径的基因作为候选基因(表3)。由GLM和MLM模型共同检测到的植酸浓度和含量显著关联的23个SNP位点的候选基因有结合盒转运蛋白MRP13 (BnaA07g07310DBnaA07g07320D)、MRP12 (BnaA07g07330D),肌醇1,3,4磷酸-5/6-激酶蛋白(BnaA10g17710D),磷脂酰肌醇磷酸酶、磷脂酶C (BnaC09g33980DBnaC09g34000DBnaC09g34010D),磷酸甘油酸变位酶家族蛋白(BnaUnng04100D)以及磷转运相关的蛋白(BnaC06g37190D)等的编码基因。此外,由GLM模型单独鉴定到的油菜种子植酸浓度和含量显著关联SNP位点的候选基因有液泡转运蛋白(BnaC07g30170D)、磷脂酰肌醇-3,4-激酶(BnaC07g30200D)、磷脂酰肌醇单磷酸(BnaC04g24700D)以及类硫酸盐转运蛋白(BnaA03g04410D)的编码基因(表3)。

        表 3  正常磷和低磷处理甘蓝型油菜种子植酸浓度和含量显著关联 SNP 的候选基因

        Table 3.  Candidate genes underlying the significant SNPs for seed phytate concentration and content of Brassica napus under sufficient and deficient phosphorus supply conditions

        性状Trait显著 SNP
        Significant SNP
        油菜基因
        Genes in Brassica napus
        模型
        Model
        基因功能注释
        Gene function annotation
        正常磷植酸含量
        PA_Cont_SP
        chrC06__35245326BnaC06g37190DGLM、MLM磷酸盐转运蛋白1;9
        Phosphate transporter 1;9 (PHT1;9)
        chrC09__37343622BnaC09g33980DGLM、MLM磷脂酰肌醇特异性磷脂酶C1
        Phospholipase C1 (PLC1)
        chrC09__37343622BnaC09g34000DGLM、MLM磷脂酰肌醇特异性磷脂酶C5
        Phospholipase C5 (PLC5)
        chrC09__37343622BnaC09g34010DGLM、MLM磷脂酰肌醇特异性磷脂酶C4
        Phospholipase C4 (PLC4)
        低磷植酸含量
        PA_Cont_DP
        chrC04__9656140BnaC04g24700DGLM磷脂酰肌醇单磷酸5-激酶-4
        Phosphatidylinositol mono phosphate5-kinase-4 (PIP5K4)
        chrC07__34857654BnaC07g30200DGLM磷脂酰肌醇-3, 4-激酶
        Phosphatidylinositol -3- and 4-kinase
        chrC07__34956361BnaC07g30170DGLM液泡膜铁转运家族蛋白
        Vacuolar iron transporter family protein (VIT)
        chrA03__19164557BnaA03g04410DMLM硫酸盐转运蛋白4;1
        Sulfate transporter 4;1 (SULTR4;1)
        chrA07__7809297BnaA07g07310DGLM、MLM多药耐药相关蛋白13
        Multidrug resistance-associated protein 13 (MRP13)
        chrA07__7809297BnaA07g07320DGLM、MLM多药耐药相关蛋白13
        Multidrug resistance-associated protein 13 (MRP13)
        chrA07__7809297BnaA07g07330DGLM、MLM多药耐药相关蛋白12
        Multidrug resistance-associated protein 12 (MRP12)
        chrA02__24207127BnaA02g34000DGLM、MLM3'(2'),5'-二磷酸核苷酸酶活性
        3'(2'),5'-bisphosphate nucleotidase activity
        chrA02__24207127BnaUnng04100DGLM、MLM磷酸甘油酸变位酶家族蛋白,参与InsP6到InsP7的转换和InsP7到InsP8的转换
        Phosphoglycerate mutase-like family protein
        低磷植酸浓度
        PA_Conc_DP
        chrA10__12775294BnaA10g17710DGLM、MLM肌醇1,3,4磷酸-5/6-激酶蛋白
        Inositol 1,3,4-trisphosphate -5/6-kinase family protein
        注(Note):PA—植酸 Phytate; Conc—浓度 Concentration; Cont—含量 Content; SP—充足磷 Sufficient P supply; DP—低磷 Deficient P supply.
      • 在A07染色体上由GLM和MLM模型同时鉴定到低磷处理植酸含量显著关联SNP位点chrA07__7809297,确定多药耐药相关蛋白的基因MRP13(BnaA07.MRP13)为该位点的候选基因(表3)。BnaA07.MRP13的CDS区域及其前面2 kb范围内共有19个SNP位点,将这19个SNP位点与低磷处理植酸含量进行候选基因关联分析。结果表明,共有4个SNP位点与低磷胁迫下植酸含量显著相关,且这4个SNP表现出很强的连锁不平衡(图3a)。对这4个SNP位点进行单倍型分析,共分成2种单倍型:高植酸单倍型“BnaA07.MRP13ContHap1”(TTTA)和低植酸单倍型“BnaA07.MRP13ContHap2”(CCCT) (图3b)。高植酸单倍型“BnaA07.MRP13ContHap1”具有218个油菜品种,其5粒种子的平均植酸含量为0.62 mg,低植酸单倍型“BnaA07.MRP13ContHap2”具有72个油菜品种,其5粒种子的平均植酸含量为0.56 mg (其他油菜品种SNP缺失,不计算在内)。含有“BnaA07.MRP13ContHap1”的品种在低磷胁迫下的植酸含量要显著低于含有“BnaA07.MRP13ContHap2”的品种(P = 0.0064) (图3b)。

        图  3  BnaA07.MRP13的候选基因关联分析(a)和单倍型分析(b)

        Figure 3.  Candidate genes association analysis (a) and haplotype analysis (b) of BnaA07.MRP13

      • 磷肥的施用量显著影响作物种子植酸浓度和含量[27-28]。研究表明,种子中植酸浓度与磷肥的施用量、磷肥在土壤中的有效性及溶解度显著相关[29]。本研究正常磷与低磷处理分别检测到42和6个与种子植酸浓度显著关联的SNP位点,分别检测到1和7个与种子植酸含量显著关联的SNP位点(表2)。但是,在正常磷和低磷胁迫下没有鉴定到相同的SNP位点。在 A07染色体上鉴定到1个低磷处理PA_Cont显著关联的SNP位点和由 3个邻近的 SNP 位点组成的显著 SNP 簇(表2)。这3个SNP位点中,chrA07__7809297和chrA07__7809300同时被GLM和MLM模型检测到,预测BnaA07.MRP12 (BnaA07g07330D) 和BnaA07.MRP13 (BnaA07g07310DBnaA07g0732D)为该位点的候选基因(表3)。拟南芥AtMRP5编码ABC 转运子家族的多药耐药相关蛋白,负责将合成的植酸运输至蛋白贮存液泡[30]。与野生型比较,水稻OsMRP5突变体植株籽粒植酸浓度降低35.8%~71.9%,无机磷浓度显著提高,且无机磷的增加量要高于植酸浓度下降量[31]BnaA07g07330DBnaA07g07310DBnaA07g0732D均属于MRP家族基因,很有可能参与油菜种子中植酸的运输(表3)。低磷处理PA_Cont显著关联SNP鉴定到硫酸盐转运蛋白家族SULTR4;1(BnaA03g04410D) (表3)。水稻OsSPDT基因编码硫酸盐转运蛋白,与野生型比较其突变体种子植酸浓度降低20%,而产量无显著降低,这有利于减少水稻种子中植酸磷通过人和动物流向环境,减轻环境污染[14]。此外,低磷处理PA_Cont显著关联SNP位点还鉴定到液泡膜铁转运家族蛋白(VIT) (BnaC07g30170D)、肌醇或磷脂酰肌醇磷酸酶(BnaA02g34000D)以及参与将肌醇-6-磷酸(InsP6)合成肌醇-7-磷酸(InsP7)的激酶(BnaUnng04100D) (表3)。这些候选基因及其功能需要通过CRISPR-Cas9等分子生物学技术进一步验证。

        植酸含量是种子植酸浓度和种子质量的乘积,因此,种子植酸含量与植酸浓度和种子质量均相关。如果种子植酸含量的差异主要由植酸浓度决定,两者的显著SNP位点可能一样;如果与两者均有关,或者主要由种子质量决定,两者的显著相关的SNP位点可能就不一样(表2)。

      • 作物种子中的磷是全球磷循环的主要驱动力,通过降低种子中植酸的浓度可以减少有限磷资源的浪费[32]。Matthus等[33]利用高效液相色谱法研究油菜种子的磷组成,发现油菜种子中的植酸浓度要显著高于其他作物。研究人员利用CRISPR-Cas9技术构建了甘蓝型油菜BnITPK1BnITPK4基因的双突变体,结果表明,与野生型比较,双突变体种子植酸浓度显著降低27.2%~35.3%;此外,双突变体种子萌发、千粒重、株高、含油量以及产量与野生型无显著差异[12-13]。这表明,通过编辑油菜中BnITPK1BnITPK4可以实现低植酸育种。在本研究中,通过低磷胁迫油菜种子植酸含量鉴定到植酸合成最后一步的MRP家族基因BnaA07.MRP13BnaA07.MRP13基因的候选基因关联分析及单倍型分析结果表明,与低磷胁迫下植酸含量显著相关的4个SNP构成2种单倍型:“BnaA07.MRP13ContHap1”(TTTA)和“BnaA07.MRP13ContHap 2”(CCCT)。其中,“BnaA07.MRP13ContHap2”单倍型属于低植酸单倍型,对于后续油菜低植酸育种研究具有重要意义。本研究油菜自然群体正常磷处理PA_Conc的变异幅度为17.52~81.32 mg/g,低磷处理PA_Conc变异范围为17.47~49.05 mg/g。这些低植酸品种也可以作为低植酸育种的种质资源。

      • 甘蓝型油菜关联分析群体正常磷和低磷处理种子植酸浓度和植酸含量均具有广泛的遗传变异。全基因组关联分析共检测到56个与种子植酸浓度和植酸含量显著关联的 SNP 位点,筛选和鉴定了14个与油菜种子植酸合成相关的候选基因,通过候选基因关联分析以及单倍型分析鉴定了油菜种子高和低植酸单倍型“BnaA07.MRP13ContHap1”(TTTA)和“BnaA07.MRP13ContHap2”(CCCT),这为解析油菜种子植酸浓度和含量的遗传机制奠定了基础,为低植酸品种的培育提供了优异种质和基因资源。但是,这些候选基因的功能及其决定高和低植酸单倍型的分子机制还不清楚,在今后的研究中,可利用 CRISPR-Cas9 等技术创制候选基因的突变体,深入开展基因功能研究,明确这些基因影响种子植酸含量的生物学机制。

    参考文献 (33)

    目录

      /

      返回文章
      返回