Skeletal muscle organoid protocol and correlation to fetal development.

(A) Brightfield images of myogenic development stages, with corresponding cytokines/growth factors protocol applications (B) Graph depicting organoid development in size (Day 0, n= 51; Day 3, n= 46; Day 5, n= 48; Day 7, n= 50; Day 11, n= 50; Day 17, n= 33; for each timepoint organoids from 3 independent derivations were measured) (C) Representative mesodermal T+, neuromesodermal T+, SOX2+ and paraxial mesodermal TBX6+ organoid expression (in live DAPI-cells) at Day 5 (D) Representative organoid PAX3+, mesodermal SOX2- and neural SOX2+ expression at Day 7; PAX3/PAX7 coexpression at Day 11 (E) Graph depicting qPCR values for anterior PAX3/UNCX/MEOX2/TBX18 and posterior TBX6/HES7 somitic mesodermal markers (F) Graph depicting qPCR values for epaxial and hypaxial dermomyotomal PAX7/EN1/SIM1/LBX1 and neural crest TFAP2A/SOX10 markers (G) Representative organoid neural plate border epithelial PAX3+ / PAX7+ / SOX2+ / TFAP2A+, paraxial mesodermal PAX3+/SOX2- and delaminating specified neural crest progenitor PAX3+ / SOX10+ expression at Day 17 (H) Representative organoid myogenic FastMyHC+, PAX7+ and neural SOX2+/TFAP2A+/SOX10+ expression at Day 23. Dashed line represents location of embryoid body embedded into Matrigel. Statistics: Values at each timepoint represent difference in mean relative expression for each gene (D7= Day 5 - Day 7, D11= Day 7 - Day 11, D17= Day 11 - Day 17, D23= Day 17 - Day 23) as derived by performing ordinary one-way ANOVA and Tukey ′s multiple comparison tests (E,F). Scale bars: 200μm (G), 100μm (A,D,H), 50μm

Neural lineage development during skeletal muscle organoid progression.

(A) Stepwise brightfield images depicting delamination / migration of progenitor population during organoid culture progression and corresponding myofiber formation (B) Gene ontology enrichment analysis comparing 4 and 8 weeks organoids attributes muscle identity at 8 weeks post differentiation and highlights muscle system development and neural lineage arrest among the top upregulated and downregulated gene ontology terms, respectively (C-E) Organoid overview at Day 35 indicates predominant expression of FastMyHC+ and PAX7+ myogenic populations, while SOX2+ neural populations demarcate SOX2 neural plate border epithelium location as observed at earlier stages (Day 16) (C); PAX7 cells are of myogenic origin (PAX7+/SOX10-), MF20+ myotubes are in their proximity and MYOD1+ cells appear at organoid periphery (D); TUJ1+ neurons are restricted to inner organoid areas and close to SOX2+ epithelium, while FastMyHC+ myofibers occupy exterior organoid areas(E) (F) Histographs based on FACS intracellular quantification depicting percentage of PAX7+ or MYOD1+ cells through differentiation protocol. For each replicate 10 organoids were pooled, n=10 Statistics: *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, ns: not significant (G) Heatmap of HOX gene cluster emphasizes organoid culture limb axial anatomical identity, by depicting transition from an initial limb bud (HOX 9-10) towards a more distal identity (HOX 11-13) at 8 weeks and 16 weeks post differentiation, respectively. Scale bars, 500μm (C), 200μm (A, D, E)

Skeletal muscle organoid characterization at single cell resolution.

