花生油脂合成相关基因的表达谱分析

许 静 潘丽娟 李昊远 王 通 陈 娜 陈明娜 王 冕 禹山林 侯艳华,* 迟晓元,*

花生油脂合成相关基因的表达谱分析

许 静1,**潘丽娟1,**李昊远2王 通1陈 娜1陈明娜1王 冕1禹山林1侯艳华2,*迟晓元1,*

1山东省花生研究所, 山东青岛 266100;2哈尔滨工业大学(威海)海洋科学与技术学院, 山东威海 264209

为研究不同发育时期花生籽仁油脂合成过程中基因表达调控模式, 本研究以高油酸、中油花生品系F18和低油酸、低油花生品种‘鲁花6号’为材料, 对下针后10、30、40、60 DAP (days after pegging)的花生种子进行表达谱芯片测序。结果表明, 130、3556、2783个基因分别在30、40、60 DAP时期差异表达。GO注释和KEGG富集结果显示, 差异表达的基因主要富集在脂肪酸合成和光合等代谢进程中, 其中、、等主要参与油酸的积累, 参与光合作用的基因均为捕光叶绿素/结合蛋白, 全部上调表达。代谢通路图结果表明, 籽仁发育的40 DAP和60 DAP时期, 脂肪酸合成途径的基因均上调表达。研究结果为花生油脂代谢的分子机制提供理论基础, 同时也为花生品质改良贡献了基因资源。

花生; 基因芯片; 差异表达分析; 油脂合成相关基因

花生(L.)是世界上重要的油料作物和经济作物, 也是重要的食用油和食用蛋白源, 具有较高的经济价值[1]。我国是世界上最大的花生消费国, 每年花生总产量的46%~48%用于榨油。花生籽仁中脂肪含量占50%左右, 其中80%为不饱和脂肪酸, 能够降低有害胆固醇, 减缓动脉粥样硬化, 有效预防冠心病等心脑血管疾病[2], 是人们理想的健康食用油。因此, 提高花生含油量是国家发展的重大需求, 市场潜力巨大, 发展前景广阔。

植物油脂的合成过程复杂, 涉及多种酶的催化作用。甘油-3-磷酸酰基转移酶(-glycerol-3-phosphate acyltransferase, GPAT)、溶血磷脂酸酰基转移酶(lysophosphatidyl acyltransferase, LPAAT)和二酰甘油酰基转移酶(diacylglycerol acyltransferase, DGAT)是参与植物油脂合成的3个关键酶。近年来, 花生油脂合成方面也取得了一定的研究进展。陈玉梅等[3]对花后20、35、50、60 d的花生籽粒进行转录组测序, 共得到135条代谢通路, 其中与油脂合成相关的有脂肪酸生物合成、不饱和脂肪酸生物合成、脂肪酸延长、脂肪酸代谢、花生四烯酸代谢等代谢通路。Chi等[4]在花生中分离得到6个基因, RT-PCR分析发现, 在花生籽仁发育时期的转录水平逐渐增加, 同时在逆境环境下,基因在叶片或根中的表达显著增加。3个基因从花生中被克隆出来,在籽仁发育的初期表达较低, 但在果针入土50 d时表达增加,在果针入土10 d和50 d高表达,在果针入土20 d和50 d表达量最高[5]。Zheng等[6]将转到烟草中, 与野生型相比, 脂肪酸含量增加14.7%~20.9%。Chi等[7]通过RT-PCR或RACE-PCR在花生中克隆得到4个基因,基因在籽粒中富集表达, 主要集中在籽粒发育初期, 然而基因在籽粒发育时期表达量下降, 主要在花中富集表达。Chen等[8]克隆了花生中的基因, RT-PCR表明在花后40 d和50 d表达量与种子中油脂的积累速率同步增加。皮广静等[9]分析高油品系农大D666及其亲本不同发育时期籽仁内三类酰基转移酶基因的表达量表明, R5期和、基因表达量之间协调性较好, 有利于油脂的合成与积累。

