LAX1调控水稻穗发育转录网络的独立验证与精修

《Plants》：Independent Validation and Refinement of the LAX1-Regulated Transcriptional Network in Rice Panicle Development

　　

摘要翻译：
水稻穗部结构对籽粒产量起决定性作用，然而调控腋生分生组织形成的调控网络仍未得到充分阐明。先前研究将LAX1鉴定为一个典型的bHLH（基本螺旋-环-螺旋）转录因子，并详细描述了其在lax1突变体中的直接靶标及全基因组转录变化。本研究利用CRISPR/Cas9技术独立构建了LAX1敲除突变体，并对幼穗进行了转录组测序。lax1-KO突变体表现出次生分枝减少和籽粒数下降。采用更严格的阈值（|log₂倍数变化（log₂FC）| ≥ 1）后，共鉴定出518个高置信度差异表达基因（DEGs），其中131个与先前报道的DEGs重叠（方向一致性为88.5%），从而支持了核心LAX1调控网络。分析揭示了PIN家族成员特异性的表达调控（OsPIN1c/1d和OsPIN2下调；OsPIN3t和OsPIN9上调），激素相关基因（OsPP2C09、OsPP2C49、OsRR3、OsRR9、OsGASR9）和OsTPP基因（OsTPP1、OsTPP4、OsTPP9）的表达变化，以及穗部调控因子（IPA1、OsTB1、DEP1、DEP3、CUC和MADS-box基因）的表达改变。基序富集分析进一步发现bHLH基序是出现频率最高的类型，而CAMTA和FRS/FRF基序在DEG启动子中具有最高的统计显著性，提示可能存在多层次转录调控结构。综上，这些发现验证并精修了LAX1调控网络，将LAX1定位为水稻穗部发育的关键协调因子。

论文解读：

水稻（Oryza sativa L.）是全球半数以上人口的主粮，其穗部结构直接决定籽粒产量。穗发育由茎端分生组织（SAM）向花序分生组织（IM）转变，进而启动腋生分生组织（AM）形成分枝和小穗。该过程受到转录因子和植物激素（尤其是生长素和细胞分裂素）复杂网络的调控。尽管已鉴定出多个参与穗发育的调控基因，如LAX PANICLE 1（LAX1）、LAX PANICLE 2（LAX2）、MONOCULM 1（MOC1）、FRIZZY PANICLE（FZP）、IDEAL PLANT ARCHITECTURE 1（IPA1）等，但该调控网络的层级组织和分子互作仍不清楚。LAX1编码bHLH转录因子，是AM起始的关键调控因子。先前研究虽已鉴定LAX1的直接靶标，但对其丧失功能后如何重塑特定激素信号分支、尤其是生长素转运蛋白网络及下游分生组织身份调控因子的理解仍不完整。为此，研究人员利用CRISPR/Cas9技术独立构建了LAX1敲除突变体，并对幼穗进行转录组测序，以期在更严格阈值下独立验证并扩展LAX1调控网络的转录组框架。

研究人员以粳稻品种中花11（ZH11）为遗传背景，通过CRISPR/Cas9靶向LAX1唯一外显子，获得了一个纯合lax1-KO系（第51 bp处单个C碱基缺失导致移码和提前终止）。脱靶分析确认前五个预测脱靶位点无突变；qRT-PCR证实LAX1表达显著降低。表型分析显示，lax1-KO株系与non-KO对照在穗长和一次枝梗数上无差异，但二次枝梗数和每穗实粒数显著减少，与LAX1促进枝梗分生组织起始的已知功能一致。随后对lax1-KO和non-KO的幼穗（≤5 mm）进行RNA-seq，每个基因型三个生物学重复。测序数据质量高，各样品Q30均超过93.21%，超过91.10%的clean reads唯一比对到水稻基因组；生物学重复间Pearson相关系数≥0.90，主成分分析（PCA）显示基因型依赖的明显聚类。以FDR < 0.05且|log₂FC| ≥ 1为阈值，共鉴定出518个DEGs，其中402个上调、116个下调。GO富集分析显示DEGs主要涉及结合、催化活性、膜组分及代谢过程；KEGG富集分析表明植物-病原互作、植物激素信号转导、内质网蛋白质加工、半乳糖代谢及淀粉和蔗糖代谢等通路显著富集。为验证RNA-seq结果，研究人员对15个基因（包括LAX2、MOC1、DEP1、DEP3、IPA1、APO2、RCN4、OsSHI1、G1等9个已知穗部调控因子和6个随机基因）进行qRT-PCR，结果与RNA-seq趋势基本一致。

