整合批量转录组与单细胞分析揭示ADPKD进展中的细胞外基质重塑与肾小管代谢衰竭

《Biology》：Integrated Bulk Transcriptomic and Single-Cell Analyses Reveal Extracellular-Matrix Remodeling and Tubular Metabolic Failure in ADPKD Progression

　　

： 背景：常染色体显性多囊肾病（autosomal dominant polycystic kidney disease，ADPKD）的进展涉及囊肿扩张、间质重塑及肾功能进行性丧失，但连接囊肿负荷与组织重塑的细胞类型解析程序尚未被完整阐明。方法：研究人员将按囊肿大小分层的人类批量转录组数据集GSE7869与时间分辨的Pkd1单细胞数据集GSE268494进行整合分析。在小囊肿、中囊肿和大囊肿组中鉴定差异表达基因，通过Mfuzz时间聚类建模，依据蛋白质-蛋白质相互作用拓扑结构和功能富集进行优先级排序，并进一步通过免疫反卷积、单细胞映射、AUCell通路评分及CellPhoneDB配体-受体推断进行评估。结果：研究人员鉴定出3842个共享差异表达基因及两个占主导的单向变化模块。进行性上调模块富集了细胞外基质组织、黏附及纤维化相关基因，包括COL1A1、COL5A2、FBN1、BGN、LAMC1、HSPG2、MMP2、LOX和ITGA5。进行性下调模块则突出肾小管转运、质子转运及线粒体/代谢相关基因，包括SLC12A3、AQP2、ATP5A1、ATP6V0C、ATP6V1B1、NDUFV1、REN和AGT。单细胞分析将基质模块主要定位至成纤维细胞/基质细胞区室，将转运/代谢模块定位至节段特异性肾小管上皮细胞群体。免疫反卷积和配体-受体推断进一步提示巨噬细胞相关的重塑以及TGFB1、CD47、CD200、PECAM1和TNF超家族相关的免疫-实质细胞通讯。结论：ADPKD进展与协调的ECM/免疫重塑及肾小管代谢衰竭相关，为未来的空间和功能验证提供了候选标志物及通讯轴。

论文解读

一、研究背景与科学问题

常染色体显性多囊肾病（autosomal dominant polycystic kidney disease，ADPKD）是最常见的遗传性囊性肾病，也是导致肾衰竭的主要单基因病因。多数病例由PKD1或PKD2基因的致病性变异引起，这两个基因编码参与初级纤毛信号、钙稳态、上皮极性和肾小管结构维持的多囊蛋白复合体。临床上，ADPKD表现为双侧肾脏进行性囊肿扩张、肾脏增大、高血压、疼痛、感染、肾结石及肾功能逐渐丧失。尽管基因型是预后的主要决定因素，但具有相似遗传背景的患者在囊肿生长速率、间质纤维化和肾功能下降方面可表现出显著差异，提示继发性组织程序在很大程度上塑造了疾病轨迹。

当前的治疗管理在临床上重要但生物学上尚不完整。托伐普坦（Tolvaptan）可在部分患者中延缓肾脏总体积扩张和估算肾小球滤过率下降，但现有疗法并不能完全阻止进展至肾衰竭，治疗选择仍主要依赖临床风险分层而非分子分期。因此，亟需定义能够反映囊肿负荷、间质重塑和肾小管功能损害的疾病相关分子模块，以支持生物标志物发现或机制导向的干预策略。

经典的ADPKD上皮模型强调钙和环磷酸腺苷信号改变、囊肿上皮增殖、平面细胞极性缺陷及跨上皮液体分泌。近年来的单细胞和多组学研究通过鉴定囊肿衬覆上皮状态、修复失败肾小管程序、促炎成纤维细胞群体、血管重塑及染色质可及性改变，拓展了这一框架。这些观察支持一种多细胞视角：肾小管上皮应激、间质重塑、血管改变和免疫激活随时间相互强化，而非作为孤立的病理区室发挥作用。

免疫和细胞外基质（extracellular matrix，ECM）重塑与ADPKD进展尤其相关。巨噬细胞极化、炎症细胞因子信号、补体激活、T细胞参与及免疫检查点通路已被认为参与囊肿生长和纤维化。ECM动力学，包括胶原沉积、基底膜重塑、基质交联和金属蛋白酶活性，可能主动改变组织硬度、上皮极性、整合素信号和机械转导，而非仅仅是终末期瘢痕形成的被动反映。与此同时，代谢重编程和线粒体功能障碍已成为连接上皮应激与囊肿扩张、转运活性受损及适应不良修复易感性的候选机制。

批量转录组研究能够捕获疾病分期相关信号，但往往难以区分主导信号究竟来自上皮细胞状态改变、间质激活、免疫浸润还是细胞组成改变。相反，单细胞数据集提供细胞分辨率，但可能受模型系统、采样深度以及将实验时间点直接对齐人类囊肿负荷的困难所限制。因此，将按囊肿大小分层的人类转录组数据与时间分辨的单细胞参考整合，为从细胞类型特异性机制中分离共享进展程序提供了实用策略。

