Comprehensive analysis of the functional microRNA–mRNA regulatory network identifies miRNA signatures associated with glioma malignant progression

Glioma is the most common and fatal primary brain tumour with poor prognosis; however, the functional roles of miRNAs in glioma malignant progression are insufficiently understood. Here, we used an integrated approach to identify miRNA functional targets during glioma malignant progression by combining the paired expression profiles of miRNAs and mRNAs across 160 Chinese glioma patients, and further constructed the functional miRNA–mRNA regulatory network. As a result, most tumour-suppressive miRNAs in glioma progression were newly discovered, whose functions were widely involved in gliomagenesis. Moreover, three miRNA signatures, with different combinations of hub miRNAs (regulations≥30) were constructed, which could independently predict the survival of patients with all gliomas, high-grade glioma and glioblastoma. Our network-based method increased the ability to identify the prognostic biomarkers, when compared with the traditional method and random conditions. Hsa-miR-524-5p and hsa-miR-628-5p, shared by these three signatures, acted as protective factors and their expression decreased gradually during glioma progression. Functional analysis of these miRNA signatures highlighted their critical roles in cell cycle and cell proliferation in glioblastoma malignant progression, especially hsa-miR-524-5p and hsa-miR-628-5p exhibited dominant regulatory activities. Therefore, network-based biomarkers are expected to be more effective and provide deep insights into the molecular mechanism of glioma malignant progression.     

Identification of Metastamirs as Metastasis-associated MicroRNAs in Clear Cell Renal Cell Carcinomas

MicroRNAs (miRNAs) play a pivotal role in cancerogenesis and cancer progression, but their specific role in the metastasis of clear cell renal cell carcinomas (ccRCC) is still limited. Based on microRNA microarray analyses from normal and cancerous samples of ccRCC specimens and from bone metastases of ccRCC patients, we identified a set of 57 differentially expressed microRNAs between these three sample groups of ccRCC. A selected panel of 33 miRNAs was subsequently validated by RT-qPCR on total 57 samples. Then, 30 of the 33 examined miRNAs were confirmed to be deregulated. A stepwise down-regulation of miRNA expression from normal, over primary tumor to metastatic tissue samples, was found to be typical. A total of 23 miRNAs (miR-10b/-19a/-19b/-20a/-29a/-29b/-29c/-100/-101/-126/-127/-130/-141/-143/-145/-148a/-192/-194/-200c/-210/-215/-370/-514) were down-regulated in metastatic tissue samples compared with normal tissue. This down-regulated expression in metastatic tissue in comparison with primary tumor tissue was also present in 21 miRNAs. In cell culture experiments with 5-aza-2''-deoxycytidine and trichostatin A, epigenetic modifications were shown as one reason of this down-regulation. The altered miRNA profiles, comprising newly identified metastasis-associated miRNAs, termed metastamir and the predicted miRNA-target interactions together with the significant correlations of miRNAs that were either lost or newly appeared in the studied sample groups, afford a solid basis for further functional analyses of individual miRNAs in RCC metastatic progression.     

Estimating survival time of patients with glioblastoma multiforme and characterization of the identified microRNA signatures

Background

Though glioblastoma multiforme (GBM) is the most frequently occurring brain malignancy in adults, clinical treatment still faces challenges due to poor prognoses and tumor relapses. Recently, microRNAs (miRNAs) have been extensively used with the aim of developing accurate molecular therapies, because of their emerging role in the regulation of cancer-related genes. This work aims to identify the miRNA signatures related to survival of GBM patients for developing molecular therapies.

Results

This work proposes a support vector regression (SVR)-based estimator, called SVR-GBM, to estimate the survival time in patients with GBM using their miRNA expression profiles. SVR-GBM identified 24 out of 470 miRNAs that were significantly associated with survival of GBM patients. SVR-GBM had a mean absolute error of 0.63 years and a correlation coefficient of 0.76 between the real and predicted survival time. The 10 top-ranked miRNAs according to prediction contribution are as follows: hsa-miR-222, hsa-miR-345, hsa-miR-587, hsa-miR-526a, hsa-miR-335, hsa-miR-122, hsa-miR-24, hsa-miR-433, hsa-miR-574 and hsa-miR-320. Biological analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway on the identified miRNAs revealed their influence in GBM cancer.

