古老杂交与系统发育冲突：探究菊科（Asteraceae）中的演化复杂性

《Applications in Plant Sciences》：Ancient hybridization and phylogenetic discordance: Exploring evolutionary complexity in Asteraceae

　　

前提：植物系统基因组学中普遍存在的冲突系统发育信号通常反映了由杂交(hybridization)、不完全谱系分选(incomplete lineage sorting, ILS)及全基因组复制(whole-genome duplication, WGD)等过程塑造的演化历史。研究人员旨在识别并评估超多样性菊科（Asteraceae）中的这些复杂过程，以揭示系统发育冲突的潜在成因。 方法：研究人员利用新的及已有的Hyb-Seq和转录组数据分析核质体不一致(cytonuclear incongruence)、WGD及网状化(reticulation)，并提供详细复杂生物信息学分析教程以提高透明度与可重现性。 结果：研究人员发现核基因树间存在广泛冲突及深层网状化事件（尤其见于南美支系）；全家族范围内检测到WGD信号，但因数据完整性差异、事件复杂性及古老起源常难以解读。 讨论：本研究及教程与日益增长的系统基因组研究共同强调了网状演化与WGD在大型多样化支系演化中的作用，同时也突显了相关挑战，预期理论方法的持续进步将进一步推动网状演化实证研究。

论文解读：古老杂交与系统发育冲突——探究菊科（Asteraceae）中的演化复杂性

研究背景与立项依据

随着系统基因组学（phylogenomics）发展，研究者发现基因树冲突（gene tree conflict）与胞核质体不一致（cytonuclear discordance，即核基因组与质体基因组系统发育拓扑结构不同）在植物类群中普遍存在。造成此类系统发育冲突（phylogenetic discordance）的生物学过程主要包括不完全谱系分选（incomplete lineage sorting, ILS）、基因流（gene flow）、古老杂交（ancient hybridization）及全基因组复制（whole-genome duplication, WGD）/多倍化（polyploidy）。菊科（Asteraceae）是被子植物中物种极丰富的超大科（>30,000种），历史上经历过多次大规模基因复制与广泛杂交事件，但其深层演化关系受冲突信号干扰而难以解析。此前Mandel等人（2019）构建了菊科骨架系统树，仍有若干节点支持率低，尤其Anthemideae、Astereae、Senecioneae等族间关系因经历快速辐射与WGD而模糊，且南美早期分化支系采样不足。为此，研究人员整合Hyb-Seq靶序列捕获与转录组数据，系统检测基因树冲突、核—质体不一致、WGD及网状演化（reticulate evolution），并配套开源生物信息学教程发表於《Applications in Plant Sciences》。

主要关键技术方法

研究人员选取菊科47个族（共50族）267个分类操作单元（OTU/taxa），外加姊妹科Calyceraceae 1种及Goodeniaceae 2种作外群；新增Feddeeae和Polymnieae各1个样本的Hyb-Seq数据（Compositae-1061探针组，Arbor Biosciences），新增18个分类群及下载51个已发表转录组（1KP项目及Zhang等）进行RNA-Seq。质体基因组从Hyb-Seq脱靶reads中用GetOrganelle组装并注释，以Bowtie 2比对至各亚科参考质体基因组，用IQ-TREE建最大似然（maximum likelihood, ML）质体系统树。核基因组用PHYLUCE流程保守预测直系同源基因，RAxML构建串联（concatenated）ML核系统树，ASTRAL-III基于单拷贝基因树推断考虑ILS的聚结物种树（coalescent-based species tree），并以调整Robinson–Foulds距离（adjusted Robinson–Foulds distance, RF_adj）比较拓扑差异。采用Quartet Sampling（Pease等）与PhyParts量化基因树—物种树冲突。用wgd v.2.0.38基于转录组旁系同源对Ks分布（intraspecific synonymous substitutions per site, K_S）拟合指数—对数正态混合模型（exponential-lognormal mixture model, ELMM）推断古老WGD时间，并以近缘类群做歧异速率校正。针对骨架关键节点，用PhyloNetworks中SNaQ（Species Networks applying Quartets）方法推断允许不同数量杂交事件（h = 0–10）的系统发育网络（phylogenetic network），并通过自举评估支持率。

研究结果

测序统计（Sequencing statistics）

新测Feddea cubensis与Polymnia canadensis分别获得约328万和476万条Hyb-Seq reads；新增18个转录组样本原始reads介于1444万～2370万条。

质体序列组装与系统重建（Plastid sequence assembly and phylogenetic reconstruction）

