环境驱动的转录组重编程调控稻米直链淀粉含量与垩白度的可塑性

《The Crop Journal》：Environment-driven transcriptomic reprogramming modulates the plasticity of amylose content and chalkiness in rice

　　

摘要译文：理解环境波动如何塑造稻米品质性状（如直链淀粉含量（AC）和垩白度（CD））的可塑性，对于培育气候韧性的优质品种至关重要。为解决这一问题，研究人员在两个可塑性品种（RD15和HXZ9）及一个稳定性品种（D50）中，于重叠的光温组合条件下进行了种植，并在籽粒灌浆期采用了时间分辨转录组学和全长异构体测序。对遗传变异、基因表达和替代剪接（AS）的整合分析揭示了一种基本的调控策略二分法：可塑性品种通过核心代谢的定向重编程和关键调控因子的条件敏感性替代剪接实现适应性，而稳定性品种则通过一组组成型表达的、参与精确RNA加工和细胞区室化的基因来维持稳健性。值得注意的是，研究人员鉴定并通过实验验证了FLO12的差异转录本使用与AC可塑性强烈相关，这为异构体水平调控与表型变异之间提供了直接的机制联系。该多组学框架从基因组变异到动态转录和转录后控制，定义了调控AC和CD环境响应性与稳定性的层级调控景观，为气候韧性优质稻米的精准育种提供了可行靶点。

论文解读：环境驱动的转录组重编程调控稻米直链淀粉含量与垩白度的可塑性

研究背景：稻米品质高度受环境条件影响，尤其是温度。灌浆期高温会增加垩白度、降低整精米率、降低表观直链淀粉含量（AC）并损害淀粉糊化特性；低温则可提高某些品种的AC并减少垩白，但常降低整精米率。太阳辐射与温度交织影响品质。表型可塑性是指同一基因型在不同环境中产生不同表型的能力，作物育种中可塑性与稳定性的目标常常对立，需采取不同策略：最小化可塑性培育稳健品种，或通过选择环境特异性等位基因培育适应特定环境的可塑品种。然而，以往研究多集中于关联定位或单个基因验证，缺乏对不同可塑性品种基因组的系统比较，也缺乏在多个灌浆阶段和环境条件下的时间分辨转录组与全长转录组整合分析。替代剪接（AS）可扩展蛋白质组多样性并快速响应环境，但AS如何在水稻籽粒灌浆期调控品质可塑性仍不清楚。因此，本研究旨在通过整合基因组变异、基因表达和AS分析，揭示AC和垩白度（CD）可塑性的层级调控机制。

研究内容：研究人员选取了三个水稻品种：RD15（优质泰国香米，AC/CD环境可塑）、HXZ9（优质高产籼稻常规品种，AC/CD可塑）和D50（热带粳稻对照品种，品质稳定）。它们在三个重叠的光温环境下种植：三亚两季（Sy1: 2023年8–11月；Sy2: 2024年1–4月）和杭州一季（Hz: 2024年6–9月）。RD15种植于Sy1和Sy2，HXZ9种植于全部三个环境，D50种植于Sy2和Hz，各有三个生物学重复。对发育中的籽粒在05、10、15、20、25 DAF（花后天数）五个阶段取样进行表型测定、RNA-seq和ONT全长转录组测序。同时采用PacBio HiFi和ONT平台对三品种进行全基因组重测序，鉴定SNP、InDel和结构变异（SV）。

关键方法：本研究主要采用了以下关键技术：①全基因组重测序（PacBio HiFi + ONT长读长测序）用于鉴定三品种间遗传变异；②时间分辨RNA-seq（短读长）分析基因表达动态；③ONT全长转录组测序构建泛转录组并检测AS事件；④DESeq2鉴定差异表达基因（DEGs）和稳定表达基因（SEGs），SUPPA2基于PSI值鉴定差异转录本使用（DTU）和稳定转录本使用（STU）；⑤WGCNA共表达网络分析和GRNBoost2预测转录因子-靶基因调控网络；⑥qRT-PCR验证关键异构体比例。样本来源为上述三个环境种植的三品种籽粒。

研究结果（按原文小标题）：

3.1 稻米品质表型可塑性与遗传变异：表型追踪显示籽粒呈双相生长模式（05–10 DAF快速膨大，后渐缓至成熟）。环境数据表明Sy2的Tmax高于Sy1，Hz温度变幅最大；Sy2的ΔT大于其他环境；Sy2和Hz的日日照时长长于Sy1。品质评估显示RD15和HXZ9的AC、CD和垩白粒率（CGR）具有环境可塑性，而D50在Sy2与Hz间保持稳定。全基因组测序显示，RD15携带最多的SNP、InDel和SV，其次是HXZ9和D50。两可塑品种（RD15与HXZ9）间的遗传变异数量大于任一与稳定品种D50间的变异。在已知AC/CD相关基因中，RD15与HXZ9间差异最大（影响130个AC基因和128个CD基因的SNP、InDel及SV），且Wx和Chalk5的变异数目呈现明显梯度。

