Genome-wide mRNA and miRNA expression profiling reveal multiple regulatory networks in colorectal cancer

CMC: Cancer miRNA Census - a list of cancer-related miRNA genes

A growing body of evidence indicates an important role of miRNAs in cancer; however, there is no definitive, convenient-to-use list of cancer-related miRNAs or miRNA genes that may serve as a reference for analyses of miRNAs in cancer. To this end, we created a list of 165 cancer-related miRNA genes called the Cancer miRNA Census (CMC). The list is based on a score, built on various types of functional and genetic evidence for the role of particular miRNAs in cancer, e.g. miRNA-cancer associations reported in databases, associations of miRNAs with cancer hallmarks, or signals of positive selection of genetic alterations in cancer. The presence of well-recognized cancer-related miRNA genes, such as MIR21, MIR155, MIR15A, MIR17 or MIRLET7s, at the top of the CMC ranking directly confirms the accuracy and robustness of the list. Additionally, to verify and indicate the reliability of CMC, we performed a validation of criteria used to build CMC, comparison of CMC with various cancer data (publications and databases), and enrichment analyses of biological pathways and processes such as Gene Ontology or DisGeNET. All validation steps showed a strong association of CMC with cancer/cancer-related processes confirming its usefulness as a reference list of miRNA genes associated with cancer.

Therapeutic targeting of the TPX2/TTK network in colorectal cancer

BACKGROUND

While the increased screening, changes in lifestyle, and recent advances in treatment regimen have decreased colorectal cancer (CRC) mortality, metastatic disease and recurrence remains a major clinical challenge. In the era of precision medicine, the identification of actionable novel therapeutic targets could ultimately offer an alternative treatment strategy for CRC.

METHODS

RNA-Seq was conducted using the illumina platform, while bioinformatics analyses were conducted using CLC genomics workbench and iDEP.951. Colony forming unit, flow cytometry, and fluorescent microscopy were used to assess cell proliferation, cell cycle distribution, and cell death, respectively. The growth potential of CRC cells under 3-dimensional (3D) conditions was assessed using Matrigel. STRING database (v11.5) and Ingenuity Pathway Analysis (IPA) tool were used for network and pathway analyses. CRISPR-Cas9 perturbational effects database was used to identify potential therapeutic targets for CRC, through integration with gene-drug interaction database. Structural modeling and molecular docking were used to assess the interaction between candidate drugs and their targets.

RESULTS

In the current study, we investigated the therapeutic potential of targeting TPX2, TTK, DDX39A, and LRP8, commonly upregulated genes in CRC identified through differential expression analysis in CRC and adjacent non-cancerous tissue. Targeted depletion of TPX2 and TTK impaired CRC proliferation, cell cycle progression, and organoid formation under 3D culture conditions, while suppression of DDX39A and LRP8 had modest effects on CRC colony formation. Differential expression analysis and bioinformatics on TPX2 and TTK-deficient cells identified cell cycle regulation as the hallmark associated with loss of TPX2 and TTK. Elevated expression of TPX2 and TTK correlated with an oncogenic state in tumor tissue from patients with colon adenocarcinoma, thus corroborating an oncogenic role for the TPX2/TTK network in the pathogenesis of CRC. Gene set enrichment and pathway analysis of TPX2high/TTKhigh CRC identified numerous additional gene targets as integral components of the TPX2/TTK network. Integration of TPX2/TTK enriched network with CRISPR-Cas9 functional screen data identified numerous novel dependencies for CRC. Additionally, gene-drug interaction analysis identified several druggable gene targets enriched in the TPX2/TTK network, including AURKA, TOP2A, CDK1, BIRC5, and many others.

CONCLUSIONS

Our data has implicated an essential role for TPX2 and TTK in CRC pathogenesis and identified numerous potential therapeutic targets and their drug interactions, suggesting their potential clinical use as a novel therapeutic strategy for patients with CRC. Video Abstract.

High Expression of EIF4G2 Mediated by the TUG1/Hsa-miR-26a-5p Axis Is Associated with Poor Prognosis and Immune Infiltration of Gastric Cancer

Objective

Eukaryotic translation initiation factor 4 gamma 2 (EIF4G2) is involved in the occurrence and development of various tumors. However, the effect of EIF4G2 in gastric cancer (GC) has not been fully explored. The purpose of this study was to explore the function and mechanism of EIF4G2 in GC.

Methods

The Tumor Immune Estimation Resource 2.0 database was used to analyze EIF4G2 expression in various cancers and the relationship between EIF4G2 expression and tumor-infiltrating immune cells. Gene Expression Profiling Interactive Analysis was utilized to assess the EIF4G2 expression level and its effect on survival in GC. UALCAN was conducted to analyze EIF4G2 expression in various subgroups of GC. The Kaplan–Meier plotter was employed for survival analysis. Receiver operator characteristic (ROC) curve analysis was applied to evaluate the diagnostic role of EIF4G2 in GC. LinkedOmics was used to identify the co-expressed genes and Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. The Tumor-Immune System Interaction database was employed to analyze the correlation between EIF4G2 expression and tumor-infiltrating lymphocytes. The starBase web platform was used to predict the upstream microRNAs and long noncoding RNAs.

Results

EIF4G2 expression was upregulated in GC tissues compared to normal controls. High expression of EIF4G2 indicated poor prognosis in GC. ROC analysis revealed that EIF4G2 had good diagnostic ability to distinguish GC from normal tissues. Immune infiltration analysis indicated that EIF4G2 expression may be involved in the modulation of tumor immune infiltration in GC. Finally, we determined that the Taurine Upregulated 1 (TUG1)/hsa-miR-26a-5p/EIF4G2 axis was the most likely regulatory pathway involved in GC development.

Conclusions

EIF4G2 was upregulated in GC and elevated expression of EIF4G2 indicated unfavorable prognosis. Moreover, EIF4G2 expression may be involved in the regulation of tumor immune cell infiltration. The TUG1/hsa-miR-26a-5p axis is a likely upstream regulatory mechanism of EIF4G2 in GC. EIF4G2 may thus serve as a prognosis biomarker and present a new therapeutic target.

Transcriptome profiling and co-expression network analysis of lncRNAs and mRNAs in colorectal cancer by RNA sequencing

Background

Long non-coding RNAs (lncRNAs) are widely involved in the pathogenesis of cancers. However, biological roles of lncRNAs in occurrence and progression of colorectal cancer (CRC) remain unclear. The current study aimed to evaluate the expression pattern of lncRNAs and messenger RNAs (mRNAs).

Methods

RNA sequencing (RNA-Seq) in CRC tissues and adjacent normal tissues from 6 CRC patients was performed and functional lncRNA-mRNA co-expression network was constructed afterwards. Gene enrichment analysis was demonstrated using DAVID 6.8 tool. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was used to validate the expression pattern of differentially expressed lncRNAs. Pearson correlation analysis was applied to evaluate the relationships between selected lncRNAs and mRNAs.

Results

One thousand seven hundred and sixteenth differentially expressed mRNAs and 311 differentially expressed lncRNAs were screened out. Among these, 568 mRNAs were up-regulated while 1148 mRNAs down-regulated, similarly 125 lncRNAs were up-regulated and 186 lncRNAs down-regulated. In addition, 1448 lncRNA–mRNA co-expression pairs were screened out from 940,905 candidate lncRNA-mRNA pairs. Gene enrichment analysis revealed that these lncRNA-related mRNAs are associated with cell adhesion, collagen adhesion, cell differentiation, and mainly enriched in ECM-receptor interaction and PI3K-Akt signaling pathways. Finally, RT-qPCR results verified the expression pattern of lncRNAs, as well as the relationships between lncRNAs and mRNAs in 60 pairs of CRC tissues.

Conclusions

In conclusion, these results of the RNA-seq and bioinformatic analysis strongly suggested that the dysregulation of lncRNA is involved in the complicated process of CRC development, and providing important insight regarding the lncRNAs involved in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09878-6.

mRNA-miRNA bipartite networks reconstruction in different tissues of bladder cancer based on gene co-expression network analysis

