IF50.5｜隐秘变异通过层级上位性驱动番茄花序形态演化

四个基因的隐秘变异组合出216种基因型，揭示层级上位性如何塑造表型空间。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

隐秘遗传变异单独效应微弱，但被认为通过上位性相互作用构成驱动性状演化的巨大遗传多样性库。泛基因组学揭示了基因家族、顺式调控区和调控网络中普遍存在的变异，但隐秘变异如何推动表型多样化仍缺乏系统解析。番茄花序结构在茄科内呈现丰富自然变异，许多突变体模拟了跨物种的自然变异，且常涉及隐秘遗传互作。此前已在番茄中鉴定出两个MADS-box转录因子J2和EJ2之间的剂量依赖冗余上位性，它们单独突变时对花序分枝无表型效应，但组合突变导致高度分枝。这一J2–EJ2旁系同源对为研究隐秘等位基因在调控网络中的上位性互作如何影响性状可演化性提供了理想系统。

科学问题

卡在哪里： 此前研究局限于单个物种或少数近缘种，且主要揭示主效变异，隐秘等位基因在自然群体和基因调控网络中的作用未被系统探索。此外，背景依赖性（可能部分源于隐秘等位基因）在进化与发育遗传学中常见，但缺乏对隐秘变异及其在表型演化中作用的系统解析。J2–EJ2旁系同源对虽提供了研究平台，但需鉴定更多网络组分才能充分发挥其潜力。

本文要回答： 利用自然和工程化的顺式调控隐秘变异，鉴定番茄花序发育调控网络中的额外冗余反式调控因子，并系统解析四个网络基因的遗传互作如何塑造花序分枝表型空间。

技术路线

作者主要用了：CRISPR–Cas9基因编辑、泛基因组变异挖掘、双荧光素酶报告、RNA-seq、F2群体数量遗传学建模。

作者首先利用番茄泛基因组数据挖掘EJ2启动子区变异，结合开放染色质和预测转录因子结合位点，筛选出两个野生种特有的顺式调控变异（EJ2Sh和EJ2Sp）。通过CRISPR–Cas9在j2突变体背景下工程化一系列EJ2启动子等位基因，验证了这些顺式调控变异的表型效应。随后，基于预测的DOF和PLT结合位点，结合表达图谱筛选候选转录因子，并通过双荧光素酶实验证实PLT3激活EJ2启动子。进一步构建plt3 plt7双突变体，发现其与j2 ej2双突变体表型相似，且PLT3/7调控EJ2和J2表达。最后，将六个EJ2pro j2系与plt3 plt7/+杂交，产生六个F2分离群体，获得216种基因型，在多个季节和环境中共计35,606个花序中量化分枝事件，构建高分辨率基因型–表型图谱。

图 1　Fig. 1。图1展示了茄属植物花序复杂性增加，以及J2和EJ2旁系同源对在控制番茄花序结构中的剂量依赖冗余关系。图中显示ej2 j2双突变体产生增殖的分生组织，而单个突变体无表型。EJ2启动子上游6 kb区域显示开放染色质、保守非编码序列、预测的TFBS和泛基因组变异。S. lycopersicum和S. habrochaites的EJ2顺式调控区比较显示一个DOF和两个PLT结合位点，重叠的DOF-PLT位点被自然变异破坏。

核心亮点

亮点 1｜EJ2启动子多个顺式调控区控制花序分枝

通过泛基因组挖掘和CRISPR–Cas9工程化，作者在EJ2启动子区鉴定出多个调控花序分枝的顺式调控元件。自然变异EJ2Sh和EJ2Sp分别来自S. habrochaites和S. pennellii，它们单独存在时几乎不引起分枝，但在j2突变背景下分别导致轻微和显著的分枝增加。工程化的七个EJ2pro等位基因中，两个约100 bp的缺失（EJ2pro7和EJ2pro8）移除了DOF-PLT结合位点，在j2背景下引起连续范围的分枝效应，而在J2功能背景下均无表型。此外，靶向另外两个启动子区域产生的缺失等位基因中，仅一个约3.3 kb的缺失导致轻度分枝。这些结果表明EJ2启动子多个区段共同正调控EJ2表达，控制花序发育。

图 2　Fig. 2。图2显示EJ2启动子153 bp靶区域，包含DOF和PLT结合位点及四个CRISPR–Cas9 gRNA。两个强表型系EJ2pro7 j2和EJ2pro8 j2携带约100 bp重叠缺失，移除了DOF-PLT结合位点。图中展示不同EJ2pro等位基因在j2背景下的分枝表型，显示连续范围的分枝效应，而在J2背景下无表型。

亮点 2｜PLT3和PLT7冗余调控EJ2及花序分枝

