Loss of Tob1 promotes muscle regeneration through muscle stem cell expansion.

Muscle stem cells (MuSCs) play an indispensable role in postnatal muscle growth and hypertrophy in adult. MuSCs also retain a highly regenerative capacity and are therefore considered a promising stem cell source for regenerative therapy of muscle diseases. In this study, we identify tumor-suppressor protein Tob1 as a Pax7-target gene that negatively controls the population expansion of MuSCs. Tob1 protein is undetectable in the quiescent state but is upregulated during activation in MuSCs. Tob1 ablation in mice accelerates MuSC population expansion and boosts muscle regeneration. Moreover, inactivation of Tob1 in MuSCs ameliorates the efficiency of MuSC transplantation in a murine muscular dystrophy model. Collectively, selective targeting of Tob1 may be a therapeutic option for the treatment of muscular diseases, including muscular dystrophy and age-related sarcopenia.

H4K20me1 and H3K27me3 are concurrently loaded onto the inactive X chromosome but dispensable for inducing gene silencing.

During X chromosome inactivation (XCI), in female placental mammals, gene silencing is initiated by the Xist long non-coding RNA. Xist accumulation at the X leads to enrichment of specific chromatin marks, including PRC2-dependent H3K27me3 and SETD8-dependent H4K20me1. However, the dynamics of this process in relation to Xist RNA accumulation remains unknown as is the involvement of H4K20me1 in initiating gene silencing. To follow XCI dynamics in living cells, we developed a genetically encoded, H3K27me3-specific intracellular antibody or H3K27me3-mintbody. By combining live-cell imaging of H3K27me3, H4K20me1, the X chromosome and Xist RNA, with ChIP-seq analysis we uncover concurrent accumulation of both marks during XCI, albeit with distinct genomic distributions. Furthermore, using a Xist B and C repeat mutant, which still shows gene silencing on the X but not H3K27me3 deposition, we also find a complete lack of H4K20me1 enrichment. This demonstrates that H4K20me1 is dispensable for the initiation of gene silencing, although it may have a role in the chromatin compaction that characterises facultative heterochromatin.

Dmrt factors determine the positional information of cerebral cortical progenitors via differential suppression of homeobox genes.

The spatiotemporal identity of neural progenitors and the regional control of neurogenesis are essential for the development of cerebral cortical architecture. Here, we report that mammalian DM domain factors (Dmrt) determine the identity of cerebral cortical progenitors. Among the Dmrt family genes expressed in the developing dorsal telencephalon, Dmrt3 and Dmrta2 show a medialhigh/laterallow expression gradient. Their simultaneous loss confers a ventral identity to dorsal progenitors, resulting in the ectopic expression of Gsx2 and massive production of GABAergic olfactory bulb interneurons in the dorsal telencephalon. Furthermore, double-mutant progenitors in the medial region exhibit upregulated Pax6 and more lateral characteristics. These ventral and lateral shifts in progenitor identity depend on Dmrt gene dosage. We also found that Dmrt factors bind to Gsx2 and Pax6 enhancers to suppress their expression. Our findings thus reveal that the graded expression of Dmrt factors provide positional information for progenitors by differentially repressing downstream genes in the developing cerebral cortex.

Transcriptome analysis of gene expression changes upon enzymatic dissociation in skeletal myoblasts.

Single-cell RNA-sequencing analysis is one of the most effective tools for understanding specific cellular states. The use of single cells or pooled cells in RNA-seq analysis requires the isolation of cells from a tissue or culture. Although trypsin or more recently cold-active protease (CAP) has been used for cell dissociation, the extent to which the gene expression changes are suppressed has not been clarified. To this end, we conducted detailed profiling of the enzyme-dependent gene expression changes in mouse skeletal muscle progenitor cells, focusing on the enzyme treatment time, amount and temperature. We found that the genes whose expression was changed by the enzyme treatment could be classified in a time-dependent manner and that there were genes whose expression was changed independently of the enzyme treatment time, amount and temperature. This study will be useful as reference data for genes that should be excluded or considered for RNA-seq analysis using enzyme isolation methods.

Modeling population size independent tissue epigenomes by ChIL-seq with single thin sections.

