Overexpression of the S100A2 protein as a prognostic marker for patients with stage II and III colorectal cancer

Protein expression of S100A2 reveals it association with patient prognosis and immune infiltration profile in colorectal cancer

PURPOSE

Colorectal cancer (CRC) is the third most diagnosed cancer worldwide. Despite a well-established knowledge of tumour development, biomarkers to predict patient outcomes are still required. S100 calcium-binding protein A2 (S100A2) has been purposed as a potential marker in many types of cancer, however, the prognostic value of S100A2 in CRC is rarely reported.

MATERIAL AND METHODS

In this study, immunohistochemistry (IHC) was performed to identify the prognostic role of S100A2 protein expression in the tumour core of the tissue microarrays (TMAs) in colorectal cancer patients (n=787). Bulk RNA transcriptomic data was used to identify significant genes compared between low and high cytoplasmic S100A2 groups. Multiplex immunofluorescence (mIF) was performed to further study and confirm the immune infiltration in tumours with low and high cytoplasmic S100A2.

RESULTS

Low cytoplasmic protein expression of S100A2 in the tumour core was associated with poor survival (HR 0.539, 95%CI 0.394-0.737, P<0.001) and other adverse tumour phenotypes. RNA transcriptomic analysis showed a gene significantly associated with the low cytoplasmic S100A2 group (AKT3, TAGLN, MYLK, FGD6 and ETFDH), which correlated with tumour development and progression. GSEA analysis identifies the enriched anti-tumour and immune activity group of genes in high cytoplasmic S100A2. Additionally, mIF staining showed that high CD3+FOXP3+ and CD163+ inversely associated with low cytoplasmic S100A2 (P<0.001, P=0.009 respectively).

CONCLUSION

Our finding demonstrates a prognostic value of S100A2 together with the correlation with immune infiltration in CRC.

Increased Expression of the RBPMS Splice Variants Inhibits Cell Proliferation in Ovarian Cancer Cells

RNA-Binding Protein with Multiple Splicing (RBPMS) is a member of family proteins that bind to nascent RNA transcripts and regulate their splicing, localization, and stability. Evidence indicates that RBPMS controls the activity of transcription factors associated with cell growth and proliferation, including AP-1 and Smads. Three major RBPMS protein splice variants (RBPMSA, RBPMSB, and RBPMSC) have been described in the literature. We previously reported that reduced RBPMS levels decreased the sensitivity of ovarian cancer cells to cisplatin treatment. However, little is known about the biological role of the RBPMS splice variants in ovarian cancer cells. We performed RT-PCR and Western blots and observed that both RBPMSA and RBPMSC are reduced at the mRNA and protein levels in cisplatin resistant as compared with cisplatin sensitive ovarian cancer cells. The mRNA and protein levels of RBPMSB were not detectable in any of the ovarian cancer cells tested. To better understand the biological role of each RBPMSA and RBPMSC, we transfected these two splice variants in the A2780CP20 and OVCAR3CIS cisplatin resistant ovarian cancer cells and performed cell proliferation, cell migration, and invasion assays. Compared with control clones, a significant reduction in the number of colonies, colony size, cell migration, and invasion was observed with RBPMSA and RBPMSC overexpressed cells. Moreover, A2780CP20-RBPMSA and A2780CP20-RBPMSC clones showed reduced senescence-associated β-galactosidase (β-Gal)-levels when compared with control clones. A2780CP20-RBPMSA clones were more sensitive to cisplatin treatment as compared with A2780CP20-RBPMSC clones. The A2780CP20-RBPMSA and A2780CP20-RBPMSC clones subcutaneously injected into athymic nude mice formed smaller tumors as compared with A2780CP20-EV control group. Additionally, immunohistochemical analysis showed lower proliferation (Ki67) and angiogenesis (CD31) staining in tissue sections of A2780CP20-RBPMSA and A2780CP20-RBPMSC tumors compared with controls. RNAseq studies revealed many common RNA transcripts altered in A2780CP20-RBPMSA and A2780CP20-RBPMSC clones. Unique RNA transcripts deregulated by each RBPMS variant were also observed. Kaplan-Meier (KM) plotter database information identified clinically relevant RBPMSA and RBPMSC downstream effectors. These studies suggest that increased levels of RBPMSA and RBPMSC reduce cell proliferation in ovarian cancer cells. However, only RBPMSA expression levels were associated with the sensitivity of ovarian cancer cells to cisplatin treatment.

Influence of S100A2 in Human Diseases

S100 proteins are a family of low-molecular-weight proteins characterized by two calcium-binding sites with a helix-loop-helix (“EF-hand-type”) domain. The S100 family of proteins is distributed across various organs and can interact with diverse molecules. Among the proteins of the S100 family, S100 calcium-binding protein A2 (S100A2) has been identified in mammary epithelial cells, glands, lungs, kidneys, and prostate gland, exhibiting various physiological and pathological actions in human disorders, such as inflammatory diseases and malignant tumors. In this review, we introduce basic knowledge regarding S100A2 regulatory mechanisms. Although S100A2 is a tumor suppressor, we describe the various influences of S100A2 on cancer and inflammatory diseases.

Transcriptional activation of S100A2 expression by HIF-1α via binding to the hypomethylated hypoxia response elements in HCC cells