(A) t-SNE visualization of color-coded clustering (n= 4.323 cells) at 12 weeks post differentiation highlights predominant presence of skeletal muscle lineage, represented by clusters corresponding to myogenic progenitors (n=1625 cells, 37% of total population) in non-dividing (n=1317 cells) and mitotic (n=308 cells) state, myoblasts (n=731 cells), myocytes (n=1147 cells) and myotubes (n=442). Additionally, mesenchymal and neural lineages are represented by 2 smaller clusters of fibroadipogenic (n=165 cells) and neural (n=213 cells) progenitors, respectively (B) Immunocytochemistry at Day 135 indicates derivation of PRDM16+ adipocyte clusters distinct from PAX7+ myogenic progenitors (C) Organoid overview at Day 98 depicts expression of PAX7+ myogenic populations and highlights SOX2+ neural populations towards organoid interior (D) Positive areas with Oil O Red staining indicate derivation of adipocytes in organoid culture at Day 135 (E) Violin plots of key markers PAX7, MYOD1, MYOG, MYH3 from each stage as in 3A depict relative expression levels and emphasize gradual transition from myogenic progenitor to myotube subcluster (F) Dot plot showing expression of representative genes related to adipogenesis and fibrogenesis across the 7 main clusters. Circle area represents percentage of gene+ cells in a cluster, color reflects average expression level (gray, low expression; blue, high expression) (G) Dot plot showing expression of representative genes related to fetal myogenesis across the 7 main clusters. Circle area represents percentage of gene+ cells in a cluster, color reflects average expression level (gray, low expression; blue, high expression). (H-J) Representative organoid overview at Day 98 indicates predominant expression of Fast MyHC fetal myofibers (H,J), positive for M-Cadherin and in PAX7+ cells proximity (I,J) (K) Representative recording (n=6) of acetylcholine (ACh)-induced changes in holding current in a single skeletal muscle cell. ACh (10 μM) was applied as indicated by the bar. Holding potential -90 mV. Downward deflections represent membrane currents in response to depolarizing voltage ramps (duration 500 ms) from -120 mV to +60 mV. Dashed line indicates zero current level. (L) I/V curve of nAChR currents obtained by subtraction of voltage ramp-induced changes of current in presence and absence of ACh (10-5 mol/L), corresponding lower-case letters b - a (M) Summarized FRET-recordings from skeletal muscle cells transfected with Twitch2B to monitor the increase in [Ca2+]i during ACh application (h) Scale bars, 1mm (C,H), 200μm (B), 100μm (I,J)

Myogenic progenitor identity and comparison to progenitors derived from fetal and adult muscle tissue.

(A) t-SNE plot visualization of color-coded clustering indicates myogenic progenitor sub-cluster with distinct molecular signatures: “dormant” PAX7high/ CHODLhigh/FBN1high, “activated” CD44high / CD98+ / MYOD1+ and “mitotic” KI-67+/CDK1+/TOP2A (B) Organoid overviews at Day 98 illustrate CD44 and PAX7 expressing cells at developing sites, which are more accessible to HGF activation signal, specificity of CD44 on MYOD+/PAX7+ progenitor expressing cells (arrows) and absence from FastMyHC+ positive myofibers is highlighted (C) FBN1+ microfibrils are located towards organoid interior (D) Pseudotime ordering for myogenic progenitors and myoblast corresponding clusters highlights distinct developmental trajectories promoting myogenic commitment and self-renewal Correlation coefficient plot for Log2 fold change (Log2 FC) values for isolated myogenic progenitors from human fetal tissue (17w) and FACS sorted CXCR4+ / ITGβ1+ organoid derived myogenic progenitors (16w). PAX7, COL1A1, COL3A1, COL4A1, COL5A1, COL15A1, FBN1 and CHODL and further extracellular matrix related genes are highlighted on the plot. Pearson’s correlation coefficient, rho=0.9 for global expression comparison and rho=0.97 for selected genes (E) UMAP color-based clustering divides non-dividing myogenic progenitors and adult satellite cells into 4 clusters with distinct molecular signatures: satellite cells characterized by SPRY1high/ZFP36high /MYF5high expression, co-clustered with dormant SPRY1high / FBN1high / CHODLhigh / PAX7high, activated CD44high / CD98high / MYOD1+ and committed, NEB high/ MYH3high / MYOD1high organoid derived myogenic progenitors. Dots correspond to adult satellite cells from GSE130646 database, triangles correspond to organoid derived myogenic progenitors. (F) Violin plots depicting relative expression levels of markers specific for quiescent PAX7 or activated MYF5 muscle stem cell state across adult satellite cells (GSE130646) and organoid derived myogenic progenitors subclusters (H) Ridge plots of developmental score distribution of myogenic progenitors across in vivo or in vitro stages based on the difference between up-regulated satellite cell and embryonic markers from human reference atlases for weeks (Wk) 5 to 18 embryonic and fetal stages, years (Yr) 7 to 42 adult satellite cells and SM (skeletal muscle) organoid (I,J) In vivo engraftment potential of human myogenic progenitors. 100.000 CD82+ sorted human cells were injected into Tibialis anterior muscle of nude mice (I) Control mice were not injected (J). Six weeks post transplantation, transverse cryosections of muscle fibers were stained with huLamin A/C (green), Dystrophin (red) and DAPI (blue). Human cells appear green and red in contrast to murine cells which only show a dystrophin positive red staining. Scale bars 200 μm in I,J

