《Frontiers in Cell and Developmental Biology》：Comprehensive profiling of immune landscape changes and oxidative stress-associated biomarkers in intervertebral disc degeneration
编辑推荐：
背景： 椎间盘退变（intervertebral disc degeneration, IDD）常被认为是腰痛的原因之一。氧化应激与免疫微环境失调已被公认为该病理过程中的重要参与者。然而，关于与氧化应激相关的核心基因及其如何参与免疫应答的整体认识仍然缺乏。
方
背景： 椎间盘退变（intervertebral disc degeneration, IDD）常被认为是腰痛的原因之一。氧化应激与免疫微环境失调已被公认为该病理过程中的重要参与者。然而，关于与氧化应激相关的核心基因及其如何参与免疫应答的整体认识仍然缺乏。
方法： 研究人员整合了来自三个IDD队列（GSE167199、GSE176205、GSE186542）的转录组数据，并进行了差异表达分析、加权基因共表达网络分析（WGCNA）、共识聚类以及CIBERSORT分析，以鉴定氧化应激相关基因、定义分子亚型并评估免疫浸润情况。研究人员进一步构建了单基因基因集富集分析（GSEA）模型、竞争性内源RNA（ceRNA）网络以及转录因子-信使RNA（TF-mRNA）相互作用网络。通过在人样本中进行定量聚合酶链反应（qPCR）和免疫组织化学（IHC）检测，并在大鼠穿刺模型中通过免疫荧光（IF）实验验证了生物标志物的有效性。
结果： 该研究揭示了两种具有不同氧化应激特征的新型IDD亚型，并鉴定出CD36、SESN3和HMOX1为关键调控因子。它们的表达与免疫浸润，特别是与静息肥大细胞、初始CD4+T细胞和巨噬细胞的浸润程度相关。GSEA分析将其与细胞周期、代谢和细胞因子通路联系起来，而ceRNA和TF网络则揭示了涉及长链非编码RNA（lncRNA）、微小RNA（miRNA）和转录因子的复杂调控机制。验证实验证实了这三种基因在人及大鼠退变椎间盘组织中的表达上调，突显了其诊断和治疗潜力。
结论： 该研究利用生物信息学方法鉴定并提出了CD36、SESN3和HMOX1作为IDD中潜在的氧化应激相关生物标志物，并强调了它们与免疫微环境重塑之间的密切关系。这些结果有助于更好地理解椎间盘退变过程中氧化过程与免疫应答之间的关系，并可能有助于构思未来具有抗氧化和免疫调节作用的治疗策略。
研究背景与科学问题
腰痛是全球范围内导致功能障碍的主要原因之一，而椎间盘退变（intervertebral disc degeneration, IDD）是其主要的病理基础。随着全球人口老龄化加速，与IDD相关的临床负担和经济成本日益增加。从结构上看，椎间盘是一种特殊的纤维软骨组织，血供稀疏。髓核细胞负责合成和维持富含II型胶原和蛋白聚糖的细胞外基质（extracellular matrix, ECM），以保持椎间盘的水合水平和机械稳定性。然而，老化、异常机械负荷和营养剥夺等因素会逐渐破坏椎间盘的内稳态，导致不可逆的退行性改变。
在IDD的众多发病机制中，氧化应激被认为是疾病起始和进展的关键驱动因素。尽管基础水平的活性氧（reactive oxygen species, ROS）对细胞内通讯是必要的，但其过度产生会扰乱氧化还原平衡，引发氧化应激。高水平的ROS不仅直接损伤脂质、蛋白质和DNA，还会激活MAPK、NF-κB和PI3K/Akt等关键信号级联反应，促进纤维环和髓核细胞的衰老与凋亡。同时，氧化应激会下调合成代谢活动，上调基质降解酶如ADAMTS家族成员和基质金属蛋白酶（matrix metalloproteinases, MMPs）的表达，从而加速ECM的分解。此外，这一过程还通过激活NLRP3炎症小体通路增强慢性炎症反应，驱动促炎介质如肿瘤坏死因子-α（tumor necrosis factor-alpha, TNF-α）、白细胞介素-6（interleukin-6, IL-6）和白细胞介素-1β（interleukin-1 beta, IL-1β）的分泌，形成氧化应激、持续性炎症和组织进行性降解的恶性循环。
近年来，越来越多的关注聚焦于IDD进程中局部免疫微环境的改变。在生理健康条件下，椎间盘被视为免疫豁免部位。但随着退变的进展，终板屏障受损，伴随病理性血管化和神经生长，促进了免疫细胞尤其是巨噬细胞和T淋巴细胞向椎间盘基质的浸润。这些浸润的免疫细胞与原生椎间盘细胞相互作用，进一步加剧退变级联反应。重要的是，氧化应激在调节免疫细胞的募集、功能活性和极化（特别是巨噬细胞）中也扮演着关键角色。ROS会促使巨噬细胞向促炎的M1表型转化，进而增加ROS和炎性细胞因子的产生，加重组织损伤。
尽管已有大量文献指出氧化应激相关基因参与IDD，但以往的研究往往是孤立地分析单个数据集或队列，样本量有限，这限制了对连接氧化应激与免疫调控的重要分子网络的阐明，也阻碍了稳健生物标志物的识别和靶向治疗策略的设计。因此，有必要对IDD背景下氧化应激框架内的分子进行详细分析，以更深入地理解该疾病的发病机制。
研究内容与主要结论
为解决上述问题，研究人员从基因表达综合数据库（Gene Expression Omnibus, GEO）下载了三个与IDD相关的转录组数据集（GSE167199、GSE176205和GSE186542），共包含21个样本（9个健康髓核组织和12个退变髓核组织）。通过整合分析，研究人员旨在发现与氧化应激相关的核心基因，探讨其与免疫细胞浸润模式的关系，并构建上游调控网络，最终通过人体和大鼠模型进行实验验证。
本研究的主要结论包括：成功鉴定了CD36、SESN3和HMOX1作为IDD中关键的氧化应激相关生物标志物；发现这些基因的表达与免疫微环境的重塑，特别是与静息肥大细胞、初始CD4+T细胞和巨噬细胞的浸润密切相关；并通过人体组织样本和大鼠模型的实验验证，确认了这些基因在退变椎间盘中的上调表达。这些发现揭示了氧化应激与免疫调控在IDD发病机制中的复杂相互作用，为未来的诊断标志物开发和治疗策略设计提供了新的方向。
关键技术方法
研究人员首先将三个GEO数据集合并，并使用ComBat算法校正批次效应。通过limma包筛选差异表达基因，并与分子特征数据库（Molecular Signatures Database, MSigDB）中的氧化应激相关基因取交集，获得候选基因。采用共识聚类分析基于氧化应激相关基因表达谱划分分子亚型，并通过t分布随机邻域嵌入（t-distributed stochastic neighbor embedding, t-SNE）验证分型效果。运用加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis, WGCNA）鉴定与氧化应激亚型高度相关的基因模块，并结合模块隶属度和基因显著性筛选枢纽基因。通过CIBERSORT反卷积算法评估22种免疫细胞的浸润比例，并进行Spearman相关性分析探索枢纽基因与免疫细胞丰度的关联。最后，通过定量聚合酶链反应（quantitative polymerase chain reaction, qPCR）和免疫组织化学（immunohistochemistry, IHC）在人髓核组织样本（9例，按Pfirrmann分级分为对照组和IDD组）中验证，并在大鼠尾椎针刺退变模型中进行免疫荧光（immunofluorescence, IF）和qPCR验证。
研究结果
3.1 鉴定IDD中氧化应激相关的差异表达基因
通过整合三个GEO数据集并进行批次校正后，主成分分析显示对照组与IDD组之间存在明显的转录差异。差异表达分析发现了大量显著上调和下调的基因。基因本体论（Gene Ontology, GO）富集分析表明这些差异表达基因与氧化应激相关生物学过程密切相关，京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路分析则显示其在MAPK等信号通路中显著富集。
3.2 共识聚类揭示不同的氧化应激相关分子亚型
基于氧化应激相关基因的表达谱，共识聚类分析将21个样本稳定地划分为两个分子亚型（Cluster 1和Cluster 2），t-SNE分析证实了这两个亚型的清晰分离。对两个亚型进行差异表达分析，并与预设的氧化应激基因列表交叉比对，获得了氧化应激相关差异表达基因。GO富集分析显示这些基因主要与氧化还原稳态、脂肪酸代谢和细胞对氧化应激的反应有关。其中，CD36、SESN3和HMOX1被确定为排名靠前的枢纽基因，在两个亚型间表达水平差异显著。
3.3 WGCNA鉴定与氧化应激亚型相关的关键模块
通过WGCNA分析，研究人员构建了基因共表达网络。在与临床性状的相关性分析中，含有2539个基因的绿松石色模块与氧化应激相关亚型表型的相关性最强。结合基因显著性大于0.6和模块隶属度大于0.7的标准，在该模块中共筛选出44个枢纽基因。
3.4 筛选枢纽生物标志物并构建调控网络
将WGCNA得到的44个枢纽基因与聚类分析得到的主要氧化应激相关差异表达基因取交集，最终确定了三个共同的候选生物标志物：CD36、SESN3和HMOX1。单基因GSEA分析揭示了这些基因与多种通路的关联，例如SESN3与病毒蛋白与细胞因子受体相互作用通路呈负相关；CD36与细胞周期、泛醌及其他萜类-醌的生物合成等通路呈负相关；HMOX1与细胞周期、果糖和甘露糖代谢等通路呈负相关。研究人员还构建了由mRNA-miRNA-lncRNA相互作用的ceRNA调控网络和TF-mRNA调控网络，揭示了这些枢纽基因可能受到复杂的转录后和转录水平的调控。
3.5 免疫浸润特征及其与枢纽基因的相关性
CIBERSORT分析显示，不同分子亚型和疾病状态下的免疫浸润模式存在差异。例如，Cluster 2中静息肥大细胞的丰度显著高于Cluster 1，而初始CD4+T细胞的丰度则较低。与对照组相比，IDD组的静息记忆CD4+T细胞比例降低。Spearman相关性分析表明，HMOX1的表达与静息肥大细胞、M0巨噬细胞和活化NK细胞呈正相关，与初始CD4+T细胞和初始B细胞呈负相关；CD36仅与静息肥大细胞呈正相关；SESN3与静息肥大细胞呈正相关，与初始CD4+T细胞呈负相关。
3.6 在人体组织中对枢纽基因进行实验验证
通过qPCR和IHC检测人髓核组织样本，结果显示与对照组（Pfirrmann II级）相比，IDD组（Pfirrmann III-V级）中CD36、SESN3和HMOX1的mRNA和蛋白表达水平均显著升高。从不同退变阶段组织中分离出的髓核细胞进行的免疫荧光分析也证实了这一结果。
3.7 在大鼠椎间盘退变模型中CD36、SESN3和HMOX1表达上调
在大鼠尾椎针刺诱导的IDD模型中，组织学染色（H&E、Masson、SOFG）证实了退变模型的成功建立。免疫荧光和qPCR检测均表明，与假手术组相比，退变组大鼠椎间盘组织中CD36、SESN3和HMOX1的表达水平显著上调。
讨论与结论
讨论总结： 本研究通过整合多个转录组数据集并应用系统的生物信息学分析方法，成功鉴定了CD36、SESN3和HMOX1作为IDD中连接氧化应激与免疫调节的关键基因。研究发现IDD中存在两种与氧化应激相关的分子亚型，支持了IDD的生物学异质性。CD36作为清道夫受体，参与脂肪酸摄取、脂质处理、炎症反应和氧化损伤，其与静息肥大细胞的正相关暗示了代谢压力与免疫激活之间的潜在联系。SESN3属于Sestrin家族，具有抗氧化保护和代谢调节功能，其与初始CD4+T细胞的负相关表明氧化还原敏感机制可能参与了局部免疫调节。HMOX1编码血红素加氧酶1，是一种经典的应激诱导的保护性酶，其与多种免疫细胞类型的广泛相关性提示它在桥接氧化应激反应与免疫适应方面发挥重要作用。研究还构建了预测性的ceRNA和TF-mRNA调控网络，为理解这些基因的上游调控机制提供了基础。本研究的优势在于使用了多个独立数据集以减少偏倚，并结合了计算筛选与人体组织和大鼠模型的实验验证。局限性包括：公共数据库汇总样本量不大；免疫浸润是通过计算方法估算而非直接测量；对枢纽基因的功能验证仅限于表达层面，缺乏直接的机制研究；ceRNA等调控网络仍为预测性结果，需实验证实。
研究结论： 该研究利用生物信息学方法鉴定并提出了CD36、SESN3和HMOX1作为IDD中潜在的氧化应激相关生物标志物，并强调了它们与免疫微环境重塑之间的密切关系。这些结果有助于更好地理解椎间盘退变过程中氧化过程与免疫应答之间的关系，并可能有助于构思未来具有抗氧化和免疫调节作用的治疗策略。
