E2F转录因子家族在宫颈癌中的表达及其与临床预后的关系

摘要/Abstract

摘要：

目的研究E2F家族在宫颈癌中的表达、预后价值以及与免疫浸润之间的关系。方法采用Oncomine数据库分析E2F转录因子家族在不同类型癌症中的表达差异。采用GEPIA 2数据库分析E2F转录因子家族在宫颈癌组织与正常组织中的表达及不同临床分期间的差异。采用Kaplan-Meier plotter数据库分析E2F转录因子家族在宫颈癌中的预后价值。使用cBioPortal数据库进行E2F转录因子家族的共表达基因和遗传学突变分析。采用WebGestalt数据库对E2F家族相关基因进行功能及通路的富集分析。采用TIMER 2.0数据库分析E2F转录因子家族与宫颈癌中免疫细胞浸润的相关性。结果宫颈癌中E2F1、E2F2、E2F3、E2F7、E2F8的表达均较正常宫颈组织高，但与临床分期无明显相关性。预后分析显示，E2F1高表达和E2F3低表达与患者更长的总生存期（overall survival，OS）和无复发生存期（recurrence free survival，RFS）明显相关。基因变异分析显示，E2F1基因变异率为5%，为E2F家族中最高。计算宫颈癌E2F家族成员之间的相关性结果表明，E2F2与E2F7、E2F8的表达，E2F3与E2F5的表达，E2F7与E2F8的表达均呈正相关；E2F2与E2F5的表达呈负相关。GO分析显示，E2F基因家族的相关功能主要指向生物调节和能量代谢；KEGG分析显示，E2F基因家族主要与细胞分裂、衰老及细胞周期密切相关。TIMER2.0数据库分析显示，部分E2F家族成员的表达水平与CD4⁺T细胞、CD8⁺T细胞浸润呈正相关，与巨噬细胞浸润呈负相关。结论E2F转录因子家族在宫颈癌中存在表达失调，家族成员E2F1，2，3，7，8是未来宫颈癌精准治疗的潜在靶点，而E2F1和E2F3可能成为评价宫颈癌预后新的生物标志物。

关键词:宫颈癌,E2F转录因子家族,预后,免疫细胞浸润,生物信息学

Abstract:

ObjectiveTo explore the expression profiles and prognostic value of E2F family members and their relations with immune infiltration and relevant functions based on bioinformatics approaches.MethodsThe expression of E2F family members in different cancer types were analyzed using Oncomine database， and their expression levels in different cervical cancer clinical stages were evaluated using GEPIA 2 database. Then the prognostic value of E2Fs expression in cervical cancer was assessed using Kaplan-Meier plotter database. The co-expression profiles and genetic alterations of E2F family members were analyzed using the cBioPortal. The function and pathway enrichment analysis of E2F family were performed using WebGestalt database and the correlation analysis of E2F family and immune cell infiltration in cervical cancer was conducted using TIMER 2.0 database.ResultsThe mRNA expression levels of E2F1， E2F2， E2F3， E2F7 and E2F8 in cervical cancer were significantly increased compared with those in normal cervical tissues， but showed no significant correlations with tumor stage. Prognostic analysis showed that high expressions of E2F1 and low expressions of E2F3 were significantly associated with longer overall survival （OS） and relapse-free survival （RFS）. Analysis of the TCGA database and cBioPortal showed that the genetic alteration rate of E2F1 was 5%， which was the highest among E2F family members. GO function enrichment analysis showed that the main functions of E2F family were oriented toward the biological processes regulation and energy metabolism. KEGG pathway enrichment analysis showed that the main functions of E2Fs were related to cell division， cell cycle and cellular senescence. The TIMER 2.0 database analysis showed that the expressions of most E2F family members were significantly positively correlated with the infiltration of CD4+ T cells， and negatively correlated with the infiltration of CD8+ T cells and macrophages.ConclusionThe molecular and expression profiles of E2Fs members are dysregulated in cervical cancer， and E2F family members （especially E2F1 and E2F3） are potential prognostic markers and new therapeutic targets， which may provide a new avenue for studying the mechanism and targeted therapy of cervical cancer.

Key words:cervical cancer,E2F transcription factors,prognosis,immune cell infiltrates,bioinformatics

陈玥, 佘婧瑶, 陈思, 曹莹, 梁春云, 卢燕, 王佩娟, 公真. E2F转录因子家族在宫颈癌中的表达及其与临床预后的关系[J]. betway必威登陆网址（betway.com ）学报, 2023, 44(3): 215-223.

