Genome-wide MyoD binding in skeletal muscle cells: a potential for broad cellular reprogramming

Molecular and cellular regulation of skeletal myogenesis

Since the seminal discovery of the cell-fate regulator Myod, studies in skeletal myogenesis have inspired the search for cell-fate regulators of similar potential in other tissues and organs. It was perplexing that a similar transcription factor for other tissues was not found; however, it was later discovered that combinations of molecular regulators can divert somatic cell fates to other cell types. With the new era of reprogramming to induce pluripotent cells, the myogenesis paradigm can now be viewed under a different light. Here, we provide a short historical perspective and focus on how the regulation of skeletal myogenesis occurs distinctly in different scenarios and anatomical locations. In addition, some interesting features of this tissue underscore the importance of reconsidering the simple-minded view that a single stem cell population emerges after gastrulation to assure tissuegenesis. Notably, a self-renewing long-term Pax7+ myogenic stem cell population emerges during development only after a first wave of terminal differentiation occurs to establish a tissue anlagen in the mouse. How the future stem cell population is selected in this unusual scenario will be discussed. Recently, a wealth of information has emerged from epigenetic and genome-wide studies in myogenic cells. Although key transcription factors such as Pax3, Pax7, and Myod regulate only a small subset of genes, in some cases their genomic distribution and binding are considerably more promiscuous. This apparent nonspecificity can be reconciled in part by the permissivity of the cell for myogenic commitment, and also by new roles for some of these regulators as pioneer transcription factors acting on chromatin state.

'From death, lead me to immortality' - mantra of ageing skeletal muscle

Skeletal muscle is a post-mitotic tissue maintained by repair and regeneration through a population of stem cell-like satellite cells. Following muscle injury, satellite cell proliferation is mediated by local signals ensuring sufficient progeny for tissue repair. Age–related changes in satellite cells as well as to the local and systemic environment potentially impact on the capacity of satellite cells to generate sufficient progeny in an ageing organism resulting in diminished regeneration. ‘Rejuvenation’ of satellite cell progeny and regenerative capacity by environmental stimuli effectors suggest that a subset of age-dependent satellite cell changes may be reversible. Epigenetic regulation of satellite stem cells that include DNA methylation and histone modifications which regulate gene expression are potential mechanisms for such reversible changes and have been shown to control organismal longevity. The area of health and ageing that is likely to benefit soonest from advances in the biology of adult stem cells is the emerging field of regenerative medicine. Further studies are needed to elucidate the mechanisms by which epigenetic modifications regulate satellite stem cell function and will require an increased understanding of stem-cell biology, the environment of the aged tissue and the interaction between the two.

A 5'-regulatory region and two coding region polymorphisms modulate promoter activity and gene expression of the growth suppressor gene ZBED6 in cattle

