Differentiation is a key cellular process in normal tissue development that is significantly altered in cancer. Although molecular signatures characterising pluripotency and multipotency exist, there is, as yet, no single quantitative mark of a cellular sample's position in the global differentiation hierarchy. Here we adopt a systems view and consider the sample's network entropy, a measure of signaling pathway promiscuity, computable from a sample's genome-wide expression profile. We demonstrate that network entropy provides a quantitative, in-silico, readout of the average undifferentiated state of the profiled cells, recapitulating the known hierarchy of pluripotent, multipotent and differentiated cell types. Network entropy further exhibits dynamic changes in time course differentiation data, and in line with a sample's differentiation stage. In disease, network entropy predicts a higher level of cellular plasticity in cancer stem cell populations compared to ordinary cancer cells. Importantly, network entropy also allows identification of key differentiation pathways. Our results are consistent with the view that pluripotency is a statistical property defined at the cellular population level, correlating with intra-sample heterogeneity, and driven by the degree of signaling promiscuity in cells. In summary, network entropy provides a quantitative measure of a cell's undifferentiated state, defining its elevation in Waddington's landscape.

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable myopathy, linked to hypomethylation of D4Z4 repeats on chromosome 4q causing expression of the DUX4 transcription factor. However, DUX4 is difficult to detect in FSHD muscle biopsies and it is debatable how robust changes in DUX4 target gene expression are as an FSHD biomarker. PAX7 is a master regulator of myogenesis that rescues DUX4-mediated apoptosis. Here, we show that suppression of PAX7 target genes is a hallmark of FSHD, and that it is as major a signature of FSHD muscle as DUX4 target gene expression. This is shown using meta-analysis of over six FSHD muscle biopsy gene expression studies, and validated by RNA-sequencing on FSHD patient-derived myoblasts. DUX4 also inhibits PAX7 from activating its transcriptional target genes and vice versa. Furthermore, PAX7 target gene repression can explain oxidative stress sensitivity and epigenetic changes in FSHD. Thus, PAX7 target gene repression is a hallmark of FSHD that should be considered in the investigation of FSHD pathology and therapy.

Skeletal muscle wasting in facioscapulohumeral muscular dystrophy (FSHD) results in substantial morbidity. On a disease-permissive chromosome 4qA haplotype, genomic and/or epigenetic changes at the D4Z4 macrosatellite repeat allows transcription of the DUX4 retrogene. Analysing transgenic mice carrying a human D4Z4 genomic locus from an FSHD-affected individual showed that DUX4 was transiently induced in myoblasts during skeletal muscle regeneration. Centromeric to the D4Z4 repeats is an inverted D4Z4 unit encoding DUX4c. Expression of DUX4, DUX4c and DUX4 constructs, including constitutively active, dominant-negative and truncated versions, revealed that DUX4 activates target genes to inhibit proliferation and differentiation of satellite cells, but that it also downregulates target genes to suppress myogenic differentiation. These transcriptional changes elicited by DUX4 in mouse have significant overlap with genes regulated by DUX4 in man. Comparison of DUX4 and DUX4c transcriptional perturbations revealed that DUX4 regulates genes involved in cell proliferation, whereas DUX4c regulates genes engaged in angiogenesis and muscle development, with both DUX4 and DUX4c modifing genes involved in urogenital development. Transcriptomic analysis showed that DUX4 operates through both target gene activation and repression to orchestrate a transcriptome characteristic of a less-differentiated cell state.

