Myogenesis in the Genomics Era

Comparative muscle transcriptome of Mali and Hampshire breeds of pigs: a preliminary study

Muscle development is an important priority of pig breeding programs. There is a considerable variation in muscularity between the breeds, but the regulation mechanisms of genes underlying myogenesis are still unclear. Transcriptome data from two breeds of pigs with divergent muscularity (Mali and Hampshire) were integrated with histology, immunofluorescence and meat yield to identify differences in myogenesis during the early growth phase. The muscle transcriptomics analysis revealed 17,721 common, 1413 and 1115 unique transcripts to Hampshire and Mali, respectively. This study identified 908 differentially expressed genes (p < 0.05; log2FC > ±1) in the muscle samples, of which 550 were upregulated and 358 were downregulated in Hampshire pigs, indicating differences in physiological process related to muscle function and development. Expression of genes related to myoblast fusion (MYMK), skeletal muscle satellite cell proliferation (ANGPT1, CDON) and growth factors (HGF, IGF1, IGF2) were higher in Hampshire than Mali, even though transcript levels of several other myogenesis-related genes (MYF6, MYOG, MSTN) were similar. The number of fibers per fascicle and the expression of myogenic marker proteins (MYOD1, MYOG and PAX7) were more in Hampshire as compared to Mali breed of pig, supporting results of transcriptome studies. The results suggest that differences in muscularity between breeds could be related to the regulation of myoblast fusion and myogenic activities. The present study will help to identify genes that could be explored for their utility in the selection of animals with different muscularities.

Single-cell RNA sequencing in skeletal muscle developmental biology

Skeletal muscle is the most extensive tissue in mammals, and they perform several functions; it is derived from paraxial mesodermal somites and undergoes hyperplasia and hypertrophy to form multinucleated, contractile, and functional muscle fibers. Skeletal muscle is a complex heterogeneous tissue composed of various cell types that establish communication strategies to exchange biological information; therefore, characterizing the cellular heterogeneity and transcriptional signatures of skeletal muscle is central to understanding its ontogeny's details. Studies of skeletal myogenesis have focused primarily on myogenic cells' proliferation, differentiation, migration, and fusion and ignored the intricate network of cells with specific biological functions. The rapid development of single-cell sequencing technology has recently enabled the exploration of skeletal muscle cell types and molecular events during development. This review summarizes the progress in single-cell RNA sequencing and its applications in skeletal myogenesis, which will provide insights into skeletal muscle pathophysiology.

Heterogeneous Skeletal Muscle Cell and Nucleus Populations Identified by Single-Cell and Single-Nucleus Resolution Transcriptome Assays

Single-cell RNA-seq (scRNA-seq) has revolutionized modern genomics, but the large size of myotubes and myofibers has restricted use of scRNA-seq in skeletal muscle. For the study of muscle, single-nucleus RNA-seq (snRNA-seq) has emerged not only as an alternative to scRNA-seq, but as a novel method providing valuable insights into multinucleated cells such as myofibers. Nuclei within myofibers specialize at junctions with other cell types such as motor neurons. Nuclear heterogeneity plays important roles in certain diseases such as muscular dystrophies. We survey current methods of high-throughput single cell and subcellular resolution transcriptomics, including single-cell and single-nucleus RNA-seq and spatial transcriptomics, applied to satellite cells, myoblasts, myotubes and myofibers. We summarize the major myonuclei subtypes identified in homeostatic and regenerating tissue including those specific to fiber type or at junctions with other cell types. Disease-specific nucleus populations were found in two muscular dystrophies, FSHD and Duchenne muscular dystrophy, demonstrating the importance of performing transcriptome studies at the single nucleus level in muscle.

Tributyrin, a Butyrate Pro-Drug, Primes Satellite Cells for Differentiation by Altering the Epigenetic Landscape

