- Open Access
Identification of aberrant microRNA expression pattern in pediatric gliomas by microarray
Diagnostic Pathologyvolume 8, Article number: 158 (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
For the study we recruited 8 patients undergoing surgery to treat astrocytomas at the XinHua Hospital. During surgery the tumor tissue and the matched adjacent noncancerous tissues were cut into small pieces and stored in liquid nitrogen. The tissues collected in XinHua Hospital were kept at −70°C until shipment to the Institute for Nutritional Sciences for RNA extraction and other experiments. Written informed consent was obtained from all patients or their representatives, and the Shanghai Committee of human rights approved the study. The general information of the tumor samples was summarized in Table 1.
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
In the microarray data analysis, we investigated more than 800 microRNAs in sum and found 40 differentially expressed microRNAs in gliomas compared to adjacent tissues. Among these, 23 up-regulated and 17 down-regulated microRNAs were involved. The information of the 40 microRNAs was summarized in Table 2. What needs illustration is that miR-24 was involved in the differential expressed microRNAs because it was with low p-value and FDR. Several recent findings suggested miR-24 played a role in the development of gliomas [23, 24]. Thus, we still listed miR-24 in Table 2 considering the fold change (1.472878) was very closed to the threshold set by us. Cluster analysis was performed of the 40 microRNAs on the 8 brain tumors and corresponding adjacent noncancerous tissues. The results are shown in Figure 1.
MicroRNA expression validation
In order to validate the microarray platform, we confirmed the expression of four randomly selected microRNAs that were strongly up- or down-regulated using qRT-PCR, using the same RNA samples that were used for the microarrays. The expression level of miR-1259, miR-21, miR-222 and miR-128 was identified. The fold change of each microRNA in tumors was calculated. The qRT-PCR results showed that miR-1259, miR-21 and miR-222 were up-regulated and that miR-128 was down-regulated, which was in line with the data generated from the microarray (Figure 2).
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.
Negative regulation of neuron apoptosis, axonogenesis, regulation of synaptic transmission, neurotransmitter secretion et al. were involved in the enriched GO categories of target genes of up-regulated microRNAs (Figure 3A), and the target genes of down-regulated microRNAs were enriched in synaptic vesicle transport, positive regulation of axon extension and positive regulation of apoptosis et al. (Figure 3B).
Pathway analysis based on the KEGG pathway database was also applied. Using a statistical method similar to that in the GO analysis, we identified the enriched pathways of the target genes. From our data, target genes of up-regulated microRNAs were involved in the ErbB signaling pathway, axon guidance, glioma, long-term potentiation, pathways in cancer, the MAPK signaling pathway et al. (Figure 4A). Target genes of down-regulated microRNAs were related to N-Glycan biosynthesis, the Wnt signaling pathway, the neurotrophin signaling pathway, the notch signaling pathway, amino acids metabolism et al. (Figure 4B).
A total of 165 target genes involved in both the enriched GO categories and KEGG pathways were assembled in a microRNA-gene network, in which the interaction between microRNAs and their corresponding target genes was displayed. The microRNAs that had crucial roles in regulating the related biological processes and pathways were identified. Among the up-regulated microRNAs, miR-24, miR-92b and miR-760 were the top three key microRNAs in the network. The top three microRNAs in the down-regulated microRNAs were miR-128, miR-218 and miR-26a. As for target genes, CACNA1E (Voltage-dependent R-type calcium channel subunit alpha-1E), FOXO3 (Forkhead box protein O3) and STX6 (Syntaxin-6) were the top three key genes in the microRNA-gene network (Figure 5).
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.
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.
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.
Pfister S, Hartmann C, Korshunov A: Histology and molecular pathology of pediatric brain tumors. J Child Neurol. 2009, 24: 1375-1386. 10.1177/0883073809339213.
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.
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 press
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.
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.
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.
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.
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 press
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.
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.
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.
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.
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.
Bilban M, Buehler LK, Head S, Desoye G, Quaranta V: Normalizing DNA microarray data. Curr Issues Mol Biol. 2002, 4: 57-64.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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 press
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.
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.
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.
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.
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.
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).
The authors declare that they have no competing interests.
JM, SYQ and YFY: conceived of the study and helped to draft the manuscript. FTL, YYX, YZ, LMT and HZ carried out part of the experiments. FTL, LMT, ZZ and YFY performed the statistical analysis. FTL wrote the manuscript. YL, GYF, BJL and LH participated in the design and coordination of the study. All authors read and approved the final manuscript.
Fatao Liu, Yuyu Xiong, Yang Zhao contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.