Modeling the ACVR1R206H mutation in human skeletal muscle stem cells

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Abnormalities in skeletal muscle repair can lead to poor function and complications such as scarring or heterotopic ossification (HO). Here, we use fibrodysplasia ossificans progressiva (FOP), a disease of progressive HO caused by ACVR1^R206H (Activin receptor type-1 receptor) mutation, to elucidate how ACVR1 affects skeletal muscle repair. Rare and unique primary FOP human muscle stem cells (Hu-MuSCs) isolated from cadaveric skeletal muscle demonstrated increased extracellular matric (ECM) marker expression, showed skeletal muscle-specific impaired engraftment and regeneration ability. Human induced pluripotent stem cell (iPSC)-derived muscle stem/progenitor cells (iMPCs) single-cell transcriptome analyses from FOP also revealed unusually increased ECM and osteogenic marker expression compared to control iMPCs. These results show that iMPCs can recapitulate many aspects of Hu-MuSCs for detailed in vitro study; that ACVR1 is a key regulator of Hu-MuSC function and skeletal muscle repair; and that ACVR1 activation in iMPCs or Hu-MuSCs may contribute to HO by changing the local tissue environment.

Editor's evaluation

The manuscript by Barruet et al., investigates an interesting and rare skeletal muscle dystrophy (FOP). They use both primary and induced (iPSC) muscle stem cells to determine how regeneration/engraftment is affected in this condition. The authors use a blend of histology and transcriptional approaches to determine ECM remodeling and myogenic capacity from both cell/tissue lines. The experiments are well conducted and use strong approaches and statistical measures to test their hypothesis. Overall, this is a quality manuscript on a rare muscle disease that will establish a novel model to study FOP and initial data to elucidate the molecular pathology of HO.

https://doi.org/10.7554/eLife.66107.sa0

Introduction

Human diseases of skeletal muscle are major medical problems. Aberrant repair after muscle injury can results in scarring, fibrosis, lipid infiltration, or heterotopic ossification (HO; bone formation in an inappropriate site). Although HO is a relatively uncommon complication, it can form after trauma, burns, spinal cord injuries, and surgical procedures such as hip arthroplasty (Agarwal et al., 2016; Matsuo et al., 2019; Meyers et al., 2019; Sullivan et al., 2013). Recently, HO has been reported in patients with severe COVID-19, suggesting a critical role for inflammation and muscle damage (Aziz et al., 2021).

To better understand the link between skeletal muscle and HO, we studied fibrodysplasia ossificans progressiva (FOP) as a prototypical condition. FOP is a congenital disease of abnormal skeletal muscle regeneration and severe HO, in which the majority of patients have the classical ACVR1^R206H (c.617> A) mutation (Shore Eileen et al., 2006) that causes hyperactivation of the BMP-SMAD signaling (Billings et al., 2008) and aberrant responses to Activin A (Hatsell et al., 2015). While previous studies suggested that multiple cell types including endothelial cells (Barruet et al., 2016; Medici et al., 2010), mesenchymal progenitors (Culbert et al., 2014; Dey et al., 2016), pericytes (Cai et al., 2015), and more recently mast cells (Brennan et al., 2018; Convente et al., 2018), macrophages (Barruet et al., 2018; Matsuo et al., 2021), and FAPs (Lees-Shepard et al., 2018) can indirectly or directly contribute to HO in FOP, how the ACVR1^R206H mutation impairs human muscle repair and muscle stem cell function has yet to be determined.

Disrupted signaling of bone morphogenetic proteins (BMPs), originally identified by their ability to induce bone formation when injected into muscle (Urist, 1965), changes muscle homeostasis (Ono et al., 2011) by controlling the proliferation and differentiation of satellite cells (SCs) (Ono et al., 2011; Stantzou et al., 2017). SCs marked by PAX7 (Paired Box 7) (Seale et al., 2000), are a prerequisite for skeletal muscle regeneration (Mauro, 1961) and are thought to be the main human muscle stem cells (Hu-MuSCs) of postnatal skeletal muscle. Upon injury, activated SCs give rise to myoblasts, which form new myofibers or fuse to existing muscle fibers to repair muscle damage (Kuang et al., 2007). A subset of SCs does not differentiate and serves to replenish the SC pool. The BMP pathway is critical for maintaining PAX7 expression in primary MuSCs and for preventing commitment to myogenic differentiation (Friedrichs et al., 2011). Abrogation of BMP signaling in SCs slowed myofiber growth (Stantzou et al., 2017), and increased BMP4 levels in Duchenne’s muscular dystrophy (DMD) can exacerbate the disease (Shi et al., 2013). However, how the BMP pathway regulates muscle repair or SC function in conditions of HO such as FOP remains unclear and is a major knowledge gap.

Since primary tissues can only and safely be obtained from deceased FOP patients because of the risks of heterotopic ossification at the surgical site, patient-derived human-induced pluripotent stem cells (hiPSCs) provide a strategy to create FOP cells for in vitro and mouse xeno-transplant studies. Interestingly, BMP signaling is a key pathway that can be manipulated to make Hu-MuSCs-like cells from hiPSCs (iMPCs) (Chal et al., 2016; Chal et al., 2018; Xi et al., 2017). Recent protocols to create Hu-MuSCs-like cells from hiPSCs (Takahashi et al., 2007) could generate skeletal myogenic lineage cells, but with limited muscle regenerative capacity (Borchin et al., 2013; Magli et al., 2017; Shelton et al., 2014; van der Wal et al., 2018; Xi et al., 2017). PAX7 cell engraftment capability was often not reported or showed limited success. Transgene-induced PAX7 or PAX3 expression in hiPSCs can increase engraftment of PAX7⁺ cells and contribution to the SC pool (Al Tanoury et al., 2020; Wu et al., 2018); however, these engineered cells may not reflect physiology due to the genetic manipulation of these master transcription factors.

In this study, we used primary human muscle stem cells isolated from deceased FOP subjects and hiPSC-derived muscle progenitor cells from FOP subjects to elucidate how ACVR1^R206H affects muscle stem cell function and regenerative ability. Primary FOP Hu-MuSCs can engraft and regenerate injured muscle of immunocompromised (NSG) mice, but at lower efficiency than Hu-MuSCs from non-FOP donors, and the source of FOP Hu-MuSCs (biceps vs. diaphragm) appears to impact their engraftment efficiency. A new non-transgenic strategy for creating human iMPCs was developed and applied to multiple hiPSC lines from control subjects and subjects with FOP. These iMPCs shared transcriptional similarities with primary Hu-MuSCs, and abnormal activation of the BMP pathway by ACVR1^R206H changed the transcription profiles of the FOP iMPCs compared to controls.

Results

Lower efficiency engraftment of primary human FOP Hu-MuSCs

Muscle tissue samples were obtained from FOP cadavers which allowed the study of markers in situ and the isolation of primary human cells with endogenous activated ACVR1 signaling. Hematoxylin and eosin and alcian blue staining of FOP primary muscle samples from two deceased FOP subjects in a region without HO showed no gross defects. FOP muscle tissue near a HO lesion showed increased ECM proteoglycan components (alcian blue staining) in the interstitial space of the muscle fibers near the HO lesion (Figure 1A). Primary Hu-MuSCs carrying the ACVR1^R206H activating mutation were isolated from unfixed muscle tissue from two FOP autopsies. Biceps brachii muscle, commonly affected by HO, and diaphragm muscle, one of the rare skeletal muscle sites spared from HO in patients with FOP, were analyzed. PAX7 staining confirmed the presence of SCs in human FOP biceps (Figure 1B, middle and right). Interestingly, a small subregion showed Collagen Type 1 expression with embedded satellite cells (Figure 1B, right) suggesting possible early HO formation despite absence of gross HO.

Figure 1 with 1 supplement see all

Download asset Open asset

Regenerative properties of primary Hu-MuSCs isolated from FOP patients.

(A) Hematoxylin and eosin staining (top) and alcian blue staining (bottom) of muscle cross sections from FOP subjects, with no heterotopic ossification (left) and with heterotopic ossification, depicted by the white dashed line (right, 100 µm scale bar). (B) Immunofluorescence staining for PAX7 and COL1 of muscle cross sections from control and FOP subjects. White arrowheads mark Hu-MuSCs, 50 µm scale bar. (C) Heat map of normalized gene expression of sorted human satellite cells from two control subjects and two different muscles of 1 FOP subject. (D) Schematic and experimental time course. Hu-MuSCs from biceps and diaphragm muscles of a deceased FOP patient were sorted and transplanted into NSG mice. (E) Immunofluorescence staining 5 weeks after xenotransplantation. White arrowheads mark Hu-MuSCs, 50 µm scale bar. (F) Quantification of human DYSTROPHIN fibers and human PAX7⁺ cells 5 weeks after transplant. (G) Mice were re-injured 5 weeks after transplantation with bupivacaine. Immunofluorescence staining was performed at week 10. White arrows mark Hu-MuSCs, 50 µm scale bar. (H) Quantification of human DYSTROPHIN fibers and human PAX7⁺ cells after re-injury at week 10. n = 1, biological replicates, n ≥ 3 technical replicates. Error bars represent mean and SD. *, p < 0.05, **, p < 0.01. Muscle specimen and transplantation details are listed in Figure 1—source data 1 and Figure 1—source data 2 .

Figure 1—source data 1 Muscle specimen information.: https://cdn.elifesciences.org/articles/66107/elife-66107-fig1-data1-v2.docx
Download elife-66107-fig1-data1-v2.docx
Figure 1—source data 2 Hu-MuSCs transplantation details.: https://cdn.elifesciences.org/articles/66107/elife-66107-fig1-data2-v2.docx
Download elife-66107-fig1-data2-v2.docx

Sufficient Hu-MuSCs (Figure 1—figure supplement 1A) were sorted for transplant and gene expression analysis. FOP Hu-MuSCs showed lower PAX7 expression compared to control Hu-MuSCs. COL1A1 and ID3 were increased in affected FOP muscle (biceps) (Figure 1C). However, ID1 was increased in both the FOP diaphragm and biceps compared to control muscles (Figure 1C). Hu-MuSCs were transplanted into the irradiated tibialis anterior (TA) muscle of immunocompromised (NSG) mice (to hinder endogenous satellite cells) previously injured with bupivacaine (Garcia et al., 2017). The bupivicaine step induces myofiber injury and promotes engraftment of donor Hu-MuSCs. Five weeks after transplantation, biceps and diaphragm FOP Hu-MuSCs had engrafted and formed human fibers (Figure 1D and E). However, the number of human DYSTROPHIN-positive fibers was significantly lower with FOP Hu-MuSCs vs. control Hu-MuSCs (Figure 1F). The number of engrafted PAX7⁺ cells was qualitatively lower with FOP vs. control Hu-MuSCs (Figure 1F and Figure 1—figure supplement 1B). Ten weeks after re-injury with bupivacaine at week 5 (Figure 1D), the number of human DYSTROPHIN fibers was significantly decreased when FOP biceps Hu-MuSCs were transplanted compared to control Hu-MuSCs, but not when FOP diaphragm Hu-MuSCs were transplanted (Figure 1G and H). No differences in the number of human PAX7 cells were identified (Figure 1G and H). No radiologic evidence of HO was found in any mice (Figure 1—figure supplement 1C).

Thus, primary Hu-MuSCs isolated from two FOP subjects can engraft and regenerate injured muscle of NSG mice. Our limited sample size suggests this may occur at lower efficiency than control MuSCs. In addition, the results suggest that the source of FOP Hu-MuSCs (biceps or diaphragm) may impact MuSC engraftment efficiency.

Human FOP iPSCs can differentiate into skeletal muscle cells

