- Open Access
MiR-1248: a new prognostic biomarker able to identify supratentorial hemispheric pediatric low-grade gliomas patients associated with progression
Biomarker Research volume 10, Article number: 44 (2022)
Pediatric low-grade gliomas (pLGGs), particularly incompletely resected supratentorial tumours, can undergo progression after surgery. However to date, there are no predictive biomarkers for progression. Here, we aimed to identify pLGG-specific microRNA signatures and evaluate their value as a prognostic tool.
We identified and validated supratentorial incompletey resected pLGG-specific microRNAs in independent cohorts from four European Pediatric Neuro-Oncology Centres.
These microRNAs demonstrated high accuracy in differentiating patients with or without progression. Specifically, incompletely resected supratentorial pLGGs with disease progression showed significantly higher miR-1248 combined with lower miR-376a-3p and miR-888-5p levels than tumours without progression. A significant (p < 0.001) prognostic performance for miR-1248 was reported with an area under the curve (AUC) of 1.00. We also highlighted a critical oncogenic role for miR-1248 in gliomas tumours. Indeed, high miR-1248 levels maintain low its validated target genes (CDKN1A (p21)/FRK/SPOP/VHL/MTAP) and consequently sustain the activation of oncogenic pathways.
Altogether, we provide a novel molecular biomarker able to successfully identify pLGG patients associated with disease progression that could support the clinicians in the decision-making strategy, advancing personalized medicine.
Pediatric low-grade gliomas (pLGGs) are a histologically heterogeneous group of tumours which include pilocytic astrocytomas (PA), the most frequent histology, and non-PA tumours, that comprise angiocentric gliomas (AG), dysembryoplastic neuroepithelial tumours (DNET), gangliogliomas (GG), and glioneuronal tumours (GNT). These tumours are classified as grade 1 and 2 tumours by the WHO and account for 40% of brain tumours under the age of 18 years [1,2,3,4].
Safely resectable pLGGs, primarily occurring in the cerebellum and in non-eloquent supratentorial locations, can often be cured by surgery with an excellent long-term prognosis . Deep supratentorial tumours and those arising in the brainstem are more difficult to resect. Of note, disease progression occurs in about 30% of the cases, and these patients receive adjuvant chemotherapy and, in selected cases, radiotherapy (e.g. SIOP‐LGG 2004 protocol and LOGGIC/LOGGIC Core) [6, 7].
The cost of survivorship is often high, since patients frequently present long-term sequelae both by the disease itself and as a consequence of the treatment, such as epilepsy, vision loss, impaired motor skills and cognitive dysfunction .
The genomic landscape of pLGGs has allowed the revision of the WHO classification resulting in the new 2021 WHO classification of the tumours of the Central Nervous System (CNS) . Moreover, recent studies have stratified patients into different risk classes on the basis of selected genetic alterations  and on histological and clinical variables .
Research groups, including ours, have investigated microRNA profiles in pLGGs. Profiling studies on PAs demonstrated the over-expression of miR-21 and miR-34 and the downregulation of miR-124 and miR-129 compared to normal brain tissue [11,12,13,14]. Furthermore, miR-487b was reported to be downregulated in these tumours [12, 15, 16]. Low levels of miR-10b-5p were reported in sporadic and NF1-associated PAs . In addition, miR-20a-5p was identified as a regulator of both the MAPK/ERK and mTORC1 pathways, in a subset of genetically driven astrocytomas , while miR-125b was under-expressed in pleomorphic xanthoastrocytomas and GG . Finally, miR-519d and miR-4758 were differentially expressed in GG compared to other histological subtypes .
In this context, we previously investigated a cohort of supratentorial PA and non-PA tumours, showing that all of them are characterized by low expression of miR-139-5p and that this deregulation contributes to tumour growth by sustaining the pro-tumoral PI3K/AKT signalling pathway .
Two studies investigated microRNAs as possible prognostic biomarkers in pLGGs. Low levels of miR-29b-3p together with a high level of its target gene, cannabinoid receptor 1 (CNR1), were suggested to be predictive of residual disease involution in a cohort of subtotally resected pLGGs . In the second study, a 7-microRNAs signature was associated with chemotherapy response .
The identification of high risk of progression in incompletely resected pLGGs remains to date an unmet clinical need. Therefore, the aim of this study was to analyse microRNA profiles in pLGG tissue samples in a large cohort of incompletely resected patients to identify potential biomarkers of progression. To this end, samples were collected from independent cohorts of patients from four European Pediatric Neuro-Oncology Centres and two different technologies were used to determine microRNAs able to differentiate patients with and without disease progression. Finally, we validated miR-1248 as a prognostic biomarker for tumours with progression.
Materials and methods
Characteristics of pLGG samples
This study was performed on 43 supratentorial incompletely resected pLGG samples (PA and non-PA) collected from four European pediatric neuro-oncology centres. The workflow of the study is depicted in Supplementary Figure S1.
The first cohort (Cohort I) derived from three independent European institutions (Bambino Gesù Children’s Hospital, A. Gemelli Hospital and Marseille Hospital) and consisted of 20 PA and 9 non-PA tumour samples. The second cohort (Cohort II) derived from one European institution, the DKFZ German Research Centre and consisted of 9 PA tumour samples.
According to the extent of surgical resection on Magnetic Resonance Imaging (MRI) after one month post-surgery or post biopsy, patients were subdivided in those with near total resection (NTR), where an amount between 50–90% of the mass was removed, and those with biopsy, where less than 50% of the tumour was removed.
Patients with any residual tumour (NTR and biopsy) were further subdivided in with or without progression according to disease progression, on the basis of MRI evaluation during the follow up and/or clinical recurrence identified by clinicians . Kaplan Meier analysis was performed using the progression free survival (PFS) data reported in Tables 1 and 2, and the results demonstrated a lower PFS for pLGG patients with progression compared to patients without progression (Supplementary Figure S2).
Histological diagnosis was performed according to WHO 2021 classification criteria  by experienced neuropathologists (FG, FDC, DFB). All patients before surgery were naïve for chemo and/or radiotherapy. Ethical approval (Rif. 5866) was obtained in accordance with the Helsinki declaration of 1964 and its later amendments. Informed written consent was obtained from the patients, parents or guardians before enrolment, according to our ethical committee guidelines.
Samples used for microRNA profiling
MicroRNA expression levels were obtained from 29 supratentorial Grade I PA and non-PA tumours belonging to Cohort I by the use of Real-Time quantitative PCR (RT-qPCR) TaqMan Low Density Array (TLDA) microfluidic cards (Human miR v3.0, Life Technologies). Clinical and pathological data of Cohort I are reported in Table 1. Mean age was 7.6 years, with a median of 7. The most frequent histology was PA (n = 20, 59%), followed by GG (n = 8, 23%), DNET (n = 4, 12%), AG (n = 1, 3%) and glioneuronal tumour (GNT, n = 1, 3%). Sixteen PAs (80%) were located in midline structures: six in optic pathways, five in the third ventricle and four in the thalamus. Only four (20%) arose in hemispheric structures. DNET and AG were hemispheric, while the GNT arose in the third ventricle. Five GG arose in the temporal lobe, while three in midline structures, specifically one in the diencephalon, one in optic pathways and one in thalamus. (62%, DNET n = 2; GG n = 6; PA n = 13) experienced tumour progression with a PFS average time of 4.5 years and tumours were equally distributed between hemispheric (n = 10) and midline (n = 11) areas. Thirteen patients (38%, AG n = 1; DNET n = 2; GG n = 2; GNT n = 1; PA n = 6) did not progress, with tumours mostly located in midline structures (9 out of 13, 69%). Of note, DNA methylation analysis was performed in 18 samples (Supplementary Table S1). On the basis of the new classifier version (v12.3) [24, 25] high score (≥ 0.84) was obtained in 12/18 (67%) patients and histological diagnosis was confirmed in 9/12 (75%) cases. Conversely, three tumours with an original pathological diagnosis of GG were classified as PA. 6/18 (33%) tumours were attributed a low methylation score with consequent inconsistent methylation classification.
MicroRNA-sequencing was performed in the 9 samples of Cohort II at the DKFZ German Research Centre to obtain microRNA expression profiles. Clinical and pathological features of this cohort are reported in Table 2. All patients underwent NTR and the mean age at diagnosis was 4.1, while the median age was 5 years old. All NTR supratentorial tumours involved midline structures: six were located in the optic pathways (60%), two in the thalamus (20%) and two in the third ventricle (20%). In the NTR subgroup, three patients (33.3%) presented progression with a mean PFS of 8.6 months.
Samples used for validation of profiling results
Validation of profiling results was performed in a series of Cohort I samples (27 out of the 29 that were used for profiling plus 5 new samples) through droplet digital PCR (ddPCR). Clinical and pathological features of samples are reported in Table 1, indicating which samples were used for profiling and validation of the results.
Genomic landscape of samples
Cohort I samples were screened for the KIAA1549:BRAF fusion gene variants (KIAA1549:BRAF exon16-exon9 (K16B9) (n = 33) and KIAA1549:BRAF exon16-exon11 (K16B11) (n = 33)) using the Applied Biosystems ViiA 7 RT-qPCR system and validated by PCR-based Sanger sequencing. The latter was also used for the detection of KIAA1549:BRAF exon15-exon9 (K15B9) fusion (n = 19). The BRAF V600E point mutation (n = 33) was analysed by RT-qPCR using TaqMan probes, as previously described  (Table 1). DNA methylation of Cohort I samples (n = 18) was also evaluated by DNA Methylation array . Results were analysed both with the old (v11b4) and new (v12.3) classifier versions (Supplementary Table S1). Samples with calibrated scores higher than 0.84 were considered informative, while low scores suggesting classifier prediction uncertainty were not taken into consideration . In addition to microRNA expression levels both RNA sequencing and whole genome sequencing data were available for seven samples of Cohort II (Table 2) .
For Cohort I, microRNA expression profiling was performed on fresh frozen (FF) (n = 19) or FFPE (n = 9) tissues using TLDA microfluidic cards (Human miR v3.0, Life Technologies), that analyse the 754 most studied microRNAs, as previously reported . For Cohort II samples microRNA sequencing, that allows the assessment of all known microRNAs (2656 mature microRNAs in the latest release of miRBase; http://www.mirbase.org/) , was performed in the DKFZ German Research Centre Heidelberg, as previously described .
Droplet Digital PCR (ddPCR)
Ten nanograms of RNA were retrotranscribed using TaqMan™ MicroRNA Assay (Life Technologies) according to manufacturer’s instructions. The resulting cDNA was diluted 1:6 and 8 μl were used to prepare a 22 μl reaction mix containing 11 μl of 2X ddPCR Supermix for Probes (Bio-Rad) and 1,1 μl 20X TaqMan miRNA PCR primer probe set (Life Technologies). The PCR mixes for each sample were loaded in a disposable cartridge (Bio-Rad) together with 70 μl of droplet generation Oil (Bio-Rad) and loaded in the QX200 droplet generator (Bio-Rad).
40 μl of droplets were then transferred in a 96 well plate and an endpoint PCR was performed using the following conditions: 95 °C for 10 min, then 45 cycles of 95 °C for 15 s and 58 °C for 1 min, and a final step at 98 °C for 10 min. Then, the 96 well plate was placed in the QX200 Droplet Reader for detection of positive droplets. The quantification of positive droplets was performed using the QuantaSoft software (Bio-Rad).
MicroRNA expression profiling data from Cohort I were processed with Statminer™ Software v 5.0 (Integromics TM) and differential expression analysis between pLGG samples with and without progression was performed, as previously described . Briefly, microRNA expression was normalized using the global normalization and differential expression was evaluated using the limma test. P values < 0.05 were considered as statistically significant.
MicroRNA sequencing data of Cohort II were analysed as described in . Samples were classified based on the presence or absence of progression, and differential expression analysis was performed using the R package DESeq2. P values < 0.05 were considered statistically significant.
Clustering and heatmaps were also generated in R using the pheatmap function. The different range of microRNA expression values in the two cohorts (Cycle Threshold and Reads per Million mapped reads values), is represented by the diverse heatmap colours. Kaplan Meier analyses were performed using PFS and OS data of all patients from cohorts I and II using GraphPad Prism software Version 8 (La Jolla, CA).
For ddPCR data, statistical analyses were performed using GraphPad Prism software Version 8 (La Jolla, CA). The t-test for unpaired data was used to analyse differences in microRNA expression between experimental groups. P values < 0.05 were considered statistically significant.
Receiver operating characteristic (ROC) curves for the microRNAs of interest were calculated using GraphPad Prism version 8 (La Jolla, California, USA).
Univariate analysis was performed in IBM SPSS Statistics version 27 (Armonk, New York, USA), using the General Linear Model (GLM) for miR-1248. In detail, the Progression Free Survival (PFS) was used as a covariate and other clinical features (site of onset, histology, BRAF status, age, sex and type of resection) as fixed factors. Univariate analysis was performed to determine the significance of the clinical features on the expression of the microRNA. The Parameter Estimates summarized the effect of each clinical factor. P < 0.05 was considered statistically significant.
MicroRNA validated target identification was performed for miR-1248 downloading the validated targets genes for Homo sapiens Release 8.0 from miRTarBase . Tumour suppressor genes were downloaded from TSGene 2.0 web tool , that reports 983 downregulated tumour suppressor genes in The Cancer Genome Atlas (TCGA) pan-cancer samples compared to normal samples. Intersection of validated target genes of miR-1248 with 983 tumour suppressor genes was performed.
Cohort I microRNA array profiling
MicroRNA profiling was performed on supratentorial PA and non-PA samples from Cohort I (Table 1) divided in tumours with or without progression (see materials and methods for details). We found a higher number of detectable microRNAs in the subgroup of patients without progression (401/754, 53%) in comparison to those with progression (331/754, 44%). Differential expression analysis resulted in 38 microRNAs, nine upregulated and 29 downregulated in tumours with progression (Supplementary Table S2). Hierarchical clustering distinguished two subgroups: one branch included predominantly tumours with progression and the second branch mixed samples (Fig. 1).
Of note, samples with progression were equally distributed between hemispheric (n = 10) and midline (n = 10) areas, whereas those who did not progress were all located in midline structures (n = 9). This result underlines that differentially expressed (DE) microRNAs reflect at least in part the distinct embryologic origin between tumours, hemispheric or midline, respectively.
Cohort II microRNA sequencing
As for Cohort I, a second independent cohort included supratentorial PA samples subdivided into tumours with or without progression (Table 2). The microRNA differential expression analysis resulted in 32 upregulated and 52 downregulated microRNAs in tumours with progression (Supplementary Tables S3). The DE microRNAs were used for hierarchical clustering analysis, which showed distinct branches for tumours that progressed from the ones that did not (Fig. 2).
All features of the samples included in the study are reported in Fig. 2 and Table 2. Samples with progression were located in the thalamus (n = 2) and in the optic pathway (n = 1), while tumours that did not progress arose in the optic pathway (n = 4) and in the third ventricle (n = 2).
Deregulated microRNAs in the two cohorts
Since our aim was to identify microRNAs able to stratify incompletely resected pLGG patients into progression risk categories, we further analysed the results obtained from Cohort I and II focusing on the deregulated microRNAs in both cohorts. Specifically, analysing data generated from two independent cohorts with two different technologies and considering all samples from both cohorts, we were able to identify three deregulated microRNAs. In detail, miR-376a-3p and miR-888-5p were downregulated, while miR-1248 was upregulated in pLGGs, both PA and non-PA, with progression (Fig. 3 and Supplementary Figure S3).
Evaluation of miR-376a-3p, miR-888-5p and miR-1248 as prognostic biomarkers
With the aim to validate the profiling results, miR-376a-3p, miR-888-5p and miR-1248 expression levels were evaluated using an absolute quantification method. Specifically, the three microRNAs were analysed in 27 out of 29 samples previously used for profiling plus 5 new cases for a total of 32 samples (see Table 1).
DdPCR allowed to confirm the downregulation of miR-376a-3p and upregulation of miR-1248 in both PA and non-PA tumours with progression, whereas miR-888-5p was not statistically significant (Fig. 4a).
ROC analysis was performed to test the specificity and sensitivity of expression of these microRNAs. MiR-376a-3p (Supplementary Figure S4) was not statistically significant, while miR-1248 resulted in an AUC = 0.815 (95% CI: 0.6682 to 0.9634) and a P value of 0.0028 (optimal miR-1248 cut-off value for pLGGs with progression > 0.31 copies/μL) (Fig. 4b). This result highlights that miR-1248 is able to distinguish both PA and non-PA with progression respect to those without.
Next, we proceeded with in-silico analysis in order to identify potential targets of miR-1248. First, miR-1248 validated target genes were queried and 180 genes were obtained. Next, we compared the 180 miR-1248 validated target genes with tumour suppressor genes downregulated in TCGA cancer samples compared to normal tissue samples. We identified 10 miR-1248 validated target genes (Table 3) and literature data confirmed the oncosuppressive role for five of them in pediatric and/or adult gliomas, namely Cyclin-dependent Kinase Inhibitor 1A (CDKN1A), Fyn-related SRC Family Tyrosine kinase (FRK), S-methyl-5'-thioadenosine phosphorylase (MTAP), Speckle-type POZ protein (SPOP) and von Hippel-Lindau disease tumour suppressor (VHL). Of interest, this can potentially explain the oncogenic properties of this microRNA in the context of pLGGs.
MiR-1248 univariate analysis
Univariate analysis was performed for miR-1248 using ddPCR data to evaluate whether it could be associated with the other clinical features of PA and non-PA patients described in Table 1. The parameter estimates resulted in statistically significant results for two clinical features, site of onset and histology (Supplementary Table S4), while the association of miR-1248 expression with patient age, sex, extent of resection, histology and BRAF status was not statistically significant.
Considering the site of onset, the observed power of effect of hemispheric site of location for miR-1248 was 0.969, whereas the effect of DNET histology was 0.987. However, the number of DNET samples evaluated in ddPCR was too small, therefore we focused our attention on the hemispheric site of location and its association with miR-1248.
MiR-1248 prognostic biomarker of progressive supratentorial hemispheric pLGGs
Based on the univariate analysis results and considering that all analysed pLGG, PA and non-PA, with progression included two distinct sites of onset, midline and hemispheric, we further analysed miR-1248 expression levels subgrouping samples into these two categories.
Notably, miR-1248 was upregulated only in hemispheric tumours in comparison with midline tumours (Fig. 5a) and this differential expression was maintained when examining all the analysed histologies (Supplementary Figure S5). Moreover, ROC analysis was performed for miR-1248 in hemispheric and midline pLGG tumours with progression and highlighted the high performance of miR-1248 as prognostic biomarker in hemispheric pLGGs with progression (AUC = 1.00, 95% CI: 1.00 to 1.00; P value = 0.0002; optimal miR-1248 cut-off value for hemispheric pLGGs with progression > 2.55 copies/uL) (Fig. 5b). These results illustrate that by combining the site of onset with the expression level of miR-1248 it is possible to identify tumours with higher risk of progression after incomplete surgical resection.
Network depicting the oncogenic role of miR-1248 with its validated target genes in gliomas
Taking into consideration the information provided by other studies on glioma tumours, we propose a model (Fig. 6) recapitulating miR-1248 possible biological role in pLGGs. Indeed, the oncogenic role of miR-1248 in pLGGs emerges by examining its previously mentioned validated targets. The first, CDKN1A (p21), is a master senescence marker thus acting as a tumour suppressor . Of note, pLGGs are characterized by oncogene-induced senescence (OIS), which is associated with growth arrest and tumour progression restriction . In addition, Bongaarts et al. described low levels of p21 in a cohort of adult and pediatric GG, however they did not take into consideration the extent of resection nor the tumour progression . Therefore, the high levels of miR-1248 in supratentorial hemispheric tumours with progression may in part reflect the loss of p21 expression in this subgroup of pLGGs. The activity of other miR-1248 validated oncosuppressor target genes, namely FRK, VHL, SPOP and MTAP, has been described by different studies. The oncosuppressor role of FRK has been reported in adult gliomas [34,35,36]. FRK inhibits glioma cell migration and invasion by interacting with the N-cadherin/β-catenin and with the JNK/c-Jun signalling pathways [34, 36]. Moreover, Wang et al. recently demonstrated that FRK inhibits the activity of the Integrin subunit β1 (ITGB1), an upstream regulator of AKT, whose functional role in pLGG has been extensively described [21, 37,38,39]. VHL is involved in HIF-1-alpha degradation, β-catenin/Tcf-4 signalling pathway and AKT activity inhibition therefore exerting a pivotal role as tumour suppressor [40,41,42,43]. MTAP and SPOP activity as tumour suppressors is less clear. In detail, SPOP downregulation in adult gliomas has been associated with disease progression and has been positively correlated with mean tumour diameter, tumour grade and histological type . MTAP expression instead has been analysed in a cohort of pediatric and adult PA samples and resulted lost in tumours arisen in the cerebral hemispheres . In addition, MTAP is located on the 9p21 chromosomal region and its loss, together with that of CDKN2A, seems to be correlated with the malignant transformation of pLGGs . Altogether, these studies summarize how high levels of miR-1248 can modulate the expression of tumour suppressor genes promoting the aberrant activation of pathways that lead to cancer maintenance and progression.
pLGGs are ideally treated by complete surgical resection. However, if not feasible, patients in clinical and/or radiological progression receive chemotherapy and, in case of older children, radiotherapy, with a PFS rate in the range of 30–40% after 5 years [47, 48]. Residual disease and/or treatments cause significant morbidity . Recent studies investigated genetic alterations proposing stratification of PA and non-PA patients into risk categories and underlining the need for multiple molecular biomarkers to better define the clinical management of these tumours [9, 10].
In this scenario, our study focused on the search for epigenetic prognostic biomarkers able to predict tumour progression in the absence of a complete surgical resection. Specifically, we focused on microRNAs since they act as pivotal players in gene regulation and are deregulated in cancer . Additionally, microRNAs present a tissue specificity that may be useful to identify different tumour types and their tissue of origin , therefore acting as potential diagnostic and prognostic biomarkers .
We conducted microRNA analyses with two highly sensitive technologies using samples collected from four different neuro-oncology centres. Despite the different histology and site of onset of the tumour samples, interestingly we were able to identify three deregulated microRNAs and a new prognostic biomarker for supratentorial hemispheric pLGGs with progression. In order to corroborate these results we employed a third highly sensitive and precise method of detection. This further validation on an extended patients’ cohort strengthens the value of miR-1248 as a prognostic biomarker in incompletely resected supratentorial hemispheric pLGGs. Interestingly, its differential expression was also significant when miR-1248 was evaluated on the basis of histology. Indeed, miR-1248 was upregulated not only in supratentorial hemispheric PA with progression compared to supratentorial midline PA with progression, but also in supratentorial hemispheric non-PA with progression when compared to midline supratentorial non-PA with progression. Therefore, we demonstrated that high levels of miR-1248 combined with the site of onset can be used as a tool to classify tumours with high risk of progression after incomplete surgical resection.
A focused literature overview allowed us to hypothesize the way in which miR-1248 exerts its oncogenic role in gliomas (Fig. 6). Some of miR-1248 validated target genes indeed act as key oncosuppressor players in gliomas, therefore low levels of these genes may substantially contribute to the major aggressive behavior of incompletely resected pLGGs with progressive disease in respect to those without progression.
Additionally, a wider study of the literature showed that miR-1248 has already been described as an oncomiR correlated with cancer progression in tumours of different tissue-origin [53,54,55,56,57,58,59,60].
High levels of miR-1248 and its target gene AGTR1 were reported in osteosarcoma (OS) and associated with chemoresistance of OS cells and poor survival of patients . In esophageal squamous cell carcinoma (ESCC), miR-1248 decreased expression by UPK1A antisense RNA 1 lncRNA sponging resulted in suppressed cellular proliferation, migration and invasion, further indicating its oncogenic role . A similar tumorigenic function of miR-1248 was reported in laryngeal squamous cell carcinoma (LSCC) where low expression of circular RNA hsa_circ_0036722 could not sponge miR-1248, thus promoting proliferation by inhibiting the tumour suppressor gene RHCG . High expression levels of miR-1248 inhibited the expression of complement anaphylatoxin C3A in lung squamous carcinoma (LSC), preventing tumour growth and reducing its aggressiveness . In addition, comparison of prostate cancer patients with and without lymphatic dissemination resulted in the identification of an 18-microRNA signature associated with lymphatic dissemination that included the highly expressed miR-1248 .
Finally, miR-1248 has also been suggested as a prognostic biomarker in patients with Wilms tumour and ESCC [59, 60], conversely its role as prognostic biomarkers for progressive supratentorial hemispheric pLGGs emerges for the first time in this study.
Altogether these evidences underline the role of miR-1248 as prognostic biomarker for all the incompletely resected supratentorial hemispheric pLGGs, filling the existing void of reliable prognostic biomarkers for patient stratification and offering a tool to guide clinicians' choices for the best treatment strategy. In conclusion, although additional research with functional experiments is required to shed light on the biological activity of miR-1248, its use as prognostic biomarker in combination with tumour location, may help clinicians’ in the decision-making strategy.
Incomplete surgical resection in pediatric low-grade glioma patients entails repeated cycles of treatment, often with lifelong clinical sequelae and sometimes mortality. The lack of prognostic biomarkers provides a disadvantage for patient management. This work reports for the first time microRNA expression patterns of incompletely resected pediatric-low grade gliomas, with miR-1248 high expression levels observed in patients with progression. Further validation of these results indicated how the use of miR-1248 expression levels along with the anatomic location can be applied for the risk stratification of incompletely resected supratentorial pediatric-low grade glioma patients.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its supplementary information files.
Pediatric low-grade gliomas
Central nervous system
Dysembryoplastic neuroepithelial tumours
Near total resection
Magnetic Resonance Imaging
Progression Free Survival
TaqMan Low Density Array
Droplet Digital PCR
Gnekow AK, Kandels D, Van Tilburg C, Azizi AA, Opocher E, Stokland T, et al. SIOP-E-BTG and GPOH guidelines for diagnosis and treatment of children and adolescents with low grade glioma. Klin Padiatr. 2019;231(03):107–35.
Sturm D, Pfister SM, Jones DT. Pediatric gliomas: current concepts on diagnosis, biology, and clinical management. J Clin Oncol. 2017;73:0242.
Jones DT, Kieran MW, Bouffet E, Alexandrescu S, Bandopadhayay P, Bornhorst M, et al. Pediatric low-grade gliomas: next biologically driven steps. Neuro Oncol. 2018;20(2):160–73.
Louis DN, Perry A, Wesseling P, Brat DJ, Cree IA, Figarella-Branger D, et al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro Oncol. 2021;23(8):1231–51.
Miklja Z, Pasternak A, Stallard S, Nicolaides T, Kline-Nunnally C, Cole B, et al. Molecular profiling and targeted therapy in pediatric gliomas: review and consensus recommendations. Neuro Oncol. 2019;21(8):968–80.
Gajjar A, Bowers DC, Karajannis MA, Leary S, Witt H, Gottardo NG. Pediatric brain tumors: innovative genomic information is transforming the diagnostic and clinical landscape. J Clin Oncol. 2015;33(27):2986–98. https://doi.org/10.1200/JCO.2014.59.9217.
KiTZ Hopp Children’s Cancer Center, Heidelberg - Translational Trial Group, INFORM-2 Trial Series and LOGGIC Project. Available from: https://www.kitz-heidelberg.de/en/for-physicians/clinical-studies/molecular-diagnostics-studies/loggic-core/.
Jones DT, Banito A, Grünewald TG, Haber M, Jäger N, Kool M, et al. Molecular characteristics and therapeutic vulnerabilities across paediatric solid tumours. Nat Rev Cancer. 2019;19(8):420–38.
Yang RR, Aibaidula A, Wang WW, Chan AKY, Shi ZF, Zhang ZY, et al. Pediatric low-grade gliomas can be molecularly stratified for risk. Acta Neuropathol. 2018;136(4):641–55.
Ryall S, Zapotocky M, Fukuoka K, Nobre L, Stucklin AG, Bennett J, et al. Integrated molecular and clinical analysis of 1,000 pediatric low-grade gliomas. Cancer Cell. 2020;37(4):569-83. e5.
Ho C-Y, Bar E, Giannini C, Marchionni L, Karajannis MA, Zagzag D, et al. MicroRNA profiling in pediatric pilocytic astrocytoma reveals biologically relevant targets, including PBX3, NFIB, and METAP2. Neuro Oncol. 2012;15(1):69–82.
Ames HM, Yuan M, Vizcaíno MA, Yu W, Rodriguez FJ. MicroRNA profiling of low-grade glial and glioneuronal tumors shows an independent role for cluster 14q32. 31 member miR-487b. Mod Pathol. 2017;30(2):204–16.
Birks DK, Barton VN, Donson AM, Handler MH, Vibhakar R, Foreman NK. Survey of MicroRNA expression in pediatric brain tumors. Pediatr Blood Cancer. 2011;56(2):211–6.
Jones TA, Jeyapalan JN, Forshew T, Tatevossian RG, Lawson AR, Patel SN, et al. Molecular analysis of pediatric brain tumors identifies microRNAs in pilocytic astrocytomas that target the MAPK and NF-κB pathways. Acta Neuropathol Commun. 2015;3(1):86.
Laddha SV, Nayak S, Paul D, Reddy R, Sharma C, Jha P, et al. Genome-wide analysis reveals downregulation of miR-379/miR-656 cluster in human cancers. Biol Direct. 2013;8(1):10.
Skalsky RL, Cullen BR. Reduced expression of brain-enriched microRNAs in glioblastomas permits targeted regulation of a cell death gene. PLoS One. 2011;6(9):e24248.
Nix JS, Yuan M, Imada EL, Ames H, Marchionni L, Gutmann DH, Rodriguez FJ. Global microRNA profiling identified miR‐10b‐5p as a regulator of neurofibromatosis 1 (NF1)‐glioma migration. Neuropathology Appl Neurobiol. 2021;47(1):96–107. https://doi.org/10.1111/nan.12641.
Bongaarts A, van Scheppingen J, Korotkov A, Mijnsbergen C, Anink JJ, Jansen FE, et al. The coding and non-coding transcriptional landscape of subependymal giant cell astrocytomas. Brain. 2020;143(1):131–49.
Yuan M, Da Silva ACA, Arnold A, Okeke L, Ames H, Correa-Cerro LS, et al. MicroRNA (miR) 125b regulates cell growth and invasion in pediatric low grade glioma. Sci Rep. 2018;8(1):1–14.
Bongaarts A, Prabowo AS, Arena A, Anink JJ, Reinten RJ, Jansen FE, et al. MicroRNA519d and microRNA4758 can identify gangliogliomas from dysembryoplastic neuroepithelial tumours and astrocytomas. Oncotarget. 2018;9(46):28103.
Catanzaro G, Besharat ZM, Miele E, Chiacchiarini M, Po A, Carai A, et al. The miR‐139‐5p regulates proliferation of supratentorial paediatric low‐grade gliomas by targeting the PI3K/AKT/mTORC1 signalling. Neuropathol Appl Neurobiol. 2018;44(7):687–706. https://doi.org/10.1111/nan.12479.
Sredni ST, Huang C-C, Suzuki M, Pundy T, Chou P, Tomita T. Spontaneous involution of pediatric low-grade gliomas: high expression of cannabinoid receptor 1 (CNR1) at the time of diagnosis may indicate involvement of the endocannabinoid system. Childs Nerv Syst. 2016;32(11):2061–7.
Tantawy M, Elzayat MG, Yehia D, Taha H. Identification of microRNA signature in different pediatric brain tumors. Genet Mol Biol. 2018;41(1):27–34.
MolecularNeuropathology.org - The platform for next generation neuropathology. Available from: https://www.molecularneuropathology.org/mnp.
Capper D, Jones DT, Sill M, Hovestadt V, Schrimpf D, Sturm D, et al. DNA methylation-based classification of central nervous system tumours. Nature. 2018;555(7697):469.
Jaunmuktane Z, Capper D, Jones DT, Schrimpf D, Sill M, Dutt M, et al. Methylation array profiling of adult brain tumours: diagnostic outcomes in a large, single centre. Acta Neuropathol Commun. 2019;7(1):1–18.
Jones DT, Hutter B, Jäger N, Korshunov A, Kool M, Warnatz H-J, et al. Recurrent somatic alterations of FGFR1 and NTRK2 in pilocytic astrocytoma. Nat Genet. 2013;45(8):927.
Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019;47(D1):D155–62.
Hovestadt V, Jones DT, Picelli S, Wang W, Kool M, Northcott PA, et al. Decoding the regulatory landscape of medulloblastoma using DNA methylation sequencing. Nature. 2014;510(7506):537–41.
Huang HY, Lin YCD, Li J, Huang KY, Shrestha S, Hong HC, et al. miRTarBase 2020: updates to the experimentally validated microRNA–target interaction database. Nucleic Acids Res. 2020;48(D1):D148–54.
Zhao M, Sun J, Zhao Z. TSGene: a web resource for tumor suppressor genes. Nucleic Acids Res. 2013;41(D1):D970–6.
Abbas T, Dutta A. p21 in cancer: intricate networks and multiple activities. Nat Rev Cancer. 2009;9(6):400–14.
Milde T, Rodriguez FJ, Barnholtz-Sloan JS, Patil N, Eberhart CG, Gutmann DH. Reimagining pilocytic astrocytomas in the context of pediatric low-grade gliomas. Neuro Oncol. 2021;23(10):1634–46.
Shi Q, Song X, Wang J, Gu J, Zhang W, Hu J, et al. FRK inhibits migration and invasion of human glioma cells by promoting N-cadherin/β-catenin complex formation. J Mol Neurosci. 2015;55(1):32–41.
Wang J, Cai C, Nie D, Song X, Sun G, Zhi T, et al. FRK suppresses human glioma growth by inhibiting ITGB1/FAK signaling. Biochem Biophys Res Commun. 2019;517(4):588–95.
Zhou X, Hua L, Zhang W, Zhu M, Shi Q, Li F, et al. FRK controls migration and invasion of human glioma cells by regulating JNK/c-Jun signaling. J Neurooncol. 2012;110(1):9–19.
Hütt-Cabezas M, Karajannis MA, Zagzag D, Shah S, Horkayne-Szakaly I, Rushing EJ, et al. Activation of mTORC1/mTORC2 signaling in pediatric low-grade glioma and pilocytic astrocytoma reveals mTOR as a therapeutic target. Neuro Oncol. 2013;15(12):1604–14.
Mueller S, Phillips J, Onar-Thomas A, Romero E, Zheng S, Wiencke JK, et al. PTEN promoter methylation and activation of the PI3K/Akt/mTOR pathway in pediatric gliomas and influence on clinical outcome. Neuro Oncol. 2012;14(9):1146–52.
Pfister S, Witt O. Pediatric Gliomas. In: von Deimling A, editor. Gliomas. Springer, Berlin Heidelberg: Berlin, Heidelberg; 2009. p. 67–81.
Chen L, Han L, Zhang K, Shi Z, Zhang J, Zhang A, et al. VHL regulates the effects of miR-23b on glioma survival and invasion via suppression of HIF-1α/VEGF and β-catenin/Tcf-4 signaling. Neuro Oncol. 2012;14(8):1026–36.
Sun X, Liu M, Wei Y, Liu F, Zhi X, Xu R, et al. Overexpression of von Hippel-Lindau tumor suppressor protein and antisense HIF-1α eradicates gliomas. Cancer Gene Ther. 2006;13(4):428–35.
Zhu W, Krishna S, Garcia C, Lin CCJ, Mitchell BD, Scott KL, et al. Daam2 driven degradation of VHL promotes gliomagenesis. Elife. 2017;6:e31926.
Lee SB, Frattini V, Bansal M, Castano AM, Sherman D, Hutchinson K, et al. An ID2-dependent mechanism for VHL inactivation in cancer. Nature. 2016;529(7585):172–7.
Ding D, Song T, Jun W, Tan Z, Fang J. Decreased expression of the SPOP gene is associated with poor prognosis in glioma. Int J Oncol. 2015;46(1):333–41.
Becker AP, Scapulatempo-Neto C, Menezes WP, Clara C, Machado HR, Oliveira RS, et al. Expression of methylthioadenosine phosphorylase (MTAP) in pilocytic astrocytomas. Pathobiology. 2015;82(2):84–9.
Frazão L, do Carmo Martins M, Nunes VM, Pimentel J, Faria C, Miguéns J, et al. BRAF V600E mutation and 9p21: CDKN2A/B and MTAP co-deletions-Markers in the clinical stratification of pediatric gliomas. BMC Cancer. 2018;18(1):1–10.
Ater JL, Zhou T, Holmes E, Mazewski CM, Booth TN, Freyer DR, et al. Randomized study of two chemotherapy regimens for treatment of low-grade glioma in young children: a report from the Children’s Oncology Group. J Clin Oncol. 2012;30(21):2641.
Wisoff JH, Sanford RA, Heier LA, Sposto R, Burger PC, Yates AJ, et al. Primary neurosurgery for pediatric low-grade gliomas: a prospective multi-institutional study from the Children’s Oncology Group. Neurosurgery. 2011;68(6):1548–55.
Fangusaro J, Bandopadhayay P. The, “Risk” in Pediatric Low-Grade Glioma. Cancer Cell. 2020;37(4):424–5.
Peng Y, Croce CM. The role of MicroRNAs in human cancer. Signal Transduct Target Ther. 2016;1(1):1–9.
Rosenfeld N, Aharonov R, Meiri E, Rosenwald S, Spector Y, Zepeniuk M, et al. MicroRNAs accurately identify cancer tissue origin. Nat Biotechnol. 2008;26(4):462–9.
Condrat CE, Thompson DC, Barbu MG, Bugnar OL, Boboc A, Cretoiu D, et al. miRNAs as biomarkers in disease: latest findings regarding their role in diagnosis and prognosis. Cells. 2020;9(2):276.
Zhao Y, Xu K, Liu P. Post-transcriptional control of angiotensin II type 1 receptor regulates osteosarcoma cell death. Cell Physiol Biochem. 2018;45(4):1581–9.
Zhang Y, Wang M, Zang X, Mao Z, Chen Y, Mao F, et al. CircHN1 affects cell proliferation and migration in gastric cancer. J Clin Lab Anal. 2020;34(10):e23433.
Du F, Guo T, Cao C. Restoration of UPK1A-AS1 expression suppresses cell proliferation, migration, and invasion in esophageal squamous cell carcinoma cells partially by sponging microRNA-1248. Cancer Manag Res. 2020;12:2653.
Guo Y, Huang Q, Zheng J, Hsueh CY, Yuan X, Heng Y, et al. Diagnostic role of dysregulated circular RNA hsa_circ_0036722 in laryngeal squamous cell carcinoma. Onco Targets Ther. 2020;13:5709.
Li G, Guo X. LncRNA STARD13-AS blocks lung squamous carcinoma cells growth and movement by targeting miR-1248/C3A. Pulm Pharmacol Ther. 2020;64:101949.
Pudova EA, Krasnov GS, Nyushko KM, Kobelyatskaya AA, Savvateeva MV, Poloznikov AA, et al. miRNAs expression signature potentially associated with lymphatic dissemination in locally advanced prostate cancer. BMC Med Genomics. 2020;13(8):1–13.
Gong Y, Zou B, Chen J, Ding L, Li P, Chen J, et al. Potential five-microRNA signature model for the prediction of prognosis in patients with Wilms tumor. Med Sci Monit. 2019;25:5435.
Zhang L, Chen J, Cheng T, Yang H, Pan C, Li H. Identification of differentially expressed genes and miRNAs associated with esophageal squamous cell carcinoma by integrated analysis of microarray data. BioMed Res Int. 2020;20:1980921.
We would like to thank the AP-HM Tumour Bank (authorization number: AC2018-31053; CRB BB-0033-00097) for providing tissue samples.
This work was supported by Association “Fondo di Giò ONLUS” and Sapienza University Grant. It was also funded by grants from patients’ associations; we thank the Association Cassandra, Imagine for Margo, the ARTCSud (Association pour la recherche sur les Tumeurs Cérébrales – Sud) and the SFCE (Société Française de lutte contre les Cancers et les leucémies de l’Enfant et de l’Adolescent).
Ethics approval and consent to participate
Ethical approval (Rif. 5866) was obtained in accordance with the Helsinki declaration of 1964 and its later amendments. Informed written consent was obtained from the patients, parents or guardians before enrolment, according to our ethical committee guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Overview of the methodological workflow of the study for the identification of miR-1248 as a progression risk stratification biomarker in pLGGs. Figure S2. Progression free survival (PFS) data for pLGG patients of cohorts I and II (p-value<0.001). Figure S3. (a) Normalized relative Ct expression levels of miR-376a-3p in pLGG without and with progression samples from cohort I (*p-value=0.025). (b) Log2 expression levels of miR-376a-3p in PA without and with progression samples from cohort II (*p-value=0.028). (c) Normalized relative Ct expression levels of miR-888-5p in pLGG without and with progression samples from cohort I (*p-value=0.007). (d) Log2 expression levels of miR-888-5p in PA without and with progression samples from cohort II (*p-value=0.007). (e) Normalized relative Ct expression levels of miR-1248 in pLGG without and with progression samples from cohort I (*p-value= 0.029). (f) Log2 expression levels of miR-1248 in PA without and with progression samples from cohort II (*p-value=0.0119). Figure S4. ROC curve of miR-376a-3p in pLGG with and without progression (AUC=0.5891; p-value=0.3986). Blue line=sensitivity, red line=identity. Figure S5. MiR-1248 levels distinguish Hemispheric PA and non-PA tumours with progression from Midline ones. (a) ddPCR copies/μL of miR-1248 in all the pLGG subgroups with progression. * p<0.05 vs PA Hemispheric with progression, ° p<0.05 vs GG Hemispheric with progression (b) ddPCR copies/μL of miR-1248 in Hemispheric PA with progression vs Midline PA with progression. ** p<0.01 vs PA Hemispheric with progression (c) ddPCR copies/μL of miR-1248 in Hemispheric non-PA with progression vs Midline non-PA with progression. The blue dots refer to hemispheric GG with progression, the fuchsia dot refers to the hemispheric DNET with progression, the turquoise dots refer to the midline GG with progression. * p<0.05 vs non-PA Hemispheric with progression.
Table S1. Methylation profiling of 18 Cohort I samples. Table S2. MicroRNAs displaying significantly upregulated and downregulated expression in supratentorial PA and non-PA tumours with and without progression (Cohort I). Table S3. MicroRNAs displaying significantly upregulated and downregulated expression in supratentorial PAs with and without progression (Cohort II). Table S4. Univariate analyses of miR-1248 in pLGG in pLGG samples of ddPCR.
About this article
Cite this article
Catanzaro, G., Besharat, Z.M., Carai, A. et al. MiR-1248: a new prognostic biomarker able to identify supratentorial hemispheric pediatric low-grade gliomas patients associated with progression. Biomark Res 10, 44 (2022). https://doi.org/10.1186/s40364-022-00389-x
- Pediatric low-grade gliomas
- Prognostic biomarker
- Tumour progression
- Risk stratification
- Personalized medicine