肾脏囊肿-透明细胞肾细胞癌共病分析揭示FOXM1+上皮通过COL6A1-CD44轴重塑SFRP2+成纤维细胞并提示isovitexin干预潜力

《Translational Oncology》：A renal cyst-clear cell renal cell carcinoma comorbidity analysis identifies a FOXM1-positive epithelial COL6A1-CD44 axis that remodels fibroblasts

透明细胞肾细胞癌（clear cell renal cell carcinoma, ccRCC）是最常见且具侵袭性的肾癌亚型，其特征是由VHL（von Hippel-Lindau）缺失驱动的广泛代谢重编程、异常血管生成、免疫逃逸以及显著的肿瘤间和肿瘤内异质性。尽管靶向治疗和免疫检查点阻断（immune checkpoint blockade, ICB）取得进展，但相当比例患者仍难以获得持久疾病控制，这凸显了深入理解控制肿瘤进展的分子程序的必要性。肾囊肿状态具有异质性，包括常见或散发性囊肿、复杂性肾囊肿病变、获得性囊性肾病、遗传性囊性疾病和囊性肾肿瘤；这些实体并不构成单一疾病连续体。因此，研究人员使用广义肾囊肿相关疾病-基因框架来识别与ccRCC共享的分子程序，而不建模特定囊性疾病或直接囊肿到癌症进展。单细胞和空间转录组学（single-cell and spatial transcriptomics）的最新进展从根本上重塑了肾癌研究格局，使该领域从体量预后特征转向机制可解释的肿瘤-基质串扰模型。里程碑式单细胞研究以前所未有分辨率解析了ccRCC细胞异质性，鉴定出恶性上皮亚群、免疫抑制生态位和活化成纤维细胞状态，这些共同塑造肿瘤微环境。互补空间转录组分析进一步证明，当表达相关性由物理细胞邻近、下游调控激活及跨多个组织切片的持续空间模式强化时，配体-受体通讯才最可靠地推断。这些方法学进展为建立多层分析框架创造了机会，其中体量风险建模鉴定临床相关基因程序，单细胞转录组学解析来源细胞状态，空间转录组学验证连接恶性上皮与其基质伙伴的组织水平回路。在本研究中，研究人员实施了整合多层分析流程，从共病基因发现开始并逐步聚焦至可成药的上皮-基质轴。首先，来自整理的疾病-基因资源的肾囊肿相关基因与TCGA-KIRC肿瘤-正常差异表达基因取交集，建立共享共病基因宇宙。第二，评估101个机器学习生存工作流以分离该共享程序的临床相关组分，产生具有稳健跨队列区分能力的随机生存森林（random survival forest, RSF）风险评分。第三，免疫解卷积和药理学敏感性分析鉴定与高风险疾病相关的候选微环境区室。第四，单细胞转录组整合将高风险信号定位至FOXM1⁺恶性上皮细胞；scTenifoldKnk虚拟FOXM1敲除鉴定出COL6A1、CD44和细胞外基质（extracellular matrix, ECM）重塑基因等显著扰动命中。第五，空间转录组学和成纤维细胞亚聚类将FOXM1⁺上皮连接至具有不良预后的CD44^high F5_SFRP2⁺ ECM重塑成纤维细胞状态。最后，通过对接和分子动力学进行药理学优先排序，鉴定isovitexin为候选FOXM1靶向化合物，体外功能实验证实其抑制ccRCC细胞增殖、迁移和侵袭的能力。总之，这些分析将FOXM1⁺上皮和COL6A1-CD44通讯与成纤维细胞重塑联系起来，并鉴定isovitexin为ccRCC中值得进一步评价的候选物。

透明细胞肾细胞癌（clear cell renal cell carcinoma, ccRCC）是最常见且具侵袭性的肾癌亚型，其分子特征包括VHL（von Hippel-Lindau）缺失驱动的代谢重编程、异常血管生成、免疫逃逸及显著的肿瘤间和肿瘤内异质性。尽管靶向治疗和免疫检查点阻断（immune checkpoint blockade, ICB）取得进展，但持久疾病控制仍具挑战。肾囊肿相关疾病具有异质性，并不构成单一疾病连续体。因此，研究人员采用广义肾囊肿相关疾病-基因框架，识别其与ccRCC共享的分子程序，而非建模特定囊性疾病或直接囊肿到癌症进展。单细胞和空间转录组学（single-cell and spatial transcriptomics）的进展使研究从体量预后特征转向机制可解释的肿瘤-基质串扰模型。研究人员实施多层分析流程：从肾囊肿相关基因与TCGA-KIRC差异表达基因交集出发，经101种机器学习生存工作流筛选出随机生存森林（random survival forest, RSF）风险评分，再结合免疫解卷积、药理学敏感性分析、单细胞转录组、scTenifoldKnk虚拟FOXM1敲除、空间转录组、成纤维细胞亚聚类、分子对接和分子动力学，最终以体外功能实验验证候选化合物。该研究得出，FOXM1⁺恶性上皮为高危来源状态，COL6A1-CD44通讯连接FOXM1⁺上皮与CD44^high F5_SFRP2⁺ ECM重塑成纤维细胞，isovitexin可抑制ccRCC细胞增殖、迁移和侵袭，作为候选治疗化合物值得进一步临床前评价。