Bladder cancer (BC) is one of the most important cancers worldwide, and if it is diagnosed early, its progression in humans can be prevented and long-term survival will be achieved accordingly. This study aimed to identify novel micro-RNA (miRNA) and gene-based biomarkers for diagnosing BC. The microarray dataset of BC tissues (GSE13507) listed in the GEO database was analyzed for this purpose. The gene expression data from three BC tissues including 165 primary bladder cancer (PBC), 58 normal looking-bladder mucosae surrounding cancer (NBMSC), and 23 recurrent non-muscle invasive tumor tissues (RNIT) were used to reconstruct gene co-expression networks. After preprocessing and normalization, deferentially expressed genes (DEGs) were obtained and used to construct the weighted gene co-expression network (WGCNA). Gene co-expression modules and low-preserved modules were extracted among BC tissues using network clustering. Next, the experimentally validated mRNA-miRNA interaction information were used to reconstruct three mRNA-miRNA bipartite networks. Reactome pathway database and Gene ontology (GO) was subsequently performed for the extracted genes of three bipartite networks and miRNAs, respectively. To further analyze the data, ten hub miRNAs (miRNAs with the highest degree) were selected in each bipartite network to reconstruct three bipartite subnetworks. Finally, the obtained biomarkers were comprehensively investigated and discussed in authentic studies. The obtained results from our study indicated a group of genes including PPARD, CST4, CSNK1E, PTPN14, ETV6, and ADRM1 as well as novel miRNAs (e.g., miR-16-5p, miR-335-5p, miR-124-3p, and let-7b-5p) which might be potentially associated with BC and could be a potential biomarker. Afterward, three drug-gene interaction networks were reconstructed to explore candidate drugs for the treatment of BC. The hub miRNAs in the mRNA-miRNA bipartite network played a fundamental role in BC progression; however, these findings need further investigation.

The role of miR-370 and miR-138 in the regulation of BMP2 suppressor gene expression in colorectal cancer: preliminary studies

PURPOSE

Colorectal cancer (CRC) is the fourth-most common cancer worldwide and the second most common cancer cause of death in the world. The components of the TGFβ-signalling pathway, which are often affected by miRNAs, are involved in the regulation of apoptosis and cell cycle. Therefore, in the current study, the expression of BMP2 gene in CRC tissues at different clinical stages compared to the non-tumour tissues has been assessed. Moreover, the plasma BMP2 protein concentration in the same group of CRC patients has been validated. Due to the constant necessity to conduct further research of the correlation between specific miRNAs and mRNAs in CRC, in silico analysis has been performed to select miRNAs that regulate BMP2 mRNA.

METHODS

The cDNA samples from tumor and non-tumor tissue were used in a qPCR reaction to determine the mRNA expression of the BMP2 gene and the expression of selected miRNAs. The concentration of BMP2 protein in plasma samples was also measured.

RESULTS

It was indicated that BMP2 was downregulated in CRC tissue. Moreover, miR-370 and miR-138 expression showed an upward trend. Decreased BMP2 with accompanied increasing miR-370 and miR-138 expression was relevant to the malignant clinicopathological features of CRC and consequently poor patient prognosis.

CONCLUSION

Our data suggest that miR-370 with its clear expression in plasma samples may be a potential diagnostic marker to determine the severity of the disease in patients at a later stage of colorectal cancer.

Interplay between K-RAS and miRNAs

K-RAS is frequently mutated in cancers, and its overactivation can lead to oncogene-induced senescence (OIS), a barrier to cellular transformation. Feedback onto K-RAS limits its signaling to avoid senescence while achieving the appropriate level of activation that promotes proliferation and survival. Such regulation could be mediated by miRNAs, as aberrant RAS signaling and miRNA activity coexist in several cancers, with miRNAs acting both up- and downstream of K-RAS. Several miRNAs both regulate and are regulated by K-RAS, suggesting a noncoding RNA-based feedback mechanism. Functional interactions between K-RAS and the miRNA machinery have also begun to unfold. This review comprehensively surveys the state of knowledge connecting K-RAS to miRNA function and proposes a model for the regulation of K-RAS signaling by noncoding RNAs.

miR-877-3p inhibits tumor growth and angiogenesis of osteosarcoma through Fibroblast Growth Factor 2 signaling

Osteosarcoma (OS) is the most common high-grade malignant bone tumor in teenagers. MicroRNAs can function as posttranscriptional regulators of gene expression, playing critical roles in cancer dev-877-3p in OS. Quantitative real-time RT-PCR was carried out for detecting miR-877-3p expression in OS. The effects of miR-877-3p on proliferation was analyzed via MTT, colony formation, and flow cytometry assays. Angiogenesis of endothelial cells were investigated by wound healing and tube formation assay. Gene profiling based on PCR array and luciferase reporter assay were conducted to determine target genes of miR-877-3p. In-vivo study was used to determine the effects of miR-877-3p on the tumor growth. The expression of miR-877-3p was markedly downregulated in OS tissues and cell lines. Low expression of miR-877-3p predicts poor prognosis of OS patients. miR-877-3p overexpression was found to inhibit the proliferation of OS cell lines. The angiogenesis assays showed that miR-877-3p attenuated the angiogenesis. Further mechanism studies showed that miR-877-3p can reduce (Fibroblast Growth Factor 2) FGF2 expression in OS cells by binding to the 3’UTR end of FGF2. Moreover, increased expression of miR-877-3p was responsible for the inhibition of tumor growth and angiogenesis. Taken together, our findings indicated that miR-877-3p might exhibit tumor suppressive role by targeting FGF2 signaling, which may serve as potential target for OS.

Identification of potential genes associated with ALDH1A1 overexpression and cyclophosphamide resistance in chronic myelogenous leukemia using network analysis

Cyclophosphamide (CP), an important alkylating agent which is used in the treatment therapy for chronic myeloid leukemia (CML). However, acquired drug resistance owing to the inactivation of its active metabolite aldophosphamide via tumoral-overexpressing aldehyde dehydrogenase (ALDH1A1) is one of the major issues with the CP therapy. However, the underlying mechanism of ALDH1A1 overexpression in cancer cells remains poorly defined. Therefore, the current study focused on analyzing the ALDH1A1-overexpressing microarray data for CP resistance and CP-sensitive CML cell lines. In this study, the microarray dataset was obtained from Gene Expression Omnibus GEO. The GEO2R tool was used to identify Differentially Expressing Genes (DEGs). Further, protein-protein interaction (PPI) network of DEGs were constructed using STRING database. Finally, Hub gene-miRNA-TFs interaction were constructed using miRNet tool. A total of 749 DEGs including 387 upregulated and 225 downregulated genes were identified from this pool of microarray data. The construction of DEGs network resulted in identification of three genes including ZEB2, EZH2, and MUC1 were found to be majorly responsible for ALDH1A1 overexpression. miRNA analysis identified that, hsa-mir-16-5p and hsa-mir-26a-5p as hub miRNA which are commonly interacting with maximum target genes. Additionally, drug-gene interaction analysis was performed to identify drugs which are responsible for ALDH1A1 expression. The entire study may provide a deeper understanding about ALDH1A1 regulatory genes responsible for its overexpression in CP resistance cancer. This understanding may be further explore for developing possible co-therapy to avoid the ALDH1A1-mediated CP resistance.

miR-26a enhances colorectal cancer cell growth by targeting RREB1 deacetylation to activate AKT-mediated glycolysis

We previously reported the inhibitory effects of microRNA-26a (miR-26a) on the conversion of pyruvate to acetyl coenzyme A in glucose metabolism by directly targeting pyruvate dehydrogenase protein X component in colorectal cancer (CRC) cells (Chen B et al., BMC Cancer 2014). Here, using microRNA in situ hybridization, we confirmed that miR-26a levels were elevated in 77 human CRC tissue samples and further investigated the key miR-26a-mediated metabolic regulation elements and signaling pathways in CRC cells through quantitative proteomic dissection combined with cancer cell biology and biochemical loss-of-function analysis. We found that AKT transcription signaling was a target pathway via miR-26a-mediated deacetylation modification of Ras-responsive element-binding protein 1 (RREB1) at the Lys-60 residue. miR-26a improved the deacetylation level of RREB1, thus contributing to RREB1 binding to the AKT1 promoter to activate AKT transcription and its related signaling pathway in glycolysis. Moreover, miR-26a promoted CRC tumorigenesis in CRC cells and subcutaneous xenograft mice. Thus, miR-26a is a key regulator of CRC tumorigenesis that mediates the deacetylation modification of RREB1 to enhance AKT1 transcription and downstream target gene expression in glycolysis for CRC growth.

Integrated whole transcriptome and small RNA analysis revealed multiple regulatory networks in colorectal cancer

Colorectal cancer (CRC) remains a global disease burden and a leading cause of cancer related deaths worldwide. The identification of aberrantly expressed messenger RNA (mRNA), long non-coding RNA (lncRNA), and microRNA (miRNA), and the resulting molecular interactions and signaling networks is essential for better understanding of CRC, identification of novel diagnostic biomarkers and potential development of therapeutic interventions. Herein, we performed microRNA (miRNA) sequencing on fifteen CRC and their non-tumor adjacent tissues and whole transcriptome RNA-Seq on six paired samples from the same cohort and identified alterations in miRNA, mRNA, and lncRNA expression. Computational analyses using Ingenuity Pathway Analysis (IPA) identified multiple activated signaling networks in CRC, including ERBB2, RABL6, FOXM1, and NFKB networks, while functional annotation highlighted activation of cell proliferation and migration as the hallmark of CRC. IPA in combination with in silico prediction algorithms and experimentally validated databases gave insight into the complex associations and interactions between downregulated miRNAs and upregulated mRNAs in CRC and vice versa. Additionally, potential interaction between differentially expressed lncRNAs such as H19, SNHG5, and GATA2-AS1 with multiple miRNAs has been revealed. Taken together, our data provides thorough analysis of dysregulated protein-coding and non-coding RNAs in CRC highlighting numerous associations and regulatory networks thus providing better understanding of CRC.

