1. 领域背景与文献
文献英文标题:RAS pathway mutations orchestrate chromatin reconfiguration and transcriptional rewiring in colorectal cancer
发表期刊:文献未明确提供;影响因子:未公开;研究领域:结直肠癌表观遗传调控,肿瘤学。
领域共识:结直肠癌是全球高发的消化道恶性肿瘤,基因组不稳定和非突变表观重编程是肿瘤恶性进展的核心特征。下一代测序技术的发展推动了多组学整合分析,逐步揭示了遗传突变、表观景观与转录程序之间的互作关系。大约40%的结直肠癌携带RAS通路突变,其中KRAS突变占比约40%中的多数,10%~20%的病例携带BRAF突变,不同突变类型的肿瘤具有 distinct的突变、转录和表观组学特征。现有研究已经构建了多种细胞和动物模型解析RAS通路突变的功能,证实RAS突变可通过上调葡萄糖转运蛋白1增强癌细胞低糖存活能力,也发现TCF7L2可能参与RAS驱动的葡萄糖代谢调控。
当前领域仍存在明确研究空白:既往CRISPR筛选仅靶向约3000个代谢基因,覆盖范围有限,TCF7L2等潜在调控因子在RAS通路突变中的作用尚未充分探索;同时HCT116和RKO两种常用结直肠癌细胞系虽同属于微卫星不稳定、CIMP+、CMS4亚型,但二者对β-连环蛋白介导存活的依赖性差异显著,RAS通路突变与WNT信号通路的互作机制在不同遗传背景下尚未阐明。针对这一问题,本研究构建了不同RAS突变状态的同基因结直肠癌细胞模型,利用转座酶可及性染色质测序(ATAC-seq)解析突变驱动的染色质重构与转录重编程,明确关键调控因子的临床意义,为结直肠癌的预后分层与靶向干预提供新的分子基础。
2. 文献综述解析
作者对现有研究的梳理按照研究方向分为两个核心维度,分别是RAS通路突变在结直肠癌中的功能研究、RAS突变与WNT信号通路的交互作用研究。现有研究已经明确,体细胞RAS通路突变是结直肠癌的常见驱动事件,不同突变类型具有不同的分子特征;基因编辑同基因模型的研究证实,敲除HCT116细胞的突变KRAS等位基因可降低细胞锚定非依赖性生长、体外增殖和体内成瘤能力,下调MYC表达,而敲除RKO细胞的突变BRAF等位基因仅降低MEK和ERK磷酸化,对细胞形态和增殖影响较小。技术层面,同基因细胞模型排除了遗传背景混杂的干扰,结果可靠性高于混合样本队列研究,既往研究也发现RAS通路突变可上调葡萄糖转运蛋白1,帮助癌细胞在低糖条件下存活,正向遗传筛选提示TCF7L2参与这一调控过程。
现有研究的局限性主要体现在三个方面:一是既往功能筛选仅覆盖约3000个代谢基因,大量非代谢调控因子的作用未被探索,TCF7L2的具体功能仍不明确;二是既往研究已经观察到不同细胞系对β-连环蛋白的依赖性存在显著差异,但RAS突变与WNT信号的互作机制在不同背景下尚未阐明;三是细胞模型发现的调控规律缺乏大样本临床队列的验证,其临床转化价值未被证实。
本研究的创新价值十分明确,相较于现有研究,本研究首次利用两种不同背景的同基因结直肠癌模型,整合ATAC-seq和转录组测序系统解析RAS通路突变驱动的全基因组染色质重构与转录重编程,明确了RUNX2和TCF7L2分别在KRAS和BRAF突变中的核心调控作用,首次证实TCF7L2表达升高是BRAF突变结直肠癌的独立预后因素,填补了RAS通路突变表观调控机制领域的研究空白,同时验证了细胞模型发现的调控规律在多临床队列中的普遍性,为结直肠癌的精准预后提供了新的分子依据。
3. 研究思路总结与详细解析
本研究的整体研究目标为解析RAS通路突变在结直肠癌中驱动的表观遗传与转录调控程序,核心科学问题是RAS通路突变如何通过重塑染色质可及性调控转录程序,哪些关键转录因子介导这一过程,这些调控改变是否与临床预后相关;整体技术路线为:构建仅RAS突变状态不同的同基因结直肠癌细胞模型→ATAC-seq检测全基因组染色质可及性,整合转录组数据→差异分析识别突变相关通路和转录因子→临床多队列验证表达规律→生存分析明确预后价值,形成了从模型探索到临床验证的完整研究闭环。
3.1 同基因细胞模型构建
实验目的:构建仅RAS突变状态不同的同基因结直肠癌细胞模型,排除遗传背景混杂对实验结果的干扰。方法细节:采用重组腺相关病毒介导的同源重组技术在内源基因位点进行敲除或敲入操作,为避免引入 unintended单核苷酸多态性,每个载体的同源臂直接从亲本细胞基因组扩增,而非使用标准参考模板,通过重叠PCR引入目标点突变,采用Gateway克隆技术组装最终靶向载体。具体构建策略为:亲本HCT116携带1个突变KRAS和1个野生型KRAS等位基因,敲除突变KRAS得到保留野生型的DK克隆,再在DK克隆基础上敲入1个突变BRAF得到MHB克隆;亲本RKO携带2个突变BRAF和1个野生型BRAF等位基因,敲除两个突变BRAF得到保留野生型的RB克隆,再分别在RB克隆基础上敲入不同热点突变KRAS得到不同突变体克隆。每个克隆制备3~4个生物学重复用于后续实验。结果解读:成功构建了遗传背景一致仅RAS突变状态不同的同基因细胞克隆,为后续实验提供了合格的研究模型。实验所用关键产品:Thermo Fisher Scientific的RIPA缓冲液(货号B9901)、BCA蛋白定量试剂盒(货号23225)、Western blot剥离缓冲液(货号46430)、SuperSignal West PICO PLUS化学发光试剂(货号34577);Sigma-Aldrich的蛋白酶磷酸酶抑制剂(货号PPC1010);Cell Signaling Technology的磷酸化ERK1/2(货号4370)、总ERK1/2抗体(货号4695);Santa Cruz Biotechnology的β-肌动蛋白抗体(货号sc-47778);Qiagen的MinElute PCR纯化试剂盒(货号28004)。
3.2 染色质可及性测序与整合分析
实验目的:检测不同RAS突变状态下的全基因组染色质可及性,识别突变相关的差异调控区域和通路。方法细节:取50000个细胞采用Omni ATAC方案构建测序文库,使用NP-40、吐温-20和洋地黄毒苷透化减少线粒体DNA污染,Tn5转座酶37℃孵育30分钟进行片段化,纯化后制备文库进行测序;下机数据经过比对、过滤低质量reads、去除PCR重复、校正Tn5偏移、调用峰、去除黑名单区域、标准化后,采用主成分分析和 Pearson相关性分析评估生物学重复重复性,使用DESeq2进行差异分析,对峰进行基因组注释后开展功能富集分析,整合chromVAR和TOBIAS分析转录因子结合活性。结果解读:主成分分析显示,同基因克隆按照RAS突变状态在表观和转录水平明确分层(图1C),亲本突变克隆和衍生野生型克隆的差异最大,不同KRAS热点突变的RKO克隆聚集性良好,说明RAS突变是驱动染色质和转录组差异的核心因素。差异分析结果显示,敲除HCT116突变KRAS得到3656个上调基因、4092个下调基因(调整后P<0.05,|log2倍数变化|>0.58,n=3~4),同时得到3327个上调开放染色质区域、1961个下调区域(调整后P<0.05,|log2倍数变化|>0.58,n=3~4);敲除RKO突变BRAF得到2694个上调基因、2933个下调基因(调整后P<0.05,|log2倍数变化|>0.58,n=3~4),同时得到4176个上调开放染色质区域、4504个下调区域(调整后P<0.05,|log2倍数变化|>0.58,n=3~4)。功能富集分析证实,WNT信号通路在RAS突变存在时显著富集,说明RAS突变在转录和表观水平共同调控WNT信号通路。
3.3 关键转录因子识别与靶基因验证
实验目的:基于差异染色质可及性识别RAS突变相关的关键转录因子,验证其对WNT信号通路靶基因的直接调控作用。方法细节:整合chromVAR的染色质变异分析和TOBIAS的转录因子足迹分析,结合转录组表达数据,筛选结合活性变化与表达变化一致的转录因子,然后将候选转录因子靶基因与WNT信号通路基因取交集,预测转录因子在靶基因启动子区的结合位点。结果解读:分析发现,KRAS突变HCT116中,RUNX2的结合活性和表达变化与突变状态高度一致,BRAF突变RKO中,TCF7L2的结合活性和表达变化与突变状态高度一致。进一步靶基因分析显示,RUNX2调控8个WNT信号通路靶基因,其中VGLL4启动子区存在明确的RUNX2结合足迹;TCF7L2调控11个WNT信号通路靶基因,其中AXIN2启动子区存在明确的TCF7L2结合足迹,证实RUNX2和TCF7L2可直接结合WNT信号通路靶基因的调控区域,介导RAS突变对WNT信号的调控。
3.4 临床队列验证与预后分析
实验目的:验证细胞模型发现的突变相关基因和转录因子在临床肿瘤中的表达规律,分析其预后价值。方法细节:纳入U-CAN、TCGA COAD、AC-ICAM三个独立临床队列,比较RAS突变型和野生型患者的基因表达差异,筛选表达调控方向与细胞模型一致的基因;采用Cox比例风险模型进行单因素和多因素生存分析,以75百分位数对表达水平分层,绘制Kaplan-Meier生存曲线;利用TCGA临床样本ATAC-seq数据,分析转录因子结合活性与突变状态的关联。结果解读:两个细胞模型共得到85个共同突变相关基因,其中8个基因在三个临床队列中均显示一致的表达变化;生存分析显示,在U-CAN队列中,仅高PPFIBP1表达与不良总生存独立相关;对RUNX2和TCF7L2的临床验证显示,TCF7L2表达在BRAF突变结直肠癌患者中持续升高,结合活性也与BRAF突变状态显著相关,多因素生存分析显示,高TCF7L2表达与不良总生存独立相关(调整后P<0.05)。
4. Biomarker研究及发现成果
Biomarker定位:本研究鉴定的核心预后生物标志物是转录因子TCF7L2,属于表达型预后生物标志物,筛选验证逻辑完整清晰:同基因细胞模型通过ATAC-seq和转录组整合分析筛选差异转录因子→两种不同背景细胞模型交叉验证→多临床队列通过转录组和ATAC-seq验证表达与结合活性→单因素和多因素生存分析明确独立预后价值。
研究过程详述:TCF7L2的表达数据来源于细胞模型转录组和三个临床队列的批量RNA测序,结合活性数据来源于TCGA临床队列的批量ATAC测序;验证方法包括差异表达分析、染色质结合活性变异分析、Cox比例风险回归模型生存分析;研究结果显示,TCF7L2在BRAF突变结直肠癌患者中的表达水平显著高于RAS野生型患者(Tukey检验,P<0.05,三个独立队列结果一致),多因素Cox分析显示,高TCF7L2表达与不良总生存独立相关(调整后P<0.05),ROC曲线相关特异性和敏感性数据文献未明确提供。
核心成果提炼:TCF7L2是BRAF突变结直肠癌相关的独立预后生物标志物,高表达提示患者总生存期更差,其预后价值独立于年龄、分期等协变量(P<0.05,三个临床队列合计样本量超过一千例,具体样本量文献未明确提供)。本研究首次在同基因结直肠癌模型中明确,RAS通路突变通过染色质重构上调TCF7L2表达,进而激活WNT信号通路,修正了既往研究认为TCF7L2是结直肠癌抑癌基因的观点,证实其在BRAF突变结直肠癌中发挥致癌作用,为BRAF突变结直肠癌的预后分层提供了新的可靠分子标志物,也为后续靶向干预提供了新的候选靶点。