朱秀瑾 郭鳳丹 夏晗 趙傳志 侯蕾



摘要:花生突變體ssm1 與其對照親本魯花11 號(LH11)相比,種子皺縮變小,百仁重、單株生產力和含油量顯著降低。 為了探究導致該突變體種子變小且含油量降低的原因,本試驗研究了種子發育關鍵轉錄因子基因LCE1 等在突變體ssm1 及其對照LH11 合子受精及種子發育不同階段的表達情況。 結果表明,AhLEC1 在花生未受精子房(DBP0)時期表達量較低,且在ssm1 和對照LH11 中無明顯差異,在受精后的子房(DBP1)和入土10 d(DAP10)和20 d(DAP20)種子中表達量上調,并且在LH11 中的表達量極顯著高于突變體ssm1、轉錄因子基因AhFUS3、AhABI3、AhAGL15 和AhWRI1 在ssm1 中表達均有不同程度下調。 其中,與脂肪酸合成和糖酵解關系密切的AhWRI1 基因在ssm1 中的表達量顯著下調。 因此推測,ssm1 中AhLEC1 基因的表達量顯著下調,影響AhFUS3、AhABI3 和AhWRI1 等基因的表達,進而導致ssm1 種子胚胎發育和油脂積累的異常。
關鍵詞:花生、小種子突變體ssm1、轉錄因子基因LEC1、基因表達、種子發育、油脂積累
中圖分類號:S565.2、Q786文獻標識號:A文章編號:1001-4942(2024)01-0010-07
油脂是油料作物種子主要的貯藏物質,為種子萌發和幼苗生長提供能量。 花生(Arachis hy ̄pogaea L.)是我國重要的油料作物之一,我國花生生產量和消費量均居世界第一。 我國食用油大量依賴進口,提高花生油產量是解決油料短缺的有效途徑[1] ,提高籽仁含油量和產量是花生育種的重要工作。 研究發現,多種轉錄因子對種子中油脂積累具有重要的調控作用。 例如,LEAFY COT ̄YLEDON1(LEC1)、LEC2、ABSCISIC ACID INSEN ̄SITIVE 3(ABI3)、FUSCA3(FUS3)、WRINKLED1(WRI1)和AGAMOUS ̄Like15(AGL15)等轉錄因子在調控胚胎發育、脂肪酸合成和糖酵解途徑等方面都起著重要作用[2-14] 。 LEC1 被認為是種子發育的關鍵調控因子,屬于核因子Y(nuclear fac ̄tor-Y, NF-Y)轉錄因子家族,編碼HAP3 復合體的一個亞基,通過與CCAAT 順式作用元件結合調控基因的表達[3] 。 在胚胎發育早期,LEC1 對維持胚柄細胞以及子葉的特性具有重要意義,而在胚胎發育后期,LEC1 在種子積累儲存脂肪和蛋白質等大分子物質以及種子獲得抗脫水能力過程發揮重要作用[4-13] 。 在油菜種子中過量表達Bn ̄LEC1 基因可使籽粒中的含油量提高2% ~ 20%,脂肪酸生物合成、蔗糖合成和運輸、糖酵解等途徑相關基因表達上調,這些生物學途徑所需碳通量增加[15] 。 在擬南芥幼苗、種子以及大豆種子中,利用芯片分析LEC1 基因的異位表達和缺失突變體,結果表明,ABI3、FUS 和LEC2 都受到LEC1 的調控[15] 。 WRI1 在擬南芥種子的油脂積累過程中發揮著極其重要的作用,擬南芥wri1 突變體中糖酵解和脂肪酸合成相關基因的表達呈下降趨勢[16] ,種子含油量顯著降低[17] 。 研究發現,WRI1 控制著多個糖酵解途徑關鍵酶的基因以及參與脂肪酸合成的基因的表達,包括PKP2、乙酰輔酶A 羧化酶(ACCase)、蔗糖合成酶(SucroseSynthase 2,SUS2)等,而WRI1 的表達受到LEC1和LEC2 的調控,因此,LEC 是調節脂肪酸合成代謝的關鍵基因[16,18] 。 AGL15 編碼一種MADS 結構域的轉錄因子,在種子發育過程中調節LEC2、ABI3 和FUS3 等基因的表達[19] 。
油料作物種子的正常生長發育與種子中油脂的積累密切相關。 植物種子油脂合成途徑已經較為清楚,但是調控花生種子油脂合成的分子機理仍不明確,因此闡明花生種子發育過程中油脂合成的調控機理,對花生高油育種有重要意義。 我們前期研究中獲得一個種子變小、種皮皺縮的花生突變體ssm1,與對照親本LH11 相比,該突變體種子含油量顯著降低,“球形胚”時期的胚胎明顯變小,且胚胎數目變少。 LH11 和ssm1 種子發育不同時期的轉錄組數據表明,在糖酵解途徑及脂肪酸生物合成和代謝途徑存在顯著富集,參與脂肪酸合成和代謝以及糖酵解途徑的眾多基因在ssm1 中表達顯著下調[20] 。 為分析種子早期發育中的關鍵調控基因,本研究利用qRT-PCR 技術分析了LEC1 等轉錄因子基因在突變體合子受精及種子發育時期的表達變化,以期為篩選導致突變表型的候選基因、解析花生發育和油脂積累的調控網絡、深入研究調控花生油脂合成的分子機理及培育高含油量花生新品種提供參考。
1 材料與方法
1.1 植物材料
以山東省農業科學院農作物種質資源研究所特色作物資源與分子技術創新團隊實驗室保存的魯花11 號(LH11)及其突變體ssm1 作為試驗材料。 將LH11 用60 Coγ 輻照誘變,對收獲的M2 代種子進行調查,獲得一個穩定遺傳的小種子株系,命名為ssm1(small seed mutant 1)。
2020 年4 月底將LH11 與突變體M4 代種子種植在山東省農業科學院試驗基地,株距20 cm,行距50 cm,單粒播種[21] 。 播種140 d 后收獲,對單株結果數、單株幼果數、百果重和百仁重等性狀進行統計分析。
于花生開花期和莢果發育不同階段,采集開花之前未受精的子房(DBP0)、受精1 d 的子房(DBP1)以及入土10 d(DAP10)、20 d(DAP20)、40 d(DAP40)膨大的未成熟種子。 采集的樣品用液氮速凍后放入-80 ℃冰箱保存。
1.2 種子營養成分含量測定
含油量及蛋白質、脂肪酸的含量測定分別參照NY/ T 1285—2007、GB 5009.5—2016、GB 5009.6—2016 規程進行,在中國科學院油料作物研究所完成。
1.3 RNA 提取與cDNA 合成
RNA 使用天根公司的RNA 提取試劑盒進行提取,然后利用艾瑞克生物公司的反轉錄試劑盒進行反轉錄合成cDNA。
1.4 qRT-PCR 檢測
參照陳娜等[22] 的方法和條件進行qRT -PCR。 qRT-PCR 的SYBR Green PreMix 試劑盒購自艾瑞克生物公司。 qRT-PCR 反應體系[23] 為20μL,包括2×SYBR? Green Pro Taq HS Premix 10μL,稀釋至8 ng/ μL 的cDNA 2 μL,正、反向引物各1 μL,ddH2O 6 μL。 反應程序:95 ℃ 30 s、95 ℃ 5 s,60 ℃ 30 s,40 個循環。 以花生Actin 基因為內參,每個時期的樣品重復3 次,采用2-ΔΔCt法[24] 計算目標基因的相對表達量,用log2FC[log2(Fold Change)] 值表示基因相對表達的差異倍數。 qRT-PCR 所用引物見表1。
2 結果與分析
2.1 花生ssm1 突變體表型分析
籽仁大小是影響花生產量的重要農藝性狀之一[25] 。 突變體ssm1 的莢果和種子明顯小于LH11,且種皮皺縮、顏色較深(圖1)。 以LH11 作為對照,對成熟后的突變體ssm1 的部分農藝性狀進行分析,發現ssm1 的單株幼果數增多,單株結果數顯著增加,分別較對照增加38.46%和88.04%、百果重、百仁重分別較對照降低64.96%、60.31%,差異達極顯著水平、公斤仁數極顯著增加,約是LH11 的2.5 倍、雖然突變體ssm1 的結果數多,但由于種子變小,幼果數多,產量大幅度降低,單株生產力較LH11 降低34.11%(表2)。
2.2 種子營養成分含量與脂肪酸的組成
對ssm1 和LH11 干種子粗蛋白含量、含油量和主要脂肪酸含量進行測定,結果(表3)表明,ssm1 的含油量顯著降低,比LH11 降低13.53%。ssm1 突變體的粗蛋白、棕櫚酸、花生酸含量稍有升高,與LH11 差異較小、硬脂酸、油酸和亞油酸含量略降低,十七碳烷酸、十七碳一烯酸分別比LH11 降低33.33%、25.00%,而亞麻酸、花生一烯酸、山崳酸、芥酸、二十四碳烷酸含量明顯升高。
2.3 LEC1 等轉錄因子基因的表達
利用qRT-PCR 分別檢測了LH11 和突變體ssm1 中AhLEC1A、AhLEC1B、AhFUS3、AhABI3、AhAGL15、AhWRI1 基因在花生合子受精及種子發育不同階段的表達情況,結果(圖2) 表明,在LH11 中,AhLEC1A、AhLEC1B、AhFUS3、AhWRI1 和AhAGL15 在種子發育早期的表達總體上均呈現先上升后下降的趨勢,而AhABI3 的表達呈一直上升趨勢。 在ssm1 中, 各基因表達水平明顯低于LH11,其中,AhLEC1A 和AhABI3 表達呈持續上調趨勢、AhWRI1 在DBP0—DAP10 時期表達量升高,之后變化不大,各時期表達量均顯著低于LH11、其余基因表達變化趨勢基本與LH11 中一致,且除AhLEC1B 在DBP0 和DAP40 以及Ah ̄FUS3 在DAP40 時期外, 表達量均顯著低于LH11。 AhLEC1A 和AhLEC1B 在LH11 和ssm1 的未受精子房( DBP0) 中差異不明顯, 但是在DBP1—DAP20 時期,即受精后子房以及胚胎發育早期,在ssm1 中的表達顯著下調,且AhLEC1B 表達量的升高具有明顯的延遲性。 其中,AhLEC1A在DAP10、DAP20 時期表達量分別比LH11 降低52.04%、45.44%、AhLEC1B 在種子胚胎發育早期(DAP10)表達較LH11 大幅下調,降低77.88%。而LEC1 下游調控的轉錄因子基因AhFUS3、AhA ̄GL15 在種子發育早期(DBP0—DAP20 時期) 的表達量逐步升高,但ssm1 的表達量顯著低于LH11、AhABI3 在ssm1 中的表達量明顯低于LH11,且除DBP1 外, 均達顯著水平。 推測在ssm1 中,AhLEC1 可能作為上游關鍵調控因子,其表達下調引起AhFUS3、AhABI3、AhWRI1 等下游轉錄因子基因表達量降低。 例如,AhWRI1 在ssm1種子發育過程中表達量顯著降低,在營養物質積累的DAP20 時期其在ssm1 中的表達量比LH11下調82.94%。
3 討論與結論
LEC1 與植物胚胎發育和種子成熟密切相關,作為先驅轉錄因子在胚胎發育早期就被激活表達,原胚期和整個胚胎發生過程中都能檢測到其表達[26] ,是調控胚胎成熟過程的關鍵因子[27] 。前期對LH11 和ssm1 的DAP10、DAP20、DAP40時期種子進行了轉錄組測序,發現與脂肪酸合成和代謝以及糖酵解相關通路中的眾多基因在ssm1 中表達顯著下調。 本研究對LEC1 等轉錄因子基因的表達模式進行分析,發現AhLEC1A 和AhLEC1B 在ssm1 合子受精和種子發育早期的表達與LH11 相比大幅度降低且表達延遲。 另外,在子房受精之前,AhLEC1A、AhLEC1B 的表達量在LH11 和ssm1 中差異并不明顯,但是受精后開始出現顯著差異,直至DAP40 時期兩者之間又無明顯差異,說明LEC1 在合子受精后開始對早期胚胎發育進行調控。
在LEC1 過表達的油菜和擬南芥中發現脂肪酸種類和脂質水平顯著增加[15] 、甘藍型油菜lec1突變體中種子含油量下降,與糖酵解和脂肪酸修飾的相關基因表達下調[28] 。 WRI1 控制著多個糖酵解途徑關鍵酶的基因以及參與脂肪酸合成的基因的表達。 擬南芥wri1 突變體中關于糖酵解和脂肪酸代謝相關酶活性降低,wri1 種子的含油量降低80%[29] 。 這些轉錄因子基因在胚胎和種子發育時期特異性表達,過量表達或抑制表達這些基因能提高含油量,證明了它們在調控種子脂肪合成方面的重要作用[15,28] 。 在本研究中,AhL ̄EC1A、AhLEC1B 正常表達的LH11 中其下游轉錄因子基因AhFUS3 和AhWRI1 在DAP10 時期后大幅度升高、但是在ssm1 中,AhLEC1A、AhLEC1B 在DBP1—DAP20 時期表達下調且顯著低于LH11,AhFUS3 和AhWRI1 的表達量同樣顯著下調且表達量的升高出現延遲性,同時,受到LEC1 調控的ABI3 和調節FUS3 等基因表達的AGL15 轉錄因子基因在ssm1 中的表達量也降低。 DAP10—DAP20 時期是養分物質積累的重要時期, 在LH11 中LEC1 以及下游轉錄因子基因的表達趨勢與油脂合成相關基因表達趨勢一致,都呈上升趨勢,與油脂積累同步、而在ssm1 中,LEC1 及其下游轉錄因子基因表達與LH11 相比在同時期顯著下調,油脂合成相關基因表達也顯著下調。 作為LEC1 靶基因的WRI1 轉錄因子在花生小種子突變體ssm1 中表達顯著降低,而受WRI1 轉錄因子調控的PKp2B(丙酮酸激酶)在ssm1 種子中表達量幾乎為零,糖酵解途徑中FBPB(果糖-1,6-二磷酸酶)和PFKB(磷酸果糖激酶)基因在ssm1中表達也顯著下調[20] ,影響了脂肪酸從頭合成所需原料乙酰CoA 的正常合成。
ssm1 表現出種子皺縮、種皮顏色較深以及含油量下降的表型,這與擬南芥wri1 突變體[29] 和甘藍型油菜lec1 突變體[28] 表型基本一致。 可以推測,AhLEC1A、AhLEC1B 作為上游轉錄因子,其表達量的顯著下調影響了AhFUS3、AhABI3 和Ah ̄WRI1 等下游轉錄因子的表達,進而影響了脂肪酸合成和代謝以及糖酵解途徑中基因的表達,導致種子含油量明顯降低。 因此推測LEC1 基因參與了早期胚胎的發育調控,并且激活了下游轉錄因子的表達,共同參與脂肪合成和糖酵解途徑。但是關于LEC1 及其下游轉錄因子如何調控脂肪合成通路還有待進一步研究。
參 考 文 獻:
[1] 黃莉,陳玉寧,羅懷勇,等. 花生種子大小相關性狀QTL 定位研究進展[J]. 作物學報,2022,48(2):280-291.
[2]Angeles ̄Nú?ez J G, Tiessen A. Mutation of the transcriptionfactor LEAFY COTYLEDON2 alters the chemical compositionof Arabidopsis seeds, decreasing oil and protein content, whilemaintaining high levels of starch and sucrose in mature seeds[J]. Journal of Plant Physiology, 2011,168(16):1891-1900.
[3] Lee H, Fischer R L, Goldberg R B, et al. Arabidopsis LEAFYCOTYLEDON1 represents a functionally specialized subunit ofthe CCAAT binding transcription factor[J]. Proceedings of theNational Academy of Sciences of the United States of America,2003,100(4):2152-2156.
[4] Pelletier J M, Kwong R W, Park S, et al. LEC1 sequentiallyregulates the transcription of genes involved in diverse develop ̄mental processes during seed development[J]. Proceedings ofthe National Academy of Sciences of the United States of Amer ̄ica, 2017,114(32):6710-6719.
[5] Zhang S B, Wong L, Meng L, et al. Similarity of expressionpatterns of knotted1 and ZmLEC1 during somatic and zygoticembryogenesis in maize (Zea mays L.)[J]. Planta, 2002,215(2):191-194.
[6] Yazawa K, Takahata K, Kamada H. Isolation of the gene enco ̄ding carrot LEAFY COTYLEDON1 and expression analysisduring somatic and zygotic embryogenesis[J]. Plant Physiologyand Biochemistry, 2004,42(3):215-223.
[7] Fambrini M, Durante C, Cionini G, et al. Characterization ofLEAFY COTYLEDON1 ̄LIKE gene in Helianthus annuus and itsrelationship with zygotic and somatic embryogenesis[J]. Devel ̄opment Genes and Evolution, 2006,216(5):253-264.
[8] Alemanno L, Devic M, Niemenak N, et al. Characterization ofleafy cotyledon1 ̄like during embryogenesis in Theobroma cacaoL.[J]. Planta, 2008,227(4):853-866.
[9] Cao S, Kumimoto R W, Siriwardana C L, et al. Identificationand characterization of NF ̄Y transcription factor families in themonocot model plant Brachypodium distachyon [ J ]. PLoSONE, 2011,6:107-126.
[10] Salvini M, Sani E, Fambrini M, et al. Molecular analysis of asunflower gene encoding an homologous of the B subunit of aCAAT binding factor[J]. Molecular Biology Reports, 2012,39(6):6449-6465.
[11] West M A L, Yee K M, Danao J, et al. LEAFY COTYLE ̄DON1 is an essential regulator of late embryogenesis and cotyle ̄don identity in Arabidopsis[J]. The Plant Cell, 1994,6(12):1731-1745.
[12] Mu J Y, Tan H L, Zheng Q, et al. LEAFY COTYLEDON1 is akey regulator of fatty acid biosynthesis in Arabidopsis[J]. PlantPhysiology, 2008,148(2):1042-1054.
[13] Junker A, M?nke G, Rutten T, et al. Elongation ̄related func ̄tions of LEAFY COTYLEDON1 during the development of Ara ̄bidopsis thaliana[J]. The Plant Journal: for Cell and Molecu ̄lar Biology, 2012,71(3):427-442.
[14] Tan H L, Yang X H, Zhang F X, et al. Enhanced seed oil pro ̄duction in canola by conditional expression of Brassica napusLEAFY COTYLEDON1 and LEC1 ̄LIKE in developing seeds[J]. Plant Physiology, 2011,156(3):1577-1588.
[15] Kagaya Y, Toyoshima R, Okuda R, et al. LEAFY COTYLE ̄DON1 controls seed storage protein genes through its regulationof FUSCA3 and ABSCISIC ACID INSENSITIVE3[J]. Plant &Cell Physiology, 2005,46(3):399-406.
[16] Ruuska S A, Girke T, Benning C, et al. Contrapuntal networksof gene expression during Arabidopsis seed filling[J]. The PlantCell, 2002,14(6):1191-1206.
[17] Cernac A, Benning C. WRINKLED1 encodes an AP2/ EREBdomain protein involved in the control of storage compound bio ̄synthesis in Arabidopsis[J]. The Plant Journal: for Cell andMolecular Biology, 2004,40(4):575-585.
[18] Baud S,Mendoza M S,To A,et al. WRINKLED1 specifies theregulatory action of LEAFY COTYLEDON2 towards fatty acidmetabolism during seed maturation in Arabidopsis [ J]. ThePlant Journal: for Cell and Molecular Biology, 2007,50(5):825-838.
[19] Zheng Y M, Ren N, Wang H, et al. Global identification oftargets of the Arabidopsis MADS domain protein AGAMOUS ̄Like15[J]. Plant Cell, 2009,21: 2563-2577.
[20] Guo F, Zhu X, Zhao C, et al. Transcriptome analysis and geneexpression profiling of the peanut small seed mutant identifiedgenes involved in seed size control[J]. International Journal ofMolecular Sciences, 2022,23(17):9726-9731.
[21] 李歡倪,仇靜靜,馬俊杰,等. 花生半矮化突變體sdm1 的表型分析與赤霉素響應研究[J]. 山東農業科學,2017,49(12):1-5.
[22] 陳娜,潘麗娟,遲曉元,等. 花生果糖-1,6-二磷酸醛縮酶基因AhFBA1 的克隆與表達[J]. 作物學報,2014,40(5):934-941.
[23] 潘麗娟,梁丹,劉風珍,等. 花生AhFUSCA3 基因的原核表達及在非生物脅迫下的表達分析[J]. 花生學報,2016,45(4):1-4,29.
[24] Livak K J, Schmittgen T D. Analysis of relative gene expressiondata using real ̄time quantitative PCR and the 2-ΔΔCt method[J]. Methods,2001,25(4):402-408.
[25] 郭建斌,李威濤,羅懷勇,等. 花生種子大小相關性狀QTL定位及與出仁率的關系[J]. 植物遺傳資源學報,2022,23(5):1465-1473.
[26] Tamar L, Masaaki O, Kelly M Y, et al. Arabidopsis LEAFYCOTYLEDON1 is sufficient to induce embryo development invegetative cells[J]. Cell,1998,93(7):5-18.
[27] 郭曉芳,嚴海燕,何驥,等. 花生ABI3 同源基因的定位分析[J]. 武漢植物學研究,2006,24(1):22-26.
[28] Thomas K, Gil S, Christiane V, et al. Regulation of storageprotein gene expression in Arabidopsis [ J ]. Development,2003,130(24):6065-6073.
[29] Focks N, Benning C. wrinkled1: a novel, low ̄seed ̄oil mutantof Arabidopsis with a deficiency in the seed ̄specific regulationof carbohydrate metabolism[J]. Plant Physiology, 1998,118(1):91-101.