Identification of prognostic markers for hepatocellular carcinoma based on miRNA expression profiles

Elucidating hepatocellular carcinoma progression: a novel prognostic miRNA-mRNA network and signature analysis

There is increasing evidence that miRNAs play an important role in the prognosis of HCC. There is currently a lack of acknowledged models that accurately predict patient prognosis. The aim of this study is to create a miRNA-based model to precisely forecast a patient's prognosis and a miRNA-mRNA network to investigate the function of a targeted mRNA. TCGA miRNA dataset and survival data of HCC patients were downloaded for differential analysis. The outcomes of variance analysis were subjected to univariate and multivariate Cox regression analyses and LASSO analysis. We constructed and visualized prognosis-related models and subsequently used violin plots to probe the function of miRNAs in tumor cells. We predicted the target mRNAs added those to the String database, built PPI protein interaction networks, and screened those mRNA using Cytoscape. The hub mRNA was subjected to GO and KEGG analysis to determine its biological role. Six of them were associated with prognosis: hsa-miR-139-3p, hsa-miR-139-5p, hsa-miR-101-3p, hsa-miR-30d-5p, hsa-miR-5003-3p, and hsa-miR-6844. The prognostic model was highly predictive and consistently performs, with the C index exceeding 0.7 after 1, 3, and 5 years. The model estimated significant differences in the Kaplan-Meier plotter and the model could predict patient prognosis independently of clinical indicators. A relatively stable miRNA prognostic model for HCC patients was constructed, and the model was highly accurate in predicting patients with good stability over 5 years. The miRNA-mRNA network was constructed to explore the function of mRNA.

MicroRNA biosensors for the detection of liver cancer

Liver cancer is one of the deadliest types worldwide and early diagnosis is highly important for successful treatment. Therefore, it is necessary to develop rapid, sensitive, simple, and inexpensive analytical tools for its detection. MicroRNAs (miRNA) represent unique biomarkers whose expression in biofluids is strongly associated with cancer in general and miR-21, -31, -122, -145, -146a, -200c, -221, -222, and -223 in liver cancer, specifically. Various biosensors for miRNA detection have been developed. These include electrochemical biosensors based on amperometric, potentiometric, conductometric and impedimetric technology. Furthermore, the use of advanced nanomaterials with enhanced chemical stability, conductivity and electrocatalytic activity have greatly increased the sensitivity and specificity of these devices. The present review focuses on recent advances in electrochemical biosensors for miRNA detection in liver cancer.

Targeting Non-Coding RNAs for the Development of Novel Hepatocellular Carcinoma Therapeutic Approaches

Hepatocellular carcinoma (HCC) remains a global health challenge, representing the third leading cause of cancer deaths worldwide. Although therapeutic advances have been made in the few last years, the prognosis remains poor. Thus, there is a dire need to develop novel therapeutic strategies. In this regard, two approaches can be considered: (1) the identification of tumor-targeted delivery systems and (2) the targeting of molecule(s) whose aberrant expression is confined to tumor cells. In this work, we focused on the second approach. Among the different kinds of possible target molecules, we discuss the potential therapeutic value of targeting non-coding RNAs (ncRNAs), which include micro interfering RNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs). These molecules represent the most significant RNA transcripts in cells and can regulate many HCC features, including proliferation, apoptosis, invasion and metastasis. In the first part of the review, the main characteristics of HCC and ncRNAs are described. The involvement of ncRNAs in HCC is then presented over five sections: (a) miRNAs, (b) lncRNAs, (c) circRNAs, (d) ncRNAs and drug resistance and (e) ncRNAs and liver fibrosis. Overall, this work provides the reader with the most recent state-of-the-art approaches in this field, highlighting key trends and opportunities for more advanced and efficacious HCC treatments.

Identification of a four-miRNA signature predicts the prognosis of papillary thyroid cancer

BACKGROUND

In recently diagnosed patients with thyroid cancer, papillary thyroid cancer (PTC), as the most common histological subtype, accounts for 90% of all cases. Although PTC is known as a relatively adolescent malignant disease, there still is a high possibility of recurrence in PTC patients with a poor prognosis. Therefore, new biomarkers are necessary to guide more effective stratification of PTC patients and personalize therapy to avoid overtreatment or inadequate treatment. Accumulating evidence demonstrates that microRNAs (miRNAs) have broad application prospects as diagnostic biomarkers in cancer.

AIM

To explore novel markers consisting of miRNA-associated signatures for PTC prognostication.

METHODS

We obtained and analyzed the data of 497 PTC patients from The Cancer Genome Atlas. The patients were randomly assigned to either a training or testing cohort.

RESULTS

We discovered 237 differentially expressed miRNAs in tumorous thyroid tissues compared with normal tissues, which contained 172 up-regulated and 65 down-regulated miRNAs. The evaluation of differently expressed miRNAs was conducted using our risk score model. We then successfully generated a four-miRNA potential prognostic signature [risk score = (-0.001 × hsa-miR-181a-2-3p) + (0.003 × hsa-miR-138-5p) + (-0.018 × hsa-miR-424-3p) + (0.284 × hsa-miR-612)], which reliably distinguished patients from high and low risk with a significant difference in the overall survival (P < 0.01) and was effective in predicting the five-year disease survival rate with the area under the receiver operating characteristic curve of 0.937 and 0.812 in the training and testing cohorts, respectively. Additionally, there was a trend indicated that high-risk patients had shorter relapse-free survival, although statistical significance was not reached (P = 0.082) in our sequencing cohort.

CONCLUSION

Our results indicated a four-miRNA signature that has a robust predictive effect on the prognosis of PTC. Accordingly, we would recommend more radical therapy and closer follow-ups for high-risk groups.

miR-139-5p and miR-451a as a Diagnostic Biomarker in LUSC

Background

Lung squamous cell carcinoma (LUSC) is a type of lung cancer that originates from segmental or subsegmental bronchial mucosa. There is evidence that miRNA plays an important role in the occurrence and progression of tumors.

Methods

In this study, plasma samples of patients with early LUSC and healthy volunteers were subjected to miRNA sequencing, and the levels of differentially expressed miRNAs (DEMs) in LUSC tissues were analyzed using R language. Cox regression and Kaplan-Meier (K-M) survival curve analyses were performed to determine the relationship between DEMs and prognosis in LUSC, and PCR method was verified for the plasma expression level of DEMs in patients with LUSC. The levels of CYFRA21-1 and SCC-Ag in plasma were measured, and area under curve (AUC) was used to evaluate the diagnostic value of the DEMs.

Results

A total of 21 DEMs were screened out by sequencing. The expression levels of DEMs in tissue samples in the TCGA database were analyzed, and four DEMs with consistent expression levels were further screened from plasma and tissue samples. Regression analysis and K-M curve were performed to select two DEMs (miR-139-5p, miR-451a) that were correlated with the prognosis. PCR verification results showed that the levels of miR-451a and miR-139-5p were low in patients, and the level of miR-139-5p in late stages III & IV with the patients of LUSC was higher than that in stages I & II. The AUC values of the four indicators (SCC-Ag, CYFRA21-1, miR-451a and miR-139-5p) in the diagnosis of LUSC, early and late cases were 0.884, 0.935 and 0.778, respectively.

Conclusion

The detection of miR-139-5p and miR-451a levels in plasma has a certain potential in the non-invasive diagnosis, especially in patients with early stages of LUSC.

The role of AURKA/miR-199b-3p in hepatocellular carcinoma cells

