Back to Journals » OncoTargets and Therapy » Volume 14

Identification of microRNA Signature and Key Genes Between Adenoma and Adenocarcinomas Using Bioinformatics Analysis

Authors Shi X, Gao GY , Shen J

Received 16 May 2021

Accepted for publication 8 August 2021

Published 4 September 2021 Volume 2021:14 Pages 4707—4720

DOI https://doi.org/10.2147/OTT.S320469

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 3

Editor who approved publication: Dr Jianmin Xu



Xinya Shi,1,* Guang Yu Gao,2,* Jiaofeng Shen2

1Department of Oncology, Changshu Second People’s Hospital, Suzhou, 215004, People’s Republic of China; 2Department of Oncology, The Second Affiliated Hospital of Soochow University, Suzhou, 215004, People’s Republic of China

*These authors contributed equally to this work

Correspondence: Xinya Shi
Department of Oncology, Changshu Second People’s Hospital, 68 Haiyu South Road, Suzhou, 215004, Jiangsu, People’s Republic of China
Email [email protected]

Background: In worldwide, colorectal cancer (CRC) is very common and the mechanisms remain unclear. This study aims to identify between adenomas with epithelial dislocation (false invasion) and adenomas with early adenocarcinoma (true invasion).
Methods: GSE41655 and GSE57965 datasets were obtained in the Gene Expression Omnibus (GEO) database. microRNA expression profiles and clinicopathological data from the TCGA (The Cancer Genome Atlas) database were downloaded to further validate the results in GEO. GEO software and the GEO2R calculation method were used to analyze two gene profiles. The co-expression of differentially expressed microRNAs (DEMs) and genes (DEGs) were identified and searched in the FunRich databases for pathway and ontology analysis. Cytoscape was utilized to construct the mRNA-microRNA network. Validation of gene expression levels was conducted by online databases and qRT-PCR and IHC experiments.
Results: In total, 6 DEMs and 34 DEGs are selected after calculating. KEGG results indicated that genes are enriched in certain tumor associated pathways. Four out of 6 microRNAs had a significant relationship with the overall survival (P < 0.05) and showed a good performance in predicting the survival risk of patients with colorectal carcinoma. Furthermore, expression levels of hsa-miR-455 and hsa-miR-125a were then verified by qRT-PCR which all target BCL2L12. IHC results showed that the expression level of BCL2L12 was higher in adenocarcinoma than in adenoma. Based on the selected gene, the top 10 small molecules were screened out as potential drugs.
Conclusion: By using microarray and bioinformatics analyses, DEMs and DEGs were selected and a complete gene network was constructed. To our knowledge, BCL2L12 and related molecules including hsa-miR-455 and hsa-miR-125a were firstly identified as potential biomarkers in the progression from adenoma to adenocarcinoma.

Keywords: microRNA, colorectal cancer, adenoma, prognostic signature

Background

Colorectal carcinoma (CRC) is the second leading cause of malignant cancer and cancer-related death worldwide. Unfortunately, 25% to 40% of people will develop recurrence despite the possibility of surgical cure. Compared with patients with earlier stage, the cure rate and overall survival of patients with advanced colorectal carcinoma are still very poor.1 Surgical treatment is still the cornerstone of the treatment of locally advanced colorectal carcinoma. At present, effective treatment for unresectable metastatic cancers is lacked, as well as tumors with the worse effect of chemotherapy and radiotherapy.2 CRC is a multi-step process from normal epithelial cells to adenomas and adenocarcinoma, and finally to different organs. In 1990, the development mode of CRC was introduced, among which APC, KRAS, TP53, and DCC were considered as the genes to promote CRC development. Since then, many studies have studied the molecular mechanism of CRC. It suggested that colorectal cancer is caused by the accumulation of genetic and genetic factors (such as Wnt, PIK3CA, and TGF - β) that changed the signaling pathway. In the pathogenesis of CRC, three main recognized pathways are chromosome instability, microsatellite instability, and CpG island methylation phenotype. Many CRC lacks the changes described above, which indicates that there are other mechanisms in the development of CRC.3 Moreover, there are significant differences in the patterns of gene methylation between the right and left colon. Notably, the prevalence of promoter methylation of the mismatch repair gene hMLH1 and the O-6-methylguanine-DNA methyltransferase (MGMT) is significantly greater in normal right colon mucosa, especially in older women, suggesting epigenetic aberrations in preneoplastic right colon mucosa that may be reflected in subsequent right-sided adenocarcinoma biology. Besides, there is differential prognosis by stage between patients with right- and left-sided CRC. Retrospective studies suggest that right-sided tumors have a slightly better prognosis in stage II CRC, but a slightly worse prognosis in stage III disease, which is probably associated with the higher prevalence of good-prognosis microsatellite unstable (MSI-high) tumors in right-sided stage II cancers.4