Conclusion

The proposed SVR-GBM using an optimal feature selection algorithm and an optimized SVR to identify the 24 miRNA signatures associated with survival of GBM patients. These miRNA signatures are helpful to uncover the individual role of miRNAs in GBM prognosis and develop miRNA-based therapies.

     

Suppression of tumorigenicity by MicroRNA-138 through inhibition of EZH2-CDK4/6-pRb-E2F1 signal loop in glioblastoma multiforme

Deregulation of microRNAs (miRNAs) is implicated in tumor progression. We attempt to indentify the tumor suppressive miRNA not only down-regulated in glioblastoma multiforme (GBM) but also potent to inhibit the oncogene EZH2, and then investigate the biological function and pathophysiologic role of the candidate miRNA in GBM. In this study, we show that miRNA-138 is reduced in both GBM clinical specimens and cell lines, and is effective to inhibit EZH2 expression. Moreover, high levels of miR-138 are associated with long overall and progression-free survival of GBM patients from The Cancer Genome Atlas dataset (TCGA) data portal. Ectopic expression of miRNA-138 effectively inhibits GBM cell proliferation in vitro and tumorigenicity in vivo through inducing cell cycles G1/S arrest. Mechanism investigation reveals that miRNA-138 acquires tumor inhibition through directly targeting EZH2, CDK6, E2F2 and E2F3. Moreover, an EZH2-mediated signal loop, EZH2-CDK4/6-pRb-E2F1, is probably involved in GBM tumorigenicity, and this loop can be blocked by miRNA-138. Additionally, miRNA-138 negatively correlates to mRNA levels of EZH2 and CDK6 among GBM clinical samples from both TCGA and our small amount datasets. In conclusion, our data demonstrate a tumor suppressive role of miRNA-138 in GBM tumorigenicity, suggesting a potential application in GBM therapy.     

LncRNA RP1-86C11.7 exacerbates the glioma progression and oncogenicity by hsa-miR-144-3p/TFRC signaling.

Glioblastoma (GBM) remains the most common and malignant tumor of the human central nervous system. Increasing evidence has highlighted that tumor cells with high transferrin receptor (TFRC) expression show advantages in growth. Long noncoding RNAs (lncRNAs) are related to glioma progression by mediating microRNAs (miRNAs). However, the underlying mechanism among TFRC, miRNA and lncRNA in GBM is limited. In the current study, we identified a new lncRNA-induced signaling mechanism that regulates the TFRC levels in GBM. The TFRC level was higher in glioma cell lines, and elevated TFRC expression promoted the proliferation and survival of glioma cells. Further study showed that hsa-miR-144a-3p bound to the 3′-UTR of TFRC mRNA and inhibited its expression, preventing the malignant properties of glioma cells, such as proliferation and survival. We also found that the lncRNA RP1-86C11.7 sponges hsa-miR-144-3p to suppress its protective role in glioma. RP1-86C11.7 overexpression in glioma cells elevated TFRC expression, increased the intracellular free iron level, and deteriorated oncogenicity, with a significant reduction in hsa-miR-144-3p. By contrast, silencing RP1-86C11.7 upregulated the hsa-miR-144-3p level, resulting in decreased TFRC expression and repressed glioma progression. However, the effect of silencing RP1-86C11.7 was reversed with simultaneous hsa-miR-144-3p inhibitor treatment: the TFRC level, intracellular iron level and proliferation in glioma cells increased. Mechanistically, our data indicated that RP1-86C11.7 exacerbates the malignant behavior of glioma through the hsa-miR-144-3p/TFRC axis. RP1-86C11.7 may be a potential biomarker or target to treat glioma in the future.     

LncRNA RP1-86C11.7 exacerbates the glioma progression and oncogenicity by hsa-miR-144-3p/TFRC signaling.

