关于儿童早期1型糖尿病中T细胞状态失衡的研究报告
本文旨在介绍一项发表于《自然·通讯》（*Nature Communications*）期刊（2025年，卷16，文章号11301）上的原创性研究。该研究由来自捷克科学院分子遗传学研究所（Institute of Molecular Genetics of the Czech Academy of Sciences）的 Veronika Niederlova、Aleš Neuwirth、Ondřej Štěpánek 及其合作者（包括来自查理大学和第二医学院、莫托尔大学医院的合作者）共同完成。
一、 研究背景与目的 1型糖尿病（Type 1 Diabetes, T1D）是一种自身免疫性疾病，其根本原因是免疫系统对产生胰岛素的胰腺β细胞失去了自我耐受。尽管遗传和环境因素已知参与其中，但T1D的具体病因，特别是免疫系统失调的详细细胞和分子机制，仍未完全阐明。T细胞，尤其是调节性T细胞（Regulatory T cells, Treg）和效应性T细胞（Effector T cells）的平衡，在维持自我耐受和自身免疫病理中扮演核心角色。既往研究提示T1D患者可能存在T细胞功能异常，但结论不一，且缺乏在疾病早期、特别是儿童发病初期的系统性、高分辨率研究。
本研究旨在利用单细胞转录组测序（single-cell RNA sequencing, scRNA-seq）技术，结合流式细胞术和公共数据集分析，深入探究新诊断T1D儿童患者外周血T细胞的特征，并与健康儿童对比。研究重点关注两个时间点：诊断时（T0）和诊断后一年（T1），以纵向追踪疾病早期的免疫动态变化。核心科学问题是：T1D儿童的T细胞在转录组和表型上是否存在系统性失调？这种失调是否涉及效应功能、干细胞样特性以及调节功能？研究目标是为理解T1D早期免疫病理机制提供新的见解。
二、 研究流程与方法 本研究是一个综合性的前瞻性队列研究，结合了scRNA-seq、T细胞受体（TCR）测序、流式细胞术验证以及公共数据再分析。整体工作流程严谨且多层次。
流程一：队列建立与样本采集 研究纳入了43名儿童供者，包括30名新诊断的T1D儿童（诊断中位时间9天）和13名年龄匹配的健康对照。对T1D儿童在诊断时（T0）和诊断后一年（T1，中位时间377天）分别采集了外周血样本。收集了详细的临床数据，包括年龄、性别、糖尿病酮症酸中毒状态、自身抗体水平、HLA单倍型、血糖控制指标（如糖化血红蛋白、C肽）和白细胞计数。从外周血中分离出外周血单个核细胞（PBMC）并进行冷冻保存，以确保后续实验批次间的一致性。
流程二：单细胞转录组测序实验设计 为了克服传统bulk测序缺乏单细胞分辨率以及流式细胞术受限于预设标记的缺点，研究团队采用了scRNA-seq（10x Genomics平台）对分选的CD4+和CD8+ T细胞进行分析。实验分为“初始实验”和“最终实验”两个阶段。 * 初始实验：对少量样本进行scRNA-seq，发现儿童供者中大部分T细胞呈初始（Naïve）表型（CD45RA+ CD45RO-）。为了更详细地研究抗原经验（非初始）T细胞，在最终实验中采用了富集策略：对每个样本，按1:5的比例分选初始（CCR7+ CD45RA+）和抗原经验T细胞，再进行scRNA-seq。这种针对性的富集策略是本方法学的关键创新点，它确保了在测序深度有限的情况下，能对相对稀少的抗原经验T细胞亚群（如效应记忆、终末分化效应细胞等）进行充分解析。最终，共处理了176个样本，平均每个样本测序约992个细胞。
流程三：单细胞数据分析流程 数据处理采用了标准化的生物信息学流程，并整合了多种先进算法以确保结果的可靠性。 1. 质量控制与预处理：使用Seurat软件包进行，去除低质量细胞、双细胞和死细胞。 2. 批次效应校正：由于文库制备分六个批次进行，研究使用STACAS算法进行整合，以消除技术批次效应，同时保留真实的生物学差异。 3. 细胞图谱构建与注释：对CD8+和CD4+ T细胞分别进行无监督聚类（UMAP降维、Louvain聚类），并依据典型标记基因（如 TCF7, LEF1, FOXP3, GZMB, NKG7 等）定义了主要的T细胞亚群。例如，CD8+ T细胞被分为初始（Naïve）、记忆（Tmem）、终末分化效应记忆（TEMRA）、增殖细胞等；CD4+ T细胞被分为初始、中央记忆（Tcm）、调节性T细胞（Treg）、Th1/Th17、Th2、TEMRA等亚群。细胞身份通过已发表的RNA-seq数据集进行了验证。 4. 差异表达与功能分析：比较了健康对照、T1D T0和T1D T1组间的基因表达差异。使用了Wilcoxon秩和检验、DESeq2（基于伪bulk表达谱）等多种统计方法。进行了基因集富集分析（GSEA），以评估特定功能通路（如效应功能、干细胞特性）的富集情况。 5. 轨迹推断与功能状态评分：对Treg细胞使用Slingshot和Condiments软件包进行伪时间轨迹分析，推断其分化状态和激活程度。使用decoupler等方法计算细胞的功能模块得分（如效应得分、干细胞样得分）。 6. TCR序列分析：从scRNA-seq数据中提取配对的TCRα和TCRβ序列，分析克隆扩增、CDR3长度和生化特性（如Boman指数、Hmoment、疏水性）。
流程四：流式细胞术验证与公共数据再分析 为了验证scRNA-seq的发现并在独立队列中确认，研究进行了两方面的补充： 1. 流式细胞术验证：使用多色流式细胞术对同一队列样本的T细胞表型（如初始/效应比例、颗粒酶B表达）和特定亚群（如Treg、γδ T细胞、CD3+CD56+ Tr3-56细胞）的频率进行了定量分析。这提供了蛋白质水平的直接证据。 2. 公共数据集再分析：为了增强结论的普适性，研究团队系统地再分析了多个公开可用的数据集，包括人类胰腺分析计划（HPAP）的scRNA-seq和流式数据，以及其他已发表的T1D相关转录组学研究。使用统一的流程处理这些数据，以验证本研究中观察到的关键特征（如效应基因下调、干细胞样特征上调）是否在其他独立队列中重现。
三、 主要研究结果 研究结果揭示了儿童早期T1D患者外周血T细胞存在深刻的、多层面的功能状态失衡。
结果一：T细胞效应程序减弱与干细胞样特征增强 这是本研究的核心发现。scRNA-seq分析显示，与健康儿童相比，新诊断的T1D儿童（T0）的CD4+和CD8+ T细胞中，与效应功能和细胞毒性相关的基因（如 GZMA, GZMB, PRF1, IFNG, *CCL5*）表达显著下调。相反，与初始、静息和干细胞样记忆T细胞状态相关的基因（如 TCF7, LEF1, IL7R, *CCR7*）表达上调，研究者将其统称为“干细胞样相关基因特征”。这种向干细胞样/初始表型的偏倚在诊断时最为明显，在诊断一年后（T1）有所缓解但未完全恢复正常。GSEA分析进一步证实，T1D儿童的T细胞富集了初始T细胞特征，而健康对照的T细胞富集了效应T细胞特征。 * 支持数据：差异表达基因热图、GSEA曲线图、基于参考数据集的decoupler评分图（均显示T1D组干细胞样特征得分高，效应特征得分低）。 * 逻辑关系：此发现提示T1D早期可能存在T细胞分化或激活障碍。为了排除scRNA-seq样本富集策略造成的偏差，研究者进行了流式细胞术验证。 * 流式验证结果：流式数据证实，在蛋白质水平上，T1D儿童T0时抗原经验T细胞比例低于健康对照，且CD8+和CD4+ T细胞内的颗粒酶B（GZMB）蛋白表达水平也较低。这与转录组数据高度一致。 * 公共数据验证：对HPAP数据库（脾脏T细胞）和其他已发表研究的再分析，均重复出了“健康对照效应特征强，T1D患者干细胞样特征强”的模式，极大地增强了这一发现的可靠性和普适性。
结果二：调节性T细胞（Treg）功能受损与成熟障碍 对Treg的深入分析揭示了其功能状态的异常。 1. Treg亚群异质性：通过重新聚类，将Treg分为四个亚群（Treg1-4）。Treg1表现出最初始/干细胞样的表型（高表达 CCR7, *TCF7*，相对低表达 FOXP3, *CTLA4*），而Treg4则表现出最活化的效应表型。 2. Treg组成偏移：在T1D儿童中，初始样Treg1亚群的比例显著高于健康对照，而活化程度高的Treg4亚群则相对较少。伪时间轨迹分析显示，T1D儿童的Treg整体上处于更不成熟、激活程度更低的状态。 3. Treg功能失调特征：T1D儿童的Treg基因表达谱与已知的功能失调Treg相似。例如，其基因特征与IPEX综合征（由 FOXP3 突变引起）患者的缺陷性Treg样细胞有重叠。此外，T1D Treg中与功能障碍相关的基因（如 CD226, IL4R, *GZMK*）表达较高，而与抑制功能相关的基因（如 *IL10RA*）表达与残留的β细胞功能（空腹C肽水平）呈正相关。 4. 独立验证：在三个独立的公共scRNA-seq数据集中，研究者均观察到了类似模式：T1D患者的Treg中初始样亚群富集，而效应亚群减少。
结果三：非常规T细胞（特别是具有调节潜能的亚群）频率与功能改变 研究者关注了此前被认为可能与T1D相关的CD3+CD56+ Tr3-56细胞。分析发现，Tr3-56表型并非一个独立的细胞群，而是主要与CD8+的黏膜相关恒定T细胞（MAIT）、γδ T细胞以及CD4+ TEMRA细胞重叠。 1. 频率降低：流式细胞术显示，T1D儿童在诊断时，外周血中γδ T细胞和CD8low（富含γδ T和CD56+ T细胞）细胞的比例低于健康对照。 2. 效应功能减弱：与常规T细胞类似，T1D儿童的这些非常规T细胞（如γδ T细胞）也表现出效应/细胞毒性基因（如 *GZMB*）表达下调的趋势。效应表型强的γδ T1亚群在健康对照中更丰富。 3. 调节潜能：具有Tr3-56表型的细胞在所有CD8+淋巴细胞亚群中，与经典Treg特征基因的相似性最高，提示其可能具有调节功能。
结果四：其他发现 1. BTN3A2表达与HLA单倍型相关：研究发现，此前有报道在T1D中上调的 BTN3A2 基因，其表达差异主要与特定的HLA单倍型相关，而非疾病状态本身，强调了对照选择在遗传关联研究中的重要性。 2. TCR受体库分析：未发现T1D特异性的“公共”克隆。同一患者在T0和T1时间点存在显著的克隆重叠，表明克隆持续性。未重复出此前报道的CDR3β长度差异，但发现了CDR3β序列生化特性（如Boman指数、Hmoment）的细微差异，需进一步大样本验证。
四、 研究结论与意义 本研究通过高分辨率的单细胞多组学分析，系统地描绘了儿童早期T1D患者外周血T细胞的转录组和表型图谱，并得出以下核心结论： 1. T1D与T细胞分化状态失衡相关：在疾病早期，患者的T细胞整体表现出效应/细胞毒性程序减弱和干细胞样/初始程序增强的特征。这种状态在诊断一年后部分正常化，但未完全恢复。这挑战了T1D solely由过度活跃的效应T细胞驱动的简单观点，提示可能存在T细胞分化或激活的内在障碍。 2. 调节性T细胞功能不成熟：T1D患者的Treg池偏向于不成熟的初始样状态，并表现出与已知功能失调Treg相似的特征，表明其调节能力可能受损。这种不成熟状态可能与无法有效控制针对β细胞的自身免疫反应有关。 3. 非常规调节性T细胞亚群减少：具有潜在调节功能的非常规T细胞（如某些γδ T和MAIT细胞）在T1D患者中频率降低且效应表型减弱，可能进一步削弱了免疫调节网络。
科学价值： * 机制新见解：提出了T1D早期免疫失调的新范式——不仅是效应T细胞的过度攻击，更是整个T细胞区室（包括效应和调节臂）的分化与功能状态失衡。干细胞样表型的偏倚可能与“卫生假说”相联系，即早期感染暴露不足导致免疫系统“不成熟”。 * 纵向动态：首次在新诊断儿童中进行了为期一年的纵向追踪，揭示了疾病早期免疫状态的动态变化，为理解疾病自然史和干预时间窗口提供了线索。 * 高分辨率图谱：生成了首个针对儿童T1D早期的高深度T细胞单细胞图谱，为领域提供了宝贵的资源。 * 方法论严谨性：通过富集策略、多批次整合、流式验证和多数据集交叉验证，确保了研究发现的稳健性。
潜在应用价值： * 生物标志物：发现的T细胞特征（如干细胞样特征评分、Treg成熟度评分）可能作为疾病活动或治疗反应的潜在生物标志物。 * 治疗靶点：研究提示，促进T细胞向效应方向正常分化，或增强Treg的功能成熟，可能成为新的免疫治疗策略。例如，靶向初始样Treg（如通过调节IL-4R或CD226信号）来促进其功能成熟。 * 预防策略：支持了早期环境因素（如微生物暴露）在塑造免疫系统、预防自身免疫中的重要性。
五、 研究亮点 1. 研究对象独特：聚焦于新诊断的儿童T1D患者及其纵向随访样本，抓住了疾病发生和早期演进的关键时间窗口，数据极具时效性和临床相关性。 2. 技术方法系统先进：结合了scRNA-seq、TCR测序、流式细胞术和公共数据再分析，形成了多层次、相互验证的研究体系。采用的抗原经验T细胞富集策略是获得深层亚群信息的关键技术亮点。 3. 发现具有颠覆性和整合性：揭示了T1D早期T细胞全局性的“低效应、高干细胞样”状态，并同时将常规Treg和非常规调节性T细胞的功能异常纳入统一框架，提供了一个更全面的免疫失调视图。 4. 结论稳健可信：所有主要发现均在内部队列通过流式细胞术验证，并在多个独立公共数据集中得到重复，显著提升了结论的可靠性和普适性。 5. 深入机制探索：不仅描述了现象，还通过轨迹分析、功能失调基因特征比对等，深入探讨了Treg功能状态异常的潜在机制。
六、 其他有价值内容 研究还澄清了领域内一个不一致的发现（BTN3A2 表达差异主要源于HLA背景而非疾病），体现了精细分析的重要性。此外，对TCR库的细致分析虽未发现强疾病关联克隆，但为未来更大规模的抗原特异性研究奠定了基础。研究者也坦诚指出了研究的局限性，如样本量相对较小、性别分布不完全匹配等，并强调本研究是探索性的，需要后续研究进一步证实和拓展。