《Biology》：Comparative Genomics of NGF, NTRK1, and NGFR Across Chordates: Joint Gene Recovery and Region-Specific Constraint
编辑推荐：
神经生长因子（Nerve Growth Factor, NGF）是一种对感觉神经元生物学和疼痛信号传导至关重要的神经营养因子。NGF与NTRK1分别是有颌类脊椎动物（gnathostome）中古老的神经营养因子家族和Trk受体家族的旁系同源基因（paralog
神经生长因子（Nerve Growth Factor, NGF）是一种对感觉神经元生物学和疼痛信号传导至关重要的神经营养因子。NGF与NTRK1分别是有颌类脊椎动物（gnathostome）中古老的神经营养因子家族和Trk受体家族的旁系同源基因（paralogs），而NGFR编码p75NTR，属于肿瘤坏死因子受体超家族。由于成熟神经营养因子是富含半胱氨酸的短配体，NTRK1具有高度相似的旁系同源基因，且NGFR位于Trk受体家族之外，因此对脊椎动物中NGF、NTRK1和NGFR的直系同源关系（orthology）评估颇具挑战性。本研究评估了三个基因在物种层面的联合基因组恢复情况，并比较了预设配体和受体区域间的进化约束。分析覆盖94种脊索动物，结合分级恢复证据与家族系统发育、基因座和结构域证据、密码子模型、氨基酸保守性及结构映射。纳入16种非有颌类物种以提供进化背景，不用于计算三个目标基因的恢复率。在78个有颌类样本中，52个物种实现了高置信度的联合恢复，8个物种获得支持的临时联合恢复，4个物种具有与可能共同恢复一致的组件级证据，14个物种尚不完整。成熟NGF的非同义替换率与同义替换率之比（dN/dS, ω）的数值估计低于前导肽区（pro-domain）（0.067对比0.352）。所有六个支架半胱氨酸均不变；在残基类别比较中，半胱氨酸支架是唯一在校正后与坐标解析背景存在显著差异的结构类别。二聚体界面、TrkA接触区和p75NTR接触区类别在校正后与该背景无显著差异。受体分析显示Trk激酶结构域和p75NTR胞内区域的ω估计值较低。综上，这些结果为目标基因的恢复和区域纯化选择（purifying constraint）提供了置信度分级的图谱，同时将非有颌类的同源基因定位于更古老的神经营养因子系统中。
论文解读：神经生长因子（NGF）、NTRK1与NGFR在脊索动物中的比较基因组学
研究背景与科学问题
神经生长因子（Nerve Growth Factor, NGF）是最早被发现的神经营养因子，在感觉神经元存活、神经生物学和疼痛信号传导中发挥核心作用。NGF主要通过由NTRK1基因编码的TrkA受体传递信号，而由NGFR基因编码的p75NTR受体则能够调节NGF与TrkA的结合及下游信号，或在不同的受体环境中介导独特的生物学效应。除了经典的神经功能外，NGF信号通路还参与损伤后的神经向内生长、骨骼修复、骨关节炎重塑、肌腱纤维化以及肝细胞癌中对乐伐替尼（lenvatinib）的耐药反应等病理生理过程。
尽管已有研究重建了神经营养素信号系统的演化历史，但这些工作主要聚焦于基因家族的演化，而非物种水平的恢复置信度。由于成熟神经营养因子是富含半胱氨酸的短配体，NTRK1需与高度相似的旁系同源基因NTRK2和NTRK3区分，而NGFR又独立于Trk受体家族之外，因此对有颌类脊椎动物中NGF、NTRK1和NGFR进行准确的直系同源关系评估极具挑战性。早期脊椎动物基因组复制事件产生了密切相关的旁系同源基因，后续的基因组重排、谱系特异性基因保留或丢失，以及不同物种间基因组组装连续性和注释质量的差异，进一步增加了比较分析的复杂性。简单地将基因标记为“存在”或“缺失”会掩盖分配的不确定性和基因家族演化史。因此，研究人员开展了这项研究，旨在通过分级证据体系，系统评估NGF、NTRK1和NGFR这三个目标基因在广泛脊索动物物种中的联合基因组恢复情况，并比较配体和受体各区域的进化约束模式。
研究方法概述
研究人员构建了一个包含94种脊索动物的分析面板，涵盖头索动物、尾索动物、圆口纲动物、软骨鱼纲动物、辐鳍鱼纲动物、非四足肉鳍鱼纲动物、非哺乳四足动物和哺乳动物。主要技术方法包括：（1）基于NCBI Datasets、UniProt、OrthoDB和Ensembl等多数据库的基因恢复清单构建，采用六种操作性证据状态对每个基因-物种组合进行分级评估；（2）利用MUSCLE和MAFFT进行多序列比对，通过IQ-TREE 2构建最大似然法系统发育树进行基因家族验证；（3）结合反向分类法（reciprocal BLASTP）、profile-HMM筛选、局部共线性（synteny）证据和结构域架构分析，对非四足动物的候选基因进行裁决；（4）使用PAML软件包中的codeml程序进行密码子替换模型分析，估算非同义替换率与同义替换率的比值（dN/dS, ω）；（5）基于PDB结构数据进行氨基酸保守性评分和结构映射分析，采用曼-惠特尼U检验比较不同结构类别的保守性分布。
研究结果
3.1 联合恢复在不同脊索动物类群中存在差异
在78个有颌类样本中，52个物种显示出NGF、NTRK1和NGFR的高置信度联合恢复，8个物种获得支持的临时联合恢复，4个物种具有与可能共同恢复一致的组件级证据，14个物种尚不完整。在软骨鱼纲中，仅两种鲨鱼获得了部分NTRK1样模型，因此被归类为可能的共同恢复而非支持的联合恢复。在辐鳍鱼纲中，20个样本中有9个呈现高置信度联合恢复。整体而言，非四足硬骨鱼类中有77.3%满足支持的联合恢复定义，四足动物中这一比例为86.0%。高置信度联合恢复在四足动物中比例更高，42/50的四足动物达到高置信度标准，而非四足硬骨鱼类中仅为10/22。此外，NTF4基因在辐鳍鱼纲和非哺乳四足动物中的恢复频率显著低于哺乳动物。
3.2 基因家族系统发育将NGF和NTRK1与其旁系同源基因分离
在未生根的神经营养因子系统发育树中，NGF与BDNF、NTF3和NTF4被一条独占边分开，71个精选NGF序列的支持率为100/100（SH-aLRT/UFBoot）。Trk受体树同样将NTRK1与NTRK2和NTRK3分开，64个NTRK1序列形成支持性的独占分支。在独立的NGFR/TNFRSF背景树中，所有65个焦点NGFR标记序列在遗传距离筛选中均最接近人类NGFR/TNFRSF16。
3.3 基因座和结构域证据完善了非四足动物候选基因的分配
对于边界案例，研究人员通过局部共线性和结构域架构分析进行了详细裁决。矛尾鱼的NGF候选基因保留了10个保守的邻近标记基因，小体鲟的NGF候选基因保留了7个，两者均被保留为临时NGF分配。雀鳝的NGFR候选基因包含四个富含半胱氨酸的结构域、一个跨膜区和一个死亡结构域，并带有11个保守的基因座标记。猫鲨的NGFR记录具有完整的p75NTR兼容架构，被高置信度分配。圆口纲动物的NGF样信号被反向分类为NTF3样，缺乏NGF的邻域特征。
3.4 成熟NGF表现出区域特异性约束
NGF前导肽区的平均氨基酸保守性为0.58，而成熟NGF为0.84。独立拟合的密码子模型显示，全编码区的ω值为0.175，前导肽区为0.352，成熟NGF为0.067。在成熟NGF内部，六个支架半胱氨酸完全保守，其保守性得分均为1.000，与42个其他具有解析坐标的成熟残基背景组相比存在显著差异。二聚体界面、TrkA接触区和p75NTR接触区类别在校正后与背景组无显著差异。
3.5 胞内受体区域显示出较低的ω估计值
Trk激酶结构域的ω值极低，NTRK1为0.026，NTRK2为0.014，NTRK3为0.014。在NGFR中，胞质非结构域区域的ω值为0.051，死亡结构域为0.063。相比之下，八个胞外受体分区的中位ω值为0.164。跨膜分区因长度不足30个密码子而未进行估计。
3.6 探索性分析显示分支水平的约束异质性
在两比率分支模型中，所有测试的前景组平均ω估计值均低于1。在27个比较中，14个在校正后仍显著，表明约束强度存在异质性，而非前景均值超过中性。MEME扫描识别出五个未校正p值≤0.1的位点，但在本杰明尼-霍赫伯格校正后无位点保持显著。
讨论与结论
此前的研究已重建了神经营养因子和Trk受体家族的演化，并将它们的扩张和多样化与早期脊椎动物基因组复制联系起来。本研究在此基础上，通过对94种脊索动物应用分级直系同源证据和旁系同源背景，建立了区分高置信度分配与临时及序列水平案例的置信度框架。在有颌类样本中，联合恢复在非四足硬骨鱼类和四足动物中都很常见，但高置信度分配在四足动物中比例更高。将高置信度与临时和序列水平案例区分开来，揭示了二元恢复汇总所无法体现的差异。
进化约束在NGF内部分布不均。成熟NGF的ω估计值低于前导肽区（0.067对比0.352）。不变的半胱氨酸构成胱氨酸结支架，其替代预计会破坏二硫键配对、折叠或分泌。Trk激酶结构域和p75NTR胞内区域的低ω估计值与强纯化选择一致，而胞外区域相对更宽松但并非不受约束。在所有分支和位点分析中，所有前景组平均ω估计值均低于1，检测到的主要进化模式是纯化选择强度的区域变异。
NGF-TrkA-p75NTR信号通路参与疼痛、修复、骨关节炎重塑、肌腱纤维化和癌症治疗反应。不变的半胱氨酸支架和低ω的胞内受体区域突出了定向功能测试的优先靶点，而更宽松的前导肽区和外周成熟NGF位点可能容纳谱系或环境依赖性变异。恢复判断受组装和注释质量影响，功能整合不能仅从联合基因组恢复推断，需要表达和生化证据支持。
研究结论：NGF、NTRK1和NGFR在大多数非四足硬骨鱼类和四足动物样本中被联合恢复，但符合本研究高置信度标准的记录在四足动物样本中更为常见。成熟NGF的估计ω值在数值上低于前导肽区；半胱氨酸支架保持不变，而二聚体界面和受体接触类别的保守性得分在校正后与其他具有解析坐标的成熟结构域残基无显著差异。Trk激酶和p75NTR胞内区域也具有低ω估计值。这些结果将NGF、NTRK1和NGFR的高置信度和支持的临时联合恢复与基于序列水平证据的可能共同恢复区分开来，并突出了三个目标基因中进化约束特别强的区域。
