LBPO.BCS02 · 生物信息与计算 · Late-Breaking
ISWI双敲除特征可对PAX3/7-FOXO1融合阳性与融合阴性横纹肌肉瘤进行分层
ISWI double-knockout signatures stratify PAX3/7-FOXO1 fusion-positive vs fusion-negative rhabdomyosarcoma
作者与单位 Authors & Affiliations
摘要 Abstract
中文摘要
背景:肉瘤融合癌基因驱动独特的转录程序,这些程序与由ATP依赖性重塑复合物塑造的染色质景观相一致。在横纹肌肉瘤(RMS)中,融合阳性(FP)肿瘤携带PAX3/7-FOXO1融合,而融合阴性(FN)肿瘤缺乏FOXO1融合。我们旨在整合融合转录本检测与基因程序评分,以量化融合状态与染色质重塑模块活性的关系,并检验实验定义的ISWI扰动特征是否在融合阳性疾病中富集。
方法:我们分析了来自公共RMS队列GSE108022的批量转录组(肿瘤 n=101;正常骨骼肌 n=5;FN n=66;FP n=35)。融合状态和融合类别(FP对FN)通过融合转录本检测和数据集注释从RNA-seq中确定。我们对一个紧凑的染色质重塑面板计算ssGSEA/GSVA富集评分,该面板捕捉核心复合物程序(如BAF、ncBAF、NuRD、ISWI)以及ISWI/SMARCA相关的靶标特征,包括两个ISWI双敲除(DKO)响应程序(Smarca1/Smarca5;ISWI ATP酶亚基SMARCA1/SNF2L和SMARCA5/SNF2H),其来源于小鼠DKO与WT的差异表达并映射至人类同源基因(ISWI_DKO_targets_UP/DOWN)。使用双侧Wilcoxon秩和检验并经Benjamini-Hochberg校正评估组间差异(肿瘤对正常;肿瘤内FP对FN)。我们还使用仅含模块的逻辑回归分类器,采用重复5折交叉验证和AUROC评估融合类别判别能力。
结果:与正常肌肉相比,RMS肿瘤显示出多个染色质重塑程序显著更高的富集,包括BAF_core、NuRD_core、ncBAF和ISWI_core(所有FDR 6.500357e-04)。ISWI_DKO_targets_UP在肿瘤中低于正常肌肉(肿瘤中位z=-0.008对正常中位z=1.569;FDR 1.936228e-03)。ISWI_DKO_targets_DOWN在肿瘤中高于正常肌肉,但校正后仅显示临界趋势(肿瘤中位z=0.108对正常中位z=-0.592;FDR 5.148731e-02)。在肿瘤内,融合类别强烈分层了机制性ISWI依赖程序:ISWI_DKO_targets_DOWN在FP中显著高于FN(p=2.521292e-07;FDR=2.521292e-06;FP中位z=0.775;FN中位z=-0.422;差值=1.197)。其他FP-FN差异包括SMARCA1_targets_UP(FDR 2.206301e-02)。一个仅含ISWI+SMARCA1/5模块的逻辑回归模型以AUROC=0.876±0.012区分FP与FN肿瘤。
结论:一个整合染色质重塑复合物核心与同源映射ISWI DKO程序的最简ssGSEA面板,可识别RMS中融合相关的染色质重塑状态。ISWI_DKO_targets_DOWN在程序水平显示出最强的FP-FN分离,支持紧凑、可重现的融合类别分层,推动对融合驱动RMS中ISWI相关依赖性进行更广泛的验证和功能随访。这一以染色质为中心的映射策略为开发候选程序水平生物标志物、以及优先排序将融合驱动的转录输出与涉及细胞状态可塑性的染色质重塑程序相联系的可检验组合假设,提供了一条实用途径。
查看英文原文 English abstract
Background: Sarcoma fusion oncogenes drive distinct transcriptional programs that are aligned with chromatin landscapes shaped by ATP-dependent remodeling complexes. In rhabdomyosarcoma (RMS), fusion-positive (FP) tumors harbor PAX3/7-FOXO1 fusions, whereas fusion-negative (FN) tumors lack FOXO1 fusions. We aimed to integrate fusion transcript detection with gene program scoring to quantify how fusion status relates to chromatin remodeler module activity, and to test whether an experimentally defined ISWI perturbation signature is enriched in fusion-positive disease.
Methods: We analyzed bulk transcriptomes from the public RMS cohort GSE108022 (Tumor n=101; Normal skeletal muscle n=5; FN n=66; FP n=35). Fusion status and fusion class (FP vs FN) were determined from RNA-seq using fusion transcript detection and dataset annotations. We computed ssGSEA/GSVA enrichment scores for a compact chromatin remodeling panel capturing core complex programs (e.g., BAF, ncBAF, NuRD, ISWI) and ISWI/SMARCA-linked target signatures, including two ISWI double-knockout (DKO) response programs (Smarca1/Smarca5; ISWI ATPase subunits SMARCA1/SNF2L and SMARCA5/SNF2H) derived from mouse DKO vs WT differential expression and mapped to human orthologs (ISWI_DKO_targets_UP/DOWN). Group differences were assessed using two-sided Wilcoxon rank-sum tests with Benjamini-Hochberg correction (Tumor vs Normal; FP vs FN within tumors). We additionally evaluated fusion-class discrimination using a module-only logistic regression classifier with repeated 5-fold cross-validation and AUROC.
Results: Compared with normal muscle, RMS tumors showed significantly higher enrichment of multiple chromatin remodeling programs, including BAF_core, NuRD_core, ncBAF, and ISWI_core (all FDR 6.500357e-04). ISWI_DKO_targets_UP was lower in tumors than normal muscle (tumor median z=-0.008 vs normal median z=1.569; FDR 1.936228e-03). ISWI_DKO_targets_DOWN was higher in tumors than normal muscle but showed only a borderline trend after correction (tumor median z=0.108 vs normal median z=-0.592; FDR 5.148731e-02). Within tumors, fusion class strongly stratified the mechanistic ISWI dependence program: ISWI_DKO_targets_DOWN was markedly higher in FP than FN (p=2.521292e-07; FDR=2.521292e-06; FP median z=0.775; FN median z=-0.422; delta=1.197). Additional FP-FN differences included SMARCA1_targets_UP (FDR 2.206301e-02). An ISWI+SMARCA1/5 module-only logistic regression model discriminated FP from FN tumors with AUROC=0.876±0.012.
Conclusions: A minimal ssGSEA panel integrating chromatin remodeling complex cores with an ortholog-mapped ISWI DKO program identifies a fusion-linked chromatin remodeling state in RMS. ISWI_DKO_targets_DOWN shows the strongest FP-FN separation at the program level and supports compact, reproducible fusion-class stratification, motivating broader validation and functional follow-up of ISWI-associated dependencies in fusion-driven RMS. This chromatin-centric mapping strategy provides a practical route to develop candidate program-level biomarkers and to prioritize testable combination hypotheses linking fusion-driven transcriptional output with chromatin remodeling programs implicated in cell-state plasticity.
利益披露 Disclosure
A. K. Aljabri, None.