Hepatocellular carcinoma (HCC) is one of the most prevalent and lethal cancers. Dysregulation of S100A2 has recently been found in many cancers including HCC. However, its regulatory mechanism in HCC remains poorly understood, especially in hypoxia. In this study, we found that S100A2 is upregulated and correlated with the clinicopathological features of HCC patients. Moreover, the elevated S100A2 showed worse overall survival. Functionally, S100A2 inhibition decreased the proliferation and migration of HepG2 cells. Interestingly, we found that HIF-1α directly binds to hypoxia response elements (HREs) of the S100A2 promoter region. S100A2 expression could be induced in an HIF-1α-dependent manner under hypoxia. Furthermore, S100A2 silencing significantly suppressed HCC cell proliferation and invasion under hypoxia. Mechanistically, pyrosequencing results showed that the hypomethylation status of CpG located in the HRE at the S100A2 promoter was correlated with S100A2 induction. Additionally, HIF-1α- mediated S100A2 activation was associated with TET2-related epigenetic inactivation. TET2 was enriched in the HRE of the S100A2 promoter in HepG2 cells. Finally, S100A2 methylation-related genes and pathways were analyzed. We found that the methylation of S100A2 is correlated with ANXA2, PPP1R15A, and FOS, which include in a hypoxia-related gene set from the GSEA database. Moreover, some EMT-related genes are associated with the methylation of S100A2 in HCC. Conclusively, our study thus uncovered a novel mechanism showing that hypoxia/HIF-1α signaling associated with DNA methylation enhances S100A2 expression in HCC. S100A2 may be useful as a target for facilitating novel diagnostic and therapeutic strategies in liver cancer.

High expression of S100A2 predicts poor prognosis in patients with endometrial carcinoma

Background

S100A2, a member of the S100 protein family, is abnormally expressed and plays a vital role in multiple cancers. However, little is known about the clinical significance of S100A2 in endometrial carcinoma.

Methods

Clinicopathological data were obtained from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), Gene Expression Omnibus (GEO), and Clinical Proteomic Tumor Analysis Consortium (CPTAC). First, the expression and prognostic value of different S100 family members in endometrial carcinoma were evaluated. Subsequently, the Kaplan–Meier plotter and Cox regression analysis were used to assess the prognostic significance of S100A2, while the association between S100A2 expression and clinical characteristics in endometrial carcinoma was also analyzed using logistic regression. A receiver operating characteristic (ROC) curve and a nomogram were constructed. The putative underlying cellular mechanisms were explored using Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and gene set enrichment analysis (GSEA).

Results

Our results revealed that S100A2 expression was significantly higher in endometrial carcinoma tissue than in non-cancerous tissue at both the mRNA and protein levels. Analysis of Kaplan–Meier plotter data revealed that patients with high S100A2 expression had shorter overall survival (OS) and disease specific survival (DSS) compared with those of patients with low S100A2 expression. Multivariate Cox analysis further confirmed that high S100A2 expression was an independent risk factor for OS in patients with endometrial carcinoma. Other clinicopathologic features found to be related to worse prognosis in endometrial carcinoma included age, clinical stage, histologic grade, and tumor invasion. Importantly, ROC analysis also confirmed that S100A2 has a high diagnostic value in endometrial carcinoma. KEGG enrichment analysis and GSEA revealed that the estrogen and IL-17 signaling pathways were significantly upregulated in the high S100A2 expression group, in which estrogen response, JAK-STAT3, K-Ras, and TNFα/NF-κB were differentially enriched.

Conclusions

S100A2 plays an important role in endometrial carcinoma progression and may represent an independent diagnostic and prognostic biomarker for endometrial carcinoma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09180-5.

The interaction between S100A2 and KPNA2 mediates NFYA nuclear import and is a novel therapeutic target for colorectal cancer metastasis

Nucleocytoplasmic transport of proteins is disrupted and dysregulated in cancer cells. Nuclear pore complexes and cargo proteins are two main transportation regulators. However, the mechanism regulating nucleocytoplasmic transport in cancer remains elusive. Here, we identified a S100A2/KPNA2 cotransport complex that transports the tumor-associated transcription factor NFYA in colorectal cancer (CRC). Through the S100A2/KNPA2 complex, depending on its interaction with S100A2, NFYA is transported to the nucleus and inhibits the transcriptional activity of E-cadherin, which in turn promotes CRC metastasis. Targeting the S100A2/KPNA2 binding sites with the specific inhibitor delanzomib is a potential therapeutic approach for CRC.

Common and mutation specific phenotypes of KRAS and BRAF mutations in colorectal cancer cells revealed by integrative -omics analysis

BACKGROUND

Genes in the Ras pathway have somatic mutations in at least 60 % of colorectal cancers. Despite activating the same pathway, the BRAF V600E mutation and the prevalent mutations in codon 12 and 13 of KRAS have all been linked to different clinical outcomes, but the molecular mechanisms behind these differences largely remain to be clarified.

METHODS

To characterize the similarities and differences between common activating KRAS mutations and between KRAS and BRAF mutations, we used genome editing to engineer KRAS G12C/D/V and G13D mutations in colorectal cancer cells that had their mutant BRAF V600E allele removed and subjected them to transcriptome sequencing, global proteomics and metabolomics analyses.

RESULTS

By intersecting differentially expressed genes, proteins and metabolites, we uncovered (i) two-fold more regulated genes and proteins when comparing KRAS to BRAF mutant cells to those lacking Ras pathway mutation, (ii) five differentially expressed proteins in KRAS mutants compared to cells lacking Ras pathway mutation (IFI16, S100A10, CD44, GLRX and AHNAK2) and 6 (CRABP2, FLNA, NXN, LCP1, S100A10 and S100A2) compared to BRAF mutant cells, (iii) 19 proteins expressed differentially in a KRAS mutation specific manner versus BRAF V600E cells, (iv) regulation of the Integrin Linked Kinase pathway by KRAS but not BRAF mutation, (v) regulation of amino acid metabolism, particularly of the tyrosine, histidine, arginine and proline pathways, the urea cycle and purine metabolism by Ras pathway mutations, (vi) increased free carnitine in KRAS and BRAF mutant RKO cells.

CONCLUSIONS

This comprehensive integrative -omics analysis confirms known and adds novel genes, proteins and metabolic pathways regulated by mutant KRAS and BRAF signaling in colorectal cancer. The results from the new model systems presented here can inform future development of diagnostic and therapeutic approaches targeting tumors with KRAS and BRAF mutations.

ELF3 Overexpression as Prognostic Biomarker for Recurrence of Stage II Colorectal Cancer