BACKGROUND

Previous studies proved that AURKA functions as an oncogene in several cancers. This article aimed to probe the miRNA-induced regulatory mechanism of AURKA in hepatocellular carcinoma (HCC).

METHODS

Differentially expressed genes in TCGA-LIHC dataset were screened by bioinformatics methods. Levels of miR-199b-3p and AURKA mRNA were examined by qRT-PCR. Western blot was utilized to evaluate protein levels of AURKA, p-AKT, and AKT. Dual-luciferase assay was introduced to explore their interaction. MTT, colony formation, scratch healing, transwell, and flow cytometry assays were introduced into cell proliferation, migration, invasion, and apoptosis assessment. The impact of miR-199b-3p/AURKA axis on HCC tumor growth was determined in a tumor xenograft model.

RESULTS

We found that AURKA was highly expressed in HCC and was coupled to poor prognosis of HCC. As manifested by cellular assays, compared to the normal cells HL-7702, AURKA presented notably high expression in HCC cell lines. Overexpressed AURKA evidently impelled the proliferation, colony formation, migration, and invasion of HCC cells while suppressing apoptosis. The regulatory gene upstream of AURKA was predicted to be miR-199b-3p by bioinformatics method, and there was a markedly negative correlation between the two. Overexpressed miR-199b-3p constrained HCC cell proliferation, migration, and invasion while fostering apoptosis, which could be counteracted by upregulating AURKA. MiR-199b-3p repressed the tumor growth in vivo by targeting AURKA and affected PI3K/AKT signaling pathway.

CONCLUSION

To summarize, this study implied the regulatory mechanism of miR-199b-3p/AURKA axis in HCC, and supplied optional therapeutic targets for HCC patients.

Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.

The Computational Analysis of Single Nucleotide Associated with MicroRNA Affecting Hepatitis B Infection

BACKGROUND

MicroRNAs (miRNAs) have a pivotal role in Hepatitis B Virus (HBV) infection and its complications by targeting the cellular transcription factors required for gene expression or directly binding to HBV transcripts. Single Nucleotide Polymorphisms (SNPs) in miRNA genes affect their expression and the regulation of target genes, clinical course, diagnosis, and therapeutic interventions of HBV infection.

METHODS

Computational assessment and cataloging of miRNA gene polymorphisms targeting mRNA transcripts straightly or indirectly through the regulation of hepatitis B infection by annotating the functional impact of SNPs on mRNA-miRNA and miRNA-RBS (miRNA binding sites) interaction were screened by applying various universally available datasets such as the miRNA SNP3.0 software.

RESULTS

2987 SNPs were detected in 139 miRNAs affecting hepatitis B infection. Among them, 313 SNPs were predicted to have a significant role in the progression of hepatitis B infection. The computational analysis also revealed that 45 out of the 313 SNPs were located in the seed region and were more important than others. Has-miR-139-3p had the largest number of SNPs in the seed region (n=6). On the other hand, proteoglycans in cancer, adherens junction, lysine degradation, NFkappa B signaling cascade, ECM-receptor binding, viral carcinogenesis, fatty acid metabolism, TGF-beta signaling pathway, p53 signaling pathway, immune evasion related pathways, and fatty acid biosynthesis were the most important pathways affected by these 139 miRNAs.

CONCLUSION

The results revealed 45 SNPs in the seed region of 25 miRNAs as the catalog in miRNA genes that regulated the hepatitis B infection. The results also showed the most important pathways regulated by these miRNAs that can be targeted for therapeutic purposes.

Exploring the underlying molecular mechanism of liver cancer cells under hypoxia based on RNA sequencing

Background

The aim of our study was to use the differentially expressed mRNAs (DEmRNAs) and differentially expressed miRNAs (DEmiRNAs) to illustrate the underlying mechanism of hypoxia in liver cancer.

Methods

In this study, a cell model of hypoxia was established, and autophagy activity was measured with western blotting and transmission electron microscopy. The effect of hypoxia conditions on the invasion of liver cancer cell was evaluated. RNA sequencing was used to identify DEmRNAs and DEmiRNAs to explore the mechanism of hypoxia in liver cancer cells.

Results

We found that autophagy activation was triggered by hypoxia stress and hypoxia might promote liver cancer cell invasion. In addition, a total of 407 shared DEmRNAs and 57 shared DEmiRNAs were identified in both HCCLM3 hypoxia group and SMMC-7721 hypoxia group compared with control group. Furthermore, 278 DEmRNAs and 24 DEmiRNAs were identified as cancer hypoxia-specific DEmRNAs and DEmiRNAs. Finally, we obtained 19 DEmiRNAs with high degree based on the DEmiRNA-DEmRNA interaction network. Among them, hsa-miR-483-5p, hsa-miR-4739, hsa-miR-214-3p and hsa-miR-296-5p may be potential gene signatures related to liver cancer hypoxia.

Conclusions

Our study may help to understand the potential molecular mechanism of hypoxia in liver cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01055-9.

5-Azacytidine Downregulates the Proliferation and Migration of Hepatocellular Carcinoma Cells In Vitro and In Vivo by Targeting miR-139-5p/ROCK2 Pathway

Simple Summary

For hepatocellular carcinoma (HCC), the second most common cause of cancer-related death, effective therapeutic approaches are lacking. As aberrant gene methylation is a major contributor to the development of HCC, demethylating drugs such as 5-azacytidine (5-Aza) have been proposed. However, despite the potential efficacy of 5-Aza in HCC, most of its mechanisms of action are still unknown. Here, we investigate the phenotypic/molecular effects of 5-Aza with a focus on miR-139-5p. Using multiple in vitro and in vivo models of HCC, we show for the first time that 5-Aza can impair HCC development via upregulation of miR-139-5p, which in turn downregulates the ROCK2/cyclin D1/E2F1/cyclin B1 pro-proliferative pathway and the ROCK2/MMP-2 pro-migratory pathway. These observations elucidate the mechanisms of action of 5-Aza in HCC, strengthen its therapeutic potential, and provide novel information about the crosstalk among ROCK2/cyclin D1/E2F1/cyclin B1/MMP-2 in HCC.

Abstract

Background: For hepatocellular carcinoma (HCC), effective therapeutic approaches are lacking. As aberrant gene methylation is a major contributor to HCC development, demethylating drugs such as 5-azacytidine (5-Aza) have been proposed. As most 5-Aza mechanisms of action are unknown, we investigated its phenotypic/molecular effects. Methods: 5-Aza effects were examined in the human HCC cell lines JHH-6/HuH-7 and in the rat cell-line N1-S1. We also employed a xenograft mouse model (HuH-7), a zebrafish model (JHH-6), and an orthotopic syngeneic rat model (N1-S1) of HCC. Results: 5-Aza downregulated cell viability/growth/migration/adhesion by upregulating miR-139-5p, which in turn downregulated ROCK2/cyclin D1/E2F1 and increased p27^kip1, resulting in G1/G0 cell accumulation. Moreover, a decrease in cyclin B1 and an increase in p27^kip1 led to G2/M accumulation. Finally, we observed a decrease in MMP-2 levels, a stimulator of HCC cell migration. Aza effects were confirmed in the mouse model; in the zebrafish model, we also demonstrated the downregulation of tumor neo-angiogenesis, and in the orthotopic rat model, we observed impaired N1-S1 grafting in a healthy liver. Conclusion: We demonstrate for the first time that 5-Aza can impair HCC development via upregulation of miR-139-5p, which in turn impairs the ROCK2/cyclin D1/E2F1/cyclin B1 pro-proliferative pathway and the ROCK2/MMP-2 pro-migratory pathway. Thus, we provide novel information about 5-Aza mechanisms of action and deepen the knowledge about the crosstalk among ROCK2/cyclin D1/E2F1/cyclin B1/p27^kip1/MMP-2 in HCC.