miRNA Clusters with Up-Regulated Expression in Colorectal Cancer

Simple Summary

As miRNAs show the capacity to be used as CRC biomarkers, we analysed experimentally validated data about frequently up-regulated miRNA clusters in CRC tissue. We identified 15 clusters that showed increased expression in CRC: miR-106a/363, miR-106b/93/25, miR-17/92a-1, miR-181a-1/181b-1, miR-181a-2/181b-2, miR-181c/181d, miR-183/96/182, miR-191/425, miR-200c/141, miR-203a/203b, miR-222/221, mir-23a/27a/24-2, mir-29b-1/29a, mir-301b/130b and mir-452/224. Cluster positions in the genome are intronic or intergenic. Most clusters are regulated by several transcription factors, and by long non-coding RNAs. In some cases, co-expression of miRNA with other cluster members or host gene has been proven. miRNA expression patterns in cancer tissue, blood and faeces were compared. The members of the selected clusters target 181 genes. Their functions and corresponding pathways were revealed with the use of Panther analysis. Clusters miR-17/92a-1, miR-106a/363, miR-106b/93/25 and miR-183/96/182 showed the strongest association with metastasis occurrence and poor patient survival, implicating them as the most promising targets of translational research.

Abstract

Colorectal cancer (CRC) is one of the most common malignancies in Europe and North America. Early diagnosis is a key feature of efficient CRC treatment. As miRNAs can be used as CRC biomarkers, the aim of the present study was to analyse experimentally validated data on frequently up-regulated miRNA clusters in CRC tissue and investigate their members with respect to clinicopathological characteristics of patients. Based on available data, 15 up-regulated clusters, miR-106a/363, miR-106b/93/25, miR-17/92a-1, miR-181a-1/181b-1, miR-181a-2/181b-2, miR-181c/181d, miR-183/96/182, miR-191/425, miR-200c/141, miR-203a/203b, miR-222/221, mir-23a/27a/24-2, mir-29b-1/29a, mir-301b/130b and mir-452/224, were selected. The positions of such clusters in the genome can be intronic or intergenic. Most clusters are regulated by several transcription factors, and miRNAs are also sponged by specific long non-coding RNAs. In some cases, co-expression of miRNA with other cluster members or host gene has been proven. miRNA expression patterns in cancer tissue, blood and faeces were compared. Based on experimental evidence, 181 target genes of selected clusters were identified. Panther analysis was used to reveal the functions of the target genes and their corresponding pathways. Clusters miR-17/92a-1, miR-106a/363, miR-106b/93/25 and miR-183/96/182 showed the strongest association with metastasis occurrence and poor patient survival, implicating them as the most promising targets of translational research.

Analysis of miRNA-mRNA Crosstalk in Radiation-Induced Mouse Thymic Lymphomas to Identify miR-486 as a Critical Regulator by Targeting IGF2BP3 mRNA

Ionizing radiation is one of the common environmental carcinogens. miRNAs play critical roles in the processes of tumor occurrence, development, metastasis. However, the relationship between radiation-induced carcinogenesis and miRNA rarely reported. This study is aimed to investigate the effect of miRNAs on radiation-induced carcinogenesis. In this study we established the radiation-induced thymic lymphoma mice model. By using miRNA array of RTL tissue and predicting for miRNAs target genes, a miRNA-mRNA crosstalk network was established. Based on this network, we identified a critical miRNA, miR-486, which was the most down-regulated in the radiation-induced carcinogenesis. Then the function of miR-486 was confirmed by using knockout mice and cellular experiments. As a result, miR-486 could inhibit proliferation of mouse lymphoma cells by targeting IGF2BP3 mRNA. The adenovirus over-expression miR-486 vector reduced tumorigenesis in vivo. MiR-486 knockout mice have a strong tendency of radiation-induced carcinogenesis. In conclusion, miR-486 inhibits the proliferation of lymphoma cells and tumorigenesis induced by radiation through targeting IGF2BP3.

MicroRNA-3148 acts as molecular switch promoting malignant transformation and adipocytic differentiation of immortalized human bone marrow stromal cells via direct targeting of the SMAD2/TGFβ pathway

MicroRNAs (miRs/miRNAs) play a key role in posttranscriptional regulation of gene expression and are implicated in a number of physiological and pathological conditions, including cellular malignant transformation. In the current study, we investigated the role of miR-3148 in regulating human stromal (mesenchymal) stem cell (hMSC) differentiation and transformation. Stable expression of miR-3148 in telomerized hMSC (hMSC-miR-3148) led to significant increase in in vitro adipocytic differentiation and suppression of osteoblastic differentiation. Concordantly, global gene expression profiling revealed significant enrichment in cholesterol biosynthesis pathway, and pathways related to enhanced cell movement and survival, whereas processes related to bone and connective tissue developments, cell death, apoptosis, and necrosis were downregulated. Global proteomic analysis using 2D-DIGE followed by mass spectrometry (MS) revealed significant changes in protein expression in hMSC-miR-3148 and enrichment in protein networks associated with carcinogenesis. Functional studies revealed that hMSC-miR-3148 exhibited enhanced in vitro cell proliferation, colony formation, migration, invasion, sphere formation, doxorubicin resistance, and increased active number of cells in S and G2/M cell cycle phases and formed sarcoma-like tumors with adipocyte infiltration when implanted into immunocompromised mice. SMAD2 was identified as bone fide gene target for miR-3148 using qRT-PCR, Western blotting, and UTR-based reporter assay. In agreement with our data, SMAD2 expression was downregulated in 47% of patients with soft tissue sarcoma. Bioinformatics analysis revealed that elevated miR-3148 expression correlates with poor prognosis in several human cancer types, including sarcoma. Our study identified miR-3148 as factor regulating hMSC differentiation and is involved in promoting malignant transformation of telomerized hMSC.

Bioinformatics Analysis Identifies the Estrogen Receptor 1 (ESR1) Gene and hsa-miR-26a-5p as Potential Prognostic Biomarkers in Patients with Intrahepatic Cholangiocarcinoma

Background

Intrahepatic cholangiocarcinoma arises from the epithelial cells of the bile ducts and is associated with poor prognosis. This study aimed to use bioinformatics analysis to identify molecular biomarkers of intrahepatic cholangiocarcinoma and their potential mechanisms.

Material/Methods

MicroRNA (miRNA) and mRNA microarrays from GSE53870 and GSE32879 were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed miRNAs (DEMs) associated with prognosis were identified using limma software and Kaplan-Meier survival analysis. Predictive target genes of the DEMs were identified using miRWalk, miRTarBase, miRDB, and TargetScan databases of miRNA-binding sites and targets. Target genes underwent Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Hub genes were analyzed by constructing the protein-protein interaction (PPI) network using Cytoscape. DEMs validated the hub genes, followed by construction of the miRNA-gene regulatory network.

Results

Twenty-five DEMs were identified. Fifteen DEMs were upregulated, and ten were down-regulated. Kaplan-Meier survival analysis identified seven upregulated DEMs and nine down-regulated DEMs that were associated with the overall survival (OS), and 130 target genes were selected. GO analysis showed that target genes were mainly enriched for metabolism and development processes. KEGG analysis showed that target genes were mainly enriched for cancer processes and some signaling pathways. Fourteen hub genes identified from the PPI network were associated with the regulation of cell proliferation. The overlap between hub genes and DEMs identified the estrogen receptor 1 (ESR1) gene and hsa-miR-26a-5p.

Conclusions

Bioinformatics analysis identified ESR1 and hsa-miR-26a-5p as potential prognostic biomarkers for intrahepatic cholangiocarcinoma.

MicroRNA Expression Profiling on Paired Primary and Lymph Node Metastatic Breast Cancer Revealed Distinct microRNA Profile Associated With LNM