BACKGROUND/AIM

Adjuvant chemotherapy for high-risk Stage II colorectal cancer (CRC) is weakly recommended; however, no consensus exists on "high-risk" definition. Prognostic biomarker identification is important for selecting patients with poor prognosis who may benefit from adjuvant chemotherapy.

MATERIALS AND METHODS

Using Microarray data analyses, ELF3 was identified as a candidate gene highly expressed in Stage II CRC with distant recurrences. ELF3 mRNA expression in 168 Stage II CRC patients was subjected to quantitative RT-PCR analysis and ELF3 protein expression in 185 patients was quantified by immunohistochemical analysis. The relationship between mRNA and protein expression levels and patient characteristics were also investigated.

RESULTS

The overall recurrence rate and relapse-free survival were significantly poorer in the ELF3 high-expression than the low-expression group at the mRNA and protein levels. High ELF3 mRNA and protein expression levels were independent poor prognostic factors.

CONCLUSION

High ELF3 expression was associated with recurrence of Stage II.

Comprehensive Profiling Reveals Distinct Microenvironment and Metabolism Characterization of Lung Adenocarcinoma

Lung adenocarcinoma has entered into an era of immunotherapy with the development of immune checkpoint inhibitors (ICIs). The identification of immune subtype is crucial to prolonging survival in patients. The tumor microenvironment (TME) and metabolism have a profound impact on prognosis and therapy. The majority of previous studies focused on only one aspect, while both of them are essential to the understanding of tumorigenesis and development. We hypothesized that lung adenocarcinoma can be stratified into immune subgroups with alterations in the TME infiltration. We aimed to explore the "TME-Metabolism-Risk" patterns in each subtypes and the mechanism behind. Glycolysis and cholesterol were selected for the analysis of metabolic states based on the first half of the study. Bioinformatic analysis was performed to investigate the transcriptomic and clinical data integrated by three lung adenocarcinoma cohorts (GSE30219, GSE31210, GSE37745, N = 415). The results were validated in an independent cohort (GSE50081, N = 127). In total, 415 lung adenocarcinoma samples were integrated and analyzed. Four major immune subtypes were indentified using bioinformatic analysis. Subtype NC1, characterized by a high level of glycolysis, with extremely low microenvironment cell infiltration. Subtype NC2, characterized by the "Silence" and "Cholesterol biosynthesis Predominant" metabolic states, with a middle degree infiltration of microenvironment cell. Subtype NC3, characterized by the lack of "Cholesterol biosynthesis Predominant" metabolic state, with abundant microenvironment cell infiltration. Subtype NC4, characterized by "Mixed" metabolic state, with a relatively low microenvironment cell infiltration. Least absolute shrinkage and selection operator (LASSO) regression and multivariate analyses were performed to calculate the risk of each sample, and we attempted to find out the potential immune escape mechanism in different subtypes. The result revealed that the lack of immune cells infiltration might contribute to the immune escape in subtypes NC1 and NC4. NC3 was characterized by the high expression of immune checkpoint molecules and fibroblasts. NC2 had defects in activation of innate immune cells. There existed an obviously survival advantage in subtype NC2. Gene set enrichment analysis (GSEA) and Gene Ontology analysis indicated that the PI3K-AKT-mTOR, TGF-β, MYC-related pathways might be correlated with this phenomenon. In addition, some differentially expressed genes (DEGs) were indentified in subtype NC3, which might be potential targets for survival phenotype transformation.

Integrative analysis identifies an immune-relevant epigenetic signature for prognostication of non-G-CIMP glioblastomas

The clinical and molecular implications of DNA methylation alterations remain unclear among the majority of glioblastomas (GBMs) without glioma-CpGs island methylator phenotype (G-CIMP); integrative multi-level molecular profiling may provide useful information. Independent cohorts of non-G-CIMP GBMs or IDH wild type (wt) lower-grade gliomas (LGGs) from local and public databases with DNA methylation and gene expression microarray data were included for discovery and validation of a multimarker signature, combined using a RISK score model. Bioinformatic and in vitro functional analyses were employed for biological validation. Using a strict multistep selection approach, we identified eight CpGs, each of which was significantly correlated with overall survival (OS) of non-G-CIMP GBMs, independent of age, the O-6-methylguanine-DNA methyltransferase (MGMT) methylation status, treatments and other identified CpGs. An epigenetic RISK signature of the 8 CpGs was developed and validated to robustly and independently prognosticate prognosis in different cohorts of not only non-G-GIMP GBMs, but also IDHwt LGGs. It also showed good discriminating value in stratified cohorts by current clinical and molecular factors. Bioinformatic analysis revealed consistent correlation of the epigenetic signature to distinct immune-relevant transcriptional profiles of GBM bulks. Functional experiments showed that S100A2 appeared to be epigenetically regulated by one identified CpG and was associated with GBM cell proliferation, apoptosis, invasion, migration and immunosuppression. The prognostic 8-CpGs RISK score signature may be of promising value for refining current glioma risk classification, and its potential links to distinct immune phenotypes make it a promising biomarker candidate for predicting response to anti-glioma immunotherapy.

The Importance of Sex in the Discovery of Colorectal Cancer Prognostic Biomarkers

Colorectal cancer (CRC) is the third leading cause of cancer deaths. Advances within bioinformatics, such as machine learning, can improve biomarker discovery and ultimately improve CRC survival rates. There are clear sex differences in CRC characteristics, but the impact of sex has not been considered with regards to CRC biomarkers. Our aim here was to investigate sex differences in the transcriptome of a normal colon and CRC, and between paired normal and tumor tissue. Next, we attempted to identify CRC diagnostic and prognostic biomarkers and investigate if they are sex-specific. We collected paired normal and tumor tissue, performed RNA-seq, and applied feature selection in combination with machine learning to identify the top CRC diagnostic biomarkers. We used The Cancer Genome Atlas (TCGA) data to identify sex-specific CRC diagnostic biomarkers and performed an overall survival analysis to identify sex-specific prognostic biomarkers. We found transcriptomic sex differences in both the normal colon tissue and in CRC. Forty-four of the top-ranked biomarkers were sex-specific and 20 biomarkers showed a sex-specific prognostic value. Our data show the importance of sex in the discovery of CRC biomarkers. We propose 20 sex-specific CRC prognostic biomarkers, including ESM1, GUCA2A, and VWA2 for males and CLDN1 and FUT1 for females.

