Front. Genet. Frontiers in Genetics Front. Genet. 1664-8021 Frontiers Media S.A. 10.3389/fgene.2021.656042 Genetics Original Research Identification of Potential Signatures and Their Functions for Acute Lymphoblastic Leukemia: A Study Based on the Cancer Genome Atlas Wang Weimin * Lyu Chunhui Wang Fei Wang Congcong Wu Feifei Li Xue Gan Silin Department of Hematology, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China

Edited by: Alessio Martino, National Research Council (CNR), Italy

Reviewed by: Michelle Lacey, Tulane University, United States; Xianghu Li, Zhejiang University, China

 *Correspondence: Weimin Wang, fccwangwm@zzu.edu.cn

This article was submitted to Statistical Genetics and Methodology, a section of the journal Frontiers in Genetics

 06 07 2021 2021 12 656042 20 01 2021 12 05 2021 Copyright © 2021 Wang, Lyu, Wang, Wang, Wu, Li and Gan. 2021 Wang, Lyu, Wang, Wang, Wu, Li and Gan

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Objective

Acute lymphoblastic leukemia (ALL) is a malignant disease most commonly diagnosed in adolescents and young adults. This study aimed to explore potential signatures and their functions for ALL.

 Methods

Differentially expressed mRNAs (DEmRNAs) and differentially expressed long non-coding RNAs (DElncRNAs) were identified for ALL from The Cancer Genome Atlas (TCGA) and normal control from Genotype-Tissue Expression (GTEx). DElncRNA–microRNA (miRNA) and miRNA–DEmRNA pairs were predicted using online databases. Then, a competing endogenous RNA (ceRNA) network was constructed. Functional enrichment analysis of DEmRNAs in the ceRNA network was performed. Protein–protein interaction (PPI) network was then constructed. Hub genes were identified. DElncRNAs in the ceRNA network were validated using Real-time qPCR.

 Results

A total of 2,903 up- and 3,228 downregulated mRNAs and 469 up- and 286 downregulated lncRNAs were identified for ALL. A ceRNA network was constructed for ALL, consisting of 845 lncRNA-miRNA and 395 miRNA–mRNA pairs. These DEmRNAs in the ceRNA network were mainly enriched in ALL-related biological processes and pathways. Ten hub genes were identified, including SMAD3, SMAD7, SMAD5, ZFYVE9, FKBP1A, FZD6, FZD7, LRP6, WNT1, and SFRP1. According to Real-time qPCR, eight lncRNAs including ATP11A-AS1, ITPK1-AS1, ANO1-AS2, CRNDE, MALAT1, CACNA1C-IT3, PWRN1, and WT1-AS were significantly upregulated in ALL bone marrow samples compared to normal samples.

 Conclusion

Our results showed the lncRNA expression profiles and constructed ceRNA network in ALL. Furthermore, eight lncRNAs including ATP11A-AS1, ITPK1-AS1, ANO1-AS2, CRNDE, MALAT1, CACNA1C-IT3, PWRN1, and WT1-AS were identified. These results could provide a novel insight into the study of ALL.

 acute lymphoblastic leukemia long non-coding RNAs functional enrichment analysis competing endogenous RNAs hub genes interaction Introduction

Acute lymphoblastic leukemia (ALL) is a malignant disease most commonly diagnosed in adolescents and young adults, especially in patients younger than 15 years. Despite significant improvements in the management of ALL, the long-term survival rate of ALL patients, especially adult patients, remains low (Jabbour et al., 2018; Richard-Carpentier et al., 2019). Therefore, it is of importance to understand the pathogenesis of ALL and identify novel diagnostic biomarkers and therapeutic targets for ALL.

LncRNA is a type of RNA longer than 200 nucleotides. Dysregulated lncRNA as tumor suppressor genes or oncogenes plays a key role in a variety of biological processes, such as cell proliferation, apoptosis, migration, and invasion. Increasing studies are focusing on the role and mechanism of lncRNA in the occurrence and development of ALL (Trimarchi et al., 2014; Arthur et al., 2017). For instance, lncRNA CASC15 could regulate SOX4 expression in RUNX1-translocated leukemia (Fernando et al., 2017). LncRNA HOTAIR is closely associated with acute leukemia patients’ poor prognosis (Zhang et al., 2016). LncRNA HOXA-AS2 induces glucocorticoid resistance by promoting ALL cell proliferation and inhibiting apoptosis (Zhao et al., 2019). Despite the fact that many studies have shown the diagnostic and prognostic values of lncRNAs in ALL, it is still required to further understand their regulatory mechanism. It has been widely accepted that lncRNAs indirectly regulate gene expression through targeted miRNAs (about 20 nucleotides) at the transcriptional or post-transcriptional level. Many miRNAs have been found to play a functional regulatory role in the development of ALL, such as miRNA-126 (Nucera et al., 2016), miRNA-155 (El-Khazragy et al., 2019), and miR-141-3p (Zhou et al., 2019). Yet, the regulatory interactions between lncRNAs and miRNAs in ALL require to be clarified.