总体而言，存在两个关键空白：第一，连接人类ADPKD囊肿负荷、间质重塑和肾小管功能障碍的可解释转录模块仍然缺乏；第二，人类ADPKD中与囊肿大小相关的转录改变究竟反映共享的进展程序还是区室特异性反应，仍不清楚。本研究旨在通过整合按囊肿大小分层的人类批量转录组数据与Pkd1时间过程单细胞数据集，重建ADPKD中分期相关的分子程序。

二、研究设计与技术方法

研究人员从GEO数据库下载人类ADPKD微阵列数据集GSE7869，选取正常对照（n=3）、小囊肿（<1 mL，n=5）、中囊肿（10–25 mL，n=5）和大囊肿（>50 mL，n=3）共16个样本进行分期分辨的转录组分析。原始探针水平数据经背景校正和分位数归一化处理。差异表达分析使用limma包，以校正后p<0.05且|log2FC|>1为阈值。共享差异表达基因通过Mfuzz软聚类建模，设置六个聚类以捕获主要的上调、下调和非单调表达趋势。进行性上调和下调聚类中的基因提交至STRING数据库构建蛋白质-蛋白质相互作用（protein-protein interaction，PPI）网络，并通过Cytoscape的cytoHubba插件以13种算法进行枢纽基因优先级排序。功能富集分析使用clusterProfiler进行基因本体（Gene Ontology，GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集。免疫浸润分析使用xCell进行免疫和基质细胞丰度推断。单细胞验证使用Pkd1时间过程数据集GSE268494，包含P66、P100和P130三个时间点共16个样本（每个时间点均含对照和Pkd1组）。单细胞分析在R中使用Seurat进行，经SoupX去除环境RNA污染，Harmony校正批次效应，DoubletFinder去除双细胞。候选基因映射至单细胞聚类并通过点图展示细胞类型特异性表达模式。AUCell用于基于MSigDB基因集计算通路评分。细胞间通讯使用CellPhoneDB v5.0.1推断配体-受体相互作用。

三、主要研究结果

3.1 囊肿大小转录组揭示分期相关的基因失调

按囊肿大小分层的转录组分析揭示了人类ADPKD肾脏中进行性、分期依赖的转录失调。小囊肿组鉴定出2594个上调基因和2757个下调基因；中囊肿组鉴定出2364个上调基因和2348个下调基因；大囊肿组鉴定出3028个上调基因和3393个下调基因。三个囊肿大小比较的交叉分析鉴定出3842个共享差异表达基因。Mfuzz软聚类将共享差异表达基因分为六个时间聚类，其中聚类3代表随囊肿增大进行性上调的基因，聚类5代表进行性下调的基因，这两个单调聚类被优先作为核心转录程序进行下游机制分析。

3.2 网络和富集分析鉴定ECM和肾小管稳态模块

聚类3的PPI网络包含高度连接的ECM、黏附和纤维化相关基因，包括COL1A1、COL5A2、FBN1、BGN、LAMC1、HSPG2、MMP2、LOX和ITGA5。功能富集分析显示上调基因富集于细胞外基质组织、细胞外结构组织、含胶原细胞外基质、细胞-底物黏附、糖胺聚糖结合、Wnt信号及基质重塑相关条目。相比之下，聚类5网络富集于上皮转运、线粒体功能和肾单位节段身份相关基因，高优先级基因包括SLC12A3、REN、NR1H4、NDUFV1、CYP3A4、ATP6V1B1、ATP6V0C、ATP5A1、AQP2和AGT。富集分析突出小分子分解代谢、有机酸代谢、顶侧和基底外侧膜区室、离子跨膜转运、质子转运及线粒体/代谢通路。cytoHubba多算法排序进一步验证了这一功能二分法：聚类3枢纽基因反复优先排序ECM和黏附基因，而聚类5枢纽基因优先排序肾小管转运和代谢基因。

3.3 免疫反卷积将候选模块与微环境重塑联系起来

免疫反卷积显示预测的免疫微环境在对照、小囊肿、中囊肿和大囊肿组间存在差异，巨噬细胞、肥大细胞和淋巴细胞特征表现出分期依赖的变化。相关性分析显示，聚类3基因尤其是MMP2、LOX、LAMC1、ITGA5和HSPG2与巨噬细胞和肥大细胞相关信号呈正相关，包括M2样巨噬细胞富集模式。聚类5基因则表现出不同的免疫相关模式，转运和代谢基因如AQP2、SLC12A3、ATP5A1和ATP6V1B1的下降可能与分化上皮功能降低和细胞组成改变有关。

3.4 单细胞映射解析批量衍生模块的细胞背景

