- Letter to the Editor
Serum microRNA profiling in patients with glioblastoma: a survival analysis
Molecular Cancervolume 16, Article number: 59 (2017)
Because circulating microRNAs (miRNAs) have drawn a great deal of attention as promising novel cancer diagnostics and prognostic biomarkers, we sought to identify serum miRNAs significantly associated with outcome in glioblastoma patients. To do this, we performed global miRNA profiling in serum samples from 106 primary glioblastoma patients. The study subjects were randomly divided into two sets: set one (n = 40) and set two (n = 66). Using a Cox regression model, 3 serum miRNAs (miR-106a-5p, miR-182, and miR-145-5p) and 5 serum miRNAs (miR-222-3p, miR-182, miR-20a-5p, miR-106a-5p, and miR-145-5p) were identified significantly associated with 2-year patient overall survival and disease-free survival (P < 0.05) in both sets and the combined set. We then created the miRNA risk scores to assess the total impact of the significant serum miRNAs on survival. The high risk scores were associated with poor patient survival (overall survival: HR = 1.92, 95% CI: 1.19, 10.23, and disease-free survival: HR = 2.03, 95%CI: 1.24, 4.28), and were independent of other clinicopathological factors. Our results suggest that serum miRNAs could serve as prognostic predictors of glioblastoma.
microRNAs (miRNAs) are small endogenous mediators of RNA interference and key regulatory components of many biological processes required for an organism’s development, cell specialization, and homeostasis [1, 2]. Many miRNAs exhibit tissue-specific patterns of expression and are deregulated in various cancers, where they can be either oncogenic (oncomirs) or tumor suppressive. In glioblastoma, the most common malignant brain tumor in adults [3–5], miRNA dysregulation plays diverse roles in various aspects of gliomagenesis, ranging from cancer stem cell regulation (miRs-7, 9, 10a/10b, 17–92, 124/137, 125a/125b, 182-5p, and 326), cell cycle control (miRs-15b, 34a, and 221/222), apoptosis (miRs-21, 26a, 101, 146b-5p, 153, 181a/181b, 196a/196b, 218, 381, 451), and angiogenesis (miRs-93 and 296), to immune modulation (miR-124, 138, 142-3p) [6–12].
miRNAs from tumors can be secreted within membrane vesicles (exosomes) [13, 14] or directly into the blood [15, 16], indicating that miRNAs probably play a key role in intracellular communication. As such, miRNAs are increasingly being profiled as biomarkers for diseases, including cancers [17–23], circulating in the bloodstream. Profiling of miRNAs circulating in serum or plasma may be particularly useful in glioblastoma patients, given the risks of sequential surgery and biopsy, and the lack of readily usable biomarkers to predict clinical outcome. Previous attempts have been made to assess circulating miRNAs in relation to glioblastoma prognosis [24–30]. However, due to small sample size and a heterogeneous patient population, prognostic impact could not be elucidated. In the current study, we perform miRNA profiling in serum samples from 106 primary glioblastoma patients from The University of Texas M. D. Anderson Cancer Center and identify a panel of miRNAs associated with patient outcome.
Study subjects were recruited from The University of Texas M. D. Anderson Cancer Center (Houston, TX) beginning in April 2013 under an Institutional Review Board approved protocol. Written informed consent was obtained from each study subject. Demographic and clinical data for the study subjects were obtained from medical record review. Serum samples were selected retrospectively at the time of analysis according to the following requirements: 1) the patients were newly diagnosed with glioblastoma defined by WHO 2007 criteria ; 2) serum was available in sufficient volume for RNA isolation; and 3) demographic, clinical, and follow-up data for the patients were available. All patients received initial standard of care treatment of resection, radiation, and chemotherapy. Case patients with recurrent glioblastoma or prior cancer (except for non-melanoma skin cancer) were excluded from the study. A total of 106 glioblastoma patients were included in the study.
Serum miRNAs were isolated using the miRNeasy kit (Qiagen), with minor modifications. miRNA expression profiling was carried out using the nCounter Human v2 miRNA Expression Assay (NanoString Technologies). For each patient, disease-free survival was calculated from the date of disease diagnosis to the date of disease recurrence/progression, or to date of death from disease, or to the last follow-up date. Overall survival was calculated from the date of disease diagnosis to the date of death from disease or to the date of last patient follow-up visit. The associations between 2-year survival and serum miRNA expression levels were estimated as hazard ratios (HRs) and 95% confidence intervals (CIs) for each miRNA using the Cox proportional hazards model, adjusted for age, sex, race, smoking status, Karnofsky Performance Scale (KPS) score, timing of blood draw, IDH1 mutation status, and body mass index (BMI). Patients with survival times over 24 months were censored at 24 months in the Cox regression analysis. The proportional hazards assumption for identified significant miRNAs was tested within each set using the Schoenfeld residuals. miRNAs with evidence for non-proportionality were deleted. The combined 3-miRNA risk score (for overall survival) and 5-miRNA risk score (for disease-free survival) for each patient was derived by linear combination of the product of the reference-normalized expression level of each miRNA times its Cox regression corresponding coefficient. All patients were dichotomized by the median risk score, and individuals with a risk score higher or lower than the median were classified as high- or low-risk groups, respectively. In all statistical analyses, a P value < 0.05 was considered significant. Then, an estimate of the false discovery rate (q value) was calculated to take into account the multiple comparisons for any miRNA significantly associated with 2-year overall survival or disease-free survival in all two sets. A low q value (q < 0.15) is an indication of high confidence in a result. miRNAs that were significantly in both two sets were deemed as the significant miRNAs.
Serum miRNAs predicting overall survival
Basic characteristics of the patient cohort at baseline are shown in Table 1. A total of 106 patients with primary glioblastoma were included in the study. We divided the study subjects into two sets: set one and set two. The set one included 40 glioblastoma patients, 20 who had survived at least 24 months and 20 who died within 24 months after diagnosis. The set two included 66 glioblastoma patients, 53 who had survived at least 24 months and 13 who died within 24 months after diagnosis. Overall, basic characteristics were very similar between two sets. Overall, the mean age was 58 years, and 66% of the study subjects were male. Most of them were Caucasians (94.3%), and nearly half of them had a history of cigarette smoking. Two-thirds of the study subjects had a KPS score of at least 90. Of 67 study subjects for whom IDH1 mutation data was available, 9 (13.4%) carried an IDH1 mutation. In terms of the timing of obtaining blood samples, nearly half of the study subjects (46.2%) had their blood drawn immediately after surgery; in 24.5%, blood was drawn after or during adjuvant chemotherapy; and in 26.4%, blood was drawn after or during chemoradiation. The median follow-up interval was 19 months. During the follow-up period, 68 study subjects (64.8%) experienced disease recurrence or progression.
Of 812 human miRNAs we were capable of detecting, 187 miRNAs were expressed in at least 80% of the patients’ serum samples. After imputation with minimum observed values for each serum miRNA, we first attempted to identify serum miRNAs associated with 2-year overall survival times (Table 2). In set one, 5 serum miRNAs (miR-630, 106a-5p, 182, 15b-5p, and 145-5p) were associated with 2-year overall survival. Among them, 2 serum miRNAs (miR-630, 106a-5p) were significantly associated with a lower probability of 2-year overall survival, and 3 serum miRNAs were significantly associated with a higher probability of 2-year overall survival (miR-182, 15b-5p and 145-5p). In set two, 5 serum miRNAs (miR-106a-5p, 20a-5p, 17-5p, 182, and 145-5p) were associated with 2-year overall survival. Among them, 3 serum miRNAs (miR-106a-5p, 20a-5p, 17-5p) were significantly associated with a lower probability of 2-year overall survival, and 2 serum miRNAs (miR-182, 145-5p) were significantly associated with a higher probability of 2-year overall survival. Three serum miRNAs were identified as significantly associated (P <0.05) with overall survival in both sets; these were miR-106a-5p, miR-182, and miR-145-5p. After adjusting for multiple comparisons, the associations were remained significant (q < 0.15). These associations with 2-year overall survival were further confirmed in a meta-analysis by combining the samples from two sets. Increased levels of miR-106a-5p and decreased levels of miR-182 and miR-145-5p were significantly associated with a lower probability of 2-year overall survival. Additionally, in the meta-analysis, we observed 2 more serum miRNAs, miR-20a-5p and miR-17-5p, whose increased levels were significantly associated with an increased hazard of death (P = 0.040 and 0.043, respectively).
Serum miRNAs predicting disease-free survival
Next we investigated the associations between serum miRNAs and 2-year disease-free survival (Table 3). In set one, 7 serum miRNAs (miR-922, 222-3p, 20a-5p, 106.5p, 182, 548a, 145-5p) were associated with 2-year disease-free survival. Among them, 4 serum miRNAs (miR-922, 222-3p, 20a-5p, 106a-5p) were significantly associated with a lower probability of 2-year disease-free survival, and 3 serum miRNAs (miR0182, 548a, 145-5p) were significantly associated with a higher probability of 2-year disease-free survival. In set two, 7 serum miRNAs (miR222-3p, 20a-5p, 106a-5p, 124a, 182, 371-5p, 145-5p) were associated with 2-year disease-free survival. Among them, 3 serum miRNAs (miR-222-3p, 20a-5p, 106a-5p) were significantly associated with a lower probability of 2-year disease-free survival, and 4 serum miRNAs (miR-124a, 182, 371a-5p, 145-5p) were significantly associated with a higher probability of 2-year disease-free survival. Five serum miRNAs (P <0.05) showed a significant association with 2-year disease-free survival in both sets: miR-222-3p, miR-182, miR-20a-5p, miR-106a-5p, and miR-145-5p. After adjusting for multiple comparisons, the associations were remained significant (q < 0.15). Their significant associations with 2-year disease-free survival were further confirmed in a meta-analysis by combining the samples from two sets. Increased levels of miR-222-3p, miR-20a-5p, and miR-106a-5p, and decreased levels of miR-182 and miR-145-5p were significantly associated with a lower probability of 2-year disease-free survival.
miRNA risk scores predicting survival outcomes
Last, we assessed the combined effects of these serum miRNAs on survival. We generated 2 risk scores, one for overall survival and the other for disease-free survival, using the linear combination of expression levels of these significant miRNAs for each individual based on the median levels of each single miRNA (Table 4). The included miRNAs were miR-106a-5p, 182 and 145-5p for overall survival, and miR-222-3p, 20a-5p, 106a-5p, 182 and 145-5p for disease-free survival. For 2-year overall survival, study subjects with high-risk scores exhibited a significantly increased hazard of death relative to those with low-risk scores (HR = 1.92, 95% CI, 1.19, 10.23). The Kaplan-Meier 2-year overall survival curve showed that patients with the high-risk scores had a median survival time (MST) of 15 months, compared with 22 months in the low-risk group (P < 0.001) (Fig. 1). For 2-year disease-free survival, study subjects with high-risk scores exhibited a significantly decreased chance for disease-free survival relative to those with low-risk scores (HR = 2.03, 95% CI, 1.24–4.28). The Kaplan-Meier 2-year disease-free survival curve showed that patients with the high-risk scores had a median survival time (MST) of 9.0 months, compared with 18 months in the low-risk group (P < 0.001) (Fig. 2).
Although the idea of using circulating biomarkers, either alone or as an adjunct to genetic and molecular profiling of the tumor, to monitor tumor progression/recurrence and clinical outcome prediction has already entered clinical practice [32, 33], little progress has been made in this context for glioblastoma patients. In the current study, we performed global miRNA profiling in serum from 106 patients with primary glioblastoma and identified a panel of serum miRNAs associated with glioblastoma prognosis.
In recent years, endogenous circulating miRNAs have attracted extensive research interest because they remain stable and are well protected even under harsh conditions. A number of studies have investigated the potential role of circulating miRNAs in glioma survival and clinical outcomes [24, 25, 27, 29, 30]. For example, Zhi et al. reported that up regulation of serum miR-20a-5p, miR-106a-5p, and miR-181b-5p was associated with advanced clinical stages of astrocytoma, and high expression of serum miR-19a-3p, miR-106a-5p, and miR-181b-5p was significantly associated with poor patient survival among astrocytoma patients . In our study, we found that elevated levels of serum miR-106a-5p and miR-20a-5p were associated with poor survival. Thus, our results support their findings. Both miR-20a and miR-106a have previously been linked to gliomagenesis. Overexpression of miR-20a and miR-106a was observed in glioma stem cells, and such overexpression could promote glioma invasiveness by directly targeting tissue inhibitor of metalloproteinases-2 (TIMP-2) . Knockdown of either miR-20a or miR-106a was observed to increase TIMP-2 expression and consequently inhibit glioma stem cell invasion. Xiao et al. found that high levels of plasma miR-182 were associated with poor overall survival (P = 0.034) and disease-free survival (P = 0.013) among glioma patients . However, in our study, we found that decreased levels of miR-182 were associated with poor overall survival (P = 0.037) and disease-free survival (P = 0.043) among glioblastoma patients. The reasons for this discrepancy are unclear, but it may be due differences in the study populations. Xiao’s study included patients with a mixture of glioma types ranging from WHO grade I to IV, whereas our study was confined to glioblastoma patients. Using oncogenomic analyses of 470 miRNAs profiled by The Cancer Genome Atlas (TCGA) program, Kouri et al. revealed that miR-182 was the only miRNA that was associated with favorable patient prognosis and temozolomide (TMZ) sensitivity, and was associated with a less aggressive oligoneural subclass of glioblastoma . In another study from the same group, miR-182 was found to play a key role in apoptosis, growth, and differentiation in glioblastoma , and overexpression of miR-182 in orthotopic glioblastoma xenografts led to reduced tumor burden and increased animal survival.
In addition to the miRNAs mentioned above, we observed that elevated serum levels of miR-222-3p and miR-17-5p, and decreased levels of miR-145-5p, were associated with a poor glioblastoma prognosis. Although similar associations have not been reported in the literature yet, these miRNAs have been linked to various phenotypes of glioblastoma. For example, overexpression of miR-222 leads to increased DNA damage, tumor growth promotion, and a poor glioblastoma prognosis . On the other hand, suppression of miR-222 expression inhibits glioma cell proliferation and invasion. Elevated levels of miR-17 were found in glioblastoma samples and were negatively related to patients’ survival . The level of miR-17 is also increased in glioblastoma spheroids, which are enriched in tumor-initiating cells (TICs) or tumor stem-like cells (TSCs) . In contrast, miR-145 is an important repressor of pluripotency in embryonic stem cells and acts as a tumor suppressor in a variety of cancers, including glioma . Lee et al. found that the low expression of miR-145 in glioblastoma was correlated with poor patient prognosis . Thus, our findings on serum miR-222-3p, miR-17-5p, and miR-145-5p are consistent with their biological functions in glioblastoma carcinogenesis.
In summary, we identified a panel of serum miRNAs associated with 2-year overall survival and disease-free survival among patients with glioblastoma. Our study provides evidence to support the role of circulating miRNAs in predicting survival with glioblastoma. Further validation from larger independent studies and further characterization of these candidate serum miRNAs are warranted.
Diederichs S, et al. The dark matter of the cancer genome: aberrations in regulatory elements, untranslated regions, splice sites, non-coding RNA and synonymous mutations. EMBO Mol Med. 2016;8:442–57.
Tomaru Y, et al. Cancer research with non-coding RNA. Cancer Sci. 2006;97:1285–90.
Chan JA, et al. MicroRNA-21 is an antiapoptotic factor in human glioblastoma cells. Cancer Res. 2005;65:6029–33.
Sasayama T, et al. MicroRNA-10b is overexpressed in malignant glioma and associated with tumor invasive factors, uPAR and RhoC. Int J Cancer. 2009;125:1407–13.
Gabriely G, et al. Human glioma growth is controlled by microRNA-10b. Cancer Res. 2011;71:3563–72.
Karsy M, et al. Current Progress on Understanding MicroRNAs in Glioblastoma Multiforme. Genes Cancer. 2012;3:3–15.
Touat M, et al. Emerging circulating biomarkers in glioblastoma: promises and challenges. Expert Rev Mol Diagn. 2015;15:1311–23.
Shea A, et al. MicroRNAs in glioblastoma multiforme pathogenesis and therapeutics. Cancer Med. 2016;5:1917–46.
Xue J, et al. miR-182-5p Induced by STAT3 Activation Promotes Glioma Tumorigenesis. Cancer Res. 2016;76:4293–304.
Wei J, et al. MiR-138 exerts anti-glioma efficacy by targeting immune checkpoints. Neuro Oncol. 2016;18:639–48.
Xu S, et al. Effect of miR-142-3p on the M2 macrophage and therapeutic efficacy against murine glioblastoma. J Natl Cancer Inst. 2014;106.
Wei J, et al. miR-124 inhibits STAT3 signaling to enhance T cell-mediated immune clearance of glioma. Cancer Res. 2013;73:3913–26.
van der Vos KE, et al. Directly visualized glioblastoma-derived extracellular vesicles transfer RNA to microglia/macrophages in the brain. Neuro Oncol. 2016;18:58–69.
Xu S, et al. Exosome-mediated microRNA transfer plays a role in radiation-induced bystander effect. RNA Biol. 2015;12:1355–63.
Fernandez-Mercado M, et al. The circulating transcriptome as a source of non-invasive cancer biomarkers: concepts and controversies of non-coding and coding RNA in body fluids. J Cell Mol Med. 2015;19:2307–23.
Igaz I, et al. Diagnostic Relevance of microRNAs in Other Body Fluids Including Urine, Feces, and Saliva. EXS. 2015;106:245–52.
Wang J, et al. Tumor-associated circulating microRNAs as biomarkers of cancer. Molecules. 2014;19:1912–38.
Schwarzenbach H, et al. Clinical relevance of circulating cell-free microRNAs in cancer. Nat Rev Clin Oncol. 2014;11:145–56.
Healy NA, et al. Systemic mirnas as potential biomarkers for malignancy. Int J Cancer. 2012;131:2215–22.
Mittra I, et al. Nucleic acids in circulation: are they harmful to the host? J Biosci. 2012;37:301–12.
Ma R, et al. Circulating microRNAs in cancer: origin, function and application. J Exp Clin Cancer Res. 2012;31:38.
Anindo MI, et al. Insights into the potential use of microRNAs as biomarker in cancer. Int J Surg. 2012;10:443–9.
Jarry J, et al. The validity of circulating microRNAs in oncology: Five years of challenges and contradictions. Mol Oncol. 2014;8:819–29.
Xiao Y, et al. Potential Diagnostic and Prognostic Value of Plasma Circulating MicroRNA-182 in Human Glioma. Med Sci Monit. 2016;22:855–62.
Herman A, et al. Analysis of Glioblastoma Patients’ Plasma Revealed the Presence of MicroRNAs with a Prognostic Impact on Survival and Those of Viral Origin. PLoS ONE. 2015;10, e0125791.
Dong L, et al. miRNA microarray reveals specific expression in the peripheral blood of glioblastoma patients. Int J Oncol. 2014;45:746–56.
Yang C, et al. Identification of seven serum microRNAs from a genome-wide serum microRNA expression profile as potential noninvasive biomarkers for malignant astrocytomas. Int J Cancer. 2013;132:116–27.
Ilhan-Mutlu A, et al. Plasma MicroRNA-21 concentration may be a useful biomarker in glioblastoma patients. Cancer Invest. 2012;30:615–21.
Zhi F, et al. Identification of 9 serum microRNAs as potential noninvasive biomarkers of human astrocytoma. Neuro Oncol. 2015;17:383–91.
Yue X, et al. Downregulation of serum microRNA-205 as a potential diagnostic and prognostic biomarker for human glioma. J Neurosurg. 2016;124:122–8.
Louis DN, et al. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol. 2007;114:97–109.
Crowley E, et al. Liquid biopsy: monitoring cancer-genetics in the blood. Nat Rev Clin Oncol. 2013;10:472–84.
Diaz Jr LA, et al. Liquid biopsies: genotyping circulating tumor DNA. J Clin Oncol. 2014;32:579–86.
Wang Z, et al. Oncogenic miR-20a and miR-106a enhance the invasiveness of human glioma stem cells by directly targeting TIMP-2. Oncogene. 2015;34:1407–19.
Kouri FM, et al. miRNA-182 and the regulation of the glioblastoma phenotype - toward miRNA-based precision therapeutics. Cell Cycle. 2015;14:3794–800.
Kouri FM, et al. miR-182 integrates apoptosis, growth, and differentiation programs in glioblastoma. Genes Dev. 2015;29:732–45.
Quintavalle C, et al. MiR-221/222 target the DNA methyltransferase MGMT in glioma cells. PLoS ONE. 2013;8, e74466.
Malzkorn B, et al. Identification and functional characterization of microRNAs involved in the malignant progression of gliomas. Brain Pathol. 2010;20:539–50.
Ernst A, et al. De-repression of CTGF via the miR-17-92 cluster upon differentiation of human glioblastoma spheroid cultures. Oncogene. 2010;29:3411–22.
Speranza MC, et al. NEDD9, a novel target of miR-145, increases the invasiveness of glioblastoma. Oncotarget. 2012;3:723–34.
Lee HK, et al. MicroRNA-145 is downregulated in glial tumors and regulates glioma cell migration by targeting connective tissue growth factor. PLoS ONE. 2013;8, e54652.
Availability of data and materials
Please contact the corresponding author for all data requests.
Conception and design: HZ and ABH; Development of methodology: HZ and JS; Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): JS, TRH and RS; Analysis and interpretation of data: HZ, JS, and ABH; Writing, review, and/or revision of the manuscript: HZ, JS, TRH, GNF, and ABH; Administrative, technical, or material support: None. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Study subjects were recruited from The University of Texas M. D. Anderson Cancer Center (Houston, TX) beginning in April 2013 under an Institutional Review Board approved protocol. Written informed consent was obtained from each study subject. Demographic and clinical data for the study subjects were obtained from medical record review.