MADS51功能性单倍型FH2介导野生稻开花期分化并促进生态物种形成

《Cell Reports》：A MADS51 allele is associated with flowering-time divergence and contributes to ecological adaptation during wild rice speciation

　　

开花期分化（flowering-time divergence）可促进生态物种形成（ecological speciation），但野生植物中的因果基因和突变仍大多未知。研究人员调查了野生稻物种 Oryza rufipogon 与 O. nivara 之间的开花期分化。通过建立野生稻转化体系（transformation system），并结合基因组扫描（genome scans）、群体基因组学（population genomics）和功能实验（functional assays），研究人员鉴定出 MADS51 是开花期分化的一个贡献基因。在 O. rufipogon 中敲除（knockout）MADS51 加速开花，而过表达（overexpression）延迟开花。一个功能性单倍型（functional haplotype 2 [FH2]）由三个非同义替换（nonsynonymous substitutions）和一个4-bp移码（4-bp frameshift）定义，在东南亚 O. nivara 中近乎固定（nearly fixed）。互补实验（complementation assays）证实，FH2 比祖先单倍型（ancestral haplotypes）促进更早开花。群体遗传分析（population genetic analyses）揭示与 FH2 相关的选择清除（selective sweep），其似乎是在 O. nivara 中新发（de novo）产生的。这些发现为理解新发突变（de novo mutation）如何促成野生稻表型分化（phenotypic differentiation）提供了见解。

研究背景与问题：物种形成是进化生物学核心议题，生殖隔离限制基因流并促进群体分化；植物中开花期分化是最常见的合子前生殖隔离（prezygotic reproductive isolation，受精前障碍）形式。尽管QTL定位和群体基因组扫描发现大量开花期候选位点，但在原生物种中经转基因功能验证者少，因果突变或等位基因更少。Oryza rufipogon 与 O. nivara 是最近缘野生稻和亚洲栽培稻直接祖先；O. nivara 适应干旱生境，近期多次从 O. rufipogon 起源，包括南亚Niv1和东南亚Niv2；早花有助于逃避季节性干旱并促进与 O. rufipogon 的生殖隔离，是平行生态物种形成的经典案例。此前杂交群体定位到十余个开花期QTL，但因果基因和突变不明，主要瓶颈是野生稻转化困难。研究人员聚焦东南亚起源，整合基因组扫描、候选基因筛选、群体基因组学和功能实验，解析开花期分化的遗传基础。

研究开展、结论与意义：研究人员建立野生稻转化体系，鉴定MADS51为开花期分化贡献基因。敲除MADS51使 O. rufipogon 早花，过表达延迟开花；FH2由3个非同义替换（nonsynonymous substitutions，导致氨基酸改变的碱基替换）和1个4-bp移码（frameshift，改变阅读框的插入缺失）定义，在东南亚 O. nivara 中近乎固定，并比祖先单倍型促进更早开花；FH2伴随选择清除（selective sweep，有利突变受选择后导致邻近多样性降低），似在 O. nivara 中新发产生。意义在于首次在原生野生稻中直接验证MADS51自然变异的功能，揭示de novo突变可促成生态适应和生殖隔离，为植物适应、物种形成及栽培稻分子改良提供线索。

关键技术方法：选择NEP0401（O. rufipogon）和NEP0233（O. nivara）自交6代降低杂合度，优化组织培养和农杆菌介导转化，建立野生稻转化体系，转化效率约75%和55%。群体分析采用已发表全基因组重测序数据：Dataset A含69个Ruf2与64个Niv2（其中39份有普通园表型：18 Ruf2、21 Niv2），Dataset B含348份代表4个遗传谱系和9个亚群。方法包括F_ST（群体分化指数）、D_XY（群体间绝对序列分歧）、XP-CLR（跨群体复合似然比）基因组扫描，314个功能已知开花期基因与654个选择区域基因交集，CRISPR/Cas9敲除、过表达、FH1/FH2互补、qRT-PCR、亚细胞定位，以及单倍型网络、Fay and Wu’s H、EHH（扩展单倍型纯合度）和est-sfs分析。

