Next Article in Journal
Myeloperoxidase PET Imaging Tracks Intracellular and Extracellular Treatment Changes in Experimental Myocardial Infarction
Next Article in Special Issue
Systemic Analyses of Cuproptosis-Related lncRNAs in Pancreatic Adenocarcinoma, with a Focus on the Molecular Mechanism of LINC00853
Previous Article in Journal
Catalytic Degradation of Ciprofloxacin in Aqueous Solution by Peroxymonosulfate Activated with a Magnetic CuFe2O4@Biochar Composite
Previous Article in Special Issue
Proteotranscriptomic Discrimination of Tumor and Normal Tissues in Renal Cell Carcinoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 6
10.3390/ijms24065705
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Expression of COA6 Is Related to Unfavorable Prognosis and Enhanced Oxidative Phosphorylation in Lung Adenocarcinoma

by
Ming Zhang
1,2,3,†,
Xiaohua Liao
1,2,†,
Guanxu Ji
1,2,
Xianming Fan
3,* and
Qiang Wu
1,2,*
1
Faculty of Chinese Medicine, Macau University of Science and Technology, Macao 999078, China
2
State Key Laboratory of Quality Research in Chinese Medicine, Macao 999078, China
3
Department of Respiratory and Critical Care Medicine, Affiliated Hospital of Southwest Medical University, Luzhou 646000, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(6), 5705; https://doi.org/10.3390/ijms24065705
Submission received: 7 November 2022 / Revised: 5 March 2023 / Accepted: 7 March 2023 / Published: 16 March 2023
(This article belongs to the Special Issue Biomarkers of Tumor Progression, Prognosis and Therapy)
Browse Figures
Versions Notes

Abstract

:
Mitochondrial metabolism plays an important role in the occurrence and development of cancers. Cytochrome C oxidase assembly factor six (COA6) is essential in mitochondrial metabolism. However, the role of COA6 in lung adenocarcinoma (LUAD) remains unknown. Here we report that the expression of COA6 mRNA and protein were upregulated in LUAD tissues compared with lung normal tissues. We found that COA6 had high sensitivity and specificity to distinguish LUAD tissues from normal lung tissues shown by a receiver operating characteristic (ROC) curve. In addition, our univariate and multivariate Cox regression analysis indicated that COA6 was an independent unfavorable prognostic factor for LUAD patients. Furthermore, our survival analysis and nomogram showed that a high expression of COA6 mRNA was related to the short overall survival (OS) of LUAD patients. Notably, our weighted correlation network analysis (WGCNA) and functional enrichment analysis revealed that COA6 may participate in the development of LUAD by affecting mitochondrial oxidative phosphorylation (OXPHOS). Importantly, we demonstrated that depletion of COA6 could decrease the mitochondrial membrane potential (MMP), nicotinamide adenine dinucleotide (NAD) + hydrogen (H) (NADH), and adenosine triphosphate (ATP) production in LUAD cells (A549 and H1975), hence inhibiting the proliferation of these cells in vitro. Together, our study strongly suggests that COA6 is significantly associated with the prognosis and OXPHOS in LUAD. Hence, COA6 is highly likely a novel prognostic biomarker and therapeutic target of LUAD.

1. Introduction

Lung cancer is a common malignant tumor worldwide. In 2020, more than two million new lung cancer cases occurred globally and the number of deaths was the highest among all cancers, up from 1.8 million cases [1]. Nonsmall cell lung cancer (NSCLC) accounts for about 85% of all lung cancers [2] and lung adenocarcinoma (LUAD) is the most common histological subtype of NSCLC [3]. Current therapies for LUAD include surgery, chemotherapy, radiotherapy, molecular targeted therapy, and immunotherapy [4]. LUAD patients usually lack symptoms in the early stage and many patients are in the middle and late stages when their symptoms appear. Hence, early diagnosis can significantly reduce the mortality of LUAD patients [5]. Although some progress has been made in the treatment of advanced patients, drug resistance often occurs, leading to treatment failure. In China, the 5-year survival rate of patients with lung cancer is less than 30% [6]. Therefore, identifying new biomarkers and therapeutic targets is important for the diagnosis and therapy of LUAD.
Mitochondrial metabolism is important for the occurrence and development of lung cancer. By using radiotracers and positron emission tomography to measure the mitochondrial membrane potential (MMP) in a mouse model of lung cancer, Momcilovic M et al. discovered that LUAD tissues possessed high MMP [7]. Analysis with data from The Cancer Genome Atlas (TCGA) indicated that the content of mitochondrial DNA in lung cancer tissues is higher than that in adjacent normal tissues [8]. Interestingly, metformin could reduce the production of cellular ATP by inhibiting mitochondrial respiratory chain complex Ⅰ [9]. In addition, in phase II clinical trial, metformin combined with an epidermal growth factor receptor-tyrosine kinase inhibitor significantly improved progression-free survival and overall survival (OS) of patients with advanced LUAD [10].
The cytochrome C oxidase assembly factor six (COA6) protein is a 14 kDa protein with a conservative motif CX9CXnCX10C. COA6 is mainly located in the inner mitochondrial membrane space. Importantly, COA6 is involved in the formation of cytochrome C oxidase (COX) [11]. Mutation or lack of COA6 can cause COX deficiency, adenosine triphosphate (ATP) synthesis disorder, and fatal neonatal cardiomyopathy [12]. These suggest that COA6 plays an important role in mitochondrial metabolism, yet the role of COA6 in LUAD remains unclear.
In this research, we found that COA6 was overexpressed in LUAD tissues. We then investigated the impact of COA6 on the diagnosis and prognosis of LUAD patients. Furthermore, we discussed the underlying mechanism of COA6 in LUAD. Finally, we verified the effect of COA6 on mitochondrial function and the proliferation of lung cancer cells in vitro.

2. Results

2.1. The Expression of COA6 Was Upregulated in LUAD

To probe whether COA6 is a potential marker of LUAD, we first examined the COA6 expression in LUAD from multiple databases. We found that the level of COA6 mRNA expression was higher in LUAD tissues than in normal lung tissues both in ONCOMINE and TCGA databases (p < 0.05) (Table 1 and Figure 1A). Our receiver operating characteristic (ROC) curve showed that COA6 had high sensitivity and specificity to distinguish LUAD and lung normal tissues in the TCGA database (area under the ROC curve (AUC) = 0.903, p < 0.0001) (Figure 1B). The level of COA6 protein expression was also higher in LUAD tissues than normal lung tissues in the Clinical Proteomic Tumor Analysis Consortium (CPTAC) database (p < 0.05) (Figure 1C). In the Human Protein Atlas (HPA), as LUAD tissues were from five patients, and lung normal tissues were from two patients, it is difficult to compare the expression of COA6 protein in LUAD tissue and normal lung tissue with statistical significance due to the small sample size. We also observed that COA6 was located in the cytoplasm and membrane. (Figure 1D,E and Table S1).

2.2. Associations of COA6 Expression Level with Clinical Parameters in LUAD

We then compared the COA6 mRNA expression in different clinical subgroups of LUAD patients in TCGA (Table S2). In age and gender subgroups, there were no significant differences in COA6 mRNA expression (Figure 2A,B). In the T (the size of the tumor) stages, the COA6 mRNA expression was higher in T2, T3, and T4 than in the T1 groups (p < 0.05) (Figure 2C). In N (lymph node metastasis) stages, the COA6 mRNA expression was higher in N1/2/3 than in N0 groups (p < 0.05) (Figure 2D). In M (distant metastasis) stages, the COA6 mRNA expression was higher in the M1 than in the M0 groups (p < 0.001) (Figure 2E). In tumor (the size of the tumor, lymph node metastasis, and distant metastasis) stages, the COA6 mRNA expression was higher in the II, III, and IV groups than the I group (p < 0.05) (Figure 2F). These data strongly suggest that the higher level of COA6 is related to the late stages of LUAD.

