郅晨陽 薛曉夢 吳 潔 李雄才 王 瑾 晏立英 王 欣 陳玉寧 康彥平 王志慧 淮東欣,* 洪彥彬 姜慧芳 雷 永,* 廖伯壽

1中國農(nóng)業(yè)科學院油料作物研究所 / 農(nóng)業(yè)農(nóng)村部遺傳育種重點實驗室, 湖北武漢 430062; 2 襄陽市農(nóng)業(yè)技術推廣中心, 湖北襄陽 441021; 3 河北省農(nóng)林科學院糧油作物研究所, 河北石家莊 050035; 4 廣東省農(nóng)業(yè)科學院作物研究所, 廣東廣州 510640

花生是世界范圍內(nèi)重要的油料作物和經(jīng)濟作物之一, 我國是世界花生生產(chǎn)和消費大國, 其產(chǎn)量和消費量持續(xù)增長[1]。近年來, 隨著人們消費需求的改變, 高品質食用花生在消費市場中嶄露頭角[2-3], 對其風味品質也提出了更高的要求。研究表明, 甜味與花生口味品質呈顯著正相關, 相關系數(shù)高達0.88[4], 甜味較高時, 花生口感較好?？扇苄蕴鞘腔ㄉ鹞兜闹饕獊碓? 其中蔗糖所占比例最高, 約占可溶性糖的90%, 與總糖含量極顯著正相關[5-6],當蔗糖含量高于5%, 花生籽仁可以品嘗到明顯的甜味[7], 表明蔗糖含量是影響花生甜味的關鍵因素,因此籽仁蔗糖含量是評價花生風味的關鍵指標。

研究發(fā)現(xiàn), 花生種質資源中蔗糖含量變異豐富,是可遺傳的數(shù)量性狀[7-9], 改良潛力大。已有研究表明, 花生蔗糖含量可能受兩對主基因+多基因調(diào)控,主要受母體效應控制[10-11]。通過定位研究, 分別在A06、A08和B06染色體上鑒定到調(diào)控花生蔗糖含量的主效QTL[12-14]。通過轉錄組-代謝組聯(lián)合分析,推測花生籽仁蔗糖含量可能來源于光合產(chǎn)物的非完全轉化[15]。雖然花生籽仁蔗糖含量研究已取得初步進展, 但其調(diào)控機制仍不明確, 因此, 開展花生籽仁蔗糖含量的遺傳研究分析, 探索其規(guī)律, 對指導花生蔗糖含量定位研究及花生品質育種具有重要意義。

本研究利用高糖品種(系) SYT5-1 (6.41%)和冀花甜1號(JHT1, 6.34%)分別構建了2個F2:3群體:NYBP×SYT5-1、19-1934×冀花甜1號, 分別獲得了467份和149份F2:3家系, 通過測定2個群體籽仁蔗糖含量, 分析蔗糖含量與含油量、蛋白質含量以及籽仁相關性狀之間的相關性, 利用SEA 2.0軟件對蔗糖含量進行了遺傳模型分析, 推測其遺傳規(guī)律,為進一步探索蔗糖含量調(diào)控機制、指導食用型花生品種培育奠定基礎。

1 材料與方法

1.1 群體材料

本研究以4個蔗糖含量存在極顯著差異的品種(系): SYT5-1 (6.41%)、南陽白皮(NYBP, 2.53%)、19-1934 (2.61%)和冀花甜1號(JHT1, 6.34%)為親本,配置2個雜交組合: NYBP×SYT5-1、19-1934×冀花甜1號, 種植于中國農(nóng)業(yè)科學院油料作物研究所陽邏基地, 按組合收獲F1、F2代種子, 分別獲得467份和149份F2:3家系的花生籽仁, 以其為試驗材料開展研究。其中, SYT5-1和19-1934分別為中國農(nóng)業(yè)科學院油料作物研究所選育的高糖品系和高油品系,南陽白皮為本單位保存的種質資源, 冀花甜1號由河北省農(nóng)林科學院育成并提供。

1.2 表型分析

