《Translational Oncology》：Shared inflammatory-metabolic transcriptional effects in tubulointerstitial fibrosis and clear cell renal cell carcinoma
编辑推荐：
慢性肾脏病（chronic kidney disease, CKD）是全球重大健康负担，肾小管间质纤维化（tubulointerstitial fibrosis）是进行性肾损伤常见的组织病理学终点。肾小管上皮细胞高度依赖氧化代谢，脂肪酸氧化（fatty-aci
慢性肾脏病（chronic kidney disease, CKD）是全球重大健康负担，肾小管间质纤维化（tubulointerstitial fibrosis）是进行性肾损伤常见的组织病理学终点。肾小管上皮细胞高度依赖氧化代谢，脂肪酸氧化（fatty-acid oxidation, FAO）受损已在实验和转录水平上与肾纤维化相关联。与此同时，透明细胞肾细胞癌（clear cell renal cell carcinoma, ccRCC）以广泛的代谢重塑为特征，包括氧感知和能量代谢的改变。大规模分子谱分析已确立代谢重编程为ccRCC中与肿瘤生物学相关的反复出现的特征。这些观察为探讨纤维化与ccRCC是否在组织病理学和临床行为截然不同的情况下仍共享可重复的转录效应提供了依据。
炎症、线粒体功能障碍、脂质处理和氧化还原调节是合理的汇聚点。铁死亡（ferroptosis）是一种由磷脂过氧化驱动的铁依赖性调节性细胞死亡，ACSL4通过塑造可供氧化的多不饱和磷脂库而成为铁死亡敏感性的重要决定因素。肾损伤模型也证明了GPX4依赖性脂质过氧化控制和铁死亡性肾小管损伤的作用。然而，铁死亡相关基因或脂质-氧化还原程序的表达不应等同于铁死亡执行的直接证据；需要功能性检测才能确立该机制。
跨疾病转录组比较因组织组成、检测平台、归一化和物种差异而复杂化。因此，本分析使用独立的人类证据层而非原始尺度合并。TCGA-KIRC作为人类ccRCC bulk组织参考，GSE12682提供人类显微切割肾小管间质纤维化队列，GSE328575仅用于ccRCC中描述性细胞类型定位。该设计避免了正式跨疾病基因水平比较中的小鼠-人类直系同源转换，也避免了RNA-seq与微阵列效应幅度的直接校准。
主要目标因此是在预设显著性和效应量标准下，识别在人类ccRCC和人类肾小管间质纤维化中均独立与疾病相关的基因和通路水平方向。研究人员进一步在TCGA-KIRC中检验了探索性临床关联，并使用ccRCC单细胞图谱确定哪些细胞群体对选定的转录信号有贡献。所得框架旨在用于证据三角验证和候选优先排序，而非因果推断。
本研究强调证据层的分析分离、队列内标准化和明确的解释边界。通过聚焦独立一致的效应，研究人员试图区分稳健的跨疾病信号与可能源于平台缩放、细胞组成或疾病特异性生物学的关联。
总之，研究人员对人类ccRCC和肾小管间质纤维化转录组进行了独立的跨疾病再分析，辅以ccRCC中的描述性单细胞定位以及HK-2、NRK-52E和786-O细胞模型中的RT-qPCR验证。计算分析识别出一组有限但统计上非随机的基因水平一致效应和一致的程序水平方向，而细胞实验提供了具有明确引物序列和RT-qPCR质量指标的物种解析mRNA验证。
论文解读
一、研究背景与问题
慢性肾脏病（CKD）已成为全球重大公共卫生负担，而肾小管间质纤维化是进行性肾损伤最常见的组织病理学终点。肾小管上皮细胞高度依赖氧化代谢，脂肪酸氧化（FAO）受损已被实验和转录水平证据与肾纤维化联系起来。与此同时，透明细胞肾细胞癌（ccRCC）以广泛的代谢重塑为特征，大规模分子谱分析已确立代谢重编程为其反复出现的肿瘤生物学特征。尽管纤维化与ccRCC在组织病理学和临床行为上截然不同，两者是否共享可重复的转录效应尚不清楚。炎症、线粒体功能障碍、脂质处理和氧化还原调节被认为是可能的汇聚点。铁死亡作为一种铁依赖性调节性细胞死亡，其敏感性由ACSL4决定，因为该酶塑造了可供氧化的多不饱和磷脂库。然而，既往研究多局限于单一疾病或单一数据类型，缺乏跨疾病、跨平台的系统比较。跨疾病转录组比较还面临组织组成、检测平台、归一化和物种差异等复杂问题。因此，有必要采用独立证据层设计，在严格控制的标准下识别两种疾病间共享的转录改变。
二、研究内容与主要结论
研究人员对TCGA-KIRC（533例原发肿瘤和72例实体组织正常样本）、GSE12682人类显微切割肾小管间质纤维化队列（13例健康与23例纤维化样本）和GSE328575 ccRCC单细胞图谱（65,770个细胞）进行了独立分析。在11,805个对齐基因中，64个基因在两种疾病中均满足预设的FDR<0.05、|log2FC|≥0.5且方向一致的共享效应标准，其中51个共享上调、13个共享下调。程序水平上，炎症评分在两种疾病中均升高（RCC +0.64，纤维化 +0.28），脂肪酸氧化评分均降低（RCC ?0.86，纤维化 ?0.13）。S100A9是原始炎症-代谢面板中唯一满足完整共享效应规则的基因，在ccRCC单细胞图谱中主要定位于单核/巨噬细胞群体。RT-qPCR实验证实，TGF-β1处理的HK-2和NRK-52E细胞以及786-O肿瘤细胞均表现为HADHA、PPARA、PPARGC1A下调而ACSL4上调的一致代谢转变。该研究的重要意义在于提供了一个保守的跨疾病汇聚模型：纤维化和ccRCC独立地表现出炎症转录增强和脂肪酸氧化减弱，共享基因富集于免疫相关信号，但两者仍保持分子层面的显著差异。
三、主要关键技术方法
本研究采用多层级独立证据设计。TCGA-KIRC作为ccRCC bulk组织参考，GSE12682作为独立人类肾小管间质纤维化队列，GSE328575仅用于描述性细胞类型定位，三个数据集分别分析、不进行原始尺度合并。基因水平整合采用人类-人类比对，无需小鼠-人类直系同源转换。GSE12682使用limma经验贝叶斯回归进行差异表达估计，以性别和纤维化状态为协变量。跨疾病共享效应要求两队列各自独立满足FDR<0.05和|log2FC|≥0.5且方向一致，候选排序采用队列内效应绝对值排名而非直接y=x校准。RT-qPCR验证使用HK-2、NRK-52E和786-O三种细胞系，以GAPDH/Gapdh为内参，采用2^-ΔΔCt法计算相对表达量，使用双侧单样本t检验和Holm校正进行统计检验。
四、研究结果
独立纤维化队列、样本结构与差异表达。 GSE12682分析保留36例显微切割肾小管间质样本，性别调整后limma火山图显示332个基因达到FDR<0.05，其中205个同时满足|log2FC|≥0.5的效应量标准。
RCC肿瘤-正常参考效应。 TCGA-KIRC中15,699个基因达到FDR<0.05，6,399个基因满足效应量标准。程序水平上，肿瘤相对正常表现为炎症升高（+0.64）、脂肪酸氧化降低（?0.86）、脂质氧化还原升高（+0.18）、糖酵解/缺氧升高（+0.60）。领先效应包括HK2（+3.34）、SCD（+3.36）、FABP1（?3.32）等。
TCGA-KIRC探索性临床背景。 530例有病理分期的肿瘤中，联合炎症-代谢评分与分期呈弱负相关（rho=?0.122）。中位切分的log-rank检验P=0.0551，联合评分单变量Cox风险比每标准差为0.833。脂肪酸氧化评分HR=0.623，炎症评分HR=1.109。HADHA表达随分期递减（rho=?0.161）。
RCC单细胞图谱背景。 GSE328575包含65,770个细胞，恶性细胞占62.2%。描述性定位显示HIF1A富集于巨噬细胞，ACSL4富集于肥大细胞，CCL2富集于成纤维细胞，IL6富集于恶性细胞。库水平上联合评分与恶性细胞比例呈正相关（rho=0.722）。
RCC单细胞程序与细胞状态背景。 炎症评分在单核细胞（1.202）、巨噬细胞（0.673）和单核/巨噬细胞（0.563）中最高，脂质氧化还原在巨噬细胞中最高（0.450），糖酵解/缺氧在恶性细胞中最高（0.154）。
次要原发-转移RCC分析。 8例原发与8例转移库比较中，脂肪酸氧化在转移中降低（效应=?0.123），NFE2L2为最大位点效应（?0.387），但无任何精选靶基因通过FDR<0.05。
独立跨疾病基因水平一致性。 11,805个对齐基因中6,690个（56.7%）方向一致，64个通过完整共享效应规则。全局效应相关性为rho=0.224（P=1.05e-133），强非随机。领先共享上调基因包括APOC1、C1QB、CD3D、MS4A6A、MYO1F、LGALS9、IL10RA，共享下调基因包括MYH14、USP2、PLA2R1。S100A9的RCC log2FC=+1.11，纤维化log2FC=+1.16。
整合证据与解释边界。 64个共享基因分为51个共享上调和13个共享下调。S100A9在ccRCC图谱中最高表达于Mono、Mono.Mac和Mac细胞群体。
体外RT-qPCR验证代谢候选基因。 HK-2细胞中TGF-β1处理后HADHA降至0.56倍、PPARA降至0.61倍、PPARGC1A降至0.48倍、ACSL4升至1.82倍。NRK-52E细胞方向一致。786-O相对HK-2表现为HADHA 0.42倍、PPARA 0.49倍、PPARGC1A 0.34倍、ACSL4 2.33倍。
五、讨论总结与结论翻译
讨论部分强调，本研究的核心发现是肾小管间质纤维化和ccRCC共享统计上非随机但定量适度的疾病相关转录改变成分。64个共享基因的有限重叠比广泛的单一共享致病程序的声明在生物学上更合理，强调两种疾病尽管存在选定的共同反应，但在分子层面仍显著不同。通路水平上最一致的模式是炎症升高伴FAO降低。S100A9作为唯一满足完整共享效应规则的精选靶基因，其信号至少部分具有组成敏感性。铁死亡相关解释应保持适当边界：本研究未直接测量脂质过氧化、GPX4蛋白丰度、铁死亡诱导剂敏感性或药理学挽救，因此观察到的模式应描述为铁死亡相关或铁死亡允许性背景而非铁死亡执行的证据。研究也未建立直接的HADHA-ACSL4调控轴，PPARα/PGC-1α协调调控模型仍属假说。探索性临床分析提供背景而非预后生物标志物的验证。GSE328575单细胞图谱存在推断局限性，细胞水平分布并非独立的患者水平重复。未来需要匹配组织区室、空间转录组学、功能性代谢测量和候选调控因子的直接扰动实验来确定哪些共享效应具有机制可操作性。
结论：对TCGA-KIRC和人类GSE12682肾小管间质纤维化队列的独立分析，在预设FDR、效应量和方向标准下识别出64个具有一致疾病相关效应的基因。程序水平上，两种疾病均表现为炎症升高和脂肪酸氧化降低，S100A9作为原始精选面板中独立共享的成员，在ccRCC单细胞图谱中主要定位于单核/巨噬细胞区室。HK-2、NRK-52E和786-O模型中的物种解析RT-qPCR实验显示HADHA、PPARA和PPARGC1A降低而ACSL4升高，为转录组分析中观察到的代谢方向提供了实验性mRNA水平支持。综合发现支持肾小管间质纤维化和ccRCC之间存在可重复但有限的炎症-代谢重叠，但不建立因果性跨疾病通路、直接铁死亡执行、优化的生存生物标志物或治疗效力。
