《Molecular Phylogenetics and Evolution》：Deciphering reticulation and reproductive trait evolution in a chasmophyte lineage: a phylogenomic approach to the evolutionary patterns within Petrocoptis
编辑推荐：
研究人员假设伊比利亚半岛岩生属 Petrocoptis的多样化由反复发生的网状演化（reticulate evolution）与繁殖性状的协同进化共同塑造，其中包括一种与种脐（strophiole）相关的种子适应性综合征，该特征可调节干旱悬崖生境下的水分吸收。
研究人员假设伊比利亚半岛岩生属 Petrocoptis的多样化由反复发生的网状演化（reticulate evolution）与繁殖性状的协同进化共同塑造，其中包括一种与种脐（strophiole）相关的种子适应性综合征，该特征可调节干旱悬崖生境下的水分吸收。研究人员分析了代表该属分布范围的63个种群，采用目标富集技术获得345个核基因座（Angiosperms353）及192个质体区域。系统基因组关系在串联法与多物种溯祖（multispecies coalescent）框架下推断，同时考虑了基因谱系冲突、细胞质–核基因组不一致（cytonuclear discordance）、不完全谱系分选（ILS）及网状演化。种子与花部性状在整合的系统发育比较框架中分析，综合评估演化模式、基因组变异与气候相关性。细胞质–核基因组不一致分析与溯祖模拟排除了ILS作为冲突的唯一驱动因素。网络分析支持至少两次网状演化事件，表明存在祖先杂交。繁殖性状演化表现为反复的协同适应性转变：花部演化显示花萼长度多次增加，并与粉色花冠及较高的基因组杂合度相关；种子形态则呈现种脐毛类型的有序转变及不同的适应性最优值——圆柱形毛对应小种子和比例较大的种脐，棒状毛对应大种子和比例较小的种脐。种脐性状与温度和降水量显著相关，构成一种气候关联的种脐适应性综合征，可能优化种子吸水调控。这些结果凸显了网状基因组演化与繁殖性状适应性重构在 Petrocoptis多样化过程中的共同作用。将系统发育网络与性状演化整合，为理解破碎化悬崖系统中基因流与生态特化如何共同推动谱系多样化提供了统一框架。
Petrocoptis是伊比利亚半岛特有的岩生植物属，生长在垂直石灰岩悬崖等高度破碎化的生境中。这类环境具有土壤稀缺、水分供应不稳定、种群间基因交流受限等特点，是植物多样化与适应性演化的理想研究对象。然而，以往基于少量分子标记的研究显示该属的系统发育关系难以解决，质体与核基因组之间存在显著的不一致，推测其演化历史可能涉及不完全谱系分选（ILS）与祖先杂交，但缺乏全基因组尺度的证据。此外，该属在花部大小、花色、种子大小及种脐形态上存在丰富的变异，这些性状的生态功能与演化模式尚不清楚。因此，研究人员开展了本项研究，旨在结合系统基因组学与比较性状分析，解析 Petrocoptis的网状演化历史及其与繁殖性状适应的关系，成果发表于《Molecular Phylogenetics and Evolution》。
研究人员对覆盖全属分布区的63个种群进行了采样，利用目标富集（target enrichment）技术捕获345个低拷贝核基因座（Angiosperms353）与192个质体区域，进行高通量测序与系统发育重建。在数据分析中，采用串联法与多物种溯祖法（ASTRAL?III、SVDquartets）推断物种树，并通过四谱系一致性（QP?IC）、溯祖模拟及四分体一致性因子（qcCF）检验基因冲突来源。为明确网状演化事件，研究人员使用 PhyloNet 与 NANUQ 算法构建物种网络。种子与花部性状的形态测量在种群水平进行，结合系统发育信号分析、祖先状态重建及适应性模型（Bayou、mvMORPH）评估性状演化模式与气候关联。此外，通过 SNP 基因型数据计算高度杂合位点数量（HH?SNP），并在系统发育混合模型中检验其与花部性状的关联。
研究结果
3.1 序列恢复与比对指标
目标富集成功恢复了345个核基因座与192个质体区域，各分类单元在矩阵中均得到良好代表，缺失数据比例低，保证了系统发育推断的稳定性。
3.2 系统发育关系
核基因组物种树在各推断方法中高度一致，支持将种群划分为五个主要分支，分别对应地理区域的谱系。质体基因组树与核基因组树存在广泛的不一致，例如 P. glaucifolia在质体树中呈现非单系，某些类群的质体谱系与其他地理群共享，提示质体捕获现象。
3.3 细胞质–核基因组及基因谱系冲突分析
溯祖模拟显示，大多数质体谱系在仅考虑 ILS 的情况下出现频率低于22%，表明 ILS 无法单独解释观察到的不一致。qcCF 分析检测到显著偏离多物种溯祖模型的四分体集中在特定节点，支持局部网状演化。θ 参数估计值整体较低，但仍伴随高冲突，进一步排除 ILS 的主导作用。
3.4 系统发育网络推断
PhyloNet 分析支持两个网状事件：一次导致 P. pardoi的起源，涉及 P. montsicciana–P. pseudoviscosa祖先与 P. guarensis祖先的基因交流；另一次涉及 P. glaucifolia内部支系的杂交。NANUQ 网络同样揭示复杂的网状结构。
3.5 连续型种子与花部性状的演化
种子面积、种脐面积、种脐相对大小及花萼长度均表现出显著的系统发育信号（Pagel’s λ = 1，Blomberg’s K 显著）。祖先状态重建显示种子面积与花萼长度在主要分支分化时发生定向变化，而种脐性状的变化更多发生在种群水平。种脐面积与相对大小均与年均最高温正相关，与年降水量负相关，种子面积与气候无显著关联。
3.6 离散性状的演化
种脐毛类型的演化最符合有序双向模型（ORD），祖先为圆柱形毛，随后在中–东比利牛斯支系中出现混合型与棒状毛，并在 P. pseudoviscosa支系中逆转回圆柱形毛。花色演化的最适模型为等速率（ER）模型，祖先状态为白色花，随后在多个支系中独立出现粉色花，并伴随少量逆转。
3.7 种子与花部性状的适应性机制转变
贝叶斯可逆跳跃模型（rjMCMC）检测到种子面积与花萼长度的多次适应性最优值转变。种子大小在某些支系中增大，在其他支系中减小；花萼长度则在所有支系中呈增加趋势。多最优值 Ornstein–Uhlenbeck（OUM）模型显示，种脐毛类型对应不同的种子–种脐形态最佳组合，花色与花萼长度也存在显著关联。
3.8 遗传多样性与花部性状的关联
HH?SNP 计数在不同谱系间差异显著，且与花萼长度和花色概率正相关。高基因组杂合度的谱系倾向于具长花萼与粉色花，低杂合度谱系多为短花萼与白色花。系统发育方差解释了大部分表型变异。
讨论与结论
研究人员指出，Petrocoptis的演化历史由快速辐射与网状基因流共同塑造，祖先杂交与 ILS 造成广泛的细胞质–核基因组不一致。种脐性状构成了一种与气候关联的适应性综合征，主要通过调节水分吸收来应对悬崖生境的干旱压力，而非单纯反映种子大小的演化。花部性状在系统发育上较为松散，多次独立演化出粉色花与较大花部，可能与传粉者选择有关。高基因组杂合度与显眼花部表型的正相关，暗示更高的基因流与异交水平有助于维持遗传多样性并促进适应性演化。研究强调了在破碎化生境中，网状演化、生态选择与微演化适应共同推动植物多样化的复杂过程。
结论部分强调，整合系统发育网络与性状演化分析为理解岩生植物的多样化提供了新视角，并建议在未来的研究中结合群体基因组学、传粉生态学与受控实验，进一步验证种脐与花部性状的功能及其遗传基础。