方法概括：研究人员从DisGeNET、GeneCards和OMIM检索肾囊肿相关基因，并与TCGA-KIRC肿瘤-正常差异表达基因取交集；使用TCGA-KIRC、GSE167573、GSE22541、E-MTAB-1980、JAVELIN_AA和IMmotion151_ATEZO进行生存及免疫治疗验证；整合GSE159115、GSE210038、GSE242299单细胞数据；采用GSE175540及Kalogirou et al.报道的转移性RCC空间免疫治疗应答研究共16个空间样本；构建RSF模型并行免疫解卷积、单细胞Harmony/inferCNV/Scissor/SCENIC/scTenifoldKnk/CellChat、空间RCTD/mistyR/联合生态位/SpaGene/核心距离、成纤维亚聚类/NicheNet/Monocle3/PROGENy、传统中医药（traditional Chinese medicine, TCM）资源分子对接与100 ns分子动力学，以及CCK-8、集落形成、划痕和Transwell实验。

结果部分保留各小标题并浓缩如下。第一，237基因肾脏囊肿-ccRCC共享程序连接上皮发育、肾脏发育和ECM重塑。研究人员通过疾病-基因资源与TCGA-KIRC差异表达基因交集获得237个共享基因，GO富集显示其集中于上皮发育、RNA聚合酶II转录调控、细胞命运决定、肾脏发育、膜极性及含胶原ECM等条目；突变分析显示114个共享基因存在体细胞突变，以错义突变为主，并总结CNV频率；两个方向孟德尔随机化分析支持肾癌遗传易感与相关肾囊肿表型之间的关联。第二，源自共享程序的RSF模型显示稳健预后并对已发表ccRCC特征具有良好基准表现。研究人员将237个基因经单变量Cox筛选至74个，再经101种生存工作流评估后选择RSF，最终67个基因进入跨队列评分；该模型在TCGA-KIRC、E-MTAB-1980、GSE22541和GSE167573中显示区分能力，并在IMmotion151_ATEZO和JAVELIN_AA中区分无进展生存。第三，高RSF风险提名富含Tregs、B细胞、浆细胞和成纤维细胞的免疫-基质微环境。研究人员以TCGA-KIRC队列特异性RSF截断值分层，IOBR整合免疫解卷积显示高风险肿瘤中调节性T细胞（regulatory T cells, Tregs）、浆细胞、B细胞谱系和成纤维细胞增加，ESTIMATE免疫评分升高，同时检查点与炎症激活并存，MHC和抗原呈递特征异质；oncoPredict/GDSC2预测药物反应提示DNA损伤反应、细胞周期、有丝分裂、RTK/血管生成、PI3K/mTOR、凋亡和细胞应激等依赖。第四，FOXM1⁺恶性上皮是高危来源状态，虚拟FOXM1敲除提名COL6A1-CD44通讯。研究人员综合cytoHubba、RSF变量重要性和表达生存分析将FOXM1列为机制枢纽；单细胞整合后，FOXM1⁺肿瘤上皮细胞较FOXM1^?上皮具有更高推断CNV负担，Scissor将体量RSF风险表型投射至单细胞，SCENIC显示Scissor⁺上皮中FOXM1及FOXM1_extended调节子活性增加；scTenifoldKnk虚拟FOXM1敲除在Scissor⁺上皮中鉴定COL6A1、CD44及LTBP1、LUM、WISP1、MMP2、MMP9、COL6A3、ADAM12、COL1A1、LOX、SERPINE1、COL5A1等ECM/重塑基因为显著命中。第五，空间转录组学支持FOXM1⁺上皮-成纤维细胞生态位和COL6A1-CD44耦合。研究人员对16个空间样本行RCTD解卷积、mistyR共定位、联合生态位划分、SpaGene和核心距离分析；联合生态位分析中FOXM1⁺肿瘤上皮与成纤维细胞共同占据生态位7、8、9，并在这些生态位中富集；SpaGene在代表样本中支持COL6A1-CD44空间耦合，核心距离分析显示成纤维细胞RCTD权重和COL6A1-CD44配体-受体评分在FOXM1⁺上皮核心或紧邻环最高并向外下降；CellChat显示FOXM1⁺肿瘤上皮向成纤维细胞存在COL6A1-CD44信号，而FOXM1^?上皮未见相同显著模式。第六，CD44^high F5_SFRP2⁺成纤维细胞代表提名通讯模型中最可能的接收状态。研究人员将成纤维细胞再聚类为F1_ADIRF、F2_LZTS1、F3_ATF3、F4_MYH11、F5_SFRP2、F6_MALAT1和F7_TOP2A；F5_SFRP2⁺成纤维细胞表达SFRP2、EFEMP1、PDGFRA、FAP、CTHRC1、LUM、NTM、THBS2、VCAN和POSTN，呈活化ECM重塑状态；NicheNet支持FOXM1⁺上皮来源COL6A1为高活性配体，F5具高CD44表达和COL6A1-CD44轴评分；F5特征在TCGA-KIRC中与较短总生存相关；Monocle3和ClusterGVis将F5样成纤维细胞置于晚期/终末重塑状态，PROGENy提示缺氧、TGF-beta、NF-kB、TNF-alpha及雄激素/雌激素、p53、MAPK、JAK-STAT和WNT信号偏倚，SCENIC提名GATA6、FOXO1、SOX6、NR2F1和RUNX2调节子。第七，TCM药理学资源和对接/分子动力学优先鉴定isovitexin为FOXM1靶向候选物。研究人员从TCMSP、HERB、TCMIP、HIT等资源筛选FOXM1相关草药活性成分，选择AKBA、honokiol和isovitexin进行分子对接和100 ns分子动力学；对接显示三者均可结合FOXM1，AKBA对接能最低，isovitexin接近，honokiol较弱；PLIP显示isovitexin形成七个氢键，涉及LEU259、LYS260、ARG297、SER306和TRP308；分子动力学中isovitexin平均氢键数最高，故被优先作为FOXM1靶向候选物。第八，Isovitexin在体外抑制ccRCC细胞增殖、迁移和侵袭。研究人员在786O和OSRC2细胞中用0、20和40 μmol/L isovitexin处理，CCK-8显示浓度和时间依赖性活力抑制，集落形成实验显示克隆能力下降，划痕实验显示迁移减弱，Transwell迁移和侵袭实验显示迁移及侵袭细胞数减少，且40 μmol/L抑制最强。

