Single-cell analysis of the muscle stem cell hierarchy identifies heterotypic communication signals involved in skeletal muscle regeneration

Muscle stem cells (MuSCs) are an essential adult stem cell population with the capacity to self-renew and regenerate muscle tissue. Functionally heterogeneous subpopulations of MuSCs have been identified based on their expression of myogenic regulatory factors and surface markers. However, a unified organization of muscle stem and progenitor cells and their subpopulations remains unresolved. Here, we performed temporal analysis of skeletal muscle regeneration using single-cell RNA-sequencing (scRNA-seq) of myotoxin-injured adult mouse hindlimb muscles. We generated over 34,000 single-cell transcriptomes spanning four muscle regeneration time-points and identified 15 distinct cell types, including a heterogeneous population of MuSCs and progenitor cells. Our analysis provides a hierarchical map of myogenic cell populations and identifies stage-specific regulatory programs that govern their contributions to muscle regeneration. In this transcriptomic atlas, we observed cell type-specific regenerative dynamics, exemplified by waves of transient amplification and diversification of multiple immune cell types and, subsequently, myogenic cells. Unbiased trajectory inference organized the myogenic cell populations within the atlas into a continuum, consisting of a hierarchy of quiescent MuSCs, cycling progenitors, committed myoblasts, and terminally differentiated myocytes. This myogenic trajectory matched prior understanding and also revealed that MuSC stages are defined by synchronous changes in regulatory factors, cell cycle-associated, and surface receptor gene expression. Lastly, we analyzed the transcriptomic atlas to identify over 100 candidate heterotypic communication signals between myogenic and non-myogenic cell populations, including many involving the fibroblast growth factor (FGF), Notch, and Syndecan receptor families and their associated ligands. Syndecan receptors were implicated in a large fraction of these cell communication interactions and were observed to exhibit transcriptional heterogeneity within the myogenic continuum. Using multiparameter mass cytometry (CyTOF), we confirmed that cycling MuSCs exhibit diversified Syndecan-1/2 expression, which suggests that dynamic alterations in Syndecan signaling interactions may coordinate stage-specific myogenic cell fate regulation. This scRNA-seq reference atlas provides a resolved hierarchical organization of myogenic subpopulations as a resource to investigate cell-cell interactions that regulate myogenic stem and progenitor cell fates in muscle regeneration.