T1D血液来源B细胞特征与桥本甲状腺炎组织的跨组织一致性：探索性空间免疫主题分析

《International Journal of Molecular Sciences》：Cross-Tissue Concordance of a T1D Blood-Derived B-Cell Signature with Hashimoto Thyroiditis Tissue: Exploratory Spatial Immune-Topic Analysis

　　

摘要专业翻译 1型糖尿病（type 1 diabetes, T1D）与桥本甲状腺炎（Hashimoto thyroiditis, HT）具有共同的自身免疫特征，但B细胞转录程序的跨组织一致性仍不确定。研究人员利用24项保留的队列–通路效应估计值，重新计算了T1D血液汇总级meta分析统计量，并另行报告了先前生成的固定14基因HT甲状腺关联。T1D meta分析得出Hedges' g = 0.51（95% CI 0.18–0.85；q = 0.022；I2 = 0%）。在HRA001684（16例HT，50例非HT）中，固定评分与HT相关（HC3校正β = 2.65，95% CI 2.37–2.94；置换P = 2 × 10?5）。单细胞数据（6512个细胞；三个HT和三个免疫检查点抑制剂相关甲状腺炎文库）缺乏健康对照，仅作描述性分析。在五个空间切片中，固定签名对比无显著差异（Δ = 0.459；精确P = 0.10；匹配基因P = 0.085），与B细胞/浆细胞主题无优先对应关系（报告P = 0.90）。独立推断的甲状腺免疫主题显示区域变异，但未定位固定评分。这些报告结果与T1D血液–HT甲状腺转录部分一致性相符，但未确立全身性HT激活、共享的细胞内在程序或经验证的空间定位。在报告的敏感性分析中，联合B细胞/浆细胞标记物调整使HT关联衰减约90%，完全调整后的置信区间包含零。与等权重及其他B细胞参考签名的比较未显示T1D来源权重具有统计学显著优势。公开GitHub工作流仅验证保留效应的汇总级合并；HT批量、空间、单细胞及标记物/权重分析的源到结果重建尚未建立。

论文解读文章

一、研究背景与问题提出

1型糖尿病（type 1 diabetes, T1D）与桥本甲状腺炎（Hashimoto thyroiditis, HT）是两种常见的自身免疫性内分泌疾病，临床上常合并存在。两者虽累及不同内分泌器官，但均涉及免疫耐受受损、HLA相关易感性、淋巴细胞活化、抗原呈递及干扰素相关炎症等共同机制。既往研究多聚焦于循环免疫细胞或蛋白水平的跨疾病比较，例如Magnusson等人的质谱流式研究报道HT患者循环浆母细胞增加，血浆蛋白分析发现CDCP1在T1D和HT中均升高、SLAMF1在长期T1D和HT中升高，这些结果确立了免疫学重叠，但均未证明完全相同的B细胞转录程序在两疾病间得以保留。此外，Yi等人对六种自身免疫或炎症性疾病（包括T1D胰岛和HT甲状腺组织）的批量及单细胞/单核数据进行分析，报道了趋同的炎症通路，但同样未回答一个预设的T1D血液来源B细胞评分是否与HT甲状腺状态相关这一具体问题。因此，研究人员提出核心科学问题：源自T1D外周血的固定B细胞转录特征是否与独立的HT甲状腺组织状态存在关联，并进一步探索该特征在空间转录组层面的分布模式。

二、研究设计与主要发现

研究人员整合了多个公开人类转录组数据集开展系统性分析。三个T1D外周血队列（GSE30210、GSE9006、GSE156035）用于通路水平meta分析和固定14基因加权签名构建；HRA001684提供HT甲状腺批量队列（16例HT，50例非HT）用于跨组织评分验证；GSE192580提供甲状腺炎单细胞数据（6512个细胞，三个HT和三个免疫检查点抑制剂相关甲状腺炎文库，无健康对照）；GSE248205提供空间转录组切片（三个HT、两个对照、三个Graves病切片，其中HT-对照直接推断仅用五个切片）；GSE230424提供外部描述性空间图谱（两个HT-only切片P3/P4和两个癌相关HT切片P1/P2）。

主要研究结果如下：

T1D外周血B细胞程序meta分析：基于24项保留的队列–通路效应估计值，随机效应meta分析显示B细胞程序汇总效应为Hedges' g = 0.51（95% CI 0.18–0.85；q = 0.022；I² = 0%）。B细胞和JAK–STAT通路具有最大的正向汇总估计，T细胞、HLA I类、I型干扰素和细胞因子/炎症程序亦为正向，而HLA II类/抗原呈递细胞（APC）估计较弱且为负向。独立标准库检查器同意全部八个保留的合并通路结果。

固定T1D签名在HT批量组织中的关联：通过逆方差合并三个T1D队列的基因水平Hedges' g并做L1归一化，构建了14基因加权签名。在HRA001684中，该固定评分与HT状态显著相关（HC3校正β = 2.6534，95% CI 2.3679–2.9389），Freedman–Lane残差置换和均值/方差匹配基因零检验均得到经验P = 2 × 10^?5。

