We used our previously developed approaches to isolate human satellite cells from fresh human muscle biopsies (Garcia et al., 2018; Garcia et al., 2017; Xu et al., 2015). Recent studies using mouse (Machado et al., 2017; Der Vartanian et al., 2019; Scaramozza et al., 2019) and human (Charville et al., 2015; Cho and Doles, 2017) cells have provided evidence that satellite cells have heterogeneous patterns of gene expression. To discern satellite cell subpopulations in normal human muscle, we performed single cell RNA sequencing on highly purified CXCR4+/CD29+/CD56+ human satellite cells (Garcia et al., 2018) from eight vasti lateralis muscle of human subjects (20–83 years old) (Figure 1a and Figure 1—figure supplement 1a). A total of 68,108 cells were analyzed (See Figure 1—figure supplement 1b,c for QC data). Utilizing the Seurat single cell analysis package (Macosko et al., 2015; Butler et al., 2018) we found 17 distinct clusters of cells, as represented in 2D uniform manifold approximation and projection (UMAP) Figure 1b. While we found the distribution of each cluster to vary among different samples we didn’t find unique clusters associated with age or sex (Figure 1c). We confirmed clusters 0–8, 10, 12 and 15 to contain the vast majority of the cells (91%) and to consist of satellite cells by virtue of the expression of PAX7 and MYF5, while clusters 11, 13, 14, 16 and 9 contained small numbers of contaminating cells that were PAX7- and expressed either a mesenchymal/fibroblastic/smooth muscle and hematopoietic pattern, or a more differentiated muscle signature (Figure 1d,e and Figure 1—figure supplement 2a). Thus, 12 transcriptionally distinct satellite cell clusters were identified to be present across vastus lateralis biopsies.

Each cluster was found to have a unique transcriptomic fingerprint with heterogeneous gene expression, Figure 1e. Each cluster was characterized by the top differentially expressed genes, shown in Figure 1e and Supplementary file 1. Cluster 0 was characterized by upregulation of genes associated with the NOTCH pathway (DLK1; Guruharsha et al., 2012), G-protein signaling (GNAS; Weinstein et al., 2004) and satellite cell quiescence (CHRDL2; Charville et al., 2015). Cluster 1 up-regulated genes were associated with cellular stress response (DNAJA1; Stark et al., 2014, ATF3; Hai et al., 1999) and cellular development (EGR1; Zhang et al., 2018). Cluster 2 had high representation of genes involved in hormone signaling and inflammation (APOE; Liu et al., 2013a; Tzioras et al., 2019, SPARCL1; Hurley et al., 2015, PLCG2; Yu et al., 2005). Cluster 3 was found to have up-regulated genes associated with response to metal ions and oxidative stress (MT1X, MT1A; Si and Lang, 2018, SOX4; Pan et al., 2017). Cluster 4 was characterized by upregulation of genes associated with response to INF gamma (CCL2; Schroder et al., 2004), TNF signaling (SOD2; Yi et al., 2017, TAGLN; Yang et al., 2009). Cluster 5 was composed of cells expressing genes implicated in extra-cellular matrix (ECM), membrane receptor, focal adhesion (COLs, CAV1; Yeh et al., 2017, ANXA1; Sheikh and Solito, 2018) and down regulation of senescence. Cluster 6 consisted up-regulated genes associated with cellular development and mesenchymal cell development and differentiation (SIX1; Wu et al., 2014, SOX8; Schmidt et al., 2003) and myogenic commitment and differentiation (MYOD1, MYOG; Dumont et al., 2015). Based on the expression of MYOD1, cells in this cluster were considered activated and or progressing toward, differentiation. Cluster 7 was found to have cells with expression patterns of genes of positive regulation of cellular development, IGF1 (Mourkioti and Rosenthal, 2005; Schiaffino and Mammucari, 2011), KRT17 (Karantza, 2011) previously not described to be expressed by satellite cells, PTGDS (Moniot et al., 2014), and quiescence (CHODL; Machado et al., 2017). Cluster 8 up-regulated genes were associated with transcriptional repression (TXNIP; Elgort et al., 2010), and cellular stress response (KDM6B; Mallaney et al., 2019, GREM1; Simeckova et al., 2019). Cluster 10 contained cells expressing MX1, recently described to be present in a subpopulation of satellite cells (Scaramozza et al., 2019), and cell and genes associated with response to IFN gamma (STAT1; Qing and Stark, 2004, ISG15; Deng et al., 2015). Cluster 12 consisted of cells with expression patterns of genes of positive regulation of cell junction assembly (CAV1; Song et al., 2007, THY1; Kumar et al., 2016, S100A10; Lee et al., 2004), collagen catabolic processes, extracellular matrix assembly/disassembly (TIMP1, MMP2; Arpino et al., 2015), response to TGF-beta (CAV1; Gvaramia et al., 2013, COLs; Biernacka et al., 2011, CLDN5; Shen et al., 2011), blood vessel development (GPC3; Ng et al., 2009, HSPG2; Martinez et al., 2018, CAV1; Yu et al., 2006, COLs; Marchand et al., 2018, THY1; Inoue et al., 2016, MMP2; Brooks et al., 1996, TIMP1; Arpino et al., 2015, B2M; Smith et al., 2015). Cluster 15 consisted of a small number of cells expressing genes associated with DNA replication and cell cycle control (TOP2A, BIRC5, MKI67) suggesting satellite cell activation. We also evaluated expression of PAX3, which has been shown to be enriched in a subset of Pax7+ mouse satellite cells (Der Vartanian et al., 2019; Scaramozza et al., 2019). We found low detectable PAX3 expression present all the myogenic clusters however cells expressing PAX3 did not form a unique cluster (Figure 1—figure supplement 2a).