Tissue mRNA for S100A4, S100A6, S100A8, S100A9, S100A11 and S100P Proteins in Colorectal Neoplasia: A Pilot Study

S100 proteins are involved in the pathogenesis of sporadic colorectal carcinoma through different mechanisms. The aim of our study was to assess tissue mRNA encoding S100 proteins in patients with non-advanced and advanced colorectal adenoma. Mucosal biopsies were taken from the caecum, transverse colon and rectum during diagnostic and/or therapeutic colonoscopy. Another biopsy was obtained from adenomatous tissue in the advanced adenoma group. The tissue mRNA for each S100 protein (S100A4, S100A6, S100A8, S100A9, S100A11 and S100P) was investigated. Eighteen biopsies were obtained from the healthy mucosa in controls and the non-advanced adenoma group (six individuals in each group) and thirty biopsies in the advanced adenoma group (ten patients). Nine biopsies were obtained from advanced adenoma tissue (9/10 patients). Significant differences in mRNA investigated in the healthy mucosa were identified between (1) controls and the advanced adenoma group for S100A6 (p = 0.012), (2) controls and the non-advanced adenoma group for S100A8 (p = 0.033) and (3) controls and the advanced adenoma group for S100A11 (p = 0.005). In the advanced adenoma group, differences between the healthy mucosa and adenomatous tissue were found in S100A6 (p = 0.002), S100A8 (p = 0.002), S100A9 (p = 0.021) and S100A11 (p = 0.029). Abnormal mRNA expression for different S100 proteins was identified in the pathological adenomatous tissue as well as in the morphologically normal large intestinal mucosa.

miR-325-3p Promotes the Proliferation, Invasion, and EMT of Breast Cancer Cells by Directly Targeting S100A2

This study was designed to investigate the precise mechanisms of miR-325-3p/S100A2 axis in breast cancer (BC). In this study, we found that the level of miR-325-3p was dramatically increased in BC tissues and cell lines, and the expression of S100A2 was significantly decreased. Also, the high level of miR-325-3p was closely associated with low expression of S100A2 in BC tissues. Moreover, introduction of miR-325-3p significantly promoted proliferation, invasion, and EMT of BC cells. Bioinformatics analysis predicted that the S100A2 was a potential target gene of miR-325-3p. Luciferase reporter assay demonstrated that miR-325-3p could directly target S100A2. In addition, miR-325-3p overexpression had similar effects with knockdown of S100A2 on BC cells. Overexpression of S100A2 in BC cells partially reversed the promoted effects of miR-325-3p mimic. Overexpression of miR-325-3p promoted cell proliferation, invasion, and EMT of BC cells by directly downregulating S100A2 expression.

S100A16 suppresses the proliferation, migration and invasion of colorectal cancer cells in part via the JNK/p38 MAPK pathway

S100 calcium binding protein A16 (S100A16) is the most recent member of the S100 calcium-binding protein family. The function of S100A16 has been associated with various types of cancer; however, its role in colorectal cancer (CRC) remains unknown. Therefore, the aim of the present study was to investigate the role of S100A16 in CRC progression. The Oncomine dataset used in the current study revealed that the expression of S100A16 was decreased in CRC compared with normal colorectal tissues. Similar results were also determined via immunohistochemistry. In addition, a negative association was identified between S100A16 expression and the prognosis of patients with CRC. Further functional experiments revealed that S100A16 knockdown promoted the proliferation, migration and invasion of HCT116 and SW480 cells, and vice versa in Lovo cells. Epithelial-mesenchymal transition (EMT) was promoted and the JNK/p38 MAPK pathway was activated in HCT116 cells following S100A16 knockdown, as determined via western blotting. Furthermore, S100A16 silencing promoted the migration and invasion of cells. EMT was also reversed when cells were treated with the JNK inhibitor (SP600125) or the p38 inhibitor (SB203580). In summary, the results of the present study demonstrated that S100A16 suppressed the proliferation, migration and invasion of CRC cells partially via the JNK/p38 MAPK signalling pathway and subsequent EMT mediation.

An Integrated Bioinformatic Analysis of the S100 Gene Family for the Prognosis of Colorectal Cancer

Background

S100 family genes exclusively encode at least 20 calcium-binding proteins, which possess a wide spectrum of intracellular and extracellular functions in vertebrates. Multiple lines of evidences suggest that dysregulated S100 proteins are associated with human malignancies including colorectal cancer (CRC). However, the diverse expression patterns and prognostic roles of distinct S100 genes in CRC have not been fully elucidated.

Methods

In the current study, we analyzed the mRNA expression levels of S100 family genes and proteins and their associations with the survival of CRC patients using the Oncomine analysis and GEPIA databases. Expressions and mutations of S100 family genes were analyzed using the cBioPortal, and protein-protein interaction (PPI) networks of S100 proteins and their mutation-related coexpressed genes were analyzed using STRING and Cytoscape.

Results

We observed that the mRNA expression levels of S100A2, S100A3, S100A9, S100A11, and S100P were higher and the level of S100B was lower in CRC tissues than those in normal colon mucosa. A high S100A10 levels was associated with advanced-stage CRC. Results from GEPIA database showed that highly expressed S100A1 was correlated with worse overall survival (OS) and disease-free survival (DFS) and that overexpressions of S100A2 and S100A11 were associated with poor DFS of CRC, indicating that S100A1, S100A2, and S100A11 are potential prognostic markers. Unexpectedly, most of S100 family genes showed no significant prognostic values in CRC.

