Skip to main content

High MICAL-L2 expression and its role in the prognosis of colon adenocarcinoma

Abstract

Background

MICAL-like protein 2 (MICAL-L2), a member of the molecules interacting with CasL (MICAL) family of proteins, is strongly associated with the malignancy of multiple types of cancer. However, the role of MICAL-L2 in colon adenocarcinoma (COAD) has not been well characterized.

Methods

In this study, we analyzed the role of MICAL-L2 in COAD using datasets available from public databases. The mRNA and protein expression of MICAL-L2 was investigated using TCGA, UALCAN, and independent immunohistochemical assays. Overall survival (OS) and disease-specific survival (DSS) of COAD patients were assessed based on the MICAL-L2 expression level using the Kaplan–Meier method. Univariate and multivariate analysis was employed to determine whether MICAL-L2 could serve as an independent prognostic indicator of OS. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA) were further utilized to explore the possible cellular mechanism underlying the role of MICAL-L2 in COAD. In addition, the correlation between MICAL-L2 expression and immune cell infiltration levels was investigated via single-sample gene set enrichment analysis (ssGSEA).

Results

Data from TCGA, HPA, and UALCAN datasets indicated that MICAL-L2 expression was significantly higher in COAD tissue than in adjacent normal tissues, and this was confirmed by immunohistochemical assays. Kaplan–Meier survival analysis revealed that patients with MICAL-L2 had shorter OS and DSS. Furthermore, multivariate Cox analysis indicated that MICAL-L2 was an independent risk factor for OS in COAD patients. ROC analysis confirmed the diagnostic value of MICAL-L2, and a prognostic nomogram involving age, M stage, and MICAL-L2 expression was constructed for OS. Functional enrichment analyses revealed that transport-related activity was closely associated with the role of MICAL-L2 in COAD. Regarding immune infiltration levels, MICAL-L2 was found to be positively associated with CD56bright NK cells.

Conclusions

Our results suggested that MICAL-L2 is a promising biomarker for determining prognosis and correlated with immune infiltration levels in COAD.

Peer Review reports

Background

Colon cancer is a commonly diagnosed malignant tumor of the digestive tract and a leading cause of cancer-related death worldwide [1]. It usually affects adults at 40–50 years of age and occurs more often in males than females. The etiology of colon cancer is mainly associated with a high-fat diet, colonic polyps, genetic make-up, and chronic inflammation [2]. Patients in the early stage of colon cancer may not present obvious clinical symptoms. As the tumor volume increases, patients may display abdominal distension and dyspepsia, and may even be able to feel a lump/mass in the abdomen. Although the 5-year survival rate of patients in the early stages of colon cancer can be higher than 90%, that of patients diagnosed at an advanced stage is lower than 20% [3,4,5]. These observations underline the need to further unravel the mechanisms underlying colon cancer progression and identify novel therapeutic targets for the treatment of this disease.

Molecules interacting with CasL (MICALs) represent an evolutionarily conserved family of proteins with roles in the regulation of cytoske leton dynamics [6]. MICAL-like protein 2 (MICAL-L2), a member of the MICAL family, has three conserved domains, namely, a calponin homology (CH) domain; a Lin11, Isl-1, and Mec-3 (LIM) domain; and a C-terminal coiled-coil (CC) domain [7]. The CH and LIM domains link MICAL-L2 to the actin cytoskeleton, while the CC domain is required for interaction with Rab GTPases. MICAL-L2 exerts its multiple biological functions primarily via processes involving cargo transportation. For example, by binding to Rab13, MICAL-L2 triggers the transportation of glucose transporter-4 (GLUT4) and mediates GLUT4-containing vesicle localization and fusion with the muscle cell membrane [8]. Rab13 and MICAL-L2 also act together in the transfer of actinin-4 from the cell body to the tips of neurites [9]. It has been well documented that MICAL-L2 is highly expressed and promotes cell migration and invasion in multiple types of cancer, including gastric cancer, ovarian cancer, and breast cancer [10,11,12]. In ovarian cancer cells, the silencing of MICAL-L2 was shown to inhibit canonical Wnt/β-catenin signaling and induce mesenchymal–epithelial transition [11]. We have previously shown that MICAL-L2 facilitates the proliferation of lung cancer cells via the de-ubiquitination of c-Myc, which blocks its degradation [13]. Recently, another MICAL family member, MICAL1, which shares sequence similarity with MICAL-L2 [14], was found to play a key role in the migration and growth of colorectal cancer cells by suppressing the ERG1/β-catenin signaling pathway [15]. However, the role of MICAL-L2 in the prognosis and possible pathogenesis of colon cancer has not been fully elucidated.