Yue CHEN, Jingyao SHE, Si CHEN, Ying Cao, Chunyun LIANG, Yan LU, Peijuan WANG, Zhen Gong. The expression for E2Fs in human cervical cancer and its relationship with clinical prognosis[J]. Journal of Shandong First Medical Unversity & Shandong Academy of Medical Sciences, 2023, 44(3): 215-223.

图/表8

图1E2F转录因子家族表达分析A：E2F家族在不同类型肿瘤与正常组织中的表达差异（Oncomine）；B：E2F家族在宫颈癌及正常宫颈组织中的表达差异（GEPIA2）

图2宫颈癌中E2F家族的相对表达水平（GEPIA2）T为肿瘤组织；N为正常组织。

表1E2F家族在不同类型宫颈癌与正常组织中表达水平的倍数变化（Oncomine）

转录因子	宫颈癌类型	倍数变化	t	P	数据来源
E2F1	宫颈鳞状细胞癌	3.640	7.886	< 0.001	Scotto Cervix 2 Statistics
E2F2	宫颈鳞状细胞癌	2.190	8.117	< 0.001	Biewenga Cervix Statistics
E2F3	宫颈鳞状细胞癌	3.296	9.625	< 0.001	Zhai Cervix Statistics
E2F3	宫颈腺癌	4.546	10.055	< 0.001	Pyeon Multi-cancer Statistics
E2F7	宫颈鳞状细胞癌	3.252	6.950	< 0.001	Scotto Cervix 2 Statistics
E2F7	宫颈腺癌	9.390	12.346	< 0.001	Pyeon Multi-cancer Statistics
E2F8	宫颈鳞状细胞癌	2.557	4.719	< 0.001	Zhai Cervix Statistics
E2F8	宫颈鳞状细胞癌	2.751	4.646	< 0.001	Scotto Cervix 2 Statistics

表1E2F家族在不同类型宫颈癌与正常组织中表达水平的倍数变化（Oncomine）

转录因子	宫颈癌类型	倍数变化	t	P	数据来源
E2F1	宫颈鳞状细胞癌	3.640	7.886	< 0.001	Scotto Cervix 2 Statistics
E2F2	宫颈鳞状细胞癌	2.190	8.117	< 0.001	Biewenga Cervix Statistics
E2F3	宫颈鳞状细胞癌	3.296	9.625	< 0.001	Zhai Cervix Statistics
E2F3	宫颈腺癌	4.546	10.055	< 0.001	Pyeon Multi-cancer Statistics
E2F7	宫颈鳞状细胞癌	3.252	6.950	< 0.001	Scotto Cervix 2 Statistics
E2F7	宫颈腺癌	9.390	12.346	< 0.001	Pyeon Multi-cancer Statistics
E2F8	宫颈鳞状细胞癌	2.557	4.719	< 0.001	Zhai Cervix Statistics
E2F8	宫颈鳞状细胞癌	2.751	4.646	< 0.001	Scotto Cervix 2 Statistics

图3E2F家族表达水平与宫颈癌临床分期的相关性（GEPIA2）

图4E2F家族与宫颈癌患者预后的关系（Kaplan-Meier plotter）

图5宫颈癌中各E2F家族成员的基因改变频率（cBioPortal）

图6E2F家族各成员之间的相关性（GEPIA2）

图7E2F家族相关分子分析（cBioPortal，Cytoscape）A为宫颈癌中190个E2F家族相关分子的蛋白相互作用网络图；B为HUB基因鉴定。

参考文献31

