Genetic and epigenetic determinants of neurogenesis and myogenesis

Multiplexed spatial mapping of chromatin features, transcriptome and proteins in tissues

The phenotypic and functional states of cells are modulated by a complex interactive molecular hierarchy of multiple omics layers, involving the genome, epigenome, transcriptome, proteome and metabolome. Spatial omics approaches have enabled the study of these layers in tissue context but are often limited to one or two modalities, offering an incomplete view of cellular identity. Here we present spatial-Mux-seq, a multimodal spatial technology that allows simultaneous profiling of five different modalities: two histone modifications, chromatin accessibility, whole transcriptome and a panel of proteins at tissue scale and cellular level in a spatially resolved manner. We applied this technology to mouse embryos and mouse brains, generating detailed multimodal tissue maps that identified more cell types and states compared to unimodal data. This analysis uncovered spatiotemporal relationships among histone modifications, chromatin accessibility, gene expression and protein levels during neuron differentiation, and revealed a radial glia niche with spatially dynamic epigenetic signals. Collectively, the spatial multi-omics approach heralds a new era for characterizing tissue and cellular heterogeneity that single-modality studies alone could not reveal.

E-box independent chromatin recruitment turns MYOD into a transcriptional repressor

MYOD is an E-box sequence-specific basic Helix-Loop-Helix (bHLH) transcriptional activator that, when expressed in non-muscle cells, induces nuclear reprogramming toward skeletal myogenesis by promoting chromatin accessibility at previously silent loci. Here, we report on the identification of a previously unrecognized property of MYOD as repressor of gene expression, via E-box-independent chromatin binding within accessible genomic elements, which invariably leads to reduced chromatin accessibility. MYOD-mediated repression requires the integrity of functional domains previously implicated in MYOD-mediated activation of gene expression. Repression of mitogen- and growth factor-responsive genes occurs through promoter binding and requires a highly conserved domain within the first helix. Repression of cell-of-origin/alternative lineage genes occurs via binding and decommissioning of distal regulatory elements, such as super-enhancers (SE), which requires the N-terminal activation domain as well as two chromatin-remodeling domains and leads to reduced strength of CTCF-mediated chromatin interactions. Surprisingly, MYOD-mediated chromatin compaction and repression of transcription do not associate with reduction of H3K27ac, the conventional histone mark of enhancer or promoter activation, but with reduced levels of the recently discovered histone H4 acetyl-methyl lysine modification (Kacme). These results extend MYOD biological properties beyond the current dogma that restricts MYOD function to a monotone transcriptional activator and reveal a previously unrecognized functional versatility arising from an alternative chromatin recruitment through E-box or non-E-box sequences. The E-box independent repression of gene expression by MYOD might provide a promiscuous mechanism to reduce chromatin accessibility and repress cell-of-origin/alternative lineage and growth factor/mitogen-responsive genes to safeguard the integrity of cell identity during muscle progenitor commitment toward the myogenic lineage.

Recruitment of homodimeric proneural factors by conserved CAT-CAT E-boxes drives major epigenetic reconfiguration in cortical neurogenesis

Proneural factors of the basic helix-loop-helix family coordinate neurogenesis and neurodifferentiation. Among them, NEUROG2 and NEUROD2 subsequently act to specify neurons of the glutamatergic lineage. Disruption of these factors, their target genes and binding DNA motifs has been linked to various neuropsychiatric disorders. Proneural factors bind to specific DNA motifs called E-boxes (hexanucleotides of the form CANNTG, composed of two CAN half sites on opposed strands). While corticogenesis heavily relies on E-box activity, the collaboration of proneural factors on different E-box types and their chromatin remodeling mechanisms remain largely unknown. Here, we conducted a comprehensive analysis using chromatin immunoprecipitation followed by sequencing (ChIP-seq) data for NEUROG2 and NEUROD2, along with time-matched single-cell RNA-seq, ATAC-seq and DNA methylation data from the developing mouse cortex. Our findings show that these factors are highly enriched in transiently active genomic regions during intermediate stages of neuronal differentiation. Although they primarily bind CAG-containing E-boxes, their binding in dynamic regions is notably enriched in CAT-CAT E-boxes (i.e. CATATG, denoted as 5'3' half sites for dimers), which undergo significant DNA demethylation and exhibit the highest levels of evolutionary constraint. Aided by HT-SELEX data reanalysis, structural modeling and DNA footprinting, we propose that these proneural factors exert maximal chromatin remodeling influence during intermediate stages of neurogenesis by binding as homodimers to CAT-CAT motifs. This study provides an in-depth integrative analysis of the dynamic regulation of E-boxes during neuronal development, enhancing our understanding of the mechanisms underlying the binding specificity of critical proneural factors.

Multiplexed spatial mapping of chromatin features, transcriptome, and proteins in tissues

The phenotypic and functional states of a cell are modulated by a complex interactive molecular hierarchy of multiple omics layers, involving the genome, epigenome, transcriptome, proteome, and metabolome. Spatial omics approaches have enabled the capture of information from different molecular layers directly in the tissue context. However, current technologies are limited to map one to two modalities at the same time, providing an incomplete representation of cellular identity. Such data is inadequate to fully understand complex biological systems and their underlying regulatory mechanisms. Here we present spatial-Mux-seq, a multi-modal spatial technology that allows simultaneous profiling of five different modalities, including genome-wide profiles of two histone modifications and open chromatin, whole transcriptome, and a panel of proteins at tissue scale and cellular level in a spatially resolved manner. We applied this technology to generate multi-modal tissue maps in mouse embryos and mouse brains, which discriminated more cell types and states than unimodal data. We investigated the spatiotemporal relationship between histone modifications, chromatin accessibility, gene and protein expression in neuron differentiation revealing the relationship between tissue organization, function, and gene regulatory networks. We were able to identify a radial glia spatial niche and revealed spatially changing gradient of epigenetic signals in this region. Moreover, we revealed previously unappreciated involvement of repressive histone marks in the mouse hippocampus. Collectively, the spatial multi-omics approach heralds a new era for characterizing tissue and cellular heterogeneity that single modality studies alone could not reveal.

SEAMoD: A fully interpretable neural network for cis-regulatory analysis of differentially expressed genes

A common way to investigate gene regulatory mechanisms is to identify differentially expressed genes using transcriptomics, find their candidate enhancers using epigenomics, and search for over-represented transcription factor (TF) motifs in these enhancers using bioinformatics tools. A related follow-up task is to model gene expression as a function of enhancer sequences and rank TF motifs by their contribution to such models, thus prioritizing among regulators. We present a new computational tool called SEAMoD that performs the above tasks of motif finding and sequence-to-expression modeling simultaneously. It trains a convolutional neural network model to relate enhancer sequences to differential expression in one or more biological conditions. The model uses TF motifs to interpret the sequences, learning these motifs and their relative importance to each biological condition from data. It also utilizes epigenomic information in the form of activity scores of putative enhancers and automatically searches for the most promising enhancer for each gene. Compared to existing neural network models of non-coding sequences, SEAMoD uses far fewer parameters, requires far less training data, and emphasizes biological interpretability. We used SEAMoD to understand regulatory mechanisms underlying the differentiation of neural stem cell (NSC) derived from mouse forebrain. We profiled gene expression and histone modifications in NSC and three differentiated cell types and used SEAMoD to model differential expression of nearly 12,000 genes with an accuracy of 81%, in the process identifying the Olig2, E2f family TFs, Foxo3, and Tcf4 as key transcriptional regulators of the differentiation process.

