编辑推荐：
本研究通过大规模野外采样和免疫反应生物测定，揭示了气候（尤其是降水量）是影响蜻蜓目宿主黑化封装免疫反应空间变异的关键驱动因素，而宿主特征和寄生虫压力作用不显著。这为宏观免疫学提供了新见解，并强调了在流行病学模型中纳入环境因素的重要性。
引言
免疫系统是生物体抵御寄生虫和病原体的首要防线。然而，对于驱动个体、种群和物种间免疫力差异的因素，我们知之甚少，这些因素也鲜少被纳入流行病学模型。宿主的特征、寄生虫暴露以及环境因素都能影响免疫力，但这些因素如何塑造免疫反应强度的空间变异仍未被充分探索。本研究旨在评估非生物环境、寄生虫压力以及宿主特征对宿主免疫力的相对影响。
材料与方法
研究系统
本研究以蜻蜓和豆娘（统称蜻蜓目昆虫）为宿主，它们是其体表寄生虫——Arrenurus属水螨的常见宿主。水螨寄生受宿主体重、水体pH值、电导率、溶解氧和气温等多种因素影响。蜻蜓目昆虫通过无脊椎动物共有的一种通用免疫防御机制——酚氧化酶级联反应来抵御水螨。该反应的终产物黑素能够包裹水螨的取食管，防止其吸食宿主的血淋巴。
采样与免疫反应测定
研究团队于2019年和2020年夏季，在加拿大东部横跨温带和北方森林生物群系的42个地点，机会性地采集了来自63个物种的1237只蜻蜓目个体。为量化免疫反应强度，研究人员将一段灭菌的尼龙单丝植入宿主胸部腹侧，模拟寄生虫的取食管。12小时后取出植入物并拍照。通过ImageJ软件将图像转换为8位灰度图，测量植入部分与未植入对照部分的平均灰度值。免疫反应指数（IRi）计算公式为：IRi= CAi- SAi，其中CAi为对照部分平均灰度值，SAi为植入部分平均灰度值。灰度值范围为0（纯黑）至255（纯白），因此IRi值越高，表明黑化封装反应越强。
研究者还记录了每个宿主的体重（干燥后称重）以及其体表附着的水螨数量。对于每个采样点，从WorldClim数据库提取了年平均温度和年平均降水量作为气候变量，并使用pH计测量了湖泊的pH值。
统计分析
所有统计分析均在R语言环境中进行。为了探究免疫反应与局部环境及寄生虫的关系，研究使用了广义加性模型（GAMs）。该模型将免疫反应指数作为响应变量，将pH值、降水量（线性项）、温度（非线性平滑项）、种群水平平均螨虫丰度以及个体水平偏离平均的螨虫负荷作为预测变量。同时，模型包含了宿主科、宿主物种、采样点以及嵌套的水域等级（5级和6级流域）作为随机效应，以解释空间和物种层面的依赖性。
此外，研究还利用结构方程模型（SEMs）来剖析局部环境、寄生虫和免疫力之间复杂的互联关系网络。该模型假设局部气候（降水、温度）影响局部水环境（pH值）和宿主体型，进而影响寄生虫感染（种群水平和个体水平），最终共同作用于黑化免疫反应。
结果
空间与系统发育信号
研究发现，所有协变量（包括黑化反应、温度、降水、pH值、体重等）均存在显著的空间自相关（所有变量的Moran's I检验p值 < 0.001）。其中，黑化反应在小于约2.2个十进制度（约270公里）的距离内呈中度正空间自相关。除体重显示出较强的系统发育保守性（ΔAIC支持系统发育模型）外，其他性状（如黑化反应、螨虫丰度等）均未检测到明显的系统发育信号，随机效应模型的支持度更高。这表明这些性状的变异主要受环境等因素驱动，而非系统发育关系的约束。
局部气候是免疫防御的关键预测因子
在对全部数据集的分析中，广义加性模型显示，黑化反应强度随采样点年平均降水量的增加而显著下降（参数回归系数 = -0.21，p < 0.0001）。然而，黑化反应与年平均温度、站点pH值、种群水平平均螨虫丰度或个体水平螨虫负荷的偏离均无显著关系。该GAM模型解释了数据中80.1%的偏差，调整后R2= 0.79。
免疫力的环境驱动因素
结构方程模型的结果进一步证实了气候的主导作用，并提供了更细致的图景。模型显示，黑化反应与降水量呈强负相关，与温度呈中等负相关，与宿主体重和种群平均螨虫丰度呈弱正相关，但与个体螨虫负荷偏离呈弱负相关。此外，pH值随降水量增加而中度下降。模型未检测到气候、pH值、螨虫寄生和免疫反应之间存在显著的间接关系网络。
讨论
本研究的结果表明，在解释蜻蜓目宿主黑化免疫反应的地理变异方面，非生物环境（尤其是气候）的影响超越了寄生虫压力和宿主特征。这与最初的假设——即局部环境和寄生虫压力共同驱动免疫防御变异——形成了对比。
降水量与免疫反应的关联机制
降水量是解释宿主个体和种群间免疫反应变异的最佳因子，在最为干旱的地区观察到了最强的免疫反应。一种可能的解释是，黑素在昆虫中除用于寄生虫防御外，还参与体温调节（通过体色变深）和保水（通过降低角质层渗透性）。在干旱环境中，增强的角质层黑素含量有助于抵抗干燥胁迫，而基于黑素的封装免疫反应可能是这种适应性进化的副产品。这一发现支持了“黑化-干燥抵抗”假说。另一种可能性是，紫外线B辐射（通常与降水量/云量呈负相关）可能驱动了观察到的模式，因为紫外线B已知会影响其他无脊椎动物的免疫防御。
寄生虫作用的复杂性
尽管在所有个体的总体分析中未发现寄生程度与黑化反应间的整体关系，但基于部分物种数据的结构方程模型提示，可能存在宿主对本地水螨寄生虫的局部适应。来自寄生虫丰度较高种群的宿主表现出更强的黑化反应，这可能是一种适应性反应。然而，重度感染本身可能消耗或抑制黑素储备，从而导致个体水平上感染偏离与免疫反应呈负相关。这种负相关关系符合黑化反应的机制逻辑：为了封堵更多的螨虫取食管，宿主需要产生更多黑素。
系统发育约束的缺失
分析显示，除体重外，螨虫丰度和免疫防御等关键性状均未表现出系统发育信号。这强化了我们的核心发现，即这些性状主要受环境调节，而非受系统发育限制。这意味着在分析蜻蜓目免疫防御时，考虑物种效应比考虑系统发育关系更为重要。
研究意义与局限
本研究基于丰富的野外数据（包括超过一千次免疫测定），强调了局部环境在调节宿主-寄生虫互作中的核心作用。将气候和环境因素如何影响免疫防御纳入分析，揭示了宿主-寄生虫互作中固有的细微差别。这一发现对于改进流行病学模型、以及理解气候变化（特别是温度和降水模式改变）情景下宿主-寄生虫关系的变化具有重要意义。
需要指出的是，本研究基于相关性数据，结构方程模型虽能检验假设的“因果”路径，但无法推断因果关系。此外，使用的寄生虫度量较为笼统（未鉴定水螨至物种水平），并且部分物种的代表性优于其他物种。未来更均衡的采样和更深入的实验研究将有助于验证和细化这些模式。
总之，本研究通过宏观生态学视角，揭示了气候条件对无脊椎动物免疫反应具有超越宿主特征和寄生虫压力的主导性影响，为宏观免疫学领域提供了关键实证，并呼吁在未来的疾病风险预测和管理模型中，更充分地纳入空间显式的环境驱动因素。
