甘藍(lán)型油菜每角粒數(shù)的全基因組關(guān)聯(lián)分析

孫程明 陳 鋒 陳 松 彭 琦 張 維 易 斌張潔夫,* 傅廷棟

研究簡報

甘藍(lán)型油菜每角粒數(shù)的全基因組關(guān)聯(lián)分析

孫程明1,2陳 鋒1陳 松1彭 琦1張 維1易 斌2,*張潔夫1,*傅廷棟2

1江蘇省農(nóng)業(yè)科學(xué)院經(jīng)濟(jì)作物研究所 / 農(nóng)業(yè)部長江下游棉花與油菜重點實驗室 / 江蘇省現(xiàn)代作物生產(chǎn)協(xié)同創(chuàng)新中心, 江蘇南京210014;2華中農(nóng)業(yè)大學(xué)植物科學(xué)技術(shù)學(xué)院 / 作物遺傳改良國家重點實驗室, 湖北武漢 430070

每角粒數(shù)是油菜重要的產(chǎn)量構(gòu)成因子, 增加每角粒數(shù)有助于提高油菜的籽粒產(chǎn)量。利用Illumina 60K SNP芯片對496份具有代表性的油菜資源進(jìn)行基因型分析, 考察該群體在2個環(huán)境中的每角粒數(shù), 利用MLM和GLM模型進(jìn)行全基因組關(guān)聯(lián)分析。結(jié)果表明, 本群體在2個環(huán)境中每角粒數(shù)的廣義遺傳力為57.7%。利用MLM和GLM模型分別檢測到9個和20個位點, 所有MLM位點均得到GLM結(jié)果的驗證。6個位點與前人定位的QTL重疊, 其中2個位點得到2次驗證, 其余14個是新位點。在7個位點附近找到了候選基因, 其中在C09染色體的位點Bn-scaff_15576_1-p74980附近找到已克隆的油菜每角粒數(shù)基因, 在其余6個位點附近找到、等已知的擬南芥每角粒數(shù)基因的同源基因。本研究結(jié)果有助于解析油菜每角粒數(shù)的遺傳基礎(chǔ)及其調(diào)控機(jī)制, 為每角粒數(shù)的遺傳改良奠定了基礎(chǔ)。

甘藍(lán)型油菜; 產(chǎn)量; 每角粒數(shù); 關(guān)聯(lián)分析; SNP標(biāo)記

油菜是我國重要的油料作物, 種植面積約690萬公頃,產(chǎn)油量占國產(chǎn)植物油總量的55%以上, 是我國最大的食用植物油來源[1]。作為食用油消費(fèi)大國, 我國食用油60%以上卻依賴進(jìn)口, 油菜的產(chǎn)量水平直接影響我國食用油的供給, 因此, 提高產(chǎn)量成為油菜育種最迫切的目標(biāo)。每角粒數(shù)是油菜重要的產(chǎn)量構(gòu)成因子, 解析其遺傳基礎(chǔ)對提高油菜產(chǎn)量具有重要意義。

油菜種子的形成起始于胚珠原基, 而胚珠原基需要經(jīng)歷珠被、珠心和胚囊的形成, 以及授粉、受精、種子發(fā)育等過程才能最終形成種子, 這一系列過程受到復(fù)雜的遺傳網(wǎng)絡(luò)調(diào)控。李永鵬等[2]研究表明, 部分胚珠原基在發(fā)育過程中敗育而最終不能形成成熟的種子, 正常發(fā)育胚珠百分比是反映品種間每角粒數(shù)差異的重要因素。Yang等[3]進(jìn)一步研究表明, 油菜每角粒數(shù)的自然變異由每子房胚珠數(shù)(30.7%)、可育胚珠比(18.2%)、可育胚珠授精率(7.1%)和授精胚珠發(fā)育成種子的比率(43.9%)的差異共同決定。