The development of transcriptome analysis and RNA sequencing technology is increasing the possibility of identifying lncRNAs that may be involved in the pathogenesis of ALL. Moreover, further studies on the function of abnormally expressed lncRNAs may help understand the pathogenesis of ALL and provide important insights for the treatment of ALL. In this study, we comprehensively analyzed DElncRNAs and DEmRNAs in bone marrow samples of ALL. A ceRNA network was constructed for ALL on the basis of DElncRNA–miRNA and miRNA–DEmRNA pairs. DEmRNAs in the ceRNA network were significantly associated with ALL-related biological processes and pathways. Among DElncRNAs in the ceRNA network, eight lncRNAs including ATP11A-AS1, ITPK1-AS1, ANO1-AS2, CRNDE, MALAT1, CACNA1C-IT3, PWRN1, and WT1-AS were validated by Real-time qPCR, which could become potential diagnostic and therapeutic targets of ALL.

 Materials and Methods ALL Data Acquisition and Differential Expression Analysis

LncRNA and mRNA RNA-seq data of 494 bone marrows with ALL (hematopoietic and reticuloendothelial systems) were retrieved from TCGA repository1, which were derived from the IlluminaHiSeq RNA-Seq platform. All the data from three phases together, including 12 cases of phase 1, 468 cases of phase 2, and 14 cases of phase 3 were enrolled in the study. There were 321 (64.98%) males, 172 (34.82%) females, and 1 unknown (0.02%). The age distribution of the ALL group is as follows: 403 cases of 0–14 years old and 91 cases of ≥ 14 years old. All data of normal tissue samples were obtained from 407 whole blood in the Genotype-Tissue Expression (GTEx) database2. There were 265 (65.11%) males and 142 (34.89%) females in the control group. The age distribution of the control group is as follows: 34 cases of 20–29 years old, 34 cases of 30–39 years old, 72 cases of 40–49 years old, 130 cases of 50–59 years old, 132 cases of 60–69 years old, and 5 cases of 70–79 years old. Complete description of the multiple ethnicity groups, the biospecimen procurement methods, and sample fixation was provided in the GTEx official annotation. Differential expression analyses between ALL samples and normal samples were carried out using the EdgeR package in R (Robinson et al., 2010). The obtained p-values were corrected by false discovery rate (FDR). mRNAs and lncRNAs with adjusted p < 0.05 and | log 2fold change (FC)| ≥ 2 were considered as DEmRNAs and DElncRNAs. Volcano plots and heatmaps were generated using the ggplot2 and packages in R, respectively.

 ceRNA Network Construction

After identification of DElncRNAs and DEmRNAs, lncRNA–miRNA pairs were predicted by miRcode3 that provides > 10,000 lncRNAs (Jeggari et al., 2012). Then, miRNAs that targeted DEmRNAs were predicted using TargetScan4 (Agarwal et al., 2015), miRDB5 (Wang, 2008), and miRTarBase database6, which provides an experimentally validated microRNA–target interactions database (Huang et al., 2020). After integration of DElncRNA–miRNA and miRNA–DEmRNA, a ceRNA network was constructed and visualized using the Cytoscape software (version 3.5.1) (Shannon et al., 2003).

 Functional Enrichment Analyses of DEmRNAs in the ceRNA Network

Gene Ontology (GO) analysis of DEmRNAs in the ceRNA network was carried out using Database for Annotation, Visualization, and Integrated Discovery (DAVID) (Dennis et al., 2003), including biological process (BP), cellular component (CC), and molecular function (MF). Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) was analyzed using the clusterProfiler in R (Yu et al., 2012). Furthermore, the KEGG results were visualized using the Cytoscape plug-in ClueGO. p < 0.05 was set as the cutoff value.

 PPI Network

The interactions between proteins were predicted using the Search Tool for the Retrieval of Interacting Genes (STRING) database7 (minimum required interaction score > 0.4) (Szklarczyk et al., 2019). Furthermore, PPI networks were embodied using the Cytoscape v3.5.0 software. In addition, we used Molecular Complex Detection (MCODE) plugin to identify the hub genes in the PPI network. The criteria were set as follows: MCODE scores > 3 and number of nodes > 4. The top 10 hub genes were identified using the ranking method of degree.

 Real-Time qPCR

Bone marrow samples were isolated from 25 ALL patients and 15 healthy participants and red blood cells were removed. Total RNA was extracted from bone marrow samples and then was stored at −80℃. Extracted samples were lysed using 1 ml of Trizol and placed for 5 min on ice. RNA concentration and purity were determined using a NanoDrop UV spectrophotometer. Then, RNA was reverse transcribed into cDNA. Primer sequences of ATP11A-AS1, ITPK1-AS1, ANO1-AS2, CRNDE, MALAT1, CACNA1C-IT3, PWRN1, and WT1-AS were designed and synthesized by Shanghai Shengong Biological Engineering Co., Ltd. (Shanghai, China). The primer sequences are listed in Table 1. PCR amplification had the following conditions: 95℃ for 3 min; 40 PCR cycle reactions (95℃ for 20 s; 60℃ for 30 s). GAPDH was used as a control. Relative expression levels of lncRNAs were calculated using 2–ΔΔCT method. Differences between the two groups were analyzed using Student’s t-test. p-value < 0.05 was considered statistically significant.