Lineage representation and organoid culture progression at early stages of differentiation protocol.

(A) qPCR analysis depicts relative expression of pluripotent (OCT4, SOX2, NANOG), mesodermal / posterior somitic (BRACHYURY,MESOGENIN,TBX6,HES7), anterior somitic (PAX3,UNCX,TBX18,MEOX2), dermomyotomal (PAX7, SIM1, EN1, LBX1) and neural tube / crest markers (SOX2, PAX6, TFAP2A, SOX10) (B) Gating strategies and FACS intracellular quantification of PAX3 and PAX7 markers at Day 10 and 16 (C) Tile scan overview at Day 23 shows neural crest SOX10+ populations close to neural plate border epithelium (D) and SOX10-/PAX3+ populations at more outer locations (E). For each replicate 6 organoids were pooled; gene expression was normalized to housekeeping gene. Statistics: *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, ns: not significant. FACS plots, yellow: unstained population, red: isotype control, gray: PAX7+ or PAX3+ population. Scale bars, 500uM (C), 200uM (D, E)

Myogenic versus neural fate during organoid development.

(A) Heatmap of limb myogenic progenitor markers emphasizes upregulation through organoid protocol development (B) Differential expression comparison highlights upregulated (n=5.919) and downregulated (n=5.005) genes between 8 vs 4 weeks organoids. (Padj< 0.001) (C) Heatmap of neural tube and neural crest markers during 4, 8 and 16 weeks post differentiation highlights neural lineage arrest from 4 weeks onwards (D) Organoid overview at Day 35 indicates presence of MF20+ myofibers in proximity of PAX7+ cells and delamination migration of mainly myogenic cells. SOX2 and TFAP2A expression is restricted towards inner organoid areas (E) Gating strategy applied to quantify MYOD1+ (committed state) and PAX7+ (progenitor state) populations at 5 and 8 weeks post differentiation. Yellow histogram: unstained population, red: isotype control, gray: PAX7+ or MYOD+ population. Scale bars, 500uM

Single cell RNA seq expression profiling and lineage representation in organoid culture at week 12

(A) Gene signatures of t-SNE described clusters based on relative expression levels of the 50 most significant markers for each of the 7 clusters (B) t-SNE plot of myogenic and fibroadipogenic markers (C) Brightfield image at Day 135 highlights presence of structures reassembling adipocytes (D) Immunocytochemistry at Day 135 emphasizing derivation of PRDM16+ cells in FastMyHC+ myofibers proximity (E) Organoid overview at Day 98 highlights muscle system development at outer and neural lineage representation at more inner sites (F) t-SNE plot of key markers PAX7, MYOD1, MYOG, MYH3 from each stage depicts relative expression levels and emphasize gradual transition from myogenic progenitor to myotube subcluster (G) Heatmap of embryonic or fetal myogenic markers depicts a fetal environment during 8 to16 weeks post differentiation (H) Two-dimensional principal component analysis based on gene expression profiling between samples separates skeletal muscle organoid cluster from pluripotent and highlights greater difference in variance between 4 and 8-16 weeks subgroups in comparison to 8 and 16 weeks subgroups. w.p.d.: weeks post differentiation. Scale bars, 1mM (D), 200uM (E, F).

Skeletal muscle organoid culture maturation and identity.

(A) Differential expression comparison highlights upregulated (n=2.856) and downregulated (n=2.369) genes between 16 and 8 weeks organoids (Padj<0.001) (B-C) Gene ontology enrichment analysis at 16 weeks depicting top statistically significant upregulated (B) and downregulated (C) gene ontology (GO) terms; (D) statistically significant upregulated GO terms further highlights skeletal muscle maturation (E) Gating strategy applied to quantify MYOD1+ (committed state) and PAX7+ (progenitor state) population at 14 w post differentiation (F) Organoid overview at Day 98 depicts high proportion of PAX7+ /MYOD+ cells towards outer areas. Yellow histogram: unstained population, Red histogram: Isotype Control, Gray Histogram: PAX7+ or MYOD+ population. Scale bars, 500uM (F)