Recent advances in genome-wide technologies have enabled analyses using small cell numbers of even single cells. However, obtaining tissue epigenomes with cell-type resolution from large organs and tissues still remains challenging, especially when the available material is limited. Here, we present a ChIL-based approach for analyzing the diverse cellular dynamics at the tissue level using high-depth epigenomic data. "ChIL for tissues" allows the analysis of a single tissue section and can reproducibly generate epigenomic profiles from several tissue types, based on the distribution of target epigenomic states, tissue morphology, and number of cells. The proposed method enabled the independent evaluation of changes in cell populations and gene activation in cells from regenerating skeletal muscle tissues, using a statistical model of RNA polymerase II distribution on gene loci. Thus, the integrative analyses performed using ChIL can elucidate in vivo cell-type dynamics of tissues.

Discriminative feature of cells characterizes cell populations of interest by a small subset of genes.

Organisms are composed of various cell types with specific states. To obtain a comprehensive understanding of the functions of organs and tissues, cell types have been classified and defined by identifying specific marker genes. Statistical tests are critical for identifying marker genes, which often involve evaluating differences in the mean expression levels of genes. Differentially expressed gene (DEG)-based analysis has been the most frequently used method of this kind. However, in association with increases in sample size such as in single-cell analysis, DEG-based analysis has faced difficulties associated with the inflation of P-values. Here, we propose the concept of discriminative feature of cells (DFC), an alternative to using DEG-based approaches. We implemented DFC using logistic regression with an adaptive LASSO penalty to perform binary classification for discriminating a population of interest and variable selection to obtain a small subset of defining genes. We demonstrated that DFC prioritized gene pairs with non-independent expression using artificial data and that DFC enabled characterization of the muscle satellite/progenitor cell population. The results revealed that DFC well captured cell-type-specific markers, specific gene expression patterns, and subcategories of this cell population. DFC may complement DEG-based methods for interpreting large data sets. DEG-based analysis uses lists of genes with differences in expression between groups, while DFC, which can be termed a discriminative approach, has potential applications in the task of cell characterization. Upon recent advances in the high-throughput analysis of single cells, methods of cell characterization such as scRNA-seq can be effectively subjected to the discriminative methods.

Regulation of RNA polymerase II activation by histone acetylation in single living cells.

In eukaryotic cells, post-translational histone modifications have an important role in gene regulation. Starting with early work on histone acetylation, a variety of residue-specific modifications have now been linked to RNA polymerase II (RNAP2) activity, but it remains unclear if these markers are active regulators of transcription or just passive byproducts. This is because studies have traditionally relied on fixed cell populations, meaning temporal resolution is limited to minutes at best, and correlated factors may not actually be present in the same cell at the same time. Complementary approaches are therefore needed to probe the dynamic interplay of histone modifications and RNAP2 with higher temporal resolution in single living cells. Here we address this problem by developing a system to track residue-specific histone modifications and RNAP2 phosphorylation in living cells by fluorescence microscopy. This increases temporal resolution to the tens-of-seconds range. Our single-cell analysis reveals histone H3 lysine-27 acetylation at a gene locus can alter downstream transcription kinetics by as much as 50%, affecting two temporally separate events. First acetylation enhances the search kinetics of transcriptional activators, and later the acetylation accelerates the transition of RNAP2 from initiation to elongation. Signatures of the latter can be found genome-wide using chromatin immunoprecipitation followed by sequencing. We argue that this regulation leads to a robust and potentially tunable transcriptional response.

Genome-wide analysis of the spatiotemporal regulation of firing and dormant replication origins in human cells.

In metazoan cells, only a limited number of mini chromosome maintenance (MCM) complexes are fired during S phase, while the majority remain dormant. Several methods have been used to map replication origins, but such methods cannot identify dormant origins. Herein, we determined MCM7-binding sites in human cells using ChIP-Seq, classified them into firing and dormant origins using origin data and analysed their association with various chromatin signatures. Firing origins, but not dormant origins, were well correlated with open chromatin regions and were enriched upstream of transcription start sites (TSSs) of transcribed genes. Aggregation plots of MCM7 signals revealed minimal difference in the efficacy of MCM loading between firing and dormant origins. We also analysed common fragile sites (CFSs) and found a low density of origins at these sites. Nevertheless, firing origins were enriched upstream of the TSSs. Based on the results, we propose a model in which excessive MCMs are actively loaded in a genome-wide manner, irrespective of chromatin status, but only a fraction are passively fired in chromatin areas with an accessible open structure, such as regions upstream of TSSs of transcribed genes. This plasticity in the specification of replication origins may minimize collisions between replication and transcription.

