欢迎访问作物学报，今天是
作物学报

作物学报 ›› 2019, Vol. 45 ›› Issue (9): 1303-1310.doi: 10.3724/SP.J.1006.2019.94021

• 作物遗传育种·种质资源·分子遗传学 • 上一篇    下一篇

甘蓝型油菜角果长度性状的全基因组关联分析

孙程明1,2,陈松1,彭琦1,张维1,易斌2,*(),张洁夫1,*(),傅廷栋2   

  1. 1 江苏省农业科学院经济作物研究所/农业部长江下游棉花与油菜重点实验室/江苏省现代作物生产协同创新中心, 江苏南京 210014
    2 华中农业大学植物科学技术学院/作物遗传改良国家重点实验室, 湖北武汉 430070
  • 收稿日期:2019-02-12 接受日期:2019-05-12 出版日期:2019-09-12 网络出版日期:2019-05-17
  • 通讯作者: 易斌,张洁夫
  • 作者简介:E-mail: suncm8331537@gmail.com
  • 基金资助:
    本研究由国家重点研发计划项目(2018YFD0100601);国家现代农业产业技术体系建设专项资助(CARS-12)

Genome-wide association study of silique length in rapeseed (Brassica napus L.)

SUN Cheng-Ming1,2,CHEN Song1,PENG Qi1,ZHANG Wei1,YI Bin2,*(),ZHANG Jie-Fu1,*(),FU Ting-Dong2   

  1. 1 Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences/Key Laboratory of Cotton and Rapeseed (Nanjing), Ministry of Agriculture/Jiangsu Collaborative Innovation Center for Modern Crop Production, Nanjing 210014, Jiangsu, China
    2 National Key Laboratory of Crop Genetic Improvement/College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, Hubei, China;
  • Received:2019-02-12 Accepted:2019-05-12 Published:2019-09-12 Published online:2019-05-17
  • Contact: Bin YI,Jie-Fu ZHANG
  • Supported by:
    The study was supported by the National Key Research and Development Program of China(2018YFD0100601);the China Agriculture Research System(CARS-12)

RichHTML

31

PDF

936

摘要：

角果长度是油菜重要的农艺性状, 适度增加角果长度有利于扩大角果库容量, 增加光合面积, 提高油菜的籽粒产量。本研究利用Illumina 60K SNP芯片对496份具有代表性的油菜资源进行基因型分析, 考察群体在4个环境中的角果长度表型, 利用MLM和GLM模型进行全基因组关联分析。结果表明, MLM模型检测到7个位点, 联合解释25.01%的表型变异; GLM模型检测到25个位点, 联合解释41.77%的表型变异。合并共同位点后得到27个位点, 其中7个与前人报道的QTL重叠, 其余20个是新鉴定的位点。效应最大的位点Bn-A09-p29991443位于A09染色体, 在MLM和GLM模型中分别解释13.89%和12.86%的表型变异, 携带其优异等位基因的材料平均角果长度增加0.89 cm。同时, 在该位点附近找到已克隆的油菜角果长度基因ARF18和BnaA9.CYP78A9。此外, 在5个位点附近发现拟南芥已知角果长度基因GID1b、FUL、EOD3、DOF4.4和GA20ox1的同源拷贝。本研究结果有助于解析角果长度的遗传基础, 为研究角果长度的调控机理, 指导角果长度的遗传改良打下基础。

关键词: 甘蓝型油菜, 产量, 角果长度, 关联分析, SNP标记

Abstract:

