植物特异性代谢产物对植物适应性及人类健康至关重要,但其生物合成的翻译调控机制尚不清楚。研究人员以拟南芥(Arabidopsis thaliana)为材料,揭示了一种由5′非翻译区(5′ untranslated region, 5′UTR)介导的翻译调控机制,该机制控制着硫代葡萄糖苷(glucosinolate)的生物合成。通过正向遗传筛选探索生长素(auxin)与硫代葡萄糖苷之间的代谢互作,研究人员鉴定了两个显性拟南芥等位基因,它们在MYB28(脂肪族硫代葡萄糖苷生物合成的关键调控因子)的5′UTR内相距13个碱基对(base pair, bp)的位置各携带一个单核苷酸替换突变。这些突变显著提高了MYB28蛋白丰度,而不影响其转录本水平,进而导致硫代葡萄糖苷产量增加。5′UTR的突变分析显示,RNA三级结构的改变影响了翻译效率,从而建立了RNA构象与代谢输出之间的联系。该研究发现揭示了植物特异性代谢中此前未被表征的转录后调控层次,并强调了5′UTR作为精准育种以增强作物性能和营养品质的潜在靶点。
植物产生数百万种特异性代谢产物,这些物质对植物生存至关重要,其中许多与人类健康密切相关。硫代葡萄糖苷便是一个典型例子,这类含氮、硫的化合物存在于十字花目(Brassicales)植物中,包括卷心菜、西兰花、萝卜、芥菜和羽衣甘蓝等重要农作物。硫代葡萄糖苷既是植物防御化合物,其水解产物对草食动物和病原体具有毒性,也是已知的促进健康化合物,具有抗癌和抗炎等活性。因此,提高硫代葡萄糖苷含量可同时增强十字花科作物的营养价值和抗病性。在拟南芥中,硫代葡萄糖苷生物合成途径已被广泛研究,其合成起始于细胞色素P450单加氧酶CYP79家族介导的氨基酸向相应醛肟(aldoxime)的转化。脂肪族硫代葡萄糖苷和吲哚硫代葡萄糖苷分别由链延伸的甲硫氨酸(Met)和色氨酸(Trp)衍生而来,分别受MYB28/MYB29/MYB76和MYB34/MYB51/MYB122等转录因子调控。特别值得注意的是,吲哚硫代葡萄糖苷途径与生长素生物合成存在代谢关联:吲哚-3-乙醛肟(indole-3-acetaldoxime, IAOx)既是吲哚硫代葡萄糖苷的中间体,也可通过替代途径转化为吲哚-3-乙酸(indole-3-acetic acid, IAA)这一主要内源生长素。当吲哚硫代葡萄糖苷生物合成下游步骤受阻时,IAOx代谢流会重定向至IAA合成,导致高生长素表型。然而,调控IAOx在防御代谢与生长素生物合成之间分配的机制仍有待阐明。
为探究IAOx与吲哚硫代葡萄糖苷生物合成之间的代谢交叉对话,研究人员以CYP83B1功能缺陷的ref5突变体为背景开展了无偏抑制子筛选。ref5突变体因IAOx积累而导致IAA水平升高,表现出典型的下胚轴伸长和叶片向下卷曲等表型。研究人员对约15万粒甲基磺酸乙酯(ethylmethanesulfonate, EMS)诱变的ref5种子进行了筛选,从约3万株M
2植株中分离出两株独立的抑制子品系#19-6和#30-3。这两株抑制子呈现出恢复的高生长素相关形态特征(restored high-auxin–associated morphological traits, rhax),其叶片形态与ref5明显不同,且IAA含量显著降低。遗传学分析表明,这两个抑制子突变均为单一位点的显性突变,且等位性测试证实它们极可能是同一基因的等位变异。
通过混池分离分析(bulked segregant analysis, BSA),研究人员将#19-6的突变定位至5号染色体上12个紧密连锁的候选单核苷酸多态性(single-nucleotide polymorphism, SNP)。尽管RNA测序(RNA sequencing, RNA-seq)显示#19-6相比ref5有广泛的转录组改变,但这12个候选基因本身的转录水平均未发生变化。进一步分析发现,AT5G61420(即MYB28)的共表达基因中有60%在#19-6中显著上调,提示该SNP可能激活了MYB28功能。测序确认#19-6在MYB28的第二转录本(MYB28.2)翻译起始密码子ATG上游152 bp处携带C-to-T突变;#30-3则在相隔13 bp的上游165 bp处携带另一个C-to-T突变。这两个突变均位于MYB28的5′UTR而非编码区,暗示其通过在翻译水平发挥调控作用。
研究人员将这两个显性等位基因分别命名为myb28-1D和myb28-2D。在野生型(wild-type, WT)背景下,myb28-1D和myb28-2D单突变体的形态与WT无差别,但脂肪族硫代葡萄糖苷主要成分3MSOP(3-甲硫基亚氧丙基硫代葡萄糖苷)、4MSOB(4-甲硫基亚氧丁基硫代葡萄糖苷)和8MSOO(8-甲硫基亚氧辛基硫代葡萄糖苷)的水平均显著高于WT,且该化学表型呈显性遗传。RNA-seq分析显示,ref5 myb28-1D中大多数脂肪族硫代葡萄糖苷生物合成基因显著上调。这些结果表明,MYB28 5′UTR中的SNPs足以独立于ref5突变增强脂肪族硫代葡萄糖苷产量,且这种增强效应不依赖转录水平改变,而可能通过促进MYB28蛋白翻译实现。
为验证上述假设,研究人员首先通过互补实验确定了MYB28.2是驱动脂肪族硫代葡萄糖苷生物合成调控的主要功能亚型:仅有MYB28.2能在myb28 myb29双突变体中恢复硫代葡萄糖苷产量,而MYB28.1不能。随后,研究人员构建了由2 kb天然启动子(含5′UTR)驱动的MYB28:HA标签融合蛋白表达载体,比较了野生型C等位基因(p2KbC::MYB28:HA)与突变型T等位基因(p2KbT::MYB28:HA)在myb28 myb29背景中的表达效果。RT-PCR显示两者转录本水平相当,但Western blot检测表明突变型T等位基因的MYB28-HA蛋白水平显著升高,相应地,脂肪族硫代葡萄糖苷水平也大幅提高。为进一步证明5′UTR本身足以影响翻译效率,研究人员构建了由35S启动子驱动、下游连接野生型或突变型MYB28 5′UTR与绿色荧光蛋白(green fluorescent protein, GFP)编码序列的报告基因体系。尽管GFP转录本水平存在变异,但携带突变型5′UTR的植株中GFP蛋白水平明确更高,证实5′UTR单独即可影响翻译效率。
关于5′UTR介导翻译调控的机制,研究人员注意到MYB28 5′UTR含有一个预测的上游开放阅读框(upstream open reading frame, uORF),而myb28-1D和myb28-2D突变均位于该uORF区域内。然而,这两个突变并未改变uORF的编码框长度,仅引起uORF编码肽的氨基酸替换(S31F和L27F),因此不太可能通过传统的uORF介导的核糖体滞留机制发挥作用。研究人员进一步利用CRISPR-Cas9系统构建了uORF下游3 bp缺失突变体myb28uORF-Δ3bp,该突变体同样表现出脂肪族硫代葡萄糖苷水平升高,且uORF翻译产物仅缩短一个氨基酸。
基于三个突变体(myb28-1D、myb28-2D和myb28uORF-Δ3bp)均不影响uORF编码肽长度的特点,研究人员推测MYB28 5′UTR的翻译调控可能依赖于RNA结构而非uORF肽序列本身。通过查询FOLDAtlas数据库中WT拟南芥的实验性硫酸二甲酯(dimethyl sulfate, DMS)探针数据,发现WT MYB28 5′UTR存在一个抑制翻译的假结(pseudoknot)结构。为验证该结构的功能,研究人员设计了两个突变体:A303G维持假结配对,A303G/U353C则通过破坏一个碱基对削弱假结。在WT拟南芥原生质体转染实验中,A303G突变体UTR的GFP报告基因翻译效率较低,而A303G/U353C突变体的翻译效率显著更高。利用RNAStructure程序进行的RNA结构预测显示,C241U(myb28-1D)和C228U(myb28-2D)突变改变了MYB28 5′UTR的整体结构并破坏了假结,有利于翻译进行;Δ3bp缺失突变同样破坏了假结结构。
为进一步巩固5′UTR结构与翻译效率之间的关联,研究人员设计了两组额外的突变体进行原生质体验证:C254G/U256A破坏假结且U256A引入提前终止密码子使uORF编码肽缺失最后三个残基;A352U则维持假结并形成不利于核糖体活性的长双链结构。结果显示,C241U具有最强的翻译效率,C254G/U256A也提高了翻译效率,两者均呈现相对松散的5′UTR RNA结构;而A352U突变因含有假结和长双链结构,翻译效率显著低于WT。特别重要的是,研究人员还构建了一个同义突变体U227G:该突变在保留uORF编码缬氨酸的前提下,改变了抑制性5′UTR结构,形成缺乏假结和长双链螺旋的整体开放构象,使GFP翻译效率提高约七倍。这一结果明确表明,翻译抑制可以独立于uORF编码肽,而通过5′UTR的结构特征实现。
在代谢互作方面,ref5 myb28-1D和ref5 myb28-2D原本因恢复ref5的高生长素表型而被筛选到。研究人员通过MYB28.2过表达证实,激活MYB28确实能够模拟ref5 myb28-1D的表型恢复效应,而MYB28.1过表达则不能。利用CRISPR-Cas9在ref5 myb28-1D背景中敲除MYB28后,rhax表型被逆转为高生长素表型,IAA水平恢复至ref5水平,证明MYB28的激活对于表型恢复是必要且充分的。进一步代谢分析显示,ref5 myb28突变体除脂肪族硫代葡萄糖苷升高外,主要吲哚硫代葡萄糖苷I3M(吲哚-3-基甲基硫代葡萄糖苷)水平也增加,同时IAOx和IAA水平降低,表明存在吲哚硫代葡萄糖苷生物合成与生长素稳态之间的代谢再平衡。RNA-seq和 quantitative PCR分析揭示,尽管CYP79B2和REF5等吲哚硫代葡萄糖苷生物合成基因未被显著诱导,REF2(CYP83A1)却被大幅上调。已知REF2虽偏好脂肪族醛肟,但也能以IAOx为底物,只是活性远低于REF5。因此,研究人员推测MYB28激活可能通过上调REF2及下游共同途径基因(如SUR1、UGT74B1),将IAOx代谢流从IAA转向吲哚硫代葡萄糖苷合成。这一假说得到遗传学验证:在ref5 myb28-1D和ref5背景中分别引入ref2突变后,IAOx和IAA水平均大幅回升且两者相当,高生长素表型也重新出现,表明REF2的功能对于MYB28激活降低生长素含量是不可或缺的。
在讨论部分,研究人员总结指出,植物特异性代谢产物的生物合成传统上被认为主要通过转录程序协调响应内外部刺激,而本研究鉴定了一种此前未知的转录后调控层次:由生物合成调控基因5′UTR中的细微序列变异介导的翻译调控。通过正向遗传筛选,研究人员发现MYB28 5′UTR中的两个显性突变通过增强翻译效率提高蛋白积累,而不改变转录本丰度,标志着翻译调控对植物特异性代谢的直接调控作用。
5′UTR通过uORFs、RNA结构和顺式调控基序等元件精细调节翻译。在MYB28案例中,尽管存在预测uORF,但突变分析表明所鉴定的突变通过破坏假结结构而非改变uORF编码肽来激活翻译。WT MYB28 5′UTR形成的假结结构抑制核糖体扫描,而多个破坏该结构的突变(myb28-1D、myb28-2D、Δ3bp、C254G/U256A、U227G)均增强翻译,维持或稳定假结的A352U突变则降低翻译效率,形成了结构-功能的明确对应关系。这种5′UTR介导的翻译激活导致MYB28蛋白大幅积累,引起脂肪族硫代葡萄糖苷过量合成,同时通过诱导REF2表达将共享中间体IAOx从IAA合成重新导向吲哚硫代葡萄糖苷生物合成,从而恢复生长素稳态。值得注意的是,REF2的底物杂泛性(promiscuity)在这种代谢分配中发挥关键作用,这种酶促灵活性可能是植物在进化过程中应对代谢分支点调控压力而形成的。
MYB28.2被确定为功能性亚型,而MYB28.1因缺失R2 DNA结合域可能仅在特定组织或条件下发挥作用。MYB28 5′UTR序列在十字花科物种中高度保守,myb28-1D和myb28-2D对应位置在所有被检物种中均保守,提示该调控机制可能具有广泛的育种应用潜力。在拟南芥中,破坏假结结构的单核苷酸替换可使脂肪族硫代葡萄糖苷产量提高达10倍,且在WT背景下不影响整体生长和形态。鉴于5′UTR介导的翻译调控提供了一条不依赖转录过表达即可控制酶丰度和激活生物合成过程的直接途径,它可能成为精准育种中增强作物防御和健康促进代谢物的高效靶标。
研究结论部分指出,本研究表明5′UTR介导的翻译调控在协调拟南芥MYB28表达与更广泛代谢网络中发挥重要作用。由MYB28翻译控制驱动的代谢流在生长(生长素)与防御(硫代葡萄糖苷)途径之间的动态再分配,展示了非编码序列如何贡献于平衡植物中相互竞争的生理需求。通过5′UTR中单碱基替换实现翻译激活的发现,为作物精准育种获取理想性状提供了一种变革性机制。