Chd2 regulates chromatin for proper gene expression toward differentiation in mouse embryonic stem cells.

Chromatin reorganization is necessary for pluripotent stem cells, including embryonic stem cells (ESCs), to acquire lineage potential. However, it remains unclear how ESCs maintain their characteristic chromatin state for appropriate gene expression upon differentiation. Here, we demonstrate that chromodomain helicase DNA-binding domain 2 (Chd2) is required to maintain the differentiation potential of mouse ESCs. Chd2-depleted ESCs showed suppressed expression of developmentally regulated genes upon differentiation and subsequent differentiation defects without affecting gene expression in the undifferentiated state. Furthermore, chromatin immunoprecipitation followed by sequencing revealed alterations in the nucleosome occupancy of the histone variant H3.3 for developmentally regulated genes in Chd2-depleted ESCs, which in turn led to elevated trimethylation of the histone H3 lysine 27. These results suggest that Chd2 is essential in preventing suppressive chromatin formation for developmentally regulated genes and determines subsequent effects on developmental processes in the undifferentiated state.

Thymine DNA glycosylase modulates DNA damage response and gene expression by base excision repair-dependent and independent mechanisms.

Thymine DNA glycosylase (TDG) is a base excision repair (BER) enzyme, which is implicated in correction of deamination-induced DNA mismatches, the DNA demethylation process and regulation of gene expression. Because of these pivotal roles associated, it is crucial to elucidate how the TDG functions are appropriately regulated in vivo. Here, we present evidence that the TDG protein undergoes degradation upon various types of DNA damage, including ultraviolet light (UV). The UV-induced degradation of TDG was dependent on proficiency in nucleotide excision repair and on CRL4CDT2 -mediated ubiquitination that requires a physical interaction between TDG and DNA polymerase clamp PCNA. Using the Tdg-deficient mouse embryonic fibroblasts, we found that ectopic expression of TDG compromised cellular survival after UV irradiation and repair of UV-induced DNA lesions. These negative effects on cellular UV responses were alleviated by introducing mutations in TDG that impaired its BER function. The expression of TDG induced a large-scale alteration in the gene expression profile independently of its DNA glycosylase activity, whereas a subset of genes was affected by the catalytic activity of TDG. Our results indicate the presence of BER-dependent and BER-independent functions of TDG, which are involved in regulation of cellular DNA damage responses and gene expression patterns.

Structure and function of human histone H3.Y nucleosome.

Histone H3.Y is a primate-specific, distant H3 variant. It is evolutionarily derived from H3.3, and may function in transcription regulation. However, the mechanism by which H3.Y regulates transcription has not been elucidated. In the present study, we determined the crystal structure of the H3.Y nucleosome, and found that many H3.Y-specific residues are located on the entry/exit sites of the nucleosome. Biochemical analyses revealed that the DNA ends of the H3.Y nucleosome were more flexible than those of the H3.3 nucleosome, although the H3.Y nucleosome was stable in vitro and in vivo Interestingly, the linker histone H1, which compacts nucleosomal DNA, appears to bind to the H3.Y nucleosome less efficiently, as compared to the H3.3 nucleosome. These characteristics of the H3.Y nucleosome are also conserved in the H3.Y/H3.3 heterotypic nucleosome, which may be the predominant form in cells. In human cells, H3.Y preferentially accumulated around transcription start sites (TSSs). Taken together, H3.Y-containing nucleosomes around transcription start sites may form relaxed chromatin that allows transcription factor access, to regulate the transcription status of specific genes.

Crystal Structure and Characterization of Novel Human Histone H3 Variants, H3.6, H3.7, and H3.8.