Breast cancer (BC) is the foremost cause of cancer-related deaths in women. BC patients are oftentimes presented with lymph node metastasis (LNM), which increases their risk of recurrence. Compelling data have recently implicated microRNAs in promoting BC metastasis. Therefore, the identification of microRNA (miRNA)-based molecular signature associated with LNM could provide an opportunity for a more personalized treatment for BC patients with high risk of LNM. In current study, we performed comprehensive miRNA profiling in matched primary breast and LNM and identified 40 miRNAs, which were differentially expressed in LNM compared to primary tumors. The expression of 14 miRNAs (Up: hsa-miR-155-5p, hsa-miR-150-5p, hsa-miR-146a-5p, hsa-miR-142-5p and down: hsa-miR-200a-3p, hsa-miR-200b-3p, hsa-miR-200c-3p, hsa-miR-205-5p, hsa-miR-210-3p, hsa-miR-214-3p, hsa-miR-141-3p, hsa-miR-127-3p, hsa-miR-125a-5p, and hsa-let-7c-5p) was subsequently validated in a second cohort of 32 breast and 32 matched LNM tumor tissues. Mechanistically, forced expression of hsa-miR-205-5p, or hsa-miR-214-3p epigenetically inhibited MDA-MB-231 cell proliferation, colony formation, and cell migration. Global gene expression profiling on MDA-MB-231 cells overexpressing hsa-miR-205-5p, or hsa-miR-214-3p in combination with in silico target prediction and ingenuity pathway analyses identified multiple bona fide targets for hsa-miR-205-5p, hsa-miR-214-3p affecting cellular proliferation and migration. Interestingly, interrogation of the expression levels of hsa-miR-205 and hsa-miR-214 in the METABRIC breast cancer dataset revealed significantly poor overall survival in patients with downregulated expression of miR-205 [HR = 0.75 (0.61–0.91)], p = 0.003 and hsa-miR-214 [HR = 0.74 (0.59–0.93) p = 0.008]. Our data unraveled the miRNA-transcriptional landscape associated with LNM and provide novel insight on the role of several miRNAs in promoting BC LNM, and suggest their potential utilization in the clinical management of BC patients.

Transgelin is a poor prognostic factor associated with advanced colorectal cancer (CRC) stage promoting tumor growth and migration in a TGFβ-dependent manner

Colorectal cancer (CRC) is the fourth most common cancer type globally. Investigating the signaling pathways that maintain cancer cell phenotype can identify new biomarkers for targeted therapy. Aberrant transforming growth factor-β (TGFβ) signaling has been implicated in CRC progression, however, the exact mechanism by which TGFβ exerts its function is still being unraveled. Herein, we investigated TAGLN expression, prognostic value, and its regulation by TGFβ in CRC. While TAGLN was generally found to be downregulated in CRC, elevated expression of TAGLN was associated with advanced CRC stage and predicted poor overall survival (hazard ratio (HR) = 1.8, log-rank test P-value = 0.014) and disease-free survival (HR = 1.6, log-rank test P-value = 0.046), hence implicating TAGLN as poor prognostic factor in CRC. Forced expression of TAGLN was associated with enhanced CRC cell proliferation, clonogenic growth, cell migration and in vivo tumor formation in immunocompromised mice, while targeted depletion of TAGLN exhibited opposing biological effects. Global gene expression profiling of TAGLN-overexpressing or TAGLN-deficient CRC cell lines revealed deregulation of multiple cancer-related genes and signaling pathways. Transmission electron microscopy (TEM) revealed ultrastructural changes due to loss of TAGLN, including disruption of actin cytoskeleton organization and aberrant actin filament distribution. Hierarchical clustering, principle component, and ingenuity pathway analyses revealed distinct molecular profile associated with TAGLN^high CRC patients with remarkable activation of a number of mechanistic networks, including SMARCA4, TGFβ1, and P38 MAPK. The P38 MAPK was the top predicted upstream regulator network promoting cell movement through regulation of several intermediate molecules, including TGFβ1. Concordantly, functional categories associated with cellular movement and angiogenesis were also enriched in TAGLN^high CRC, supporting a model for the molecular mechanisms linking TGFβ-induced upregulation of TAGLN and CRC tumor progression and suggesting TAGLN as potential prognostic marker associated with advanced CRC pathological stage.

Comparative bioinformatical analysis of pancreatic head cancer and pancreatic body/tail cancer

This study is to analyze differentially expressed genes (DEGs) and mutation signatures of pancreatic head cancer and pancreatic body/tail cancer. Pancreatic Adenocarcinoma (PAAD) RNA-seq data, mutation data and clinical data were downloaded and collected from The Cancer Genome Atlas (TCGA), FireHose and CBioPortal. According to the anatomic location, the patients were divided into 146 cases of pancreatic head cancer and 28 cases of pancreatic body/tail cancer. Then survival analysis was performed by Kaplan-Meier and log-rank test. Furthermore, DEGs were screened by R package Deseq2. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein-protein interaction (PPI) were then carried out by DAVID and String. Online tool TIMER was used to analyze the immune cells infiltration. R package maftools and GenVisR were applied to analyze frequently mutated genes and mutant-allele tumor heterogeneity (MATH) of PAAD. Survival of patients with pancreatic body/tail cancer was better than those with pancreatic head cancer (median survival, 24.05 vs 19.45 months, p = 0.048). And 496 significant DEGs (|log2 FoldChange| > 1.5,false discovery rate (FDR) < 0.05) were identified, including 253 downregulated genes and 243 upregulated genes. And there were 13 Go terms (4 biological processes, 6 cellular components and 3 molecular functions) and 3 KEGG pathways (Pancreatic secretion, Fat digestion and absorption, Protein digestion and absorption) (FDR < 0.05). B cells and CD4 + T cells infiltration were more significant in pancreatic head cancer. MATH scores of pancreatic body/tail cancer were higher than pancreatic head cancer, while χ² test of top 10 frequently mutated genes showed little difference between them. There were prognostic and genetic differences between pancreatic head cancer and pancreatic body/tail cancer. PAAD originated from different location may have different biology natures and should not be treated with same strategy.

Effect of miR-181a-3p on osteogenic differentiation of human bone marrow-derived mesenchymal stem cells by targeting BMP10

Objective: To explore the regulation relationship between miR-181a-3p and BMP10, and their mechanism of osteogenic differentiation of human bone marrow-derived mesenchymal stem cells (MSCs).Methods: After osteogenic induction of MSCs, the ALP activity was detected by ELISA. The expression of miRNA-181a-3p and BMP10 was detected by RT-qPCR, and the protein levels of BMP10 and osteogenic differentiation marker proteins ALK and RUNX2 were detected by Western blot. The TargetScan online website was used to predict the putative target of miR-181a-3p, and dual luciferase reporter assay was performed to validate the targeting relationship between miR-181a-3p and BMP10.Results: In osteogenic differentiation of MSCs, ALP activity, the level of ALK and RUNX2 was evidently increased (p < .05), and the expression of miR-181a-3p was significantly downregulated (p < .05). Moreover, overexpression of miR-181a-3p obviously decreased the expression of BMP10 (p < .05), miR-181a-3p knockdown increased the expression of BMP10 prominently (p < .05). The transfection of miR-181a-3p mimics resulted in significantly downregulation of ALP activity and RUNX2 protein expression in MSCs (p < .05). In addition, overexpression of BMP10 could reverse the inhibitory effect of miR-181a-3p on osteogenic differentiation (p < .05).Conclusions: In conclusion, we found that miR-181a-3p inhibited osteogenic differentiation of MCSs by targeting BMP10.

MicroRNA Expression Profiling on Paired Primary and Lymph Node Metastatic Breast Cancer Revealed Distinct microRNA Profile Associated With LNM

Breast cancer (BC) is the foremost cause of cancer-related deaths in women. BC patients are oftentimes presented with lymph node metastasis (LNM), which increases their risk of recurrence. Compelling data have recently implicated microRNAs in promoting BC metastasis. Therefore, the identification of microRNA (miRNA)-based molecular signature associated with LNM could provide an opportunity for a more personalized treatment for BC patients with high risk of LNM. In current study, we performed comprehensive miRNA profiling in matched primary breast and LNM and identified 40 miRNAs, which were differentially expressed in LNM compared to primary tumors. The expression of 14 miRNAs (Up: hsa-miR-155-5p, hsa-miR-150-5p, hsa-miR-146a-5p, hsa-miR-142-5p and down: hsa-miR-200a-3p, hsa-miR-200b-3p, hsa-miR-200c-3p, hsa-miR-205-5p, hsa-miR-210-3p, hsa-miR-214-3p, hsa-miR-141-3p, hsa-miR-127-3p, hsa-miR-125a-5p, and hsa-let-7c-5p) was subsequently validated in a second cohort of 32 breast and 32 matched LNM tumor tissues. Mechanistically, forced expression of hsa-miR-205-5p, or hsa-miR-214-3p epigenetically inhibited MDA-MB-231 cell proliferation, colony formation, and cell migration. Global gene expression profiling on MDA-MB-231 cells overexpressing hsa-miR-205-5p, or hsa-miR-214-3p in combination with in silico target prediction and ingenuity pathway analyses identified multiple bona fide targets for hsa-miR-205-5p, hsa-miR-214-3p affecting cellular proliferation and migration. Interestingly, interrogation of the expression levels of hsa-miR-205 and hsa-miR-214 in the METABRIC breast cancer dataset revealed significantly poor overall survival in patients with downregulated expression of miR-205 [HR = 0.75 (0.61-0.91)], p = 0.003 and hsa-miR-214 [HR = 0.74 (0.59-0.93) p = 0.008]. Our data unraveled the miRNA-transcriptional landscape associated with LNM and provide novel insight on the role of several miRNAs in promoting BC LNM, and suggest their potential utilization in the clinical management of BC patients.