1.2.1 花生籽仁蔗糖含量及蛋白質含量測定 利用近紅外模型測定花生籽仁蔗糖含量和蛋白質含量。每份材料選取15～20粒飽滿無病變的花生種子,將其放置于小樣品杯中, 利用Sperctra Star XL近紅外光譜儀(美國Unity科技公司)采集光譜并預測蔗糖含量和蛋白質含量, 每份材料重復3次, 分別計算蔗糖含量和蛋白含量的平均值。

1.2.2 花生籽仁含油量測定 利用脈沖核磁共振儀分析法, 使用核磁共振儀(PQ001)對花生籽仁含油量進行測定。每份材料選用10粒飽滿無病變的花生籽仁稱取其質量, 利用NMR (Nuclear Magnetic Resonance)技術得到的信號量與含油量的線性方程,計算樣品含油量, 每份材料重復3次, 以平均值作為該材料的含油量。

1.2.3 花生籽仁長、籽仁寬及百仁重測定 花生籽仁產(chǎn)量性狀統(tǒng)計: 選取15粒均勻飽滿的花生籽仁,利用萬深SC-G型自動考種儀及千粒重儀, 分別獲得籽仁長、籽仁寬和百仁重, 每份材料重復3次, 取平均值。

1.2.4 數(shù)據(jù)分析 利用SPSS 19統(tǒng)計F2:3群體花生蔗糖含量、含油量、蛋白質含量和籽仁性狀數(shù)據(jù),并計算各性狀間的相關系數(shù)。

應用SEA 2.0軟件程序[16]對花生蔗糖含量進行主基因+多基因遺傳模型分析。依據(jù)AIC準則和極大似然值篩選AIC值最小及極大似然值較大的遺傳模型作為備選模型, 并對獲得的所有模型進行適合性檢驗, 包括: 均勻性檢驗(U12、U22、U32)、Smirnov檢驗(nW2)、Kolmogorov檢驗(DN), 篩選適合性檢驗參數(shù)均未達到顯著或極顯著水平、具有其統(tǒng)計學意義的模型進行遺傳參數(shù)分析。使用最小二乘法對蔗糖含量遺傳模型進行遺傳參數(shù)估算確定最佳遺傳模型。

2 結果與分析

2.1 2個群體F2蔗糖含量頻數(shù)分布特征

對兩群體花生籽仁蔗糖含量頻數(shù)分布進行統(tǒng)計(圖1)發(fā)現(xiàn), 在NYBP×SYT5-1中, 蔗糖含量的變異范圍為0.19%～8.20%, 呈現(xiàn)連續(xù)分布的特點, 且群體蔗糖含量偏度與峰度絕對值均小于1 (附表1), 總體呈現(xiàn)近似正態(tài)分布, 表明蔗糖含量可能受多基因調(diào)控; 19-1934×JHT1群體中, 蔗糖含量變異范圍更大, 為0.09%～9.51% (附表1), 同樣呈現(xiàn)連續(xù)分布,但是表現(xiàn)明顯的雙峰分布趨勢, 表明蔗糖含量可能由少量主基因調(diào)控。在2個群體中都觀察到明顯的超親分離現(xiàn)象, 其中在NYBP×SYT5-1中蔗糖含量高于親本SYT5-1的家系共27個, 在19-1934×JHT1群體中蔗糖含量高于親本JHT1的家系共9個, 這些優(yōu)異家系為花生高糖育種提供了材料基礎。

附表1 F2:3群體蔗糖含量表型統(tǒng)計Table S1 Phenotypic statistics of sucrose content in F2:3 population

表1 花生籽仁蔗糖含量遺傳模型分析Table 1 Genetic model analysis of sucrose content in peanut seeds

圖1 花生籽仁蔗糖含量頻數(shù)分布圖Fig. 1 Frequency distribution of sucrose content in peanut kernel

2.2 花生籽仁蔗糖含量與含油量、蛋白質含量及產(chǎn)量相關性狀之間的相關性分析

