Toxicity in hematopoietic stem cells from bone marrow and peripheral blood in mice after benzene exposure: Single-cell transcriptome sequencing analysis

Benzene is a ubiquitous, occupational, and environmental hematotoxic and leukemogen. Damage to hematopoietic stem cells (HSCs) induced by benzene and its metabolites is a key event in bone marrow (BM) depression and leukemogenesis. There are no reports on transcriptome profiles of HSCs following benzene exposure. Here, Smart-seq2 single-cell transcriptome sequencing was used to detect transcriptomic alternations in BM HSCs and peripheral blood HSCs (PBSCs) in male C57B/6 mice exposed to benzene. We found that benzene caused hematotoxicity which was confirmed by routine blood test, pathological examination, and HSCs percentage analysis. A total of 1514 differentially expressed genes (DEGs) in BM HSCs and 1703 DEGs in PBSCs were screened after treatment with benzene. Weighted gene correlation network analysis revealed that pathways in cancer, transcriptional misregulation in cancer, and hematopoietic cell lineage are vital pathways involved in benzene-induced toxicity in BM HSCs, whereas hematopoietic cell lineage and leukocyte transendothelial migration are critical pathways in PBSCs. Of note, there were 164 common DEGs in both HSCs, out of which 53 genes were co-regulated in both types of HSCs. Subsequent pathway analysis of these 53 genes indicated that the most relevant pathways involved neutrophil degranulation and CD93 localized in the core of the network of the 53 genes, which are known to regulate leukemia stem cell self-renewal and quiescence. Our results could enhance our understanding of HSC responses to benzene, facilitate the identification of potential molecular biomarkers and future studies on its mechanism of toxicity toward HSCs. •Exposure to benzene induced hematological damage in mice.•Genes and pathways involved in benzene-caused toxicity in bone marrow HSCs and peripheral blood HSCs were identified.•Fifty-three genes were co-regulated in both HSCs in mice and CD93 localized in the core of the network of 53 genes.