MyoD-Induced Trans-Differentiation: A Paradigm for Dissecting the Molecular Mechanisms of Cell Commitment, Differentiation and Reprogramming

The discovery of the skeletal muscle-specific transcription factor MyoD represents a milestone in the field of transcriptional regulation during differentiation and cell-fate reprogramming. MyoD was the first tissue-specific factor found capable of converting non-muscle somatic cells into skeletal muscle cells. A unique feature of MyoD, with respect to other lineage-specific factors able to drive trans-differentiation processes, is its ability to dramatically change the cell fate even when expressed alone. The present review will outline the molecular strategies by which MyoD reprograms the transcriptional regulation of the cell of origin during the myogenic conversion, focusing on the activation and coordination of a complex network of co-factors and epigenetic mechanisms. Some molecular roadblocks, found to restrain MyoD-dependent trans-differentiation, and the possible ways for overcoming these barriers, will also be discussed. Indeed, they are of critical importance not only to expand our knowledge of basic muscle biology but also to improve the generation skeletal muscle cells for translational research.

Epigenetic genes and epilepsy - emerging mechanisms and clinical applications

An increasing number of epilepsies are being attributed to variants in genes with epigenetic functions. The products of these genes include factors that regulate the structure and function of chromatin and the placing, reading and removal of epigenetic marks, as well as other epigenetic processes. In this Review, we provide an overview of the various epigenetic processes, structuring our discussion around five function-based categories: DNA methylation, histone modifications, histone-DNA crosstalk, non-coding RNAs and chromatin remodelling. We provide background information on each category, describing the general mechanism by which each process leads to altered gene expression. We also highlight key clinical and mechanistic aspects, providing examples of genes that strongly associate with epilepsy within each class. We consider the practical applications of these findings, including tissue-based and biofluid-based diagnostics and precision medicine-based treatments. We conclude that variants in epigenetic genes are increasingly found to be causally involved in the epilepsies, with implications for disease mechanisms, treatments and diagnostics.

Disturbance of calcium homeostasis and myogenesis caused by TET2 deletion in muscle stem cells

Skeletal muscle myogenesis is a sophisticated process controlled by genetic and epigenetic regulators. In animals, one of the key enzymes for the DNA demethylation of 5-methylcytosine is TET2. Although TET2 is essential for muscle development, the mechanisms by which TET2 regulates myogenesis, particularly the implication for muscle stem cells, remains unclear. In the present study, we employed the TET2 knockout mouse model to investigate the function of TET2 in muscle development and regeneration. We observed that TET2 deficiency caused impaired muscle stem cell proliferation and differentiation, resulting in the reduction in both myofiber number and muscle tissue size. Specifically, TET2 maintains calcium homeostasis in muscle stem cells by controlling the DNA methylation levels of the calcium pathway genes. Forced expression of the sodium/calcium exchanger protein SLC8A3 could rescue the myogenic defects in TET2 knockout cells. Our data not only illustrated the vital function of TET2 during myogenesis but also identified novel targets that contribute to calcium homeostasis for enhancing muscle function.

MyoD is a 3D genome structure organizer for muscle cell identity

The genome exists as an organized, three-dimensional (3D) dynamic architecture, and each cell type has a unique 3D genome organization that determines its cell identity. An unresolved question is how cell type-specific 3D genome structures are established during development. Here, we analyzed 3D genome structures in muscle cells from mice lacking the muscle lineage transcription factor (TF), MyoD, versus wild-type mice. We show that MyoD functions as a “genome organizer” that specifies 3D genome architecture unique to muscle cell development, and that H3K27ac is insufficient for the establishment of MyoD-induced chromatin loops in muscle cells. Moreover, we present evidence that other cell lineage-specific TFs might also exert functional roles in orchestrating lineage-specific 3D genome organization during development.

Pioneer transcription factors (TFs) have been proposed to act as protein anchors to orchestrate cell type-specific 3D genome architecture. MyoD is a pioneer TF for myogenic lineage specification. Here the authors provide further support for the role of MyoD in 3D genome architecture in muscle stem cells by comparing MyoD knockout and wild-type mice.

Programmed genomic instability regulates neural transdifferentiation of human brain microvascular pericytes

BACKGROUND

Transdifferentiation describes transformation in vivo of specialized cells from one lineage into another. While there is extensive literature on forced induction of lineage reprogramming in vitro, endogenous mechanisms that govern transdifferentiation remain largely unknown. The observation that human microvascular pericytes transdifferentiate into neurons provided an opportunity to explore the endogenous molecular basis for lineage reprogramming.

RESULTS

We show that abrupt destabilization of the higher-order chromatin topology that chaperones lineage memory of pericytes is driven by transient global transcriptional arrest. This leads within minutes to localized decompression of the repressed competing higher-order chromatin topology and expression of pro-neural genes. Transition to neural lineage is completed by probabilistic induction of R-loops in key myogenic loci upon re-initiation of RNA polymerase activity, leading to depletion of the myogenic transcriptome and emergence of the neurogenic transcriptome.

CONCLUSIONS

These findings suggest that the global transcriptional landscape not only shapes the functional cellular identity of pericytes, but also stabilizes lineage memory by silencing the competing neural program within a repressed chromatin state.

Mechanisms of Binding Specificity among bHLH Transcription Factors

The transcriptome of every cell is orchestrated by the complex network of interaction between transcription factors (TFs) and their binding sites on DNA. Disruption of this network can result in many forms of organism malfunction but also can be the substrate of positive natural selection. However, understanding the specific determinants of each of these individual TF-DNA interactions is a challenging task as it requires integrating the multiple possible mechanisms by which a given TF ends up interacting with a specific genomic region. These mechanisms include DNA motif preferences, which can be determined by nucleotide sequence but also by DNA’s shape; post-translational modifications of the TF, such as phosphorylation; and dimerization partners and co-factors, which can mediate multiple forms of direct or indirect cooperative binding. Binding can also be affected by epigenetic modifications of putative target regions, including DNA methylation and nucleosome occupancy. In this review, we describe how all these mechanisms have a role and crosstalk in one specific family of TFs, the basic helix-loop-helix (bHLH), with a very conserved DNA binding domain and a similar DNA preferred motif, the E-box. Here, we compile and discuss a rich catalog of strategies used by bHLH to acquire TF-specific genome-wide landscapes of binding sites.

Myod1 and GR coordinate myofiber-specific transcriptional enhancers