硬脂酰-ACP去饱和酶(Stearoyl-ACP desaturase, FAB2)是实现硬脂酸转化为不饱和脂肪酸的关键酶, 微粒体Δ12去饱和酶(Microsomal Δ12 desaturase, FAD2)是调控花生油酸、亚油酸含量和油亚比(O/L)的关键酶。在酵母中异源表达和发现了亚油酸的形成, 证实了两基因的功能[10]。此外, Chi等[11]克隆了一系列的脂肪酸去饱和酶, RT-PCR分析发现和在籽粒中较其他组织表达量高, 在籽粒发育时期,、、表达量逐渐上升达到最高水平, 然后降低。转录因子WRINKLED(WRI1)是在种子皱缩的拟南芥突变体中发现并命名, 研究发现该基因会导致突变体种子含油量降低80%[12]。随后, WRI1转录因子常被用于油料作物和模式植物油分改良相关的研究。Pouvreau等[13]研究表明, 过表达基因, 能使转基因玉米种子中含油量增加, 同时还发现基因具有调节细胞质中三酰甘油主链生成的作用。孙金波等[14]从二倍体花生野生种和四倍体栽培种克隆到20个WRI1转录因子, 转录组分析表明在种子中的表达高于其他组织, RT-PCR结果进一步表明,表达随种子成熟水平逐渐增加。为进一步研究花生油脂合成的调控机制, 本研究选用油酸和油脂含量差异显著的2个花生品种(系)。通过基因芯片测序技术, 研究花生种子不同发育阶段的全基因组表达谱, 并分离克隆出与种子发育有关, 特别是与油脂或脂肪酸合成和调控有关的特异表达基因, 以期为花生品质的遗传改良奠定基础。

1 材料与方法

1.1 材料

花生(L.)材料选择山东省花生研究所杂交育种课题组培育的高油酸、中油花生品系F18和低油酸、低油花生品种‘鲁花6号’。

1.2 方法

1.2.1 样品采集 分别在10、30、40、60 DAP选取无病害、长势一致的花生荚果, 剥去花生壳后迅速将籽仁放入液氮中速冻,-80℃保存, 作为芯片测序的样品。另取各个时期花生荚果适量, 作为油脂和脂肪酸含量测定的样品。

1.2.2 花生籽仁油脂和脂肪酸含量测定 采用索氏提取法[15]提取油脂, 脂肪酸的提取及气相色谱分析方法具体参考迟晓元等[16]的方法。样品采用3 g左右研磨成粉末的花生种子, 每个样品重复3次。

1.2.3 RNA提取和纯化 采用北京天根生化科技有限公司试剂盒(RNeasy Mini Kit)提取样品总RNA。使用NucleoSpin RNA纯化试剂盒(MAC HEREYNAGEL, 德国) 进行RNA样品的纯化。用甲醛琼脂糖凝胶电泳对RNA质检, 然后用Nanovue Plus分光光度计(General Electric Company, USA)对RNA进行定量。所有RNA样品均保存于-80℃冰箱备用。

1.2.4 基因芯片制备 芯片测序工作由厦门博奥生物科技有限公司完成。采用Roche公司的NimbleGen花生表达谱芯片, 芯片包含36,158个探针组。Roche NimbleGen表达谱芯片试验流程按照博奥生物公司的规定流程进行, 具体方法参考Yu等[17]的方法。

1.2.5 数据分析 利用NimbleScan v2.6软件将杂交芯片的图像信号转换为数字信号。采用Genespring v13.0软件分位数法进行数据标准化处理, 以上工作均由博奥生物有限公司完成。根据统计学分析, 差异基因的筛选条件为: 校正后的-value≤0.05且2个样品的差异倍数|fold change|≥2。

