Expression of the human FSHD-linked DUX4 gene induces neurogenesis during differentiation of murine embryonic stem cells

The FSHD muscle-blood biomarker: a circulating transcriptomic biomarker for clinical severity in facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable skeletal myopathy. Clinical trials for FSHD are hindered by heterogeneous biomarkers poorly associated with clinical severity, requiring invasive muscle biopsy. Macroscopically, FSHD presents with slow fatty replacement of muscle, rapidly accelerated by inflammation. Mis-expression of the transcription factor DUX4 is currently accepted to underlie pathogenesis, and mechanisms including PAX7 target gene repression have been proposed. Here, we performed RNA-sequencing on MRI-guided inflamed and isogenic non-inflamed muscle biopsies from the same clinically characterized FSHD patients (n = 24), alongside isogenic peripheral blood mononucleated cells from a subset of patients (n = 13) and unaffected controls (n = 11). Multivariate models were employed to evaluate the clinical associations of five published FSHD transcriptomic biomarkers. We demonstrated that PAX7 target gene repression can discriminate control, inflamed and non-inflamed FSHD muscle independently of age and sex (P < 0.013), while the discriminatory power of DUX4 target genes was limited to distinguishing FSHD muscle from control. Importantly, the level of PAX7 target gene repression in non-inflamed muscle associated with clinical assessments of FSHD severity (P = 0.04). DUX4 target gene biomarkers in FSHD muscle showed associations with lower limb fat fraction and D4Z4 array length but not clinical assessment. Lastly, PAX7 target gene repression in FSHD muscle correlated with the level in isogenic peripheral blood mononucleated cells (P = 0.002). A refined PAX7 target gene biomarker comprising 143/601 PAX7 target genes computed in peripheral blood (the FSHD muscle-blood biomarker) associated with clinical severity in FSHD patients (P < 0.036). Our new circulating biomarker validates as a classifier of clinical severity in an independent data set of 54 FSHD patient and 29 matched control blood samples, with improved power in older patients (P = 0.03). In summary, we present the minimally invasive FSHD muscle-blood biomarker of FSHD clinical severity valid in patient muscle and blood, of potential use in routine disease monitoring and clinical trials.

Canine DUXC: Implications for DUX4 retrotransposition and preclinical models of FSHD

Mis-expression of DUX4 in skeletal muscle causes facioscapulohumeral muscular dystrophy (FSHD). Human DUX4 and mouse Dux are retrogenes derived from retrotransposition of the mRNA from the parental DUXC gene. Primates and rodents have lost the parental DUXC gene and it is unknown whether DUXC had a similar role in driving an early pluripotent transcriptional program. Dogs and other Laurasiatherians have retained DUXC, providing an opportunity to determine the functional similarity to the retrotransposed DUX4 and Dux. Here we identify expression of two isoforms of DUXC mRNA in canine testis tissues, one encoding the canonical double homeodomain protein (DUXC), similar to DUX4/Dux, and a second that includes an in-frame alternative exon that disrupts the conserved amino-acid sequence of the first homeodomain (DUXC-ALT). Expression of DUXC in canine cells induces a pluripotent program similar to DUX4 and Dux, and induces expression of a similar set of retrotransposons of the ERV/MaLR and LINE-1 families, as well as pericentromeric satellite repeats; whereas DUXC-ALT did not robustly activate gene expression in these assays. Important for preclinical models of FSHD, human DUX4 and canine DUXC show higher conservation of their homeodomains and corresponding binding motifs compared to the conservation between human DUX4 and mouse Dux, and human DUX4 activates a highly similar transcriptional program in canine cells. Together, these findings show that retrotransposition resulted in the loss of an alternatively spliced isoform and that DUXC containing mammals might be good candidates for certain pre-clinical models of FSHD.

A SARS-CoV-2 host infection model network based on genomic human Transcription Factors (TFs) depletion

In December 2019 a new beta-coronavirus was isolated and characterized by sequencing samples from pneumonia patients in Wuhan, Hubei Province, China. Coronaviruses are positive-sense RNA viruses widely distributed among different animal species and humans in which they cause respiratory, enteric, liver and neurological symptomatology. Six species of coronavirus have been described (HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1) that cause cold-like symptoms in immunocompetent or immunocompromised subjects and two strains of sometimes fatal zoonotic origin that cause severe acute respiratory syndrome (SARS-CoV and MERS-CoV). The SARS-CoV-2 strain is the emerging seventh member of the coronavirus family, which is actually determining a global emergency.