The differentiated muscle cluster (9), hematopoietic cluster (14) and fibrogenic/mesenchymal clusters (13,16) were confirmed to have expression profiles with up-regulated genes consistent with their identities respectively (Figure 1e and Supplementary file 1). Cluster 11 was made of cells expressing high levels of DNA damage expressing genes (CYCS; Pal et al., 2010, H2AFZ; Flint et al., 2007, DDX5; Nicol et al., 2013) and ribosomal genes suggesting low quality cells. We also analyzed human satellite cells from the rectus femoris (a different muscle from the quadriceps group) of an 84 year old male (Figure 1—figure supplement 2b). From this individual sample, a total of 4,575 cells were analyzed (Figure 1—figure supplement 1c). Similar to the analysis of eight combined vasti, we confirmed clusters to consist of satellite cells by the expression of PAX7 and MYF5 (Clusters 0–3), while clusters 4 and 5 contained a small number of cells that were PAX7- and expressed either a mesenchymal/fibroblastic/endothelial pattern, or a more differentiated muscle signature, respectively (Figure 1—figure supplement 2c). Analogous cluster markers to the vasti samples were found in this rectus femoris sample (i.e. clusters expressing genes associated with NOTCH signaling (0), cellular stress (1), response to IFN gamma (cluster 2-CCL2, MYC, ICAM1, IRF1) (Schroder et al., 2004), ECM and cell adhesion (cluster 3-CAV1, COLs, VCAM1) (Figure 1—figure supplement 2d and Supplementary file 2). Taken together, these scRNAseq results show that satellite cells in human uninjured muscle, can be separated into transcriptionally distinct subpopulations (Figure 1f).

We selected highly expressed markers to validate the predicted human satellite cell subpopulations at the protein level, by immunofluorescence. We mapped five of the top differentially expressed genes (CYCS, DLK1, ICAM1 and VCAM1) onto the 2D UMAP space, (Figure 2a). These genes are heterogeneously expressed across the eight vasti lateralis (Figure 1e), the rectus femoris sample, two pooled recti abdominis and two pooled pectoralis major samples (Figure 2a and Figure 2—figure supplement 1a,b,c). We then assessed proportions of cells with each marker detectable in human muscles by both immunofluorescence staining and flow cytometry (Figure 2b and Figure 2—figure supplement 1d and Figure 2—figure supplement 2a–e). Muscle sections from human tissue were stained with PAX7 and LAMININ to identify PAX7 expressing satellite cells. Sections were also stained with each cluster marker, as shown in Figure 2b and Figure 2—figure supplement 2e. With all five markers, we found PAX7 satellite cells in individual muscles that were either negative or positive for expression of each respective protein, confirming in vivo heterogeneity of sublaminar satellite cells. We also confirmed the detection of the surface markers ICAM1, DLK1 and VCAM1 on subsets of satellite cells isolated from human muscle using flow cytometry (Figure 2—figure supplement 2a–e). This utilization of flow cytometry for validation of surface markers demonstrates the feasibility of separating and isolating heterogeneous subpopulations of human satellite cells for downstream use and experimentation. The proportions of satellite cells expressing each cluster marker were quantified and compared among data from immunofluorescence staining, flow cytometry, and the number of expressing cells within the in silico data (Figure 2c). We discovered the proportion of PAX7+ satellite cells expressing DLK1, ICAM1, CYCS, and VCAM1 by immunofluorescence staining to be 35.0 ± 1.7%, 38.8 ± 1.9%, 49.7 ± 1.7%, and 10.6 ± 2.8%, respectively. In addition, we also found heterogeneous co-expression of those markers (DLK1+/ICAM1+, 28%; ICAM1+/VCAM1+, 21% and DLK1+/VCAM1+, 16%) in PAX7+ cells (Figure 2—figure supplement 2d, e). The portion of CXCR4+/CD29+/CD56+ cells expressing ICAM1, DLK1, and VCAM1 measured by flow cytometry was found to be 63.9 ± 13.8%, 37.0 ± 0.6%, and 13.5 ± 4.3%, respectively. Our results confirmed that differentially expressed genes (DLK1, ICAM1, CYCS, and VCAM1) identified via scRNAseq are expressed heterogeneously by human satellite cells in vivo.

We next evaluated the transcriptional signatures associated with stemness, quiescence and activation in the previously identified satellite cell clusters 0–8, 10, 12 and 15, as well as differentiated myogenic cells from cluster 9 (Figure 3a). We evaluated the expression of cell cycle markers previously described as part of the quiescent stem cell gene signature of satellite cells (Cheung and Rando, 2013). We found one cluster (15) composed of fully activated satellite cells (0.52%) with increased expression of G2/M/S phase markers such as ANLN, BIRC5, CCNA2, CCNB1 or CCNE2 while in the rest of the satellite cell clusters these were downregulated (Figure 3b). This was also confirmed by the cell cycle scoring vignette from Seurat and the increased RNA counts in cluster 15 (Figure 3—figure supplement 1a,b). Moreover, Ki67 was only detected in cluster 15 providing additional evidence that the majority of satellite cells were not activated. We noticed an increase of the RNA count in cluster eight despite low levels of activation-associated cell cycle genes although this cluster also showed increased level expression of MEF2C compared to the rest of satellite cell clusters (Figure 3c). As expected, the cluster consisted of differentiated myocytes and showed characteristics of non-proliferative cells (e.g. low expression of G2/M/S phase markers (Figure 3b), low RNA count (Figure 3—figure supplement 1b).

We also evaluated previously used markers of stemness (KLF4, MYC) (Takahashi et al., 2007), quiescence (NDRG2, DAG1 and CHRNA1 (Charville et al., 2015), SPRY1 (Bigot et al., 2015; Chakkalakal et al., 2012), Shea et al., 2010), the NOTCH pathway targets HEY1 (Jiang et al., 2014), EGR1 (Min et al., 2008), HES1 (Jiang et al., 2014), CD82 (Alexander et al., 2016), and activation (MEF2C; Liu et al., 2014) as shown in Figure 3c. Each satellite cell cluster demonstrated a mixed pattern of expression of those markers. We found differential expression among satellite cell clusters including KLF4, MYC, and EGR1 (cluster 1); HEY1 (cluster 12) and HES1 (cluster 6, 9, 12); CD82 (cluster 5, 12, 15)); DAG1 (6, 15, 4) and CHRNA1 (cluster 6, 12). Cluster six in which MYOD1 and MYOG were differentially expressed also had HES1, DAG1 and CHRNA1 significantly upregulated. SPRY1 was expressed in all of the satellite clusters at varying levels, while MEF2C was differentially upregulated in the myocyte cluster 9. In agreement with in silico analysis, immunofluorescence of human muscle tissue revealed that SPRY1 was detectable in a subset (~33%) of satellite cells in vivo (Figure 3d). We also found SOX8, a previously described satellite cell marker (Schmidt et al., 2003), to be associated with activation in human satellite cells (Figure 3c). Independent analysis of the rectus femoris resulted in similar findings (Figure 3—figure supplement 1c). Taken together, while we observed evidence of variable profiles consistent with quiescence, priming, or early activation, and notwithstanding a presumed effect of the purification and preparation process (i.e, cluster 6, 8, 15) the majority of satellite cells analyzed are best characterized by pre-activation, relatively quiescent states.

In order to estimate the lineage relationships between the satellite cell clusters, we performed pseudotime analysis of our single cell data on all myogenic clusters utilizing the R package Monocle (Qiu et al., 2017a; Qiu et al., 2017b; Trapnell et al., 2014; Figure 3e). We found the cells ordered from proximal to distal in an arrangement compatible with quiescent satellite cells transitioning towards myogenic differentiation with several branching points (Figure 3e). Cells from cluster 9, which express genes of myogenic commitment and terminal differentiation were located distal to the satellite cell clusters at the end point of pseudotime while satellite cell clusters were distributed more proximally along sub-branches. The gene expression pattern across pseudotime suggests that two main populations of satellite cells diverge in accordance with expression of PAX7, SPRY1, HEY1, DLK1 and CAV1 (Figure 3f, dashed line). Analysis of the rectus femoris was confirmatory, with satellite cells diverging in multiple branches while the most distal branch consisted of differentiated myogenic cells. The branch represented by the solid line in Figure 3—figure supplement 1d, e) was driven by the expression of PAX7, MYF5, HEY1 and CAV1 (Figure 3—figure supplement 1e). Taken together, the pseudotime analysis shows that within resting human muscle, satellite cells can be ordered in various states of transition from stem-like cell to more differentiated progenitors. Transcription factors associated with more stem-like and quiescent, non-activated satellite cell states early in pseudotime (PAX7, SPRY1, HEY1) were associated with surface proteins DLK1 and CAV1, indicating a potential avenue for a physical sorting strategy.

