《Molecular Genetics and Genomics》：Single-cell transcriptomic profiling combined with Mendelian randomization illuminates molecular drivers of bladder cancer
编辑推荐：
本文针对膀胱癌细胞高度异质性、分子驱动因素难以明确的问题，整合单细胞RNA测序与孟德尔随机化技术，系统性筛选出与膀胱癌发生发展因果相关的关键基因。研究发现，B细胞亚群在疾病中变化最显著，并通过MR分析锁定了六个潜在驱动基因。其中，功能验证证实ARHGEF18是风险基因，而YPEL5具有保护作用。该研究为理解膀胱癌的细胞类型特异性机制提供了新见解，并展示了一种可扩展的、连接遗传风险与功能机制的整合性研究策略。
膀胱癌是全球范围内最常见的泌尿系统恶性肿瘤之一，其发生发展受到遗传风险、细胞异质性和肿瘤-免疫相互作用的共同影响。尽管在手术、化疗和免疫治疗方面取得了一定进展，但对于晚期，尤其是肌层浸润性和转移性膀胱癌的管理，仍然是一个重大挑战。这些挑战的根源，很大程度上在于肿瘤细胞的高度异质性，以及对其起始和进展的潜在分子过程认知有限。传统的研究方法，如全基因组关联研究，虽然鉴定出了一些与疾病风险相关的基因位点，但难以确定其确切的因果基因及其在特定细胞类型中的功能。同样，高通量测序技术发现的体细胞突变，也常常忽略了非编码基因组元件的影响以及肿瘤微环境内复杂的细胞组成。这种知识与技术上的鸿沟，凸显了阐明膀胱癌中细胞类型特异性调控机制的迫切性。近年来，单细胞RNA测序技术与孟德尔随机化分析方法的结合，为解决这一难题提供了强大的工具。前者能以无与伦比的分辨率描绘肿瘤的细胞图谱，后者则能为遗传变异与基因表达乃至疾病结局之间的因果关系提供推断依据。将两者结合，有望从细胞层面精准定位那些既有遗传基础又有功能影响的“分子驱动者”。
本研究正是运用了这一整合策略，旨在系统地识别在特定细胞环境中与膀胱癌相关的关键基因，以加深对膀胱癌生物学的理解，并为开发更精准的诊断和治疗策略奠定基础。该论文已发表在《Molecular Genetics and Genomics》期刊上。
为了开展这项研究，研究人员主要运用了以下几项关键技术：首先，利用公共数据库获取了膀胱癌的单细胞转录组数据（GSE222315，含13个样本）以及多个批量转录组数据集用于验证。其次，通过Seurat和Harmony等生物信息学工具对单细胞数据进行严格的质控、标准化、整合、降维和细胞亚群注释。接着，基于B细胞亚群中筛选出的差异表达基因，利用来自eQTLGen联盟的血液表达数量性状位点数据和来自IEU OpenGWAS数据库的膀胱癌全基因组关联研究汇总数据，进行了多方法的孟德尔随机化分析，以推断基因表达与疾病风险的因果关系。此外，研究还综合运用了免疫浸润分析、基因集富集分析、转录因子调控网络构建、细胞伪时间轨迹分析等多种生物信息学方法。最后，通过免疫组化、实时荧光定量PCR、蛋白免疫印迹以及包括CCK-8、划痕愈合、Transwell侵袭和流式细胞术在内的体外细胞功能实验，对候选基因进行了多层次验证。其中，用于免疫组化验证的临床样本来自南充中心医院在2021年3月至2025年3月期间存档的根治性膀胱切除术标本。
研究结果
膀胱癌中的细胞亚群与功能特征
通过对13个膀胱癌相关组织样本的单细胞转录组数据分析，研究者最终获得了72,664个高质量细胞，并将其注释为11种细胞类型，包括上皮细胞、CD4+T细胞、调节性T细胞、CD8+T细胞、B细胞、浆细胞、肥大细胞、巨噬细胞、成纤维细胞、内皮细胞和一个未分类群。通过评估细胞丰度变化和转录谱变化，研究者发现B细胞对膀胱癌病理的贡献最为突出。对B细胞特异性标志基因的通路分析显示，这些基因显著富集于免疫反应上调、利什曼病相关通路和细胞质蛋白质合成等生物学事件中。
多方法验证的孟德尔随机化分析鉴定膀胱癌驱动候选基因
为了从B细胞的标志基因中进一步筛选关键基因，研究者利用孟德尔随机化方法分析了六个标志基因的表达与膀胱癌风险之间的因果关系。分析结果强调了六个与膀胱癌相关的基因：ARHGEF18、HLA-DRB5、ISG20、NCF1、RPL13和YPEL5。其中，ARHGEF18、HLA-DRB5、ISG20、NCF1和RPL13与膀胱癌风险升高相关，而YPEL5则显示出保护性关联。留一法敏感性分析证实了这些因果关联的稳健性。因此，这六个基因被确定为后续研究的关键靶点。
肿瘤微环境中的基因免疫相关性
对肿瘤微环境中免疫浸润模式的分析显示，与对照组相比，疾病样本中M1型巨噬细胞等免疫亚群的浸润显著更高。相关性分析进一步揭示了关键基因与特定免疫细胞亚群之间的关联。例如，ARHGEF18与活化的肥大细胞和静息的自然杀伤细胞呈正相关，YPEL5与嗜酸性粒细胞呈强正相关。这些基因还与一系列免疫调节剂（如趋化因子、免疫刺激分子）密切相关，表明它们可能对免疫微环境的组成具有关键影响。
基因富集分析及其与膀胱癌调控和转录因子的关联
基因集富集分析揭示了关键基因参与的信号通路。例如，ARHGEF18主要富集于核糖体功能和氧化磷酸化通路；HLA-DRB5、ISG20和NCF1与NF-κB、Toll样受体和趋化因子信号通路相关；YPEL5则参与cGMP-PKG和NF-κB信号网络。基因集变异分析也得出了类似的富集趋势。此外，研究者还构建了这些关键基因的转录因子调控网络，发现了可能调控它们表达的潜在上游转录因子。
关键基因在单细胞数据中的作用
在单细胞转录组数据中，六个关键基因的表达展现出细胞类型特异性分布。研究者特别聚焦于B细胞，利用伪时间轨迹分析揭示了B细胞在分化过程中的动态变化，并展示了关键基因表达沿发育路径的变化。此外，通过AUCell方法评估发现，这些关键基因在特定的免疫代谢通路中表现出不同的活性，例如YPEL5在MYC靶点v1和氧化磷酸化通路中活性较高，而HLA-DRB5在干扰素γ反应和移植物排斥通路中活性显著。
基因表达差异及其与临床病理特征的关系
在多个独立的转录组数据集和TCGA-BLCA队列中，验证了关键基因的表达差异。一致的趋势显示，ARHGEF18和ISG20在肿瘤组织中表达上调，而YPEL5表达下调。免疫组化实验在38例膀胱癌患者的肿瘤组织和配对癌旁正常组织中也证实了这些表达趋势。分子亚型分析提示，ARHGEF18高表达组中管腔乳头状亚型居多，而YPEL5低表达组中基底鳞状亚型更常见。
ARHGEF18和YPEL5在膀胱癌细胞中的体外分析**
在细胞水平上，与正常尿路上皮细胞相比，膀胱癌细胞系T24和5637中ARHGEF18的表达显著升高，而YPEL5的表达降低。功能实验表明，敲低ARHGEF18或过表达YPEL5，均能导致T24和5637细胞的增殖、迁移和侵袭能力下降，并增加细胞凋亡。这直接证实了ARHGEF18的促癌作用和YPEL5的抑癌作用。
研究结论与讨论
本研究通过整合单细胞转录组学、孟德尔随机化及实验验证，成功识别出ARHGEF18和YPEL5作为膀胱癌中具有遗传和转录支持的关键调控因子。这一整合策略的优势在于，它将遗传风险位点与细胞分辨率的转录组数据相叠加，不仅提供了统计上的稳健性，还赋予了生物学背景，有助于过滤高通量分析中常见的假阳性信号，并加强了生物学解释。
YPEL5作为一个新兴的肿瘤抑制因子，在本研究和既往报道中均显示其下调与多种癌症进展相关，并可能通过调控细胞周期和凋亡发挥作用。本研究还发现其表达与嗜酸性粒细胞浸润正相关，暗示了其潜在的免疫调节角色。而ARHGEF18作为一种Rho鸟嘌呤核苷酸交换因子，其过表达与肿瘤细胞侵袭、迁移以及不良预后相关。本研究不仅支持了其在细胞骨架重塑和运动中的经典功能，其与免疫细胞（如活化肥大细胞）的相关性也提示它可能参与塑造免疫抑制性肿瘤微环境。
这项研究也存在一些局限性，例如用于单细胞分析的初始样本量有限，孟德尔随机化使用的是血液eQTL数据可能无法完全反映组织特异性调控，以及临床样本规模限制了深入的生存分析和亚型分层分析。然而，它所展示的整合遗传因果推断与单细胞分辨率的研究框架，为在膀胱癌这类高度异质性的肿瘤中，超越描述性关联、实现基于生物学的关键驱动基因优先排序，提供了一个可扩展且强有力的策略。该工作加深了人们对膀胱癌细胞类型特异性机制的理解，并为未来开发针对ARHGEF18和YPEL5通路的新型治疗干预措施提供了潜在靶点。