在激素相关基因方面，生长素信号组分在差异表达转录本中占显著比例。OsPID（水稻中拟南芥PINOID激酶直系同源物）在lax1-KO中下调，与LAX1作为OsPID转录激活因子的已知作用一致。研究人员进一步检测了PIN家族成员的表达，发现OsPIN1c、OsPIN1d和OsPIN2显著下调，而OsPIN3t和OsPIN9上调；OsPIN1a、OsPIN1b、OsPIN5b和OsPIN8无显著变化。值得注意的是，四个OsPIN1旁系同源物中仅在幼穗中特异性表达的OsPIN1c/1d下调，而主要参与根和茎发育的OsPIN1a/1b不变。此外，OsGH3.1和OsGH3.3下调，OsIAA14、OsIAA23和OsIAA24也发生显著变化；ABA相关基因OsPP2C09和OsPP2C49、细胞分裂素响应因子OsRR3和OsRR9、赤霉素诱导基因OsGASR9均表现出显著表达变化。在淀粉和蔗糖代谢通路中，多个基因在lax1-KO中下调，包括OsTPP4、OsTPP1、OsTPP9、OsHXK7、OsTPS8及若干葡聚糖内切-1,3-β-葡萄糖苷酶基因；其中OsTPP4和OsHXK7在本研究及先前报道中均一致下调，提示LAX1与海藻糖-6-磷酸（Tre6P）信号之间存在保守联系。

在已知穗部调控因子方面，研究人员检测了41个已表征穗发育相关基因的表达，发现8个基因（IPA1、OsSHI1、LAX1、DEP1、DEP3、OsTB1、G1和APO2/RFL）发生差异表达。qRT-PCR验证确认IPA1和OsTB1下调，DEP1和DEP3上调；OsCUC1和OsCUC3显著下调；OsMADS4、OsMADS51和OsMADS56表达降低；APO2/RFL、FZP和APO1也下调；而OsbHLH067、OsbHLH068、OsbHLH069和LAX2表达无显著变化。为探究可能的合作转录因子，研究人员对518个DEGs的启动子进行了基序富集分析。在763个植物TF结合基序中，共有173个基序显著富集。最显著富集的基序包括CAMTA2、CAMTA3、FAR1、FHY3和CAMTA1，属于钙调蛋白结合转录激活因子（CAMTA）和FRS/FRF转录因子家族；其余前列基序主要与BES/BZR和bHLH家族相关。bHLH家族是出现频率最高的类型，而CAMTA和FRS/FRF家族基序统计显著性最高。配对相似性分析显示，前13个位置权重矩阵（PWM）存在高度冗余的DNA结合偏好。

讨论部分指出，lax1-KO表型与中等至强等位基因一致，确认其功能缺失状态。本研究鉴定的DEG数（518）少于Fu等人报道的2100个，可能源于取样时期较晚（≤5 mm vs ≤2 mm）、更严格的倍数变化阈值（2倍 vs 1.5倍）、不同突变等位基因以及更保守的统计校正。与Fu等数据集的DEG重叠分析显示，518个DEGs中有131个（25.3%）重叠，方向一致性88.5%，两者在核心转录趋势（OsPID下调和激素及代谢通路富集）上高度一致，支持核心LAX1调控网络的稳健性。PIN基因成员特异性表达变化与功能相关，ospin1c/1d双突变体可模拟lax1穗部缺陷；OsPIN2下调可能影响分蘖和株高，OsPIN3t和OsPIN9上调可能代表对生长素运输破坏的补偿反应。Tre6P通路相关基因的保守下调提示LAX1可能通过影响糖信号调节分生组织活性。IPA1和OsTB1的协同下调与IPA1直接调控OsTB1的已知关系一致，而DEP1上调可能涉及其他因子或阶段特异性效应。APO2/RFL和CUC基因的下调可能反映AM起始受损的反馈响应。MADS-box基因下调中，OsMADS51已被证明是LAX1直接靶标。基序富集分析支持多层次转录架构：直接层由LAX1的bHLH结构域介导，间接层通过与其他TF家族（包括bHLH和CAMTA）合作实现。但间接调控关系仍属假设，需通过ChIP-qPCR、酵母双杂交和遗传上位性实验等直接证据加以验证。

研究结论：本研究独立验证并精修了LAX1调控水稻穗发育的转录网络。518个DEGs揭示多个激素信号通路相关表达变化，尤其是PIN基因成员特异性表达改变；OsTPP基因（含OsTPP1和OsTPP9）下调提示LAX1可能与Tre6P信号相关；IPA1、OsTB1、DEP1、DEP3、CUC和MADS-box等分生组织身份基因表达改变。基序富集分析表明bHLH基序最为常见，CAMTA和FRS/FRF基序统计显著性最高，符合可能的多层次转录调控结构。这些发现将LAX1定位为通过多级转录级联协调穗发育多方面过程的重要调控因子，为未来解析穗部结构组合调控密码提供了基础。