The rarity of FOP disease and difficulty obtaining human tissue samples from patients with FOP makes it difficult to obtain a reliable source of muscle stem cells. Therefore, we used established and fully characterized control hiPSCs (Wtc11, 1323–2, and BJ2) and ACVR1^R206H hiPSCs (F1-1, F2-3, F3-2) lines previously derived from patients with FOP (Matsumoto et al., 2013). We first asked whether our lines were able to differentiate into myogenic cells and specifically into muscle stem cell-like cells using an adapted version of recently published protocols (Hicks et al., 2018; Shelton et al., 2014; Figure 2A). The control hiPSCs yielded PAX7, MYOGENIN-expressing cells (Figure 2B and C and Figure 2—figure supplement 1A), expressed DYSTROPHIN (Figure 2B), and formed contractile myotubes (Figure 2C; Video 1 and Video 2).

Figure 2 with 1 supplement see all

Download asset Open asset

Skeletal muscle differentiation of human control and FOP hiPSC lines.

(A) Differentiation schematic of hiPSCs into skeletal muscle stem cell-like cells (iMPCs). (B) Immunofluorescence staining for PAX7, DYSTROPHIN and nuclear stain DAPI of iMPCs at day 50, 50 µm scale bar. (C) Immunofluorescence staining of PAX7 and MYOGENIN expressing cells. Control and FOP hiPSCs can form contractile myotubes, 100 µm scale bar. (D) *PAX7*, *MYF5*, *MYOD1,* and *MYOGENIN* gene expression at day 50 of differentiation (n = 3 biological replicates and n ≥ 3 technical replicates). Error bars represent mean ± SD. No statistically significant differences were detected. (E) Schematic describing the addition of the BMP pathway inhibitor LDN193189 (LDN). (F) Representative immunofluorescence staining of 2 control and 2 FOP hiPSC lines differentiated and stained for MHC (Myosin Heavy Chain) and nuclear stain DAPI at day 50±LDN addition, 200 µm scale bar.

Video 1

Download asset

posterframe for video — Myotube contraction of differentiated Control (1323-2) hiPSCs.

Video 2

Download asset

Since BMPs control skeletal muscle differentiation from hiPSCs (Chal et al., 2016; Xi et al., 2017), we investigated if genetic activation of the BMP pathway via ACVR1^R206H could alter the myogenic differentiation of hiPSCs (F1-1, F2-3, F3-2) derived from patients with FOP (Matsumoto et al., 2013). All three FOP iPSC lines formed contractile myotubes with cells expressing PAX7, MYOGENIN, and DYSTROPHIN (Figure 2B and C, Figure 2—figure supplement 1A). PAX7, MYF5, MYOD1, and MYOGENIN (Figure 2D) gene expression showed some heterogeneity among the different hiPSC lines but these differences were not statistically significant. The heterogeneity seen is likely due to the presence of other cell types in our differentiation such as neuronal cells (Figure 2—figure supplement 1C). Adding a SMAD inhibitor of the BMP pathway (LDN193189) into the differentiation protocol (Figure 2E) improved the differentiation of the control hiPSC lines. LDN also improved differentiation of the FOP F3-2 (Figure 2F) line, which had shown lower PAX7 and MYOGENIN expression (Figure 2D). These findings show that individual hiPSC lines are heterogeneous in differentiation to skeletal muscle lineages, similar to other protocols (Volpato and Webber, 2020), that FOP hiPSCs can form skeletal muscle cells and Hu-MuSC-like cells expressing PAX7 despite upregulation of the BMP pathway by ACVR1^R206H, and that chemical blockade of the BMP pathway can improve the formation of Hu-MuSC-like cells from the FOP iPSC line that showed the lowest efficiency.

PAX7 expressing cells can be isolated from myogenic differentiation

To test the regenerative properties of PAX7-expressing MuSCs, we used FACS to purify HNK1^-CD45^-CD31^- cells co-expressing CD29, CXCR4, and CD56 markers present on human PAX7⁺ cells (Garcia Steven et al., 2018; Figure 3A and Figure 3—figure supplement 1A–C). FACS analysis identified intermediate CD56 cells expressing high PAX7 and low MYOGENIN (Figure 3B and Figure 3—figure supplement 1C,D), consistent with Hu-MuSC expression profiles. All HNK1^-CD45^-CD31^- CXCR4⁺CD29⁺CD56^dim cells formed myotubes expressing MHC (Figure 3—figure supplement 1E) when cultured in terminal differentiation media demonstrating isolation of functional iMPCs with satellite cell-like characteristics from the cultures.

Figure 3 with 1 supplement see all

Download asset Open asset

Isolation and transplantation of PAX7-expressing cells from hiPSC muscle differentiation culture.

hiPSCs were differentiated into iMPCs until day 50 and sorted via flow cytometry. (A) Gating strategy. (B) Myogenic gene expression of CD56^dim and CD56^hi cells. CD56^dim vs CD56^hi, n = 6 (3 Ctrl and 3 FOP lines) biological replicates, ** p < 0.01, *** p < 0.001, **** p < 0.0001. No significant differences were found between control and FOP lines. (C) Representative human DYSTROPHIN (top, 200 µm scale bar), and human LAMIN A/C, human SPECTRIN, and PAX7 (bottom, 100 µm scale bar) immunohistochemistry of NSG mice anterior tibialis, where sorted iMPCs were transplanted. White arrows show engrafted hiPSC-derived muscle stem cells. (D) Quantification of human DYSTROPHIN fibers and human PAX7 cells at week 5 after transplant (n = 3 biological replicates). Transplantation details are in Figure 3—source data 1.

Figure 3—source data 1 hiPSC-derived HNK1^-CD45^-CD31^- CXCR4⁺CD29⁺CD56^dim cell transplants.: https://cdn.elifesciences.org/articles/66107/elife-66107-fig3-data1-v2.docx
Download elife-66107-fig3-data1-v2.docx

Isolated iMPCs can regenerate injured mouse muscle and form human fibers

To assess iMPC regenerative capacity in vivo, we injected 1000–10,000 iMPCs derived from control or FOP hiPSCs (Figure 3—source data 1) into the irradiated TA muscle of NSG mice. New fibers expressing human DYSTROPHIN and PAX7 cells were found after 5 weeks, showing that iMPCs could engraft and promote muscle regeneration (Figure 3C). However, the number of human fibers and human PAX7⁺ cells remained low (Figure 3D) compared to primary Hu-MuSCs. While some iMPC transplants yielded up to 60 new human fibers, some did not yield any human fibers. By comparison, 2000 primary non-FOP Hu-MuSCs resulted in an average of 155 human fibers based on prior assessments using the same assay (Garcia Steven et al., 2018). No significant differences between control and FOP iMPCs were identified, though some individual FOP samples showed higher engraftment (Figure 3D).

These results showed that iMPCs have engraftment potential into a muscle injury site in mice, but engraftment efficiency may be lower than primary Hu-MuSCs or be the result of differences in experimental conditions. Also, ACVR1^R206H did not significantly impact muscle fiber regeneration in this assay.

Transcriptional profiling of iMPCs

The lower engraftment of iMPCs compared to primary Hu-MuSCs suggested that the FACS-purified population was still heterogeneous or that iMPCs do not fully recapitulate adult primary Hu-MuSCs. Single cell RNA sequencing (scRNAseq) from control (1323–2) and FOP (F3-2) iMPCs (these lines were selected based on their lower intra-line variability, Figure 2) were analyzed (Figure 4A and Figure 4—source data 1). Cell populations for both samples were defined by the dimension reduction technique of uniform manifold approximation and projection (UMAP) (Becht et al., 2018) and unsupervised clustering with Seurat v3 package (Stuart et al., 2019; Figure 4—figure supplement 1A–B). Both control (Figure 4—figure supplement 1A) and FOP (Figure 4—figure supplement 1B) samples had clusters expressing myogenic genes (PAX7 and MYOD); mesenchymal genes (PDGFRA); and neuronal genes (SOX2) (Figure 4—figure supplement 1C). Merged analysis to allow direct comparison identified 13 distinct clusters (Figure 4B and Figure 4—figure supplement 1D). Cells expressing muscle markers (PAX7, MYF5, MYOD, and MYOG [MYOGENIN]) were found in clusters 0–2 (Figure 4B–C and Figure 4—figure supplement 1E). Mesenchymal genes (PDGFRA, ASPN, and COL1A1) were expressed in cluster 3. Cluster four consisted cells expressing retinal (OTX2) and neuroprogenitor cell (NPC) markers (SOX2, MAP2). Neuroprogenitor markers (SOX2, MAP2) and HES6 were upregulated in cluster 5. Clusters 6 and 7 were made of cycling neuroprogenitor cells (MKI67, TOP2A). We also found three clusters (8-11) to be neuroepithelial cells (NECs), with cluster 10 having and increased in expression of mTor gene downstream target. Glial cells (TNC, SLC1A2) were comprised in cluster 12 (Figure 4C and D and Figure 4—figure supplement 1E). The frequency of muscle cells (clusters 0–2) was higher in FOP (14%, 2.5%, and 1.7%) compared to controls (4.2%, 0.5%, and 0.5%) (Figure 4F) suggesting that FOP hiPSCs may be more efficient at making muscle progenitor cells. Thus, FACS-purified iMPC cultures contain muscle stem/progenitor cells, but other cell types such as mesenchymal, neuronal progenitor cells, and myoblasts persist in the culture. Furthermore, the frequency of iMPCs appears to be higher in FOP vs. control cell cultures.

Figure 4 with 1 supplement see all

Download asset Open asset

Single-cell RNA sequencing of HNK1^-CD45^-CD31^-CXCR4⁺CD29⁺CD56^dim cells.

(A) Schematic of differentiating hiPSCs, sorting of iMPCs, and scRNAseq. (B) UMAP visualization plots of cells combined from both control and FOP samples. Cell types were assigned based on the expression of marker genes (Figure 4—figure supplement 1), MPC (muscle progenitor cells), NPCs (neuroprogenitor cells) and NECs (neuroepithelial cells). (C) Feature expression plots showing the localization of cells expressing myogenic markers (*PAX7*, *MYF5*), mesenchymal (*COL1A1*), and neural cell marker (*SOX2*). (D) Dot plot displaying expression genes associated with myogenesis, mesenchymal and neurogenesis markers. (E) Proportion of cells per cluster for each sample.

Figure 4—source data 1 Quality control information for each sample.: https://cdn.elifesciences.org/articles/66107/elife-66107-fig4-data1-v2.docx
Download elife-66107-fig4-data1-v2.docx

FACS-sorted iMPC transcriptome is heterogeneous

Typical muscle progenitor cultures are expected to contain cells undergoing expansion, differentiation, and maturation. Sub-clustering (Figure 5A and B) identified five new myogenic subpopulations (Figure 5A and B). Clusters 0–3 expressed higher levels of PAX7 and MYF5 (markers of more quiescent MuSCs) while cluster four expressed higher levels of MYOD (or MYOD1), MYOG, SOX8, and MEF2C (markers of differentiated MuSCs/myoblasts/myocytes). Muscle stem-like cells expressing APOE, KRT17, and CAV1 defined cluster 0. Quiescence markers (SPRY1, HEY1, HES1) were also highly expressed in this cluster. CDH15, a niche regulator of SCs quiescence (Goel et al., 2017) was enriched in cluster 1. Cluster 1 also had increased expression of DES, CHRNA1, and MYOD1 suggesting that these cells have a transcriptome resembling the profile of progenitor cells. Cluster 2 consisted of muscle stem-like cells expressing high levels of FOS (Figure 5B and Figure 5—figure supplement 1A). Cell cycle (TOP2A or KI67) markers were increased in cluster 3 (Figure 5B and Figure 5—figure supplement 1B (left)). Cluster 3 had a higher proportion of cells in G2M and S phase (Figure 5—figure supplement 1B (left)). The cell cycle distribution was similar in control and FOP (Figure 5—figure supplement 1B (right)). The proportion of cells in clusters 3 and 4 were similar in both control (10% and 12.7%) and FOP (14.3% and 9%) samples. The proportion of cells in clusters 0 and 3 was higher in FOP (31.4% and 23.4%) compared to control (20% and 5.5%), while the proportion of cells in cluster two was increased in control (51.8% vs 21.9%) (Figure 5—figure supplement 1D).