We were interested in separating functionally distinct satellite cell subpopulations from normal adult muscle. We therefore decided to analyze the transcriptome data based on previously established markers of satellite cell quiescence: Sprouty1 (SPRY1) (Shea et al., 2010) and Hairy/enhancer-of-split related with YRPW motif protein 1 (HEY1) (Fukada et al., 2011). To do this, we performed an in silico grouping of the satellite cell populations in the adult sample dataset (clusters 0–8, 10, 12 and 15) excluding the rare contaminating and differentiated non-satellite cell populations. Cells were then sorted into two groups in silico by expression of SPRY1 and HEY1: SPRY1/HEY1 high expressing cells (SPRY1^hi/HEY1^hi) and SPRY1/HEY1 low/negative expressing cells (SPRY1^low/neg/HEY1^low/neg) (Figure 3g). The two groups were analyzed for differentially expressed genes with greater than 1.5-fold upregulation. In agreement with known MRF expression in quiescent and activated satellite cells, we found that both PAX7 and MYF5 were upregulated in the SPRY1^hi/HEY1^hi group, while MYOD1, MYOG, and MYF6 were upregulated in the SPRY1^low/neg/HEY1^low/neg group (Figure 3h). Indeed, several other genes associated with quiescence and activation were associated with the SPRY1^hi/HEY1^hi group (Fibronectin 1, quiescence marker CHRNA1, DAG1, DLK1 and CHRDL2) and the SPRY1^low/neg/HEY1^low/neg group (CYCS, ElF4A1, MTIX, PDK4, SOD2 and CYR61) respectively (Supplementary file 3). Moreover, caveolae scaffolding protein, Caveolin-1 (CAV1), which was associated with early pseudotime satellite cells, showed a larger probability density in the SPRY1^hi/HEY1^hi cells. CAV1 has been shown previously to be expressed by mouse satellite cells (Kann and Krauss, 2019; Gnocchi et al., 2009) and expression is reported to be lost with activation (Volonte et al., 2005). This prompted us to perform in silico sorting of satellite cells by expression of CAV1: CAV1 high expressing cells (CAV1^high) vs CAV1 low expressing cells (CAV1^low) (Figure 4a). The average expression of quiescence marker SPRY1, HEY1, CD82, DAG1 and CHRNA1 was increased in CAV1^high expressing satellite cells while the expression KLF4 and MYC were decreased (Figure 4b) further suggesting that CAV1 may be associated with satellite cells in a more quiescent state. The two groups were analyzed for differentially expressed genes. Differentially up- and down-regulated genes were subjected to Gene Ontology (GO) analysis and Pathway analysis (Figure 4c,d and Supplementary file 4). The GO analysis revealed that up-regulated genes were predominantly associated with extracellular matrix organization and cell-cell adhesion regulation while down-regulated genes were associated with myogenesis and stress responses. Pathway analysis showed that up-regulated genes in the CAV1^high cells were associated with extracellular matrix, membrane receptors, VEGF signaling and focal adhesion. Down-regulated genes were associated with TGF-beta signaling. These gene expression analyses together suggest that CAV1 expression correlates with both quiescence and cell adhesion characteristics in human satellite cells.