Primer sequence information for Real-time qPCR.

Gene symbol Primer sequence (5′–3′)
Human GAPDH 5′-CGGAGTCAACGGATTTGGTCGTAT-3′ (forward) 5′-AGCCTTCTCCATGGTGGTGAAGAC-3′ (reverse)
Human ANO1-AS2 5′-CCGGAACAAGAACCTCGCTC-3′ (forward) 5′-GGTCCTCGCCTACCATCCAA-3′ (reverse)
Human PWRN1 5′-ACATTCGAAACCCAGGTGCC-3′ (forward) 5′-GGAAGTGGATGCTGACGCTC-3′ (reverse)
Human MALAT1 5′-GGTTCAGAAGGTCTGAAGCTC-3′ (forward) 5′-CCCAGAAGTGTTTACACTGCT-3′ (reverse)
Human CACNA1C-IT3 5′-GCCAGGACCAAGACACCAAGAC-3′ (forward) 5′-TTGGGCAGGGCTCGGTTCC-3′ (reverse)
Human ITPK1-AS1 5′-AATCCTGTGCGCTGTCATCC-3′ (forward) 5′-GATTGCTCTTGGCTGTGCCT-3′ (reverse)
Human ATP11A-AS2 5′-ACAGTCCCTTCCCTTACGCT-3′ (forward) 5′-TGAACGCTGCACTTGTGGAC-3′ (reverse)
Human CRNDE 5′-GAGGACGTGCTGGGGCT-3′ (forward)
5′-CTGAGTCCATGTCCCGAATC-3′(reverse)
Human WT1-AS 5′-GCCTCTCTGTCCTCTTCTTTGT-3′ (forward)
5′-GCTGTGAGTCCTGGTGCTTAG-3′ (reverse)
Results Identification of DElncRNAs and DEmRNAs for ALL

The workflow of this study is shown in Figure 1. According to adjusted p-value < 0.05 and | log2 fold change (FC)| ≥ 2, 2903 up− and 3228 downregulated mRNAs were identified for ALL, as shown in the volcano plot (Figure 2A and Supplementary Material 1). Furthermore, there were 469 up− and 286 downregulated lncRNAs for ALL (Figure 2B and Supplementary Material 2). Heatmaps depicted the differences in expression patterns of all DEmRNAs (Figure 2C) and DElncRNAs (Figure 2D) between ALL bone marrow samples and normal samples.

The workflow of this study.

Identification of DEmRNAs and DElncRNAs for ALL. Volcano plot showing DEmRNAs (A) and DElncRNAs (B) for ALL. As shown in heatmaps, the differences in expression patterns of DEmRNAs (C) and DElncRNAs (D) between ALL bone marrow samples and normal samples. Red represents upregulation and green represents downregulation. DEmRNAs: differentially expressed mRNAs; DElncRNAs: Differentially expressed lncRNAs; ALL, Acute lymphoblastic leukemia.

 Construction of ceRNA Network for ALL

The miRNAs that targeted DEmRNAs were predicted using TargetScan, miRDB, and miRTarBase databases. After integration of prediction results from the three databases, 297 DEmRNAs were intersected and identified for the construction of ceRNA network (Figure 3). Furthermore, DElncRNA–miRNA relationships were predicted using miRcode database. By comprehensively analyzing DElncRNA–miRNA and miRNA–DEmRNA pairs, a ceRNA network was constructed for ALL (Figure 4). There were 845 lncRNA–miRNA pairs (Supplementary Material 3) and 395 miRNA–mRNA pairs (Supplementary Material 4) in the ceRNA network.

Venn diagram showing 297 differentially expressed mRNAs targeted by miRNAs via intersection of prediction results of TargetScan, miRDB, and miRTarBase database.

A ceRNA network construction for acute lymphoblastic leukemia. Blue rhombus represents lncRNAs; green circle represents miRNAs and red triangle represents mRNAs.

 Functional Enrichment Analysis of DEmRNAs in the ceRNA Network

As depicted in heatmaps, there were obvious differences in the expression patterns of all DEmRNAs in the ceRNA network between ALL bone marrow samples and normal samples (Figure 5A). Bubble diagrams showed the top 40 GO enrichment analysis results enriched by DEmRNAs in the ceRNA network (Figure 5B). We found that these mRNAs were mainly enriched in ALL-related biological processes such as transcription, programmed cell death, apoptosis, cell cycle, proliferation, and so on. Figure 6 depicted the relationships between DEmRNAs and enriched biological processes including morphogenesis of an epithelium, kidney epithelium development, ureteric bud development, mesonephric epithelium development, and mesonephric tubule development. As for KEGG pathway enrichment analysis results, these DEmRNAs were mainly enriched in pathways in cancer, cell cycle, small cell lung cancer, p53 signaling pathway, Wnt signaling pathway, pentose phosphate pathway, and non-small cell lung cancer (Figures 7A,B).