2.3. Prognostic Impact of COA6 Expression in LUAD

To investigate whether COA6 expression can be a potential prognostic marker in LUAD, we performed bioinformatic analysis using datasets from TCGA (485 patients) and GSE31210 (204 patients) (Table 2). Our univariate and multiple Cox regression analysis of these two databases indicated that COA6 and tumor stage were independent unfavorable prognostic factors of LUAD patients (Figure 3A–D). Notably, our survival analysis showed that patients with a higher level of COA6 had a shorter OS (Figure 3E,F). In TCGA, we built a nomogram based on the two independent prognostic factors (COA6 and tumor stage) for OS of LUAD patients and found that the increase of COA6 is related to a reduction of the survival rate of patients (Figure 4A). Therefore, our ROC curve showed certain specificity and sensitivity of the nomogram for prognosis (AUC of 1, 2, 3-year survival: 0.732, 0.705, 0.701) (Figure 4B). Furthermore, our calibration curve showed a good prediction accuracy of the nomogram for one, three, and five-year survival (Figure 4C). The C-index of the nomogram based on COA6 and tumor stage was 0.684 (0.661–0.706), and the c-index of a nomogram based only on tumor stage was 0.664 (0.643–0.684). Taken together, we suggest that an increased COA6 level could be used to predict the prognosis of LUAD patients.

2.4. Difference of Biological Pathways between High and Low COA6 Groups in LUAD

Next, to identify the difference between biological pathways of high and low COA6 groups, we compared global gene expression differences between the two groups in TCGA. We found that there were a total of 19,101 genes with log2foldchange (Figure 5A,B and Table S3). We then performed a GSEA with the rank of log2foldchange. We discovered that there were 19 activated and 6 suppressed hallmark pathways in the high COA6 group if compared with the low COA6 group, (normalized enrichment score (NES)|>1, p.adjust < 0.05) (Figure 5C,D and Table S4). Notably, in the high COA6 group, the main activated pathways were OXPHOS, MYC targets, and mTORC1 signaling, whereas, the main suppressed pathways were UV RESPONSE DN, myogenesis, and IL6-JAK-STAT3 signaling.

2.5. Biological Pathways Correlated with COA6 in LUAD

To find key module eigengenes related to COA6, we performed a WGCNA with the data of 485 LUAD patients from TCGA. We identified 4802 genes (top 25% of all genes according to variance) extracted by WGCNA (Table S5). We used clinical parameters and COA6 to construct a trait heatmap (Figure 6A). The soft threshold power was set as six (scale R^2 = 0.92) (Figure 6B) and mergeCutHeight was set as 0.25 (Figure 6C). We identified twenty modules (Figure 6D), in which the brown modules were most correlated with COA6 (Pearson coefficient = 0.47 and p = 9 × 10−27) (Figure 6E). We found 341 brown module eigengenes in a scatter plot (Figure 6F and Table S6). We also performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis with the brown module eigengenes. We found that these eigengenes were mainly correlated with the cellular component (CC) of the ribosome and mitochondrial inner membrane. Their molecular functions (MF) were focused on related electron transfer activity, cytochrome-c oxidase activity, and oxidoreductase activity. Moreover, the biological processes (BP) in which they participate included ribosome biogenesis and oxidative phosphorylation. Finally, the KEGG pathway closely associated with them contained ribosome, oxidative phosphorylation, chemical carcinogenesis-reactive oxygen species, and so on (Figure 7A–D and Tables S7 and S8). These data clearly demonstrate that COA6 plays essential roles in mitochondrial metabolism.

2.6. Down-Regulation of COA6 Inhibited the Proliferation of LUAD

To validate that COA6 is crucial in mitochondrial metabolism and LUAD, we used siRNA against COA6 to knock down the expression of COA6 in LUAD cells (A549 and H1975). Our qRT-PCR results showed that the COA6 mRNA level was significantly decreased by COA6 siRNA (p < 0.001) (Figure 8A). Our Western blot further confirmed that the protein levels of COA6, together with COX2, were decreased by COA6 siRNA in the LUAD cells (Figure 8B,C). In addition, our immunofluorescence affirmed that the expression of COA6 protein, which was mainly expressed in the cytoplasm in LUAD cells, was decreased upon COA6 RNAi (Figure 8D). Importantly, our MTT analysis indicated that the knockdown of COA6 significantly reduced the proliferation rate of LUAD cells (Figure 8E). Our colony formation assays also showed that the knockdown of COA6 could inhibit the formation of a colony of LUAD cells (Figure 8F,G). These results suggest that the knockdown of COA6 could decrease COX2 protein level and inhibit the proliferation of LUAD cells in vitro.

2.7. Down Regulation of COA6 Affected the MMP, Reactive Oxygen Species (ROS), and Nicotinamide Adenine Dinucleotide (NAD) + Hydrogen (H) (NADH) of LUAD

To further investigate the role of COA6 in mitochondrial metabolism, we measured mitochondrial membrane potential and intracellular ROS by the fluorescent probe. Our results indicated that knockdown of COA6 significantly decreased the MMP of LUAD cells (p < 0.05) (Figure 9A,C,D) and also inhibited the production of ROS of H1975 cells (p < 0.05) since the ROS level in A549 cells was lower in the COA6 siRNA group than in the control siRNA group, although without significant difference between the two groups (Figure 9B,E,F). We also observed that the NADH/NAD+ ratio in the COA6 siRNA group was lower than in the control siRNA group (p < 0.001) (Figure 9B,G,H). These suggest that the knockdown of COA6 can significantly reduce NADH production in LUAD cells.

2.8. Down Regulation of COA6 Deceased the Production of ATP in LUAD

We performed a real-time ATP rate assay to uncover whether COA6 has a role in ATP production, We found that the total ATP and the ATP produced by mitochondrial OXPHOS were lower in the COA6 siRNA group than in the control siRNA group (p < 0.05), whereas there was no significant difference of ATP produced by glycolysis between the two groups (Figure 10A–D). In contrast, in H1975 cells, the total ATP, ATP produced by mitochondrial OXPHOS, and ATP produced by glycolysis were all lower in the COA6 siRNA group than in the control siRNA group (p < 0.05) (Figure 10E–H). These results suggest that the knockdown of COA6 could decrease the production of ATP in LUAD cells in vitro.

3. Discussion

