白介素10启动子变异与精神病夸大症状的机器学习基因集研究

原文：Interleukin pathway genotype and the grandiosity dimension of psychosis: a machine learning gene set analysis with unsupervised anomaly detection controls.
期刊：Frontiers in immunology（影响因子 7.0） |  PubMed | PMC 全文（含 PDF）

在精神分裂症等精神病性障碍的研究中，免疫系统参与疾病发生早已不是新鲜话题。流行病学证据显示，母亲孕期感染会增加后代罹患精神病的风险；遗传学上，精神分裂症最强的常见变异信号位于主要组织相容性复合体，并由补体成分4（C4）的表达差异所解释——风险等位基因提高C4表达，而C4会标记突触以供小胶质细胞清除。成年患者外周血细胞因子也呈现可重复的改变：一部分在急性期升高、随治疗恢复正常，另一部分则不论临床状态如何都持续偏高，形成所谓“状态与特质并存”的免疫表型。

然而，既往免疫研究大多聚焦于阴性症状、缺陷症状和偏执症状，对“夸大（grandiosity）”这一维度几乎无人问津。夸大信念——对自身力量、知识、身份或价值的膨胀感——在约一半的精神病患者中出现，在躁狂患者中更可达三分之二，但其生物学基础基本未知。更关键的是，关于炎症如何驱动精神行为，学界存在两种对立假说：一种认为促炎状态驱动“激活”，从而应增加夸大；另一种认为炎症抑制纹状体多巴胺信号，导致快感缺失和精神运动迟滞，从而应减少夸大。阴性症状无法区分这两种假说，而夸大恰好是一个判别性症状。

近期发表在 Frontiers in Immunology（影响因子 7.0）上的一项研究，正是利用波兰波兹南（Poznań）精神病遗传学队列，对489例精神分裂症患者进行了白介素通路基因集与夸大维度的机器学习分析，并引入了无监督异常检测作为对照。论文题为《Interleukin pathway genotype and the grandiosity dimension of psychosis: a machine learning gene set analysis with unsupervised anomaly detection controls》。

研究背景：为什么盯住“夸大”和白介素通路

研究者组装了一个覆盖促炎与抗炎两类细胞因子及其受体、拮抗剂的基因组合，共8个白介素通路基因中的26个功能性变异：包括 IL1A、IL1B、IL1RN、IL6、IL6R、IL10、IL10RA 和 TGFB1。这些位点均为启动子、编码或调控区功能变异，用以反映个体稳定的细胞因子产生“设定点”。其中，白介素-10（IL-10）是该网络中的主要抗炎细胞因子，其启动子多态性在既往荟萃分析中被认为轻微增加精神分裂症风险。

夸大症状的操作化定义来自 OPCRIT（Operational Criteria Checklist）量表：将“ inflated self-esteem（ inflated 自尊，第56项）”与“ grandiose delusions（夸大妄想，第57项）”两项相加，每项0–2分，总分0–4分。489名患者中，73%的人得分为0，即终生从未出现可观察的夸大表现；分析主要使用“是否出现过任何夸大（得分>0）”以及经秩逆正态变换后的连续分数两种形式。

核心发现：整体基因面板无效，但 IL10 启动子位点显著

研究首先把整个26标记面板作为整体，用岭回归在五折交叉验证中预测夸大分数。结果显示，面板整体的样本外 Spearman 相关仅为 rho = 0.041，置换检验 p = 0.159，即不能预测夸大。多种监督分类器（逻辑回归、弹性网络、随机森林、梯度提升、支持向量机）的汇总ROC也未见有效预测能力。图1顶部展示了面板水平检验：按夸大着色的前两个基因型主成分中，夸大患者并未形成可分离的聚类。

随后，研究者以基因为单位，用 SKAT（sequence kernel association test）进行方差成分检验，并控制错误发现率（FDR）及进行选择校正的置换检验。八个基因中，只有 IL10 通过阈值：SKAT p = 0.0028，q = 0.023；选择校正后 p = 0.025。进一步定位发现，信号局限于 IL10 的 -819/-592 启动子位点（由两个近乎完全连锁不平衡的标签 SNP rs1800871 与 rs1800872 构成，剂量相关 r = 0.994，被合并为单一单倍型剂量）。

该位点呈现隐性效应：携带两份拷贝（即 -819/-592 共享单倍型纯合）的患者，用 Firth 惩罚逻辑回归估计的出现夸大的优势比（OR）高达 5.28（95% CI 2.37–12.48，p = 4.5×10⁻⁵）。在26名两份拷贝携带者中，65% 有夸大表现；而非携带者中仅 24% 有夸大。效应特异于夸大维度（标准化 beta = 0.19，p = 2×10⁻⁵），在随机分半样本中中等可重复（r = 0.36；以躁狂簇为对照时 r = 0.03）。该效应不依赖性別、评估年龄或发病年龄，在替换 SKAT 核、加入面板主成分、调整神经阻滞剂反应（OPCRIT 第89项）后均不变。