Figure 5 with 1 supplement see all

Download asset Open asset

Transcriptional profile of iMPCs.

(A) Identified myogenic clusters (in red) were sub-clustered from the other clusters (neurogenic and mesenchymal) and re-analyzed. UMAP of the new myogenic subclusters was generated. (B) Dot plot displaying expression of cluster defining genes. (C) Feature expression plots of myogenic markers. (D) Myogenic cells were ordered as a function of pseudotime using the Monocle package. (E) Pseudotime trajectory plot of myogenic clusters.Arrows represents the direction and major branches of the pseudotime. Branches A-C are marked with letters along the trajectories. (F) Gene plots displaying the expression of specific myogenic genes as a function of pseudotime. (**G–H**) Violin plots of myogenic (G) and chondro/osteogenic genes (H) that were significantly differentially expressed. Violin plot width depicts the larger probability density of cells expressing each particular gene at the indicated expression level. *, significantly different following differential expression testing using the Wilcoxon rank sum test per cluster. (I) SMAD and P38MAPK pathway markers significantly differentially expressed between control and FOP myogenic cells. p values for (**G–H**) are in Figure 5—source data 1.

Figure 5—source data 1 Differential expression analysis of control and FOP myogenic cells.: https://cdn.elifesciences.org/articles/66107/elife-66107-fig5-data1-v2.xlsx
Download elife-66107-fig5-data1-v2.xlsx

We ordered the myogenic cells into three major branches using the Monocle analysis package (Trapnell et al., 2014) based on genes that differ between clusters and constructed pseudotime differentiation trajectories (Figure 5D). The pseudotime ordering of the cells (Figure 5D) showed that cycling MuSCs (cluster 3) were located at the tip of the tree branch A while cluster 0 comprised of the most quiescent MuSCs were mainly distributed along branch A. Cells from the FOS upregulated MuSCs (cluster 2) were distributed along branch B. Cluster 1 (muscle progenitor cells) was distributed in the more proximal parts of branches B and C. Distal parts of branches B and C, were notably comprised of more mature cells from cluster 4 (myoblasts) (Figure 5E). PAX7 and MYF5 were upregulated in branch A but downregulated in the tip of branches B and C (Figure 5C and E) confirming the correct cell ordering. MYOD and MYOG showed higher expression patterns early in branch B and late in branch C (Figure 5E and F).

Thus, hiPSC differentiation cultures contain subpopulations of iMPCs showing the expected spectrum of quiescence, activation, and differentiation with FOP cultures having a higher proportion of cells in the stem cell/progenitor and proliferating phases and fewer mature myoblasts.

FOP iMPCs cells express increased chondro/osteogenic and ECM markers

Differential expression analysis on the transcriptional profiles of the sub-clustered myogenic cells (Figure 5A and B and Figure 5—figure supplement 1A-E) was used to see if ACVR1^R206H altered transcriptional signatures. While PAX7 was significantly increased in control vs. FOP cells in cluster three only, MYF5 was significantly increased in FOP cells from cluster 1 compared to control cells (Figure 5G, Figure 5—source data 1). Since Hu-MuSCs show heterogeneous levels of PAX7 and MYF5 expression (Kuang et al., 2007), this suggests the ACVR1^R206H mutation may favor one sub-population over another. Interestingly, SPRY1, a known regulator of quiescence (Shea et al., 2010) which decreases with age (Bigot et al., 2015), was significantly downregulated in FOP cells in all the clusters, while DLK1, which act as a muscle regeneration inhibitor (Andersen et al., 2013) was significantly increased in FOP cells (Figure 5G and Figure 5—source data 1).

Since the ACVR1^R206H mutation increases BMP pathway activity and expression of chondrogenic and osteogenic markers in multiple lineages (Barruet et al., 2016; Culbert et al., 2014; Matsumoto et al., 2013), we examined the expression levels of extracellular matrix, fibrogenic, chondrogenic, and osteogenic genes. ITM2A, COL1A1, DCN, CD151, SPARC, and TIMP1 were significantly increased in FOP cells (Figure 5H, Figure 5—source data 1). ECM proteoglycans known to be involved in inflammation, including BGN (Nastase et al., 2012) and LUM (Nikitovic et al., 2014), TAGLN which regulates osteogenic differentiation (Elsafadi et al., 2016), and IGBP5 which is increased in aged satellite cells (Soriano-Arroquia et al., 2016), were also increased in FOP cells (Figure 5—figure supplement 1E). We also investigated chondrogenic and osteogenic genes that have been described to be increased in other FOP cell types (Culbert et al., 2014; Hino et al., 2015; Matsumoto et al., 2013) such as SOX9, SP7, RUNX2, ACAN, and BGLAP. Overall, those markers had very low level of expression in addition to have few cells expressing them. However, a higher number of cells expressed SOX9 in FOP than control cells (Figure 5—figure supplement 1F-G). SP7 was not found to be expressed in any of the samples or clusters.

Finally, expression of target genes of the BMP pathway (ID1, ID3, BMPs, and SMADs) and the p38MAPK pathway (Figure 5I and Figure 5—figure supplement 1H) was assessed to see if activated ACVR1 altered these pathways. ID1 was significantly higher in clusters 0 (MuSCs 1), 1 (MPCs), 3 (Cycling MuSCs), and 4 (Myoblasts) in control cells, while ID3 was significantly higher in FOP cells clusters 1 (MPCs) and 2 (MuSC 2). In addition, the BMP/TGFβ pathway downstream target gene SMAD5 was significantly higher in cluster 2 (MuSCs 2) of control cells. The p38 pathway components JUN (clusters 1, 3), JUNB (clusters 0–2), and FOS (clusters 1, 3) were significantly increased in FOP cells. Within the known ACVR1 co-receptors, BMPR2 (clusters 3, 4) and ACVR2B (clusters 0, 2, 3) expression were significantly higher in control cells while ACVR2A expression was significantly higher in FOP cells (clusters 1, 4). ACVR2B was significantly higher in control cells (clusters 0, 1, 4) (Figure 5I). Similar to the primary FOP Hu-MuSCs (Figure 1), these results suggest that FOP iMPCs may have a chondrogenic/osteogenic signature, increased ID3 expression, and also showed higher p38 pathway activity and higher levels of the ACVR2A co-receptor at different stages of myogenic differentiation.

iMPC transcriptome shows similarities to primary Hu-MuSCs

Comparing the iMPC scRNAseq to primary Hu-MuSCs data of sorted satellite cells from a human vastus lateralis muscle (Barruet et al., 2020) was used to identify if their lower engraftment efficiency was due to transcriptional differences. The merged data UMAP (Figure 6A) showed that clusters 0–5 contained myogenic cells (Figure 6B and C). PAX7⁺ + were identified in clusters 0–4 (Figure 6C). Myocyte contaminants present in the primary sorted cells constituted cluster 5. Although iMPCs expressed PAX7 and MYF5, the gene expression levels were higher in primary Hu-MuSCs. In contrast, MYOD was higher in iMPCs (Figure 6C and Figure 6—figure supplement 1). Mesenchymal cells were identified in cluster six while NPCs, NECs and Glial cells were identified in clusters 7–9, 10–12 and 13, respectively (Figure 6B and C).

Figure 6 with 1 supplement see all

Download asset Open asset

iMPC transcriptional signature compared to human primary muscle stem cells.

(**A–B**) UMAP of cells combined from human primary muscle stem cells (vastus muscle) and the control and FOP samples. (A) UMAP showing the distribution of cells per sample, and (B) with clusters labeled. (C) Dot plot displaying expression of cluster defining genes. (D) Myogenic cells (in red) are comprised in cluster 0–5. UMAP of the myogenic clusters. (**E–G**) Pseudotime trajectory plot generated via Monocle analysis depicting all myogenic clusters (E) and samples (F). (G) Level of expression of ECM/osteogenic genes along the cell trajectories. (H) Heatmap representing genes that are significantly branch dependent using the BEAM analysis (Figure 6—source data 1) and also genes that have similar lineage-dependent expression patterns. Branches are shown in the upper right panel.

Figure 6—source data 1 BEAM analysis.: https://cdn.elifesciences.org/articles/66107/elife-66107-fig6-data1-v2.xlsx
Download elife-66107-fig6-data1-v2.xlsx

Detailed analysis of the myogenic cell subset (cluster 0–5, Figure 6D) was performed using pseudotime trajectory analysis to elucidate the states of the iMPCs with respect to primary Hu-MuSCs. Cells from clusters 0, 1, and 4 were distributed along branches B and C. Myocytes ordered at the distal end of branch A. iMPCs ordered away from primary Hu-MuSCs (PAX7+) and primary myocytes (MYL1+) (Figure 6E–G). Branch expression analysis modeling (BEAM) allowed us to investigate significant gene that are branch-dependent in their expression (Figure 6—source data 1). Branch A was mainly constituted of primary Hu-MuSCs while branch B consisted of iMPCs and subset of primary Hu-MuSCs, expressed significantly higher levels of genes associated with mesenchymal, fibrogenic, chondrogenic, osteogenic lineages and extracellular matrix (Figure 6G and H). Thus, iMPCs retained strong bi-potency compared to primary Hu-MuSCs, suggesting that iMPCs may not be as committed to the muscle lineage as primary adult muscle stem cells (Xi et al., 2020).

Discussion

Developing optimal strategies for skeletal muscle regeneration and repair requires a detailed understanding of how these processes are regulated. Our study uses the abnormal ACVR1 activity of FOP to dissect how the BMP/ACVR1 pathways impact skeletal muscle repair (Figure 7). In the case of FOP, muscle injury can trigger disease progression or severe complications. Thus, generating large numbers of human iMPCs using transgene-free protocols holds promise for understanding pathological mechanisms and finding new therapeutic targets. The findings in FOP are also likely to help inform why HO sometimes forms in other conditions of muscle injury or inflammation, including burns, CNS injury, rheumatologic disorders, trauma, and surgery (Matsuo et al., 2019).

Figure 7

Download asset Open asset

Testing the role of activated ACVR1 signaling in muscle repair.

FOP iMPCs and FOP primary human muscle stem cells are used for modeling muscle stem cell engraftment and regeneration properties.

This study revealed several novel findings, including deficiencies in the muscle repair capacity of FOP Hu-MuSCs obtained from two different subjects, abnormal production of extracellular matrix components by iMSCs and FOP Hu-MuSCs that may contribute to FOP disease progression, and important differences between primary Hu-MuSCs and iMPCs that may impact other studies. In addition, our results suggest that there are muscle specific differences in FOP Hu-MuSC repair capacity, potentially explaining why some skeletal muscles appear to be protected from developing HO.

The iMPC and primary Hu-MuSC single cell transcriptomes revealed clear subsets of cells at all stages of muscle differentiation (quiescent, activated, and differentiated). Satellite cell subtype markers e.g. COLs, DLK1, ID3, and HES1, (Barruet et al., 2020) were also highly expressed in iMPCs suggesting that in vitro differentiation may favor specific subtypes of satellite cells. iMPCs also expressed high levels of mesenchymal/ECM markers, similar to that in human fetal muscle progenitor cells (Xi et al., 2020). These results suggest that iMPCs may retain a more progenitor-like phenotype as compared to Hu-MuSCs. In addition, we showed that PAX7⁺ iMPCs can engraft in muscle injury models in mice but at lower efficiency.

