Abstract
ETV6/RUNX1 gene fusion is the most common chromosomal translocation abnormality occurred in pediatric B-cell acute lymphoblastic leukemia (B-ALL). Compared with ETV6-RUNX1-negative patients, ETV6-RUNX1-positive patients possess more improved treatment strategies but higher risk to relapse. In this research, the potential gene interaction networks were constructed intending for elucidating the pathogenesis of B-ALL. We performed the weighted gene co-expression network analysis (WGCNA) to assess the involvement of lncRNA-mRNA pairs in B-ALL patients consisting of 24 ETV6-RUNX1-positive patients and 18 ETV6-RUNX1-negative patients and found a module that was significantly associated with positive/negative trait. Gene Ontology analysis showed that mRNAs in this module were enriched in the positive regulation of MAPK cascade, positive regulation of JNK cascade, and myeloid cell differentiation pathway. To further investigate the relationship between lncRNAs and mRNAs in this significant module, we constructed the lncRNA-mRNA co-expression network. 3 lncRNAs (RP11-170J3.2, RP11-135F9.1 and RP1-151B14.9) were found at the core of the lncRNA-mRNA co-expression network, which had the most co-expression connections with mRNAs. And several related mRNAs (ACTN1, TNFRSF21 and NLRP3) had a significant correlation with the patient survival prediction. Our findings may explicate the pathogenesis of B-ALL, and the disease-associated genes could provide clues to find novel biomarkers for prognosis.