《Biology》：Integrated Transcriptomic and Metabolomic Analysis Reveals Key Genes in Fatty Acid and Flavonoid Biosynthesis in Macadamia Nut (Macadamia integrifolia)
编辑推荐：
bHLH转录因子在植物发育及脂肪酸和类黄酮代谢中处于核心地位。通过对两个脂肪酸含量存在显著差异的澳洲坚果品种‘桂热1号’（GR1）和‘A4’进行转录组学、代谢组学以及植物激素/元素谱分析，研究人员鉴定出142个MibHLH转录因子，将其划分为21个亚家族，其中
bHLH转录因子在植物发育及脂肪酸和类黄酮代谢中处于核心地位。通过对两个脂肪酸含量存在显著差异的澳洲坚果品种‘桂热1号’（GR1）和‘A4’进行转录组学、代谢组学以及植物激素/元素谱分析，研究人员鉴定出142个MibHLH转录因子，将其划分为21个亚家族，其中14个亚家族与拟南芥AtbHLH成员表现出共线性关系。MibHLH46呈现显著差异表达，其表达与脂肪酸生物合成基因的表达高度相关。启动子分析进一步揭示，脂肪酸和类黄酮生物合成基因的启动子区域富含MYC和G-box顺式作用元件。比较代谢组学分析显示，‘GR1’比‘A4’积累更高水平的多种矿质元素（Fe、Ca、Mg）和脱落酸（ABA）。qRT-PCR验证了转录组学趋势，确认了‘GR1’中与脂肪酸氧化、α-亚麻酸代谢及类黄酮途径相关基因的上调表达。研究人员提出，MibHLH46作为一个候选调控因子，通过转录调控关键生物合成基因参与核仁脂肪酸和类黄酮代谢，这为鉴定澳洲坚果脂质代谢中的调控基因奠定了基础。
**论文解读：澳洲坚果脂肪酸与类黄酮生物合成的转录调控机制研究**
**一、研究背景与科学问题**
澳洲坚果（*Macadamia integrifolia*）以其高油含量和丰富的有益不饱和脂肪酸而闻名，其核仁中单不饱和脂肪酸（MUFAs）尤其是油酸（C18:1）的比例极高，赋予其显著的营养与经济价值。然而，尽管脂肪酸组成对坚果品质至关重要，其合成的遗传调控网络，特别是转录水平上的精细调控机制，尚未被完全阐明。类黄酮作为植物重要的次生代谢产物，不仅参与坚果壳的形成，也可能与核仁品质相关。已有研究表明，碱性螺旋-环-螺旋（basic helix-loop-helix, bHLH）转录因子家族是真核生物中最大的转录因子家族之一，在植物发育、胁迫响应以及次生代谢，尤其是类黄酮和脂肪酸生物合成中扮演核心角色。然而，在澳洲坚果中，bHLH家族成员的全基因组鉴定、进化分析及其在脂肪酸和类黄酮代谢中的具体调控作用，此前缺乏系统研究。因此，有必要通过整合转录组与代谢组分析，筛选控制澳洲坚果脂肪酸合成的关键调控基因，为高品质品种的分子育种提供基因资源。
**二、研究内容与技术方法**
该研究发表于《Biology》期刊。研究人员选用来自中国广西亚热带农业科学研究所澳洲坚果种质资源圃（106°79′22″E, 22°34′13″N）的两个脂肪酸含量差异显著的澳洲坚果品种‘GR1’（高脂肪酸和高类黄酮含量）和‘A4’（低含量）作为材料，开展系统比较分析。研究主要采用了以下关键技术方法：
1. **转录组测序（RNA-Seq）**：利用Illumina平台对两个品种的核仁进行转录组测序，获取全局基因表达谱，并通过HISAT2比对到澳洲坚果参考基因组（GCA_013358625.1），使用StringTie进行转录本组装，鉴定差异表达基因（DEGs）。
2. **非靶向代谢组学分析**：采用液相色谱-质谱联用（LC-MS）技术对核仁代谢物进行提取、鉴定与定量，结合OPLS-DA模型筛选差异代谢物并富集KEGG通路。
3. **植物激素与矿质元素定量**：使用电感耦合等离子体光谱法（ICP）测定矿质元素，利用液相色谱-串联质谱法（LC-MS/MS）定量内源植物激素。
4. **生物信息学分析**：利用MUSCLE进行多序列比对，IQ-TREE构建最大似然系统发育树，PlantTFDB数据库进行转录因子注释，MEME进行保守基序分析，并开展启动子区顺式作用元件分析。
5. **qRT-PCR验证**：以澳洲坚果苹果酸脱氢酶基因*MiMADH*为内参，采用2
-ΔΔCT法对关键候选基因的表达进行独立验证。
**三、主要研究结果**
**3.1 两品种表型与代谢特征比较**
营养分析显示，‘GR1’核仁的粗蛋白（8.76 g 100 g
-1 DW）、粗脂肪（79.16% DW）和总糖（32.94 mg g
-1 DW）含量均显著高于‘A4’（分别为6.66 g 100 g
-1 DW、72.14% DW和30.07 mg g
-1 DW），而灰分含量‘GR1’（1.14%）显著低于‘A4’（1.23%）。元素谱分析表明，‘GR1’的铁（Fe, 26.74 mg kg
-1）、钙（Ca）、镁（Mg）、磷（P）含量分别是‘A4’的2.34倍、1.30倍、1.14倍和1.09倍。植物激素分析显示，‘GR1’中脱落酸（ABA）含量高达1279.07 ng g
-1，约为‘A4’的48.95倍，差异极为显著。此外，‘GR1’中多种细胞分裂素和赤霉素（GA3、GA4、GA7）含量也显著高于‘A4’。
**3.2 转录组与代谢组差异图谱**
转录组PCA分析显示，两个品种沿第一主成分（解释96.52%变异）明显分离。共鉴定到3465个DEGs（1667个上调，1798个下调）。GO富集显示DEGs富集于“蛋白二聚化活性”、“磷酸酯水解酶活性”等分子功能，以及“核糖核苷单磷酸生物合成过程”等生物学过程。KEGG富集显示“脂肪酸降解”、“糖酵解/糖异生”、“丙酮酸代谢”等脂质代谢相关通路显著富集。代谢组学共检测到3667种代谢物，其中494种上调、965种下调。KEGG分析显示苯丙烷生物合成、类黄酮生物合成等次级代谢通路显著差异。
**3.3 MibHLH转录因子家族全基因组鉴定**
从澳洲坚果基因组中鉴定出142个MibHLH转录因子，分为21个亚家族。基因结构分析显示这些基因具有多样的外显子-内含子组织和保守蛋白基序。重复事件分析表明，WGD/片段重复占16.9%，串联重复占8.5%。共线性分析显示14个MibHLH成员与拟南芥AtbHLH存在共线性关系。
**3.4 IIIf亚家族成员MibHLH46与类黄酮途径的关系**
表达分析显示，*MibHLH46*在‘A4’中高表达而在‘GR1’中低表达，差异显著。启动子分析表明，*MibHLH46*启动子主要含有光响应和SA响应元件。而脂肪酸和类黄酮生物合成关键基因（如*PAL*、*4CL*、*CHS*、*F3H*、*FLS*、*UGT*）的启动子区域富含MYC和G-box等bHLH结合元件，提示MibHLH蛋白可能通过结合这些元件调控靶基因转录。
**3.5 qRT-PCR验证**
qRT-PCR验证了‘GR1’中脂肪酸氧化相关基因（*FAO-MFP-a-like*、*FAO-MFP2-like*、*FAO-MFP2-like2*、*FALDH*家族等）、α-亚麻酸代谢基因（*ALA-βOx-pMFP*家族、*JMT2*、*LA-9LOX5*）的显著上调。类黄酮上游基因*PAL*、*PAL2*和*4CL*在‘GR1’中分别上调1.7倍、2.5倍和4.9倍，但关键酶*CHS*和*FLS*无显著差异，表明该途径并非在所有节点协同上调。
**3.6 bHLH介导的脂肪酸生物合成分子通路**
综合转录组和代谢组数据，构建了涉及苯丙氨酸起始的类黄酮生物合成与脂肪酸代谢之间交叉调控的分子网络，*PAL*、*4CL*、*CHS*、*F3H*、*FLS*等基因在两个品种间表现出显著不同的表达模式。
**四、讨论与结论**
讨论部分指出，bHLH转录因子可能通过与MYB和WD40蛋白形成MBW复合物，直接调控脂肪酸和类黄酮途径基因的表达。‘GR1’中显著升高的Fe、Ca、Mg可作为脂质代谢酶的辅因子，而极高的ABA水平可能与种子成熟和胁迫响应相关，共同参与调节脂质生物合成通量。*MibHLH46*启动子虽以光响应和SA响应元件为主，但E-box基序属于G-box家族，可能间接参与ABA调控。研究同时指出，仅依赖‘GR1’和‘A4’两个品种的比较，难以完全排除遗传背景差异的影响，未来需要通过多品种比较或基因过表达/敲除实验进行功能验证。
研究结论认为，共鉴定到142个MibHLH转录因子，分属21个亚家族。*MibHLH46*呈现显著差异表达，其表达与脂肪酸合成相关基因密切相关，并通过qRT-PCR得到验证。研究人员提出*MibHLH46*可能通过转录调控关键生物合成基因（如*MiPAL*）参与类黄酮和脂肪酸合成，尽管其功能仍需进一步验证。这一发现为澳洲坚果营养品质改良的分子育种提供了重要候选基因和理论依据。
