Search

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.

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.

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.

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.

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.

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.

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.

© 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.

Summary Adaptation to mechanical load, leading to enhanced force and power output, is a characteristic feature of skeletal muscle. Formation of new myonuclei required for efficient muscle hypertrophy relies on prior activation and proliferation of muscle stem cells (MuSCs). However, the mechanisms controlling MuSC expansion under conditions of increased load are not fully understood. Here we demonstrate that interstitial mesenchymal progenitors respond to mechanical load and stimulate MuSC proliferation in a surgical mouse model of increased muscle load. Mechanistically, transcriptional activation of Yes-associated protein 1 (Yap1)/transcriptional coactivator with PDZ-binding motif (Taz) in mesenchymal progenitors results in local production of thrombospondin-1 (Thbs1), which, in turn, drives MuSC proliferation through CD47 signaling. Under homeostatic conditions, however, CD47 signaling is insufficient to promote MuSC proliferation and instead depends on prior downregulation of the Calcitonin receptor. Our results suggest that relayed signaling between mesenchymal progenitors and MuSCs through a Yap1/Taz-Thbs1-CD47 pathway is critical to establish the supply of MuSCs during muscle hypertrophy.

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.

In multicellular organisms, individual cells are characterized by their gene expression profiles and the spatial interactions among cells enable the elaboration of complex functions. Expression profiling in spatially defined regions is crucial to elucidate cell interactions and functions. Here, we established a transcriptome profiling method coupled with photo-isolation chemistry (PIC) that allows the determination of expression profiles specifically from photo-irradiated regions of whole tissues. PIC uses photo-caged oligodeoxynucleotides for in situ reverse transcription. After photo-irradiation of limited areas, gene expression was detected from at least 10 cells in the tissue sections. PIC transcriptome analysis detected genes specifically expressed in small distinct areas of the mouse embryo. Thus, PIC enables transcriptome profiles to be determined from limited regions at a spatial resolution up to the diffraction limit.

We previously demonstrated that CRM1, a major nuclear export factor, accumulates at Hox cluster regions to recruit nucleoporin-fusion protein Nup98HoxA9, resulting in robust activation of Hox genes (Oka et al., 2016). However, whether this phenomenon is general to other leukemogenic proteins remains unknown. Here, we show that two other leukemogenic proteins, nucleoporin-fusion SET-Nup214 and the NPM1 mutant, NPM1c, which contains a nuclear export signal (NES) at its C-terminus and is one of the most frequent mutations in acute myeloid leukemia, are recruited to the HOX cluster region via chromatin-bound CRM1, leading to HOX gene activation in human leukemia cells. Furthermore, we demonstrate that this mechanism is highly sensitive to a CRM1 inhibitor in leukemia cell line. Together, these findings indicate that CRM1 acts as a key molecule that connects leukemogenic proteins to aberrant HOX gene regulation either via nucleoporin-CRM1 interaction (for SET-Nup214) or NES-CRM1 interaction (for NPM1c).

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.

Nucleosomes in contact In eukaryotic cells, genomic DNA must be compacted to fit inside the nucleus. A key player in DNA packaging is the nucleosome, which comprises a segment of DNA wrapped around an octamer of histone proteins. During replication and transcription, nucleosomes must reposition themselves on the DNA. In this process, nucleosomes can collide to form a dinucleosome. Kato et al. report a high-resolution crystal structure of a dinucleosome. One of the octamers has lost a histone dimer so that the dinucleosome comprises an octamer and a hexamer. The structure may represent an intermediate during chromatin remodeling. Science , this issue p. 205

Nucleosomes are dynamic entities that are repositioned along DNA by chromatin remodeling processes. A nucleosome repositioned by the switch-sucrose nonfermentable (SWI/SNF) remodeler collides with a neighbor and forms the intermediate “overlapping dinucleosome.” Here, we report the crystal structure of the overlapping dinucleosome, in which two nucleosomes are associated, at 3.14-angstrom resolution. In the overlapping dinucleosome structure, the unusual “hexasome” nucleosome, composed of the histone hexamer lacking one H2A-H2B dimer from the conventional histone octamer, contacts the canonical “octasome” nucleosome, and they intimately associate. Consequently, about 250 base pairs of DNA are left-handedly wrapped in three turns, without a linker DNA segment between the hexasome and octasome moieties. The overlapping dinucleosome structure may provide important information to understand how nucleosome repositioning occurs during the chromatin remodeling process.

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.

