- Open Access
Identification of COX4I2 as a hypoxia-associated gene acting through FGF1 to promote EMT and angiogenesis in CRC
Cellular & Molecular Biology Letters volume 27, Article number: 76 (2022)
Current evidence suggests that the hypoxic tumor microenvironment further aggravates tumor progression, leading to poor therapeutic outcomes. There is as yet no biomarker capable of evaluating the hypoxic state of the tumor. The cytochrome c oxidase (COX) subunit is crucial to the mitochondrial respiratory chain.
We investigated the potential oncogenic role of COX subunit 4 isoform 2 gene (COX4I2) in colorectal cancer (CRC) by least absolute shrinkage and selection operator (LASSO) and COX regression analysis to examine whether COX4I2 overexpression can predict colorectal cancer (CRC) prognosis. The association of COX4I2 levels with clinical features and its biological actions were evaluated both in vitro and in vivo.
Our analysis showed that elevated COX4I2 levels were correlated with poor clinical outcomes. We also observed that that COX4I2 may be involved in epithelial-mesenchymal transition, activation of cancer-related fibroblasts and angiogenesis in relation to fibroblast growth factor 1.
The COX4I2 level may be a predictor of outcome in CRC and may represent a novel target for treatment development.
Colorectal cancer (CRC) is a common digestive neoplasm, ranking third in incidence among malignancies and second in lethality among cancers . Its pathogenesis is complex, involving numerous predisposing factors [2,3,4], and has been linked to genetics, diet, and lifestyle habits. Statistics show that 30–50% of CRC tumors recur with distant metastases after surgery, especially to the liver, lung, bone and ovaries, and 5-year survival often does not exceed 10% .
Tumor growth and metastasis require vascular nourishment, and tumor neovascularization is a feature of the tumor microenvironment . The hypoxic microenvironment disrupts the balance of vascular-related factors, activates endothelial cells, and increases vascular permeability, thereby promoting the delivery of nutrients to cancer cells to enhance migration and proliferation, resulting in tumor progression [7, 8]. In recent years, anti-angiogenic therapy, such as bevacizumab and regorafenib, has been widely used in the treatment of CRC [9, 10]. However, the formation of tumor-associated blood vessels involves the action of multiple genes and signaling pathways, with wide-ranging and complex regulatory mechanisms, and relatively few patients can benefit from anti-angiogenic therapy. Thus, it is particularly important to screen for suitable target populations and prognostic markers associated with anti-angiogenic therapy.
Cytochrome c oxidase (COX) is the last enzyme in the respiratory chain and is responsible for the catalytic transfer of electrons to oxygen (O2) . The mitochondrial COX gene (COX) encodes the catalytic subunit that is complexed with nuclear gene-encoded structural subunits . Subunit 4 (COX4) is the largest of these and is involved in the mitochondrial respiratory chain, optimizing respiratory efficiency at different oxygen concentrations. As a result, mammalian cells respond to hypoxic environments by altering the composition of the COX subunits . Levels of the protein-coding genes COX4I1 and COX4I2 are regulated by oxygen concentration: under normal concentrations, COX4I1 levels rise, while those of COX4I2 decrease. In hypoxic conditions, hypoxia-inducible factor (HIF) induces the expression of COX4I2, thereby raising COX4I2 levels, altering the membrane potential of the mitochondrion and elevating levels of ATP and reactive oxygen species (ROS) . ROS levels have been linked to angiogenesis, including activities such as stimulating growth and migration in endothelial cells and modulating the levels of vascular endothelial growth factor (VEGF), its receptor, nuclear factor kappa-B (NF-κB), mitogen-activated protein kinase (MAPK), and matrix metalloproteinases (MMPs) [15,16,17].
It is well-known that hypoxia is related to tumor progression and poor prognosis. However, the action of COX4I2, which is upregulated under hypoxic conditions, in tumor progression is poorly understood. Here, we demonstrated a novel molecular pathway whereby activation of COX4I2 significantly induces fibroblast growth factor 1 (FGF1) expression and thus promotes tumor-associated angiogenesis and the activation of cancer-associated fibroblasts (CAFs).
Materials and methods
The workflow and mechanism diagram of this study is provided in the Pictorial Abstract.
Reagents and antibodies
Information on the reagents and antibodies is given in Additional file 1: Table S1. The concentrations used were based on previous studies, as recommended in the protocols provided. Further details are given in the Additional file.
The human CRC cell lines SW620 (Cat. no. TCHu101), SW480 (Cat. no. SCSP-5033), RKO (Cat. no. TCHu117), HCT116 (Cat. no CL-0096), HT-29 (Cat. no. TCHu103), and LoVo (Cat. no. SCSP-514) were acquired from the cell bank of the Chinese Academy of Sciences (Shanghai, China). The human normal colon epithelial cell (NCM460) (Cat no. CP-H040), human mesenchymal stromal cells (MSCs) (Cat. no. CP-H166), and vascular endothelial cells (human umbilical vein endothelial cells [HUVECs]; Cat. no. CL-0122) were obtained from Cyagen Biosciences (Guangzhou, China). SW620 cells were grown in Leibovitz’s L-15 medium with 10% fetal bovine serum (FBS), while the other lines were grown in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% FBS under 5% CO2 at 37 °C. All media were supplemented with 1% penicillin/streptomycin. Cells passed regular mycoplasma contamination testing. All cells were recently authenticated by short tandem repeat (STR) profiling.
Ethics and sample collection
Thirty paired CRC tissues and adjoining normal tissues (margin: 5 cm) were obtained intra-operatively from untreated patients at the Jiangsu Province Hospital of Chinese Medicine. All experiments were conducted in accordance with the Declaration of Helsinki and in accordance with the relevant designated guidelines and regulations . Tumors were scored using the 8th edition of the American Joint Committee on Cancer tumor-node-metastasis (TNM) system . The samples were rinsed with cold phosphate buffered saline, quick-frozen in liquid nitrogen, and were stored at − 80 °C until analysis. Serum samples obtained preoperatively were stored at – 80 °C.
Immunohistochemistry (IHC) was performed according to a previously published protocol . Micrographs were obtained using a Nikon Eclipse Ni-E microscope (Nikon Corp., Tokyo, Japan) (original magnification: ×400). The H-SCORE (range: 0–300, with higher scores indicating stronger positive staining) was determined according to the published protocol.
Establishment of hypoxia model
A stock solution of CoCl2 was prepared by dissolving solid CoCl2 in serum-free DMEM. Using a previously published method , we pretreated cells with 300 µmol/L CoCl2 for 30 h. CoCl2 induces hypoxia through blocking the degradation of HIF-α, thereby activating the hypoxia cascade .
The western blotting method is described in . The protein bands (including β-actin as the loading control) were visualized with a gel imaging system (ChemiDoc XRS+; Bio-Rad Laboratories, Hercules, CA, USA), and relative concentrations were determined.
RNA interference and overexpression of plasmid construction and transfection
SW480 and RKO cells, both of which express relatively high levels of COX4I2, were selected for subsequent experiments. The plasmids described here were all constructed by the commercial hi-tech company GeneChem (Shanghai, China). The procedures for constructing overexpression and RNA interference (RNAi) plasmids are described in Additional file. Three short hairpin interfering RNAs (shRNAs) targeting COX4I2 were constructed, and the one with the highest inhibitory efficiency was selected for subsequent experiments. The shRNA-COX4I2 plasmid (sh-COX4I2), a control non-targeting plasmid (NC), and the overexpressing plasmid (oe-COX4I2) were transfected into 70%-confluent cells using Lipofectamine 3000, in accordance with the described protocol. Overexpression, knockdown, and transfection efficiency were evaluated using western blotting and green fluorescent protein (GFP) expression.
Colony-forming assays were conducted according to published protocols . Colony numbers were determined under a compound light microscope (model BX53; Olympus Corp., Tokyo, Japan).
Cell Counting Kit-8 assay
Viability was assessed with a Cell Counting Kit-8 (CCK8) assay kit following published methodology . Cells were grown in 96-well plates (5 × 103 cells/well) in 100 μL serum-free medium for 12 or 24 h at 37 °C. CCK-8 solution (100 μL/well) was added for 2 h, following which the cells were grown for 12 or 24 h and absorbance read at 450 nm using a microplate reader (Synergy HT; Biotek Instruments, Winooski, VT, USA ).
For COX4I2 overexpression rescue experiments, the COX4I2-overexpressing cell lines were pretreated with 5 μM PD166866 for 12 h. PD166866 is a synthetic small molecule that inhibits the FGF-1 receptor tyrosine kinase (FGFR1); it is highly selective for FGFR1 with potential use as anti-proliferative/antiangiogenic agent for such therapeutic targets as tumor growth and neovascularization [26, 27]. Optimal concentration selection of PD166866 was evaluated with CCK8 kits (Additional file 1: Section：Concentration screening for PD-166866), with concentrations determined with reference to published studies.
Wound healing assay
The wound healing assay was conducted according to a published protocol . An inverted fluorescence microscope (CKX-41; Olympus Corp.) was used to image migration. GFP levels in the cells on the lower surface of the membrane were expressed as the mean ± standard error of the mean (SEM) from three separate experiments.
Cell invasion was examined using a Transwell assay as previously described . HUVECs or CAFs were seeded in 500 μL serum-free DMEM in the lower chamber of the 24-well Transwell apparatuses. HUVECs (6 × 103/well) or CAFs (5 × 103/well) were suspended in 200 μl 10% FBS medium and seeded in the upper chamber. After co-culturing for 24 h, the migrated cells were fixed in methanol and stained with 0.1% crystal violet for 15 min. The cells on the membrane were imaged under a light microscope (model BX53; Olympus Corp.) and counted with Image J software.
Immunofluorescence staining followed the method described in our previous study . Staining was examined by epifluorescence microscopy (model BX60-32FB2-A03; Olympus Corp.), and different filters were used for capturing images using an Olympus DP50 camera.
Enzyme-linked immunosorbent assay
Proteins in the culture supernatants or mouse/patient sera were measured using enzyme-linked immunorobent assay ( (ELISA)) kits, as per the manufacturer’s instructions. Absorbances (450 nm) were determined using a microplate reader (Synergy HT; Biotek Instruments) (Additional File 1; Table S1).
Endothelial cell tube formation assay
Tube formation assays are useful for assessing the angiogenic capacity of cells . In the present study we used the branch point method. In brief, CRC cells were grown in the upper chambers of 24-well Transwell apparatuses (pore size: 0.4 µm) and HUVEC cells (5 × 104/well) were grown in the lower chambers. The plates were pre-coated with Matrigel (50 µL) and cultured at 37 °C in a 5% CO2 incubator. After 12 h, the HUVECs were incubated with 1 mM calcein-acetoxymethyl ester (calcein-AM) for 1 h at 37 °C and examined under a fluorescence microscope (model CKX-41; Olympus Corp.) (Additional File 1; Table S1).
Establishment of a CRC cell/MSC co-culture unit
A non-contact MSC and CRC cell co-culture system was constructed using a co-culture Transwell system (pore size: 0.4 μm) (upper chamber, CRC cells; lower chamber, MSCs) in which the culture medium was diffusible but the cells could not permeate. The medium was replaced at 48-h intervals. After 5 days, the cells in the lower chamber were collected .
Xenograft tumor model
All of the animal experiments were performed according to the “Guide for the Care and Use of Laboratory Animals” (NIH Publication No. 80-23, revised 1996; National Academies Press, Washington DC, USA) and approved by the ethics committee of the Jiangsu Province Hospital of Chinese Medicine (no. 2021-10-029). Twenty-four 4-week-old male BALB/c nude mice were obtained from the Beijing Institute of Biomedicine (Beijing, China) (Certificate No. SYXK2019-0010). RKO cells transfected with sh-COX4I2, oe-COX4I2, as well as NC cells (1 × 107 cell/mouse) were injected subcutaneously into the right armpit region (n = 6 per group). The presence of tumors was visible after 7 days. Tumor diameters (maximum and minimum) were measured twice a week. The mice were euthanized under CO2 on day 35, and the sera and tumors were excised. Euthanasia was fully consistent with the recommendations of the Guidelines on Euthanasia of the American Veterinary Medical Association . Tumor volumes were calculated as V = 1/2ab2, and tumor growth curves were drawn.
Transcriptomic expression and survival analysis
The differential expression of COX4I2 and COX19 (cytochrome c oxidase assembly factor COX19)in CRC tissue was investigated using The Cancer Genome Atlas (TCGA)-Colon Adenocarcinoma (COAD) cohort and the GSE10950  (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE10950) and GSE37182  (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE37182) datasets. Expression levels were analyzed by T (Tumor) stage, N (Node) stage, and M (Metastasis) according to the TCGA-COAD data (https://portal.gdc.cancer.gov/projects/TCGA-COAD) . To investigate the effects of COX4I2 and COX19 on patient outcomes, survival analyses, including overall survival (OS), disease-free survival (DFS), and progression-free interval (PFI) were conducted using log-rank tests on the TCGA-COAD data.
Least absolute shrinkage and selection operator and COX regression
Raw RNA-sequencing (RNA-Seq) data and clinical information on COX-related genes were downloaded from the TCGA dataset (https://portal.gdc.cancer.gov/) in January 2020. Least absolute shrinkage and selection operator (LASSO) analysis was performed with the R package ‘glmnet’ (R Foundation for Statistical Computing, Vienna, Austria). Univariate and multivariate Cox regression was conducted to identify the prognostic risk of COX19 and COX4I2. Forest plots were determined using the “forestplot” R package.
COX4I2-correlated gene enrichment analysis
High and low COX4I2 expression groups were defined by the median cutoff in the TCGA-COAD dataset. The “DESeq” R package  was used to identify COX4I2-associated genes. The “clusterProfiler” R package  was used for functional enrichment analysis.
Gene set enrichment analysis (GSEA) was performed using the Broad Institute GSEA software 3.0 . The gene set “subset of Gene Ontology (GO)” was acquired from the Molecular Signatures Databases (http://www.gsea-msigdb.org/gsea/msigdb/index.jsp) and was used for GO enrichment analysis. The false discovery rate (FDR) < 0.1 was considered to be statistically significant.
We downloaded a single-cell RNA-Seq dataset named CRC-016-01-1A (https://ngdc.cncb.ac.cn/cancerscem/browse) from the Database of Single-cell Expression Map (CancerSCEM), which consists of one CRC patient’s cells. All samples were imported into Seurat V3 and visualized in Uniform Manifold Approximation and Projection (UMAP) after a standardized quality control process. Next, the “FindAllMarkers” function (min.pct = 0.25, logfc.threshold = 2, tes.use = “wilcox”) was set to find markers for all cell types, and the “AddModuleScore” function was conducted to visualize gene set enrichment results. The “Dotplot” function was used to plot dot plots, and the “Vlnplot” function was used to plot violin plots.
The Tumor Immune Estimation Resource (TIMER) 2.0 web tool (http://timer.cistrome.org) was used to evaluate the relationship between COX4I2 and fibroblasts . Spearman’s rank correlation coefficient was calculated for pairwise correlation comparisons and P < 0.05 was considered to be statistically significant.
Data are expressed as the mean ± SEM. Comparisons between two groups and multiple groups were assessed by t-tests and one-way analysis of variance, respectively. All data were analyzed with SPSS version 29.0 software (SPSS IBM Corp., Armonk, NY, USA) and illustrated using GraphPad Prism 8.0 (GraphPad Software, Inc., San Diego, CA, USA). All experiments were carried out at least three times. Statistical significance was determined as: ***P < 0.001, **P < 0.01 and *P < 0.05.
Identification of the research targets
A total of 42 COX-related genes were used for LASSO regression to identify robust markers for follow-up studies. Cross-validation was performed in 10 rounds to prevent overfitting (Fig. 1a, b). The results of sequential univariate and multivariate Cox regression analyses showed that both COX4I2 and COX19 were associated with CRC prognosis, identifying these genes as research targets (Fig. 1c–h).
Expression and survival analysis
The expression of COX19 and COX4I2 in various tumors was determined using the TIMER web tool. As shown in Fig. 2a, b, both of these genes were highly expressed in COAD. To further explore their expression in CRC tissues, data from the Gene Expression Omnibus (GEO) database were analyzed. As shown in Fig. 2c–f, the data extracted from GSE37182 and GSE10950 indicated that there was a significant difference between COX19 and COX4I2 expression in tumor tissues and controls (P < 0.05). The expression of the genes was further analyzed through TCGA, and the results showed that expression levels of both COX19 and COX4I2 were increased in COAD tissues compared with controls (Fig. 2g) (P < 0.01).
The association between COX19 and COX4I2 levels and clinicopathological characteristics is presented in Fig. 2h (factors including T/N/M stage). The mRNA expressions of COX19 correlate with the N1, N2, and M1 stage, while COX4I2 expression correlates with the N2 stage (P < 0.05).
Results from the Gene Expression Profiling Interactive Analysis (GEPIA) database showed both lower overall survival (OS; log-rank test, P = 0.011) and disease-free survival (DFS; log-rank test, P = 0.026) in the high COX4I2 expression group, while high COX19 levels were only associated with OS (log = rank test, P = 0.04) (Fig.2I, J). In addition, Kaplan–Meier (K–M) survival curves based on TCGA-COAD confirmed the relationship between elevated COX4I2 and OS, DFS, and progression-free interval (PFI) (Fig. 2k, l).
Functional enrichment analysis
A functional network based on functional predictions for COX19 and COX4I2 was constructed with GeneMANIA (Fig. 3a). The differentially expressed genes (DEGs) are shown in Fig. 3b, c. Enrichment analysis of the DEGs using METASCAPE suggested that COX4I2 may be more involved in fibrosis-related pathways, such as elastic fiber, extracellular matrix, integrin, and collagen, compared to COX19. Consequently, COX4I2 was selected for further study. The protein–protein interaction (PPI) network for COX4I2 is shown in Fig. 3e. GO and KEGG analyses indicated the involvement of COX4I2 in angiogenesis and stromal remodeling (Fig. 3f, g). The key subnetwork (Fig. 3h) obtained by the molecular complex detection (MCODE) plug-in showed a stronge association with angiogenesis (Fig. 3i). GSEA was then applied to further identify the possible mechanisms involved, and we observed that FGF may be closely associated with COX4I2 (Fig. 3j). Finally, we calculated the Pearson correlation coefficient between COX4I2 and all FGFs (Fig. 3k) at the transcriptomic level based on the TCGA-COAD data (Fig. 3l), from which we observed a significant correlation between FGF1 and COX4I2, with the results from the CRC patients’ sera supporting this (Fig. 3m). Finally, we found a high positive correlation between COX4I2 level and fibroblast abundance by MCPCOUNTER (R = 0.530, P = 4.87e−22; Fig. 3n) and EPIC methods (R = 0.508, P = 2.2e−19; Fig. 3o). These results implied that CXO4I2 was closely associated with tumor fibrosis in CRC.
Single-cell profiling enables a better understanding of the composition of the CRC microenvironment. We investigated the association between COX4I2 and FGFs using single-cell datasets. We identified additional genes that strongly and specifically marked each major cell population, as demonstrated in Fig. 4a. Figure 4b shows the cell markers of each cell population using UMAP. As seen in Fig. 4c, d, we observed, using UMAP and violin plots with unique markers for each cell type, that COX4I2 was mainly expressed in cluster 7. As shown in Fig. 4e–g, functional enrichment analysis suggested that cells expressing COX4I2 tended to express genes related to “HYPOXIA,” “EPITHELIAL-INTERSTITIAL-TRANSFORMATION,” and “ANGIOGENESIS”. The violin plot (Fig. 4h–j) shows that scores of “HYPOXIA,” “EPITHELIAL-INTERSTITIAL-TRANSFORMATION,” and “ANGIOGENESIS” were significantly higher in cluster 7 than in other clusters. Additional file 1: Fig.S1 also confirms an association between COX4I2 and hypoxia-related gene expression and may contribute to an understanding of the potential mechanisms of COX4I2.
COX4I2 promotes malignant CRC phenotypes
COX4I2 expression in CRC tissue was further explored using western blotting and IHC, which indicated high COX4I2 levels in both CRC cells and tissues. The mean (± SEM) H-SCOREs for COX4I2 expression in CRC and the paired control tissues were 95.53 ± 15.96 and 19.14 ± 6.65, respectively (Fig. 5a, b). COX4I2 levels were highest in the RKO and SW480 cell lines and were induced and activated by hypoxia in a time-dependent manner, consistent with HIF-1α (Fig. 5c). The efficiency of transfection was confirmed by GFP expression and western blotting (Fig. 5d). Silencing of COX4I2 reduced clone formation (Fig. 5e). Moreover, stable COX4I2 overexpression resulted in successful formation of xenograft tumors in mice (Fig. 5f–h). Ki-67 is a known prognosis-related factor for CRC, with strong expression linked to aggressive tumors and poor outcomes, and it is also a proliferation marker [41, 42]. We observed that CRC cases and tumors of nude mice with high COX4I2 expression also exhibited strong Ki-67 staining signals (Fig. 5i, j). These results indicated that COX4I2 may be correlated with Ki-67.
Based on the results of the enrichment analysis and observed correlation between COX4I2 and FGF1 in CRC patient sera, we speculated that FGF1 may be involved in some of the cancer-associated behaviors regulated by COX4I2. Subsequently, we further verified the correlation of COX4I2 with FGF1 using culture supernatants (Fig. 5k), mouse sera (Fig. 5l), mouse tumors (Fig. 5m), and CRC tissues (Fig. 5n, o). COX4I2 expression was found to be strongly correlated with FGF1 expression in all cases.
COX4I2 overexpression promoted the epithelial–mesenchymal transition process
GSEA and single-cell RNA-Seq analyses showed that COX4I2 expression was closely correlated with the epithelial–mesenchymal transition (EMT) process. Increases in EMT-related phenotypes were found to be linked to COX4I2 overexpression. This influence decreased with treatment with the integrin FGF1-specific inhibitor PD-166866 (Fig. 6a, b). Addition of FGF1 (10 ng/mL) rescued COX4I2 knockdown CRC cells EMT phenotype (Additional file 1: Fig. S2). The cell viability assays showed that cell proliferation was not influenced at the chosen PD-166866 concentration (Additional file 1: Fig.S3; Additional file 1: Fig.S4). Correlation analysis (Fig. 6c) showed that COX4I2 expression was positively associated with CDH2 (R = 0.320, P < 0.001), MMP2 (R = 0.450, P < 0.001), MMP9 (R = 0.400, P < 0.001), SNAI2 (R = 0.380, P < 0.001), and SNAI1 (R = 0.450, P < 0.001) expression and negatively linked with CDH1 expression (R = − 0.070, P = 0.125), in agreement with the in vitro observations. These results were confirmed in vivo at the protein level (Fig. 6d).
COX4I2 promotes angiogenesis and CAF activation
The TFacts tool (http://www.tfacts.org/) was employed to predict transcription factors that regulate COX4I2-related gene expression in CRC, using all genes that were found to be positively correlated with COX4I2 in the COAD dataset of TCGA. This analysis showed that transcription factors involved in angiogenesis and hypoxia-related signaling pathways were significantly enriched, such as SMAD2/7 and NF-κB1. The enrichment analysis also suggested that these transcription factors are closely associated with angiogenesis and CRC (Additional file 1: Fig.S5).
Considering these results, and also to verify the GSEA and single-cell RNA-Seq analyses of angiogenesis and fibroblast activation, we established several in vitro co-culture models. We first noted that overexpression of COX4I2 in CRC promoted the migration of HUVECs (Fig. 7a–c) as well as their ability to form blood vessels (Fig. 7d–f). This malignant phenotype could be rescued by inhibitors of FGF1. CAFs, which are known for their properties in forming the tumor microenvironment, originate from circulating MSCs. The non-contact co-culture model (Fig. 7g) showed that COX4I2 overexpression could promote the differentiation of MSCs to CAFs, which was mainly reflected by a significant increase in markers of CAF activation in the co-culture system (Fig. 7h–l). Fibroblast activation protein (FAP), fibroblast-specific protein-1 (S100A4), and α-smooth muscle actin (α-SMA) are all commonly used to label activated CAFs. As we have previously demonstrated, these observations also showed that pretreatment with an FGF1 inhibitor reduced MSC differentiation into CAFs in the co-culture unit. The invasion of cancer cells not only refers to the tumor cells themselves, but also encompasses the recruitment of endothelial cells and fibroblasts in the surrounding stromal microenvironment. Migration assays confirmed that COX4I2 overexpression increased CAF recruitment by the cancer cells (Fig. 7m–o). Additional addition of FGF1 (10 ng/mL) rescued the phenotype after knockdown of COX4I2, including angiogenesis and fibroblast activation (Additional file 1: Fig.S6).
Together, these findings indicated that COX4I2 promotes CRC progression by stimulating angiogenesis and fibroblast activation.
Colorectal cancer is a highly invasive and heterogeneous disease . There are several aspects to its heterogeneity. On the one hand, the disordered proliferation of the cancer cell leads to increased oxygen consumption and hypoxia in the local microenvironment, stimulating angiogenesis and the formation of abnormal tumor-associated blood vessels . On the other hand, tumor-associated stromal cells are exposed to the same hypoxic microenvironment as the tumor parenchymal cells and contribute to tumor progression . As one of the most striking features of the tumor microenvironment, hypoxia represents an adaptive mechanism for tumor development that drives tumor progression and adversely affects treatment outcome [46, 47]. There are currently no methods to quantify or determine the degree of tumor hypoxia. Therefore, the development of hypoxia-related biochemical markers could help improve our understanding of the biological mechanisms of CRC and the management of CRC therapy.
Mammals rely heavily on aerobic respiration for their energy needs, with most of the inhaled oxygen being used by COX  which, as the last enzyme in the respiratory chain, must cope with fluctuating oxygen concentrations . Similar response mechanisms also exist for cancer cells, and these are closely linked to the malignant properties of cancer cells, of which one of the best known aspects is the ability of increased ROS production through the mitochondrial respiratory respiratory chain to strongly drive cancer cell invasion and metastasis [50, 51]. This activity indicates the likely involvement of respiratory chain complexes in remodeling the hypoxic tumor microenvironment. We therefore performed a LASSO-COX regression based on TCGA-COAD data and identified COX4I2 as well as COX19 in respiratory chain complexes as potential cancer risk factors. There is little documentation of these two genes in relation to tumors, and the topic is worthy of investigation. We then compared the potential biological functions of both genes by functional enrichment analysis, which indicated the likely involvement of COX4I2 with multiple fibrosis-associated pathways, suggesting that this gene would be a noval target of further investigation. The induction of COX4I2 expression under hypoxic conditions has been reported in numerous cell lines, and it has been proposed that the COX4 isoform is specifically linked to stimulating mitochondrial respiration in response to reduced oxygen. Our findings showed that COX4I2 is a significant risk factor in CRC, and we speculate that it may play a similar role in oxygen adaptation in CRC.
The analysis of survival data showed that elevated COX4I2 expression was related to unfavorable CRC prognosis and that COX4I2 levels were significantly elevated in CRC tissues at both the transcriptomic and proteomic levels. Ki67 is involved in early rRNA synthesis, has been linked to a variety of signaling pathways, and is induced by hypoxia [52, 53]. We used Ki67 as a positive control, although Ki67 expression does not show a homogeneous pattern in all tissues (related to the spatial heterogeneity and temporal variability of hypoxia in tumors), We found that COX4I2 expression was highly consistent with Ki67 expression and that COX4I2 overexpression significantly stimulated the proliferative capacity of CRC cells. We hypothesize that COX4I2 may be a downstream factor of hypoxia-triggered signaling cascades involved in the onset of cancer proliferation.
Enrichment analysis suggested that COX4I2 may be associated with the activation of FGFs, and evidence at the single-cell level implied a close link between COX4I2 and hypoxia, EMT, and angiogenesis. Calculation of the correlations between COX4I2 and all FGF family members based on TCGA-COAD revealed that the closest link may exist between COX4I2 and FGF1. FGF1 has oncogenic activity and has been found to modulate numerous cellular processes, including proliferation, invasion, and survival . FGF1 levels are elevated in a variety of tumors and the gene is implicated in malignant progression . In tumors, FGF1 may be derived from both tumor and surrounding stromal cells and may have both paracrine or autocrine activity . Secreted FGF1 binds and activates specific FGF receptors (FGFRs) which have tyrosine kinase activity, with intrinsic tyrosine kinase activity on the cell surface initiating signal transduction cascades with various effects, including the promotion of the EMT . It has been shown that FGF1 can promote angiogenesis and fibroblast differentiation [58, 59]. Activated fibroblasts can further secrete FGF1, contributing to a vicious cycle of malignant phenotypes . Evidence from clinical samples as well as in vivo experiments both suggested that COX4I2 acts as a powerful FGF simulator for FGF1, a fact that we also consider being a likely reason for COX4I2’s link to poor prognosis.
It has been well established that FGF1 promotes the metastasis of various tumors through the EMT . The role of the EMT in neovascularization is more complex [62, 63]. Because the process of tumor metastasis is a continuation of the cascade response, it can be envisaged that the end of one cascade will feed into the beginning of the next . For example, SNAIL1-mediated EMT may contribute to neoangiogenesis in the case of metastases that rely on an angiogenic pattern rather than the EMT . In our study, the Transwell and wound-healing assays confirmed the ability of COX4I2 to activate the EMT via the FGF1 pathway, a conclusion reinforced by calculations based on the TCGA-COAD data and in vitro experiments.
FGF1 is a non-classically released growth factor which is essential for vascular endothelial cell proliferation and is an even more potent angiogenic factor than VEGF or platelet-derived growth factors [65, 66]. We clarified that COX4I2 can promote angiogenesis by regulating the level of FGF1. CAFs originate from circulating MSCs  and are involved in the formation of tumor-associated stroma . The differentiation of tumor-infiltrating MSCs to CAFs can be facilitated by FGF1. CAFs can further secrete FGF1 to actively participate in tumor-stromal crosstalk and further drive tumor progression . The results of different algorithms in our study all indicated a significant relationship between COX4I2 levels and the degree of fibroblast infiltration. Also, in the co-culture system we established, we obtained direct evidence that elevated COX4I2 expression enhanced the recruitment and activation of fibroblasts.
This study represents a preliminary exploration of the possible role of COX4I2 in CRC, with the results revealing the complexity of its underlying mechanisms. Although there are still many questions about how COX4I2 affects FGF1, we present evidence that FGF1 is an important molecule in the signaling network in which COX4I2 is involved. In addition, how tumor cells adapt to hypoxia through COX in the tumor microenvironment is an interesting and significant topic, and our subsequent work will focus on these two aspects to try to elucidate the exact pro-cancer mechanism of COX4I2.
Availability of data and materials
Jie-pin Li can be contacted (firstname.lastname@example.org) regarding the availability of data and materials.
Cytochrome c oxidase
- COX19 :
Gene encoding cytochrome c oxidase assembly factor COX19
Cytochrome c oxidase subunit 4 isoform 2 encoded by the COX4I2 gene
Dulbecco’s Modified Eagle Medium
Enzyme-linked immunosorbent assay
Fetal bovine serum
False discovery rate
Fibroblast growth factor 1
Human umbilical vein endothelial cells
Gene expression profiling interactive analysis
Gene set enrichment analysis
Kyoto Encyclopedia of Genes and Genomes
Least absolute shrinkage and selection operator
Mitogen-activated protein kinase
Mesenchymal stromal cells
Nuclear factor kappa-B
Platelet-derived growth factors
Reactive oxygen species
Short hairpin interfering RNA
Short tandem repeat
The Cancer Genome Atlas
Tumor Immune Estimation Resource
Uniform manifold approximation and projection
Vascular endothelial growth factor
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33.
Dekker E, Tanis PJ, Vleugels JLA, Kasi PM, Wallace MB. Colorectal cancer. Lancet. 2019;394(10207):1467–80.
Picard E, Verschoor CP, Ma GW, Pawelec G. Relationships between immune landscapes, genetic subtypes and responses to immunotherapy in colorectal cancer. Front Immunol. 2020;11:369.
Thanikachalam K, Khan G. Colorectal cancer and nutrition. Nutrients. 2019;11(1):164.
Yu H, Hemminki K. Genetic epidemiology of colorectal cancer and associated cancers. Mutagenesis. 2020;35(3):207–19.
Jiang SH, Hu LP, Wang X, Li J, Zhang ZG. Neurotransmitters: emerging targets in cancer. Oncogene. 2020;39(3):503–15.
Catalano V, Turdo A, Di Franco S, Dieli F, Todaro M, Stassi G. Tumor and its microenvironment: a synergistic interplay. Semin Cancer Biol. 2013;23(6):522–32.
Maishi N, Hida K. Tumor endothelial cells accelerate tumor metastasis. Cancer Sci. 2017;108(10):1921–6.
Heinemann V, von Weikersthal LF, Decker T, Kiani A, Kaiser F, Al-Batran SE, et al. FOLFIRI plus cetuximab or bevacizumab for advanced colorectal cancer: final survival and per-protocol analysis of FIRE-3, a randomised clinical trial. Br J Cancer. 2021;124(3):587–94.
Grothey A, Van Cutsem E, Sobrero A, Siena S, Falcone A, Ychou M, et al. Regorafenib monotherapy for previously treated metastatic colorectal cancer (CORRECT): an international, multicentre, randomised, placebo-controlled, phase 3 trial. Lancet. 2013;381(9863):303–12.
Pajuelo Reguera D, Čunátová K, Vrbacký M, Pecinová A, Houštěk J, Mráček T, et al. Cytochrome c oxidase subunit 4 isoform exchange results in modulation of oxygen affinity. Cells. 2020;9(2):443.
Tsukihara T, Aoyama H, Yamashita E, Tomizaki T, Yamaguchi H, Shinzawa-Itoh K, et al. The whole structure of the 13-subunit oxidized cytochrome c oxidase at 2.8 A. Science. 1996;272(5265):1136–44.
Hüttemann M, Kadenbach B, Grossman LI. Mammalian subunit IV isoforms of cytochrome c oxidase. Gene. 2001;267(1):111–23.
Fukuda R, Zhang H, Kim JW, Shimoda L, Dang CV, Semenza GL. HIF-1 regulates cytochrome oxidase subunits to optimize efficiency of respiration in hypoxic cells. Cell. 2007;129(1):111–22.
Prasad S, Gupta SC, Tyagi AK. Reactive oxygen species (ROS) and cancer: role of antioxidative nutraceuticals. Cancer Lett. 2017;387:95–105.
Zou J, Fei Q, Xiao H, Wang H, Liu K, Liu M, et al. VEGF—a promotes angiogenesis after acute myocardial infarction through increasing ROS production and enhancing ER stress-mediated autophagy. J Cell Physiol. 2019;234(10):17690–703.
Yang G, Chang CC, Yang Y, Yuan L, Xu L, Ho CT, et al. Resveratrol alleviates rheumatoid arthritis via reducing ROS and inflammation, inhibiting MAPK signaling pathways, and suppressing angiogenesis. J Agric Food Chem. 2018;66(49):12953–60.
World Medical Association Declaration of Helsinki. Ethical principles for medical research involving human subjects. JAMA. 2013;310(20):2191–4.
Weiser MR. AJCC 8th edition: colorectal cancer. Ann Surg Oncol. 2018;25(6):1454–5.
Paschalis A, Sheehan B, Riisnaes R, Rodrigues DN, Gurel B, Bertan C, et al. Prostate-specific membrane antigen heterogeneity and DNA repair defects in prostate cancer. Eur Urol. 2019;76(4):469–78.
Wu H, Huang S, Chen Z, Liu W, Zhou X, Zhang D. Hypoxia-induced autophagy contributes to the invasion of salivary adenoid cystic carcinoma through the HIF-1α/BNIP3 signaling pathway. Mol Med Rep. 2015;12(5):6467–74.
Rana NK, Singh P, Koch B. CoCl(2) simulated hypoxia induce cell proliferation and alter the expression pattern of hypoxia associated genes involved in angiogenesis and apoptosis. Biol Res. 2019;52(1):12.
Kuriakose S, Onyilagha C, Singh R, Olayinka-Adefemi F, Jia P, Uzonna JE. TLR-2 and MyD88-dependent activation of MAPK and STAT proteins regulates proinflammatory cytokine response and immunity to experimental Trypanosoma congolense infection. Front Immunol. 2019;10:2673.
Peng HH, Wang JN, Xiao LF, Yan M, Chen SP, Wang L, et al. Elevated serum FGG levels prognosticate and promote the disease progression in prostate cancer. Front Genet. 2021;12: 651647.
Fang H, Jin J, Huang D, Yang F, Guan X. PAI-1 induces Src inhibitor resistance via CCL5 in HER2-positive breast cancer cells. Cancer Sci. 2018;109(6):1949–57.
Calandrella N, Risuleo G, Scarsella G, Mustazza C, Castelli M, Galati F, et al. Reduction of cell proliferation induced by PD166866: an inhibitor of the basic fibroblast growth factor. J Exp Clin Cancer Res. 2007;26(3):405–9.
Chen Y, Xie X, Li X, Wang P, Jing Q, Yue J, et al. FGFR antagonist induces protective autophagy in FGFR1-amplified breast cancer cell. Biochem Biophys Res Commun. 2016;474(1):1–7.
Han JH, Yoon JS, Chang DY, Cho KG, Lim J, Kim SS, et al. CXCR4-STAT3 axis plays a role in tumor cell infiltration in an orthotopic mouse glioblastoma model. Mol Cells. 2020;43(6):539–50.
Cao R, Wang G, Qian K, Chen L, Qian G, Xie C, et al. Silencing of HJURP induces dysregulation of cell cycle and ROS metabolism in bladder cancer cells via PPARγ-SIRT1 feedback loop. J Cancer. 2017;8(12):2282–95.
Liu YJ, Zeng SH, Hu YD, Zhang YH, Li JP. Overexpression of NREP promotes migration and invasion in gastric cancer through facilitating epithelial–mesenchymal transition. Front Cell Dev Biol. 2021;9:746194.
Yang Z, Li H, Luo P, Yan D, Yang N, Zhang Y, et al. UNC5B promotes vascular endothelial cell senescence via the ROS-mediated P53 pathway. Oxid Med Cell Longev. 2021;2021:5546711.
Fan H, Atiya HI, Wang Y, Pisanic TR, Wang TH, Shih IM, et al. Epigenomic reprogramming toward mesenchymal–epithelial transition in ovarian-cancer-associated mesenchymal stem cells drives metastasis. Cell Rep. 2020;33(10): 108473.
Kilkenny C, Browne WJ, Cuthill IC, Emerson M, Altman DG. Improving bioscience research reporting: the ARRIVE guidelines for reporting animal research. PLoS Biol. 2010;8(6):e1000412.
Ebadfardzadeh J, Kazemi M, Aghazadeh A, Rezaei M, Shirvaliloo M, Sheervalilou R. Employing bioinformatics analysis to identify hub genes and microRNAs involved in colorectal cancer. Med Oncol. 2021;38(9):114.
Sharma A, Yadav D, Rao P, Sinha S, Goswami D, Rawal RM, et al. Identification of potential therapeutic targets associated with diagnosis and prognosis of colorectal cancer patients based on integrated bioinformatics analysis. Comput Biol Med. 2022;146: 105688.
Shi X, Kaller M, Rokavec M, Kirchner T, Horst D, Hermeking H. Characterization of a p53/miR-34a/CSF1R/STAT3 feedback loop in colorectal cancer. Cell Mol Gastroenterol Hepatol. 2020;10(2):391–418.
Bao Y, Wang L, Shi L, Yun F, Liu X, Chen Y, et al. Transcriptome profiling revealed multiple genes and ECM-receptor interaction pathways that may be associated with breast cancer. Cell Mol Biol Lett. 2019;24:38.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Reimand J, Isserlin R, Voisin V, Kucera M, Tannus-Lopes C, Rostamianfar A, et al. Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, Cytoscape and Enrichment Map. Nat Protoc. 2019;14(2):482–517.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.
Ishida H, Miwa H, Tatsuta M, Masutani S, Imamura H, Shimizu J, et al. Ki-67 and CEA expression as prognostic markers in Dukes’ C colorectal cancer. Cancer Lett. 2004;207(1):109–15.
Menon SS, Guruvayoorappan C, Sakthivel KM, Rasmi RR. Ki-67 protein as a tumour proliferation marker. Clin Chim Acta. 2019;491:39–45.
Sagaert X, Vanstapel A, Verbeek S. Tumor heterogeneity in colorectal cancer: what do we know so far? Pathobiol J Immunopathol Mol Cell Biol. 2018;85(1–2):72–84.
McDonald PC, Chafe SC, Dedhar S. Overcoming hypoxia-mediated tumor progression: combinatorial approaches targeting pH regulation, angiogenesis and immune dysfunction. Front Cell Dev Biol. 2016;4:27.
Semenza GL. Cancer-stromal cell interactions mediated by hypoxia-inducible factors promote angiogenesis, lymphangiogenesis, and metastasis. Oncogene. 2013;32(35):4057–63.
Jing X, Yang F, Shao C, Wei K, Xie M, Shen H, et al. Role of hypoxia in cancer therapy by regulating the tumor microenvironment. Mol Cancer. 2019;18(1):157.
Erin N, Grahovac J, Brozovic A, Efferth T. Tumor microenvironment and epithelial mesenchymal transition as targets to overcome tumor multidrug resistance. Drug Resist Updat. 2020;53: 100715.
Little AG, Lau G, Mathers KE, Leary SC, Moyes CD. Comparative biochemistry of cytochrome c oxidase in animals. Comp Biochem Physiol B Biochem Mol Biol. 2018;224:170–84.
Yoshikawa S, Shimada A. Reaction mechanism of cytochrome c oxidase. Chem Rev. 2015;115(4):1936–89.
Vogt S, Rhiel A, Weber P, Ramzan R. Revisiting Kadenbach: electron flux rate through cytochrome c-oxidase determines the ATP-inhibitory effect and subsequent production of ROS. BioEssays. 2016;38(6):556–67.
Chen ZX, Pervaiz S. Involvement of cytochrome c oxidase subunits Va and Vb in the regulation of cancer cell metabolism by Bcl-2. Cell Death Differ. 2010;17(3):408–20.
Huang W, Zhang Q, Wu G, Chen PP, Li J, McCabe Gillen K, et al. DCE-MRI quantitative transport mapping for noninvasively detecting hypoxia inducible factor-1α, epidermal growth factor receptor overexpression, and Ki-67 in nasopharyngeal carcinoma patients. Radiother Oncol. 2021;164:146–54.
Riffle S, Pandey RN, Albert M, Hegde RS. Linking hypoxia, DNA damage and proliferation in multicellular tumor spheroids. BMC Cancer. 2017;17(1):338.
Manousakidi S, Guillaume A, Pirou C, Bouleau S, Mignotte B, Renaud F, et al. FGF1 induces resistance to chemotherapy in ovarian granulosa tumor cells through regulation of p53 mitochondrial localization. Oncogenesis. 2018;7(2):18.
Bai YP, Shang K, Chen H, Ding F, Wang Z, Liang C, et al. FGF-1/-3/FGFR4 signaling in cancer-associated fibroblasts promotes tumor progression in colon cancer through Erk and MMP-7. Cancer Sci. 2015;106(10):1278–87.
Sun Y, Fan X, Zhang Q, Shi X, Xu G, Zou C. Cancer-associated fibroblasts secrete FGF-1 to promote ovarian proliferation, migration, and invasion through the activation of FGF-1/FGFR4 signaling. Tumour Biol. 2017;39(7):1010428317712592.
Sun J, Stathopoulos A. FGF controls epithelial-mesenchymal transitions during gastrulation by regulating cell division and apicobasal polarity. Development. 2018;145(19):dev161927.
Matsuki M, Adachi Y, Ozawa Y, Kimura T, Hoshi T, Okamoto K, et al. Targeting of tumor growth and angiogenesis underlies the enhanced antitumor activity of lenvatinib in combination with everolimus. Cancer Sci. 2017;108(4):763–71.
Takahashi M, Okubo N, Chosa N, Takahashi N, Ibi M, Kamo M, et al. Fibroblast growth factor-1-induced ERK1/2 signaling reciprocally regulates proliferation and smooth muscle cell differentiation of ligament-derived endothelial progenitor cell-like cells. Int J Mol Med. 2012;29(3):357–64.
Erdogan B, Webb DJ. Cancer-associated fibroblasts modulate growth factor signaling and extracellular matrix remodeling to regulate tumor metastasis. Biochem Soc Trans. 2017;45(1):229–36.
Jiao J, Zhao X, Liang Y, Tang D, Pan C. FGF1-FGFR1 axis promotes tongue squamous cell carcinoma (TSCC) metastasis through epithelial-mesenchymal transition (EMT). Biochem Biophys Res Commun. 2015;466(3):327–32.
Maroufi NF, Amiri M, Dizaji BF, Vahedian V, Akbarzadeh M, Roshanravan N, et al. Inhibitory effect of melatonin on hypoxia-induced vasculogenic mimicry via suppressing epithelial-mesenchymal transition (EMT) in breast cancer stem cells. Eur J Pharmacol. 2020;881:173282.
Yang F, Shao C, Wei K, Jing X, Qin Z, Shi Y, et al. miR-942 promotes tumor migration, invasion, and angiogenesis by regulating EMT via BARX2 in non-small-cell lung cancer. J Cell Physiol. 2019;234(12):23596–607.
Cabrerizo-Granados D, Peña R, Palacios L, Carrillo-Bosch L, Lloreta-Trull J, Comerma L, et al. Snail1 expression in endothelial cells controls growth, angiogenesis and differentiation of breast tumors. Theranostics. 2021;11(16):7671–84.
Boilly B, Vercoutter-Edouart AS, Hondermarck H, Nurcombe V, Le Bourhis X. FGF signals for cell proliferation and migration through different pathways. Cytokine Growth Factor Rev. 2000;11(4):295–302.
Minina E, Kreschel C, Naski MC, Ornitz DM, Vortkamp A. Interaction of FGF, Ihh/Pthlh, and BMP signaling integrates chondrocyte proliferation and hypertrophic differentiation. Dev Cell. 2002;3(3):439–49.
Timaner M, Tsai KK, Shaked Y. The multifaceted role of mesenchymal stem cells in cancer. Semin Cancer Biol. 2020;60:225–37.
Quante M, Tu SP, Tomita H, Gonda T, Wang SS, Takashi S, et al. Bone marrow-derived myofibroblasts contribute to the mesenchymal stem cell niche and promote tumor growth. Cancer Cell. 2011;19(2):257–72.
Bhattacharyya S, Oon C, Kothari A, Horton W, Link J, Sears RC, et al. Acidic fibroblast growth factor underlies microenvironmental regulation of MYC in pancreatic cancer. J Exp Med. 2020;217(8):e20191805.
The authors gratefully acknowledge contributions from the GEO and TCGA databases.
The present study was supported by the Youth Science and Technology Project of Suzhou (No. KJXW2019059); the Suzhou Science and Technology Development Plan (No. SYSD2019006); Zhangjiagang Health Youth Science and Technology Project (No. ZJGQNKJ202036); Peak Academic Talent Project of Jiangsu Provincial Hospital of Traditional Chinese Medicine (y2021rc19); Advantageous Disciplines Program of Nanjing University of Chinese Medicine (ZYX03KF022 and ZYX03KF019); Science and Technology Project of Affiliated Hospital of Nanjing University of Chinese Medicine (Y2020CX62); State Administration of Chinese Medicine Project (20085-9-3); Jiangsu Provincial Science and Technology Department Project (BE2019771); Jiangsu Province Postgraduate Research Innovation Program Project (KYCX21-1677).
Ethics approval and consent to participate
This study was conducted according to the ethical standards contained in the Declaration of Helsinki, and in national and international guidelines. The study’s protocol was approved by the ethics committee of the Jiangsu Province Hospital of Chinese Medicine, and informed consent was obtained from clinicians and patients (2020NL-107–01). Written information was provided and written consent was given to all participants prior to collection of specimens for research purposes. All mouse experiments were performed in accordance with the Guide for Care and Use of Laboratory Animals and approved by the Animal Ethics Committee of Jiangsu Province Hospital of Chinese Medicine (2021-10-029).
Adherence to ARRIVE guidelines
Our study complies with the ARRIVE guidelines (https://arriveguidelines.org).
Consent for publication
All authors have given consent for publication.
The authors have declared no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, Jp., Liu, Yj., Zeng, Sh. et al. Identification of COX4I2 as a hypoxia-associated gene acting through FGF1 to promote EMT and angiogenesis in CRC. Cell Mol Biol Lett 27, 76 (2022). https://doi.org/10.1186/s11658-022-00380-2
- Colorectal cancer
- Fibroblast growth factor 1
- Epithelial–mesenchymal transition
- Cancer-associated fibroblasts