《The Plant Genome》：Traversing the effects of ploidy changes in different Eragrostis curvula genotypes through high-throughput RNA sequencing
编辑推荐：
本研究旨在探究基因组加倍（多倍体化）对画眉草（Eragrostis curvula）转录组的影响及其对牧草品质的调控机制。研究人员对从二倍体到八倍体的10个不同倍性基因型的叶片组织进行了高通量RNA测序分析。研究发现，随着倍性水平升高，差异表达基因（DEGs）数量增加，基因调控网络发生重塑。在高倍性水平下，多个与胁迫耐受、表观遗传调控、木质素生物合成、细胞壁重塑以及多糖代谢相关的基因表达发生改变。研究结果表明，多倍体化可能通过精细的转录调控，在增强植株胁迫耐受性的同时，可能降低了牧草消化性。该研究为理解多倍体化的分子机制提供了新见解，对培育平衡胁迫耐受性和生物质消化性的牧草品种具有指导意义。
在许多植物中，拥有额外染色体组（多倍体化）是一种普遍现象，它被视为基因组进化和物种形成的主要驱动力。多倍体植物通常展现出更强的适应性、活性和对不良环境的耐受性，这些特性使其在农业和生态学上具有重要价值。然而，基因组加倍这把“双刃剑”在带来优势的同时，也可能引发复杂的转录组重塑，影响植物的其他重要农艺性状，例如作为饲料作物的营养价值。画眉草（Eragrostis curvula），一种多年生C4禾本科牧草，因其在自然界中天然存在从二倍体到八倍体的不同倍性水平，成为研究多倍体化效应的绝佳模型。尽管画眉草在干旱、贫瘠土壤中表现出色，但其饲料品质，特别是消化性，却受到高纤维含量等因素的限制。那么，基因组加倍究竟如何重塑画眉草的转录图谱？这种重塑又如何影响其胁迫耐受性和作为饲料的消化性？为了解开这些谜团，一支由D. F. Santoro、J. Carballo等人组成的研究团队开展了一项深入的研究，相关成果发表在了《The Plant Genome》期刊上。
为了系统探究倍性变化对画眉草转录组的影响，研究人员采用了高通量RNA测序（RNA-Seq）技术作为核心方法。他们选取了10个具有不同倍性水平（二倍体到八倍体）和繁殖模式的画眉草基因型，每个基因型设置三个生物学重复，共计30个叶片样本。提取总RNA并质检合格后，使用Illumina NovaSeq X Plus平台进行双端测序。获得的序列数据与画眉草参考基因组进行比对，并利用DESeq2软件进行差异表达基因（DEGs）分析。此外，研究还综合运用了多种生物信息学方法，包括基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析、转录因子（TF）家族鉴定、加权基因共表达网络分析（WGCNA）以及SWItchMiner（SWIM）网络分析，以全面解析多倍体化引起的转录组变化、关键通路、调控网络及核心调控基因。
3.1 测序与定量
研究共获得了高质量的测序数据，主成分分析（PCA）显示样本能按倍性水平清晰聚类，表明转录组存在倍性依赖性的差异。
3.2 差异表达基因的鉴定
差异表达分析揭示，相对于二倍体对照，随着倍性水平（4x, 6x, 7x, 8x）升高，独特差异表达基因（DEGs）的数量显著增加，从4x的42个激增至8x的1518个，表明转录组随基因组加倍发生了大规模重编程。在最高的8x水平，下调基因数量超过上调基因，暗示存在部分剂量补偿效应。研究人员还鉴定出一个由433个基因组成的核心集合（Ploidy_vs_2x DEGs），它们在所有多倍体基因型中均呈现一致的差异表达模式。
3.3 多倍体化影响的GO功能分类
基因本体富集分析显示，不同倍性水平激活了独特的生物学过程。例如，4x水平独特富集了响应渗透胁迫等术语；6x水平与磷酸化、镁离子结合等相关；7x水平涉及DNA结合、响应脱落酸和寒冷；8x水平则与高尔基体、液泡、蛋白质泛素化等功能联系紧密。核心的433个基因则与线粒体、响应生长素和脱落酸等通路相关。
3.4 与画眉草倍性水平相关的代谢通路
KEGG通路分析表明，随着倍性升高，涉及的代谢通路愈发复杂。在6x和7x水平，“二苯乙烯、二芳基庚烷类和姜辣素生物合成”通路被富集；8x水平独特富集了“错配修复”、“同源重组”和“DNA复制”等与基因组稳定性相关的通路。核心基因集则显著富集“植物-病原体互作”通路。
3.5 转录因子家族的鉴定
研究发现，特定转录因子家族的表达随倍性变化而改变。例如，FAR1家族成员在多个高倍性水平上调。在8x水平，与胚胎发育和分生组织维持相关的转录因子（如PDF2、KNAT1、KNAT6、YABBY1）表达上调。此外，与染色质重塑（如HAC1、SUVH1）和泛素化过程（如BPM1）相关的转录因子也被激活。统计分析显示，GRAS转录因子家族在核心的倍性敏感基因集中被显著富集。
3.6 共表达网络分析
通过加权基因共表达网络分析（WGCNA），研究识别出49个基因共表达模块，其中多个模块与倍性水平显著相关。这些模块中富含与木质素生物合成、细胞壁重塑、胁迫响应以及表观遗传调控相关的基因。例如，与细胞壁修饰相关的基因（如HST、CAD7、XTH32）上调，而与果胶降解相关的基因（如F21M12.28）下调，这暗示高倍性可能导致生物质更难以消化。同时，多种胁迫响应基因（如RBOHC、GSTU19）和表观遗传调控基因（如NRPD1、FDM1、EZA1）也特异性上调。进一步的SWIM网络分析预测了潜在的“开关基因”，其中内切葡聚糖酶GH9B16和漆酶LAC7的下调，可能与生物质糖化效率降低有关。
本研究通过对不同倍性画眉草基因型的转录组进行系统性分析，揭示了多倍体化引起的复杂而有序的分子重塑网络。主要结论如下：首先，转录组重编程的程度与倍性水平正相关，高倍性引发更广泛的基因表达变化，并在最高倍性（8x）出现剂量补偿迹象。其次，多倍体化差异性地激活了与胁迫响应、细胞壁生物合成与重塑、激素信号、表观遗传修饰等相关的特定通路和转录因子家族。尤为重要的是，研究发现了一种潜在的“权衡”机制：高倍性在上调一系列胁迫耐受相关基因（如抗氧化、病原防御、泛素-蛋白酶体系系统相关基因）和增强细胞壁结构（如木质素和细胞壁重塑基因上调）的同时，可能导致了牧草消化性和生物质糖化效率的下降。这种权衡通过共表达网络分析和关键“开关基因”（如GH9B16）的鉴定得到了进一步支持。
这项研究的意义深远。它不仅深化了我们对植物多倍体化这一重要进化事件分子机制的理解，揭示了其在增强环境适应性方面的代价，而且为画眉草及其他饲料作物的分子育种提供了宝贵的遗传资源和明确的靶点。未来，通过基因编辑或分子标记辅助选择，针对本研究所鉴定的关键基因（如那些参与木质素合成与降解、细胞壁可及性调控的基因）进行定向改良，有望打破这种“耐受性-消化性”的权衡，培育出既抗逆高产、又优质易消化的新一代牧草品种，对保障畜牧业可持续发展和生物能源开发具有重要的应用价值。