多位研究者采用連鎖分析定位了每角粒數(shù)QTL, Shi等[4]利用TN DH及F2群體在13條染色體檢測到69個每角粒數(shù)QTL, 解釋的表型變異范圍是3.2%~15.5%; Luo等[5]利用60K SNP芯片再次對TN DH群體進(jìn)行QTL定位, 共檢測到64個每角粒數(shù)QTL, 其中貢獻(xiàn)率最大的2個QTL和均位于A01染色體, 分別解釋20.4%和24.1%的表型變異; 漆麗萍[6]在A01、A06~A08、C01、C03和C05染色體檢測到15個QTL, 其中貢獻(xiàn)率最大的2個QTL和均位于A01染色體, 分別解釋26.3%和26.6%的表型變異; Cai等[7]在A01~A04、A06、A08~A10、C01、C03、C05、C06和C08染色體檢測到17個QTL, 可解釋5.6%~26.4%的表型變異; Yang等[8]利用一個RIL群體在A02、A06、C01、C04和C06染色體檢測到5個QTL, 可解釋5.3%~24.8%的表型變異, 其中位于A06的QTL在多個環(huán)境中被重復(fù)檢測到, 解釋20.1%的表型變異。

近年來, 隨著高密度SNP基因分型芯片和全基因組測序等技術(shù)的發(fā)展, 全基因組關(guān)聯(lián)分析(genome-wide association study, GWAS)已廣泛應(yīng)用于油菜復(fù)雜性狀的遺傳結(jié)構(gòu)解析[9-12]。本研究以496份具有代表性的油菜資源為關(guān)聯(lián)群體, 利用60K SNP芯片對群體進(jìn)行基因型分析, 并對每角粒數(shù)進(jìn)行全基因組關(guān)聯(lián)分析, 旨在挖掘顯著的SNP位點, 分析其候選基因, 為油菜每角粒數(shù)的遺傳改良奠定基礎(chǔ)。

1 材料與方法

1.1 試驗材料

用于關(guān)聯(lián)分析的496份甘藍(lán)型油菜包括國內(nèi)外的地方品種、育成品種及高世代育種材料(附表1)。其中國內(nèi)資源444份, 主要來自湖北、重慶、江蘇、湖南、四川、陜西等油菜主產(chǎn)省市; 國外資源52份, 主要來自德國、瑞典、朝鮮、加拿大等國家。所有材料均由華中農(nóng)業(yè)大學(xué)國家油菜工程技術(shù)研究中心提供。

1.2 田間試驗與表型調(diào)查

于2015年和2016年在江蘇泰州(15TZ和16TZ)種植自然群體。采用完全隨機(jī)區(qū)組設(shè)計, 設(shè)置2個重復(fù), 單個材料每重復(fù)種2行, 每行15株, 行寬1.5 m, 行距40 cm。每年10月上旬大田直播, 田間管理按當(dāng)?shù)爻Ｒ?guī)方式進(jìn)行。在油菜成熟期, 取每小區(qū)8株長勢一致的植株, 選主花序中部10個角果測量每角粒數(shù)。利用R軟件對各環(huán)境的每角粒數(shù)表型進(jìn)行統(tǒng)計分析和相關(guān)性分析, 并計算廣義遺傳力和最佳線性無偏預(yù)測值(best linear unbiased prediction, BLUP)[13-14]。

1.3 基因型檢測與分析

利用Illumina 60K SNP芯片對496份甘藍(lán)型油菜資源進(jìn)行基因型分析。在Genome Studio軟件中, 剔除Call Freq參數(shù) < 0.75, 稀有等位基因頻率 < 0.05, AA、BB頻率 < 0.03及GenTrain值 < 0.50的SNP, 基因型雜合的SNP按缺失值處理。將過濾后的33,218個SNP序列比對至油菜Darmor-基因組, e-value閾值設(shè)為10-12, 保留19,167個單拷貝、位置清楚的SNP標(biāo)記用于后續(xù)分析。

1.4 全基因組關(guān)聯(lián)分析

用Structure 2.3.3進(jìn)行群體結(jié)構(gòu)Q矩陣的計算[15], 用SPAGeDi計算材料間的親緣關(guān)系K矩陣[16]。將基因型、表型、Q矩陣和K矩陣導(dǎo)入Tassel 4.0后, 以Q矩陣和Q+K矩陣作協(xié)變量, 分別進(jìn)行基于一般線性模型(General Linear Model, GLM)和混合線性模型(Mixed Linear Model, MLM)的全基因組關(guān)聯(lián)分析[17]。參考前人文獻(xiàn)[9-12], 關(guān)聯(lián)分析的顯著性閾值 (Bonferroni threshold) 定為1/總標(biāo)記數(shù), 即-lg() = 4.28。當(dāng)1 Mb區(qū)間內(nèi)存在多個顯著SNP時, 如果兩兩間的2≥0.1, 則將這些SNP歸為一個關(guān)聯(lián)位點, 以最小值的SNP作為代表。用R軟件包qqman繪制曼哈頓圖和QQ (Quantile- Quantile)圖[18]。