A novel nine-microRNA-based model to improve prognosis prediction of renal cell carcinoma

Background

With the improved knowledge of disease biology and the introduction of immune checkpoints, there has been significant progress in treating renal cell carcinoma (RCC) patients. Individual treatment will differ according to risk stratification. As the clinical course varies in RCC, it has developed different predictive models for assessing patient’s individual risk. However, among other prognostic scores, no transparent preference model was given. MicroRNA as a putative marker shown to have prognostic relevance in RCC, molecular analysis may provide an innovative benefit in the prophetic prediction and individual risk assessment. Therefore, this study aimed to establish a prognostic-related microRNA risk score model of RCC and further explore the relationship between the model and the immune microenvironment, immune infiltration, and immune checkpoints. This practical model has the potential to guide individualized surveillance protocols, patient counseling, and individualized treatment decision for RCC patients and facilitate to find more immunotherapy targets.

Methods

Downloaded data of RCC from the TCGA database for difference analysis and divided it into a training set and validation set. Then the prognostic genes were screened out by Cox and Lasso regression analysis. Multivariate Cox regression analysis was used to establish a predictive model that divided patients into high-risk and low-risk groups. The ENCORI online website and the results of the RCC difference analysis were used to search for hub genes of miRNA. Estimate package and TIMER database were used to evaluate the relationship between risk score and tumor immune microenvironment (TME) and immune infiltration. Based on Kaplan-Meier survival analysis, search for immune checkpoints related to the prognosis of RCC.

Results

There were nine miRNAs in the established model, with a concordance index of 0.702 and an area under the ROC curve of 0.701. Nine miRNAs were strongly correlated with the prognosis (P < 0.01), and those with high expression levels had a poor prognosis. We found a common target gene PDGFRA of hsa-miR-6718, hsa-miR-1269b and hsa-miR-374c, and five genes related to ICGs (KIR2DL3, TNFRSF4, LAG3, CD70 and TNFRSF9). The immune/stromal score, immune infiltration, and immune checkpoint genes of RCC were closely related to its prognosis and were positively associated with a risk score.

Conclusions

The established nine-miRNAs prognostic model has the potential to facilitate prognostic prediction. Moreover, this model was closely related to the immune microenvironment, immune infiltration, and immune checkpoint genes of RCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09322-9.

miR-1269a and miR-1269b: Emerging Carcinogenic Genes of the miR-1269 Family

miRNAs play an important role in the occurrence and development of human cancer. Among them, hsa-mir-1269a and hsa-mir-1269b are located on human chromosomes 4 and 17, respectively, and their mature miRNAs (miR-1269a and miR-1269b) have the same sequence. miR-1269a is overexpressed in 9 cancers. The high expression of miR-1269a not only has diagnostic significance in hepatocellular carcinoma and non-small cell lung cancer but also is related to the poor prognosis of cancer patients such as esophageal cancer, hepatocellular carcinoma, and glioma. miR-1269a can target 8 downstream genes (CXCL9, SOX6, FOXO1, ATRX, RASSF9, SMAD7, HOXD10, and VASH1). The expression of miR-1269a is regulated by three non-coding RNAs (RP11-1094M14.8, LINC00261, and circASS1). miR-1269a participates in the regulation of the TGF-β signaling pathway, PI3K/AKT signaling pathway, p53 signaling pathway, and caspase-9-mediated apoptotic pathway, thereby affecting the occurrence and development of cancer. There are fewer studies on miR-1269b compared to miR-1269a. miR-1269b is highly expressed in hepatocellular carcinoma, non-small cell lung cancer, oral squamous cell carcinoma, and pharyngeal squamous cell carcinoma, but miR-1269b is low expressed in gastric cancer. miR-1269b can target downstream genes (METTL3, CDC40, SVEP1, and PTEN) and regulate the PI3K/AKT signaling pathway. In addition, sequence mutations on miR-1269a and miR-1269b can affect their regulation of cancer. The current studies have shown that miR-1269a and miR-1269b have the potential to be diagnostic and prognostic markers for cancer. Future research on miR-1269a and miR-1269b can focus on elucidating more of their upstream and downstream genes and exploring the clinical application value of miR-1269a and miR-1269b.At present, there is no systematic summary of the research on miR-1269a and miR-1269b. This paper aims to comprehensively analyze the abnormal expression, diagnostic and prognostic value, and molecular regulatory pathways of miR-1269a and miR-1269b in multiple cancers. The overview in our work can provide useful clues and directions for future related research.

Identification of prognostic metabolic genes in adrenocortical carcinoma and establishment of a prognostic nomogram: A bioinformatic study

ABSTRACT

Adrenocortical carcinoma is an invasive malignancy with poor prognosis, high recurrence rate and limited therapeutic options. Therefore, it is necessary to establish an effective method to diagnose and evaluate the prognosis of patients, so as to realize individualized treatment and improve their survival rate.This study investigated metabolic genes that may be potential therapeutic targets for Adrenocortical carcinoma (ACC). Level 3 gene expression data from the ACC cohort and the relevant clinical information were obtained from The Cancer Genome Atlas (TCGA) database. To verify, other ACC datasets (GSE76021, GSE19750) were downloaded from the Gene Expression Omnibus (GEO) database. The ACC datasets from TCGA and GEO were used to screen metabolic genes through the Molecular Signatures Database using gene set enrichment analysis. Then, the overlapping metabolic genes of the 2 datasets were identified.A signature of five metabolic genes (CYP11B1, GSTM2, IRF9, RPL31, and UBE2C) was identified in patients with ACC. The signature could be used to divide the patients with ACC into high- and low-risk groups based on their median risk score. Multivariate Cox regression analysis was performed to determine the independent prognostic factors of ACC. Time-dependent receiver operating characteristic (ROC) curve analysis was conducted to assess the prediction accuracy of the prognostic signature. Last, a nomogram was established to assess the individualized prognosis prediction model.The results indicated that the signature of 5 metabolic genes had excellent predictive value for ACC. These findings might help improve personalized treatment and medical decisions.

Genome wide profiling of miRNAs relevant to the DNA damage response induced by hexavalent chromium exposure (DDR-related miRNAs in response to Cr (VI) exposure)

AIM

We aimed to explore the expression of miRNAs and their potential roles in the DNA damage response (DDR) induced by Cr (VI) exposure in human B lymphoblast cells (HMy2.CIR cells) and in a population of Cr (VI)-exposed humans.

METHODS

Differentially expressed miRNAs were found by a combination of miRNA sequencing and RT-qPCR validation in HMy2.CIR cells treated with K₂Cr₂O₇. Differentially expressed miRNAs related to DDR were selected for functional study. The expression levels of differential miRNAs were also investigated in chromate workers.

RESULTS

A total of 214 differentially expressed miRNAs were identified by sequencing, and the expression of 5 miRNAs among 25 associated with DDR was found to be consistent between sequencing and validation studies.Functional studies showed that miR-148a-3p, miR-21-5p, and miR-424-3p might be related to Cr (VI)-induced cell apoptosis, and miR-221-3p might participate in Cr (VI)-induced DDR. We also found that the expression of miR-21-5p and miR-424-3p was upregulated in chromate workers.