Non-allelic histone variants are considered as epigenetic factors that regulate genomic DNA functions in eukaryotic chromosomes. In this study, we identified three new human histone H3 variants (named H3.6, H3.7, and H3.8), which were previously annotated as pseudogenes. H3.6 and H3.8 conserve the H3.3-specific amino acid residues, but H3.7 shares the specific amino acid residues with H3.1. We successfully reconstituted the nucleosome containing H3.6 in vitro and determined its crystal structure. In the H3.6 nucleosome, the H3.6-specific Val62 residue hydrophobically contacts the cognate H4 molecule, but its contact area is smaller than that of the corresponding H3.3 Ile62 residue. The thermal stability assay revealed that the H3.6 nucleosome is substantially unstable, as compared to the H3.3 nucleosome. Interestingly, mutational analysis demonstrated that the H3.6 Val62 residue is fully responsible for the H3.6 nucleosome instability, probably because of the weakened hydrophobic interaction with H4. We also reconstituted the nucleosome containing H3.8, but its thermal stability was quite low. In contrast, purified H3.7 failed to form nucleosomes in vitro. The identification and characterization of these novel human histone H3 variants provide important new insights into understanding the epigenetic regulation of the human genome.

Incorporation of histone H3.1 suppresses the lineage potential of skeletal muscle.

Lineage potential is triggered by lineage-specific transcription factors in association with changes in the chromatin structure. Histone H3.3 variant is thought to play an important role in the regulation of lineage-specific genes. To elucidate the function of H3.3 in myogenic differentiation, we forced the expression of GFP-H3.1 to alter the balance between H3.1 and H3.3 in mouse C2C12 cells that could be differentiated into myotubes. GFP-H3.1 replaced H3.3 in the regulatory regions of skeletal muscle (SKM) genes and induced a decrease of H3K4 trimethylation (H3K4me3) and increase of H3K27 trimethylation (H3K27me3). Similar results were obtained by H3.3 knockdown. In contrast, MyoD-dependent H3.3 incorporation into SKM genes in fibroblasts induced an increase of H3K4me3 and H3K27me3. In mouse embryos, a bivalent modification of H3K4me3 and H3K27me3 was formed on H3.3-incorporated SKM genes before embryonic skeletal muscle differentiation. These results suggest that lineage potential is established through a selective incorporation of specific H3 variants that governs the balance of histone modifications.

Cdt1-binding protein GRWD1 is a novel histone-binding protein that facilitates MCM loading through its influence on chromatin architecture.

Efficient pre-replication complex (pre-RC) formation on chromatin templates is crucial for the maintenance of genome integrity. However, the regulation of chromatin dynamics during this process has remained elusive. We found that a conserved protein, GRWD1 (glutamate-rich WD40 repeat containing 1), binds to two representative replication origins specifically during G1 phase in a CDC6- and Cdt1-dependent manner, and that depletion of GRWD1 reduces loading of MCM but not CDC6 and Cdt1. Furthermore, chromatin immunoprecipitation coupled with high-throughput sequencing (Seq) revealed significant genome-wide co-localization of GRWD1 with CDC6. We found that GRWD1 has histone-binding activity. To investigate the effect of GRWD1 on chromatin architecture, we used formaldehyde-assisted isolation of regulatory elements (FAIRE)-seq or FAIRE-quantitative PCR analyses, and the results suggest that GRWD1 regulates chromatin openness at specific chromatin locations. Taken together, these findings suggest that GRWD1 may be a novel histone-binding protein that regulates chromatin dynamics and MCM loading at replication origins.

Chd5 Regulates MuERV-L/MERVL Expression in Mouse Embryonic Stem Cells Via H3K27me3 Modification and Histone H3.1/H3.2.

Chd5 is an essential factor for neuronal differentiation and spermatogenesis and is a known tumor suppressor. H3K27me3 and H3K4un are modifications recognized by Chd5; however, it remains unclear how Chd5 remodels chromatin structure. We completely disrupted the Chd5 locus using the CRISPR-Cas9 system to generate a 52 kbp long deletion and analyzed Chd5 function in mouse embryonic stem cells. Our findings show that Chd5 represses murine endogenous retrovirus-L (MuERV-L/MERVL), an endogenous retrovirus-derived retrotransposon, by regulating H3K27me3 and H3.1/H3.2 function.

agplus: a rapid and flexible tool for aggregation plots.