从267个分类群脱靶数据中成功组装并注释22个完整环状质体基因组，涵盖Calyceraceae及菊科7个亚科；另从GenBank下载7个质体基因组作参考。质体比较显示Asteraceae具大单拷贝区（large single-copy region, LSC）特征性大片段倒位，Schlechtendalia luzulifolia（Barnadesieae族）具两处大规模倒位区别于其他菊科。基于高覆盖度代表分类群构建的质体ML树显示Calyceraceae（Gamocarpha macrocephala）与Barnadesieae族（Fulcaldea stuessyi）形成姊妹关系，整体质体拓扑与串联核基因组及聚结核基因组树均存在明显不一致，证实胞核质体冲突。

核系统发育分析（Nuclear phylogenetic analyses）

PHYLUCE保留1047个靶位点（原1061），各分类群回收位点数16～584不等。串联ML核树中几乎所有族为单系，Cyclolepis genistoides在本研究中恢复为Gochnatieae的姊妹而非Wunderlichieae。ASTRAL-III聚结物种树与串联树拓扑差异显著（RF_adj= 0.303），部分支系关系不稳定，反映ILS影响。转录组BUSCO完整性均值约66.9%，Famatinanthus decussatus因来自基因组数据仅10.7%完整。

鉴定系统发育冲突（Identifying phylogenetic discordance）

Quartet Sampling显示串联核树中菊科主干若干节点（Barnadesioideae与其余菊科分离节点、Stifftioideae＋Mutisioideae与其他支系分离节点、Gochnatieae—Hecastocleideae—Oldenburgieae＋Tarchonantheae相关连续节点）具低四基因合意值（Quartet Concordance, QC）与高冲突；聚结树中Calyceraceae＋Barnadesieae与其内部相互分离的节点呈高冲突。PhyParts同样显示主干基因树高度不一致，尤在Barnadesioideae之后分支。

推断全基因组复制（Inferring whole-genome duplications）

对69份转录组旁系同源K_S分布拟合显示Asteraceae多数分类群在K_S≈ 0.5、2.25、3.75附近出现峰值；Calyceraceae与Goodeniaceae也出现近似K_S≈ 0.5–0.55峰，Campanulaceae无此峰。混合旁系/直系同源K_S比较表明Calyceraceae与Asteraceae共享K_S峰（≈0.52–0.55），Goodeniaceae直系同源峰偏高（≈0.64），提示Goodeniaceae分歧早于Calyceraceae–Asteraceae共有WGD事件。Barnadesieae、Cardueae、Cichorieae、Gnaphalieae、Senecioneae、Anthemideae、Helenieae、Neurolaeneae、Chaenectideae及Eupatorieae部分分类群在近K_S较小区间显示较近期WGD信号，但因部分转录组不完整致个别推断存疑。

检验网状演化证据（Testing for evidence of reticulate evolution）

对六个重点分类取样子集运行SNaQ，最佳网络模型允许杂交事件数h = 1（南美子集SA）至h = 7（整体骨架Overall）。排除涉及外群的虚假杂交后，多个子网重复捕获Mutisieae、Barnadesieae及Pertyeae参与网状化；其中三个子网显示Mutisieae参与杂交（亲本涉及Nassauvieae与Onoserideae），四个祖先分支获≥3个独立分析支持含杂交事件。Famatinanthus decussatus分化前及Hecastocleis shockleyi分化后均检出杂交信号，暗示古老基因流可能促成独特征状或适应性。

讨论与结论翻译（Conclusions）

研究人员发现，全基因组复制、古老杂交及基因渐渗在菊科早期多样化中发挥重要作用，可能是形态创新、生态适应及谱系分化的催化剂；灭绝事件亦可能造成主干冲突但无法取样验证。扩充菊科基因组资源与更广分类采样将深化对基因流、系统关系及基因组演化的解析，高分辨率基因组数据与系统发育网络（phylogenetic network）工具对阐明古老网状化、渐渗及分化后果至关重要。本研究随文发布基于R Markdown的生物信息学教程，期望促进后续研究纳入网状分析以拓展对菊科杂交事件的认知。尽管演化复杂性早已被认识，基因组数据增长与计算手段进步使我们有能力剖析并更精确记录贡献于"生命之网"（network of life）——而非单纯"生命之树"——的复杂演化历史。

数据可用性：新Hyb-Seq与转录组原始数据存于NCBI SRA（BioProject PRJNA1246973），其余见PRJNA540287与先前BioProjects；完整质体基因组、树文件及比对结果见Zenodo；分析教程与代码发布于GitHub（https://github.com/erika-r-moore/Ellestad_etal_2025_APPS_Hybridizations）。