In silico analysis is a promising approach for understanding biological events in complex diseases and due to serious worldwide emergency and serious threat to global health, it is extremely important to use bioinformatics methods able to study an emerging pathogen like SARS-CoV-2.

Herein, we report on in silico comparative analysis between complete genome of SARS-CoV, MERS-CoV, HCoV-OC43 and SARS-CoV-2 strains, to identify the occurrence of specific conserved motifs on viral genomic sequences which should be able to bind and therefore induce a subtraction of host's Transcription Factors (TFs) which lead to a depletion, an effect comparable to haploinsufficiency (a genetic dominant condition in which a single copy of wild-type allele at a locus, in heterozygous combination with a variant allele, is insufficient to produce the correct quantity of transcript and, therefore, of protein, for a correct standard phenotypic expression).

In this competitive scenario, virus versus host, the proposed in silico protocol identified the TFs same as the distribution of TFBSs (Transcription Factor Binding Sites) on analyzed viral strains, potentially able to influence genes and pathways with biological functions confirming that this approach could brings useful insights regarding SARS-CoV-2. According to our results obtained by this in silico approach it is possible to hypothesize that TF-binding motifs could be of help in the explanation of the complex and heterogeneous clinical presentation in SARS-CoV-2 and subsequently predict possible interactions regarding metabolic pathways, and drug or target relationships.

Mouse Dux is myotoxic and shares partial functional homology with its human paralog DUX4

D4Z4 repeats are present in at least 11 different mammalian species, including humans and mice. Each repeat contains an open reading frame encoding a double homeodomain (DUX) family transcription factor. Aberrant expression of the D4Z4 ORF called DUX4 is associated with the pathogenesis of Facioscapulohumeral muscular dystrophy (FSHD). DUX4 is toxic to numerous cell types of different species, and over-expression caused dysmorphism and developmental arrest in frogs and zebrafish, embryonic lethality in transgenic mice, and lesions in mouse muscle. Because DUX4 is a primate-specific gene, questions have been raised about the biological relevance of over-expressing it in non-primate models, as DUX4 toxicity could be related to non-specific cellular stress induced by over-expressing a DUX family transcription factor in organisms that did not co-evolve its regulated transcriptional networks. We assessed toxic phenotypes of DUX family genes, including DUX4, DUX1, DUX5, DUXA, DUX4-s, Dux-bl and mouse Dux. We found that DUX proteins were not universally toxic, and only the mouse Dux gene caused similar toxic phenotypes as human DUX4. Using RNA-seq, we found that 80% of genes upregulated by Dux were similarly increased in DUX4-expressing cells. Moreover, 43% of Dux-responsive genes contained ChIP-seq binding sites for both Dux and DUX4, and both proteins had similar consensus binding site sequences. These results suggested DUX4 and Dux may regulate some common pathways, and despite diverging from a common progenitor under different selective pressures for millions of years, the two genes maintain partial functional homology.

Dux4 controls migration of mesenchymal stem cells through the Cxcr4-Sdf1 axis

We performed transcriptome profiling of human immortalized myoblasts (MB) transiently expressing double homeobox transcription factor 4 (DUX4) and double homeobox transcription factor 4 centromeric (DUX4c) and identified 114 and 70 genes differentially expressed in DUX4- and DUX4c-transfected myoblasts, respectively. A significant number of differentially expressed genes were involved in inflammation, cellular migration and chemotaxis suggesting a role for DUX4 and DUX4c in these processes. DUX4 but not DUX4c overexpression resulted in upregulation of the CXCR4 (C-X-C motif Receptor 4) and CXCL12 (C-X-C motif ligand 12 also known as SDF1) expression in human immortalized myoblasts. In a Transwell cell migration assay, human bone marrow-derived mesenchymal stem cells (BMSCs) were migrating more efficiently towards human immortalized myoblasts overexpressing DUX4 as compared to controls; the migration efficiency of DUX4-transfected BMSCs was also increased. DUX4c overexpression in myoblasts or in BMSCs had no impact on the rate of BMSC migration. Antibodies against SDF1 and CXCR4 blocked the positive effect of DUX4 overexpression on BMSC migration. We propose that DUX4 controls the cellular migration of mesenchymal stem cells through the CXCR4 receptor.