UNLABELLED: Aggregation plots are frequently used to evaluate signal distributions at user-interested points in ChIP-Seq data analysis. agplus, a new and simple command-line tool, enables rapid and flexible generation of text tables tailored for aggregation plots from which users can easily design multiple groups based on user-definitions such as regulatory regions or transcription initiation sites. AVAILABILITY AND IMPLEMENTATION: This software is implemented in Ruby, supported on Linux and Mac OSX, and freely available at http://github.com/kazumits/agplus CONTACT: yohkawa@epigenetics.med.kyushu-u.ac.jp.

Distribution of histone H4 modifications as revealed by a panel of specific monoclonal antibodies.

Post-translational histone modifications play a critical role in genome functions such as epigenetic gene regulation and genome maintenance. The tail of the histone H4 N-terminus contains several amino acids that can be acetylated and methylated. Some of these modifications are known to undergo drastic changes during the cell cycle. In this study, we generated a panel of mouse monoclonal antibodies against histone H4 modifications, including acetylation at K5, K8, K12, and K16, and different levels of methylation at K20. Their specificity was evaluated by ELISA and immunoblotting using synthetic peptide and recombinant proteins that harbor specific modifications or amino acid substitutions. Immunofluorescence confirmed the characteristic distributions of target modifications. An H4K5 acetylation (H4K5ac)-specific antibody CMA405 reacted with K5ac only when the neighboring K8 was unacetylated. This unique feature allowed us to detect newly assembled H4, which is diacetylated at K5 and K12, and distinguish it from hyperacetylated H4, where K5 and K8 are both acetylated. Chromatin immunoprecipiation combined with deep sequencing (ChIP-seq) revealed that acetylation of both H4K8 and H4K16 were enriched around transcription start sites. These extensively characterized and highly specific antibodies will be useful for future epigenetics and epigenome studies.

SraTailor: graphical user interface software for processing and visualizing ChIP-seq data.

Raw data from ChIP-seq (chromatin immunoprecipitation combined with massively parallel DNA sequencing) experiments are deposited in public databases as SRAs (Sequence Read Archives) that are publically available to all researchers. However, to graphically visualize ChIP-seq data of interest, the corresponding SRAs must be downloaded and converted into BigWig format, a process that involves complicated command-line processing. This task requires users to possess skill with script languages and sequence data processing, a requirement that prevents a wide range of biologists from exploiting SRAs. To address these challenges, we developed SraTailor, a GUI (Graphical User Interface) software package that automatically converts an SRA into a BigWig-formatted file. Simplicity of use is one of the most notable features of SraTailor: entering an accession number of an SRA and clicking the mouse are the only steps required to obtain BigWig-formatted files and to graphically visualize the extents of reads at given loci. SraTailor is also able to make peak calls, generate files of other formats, process users' own data, and accept various command-line-like options. Therefore, this software makes ChIP-seq data fully exploitable by a wide range of biologists. SraTailor is freely available at http://www.devbio.med.kyushu-u.ac.jp/sra_tailor/, and runs on both Mac and Windows machines.

Human TREX component Thoc5 affects alternative polyadenylation site choice by recruiting mammalian cleavage factor I.

The transcription-export complex (TREX) couples mRNA transcription, processing and nuclear export. We found that CFIm68, a large subunit of a heterotetrameric protein complex mammalian cleavage factor I (CFIm), which is implicated in alternative polyadenylation site choice, co-purified with Thoc5, a component of human TREX. Immunoprecipitation using antibodies against different components of TREX indicated that most likely both complexes interact via an interaction between Thoc5 and CFIm68. Microarray analysis using human HeLa cells revealed that a subset of genes was differentially expressed on Thoc5 knockdown. Notably, the depletion of Thoc5 selectively attenuated the expression of mRNAs polyadenylated at distal, but not proximal, polyadenylation sites, which phenocopied the depletion of CFIm68. Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) indicated that CFIm68 preferentially associated with the 5' regions of genes; strikingly, the 5' peak of CFIm68 was significantly and globally reduced on Thoc5 knockdown. We suggest a model in which human Thoc5 controls polyadenylation site choice through the co-transcriptional loading of CFIm68 onto target genes.

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.