Skeletal muscle is a dynamic tissue the size of which can be remodeled through the concerted actions of various cues. Here, we investigated the skeletal muscle transcriptional program and identified key tissue-specific regulatory genetic elements. Our results show that Myod1 is bound to numerous skeletal muscle enhancers in collaboration with the glucocorticoid receptor (GR) to control gene expression. Remarkably, transcriptional activation controlled by these factors occurs through direct contacts with the promoter region of target genes, via the CpG-bound transcription factor Nrf1, and the formation of Ctcf-anchored chromatin loops, in a myofiber-specific manner. Moreover, we demonstrate that GR negatively controls muscle mass and strength in mice by down-regulating anabolic pathways. Taken together, our data establish Myod1, GR and Nrf1 as key players of muscle-specific enhancer-promoter communication that orchestrate myofiber size regulation.

Fusion of single-cell transcriptome and DNA-binding data, for genomic network inference in cortical development

BACKGROUND

Network models are well-established as very useful computational-statistical tools in cell biology. However, a genomic network model based only on gene expression data can, by definition, only infer gene co-expression networks. Hence, in order to infer gene regulatory patterns, it is necessary to also include data related to binding of regulatory factors to DNA.

RESULTS

We propose a new dynamic genomic network model, for inferring patterns of genomic regulatory influence in dynamic processes such as development. Our model fuses experiment-specific gene expression data with publicly available DNA-binding data. The method we propose is computationally efficient, and can be applied to genome-wide data with tens of thousands of transcripts. Thus, our method is well suited for use as an exploratory tool for genome-wide data. We apply our method to data from human fetal cortical development, and our findings confirm genomic regulatory patterns which are recognised as being fundamental to neuronal development.

CONCLUSIONS

Our method provides a mathematical/computational toolbox which, when coupled with targeted experiments, will reveal and confirm important new functional genomic regulatory processes in mammalian development.

Intronic enhancer region governs transcript-specific Bdnf expression in rodent neurons

Brain-derived neurotrophic factor (BDNF) controls the survival, growth, and function of neurons both during the development and in the adult nervous system. Bdnf is transcribed from several distinct promoters generating transcripts with alternative 5' exons. Bdnf transcripts initiated at the first cluster of exons have been associated with the regulation of body weight and various aspects of social behavior, but the mechanisms driving the expression of these transcripts have remained poorly understood. Here, we identify an evolutionarily conserved intronic enhancer region inside the Bdnf gene that regulates both basal and stimulus-dependent expression of the Bdnf transcripts starting from the first cluster of 5' exons in mouse and rat neurons. We further uncover a functional E-box element in the enhancer region, linking the expression of Bdnf and various pro-neural basic helix–loop–helix transcription factors. Collectively, our results shed new light on the cell-type- and stimulus-specific regulation of the important neurotrophic factor BDNF.

Disruption of NEUROD2 causes a neurodevelopmental syndrome with autistic features via cell-autonomous defects in forebrain glutamatergic neurons

While the transcription factor NEUROD2 has recently been associated with epilepsy, its precise role during nervous system development remains unclear. Using a multi-scale approach, we set out to understand how Neurod2 deletion affects the development of the cerebral cortex in mice. In Neurod2 KO embryos, cortical projection neurons over-migrated, thereby altering the final size and position of layers. In juvenile and adults, spine density and turnover were dysregulated in apical but not basal compartments in layer 5 neurons. Patch-clamp recordings in layer 5 neurons of juvenile mice revealed increased intrinsic excitability. Bulk RNA sequencing showed dysregulated expression of many genes associated with neuronal excitability and synaptic function, whose human orthologs were strongly associated with autism spectrum disorders (ASD). At the behavior level, Neurod2 KO mice displayed social interaction deficits, stereotypies, hyperactivity, and occasionally spontaneous seizures. Mice heterozygous for Neurod2 had similar defects, indicating that Neurod2 is haploinsufficient. Finally, specific deletion of Neurod2 in forebrain excitatory neurons recapitulated cellular and behavioral phenotypes found in constitutive KO mice, revealing the region-specific contribution of dysfunctional Neurod2 in symptoms. Informed by these neurobehavioral features in mouse mutants, we identified eleven patients from eight families with a neurodevelopmental disorder including intellectual disability and ASD associated with NEUROD2 pathogenic mutations. Our findings demonstrate crucial roles for Neurod2 in neocortical development, whose alterations can cause neurodevelopmental disorders including intellectual disability and ASD.

Target Molecules of STIM Proteins in the Central Nervous System