Satellite cells (SC) are a population of muscle resident stem cells that are responsible for postnatal muscle growth and repair. With investigation into the genomic regulation of SC fate, the role of the epigenome in governing SC myogenesis is becoming clearer. Histone deacetylase (HDAC) inhibitors have been demonstrated to be effective at enhancing the myogenic program of SC, but their role in altering the epigenetic landscape of SC remains undetermined. Our objective was to determine how an HDAC inhibitor, butyrate, promotes myogenic differentiation. SC from tributyrin treated neonatal piglets showed a decrease relative to SC from control animals in the expression of enhance of zeste homologue-2 (EZH2), a chromatin modifier, ex vivo. Chromatin Immunoprecipitation-Sequencing (ChIP-Seq) analysis of SC isolated from tributyrin treated pigs showed a global reduction of the tri-methylation of lysine 27 of histone H3 (H3K27me3) repressive chromatin mark. To determine if reductions in EZH2 was the primary mechanism through which butyrate affects SC behavior, SC were transfected with siRNA targeting EZH2, treated with 0.5 mM butyrate, or both. Treatment with butyrate reduced paired-box-7 (Pax7) and myogenic differentiation-1 (MyoD) gene expression, while siRNA caused reductions in EZH2 had no effect on their expression. EZH2 depletion did result in an increase in differentiating SC, but not in myotube hypertrophy. These results indicate that while EZH2 reduction may force myogenic differentiation, butyrate may operate through a parallel mechanism to enhance the myogenic program.

Incomplete MyoD-induced transdifferentiation is associated with chromatin remodeling deficiencies

Our current understanding of cellular transdifferentiation systems is limited. It is oftentimes unknown, at a genome-wide scale, how much transdifferentiated cells differ quantitatively from both the starting cells and the target cells. Focusing on transdifferentiation of primary human skin fibroblasts by forced expression of myogenic transcription factor MyoD, we performed quantitative analyses of gene expression and chromatin accessibility profiles of transdifferentiated cells compared to fibroblasts and myoblasts. In this system, we find that while many of the early muscle marker genes are reprogrammed, global gene expression and accessibility changes are still incomplete when compared to myoblasts. In addition, we find evidence of epigenetic memory in the transdifferentiated cells, with reminiscent features of fibroblasts being visible both in chromatin accessibility and gene expression. Quantitative analyses revealed a continuum of changes in chromatin accessibility induced by MyoD, and a strong correlation between chromatin-remodeling deficiencies and incomplete gene expression reprogramming. Classification analyses identified genetic and epigenetic features that distinguish reprogrammed from non-reprogrammed sites, and suggested ways to potentially improve transdifferentiation efficiency. Our approach for combining gene expression, DNA accessibility, and protein-DNA binding data to quantify and characterize the efficiency of cellular transdifferentiation on a genome-wide scale can be applied to any transdifferentiation system.

Temporal regulation of chromatin during myoblast differentiation

The commitment to and execution of differentiation programmes involves a significant change in gene expression in the precursor cell to facilitate development of the mature cell type. In addition to being regulated by lineage-determining and auxiliary transcription factors that drive these changes, the structural status of the chromatin has a considerable impact on the transcriptional competence of differentiation-specific genes, which is clearly demonstrated by the large number of cofactors and the extraordinary complex mechanisms by which these genes become activated. The terminal differentiation of myoblasts to myotubes and mature skeletal muscle is an excellent system to illustrate these points. The MyoD family of closely related, lineage-determining transcription factors directs, largely through targeting to chromatin, a cascade of cooperating transcription factors and enzymes that incorporate or remove variant histones, post-translationally modify histones, and alter nucleosome structure and positioning via energy released by ATP hydrolysis. The coordinated action of these transcription factors and enzymes prevents expression of differentiation-specific genes in myoblasts and facilitates the transition of these genes from transcriptionally repressed to activated during the differentiation process. Regulation is achieved in both a temporal as well as spatial manner, as at least some of these factors and enzymes affect local chromatin structure at myogenic gene regulatory sequences as well as higher-order genome organization. Here we discuss the transition of genes that promote myoblast differentiation from the silenced to the activated state with an emphasis on the changes that occur to individual histones and the chromatin structure present at these loci.