Silique length is a key agronomic trait of rapeseed. Moderately increasing silique length is conducive to high seed yield by enlarging photosynthetic area and seed volume. A collection of 496 representative rapeseed accessions was genotyped by the Illumina 60K SNP array, and phenotyped for silique length in four environments. The genome-wide association study (GWAS) of silique length was performed via the MLM (Mixed linear model) and GLM (General linear model). Seven loci and 25 loci were detected with MLM and GLM, which explained 25.01% and 41.77% of the phenotypic variance, respectively. Combining the common loci between two models, we finally got 27 unique loci, of which seven were overlapped with reported QTLs, and 20 were new one. Bn-A09-p29991443, the most effective locus, was located on chromosome A09, accounting for 13.89% and 12.86% of the phenotypic variance in MLM and GLM, respectively. Silique length of accessions with the favorable allele of Bn-A09-p29991443 was averagely 0.89 cm longer than that with the unfavorable allele. The cloned silique length genes ARF18 and BnaA9.CYP78A9 in rapeseed was found to be colocalized with Bn-A09-p29991443. Besides, five candidates including GID1b, FUL, EOD3, DOF4.4 and GA20ox1, orthologous to documented Arabidopsis silique length genes, were found near our GWAS loci. The results provide an insight into the genetic basis of silique length and lay a foundation for further mechanism exploration and breeding for this trait in B. napus.

Key words: Brassica napus L, yield, silique length, GWAS, SNP

表1

关联群体角果长度性状的统计分析"

环境
Environment		 最小值
Min.		 最大值
Max.		 平均值±标准差
Mean ± SD		 变异系数
CV
2013/2014 Nanjing 3.42 11.37 5.40±1.01 0.19
2013/2014 Nanjing 3.28 15.80 6.39±1.25 0.20
2014/2015 Taizhou 3.31 10.18 5.79±0.88 0.15
2015/2016 Taizhou 4.01 9.94 5.90±0.77 0.13

图1

关联群体在4个环境的角果长度分布 A和B分别指群体在南京和泰州2年的角果长度频数分布图。"

表2

4个环境的角果长度相关系数"

2012/2013 Nanjing 2013/2014
Nanjing		 2014/2015 Taizhou
2013/2014 Nanjing 0.57***
2014/2015 Taizhou 0.64*** 0.70***
2015/2016 Taizhou 0.56*** 0.57*** 0.75***

表3

MLM角果长度显著关联位点(BLUP) "

标记
Marker		 染色体
Chr.		 位置
Position		 -lg (P) 表型变异
R2 (%)		 环境
Environment		 已报道QTL
Reported QTL
Bn-A07-p8572785 A07 10,018,929 5.80 5.10 13NJ/15TZ
Bn-A09-p28925363 A09 26,858,796 5.15 4.92 13NJ/15TZ [8, 22]
Bn-A09-p29991443 A09 27,815,620 13.59 13.89 13NJ/14NJ/15TZ [6, 8, 22, 23]
Bn-scaff_15712_2-p492440 C02 38,697,584 4.67 3.54
Bn-scaff_17869_1-p813291 C04 9,884,407 5.41 4.73 14NJ /15TZ
Bn-scaff_20817_1-p60579 C04 48,635,772 4.56 4.27 15TZ
Bn-scaff_15576_1-p68473 C09 41,121,951 4.29 3.21 13NJ/15TZ

图2

油菜角果长度全基因组关联分析(BLUP) A: 角果长度MLM曼哈顿图; B: 角果长度MLM QQ图; C: 角果长度GLM曼哈顿图; D: 角果长度GLM QQ图。水平线代表Bonferroni阈值。"

表4

GLM角果长度显著关联位点(BLUP) "

