Search

PNPLA8, one of the calcium-independent phospholipase A2 enzymes, is involved in various physiological processes through the maintenance of membrane phospholipids. However, little is known about its role in brain development. Here, we report 12 individuals from 10 unrelated families with biallelic ultra-rare variants inPNPLA8presenting with a wide spectrum of clinical features ranging from developmental and epileptic-dyskinetic encephalopathy (DEDE) to progressive movement disorders. Complete loss of PNPLA8 was associated with the severe end of the spectrum, showing DEDE manifestations and congenital or progressive microcephaly. Using cerebral organoids generated from human induced pluripotent stem cells, we found that loss of PNPLA8 reduced the number of basal radial glial cells (bRGCs) and upper-layer neurons. By spatial transcriptomic analysis targeting apical radial glial cells (aRGCs), we found the downregulation of bRGC-related gene sets in patient-derived cerebral organoids. Lipidomic analysis revealed a decrease in the amount of lysophosphatidic acid, lysophosphatidylethanolamine, and phosphatidic acid, indicative of the disturbed phospholipid metabolism inPNPLA8knockout neural progenitor cells. Our data suggest that PNPLA8 has a critical role in the bRGC-mediated expansion of the developing human cortex by regulating the fate commitment of aRGCs.

Neuropathic pain is often caused by injury and diseases that affect the somatosensory system. Although pain development has been well studied, pain recovery mechanisms remain largely unknown. Here, we found that CD11c-expressing spinal microglia appear after the development of behavioral pain hypersensitivity following nerve injury. Nerve-injured mice with spinal CD11c + microglial depletion failed to recover spontaneously from this hypersensitivity. CD11c + microglia expressed insulin-like growth factor-1 (IGF1), and interference with IGF1 signaling recapitulated the impairment in pain recovery. In pain-recovered mice, the depletion of CD11c + microglia or the interruption of IGF1 signaling resulted in a relapse in pain hypersensitivity. Our findings reveal a mechanism for the remission and recurrence of neuropathic pain, providing potential targets for therapeutic strategies.

Liver cancer is an aggressive cancer associated with a poor prognosis. Development of therapeutic strategies for liver cancer requires fundamental research using suitable experimental models. Recent progress in direct reprogramming technology has enabled the generation of many types of cells that are difficult to obtain and provide a cellular resource in experimental models of human diseases. In this study, we aimed to establish a simple one-step method for inducing cells that can form malignant human liver tumors directly from healthy endothelial cells using nonintegrating episomal vectors. To screen for factors capable of inducing liver cancer-forming cells (LCCs), we selected nine genes and one short hairpin RNA that suppresses tumor protein p53 (TP53) expression and introduced them into human umbilical vein endothelial cells (HUVECs), using episomal vectors. To identify the essential factors, we examined the effect of changing the amounts and withdrawing individual factors. We then analyzed the proliferation, gene and protein expression, morphologic and chromosomal abnormality, transcriptome, and tumor formation ability of the induced cells. We found that a set of six factors, forkhead box A3 (FOXA3), hepatocyte nuclear factor homeobox 1A (HNF1A), HNF1B, lin-28 homolog B (LIN28B), MYCL proto-oncogene, bHLH transcription factor (L-MYC), and Kruppel-like factor 5 (KLF5), induced direct conversion of HUVECs into LCCs. The gene expression profile of these induced LCCs (iLCCs) was similar to that of human liver cancer cells, and these cells effectively formed tumors that resembled human combined hepatocellular-cholangiocarcinoma following transplantation into immunodeficient mice. Conclusion: We succeeded in the direct induction of iLCCs from HUVECs by using nonintegrating episomal vectors. iLCCs generated from patients with cancer and healthy volunteers will be useful for further advancements in cancer research and for developing methods for the diagnosis, treatment, and prognosis of liver cancer.

In the field of genomic medical research, the amount of large-scale information continues to increase due to advances in measurement technologies, such as high-performance sequencing and spatial omics, as well as the progress made in genomic cohort studies involving more than one million individuals. Therefore, researchers require more computational resources to analyze this information. Here, we introduce a hybrid cloud system consisting of an on-premise supercomputer, science cloud, and public cloud at the Kyoto University Center for Genomic Medicine in Japan as a solution. This system can flexibly handle various heterogeneous computational resource-demanding bioinformatics tools while scaling the computational capacity. In the hybrid cloud system, we demonstrate the way to properly perform joint genotyping of whole-genome sequencing data for a large population of 11,238, which can be a bottleneck in sequencing data analysis. This system can be one of the reference implementations when dealing with large amounts of genomic medical data in research centers and organizations.

The analysis of gene expression regulation, or the epigenome analysis, at the single-cell level is at the forefront of genomics research. To elucidate the mechanisms that regulate gene expression, chromatin immunoprecipitation has been conventionally used for determining the binding sites of DNA-binding proteins, such as histones and transcription factors. Now several new approaches have been emerged to reveal epigenome states at the single-cell level. Instead of using immunoprecipitation of fragmented chromatin, in situ reactions using cells or nuclei, combining with transposase tagging and other methods, have enabled single-cell analysis. Furthermore, single-cell multiomics techniques to simultaneously profiling transcriptome and open chromatin or histone modification have been developed. These single-cell analyses have the potential to identify different cell types in a cell population and reveal the dynamic changes of gene regulation, although those technologies have not yet reached a level for general application.

SummaryTranscription is a dynamic process that stochastically switches between the ON and OFF states. To detect the dynamic relationship among protein clusters of RNA polymerase II (RNAPII) and coactivators, gene loci, and transcriptional activity, we inserted an MS2 repeat, a TetO repeat, and inteins with a selection marker just downstream of the transcription start site (TSS). By optimizing the individual elements, we have developed the Spliced TetO REpeAt, MS2 repeat, and INtein sandwiched reporter Gene tag (STREAMING-tag) system. Clusters of RNAPII and BRD4 were observed proximally to the TSS of Nanog when the gene was transcribed in mouse embryonic stem cells. In contrast, clusters of MED19 and MED22 Mediator subunits were constitutively located near the TSS. Thus, the STREAMING-tag system revealed the spatiotemporal relationships between transcriptional activity and protein clusters near the gene. This powerful tool is useful for quantitatively understanding dynamic transcriptional regulation in living cells.

