A MiRNA Signature for Defining Aggressive Phenotype and Prognosis in Gliomas

Gliomas represent a disparate group of tumours for which there are to date no cure. Thus, there is a recognized need for new diagnostic and therapeutic approaches based on increased understanding of their molecular nature. We performed the comparison of the microRNA (miRNA) profile of 8 WHO grade II gliomas and 24 higher grade tumours (2 WHO grade III and 22 glioblastomas) by using the Affymetrix GeneChip miRNA Array v. 1.0. A relative quantification method (RT-qPCR) with standard curve was used to confirm the 22 miRNA signature resulted by array analysis. The prognostic performances of the confirmed miRNAs were estimated on the Tumor Cancer Genome Atlas (TCGA) datasets. We identified 22 miRNAs distinguishing grade II gliomas from higher grade tumours. RT-qPCR confirmed the differential expression in the two patients'' groups for 13 out of the 22 miRNAs. The analysis of the Glioblastoma Multiforme (GBM) and Lower Grade Glioma (LGG) datasets from TCGA demonstrated the association with prognosis for 6 of those miRNAs. Moreover, in the GBM dataset miR-21 and miR-210 were predictors of worse prognosis in both univariable and multivariable Cox regression analyses (HR 1.19, p = 0.04, and HR 1.18, p = 0.029 respectively). Our results support a direct contribution of miRNAs to glioma cancerogenesis and suggest that miR-21 and miR-210 may play a role in the aggressive clinical behaviour of glioblastomas.     

Involvement of MicroRNAs in Infection of Silkworm with Bombyx mori Cytoplasmic Polyhedrosis Virus (BmCPV)

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is one of the most important pathogens of silkworm. MicroRNAs (miRNAs) have been demonstrated to play key roles in regulating host-pathogen interaction. However, there are limited reports on the miRNAs expression profiles during insect pathogen challenges. In this study, four small RNA libraries from BmCPV-infected midgut of silkworm at 72 h post-inoculation and 96 h post-inoculation and their corresponding control midguts were constructed and deep sequenced. A total of 316 known miRNAs (including miRNA*) and 90 novel miRNAs were identified. Fifty-eight miRNAs displayed significant differential expression between the infected and normal midgut (P value < = 0.01 and fold change > = 2.0 or < = 0.5), among which ten differentially expressed miRNA were validated by qRT-PCR method. Further bioinformatics analysis of predicted target genes of differentially expressed miRNAs showed that the miRNA targets were involved in stimulus and immune system process in silkworm.     

High Bone Sialoprotein (BSP) Expression Correlates with Increased Tumor Grade and Predicts a Poorer Prognosis of High-Grade Glioma Patients

Objectives

To investigate the expression and prognostic value of bone sialoprotein (BSP) in glioma patients.

Methods

We determined the expression of BSP using real-time RT-PCR and immunohistochemistry in tissue microarrays containing 15 normal brain and 270 glioma samples. Cumulative survival was calculated by the Kaplan-Meier method and analyzed by the log-rank test. Univariate and multivariate analyses were performed by the stepwise forward Cox regression model.

Results

Both BSP mRNA and protein levels were significantly elevated in high-grade glioma tissues compared with those of normal brain and low-grade glioma tissues, and BSP expression positively correlated with tumor grade (P<0.001). Univariate and multivariate analysis showed high BSP expression was an independent prognostic factor for a shorter progression-free survival (PFS) and overall survival (OS) in both grade III and grade IV glioma patients [hazard ratio (HR) = 2.549 and 3.154 for grade III glioma, and HR = 1.637 and 1.574 for grade IV glioma, respectively]. Patients with low BSP expression had a significantly longer median OS and PFS than those with high BSP expression. Small extent of resection and lineage of astrocyte served as independent risk factors of both shorter PFS and OS in grade III glioma patients; GBM patients without O⁶-methylguanine (O⁶-meG) DNA methyltransferase (MGMT) methylation and Karnofsky performance score (KPS) less than 70 points were related to poor prognosis. Lack of radiotherapy related to shorter OS but not affect PFS in both grade III and grade IV glioma patients.

Conclusion

High BSP expression occurs in a significant subset of high-grade glioma patients and predicts a poorer outcome. The study identifies a potentially useful molecular marker for the categorization and targeted therapy of gliomas.     

High Cytoplasmic FOXO1 and pFOXO1 Expression in Astrocytomas Are Associated with Worse Surgical Outcome