标记
Marker		 染色体
Chr.		 位置
Position		 -lg (P) 表型变异
R2 (%)		 环境
Environment		 已报道QTL
Reported QTL
Bn-A02-p11454573 A02 8,166,563 4.61 3.72 13NJ/14NJ/15TZ
Bn-A03-p4408895 A03 3,938,915 4.39 2.91
Bn-scaff_16092_1-p866917 A03 18,418,252 4.74 4.66 13NJ/15TZ
Bn-scaff_27914_1-p34836 A06 15,544,511 6.21 5.30 15TZ/16TZ
Bn-A06-p23042849 A06 22,092,944 4.72 3.85 15TZ/16TZ
Bn-A07-p1228602 A07 860,252 4.51 3.76 16TZ
Bn-A10-p12099059 A07 2,762,935 5.02 3.38 15TZ/16TZ
Bn-A07-p7606228 A07 9,162,030 6.03 4.90 13NJ/15TZ/16TZ [23]
Bn-A07-p8572785 A07 10,018,929 9.38 7.76 13NJ/15TZ/16TZ
Bn-A07-p10340211 A07 11,509,208 4.45 3.70 15TZ/16TZ
Bn-A07-p13957160 A07 15,884,413 5.18 4.20 15TZ/16TZ [23]
Bn-A07-p18319586 A07 20,221,220 5.14 4.33 14NJ/15TZ [6, 23]
Bn-A09-p3051349 A09 2,971,335 4.31 3.32 15TZ
Bn-A09-p28925363 A09 26,858,796 6.08 4.89 13NJ/15TZ/16TZ [8, 22]
Bn-A09-p29991443 A09 27,815,620 16.27 12.86 13NJ/14NJ/15TZ/16TZ [6, 8, 22, 23]
Bn-A10-p3966740 A10 885,133 6.47 4.46 13NJ/14NJ/15TZ/16TZ [6]
Bn-A10-p15793623 A10 15,756,213 4.60 3.65
Bn-scaff_17177_1-p381225 C02 44,843,894 5.12 4.13 15TZ
Bn-A05-p4304172 C04 6,449,454 6.96 5.75 14NJ/15TZ/16TZ
Bn-scaff_27914_1-p9919 C04 28,255,420 7.06 6.09 13NJ/14NJ/15TZ/16TZ
Bn-scaff_20817_1-p60579 C04 48,635,772 7.98 6.89 13NJ/14NJ/15TZ/16TZ
Bn-scaff_16069_1-p1668600 C07 38,086,007 6.38 5.29 15TZ/16TZ
Bn-scaff_16069_1-p3731985 C07 40,142,298 5.92 4.63 13NJ/14NJ/15TZ/16TZ
Bn-scaff_16361_1-p241830 C08 27,753,144 4.34 3.48 13NJ [9]
Bn-scaff_15576_1-p68473 C09 41,121,951 5.47 3.72 14NJ/15TZ/16TZ

表5

角果长度关联位点候选基因信息"