1.5 候選基因挖掘

以2= 0.1作為衰減閾值, 用Tassel計算顯著關(guān)聯(lián)位點的LD衰減距離作為其置信區(qū)間, 提取置信區(qū)間內(nèi)的基因CDS序列, 與擬南芥的基因CDS進(jìn)行BLAST比對, 設(shè)e-value閾值為10-10, 以相似度最高的擬南芥基因信息來注釋油菜基因。將落在置信區(qū)間的油菜每角粒數(shù)基因和擬南芥每角粒數(shù)基因的同源基因作為候選基因。

2 結(jié)果與分析

2.1 每角粒數(shù)表型分析

496份油菜資源的每角粒數(shù)在2個環(huán)境中具有廣泛的變異, 單個環(huán)境的每角粒數(shù)極差范圍是15.24~18.90。單個環(huán)境表型平均值范圍是(20.56±1.99)~(21.45±2.36), 變異系數(shù)范圍是0.10~0.11 (表1和圖1)。2個環(huán)境間的相關(guān)系數(shù)為0.55 (≤ 0.001), 達(dá)到極顯著水平, 這表明本研究考察的表型數(shù)據(jù)具有較高的可靠性和重復(fù)性。每角粒數(shù)在2個環(huán)境的廣義遺傳力為57.7%。

2.2 每角粒數(shù)全基因組關(guān)聯(lián)分析

利用MLM模型對各環(huán)境的每角粒數(shù)表型和BLUP值進(jìn)行全基因組關(guān)聯(lián)分析, 當(dāng)閾值為?lg() = 4.28時, 在15TZ、16TZ和基于BLUP分別檢測到6、2和2個顯著位點(表2和圖2)。位于C02的位點Bn-scaff_15712_6- p1336179在15TZ和16TZ兩環(huán)境中被重復(fù)檢測到。合并重疊的位點后得到9個顯著位點, 分布在A01、A04、A07、A08、C01~C04染色體, 可解釋3.26%~4.61%的表型變異。

表1 關(guān)聯(lián)群體每角粒數(shù)性狀的統(tǒng)計分析

圖1 關(guān)聯(lián)群體在2個環(huán)境的每角粒數(shù)分布

圖2 油菜每角粒數(shù)全基因組關(guān)聯(lián)分析(MLM)

A: 每角粒數(shù)BLUP值MLM曼哈頓圖; B: 15TZ每角粒數(shù)MLM曼哈頓圖; C: 16TZ每角粒數(shù)MLM曼哈頓圖。水平線代表Bonferroni閾值。

A: Manhattan plot of MLM based on BLUP value; B: Manhattan plot of MLM in15TZ; C: Manhattan plot of MLM in 16TZ. The dashed horizontal line depicts the Bonferroni significance threshold.

表2 MLM每角粒數(shù)顯著關(guān)聯(lián)位點

利用GLM模型對各環(huán)境的每角粒數(shù)表型和BLUP值進(jìn)行全基因組關(guān)聯(lián)分析, 當(dāng)閾值為?lg() = 4.28時, 在15TZ、16TZ和基于BLUP分別檢測到10、6和5個顯著位點。位于C04的位點Bn-scaff_23907_1-p3780基于BLUP和在15TZ均被重復(fù)檢測到, 合并重疊的位點后得到20個顯著位點, 分布在A01、A04、A07、A08、C01~C07和C09染色體, 可解釋2.90%~4.83%的表型變異(表3和圖3)。比較2個關(guān)聯(lián)模型的結(jié)果, 所有MLM位點均得到GLM結(jié)果的驗證, 合并2個模型的結(jié)果共得到20個顯著位點。

圖3 油菜每角粒數(shù)全基因組關(guān)聯(lián)分析(GLM)

A: 每角粒數(shù)BLUP值GLM曼哈頓圖; B: 15TZ每角粒數(shù)GLM曼哈頓圖; C: 16TZ每角粒數(shù)GLM曼哈頓圖。水平線代表Bonferroni閾值。

A: Manhattan plot of GLM based on BLUP value; B: Manhattan plot of GLM in 15TZ; C: Manhattan plot of GLM in 16TZ. The dashed horizontal line depicts the Bonferroni significance threshold.

表3 GLM每角粒數(shù)顯著關(guān)聯(lián)位點

2.3 候選基因挖掘

