An in silico FSHD muscle fibre for modelling DUX4 dynamics and predicting the impact of therapy

  1. Matthew V Cowley
  2. Johanna Pruller
  3. Massimo Ganassi
  4. Peter S Zammit
  5. Christopher RS Banerji  Is a corresponding author
  1. University of Bath, United Kingdom
  2. King's College London, United Kingdom
  3. The Alan Turing Institute, United Kingdom

Abstract

Facioscapulohumeral muscular dystrophy (FSHD) is an incurable myopathy linked to over-expression of the myotoxic transcription factor DUX4. Targeting DUX4 is the leading therapeutic approach, however it is only detectable in 0.1-3.8% of FSHD myonuclei. How rare DUX4 drives FSHD and the optimal anti-DUX4 strategy is unclear. We combine stochastic gene expression with compartment models of cell states, building a simulation of DUX4 expression and consequences in FSHD muscle fibres. Investigating iDUX4 myoblasts, scRNAseq and snRNAseq of FSHD muscle we estimate parameters including DUX4 mRNA degradation, transcription and translation rates and DUX4 target gene activation rates. Our model accurately recreates the distribution of DUX4 and target gene positive cells seen in scRNAseq of FSHD myocytes. Importantly we show DUX4 drives significant cell death despite expression in only 0.8% of live cells. Comparing scRNAseq of unfused FSHD myocytes to snRNAseq of fused FSHD myonuclei, we find evidence of DUX4 protein syncytial diffusion and estimate its rate via genetic algorithms. We package our model into freely available tools, to rapidly investigate consequences of anti-DUX4 therapy.

Data availability

All data generated or analysed during this study are publicly available or included in the manuscript, all code employed is published as part of our shiny app at 3 public domain URLs listed in the manuscript.

The following previously published data sets were used

Article and author information

Author details

  1. Matthew V Cowley

    Department of Chemistry, University of Bath, Bath, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5258-8024
  2. Johanna Pruller

    Randall Centre for Cell and Molecular Biophysics, King's College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  3. Massimo Ganassi

    Randall Centre for Cell and Molecular Biophysics, King's College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3163-9707
  4. Peter S Zammit

    Randall Centre for Cell and Molecular Biophysics, King's College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9562-3072
  5. Christopher RS Banerji

    The Alan Turing Institute, London, United Kingdom
    For correspondence
    cbanerji@turing.ac.uk
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4373-7657

Funding

EPSRC Centre for Doctoral Training in Sustainable Chemical Technologies (EP/L016354/1)

  • Matthew V Cowley

Friends of FSH Research

  • Matthew V Cowley

Muscular Dystrophy UK (19GRO-PG12-0493)

  • Johanna Pruller

FSHD Society (FSHD-Winter2021-4491649104)

  • Johanna Pruller

Medical Research Council (MR/S002472/1)

  • Massimo Ganassi

Association Francaise contre les Myopathies

  • Peter S Zammit

SOLVE FSHD

  • Massimo Ganassi

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

Reviewing Editor

  1. Murim Choi, Seoul National University, Republic of Korea

Version history

  1. Preprint posted: December 12, 2022 (view preprint)
  2. Received: April 13, 2023
  3. Accepted: May 14, 2023
  4. Accepted Manuscript published: May 15, 2023 (version 1)
  5. Version of Record published: June 22, 2023 (version 2)

Copyright

© 2023, Cowley et al.

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

Metrics

  • 914
    views
  • 133
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Matthew V Cowley
  2. Johanna Pruller
  3. Massimo Ganassi
  4. Peter S Zammit
  5. Christopher RS Banerji
(2023)
An in silico FSHD muscle fibre for modelling DUX4 dynamics and predicting the impact of therapy
eLife 12:e88345.
https://doi.org/10.7554/eLife.88345

Share this article

https://doi.org/10.7554/eLife.88345

Further reading

    1. Cell Biology
    2. Neuroscience
    Marcos Moreno-Aguilera, Alba M Neher ... Carme Gallego
    Research Article Updated

    Alternative RNA splicing is an essential and dynamic process in neuronal differentiation and synapse maturation, and dysregulation of this process has been associated with neurodegenerative diseases. Recent studies have revealed the importance of RNA-binding proteins in the regulation of neuronal splicing programs. However, the molecular mechanisms involved in the control of these splicing regulators are still unclear. Here, we show that KIS, a kinase upregulated in the developmental brain, imposes a genome-wide alteration in exon usage during neuronal differentiation in mice. KIS contains a protein-recognition domain common to spliceosomal components and phosphorylates PTBP2, counteracting the role of this splicing factor in exon exclusion. At the molecular level, phosphorylation of unstructured domains within PTBP2 causes its dissociation from two co-regulators, Matrin3 and hnRNPM, and hinders the RNA-binding capability of the complex. Furthermore, KIS and PTBP2 display strong and opposing functional interactions in synaptic spine emergence and maturation. Taken together, our data uncover a post-translational control of splicing regulators that link transcriptional and alternative exon usage programs in neuronal development.

    1. Cell Biology
    Ang Li, Jianxun Yi ... Jingsong Zhou
    Research Article

    Amyotrophic lateral sclerosis (ALS) is a fatal neuromuscular disorder characterized by progressive weakness of almost all skeletal muscles, whereas extraocular muscles (EOMs) are comparatively spared. While hindlimb and diaphragm muscles of end-stage SOD1G93A (G93A) mice (a familial ALS mouse model) exhibit severe denervation and depletion of Pax7+satellite cells (SCs), we found that the pool of SCs and the integrity of neuromuscular junctions (NMJs) are maintained in EOMs. In cell sorting profiles, SCs derived from hindlimb and diaphragm muscles of G93A mice exhibit denervation-related activation, whereas SCs from EOMs of G93A mice display spontaneous (non-denervation-related) activation, similar to SCs from wild-type mice. Specifically, cultured EOM SCs contain more abundant transcripts of axon guidance molecules, including Cxcl12, along with more sustainable renewability than the diaphragm and hindlimb counterparts under differentiation pressure. In neuromuscular co-culture assays, AAV-delivery of Cxcl12 to G93A-hindlimb SC-derived myotubes enhances motor neuron axon extension and innervation, recapitulating the innervation capacity of EOM SC-derived myotubes. G93A mice fed with sodium butyrate (NaBu) supplementation exhibited less NMJ loss in hindlimb and diaphragm muscles. Additionally, SCs derived from G93A hindlimb and diaphragm muscles displayed elevated expression of Cxcl12 and improved renewability following NaBu treatment in vitro. Thus, the NaBu-induced transcriptomic changes resembling the patterns of EOM SCs may contribute to the beneficial effects observed in G93A mice. More broadly, the distinct transcriptomic profile of EOM SCs may offer novel therapeutic targets to slow progressive neuromuscular functional decay in ALS and provide possible ‘response biomarkers’ in pre-clinical and clinical studies.