FOXO1 is at a convergence point of receptor tyrosine kinase (RTK) signaling, which is one of the three core pathways implicated in glioblastoma. It was recently shown that FOXO1 can effectively induce glioma cell death and inhibit tumor growth through cell cycle arrest and apoptosis. We therefore evaluated FOXO1 and pFOXO1 protein expression in 181 primary astrocytoma samples and 16 normal brain samples. Astrocytoma samples expressed higher cytoplasmic FOXO1 and pFOXO1 than normal brain samples. Nuclear pFOXO1 level was significantly higher than nuclear FOXO1 in astrocytomas. High cytoplasmic FOXO1 expression was associated with older onset age (P = 0.001) and higher WHO grade (P = 0.001). The trend was also observed between cytoplasmic pFOXO1 expression and WHO grade although not significant. Univariate survival analysis showed that both high cytoplasmic FOXO1 and pFOXO1 expression indicated a significantly shorter median overall survival and progression-free survival. Multivariate survival analysis revealed cytoplasmic FOXO1 expression, cytoplasmic pFOXO1 expression, WHO grade, gender, extent of resection and radiotherapy to be independent prognostic factors for overall survival and progression-free survival. Thus, our data suggested that cytoplasmic FOXO1 and pFOXO1 expression may serve as valuable prognostic variables in astrocytomas and may have significant implications for the development and application of targeted therapy.     

Integrated miRNA-mRNA Analysis Revealing the Potential Roles of miRNAs in Chordomas

Introduction

Emerging evidence suggests that microRNAs (miRNAs) are crucially involved in tumorigenesis and that paired expression profiles of miRNAs and mRNAs can be used to identify functional miRNA-target relationships with high precision. However, no studies have applied integrated analysis to miRNA and mRNA profiles in chordomas. The purpose of this study was to provide insights into the pathogenesis of chordomas by using this integrated analysis method.

Methods

Differentially expressed miRNAs and mRNAs of chordomas (n = 3) and notochord tissues (n = 3) were analyzed by using microarrays with hierarchical clustering analysis. Subsequently, the target genes of the differentially expressed miRNAs were predicted and overlapped with the differentially expressed mRNAs. Then, GO and pathway analyses were performed for the intersecting genes.

Results

The microarray analysis indicated that 33 miRNAs and 2,791 mRNAs were significantly dysregulated between the two groups. Among the 2,791 mRNAs, 911 overlapped with putative miRNA target genes. A pathway analysis showed that the MAPK pathway was consistently enriched in the chordoma tissue and that miR-149-3p, miR-663a, miR-1908, miR-2861 and miR-3185 likely play important roles in the regulation of MAPK pathways. Furthermore, the Notch signaling pathway and the loss of the calcification or ossification capacity of the notochord may also be involved in chordoma pathogenesis.

Conclusion

This study provides an integrated dataset of the miRNA and mRNA profiles in chordomas, and the results demonstrate that not only the MAPK pathway and its related miRNAs but also the Notch pathway may be involved in chordoma development. The occurrence of chordoma may be associated with dysfunctional calcification or ossification of the notochord.     

Comprehensive analysis of miRNA-gene regulatory network with clinical significance in human cancers

microRNAs(miRNAs), particularly the exosomal miRNAs have been widely used as biomarkers and promising therapeutic targets in cancer. However, a comprehensive analysis of miRNA-gene regulatory network with clinical significance remains scarce. The emergence of high-throughput multi-omics data over large, well-characterized patient cohorts provides an unprecedented opportunity to address this problem. Herein, we performed a clinic-centered analysis to identify cancer-associated miRNAs, miRNA-target axis. We first calculated the correlation among miRNA, mRNA and 75 unique clinico-pathological characteristics(CPCs) in 26 cancer types, and established an online resource(4CR). Interestingly, we found that the high expression of several DNA methylation-related enzymes was associated with adverse outcomes of cancer patients, and these genes were regulated by a cluster of miRNAs. Furthermore, by integrating exosomal miRNA and m RNA databases, we identified exosomal miRNA biomarkers for non-invasive cancer surveillance and therapy monitoring. Finally, we explored the role of CPC-related miRNAs for therapeutic effect prediction of drugs based on their shared targets. Our analysis pipeline illustrated the significance of clinic-centered analysis in miRNA-gene pair identification and provided helpful clues for future cancer studies.     

Glioma IL13Rα2 Is Associated with Mesenchymal Signature Gene Expression and Poor Patient Prognosis