基于油菜基因組注釋信息, 在C09染色體Bn-scaff_ 15576_1-p74980位點下游82 kb找到已克隆的油菜每角粒數(shù)基因, 該基因是擬南芥的同源基因(表4)。參與調(diào)控油菜大孢子母細(xì)胞的減數(shù)分裂, 該基因的自然缺失導(dǎo)致部分雌配子體發(fā)育異常, 無法形成正常功能的胚珠, 最終導(dǎo)致油菜每角果粒數(shù)顯著降低[19]。此外, 在A04染色體Bn-A04-p10393460位點上游522 kb找到候選基因, 該基因的擬南芥同源基因編碼一個短富含甘氨酸結(jié)構(gòu)域的蛋白, 該基因突變導(dǎo)致植株角果長度縮短、每角粒數(shù)減少、粒重降低[20]。在C03染色體Bn-scaff_23954_1-p635109位點下游162 kb找到候選基因, 該基因的擬南芥同源基因編碼一個E3泛素連接酶, 其功能缺失突變體每角粒數(shù)增加、粒重增加[21]。此外, 本研究還找到其他候選基因, 包括擬南芥已知每角粒數(shù)調(diào)控基因[22]、[23][24]和[25]在油菜中的同源拷貝。

3 討論

每角粒數(shù)是油菜重要的產(chǎn)量構(gòu)成因子, 增加每角粒數(shù)有利于提高單株籽粒產(chǎn)量, 從而提高群體的籽?？偭俊１狙芯繉?96份油菜資源的每角粒數(shù)進(jìn)行2個環(huán)境的考察, 分析表明本群體的每角粒數(shù)廣義遺傳力為57.7%。利用2個關(guān)聯(lián)模型MLM和GLM, 檢測到20個顯著位點, 分布于12條染色體。本研究中每角粒數(shù)關(guān)聯(lián)位點的效應(yīng)較小(< 5%), 在2個環(huán)境間能被重復(fù)檢測到的位點少, 表明該群體中每角粒數(shù)主要受微效多基因調(diào)控, 且基因效應(yīng)受環(huán)境影響較大。

本研究獲得的每角粒數(shù)QTL有6個位點與前人結(jié)果一致(表2和表3), 其中位于C01染色體的位點Bn-scaff_15936_1-p270915, 物理位置36.41 Mb, 與Luo等[5]檢測到的每角粒數(shù)QTL(范圍: 36.53~37.48 Mb)和漆麗萍[6]檢測到的QTL(范圍: 35.69~36.81 Mb)重疊。位于C09染色體的位點Bn-scaff_15576_ 1-p74980物理位置是41.13 Mb, 與Li等[19]已克隆的每角粒數(shù)基因(位置: 41.21 Mb)和Luo等[5]檢測到的QTL(范圍: 41.28~41.67 Mb)重疊。這些結(jié)果證實了本研究結(jié)果的可靠性。此外, 本研究還檢測到14個新位點, 如位于A04染色體的Bn-A04-p10393460和C04染色體的Bn-scaff_23907_1-p3780, 表型貢獻(xiàn)率均相對較高; 位于C02染色體的位點Bn-scaff_15712_6- p1336179在2個環(huán)境中被重復(fù)檢測到, 這些位點值得進(jìn)一步驗證和研究。

目前油菜中已報道的每角粒數(shù)基因較少, Li等[19]利用圖位克隆的方法, 在C09染色體克隆了每角粒數(shù)基因; Yang等[26]應(yīng)用CRISPR/Cas9系統(tǒng)敲除在油菜A04和C04染色體的2個拷貝, 雙突變植株葉片增多, 每角粒數(shù)和粒重均增加; Shah等[27]發(fā)現(xiàn)擬南芥在油菜A02染色體的拷貝突變, 導(dǎo)致了油菜株高、分枝高度、分枝數(shù)、每角粒數(shù)和單株產(chǎn)量等性狀的改變。本研究在Bn-scaff_15576_1-p74980位點下游82 kb找到(表4)。此外, 還在另外6個位點附近找到了候選基因, 這些基因的擬南芥同源基因均參與每角粒數(shù)的調(diào)控, 部分基因還有一因多效的作用, 如[20]和[22]同時影響角果長、粒型和粒重等性狀。后續(xù)研究可通過轉(zhuǎn)基因或CRISPR/Cas9敲除, 進(jìn)一步驗證這些候選基因的功能。

表4 每角粒數(shù)關(guān)聯(lián)位點候選基因信息

