《Genes》：Transcriptional Response of Rice Mesocotyl Elongation to Sowing Depth and Identification of Key Regulatory Factors
Ya Wang,
Dong Liu,
Mengjuan Ma,
Ming Li,
Jing Fu,
Fengjiang Yu,
Qiulin Li,
Yuetao Wang,
Fuhua Wang and
Haiqing Yin
+ 1 author
编辑推荐：
本研究针对直播水稻中胚轴（Mesocotyl）长度（ML）与播种深度的响应关系，利用转录组学（RNA-seq）、加权基因共表达网络分析（WGCNA）和层次聚类分析（HCA）等技术，鉴定出苯丙烷生物合成、丝裂原活化蛋白激酶（MAPK）信号通路和二萜生物合成等关键通路及多个候选基因（如OsCCR1、LOC_Os07g03319、OsCYP71Z21），为解析水稻中胚轴伸长分子机制和培育耐深播水稻品种提供了理论基础。
背景：耐深播——水稻直播技术推广的关键“卡脖子”难题
在全球人口增长和环境变化的双重压力下，水资源节约、劳力集约的水稻直播技术正受到越来越多的关注。然而，与传统的育苗移栽相比，直播技术，特别是旱直播，在实际生产中面临着一大难题：播种过深会导致出苗率低、出苗不齐、幼苗活力差，严重制约了其规模化推广。这背后的核心生理学原因，与水稻幼苗的“中胚轴”息息相关。中胚轴是连接幼芽和种子根之间的一个关键胚胎组织，其伸长是推动幼苗顶出土壤、成功出苗的主要动力。研究表明，水稻的出苗能力与中胚轴的长度呈显著正相关。换句话说，拥有更长中胚轴的水稻品种，能更好地适应较深的播种深度，实现更高的出苗率和更健壮的幼苗。因此，如何从分子层面揭示水稻中胚轴响应播种深度、调控其自身伸长的机制，成为了作物科学家亟待破解的关键科学问题，也是培育耐深播水稻品种、保障国家粮食安全的重要突破口。
此前的研究已揭示了植物激素（如赤霉素、乙烯、油菜素内酯等）和一些关键基因（如SD1、GY1、OsGSK2等）在中胚轴伸长调控中的核心作用。然而，面对从浅播到深播的连续梯度变化，驱动中胚轴伸长动态调控的核心转录组图谱、关键生物合成与信号通路，以及其中的枢纽基因仍不甚清晰。为此，来自国内科研团队的研究人员，以粳稻品种郑稻209为材料，通过设置3厘米、5厘米和7厘米三个播种深度梯度，巧妙地构建了中胚轴长度显著递变的实验模型。他们利用高通量转录组测序（RNA-Seq）技术，结合生物信息学分析方法，系统绘制了水稻中胚轴响应播种深度的动态基因表达谱，旨在鉴定调控中胚轴伸长、从而影响水稻耐深播能力的关键通路和基因。这项研究的结果最终以“Transcriptional Response of Rice Mesocotyl Elongation to Sowing Depth and Identification of Key Regulatory Factors”为题，发表于学术期刊《Genes》。
关键研究方法
研究采用了粳稻品种郑稻209作为实验材料，在生长室中设置3、5、7厘米三个播种深度进行培养。主要技术方法包括：
- 1.
表型与细胞学观测：测量了播种8天后幼苗的中胚轴长度，并利用显微镜观察了中胚轴细胞长度。
- 2.
转录组测序与差异分析：收集8天龄幼苗的中胚轴组织，提取总RNA，构建cDNA文库，在Illumina NovaSeq 6000平台上进行RNA-Seq。对获得的序列数据进行质量控制，并比对到水稻参考基因组。通过DESeq2软件，以|log2FC| > 1和FDR < 0.05为标准，识别了不同播种深度间的差异表达基因。
- 3.
生物信息学与共表达分析：对差异表达基因进行了基因本体（Gene Ontology, GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路富集分析。使用加权基因共表达网络分析，构建了基因共表达网络，以鉴定与中胚轴长度高度相关的基因模块。通过层次聚类分析，对差异表达基因的表达模式进行了无监督聚类。
- 4.
实时荧光定量PCR验证：针对转录组分析筛选出的关键候选基因，设计了特异性引物，利用实时荧光定量PCR技术，验证了其在不同播种深度下的表达趋势，以确认RNA-Seq数据的可靠性。
研究结果
3.1. 中胚轴长度随播种深度增加而增加
研究发现，郑稻209的中胚轴长度随着播种深度（从3厘米到7厘米）的增加而显著增加。在7厘米深度时达到最长，最长可达3.4厘米。细胞学观察也表明，中胚轴细胞长度也呈现相同的变化趋势。这证实了在较深的播种压力下，中胚轴及其细胞的伸长是水稻幼苗为适应深播环境而启动的一种关键形态建成策略。
3.2. 转录组测序与差异基因鉴定
RNA-Seq分析产生了高质量的数据。通过比较不同深度组（W3 vs. W5, W3 vs. W7, W5 vs. W7），共识别出上千个差异表达基因。其中，三个比较组中共有60个共同差异表达基因。富集分析显示，这些差异基因主要富集在苯丙烷生物合成、角质、木栓质和蜡质生物合成、植物丝裂原活化蛋白激酶信号通路、二萜生物合成、氰基氨基酸代谢、光合生物的碳固定、类黄酮生物合成和谷胱甘肽代谢等通路。
3.3. 共表达网络与关键通路分析
WGCNA分析将差异表达基因聚类为不同的模块，并发现蓝色模块和棕色模块的基因与中胚轴长度的变化高度相关。对这两个模块的基因进行KEGG富集分析，进一步将核心调控通路聚焦于苯丙烷生物合成、碳代谢、光合作用天线蛋白和植物-病原体互作。HCA进一步将差异表达基因分为8个表达谱，揭示了基因表达随播种深度变化的动态模式。其中，在浅播条件下高表达、在深播条件下被抑制的基因（Profile 0和Profile 1）富集在光合作用天线蛋白和光合作用通路，表明浅播下中胚轴的能量供给充足。而在中等深度（5厘米）表达达到峰值的基因（Profile 3）则富集在苯丙烷生物合成通路，这可能是在适度深播下平衡细胞壁机械强度与延展性的关键。
3.4. 核心通路中关键基因的表达模式
研究者对富集到的核心通路中的关键基因进行了具体分析。在苯丙烷生物合成相关基因中，随着播种深度增加，木质素合成关键基因OsCCR1的表达下调，而编码腐胺羟基肉桂酰基转移酶的OsPHT3和LOC_Os04g59260基因的表达则上调。在MAPK信号通路相关基因中，LOC_Os07g03319和LOC_Os07g03580的表达上调，而LOC_Os07g03409的表达显著下调。在二萜生物合成相关基因中，随着播种深度增加，OsCYP76M5和OsCYP71Z2的表达下调，而OsCYP71Z21的表达上调。这些基因的表达趋势通过qRT-PCR实验得到了验证，表明它们可能是响应播种深度、调控中胚轴伸长的关键候选基因。
结论与意义
本研究通过系统的转录组学分析，构建了水稻中胚轴响应播种深度梯度的动态基因表达图谱，并揭示了调控中胚轴伸长的关键分子网络。研究表明，中胚轴的伸长是一个由多种代谢和信号通路协同调控的复杂过程，涉及细胞壁修饰、能量代谢、激素信号整合等多个层面。其中，苯丙烷生物合成通路（特别是木质素代谢分支）、MAPK信号通路和二萜生物合成通路（核心是赤霉素生物合成分支）是三个关键的调控枢纽。深播诱导的黑暗信号可能通过调节MAPK级联反应，进而影响植物激素信号（如乙烯、赤霉素、生长素等），并重塑二萜代谢谱，最终精细调控赤霉素的生物合成。同时，深播条件下木质素合成相关基因（如OsCCR1）的下调，可能降低了细胞壁的刚性，增加了细胞壁的延展性，从而促进了细胞的伸长。此外，研究还鉴定出了一批具有不同深度响应模式的转录因子，为后续功能验证和调控网络解析提供了丰富的候选靶点。
这项研究的重要意义在于，它不仅为理解水稻中胚轴适应深播胁迫的分子机制提供了新见解，更重要的是，鉴定出的多个候选基因（如OsCCR1、LOC_Os07g03319、OsCYP71Z21等）为后续的基因功能验证和分子设计育种提供了直接靶点。通过分子标记辅助选择或基因编辑技术，利用这些关键基因的优异等位变异，有望快速培育出具有更长中胚轴、更强耐深播能力的水稻新品种，从而解决水稻直播技术推广中的“出苗关”瓶颈，为实现节水增效、机械化种植的现代稻作生产模式提供重要的种质资源与理论支持。