Conclusions

Our findings, though still need to be ascertained, offer novel insights into the prognostic implications of the S100 family in CRC and will inspire more clinical trials to explore potential S100-targeted inhibitors for the treatment of CRC.

Comprehensive analysis of the expression and prognosis for S100 in human ovarian cancer: A STROBE study

S100 family members are frequently deregulated in human malignancies, including ovarian cancer. However, the prognostic roles of each individual S100 family member in ovarian cancer (OC) patients remain elusive. In the present study, we assessed the prognostic roles and molecular function of 20 individual members of the S100 family in OC patients using GEPIA, Kaplan-Meier plotter, SurvExpress, GeneMANIA and Funrich database. Our results indicated that the mRNA expression levels of S100A1, S100A2, S100A4, S100A5, S100A11, S100A14, and S100A16 were significantly upregulated in patients with OC, and high mRNA expression of S100A1, S100A3, S100A5, S100A6, and S100A13 were significantly correlated with better overall survival, while increased S100A2, S100A7A, S100A10, and S100A11 mRNA expressions were associated with worse prognosis in OC patients. In stratified analysis, the trends of high expression of individual S100 members were nearly the same in different pathological grade, clinical stage, TP53 mutation status, and treatment. More importantly, S100 family signatures may be useful potential prognostic markers for OC. These findings suggest that S100 family plays a vital role in prognostic value and could potentially be an S100-targeted inhibitors for OC patients.

Identifying four DNA methylation gene sites signature for predicting prognosis of osteosarcoma

BACKGROUND

Osteosarcoma (OS) is a common malignant bone tumor in children and adolescents. DNA methylation plays a crucial role in the prognosis prediction of cancer. Identification of novel DNA methylation sites biomarkers could be beneficial for the prognosis of OS patients. In this study, we aim to find an efficient methylated site model for predicting survival in OS.

METHODS

DNA methylation data were downloaded from the Cancer Genome Atlas database (TCGA) and the GEO database. Cox proportional hazard regression and random survival forest algorithm (RSFVH) were applied to identify DNA methylated site signature in the samples randomly assigned to the training subset and the other samples as the test subset. By randomizing 71 clinical samples into two individual groups and a series of statistical analyses between the two groups, a DNA methylation signature is verified.

RESULTS

This signature comprises four methylation sites (cg04533248, cg12401425, cg13997435, and cg15075357) associated with the patient training group from the univariate Cox proportional hazards regression analysis, RSFVH, and multivariate Cox regression analysis. Kaplan-Meier survival curves showed the OS patients in the high-risk group have a poor 5-year overall survival compared with the low-risk group, and this finding was identified in the test data set. A ROC analysis was performed in the current research. The results revealed that this signature was an independent predictor of patient survival by investigating the AUC of the four methylation sites signature in the training data set (AUC =0.861) and test data set, respectively (AUC =0.920). The nomogram described in the current study placed a great guiding value for predicting 1-, 2-, 3-year survival of the OS by combining age, gender, grade, and TNM stage as covariates with the RS of patients' methylation related signatures.

CONCLUSIONS

Our study proved that this signature might be a powerful prognostic tool for survival rate evaluation and guide tailored therapy for OS patients.

S100A2 promotes glycolysis and proliferation via GLUT1 regulation in colorectal cancer

The deregulation of S100A2 has been implicated in the pathogenesis of several types of cancers. However, the molecular mechanisms underlying the protumorigenic capacities of S100A2 have not been fully elucidated. Here, we demonstrated the molecular mechanisms underlying the roles of S100A2 in glycolysis reprogramming and proliferation of colorectal cancer (CRC) cells. The results indicated that S100A2 overexpression raises glucose metabolism and proliferation. Mechanistically, S100A2 activated the PI3K/AKT signaling pathway, upregulated GLUT1 expression, induced glycolytic reprogramming, and consequently increased proliferation. Clinical data showed significantly increased S100A2 levels in CRC tissues and the Oncomine database. In addition, analysis revealed a positive correlation between S100A2 and GLUT1 mRNA expression in CRC tissues. Together, these results demonstrate that the S100A2/GLUT1 axis can promote the progression of CRC by modulating glycolytic reprogramming. Our results further suggest that targeting S100A2 could present a promising therapeutic avenue for the prevention of colorectal cancer progression.

Identification of EMT-related high-risk stage II colorectal cancer and characterisation of metastasis-related genes

Background

Our laboratory previously reported an individual-level prognostic signature for patients with stage II colorectal cancer (CRC). However, this signature was not applicable for RNA-sequencing datasets. In this study, we constructed a robust epithelial-to-mesenchymal transition (EMT)- related gene pair prognostic signature.

Methods

Based on EMT-related genes, metastasis-associated gene pairs were identified between metastatic and non-metastatic samples. Then, we selected prognosis-associated gene pairs, which were significantly correlated with disease-free survival of stage II CRC using multivariate Cox regression model, as the EMT-related prognosis signature.

Results

An EMT-related signature composed of fifty-one gene pairs (51-GPS) for prediction-relapse risk of patients with stage II CRC was developed, whose prognostic efficiency was validated in independent datasets. Moreover, 51-GPS achieved better predictive performance than other reported signatures, including a commercial signature Oncotype Dx colon cancer and an immune-related gene pair signature. Besides, EMT-related functional gene sets achieved high enrichment scores in high-risk samples. Especially, loss-of-function antisense approach showed that DEGs between the predicted two clusters were metastasis-related.

Conclusions

The EMT-related gene pair signature can identify the high relapse-risk patients with stage II CRC, which can facilitate individualised management of patients.

Serum S100A6, S100A8, S100A9 and S100A11 proteins in colorectal neoplasia: results of a single centre prospective study

