单细胞与批量测序数据的整合分析描绘P2ry12在缺血再灌注损伤后小胶质细胞中的表达与功能

Integrative Analysis of Single-Cell and Bulk Sequencing Data Depicting the Expression and Function of P2ry12 in Microglia Post Ischemia-Reperfusion Injury.

作者信息Chenglong Wang, Li Peng, Yuan Wang, Ying Xue, Tianyi Chen, Yanyan Ji, Yishan Li, Yong Zhao, Shanshan Yu
PMID37047745
发布时间2023-04-05
DOI10.3390/ijms24076772

实验完整度

整合scRNA-seq和bulk-seq数据鉴定MSPGs,结合轨迹分析、细胞互作分析,并通过qPCR、Western blot、免疫荧光等体内实验验证P2ry12表达与功能。

主要模型

小鼠tMCAO模型(GSE167593、GSE174574) 大鼠tMCAO模型(Sprague-Dawley大鼠) 小胶质细胞(离体微胶质细胞排序数据)

重点核对

tMCAO模型:成年雄性Sprague-Dawley大鼠(250-280g),缺血1小时再灌注 P2ry12拮抗剂PSB0739:剂量0.3 mg/kg,脑室内注射,tMCAO前2小时 qPCR检测P2ry12 mRNA:时间点为再灌注后1小时、6小时、24小时 Western blot检测P2ry12蛋白:再灌注后24小时,以β-actin为对照 免疫荧光双标:Iba1(1:300,Abcam)和CD34(1:100,ABclonal)

摘要

P2ry12是小胶质细胞标志基因。最近,越来越多的证据表明其表达水平可随不同的中枢神经系统疾病而变化,并可影响小胶质细胞的功能,如极化、可塑性和迁移。然而,缺血再灌注损伤(IRI)期间小胶质细胞中P2ry12的表达和功能仍不清楚。在此,我们开发了一种计算方法,利用与IRI相关的测序数据获得小胶质细胞特异性P2ry12基因(MSPGs)。我们评估了IRI期间MSPGs综合表达水平的变化,并将其与P2ry12的表达进行比较以确定相似性。随后,通过生物信息学方法利用MSPGs探索小胶质细胞中P2ry12的功能。此外,还进行了多项动物实验以确认结果的可靠性。观察到P2ry12的表达在损伤后24小时内逐渐下降。作为响应,P2ry12表达降低的小胶质细胞中一个受体编码基因(Flt1)和三个配体编码基因(Nampt、Igf1和Cxcl2)的表达增加。此外,双标记免疫荧光染色显示,抑制P2ry12在IRI期间阻断了小胶质细胞向血管的迁移。总的来说,我们采用计算与实验相结合的方法成功探索了IRI期间小胶质细胞中P2ry12的表达和功能。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
缺血再灌注损伤(IRI)期间,小胶质细胞中P2ry12的表达和功能如何变化?
核心机制
低表达P2ry12的小胶质细胞激活,通过上调配体(Nampt、Cxcl2、Igf1)和受体(Flt1)介导细胞间通讯,促进小胶质细胞向血管迁移。
主要证据
整合scRNA-seq和bulk-seq数据鉴定MSPGs,MSPGs评分显示tMCAO组显著降低;qPCR和Western blot验证P2ry12 mRNA和蛋白水平在IRI后下降;免疫荧光显示PSB0739抑制P2ry12后,血管壁小胶质细胞密度降低。
研究意义
研究为理解IRI中小胶质细胞功能提供新见解,并提示P2ry12抑制剂的给药时机可能影响IRI结局;所开发的计算方法可用于单基因表达与功能分析。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

scRNA-seq数据获取与预处理

获取IRI小鼠大脑皮层单细胞转录组数据,并进行质量控制、去噪和数据整合。

从GEO获取两个scRNA-seq数据集(GSE167593、GSE174574),按标准进行细胞和基因筛选,去除线粒体污染,并使用ALRA算法进行技术丢失数据的插补,使用Harmony算法进行数据整合。

2

细胞类型注释

识别并注释各细胞群,确定小胶质细胞群。

使用SingleR、scmap-cell和SCINA三种自动注释工具,综合多数标签确定细胞类型。

3

鉴定MSPGs

通过整合scRNA-seq和bulk-seq数据,获得与小胶质细胞中P2ry12表达强相关的基因集。

计算scRNA-seq数据中各基因与P2ry12的Pearson相关系数,使用meta分析整合两个数据集;通过最小超几何测试确定最优相关系数阈值(0.74),最终得到88个MSPGs。

4

MSPGs评分验证

验证MSPGs表达能否替代P2ry12表达。

使用AUCell和AddModuleScore计算每个细胞的MSPGs诱导评分,比较不同细胞群及sham和tMCAO组的评分差异。

5

P2ry12表达变化验证

验证IRI后P2ry12表达是否下降。

在多个bulk-seq数据集中计算MSPGs平均表达,并在大鼠tMCAO模型中通过qPCR和Western blot检测P2ry12 mRNA和蛋白水平。

6

轨迹分析与差异基因分析

探索P2ry12表达与小胶质细胞激活状态的关系。

使用Monocle R包基于MSPGs构建小胶质细胞轨迹,将细胞分为三个状态,并进行差异基因表达分析。

7

细胞间通讯分析

识别低P2ry12小胶质细胞与其他细胞类型的相互作用。

使用CellChat R包进行细胞间通讯分析,比较低/高P2ry12小胶质细胞中配体和受体的表达差异,并在bulk-seq数据中验证。

8

免疫荧光验证P2ry12调控小胶质细胞迁移

验证P2ry12是否调控IRI中小胶质细胞向血管的迁移。

建立大鼠tMCAO模型,给予PSB0739抑制P2ry12,通过双标记免疫荧光染色检测血管壁附着的小胶质细胞密度。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

单细胞RNA测序数据分析
批量RNA测序数据分析
基因相关性分析
基因组评分
细胞间相互作用分析
基因本体富集分析

产品清单

实验环节名称品牌货号
尼龙线栓Cinontech--
激光多普勒血流仪PerimedPeriflux System 5000
PSB0739GLPBIO--
2,3,5-三苯基氯化四氮唑Sigma--
甲苯胺蓝----
Iba1抗体Abcamab178846
CD34抗体ABclonalA10796
488偶联山羊抗大鼠IgGAbbkineA23220
cy3偶联山羊抗大鼠IgGAbbkineA22220
含DAPI的Vectashield封片剂----
RNAiso PlusTaKaRa Biotechnology--
PrimeScript逆转录试剂盒TaKaRa Biotechnology--
SYBR Premix Ex Taq IITaKaRa Biotechnology--
P2ry12引物Sangon Biotech--
CFX96 Touch实时荧光定量PCR检测系统Bio-Rad--
PVDF膜Millipore--
P2ry12抗体Signalway Antibody31255
β-肌动蛋白抗体ABclonalAC004
成像密度计Bio-Rad--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
tMCAO模型建立
动物品系及体重、麻醉剂量、缺血时间、再灌注时间点
阅读提示:Methods 4.10
PSB0739给药
剂量、给药途径、给药时间
阅读提示:Methods 4.10
qPCR
引物序列、所用试剂盒及仪器型号
阅读提示:Methods 4.13
Western blot
蛋白上样量、抗体稀释比例、抗体品牌及货号
阅读提示:Methods 4.14
免疫荧光
抗体及稀释比例、封片剂、成像显微镜
阅读提示:Methods 4.12