Zinc finger, BED-type containing 6 (ZBED6) is an important transcription factor in placental mammals, affecting development, cell proliferation and growth. Polymorphisms in its promoter and coding regions are likely to impact ZBED6 transcription and growth traits. In this study, rapid amplification of 5’ cDNA ends (5'-RACE) analysis revealed two transcription start sites (TSS) for the bovine ZBED6 starting within exon 1 of the ZC3H11A gene (TSS-1) and upstream of the translation start codon of the ZBED6 gene (TSS-2). There was one SNP in the promoter and two missense mutations in the coding region of the bovine ZBED6 by sequencing of the pooled DNA samples (Pool-Seq, n = 100). The promoter and coding region are the key regions for gene function; polymorphisms in these regions can alter gene expression. Quantitative real-time PCR (qPCR) analysis showed that ZBED6 has a broad tissue distribution in cattle and is highly expressed in skeletal muscle. Eleven promoter-detection vectors were constructed, which enabled the cloning of putative promoter sequences and analysis of ZBED6 transcriptional activity by luciferase reporter gene assays. The core region of the basal promoter of bovine ZBED6 is located within region -866 to -556. The activity of WT-826G-pGL3 in driving reporter gene transcription is significantly higher than that of the M-826A-pGL3 construct (P < 0.01). Analysis of gene expression patterns in homozygous full-sibling Chinese Qinchuan cattle showed that the mutant-type Hap-AGG exhibited a lower mRNA level than the wild-type Hap-GCA (P < 0.05) in longissimus dorsi muscle (LDM). Moreover, ZBED6 mRNA expression was low in C2C12 cells overexpressing the mutant-type ZBED6 (pcDNA3.1⁺-Hap-GG) (P < 0.01). Our results suggest that the polymorphisms in the promoter and coding regions may modulate the promoter activity and gene expression of bovine ZBED6 in the skeletal muscles of these cattle breeds.

Genome-wide binding of the basic helix-loop-helix myogenic inhibitor musculin has substantial overlap with MyoD: implications for buffering activity

Background

Musculin (MSC) is a basic helix-loop-helix transcription factor that inhibits myogenesis during normal development and contributes to the differentiation defect in rhabdomyosarcoma. As one of many transcription factors that impede myogenesis, its binding on a genome-wide scale relative to the widespread binding of the myogenic factor MyoD is unknown.

Methods

Chromatin immunoprecipitation coupled to high-throughput sequencing was performed for endogenous MSC in rhabdomyosarcoma cells and its binding was compared to that of MyoD in the same type of cells.

Results

MSC binds throughout the genome, in a pattern very similar to MyoD. Its binding overlaps strongly with regions enriched for acetylated histone H4, as well as regions that score high for DNase hypersensitivity in human myoblasts. In contrast to MyoD, MSC has a more relaxed binding sequence preference in the nucleotides that flank the core E-box motif.

Conclusions

The myogenic inhibitor MSC binds throughout the genome of rhabdomyosarcoma cells, in a pattern highly similar to that of MyoD, suggesting a broad role in buffering the activity of MyoD in development and rhabdomyosarcomas.

Genome wide analysis reveals Zic3 interaction with distal regulatory elements of stage specific developmental genes in zebrafish

Zic3 regulates early embryonic patterning in vertebrates. Loss of Zic3 function is known to disrupt gastrulation, left-right patterning, and neurogenesis. However, molecular events downstream of this transcription factor are poorly characterized. Here we use the zebrafish as a model to study the developmental role of Zic3 in vivo, by applying a combination of two powerful genomics approaches – ChIP-seq and microarray. Besides confirming direct regulation of previously implicated Zic3 targets of the Nodal and canonical Wnt pathways, analysis of gastrula stage embryos uncovered a number of novel candidate target genes, among which were members of the non-canonical Wnt pathway and the neural pre-pattern genes. A similar analysis in zic3-expressing cells obtained by FACS at segmentation stage revealed a dramatic shift in Zic3 binding site locations and identified an entirely distinct set of target genes associated with later developmental functions such as neural development. We demonstrate cis-regulation of several of these target genes by Zic3 using in vivo enhancer assay. Analysis of Zic3 binding sites revealed a distribution biased towards distal intergenic regions, indicative of a long distance regulatory mechanism; some of these binding sites are highly conserved during evolution and act as functional enhancers. This demonstrated that Zic3 regulation of developmental genes is achieved predominantly through long distance regulatory mechanism and revealed that developmental transitions could be accompanied by dramatic changes in regulatory landscape.

The Zic3 transcription factor regulates early embryonic patterning, and the loss of its function leads to defects in left-right body asymmetry. Previous studies have only identified a small number of Zic3 targets, which renders the molecular mechanism underlying its activity insufficiently understood. Utilizing two genomics technologies, next generation sequencing and microarray, we profile the genome-wide binding sites of Zic3 and identified its target genes in the developing zebrafish embryo. Our results show that Zic3 regulates its target genes predominantly through regulatory elements located far from promoters. Among the targets of Zic3 are the Nodal and Wnt pathways known to regulate gastrulation and left-right body asymmetry, as well as neural pre-pattern genes regulating proliferation of neural progenitors. Using enhancer activity assay, we further show that genomic regions bound by Zic3 function as enhancers. Our study provides a genome-wide view of the regulatory landscape of Zic3 and its changes during vertebrate development.

Discriminative motif analysis of high-throughput dataset

MOTIVATION

High-throughput ChIP-seq studies typically identify thousands of peaks for a single transcription factor (TF). It is common for traditional motif discovery tools to predict motifs that are statistically significant against a naïve background distribution but are of questionable biological relevance.

RESULTS

We describe a simple yet effective algorithm for discovering differential motifs between two sequence datasets that is effective in eliminating systematic biases and scalable to large datasets. Tested on 207 ENCODE ChIP-seq datasets, our method identifies correct motifs in 78% of the datasets with known motifs, demonstrating improvement in both accuracy and efficiency compared with DREME, another state-of-art discriminative motif discovery tool. More interestingly, on the remaining more challenging datasets, we identify common technical or biological factors that compromise the motif search results and use advanced features of our tool to control for these factors. We also present case studies demonstrating the ability of our method to detect single base pair differences in DNA specificity of two similar TFs. Lastly, we demonstrate discovery of key TF motifs involved in tissue specification by examination of high-throughput DNase accessibility data.

AVAILABILITY

The motifRG package is publically available via the bioconductor repository.

CONTACT

yzizhen@fhcrc.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Hierarchical mechanisms for direct reprogramming of fibroblasts to neurons

Direct lineage reprogramming is a promising approach for human disease modeling and regenerative medicine, with poorly understood mechanisms. Here, we reveal a hierarchical mechanism in the direct conversion of fibroblasts into induced neuronal (iN) cells mediated by the transcription factors Ascl1, Brn2, and Myt1l. Ascl1 acts as an "on-target" pioneer factor by immediately occupying most cognate genomic sites in fibroblasts. In contrast, Brn2 and Myt1l do not access fibroblast chromatin productively on their own; instead, Ascl1 recruits Brn2 to Ascl1 sites genome wide. A unique trivalent chromatin signature in the host cells predicts the permissiveness for Ascl1 pioneering activity among different cell types. Finally, we identified Zfp238 as a key Ascl1 target gene that can partially substitute for Ascl1 during iN cell reprogramming. Thus, a precise match between pioneer factors and the chromatin context at key target genes is determinative for transdifferentiation to neurons and likely other cell types.

Emerin and histone deacetylase 3 (HDAC3) cooperatively regulate expression and nuclear positions of MyoD, Myf5, and Pax7 genes during myogenesis

The spatial organization of chromatin is critical in establishing cell-type dependent gene expression programs. The inner nuclear membrane protein emerin has been implicated in regulating global chromatin architecture. We show emerin associates with genomic loci of muscle differentiation promoting factors in murine myogenic progenitors, including Myf5 and MyoD. Prior to their transcriptional activation Myf5 and MyoD loci localized to the nuclear lamina in proliferating progenitors and moved to the nucleoplasm upon transcriptional activation during differentiation. The Pax7 locus, which is transcribed in proliferating progenitors, localized to the nucleoplasm and Pax7 moved to the nuclear lamina upon repression during differentiation. Localization of Myf5, MyoD, and Pax7 to the nuclear lamina and proper temporal expression of these genes required emerin and HDAC3. Interestingly, activation of HDAC3 catalytic activity rescued both Myf5 localization to the nuclear lamina and its expression. Collectively, these data support a model whereby emerin facilitates repressive chromatin formation at the nuclear lamina by activating the catalytic activity of HDAC3 to regulate the coordinated spatiotemporal expression of myogenic differentiation genes.

DNA damage-activated ABL-MyoD signaling contributes to DNA repair in skeletal myoblasts

Previous works have established a unique function of MyoD in the control of muscle gene expression during DNA damage response in myoblasts. Phosphorylation by DNA damage-activated ABL tyrosine kinase transiently inhibits MyoD-dependent activation of transcription in response to genotoxic stress. We show here that ABL-MyoD signaling is also an essential component of the DNA repair machinery in myoblasts exposed to genotoxic stress. DNA damage promoted the recruitment of MyoD to phosphorylated Nbs1 (pNbs1)-containing repair foci, and this effect was abrogated by either ABL knockdown or the ABL kinase inhibitor imatinib. Upon DNA damage, MyoD and pNbs1 were detected on the chromatin to MyoD target genes without activating transcription. DNA damage-mediated tyrosine phosphorylation was required for MyoD recruitment to target genes, as the ABL phosphorylation-resistant MyoD mutant (MyoD Y30F) failed to bind the chromatin following DNA damage, while retaining the ability to activate transcription in response to differentiation signals. Moreover, MyoD Y30F exhibited an impaired ability to promote repair in a heterologous system, as compared with MyoD wild type (WT). Consistently, MyoD-null satellite cells (SCs) displayed impaired DNA repair that was rescued by reintroduction of MyoD WT but not by MyoD Y30F. In addition, inhibition of ABL kinase prevented MyoD WT-mediated rescue of DNA repair in MyoD-null SCs. These results identify an unprecedented contribution of MyoD to DNA repair and suggest that ABL-MyoD signaling coordinates DNA repair and transcription in myoblasts.

Divergent transcriptional regulatory logic at the intersection of tissue growth and developmental patterning

The Yorkie/Yap transcriptional coactivator is a well-known regulator of cellular proliferation in both invertebrates and mammals. As a coactivator, Yorkie (Yki) lacks a DNA binding domain and must partner with sequence-specific DNA binding proteins in the nucleus to regulate gene expression; in Drosophila, the developmental regulators Scalloped (Sd) and Homothorax (Hth) are two such partners. To determine the range of target genes regulated by these three transcription factors, we performed genome-wide chromatin immunoprecipitation experiments for each factor in both the wing and eye-antenna imaginal discs. Strong, tissue-specific binding patterns are observed for Sd and Hth, while Yki binding is remarkably similar across both tissues. Binding events common to the eye and wing are also present for Sd and Hth; these are associated with genes regulating cell proliferation and “housekeeping” functions, and account for the majority of Yki binding. In contrast, tissue-specific binding events for Sd and Hth significantly overlap enhancers that are active in the given tissue, are enriched in Sd and Hth DNA binding sites, respectively, and are associated with genes that are consistent with each factor's previously established tissue-specific functions. Tissue-specific binding events are also significantly associated with Polycomb targeted chromatin domains. To provide mechanistic insights into tissue-specific regulation, we identify and characterize eye and wing enhancers of the Yki-targeted bantam microRNA gene and demonstrate that they are dependent on direct binding by Hth and Sd, respectively. Overall these results suggest that both Sd and Hth use distinct strategies – one shared between tissues and associated with Yki, the other tissue-specific, generally Yki-independent and associated with developmental patterning – to regulate distinct gene sets during development.

The Hippo tumor suppressor pathway controls proliferation in a tissue-nonspecific fashion in Drosophila epithelial progenitor tissues via the transcriptional coactivator Yorkie (Yki). However, despite the tissue-nonspecific role that Yki plays in tissue growth, the transcription factors that recruit Yki to DNA, most notably Scalloped (Sd) and Homothorax (Hth), are important regulators of developmental patterning with many tissue-specific functions. Thus, these three transcriptional regulators – Yki, Sd, and Hth – provide a model for exploring the properties of protein-DNA interactions that regulate both tissue-shared and tissue-specific functions. With this goal in mind, we identified the positions in the fly genome that are bound by Yki, Sd, and Hth in the progenitors of the wing and eye-antenna structures of the fly. These data not only provide a global view of the Yki gene regulatory network, they reveal an unusual amount of tissue specificity in the genomic regions targeted by Sd and Hth, but not Yki. The data also reveal that tissue-specific binding is very likely to overlap tissue-specific enhancer regions, provide important clues for how tissue-specific Sd and Hth binding occurs, and support the idea that gene regulatory networks are plastic, with spatial differences in binding significantly impacting network structures.

The Igf2/H19 muscle enhancer is an active transcriptional complex

In eukaryotic cells, gene expression is mediated by enhancer activation of RNA polymerase at distant promoters. Recently, distinctions between enhancers and promoters have been blurred by the discovery that enhancers are associated with RNA polymerase and are sites of RNA synthesis. Here, we present an analysis of the insulin-like growth factor 2/H19 muscle enhancer. This enhancer includes a short conserved core element that is organized into chromatin typical of mammalian enhancers, binds tissue-specific transcription factors and functions on its own in vitro to activate promoter transcription. However, in a chromosomal context, this element is not sufficient to activate distant promoters. Instead, enhancer function also requires transcription in cis of a long non-coding RNA, Nctc1. Thus, the insulin-like growth factor 2/H19 enhancer is an active transcriptional complex whose own transcription is essential to its function.

Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

eRNAs promote transcription by establishing chromatin accessibility at defined genomic loci

Transcription factors and DNA regulatory binding motifs are fundamental components of the gene regulatory network. Here, by using genome-wide binding profiling, we show extensive occupancy of transcription factors of myogenesis (MyoD and Myogenin) at extragenic enhancer regions coinciding with RNA synthesis (i.e., eRNA). In particular, multiple regions were transcribed to eRNA within the regulatory region of MYOD1, including previously characterized distal regulatory regions (DRR) and core enhancer (CE). While (CE)RNA enhanced RNA polymerase II (Pol II) occupancy and transcription at MYOD1, (DRR)RNA acted to activate the downstream myogenic genes. The deployment of transcriptional machinery to appropriate loci is contingent on chromatin accessibility, a rate-limiting step preceding Pol II assembly. By nuclease sensitivity assay, we found that eRNAs regulate genomic access of the transcriptional complex to defined regulatory regions. In conclusion, our data suggest that eRNAs contribute to establishing a cell-type-specific transcriptional circuitry by directing chromatin-remodeling events.

Genome-wide survey by ChIP-seq reveals YY1 regulation of lincRNAs in skeletal myogenesis

Skeletal muscle differentiation is orchestrated by a network of transcription factors, epigenetic regulators, and non-coding RNAs. The transcription factor Yin Yang 1 (YY1) silences multiple target genes in myoblasts (MBs) by recruiting Ezh2 (Enhancer of Zeste Homologue2). To elucidate genome-wide YY1 binding in MBs, we performed chromatin immunoprecipitation (ChIP)-seq and found 1820 specific binding sites in MBs with a large portion residing in intergenic regions. Detailed analysis demonstrated that YY1 acts as an activator for many loci in addition to its known repressor function. No significant co-occupancy was found between YY1 and Ezh2, suggesting an additional Ezh2-independent function for YY1 in MBs. Further analysis of intergenic binding sites showed that YY1 potentially regulates dozens of large intergenic non-coding RNAs (lincRNAs), whose function in myogenesis is underexplored. We characterized a novel muscle-associated lincRNA (Yam-1) that is positively regulated by YY1. Yam-1 is downregulated upon differentiation and acts as an inhibitor of myogenesis. We demonstrated that Yam-1 functions through in cis regulation of miR-715, which in turn targets Wnt7b. Our findings not only provide the first genome-wide picture of YY1 association in muscle cells, but also uncover the functional role of lincRNA Yam-1.

DNA methylation and differentiation: HOX genes in muscle cells

BACKGROUND

Tight regulation of homeobox genes is essential for vertebrate development. In a study of genome-wide differential methylation, we recently found that homeobox genes, including those in the HOX gene clusters, were highly overrepresented among the genes with hypermethylation in the skeletal muscle lineage. Methylation was analyzed by reduced representation bisulfite sequencing (RRBS) of postnatal myoblasts, myotubes and adult skeletal muscle tissue and 30 types of non-muscle-cell cultures or tissues.

RESULTS

In this study, we found that myogenic hypermethylation was present in specific subregions of all four HOX gene clusters and was associated with various chromatin epigenetic features. Although the 3' half of the HOXD cluster was silenced and enriched in polycomb repression-associated H3 lysine 27 trimethylation in most examined cell types, including myoblasts and myotubes, myogenic samples were unusual in also displaying much DNA methylation in this region. In contrast, both HOXA and HOXC clusters displayed myogenic hypermethylation bordering a central region containing many genes preferentially expressed in myogenic progenitor cells and consisting largely of chromatin with modifications typical of promoters and enhancers in these cells. A particularly interesting example of myogenic hypermethylation was HOTAIR, a HOXC noncoding RNA gene, which can silence HOXD genes in trans via recruitment of polycomb proteins. In myogenic progenitor cells, the preferential expression of HOTAIR was associated with hypermethylation immediately downstream of the gene. Other HOX gene regions also displayed myogenic DNA hypermethylation despite being moderately expressed in myogenic cells. Analysis of representative myogenic hypermethylated sites for 5-hydroxymethylcytosine revealed little or none of this base, except for an intragenic site in HOXB5 which was specifically enriched in this base in skeletal muscle tissue, whereas myoblasts had predominantly 5-methylcytosine at the same CpG site.

CONCLUSIONS

Our results suggest that myogenic hypermethylation of HOX genes helps fine-tune HOX sense and antisense gene expression through effects on 5' promoters, intragenic and intergenic enhancers and internal promoters. Myogenic hypermethylation might also affect the relative abundance of different RNA isoforms, facilitate transcription termination, help stop the spread of activation-associated chromatin domains and stabilize repressive chromatin structures.

Identification of cis-regulatory modules in promoters of human genes exploiting mutual positioning of transcription factors

In higher organisms, gene regulation is controlled by the interplay of non-random combinations of multiple transcription factors (TFs). Although numerous attempts have been made to identify these combinations, important details, such as mutual positioning of the factors that have an important role in the TF interplay, are still missing. The goal of the present work is in silico mapping of some of such associating factors based on their mutual positioning, using computational screening. We have selected the process of myogenesis as a study case, and we focused on TF combinations involving master myogenic TF Myogenic differentiation (MyoD) with other factors situated at specific distances from it. The results of our work show that some muscle-specific factors occur together with MyoD within the range of ±100 bp in a large number of promoters. We confirm co-occurrence of the MyoD with muscle-specific factors as described in earlier studies. However, we have also found novel relationships of MyoD with other factors not specific for muscle. Additionally, we have observed that MyoD tends to associate with different factors in proximal and distal promoter areas. The major outcome of our study is establishing the genome-wide connection between biological interactions of TFs and close co-occurrence of their binding sites.

Epigenetic control of skeletal muscle regeneration: Integrating genetic determinants and environmental changes

During embryonic development, pluripotent cells are genetically committed to specific lineages by the expression of cell-type-specific transcriptional activators that direct the formation of specialized tissues and organs in response to developmental cues. Chromatin-modifying proteins are emerging as essential components of the epigenetic machinery, which establishes the nuclear landscape that ultimately determines the final identity and functional specialization of adult cells. Recent evidence has revealed that discrete populations of adult cells can retain the ability to adopt alternative cell fates in response to environmental cues. These cells include conventional adult stem cells and a still poorly defined collection of cell types endowed with facultative phenotype and functional plasticity. Under physiological conditions or adaptive states, these cells cooperate to support tissue and organ homeostasis, and to promote growth or compensatory regeneration. However, during chronic diseases and aging these cells can adopt a pathological phenotype and mediate maladaptive responses, such as the formation of fibrotic scars and fat deposition that progressively replaces structural and functional units of tissues and organs. The molecular determinants of these phenotypic transitions are only emerging from recent studies that reveal how dynamic chromatin states can generate flexible epigenetic landscapes, which confer on cells the ability to retain partial pluripotency and adapt to environmental changes. This review summarizes our current knowledge on the role of the epigenetic machinery as a 'filter' between genetic commitment and environmental signals in cell types that can alternatively promote skeletal muscle regeneration or fibro-adipogenic degeneration.

Chromatin tandem affinity purification sequencing

Chromatin immunoprecipitation coupled with ultra-high-throughput sequencing (ChIP-seq) is a widely used method for mapping the interactions of proteins with DNA. However, the requirements for ChIP-grade antibodies impede wider application of this method, and variations in results can be high owing to differences in affinity and cross-reactivity of antibodies. Therefore, we developed chromatin tandem affinity purification (ChTAP) as an effective alternative to ChIP. Through the use of affinity tags and reagents that are identical for all proteins investigated, ChTAP enables one to directly compare the binding between different transcription factors and to directly assess the background in control experiments. Thus, ChTAP-seq can be used to rapidly map the genome-wide binding of multiple DNA-binding proteins in a wide range of cell types. ChTAP can be completed in 3-4 d, starting from cross-linking of chromatin to purification of ChIP DNA.

Massively parallel in vivo enhancer assay reveals that highly local features determine the cis-regulatory function of ChIP-seq peaks

Transcription factors (TFs) recognize short sequence motifs that are present in millions of copies in large eukaryotic genomes. TFsmust distinguish their target binding sites from a vast genomic excess of spurious motif occurrences; however, it is unclear whether functional sites are distinguished from nonfunctional motifs by local primary sequence features or by the larger genomic context in which motifs reside. We used a massively parallel enhancer assay in living mouse retinas to compare 1,300 sequences bound in the genome by the photoreceptor transcription factor Cone-rod homeobox (Crx), to 3,000 control sequences. We found that very short sequences bound in the genome by Crx activated transcription at high levels, whereas unbound genomic regions with equal numbers of Crx motifs did not activate above background levels, even when liberated from their larger genomic context. High local GC content strongly distinguishes bound motifs from unbound motifs across the entire genome. Our results show that the cis-regulatory potential of TF-bound DNA is determined largely by highly local sequence features and not by genomic context.

Six1 regulates MyoD expression in adult muscle progenitor cells

Quiescent satellite cells are myogenic progenitors that enable regeneration of skeletal muscle. One of the early events of satellite cell activation following myotrauma is the induction of the myogenic regulatory factor MyoD, which eventually induces terminal differentiation and muscle function gene expression. The purpose of this study was to elucidate the mechanism by which MyoD is induced during activation of satellite cells in mouse muscle undergoing regeneration. We show that Six1, a transcription factor essential for embryonic myogenesis, also regulates MyoD expression in muscle progenitor cells. Six1 knock-down by RNA interference leads to decreased expression of MyoD in myoblasts. Chromatin immunoprecipitation assays reveal that Six1 binds the Core Enhancer Region of MyoD. Further, transcriptional reporter assays demonstrate that Core Enhancer Region reporter gene activity in myoblasts and in regenerating muscle depends on the expression of Six1 and on Six1 binding sites. Finally, we provide evidence indicating that Six1 is required for the proper chromatin structure at the Core Enhancer Region, as well as for MyoD binding at its own enhancer. Together, our results reveal that MyoD expression in satellite cells depends on Six1, supporting the idea that Six1 plays an important role in adult myogenesis, in addition to its role in embryonic muscle formation.

Signal-dependent fra-2 regulation in skeletal muscle reserve and satellite cells

Activator protein-1 (AP-1) is a ubiquitous transcription factor that paradoxically also has some tissue-specific functions. In skeletal muscle cells, we document that the AP-1 subunit, Fra-2, is expressed in the resident stem cells (Pax7-positive satellite cells) and also in the analogous undifferentiated ‘reserve' cell population in myogenic cultures, but not in differentiated myofiber nuclei. Silencing of Fra-2 expression enhances the expression of differentiation markers such as muscle creatine kinase and myosin heavy chain, indicating a possible role of Fra-2 in undifferentiated myogenic progenitor cells. We observed that Fra-2 is a target of cytokine-mediated extracellular signal-regulated kinase-1/2 signaling in cultured muscle cells, and extensive mass spectrometry and mutational analysis identified S320 and T322 as regulators of Fra-2 protein stability. Interestingly, Fra-2 S320 phosphorylation occurs transiently in activated satellite cells and is extinguished in myogenin-positive differentiating cells. Thus, cytokine-mediated Fra-2 expression and stabilization is linked to regulation of myogenic progenitor cells having implications for the molecular regulation of adult muscle stem cells and skeletal muscle regeneration.

The role of MyoD1 and histone modifications in the activation of muscle enhancers

MyoD1 is a key regulator that orchestrates skeletal muscle differentiation through the regulation of gene expression. Although many studies have focused on its role in transcriptional control at gene promoters, less is known regarding the role of MyoD1 in the assembly of active enhancers. Here, we discuss novel data that point to the ability of MyoD1 to mediate the assembly of active enhancers that augment the transcription of genes essential for muscle development and lineage specification. Based on genome-wide studies of epigenetic marks that typify active enhancers, we recently identified the compendium of distal regulatory elements that dictate transcriptional programs during myogenesis. Superimposition of MyoD1 binding sites upon the locations of muscle enhancers revealed its unequivocal binding to a core region of nearly a third of condition-specific muscle enhancers. Further studies exploring deposition of enhancer-related epigenetic marks in myoblasts lacking MyoD1 demonstrate the dependence of muscle enhancer assembly on the presence of MyoD1. We propose a model wherein MyoD1 mediates recruitment of Set7, H3K4me1, H3K27ac, p300, and RNAP II to MyoD1-bound enhancers to establish condition-specific activation of muscle genes. Moreover, muscle enhancers are modulated through coordinated binding of transcription factors, including c-Jun, Jdp2, Meis, and Runx1, which are recruited to muscle enhancers in a MyoD1-dependent manner. Thus, MyoD1 and enhancer-associated transcription factors function coordinately to assemble and regulate enhancers, thereby augmenting expression of muscle-related genes.

Chromatin modifications as determinants of muscle stem cell quiescence and chronological aging

The ability to maintain quiescence is critical for the long-term maintenance of a functional stem cell pool. To date, the epigenetic and transcriptional characteristics of quiescent stem cells and how they change with age remain largely unknown. In this study, we explore the chromatin features of adult skeletal muscle stem cells, or satellite cells (SCs), which reside predominantly in a quiescent state in fully developed limb muscles of both young and aged mice. Using a ChIP-seq approach to obtain global epigenetic profiles of quiescent SCs (QSCs), we show that QSCs possess a permissive chromatin state in which few genes are epigenetically repressed by Polycomb group (PcG)-mediated histone 3 lysine 27 trimethylation (H3K27me3), and a large number of genes encoding regulators that specify nonmyogenic lineages are demarcated by bivalent domains at their transcription start sites (TSSs). By comparing epigenetic profiles of QSCs from young and old mice, we also provide direct evidence that, with age, epigenetic changes accumulate and may lead to a functional decline in quiescent stem cells. These findings highlight the importance of chromatin mapping in understanding unique features of stem cell identity and stem cell aging.

Skeletal muscle programming and re-programming

The discovery of the transcription factor MyoD and its ability to induce muscle differentiation was the first demonstration of genetically programmed cell transdifferentiation. MyoD functions by activating a feed-forward circuit to regulate muscle gene expression. This requires binding to specific E-boxes throughout the genome, followed by recruitment of chromatin modifying complexes and transcription machinery. MyoD binding can be modified by both cooperative factors and inhibitors, including microRNAs that may serve as important developmental switches. Recent studies indicate that epigenetic regulation of MyoD binding sites is another important mechanism for controlling MyoD activity, which may ultimately limit its ability to induce transdifferentiation to cells with permissive epigenetic 'landscapes.'

ASCL1 reprograms mouse Muller glia into neurogenic retinal progenitors

Non-mammalian vertebrates have a robust ability to regenerate injured retinal neurons from Müller glia (MG) that activate the gene encoding the proneural factor Achaete-scute homolog 1 (Ascl1; also known as Mash1 in mammals) and de-differentiate into progenitor cells. By contrast, mammalian MG have a limited regenerative response and fail to upregulate Ascl1 after injury. To test whether ASCL1 could restore neurogenic potential to mammalian MG, we overexpressed ASCL1 in dissociated mouse MG cultures and intact retinal explants. ASCL1-infected MG upregulated retinal progenitor-specific genes and downregulated glial genes. Furthermore, ASCL1 remodeled the chromatin at its targets from a repressive to an active configuration. MG-derived progenitors differentiated into cells that exhibited neuronal morphologies, expressed retinal subtype-specific neuronal markers and displayed neuron-like physiological responses. These results indicate that a single transcription factor, ASCL1, can induce a neurogenic state in mature MG.

Sources for skeletal muscle repair: from satellite cells to reprogramming

Skeletal muscle regeneration is the process that ensures tissue repair after damage by injury or in degenerative diseases such as muscular dystrophy. Satellite cells, the adult skeletal muscle progenitor cells, are commonly considered to be the main cell type involved in skeletal muscle regeneration. Their mechanism of action in this process is extensively characterized. However, evidence accumulated in the last decade suggests that other cell types may participate in skeletal muscle regeneration. Although their actual contribution to muscle formation and regeneration is still not clear; if properly manipulated, these cells may become new suitable and powerful sources for cell therapy of skeletal muscle degenerative diseases. Mesoangioblasts, vessel associated stem/progenitor cells with high proliferative, migratory and myogenic potential, are very good candidates for clinical applications and are already in clinical experimentation. In addition, pluripotent stem cells are very promising sources for regeneration of most tissues, including skeletal muscle. Conditions such as muscle cachexia or aging that severely alter homeostasis may be counteracted by transplantation of donor and/or recruitment and activation of resident muscle stem/progenitor cells. Advantages and limitations of different cell therapy approaches will be discussed.

Six homeoproteins directly activate Myod expression in the gene regulatory networks that control early myogenesis

In mammals, several genetic pathways have been characterized that govern engagement of multipotent embryonic progenitors into the myogenic program through the control of the key myogenic regulatory gene Myod. Here we demonstrate the involvement of Six homeoproteins. We first targeted into a Pax3 allele a sequence encoding a negative form of Six4 that binds DNA but cannot interact with essential Eya co-factors. The resulting embryos present hypoplasic skeletal muscles and impaired Myod activation in the trunk in the absence of Myf5/Mrf4. At the axial level, we further show that Myod is still expressed in compound Six1/Six4:Pax3 but not in Six1/Six4:Myf5 triple mutant embryos, demonstrating that Six1/4 participates in the Pax3-Myod genetic pathway. Myod expression and head myogenesis is preserved in Six1/Six4:Myf5 triple mutant embryos, illustrating that upstream regulators of Myod in different embryonic territories are distinct. We show that Myod regulatory regions are directly controlled by Six proteins and that, in the absence of Six1 and Six4, Six2 can compensate.

The onset of skeletal muscle formation is controlled by complex gene regulatory networks. By manipulation of these genetic pathways in the mouse embryo, we have examined the interplay between genes encoding the transcriptional regulator Pax3; the major myogenic determination proteins Myf5, Mrf4, and Myod; as well as genes encoding homeodomain proteins Six1 and Six4. In the absence of Myf5 and Six1/4, Myod expression is compromised. We demonstrate that key regulatory elements of the Myod gene are directly targeted by Six factors, including Six2, which is unexpectedly upregulated in the absence of Six1 and Six4. This work therefore reveals new aspects of the gene regulatory networks that control myogenesis.

The Scaffold attachment factor b1 (Safb1) regulates myogenic differentiation by facilitating the transition of myogenic gene chromatin from a repressed to an activated state

The regulation of skeletal muscle gene expression during myogenesis is mediated by lineage-specific transcription factors in combination with numerous cofactors, many of which modify chromatin structure. However, the involvement of scaffolding proteins that organize chromatin and chromatin-associated regulatory proteins has not extensively been explored in myogenic differentiation. Here, we report that Scaffold attachment factor b1 (Safb1), primarily associated with transcriptional repression, functions as a positive regulator of myogenic differentiation. Knockdown of Safb1 inhibited skeletal muscle marker gene expression and differentiation in cultured C2C12 myoblasts. In contrast, over-expression resulted in the premature expression of critical muscle structural proteins and formation of enlarged thickened myotubes. Safb1 co-immunoprecipitated with MyoD and was co-localized on myogenic promoters. Upon Safb1 knockdown, the repressive H3K27me3 histone mark and binding of the Polycomb histone methyltransferase Ezh2 persisted at differentiation-dependent gene promoters. In contrast, the appearance of histone marks and regulators associated with myogenic gene activation, such as myogenin and the SWI/SNF chromatin remodelling enzyme ATPase, Brg1, was blocked. These results indicate that the scaffold protein Safb1 contributes to the activation of skeletal muscle gene expression during myogenic differentiation by facilitating the transition of promoter sequences from a repressive chromatin structure to one that is transcriptionally permissive.

Pbx and Prdm1a transcription factors differentially regulate subsets of the fast skeletal muscle program in zebrafish

The basic helix–loop–helix factor Myod initiates skeletal muscle differentiation by directly and sequentially activating sets of muscle differentiation genes, including those encoding muscle contractile proteins. We hypothesize that Pbx homeodomain proteins direct Myod to a subset of its transcriptional targets, in particular fast-twitch muscle differentiation genes, thereby regulating the competence of muscle precursor cells to differentiate. We have previously shown that Pbx proteins bind with Myod on the promoter of the zebrafish fast muscle gene mylpfa and that Pbx proteins are required for Myod to activate mylpfa expression and the fast-twitch muscle-specific differentiation program in zebrafish embryos. Here we have investigated the interactions of Pbx with another muscle fiber-type regulator, Prdm1a, a SET-domain DNA-binding factor that directly represses mylpfa expression and fast muscle differentiation. The prdm1a mutant phenotype, early and increased fast muscle differentiation, is the opposite of the Pbx-null phenotype, delayed and reduced fast muscle differentiation. To determine whether Pbx and Prdm1a have opposing activities on a common set of genes, we used RNA-seq analysis to globally assess gene expression in zebrafish embryos with single- and double-losses-of-function for Pbx and Prdm1a. We find that the levels of expression of certain fast muscle genes are increased or approximately wild type in pbx2/4-MO;prdm1a−/− embryos, suggesting that Pbx activity normally counters the repressive action of Prdm1a for a subset of the fast muscle program. However, other fast muscle genes require Pbx but are not regulated by Prdm1a. Thus, our findings reveal that subsets of the fast muscle program are differentially regulated by Pbx and Prdm1a. Our findings provide an example of how Pbx homeodomain proteins act in a balance with other transcription factors to regulate subsets of a cellular differentiation program.

Comparison of endogenous and overexpressed MyoD shows enhanced binding of physiologically bound sites

BACKGROUND

Transcription factor overexpression is common in biological experiments and transcription factor amplification is associated with many cancers, yet few studies have directly compared the DNA-binding profiles of endogenous versus overexpressed transcription factors.

METHODS

We analyzed MyoD ChIP-seq data from C2C12 mouse myotubes, primary mouse myotubes, and mouse fibroblasts differentiated into muscle cells by overexpression of MyoD and compared the genome-wide binding profiles and binding site characteristics of endogenous and overexpressed MyoD.

RESULTS

Overexpressed MyoD bound to the same sites occupied by endogenous MyoD and possessed the same E-box sequence preference and co-factor site enrichments, and did not bind to new sites with distinct characteristics.

CONCLUSIONS

Our data demonstrate a robust fidelity of transcription factor binding sites over a range of expression levels and that increased amounts of transcription factor increase the binding at physiologically bound sites.

miR-411 is up-regulated in FSHD myoblasts and suppresses myogenic factors

BACKGROUND

Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal dominant muscle disorder, which is linked to the contraction of the D4Z4 array at chromosome 4q35. Recent studies suggest that this shortening of the D4Z4 array leads to aberrant expression of double homeobox protein 4 (DUX4) and causes FSHD. In addition, misregulation of microRNAs (miRNAs) has been reported in muscular dystrophies including FSHD. In this study, we identified a miRNA that is differentially expressed in FSHD myoblasts and investigated its function.

METHODS

To identify misregulated miRNAs and their potential targets in FSHD myoblasts, we performed expression profiling of both miRNA and mRNA using TaqMan Human MicroRNA Arrays and Affymetrix Human Genome U133A plus 2.0 microarrays, respectively. In addition, we over-expressed miR-411 in C₂C₁₂ cells to determine the effect of miR-411 on myogenic markers.

RESULTS

Using miRNA and mRNA expression profiling, we identified 8 miRNAs and 1,502 transcripts that were differentially expressed in FSHD myoblasts during cell proliferation. One of the 8 differentially expressed miRNAs, miR-411, was validated by quantitative RT-PCR in both primary (2.1 fold, p<0.01) and immortalized (2.7 fold, p<0.01) myoblasts. In situ hybridization showed cytoplasmic localization of miR-411 in FSHD myoblasts. By analyzing both miRNA and mRNA data using Partek Genomics Suite, we identified 4 mRNAs potentially regulated by miR-411 including YY1 associated factor 2 (YAF2). The down-regulation of YAF2 in immortalized myoblasts was validated by immunoblotting (-3.7 fold, p<0.01). C₂C₁₂ cells were transfected with miR-411 to determine whether miR-411 affects YAF2 expression in myoblasts. The results showed that over-expression of miR-411 reduced YAF2 mRNA expression. In addition, expression of myogenic markers including Myod, myogenin, and myosin heavy chain 1 (Myh1) were suppressed by miR-411.

CONCLUSIONS

The study demonstrated that miR-411 was differentially expressed in FSHD myoblasts and may play a role in regulating myogenesis.

Intrinsic epigenetic regulation of the D4Z4 macrosatellite repeat in a transgenic mouse model for FSHD

Facioscapulohumeral dystrophy (FSHD) is a progressive muscular dystrophy caused by decreased epigenetic repression of the D4Z4 macrosatellite repeats and ectopic expression of DUX4, a retrogene encoding a germline transcription factor encoded in each repeat. Unaffected individuals generally have more than 10 repeats arrayed in the subtelomeric region of chromosome 4, whereas the most common form of FSHD (FSHD1) is caused by a contraction of the array to fewer than 10 repeats, associated with decreased epigenetic repression and variegated expression of DUX4 in skeletal muscle. We have generated transgenic mice carrying D4Z4 arrays from an FSHD1 allele and from a control allele. These mice recapitulate important epigenetic and DUX4 expression attributes seen in patients and controls, respectively, including high DUX4 expression levels in the germline, (incomplete) epigenetic repression in somatic tissue, and FSHD–specific variegated DUX4 expression in sporadic muscle nuclei associated with D4Z4 chromatin relaxation. In addition we show that DUX4 is able to activate similar functional gene groups in mouse muscle cells as it does in human muscle cells. These transgenic mice therefore represent a valuable animal model for FSHD and will be a useful resource to study the molecular mechanisms underlying FSHD and to test new therapeutic intervention strategies.

Facioscapulohumeral dystrophy (FSHD) is a progressive muscle disorder that is associated with contraction and chromatin relaxation of the D4Z4 macrosatellite repeat on chromosome 4q. Each unit of the repeat contains a copy of the primate-specific DUX4 retrogene, encoding a germline transcription factor that is repressed in somatic tissue. In FSHD, somatic repression of the DUX4 gene is compromised, leading to a variegated expression pattern of DUX4 in muscle cells. The complex (epi)genetic etiology of FSHD has long hampered the generation of a faithful animal model, and thus far the role of FSHD candidate genes has only been studied in model organisms by overexpression approaches. Here we present two transgenic mouse models containing either patient- or control-sized D4Z4 repeats. In our mice, the regulation of the FSHD locus is preserved in both lines, and only in the disease model somatic derepression and variegated expression of DUX4 is observed. These mice thus reflect many aspects of the complex regulation of DUX4 expression in humans. These models may therefore become valuable tools in understanding the in vivo regulation and function of DUX4, its role in FSHD, and the evaluation of therapeutic strategies.

Proteomic and functional analyses reveal the role of chromatin reader SFMBT1 in regulating epigenetic silencing and the myogenic gene program

SFMBT1 belongs to the malignant brain tumor domain-containing chromatin reader family that recognizes repressive histone marks and represses transcription. The biological functions and molecular basis underlying SFMBT1-mediated transcriptional repression are poorly elucidated. Here, our proteomic analysis revealed that SFMBT1 is associated with multiple transcriptional corepressor complexes, including CtBP/LSD1/HDAC complexes, polycomb repressive complexes, and malignant brain tumor family proteins, that collectively contribute to SFMBT1 repressor activity. During myogenesis, Sfmbt1 represses myogenic differentiation of cultured and primary myoblasts. Mechanistically, Sfmbt1 interacts with MyoD and mediates epigenetic silencing of MyoD target genes via recruitment of its associated corepressors and subsequent induction of epigenetic modifications and chromatin compaction. Therefore, our study identified novel mechanisms accounting for SFMBT1-mediated transcription repression and revealed an essential role of Sfmbt1 in regulating MyoD-mediated transcriptional silencing that is required for the maintenance of undifferentiated states of myogenic progenitor cells.

Satellite cells and the muscle stem cell niche

Adult skeletal muscle in mammals is a stable tissue under normal circumstances but has remarkable ability to repair after injury. Skeletal muscle regeneration is a highly orchestrated process involving the activation of various cellular and molecular responses. As skeletal muscle stem cells, satellite cells play an indispensible role in this process. The self-renewing proliferation of satellite cells not only maintains the stem cell population but also provides numerous myogenic cells, which proliferate, differentiate, fuse, and lead to new myofiber formation and reconstitution of a functional contractile apparatus. The complex behavior of satellite cells during skeletal muscle regeneration is tightly regulated through the dynamic interplay between intrinsic factors within satellite cells and extrinsic factors constituting the muscle stem cell niche/microenvironment. For the last half century, the advance of molecular biology, cell biology, and genetics has greatly improved our understanding of skeletal muscle biology. Here, we review some recent advances, with focuses on functions of satellite cells and their niche during the process of skeletal muscle regeneration.

Myod and H19-Igf2 locus interactions are required for diaphragm formation in the mouse

The myogenic regulatory factor Myod and insulin-like growth factor 2 (Igf2) have been shown to interact in vitro during myogenic differentiation. In order to understand how they interact in vivo, we produced double-mutant mice lacking both the Myod and Igf2 genes. Surprisingly, these mice display neonatal lethality due to severe diaphragm atrophy. Alteration of diaphragm muscle development occurs as early as 15.5 days post-coitum in the double-mutant embryos and leads to a defect in the terminal differentiation of muscle progenitor cells. A negative-feedback loop was detected between Myod and Igf2 in embryonic muscles. Igf2 belongs to the imprinted H19-Igf2 locus. Molecular analyses show binding of Myod on a mesodermal enhancer (CS9) of the H19 gene. Chromatin conformation capture experiments reveal direct interaction of CS9 with the H19 promoter, leading to increased H19 expression in the presence of Myod. In turn, the non-coding H19 RNA represses Igf2 expression in trans. In addition, Igf2 also negatively regulates Myod expression, possibly by reducing the expression of the Srf transcription factor, a known Myod activator. In conclusion, Igf2 and Myod are tightly co-regulated in skeletal muscles and act in parallel pathways in the diaphragm, where they affect the progression of myogenic differentiation. Igf2 is therefore an essential player in the formation of a functional diaphragm in the absence of Myod.

Differentiation-associated microRNAs antagonize the Rb-E2F pathway to restrict proliferation

Transcriptional regulation by Rb–E2F and posttranscriptional regulation by microRNAs control the expression of cell cycle and DNA replication genes and restrict cellular proliferation.

The cancer-associated loss of microRNA (miRNA) expression leads to a proliferative advantage and aggressive behavior through largely unknown mechanisms. Here, we exploit a model system that recapitulates physiological terminal differentiation and its reversal upon oncogene expression to analyze coordinated mRNA/miRNA responses. The cell cycle reentry of myotubes, forced by the E1A oncogene, was associated with a pattern of mRNA/miRNA modulation that was largely reciprocal to that induced during the differentiation of myoblasts into myotubes. The E1A-induced mRNA response was preponderantly Retinoblastoma protein (Rb)-dependent. Conversely, the miRNA response was mostly Rb-independent and exerted through tissue-specific factors and Myc. A subset of these miRNAs (miR-1, miR-34, miR-22, miR-365, miR-29, miR-145, and Let-7) was shown to coordinately target Rb-dependent cell cycle and DNA replication mRNAs. Thus, a dual level of regulation—transcriptional regulation via Rb–E2F and posttranscriptional regulation via miRNAs—confers robustness to cell cycle control and provides a molecular basis to understand the role of miRNA subversion in cancer.

Genome-wide identification of enhancers in skeletal muscle: the role of MyoD1

To identify the compendium of distal regulatory elements that govern myogenic differentiation, we generated chromatin state maps based on histone modifications and recruitment of factors that typify enhancers in myoblasts and myotubes. We found a striking concordance between the locations of these newly defined enhancers, MyoD1-binding events, and noncoding RNA transcripts. These enhancers recruit several sequence-specific transcription factors in a spatially constrained manner around MyoD1-binding sites. Remarkably, MyoD1-null myoblasts show a wholesale loss of recruitment of these factors as well as diminished monomethylation of H3K4 (H3K4me1) and acetylation of H3K27 (H3K27ac) and reduced recruitment of Set7, an H3K4 monomethylase. Surprisingly, we found that H3K4me1, but not H3K27ac, could be restored by re-expression of MyoD1 in MyoD1(-/-) myoblasts, although re-expression of this factor in MyoD1-null myotubes restored both histone modifications. Our studies identified a role for MyoD1 in condition-specific enhancer assembly through recruitment of transcription factors and histone-modifying enzymes that shape muscle differentiation.

Transcription factor and chromatin features predict genes associated with eQTLs

Cell type-specific gene expression in humans involves complex interactions between regulatory factors and DNA at enhancers and promoters. Mapping studies for expression quantitative trait loci (eQTLs), transcription factors (TFs) and chromatin markers have become widely used tools for identifying gene regulatory elements, but prediction of target genes remains a major challenge. Here, we integrate genome-wide data on TF-binding sites, chromatin markers and functional annotations to predict genes associated with human eQTLs. Using the random forest classifier, we found that genomic proximity plus five TF and chromatin features are able to predict >90% of target genes within 1 megabase of eQTLs. Despite being regularly used to map target genes, proximity is not a good indicator of eQTL targets for genes 150 kilobases away, but insulators, TF co-occurrence, open chromatin and functional similarities between TFs and genes are better indicators. Using all six features in the classifier achieved an area under the specificity and sensitivity curve of 0.91, much better compared with at most 0.75 for using any single feature. We hope this study will not only provide validation of eQTL-mapping studies, but also provide insight into the molecular mechanisms explaining how genetic variation can influence gene expression.

Myogenin recruits the histone chaperone facilitates chromatin transcription (FACT) to promote nucleosome disassembly at muscle-specific genes

Facilitates chromatin transcription (FACT) functions to reorganize nucleosomes by acting as a histone chaperone that destabilizes and restores nucleosomal structure. The FACT complex is composed of two subunits: SSRP1 and SPT16. We have discovered that myogenin interacts with the FACT complex. Transfection of FACT subunits with myogenin is highly stimulatory for endogenous muscle gene expression in 10T1/2 cells. We have also found that FACT subunits do not associate with differentiation-specific genes while C2C12 cells are proliferating but are recruited to muscle-specific genes as differentiation initiates and then dissociate as differentiation proceeds. The recruitment is dependent on myogenin, as knockdowns of myogenin show no recruitment of the FACT complex. These data suggest that FACT is involved in the early steps of gene activation through its histone chaperone activities that serve to open the chromatin structure and facilitate transcription. Consistent with this hypothesis, we find that nucleosomes are depleted at muscle-specific promoters upon differentiation and that this activity is dependent on the presence of FACT. Our results show that the FACT complex promotes myogenin-dependent transcription and suggest that FACT plays an important role in the establishment of the appropriate transcription profile in a differentiated muscle cell.

Epigenetic modulation by TFII-I during embryonic stem cell differentiation

TFII-I transcription factors play an essential role during early vertebrate embryogenesis. Genome-wide mapping studies by ChIP-seq and ChIP-chip revealed that TFII-I primes multiple genomic loci in mouse embryonic stem cells and embryonic tissues. Moreover, many TFII-I-bound regions co-localize with H3K4me3/K27me3 bivalent chromatin within the promoters of lineage-specific genes. This minireview provides a summary of current knowledge regarding the function of TFII-I in epigenetic control of stem cell differentiation.

DNA methylation and differentiation: silencing, upregulation and modulation of gene expression

Differentiation-related DNA methylation is receiving increasing attention, partly owing to new, whole-genome analyses. These revealed that cell type-specific differential methylation in gene bodies is more frequent than in promoters. We review new insights into the functionality of DNA methylation during differentiation, with emphasis on the methylomes of myoblasts, myotubes and skeletal muscle versus non-muscle samples. Biostatistical analyses of data from reduced representation bisulfite sequencing are discussed. Lastly, a model is presented for how promoter and intragenic DNA hypermethylation affect gene expression, including increasing the efficiency of polycomb silencing at some promoters, downmodulating other promoters rather than silencing them, counteracting enhancers with heterologous specificity, altering chromatin conformation by inhibiting the binding of CTCF, modulating mRNA transcript levels by inhibiting overlapping promoters of noncoding RNA genes or by regulating the use of alternative mRNA promoters, modulating transcription termination, regulating alternative splicing and acting as barriers to the spread of activating chromatin.

Molecular roadblocks for cellular reprogramming

During development, diverse cellular identities are established and maintained in the embryo. Although remarkably robust in vivo, cellular identities can be manipulated using experimental techniques. Lineage reprogramming is an emerging field at the intersection of developmental and stem cell biology in which a somatic cell is stably reprogrammed into a distinct cell type by forced expression of lineage-determining factors. Lineage reprogramming enables the direct conversion of readily available cells from patients (such as skin fibroblasts) into disease-relevant cell types (such as neurons and cardiomyocytes) or into induced pluripotent stem cells. Although remarkable progress has been made in developing novel reprogramming methods, the efficiency and fidelity of reprogramming need to be improved in order increase the experimental and translational utility of reprogrammed cells. Studying the mechanisms that prevent successful reprogramming should allow for improvements in reprogramming methods, which could have significant implications for regenerative medicine and the study of human disease. Furthermore, lineage reprogramming has the potential to become a powerful system for dissecting the mechanisms that underlie cell fate establishment and terminal differentiation processes. In this review, we will discuss how transcription factors interface with the genome and induce changes in cellular identity in the context of development and reprogramming.

Comparison of genome-wide binding of MyoD in normal human myogenic cells and rhabdomyosarcomas identifies regional and local suppression of promyogenic transcription factors

Rhabdomyosarcoma is a pediatric tumor of skeletal muscle that expresses the myogenic basic helix-loop-helix protein MyoD but fails to undergo terminal differentiation. Prior work has determined that DNA binding by MyoD occurs in the tumor cells, but myogenic targets fail to activate. Using MyoD chromatin immunoprecipitation coupled to high-throughput sequencing and gene expression analysis in both primary human muscle cells and RD rhabdomyosarcoma cells, we demonstrate that MyoD binds in a similar genome-wide pattern in both tumor and normal cells but binds poorly at a subset of myogenic genes that fail to activate in the tumor cells. Binding differences are found both across genomic regions and locally at specific sites that are associated with binding motifs for RUNX1, MEF2C, JDP2, and NFIC. These factors are expressed at lower levels in RD cells than muscle cells and rescue myogenesis when expressed in RD cells. MEF2C is located in a genomic region that exhibits poor MyoD binding in RD cells, whereas JDP2 exhibits local DNA hypermethylation in its promoter in both RD cells and primary tumor samples. These results demonstrate that regional and local silencing of differentiation factors contributes to the differentiation defect in rhabdomyosarcomas.

DNA regions bound at low occupancy by transcription factors do not drive patterned reporter gene expression in Drosophila

In animals, each sequence-specific transcription factor typically binds to thousands of genomic regions in vivo. Our previous studies of 20 transcription factors show that most genomic regions bound at high levels in Drosophila blastoderm embryos are known or probable functional targets, but genomic regions occupied only at low levels have characteristics suggesting that most are not involved in the cis-regulation of transcription. Here we use transgenic reporter gene assays to directly test the transcriptional activity of 104 genomic regions bound at different levels by the 20 transcription factors. Fifteen genomic regions were selected based solely on the DNA occupancy level of the transcription factor Kruppel. Five of the six most highly bound regions drive blastoderm patterns of reporter transcription. In contrast, only one of the nine lowly bound regions drives transcription at this stage and four of them are not detectably active at any stage of embryogenesis. A larger set of 89 genomic regions chosen using criteria designed to identify functional cis-regulatory regions supports the same trend: genomic regions occupied at high levels by transcription factors in vivo drive patterned gene expression, whereas those occupied only at lower levels mostly do not. These results support studies that indicate that the high cellular concentrations of sequence-specific transcription factors drive extensive, low-occupancy, nonfunctional interactions within the accessible portions of the genome.

BAF60 A, B, and Cs of muscle determination and renewal

Developmental biologists have defined many of the diffusible and transcription factors that control muscle differentiation, yet we still have only rudimentary knowledge of the mechanisms that dictate whether a myogenic progenitor cell forms muscle versus alternate lineages, including those that can be pathological in a state of disease or degeneration. Clues about the molecular basis for lineage determination in muscle progenitors are only now emerging from studies of chromatin modifications that avail myogenic genes for transcription, together with analysis of the composition and activities of the chromatin-modifying complexes themselves. Here we review recent progress on muscle determination and explore a unifying theme that environmental cues from the stem or progenitor niche control the selection of specific subunit variants of the switch/sucrose nonfermentable (SWI/SNF) chromatin-modifying complex, creating a combinatorial code that dictates whether cells adopt myogenic versus nonmyogenic cell fates. A key component of the code appears to be the mutually exclusive usage of the a, b, and c variants of the 60-kD structural subunit BAF60 (BRG1/BRM-associated factor 60), of which BAF60c is essential to activate both skeletal and cardiac muscle programs. Since chromatin remodeling governs myogenic fate, the combinatorial assembly of the SWI/SNF complex might be targeted to develop drugs aimed at the therapeutic reduction of compensatory fibrosis and fatty deposition in chronic muscular disorders.

A co-localization model of paired ChIP-seq data using a large ENCODE data set enables comparison of multiple samples

Deep sequencing approaches, such as chromatin immunoprecipitation by sequencing (ChIP-seq), have been successful in detecting transcription factor-binding sites and histone modification in the whole genome. An approach for comparing two different ChIP-seq data would be beneficial for predicting unknown functions of a factor. We propose a model to represent co-localization of two different ChIP-seq data. We showed that a meaningful overlapping signal and a meaningless background signal can be separated by this model. We applied this model to compare ChIP-seq data of RNA polymerase II C-terminal domain (CTD) serine 2 phosphorylation with a large amount of peak-called data, including ChIP-seq and other deep sequencing data in the Encyclopedia of DNA Elements (ENCODE) project, and then extracted factors that were related to RNA polymerase II CTD serine 2 in HeLa cells. We further analyzed RNA polymerase II CTD serine 7 phosphorylation, of which their function is still unclear in HeLa cells. Our results were characterized by the similarity of localization for transcription factor/histone modification in the ENCODE data set, and this suggests that our model is appropriate for understanding ChIP-seq data for factors where their function is unknown.

Notch3 and Mef2c proteins are mutually antagonistic via Mkp1 protein and miR-1/206 microRNAs in differentiating myoblasts

BACKGROUND

Notch3 is expressed in myogenic precursors, but its function is not well known.

RESULTS

Notch3 represses the activity of Mef2c and is in turn inhibited by the microRNAs-1 and -206.

CONCLUSION

Notch3 serves as a regulator for preventing premature myogenic differentiation.

SIGNIFICANCE

Understanding how precocious differentiation is prevented is critical for designing therapy for skeletal muscle regeneration. The Notch signaling pathway is a well known regulator of skeletal muscle stem cells known as satellite cells. Loss of Notch1 signaling leads to spontaneous myogenic differentiation. Notch1, normally expressed in satellite cells, is targeted for proteasomal degradation by Numb during differentiation. A homolog of Notch1, Notch3, is also expressed in these cells but is not inhibited by Numb. We find that Notch3 is paradoxically up-regulated during the early stages of differentiation by an enhancer that requires both MyoD and activated Notch1. Notch3 itself strongly inhibits the myogenic transcription factor Mef2c, most likely by increasing the p38 phosphatase Mkp1, which inhibits the Mef2c activator p38 MAP kinase. Active Notch3 decreases differentiation. Mef2c, however, induces microRNAs miR-1 and miR-206, which directly down-regulate Notch3 and allow differentiation to proceed. Thus, the myogenic differentiation-induced microRNAs miR-1 and -206 are important for differentiation at least partly because they turn off Notch3. We suggest that the transient expression of Notch3 early in differentiation generates a temporal lag between myoblast activation by MyoD and terminal differentiation into myotubes directed by Mef2c.

EZH2: a pivotal regulator in controlling cell differentiation

Epigenetic regulation plays an important role in stem cell self-renewal, maintenance and lineage differentiation. The epigenetic profiles of stem cells are related to their transcriptional signature. Enhancer of Zeste homlog 2 (EZH2), a catalytic subunit of epigenetic regulator Polycomb repressive complex 2 (PRC2), has been shown to be a key regulator in controlling cellular differentiation. EZH2 is a histone methyltransferase that not only methylates histone H3 on Lys 27 (H3K27me3) but also interacts with and recruits DNA methyltransferases to methylate CpG at certain EZH2 target genes to establish firm repressive chromatin structures, contributing to tumor progression and the regulation of development and lineage commitment both in embryonic stem cells (ESCs) and adult stem cells. In addition to its well-recognized epigenetic gene silencing function, EZH2 also directly methylates nonhistone targets such as the cardiac transcription factor, GATA4, resulting in attenuated GATA4 transcriptional activity and gene repression. This review addresses recent progress toward the understanding of the biological functions and regulatory mechanisms of EZH2 and its targets as well as their roles in stem cell maintenance and cell differentiation.