Stromal interaction molecules (STIMs), including STIM1 and STIM2, are single-pass transmembrane proteins that are located predominantly in the endoplasmic reticulum (ER). They serve as calcium ion (Ca²⁺) sensors within the ER. In the central nervous system (CNS), they are involved mainly in Orai-mediated store-operated Ca²⁺ entry (SOCE). The key molecular components of the SOCE pathway are well-characterized, but the molecular mechanisms that underlie the regulation of this pathway need further investigation. Numerous intracellular target proteins that are located in the plasma membrane, ER, cytoskeleton, and cytoplasm have been reported to play essential roles in concert with STIMs, such as conformational changes in STIMs, their translocation, the stabilization of their interactions with Orai, and the activation of other channels. The present review focuses on numerous regulators, such as Homer, SOCE-associated regulatory factor (SARAF), septin, synaptopodin, golli proteins, partner of STIM1 (POST), and transcription factors and proteasome inhibitors that regulate STIM-Orai interactions in the CNS. Further we describe novel roles of STIMs in mediating Ca²⁺ influx via other than Orai pathways, including TRPC channels, VGCCs, AMPA and NMDA receptors, and group I metabotropic glutamate receptors. This review also summarizes recent findings on additional molecular targets of STIM proteins including SERCA, IP₃Rs, end-binding proteins (EB), presenilin, and CaMKII. Dysregulation of the SOCE-associated toolkit, including STIMs, contributes to the development of neurodegenerative disorders (e.g., Alzheimer's disease, Parkinson's disease, and Huntington's disease), traumatic brain injury, epilepsy, and stroke. Emerging evidence points to the role of STIM proteins and several of their molecular effectors and regulators in neuronal and glial physiology and pathology, suggesting their potential application for future therapeutic strategies.

Adaptation of the carbamoyl-phosphate synthetase enzyme in an extremophile fish

Tetrapods and fish have adapted distinct carbamoyl-phosphate synthase (CPS) enzymes to initiate the ornithine urea cycle during the detoxification of nitrogenous wastes. We report evidence that in the ureotelic subgenus of extremophile fish Oreochromis Alcolapia, CPS III has undergone convergent evolution and adapted its substrate affinity to ammonia, which is typical of terrestrial vertebrate CPS I. Unusually, unlike in other vertebrates, the expression of CPS III in Alcolapia is localized to the skeletal muscle and is activated in the myogenic lineage during early embryonic development with expression remaining in mature fish. We propose that adaptation in Alcolapia included both convergent evolution of CPS function to that of terrestrial vertebrates, as well as changes in development mechanisms redirecting CPS III gene expression to the skeletal muscle.

Poly(ADP-ribose) Polymerase 1 (PARP1) restrains MyoD-dependent gene expression during muscle differentiation

The myogenic factor MyoD regulates skeletal muscle differentiation by interacting with a variety of chromatin-modifying complexes. Although MyoD can induce and maintain chromatin accessibility at its target genes, its binding and trans-activation ability can be limited by some types of not fully characterized epigenetic constraints. In this work we analysed the role of PARP1 in regulating MyoD-dependent gene expression. PARP1 is a chromatin-associated enzyme, playing a well recognized role in DNA repair and that is implicated in transcriptional regulation. PARP1 affects gene expression through multiple mechanisms, often involving the Poly(ADP-ribosyl)ation of chromatin proteins. In line with PARP1 down-regulation during differentiation, we observed that PARP1 depletion boosts the up-regulation of MyoD targets, such as p57, myogenin, Mef2C and p21, while its re-expression reverts this effect. We also found that PARP1 interacts with some MyoD-binding regions and that its presence, independently of the enzymatic activity, interferes with MyoD recruitment and gene induction. We finally suggest a relationship between the binding of PARP1 and the loss of the activating histone modification H3K4me3 at MyoD-binding regions. This work highlights not only a novel player in the epigenetic control of myogenesis, but also a repressive and catalytic-independent mechanisms by which PARP1 regulates transcription.

Functional analysis of two MyoDs revealed their role in the activation of myomixer expression in yellowfin seabream (Acanthopagrus latus) (Hottuyn, 1782)

Myoblast determination protein (MyoD), a muscle-specific basic helix-loop-helix (bHLH) transcription factor, plays a pivotal role in regulating skeletal muscle growth and development. However, the regulation mechanism of MyoD has not been determined in marine fishes. In the present study, we isolated the MyoD1 (AlMyoD1) and MyoD2 (AlMyoD2) genomic sequences and analyzed the expression patterns in different tissues of yellowfin seabream (Acanthopagrus latus). The open reading frame (ORF) sequences of AlMyoD1 and AlMyoD2 encoded 297 and 271 amino acids possessing three common characteristic domains, respectively, containing a myogenic basic domain, a bHLH domain, and a ser-rich region (helix III). Phylogenetic and genome structure analyses exhibited classic phylogeny and highly conserved exon/intron architecture. Furthermore, the AlMyoD1 and AlMyoD2 transcription levels were higher in white muscle than in the other tissues. In order to further study AlMyoD function in muscle, promoter sequence analysis found that several E-box binding sites were present. Additionally, binding sites of Almyomixer involved in mammal myoblast fusion, which expression was also the highest in white muscle, were found in the promoter of AlMyoD. Pomoter activity assays further confirmed that both AlMyoD1 and AlMyoD2 can dramatically activate Almyomixer expression, and the AlMyoD1 M2 and AlMyoD2 M5 E-box binding sites were functionally important for Almyomixer transcription based on mutation analysis and electrophoretic mobile shift assays (EMSA). In summary, two MyoDs play a core role in Almyomixer regulation and may promote myofibre formation during muscle development and growth by regulating Almyomixer expression.

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system1. Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis-element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

Dissecting myogenin-mediated retinoid X receptor signaling in myogenic differentiation

Deciphering the molecular mechanisms underpinning myoblast differentiation is a critical step in developing the best strategy to promote muscle regeneration in patients suffering from muscle-related diseases. We have previously established that a rexinoid x receptor (RXR)-selective agonist, bexarotene, enhances the differentiation and fusion of myoblasts through a direct regulation of MyoD expression, coupled with an augmentation of myogenin protein. Here, we found that RXR signaling associates with the distribution of myogenin at poised enhancers and a distinct E-box motif. We also found an association of myogenin with rexinoid-responsive gene expression and identified an epigenetic signature related to histone acetyltransferase p300. Moreover, RXR signaling augments residue-specific histone acetylation at enhancers co-occupied by p300 and myogenin. Thus, genomic distribution of transcriptional regulators is an important designate for identifying novel targets as well as developing therapeutics that modulate epigenetic landscape in a selective manner to promote muscle regeneration.

Sequence and chromatin determinants of transcription factor binding and the establishment of cell type-specific binding patterns

Transcription factors (TFs) selectively bind distinct sets of sites in different cell types. Such cell type-specific binding specificity is expected to result from interplay between the TF's intrinsic sequence preferences, cooperative interactions with other regulatory proteins, and cell type-specific chromatin landscapes. Cell type-specific TF binding events are highly correlated with patterns of chromatin accessibility and active histone modifications in the same cell type. However, since concurrent chromatin may itself be a consequence of TF binding, chromatin landscapes measured prior to TF activation provide more useful insights into how cell type-specific TF binding events became established in the first place. Here, we review the various sequence and chromatin determinants of cell type-specific TF binding specificity. We identify the current challenges and opportunities associated with computational approaches to characterizing, imputing, and predicting cell type-specific TF binding patterns. We further focus on studies that characterize TF binding in dynamic regulatory settings, and we discuss how these studies are leading to a more complex and nuanced understanding of dynamic protein-DNA binding activities. We propose that TF binding activities at individual sites can be viewed along a two-dimensional continuum of local sequence and chromatin context. Under this view, cell type-specific TF binding activities may result from either strongly favorable sequence features or strongly favorable chromatin context.

The regulatory role of Myomaker and Myomixer-Myomerger-Minion in muscle development and regeneration

Skeletal muscle plays essential roles in motor function, energy, and glucose metabolism. Skeletal muscle formation occurs through a process called myogenesis, in which a crucial step is the fusion of mononucleated myoblasts to form multinucleated myofibers. The myoblast/myocyte fusion is triggered and coordinated in a muscle-specific way that is essential for muscle development and post-natal muscle regeneration. Many molecules and proteins have been found and demonstrated to have the capacity to regulate the fusion of myoblast/myocytes. Interestingly, two newly discovered muscle-specific membrane proteins, Myomaker and Myomixer (also called Myomerger and Minion), have been identified as fusogenic regulators in vertebrates. Both Myomaker and Myomixer-Myomerger-Minion have the capacity to directly control the myogenic fusion process. Here, we review and discuss the latest studies related to these two proteins, including the discovery, structure, expression pattern, functions, and regulation of Myomaker and Myomixer-Myomerger-Minion. We also emphasize and discuss the interaction between Myomaker and Myomixer-Myomerger-Minion, as well as their cooperative regulatory roles in cell-cell fusion. Moreover, we highlight the areas for exploration of Myomaker and Myomixer-Myomerger-Minion in future studies and consider their potential application to control cell fusion for cell-therapy purposes.

Pro-neuronal activity of Myod1 due to promiscuous binding to neuronal genes

The on-target pioneer factors Ascl1 and Myod1 are sequence-related but induce two developmentally unrelated lineages-that is, neuronal and muscle identities, respectively. It is unclear how these two basic helix-loop-helix (bHLH) factors mediate such fundamentally different outcomes. The chromatin binding of Ascl1 and Myod1 was surprisingly similar in fibroblasts, yet their transcriptional outputs were drastically different. We found that quantitative binding differences explained differential chromatin remodelling and gene activation. Although strong Ascl1 binding was exclusively associated with bHLH motifs, strong Myod1-binding sites were co-enriched with non-bHLH motifs, possibly explaining why Ascl1 is less context dependent. Finally, we observed that promiscuous binding of Myod1 to neuronal targets results in neuronal reprogramming when the muscle program is inhibited by Myt1l. Our findings suggest that chromatin access of on-target pioneer factors is primarily driven by the protein-DNA interaction, unlike ordinary context-dependent transcription factors, and that promiscuous transcription factor binding requires specific silencing mechanisms to ensure lineage fidelity.

Reprogramming of DNA methylation at NEUROD2-bound sequences during cortical neuron differentiation

The characteristics of DNA methylation changes that occur during neurogenesis in vivo remain unknown. We used whole-genome bisulfite sequencing to quantitate DNA cytosine modifications in differentiating neurons and their progenitors isolated from mouse brain at the peak of embryonic neurogenesis. Localized DNA hypomethylation was much more common than hypermethylation and often occurred at putative enhancers within genes that were upregulated in neurons and encoded proteins crucial for neuronal differentiation. The hypomethylated regions strongly overlapped with mapped binding sites of the key neuronal transcription factor NEUROD2. The 5-methylcytosine oxidase ten-eleven translocation 2 (TET2) interacted with NEUROD2, and its reaction product 5-hydroxymethylcytosine accumulated at the demethylated regions. NEUROD2-targeted differentially methylated regions retained higher methylation levels in Neurod2 knockout mice, and inducible expression of NEUROD2 caused TET2-associated demethylation at its in vivo binding sites. The data suggest that the reorganization of DNA methylation in developing neurons involves NEUROD2 and TET2-mediated DNA demethylation.

Transcription Factor-Directed Re-wiring of Chromatin Architecture for Somatic Cell Nuclear Reprogramming toward trans-Differentiation

MYOD-directed fibroblast trans-differentiation into skeletal muscle provides a unique model to investigate how one transcription factor (TF) reconfigures the three-dimensional chromatin architecture to control gene expression, which is otherwise achieved by the combinatorial activities of multiple TFs. Integrative analysis of genome-wide high-resolution chromatin interactions, MYOD and CTCF DNA-binding profile, and gene expression, revealed that MYOD directs extensive re-wiring of interactions involving cis-regulatory and structural genomic elements, including promoters, enhancers, and insulated neighborhoods (INs). Re-configured INs were hot-spots of differential interactions, whereby MYOD binding to highly constrained sequences at IN boundaries and/or inside INs led to alterations of promoter-enhancer interactions to repress cell-of-origin genes and to activate muscle-specific genes. Functional evidence shows that MYOD-directed re-configuration of chromatin interactions temporally preceded the effect on gene expression and was mediated by direct MYOD-DNA binding. These data illustrate a model whereby a single TF alters multi-loop hubs to drive somatic cell trans-differentiation.

Cell Reprogramming: The Many Roads to Success

Cellular reprogramming experiments from somatic cell types have demonstrated the plasticity of terminally differentiated cell states. Recent efforts in understanding the mechanisms of cellular reprogramming have begun to elucidate the differentiation trajectories along the reprogramming processes. In this review, we focus mainly on direct reprogramming strategies by transcription factors and highlight the variables that contribute to cell fate conversion outcomes. We review key studies that shed light on the cellular and molecular mechanisms by investigating differentiation trajectories and alternative cell states as well as transcription factor regulatory activities during cell fate reprogramming. Finally, we highlight a few concepts that we believe require attention, particularly when measuring the success of cell reprogramming experiments.

Cis-regulatory determinants of MyoD function

Muscle-specific transcription factor MyoD orchestrates the myogenic gene expression program by binding to short DNA motifs called E-boxes within myogenic cis-regulatory elements (CREs). Genome-wide analyses of MyoD cistrome by chromatin immnunoprecipitation sequencing shows that MyoD-bound CREs contain multiple E-boxes of various sequences. However, how E-box numbers, sequences and their spatial arrangement within CREs collectively regulate the binding affinity and transcriptional activity of MyoD remain largely unknown. Here, by an integrative analysis of MyoD cistrome combined with genome-wide analysis of key regulatory histones and gene expression data we show that the affinity landscape of MyoD is driven by multiple E-boxes, and that the overall binding affinity—and associated nucleosome positioning and epigenetic features of the CREs—crucially depend on the variant sequences and positioning of the E-boxes within the CREs. By comparative genomic analysis of single nucleotide polymorphism (SNPs) across publicly available data from 17 strains of laboratory mice, we show that variant sequences within the MyoD-bound motifs, but not their genome-wide counterparts, are under selection. At last, we show that the quantitative regulatory effect of MyoD binding on the nearby genes can, in part, be predicted by the motif composition of the CREs to which it binds. Taken together, our data suggest that motif numbers, sequences and their spatial arrangement within the myogenic CREs are important determinants of the cis-regulatory code of myogenic CREs.

Causal Gene Regulatory Network Modeling and Genomics: Second-Generation Challenges

Gene regulatory network modeling has played a major role in advancing the understanding of developmental systems, by crystallizing structures of relevant extant information, by formally posing hypothetical functional relationships between network elements, and by offering clear predictive tests to improve understanding of the mechanisms driving developmental progression. Both ordinary differential equation (ODE)-based and Boolean models have also been highly successful in explaining dynamics within subcircuits of more complex processes. In a very small number of cases, gene regulatory network models of much more global scope have been proposed that successfully predict the dynamics of the processes establishing most of an embryonic body plan. Can such successes be expanded to very different developmental systems, including post-embryonic mammalian systems? This perspective discusses several problems that must be solved in more quantitative and predictive theoretical terms, to make this possible. These problems include: the effects of cellular history on chromatin state and how these affect gene accessibility; the dose dependence of activities of many transcription factors (a problem for Boolean models); stochasticity of some transcriptional outputs (a problem for simple ODE models); response timing delays due to epigenetic remodeling requirements; functionally different kinds of repression; and the regulatory syntax that governs responses of genes with multiple enhancers.

Twist2 amplification in rhabdomyosarcoma represses myogenesis and promotes oncogenesis by redirecting MyoD DNA binding

Li et al. show that TWIST2 shapes the epigenetic landscape to drive chromatin opening at oncogenic loci and chromatin closing at myogenic loci. These epigenetic changes redirect MyoD binding from myogenic genes toward oncogenic, metabolic, and growth genes.

Rhabdomyosarcoma (RMS) is an aggressive pediatric cancer composed of myoblast-like cells. Recently, we discovered a unique muscle progenitor marked by the expression of the Twist2 transcription factor. Genomic analyses of 258 RMS patient tumors uncovered prevalent copy number amplification events and increased expression of TWIST2 in fusion-negative RMS. Knockdown of TWIST2 in RMS cells results in up-regulation of MYOGENIN and a decrease in proliferation, implicating TWIST2 as an oncogene in RMS. Through an inducible Twist2 expression system, we identified Twist2 as a reversible inhibitor of myogenic differentiation with the remarkable ability to promote myotube dedifferentiation in vitro. Integrated analysis of genome-wide ChIP-seq and RNA-seq data revealed the first dynamic chromatin and transcriptional landscape of Twist2 binding during myogenic differentiation. During differentiation, Twist2 competes with MyoD at shared DNA motifs to direct global gene transcription and repression of the myogenic program. Additionally, Twist2 shapes the epigenetic landscape to drive chromatin opening at oncogenic loci and chromatin closing at myogenic loci. These epigenetic changes redirect MyoD binding from myogenic genes toward oncogenic, metabolic, and growth genes. Our study reveals the dynamic interplay between two opposing transcriptional regulators that control the fate of RMS and provides insight into the molecular etiology of this aggressive form of cancer.

Mechanistic Insights in NeuroD Potentiation of Mineralocorticoid Receptor Signaling

Mineralocorticoid receptor (MR)-mediated signaling in the brain has been suggested as a protective factor in the development of psychopathology, in particular mood disorders. We recently identified genomic loci at which either MR or the closely related glucocorticoid receptor (GR) binds selectively, and found members of the NeuroD transcription factor family to be specifically associated with MR-bound DNA in the rat hippocampus. We show here using forebrain-specific MR knockout mice that GR binding to MR/GR joint target loci is not affected in any major way in the absence of MR. Neurod2 binding was also independent of MR binding. Moreover, functional comparison with MyoD family members indicates that it is the chromatin remodeling aspect of NeuroD, rather than its direct stimulation of transcription, that is responsible for potentiation of MR-mediated transcription. These findings suggest that NeuroD acts in a permissive way to enhance MR-mediated transcription, and they argue against competition for DNA binding as a mechanism of MR- over GR-specific binding.

Helix-loop-helix proteins and the advent of cellular diversity: 30 years of discovery

In this review from Murre, the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

Helix–loop–helix (HLH) proteins are dimeric transcription factors that control lineage- and developmental-specific gene programs. Genes encoding for HLH proteins arose in unicellular organisms >600 million years ago and then duplicated and diversified from ancestral genes across the metazoan and plant kingdoms to establish multicellularity. Hundreds of HLH proteins have been identified with diverse functions in a wide variety of cell types. HLH proteins orchestrate lineage specification, commitment, self-renewal, proliferation, differentiation, and homing. HLH proteins also regulate circadian clocks, protect against hypoxic stress, promote antigen receptor locus assembly, and program transdifferentiation. HLH proteins deposit or erase epigenetic marks, activate noncoding transcription, and sequester chromatin remodelers across the chromatin landscape to dictate enhancer–promoter communication and somatic recombination. Here the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

Loci-specific histone acetylation profiles associated with transcriptional coactivator p300 during early myoblast differentiation

Molecular regulation of stem cell differentiation is exerted through both genetic and epigenetic determinants over distal regulatory or enhancer regions. Understanding the mechanistic action of active or poised enhancers is therefore imperative for control of stem cell differentiation. Based on the genome-wide co-occurrence of different epigenetic marks in committed proliferating myoblasts, we have previously generated a 14-state chromatin state model to profile rexinoid-responsive histone acetylation in early myoblast differentiation. Here, we delineate the functional mode of transcription regulators during early myogenic differentiation using genome-wide chromatin state association. We define a role of transcriptional coactivator p300, when recruited by muscle master regulator MyoD, in the establishment and regulation of myogenic loci at the onset of myoblast differentiation. In addition, we reveal an enrichment of loci-specific histone acetylation at p300 associated active or poised enhancers, particularly when enlisted by MyoD. We provide novel molecular insights into the regulation of myogenic enhancers by p300 in concert with MyoD. Our studies present a valuable aptitude for driving condition-specific chromatin state or enhancers pharmacologically to treat muscle-related diseases and for the identification of additional myogenic targets and molecular interactions for therapeutic development.

Abbreviations: MRF: Muscle regulatory factor; HAT: Histone acetyltransferase; CBP: CREB-binding protein; ES: Embryonic stem; ATCC: American type culture collection; DM: Differentiation medium; DMEM: Dulbecco’s Modified Eagle Medium; GM: Growth medium; GO: Gene ontology; GREAT: Genomic regions enrichment of annotations tool; FPKM: Fragments per kilobase of transcript per million; GEO: Gene expression omnibus; MACS: Model-based analysis for ChIP-seq

The muscle regulatory transcription factor MyoD participates with p53 to directly increase the expression of the pro-apoptotic Bcl2 family member PUMA

The muscle regulatory transcription factor MyoD is a master regulator of skeletal myoblast differentiation. We have previously reported that MyoD is also necessary for the elevated expression of the pro-apoptotic Bcl2 family member PUMA, and the ensuing apoptosis, that occurs in a subset of myoblasts induced to differentiate. Herein, we report the identification of a functional MyoD binding site within the extended PUMA promoter. In silico analysis of the murine PUMA extended promoter revealed three potential MyoD binding sites within 2 kb of the transcription start site. Expression from a luciferase reporter construct containing this 2 kb fragment was enhanced by activation of MyoD in both myoblasts and fibroblasts and diminished by silencing of MyoD in myoblasts. Experiments utilizing truncated versions of this promoter region revealed that the potential binding site at position - 857 was necessary for expression. Chromatin immunoprecipitation (ChIP) analysis confirmed binding of MyoD to the DNA region encompassing position - 857. The increase in MyoD binding to the PUMA promoter as a consequence of culture in differentiation media (DM) was comparable to the increase in MyoD binding at the myogenin promoter and was diminished in myoblasts silenced for MyoD expression. Finally, ChIP analysis using an antibody specific for the transcription factor p53 demonstrated that, in myoblasts silenced for MyoD expression, p53 binding to the PUMA promoter was diminished in response to culture in DM. These data indicate that MyoD plays a direct role in regulating PUMA expression and reveal functional consequences of MyoD expression on p53 mediated transcription of PUMA.

Shaping Gene Expression by Landscaping Chromatin Architecture: Lessons from a Master

Since its discovery as a skeletal muscle-specific transcription factor able to reprogram somatic cells into differentiated myofibers, MyoD has provided an instructive model to understand how transcription factors regulate gene expression. Reciprocally, studies of other transcriptional regulators have provided testable hypotheses to further understand how MyoD activates transcription. Using MyoD as a reference, in this review, we discuss the similarities and differences in the regulatory mechanisms employed by tissue-specific transcription factors to access DNA and regulate gene expression by cooperatively shaping the chromatin landscape within the context of cellular differentiation.

Loss of myogenic potential and fusion capacity of muscle stem cells isolated from contractured muscle in children with cerebral palsy

Cerebral palsy (CP) is the most common cause of pediatric neurodevelopmental and physical disability in the United States. It is defined as a group of motor disorders caused by a nonprogressive perinatal insult to the brain. Although the brain lesion is nonprogressive, there is a progressive, lifelong impact on skeletal muscles, which are shorter, spastic, and may develop debilitating contractures. Satellite cells are resident muscle stem cells that are indispensable for postnatal growth and regeneration of skeletal muscles. Here we measured the myogenic potential of satellite cells isolated from contractured muscles in children with CP. When compared with typically developing (TD) children, satellite cell-derived myoblasts from CP differentiated more slowly (slope: 0.013 (SD 0.013) CP vs. 0.091 (SD 0.024) TD over 24 h, P < 0.001) and fused less (fusion index: 21.3 (SD 8.6) CP vs. 81.3 (SD 7.7) TD after 48 h, P < 0.001) after exposure to low-serum conditions that stimulated myotube formation. This impairment was associated with downregulation of several markers important for myoblast fusion and myotube formation, including DNA methylation-dependent inhibition of promyogenic integrin-β 1D (ITGB1D) protein expression levels (-50% at 42 h), and ~25% loss of integrin-mediated focal adhesion kinase phosphorylation. The cytidine analog 5-Azacytidine (5-AZA), a demethylating agent, restored ITGB1D levels and promoted myogenesis in CP cultures. Our data demonstrate that muscle contractures in CP are associated with loss of satellite cell myogenic potential that is dependent on DNA methylation patterns affecting expression of genetic programs associated with muscle stem cell differentiation and muscle fiber formation.

Intrinsic DNA binding properties demonstrated for lineage-specifying basic helix-loop-helix transcription factors

During development, transcription factors select distinct gene programs, providing the necessary regulatory complexity for temporal and tissue-specific gene expression. How related factors retain specificity, especially when they recognize the same DNA motifs, is not understood. We address this paradox using basic helix-loop-helix (bHLH) transcription factors ASCL1, ASCL2, and MYOD1, crucial mediators of lineage specification. In vivo, these factors recognize the same DNA motifs, yet bind largely different genomic sites and regulate distinct transcriptional programs. This suggests that their ability to identify regulatory targets is defined either by the cellular environment of the partially defined lineages in which they are endogenously expressed, or by intrinsic properties of the factors themselves. To distinguish between these mechanisms, we directly compared the chromatin binding properties of this subset of bHLH factors when ectopically expressed in embryonic stem cells, presenting them with a common chromatin landscape and cellular components. We find that these factors retain distinct binding sites; thus, specificity of binding is an intrinsic property not requiring a restricted landscape or lineage-specific cofactors. Although the ASCL factors and MYOD1 have some distinct DNA motif preference, it is not sufficient to explain the extent of the differential binding. All three factors can bind inaccessible chromatin and induce changes in chromatin accessibility and H3K27ac. A reiterated pattern of DNA binding motifs is uniquely enriched in inaccessible chromatin at sites bound by these bHLH factors. These combined properties define a subclass of lineage-specific bHLH factors and provide context for their central roles in development and disease.

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.

Function of the myogenic regulatory factors Myf5, MyoD, Myogenin and MRF4 in skeletal muscle, satellite cells and regenerative myogenesis

Discovery of the myogenic regulatory factor family of transcription factors MYF5, MYOD, Myogenin and MRF4 was a seminal step in understanding specification of the skeletal muscle lineage and control of myogenic differentiation during development. These factors are also involved in specification of the muscle satellite cell lineage, which becomes the resident stem cell compartment inadult skeletal muscle. While MYF5, MYOD, Myogenin and MRF4 have subtle roles in mature muscle, they again play a crucial role in directing satellite cell function to regenerate skeletal muscle: linking the genetic control of developmental and regenerative myogenesis. Here, I review the role of the myogenic regulatory factors in developing and mature skeletal muscle, satellite cell specification and muscle regeneration.

Requirement of the fusogenic micropeptide myomixer for muscle formation in zebrafish

Skeletal muscle formation requires fusion of mononucleated myoblasts to form multinucleated myofibers. The muscle-specific membrane proteins myomaker and myomixer cooperate to drive mammalian myoblast fusion. Whereas myomaker is highly conserved across diverse vertebrate species, myomixer is a micropeptide that shows relatively weak cross-species conservation. To explore the functional conservation of myomixer, we investigated the expression and function of the zebrafish myomixer ortholog. Here we show that myomixer expression during zebrafish embryogenesis coincides with myoblast fusion, and genetic deletion of myomixer using CRISPR/Cas9 mutagenesis abolishes myoblast fusion in vivo. We also identify myomixer orthologs in other species of fish and reptiles, which can cooperate with myomaker and substitute for the fusogenic activity of mammalian myomixer. Sequence comparison of these diverse myomixer orthologs reveals key amino acid residues and a minimal fusogenic peptide motif that is necessary for promoting cell-cell fusion with myomaker. Our findings highlight the evolutionary conservation of the myomaker-myomixer partnership and provide insights into the molecular basis of myoblast fusion.

Myomerger induces fusion of non-fusogenic cells and is required for skeletal muscle development

Despite the importance of cell fusion for mammalian development and physiology, the factors critical for this process remain to be fully defined, which has severely limited our ability to reconstitute cell fusion. Myomaker (Tmem8c) is a muscle-specific protein required for myoblast fusion. Expression of myomaker in fibroblasts drives their fusion with myoblasts, but not with other myomaker-expressing fibroblasts, highlighting the requirement of additional myoblast-derived factors for fusion. Here we show that Gm7325, which we name myomerger, induces the fusion of myomaker-expressing fibroblasts. Thus, myomaker and myomerger together confer fusogenic activity to otherwise non-fusogenic cells. Myomerger is skeletal muscle-specific and genetic deletion in mice results in a paucity of muscle fibres demonstrating its requirement for normal muscle formation. Myomerger deficient myocytes differentiate and harbour organized sarcomeres but are fusion-incompetent. Our findings identify myomerger as a fundamental myoblast fusion protein and establish a system that begins to reconstitute mammalian cell fusion.

Cellular fusion is fundamental for skeletal muscle development. Here the authors show that myomerger is expressed in myoblasts, is essential for myoblast fusion in mice, and in co-operation with myomaker confers fusogenic ability to non-fusogenic cells.

Destabilization of the TWIST1/E12 complex dimerization following the R154P point-mutation of TWIST1: an in silico approach

Background

The bHLH transcription factor TWIST1 plays a key role in the embryonic development and in tumorigenesis. Some loss-of-function mutations of the TWIST1 gene have been shown to cause an autosomal dominant craniosynostosis, known as the Saethre-Chotzen syndrome (SCS). Although the functional impacts of many TWIST1 mutations have been experimentally reported, little is known on the molecular mechanisms underlying their loss-of-function. In a previous study, we highlighted the predictive value of in silico molecular dynamics (MD) simulations in deciphering the molecular function of TWIST1 residues.

Results

Here, since the substitution of the arginine 154 amino acid by a glycine residue (R154G) is responsible for the SCS phenotype and the substitution of arginine 154 by a proline experimentally decreases the dimerizing ability of TWIST1, we investigated the molecular impact of this point mutation using MD approaches. Consistently, MD simulations highlighted a clear decrease in the stability of the α-helix during the dimerization of the mutated R154P TWIST1/E12 dimer compared to the wild-type TE complex, which was further confirmed in vitro using immunoassays.

Conclusions

Our study demonstrates that MD simulations provide a structural explanation for the loss-of-function associated with the SCS TWIST1 mutation and provides a proof of concept of the predictive value of these MD simulations. This in silico methodology could be used to determine reliable pharmacophore sites, leading to the application of docking approaches in order to identify specific inhibitors of TWIST1 complexes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0076-x) contains supplementary material, which is available to authorized users.

