Skip to main content
Download PDF

Evaluation of clinical value and potential mechanism of MTFR2 in lung adenocarcinoma via bioinformatics

BMC Cancer volume 21, Article number: 619 (2021) Cite this article

Abstract

Background

Mitochondrial fission regulator 2 (MTFR2) was involved in the progression and development of various cancers. However, the relationship between MTFR2 with lung adenocarcinoma (LUAD) had not been reported. Herein, this study analyzed the clinical significance and potential mechanisms of MTFR2 in LUAD via bioinformatics tools.

Results

We found that the level of MTFR2 was increased, and correlated with sex, age, smoking history, neoplasm staging, histological subtype and TP53 mutation status in LUAD patients. Kaplan-Meier survival analysis showed LUAD patients with increased MTFR2 had a poor prognosis. In addition, univariate COX regression analysis showed neoplasm staging, T stage, distant metastasis and MTFR2 level were risk factors for the prognosis of LUAD. A total of 1127 genes were coexpressed with MTFR2, including 840 positive and 208 negative related genes. KEGG and GSEA found that MTFR2 participated in the progression of LUAD by affecting cell cycle, DNA replication, homologous recombination, p53 signaling pathway and other mechanisms. The top 10 coexpressed genes, namely CDK1, CDC20, CCNB1, PLK1, CCNA2, AURKB, CCNB2, BUB1B, MAD2L1 and BUB1 were highly expressed, and were associated with poor prognosis in LUAD.

Conclusions

Consequently, we elucidated MTFR2 was a biomarker for diagnosis and poor prognosis in LUAD, and might participate in the progression of LUAD via affecting cell cycle, DNA replication, homologous recombination and p53 signaling pathway.

Peer Review reports

Introduction

Non-small cell lung cancer (NSCLC) was one of the most common malignant tumors in the world, and lung adenocarcinoma (LUAD) was one of the common subtype of NSCLC [1, 2]. LUAD patients in early stage got improved long-term prognosis by surgery and neoadjuvant chemotherapy, while due to incomplete excision, primary and secondary drug resistance and other reasons, recurrence and metastasis happened. The 5-year survival rate of the patients was still very low, and its treatment was extremely challenging [2, 3]. In recent years, targeted therapy had greatly improved the prognosis of patients with LUAD, but the mortality of patients with LUAD remained at a high level [4,5,6]. Therefore, it was of great significance to further discover the novel molecules and explore their regulatory mechanisms involved in the occurrence and development of LUAD resulting in finding novel therapeutic targets.

Mitochondrial fission regulator 2 (MTFR2), also known as FAM54A, was located on the 6q23.2 chromosome and belonged to the MTFR family. Abnormal mitosis of mitochondria was related to the pathogenesis of many diseases, as well as occurrence and prognosis of several tumors [7,8,9,10,11]. MTFR2 was up-regulated in breast cancer and was related to age, lymph node metastasis, and prognosis in HER2 positive breast cancer patients. Knockout of MTFR2 prohibited the proliferation, migration and invasion of breast cancer cells, and blocked epithelial-stromal transformation (EMT) [8]. MTFR2 affected the transcription of TTK by activating the transcriptional promoter. MTFR2-dependent TTK regulation played a key role in maintaining GSC (Glioma stem-like cells) in glioma, and may be a target molecule of new drugs for glioma patients [12]. Whereas, at present, rare researches were focused on MTFR2 in other tumors including LUAD, except for those cancers listed above. Herein, the present study explored the clinical value and potential regulatory mechanism of MTFR2 in the progression of LUAD through multiple databases.

Results

MTFR2 was upregulated and correlated with sex, age, smoking history, cancer stage, histological subtype and TP53 mutation status in LUAD patients

MTFR2 was abnormally expressed in pan-cancerous tissues, and was mainly over-expressed in a variety of tumor tissues (Fig. 1). In detail, in the Oncomine database, MTFR2 was highly expressed in breast cancer, colorectal cancer, gastric cancer, lung cancer and other tumor tissues, while lowly- expressed in brain nerve tumor, breast cancer, leukemia and other tumor tissues (Fig. 1A). In the Timer database, MTFR2 was highly expressed in BLCA, BRCA, CHOL, LUSC, LUAD and other tumor tissues (Fig. 1B).