1	Sung H， Ferlay J， Siegel RL， et al. Global cancer statistics 2020： GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries［J］. CA Cancer J Clin， 2021， 71（3）： 209.
2	Bray F， Lortet-Tieulent J， Znaor A， et al. Patterns and trends in human papillomavirus-related diseases in Central and Eastern Europe and Central Asiaa［J］. Vaccine， 2013， 31 ： H32.
3	Lees BF， Erickson BK， Huh WK. Cervical cancer screening： evidence behind the guidelines［J］. Am J Obstet Gynecol， 2016， 214（4）： 438.
4	Wang W， Jiao Y， Zhang L， et al. Multiparametric MRI-based radiomics analysis： differentiation of subtypes of cervical cancer in the early stage［J］. Acta radiol， 2022， 63（6）： 847.
5	Chen HZ， Ouseph MM， Li J， et al. Canonical and atypical E2Fs regulate the mammalian endocycle［J］. Nat Cell Biol， 2012， 14（11）： 1192.
6	Rowland BD， Bernards R. Re-evaluating cell-cycle regulation by E2Fs［J］. Cell， 2006， 127（5）： 871.
7	Lammens T， Li J， Leone G， et al. Atypical E2Fs： new players in the E2F transcription factor family［J］. Trends Cell Biol， 2009， 19（3）： 111.
8	Zhang C， Cui Y， Wang G， et al. Comprehensive analysis of the expression and prognosis for E2Fs in human clear cell renal cell carcinoma［J］. J Healthc Eng， 2021， 2021： 5790416.
9	Liu J， Xia L， Wang S， et al.E2F4promotes the proliferation of hepatocellular carcinoma cells through upregulation ofCDCA3［J］. J Cancer， 2021， 12（17）： 5173.
10	Hu D， Meng N， Lou X， et al. Prognostic values of E2F1/2 transcriptional expressions in chromophobe renal cell carcinoma patients： evidence from bioinformatics analysis［J］. Int J Gen Med， 2021， 14： 3593.
11	Cheng Y， Han X， Mo F， et al. Apigenin inhibits the growth of colorectal cancer through down-regulation of E2F1/3 by miRNA-215-5p［J］. Phytomedicine， 2021， 89： 153603.
12	Emmrich S， Pützer BM. Checks and balances： E2F-microRNA crosstalk in cancer control［J］. Cell Cycle， 2010， 9（13）： 2555.
13	Yang C， Zhang ZC， Liu TB， et al. E2F1/2/7/8 as independent indicators of survival in patients with cervical squamous cell carcinoma［J］. Cancer Cell Int， 2020， 20： 500.
14	Rhodes DR， Kalyana-Sundaram S， Mahavisno V， et al. Oncomine 3.0： genes， pathways， and networks in a collection of 18，000 cancer gene expression profiles［J］. Neoplasia， 2007， 9（2）： 166.
15	Tang Z， Kang B， Li C， et al. GEPIA2： an enhanced web server for large-scale expression profiling and interactive analysis［J］. Nucleic Acids Res， 2019， 47（W1）： W556.
16	Nagy Á， Lánczky A， Menyhárt O， et al. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets［J］. Sci Rep， 2018， 8（1）： 9227.
17	Gao J， Aksoy BA， Dogrusoz U， et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal［J］. Sci Signal， 2013， 6（269）： pl1.
18	Szklarczyk D， Gable AL， Nastou KC， et al. The STRING database in 2021： customizable protein-protein networks， and functional characterization of user-uploaded gene/measurement sets［J］. Nucleic Acids Res， 2021， 49（D1）： D605.
19	Wang J， Vasaikar S， Shi Z， et al. WebGestalt 2017： a more comprehensive， powerful， flexible and interactive gene set enrichment analysis toolkit［J］. Nucleic Acids Res， 2017， 45（W1）： W130.
20	Li T， Fu J， Zeng Z， et al. TIMER2.0 for analysis of tumor-infiltrating immune cells［J］. Nucleic Acids Res， 2020， 48（W1）： W509.
21	Minnaar CA， Maposa I， Kotzen JA， et al. Effects of modulated electro-hyperthermia （mEHT） on two and three year survival of locally advanced cervical cancer patients［J］. Cancers （Basel）， 2022， 14（3）： 656.
22	Xu J， Liu H， Yang Y， et al. Genome-wide profiling of cervical RNA-binding proteins identifies human papillomavirus regulation of RNASEH2A expression by viral E7 and E2F1［J］. mBio， 2019， 10（1）： e02687.
23	Liu Z， Mo H， Liu R， et al. Matrix stiffness modulates hepatic stellate cell activation into tumor-promoting myofibroblasts via E2F3-dependent signaling and regulates malignant progression［J］. Cell Death Dis， 2021， 12（12）： 1134.
24	Wu X， Kong C， Wu Y. Long intergenic non-protein coding RNA 1094 （LINC01094） promotes the progression of breast cancer （BC） by regulating the microRNA-340-5p （miR-340-5p）/E2F transcription factor 3 （E2F3） axis［J］. Bioengineered， 2021， 12（1）： 9046.
25	Zong S， Liu X， Zhou N， et al. E2F7， EREG， miR-451a and miR-106b-5p are associated with the cervical cancer development［J］. Arch Gynecol Obstet， 2019， 299（4）： 1089.
26	Ye L， Guo L， He Z， et al. Upregulation of E2F8 promotes cell proliferation and tumorigenicity in breast cancer by modulating G1/S phase transition［J］. Oncotarget， 2016， 7（17）： 23757.
27	Kim LK， Park SA， Eoh KJ， et al. E2F8 regulates the proliferation and invasion through epithelial-mesenchymal transition in cervical cancer［J］. Int J Biol Sci， 2020， 16（2）： 320.
28	Peng H， Shen J， Long X， et al. Local release of TGF-β inhibitor modulates tumor-associated neutrophils and enhances pancreatic cancer response to combined irreversible electroporation and immunotherapy［J］. Adv Sci （Weinh）， 2022， 9（10）： e2105240.
29	Zhang S， Xiong H， Yang J， et al. Pan-cancer analysis reveals the multidimensional expression and prognostic and immunologic roles ofVSTM2Lin cancer［J］. Front Mol Biosci， 2022， 8： 792154.
30	Waldman AD， Fritz JM， Lenardo MJ. A guide to cancer immunotherapy： from T cell basic science to clinical practice［J］. Nat Rev Immunol， 2020， 20（11）： 651.
31	Liu J， Yang J， Gao F， et al. A microRNA-messenger RNA regulatory network and its prognostic value in cervical cancer［J］. DNA Cell Biol， 2020， 39（7）： 1328.