The Role of TGF-β Signaling Regulatory MicroRNAs in the Pathogenesis of Colorectal Cancer

Colorectal cancer (CRC) is one of the most common cancers globally and is associated with a high mortality rate. The transforming growth factor beta (TGF-β) signaling pathway plays an important role in normal intestinal tissue function, but has also been implicated in the development of CRC. MicroRNAs (miRNAs) have also recently emerged as important regulators of cancer development and progression. They act by targeting multiple signaling pathways including the TGF-β signaling pathway. There is growing evidence demonstrating that miRNAs target various components of the TGF-β signaling pathway, including TGF-β1, TGF-β2, regulatory SMADs (SMAD1, 2, 3, 5 and 9), co-mediator SMAD4, inhibitory SMADs (SMAD6 and 7) and the TGF-β receptors, and thereby alter the proliferation and migration of CRC cells. In this review, we summarize the data concerning the interaction between TGF-β signaling pathway and miRNAs with the aim to better understanding the CRC molecular mechanisms and hence better management of this disease.

Integrated Transcriptome and Pathway Analyses Revealed Multiple Activated Pathways in Breast Cancer

Breast cancer (BC) is the leading cause of cancer-related death in women. Therefore, a better understanding of BC biology and signaling pathways might lead to the development of novel biomarkers and targeted therapies. Although a number of transcriptomic studies have been performed on breast cancer patients from various geographic regions, there are almost no such comprehensive studies performed on breast cancer from patients in the gulf region. This study aimed to provide a better understanding of the altered molecular networks in BC from the gulf region. Herein, we compared the transcriptome of BC to adjacent normal tissue from six BC patients and identified 1,108 upregulated and 518 downregulated transcripts. A selected number of genes from the RNA-Seq analysis were subsequently validated using qRT-PCR. Differentially expressed (2.0-fold change, adj. p < 0.05) transcripts were subjected to ingenuity pathway analysis, which revealed a myriad of affected signaling pathways and functional categories. Activation of ERBB2, FOXM1, ESR1, and IGFBP2 mechanistic networks was most prominent in BC tissue. Additionally, BC tissue exhibited marked enrichment in genes promoting cellular proliferation, migration, survival, and DNA replication and repair. The presence of genes indicative of immune cell infiltration and activation was also observed in BC tissue. We observed high concordance [43.5% (upregulated) and 62.1% (downregulated)] between differentially expressed genes in our study group and those reported for the TCGA BC cohort. Our data provide novel insight on BC biology and suggest common altered molecular networks in BC in this geographic region. Our data suggest future development of therapeutic interventions targeting those common signaling pathways.

A functional variant in the flanking region of pri-let-7f contributes to colorectal cancer risk in a Chinese population

Let-7f was reported to be downregulated in patients with colorectal cancer (CRC). However, little is known about the role of let-7f in CRC carcinogenesis. The aim of this study was to investigate the correlation between genetic polymorphisms in the flanking region of pri-let-7f and CRC risk, as well as the potential role of let-7f in CRC cell migration and invasion. The pri-let-7f-1 rs10739971 and pri-let-7f-2 rs17276588 were genotyped using TaqMan (Applied Biosystems, Foster City, CA) assay. The luciferase activity was detected using Dual-Luciferase Reporter Assay. CRC cell migration and invasion were evaluated using transwell chamber assay. The rs17276588 AG and AG/AA genotypes had a significantly increased CRC risk (AG vs. GG: adjusted odds ratio [OR] = 1.48, 95% confidence interval [CI] = 1.19-1.83, p < 0.001; AG/AA vs. GG: adjusted OR = 1.43, 95% CI = 1.17-1.75, p < 0.001). Stratification analyses showed that the increased risk was observed in CRC patients with well-moderately differential status, patients with clinical Stages I-II, and patients without lymph node metastasis. The rs17276588A allele displayed a decreased transcriptional activity and low levels of let-7f. Moreover, let-7f inhibited migration and invasion in Caco-2 and Lovo cells. These findings indicate that the rs17276588 AG/AA genotypes increased CRC risk by reducing the expression of tumor suppressor let-7f.

MicroRNAs as potential markers of parenteral nutrition associated liver disease in adult patients

Parenteral nutrition-associated liver disease (PNALD) is a severe complication in patients completely dependent on parenteral nutrition (PN). The gold diagnostic standard, liver biopsy, is associated with significant health risk and therefore its use is limited. MicroRNAs (miRNAs) are small non-coding regulatory RNA molecules with highly tissue-specific expression and the secreted miRNAs may serve as non-invasive diagnostic biomarkers. The aim of this study was to evaluate the expression of a panel of specific miRNAs associated with liver diseases of different origin in PN-dependent adult patients in order to design miRNA panel enabling to precise monitoring of PNALD progression. Twelve PN-dependent patients with short bowel syndrome (SBS) were monitored on three/four-month basis for up to 24 months. Forty-five age- and sex-matched subjects without any known liver pathology served as controls. Specific miRNAs expression was determined by RT-qPCR using TaqMan probes (Thermofisher). Liver function test parameters were determined in certified clinical laboratories. Six of the tested miRNAs exhibited significantly altered expression compared with healthy controls, three of them (MIR122, MIR1273g, and MIR500a) were upregulated while three were down-regulated (MIR505, MIR199a, MIR139). MIR122 positively correlated with serum AST and ALT activities while MIR1273g positively correlated with serum CRP concentration and GGT activity. MIR505, MIR199a, and MIR139 negatively correlated with serum GGT activity. Fluctuation of these parameters well paralleled serum miRNA concentrations in all patients throughout the whole observation period. We identified six miRNAs whose serum concentrations are significantly altered in PN-dependent patients with PNALD and correlate with markers of inflammation, cholestasis or hepatic injury. Their reliability as markers of PNALD progression needs to be further evaluated.

m1A Regulated Genes Modulate PI3K/AKT/mTOR and ErbB Pathways in Gastrointestinal Cancer

BACKGROUND: Gene expression can be posttranscriptionally regulated by a complex network of proteins. N1-methyladenosine (m1A) is a newly validated RNA modification. However, little is known about both its influence and biogenesis in tumor development. METHODS: This study analyzed TCGA data of patients with five kinds of gastrointestinal (GI) cancers. Using data from cBioPortal, molecular features of the nine known m1A-related enzymes in GI cancers were investigated. Using a variety of bioinformatics approach, the impact of m1A regulators on its downstream signaling pathway was studied. To further confirm this regulation, the effect of m1A writer ALKBH3 knockdown was studied using RNA-seq data from published database. RESULTS: Dysregulation and multiple types of genetic alteration of putative m1A-related enzymes in tumor samples were observed. The ErbB and mTOR pathways with ErbB2, mTOR, and AKT1S1 hub genes were identified as being regulated by m1A-related enzymes. The expression of both ErbB2 and AKT1S1 was decreased after m1A writer ALKBH3 knockdown. Furthermore, Gene Ontology analysis revealed that m1A downstream genes were associated with cell proliferation, and the results showed that m1A genes are reliably linked to mTOR. CONCLUSION: This study demonstrated for the first time the dysregulation of m1A regulators in GI cancer and its signaling pathways and will contribute to the understanding of RNA modification in cancer.

SOX4: Epigenetic regulation and role in tumorigenesis

Sex-determining region Y-related (SRY) high-mobility group box 4 (SOX4) is a member of the group C subfamily of SOX transcription factors and promotes tumorigenesis by endowing cancer cells with survival, migratory, and invasive capacities. Emerging evidence has highlighted an unequivocal role for this transcription factor in mediating various signaling pathways involved in tumorigenesis, epithelial-to-mesenchymal transition (EMT), and tumor progression. During the last decade, numerous studies have highlighted the epigenetic interplay between SOX4-targeting microRNAs (miRNAs), long noncoding RNAs (lncRNAs) and SOX4 and the subsequent modulation of tumorigenesis, invasion and metastasis. In this review, we summarize the current state of knowledge about the role of SOX4 in cancer development and progression, the epigenetic regulation of SOX4, and the potential utilization of SOX4 as a diagnostic and prognostic biomarker and its depletion as a therapeutic target.

Circulating miR-141-3p, miR-143-3p and miR-200c-3p are differentially expressed in colorectal cancer and advanced adenomas