In overloaded and regenerating muscle, the generation of new myonuclei depends on muscle satellite cells (MuSCs). Because MuSC behaviors in these two environments have not been considered separately, MuSC behaviors in overloaded muscle remain unexamined. Here, we show that most MuSCs in overloaded muscle, unlike MuSCs in regenerating muscle, proliferate in the absence of MyoD expression. Mechanistically, MuSCs in overloaded muscle sustain the expression of Heyl, a Notch effector gene, to suppress MyoD expression, which allows effective MuSC proliferation on myofibers and beneath the basal lamina. Although Heyl-knockout mice show no impairment in an injury model, in a hypertrophy model, their muscles harbor fewer new MuSC-derived myonuclei due to increased MyoD expression and diminished proliferation, which ultimately causes blunted hypertrophy. Our results show that sustained HeyL expression is critical for MuSC proliferation specifically in overloaded muscle, and thus indicate that the MuSC-proliferation mechanism differs in overloaded and regenerating muscle.

MCF7 cells acquire estrogen-independent proliferation after long-term estrogen deprivation (LTED), which recapitulates endocrine therapy resistance. LTED cells can become primed for apoptosis, but the underlying mechanism is largely unknown. We previously reported that Eleanor non-coding RNAs (ncRNAs) upregulate the ESR1 gene in LTED cells. Here, we show that Eleanors delineate the topologically associating domain (TAD) of the ESR1 locus in the active nuclear compartment of LTED cells. The TAD interacts with another transcriptionally active TAD, which is 42.9 Mb away from ESR1 and contains a gene encoding the apoptotic transcription factor FOXO3. Inhibition of a promoter-associated Eleanor suppresses all genes inside the Eleanor TAD and the long-range interaction between the two TADs, but keeps FOXO3 active to facilitate apoptosis in LTED cells. These data indicate a role of ncRNAs in chromatin domain regulation, which may underlie the apoptosis-prone nature of therapy-resistant breast cancer cells and could be good therapeutic targets., Long term estrogen deprivation can result in apoptosis in breast cancer cells. Here, the authors show that this apoptosis is induced by the long-range chromatin interaction of loci containing the ESR1 and FOXO3 genes, resulting in FOXO3-mediated apoptosis.

Histone proteins, including their variants, play indispensable roles in chromatin structure and function. The incorporation of specific histones is the key to epigenetic features and contributes to cell differentiation and development. To understand how each histone is targeted to and incorporated into a specific chromatin region, we introduced epitope-tagged histone complexes into permeabilized cells and allowed them to become incorporated into the chromatin. The newly incorporated histone complex is easily distinguished from the parental histones by the epitope-tag. We analyzed the correlation between the histone incorporation and the chromatin structure, using the H2A-H2B and H2A.Z-H2B complexes. The histone incorporation mainly occurred at less condensed chromatin (open), suggesting that the condensed chromatin (closed) is a barrier for histone incorporation. To overcome this barrier, H2A, but not H2A.Z, uses a replication-coupled deposition mechanism. This leads to the genome-wide distribution of H2A and the exclusion of H2A.Z from the closed chromatin. Moreover, H2A.Z-H2B is specifically incorporated into the chromatin around the transcription start sites. However, the H2A.Z mutant, in which the H2A.Z-specific region (M6) is swapped with the corresponding H2A residues, is incorporated into the closed chromatin in a replication-coupled manner. The mutant no longer functions as H2A.Z, in terms of the deposition. This may be the reason why the mutant does not compensate for the lethality of the H2A.Z knockout in flies. Our results show that histone incorporations are regulated by the chromatin structures, which may be crucial for maintaining the epigenetic chromatin states.