CAV1 was detected in 55% of cells, however expression was significantly upregulated in clusters 4, 5, 12 and 15 across the eight vasti lateralis and rectus femoris, rectus abdominis and pectoralis major samples (Figure 5a). CAV1 protein expression was also heterogeneous when assayed in muscle sections and by flow cytometry of satellite cells (Figure 5b and c). We detected CAV1 expression on 10.65 ± 3.6% (n = 19) of satellite cells and 0.7 ± 0.3% of the whole cell population by flow cytometry (Figure 5c). The less frequent detection of CAV1 by flow cytometry of live cells reflects either the various intracellular localization possibilities of CAV1 (Boscher and Nabi, 2012), which are detected by immunofluorescence of sections, whereas flow cytometry will only detect CAV1 present on the cell surface (Wu and Terrian, 2002). Alternatively, CAV1 protein may be degraded rapidly upon SC isolation from the niche.

To physically separate CAV1 human satellite cells, CXCR4+/CD29+/CD56+ cells were sorted based on CAV1 expression (Figure 5d), marking human satellite cells with differential surface CAV1 expression (CAV1+ and CAV-). (Figure 5f) Upon back-gating analysis of the CAV1+ and CAV1- populations (Figure 5e) we found overlap in regard to cell size, cell granularity, and CXCR4, CD29 and CD56 expression, indicating that separation of these two populations is not attributable to other commonly assessed satellite cell characteristics. CAV1+ sorted cells expressed significantly higher transcriptional levels of CAV1 compared to CAV1- (Figure 5f) confirming that flow cytometry captures CAV1^high expressing satellite cells identified by scRNA sequencing. Moreover, both CAV1+ and CAV1- sorted satellite cells expressed PAX7 by immunofluorescence and qPCR (Figure 5—figure supplement 1a,b), ruling out the possibility that CAV1- cells are non-satellite cell contaminants. This possibility is also contrary to the transcriptome data, which show that a large proportion of satellite cells are CAV1- or CAV1^low.