由圖2可知, 在NYBP×SYT5-1群體中, 花生籽仁蔗糖含量與含油量呈極顯著負相關(r= -0.57**),與蛋白質含量顯著正相關(r= 0.52**), 但與籽仁長、籽仁寬和百粒重間無明顯相關關系。在19-1934 ×JHT1中, 籽仁蔗糖含量與含油量同樣顯著負相關(r= -0.80**), 與蛋白質含量也顯著正相關(r=0.75**),而與籽仁長、籽仁寬和百仁重間都表現(xiàn)為顯著負相關, 相關系數(shù)分別為-0.41**、-0.39**、-0.42**。同時發(fā)現(xiàn), 在2個群體中, 籽仁含油量與蛋白質含量顯著負相關, 相關系數(shù)在2個群體間分別為-0.47**、-0.62**。

圖2 花生籽仁蔗糖含量與其他性狀之間相關性Fig. 2 Correlations between sucrose content and other traits in peanut kernel

綜上, 在兩群體間, 籽仁蔗糖含量與含油量均呈顯著負相關, 與蛋白質含量均顯著正相關, 同時發(fā)現(xiàn), 籽仁含油量與蛋白質含量間也存在極顯著的負相關。但在2個群體中, 蔗糖含量與籽仁長、籽仁寬和百仁重3個性狀間的相關性不一致, 在19-1934 × JHT1中蔗糖含量與三者顯著負相關, 但在NYBP×SYT5-1群體中不存在相關關系。

2.3 花生籽仁蔗糖含量的遺傳模型分析

使用SEA 2.0軟件中SEA-F23程序對花生籽仁蔗糖含量進行主基因+多基因模型分析發(fā)現(xiàn), NYBP×SYT5-1與19-1934×JHT1群體均獲得11種遺傳模型(表1), 模型適合性檢驗結果顯示, 所有模型均為達到顯著或極顯著水平, 具有其統(tǒng)計學意義(表2)。

表2 花生籽仁蔗糖含量候選模型的適合性檢驗Table 2 Suitability test of candidate models for sucrose content in peanut seeds

依據(jù)極大似然值最大及AIC值最小原則對2個組合間所有遺傳模型進行篩選(表1)發(fā)現(xiàn), 在NYBP×SYT5-1群體中, 2MG-AD極大似然值為–847.85,AIC值為1707.69, 是最優(yōu)候選模型; 在19-1934×JHT1中, 2MG-AD同樣具有最大的極大似然值及最小AIC值, 分別為-298.02與608.03, 是最優(yōu)候選模型。綜上, 兩群體最優(yōu)候選模型均為2MG-AD, 即兩對主基因控制加性顯性效應是最優(yōu)遺傳模型。

一階參數(shù)結果表明(表3), 在NYBP×SYT5-1群體中, 第1對主基因與第2對主基因加性效應分別為da=1.6676、db=1.2171, 顯性效應分別為ha=-0.6996、hb= -0.3973; 19-1934×JHT1群體中, 兩對主基因的加性效應分別為da=1.7948、db=0.7131, 顯性效應分別為ha=1.9106、hb= -0.1276。

表3 花生籽仁蔗糖含量一階參數(shù)估計值Table 3 Estimates of first order genetic parameters for sugar content in peanut kernels

兩群體間, 兩對主基因的加性效應均表現(xiàn)為|da|>|db|, 即表明兩對主基因調(diào)控蔗糖含量表型間,以其中一對主基因的加性效應為主; 顯性效應|ha|>|hb|, 表明組合中控制籽仁蔗糖含量的其中一對主基因顯性效應大于另一對主基因, 顯性效應更加明顯。因此, 調(diào)控花生籽仁蔗糖含量的主基因間存在明顯的互作, 且其中一對主基因的加性效應和顯性效應大于另一對主基因, 對控制花生籽仁蔗糖含量的效應更加明顯。

