Heterogeneity analysis of astrocytes following spinal cord injury at single‐cell resolution

Abstract Astrocytes play many important functions in response to spinal cord injury (SCI) in an activated manner, including clearance of necrotic tissue, formation of protective barrier, maintenance of microenvironment balance, interaction with immune cells, and formation of the glial scar. More and more studies have shown that the astrocytes are heterogeneous, such as inflammatory astrocyte 1 (A1) and neuroprotective astrocyte 2 (A2) types. However, the subtypes of astrocyte resulting from SCI have not been clearly defined. In this study, using single‐cell RNA sequencing, we constructed the transcriptomic profile of astrocytes from uninjured spinal cord tissue and injured tissue nearby the lesion epicenter at 0.5, 1, 3, 7, 14, 60, and 90 days after mouse hemisection spinal cord surgery. Our analysis uncovered six transcriptionally distinct astrocyte states, including Atp1b2 + , S100a4 + , Gpr84 + , C3 + /G0s2 + , GFAP + /Tm4sf1 + , and Gss + /Cryab + astrocytes. We used these new signatures combined with canonical astrocyte markers to determine the distribution of morphologically and physiologically distinct astrocyte population at injured sites by immunofluorescence staining. Then we identified the dynamic evolution process of each astrocyte subtype following SCI. Finally, we also revealed the evolution of highly expressed genes in these astrocyte subtypes at different phases of SCI. Together, we provided six astrocyte subtypes at single‐cell resolution following SCI. These data not only contribute to understand the heterogeneity of astrocytes during SCI but also help to find new astrocyte subtypes as a target for SCI repair.