This iMPC system revealed that activation of the BMP pathway by the FOP ACVR1^R206H mutation induced changes that could contribute to abnormal muscle healing. As expected, the classical SMAD pathways were active in FOP cells; however, it was unexpected to find that the p38 pathway via JUN/FOS was more active in the FOP iMPCs, considering that FOP cultures had a higher proportion of cells in the stem cell/progenitor and proliferating phases. While increased SMAD activity (Billings et al., 2008) or misinterpretation of the Activin A ligand (Hatsell et al., 2015) occur with the ACVR1^R206H mutation, abnormal p38 signaling may be critical in some cell types like macrophages (Barruet et al., 2018). Since p38 is a major regulator of Hu-MuSC function (Segalés et al., 2016) and has been associated with inflammation and ECM accumulation in aged muscle (Cosgrove et al., 2014), further studies are needed to determine how this p38 signaling contributes to satellite cell subtype specification (Barruet et al., 2020) and affect healing in FOP.

Importantly, our iMPC and Hu-MuSC transplant studies showed no HO in the recipient mice, consistent with prior FOP genetic studies (Dey et al., 2016; Lees-Shepard et al., 2018) describing that other cell types may be directly responsible for bone formation. Indeed, several reports have suggested fibro/adipogenic progenitors (FAPs) as a candidate cell population to FOP pathogenesis. Interestingly, there are functional cross-talk between satellite cells and FAPs which is critical for muscle repair in homeostasis, disease, and aging contexts (Biferali et al., 2019; Farup et al., 2015; Lukjanenko et al., 2019; Wosczyna and Rando, 2018). However, FOP iMPCs showed increased chondrogenic/osteogenic and ECM gene expression, a feature seen in other bone-related cell types in FOP (Barruet et al., 2016; Culbert et al., 2014; Lees-Shepard et al., 2018). Thus, one possibility would be for FOP Hu-MuSCs to further contribute to HO formation indirectly by modulation of the osteogenic environment such as contributing to changes in muscle stiffness, as previously reported in FOP mice (Haupt et al., 2019; Stanley et al., 2019). This has also been reported in aged muscle, which show increased myofiber stiffness due to increased ECM deposition (Hwang and Brack, 2018). In addition, while our results suggest no major differences or possibly a slight increase in the ability to form skeletal muscle progenitors, activated ACVR1 by the R206H mutation decreased (but did not abrogate) in vitro formation of mature myoblasts of FOP hiPSCs and in vivo muscle repair after transplant of primary FOP Hu-MuSCs.

Comparing primary FOP Hu-MuSCs from muscles that develop HO (biceps) and non-affected muscle (diaphragm) suggests that the source of the Hu-MuSCs may impact engraftment efficiency. The re-injury model showed that engraftment of Hu-MuSCs from FOP biceps, but not from diaphragm, remained significantly decreased. Thus, it is intriguing to consider that the clinical sparing of the diaphragm from HO in patients with FOP may result from a less impaired or unimpaired muscle repair process in diaphragm satellite cells. Further delineation of muscle-specific Hu-MuSC properties will be revealing to understand this observation. Our finding that primary FOP Hu-MuSCs have lower engraftment ability provides a potential explanation for the poor skeletal muscle repair observed in patients with FOP (Shore, 2012). However, the FOP and control iMPCs showed no major differences in engraftment, possibly due to decreased assay sensitivity from the lower engraftment efficiency of iMPCs in general, which could be explained by iMPCs being more immature than primary Hu-MuSCs, the presence of cell contaminants after isolation, or that iMPCs represent a subtype of cells that may be more reflective of non-ossifying skeletal muscle like diaphragm.

This study has several limitations. Our transplant experiments used immunocompromised mice where mature B and T cells are absent and macrophages are defective (Shultz et al., 2005). This may dampen the engraftment/regenerative phenotypes, particularly as immune cells are important in muscle regeneration (Furrer and Handschin, 2017) and there is growing awareness of the role of the immune system in the pathogenesis of FOP (Barruet et al., 2018; Convente et al., 2018). Also, while our FACS strategy can isolate Hu-MuSCs and iMPCs without the need of transcription factors or markers, we noted that the engraftment efficiency still varied among the iMPC lines despite FACS purification. This heterogeneity is an ongoing problem among all published iMPC protocols to date. In our case, this was counter-balanced by our consistent findings across multiple lines and between iMPCs and the availability of primary Hu-MuSCs despite the absence of isogenic control and FOP lines. Although our primary cell studies were limited by the rarity of FOP (estimated at 1 in 1.4 million people) and even rarer suitable cadaveric samples where we were only able to obtain Hu-MSCs from two FOP subjects, the combination of iPSC-derived lineages with the rare primary samples provided multiple avenues for supporting our conclusions. Finally, the ACVR1^R206H mutation increases BMP signaling but also introduces neofunction to Activin A. Our studies are not able to distinguish between these two contributing pathways. Our finding of increased p38 activity in the FOP iMPCs was also seen in FOP subject monocyte-derived macrophages (Barruet et al., 2018), suggesting that this alternate signaling pathway by ACVR1 should also be investigated further. Future studies, including the use of isogenic control lines, testing of neutralizing antibodies to directly target the ACVR1^WT or ACVR1^R206H receptors, and the use of other small molecule inhibitors of pathways downstream of ACVR1 (e.g. p38 or SMAD), will help elucidate the different factors that contribute to the hiPSC-line specific effects, including potential roles for disease modifier genes or BMP pathway modulators.