A major challenge for successful immunotherapy against glioma is the identification and characterization of validated targets. We have taken a bioinformatics approach towards understanding the biological context of IL-13 receptor α2 (IL13Rα2) expression in brain tumors, and its functional significance for patient survival. Querying multiple gene expression databases, we show that IL13Rα2 expression increases with glioma malignancy grade, and expression for high-grade tumors is bimodal, with approximately 58% of WHO grade IV gliomas over-expressing this receptor. By several measures, IL13Rα2 expression in patient samples and low-passage primary glioma lines most consistently correlates with the expression of signature genes defining mesenchymal subclass tumors and negatively correlates with proneural signature genes as defined by two studies. Positive associations were also noted with proliferative signature genes, whereas no consistent associations were found with either classical or neural signature genes. Probing the potential functional consequences of this mesenchymal association through IPA analysis suggests that IL13Rα2 expression is associated with activation of proinflammatory and immune pathways characteristic of mesenchymal subclass tumors. In addition, survival analyses indicate that IL13Rα2 over-expression is associated with poor patient prognosis, a single gene correlation ranking IL13Rα2 in the top ~1% of total gene expression probes with regard to survival association with WHO IV gliomas. This study better defines the functional consequences of IL13Rα2 expression by demonstrating association with mesenchymal signature gene expression and poor patient prognosis. It thus highlights the utility of IL13Rα2 as a therapeutic target, and helps define patient populations most likely to respond to immunotherapy in present and future clinical trials.     

The Expression Pattern of microRNAs in Granulosa Cells of Subordinate and Dominant Follicles during the Early Luteal Phase of the Bovine Estrous Cycle

This study aimed to investigate the miRNA expression patterns in granulosa cells of subordinate (SF) and dominant follicle (DF) during the early luteal phase of the bovine estrous cycle. For this, miRNA enriched total RNA isolated from granulosa cells of SF and DF obtained from heifers slaughtered at day 3 and day 7 of the estrous cycle was used for miRNAs deep sequencing. The results revealed that including 17 candidate novel miRNAs, several known miRNAs (n = 291–318) were detected in SF and DF at days 3 and 7 of the estrous cycle of which 244 miRNAs were common to all follicle groups. The let-7 families, bta-miR-10b, bta-miR-26a, bta-miR-99b and bta-miR-27b were among abundantly expressed miRNAs in both SF and DF at both days of the estrous cycle. Further analysis revealed that the expression patterns of 16 miRNAs including bta-miR-449a, bta-miR-449c and bta-miR-222 were differentially expressed between the granulosa cells of SF and DF at day 3 of the estrous cycle. However, at day 7 of the estrous cycle, 108 miRNAs including bta-miR-409a, bta-miR-383 and bta-miR-184 were differentially expressed between the two groups of granulosa cell revealing the presence of distinct miRNA expression profile changes between the two follicular stages at day 7 than day 3 of the estrous cycle. In addition, unlike the SF, marked temporal miRNA expression dynamics was observed in DF groups between day 3 and 7 of the estrous cycle. Target gene prediction and pathway analysis revealed that major signaling associated with follicular development including Wnt signaling, TGF-beta signaling, oocyte meiosis and GnRH signaling were affected by differentially expressed miRNAs. Thus, this study highlights the miRNA expression patterns of granulosa cells in subordinate and dominant follicles that could be associated with follicular recruitment, selection and dominance during the early luteal phase of the bovine estrous cycle.     

Evidence for X-Chromosomal Schizophrenia Associated with microRNA Alterations

Background

Schizophrenia is a severe disabling brain disease affecting about 1% of the population. Individual microRNAs (miRNAs) affect moderate downregulation of gene expression. In addition, components required for miRNA processing and/or function have also been implicated in X-linked mental retardation, neurological and neoplastic diseases, pointing to the wide ranging involvement of miRNAs in disease.

Methods and Findings

To explore the role of miRNAs in schizophrenia, 59 microRNA genes on the X-chromosome were amplified and sequenced in males with (193) and without (191) schizophrenia spectrum disorders to test the hypothesis that ultra-rare mutations in microRNA collectively contribute to the risk of schizophrenia. Here we provide the first association of microRNA gene dysfunction with schizophrenia. Eight ultra-rare variants in the precursor or mature miRNA were identified in eight distinct miRNA genes in 4% of analyzed males with schizophrenia. One ultra-rare variant was identified in a control sample (with a history of depression) (8/193 versus 1/191, p = 0.02 by one-sided Fisher''s exact test, odds ratio = 8.2). These variants were not found in an additional 7,197 control X-chromosomes.

Conclusions

Functional analyses of ectopically expressed copies of the variant miRNA precursors demonstrate loss of function, gain of function or altered expression levels. While confirmation is required, this study suggests that microRNA mutations can contribute to schizophrenia.     

MiRNA Expression Profile of Human Subcutaneous Adipose and during Adipocyte Differentiation

Background

Potential regulators of adipogenesis include microRNAs (miRNAs), small non-coding RNAs that have been recently shown related to adiposity and differentially expressed in fat depots. However, to date no study is available, to our knowledge, regarding miRNAs expression profile during human adipogenesis. Thereby, the aim of this study was to investigate whether miRNA pattern in human fat cells and subcutaneous adipose tissue is associated to obesity and co-morbidities and whether miRNA expression profile in adipocytes is linked to adipogenesis.