基于顺式调控变异破坏的DOF和PLT结合位点，作者聚焦于两个PLT旁系同源基因PLT3和PLT7。双荧光素酶实验显示PLT3强烈激活完整EJ2启动子，但不能激活EJ2pro8、EJ2pro-Sh和EJ2pro-Sp等突变启动子。CRISPR突变体分析表明plt7单突变体与野生型相似，plt3单突变体产生少量分枝，而plt3 plt7双突变体表现出极端的茎端分生组织过度增殖和分枝。这些突变等位基因也显示剂量依赖冗余：plt3/+ plt7基因型显示弱分枝，plt3 plt7/+基因型显示中等分枝。RNA-seq分析显示J2和EJ2在plt3 plt7分生组织中表达极低，而PLT3和PLT7在j2 ej2分生组织中仍表达，支持PLT调控SEP基因。

图 3　Fig. 3。图3显示EJ2、J2、PLT3和PLT7在五个发育阶段的表达动态，PLT3和PLT7在花转变后与J2和EJ2相似表达。PLT3和PLT7蛋白的系统发育树显示它们与拟南芥AtPLT3和AtPLT7是直系同源，但来自独立复制。双荧光素酶实验显示PLT3激活完整EJ2启动子，但不激活突变启动子。plt3 plt7双突变体表现出极端分枝，且plt3/+ plt7和plt3 plt7/+基因型显示剂量依赖效应。RNA-seq PCA显示plt3 plt7和j2 ej2双突变体聚集在花分生组织成熟阶段。

亮点 3｜旁系同源对内部存在超乘性协同上位性

通过构建216种基因型的F2群体并量化35,606个花序的分枝事件，作者发现突变在不同位点间倾向于乘性组合（对数链接的加性模型解释59.40%偏差，优于同一性链接的加性模型46.05%）。进一步拟合成对互作模型显著优于乘性模型（解释83.11%偏差，P<10^{-16}），检测到普遍的上位性，43/80个上位性项P<0.05。其中plt3与plt7之间以及EJ2pro等位基因与j2之间存在强烈的超乘性协同互作，组合突变产生的平均分枝事件数比单个突变效应相乘的预期高出3至9倍。

图 4　Fig. 4。图4展示了PLT3/7–EJ2/J2遗传和分子网络的简化架构，以及六个F2群体的基因型–表型图谱。图中显示216种基因型的花序分枝事件分布，以及不同plt3 plt7背景下j2和EJ2pro突变组合的对数表型关系，斜率随PLT突变累积而降低。层级上位性模型的响应面显示旁系同源对内部的协同效应和对间的掩蔽效应，例如plt3 plt7/+在WT背景下有81.30倍效应，但在EJ2pro8 j2背景下仅2.49倍。

亮点 4｜旁系同源对之间存在剂量依赖的掩蔽互作

作者发现非旁系同源对之间的加性×加性上位性项中11/14为负且P<0.05，表明存在系统性掩蔽效应。通过绘制不同plt3 plt7突变背景下j2和EJ2pro突变组合的对数表型关系，发现斜率随PLT3/PLT7突变累积而降低，几乎降至零。同样，在不同j2和EJ2pro8背景下分析plt3 plt7突变组合，斜率随EJ2 J2扰动强度增加而降低。这种模式在所有六个EJ2pro等位基因中一致，表明一个旁系同源对中的突变线性地重新标度另一个旁系同源对中突变的效应，即一个对的突变效应在另一个对高度突变时被削弱。

亮点 5｜层级上位性模型以单参数捕获掩蔽互作

作者提出层级上位性模型，将旁系同源对内部的任意显性和上位性模式与对间的多线性互作相结合，仅用一个参数描述一个旁系同源对突变对另一个对效应的掩蔽程度。该模型解释了82.14%的偏差，与成对互作模型（83.11%）相当，但参数数量大幅减少（AIC差值-19.64，层级模型约18,358倍更可能）。留一季交叉验证的平均R²为0.89，与成对模型相同，表明层级模型在保持高预测力的同时更简洁。响应面可视化显示，随着一个旁系同源对的表型效应增加，另一个对的突变效应逐渐减弱，例如plt3 plt7/+在WT EJ2 J2背景下有81.30倍效应，但在EJ2pro8 j2背景下仅2.49倍。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： CRISPR–Cas9基因编辑、泛基因组变异挖掘、双荧光素酶报告、RNA-seq、F2群体数量遗传学建模
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事数量遗传学、进化发育生物学或基因调控网络建模的研究者，尤其是关注隐秘变异和上位性的人。

带走一句

层级上位性模型用极简参数捕获复杂遗传互作，提示在分析多基因性状时，应考虑旁系同源对内部协同与对间掩蔽的嵌套结构，而非仅依赖成对互作模型。