Inflammation-20240417:基于 hdWGCNA 和单细胞生物信息学分析的免疫细胞浸润调控网络在椎间盘退化中的作用【hdWGCNA应用实例】

  • 免疫浸润在椎间盘退化 (IDD) 中起着至关重要的作用。在这项研究中,我们通过单细胞生物信息学分析探讨了 IDD 的免疫微环境。本研究整合了三个单细胞数据集。根据特征基因将髓核细胞 (NPCs) 分为亚组,通过伪时间轨迹分析分析每个亚组在 IDD 过程中的作用。使用 hdWGCNA 获得枢纽基因,并通过批量数据集和伪时间序列进一步鉴定。枢纽基因的表达定义了与免疫浸润相关的 NPC,并探讨了这些 NPC 与免疫细胞之间的相互作用。NPC 分为四个亚组:储备 NPC、HCL-NPC、反应 NPC 和支持 NPC,它们分别主导 IDD 的四个过程:非、轻度、中度和重度变性。SPP1 和 ICAM1 被鉴定为髓核免疫浸润枢纽基因。巨噬细胞和中幼粒细胞分别通过 SPP1-CD44 通路和 ICAM1-ITGB2 配体-受体通路在 SPP1-ICAM 双通 NPC 组中发挥促炎作用。同时,双方 NPC 都通过 ANXA1-FPR1 通路寻求中性粒细胞的自助炎症缓解。对分化和免疫浸润景观的系统分析有助于了解 IDD 的整体发展过程。我们的数据表明,SPP1 和 ICAM1 可能是治疗 IDD 炎症浸润的新靶点。

直接上结果

1. 本研究的流程图。

image.png

2. 单细胞亚群

集成单细胞数据集的 NPC 亚群的特征。a 所有细胞的主成分分析 (PCA)。b 批量移除后的单细胞分布图。c 来自三个单细胞数据集的 17 个髓核样品的 RNA 质量的小提琴图。d 所有 NPC 的 T-SNE 图,颜色表示不同的亚组。e NPC 亚群基因表达特征热图。f 四个 NPC 亚组的 GO 富集分析的气泡图。

3. 伪时间分析阐明了 IDD 进展的特点

本研究观察到四种类型 NPC 的比例在 IDD 变性过程中有所不同。后备 NPC 在前期占据了主要份额,然后又分为两个分支,其中一个分支仍然以后备 NPC 为主(超过一半),而另一个主要分化为 HCL NPC。随后,它到达了关键的分叉点(红色圆圈),其中一个分支由绝大多数 HCL NPC 组成,而另一个分支继续退化并从 HCL 分化为反应 NPC。大多数响应 NPC 继续退化为辅助 NPC(图 D)。3a-c)。在组织样本水平上进一步检查细胞类型的分布,得出的结果与伪时间分析相似但不完全一致(图 D)。3d)。典型的 IDD 标记基因(如 COL1A1 和 MMP2)的表达呈具有伪时间序列的波浪状,而包括 ACAN 和 COL2A1 在内的健康 NPC 标记基因的表达逐渐降低(图 D)。3e)。这些结果证实了该序列与 IDD 的发展方向一致。

图 3: 伪时间序列揭示了 NPC 亚群分化与 IDD 进展之间的关系。a NPC 的伪时间分析,红色圆圈表示时间轨迹中的关键分叉点。b NPC 的四个亚组的比例分布随序列的进展而变化。c 伪时间序列的每个时期亚群细胞的百分比图。d 显示了不同 IDD 水平的组织水平 NPC 亚组的百分比。e 典型 IDD 相关基因在伪时间序列中的表达变化。f 序列过程中基因表达热图和主要富集 KEGG 通路的聚类。

4. IDD Hub 模块和 Hub 基因的鉴定【这儿用的hdWGCNA】

响应 NPC 在时间轨迹的关键分叉点(图 1 中的红色圆圈)分支出来。3a),这是 HCL NPC 分化为支持 NPC 的过渡细胞亚群。因此,进一步评估了响应 NPC 中的特征模块。使用和声去除批处理效果后,确定了五个特征模块(软阈值 = 4)(图 D)。4a 和 b)。模块 1 (M1) 的 KEGG 通路富集分析显示 TNF-α 信号通路NF-κB 信号通路的富集(图 1)。4c)。响应 NPC 中 M1 的和谐模块特征基因 (hME) 值与其他细胞类型显著不同(图 D)。4d),并且 M1 与其他三个模块没有显著相关性(图 D)。4e),表明 M1 在反应 NPC 的炎症反应中是必不可少的。因此,从 M1 中鉴定 IDD 免疫浸润中心基因是合理的,因为ICAM1 和 NFKB1 等炎症标志物都存在于 M1 基因网络中(图 D)。4f)。