GO enrichment analysis of DEmRNAs in the ceRNA network. (A) Heatmaps showing differences in the expression patterns of all DEmRNAs in the ceRNA network between ALL bone marrow samples and normal samples. Red stands for upregulation and green stands for downregulation. (B) The top 40 GO enrichment analysis results including biological process, cellular component, and molecular function. DEmRNAs: Differentially expressed mRNAs.

The top five biological processes enriched by DEmRNAs in the ceRNA network. DEmRNAs: Differentially expressed mRNAs.

KEGG pathway enrichment analysis of DEmRNAs in the ceRNA network. (A) Seven enriched KEGG pathways. (B) Visualization of KEGG enrichment analysis results. DEmRNAs: Differentially expressed mRNAs.

 Identification of Hub Genes in the PPI Network

The DEmRNAs in the ceRNA network were imported into STRING database. Then, a PPI network was constructed for ALL (Figure 8A). Two PPI subnetworks were then constructed (Figures 8B,C). Ten hub genes were identified for ALL, including SMAD3, SMAD7, SMAD5, ZFYVE9, FKBP1A, FZD6, FZD7, LRP6, WNT1, and SFRP1.

Identification of hub genes in the PPI network. (A) Construction of PPI network based on the DEmRNAs in the ceRNA network. (B,C) Two PPI subnetworks for ALL. Red represents upregulation and green represents downregulation. DEmRNAs: Differentially expressed mRNAs.

 Correlation Between Hub Genes and DElncRNAs

Correlation analysis between hub genes and DElncRNAs was performed by corrplot package. The significant correlations between DElncRNAs and hub genes are shown in Figure 9 and Supplementary Material 5. There was strong correlation between WT1-AS and FZD7 (r = 0.751907203; p < 0.0001). Furthermore, PWRN1 and SMAD3 were significantly correlated (r = 0.521493415 and p = 4.32E−08).

Heatmaps showing the correlation between hub genes and DElncRNAs. DEmRNAs: Differentially expressed mRNAs. The right bar indicates the color legend of Pearson correlation values.

 Validation of Eight lncRNAs in Bone Marrow of ALL

Among all DElncRNAs in the ceRNA network, the most significant difference between 10 lncRNAs (AC009093.1, C17orf77, ATP11A-AS1, ITPK1-AS1, ANO1-AS2, ITCH-IT1, CRNDE, MALAT1, CACNA1C-IT3, and PWRN1) in the ceRNA network and WT1-AS (which was closely related to hub gene FZD7) was selected for verification. However, as the primers of AC009093.1, C17orf77, and ITCH-IT1 for RQ-PCR were not ideal, the remaining eight lncRNAs were validated. As Figure 10 shows, these eight lncRNAs were significantly upregulated in ALL bone marrow samples (n = 25) compared to normal samples (n = 15) by Real-time qPCR.

Validation of eight lncRNAs in ALL bone marrow samples using Real-time qPCR. (A) ATP11A-AS1; (B) ITPK1-AS1; (C) ANO1-AS2; (D) CRNDE; (E) MALAT1; (F) CACNA1C-IT3; (G) PWRN1; and (H) WT1-AS. Control: n = 15; ALL: n = 25.

 Discussion

In this study, we constructed a ceRNA network for ALL based on DElncRNA–miRNA and miRNA–DEmRNA relationships. Among all DElncRNAs in the ceRNA network, eight lncRNAs were validated in ALL bone marrow samples using Real-time qPCR. These lncRNAs might become potential biomarkers for ALL.

To explore potential functions of DEmRNAs in the ceRNA network, we performed functional enrichment analysis. We found that these mRNAs were mainly enriched in ALL-related biological processes such as transcription (Gocho and Yang, 2019), programmed cell death (Hass et al., 2016), apoptosis, cell cycle (Jing et al., 2018), and proliferation (Sun et al., 2019). The DEmRNAs in these biological processes could modulate the development of ALL. Furthermore, these DEmRNAs were significantly associated with pathways in cancer, cell cycle, p53 signaling pathway, Wnt signaling pathway, and pentose phosphate pathway. It has been widely accepted that the p53 signaling pathway is a promising drug target in ALL (Trino et al., 2016). In particular, alterations of the tumor suppressor gene TP53 were frequently found in pediatric ALL (Demir et al., 2020). As for the Wnt signaling pathway, it was significantly correlated with the pathogenesis of ALL (Montano et al., 2018). Recent findings reported that inhibiting Wnt/β catenin could reverse multidrug resistance in children ALL (Fu et al., 2019). Moreover, the pathway is regulated by many factors. For example, miR-181a-5p could promote ALL cell proliferation via targeting the Wnt pathway (Lyu et al., 2017). Our results indicated that the DEmRNAs in the ceRNA network could be involved in the pathogenesis of ALL.

