Single-cell spatial transcriptomics reveals a dystrophic trajectory following a developmental bifurcation of myoblast cell fates in facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) is linked to abnormal derepression of the transcription activator DUX4. This effect is localized to a low percentage of cells, requiring single-cell analysis. However, single-cell/nucleus RNA-seq cannot fully capture the transcriptome of multinucleated large myotubes. To circumvent these issues, we use multiplexed error-robust fluorescent in situ hybridization (MERFISH) spatial transcriptomics that allows profiling of RNA transcripts at a subcellular resolution. We simultaneously examined spatial distributions of 140 genes, including 24 direct DUX4 targets, in in vitro differentiated myotubes and unfused mononuclear cells (MNCs) of control, isogenic D4Z4 contraction mutant and FSHD patient samples, as well as the individual nuclei within them. We find myocyte nuclei segregate into two clusters defined by the expression of DUX4 target genes, which is exclusively found in patient/mutant nuclei, whereas MNCs cluster based on developmental states. Patient/mutant myotubes are found in "FSHD-hi" and "FSHD-lo" states with the former signified by high DUX4 target expression and decreased muscle gene expression. Pseudotime analyses reveal a clear bifurcation of myoblast differentiation into control and FSHD-hi myotube branches, with variable numbers of DUX4 target-expressing nuclei found in multinucleated FSHD-hi myotubes. Gene coexpression modules related to extracellular matrix and stress gene ontologies are significantly altered in patient/mutant myotubes compared with the control. We also identify distinct subpathways within the DUX4 gene network that may differentially contribute to the disease transcriptomic phenotype. Taken together, our MERFISH-based study provides effective gene network profiling of multinucleated cells and identifies FSHD-induced transcriptomic alterations during myoblast differentiation.

Engineered FSHD mutations results in D4Z4 heterochromatin disruption and feedforward DUX4 network activation.

Facioscapulohumeral dystrophy (FSHD) is linked to contraction of D4Z4 repeats on chromosome 4q with SMCHD1 mutations acting as a disease modifier. D4Z4 heterochromatin disruption and abnormal upregulation of the transcription factor DUX4, encoded in the D4Z4 repeat, are the hallmarks of FSHD. However, defining the precise effect of D4Z4 contraction has been difficult because D4Z4 repeats are primate-specific and DUX4 expression is very rare in highly heterogeneous patient myocytes. We generated isogenic mutant cell lines harboring D4Z4 and/or SMCHD1 mutations in a healthy human skeletal myoblast line. We found that the mutations affect D4Z4 heterochromatin differently, and that SMCHD1 mutation or disruption of DNA methylation stabilizes otherwise variegated DUX4 target activation in D4Z4 contraction mutant cells, demonstrating the critical role of modifiers. Our study revealed amplification of the DUX4 signal through downstream targets, H3.X/Y and LEUTX. Our results provide important insights into how rare DUX4 expression leads to FSHD pathogenesis.