标记
Marker		 油菜基因
Rapeseed gene		 染色体
Chr.		 位置
Position		 拟南芥同源基因
Ar. homolog		 参考基因组
Ref. genome
Bn-A07-p13957160 BnaA07g19530 A07 15,590,253 GID1b Darmor-bzh
Bn-A09-p3051349 BnaA09g05500 A09 2,718,832 FUL Darmor-bzh
Bn-scaff_20817_1-p60579 BnaC04g50960 C04 48,343,005 EOD3 Darmor-bzh
Bn-scaff_16069_1-p1668600 BnaC07g36530 C07 38,575,790 DOF4.4 Darmor-bzh
Bn-scaff_16069_1-p3731985 BnaC07g39650 C07 40,392,394 GA20ox1 Darmor-bzh
Bn-A09-p29991443 BnA09g0377700 A09 36,712,502 ARF18 ZS11
Bn-A09-p29991443 BnA09g0377760 A09 36,740,578 CYP78A9 ZS11
[1] 王汉中 . 我国油菜产业发展的历史回顾与展望. 中国油料作物学报, 2010,32:300-302.
Wang H Z . Review and future development of rapeseed industry in China. Chin J Oil Crop Sci, 2010,32:300-302 (in Chinese with English abstract).
[2] 冷锁虎, 唐瑶, 李秋兰, 左青松, 杨萍 . 油菜的源库关系研究: I. 角果大小对油菜后期源库的调节. 中国油料作物学报, 2005,27(3):37-40.
Leng S H, Tang Y, Li Q L, Zuo Q S, Yang P . Studies on source and sink of rapeseed: I. Regulation of pod size on source and sink in rapeseed after flowering. Chin J Oil Crop Sci, 2005,27(3):37-40 (in Chinese with English abstract).
[3] 王春丽, 海江波, 田建华, 杨建利, 赵晓光 . 油菜终花后角果和叶片光合对籽粒产量和品质的影响. 西北植物学报, 2014,34:1620-1626.
Wang C L, Hai J B, Tian J H, Yang J L, Zhao X G . Influence of silique and leaf photosynthesis on yield and quality of seed of oilseed rape (Brassica napus L.) after flowering. Acta Bot Boreal-Occident Sin, 2014,34:1620-1626 (in Chinese with English abstract).
[4] Udall J A, Quijada P A, Lambert B, Osborn T C . Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 2. Identification of alleles from unadapted germplasm. Theor Appl Genet, 2006,113:597-609.
[5] Yang P, Shu C, Chen L, Xu J, Wu J, Liu K . Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L.). Theor Appl Genet, 2012,125:285-296.
[6] 漆丽萍 . 甘蓝型油菜株型与角果相关性状的QTL分析. 华中农业大学博士学位论文, 湖北武汉, 2014.
Qi L P . QTL Analysis for the Traits Associated with Plant Architecture and Silique in Brassica napus L. PhD Dissertation of Huazhong Agricultural University, Wuhan, Hubei, China, 2014. (in Chinese with English abstract).
[7] Wang X D, Chen L, Wang A N, Wang H, Tian J H, Zhao X P, Chao H B, Zhao Y J, Zhao W G, Xiang J, Gan J P, Li M T . Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus. BMC Plant Biol, 2016,16:71.
[8] Fu Y, Wei D Y, Dong H L, He Y J, Cui Y X, Mei J Q, Wan H F, Li J N, Snowdon R, Friedt W, Li X R, Qian W . Comparative quantitative trait loci for silique length and seed weight in Brassica napus. Sci Rep, 2015,5:14407
[9] Yang Y, Shen Y S, Li S D, Ge X H, Li Z Y . High density linkage map construction and QTL detection for three silique-related traits in Orychophragmus violaceus derived Brassica napus population. Front Plant Sci, 2017,8:1512
[10] Li H, Peng Z Y, Yang X H, Wang W D, Fu J J, Wang J H, Han Y J, Chai Y C, Guo T T, Yang N, Liu J, Warburton M, Cheng Y B, Hao X M, Zhang P, Zhao J Y, Liu Y J, Wang G Y, Li J S, Yan J B . Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat Genet, 2013,45:43-52.
[11] Wen W W, Li D, Li X, Gao Y Q, Li W Q, Li H H, Liu J, Liu H J, Chen W, Luo J, Yan J B . Metabolome-based genome-wide association study of maize kernel leads to novel biochemical insights. Nat Commun, 2014,5:3438
[12] Dong H J, Zhao H, Xie W B, Han Z M, Li G W, Yao W, Bai X F, Hu Y, Guo Z L, Lu K, Yang L, Xing Y Z . A novel tiller angle gene, TAC3, together with TAC1 and D2 largely determine the natural variation of tiller angle in rice cultivars. PLoS Genet, 2016,12:e1006412.
[13] Huang X H, Kurata N, Wei X H, Wang Z X, Wang A H, Zhao Q, Zhao Y, Liu K Y, Lu H Y, Li W J, Guo Y L, Lu Y Q, Zhou C C, Fan D L, Weng Q J, Zhu C R, Huang T, Zhang L, Wang Y C, Feng L, Furuumi H, Kubo To, Miyabayashi T, Yuan X P, Xu Q, Dong G J, Zhan Q L, Li C Y, Fujiyama A, Toyoda A, Lu T T, Feng Q, Qian Q, Li J Y, Han B . A map of rice genome variation reveals the origin of cultivated rice. Nature, 2012,490:497-501.
[14] Sun C M, Wang B Q, Yan L, Hu K N, Liu S, Zhou Y M, Guan C Y, Zhang Z Q, Li J N, Zhang J F, Chen S, Wen J, Ma C Z, Tu J X, Shen J X, Fu T D, Yi B . Genome-wide association study provides insight into the genetic control of plant height in rapeseed (Brassica napus L.). Front Plant Sci, 2016,7:1102.
[15] Sun C M, Wang B Q, Wang X H, Hu K N, Li K D, Li Z Y, Li S, Yan L, Guan C Y, Zhang J F, Zhang Z Q, Chen S, Wen J, Tu J X, Shen J X, Fu T D, Yi B . Genome-wide association study dissecting the genetic architecture underlying the branch angle trait in rapeseed (Brassica napus L.). Sci Rep, 2016,6:33673.
[16] Ihaka R, Gentleman R . R: a language for data analysis and graphics. J Comput Graph Stat, 1996,5:299-314.
[17] Merk H L, Yarnes S C, Van Deynze A, Tong N, Menda N, Mueller L A, Mutschler M A, Loewen S A, Myers J R, Francis D M . Trait diversity and potential for selection indices based on variation among regionally adapted processing tomato germplasm. J Am Soc Hortic Sci, 2012,13:427-437.
[18] Evanno G, Regnaut S, Goudet J . Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol, 2005,14:2611-2620.
[19] Hardy O J, Vekemans X . SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes, 2002,2:618-620.
[20] Bradbury P J, Zhang Z, Kroon D E, Casstevens T M, Ramdoss Y, Buckler E S . TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics, 2007,23:2633-2635.
[21] Turner S D . qqman: an R package for visualizing GWAS results using QQ and manhattan plots. BioRxiv, 2014. doi: 10.1101/005165.
[22] 周庆红, 周灿, 郑伟, 付东辉 . 甘蓝型油菜角果长度全基因组关联分析. 中国农业科学, 2017,50:228-239.
Zhou Q H, Zhou C, Zheng W, Fu D H . Genome wide association analysis of silique length in Brassica napus L. Sci Agric Sin, 2017,50:228-239 (in Chinese with English abstract).
[23] Dong H, Tan C, Li Y, He Y, Wei S, Cui Y, Chen Y, Wei D, Fu Y, He Y, Wan H, Liu Z, Xiong Q, Lu K, Li J, Qian W . Genome-wide association study reveals both overlapping and independent genetic loci to control seed weight and silique length in Brassica napus. Front Plant Sci, 2018,9:921.
[24] Liu J, Hua W, Hu Z, Yang H, Zhang L, Li R, Deng L, Sun X, Wang X, Wang H . Natural variation in ARF18 gene simultaneously affects seed weight and silique length in polyploid rapeseed. Proc Natl Acad Sci USA, 2015,112:5123-5132.
[25] Shi L, Song J, Guo C, Wang B, Guan Z, Yang P, Chen X, Zhang Q, King G J, Wang J, Liu K . A CACTA-like transposable element in the upstream region of BnaA9. CYP 78A9 acts as an enhancer to increase silique length and seed weight in rapeseed. Plant J, 2019. doi: 10.1111/tpj.14236.
[26] Gu Q, Ferrándiz C, Yanofsky M F, Martienssen R . The FRUITFULL MADS-box gene mediates cell differentiation during Arabidopsis fruit development. Development, 1998,125:1509-1517.
[27] Fang W, Wang Z, Cui R, Li J, Li Y . Maternal control of seed size by EOD3/CYP78A6 inArabidopsis thaliana. Plant J, 2012,70:929-939.
[28] Zou H F, Zhang Y Q, Wei W, Chen H W, Song Q X, Liu Y F, Zhao M Y, Wang F, Zhang B C, Lin Q . The transcription factor AtDOF4.2 regulates shoot branching and seed coat formation in Arabidopsis. Biochem J, 2013,449:373-388.
[29] Griffiths J, Murase K, Rieu I, Zentella R, Zhang Z L, Powers S J, Gong F, Phillips A L, Hedden P, Sun T P . Genetic characterization and functional analysis of the GID1 gibberellin receptors in Arabidopsis. Plant Cell, 2006,18:3399-3414.
[30] Rieu I, Ruiz-Rivero O, Fernandez-Garcia N, Griffiths J, Powers S J, Gong F, Linhartova T, Eriksson S, Nilsson O, Thomas S G . The gibberellin biosynthetic genes AtGA20ox1 and AtGA20ox2 act, partially redundantly, to promote growth and development throughout the Arabidopsis life cycle. Plant J, 2008,53:488-504.
[1] 习千辉, 徐梓瑗, 刘梦梦, 王宏艺, 郎凯琳, 井震海, 陈锋, 赵磊. 小麦籽粒铜含量的全基因组关联分析及候选基因预测[J]. 作物学报, 2026, 52(6): 1604-1617.
[2] 毛嘉琦, 黄朋雨, 赵佳佳, 郑兴卫, 武棒棒, 郝宇琼, 屈非, 刘成, 马朋涛, 郑军. 山西小麦品种白粉病抗性评价及抗病基因分子检测[J]. 作物学报, 2026, 52(6): 1669-1681.
[3] 胡川, 赵凯男, 黄修利, 吴金芝, 任开明, 王贺正, 付国占, 黄明, 李友军. 一次灌溉下耕作方式和氮肥用量对旱地小麦产量和品质的影响[J]. 作物学报, 2026, 52(6): 1830-1846.
[4] 马胜乾, 王志平, 陈浩天, 窦淑贤, 张燕, 邓艾兴, 张卫建, 原向阳, 宋振伟. 秸秆还田下耕作方式与氮肥施用量对东北玉米产量及土壤团聚体的影响[J]. 作物学报, 2026, 52(6): 1802-1816.
[5] 唐宽强, 李公允, 宋美毅, 赵雪, 常春玲. 大豆株高性状全基因组关联分析及预测模型构建[J]. 作物学报, 2026, 52(6): 1743-1756.
[6] 张思思, 赵向辉, 周洋, 姚云凤, 朱荣昱, 董元杰, 胡国庆, 徐通, 刘兆新. 冬闲期翻耕和绿肥还田对连作花生田土壤理化性质和产量的影响[J]. 作物学报, 2026, 52(5): 1472-1486.
[7] 张宁宁, 滕雨菲, 任娜娜, 魏兴卓, 闫书豪, 樊可心, 王永宏, 陈文康, 张兴华, 朱万超, 徐淑兔, 薛吉全. 201份玉米自交系抗旱表型评价及可塑性分析[J]. 作物学报, 2026, 52(5): 1309-1325.
[8] 王壮壮, 武紫君, 张永新, 张芯源, 袁丽雪, 陈如雪, 刘世举, 段剑钊, 冯伟, 王同朝, 王永华. 豫东南黏壤潮土区水氮优化协同提高冬小麦产量和氮素利用效率[J]. 作物学报, 2026, 52(5): 1501-1521.
[9] 杨欣雨, 崔文涛, 迪力尼格尔·阿力木, 汪凯翔, 吴鹏昊, 任姣姣. 玉米穗上叶片数全基因组关联分析和全基因组选择[J]. 作物学报, 2026, 52(5): 1573-1590.
[10] 张振, 冯连杰, 石玉, 于振文, 张永丽. 节水补灌下不同穗型小麦产量形成差异研究[J]. 作物学报, 2026, 52(5): 1522-1535.
[11] 闫安, 蒋昆炜, 王蓉圆, 田林, 张璐, 王韵, 徐建龙. 水稻剑叶小维管束数基因SVN7的鉴定与克隆[J]. 作物学报, 2026, 52(5): 1364-1372.
[12] 刘昕萌, 任昊, 张继波, 张吉旺, 赵斌, 任佰朝, 刘鹏, 王洪章. 茉莉酸甲酯(MeJA)缓解高温影响玉米雌穗分化的生理机制[J]. 作物学报, 2026, 52(5): 1561-1572.
[13] 张鸿蓉, 王菲儿, 李盼, 仇海龙, 朱静, 赵连豪, 南运有, 何蔚, 樊志龙, 胡发龙, 柴强, 殷文. 减量20%灌水与25%有机肥替代化肥提高青贮玉米产量的光合特性[J]. 作物学报, 2026, 52(5): 1487-1500.
[14] 王宇诚, 张露, 刘阿康, 黄见良, 彭少兵, 袁珅. 基于产量差的作物大面积单产提升策略与展望[J]. 作物学报, 2026, 52(5): 1279-1290.
[15] 赵佳雪, 周龙昊, 郭岂源, 尚伦霄, 王涵, 刘志涛, 陈曦, 张晓佩, 宋宪亮, 毛丽丽. 长期秸秆还田与深松通过改善土壤环境与棉花光合特性提高滨海盐碱地棉花产量[J]. 作物学报, 2026, 52(5): 1548-1560.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
京ICP备15008789号-2