MicroRNAs are a class of non-coding single-stranded RNA molecules, which are about 22 nucleotides in length and encoded by endogenous genes.5 In the whole development process of organisms, microRNA may have biological functions such as regulating early cell development, participating in cell differentiation and tissue development, and regulating gene expression, the main role of which is to regulate gene expression.6 In particular, Wang et al found that microRNA-217 inhibited high glucose induced proliferation and migration of vascular smooth muscle cells by targeting ROCK1.7 Ma et al found that microRNA-214 was reduced in cutaneous squamous cell carcinoma tissues and cells, and patients with microRNA-214 overexpression had a higher survival rate.8 Yan et al also reported that microRNA-224, microRNA-147b, and microRNA-31 are related to lymph node metastasis and overall survival for lung cancer by regulating PRPF4B, WDR82, or NR3C2.9 High expression of microRNA-21 has been correlated with shorter disease-free-survival in stage IIcolorectal cancer patients. Apart from microRNA-21, the classifiers based on microRNAs including microRNA-20a-5p, microRNA-103a-3p, microRNA-106a-5p, and microRNA-143-5p have been reported as novel predictive markers for the recurrence of stage II disease. Furthermore, microRNA-320e has also been recognized as a prognostic biomarker; high level of microRNA-320e is correlated with more advanced disease, recurrence, and dismal prognosis in those with stage II and III colorectal cancer. Moreover, there is some evidence for contribution of microRNA-148a to the carcinogenesis of CRC. It has also been reported that postoperative plasma microRNA-31, microRNA-141, and microRNA-16 are suggested biomarkers of disease recurrence after the surgical resection. Finally, microRNA-429 expression is upregulated in colorectal cancer tissue and as such closely related to tumor size, lymph node involvement, and distant metastases, whereas it leads to shorter survival.10 However, the regulations of microRNAs in the development of adenoma to carcinoma remain unknown.

In this article, microarray data from the GEO database and colorectal cancer sample data in the TCGA database were used for identifying differently expressed microRNAs (DEMs) between adenoma tissues and adenocarcinoma tissues. By using a variety of mRNA and microRNA-related functional databases and performing verification experiments, several genes associated with adenoma progression to cancer and subsequent pathways were identified. For the study, we selected 2 different sets of data, including adenoma and cancer, 1 of which was not published.

Methods

Microarray Data

GEO database, the full name of gene expression omnibus, is a gene expression database created and maintained by National Biotechnology Information Center NCBI. It was constructed in 2000 and contains high-throughput gene expression data submitted by research institutions around the world. In our study, GSE41655 and GSE57965 datasets were downloaded from GEO.

GSE41655, including 20 healthy colonics, 20 colorectal adenomas, and 20 colorectal cancer tissue. MicroRNA expression profiling analysis of these tissues was constructed on Affymetrix Human Transcriptome Array 2.0 (GPL17586). Dataset GSE57965, including 5 adenomas and 9 adenocarcinomas, was processed by Affymetrix Human Gene Expression Array (PrimeView) (GPL15207).

DEMs Analysis

GEO2R, an R-associated web application, was applied to filtrate DEMs between adenoma samples and adenocarcinoma samples. We also used R software to analyze two sets of data. The p < 0.05 and| log FC| ≥ 1 were considered as cutoff criterion.11

Functional and Pathway Enrichment Analysis

GO functional analysis and KEGG pathway analysis were performed to predict the potential functions of the DEMs by using the Database for Annotation, Visualization, and Integrated Discovery (FunRich; http://www.funrich.org). Upregulated and downregulated DEMs were submitted to the FunRich online program. The top 10 items of the cellular component (CC), biological process (BP), and molecular function (MF) categories were then sorted and presented in the form of pie graphs. Cytoscape software was utilized to perform KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis. ClueGO is a Cytoscape plug-in, which visualizes non-redundant biological terms for a large number of gene clusters in functional packet networks. By uploading target genes, the network graph according to kappa statistics was established.

Prediction of Potential DEMs Target mRNAs and microRNA-mRNA Regulatory Network

MicroRNAs inhibit the expression of target genes mainly by binding to target mRNA, or promoting the degradation of mRNA, or hindering its translation. Accurate and rapid prediction of microRNA target genes by bioinformatics methods can provide clues for the study of microRNA function. Upregulated and downregulated DEMs were submitted to the FunRich online program to achieve target genes. Besides, GSE57965 was obtained from the GEO database. By combining the target genes from FunRich software and differential expression analysis of GSE57965, the intersection genes between the two results were selected and the microRNA-mRNA regulatory network was constructed by utilizing Cytoscape software.

Analysis of the microRNAs and Their Association with Prognosis of Colorectal Cancer Patients

The Kaplan–Meier plotter can evaluate the effect of genes (mRNA, microRNA, protein) on the survival rate of 21 cancer types (including breast cancer, ovarian cancer, lung cancer, and gastric cancer). The sources of the database include GEO, EGA, and TCGA.12 Each identified microRNA will be entered into this online tool to assess the survival rate of colorectal cancer patients according to the Kaplan Meier curve.

Independent Prognostic Ability of the microRNA Signature

To explore the effect of DEMs on the prognosis of patients with colorectal cancer, univariate and multivariate Cox proportional hazards regression analysis was performed on the DEMs. The microRNAs associated with the prognosis of colorectal cancer were identified, and a risk model based on TCGA information was built. Five hundred and sixty-nine patients with colorectal cancer obtained from the TCGA database. R software was utilized to divide the samples with complete survival information and DEMs expression profiles into two groups (train group and test group). To reduce the number of microRNAs with similar expression, stepwise Cox regression was performed for microRNAs with P < 0.005 to establish a prognostic model. In multivariate Cox regression analysis, we use the function of “Coxph” and “direction = both” in R language survival package.13,14 According to multivariate Cox regression analysis, these microRNAs were further researched to select important targets and constructed a risk linear model.

To further learn about the association between the identified microRNAs and the prognosis of colorectal cancer patients, we constructed a risk prediction model. The area under the curve (AUC) and receiver operating characteristic (ROC) of 3 and 5 years were calculated by the “survivalROC” software package to evaluate the predictive ability of microRNAs.

Quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR)