S100 proteins are involved in biological events related to colorectal carcinogenesis. Aim of this prospective study was to assess serum concentration of S100A6, A8, A9 and A11 proteins in patients with colorectal neoplasia. Eighty-four subjects were enrolled: 20 controls (average risk population with normal findings on colonoscopy; 7 men, 13 women, age 23-74, mean 55 ± 14), 20 patients with non-advanced colorectal adenoma (non-AA, 10 men, 10 women, age 41-82, mean 62 ± 11), 22 with advanced colorectal adenoma (AA, 15 men, 7 women, age 49-80, mean 64 ± 8) and 22 with colorectal cancer (CRC, 12 men, 10 women, age 49-86, mean 69 ± 10). Peripheral venous blood was obtained. Serum S100 proteins were investigated by enzyme immunoassay technique. Serum S100A6 was significantly lower in CRC (mean 8530 ± 4743 ng/L), p = .035 compared to controls (mean 11308 ± 2968 ng/L). Serum S100A8 was significantly higher in AA (median 11955 ng/L, IQR 2681-34756 ng/L), p = .009 and in CRC (median 27532 ng/L, IQR 6794-35092 ng/L), p < .001 compared to controls (median 2513 ng/L, IQR 2111-4881 ng/L). Serum S100A9 concentrations did not differ between any tested group and controls, p > .05. Serum concentration of S100A11 was significantly lower in non-AA (mean 3.5 ± 2.4 μg/L), p = .004 and in CRC (mean 3.4 ± 2.4 μg/L), p = .002 compared to controls (mean 5.9 ± 2.5 μg/L). Sensitivity and specificity for S100A8 protein in patients with CRC were 94% and 73%; positive predictive value 68% and negative predictive value 95%. Patients with colorectal neoplasia have significantly lower serum S100A6 and S100A11 levels, significantly higher S100A8 and unaltered serum S100A9 levels.

Dysregulation of pseudogene/lncRNA-hsa-miR-363-3p-SPOCK2 pathway fuels stage progression of ovarian cancer

Objective: Ovarian cancer is one of the most common and lethal cancer types in women. The molecular mechanism of ovarian cancer progression is still unclear.

Results: Here, we first reported that expression levels of three genes, GJB2, S100A2 and SPOCK2, were significantly higher in advanced stage than that in early stage of ovarian cancer, and upregulation of them indicated poor prognosis of patients with ovarian cancer. Subsequently, 8, 6 and 20 miRNAs were predicted to target GJB2, S100A2 and SPOCK2, respectively. Among these miRNA-mRNA pairs, hsa-miR-363-3p-SPOCK2 axis was the most potential in suppressing progression of ovarian cancer. Mechanistically, we found that hsa-miR-363-3p-SPOCK2 axis was involved in regulation of actin cytoskeleton. Moreover, 6 pseudogenes and 8 lncRNAs were identified to potentially inhibit hsa-miR-363-3p-SPOCK2 axis in ovarian cancer.

Conclusions: Collectively, we elucidate a regulatory role of pseudogene/lncRNA-hsa-miR-363-3p-SPOCK2 pathway in progression of ovarian cancer, which may provide effective therapeutic approaches and promising prognostic biomarkers for ovarian cancer.

Materials and methods: Differentially expressed genes (DEGs) in ovarian cancer were first screened using GSE12470, after which DEGs expression were validated using GEPIA. Kaplan-Meier analysis was employed to assess the prognostic values. Potential miRNAs were predicted by seven target prediction databases, and upstream lncRNAs and pseudogenes of hsa-miR-363-3p were forecasted through miRNet or starBase. UALCAN and starBase were used to obtain the co-expressed genes of SPOCK. Enrichment analysis for these co-expressed genes was performed by Enrichr.

A network-based predictive gene expression signature for recurrence risks in stage II colorectal cancer

The current criteria for defining the recurrence risks of stage II colorectal cancer (CRC) are not robust; therefore, we aimed to explore novel gene signatures to predict recurrence risks and to reveal the underlying mechanisms of stage II CRC. First, the gene expression profiles of 124 patients with stage II CRC from The Cancer Genome Atlas (TCGA) database were obtained to screen differentially expressed genes (DEGs). A total of 202 DEGs, including 128 upregulated and 74 downregulated, were identified in the recurrence group (n = 24) compared to the nonrecurrence group (n = 100). Furthermore, the top 5 DEGs (ZNF561, WFS1, SLC2A1, MFI2, and PTGR1) were identified by random forest variable hunting, and four (ZNF561, WFS1, SLC2A1, and PTGR1) were selected to create a four‐gene recurrent model (GRM), with an area under the curve (AUC) of 0.882 according to the receiver operating characteristic curve, and the robust diagnostic effectiveness of the GRM was further validated with another gene expression profiling dataset (GSE12032), with an AUC of 0.943. The diagnostic effectiveness of the GRM regarding recurrence was associated with poor disease‐free survival in all stages of CRC. In addition, gene ontology functional annotation and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses revealed 18 enriched functions and 6 enriched pathways. Four genes, ABCG2, CACNA1F, CYP19A1, and TF, were identified as hub genes by the protein‐protein interaction network, which further validated that these genes were correlated with a poor pathologic stage and overall survival in all stages of CRC. In conclusion, the GRM can effectively classify stage II CRC into groups of high and low risks of recurrence, thereby making up for the prognostic value of the traditional clinicopathological risk factors defined by the National Comprehensive Cancer Network guidelines. The hub genes may be useful therapeutic targets for recurrence. Thus, the GRM and hub genes could offer clinical value in directing individualized and precision therapeutic regimens for stage II CRC patients.

Expression Profiling Reveals Involvement of WNT Pathway in the Malignant Progression of Sessile Serrated Adenomas

