《Translational Oncology》：Integrated bulk and single-cell transcriptomic analysis of ferroptosis-related programs in high-grade serous ovarian carcinoma
编辑推荐：
中文摘要翻译
铁死亡（ferroptosis）是一种由铁依赖性脂质过氧化物积累驱动的调节性细胞死亡形式，其在高级别浆液性卵巢癌（high-grade serous ovarian carcinoma，HGSOC）中的转录特征及临床意义尚不明确。研究人员整合了T
中文摘要翻译
铁死亡（ferroptosis）是一种由铁依赖性脂质过氧化物积累驱动的调节性细胞死亡形式，其在高级别浆液性卵巢癌（high-grade serous ovarian carcinoma，HGSOC）中的转录特征及临床意义尚不明确。研究人员整合了TCGA-OV队列（n=421例原发肿瘤）的bulk RNA-seq数据与治疗初治的GSE300897单细胞队列（9例供体，16,939个细胞）数据，对25个铁死亡相关基因（ferroptosis-related genes，FRGs）进行了跨平台分析。研究结果显示：在TCGA-OV队列中，聚合FRG评分与总生存期（overall survival，OS）无显著关联（中位数分层log-rank p=0.80；多变量HR=1.46，95% CI 0.81–2.63，p=0.21），但该评分与细胞毒性T细胞（ρ=0.38）、CD8+ T细胞（ρ=0.37）及巨噬细胞/单核细胞（ρ=0.30）浸润评分显著正相关，且FRG高评分肿瘤中HAVCR2、TIGIT、PDCD1和CTLA4等免疫检查点基因表达升高。单细胞分析显示，髓系细胞具有最高的FRG评分（均值0.83），并优先表达铁代谢相关基因FTH1、FTL、TFRC和NCOA4。化疗难治性肿瘤（n=4）与化疗敏感性肿瘤（n=5）相比，髓系细胞比例升高2.5倍（27.2%对10.7%），成纤维细胞和浆细胞比例降低；在恶性细胞中，GLS（log2FC=+0.88）、CS（+0.68）和TFRC（+0.58）趋向升高，而抗氧化臂基因SLC7A11（?1.32）、HMOX1（?1.08）和GPX4（?0.82）趋向降低，但均未达到BH校正显著性。跨平台相关性分析显示，FRG评分与缺氧（TCGA-OV ρ=0.39；GSE300897 ρ=0.44）、糖酵解（0.44对0.37）和NF-κB（0.17对0.28）程序在两种平台上均呈正相关，而氧化磷酸化（OXPHOS）仅在bulk水平呈正相关（ρ=0.41），单细胞水平无相关性（ρ=?0.02）。qRT-PCR验证实验确认SKOV3和OVCAR3细胞中FTH1、FTL、GPX4和SLC7A11的相对表达量均高于IOSE80正常卵巢上皮细胞。该研究提供了细胞类型分辨的、基于关联的铁死亡相关转录程序特征，为后续功能验证奠定了假设基础。
论文解读
一、研究背景与问题提出
卵巢癌是妇科恶性肿瘤中致死率最高的类型，在发达国家女性癌症相关死亡原因中位居第五。高级别浆液性卵巢癌（high-grade serous ovarian carcinoma，HGSOC）是其最主要的组织学亚型，多数患者确诊时已处于FIGO III–IV期。尽管HGSOC对铂类-紫杉醇联合化疗初始敏感，但约70–80%的患者最终复发并产生耐药，后续治疗基本转为姑息性。因此，阐明卵巢癌进展和治疗耐药的分子机制是当前研究的迫切任务。
铁死亡（ferroptosis）是2012年首次系统描述的一种非凋亡性调节性细胞死亡形式，其特征为铁依赖性脂质氢过氧化物在细胞膜上的致死性积累。GPX4作为铁死亡的主要抑制因子，通过将磷脂氢过氧化物还原为无毒羟基产物发挥保护作用；SLC7A11（xCT）通过导入胞外胱氨酸维持谷胱甘肽合成以支持GPX4活性；而ACSL4通过将多不饱和脂肪酸掺入磷脂促进铁死亡易感性，NCOA4介导的铁蛋白自噬则从铁蛋白储存（FTL/FTH1）中释放不稳定铁以放大氧化应激。
尽管已有研究报道GPX4和SLC7A11在卵巢癌细胞中高表达并与铂类反应相关，但在HGSOC中尚缺乏聚焦的跨平台分析来定义铁死亡相关转录程序及其与免疫和代谢状态的关系。具体而言，三个关键问题尚未解决：第一，紧凑的铁死亡相关基因（ferroptosis-related genes，FRGs）程序在bulk肿瘤水平是否携带临床预后信息，抑或主要反映细胞组成差异；第二，HGSOC中哪些细胞类型对铁死亡相关转录信号贡献最大；第三，铁死亡相关转录在bulk和单细胞平台上与免疫、炎症和代谢程序的共变程度如何。此外，肿瘤免疫微环境（tumor immune microenvironment，TIME）已被证实是卵巢癌预后和治疗反应的关键决定因素，但不同细胞区室中铁死亡基因特征是否参与免疫调节仍不清楚。
二、研究内容与主要结论
研究人员整合了TCGA-OV bulk RNA-seq数据（421例原发肿瘤）与治疗初治的GSE300897单细胞队列（9例HGSOC患者，含4例化疗难治性和5例化疗敏感性肿瘤），围绕三个核心问题展开分析：25基因FRG panel在HGSOC肿瘤中的表达与共调控模式；聚合FRG评分与生存及免疫/基质组成的关系；哪些细胞区室和转录程序解释铁死亡相关异质性，包括化疗难治性与敏感性供体间的差异。
研究的主要发现包括：聚合FRG评分在TCGA-OV中不具预后价值，但该评分与免疫浸润和免疫检查点表达显著相关；单细胞分析揭示髓系细胞具有最高的FRG评分并优先表达铁代谢基因；化疗难治性肿瘤以组成重塑为主，恶性细胞中抗氧化臂基因趋向降低；FRG评分与缺氧、糖酵解和NF-κB程序在两种平台上均呈正相关，而OXPHOS耦合具有平台依赖性。该研究的主要贡献在于提供了跨平台、细胞类型分辨的铁死亡相关转录特征，并生成可供功能验证的假设。
三、主要关键技术方法
研究人员采用了以下关键技术方法：TCGA-OV bulk转录组数据来自GDC门户，单细胞数据来自GEO accession GSE300897（DECIDER观察性试验NCT04846933，9例治疗初治HGSOC活检）。bulk数据以log2(CPM+1)标准化，25基因FRG panel按功能模块分组（铁代谢、抗氧化系统、脂代谢、调控、代谢/其他），采用Scanpy score_genes实现计算per-sample FRG评分。生存分析采用Kaplan–Meier估计、log-rank检验和时间依赖性ROC曲线；免疫微环境通过13个谱系marker基因特征评分量化。单细胞数据经质控后保留16,939个细胞和32,300个基因，采用Leiden聚类（resolution 0.5）和UMAP可视化，细胞身份采用原始注释（9种亚型）。通路评分覆盖I型干扰素（IFN-I）、缺氧、NF-κB、JAK–STAT、糖酵解、氧化磷酸化（OXPHOS）和脂肪酸氧化（FAO）。化疗难治性与敏感性比较采用供体水平pseudobulk分析（Welch t检验，BH校正）。qRT-PCR验证在IOSE80、SKOV3和OVCAR3细胞系中进行，采用2^-ΔΔCt方法。
四、研究结果
铁死亡基因在TCGA-OV中广泛表达且弱协调。 25个FRGs在421例TCGA-OV肿瘤中的中位表达跨度超过一个数量级，铁储存基因FTL（中位11.2 log2(CPM+1)）和FTH1（10.4）及核糖体基因RPL8（11.2）表达最高，而促铁死亡基因NOX1（0.6）、SLC7A11（3.6）和LPCAT3（4.3）处于低表达端。相关性分析显示最强的正相关对为LPCAT3–BECN1（r=0.49）、FTL–HMOX1（r=0.48）和FTL–GPX4（r=0.47）。主成分分析中PC1和PC2分别捕获panel方差的14.8%和12.0%。
bulk FRG评分不具预后价值，但个别基因显示名义关联。 中位数分层（log-rank p=0.80）和三分位分层（log-rank p=0.24）均未显示生存差异；时间依赖性ROC的1、3、5年AUC分别为0.52、0.54和0.46。单基因Cox分析发现SLC7A11（HR=0.86，p=0.004）和PCBP2（HR=0.77，p=0.013）与较长生存相关，而LPCAT3（HR=1.24，p=0.011）与较短生存相关。多变量模型中FRG评分无独立关联（HR=1.46，p=0.21）。
FRG评分追踪TCGA-OV中炎症性、细胞毒性免疫微环境。 FRG评分与细胞毒性（ρ=0.38）、CD8+ T细胞（ρ=0.37）、巨噬细胞/单核细胞（ρ=0.30）、Treg（ρ=0.27）和树突状细胞（ρ=0.24）评分显著正相关。FRG高评分肿瘤中HAVCR2（Δ=+0.49，p<0.001）、TIGIT（Δ=+0.22，p=0.002）、PDCD1（Δ=+0.21，p=0.005）和CTLA4（Δ=+0.16，p=0.003）表达升高。CD8A（r=0.23）和GZMB（r=0.29）随FRG评分单调增加。
化疗难治性与敏感性HGSC的单细胞景观。 质控后保留16,939个细胞（88.4%），UMAP分离出九种细胞亚型。难治性肿瘤中髓系细胞比例显著升高（27.2%对10.7%），成纤维细胞（9.5%对23.3%）和浆细胞（1.0%对5.3%）比例降低，T细胞比例相当（34.2%对34.8%）。
铁死亡程序在髓系区室中最强。 髓系细胞具有最高FRG评分（均值0.83），其次为肥大细胞（0.66）、树突状细胞（0.59）、成纤维细胞（0.56）和卵巢癌细胞（0.48），淋巴细胞最低。铁代谢模块基因FTH1、FTL、TFRC和NCOA4优先在髓系细胞表达。FRG评分与测序深度仅弱相关（ρ=0.12），排除文库大小伪影。
化疗难治性重塑组成并改变癌细胞中抗氧化铁死亡臂。 组成差异超过转录差异。恶性细胞中GLS（log2FC=+0.88）、CS（+0.68）和TFRC（+0.58）趋向升高，而抗氧化臂SLC7A11（?1.32）、HMOX1（?1.08）和GPX4（?0.82）趋向降低，但均未达到BH校正显著性。
难治性肿瘤结合较低的IFN-I活性和髓系炎症偏移。 IFN-I活性在难治性肿瘤中平均降低0.15分，内皮细胞中降幅最大（Δ=?0.32）；缺氧（+0.04）和NF-κB（+0.03）略升高。FRG与IFN-I评分在单细胞水平几乎不相关（ρ=0.02）。FTH1与IL6（ρ=0.15）、FTL与IL1B的相关性为正向但较弱。
程序共激活跨平台复制，而单细胞OXPHOS耦合未复制。 FRG评分与缺氧（TCGA-OV ρ=0.39；GSE300897 ρ=0.44）、糖酵解（0.44对0.37）和NF-κB（0.17对0.28）在两种平台均呈正相关；OXPHOS在bulk中正相关（ρ=0.41）但单细胞中无相关性（ρ=?0.02）。FTH1–LDHA共表达弱但显著（ρ=0.17，p=8.4×10?21）。SC-derived hub模块（FTH1/LDHA/GPX4/NCOA4）在TCGA-OV中不区分生存（log-rank p=0.62）。
qRT-PCR验证。 FTH1、FTL、GPX4和SLC7A11在SKOV3和OVCAR3中的相对表达量均高于IOSE80。FTH1在IOSE80、SKOV3和OVCAR3中分别为1.00±0.06、1.84±0.12和2.21±0.17；FTL为1.00±0.05、1.62±0.11和1.95±0.13；GPX4为1.00±0.05、1.48±0.10和1.71±0.12；SLC7A11为1.00±0.07、2.05±0.16和2.46±0.19。
五、讨论与结论总结
讨论部分指出，该研究不支持将聚合FRG评分作为预后生物标志物，这一阴性结果具有信息价值，表明无方向性的促/抗铁死亡基因复合物不应被解读为肿瘤侵袭性或铁死亡敏感性的简单度量。在bulk水平，FRG评分追踪免疫浸润表型，但这一关联可能源于铁处理、抗氧化防御、代谢状态或FRG富集细胞群体丰度的变化，而非铁死亡本身驱动免疫浸润。单细胞结果支持组成性解释：髓系细胞是铁死亡相关转录信号的主要来源，但高FRG评分不应等同于铁死亡执行增加。化疗难治性肿瘤的组成重塑超过基因表达差异，恶性细胞中抗氧化臂的趋向性降低值得未来功能实验验证。跨平台可复制的缺氧/糖酵解/NF-κB正相关识别出可优先进行功能测试的共激活转录程序。单细胞hub网络将FTH1、LDHA和RELA定位在铁处理、代谢和炎症信号的交界面，但这些观察仅为假设生成性质。
研究结论为：通过整合TCGA-OV bulk转录组与供体感知的单细胞分析，该研究识别出一个不具独立预后价值但与免疫组成及特定炎症/代谢状态密切相关的铁死亡相关转录程序。髓系细胞显示最高聚合FRG评分并优先表达铁处理基因，化疗难治性肿瘤表现出显著的组成重塑和恶性细胞中若干抗氧化FRGs的探索性降低。缺氧、糖酵解和NF-κB程序与FRG评分在跨平台中显示可复现的正相关，而OXPHOS耦合具有平台依赖性。这些发现定义了细胞类型分辨的、基于关联的假设，供未来功能验证，而非已确立的铁死亡机制或治疗靶点。