NeuroD Factors Discriminate Mineralocorticoid From Glucocorticoid Receptor DNA Binding in the Male Rat Brain

In the limbic brain, mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) both function as receptors for the naturally occurring glucocorticoids (corticosterone/cortisol) but mediate distinct effects on cellular physiology via transcriptional mechanisms. The transcriptional basis for specificity of these MR- vs GR-mediated effects is unknown. To address this conundrum, we have identified the extent of MR/GR DNA-binding selectivity in the rat hippocampus using chromatin immunoprecipitation followed by sequencing. We found 918 and 1450 nonoverlapping binding sites for MR and GR, respectively. Furthermore, 475 loci were co-occupied by MR and GR. De novo motif analysis resulted in a similar binding motif for both receptors at 100% of the target loci, which matched the known glucocorticoid response element (GRE). In addition, the Atoh/NeuroD consensus sequence was found in co-occurrence with all MR-specific binding sites but was absent for GR-specific or MR-GR overlapping sites. Basic helix-loop-helix family members Neurod1, Neurod2, and Neurod6 showed hippocampal expression and were hypothesized to bind the Atoh motif. Neurod2 was detected at rat hippocampal MR binding sites but not at GR-exclusive sites. All three NeuroD transcription factors acted as DNA-binding-dependent coactivators for both MR and GR in reporter assays in heterologous HEK293 cells, likely via indirect interactions with the receptors. In conclusion, a NeuroD family member binding to an additional motif near the GRE seems to drive specificity for MR over GR binding at hippocampal binding sites.