We constructed a PPI network for B-ALL on the basis of DEmRNAs in the ceRNA network. Ten hub genes were identified for ALL, including SMAD3, SMAD7, SMAD5, ZFYVE9, FKBP1A, FZD6, FZD7, LRP6, WNT1, and SFRP1. Among them, the loss of the Smad3 protein has been identified as a key feature of acute T-cell lymphoblastic leukemia (Wolfraim et al., 2004). Smad7 is a promising therapeutic target for B-cell ALL (Guo et al., 2018). Furthermore, microRNA-181a might regulate its expression for pediatric ALL (Nabhan et al., 2017). Wnt signaling pathway can enhance hematopoietic cell proliferation (Doubravska et al., 2008). It could mediate growth and prognosis of B-cell progenitor ALL, which could be a potential treatment strategy in ALL (Khan et al., 2007; Mochmann et al., 2011). In the pathway, FZD6, FZD7, LRP6, and WNT1 were marker proteins. LRP6 has been reported to be a candidate tumor suppressor gene in pre-B ALL (Montpetit et al., 2004). Furthermore, low expression of SFRP1 was significantly associated with clinical outcomes of patients with Philadelphia-positive ALL (Martin et al., 2008).

Consistently with differential expression analysis results, eight lncRNAs including ATP11A-AS1, ITPK1-AS1, ANO1-AS2, CRNDE, MALAT1, CACNA1C-IT3, PWRN1, and WT1-AS were significantly upregulated in ALL bone marrow, indicating that these abnormally expressed lncRNAs could be involved in the development of ALL. Among them, CRNDE was upregulated in the bone marrow of B-cell precursor acute lymphoblastic leukemia (BCP-ALL) patients and BCP-ALL cell lines (NALM-6 and RS4;11). Functionally, CRNDE upregulated CREB expression by suppressing miR-345-5p, thus promoting cell proliferation and reducing cell apoptosis in BCP-ALL (Wang W. et al., 2020). A large amount of research has reported that aberrantly expressed MALAT1 was involved in a variety of cancers, such as breast cancer metastasis (Kim et al., 2018), colon cancer (Wu et al., 2018), and non-small cell lung cancer (Li et al., 2018). Abnormally expressed is in significant association with poor prognosis in childhood ALL (Pouyanrad et al., 2019). Furthermore, miR-125b in combination with miR-99a and/or miR-100 could inhibit the expression of MALAT1 in vincristine-resistant children ALL cells (Moqadam et al., 2013). PWRN1 was significantly underexpressed in gastric cancer tissues and cells (Chen et al., 2018). Overexpressed PWRN1 could inhibit the proliferation and metastasis of gastric cancer cells and tumor growth. Furthermore, PWRN1 may regulate miR-425-5p expression by acting as its sponge in gastric cancer cells. ITPK1-AS1 expression could predict gastric cancer patients’ survival (Hu et al., 2019). WT1-AS has been characterized as a tumor-suppressive lncRNA in several cancers including cervical squamous cell carcinoma (Zhang et al., 2019), gastric cancer (Du et al., 2016), papillary thyroid carcinoma (Le et al., 2020), non-small cell lung cancer cell (Jiang et al., 2020), and hepatocellular carcinoma (Lv et al., 2015). Besides, WT1-AS can regulate WT1 on oxidative stress injury and apoptosis of neurons in Alzheimer’s disease via inhibition of the miR-375/SIX4 axis (Wang Q. et al., 2020). However, other lncRNAs have not been reported yet. According to our results, these lncRNAs deserve more research on ALL.

However, there are several limitations in this study. First, since there was no normal control of ALL in TCGA database, data of 407 whole blood in the GTEx database were obtained as control. Given that ALL primarily affects younger individuals, the age distribution of the control group is not ideally matched with the ALL from TCGA database, which may cause confounding. Second, the sample size of this study is small, and larger clinical samples should be used to verify these lncRNAs. In addition, this study lacks functional experiments. In future research, we will further the function and clinical value of these lncRNAs in ALL.

 Conclusion

In our study, a ceRNA network was constructed for ALL. Among all DElncRNAs in the ceRNA network, eight lncRNAs were validated in ALL bone marrow samples using Real-time qPCR, which might provide a novel insight into the further study of ALL.

 Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

 Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Research Ethics Committee of the First Affiliated Hospital of Zhengzhou University (2019-KY-194). The study was conducted in accordance with the Declaration of Helsinki. The patients/participants provided their written informed consent to participate in this study.

 Author Contributions