Glioblastoma (GBM) remains the most common and malignant tumor of the human central nervous system. Increasing evidence has highlighted that tumor cells with high transferrin receptor (TFRC) expression show advantages in growth. Long noncoding RNAs (lncRNAs) are related to glioma progression by mediating microRNAs (miRNAs). However, the underlying mechanism among TFRC, miRNA and lncRNA in GBM is limited. In the current study, we identified a new lncRNA-induced signaling mechanism that regulates the TFRC levels in GBM. The TFRC level was higher in glioma cell lines, and elevated TFRC expression promoted the proliferation and survival of glioma cells. Further study showed that hsa-miR-144a-3p bound to the 3′-UTR of TFRC mRNA and inhibited its expression, preventing the malignant properties of glioma cells, such as proliferation and survival. We also found that the lncRNA RP1-86C11.7 sponges hsa-miR-144-3p to suppress its protective role in glioma. RP1-86C11.7 overexpression in glioma cells elevated TFRC expression, increased the intracellular free iron level, and deteriorated oncogenicity, with a significant reduction in hsa-miR-144-3p. By contrast, silencing RP1-86C11.7 upregulated the hsa-miR-144-3p level, resulting in decreased TFRC expression and repressed glioma progression. However, the effect of silencing RP1-86C11.7 was reversed with simultaneous hsa-miR-144-3p inhibitor treatment: the TFRC level, intracellular iron level and proliferation in glioma cells increased. Mechanistically, our data indicated that RP1-86C11.7 exacerbates the malignant behavior of glioma through the hsa-miR-144-3p/TFRC axis. RP1-86C11.7 may be a potential biomarker or target to treat glioma in the future.     

MicroRNA-21 Promotes Glioblastoma Tumorigenesis by Down-regulating Insulin-like Growth Factor-binding Protein-3 (IGFBP3)

Despite advances in surgery, imaging, chemotherapy, and radiation, patients with glioblastoma multiforme (GBM), the most common histological subtype of glioma, have an especially dismal prognosis; >70% of GBM patients die within 2 years of diagnosis. In many human cancers, the microRNA miR-21 is overexpressed, and accumulating evidence indicates that it functions as an oncogene. Here, we report that miR-21 is overexpressed in human GBM cell lines and tumor tissue. Moreover, miR-21 expression in GBM patient samples is inversely correlated with patient survival. Knockdown of miR-21 in GBM cells inhibited cell proliferation in vitro and markedly inhibited tumor formation in vivo. A number of known miR-21 targets have been identified previously. By microarray analysis, we identified and validated insulin-like growth factor (IGF)-binding protein-3 (IGFBP3) as a novel miR-21 target gene. Overexpression of IGFBP3 in glioma cells inhibited cell proliferation in vitro and inhibited tumor formation of glioma xenografts in vivo. The critical role that IGFBP3 plays in miR-21-mediated actions was demonstrated by a rescue experiment, in which IGFBP3 knockdown in miR-21KD glioblastoma cells restored tumorigenesis. Examination of tumors from GBM patients showed that there was an inverse relationship between IGFBP3 and miR-21 expression and that increased IGFBP3 expression correlated with better patient survival. Our results identify IGFBP3 as a novel miR-21 target gene in glioblastoma and suggest that the oncogenic miRNA miR-21 down-regulates the expression of IGFBP3, which acts as a tumor suppressor in human glioblastoma.     

Comparative microarray analysis of microRNA expression profiles in primary cutaneous malignant melanoma, cutaneous malignant melanoma metastases, and benign melanocytic nevi

