Identification of aberrant microRNA expression pattern in pediatric gliomas by microarray
- Fatao Liu†1,
- Yuyu Xiong†2,
- Yang Zhao†3,
- Liming Tao1,
- Zhou Zhang1,
- Hong Zhang4,
- Yun Liu4, 5,
- Guoyin Feng2,
- Baojie Li2,
- Lin He1, 2, 4,
- Jie Ma3Email author,
- Shengying Qin2Email author and
- Yifeng Yang1Email author
© Liu et al.; licensee BioMed Central Ltd. 2013
Received: 20 August 2013
Accepted: 9 September 2013
Published: 20 September 2013
Brain tumor remains the leading cause of disease-related death in children. Many studies have focused on the complex biological process involved in pediatric brain tumors but little is know about the possible role of microRNAs in the genesis of these tumors.
In this study, we used a microRNA microarray assay to study the expression pattern of microRNAs in pediatric gliomas and matched normal tissues.
We found 40 differentially expressed microRNAs, among which miR-1321, miR-513b, miR-769-3p were found be related to cancer genesis for the first time. The expression of selected microRNAs were then confirmed by qRT-PCR. Furthermore, GO and pathway analysis showed that the target genes of the 40 differentially expressed microRNAs were significantly enriched in nervous system-related and tumor-related biological processes and signaling pathways. Additionally, an apoptosis-related network of microRNA–mRNA interaction, representing the critical microRNAs and their targets, was constructed based on microRNA status.
In the present study we identified the changed expression pattern of microRNAs in pediatric gliamas. Our study also provides a better understanding of pediatric brain tumor biology and may assist in the development of less toxic therapies and in the search for better markers for disease stratification.
The virtual slide(s) for this article can be found here: http://www.diagnosticpathology.diagnomx.eu/vs/1323049861105720
Brain tumor, together with leukaemia, remains the leading cause of disease-related death in children . According to a population-based study by Kaatsch and colleagues , approximately 60% of pediatric brain tumors are gliomas. Glioma is the most common type of primary brain tumor. Brain glioma can cause headaches, nausea and vomiting, seizures, and cranial nerve disorders as a result of increased intracranial pressure. Based on the observations that different gliomas share morphological similarities in different lineages of glial cells, the respective tumors have been classified as astrocytomas, oligodendrogliomas, and ependymomas . However, because the pathogenesis of these tumors is unclear, treatment is particularly complex. Many children who have been treated for brain tumors experience significant long-term problems, such as changes in intellectual and motor function . A better understanding of pediatric brain tumor pathogenesis is necessary to provide better markers for disease stratification and to assist in the development of less toxic therapies.
Aberrant microRNA expression has been found to be associated with a wide variety of human tumors, including lung cancer , cervical cancer , bladder cancer , esophageal adenocarcinoma  and pituitary adenomas  et al. MicroRNAs also play a significant role in brain tumor pathology by regulating target gene expression, apoptosis and autophagy . The expression of miR-21 has been found to be increased between 5 and 100-fold in human glioblastoma tissues compared to control non-neoplastic brains . A set of brain-enriched microRNAs, miR-128, miR-181a, miR-181b and miR-181c have been found to be down-regulated in glioblastoma . Bottoni and his colleagues, using Northern blot, found that two microRNAs, miR-15a and miR-16-1 had a reduced expression in pituitary adenomas as compared to normal pituitary tissue . Upregulation of mir-372 has also been found to be related with poor prognosis in glioma .
Large profiling studies using solid tissue and hematological tumors have established the usefulness of microRNA profiling for diagnosis and prognosis . We therefore sought to determine the expression profiles of microRNAs in pediatric gliomas and matched normal tissues using microRNAs microarrays. We also performed Gene Ontology (GO) and pathway analysis to investigate the changed biological processes and signaling pathways involved in pediatric gliomas.
Materials and methods
General information of the tumor samples involved in the present study
WHO grade II
WHO grade II-III
WHO grade II
WHO grade I-II
WHO grade IV
WHO grade IV
WHO grade II-III
WHO grade II
MicroRNA microarray assay
Total RNA extraction and amplification
All of the RNA samples were extracted from tissues using the Trizol (Life technology) method. DNase I (New England Biolabs) was then added to digest the residual RNA. After purification, the concentration was measured using Nanodrop2000, all total RNA was reversely transcribed to cDNA using a MessageAmp™ II aRNA Amplification Kit (Ambion). T7 Enzyme Mix (Life technology) was then used to perform transcription of cDNAs to aRNAs in vitro. The linear amplified aRNAs were labeled using dissolved Mono-functional CyDye (Cy5) and purified by resin column.
Fluorescence labeling, hybridization and scanning
Firstly, microRNA probes (a total of 866 human microRNAs; miRBase, release 12.0) printed on the microRNA microarray (LC Sciences) were denatured at 95°C for 5 minutes. The microarray and probes were then hybridized with the Cy5-labeled aRNAs in the last step for 16 hours at 50°C. Next, 2X SSC + 0.2% SDS, 0.1% × SSC + 0.2% SDS, 0.1% × SSC (Amresco) were each used consecutively for ten minutes to wash the microarray. The microarray was then scanned and fluorescence signals were assessed.
Microarray data analysis
Data generated from the microarray was imported to Microsoft Excel. After normalizing the signal of each microRNA by using global average normalization as described by Bilban et al. , the expression level of each microRNA was calculated and student’s t test was performed to estimate between-group differences. For each microRNA, the difference between brain tumor and the matched adjacent noncancerous tissues was set to be significant if the fold change > 1.5 or < 0.67 and the p-value < 0.05.
Quantitative RT-PCR analysis
Three of the up-regulated microRNAs, miR-1259, miR-21 and miR-222 and one of the down-regulated microRNAs, miR-128, were selected at representatives for verification. MiR-1259 and miR-128 exhibited lowest p-values, miR-21 was a well-known Onco-miR and miR-222 expression was also proved to be related with the development of a variety of tumors [17–20]. Expression of these microRNAs was assayed using stem-loop RT followed by PCR analysis as previously described (qRT-PCR) . Real-time PCR was performed using TaqMan PCR kit (Life technology) on an Applied Biosystems 7900 Fast Real-Time PCR System (Applied Biosystems). In our study, qRT-PCR was performed in triplicate for each sample. The relative amount of microRNAs was normalized against U6 snRNA, the stability value of which was 0.046 according to Normfinder (http://www.multid.se/genex/hs410.htm). The fold change for each microRNA was calculated using the 2-[delta][delta] Ct method .
We performed cluster analysis on the microRNAs that showed differential expression in brain tumors using Cluster (version 3.0, http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm) software. The TreeView tool (http://jtreeview.sourceforge.net/) was used to show the results generated by Cluster 3.0.
MicroRNA target gene prediction
The sequences of microRNAs with differential expression were obtained from the Sanger microRNA Registry (http://www.sanger.ac.uk/software/Rfam/mirna). The sequences were used as query sequences to predict target genes on the UCSC human genome sequences and refgenes database (http://genome.ucsc.edu/). Only target genes that were expressed in the brain were retained.
Gene Ontology (GO) analysis
We performed GO analysis on target genes of microRNAs with differential expression based on the Gene ontology database (http://www.geneontology.org/). Fisher’s two-side exact test and the Chi-square test were used to classify the GO categories, and the false discovery rate (FDR) was calculated to correct the P-value. We chose only GOs that had a p-value of <0.05 and a FDR of < 0.05.
Pathway analysis was also performed on target genes of the differentially expressed microRNAs based on the KEGG database (http://www.genome.jp/kegg/). As in the GO analysis, two-side Fisher’s exact test and Chi-square test were used to classify the KEGG pathway categories, and the false discovery rate (FDR) was also calculated. Using a p-value cutoff of 0.05 and FDR cutoff of 0.05, we determined the enriched pathways.
MicroRNA microarray analysis
Characteristics of microRNAs that were significantly up- or down-regulated in pediatric brain tumors
MicroRNA expression validation
Character analysis of MicroRNA target genes
We determined the target genes of each microRNA based on the UCSC Genome Browser Database. 2276 target genes of up-regulated microRNAs and 2403 target genes of down-regulated microRNAs were found to be expressed in the human brain. We then performed GO analysis on target genes of up- and down-regulated microRNAs, respectively. We found 143 significant GO categories for target genes of up-regulated microRNAs and 168 significant GO categories for target genes of down-regulated microRNAs.
In the present study we used a microRNA microarray assay to study the differential expression pattern of microRNA in pediatric gliomas and the matched adjacent noncancerous tissues. Among the microRNAs detected in the microarray, we identified 40 microRNAs showing significantly higher or lower expression levels in tumors compared to the matched normal tissues. Results of GO analysis and KEGG pathway analysis suggested that target genes were closely associated with nervous system-related and tumor-related biological processes and signaling pathways.
Some of the microRNAs that showed differential expression between brain tumors and normal tissues in our study had previously been studied or found to be associated with gliomas, including pediatric gliomas. MiR-21, which was discovered to be upregulated by us, was also found and validated to be with increased expression in pediatric pilocytic astrocytoma (PA), a World Health Organization grade I pediatric glioma, in the study of Ho CY et al. . Overexpression of miR-222 in gliomas was also observed by the research of Li Q et al. . They proved that miR-222 could regulate Wnt/β-catenin signaling pathway and promote glioma genesis by using RNA interference and western blot technology. Furthermore, downregulation of miR-204, which was found to be associated with gliomas in our study, was also proved to contribute to glioma migration by targeting the migration-promoting receptor EphB2 . In the present study miR-218 was discovered to be with decreased expression in pediatric gliomas, indicating the reverse relationship between miR-218 expression and development of gliomas. This was consistent with the findings of Tu Y et al. . They found that upregulation of miR-218 reduced the migration, invasion and proliferation of glioma cells dramatically by regulating a wide range of genes and pathways. A new target gene of miR-218, which was noticed to be down-regulated in glioma in our study, has recently been found in the study of Shi ZM et al. . They identified p70S6K1 as a novel direct target of miR-128 and overexpression of p70S6K1 can partly rescue the inhibitory effect of miR-128 in the glioma cells. This result further confirms our findings.
As expected, some of the microRNAs that showed differential expression in gliomas in our data were novel. There is little evidence about the relationship between miR-1321, miR-513b, miR-769-3p and cancer genesis. In addition, the upregulated microRNAs including miR-424*, miR-760, miR-513b, miR-361-5p, miR-1259, miR-363, miR-199a, miR-1827, miR-423-5p, miR-1308, miR-1274, miR-1224, miR-513b and the downregulated microRNAs including miR-885, miR-769, miR-1296, miR-192, miR-331-3p, miR-484, miR-99b were found to be related with gliomas for the first time. This may provide new clues for gliomas research.
Aberrant expression of microRNAs may contribute to the induction of pediatric gliomas by regulating genes involved in nervous system-related and tumor-related biological processes and signaling pathways. According to the result of GO analysis, the target genes were involved in negative regulation of neuron apoptosis, axonogenesis, regulation of synaptic transmission, neurotransmitter secretion and synaptic vesicle transport et al.. These biological processes are crucial for maintaining normal functioning of the nervous system. Alteration of the processes induced by aberrant expression of certain microRNAs may favor glioma-genesis. As for biological pathways, several cancer-related pathways, including glioma, endometrial cancer, the TGF-beta signaling pathway, the MAPK signaling pathway, the wnt signaling pathway and the notch signaling pathway et al. were among the enriched pathways of the target genes. TGF-beta signaling has long been considered to contribute to glioma pathogenesis by direct support of tumor growth, self-renewal of glioma initiating stem cells and the inhibiting of anti-tumor immunity . Wnt signaling and MAPK signaling also play significant roles in almost all kinds of tumors, including glioma . In recent studies, notch signaling has been proved to be dysregulated in brain tumors and to contribute to the malignant potential of these tumors . Peng Xu et al. found out the different roles of Notch1 and Notch2. They proved that both upregulating of Notch1 and knocking down Notch2 had the effect of suppressing glioma cell growth and invasion as well as inducing apoptosis .
In conclusion, we investigated the changed expression patterns of microRNAs in pediatric gliomas and identified the possible role of microRNA and corresponding target genes in brain tumor pathology. However, due to the limitation of the microarray assay, only 866 known human microRNAs were detected. The role of more microRNAs in brain tumors needs to be investigated using more advanced techniques, such as second-generation sequencing and higher capacity microarrays. The present study does at least provide new insights into pediatric brain tumor biology and may assist in finding new diagnostic and therapeutic tools for these tumors.
The changed expression pattern of microRNAs in pediatric gliomas was investigated and it may assist in finding new diagnostic and therapeutic strategies for these tumors. The specific microRNA expression pattern provides biomarkers for the diagnosis, staging and prognosis of gliomas. The differential expressed microRNA may also provide targets for drug development. MicroRNA mimics and inhibitors can strengthen and weaken the regulation of target genes respectively and may be novel antitumor candidates.
This work was supported by 973 Program (2011CB504000, 2010CB945203, 2010CB529600, 2010CB945203), the National Key Technology R & D Program (2012BAI01B09), the National Natural Science Foundation (81121001, 81271382, 81302005), the Shanghai Leading Academic Discipline Project (B205), the Shanghai Municipality Science & Technology Commission (09DJ1400601, 2008ZD001, 07DJ14005, 09XD1403300, 13ZR1438500) and the Chinese Academy of Sciences (2009KIP305).
- Turner CD, Rey-Casserly C, Liptak CC, Chordas C: Late effects of therapy for pediatric brain tumor survivors. J Child Neurol. 2009, 24: 1455-1463. 10.1177/0883073809341709.View ArticlePubMedGoogle Scholar
- Kaatsch P, Rickert CH, Kuhl J, Schuz J, Michaelis J: Population-based epidemiologic data on brain tumors in German children. Cancer. 2001, 92: 3155-3164. 10.1002/1097-0142(20011215)92:12<3155::AID-CNCR10158>3.0.CO;2-C.View ArticlePubMedGoogle Scholar
- Pfister S, Hartmann C, Korshunov A: Histology and molecular pathology of pediatric brain tumors. J Child Neurol. 2009, 24: 1375-1386. 10.1177/0883073809339213.View ArticlePubMedGoogle Scholar
- Achanta P, Sedora Roman NI, Quinones-Hinojosa A: Gliomagenesis and the use of neural stem cells in brain tumor treatment. Anticancer Agents Med Chem. 2010, 10: 121-130. 10.2174/187152010790909290.PubMed CentralView ArticlePubMedGoogle Scholar
- Peng Y, Dai Y, Hitchcock C, Yang X, Kassis ES, Liu L, Luo Z, Sun HL, Cui R, Wei H, et al.: Insulin growth factor signaling is regulated by microRNA-486, an underexpressed microRNA in lung cancer. Proc Natl Acad Sci USA. 2013, in pressGoogle Scholar
- Shen SN, Wang LF, Jia YF, Hao YQ, Zhang L, Wang H: Upregulation of microRNA-224 is associated with aggressive progression and poor prognosis in human cervical cancer. Diagn Pathol. 2013, 8: 69-10.1186/1746-1596-8-69.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang S, Xue S, Dai Y, Yang J, Chen Z, Fang X, Zhou W, Wu W, Li Q: Reduced expression of microRNA-100 confers unfavorable prognosis in patients with bladder cancer. Diagn Pathol. 2012, 7: 159-10.1186/1746-1596-7-159.PubMed CentralView ArticlePubMedGoogle Scholar
- Luzna P, Gregar J, Uberall I, Radova L, Prochazka V, Ehrmann J: Changes of microRNAs-192, 196a and 203 correlate with Barrett's esophagus diagnosis and its progression compared to normal healthy individuals. Diagn Pathol. 2011, 6: 114-10.1186/1746-1596-6-114.PubMed CentralView ArticlePubMedGoogle Scholar
- Mao ZG, He DS, Zhou J, Yao B, Xiao WW, Chen CH, Zhu YH, Wang HJ: Differential expression of microRNAs in GH-secreting pituitary adenomas. Diagn Pathol. 2010, 5: 79-10.1186/1746-1596-5-79.PubMed CentralView ArticlePubMedGoogle Scholar
- Palumbo S, Miracco C, Pirtoli L, Comincini S: Emerging roles of microRNA in modulating cell-death processes in malignant glioma. J Cell Physiol. 2013, in pressGoogle Scholar
- Chan JA, Krichevsky AM, Kosik KS: MicroRNA-21 is an antiapoptotic factor in human glioblastoma cells. Cancer Res. 2005, 65: 6029-6033. 10.1158/0008-5472.CAN-05-0137.View ArticlePubMedGoogle Scholar
- Ciafre SA, Galardi S, Mangiola A, Ferracin M, Liu CG, Sabatino G, Negrini M, Maira G, Croce CM, Farace MG: Extensive modulation of a set of microRNAs in primary glioblastoma. Biochem Biophys Res Commun. 2005, 334: 1351-1358. 10.1016/j.bbrc.2005.07.030.View ArticlePubMedGoogle Scholar
- Bottoni A, Piccin D, Tagliati F, Luchin A, Zatelli MC, degli Uberti EC: miR-15a and miR-16-1 down-regulation in pituitary adenomas. J Cell Physiol. 2005, 204: 280-285. 10.1002/jcp.20282.View ArticlePubMedGoogle Scholar
- Li G, Zhang Z, Tu Y, Jin T, Liang H, Cui G, He S, Gao G: Correlation of microRNA-372 upregulation with poor prognosis in human glioma. Diagn Pathol. 2013, 8: 1-10.1186/1746-1596-8-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicoloso MS, Calin GA: MicroRNA involvement in brain tumors: from bench to bedside. Brain Pathol. 2008, 18: 122-129. 10.1111/j.1750-3639.2007.00119.x.View ArticlePubMedGoogle Scholar
- Bilban M, Buehler LK, Head S, Desoye G, Quaranta V: Normalizing DNA microarray data. Curr Issues Mol Biol. 2002, 4: 57-64.PubMedGoogle Scholar
- Hwang MS, Yu N, Stinson SY, Yue P, Newman RJ, Allan BB, Dornan D: miR-221/222 targets adiponectin receptor 1 to promote the epithelial-to-mesenchymal transition in breast cancer. PloS one. 2013, 8: e66502-10.1371/journal.pone.0066502.PubMed CentralView ArticlePubMedGoogle Scholar
- Walter BA, Valera VA, Pinto PA, Merino MJ: Comprehensive microRNA Profiling of Prostate Cancer. J Cancer Educ. 2013, 4: 350-357. 10.7150/jca.6394.View ArticleGoogle Scholar
- Li BL, Lu W, Lu C, Qu JJ, Yang TT, Yan Q, Wan XP: CpG island hypermethylation-associated silencing of microRNAs promotes human endometrial cancer. Cancer Cell Int. 2013, 13: 44-10.1186/1475-2867-13-44.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Zhang H, He L, Dong W, Li J, Shan Z, Teng W: Association between the expression of four upregulated miRNAs and extrathyroidal invasion in papillary thyroid carcinoma. OncoTargets Ther. 2013, 6: 281-287.View ArticleGoogle Scholar
- Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, et al.: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-10.1093/nar/gni178.PubMed CentralView ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Chen L, Zhang A, Li Y, Zhang K, Han L, Du W, Yan W, Li R, Wang Y, Wang K, et al.: MiR-24 regulates the proliferation and invasion of glioma by ST7L via beta-catenin/Tcf-4 signaling. Cancer Lett. 2013, 329: 174-180. 10.1016/j.canlet.2012.10.025.View ArticlePubMedGoogle Scholar
- Xu W, Liu M, Peng X, Zhou P, Zhou J, Xu K, Xu H, Jiang S: miR-24-3p and miR-27a-3p promote cell proliferation in glioma cells via cooperative regulation of MXI1. Int J Oncol. 2013, 42: 757-766.PubMedGoogle Scholar
- Ho CY, Bar E, Giannini C, Marchionni L, Karajannis MA, Zagzag D, Gutmann DH, Eberhart CG, Rodriguez FJ: MicroRNA profiling in pediatric pilocytic astrocytoma reveals biologically relevant targets, including PBX3, NFIB, and METAP2. Neuro-oncology. 2013, 15: 69-82. 10.1093/neuonc/nos269.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Q, Shen K, Zhao Y, He X, Ma C, Wang L, Wang B, Liu J, Ma J: MicroRNA-222 promotes tumorigenesis via targeting DKK2 and activating the Wnt/beta-catenin signaling pathway. FEBS Lett. 2013, 587: 1742-1748. 10.1016/j.febslet.2013.04.002.View ArticlePubMedGoogle Scholar
- Ying Z, Li Y, Wu J, Zhu X, Yang Y, Tian H, Li W, Hu B, Cheng SY, Li M: Loss of miR-204 expression enhances glioma migration and stem cell-like phenotype. Cancer Res. 2013, 73: 990-999. 10.1158/0008-5472.CAN-12-2895.PubMed CentralView ArticlePubMedGoogle Scholar
- Tu Y, Gao X, Li G, Fu H, Cui D, Liu H, Jin W, Zhang Y: MicroRNA-218 Inhibits Glioma Invasion, Migration, Proliferation and Cancer Stem-like Cell Self-renewal by Targeting the Polycomb Group Gene Bmi1. Cancer Res. 2013, in pressGoogle Scholar
- Shi ZM, Wang J, Yan Z, You YP, Li CY, Qian X, Yin Y, Zhao P, Wang YY, Wang XF, et al.: MiR-128 inhibits tumor growth and angiogenesis by targeting p70S6K1. PloS One. 2012, 7: e32709-10.1371/journal.pone.0032709.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaminska B, Kocyk M, Kijewska M: TGF beta signaling and its role in glioma pathogenesis. Adv Exp Med Biol. 2013, 986: 171-187. 10.1007/978-94-007-4719-7_9.View ArticlePubMedGoogle Scholar
- Gong A, Huang S: FoxM1 and Wnt/beta-catenin signaling in glioma stem cells. Cancer Res. 2012, 72: 5658-5662. 10.1158/0008-5472.CAN-12-0953.PubMed CentralView ArticlePubMedGoogle Scholar
- Stockhausen MT, Kristoffersen K, Poulsen HS: Notch signaling and brain tumors. Adv Exp Med Biol. 2012, 727: 289-304. 10.1007/978-1-4614-0899-4_22.View ArticlePubMedGoogle Scholar
- Xu P, Zhang A, Jiang R, Qiu M, Kang C, Jia Z, Wang G, Han L, Fan X, Pu P: The different role of Notch1 and Notch2 in astrocytic gliomas. PloS One. 2013, 8: e53654-10.1371/journal.pone.0053654.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.