《Stress Biology》：Integrative bioinformatics analysis of potential convergent evolution in tardigrade stress-response proteins
编辑推荐：
为探究缓步动物抵御极端环境能力的分子机制，研究人员聚焦于SAHS、CAHS、MAHS和Dsup等关键应激反应蛋白家族，通过对四种缓步动物物种的基因结构、同源关系、保守基序及三维结构进行整合分析，揭示了这些蛋白家族在物种间存在显著的拷贝数变异和结构多样性，并鉴定出一个新的Dsup蛋白(H.Henanensis.Chr5.66)，其结构分析显示跨物种保守的α-螺旋核心。该研究深入阐释了缓步动物应激耐受性的分子进化基础，并为相关生物技术应用提供了新见解。
缓步动物，这种仅有毫米大小、外表憨态可掬的微型水生生物，因其拥有“地表最强”的生存能力而闻名遐迩。它们能在几乎完全脱水、接近绝对零度、沸水高温、强辐射甚至太空真空的极端恶劣环境中生存，被誉为自然界的“终极生存大师”。这种令人难以置信的生存韧性，其背后隐藏的分子密码，一直是科学家们孜孜以求探寻的奥秘。目前，研究普遍认为，缓步动物体内一系列独特的应激反应蛋白，如分泌性丰富热溶性蛋白(Secretory Abundant Heat-Soluble protein, SAHS)、胞质丰富热溶性蛋白(Cytoplasmic Abundant Heat-Soluble protein, CAHS)、线粒体丰富热溶性蛋白(Mitochondrial Abundant Heat-Soluble protein, MAHS)以及损伤抑制蛋白(Damage Suppressor protein, Dsup)，构成了其极端耐受性的核心防线。然而，以往的研究大多局限于少数几个模式物种，例如Ramazzottius varieornatus和 Hypsibius exemplaris。这些关键的应激蛋白家族在其他缓步动物谱系中的分类广度、基因拷贝数变异以及结构多样性究竟如何？它们是否共享着某些趋同进化的分子特征，从而共同铸就了缓步动物的“超能力”？为了回答这些问题，研究人员利用日益丰富的缓步动物基因组资源，开展了一项跨物种的整合生物信息学分析。
为了全面解析缓步动物应激反应蛋白的进化与结构特征，研究人员从公开数据库中获取了四种缓步动物（R. varieornatus、H. exemplaris、Paramacrobiotus metropolitanus和 H. henanensis）的高质量、完整注释的基因组及蛋白质组数据。研究采用了整合的生物信息学方法，包括本地BLAST搜索以鉴定同源蛋白序列，利用GSDS 2.0进行基因结构（外显子-内含子）分析，通过OrthoVenn3平台的OrthoMCL算法进行同源簇分析以识别直系同源和旁系同源关系。系统发育分析则通过MAFFT或MUSCLE进行多序列比对，并使用IQ-TREE2基于最大似然法构建分子进化树。保守基序通过MEME Suite进行识别，蛋白质三维结构则利用AlphaFold Server进行预测建模。此外，研究还使用ConDor工具对缓步动物特有紊乱蛋白进行了潜在的进化趋同位点分析。
研究结果
缓步动物应激反应蛋白的基因结构分析
研究人员首先分析了四种主要应激反应蛋白家族（SAHS、CAHS、MAHS、Dsup）在四个缓步动物物种中的基因结构。结果显示，不同家族和物种间的外显子和内含子长度存在显著变异，反映了物种特异性的适应性分化。例如，SAHS家族的外显子长度跨度最大，而MAHS基因则相对较短。在基因拷贝数上，也观察到明显的物种差异：SAHS和CAHS家族在R. varieornatus、H. exemplaris和 P. metropolitanus中具有较多拷贝，而H. henanensis的拷贝数略少。MAHS家族在所有物种中均严格保守，仅有一个拷贝。尤为值得注意的是，Dsup样序列仅在R. varieornatus、H. exemplaris和 H. henanensis中被发现，在P. metropolitanus中则完全缺失，这可能暗示了该物种采用了替代性的DNA损伤保护策略。
SAHS和CAHS的同源簇分析
通过对SAHS和CAHS蛋白进行同源簇分析，研究人员揭示了这些蛋白在物种间共享和特有的进化模式。SAHS蛋白被分为7个簇和16个单例，其中单例蛋白比例较高，提示可能存在谱系特异性的适应。CAHS蛋白被分为13个簇和16个单例。系统发育树显示，H. exemplaris和 H. henanensis的蛋白倾向于聚在一起，而R. varieornatus和 P. metropolitanus的蛋白形成另一支。与基于保守蛋白构建的物种树相比，SAHS和CAHS的基因树在整体上遵循物种拓扑结构，但部分同源簇显示出不一致的聚类模式，这可能与基因复制和谱系特异性进化事件有关。
应激反应蛋白的分子系统发育分析、保守序列基序和3D结构建模
最大似然法系统发育分析进一步支持了同源簇分析的结果，并显示出与物种亲缘关系一致的聚类模式。保守基序分析在SAHS、CAHS和MAHS家族中分别识别出三个、两个和三个高度保守的基序。对Dsup蛋白的多序列比对也揭示了一些跨物种保守的短基序。蛋白质三维结构建模显示，SAHS1蛋白（以R. varieornatus为模型）具有独特的β-折叠区域，而其他大多数应激反应蛋白则表现出高度的内在无序性。这些结构特征的差异可能对应着不同的应激适应机制。
TDPs中的候选趋同位点
利用ConDor工具对缓步动物特有紊乱蛋白进行的趋同进化分析，在SAHS和CAHS家族中分别识别出三个具有统计学显著性的候选趋同位点。例如，SAHS1参考序列中的K165位点显示出向赖氨酸(K)的十次独立替换事件。重要的是，这些候选位点在结构模型中都定位于之前识别出的保守基序区域内，这增加了它们具有功能相关性的生物学合理性，暗示这些位点的重复替换可能是适应极端环境压力的结果。
Dsup蛋白家族的比较结构和分子分析
本研究的一个重要发现是，在H. henanensis中鉴定出一个新的Dsup蛋白成员（H.Henanensis.Chr5.66），从而将已知的Dsup蛋白多样性扩展至三个。尽管这三个Dsup蛋白在序列水平上相似性较低，但结构建模和比对揭示它们共享一个保守的α-螺旋核心区域，其结构叠加的均方根偏差值在3.8-4.8 ?之间，表明存在中等的折叠保守性。进一步的比较分析显示，这三个蛋白在疏水性分布、电荷分布和离子化特性上也高度相似，特别是都具有一个带正电荷的C末端结构域。这些保守的结构和理化特性暗示了Dsup蛋白在保护DNA免受应激损伤方面可能具有普遍的作用机制。生物信息学预测也支持所有三个Dsup蛋白定位于细胞核。
研究结论与讨论
本研究通过整合生物信息学方法，深入揭示了缓步动物关键应激反应蛋白家族的进化与结构蓝图。研究发现，SAHS和CAHS家族在基因拷贝数上存在扩张和物种特异性变异，这可能是通过基因复制适应不同环境压力的结果。相比之下，MAHS家族在所有分析的物种中严格单拷贝，突显了其在线粒体保护中可能具有不可替代的必需功能。尤其值得注意的是在H. henanensis中发现的新Dsup同源蛋白，以及Dsup在P. metropolitanus中的缺失，这提示不同的缓步动物谱系可能进化出了多样化的基因组稳定性维持策略。
尽管序列多样，但Dsup蛋白在结构上保守的α-螺旋核心、疏水区模式以及带正电的C末端，强烈暗示这些特征对其发挥DNA保护功能至关重要。结合预测的核定位信号，进一步支持了Dsup通过直接与染色质相互作用来稳定基因组、抵御辐射等损伤的假说。对SAHS和CAHS蛋白的候选趋同位点分析，为这些蛋白在应对极端环境压力时可能经历了趋同进化提供了计算证据。
这项研究的意义不仅在于增进了对缓步动物这一神奇生物极端耐受性分子基础的理解，更在于为未来的生物技术应用铺平了道路。Dsup蛋白在辐射防护、CAHS和MAHS蛋白在生物分子干燥保存等方面的潜力，使其在生物医药、农业、太空探索等领域具有广阔的应用前景。例如，利用Dsup增强人类细胞的辐射抗性，或利用CAHS蛋白稳定疫苗和酶制剂。当然，本研究的所有功能推断均基于计算分析，需要后续的异源表达、基因敲除、生化与细胞实验来验证这些候选蛋白和位点的具体机制与生理效应。总的来说，该研究为理解生命极限与分子适应提供了重要见解，并将激发更多针对缓步动物独特蛋白质资源的探索与利用。相关研究成果已发表在期刊《Stress Biology》上。