二階遺傳參數(shù)結果顯示(表4), NYBP×SYT5-1中, 最優(yōu)模型2MG-AD的主基因方差與主基因遺傳率最高, 分別為2.2810和85.9102; 19-1934×JHT1組合中, 主基因方差與遺傳率為4.0724和97.7852。在兩群體間蔗糖含量的主基因遺傳率很高, 表明蔗糖含量主要受遺傳效應控制, 受環(huán)境影響較小。

表4 花生籽仁蔗糖含量二階遺傳參數(shù)估計值Table 4 Estimated values of genetic parameters of sucrose content in peanut kernels

3 討論

在我國, 花生是主要的油料和經(jīng)濟作物之一,我國花生總產(chǎn)的約40%用于食用加工, 食用花生的消費量也在逐年增長[17]。蔗糖含量是影響花生及其加工食品風味的重要因素, 已成為食用花生品質遺傳改良的重要方向, 研究蔗糖含量的遺傳規(guī)律將為開展高糖花生育種奠定重要的理論基礎。本研究利用2個高糖材料SYT5-1和JHT1, 分別構建了F2:3群體, 雜交后代中蔗糖含量呈連續(xù)分布, 變異范圍分別為0.19%～8.20%和0.09%～9.06% (圖1), 與Bishi等在60份珍珠豆型資源材料及41份印度主栽培品種中檢測的蔗糖含量(2.44%～7.61%和2.61%～6.50%)相比[7-8], 雜交群體中蔗糖變異范圍更大, 創(chuàng)造了更多的變異。同時發(fā)現(xiàn), 本研究2個群體間蔗糖含量表現(xiàn)出顯著不同的連續(xù)分布特征, NYBP×SYT5-1群體蔗糖含量呈正態(tài)分布(圖1-A), 但19-1934×JHT1群體的雙峰分布趨勢明顯(圖1-B)。對比前人研究結果,同樣發(fā)現(xiàn)不同材料構建的群體間, 蔗糖含量的分布特征存在差異, NYBP×SYT5-1群體與張曉軍等[12]構建的分離群體中, 蔗糖含量均呈近似正態(tài)分布, 表明花生蔗糖含量在這些群體中可能受多基因調(diào)控[18];而本研究19-1934×JHT1群體與秦利等[11]、Guo等[14]的試驗結果均表現(xiàn)出明顯的雙峰分布, 可能受少量主基因調(diào)控。造成這種差異的原因可能與選用的雜交親本攜帶不同的蔗糖含量主效位點相關, 加上不同組合群體間親本材料遺傳背景存在差異, 導致不同群體間蔗糖含量分布特征存在差異。

本研究的2個雜交分離群體中, 花生蔗糖含量與含油量均呈極顯著負相關(r= -0.57**,r= –0.80**)(圖2)。在大豆、油菜等油料作物中, 蔗糖含量與含油量也表現(xiàn)出類似的負相關性[19-21], 這種負相關性可能與油料作物籽仁中的油脂來源于蔗糖相關, 多項研究結果表明, 油料作物的光合產(chǎn)物以蔗糖形式運輸?shù)阶讶? 并被轉化分解, 為脂肪酸和三酰甘油的合成積累提供最主要碳源[22-24]。在油菜種子胚中,隨著胚發(fā)育, 蔗糖含量逐步降低, 而油脂含量不斷積累增加[24]。同時, 蔗糖還作為一種重要的信號物質調(diào)控糖酵解途徑, 影響該途徑中多種基因的表達,對下游油脂合成途徑產(chǎn)生影響[25]。蓖麻種子中糖代謝相關基因表達量下調(diào)引起油脂積累相關基因表達量的明顯提高[26]。在以淀粉作為主要儲藏物質的玉米、水稻等作物相關研究中, 與糖代謝相關的基因的突變造成糖分積累并引起淀粉含量的下降或淀粉性質的變化[27-29],OsSWEET11突變體中糖-淀粉代謝途徑明顯變化, 對水稻穎果發(fā)育產(chǎn)生了重要影響[30]。因此推測, 花生籽仁中蔗糖含量的積累可能是油脂合成受阻引起。