Perturbations in microRNA (miRNA) expression profiles have been reported for cutaneous malignant melanoma (CMM) predominantly when examined in cell lines. Despite the rapidly growing number of newly discovered human miRNA sequences, the availability of up-to-date miRNA expression profiles for clinical samples of primary cutaneous malignant melanoma (PCMM), cutaneous malignant melanoma metastases (CMMM), and benign melanocytic nevi (BMN) is limited. Specimens excised from the center of tumors (lesional) from patients with PCMM (n=9), CMMM (n=4), or BMN (n=8) were obtained during surgery. An exploratory microarray analysis was performed by miRNA expression profiling based on Agilent platform screening for 1205 human miRNAs. The results from the microarray analysis were validated by TaqMan quantitative real-time polymerase chain reaction. In addition to several miRNAs previously known to be associated with CMM, 19 unidentified miRNA candidates were found to be dysregulated in CMM patient samples. Among the 19 novel miRNA candidates, the genes hsa-miR-22, hsa-miR-130b, hsa-miR-146b-5p, hsa-miR-223, hsa-miR-301a, hsa-miR-484, hsa-miR-663, hsa-miR-720, hsa-miR-1260, hsa-miR-1274a, hsa-miR-1274b, hsa-miR-3663-3p, hsa-miR-4281, and hsa-miR-4286 were upregulated, and the genes hsa-miR-24-1*, hsa-miR-26a, hsa-miR-4291, hsa-miR-4317, and hsa-miR-4324 were downregulated. The results of this study partially confirm previous CMM miRNA profiling studies identifying miRNAs that are dysregulated in CMM. However, we report several novel miRNA candidates in CMM tumors; these miRNA sequences require further validation and functional analysis to evaluate whether they play a role in the pathogenesis of CMM.     

C-Myc negatively controls the tumor suppressor PTEN by upregulating miR-26a in glioblastoma multiforme cells

The c-Myc oncogene is amplified in many tumor types. It is an important regulator of cell proliferation and has been linked to altered miRNA expression, suggesting that c-Myc-regulated miRNAs might contribute to tumor progression. Although miR-26a has been reported to be upregulated in glioblastoma multiforme (GBM), the mechanism has not been established. We have shown that ectopic expression of miR-26a influenced cell proliferation by targeting PTEN, a tumor suppressor gene that is inactivated in many common malignancies, including GBM. Our findings suggest that c-Myc modulates genes associated with oncogenesis in GBM through deregulation of miRNAs via the c-Myc–miR-26a–PTEN signaling pathway. This may be of clinical relevance.     

miRNA-182 and the regulation of the glioblastoma phenotype - toward miRNA-based precision therapeutics

Glioblastoma (GBM) is an incurable cancer, with survival rates of just 14-16 months after diagnosis.¹ Functional genomics have identified numerous genetic events involved in GBM development. One of these, the deregulation of microRNAs (miRNAs), has been attracting increasing attention due to the multiple biologic processes that individual miRNAs influence. Our group has been studying the role of miR-182 in GBM progression, therapy resistance, and its potential as GBM therapeutic. Oncogenomic analyses revealed that miR-182 is the only miRNA, out of 470 miRNAs profiled by The Cancer Genome Atlas (TCGA) program, which is associated with favorable patient prognosis, neuro-developmental context, temozolomide (TMZ) susceptibility, and most significantly expressed in the least aggressive oligoneural subclass of GBM. miR-182 sensitized glioma cells to TMZ-induced apoptosis, promoted glioma initiating cell (GIC) differentiation, and reduced tumor cell proliferation via knockdown of Bcl2L12, c-Met and HIF2A.² To deliver miR-182 to intracranial gliomas, we have characterized Spherical Nucleic Acids covalently functionalized with miR-182 sequences (182-SNAs). Upon systemic administration, 182-SNAs crossed the blood-brain/blood-tumor barrier (BBB/BTB), reduced tumor burden, and increased animal subject survival.^2-4 Thus, miR-182-based SNAs represent a tool for systemic delivery of miRNAs and a novel approach for the precision treatment of malignant brain cancers.     

The mRNA–miRNA–lncRNA Regulatory Network and Factors Associated with Prognosis Prediction of Hepatocellular Carcinoma