Functional properties of organoid derived skeletal muscle myofibers.

(A) Second Harmonic Generation imaging (SHG) from areas with dense fiber network reveals myofibril formation with distinct sarcomeres (z plane of stack right) (B) Electron microscopy depicts myofibrils with well-developed sarcomeres (C,D) Skeletal muscle organoid derived myofibers expressed LAMININ and DYSTROPHIN (E) Representative recording of acetylcholine (Ach)-induced changes in holding current in a single skeletal muscle cell from organoid derived skeletal muscle cells from iPSC lines with Duchenne (DMD_IPS1, DMD_IPSCORE_65_1) and Becker (BMD_IPS1) muscular dystrophy patients. ACh (10 μM) was applied as indicated by the bar. Holding potential -90 mV. Downward deflections represent membrane currents in response to depolarizing voltage ramps (duration 500 ms) from -120 mV to +60 mV. Dashed line indicates zero current level. (F) Summarized FRET-recordings from organoid derived skeletal muscle cells from iPSC lines with Duchenne (DMD_IPS1, DMD_IPSCORE_65_1), Becker (BMD_IPS1) transfected with Twitch2B to monitor the increase in [Ca2+]i during ACh application. Scale Bars, 100uM (C,D), 20uM in (A), 500nM (B)

Subclustering of myogenic progenitors and NOTCH signaling.

(A) Dot plot showing expression of representative genes related to satellite cell identity across the 7 main clusters. (B) Dot plot showing expression of representative genes related to NOTCH signaling, which are identified as regulators of satellite cell quiescence across the 7 main clusters (C) Organoid immunohistochemistry at Day 98 indicates PAX7 myogenic progenitors that are Ki-67-/EdU-(D) Gating strategy applied to quantify EdU+ proliferating cells in organoid culture at Day 98 (14 weeks); Histograph depicting percentage of EdU+ cells; Dot plot, circle area represents percentage of gene+ cells in a cluster, color reflects average expression level (gray, low expression; blue, high expression). Scale Bars, 100uM in (C).

Pseudotime ordering of myogenic progenitor revealing distinct states and cell fate decisions.

