Elucidating the changes in the heterogeneity and function of radiation-induced cardiac macrophages using single-cell RNA sequencing

Abstract

Purpose: A mouse model of irradiation (IR)-induced heart injury was established to investigate the early changes in cardiac function after radiation and the role of cardiac macrophages in this process. Methods: Cardiac function was evaluated by heart-to-tibia ratio, lung-to-heart ratio and echocardiography. Immunofluorescence staining and flow cytometry analysis were used to evaluate the changes of macrophages in the heart. Immune cells from heart tissues were sorted by magnetic beads for single-cell RNA sequencing, and the subsets of macrophages were identified and analyzed. Trajectory analysis was used to explore the differentiation relationship of each macrophage subset. The differentially expressed genes (DEGs) were compared, and the related enriched pathways were identified. Single-cell regulatory network inference and clustering (SCENIC) analysis was performed to identify the potential transcription factors (TFs) which participated in this process. Results: Cardiac function temporarily decreased on Day 7 and returned to normal level on Day 35, accompanied by macrophages decreased and increased respectively. Then, we identified 7 clusters of macrophages by single-cell RNA sequencing and found two kinds of stage specific macrophages: senescence-associated macrophage (Cdkn1ahighC5ar1high) on Day 7 and interferon-associated macrophage (Ccr2highIsg15high) on Day 35. Moreover, we observed cardiac macrophages polarized over these two-time points based on M1/M2 and CCR2/major histocompatibility complex II (MHCII) expression. Finally, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses suggested that macrophages on Day 7 were characterized by an inflammatory senescent phenotype with enhanced chemotaxis and inflammatory factors, while macrophages on Day 35 showed enhanced phagocytosis with reduced inflammation, which was associated with interferon-related pathways. SCENIC analysis showed AP-1 family members were associated with IR-induced macrophages changes. Conclusion: We are the first study to characterize the diversity, features, and evolution of macrophages during the early stages in an IR-induced cardiac injury animal model.