EWSR1过表达调节结直肠癌中增殖-迁移权衡：一项多队列转录组分析

《Journal of Bio-X Research》：EWSR1 Overexpression Modulates the Proliferation–Migration Trade-off in Colorectal Cancer: A Multicohort Transcriptomic Analysis

　　

背景：结直肠癌（colorectal cancer, CRC）的进展涉及增殖状态与迁移肿瘤状态之间的动态转换，这通常被描述为增殖-迁移权衡。然而，肿瘤细胞可能表现出反映表型可塑性的混合表型。RNA结合蛋白Ewing sarcoma breakpoint region 1（EWSR1）调节转录和RNA加工，但其在融合阴性CRC中协调这些表型的作用仍不清楚。方法：研究人员对来自3个队列（The Cancer Genome Atlas、Clinical Proteomic Tumor Analysis Consortium和Sidra–Leiden University Medical Center）的988例CRC肿瘤进行了整合转录组分析。使用单样本基因集富集分析（single-sample gene set enrichment analysis, ssGSEA）并经过ComBat批次校正后，量化了增殖、上皮-间质转化（epithelial–mesenchymal transition, EMT）和基质程序。通过将EMT评分对基质特征进行残差化处理，得到了肿瘤内在迁移信号。使用相关性分析、多变量回归和交互作用模型评估关联。结果：原始EMT评分与基质特征强相关（r = 0.964, P < 0.001）。调整后，增殖与迁移之间出现了清晰的负相关关系（r = ?0.226, P < 0.001）。EWSR1表达与增殖强相关（r = 0.341, P < 0.001），与迁移弱相关（r = 0.076, P = 0.017），但与基质含量无关。在多变量模型中，增殖负向预测迁移（β = ?0.286, P < 0.001），而EWSR1正向预测迁移（β = 0.174, P < 0.001）。交互作用不显著（P = 0.630），但分层分析显示，在EWSR1高表达肿瘤中，该权衡有所减弱。结论：EWSR1表达与迁移相关活性相关，而未耦合增殖和迁移；相反，EWSR1表达、迁移和增殖之间的关系表明，EWSR1通过部分放松CRC中增殖-迁移权衡而参与表型可塑性。

**研究背景**
结直肠癌（CRC）是全球癌症发病和死亡的主要原因之一，其进展涉及遗传改变和转录重编程的复杂相互作用。经典驱动突变（如APC、KRAS、TP53、BRAF）已被广泛研究，近年关注转向调控动态肿瘤表型（包括增殖、侵袭和上皮可塑性）的转录和转录后调控因子。肿瘤细胞并非静态，而是表现出显著可塑性，根据内在和微环境信号在增殖和迁移状态之间转换。肿瘤生物学中的核心概念是“go-or-grow”假说，即细胞增殖与迁移之间存在互斥关系。然而，新证据表明这种二分法并非绝对，在某些调控条件下，肿瘤细胞可能呈现同时支持生长和侵袭潜能的混合表型。识别协调或耦合这些过程的分子调控因子仍是癌症生物学的关键挑战。RNA结合蛋白EWSR1（Ewing sarcoma breakpoint region 1）在尤因肉瘤中因参与致癌融合蛋白而闻名，但其在融合阴性CRC中的过表达作用尚不明确。EWSR1作为转录和RNA加工的多功能调节因子，影响替代剪接、mRNA运输以及与细胞生长、分化和细胞骨架组织相关的基因表达程序。已有证据提示EWSR1可能通过调节上皮-间质转化（EMT）、基质相互作用和应激反应途径促进癌症进展。然而，EWSR1是否独立关联这些表型或调节它们之间的关系，从而促成协调或耦合的肿瘤状态，仍不清楚。因此，研究人员开展本研究，假设EWSR1表达调节CRC中增殖与迁移相关转录程序之间的关系，旨在区分两种模型：（a）表型协调——EWSR1独立影响增殖和迁移；（b）功能耦合——增殖与迁移之间的关系本身依赖于EWSR1表达。该研究发表在《Journal of Bio-X Research》。

**研究方法**
研究人员采用多队列转录组驱动分析框架，整合来自3个独立CRC队列（The Cancer Genome Atlas [TCGA]、Clinical Proteomic Tumor Analysis Consortium [CPTAC-2]和Sidra–Leiden University Medical Center [Sidra-LUMC]）的988例肿瘤RNA测序数据。数据经基因符号统一化后合并，并使用ComBat算法进行批次校正。使用单样本基因集富集分析（ssGSEA）量化转录表型：增殖评分基于Hallmark E2F靶点和G2/M检查点基因集；迁移相关活性通过Hallmark EMT基因集近似量化；基质含量使用细胞外基质富集基因特征量化。为分离肿瘤内在迁移信号，通过线性回归对EMT评分进行基质混杂调整，残差作为基质校正后的EMT评分（EMT_adjusted）。所有连续变量在后续分析前进行z-score标准化。主要技术方法：1）ssGSEA（单样本基因集富集分析）计算每个样本的富集评分；2）ComBat批次校正消除技术变异性；3）线性回归残差化处理去除基质混杂；4）Pearson相关性分析评估变量间基线关系；5）多变量线性回归模型评估增殖和EWSR1对迁移的独立贡献；6）交互作用建模（增殖×EWSR1）检验EWSR1是否修饰增殖-迁移关系。所有分析在SPSS 29和R 4.5.2中进行，双尾P<0.05视为显著。

