- Primary research
- Open Access
Cofilin-1, LIMK1 and SSH1 are differentially expressed in locally advanced colorectal cancer and according to consensus molecular subtypes
Cancer Cell International volume 21, Article number: 69 (2021)
Colorectal cancer (CRC) is among the deadliest cancers, wherein early dissemination of tumor cells, and consequently, metastasis formation, are the main causes of mortality and poor prognosis. Cofilin-1 (CFL-1) and its modulators, LIMK1/SSH1, play key roles in mediating the invasiveness by driving actin cytoskeleton reorganization in various cancer types. However, their clinical significance and prognostic value in CRC has not been fully explored. Here, we evaluated the clinical contribution of these actin regulators according to TNM and consensus molecular subtypes (CMSs) classification.
CFL-1, LIMK1 and SSH1 mRNA/protein levels were assessed by real-time PCR and immunohistochemical analyses using normal adjacent and tumor tissues obtained from a clinical cohort of CRC patients. The expression levels of these proteins were associated with clinicopathological features by using the chi square test. In addition, using RNA-Seq data of CRC patients from The Cancer Genome Atlas (TCGA) database, we determine how these actin regulators are expressed and distributed according to TNM and CMSs classification. Based on gene expression profiling, Kaplan–Meier survival analysis was used to evaluated overall survival.
Bioinformatic analysis revealed that LIMK1 expression was upregulated in all tumor stages. Patients with high levels of LIMK1 demonstrated significantly lower overall survival rates and exhibited greater lymph node metastatic potential in a clinical cohort. In contrast, CFL-1 and SSH1 have expression downregulated in all tumor stages. However, immunohistochemical analyses showed that patients with high protein levels of CFL-1 and SSH1 exhibited greater lymph node metastatic potential and greater depth of local invasion. In addition, using the CMSs classification to evaluate different biological phenotypes of CRC, we observed that LIMK1 and SSH1 genes are upregulated in immune (CMS1) and mesenchymal (CMS4) subtypes. However, patients with high levels of LIMK1 also demonstrated significantly lower overall survival rates in canonical (CMS2), and metabolic (CMS3) subtypes.
We demonstrated that CFL-1 and its modulators, LIMK1/SSH1, are differentially expressed and associated with lymph node metastasis in CRC. Finally, this expression profile may be useful to predict patients with aggressive signatures, particularly, the immune and mesenchymal subtypes of CRC.
Colorectal cancer (CRC) it is the third most commonly diagnosed cancer and the second most common cause of cancer related deaths worldwide. Additionally, metastasis is the main cause of death in the majority of patients, and effective biomarkers for risk stratification of early metastatic disease are still required . Decisions regarding clinical treatment are almost exclusively based on clinicopathological parameters, since CRC is a highly heterogeneous disease molecular classifications have been proposed based on the gene expression signatures and tumor biology features, which could lead to improved therapeutic outcomes [2, 3]. Later, one research group performed a more detailed classification of primary colorectal tissue, in which a large dataset (n = 4151 patients) was analyzed and four consensus molecular subtypes (CMS) for CRC were proposed with distinguishing features: CMS1 (microsatellite instability and immunity), CMS2 (epithelial and canonical), CMS3 (epithelial and metabolic), and CMS4 (mesenchymal) . Hence, research has attempted to identify new phenotypic signatures in each subgroup to design an innovative strategy to predict patient stratification and to translate this information into a more homogeneous and effective drug response.
Particularly, the CMS4 displays a lower overall survival rate and relapse-free survival rate, exhibiting a high stromal content and activation of pathways related to epithelial–mesenchymal transition (EMT), mediated by transforming growth factor-β (TGFβ) signaling . EMT is a reversible cellular program in which epithelial cells transiently lose proteins of the junctional complex and acquire mesenchymal proteins which develop quasi-mesenchymal cell states, leading to more aggressive phenotypes [4, 5]. During EMT, the cellular actin cytoskeleton is rearranged and distinct actin structures are stimulated, allowing cell migration and invasion . This process is regulated by various actin-binding proteins that mediate the construction of protrusive and contractile structures, associated with rapid polymerization/depolymerization of filamentous actin (F-actin) [6, 7]. The members of cofilin-1 (CFL-1) and ADF protein family are the most critical components of actin dynamic regulation. CFL-1 promotes actin filament disassembly by severing F-actin and rapidly increasing depolymerization of the filament. CFL-1 activity is tightly regulated by LIM domain kinase 1 (LIMK1) and slingshot protein phosphatase 1 (SSH1) proteins . The precise balance between expression and subcellular localization of CFL-1, LIMK1, and SSH1 is pivotal for small changes in the dynamics of the actin cytoskeleton, and may augment features of tumor aggressiveness, including tumor dissemination . In various cancer types, including gastric, prostate, urothelial, breast, and vulvar, elevated expression of CFL-1 and its regulators, LIMK1/SSH1, have been implicated in tumor progression and aggressiveness, as well as in poor survival rate [10,11,12,13,14,15,16]. In the case of CRC, CFL-1 gene expression information is limited and divergent. Some studies have demonstrated elevated gene and protein expression, while others showed decreased expression in relation to normal tissues [17,18,19,20,21].
Since a significant proportion of patients present dissemination of tumor cells with subsequent metastasis growth, both innovative therapies and biomarkers to detect patients with lymph node metastasis potential are needed in CRC treatment. Therefore, an analysis of the profile of gene and protein expression of CFL-1 and its regulators, LIMK1/SSH1, may provide important insights into the local invasion processes of primary tumor cells and could help in disease stratification in routine pathology. In this study, using bioinformatic approaches and a clinical cohort, we aimed to determine the potential role of expression levels, as well as clinical implications, of CFL-1 and its regulators, LIMK1/SSH1, in CRC.
Materials and methods
Tumor tissue samples
Primary tumor specimens and their corresponding normal adjacent tissue (at least 5 cm away from the lesion) were obtained from patients with CRC who underwent colectomy and rectosigmoidectomy at the surgical center of the Hospital do Câncer I—Instituto Nacional de Câncer (INCA) (Rio de Janeiro, Brazil). Fresh samples were collected from 2008 to 2018 and were immediately stored in RNA-later (Invitrogen, Carlsbad, CA, USA) at -20 °C for RT-qPCR analysis (n = 38). Formalin-fixed paraffin-embedded CRC tissues (n = 59) from 2009 to 2018 were obtained from the Department of Pathology of INCA and were used for immunohistochemical (IHC) analysis after identification (primary tumors n = 49 and liver metastasis tissue n = 10). None of the patients enrolled had received preoperative radiotherapy or chemotherapy. Details of the clinicopathological description of the patients included in this study are available in Table 1. This study was approved by the Instituto Nacional de Cancer Ethics Research Committee (CEP-INCA), following all relevant guidelines and regulations. Informed consent was obtained from all participants.
Transcription polymerase chain reaction (RT-qPCR)
RNA was isolated by extraction using Trizol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol. The RNA integrity was evaluated by conventional PCR and agarose gel electrophoresis. RNA (1 µg) was treated with DNase I, RNase-free, (Thermo Scientific, Wilmington, NC, USA) to eliminate genomic DNA contamination. cDNA was synthesized from the specimens as described previously . Quantitative PCR was performed using the SYBR Green® PCR Master Mix (Invitrogen, Carlsbad, CA—USA) and Applied Biosystems 7500 Real-Time PCR System, using the following primers: LIMK1 F 5′-CAAGGGACTGGTTATGGTGGC-3′; LIMK1 R 5′-CCCCGTCACCGATAAAGGTC-3′; CFL-1 F 5′-GATAAGGACTGCCGCTATGC-3′; CFL-1 R 5′-GCTTGATCCCTGTCAGCTTC-3′; SSH1 F 5′-ACACCGAGGAGAATATCTTGC-3′; SSH1 R 5′-TGAACCCACCATCTCCATCAAG-3′; β-Actin F 5′-TACAATGAGCTGCGTGTGG-3′; and β-Actin R 5′-TAGCACAGCCTGGATAGCAA-3′. We performed the reaction using 30 s at 95 °C, 30 s at 60 °C, and 30 s at 72 °C for CFL-1 primers, 30 s at 95 °C, 30 s at 57 °C, and 30 s at 72 °C for SSH1 primers, and 30 s at 95 °C and 45 s at 65 °C for LIMK1 primers, all using 40 cycles [23,24,25]. Negative template controls were run in all qPCR experiments to evaluate whether reaction was contaminated with exogenous DNA or to detect primer dimers formation. The differential expression of the selected genes was calculated by the 2−ΔΔCT or 2−ΔCT method . Tumor samples were classified into high or low groups according to the median value.
Immunohistochemistry and pathological analysis
The histological analysis of cancer specimens was conducted according to the guidelines contained in the 8th edition of the Cancer Staging Manual, edited by the American Joint Committee on Cancer, prior to immunohistochemical analysis by pathologists of INCA . The patient cohort selected included cases of all tumor stages and pathological grades (I–IV). Immunohistochemistry analyses of CFL-1, LIMK1, and SSH1 proteins were performed using formalin-fixed paraffin-embedded CRC blocks from 59 patients (49 with primary tumors and 10 with liver metastases) stored in the Division of Pathology of INCA. Briefly, 3-micron slices were de-paraffinized and antigen retrieval was performed in Trilogy Buffer (Cell Marque, Sigma-Aldrich, Rocklin, CA, USA), at 98 °C, using the steam process for 30 min. Endogenous peroxidases were blocked with a NovoLink Max Polymer Detection kit (Leica Microsystems, Wetzlar, Germany) for 5 min. Anti-LIMK1 (Thermo Fisher—#PA5-14938, 1:1400 dilution), anti-SSH1 (Sigma—#HPA019845, 1:600 dilution), and anti-CFL-1 (Cell Signaling Technology—#5175, 1:2500 dilution) antibodies were incubated overnight at 4 °C. After incubation, the post-primary antibody and the polymer (Novolink, Newcastle upon Tyne, United Kingdom) were added and incubated for 30 min, rinsed, and exposed to a solution of diaminobenzidine for 3 min. Next, the samples were dehydrated with alcohol, cleared in xylene, and mounted. Negative controls were acquired using the same protocol described above, with the omission of the primary antibody. Immunoreactivity (IR) was determined under an optical microscope based on evaluation of staining intensity (on a scale of 1–3) and percentage of tumor cells with positive IR staining (0 to 100%) for each antibody on each slide. The value was calculated by a pathologist (unfamiliar with the experimental groups) by multiplying staining intensity by percentage of cells with positive IR. The median score among all patients was used for separation between high and low expression.
CMS classification and gene expression analysis
We downloaded the RNA-Seq data from human CRC samples (n = 622, TCGA-COAD, and TCGA-READ) and normal colorectal samples (n = 51) databases available through the Cancer Genome Atlas (TCGA) project using the R package TCGA Biolinks . Patients were selected when they had RNA-Seq data available online through the TCGA Biolinks package, and the exclusion criterion was not having been classified with the CRC molecular subtype classifier tool  or not having the staging described in the clinical-pathological information table. Differences in CFL-1, SSH1, and LIMK1 gene expression among the CMS groups were evaluated after samples were classified using the CRC molecular subtype classifier . We also analyzed gene expression according to the tumor stages of the samples when clinicopathological information was available. Correlations between gene pairs were evaluated as well as their impact on patient survival. All analyses were performed in the R environment, and plots were constructed using the ggplot2 package .
Patient cohort groups were analyzed using P values calculated from Student’s t-test where two groups were compared. For comparison between more than two groups, ANOVA followed by the Bonferroni post-test were used. The association between CFL-1, SSH1, and LIMK1 expression and clinicopathological parameters was analyzed by χ2 or Fisher’s exact tests. Data were analyzed using GraphPad Prism 5.0 (GraphPad Software, San Diego, CA, USA). Lymph node metastasis risk was analyzed using the binary logistic regression method, using univariate and multivariate analysis. Data were analyzed using SPSS version 20.0 software (SPSS; Chicago, IL, USA). Differences in expression levels from TCGA data were analyzed using the Kruskal–Wallis test to compare differences between CMSs and tumor stages. This was followed by post-hoc pairwise comparisons using the Dunn’s test to determine significant differences between each subgroup. The correlation between CFL-1 and SSH1 or LIMK1 was evaluated using the Spearman correlation coefficient test. Overall survival curves were plotted by the Kaplan–Meier method using the R packages for survival analysis “survival” , “survminer”, and the log-rank test was used to compare survival estimates across different groups. Hazard ratios (HRs) were adjusted using Cox regression, with a 95% confidence interval. To overall survival analysis, samples expressing the gene above or equal the median were classified as high expression and the others as low expression. Differences were considered statistically significant when *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001.
Overexpression of LIMK1 predicts worse overall survival in CRC patients
To investigate the role that CFL-1 and its regulators, LIMK1 and SSH1, play in human CRC we analyzed the mRNA expression patterns of their coding genes using an in-silico analysis of RNAseq data from patients with CRC (TCGA database). This data comprised of both normal (n = 51) and tumor (n = 622) samples. First, we compared the gene expression levels of each gene among different tumor stages (I (n = 105), II (n = 227), III (n = 179), and IV (n = 88). Twenty-three samples were excluded from this analysis due to a lack of staging information. LIMK1 expression was upregulated (P < 0.0001) whereas CFL-1 and SSH1 were downregulated (P < 0.0001) in all tumor stages when compared with normal tissues (Fig. 1a). However, no significant differences were observed comparing expression levels of each individual gene among different tumor stages (P > 0.05). Since LIMK1 was overexpressed in all tumor stages, we analyzed the impact of these overexpression on overall survival of the patients by Kaplan-Meier curves analysis using TCGA database. Patients were separated according to all tumor stages (I-IV) (Fig. 1b), early stage (I-II) (Fig. 1c), and late-stage (III-IV) (Fig. 1d). A significant worse overall survival rate was observed for patients with high levels of LIMK1 in all tumor stages, early stage, and late stage (P < 0.001, P < 0.01, and P < 0.01, respectively), indicating that high level LIMK1 was an unfavorable prognostic indicator in patients with CRC. However, there were no significant differences in overall survival among subgroups of patients and expression levels of SSH1 (Additional file 1: Figure S1A-B) and to CFL-1 in early stages (Additional file 2: Figure S2A). Patients with low levels of CFL-1 (late-stage III–IV) exhibited lower overall survival (P < 0.05) (Additional file 2: Figure S2B).
Overexpression of LIMK1 is associated positive lymph node metastasis in CRC patients
To confirm the results obtained by bioinformatic analysis, we evaluated CFL-1, LIMK1, and SSH1 mRNA levels by RT-qPCR using tumor tissues from patients with CRC paired with adjacent normal tissues in a clinical cohort. Heterogeneous expression for CFL-1, SSH1 and LIMK1 were detected when tumor tissues were paired with adjacent normal samples (Fig. 2a–c). However, when we separate tumor samples in groups with high or low expression according to the median value obtained by RT-qPCR analysis, only those with higher levels of LIMK1 were associated with lymph node metastasis (N1–N2) (χ2, 8.081; P < 0.005) (Fig. 2d). These results suggest that high mRNA level of LIMK1 is associated with positive regional lymph node metastasis in CRC patients.
CFL-1 and SSH1 protein levels are associated with positive lymph node metastasis in CRC patients
We also explored the protein expression patterns in 43 normal colorectal tissues, 49 primary tumors (stages I–IV), and 10 liver metastases by immunohistochemistry (IHC) analysis. CRC samples displayed CFL-1 immunostaining localized in the cytoplasm and membrane, a small number of cells also exhibited nuclear localization. Adjacent normal epithelial tissues displayed either membrane and/or cytoplasmic CFL-1 in the glands of the intestinal crypt (Fig. 3a). Quantitative analysis of CFL-1 staining revealed a higher level (2.7 median score) of CFL-1 in stage III CRC. Stage I patients’ tissues had low-level expression of CFL-1, and exhibited scores of less than 1 (10/13 cases). In patients in stages II, IV, and in liver metastases tissues, CFL-1 expression varied widely, although some isolated cells clearly showed positive staining, even in tissues with low expression levels (Fig. 3b).
To further analyze the association between the data of CFL-1 obtained by IHC score with clinicopathological parameters, we separate tumor samples into high or low expression groups according to the median value obtained in the IHC score. CFL-1 score was significantly related to depth of invasion and metastatic lymph nodes (P = 0.001 and P = 0.02, respectively) (Table 2). Patients with locally advanced invasion (through the muscularis propria into pericolorectal tissues—T3), exhibited a higher CFL-1 score (81.3%). Furthermore, patients with lymph node metastasis also exhibited a higher CFL-1 score (73.7%). However, no significant association between sex, age, localization, stage, lymphovascular and perineural invasion, liver metastasis, tumor size, tumor type, and tumor grade with CFL-1 IHC score was identified (Table 2).
Similarly, IHC staining for SSH1 displayed localization in the cytoplasm and membrane in the normal adjacent epithelial layer. IHC quantification revealed high levels (0.8 median score) of SSH1 in stage III tumor tissues, while tissues in stages I, II, IV and in metastatic tissue revealed weak focally positive staining intensity in groups of tumor cells (Fig. 4a, b). When tumor samples were separate into high or low SSH1 expression groups, the SSH1 IHC score was significantly related to lymph node metastasis, similar to CFL-1 data (Table 3). In fact, patients expressing higher levels of SSH1 exhibited a high percentage of T3 grade (68.7%), and lymph node tumor positivity was associated with the SSH1 IHC score (P = 0.02). However, a non-significant association between other clinicopathological parameters and the SSH1 IHC score was observed (Table 3). These data suggest that both CFL-1 and SSH1 play an important role in promoting tumor migration and invasion, leading to dissemination of tumor cells.
LIMK1 immunostaining also occurred in the glands of the intestinal crypt with cytoplasmic localization in adjacent normal tissues. IHC quantitative analysis of LIMK1 did not reach statistical significance because the immunoreactivity levels varied widely among tissue samples (Fig. 5a, b). Furthermore, no association was observed between the LIMK1 IHC score with clinicopathological parameters (Table 4).
To address whether CFL-1 and SSH1 protein levels could be prognostic factors for lymph node metastasis, we performed univariate and multivariate binary logistic regression analyses. In univariate analysis, parameters such as T grade (P = 0.011), lymphovascular invasion (P = 0.012), and high levels of CFL-1 and SSH1 (P = 0.0125) were significantly correlated with lymph node metastasis risk (Table 5). However, multivariate analysis did not support high levels of CFL-1 and SSH1 as independent prognostic factors for CRC (Additional file 3: Table S1). These results suggest that between the clinicopathological factors, depth of invasion, and lymphovascular invasion, the protein levels of CFL-1 and SSH1 could be important factors in predicting lymph node metastases and local advanced disease in CRC.
LIMK1 and SSH1 are differentially expressed among the consensus molecular subtypes (CMS)
Next, we sought to determine how CFL-1 and its regulators, LIMK1 and SSH1, are expressed and distributed according to the CMS classification. RNA-Seq data of CRC patients from the TCGA database were used to classify each sample according to the CMS system: 68 samples were classified as CMS1, 155 as CMS2, 84 as CMS3, and 161 as CMS4. One hundred fifty-four samples were not classified. LIMK1 and SSH1 genes, but not CFL-1 gene showed significant differential expression among the CMSs [Kruskal–Wallis P value < 0.05 for LIMK1 and SSH1and P value > 0.05, to CFL-1 (Fig. 6a)]. Notably, LIMK1 expression was highest in CMS4 (median of 10.78) compared to CMS2 and CMS3 (median of 10.60 and 10.43 respectively, P < 0.0001 for both), but was similar to CMS1. In addition, the expression level of SSH1 was similar in CMS2 and CMS3 (median of 10.31 and 10.30 respectively), while CMS1 (median of 10.61) significantly express more SSH1 when compared with CMS2 and CMS3 (P < 0.0001 for both). CMS4 (median of 10.77) exhibited also higher levels of SSH1 expression when compared to the CMS3 and CMS2 subtypes (P < 0.0001 for both). These results revealed that LIMK1 and SSH1 genes are upregulated in immune and mesenchymal subtypes, suggesting a distinct role in the actin dynamic regulation during tumor progression of CRC.
Finally, we analyze the correlation between CFL-1 levels and the regulators SSH1 or LIMK1 according to the CMS classification. We observed that all correlations between CFL-1 and SSH1 mRNA levels were negative (Fig. 6b). In addition, the CMS4 subtype exhibited the highest correlation value (P < 0.0001 and R = − 0.4) among the subtypes. CMS1 (P < 0.05) and CMS2 (P < 0.0001) also exhibited significant correlation between CFL-1 and SSH1 mRNA levels (Fig. 6b). In contrast, all correlations between CFL-1 and LIMK1 mRNA levels were positive. Furthermore, CMS3 and CMS4 exhibited significant correlations between CFL-1 and LIMK1 mRNA levels (P < 0.0001 for CMS3; P < 0.05 for CMS4), and CMS3 exhibited the highest correlation value (R = 0.47) among the subtypes (Fig. 6c). Interestingly, high levels of LIMK1 expression were significantly correlated with lower overall survival rates in the CMS2 and CMS3 subtypes (P < 0.05, Fig. 6d). These results revealed that overexpression of LIMK1 could be useful to detect patients with poor prognosis in CMS2 and CMS3 subtypes. However, there were no significant differences in overall survival among CMSs expression levels of CFL-1 and SSH1 (Additional file 4: Figure S3 and Additional file 5: Figure S4).
Invasion of tumor cells is a key trait of the aggressive phenotype of cancers, and greatly affects the survival of patients due to metastatic outgrowth of disseminated tumor cells . There is considerable evidence indicating that CFL-1 and its regulators, LIMK1 and SSH1, are convergent points of cell signaling networks that play a crucial role in modulating actin cytoskeleton reorganization during malignant progression [8, 30]. However, the clinical significance of these proteins in CRC is not yet fully understood.
We found that LIMK1 mRNA was upregulated in all tumor stages, and patients with high levels exhibited lower overall survival. Although tumor tissue data from our clinical cohort showed heterogeneous expression of LIMK1 mRNA, positive regional lymph node metastasis were positively associated with high levels of LIMK1 expression. These data indicate that LIMK1 could be associated with local invasion during CRC progression and could participate in regional lymph node metastasis. These findings are consistent with other studies in CRC showing that LIMK1 mRNA and protein levels were associated with poor prognosis, including depth of invasion, lymph node metastasis, tumor stage, distal metastases, and reduced survival rate [19, 31, 32]. The mechanism by which LIMK1 can influence poor survival in CRC patients is not known. However, it is possible to suggest that LIMK1 could coordinate proliferation and dissemination of tumor cells by both microtubule disassembly and actin polymerization [33, 34]. Reinforcing this suggestion, it was reported that LIMK1 is crucial for invasiveness and metastatic activity, and regulating the phosphorylation of CFL-1 during mitosis contributes to appropriate cytokinesis, proliferation, and extension of the tumor [34, 35].
Regarding the clinical and biological significance of CFL-1 and SSH1 in CRC, we found that CFL-1 and SSH1 mRNA levels were downregulated in all tumor stages when analyzed by bioinformatic. Indeed, in our clinical cohort, some patients showed low expression of CFL-1 and SSH1 mRNA. Conversely, the upregulation of SSH1 mRNA levels has been observed in CRC tumor tissues when compared with control tissues, and was also associated with tumor stage, lymph node metastasis, and poor prognosis . In addition, SSH1 upregulation was also observed in breast and pancreatic tumors. In breast tumors, SSH1 upregulation was associated with increased metastasis and mortality [15, 36]. In contrast with our data, CFL-1 was upregulated in another small cohort (30 patients) of CRC patients; however, no correlation was found between CFL-1 mRNA levels and clinicopathological features of patients . Interestingly, a study demonstrated progressively reduced levels of phosphorylation, but not mRNA expression, of CFL-1 in CRC patient tissues . Indeed, the low expression levels of CFL-1 and SSH1 in CRC may be explained by the number of patients in each study and particularly the level of methylation in the promoter region that regulates gene expression. In according, we found increased methylation level of the CFL-1 and SSH1 promoter regions in colon and rectum adenocarcinoma tissues from the TCGA database, suggesting that differences in the expression of CFL-1 and SSH1 in CRC may be related to epigenetic modifications (Additional file 6: Figure S5). Supporting this idea, LIMK2, another actin regulator, was downregulated in CRC tissues due to increased promoter methylation .
Our IHC profiling identified a direct association between CFL-1 and SSH1 levels with lymph node metastasis in CRC tissues, which was significantly correlated with lymph node metastasis risk by univariate analysis, supporting their role in local tumor dissemination. The results obtained using the multivariate analysis were not significative possibly due to small number of patients enrolled in this study. Further analysis is needed to assess their role in lymph node metastasis risk. However, CFL-1 levels association with dedifferentiation, infiltration depth, pathological staging, and presence of lymph node metastasis have also been found in other cancer types [11,12,13, 37, 38]. It is possible to speculate that the invasive margin may be region within the tumor that exhibits the highest levels of CFL-1, which increases aggressiveness and leads to local dissemination. However, further investigation is required to determine the relationship between CFL-1 and SSH1 expression levels and their role in different areas of the tumor.
Our bioinformatic analysis showed that the CMS1 and CMS4 subtypes exhibited the highest levels of LIMK1 and SSH1 expression among the analyzed subtypes, while CFL-1 levels were homogeneous according to CMSs. It is important point out that CMS1 and CMS4 subtypes are characterized by an extensive prominent stromal invasion, fibroblast-rich content, angiogenesis, and poor outcome among the CMSs. CMS1 subtype have the worst survival after relapse while CMS4 exhibits the worst prognosis subtype . In agreement with our results, previous reports have shown a relevant association between these proteins and EMT markers [E-cadherin, β-catenin (nuclear), ZEB1, and SNAIL] in CRC tissues [19, 21], indicating their central role during EMT in CRC. We previously demonstrated, using an in vitro model, the dynamic role of CFL-1 and LIMK2 during EMT in CRC . In accordance with these data, studies using bladder, prostate, and gastric models have also correlated CFL-1 and LIMK1 expression pathways with the EMT program [14, 39, 40], as well as tumor initiation and metastatic colonization in breast model . Notably, these gene expression could not be constant throughout the different areas of the tumor probably due to the high intratumor heterogeneity of CRC . One study using CRC tissues showed that tumor buddings exhibited differential gene expression related to the actin cytoskeleton remodeling pathway, including EMT signatures when compared with tumor bulks, and also alternated from CMS2 in the tumor bulk to CMS4 in budding cells . Furthermore, we noticed a significant correlation between CFL-1 and LIMK1/SSH1 expression according to the CMS classification suggesting that the modulation of CFL-1 (activation/non-activation) can be more complex and regulated to different degrees depending on the CMS. Aggressive subtypes, such as mesenchymal (CMS4), may have a higher degree of dependency on CFL-1 regulators, and could impact on higher invasion rates.
Moreover, we observed that high levels of LIMK1 mRNA significantly correlated with lower overall survival in canonical and metabolic subtypes. The major pathway related to the canonical subtype is Wnt signaling. Interestingly, evidence has shown that LIMK1 can bind β-catenin protein and translocate into the nucleus, enhancing CRC cell metastasis . Regarding metabolic LIMK1 functions, some studies have suggested that insulin and high glucose levels may stimulate actin cytoskeleton remodeling, with consequent phosphorylation of LIMK1 mediating GLUT4 translocation to the cell surface and glucose uptake [43,44,45]. However, further investigation is required to verify the role of LIMK1 in canonical and metabolic subtypes in CRC since most of these studies were evaluated using in vitro models. Here, we used bioinformatics approaches to show the clinical significance of these proteins and their relationship with distinct biological phenotypes, according to the CMS classification. Additionally, since these proteins could mediate the dissemination of tumor cells, their mRNA and correlation analyses may be useful for stratifying patients with potential risk for metastatic disease in each CMS.
In conclusion, our data suggest that CFL-1 and its regulators LIMK1 and SSH1, are differentially expressed in CRC. Since patients with lymph node metastasis exhibited RNA and/or protein upregulation, these might indicate their crucial role in regional cancer cell dissemination and also might be useful for predicting metastatic disease in CRC. A better and more detailed clinical stratification by CMS, as well as new treatment approaches using these proteins as targets, might provide improvements in therapeutic outcomes for CRC patients, avoiding early tumor dissemination, particularly in immune and mesenchymal subtypes.
Adjacent normal tissue
The American Joint Committee on Cancer
Consensus Molecular Subtypes
LIM domain kinase 1
Slingshot protein phosphatase 1
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Guinney J, Dienstmann R, Wang X, de Reyniès A, Schlicker A, Soneson C, et al. The consensus molecular subtypes of colorectal cancer. Nat Med. 2015;21:1350–6.
Isella C, Brundu F, Bellomo SE, Galimi F, Zanella E, Porporato R, et al. Selective analysis of cancer-cell intrinsic transcriptional traits defines novel clinically relevant subtypes of colorectal cancer. Nat Commun. 2017;8:15107.
Dongre A, Weinberg RA. New insights into the mechanisms of epithelial-mesenchymal transition and implications for cancer. Nat Rev Mol Cell Biol. 2019;20:69–84.
Yang J, Antin P, Berx G, Blanpain C, Brabletz T, Bronner M, et al. Guidelines and definitions for research on epithelial-mesenchymal transition. Nat Rev Mol Cell Biol. 2020;21:341–52.
Kristó I, Bajusz I, Bajusz C, Borkúti P, Vilmos P. Actin, actin-binding proteins, and actin-related proteins in the nucleus. Histochem Cell Biol. 2016;145:373–88.
Hodge RG, Ridley AJ. Regulating Rho GTPases and their regulators. Nat Rev Mol Cell Biol. 2016;17:496–510.
Kanellos G, Frame MC. Cellular functions of the ADF/cofilin family at a glance. J Cell Sci. 2016;129:3211–8.
Wang W, Mouneimne G, Sidani M, Wyckoff J, Chen X, Makris A, et al. The activity status of cofilin is directly related to invasion, intravasation, and metastasis of mammary tumors. J Cell Biol. 2006;173:395–404.
You T, Gao W, Wei J, Jin X, Zhao Z, Wang C, et al. Overexpression of LIMK1 promotes tumor growth and metastasis in gastric cancer. Biomed Pharmacother. 2015;69:96–101.
Lu LI, Fu NI, Luo XU, Li X-Y, Li X-P. Overexpression of cofilin 1 in prostate cancer and the corresponding clinical implications. Oncol Lett. 2015;9:2757–61.
Maimaiti Y, Liu Z, Tan J, Abudureyimu K, Huang B, Liu C, et al. Dephosphorylated cofilin expression is associated with poor prognosis in cases of human breast cancer: a tissue microarray analysis. Onco Targets Ther. 2016;9:6461–6.
Wu Q, Jiang Y, Cui S, Wang Y, Wu X. The role of cofilin-l in vulvar squamous cell carcinoma: a marker of carcinogenesis, progression and targeted therapy. Oncol Rep. 2016;35:2743–54.
Hensley PJ, Zetter D, Horbinski CM, Strup SE, Kyprianou N. Association of epithelial-mesenchymal transition and nuclear cofilin with advanced urothelial cancer. Hum Pathol. 2016;57:68–77.
Chen C, Maimaiti Y, Zhijun S, Zeming L, Yawen G, Pan Y, et al. Slingshot-1L, a cofilin phosphatase, induces primary breast cancer metastasis. Oncotarget. 2017;8:66195–203.
Maimaiti Y, Maimaitiming M, Li Y, Aibibula S, Ainiwaer A, Aili A, et al. SSH1 expression is associated with gastric cancer progression and predicts a poor prognosis. BMC Gastroenterol. 2018;18:12.
Lourenço FC, Munro J, Brown J, Cordero J, Stefanatos R, Strathdee K, et al. Reduced LIMK2 expression in colorectal cancer reflects its role in limiting stem cell proliferation. Gut. 2014;63:480–93.
Zhang Y, Li A, Shi J, Fang Y, Gu C, Cai J, et al. Imbalanced LIMK1 and LIMK2 expression leads to human colorectal cancer progression and metastasis via promoting β-catenin nuclear translocation. Cell Death Dis. 2018;9:749.
Aggelou H, Chadla P, Nikou S, Karteri S, Maroulis I, Kalofonos HP, et al. LIMK/cofilin pathway and Slingshot are implicated in human colorectal cancer progression and chemoresistance. Virchows Arch. 2018;472:727–37.
Mousavi S, Safaralizadeh R, Hosseinpour-Feizi M, Azimzadeh-Isfanjani A, Hashemzadeh S. Study of cofilin 1 gene expression in colorectal cancer. J Gastrointest Oncol. 2018;9:791–6.
Song X, Xie D, Xia X, Tan F, Pei Q, Li Y, et al. Role of SSH1 in colorectal cancer prognosis and tumor progression. J Gastroenterol Hepatol. 2020. https://doi.org/10.1111/jgh.15001.
Bastos LG, dos R, de Marcondes, de-Freitas-Junior PG, Leve JCM, Mencalha F, de Souza AL. WF, et al. Progeny from irradiated colorectal cancer cells acquire an EMT-like phenotype and activate Wnt/β-catenin pathway. J Cell Biochem. 2014;115:2175–87.
Rubie C, Kempf K, Hans J, Su T, Tilton B, Georg T, et al. Housekeeping gene variability in normal and cancerous colorectal, pancreatic, esophageal, gastric and hepatic tissues. Mol Cell Probes. 2005;19:101–9.
Li D, Wang H, Song H, Xu H, Zhao B, Wu C, et al. The microRNAs miR-200b-3p and miR-429-5p target the LIMK1/CFL1 pathway to inhibit growth and motility of breast cancer cells. Oncotarget. 2017;8:85276–89.
Sousa-Squiavinato ACM, Rocha MR, Barcellos-de-Souza P, de Souza WF, Morgado-Diaz JA. Cofilin-1 signaling mediates epithelial-mesenchymal transition by promoting actin cytoskeleton reorganization and cell-cell adhesion regulation in colorectal cancer cells. Biochim Biophys Acta Mol Cell Res. 2019;1866:418–29.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3:1101–8.
Weiser MR. AJCC 8th Edition: Colorectal Cancer. Ann Surg Oncol. 2018;25:1454–5.
Colaprico A, Silva TC, Olsen C, Garofano L, Cava C, Garolini D, et al. TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016;44:e71.
Therneau TM until. 2009) TL (original S->R port and R maintainer, Elizabeth A, Cynthia C. Survival: survival analysis. 2020. https://CRAN.R-project.org/package=survival. Accessed 25 Jun 2020.
Mizuno K. Signaling mechanisms and functional roles of cofilin phosphorylation and dephosphorylation. Cell Signal. 2013;25:457–69.
Liao Q, Li R, Zhou R, Pan Z, Xu L, Ding Y, et al. LIM kinase 1 interacts with myosin-9 and alpha-actinin-4 and promotes colorectal cancer progression. Br J Cancer. 2017;117:563–71.
Chen Y, Chen G, Zhang B, Liu C, Yu Y, Jin Y. miR-27b-3p suppresses cell proliferation, migration and invasion by targeting LIMK1 in colorectal cancer. Int J Clin Exp Pathol. 2017;10:9251–61.
Gorovoy M, Niu J, Bernard O, Profirovic J, Minshall R, Neamu R, et al. LIM kinase 1 coordinates microtubule stability and actin polymerization in human endothelial cells. J Biol Chem. 2005;280:26533–42.
Lee M-H, Kundu JK, Chae J-I, Shim J-H. Targeting ROCK/LIMK/cofilin signaling pathway in cancer. Arch Pharm Res. 2019;42:481–91.
Sumi T, Matsumoto K, Nakamura T. Mitosis-dependent phosphorylation and activation of LIM-kinase 1. Biochem Biophys Res Commun. 2002;290:1315–20.
Wang Y, Kuramitsu Y, Kitagawa T, Baron B, Yoshino S, Maehara S-I, et al. Cofilin-phosphatase slingshot-1L (SSH1L) is over-expressed in pancreatic cancer (PC) and contributes to tumor cell migration. Cancer Lett. 2015;360:171–6.
Wang W-S, Zhong H-J, Xiao D-W, Huang X, Liao L-D, Xie Z-F, et al. The expression of CFL1 and N-WASP in esophageal squamous cell carcinoma and its correlation with clinicopathological features. Dis Esophagus. 2010;23:512–21.
Rangel MP, Antonangelo L, Acencio MMP, Faria CS, de Sá VK, Leão PS, et al. Detection of sputum cofilin-1 as indicator of malignancy. Braz J Med Biol Res. 2018;51:e7138.
Collazo J, Zhu B, Larkin S, Martin SK, Pu H, Horbinski C, et al. Cofilin drives cell-invasive and metastatic responses to TGF-β in prostate cancer. Cancer Res. 2014;74:2362–73.
Wang H, Tao L, Jin F, Gu H, Dai X, Ni T, et al. Cofilin 1 induces the epithelial-mesenchymal transition of gastric cancer cells by promoting cytoskeletal rearrangement. Oncotarget. 2017;8:39131–42.
Shibue T, Brooks MW, Weinberg RA. An integrin-linked machinery of cytoskeletal regulation that enables experimental tumor initiation and metastatic colonization. Cancer Cell. 2013;24:481–98.
De Smedt L, Palmans S, Andel D, Govaere O, Boeckx B, Smeets D, et al. Expression profiling of budding cells in colorectal cancer reveals an EMT-like phenotype and molecular subtype switching. Br J Cancer. 2017;116:58–65.
Ishibashi F. High glucose increases phosphocofilin via phosphorylation of LIM kinase due to Rho/Rho kinase activation in cultured pig proximal tubular epithelial cells. Diabetes Res Clin Pract. 2008;80:24–33.
Palanivel R, Ganguly R, Turdi S, Xu A, Sweeney G. Adiponectin stimulates Rho-mediated actin cytoskeleton remodeling and glucose uptake via APPL1 in primary cardiomyocytes. Metab Clin Exp. 2014;63:1363–73.
Zhang H-S, Zhao J-W, Wang H, Zhang H-Y, Ji Q-Y, Meng L-J, et al. LIM kinase 1 is required for insulin–dependent cell growth of osteosarcoma cell lines. Mol Med Rep. 2014;9:103–8.
The authors are grateful to the Banco Nacional de tumores of Instituto Nacional de Câncer (BNT-INCA) and the Dipat-INCA for provide the CRC specimens, and for their assistance in the RNA extraction and IHC technique. The authors thank also all sample donors. A.C.M.S.S. is the recipient of a postdoctoral fellow from Instituto Nacional de Câncer—Ministério da Saúde.
This research was supported by grants of Conselho Nacional de Desenvolvimento Científico e Tecnológico, grant number 309259/2018-5; Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, grant number E-26/203.325/2017, E-26/010.101072/2018 and Instituto Nacional de Câncer—Ministério da Saúde.
Ethics approval and consent to participate
The analyses involving human participants were reviewed and approved by Ethics Research Committee of Brazilian National Cancer Institute (CEP-INCA), number 84/04 and updated on 13/07/2017. The patients/participants provided their written informed consent to participate in this study.
Consent for publication
There is no competing interesting to declaim.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Analysis of overall survival according to SSH1 expression according to tumor stage. Kaplan–Meier curves depicting the overall survival using CRC patient’s data from TCGA Data Bank were stratified based on SSH1 expression level according to tumor stage (early/late). (a) All tumor stages I–IV (High n=266; Low n=264), (b) early stage I-II (High n=174; Low n=173 and (c) late stage III-IV (High n=137; Low n=136). Samples expressing the gene above or equal the median were classified as high expression and the others as low expression. The P values were derived from the log-rank test, and are indicated.
Analysis of overall survival according to CFL-1 expression according to tumor stage. Kaplan–Meier curves depicting the overall survival using CRC patient’s data from TCGA Data Bank were stratified based on CFL-1 expression level according to tumor stage (early/late). (a) Early stage I-II (High n=174; Low n=173 and (b) late stage III-IV (High n=137; Low n=136). Samples expressing the gene above or equal the median were classified as high expression and the others as low expression. The P values were derived from the log-rank test, and are indicated.
Multivariate analysis for the lymph node metastatic risk.
Analysis of overall survival according to CFL-1 expression in CRC tissues according to CMS classification. Kaplan–Meier curves depicting the overall survival were generated using CRC patient’s data from TCGA Data Bank. (a) CMS1 immune (High n=42; Low n=41), (b) CMS2 canonical (High n=113; Low n=112), (c) CMS3 metabolic (High n=36; Low n=36), and (d) CMS4 mesenchymal (High n=75; Low n=75). The P values were derived from the log-rank test.
Analysis of overall survival according to SSH1 expression in CRC tissues according to CMS classification. Kaplan–Meier curves depicting the overall survival were generated using CRC patient’s data from TCGA Data Bank. (a) CMS1 immune (High n=42; Low n=41), (b) CMS2 canonical (High n=113; Low n=112), (c) CMS3 metabolic (High n=36; Low n=36), and (d) CMS4 mesenchymal (High n=75; Low n=75). The P values were derived from the log-rank test.
Analysis of methylation levels in the promoter region of CFL-1 (a) and SSH1 (b) using the UALCAN approach. Colon adenocarcinoma (COAD) and rectal adenocarcinoma (READ) from the TCGA database were used. *P<0.05; **P<0.01; ***P<0.001 and ****P<0.0001. Not significant (n.s.).
About this article
Cite this article
Sousa-Squiavinato, A.C.M., Vasconcelos, R.I., Gehren, A.S. et al. Cofilin-1, LIMK1 and SSH1 are differentially expressed in locally advanced colorectal cancer and according to consensus molecular subtypes. Cancer Cell Int 21, 69 (2021). https://doi.org/10.1186/s12935-021-01770-w
- Colorectal cancer
- Consensus molecular subtypes
- Lymph node metastasis