Ret function in muscle stem cells points to tyrosine kinase inhibitor therapy for facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) involves sporadic expression of DUX4, which inhibits myogenesis and is pro-apoptotic. To identify target genes, we over-expressed DUX4 in myoblasts and found that the receptor tyrosine kinase Ret was significantly up-regulated, suggesting a role in FSHD. RET is dynamically expressed during myogenic progression in mouse and human myoblasts. Constitutive expression of either RET9 or RET51 increased myoblast proliferation, whereas siRNA-mediated knockdown of Ret induced myogenic differentiation. Suppressing RET activity using Sunitinib, a clinically-approved tyrosine kinase inhibitor, rescued differentiation in both DUX4-expressing murine myoblasts and in FSHD patient-derived myoblasts. Importantly, Sunitinib also increased engraftment and differentiation of FSHD myoblasts in regenerating mouse muscle. Thus, DUX4-mediated activation of Ret prevents myogenic differentiation and could contribute to FSHD pathology by preventing satellite cell-mediated repair. Rescue of DUX4-induced pathology by Sunitinib highlights the therapeutic potential of tyrosine kinase inhibitors for treatment of FSHD.

DOI:http://dx.doi.org/10.7554/eLife.11405.001

DNA-binding sequence specificity of DUX4

Background

Misexpression of the double homeodomain transcription factor DUX4 results in facioscapulohumeral muscular dystrophy (FSHD). A DNA-binding consensus with two tandem TAAT motifs based on chromatin IP peaks has been discovered; however, the consensus has multiple variations (flavors) of unknown relative activity. In addition, not all peaks have this consensus, and the Pitx1 promoter, the first DUX4 target sequence mooted, has a different TAAT-rich sequence. Furthermore, it is not known whether and to what extent deviations from the consensus affect DNA-binding affinity and transcriptional activation potential.

Results

Here, we take both unbiased and consensus sequence-driven approaches to determine the DNA-binding specificity of DUX4 and its tolerance to mismatches at each site within its consensus sequence. We discover that the best binding and the greatest transcriptional activation are observed when the two TAAT motifs are separated by a C residue. The second TAAT motif in the consensus sequence is actually (T/C)AAT. We find that a T is preferred here. DUX4 has no transcriptional activity on “half-sites”, i.e., those bearing only a single TAAT motif. We further find that DUX4 does not bind to the TAATTA motif in the Pitx1 promoter, that Pitx1 sequences have no competitive band shift activity, and that the Pitx1 sequence is transcriptionally inactive, calling into question PITX1 as a DUX4 target gene. Finally, by multimerizing binding sites, we find that DUX4 transcriptional activation demonstrates tremendous synergy and that at low DNA concentrations, at least two motifs are necessary to detect a transcriptional response.

Conclusions

These studies illuminate the DNA-binding sequence preferences of DUX4.

Electronic supplementary material

The online version of this article (doi:10.1186/s13395-016-0080-z) contains supplementary material, which is available to authorized users.

Homologous Transcription Factors DUX4 and DUX4c Associate with Cytoplasmic Proteins during Muscle Differentiation