Forty participants, including 20 patients with adenomas and 20 patients with adenocarcinomas, were recruited from The Second Affiliated Hospital of Soochow University (Table 1). Before the study began, all the patients agreed. Tissue samples were collected from each participant before treatment began. After that, Total RNAs of tissues were isolated by TRIzol (Invitrogen, Carlsbad, CA, USA). Reverse transcription was performed to synthesize cDNA using Primescript RT Reagent (TaKaRa, Otsu, Shiga, Japan). qRT-PCR reactions were performed using SYBR reagent (TaKaRa, Otsu, Shiga, Japan) by Real-time PCR System (Applied Biosystems, Foster City, CA, USA). The following primers were used for qRT-PCR reactions: hsa-miR-455: 5ʹ-GCAGUCCAUGGGCAUAUACAC-3ʹ and hsa-miR-125a: 5ʹ-CCGTCCCTGAGACCCTTTAAC-3ʹ.

Table 1 Clinicopathological Characteristics of the Whole Series

Immunohistochemical Staining

Genes expression in adenocarcinoma tissues and normal tissues was extracted from the human protein atlas (www.proteinatlas.org). The Human Protein Atlas is a Swedish-based program initiated in 2003 to map all the human proteins in cells, tissues, and organs using an integration of various omics technologies, including antibody-based imaging, mass spectrometry-based proteomics, transcriptomics, and systems biology.

Identification of Candidate Small Molecule Drugs