Approximately 15% to 20% of colorectal cancers are developed through the serrated pathway of tumorigenesis, which is associated with BRAF mutation, CpG island methylation phenotype, and MLH1 methylation. However, the detailed process of progression from sessile serrated adenoma (SSA) to dysplasia and carcinoma has not been elucidated. To further characterize mechanisms involved in the dysplastic progression of SSA, we investigated differential expressions of mRNAs between areas with and without dysplasia within the same SSA polyps. Significantly dysregulated genes in paired samples were applied for functional annotation and biological significance. The same lysates from a subset of matched samples were subjected for miRNA expression profiling. Differentially expressed miRNAs were determined, and their targeted mRNAs were compared in parallel to the list of differentially expressed mRNAs from an RNA sequencing study. Fourteen common mRNA targets were identified, which include AXIN2, a known indicator of WNT/β-catenin pathway activation. Together, in this study, different genes, pathways, and biological processes involved in the initiation and progression of dysplasia in the serrated pathway are documented. One of the most significant findings is the involvement of the WNT/β-catenin pathway in the dysplastic progression of SSAs with different genes being targeted in early versus advanced dysplasia.

Prognostic values of S100 family members in ovarian cancer patients

Objective

Exhibiting high consistence in sequence and structure, S100 family members are interchangeable in function and they show a wide spectrum of biological processes, including proliferation, apoptosis, migration, inflammation and differentiation and the like. While the prognostic value of each individual S100 in ovarian cancer is still elusive. In current study, we investigated the prognostic value of S100 family members in the ovarian cancer.

Methods

We used the Kaplan Meier plotter (KM plotter) database, in which updated gene expression data and survival information are from 1657 ovarian cancer patients, to assess the relevance of individual S100 family mRNA expression to overall survival in various ovarian cancer subtypes and different clinicopathological features.

Results

It was found that high expression of S100A2 (HR = 1.18, 95%CI: 1.04–1.34, P = 0.012), S100A7A (HR = 1.3, 95%CI: 1.04–1.63, P = 0.02),S100A10 (HR = 1.2, 95%CI: 1.05–1.38, P = 0.0087),and S100A16 (HR = 1.23, 95%CI: 1–1.51, P = 0.052) were significantly correlated with worse OS in all ovarian cancer patients, while the expression of S100A1 (HR = 0.87, 95%CI: 0.77–0.99, P = 0.039), S100A3 (HR = 0.83, 95%CI: 0.71–0.96, P = 0.0011), S100A5 (HR = 0.84, 95%CI: 0.73–0.97, P = 0.017), S100A6 (HR = 0.84, 95%CI: 0.72–0.98, P = 0.024), S100A13 (HR = 0.85, 95%CI:0.75–0.97, P = 0.014) and S100G (HR = 0.86, 95%CI: 0.74–0.99, P = 0.041) were associated with better prognosis. Furthermore, we assessed the prognostic value of S100 expression in different subtypes and the clinicopathological features, including pathological grades, clinical stages and TP53 mutation status, of ovarian cancer patients.

Conclusion

Comprehensive understanding of the S100 family members may have guiding significance for the diagnosis and outcome of ovarian cancer patients.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-5170-3) contains supplementary material, which is available to authorized users.

Multiple S100 protein isoforms and C-terminal phosphorylation contribute to the paralog-selective regulation of nonmuscle myosin 2 filaments

Nonmuscle myosin 2 (NM2) has three paralogs in mammals, NM2A, NM2B, and NM2C, which have both unique and overlapping functions in cell migration, formation of cell-cell adhesions, and cell polarity. Their assembly into homo- and heterotypic bipolar filaments in living cells is primarily regulated by phosphorylation of the N-terminally bound regulatory light chain. Here, we present evidence that the equilibrium between these filaments and single NM2A and NM2B molecules can be controlled via S100 calcium-binding protein interactions and phosphorylation at the C-terminal end of the heavy chains. Furthermore, we show that in addition to S100A4, other members of the S100 family can also mediate disassembly of homotypic NM2A filaments. Importantly, these proteins can selectively remove NM2A molecules from heterotypic filaments. We also found that tail phosphorylation (at Ser-1956 and Ser-1975) of NM2B by casein kinase 2, as well as phosphomimetic substitutions at sites targeted by protein kinase C (PKC) and transient receptor potential cation channel subfamily M member 7 (TRPM7), down-regulates filament assembly in an additive fashion. Tail phosphorylation of NM2A had a comparatively minor effect on filament stability. S100 binding and tail phosphorylation therefore preferentially disassemble NM2A and NM2B, respectively. These two distinct mechanisms are likely to contribute to the temporal and spatial sorting of the two NM2 paralogs within heterotypic filaments. The existence of multiple NM2A-depolymerizing S100 paralogs offers the potential for diverse regulatory inputs modulating NM2A filament disassembly in cells and provides functional redundancy under both physiological and pathological conditions.

A multiomics analysis of S100 protein family in breast cancer

The S100 gene family is the largest subfamily of calcium binding proteins of EF-hand type, expressed in tissue and cell-specific manner, acting both as intracellular regulators and extracellular mediators. There is a growing interest in the S100 proteins and their relationships with different cancers because of their involvement in a variety of biological events closely related to tumorigenesis and cancer progression. However, the collective role and the possible coordination of this group of proteins, as well as the functional implications of their expression in breast cancer (BC) is still poorly known. We previously reported a large-scale proteomic investigation performed on BC patients for the screening of multiple forms of S100 proteins. Present study was aimed to assess the functional correlation between protein and gene expression patterns and the prognostic values of the S100 family members in BC. By using data mining, we showed that S100 members were collectively deregulated in BC, and their elevated expression levels were correlated with shorter survival and more aggressive phenotypes of BC (basal like, HER2 enriched, ER-negative and high grading). Moreover a multi-omics functional network analysis highlighted the regulatory effects of S100 members on several cellular pathways associated with cancer and cancer progression, expecially immune response and inflammation. Interestingly, for the first time, a pathway analysis was successfully applied on different omics data (transcriptomics and proteomics) revealing a good convergence between pathways affected by S100 in BC. Our data confirm S100 members as a promising panel of biomarkers for BC prognosis.

Regulatory activity based risk model identifies survival of stage II and III colorectal carcinoma

