单细胞荟萃分析揭示胰腺内皮细胞转录组特征并确定NKX2-3在PLVAP表达中的关键作用

Single-Cell Meta-Analysis Uncovers the Pancreatic Endothelial Cell Transcriptomic Signature and Reveals a Key Role for NKX2-3 in PLVAP Expression

作者信息Safwat T Khan, Neha Ahuja, Sonia Taïb, Shabana Vohra, Ondine Cleaver, Sara S Nunes
PMID39445426
发布时间2024-12
DOI10.1161/ATVBAHA.124.321781

实验完整度

包含多组学数据整合、独立数据验证、空间转录组学、原位RNAscope/免疫荧光以及体外转染功能验证,证据层级多元。

主要模型

人脐静脉内皮细胞(HUVEC) 胚胎期18.5天小鼠胰腺组织 人类胎儿、新生儿和成人胰腺单核/单细胞RNA测序数据 人类胰岛和外分泌内皮细胞批量RNA测序样本

重点核对

转基因NKX2-3-GFP在HUVEC中的转染效率和表达时间点(4小时约80%,24小时约20%) NKX2-3转染后HUVEC中PLVAP、SPARCL1、SMAD7、EDN1表达上调,CLEC14A、HSPG2表达下调 E18.5小鼠胰岛与内分泌区域PLVAP免疫荧光强度、SPARC和NKX2-3 RNAscope探针结合面积的差异 空间转录组中胰岛与外分泌区域内皮亚群基因特征的Fisher z检验得分

摘要

背景:胰腺血管系统表现出组织特异性的生理和功能适应,以支持β细胞的快速胰岛素反应。然而,消化酶使得表征胰腺内皮细胞(EC)变得困难,导致对胰腺EC特化的理解不足。方法:挖掘可用的单核/单细胞RNA测序数据集,以鉴定胰腺EC富集的标志基因,并构建人类胰腺EC的综合图谱。我们使用独立的单核/单细胞RNA测序数据、分离的EC的批量RNA测序数据、空间转录组学数据、免疫荧光和所选标志物的RNAScope验证了这些发现。通过基因转染在人脐静脉内皮细胞(HUVEC)中表达NK2同源盒3(NKX2-3)转录因子(TF),并通过RT-qPCR评估胰腺EC富集标志基因的表达。结果:我们定义了一个跨物种和发育阶段保守的胰腺EC富集基因特征,包括参与ECM(细胞外基质)组成(COL15A1和COL4A1)、通透性和屏障功能(PLVAP、EHD4、CAVIN3、HSPG2、ROBO4、HEG1和CLEC14A)以及关键信号通路(S1P [鞘氨醇-1-磷酸]、TGF-β [转化生长因子-β]、RHO/RAC GTPase [三磷酸鸟苷酶]、PI3K/AKT [磷脂酰肌醇3-激酶/蛋白激酶B]和PDGF [血小板衍生生长因子])的基因。整合图谱揭示了胰腺内的血管层级。我们鉴定并验证了一个特化的胰岛毛细血管亚群,其特征在于参与通透性(PLVAP和EHD4)、免疫调节(FABP5、HLA-C和B2M)、ECM组成(SPARC和SPARCL1)、IGF(胰岛素样生长因子)信号传导(IGFBP7)和膜转运(SLCO2A1、SLC2A3和CD320)的基因。重要的是,我们确定了NKX2-3是胰腺EC中富集的关键TF。DNA结合基序分析发现NKX2-3基序存在于约40%的标志基因中。在HUVEC中诱导NKX2-3促进了胰岛毛细血管EC富集基因PLVAP和SPARCL1的表达。结论:我们定义了一个经过验证的胰腺EC转录组特征,并揭示了其组织内转录组异质性。我们证明NKX2-3作用于PLVAP的上游,并提供了一个可供社区进一步探索的单细胞在线资源:https://vasconcelos.shinyapps.io/pancreatic_endothelial/。

实验结论

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

研究问题
胰腺内皮细胞的转录组特征及其组织内异质性如何,以及哪些转录因子调控这些特征?
核心机制
NKX2-3作为转录因子直接或间接调控胰腺内皮细胞特异性基因的表达,包括上调PLVAP、SPARCL1、SMAD7和EDN1,下调CLEC14A和HSPG2,从而参与内皮细胞表型的特化。
主要证据
通过多组学数据整合和独立验证,发现胰腺EC特征基因在胰岛毛细血管EC中显著富集;体外转染实验证明NKX2-3诱导改变多个特征基因表达。
研究意义
研究为理解胰腺内皮细胞在葡萄糖稳态和糖尿病病理生理中的作用提供了资源,并可能有助于未来在体外模型中生成胰岛内皮细胞。