Control of muscle formation by the fusogenic micropeptide myomixer

Skeletal muscle formation occurs through fusion of myoblasts to form multinucleated myofibers. From a genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) loss-of-function screen for genes required for myoblast fusion and myogenesis, we discovered an 84-amino acid muscle-specific peptide that we call Myomixer. Myomixer expression coincides with myoblast differentiation and is essential for fusion and skeletal muscle formation during embryogenesis. Myomixer localizes to the plasma membrane, where it promotes myoblast fusion and associates with Myomaker, a fusogenic membrane protein. Myomixer together with Myomaker can also induce fibroblast-fibroblast fusion and fibroblast-myoblast fusion. We conclude that the Myomixer-Myomaker pair controls the critical step in myofiber formation during muscle development.

NEUROD2 Regulates Stim1 Expression and Store-Operated Calcium Entry in Cortical Neurons

Calcium signaling controls many key processes in neurons, including gene expression, axon guidance, and synaptic plasticity. In contrast to calcium influx through voltage- or neurotransmitter-gated channels, regulatory pathways that control store-operated calcium entry (SOCE) in neurons are poorly understood. Here, we report a transcriptional control of Stim1 (stromal interaction molecule 1) gene, which is a major sensor of endoplasmic reticulum (ER) calcium levels and a regulator of SOCE. By using a genome-wide chromatin immunoprecipitation and sequencing approach in mice, we find that NEUROD2, a neurogenic transcription factor, binds to an intronic element within the Stim1 gene. We show that NEUROD2 limits Stim1 expression in cortical neurons and consequently fine-tunes the SOCE response upon depletion of ER calcium. Our findings reveal a novel mechanism that regulates neuronal calcium homeostasis during cortical development.