The ability to separate CAV1^high expressing cells permitted us to test phenotypic and functional characteristics. CAV1+ and CAV1- satellite cells were sorted and cultured in growth-promoting conditions. The morphological phenotype differed significantly, with CAV1+ cells adopting a round shape compared to the more spindle and elongated morphology of CAV1- satellite cells (Figure 5g). We next used time lapse video microscopy to evaluate time to first division (Kuang et al., 2007; Marti et al., 2013; Siegel et al., 2009). Three separate experiments using satellite cells from four unique human muscles were performed by comparing CAV1+ and CAV1- satellite cells from the same muscle. Cells were sorted and then placed into culture and live stained with CellTracker Green for live visualization for 6 days. While none of the CAV1+ Hu-MuSC divided during this time, 5.4% (range 2.5–11%) of CAV1- Hu-MuSC completed their first cellular division in the first six days in vitro. The average time for CAV1 negative Hu-MuSC to complete their first division was 4.4 days (range 2.8–5.7 days). The time to first division was significantly longer in CAV1+ cells compared to CAV1- (Figure 5h). In support of this finding of slower division, we found that Ki67 expression is expressed in fewer CAV1+ Hu-MuSCs at day 3 of culture (Figure 5i). Finally, CAV1+ and CAV1- satellite cells were cultured in growth media for 3 days and MYOD expression was evaluated. CAV1+ Hu-MuSC expressed significantly less MYOD than their CAV- counterparts (Figure 5j,k). These in vitro assays are consistent with the previously discussed transcriptome data and indicate that CAV1+ cells differ from other human satellite cells in canonical assays of satellite cell activation and proliferation, collectively characterized by resistance to activation.

The capacity to separate a subpopulation of satellite cells based on CAV1 expression enabled us to assay in vivo functional heterogeneity. From single muscle samples, 500 CAV1- and CAV1+ satellite cells were transplanted into the tibialis anterior (TA) of pre-irradiated NOD scid gamma (NSG) mice (Garcia et al., 2018; Xu et al., 2015). We (Garcia et al., 2018; Xu et al., 2015) and others have (Arpke and Kyba, 2016; Gayraud-Morel et al., 2012) previously demonstrated that the use of small number of cells are effective for transplantation studies. Mice were sacrificed 5 weeks later and evaluated for engraftment and myogenic differentiation using human specific DYSTROPHIN antibody (Figure 6a). The engraftment of human muscle was quantified by counting the maximum number of human specific DYSTROPHIN fibers per cross-section in each experimental group. We found that transplantation with CAV1+ human satellite cells led to robust engraftment in contrast to transplantation with equal numbers of CAV1- cells (Figure 6b). The CAV1+ transplants resulted in 4-fold higher engraftment over the CAV1- group (69.3 ± 18.4* vs 14.4 ± 6.3* human fibers respectively (p=0.013)), and the efficiency of human fibers corresponded to roughly 1 fiber per seven satellite cells transplanted. Repopulation of the satellite cell niche by human PAX7 cells was also significantly increased in the CAV1+ transplants when evaluated by counting human-derived PAX7 sublaminar cells on cross sections shown in Figure 6c, and quantified in Figure 6d. Finally, using the re-isolation assay we previously developed (Garcia et al., 2018) we were able to re-purify human satellite cells from mice transplanted with CAV1+ but not CAV1- cells (Figure 6e). In conclusion, human CAV1+ satellite cells engraft, differentiate, repopulate the satellite cell niche and retain satellite cell phenotypes efficiently and to a greater extent than the more numerous CAV1- satellite cells.

This study was conducted under the approval of the Institutional Review Board at The University of California San Francisco (UCSF). Biopsies were obtained from individuals undergoing surgery at UCSF. Written informed consent was obtained from all subjects. All types of muscle used for each experiment are listed in Supplementary file 5.

