-
摘要:目的
探讨白蛋白与纤维蛋白原比值(albumin to fibrinogen ratio, AFR)与乳腺浸润性导管癌患者总生存期的关系。
方法收集2009年1月—2012年4月河北医科大学第四医院乳腺中心收治的230例乳腺浸润性导管癌患者血液及临床资料,回顾性分析患者术前临床资料,计算AFR,Kaplan-Meier法及Cox比例回归风险模型进行生存分析,采用二分类方式,分为AFR < 13.1与AFR≥13.1两组,采用单因素及多因素分析其与患者的生存关系,并建立预后模型。
结果AFR≥13.1时患者生存率降低,AFR是改善总生存期(OS)的独立预后因素,将多因素分析里P < 0.05的变量,即将AFR、白细胞(WBC)、雌激素受体(ER)孕激素受体(PR)、TNM分期、核分级、p53及脉管瘤栓7个指标建立OS预测模型,其受试者工作特征曲线(ROC)面积为0.882。
结论AFR低表达是改善乳腺癌OS的独立预后因素,AFR、WBC、ER或PR、TNM、核分级、p53及脉管瘤栓7个指标建立OS预测模型,具有很好预测乳腺浸润性导管癌患者生存率的作用。
-
关键词:
- 白蛋白与纤维蛋白原比值 /
- 乳腺浸润性导管癌 /
- 总生存期
Abstract:ObjectiveTo investigate the prognostic significance of albumin to fibrinogen ratio (AFR) in breast infiltrating ductal carcinoma.
MethodsWe collected 230 blood and clinical data of breast invasive ductal carcinoma admitted to the Breast Center of the Fourth Hospital of Hebei Medical University from January 2009 to April 2012. A retrospective study was conducted to collect fibrinogen and albumin indicators from preoperative patients. Survival analysis was performed using Kaplan-Meier method and Cox proportional regression risk model to explore the relation between AFR expression and prognosis. AFR was divided into two groups by two classifications: AFR < 13.1 and AFR≥13.1. The univariate and multivariate analyses were used to analysis the relation between AFR and patients' survival, then a prognosis model was established.
ResultsThe survival rate of patients was reduced when AFR≥13.1. In the multivariate analysis, the variables with P < 0.05 were used to create a model. AFR, white blood cell (WBC), estrogen receptor(ER) or progesterone receptor (PR), TNM, nuclear grading, p53, vascular tumor thrombus and other OS indicators were used to predict the survival time of patients. The receiver operating characteristic(ROC) curve area was 0.882.
ConclusionLow AFR expression is an independent prognostic factor for breast cancer. Seven indicators, including AFR, WBC, ER/PR, TNM, nuclear grading, p53 and vascular tumor thrombus, were used to create an OS predictive model, which could predict the prognosis of breast infiltrating ductal carcinoma patients.
-
0 引言
铁死亡是一种铁离子依赖的脂质过氧化物累积而导致的非凋亡形式的新型细胞死亡方式,2012年由Brent R. Stockwell提出[1],研究人员发现,小分子化合物Erastin可以特异性诱导Ras突变细胞发生铁离子依赖的、非凋亡性的细胞死亡,将这一过程称为铁死亡[2]。铁死亡的形态特征包括线粒体体积缩小、膜密度增高、线粒嵴减少或消失[3];其生化特点表现为细胞内铁和活性氧(reactive oxygen species, ROS) 的大量积累,产生过度的氧化应激以及脂质过氧化[4]。研究发现铁死亡与肿瘤放射治疗增敏、免疫治疗及耐药性方面密切相关[5-6]。TCGA数据分析发现高表达铁死亡相关转运体-SLC7A11的神经胶质瘤的患者,在接受放射治疗后肿瘤复发时间和总体生存期更短,提示铁死亡与神经胶质瘤患者放疗预后不良有关[7]。对传统疗法有抵抗力或具有高转移倾向的癌细胞可能特别容易发生铁死亡敏感,铁死亡诱导剂还可以与更传统的药物(如顺铂)协同作用,抑制肿瘤生长[8]。
前期研究中,我们构建与分析了肺腺癌(lung adenocarcinoma, LUAD)中与铁死亡有关的四个基因,其中肝型谷氨酰胺酶(GLS2)是LUAD预后的保护性因素,促进铁死亡。GLS2是谷氨酰胺酶的成员,它的活性在一些癌细胞中明显失调,GLS2是抑癌基因或癌基因取决于肿瘤的类型[9-10]。在乳腺癌中,GLS2与增强的体外细胞迁移和侵袭、间充质标志物和体内肺转移有关[11]。在胶质母细胞瘤中,GLS2可以抑制胶质母细胞瘤的生长[12]。GLS2介导的治疗在大规模临床环境中的应用仍处于初级阶段,其作为人类癌症的治疗及疾病进展标志仍存在争议。
本研究运用多个在线生物信息学平台和工具,评估了铁死亡相关基因GLS2在不同癌组织中的表达及其与癌症的可能联系,全面研究了GLS2在泛癌中的表达情况、预后和功能,分析GLS2是否可能成为癌症治疗及疾病进展的潜在靶点。
1 材料与方法
1.1 数据收集和处理
收集癌症基因组图谱(TCGA)数据库和广泛研究所癌症细胞系百科全书(CCLE)数据库的泛癌测序数据(Illumina平台),以及国际癌症基因组联盟(ICGC)数据库中的数据,所有数据均通过其门户网站进行分析。使用R包中的RMA功能(R studio版本:1.2.1335,R版本:3.6.1),过滤整个数据集,删除缺失和重复的结果,并通过log2(TPM+1)转换。通过门户网站提取相应的临床信息,包括患者的年龄、性别、肿瘤分期和临床分期。此外,仅从TCGA数据库下载肿瘤突变负荷(tumor mutational burden, TMB)和微卫星不稳定性(microsatellite instability, MSI)等信息。TMB计算为每百万碱基对的总突变发生率,MSI计算为基因重复序列中发生的插入或缺失事件数。
1.2 Cox回归分析和生存分析
Cox回归分析用于从R环境中的ICGC和TCGA数据库中检测GLS2的表达与每种癌症类型患者的总体生存率(overall survival, OS)、疾病特异性生存率(disease–specific survival, DSS)和无进展间隔期(progression–free interval, PFI)之间的相关性。根据最佳分离方法将患者分为高和低GLS2表达组后,采用Kaplan–Meier方法构建每种癌症类型患者的生存曲线。通过在R软件包(rdocumentation.org/packages/survival)中的survival ROC和survival分析,探究了生存的特异性和时间依赖性敏感性。使用对数秩检验检验曲线之间的差异,P<0.05为差异具有统计学意义。
1.3 免疫细胞浸润富集
肿瘤免疫评估资源数据库(TIMER)是一个用于免疫细胞浸润计算的数据库,可用于分析包括B细胞、CD4+T细胞、CD8+T细胞、巨噬细胞、中性粒细胞和树突状细胞的六种常见类型免疫细胞的浸润分数。TCGA数据库中泛癌数据的免疫细胞浸润分数用于分析其与GSL2表达的相关性。
1.4 统计学方法
Spearman相关检验用于评估GLS2表达与相关目标之间的相关性,包括免疫细胞浸润分数、TMB、MSI、错配修复(mismatch repair, MMR)基因和甲基化转移酶基因。根据样本是否配对,采用配对测试或T检验对各组之间或肿瘤与正常组织之间的GLS2表达水平进行比较。P<0.05为差异有统计学意义。所有图通过ggplot2和forestplot的R包生成。
2 结果
2.1 GLS2基因在正常组织和癌组织中的表达
利用GTEx数据库分析GLS2基因在各个正常组织中的表达水平,结果显示:GLS2在31个正常组织中的表达水平均相似,见图1A。从CCLE数据库下载各个肿瘤细胞系的数据,发现不仅GLS2表达水平相较于正常组织普遍升高,而且与正常组织中的表达范围相比,GLS2的表达范围更窄,见图1B。进一步从TCGA数据库获取GLS2在各个肿瘤样本中癌与癌旁的差异,发现GLS2在大多数癌与癌旁组织中表达差异具有统计学意义(均P<0.05),见图1C,结合TCGA及GTEx的数据,发现27种癌症中,GLS2在25种癌组织和癌旁组织的表达差异存在统计学意义(均P<0.05),见图1D。
图 1 不同组织和肿瘤中GLS2的mRNA表达水平Figure 1 mRNA expression levels of GLS2 in different tissue and tumorsA: normal mRNA expression levels of GLS2 in different tissues from GTEx database; B: mRNA expression levels of GLS2 in various tumor cell lines from CCLE database; C: mRNA expression difference of GLS2 between tumor and peri-tumor samples from TCGA database; D: mRNA expression difference of GLS2 between normal, peri-tumor and tumor samples, combining data from TCGA and GTEx databases; *: P<0.05, **: P<0.01, ***: P <0.001.2.2 GLS2基因表达水平与泛癌患者预后的相关性分析
我们利用TCGA数据库的数据评估了GLS2在不同癌症类型中与OS的相关性,其中GLS2在肾上腺皮质癌(ACC)中的风险效应最高(HR=4.37,P=0.0052),见图2。使用患者数据对每种癌症类型的中位表达值采用二分类进行生存率分析,结果显示高表达的GLS2不利于ACC患者预后(P<0.05),而在肾透明细胞癌(KIRC)、肺腺癌(LUAD)及胰腺癌(PAAD)中,高表达的GLS2与较好的预后相关(P<0.05),见图3。因此,GLS2可能作为ACC、KIRC、LUAD及PAAD患者预后的候选标志物。
考虑到随访期间,OS可能会受到非癌症相关死亡的影响,我们采用单变量Cox回归重新分析不同癌种中DSS和GLS2表达水平之间的相关性,分析结果与OS相关的结果相似,GLS2在ACC中的风险效应最高(HR=4.75,P=0.0013),见图4。在生存分析中,GLS2高表达在ACC中再次表现出更差的预后,见图5。最后,通过Cox回归分析GLS2表达和PFI之间的相关性(图略,请扫描OSID码),结果与上述OS及DSS分析结果相同。上述数据表明,GLS2的表达高低与不同肿瘤预后之间存在相关性。
2.3 GLS2基因可能通过影响不同癌种中的免疫细胞浸润来调节肿瘤微环境
为了研究GLS2基因表达与肿瘤微环境(TME)之间的关系,我们利用来自TCGA的TIMER数据库中6种免疫细胞类型(B细胞、CD4+T细胞、CD8+T细胞、中性粒细胞、巨噬细胞和树突状细胞)的浸润分数进行分析,结果发现在多种肿瘤中该基因表达与免疫细胞浸润水平存在显著相关性,其中排名前三位的肿瘤队列是低级别胶质瘤(LGG)、前列腺癌(PRAD)和肉瘤(SARC)。在LGG中,GLS2表达水平与CD8+T细胞浸润分数呈正相关,与B细胞、CD4+T细胞、中性粒细胞、巨噬细胞和树突状细胞浸润分数呈负相关。在PRAD中,GLS2表达水平与B细胞、CD4+T细胞、CD8+T细胞、中性粒细胞、巨噬细胞和树突状细胞浸润水平呈正相关,而在SARC中,GLS2表达水平与上述免疫细胞浸润水平呈负相关。上述结果表明,GLS2表达与LGG、PARP及SARC中免疫细胞的浸润水平有关,见图6。
2.4 GLS2表达与不同癌种的免疫检查点基因表达的相关性
通过mRNA序列数据库,我们评估了GLS2的表达与免疫检查点基因的表达关系。研究表明,在大部分肿瘤中,GLS2与多种免疫检查点相关基因高度相关(P<0.05),见图7。
图 7 来自TCGA数据库的多个肿瘤中GLS2 mRNA表达水平与免疫检查点mRNA表达之间的相关性Figure 7 Correlation between GLS2 mRNA expression levels and immune checkpoints’ mRNA expression in multiple tumors from the TCGA databaseThe lower triangle in each tile indicates coefficients calculated by Pearson’s correlation test, and the upper triangle indicates log10 transformed P value. *: P<0.05,**: P<0.01,***: P<0.001.2.5 GLS2的表达与某些癌症TMB和MSI相关
在32种癌症类型中,GLS2的表达与TMB之间的相关性在12种中具有显著性(P<0.05),见图8A,6种中监测到GLS2的表达与MSI具有相关性(P<0.05),见图8B。
2.6 GLS2表达与不同癌症中的MMR缺陷密切相关,并可能干扰转录后的甲基化
我们分析了GLS2表达与一些成熟的MMR基因(MLH1、MSH2、MSH6、PMS2和EPCAM)之间的关系。结果表明,在泛癌中,GLS2的表达与MMR基因的表达显著且高度相关,见图9。关于GLS2的表达与四种甲基化转移酶(DNMT1、DNMT2、DNMT3A、DNMT3B)表达的相关性,结果与上述结果一致,在特定癌症类型中存在共表达(图略,请扫描OSID码)。
3 讨论
GLS2是谷氨酰胺酶家族成员,研究发现,谷氨酰胺酶家族有助于各种人类癌症的生长。此外,许多研究表明谷氨酰胺酶家族介导的靶向治疗是有希望抑制癌症发展的分子药物[13]。我们之前的研究中,GLS2被证实为与肺腺癌生存有关的铁死亡基因。然而,GLS2对多种癌症进展的影响仍不清楚。为了系统地了解GLS2在其他癌症进展中的作用,我们对基因表达数据集进行了比较数据挖掘。
基于肿瘤的表达分析,发现GLS2在不同组织中广泛表达,在癌症组织与正常组织中的表达存在差异,其表达程度也因组织类型而异,并且该表达与这些癌症类型中的生存相关,其中ACC的OS、DSS、PFI的HR值最高,生存曲线也显示高GLS2表达患者具有较差预后。GLS2在不同类型的肿瘤中表现出不同的表达模式和功能,并根据特定的肿瘤类型抑制或促进肿瘤的发展[14]。这与我们的研究结果相似,这里发现的GLS2表达与各种组织中的癌症之间的相关性,包括与癌症生存和预后的联系,使其成为癌症监测和进一步研究的潜在生物标志物。
TME在肿瘤生长、迁徙过程中起着多方面的作用,它可以产生免疫激活或抑制效应[15-16]。为了探究GLS2与TMB的相关性,我们测试了GLS2的表达与每种癌症类型的免疫细胞浸润水平之间的相关性,发现GLS2在多种癌症(如BLCA、BRCA等)中与免疫检查点mRNA表达之间具有相关性。在Qian之前的研究中,GLS2与LUAD的预后和免疫治疗相关[17]。Zou的研究发现,谷氨酰胺酶家族在多种癌症中与肿瘤的发展和TME有关[18]。本研究与上述研究结果一致。
TMB和MSI是多种肿瘤有效的预后生物标志物和免疫治疗反应指标[19-21]。MMR为细胞内的错配修复机制,该机制的关键基因功能丧失会导致DNA复制错误无法被修复,进而导致较高的体细胞突变的产生,MMR缺陷的肿瘤对免疫治疗有较好的反应[22]。甲基化为DNA化学修饰的一种形式,能够在不改变DNA序列的前提下,改变遗传表现,控制基因表达[23]。我们研究证明了GLS2的表达与某些癌症类型的TMB和MSI之间存在联系,并且某些肿瘤也表现出GLS2与主要MMR基因和甲基化转移酶的共同表达。先前的研究将TMB和MSI与患者的药物反应相关联,尤其是针对免疫检查点抑制剂的药物,如TGF-β拮抗剂和PD-1抑制剂[24-25]。提示GLS2可以用于癌症患者给药后免疫治疗评估,可能是潜在的药物反应良好指标。
本研究中,我们使用不同的在线生物信息学平台和工具,系统分析了GLS2在各种癌症中的表达水平及其与免疫细胞浸润、TMB、MSI、MMR基因和甲基化转移酶的关系。我们的分析显示GLS2在癌症发展中发挥着不同的作用,并对癌症的临床结果进行了不同的调节。GLS2可作为肾上腺皮质癌、肾透明细胞癌和胰腺癌患者的癌症治疗靶点及预后标志物。总之,本研究结果为不同类型癌症的分子和临床特征提供了一些见解,因此可以用来帮助将基因组知识转化为癌症治疗。
然而,完全依赖于开放存取数据库中可用的信息,没有经过实验的确认,可能出现过度拟合和欠拟合,此外,GLS2在不同癌症中的具体作用分子机制和信号通路尚不清晰,因此,这项研究需要更多的理论、实验和临床研究来验证和进一步挖掘。
作者贡献郑丽华:论文撰写及修改刘峰、李卫:数据统计赵亚恒、杨艳、吴尚:收集数据刘运江:指导文章撰写 -
表 1 230例原发性乳腺浸润性导管癌患者一般资料
Table 1 General information of 230 patients with primary breast invasive ductal carcinoma
表 2 纤维蛋白原与TNM分期的关系
Table 2 Relation between fbrinogen and TNM staging
表 3 AFR二分类后(即AFR < 13.1和AFR≥13.1组)原发性乳腺浸润性导管癌患者临床特征
Table 3 Clinicopathological features of breast invasive ductal carcinoma patients after AFR second classification (AFR < 13.1 and AFR≥13.1)
表 4 乳腺浸润性导管癌患者的单因素分析
Table 4 Univariate analysis of breast invasive ductal carcinoma patients
表 5 影响乳腺浸润性导管癌患者生存的多因素分析
Table 5 Multivariate analysis of breast invasive ductal carcinoma patients
-
[1] Bugada D, Allegri M, Lavand'homme P, et al. Inflamation-based scores: a new method for patient-targeted strategies and improved perioperative outcome in cancer patients[J]. Biomed Res Int, 2014, 2014: 142425.
[2] Sun J, Chen X, Gao P, et al. Can the neutrophil to lymphocyte ratio be used to determine gastric cancer treatment outcomes? A systematic review and meta-analysis[J]. Dis Markers, 2016, 2016: 7862469. http://cn.bing.com/academic/profile?id=0735fc664d72d8ba86fe2274104680ac&encoded=0&v=paper_preview&mkt=zh-cn
[3] Templeton AJ, Ace O, McNamara MG, et al. Prognostic role of platelet to lymphocyte ratio in solid tumors: a systematic review and meta-analysis[J]. Cancer Epidemiol Biomarkers Prev, 2014, 23(7): 1204-1212. doi: 10.1158/1055-9965.EPI-14-0146
[4] Gu L, Li H, Chen L, et al. Prognostic role of lymphocyte to monocyte ratio for patients with cancer: evidence from a systematic review and meta-analysis[J]. Oncotarget, 2016, 7(22): 31926-31942. http://cn.bing.com/academic/profile?id=50d2b126f9a5cf6acb9ca93e5e651b25&encoded=0&v=paper_preview&mkt=zh-cn
[5] 樊宁波, 陈冬妮, 温浙盛, 等.术前血浆白蛋白与纤维蛋白原比值对食管鳞癌根治术后生存的预测价值[J].中华消化外科杂志, 2019, 18(6): 563-569. doi: 10.3760/cma.j.issn.1673-9752.2019.06.010 Fan NB, Chen DN, Wen ZS, et al. Predictive value of preoperative albumin-to-fibrinogen ratio for postoperative survival of esophageal squamous cell carcinoma after radical esophagectomy[J]. Zhonghua Xiao Hua Wai Ke Za Zhi, 2019, 18(6): 563-569. doi: 10.3760/cma.j.issn.1673-9752.2019.06.010
[6] Ward DG, Wei W, Buckels J, et al. Detection of pancreatic adenocarcinoma using circulating fragments of fibrinogen[J]. Eur J Gastroenterol Hepatol, 2010, 22(11): 1358-1363. doi: 10.1097/MEG.0b013e32833b6d41
[7] Bootle-Wilbraham CA, Tazzyman S, Marshall JM, et al. Fibrinogen E-fragment Inhibits the Migration and Tubule Formation of Human Dermal Microvascular Endothelial Cells in vitro[J]. Cancer Res, 2000, 60(17): 4719-4724. http://cn.bing.com/academic/profile?id=62c4d8f55e24d2fb8293ed0e27e8ff28&encoded=0&v=paper_preview&mkt=zh-cn
[8] Jensen T, Kierulf P, Sandset PM, et al. Skjonsberg, Fibrinogen and fifibrin induce synthesis of proinflflammatory cytokines from isolated peripheral blood mononuclear cells, Thromb[J]. Thromb Haemost, 2007, 97(5): 822-829. doi: 10.1160/TH07-01-0039
[9] Redman CM, Xia H. Fibrinogen biosynthesis. Assembly, intracellular degradation, and association with lipid synthesis and secretion[J]. Ann N Y Acad Sci, 2001, 936: 480-495. http://cn.bing.com/academic/profile?id=607e28eeadeb1d726cb0ac0460ed4eea&encoded=0&v=paper_preview&mkt=zh-cn
[10] Palumbo JS, Kombrinck KW, Drew AF, et al. Fibrinogen is an important determinant of the metastatic potential of circulating tumor cells[J]. Blood, 2000, 96(10): 3302-3309. doi: 10.1182/blood.V96.10.3302
[11] Palumbo JS, Talmage KE, Massari JV, et al. Platelets and fibrin(ogen) increase metastatic potential by impeding natural killer cell-mediated elimination of tumor cells[J]. Blood, 2005, 105(1): 178-185. doi: 10.1182/blood-2004-06-2272
[12] Sahni A, Francis CW. Vascular endothelial growth factor binds to fibrinogen and fibrin and stimulates endothelial cell proliferation[J]. Blood, 2000, 96(12): 3772-3778. doi: 10.1182/blood.V96.12.3772
[13] Menter DG, Tucker SC, Kopetz S, et al. Platelets and cancer: a casual or causal relationship: revisited[J]. Cancer Metastasis Rev, 2014, 33(1): 231-269. doi: 10.1007/s10555-014-9498-0
[14] 吴菲, 莫淼, 覃肖潇, 等.上海市女性乳腺癌不同筛查模式的成本效果分析[J].中华流行病学杂志, 2017, 38(12): 1665-1671. doi: 10.3760/cma.j.issn.0254-6450.2017.12.017 Wu F, Mo M, Tan XX, et al. Cost-effectiveness of multiple screening modalities on breast cancer in Chinese women from Shanghai[J]. Zhonghua Liu Xing Bing Xue Za Zhi, 2017, 38(12): 1665-1671. doi: 10.3760/cma.j.issn.0254-6450.2017.12.017
[15] Liang Y, Wang W, Que Y, et al. Prognostic value of the fibrinogen/albumin ratio (FAR) in patients with operable soft tissue sarcoma[J]. BMC Cancer, 2018, 18(1): 942. doi: 10.1186/s12885-018-4856-x
[16] You X, Zhou Q, Song J, et al. Preoperative albumin-to-fibrinogen ratio predicts severe postoperative complications in elderly gastric cancer subjects after radical laparoscopic gastrectomy[J]. BMC Cancer, 2019, 19(1): 931. doi: 10.1186/s12885-019-6143-x
[17] Wang YY, Liu ZZ, Xu D, et al. Fibrinogen-Albumin Ratio Index (FARI): A More Promising Inflammation-Based Prognostic Marker for Patients Undergoing Hepatectomy for Colorectal Liver Metastases[J]. Ann Surg Oncol, 2019, 26(11): 3682-3692. doi: 10.1245/s10434-019-07586-3
[18] Klar M, Jochmann A, Foeldi M, et al. The MSKCC nomogram for prediction the likelihood of non-sentinel node involvement in a German breast cancer population[J]. Breast Cancer Res Treat, 2008, 112(3): 523-531. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=20d4f0ad6f5aa631427488b2ae75b8f0
-
期刊类型引用(1)
1. 吴晓昀,荆林林,齐夏丽. 2型糖尿病合并乳腺癌患者预后高危因素分析. 实用癌症杂志. 2024(08): 1356-1358 . 百度学术
其他类型引用(0)