(A) Gene signatures of t-SNE described clusters based on relative expression levels of the 50 most significant markers for each of the 3 clusters (B) Expression of selected genes along pseudotime. Group of genes selected for myogenic progenitors: PAX7, SPRY1, CHODL (dormant state), CD44, CD98, MYOD1 (activated state), TOP2A (mitotic state), and for myoblasts: MYOD1 (C, D) UMAP feature plots depicting relative expression of extracellular matrix proteins COL1A2, COL5A2, COL5A3, FN1 (C) and selected transcriptional regulators that define dormant (PAX7, FBN1, CHODL, SPRY1), activated (CD44, CD98, MYOD1, VEGFA) and mitotic (TOP2A) state of myogenic progenitor and myoblast (MYOD1) clusters (D). (E,F) UMAP feature plots depicting relative expression of genes regulating asymmetric divisions and self-renewal (PARD3, p38a/b, CD9, NOTCH3) (E) and extracellular matrix collagens (COL1A1, COL3A1, COL4A1, COL5A1, COL6A1, COL15A1; COL4A2, COL6A2 (F).

Organoid derived myogenic progenitors and correlation to fetal muscle progenitors.

(A) Gating strategy applied for FACS sorting CD29+ / CXCR4+ cells from 15 -16 weeks skeletal muscle organoids (B,C) Re-plating CD29+ / CXCR4+ cells and culturing for 14 days highlights fetal myogenic potential, illustrated by brightfield and immunocytochemistry images for Fast MyHC+, TITIN+ and PAX7+ populations (D) t-SNE feature plots for ITGβ1 and CXCR4 demarcate expression of FACS isolated calls into myogenic progenitor subcluster (E) Gating strategy from (A) together with unstained population (yellow), isotype control (red) and CD29+ / CXCR4+ (gray) population (F) Correlation coefficient plot for Log2 fold change (Log2 FC) values for isolated myogenic progenitors from human fetal tissue (17w) and FACS sorted CXCR4+ / ITGβ1+ organoid derived myogenic progenitors (16w). PAX7, MYF5, MYOD1, MYOG are highlighted on the plot (G) Correlation coefficient plot for Log2 fold change (Log2 FC) values for FACS sorted CXCR4+ / ITGβ1+ organoid derived myogenic progenitors (16 w) from CB CD34+, DMD_IPS1, BMD_iPS1 and iPSCORE_65_1 pluripotent lines. Scale bars, 200uM (B), 100uM (C).

Organoid derived myogenic progenitors and correlation to adult human satellite cells.

(A) Violin plots highlighting developmental status for organoid derived myogenic progenitors and satellite cells through relative expression of key signature markers ZFP36, CD44, FBN1, NEB and SPRY1 for each subcluster (B) Enhanced Volcano plot from comparing transcript levels between all adult satellite cells (GSE130646) and organoid derived myogenic progenitors. Log2 fold-change in normalized gene expression vs. - Log10 adjusted p-value is plotted. Differentially expressed genes: blue, adjusted p-value < 0.001; green, adjusted p-value > 0.001 and log2 fold-change > 1;,red, log2 fold-change > 1 and adjusted p value < 0.001; gray, no significance. (C) Pseudotime ordering for organoid derived myogenic progenitors and adult satellite cells together with expression of selected genes (D) Expression of genes selected for dormant myogenic progenitors (PAX7, NOTCH3, RBPJ) and activated satellite cells (MYF5, JUN, FOS). (E) UMAP feature plots depicting expression of PAX7, CD82+, MYF5 and NEB within the myogenic progenitor and adult satellite cell clusters corresponding to Figure 4F.

Characterization of cell-cell communication network for all clusters at week 12 of human skeletal muscle organoid development

(A) Circle plot illustrates the aggregated cell-cell communication network for all clusters at week 12 of human skeletal muscle organoids development. Circle sizes are proportional to the number of cells in each cell group and edge width represents the communication probability. (B) Heatmap indicating signals contributing most to the outgoing or incoming signalling among activated, mitotic, dormant/specification resistant myogenic progenitors and fibroadipogenic progenitors, neural progenitors and myofibers related clusters at week 12 of human skeletal muscle organoid development.

Reproducibility of organoid culture at early and mature stages.

(A) Pseudotemporal ordering of organoids from Day 2, 7 and 11 based on qPCR expression profiling of selected 32 genes indicates robust homogeneous development of organoids during WNT activation and BMP Inhibition, while the subsequent development of paraxial mesoderm and neural tube related lineages via bFGF, Retinoic Acid (Day 7) and WNT1a, hSHH (Day 11) stimulation, introduced small variations to the organoid culture system. (B) Integrative analysis between datasets of organoids from different lines demonstrate highly conserved lineage representation for all datasets at mature stages (week 12) of human skeletal muscle organoid development. Feature plots on PAX7, CD44, FBN1, MYH3, POSTN, SOX2 as well as NOTCH3, MYOG, TOP2A indicate presence of dormant, activated and mitotic myogenic progenitors together with skeletal muscle myofibers, fibroadipogenic progenitors and neural progenitors related clusters.

Comparison between 2D in vitro myogenic differentiation protocols and in vivo staging.

Violin plots highlighting expression of key myogenesis markers: Comparison between myogenic progenitors derived in 2D differentiation protocols (HX protocol, Xi et al., 2017; JC protocol, Chal et al., 2015; MS protocol, Shelton et al., 2014) and in vivo myogenic progenitors (5 stages: 1,2: embryonic; 3,4: fetal; 5: postnatal). Data were acquired from Pyle ′s LAB UCLA website, (https://aprilpylelab.com/datasets/in-vivo-and-in-vitro-multiple-protocols-smpcs-and-scs/) and Xi et al. (2020).

qPCR primer pairs applied to detect relative expression of key markers during skeletal muscle organoid development.

qPCR primer pairs applied for diffusion map analysis of early skeletal muscle organoid development.