Fig. 1
figure 1

MTFR2 expressions in cancer tissues. A Oncomine database; B Timer database. Note: A Red, increased expression of MTFR2; blue, decreased expression of MTFR2; B red, cancer tissues; blue, normal tissues. *, p < 0.05; **, p < 0.01; ***, p < 0.001

Full size image

In addition, we found consistent results in the Ualcan database, where MTFR2 was highly expressed in LUAD tissues (Table 1). In addition, in the TCGA database, the expression of MTFR2 was increased in unpaired LUAD patients (Fig. 2A) and in 57 pairs of LUAD patients (Fig. 2B).

Table 1 The expression of MTFR2 in Ualcan database was correlated to the clinicopathological characteristics in LUAD patients
Full size table
Fig. 2
figure 2

The expression of MTFR2 was correlated to diagnosis and prognosis in LUAD patients. A Expression of MTFR2 in unpaired LUAD tissues; B Expression of MTFR2 in 57 paired LUAD tissues; C The diagnostic value of MTFR2 in unpaired LUAD tissues; D The diagnostic value of MTFR2 in 57 paired LUAD tissues; E The expression of MTFR2 in TCGA database was correlated with the OS of LUAD patients; F-G The levels of MTFR2 (228069_at) and MTFR2 (234944_at) were correlated to the OS of LUAD patients with Kaplan-Meier Plotter analysis; H-I The levels of MTFR2 (228069_at) and MTFR2 (234944_at) were correlated to the early progression of LUAD patients with Kaplan-Meier Plotter analysis. Note: Normal, normal lung tissues; LUAD, lung adenocarcinoma tissues; AUC, Area under curve

Full size image

In the Ualcan database, we found that the expression of MTFR2 was related to the clinicopathological features of LUAD patients (Table 1). In detail, the level of MTFR2 was related to the sex (Male vs Female), Age (Age(41-60Yrs) vs Age (61-80Yrs)), Smoking (Non smoker vs Smoker; Non smoker vs Reformed smoker2; Smoker vs Reformed smoker1; Smoker vs Reformed smoker2; Reformed smoker1 vs Reformed smoker2), histological subtype (NOS vs Mixed; NOS vs LBC-Non Mucinous; NOS vs Solid Pattern Predominant; NOS vs Acinar; NOS vs LBC-Mucinous; NOS vs Papillary; Mixed vs Solid Pattern Predominant; Mixed vs LBC-Mucinous; Mixed vs Papillary; Clear Cell vs LBC-Non Mucinous; Clear Cell vs Papillary; LBC-Non Mucinous vs Solid Pattern Predominant; LBC-Non Mucinous vs LBC-Mucinous; Solid Pattern Predominant vs LBC-Mucinous; Solid Pattern Predominant vs Mucinous carcinoma; Solid Pattern Predominant vs Papillary; Acinar vs LBC-Mucinous; LBC-Mucinous vs Papillary) and TP53 mutation (TP53-Mutant vs TP53-Non Mutant) of LUAD patients.

In TCGA database, the area under the curve (AUC) showed that the level of MTFR2 had a diagnostic value in LUAD patients via the ROC analysis (Fig. 2C and D). Further Kaplan-Meier analysis showed that the LUAD patients with increased MTFR2 had a poor prognosis (Fig. 2E). In addition, in the Kaplan-Meier Plotter database, the levels of MTFR2 (228069_at) and MTFR2 (234944_at) were associated with the overall survival and early disease progression of LUAD patients (Fig. 2F-I). Univariate Cox analysis of complete clinical data of LUAD patients in TCGA database showed that cancer stage, T stage, lymphatic metastasis as well as MTFR2 expression were risk factors for the prognosis of LUAD patients (Table 2). Further multivariate Cox analysis showed that cancer stage was an independent risk factor for the prognosis of LUAD patients (Table 2).

Table 2 COX regression analysis of the correlation between clinicopathological features with prognosis
Full size table

Screening of MTFR2 co-expressed genes

We screened 1127 MTFR2 co-expressed genes in TCGA database by screening conditions, most of which were positive related genes. Among them, there were 840 positive and 208 negative related genes (S1 Table). The top 10 genes that were positively and negatively related to MTFR2 (Table 3). The expression of the top 10 positive and negative related genes of MTFR2 in LUAD tissues from TCGA were abnormally expressed (Fig. 3).