3.2 籽粒发育过程中的转录组动态与功能可塑性：转录组分析显示所有品种均存在一致的早期（05–15 DAF）到晚期（20–25 DAF）的阶段性转换。环境驱动的转录组分歧主要表现为20 DAF后同环境样本聚类。Mfuzz聚类获得五个表达簇，C1和C2在各品种/环境中高度保守（富集异源物质转运、蛋白泛素化、细胞定位、碳水化合物代谢等），C3、C4、C5则显示品种与环境特异的功能富集，C5几乎每个品种-环境组合都有独特组合。

3.3 水稻籽粒发育期间泛转录组动态景观：构建的泛转录组包含76,596个基因和112,053个转录本，其中77.4%为新转录本。泛转录组转录本平均长度短于参考注释，多异构体基因比例更高。AS分析检测到37,760个事件，内含子保留（RI）占比最高（36.59%），其次为A3SS（24.89%）和A5SS（15.26%）。各品种在不同环境下有独特的环境特异转录本集，并富集不同通路（如RD15的Sy1特异转录本富集翻译，Sy2特异转录本富集RNA加工）。跨环境共有转录本核心为37,698（Sy1）至45,927（Sy2）个。大量转录本显示严格时期特异性（41,750仅在单一时期表达），其中05和25 DAF特异转录本占比超70%。

3.4 品种AC和CD可塑性与稳定性的不同转录策略：RD15在Sy1 vs Sy2中DEGs最少（7857上调、6848下调）但SEGs最多（4696）；D50在Sy2 vs Hz中DEGs最多（21,250总数）但SEGs较少（1406）。HXZ9与D50在Sy2 vs Hz比较中DEG和SEG比例无显著差异。GO富集显示RD15和HXZ9的DEGs富集羧酸/有机酸代谢，D50的DEGs富集磷脂代谢和DNA损伤响应；D50的SEGs富集甲基化、RNA加工/修饰、胞内运输、蛋白定位和去磷酸化，而RD15/HXZ9的SEGs富集RNA修饰和定位。相关性分析（|r|>0.7, P<0.05）分别鉴定出RD15的910个AC相关和1105个CD相关可塑性基因，HXZ9的930个AC相关（Sy2 vs Hz）等，D50的62个AC相关和118个CD相关稳定性基因。

3.5 品种AC和CD可塑性与稳定性的不同转录后策略：HXZ9在Sy1 vs Hz中DTU事件最少（2534上调、2761下调）但STU最多（31,080）；D50在Sy2 vs Hz中DTU最多（9976），同时维持36,143个STU。HXZ9与D50在Sy2 vs Hz比较中，HXZ9的DTU比例显著低于D50（OR=0.961, P=0.012），STU比例无显著差异。DTU基因（DTUGs）富集于羧酸代谢、有机酸代谢、RNA加工、剪接、胞内蛋白运输和蛋白定位；STU基因（STUGs）富集于翻译、酰胺生物合成、氨基酸代谢和胞内运输。D50贡献了最大的稳定性相关转录本集（3073个AC相关、3100个CD相关稳定转录本），包含多个已知品质基因异构体。

3.6 品质可塑性与稳定性的调控网络及候选基因：整合基因和转录本水平数据，RD15获得1248个AC可塑性和1578个CD可塑性基因；HXZ9获得1272个AC可塑性（Sy2 vs Hz）和1694/2006个CD可塑性基因；D50获得2168个AC稳定性和2280个CD稳定性基因。WGCNA鉴定出AC关键模块：RD15的MEturquoise（547基因，富集蛋白翻译和能量代谢）、HXZ9的MEblue（335基因，富集嘧啶化合物代谢）、D50的MEturquoise（943基因，富集细胞组分组装和肽生物合成）。TF调控网络显示AC可塑性网络共享flo4。FLO12作为HXZ9与D50共有的AC基因，在5和10 DAF高表达（>500 TPM），其中一个特异转录本在HXZ9的15 DAF表现为DTU，而D50在各阶段均为STU。FLO12在HXZ9与D50间存在2个内含子SNP和3个5′UTR InDel。qRT-PCR验证了该特异异构体比例与RNA-seq PSI值一致，且与AC表型显著正相关（r=0.71, P=0.005）。CD关键模块中，RD15和HXZ9的CD可塑性模块（Sy1 vs Sy2）富集氨基酸代谢和翻译，D50的CD稳定性模块富集RNA加工。TF网络显示OsPK1基因在10 DAF表达与CD显著负相关（r=?0.79, P=0.03），但因其存在大量遗传变异而被排除为候选基因。

讨论与结论：研究发现，可塑性品种通过定向重编程核心能量代谢和关键调控因子（如FLO12）的条件敏感性AS实现适应性，而稳定品种D50则通过组成型稳定表达的基因套组（涉及表观遗传调控、RNA加工和细胞区室化）辅以主动的代偿性转录响应来维持稳定性。FLO12的差异转录本使用与AC可塑性之间的直接联系，提供了异构体水平调控决定品质性状变异的机制证据。该多组学框架揭示了从基因组变异到动态转录和转录后控制的层级调控景观，为培育气候韧性优质稻米提供了可行靶点。研究也存在局限性，如实验设计不对称、缺乏多年环境重复、相关分析不能证明因果关系、调控网络为计算预测等，未来需通过CRISPR异构体特异性敲除、近等基因系或多亲本群体进行验证。研究发表在《The Crop Journal》。