CONCLUSIONS

Cr (VI) exposure could significantly impact miRNAs expression in vitro and in chromate workers. Functional studies showed that miR-148a-3p, miR-21-5p and miR-221-3p might take a crucial role in the cellular DDR induced by Cr (VI) exposure.

Resveratrol inhibits the migration, invasion and epithelial-mesenchymal transition in liver cancer cells through up- miR-186-5p expression

To investigate the molecular mechanism of resveratrol inhibiting the metastasis of liver cancer . HepG2 and Huh7 cells were treated with different concentrations of resveratrol, and the cell viability was determined by CCK-8 assay to determine the optimal concentration of resveratrol for subsequent experiments. The expressions of miR-186-5p in liver cancer tissues and liver cancer cells were determined by quantitative real-time RT-PCR. The migration and invasion of HepG2 and Huh7 cells were detected by wound healing assay and Transwell assay, and the expression levels of epithelial-mesenchymal transition (EMT) related proteins were determined by Western blotting. Resveratrol with concentration of had no effect on the viability of HepG2 and Huh7 cells, so the concentration of resveratrol in subsequent experiments was 6.25 μmol/L. Resveratrol inhibited the wound healing and invasion of liver cancer cells; increased the expression of E-cadherin, and decreased the expression of vimentin and Twist1. The expression of miR-186-5p was significantly down-regulated in liver cancer tissues and cells compared with the adjacent tissues and normal liver cells (both <0.05). Furthermore, resveratrol induced the expression of miR-186-5p in liver cancer cells (both <0.01). Overexpression of miR-186-5p suppressed the migration, invasion and EMT of liver cancer cells. Knockdown of miR-186-5p blocked the inhibition effects of resveratrol on the migration, invasion and EMT of liver cancer cells. Resveratrol could inhibit the metastasis of liver cancer , which might be associated with up-regulating miR-186-5p.

Genome-Wide Expression Difference of MicroRNAs in Basal Cell Carcinoma

Distinct expression of the miRNAs has rarely been explored in basal cell carcinoma (BCC) of skin, and the regulatory role of miRNAs in BCC development remains quite opaque. Here, we collected control tissues from adjacent noncancerous skin (n = 15; control group) and tissues at tumor centers from patients with cheek BCC (n = 15; BCC group) using punch biopsies. After six small RNA sequencing- (sRNA-seq-) based miRNA expression profiles were generated for both BCC and controls, including three biological replicates, we conducted comparative analysis on the sRNA-seq dataset, discovering 181 differentially expressed miRNAs (DEMs) out of the 1,873 miRNAs in BCCs. In order to validate the sRNA-seq data, expression of 15 randomly selected DEMs was measured using the TaqMan probe-based quantitative real-time PCR. Functional analysis of predicted target genes of DEMs in BCCs shows that these miRNAs are primarily involved in various types of cancers, immune response, epithelial growth, and morphogenesis, as well as energy production and metabolism, indicating that BCC development is caused, at least in part, by changes in miRNA regulation for biological and disease processes. In particular, the “basal cell carcinoma pathways” were found to be enriched by predicted DEM targets, and regulatory relationships between DEMs and their targeted genes in this pathway were further uncovered. These results revealed the association between BCCs and abundant miRNA molecules that regulate target genes, functional modules, and signaling pathways in carcinogenesis.

hsa-miR-199b-3p Prevents the Epithelial-Mesenchymal Transition and Dysfunction of the Renal Tubule by Regulating E-cadherin through Targeting KDM6A in Diabetic Nephropathy

Diabetic nephropathy (DN) is the leading cause of end-stage renal disease. The association between epithelial-mesenchymal transition (EMT) and fibrosis is quite ascertained, but its link to eventual tubule dysfunction is missing. Here, we show that human microRNA- (hsa-miR-) 199b-3p protects renal tubules from diabetic-induced injury by repressing KDM6A, a histone lysine demethylase regulating E-cadherin expression. Lower E-cadherin expression is related to a higher level of KDM6A, while E-cadherin is promoted upon treatment with the KDM6A inhibitor GSK-J4 in both high glucose- (HG-) induced HK2 cells and the kidneys from streptozotocin- (STZ-) induced type 1 diabetic mice. However, overexpression or RNA silencing of E-cadherin fails to alter KDM6A expression. We also show that the upregulation of KDM6A is associated with the increased methylation level of the E-cadherin promoter. Then, the target prediction results and a dual-luciferase assay show that hsa-miR-199b-3p is a new miRNA that targets KDM6A. Overexpression of hsa-miR-199b-3p increases E-cadherin expression and prevents EMT through repressing KDM6A expression in HG-induced HK2 cells. In contrast, inhibitor-induced hsa-miR-199b-3p knockdown has opposite effects, as it decreases E-cadherin level and worsens EMT, accompanied by increased levels of KDM6A. Besides, Mir199b-knockout mice without mmu-miR-119b-3p expression exhibit more renal tubule dysfunction and more serious kidney tissue damage upon treatment with STZ. These results demonstrate that hsa-miR-199b-3p improves E-cadherin expression and prevents the progression of DN through targeting KDM6A. miR-199b-3p could be a future biomarker or target for the diagnosis or treatment of DN.

Identification of hub lncRNAs in head and neck cancer based on weighted gene co-expression network analysis and experiments

Head and neck squamous cell carcinoma (HNSCC) ranks as the sixth most common cancer among systemic malignant tumors, with 600 000 new cases occurring every year worldwide. Since HNSCC has high heterogeneity and complex pathogenesis, no effective prognostic indicator has yet been identified. Here, we aimed to identify a lncRNA signature associated with the prognosis of HNSCC as a potential new biomarker. LncRNA expression data were downloaded from The Cancer Genome Atlas database. A polygenic risk score model was constructed by using Lasso-Cox regression analysis. Weighted gene co-expression network analysis (WGCNA) was applied to analyze the co-expression modules of lncRNAs associated with the prognosis of HNSCC. The robustness of the signature was validated in testing and external cohorts. Polymerase chain reaction was performed to detect the expression levels of identified lncRNAs in cancer and adjacent tissues. We constructed an 8-lncRNA signature (LINC00567, LINC00996, MTOR-AS1, PRKG1-AS1, RAB11B-AS1, RPS6KA2-AS1, SH3BP5-AS1, ZNF451-AS1) that could be used as an independent prognostic factor of HNSCC. The signature showed strong robustness and had stable prediction performance in different cohorts. WGCNA results showed that modules related to risk score mainly participated in biological processes such as blood vessel development, positive regulation of catabolic processes, and regulation of growth. The prognostic risk score model based on lncRNA for HNSCC may help clinicians conduct individualized treatment plans.

Mir-139-5p inhibits glioma cell proliferation and progression by targeting GABRA1