[1]	戚孟琪, 吕玉芹, 高蕾, 张敬军.基于生物信息学方法探究布瓦西坦抗癫痫机制[J]. betway必威登陆网址（betway.com ）学报, 2023, 44(6): 401-407.
[2]	周海英, 郭冉, 赵聪, 吴培培, 陈卫东, 孙宗文.血清NLR、PLR与非小细胞肺癌预后相关性研究进展[J]. betway必威登陆网址（betway.com ）学报, 2023, 44(2): 147-151.
[3]	苏颖, 姚海平, 刘佩.手术治疗新生儿先天性膈疝临床疗效及预后的相关预测因素[J]. betway必威登陆网址（betway.com ）学报, 2023, 44(2): 121-124.
[4]	王朝亚, 潘淳, 贾犇黎, 葛亮.RCOR3基因在多形性胶质母细胞瘤中的表达和预后[J]. betway必威登陆网址（betway.com ）学报, 2023, 44(2): 93-100.
[5]	黄永胜, 黄彩娜, 董学岭, 卓秀丽, 房娟娟, 宋文霞, 张雨露, 阎磊, 陈刚, 吕仁广.膀胱尿路上皮癌蛋白组学相关个体化预后特征的推导与验证[J]. betway必威登陆网址（betway.com ）学报, 2023, 44(1): 15-23.
[6]	徐汝飞.高压氧辅助神经内镜经侧裂入路治疗高血压脑出血的临床疗效[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(7): 525-528.
[7]	陈荣华, 王艳梅, 石瑞平.EZH2在AML1/ETO阳性急性髓系白血病中的表达及临床意义[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(5): 360-364.
[8]	徐雯, 吕玉芹, 高蕾, 张敬军.基于多组学数据研究奥卡西平抗痫机制[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(5): 327-334.
[9]	崔杨, 金洋, 刘玥.二元应对护理干预在宫颈癌放化疗患者及配偶中的应用效果[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(11): 839-843.
[10]	卜文超, 曾宪智, 关云茜, 谢思源, 陈世新, 汤挺兵, 曹明国.SOCS基因家族在肝细胞癌中的表达及生物信息学分析[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(10): 740-747.
[11]	安鑫.应激性高血糖对非糖尿病急性心肌梗死患者PCI术后心肌灌注、预后的影响[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(1): 41-44.
[12]	接智慧, 周岩, 杨海燕, 杨洪川, 程俊兰, 许文.miR-17-92基因簇在子宫内膜病变组织中的表达及临床意义[J]. betway必威登陆网址（betway.com ）学报, 2022, 43(1): 16-19.
[13]	吕玉芹, 高蕾, 徐雯, 戚孟琪, 张敬军.颞叶癫痫差异发作频率相关生物通路及蛋白-蛋白相互作用网络分析[J]. betway必威登陆网址（betway.com ）学报, 2021, 42(7): 508-515.
[14]	戴文玲.急性胰腺炎患者血清Ghrelin及NGAL水平变化与临床预后的关系[J]. betway必威登陆网址（betway.com ）学报, 2021, 42(6): 454-456.
[15]	陈蕊, 张俊勇, 赵琪.炎症标志物与消化系统恶性肿瘤相关性的研究进展[J]. betway必威登陆网址（betway.com ）学报, 2021, 42(4): 317-320.