WW designed the project, proposed the research concept, and wrote the manuscript. WW, CL, and FWa constructed the bioinformatic analysis, constructed the graphic images and data charts, and performed the statistical processing. CW, FWu, XL, and SG performed the experiments. All authors read and approved the manuscript and agreed to be accountable for all aspects of the research in ensuring that the accuracy and integrity of any part of the work are appropriately investigated and resolved.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding. This study was supported by grants from the National Natural Science Foundation of China (Grant No. 81700138 to WW) and the Medical Science and Technique Foundation of Health Commission of Henan Provincial (Grant No. SBGJ202003041 to WW).

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.656042/full#supplementary-material

Differentially expressed mRNAs.

Differentially expressed long non-coding RNAs.

845 lncRNA–miRNA pairs in the ceRNA network.

395 miRNA–mRNA pairs in the ceRNA network.

The specific p-values of correlation analysis between hub genes and DElncRNAs.

 References Agarwal V. Bell G. W. Nam J. Bartel D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. Elife 4:e05005. 10.7554/eLife.05005 26267216 Arthur G. Almamun M. Taylor K. (2017). Hypermethylation of antisense long noncoding RNAs in acute lymphoblastic leukemia. Epigenomics 9 635–645. 10.2217/epi-2016-0156 28093925 Chen Z. Ju H. Yu S. Zhao T. Jing X. Li P. (2018). Prader-Willi region non-protein coding RNA 1 suppressed gastric cancer growth as a competing endogenous RNA of miR-425-5p. Clin. Sci. 132 1003–1019. 10.1042/CS20171588 29535266 Demir S. Boldrin E. Sun Q. Hampp S. Tausch E. Eckert C. (2020). Therapeutic targeting of mutant p53 in pediatric acute lymphoblastic leukemia. Haematologica 105 170–181. 10.3324/haematol.2018.199364 31073076 Dennis G. Sherman B. T. Hosack D. A. Yang J. Gao W. Lane H. C. (2003). DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 4:R609. 10.1186/gb-2003-4-9-r60 Doubravska L. Simova S. Cermak L. Valenta T. Korinek V. Andera L. (2008). Wnt-expressing rat embryonic fibroblasts suppress Apo2L/TRAIL-induced apoptosis of human leukemia cells. Apoptosis 13 573–587. 10.1007/s10495-008-0191-z 18347988 Du T. Zhang B. Zhang S. Jiang X. Zheng P. Li J. (2016). Decreased expression of long non-coding RNA WT1-AS promotes cell proliferation and invasion in gastric cancer. Biochim. Biophys. Acta 1862 12–19. 10.1016/j.bbadis.2015.10.001 26449525 El-Khazragy N. Noshi M. A. Abdel-Malak C. Zahran R. F. Swellam M. (2019). miRNA-155 and miRNA-181a as prognostic biomarkers for pediatric acute lymphoblastic leukemia. J. Cell. Biochem. 120 6315–6321. 10.1002/jcb.27918 30362157 Fernando T. R. Contreras J. R. Zampini M. Rodriguez-Malave N. I. Alberti M. O. Anguiano J. (2017). The lncRNA CASC15 regulates SOX4 expression in RUNX1-rearranged acute leukemia. Mol. Cancer 16:126. 10.1186/s12943-017-0692-x 28724437 Fu J. Si L. Zhuang Y. Zhang A. Sun N. Li D. (2019). Wnt/betacatenin inhibition reverses multidrug resistance in pediatric acute lymphoblastic leukemia. Oncol. Rep. 41 1387–1394. 10.3892/or.2018.6902 30535484 Gocho Y. Yang J. J. (2019). Genetic defects in hematopoietic transcription factors and predisposition to acute lymphoblastic leukemia. Blood 134 793–797. 10.1182/blood.2018852400 31311817 Guo Y. Fang Q. Ma D. Yu K. Cheng B. Tang S. (2018). Up-regulation of HO-1 promotes resistance of B-cell acute lymphocytic leukemia cells to HDAC4/5 inhibitor LMK-235 via the Smad7 pathway. Life Sci. 207 386–394. 10.1016/j.lfs.2018.06.004 29886060 Hass C. Belz K. Schoeneberger H. Fulda S. (2016). Sensitization of acute lymphoblastic leukemia cells for LCL161-induced cell death by targeting redox homeostasis. Biochem. Pharmacol. 105 14–22. 10.1016/j.bcp.2016.01.004 26774450 Hu Z. Yang D. Tang Y. Zhang X. Wei Z. Fu H. (2019). Five-long non-coding RNA risk score system for the effective prediction of gastric cancer patient survival. Oncol. Lett. 17 4474–4486. 10.3892/ol.2019.10124 30988816 Huang H. Lin Y. Li J. Huang K. Shrestha S. Hong H. (2020). miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 48 D148–D154. 10.1093/nar/gkz896 31647101 Jabbour E. Pui C. Kantarjian H. (2018). Progress and innovations in the management of adult acute lymphoblastic leukemia. Jama Oncol. 4 1413–1420. 10.1001/jamaoncol.2018.1915 29931220 Jeggari A. Marks D. S. Larsson E. (2012). miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics 28 2062–2063. 10.1093/bioinformatics/bts344 22718787 Jiang X. Wang J. Fang L. (2020). LncRNA WT1-AS over-expression inhibits non-small cell lung cancer cell stemness by down-regulating TGF-beta1. BMC Pulm. Med. 20:113. 10.1186/s12890-020-1146-6 32349718 Jing D. Huang Y. Liu X. Sia K. C. S. Zhang J. C. Tai X. (2018). Lymphocyte-specific chromatin accessibility pre-determines glucocorticoid resistance in acute lymphoblastic leukemia. Cancer Cell 34:906. 10.1016/j.ccell.2018.11.002 30537513 Khan N. I. Bradstock K. F. Bendall L. J. (2007). Activation of Wnt/beta-catenin pathway mediates growth and survival in B-cell progenitor acute lymphoblastic leukaemia. Br. J. Haematol. 138 338–348. 10.1111/j.1365-2141.2007.06667.x 17614820 Kim J. Piao H. Kim B. Yao F. Han Z. Wang Y. (2018). Long noncoding RNA MALAT1 suppresses breast cancer metastasis. Nat. Genet. 50:1705. 10.1038/s41588-018-0252-3 30349115 Le F. Luo P. Ouyang Q. Zhong X. (2020). LncRNA WT1-AS downregulates survivin by upregulating miR-203 in papillary thyroid carcinoma. Cancer Manag. Res. 12 443–449. 10.2147/CMAR.S232294 32021456 Li S. Mei Z. Hu H. Zhang X. (2018). The lncRNA MALAT1 contributes to non-small cell lung cancer development via modulating miR-124/STAT3 axis. J. Cell. Physiol. 233 6679–6688. 10.1002/jcp.26325 29215698 Lv L. Chen G. Zhou J. Li J. Gong J. (2015). WT1-AS promotes cell apoptosis in hepatocellular carcinoma through down-regulating of WT1. J. Exp. Clin. Cancer Res. 34:119. 10.1186/s13046-015-0233-7 26462627 Lyu X. Li J. Yun X. Huang R. Deng X. Wang Y. (2017). miR-181a-5p, an inducer of Wnt-signaling, facilitates cell proliferation in acute lymphoblastic leukemia. Oncol. Rep. 37 1469–1476. 10.3892/or.2017.5425 28184923 Martin V. Agirre X. Jimenez-Velasco A. Jose-Eneriz E. S. Cordeu L. Garate L. (2008). Methylation status of Wnt signaling pathway genes affects the clinical outcome of Philadelphia-positive acute lymphoblastic leukemia. Cancer Sci. 99 1865–1868. 10.1111/j.1349-7006.2008.00884.x 18549404 Mochmann L. H. Bock J. Ortiz-Tanchez J. Schlee C. Bohne A. Neumann K. (2011). Genome-wide screen reveals WNT11, a non-canonical WNT gene, as a direct target of ETS transcription factor ERG. Oncogene 30 2044–2056. 10.1038/onc.2010.582 21242973 Montano A. Forero-Castro M. Marchena-Mendoza D. Benito R. Maria Hernandez-Rivas J. (2018). New challenges in targeting signaling pathways in acute lymphoblastic leukemia by NGS approaches: an update. Cancers 10:1104. 10.3390/cancers10040110 29642462 Montpetit A. Larose J. Boily G. Langlois S. Trudel N. Sinnett D. (2004). Mutational and expression analysis of the chromosome 12p candidate tumor suppressor genes in pre-B acute lymphoblastic leukemia. Leukemia 18 1499–1504. 10.1038/sj.leu.2403441 15284860 Moqadam F. A. Lange-Turenhout E. A. M. Aries I. M. Pieters R. den Boer M. L. (2013). MiR-125b, miR-100 and miR-99a co-regulate vincristine resistance in childhood acute lymphoblastic leukemia. Leukemia Res. 37 1315–1321. 10.1016/j.leukres.2013.06.027 23915977 Nabhan M. Louka M. L. Khairy E. Tash F. Ali-Labib R. El-Habashy S. (2017). MicroRNA-181a and its target Smad 7 as potential biomarkers for tracking child acute lymphoblastic leukemia. Gene 628 253–258. 10.1016/j.gene.2017.07.052 28732737 Nucera S. Giustacchini A. Boccalatte F. Calabria A. Fanciullo C. Plati T. (2016). miRNA-126 orchestrates an oncogenic program in b cell precursor acute lymphoblastic leukemia. Cancer Cell 29 905–921. 10.1016/j.ccell.2016.05.007 27300437 Pouyanrad S. Rahgozar S. Ghodousi E. S. (2019). Dysregulation of miR-335-3p, targeted by NEAT1 and MALAT1 long non-coding RNAs, is associated with poor prognosis in childhood acute lymphoblastic leukemia. Gene 692 35–43. 10.1016/j.gene.2019.01.003 30639603 Richard-Carpentier G. Kantarjian H. Jabbour E. (2019). Recent advances in adult acute lymphoblastic leukemia. Curr. Hematol. Malig. Rep. 14 106–118. 10.1007/s11899-019-00503-1 30879177 Robinson M. D. McCarthy D. J. Smyth G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26 139–140. 10.1093/bioinformatics/btp616 19910308 Shannon P. Markiel A. Ozier O. Baliga N. S. Wang J. T. Ramage D. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13 2498–2504. 10.1101/gr.1239303 14597658 Sun H. Zhang Z. Luo W. Liu J. Lou Y. Xia S. (2019). NET1 enhances proliferation and chemoresistance in acute lymphoblastic leukemia cells. Oncol. Res. 27 935–944. 10.3727/096504019X15555388198071 31046876 Szklarczyk D. Gable A. L. Lyon D. Junge A. Wyder S. Huerta-Cepas J. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 D607–D613. 10.1093/nar/gky1131 30476243 Trimarchi T. Bilal E. Ntziachristos P. Fabbri G. Dalla-Favera R. Tsirigos A. (2014). Genome-wide mapping and characterization of notch-regulated long noncoding RNAs in acute leukemia. Cell 158 593–606. 10.1016/j.cell.2014.05.049 25083870 Trino S. De Luca L. Laurenzana I. Caivano A. Del Vecchio L. Martinelli G. (2016). P53-MDM2 pathway: evidences for a new targeted therapeutic approach in b-acute lymphoblastic leukemia. Front. Pharmacol. 7:491. 10.3389/fphar.2016.00491 28018226 Wang Q. Ge X. Zhang J. Chen L. (2020). Effect of lncRNA WT1-AS regulating WT1 on oxidative stress injury and apoptosis of neurons in Alzheimer’s disease via inhibition of the miR-375/SIX4 axis. Aging (Albany NY) 12 23974–23995. 10.18632/aging.104079 33234729 Wang W. Wu F. Ma P. Gan S. Li X. Chen L. (2020). LncRNA CRNDE promotes the progression of B-cell precursor acute lymphoblastic leukemia by targeting the miR-345-5p/CREB axis. Mol. Cells 43 718–727. 10.14348/molcells.2020.0065 32868489 Wang X. (2008). miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA 14 1012–1017. 10.1261/rna.965408 18426918 Wolfraim L. A. Fernandez T. M. Mamura M. Fuller W. L. Kumar R. Cole D. E. (2004). Loss of Smad3 in acute T-cell lymphoblastic leukemia. N. Engl. J. Med. 351 552–559. 10.1056/NEJMoa031197 15295048 Wu Q. Meng W. Jie Y. Zhao H. (2018). LncRNA MALAT1 induces colon cancer development by regulating miR-129-5p/HMGB1 axis. J. Cell. Physiol. 233 6750–6757. 10.1002/jcp.26383 29226325 Yu G. Wang L. Han Y. He Q. (2012). clusterProfiler: an r package for comparing biological themes among gene clusters. Omics 16 284–287. 10.1089/omi.2011.0118 22455463 Zhang Y. Huang S. Zhou H. Chen C. Tian L. Shen J. (2016). Role of HOTAIR in the diagnosis and prognosis of acute leukemia. Oncol. Rep. 36 3113–3122. 10.3892/or.2016.5147 27748863 Zhang Y. Na R. Wang X. (2019). LncRNA WT1-AS up-regulates p53 to inhibit the proliferation of cervical squamous carcinoma cells. BMC Cancer 19:1052. 10.1186/s12885-019-6264-2 31694597 Zhao Q. Zhao S. Li J. Zhang H. Qian C. Wang H. (2019). TCF7L2 activated HOXA-AS2 decreased the glucocorticoid sensitivity in acute lymphoblastic leukemia through regulating HOXA3/EGFR/Ras/Raf/MEK/ERK pathway. Biomed. Pharmacother. 109 1640–1649. 10.1016/j.biopha.2018.10.046 30551418 Zhou R. Mo W. Wang S. Zhou W. Chen X. Pan S. (2019). miR-141-3p and TRAF5 network contributes to the progression of t-cell acute lymphoblastic leukemia. Cell Transplan. 281(1Suppl) 59S–65S. 10.1177/0963689719887370 31722554 Abbreviations ALL

acute lymphoblastic leukemia

 DEmRNAs

differentially expressed mRNAs

 DElncRNAs

differentially expressed long non-coding RNAs

 TCGA

The Cancer Genome Atlas

 GTEx

Genotype-Tissue Expression

 ceRNA

competing endogenous RNA

 PPI

protein–protein interaction

 FC

fold change

 GO

Gene Ontology

 BP

biological process

 CC

cellular component

 MF

molecular function

 KEGG

Kyoto Encyclopedia of Genes and Genomes.

https://portal.gdc.cancer.gov/

https://gtexportal.org/home/datasets

http://www.mircode.org/

http://www.targetscan.org/

http://www.mirdb.org/

http://mirtarbase.mbc.nctu.edu.tw/

http://string-db.org/