In this research, we first reported that both COA6 mRNA and protein expression were overexpressed in LUAD and COA6 had high sensitivity and specificity to distinguish LUAD tissues and normal lung tissues. In addition, the expression of COA6 was correlated with the clinical parameters of LUAD patients. Our analysis indicated that COA6 is an independent risk factor for the prognosis of LUAD patients. Patients with a higher expression of COA6 had a shorter OS. These strongly suggest that COA6 is an important gene for the development of LUAD.
As a cytochrome C oxidase assembly factor, COA6 can interact with the copper-containing catalytic domain to maintain the stability of the newly synthesized COX subunit two (COX2) structure [13]. Loss of COA6 can lead to the rapid degradation of newly synthesized COX2. Subsequently, the level of COX decreases and the mitochondrial metabolism is impaired [14]. Since COX is essential for the process of OXPHOS as a terminal enzyme in the respiratory chain, a lack of COX usually leads to severe early-onset neuromuscular diseases and cardiomyopathy [15,16,17,18]. Current research on COA6 has mainly been concentrated on its role in mitochondrial metabolic diseases. Calvo SE et al. found that the mutation of COA6 was complicated by combined deficiency of mitochondrial complexes I and IV in the myocardial tissue of a patient with hypertrophic cardiomyopathy [19]. Subsequent studies confirmed that the deficiency of COA6 can disrupt COX assembly, reduce mitochondrial OXPHOS, and cause neonatal hypertrophic cardiomyopathy. Copper treatment can repair COX defects partially [12,20]. Recently, it has been proved that COA6 mainly acted as a thiol reductase in maintaining the structure’s stability of COX2, reducing the disulfide bond of the key cysteine residues in the copper metal chaperones (especially Sco1 and Sco2), and participating in the transport of copper [21,22]. Consistent with these findings, our results showed that the genes correlated with COA6 were mainly concentrated in the mitochondria, related to the activities of oxidoreductase, NADH dehydrogenase, and COX. In KEGG enrichment analysis, these genes were mainly related to oxidative phosphorylation. These strongly indicate that COA6 also has a crucial role in mitochondrial metabolism in LUAD.
Glycolysis was once considered to be the main metabolic pathway for the rapid proliferation of cancer cells, while mitochondrial metabolism was dispensable [23]. However, in recent years, research has proven that mitochondrial metabolism plays an important role in the occurrence and development of tumors. In mouse lung cancer driven by the KRAS gene, ROS, which is required for the growth of cancer, is mainly produced in mitochondrial metabolism and the deletion of mitochondrial transcription factor A could disrupt mitochondrial function and decrease tumor incidence [24]. In patients with lung cancer, tumor tissues had higher levels of glucose oxidation and tricarboxylic acid cycle than adjacent normal tissues [25]. In mouse lung cancer models, tumor tissues had a high MMP level, suggesting a high level of OXPHOS [7]. Decreased expression of heat shock protein family D member one could reduce mitochondrial OXPHOS, leading to inhibition of the growth of A549 cells [26]. Consistently, our research found that higher levels of COA6 mRNA and protein in LUAD tissues and a large number of genes positively correlated with COA6 were enriched in the OXPHOS gene set, suggesting that COA6 is closely related to OXPHOS in LUAD. The OXPHOS pathway was significantly activated in the high COA6 group. These may explain why the expression of COA6 was associated with tumor stage and T stage in LUAD and patients with a higher expression of COA6 were of shorter OS. Our in vitro data also affirmed that the depletion of COA6 could decrease the mitochondrial membrane potential (MMP), NADH, and ATP production of LUAD cells, thus inhibiting the proliferation of these cells. Together, we propose that the uplifted COA6 level is closely related to the enhanced OXPHOS, which is beneficial to tumor growth and leads to the poor prognosis of LUAD patients.
Our gene-set enrichment analysis (GSEA) also showed COA6 is associated with MYC targets, ROS, mTORC1 signaling, glycolysis, E2F targets, DNA repair, and G2M checkpoints. As a transcription factor, MYC can enhance the growth and survival of lung cancer and other tumors [27,28,29]. Interestingly, MYC is important to activate mitochondrial biogenesis in cancer and could promote the growth of LUAD cells by regulating the structure and function of mitochondria [30,31]. Thus, we propose that COA6 is possibly involved in MYC-mediated tumor promotion. Intracellular ROS are mainly produced in mitochondrial metabolism and there are eight sites in the mitochondria that can produce superoxide [32,33]. ROS can activate the PI3K pathway by inhibiting PTEN (a PtdIns(3,4,5)P3 phosphatase) and increase the activity of AKT. Activation of the PI3K–AKT–mTOR pathway affects many cell biological processes, such as cell growth, proliferation, and differentiation [34,35]. For instance, activation of the mTOR pathway is related to the clinical outcome, invasiveness, and metastasis of LUAD [36]. Under hypoxia, ROS could activate hypoxia-inducible factors (HIFs) which can increase glycolysis and activate angiogenesis genes, and then promote the growth of the tumor [37,38]. In addition, E2F targets, DNA repair, and G2M checkpoints are important signal pathways in tumorigenesis and the development of LUAD [39,40,41]. Our result also showed that depletion of COA6 could decrease glycolysis-produced ATP and ROS in H1975 cells. Therefore, we suggest that COA6 may promote the proliferation of LUAD cells by multiple tumor promotion pathways including MYC, PI3K–AKT–mTOR, and ROS.
In conclusion, we demonstrated that COA6 is important in regulating OXPHOS. LUAD was significantly affected by COA6 and depletion of COA6 can inhibit OXPHOS and decrease ATP production, leading to inhibited proliferation of LUAD cells. Further clinical research can be carried out to affirm the functions of COA6 in mitochondrial metabolism in LUAD. Nonetheless, our COA6 study has provided novel insights into the roles of COA6 in cancer progression.

4. Materials and Methods

4.1. Bioinformatics Analysis