Glioma is an extremely aggressive malignant neoplasm of the central nervous system. MicroRNA (miRNA) are known to bind to specific target mRNA to regulate post-transcriptional gene expression and are, therefore, currently regarded as promising biomarkers for glioma diagnosis and prognosis. The aim of the present study was to examine the pathogenesis and potential molecular markers of glioma by comparing the differential expression of miRNA and mRNA between glioma tissue and peritumor brain tissue. We explored the impact of screened core miRNA and mRNA on cell proliferation, invasion, and migration of glioma. An miRNA expression profile dataset (GSE90603) and a transcriptome profile dataset (GSE90598) were downloaded from combined miRNA-mRNA microarray chips in the Gene Expression Omnibus (GEO) database. Overall, 59 differentially expressed miRNAs (DEMs) and 419 differentially expressed genes (DEGs) were identified using the R limma software package. FunRich software was used to predict DEM target genes and miRNA-gene pairs, and Perl software was used to find overlapping genes between DEGs and DEM target genes. There were 129 overlapping genes regulated by nine miRNAs between target genes of the DEMs and DEGs. The Chinese Glioma Genome Atlas(CGGA) was analyzed in order to identify miRNAs with diagnostic and prognostic significance. MiR-139-5p, miR-137, and miR-338-3p were validated to be significantly linked to prognosis in glioma patients. Finally, we validated that miR-139-5p affected glioma malignant biological behavior via targeting gamma-aminobutyric acid A receptor alpha 1(GABRA1) through rescue experiments. Low miR-139-5p expression was correlated with survival probability and World Health Organization (WHO) grade. MiR-139-5p overexpression inhibited cell proliferation, migration, and invasion of glioma in vitro. GABRA1 was identified as a functional downstream target of miR-139-5p. Decreased GABRA1 expression was related to similar biological roles as miR-139-5p overexpression while upregulation of GABRA1 effectively reversed the inhibition effects of miR-139-5p. These results demonstrate a novel axis for miR-139-5p/GABRA1 in glioma progression and provide potential prognostic predictors and therapeutic target for glioma patients.

Exosomes function as nanoparticles to transfer miR-199a-3p to reverse chemoresistance to cisplatin in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a frequently seen malignant tumor globally. The occurrence of cisplatin (DDP) resistance is one of the main reasons for the high mortality of HCC patients. Therefore, it is of great theoretical significance and application value to explore the mechanism of chemotherapy resistance. Drug resistance can be modulated by exosomes containing mRNAs, micro RNAs (miRNAs) and other non-coding RNA (ncRNAs). Exosomal miR-199a-3p (Exo-miR-199a-3p) was subjected to extraction and verification. Whether exo-miR-199a-3p could make HCC cells sensitive to DDP in vitro was verified via flow cytometry, Cell Counting Kit-8 (CCK-8) assay, immunofluorescence assay and Transwell assay. Intravenous injection of exo-miR-199a-3p and intraperitoneal injection of DDP were carried out in vivo. Moreover, the possible targets of miR-199a-3p were screened through bioinformatics analysis, which were ascertained by Western blotting (WB). Then, miR-199a-3p levels in human normal liver epithelial cell line HL-7702 and HCC cell lines HuH7 and HuH7/DDP were elevated in a concentration-dependent manner. Exo-miR-199a-3p has abilities to adjust underlying targets and conjugate cells, to repress cells to invade, stimulate their apoptosis and abate their ability. Additionally, the caudal injection of exo-miR-199a-3p reversed the chemoresistance of tumors and slowed down their growth in the body owing to the up-regulation of miR-199a-3p and down-regulation of underlying target proteins in tumors. Finally, exo-miR-199a-3p was found to overturn the HCC’s resistance to DDP, and it may function in DDP-refractory HCC therapy as an underlying option in the future.

miR-369-3p serves as prognostic factor and regulates cancer progression of hepatocellular carcinoma

Background: The aim of this study was to investigate the role of miR-369-3p in hepatocellular carcinoma (HCC). Materials & methods: The expression levels of miR-369-3p were detected using the quantitative real-time reverse transcription-PCR analysis. The cell counting kit-8 and transwell assays were used to explore the effects of miR-369-3p on cell proliferation, migration and invasion of HCC cells. Results: The miR-369-3p expression was downregulated in HCC tissues and cell lines, in comparison to the normal controls, respectively. In vitro, overexpression of miR-369-3p in Hep 3B and Huh7 cells inhibited cell proliferation, migration and invasion. SOX4 was a direct target of miR-369-3p. Conclusion: Our results suggested that miR-369-3p may be a tumor suppressor in HCC by targeting SOX4.

Analysis and Validation of circRNA-miRNA Network in Regulating m6A RNA Methylation Modulators Reveals CircMAP2K4/miR-139-5p/YTHDF1 Axis Involving the Proliferation of Hepatocellular Carcinoma

The m⁶A RNA methylation modulators play a crucial role in regulating hepatocellular carcinoma (HCC) progression. The circular RNA (circRNA) regulatory network in regulating m⁶A RNA methylation modulators in HCC remains largely unknown. In this study, 5 prognostic m⁶A RNA methylation modulators in HCC were identified from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) projects. The differentially expressed microRNAs (DEmiRNAs) and circRNAs (DEcircRNAs) between paired tumor and normal tissues were screened out from TCGA and or Gene Expression Omnibus (GEO) database to construct the circRNA-miRNA- m⁶A RNA methylation modulator regulatory network, which included three m⁶A RNA methylation modulators (HNRNPC, YTHDF1, and YTHDF2), 11 DEmiRNAs, and eight DEcircRNAs. Among the network, hsa-miR-139-5p expression was negatively correlated with YTHDF1. Hsa-miR-139-5p low or YTHDF1 high expression was correlated with high pathological grade, advanced stage and poor survival of HCC. Additionally, cell cycle, base excision repair, and homologous recombination were enriched in YTHDF1 high expression group by GSEA. A hub circRNA regulatory network was constructed based on hsa-miR-139-5p/YTHDF1 axis. Furthermore, hsa_circ_0007456(circMAP2K4) was validated to promote HCC cell proliferation by binding with hsa-miR-139-5p to promote YTHDF1 expression. Taken together, we identified certain circRNA regulatory network related to m⁶A RNA methylation modulators and provided clues for mechanism study and therapeutic targets for HCC.

Bioinformatics Analyses of Potential miRNA-mRNA Regulatory Axis in HBV-related Hepatocellular Carcinoma

Aims: We aimed to explore the crucial miRNA-mRNA axis through bioinformatics analysis and provide evidences for the development of pathophysiological mechanisms and new therapies for HBV-related HCC. Methods: MiRNA (GSE76903) and mRNA (GSE77509) dataset were used to screen differentially expressed miRNAs (DE-miRNAs) and differentially expressed mRNAs (DE-mRNAs) using R software. Overlapping genes between DE-mRNAs and target genes of DE-miRNAs were identified as candidate genes. Hub genes were obtained via cytohubba analysis. The expression at protein and mRNA levels and prognostic value of hub genes were evaluated based on The Cancer Genome Atlas (TCGA) data. Key miRNA-mRNA axes were constructed according to predicted miRNA-mRNA pairs. MiRNA expression and prognostic role were respectively identified using starBase v3.0 and Kaplan-Meier plotter database. Real-time PCR was performed to verify the expression of crucial miRNAs and mRNAs. Coexpression of crucial miRNA and mRNA were analyzed using starBase v3.0. Results: CDK1, CCNB1, CKS2 and CCNE1 were screened as hub genes, which were significantly upregulated at protein and mRNA levels. These up-regulated hub genes were also significantly associated with poor prognosis. Hsa-mir-195-5p/CDK1, hsa-mir-5589-3p/CCNB1 and hsa-let-7c-3p/CKS2 were screened as critical miRNA-mRNA axes. Critical miRNAs were decreased in HCC, which indicates unfavourable prognosis. QPCR results showed that crucial miRNAs were decreased, whereas critical mRNAs were increased in HBV-related HCC. A reverse relationship between miRNA and mRNA in crucial axis was further verified. Conclusion: This study identified several miRNA-mRNA axes in HBV-related HCC. Hsa-mir-195-5p/CDK1, hsa-mir-5589-3p/CCNB1 and hsa-let-7c-3p/CKS2 might serve as potential prognostic biomarkers and therapeutic targets for HBV-related HCC.