Nuclear bodies reorganize during myogenesis in vitro and are differentially disrupted by expression of FSHD-associated DUX4

Background

Nuclear bodies, such as nucleoli, PML bodies, and SC35 speckles, are dynamic sub-nuclear structures that regulate multiple genetic and epigenetic processes. Additional regulation is provided by RNA/DNA handling proteins, notably TDP-43 and FUS, which have been linked to ALS pathology. Previous work showed that mouse cell line myotubes have fewer but larger nucleoli than myoblasts, and we had found that nuclear aggregation of TDP-43 in human myotubes was induced by expression of DUX4-FL, a transcription factor that is aberrantly expressed and causes pathology in facioscapulohumeral dystrophy (FSHD). However, questions remained about nuclear bodies in human myogenesis and in muscle disease.

Methods

We examined nucleoli, PML bodies, SC35 speckles, TDP-43, and FUS in myoblasts and myotubes derived from healthy donors and from patients with FSHD, laminin-alpha-2-deficiency (MDC1A), and alpha-sarcoglycan-deficiency (LGMD2D). We further examined how these nuclear bodies and proteins were affected by DUX4-FL expression.

Results

We found that nucleoli, PML bodies, and SC35 speckles reorganized during differentiation in vitro, with all three becoming less abundant in myotube vs. myoblast nuclei. In addition, though PML bodies did not change in size, both nucleoli and SC35 speckles were larger in myotube than myoblast nuclei. Similar patterns of nuclear body reorganization occurred in healthy control, MDC1A, and LGMD2D cultures, as well as in the large fraction of nuclei that did not show DUX4-FL expression in FSHD cultures. In contrast, nuclei that expressed endogenous or exogenous DUX4-FL, though retaining normal nucleoli, showed disrupted morphology of some PML bodies and most SC35 speckles and also co-aggregation of FUS with TDP-43.

Conclusions

Nucleoli, PML bodies, and SC35 speckles reorganize during human myotube formation in vitro. These nuclear body reorganizations are likely needed to carry out the distinct gene transcription and splicing patterns that are induced upon myotube formation. DUX4-FL-induced disruption of some PML bodies and most SC35 speckles, along with co-aggregation of TDP-43 and FUS, could contribute to pathogenesis in FSHD, perhaps by locally interfering with genetic and epigenetic regulation of gene expression in the small subset of nuclei that express high levels of DUX4-FL at any one time.

Genome-wide association between Six4, MyoD, and the histone demethylase Utx during myogenesis

Adult skeletal muscles can regenerate after injury, due to the presence of satellite cells, a quiescent population of myogenic progenitor cells. Once activated, satellite cells repair the muscle damage by undergoing myogenic differentiation. The myogenic regulatory factors (MRFs) coordinate the process of progenitor differentiation in cooperation with other families of transcription factors (TFs). The Six1 and Six4 homeodomain TFs are expressed in developing and adult muscle and Six1 is critical for embryonic and adult myogenesis. However, the lack of a muscle developmental phenotype in Six4-null mice, which has been attributed to compensation by other Six family members, has discouraged further assessment of the role of Six4 during adult muscle regeneration. By employing genome-wide approaches to address the function of Six4 during adult skeletal myogenesis, we have identified a core set of muscle genes coordinately regulated in adult muscle precursors by Six4 and the MRF MyoD. Throughout the genome of differentiating adult myoblasts, the cooperation between Six4 and MyoD is associated with chromatin repressive mark removal by Utx, a demethylase of histone H3 trimethylated at lysine 27. Among the genes coordinately regulated by Six4 and MyoD are several genes critical for proper in vivo muscle regeneration, implicating a role of Six4 in this process. Using in vivo RNA interference of Six4, we expose an uncompensated function of this TF during muscle regeneration. Together, our results reveal a role for Six4 during adult muscle regeneration and suggest a widespread mechanism of cooperation between Six4 and MyoD.