讨论部分总结：该研究提出一个肾脏囊肿相关共享程序，连接FOXM1⁺上皮、COL6A1-CD44通讯和CD44^high F5_SFRP2⁺成纤维细胞重塑。其核心前提是，共病驱动基因发现与多组学及空间分辨率整合，可揭示传统肿瘤-only分析可能忽略的可操作肿瘤-基质回路。高危组混合免疫-基质状态与既往KIRC炎症相关预后特征及免疫细胞组成、免疫检查点表达研究一致。分析提示FOXM1可能除增殖相关转录因子作用外，还参与ccRCC上皮-基质通讯；单细胞将FOXM1定位于高CNV负担恶性上皮，SCENIC显示FOXM1调节子在Scissor⁺高危上皮中富集，scTenifoldKnk虚拟敲除提名COL6A1和CD44等ECM重塑基因，支持可测试的FOXM1-COL6A1-CD44通讯假说，但未建立FOXM1对COL6A1或CD44的直接调控。空间转录组提供组织水平验证：16个ccRCC切片中FOXM1⁺上皮与成纤维细胞反复共定位并共享生态位，COL6A1-CD44信号富集于FOXM1⁺上皮-成纤维细胞界面，径向梯度支持局部旁分泌通讯模型。F5_SFRP2⁺成纤维细胞是最可能接收状态，具高CD44表达、不良预后影响和肿瘤相关肌成纤维细胞重塑、缺氧、TGF-beta及NF-kB偏倚，伪时间分析将其置于终末ECM重塑状态，SCENIC提示FOXO1和RUNX2等候选机制连接CD44 engagement与成纤维细胞活化。药理学分析通过对接和100 ns分子动力学优先选择isovitexin，体外实验证实其浓度依赖性抑制增殖、克隆能力、迁移和侵袭。限制包括：肾囊肿组分依赖整理的疾病-基因资源和群体遗传工具，因此识别共享分子程序而非证明直接囊肿到癌症进展；空间共定位支持物理邻近和配体-受体耦合，但正式功能依赖需通过上皮特异性FOXM1敲除等模型强化因果推断；isovitexin验证在已建细胞系中完成，需扩展至患者来源类器官或体内模型评价生理相关背景中的疗效。结论翻译：一个肾脏囊肿相关框架鉴定出237个共享基因和一个RSF风险状态。单细胞和空间分析将FOXM1⁺上皮与成纤维细胞富集的生态位联系起来，并支持一个COL6A1-CD44通讯轴。F5_SFRP2⁺成纤维细胞是最可能的接收状态，并显示ECM重塑表型。Isovitexin抑制ccRCC细胞中的恶性表型，并作为候选治疗化合物值得进一步临床前评价。