Facioscapulohumeral muscular dystrophy (FSHD) is an incurable disease, characterized by skeletal muscle weakness and wasting. Genetically, FSHD is characterized by contraction or hypomethylation of repeat D4Z4 units on chromosome 4, which causes aberrant expression of the transcription factor DUX4 from the last repeat. Many genes have been implicated in FSHD pathophysiology, but an integrated molecular model is currently lacking. We developed a novel differential network methodology, Interactome Sparsification and Rewiring (InSpiRe), which detects network rewiring between phenotypes by integrating gene expression data with known protein interactions. Using InSpiRe, we performed a meta-analysis of multiple microarray datasets from FSHD muscle biopsies, then removed secondary rewiring using non-FSHD datasets, to construct a unified network of rewired interactions. Our analysis identified b-catenin as the main coordinator of FSHD-associated protein interaction signalling, with pathways including canonical Wnt, HIF1-a and TNF-a clearly perturbed. To detect transcriptional changes directly elicited by DUX4, gene expression profiling was performed using microarrays on murine myoblasts. This revealed that DUX4 significantly modified expression of the genes in our FSHD network. Furthermore, we experimentally confirmed that Wnt/ b-catenin signalling is affected by DUX4 in murine myoblasts. Thus, we provide the first unified molecular map of FSHD signalling, capable of uncovering pathomechanisms and guiding therapeutic development.

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable myopathy, linked to epigenetic derepression of D4Z4 repeats on chromosome 4q, leading to ectopic DUX4 expression. FSHD patient myoblasts have defective myogenic differentiation, forming smaller myotubes with reduced myosin content. However, molecular mechanisms driving such disrupted myogenesis in FSHD are poorly understood. We performed high-throughput morphological analysis describing FSHD and control myogenesis, revealing altered myogenic differentiation results in hypotrophic myotubes. Employing polynomial models and an empirical Bayes approach, we established eight critical time points during which human healthy and FSHD myogenesis differ. RNA-sequencing at these eight nodal time points in triplicate, provided temporal depth for a multivariate regression analysis, allowing assessment of interaction between progression of differentiation and FSHD disease status. Importantly, the unique size and structure of our data permitted identification of many novel FSHD pathomechanisms undetectable by previous approaches. For further analysis here, we selected pathways that control mitochondria: of interest considering known alterations in mitochondrial structure and function in FSHD muscle, and sensitivity of FSHD cells to oxidative stress. Notably, we identified suppression of mitochondrial biogenesis, in particular via peroxisome proliferator-activated receptor gamma coactivator 1-α (PGC1α), the cofactor and activator of oestrogen-related receptor α (ERRα). PGC1α knock-down caused hypotrophic myotubes to form from control myoblasts. Known ERRα agonists and safe food supplements biochanin A, daidzein or genistein, each rescued the hypotrophic FSHD myotube phenotype. Together our work describes transcriptomic changes in high resolution that occur during myogenesis in FSHD ex vivo, identifying suppression of the PGC1α-ERRα axis leading to perturbed myogenic differentiation, which can effectively be rescued by readily available food supplements.

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable skeletal myopathy. The condition is linked to hypomethylation of the D4Z4 macrosatellite repeat at chromosome 4q35, leading to epigenetic derepression of the transcription factor DUX4; coupled with a permissive 4qA haplotype supplying a poly(A) signal. DUX4 may drive FSHD pathology via both induction of target genes and inhibition of the function of the myogenic master regulator PAX7. Biomarkers for FSHD have focused on DUX4 target gene expression. We have, however, reported that PAX7 target gene repression is a hallmark of FSHD skeletal muscle. Here we demonstrate that PAX7 target gene repression is an equivalent biomarker to DUX4 target gene expression when considering RNA-Sequencing data from magnetic resonance imaging-guided muscle biopsies. Moreover, PAX7 target gene repression correlates with active disease, independent to DUX4 target gene expression. PAX7 target genes are also repressed in RNA-Sequencing data from single cells, representing a significantly better biomarker of FSHD cells than DUX4 target gene expression. Importantly, PAX7 target gene repression is a significant biomarker in the majority of FSHD cells that are DUX4 target gene negative, and on which the DUX4 biomarker is indiscriminate. To facilitate the evaluation of validated biomarkers we provide a simple tool that outputs biomarker values from a normalized expression data matrix. In summary, PAX7 target gene repression in FSHD correlates with disease severity, independently of DUX4 target gene expression. At the single-cell level, PAX7 target gene repression can efficiently discriminate FSHD cells, even when no DUX4 target genes are detectable.

One of the key characteristics of cancer cells is an increased phenotypic plasticity, driven by underlying genetic and epigenetic perturbations. However, at a systems-level it is unclear how these perturbations give rise to the observed increased plasticity. Elucidating such systems-level principles is key for an improved understanding of cancer. Recently, it has been shown that signaling entropy, an overall measure of signaling pathway promiscuity, and computable from integrating a sample's gene expression profile with a protein interaction network, correlates with phenotypic plasticity and is increased in cancer compared to normal tissue. Here we develop a computational framework for studying the effects of network perturbations on signaling entropy. We demonstrate that the increased signaling entropy of cancer is driven by two factors: (i) the scale-free (or near scale-free) topology of the interaction network, and (ii) a subtle positive correlation between differential gene expression and node connectivity. Indeed, we show that if protein interaction networks were random graphs, described by Poisson degree distributions, that cancer would generally not exhibit an increased signaling entropy. In summary, this work exposes a deep connection between cancer, signaling entropy and interaction network topology.

Facioscapulohumeral muscular dystrophy (FSHD) is characterised by progressive skeletal muscle weakness and wasting. FSHD is linked to epigenetic derepression of the subtelomeric D4Z4 macrosatellite at chromosome 4q35. Epigenetic derepression permits the distal‐most D4Z4 unit to transcribe DUX4, with transcripts stabilised by splicing to a poly(A) signal on permissive 4qA haplotypes. The pioneer transcription factor DUX4 activates target genes that are proposed to drive FSHD pathology. While this toxic gain‐of‐function model is a satisfying “bottom‐up” genotype‐to‐phenotype link, DUX4 is rarely detectable in muscle and DUX4 target gene expression is inconsistent in patients. A reliable biomarker for FSHD is suppression of a target gene score of PAX7, a master regulator of myogenesis. However, it is unclear how this “top‐down” finding links to genomic changes that characterise FSHD and to DUX4. Here, we explore the roles and interactions of DUX4 and PAX7 in FSHD pathology and how the relationship between these two transcription factors deepens understanding via the immune system and muscle regeneration. Considering how FSHD pathomechanisms are represented by “DUX4opathy” models has implications for developing therapies and current clinical trials.