Identifying miRNA modules associated with progression of keloids through weighted gene co-expression network analysis and experimental validation in vitro

Keloid is a type of skin fibroproliferative disease, characterized by excessive deposition of collagen in the extracellular matrix, myofibroblast activation and invasive growth to the surrounding normal skin tissue. However, the specific pathogenesis of keloids is not yet fully understood and existing treatment strategies are unsatisfied. It is therefore urgent to explore new biomarkers associated with its progression for keloids. In this study, the microarray dataset GSE113620 was downloaded from the Gene Expression Omnibus (GEO) database to screen out the differential expression of miRNAs (DEMs). The DEMs with large variance were applied to construct a weighted gene co-expression network to identify miRNA modules that are closely relevant to keloid progression. It is worth noting that miR-424-3p in the blue module (r = 0.98, p = 1e-18) is considered to be the ultimate target most relevant to keloid progression through co-expressed network analysis. Subsequently, the results of molecular biology experiments determine that miR-424-3p targeting Smad7 significantly enhanced the ability of cell proliferation, migration and collagen secretion after transfection with miR-424-3p mimic, while the apoptosis rate was significantly reduced. On the contrary, the miR-424-3p inhibitor performs the exact opposite function.

Identification of modules and novel prognostic biomarkers in liver cancer through integrated bioinformatics analysis

Liver cancer is a common malignant tumor with poor prognosis. Due to the lack of specific clinical manifestations at early stages, most patients are already at advanced stages of the disease by the time of diagnosis. Identification of novel biomarkers for liver cancer may thus enable earlier detection, improving outcome. MicroRNAs (miRNAs) are small endogenous noncoding RNAs of 18–22 nucleotides in length, which have a regulatory role in the expression of target proteins. Increased evidence suggests that miRNAs are abnormally expressed in a variety of cancer malignancies. Here, we combined RNA sequencing data and clinical information from The Cancer Genome Atlas Liver Hepatocellular Carcinoma database for weighted gene coexpression network analysis to identify potential miRNA prognostic biomarkers. We constructed nine coexpression modules, allowing us to identify that miR‐105‐5p, miR‐767‐5p, miR‐1266‐5p, miR‐4746‐5p, miR‐500a‐3p, miR‐1180‐3p and miR‐139‐5p are significantly associated with liver cancer prognosis. We found that these miRNAs exhibit significant association with prognosis of patients with liver cancer and confirmed the expression of these miRNAs in liver cancer tissues. Multivariate Cox regression analysis showed that miR‐105‐5p and miR‐139‐5p may be considered as independent factors. In summary, here we report that seven miRNAs have potential value as prognostic biomarkers of liver cancer.

Construction and analysis of macrophage infiltration related circRNA-miRNA-mRNA regulatory networks in hepatocellular carcinoma

BACKGROUND

Macrophage play a crucial role in regulating tumor progression. This study intended to investigate the circular RNA (circRNA) regulatory network associated with macrophage infiltration in hepatocellular carcinoma (HCC).

METHODS

The immune cell fractions of HCC from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium were calculated by Estimation of the Proportion of Immune and Cancer cells algorithm. The differentially expressed mRNAs (DEmRNAs), microRNAs (DEmiRNAs) and circRNAs (DEcircRNAs) were identified from HCC and adjacent non-tumor cases of TCGA or Gene Expression Omnibus database. The DEmRNAs related to macrophage were selected by weighted gene co-expression network analysis and then utilized to generate the circRNA-miRNA-mRNA network. A hub circRNA regulatory network was established based on the co-expressed DEmiRNAs and DEmRNAs owning contrary correlation with the clinical characteristics, survival and macrophage infiltration level. A gene signature based on the DEmRNAs in hub network was also generated for further evaluation. The circRNA binding bite for miRNA was detected by luciferase assay.

RESULTS

High macrophage fraction predicted good survival for HCC. A circRNA-miRNA-mRNA network was constructed by 27 macrophage related DEmRNAs, 21 DEmiRNAs, and 15 DEcircRNAs. Among this network, the expression of hsa-miR-139-5p was negatively correlated with CDCA8, KPNA2, PRC1 or TOP2A. Hsa-miR-139-5p low or targeted DEmRNA high expression was associated with low macrophage infiltration, high grade, advanced stage and poor prognosis of HCC. Additionally, the risk score generated by 4-DEmRNA signature could reflect the macrophage infiltration status and function as an independent prognostic factor for HCC. Finally, hsa_circ_0007456 acting on hsa-miR-139-5p related network was viewed as the hub circRNA regulatory network. Taken together, some circRNA regulatory networks may be associated with macrophage infiltration, which provides clues for mechanism study and therapeutic strategies of HCC.

Hsa_circ_0003645 Promotes Breast Cancer Progression by Regulating miR-139-3p/HMGB1 Axis

Background

The aim of the present study was to investigate the effect of over-expressing circular RNA (circ_0003645) on cell functions and its molecular mechanism in breast cancer.

Methods

The expression profile of circ_0003645, breast cancer cell lines, and the transcription levels of circular RNA, miRNA and HMGB1 gene were detected by qRT-PCR. Flow cytometry analysis was manipulated to evaluate cancer cell proliferation and cell apoptosis. The correlation between miR-139p-3p and circular_0003645 or HMGB1 was predicted by GEO, and TCGA was confirmed using the dual-luciferase reporter assay.

Results

Circ_0003645 expression was conspicuously increased in both the breast cancer tissues and cell lines. Circ_0003645 knockdown inhibited cell proliferation and induced the apoptosis of breast cancer cells in vitro and in vivo. By sponging miR-139-3p, circ_0003645 promoted the breast cancer cells progression and positively regulated HMGB1 gene.

Conclusion

Circ_0003645 functions as a ceRNA for miR-139-3p, which could upregulate HMGB1 and further promote cell proliferation in breast cancer.

Identification and characterization of miRNA expression profiles across five tissues in giant panda

microRNA (miRNA) is a small endogenous noncoding RNA molecule that plays multiple roles in regulating most biological processes. However, for China's national treasure giant panda, a world-famous rare and protected species, reports of its miRNA have been found only in blood and breast milk. To explore the miRNA expression differences between different giant panda tissues, here, we generated the miRNA profiles of five tissues (heart, liver, spleen, lung, and kidney) from four giant pandas with Illumina Hiseq 2500 platform, and filtered the differentially expressed miRNAs (DEmiRs) in each tissue, predicted the target genes of miRNA from each tissue based on the DEmiRs. Then, the GO and KEGG enrichment analysis were conducted using the target genes predicted from DEmiRs in each tissue. The RNA-seq generated an average of 0.718 GB base per sample. A total of 1,256 known miRNAs and 12 novel miRNAs were identified, and there were 215, 131, 185, 83, and 126 tissue-specific DEmiRs filtered in the heart, liver, spleen, lung, and kidney, respectively, including miR-1b-5p, miR-122-5p, miR-143, miR-126-5p, and miR-10b-5p, respectively. The predicted target genes, including MYL2, LRP5, MIF, CFD, and PEBP1 in the heart, liver, spleen, lung, and kidney, respectively, were closely associated with tissue-specific biological functions. The enrichment analysis results of target genes showed tissue-specific characteristics, such as the significantly enriched GO terms extracellular matrix in the heart and insulin-like growth factor binding in the liver. The miRNA profiles of the heart, liver, spleen, lung, and kidney of giant panda have been reported in this study, it reveals the miRNA expression differences between different tissues of the giant panda, and provides valuable genetic resources for the further related molecular genetic research of the rare and protected species giant panda and other mammals.

