《Ecology and Evolution》：Feeding Strategy Shapes the Evolution of the Hyaluronidase Gene Family in Medicinal Leeches
编辑推荐：
透明质酸酶（hyaluronidase）在水蛭吸血过程中作为“扩散因子”发挥关键作用，然而该基因家族的分子进化仍知之甚少。在本研究中，研究人员结合了系统发育重建、选择压力分析和基于转录组的表达谱分析，调查了八个医用水蛭物种（代表吸血和非吸血谱系）中透明质酸酶基
透明质酸酶（hyaluronidase）在水蛭吸血过程中作为“扩散因子”发挥关键作用，然而该基因家族的分子进化仍知之甚少。在本研究中，研究人员结合了系统发育重建、选择压力分析和基于转录组的表达谱分析，调查了八个医用水蛭物种（代表吸血和非吸血谱系）中透明质酸酶基因家族的进化特征。研究结果表明，所有吸血物种均保留了五个完整的透明质酸酶基因家族成员（hya1–hya5），而非吸血物种缺乏hya5。系统发育分析显示，每个基因家族成员形成一个支持度良好的单系分支，表明祖先基因复制发生在物种分化之前。尽管所有透明质酸酶基因均处于强纯化选择（ω?**论文解读文章**
**一、研究背景与问题**
水蛭（leeches）作为古老环节动物，在东西方医学中均有悠久应用历史，其药用价值主要源于唾液中的多种生物活性成分，如抗凝血剂hirudin、弹性蛋白酶抑制剂Eglin C、凝血因子Xa抑制剂antistasin及血小板受体拮抗剂decorsin等。然而，现有研究多聚焦于直接抗血栓物质，对另一类关键分子——透明质酸酶（hyaluronidase）关注不足。透明质酸酶通过降解宿主细胞外基质（extracellular matrix, ECM）中的透明质酸，作为“扩散因子”加速抗凝血等生物活性分子的渗透，对吸血效率至关重要。尽管已有研究报道了部分物种的透明质酸酶三维结构，但其基因家族在不同水蛭谱系间的分子进化特征、与摄食策略的关联仍不清楚。因此，研究人员选取了8个代表不同摄食策略（吸血与非吸血）的医用水蛭物种，旨在阐明透明质酸酶基因家族的进化动态及其影响因素。该研究发表在《Ecology and Evolution》。
**二、主要技术方法**
研究人员整合了比较基因组学与转录组学手段。首先，对野外采集的6个物种（Hirudinaria manillensis、Hirudo nipponia、Hirudo tianjinensis、Whitmania acranulata、Whitmania laevis、Whitmania pigra，分别来自海南定安和天津宝坻）进行头部组织转录组测序（RNA-Seq），并补充公共数据库中Hirudo medicinalis和Hirudinaria bpling的基因组与转录组数据。通过BLAST同源搜索、Exonerate验证外显子-内含子结构，鉴定出具有转录支持的透明质酸酶基因。随后，利用最大似然法（IQ-TREE）构建系统发育树（以Haemadipsa yanyuanensis为外群）；采用PAML的two-ratio分支模型分析选择压力（ω=dN/dS）；使用Salmon量化表达水平（TPM），并通过非参数检验和系统发育方差分析（phylogenetic ANOVA）比较吸血与非吸血组间的表达差异。
**三、研究结果**
**3.1 序列特征与结构比较**
通过多序列比对和基因结构分析，发现所有8个物种中共鉴定出5个透明质酸酶基因家族成员（hya1–hya5）。吸血物种（Hirudinaria和Hirudo属）均保留全部5个成员，而非吸血物种（Whitmania属）缺失hya5。各基因的CDS长度（1455–1566 bp）和编码蛋白长度（485–522 aa）高度保守，且hya1–hya3含10个外显子，hya4–hya5含11个外显子，外显子结构在物种间一致。序列一致性分析表明，近缘物种间序列一致性超过95%，远缘物种间约72%–90%。
**3.2 系统发育关系**
基于CDS构建的系统发育树（图1）显示，每个基因家族成员（hya1–hya5）均形成支持度良好的单系分支（bootstrap值>75%），表明祖先基因复制发生在现存物种分化之前。吸血物种在所有5个分支中均有对应拷贝，而非吸血物种在hya5分支中缺失代表序列。此外，H. tianjinensis的hya4基因检测到提前终止密码子，可能为假基因。
**3.3 选择压力分析**
采用two-ratio分支模型（吸血谱系为前景支，非吸血谱系为背景支）进行选择压力分析。整体上，吸血谱系的平均ω值（0.198）显著高于背景支（0.119）（p<0.001），但所有ω均远小于1，表明纯化选择是主要进化力量。基因特异性分析显示，hya2、hya3、hya4在吸血谱系中ω显著更高（p<0.05），而hya1无显著差异。这表明吸血谱系中透明质酸酶基因经历了相对放松的选择约束，但未检测到正选择信号。
**3.4 不同摄食策略水蛭中透明质酸酶基因的表达比较**
转录组表达分析（TPM）表明，吸血水蛭透明质酸酶基因的总表达量显著高于非吸血水蛭（Mann-Whitney U检验，p<0.0001；系统发育ANOVA，p=0.015）。其中，hya1、hya2和hya5在吸血谱系中表达显著更高（p<0.05）。层级聚类和主成分分析（PCA）清晰地将吸血与非吸血物种分开，吸血物种聚为一类且表达水平较高，非吸血物种聚为另一类且表达较低。
**3.5 表达差异与调控模式**
进一步分析显示，透明质酸酶基因家族成员间表达存在显著差异（Kruskal-Wallis检验，p<0.0001），hya1和hya2表达水平在所有成员中最高，且在吸血物种中显著高于非吸血物种。种内比较表明，除W. pigra外，其余物种内各成员表达差异显著（Friedman检验，p<0.001）。吸血物种维持多个成员的高表达，而非吸血物种缺少hya5，且整体表达较低。
**四、讨论与结论总结**
讨论部分指出，透明质酸酶基因家族的进化与摄食策略密切相关。首先，祖先基因复制后，吸血物种通过长期保留更多祖先拷贝（剂量效应和亚功能分化）增强ECM降解能力，而非吸血物种因摄食改变而丢失部分成员。其次，尽管所有基因受强纯化选择以维持核心催化功能，但吸血谱系中ω值升高表明选择约束相对放松，可能为功能适应性提供进化灵活性。最后，表达分析显示吸血物种中透明质酸酶基因高表达（尤其是hya1和hya2）是适应吸血行为的关键转录调控策略；而H. tianjinensis中hya4的低表达提示近缘物种间存在谱系特异性表达分歧。研究也指出存在局限性，如样本量有限、有效种群大小未知、地理来源单一等，可能影响推断强度。
**研究结论翻译：**
“本研究系统分析了具有不同摄食策略的水蛭中透明质酸酶基因家族，检查了基因拷贝数、选择压力和表达模式，从而揭示了与吸血相关的多个进化模式。研究结果表明，每个透明质酸酶基因家族成员形成一个支持度良好的单系群，这些基因在物种间高度保守，表明这些基因在物种分化前已发生分化。吸血水蛭并未通过近期基因复制获得额外的透明质酸酶拷贝；相反，祖先多拷贝基因的长期保留可能有助于维持与吸血相关的ECM降解过程。在分子进化水平上，透明质酸酶基因家族主要受纯化选择（ω<1）塑造，反映了其关键酶学功能的强功能约束。然而，吸血谱系表现出相对升高的ω值，表明选择约束适度放松而非强正选择，意味着吸血物种中有限但可能具有功能的进化灵活性。表达分析进一步证明，透明质酸酶基因在吸血水蛭中一致地更高表达，其中hya1和hya2表现出特别高的表达水平，表明转录调控可能在透明质酸酶功能与吸血策略的关联中起重要作用。值得注意的是，H. tianjinensis中hya4表达的变化表明，在近缘物种中，单个基因拷贝可能表现出谱系特异性表达差异。总体而言，透明质酸酶基因家族的进化受本研究所检查的Arhynchobdellida目水蛭物种中基因保留、序列水平约束和表达调控的共同影响。在包括Rhynchobdellida目在内的其他水蛭谱系中是否出现类似的基因保留、选择约束和表达调控模式，有待未来研究检验。”