**研究结果**
**1. 基质混杂主导CRC中原始EMT信号**：首先评估增殖、迁移（EMT）和基质程序之间的成对关系。增殖与EMT评分之间呈弱但显著的负相关（r=?0.098，P=0.002）。但EMT评分与基质特征强相关（r=0.964，P<0.001，R²≈0.93），表明基质含量解释了EMT评分的大部分方差。经基质调整后，EMT信号与基质评分无关联（r=0.003，P=0.937），确认了混杂的有效去除。使用基质校正后EMT评分时，增殖与迁移之间负相关更加显著（r=?0.226，P<0.001），支持CRC中存在内源性“go-or-grow”约束。

**2. EWSR1表达与增殖肿瘤状态强相关，与迁移弱相关**：相关性分析显示EWSR1表达与增殖强正相关（r=0.341，P<0.001），与迁移活性在基质调整前后均呈弱正相关（调整前r=0.063，P=0.046；调整后r=0.076，P=0.017）。EWSR1表达与基质评分无关联（r=0.045，P=0.159），表明其反映肿瘤内在生物学而非微环境组成。

**3. 增殖和EWSR1独立预测迁移相关活性**：以基质校正后EMT评分（EMT_adjusted）为因变量的多变量线性回归显示，增殖与迁移强负相关（β=?0.286，P<0.001），EWSR1表达与迁移正相关（β=0.174，P<0.001），基质评分不显著（P=0.591）。反向模型（以增殖为因变量）中，EMT_adjusted（β=?0.253，P<0.001）和EWSR1表达（β=0.360，P<0.001）均为显著独立预测因子。结果确认了增殖与迁移的负相关，且EWSR1与两个过程独立关联，但其对迁移的预测方向为正，表明EWSR1可能部分削弱增殖-迁移权衡而非直接耦合。

**4. 交互作用分析显示EWSR1的不对称情境依赖性调节**：以EMT_adjusted为因变量的模型中，增殖×EWSR1交互项不显著（β=?0.016，P=0.630），但增殖和EWSR1的主效应显著。以增殖为因变量的反向模型中，EWSR1×EMT_adjusted交互项显著（β=0.103，P<0.001），表明高EWSR1表达减弱了迁移活性与增殖之间的负关联。这说明EWSR1通过影响迁移活性对增殖的抑制作用而发挥方向依赖性（不对称）调节。

**5. EWSR1表达与增殖-迁移关系减弱相关**：基于中位分层的分析显示，EWSR1低表达肿瘤中增殖与EMT_adjusted呈强负相关（r=?0.310，P<0.001），而高表达肿瘤中相关性减弱（r=?0.153，P<0.001），Fisher’s r-to-z转换证实组间差异显著（Z=?2.61，P=0.009）。采用低三分位数阈值时结果一致（r=?0.340 vs ?0.184，Z=?2.55，P=0.011）。协方差分析（ANCOVA）也显示增殖×EWSR1分组的交互作用显著（F=4.439，P=0.035），支持EWSR1以阈值依赖方式削弱增殖-迁移二分法。

**6. EWSR1表达与基因组不稳定性和分子亚型相关**：非参数分析显示EWSR1表达与BRAF突变（P=1.99×10^?6）、微卫星不稳定性（MSI）状态（P=1.10×10^?6）、肿瘤突变负荷（P=5.79×10^?6）、非整倍体（P=1.08×10^?5）和基因组改变分数（P=2.44×10^?6）显著相关。分子亚型间差异显著（P=3.51×10^?7），高突变/MSI亚型与间质型和上皮型亚型均不同。EWSR1表达与总生存期（OS）或无病生存期（DFS）无显著关联（HR=0.945，P=0.317），表明其反映表型转录状态而非直接预后。

**7. EWSR1独立于MSI状态削弱增殖-迁移二分法**：在总队列和MSI分层子集中，多变量回归显示EWSR1×EMT_adjusted交互项在调整MSI后仍显著（β=0.117 vs 0.125，均P<0.001），MSI本身不改变效应大小。在MSI和微卫星稳定（MSS）子集中，EWSR1表达与增殖的正相关强度相似（MSI: β=0.351; MSS: β=0.358），表明EWSR1的调节效应独立于MSI。

**总结讨论**：研究表明EWSR1表达调节CRC中增殖与迁移转录程序之间的关系，支持表型可塑性模型而非功能耦合。原始EMT信号受基质严重混杂，调整后出现清晰的增殖-迁移拮抗。EWSR1与增殖强相关，并独立促进迁移活性，但不耦合两者。交互作用分析显示EWSR1不对称地削弱增殖-迁移约束，且该效应独立于MSI状态。EWSR1在基因组不稳定和高突变亚型中富集，提示其与转录可塑性相关。研究局限性包括：使用bulk转录组数据需计算校正基质，ssGSEA提供相对富集而非直接功能测量，缺乏实验验证。未来需单细胞转录组和功能扰动模型验证。

**研究结论**：EWSR1表达与CRC中增殖肿瘤状态相关，并独立关联迁移相关活性。但它并未功能性地耦合增殖和迁移；相反，它通过部分放松增殖-迁移拮抗而促进可塑性肿瘤表型。这些发现将EWSR1突出为肿瘤适应性的潜在生物标志物，并强调表型可塑性是CRC中临床相关的靶点。