研究路径

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

1

识别胰腺内皮细胞特征基因

从多器官单细胞数据中鉴定胰腺内皮细胞特异性富集的基因,排除环境RNA污染。

分析Descartes和Tabula Muris EC图谱,比较胰腺EC与其他组织EC的差异表达,并与多个胰腺数据集中EC富集的基因求交集,得到65个特征基因。

2

验证特征基因的胰腺特异性

在独立的多器官图谱和原位实验中验证特征基因的表达特异性。

在Tabula Sapiens成人人类图谱中评估特征基因的模块得分,并进行聚类;通过RNAscope在E18.5小鼠胰腺组织中检测Nkx2-3、Plcb1、Clec14a和Heg1的表达。

3

构建人类胰腺内皮细胞整合图谱

描绘胰腺内血管层级和内皮细胞亚群。

整合三个数据集的胰腺EC,使用Harmony去除批次效应,聚类得到6个亚群,包括动脉、毛细血管和静脉,并通过轨迹分析推断分化关系。

4

验证胰岛毛细血管亚群

确认胰岛毛细血管EC亚群的特性和空间位置。

分析bulk RNAseq数据(Jonsson等人),比较胰岛和外分泌EC的基因表达;在胎儿胰腺空间转录组数据中评估亚群特征;通过免疫荧光和RNAscope验证PLVAP、SPARC和NKX2-3的表达。

5

验证NKX2-3的调控作用

评估NKX2-3对胰腺内皮细胞特征基因表达的调控作用。

进行DNA结合基序分析,发现NKX2-3基序存在于37%的特征基因中;在HUVEC中瞬时转染NKX2-3-GFP,通过RT-qPCR检测目标基因表达变化。

研究方法

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

数据整合与差异表达分析
基因集富集分析
转录因子基序分析
空间转录组分析
体外功能验证

产品清单

实验环节名称品牌货号
胰岛素抗体Cell Signaling Technologycatalog No. 4590
胰高血糖素抗体Milliporecatalog No. 4030-01F
VE-钙黏蛋白抗体Santa Cruzcatalog No. sc-6458
Pecam1抗体Santa Cruzcatalog No. c-1506
Plvap抗体BD Pharmingencatalog No. 550563
Endomucin抗体Santa Cruzcatalog No. sc-65495
GFP抗体Thermo Fisher Scientificcatalog No. A10262
NKX2-3抗体European Monoclonal Antibodies Network--
NKX2-3-GFP质粒OriGeneRG206667
GFP质粒OriGenePS100010
人脐静脉内皮细胞ScienCellcatalog No. 8000
内皮细胞培养基ScienCellcatalog No. 1001
HUVEC核转染试剂盒LonzaVPB-1002
RNeasy Micro试剂盒QIAgen74004
iScript cDNA合成试剂盒Bio-Rad1708890
SsoAdvanced通用SYBR Green预混液Bio-Rad1725271
罗氏LightCycler 480 II仪器Roche--
RNAscope HiPlex v2检测Advanced Cell Diagnostics--
尼康A1R共聚焦显微镜Nikon--
BioTek Lionheart FX自动显微镜BioTek--

关键环节

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

环节核对要点
数据获取与分析
使用的数据集名称、来源、版本和处理方式(如过滤阈值、标准化方法)需与原文一致。
阅读提示:参见Materials and Methods中Data Acquisition和Seurat Analysis部分
NKX2-3转染实验
HUVEC的代数(2-5代)、转染的质粒量和细胞数量、转染后培养时间(24小时)、转染试剂盒型号、质粒来源
阅读提示:参见Materials and Methods中Transfection部分
RT-qPCR
RNA提取试剂盒、逆转录试剂盒、qPCR预混液、内参基因(ACTB和GAPDH)、基因特异性引物序列(见补充材料)
阅读提示:参见Materials and Methods中RT-qPCR部分和补充材料Major Resources Table
原位验证
小鼠品系(CD1)、胚胎日龄(E18.5)、组织处理步骤(固定、渗透、封闭)、抗体稀释度(1:100或1:200)、RNAscope探针浓度
阅读提示:参见Materials and Methods中Immunofluorescence Staining、In Situ RNAscope HiPlex v2 Assay部分
空间转录组分析
空间转录组数据的来源(GSE197317)、处理流程(SCTransform)、用于计算亚群特征得分的基因列表
阅读提示:参见Materials and Methods中Spatial Transcriptomics Analysis和Results中Spatial Transcriptomics部分