Lineage-specific stem cells, signals and asymmetries during stomatal development

Stomata are dispersed pores found in the epidermis of land plants that facilitate gas exchange for photosynthesis while minimizing water loss. Stomata are formed from progenitor cells, which execute a series of differentiation events and stereotypical cell divisions. The sequential activation of master regulatory basic-helix-loop-helix (bHLH) transcription factors controls the initiation, proliferation and differentiation of stomatal cells. Cell-cell communication mediated by secreted peptides, receptor kinases, and downstream mitogen-activated kinase cascades enforces proper stomatal patterning, and an intrinsic polarity mechanism ensures asymmetric cell divisions. As we review here, recent studies have provided insights into the intrinsic and extrinsic factors that control stomatal development. These findings have also highlighted striking similarities between plants and animals with regards to their mechanisms of specialized cell differentiation.

Distinct Activities of Myf5 and MyoD Indicate Separate Roles in Skeletal Muscle Lineage Specification and Differentiation

Most transcription factor families contain highly related paralogs generated by gene duplication, and functional divergence is generally accomplished by activation of distinct sets of genes by each member. Here we compare the molecular functions of Myf5 and MyoD, two highly related bHLH transcription factors that regulate skeletal muscle specification and differentiation. We find that MyoD and Myf5 bind the same sites genome-wide but have distinct functions: Myf5 induces histone acetylation without Pol II recruitment or robust gene activation, whereas MyoD induces histone acetylation, recruits Pol II, and robustly activates gene transcription. Therefore, the initial specification of the muscle lineage by Myf5 occurs without significant induction of gene transcription. Transcription of the skeletal muscle program is then achieved by the subsequent expression of MyoD, which binds to the same sites as Myf5, indicating that each factor regulates distinct steps in gene initiation and transcription at a shared set of binding sites.

