通过选择多个位点塑造金鱼草（Antirrhinum）中的发育梯度

《SCIENCE ADVANCES》：Shaping of developmental gradients through selection on multiple loci in Antirrhinum

　　

发育依赖于分子梯度的精确塑造，但自然选择如何建立这种精确性尚不清楚。本研究分析了控制金鱼草（Antirrhinum）两个变种之间黄色花颜色梯度差异的基因。研究人员表明，这些差异部分取决于色素生物合成基因FLAVIA（FLA）的顺式调控变异（cis-regulatory variation）。FLA与另外三个位点相乘性地相互作用，其中一个是FLA的反式作用调控因子（trans-acting regulator），进一步塑造了黄色梯度。所有位点在杂交带（hybrid zone）中表现出渐变（cline），其宽度与表型效应相关，显示了选择如何通过作用于多个位点的顺式和反式变异（cis and trans variation）以显著精度打磨梯度形状。

论文解读文章

**研究背景**
发育过程中模式形成的关键在于分子梯度的精确塑造，例如早期胚胎中区域特性的建立或成年个体中着色模式的生成。尽管梯度的可重复性和精确性在多种生物系统中被广泛观察，但自然选择如何通过调整梯度特性（如陡度、范围）来实现这种精确性，仍是一个未解之谜。杂交带（hybrid zone）——物种或变种之间发生长期基因交流的区域——为研究这一问题提供了机会。在杂交带中，受选择的位点可通过等位基因频率的陡峭渐变（cline）被识别。例如，蝴蝶翅膀的警戒色或花的引导图案相关基因已被发现受选择影响。然而，这些位点如何相互作用以精细调控颜色梯度，尚缺乏系统研究。本研究以金鱼草（Antirrhinum）的两个变种——A.m.m. var. pseudomajus和A.m.m. var. striatum——为模型，探究了黄色花颜色梯度差异的遗传基础及其选择机制。这两个变种分布在比利牛斯山附近的杂交带两侧，具有互补的花色模式：前者在洋红色背景上具有黄色斑块，后者在黄色背景上具有洋红色斑块，且黄色梯度陡度显著不同。研究人员旨在揭示自然选择如何通过多个位点的顺式和反式调控变异来塑造这种发育梯度。

**研究内容与结论**
研究人员通过全基因组F_ST扫描、RNA测序（RNA-seq）、原位杂交（in situ hybridization）以及渐变群拟合（cline fitting）等技术，鉴定了控制黄色梯度的关键基因及其相互作用方式。主要结论如下：黄色梯度的陡度差异部分由FLAVIA（FLA）基因的顺式调控变异决定，该基因编码查尔酮4′-O-葡萄糖基转移酶（chalcone 4′-O-glucosyltransferase, 4′CGT），参与黄色色素（aurone）的合成。FLA与SULFUREA（SULF）基因相乘性地相互作用，后者通过小RNA在反式水平抑制FLA表达。此外，AURINA（AUN）和CREMOSA（CRE）两个位点也以类似方式参与调控。所有四个位点（FLA、SULF、AUN、CRE）的表型效应在杂交带中与渐变宽度显著相关，效应越大的位点渐变越窄，表明选择强度与表型效应大小成正比。总体而言，黄色梯度是通过多个位点相乘性相互作用（multiplicative interaction）被精确塑造的，而非单个位点的修饰。

**重要性与发表刊物**
该研究揭示了自然选择如何通过作用于多个位点的顺式和反式变异，以极高精度打磨发育梯度的形状，为理解进化中的“修补”（evolutionary tinkering）机制提供了新证据。论文发表在《SCIENCE ADVANCES》。

**主要关键技术与方法**
研究人员利用全基因组F_ST扫描识别受选择位点，通过RNA-seq比较两个变种及基因渗入系中的转录水平差异，采用RNA原位杂交可视化FLA表达的梯度分布。表型分析基于ImageJ测量花瓣上黄色强度（绿光减去蓝光反射率）沿近端-远端纵轴的变化。渐变群拟合采用对称S形模型（sigmoid cline）估计等位基因频率的渐变宽度和中心，并使用模拟退火算法优化参数。样本来源包括西班牙比利牛斯山杂交带（Planoles）的野生种群，以及实验室构建的F₂、F₃和F₄群体。

**研究结果**

**Yellow gradient steepness varies between varieties**
通过测量黄色强度沿花瓣纵轴的变化，发现A.m.m. var. pseudomajus的黄色梯度陡峭，而A.m.m. var. striatum的梯度平缓。这种差异与洋红色素的分布有关，陡峭梯度可减少黄色与洋红色的重叠，增强视觉对比。

**An island of high relative divergence is located on chromosome 2**
全基因组F_ST扫描在染色体2上发现一个高分化岛（F_ST island），其大小约为789 kb，远大于已知的ROS/EL位点岛。该岛位于着丝粒附近，重组率极低（约0.1 cM/Mb），因而范围较大。该岛不包含已知花色基因，提示存在新位点。

**The chromosome 2 island harbors a novel flower color locus, FLA**
通过F₂和F₄群体的表型排名与基因型关联分析，发现该岛上的SNP与黄色强度显著相关。将striatum来源的等位基因命名为FLAs，pseudomajus来源的命名为fla^p。FLAs在黄色高表达组中富集，表明该位点控制黄色色素分布。

**Both cis and trans variation controls the gradient of FLA expression**
RNA-seq和原位杂交显示，FLA编码4′CGT，其表达在striatum中呈平缓梯度，在pseudomajus中呈陡峭梯度。通过基因渗入实验分离顺式和反式效应：fla^p在顺式水平限制FLA表达并增强梯度陡度，而SULF^p在反式水平通过小RNA抑制FLA，两者相加效应更强。

**FLA and SULF interact multiplicatively to shape the yellow gradient**
对比不同基因型组合的黄色强度纵轴曲线，发现fla^p和SULF^p各自减小黄色强度的效应沿纵轴呈下降梯度，且二者相乘预测的效应与实际双突变体效应高度吻合，证明它们相乘性地相互作用以陡化梯度。

**Selection intensity correlates with effects of FLA and SULF alleles on the yellow gradient**
在杂交带中，FLA和SULF均呈现等位基因频率渐变。SULF^p的渐变宽度约1 km，fla^p约3 km，表明选择系数相差4–16倍。进一步分析重组等位基因FLA^sp（具有striatum上游和pseudomajus下游区域），其渐变宽度约9 km，表型效应更小，验证了选择强度与表型效应正相关。

**Two additional loci shape the yellow gradient through multiplicative interactions**
AUN和CRE位点对黄色梯度的影响较小，但同样产生负斜率，陡化梯度。将四个黄色位点（AUN、CRE、FLA、SULF）的pseudomajus等位基因相乘预测的黄色梯度与实际pseudomajus梯度在误差范围内一致，而加性模型在远端偏离，支持相乘性模型。

**总结讨论**
研究人员翻译了研究结论部分：
我们的结果表明，黄色发育梯度的精确打磨是通过选择作用于多个位点而非单个位点的修饰来实现的，这与进化修补（evolutionary tinkering）的概念一致。单个等位基因的表型效应可以非常微小，但足以产生选择所维持的陡峭渐变，从而在杂交种群之间保持差异。适应度差异依赖于与其他位点（尤其是控制洋红色的位点）的上位相互作用（epistatic interactions）。因此，根据遗传背景的不同，不同的梯度可能被选择所青睐。