All mice were bred and housed in a pathogen-free facility at UCSF. All procedures were approved and performed in accordance with the UCSF Institutional Animal Care and Use Committee. All experiments were unblinded and performed in 8–12 week-old NSG. Mice were randomized to all experimental groups by sex and littermates and were pretreated with 18 gamma (Gy) on the day before transplantation. Hu-MuSCs were injected along with 50 µl 0.5% bupivacaine directly into the TA muscle of one leg as indicated (Garcia et al., 2017). Detailed information can be found in the Supplemental Experimental Procedures section.

Freshly harvested human muscle was either immediately digested or stored in DMEM with 30% FBS at 4°C overnight. Muscle was digested, erythrocytes were lysed, and hematopoietic and endothelial cells were depleted with magnetic column depletion (Miltenyi Biotech). Viable cells were depleted for CD31, CD34, and CD45 expressing cells. Cells that remained after depletion were sorted for CXCR4+/CD29+/CD56+ and collected for further experimentation (Garcia et al., 2018; Garcia et al., 2017).

To capture individual cells, we utilized the Chromium Single Cell 3' Reagent Version one and Version 3 Kit from 10X Genomics (Zheng et al., 2017). For all samples (vasti, rectus femoris, recti abdominis and pectoralis major) 18,000–30,000 satellite cells isolated as in Garcia et al. (2018) were loaded onto one well of the 10X chip to produce Gel Bead-in-Emulsions (GEMs). GEMs underwent reverse transcription to barcode RNA before cleanup and cDNA amplification. Libraries were prepared with the Chromium Single Cell 3' Reagent Version 1and 3 Kit (see Figure 1—figure supplement 1b,c). Each sample was sequenced on 1 lane of the HiSeq2500 (Illumina) run in Rapid Run Mode with paired-end sequencing parameters or 1 lane of the NovaSeq 6000 S4. The estimated number of cells, mean reads per cell, median genes per cells, median UMI (Unique Molecular Identifier) counts per cells as well as other quality control information are summarized in Figure 1—figure supplement 1b,c. Gene-barcoded matrices were analyzed with the R package Seurat v3.1, (Satija et al., 2015; R Development Core Team, 2014; Zheng et al., 2017). Cells with fewer than 500 genes, greater than 6000 genes and genes expressed in fewer than 5 cells were not included in the downstream analyses. We also filtered cells that had more than 10% mitochondrial counts. In all samples UMI counts (or RNA counts) were normalized with NormalizeData using default settings. The FindVariableFeatures function was used to determine subset of feature that exhibit high cell-to-cell variation in each dataset based on a variance stabilizing transformation (‘vst’). We used the default setting returning 2000 feature per dataset. These were used for downstream analysis. In the case of the merged data analysis samples were combined utilizing the FindIntegrationAnchors function with the ‘dimensionality’ set at 30. Then, we ran these ‘anchors’ to the IntegratData function for batch correction for all cells enabling them to be jointly analyzed. The resulting outputs were scaled and UMI counts and mitochondrial contamination regressed out with the ScaleData function. We didn’t regress out heterogeneity associated with cell cycle stage since it is an important factor in determining the state of quiescence of our sorted human satellite cells. PCA was performed with RunPCA, and significant PCs determined based on the Scree plot utilizing the function PCElbowPlot. The resolution parameter in FindClusters was adjusted to 0.5. Clusters were visualized by UMAP with Seurat’s RunUMAP function. Differentially expressed genes were determined with the FindAllMarkers function. We performed differential gene-expression utilizing Seurat v3’s FindMarkers function with the Model-based Analysis of Single-Cell Transctiptomics (MAST) that treates cellular detection as a covariate to calculate adjusted p values for multiple comparisons. Lists of differentially expressed genes for individual analyses are provided in Supplementary file 1 and 2. The cell-cycle scoring vignette from Seurat v3 was used to calculate cell cycle phase scores for each cell based on its expression of G2/M and S phase markers. Cells that didn’t express G2/M and S phase markers were scored as not cycling cells in G0/G1 phase. Scores were assigned using the CellCycleScoring function and visualized in a barplot for each cluster. The in silico FACS were done using the subset function for each gene of interest. The differential expression was performed as described earlier.

We utilized Monocle 2.12.0 to order cells in pseudotime based on their transcriptomic similarity (Qiu et al., 2017b). Variable genes from Seurat analysis were used as input and clusters were projected onto the minimum spanning tree after ordering. For computing power purposes, the combined vasti object was downsampled to 500 cells per cluster. Gene expression patterns were plotted with plot_genes_branched_pseudotime, and plot_multiple_branches_pseudotime function.