附表 請見網(wǎng)絡(luò)版: 1) 本刊網(wǎng)站http://zwxb. chinacrops.org/; 2) 中國知網(wǎng)http://www.cnki.net/; 3) 萬方數(shù)據(jù)http://c.wanfangdata.com.cn/Periodical-zuowxb.aspx。

[1] 王漢中. 我國油菜產(chǎn)業(yè)發(fā)展的歷史回顧與展望. 中國油料作物學(xué)報, 2010, 32: 300–302. Wang H Z. Review and future development of rapeseed industry in China., 2010, 32: 300–302 (in Chinese with English abstract).

[2] 李永鵬, 程焱, 蔡光勤, 范楚川, 周永明. 油菜每角果粒數(shù)差異的細(xì)胞學(xué)基礎(chǔ)和分子機(jī)理. 中國科學(xué): 生命科學(xué), 2014, 44: 822–831. Li Y P, Cheng Y, Cai G Q, Fan C C, Zhou Y M. Cytological basis and molecular mechanism of variation in number of seeds per pod in., 2014, 44: 822–831 (in Chinese with English abstract).

[3] Yang Y, Wang Y, Zhan J, Shi J, Wang X, Liu G, Wang H. Genetic and cytological analyses of the natural variation of seed number per pod in rapeseed (L.)., 2017, 8: 1890. doi: 10.3389/fpls.2017.01890.

[4] Shi J, Li R, Qiu D, Jiang C, Long Y, Morgan C, Bancroft I, Zhao J, Meng J. Unraveling the complex trait of crop yield with quantitative trait loci mapping in., 2009, 182: 851–861.

[5] Luo Z, Wang M, Long Y, Huang Y, Shi L, Zhang C, Liu X, Fitt B D, Xiang J, Mason A S. Incorporating pleiotropic quantitative trait loci in dissection of complex traits: seed yield in rapeseed as an example., 2017, 130: 1569–1585.

[6] 漆麗萍. 甘藍(lán)型油菜株型與角果相關(guān)性狀的QTL分析. 華中農(nóng)業(yè)大學(xué)博士學(xué)位論文, 湖北武漢, 2014. Qi L P. QTL Analysis for the Traits Associated with Plant Architecture and Silique inL. PhD Dissertation of Huazhong Agricultural University, Wuhan, Hubei, China, 2014 (in Chinese with English abstract).

[7] Cai G, Yang Q, Chen H, Yang Q, Zhang C, Fan C, Zhou Y. Genetic dissection of plant architecture and yield-related traits in, 2016, 6: 21625. doi: 10.1038/srep 21625.

[8] Yang Y, Shi J, Wang X, Liu G, Wang H. Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis., 2016, 6: 24124. doi: 10.1038/srep24124.

[9] Sun C, Wang B, Yan L, Hu K, Liu S, Zhou Y, Guan C, Zhang Z, Li J, Zhang J, Chen S, Wen J, Ma C, Tu J, Shen J, Fu T, Yi B. Genome-wide association study provides insight into the genetic control of plant height in rapeseed (L.)., 2016, 7: 1102. doi: 10.3389/fpls.2016.01102.

[10] Lu K, Wei L, Li X, Wang Y, Wu J, Liu M, Zhang C, Chen Z, Xiao Z, Jian H. Whole-genome resequencing revealsorigin and genetic loci involved in its improvement., 2019, 10: 1154. doi: 10.1038/s41467-019-09134-9.

[11] Chen L, Wan H, Qian J, Guo J, Sun C, Wen J, Yi B, Ma C, Tu J, Song L. Genome-wide association study of cadmium accumulation at the seedling stage in rapeseed (L.)., 2018, 9: 375. doi: 10.3389/fpls.2018.00375.

[12] Xu L, Hu K, Zhang Z, Guan C, Chen S, Hua W, Li J, Wen J, Yi B, Shen J. Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (L.)., 2015, 23: 43–52.

[13] 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., 2012, 137: 427–437.

[14] Ihaka R, Gentleman R. R: a language for data analysis and graphics., 1996, 5: 299–314.

[15] Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study., 2005, 14: 2611–2620.

[16] Hardy O J, Vekemans X. SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels., 2002, 2: 618–620.

[17] 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., 2007, 23: 2633–2635.

[18] Turner S D. qqman: an R package for visualizing GWAS results using QQ and manhattan plots., 2014, 1: 005165. doi: http://dx.doi.org/10.1101/005165.