In vertebrates, female receptivity to male courtship is highly dependent on ovarian secretion of estrogens and prostaglandins. We recently identified female-specific neurons in the medaka (Oryzias latipes) preoptic area that express Npba, a neuropeptide mediating female sexual receptivity, in response to ovarian estrogens. Here we show by transcriptomic analysis that these neurons express a multitude of neuropeptides, in addition to Npba, in an ovarian-dependent manner, and we thus termed them female-specific, sex steroid-responsive peptidergic (FeSP) neurons. Our results further revealed that FeSP neurons express a prostaglandin E2 receptor gene, ptger4b, in an ovarian estrogen-dependent manner. Behavioral and physiological examination of ptger4b-deficient female medaka found that they exhibit increased sexual receptivity while retaining normal ovarian function and that their FeSP neurons have reduced firing activity and impaired neuropeptide release. Collectively, this work provides evidence that prostaglandin E2/Ptger4b signaling mediates the estrogenic regulation of FeSP neuron activity and female sexual receptivity.

Leydig cells in fetal testes play crucial roles in masculinizing fetuses through androgen production. Gene knockout studies have revealed that growth factors are implicated in fetal Leydig cell (FLC) differentiation, but little is known about the mechanisms regulating this process. We investigate this issue by characterizing FLC progenitor cells using single-cell RNA sequencing. The sequence datasets suggest that thymosin β10 (Tmsb10) is transiently upregulated in the progenitors. While studying the function of Tmsb10, we reveal that platelet-derived growth factor (PDGF) regulates ciliogenesis through the RAS/ERK and PI3K/AKT pathways, and thereby promotes desert hedgehog (DHH)-dependent FLC differentiation. Tmsb10 expressed in the progenitor cells induces their differentiation into FLCs by suppressing the RAS/ERK pathway. Through characterizing the transiently expressed Tmsb10 in the FLC progenitors, this study unveils the molecular process of FLC differentiation and shows that it is cooperatively induced by DHH and PDGF.

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.

Photo-isolation chemistry (PIC) enables isolation of transcriptome information from locally defined areas by photo-irradiation. Here, we present an optimized PIC protocol for formalin-fixed frozen and paraffin mouse sections and fresh-frozen mouse sections. We describe tissue section preparation and permeabilization, followed by in situ reverse transcription using photo-caged primers. We then detail immunostaining and UV-mediated uncaging to the target areas, followed by linear amplification of uncaged cDNAs, library preparation, and quantification. This protocol can be applied to various animal tissue types.

Nucleoporins NUP98 and NUP214 form chimeric fusion proteins that assemble into phase-separated nuclear bodies. However, the function and physiological significance of these nuclear bodies remain largely unknown. Previously, we reported that both NUP98-HOXA9 and SET-NUP214 are recruited to HOX cluster regions via chromatin-bound CRM1, a nuclear export receptor (Oka et al., 2019). Here, we show that these nuclear bodies promote the condensation of mixed lineage leukemia 1 (MLL1), a histone methyltransferase which is essential for the maintenance of HOX gene expression. Our analysis revealed that SET-NUP214 and CRM1 robustly associate with MLL1 to form nuclear bodies and are colocalized on chromatin. We also showed that MLL1 and CRM1 are recruited to the nuclear bodies of NUP98-HOXA9 and that the NUP98-HOXA9/CRM1/MLL1 complex accumulates on its target gene loci, including HOX clusters and MEIS1. These phenomena were not observed in phase-separation–deficient mutants or non-DNA-binding mutants of NUP98-HOXA9. Collectively, these results show that both phase separation and proper targeting of nucleoporin fusions to specific sites could enhance the activation of a wide range of target genes by promoting the condensation of MLL1 and CRM1.

Recent advances have enabled the direct induction of human tissue-specific stem and progenitor cells from differentiated somatic cells. However, it is not known whether human hepatic progenitor cells (hHepPCs) can be generated from other cell types by direct lineage reprogramming with defined transcription factors. Here, we show that a set of three transcription factors, FOXA3, HNF1A, and HNF6, can induce human umbilical vein endothelial cells to directly acquire the properties of hHepPCs. These induced hHepPCs (hiHepPCs) propagate in long-term monolayer culture and differentiate into functional hepatocytes and cholangiocytes by forming cell aggregates and cystic epithelial spheroids, respectively, under three-dimensional culture conditions. After transplantation, hiHepPC-derived hepatocytes and cholangiocytes reconstitute damaged liver tissues and support hepatic function. The defined transcription factors also induce hiHepPCs from endothelial cells circulating in adult human peripheral blood. These expandable and bipotential hiHepPCs may be useful in the study and treatment of human liver diseases., The conditions to induce human hepatic progenitor cells from other cell types are unclear. Here, the authors reprogram human endothelial cells to hepatic progenitor cells by expressing FOXA3, HNF1A and HNF6, capable of giving rise to hepatocytes and cholangiocytes that reconstitute damaged liver tissues on transplantation.

To elucidate dynamic changes in native BCR‐ABL and alternatively spliced tyrosine kinase inhibitor (TKI)‐resistant but function‐dead BCR‐ABLIns35bp variant, following commencement or discontinuation of TKI therapy, each transcript was serially quantified in patients with chronic myeloid leukemia (CML) by deep sequencing. Because both transcripts were amplified together using conventional PCR system for measuring International Scale (IS), deep sequencing method was used for quantifying such BCR‐ABL variants. At the initial diagnosis, 7 of 9 patients presented a small fraction of cells possessing BCR‐ABLIns35bp, accounting for 0.8% of the total IS BCR‐ABL, corresponding to actual BCR‐ABLIns35bp value of 1.1539% IS. TKI rapidly decreased native BCR‐ABL but not BCR‐ABLIns35bp, leading to the initial increase in the proportion of BCR‐ABLIns35bp. Thereafter, both native BCR‐ABL and BCR‐ABLIns35bp gradually decreased in the course of TKI treatment, whereas small populations positive for TKI‐resistant BCR‐ABLIns35bp continued fluctuating at low levels, possibly underestimating the molecular response (MR). Following TKI discontinuation, sequencing analysis of 54 patients revealed a rapid relapse, apparently derived from native BCR‐ABL + clones. However, IS fluctuating at low levels around MR4.0 marked a predominant persistence of cells expressing function‐dead BCR‐ABLIns35bp, suggesting that TKI resumption was unnecessary. We clarified the possible mechanism underlying mis‐splicing BCR‐ABLIns35bp, occurring at the particular pseudo‐splice site within intron8, which can be augmented by TKI treatment through inhibition of RNA polymerase II phosphorylation. No mutations were found in spliceosomal genes. Therefore, monitoring IS functional BCR‐ABL extracting BCR‐ABLIns35bp would lead us to a correct evaluation of MR status, thus determining the adequate therapeutic intervention., Dynamic changes were observed in BCR‐ABL and BCR‐ABLIns35bp after TKI discontinuation. BCR‐ABL+ cells responsible for early relapse after cessation of TKI treatment might derive from the residual native BCR‐ABL addicted clones. In contrast, some patients whose MRD mainly consisted of BCR‐ABLIns35bp transiently lost DMR during a relatively late phase after TKI discontinuation and thereafter regained DMR spontaneously without TKI resumption.

