- Primary research
- Open Access
Overexpression of DDIT4 and TPTEP1 are associated with metastasis and advanced stages in colorectal cancer patients: a study utilizing bioinformatics prediction and experimental validation
Cancer Cell International volume 21, Article number: 303 (2021)
Various diagnostic and prognostic tools exist in colorectal cancer (CRC) due to multiple genetic and epigenetic alterations causing the disease. Today, the expression of RNAs is being used as prognostic markers for cancer.
In the current study, various dysregulated RNAs in CRC were identified via bioinformatics prediction. Expression of several of these RNAs were measured by RT-qPCR in 48 tissues from CRC patients as well as in colorectal cancer stem cell-enriched spheroids derived from the HT-29 cell line. The relationships between the expression levels of these RNAs and clinicopathological features were analyzed.
Our bioinformatics analysis determined 11 key mRNAs, 9 hub miRNAs, and 18 lncRNAs which among them 2 coding RNA genes including DDIT4 and SULF1 as well as 3 non-coding RNA genes including TPTEP1, miR-181d-5p, and miR-148b-3p were selected for the further investigations. Expression of DDIT4, TPTEP1, and miR-181d-5p showed significantly increased levels while SULF1 and miR-148b-3p showed decreased levels in CRC tissues compared to the adjacent normal tissues. Positive relationships between DDIT4, SULF1, and TPTEP1 expression and metastasis and advanced stages of CRC were observed. Additionally, our results showed significant correlations between expression of TPTEP1 with DDIT4 and SULF1.
Our findings demonstrated increased expression levels of DDIT4 and TPTEP1 in CRC were associated with more aggressive tumor behavior and more advanced stages of the disease. The positive correlations between TPTEP1 as non-coding RNA and both DDIT4 and SULF1 suggest a regulatory effect of TPTEP1 on these genes.
Colorectal cancer (CRC) is the second most common cancer and leading cause of cancer-related deaths in the world . CRC is now known to be a heterogeneous disease due to the various genetic and epigenetic alterations causing the disease . The existence of a subset of cancer cells named cancer stem cells (CSCs) also leads to tumor heterogeneity by utilizing self-renewal and multi-lineage differentiation features in the tumor . These alterations and CSCs play important roles in development and progression of CRC [4, 5]. CRC is typically classified according to the pathological and clinical features of the American Joint Committee on Cancer (AJCC) and the staging system is used to evaluate prognosis and guide treatment strategies [6, 7]. There are some genetic biomarkers which can aid in estimating prognosis and in guiding treatment selection in CRC patients such as 18q loss of heterozygosity (LOH), p27 Kip1, DNA microsatellite instability , K-RAS mutation  and RNA expression profile [9,10,11]. It is important to find sensitive and specific biomarkers to best guide early and appropriate treatment before disease progression .
Bioinformatics can serve as a very useful tool to investigate the complexity of big datasets, discover novel biomarkers and analyze their validation in clinical studies .
Nowadays, some RNA expression panels are used in clinical cancer such as PAM50 . Although the main focus is on transcripts of coding RNA genes, there are some evidence that non-coding RNAs (ncRNAs) are also involved in hallmarks and pathological processes of cancer [15, 16]. The recent discovery from whole genomes sequencing has revealed that 98% of the human transcriptome contain ncRNAs . Evidence shows that the biological functions of many ncRNAs that are involved in the diseases are unknown. The biology of microRNAs (miRNAs) as the abundant small ncRNAs has been better understood [18, 19]. They can interfere in tumorigenesis by regulating oncogenes and tumor suppressor genes . Small ncRNAs that regulate mRNAs can be predicted by numerous in-silico computational programs [21, 22]. Long non-coding RNAs (lncRNAs) are another type of ncRNAs which are expressed in tissue-specific pattern and dysregulated in cancer  and play important functions in cellular processes such as cell proliferation, motility, and apoptosis . Some reports have demonstrated that the levels of some lncRNAs, miRNAs and mRNAs are controlled and regulated by each other in cancer . Identification of the interacting target RNAs of each lncRNA is an important step in understanding lncRNA functions which can be done through computational prediction of lncRNA–RNA interactions . In the current study, by getting help from bioinformatics analysis and computational algorithms, we selected several genes for further investigation of RNA levels in our CRC patients. These genes included DNA-damage-inducible transcript 4 (DDIT4), sulfatase 1 (SULF1) as coding RNA genes and miR-181d-5p, miR-148b-3p and TPTE Pseudogene 1 (TPTEP1) as ncRNA genes.
DDIT4 also known as REDD1 or RTP801, is expressed in response to diverse stress conditions and its abnormal expression is linked to cancer via the effects on PI3K/Akt/mTOR signaling [27, 28]. In-Silico evaluation has shown dysregulation in RNA expression levels of DDIT4 in several cancers which may be used as a poor prognostic factor in colon cancer .
SULF1 is a sulfatase that selectively remove 6-O-sulfate groups from heparan sulfate (HS). Alternation of HS chains is important in signaling events because heparan sulfate proteoglycans (HSPGs) are released into the extracellular matrix and act as co-receptors which contributes to regulation of cellular processes . Some studies have reported dysregulation of SULF1 expression in CRC [31,32,33].
It was described in a meta-analysis report that dysregulation of miRNA-181d family membrane can be used as prognostic marker in different cancers . Also, dysregulation of miRNA-148b-3p expression was reported in numerous cancers including breast , thyroid , prostate , colorectal  and gastric cancer . Until now, dysregulated expression of TPTEP1 has been reported mainly in patients with human lung  and liver cancer .
In this study, we explored effector networks of mRNAs, miRNAs, and lncRNAs in CRC based on predicted relationships of these RNAs via bioinformatics tools. DDIT4, SULF1, miR-181d-5p, miR-148b-3p and TPTEP1 were selected as the potential biomarkers in CRC patients and their expression levels were measured by RT-qPCR. Also, we investigated the association between these RNA expression levels and clinicopathological features in CRC tissue samples. Although some of these RNAs were reviewed in CRC in the past, there is no data about RNA expression levels of DDIT4 and TPTEP1 and their clinical significance in CRC patients as well as in the colorectal CSC-enriched spheroids. Therefore, based on our knowledge, our study is the first to report these data, and also to explore the correlations of these RNA expression levels amongst each based on our prediction analysis via bioinformatics tools.
Bioinformatics prediction study
Data sources and network construction
In our previous study, we detected differentially expressed genes (DEGs) in total 231 CRC patients obtained from merged five data series on Gene Expression Omnibus (GEO), including GSE41011, GSE62932, GSE63624, GSE77953, and GSE78248. Up-regulated genes with score > 3 were included in the current study, the mean score of up-regulated genes was used as cut-off criteria (Additional file 1: Table S1). The score was obtained from the merging series mentioned above as we described previously . In the following step, the genes which were highly associated with carcinogenesis and colorectal cancer diseases (P < 0.0001) were screened among the up-regulated genes according to the DisGeNET library  on Enrichr . Protein–protein interaction (PPI) network was found using STRING database with stringApp (confidence score > 0.4)  in Cytoscape software . K-means algorithm was used for clustering of the STRING database, and the genes in the largest cluster were selected as the entry criteria for subsequent analysis. Workflow of bioinformatics analysis steps were descripted in Fig. 1.
Network enrichment analysis
The pathway enrichment analysis and Gene Ontology (GO) were done using Enrichr in order to better understanding the biological process and functions of genes. Enrichr is a powerful enrichment analysis online tool which is linked to mammalian gene sets libraries and pathway databases . We used KEGG , Reactome , BioPlanet , and WikiPathways  which are important databases and store a lot of data on biological pathways, for our pathway analysis on Enrichr. The key genes for the present study were selected according to pathway and GO analysis. To visualize results of pathway and GO analysis for key genes, ClueGO plug-in by Cytoscape software was used .
mRNA–miRNA network and lncRNA predictions
Prediction of MicroRNAs was performed for key genes using Target scan , DIANA , miR2Disease , miRWalk , miRNet , and microRNA.org . mRNA–miRNA bipartite network was constructed for experimented and high-quality predicted mRNA–miRNA data. The highest degree and most common miRNAs for key genes were selected as hub miRNAs of mRNA–miRNA bipartite network. Finally, in order to investigate the communication of the miRNAs with other RNAs, lncRNA predictions were performed through the hub miRNAs on the lncRNA SNP2 , miRwalk , LNCeDB  databases.
Selecting genes process among the key genes for experimental study
Enrichment analysis and literature review led to selecting two mRNAs and miRNAs among numerous key genes and hub miRNAs related to CSCs for our experimental study. Also, we selected one of topmost lncRNAs that has putative target site interaction found via computational prediction of lncRNA–mRNA on http://rtools.cbrc.jp/cgi-bin/RNARNA/index.pl  based on RactIP  and IntaRNA databases .
Experimental studies in CRC tissues and colorectal CSC-enriched spheroids
Tissue specimens and clinical data collection
Forty-eight fresh tissue samples (tumor and adjacent normal tissues as control) from the patients with CRC, who have not received any preoperative radiotherapy and other antitumor therapy, were harvested through the surgery at the Firoozgar and Bahman hospitals (Tehran, Iran) between April 2017 and May 2018. All samples were transferred to RNA later (EURx, Poland) immediately after resection and placed into prepared cryogenic vials, and frozen in liquid nitrogen to avoid RNA degradation. The diagnosis of CRC was made by postoperative pathological examination according to the diagnostic criteria from the AJCC . Clinicopathological data from the patients were collected from their electronic medical record system. Clinicopathological features for tumor included: tumor size, vascular invasion, perineural invasion, TNM stage, metastasis, and histologic grade (tumor differentiation) in addition to the sex and age of the patients.
Isolation and confirmation of colorectal CSC-enriched spheroids
Colorectal CSC-enriched spheroids were isolated from the HT-29 cancer cell line which were purchased from Iranian Biological Resource Center (IBRC; C10097, RRID: CVCL_0320). The HT-29 cancer cells were cultured in DMEM/High glucose medium (Gibco, Germany) supplemented with 10% fetal bovine serum (Gibco, Germany), 1% non‐essential amino acids (Gibco, Germany), 2 mM l-glutamine (Gibco, Germany) and 1% Penicillin–Streptomycin (Biowest, France) at 37 °C in 5% CO2 and 95% humidified incubator. Colorectal CSC-enriched spheroids were generated using hanging droplet technique of HT-29 cancer cells descripted earlier . Spheroid culture medium included DMEM/F12 (Gibco, Germany), supplemented with 2% B27 (Gibco, Germany), 10 ng/ml of basic fibroblast growth factor (bFGF), 20 ng/ml epidermal growth factor (EGF) (PeproTech, USA), 1% nonessential amino acids, 2 mM l-glutamine, and 1% Penicillin–Streptomycin. In brief, 70%‐80% confluent HT‐29 cells were detached with 0.05% trypsin/EDTA (Gibco, Germany) and were washed twice with PBS and serum‐free medium. Then, 25 μL droplets containing 10,000 cells re-suspended in spheroid culture medium on the lid of Petri dishes containing 5 ml PBS at 37 °C incubator with 5% CO2 and 95% humidified incubator for 96 h. Drops containing spheroids were harvested by washing with a gentle shaking of media and transferred onto nonattachment flasks (coated flasks with 1.2% poly-HEMA (Sigma, USA)) for 6 days in spheroid culture medium conditions as described above. After observation of sphere morphology by microscope, spheroids were evaluated for CSC features based on RNA expression of stemness genes (OCT4, SOX2, C-MYC, KLF4 and NANOG), ABC transporter genes (ABCB1, ABCG2, and ABCC1) and epithelial-mesenchymal transition (EMT) genes (TWIST1, SNAIL1, Vimentin, and ZEB1) .
RNA extraction and cDNA synthesis
Total RNA was extracted from frozen tissues and cells (HT-29 cell line and colorectal CSC-enriched spheroids), using miRNeasy mini kit (QIAGEN GmbH-Germany) according to the manufacturer’s instructions. RNA samples were separated by agarose gel electrophoresis and their concentration was measured by optical absorbance at 260/280 nm. Complementary DNA (cDNA) was synthesized from extracted RNA using cDNA synthesis kit (TaKaRa Bio, Shiga, Japan) and miRNA cDNA synthesis kit (Bon Yakhteh, Iran).
Real time-quantitative polymerase chain reaction (RT‐qPCR)
The specific primers for amplification with RT-qPCR were designed using Primer-BLAST  and OligoAnalyzer 3.1 software (Integrated DNA Technologies) (Table 1). RT-qPCR was performed to find the expression levels of selected genes from bioinformatics analysis and stemness, ABC transporter and EMT genes that were used for validating colorectal CSC-enriched spheroids. RT-qPCR reactions were performed by SYBR Green PCR Master Mix (Takara, Japan) on Real-Time PCR System (Rotor-Gene Q MDx, Germany). The expression of miRNA genes and TPTEP1 was normalized to internal control of kit and RNU6 (U6) expression levels, respectively. To normalize other mRNAs expression, GAPDH gene was used as an internal control gene. The relative expression levels of the genes were calculated by 2−ΔΔCt method .
Statistical analysis was performed using SPSS 21.0 software (SPSS Inc, Chicago, IL). All data in statistical analyses were expressed as median of RNA expression levels. Significance differences in expression levels of candidate genes between tumor and adjacent normal tissue samples as well as between colorectal CSC-enriched spheroids and HT-29 cell line were analyzed using nonparametric test (Mann–Whitney U test). For comparisons of quantitative values between more than two groups, Kruskal–Wallis test was used. The Spearman's test was applied to evaluate the association between expression levels of these RNAs amongst each other and clinicopathological features. P-value less than 0.05 was considered as statistically significant. GraphPad Prism version 8 software (GraphPad Software, La Jolla, CA) was used for making the boxplots, heat map graph and scatterplots.
Bioinformatics analysis and selecting target genes
Network analysis and clustering genes
Three hundred and seventy up-regulated genes were included in network analysis based on score > 3 which found to be involved in carcinogenesis or colorectal cancer (P < 0.0001) on DisGeNET (Additional file 1: Table S1). PPI network analysis explored the interactions of these up-regulated genes amongst each other. Five main clusters were obtained from k-means algorithm for genes with confidence ≥ 0.4 (Additional file 2: Figure S1). In order to limit the number of genes, largest cluster covering 167 genes (Additional file 1: Table S1) were selected for subsequent analysis in which their PPI network is shown in Fig. 2.
Pathway and GO enrichment analysis
To find better characteristics of the 167 genes, pathway and GO enrichment analysis were performed using the Enrichr tool. The top 10 results of pathway and GO annotation analysis (P < 0.05) were shown in the Additional file 2: Figure S2 and S3. Enrichment analysis displayed the “spliceosome”, “miRNA biogenesis”, “P53 signaling pathway”,” DNA repair”, “MAPK signaling pathway” and “gene expression” are parts of the top 10 pathways.
A closer check of GO and pathway analysis indicated that some of the genes participate in “microRNAs in cancer”, “proteoglycans in cancer”, “apoptosis” and “cell cycle” pathways. These genes contribute to several key biological processes including “extracellular matrix organization”, “regulation of cell migration” and “positive regulation of cell proliferation” based on GO analysis which their disorder was reported in cancer . Information of functional characteristics of genes led to restricting these genes to 11 key genes (TWIST1, DDIT4, LAMC2, SULF1, REG1A, REG3A, VSNL1, BNIP3, GPSM2, GTF3A and SMNDC1). The information pathways and common features of 11 key genes are summarized in Fig. 3.
mRNA–miRNA network and prediction of lncRNAs
One mRNA-miRNA bipartite network was created from key genes and miRNAs related to them (Fig. 4 and Additional file 3: Table S2). Nine miRNAs with the highest degree and most common for 11 key genes (hub miRNAs) were selected as most effective miRNAs on the network of CRC (hsa-miR-1, hsa-miR-125a-5p, hsa-miR-129-5p, hsa-miR-1297, hsa-miR-137, hsa-miR-145, hsa-miR-148b, hsa-miR-181d and hsa-miR-185). To reduce analysis complexity, lncRNAs for hub miRNAs were predicted using miRNA–lncRNA target algorithms in several databases (Additional file 4: Table S3). The Topmost of lncRNAs (18 lncRNAs), those experimentally have been supported are summarized in Table 2.
Finally, reviewing the literature led to selecting genes of DDIT4, SULF1, miR-181d-5p and miR-148b-3p that are involved in CSCs amongst the key genes and hub miRNAs for our experimental validation [68,69,70,71]. Also, TPTEP1 was selected on topmost of the predicted lncRNAs for our experimental study that has target and interaction sites in untranslated region (UTR) and coding sequence (CDS) region for DDIT4 and SULF1, respectively, as found by prediction algorithms (Fig. 5).
Experimental studies in CRC tissues and colorectal CSC-enriched spheroids
Out of 48 patients with CRC, the number of males and females were 29 (60.4%) and 19 (39.6%) respectively. The patients were in the age group of 20–87 years with the mean age of 59 ± 13.7 (mean ± SD) years. Twenty-five (52.1%) of the cases were in early stages (I-II) while 23 (47.9%) were in advanced stages (III-IV) of the tumor. Twenty-one patients (43.8%) had well differentiated cancer cells (Grade 1), followed by 19 cases (39.6%) with moderately differentiated (Grade 2) and 8 cases (16.6%) with poorly differentiated cancer cells (Grade 3). Patients were classified into 3 groups: 25 patients (52.1%) didn’t have any metastasis, 14 patients (29.1%) had only lymph node metastasis and 9 patients (18.7%) had distant metastasis. Vascular and perineural invasion were observed in 14 (29.1%) and 18 (37.5%) patients respectively.
RNA expression levels of selected genes in CRC tissues and the relationship with clinicopathological features
The mRNAs and miRNAs expression levels were evaluated in 48 tumor tissues and their adjacent normal tissues of CRC patients by RT-qPCR. The analysis of the RT-qPCR data using Mann–Whitney U test demonstrated median expression levels of DDIT4 (P = 0.007), TPTEP1 (P = 0.035) and miR-181d-5p (P = 0.020) were significantly higher in CRC tissues compared to the adjacent normal tissues (Table 3) (frame A, C and D of Fig. 6, respectively). In contrast, the expression levels of SULF1 (P = 0.032) and miRNA-148b-3p (P < 0. 001) were significantly lower in CRC tissues compared to the adjacent normal tissues (Table 3) (frame B and E of Fig. 6, respectively). Figure 6F summaries the data of frame A-E in a heat map graph showing the pattern of these gene expression from the up to down-regulated expression levels in CRC tissues compared to the adjacent normal tissues in order, where the DDIT4 shows the highest up-regulated and SULF1 the lowest down-regulated expression levels in CRC tissues compared to the adjacent normal tissues.
Additional analyses were performed to find any association between the median expression of the selected genes and the clinicopathological features of the CRC patients (Table 4). Results displayed significant relationship between TNM stage and expression of SULF1 (P = 0.023) and TPTEP1 (P < 0.01) in CRC patients. Median expression of SULF1 and TPTEP1 showed significantly increased levels in tumor tissues obtained from CRC patients with the advanced stages (frame B and C of Fig. 7). Kruskal–Wallis test showed that expression levels of DDIT4 (P = 0.048), SULF1 (P = 0.009) and TPTEP1 (P = 0.035) were significantly related to metastasis (Table 4). Our results demonstrated median expression of DDIT4 (P = 0.029) and SULF1 (P < 0.001) were significantly higher in patients with distant metastasis than patients without metastasis (frame D and E of Fig. 7) while median expression of TPTEP1 was higher in patients with only lymph node metastasis than patients without metastasis (P = 0.017) (frame F of Fig. 7). The levels of miR-181d-5p expression were found significantly related to the histologic grading of the CRC tumor (P = 0.016) (Table 4). As shown in frame G of Fig. 7, the median expression levels of miR-181d-5p was significantly higher in CRC patients with Grade 3 than patients with Grade 2 (P = 0.006). A significant relationship between median expression levels of miR-148b-3p and presence of perineural invasion was observed, where the median expression levels of miR-148b-3p in patients with the present of perineural invasion were significantly lower compared to that in patients without the perineural invasion (P = 0.021) (Table 4 and frame H of Fig. 7).
There were numerous significant correlations between the RNA expression levels of the selected genes amongst each other both in CRC tissues (Table 5) and adjacent normal tissues (Table 6) as shown by spearman correlation values. The expression levels of TPTEP1 were positively correlated with the expression levels of DDIT4 both in CRC (P = 0.019, rs: 0.37) and adjacent normal tissues (P < 0.001, rs: 0.53) (frame A of Fig. 8). Also, a significant positive correlation was found between TPTEP1 and SULF1 expression levels both in CRC tissues (P < 0.001, rs: 0.65) and to a lesser extent in adjacent normal tissues (P = 0.046, rs: 0.37) (frame B of Fig. 8). There was a negative correlation between DDIT4 and miR-148b-3p expression levels in CRC tissues (P = 0.037, rs: − 0.32) (frame C of Fig. 8). Moreover, a significant negative correlation between SULF1 and miR-148b-3p expression levels in CRC tissues (P = 0.010, rs: − 0.40) and to a greater extent, but positively, in adjacent normal tissues (P = 0.004, rs: 0.47) was found (frame D of Fig. 8). The expression levels of miR-148b-3p were positively and strongly correlated with the expression levels of miRNA-181d-5p both in CRC (P < 0.001, rs: 0.88) and adjacent normal tissues (P < 0.001, rs: 0.64) (frame E of Fig. 8).
Validation of CSC marker genes and expression levels of selected genes in CSC-enriched spheroids compared to HT-29 cancer cells
We evaluated stemness, ABC transporter, and EMT marker genes as CSC features in CSC-enriched spheroids derived from HT-29 as determined by RT-qPCR after observation of sphere morphology under the microscope (frame A and B of Fig. 9). Our results showed significantly higher expression levels of stemness genes (OCT4, SOX2, C-MYC, and KLF4), ABC transporter genes (ABCB1, ABCG2, and ABCC1) and EMT genes (TWIST1, SNAIL1, and ZEB1) in colorectal CSC-enriched spheroids compared to HT-29 cancer cells (control) (frame C-E of Fig. 9).
After detection of CSC features, the RNA expression levels of DDIT4, SULF1, TPTEP1 and miRNAs (miR-181d-5p and miR-148b-3p) were measured using RT‐qPCR in colorectal CSC-enriched spheroids and HT-29 cancer cells. RNA expression levels of DDIT4 (P = 0.042), SULF1 (P = 0.032) and TPTEP1 (P = 0.021) were significantly higher in colorectal CSC-enriched spheroids compared to the HT-29 cancer cells (frame A-C of Fig. 10). No significant difference was found in miRNAs expression levels (miR-181d-5p and miR-148b-3p) between colorectal CSC-enriched spheroids and HT-29 cancer cells.
Numerous RNAs (mRNAs, lncRNAs and miRNAs) can be used as potential biomarkers for diagnosis, prognosis and treatment in various cancers and their dysregulation shown to be associated with the development of different cancers. RNA biomarkers provide dynamics insights into cell regulation and processes compared to DNA biomarkers. They have more sensitivity and specificity than protein biomarkers . The biological roles of RNAs make them as important as the functions of proteins . The use of RNA studies in medicine has led to attracting numerous companies to develop new RNA-based diagnostic, prognostic tools, and drugs . The panels such as ThyraMIR/ThyGENX and approval of the first RNAi drug Onpattro have made the RNAs, especially ncRNAs, studies important [74, 75].
In the current study, to identify RNA biomarkers in CRC, bioinformatics analysis was applied to detect DEGs in CRC microarray data, in which 11 key genes were selected for further analysis for predictions of miRNAs and lncRNAs. Then, we evaluated the expression levels of DDIT4, SULF1, TPTEP1, miR-181d-5p and miR-148b-3p, as potential biomarkers, in CRC patients and their association with clinicopathological features.
DDIT4 is a suppressor for mammalian target of rapamycin (mTOR) signaling pathway which is induced in various cellular stress conditions such as hypoxia and DNA damage [28, 76, 77]. Moreover, DDIT4 gene has a p53 transcription-factor binding site which can play a key role in the p53-dependent tumorigenesis . Despite its repressive role on mTOR signaling pathway, up-regulation of DDIT4 has been shown to promote cell proliferation and reduce apoptotic rate in various cell types [79, 80]. DDIT4 was suitable candidate gene in order to predict miRNAs in CRC because of not only its participation in the PI3K/Akt/mTOR signaling pathway  but also its involvement in “miRNAs in cancer” pathway as shown by KEGG pathway analysis. We found RNA expression levels of DDIT4 were significantly higher in CRC tissues compared to the adjacent normal tissues. This result is in line with the study on gastric cancer displaying up-regulation of DDIT4 expression in the tumor tissues compared to the adjacent normal tissues as found by RT-qPCR and immunohistochemically staining . Since earlier findings reported that inhibition of mTOR pathway is leading to enrichment of cancer stem cells, high DDIT4 expression could be related to expression of stem-cells markers [68, 81]. As expected, up-regulation of DDIT4 expression was observed in colorectal CSC-enriched spheroids compared to HT-29 cancer cells in our study. This result supports higher expression levels of DDIT4 in the CRC patients with metastasis because CSCs play role in tumor metastasis , and are key drivers in tumor progression . These findings indicate that this higher RNA expression of DDIT4 is significantly associated with more aggressive tumor behavior. Our report is the first study to show high mRNA levels of DDIT4 expression and its clinical significance in CRC tissues as well as in colorectal CSC-enriched spheroids.
SULF1, another candidate gene, is a subtype of proteinase released by various cells in the extra cell matrix (ECM) and alters its function by modifying HS. This alteration affects several signaling molecules toward the development and spread of cancer in the microenvironment . Several experimental studies reported SULF1 as a tumor suppressor effector and its down-regulation levels related to several cancers such as pancreatic, ovarian and gastric cancer [85,86,87]. While, some other studies have shown up-regulation of SULF1 expression in gastric, colorectal and bladder cancer [31, 88, 89]. Our bioinformatics analysis showed up-regulation of SULF1 expression levels. This data is consistent with the previous results on ONCOMINE database showing increased expression levels of SULF1 in CRC tissues compared to the adjacent normal tissues . In our study, RT-qPCR data showed down-regulation of SULF1 expression levels in CRC tissues compared to the adjacent normal tissues, although its expression levels showed significantly increased levels in patients with more advanced stage and metastasis of the tumor. This result is in line with the previous findings observing down-regulating of SULF1 in early stage of ovarian tumors [86, 91]. The increased SULF1 expression levels have been also reported at the later stages of malignancy progression in CRC patients [32, 33]. It has been described that SULF1 has ambivalent functions and there is insufficient information to understand the conflicting results regarding the role of SULF1 in cancer [92, 93]. The tumor suppressor effect of SULF1 was described under hypoxic conditions in solid tumors. The level reduction of SULF1 in such environments causes increasing in 6-O-sulfate on HSPGs which subsequently leads to increasing of the fibroblast growth factor (FGF) signaling and cancer progression . Besides, the oncogenic effect of SULF1 was proposed due to the high-affinity of HS-Wnt complex. In fact, extracellular removal of the 6-O-sulfate on the HSPGs by SULF1 allows initiation of the Wnt signaling [94, 95]. Evidence suggests that overexpression of SULF1 is related to expression of EMT genes and can promote EMT in human hepatocellular carcinoma . In this regard, we measured SULF1 expression in colorectal CSC-enriched spheroids with increased EMT gene expression and observed significantly higher expression levels of SULF1 in colorectal CSC-enriched spheroids compared to the HT-29 cells. This finding is in line with the previous data showing up-regulation of SULF1 expression levels in high metastatic colorectal cancer cell lines  and breast CSCs . Increased SULF1 expression levels in patients with distant metastasis, advanced stages of the tumor, as well as colorectal CSC-enriched spheroids indicate that the levels of SULF1 expression is being increased by tumor progression in CRC. Despite dysregulation of SULF1 expression levels in CRC, such challenging observations make it difficult to offer SULF1 as a “biomarker”. Therefore, more studies are needed to reveal the dual roles of SULF1 and its expression pattern in cancer patients.
In the present study, we also investigated expression of some ncRNAs including miRNAs (miR-181d-5p and miR-148b-3p) and TPTEP1 in CRC tissues. Previous reports indicated that miR-181d contributes in regulation of Akt pathway in breast cancer and CRC cell glycolysis which acts as an oncomiR [97, 98]. We demonstrated that miR-181d-5p is significantly up-regulated in tumor tissues compared with the adjacent normal tissues in CRC patients. This data is in good agreement with the previous data in CRC patients . Moreover, the association between overexpression of miR-181d-5p and high-grade tumor cells may indicate a possible influence of increased miR-181d-5p expression in the progression of cancer. Despite what was previously described in the breast cancer cells and CRC patients about association between high expression of miR-181d-5p and increased invasion and migration of the tumors [97, 99], our data analysis didn't show any significant difference of miR-181d-5p expression levels with various metastasis groups in CRC patients. In line with this data, our results didn’t show any significant different in expression levels of miR-181d-5p between HT-29 cancer cells and colorectal CSC-enriched spheroids. While up-regulation of miR-181 family has been previously observed in the liver cancer stem/progenitor cells .
Cancer reports displayed that miR-148b, especially miR-148b-3p, plays an important role as a tumor suppressor by influencing on cell growth and proliferation , apoptosis , metastasis dissemination and cancer therapy responses . Our result demonstrated a lower expression of miR-148b-3p in CRC tissues compared with the adjacent normal tissues which is in line with previous result in CRC patients . Also, we observed down-regulation of miR-148b-3p in patients with vascular invasion compared to those without this invasion. No significant difference in miR-148b-3p expression was revealed between colorectal CSC-enriched spheroids and HT-29 cancer cells, nor between metastasis groups of CRC patients. In contrast to our findings, decreased expression and suppressor role of miR-148b-3p has been previously reported in the hepatic CSCs [69, 102].
Our investigation, for the first time, identified dysregulated expression of TPTEP1 in CRC patients. Contrary to the lung  and liver  cancer studies, the expression of TPTEP1 showed an up-regulation pattern in our CRC tissues compared to the adjacent normal tissues. Moreover, we observed higher expression of TPTEP1 in colorectal CSC-enriched spheroids than HT-29 cancer cells. As expected, based on predictions, DDIT4 and SULF1 expression levels were significantly correlated with the TPTEP1 expression levels in CRC. This result may be related to interaction between these RNAs amongst each other and can support findings about predicted binding sites based on bioinformatics algorithms for TPTEP1. According to the predicted binding site, DDIT4 in RNA level from 3′UTR region interacts with TPTEP1. It is remarkable that 3' UTRs play critical roles in gene expression regulation through bindings ncRNAs . Also, RNA expression levels of DDIT4 and SULF1 were significantly correlated negatively with miR-148b-3p expression levels in CRC tissues. These correlations may be explained by regulatory effects of miR-148b-3p expression on these RNAs as predicted based on the mRNA-miRNA network. We aware that our research has limitations to describe in details these regulatory effects and further studies are warranted to understand the relationship between these RNAs.
Overexpression of DDIT4 and TPTEP1 in CRC patients with metastasis and advanced stages as well as in colorectal CSC-enriched spheroids indicates that increased RNA expression of these markers may be useful indicators of more aggressive tumor behavior and further disease progression in CRC patients. Moreover, correlations and predicted interactions of TPTEP1 and miR-148b-3p with DDIT4 and SULF1 in mRNA level might be due to the regulatory effects of these RNAs amongst each other. According to the expression differences of DDIT4, SULF1, TPTEP1, miR-181d-5p, and miR-148b-3p in CRC tissues compared to the adjacent normal tissues, we believe our results provide a valuable resource in order to find biomarkers clinicopathologically relevant to CRC patients. From these findings, we are able to conclude that analysis of hub mRNA–miRNA genes can help to predict some important lncRNAs which are dysregulated in CRC patients.
Availability of data and materials
The analyzed data during the current study are available from the corresponding authors on reasonable request.
Cancer stem cells
Long noncoding RNA
Heparan sulfate proteoglycans
Differentially expressed genes
Extra cell matrix
Fibroblast growth factors
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(6):394–424.
Choong MK, Tsafnat G. Genetic and epigenetic biomarkers of colorectal cancer. Clin Gastroenterol Hepatol. 2012;10(1):9–15.
Ayob AZ, Ramasamy TS. Cancer stem cells as key drivers of tumour progression. J Biomed Sci. 2018;25(1):20.
Zhou Y, Xia L, Wang H, Oyang L, Su M, Liu Q, et al. Cancer stem cells in progression of colorectal cancer. Oncotarget. 2017;9(70):33403–15.
Anderson EC, Hessman C, Levin TG, Monroe MM, Wong MH. The role of colorectal cancer stem cells in metastatic disease and therapeutic response. Cancers. 2011;3(1):319–39.
Wolpin BM, Mayer RJ. Systemic treatment of colorectal cancer. Gastroenterology. 2008;134(5):1296-310.e1.
Weiser MR. AJCC 8th Edition: colorectal cancer. Ann Surg Oncol. 2018;25(6):1454–5.
Arrington AK, Heinrich EL, Lee W, Duldulao M, Patel S, Sanchez J, et al. Prognostic and predictive roles of KRAS mutation in colorectal cancer. Int J Mol Sci. 2012;13(10):12153–68.
Salazar R, Roepman P, Capella G, Moreno V, Simon I, Dreezen C, et al. Gene expression signature to improve prognosis prediction of stage II and III colorectal cancer. J Clin Oncol. 2011;29(1):17–24.
Sanz-Pamplona R, Berenguer A, Cordero D, Riccadonna S, Sole X, Crous-Bou M, et al. Clinical value of prognosis gene expression signatures in colorectal cancer: a systematic review. PLoS ONE. 2012;7(11):e48877.
Nugent M, Miller N, Kerin MJ. Circulating miR-34a levels are reduced in colorectal cancer. J Surg Oncol. 2012;106(8):947–52.
Alves Martins BA, de Bulhoes GF, Cavalcanti IN, Martins MM, de Oliveira PG, Martins AMA. Biomarkers in colorectal cancer: the role of translational proteomics research. Front Oncol. 2019;9:1284.
Anandaram H. A review on application of biomarkers in the field of bioinformatics & nanotechnology for individualized cancer treatment. MOJ Proteom Bioinform. 2017;5(6):179–84.
Supplitt S, Karpinski P, Sasiadek M, Laczmanska I. Current achievements and applications of transcriptomics in personalized cancer medicine. Int J Mol Sci. 2021;22(3):1422.
Iyer MK, Niknafs YS, Malik R, Singhal U, Sahu A, Hosono Y, et al. The landscape of long noncoding RNAs in the human transcriptome. Nat Genet. 2015;47(3):199–208.
Berindan-Neagoe I, Monroig Pdel C, Pasculli B, Calin GA. MicroRNAome genome: a treasure for cancer diagnosis and therapy. CA Cancer J Clin. 2014;64(5):311–36.
Yamamura S, Imai-Sumida M, Tanaka Y, Dahiya R. Interaction and cross-talk between non-coding RNAs. Cell Mol Life Sci. 2018;75(3):467–84.
Okugawa Y, Toiyama Y, Goel A. An update on microRNAs as colorectal cancer biomarkers: where are we and what’s next? Expert Rev Mol Diagn. 2014;14(8):999–1021.
Huang T, Alvarez A, Hu B, Cheng SY. Noncoding RNAs in cancer and cancer stem cells. Chin J Cancer. 2013;32(11):582–93.
Macfarlane LA, Murphy PR. MicroRNA: biogenesis, function and role in cancer. Curr Genomics. 2010;11(7):537–61.
Fukunaga T, Hamada M. RIblast: an ultrafast RNA-RNA interaction prediction system based on a seed-and-extension approach. Bioinformatics. 2017;33(17):2666–74.
Dweep H, Sticht C, Gretz N. In-silico algorithms for the screening of possible microRNA binding sites and their interactions. Curr Genomics. 2013;14(2):127–36.
Ren X. Genome-wide analysis reveals the emerging roles of long non-coding RNAs in cancer. Oncol Lett. 2020;19(1):588–94.
Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011;12(12):861–74.
Lopez-Urrutia E, Bustamante Montes LP, Ladron de Guevara Cervantes D, Perez-Plasencia C, Campos-Parra AD. Crosstalk between long non-coding RNAs, micro-RNAs and mRNAs: deciphering molecular mechanisms of master regulators in cancer. Front Oncol. 2019;9(669):669.
Terai G, Iwakiri J, Kameda T, Hamada M, Asai K. Comprehensive prediction of lncRNA-RNA interactions in human transcriptome. BMC Genomics. 2016;17(Suppl 1):12.
Brugarolas J, Lei K, Hurley RL, Manning BD, Reiling JH, Hafen E, et al. Regulation of mTOR function in response to hypoxia by REDD1 and the TSC1/TSC2 tumor suppressor complex. Genes Dev. 2004;18(23):2893–904.
Tirado-Hurtado I, Fajardo W, Pinto JA. DNA damage inducible transcript 4 gene: the switch of the metabolism as potential target in cancer. Front Oncol. 2018;8(106):106.
Pinto JA, Rolfo C, Raez LE, Prado A, Araujo JM, Bravo L, et al. In silico evaluation of DNA damage inducible Transcript 4 gene (DDIT4) as prognostic biomarker in several malignancies. Sci Rep. 2017;7(1):1526.
Sasisekharan R, Venkataraman G. Heparin and heparan sulfate: biosynthesis, structure and function. Curr Opin Chem Biol. 2000;4(6):626–31.
Vicente CM, Lima MA, Yates EA, Nader HB, Toma L. Enhanced tumorigenic potential of colorectal cancer cells by extracellular sulfatases. Mol Cancer Res. 2015;13(3):510–23.
Suhovskih AV, Aidagulova SV, Kashuba VI, Grigorieva EV. Proteoglycans as potential microenvironmental biomarkers for colon cancer. Cell Tissue Res. 2015;361(3):833–44.
Babel I, Barderas R, Diaz-Uriarte R, Moreno V, Suarez A, Fernandez-Acenero MJ, et al. Identification of MST1/STK4 and SULF1 proteins as autoantibody targets for the diagnosis of colorectal cancer by using phage microarrays. Mol Cell Proteomics. 2011;10(3):M110.
Pop-Bica C, Pintea S, Cojocneanu-Petric R, Del Sal G, Piazza S, Wu ZH, et al. MiR-181 family-specific behavior in different cancers: a meta-analysis view. Cancer Metastasis Rev. 2018;37(1):17–32.
Dai W, He J, Zheng L, Bi M, Hu F, Chen M, et al. miR-148b-3p, miR-190b, and miR-429 regulate cell progression and act as potential biomarkers for breast cancer. J Breast Cancer. 2019;22(2):219–36.
Celano M, Rosignolo F, Maggisano V, Pecce V, Iannone M, Russo D, et al. MicroRNAs as biomarkers in thyroid carcinoma. Int J Genomics. 2017;2017:6496570.
Arambula-Meraz E, Bergez-Hernandez F, Leal-Leon E, Romo-Martinez E, Picos-Cardenas V, Luque-Ortega F, et al. Expression of miR-148b-3p is correlated with overexpression of biomarkers in prostate cancer. Genet Mol Biol. 2020;43(1):e20180330.
Song Y, Xu Y, Wang Z, Chen Y, Yue Z, Gao P, et al. MicroRNA-148b suppresses cell growth by targeting cholecystokinin-2 receptor in colorectal cancer. Int J Cancer. 2012;131(5):1042–51.
Song YX, Yue ZY, Wang ZN, Xu YY, Luo Y, Xu HM, et al. MicroRNA-148b is frequently down-regulated in gastric cancer and acts as a tumor suppressor by inhibiting cell proliferation. Mol Cancer. 2011;10:1.
Cao F, Wang Z, Feng Y, Zhu H, Yang M, Zhang S, et al. lncRNA TPTEP1 competitively sponges miR3285p to inhibit the proliferation of nonsmall cell lung cancer cells. Oncol Rep. 2020;43(5):1606–18.
Ding H, Liu J, Zou R, Cheng P, Su Y. Long non-coding RNA TPTEP1 inhibits hepatocellular carcinoma progression by suppressing STAT3 phosphorylation. J Exp Clin Cancer Res. 2019;38(1):189.
Fattahi F, Kiani J, Khosravi M, Vafaei S, Mohammadi A, Madjd Z, et al. Enrichment of up-regulated and down-regulated gene clusters using gene ontology, miRNAs and lncRNAs in colorectal cancer. Comb Chem High Throughput Screen. 2019;22(8):534–45.
Pinero J, Bravo A, Queralt-Rosinach N, Gutierrez-Sacristan A, Deu-Pons J, Centeno E, et al. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017;45(D1):D833–9.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.
Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: network analysis and visualization of proteomics data. J Proteome Res. 2019;18(2):623–32.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Sidiropoulos K, Viteri G, Sevilla C, Jupe S, Webber M, Orlic-Milacic M, et al. Reactome enhanced pathway visualization. Bioinformatics. 2017;33(21):3461–7.
Huang R, Grishagin I, Wang Y, Zhao T, Greene J, Obenauer JC, et al. The NCATS BioPlanet - an integrated platform for exploring the universe of cellular signaling pathways for toxicology, systems biology, and chemical genomics. Front Pharmacol. 2019;10:445.
Slenter DN, Kutmon M, Hanspers K, Riutta A, Windsor J, Nunes N, et al. WikiPathways: a multifaceted pathway database bridging metabolomics to other omics research. Nucleic Acids Res. 2018;46(D1):D661–7.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. ELife. 2015;4:e05005.
Paraskevopoulou MD, Georgakilas G, Kostoulas N, Vlachos IS, Vergoulis T, Reczko M, et al. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013;41:W169-73.
Jiang Q, Wang Y, Hao Y, Juan L, Teng M, Zhang X, et al. miR2Disease: a manually curated database for microRNA deregulation in human disease. Nucleic Acids Res. 2009;37:D98-104.
Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.
Fan Y, Siklenka K, Arora SK, Ribeiro P, Kimmins S, Xia J. miRNet - dissecting miRNA-target interactions and functional associations through network-based visual analysis. Nucleic Acids Res. 2016;44(W1):W135–41.
Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNAorg resource: targets and expression. Nucleic Acids Res. 2008;36:D149-53.
Miao YR, Liu W, Zhang Q, Guo AY. lncRNASNP2: an updated database of functional SNPs and mutations in human and mouse lncRNAs. Nucleic Acids Res. 2018;46(D1):D276–80.
Das S, Ghosal S, Sen R, Chakrabarti J. lnCeDB: database of human long noncoding RNA acting as competing endogenous RNA. PLoS ONE. 2014;9(6):e98965.
Iwakiri J, Terai G, Hamada M. Computational prediction of lncRNA-mRNA interactions by integrating tissue specificity in human transcriptome. Biol Direct. 2017;12(1):15.
Kato Y, Sato K, Hamada M, Watanabe Y, Asai K, Akutsu T. RactIP: fast and accurate prediction of RNA-RNA interaction using integer programming. Bioinformatics. 2010;26(18):i460–6.
Mann M, Wright PR, Backofen R. IntaRNA 2.0: enhanced and customizable prediction of RNA-RNA interactions. Nucleic Acids Res. 2017;45(W1):W435–9.
Naseri M, Zoller M, Hadjati J, Ghods R, Ranaei Pirmardan E, Kiani J, et al. Dendritic cells loaded with exosomes derived from cancer stem cell-enriched spheroids as a potential immunotherapeutic option. J Cell Mol Med. 2021;25(7):3312–26.
Eini L, Naseri M, Karimi-Busheri F, Bozorgmehr M, Ghods R, Madjd Z. Primary colonospheres maintain stem cell-like key features after cryopreservation. J Cell Physiol. 2020;235(3):2452–63.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. 2012;13:134.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Pinto JA, Bravo L, Chirinos LA, Vigil CE. Expression of DDIT4 is correlated with NOTCH1 and high molecular risk in acute myeloid leukemias. Blood. 2016;128(22):5254.
Liu Q, Xu Y, Wei S, Gao W, Chen L, Zhou T, et al. miRNA-148b suppresses hepatic cancer stem cell by targeting neuropilin-1. Biosci Rep. 2015;35(4):e00229.
Ji J, Yamashita T, Budhu A, Forgues M, Jia H-L, Li C, et al. Identification of microRNA-181 by genome-wide screening as a critical player in EpCAM-positive hepatic cancer stem cells. Hepatology. 2009;50(2):472–80.
Stebbins NW. Decoding structure-function relationships of glycans [dissertation on the Internet]. Massachusetts Institute of Technology; 2017. http://hdl.handle.net/1721.1/110887.
Xi X, Li T, Huang Y, Sun J, Zhu Y, Yang Y, et al. RNA biomarkers: frontier of precision medicine for cancer. Noncoding RNA. 2017;3(1):9.
Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, et al. Landscape of transcription in human cells. Nature. 2012;489(7414):101–8.
Bonneau E, Neveu B, Kostantin E, Tsongalis GJ, De Guire V. How close are miRNAs from clinical practice? A perspective on the diagnostic and therapeutic market. EJIFCC. 2019;30(2):114–27.
Slack FJ, Chinnaiyan AM. The role of non-coding RNAs in oncology. Cell. 2019;179(5):1033–55.
Shoshani T, Faerman A, Mett I, Zelin E, Tenne T, Gorodin S, et al. Identification of a novel hypoxia-inducible factor 1-responsive gene, RTP801, involved in apoptosis. Mol Cell Biol. 2002;22(7):2283–93.
Gordon BS, Steiner JL, Williamson DL, Lang CH, Kimball SR. Emerging role for regulated in development and DNA damage 1 (REDD1) in the regulation of skeletal muscle metabolism. Am J Physiol Endocrinol Metab. 2016;311(1):E157–74.
Wei CL, Wu Q, Vega VB, Chiu KP, Ng P, Zhang T, et al. A global map of p53 transcription-factor binding sites in the human genome. Cell. 2006;124(1):207–19.
Du F, Sun L, Chu Y, Li T, Lei C, Wang X, et al. DDIT4 promotes gastric cancer proliferation and tumorigenesis through the p53 and MAPK pathways. Cancer Commun. 2018;38(1):45.
Chang B, Liu G, Yang G, Mercado-Uribe I, Huang M, Liu J. REDD1 is required for RAS-mediated transformation of human ovarian epithelial cells. Cell Cycle. 2009;8(5):780–6.
Bhola NE, Jansen VM, Koch JP, Li H, Formisano L, Williams JA, et al. Treatment of triple-negative breast cancer with TORC1/2 inhibitors sustains a drug-resistant and notch-dependent cancer stem cell population. Cancer Res. 2016;76(2):440–52.
Shiozawa Y, Nie B, Pienta KJ, Morgan TM, Taichman RS. Cancer stem cells and their role in metastasis. Pharmacol Ther. 2013;138(2):285–93.
Liao W-T, Ye Y-P, Deng Y-J, Bian X-W, Ding Y-Q. Metastatic cancer stem cells: from the concept to therapeutics. Am J Stem Cells. 2014;3(2):46–62.
Hammond E, Khurana A, Shridhar V, Dredge K. The role of heparanase and sulfatases in the modification of heparan sulfate proteoglycans within the tumor microenvironment and opportunities for novel cancer therapeutics. Front Oncol. 2014;4:195.
Lai JP, Sandhu DS, Shire AM, Roberts LR. The tumor suppressor function of human sulfatase 1 (SULF1) in carcinogenesis. J Gastrointest Cancer. 2008;39(1–4):149–58.
Lai J, Chien J, Staub J, Avula R, Greene EL, Matthews TA, et al. Loss of HSulf-1 up-regulates heparin-binding growth factor signaling in cancer. J Biol Chem. 2003;278(25):23107–17.
Narita K, Staub J, Chien J, Meyer K, Bauer M, Friedl A, et al. HSulf-1 inhibits angiogenesis and tumorigenesis in vivo. Cancer Res. 2006;66(12):6025–32.
Junnila S, Kokkola A, Mizuguchi T, Hirata K, Karjalainen-Lindsberg ML, Puolakkainen P, et al. Gene expression analysis identifies over-expression of CXCL1, SPARC, SPP1, and SULF1 in gastric cancer. Genes Chromosomes Cancer. 2010;49(1):28–39.
Lee HY, Yeh BW, Chan TC, Yang KF, Li WM, Huang CN, et al. Sulfatase-1 overexpression indicates poor prognosis in urothelial carcinoma of the urinary bladder and upper tract. Oncotarget. 2017;8(29):47216–29.
Rosen SD, Lemjabbar-Alaoui H. Sulf-2: an extracellular modulator of cell signaling and a cancer target candidate. Expert Opin Ther Targets. 2010;14(9):935–49.
Khurana A, Beleford D, He X, Chien J, Shridhar V. Role of heparan sulfatases in ovarian and breast cancer. Am J Cancer Res. 2013;3(1):34–45.
Vives RR, Seffouh A, Lortat-Jacob H. Post-synthetic regulation of HS structure: The Yin and Yang of the sulfs in cancer. Front Oncol. 2014;3:331.
Pascale RM, Calvisi DF, Feo F. Sulfatase 1: a new Jekyll and Hyde in hepatocellular carcinoma? Transl Gastroenterol Hepatol. 2016;1:43.
Ai X, Do AT, Lozynska O, Kusche-Gullberg M, Lindahl U, Emerson CP Jr. QSulf1 remodels the 6-O sulfation states of cell surface heparan sulfate proteoglycans to promote Wnt signaling. J Cell Biol. 2003;162(2):341–51.
Nawroth R, van Zante A, Cervantes S, McManus M, Hebrok M, Rosen SD. Extracellular sulfatases, elements of the Wnt signaling pathway, positively regulate growth and tumorigenicity of human pancreatic cancer cells. PLoS ONE. 2007;2(4):e392.
Dhanasekaran R, Nakamura I, Hu C, Chen G, Oseini AM, Seven ES, et al. Activation of the transforming growth factor-beta/SMAD transcriptional pathway underlies a novel tumor-promoting role of sulfatase 1 in hepatocellular carcinoma. Hepatology. 2015;61(4):1269–83.
Guo X, Zhu Y, Hong X, Zhang M, Qiu X, Wang Z, et al. miR-181d and c-myc-mediated inhibition of CRY2 and FBXL3 reprograms metabolism in colorectal cancer. Cell Death Dis. 2017;8(7):e2958.
Strotbek M, Schmid S, Sanchez-Gonzalez I, Boerries M, Busch H, Olayioye MA. miR-181 elevates Akt signaling by co-targeting PHLPP2 and INPP4B phosphatases in luminal breast cancer. Int J Cancer. 2017;140(10):2310–20.
Wang H, Wei H, Wang J, Li L, Chen A, Li Z. MicroRNA-181d-5p-containing exosomes derived from CAFs promote EMT by regulating CDX2/HOXA5 in breast cancer. Mol Ther Nucleic Acids. 2020;19:654–67.
Wu Y, Liu GL, Liu SH, Wang CX, Xu YL, Ying Y, et al. MicroRNA-148b enhances the radiosensitivity of non-Hodgkin’s Lymphoma cells by promoting radiation-induced apoptosis. J Radiat Res. 2012;53(4):516–25.
Cimino D, De Pitta C, Orso F, Zampini M, Casara S, Penna E, et al. miR148b is a major coordinator of breast cancer progression in a relapse-associated microRNA signature by targeting ITGA5, ROCK1, PIK3CA, NRAS, and CSF1. FASEB J. 2013;27(3):1223–35.
Xu Y, Xie Y, Wang X, Chen X, Liu Q, Ying M, et al. Identification of cancer stem cells from hepatocellular carcinoma cell lines and their related microRNAs. Oncol Rep. 2013;30(5):2056–62.
Szczesniak MW, Makalowska I. lncRNA-RNA Interactions across the Human Transcriptome. PLoS ONE. 2016;11(3):e0150353.
We are grateful to Prof. Carlos Cordon-Cardo (University of Icahn School of Medicine) and Dr. Fatemeh Fattahi (University of Michigan Medical School) for their feedback and pre-reviewing our article.
This research was supported by a grant of the Iran University of Medical Sciences (Grant: 95-04-87-30106).
Ethics approval and consent to participate
The study protocol was approved by the Iran University of Medical Sciences Research (IUMS) Ethics Committee (IR.IUMS.REC 1395.95-04-87-30106). Informed consent was obtained from all individual participants at the time of sample collection with routine consent forms.
Consent for publication
Patients signed informed consent regarding publishing their data.
All authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Information obtained in the process of identification and selection of key genes including up-regulated genes with combined gene score > 3 of previous analysis (PMID: 31654507), up-regulated genes in carcinogenesis & colorectal cancer of DisGeNET with pvalue < 0.0001, the largest cluster of k-means clustering in STRING (confidence ≥ 0.4), and key genes.
Additional file 2: Figure S1.
Protein–protein interaction (PPI) network analysis. PPI network explored the interactions between the 370 up-regulated genes with confidence ≥ 0.4. Five main clusters were obtained of the k-means algorithm that five colors were applied for indication gene clusters (each color indicates a cluster gene). Figure S2. Pathway analysis for the largest cluster covering 167 genes on Enrichr. Top ten results of pathway analysis that was performed based on BioPlanet, KEGG, WikiPathways, and Reactome libraries. Figure S3. Gene ontology analysis for the largest cluster covering 167 genes on Enrichr. Top ten results of gene ontology (GO) based on p-value. Results of GO analysis was contained cellular component (CC), biological process (BP), molecular function (MF) and Jensen diseases.
Additional file 3: Table S2.
Information of miRNA-mRNA network for 11 key genes.
Additional file 4: Table S3.
All of the predicted lncRNAs for 9 hub miRNAs.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Fattahi, F., Kiani, J., Alemrajabi, M. et al. Overexpression of DDIT4 and TPTEP1 are associated with metastasis and advanced stages in colorectal cancer patients: a study utilizing bioinformatics prediction and experimental validation. Cancer Cell Int 21, 303 (2021). https://doi.org/10.1186/s12935-021-02002-x
- Colorectal cancer (CRC)
- Bioinformatics analysis