MDEHT: a multivariate approach for detecting differential expression of microRNA isoform data in RNA-sequencing studies

Motivation

miRNA isoforms (isomiRs) are produced from the same arm as the archetype miRNA with a few nucleotides different at 5 and/or 3 termini. These well-conserved isomiRs are functionally important and have contributed to the evolution of miRNA genes. Accurate detection of differential expression of miRNAs can bring new insights into the cellular function of miRNA and a further improvement in miRNA-based diagnostic and prognostic applications. However, very few methods take isomiR variations into account in the analysis of miRNA differential expression.

Results

To overcome this challenge, we developed a novel approach to take advantage of the multidimensional structure of isomiR data from the same miRNAs, termed as a multivariate differential expression by Hotelling’s T² test (MDEHT). The utilization of the information hidden in isomiRs enables MDEHT to increase the power of identifying differentially expressed miRNAs that are not marginally detectable in univariate testing methods. We conducted rigorous and unbiased comparisons of MDEHT with seven commonly used tools in simulated and real datasets from The Cancer Genome Atlas. Our comprehensive evaluations demonstrated that the MDEHT method was robust among various datasets and outperformed other commonly used tools in terms of Type I error rate, true positive rate and reproducibility.

Availability and implementation

The source code for identifying and quantifying isomiRs and performing miRNA differential expression analysis is available at https://github.com/amanzju/MDEHT.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comprehensive analysis of differentially expressed microRNAs and mRNAs involved in diabetic corneal neuropathy

AIMS

Corneal nerve fibers are derived from the ophthalmic division of the trigeminal ganglion (TG). Here, by sequencing of microRNAs (miRNAs) and messenger RNAs (mRNAs) from diabetic and normal TG tissues, we aimed to uncover potential miRNAs, mRNAs, and the network of their interactions involved in the pathogenesis of diabetic corneal neuropathy.

MAIN METHODS

We performed RNA sequencing to systematically screen out differentially expressed miRNAs and mRNAs in TG tissues from diabetic and normal mice. Functional enrichment analyses were performed to illustrate the biological functions of differentially expressed mRNAs (DEmRNAs). Following this, miRNA-mRNA regulatory networks were built by means of bioinformatics methods to suggest regulatory role for miRNAs in the pathogenesis of diabetic corneal neuropathy. Finally, the credibility of the sequencing-based results was validated using qRT-PCR.

KEY FINDINGS

Sequencing analyses disclosed that 68 miRNAs and 114 mRNAs were differentially expressed in diabetic TG tissues compared with normal TG samples. The functional analyses showed that DEmRNAs participated in diabetes-related biological processes. After applying an optimized approach to predict miRNA-mRNA pairs, a miRNA-mRNA interacting network was inferred. Subsequently, the expression and correlation of miR-350-5p and Mup20, miR-592-5p and Angptl7 as well as miR-351-5p and Elovl6 were preliminarily validated.

SIGNIFICANCE

Our study provides a systematic characterization of miRNA and mRNA expression in the TG during diabetic corneal neuropathy and will contribute to the development of clinical diagnostic and therapeutic strategies for diabetic corneal neuropathy.

ExoceRNA atlas: A database of cancer ceRNAs in human blood exosomes

AIMS

Competing endogenous RNAs (ceRNAs) play essential roles in cancer pathogenesis and those in exosomes have been the promising biomarkers for cancer diagnose and therapy. We aim to identify potential active ceRNA pairs in cancer blood exosomes by combining TCGA and exoRBase.

MAIN METHODS

Two strict screening criteria were implemented, including hypergeometric test on the targets predicted by RNA22 for differential miRNAs and Pearson test on the candidate mRNAs and lncRNAs for each cancer. Then2638292, 4925485 and 70669 ceRNAs in blood exosomes are available for colorectal cancer (CRC), hepatocellular carcinoma (HCC) and pancreatic adenocarcinoma (PAAD), respectively.

KEY FINDINGS