Leydig cells in fetal testes play crucial roles in masculinizing fetuses through androgen production. Gene knockout studies have revealed that growth factors are implicated in fetal Leydig cell (FLC) differentiation, but little is known about the mechanisms regulating this process. We investigated this issue by characterizing FLC progenitor cells using single-cell RNA sequencing. The sequence datasets suggested that thymosin β10 (Tmsb10) was transiently upregulated in the progenitors. While studying the function of Tmsb10, we revealed that platelet-derived growth factor (PDGF) regulated ciliogenesis through the RAS/ERK and PI3K/AKT pathways, and thereby promoted desert hedgehog (DHH)-dependent FLC differentiation. Tmsb10 expressed in the progenitor cells induced their differentiation into FLCs by suppressing the RAS/ERK pathway. Through characterizing the transiently expressed Tmsb10 in the FLC progenitors, this study unveiled the molecular process of FLC differentiation and showed that it is cooperatively induced by DHH and PDGF.

In vertebrates, female receptivity to male courtship is highly dependent on ovarian secretion of estrogens and prostaglandins. We recently identified female-specific neurons in the medaka (Oryzias latipes) preoptic area that express Npba, a neuropeptide mediating female sexual receptivity, in response to ovarian estrogens. Here we show by transcriptomic analysis that these neurons express a multitude of neuropeptides, in addition to Npba, in an ovarian-dependent manner, and we thus termed them female-specific, sex steroid-responsive peptidergic (FeSP) neurons. Our results further revealed that FeSP neurons express a prostaglandin E

