Xiangwei Zhang
Wei Dong
Jishuai Zhang2
Duozhi Shi
Wei Ma- 1Department of General Thoracic, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
- 2Department of General Thoracic, Feicheng Hospital Affiliated to Shandong First Medical University, Feicheng, China
- 3Department of General Thoracic, Shenxian County People’s Hospital of Shandong Provincial Group, Liaocheng, China
- 4Lifehealthcare Clinical Laboratories, Hangzhou, China
Background: Lung cancer is the leading cause of cancer-related death worldwide, of which lung adenocarcinoma (LUAD) is one of the main histological subtypes. Mitochondria are vital for maintaining the physiological function, and their dysfunction has been found to be correlated with tumorigenesis and disease progression. Although, some mitochondrial-related genes have been found to correlate with the clinical outcomes of multiple tumors solely. The integrated relationship between nuclear mitochondrial genes (NMGs) and the prognosis of LUAD remains unclear.
Methods: The list of NMGs, gene expression data, and related clinical information of LUAD were downloaded from public databases. Bioinformatics methods were used and obtained 18 prognostic related NMGs to construct a risk signature.
Results: There were 18 NMGs (NDUFS2, ATP8A2, SCO1, COX14, COA6, RRM2B, TFAM, DARS2, GARS, YARS2, EFG1, GFM1, MRPL3, MRPL44, ISCU, CABC1, HSPD1, and ETHE1) identified by LASSO regression analysis. The mRNA expression of these 18 genes was positively correlated with their relative linear copy number alteration (CNA). Meanwhile, the established risk signature could effectively distinguish high- and low-risk patients, and its predictive capacity was validated in three independent gene expression omnibus (GEO) cohorts. Notably, a significantly lower prevalence of actionable EGFR alterations was presented in patients with high-risk NMGs signature but accompanied with a more inflame immune tumor microenvironment. Additionally, multicomponent Cox regression analysis showed that the model was stable when risk score, tumor stage, and lymph node stage were considered, and the 1-, 3-, and 5-year AUC were 0.74, 0.75, and 0.70, respectively.
Conclusion: Together, this study established a signature based on NMGs that is a prognostic biomarker for LUAD patients and has the potential to be widely applied in future clinical settings.
Introduction
Mitochondria are complex organelles of bioenergetic, biosynthetic, and signaling that are correlated to several diseases, including cardiovascular diseases, neurological disorders, and metabolism disorders (Vyas et al., 2016; Genovese et al., 2020). In addition, recently, mitochondria are demonstrated to participate in multiple aspects of tumor formation and progression (Hsu et al., 2016; Bonora et al., 2021). The mechanisms of mitochondria involved in tumorigenesis have been extensively studied (Nunes et al., 2015), and some specific nuclear mitochondrial genes (NMGs) were regarded as the potential targets for the development of the next generation of cancer therapeutics (Su et al., 2016; Zong et al., 2016).
mtDNA is the own genome of mammalian mitochondria, 16.5- kb double-stranded circular DNA, which encodes a total of 13 proteins that are all core components of oxidative phosphorylation (OXPHOS). However, all remaining mitochondrial functions proteins (∼1,300) are encoded in the nuclear DNA (nDNA). Following their translation, they are imported into the organelle through the import machinery instead (Anderson et al., 1981; DiMauro and Schon, 2003; Schon and Przedborski, 2011). Thus, the mitochondrial function requires cooperation between the nuclear and mitochondrial genomes. Mitochondrial genome instability and mitochondrial dysfunction are novel markers for cancer (Wang et al., 2015; Idaghdour and Hodgkinson, 2017). Mitochondrial dysfunction, caused by smoking or other factors, plays a vital role in the carcinogenesis of multiple cancers, especially lung cancer (Cloonan et al., 2020). The dysfunction of mitochondria may shift the energy production from mitochondria-mediated oxidative to glycolysis, which in turn increases the cell metabolism and activity (Chen et al., 2021). Meanwhile, mitochondrial dysfunction causes ROS outbreak, which mediates epithelial-mesenchymal transition and cell invasion in lung cancer via coordinate intracellular signal transduction (Nunes et al., 2015; He et al., 2016; Han et al., 2018). mtDNA variations have been frequently suggested to serve as initiators for a variety of tumors, such as prostate cancer (Hopkins et al., 2017), LUAD (Yuan et al., 2015), breast cancer (Weerts et al., 2018), pancreatic cancer (Lam et al., 2012), gynecological malignancies (Guerra et al., 2014), and acute myeloid leukemia (Reznik et al., 2016; Tyagi et al., 2018). Contrary to the increased study on the association of mtDNA alterations and cancer, fewer studies have shed light on the role of NMGs on cancer. The germline or somatic alteration in nuclear mitochondrial complex II genes encoding SDH subunits (SDHB, C, and D) were found associated with pheochromocytomas and paragangliomas (Gimm et al., 2000; van Nederveen et al., 2007; Bardella et al., 2011). Similarly, one study found that the NMGs NDUFS1 and NDUFS8 (encoding subunits of mitochondrial complex 1) had significant prognostic power in the patients with non-small cell lung cancer (NSCLC) by analyzed immunohistochemical staining and RNA expression data (Su et al., 2016). Alterations in the genes mentioned above have also been detected in patients with mitochondrial disease, such as mitochondrial complex I deficiency and mitochondrial respiratory chain deficiency (Bourgeron et al., 1995; Loeffen et al., 1998; Björkman et al., 2015). Therefore, we examined the hypothesis that the NMGs may have significant predictive value in the identification of high-risk cancer patients with poor overall clinical survival.
Lung cancer is the leading cause of cancer-related mortality worldwide causing over 1.7 million deaths annually (Bray et al., 2018), of which LUAD are the major subtypes. In the present study, we investigate the NMGs associated with LUAD survival by bioinformatics analysis, following univariate Cox regression analyses to develop an NMGs signature to provide new clues for improving the diagnosis and treatment of patients with LUAD.
Materials and Methods
Data Sources
The data of gene expression and corresponding clinical of LUAD were retrieved from cbioportal (TCGA, Firehose Legacy) and Gene Expression Omnibus (GEO) (GSE13213, GSE42127, and GSE72094) database. TCGA-LUAD dataset was chosen as the training dataset, and the other three datasets were selected as external validation. The last of NMGs was downloaded from MITOMAP (Last update: January 2018)1 and was detailed in Supplementary Table 1.
Construction and Evaluation of a Survival-Related Nuclear Mitochondrial Genes Signature
First, univariate Cox regression analysis was used to identify prognosis-related NMGs (P < 0.05). Then, the least absolute shrinkage and selection operator (LASSO)-penalized Cox regression analysis was performed to build a prognostic model with the “glmnet” R package. The risk score calculation was based on the expression level of the normalized gene and regression coefficient of the corresponding gene, which was as follows: risk score = sum (expression level of each gene × coefficient of corresponding genes).
The patients in the training and validation cohorts were grouped into high- or low-risk groups based on the median risk score. Kaplan-Meier (KM) analysis was performed to compare the survival differences between the high- and low-risk groups using the R package. In addition, the area under the curve (AUC) was calculated to measure the prognostic capability of the NMGs signature, and the nomogram was drawn by the R package.
Multivariate Cox regression analyses were performed to determine the prognostic values for the signature and some clinical features. R package “rms” was used to draw nomograms and calibration curves.
Gene Set Enrichment Analysis
To explore biological processes in the high-risk and low-risk groups, GSEA was performed by ClusterProfiler package in R studio. The HALLMARK gene sets and Kyoto Encyclopedia of Genes and Genomes (KEGG) gene sets from the Molecular Signatures Database (MSigDB)2 were used.
Genomic Mutations and Tumor Microenvironment Analysis
To compare the differences in genomic mutations between high and low risk groups, mutation profiles were analyzed and visualized using the “maftools” R. Meanwhile, transcriptome profiles from the TCGA cohort were used to identify the immune cell fractions of 22 distinct leukocyte subsets by the CIBERSORT tool.
Statistical Analyses
Fisher’s exact test was applied to compare the difference in proportions between the high- and low-risk groups. R software was used to perform statistical analysis and all statistical results with p-values < 0.05 were considered statistically significant.
Results
Baseline Characteristics of Lung Adenocarcinoma Patients Included in This Study
In this study, TCGA-LUAD, which included 502 samples, was applied as the discovery dataset. In addition, GSE13213, GSE42127, and GSE72094 databases, which contained 117, 133, and 420 LUAD patients, were used as the validation cohort. The details of baseline characteristics are shown in Table 1.
Table 1. Baseline characteristics of the patients in this study.
Construction of a Nuclear Mitochondrial Genes Signature
The patients from TCGA were used as a training cohort to identify an NMGs signature. A univariate Cox regression analysis was conducted, and 26 survival-related NMGs (Figure 1A) from a total of 146 NMGs were identified (P < 0.05) (Supplementary Table 1). Among the 26 genes, the hazard ratio (HR) of seven genes (ATP8A2, CABC1, TK2, ANT1, RRM2B, COX14, and ISCU) was less than 1, while the HR of the remaining 19 genes was more than 1, which indicated that the 7 genes were associated with a poor prognosis in LUAD patients.
Figure 1. Construction of a survival-related signature based on nuclear mitochondrial genes (NMGs). (A) The forest plot of the associations between the expression levels of 26 genes and overall survival in LUAD patients. Hazard ratios (HR), P-value, and corresponding 95% confidence intervals were calculated by univariate Cox regression analysis. (B) The result of LASSO regression analysis. The vertical dotted lines were drawn at the optimal values by using the maximum and lambda.min criteria. (C) Kaplan-Meier OS curves for patients with high- and low-risk NMGs signature in the TCGA cohort. (D) The ROC curves for 1-, 3-, and 5-year survival of established NMGs signature. OS, overall survival; LASSO, the least absolute shrinkage and selection operator; ROC, receiver operating characteristic. *P < 0.05.
Then we established a risk score model with these identified 26 survival-related NMGs using LASSO-based Cox regression analysis and 18 NMGs were selected (Figure 1B). The risk score was calculated using the following formula:
Using the median risk score, we classified the patients from the TCGA cohort into high-risk and low-risk groups. The Kaplan-Meier curve indicated that the patients in the low-risk group exhibited a longer survival time than those in the high-risk group [median overall survival (OS], 60.12 vs. 38.47 months, P < 0.0001, Figure 1C). Furthermore, the 1-, 3-, and 5-year OS predictions for LUAD patients and all the calibrated curves were well-fitted with AUCs of 0.72, 0.68, and 0.64, respectively (Figure 1D).
We analyzed the expression of the 18 NMGs in LUAD samples compared to normal samples using GEPIA to compile Genome Tissue Expression (GTEx) and TCGA datasets. The 18 genes were found to be expressed in both tumor and normal tissues, of which COA6, DARS2, and MRPL3 were significantly more expressed in tumor tissues (P < 0.05, Figure 2).
Figure 2. Expression analyses of the 18 genes (NDUFS2, ATP8A2, SCO1, COX14, COA6, RRM2B, TFAM, DARS2, GARS, YARS2, EFG1, GFM1, MRPL3, MRPL44, ISCU, CABC1, HSPD1, ETHE1) in LUAD samples and normal samples. *P < 0.05.
18 Nuclear Mitochondrial Genes Signature Validation in Three Gene Expression Omnibus Cohorts
We assessed the performance of established NMG’s signature using the three independent validation cohorts (GSE13213, GSE42127, and GSE72094) from the GEO. Based on the median value of risk score calculated as described above, the patients in the GSE13213, GSE42127, and GSE72094 were assigned to the high-risk and low-risk groups, respectively. The results were satisfactory that the patients in the low-risk group in all three validation cohorts had better survival than those in the high-risk group (all P < 0.05, Figures 3A,C,E). The receiver operating characteristic (ROC) curve for 1-, 3-, and 5-year OS predictions suggested that the model possessed predictive accuracy with AUCs above 0.600 (Figures 3B,D,F).
Figure 3. Validation of the established NMGs signature. Kaplan-Meier OS curves for patients assigned to high- and low-risk groups based on the risk score in the GSE13213 (A), GSE42127 (C), and GSE72094 (E) cohort. The ROC curves of 1-, 3-, and 5-year OS in the GSE13213 (B), GSE42127 (D), and GSE72094 (F) cohort. NMGs: nuclear mitochondrial genes; OS: overall survival; ROC: receiver operating characteristic.
Multivariate COX Regression Analyses Regarding Overall Survival in the TCGA-LUAD Cohort
As depicted in Supplementary Figure 1, the risk model was associated with Fev1 Fvc ratio of prebroncholiator, neoplasm disease stage, tumor (T) stage, and lymph node (N) stage (all P < 0.05). Meanwhile, it was unrelated with Fev1 Fvc ratio of postbroncholiator, age, long-distant metastasis (M) stage, adjuvant postoperative targeted therapy, and smoking history. Therefore, the nomogram was constructed based on the variables of risk score, Fev1 Fvc ratio of prebroncholiator, neoplasm disease stage, T stage, and N stage to predict the OS of patients with LUAD (Figure 4A). The concordant index (C-index) of the nomogram model to predict OS was 0.65 and the calibration curves of the 1-, 3-, and 5-year survival prediction were all close to the ideal performance of Figure 4B, exhibiting good predictive accuracy. Subsequently, these four variables (Fev1 Fvc ratio of prebroncholiator, neoplasm disease stage, T stage, and N stage) and risk score were selected for multivariate Cox regression analysis. Risk score was independently associated with a worse survival (HR = 3.56, 95% CI = 1.70–7.46, P < 0.001); on the contrary, low T stage (HR = 0.35, 95% CI = 0.15–0.84, P = 0.018) and low N stage (HR = 0.43, 95% CI = 0.21–0.86, P = 0.017) were significantly associated with better OS (Figure 4C). Then, the AUC for the 1-, 3-, and 5-year survival of the constructed nomogram were 0.74, 0.75, and 0.70, respectively (Figure 4D), indicating that the model was stable and robust for predicting the survival of patients with LUAD.
Figure 4. Development of prognostic NMGs nomogram for LUAD patients. (A) The nomogram for predicting 1-, 3-, and 5-year OS probabilities of LUAD patients. (B) Calibration plot of the established nomogram for predicting probabilities of 1-, 3-, and 5-year OS. (C) Multivariate Cox regression analysis of risk factors (including risk score, Fev1 Fvc ratio of prebroncholiator, neoplasm disease stage, T stage, and N stage). (D) The ROC curves of the nomogram for predicting 1-, 3-, and 5-year OS. OS, overall survival; ROC, receiver operating characteristic; T stage, tumor stage; N stage, lymph node stage.
Biological Processes in the High-Risk and Low-Risk Subgroups
HALLMARK enriched results revealed that the E2F targets, MYC targets, DNA repair, and G2M checkpoint signaling pathways were highly enriched in the high-risk subgroups (P < 0.05, Figure 5A). In other words, overexpression of the genes in the pathways mentioned above was positively correlated with the risk score. Meanwhile, we conducted the KEGG pathway analysis and found that NMGs were involved in the cell cycle, proteasome, and DNA replication signaling pathways (Figure 5B). Interestingly, both HALLMARK and KEGG analyses indicated that NMGs were enriched in metabolism-related signaling pathways such as glycolysis and oxidative phosphorylation signaling pathway in the high-risk subgroups.
Figure 5. HALLMARK (A) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (B) enrichment analysis in the high-risk and low-risk subgroups. The enriched items were selected with a corrected q-value 0.01; the length and color of the bar represent the absolute value of NES and the q-values, respectively. NES: normalized enrichment score.
Genomic Mutations and the Immune Microenvironment in the High-Risk and Low-Risk Groups
The top 20 most frequently mutated genes in the high-risk and low-risk groups are shown in Figures 6A,B. The mutation count in the high-risk group was significantly higher than that in the low-risk group (Figure 6C). In addition, by comparison the prevalence of actionable alterations in LUAD between two groups, we found that there were more patients in the low-risk group having actionable EGFR alterations (Figure 6D).
Figure 6. Genomic mutations and tumor microenvironment analysis in the high-risk and low-risk subgroups. The top 20 most frequently mutated genes in high-risk (A) and low-risk (B) subgroups. (C) Comparison of mutation count between the two risk subgroups. (D) Comparison of the prevalence of actionable genes in LUAD between two risk subgroups. (E) Comparison of immune cell infiltration between two risk subgroups.
Meanwhile, analysis of the tumor-infiltrated cells of TCGA-LUAD indicated the difference in immune cell infiltration between the high-risk and low-risk subgroups. The low-risk group had a significantly higher abundance of memory B cells, resting CD4 memory T cell, regulatory T cell (Tregs), monocyte, resting myeloid dendritic cell, and activated mast cell. On the other hand, the level of CD8 T cell, activated CD4 memory T cell, helper follicular T cell, activated NK cell, M0 macrophage, M1 macrophage, and resting mast cell in the high-risk group were higher than that of the low-risk group (Figure 6E).
Genetic Alteration in Nuclear Mitochondrial Genes
In the TCGA-LUAD cohort, 10.87% of patients have genetic alterations among 18 NMGs, of which missense mutation accounts for the most prevalent type. The most prevalent altered gene is ATP8A2 (6%), whereas other genes were rarely altered (Figure 7A). In addition, we found the mRNA expression of 18 genes was positively correlated with their relative linear copy number alteration (CNA) value, especially EFG1, GFM1, NDUFS2, SCO1, and YARS2 (Figure 7B).
Figure 7. Genetic alteration and copy number alterations (CNAs) of the 18 genes in LUAD samples. (A) Genetic alteration percentage of 18 genes. (B) Correlation of each 18 genes mRNA expression level and linearized copy-number alteration.
The Therapeutic Benefit of the Risk Score Value
Pathway enrichment analysis showed that overexpression of genes in the DNA repair pathway was positively correlated with risk score, which suggests that a higher risk score indicated an increased ability of DNA repair. Thus, we speculate that the patients in the high-risk group exhibited a weaker response to radiation therapy compared with the low-risk group. To confirm our hypothesis, the patients who received radiation therapy (N = 61) from the TCGA cohort were separated into a high-risk and low-risk group based on the median risk score, and the patients with higher risk scores exhibited significantly worse prognosis (p = 0.038) (Supplementary Figure 2), which indicated NMGs signature might serve as a potential indicator of patients with LUAD response to radiotherapy.
Discussion
The pivotal roles of mitochondria in maintaining normal cell function as well as in cancer development are widely acknowledged. Dysfunction in some specific members of mitochondrial genes, especially FH, SDH family, and IDH1/2, have been demonstrated as the milestone event for the development of multiple cancers like gliomas (Hartmann et al., 2009), glioblastomas (Balss et al., 2008), acute myeloid leukemia (Paschka et al., 2010), pheochromocytoma (Astuti et al., 2001), paraganglioma (Niemann and Müller, 2000; Astuti et al., 2001), and papillary renal cell cancer (Tomlinson et al., 2002; Luo et al., 2020). However, the integrated role of NMGs in the prediction for the prognosis of lung cancer patients remains undefined.
In this study, we first identified 26 NMGs whose expression were correlated with LUAD patients’ survival. Among those genes, only seven genes (ATP8A2, CABC1, TK2, ANT1, RRM2B, COX14, and ISCU) were associated with a favorable outcome in LUAD patients. Although these genes were previously reported in some cancer types, their specific role in lung cancer is unclarified except for ATP8A2 and RRM2B. ATP8A2, encoded ATPase Phospholipid Transporting 8A2, belongs to the P4-ATPase family that actively flips phosphatidylserine and phosphatidylethanolamine from the exoplasmic to the cytoplasmic leaflet of cell membranes to generate and maintain phospholipid asymmetry (Choi et al., 2019). As this gene involves neurite elongation and neuron survival, various genomic variants of ATP8A2 are identified in patients with neurological disorders (Guissart et al., 2020). Though the relationship between ATP8A2 and lung cancer is poorly investigated, a recent study found the alteration in its expression was associated with the prognosis of LUAD (Wang et al., 2020). ATP8A2 was hypermethylated in tumor tissues and, therefore, poorly expressed in the tumor tissues compared with the normal tissues (Wang et al., 2020). The same as the previously published result by Wang et al. (2020), we found ATP8A2 was associated with a better prognosis in LUAD. RRM2B, which has been widely recognized for its important role in maintaining genomic stability, was widely overexpressed in many types of tumors. However, a previous study suggested the distinct roles of RRM2B in tumor progression: on one hand, its amplification was related to worse outcomes in breast cancers (Chae et al., 2016; Iqbal et al., 2021); on the other hand, in vitro analysis found its overexpression could inhibit the proliferation of lung cancer cells by the regulation work of FOXO3 (Cho et al., 2014). The latter result was also supported by our results as RRM2B was significantly associated with better survival in LUAD patients. Meanwhile, 18 NMGs significantly correlated with unfavorable prognosis in LUAD patients were identified. Among them, SCO1, HSPD1, IARS2, MRPL44, TFAM, NDUFS2, MRPL13, OPA1, and YARS have been previously identified as oncogenes in lung cancer or other types of cancers, serving as unfavorable prognosis biomarkers individually (Telang et al., 2012; Lee et al., 2017; Sotgia and Lisanti, 2017; Di et al., 2019; Dunham-Snary et al., 2019; Liu et al., 2019; Witherspoon et al., 2019; Li et al., 2020; Zhang et al., 2020; Cai et al., 2021). However, it is noteworthy that nearly all these findings were established and assessed in vitro, and all these NMGs were analyzed separately instead of integrating. Our novel identified unfavorable prognosis predictors, including GFM1, KARS, DARS2, GARS, DLP1, MRPL12, COA6, SCO1, and C10ORF have not been functionally analyzed previously. We are the first study to confirm their association with the progression of LUAD patients. Meanwhile, previous research found COA6, an OXPHOS complexes assembly factor, negatively correlated with vimentin levels, which was suggested as a favorable predictor for cancer prognosis (Cruz-Bermúdez et al., 2019). On the contrary, our findings suggested its association with unfavorable survival in LUAD patients. Furthermore, MRPL44 and MRPL12 were both members of mitochondrial ribosomal proteins, whose abnormal expression was found associated with the tumorigenesis and development of lung cancer (Huang et al., 2020). Glycolysis is a dominant metabolism for cancer cells to produce energy, even in the presence of oxygen. Research has demonstrated that a higher glycolytic rate in tumor cells was shown to promote resistance to chemotherapeutics (Ganapathy-Kanniappan and Geschwind, 2013). On the contrary, it has been revealed that mitochondrial OXPHOS is also utilized by cancer cells. Meanwhile, some evidence suggests that OXPHOS contributes to cancer metastasis (Viale et al., 2015; Yu et al., 2017; Rao et al., 2019). In the present study, glycolysis and oxidative phosphorylation signaling pathways were enriched in the high-risk subgroup in HALLMARK and KEGG enrichment analyses. These signaling pathways may play important roles in tumorigenesis and the progression of LUAD. Nowadays, the treatment landscape of NSCLC includes target therapy, immunotherapy, and combination therapy, which is mainly decided by the patients’ driver genes status, especially EGFR and ALK (Grant et al., 2021). Interestingly, our findings showed that LUAD patients with high or low NMGs signature may have distinct treatment choices as the differences in the tumor microenvironment and driver alterations. LUAD patients with low NMGs signature have not only improved outcomes but also significantly higher prevalence of EGFR alterations, which may contribute to the response to EGFR tyrosine-kinase inhibitors. On the contrary, patients with high NMGs have more inflamed immune features, including higher mutation counts, high presence of immune-prone TILs (such as CD8 positive T cells, activated CD4 memory T cell and NK cell), and lower presence of immunosuppressive TILs (such as Treg cells). Previous studies have found the tumor-infiltrated CD8 positive T cells prior to immune checkpoint inhibitors could serve as an indicator for the response of anti-PD-1/L1 therapy (Tumeh et al., 2014; McDermott et al., 2018). Based on these features, it is indicated that LUAD patients with high NMGs signature may benefit more from immune checkpoint inhibitors (Lucibello et al., 2021). In concordance with a more inflamed tumor microenvironment, tumor mutation counts (also known as tumor mutation burden, TMB) level was significantly higher in the high NMGs signature than in the low group, which has already been granted as an effective biomarker for selecting patients who may respond to pembrolizumab by FDA (Marcus et al., 2021).
There were some limitations in our study. Although with ample validation, the signature was established and validated based on the public datasets. Further validation in a local cohort using prospective clinical samples and data is merited. Though we found the association between NMGs signature and tumor microenvironment or driver gene alterations (especially EGFR), which may affect the decision-making for clinical management for LUAD patients, the validation of the indication that the difference in the efficacy of immune checkpoint inhibitor and targeted therapy on patients with high- or low-risk signature is worthy to explore in further study, as well as the underlying regulatory mechanism in vivo or in vitro. Based on the above limitations, we have designed a local LUAD cohort from our hospital and performed testing of the involved NMGs based on real-time PCR to further verify the results of this study.
In conclusion, we conducted an integrative analysis of the role of NMGs in the prognosis of LUAD and successfully established a robust and stable prognosis predicting signature with sufficient validation. Except for prognostic function, this established signature could also distinct LUAD patients who may have more clinical benefits from radiotherapy, targeted therapy, and immune checkpoint inhibitors, which may give clinicians more decision-making help to manage patients.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories is cBioPortal and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
WM, XZ, and WD: conception and design, data analysis and interpretation (statistical analysis), and manuscript writing and revision. JZ, WL, DS, and JY: data acquisition and manuscript writing. All authors read and approved the final manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (81802282), Natural Science Foundation of Shandong Province (ZR2020MH233), Jinan Science and Technology Development Program (201907112), and Clinical Medicine Technology and Innovations Plan of Jinan City (201907066).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.740487/full#supplementary-material
Supplementary Figure 1 | Comparison of clinical characteristics between the high-risk and low-risk groups.
Supplementary Figure 2 | The therapeutic benefit of the risk score value. Kaplan-Meier OS curves for patients who received radiation therapy assigned to high- and low-risk groups in the TCGA cohort.
Footnotes
References
Anderson, S., Bankier, A. T., Barrell, B. G., de Bruijn, M. H., Coulson, A. R., Drouin, J., et al. (1981). Sequence and organization of the human mitochondrial genome. Nature 290, 457–465. doi: 10.1038/290457a0
Astuti, D., Latif, F., Dallol, A., Dahia, P. L., Douglas, F., George, E., et al. (2001). Gene mutations in the succinate dehydrogenase subunit SDHB cause susceptibility to familial pheochromocytoma and to familial paraganglioma. Am. J. Hum. Genet. 69, 49–54. doi: 10.1086/321282
Balss, J., Meyer, J., Mueller, W., Korshunov, A., Hartmann, C., and von Deimling, A. (2008). Analysis of the IDH1 codon 132 mutation in brain tumors. Acta Neuropathol. 116, 597–602. doi: 10.1007/s00401-008-0455-2
Bardella, C., Pollard, P. J., and Tomlinson, I. (2011). SDH mutations in cancer. Biochim. Biophys. Acta 1807, 1432–1443. doi: 10.1016/j.bbabio.2011.07.003
Björkman, K., Sofou, K., Darin, N., Holme, E., Kollberg, G., Asin-Cayuela, J., et al. (2015). Broad phenotypic variability in patients with complex I deficiency due to mutations in NDUFS1 and NDUFV1. Mitochondrion 21, 33–40. doi: 10.1016/j.mito.2015.01.003
Bonora, M., Missiroli, S., Perrone, M., Fiorica, F., Pinton, P., and Giorgi, C. (2021). Mitochondrial control of genomic instability in cancer. Cancers (Basel) 13:1914. doi: 10.3390/cancers13081914
Bourgeron, T., Rustin, P., Chretien, D., Birch-Machin, M., Bourgeois, M., Viegas-Péquignot, E., et al. (1995). Mutation of a nuclear succinate dehydrogenase gene results in mitochondrial respiratory chain deficiency. Nat. Genet. 11, 144–149. doi: 10.1038/ng1095-144
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492
Cai, M., Li, H., Chen, R., and Zhou, X. (2021). MRPL13 promotes tumor cell proliferation, migration and EMT process in breast cancer through the PI3K-AKT-mTOR PATHWAY. Cancer Manag. Res. 13, 2009–2024. doi: 10.2147/cmar.S296038
Chae, Y., Anker, J., Carneiro, B., Chandra, S., Kaplan, J., Kalyan, A., et al. (2016). Genomic landscape of DNA repair genes in cancer. Oncotarget 7, 23312–23321. doi: 10.18632/oncotarget.8196
Chen, X., Hao, B., Li, D., Reiter, R. J., Bai, Y., Abay, B., et al. (2021). Melatonin inhibits lung cancer development by reversing the Warburg effect via stimulating the SIRT3/PDH axis. J. Pineal Res. 71:e12755. doi: 10.1111/jpi.12755
Cho, E., Kuo, M., Liu, X., Yang, L., Hsieh, Y., Wang, J., et al. (2014). Tumor suppressor FOXO3 regulates ribonucleotide reductase subunit RRM2B and impacts on survival of cancer patients. Oncotarget 5, 4834–4844. doi: 10.18632/oncotarget.2044
Choi, H., Andersen, J. P., and Molday, R. S. (2019). Expression and functional characterization of missense mutations in ATP8A2 linked to severe neurological disorders. Hum. Mutat. 40, 2353–2364. doi: 10.1002/humu.23889
Cloonan, S. M., Kim, K., Esteves, P., Trian, T., and Barnes, P. J. (2020). Mitochondrial dysfunction in lung ageing and disease. Eur. Respir. Rev. 29:200165. doi: 10.1183/16000617.0165-2020
Cruz-Bermúdez, A., Laza-Briviesca, R., Vicente-Blanco, R., García-Grande, A., Coronado, M., Laine-Menéndez, S., et al. (2019). Cancer-associated fibroblasts modify lung cancer metabolism involving ROS and TGF-β signaling. Free Rad. Biol. Med. 130, 163–173. doi: 10.1016/j.freeradbiomed.2018.10.450
Di, X., Jin, X., Ma, H., Wang, R., Cong, S., Tian, C., et al. (2019). The oncogene IARS2 promotes non-small cell lung cancer tumorigenesis by activating the AKT/MTOR pathway. Front. Oncol. 9:393. doi: 10.3389/fonc.2019.00393
DiMauro, S., and Schon, E. A. (2003). Mitochondrial respiratory-chain diseases. N. Engl. J. Med. 348, 2656–2668. doi: 10.1056/NEJMra022567
Dunham-Snary, K., Wu, D., Potus, F., Sykes, E., Mewburn, J., Charles, R., et al. (2019). Ndufs2, a core subunit of mitochondrial complex I, is essential for acute oxygen-sensing and hypoxic pulmonary vasoconstriction. Circ. Res. 124, 1727–1746. doi: 10.1161/circresaha.118.314284
Ganapathy-Kanniappan, S., and Geschwind, J. F. (2013). Tumor glycolysis as a target for cancer therapy: progress and prospects. Mol. Cancer 12:152. doi: 10.1186/1476-4598-12-152
Genovese, I., Vezzani, B., Danese, A., Modesti, L., Vitto, V. A. M., Corazzi, V., et al. (2020). Mitochondria as the decision makers for cancer cell fate: from signaling pathways to therapeutic strategies. Cell Calcium 92:102308. doi: 10.1016/j.ceca.2020.102308
Gimm, O., Armanios, M., Dziema, H., Neumann, H. P., and Eng, C. (2000). Somatic and occult germ-line mutations in SDHD, a mitochondrial complex II gene, in nonfamilial pheochromocytoma. Cancer Res. 60, 6822–6825.
Grant, M., Herbst, R., and Goldberg, S. (2021). Selecting the optimal immunotherapy regimen in driver-negative metastatic NSCLC. Nat. Rev. Clin. Oncol. 18, 625–644. doi: 10.1038/s41571-021-00520-1
Guerra, F., Girolimetti, G., Perrone, A. M., Procaccini, M., Kurelac, I., Ceccarelli, C., et al. (2014). Mitochondrial DNA genotyping efficiently reveals clonality of synchronous endometrial and ovarian cancers. Mod. Pathol. 27, 1412–1420. doi: 10.1038/modpathol.2014.39
Guissart, C., Harrison, A. N., Benkirane, M., Oncel, I., Arslan, E. A., Chassevent, A. K., et al. (2020). ATP8A2-related disorders as recessive cerebellar ataxia. J. Neurol. 267, 203–213. doi: 10.1007/s00415-019-09579-4
Han, S. Y., Jeong, Y. J., Choi, Y., Hwang, S. K., Bae, Y. S., and Chang, Y. C. (2018). Mitochondrial dysfunction induces the invasive phenotype, and cell migration and invasion, through the induction of AKT and AMPK pathways in lung cancer cells. Int. J. Mol. Med. 42, 1644–1652. doi: 10.3892/ijmm.2018.3733
Hartmann, C., Meyer, J., Balss, J., Capper, D., Mueller, W., Christians, A., et al. (2009). Type and frequency of IDH1 and IDH2 mutations are related to astrocytic and oligodendroglial differentiation and age: a study of 1,010 diffuse gliomas. Acta Neuropathol. 118, 469–474. doi: 10.1007/s00401-009-0561-9
He, K., Guo, X., Liu, Y., Li, J., Hu, Y., Wang, D., et al. (2016). TUFM downregulation induces epithelial—mesenchymal transition and invasion in lung cancer cells via a mechanism involving AMPK-GSK3β signaling. Cell Mol. Life Sci. 73, 2105–2121. doi: 10.1007/s00018-015-2122-9
Hopkins, J. F., Sabelnykova, V. Y., Weischenfeldt, J., Simon, R., Aguiar, J. A., Alkallas, R., et al. (2017). Mitochondrial mutations drive prostate cancer aggression. Nat. Commun. 8:656. doi: 10.1038/s41467-017-00377-y
Hsu, C. C., Tseng, L. M., and Lee, H. C. (2016). Role of mitochondrial dysfunction in cancer progression. Exp. Biol. Med. (Maywood) 241, 1281–1295. doi: 10.1177/1535370216641787
Huang, G., Li, H., and Zhang, H. (2020). Abnormal expression of mitochondrial ribosomal proteins and their encoding genes with cell apoptosis and diseases. Int. J. Mol. Sci. 21:8879. doi: 10.3390/ijms21228879
Idaghdour, Y., and Hodgkinson, A. (2017). Integrated genomic analysis of mitochondrial RNA processing in human cancers. Genome Med. 9:36. doi: 10.1186/s13073-017-0426-0
Iqbal, W., Demidova, E., Serrao, S., ValizadehAslani, T., Rosen, G., and Arora, S. (2021). RRM2B is frequently amplified across multiple tumor types: implications for DNA repair, cellular survival, and cancer therapy. Front. Genet. 12:628758. doi: 10.3389/fgene.2021.628758
Lam, E. T., Bracci, P. M., Holly, E. A., Chu, C., Poon, A., Wan, E., et al. (2012). Mitochondrial DNA sequence variation and risk of pancreatic cancer. Cancer Res. 72, 686–695. doi: 10.1158/0008-5472.can-11-1682
Lee, W., Na, H., Lee, S., Lim, W., Kim, N., Lee, J., et al. (2017). Transcriptomic analysis of mitochondrial TFAM depletion changing cell morphology and proliferation. Sci. Rep. 7:17841. doi: 10.1038/s41598-017-18064-9
Li, M., Wang, L., Wang, Y., Zhang, S., Zhou, G., Lieshout, R., et al. (2020). Mitochondrial fusion Via OPA1 and MFN1 supports liver tumor cell metabolism and growth. Cells 9:121. doi: 10.3390/cells9010121
Liu, L., Qi, L., Knifley, T., Piecoro, D., Rychahou, P., Liu, J., et al. (2019). S100A4 alters metabolism and promotes invasion of lung cancer cells by up-regulating mitochondrial complex I protein NDUFS2. J. Biol. Chem. 294, 7516–7527. doi: 10.1074/jbc.RA118.004365
Loeffen, J., Smeitink, J., Triepels, R., Smeets, R., Schuelke, M., Sengers, R., et al. (1998). The first nuclear-encoded complex I mutation in a patient with Leigh syndrome. Am. J. Hum. Genet. 63, 1598–1608. doi: 10.1086/302154
Lucibello, G., Mograbi, B., Milano, G., Hofman, P., and Brest, P. (2021). PD-L1 regulation revisited: impact on immunotherapeutic strategies. Trends Mol. Med. 27, 868–881. doi: 10.1016/j.molmed.2021.06.005
Luo, Y., Ma, J., and Lu, W. (2020). The significance of mitochondrial dysfunction in cancer. Int. J. Mol. Sci. 21:5598. doi: 10.3390/ijms21165598
Marcus, L., Fashoyin-Aje, L., Donoghue, M., Yuan, M., Rodriguez, L., Gallagher, P., et al. (2021). FDA approval summary: pembrolizumab for the treatment of tumor mutational burden-high solid tumors. Clin. Cancer Res. 27, 4685–4689. doi: 10.1158/1078-0432.Ccr-21-0327
McDermott, D., Huseni, M., Atkins, M., Motzer, R., Rini, B., Escudier, B., et al. (2018). Clinical activity and molecular correlates of response to atezolizumab alone or in combination with bevacizumab versus sunitinib in renal cell carcinoma. Nat. Med. 24, 749–757. doi: 10.1038/s41591-018-0053-3
Niemann, S., and Müller, U. (2000). Mutations in SDHC cause autosomal dominant paraganglioma, type 3. Nat. Genet. 26, 268–270. doi: 10.1038/81551
Nunes, J. B., Peixoto, J., Soares, P., Maximo, V., Carvalho, S., Pinho, S. S., et al. (2015). OXPHOS dysfunction regulates integrin-β1 modifications and enhances cell motility and migration. Hum. Mol. Genet. 24, 1977–1990. doi: 10.1093/hmg/ddu612
Paschka, P., Schlenk, R. F., Gaidzik, V. I., Habdank, M., Krönke, J., Bullinger, L., et al. (2010). IDH1 and IDH2 mutations are frequent genetic alterations in acute myeloid leukemia and confer adverse prognosis in cytogenetically normal acute myeloid leukemia with NPM1 mutation without FLT3 internal tandem duplication. J. Clin. Oncol. 28, 3636–3643. doi: 10.1200/jco.2010.28.3762
Rao, S., Mondragón, L., Pranjic, B., Hanada, T., Stoll, G., Köcher, T., et al. (2019). AIF-regulated oxidative phosphorylation supports lung cancer development. Cell Res. 29, 579–591. doi: 10.1038/s41422-019-0181-4
Reznik, E., Miller, M. L., Şenbabaoğlu, Y., Riaz, N., Sarungbam, J., Tickoo, S. K., et al. (2016). Mitochondrial DNA copy number variation across human cancers. Elife 5:e10769. doi: 10.7554/eLife.10769
Schon, E. A., and Przedborski, S. (2011). Mitochondria: the next (neurode)generation. Neuron 70, 1033–1053. doi: 10.1016/j.neuron.2011.06.003
Sotgia, F., and Lisanti, M. (2017). Mitochondrial markers predict survival and progression in non-small cell lung cancer (NSCLC) patients: use as companion diagnostics. Oncotarget 8, 68095–68107. doi: 10.18632/oncotarget.19677
Su, C. Y., Chang, Y. C., Yang, C. J., Huang, M. S., and Hsiao, M. (2016). The opposite prognostic effect of NDUFS1 and NDUFS8 in lung cancer reflects the oncojanus role of mitochondrial complex I. Sci. Rep. 6:31357. doi: 10.1038/srep31357
Telang, S., Nelson, K., Siow, D., Yalcin, A., Thornburg, J., Imbert-Fernandez, Y., et al. (2012). Cytochrome c oxidase is activated by the oncoprotein Ras and is required for A549 lung adenocarcinoma growth. Mol. Cancer 11:60. doi: 10.1186/1476-4598-11-60
Tomlinson, I. P., Alam, N. A., Rowan, A. J., Barclay, E., Jaeger, E. E., Kelsell, D., et al. (2002). Germline mutations in FH predispose to dominantly inherited uterine fibroids, skin leiomyomata and papillary renal cell cancer. Nat. Genet. 30, 406–410. doi: 10.1038/ng849
Tumeh, P., Harview, C., Yearley, J., Shintaku, I., Taylor, E., Robert, L., et al. (2014). PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature 515, 568–571. doi: 10.1038/nature13954
Tyagi, A., Pramanik, R., Vishnubhatla, S., Ali, S., Bakhshi, R., Chopra, A., et al. (2018). Pattern of mitochondrial D-loop variations and their relation with mitochondrial encoded genes in pediatric acute myeloid leukemia. Mutat. Res. 810, 13–18. doi: 10.1016/j.mrfmmm.2018.05.002
van Nederveen, F. H., Korpershoek, E., Lenders, J. W., de Krijger, R. R., and Dinjens, W. N. (2007). Somatic SDHB mutation in an extraadrenal pheochromocytoma. N. Engl. J. Med. 357, 306–308. doi: 10.1056/NEJMc070010
Viale, A., Corti, D., and Draetta, G. F. (2015). Tumors and mitochondrial respiration: a neglected connection. Cancer Res. 75, 3685–3686. doi: 10.1158/0008-5472.can-15-0491
Vyas, S., Zaganjor, E., and Haigis, M. C. (2016). Mitochondria and Cancer. Cell 166, 555–566. doi: 10.1016/j.cell.2016.07.002
Wang, X., Shi, D., Zhao, D., and Hu, D. (2020). Aberrant methylation and differential expression of SLC2A1, TNS4, GAPDH, ATP8A2, and CASZ1 are associated with the prognosis of lung adenocarcinoma. Biomed. Res. Int. 2020:1807089. doi: 10.1155/2020/1807089
Wang, Z., Choi, S., Lee, J., Huang, Y. T., Chen, F., Zhao, Y., et al. (2015). Mitochondrial variations in non-small cell lung cancer (NSCLC) survival. Cancer Inform. 14, 1–9.
Weerts, M. J. A., Smid, M., Foekens, J. A., Sleijfer, S., and Martens, J. W. M. (2018). Mitochondrial RNA expression and single nucleotide variants in association with clinical parameters in primary breast cancers. Cancers (Basel) 10:500. doi: 10.3390/cancers10120500
Witherspoon, M., Sandu, D., Lu, C., Wang, K., Edwards, R., Yeung, A., et al. (2019). ETHE1 overexpression promotes SIRT1 and PGC1α mediated aerobic glycolysis, oxidative phosphorylation, mitochondrial biogenesis and colorectal cancer. Oncotarget 10, 4004–4017. doi: 10.18632/oncotarget.26958
Yu, L., Lu, M., Jia, D., Ma, J., Ben-Jacob, E., Levine, H., et al. (2017). Modeling the genetic regulation of cancer metabolism: interplay between glycolysis and oxidative phosphorylation. Cancer Res. 77, 1564–1574. doi: 10.1158/0008-5472.can-16-2074
Yuan, Y., Wang, W., Li, H., Yu, Y., Tao, J., Huang, S., et al. (2015). Nonsense and missense mutation of mitochondrial ND6 gene promotes cell migration and invasion in human lung adenocarcinoma. BMC Cancer 15:346. doi: 10.1186/s12885-015-1349-z
Zhang, C., Lin, X., Zhao, Q., Wang, Y., Jiang, F., Ji, C., et al. (2020). YARS as an oncogenic protein that promotes gastric cancer progression through activating PI3K-Akt signaling. J. Cancer Res. Clin. Oncol. 146, 329–342. doi: 10.1007/s00432-019-03115-7
Keywords: lung adenocarcinoma, nuclear mitochondrial genes, risk score, overall survival, signature
Citation: Zhang X, Dong W, Zhang J, Liu W, Yin J, Shi D and Ma W (2021) A Novel Mitochondrial-Related Nuclear Gene Signature Predicts Overall Survival of Lung Adenocarcinoma Patients. Front. Cell Dev. Biol. 9:740487. doi: 10.3389/fcell.2021.740487
Received: 13 July 2021; Accepted: 28 September 2021;
Published: 25 October 2021.
Edited by:
Giampaolo Morciano, University of Ferrara, ItalyReviewed by:
Daniela Ramaccini, University of Ferrara, ItalyNataliia Naumova, University of Padua, Italy
Copyright © 2021 Zhang, Dong, Zhang, Liu, Yin, Shi and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Ma, Mw5088@163.com
†These authors have contributed equally to this work and share first authorship