Hundreds of double homeobox (DUX) genes map within 3.3-kb repeated elements dispersed in the human genome and encode DNA-binding proteins. Among these, we identified DUX4, a potent transcription factor that causes facioscapulohumeral muscular dystrophy (FSHD). In the present study, we performed yeast two-hybrid screens and protein co-purifications with HaloTag-DUX fusions or GST-DUX4 pull-down to identify protein partners of DUX4, DUX4c (which is identical to DUX4 except for the end of the carboxyl terminal domain) and DUX1 (which is limited to the double homeodomain). Unexpectedly, we identified and validated (by co-immunoprecipitation, GST pull-down, co-immunofluorescence and in situ Proximal Ligation Assay) the interaction of DUX4, DUX4c and DUX1 with type III intermediate filament protein desmin in the cytoplasm and at the nuclear periphery. Desmin filaments link adjacent sarcomere at the Z-discs, connect them to sarcolemma proteins and interact with mitochondria. These intermediate filament also contact the nuclear lamina and contribute to positioning of the nuclei. Another Z-disc protein, LMCD1 that contains a LIM domain was also validated as a DUX4 partner. The functionality of DUX4 or DUX4c interactions with cytoplasmic proteins is underscored by the cytoplasmic detection of DUX4/DUX4c upon myoblast fusion. In addition, we identified and validated (by co-immunoprecipitation, co-immunofluorescence and in situ Proximal Ligation Assay) as DUX4/4c partners several RNA-binding proteins such as C1QBP, SRSF9, RBM3, FUS/TLS and SFPQ that are involved in mRNA splicing and translation. FUS and SFPQ are nuclear proteins, however their cytoplasmic translocation was reported in neuronal cells where they associated with ribonucleoparticles (RNPs). Several other validated or identified DUX4/DUX4c partners are also contained in mRNP granules, and the co-localizations with cytoplasmic DAPI-positive spots is in keeping with such an association. Large muscle RNPs were recently shown to exit the nucleus via a novel mechanism of nuclear envelope budding. Following DUX4 or DUX4c overexpression in muscle cell cultures, we observed their association with similar nuclear buds. In conclusion, our study demonstrated unexpected interactions of DUX4/4c with cytoplasmic proteins playing major roles during muscle differentiation. Further investigations are on-going to evaluate whether these interactions play roles during muscle regeneration as previously suggested for DUX4c.

The Role of D4Z4-Encoded Proteins in the Osteogenic Differentiation of Mesenchymal Stromal Cells Isolated from Bone Marrow

Facioscapulohumeral muscular dystrophy (FSHD) is associated with an activation of the double homeobox 4 (DUX4) gene, which we previously identified within the D4Z4 repeated elements in the 4q35 subtelomeric region. The pathological DUX4 mRNA is derived from the most distal D4Z4 unit and extends unexpectedly within the flanking pLAM region, which provides an intron and polyadenylation signal. The conditions that are required to develop FSHD are a permissive allele providing the polyadenylation signal and hypomethylation of the D4Z4 repeat array compared with the healthy muscle. The DUX4 protein is a 52-kDa transcription factor that initiates a large gene deregulation cascade leading to muscle atrophy, inflammation, differentiation defects, and oxidative stress, which are the key features of FSHD. DUX4 is a retrogene that is normally expressed in germline cells and is submitted to repeat-induced silencing in adult tissues. Since DUX4 mRNAs have been detected in human embryonic and induced pluripotent stem cells, we investigated whether they could also be expressed in human mesenchymal stromal cells (hMSCs). We found that DUX4 mRNAs were induced during the differentiation of hMSCs into osteoblasts and that this process involved DUX4 and new longer protein forms (58 and 70 kDa). A DUX4 mRNA with a more distant 5' start site was characterized that presented a 60-codon reading frame extension and encoded the 58-kDa protein. Transfections of hMSCs with an antisense oligonucleotide targeting DUX4 mRNAs decreased both the 52- and 58-kDa protein levels and confirmed their identity. Gain- and loss-of-function experiments in hMSCs suggested these DUX4 proteins had opposite roles in osteogenic differentiation as evidenced by the alkaline phosphatase activity and calcium deposition. Differentiation was delayed by the 58-kDa DUX4 expression and it was increased by 52-kDa DUX4. These data indicate a role for DUX4 protein forms in the osteogenic differentiation of hMSCs.

Facioscapulohumeral dystrophy: the path to consensus on pathophysiology

Although the pathophysiology of facioscapulohumeral dystrophy (FSHD) has been controversial over the last decades, progress in recent years has led to a model that incorporates these decades of findings and is gaining general acceptance in the FSHD research community. Here we review how the contributions from many labs over many years led to an understanding of a fundamentally new mechanism of human disease. FSHD is caused by inefficient repeat-mediated epigenetic repression of the D4Z4 macrosatellite repeat array on chromosome 4, resulting in the variegated expression of the DUX4 retrogene, encoding a double-homeobox transcription factor, in skeletal muscle. Normally expressed in the testis and epigenetically repressed in somatic tissues, DUX4 expression in skeletal muscle induces expression of many germline, stem cell, and other genes that might account for the pathophysiology of FSHD. Although some disagreements regarding the details of mechanisms remain in the field, the coalescing agreement on a central model of pathophysiology represents a pivot-point in FSHD research, transitioning the field from discovery-oriented studies to translational studies aimed at developing therapies based on a sound model of disease pathophysiology.