Colorectal cancer (CRC) is one of the prominent causes of cancer related deaths because, in part, there is not an early, non-invasive, effective detection strategy. Circulating microRNAs (miRNAs) have been proposed as potential non-invasive biomarkers for CRC. In this study, we evaluated the miRNA profile in sixteen CRC tissues by Next-Generation-Sequencing and compared the circulating expression levels of 22 miRNAs among 45 CRC, 14 hyperplastic polyps, 11 advanced adenoma patients and 45 control subjects, by reverse transcription-quantitative PCR, to search for miRNAs which could be potential biomarkers. In total, nine of them represented 70% of total read counts (miR-10a-5p, miR-192-5p, miR-10b-5p, miR-22-3p, miR-26a-5p, miR-148a-3p, miR-181a-5p, miR-92a-3p and miR-143-5p). In silico analysis found eight candidates to mature miRNAs. With respect to circulating miRNA, we found higher serum expression levels of miR-143-3p, miR-141-3p and miR-200c-3p in the CRC and adenoma groups compared with controls (P<0.002), and we also found significant higher levels of miR-141-3p and miR-200c-3p in serum of adenoma patients compared with the CRC group. In conclusion, the measurement of miRNAs in the blood could complement current screening methods for CRC and might provide new insights into mechanisms of tumorigenesis. miR-143-3p, miR-141-3p and miR-200c-3p could be interesting miRNAs to study as potential biomarkers for CRC.

Neoplastic Transformation of Human Mesenchymal Stromal Cells Mediated via LIN28B

Bone marrow stromal (Mesenchymal) stem cells (MSCs) are multipotent bone cells capable of differentiating into mesoderm-type cells, such as osteoblasts and adipocytes. Existing evidence suggests that transformation of MSCs gives rise to sarcoma. In order to identify the molecular mechanism leading to spontaneous transformation of human bone marrow MSCs (hBMSCs), we performed comprehensive microRNA (miRNA) and mRNA profiling in the transformed hBMSC-Tum line compared to the parental clone. As a result, we identified multiple dysregulated molecular networks associated with the hBMSC transformed phenotype. LIN28B was upregulated 177.0-fold in hBMSC-Tum, which was associated with marked reduction in LET-7 expression and upregulated expression of its target HMGA2. Targeted depletion of LIN28B or exogenous expression of LET-7b suppressed hBMSC-Tum proliferation, colony formation, and migration. On the other hand, forced expression of LIN28B promoted malignant transformation of parental hBMSC cells as shown by enhanced in vitro colony formation, doxorubicin resistance, and in vivo tumor formation in immunocompromised mice. Analysis of LIN28B and HMGA2 expression levels in cohorts from The Cancer Genome Atlas sarcoma dataset revealed a strong inverse-relationship between elevated expression and overall survival (OS) in 260 patients (p = 0.005) and disease-free survival (DFS) in 231 patients (p = 0.02), suggesting LIN28B and HMGA2 are important regulators of sarcoma biology. Our results highlight an important role for the LIN28B/LET-7 axis in human sarcoma pathogenesis and suggest that the therapeutic targeting of LIN28B may be relevant for patients with sarcoma.

DNA methylation and chromatin modifiers in colorectal cancer

Colorectal carcinogenesis is a multistep process involving the accumulation of genetic alterations over time that ultimately leads to disease progression and metastasis. Binding of transcription factors to gene promoter regions alone cannot explain the complex regulation pattern of gene expression during this process. It is the chromatin structure that allows for a high grade of regulatory flexibility for gene expression. Posttranslational modifications on histone proteins such as acetylation, methylation, or phosphorylation determine the accessibility of transcription factors to DNA. DNA methylation, a chemical modification of DNA that modulates chromatin structure and gene transcription acts in concert with these chromatin conformation alterations. Another epigenetic mechanism regulating gene expression is represented by small non-coding RNAs. Only very recently epigenetic alterations have been included in molecular subtype classification of colorectal cancer (CRC). In this chapter, we will provide examples of the different epigenetic players, focus on their role for epithelial-mesenchymal transition and metastatic processes and discuss their prognostic value in CRC.

Joint bioinformatics analysis of underlying potential functions of hsa-let-7b-5p and core genes in human glioma

Background

Glioma accounts for a large proportion of cancer, and an effective treatment for this disease is still lacking because of the absence of specific driver molecules. Current challenges in the treatment of glioma are the accurate and timely diagnosis of brain glioma and targeted treatment plans. To investigate the diagnostic biomarkers and prospective role of miRNAs in the tumorigenesis and progression of glioma, we analyzed the expression of miRNAs and key genes in glioma based on The Cancer Genome Atlas database.

Methods

Of the 701 cases that were downloaded, five were normal and 696 were glioma. Then, 1626 differentially expressed genes were identified, and 173 aberrantly expressed miRNAs were calculated by edgeR. GO and KEGG pathway enrichment analyses were performed using Cytoscape software. A coexpression network was built by weighted correlation network analysis (WGCNA). A cell scratch test and transwell, cell apoptosis and cell cycle assays were performed to validate the function of hsa-let-7b-5p.

Results

Based on crosstalk genes in the KEGG, PPI network, and WGCNA analyses, PLK1, CCNA2, cyclin B2 (CCNB2), and AURKA were screened as candidate diagnostic marker genes. The survival analysis revealed that high mRNA expression of PLK1, CCNA2, and AURKA was significantly associated with poor overall survival. Furthermore, hsa-let-7b-5p was identified as a core miRNA in the regulation of candidate genes involved in glioma development. We confirmed that hsa-let-7b-5p could inhibit the migration, invasion, and cell cycle of glioma cells.

Conclusions

This study provides four potential biomarkers for the diagnosis of glioma, offers a potential explanation of its pathogenesis, and proposes hsa-let-7b-5p as a therapeutic target.

Electronic supplementary material

The online version of this article (10.1186/s12967-019-1882-7) contains supplementary material, which is available to authorized users.

CXCR7 signaling promotes breast cancer survival in response to mesenchymal stromal stem cell-derived factors

The interaction between cancer cells and molecular cues provided by tumor stromal cells plays a crucial role in cancer growth and progression. We have recently reported that the outcome of interaction between tumor cells and stromal cells is dependent on the gene expression signature of tumor cells. In the current study, we observed that several cancer cell lines, e.g., MCF7 breast cancer line, exhibited growth advantage when cultured in the presence of conditioned media (CM) derived from human bone marrow stromal stem cells (hBMSCs). Regarding the underlying molecular mechanism, we have identified CXCR7 as highly expressed by MCF7 cells and that it mediated the enhanced growth in response to hBMSC CM. Regarding the clinical relevance, we found an inverse correlation between the level of tumor gene expression of CXCR7 in bladder, breast, cervical, kidney, liver, lung, pancreatic, stomach, and uterine cancers, and patients' overall survival. Interestingly, significant positive correlation between CXCR7 and CXCL12 gene expression (Pearson = 0.3, p = 2.0 × 10^-16) was observed in breast cancer patients, suggesting a biological role for the CXCR7/CXCL12 genetic circuit in breast cancer biology. Our data provide insight into the molecular mechanisms by which stromal-derived microenvironmental cues mediate CXCR7 signaling and growth enhancement of breast cancer cells. Therapeutic targeting of this circuit might provide novel therapeutic opportunity for breast cancer.

Differential expression profiles and functional analysis of plasma miRNAs associated with chronic myeloid leukemia phases

Aim: This study was aimed to investigate the expression profiles and biological function of plasma miRNAs at different phases of chronic myeloid leukemia (CML). Materials & methods: Differentially expressed miRNAs were identified by microarray. The candidate miRNAs were validated by quantitative real-time PCR at chronic phase, accelerated phase and blast crisis. The functional analysis of miRNAs was carried out by using DAVID. Results: The putative targets of dysregulated miRNAs were involved in important signaling pathways. Plasma let-7b-5p and miR-451a expression was lower in CML patients, and plasma miR-451a gradually decreased from chronic phase to accelerated phase and blast crisis. Conclusion: Dysregulated plasma miRNAs maybe play regulatory roles in pathogenesis of CML. Let-7b-5p and miR-451a can be used as potential biomarkers for the diagnosis and prognosis of CML.

Emerging microRNA biomarkers for colorectal cancer diagnosis and prognosis

MicroRNAs (miRNAs) are one abundant class of small, endogenous non-coding RNAs, which regulate various biological processes by inhibiting expression of target genes. miRNAs have important functional roles in carcinogenesis and development of colorectal cancer (CRC), and emerging evidence has indicated the feasibility of miRNAs as robust cancer biomarkers. This review summarizes the progress in miRNA-related research, including study of its oncogene or tumour-suppressor roles and the advantages of miRNA biomarkers for CRC diagnosis, treatment and recurrence prediction. Along with analytical technique improvements in miRNA research, use of the emerging extracellular miRNAs is feasible for CRC diagnosis and prognosis.

Computational Analysis of miRNA and their Gene Targets Significantly Involved in Colorectal Cancer Progression

BACKGROUND