Clinical and pathological indicators are inadequate for prognosis of stage II and III colorectal carcinoma (CRC). In this study, we utilized the activity of regulatory factors, univariate Cox regression and random forest for variable selection and developed a multivariate Cox model to predict the overall survival of Stage II/III colorectal carcinoma in GSE39582 datasets (469 samples). Patients in low-risk group showed a significant longer overall survival and recurrence-free survival time than those in high-risk group. This finding was further validated in five other independent datasets (GSE14333, GSE17536, GSE17537, GSE33113, and GSE37892). Besides, associations between clinicopathological information and risk score were analyzed. A nomogram including risk score was plotted to facilitate the utilization of risk score. The risk score model is also demonstrated to be effective on predicting both overall and recurrence-free survival of chemotherapy received patients. After performing Gene Set Enrichment Analysis (GSEA) between high and low risk groups, we found that several cell-cell interaction KEGG pathways were identified. Funnel plot results showed that there was no publication bias in these datasets. In summary, by utilizing the regulatory activity in stage II and III colorectal carcinoma, the risk score successfully predicts the survival of 1021 stage II/III CRC patients in six independent datasets.

Current Tissue Molecular Markers in Colorectal Cancer: A Literature Review

BACKGROUND

Colorectal cancer (CRC) is one of the most spread neoplasia types all around the world, especially in western areas. It evolves from precancerous lesions and adenomatous polyps, through successive genetic and epigenetic mutations. Numerous risk factors intervene in its development and they are either environmental or genetic.

AIM OF THE REVIEW

Alongside common screening techniques, such as fecal screening tests, endoscopic evaluation, and CT-colonography, we have identified the most important and useful biomarkers and we have analyzed their role in the diagnosis, prevention, and prognosis of CRC.

CONCLUSION

Biomarkers can become an important tool in the diagnostic and therapeutic process for CRC. But further studies are needed to identify a noninvasive, cost-effective, and highly sensible and specific screening test for their detection and to standardize their use in clinical practice.

Standardising RNA profiling based biomarker application in cancer-The need for robust control of technical variables

Histopathology-based staging of colorectal cancer (CRC) has utility in assessing the prognosis of patient subtypes, but as yet cannot accurately predict individual patient's treatment response. Transcriptomics approaches, using array based or next generation sequencing (NGS) platforms, of formalin fixed paraffin embedded tissue can be harnessed to develop multi-gene biomarkers for predicting both prognosis and treatment response, leading to stratification of treatment. While transcriptomics can shape future biomarker development, currently <1% of published biomarkers become clinically validated tests, often due to poor study design or lack of independent validation. In this review of a large number of CRC transcriptional studies, we identify recurrent sources of technical variability that encompass collection, preservation and storage of malignant tissue, nucleic acid extraction, methods to quantitate RNA transcripts and data analysis pipelines. We propose a series of defined steps for removal of these confounding issues, to ultimately aid in the development of more robust clinical biomarkers.

Predictive and prognostic biomarkers in colorectal cancer: A systematic review of recent advances and challenges

BACKGROUND

Colorectal cancer (CRC) is one of the leading cause of cancer deaths worldwide. Since CRC is largely asymptomatic until alarm features develop to advanced stages, the implementation of the screening programme is very much essential to reduce cancer incidence and mortality rates. CRC occurs predominantly from accumulation of genetic and epigenetic changes in colon epithelial cells, which later gets transformed into adenocarcinomas.

SCOPE OF REVIEW

The current challenges of screening paradigm and diagnostic ranges are from semi-invasive methods like colonoscopy to non-invasive stool-based test, have resulted in over-diagnosis and over-treatment of CRC. Hence, new screening initiatives and deep studies are required for early diagnosis of CRC. In this regard, we not only summarise current predictive and prognostic biomarkers with their potential for diagnostic and therapeutic applications, but also describe current limitations, future perspectives and challenges associated with the progression of CRC.

MAJOR CONCLUSIONS

Currently many potential biomarkers have already been successfully translated into clinical practice eg. Fecal haemoglobin, Carcinoembryonic antigen (CEA) and CA19.9, although these are not highly promising diagnostic target for personalized medicine. So there is a critical need for reliable, minimally invasive, highly sensitive and specific genetic markers of an individualised and optimised patient treatment at the earliest disease stage possible.

GENERAL SIGNIFICANCE

Identification of a new biomarker, or a set of biomarkers to the development of a valid, and clinical sensible assay that can be served as an alternative tool for early diagnosis of CRC and open up promising new targets in therapeutic intervention strategies.

Large-Scale Analysis of Gene Expression Data Reveals a Novel Gene Expression Signature Associated with Colorectal Cancer Distant Recurrence

Colorectal cancer (CRC) is the fourth-ranked cause of cancer-related deaths worldwide. Despite recent advances in CRC management, distant recurrence (DR) remains the major cause of mortality in patients with preoperative chemotherapy and radiotherapy, underscoring a need to precisely identify novel gene signatures for predicting the risk of systemic relapse. Herein, we integrated two independent CRC gene expression datasets: the GSE71222 dataset, including 26 patients who developed DR and 126 patients who did not develop DR, and the GSE21510 dataset, including 23 patients who developed DR and 76 patients who did not develop DR. Our data revealed 37 common upregulated genes (fold change (FC) ≥ 1.5, P < 0.05) and three common downregulated genes (FC ≤ 1.5, P < 0.05) between DR and non-recurrent patients from the two datasets. We subsequently validated the upregulated gene panel in the Cancer Genome Atlas CRC datasets (379 patients), which identified a five-gene signature (S100A2, VIP, HOXC6, DACT1, KIF26B) associated with poor overall survival (OS, log-rank test P-value: 1.19 × 10⁻⁴) and poor disease-free survival (DFS, log-rank test P-value: 0.002). In a Cox proportional hazards multiple regression model, the five-gene signature and tumor stage retained their significance as independent prognostic factors for CRC DFS and OS. Therefore, our data identified a novel DR gene expression signature associated with worse prognosis in CRC.