《HUMAN MUTATION》：Integrative Mendelian Randomization and Single-Cell Pseudotime Analysis Reveal DKK3 as a PI3K–AKT-Modulated Driver of Esophageal Squamous Cell Carcinoma
编辑推荐：
食管鳞癌（ESCC）的恶性进展机制尚不清晰，Dickkopf-3 (DKK3)的作用存在争议。为此，研究者整合了cis-eQTL孟德尔随机化、TCGA差异表达、单细胞RNA测序（scRNA-seq）与体内外功能实验，系统探究了DKK3的因果作用。研究将DKK3鉴定为ESCC的一个因果性致癌调控因子，其通过激活PI3K–AKT信号通路驱动肿瘤增殖、迁移与克隆形成。该发现支持DKK3作为ESCC的潜在生物标志物和治疗靶点。
食管鳞状细胞癌（Esophageal Squamous Cell Carcinoma, ESCC）是全球范围内，尤其是在东亚地区，一种极具侵袭性的恶性肿瘤。尽管在手术、放化疗等方面不断取得进展，但患者的预后依然不理想，生存率有待提高。这背后一个关键瓶颈在于，我们对驱动ESCC发生和发展的深层分子机制，特别是那些具有因果关系的核心驱动因子，仍然缺乏清晰、完整的认识。传统的关联性研究，包括基于批量组织测序的分析，往往难以区分因果与伴随现象，也无法精确解析肿瘤内部复杂的细胞异质性和动态演变过程。为了破解这个难题，并发现新的、有效的治疗靶点，寻找并验证具有因果作用的致癌基因已成为领域内的迫切需求。近期，一项发表在《HUMAN MUTATION》的研究，通过整合遗传流行病学、多组学数据和前沿功能实验，为我们提供了一个从因果推断到功能验证的完整研究范例，并鉴定出了一个关键基因——Dickkopf-3 (DKK3)。
为了系统性探索与ESCC风险存在因果关联的基因，并阐明其在肿瘤发展中的动态角色，研究者采用了一套整合性的创新研究策略。其主要技术方法包括：1）利用eQTLGen联盟（31,684人）的顺式表达数量性状位点（cis-eQTL）和大型GWAS汇总数据进行孟德尔随机化（Mendelian Randomization, MR）分析，以筛选与ESCC风险存在潜在因果关联的基因。2）整合癌症基因组图谱（The Cancer Genome Atlas, TCGA）数据库中ESCC肿瘤与正常组织的转录组数据，进行差异表达基因（DEGs）分析。3）从基因表达综合数据库（GEO）获取并处理食管鳞癌的单细胞RNA测序（scRNA-seq）数据，使用Seurat工具包进行细胞聚类、注释和拟时序（Pseudotime）分析，以解析细胞异质性和发育轨迹。4）在体外利用食管鳞癌细胞系（TE-1和ECA109）进行功能验证，包括蛋白质印迹（Western Blot）、CCK-8细胞增殖实验、划痕愈合实验和克隆形成实验。5）在体内通过构建裸鼠皮下移植瘤模型，评估目标基因对肿瘤生长的作用。
3.1. 差异表达分析与孟德尔随机化
研究首先对TCGA数据库的7例正常组织和184例ESCC肿瘤组织进行了差异表达分析，共鉴定出1186个显著下调基因和5947个显著上调基因。接着，利用基于cis-eQTL的孟德尔随机化分析，筛选出231个与ESCC风险存在因果关联的基因。将两组基因取交集后，最终聚焦于22个候选基因，其中3个基因在肿瘤中下调且与ESCC风险呈负相关，19个基因在肿瘤中上调且与ESCC风险呈正相关。这一步从大量基因中初步锁定了可能与ESCC存在因果关系的核心基因集。
3.2. 单细胞分析中的质量控制与降维
为了在细胞层面深入理解ESCC的复杂性，研究者对来自GEO数据库的scRNA-seq数据进行了严格的质量控制和降维聚类分析。通过主成分分析（PCA）和均匀流形近似与投影（UMAP）算法，成功将细胞分为B细胞、T细胞、自然杀伤细胞、成纤维细胞、上皮细胞、肥大细胞、线粒体相关细胞和巨噬细胞等八大类，并可视化比较了ESCC细胞与正常细胞在UMAP空间中的分布差异。这为后续在特定细胞类型中探究候选基因的表达模式奠定了基础。
3.3. MR基因评分与可视化
研究者进一步将目光聚焦于与肿瘤发生最密切相关的上皮细胞。通过分析上皮细胞在ESCC与正常组织中的差异表达基因，并与之前筛选出的22个MR候选基因进行比对，最终确定了三个在ESCC上皮细胞中高表达且与MR结果一致的基因：DKK3、NINJ2和SAPCD2。通过计算这三个基因的MR评分并定位到单细胞层面，发现聚类0和12是MR评分最高的上皮细胞亚群，提示这些细胞可能富含与ESCC发生发展密切相关的分子特征。
3.4. 拟时序分析与特征基因鉴定
为了解这三个候选基因在肿瘤演进过程中的动态变化，研究者对上皮细胞进行了拟时序分析。结果显示，在DKK3、NINJ2和SAPCD2这三个基因中，只有DKK3的表达水平随着上皮细胞拟时序（即假定的恶性演进轨迹）的推进而呈现持续、显著的升高趋势。这一关键发现表明，DKK3的表达变化与上皮细胞的恶性转化进程紧密相关，因此被选为进一步机制验证的核心目标。
3.5. 特征基因的敲低与验证
研究者首先在蛋白质水平验证了DKK3在ESCC细胞系（TE-1, ECA109）中高表达，而在正常食管上皮细胞系（HET-1A）中低表达。随后，利用小干扰RNA（siRNA）成功敲低了TE-1和ECA109细胞中的DKK3表达。功能实验表明，敲低DKK3能显著抑制ESCC细胞的增殖能力、迁移能力以及克隆形成能力。这些体外实验结果首次在功能层面证实了DKK3是促进ESCC细胞恶性表型的关键因子。
3.6. 特征基因在动物模型中的作用
为了在更接近人体的环境中验证DKK3的功能，研究者在裸鼠体内构建了移植瘤模型。将对照（si-NC）或DKK3敲低（si-DKK3-1）的ECA109细胞皮下注射到裸鼠体内。活体成像、肿瘤体积测量和重量分析均一致表明，与对照组相比，DKK3敲低显著抑制了移植瘤的生长。这项体内实验有力地证实了DKK3是驱动ESCC肿瘤体内生长所必需的。
3.7. 富集通路分析与验证
为了揭示DKK3促进ESCC进展的下游分子机制，研究者对TCGA数据中与DKK3表达相关的基因进行了基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析。分析结果强烈提示，DKK3可能通过调节PI3K-AKT信号通路发挥作用。为了验证这一生物信息学预测，研究者检测了DKK3敲低后ESCC细胞中关键信号蛋白的磷酸化水平。Western Blot结果显示，敲低DKK3显著降低了磷酸化PI3K（p-PI3K）和磷酸化AKT（p-AKT）的蛋白水平。这从分子机制上证实，DKK3的确是通过激活PI3K-AKT信号通路来发挥其促癌作用的。
3.8. DKK3的孟德尔随机化研究结果
在研究的最后，研究者回归到最初的因果推断问题，对DKK3与ESCC的因果关系进行了详尽的孟德尔随机化验证。通过森林图、散点图、漏斗图和留一法敏感性分析等多种方法，均得到一致且稳健的结果。这些分析共同强有力地支持了“DKK3基因表达水平升高是导致ESCC风险增加的原因”这一因果假设，排除了反向因果关系和明显的发表偏倚，为整个研究的因果推断链条提供了坚实的人口遗传学证据。
结论与讨论
本研究的结论清晰而有力：通过整合孟德尔随机化、单细胞转录组学和功能实验，首次将DKK3鉴定为ESCC的一个因果性致癌驱动因子。DKK3在ESCC中高表达，其表达量随着上皮细胞恶性进展而增加；在功能上，DKK3是维持ESCC细胞增殖、迁移、克隆形成和体内肿瘤生长所必需的；在机制上，DKK3至少部分通过激活PI3K-AKT信号通路来发挥其促癌作用。这项研究的意义重大。首先，在方法论上，它成功示范了如何将基于人群的遗传因果推断（MR）与高分辨率的单细胞肿瘤图谱（scRNA-seq）及体内外功能验证相结合，构建了一个从基因到表型、从群体到单细胞、从因果关联到机制阐释的完整研究范式。这为未来在复杂疾病中鉴定和验证新的治疗靶点提供了可借鉴的蓝图。其次，在科学发现上，该研究明确了DKK3在ESCC中的致癌作用及其机制，解决了先前关于该基因在ESCC中功能角色的争议，为将DKK3开发为ESCC的潜在生物标志物、预后指标乃至药物治疗靶点奠定了坚实的理论基础。最后，对PI3K-AKT通路的牵连也为探索联合靶向治疗提供了新思路。总之，这项研究不仅加深了我们对食管鳞癌发病机制的理解，也为未来的精准医疗策略带来了新的希望。