Globally, colorectal cancer (CRC) is the third most common cancer in women and the fourth most common cancer in men. Dysregulation of small non-coding miRNAs have been correlated with colon cancer progression. Since there are increasing reports of candidate miRNAs as potential biomarkers for CRC, this makes it important to explore common miRNA biomarkers for colon cancer. As computational prediction of miRNA targets is a critical initial step in identifying miRNA: mRNA target interactions for validation, we aim here to construct a potential miRNA network and its gene targets for colon cancer from previously reported candidate miRNAs, inclusive of 10 up- and 9 down-regulated miRNAs from tissues; and 10 circulatory miRNAs.

METHODS

The gene targets were predicted using DIANA-microT-CDS and TarBaseV7.0 databases. Each miRNA and its targets were analyzed further for colon cancer hotspot genes, whereupon DAVID analysis and mirPath were used for KEGG pathway analysis.

RESULTS

We have predicted 874 and 157 gene targets for tissue and serum specific miRNA candidates, respectively. The enrichment of miRNA revealed that particularly hsa-miR-424-5p, hsa-miR-96-5p, hsa-miR-1290, hsa-miR-224, hsa-miR-133a and has-miR-363-3p present possible targets for colon cancer hallmark genes, including BRAF, KRAS, EGFR, APC, amongst others. DAVID analysis of miRNA and associated gene targets revealed the KEGG pathways most related to cancer and colon cancer. Similar results were observed in mirPath analysis. A new insight gained in the colon cancer network pathway was the association of hsa-mir-133a and hsa-mir-96-5p with the PI3K-AKT signaling pathway. In the present study, target prediction shows that while hsa-mir-424-5p has an association with mostly 10 colon cancer hallmark genes, only their associations with MAP2 and CCND1 have been experimentally validated.

CONCLUSION

These miRNAs and their targets require further evaluation for a better understanding of their associations, ultimately with the potential to develop novel therapeutic targets.

MicroRNA-4458 suppresses migration and epithelial-mesenchymal transition via targeting HMGA1 in non-small-cell lung cancer cells

Purpose

Increasing studies have shown that microRNA-4458 (miR-4458) is associated with human cancer progression. However, the molecular mechanism of miR-4458 in non-small-cell lung cancer (NSCLC) remains largely unknown. This study aims to reveal the biological function of miR-4458 in NSCLC.

Materials and methods

The expression of miR-4458 in NSCLC cells was evaluated by qRT-PCR. Cell proliferation and migration assay were carried out in vitro after transfection. A luciferase reporter and Western blot assay were performed to identify the functional target of miR-4458.

Results

The study indicated that miR-4458 was markedly downregulated in NSCLC cells. Overexpression of miR-4458 strongly reduced the proliferation and migration in NSCLC cell lines. In addition, miR-4458 inhibited the progression of migration and epithelial–mesenchymal transition (EMT) through the PI3K/AKT pathway. Luciferase report assay demonstrated that HMGA1 was a target gene for miR-4458.

Conclusion

The results indicate that miR-4458 participated in the process of migration and EMT via directly targeting HMGA1 and miR-4458 might be a potential novel therapeutic target in NSCLC.

Identification of a 6-gene signature predicting prognosis for colorectal cancer

BACKGROUND

An accurate and robust gene signature is of the utmost importance in assisting oncologists to make a more accurate evaluation in clinical practice. In our study, we extracted key mRNAs significantly related to colorectal cancer (CRC) prognosis and we constructed an expression-based gene signature to predict CRC patients' survival.

METHODS

mRNA expression profiles and clinicopathological data of colon adenocarcinoma (COAD) cases and rectum adenocarcinoma (READ) were collected from The Cancer Genome Atlas database to investigate gene expression alteration associated to the prognosis of CRC. Differentially expressed mRNAs (DEMs) were detected between COAD/READ and normal tissue samples. Relying on a univariate and multivariate Cox regression analyses, a mRNA panel signature was established and used for predicting the overall survival (OS) in CRC patients. Receiver operating characteristic curve was used to evaluate the prognosis performance of our model through calculating the AUC values corresponding to the 3-year and 5-year survival. To assess the performance of gene signature in the given cancer subgroups (CRC entire cohort, COAD cohort, and READ cohort), a stratified analysis was carried out according to clinical factors.

RESULTS

A total of 5341 and 5594 DEMs were collected from COAD vs. normal tissue samples, and READ vs. normal samples respectively. A univariate regression analysis for the common DEMs between COAD and READ cohorts resulted in 14 common mRNAs related to OS. The multivariate Cox regression analysis revealed that 6 of these mRNAs (EPHA6, TIMP1, IRX6, ART5, HIST3H2BB, and FOXD1) had significant prognostic value allowing the discrimination between high- and low-risk patients, implying poor and good outcomes, respectively. The stratified analysis identified 6-gene signature as an independent prognostic signature in predicting CRC patients' survival.

CONCLUSIONS

The 6-gene signature could act as an independent biomarker for survival prediction of CRC patients.

An Integrated Approach for Identification of Functionally Similar MicroRNAs in Colorectal Cancer

Colorectal cancer (CRC) is one of the most prevalent cancers around the globe. However, the molecular reasons for pathogenesis of CRC are still poorly understood. Recently, the role of microRNAs or miRNAs in the initiation and progression of CRC has been studied. MicroRNAs are small, endogenous noncoding RNAs found in plants, animals, and some viruses, which function in RNA silencing and posttranscriptional regulation of gene expression. Their role in CRC development is studied and they are found to be potential biomarkers in diagnosis and treatment of CRC. Therefore, identification of functionally similar CRC related miRNAs may help in the development of a prognostic tool. In this regard, this paper presents a new algorithm, called μSim. It is an integrative approach for identification of functionally similar miRNAs associated with CRC. It integrates judiciously the information of miRNA expression data and miRNA-miRNA functionally synergistic network data. The functional similarity is calculated based on both miRNA expression data and miRNA-miRNA functionally synergistic network data. The effectiveness of the proposed method in comparison to other related methods is shown on four CRC miRNA data sets. The proposed method selected more significant miRNAs related to CRC as compared to other related methods.

Whole genome mRNA expression profiling revealed multiple deregulated pathways in stromal vascular fraction from erectile dysfunction patients

Background: Stem-cell-based therapies have recently been explored in the field of erectile dysfunction (ED). However, the cellular and molecular phenotype of adipose derived stem cells (ADSCs) stromal vascular fraction (SVF) from ED patients remains largely unknown. Herein we compared the global gene expression profile in the SVF from ED patients and healthy individuals and identified altered signaling pathways between the two groups.

Methods: Samples (2–5 g) of abdominal adipose tissue from ED patients (n = 6) and healthy individual controls (n = 3) undergoing elective cosmetic liposuction were collected. Immediately after removal, SVF was separated using Collagenase type I and type IV protocol. RNA was isolated and microarray experiments were conducted using the Agilent platform. Data were normalized and pathway analyses were performed using GeneSpring software.

Results: Our data revealed multiple differentially expressed genes between the ED and control group. Hierarchical clustering based on differentially expressed mRNAs revealed clear separation of the two groups. The distribution of the top enriched pathways for the up-regulated genes indicated enrichment in inflammatory response and T-cell receptor signaling, while pathway analysis performed on the down-regulated genes revealed enrichment in mitogen-activated protein kinase, TGF-β, senescence, FAK, adipogenesis, androgen receptor, and EGF–EGFR signaling pathways in SVF from ED patient.

Conclusion: Our data revealed the existence of multiple altered signaling pathways in the SVF from ED patients, which could potentially play a role in the etiology of this disease. Therefore, therapeutic strategies targeting these pathways might provide novel therapeutic opportunity for ED patients.

Independent prognostic genes and mechanism investigation for colon cancer

PROPOSE

We aimed to explore the potential molecular mechanism and independent prognostic genes for colon cancer (CC).

METHODS

Microarray datasets GSE17536 and GSE39582 were downloaded from Gene Expression Omnibus. Meanwhile, the whole CC-related dataset were downloaded from The Cancer Genome Atlas (TCGA) database. Differentially expressed mRNA (DEMs) were identified between cancer tissue samples and para-carcinoma tissue samples in TCGA dataset, followed by the KEGG pathway and GO function analyses. Furthermore, the clinical prognostic analysis including overall survival (OS) and disease-free survival (DFS) were performed in all three datasets.

RESULTS

A total of 633 up- and 321 down-regulated mRNAs were revealed in TCGA dataset. The up-regulated mRNAs were mainly assembled in functions including extracellular matrix and pathways including Wnt signaling. The down-regulated mRNAs were mainly assembled in functions like Digestion and pathways like Drug metabolism. Furthermore, up-regulation of UL16-binding protein 2 (ULBP2) was associated with OS in CC patients. A total of 12 DEMs including Surfactant Associated 2 (SFTA2) were potential DFS prognostic genes in CC patients. Meanwhile, the GRP and Transmembrane Protein 37 (TMEM37) were two outstanding independent DFS prognostic genes in CC.

CONCLUSIONS