ONCOMINE [42] was used to compare the COA6 mRNA expression of LUAD tissues and lung normal tissues. There were four datasets containing gene expression data of LUAD patients, respectively, from Stanford (Garber Lung [43]), Erasmus (Hou Lung [44]), Tokyo (Okayama Lung [45]), and Southern California (Selamat Lung [46]). UALCAN [47] was also used to compare their COA6 protein expression of them. The HPA [48] was used to obtain the immunohistochemical results of COA6 protein in LUAD and lung normal tissues.
To further analyze the relationship between COA6 and LUAD, we downloaded the RNA-seq and clinical data of LUAD patients from TCGA (https://portal.gdc.cancer.gov, accessed on 20 October 2021) and GSE31210 (http://www.prognoscan.org, accessed on 20 April 2022). RNA-seq data of 519 LUAD tissues and 58 adjacent normal tissues from TCGA were used to verify the results in ONCOMINE and plot the ROC curve. Patients’ inclusion criteria for clinically relevant analysis: the patients with RNA-seq data, complete clinical stages (T stage, N stage, M stage, and tumor stage), and survival data (OS and survival status). Then we compared the COA6 mRNA expression of different clinical subgroups in 485 LUAD patients from TCGA.
In TCGA and GSE31210, the impact of COA6 on patients’ OS was determined by univariate and multivariate Cox regression analysis. The patients were divided into a low COA6 group and a high COA6 group (the median of COA6 mRNA expression was set as the group cutoff). Then, the OS of the two groups were compared by Kaplan–Meier analysis. A prognostic model was built using the independent prognostic factors of LUAD patients in TCGA through COX multivariate regression. The model was displayed using a nomogram. The specificity and sensitivity of the model were shown by the ROC curve. The accuracy of the model was evaluated by calibration curve and C-index.
Differentially expressed genes between the low COA6 and high COA6 groups from TCGA were used for GSEA, and hallmark gene sets (h.all.v7.4.symbols.gmt) were downloaded from the Molecular Signatures Database [49]. To find the module eigengenes most correlated with COA6, the data of 485 LUAD patients from TCGA were used for WGCNA. Then, functional enrichment analysis of these module eigengenes was performed.

4.2. Cell Culture and Transfection

LUAD cells (A549, H1975) were obtained from the cell bank of the Chinese Academy of Sciences. The cells were cultured in 1640 medium (GIBCO, Waltham, MA, USA), added with 10% fetal bovine serum (GIBCO, Waltham, MA, USA), and cultured in a humidified 5% CO2 incubator at 37 °C. lipornai™ transfection reagent (Beyotime, Shanghai, China) was used to transfect the control and COA6 siRNA (Santa Cruz, Dallas, TX, USA) into cells. After 48 h from transfection, the effect of transfection was determined by qRT-PCR, western blot analysis, and immunofluorescence analysis.

4.3. RNA Isolation and qRT-PCR

Total RNA isolation, reverse transcription, and qRT-PCR were performed according to the product manual. The main kits included TransZol Up Plus RNA Kit (Transgen Biotech, Beijing, China), All-in-one First-Strand cDNA Synthesis SuperMix for qPCR (Transgen Biotech, Beijing, China), and Green qPCR SuperMix (Transgen Biotech, Beijing, China). The primer sequences for qRT-PCR were as follows: COA6: forward (5′-CAGTAGGAATGGCAGCCCCATCT-3′) and reverse (5′-CCCCCAGCAGACCTGTCTTTCCTT-3′); β-actin: forward (5′-CAGGGCGTGATGGTGGGCAT-3′) and reverse (5′-GATGCCGTGCTCGATGGGGT-3′).

4.4. Protein Extraction and Western Blot Analysis

Total protein was extracted using radioimmunoassay lysis buffer (RIPA, Thermo Fisher Scientific, Waltham, MA, USA). Protein assay dye reagent concentration (Bio Rad, Hercules, CA, USA) was used to measure the protein concentration. Western blot analysis was performed previously as described [50]. The primary antibodies were as follows: anti-COA6 (24209-1-AP, Proteintech, Wuhan, China), anti-COX2 (55070-1-AP, Proteintech, Wuhan, China), and anti-GAPDH (60004-1-Ig, Proteintech, Wuhan, China). Ultrasignal hypersensitive ECL chemiluminescence substrate (4A biotech, Suzhou, China) was used for chemiluminescence.

4.5. Cellular Immunofluorescence Assay

Cellular immunofluorescence assay was performed previously as described [51]. The cells were fixed with 4% formaldehyde and then infiltrated with 0.1% Triton X-100. Immunofluorescence stains of cells were performed with anti-COA6 antibody (24209-1-ap, Proteintech) and Alexa fluor 594 Goat antirabbit IgG (Thermo Fisher Scientific, Waltham, MA, USA). The nuclei were stained with 4′,6-diamidino-2-phenylindole.

4.6. MTT Assay and Colony Formation Experiment

Cells (1500 cells/well) were cultured in 96 well plates and transfected with siRNA. After 24, 48, 72, and 96 h, the cells were incubated with MTT solution for another 4 h. Cell viability was reflected by optical density (OD) at 490 nm. The transfected A549 and H1975 cells (1000/well) were seeded into 6-well plates. After 10 days, the cells were stained with crystal violet and the number of colonies was counted.

4.7. Detection of MMP, Intracellular ROS, and NADH/NAD+

Cells were cultured in 6 well plates and transfected with siRNA. After 72 h, MMP was labeled with Mito-Tracker Red CMXRos (Beyotime, Shanghai, China) and the intracellular ROS was labeled with DHE-ROS fluorescent probe (Bestbio, Shanghai, China). The cells were fixed with 4% formaldehyde and then infiltrated with 0.1% Triton X-100. The nuclei were stained with 4′,6-diamidino−2−phenylindole [52]. Cell number and fluorescence intensity were measured by Image-pro-plus. The intracellular NADH and NAD+ were detected with NAD+/NADH test kit (WST−8) (Beyotime, Shanghai, China). The concentration of NADH and NAD+ was reflected by optical density (OD) at 450 nm.

4.8. Real-Time ATP Rate Assay

Cells were cultured in 6 well plates and transfected with siRNA. After 48 h, plated A549 (12,000 cells/well) and H1975 (18,000 cells/well) in the Seahorse XFp Microplate, and 24 h later, Oligomycin and Rotenone/Antimycin A were used to affect the production of ATP. ATP rate was detected by Agilent Seahorse XFp according to the user manual. Cell number was estimated using Sulforhodamine B colorimetric assay [53].

4.9. Statistical Analysis

R (4.1.1) software and SPSS (25.0.0.0) software were used for statistical analyses, the main R packages used in the study were as follows: “ROCR”, “survival”, “survminer”, “DESeq2”, “clusterprofiler” and “WGCNA”. All in vitro experiments were repeated three times. A t-test and analysis of variance were used in comparing the means of two groups and multiple groups respectively, p < 0.05 means statistical significance.

5. Conclusions

In LUAD, COA6 is significantly correlated with the prognosis of patients and is important for mitochondrial metabolism. It is highly possible that COA6 is a novel prognostic biomarker and therapeutic target.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24065705/s1.

Author Contributions

Conceptualization, M.Z. and Q.W.; methodology, M.Z., X.L. and G.J.; software, M.Z.; validation, M.Z. and Q.W.; formal analysis, M.Z.; investigation, M.Z. and X.L.; resources, M.Z.; data curation, M.Z. and X.L.; writing—original draft preparation, M.Z. and X.L.; writing—review and editing, M.Z., X.F. and Q.W.; visualization, M.Z. and X.L.; supervision, X.F. and Q.W.; project administration, Q.W.; funding acquisition, Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Science and Technology Development Fund (FDCT Macau) (0072/2019/A2 and SKL-QRCM (MUST)-2023-2025).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Yao, G.; Chen, K.; Qin, Y.; Niu, Y.; Zhang, X.; Xu, S.; Zhang, C.; Feng, M.; Wang, K. Long Non-coding RNA JHDM1D-AS1 Interacts with DHX15 Protein to Enhance Non-Small-Cell Lung Cancer Growth and Metastasis. Mol. Ther. Nucleic. Acids 2019, 18, 831–840. [Google Scholar] [CrossRef] [Green Version]
  3. Mao, H. Clinical relevance of mutant-allele tumor heterogeneity and lung adenocarcinoma. Ann. Transl. Med. 2019, 7, 432. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, Z.; Yang, S.; Ma, Y.; Zhou, H.; Wu, X.; Han, J.; Hou, J.; Hao, L.; Spicer, J.D.; Koh, Y.W.; et al. Consistency of recommendations for the diagnosis and treatment of non-small cell lung cancer: A systematic review. Transl. Lung Cancer Res. 2021, 10, 2715–2732. [Google Scholar] [CrossRef]
  5. National Lung Screening Trial Research, T.; Church, T.R.; Black, W.C.; Aberle, D.R.; Berg, C.D.; Clingan, K.L.; Duan, F.; Fagerstrom, R.M.; Gareen, I.F.; Gierada, D.S.; et al. Results of initial low-dose computed tomographic screening for lung cancer. N. Engl. J. Med. 2013, 368, 1980–1991. [Google Scholar] [CrossRef] [Green Version]
  6. Allemani, C.; Matsuda, T.; Di Carlo, V.; Harewood, R.; Matz, M.; Niksic, M.; Bonaventure, A.; Valkov, M.; Johnson, C.J.; Esteve, J.; et al. Global surveillance of trends in cancer survival 2000-14 (CONCORD-3): Analysis of individual records for 37 513 025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries. Lancet 2018, 391, 1023–1075. [Google Scholar] [CrossRef] [Green Version]
  7. Momcilovic, M.; Jones, A.; Bailey, S.T.; Waldmann, C.M.; Li, R.; Lee, J.T.; Abdelhady, G.; Gomez, A.; Holloway, T.; Schmid, E.; et al. In vivo imaging of mitochondrial membrane potential in non-small-cell lung cancer. Nature 2019, 575, 380–384. [Google Scholar] [CrossRef] [PubMed]
  8. Reznik, E.; Miller, M.L.; Senbabaoglu, Y.; Riaz, N.; Sarungbam, J.; Tickoo, S.K.; Al-Ahmadie, H.A.; Lee, W.; Seshan, V.E.; Hakimi, A.A.; et al. Mitochondrial DNA copy number variation across human cancers. Elife 2016, 5, e10769. [Google Scholar] [CrossRef]
  9. El-Mir, M.Y.; Nogueira, V.; Fontaine, E.; Averet, N.; Rigoulet, M.; Leverve, X. Dimethylbiguanide inhibits cell respiration via an indirect effect targeted on the respiratory chain complex I. J. Biol. Chem. 2000, 275, 223–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Arrieta, O.; Barron, F.; Padilla, M.S.; Aviles-Salas, A.; Ramirez-Tirado, L.A.; Arguelles Jimenez, M.J.; Vergara, E.; Zatarain-Barron, Z.L.; Hernandez-Pedro, N.; Cardona, A.F.; et al. Effect of Metformin Plus Tyrosine Kinase Inhibitors Compared With Tyrosine Kinase Inhibitors Alone in Patients with Epidermal Growth Factor Receptor-Mutated Lung Adenocarcinoma: A Phase 2 Randomized Clinical Trial. JAMA Oncol. 2019, 5, e192553. [Google Scholar] [CrossRef]
  11. Vogtle, F.N.; Burkhart, J.M.; Rao, S.; Gerbeth, C.; Hinrichs, J.; Martinou, J.C.; Chacinska, A.; Sickmann, A.; Zahedi, R.P.; Meisinger, C. Intermembrane space proteome of yeast mitochondria. Mol. Cell Proteomics. 2012, 11, 1840–1852. [Google Scholar] [CrossRef] [Green Version]
  12. Baertling, F.; AM van den Brand, M.; Hertecant, J.L.; Al-Shamsi, A.; P van den Heuvel, L.; Distelmaier, F.; Mayatepek, E.; Smeitink, J.A.; Nijtmans, L.G.; Rodenburg, R.J. Mutations in COA6 cause cytochrome c oxidase deficiency and neonatal hypertrophic cardiomyopathy. Hum. Mutat. 2015, 36, 34–38. [Google Scholar] [CrossRef]
  13. Stroud, D.A.; Maher, M.J.; Lindau, C.; Vogtle, F.N.; Frazier, A.E.; Surgenor, E.; Mountford, H.; Singh, A.P.; Bonas, M.; Oeljeklaus, S.; et al. COA6 is a mitochondrial complex IV assembly factor critical for biogenesis of mtDNA-encoded COX2. Hum. Mol. Genet. 2015, 24, 5404–5415. [Google Scholar] [CrossRef] [Green Version]
  14. Pacheu-Grau, D.; Bareth, B.; Dudek, J.; Juris, L.; Vogtle, F.N.; Wissel, M.; Leary, S.C.; Dennerlein, S.; Rehling, P.; Deckers, M. Cooperation between COA6 and SCO2 in COX2 maturation during cytochrome c oxidase assembly links two mitochondrial cardiomyopathies. Cell Metab. 2015, 21, 823–833. [Google Scholar] [CrossRef] [Green Version]
  15. Bi, R.; Zhang, W.; Zhang, D.F.; Xu, M.; Fan, Y.; Hu, Q.X.; Jiang, H.Y.; Tan, L.; Li, T.; Fang, Y.; et al. Genetic association of the cytochrome c oxidase-related genes with Alzheimer’s disease in Han Chinese. Neuropsychopharmacology 2018, 43, 2264–2276. [Google Scholar] [CrossRef] [Green Version]
  16. Brischigliaro, M.; Zeviani, M. Cytochrome c oxidase deficiency. Biochim. Biophys. Acta Bioenerg. 2021, 1862, 148335. [Google Scholar] [CrossRef] [PubMed]
  17. Visuttijai, K.; Hedberg-Oldfors, C.; Lindgren, U.; Nordstrom, S.; Eliasdottir, O.; Lindberg, C.; Oldfors, A. Progressive external ophthalmoplegia associated with novel MT-TN mutations. Acta Neurol. Scand. 2021, 143, 103–108. [Google Scholar] [CrossRef] [PubMed]
  18. Wintjes, L.T.M.; Kava, M.; van den Brandt, F.A.; van den Brand, M.A.M.; Lapina, O.; Bliksrud, Y.T.; Kulseth, M.A.; Amundsen, S.S.; Selberg, T.R.; Ybema-Antoine, M.; et al. A novel variant in COX16 causes cytochrome c oxidase deficiency, severe fatal neonatal lactic acidosis, encephalopathy, cardiomyopathy, and liver dysfunction. Hum. Mutat. 2021, 42, 135–141. [Google Scholar] [CrossRef] [PubMed]
  19. Calvo, S.E.; Compton, A.G.; Hershman, S.G.; Lim, S.C.; Lieber, D.S.; Tucker, E.J.; Laskowski, A.; Garone, C.; Liu, S.; Jaffe, D.B.; et al. Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing. Sci. Transl. Med. 2012, 4, 118ra10. [Google Scholar] [CrossRef] [Green Version]
  20. Ghosh, A.; Trivedi, P.P.; Timbalia, S.A.; Griffin, A.T.; Rahn, J.J.; Chan, S.S.; Gohil, V.M. Copper supplementation restores cytochrome c oxidase assembly defect in a mitochondrial disease model of COA6 deficiency. Hum. Mol. Genet. 2014, 23, 3596–3606. [Google Scholar] [CrossRef]
  21. Soma, S.; Morgada, M.N.; Naik, M.T.; Boulet, A.; Roesler, A.A.; Dziuba, N.; Ghosh, A.; Yu, Q.; Lindahl, P.A.; Ames, J.B.; et al. COA6 Is Structurally Tuned to Function as a Thiol-Disulfide Oxidoreductase in Copper Delivery to Mitochondrial Cytochrome c Oxidase. Cell Rep. 2019, 29, 4114–4126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Pacheu-Grau, D.; Wasilewski, M.; Oeljeklaus, S.; Gibhardt, C.S.; Aich, A.; Chudenkova, M.; Dennerlein, S.; Deckers, M.; Bogeski, I.; Warscheid, B.; et al. COA6 Facilitates Cytochrome c Oxidase Biogenesis as Thiol-reductase for Copper Metallochaperones in Mitochondria. J. Mol. Biol. 2020, 432, 2067–2079. [Google Scholar] [CrossRef] [PubMed]
  23. Warburg, O. On the origin of cancer cells. Science 1956, 123, 309–314. [Google Scholar] [CrossRef] [PubMed]
  24. Weinberg, F.; Hamanaka, R.; Wheaton, W.W.; Weinberg, S.; Joseph, J.; Lopez, M.; Kalyanaraman, B.; Mutlu, G.M.; Budinger, G.R.; Chandel, N.S. Mitochondrial metabolism and ROS generation are essential for Kras-mediated tumorigenicity. Proc. Natl. Acad. Sci. USA 2010, 107, 8788–8793. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hensley, C.T.; Faubert, B.; Yuan, Q.; Lev-Cohain, N.; Jin, E.; Kim, J.; Jiang, L.; Ko, B.; Skelton, R.; Loudat, L.; et al. Metabolic Heterogeneity in Human Lung Tumors. Cell 2016, 164, 681–694. [Google Scholar] [CrossRef] [Green Version]
  26. Parma, B.; Ramesh, V.; Gollavilli, P.N.; Siddiqui, A.; Pinna, L.; Schwab, A.; Marschall, S.; Zhang, S.; Pilarsky, C.; Napoli, F.; et al. Metabolic impairment of non-small cell lung cancers by mitochondrial HSPD1 targeting. J. Exp. Clin. Cancer Res. 2021, 40, 248. [Google Scholar] [CrossRef]
  27. Chanvorachote, P.; Sriratanasak, N.; Nonpanya, N. C-myc Contributes to Malignancy of Lung Cancer: A Potential Anticancer Drug Target. Anticancer Res. 2020, 40, 609–618. [Google Scholar] [CrossRef]
  28. Elbadawy, M.; Usui, T.; Yamawaki, H.; Sasaki, K. Emerging Roles of C-Myc in Cancer Stem Cell-Related Signaling and Resistance to Cancer Chemotherapy: A Potential Therapeutic Target Against Colorectal Cancer. Int. J. Mol. Sci. 2019, 20, 2340. [Google Scholar] [CrossRef] [Green Version]
  29. Kalkat, M.; De Melo, J.; Hickman, K.A.; Lourenco, C.; Redel, C.; Resetca, D.; Tamachi, A.; Tu, W.B.; Penn, L.Z. MYC Deregulation in Primary Human Cancers. Genes 2017, 8, 151. [Google Scholar] [CrossRef] [Green Version]
  30. Li, F.; Wang, Y.; Zeller, K.I.; Potter, J.J.; Wonsey, D.R.; O’Donnell, K.A.; Kim, J.W.; Yustein, J.T.; Lee, L.A.; Dang, C.V. Myc stimulates nuclearly encoded mitochondrial genes and mitochondrial biogenesis. Mol. Cell Biol. 2005, 25, 6225–6234. [Google Scholar] [CrossRef] [Green Version]
  31. Graves, J.A.; Wang, Y.; Sims-Lucas, S.; Cherok, E.; Rothermund, K.; Branca, M.F.; Elster, J.; Beer-Stolz, D.; Van Houten, B.; Vockley, J.; et al. Mitochondrial structure, function and dynamics are temporally controlled by c-Myc. PLoS ONE 2012, 7, e37699. [Google Scholar] [CrossRef] [Green Version]
  32. Brand, M.D. The sites and topology of mitochondrial superoxide production. Exp. Gerontol. 2010, 45, 466–472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Quinlan, C.L.; Treberg, J.R.; Perevoshchikova, I.V.; Orr, A.L.; Brand, M.D. Native rates of superoxide production from multiple sites in isolated mitochondria measured using endogenous reporters. Free Radic. Biol. Med. 2012, 53, 1807–1817. [Google Scholar] [CrossRef] [Green Version]
  34. Leslie, N.R.; Bennett, D.; Lindsay, Y.E.; Stewart, H.; Gray, A.; Downes, C.P. Redox regulation of PI 3-kinase signalling via inactivation of PTEN. EMBO J. 2003, 22, 5501–5510. [Google Scholar] [CrossRef]
  35. Tamaddoni, A.; Mohammadi, E.; Sedaghat, F.; Qujeq, D.; As’Habi, A. The anticancer effects of curcumin via targeting the mammalian target of rapamycin complex 1 (mTORC1) signaling pathway. Pharmacol. Res. 2020, 156, 104798. [Google Scholar] [CrossRef]
  36. Krencz, I.; Sebestyen, A.; Fabian, K.; Mark, A.; Moldvay, J.; Khoor, A.; Kopper, L.; Papay, J. Expression of mTORC1/2-related proteins in primary and brain metastatic lung adenocarcinoma. Hum. Pathol. 2017, 62, 66–73. [Google Scholar] [CrossRef]
  37. Chandel, N.S.; McClintock, D.S.; Feliciano, C.E.; Wood, T.M.; Melendez, J.A.; Rodriguez, A.M.; Schumacker, P.T. Reactive oxygen species generated at mitochondrial complex III stabilize hypoxia-inducible factor-1alpha during hypoxia: A mechanism of O2 sensing. J. Biol. Chem. 2000, 275, 25130–25138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Semenza, G.L. Targeting HIF-1 for cancer therapy. Nat. Rev. Cancer 2003, 3, 721–732. [Google Scholar] [CrossRef] [PubMed]
  39. Hu, B.; Liu, D.; Liu, Y.; Li, Z. DNA Repair-Based Gene Expression Signature and Distinct Molecular Subtypes for Prediction of Clinical Outcomes in Lung Adenocarcinoma. Front. Med. Lausanne 2020, 7, 615981. [Google Scholar] [CrossRef]
  40. Zhang, Q.; Shi, R.; Bai, Y.; Meng, L.; Hu, J.; Zhu, H.; Liu, T.; De, X.; Wang, S.; Wang, J.; et al. Meiotic nuclear divisions 1 (MND1) fuels cell cycle progression by activating a KLF6/E2F1 positive feedback loop in lung adenocarcinoma. Cancer Commun. Lond. 2021, 41, 492–510. [Google Scholar] [CrossRef]
  41. Halike, X.; Li, J.; Yuan, P.; Yasheng, K.; Chen, M.; Xia, L.; Li, J. The petroleum ether extract of Brassica rapa L. induces apoptosis of lung adenocarcinoma cells via the mitochondria-dependent pathway. Food Funct. 2021, 12, 10023–10039. [Google Scholar] [CrossRef]
  42. Rhodes, D.R.; Yu, J.; Shanker, K.; Deshpande, N.; Varambally, R.; Ghosh, D.; Barrette, T.; Pandey, A.; Chinnaiyan, A.M. ONCOMINE: A cancer microarray database and integrated data-mining platform. Neoplasia 2004, 6, 1–6. [Google Scholar] [CrossRef] [Green Version]
  43. Garber, M.E.; Troyanskaya, O.G.; Schluens, K.; Petersen, S.; Thaesler, Z.; Pacyna-Gengelbach, M.; van de Rijn, M.; Rosen, G.D.; Perou, C.M.; Whyte, R.I.; et al. Diversity of gene expression in adenocarcinoma of the lung. Proc. Natl. Acad. Sci. USA 2001, 98, 13784–13789. [Google Scholar] [CrossRef] [Green Version]
  44. Hou, J.; Aerts, J.; den Hamer, B.; van Ijcken, W.; den Bakker, M.; Riegman, P.; van der Leest, C.; van der Spek, P.; Foekens, J.A.; Hoogsteden, H.C.; et al. Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS ONE 2010, 5, e10312. [Google Scholar] [CrossRef] [PubMed]
  45. Okayama, H.; Kohno, T.; Ishii, Y.; Shimada, Y.; Shiraishi, K.; Iwakawa, R.; Furuta, K.; Tsuta, K.; Shibata, T.; Yamamoto, S.; et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012, 72, 100–111. [Google Scholar] [CrossRef] [Green Version]
  46. Selamat, S.A.; Chung, B.S.; Girard, L.; Zhang, W.; Zhang, Y.; Campan, M.; Siegmund, K.D.; Koss, M.N.; Hagen, J.A.; Lam, W.L.; et al. Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 2012, 22, 1197–1211. [Google Scholar] [CrossRef] [Green Version]
  47. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef] [PubMed]
  48. Uhlen, M.; Zhang, C.; Lee, S.; Sjostedt, E.; Fagerberg, L.; Bidkhori, G.; Benfeitas, R.; Arif, M.; Liu, Z.; Edfors, F.; et al. A pathology atlas of the human cancer transcriptome. Science 2017, 357, eaan2507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Liberzon, A.; Birger, C.; Thorvaldsdottir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015, 1, 417–425. [Google Scholar] [CrossRef] [Green Version]
  50. Ji, G.; Xiao, X.; Huang, M.; Wu, Q. Jmjd6 regulates ES cell homeostasis and enhances reprogramming efficiency. Heliyon 2022, 8, e09105. [Google Scholar] [CrossRef]
  51. Ma, L.; Huang, M.; Liao, X.; Cai, X.; Wu, Q. NR2F2 Regulates Cell Proliferation and Immunomodulation in Whartons’ Jelly Stem Cells. Genes 2022, 13, 1458. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, H.; Liu, S.; Cai, Z.; Dong, W.; Ye, J.; Cai, Z.; Han, Z.; Liang, Y.; Zhuo, Y.; Luo, Y.; et al. Down-regulation of ACACA suppresses the malignant progression of Prostate Cancer through inhibiting mitochondrial potential. J. Cancer 2021, 12, 232–243. [Google Scholar] [CrossRef] [PubMed]
  53. Vichai, V.; Kirtikara, K. Sulforhodamine B colorimetric assay for cytotoxicity screening. Nat. Protoc. 2006, 1, 1112–1116. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 05705 g001 550