值得注意的是，在466例患者与482例对照之间，两份拷贝携带频率并无差异，说明该 IL10 启动子变异并非精神分裂症病例-对照风险因子，而是与“夸大”这一症状维度特异相关。图2以交叉拟合方式展示了每位患者在白介素面板上的贡献：489行按夸大分数从高到低排列，可直观看到高夸大患者集中在 IL10 相关标记的特定贡献模式上。

方法亮点：从聚合预测到基因集检验，再到无监督异常检测

这项研究在方法上有一套严格的“门控”流程。由于面板中预期没有单个标记能携带可检测效应，证据先在基因内聚合：SKAT 作为方差成分得分检验，能把一个基因内多个变异（甚至方向相反者）汇总为单一统计量。作者采用线性核与等权重，因为面板是常见功能变异，而非罕见变异测序数据；若使用针对罕见变异的 Beta(1,25) 权重，反而会把检验集中到最罕见标记上，而此处两份拷贝携带者其实是常见等位基因（次要等位基因频率 0.23）的纯合子。

为控制多重检验，八个基因采用 Benjamini–Hochberg FDR；由于报告的是八者中 p 值最小者，还进行了“选择校正置换”：把夸大表型打乱，重跑整个八基因筛选，比较观察到的最小 p 值与每次打乱中的最小值。此外，用重复“发现—留出”分半、刀切法（jackknife over carriers）、脆性分析（fragility analysis：需要多少名夸大携带者被重编码或移除才会使 p>0.05）、以及300次随机留 k 出子集来检验稳定性。

无监督方向是本文另一创新。研究者应用了四种异常检测器：隔离森林（isolation forest）、单类支持向量机（径向基核）、全连接神经自编码器（26-16-4-16-26，瓶颈4单元，Adam优化，L2惩罚，3000次迭代，以均方重建误差评分），以及另一集成方法。目的是在不参考夸大表型的情况下，先找出整体基因型谱“异常”的患者，再看他们是否在夸大上不同——这对应工程状态监测中的异常检测逻辑。图3顶部展示了每位患者的集成异常分数，以及各检测器在包含 IL10 标记时异常分数与夸大的相关性。

由于该队列没有全基因组基因分型，无法计算祖先主成分，作者用面板自身允许的两个检查替代：比较病例与对照间26标记基因型矩阵的前五个主成分；以及在加入24个非 IL10 标记的前五个主成分作为协变量后重新估计 IL10 基因集检验和隐性 OR。所有模型超参数事先固定，不根据结果调整；缺失剂量（12,714个基因型细胞中169个，1.33%）在拟合任何模型前一次性均值填补，且不使用表型信息，因而不会在交叉验证折间泄漏结果信息。

应用意义：为“抗炎—激活”假说提供可检验的候选机制

这项研究的意义在于，它把免疫遗传学从“是否患病”推进到“何种症状维度”。整体面板不能预测夸大，说明常见的白介素功能变异并不以多基因聚合方式普遍影响该症状；但 IL10 -819/-592 启动子位点的隐性效应，则提示抗炎细胞因子 IL-10 的遗传设定点可能特异调控夸大这一“激活”行为表达。

由于两份拷贝携带者在病例与对照间频率相同，该变异不是精神分裂症易感标记，而可能是症状修饰因子：在已患病者中，IL-10 低产生（或特定启动子单倍型）相关的隐性遗传背景，使个体更易表现夸大。作者谨慎地将其作为“待独立重复的假设”提出，而非既定结论——样本仅489例、携带者仅26人，且队列为单一波兰站点，需要其他独立样本验证。

对科研而言，该工作示范了如何把机器学习预测、基因集检验、惩罚回归与无监督异常检测组合成一个可重复、防过拟合的分析流水线；所有数字与图形均可通过固定随机种子的 Jupyter 笔记本复现，并遵循 STREGA 与 TRIPOD+AI 报告规范。对临床而言，若未来重复成功，IL10 启动子基因型或有助于识别精神病患者中更易出现夸大、激活症状的亚群，从而为理解炎症—行为关系及潜在免疫调节治疗分层提供线索。

总之，这项研究在“免疫参与精神病”的版图上补上了一块长期缺失的拼图：不是阴性或偏执症状，而是夸大维度，可能藏着白介素-10 通路与大脑激活状态之间的特异性联系。

本文图片转载自 PMC 开放获取全文，仅供学习参考，请以原文为准。