Tissue regeneration requires coordination between resident stem cells and local niche cells1,2. Here we identify that senescent cells are integral components of the skeletal muscle regenerative niche that repress regeneration at all stages of life. The technical limitation of senescent-cell scarcity3 was overcome by combining single-cell transcriptomics and a senescent-cell enrichment sorting protocol. We identified and isolated different senescent cell types from damaged muscles of young and old mice. Deeper transcriptome, chromatin and pathway analyses revealed conservation of cell identity traits as well as two universal senescence hallmarks (inflammation and fibrosis) across cell type, regeneration time and ageing. Senescent cells create an aged-like inflamed niche that mirrors inflammation associated with ageing (inflammageing4) and arrests stem cell proliferation and regeneration. Reducing the burden of senescent cells, or reducing their inflammatory secretome through CD36 neutralization, accelerates regeneration in young and old mice. By contrast, transplantation of senescent cells delays regeneration. Our results provide a technique for isolating in vivo senescent cells, define a senescence blueprint for muscle, and uncover unproductive functional interactions between senescent cells and stem cells in regenerative niches that can be overcome. As senescent cells also accumulate in human muscles, our findings open potential paths for improving muscle repair throughout life. We thank M. Jardí, A. Navarro, J. M. Ballestero, K. Slobodnyuk, M. González, J. López and M. Raya for their technical contributions; A. Harada and K. Tanaka for assistance in ATAC-seq; all of the members of the P.M.-C. laboratory for discussions; J. Campisi for p16-3MR mice; J. A. Fernández-Blanco (PRBB Animal Facility); O. Fornas (UPF/CRG FACS Facility); E. Rebollo (IBMB Molecular Imaging Platform); V. A. Raker for manuscript editing; and the members of the Myoage network (A. Maier) for human material. We acknowledge funding from MINECO-Spain (RTI2018-096068, to P.M.-C. and E.P.); ERC-2016-AdG-741966, LaCaixa-HEALTH-HR17-00040, MDA, UPGRADE-H2020-825825, AFM, DPP-Spain, Fundació La MaratóTV3-80/19-202021 and MWRF to P.M.-C.; Fundació La MaratóTV3-137/38-202033 to A.L.S.; María-de-Maeztu Program for Units of Excellence to UPF (MDM-2014-0370) and Severo-Ochoa Program for Centers of Excellence to CNIC (SEV-2015-0505). This work was also supported by JST-CREST JPMJCR16G1 and MEXT/JSPS JP20H00456/18H05527 to Y.O.; the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA16030502) to M.A.E.; V.M. and A.C. were supported by FPI and Maria-de-Maeztu predoctoral fellowships, respectively, and V.S. by a Marie Skłodowska-Curie individual fellowship. Parts of the figures were drawn using pictures from Servier Medical Art. Servier Medical Art by Servier is licensed under a Creative Commons Attribution 3.0 Unported License (https://creativecommons.org/licences/by/3.0/).

Uchino et al. have developed a genetically encoded probe to specifically detect the Ser2-phosphorylated, elongating form of RNA polymerase II in living cells. Ser2-phosphorylated RNA polymerase foci exhibit more dynamic motion than chromatin replication domains and p300-enriched foci., In eukaryotic nuclei, most genes are transcribed by RNA polymerase II (RNAP2), whose regulation is a key to understanding the genome and cell function. RNAP2 has a long heptapeptide repeat (Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7), and Ser2 is phosphorylated on an elongation form. To detect RNAP2 Ser2 phosphorylation (RNAP2 Ser2ph) in living cells, we developed a genetically encoded modification-specific intracellular antibody (mintbody) probe. The RNAP2 Ser2ph-mintbody exhibited numerous foci, possibly representing transcription “factories,” and foci were diminished during mitosis and in a Ser2 kinase inhibitor. An in vitro binding assay using phosphopeptides confirmed the mintbody’s specificity. RNAP2 Ser2ph-mintbody foci were colocalized with proteins associated with elongating RNAP2 compared with factors involved in the initiation. These results support the view that mintbody localization represents the sites of RNAP2 Ser2ph in living cells. RNAP2 Ser2ph-mintbody foci showed constrained diffusional motion like chromatin, but they were more mobile than DNA replication domains and p300-enriched foci, suggesting that the elongating RNAP2 complexes are separated from more confined chromatin domains.

Chromatin profiling in single cells has been extremely challenging and almost exclusively limited to histone proteins. In cases where single-cell methods have shown promise, many require highly specialized equipment or cell type–specific protocols and are relatively low throughput. Here, we combine the advantages of tagmentation, linear amplification, and combinatorial indexing to produce a high-throughput single-cell DNA binding site mapping method that is simple, inexpensive, and capable of multiplexing several independent samples per experiment. Targeted insertion of promoters sequencing (TIP-seq) uses Tn5 fused to proteinA to insert a T7 RNA polymerase promoter adjacent to a chromatin protein of interest. Linear amplification of flanking DNA with T7 polymerase before sequencing library preparation provides ∼10-fold higher unique reads per single cell compared with other methods. We applied TIP-seq to map histone modifications, RNA polymerase II (RNAPII), and transcription factor CTCF binding sites in single human and mouse cells.

Skeletal muscles display sexually dimorphic features. Biochemically, glycolysis and fatty acid β-oxidation occur preferentially in the muscles of males and females, respectively. However, the mechanisms of the selective utilization of these fuels remains elusive. Here, we obtain transcriptomes from quadriceps type IIB fibers of untreated, gonadectomized, and sex steroid-treated mice of both sexes. Analyses of the transcriptomes unveil two genes, Pfkfb3 (phosphofructokinase-2) and Pdk4 (pyruvate dehydrogenase kinase 4), that may function as switches between the two sexually dimorphic metabolic pathways. Interestingly, Pfkfb3 and Pdk4 show male-enriched and estradiol-enhanced expression, respectively. Moreover, the contribution of these genes to sexually dimorphic metabolism is demonstrated by knockdown studies with cultured type IIB muscle fibers. Considering that skeletal muscles as a whole are the largest energy-consuming organs, our results provide insights into energy metabolism in the two sexes, during the estrus cycle in women, and under pathological conditions involving skeletal muscles., Baba et al. analyzed the transcriptomes from quadriceps type IIB fibers of untreated, gonadectomized, and sex steroid-treated mice of both sexes and identified Pfkfb3 and Pdk4 as differentially regulated genes between males and diestrus females. The authors found that Pfkfb3 and Pdk4 may act as metabolic switches, showed male-enriched and estradiol-enhanced expression, respectively and contributed to sexually dimorphic metabolism.

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., Author summary Statistical methods for detecting differences in individual gene expression are indispensable for understanding cell types. However, conventional statistical methods, such as differentially expressed gene (DEG)-based analysis, have faced difficulties associated with the inflation of P-values because of both the large sample size and selection bias introduced by exploratory data analysis such as single-cell transcriptomics. 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 the discrimination of 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 it 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 differentially expressed gene-based methods for interpreting large data sets.

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., “ChIL for tissues” allows obtaining epigenome information from a single‐thin tissue section with high sensitivity. A statistical model of RNApolII distribution is combined with ChIL to analyze cell‐type‐specific transcriptional states in a tissue.

Histone H3mm18 is a non-allelic H3 variant expressed in skeletal muscle and brain in mice. However, its function has remained enigmatic. We found that H3mm18 is incorporated into chromatin in cells with low efficiency, as compared to H3.3. We determined the structures of the nucleosome core particle (NCP) containing H3mm18 by cryo-electron microscopy, which revealed that the entry/exit DNA regions are drastically disordered in the H3mm18 NCP. Consistently, the H3mm18 NCP is substantially unstable in vitro. The forced expression of H3mm18 in mouse myoblast C2C12 cells markedly suppressed muscle differentiation. A transcriptome analysis revealed that the forced expression of H3mm18 affected the expression of multiple genes, and suppressed a group of genes involved in muscle development. These results suggest a novel gene expression regulation system in which the chromatin landscape is altered by the formation of unusual nucleosomes with a histone variant, H3mm18, and provide important insight into understanding transcription regulation by chromatin.

DNA replication is a key step in initiating cell proliferation. Loading hexameric complexes of minichromosome maintenance (MCM) helicase onto DNA replication origins during the G1 phase is essential for initiating DNA replication. Here, we examined MCM hexamer states during the cell cycle in human hTERT-RPE1 cells using multicolor immunofluorescence-based, single-cell plot analysis, and biochemical size fractionation. Experiments involving cell-cycle arrest at the G1 phase and release from the arrest revealed that a double MCM hexamer was formed via a single hexamer during G1 progression. A single MCM hexamer was recruited to chromatin in the early G1 phase. Another single hexamer was recruited to form a double hexamer in the late G1 phase. We further examined relationship between the MCM hexamer states and the methylation levels at lysine 20 of histone H4 (H4K20) and found that the double MCM hexamer state was correlated with di/trimethyl-H4K20 (H4K20me2/3). Inhibiting the conversion from monomethyl-H4K20 (H4K20me1) to H4K20me2/3 retained the cells in the single MCM hexamer state. Non-proliferative cells, including confluent cells or Cdk4/6 inhibitor-treated cells, also remained halted in the single MCM hexamer state. We propose that the single MCM hexamer state is a halting step in the determination of cell cycle progression.

The musculoskeletal system is integrated by tendons that are characterized by the expression of scleraxis (Scx), a functionally important transcription factor. Here, we newly developed a tenocyte induction method using induced pluripotent stem cells established from ScxGFP transgenic mice by monitoring fluorescence, which reflects a dynamic differentiation process. Among several developmentally relevant factors, transforming growth factor-beta 2 (TGF-β2) was the most potent inducer for differentiation of tenomodulin-expressing mature tenocytes. Single-cell RNA sequencing (scRNA-seq) revealed 11 distinct clusters, including mature tenocyte population and tenogenic differentiation trajectory, which recapitulated the in vivo developmental process. Analysis of the scRNA-seq dataset highlighted the importance of retinoic acid (RA) as a regulatory pathway of tenogenic differentiation. RA signaling was shown to have inhibitory effects on entheseal chondrogenic differentiation as well as TGF-β2-dependent tenogenic/fibrochondrogenic differentiation. The collective findings provide a new opportunity for tendon research and further insight into the mechanistic understanding of the differentiation pathway to a tenogenic fate.

The precursor of heme, protoporphyrin IX (PPIX), accumulates abundantly in the uteri of birds, such as Japanese quail, Coturnix japonica, which has brown-speckled eggshells; however, the molecular basis of PPIX production in the uterus remains largely unknown. Here, we investigated the cause of low PPIX production in a classical Japanese quail mutant exhibiting white eggshells by comparing its gene expression in the uterus with that of the wild type using transcriptome analysis. We also performed genetic linkage analysis to identify the causative genomic region of the white eggshell phenotype. We found that 11 genes, including 5’-aminolevulinate synthase 1 (ALAS1) and hephaestin-like 1 (HEPHL1), were specifically upregulated in the wild-type uterus and downregulated in the mutant. We mapped the 172 kb candidate genomic region on chromosome 6, which contains several genes, including a part of the paired-like homeodomain 3 (PITX3), which encodes a transcription factor. ALAS1, HEPHL1, and PITX3 were expressed in the apical cells of the luminal epithelium and lamina propria cells of the uterine mucosa of the wild-type quail, while their expression levels were downregulated in the cells of the mutant quail. Biochemical analysis using uterine homogenates indicated that the restricted availability of 5’-aminolevulinic acid is the main cause of low PPIX production. These results suggest that uterus-specific transcriptional regulation of heme-biosynthesis-related genes is an evolutionarily acquired mechanism of eggshell pigment production in Japanese quail. Based on these findings, we discussed the molecular basis of PPIX production in the uteri of Japanese quails.

The CCAAT motif-binding factor NF-Y consists of three different subunits, NF-YA, NF-YB, and NF-YC. Although it is suggested that NF-Y activity is essential for normal tissue homeostasis, survival, and metabolic function, its precise role in lipid metabolism is not clarified yet. In Drosophila, eye disc specific knockdown of Drosophila NF-YA (dNF-YA) induced aberrant morphology of the compound eye, the rough eye phenotype in adults and mutation of the lipase 4 (lip4) gene suppressed the rough eye phenotype. RNA-seq analyses with dNF-YA knockdown third instar larvae identified the lip4 gene as one of the genes that are up-regulated by the dNF-YA knockdown. We identified three dNF-Y-binding consensuses in the 5'flanking region of the lip4 gene, and a chromatin immunoprecipitation assay with the specific anti-dNF-YA IgG demonstrated dNF-Y binding to this genomic region. The luciferase transient expression assay with cultured Drosophila S2 cells and the lip4 promoter-luciferase fusion genes with and without mutations in the dNF-Y-binding consensuses showed that each of the three dNF-Y consensus sequences negatively regulated lip4 gene promoter activity. Consistent with these results, qRT-PCR analysis with the dNF-YA knockdown third instar larvae revealed that endogenous lip4 mRNA levels were increased by the knockdown of dNF-YA in vivo. The specific knockdown of dNF-YA in the fat body with the collagen-GAL4 driver resulted in smaller oil droplets in the fat body cells. Collectively, these results suggest that dNF-Y is involved in lipid storage through its negative regulation of lip4 gene transcription.

To gain insights into tissue-specific gene expression in multicellular systems, gene expression profiles are required to be precisely linked with spatial information. Here, we establish a hight-depth spatial transcriptomics method, photo-isolation chemistrty (PIC), which is able to isolate gene expression profiles only from UV-irradiatied region out of whole tissue section. This method performs reverse transcription on tissue sections using photocaged oligo DNAs. After the UV irradiation, the cDNAs in the irradiated regions are allowed to be amplified and sequenced, thereby providing gene expression profiles linked with spatial information.

In multicellular organisms, expression profiling in spatially defined regions is crucial to elucidate cell interactions and functions. Here, we establish a transcriptome profiling method coupled with photo-isolation chemistry (PIC) that allows the determination of expression profiles specifically from photo-irradiated regions of interest. PIC uses photo-caged oligodeoxynucleotides for in situ reverse transcription. PIC transcriptome analysis detects genes specifically expressed in small distinct areas of the mouse embryo. Photo-irradiation of single cells demonstrated that approximately 8, 000 genes were detected with 7 × 10⁴ unique read counts. Furthermore, PIC transcriptome analysis is applicable to the subcellular and subnuclear microstructures (stress granules and nuclear speckles, respectively), where hundreds of genes can be detected as being specifically localised. The spatial density of the read counts is higher than 100 per square micrometre. Thus, PIC enables high-depth transcriptome profiles to be determined from limited regions up to subcellular and subnuclear resolutions., 光照射を用いた超高解像度な遺伝子解析技術の開発に成功 --組織内に潜むがん細胞の病理診断などに応用可能--. 京都大学プレスリリース. 2021-07-27.

During neocortical development, tight regulation of neurogenesis-to-astrogenesis switching of neural precursor cells (NPCs) is critical to generate a balanced number of each neural cell type for proper brain functions. Accumulating evidence indicates that a complex array of epigenetic modifications and the availability of extracellular factors control the timing of neuronal and astrocytic differentiation. However, our understanding of NPC fate regulation is still far from complete. Bone morphogenetic proteins (BMPs) are renowned as cytokines that induce astrogenesis of gliogenic late-gestational NPCs. They also promote neurogenesis of mid-gestational NPCs, although the underlying mechanisms remain elusive. By performing multiple genome-wide analyses, we demonstrate that Smads, transcription factors that act downstream from BMP signaling, target dramatically different genomic regions in neurogenic and gliogenic NPCs. We found that histone H3K27 trimethylation and DNA methylation around Smad-binding sites change rapidly as gestation proceeds, strongly associated with the alteration of accessibility of Smads to their target binding sites. Furthermore, we identified two lineage-specific Smad-interacting partners—Sox11 for neurogenic and Sox8 for astrocytic differentiation—that further ensure Smad-regulated fate-specific gene induction. Our findings illuminate an exquisite regulation of NPC property change mediated by the interplay between cell-extrinsic cues and -intrinsic epigenetic programs during cortical development.

The ovulatory luteinizing hormone (LH) surge induces rapid changes of gene expression and cellular functions in granulosa cells (GCs) undergoing luteinization. However, it remains unclear how the changes in genome-wide gene expression are regulated. H3K4me3 histone modifications are involved in the rapid alteration of gene expression. In this study, we investigated genome-wide changes of transcriptome and H3K4me3 status in mouse GCs undergoing luteinization. GCs were obtained from mice treated with equine chorionic gonadotropin (hCG) before, 4 hours, and 12 hours after human chorionic gonadotropin injection. RNA-sequencing identified a number of upregulated and downregulated genes, which could be classified into 8 patterns according to the time-course changes of gene expression. Many genes were transiently upregulated or downregulated at 4 hours after hCG stimulation. Gene Ontology terms associated with these genes included steroidogenesis, ovulation, cumulus-oocyte complex (COC) expansion, angiogenesis, immune system, reactive oxygen species (ROS) metabolism, inflammatory response, metabolism, and autophagy. The cellular functions of DNA repair and cell growth were newly identified as being activated during ovulation. Chromatin immunoprecipitation–sequencing revealed a genome-wide and rapid change in H3K4me3 during ovulation. Integration of transcriptome and H3K4me3 data identified many H3K4me3-associated genes that are involved in steroidogenesis, ovulation, COC expansion, angiogenesis, inflammatory response, immune system, ROS metabolism, lipid and glucose metabolism, autophagy, and regulation of cell size. The present results suggest that genome-wide changes in H3K4me3 after the LH surge are associated with rapid changes in gene expression in GCs, which enables GCs to acquire a lot of cellular functions within a short time that are required for ovulation and luteinization.

Adult muscles and their associated stem cells retain functional positional memory based on developmental origin., Muscle stem cells (satellite cells) are distributed throughout the body and have heterogeneous properties among muscles. However, functional topographical genes in satellite cells of adult muscle remain unidentified. Here, we show that expression of Homeobox-A (Hox-A) cluster genes accompanied with DNA hypermethylation of the Hox-A locus was robustly maintained in both somite-derived muscles and their associated satellite cells in adult mice, which recapitulates their embryonic origin. Somite-derived satellite cells were clearly separated from cells derived from cranial mesoderm in Hoxa10 expression. Hoxa10 inactivation led to genomic instability and mitotic catastrophe in somite-derived satellite cells in mice and human. Satellite cell–specific Hoxa10 ablation in mice resulted in a decline in the regenerative ability of somite-derived muscles, which were unobserved in cranial mesoderm–derived muscles. Thus, our results show that Hox gene expression profiles instill the embryonic history in satellite cells as positional memory, potentially modulating region-specific pathophysiology in adult muscles.

Morphogen signalling forms an activity gradient and instructs cell identities in a signalling strength-dependent manner to pattern developing tissues. However, developing tissues also undergo dynamic morphogenesis, which may produce cells with unfit morphogen signalling and consequent noisy morphogen gradients. Here we show that a cell competition-related system corrects such noisy morphogen gradients. Zebrafish imaging analyses of the Wnt/β-catenin signalling gradient, which acts as a morphogen to establish embryonic anterior-posterior patterning, identify that unfit cells with abnormal Wnt/β-catenin activity spontaneously appear and produce noise in the gradient. Communication between unfit and neighbouring fit cells via cadherin proteins stimulates apoptosis of the unfit cells by activating Smad signalling and reactive oxygen species production. This unfit cell elimination is required for proper Wnt/β-catenin gradient formation and consequent anterior-posterior patterning. Because this gradient controls patterning not only in the embryo but also in adult tissues, this system may support tissue robustness and disease prevention., Gradients of morphogens such as Wnt provide instructive cues for cell identities during development. Here, the authors report that in the developing zebrafish embryo, cell competition and elimination of unfit cells are required for proper Wnt gradient formation.

Although the postoperative recurrence rate for pancreatic neuroendocrine tumors (PNETs) is reported to be 13.5%‐30%, the paucity of valuable biomarkers to predict recurrence poses a problem for the early detection of relapse. Hence, this study aimed to identify new biomarkers to predict the recurrence of PNETs. We performed RNA sequencing (RNA‐Seq) on RNA isolated from frozen primary tumors sampled from all localized G1/G2 PNETs resected curatively from 1998 to 2015 in our institution. We calculated differentially expressed genes (DEGs) in tumor with and without recurrence (≥3 years) for the propensity‐matched cohort. Gene ontology analysis for the identified DEGs was also performed. Furthermore, we evaluated the expression levels of candidate genes as recurrence predictors via immunostaining. Comparison of transcriptional levels in tumors with and without recurrence identified 166 DEGs. Up‐ and downregulated genes with high significance in these tumors were mainly related to extracellular organization and cell adhesion, respectively. We observed the top three upregulated genes, C‐type lectin domain family 3 member A (CLEC3A), matrix metalloproteinase‐7 (MMP7), and lipocalin2 (LCN2) immunohistochemically and compared their levels in recurrent and nonrecurrent tumors. Significantly higher recurrence rate was shown in patients with positive expression of CLEC3A (P = 0.028), MMP7 (P = 0.003), and LCN2 (P = 0.040) than that with negative expression. We identified CLEC3A, MMP7, and LCN2 known to be associated with the phosphatidylinositol‐3‐kinase/Akt pathway, as potential novel markers to predict the postoperative recurrence of PNETs.

A set of key developmental genes is essential for skeletal growth from multipotent progenitor cells at weaning. Polycomb group proteins, which regulate such genes contributes to the cell lineage commitment and subsequent differentiation via epigenetic chromatin modification and remodeling. However, it is unclear which cell lineage and gene sets are targeted by polycomb proteins during skeletal growth. We now report that mice deficient in a polycomb group gene Cbx2cterm/cterm exhibited skeletal hypoplasia in the tibia, femur, and cranium. Long bone cavities in these mice contained fewer multipotent mesenchymal stromal cells. RNA-sequencing of bone marrow cells showed downregulation and upregulation of osteoblastic and adipogenic genes, respectively. Furthermore, the expression levels of genes specifically expressed in B-cell precursors were decreased. Forced expression of Cbx2 in Cbx2cterm/cterm bone marrow stromal cell recovered fibroblastic colony formation and suppressed adipogenic differentiation. Collectively, our results suggest that Cbx2 controls the maintenance and adipogenic differentiation of mesenchymal stromal cells in the bone marrow.

Background. We previously demonstrated that step training leads to reorganization of neuronal networks in the lumbar spinal cord of rodents after a hemisection (HX) injury and step training, including increases excitability of spinally evoked potentials in hindlimb motor neurons. Methods. In this study, we investigated changes in RNA expression and synapse number using RNA-Seq and immunohistochemistry of the lumbar spinal cord 23 days after a mid-thoracic HX in rats with and without post-HX step training. Results. Gene Ontology (GO) term clustering demonstrated that expression levels of 36 synapse-related genes were increased in trained compared with nontrained rats. Many synaptic genes were upregulated in trained rats, but Lrrc4 (coding NGL-2) was the most highly expressed in the lumbar spinal cord caudal to the HX lesion. Trained rats also had a higher number of NGL-2/synaptophysin synaptic puncta in the lumbar ventral horn. Conclusions. Our findings demonstrate clear activity-dependent regulation of synapse-related gene expression post-HX. This effect is consistent with the concept that activity-dependent phenomena can provide a mechanistic drive for epigenetic neuronal group selection in the shaping of the reorganization of synaptic networks to learn the locomotion task being trained after spinal cord injury.

Despite recent therapeutic advances for multiple myeloma (MM), relapse is very common. Here, we conducted longitudinal single-cell transcriptome sequencing (scRNA-seq) of MM cells from a patient with relapsed MM, treated with multiple anti-myeloma drugs. We observed five subclusters of MM cells, which appeared and/or disappeared in response to the therapeutic pressure, and identified cluster 3 which emerged during lenalidomide treatment and disappeared after proteasome inhibitor (PI) treatment. Among the differentially expressed genes in cluster 3, we found a candidate drug-response gene; pellino E3 ubiquitin-protein ligase family member 2 (

Chromatin plays a crucial role in gene regulation, and chromatin immunoprecipitation followed by sequencing (ChIP–seq) has been the standard technique for examining protein–DNA interactions across the whole genome. However, it is difficult to obtain epigenomic information from limited numbers of cells by ChIP–seq because of sample loss during chromatin preparation and inefficient immunoprecipitation. In this study, we established an immunoprecipitation-free epigenomic profiling method named chromatin integration labelling (ChIL), which enables the amplification of genomic sequences closely associated with the target molecules before cell lysis. Using ChIL followed by sequencing (ChIL–seq), we reliably detected the distributions of histone modifications and DNA-binding factors in 100–1,000 cells. In addition, ChIL–seq successfully detected genomic regions associated with histone marks at the single-cell level. Thus, ChIL–seq offers an alternative method to ChIP–seq for epigenomic profiling using small numbers of cells, in particular, those attached to culture plates and after immunofluorescence. Harada et al. develop a chromatin integration labelling (ChIL) method to map distributions of histone modifications and DNA-binding factors at low-input or even single-cell levels.

In eukaryotes, histone variant distribution within the genome is the key epigenetic feature. To understand how each histone variant is targeted to the genome, we developed a new method, the RhIP (Reconstituted histone complex Incorporation into chromatin of Permeabilized cell) assay, in which epitope-tagged histone complexes are introduced into permeabilized cells and incorporated into their chromatin. Using this method, we found that H3.1 and H3.3 were incorporated into chromatin in replication-dependent and -independent manners, respectively. We further found that the incorporation of histones H2A and H2A.Z mainly occurred at less condensed chromatin (open), suggesting that condensed chromatin (closed) is a barrier for histone incorporation. To overcome this barrier, H2A, but not H2A.Z, uses a replication-coupled deposition mechanism. Our study revealed that the combination of chromatin structure and DNA replication dictates the differential histone deposition to maintain the epigenetic chromatin states.

In eukaryotic nuclei, most genes are transcribed by RNA polymerase II (RNAP2). How RNAP2 transcription is regulated in the nucleus is a key to understanding the genome and cell function. The largest subunit of RNAP2 has a long heptapeptide repeat (Tyr1-Ser2-Pro3-Thr4-Ser5- Pro6-Ser7) at the C-terminal domain and Ser2 is phosphorylated on an elongation form of RNAP2. To detect RNAP2 Ser2 phosphorylation (RNAP2 Ser2ph) in living cells, we developed a genetically encoded modification-specific intracellular antibody (mintbody) probe. The RNAP2 Ser2ph-mintbody probe exhibited numerous foci, possibly representing transcription “factories” in living HeLa cells, and foci were diminished when cells were treated with triptolide to induce RNAP2 degradation and with flavopiridol to inhibit Ser2ph. An in vitro binding assay using phospho-peptides confirmed the Ser2ph-specific binding of the mintbody. These results support the view that mintbody localization represents the sites of RNAP2 Ser2ph in living cells. RNAP2 Ser2ph-mintbody foci were colocalized with proteins associated with elongating RNAP2, such as the CDK12 and Paf1 complex component, compared to factors involved in transcription activation around the transcription start sites, such as CDK9 and BRD4. Tracking analysis revealed that RNAP2 Ser2ph-mintbody foci showed constrained diffusional motion like chromatin, but was more mobile compared to euchromatin domains, suggesting that the elongating RNAP2 complexes are separated from the more confined initiating clusters.SummaryThe authors developed a genetically encoded probe to specifically detect the Ser2- phosphorylated, elongating form of RNA Polymerase II in living cells. The motion of Ser2- phosphorylated polymerase foci was more dynamic than chromatin domains, suggesting that the elongating complexes are separated from the more confined initiating clusters.

DNA methylation is an essential form of epigenetic regulation responsible for cellular identity. In muscle stem cells, termed satellite cells, DNA methylation patterns are tightly regulated during differentiation. However, it is unclear how these DNA methylation patterns affect the function of satellite cells. We demonstrate that a key epigenetic regulator, ubiquitin like with PHD and RING finger domains 1 (Uhrf1), is activated in proliferating myogenic cells but not expressed in quiescent satellite cells or differentiated myogenic cells in mice. Ablation of Uhrf1 in mouse satellite cells impairs their proliferation and differentiation, leading to failed muscle regeneration.

Summary An effective combination of multi-omic datasets can enhance our understanding of complex biological phenomena. To build a context-dependent network with multiple omic layers, i.e., a trans-omic network, we perform phosphoproteomics, transcriptomics, proteomics, and metabolomics of murine liver for 4 h after insulin administration and integrate the resulting time series. Structural characteristics and dynamic nature of the network are analyzed to elucidate the impact of insulin. Early and prominent changes in protein phosphorylation and persistent and asynchronous changes in mRNA and protein levels through non-transcriptional mechanisms indicate enhanced crosstalk between phosphorylation-mediated signaling and protein expression regulation. Metabolic response shows different temporal regulation with transient increases at early time points across categories and enhanced response in the amino acid and nucleotide categories at later time points as a result of process convergence. This extensive and dynamic view of the trans-omic network elucidates prominent regulatory mechanisms that drive insulin responses through intricate interlayer coordination.

The SRY gene induces testis development even in XX individuals. However, XX/Sry testes fail to produce mature sperm, due to the absence of Y chromosome carrying genes essential for spermatogenesis. XX/Sry Sertoli cells show abnormalities in the production of lactate and cholesterol required for germ cell development. Leydig cells are essential for male functions through testosterone production. However, whether XX/Sry adult Leydig cells (XX/Sry ALCs) function normally remains unclear. In this study, the transcriptomes from XY and XX/Sry ALCs demonstrated that immediate early and cholesterogenic gene expressions differed between these cells. Interestingly, cholesterogenic genes were upregulated in XX/Sry ALCs, although downregulated in XX/Sry Sertoli cells. Among the steroidogenic enzymes, CYP17A1 mediates steroid 17α-hydroxylation and 17,20-lyase reaction, necessary for testosterone production. In XX/Sry ALCs, the latter reaction was selectively decreased. The defects in XX/Sry ALCs, together with those in the germ and Sertoli cells, might explain the infertility of XX/Sry testes.

© 2021 EMBO. This is an open access article under the terms of the Creative Commons Attribution License,which permits use, distribution and reproduction in any medium, provided the original work is properly cited., 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., This work was supported by Fundação para a Ciência e Tecnologia (S.T.d.R), project grants PTDC/BIA‐ MOL/29320/2017 IC&DT (A. C. R. & S.T.d.R), CEECUIND/01234/207 (S.T.d.R), and SFRH/BD/137099/2018 (A.C.R.), by an ERC Advanced Investigator award ERC‐ADG‐2014 671027 attributed to E.H., Sir Henry Wellcome Postdoctoral Fellowship (J.J.Z.), Japan Society for the Promotion of Science KAKENHI grants (JP17KK0143 and JP20K06484 to Y.S., JP19H04970, JP19H03158 and JP20H05393 to K.M., JP17K17719 to T.H., JP18H05534 to H.Ku, JP18H05527 and JP20H00456 to Y.O., JP17H01417 and JP18H05527 to H.Ki), and Japan Science and Technology Agency (JST) CREST JPMJCR16G1 to T.K., H.Ku, Y.O. and H.Ki, PREST JPMJPR2026 to K.M., and ERATO JPMJER1901 to H.Ku. J.J.Z. is supported by core funding of The Novo Nordisk Foundation Center for Stem Cell Biology (Novo Nordisk Foundation grant number NNF17CC0027852). Open Access funding enabled and organized by Projekt DEAL.

Recent advances in omics studies have enabled analysis at the single-cell level; however, methods for analyzing the whole cell of large organs and tissues remain challenging. Here, we developed a method named tsChIL to understand the diverse cellular dynamics at the tissue level using high-depth epigenomic data. tsChIL allowed the analysis of a single tissue section and could reproducibly acquire epigenomic profiles from several types of tissues, 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 of cells in regenerating skeletal muscle tissues, using a statistical model of RNA polymerase II distribution on gene loci. Thus, the integrative analysis by tsChIL can elucidate in vivo cell-type dynamics of tissues.

MyoD, a myogenic differentiation protein, has been studied for its critical role in skeletal muscle differentiation. MyoD-expressing myoblasts have a potency to be differentiated with proliferation of ectopic cells. However, little is known about the effect on chromatin structure of MyoD binding in proliferative myoblasts. In this study, we evaluated the chromatin structure around MyoD-bound genome regions during the cell cycle by chromatin immunoprecipitation sequencing. Genome-wide analysis of histone modifications was performed in proliferative mouse C2C12 myoblasts during three phases (G1, S, G2/M) of the cell cycle. We found that MyoD-bound genome regions had elevated levels of active histone modifications, such as H3K4me1/2/3 and H3K27ac, compared with MyoD-unbound genome regions during the cell cycle. We also demonstrated that the elevated H3K4me2/3 modification level was maintained during the cell cycle, whereas the H3K27ac and H3K4me1 modification levels decreased to the same level as MyoD-unbound genome regions during the later phases. Immunoblot analysis revealed that MyoD abundance was high in the G1 phase then decreased in the S and G2/M phases. Our results suggest that MyoD binding formed selective epigenetic memories with H3K4me2/3 during the cell cycle in addition to myogenic gene induction via active chromatin formation coupled with transcription.

During X chromosome inactivation (XCI), in female placental mammals, gene silencing is initiated by the Xist long-noncoding 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 molecular mechanism allowing for H4K20me1 enrichment. To follow XCI dynamics in living cells, we developed a genetically-encoded, H3K27me3-specific intracellular antibody, or H3K27me3-mintbody. By combining it with live-imaging of H4K20me1, the X chromosome and Xist RNA we uncover similarities in the initial accumulation dynamics of H3K27me3 and H4K20me1. Further ChIP-seq analysis confirmed concurrent accumulation of both marks during XCI albeit with distinct genomic distributions. Using a Xist B and C repeat mutant, which can silence the X but does not allow for H3K27me3 deposition, we also found a lack of H4K20me1 enrichment. Thus, these two marks accumulate at the X thanks to the same region of Xist and H4K20me1 in particular may have a role in the chromatin compaction that characterises facultative heterochromatin.