Methodology/Principal Findings

We performed a global miRNA expression microarray of 723 human and 76 viral mature miRNAs in human adipocytes during differentiation and in subcutaneous fat samples from non-obese (n = 6) and obese with (n = 9) and without (n = 13) Type-2 Diabetes Mellitus (DM-2) women. Changes in adipogenesis-related miRNAs were then validated by RT-PCR. Fifty of 799 miRNAs (6.2%) significantly differed between fat cells from lean and obese subjects. Seventy miRNAs (8.8%) were highly and significantly up or down-regulated in mature adipocytes as compared to pre-adipocytes. Otherwise, 17 of these 799 miRNAs (2.1%) were correlated with anthropometrical (BMI) and/or metabolic (fasting glucose and/or triglycerides) parameters. We identified 11 miRNAs (1.4%) significantly deregulated in subcutaneous fat from obese subjects with and without DM-2. Interestingly, most of these changes were associated with miRNAs also significantly deregulated during adipocyte differentiation.

Conclusions/Significance

The remarkable inverse miRNA profile revealed for human pre-adipocytes and mature adipocytes hints at a closely crosstalk between miRNAs and adipogenesis. Such candidates may represent biomarkers and therapeutic targets for obesity and obesity-related complications.     

A Seven-microRNA Expression Signature Predicts Survival in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the fifth common cancer. The differential expression of microRNAs (miRNAs) has been associated with the prognosis of various cancers. However, limited information is available regarding genome-wide miRNA expression profiles in HCC to generate a tumor-specific miRNA signature of prognostic values. In this study, the miRNA profiles in 327 HCC patients, including 327 tumor and 43 adjacent non-tumor tissues, from The Cancer Genome Atlas (TCGA) Liver hepatocellular carcinoma (LIHC) were analyzed. The associations of the differentially expressed miRNAs with patient survival and other clinical characteristics were examined with t-test and Cox proportional regression model. Finally, a tumor-specific miRNA signature was generated and examined with Kaplan–Meier survival, univariate\multivariate Cox regression analyses and KEGG pathway analysis. Results showed that a total of 207 miRNAs were found differentially expressed between tumor and adjacent non-tumor HCC tissues. 78 of them were also discriminatively expressed with gender, race, tumor grade and AJCC tumor stage. Seven miRNAs were significantly associated with survival (P value <0.001). Among the seven significant miRNAs, six (hsa-mir-326, hsa-mir-3677, hsa-mir-511-1, hsa-mir-511-2, hsa-mir-9-1, and hsa-mir-9-2) were negatively associated with overall survival (OS), while the remaining one (hsa-mir-30d) was positively correlated. A tumor-specific 7-miRNAs signature was generated and validated as an independent prognostic predictor. Collectively, we have identified and validated an independent prognostic model based on the expression of seven miRNAs, which can be used to assess patients’ survival. Additional work is needed to translate our model into clinical practice.     

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.     

Pulmonary microRNA profiling in a mouse model of ventilator-induced lung injury

The aim of this study was to investigate the changes induced by high tidal volume ventilation (HVTV) in pulmonary expression of micro-RNAs (miRNAs) and identify potential target genes and corresponding miRNA-gene networks. Using a real-time RT-PCR-based array in RNA samples from lungs of mice subjected to HVTV for 1 or 4 h and control mice, we identified 65 miRNAs whose expression changed more than twofold upon HVTV. An inflammatory and a TGF-β-signaling miRNA-gene network were identified by in silico pathway analysis being at highest statistical significance (P = 10(-43) and P = 10(-28), respectively). In the inflammatory network, IL-6 and SOCS-1, regulated by miRNAs let-7 and miR-155, respectively, appeared as central nodes. In TGF-β-signaling network, SMAD-4, regulated by miR-146, appeared as a central node. The contribution of miRNAs to the development of lung injury was evaluated in mice subjected to HVTV treated with a precursor or antagonist of miR-21, a miRNA highly upregulated by HVTV. Lung compliance was preserved only in mice treated with anti-miR-21 but not in mice treated with pre-miR-21 or negative-control miRNA. Both alveolar-arterial oxygen difference and protein levels in bronchoalveolar lavage were lower in mice treated with anti-miR-21 than in mice treated with pre-miR-21 or negative-control miRNA (D(A-a): 66 ± 27 vs. 131 ± 22, 144 ± 10 mmHg, respectively, P < 0.001; protein concentration: 1.1 ± 0.2 vs. 2.3 ± 1, 2.1 ± 0.4 mg/ml, respectively, P < 0.01). Our results show that HVTV induces changes in miRNA expression in mouse lungs. Modulation of miRNA expression can affect the development of HVTV-induced lung injury.