[19] Li S, Chen L, Zhang L, Li X, Liu Y, Wu Z, Dong F, Wan L, Liu K, Hong D. BnaC9. SMG7b functions as a positive regulator of the number of seeds per silique inby regulating the formation of functional female gametophytes., 2015, 169: 2744–2760.

[20] Rodríguez-Hernández A A, Muro-Medina C V, Ramírez-Alonso J I, Jiménez-Bremont J F. Modification of AtGRDP1 gene expression affects silique and seed development in., 2017, 486: 252–256.

[21] Xia T, Li N, Dumenil J, Li J, Kamenski A, Bevan M W, Gao F, Li Y. The ubiquitin receptor DA1 interacts with the E3 ubiquitin ligase DA2 to regulate seed and organ size in Arabidopsis., 2013, 25: 3347–3359.

[22] Braud C, Zheng W, Xiao W. LONO1 encoding a nucleoporin is required for embryogenesis and seed viability in Arabidopsis., 2012, 160: 823–836.

[23] Groszmann M, Paicu T, Smyth D R. Functional domains of SPATULA, a bHLH transcription factor involved in carpel and fruit development in Arabidopsis., 2008, 55: 40–52.

[24] Chen N N, Chen H R, Yeh S Y, Vittore G, Ho H D. Autophagy is enhanced and floral development is impaired in AtHVA22d RNA interference Arabidopsis., 2009, 149: 1679–1689.

[25] Leroy O, Hennig L, Breuninger H, Laux T, K?hler C. Polycomb group proteins function in the female gametophyte to determine seed development in plants., 2007, 134: 3639–3648.

[26] Yang Y, Zhu K, Li H, Han S, Meng Q, Khan S U, Fan C, Xie K, Zhou Y. Precise editing of CLAVATA genes inL. regulates multilocular silique development., 2018, 16: 1322–1335.

[27] Shah S, Karunarathna N L, Jung C, Emrani N. An APETALA1 ortholog affects plant architecture and seed yield component in oilseed rape (L.)., 2018, 18: 380. doi: 10.1186/s12870-018-1606-9.

Genome-wide association study of seed number per silique in rapeseed (L.)

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

1Institute 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;2National Key Laboratory of Crop Genetic Improvement / College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, Hubei, China

Seed number per silique (SSN) is a key component of seed yield in rapeseed, increasing SSN can improve the seed yield of plants. A collection of 496 representative rapeseed accessions was genotyped by the Illumina 60K SNP array and phenotyped for SSN in two environments. The genome-wide association study (GWAS) of SSN was performed via the MLM (Mixed linear model) and GLM (General linear model). The broad-sense heritability of SSN was 57.7%. Nine and twenty loci were detected with MLM and GLM, respectively, and all loci detected by MLM were included those by GLM. Six loci were overlapped with reported QTLs, and two of them were validated by two independent researches, and the rest 14 loci were new. We identified plausible candidate genes nearby seven loci, and the reported rapeseed SSN genewas found near the locus Bn-scaff_15576_1-p74980 on C09 chromosome detected in this study. Besides, six candidates orthologous to documentedSSN genes, like,,,and, were found near our GWAS loci. The results provide an insight into the genetic basis of seed number per silique and lay a foundation for further mechanism exploration and breeding for this trait in.

L.; yield; seed number per silique; GWAS; SNP

2019-04-15;

2019-08-09;

2019-09-11.

10.3724/SP.J.1006.2020.94060

本研究由國家重點研發(fā)計劃項目(2018YFD0100602), 國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(CARS-12), 江蘇省農(nóng)業(yè)科技自主創(chuàng)新基金(CX(19)3055), 江蘇省基礎(chǔ)研究計劃(自然科學(xué)基金)項目(BK20190260)和中央高?；究蒲袠I(yè)務(wù)費(fèi)專項資金資助項(2662016PY063)資助。

The study was supported by the National Key Research and Development Program of China (2018YFD0100602), the Earmarked Fund for China Agriculture Research System (CARS-12), the Jiangsu Agriculture Science and Technology Innovation Fund (CX(19)3055), the Natural Fund Project of Jiangsu Basic Research Program (BK20190260), and the Fundamental Research Funds for the Central Universities (2662016PY063).

張潔夫, E-mail: jiefu_z@163.com; 易斌, E-mail:yibin@mail.hzau.edu.cn

E-mail: suncm8331537@gmail.com

URL:http://kns.cnki.net/kcms/detail/11.1809.S.20190910.1542.018.htm