图4: 基于hdWCGNA的关键模块识别。 a 响应 NPC 的树状图。 b hdWCGNA 识别的响应 NPC 的 Hub 模块。 c 基于 KEGG 的富集路径气泡图。 d 每个NPC子组中关键模块的表达。 e 模块之间的关联性。 f 来自response-M1模块的代表性基因的网络。

与 M1 模块中心基因、反应 NPC 特征标记基因和 2499 个免疫相关基因交叉后,鉴定出 7 个共表达基因(GBP2、ICAM1、NFKB1、PTX3、SAA1、SPP1 和 TNFAIP3)。通过伪时间分析验证了这 7 个基因,其中 4 个(GBP2、ICAM1、NFKB1 和 SPP1)显着表达。 GBP2、SPP1、ICAM1和NFKB1的表达从中期IDD开始增加,这与反应NPC比例的增加相一致(图5a)。此外,这七个中心基因的表达水平在 GSE70362 RNA-seq 数据集中得到了进一步验证,该数据集包括 24 个不同 IDD 程度的 NP 组织样本(图5b)。仅SPP1和ICAM1的表达模式符合IDD中期增加、晚期减少的特点。随后,我们检测了不同变性程度的原代人髓核细胞中SPP1和ICAM1 mRNA的相对表达水平,发现它们的趋势与数据集一致(图5c)。 Ⅲ级SPP1表达量显着高于Ⅱ、IV级(P < 0.05,P < 0.01); III级ICAM1表达量显着高于II、IV、V级(P < 0.05,P < 0.001,P < 0.01)。与软基质对照组相比,硬细胞外基质促炎模型也可以增加这两个基因的表达(图5d)。此外,WB结果显示NP组织中SPP1和ICAM1的蛋白表达水平与mRNA水平一致(图5e)。因此,最终选择SPP1和ICAM1作为进一步研究的中心基因。

图5: 鼻咽癌免疫相关枢纽基因的筛选与验证。 a 伪时间分析显示hub基因的表达趋势。 b 批量数据集 GSE70362 验证了不同 IDD 水平下 hub 基因的表达水平。 c 以GAPDH为内参,人原代髓核细胞中SPP1和ICAM1的相对mRNA表达水平。 d, e 基于 GSEA 的双上 NPC 与双下 NPC 的富集途径和交互网络。 *P < 0.05,**P < 0.01,***P < 0.001,****P < 0.0001。

5. SPP1和ICAM1双向NPC主要与巨噬细胞通讯

细胞间通讯在决定各种细胞分化命运中起着至关重要的作用。对免疫细胞和 NPC 之间细胞间通讯的分析表明,秘密信号传导和细胞间接触途径是主要的通讯模式(图6a)。对免疫细胞与NPC之间细胞间通讯的分析表明,巨噬细胞(单核细胞)是细胞通讯的关键枢纽,与免疫细胞中的中性粒细胞相互作用密切且最为密切(图6b)。在 NPC 中,双向组细胞主要通过 SPP1 途径进行自分泌通讯。此外,巨噬细胞通过SPP1-CD44配体-受体对通过SPP1途径与双向NPC进行通信(图6c-e)。骨髓细胞通过 ICAM1-ITGB2 通路与 Both-up NPC 进行通信(图6f)。相反,当双上NPC为源细胞时,它们主要通过ICAM途径与骨髓细胞相互作用,其中ICAM1 − (ITGAX + ITGB2)配体-受体对是最重要的(图6 g)。当靶向中性粒细胞时,both-up NPC 主要依赖膜联蛋白途径,其中 ANXA1-FPR1 配体-受体对最为重要(图6h)。

图6: ICAM1 和 SPP1 参与免疫细胞和 NPC 之间的通讯。 a NPC 和免疫细胞之间的细胞-细胞接触网络。 b NPC与免疫细胞相互作用的微环境图。 c–e SPP1信号通路网络。 f ICAM信号通路网络。 g ITGB2信号通路网络。 h ANNEXIN 信号通路网络。
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 194,491评论 5 459
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 81,856评论 2 371
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 141,745评论 0 319
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,196评论 1 263
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 61,073评论 4 355
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,112评论 1 272
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,531评论 3 381
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,215评论 0 253
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,485评论 1 290
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,578评论 2 309
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,356评论 1 326
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,215评论 3 312
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,583评论 3 299
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 28,898评论 0 17
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,174评论 1 250
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,497评论 2 341
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,697评论 2 335

推荐阅读更多精彩内容