Cells were plated immediately after sorting on Matrigel coated chamber slides. 3 hr after plating cells were stained for PAX7 (DSHB). CD56-CD29-CXCR4- cells were used as controls. MYOD and Ki67 protein expression were assessed 3 days after plating. Cryosection slides or sorted human satellite cells were fixed with 4% PFA at room temperature for 10 min, washed in PBST (Phosphate Buffered Saline Tween20 0.1%), permeabilized with 0.1%Triton-100X (Sigma-Aldrich) and then blocked with protein-free serum block (DAKO) or 2% goat serum and incubated at room temperature overnight with primary antibodies (Supplemental Experimental Procedures and Key Resources Table). After PBST wash the corresponding secondary antibodies were applied for 1 hr at room temperature. Sections were mounted with VECTASHIELD mounting medium with DAPI (Vector Laboratories) and all samples were examined using a Leica upright or DMi8 Leica microscope.

We used the HEK293 cell line, ATCC Cat# PTA-4488, as a control in Figure 5—figure supplement 1. The cell line is tested periodically for mycoplasma and is negative to date.

Tissues were collected in RLT buffer (Qiagen), total RNA was isolated using the RNAeasy isolation kit (Qiagen). RNA was transcribed into cDNA with High-Capacity cDNA Reverse Transcription kit (ThermoFisher Scientific). cDNA was then pre-amplified with GE PreAmp Master Mix (Fluidigm Inc). Real-time quantitative PCR was performed in triplicated with either Taqman Universal PCR Master Mix (Life Technologies) on either a Viia7 thermocycler (Life Technologies). Taqman primers are listed in the Key Resources Table. Beta actin was used for normalization as endogenous control.

Hu-MuSCs were sorted and plated at a density of 500 to 1000 cells per well on a 48 glass well plate (Mattek) precoated with Matrigel (Corning). Hu-MuSCs were grown in Growth media: DMEM high glucose phenol free media, 20% FBS and 1% pen/strep (Gibco). The following day, satellite cells were incubated with 10 mM final concentration of CellTracker Green CMFDA dye (Thermofisher Scientific) in phenol free media for 45 min at 37°C to track cell division. After a media wash, fresh growth media was added for subsequent time lapse experiments, cells were imaged using Zeiss Confocal Microscope. Images of Hu-MuSC were taken every 15 min for 6 days. Images and videos were analyzed using Zeiss Zen microscope software. Statistical analysis was done using GraphPad Prism.

Normality of the data was checked utilizing the Shapiro-Wilk normality test in GraphPad Prism. Means between or across groups were compared using two-tailed t-tests for experiments involving two groups, or one-way ANOVA with post hoc Tukey multiple comparisons when comparisons were made across three or more groups to determine significance (p<0.05) between test conditions and controls, and multiple groups. Multivariate regression was utilized as indicated for comparing satellite cell yield per gram controlling for age, gender, and muscle type. All human muscle samples collected over the past one year and processed within 12 hr after biopsy were used for data analyses in Figure 1. At least three mice per group were used for all transplantation experiments. At least three biological replicates (three different muscle source) for each experiment were performed unless otherwise noted, with exact n values listed in each figure legend. For CAV1+/- statistical analysis paired or unpaired t-tests were used. All error bars are depicted as s.e.m. p-values are indicated with asterisks (*p<0.05, **p<0.01, ***p<0.001).