Epigenetics and developmental programming of welfare and production traits in farm animals

The concept that postnatal health and development can be influenced by events that occur in utero originated from epidemiological studies in humans supported by numerous mechanistic (including epigenetic) studies in a variety of model species. Referred to as the 'developmental origins of health and disease' or 'DOHaD' hypothesis, the primary focus of large-animal studies until quite recently had been biomedical. Attention has since turned towards traits of commercial importance in farm animals. Herein we review the evidence that prenatal risk factors, including suboptimal parental nutrition, gestational stress, exposure to environmental chemicals and advanced breeding technologies, can determine traits such as postnatal growth, feed efficiency, milk yield, carcass composition, animal welfare and reproductive potential. We consider the role of epigenetic and cytoplasmic mechanisms of inheritance, and discuss implications for livestock production and future research endeavours. We conclude that although the concept is proven for several traits, issues relating to effect size, and hence commercial importance, remain. Studies have also invariably been conducted under controlled experimental conditions, frequently assessing single risk factors, thereby limiting their translational value for livestock production. We propose concerted international research efforts that consider multiple, concurrent stressors to better represent effects of contemporary animal production systems.

MyoD reprogramming requires Six1 and Six4 homeoproteins: genome-wide cis-regulatory module analysis

Myogenic regulatory factors of the MyoD family have the ability to reprogram differentiated cells toward a myogenic fate. In this study, we demonstrate that Six1 or Six4 are required for the reprogramming by MyoD of mouse embryonic fibroblasts (MEFs). Using microarray experiments, we found 761 genes under the control of both Six and MyoD. Using MyoD ChIPseq data and a genome-wide search for Six1/4 MEF3 binding sites, we found significant co-localization of binding sites for MyoD and Six proteins on over a thousand mouse genomic DNA regions. The combination of both datasets yielded 82 genes which are synergistically activated by Six and MyoD, with 96 associated MyoD+MEF3 putative cis-regulatory modules (CRMs). Fourteen out of 19 of the CRMs that we tested demonstrated in Luciferase assays a synergistic action also observed for their cognate gene. We searched putative binding sites on these CRMs using available databases and de novo search of conserved motifs and demonstrated that the Six/MyoD synergistic activation takes place in a feedforward way. It involves the recruitment of these two families of transcription factors to their targets, together with partner transcription factors, encoded by genes that are themselves activated by Six and MyoD, including Mef2, Pbx-Meis and EBF.