A comprehensive functional analysis on differential miRNAs in cancer blood exosomes indicates that they play important roles in development of cancer by degrading or inhibiting the post-transcription translation level of mRNA or by acting as mediators to regulate the expression of mRNA. Topological and biological functional analysis of ceRNA networks demonstrate that hub ceRNAs involve in cancer-related biological pathways and processes, so as to influence the occurrence and development of cancer and would be the potential biomarkers for three cancers. Finally, we designed a web-accessible database, ExoceRNA Atlas (https://www.exocerna-atlas.com/exoceRNA#/) as a repository of ceRNAs in blood exosomes. It can friendly search, browse and visualize ceRNA networks of the query genes along with giving the detailed functional analysis results. The entire ceRNA data can also be freely downloaded.

SIGNIFICANCE

ExoceRNA Atlas will serve as a powerful public resource for identifying ceRNAs and greatly deepen our understanding their functions in cancer exosomes.

Comprehensive Transcriptomic Analysis and Experimental Validation Identify lncRNA HOXA-AS2/miR-184/COL6A2 as the Critical ceRNA Regulation Involved in Low-Grade Glioma Recurrence

Purpose

The recurrence and metastasis of glioma are closely related to complex regulatory networks among protein-coding genes, lncRNAs and microRNAs. The aim of this study was to investigate core genes, lncRNAs, miRNAs and critical ceRNA regulatory mechanisms, which are involved in lower-grade glioma (LGG) recurrence.

Materials and Methods

We employed multiple datasets from Chinese Glioma Genome Atlas (CGGA) database and The Cancer Genome Atlas (TCGA) to perform comprehensive transcriptomic analysis. Further in vitro experiments including cell proliferation assay, luciferase reporter assay, and Western blot were performed to validate our results.

Results

Recurrent LGG and glioblastoma (GBM) showed different transcriptome characteristics with less overlap of differentially expressed protein-coding genes (DEPs), lncRNAs (DELs) and miRNAs (DEMs) compared with primary samples. There were no overlapping gene in ontology (GO) terms related to GBM recurrence in the TCGA and CGGA databases, but there were overlaps associated with LGG recurrence. GO analysis and protein–protein interaction (PPI) network analysis identified three core genes: TIMP1, COL1A1 and COL6A2. By hierarchical cluster analysis of them, LGGs could be clustered as Low_risk and High_risk group. The High_risk group with high expression of TIMP1, COL1A1, and COL6A2 showed worse prognosis. By coexpression networks analysis, competing endogenous RNA (ceRNA) network analysis, cell proliferation assay and luciferase reporter assay, we confirmed that lncRNA HOXA-AS2 functioned as a ceRNA for miR-184 to regulate expression of COL6A2, which induced cell proliferation of low-grade glioma.

Conclusion

In this study, we revealed a 3-hub protein-coding gene signature to improve prognostic prediction in LGG, and identified a critical ceRNA regulation involved in LGG recurrence.

Diagnostic, prognostic, and therapeutic significance of miR-139-5p in cancers

miRNAs are a group of non-coding RNAs that have regulatory functions in post-transcriptional gene expression. These molecules play a fundamental role in cellular processes, for instance cell proliferation, apoptosis, migration, and invasion. Scientific investigations have previously established that miRNAs can either promote or suppress tumor development by mediating different signaling pathways. miR-139-5p, located on chromosome 11q13.4, has been examined extensively in cancers. Studies have demonstrated that miR-139-5p might be an attractive cancer biomarker. Herein, we will review how miR-139-5p acts in cancer diagnosis, prognosis, and therapy, as well as elucidating its major target genes and associated signaling pathways.

The novel miR-1269b-regulated protein SVEP1 induces hepatocellular carcinoma proliferation and metastasis likely through the PI3K/Akt pathway

Decreased intercellular adhesion is a key step in the metastasis and recurrence of many cancers, including hepatocellular carcinoma (HCC). SVEP1 is an important cell adhesion molecule that plays a key role in regulating intercellular adhesion and embryonic lymphatic development. However, the expression patterns and roles of SVEP1 in HCC are still largely unknown. We identified SVEP1 expression by analyzing 220 HCC samples from our cancer center. TCGA and GEO online-databases were used for data calibration and validation. SVEP1 was differentially expressed in two groups of HCCs with different risks of recurrence and was deemed as an independent risk factor for the prognosis of HCC. The expression of SVEP1 is negatively related to the proliferation and metastasis of HCC. Downregulation of SVEP1 expression promoted in vitro HCC cell migration, chemotaxis, invasion and proliferation, as well as in vivo tumor growth, local invasion and metastasis in a mouse model. Bioinformatic analysis and RT-PCR results showed that miR-1269b expression is negatively correlated with the SVEP1 expression and the prognosis of HCC patients. Further experiments showed that miR-1269b directly targets and downregulates the expression of SVEP1, which further induces the phosphorylation of Akt at thr308. These regulatory effects ultimately mediate the proliferation and metastasis of HCC cells. SVEP1 could serve as a promising prognostic marker of HCC. MiR-1269b downregulates SVEP1 expression and promotes HCC proliferation and metastasis likely through the PI3k/Akt signaling pathway.

miR-487b and TRAK2 that form an axis to regulate the aggressiveness of osteosarcoma, are potential therapeutic targets and prognostic biomarkers

To investigate the effect of microRNA-487b (miR-487b) as well as the underlying mechanism in osteosarcoma (OS). Data downloaded from the Gene Expression Omnibus (GEO) database were used to analyze the expression and prognostic value of miR-487b/TRAK2. Cell counting kit-8, colony formation, and transwell assays were performed to investigate the biological functions of miR-487b and TRAK2. Luciferase reporter assay was applied to confirm the interactions between miR-487b and TRAK2. miR-487b was overexpressed in OS tissues and was inversely associated with the prognosis of OS patients. We discovered that miR-487b could contribute to the proliferative, clonogenic, invasive, and migratory capabilities of OS cells. Through target prediction using miRWalk and differential expression analysis based on the GEO data set, trafficking kinesin protein 2 (TRAK2) was recognized as a potential target of miR-487b, which was further verified by luciferase reporter assay. The expression of TRAK2 was decreased in OS tissues compared with normal tissues and was positively correlated with the prognosis of OS patients. A negative relevance was presented between the expression of miR-487b and TRAK2 in OS cells. Of note, further mechanistic analyses indicated that TRAK2 was implicated in the regulatory effect of miR-487b on the cell malignant behaviors in OS. To sum up, these results demonstrated that miR-487b played an oncogenic role in OS progression via directly targeting TRAK2, which could advance the development of cancer treatment.

circKIF4A promotes tumorogenesis of glioma by targeting miR-139-3p to activate Wnt5a signaling

Background

Glioma has the characteristics of high incidence and mortality, and is a common malignant tumor of the central nervous system. Circular RNAs (circRNAs) have been reported to play vital roles in progression of cancer including glioma, and circKIF4A is up-regulated in glioma tissues. However, its role and mechanisms in gliomas are unclear.

Methods

circKIF4A and miR-139-3p were determined by qRT-PCR. Transwell assay, wound-healing assay, cell colony formation and flow cytometry were performed to measure cell invasion, migration, proliferation and apoptosis. Western blotting was used to evaluate Wnt/β-catenin pathway-related protein. Luciferase reporter assays confirmed the relationship among circKIF4A, miR-139-3p and Wnt5a. Sphere formation was performed to measure the ability of glioma-initiating cells (GICs) spheroid formation. A nude mouse xenograft model was established and immunohistochemical staining was used to detect Ki-67 and Wnt5a levels.

Results

circKIF4A and Wnt5a were up-regulated and miR-139-3p was down-regulated in both glioma cells and tissues. circKIF4A promoted Wnt5a expression by sponging miR-139-3p. Knockdown of circKIF4A inhibited the colony formation ability, migration and invasion, and promoted the apoptosis of glioma cells by regulating miR-139-3p. Knockdown of circKIF4A inhibited Wnt/β-catenin signaling pathway and proliferation-related signal via miR-139-3p. Furthermore, knockdown of circKIF4A or overexpression of miR-139 suppressed the ability of sphere formation of GICs and inhibitd Wnt/β-catenin signaling pathway and proliferation-related signal in GICs. Additionally, depletion of circKIF4A decreased the expression level of Wnt5a and Ki-67, inhibited tumorigenesis in xenograft modes.

Conclusion

circKIF4A was overexpressed in glioma, and knockdown of circKIF4A suppressed glioma progression via miR-139-3p/Wnt5a axis. The results indicated that circKIF4A may be a potential target for clinical treatment of glioma.

Transcriptome analysis revealed key prognostic genes and microRNAs in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide. However, the molecular mechanisms involved in HCC remain unclear and are in urgent need of elucidation. Therefore, we sought to identify biomarkers in the prognosis of HCC through an integrated bioinformatics analysis.

Methods

Messenger RNA (mRNA) expression profiles were obtained from the Gene Expression Omnibus database and The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) for the screening of common differentially expressed genes (DEGs). Function and pathway enrichment analysis, protein-protein interaction network construction and key gene identification were performed. The significance of key genes in HCC was validated by overall survival analysis and immunohistochemistry. Meanwhile, based on TCGA data, prognostic microRNAs (miRNAs) were decoded using univariable and multivariable Cox regression analysis, and their target genes were predicted by miRWalk.

Results

Eleven hub genes (upregulated ASPM, AURKA, CCNB2, CDC20, PRC1 and TOP2A and downregulated AOX1, CAT, CYP2E1, CYP3A4 and HP) with the most interactions were considered as potential biomarkers in HCC and confirmed by overall survival analysis. Moreover, AURKA, PRC1, TOP2A, AOX1, CYP2E1, and CYP3A4 were considered candidate liver-biopsy markers for high risk of developing HCC and poor prognosis in HCC. Upregulation of hsa-mir-1269b, hsa-mir-518d, hsa-mir-548aq, hsa-mir-548f-1, and hsa-mir-6728, and downregulation of hsa-mir-139 and hsa-mir-4800 were determined to be risk factors of poor prognosis, and most of these miRNAs have strong potential to help regulate the expression of key genes.

Conclusions

This study undertook the first large-scale integrated bioinformatics analysis of the data from Illumina BeadArray platforms and the TCGA database. With a comprehensive analysis of transcriptional alterations, including mRNAs and miRNAs, in HCC, our study presented candidate biomarkers for the surveillance and prognosis of the disease, and also identified novel therapeutic targets at the molecular and pathway levels.