在本研究兩群體中, 花生籽仁蔗糖與蛋白質含量均呈顯著正相關(r= 0.52**,r= 0.75**), 含油量與蛋白質含量顯著負相關(r= –0.47**,r= –0.62**, 圖2),大豆相關研究存在相似結論[31]。籽仁中油脂和蛋白質的合成主要以蔗糖的分解產(chǎn)物丙酮酸為共同底物[32], 二者存在底物競爭關系, 油脂含量與蛋白質含量通常表現(xiàn)為顯著負相關[33-35]。因此推測, 花生籽仁內(nèi)蔗糖含量與蛋白質含量顯著正相關可能與油脂合成受阻、碳源流向蛋白質合成有關。

在19-1934×JHT1群體中, 蔗糖含量與籽仁長、籽仁寬、百仁重顯著負相關(圖2-B), 這與油菜、大豆中的研究結果相似[36-38], 籽仁中蔗糖含量變化可以引起每角果粒數(shù)(莢粒數(shù))和粒重等相關性狀的變化。蔗糖含量高低總體體現(xiàn)了植物體內(nèi)碳源供應水平, 除了為種子油脂積累提供所需的碳源, 蔗糖也為胚組織細胞的發(fā)育和分裂提供養(yǎng)分[39-40], 進而影響了籽仁發(fā)育。在玉米、水稻、大豆等研究中表明,蔗糖轉運蛋白SWEET相關基因在作物馴化過程中受到選擇, 對產(chǎn)量和品質都具有重要影響, SWEET基因的突變造成蔗糖轉運受阻, 引起胚乳發(fā)育不良,籽粒大小及粒重均顯著降低[27,41]。膜轉運蛋白MFS家族也可通過蛋白互作方式影響水稻光合產(chǎn)物的分布平衡, 調(diào)控灌漿、籽粒大小及粒數(shù)等重要性狀[42]。在本研究2個群體中, 蔗糖含量與籽仁大小相關性狀的相關性不一致, 也可能與親本材料具有不同的遺傳背景相關。

本研究遺傳分析結果表明, 蔗糖含量具有很高的遺傳率, 分別為85.91%與97.79% (表4), 與秦利等[11]、胡美玲等[10]研究結果一致?；ㄉ讶收崽呛窟z傳率高, 表明其主要受基因型調(diào)控, 環(huán)境因素對其影響較小, 因此在育種群體低世代即可進行蔗糖含量的選擇, 以加快育種進程, 提高育種效率。遺傳模型分析表明, 花生籽仁蔗糖含量主要受兩對主基因的加性顯性效應調(diào)控(表3), 與秦利等[11]研究結果一致。兩對主基因的加性效應|da|>|db|、顯性效應|ha|>|hb|, 表明主基因間存在明顯互作, 且以其中一對主基因的加性顯性效應為主。同時, 兩群體中加性效應da>db>0, 兩對主基因的加性效應均為正效應,即基因間的加性效應可以增強, 具有累加效應?；ㄉ讶收崽呛窟z傳率高且受兩對主基因控制, 表明蔗糖含量可以通過開發(fā)與蔗糖含量相關的分子標記進行快速選擇育種, 利用主基因的累加效應開展超親育種, 進行種質創(chuàng)新和高蔗糖品種育種具有充分的可行性, 為培育不同蔗糖含量的食用花生品種提供了豐富的親本材料和選擇空間, 對花生籽仁蔗糖含量的改良具有重要指導意義。

4 結論

本研究利用高糖材料SYT5-1和JHT1分別構建了2個F2:3群體, 花生籽仁蔗糖含量在兩群體中都呈現(xiàn)連續(xù)分布特點, 超親分離現(xiàn)象明顯。在兩群體中, 花生籽仁的蔗糖含量均與含油量顯著負相關,與蛋白含量顯著正相關; 但與籽仁長、籽仁寬和百仁重的相關性不一致?；ㄉ崽呛恐饕軆蓪χ骰蚩刂? 兩對主基因間存在明顯互作, 累加效應明顯, 且其中一對主基因的加性效應和顯性效應大于另一對主基因, 對花生蔗糖含量的調(diào)控效果更加顯著。