1.2.6 GO功能注释及富集分析 首先利用blast+ v2.2.30将基因序列与NCBI中的NR、Swissprot、TrEm1等数据库进行比对。然后将比对结果导入Blast2GO中用于GO注释, 主要包括分子功能、生物过程和细胞组分信息注释。最后利用WEGO对GO功能进行分类统计(http://wego.genomics.org.cn/cgi-bin/wego/index.pl)。

1.2.7 KEGG Pathway分类与富集分析 利用KEGG网站(http://www.kegg.jp/)对获得的序列进行代谢通路分类, 同一个基因序列可能注释到多个代谢通路中。使用R软件中ggplot2软件包进行富集分析。

1.2.8 蛋白互作网络分析 将F18和‘鲁花6号’40 DAP和60 DAP 2个时期的所有差异表达的基因序列比对结果上传至STRING数据库(https://string-db.org/)中, 分析蛋白之间的互作关系。然后导入到cytoscape中进行编辑作图。

1.2.9 差异表达基因的RT-PCR验证 参照陈娜等[18]的方法进行荧光定量PCR的反应。每个样品3次重复, 内参基因选择, 采用2–DDCt方法分析数据。荧光定量PCR所用引物见表1。

2 结果与分析

2.1 不同发育时期花生籽仁油脂和脂肪酸含量分析

F18和‘鲁花6号’不同发育时期花生籽仁含油量和脂肪酸组分含量如表2, 两品种(系)含油量和油酸含量随着籽仁的发育不断积累。除了30 DAP发育时期, F18在各个发育时期的含油量均高于‘鲁花6号’, 两品种(系)在60 DAP发育时期含油量均最高, 分别为50.81%和45.44%。F18在4个发育时期的油酸含量均高于‘鲁花6号’, 在40 DAP发育时期达76.06%, 而‘鲁花6号’在60 DAP达到最高, 为43.98%。亚油酸含量在两品种(系)中随着籽仁发育呈下降趋势, F18和‘鲁花6号’均在60 DAP发育时期达到最低, 为1.55%和34.28%, 油亚比均最高为45.65和1.28。

表1 RT-PCR验证基因所用的引物

2.2 不同发育时期花生籽仁差异表达基因分析

根据芯片杂交结果中不同发育时期籽仁基因表达水平, 筛选出显著差异表达的基因。10 DAP时期两品种(系)基因表达水平差异不显著(图1-A); 30 DAP时期在两品种(系)中差异表达的基因共130个, 与‘鲁花6号’相比, F18上调表达基因20个, 下调表达基因110个(图1-B); 两品种(系)在40 DAP时期差异表达的基因数目最多, 达3556个, 其中上调的差异表达基因有2056个, 下调的有1500个(图1-C); 60 DAP时期两品种(系)差异表达的基因共2783个, 相对于‘鲁花6号’, F18上调表达的基因有1746个, 下调表达的基因有1037个(图1-D)。

表2 不同发育时期两品种(系)油脂和脂肪酸含量

A: 10 DAP差异表达基因; B: 30 DAP差异表达基因; C: 40 DAP差异表达基因; D: 60 DAP差异表达基因。红色表示F18相对于‘鲁花6号’差异表达的基因。

A: different expression genes of 10 DAP; B: different expression genes of 30 DAP; C: different expression genes of 40 DAP; D: different expression genes of 60 DAP. Red blots represent the different expression genes of F18 compared with ‘Luhua 6’.

根据筛选的差异表达基因分析发现, 一些基因在花生籽仁发育的不同时期均有表达。由图2可知, 与‘鲁花6号’相比, F18在30、40、60 DAP 3个时期同时表达上调的基因有6个, 下调表达的基因为50个; 在40、60 DAP 2个时期共同差异表达上调的有712个, 差异表达下调的有362个。

A: F18相比于‘鲁花6号’上调表达基因; B: F18相比于‘鲁花6号’下调表达基因。

A: the up-regulated expression genes of F18 compared with ‘Luhua 6’; B: the down-regulated expression genes of F18 compared with ‘Luhua 6’.

2.3 不同发育时期差异表达基因的GO功能注释分析

对筛选到的差异表达基因进行GO功能注释, 3个不同发育时期分别注释22、1914和1559个差异基因, 主要包括生物过程、细胞组分和分子功能三大类别。如图3-A所示, 40 DAP和60 DAP时期差异表达基因主要参与代谢调控(metabolic process GO: 0008152)和催化进程(catalytic activity GO: 0003824)。有趣的是, 参与光合作用的基因全部上调(53个)而没有下调, 这些参与光合作用的基因大量上调导致F18固定空气中更多的CO2, 并合成为自身碳源, 最终流向脂质合成。此外, 还发现前体代谢物和能量的生成上调的基因有42个, 而下调仅有1个。这些光合作用、糖、蛋白等代谢也对油脂合成和积累有一定影响。进一步对40 DAP和60 DAP 2个时期差异表达的1074个基因进行富集发现, 共有598个基因富集到11个GO类别中(图3-B)。其中代谢进程富集的基因最多, 这表明‘鲁花6号’和F18参与代谢途径的基因的表达有显著差异。

2.4 不同发育时期差异表达基因的KEGG Pathway富集分析

根据KEGG数据库, 本研究对花生籽仁40 DAP和60 DAP 2个发育时期1074个差异表达的基因(上调表达基因712个, 下调表达基因362个)进行KEGG Pathway富集分析。上调表达的基因主要富集在光合、氰氨基酸代谢、脂肪酸合成和淀粉蔗糖代谢4个代谢途径中(图4-A), 注释的基因分别为25、12、21和9个。下调表达的基因只富集在柠檬烯和蒎烯降解途径中, 没有富集到与脂肪酸代谢相关的途径(图4-B)。这表明, 与‘鲁花6号’相比, F18中参与脂肪酸合成的基因在花生籽仁发育的后期表达量普遍较高。另外, 本研究发现, 光合途径差异表达基因编码的都是捕光叶绿素/结合蛋白, 这些蛋白复合体能够捕获并传递光能, 调节光能分配, 而光合作用是植物一切物质积累的源头, 这一发现反应了2个品种(系)自身的遗传差异。

针对我们感兴趣的脂肪酸合成途径, 绘制了代谢通路图(图5)。结果显示, 在F18和‘鲁花6号’后2个发育时期脂肪酸合成途径中没有找到下调的基因, 说明在这一时期中, F18的脂肪酸合成代谢比‘鲁花6号’更旺盛, 产生更多的脂肪酸。在代谢通路的最开始, 第1个上调的是ACCase (EC: 6.4.1.2)基因, 催化脂肪酸生物合成的起始; 而后在脂肪酸合成循环的通路中, β-酮脂酰ACP合酶(FabF)、β-酮脂酰ACP还原酶(FabG)、β-羟酰基ACP脱水酶(FabZ)和烯酰-ACP还原酶(FabI)均上调, 他们共同催化脂肪酸的缩合反应; 在途径的末端, 长链酰基合成酶(EC:6.2.1.3)基因上调表达, 合成长链的脂肪酸。

2.5 油脂代谢相关基因的筛选

根据GO注释的结果, 我们分别在40 DAP和60 DAP籽仁发育时期筛选到64个和62个参与油脂代谢相关的基因, 其中分别鉴定到25个和27个油脂合成的关键基因(表3和表4)。硬脂酰ACP去饱和酶(FAB2)是催化硬脂酸产生油酸的关键酶, ω6脂肪酸去饱和酶(FAD2)是催化油酸C12位去饱和酶形成亚油酸的一类酶。花生籽仁发育的40 DAP时期, F18花生品系中两基因均上调表达(和), 比‘鲁花6号’花生品种中的表达量分别高3.48倍和2.06倍。在60 DAP发育时期,起关键作用, 与‘鲁花6号’相比, F18中()表达量显著增加, 高达29.89倍。

2.6 WRI1蛋白互作网络分析

在花生籽仁发育后期的差异表达基因中筛选到WRI1转录因子, 以WRI1为核心进行蛋白互作网络分析, 共发现9个与WRI1互作的蛋白(图6)。其中KASI (3-酮酰基酰基载体蛋白合酶)、MOD1 (烯酰基载体蛋白还原酶)和GLB1 (氮调控蛋白)与脂质代谢有关, SUC5 (蔗糖转运蛋白)和UGT84B1 (葡萄糖基转移酶)参与糖代谢, 进一步使碳源从糖代谢流向脂肪酸合成, 促进脂肪酸合成。

2.7 不同发育时期差异表达基因的RT-PCR验证

为验证芯片测序结果中基因表达量的准确性, 随机挑选了不同发育时期的16个基因进行RT-PCR验证, 16个差异表达基因及内参基因引物如表1所示。将RT-PCR结果与芯片测序数据比对发现, 16个差异表达基因与芯片数据基本吻合(图7和表5), 表明芯片测序数据的可靠性, 可以用来分析不同发育时期花生籽仁差异表达的基因。

红色代表相对于‘鲁花6号’, F18上调表达的基因。

Red represents the up-regulated expression genes of F18 compared with ‘Luhua 6’.

表3 40 DAP发育时期油脂代谢相关基因

(续表3)

(续表3)

a表示F18相对于‘鲁花6号’的上调和下调。

arepresents the up- and down-regulated expression genes of F18 compared with ‘Luhua 6’.

表4 60 DAP发育时期油脂代谢相关基因

(续表4)

a表示F18相对于‘鲁花6号’的上调和下调。

arepresents the up- and down-regulated expression genes of F18 compared with ‘Luhua 6’.

3 讨论

油酸是有益脂肪酸, 能够减少心脑血管的发生, 维持身体健康[19]。最近几年已有多家单位开展高油酸花生育种工作: 姜慧芳等[20]发现, 远缘杂交后脂肪酸的变异远超过亲本间的差异, 并且获得了6份油酸含量达64.0%以上且棕榈酸含量在8.5%以下的新种质; 于明洋等[21]将‘花育22’与‘开农1716’杂交并不断回交, 获得了24个油酸含量高于70%的花生新品系。然而, 高油酸代谢的分子机制研究相对较少, 尤其是不同花生品种(系)之间。本研究对‘鲁花6号’和F18不同发育时期的花生籽仁进行芯片测序, 分别在10 DAP、30、40、60发育时期筛选到130、3556、2783个差异表达的基因(10 DAP发育时期无差异表达基因)。本研究发现在发育初期(10 DAP和30 DAP)两品种(系)之间差异表达基因较少, 而在40 DAP发育时期差异表达基因数目最多, 60 DAP发育时期相对减少, 这与植物种子油脂积累呈现慢—快—慢的变化规律相一致[22]。

黑色和灰色分别代表相对于‘鲁花6号’, F18上调和下调表达的基因。

Black and grey represent the up- and down-regulated expression genes of F18 compared with ‘Luhua 6’.

表5 差异表达基因的RT-PCR验证结果

(续表5)

a表示F18相对于‘鲁花6号’的上调和下调。

arepresents the up- and down-regulated expression genes of F18 compared with ‘Luhua 6’.

GO注释及富集分析发现, 40 DAP和60 DAP 2个时期共有416个差异表达的基因富集到代谢进程(metabolic process)。其中相比于‘鲁花6号’, F18在脂质代谢进程中上调表达基因居多, 表明F18在花生籽仁发育后期, 脂质代谢相关基因活跃, 油脂合成加快, 这可能是造成两品种(系)间含油量差异的原因。GO富集分析发现, F18中参与光合作用的基因表达量均比‘鲁花6号’高, 这将促使F18固定更多的碳源并流向油脂合成途径。此外, F18中前体代谢和能量生成途径上的基因上调表达。这些光合作用、糖、蛋白等代谢进程的活跃将促进油脂合成和积累, 与前人研究相符[23]。

已有研究证明和在改善油料作物脂肪酸组分方面发挥了作用。与野生型相比, 转基因的烟草植株不饱和脂肪酸(尤其是油酸)含量更高[24]。FAD2是调控花生油酸、亚油酸含量和油亚比(O/L)的关键酶,酶活性的丧失是产生高O/L比值花生的主要原因之一[25]。本研究发现, 60 DAP花生发育时期, F18中两基因均上调表达, 但的表达量差异倍数(29.89)远高于(3.82), 这表明F18品系花生籽仁发育后期,基因起主要作用, 因此导致F18的油酸含量高于‘鲁花6号’。这与官梅等[26]利用芯片测序的结果一致。

转录因子与基因启动子区相结合调控基因的表达, 有研究发现, 一些转录因子能够参与植物油脂代谢[27]。WRI1 (WRINKLED1)转录因子是AP2/EREBP家族的一员, Maeo等[28]发现, 超表达基因能够上调脂肪酸合成途径的一系列基因, 包括、、、等。Fukuda等[29]的研究表明,能够与脂肪酸合成相关基因的上游序列AW-box结合, 调控乙酰辅酶A羧化酶基因的表达, 进而调控油脂代谢。本研究在花生籽仁发育的后期筛选到了转录因子, 其在F18中的表达量比在‘鲁花6号’高3.99倍, 但是本文中筛选到的与油脂合成相关的差异表达基因是否受到该转录因子的调控有待进一步研究。

以上数据表明, 花生籽仁油脂的合成是一个复杂的过程, 受许多基因和转录因子的共同调控。研究结果可为花生籽仁油脂合成机理研究提供理论基础, 为花生品质育种提供新的基因资源, 但是相关基因复杂的调控网络及功能有待进一步验证。

4 结论

分别在30、40、60 DAP时期检测到130、3556、2783个差异表达的基因。GO注释、KEGG富集和代谢通路图表明, 差异表达的基因主要参与脂肪酸合成和光合等代谢进程, 其中硬脂酰-ACP去饱和酶(FAB2)、微粒体Δ12去饱和酶(FAD2)、WRI1转录因子主要参与油酸的积累。

[1] Sarvamangala C, Gowda M, Varshney R. Identification of quantitative trait loci for protein content, oil content and oil quality for groundnut (L.)., 2011, 122: 49–59.

[2] Huth P J, Fulgoni III V L, Larson B T. A systematic review of high-oleic vegetable oil substitutions for other fats and oils on cardiovascular disease risk factors: implications for novel high-oleic soybean oils., 2015, 6: 674–693.

[3] 陈玉梅, 李璐璐, 陈锦玲, 徐媛, 李惠敏, 秦新民. 基于转录组测序的花生籽粒不同发育时期油脂合成相关基因差异表达分析. 河南农业科学, 2019, 48(7): 24–37. Chen Y M, Li L L, Chen J L, Xu Y, Li H M, Qin X M. Differential expression analysis of genes related to lipid synthesis through transcriptome sequencing during different developmental stages in peanut seed., 2019, 48(7): 24–37 (in Chinese with English abstract).

[4] Chi X Y, Yang Q L, Pan L J, Chen N, Wang T, Wang M, Yang Z, Guan X, Yu S L. Isolation and expression analysis of glycerol-3-phosphate acyltransferase genes from peanuts (L.)., 2015, 66: e093.

[5] Chi X Y, Hu R B, Zhang X W, Chen M N, Chen N, Pan L J, Wang T, Wang M, Yang Z, Wang Q F, Yu S L. Cloning and functional analysis of three diacylglycerol acyltransferase genes from peanut (L.)., 2014, 9: e105834.

[6] Zheng L, Shockey J, Guo F, Shi L M, Li X G, Shan L, Wan S B, Peng Z Y. Discovery of a new mechanism for regulation of plant triacylglycerol metabolism: the peanutgene family transcriptome is highly enriched in alternative splicing variants., 2017, 219: 62–70.

[7] Chi X Y, Dong F, Yang Q L, Chen M N, Chen N, Pan L J, Wang T, Wang M, Yang Z, He Y N, Yu S L. Expression and characterization ofgenes from peanut (L.)., 2014, 15: 141–153.

[8] Chen S L, Huang J Q, Lei Y, Zhang Y T, Ren X P, Chen Y N, Jiang H F, Yan L Y, Li Y R, Liao B S. Identification and characterization of a gene encoding a putative lysophosphatidyl acyltransferase from., 2012, 37: 1029–1039.

[9] 皮广静, 刘风珍, 万勇善, 张昆, 吕玉英, 张秀荣. 花生高油品系农大D666及其亲本油脂合成酰基转移酶基因的表达分析. 分子植物育种, 2018, 16: 1057–1065. Pi G J, Liu F Z, Wan Y S, Zhang K, Lyu Y Y, Zhang X R. Expression analysis of Acyltransferase genes involved in oil biosynthesis in high-oil peanut line Nongda D666 and parents., 2018, 16: 1057–1065 (in Chinese with English abstract).

[10] Chi X Y, Yang Q L, Pan L J, Chen M N, He Y N, Yang Z, Yu S L. Isolation and characterization of fatty acid desaturase genes from peanut (L.)., 2011, 30: 1393–1404.

[11] Chi X Y, Zhang Z M, Chen N, Zhang X W, Wang M, Chen M N, Wang T, Pan L J, Chen J, Yang Z, Guan X Y, Yu S L. Isolation and functional analysis of fatty acid desaturase genes from peanut (L.)., 2017, 12: e0189759.

[12] Focks N, Benning C. Wrinkled1: a novel, low-seed-oil mutant ofwith a deficiency in the seed-specific regulation of carbohydrate metabolism., 1998, 118: 91–101.

[13] Pouvreau B, Baud B, Vernoud V, Morin V, Py C, Gendrot G, Pichon J P, Rouster J, Paul W, Rogowsky P M. Duplicate maizetranscription factors activate target genes involved in seed oil biosynthesis., 2011, 156: 674–686.

[14] 孙金波, 石素华, 杨利, 李凤丽, 王兴军, 赵术珍. 花生基因家族的全基因组与表达谱分析. 花生学报, 2020, 49(1): 9–18. Sun J B, Shi S H, Yang L, Li F L, Wang X J, Zhao S Z. Genome-wide analysis ofgene family and their expression profiles in peanut., 2020, 49(1): 9–18 (in Chinese with English abstract).

[15] Harwood H J. Oleochemicals as a fuel: mechanical and economic feasibility., 1984, 61: 315–324.

[16] 迟晓元, 郝翠翠, 潘丽娟, 陈娜, 陈明娜, 王通, 王冕, 杨珍, 梁成伟. 不同花生品种脂肪酸组成及其积累规律的研究. 花生学报, 2016, 45(3): 32–36. Chi X Y, Hao C C, Pan L J, Chen N, Chen M N, Wang T, Wang M, Yang Z, Liang C W. Fatty acid accumulation pattern in different types of peanut., 2016, 45: 32–36 (in Chinese with English abstract).

[17] Yu P, Wang C H, Xu Q, Feng Y, Yuan X P, Yu H Y, Wang Y P, Tang S X, Wei X H. Detection of copy number variations in rice using array-based comparative genomic hybridization., 2011, 12: 372.

[18] 陈娜, 迟晓元, 程果, 潘丽娟, 陈明娜, 王通, 王冕, 杨珍, 禹山林. 花生中低温胁迫相关转录因子基因的筛选. 核农学报, 2016, 30(1): 19–27. Chen N, Chi X Y, Cheng G, Pan L J, Chen M N, Wang T, Wang M, Yang Z, Yu S L. Profiling of genes encoding cold stress-related transcription factors in peanut., 2016, 30: 19–27 (in Chinese with English abstract).

[19] Barbour J A, Howe P R, Buckley J D, Bryan J, Coates A M. Cerebrovascular and cognitive benefits of high-oleic peanut consumption in healthy overweight middle-aged adults., 2017, 20: 555–562.

[20] 姜慧芳, 任小平, 黄家权, 雷永, 廖伯寿. 野生花生脂肪酸组成的遗传变异及远缘杂交创造高油酸低棕榈酸花生新种质. 作物学报, 2009, 35: 25–32. Jiang H F, Ren X P, Huang J Q, Lei Y, Liao B S. Genetic variation of fatty acid components inspecies and development of interspecific hybrids with high oleic and low palmitic acids., 2009, 35: 25–32 (in Chinese with English abstract).

[21] 于明洋, 孙明明, 郭悦, 姜平平, 雷永, 黄冰艳, 冯素萍, 郭宝珠, 隋炯明, 王晶珊. 利用回交法快速选育高油酸花生新品系. 作物学报, 2017, 43: 855–861. Yu M Y, Sun M M, Guo Y, Jiang P P, Lei Y, Huang B Y, Feng S P, Guo B Z, Sui J M, Wang J S. Breeding new peanut line with high oleic acid content using backcross method., 2017, 43: 855–861 (in Chinese with English abstract).

[22] 陈四龙. 花生油脂合成相关基因的鉴定与功能研究. 中国农业科学院研究生院博士学位论文, 北京, 2012, pp 60–61. Chen S L. Identification and Functional Analysis of Lipid Biosynthesis Related Genes in Peanut (L.). PhD Dissertation of Graduate School of Chinese Academy of Agricultural Sciences, Beijing, China, 2012, pp 60–61 (in Chinese with English abstract).

[23] Li F P, Ma C Z, Wang X, Gao C B, Zhang J F, Wang Y Y, Cong N, Li X H, Wen J, Yi B, Shen J X, Tu J X, Fu T D. Characterization of sucrose transporter alleles and their association with seed yield-related traits inL., 2011, 11: 168.

[24] Zaborowska Z, Starzycki M, Femiak I, Swiderski M, Legocki A B. Yellow lupine gene encoding stearoyl-ACP desaturase: organization, expression and potential application., 2002, 49: 29–42.

[25] Bruner A C, Jung S, Abbott A G, Powell G L. The naturally occurring high oleate oil character in some peanut varieties results from reduced oleoyl-PC desaturase activity from mutation of aspartate 150 to asparagine., 2001, 41: 522–526.

[26] 官梅, 李栒, 官春云. 利用基因芯片技术研究甘蓝型油菜油酸合成中差异表达基因. 作物学报, 2010, 36: 968–978. Guan M, Li X, Guan C Y. Differentially expressed genes in oleic acid synthesis ofby detected gene chip., 2010, 36: 968–978 (in Chinese with English abstract).

[27] 李玉兰, 孙勤富, 王幼平. 植物油脂合成的转录调控研究进展. 分子植物育种, 2016, 14: 2509–2518. Li Y L, Sun Q F, Wang Y P. Research advance in transcriptional regulation of lipid synthesis and accumulation in plant., 2016, 14: 2509–2518 (in Chinese with English abstract).

[28] Maeo K, Tokuda T, Ayame A, Mitsui N, Kawai T, Tsukagoshi H, Ishiguro S, Nakamura K. An AP2-type transcription factor, WRINKLED1, ofbinds to the AW-box sequence conserved among proximal upstream regions of genes involved in fatty acid synthesis., 2009, 60: 476–487.

[29] Fukuda N, Ikawa Y, Aoyagi T, Kozaki A. Expression of the genes coding for plastidic acetyl-CoA carboxylase subunits is regulated by a location-sensitive transcription factor binding site., 2013, 82: 473–483.

Expression pattern analysis of genes related to lipid synthesis in peanut

XU Jing1,**, PAN Li-Juan1,**, LI Hao-Yuan2, WANG Tong1, CHEN Na1, CHEN Ming-Na1, WANG Mian1, YU Shan-Lin1, HOU Yan-Hua2,*, and CHI Xiao-Yuan1,*

1Shandong Peanut Research Institute, Qingdao 266100, Shandong, China;2School of Marine and Technology, Harbin Institute of Technology, Weihai 264209, Shandong, China

In order to survey the regulation patterns of genes expression in the synthesis of oil at different developmental stages in peanut seed, F18 (high-oleic medium oil peanut variety) and ‘Luhua 6’ (low-oleic acid low oil variety) were used as research materials. The expression pattern analysis of genes was performed for peanut seeds on the 10, 30, 40, and 60 DAP, respectively. The results indicated that 130, 3556, and 2783 genes were significantly differentially expressed in the two varieties (lines) on the 30, 40, and 60 DAP, respectively. GO annotation and KEGG enrichment showed that DEGs were mainly enriched in the fatty acid synthesis and photosynthetic pathways, such as,,genes, which were involved in the accumulation of oleic acid. All the genes involved in photosynthesis pathway were photochlorophyll binding proteins and all of them were up-regulated. KEGG pathway indicated that all the genes involved in fatty acid biosynthesis pathway on the 40 DAP and 60 DAP were up-regulated. In summary, these results provide a theoretical basis for molecular study of fatty acid synthesis in the development stage of peanut seeds and offer some candidate genes as gene resources of quality improvement in peanut breeding.

peanut; microarray; differentially expressed genes; genes involved in lipid synthesis

10.3724/SP.J.1006.2021.04150

本研究由国家现代农业产业技术体系(花生)建设专项(CARS-13), 青岛市科技惠民示范引导专项(20-3-4-26-nsh), 山东省农业科学院农业科技创新工程(CXGC2016B02, CXGC2018E21), 山东省重大科技创新工程项目(2019JZZY010702), 广东省重点领域研发计划项目(2020B020219003), 山东省良种工程(2017LZGC003), 临沂市重点研发计划项目(2019YD009)和国家自然科学基金项目(31701464)资助。

This study was supported by the China Agriculture Research System (CARS-13), the Qingdao People’s livelihood Science and Technology Program (20-3-4-26-nsh), the Agricultural Scientific and Technological Innovation Project of Shandong Academy of Agricultural Sciences (CXGC2016B02, CXGC2018E21), the Breeding Project from Department Science & Technology of Shandong Province (2019JZZY010702), the Key-Area Research and Development Program of Guangdong Province (2020B020219003), the Improved Variety Engineering Project of Shandong Province (2017LZGC003), the Key Research and Development Program of Linyi (2019YD009), and the National Natural Science Foundation of China (31701464).

迟晓元, E-mail: chi000@126.com; 侯艳华, E-mail: houyanhua2006@163.com

**同等贡献(Contributed equally to this work)

许静, E-mail: xu_jing_cool@yeah.net

2020-07-10;

2020-11-13;

2020-12-23.

URL: https://kns.cnki.net/kcms/detail/11.1809.S.20201223.1136.006.html