Table 3 The top 10 genes with positive and negative correlation to MTFR2
Full size table
Fig. 3
figure 3

The expressions of the top 10 positive and negative related genes of MTFR2 in LUAD tissues. A The top 10 positively correlated genes; B The top 10 negatively correlated genes. Note: ***, p < 0.001; N, normal lung tissue; T, lung adenocarcinoma tissue

Full size image

Analysis of MTFR2 function and its coexpressed genes by GO, KEGG and GSEA

In order to further understand the potential function of MTFR2 in the development of LUAD, we carried out GO and KEGG analysis of MTFR2 co-expressed genes. MTFR2 coexpressed genes were involved in BP, CC and MF (Fig. 4A-C), and in details, were mainly involved in DNA replication, chromosome segregation, mitotic nuclear division, mitotic cell cycle phase transition and other processes (S2 Table). KEGG pathway analysis indicated MTFR2 co-expressed genes were mainly involved in cell cycle, DNA replication, proteasome, homologous recombination, spliceosome, nucleotide excision repair, human T-cell leukemia virus 1 infection, p53 signaling pathway and other signal pathways (Fig. 4D and Table 4). What’s more, in GSEA database, we found that cell cycle, DNA_replication, homologous recombination, p53 signal pathway, oocyte meiosis and base_excision repair signaling pathways were highly enriched in MTFR2 highly expressed group (Fig. 5 and Table 5). To sum up, these results suggested that MTFR2 might regulate the progress of LUAD through cell cycle, DNA replication, homologous recombination and p53 signaling pathway.

Fig. 4
figure 4

Analysis of MTFR2 co-expressed genes by GO and KEGG. A BP; B CC; C MF; D KEGG

Full size image
Table 4 MTFR2 co-expressed genes by KEGG
Full size table
Fig. 5
figure 5

GSEA displayed the main signal pathways derived from high expression of MTFR2 enrichment. A Cell cycle; B DNA replication; C Homologous recombination; D p53 signaling pathway; E Oocyte meiosis; F Base_excision_ repair

Full size image
Table 5 GSEA displayed the main signal pathways derived from high expression of MTFR2 enrichment
Full size table

Construction of PPI network and the expressions of hub genes and their prognostic value

By identifying the function of MTFR2 co-expressed genes, the potential biological function of MTFR2 was also inferred. The top 10 Hub genes in the PPI network were selected by the CytoHubba plug-in as CDK1, CDC20, CCNB1, PLK1, CCNA2, AURKB, CCNB2, BUB1B, MAD2L1 and BUB1 (S1 Fig and Table 6). The expression of MTFR2 was correlated with that of CDK1, CDC20, CCNB1, PLK1, CCNA2, AURKB, CCNB2, BUB1B, MAD2L1 and BUB1 (Fig. 6). In addition, GEPIA database was used to analyze the expressions of Hub genes and their prognostic values. We found that the expressions of CDK1, CDC20, CCNB1, PLK1, CCNA2, AURKB, CCNB2, BUB1B, MAD2L1 and BUB1 genes were significantly increased in LUAD tissues (Fig. 7). The top10 Hub genes were associated with OS of LUAD patients, and CDK1, CCNB1, PLK1, AURKB, CCNB2, BUB1B and BUB1 were associated with DFS in LUAD patients with disease-free progression (Fig. 8).

Table 6 The top 10 Hub genes in PPI Network
Full size table
Fig. 6
figure 6

The expression of MTFR2 was correlated to the top 10 Hub genes. A CDK1; B CDC20; C CCNB1; D PLK1; E CCNA2; F AURKB; G CCNB2; H BUB1B; I MAD2L1; J BUB1

Full size image
Fig. 7
figure 7

The expressions of the top10 Hub genes in LUAD tissues from TCGA Database. A CDK1; B CDC20; C CCNB1; D PLK1; E CCNA2; F AURKB; G CCNB2; H BUB1B; I MAD2L1; (J) BUB1. Note: T, LUAD tissues; N, normal tissues; red, LUAD tissues; gray, normal lung tissues; *, p < 0.05