Figure 1. Expression of COA6 in LUAD and lung normal tissues. (A) COA6 mRNA expression in LUAD shown in TCGA; (B) ROC curve showing the diagnostic value of COA6 in LUAD shown in TCGA; (C) COA6 protein expression of LUAD in the University of Alabama at Birmingham cancer data analysis Portal (UALCAN); (D,E) Immunohistochemical staining of COA6 in normal lung tissues and LUAD tissues in HPA. COA6 antibody was labeled with DAB (3,3′-diaminobenzidine), the brown color indicates the expression of COA6 protein, magnification: 50× (upper panel), 200× (lower panel). t-test of independent samples was used in comparing the means of two groups, *** p < 0.001.
Figure 1. Expression of COA6 in LUAD and lung normal tissues. (A) COA6 mRNA expression in LUAD shown in TCGA; (B) ROC curve showing the diagnostic value of COA6 in LUAD shown in TCGA; (C) COA6 protein expression of LUAD in the University of Alabama at Birmingham cancer data analysis Portal (UALCAN); (D,E) Immunohistochemical staining of COA6 in normal lung tissues and LUAD tissues in HPA. COA6 antibody was labeled with DAB (3,3′-diaminobenzidine), the brown color indicates the expression of COA6 protein, magnification: 50× (upper panel), 200× (lower panel). t-test of independent samples was used in comparing the means of two groups, *** p < 0.001.
Ijms 24 05705 g001
Ijms 24 05705 g002 550
Figure 2. COA6 mRNA expression in different clinical subgroups of LUAD patients in TCGA. (A) Age; (B) Gender; (C) T (the size of the tumor) stage; (D) N (lymph node metastasis) stage; (E) M (distant metastasis) stage; (F) Tumor (the size of the tumor, lymph node metastasis, and distant metastasis) stage. t−test of independent samples was used in comparing the means of two groups, * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 2. COA6 mRNA expression in different clinical subgroups of LUAD patients in TCGA. (A) Age; (B) Gender; (C) T (the size of the tumor) stage; (D) N (lymph node metastasis) stage; (E) M (distant metastasis) stage; (F) Tumor (the size of the tumor, lymph node metastasis, and distant metastasis) stage. t−test of independent samples was used in comparing the means of two groups, * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 24 05705 g002
Ijms 24 05705 g003 550
Figure 3. Impact of COA6 on the prognosis of LUAD patients. (A,B) Univariate Cox regression analysis in TCGA and GSE31210; (C,D) Multivariate Cox regression analysis in TCGA and GSE31210. (E,F) Kaplan–Meier survival curve for impacts of COA6 on overall survival of LUAD patients in TCGA and GSE31210.
Figure 3. Impact of COA6 on the prognosis of LUAD patients. (A,B) Univariate Cox regression analysis in TCGA and GSE31210; (C,D) Multivariate Cox regression analysis in TCGA and GSE31210. (E,F) Kaplan–Meier survival curve for impacts of COA6 on overall survival of LUAD patients in TCGA and GSE31210.
Ijms 24 05705 g003
Ijms 24 05705 g004 550
Figure 4. Prognostic model based on COA6 and tumor stage for LUAD patients in TCGA. (A) Nomogram based on COA6 and tumor stage. (B) ROC curve respecting specificity and sensitivity of nomogram for prognosis. (C) Calibration curve of the nomogram for overall survival at 1 year, 3 years, and 5 years.
Figure 4. Prognostic model based on COA6 and tumor stage for LUAD patients in TCGA. (A) Nomogram based on COA6 and tumor stage. (B) ROC curve respecting specificity and sensitivity of nomogram for prognosis. (C) Calibration curve of the nomogram for overall survival at 1 year, 3 years, and 5 years.
Ijms 24 05705 g004
Ijms 24 05705 g005 550
Figure 5. GSEA of DEGs between high and low COA6 groups. (A) Volcano map of DEGs; (B) Heatmap of DEGs with significant differences (p.adjust < 0.01, log2foldchange > 3). (C) The top 5 activated hallmark pathways in the high COA6 group; (D) The top 5 suppressed hallmark pathways in the high COA6 group.
Figure 5. GSEA of DEGs between high and low COA6 groups. (A) Volcano map of DEGs; (B) Heatmap of DEGs with significant differences (p.adjust < 0.01, log2foldchange > 3). (C) The top 5 activated hallmark pathways in the high COA6 group; (D) The top 5 suppressed hallmark pathways in the high COA6 group.
Ijms 24 05705 g005
Ijms 24 05705 g006 550
Figure 6. Screening for key modules related to COA6 in TCGA- LUAD through WGCNA. (A) Sample dendrogram and trait heatmap; (B) Calculation of the scale-free fit index of various soft-thresholding powers; (C) Clustering of module eigengenes. The red line indicates mergeCutHeight (0.25); (D) Clustering dendrograms; (E) Correlation heatmap between module eigengenes and clinical parameters and COA6; (F) Scatter plot of brown module eigengenes.
Figure 6. Screening for key modules related to COA6 in TCGA- LUAD through WGCNA. (A) Sample dendrogram and trait heatmap; (B) Calculation of the scale-free fit index of various soft-thresholding powers; (C) Clustering of module eigengenes. The red line indicates mergeCutHeight (0.25); (D) Clustering dendrograms; (E) Correlation heatmap between module eigengenes and clinical parameters and COA6; (F) Scatter plot of brown module eigengenes.
Ijms 24 05705 g006
Ijms 24 05705 g007 550
Figure 7. Function enrichment plots of genes in the brown module. (A) The top 10 enriched GO−CC terms; (B) The top 10 enriched GO−MF terms; (C) The top 10 enriched GO−BP terms; (D) The top 10 enriched KEGG pathways.
Figure 7. Function enrichment plots of genes in the brown module. (A) The top 10 enriched GO−CC terms; (B) The top 10 enriched GO−MF terms; (C) The top 10 enriched GO−BP terms; (D) The top 10 enriched KEGG pathways.
Ijms 24 05705 g007
Ijms 24 05705 g008 550
Figure 8. COA6 regulates the proliferation of LUAD cells (A549 and H1975). (A) Histogram showing the known down of COA6 mRNA by siRNA (mean ± SD, n = 3); (B,C) Western blot showed that the decreased protein level of COA6 and COX2 by the knockdown of COA6 (mean ± SD, n = 3); (D) Immunofluorescence staining showing the inhibition of the COA6 protein expression, Magnification: 400×; (E) MTT showing the inhibited proliferation of cells by the knockdown of COA6, the optical density (OD) values of two groups were compared at the 96th hour (mean ± SD, n = 3); (F,G) Colony formation assay showing the reduced number of colonies by the knockdown of COA6 (mean ± SD, n = 3). t−test of independent samples was used in comparing the means of two groups, * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 8. COA6 regulates the proliferation of LUAD cells (A549 and H1975). (A) Histogram showing the known down of COA6 mRNA by siRNA (mean ± SD, n = 3); (B,C) Western blot showed that the decreased protein level of COA6 and COX2 by the knockdown of COA6 (mean ± SD, n = 3); (D) Immunofluorescence staining showing the inhibition of the COA6 protein expression, Magnification: 400×; (E) MTT showing the inhibited proliferation of cells by the knockdown of COA6, the optical density (OD) values of two groups were compared at the 96th hour (mean ± SD, n = 3); (F,G) Colony formation assay showing the reduced number of colonies by the knockdown of COA6 (mean ± SD, n = 3). t−test of independent samples was used in comparing the means of two groups, * p < 0.05, ** p < 0.01, and *** p < 0.001.
Ijms 24 05705 g008
Ijms 24 05705 g009 550
Figure 9. COA6 regulates the MMP, ROS, and NADH of LUAD cells (A549 and H1975). (A) Fluorescence staining showing the MMP affected by siRNA, Magnification: 200×; (B) Fluorescence staining showing the ROS affected by siRNA, Magnification: 200×; (C,D) The histograms showing that knockdown of COA6 decreased MMP of LUAD cells (mean ± SD, C: ncontrol = 11, ncoa6 = 13, D: ncontrol = 17, ncoa6 = 11); (E,F) The histograms showing that knockdown of COA6 did not decrease the ROS in A549 cells significantly, but decreased the ROS in H1975 cells significantly (mean ± SD, C: ncontrol = 30, ncoa6 = 32, D: ncontrol = 22, ncoa6 = 13); (G,H) The histograms showing that knockdown of COA6 decreased NADH in LUAD cells (mean ± SD, n = 3). Independent samples t−test was used in comparing the means of two groups, * p < 0.05 and *** p < 0.001.
Figure 9. COA6 regulates the MMP, ROS, and NADH of LUAD cells (A549 and H1975). (A) Fluorescence staining showing the MMP affected by siRNA, Magnification: 200×; (B) Fluorescence staining showing the ROS affected by siRNA, Magnification: 200×; (C,D) The histograms showing that knockdown of COA6 decreased MMP of LUAD cells (mean ± SD, C: ncontrol = 11, ncoa6 = 13, D: ncontrol = 17, ncoa6 = 11); (E,F) The histograms showing that knockdown of COA6 did not decrease the ROS in A549 cells significantly, but decreased the ROS in H1975 cells significantly (mean ± SD, C: ncontrol = 30, ncoa6 = 32, D: ncontrol = 22, ncoa6 = 13); (G,H) The histograms showing that knockdown of COA6 decreased NADH in LUAD cells (mean ± SD, n = 3). Independent samples t−test was used in comparing the means of two groups, * p < 0.05 and *** p < 0.001.
Ijms 24 05705 g009
Ijms 24 05705 g010 550
Figure 10. COA6 regulates the production of ATP in LUAD cells (A549 and H1975). (A,E) Line graphs showing the impact on OCR by OLI and Rtn/AA in different groups; (B,F) The histograms showing that knockdown of COA6 decreased ATP produced by mitochondrial OXPHOS (mean ± SD, n = 3); (C,G) The histograms showing that knockdown of COA6 did not decrease ATP produced by glycolysis significantly in A549 cells, but decreased ATP produced by glycolysis significantly in H1975 cells (mean ± SD, n = 3); (D,H) The histograms showing that knockdown of COA6 decreased total ATP in LUAD cells (mean ± SD, n = 3). OCR, oxygen consumption rate; OLI, oligomycin; Rtn/AA, Rotenone/Antimycin A. Independent samples t-test was used in comparing the means of two groups, * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 10. COA6 regulates the production of ATP in LUAD cells (A549 and H1975). (A,E) Line graphs showing the impact on OCR by OLI and Rtn/AA in different groups; (B,F) The histograms showing that knockdown of COA6 decreased ATP produced by mitochondrial OXPHOS (mean ± SD, n = 3); (C,G) The histograms showing that knockdown of COA6 did not decrease ATP produced by glycolysis significantly in A549 cells, but decreased ATP produced by glycolysis significantly in H1975 cells (mean ± SD, n = 3); (D,H) The histograms showing that knockdown of COA6 decreased total ATP in LUAD cells (mean ± SD, n = 3). OCR, oxygen consumption rate; OLI, oligomycin; Rtn/AA, Rotenone/Antimycin A. Independent samples t-test was used in comparing the means of two groups, * p < 0.05, ** p < 0.01, and *** p < 0.001.
Ijms 24 05705 g010
Table
Table 1. Significant difference between COA6 mRNA expression in LUAD tissues and normal lung tissues (ONCOMINE).
Table 1. Significant difference between COA6 mRNA expression in LUAD tissues and normal lung tissues (ONCOMINE).
DatasetLUAD CasesNormal CasesFold Changep-Valuet-Test
Garber Lung6401.830.0043.395
Hou Lung65451.6985.33 × 10−138.652
Okayama Lung202261.6733.08× 10−1210.077
Selamat Lung58581.1670.0022.954
Table
Table 2. Clinical data of the LUAD patients in TCGA and test GSE31210.
Table 2. Clinical data of the LUAD patients in TCGA and test GSE31210.
Clinical FactorTCGA (n = 485)GSE31210 (n = 204)
Survival time (days)868.99 ± 872.111764.70 ± 683.65
Status
Live307172
Dead17832
Age
<65216145
≥6526959
Gender
Female262109
Male22395
Tumor stage
I263162
II11642
III800
IV260
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, M.; Liao, X.; Ji, G.; Fan, X.; Wu, Q. High Expression of COA6 Is Related to Unfavorable Prognosis and Enhanced Oxidative Phosphorylation in Lung Adenocarcinoma. Int. J. Mol. Sci. 2023, 24, 5705. https://doi.org/10.3390/ijms24065705

AMA Style

Zhang M, Liao X, Ji G, Fan X, Wu Q. High Expression of COA6 Is Related to Unfavorable Prognosis and Enhanced Oxidative Phosphorylation in Lung Adenocarcinoma. International Journal of Molecular Sciences. 2023; 24(6):5705. https://doi.org/10.3390/ijms24065705

Chicago/Turabian Style

Zhang, Ming, Xiaohua Liao, Guanxu Ji, Xianming Fan, and Qiang Wu. 2023. "High Expression of COA6 Is Related to Unfavorable Prognosis and Enhanced Oxidative Phosphorylation in Lung Adenocarcinoma" International Journal of Molecular Sciences 24, no. 6: 5705. https://doi.org/10.3390/ijms24065705

APA Style

Zhang, M., Liao, X., Ji, G., Fan, X., & Wu, Q. (2023). High Expression of COA6 Is Related to Unfavorable Prognosis and Enhanced Oxidative Phosphorylation in Lung Adenocarcinoma. International Journal of Molecular Sciences, 24(6), 5705. https://doi.org/10.3390/ijms24065705

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop