硬骨鱼多基因趋同进化揭示适应性变异的分子机制

《Molecular Biology and Evolution》：Estimates of molecular convergence reveal multiple genes with adaptive variation across teleost fish

本研究为揭示硬骨鱼适应性进化的分子基础，采用无监督方法在143个物种中检测到89个蛋白编码基因家族存在分子趋同现象。通过整合蛋白质特性、跨物种组织基因表达、斑马鱼胚胎单细胞RNA测序及基因扰动实验，发现趋同基因具有更多拷贝数，涉及胚胎发育、组织形态发生、代谢和热应激反应等过程。趋同替代比非趋同替代更具激进性，且多数基因具有多组织表达特征和多效性，为理解硬骨鱼类谱系特异性适应的进化机制提供了重要见解。

　　

在生命演化的长河中，面对相似的环境挑战，不同物种往往会独立演化出相似的表型特征，这种现象被称为趋同进化。从蝙蝠与海豚分别演化出回声定位能力，到食肉植物独立发展出捕虫机制，自然界充满了这类令人惊叹的进化案例。然而，这些宏观表型背后的分子机制是什么？是否在不同谱系中存在着相同的遗传解决方案？这些问题是进化生物学领域的核心议题。

传统研究通常基于预先设定的表型假设来探索分子趋同现象，这种方法可能遗漏许多未被认识的适应性变化。随着基因组数据的爆炸式增长，科学家们现在能够以更系统、更无偏的方式在全基因组范围内搜索分子趋同的信号。硬骨鱼类作为脊椎动物中物种最丰富的类群，其丰富的表型多样性和生境适应性为研究趋同进化提供了理想模型。

在《Molecular Biology and Evolution》上发表的最新研究中，科学家们开展了一项大规模研究，旨在揭示硬骨鱼类中分子趋同的普遍性和功能意义。研究人员收集了143种硬骨鱼的染色体级别基因组数据，覆盖40个目85个科，利用OrthoFinder进行直系同源基因分析，最终获得了9,273个直系同源群组，包含约200万个直系同源基因。研究采用CSUBST（Combinatorial SUBSTitutions）算法——一种由Fukushima和Pollock（2023年）开发的无监督方法，来检测蛋白质编码基因中的趋同替代现象。

CSUBST算法的创新之处在于其能够评估组合替代的相对速率，通过ω_C（dN_C/dS_C）这一指标来量化趋同进化程度。类似于经典的正选择指标ω（dN/dS）>1，ω_C>1表示超出了中性期望的定向趋同。研究人员设定了严格的阈值（ω_C≥ 5且OCN ≥ 5）来鉴定可靠的趋同事件，最终在89个直系同源群组（274个基因）中发现了显著的趋同证据。

为开展研究，团队整合了多组学数据和分析方法：从NCBI和Ensembl数据库获取143种硬骨鱼基因组并进行直系同源分析；应用CSUBST算法检测分子趋同；利用Bgee数据库进行基因表达和组织特异性分析；通过ZFIN数据库获取基因扰动表型数据；结合单细胞RNA测序（斑马鱼胚胎发育图谱）揭示表达动态；并借助蛋白结构映射和基因本体富集分析解析功能意义。

数据收集与直系同源分析

研究人员从NCBI和Ensembl基因组数据库中筛选了具有染色体级别注释基因组的硬骨鱼物种，最终数据集包含143个物种，覆盖40个目（以鲈形目最多，17种）和85个科（以鲑科最多，10种）。使用斑点雀鳝（Lepisosteus oculatus）作为外群，因为它没有经历硬骨鱼特异的全基因组复制事件。通过OrthoFinder分析，获得了40,940个直系同源群组，包含超过400万个蛋白编码直系同源基因。经过过滤后，最终数据集包含9,273个直系同源群组，约200万个基因，聚焦于跨物种保守的基因集。

硬骨鱼中有多少基因发生了趋同替代？

研究发现大多数趋同事件涉及3个分支，少数案例涉及4、5或7个分支。在分类分布上，鲈形目、刺尾鱼目和鳉形目以及慈鲷科、鲈科和四齿鲀科具有最多的趋同直系同源群组。生态特征分析显示，热带海水鱼中趋同直系同源群的占比最高，但未发现趋同与气候、生态系统类型或食物项等宏观生态特征的整体显著关联。

评估氨基酸替代的功能后果

蛋白结构映射显示，大多数趋同位点位于蛋白质表面残基上，可能与蛋白质相互作用有关。BLOSUM62评分比较表明，趋同替代的平均得分显著低于非趋同替代，表明这些变化在进化中较为罕见，更可能改变蛋白质功能。理化性质分析发现，95%的趋同直系同源群（85个中的89个）在至少一种理化性质上发生了显著变化，这些变化被归类为"激进"氨基酸变化。最显著的趋同变化影响了大体积性质、侧链大小和α转角倾向性等属性。

趋同基因参与哪些生物过程？

基因本体富集分析显示，趋同基因显著富集于生物分子代谢、激素或刺激（如热）响应、胚胎发育和组织形态发生等过程。同时发现了与异生物处理、氧气运输和相机型眼视网膜形态发生相关的富集术语。语义相似性分析表明，趋同直系同源群的GO术语语义相似性低于背景基因集，提示它们与更多样化的生物过程相关。

组织特异性分析发现，约三分之一（28/89）的趋同直系同源群是组织特异性的，其中肝脏中数量最多。与背景基因集相比，趋同基因在组织特异性分布上存在显著差异，中间值组织特异性的频率更高。斑马鱼基因扰动数据显示，趋同基因可以影响多个组织的结构，表现出多效性特征。

胚胎发育过程中趋同基因的单细胞表达动态

对斑马鱼胚胎发育单细胞RNA测序数据的分析揭示了趋同基因的三种表达模式：第一类基因在胚胎发育过程中 ubiquitous 表达，涉及几乎所有细胞类型，如adam17a、eef2b和eef1a1|1；第二类基因在整个胚胎发育过程中表达，但限于特定细胞类型，如alcam-a和pax3a；第三类基因主要在胚胎发育后期表达，如slc22a21和hbae5，分别与排泄血管系统和造血细胞相关。

位点特异性功能的潜在案例

研究发现几个值得关注的趋同位点功能案例。在胚胎血红蛋白α5蛋白（hbae5）中，趋同和分歧替代均位于原卟啉IX含铁配体结合位点附近，这些变化可能与缺氧耐受性适应相关。在溶质载体蛋白slc22a21中，趋同替代位于桶状结构内部通道，靠近1-甲基-4-苯基吡啶鎓配体结合位点，可能影响有机离子稳态。N-乙酰神经氨酸合成酶（nansa）中的趋同替代位于蛋白表面和内部功能配体结合位点，与骨骼发育特别是头部形成密切相关。ADP-核糖基精氨酸水解酶（adprhl1）的趋同替代位于腺苷-5'-二磷酸配体位点，可能与蛋白质翻译后修饰有关。

多效性与适应性进化

研究表明，许多趋同基因具有多功能性和多组织表达特征，可能在生态转变时期产生有益的多效性效应。与传统观点认为多效性会限制适应不同，中等程度的多效性可能通过增加有益等位基因的固定而增强适应潜力。研究发现趋同直系同源群的基因拷贝数分布高于随机选择的非趋同直系同源群，这种模式在86%的目级和83%的科级分类群中一致存在。

基因复制与旁系同源基因表达模式的差异可能通过分割遗传变异效应，帮助克服遗传约束，使旁系同源基因间的变异以隐性分歧形式持续存在。当环境发生变化时，这些隐性变异可能显现，产生适应性表型。

研究结论与意义

这项研究系统鉴定了硬骨鱼类中的分子趋同事件，揭示了趋同进化在硬骨鱼类适应性进化中的普遍性和功能重要性。研究发现趋同基因通常具有多效性特征，在多种组织和发育阶段表达，参与多样的生物学过程。这些基因中的趋同替代往往是激进的氨基酸变化，可能对蛋白质功能产生显著影响。

研究结果支持了多效性基因在适应性进化中的重要作用，挑战了多效性会限制适应的传统观点。相反，中等程度的多效性可能通过同时影响多个性状而促进适应，特别是在生态转变时期。

该研究为理解重复进化提供了新视角，表明分子趋同可能是适应性进化的重要机制。未来研究可以进一步探索这些趋同基因在特定适应性性状中的作用，以及它们在不同生态背景下如何贡献于表型进化。随着更多鱼类基因组和功能数据的积累，我们将能更全面理解分子趋同在生物多样性形成中的贡献。