研究结果：
1. 建立O. rufipogon/O. nivara转化体系：通过筛选8份代表材料、选择NEP0401和NEP0233并自交6代，全基因组重测序确认杂合度0.347%和0.343%；优化萌发与组织培养，获得约75%和55%转化效率；为在原生物种中验证基因功能提供工具。
2. 两个野生稻物种开花期分化与群体遗传结构：133份东南亚材料普通园实验中，Niv2于7月中至9月中开花，Ruf2于10月底至11月底开花，超过1个月，近乎完全合子前隔离；ADMIXTURE和邻接树形成两个簇，F_ST=0.384，Niv2的π（nucleotide diversity，核苷酸多样性）=0.0020略低于Ruf2的0.0023；支持Niv2起源于Ruf2且开花期高度分化。
3. 全基因组选择清除检测与候选开花期基因鉴定：F_ST与D_XY初筛得1,282个区域25.8 Mb（8.29%），XP-CLR得686个区域11.9 Mb（3.83%），交集485个高置信选择清除区域5.784 Mb（1.55%）含654个注释基因；GO富集生殖发育、碳水化合物合成和金属离子跨膜运输等；将314个功能已知开花期基因与654个基因交集得11个高优先候选，FTL4、HUB1和MADS51同时具π降低、Niv2中Fay and Wu’s H显著负及基因F_ST居全基因组前5%的信号；MADS51含5个2-kb上游/UTR变异和4个编码变异（3个非同义替换和1个4-bp移码），分化强，故选其验证。
4. MADS51功能验证：在O. rufipogon NEP0401中CRISPR/Cas9敲除exon 1产生提前终止等位基因，田间自然生长季中敲除株比野生型早花>20天，过表达株晚花5-11天，表明MADS51抑制开花；敲除株茎秆更直立、每穗小穗数增加，趋向O. nivara表型；栽培稻Hwayoung转移DNA插入突变体在长日照延迟10天、短日照延迟7天，显示背景依赖。
5. MADS51序列变异及其对开花期影响：基于de novo组装和群体重测序，MADS51有23个SNP（17个在2-kb上游/UTR，6个编码）和6个indel（5和1）；9个显著分化变异中5个调控、4个编码；调控变异未与预测转录因子结合位点重叠，叶片MADS51表达无显著种间差异；编码区π降低且Fay and Wu’s H显著负；四个编码变异为S61C（位于保守MADS-box结构域，与DNA结合和转录调控相关）、G65D、S97Y和4-bp缺失（Wei et al. (2021)报道与水稻抽穗期相关），在Niv2中独特衍生；4-bp缺失在终止密码子上游造成11个氨基酸C端延伸但不影响核定位；四变异紧密连锁，个体贡献难分。
6. 单倍型分析与正选择检验：基于4个候选功能变异定义9个功能性单倍型（FH），Niv2主要固定FH2（Ruf2中仅2.2%），FH1和FH3仅见于Ruf2；普通园中FH1和FH3分别221和216天开花，FH2为125天；互补实验中FH2平均104天早于FH1的114天，确认FH2促早花，但10天差异远小于自然约3个月，提示MADS51效应中等；FH2组MADS51周围多样性降低（S=186，π=0.0008 vs FH1组S=238，π=0.0015），衍生等位基因频谱呈U形，EHH衰减更慢，支持选择清除和FH2为衍生适应等位基因。
7. MADS51等位基因起源与扩散：用348份Dataset B构建单倍型网络，5个携带FH2的单倍型形成独立且高支持组，主要分布于东南亚Niv2（91.06%），仅Hap1含2个Ruf1、3个Ruf2和1个Niv1；该组与最近祖先Hap6至少相隔5个突变步；FH2主要见于Niv2a和Niv2c；支持FH2在Niv2群体中新发产生而非来自祖先Ruf2，并促进O. nivara生态适应。

讨论总结与结论翻译：讨论指出，MADS51变异贡献O. rufipogon与O. nivara开花期分化；转化体系突破野生稻功能验证瓶颈，首次提供MADS51调控野生稻自然开花期变异的直接证据。敲除MADS51早花>20天并产生O. nivara样直立茎和更多小穗，过表达晚花5-11天；栽培稻中功能丧失反而延迟开花，说明功能受遗传背景和互作网络影响。MADS51仅解释自然种间差异的小部分，开花期是多基因性状，先前QTL研究分别检测到5和11个开花期QTL。FH2在Niv2固定并促早花，互补实验显示104 vs 114天；强选择信号支持其为适应衍生等位基因；MADS51多效性可能促成开花期、株型和小穗数相关变化，连锁亦可能参与。起源分析支持FH2在东南亚Niv2中新发，证据包括主要限于Niv2（91.1%）、与最近祖先相差5步、周围强选择；这提示生态物种形成适应等位基因并非均来自standing genetic variation（祖先已有遗传变异），de novo突变在强定向选择下可快速适应和成种。局限为MADS51仅一个组分、四变异个体贡献未定、背景依赖机制未明。研究结论可译为：这些发现为理解新发突变（de novo mutation）如何促成野生稻表型分化（phenotypic differentiation）提供见解；总体而言，研究人员的发现支持一个模型，即在东南亚O. nivara中新发产生的早花MADS51等位基因促进O. nivara对新环境的生态适应。