Freshly harvested human muscle was either immediately digested or stored in DMEM with 30% FBS at 4°C. Muscle was trimmed of excess fat, tendon, connective tissue, and fascia and mechanically minced. The tissue was then digested in 1 mg/ml collagenase XI (Sigma-Aldrich) in Dulbecco's Modified Eagle Medium (DMEM) with high glucose, 10% FBS and 1% Penicillin/Streptomycin at 37°C for 70 min with intermittent manual needle trituration, performed slowly with an 18-gauge needle. Digests were washed with PBS and further digested with 0.25% trypsin at 37°C for 12–15 min. Suspensions were passed through 40 µm nylon mesh, erythrocytes were lysed with ACK lysing buffer (ThermoFisher) for 5–7 min on ice, and washed with PBS. Magnetic column depletion of hematopoietic and endothelial cells was performed after cells were stained with anti-CD45 and anti-CD31 magnetic beads (Miltenyi Biotec). This step has the added benefit of removing small fiber fragments and facial tissue, which are a cause of high background on the flow cytometer. Unbound cells were washed and stained with anti-CD29-488 or 647 (eBioscience), anti-CD31-450 (eBioscience), anti-CD34-450 (eBbioscience), anti-CD45-450 (eBbioscience), anti-CD56-APC-vio-770 (Miltenyi Biotec), and anti-CXCR4-PE or APC (eBbioscience) (Note for the reisolation of Hu-MuSCs from transplanted mice: mouse muscle was processed as stated for human muscle, stained with the following antibodies: anti-human CD29-488 or 647 (Ebioscience), anti-human CD31-450 (eBbioscience), anti-human CD45-450 (eBbioscience), anti-human CD56-APC-vio-770 (Miltenyi Biotec), anti-human CXCR4-PE or APC (eBioscience), anti-mouse CD31-450 (eBioscience), anti-mouse CD45-450 (eBbioscience), and anti-mouse Sca1-450 (eBbioscience)). Cells were washed and resuspended in flow cytometry buffer with 1:1000 sytox blue (Life Technologies). Flow cytometry antibodies listed in the Key Resources Table. Flow cytometry analysis and cell sorting were performed at the University of California San Francisco Flow Cytometry Core with the BD FACSAria2 operated using FACSDiva software. Viable cells were depleted for CD31, CD34, and CD45 expressing cells. Cells that remained after depletion were sorted for CXCR4+/CD29+/CD56+ and collected for further experimentation. We have previously published FMO controls for CD56 and CD29 use in Hu-MuSC isolation (Xu et al., 2015). Cells were sorted in 20% FBS in DMEM supplemented with 10 μM Rho-associated protein kinase inhibitor (ROCKi) (Y-27632 2HCl, Selleck Chemicals). See (Garcia et al., 2017; Xu et al., 2015) for details of the authors’ prior muscle digestion and Hu-MuSC isolation protocol. Flow cytometry isolations were analyzed with FACSDiva and FlowJo software. Satellite cell subpopulations were analyzed and/or sorted with the following antibodies: anti-human CAV1-PE (abcam), ICAM1-PE (abcam), VCAM1 (eBioscience) and DLK1 (R and Dsytems).

All glass slides were removed from −80°C and warmed at room temperature for 10 min. For human DYSTROPHIN immunostaining, sections were fixed in 4% PFA for 10 min at room temperature and then washed in PBST (PBS with 0.1% Tween-20 (Sigma-Aldrich). The sections were blocked with 10% goat serum in PBS for 10 min at room temperature. The sections were then incubated overnight at room temperature with mouse monoclonal anti-human DYSTROPHIN (1:10 DSHB), human specificity of which was previously confirmed (Xu et al., 2015). The sections were then washed in PBST followed by 1 hr of incubation at room temperature with Alexa Fluor 594 goat anti-mouse IgG (1:500 Thermo) in 10% normal goat serum in PBS. Sections were mounted with VECTASHIELD mounting medium with DAPI (Vector Laboratories) and all samples were examined using a Leica upright microscope. Human-derived fibers (e.g. hDYSTROPHIN positive) were quantified by counting the number of positively stained fibers in the section with the most positive fibers after analyzing sections along the length of the muscle as has been previously reported (Rozkalne et al., 2014; Xu et al., 2015). For all other immunostainings, the slides were fixed in 4% PFA at room temperature for 10 min, washed in PBST, and then blocked with protein-free serum block (DAKO) and incubated at room temperature overnight with the following primary antibodies: mouse monoclonal IgG1 anti-PAX7 (1:10 DSHB), rabbit polyclonal anti-LAMININ (1:250 Sigma-Aldrich), mouse monoclonal IgG2b anti-human SPECTRIN (Leica Microsystems), mouse monoclonal IgG2b anti-human LAMIN A/C (Vector Laboratories), DLK1(1:50 abcam), ICAM1(1:50 Santa Cruz Biotechnology), CYCS (1:100 LifeSpan Bioscience), VCAM1 (1:75 abcam), CAV1 (1:500 abcam) and Sprouty1 (1:50). After PBST wash the following corresponding secondary antibodies were applied for 1 hr at room temperature: FITC donkey anti-mouse (1:500 Jackson Immunology), Cy3 goat anti-mouse (1:500 Jackson Immunology), Cy5 donkey anti-mouse (1:500 Jackson Immunology), Cy5 donkey anti-rabbit (1:300 Jackson Immunology), Alexa Fluor 488 goat anti-mouse IgG1 (1:500 Thermo), Alexa Fluor 594 goat anti-mouse IgG1 (1:500 Thermo), Alexa Fluor 488 goat anti-mouse IgG2b (1:500 Thermo), Alexa Fluor 594 goat anti-mouse IgG2b (1:500 Thermo). Sections were mounted with VECTASHIELD mounting medium with DAPI (Vector Laboratories) and all samples were examined using a Leica upright microscope.

Single cell gene expression data have been deposited and can be found here: https://doi.org/10.7272/Q65X273X (Pomerantz and Barruet, 2020). Detailed scripts for each analysis are in Source code 1–5.