Full size image
Fig. 8
figure 8

The prognosis values of Hub genes in LUAD from GEPIA. A CDK1; B CDC20; C CCNB1; D PLK1; E CCNA2; F AURKB; G CCNB2; H BUB1B; I MAD2L1; J BUB1

Full size image

Discussion

Mitochondria, via regulating the membrane potential, participated in programmed cell death as well as cell proliferation [13, 14]. MTFR2 played an vital role in promoting the division of mitochondria in cells. Under pathological condition, abnormal mitosis of mitochondria was related to the development of tumors [10, 11]. At present, studies found that the expression of MTFR2 in breast cancer tissues were increased, and was correlated to the clinicopathological features and poor prognosis of patients. Interfering with MTFR2 inhibited the growth and migration of breast cancer cells, indicating its biological role as oncogene and with prognostic value [8]. However, the role of MTFR2 in LUAD was still not been elucidated. In present study, it was found for the first time that MTFR2 was highly expressed in LUAD tissues in multiple databases, and the level of MTFR2 was related to sex, age, smoking history, cancer stage, histological subtypes and TP53 mutation status in LUAD patients. ROC analysis showed that the expression of MTFR2 had diagnostic value in LUAD. Kaplan-Meier analysis showed that the overall survival time of LUAD patients with increased MTFR2 expression was related to the early disease progression. Univariate Cox analysis showed that cancer stage, T stage, distant metastasis and the level of MTFR2 were risk factors for the prognosis in LUAD. These results indicated that MTFR2 played an important role in the development of LUAD.

Abnormal cell cycle leaded to abnormal regulation of cell growth and induced overgrowth of cells, resulting in poor prognosis of tumor patients. Therefore, through a detailed understanding of the mechanism of tumor growth and inhibition of tumor occurrence and development, the prognosis of patients could be improved [15,16,17]. Nichols et al. reported that high mutation burden in LUAD patients with increased HORMAD1 expression decreased the survival rate of LUAD patients. HORMAD1 was one of the key genes in the process of homologous recombination (HR), and promoted the formation of RAD51 filaments to participate in the process of homologous recombination. Interfering with HORMAD1 reduced the growth of tumor cells [18]. FAM111B was highly expressed in LUAD tissues. LUAD patients with high expression of FAD111B had low relapse-free survival (RFS) and overall survival (OS). Interfering with FAM111B inhibited the proliferation, cell cycle arrest and migration of LUAD cells and the ability of tumor formation in nude mice, but induced apoptosis. FAM111B might regulate the growth of LUAD through p53 signaling pathway, inhibited the expression of BAG3 and BCL2, and played an important role in cell cycle and apoptosis [19]. In present study, we found that MTFR2 coexpressed genes were mainly involved in cell cycle, DNA replication, homologous recombination, p53 signaling pathway and other signaling pathways by KEGG. Secondly, in GSEA, we also found that cell cycle, DNA replication, homologous recombination and p53 signaling pathways were highly enriched in MTFR2 highly-expressed group. Up to now, it was found that interfering with MTFR2 affected the process of EMT [8], but the role of MTFR2 in LUAD was still not been elucidated. However, MTFR2 promoted mitochondrial division in cells, which was related to tumor progression. Therefore, the interaction between MTFR2 with cell cycle, DNA replication, homologous recombination and p53 signaling pathway was worth in-depth study in LUAD.

