Identification of immunity-related lncRNAs and construction of a ceRNA network of potential prognostic biomarkers in acute myeloid leukemia

Abstract

Objective: Using bioinformatics analyses, this study aimed to identify lncRNAs related to the immune status of acute myeloid leukemia (AML) patients and ascertain the potential impact in immunity-related competing endogenous RNA (ceRNA) networks on AML prognosis.Methods: AML-related RNA-seq FPKM data, AML-related miRNA expression microarray data, and gene sets associated with immunity-related pathways were, respectively, obtained from the TCGA, GEO, and ImmReg databases. An immunity-related ceRNA network was then constructed according to the predicted interactions between AML-related mRNAs, lncRNAs, and miRNAs. After performing LASSO and multivariate Cox regression analyses, lncRNAs in the ceRNA network were used to establish an AML prognostic model. According to mutual regulatory relationships and consistent trends of expression among candidate ceRNAs, two ceRNA subnetworks related to the AML prognostic model were determined. Finally, the correlation between the expression levels of mRNAs, lncRNAs, and miRNAs in each ceRNA subnetwork and immune cell infiltration (assessed by combining the ESTIMATE and CIBERSORT methods and ssGSEA) was analyzed.Results: A total of 424 immunity-related differentially expressed (IR-DE) mRNAs (IR-DEmRNAs), 191 IR-DElncRNAs, and 69 IR-DEmiRNAs were obtained, and a ceRNA network of 20 IR-DElncRNAs, 6 IR-DEmRNAs, and 3 IR-DEmiRNAs was established. Univariate Cox regression analysis was conducted on 20 IR-DElncRNAs, and 7 of these were identified to be significantly correlated with the overall survival (OS) time in AML patients. Then, two IR-DElncRNAs (MEG3 and HCP5) were screened as independent OS-related factors by LASSO and multivariable Cox regression analyses, and a prognostic model was constructed to evaluate the survival risk in AML patients. Survival analyses indicated that the OS of patients was often poor in the high-risk group. Additionally, from this model, two ceRNA regulatory pathways, namely, MEG3/miR-125a-5p/SEMA4C and HCP5/miR-125b-5p/IL6R, which were potentially involved in the immune regulation of AML prognosis were identified.Conclusion: lncRNAs HCP5 and MEG3 may act as key ceRNAs in the pathogenesis in AML by regulating immune cell representation as part of the regulatory lncRNA-miRNA-mRNA axes. The candidate mRNAs, lncRNAs, and miRNAs included in the ceRNA network identified here may serve as useful prognostic biomarkers and immunotherapeutic targets for AML.