ULBP2 might be a potential novel OS prognostic biomarker in CC, while GRP and TMEM37 could be served as the independent DFS prognostic genes in CC. Furthermore, functions including extracellular matrix and digestion, as well as pathways including Wnt signaling and drug metabolism might play important roles in the process of CC.

miR-223-3p Inhibits Human Osteosarcoma Metastasis and Progression by Directly Targeting CDH6

Cadherin-6 (CDH6) is aberrantly expressed in cancer and closely associated with tumor progression. However, the functions of CDH6 in human osteosarcoma and the molecular mechanisms underlying CDH6 in osteosarcoma oncogenesis remain poorly understood. In this work, we assessed the role of CDH6 in human osteosarcoma and identified that the expression of CDH6 was closely related with the overall survival and poor prognosis of osteosarcoma patients. MicroRNAs (miRNAs) have been implicated as important epigenetic regulators during the progression of osteosarcoma. Using dual-luciferase reporter assays, we showed that miR-223-3p suppresses CDH6 expression by directly binding to the 3' UTR of CDH6. miR-223-3p overexpression significantly inhibited cell invasion, migration, growth, and proliferation by suppressing the CDH6 expression in vivo and in vitro. Besides, CDH6 overexpression in the miR-223-3p-transfected osteosarcoma cells effectively rescued the inhibition of cell invasion, migration, growth, and proliferation mediated by miR-223-3p. Additionally, Kaplan-Meier analysis suggests that the expression of miR-223-3p predicts favorable clinical outcomes for osteosarcoma patients. Moreover, the expression of miR-223-3p was downregulated in osteosarcoma patients and was negatively associated with the expression of CDH6. Collectively, these data highlight that miR-223-3p/CDH6 axis is an important novel pleiotropic regulator and could early predict the metastatic potential in human osteosarcoma treatments.

The Application of Gene Expression Profiling in Predictions of Occult Lymph Node Metastasis in Colorectal Cancer Patients

A key factor in determining the likely outcome for a patient with colorectal cancer is whether or not the tumour has metastasised to the lymph nodes-information which is also important in assessing any possibilities of lymph node resection so as to improve survival. In this review we perform a wide-range assessment of literature relating to recent developments in gene expression profiling (GEP) of the primary tumour, to determine their utility in assessing node status. A set of characteristic genes seems to be involved in the prediction of lymph node metastasis (LNM) in colorectal patients. Hence, GEP is applicable in personalised/individualised/tailored therapies and provides insights into developing novel therapeutic targets. Not only is GEP useful in prediction of LNM, but it also allows classification based on differences such as sample size, target gene expression, and examination method.

Molecular profiling of ALDH1+ colorectal cancer stem cells reveals preferential activation of MAPK, FAK, and oxidative stress pro-survival signalling pathways

Tumour heterogeneity leads to variable clinical response and inaccurate diagnostic and prognostic assessment. Cancer stem cells (CSCs) represent a subpopulation responsible for invasion, metastasis, therapeutic resistance, and recurrence in many human cancer types. However, the true identity of colorectal cancer (CRC) SCs remains elusive. Here, we aimed to characterize and define the gene expression portrait of CSCs in CRC-model SW403 cells. We found that ALDH⁺ positive cells are clonogenic and highly proliferative; their global gene expression profiling-based molecular signature revealed gene enrichment related to DNA damage, MAPK, FAK, oxidative stress response, and Wnt signalling. ALDH⁺ cells showed enhanced ROS stress resistance, whereas MAPK/FAK pathway pharmacologic inhibition limited their survival. Conversely, 5-fluorouracil increased the ALDH⁺ cell fraction among the SW403, HCT116 and SW620 CRC models. Notably, analysis of ALDH1A1 and POU5F1 expression levels in cohorts of 462 or 420 patients for overall (OS) or disease-free (DFS) survival, respectively, obtained from the Cancer Genome Atlas CRC dataset, revealed strong association between elevated expression and poor OS (p = 0.006) and poor DFS (p = 0.05), thus implicating ALDH1A1 and POU5F1 in CRC prognosis. Our data reveal distinct molecular signature of ALDH⁺ CSCs in CRC and suggest pathways relevant for successful targeted therapies and management of CRC.

MicroRNA‑195 is associated with regulating the pathophysiologic process of human laryngeal squamous cell carcinoma

MicroRNAs (miRNAs) have been reported to be associated with the modulation of tumor development, including alterations associated with the development of human laryngeal squamous cell carcinoma (LSCC). The present study was designed to investigate whether miRNA‑195 was associated with the pathophysiologic process of human LSCC and to identify its potential roles and underlying molecular mechanisms. To determine whether miRNA‑195 serves a role in LSCC, reverse transcription‑quantitative polymerase chain reaction was used to detect miRNA‑195 expression in LSCC tissues. The tumor‑suppressive effect of miRNA‑195 was determined by in vitro assays. Gain‑of‑function studies using miRNA‑195 mimics were performed to investigate cell viability, migration and invasion, and apoptosis in the AMC‑HN‑8 cell line. Western blotting was performed to reveal the molecular mechanisms of miRNA‑195 and its downstream signaling pathways in the LSCC AMC‑HN‑8 cell line. The present study demonstrated that miRNA‑195 is downregulated in primary LSCC tumors. Upregulating miRNA‑195 in vitro suppressed cell viability, migration and invasion in AMC‑HN‑8 cells. Overexpression of miRNA‑195 alone in AMC‑HN‑8 cells was sufficient to induce cell apoptosis, as identified by terminal deoxynucleotidyl transferase dUTP nick end labeling assay. Compared with the high expression of miRNA‑195 in AMC‑HN‑8 cells, the expression levels of vascular endothelial growth factor receptor‑II protein and downstream signaling pathway proteins, which were associated with cell viability, migration, invasion and apoptosis, were markedly decreased compared with control or miRNA‑195 negative control treatment group. Together, these data suggest the therapeutic potential of miRNA‑195 in modulating cell growth, migration and apoptosis during the pathophysiological progression of LSCC and that miRNA‑195 may serve as a potential therapeutic target in human LSCC.

Analysis of clock gene-miRNA correlation networks reveals candidate drivers in colorectal cancer

Altered functioning of the biological clock is involved in cancer onset and progression. MicroRNAs (miRNAs) interact with the clock genes modulating the function of genetically encoded molecular clockworks. Collaborative interactions may take place within the coding-noncoding RNA regulatory networks. We aimed to evaluate the cross-talk among miRNAs and clock genes in colorectal cancer (CRC). We performed an integrative analysis of miRNA-miRNA and miRNA-mRNA interactions on high-throughput molecular profiling of matched human CRC tissue and non-tumor mucosa, pinpointing core clock genes and their targeting miRNAs. Data obtained in silico were validated in CRC patients and human colon cancer cell lines. In silico we found severe alterations of clock gene–related coding-noncoding RNA regulatory networks in tumor tissues, which were later corroborated by the analysis of human CRC specimens and experiments performed in vitro. In conclusion, specific miRNAs target and regulate the transcription/translation of clock genes and clock gene-related miRNA-miRNA as well as mRNA-miRNA interactions are altered in colorectal cancer. Exploration of the interplay between specific miRNAs and genes, which are critically involved in the functioning of the biological clock, provides a better understanding of the importance of the miRNA-clock genes axis and its derangement in colorectal cancer.

Double-strand break repair and colorectal cancer: gene variants within 3' UTRs and microRNAs binding as modulators of cancer risk and clinical outcome

Genetic variations in 3' untranslated regions of target genes may affect microRNA binding, resulting in differential protein expression. microRNAs regulate DNA repair, and single-nucleotide polymorphisms in miRNA binding sites (miRSNPs) may account for interindividual differences in the DNA repair capacity. Our hypothesis is that miRSNPs in relevant DNA repair genes may ultimately affect cancer susceptibility and impact prognosis.In the present study, we analysed the association of polymorphisms in predicted microRNA target sites of double-strand breaks (DSBs) repair genes with colorectal cancer (CRC) risk and clinical outcome. Twenty-one miRSNPs in non-homologous end-joining and homologous recombination pathways were assessed in 1111 cases and 1469 controls. The variant CC genotype of rs2155209 in MRE11A was strongly associated with decreased cancer risk when compared with the other genotypes (OR 0.54, 95% CI 0.38-0.76, p = 0.0004). A reduced expression of the reporter gene was observed for the C allele of this polymorphism by in vitro assay, suggesting a more efficient interaction with potentially binding miRNAs. In colon cancer patients, the rs2155209 CC genotype was associated with shorter survival while the TT genotype of RAD52 rs11226 with longer survival when both compared with their respective more frequent genotypes (HR 1.63, 95% CI 1.06-2.51, p = 0.03 HR 0.60, 95% CI 0.41-0.89, p = 0.01, respectively).miRSNPs in DSB repair genes involved in the maintenance of genomic stability may have a role on CRC susceptibility and clinical outcome.