Differential expression of genes in the process of tumorigenesis and development regulated the progression of cancers and affected the prognosis [20,21,22,23]. Polo-like kinase1 (PLK1) was a key regulatory factor in cell cycle and apoptosis. PLK1 was overexpressed in many kinds of Cutaneous T-cell lymphomas (CTCL) cells. Down-regulation of PLK1 inhibited cell growth, stability and proliferation, and also leaded to cell cycle arrest, mitotic protein changes, apoptosis and mitotic abnormalities [20]. CCNA2 was highly expressed, and was correlated to poor prognosis in NSCLC patients [22]. Topoisomerase IIA (TOP2A) was highly expressed, and was associated with poor prognosis in LUAD, indicating that TOP2A was a prognostic marker. Knockdown of TOP2A in A549 and GLC82 cells inhibited cell proliferation, migration and invasion, and decreased the expression of CCNB1 and CCNB2, indicating that it promoted the progress of LUAD by targeting CCNB1 and CCNB2 [23]. In this study, we found that the hub genes were involved in a variety of signal pathways including those in the occurrence and development of tumors. PPI network analyzed by GEPIA showed that the expressions of CDK1, CDC20, CCNB1, PLK1, CCNA2, AURKB, CCNB2, BUB1B, MAD2L1 and BUB1 were increased and related to the overall survival (OS), and CDK1, CCNB1, PLK1, AURKB, CCNB2, BUB1B and BUB1 were associated with disease-free survival (DFS) in LUAD patients. The expression of MTFR2 was also closely related to that of CDK1, CDC20, CCNB1, PLK1, CCNA2, AURKB, CCNB2, BUB1B, MAD2L1 and BUB1, which further indicated that MTFR2 and Hub genes played an vital role in the diagnosis, treatment and disease progression of LUAD.

Generally speaking, we found that the high level of MTFR2 was related to the clinicopathological features, diagnosis and prognosis, which might be used as a potential target for diagnosis and prognosis, and MTFR2 promoted the occurrence and development through cell cycle, DNA replication, homologous recombination and p53 signal pathway in LUAD.

Materials and methods

Oncomine and Timer database