The Connectivity Map (CMap) (http://www.complement.us/labweb/cmap/) is a database used for predicting potential drugs that may induce or reverse gene expression. The identified genes were input into the CMap database to Search for drugs against the transformation of adenoma to adenocarcinoma, and the enrichment scores were also calculated.

Results

Identification of the microRNAs Between Adenoma Tissues and Adenocarcinoma Tissues

R software was utilized to study the microRNA and mRNA expression profiles from the GSE41655 and GSE57965. Through the cut-off criteria (P < 0.05 and |log2FC|≥1), 27 DEMs and 692 DEGs were identified. A volcano plot and a heatmap were performed to show upregulated (red) and downregulated (blue) genes between adenomas and adenocarcinomas, respectively (Figure 1).

Figure 1 Heat map and volcano plot of differentially expressed microRNAs.

Abbreviations: N, adenoma; T, adenocarcinoma.

Screening of Potential Transcription Factors and Enrichment Analysis

For the reason that transcription factors are crucial in microRNA, FunRich is used to research the top 10 enriched transcription factors, namely EGR1, SP1, SOX1, POU2F1, RREB1, TCF3, RORA, FOXA1, SP4, and TEAD1. To further research the function of these microRNAs, we submitted them to FunRich to perform gene ontology analysis. The result showed that DEMs were most enriched in the nucleus, transcription factor activity, cell communication, and signal transduction (Figure 2). KEGG pathway analysis showed that these potential target genes were mainly enriched in 6 pathways including glycosphingolipid biosynthesis, retinol metabolism, AGE-RAGE signaling pathway in diabetic complications, pancreatic secretion, inflammatory mediator regulation of TRP channels, and renin secretion (Figure 3) (Table 2).

Table 2 KEGG Enrichment Analysis of Potential Target mRNAs

Figure 2 Screening of potential TFs and target mRNAs of differentially expressed microRNAs. (A) Identification of the potential TFs of differentially expressed microRNAs by FunRich software. (B) The top 10 of biological process, cellular component, and molecular function of the target genes of microRNAs.

Figure 3 KEGG pathway enriched by potential target mRNAs.

microRNA-mRNA Regulatory Network

By utilizing FunRich software, 1156 potential target genes were obtained and only 34 of them differentially expressed in GSE57965. The Venn Diagram showed the results we screened (Figure 4A). To show the composition and relationship of microRNA and mRNA more clearly, we built a microRNA-mRNA regulatory network by utilizing Cytoscape software. Finally, 41 essential microRNA-mRNA pairs were selected which may play a key role in the transformation of adenoma to adenocarcinoma (Figure 4B).

Figure 4 Identified target mRNAs and microRNA-mRNA regulatory network. (A) Venn Diagram of GSE41655 and GSE57965. (B) microRNA-mRNA regulatory network.

Analysis of the DEMs and Their Association with Prognosis of CRC Patients

Kaplan–Meier Plotter was used to analyze the prognosis of patients with colorectal carcinoma. After uploading six microRNAs, we got the corresponding survival curve. The results showed that overexpression of microRNA-199a, microRNA-199b, microRNA-127, and microRNA-125a (Figure 5) were associated with poor overall survival in patients with colorectal carcinoma. However, the expression level of microRNA-455 and microRNA-146b may have no obvious association with the overall survival. This indicated that the identified microRNAs may be potential targets.

Figure 5 The association between microRNAs and colorectal cancer prognosis.

Construction of Prognostic Risk Model and Predictability Evaluation

To determine the best prognostic microRNAs, 20 microRNAs were analyzed by LASSO Cox regression via glmnet package in R software. Nine microRNAs were identified to construct risk markers according to the minimum standard (Figure 6A). Then, multivariate Cox proportional risk regression analysis was performed to evaluate the independent prognostic value of 9 candidates prognostic microRNAs. Through Cox model, 7 candidate microRNAs (microRNA-125a-5p, microRNA-377, microRNA-376c, microRNA-455-3p, microRNA-126, microRNA-99a, and microRNA-193b) were identified as independent prognostic factors. Then, we combined the 7 microRNAs to construct a model to predict patient outcomes. The AUC of 3-year survival and 5-year survival of 7-microRNA was 0.809 and 0.981 respectively, which indicated that the model had a good effect in predicting the survival risk of colorectal carcinoma patients (Figure 6B). Through the risk model, CRC patients were divided into a high-risk group and a low-risk group. Results indicated that the model can predict the clinical outcomes of CRC patients. The survival status, risk score, and distribution of 7 microRNAs expressions in every CRC patient were also studied. The survival results indicated that patients with a low-risk score had a greater survival time than patients with a high-risk score (P < 0.0001) (Figure 6C and D).

Figure 6 Identification of the signature significantly associated with the survival of patients with CRC in the training group. (A) LASSO Cox regression algorithm was used to reduce the scope. (B) Time-dependent ROC curves analysis. (C and D) Risk score distribution and survival status for patients in high- and low-risk groups by the signature. LASSO, least absolute shrinkage and selection operator; ROC, receiver operating characteristic.

Verification of Potential Biomarker Expression by qRT-PCR and IHC

According to the combination of the GEO and TCGA analysis results, 2 microRNAs were identified for further research (Figure 7A). Then, the selected biomarkers including microRNA-125a-5p and microRNA-455-3p were validated in colorectal cancer tissue samples using qRT-PCR analysis. Consistent with the prediction, the results showed that the expression levels of microRNA-125a (P-value = 0.015) and microRNA-455 (P-value = 0.013) in the tissue of adenoma patients were lower than that of adenocarcinoma patients (Figure 7B). By using FunRich software, we also found that microRNA-125a-5p and microRNA-455-3p had high reliability, which all target BCL2L12. The results of immunohistochemistry also revealed that BCL2L12 was downregulated in adenoma tissues compared with adenocarcinoma tissues (Figure 7C). This mRNA expression was detected by the polyclonal antibody using immuno-histochemistry and sections were counterstained lightly with hematoxylin. Statistics analysis have shown that higher BCL2L12 protein expression in colorectal carcinomas than overall colorectal adenoma (Table 3).

Table 3 Expression of BCL2L12 Protein in Colorectal Adenoma and Carcinoma

Figure 7 Genes expression in human colorectal carcinoma specimens and adenomas. (A) Venn Diagram of GSE41655 and TCGA. (B) QRT-PCR results indicate that the expression levels of microRNA-455 and microRNA-125a in tissues of adenoma patients and cancer patients. (C) BCL2L12 expression in adenoma patients and colorectal carcinoma patients.

Screening of Small Molecule Drugs

By using CMap software, Top ten molecular drugs significantly correlated molecules with BCL2L12 were achieved including navitoclax, maraviroc, HIC, and imiquimod. The result is shown in Table 4. The molecular chemical structure diagram of navitoclax is also shown in Figure 8.

Table 4 Results of Connectivity Map Analysis

Figure 8 Chemical structure diagram of navitoclax.

Discussion

Colorectal adenoma is closely related to colorectal cancer. Studies suggest that at least 80% of colorectal cancer evolved from colorectal adenoma, which lasted for more than 5 years, with an average of 10–15 years. Early detection and complete endoscopic resection of adenoma are the keys to survival, and there is almost no chance of cancer development.15–17 In this present study, GSE41655 and GSE57965 were obtained from the GEO. Five hundred and sixty-nine colorectal carcinoma patients’ information was downloaded from the TCGA database. Twenty-seven CRC-related microRNAs were identified and may play a vital role in promoting the development of adenomas into adenocarcinomas. To further study the regulatory mechanism of the 27 microRNAs in colorectal cancer, we used FunRich for further study. Gene ontology enrichment analysis showed that these microRNAs were primarily associated with the nucleus, transcription factor activity, cell communication, and signal transduction. This is consistent with the recognition that cell cycle and cell proliferation regulator function defects are the key reasons for cancer occurrence and development.18,19 The ion transport in tumor cells is essentially different from that in normal cells.20 KEGG enrichment pathway analysis showed that these selected microRNAs were mainly enriched in 6 pathways including glycosphingolipid biosynthesis, retinol metabolism, AGE-RAGE signaling pathway in diabetic complications, pancreatic secretion, inflammatory mediator regulation of TRP channels, and renin secretion. Several studies on cancer-related glycosylation have shown that abnormal glycosylation is a common feature of various stages of malignant transformation and tumor progression.21,22 It is important to note that the glycosylation changes observed so far are relatively specific to the type and stage of cancer, making glycans potential biomarkers and targets for drug therapy. As for the AGE-RAGE signaling pathway in diabetic complications, previous studies have shown that glycine may protect diabetic macrovascular complications by improving Glo1 function, inhibiting AGE/RAGE pathway and subsequent oxidative stress.23 TRP channels are regulated by proinflammatory mediators, neuropeptides, and cytokines. Antagonists or agonists targeting these receptors have made great progress in the treatment of pain.24 However, the regulation mechanism of these pathways in colorectal cancer has not been discovered yet.

The microRNA-mRNA regulatory network was constructed based on FunRich and Cytoscape. Six DEMs (microRNA-199a, microRNA-199b, microRNA-127, microRNA-125a, microRNA-455, and microRNA-146b) and 34 potential genes were selected. Previous studies have shown that they are involved in the development of cancer. For instance, a previous study reported down-regulation of microRNA-199a-5p in NSCLC altered the expression level of GRP78 and spliced XBP1 through UPR pathway.25 As for hsa-miR-127, a previous study reported it affects the survival of glioma patients by regulating replication initiation factor 1 and promotes the proliferation, migration, and invasion of tumor cells. The hsa-mir-127/REPIN1 pathway is involved in gliomas and could be a potential therapeutic target.26 Onnis et al also demonstrated that combined expression of Epstein–Barr nuclear antigen 1 and microRNA-127 affects the expression of major B cell regulatory factors in memory B cells.27 As for hsa-miR-125a, a previous study reported upregulated microRNA-125a obviously enhanced cell proliferation, migration and invasion in HSCC, with upregulation of C–C Chemokine Receptor Type 7.28 It also activates p53 and Induces Apoptosis in Lung Cancer Cells.29 Besides, microRNA-125a may play an anti-cancer role by regulating BRCA1 signaling pathway, and reintroduction of microRNA-125a analogues may be a potential adjuvant therapy for advanced/chemoresistant breast carcinoma.30 Survival analysis indicated that overexpression of microRNA-199a, microRNA-199b, microRNA-127, and microRNA-125a were associated with poor OS in patients with colorectal cancer by utilizing km-plot software.

Moreover, 34 potential genes including AQP11, ATP13A2, KLHL23, BCL2L12, DPF2, EVA1A, FOXQ1 were further investigated in our research. KLHL23 (Kelch-like Family Member 23) has been studied a lot in recent decades. Overexpression of KLHL23 may affect the treatment of gastric cancer.31 DPF2 (Double plant homeodomain finger 2) is a highly conserved member of the D4 protein family, which is widely expressed in human tissues. Recently, it has been found that DPF2 can inhibit myeloid differentiation of hematopoietic stem/progenitor cells and acute myeloid leukemia cells.32 FOXQ1, forkhead box q1, recent studies have found that FOXQ1 and other gene changes may be involved in the mechanism of secondary resistance to anti EGFR antibody therapy in colorectal carcinoma metastasis.33 It has also been found to be expressed in breast cancer, ovarian cancer, and endometrial cancer, in which the receptor, through its ligand, is produced by tumor cells or matrix elements to stimulate the invasion and metastasis of tumor cells.34 Besides, another study also reported that FOXQ1 is up-regulated in colorectal cancer cells and clinical specimens, which may be of great significance for the diagnosis of colorectal carcinoma.

Besides, a 7-microRNAs signature (microRNA-125a-5p, microRNA-377, microRNA-376c, microRNA-455-3p, microRNA-126, microRNA-99a, and microRNA-193b) of 569 patients with colorectal cancer were constructed by using single variable and multivariate Cox and risk scoring methods. CRC patients were divided into a high-risk group and a low-risk group. The results showed that the signal had good repeatability and robustness in predicting the prognosis of CRC patients.

According to the combination of the GEO and TCGA analysis results, 2 microRNAs (microRNA-455 and microRNA-125a) were identified for further research. After using qRT-PCR analysis, we found that microRNA-455 and microRNA-125a were differently expressed in the tissue of adenoma patients and cancer patients, and we also found that microRNA-455 and microRNA-125a had strong effects, which all target BCL2L12. By using IHC, we also found it has a positive strong expression in colorectal cancer tissues and a negative weak expression in adenoma tissues. The BCL2L12 gene, a member of the BCL2 family, was discovered and cloned in 2001.35

The BCL2L12 gene is located on chromosome 19q13.3, between IRF3 and PRMT1/hrmt1l2. It is composed of 7 coding exons and 6 inserted introns. It spans 8.8kb genome region and encodes a proline rich protein with a BH2 and a putative BH3 domain. It is obviously related to apoptosis. BCL2L12 mRNA was expressed in many tissues.36 In addition to the BH2 domain, BCL2L12 protein contains five consistent PXXP tetrapeptide motifs and a proline rich region, indicating that it interacts with protein tyrosine kinases such as phospholipase, RAS guanosine triphosphatase activating protein and SRC like tyrosine kinases.37 The prognostic significance of BCL2L12 expression in several tumor types has been assessed.38,39 Many preclinical epithelial ovarian cancer studies correlate Bcl-2-regulated apoptosis to metformin’s chemo-sensitizing effects. The chemo-sensitizing effect of metformin seems to be correlated with p53 function. In the presence of p53, metformin suppresses hexokinase II (glycolytic enzyme) and pyruvate dehydrogenase kinase (anti-apoptotic serine/threonine kinase). As a result, epithelial ovarian cancer cells are sensitized to metformin.40 Besides, BCL2L12 protein expression could be utilized as a significant prognostic tissue biomarker in patients with primary advanced-stage Laryngeal squamous cell carcinoma.41 For colorectal cancer, the expression of BCL2L12 mRNA gradually increased during the treatment, similar to the expression of other BCL-2 family genes that are beneficial to apoptosis or specific pro-apoptotic transcripts, so it is suggested that BCL2L12 can promote apoptosis in colorectal cancer cells treated with chemotherapy.42 By using CMap software, several small molecules with potential anti-BCL2L12 effects were identified. Navitoclax, a targeted high-affinity inhibitor of BCL-2, was found to enhance the activity of AML patients combined with other chemotherapy drugs in early clinical trials.43 However, researches about navitoclax Inhibiting the transformation of adenoma to adenocarcinoma has not been performed yet, and more clinical trials are needed to understand this relationship.

At present, with the development of cancer treatment, individual difference therapy has been paid more and more attention. Therefore, it is important to discover novel therapeutic targets and methods for cancer patients. Our conclusions demonstrated that many DE mRNAs and microRNAs taken part in the transformation from adenoma to adenocarcinoma and had prognostic worth. Because all of our data are obtained from the GEO and TCGA database through R software, further data analysis and basic experiments are needed to verify.

Conclusion

Our research summarized some mechanisms of the progression of colorectal cancer. Plenty of DEMs and DEGs were selected between adenomas and adenocarcinomas by utilizing bioinformatics methods. Also, microRNA-455, microRNA-125a and BCL2L12 were identified as potential targets for inhibiting the transformation of adenoma to adenocarcinoma. However, these findings may be considered for future investigation.

Data Sharing Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethics Approval and Consent to Participate

The present study was conducted by the Declaration of Helsinki. The Ethics Committee of the Second Affiliated Hospital of Soochow University approved the study. The participants approved the use of clinical samples by providing written informed consent.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Disclosure

The authors have declared that no competing interest exists.

References

1. Steele SR, Chang GJ, Hendren S, et al. Practice guideline for the surveillance of patients after curative treatment of colon and rectal cancer. Dis Colon Rectum. 2015;58(8):713–725. doi:10.1097/DCR.0000000000000410

2. Vogel JD, Eskicioglu C, Weiser MR, Feingold DL, Steele SR. The American society of colon and rectal surgeons clinical practice guidelines for the treatment of colon cancer. Dis Colon Rectum. 2017;60(10):999–1017. doi:10.1097/DCR.0000000000000926

3. Balch C, Ramapuram JB, Tiwari AK. The Epigenomics of embryonic pathway signaling in colorectal cancer. Front Pharmacol. 2017;8:267. doi:10.3389/fphar.2017.00267

4. Zarkavelis G, Boussios S, Papadaki A, Katsanos KH, Christodoulou DK, Pentheroudakis G. Current and future biomarkers in colorectal cancer. Ann Gastroenterol. 2017;30(6):613–621.

5. Ma Y, Pan X, Xu P, et al. Plasma microRNA alterations between EGFR-activating mutational NSCLC patients with and without primary resistance to TKI. Oncotarget. 2017;8(51):88529–88536. doi:10.18632/oncotarget.19874

6. Ren P, Gong F, Zhang Y, Jiang J, Zhang H. MicroRNA-92a promotes growth, metastasis, and chemoresistance in non-small cell lung cancer cells by targeting PTEN. Tumour Biol. 2016;37(3):3215–3225. doi:10.1007/s13277-015-4150-3

7. Sin TK, Wang F, Meng F, et al. Implications of MicroRNAs in the treatment of gefitinib-resistant non-small cell lung cancer. Int J Mol Sci. 2016;17(2):237. doi:10.3390/ijms17020237

8. Ma X, Liang AL, Liu YJ. Research progress on the relationship between lung cancer drug-resistance and microRNAs. J Cancer. 2019;10(27):6865–6875. doi:10.7150/jca.31952

9. Kania EE, Carvajal-Moreno J, Hernandez VA, et al. hsa-miR-9-3p and hsa-miR-9-5p as Post-Transcriptional Modulators of DNA Topoisomerase IIalpha in Human Leukemia K562 Cells with Acquired Resistance to Etoposide. Mol Pharmacol. 2020;97(3):159–170.

10. Boussios S, Ozturk MA, Moschetta M, et al. The developing story of predictive biomarkers in colorectal cancer. J Pers Med. 2019;9(1):12. doi:10.3390/jpm9010012

11. Li L, Wang G, Li N, Yu H, Si J, Wang J. Identification of key genes and pathways associated with obesity in children. Exp Ther Med. 2017;14(2):1065–1073. doi:10.3892/etm.2017.4597

12. Ikeda N, Nakajima Y, Tokuhara T, et al. Clinical significance of aminopeptidase N/CD13 expression in human pancreatic carcinoma. Clin Cancer Res. 2003;9(4):1503–1508.

13. Veenman CJ, Tax DM. LESS: a model-based classifier for sparse subspaces. IEEE Trans Pattern Anal Mach Intell. 2005;27(9):1496–1500. doi:10.1109/TPAMI.2005.182

14. Pollock BE, Storlie CB, Link MJ, Stafford SL, Garces YI, Foote RL. Comparative analysis of arteriovenous malformation grading scales in predicting outcomes after stereotactic radiosurgery. J Neurosurg. 2017;126(3):852–858. doi:10.3171/2015.11.JNS151300

15. Kitahara O, Furukawa Y, Tanaka T, et al. Alterations of gene expression during colorectal carcinogenesis revealed by cDNA microarrays after laser-capture microdissection of tumor tissues and normal epithelia. Cancer Res. 2001;61(9):3544–3549.

16. Lechner S, Muller-Ladner U, Renke B, Scholmerich J, Ruschoff J, Kullmann F. Gene expression pattern of laser microdissected colonic crypts of adenomas with low grade dysplasia. Gut. 2003;52(8):1148–1153. doi:10.1136/gut.52.8.1148

17. Notterman DA, Alon U, Sierk AJ, Levine AJ. Transcriptional gene expression profiles of colorectal adenoma, adenocarcinoma, and normal tissue examined by oligonucleotide arrays. Cancer Res. 2001;61(7):3124–3130.

18. Perez R, Wu N, Klipfel AA, Beart RW. A better cell cycle target for gene therapy of colorectal cancer: cyclin G. J Gastrointest Surg. 2003;7(7):884–889. doi:10.1007/s11605-003-0034-8

19. Tominaga O, Nita ME, Nagawa H, Fujii S, Tsuruo T, Muto T. Expressions of cell cycle regulators in human colorectal cancer cell lines. Jpn J Cancer Res. 1997;88(9):855–860. doi:10.1111/j.1349-7006.1997.tb00461.x

20. Djamgoz MB, Coombes RC, Schwab A. Ion transport and cancer: from initiation to metastasis. Philos Trans R Soc Lond B Biol Sci. 2014;369(1638):20130092. doi:10.1098/rstb.2013.0092

21. Hakomori S. Tumor-associated carbohydrate antigens defining tumor malignancy: basis for development of anti-cancer vaccines. Adv Exp Med Biol. 2001;491:369–402.

22. Yue T, Goldstein IJ, Hollingsworth MA, Kaul K, Brand RE, Haab BB. The prevalence and nature of glycan alterations on specific proteins in pancreatic cancer patients revealed using antibody-lectin sandwich arrays. Mol Cell Proteomics. 2009;8(7):1697–1707. doi:10.1074/mcp.M900135-MCP200

23. Wang Z, Zhang J, Chen L, Li J, Zhang H, Guo X. Glycine Suppresses AGE/RAGE Signaling pathway and subsequent oxidative stress by restoring Glo1 function in the aorta of diabetic rats and in HUVECs. Oxid Med Cell Longev. 2019;2019:4628962.

24. Premkumar LS, Abooj M. TRP channels and analgesia. Life Sci. 2013;92(8–9):415–424. doi:10.1016/j.lfs.2012.08.010

25. Ahmadi A, Khansarinejad B, Hosseinkhani S, Ghanei M, Mowla SJ. miR-199a-5p and miR-495 target GRP78 within UPR pathway of lung cancer. Gene. 2017;620:15–22. doi:10.1016/j.gene.2017.03.032

26. Wang Y, Lin Y. Hsa-mir-127 impairs survival of patients with glioma and promotes proliferation, migration and invasion of cancerous cells by modulating replication initiator 1. Neuroreport. 2018;29(14):1166–1173. doi:10.1097/WNR.0000000000001089

27. Onnis A, Navari M, Antonicelli G, et al. Epstein-Barr nuclear antigen 1 induces expression of the cellular microRNA hsa-miR-127 and impairing B-cell differentiation in EBV-infected memory B cells. New insights into the pathogenesis of Burkitt lymphoma. Blood Cancer J. 2012;2:e84. doi:10.1038/bcj.2012.29

28. Jin S, Liu MD, Wu H, et al. Overexpression of hsa-miR-125a-5p enhances proliferation, migration and invasion of head and neck squamous cell carcinoma cell lines by upregulating C-C chemokine receptor type 7. Oncol Lett. 2018;15(6):9703–9710.

29. Jiang L, Chang J, Zhang Q, Sun L, Qiu X. MicroRNA hsa-miR-125a-3p activates p53 and induces apoptosis in lung cancer cells. Cancer Invest. 2013;31(8):538–544. doi:10.3109/07357907.2013.820314

30. Xu X, Lv YG, Yan CY, Yi J, Ling R. Enforced expression of hsa-miR-125a-3p in breast cancer cells potentiates docetaxel sensitivity via modulation of BRCA1 signaling. Biochem Biophys Res Commun. 2016;479(4):893–900. doi:10.1016/j.bbrc.2016.09.087

31. Choi ES, Lee H, Lee CH, Goh SH. Overexpression of KLHL23 protein from read-through transcription of PHOSPHO2-KLHL23 in gastric cancer increases cell proliferation. FEBS Open Bio. 2016;6(11):1155–1164. doi:10.1002/2211-5463.12136

32. Huber FM, Greenblatt SM, Davenport AM, et al. Histone-binding of DPF2 mediates its repressive role in myeloid differentiation. Proc Natl Acad Sci U S A. 2017;114(23):6016–6021. doi:10.1073/pnas.1700328114

33. Abba M, Patil N, Rasheed K, et al. Unraveling the role of FOXQ1 in colorectal cancer metastasis. Mol Cancer Res. 2013;11(9):1017–1028. doi:10.1158/1541-7786.MCR-13-0024

34. Peng X, Luo Z, Kang Q, et al. FOXQ1 mediates the crosstalk between TGF-beta and Wnt signaling pathways in the progression of colorectal cancer. Cancer Biol Ther. 2015;16(7):1099–1109. doi:10.1080/15384047.2015.1047568

35. Scorilas A, Kyriakopoulou L, Yousef GM, Ashworth LK, Kwamie A, Diamandis EP. Molecular cloning, physical mapping, and expression analysis of a novel gene, BCL2L12, encoding a proline-rich protein with a highly conserved BH2 domain of the Bcl-2 family. Genomics. 2001;72(2):217–221. doi:10.1006/geno.2000.6455

36. Scorilas A, Kyriakopoulou L, Yousef GM, Ashworth LK, Kwamie A, Diamandis EP. Molecular cloning, physical mapping, and expression analysis of a novel gene, BCL2L12, encoding a proline-rich protein with a highly conserved BH2 domain of the Bcl-2 family. Genomics. 2001;72(2):217–221.

37. Koch CA, Anderson D, Moran MF, Ellis C, Pawson T. SH2 and SH3 domains: elements that control interactions of cytoplasmic signaling proteins. Science. 1991;252(5006):668–674. doi:10.1126/science.1708916

38. Florou D, Papadopoulos IN, Scorilas A. Molecular analysis and prognostic impact of the novel apoptotic gene BCL2L12 in gastric cancer. Biochem Biophys Res Commun. 2010;391(1):214–218. doi:10.1016/j.bbrc.2009.11.034

39. Talieri M, Diamandis EP, Katsaros N, Gourgiotis D, Scorilas A. Expression of BCL2L12, a new member of apoptosis-related genes, in breast tumors. Thromb Haemost. 2003;89(6):1081–1088. doi:10.1055/s-0037-1613411

40. Boussios S, Mikropoulos C, Samartzis E, et al. Wise Management of Ovarian Cancer: on the Cutting Edge. J Pers Med. 2020;10(2):41. doi:10.3390/jpm10020041

41. Giotakis AI, Lazaris AC, Kataki A, Kontos CK, Giotakis EI. Positive BCL2L12 expression predicts favorable prognosis in patients with laryngeal squamous cell carcinoma. Cancer Biomark. 2019;25(2):141–149. doi:10.3233/CBM-181772

42. Kontos CK, Avgeris M, Vassilacopoulou D, Ardavanis A, Scorilas A. Molecular Effects of Treatment of Human Colorectal Cancer Cells with Natural and Classical Chemotherapeutic Drugs: alterations in the Expression of Apoptosis-related BCL2 Family Members, Including BCL2L12. Curr Pharm Biotechnol. 2018;19(13):1064–1075. doi:10.2174/1389201019666181112101410

43. Kivioja JL, Thanasopoulou A, Kumar A, et al. Dasatinib and navitoclax act synergistically to target NUP98-NSD1(+)/FLT3-ITD(+) acute myeloid leukemia. Leukemia. 2019;33(6):1360–1372. doi:10.1038/s41375-018-0327-2

Creative Commons License © 2021 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.