长期疟疾暴露如何塑造健康成人的基线免疫状态？

肯尼亚四地构成国内对照

这篇论文要回答的问题是：长期居住在高疟疾传播环境中，是否会在当前没有患病的成人身上留下可测量的基线免疫印记。为减少跨国比较（例如非洲人与欧洲人比较）带来的混杂，作者把全部60名参与者都招募在肯尼亚境内，来自按恶性疟原虫传播强度选出的四个研究点：Ngerenya（低至可忽略传播）、Nairobi（极低至无传播）、Junju（中等传播）和Ahero（高传播）。每个研究点15名成人，均临床健康，且在当地居住至少10年。队列以男性为主（68%），中位年龄26岁。

暴露并不是由单一检测定义的。作者把地区传播史、抗恶性疟原虫IgG应答，以及对另外八种常见病毒和寄生虫病原体的血清学谱结合起来，这些病原体包括EB病毒、巨细胞病毒、HSV-1、腺病毒、蛔虫、曼氏血吸虫、粪类圆线虫和溶组织内阿米巴。这一点关系到如何理解后面每一项结果：比较的是低与高累积病原体暴露，疟疾传播是组织这一梯度的主线，而不是被单独分离出来的变量。

抗体谱区分出不同暴露环境

研究点层面的热图显示，与Nairobi和Ngerenya相比，疟疾传播较高的Junju和Ahero在多种病原体特异性IgG应答上普遍升高。对IgG谱做主成分分析后，低暴露与高暴露研究点出现部分分组，两类之间仍有重叠。恶性疟原虫特异性IgG在不同地区存在差异，其中Ahero与Nairobi的对比最大（FDR = 4.37 × 10⁻⁷）。

这种重叠本身就是结果的一部分，而不是可以跳过的缺陷。随后参与者被分为低暴露组（n = 30）和高暴露组（n = 30），血清学显示这一分组是累积暴露的梯度，而不是两个完全分离的人群。作者还指出，队列中相对较低的蠕虫血清阳性率与肯尼亚在持续大规模给药等控制措施后蠕虫传播下降的近期流行病学数据一致，因此所测病原体组合并未覆盖研究点之间所有可能的暴露差异。

658个基因随暴露升高而改变

核心分子测量是对外周血单个核细胞（PBMC，即血液中的混合白细胞组分）做批量RNA测序，并结合65重血浆免疫介质检测。过滤后，14,576个基因进入edgeR差异表达分析，采用Benjamini-Hochberg错误发现率低于0.05的阈值，并把文库制备批次作为协变量。主模型比较高暴露与低暴露，识别出658个差异表达基因。

年龄和性别没有进入主模型：性别分布在暴露组间没有差异，年龄被视为累积暴露的替代指标。这是一个有代价的主动选择。加入年龄、性别或两者同时加入的敏感性分析减少了显著基因的数量，但大体保留了与暴露相关的核心转录特征，因此作者仍以未调整模型为主分析。读者应把658这个数字理解为未调整信号的规模，而不是在所有调整下都固定不变的效应量。

适应性免疫活化伴随炎症程序下降

变化方向是这篇论文最有辨识度的结果。较高暴露与B细胞活化、抗原呈递、免疫调节和转录后RNA加工相关通路的表达升高相关。与此同时，炎症、趋化和血管基因程序降低。作者把这种状态描述为适应性免疫活化加上受抑制的炎症转录，而不是整体性的免疫激活。

加权基因共表达网络分析（WGCNA）把基因归入模块，并检验模块与暴露状态的关联。三个模块说明了这一模式：GM17以B细胞活化和受体信号为中心，围绕CD19，并包含B细胞枢纽基因BANK1和MS4A1；GM20围绕HNRNPA1，涵盖RNA加工和剪接机制；GM30富集以IL1B和IL6为中心的炎症介质。模块结构是论文用来证明这些是协同程序而非零散基因变化的证据，但它只是与暴露组的统计关联，并未证明暴露导致了这种改变。

循环介质升高而炎症转录降低

读者可能会自然地预期：炎症基因表达降低，循环炎症蛋白也应降低。论文报告的却是相反的模式：尽管炎症基因程序降低，循环免疫介质却轻度升高。作者称之为全身介质水平与细胞转录之间的部分分离。

这是一个值得单独保留的测量差异。基因表达读数来自PBMC转录本，而介质读数来自用Luminex平台测得的血浆蛋白浓度（pg/mL）。它们回答的是关于不同组分、不同分子层次的不同问题，因此不能互相替代。论文的解释是：循环细胞中炎症转录基线降低，可以与血浆介质轻度升高同时存在，这也是作者把该状态称为重新校准而非简单抑制的原因。

整合评分区分两个暴露组

为检验各项独立读数是否构成一个连贯状态，作者把选定的基因模块和计算推断的免疫细胞群体标准化为Z分数，再用LASSO逻辑回归整合，以低暴露对高暴露作为二分类结局。得到的整合免疫状态评分在两组间存在差异，双侧Wilcoxon秩和检验P值为1.08 × 10⁻¹¹。

这个数字旁边需要放两点提醒。第一，评分是在同一批60名参与者（每组30人）上构建并评估的，因此它展示的是本队列内的分离，而不是在新人群中验证过的预测性能。第二，免疫细胞比例是用反卷积方法从批量RNA测序中计算推断的（ENIGMA配合LM22特征矩阵，并用独立的ABIS方法重复），并非用流式细胞术实测。这个评分最合理的解读是：转录、细胞和介质特征会一起变化，论文用整合特征之间的Spearman相关矩阵支持这一点。

疟疾处于更广泛的流行暴露之中

研究设计围绕疟疾传播构建，但作者明确指出，这一梯度也反映共流行病原体暴露的更广泛生态差异。由于疟疾传播是研究点选择的基础，也是最显著的暴露梯度，作者认为疟疾很可能是所观察免疫变异的关键贡献者，但它是在对多种流行病原体累积暴露的背景下发挥作用。因此论文把结果描述为环境性免疫适应，而不是疟疾特异性效应。

这一范围条件改变了结果的使用方式。658个基因、模块方向和整合评分描述的是这四个肯尼亚环境中的成人；它们不能确定哪种病原体驱动了哪条通路，也不能证明单是疟疾就产生了这种状态。作者把这一模式与既往报道的高流行环境中适应性免疫扩增、以及疟疾暴露成人中炎症单核细胞表型受抑联系起来，但这些比较只是解释背景，不是本研究新产生的证据。

基线状态可能影响后续应答

论文的前瞻性主张是有条件的：这种独特的基线免疫状态可能影响对后续感染和疫苗接种的应答。这是关于后果的假设，不是实测结局。研究没有检验这些参与者对感染或疫苗如何应答，因此不应把结果读成高暴露成人应答更差或更好的证据。

最有用的收获是一个比较关系，而不是单个数字。在这60名健康成人中，较高的累积病原体暴露与B细胞和抗原呈递程序上调、炎症趋化和血管转录下调相伴随，而血浆介质的变化方向与细胞炎症信号轻度相反。因此，任何评估流行地区免疫基线数据的人，都应预期适应性活化和炎症抑制同时出现，也不应把低炎症转录谱当作全身介质水平低的证据。

本文引用的文献