Identification of ultrasound-associated gene candidates in myeloid cells and construction of a prognostic risk model for acute myeloid leukemia.

BACKGROUND

Incorporating ultrasound (US) treatment sensitivity analysis may improve the treatment of acute myeloid leukemia (AML).

METHODS

This study integrated single-cell and bulk datasets for analysis. Differential expression analysis between US-treated and control samples was performed using limma package. The AUCell package was used to calculate US-associated scores in the single-cell dataset. Differentially expressed genes (DEGs) between the specific groups were identified, followed by intersection analysis with previously identified DEGs. Univariate regression, Least Absolute Shrinkage and Selection Operator (LASSO) analysis (using the glmnet package), and stepwise multivariate regression (using the MASS package) were used to refine the candidate genes and to construct a risk model. The model genes were validated using in vitro experiments. Enrichment analysis was conducted using gene set enrichment analysis (GSEA), and immune infiltration was evaluate by single-sample GSEA (ssGSEA) and ESTIMATE algorithms. The correlations between RiskScores and drug sensitivity were analyzed by oncoPredict package. Finally, tumor mutational burden (TMB) and genomic mutations were compared between the risk groups.

RESULTS

Nine prognostic signatures (SPINK2, HNRNPAB, SH3BGRL3, CLEC11A, ITGA4, RPL39L, MX1, HEXIM1, and MAP4K4) were identified. Particularly, low expression of SPINK2 attenuated the activity and invasion of AML cells. High-risk group had higher immune cell infiltration. Eight drugs were predicted to be correlated with the RiskScore model. DNMT3A and RUNX1 showed higher mutation frequencies in the high-risk group, whereas KIT and MUC16 showed higher mutation frequencies in the low-risk group.

CONCLUSION

The RiskScore model established in this study provides a theoretical basis for clinically screening responsive populations and optimizing treatment strategies.