Colon adenocarcinoma (COAD) is one the most common type of colon cancer. In this study, several informatics tools were used to evaluate the expression profile and the prognostic significance of MICAL-L2 in COAD. Moreover, the correlation between MICAL-L2 expression and immune infiltration, and the putative mechanisms underlying the role of MICAL-L2 in COAD, were also investigated. This is the first comprehensive study of the association between MICAL-L2 expression and its clinical characteristics in COAD and our findings may contribute to our understanding of MICAL-L2-related processes in this cancer.

Methods

Ethics statement

All immunohistochemical assays with human tumor specimens were conducted according to the institutional guidelines of Jiangsu Province.

MICAL-L2 mRNA expression and analysis of prognosis

The mRNA expression of MICAL-L2 in COAD and the corresponding clinical information data were downloaded from The Cancer Genome Atlas (TCGA) database (https://tcga-data.nci.nih.gov/tcga/) [16]. MICAL-L2 mRNA expression and its association with overall survival (OS) and disease-specific survival (DSS) of patients with COAD were also analyzed using the TCGA–COAD dataset. The expression of MICAL-L2 was assessed in 456 COAD and 41 adjacent normal tissue samples from the TCGA database. According to the median values of mRNA expression, patients with COAD were divided into high and low expression groups. Data were collected and analyzed using R3.6.3 software [17].

MICAL-L2 protein expression analysis

The Human Protein Atlas (HPA) database (https://www.proteinatlas.org) and the University of Alabama Cancer Database (UALCAN) (http://ualcan.path.uab.edu/index.html) were used to compare MICAL-L2 protein expression between normal and COAD tissues.

Immunohistochemistry

Immunohistochemistry was performed as previously described [18]. COAD tissue microarrays were purchased from Outdo Biotech (Shanghai, China). Thirty paired COAD and paracancerous tissue samples were used for MICAL-L2 immunohistochemical assays. After dewaxing and hydration, the microarray was incubated with 3% H2O2 for 30 min, subjected to antigen retrieval with citric acid at 95 °C for 20 min, blocked for 2 h at room temperature, incubated with primary antibody against MICAL-L2 at 4 °C overnight, and then with a species-matched secondary antibody for 2 h at room temperature. DAB staining was employed to detect the expression of MICAL-L2, with hematoxylin serving as the counterstain. Images were captured using an Olympus BX51 microscope. The immunoreactivity score (IRS) was obtained by multiplying the percentage of stained cells by the staining intensity scores of MICAL-L2, as previously described [19, 20].

Enrichment analysis for MICAL-L2 function

An ordered list of genes was generated based on the correlation between all genes and MICAL-L2 expression. Enriched pathways were determined using Gene Ontology (GO) [21, 22], KEGG [23,24,25], and GSEA [26, 27]. In the KEGG analysis, genes were determined to be differentially expressed based on a log2 fold-change of > 1.0 and an adjusted P-value < 0.05 (www.kegg.jp/kegg/kegg1.html), as previously reported [28]. “GSEA is a computational method that determines whether an a priori defined set of genes shows statistically significant, concordant differences between two biological states” [26, 27]. In this study, the predefined gene set was obtained from the MSigDB database (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp). STRING (https://cn.string-db.org/) and Cytoscape were used to predict and display the protein-protein interaction network of MICAL-L2 co-expressed genes.

Immune cell infiltration analysis using single-sample GSEA

Immune infiltration analysis of COAD tissue was performed using single-sample GSEA (ssGSEA) [27, 29]. The infiltration levels of 24 immune cell types were quantified from gene expression profiles, as previously described [30]. In addition, Spearman’s correlation was used to investigate the association between MICAL-L2 expression and immune cell infiltration.

Statistical analysis

SPSS 22.0 software was used for statistical analysis. The chi-square test was used to analyze and compare the clinical and pathological conditions of the two groups. The Kaplan–Meier method was used to evaluate the survival of patients and the log rank test was used to test the significance. A Cox proportional hazards regression model was used to identify significant and independent prognostic factors for COAD patients. Finally, R language was used to draw a nomogram and build a prediction model. P < 0.05 indicates significance (two-tailed).

Results

MICAL-L2 is highly expressed in COAD samples

Data mining in TCGA database showed that the mRNA expression of MICAL-L2 was elevated in most types of cancer (Fig. 1A). Focusing on COAD, a common histological subtype of colon cancer, we then examined the expression of MICAL-L2 in 456 COAD samples and 41 adjacent normal tissue samples from TCGA. We found that the mRNA expression of MICAL-L2 was significantly upregulated in COAD tissues compared with that in adjacent normal tissues (P < 0.001) (Fig. 1B). Similarly, in 41 paired cancerous and adjacent normal tissues, MICAL-L2 mRNA expression was also markedly higher in the cancer samples (P < 0.001) (Fig. 1C). Receiver operating characteristic (ROC) curve analysis was also applied to evaluate the diagnostic value of MICAL-L2 expression levels in COAD, and the area under the curve (AUC) was found to be 0.755 (95% CI = 0.691–0.819) (Fig. 1D).

Fig. 1
figure 1

mRNA expression data for MICAL-L2 in colon adenocarcinoma (COAD) obtained from The Cancer Genome Atlas (TCGA) database. A MICAL-L2 expression levels in different human tumor types according to TCGA data. B The differential expression of MICAL-L2 between COAD tissues and adjacent normal tissues. C The differential expression of MICAL-L2 between COAD tissues and paired adjacent normal tissues. D The receiver operating characteristic (ROC) curve for MICAL-L2 shows promising discrimination power between normal and COAD samples. **: P < 0.01, ***: P < 0.001

Analysis of the HPA and UALCAN data showed that the protein expression level of MICAL-L2 was higher in COAD tissues than in normal adjacent tissues (Fig. 2A, B). MICAL-L2 protein levels were also analyzed in a tissue microarray containing COAD and paracancerous tissues. Although some signal was lost during sample preparation, the immunohistochemical analysis nevertheless showed that MICAL-L2 protein levels were significantly higher in COAD tissues than in paracancerous normal tissues (Fig. 2C). Combined, these results indicated that MICAL-L2 is highly expressed in COAD at both the mRNA and protein levels.

Fig. 2
figure 2

The protein expression of MICAL-L2 in colon adenocarcinoma (COAD) tissues. A Differential protein expression of MICAL-L2 between COAD and adjacent normal tissues. B Representative images of MICAL-L2 expression in COAD. C Representative images of MICAL-L2 staining in COAD tissues. A scatterplot showing correlations between protein levels in COAD or paracancerous tissue as determined by immunoreactivity scores (IRS). ***: P < 0.001

Correlation between MICAL-L2 expression and clinicopathological features

The characteristics of 454 patients with COAD, including gene expression and clinical data, were collected from TCGA database. The patients were divided into high and low MICAL-L2 expression groups based on the mean value of MICAL-L2 expression (Table 1), following which putative correlations between MICAL-L2 expression and clinical characteristics were evaluated using logistic regression analysis. The results showed that MICAL-L2 mRNA expression was significantly associated with lymphatic invasion and primary therapy outcome (progressive disease [PD] + stable disease [SD] + partial response [PR] vs complete response [CR]) (Table 2).

Table 1 Association between MICAL-L2 expression and clinicopathologic features in the validation cohort
Table 2 Association between MICAL-L2 expression and clinicopathologic characteristics (logistic regression analysis)

Prognostic value of MICAL-L2 in COAD patients

We next determined the prognostic value of MICAL-L2 in COAD. For this, we evaluated the relationship between MICAL-L2 expression and clinical follow-up data using Kaplan–Meier analysis. Significance was assessed using the log rank test. The results showed that high MICAL-L2 expression was negatively correlated with OS (n = 453, P = 0.006; Fig. 3A) and DSS (n = 437, P = 0.028; Fig. 3B), indicating that MICAL-L2 expression levels were significantly associated with the prognosis of COAD patients. High MICAL-L2 expression was associated with poor OS in COAD patients who were over the age of 65, had stage T3 and T4 disease, or were female (Fig. 3C–H).

Fig. 3
figure 3

Prognostic value of MICAL-L2 expression for clinical outcomes in colon adenocarcinoma (COAD) patients. A&B Kaplan–Meier analysis of overall survival (OS) (A) and disease-specific survival (DSS) (B) in patients with COAD. CH Kaplan–Meier analysis of subgroup OS in patients with COAD. (C) Age ≤ 65, (D) stages T1 and T2, (E) male, (F) age > 65, (G) stages T3 and T4, and (H) female

To further identify the risk factors associated with OS in patients with COAD, univariate and multivariate analyses were performed using TCGA–COAD dataset. Univariate analysis showed that T stage, N stage, M stage, age, lymphatic invasion, and MICAL-L2 expression were the factors influencing OS. The multivariate analysis showed that MICAL-L2 expression (P = 0.032), age, and M stage were independent risk factors for OS (Table 3 & Fig. 4). Combined, these data suggested that MICAL-L2 may serve as a biomarker for the prediction of OS among COAD patients.

Table 3 Univariate and multivariate Cox proportional hazards analysis of the correlation between overall survival (OS) and MICAL-L2 expression levels
Fig. 4
figure 4

Forest plot of the multivariate Cox regression analysis for overall survival (OS) in colon adenocarcinoma (COAD) patients

Based on multivariate Cox regression analysis for OS, a nomogram was generated for internal validation. Prediction models were constructed for 1-, 3-, and 5-year OS in patients with COAD (Fig. 5A) while calibration plots to validate the efficiency of the nomograms for predicting OS were also generated. As shown in Fig. 5B, the bias-corrected line in the calibration plot was close to the ideal curve, indicative of an intimate relationship between the observed and predicted values.

Fig. 5
figure 5

Construction and validation of a nomogram to show the relationship between MICAL-L2 and other clinical factors with survival probability. A A nomogram for predicting the probability of 1-, 3-, and 5-year overall survival (OS) in colon adenocarcinoma (COAD) patients. B Calibration plots validating the efficiency of the nomograms for OS

Function enrichment analysis of MICAL-L2 in COAD

As we found that COAD patients with high levels of MICAL-L2 expression have worse OS and DSS than those with low MICAL-L2 expression, we explored the possible underlying cellular mechanism through KEGG and GSEA. As shown in Fig. 6A, 434 differentially expressed genes (DEGs) (|logFC| > 1, adjusted P-value < 0.05) were identified between the high and low MICAL-L2 expression groups, including 338 that were upregulated and 96 that were downregulated. The 10 genes showing the greatest positive or negative correlation with MICAL-L2 expression are shown in Fig. 6B. A network of potential co-expressed genes of MICAL-L2 in COAD are shown in Fig. S1.

Fig. 6
figure 6

Function enrichment analysis of MICAL-L2 in colon adenocarcinoma (COAD) using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). A Volcano plot of differentially expressed genes (DEGs) between the high and low MICAL-L2 expression groups. B Heat map showing the top 10 DEGs between the high and low MICAL-L2 expression groups. C GO and KEGG analysis of the DEGs. D Enrichment plots from GSEA. Hallmark_kras_signaling_DN (left) and hallmark_apical_junction (right)

The identified DEGs were submitted to GO term and KEGG pathway enrichment analysis. The following biological processes were found to be significantly affected: Chylomicron assembly, triglyceride-rich lipoprotein particle remodeling, and regulation of sensory perception of pain. The most enriched cellular component terms were apical plasma membrane, apical part of cell, and chylomicron. For molecular function, the most enriched terms were passive transmembrane transporter activity, channel activity, substrate-specific channel activity. The most enriched KEGG terms were cholesterol metabolism, neuroactive ligand–receptor interaction, and bile secretion (Fig. 6C). The GSEA results indicated that the co-expressed genes were mainly associated with the hallmark_kras_signaling_DN and hallmark_apical_junction pathways (Fig. 6D). We will further explore these pathways in future studies to better understand the function of MICAL-L2 in COAD.

Correlation between immune cell infiltration and MICAL-L2 expression levels in TCGA

The correlation between MICAL-L2 expression and immune infiltrate abundance in COAD was evaluated by ssGSEA using Spearman’s correlation tests (Fig. 7A). As shown in Fig. 7A, CD56bright natural killer (NK) cells, regulatory T cells (Tregs), and NK cells were all positively correlated with MICAL-L2 expression, whereas the opposite was seen for T-helper (Th) cells, gamma delta T (Tgd) cells, and Th2 cells. We further evaluated the infiltration levels of CD56bright NK cells, which displayed the greatest positive correlation with MICAL-L2 expression. The results showed that MICAL-L2 was significantly and positively correlated with the infiltration levels of CD56bright NK cells (P < 0.01, Fig. 7B, C).

Fig. 7
figure 7

Analysis of the association between MICAL-L2 expression and immune cell infiltration in colon adenocarcinoma (COAD). A Association between MICAL-L2 expression and 24 immune infiltrating cells. B Enrichment of CD56bright NK cells between high and low MICAL-L2 expression groups. C Correlation between MICAL-L2 and CD56bright NK cells. aDCs, activated DCs; DCs, dendritic cells; iDCs, immature DCs; NK, natural killer; pDCs, plasmacytoid DCs; Tcm, T central memory; Tem, T effector memory; Tfh, T follicular helper; Tgd, T gamma delta; Th, T helper cells; Th1, type 1 Th cells; Th2, type 2 Th cells; Th17, type 17 Th cells; Treg, regulatory T cell

Discussion

While there is only one MICAL-encoding gene in Drosophila, vertebrate genomes express genes encoding three MICAL (MICAL1–3) and two MICAL-like (MICAL-L1, MICAL-L2) isoforms. The disruption of MICAL1 activity was shown to impair cytoskeleton organization and breast tumor growth in an orthotopic model [31]. Additionally, high MICAL2 expression has been associated with lymphatic metastasis and shorter OS in lung cancer patients [32]. The three MICAL isoforms (MICAL1–3) contain a FAD domain and exhibit flavoprotein monooxygenase catalytic activity. Of note, MICAL1 exerts its effect on proliferation via reactive oxygen species (ROS)-sensitive PI3K/AKT/ERK signaling in breast cancer cells [33]. Similarly, MICAL2-induced ROS generation has also been reported to enhance the migratory potential of gastric cancer cells [34]. MICAL-L2 lacks the FAD domain and cannot generate ROS [35], and although several studies have suggested that MICAL-L2 may positively influence cancer progression [10, 13, 36], whether and how MICAL-L2 may be involved in this process remains unclear.

MICAL-L2 has been shown to be significantly upregulated in ovarian cancer tissues in a FIGO stage-dependent manner and has also been associated with histologic subgroups of ovarian cancer [11]. Consistent with these observations, our results revealed that, compared with adjacent normal tissues, MICAL-L2 expression was significantly upregulated in COAD tissues at both the mRNA and protein levels. ROC curve analysis also confirmed the diagnostic value of MICAL-L2. These findings strongly suggested that MICAL-L2 may play an oncogenic role in COAD. Accordingly, we assessed the prognostic value of MICAL-L2 in COAD using Kaplan–Meier survival analysis and found that patients with high MICAL-L2 expression have shorter OS and DSS compared with those with low MICAL-L2 expression. Univariate and multivariate analysis further revealed that high MICAL-L2 expression was an independent risk factor for OS in individuals with COAD. Collectively, these results indicated that MICAL-L2 may predict the prognosis of COAD and may represent a promising therapeutic target for the treatment of this cancer.

Abundant evidence supports that MICAL-L2 serves as a regulator of actin cytoskeleton organization, affecting processes such as cell vesicle trafficking and cellular morphology, among many other cytological behaviors [37,38,39]. For example, the interaction between MICAL-L2 and Rab8 and Rab13 was found to coordinate tight junction and adherens junction assembly [40, 41]. As expected, GO, KEGG, and GSEA indicated that MICAL-L2 localized to the apical part of the cell and regulated multiple types of transportation events. MICAL-L2 has been reported to mediate the endocytic recycling of occludin [42], while we have also recently shown that MICAL-L2 promotes the migration of gastric cancer cells via inhibiting EGFR transportation and degradation [10]. Although these findings demonstrated that MICAL-L2 might play a critical role in COAD progression through its transport-related activity, how MICAL-L2 precisely regulates cellular trafficking and then promotes COAD progression requires further exploration.

Over the past few years, the treatments used for COAD have transited from traditional chemical remedies to the use of targeted or immunotherapeutic drugs [43, 44]. In this study, ssGSEA was used to explore the association between MICAL-L2 expression and immune cell infiltration in COAD. Among the immune cell subpopulations, CD56bright NK cells showed the most enrichment in the high-MICAL-L2 expression group. NK cells are effective at killing tumors, and are commonly divided into CD56bright and CD56dim subtypes. Until recently, CD56bright NK cells were thought to exhibit potent antitumor activity [45, 46]. However, it has since been shown that CD56bright NK cells, which are enriched in human non-small-cell lung cancer infiltrate, display an impaired capability to kill tumor cells [47]. Further investigation, especially involving reciprocal activating crosstalk between immune cells and MICAL-L2, is necessary to delineate the regulatory immune mechanisms associated with MICAL-L2.

In conclusion, the results obtained in this study provide promising clues for a new mechanistic connection between MICAL-L2 expression and prognosis in COAD patients. Our findings indicated that MICAL-L2 may serve as an independent prognostic factor for patients with COAD, and further suggested that MICAL-L2 regulates cellular trafficking and promotes immune cell infiltration in COAD. However, how MICAL-L2 precisely regulates COAD progression remains to be characterized.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.

Abbreviations

AUC:

Area Under Curve

CI:

Confidence Interval

CH:

Calponin homology

CC:

C-terminal coiled-coil

COAD:

Colon Adenocarcinoma

DEGs:

Differentially expressed genes

DSS:

Disease Specific Survival

GSEA:

Gene Set Enrichment Analysis

HPA:

Human Protein Atlas Database

KEGG:

Kyoto Encyclopedia of Genes and Genomes

HR:

Hazard Ratio

LIM:

Lin11, Isl-1 and Mec-3

IRS:

Immunoreactivity score

MICALs:

Molecules Interacting With CasLs

NCBI:

National Center for Biotechnology Information

NES:

Normalized Enrichment Score

NK:

Natural killer

OR:

Odds Ratio

OS:

Overall Survival

ROC:

Receiver Operating Characteristic

ROS:

Reactive Oxygen Species

ssGSEA:

Single-sample Gene Set Enrichment Analysis

Tregs:

Regulatory T Cells

TNM:

Tumor Node Metastasis

TCGA:

The Cancer Genome Atlas

UALCAN:

University of Alabama Cancer Database

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Rubin DC, Shaker A, Levin MS. Chronic intestinal inflammation: inflammatory bowel disease and colitis-associated colon cancer. Front Immunol. 2012;3:107.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Akinkuotu AC, Maduekwe UN, Hayes-Jordan A. Surgical outcomes and survival rates of colon cancer in children and young adults. Am J Surg. 2021;221(4):718–24.

    Article  PubMed  Google Scholar 

  4. Araghi M, Arnold M, Rutherford MJ, Guren MG, Cabasag CJ, Bardot A, et al. Colon and rectal cancer survival in seven high-income countries 2010-2014: variation by age and stage at diagnosis (the ICBP SURVMARK-2 project). Gut. 2021;70(1):114–26.

    Article  PubMed  Google Scholar 

  5. Qaderi SM, Dickman PW, de Wilt JHW, Verhoeven RHA. Conditional survival and cure of patients with Colon or rectal Cancer: a population-based study. J Natl Compr Cancer Netw. 2020;18(9):1230–7.

    Article  Google Scholar 

  6. Yoon J, Terman JR. MICAL redox enzymes and actin remodeling: new links to classical tumorigenic and cancer pathways. Mol Cell Oncol. 2018;5(1):e1384881.

    Article  PubMed  CAS  Google Scholar 

  7. Giridharan SS, Caplan S. MICAL-family proteins: complex regulators of the actin cytoskeleton. Antioxid Redox Signal. 2014;20(13):2059–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Sun Y, Jaldin-Fincati J, Liu Z, Bilan PJ, Klip A. A complex of Rab13 with MICAL-L2 and alpha-actinin-4 is essential for insulin-dependent GLUT4 exocytosis. Mol Biol Cell. 2016;27(1):75–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Sakane A, Honda K, Sasaki T. Rab13 regulates neurite outgrowth in PC12 cells through its effector protein, JRAB/MICAL-L2. Mol Cell Biol. 2010;30(4):1077–87.

    Article  CAS  PubMed  Google Scholar 

  10. Min P, Zhao S, Liu L, Zhang Y, Ma Y, Zhao X, et al. MICAL-L2 potentiates Cdc42-dependent EGFR stability and promotes gastric cancer cell migration. J Cell Mol Med. 2019;23(6):4475–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Zhu LY, Zhang WM, Yang XM, Cui L, Li J, Zhang YL, et al. Silencing of MICAL-L2 suppresses malignancy of ovarian cancer by inducing mesenchymal-epithelial transition. Cancer Lett. 2015;363(1):71–82.

    Article  CAS  PubMed  Google Scholar 

  12. Ioannou MS, Bell ES, Girard M, Chaineau M, Hamlin JN, Daubaras M, et al. DENND2B activates Rab13 at the leading edge of migrating cells and promotes metastatic behavior. J Cell Biol. 2015;208(5):629–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Min P, Zhang L, Wang Y, Qi C, Song Y, Bibi M, et al. MICAL-L2 is essential for c-Myc Deubiquitination and stability in non-small cell lung Cancer cells. Front Cell Dev Biol. 2020;8:575903.

    Article  PubMed  Google Scholar 

  14. Vanoni MA. Structure-function studies of MICAL, the unusual multidomain flavoenzyme involved in actin cytoskeleton dynamics. Arch Biochem Biophys. 2017;632:118–41.

    Article  CAS  PubMed  Google Scholar 

  15. Gu H, Li Y, Cui X, Cao H, Hou Z, Ti Y, et al. MICAL1 inhibits colorectal cancer cell migration and proliferation by regulating the EGR1/beta-catenin signaling pathway. Biochem Pharmacol. 2022;195:114870.

    Article  CAS  PubMed  Google Scholar 

  16. Blum A, Wang P, Zenklusen JC. SnapShot: TCGA-analyzed tumors. Cell. 2018;173(2):530.

    Article  CAS  PubMed  Google Scholar 

  17. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Qi C, Min P, Wang Q, Wang Y, Song Y, Zhang Y, et al. MICAL2 contributes to gastric Cancer cell proliferation by promoting YAP Dephosphorylation and nuclear translocation. Oxidative Med Cell Longev. 2021;2021:9955717.

    Article  CAS  Google Scholar 

  19. Deng W, Gu L, Li X, Zheng J, Zhang Y, Duan B, et al. CD24 associates with EGFR and supports EGF/EGFR signaling via RhoA in gastric cancer cells. J Transl Med. 2016;14:32.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Dumitru CA, Bankfalvi A, Gu X, Zeidler R, Brandau S, Lang S. AHNAK and inflammatory markers predict poor survival in laryngeal carcinoma. Plos One. 2013;8(2):e56420.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Gene Ontology C. The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. 2021;49(D1):D325–34.

    Article  CAS  Google Scholar 

  23. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  PubMed  Google Scholar 

  26. Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.

    Article  CAS  PubMed  Google Scholar 

  27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang Q, Xia T, Qi C, Du J, Ye C. High expression of S100A2 predicts poor prognosis in patients with endometrial carcinoma. BMC Cancer. 2022;22(1):77.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462(7269):108–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yan Y, Nie K, Zheng J, Jiang X, Huang Y, Zheng Z, et al. High endothelin receptor type a expression as an independent prognostic biomarker and correlated with immune infiltrates in stomach adenocarcinoma. Cancer Manag Res. 2021;13:5013–26.

    Article  PubMed  PubMed Central  Google Scholar 

  31. McGarry DJ, Armstrong G, Castino G, Mason S, Clark W, Shaw R, et al. MICAL1 regulates actin cytoskeleton organization, directional cell migration and the growth of human breast cancer cells as orthotopic xenograft tumours. Cancer Lett. 2021;519:226–36.

    Article  CAS  PubMed  Google Scholar 

  32. Zhou W, Liu Y, Gao Y, Cheng Y, Chang R, Li X, et al. MICAL2 is a novel nucleocytoplasmic shuttling protein promoting cancer invasion and growth of lung adenocarcinoma. Cancer Lett. 2020;483:75–86.

    Article  CAS  PubMed  Google Scholar 

  33. Deng W, Wang Y, Zhao S, Zhang Y, Chen Y, Zhao X, et al. MICAL1 facilitates breast cancer cell proliferation via ROS-sensitive ERK/cyclin D pathway. J Cell Mol Med. 2018;22(6):3108–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wang Y, Min P, Qi C, Zhao S, Yu M, Zhang Y, et al. MICAL2 facilitates gastric Cancer cell migration via MRTF-A-mediated CDC42 activation. Front Mol Biosci. 2021;8:568868.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Giridharan SS, Rohn JL, Naslavsky N, Caplan S. Differential regulation of actin microfilaments by human MICAL proteins. J Cell Sci. 2012;125(Pt 3):614–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Sakane A, Yoshizawa S, Yokota H, Sasaki T. Dancing styles of collective cell migration: image-based computational analysis of JRAB/MICAL-L2. Front Cell Dev Biol. 2018;6:4.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Miyake K, Sakane A, Tsuchiya Y, Sagawa I, Tomida Y, Kasahara J, et al. Actin cytoskeletal reorganization function of JRAB/MICAL-L2 is fine-tuned by intramolecular interaction between first LIM zinc finger and C-terminal coiled-coil domains. Sci Rep. 2019;9(1):12794.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Sakane A, Abdallah AA, Nakano K, Honda K, Ikeda W, Nishikawa Y, et al. Rab13 small G protein and junctional Rab13-binding protein (JRAB) orchestrate actin cytoskeletal organization during epithelial junctional development. J Biol Chem. 2012;287(51):42455–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Rahajeng J, Giridharan SS, Cai B, Naslavsky N, Caplan S. Important relationships between Rab and MICAL proteins in endocytic trafficking. World J Biol Chem. 2010;1(8):254–64.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Yamamura R, Nishimura N, Nakatsuji H, Arase S, Sasaki T. The interaction of JRAB/MICAL-L2 with Rab8 and Rab13 coordinates the assembly of tight junctions and adherens junctions. Mol Biol Cell. 2008;19(3):971–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Nakatsuji H, Nishimura N, Yamamura R, Kanayama HO, Sasaki T. Involvement of actinin-4 in the recruitment of JRAB/MICAL-L2 to cell-cell junctions and the formation of functional tight junctions. Mol Cell Biol. 2008;28(10):3324–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Terai T, Nishimura N, Kanda I, Yasui N, Sasaki T. JRAB/MICAL-L2 is a junctional Rab13-binding protein mediating the endocytic recycling of occludin. Mol Biol Cell. 2006;17(5):2465–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Andre T, Shiu KK, Kim TW, Jensen BV, Jensen LH, Punt C, et al. Pembrolizumab in microsatellite-instability-high advanced colorectal Cancer. N Engl J Med. 2020;383(23):2207–18.

    Article  CAS  PubMed  Google Scholar 

  44. Overman MJ, McDermott R, Leach JL, Lonardi S, Lenz HJ, Morse MA, et al. Nivolumab in patients with metastatic DNA mismatch repair-deficient or microsatellite instability-high colorectal cancer (CheckMate 142): an open-label, multicentre, phase 2 study. Lancet Oncol. 2017;18(9):1182–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wagner JA, Rosario M, Romee R, Berrien-Elliott MM, Schneider SE, Leong JW, et al. CD56bright NK cells exhibit potent antitumor responses following IL-15 priming. J Clin Invest. 2017;127(11):4042–58.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Hwang S, Han J, Baek JS, Tak E, Song GW, Lee SG, et al. Cytotoxicity of human hepatic Intrasinusoidal CD56(bright) natural killer cells against hepatocellular carcinoma cells. Int J Mol Sci. 2019;20(7):1564.

    Article  CAS  PubMed Central  Google Scholar 

  47. Carrega P, Morandi B, Costa R, Frumento G, Forte G, Altavilla G, et al. Natural killer cells infiltrating human nonsmall-cell lung cancer are enriched in CD56 bright CD16(−) cells and display an impaired capability to kill tumor cells. Cancer. 2008;112(4):863–75.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (82073226), “Shilv” Plan founded by School of Basic Medical Sciences of NJMU.

Author information

Authors and Affiliations

Authors

Contributions

JD designed the study. YY, FY performed the statistical analysis. TX performed the immunohistochemistry assays. JD, QW,YZ drafted the manuscript. JD supervised the experimental work. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yujie Zhang or Jun Du.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with the Declaration of Helsinki. No ethics approval was required for this work. All utilized public data sets were generated by others who obtained ethical approval.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Figure S1.

Network of co-expressed genes of MICAL-L2.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, Y., Ye, F., Xia, T. et al. High MICAL-L2 expression and its role in the prognosis of colon adenocarcinoma. BMC Cancer 22, 487 (2022). https://doi.org/10.1186/s12885-022-09614-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09614-0

Keywords