The level of MTFR2 in pan-cancer tissues was analyzed by Oncomine database (https:// www.oncomine.org/resource/login) [24] and Timer database (https://cistrome.shinyapps.io/timer/) [25]. Screening criteria for Oncomine database: (1) Gene: FAM54A; (2) Analysis Type: Cancer vs Normal; (3) Date Type: mRNA; (4) P < 0.05; (5) Fold Chang ≥1.5. Search for FAM54A in the gene module of Timer database to analyze its level in pan-cancer tissues.

Ualcan database

Ualcan database (http://ualcan.path.uab.edu) [26] was used to explore the expression of MTFR2, and the correlation between the level of MTFR2 with the clinicopathological parameters (including sex, age, race, smoking history, stage, histological subtype, etc.) in LUAD.

TCGA and Kaplan-Meier plotter database

The gene expression data of 594 cases of LUAD with HTSeq-FPKM were downloaded from TCGA (https://portal.gdc.cancer.gov/projects/) [27], including 59 cases of normal lung tissues and 535 cases of LUAD tissues. Among them, the gene expression data of 57 normal lung tissues and 57 LUAD tissues came from the same patients. TCGA data were used to analyze the expression of MTFR2 and its diagnostic value via the ROC analysis.

The clinical data of 522 patients with LUAD were downloaded and screened. The patients with unknown or incomplete clinical characteristic information and lack of prognostic follow-up data were excluded, and the rest were retained for the survival analysis, univariate Cox analysis and multivariate Cox analysis. The mRNA data for Kaplan-Meier Plotter (http://kmplot.com/analysis/) analysis [28] were derived from TCGA and GEO databases. According to the median value of MTFR2 mRNA level, the data were divided into two groups: high and low expression groups to explore the correlation with prognosis.

Screening of MTFR2 co-expressed genes

Pearson coefficient (r) was applied to indicate the correlation between genes and reflect the biological relationship between the two genes. MTFR2 co-expressed genes were screened in 535 cases of LUAD tissues from TCGA via the R (version 3.6.1) [29]. Screening criteria: P < 0.001 and Pearson coefficient (r > 0.4 or r < − 0.4), were defined as moderate or above.

GO, KEGG and GSEA

The co-expressed genes of MTFR2 were analyzed by GO and KEGG [30] using R clusterProfiler package to explore biological function and relevant signal pathways in LUAD [27]. The gene expression data of 535 patients with LUAD from TCGA were divided into high and low expression groups according to the median MTFR2 level. The genes were arranged according to the degree of differential expression by GSEA (version 4.0.1) to explore the effect of MTFR2 on each gene. Each analysis was repeated 1000 times [29, 31]. Screening criteria: NOM p < 0.05.

PPI network construction

STRING (https://string-db.org/) [32] database was used to analyze the protein-protein interaction (PPI) network. Screening condition: combined score > 0.9 was considered to be statistically significant. The obtained PPI network was introduced into Cytoscape 3.6.1 software, and the top 10 genes with high connectivity were screened by CytoHubba plug-in, and were defined as Hub genes [33].

GEPIA database

The mRNA data from TCGA and GTEx, and data of normal lung tissues from GEPIA database were selected to verify the expression of Hub genes in LUAD compared with normal lung tissues. According to the median value of Hub gene expressions, the data were divided into two groups: high and low expression groups, and the relationships between their levels and OS and DFS were analyzed.

Availability of data and materials

Publicly available datasets were analyzed in this study. The datasets generated for this study could be found in TCGA (https://portal.gdc.cancer.gov/projects/TCGA-LUAD (HTSeq-FPKM)).

Abbreviations

BP:

Biological processes

CC:

Cellular component

MF:

Molecular function

OS:

Overall survival

DFS:

Disease-free survival

GEPIA:

Gene Expression Profiling Interactive Analysis

MTFR2:

Mitochondrial fission regulator 2

LUAD:

Lung adenocarcinoma

NSCLC:

Non-small cell lung cancer

EMT:

Epithelial-stromal transformation

GSC:

Glioma stem-like cells

RPA1:

Replication protein A1

HR:

Homologous recombination

RFS:

Relapse-free survival

PLK1:

Polo-likekinase1

CTCL:

Cutaneous T-cell lymphomas

CDK:

Cyclin-dependent kinase

References

  1. Song Y, Chen D, Zhang X, Luo Y, Li S. Integrating genetic mutations and expression profiles for survival prediction of lung adenocarcinoma. Thoracic Cancer. 2019;10(5):1220–8. https://doi.org/10.1111/1759-7714.13072.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Chen QY, Zhong Q, Zhou JF, Qiu XT, Dang XY, Cai LS, et al. Conditional survival and recurrence of remnant gastric cancer after surgical resection: a multi-institutional study. Cancer Sci. 2020;111(2):502–12. https://doi.org/10.1111/cas.14231.

    Article  CAS  PubMed  Google Scholar 

  3. Mizutani K, Guo X, Shioya A, Zhang J, Zheng J, Kurose N, et al. The impact of PRDX4 and the EGFR mutation status on cellular proliferation in lung adenocarcinoma. Int J Med Sci. 2019;16(9):1199–206. https://doi.org/10.7150/ijms.36071.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Okita R, Maeda A, Shimizu K, Nojima Y, Saisho S, Nakata M. PD-L1 overexpression is partially regulated by EGFR/HER2 signaling and associated with poor prognosis in patients with non-small-cell lung cancer. Cancer Immunol Immunother. 2017;66(7):865–76. https://doi.org/10.1007/s00262-017-1986-y.

    Article  CAS  PubMed  Google Scholar 

  5. Khunger M, Jain P, Rakshit S, Pasupuleti V, Hernandez AV, Stevenson J, et al. Safety and efficacy of PD-1/PD-L1 inhibitors in treatment-naive and chemotherapy-refractory patients with non-small-cell lung Cancer: a systematic review and meta-analysis. Clin Lung Cancer. 2018;19(3):e335–48. https://doi.org/10.1016/j.cllc.2018.01.002.

    Article  CAS  PubMed  Google Scholar 

  6. Jan YH, Lai TC, Yang CJ, Huang MS, Hsiao M. A co-expressed gene status of adenylate kinase 1/4 reveals prognostic gene signature associated with prognosis and sensitivity to EGFR targeted therapy in lung adenocarcinoma. Sci Rep. 2019;9(1):12329. https://doi.org/10.1038/s41598-019-48243-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Feng J, Ma Y, Chen Z, Hu J, Yang Q, Ding G. Mitochondrial pyruvate carrier 2 mediates mitochondrial dysfunction and apoptosis in high glucose-treated podocytes. Life Sci. 2019;237:116941. https://doi.org/10.1016/j.lfs.2019.116941.

    Article  CAS  PubMed  Google Scholar 

  8. Lu G, Lai Y, Wang T, Lin W, Lu J, Ma Y, et al. Mitochondrial fission regulator 2 (MTFR2) promotes growth, migration, invasion and tumour progression in breast cancer cells. Aging. 2019;11(22):10203–19. https://doi.org/10.18632/aging.102442.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Wu Y, Chen M, Jiang J. Mitochondrial dysfunction in neurodegenerative diseases and drug targets via apoptotic signaling. Mitochondrion. 2019;49:35–45. https://doi.org/10.1016/j.mito.2019.07.003.

    Article  CAS  PubMed  Google Scholar 

  10. Zhou Q, Li H, Li Y, Tan M, Fan S, Cao C, et al. Inhibiting neddylation modification alters mitochondrial morphology and reprograms energy metabolism in cancer cells. JCI Insight. 2019;4(4):e121582.

  11. Zeng Y, Guo Z, Hu Z, Liu M, Chen Y, Chen S, et al. FGD1 exhibits oncogenic properties in hepatocellular carcinoma through regulating cell morphology, autophagy and mitochondrial function. Biomed Pharmacother. 2020;125:110029.

    Article  CAS  PubMed  Google Scholar 

  12. Wang J, Xie Y, Bai X, Wang N, Yu H, Deng Z, et al. Targeting dual specificity protein kinase TTK attenuates tumorigenesis of glioblastoma. Oncotarget. 2018;9(3):3081–8. https://doi.org/10.18632/oncotarget.23152.

    Article  PubMed  Google Scholar 

  13. Yang M, Wang B, Gao J, Zhang Y, Xu W, Tao L. Spinosad induces programmed cell death involves mitochondrial dysfunction and cytochrome C release in Spodoptera frugiperda Sf9 cells. Chemosphere. 2017;169:155–61. https://doi.org/10.1016/j.chemosphere.2016.11.065.

    Article  CAS  PubMed  Google Scholar 

  14. Yang S, Liu Y, Guo Y, Liu R, Qi F, Li X, et al. Circadian gene clock participates in mitochondrial apoptosis pathways by regulating mitochondrial membrane potential, mitochondria out membrane permeablization and apoptosis factors in AML12 hepatocytes. Mol Cell Biochem. 2020;467(1–2):65–75. https://doi.org/10.1007/s11010-020-03701-1.

    Article  CAS  PubMed  Google Scholar 

  15. Pan B, Zhong W, Deng Z, Lai C, Chu J, Jiao G, et al. Inhibition of prostate cancer growth by solanine requires the suppression of cell cycle proteins and the activation of ROS/P38 signaling pathway. Cancer Medicine. 2016;5(11):3214–22. https://doi.org/10.1002/cam4.916.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Yue M, Li S, Yan G, Li C, Kang Z. Paeoniflorin inhibits cell growth and induces cell cycle arrest through inhibition of FoxM1 in colorectal cancer cells. Cell cycle (Georgetown, Tex). 2018;17(2):240–9.

    Article  CAS  Google Scholar 

  17. Li W, Zheng G, Xia J, Yang G, Sun J, Wang X, et al. Cell cycle-related and expression-elevated protein in tumor overexpression is associated with proliferation behaviors and poor prognosis in non-small-cell lung cancer. Cancer Sci. 2018;109(4):1012–23. https://doi.org/10.1111/cas.13524.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Nichols BA, Oswald NW, McMillan EA, McGlynn K, Yan J, Kim MS, et al. HORMAD1 is a negative prognostic Indicator in lung adenocarcinoma and specifies resistance to oxidative and genotoxic stress. Cancer Res. 2018;78(21):6196–208. https://doi.org/10.1158/0008-5472.CAN-18-1377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Sun H, Liu K, Huang J, Sun Q, Shao C, Luo J, et al. FAM111B, a direct target of p53, promotes the malignant process of lung adenocarcinoma. OncoTargets Ther. 2019;12:2829–42. https://doi.org/10.2147/OTT.S190934.

    Article  CAS  Google Scholar 

  20. Nihal M, Stutz N, Schmit T, Ahmad N, Wood GS. Polo-like kinase 1 (Plk1) is expressed by cutaneous T-cell lymphomas (CTCLs), and its downregulation promotes cell cycle arrest and apoptosis. Cell Cycle (Georgetown, Tex). 2011;10(8):1303–11.

    Article  CAS  Google Scholar 

  21. Raghavan P, Tumati V, Yu L, Chan N, Tomimatsu N, Burma S, et al. AZD5438, an inhibitor of Cdk1, 2, and 9, enhances the radiosensitivity of non-small cell lung carcinoma cells. Int J Radiat Oncol Biol Phys. 2012;84(4):e507–14. https://doi.org/10.1016/j.ijrobp.2012.05.035.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Ko E, Kim Y, Cho EY, Han J, Shim YM, Park J, et al. Synergistic effect of Bcl-2 and cyclin A2 on adverse recurrence-free survival in stage I non-small cell lung cancer. Ann Surg Oncol. 2013;20(3):1005–12. https://doi.org/10.1245/s10434-012-2727-2.

    Article  PubMed  Google Scholar 

  23. Kou F, Sun H, Wu L, Li B, Zhang B, Wang X, et al. TOP2A promotes lung adenocarcinoma Cells' malignant progression and predicts poor prognosis in lung adenocarcinoma. J Cancer. 2020;11(9):2496–508. https://doi.org/10.7150/jca.41415.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80. https://doi.org/10.1593/neo.07112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58. https://doi.org/10.1016/j.neo.2017.05.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Guo Q, Ke XX, Liu Z, Gao WL, Fang SX, Chen C, et al. Evaluation of the prognostic value of STEAP1 in lung adenocarcinoma and insights into its potential molecular pathways via Bioinformatic analysis. Front Genet. 2020;11:242. https://doi.org/10.3389/fgene.2020.00242.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Győrffy B, Surowiak P, Budczies J, Lánczky A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS One. 2013;8(12):e82241. https://doi.org/10.1371/journal.pone.0082241.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Li X, Sun L, Wang X, Wang N, Xu K, Jiang X, et al. A five immune-related lncRNA signature as a prognostic target for glioblastoma. Front Mol Biosci. 2021;8:632837. https://doi.org/10.3389/fmolb.2021.632837.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wu H, Zhang J. Decreased expression of TFAP2B in endometrial cancer predicts poor prognosis: a study based on TCGA data. Gynecol Oncol. 2018;149(3):592–7. https://doi.org/10.1016/j.ygyno.2018.03.057.

    Article  CAS  PubMed  Google Scholar 

  32. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49(D1):D605–12. https://doi.org/10.1093/nar/gkaa1074.

    Article  CAS  PubMed  Google Scholar 

  33. Ma H, He Z, Chen J, Zhang X, Song P. Identifying of biomarkers associated with gastric cancer based on 11 topological analysis methods of CytoHubba. Sci Rep. 2021;11(1):1331. https://doi.org/10.1038/s41598-020-79235-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

Thanks to Doctor Qing Luo for his nice help in the application of Databases.

Funding

This work was supported by Grants from the National Natural Science Foundation of China (No. 81960532).

Author information

Authors and Affiliations

  1. Department of Thoracic Surgery, Affiliated Hospital of Zunyi Medical University, Zunyi, 563000, Guizhou, China

    Cheng Chen, Yang Tang, Wen-Dong Qu, Xu Han, Jie-Bin Zuo, Qing-Yong Cai, Gang Xu, Yong-Xiang Song & Xi-Xian Ke

Authors
  1. Cheng Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yang Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Wen-Dong Qu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xu Han
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jie-Bin Zuo
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Qing-Yong Cai
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Gang Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Yong-Xiang Song
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Xi-Xian Ke
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

This research was finished by all the authors. SYX conceived the research topic, made the research plan and directed the implementation of the whole research. KXX drafted the manuscript and processed the data. CC, XG, CQY and TY assisted to collect and analyze the data and drafted the manuscript. QWD, HX and ZJB were mainly involved in the management of data collection. All authors approved the final manuscript.

Corresponding authors

Correspondence to Yong-Xiang Song or Xi-Xian Ke.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflicts of interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. Co-expressed gene of MTFR2.

Additional file 2: Table S2

. GO analysis of MTFR2 co-expressed genes.

Additional file 3: Figure S1

. PPI network displayed the relationship among MTFR2 coexpressed genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, C., Tang, Y., Qu, WD. et al. Evaluation of clinical value and potential mechanism of MTFR2 in lung adenocarcinoma via bioinformatics. BMC Cancer 21, 619 (2021). https://doi.org/10.1186/s12885-021-08378-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-021-08378-3

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us