探索性RNA优先空间免疫主题表征：从GSE248205八个切片推断出12个转录组主题，其中五个免疫相关主题（T1、T2、T7、T8、T9）显示重叠但不同的B细胞、浆细胞、HLA、抗原呈递、炎症和补体特征。主题评分在发现集和外部切片间存在变异。外部GSE230424的HT-only P3/P4切片用于说明性图谱。APC/HLA-II相关主题7评分在HT1、P3和P4内显示空间异质性。H&E关联比较显示数值上多数切片主题7相关性更高，但变化异质且存在一个负外部相关。

固定签名空间投影与单细胞描述性背景：固定评分在GSE248205五个切片中的HT–对照差异不显著（Δ = 0.459；精确切片标签P = 0.10；匹配基因零P = 0.085），与预设B细胞/浆细胞主题的对应关系亦未获支持（P = 0.90）。GSE192580单细胞数据描述了15个原始Seurat聚类中的谱系标记物表达，B细胞标记物见于聚类10和12，浆细胞标记物见于聚类11，但该数据集缺乏健康对照，仅作描述性解释。

B细胞标记物调整与签名权重敏感性分析：原始签名在调整年龄、性别和文库制备批次后HT系数为β = 2.65。B细胞标记物代理调整使系数降至β = 0.45，浆细胞代理单独调整产生较小降低，联合调整进一步衰减，完全调整后系数接近零且置信区间包含零。报告的自举分布显示联合调整后约90%衰减。随机权重基准显示残差偏R²值存在实质变异，原始固定权重为0.0187，等权重版本为0.0111，Travaglini肺B细胞参考签名为0.037，Aizarani肝B细胞参考签名为0.003。这些比较未显示T1D来源权重具有统计学显著优势。

三、关键技术方法

研究人员主要采用以下关键技术方法：跨队列基因水平效应合并采用固定效应逆方差加权meta分析，通路水平采用DerSimonian–Laird随机效应合成并辅以Benjamini–Hochberg校正；固定签名评分通过基因稳健标准化后加权求和实现；HT批量关联采用协变量调整普通最小二乘回归结合HC3稳健协方差推断；显著性检验采用Freedman–Lane残差置换和均值/方差匹配基因竞争性零检验；空间分析采用RNA优先主题推断、切片水平标签置换检验及余弦相似度跨数据集主题匹配；单细胞分析采用Seurat聚类和谱系标记物描述。所有分析基于公开数据集，样本队列来源包括GSE30210、GSE9006、GSE156035、HRA001684、GSE192580、GSE248205和GSE230424。

四、讨论与结论总结

讨论部分指出，三个T1D血液队列中B细胞相关程序显示一致的汇总效应，固定14基因评分与HT批量甲状腺状态相关，这是主要的生物学观察，但需要独立重建固定权重、表达标准化和样本水平HT评分才能实现完全计算可重复性。单细胞和RNA优先空间分析表征了相关免疫特征，但均未直接检验转移评分的细胞来源或空间定位。标记物代理调整使HT关联显著衰减，可用基准结果未显示T1D来源权重的优越性。研究人员强调，该关联可能反映趋同炎症通路、局部B细胞丰度或细胞内在激活，当前批量及描述性空间数据无法区分这些解释。研究未在HT外周血、配对血液–甲状腺样本或T1D胰岛转录组中分析固定签名，也未进行留一基因或留一队列稳健性分析。空间推断仅用三个HT和两个对照切片，外部说明性图谱来自含癌相关HT的集合子集。公开代码仓库仅支持保留效应meta分析的汇总级验证，图3、图4和图5的端到端独立复现尚未建立。

结论部分总结：源自T1D血液队列的固定14基因签名与HT批量甲状腺状态相关。甲状腺炎单细胞分析在无健康对照条件下描述了免疫谱系标记物，独立推断的空间主题显示区域异质性。五切片固定评分对比不显著（精确P = 0.10；匹配基因P = 0.085），与B/浆细胞主题的优先对应关系未获支持（P = 0.90）。这些报告观察与部分血液–甲状腺转录一致性相符，但未确立全身性HT B细胞激活、转移评分定位、共享因果机制或HT关联的独立端到端复现。标记物代理敏感性分析中，联合调整使HT系数衰减约90%，完全调整后置信区间包含零。与等权重及外部B细胞参考签名的比较未显示T1D来源签名权重的优越性。对于生物医学研究者，固定评分在签名权重、基因归一化和完整可执行评分工作流经验证并发布后，可作为外部队列评估的假设。临床导向研究可在配对或纵向血液–甲状腺队列中前瞻性评估，但独立验证和增量临床价值证据是诊断、风险分层、预后或治疗选择考量之前的必要条件。