Pkd1时间过程单细胞数据集鉴定出主要的肾上皮、基质、血管和免疫区室，包括近端肾小管节段、亨勒袢群体、远曲小管、连接小管、集合管主细胞和闰细胞、足细胞、成纤维细胞、内皮细胞、B细胞、T细胞和髓系细胞。候选基因映射显示ECM重塑基因主要富集于成纤维细胞/基质区室，而肾小管转运/代谢基因表现出肾单位节段限制性上皮表达：Slc12a3富集于远曲小管细胞，Aqp2富集于集合管主细胞，Atp6v1b1和Atp6v0c富集于闰细胞相关区室。分期分层单细胞分析显示人类批量囊肿大小谱与小鼠Pkd1疾病模型之间的表达趋势一致。

3.5 AUCell分析揭示细胞类型特异性通路重编程

AUCell通路评分显示通路重塑是分期和细胞类型依赖性的。显著改变的通路数量随疾病分期增加，并在后期阶段的皮区室中突出。模式分类将通路程序分为早期、中期、晚期、瞬时、进行性、持续或双向类别。热图、模块汇总、轨迹图和提琴/密度图显示纤维化、缺氧、炎症、凋亡、上皮-间质转化和干扰素相关通路激活，同时氧化磷酸化和脂肪酸代谢程序在选定肾小管区室中减弱。模块相关性分析进一步将批量ECM和肾小管模块与细胞类型特异性通路活性联系起来。

3.6 配体-受体推断提示免疫-实质细胞通讯改变

CellPhoneDB配体-受体推断在Pkd1缺陷肾脏中鉴定出免疫细胞与肾实质细胞之间重连的细胞间信号。B细胞相关相互作用包括CD200–CD200R1、LTB–LTBR、TGFB1受体复合物、TGFB1–整合素-αVβ6、PECAM1–CD38及胆固醇相关受体轴。T细胞相关通讯涉及CD226–NECTIN2、CD96–NECTIN1、CD40LG–整合素、TNFSF10–受体、PECAM1–CD38及TGFB1相关相互作用。髓系细胞相互作用包括TGFB1受体复合物、SIRPA/SIRPG–CD47、TNFSF13B–TNFRSF13B、CD93–IFNGR1及TYROBP–CD44。促纤维化的TGFB1信号可能桥接免疫细胞与成纤维细胞和基质激活；CD47、CD200和PECAM1相互作用可能介导免疫调节和黏附依赖性信号；TNF超家族信号可能维持炎症反应。

四、讨论与结论总结

讨论部分指出，本研究支持ADPKD进展的分期化、多细胞模型，其中囊肿扩大与ECM重塑、免疫激活和肾小管转运/代谢身份丧失相耦合。进行性上调模块以ECM组织、含胶原基质、黏附和基质重塑条目为主，ECM重塑可能不仅是被动的晚期瘢痕形成后果，还可通过改变组织硬度、配体可用性和囊周结构，修饰整合素信号、上皮极性、机械转导、成纤维细胞激活和巨噬细胞募集。ECM相关基因定位于成纤维细胞和基质区室，提示聚类3程序反映成纤维细胞/基质细胞内基质相关基因表达增加。进行性下调模块则突出分化肾小管上皮功能的丧失，转运和代谢基因的下降提示囊肿进展伴随上皮功能去分化。免疫结果增加了第三层：xCell反卷积将ECM枢纽基因与巨噬细胞和肥大细胞相关信号联系起来，配体-受体推断提示通过TGFB1相关、CD47相关、TNF超家族、CD200–CD200R1、PECAM1–CD38和整合素相关轴的通讯改变。AUCell结果进一步表明通路重塑是分期和细胞类型特异性的，ADPKD进展并非统一的肾脏广泛应激反应，而是不同细胞群体随时间进入不同的适应不良状态。从转化角度看，候选基因应被视为实验验证的优先标志物而非直接治疗靶点。

研究结论翻译如下：总之，对按囊肿大小分层的人类转录组和细胞类型分辨单细胞数据的整合分析，在ADPKD中鉴定出两个相反的进展程序：一个以COL1A1、COL5A2、FBN1、BGN、LAMC1、HSPG2、MMP2、LOX和ITGA5为中心的进行性上调ECM/基质模块，以及一个以SLC12A3、AQP2、ATP5A1、ATP6V0C、ATP6V1B1、NDUFV1、REN和AGT为中心的进行性下调肾小管转运/代谢模块。免疫反卷积、AUCell评分和配体-受体推断进一步将这些模块与巨噬细胞相关重塑、纤维化/缺氧/炎症通路激活以及涉及TGFB1、CD47、CD200–CD200R1、PECAM1–CD38、整合素相关信号和TNF超家族相互作用的免疫-实质细胞通讯轴联系起来。这些发现为ADPKD进展提供了分期分辨的分子框架，并提名了特定的ECM、肾小管代谢和免疫-基质候选分子，供未来空间验证和功能测试。