《The Plant Genome》：Leaf transcriptome differences between diploid and tetraploid bahiagrass
编辑推荐：
本研究通过对亚热带牧草巴哈雀稗（Paspalum notatum）的二倍体与四倍体性系植株进行叶片转录组比较，系统揭示了全基因组复制（WGD）引起的广泛转录组重塑。研究发现，多倍体化影响了发育、氧化还原稳态、光合作用等关键通路，鉴定出大量差异表达转录本（DET），包括转录因子（TF）和激素相关基因等关键育种靶点。共表达网络分析进一步鉴定出532个受基因组加倍影响的主调控基因（switch genes），并发现其在基因组特定区域（如5号、10号染色体）存在非随机聚集的“热点”。该研究为理解多倍体化的分子后果提供了新视角，所建立的转录本目录及鉴定的候选基因，为旨在提高胁迫耐受性、活力等农艺性状的分子育种项目提供了宝贵资源。
引言：多倍体与巴哈雀稗研究模型
多倍体，即全基因组复制（WGD）导致的多套染色体组现象，是植物进化中的常见事件，在增强遗传多样性和环境适应性方面扮演关键角色。亚热带牧草巴哈雀稗是研究多倍体、生殖模式与胁迫关联的理想模型。自然界中，该物种的二倍体（2n = 2x = 20）基因型通过自交不亲和的性生殖繁殖，而主要的四倍体（2n = 4x = 40）基因型则通过假配子无孢子生殖（apomixis）进行无性繁殖。尽管此前有研究关注花器官中与生殖相关的基因表达变化，但叶片对多倍体化的转录组响应尚未被专门探究。本研究旨在通过比较性系二倍体与四倍体巴哈雀稗的叶片转录组，排除无孢子生殖控制位点（ACL）的干扰，明确鉴定出与倍性水平（如胁迫耐受性、巨型效应）特异性相关的候选基因，为基于分子的育种策略提供支持。
材料与方法：转录组构建与分析流程
研究使用了来自NCBI数据库的转录组测序原始数据，涉及两个性系二倍体（2xsex）和一个性系四倍体（4xsex）基因型，每个基因型包含三个生物学重复。首先，利用已有的巴哈雀稗R1参考基因组进行参考转录本组装，构建了二倍体和四倍体叶片转录组。通过主成分分析（PCA）和聚类分析验证了样本按倍性水平清晰分离。差异表达（DE）分析基于先前构建的参考叶片转录本组，使用Kallisto进行表达量定量，并采用DESeq2和limma-voom流程（后者结合duplicateCorrelation功能以考虑基因型内相关性）进行差异表达鉴定，设定错误发现率（FDR）< 0.05且|log2FC| > 2为标准。对鉴定出的差异表达转录本（DET）进行了基因本体（GO）富集分析、KEGG通路映射、转录因子（TF）和激素相关基因分类。此外，利用SWItch Miner（SWIM）工具构建基因共表达网络，识别受基因组加倍影响的主调控基因（switch genes），并通过泊松富集检验分析其在基因组上的分布。
结果
1. 样本聚类与转录组组装
PCA和聚类分析显示，所有生物学重复紧密聚集，并清晰地按倍性（二倍体 vs. 四倍体）分组。组装得到的二倍体叶片转录本组包含76,681个转录本，四倍体转录本组包含69,996个转录本。BUSCO评估显示转录本组完整性良好。
2. 二倍体与四倍体叶片中的差异表达转录本
初步DESeq2分析在二倍体与四倍体样本间鉴定出40,844个具有显著差异表达的转录本（DET），其中18,785个在二倍体中过表达，22,059个在四倍体中过表达。为验证倍性效应并减少假阳性，后续采用包含更多基因型的分析以及limma-voom流程进行验证，最终获得28,054个重叠的“真实DET”用于后续分析。
3. 受倍性增加影响的分子通路
GO富集分析显示，在细胞组分层面，最富集的类别包括质体膜、外膜、叶绿体类囊体等，表明多倍体化对光合装置有显著影响。在生物过程层面，受影响最大的通路包括蛋白质运输、非编码RNA加工、单羧酸代谢过程和叶片发育等。分子功能类别则富集了蛋白质复合物结合、ATP水解酶活性等。KEGG通路富集分析表明，卡尔文循环（M00165）、糖异生（M00003）、单木质醇生物合成（M00039）等中心碳水化合物代谢和木质素合成通路在多倍体植物中被显著重组。MapMan分析进一步证实，全基因组复制（WGD）影响了蛋白质稳态、RNA生物合成、染色质组织、脂质代谢和光合作用等主要类别。与胁迫响应和信号传导相关的类别，如植物激素作用、氧化还原稳态和DNA损伤反应也显著富集。特别是在氧化还原稳态类别中，多个活性氧清除酶（如超氧化物歧化酶、过氧化氢酶、谷胱甘肽过氧化物酶等）的表达发生失调。
4. 二倍体与多倍体植物中差异表达的转录因子和激素相关转录本
转录因子分析显示，受WGD影响最显著的家族包括bHLH、bZIP、C2H2、MYB-related、NAC和WRKY等。其中，一些TF家族（如B3、Dof、TCP等）在四倍体中持续上调，而另一些（如DBB、YABBY）则在二倍体中过表达。激素相关转录本分析表明，脱落酸（ABA）、生长素（Auxin）和茉莉酸（JA）是受影响最显著的类别。热图清晰展示了与赤霉素、茉莉酸、脱落酸、乙烯、细胞分裂素和生长素相关的多个转录本群的差异表达模式。
5. 共表达网络分析
SWIM相关性网络分析从2835个节点中鉴定出532个主调控基因（switch genes）。网络稳健性通过留一法（LOO）和80%子抽样分析得到证实。在这些主调控基因中，有17个编码转录因子（如bHLH、AP2、NAC等家族成员），14个编码蛋白激酶，10个与光合装置特异性相关（如PsbP），35个被鉴定为胁迫响应基因（如谷胱甘肽S-转移酶GST）。这支持了多倍体在胁迫耐受和稳态维持中的作用。
6. 巴哈雀稗基因组中的主调控基因热点
将532个主调控基因转录本映射到二倍体参考基因组上，泊松概率分布分析揭示了它们在特定基因组区域存在统计上显著的聚集（热点）。主要热点位于5号染色体（10-15 Mb区域）和10号染色体（35-40 Mb的着丝粒周围区域）。其他非随机热点分布在1-4、7和9号染色体上。这表明倍性变化可能影响包含主要调控基因的特定基因组区域，从而大规模改变表达网络。
讨论与意义
巴哈雀稗是南美天然牧场的主要草种，其天然的无孢子生殖特性已催生了新颖、快速的育种计划。本研究通过比较性系二倍体与四倍体基因型的叶片转录组，提供了一个全面的叶片表达基因目录，并鉴定了在基因组复制后上调或下调的基因。本体论分析揭示了与多倍体表型典型相关的类别富集，如光合作用、细胞周期和胁迫响应。转录因子和激素相关转录本的鉴定为靶向育种提供了关键候选目标。值得注意的是，主调控基因在基因组中并非随机分布，而是聚集在特定热点区域，这些区域可能是协调从2倍体到4倍体转变过程中复杂生物途径的基因组枢纽。
综上所述，本研究的结果深化了对巴哈雀稗基因组复制转录组后果的理解，并突出了一系列可用于未来分子育种或基因组编辑计划的候选基因和调控热点，旨在改善该物种及其他重要农艺物种的胁迫耐受性、生产力和恢复力等理想性状。结合已有的参考转录本组、参考基因组和成熟的转化平台，这些资源将助力设计靶向策略，推动分子辅助育种和遗传工程的发展。
