《Frontiers in Immunology》：Changes in gene regulation are associated with the evolution of resistance to a novel parasite
编辑推荐：
本研究通过对比不同种群三刺棘鱼在明矾佐剂和绦虫蛋白刺激下的转录组应答，揭示了宿主在适应新寄生虫过程中基因调控网络的重塑机制。重点阐明了纤维化（fibrosis）这一古老生理过程的征用（co-option）及其在抗绦虫免疫中的关键作用，强调了基因表达调控和负反馈增强在进化新防御功能中的重要性。
引言
宿主与寄生虫的相互作用普遍存在，是驱动宿主多样化和进化的重要力量。免疫系统常成为寄生虫驱动选择的目标，其快速进化通常被描述为宿主与寄生虫物种之间持续的“军备竞赛”。然而，免疫进化也常由获得新的、不熟悉的寄生虫所驱动。例如，当海洋种群的三刺棘鱼（Gasterosteus aculeatus）在大约12,000年前殖民淡水湖时，它们遇到了淡水限制性绦虫Schistocephalus solidus，并进化出了抗性。本研究旨在探究宿主适应新寄生虫过程中，是涉及现有免疫通路的调整、旧通路的丢失，还是现有通路被征用于新目的。
方法
本研究比较了实验室饲养的、对绦虫易感性不同的三个三刺棘鱼种群（祖先海洋型SAY、相对抗性湖种群RSL、相对耐受湖种群GOS）在暴露于几种免疫刺激物（明矾、绦虫蛋白或对照注射）后的转录组反应。实验鱼被腹腔注射磷酸盐缓冲盐水（PBS，对照）、明矾（一种佐剂）或绦虫蛋白匀浆。在注射后第1、10、42和90天收集头肾组织进行RNA提取和3‘-TagSeq测序。使用准二项广义线性模型（GLM）分析基因表达数据，评估处理、时间、种群及其交互作用对基因表达的影响，并分析与纤维化评分的关联。
结果
差异表达结果摘要
差异表达建模显示，明矾注射对头肾转录组组成有显著影响，1,365个基因（约10.5%）因明矾主效应而差异表达，远超1%的假发现率预期。这些基因中绝大多数（1,347个）仅受明矾主效应影响，表明明矾的影响独立于种群或时间。明矾×时间交互作用影响195个基因。相比之下，绦虫蛋白处理的主效应仅影响75个基因（0.58%），但处理×时间交互作用影响377个基因，表明绦虫蛋白诱导了随时间变化的动态反应。此外，2,888个基因的表达与纤维化评分相关，389个基因显著依赖于纤维化与种群的交互作用。
对明矾和绦虫蛋白注射的共享与独特反应
比较对明矾和绦虫蛋白的反应，虽然仅有12个基因同时对两种处理差异表达，但两种处理的效应大小在多个水平上呈显著正相关，表明存在共享的转录组反应基础。特别值得注意的是，对明矾有持续反应的基因集与对绦虫蛋白有瞬时反应的基因集存在显著重叠。其中8个免疫相关基因在响应明矾时持续上调，而在绦虫蛋白处理的鱼中，这些基因在第一天上调，但到第42天下调，表明存在绦虫蛋白反应特有的补偿性下调。这些基因大多与纤维化相关，并且在RSL种群中组成型表达更高。
在独特的明矾反应中，大量基因与T细胞功能和抗原呈递相关。近一半（46/110）的独特免疫基因与T细胞有关，其中35个基因因明矾主效应而差异表达。这些T细胞相关基因的表达变化与纤维化密切相关，效应大小呈正相关，且许多在RSL种群中组成型高表达。一些T细胞基因还表现出随时间变化的动态调节模式，初期上调，后期被抑制。
在独特的绦虫蛋白反应中，NF-κB通路相关基因富集。13个受绦虫蛋白主效应或交互效应影响的免疫基因与NF-κB功能相关，其中7个基因的表达随时间动态变化，呈现初期激活、后期抑制的模式。此外，一些炎症相关基因的反应具有种群特异性，特别是在GOS鱼中与其他种群（RSL和SAY）的反应存在差异。
细胞类型特异性标记物的变化
利用已鉴定的头肾细胞类型标记物分析表明，差异表达基因的某些模式可能反映了细胞群丰度的变化。响应明矾主效应的基因显著富集抗原呈递细胞（APC）、B细胞、造血细胞（HC）和自然杀伤细胞（NKC）的标记物，其中APC、B细胞和HC标记物基因主要上调，而NKC标记物基因主要下调。绦虫蛋白×时间交互作用基因则只富集APC标记物，且这些基因均下调。与纤维化主效应相关的基因则富集APC、B细胞、HC和中性粒细胞标记物，且大多与纤维化呈正相关。
种群特异性反应
转录组范围的比较显示，除了GOS与RSL种群对绦虫蛋白主效应的反应呈负相关外，其他种群间对明矾和绦虫蛋白的反应总体上保守。轨迹分析揭示了种群间在应对处理时的差异，尤其是在恢复期。在实验的前42天，各种群和处理类型的轨迹大致相似，且均未观察到回到PBS对照状态。但当纳入第90天数据时，发现种群间恢复轨迹不同：RSL鱼在处理后表现出向基线回归的趋势，而SAY和GOS鱼则继续远离对照状态，这与RSL鱼能从纤维化中部分恢复的观察一致。
讨论
寄生虫对宿主施加巨大的选择压力，是宿主多样化和进化新颖性的关键驱动力。本研究通过比较野生宿主种群的转录组，揭示了宿主在适应新寄生虫过程中防御反应的进化机制。重点在于阐明了一种普遍生理反应（纤维化）被征用以防御特定寄生虫（S. solidus）的机制。
结果表明，宿主适应新寄生虫涉及对古老遗传通路的大规模征用，这些通路本身保持完整，但其表达水平和时序在适应了绦虫的淡水种群中发生了改变。这为越来越多的关于野生种群免疫基因表达快速进化以适应寄生虫群落变化的研究增添了新证据。
对明矾的一般纤维化反应由T细胞相关基因的变化介导，并伴有自我调节模式
对独特明矾反应的转录组分析显示出T细胞活性复杂变化的显著特征，这些变化与纤维化表型密切相关。本研究发现明矾引起了多种T细胞类型相关基因的广泛变化，这与哺乳动物中明矾偏向诱导Th2型反应的传统认知有所不同，表明鱼类中明矾与T细胞活性的相互作用可能更为复杂。值得注意的是，这些T细胞基因活性的变化与纤维化的基因表达模型有显著重叠。许多与纤维化相关的T细胞基因在抗性最强的RSL种群中组成型高表达，但在SAY和GOS种群中，其表达与纤维化评分呈正相关，而在RSL中则不相关或呈负相关。这提示在RSL种群中，与纤维化相关的T细胞通路可能被“预激活”或组成型上调，从而实现对寄生虫更快、更有针对性的防御部署，这种从诱导性到组成型表达的转变是抗性进化中的一个关键步骤。
绦虫特异性反应涉及NF-κB和炎症基因的动态调节
绦虫蛋白诱导了强烈的NF-κB通路基因反应。NF-κB是先天免疫的核心信号通路，在许多动物中功能高度保守。对绦虫蛋白的独特反应涉及多个调控NF-κB活性的基因，这些基因的表达呈现动态变化：在感染初期上调以启动免疫反应，随后被下调以防止过度炎症和免疫病理。这种激活后抑制的模式对于平衡免疫效益与成本至关重要。研究发现，一些NF-κB调控基因在RSL种群中组成型高表达，但在SAY和GOS种群中，其表达与纤维化更密切相关。这表明，在适应过程中，关键的免疫调节基因可能从需要诱导的状态转变为“随时待命”的组成型高表达状态，从而实现对威胁的更快速反应。
种群差异揭示了抗性进化的不同策略
转录组范围的比较和轨迹分析揭示了种群间免疫应答策略的差异。尽管早期反应在许多方面保守，但在恢复期出现了分歧。RSL种群（最具抗性）在处理后表现出部分恢复的迹象，而SAY（祖先型/易感）和GOS（耐受）种群则持续偏离基线状态。这种差异与RSL种群能从纤维化中恢复，而其他种群不能的观察相一致。此外，GOS与RSL种群对绦虫蛋白的反应在基因水平上呈负相关，暗示这两个湖泊种群可能进化出了不同的免疫策略（例如，抗性与耐受）。这些种群差异凸显了免疫调节系统在微进化时间尺度（约10,000代）上发生实质性变化的能力。
结论
本研究结果提供了强有力的证据，表明基因调控的变化以及为减轻免疫病理而增加的负反馈，是现有通路被进化征用以防御新寄生虫感染的关键步骤。三刺棘鱼对绦虫S. solidus的抗性进化涉及对古老、保守的通路（如纤维化和NF-κB信号）的重新利用，并通过改变这些通路的调控方式（如从诱导性到组成型表达，以及增强反馈抑制）来实现。这项工作增进了我们对宿主免疫系统如何通过调节现有通路来快速适应新寄生虫威胁的理解，强调了基因表达调控在免疫进化中的核心作用。
