Integrated GSE214207, GSE181989 and GSE145668 (130,159 cells/35 untreated bone marrow and PBMC samples) to map immune and hematopoietic alterations in aplastic anemia. The workflow included Harmony integration and lineage-resolved annotation with AUCell pathway scoring, MiloR/edgeR differential abundance and expression, Monocle 3 trajectory inference, and CellChat ligand–receptor analysis to identify reproducible disease-associated immune states, progenitor dysfunction, altered differentiation, and disrupted intercellular signaling across cohorts.