This study shows that human iPSC-derived muscle progenitor cells can be a valuable tool to model musculoskeletal diseases of skeletal muscle injury and repair. Correlating the findings in iMPCs with primary Hu-MuSCs revealed an indirect role for skeletal muscle progenitors in HO formation, as well as potential subtypes of Hu-MuSCs that could contribute to skeletal muscle specific regenerative capacity (Figure 7). These studies highlight the importance of skeletal muscle regeneration in disease pathogenesis and establish a foundation for understanding how skeletal muscle repair and osteogenesis are linked.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Strain, strain background (Mus musculus)	NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ	https://www.jax.org/strain/005557	MGI:3577020	8–12 weeks old
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	ACTBHs01060665_g1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	DYSTROPHINHs007758098_m1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	PAX7Hs00242962_m1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	CD56Hs00941830_m1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	MYOD1Hs00159528_m1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	MYF5Hs00929416_g1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	MYOGENIN Hs01072232_m1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	COL1A1Hs01076780_g1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	ID1Hs03676575_s1
Sequence-based reagent	Human RT-PCR Primers	Applied Biosystems Taqman Assays	ID3Hs00954037_g1
Antibody	Anti-Human DYSTROPHIN (mouse monoclonal)	DSHB	AB_2618157	IF(1:10)
Antibody	Anti-Human/Mouse DYSTROPHIN (mouse monoclonal)	Thermofisher	AB_10978300	IF(1:500)
Antibody	Anti-Human/Mouse PAX7 (mouse monoclonal)	DSHB	AB_528428	IF(1:10)
Antibody	Anti-Human/Mouse MHC (mouse monoclonal)	DSHB	AB_2147781	IF(1:100)
Antibody	Anti-LAMININ (rabbit polyclonal)	Sigma-Aldrich	AB_477163	IF(1:250)
Antibody	Anti-Human SPECTRIN (mouse monoclonal)	Leica Microsystems	AB_442135	IF(1:100)
Antibody	Anti-Human LAMIN A/C (mouse monoclonal)	Vector Laboratories	AB_2336546	IF(1:100)
Antibody	Anti-Human CD31 (Beads) (Mouse monoclonal)	Miltenyi Biotec	130-091-935	2 μl
Antibody	Anti-Human CD45 (Beads, mouse monoclonal)	Miltenyi Biotec	AB_2783001	5 μl
Antibody	Anti-Human CD31 AF450(WM-59, mouse monoclonal)	Ebioscience	AB_10854276	5 μl
Antibody	Anti-Human CD34 eFluor450(4H11, mouse monoclonal)	Ebioscience	AB_10733282	5 μl
Antibody	Anti-Human CD45 AF450(30-F11, mouse monoclonal)	Ebioscience	AB_1518806	5 μl
Antibody	Anti-Human CD29 FITC(TS2/16, mouse monoclonal)	Ebioscience	AB_2043830	8 μl
Antibody	Recombinant human anti-CD56 APC-vio-770 (REA196)	Miltenyi Biotec	AB_2733136	8 μl
Antibody	Anti-Human CXCR4 PE(12G5, mouse monoclonal)	Ebioscience	AB_10669164	8 μl
Antibody	Anti-Human HNK1 PE (TB01, mouse monoclonal)	Ebioscience	AB_10804531	5 μl
Antibody	Anti-Human CD45 PE(30-F11, mouse monoclonal)	Ebioscience	AB_465668	5 μl
Antibody	Anti-Human CD31 PE(390, mouse monoclonal)	Ebioscience	AB_465632	5 μl
Antibody	Anti-Human CD56 APC (CMSSB, mouse monoclonal)	Ebioscience	AB_10854573	5 μl
Antibody	Anti-Human CXCR4 PE-Cy7 (12G5, mouse monoclonal)	Ebioscience	AB_1659706	5 μl
Antibody	FcR block	Miltenyi Biotec	AB_2892112	5 μl
Software, algorithm	GraphPad Prism	GraphPad Prism (https://graphpad.com)	SCr_002798
Software, algorithm	Seurat (3.1.5)	https://satijalab.org/seurat/	SCR_007322	SCR_021002
Software, algorithm	Monocle (2.12.0)	http://cole-trapnell-lab.github.io/monocle-release/	SCR_016339
Software, algorithm	cellranger	https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/feature-bc	SCR_021002
Software, algorithm	FlowJo	https://www.flowjo.com	SCR_008520

Cell culture and differentiation

Request a detailed protocol

Pluripotent hiPSC lines derived from control and FOP fibroblasts (Matsumoto et al., 2013; Spencer et al., 2014) were cultured in mTeSR1 medium (StemCell Technologies) on irradiated SNL feeder cells (McMahon and Bradley, 1990) as described previously. hiPSCs were passaged at least once on Matrigel (Corning)-coated plates (150–300 µg/ml) to remove the SNLs before use in differentiation assays. ROCK inhibitor Y-27632 (10 µM, StemCell Technologies) was added to mTeSR1 when cells were split and removed the following day.

hiPSC lines were differentiated into skeletal muscle cells using modifications based on prior protocols (Chal et al., 2016; Shelton et al., 2014). Our hiPSC lines differentiated better with a lower cell number seeding and a longer time of recovery between the seeding and the start of the differentiation (2 days, data not shown). Cells were seeded at 7.5 × 10⁵ cells per well of a 12-well plate on Matrigel two days before the differentiation medium (E6 medium, supplemented with either 10 µM CHIR99021 (Tocris) for 2 days or with 3 µM CHIR99021 and 0.5 µM LDN193189 for six days). Cells were then grown in un-supplemented E6 media until day 12, then changed to StemPro-34 media supplemented with 10 ng/ml bFGF until day 20. The medium was then replaced by E6 medium until day 35, when DMEM/F12 supplemented with N2 (Gibco) and Insulin-Transferrin-Selenium (ITS-A, 100 X Gibco) was added. Media was changed daily until harvest at day 50. After sorting, cells were plated in satellite cell media [DMEM/F12 (Gibco), 20 % FBS (Hyclone), 1 X ITS (Gibco), 1 X Penicillin/ Streptomycin (Gibco)] for further functional assays (Figure 2A and F). Once cells reached confluence, cells were cultured in differentiation media (DMEM, 2 % horse serum, 1 X Penicillin/ Streptomycin; Gibco).

Cell lines

Request a detailed protocol

Pluripotent hiPSC lines were previously derived from the following control and FOP fibroblasts (Matsumoto et al., 2013; Spencer et al., 2014): FF-WT-BJ (Foreskin, Stemgent 08–0027), HDF-WT-1323 (Fibroblast, Cell Applications 1323), HDF-WTc (Skin biopsies obtained from donors), HDF-FOP1 (Skin Fibroblast, Corrielle GM00513), HDF-FOP2 (Skin Fibroblast, Corrielle GM00783), and HDF-FOP3 (Skin biopsies obtained from donors). All hiPSCs are tested periodically for mycoplasma and have been negative to date. Cell lines are yearly authenticated using DNA fingerprinting.

Flow cytometry

Request a detailed protocol

HNK1^-CD45^-CD31^-CXCR4⁺CD29⁺CD56^dim cells were sorted from skeletal muscle differentiation of control and FOP hiPSCs. Human primary satellite cells were isolated and sorted as described (Garcia Steven et al., 2018; Garcia et al., 2017). Human muscle was freshly harvested and stored in DMEM with 30 % FBS at 4 °C overnight or for two extra days (delay due to shipping). Muscle samples were digested, erythrocytes were lysed, and hematopoietic and endothelial cells were depleted with magnetic column depletion using CD31, CD34, and CD45 (eBioscience). Cells were further gated as described in Figure 1—figure supplement 1A and sorted for CXCR4⁺/CD29⁺/CD56⁺ and collected for subsequent experimentation.

For the staining, cells were treated with Accutase for 20 min at 37 C, washed with FACS buffer, and stained for HNK1-PE, CD45-PE, CD31-PE, CD56-APC, CXCR4-PE-Cy7, and CD29-FITC (eBioscience). HNK1^-CD45^-CD31^- cells [to select against neuronal cells (HNK1, Human Natural Killer-1) (Choi et al., 2016; Hicks et al., 2018), hematopoietic cells (CD45), and endothelial cells (CD31)], co-expressing CD29, CXCR4, and intermediate CD56, markers present on human PAX7⁺ cells (Garcia Steven et al., 2018; Garcia et al., 2017; Xu et al., 2015) were sorted with a FACSAriaIII (BD Biosciences) and Sytox Blue (Life Technologies) was used as a viability marker. Alternatively, cells were permeabilized and fixed (Fix/Perm Buffer Set, BioLegend). Fixed cells were first incubated with primary antibodies PAX7 (DSHB), CD29 (BD Biosciences), and CD56 (BD Biosciences) following by secondary antibodies (Alexa350-conjugated goat anti-mouse IgG, Alexa488-conjugated goat anti-rat IgG and Alexa546-conjugated donkey anti-goat IgG, Life Technologies).

Animal care and transplantation studies

Request a detailed protocol

All mouse studies were performed using protocols approved by the UCSF Institutional Animal Care and Use Committee. Mice were either bred and housed in a pathogen-free facility at UCSF or purchased from The Jackson Laboratory. Eight to 12-week-old NSG (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ) mice were randomized to all experimental groups by sex and littermates. Each mouse was irradiated with 18 Gy before transplantation. Isolated primary human satellite cells (Hu-MuSCs) or iMPCs were injected with 50 µl 0.5 % bupivacaine directly into the tibialis anterior (TA) muscle of one leg as described (Garcia Steven et al., 2018; Garcia et al., 2017). The TA for each mouse was harvested at week 5 or week 10 after transplantation and frozen in O.C.T. compound in 2-methylbutane chilled in liquid nitrogen. Serial 6 µm transverse frozen sections were analyzed or stored at –80 °C.

Collected tibialis anterior frozen cross sections were fixed in 4 % PFA for 10 min at room temperature, washed with PBST PBS with 0.1 % Tween-20 (Sigma-Aldrich), and blocked in PBS with 10 % goat serum for 1 hr at room temperature. Slides were then incubated 4 hrs at room temperature with the following: mouse monoclonal anti-human DYSTROPHIN (DSHB), mouse monoclonal IgG1 anti-PAX7 (DSHB), rabbit polyclonal anti-Laminin (Sigma-Aldrich), mouse monoclonal IgG2b anti-human SPECTRIN (Leica Microsystems), and mouse monoclonal IgG2b anti-human LAMIN A/C (Vector Laboratories). After PBST wash, slides were incubated with the following secondary antibodies: Alexa Fluor 555 goat anti-mouse IgG, Alexa Fluor 594 goat anti-mouse IgG1, Alexa Fluor 488 goat anti-mouse IgG2b, and Alexa Fluor 647 goat anti-rabbit (Life Technologies). Finally, sections were mounted with VECTASHIELD mounting media with DAPI (Vector Laboratories). All samples were examined using a Leica upright or DMi8 Leica microscope. Sections with the most human fibers were used for human DYSTROPHIN and PAX7 quantification for each condition, as described below.

Cell Immunostaining and NSG Tibialis Anterior Analysis iMPCs were fixed with 4%PFA/PBS for 10 min at room temperature, permeabilized with 0.1%Triton-100X (Sigma-Aldrich), and blocked with 5 % BSA (Sigma-Aldrich). Cells were stained overnight with primary antibodies for PAX7, MYOGENIN, DYSTROPHIN, and MHC. Cells were then incubated for 1 hr at room temperature in the dark with secondary antibodies Alexa488-conjugated goat anti-mouse IgG and Alexa546-conjugated goat anti-mouse IgG (Invitrogen). Nuclei were stained with DAPI (Sigma-Aldrich). Images were taken on a Nikon Eclipse E800 or Leica DMI 4000B.

Immunohistochemistry and immunofluorescence of human muscle samples

Request a detailed protocol

Human muscle samples were fixed in neutral buffered formalin for 24 hr and then placed in 70 % ethanol for at least 24 hr. The sample with heterotopic bone was decalcified in 10 % EDTA (pH 7.2–7.4) before paraffin embedding and sectioning. Sections were stained with hematoxylin and eosin (J. David Gladstone Institutes Histology Core) or for alcian blue (pH 1.0) for cartilage and nuclear red stain for nuclei.

Freshly harvested human muscle was stored in DMEM with 30 % FBS at 4 °C, or snap in frozen in O.C.T. compound in 2-methylbutane chilled in liquid nitrogen. Serial 6 µm transverse frozen sections were analyzed or stored at –80 °C and processed similarly to the mouse TA samples above. Sections were stained with PAX7 (DSHB) and mouse monoclonal anti-Collagen Type I (Millipore-Sigma). Details about specimens are in the Figure 1—source data 1.

RT-PCR and quantitative analysis

Request a detailed protocol

Tissues were collected in TRI Reagent (Sigma-Aldrich) to isolate total RNA using the Arcturus PicoPure RNA kit (Applied Biosystems) as previously described for small samples (Schepers et al., 2012). A total of 0.2–0.5 µg of RNA were transcribed into cDNA with VeriScript cDNA synthesis kit (Affymetrix). cDNA was then pre-amplified with GE PreAmp Master Mix (Fluidigm Inc). Real-time quantitative PCR was performed in triplicated with either VeriQuest Probe qPCR Master Mix (Affymetrix) or Taqman Universal PCR Master Mix (Life Technologies) on either a Viia7 thermocycler (Life Technologies) or on a BioMark 48.48 dynamic array nanofluidic chip (Fluidigm, Inc) according to manufacturers’ instructions. ACTB was used for normalization as endogenous control.

Single cell RNA Sequencing and Analysis scRNAseq was performed using the Chromium Single Cell 3' Reagent Version 2 Kit from 10 X Genomics. 45,000 (FOP) and 30,000 (control) HNK1^-CD45^-CD31^- CXCR4⁺CD29⁺CD56^dim cells isolated from the iMPC differentiations were loaded onto one well of a 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 2 Kit. Each sample was sequenced on 1 lane of the NovaSeq 6,000 S4. Sequencing reads were processed with Cell Ranger version 2.0.0. using the human reference transcriptome GRCh38. 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 4—source data 1. Gene-barcoded matrices were analyzed with the R package Seurat v3.1.5 (Satija et al., 2015; Stuart et al., 2019; Team, 2014; Zheng et al., 2017). Gene core matrices from single cell RNA sequencing of primary human satellite cells isolated from a vastus muscle (Barruet et al., 2020) was used when comparing the transcriptional profile of hiPS-derived HNK1^-CD45^-CD31^- CXCR4⁺CD29⁺CD56^dim cells. For the comparison with primary Hu-MuSCs, hiPSC-derived cell sequencing reads were re-aligned using the human reference transcriptome hg19. Cells with fewer than 500 genes, greater than 5,000 genes and genes expressed in fewer than five cells were not included in the downstream analyses. Cells with more than 10 % mitochondrial counts were filtered out. Samples 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 mitochondrial contamination regressed out with the ScaleData function. In addition, while 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 muscle stem cells, we regressed out differences between G2/M and S cell cycle stage. 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. We performed differential gene-expression utilizing Seurat v3’s FindMarkers function with default settings which utilizes the Wilcoxon rank-sum test to calculate adjusted p values for multiple comparisons. We used the CellCycleScoring function to assign score based on the expression of G2/M and S phase markers (Regev et al., 2017). Myogenic cells were further analyzed by sub-clustering using the subset function. We then use the FindNeighbors (dims = 15) and FindClusters (resolution = 0.4) functions on the myogenic cell subset of the merged hiPSCs samples only to identify sub-clusters corresponding to different myogenic states. To order the cells in pseudotime based on their transcriptional similarity we used Monocle 2.12. Variable genes from Seurat analysis were used as input and clusters were projected onto the minimum spanning tree after ordering. Gene expression patterns were plotted with plot_genes_branched_heatmap, plot_genes_branched_pseudotime, and plot_multiple_branches_pseudotime. The BEAM (branch expression analysis modeling) function was used to score gene significance in a branch-dependent manner. Cells were re-ordered using the orderCells function to set branch A (myocytes) in Figure 6H as the ‘root-state’. This allowed us to determine genes that were significantly branch dependent in branch A (mainly Hu-MuSCs) vs branch B (hiPS-derived cells) using the BEAM analysis.

Statistical analysis

Request a detailed protocol

The data were analyzed utilizing GraphPad Prism v.7 software (GraphPad) using one-way (transplant) and two-way ANOVA with post hoc Tukey’s or Sidak’s multiple comparison test (gene expression). The Sidak test was used when comparing means between control and FOP, and the Tukey test was used when means of both control and FOP were compared together with other groups for the gene expression data. For the transplantation studies, at least three mice were used per group. At least three biological replicates were performed for each experiment unless indicated otherwise. All error bars are depicted as standard deviation, p-values are (*p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).

Human specimen procurement

Request a detailed protocol

Human samples were collected through the UCSF Biospecimens and Skeletal Tissues for Rare and Orphan Disease Genetics (BSTROnG) Biobank, using protocols approved by the UCSF Institutional Review Board. All participants provided written consent. For the transplantation studies (Figure 1D–H and Figure 1—figure supplement 1), biopsy of the control subject was obtained from a 44 yo female healthy individual undergoing surgery at UCSF, and the muscle from the FOP patient was obtained at autopsy from a 55 yo female. Written informed consent was obtain from all subjects or their families.

Data availability

Sequencing data have been deposited in GEO under accession codes GSE151918. All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 4,5,6. The dataset used for the primary Hu-MuSCs can be found here: https://doi.org/10.7272/Q65X273X. Detailed scripts can be found at https://github.com/EmilieB12/FOP_muscle/tree/main (copy archived at https://archive.softwareheritage.org/swh:1:rev:2c30ff4ef2608c0a44360047e54b9ad2d3017d0c) or in Source code 1–4.

The following data sets were generated

(2020) NCBI Gene Expression Omnibus
ID GSE151918. Control and FOP Human iPSC derived myogenic cell single cell RNA sequencing.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE151918

The following previously published data sets were used

1. Pomerantz JH
2. Barruet E
(2020) Dryad Digital Repository
Functionally heterogeneous human satellite cells identified by single cell RNA sequencing.

https://doi.org/10.7272/Q65X273X

References

1. Agarwal S
2. Loder S
3. Brownley C
4. Cholok D
5. Mangiavini L
6. Li J
7. Breuler C
8. Sung HH
9. Li S
10. Ranganathan K
11. Peterson J
12. Tompkins R
13. Herndon D
14. Xiao W
15. Jumlongras D
16. Olsen BR
17. Davis TA
18. Mishina Y
19. Schipani E
20. Levi B
(2016) Inhibition of Hif1α prevents both trauma-induced and genetic heterotopic ossification
PNAS 113:E338–E347.

https://doi.org/10.1073/pnas.1515397113
- PubMed
- Google Scholar
1. Al Tanoury Z
2. Rao J
3. Tassy O
4. Gobert B
5. Gapon S
6. Garnier J-M
7. Wagner E
8. Hick A
9. Hall A
10. Gussoni E
11. Pourquié O
(2020) Differentiation of the human PAX7-positive myogenic precursors/satellite cell lineage in vitro
Development 147:e12.

https://doi.org/10.1242/dev.187344
- PubMed
- Google Scholar
1. Andersen DC
2. Laborda J
3. Baladron V
4. Kassem M
5. Sheikh SP
6. Jensen CH
(2013) Dual role of delta-like 1 homolog (DLK1) in skeletal muscle development and adult muscle regeneration
Development 140:3743–3753.

https://doi.org/10.1242/dev.095810
- PubMed
- Google Scholar
(2021) Heterotopic ossification post COVID-19: Report of two cases
Radiology Case Reports 16:404–409.

https://doi.org/10.1016/j.radcr.2020.12.002
- PubMed
- Google Scholar
1. Barruet E
2. Morales BM
3. Lwin W
4. White MP
5. Theodoris CV
6. Kim H
7. Urrutia A
8. Wong SA
9. Srivastava D
10. Hsiao EC
(2016) The ACVR1 R206H mutation found in fibrodysplasia ossificans progressiva increases human induced pluripotent stem cell-derived endothelial cell formation and collagen production through BMP-mediated SMAD1/5/8 signaling
Stem Cell Research & Therapy 7:e115.

https://doi.org/10.1186/s13287-016-0372-6
- PubMed
- Google Scholar
1. Barruet E
2. Morales BM
3. Cain CJ
4. Ton AN
5. Wentworth KL
6. Chan TV
7. Moody TA
8. Haks MC
9. Ottenhoff TH
10. Hellman J
11. Nakamura MC
12. Hsiao EC
(2018) NF-κB/MAPK activation underlies ACVR1-mediated inflammation in human heterotopic ossification
JCI Insight 3:e22.

https://doi.org/10.1172/jci.insight.122958
- PubMed
- Google Scholar
1. Barruet E
2. Garcia SM
3. Striedinger K
4. Wu J
5. Lee S
6. Byrnes L
7. Wong A
8. Xuefeng S
9. Tamaki S
10. Brack AS
11. Pomerantz JH
(2020) Functionally heterogeneous human satellite cells identified by single cell RNA sequencing
eLife 9:e51576.

https://doi.org/10.7554/eLife.51576
- PubMed
- Google Scholar
1. Becht E
2. McInnes L
3. Healy J
4. Dutertre CA
5. Kwok IWH
6. Ng LG
7. Ginhoux F
8. Newell EW
(2018) Dimensionality reduction for visualizing single-cell data using UMAP
Nature Biotechnology 1:e4314.

https://doi.org/10.1038/nbt.4314
- PubMed
- Google Scholar
(2019) Fibro-Adipogenic Progenitors Cross-Talk in Skeletal Muscle: The Social Network
Frontiers in Physiology 10:1074.

https://doi.org/10.3389/fphys.2019.01074
- PubMed
- Google Scholar
1. Bigot A
2. Duddy WJ
3. Ouandaogo ZG
4. Negroni E
5. Mariot V
6. Ghimbovschi S
7. Harmon B
8. Wielgosik A
9. Loiseau C
10. Devaney J
11. Dumonceaux J
12. Butler-Browne G
13. Mouly V
14. Duguez S
(2015) Age-Associated Methylation Suppresses SPRY1, Leading to a Failure of Re-quiescence and Loss of the Reserve Stem Cell Pool in Elderly Muscle
Cell Reports 13:1172–1182.

https://doi.org/10.1016/j.celrep.2015.09.067
- PubMed
- Google Scholar
1. Billings PC
2. Fiori JL
3. Bentwood JL
4. O’Connell MP
5. Jiao X
6. Nussbaum B
7. Caron RJ
8. Shore EM
9. Kaplan FS
(2008) Dysregulated BMP signaling and enhanced osteogenic differentiation of connective tissue progenitor cells from patients with fibrodysplasia ossificans progressiva (FOP)
Journal of Bone and Mineral Research 23:305–313.

https://doi.org/10.1359/jbmr.071030
- PubMed
- Google Scholar
(2013) Derivation and FACS-mediated purification of PAX3+/PAX7+ skeletal muscle precursors from human pluripotent stem cells
Stem Cell Reports 1:620–631.

https://doi.org/10.1016/j.stemcr.2013.10.007
- PubMed
- Google Scholar
(2018) Mast cell inhibition as a therapeutic approach in fibrodysplasia ossificans progressiva (FOP)
Bone 109:259–266.

https://doi.org/10.1016/j.bone.2017.08.023
- PubMed
- Google Scholar
1. Cai J
2. Orlova VV
3. Cai X
4. Eekhoff EMW
5. Zhang K
6. Pei D
7. Pan G
8. Mummery CL
9. Ten Dijke P
(2015) Induced Pluripotent Stem Cells to Model Human Fibrodysplasia Ossificans Progressiva
Stem Cell Reports 5:963–970.

https://doi.org/10.1016/j.stemcr.2015.10.020
- PubMed
- Google Scholar
1. Chal J
2. Al Tanoury Z
3. Hestin M
4. Gobert B
5. Aivio S
6. Hick A
7. Cherrier T
8. Nesmith AP
9. Parker KK
10. Pourquié O
(2016) Generation of human muscle fibers and satellite-like cells from human pluripotent stem cells in vitro
Nature Protocols 11:1833–1850.

https://doi.org/10.1038/nprot.2016.110
- PubMed
- Google Scholar
1. Chal J
2. Al Tanoury Z
3. Oginuma M
4. Moncuquet P
5. Gobert B
6. Miyanari A
7. Tassy O
8. Guevara G
9. Hubaud A
10. Bera A
11. Sumara O
12. Garnier J-M
13. Kennedy L
14. Knockaert M
15. Gayraud-Morel B
16. Tajbakhsh S
17. Pourquié O
(2018) Recapitulating early development of mouse musculoskeletal precursors of the paraxial mesoderm in vitro
Development 145:dev157339.

https://doi.org/10.1242/dev.157339
- PubMed
- Google Scholar
1. Choi IY
2. Lim H
3. Estrellas K
4. Mula J
5. Cohen TV
6. Zhang Y
7. Donnelly CJ
8. Richard J-P
9. Kim YJ
10. Kim H
11. Kazuki Y
12. Oshimura M
13. Li HL
14. Hotta A
15. Rothstein J
16. Maragakis N
17. Wagner KR
18. Lee G
(2016) Concordant but Varied Phenotypes among Duchenne Muscular Dystrophy Patient-Specific Myoblasts Derived using a Human iPSC-Based Model
Cell Reports 15:2301–2312.

https://doi.org/10.1016/j.celrep.2016.05.016
- PubMed
- Google Scholar
1. Convente MR
2. Chakkalakal SA
3. Yang E
4. Caron RJ
5. Zhang D
6. Kambayashi T
7. Kaplan FS
8. Shore EM
(2018) Depletion of Mast Cells and Macrophages Impairs Heterotopic Ossification in an Acvr1R206H Mouse Model of Fibrodysplasia Ossificans Progressiva
Journal of Bone and Mineral Research 33:269–282.

https://doi.org/10.1002/jbmr.3304
- PubMed
- Google Scholar
1. Cosgrove BD
2. Gilbert PM
3. Porpiglia E
4. Mourkioti F
5. Lee SP
6. Corbel SY
7. Llewellyn ME
8. Delp SL
9. Blau HM
(2014) Rejuvenation of the muscle stem cell population restores strength to injured aged muscles
Nature Medicine 20:255–264.

https://doi.org/10.1038/nm.3464
- PubMed
- Google Scholar
(2014) Alk2 regulates early chondrogenic fate in fibrodysplasia ossificans progressiva heterotopic endochondral ossification
Stem Cells 32:1289–1300.

https://doi.org/10.1002/stem.1633
- PubMed
- Google Scholar
1. Dey D
2. Bagarova J
3. Hatsell SJ
4. Armstrong KA
5. Huang L
6. Ermann J
7. Vonner AJ
8. Shen Y
9. Mohedas AH
10. Lee A
11. Eekhoff EMW
12. van Schie A
13. Demay MB
14. Keller C
15. Wagers AJ
16. Economides AN
17. Yu PB
(2016) Two tissue-resident progenitor lineages drive distinct phenotypes of heterotopic ossification
Science Translational Medicine 8:366ra163.

https://doi.org/10.1126/scitranslmed.aaf1090
- PubMed
- Google Scholar
1. Elsafadi M
2. Manikandan M
3. Dawud RA
4. Alajez NM
5. Hamam R
6. Alfayez M
7. Kassem M
8. Aldahmash A
9. Mahmood A
(2016) Transgelin is a TGFβ-inducible gene that regulates osteoblastic and adipogenic differentiation of human skeletal stem cells through actin cytoskeleston organization
Cell Death & Disease 7:e2321.

https://doi.org/10.1038/cddis.2016.196
- PubMed
- Google Scholar
1. Farup J
2. Madaro L
3. Puri PL
4. Mikkelsen UR
(2015) Interactions between muscle stem cells, mesenchymal-derived cells and immune cells in muscle homeostasis, regeneration and disease
Cell Death & Disease 6:e1830.

https://doi.org/10.1038/cddis.2015.198
- PubMed
- Google Scholar
(2011) BMP signaling balances proliferation and differentiation of muscle satellite cell descendants
BMC Cell Biology 12:26.

https://doi.org/10.1186/1471-2121-12-26
- PubMed
- Google Scholar
Book
1. Furrer R
2. Handschin C
(2017)
Optimized Engagement of Macrophages and Satellite Cells in the Repair and Regeneration of Exercised Muscle

In: Spiegelman B, editors. Hormones, Metabolism and the Benefits of Exercise. Chamcham: Furrer R. pp. 57–66.
- Google Scholar
(2017) Human Satellite Cell Isolation and Xenotransplantation
Methods in Molecular Biology 1668:105–123.

https://doi.org/10.1007/978-1-4939-7283-8_8
- PubMed
- Google Scholar
1. Garcia Steven M
2. Tamaki S
3. Lee S
4. Wong A
5. Jose A
6. Dreux J
7. Kouklis G
8. Sbitany H
9. Seth R
10. Knott PD
11. Heaton C
12. Ryan WR
13. Kim EA
14. Hansen SL
15. Hoffman WY
16. Pomerantz JH
(2018) High-Yield Purification, Preservation, and Serial Transplantation of Human Satellite Cells
Stem Cell Reports 10:1160–1174.

https://doi.org/10.1016/j.stemcr.2018.01.022
- PubMed
- Google Scholar
1. Goel AJ
2. Rieder MK
3. Arnold HH
4. Radice GL
5. Krauss RS
(2017) Niche Cadherins Control the Quiescence-to-Activation Transition in Muscle Stem Cells
Cell Reports 21:2236–2250.

https://doi.org/10.1016/j.celrep.2017.10.102
- PubMed
- Google Scholar
1. Hatsell SJ
2. Idone V
3. Wolken DMA
4. Huang L
5. Kim HJ
6. Wang L
7. Wen X
8. Nannuru KC
9. Jimenez J
10. Xie L
11. Das N
12. Makhoul G
13. Chernomorsky R
14. D’Ambrosio D
15. Corpina RA
16. Schoenherr CJ
17. Feeley K
18. Yu PB
19. Yancopoulos GD
20. Murphy AJ
21. Economides AN
(2015) ACVR1R206H receptor mutation causes fibrodysplasia ossificans progressiva by imparting responsiveness to activin A
Science Translational Medicine 7:303ra137.

https://doi.org/10.1126/scitranslmed.aac4358
- PubMed
- Google Scholar
1. Haupt J
2. Stanley A
3. McLeod CM
4. Cosgrove BD
5. Culbert AL
6. Wang L
7. Mourkioti F
8. Mauck RL
9. Shore EM
(2019) ACVR1R206H FOP mutation alters mechanosensing and tissue stiffness during heterotopic ossification
Molecular Biology of the Cell 30:17–29.

https://doi.org/10.1091/mbc.E18-05-0311
- PubMed
- Google Scholar
1. Hicks MR
2. Hiserodt J
3. Paras K
4. Fujiwara W
5. Eskin A
6. Jan M
7. Xi H
8. Young CS
9. Evseenko D
10. Nelson SF
11. Spencer MJ
12. Handel BV
13. Pyle AD
(2018) ERBB3 and NGFR mark a distinct skeletal muscle progenitor cell in human development and hPSCs
Nature Cell Biology 20:46–57.

https://doi.org/10.1038/s41556-017-0010-2
- PubMed
- Google Scholar
1. Hino K
2. Ikeya M
3. Horigome K
4. Matsumoto Y
5. Ebise H
6. Nishio M
7. Sekiguchi K
8. Shibata M
9. Nagata S
10. Matsuda S
11. Toguchida J
(2015) Neofunction of ACVR1 in fibrodysplasia ossificans progressiva
PNAS 112:15438–15443.

https://doi.org/10.1073/pnas.1510540112
- PubMed
- Google Scholar
1. Hwang AB
2. Brack AS
(2018) Muscle Stem Cells and Aging
Current Topics in Developmental Biology 126:299–322.

https://doi.org/10.1016/bs.ctdb.2017.08.008
- PubMed
- Google Scholar
(2007) Asymmetric self-renewal and commitment of satellite stem cells in muscle
Cell 129:999–1010.

https://doi.org/10.1016/j.cell.2007.03.044
- PubMed
- Google Scholar
(2018) Activin-dependent signaling in fibro/adipogenic progenitors causes fibrodysplasia ossificans progressiva
Nature Communications 9:471.

https://doi.org/10.1038/s41467-018-02872-2
- PubMed
- Google Scholar
1. Lukjanenko L
2. Karaz S
3. Stuelsatz P
4. Gurriaran-Rodriguez U
5. Michaud J
6. Dammone G
7. Sizzano F
8. Mashinchian O
9. Ancel S
10. Migliavacca E
11. Liot S
12. Jacot G
13. Metairon S
14. Raymond F
15. Descombes P
16. Palini A
17. Chazaud B
18. Rudnicki MA
19. Bentzinger CF
20. Feige JN
(2019) Aging Disrupts Muscle Stem Cell Function by Impairing Matricellular WISP1 Secretion from Fibro-Adipogenic Progenitors
Cell Stem Cell 24:433–446.

https://doi.org/10.1016/j.stem.2018.12.014
- PubMed
- Google Scholar
1. Magli A
2. Incitti T
3. Kiley J
4. Swanson SA
5. Darabi R
6. Rinaldi F
7. Selvaraj S
8. Yamamoto A
9. Tolar J
10. Yuan C
11. Stewart R
12. Thomson JA
13. Perlingeiro RCR
(2017) PAX7 Targets, CD54, Integrin α9β1, and SDC2, Allow Isolation of Human ESC/iPSC-Derived Myogenic Progenitors
Cell Reports 19:2867–2877.

https://doi.org/10.1016/j.celrep.2017.06.005
- PubMed
- Google Scholar
1. Matsumoto Y
2. Hayashi Y
3. Schlieve CR
4. Ikeya M
5. Kim H
6. Nguyen TD
7. Sami S
8. Baba S
9. Barruet E
10. Nasu A
11. Asaka I
12. Otsuka T
13. Yamanaka S
14. Conklin BR
15. Toguchida J
16. Hsiao EC
(2013) Induced pluripotent stem cells from patients with human fibrodysplasia ossificans progressiva show increased mineralization and cartilage formation
Orphanet Journal of Rare Diseases 8:190.

https://doi.org/10.1186/1750-1172-8-190
- PubMed
- Google Scholar
1. Matsuo K
2. Chavez RD
3. Barruet E
4. Hsiao EC
(2019) Inflammation in Fibrodysplasia Ossificans Progressiva and Other Forms of Heterotopic Ossification
Current Osteoporosis Reports 17:387–394.

https://doi.org/10.1007/s11914-019-00541-x
- PubMed
- Google Scholar
1. Matsuo Koji
2. Lepinski A
3. Chavez RD
4. Barruet E
5. Pereira A
6. Moody TA
7. Ton AN
8. Sharma A
9. Hellman J
10. Tomoda K
11. Nakamura MC
12. Hsiao EC
(2021) ACVR1R206H extends inflammatory responses in human induced pluripotent stem cell-derived macrophages
Bone 153:e116129.

https://doi.org/10.1016/j.bone.2021.116129
- PubMed
- Google Scholar
1. Mauro A
(1961) Satellite cell of skeletal muscle fibers
The Journal of Biophysical and Biochemical Cytology 9:493–495.

https://doi.org/10.1083/jcb.9.2.493
- PubMed
- Google Scholar
1. McMahon AP
2. Bradley A
(1990) The Wnt-1 (int-1) proto-oncogene is required for development of a large region of the mouse brain
Cell 62:1073–1085.

https://doi.org/10.1016/0092-8674(90)90385-r
- PubMed
- Google Scholar
1. Medici D
2. Shore EM
3. Lounev VY
4. Kaplan FS
5. Kalluri R
6. Olsen BR
(2010) Conversion of vascular endothelial cells into multipotent stem-like cells
Nature Medicine 16:1400–1406.

https://doi.org/10.1038/nm.2252
- PubMed
- Google Scholar
1. Meyers C
2. Lisiecki J
3. Miller S
4. Levin A
5. Fayad L
6. Ding C
7. Sono T
8. McCarthy E
9. Levi B
10. James AW
(2019) Heterotopic Ossification: A Comprehensive Review
JBMR Plus 3:e10172.

https://doi.org/10.1002/jbm4.10172
- PubMed
- Google Scholar
(2012) Biglycan: a multivalent proteoglycan providing structure and signals
The Journal of Histochemistry and Cytochemistry 60:963–975.

https://doi.org/10.1369/0022155412456380
- PubMed
- Google Scholar
(2014) Lumican affects tumor cell functions, tumor-ECM interactions, angiogenesis and inflammatory response
Matrix Biology 35:206–214.

https://doi.org/10.1016/j.matbio.2013.09.003
- PubMed
- Google Scholar
1. Ono Y
2. Calhabeu F
3. Morgan JE
4. Katagiri T
5. Amthor H
6. Zammit PS
(2011) BMP signalling permits population expansion by preventing premature myogenic differentiation in muscle satellite cells
Cell Death and Differentiation 18:222–234.

https://doi.org/10.1038/cdd.2010.95
- PubMed
- Google Scholar
1. Regev A
2. Teichmann SA
3. Lander ES
4. Amit I
5. Benoist C
6. Birney E
7. Bodenmiller B
8. Campbell P
9. Carninci P
10. Clatworthy M
11. Clevers H
12. Deplancke B
13. Dunham I
14. Eberwine J
15. Eils R
16. Enard W
17. Farmer A
18. Fugger L
19. Göttgens B
20. Hacohen N
21. Haniffa M
22. Hemberg M
23. Kim S
24. Klenerman P
25. Kriegstein A
26. Lein E
27. Linnarsson S
28. Lundberg E
29. Lundeberg J
30. Majumder P
31. Marioni JC
32. Merad M
33. Mhlanga M
34. Nawijn M
35. Netea M
36. Nolan G
37. Pe’er D
38. Phillipakis A
39. Ponting CP
40. Quake S
41. Reik W
42. Rozenblatt-Rosen O
43. Sanes J
44. Satija R
45. Schumacher TN
46. Shalek A
47. Shapiro E
48. Sharma P
49. Shin JW
50. Stegle O
51. Stratton M
52. Stubbington MJT
53. Theis FJ
54. Uhlen M
55. van Oudenaarden A
56. Wagner A
57. Watt F
58. Weissman J
59. Wold B
60. Xavier R
61. Yosef N
62. Human Cell Atlas Meeting Participants
(2017) The Human Cell Atlas
eLife 6:e27041.

https://doi.org/10.7554/eLife.27041
- PubMed
- Google Scholar
1. Satija R
2. Farrell JA
3. Gennert D
4. Schier AF
5. Regev A
(2015) Spatial reconstruction of single-cell gene expression data
Nature Biotechnology 33:495–502.

https://doi.org/10.1038/nbt.3192
- PubMed
- Google Scholar
1. Schepers K
2. Hsiao EC
3. Garg T
4. Scott MJ
5. Passegué E
(2012) Activated Gs signaling in osteoblastic cells alters the hematopoietic stem cell niche in mice
Blood 120:3425–3435.

https://doi.org/10.1182/blood-2011-11-395418
- PubMed
- Google Scholar
(2000) Pax7 is required for the specification of myogenic satellite cells
Cell 102:777–786.

https://doi.org/10.1016/s0092-8674(00)00066-0
- PubMed
- Google Scholar
(2016) Regulation of Muscle Stem Cell Functions: A Focus on the p38 MAPK Signaling Pathway
Frontiers in Cell and Developmental Biology 4:e91.

https://doi.org/10.3389/fcell.2016.00091
- PubMed
- Google Scholar
1. Shea KL
2. Xiang W
3. LaPorta VS
4. Licht JD
5. Keller C
6. Basson MA
7. Brack AS
(2010) Sprouty1 regulates reversible quiescence of a self-renewing adult muscle stem cell pool during regeneration
Cell Stem Cell 6:117–129.

https://doi.org/10.1016/j.stem.2009.12.015
- PubMed
- Google Scholar
1. Shelton M
2. Metz J
3. Liu J
4. Carpenedo RL
5. Demers SP
6. Stanford WL
7. Skerjanc IS
(2014) Derivation and expansion of PAX7-positive muscle progenitors from human and mouse embryonic stem cells
Stem Cell Reports 3:516–529.

https://doi.org/10.1016/j.stemcr.2014.07.001
- PubMed
- Google Scholar
(2013) Overactive bone morphogenetic protein signaling in heterotopic ossification and Duchenne muscular dystrophy
Cellular and Molecular Life Sciences 70:407–423.

https://doi.org/10.1007/s00018-012-1054-x
- PubMed
- Google Scholar
1. Shore E. M
(2012) Fibrodysplasia ossificans progressiva: a human genetic disorder of extraskeletal bone formation, or--how does one tissue become another?
Wiley Interdisciplinary Reviews. Developmental Biology 1:153–165.

https://doi.org/10.1002/wdev.9
- PubMed
- Google Scholar
1. Shore Eileen M
2. Xu M
3. Feldman GJ
4. Fenstermacher DA
5. Cho T-J
6. Choi IH
7. Connor JM
8. Delai P
9. Glaser DL
10. LeMerrer M
11. Morhart R
12. Rogers JG
13. Smith R
14. Triffitt JT
15. Urtizberea JA
16. Zasloff M
17. Brown MA
18. Kaplan FS
(2006) A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva
Nature Genetics 38:525–527.

https://doi.org/10.1038/ng1783
- PubMed
- Google Scholar
1. Shultz LD
2. Lyons BL
3. Burzenski LM
4. Gott B
5. Chen X
6. Chaleff S
7. Kotb M
8. Gillies SD
9. King M
10. Mangada J
11. Greiner DL
12. Handgretinger R
(2005) Human lymphoid and myeloid cell development in NOD/LtSz-scid IL2R gamma null mice engrafted with mobilized human hemopoietic stem cells
Journal of Immunology 174:6477–6489.

https://doi.org/10.4049/jimmunol.174.10.6477
- PubMed
- Google Scholar
(2016) Age-related changes in miR-143-3p:Igfbp5 interactions affect muscle regeneration
Aging Cell 15:361–369.

https://doi.org/10.1111/acel.12442
- PubMed
- Google Scholar
1. Spencer CI
2. Baba S
3. Nakamura K
4. Hua EA
5. Sears MAF
6. Fu C
7. Zhang J
8. Balijepalli S
9. Tomoda K
10. Hayashi Y
11. Lizarraga P
12. Wojciak J
13. Scheinman MM
14. Aalto-Setälä K
15. Makielski JC
16. January CT
17. Healy KE
18. Kamp TJ
19. Yamanaka S
20. Conklin BR
(2014) Calcium transients closely reflect prolonged action potentials in iPSC models of inherited cardiac arrhythmia
Stem Cell Reports 3:269–281.

https://doi.org/10.1016/j.stemcr.2014.06.003
- PubMed
- Google Scholar
1. Stanley A
2. Heo SJ
3. Mauck RL
4. Mourkioti F
5. Shore EM
(2019) Elevated BMP and Mechanical Signaling Through YAP1/RhoA Poises FOP Mesenchymal Progenitors for Osteogenesis
Journal of Bone and Mineral Research 34:1894–1909.

https://doi.org/10.1002/jbmr.3760
- PubMed
- Google Scholar
1. Stantzou A
2. Schirwis E
3. Swist S
4. Alonso-Martin S
5. Polydorou I
6. Zarrouki F
7. Mouisel E
8. Beley C
9. Julien A
10. Le Grand F
11. Garcia L
12. Colnot C
13. Birchmeier C
14. Braun T
15. Schuelke M
16. Relaix F
17. Amthor H
(2017) BMP signaling regulates satellite cell-dependent postnatal muscle growth
Development 144:2737–2747.

https://doi.org/10.1242/dev.144089
- PubMed
- Google Scholar
1. Stuart T
2. Butler A
3. Hoffman P
4. Hafemeister C
5. Papalexi E
6. Mauck WM
7. Hao Y
8. Stoeckius M
9. Smibert P
10. Satija R
(2019) Comprehensive Integration of Single-Cell Data
Cell 177:1888–1902.

https://doi.org/10.1016/j.cell.2019.05.031
- PubMed
- Google Scholar
1. Sullivan MP
2. Torres SJ
3. Mehta S
4. Ahn J
(2013) Heterotopic ossification after central nervous system trauma: A current review
Bone & Joint Research 2:51–57.

https://doi.org/10.1302/2046-3758.23.2000152
- PubMed
- Google Scholar
1. Takahashi K
2. Tanabe K
3. Ohnuki M
4. Narita M
5. Ichisaka T
6. Tomoda K
7. Yamanaka S
(2007) Induction of pluripotent stem cells from adult human fibroblasts by defined factors
Cell 131:861–872.

https://doi.org/10.1016/j.cell.2007.11.019
- PubMed
- Google Scholar
Software
1. Team RC
(2014)
A language and environment for statistical computing, version 2.6.2

R:Development Core Team.
1. Trapnell C
2. Cacchiarelli D
3. Grimsby J
4. Pokharel P
5. Li S
6. Morse M
7. Lennon NJ
8. Livak KJ
9. Mikkelsen TS
10. Rinn JL
(2014) The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells
Nature Biotechnology 32:381–386.

https://doi.org/10.1038/nbt.2859
- PubMed
- Google Scholar
1. Urist MR
(1965) Bone: formation by autoinduction
Science 150:893–899.

https://doi.org/10.1126/science.150.3698.893
- PubMed
- Google Scholar
(2018) Large-Scale Expansion of Human iPSC-Derived Skeletal Muscle Cells for Disease Modeling and Cell-Based Therapeutic Strategies
Stem Cell Reports 10:1975–1990.

https://doi.org/10.1016/j.stemcr.2018.04.002
- PubMed
- Google Scholar
1. Volpato V
2. Webber C
(2020) Addressing variability in iPSC-derived models of human disease: guidelines to promote reproducibility
Disease Models & Mechanisms 13:dmm042317.

https://doi.org/10.1242/dmm.042317
- PubMed
- Google Scholar
1. Wosczyna MN
2. Rando TA
(2018) A Muscle Stem Cell Support Group: Coordinated Cellular Responses in Muscle Regeneration
Developmental Cell 46:135–143.

https://doi.org/10.1016/j.devcel.2018.06.018
- PubMed
- Google Scholar
1. Wu J
2. Matthias N
3. Lo J
4. Ortiz-Vitali JL
5. Shieh AW
6. Wang SH
7. Darabi R
(2018) A Myogenic Double-Reporter Human Pluripotent Stem Cell Line Allows Prospective Isolation of Skeletal Muscle Progenitors
Cell Reports 25:1966–1981.

https://doi.org/10.1016/j.celrep.2018.10.067
- PubMed
- Google Scholar
1. Xi H
2. Fujiwara W
3. Gonzalez K
4. Jan M
5. Liebscher S
6. Van Handel B
7. Schenke-Layland K
8. Pyle AD
(2017) In Vivo Human Somitogenesis Guides Somite Development from hPSCs
Cell Reports 18:1573–1585.

https://doi.org/10.1016/j.celrep.2017.01.040
- PubMed
- Google Scholar
1. Xi H
2. Langerman J
3. Sabri S
4. Chien P
5. Young CS
6. Younesi S
7. Hicks M
8. Gonzalez K
9. Fujiwara W
10. Marzi J
11. Liebscher S
12. Spencer M
13. Van Handel B
14. Evseenko D
15. Schenke-Layland K
16. Plath K
17. Pyle AD
(2020) A Human Skeletal Muscle Atlas Identifies the Trajectories of Stem and Progenitor Cells across Development and from Human Pluripotent Stem Cells
Cell Stem Cell 27:158–176.

https://doi.org/10.1016/j.stem.2020.04.017
- PubMed
- Google Scholar
1. Xu X
2. Wilschut KJ
3. Kouklis G
4. Tian H
5. Hesse R
6. Garland C
7. Sbitany H
8. Hansen S
9. Seth R
10. Knott PD
11. Hoffman WY
12. Pomerantz JH
(2015) Human Satellite Cell Transplantation and Regeneration from Diverse Skeletal Muscles
Stem Cell Reports 5:419–434.

https://doi.org/10.1016/j.stemcr.2015.07.016
- PubMed
- Google Scholar
1. Zheng GXY
2. Terry JM
3. Belgrader P
4. Ryvkin P
5. Bent ZW
6. Wilson R
7. Ziraldo SB
8. Wheeler TD
9. McDermott GP
10. Zhu J
11. Gregory MT
12. Shuga J
13. Montesclaros L
14. Underwood JG
15. Masquelier DA
16. Nishimura SY
17. Schnall-Levin M
18. Wyatt PW
19. Hindson CM
20. Bharadwaj R
21. Wong A
22. Ness KD
23. Beppu LW
24. Deeg HJ
25. McFarland C
26. Loeb KR
27. Valente WJ
28. Ericson NG
29. Stevens EA
30. Radich JP
31. Mikkelsen TS
32. Hindson BJ
33. Bielas JH
(2017) Massively parallel digital transcriptional profiling of single cells
Nature Communications 8:14049.

https://doi.org/10.1038/ncomms14049
- PubMed
- Google Scholar

Article and author information

Author details

Emilie Barruet
1. Division of Plastic and Reconstructive Surgery, Departments of Surgery and Orofacial Sciences, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States
2. Division of Endocrinology and Metabolism, Department of Medicine, the Institute for Human Genetics, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States
3. Division of Plastic and Reconstructive Surgery, Department of Surgery, University of California - San Francisco, San Francisco, United States
Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft, Writing – review and editing

For correspondence
emilie.barruet@ucsf.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4593-024X
Steven M Garcia

Division of Plastic and Reconstructive Surgery, Departments of Surgery and Orofacial Sciences, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States

Contribution
Investigation, Validation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-7833-6677
Jake Wu

Division of Plastic and Reconstructive Surgery, Departments of Surgery and Orofacial Sciences, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States

Contribution
Investigation
Blanca M Morales

Division of Endocrinology and Metabolism, Department of Medicine, the Institute for Human Genetics, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States

Contribution
Investigation

Competing interests
No competing interests declared
Stanley Tamaki

Division of Plastic and Reconstructive Surgery, Departments of Surgery and Orofacial Sciences, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States

Contribution
Investigation

Competing interests
No competing interests declared
Tania Moody

Division of Endocrinology and Metabolism, Department of Medicine, the Institute for Human Genetics, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States

Contribution
Data curation, Software

Competing interests
No competing interests declared
Jason H Pomerantz
1. Division of Plastic and Reconstructive Surgery, Departments of Surgery and Orofacial Sciences, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States
2. Division of Plastic and Reconstructive Surgery, Department of Surgery, University of California - San Francisco, San Francisco, United States
Contribution
Clinical samples, Clinical samples, Conceptualization, Funding acquisition, Resources, Supervision, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5107-1883
Edward C Hsiao

Division of Endocrinology and Metabolism, Department of Medicine, the Institute for Human Genetics, the Program in Craniofacial Biology, and the Eli and Edythe Broad Center of Regeneration Medicine, University of California - San Francisco, San Francisco, United States

Contribution
Clinical samples, Clinical samples, Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review and editing

For correspondence
Edward.Hsiao@ucsf.edu

Competing interests
ECH receives clinical trial research funding from Clementia Pharmaceuti-cals, an Ipsen company, and Neurocrine Biosciences, Inc., through his institution. ECH received prior funding from Regeneron Pharmaceuticals, through his institution. ECH serves in an unpaid capacity on the international FOP Association Medical Registry Advisory Board, on the International Clinical Council on FOP, and on the Fibrous Dysplasia Foundation Medical Advisory Board. These activities pose no conflicts for the presented research

"This ORCID iD identifies the author of this article:" 0000-0001-8924-106X

Funding

National Institute of Arthritis and Musculoskeletal and Skin Diseases (R01AR066735)

Edward C Hsiao

French Muscular Dystrophy Association (Trampoline grant)

Emilie Barruet
Edward C Hsiao

Radiant Hope Foundation

Edward C Hsiao

University of California, San Francisco (UCSF Cohort Development Grant)

Edward C Hsiao

California Institute for Regenerative Medicine (TG2-01153)

Emilie Barruet

University of California, San Francisco (UCSF Program for Breakthrough Biomedical Research)

Emilie Barruet

National Institute of Arthritis and Musculoskeletal and Skin Diseases (R01AR072638-03)

Jason H Pomerantz

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors thank Kelly Wentworth and Samuel Kou for their assistance collecting the autopsy samples and Francesco Tedesco for his technical support and discussion on the AFM grant. This work was supported by a NIH/NIAMS R01AR066735 to ECH, a French Muscular Association (AFM-Telethon) Trampoline grant to ECH and EB, the Radiant Hope Foundation to ECH, and the UCSF Cohort Development Grant to ECH; the California Institute for Regenerative Medicine Fellowship Program to UCSF (TG2-01153) to EB, and the UCSF Program for Breakthrough Biomedical Research (PBBR) to EB; and the NIH R01AR072638-03 to JHP. Finally, the authors would like to thank the patients and their families for their generous specimen donations.

Ethics

Human subjects: Human samples were collected through the University of California - San Francisco Biospecimens and Skeletal Tissues for Rare and Orphan Disease Genetics (BSTROnG) Biobank, using protocols approved by the UCSF Institutional Review Board (10-03053 and 11-06711). All participants provided written consent, which includes consent to collect, use, and publish research data.

All mouse studies were performed using protocols approved by the UCSF Institutional Animal Care and Use Committee. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (AN181101-02A) of the University of California, San Francisco.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.