The aim of this study was to identify novel prognostic mRNA and microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) using methods in systems biology. Differentially expressed mRNAs, miRNAs, and long non-coding RNAs (lncRNAs) were compared between HCC tumor tissues and normal liver tissues in The Cancer Genome Atlas (TCGA) database. Subsequently, a prognosis-associated mRNA co-expression network, an mRNA–miRNA regulatory network, and an mRNA–miRNA–lncRNA regulatory network were constructed to identify prognostic biomarkers for HCC through Cox survival analysis. Seven prognosis-associated mRNA co-expression modules were obtained by analyzing these differentially expressed mRNAs. An expression module including 120 mRNAs was significantly correlated with HCC patient survival. Combined with patient survival data, several mRNAs and miRNAs, including CHST4, SLC22A8, STC2, hsa-miR-326, and hsa-miR-21 were identified from the network to predict HCC patient prognosis. Clinical significance was investigated using tissue microarray analysis of samples from 258 patients with HCC. Functional annotation of hsa-miR-326 and hsa-miR-21-5p indicated specific associations with several cancer-related pathways. The present study provides a bioinformatics method for biomarker screening, leading to the identification of an integrated mRNA–miRNA–lncRNA regulatory network and their co-expression patterns in relation to predicting HCC patient survival.     

The mRNA–miRNA–lncRNA Regulatory Network and Factors Associated with Prognosis Prediction of Hepatocellular Carcinoma

The aim of this study was to identify novel prognostic mRNA and microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) using methods in systems biology. Differentially expressed mRNAs, miRNAs, and long non-coding RNAs (lncRNAs) were compared between HCC tumor tissues and normal liver tissues in The Cancer Genome Atlas (TCGA) database. Subsequently, a prognosis-associated mRNA co-expression network, an mRNA–miRNA regulatory network, and an mRNA–miRNA–lncRNA regulatory network were constructed to identify prognostic biomarkers for HCC through Cox survival analysis. Seven prognosis-associated mRNA co-expression modules were obtained by analyzing these differentially expressed mRNAs. An expression module including 120 mRNAs was significantly correlated with HCC patient survival. Combined with patient survival data, several mRNAs and miRNAs, including CHST4, SLC22A8, STC2, hsa-miR-326, and hsa-miR-21 were identified from the network to predict HCC patient prognosis. Clinical significance was investigated using tissue microarray analysis of samples from 258 patients with HCC. Functional annotation of hsa-miR-326 and hsa-miR-21-5p indicated specific associations with several cancer-related pathways. The present study provides a bioinformatics method for biomarker screening, leading to the identification of an integrated mRNA–miRNA–lncRNA regulatory network and their co-expression patterns in relation to predicting HCC patient survival.     

The Epstein-Barr virus (EBV)-induced tumor suppressor microRNA MiR-34a is growth promoting in EBV-infected B cells

Epstein-Barr virus (EBV) infection of primary human B cells drives their indefinite proliferation into lymphoblastoid cell lines (LCLs). B cell immortalization depends on expression of viral latency genes, as well as the regulation of host genes. Given the important role of microRNAs (miRNAs) in regulating fundamental cellular processes, in this study, we assayed changes in host miRNA expression during primary B cell infection by EBV. We observed and validated dynamic changes in several miRNAs from early proliferation through immortalization; oncogenic miRNAs were induced, and tumor suppressor miRNAs were largely repressed. However, one miRNA described as a p53-targeted tumor suppressor, miR-34a, was strongly induced by EBV infection and expressed in many EBV and Kaposi's sarcoma-associated herpesvirus (KSHV)-infected lymphoma cell lines. EBV latent membrane protein 1 (LMP1) was sufficient to induce miR-34a requiring downstream NF-κB activation but independent of functional p53. Furthermore, overexpression of miR-34a was not toxic in several B lymphoma cell lines, and inhibition of miR-34a impaired the growth of EBV-transformed cells. This study identifies a progrowth role for a tumor-suppressive miRNA in oncogenic-virus-mediated transformation, highlighting the importance of studying miRNA function in different cellular contexts.     

Identification of MicroRNAs in Helicoverpa armigera and Spodoptera litura Based on Deep Sequencing and Homology Analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.