Models of global gene expression define major domains of cell type and tissue identity

Rapid retinoic acid-induced trophoblast cell model from human induced pluripotent stem cells

A limited number of accessible and representative models of human trophoblast cells currently exist for the study of placentation. Current stem cell models involve either a transition through a naïve stem cell state or precise dynamic control of multiple growth factors and small-molecule cues. Here, we demonstrated that a simple five-day treatment of human induced pluripotent stem cells with two small molecules, retinoic acid (RA) and Wnt agonist CHIR 99021 (CHIR), resulted in rapid, synergistic upregulation of CDX2. Transcriptomic analysis of RA + CHIR-treated cells showed high similarity to primary trophectoderm cells. Multipotency was verified via further differentiation towards cells with syncytiotrophoblast or extravillous trophoblast features. RA + CHIR-treated cells were also assessed for the established criteria defining a trophoblast cell model, and they possess all the features necessary to be considered valid. Collectively, our data demonstrate a facile, scalable method for generating functional trophoblast-like cells in vitro to better understand the placenta.

In vitro generation of mouse morula-like cells

Generating cells with the molecular and functional properties of embryo cells and with full developmental potential is an aim with fundamental biological significance. Here we report the in vitro generation of mouse transient morula-like cells (MLCs) via the manipulation of signaling pathways. MLCs are molecularly distinct from embryonic stem cells (ESCs) and cluster instead with embryo 8- to 16-cell stage cells. A single MLC can generate a blastoid, and the efficiency increases to 80% when 8-10 MLCs are used. MLCs make embryoids directly, efficiently, and within 4 days. Transcriptomic analysis shows that day 4-5 MLC-derived embryoids contain the cell types found in natural embryos at early gastrulation. Furthermore, MLCs introduced into morulae segregate into epiblast (EPI), primitive endoderm (PrE), and trophectoderm (TE) fates in blastocyst chimeras and have a molecular signature indistinguishable from that of host embryo cells. These findings represent the generation of cells that are molecularly and functionally similar to the precursors of the first three cell lineages of the embryo.

c-JUN is a barrier in hESC to cardiomyocyte transition

Loss of c-JUN leads to early mouse embryonic death, possibly because of a failure to develop a normal cardiac system. How c-JUN regulates human cardiomyocyte cell fate remains unknown. Here, we used the in vitro differentiation of human pluripotent stem cells into cardiomyocytes to study the role of c-JUN. Surprisingly, the knockout of c-JUN improved cardiomyocyte generation, as determined by the number of TNNT2+ cells. ATAC-seq data showed that the c-JUN defect led to increased chromatin accessibility on critical regulatory elements related to cardiomyocyte development. ChIP-seq data showed that the knockout c-JUN increased RBBP5 and SETD1B expression, leading to improved H3K4me3 deposition on key genes that regulate cardiogenesis. The c-JUN KO phenotype could be copied using the histone demethylase inhibitor CPI-455, which also up-regulated H3K4me3 levels and increased cardiomyocyte generation. Single-cell RNA-seq data defined three cell branches, and knockout c-JUN activated more regulons that are related to cardiogenesis. In summary, our data demonstrated that c-JUN could regulate cardiomyocyte cell fate by modulating H3K4me3 modification and chromatin accessibility and shed light on how c-JUN regulates heart development in humans.

Nuclear localization of TET2 requires β-catenin activation and correlates with favourable prognosis in colorectal cancer

Mutation-induced malfunction of ten-eleven translocation methylcytosine dioxygenase 2 (TET2) is widely reported in haematological malignancies. However, the role of TET2 in solid cancers, including colorectal cancer (CRC), is unclear. Here, we found that TET2 malfunction in CRC is mostly due to decreased nuclear localization and that nuclear localization of TET2 is correlated with better survival of patients. To explore the underlying mechanisms, 14 immortalized solid tumour cell lines and 12 primary CRC cell lines were used. TET2 was mostly detected in the nucleus, and it induced significant DNA demethylation and suppressed cell growth by demethylating RORA and SPARC in cell lines like SW480. While in cell lines like SW620, TET2 was observed in the cytosol and did not affect DNA methylation or cell growth. Further examination with immunoprecipitation-mass spectrometry illustrated that β-catenin activation was indispensable for the nuclear localization and tumour suppression effects of TET2. In addition, the β-catenin pathway activator IM12 and the TET2 activator vitamin C were used simultaneously to enhance the effects of TET2 under low-expression conditions, and synergistic inhibitory effects on the growth of cancer were observed both in vitro and in vivo. Collectively, these data suggest that β-catenin-mediated nuclear localization of TET2 is an important therapeutic target for solid tumours.

H3K36 methylation maintains cell identity by regulating opposing lineage programmes

The epigenetic mechanisms that maintain differentiated cell states remain incompletely understood. Here we employed histone mutants to uncover a crucial role for H3K36 methylation in the maintenance of cell identities across diverse developmental contexts. Focusing on the experimental induction of pluripotency, we show that H3K36M-mediated depletion of H3K36 methylation endows fibroblasts with a plastic state poised to acquire pluripotency in nearly all cells. At a cellular level, H3K36M facilitates epithelial plasticity by rendering fibroblasts insensitive to TGFβ signals. At a molecular level, H3K36M enables the decommissioning of mesenchymal enhancers and the parallel activation of epithelial/stem cell enhancers. This enhancer rewiring is Tet dependent and redirects Sox2 from promiscuous somatic to pluripotency targets. Our findings reveal a previously unappreciated dual role for H3K36 methylation in the maintenance of cell identity by integrating a crucial developmental pathway into sustained expression of cell-type-specific programmes, and by opposing the expression of alternative lineage programmes through enhancer methylation.

Restricting epigenetic activity promotes the reprogramming of transformed cells to pluripotency in a line-specific manner

Somatic cell reprogramming and oncogenic transformation share surprisingly similar features, yet transformed cells are resistant to reprogramming. Epigenetic barriers must block transformed cells from reprogramming, but the nature of those barriers is unclear. In this study, we generated a systematic panel of transformed mouse embryonic fibroblasts (MEFs) using oncogenic transgenes and discovered transformed cell lines compatible with reprogramming when transfected with Oct4/Sox2/Klf4/Myc. By comparing the reprogramming-capable and incapable transformed lines we identified multiple stages of failure in the reprogramming process. Some transformed lines failed at an early stage, whilst other lines seemed to progress through a conventional reprogramming process. Finally, we show that MEK inhibition overcomes one critical reprogramming barrier by indirectly suppressing a hyperacetylated active epigenetic state. This study reveals that diverse epigenetic barriers underly resistance to reprogramming of transformed cells.

The pluripotency factor Tex10 finetunes Wnt signaling for PGC and male germline development

Testis-specific transcript 10 (Tex10) is a critical factor for pluripotent stem cell maintenance and preimplantation development. Here, we dissect its late developmental roles in primordial germ cell (PGC) specification and spermatogenesis using cellular and animal models. We discover that Tex10 binds the Wnt negative regulator genes, marked by H3K4me3, at the PGC-like cell (PGCLC) stage in restraining Wnt signaling. Depletion and overexpression of Tex10 hyperactivate and attenuate the Wnt signaling, resulting in compromised and enhanced PGCLC specification efficiency, respectively. Using the Tex10 conditional knockout mouse models combined with single-cell RNA sequencing, we further uncover critical roles of Tex10 in spermatogenesis with Tex10 loss causing reduced sperm number and motility associated with compromised round spermatid formation. Notably, defective spermatogenesis in Tex10 knockout mice correlates with aberrant Wnt signaling upregulation. Therefore, our study establishes Tex10 as a previously unappreciated player in PGC specification and male germline development by fine-tuning Wnt signaling.

Comparative functional genomics identifies unique molecular features of EPSCs

The authors provide a comprehensive resource on proteomics, transcriptomic, and epigenetic level details of EPSCs to shed light on possible molecular pathways regulating their expanded pluripotency potential.

Extended pluripotent or expanded potential stem cells (EPSCs) possess superior developmental potential to embryonic stem cells (ESCs). However, the molecular underpinning of EPSC maintenance in vitro is not well defined. We comparatively studied transcriptome, chromatin accessibility, active histone modification marks, and relative proteomes of ESCs and the two well-established EPSC lines to probe the molecular foundation underlying EPSC developmental potential. Despite some overlapping transcriptomic and chromatin accessibility features, we defined sets of molecular signatures that distinguish EPSCs from ESCs in transcriptional and translational regulation as well as metabolic control. Interestingly, EPSCs show similar reliance on pluripotency factors Oct4, Sox2, and Nanog for self-renewal as ESCs. Our study provides a rich resource for dissecting the regulatory network that governs the developmental potency of EPSCs and exploring alternative strategies to capture totipotent stem cells in culture.

TNFAIP6 defines the MSC subpopulation with enhanced immune suppression activities

Background

Mesenchymal stromal/stem cells (MSCs) have been intensively investigated in both pre-clinical and clinical studies. However, the therapeutic efficacy varies resulting from the heterogenicity of MSCs. Therefore, purifying the specific MSC subpopulation with specialized function is necessary for their therapeutic applications.

Methods

The large-scale RNA sequencing analysis was performed to identify potential cell markers for the mouse MSCs. Then, the immune suppression activities of the purified MSC subpopulation were assessed in vitro and in vivo.

Results

The TNFAIP6 (tumor necrosis factor alpha-induced protein 6) has been identified as a potential cell marker for mouse MSCs, irrespective of tissue origin and laboratory origin. The TNFAIP6⁺ mouse MSCs showed enhanced immune suppression activities and improved therapeutic effects on the mouse model of acute inflammation, resulting from faster response to immune stimulation.

Conclusions

Therefore, we have demonstrated that the TNFAIP6⁺ MSC subpopulation has enhanced immune suppression capabilities.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-03176-5.

Metabolic and epigenetic dysfunctions underlie the arrest of in vitro fertilized human embryos in a senescent-like state

Around 60% of in vitro fertilized (IVF) human embryos irreversibly arrest before compaction between the 3- to 8-cell stage, posing a significant clinical problem. The mechanisms behind this arrest are unclear. Here, we show that the arrested embryos enter a senescent-like state, marked by cell cycle arrest, the down-regulation of ribosomes and histones and down-regulation of MYC and p53 activity. The arrested embryos can be divided into 3 types. Type I embryos fail to complete the maternal-zygotic transition, and Type II/III embryos have low levels of glycolysis and either high (Type II) or low (Type III) levels of oxidative phosphorylation. Treatment with the SIRT agonist resveratrol or nicotinamide riboside (NR) can partially rescue the arrested phenotype, which is accompanied by changes in metabolic activity. Overall, our data suggests metabolic and epigenetic dysfunctions underlie the arrest of human embryos.

eIF4A2 targets developmental potency and histone H3.3 transcripts for translational control of stem cell pluripotency

Translational control has emerged as a fundamental regulatory layer of proteome complexity that governs cellular identity and functions. As initiation is the rate-limiting step of translation, we carried out an RNA interference screen for key translation initiation factors required to maintain embryonic stem cell (ESC) identity. We identified eukaryotic translation initiation factor 4A2 (eIF4A2) and defined its mechanistic action through ribosomal protein S26-independent and -dependent ribosomes in translation initiation activation of messenger RNAs (mRNAs) encoding pluripotency factors and the histone variant H3.3 with demonstrated roles in maintaining stem cell pluripotency. eIF4A2 also mediates translation initiation activation of Ddx6, which acts together with eIF4A2 to restrict the totipotent two-cell transcription program in ESCs through Zscan4 mRNA degradation and translation repression. Accordingly, knockdown of eIF4A2 disrupts ESC proteome, causing the loss of ESC identity. Collectively, we establish a translational paradigm of the protein synthesis of pluripotency transcription factors and epigenetic regulators imposed on their established roles in controlling pluripotency.

MYOCD is Required for Cardiomyocyte-like Cells Induction from Human Urine Cells and Fibroblasts Through Remodeling Chromatin

Despite direct reprogramming of human cardiac fibroblasts into induced cardiomyocytes (iCM) holds great potential for heart regeneration, the mechanisms are poorly understood. Whether other human somatic cells could be reprogrammed into cardiomyocytes is also unknown. Here, we report human urine cells (hUCs) could be converted into CM-like cells from different donors and the related chromatin accessibility dynamics (CAD) by assay for transposase accessible chromatin(ATAC)-seq. hUCs transduced by MEF2C, TBX5, MESP1 and MYOCD but without GATA4 expressed multiple cardiac specific genes, exhibited Ca²⁺ oscillation potential and sarcomeric structures, and contracted synchronously in coculture with mouse CM. Additionally, we found that MYOCD is required for both closing and opening critical loci, mainly by hindering the opening of loci enriched with motifs for the TEAD and AP1 family and promoting the closing of loci enriched with ETS motifs. These changes differ partially from CAD observed during iCM induction from human fibroblasts. Collectively, our study offers one practical platform for iCM generation and insights into mechanisms for iCM fate determination.

Genome-wide annotation of protein-coding genes in pig

BACKGROUND

There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs.

RESULTS

An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented.

CONCLUSIONS

Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

Phenotype-tissue expression and exploration (PTEE) resource facilitates the choice of tissue for RNA-seq-based clinical genetics studies

BACKGROUND

RNA-seq emerges as a valuable method for clinical genetics. The transcriptome is "dynamic" and tissue-specific, but typically the probed tissues to analyze (TA) are different from the tissue of interest (TI) based on pathophysiology.

RESULTS

We developed Phenotype-Tissue Expression and Exploration (PTEE), a tool to facilitate the decision about the most suitable TA for RNA-seq. We integrated phenotype-annotated genes, used 54 tissues from GTEx to perform correlation analyses and identify expressed genes and transcripts between TAs and TIs. We identified skeletal muscle as the most appropriate TA to inquire for cardiac arrhythmia genes and skin as a good proxy to study neurodevelopmental disorders. We also explored RNA-seq limitations and show that on-off switching of gene expression during ontogenesis or circadian rhythm can cause blind spots for RNA-seq-based analyses.

CONCLUSIONS

PTEE aids the identification of tissues suitable for RNA-seq for a given pathology to increase the success rate of diagnosis and gene discovery. PTEE is freely available at https://bioinf.eva.mpg.de/PTEE/.

Identification of New Transcription Factors that Can Promote Pluripotent Reprogramming

BACKGROUND

Four transcription factors, Oct4, Sox2, Klf4, and c-Myc (the Yamanka factors), can reprogram somatic cells to induced pluripotent stem cells (iPSCs). Many studies have provided a number of alternative combinations to the non-Yamanaka factors. However, it is clear that many additional transcription factors that can generate iPSCs remain to be discovered.

METHODS

The chromatin accessibility and transcriptional level of human embryonic stem cells and human urine cells were compared by Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing (RNA-seq) to identify potential reprogramming factors. Selected transcription factors were employed to reprogram urine cells, and the reprogramming efficiency was measured. Urine-derived iPSCs were detected for pluripotency by Immunofluorescence, quantitative polymerase chain reaction, RNA sequencing and teratoma formation test. Finally, we assessed the differentiation potential of the new iPSCs to cardiomyocytes in vitro.

RESULTS

ATAC-seq and RNA-seq datasets predicted TEAD2, TEAD4 and ZIC3 as potential factors involved in urine cell reprogramming. Transfection of TEAD2, TEAD4 and ZIC3 (in the presence of Yamanaka factors) significantly improved the reprogramming efficiency of urine cells. We confirmed that the newly generated iPSCs possessed pluripotency characteristics similar to normal H1 embryonic stem cells. We also confirmed that the new iPSCs could differentiate to functional cardiomyocytes.

CONCLUSIONS

In conclusion, TEAD2, TEAD4 and ZIC3 can increase the efficiency of reprogramming human urine cells into iPSCs, and provides a new stem cell sources for the clinical application and modeling of cardiovascular disease.

Transposable element sequence fragments incorporated into coding and noncoding transcripts modulate the transcriptome of human pluripotent stem cells

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

Histone H2Bub1 deubiquitylation is essential for mouse development, but does not regulate global RNA polymerase II transcription

Co-activator complexes dynamically deposit post-translational modifications (PTMs) on histones, or remove them, to regulate chromatin accessibility and/or to create/erase docking surfaces for proteins that recognize histone PTMs. SAGA (Spt-Ada-Gcn5 Acetyltransferase) is an evolutionary conserved multisubunit co-activator complex with modular organization. The deubiquitylation module (DUB) of mammalian SAGA complex is composed of the ubiquitin-specific protease 22 (USP22) and three adaptor proteins, ATXN7, ATXN7L3 and ENY2, which are all needed for the full activity of the USP22 enzyme to remove monoubiquitin (ub1) from histone H2B. Two additional USP22-related ubiquitin hydrolases (called USP27X or USP51) have been described to form alternative DUBs with ATXN7L3 and ENY2, which can also deubiquitylate H2Bub1. Here we report that USP22 and ATXN7L3 are essential for normal embryonic development of mice, however their requirements are not identical during this process, as Atxn7l3^-/- embryos show developmental delay already at embryonic day (E) 7.5, while Usp22^-/- embryos are normal at this stage, but die at E14.5. Global histone H2Bub1 levels were only slightly affected in Usp22 null embryos, in contrast H2Bub1 levels were strongly increased in Atxn7l3 null embryos and derived cell lines. Our transcriptomic analyses carried out from wild type and Atxn7l3^-/- mouse embryonic stem cells (mESCs), or primary mouse embryonic fibroblasts (MEFs) suggest that the ATXN7L3-related DUB activity regulates only a subset of genes in both cell types. However, the gene sets and the extent of their deregulation were different in mESCs and MEFs. Interestingly, the strong increase of H2Bub1 levels observed in the Atxn7l3^-/- mESCs, or Atxn7l3^-/- MEFs, does not correlate with the modest changes in RNA Polymerase II (Pol II) occupancy and lack of changes in Pol II elongation observed in the two Atxn7l3^-/- cellular systems. These observations together indicate that deubiquitylation of histone H2Bub1 does not directly regulate global Pol II transcription elongation.

Prioritizing and characterizing functionally relevant genes across human tissues

Knowledge of genes that are critical to a tissue's function remains difficult to ascertain and presents a major bottleneck toward a mechanistic understanding of genotype-phenotype links. Here, we present the first machine learning model-FUGUE-combining transcriptional and network features, to predict tissue-relevant genes across 30 human tissues. FUGUE achieves an average cross-validation auROC of 0.86 and auPRC of 0.50 (expected 0.09). In independent datasets, FUGUE accurately distinguishes tissue or cell type-specific genes, significantly outperforming the conventional metric based on tissue-specific expression alone. Comparison of tissue-relevant transcription factors across tissue recapitulate their developmental relationships. Interestingly, the tissue-relevant genes cluster on the genome within topologically associated domains and furthermore, are highly enriched for differentially expressed genes in the corresponding cancer type. We provide the prioritized gene lists in 30 human tissues and an open-source software to prioritize genes in a novel context given multi-sample transcriptomic data.

MIMIC: an optimization method to identify cell type-specific marker panel for cell sorting

Multi-omics data allow us to select a small set of informative markers for the discrimination of specific cell types and study of cellular heterogeneity. However, it is often challenging to choose an optimal marker panel from the high-dimensional molecular profiles for a large amount of cell types. Here, we propose a method called Mixed Integer programming Model to Identify Cell type-specific marker panel (MIMIC). MIMIC maintains the hierarchical topology among different cell types and simultaneously maximizes the specificity of a fixed number of selected markers. MIMIC was benchmarked on the mouse ENCODE RNA-seq dataset, with 29 diverse tissues, for 43 surface markers (SMs) and 1345 transcription factors (TFs). MIMIC could select biologically meaningful markers and is robust for different accuracy criteria. It shows advantages over the standard single gene-based approaches and widely used dimensional reduction methods, such as multidimensional scaling and t-SNE, both in accuracy and in biological interpretation. Furthermore, the combination of SMs and TFs achieves better specificity than SMs or TFs alone. Applying MIMIC to a large collection of 641 RNA-seq samples covering 231 cell types identifies a panel of TFs and SMs that reveal the modularity of cell type association networks. Finally, the scalability of MIMIC is demonstrated by selecting enhancer markers from mouse ENCODE data. MIMIC is freely available at https://github.com/MengZou1/MIMIC.

Identifying transposable element expression dynamics and heterogeneity during development at the single-cell level with a processing pipeline scTE

Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity in unclear ways. Single-cell sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE expression has not been addressed. Here, we develop a single-cell TE processing pipeline, scTE, and report the expression of TEs in single cells in a range of biological contexts. Specific TE types are expressed in subpopulations of embryonic stem cells and are dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs are expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we apply scTE to single-cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results classify the dynamic patterns of TEs in single cells and their contributions to cell heterogeneity.

How transposable elements (TE) contribute to cell fate changes is unclear. Here, the authors generate a pipeline to quantify TE expression from single cell data. They show the dynamic expression of TEs from gastrulation to somatic cell reprogramming and human disease

Proteins in DNA methylation and their role in neural stem cell proliferation and differentiation

BACKGROUND

Epigenetic modifications, namely non-coding RNAs, DNA methylation, and histone modifications such as methylation, phosphorylation, acetylation, ubiquitylation, and sumoylation play a significant role in brain development. DNA methyltransferases, methyl-CpG binding proteins, and ten-eleven translocation proteins facilitate the maintenance, interpretation, and removal of DNA methylation, respectively. Different forms of methylation, including 5-methylcytosine, 5-hydroxymethylcytosine, and other oxidized forms, have been detected by recently developed sequencing technologies. Emerging evidence suggests that the diversity of DNA methylation patterns in the brain plays a key role in fine-tuning and coordinating gene expression in the development, plasticity, and disorders of the mammalian central nervous system. Neural stem cells (NSCs), originating from the neuroepithelium, generate neurons and glial cells in the central nervous system and contribute to brain plasticity in the adult mammalian brain.

MAIN BODY

Here, we summarized recent research in proteins responsible for the establishment, maintenance, interpretation, and removal of DNA methylation and those involved in the regulation of the proliferation and differentiation of NSCs. In addition, we discussed the interactions of chemicals with epigenetic pathways to regulate NSCs as well as the connections between proteins involved in DNA methylation and human diseases.

CONCLUSION

Understanding the interplay between DNA methylation and NSCs in a broad biological context can facilitate the related studies and reduce potential misunderstanding.

Chromatin and Epigenetic Rearrangements in Embryonic Stem Cell Fate Transitions

A major event in embryonic development is the rearrangement of epigenetic information as the somatic genome is reprogrammed for a new round of organismal development. Epigenetic data are held in chemical modifications on DNA and histones, and there are dramatic and dynamic changes in these marks during embryogenesis. However, the mechanisms behind this intricate process and how it is regulating and responding to embryonic development remain unclear. As embryos develop from totipotency to pluripotency, they pass through several distinct stages that can be captured permanently or transiently in vitro. Pluripotent naïve cells resemble the early epiblast, primed cells resemble the late epiblast, and blastomere-like cells have been isolated, although fully totipotent cells remain elusive. Experiments using these in vitro model systems have led to insights into chromatin changes in embryonic development, which has informed exploration of pre-implantation embryos. Intriguingly, human and mouse cells rely on different signaling and epigenetic pathways, and it remains a mystery why this variation exists. In this review, we will summarize the chromatin rearrangements in early embryonic development, drawing from genomic data from in vitro cell lines, and human and mouse embryos.

JMJD3 acts in tandem with KLF4 to facilitate reprogramming to pluripotency

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3’s catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription.

Previous work suggested that histone demethylase JMJD3 is detrimental to somatic cell reprogramming. Here, the authors show that while JMJD3 has a context-independent detrimental effect on early stages of reprogramming, during late stages it activates epithelial and pluripotency genes together with Klf4.

β-Catenin safeguards the ground state of mousepluripotency by strengthening the robustness of the transcriptional apparatus

Mouse embryonic stem cells cultured with MEK (mitogen-activated protein kinase kinase) and GSK3 (glycogen synthase kinase 3) inhibitors (2i) more closely resemble the inner cell mass of preimplantation blastocysts than those cultured with SL [serum/leukemia inhibitory factor (LIF)]. The transcriptional mechanisms governing this pluripotent ground state are unresolved. Release of promoter-proximal paused RNA polymerase II (Pol2) is a multistep process necessary for pluripotency and cell cycle gene transcription in SL. We show that β-catenin, stabilized by GSK3 inhibition in medium with 2i, supplies transcriptional coregulators at pluripotency loci. This selectively strengthens pluripotency loci and renders them addicted to transcription initiation for productive gene body elongation in detriment to Pol2 pause release. By contrast, cell cycle genes are not bound by β-catenin, and proliferation/self-renewal remains tightly controlled by Pol2 pause release under 2i conditions. Our findings explain how pluripotency is reinforced in the ground state and also provide a general model for transcriptional resilience/adaptation upon network perturbation in other contexts.

Metabolic switch and epithelial-mesenchymal transition cooperate to regulate pluripotency

Both metabolic switch from oxidative phosphorylation to glycolysis (OGS) and epithelial-mesenchymal transition (EMT) promote cellular reprogramming at early stages. However, their connections have not been elucidated. Here, when a chemically defined medium was used to induce early EMT during mouse reprogramming, a facilitated OGS was also observed at the same time. Additional investigations suggested that the two events formed a positive feedback loop via transcriptional activation, cooperated to upregulate epigenetic factors such as Bmi1, Ctcf, Ezh2, Kdm2b, and Wdr5, and accelerated pluripotency induction at the early stage. However, at late stages, by over-inducing glycolysis and preventing the necessary mesenchymal-epithelial transition, the two events trapped the cells at a new pluripotency state between naïve and primed states and inhibited further reprogramming toward the naïve state. In addition, the pluripotent stem cells at the new state have high similarity to epiblasts from E4.5 and E5.5 embryos, and have distinct characteristics from the previously reported epiblast-like or formative states. Therefore, the time-dependent cooperation between OGS and EMT in regulating pluripotency should extend our understanding of related fields.

DPre: computational identification of differentiation bias and genes underlying cell type conversions

SUMMARY

Cells are generally resistant to cell type conversions, but can be converted by the application of growth factors, chemical inhibitors and ectopic expression of genes. However, it remains difficult to accurately identify the destination cell type or differentiation bias when these techniques are used to alter cell type. Consequently, there is demand for computational techniques that can help researchers understand both the cell type and differentiation bias. While advanced tools identifying cell types exist for single cell data and the deconvolution of mixed cell populations, the problem of exploring partially differentiated cells of indeterminate transcriptional identity has not been addressed. To fill this gap, we developed driver-predictor, which relies on scoring per gene transcriptional similarity between RNA-Seq datasets to reveal directional bias of differentiation. By comparing against large cell type transcriptome libraries or a desired target expression profile, the tool enables the user to visualize both the changes in transcriptional identity as well as the genes accounting for the cell type changes. This software will be a powerful tool for researchers to explore in vitro experiments that involve cell type conversions.

AVAILABILITY AND IMPLEMENTATION

Source code is open source under the MIT license and is freely available on https://github.com/LoaloaF/DPre.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transgenerational Self-Reconstruction of Disrupted Chromatin Organization After Exposure To An Environmental Stressor in Mice

Exposure to environmental stressors is known to increase disease susceptibility in unexposed descendants in the absence of detectable genetic mutations. The mechanisms mediating environmentally-induced transgenerational disease susceptibility are poorly understood. We showed that great-great-grandsons of female mice exposed to tributyltin (TBT) throughout pregnancy and lactation were predisposed to obesity due to altered chromatin organization that subsequently biased DNA methylation and gene expression. Here we analyzed DNA methylomes and transcriptomes from tissues of animals ancestrally exposed to TBT spanning generations, sexes, ontogeny, and cell differentiation state. We found that TBT elicited concerted alterations in the expression of “chromatin organization” genes and inferred that TBT-disrupted chromatin organization might be able to self-reconstruct transgenerationally. We also found that the location of “chromatin organization” and “metabolic” genes is biased similarly in mouse and human genomes, suggesting that exposure to environmental stressors in different species could elicit similar phenotypic effects via self-reconstruction of disrupted chromatin organization.

Computational Methods for Mapping, Assembly and Quantification for Coding and Non-coding Transcripts

The measurement of gene expression has long provided significant insight into biological functions. The development of high-throughput short-read sequencing technology has revealed transcriptional complexity at an unprecedented scale, and informed almost all areas of biology. However, as researchers have sought to gather more insights from the data, these new technologies have also increased the computational analysis burden. In this review, we describe typical computational pipelines for RNA-Seq analysis and discuss their strengths and weaknesses for the assembly, quantification and analysis of coding and non-coding RNAs. We also discuss the assembly of transposable elements into transcripts, and the difficulty these repetitive elements pose. In summary, RNA-Seq is a powerful technology that is likely to remain a key asset in the biologist's toolkit.

Functional Landscape of PCGF Proteins Reveals Both RING1A/B-Dependent-and RING1A/B-Independent-Specific Activities

Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) control cell identity by establishing facultative heterochromatin repressive domains at common sets of target genes. PRC1, which deposits H2Aub1 through the E3 ligases RING1A/B, forms six biochemically distinct subcomplexes depending on the assembled PCGF protein (PCGF1–PCGF6); however, it is yet unclear whether these subcomplexes have also specific activities. Here we show that PCGF1 and PCGF2 largely compensate for each other, while other PCGF proteins have high levels of specificity for distinct target genes. PCGF2 associates with transcription repression, whereas PCGF3 and PCGF6 associate with actively transcribed genes. Notably, PCGF3 and PCGF6 complexes can assemble and be recruited to several active sites independently of RING1A/B activity (therefore, of PRC1). For chromatin recruitment, the PCGF6 complex requires the combinatorial activities of its MGA-MAX and E2F6-DP1 subunits, while PCGF3 requires an interaction with the USF1 DNA binding transcription factor.

•

PRC1 complexes retain high target specificity with little compensatory functions

•

PCGF2 defines repressive signatures and PCGF3 and PCGF6 transcription active states

•

RING1A/B is dispensable for the assembly and recruitment of PCGF3 and PCGF6 complexes

•

PCGF3 and PCGF6 complexes bind target genes through specific DNA binding activities

Single cell RNA-seq identifies the origins of heterogeneity in efficient cell transdifferentiation and reprogramming

Forced transcription factor expression can transdifferentiate somatic cells into other specialised cell types or reprogram them into induced pluripotent stem cells (iPSCs) with variable efficiency. To better understand the heterogeneity of these processes, we used single-cell RNA sequencing to follow the transdifferentation of murine pre-B cells into macrophages as well as their reprogramming into iPSCs. Even in these highly efficient systems, there was substantial variation in the speed and path of fate conversion. We predicted and validated that these differences are inversely coupled and arise in the starting cell population, with Myc^high large pre-BII cells transdifferentiating slowly but reprogramming efficiently and Myc^low small pre-BII cells transdifferentiating rapidly but failing to reprogram. Strikingly, differences in Myc activity predict the efficiency of reprogramming across a wide range of somatic cell types. These results illustrate how single cell expression and computational analyses can identify the origins of heterogeneity in cell fate conversion processes.

Hemi-methylated CpG sites connect Dnmt1-knockdown-induced and Tet1-induced DNA demethylation during somatic cell reprogramming

The relationship between active DNA demethylation induced by overexpressing Tet1 and passive DNA demethylation induced by suppressing Dnmt1 remains unclear. Here, we found that DNMT1 preferentially methylated, but TET1 preferentially demethylated, hemi-methylated CpG sites. These phenomena resulted in a significant overlap in the targets of these two types of DNA demethylation and the counteractions of Dnmt1 and Tet1 during somatic cell reprogramming. Since the hemi-methylated CpG sites generated during cell proliferation were enriched at core pluripotency loci, DNA demethylation induced by Tet1 or sh-RNA against Dnmt1 (sh-Dnmt1) was enriched in these loci, which, in combination with Yamanaka factors, led to the up-regulation of these genes and promoted somatic cell reprogramming. In addition, since sh-Dnmt1 induces DNA demethylation by impairing the further methylation of hemi-methylated CpG sites generated during cell proliferation, while Tet1 induced DNA demethylation by demethylating these hemi-methylated CpG sites, Tet1-induced DNA demethylation, compared with sh-Dnmt1-induced DNA demethylation, exhibited a higher ability to open the chromatin structure and up-regulate gene expression. Thus, Tet1-induced but not sh-Dnmt1-induced DNA demethylation led to the up-regulation of an additional set of genes that can promote the epithelial-mesenchymal transition and impair reprogramming. When vitamin C was used to further increase the demethylation ability of TET1 during reprogramming, Tet1 induced a larger up-regulation of these genes and significantly impaired reprogramming. Therefore, the current studies provide additional information regarding DNA demethylation during somatic cell reprogramming.

Transposable elements are regulated by context-specific patterns of chromatin marks in mouse embryonic stem cells

The majority of mammalian genomes are devoted to transposable elements (TEs). Whilst TEs are increasingly recognized for their important biological functions, they are a potential danger to genomic stability and are carefully regulated by the epigenetic system. However, the full complexity of this regulatory system is not understood. Here, using mouse embryonic stem cells, we show that TEs are suppressed by heterochromatic marks like H3K9me3, and are also labelled by all major types of chromatin modification in complex patterns, including bivalent activatory and repressive marks. We identified 29 epigenetic modifiers that significantly deregulated at least one type of TE. The loss of Setdb1, Ncor2, Rnf2, Kat5, Prmt5, Uhrf1, and Rrp8 caused widespread changes in TE expression and chromatin accessibility. These effects were context-specific, with different chromatin modifiers regulating the expression and chromatin accessibility of specific subsets of TEs. Our work reveals the complex patterns of epigenetic regulation of TEs.

Transposable elements (TEs) fulfill essential but poorly understood roles in genome organization and gene expression control. Here the authors show that the regulation of TEs occurs through overlapping epigenetic mechanisms that control the expression and chromatin signatures at TEs.

Chromatin Accessibility Dynamics during iPSC Reprogramming

Cell-fate decisions remain poorly understood at the chromatin level. Here, we map chromatin remodeling dynamics during induction of pluripotent stem cells. ATAC-seq profiling of MEFs expressing Oct4-Sox2-Klf4 (OSK) reveals dynamic changes in chromatin states shifting from open to closed (OC) and closed to open (CO), with an initial burst of OC and an ending surge of CO. The OC loci are largely composed of genes associated with a somatic fate, while the CO loci are associated with pluripotency. Factors/conditions known to impede reprogramming prevent OSK-driven OC and skew OC-CO dynamics. While the CO loci are enriched for OSK motifs, the OC loci are not, suggesting alternative mechanisms for chromatin closing. Sap30, a Sin3A corepressor complex component, is required for the OC shift and facilitates reduced H3K27ac deposition at OC loci. These results reveal a chromatin accessibility logic during reprogramming that may apply to other cell-fate decisions.

Diversity among POU transcription factors in chromatin recognition and cell fate reprogramming

The POU (Pit-Oct-Unc) protein family is an evolutionary ancient group of transcription factors (TFs) that bind specific DNA sequences to direct gene expression programs. The fundamental importance of POU TFs to orchestrate embryonic development and to direct cellular fate decisions is well established, but the molecular basis for this activity is insufficiently understood. POU TFs possess a bipartite 'two-in-one' DNA binding domain consisting of two independently folding structural units connected by a poorly conserved and flexible linker. Therefore, they represent a paradigmatic example to study the molecular basis for the functional versatility of TFs. Their modular architecture endows POU TFs with the capacity to accommodate alternative composite DNA sequences by adopting different quaternary structures. Moreover, associations with partner proteins crucially influence the selection of their DNA binding sites. The plentitude of DNA binding modes confers the ability to POU TFs to regulate distinct genes in the context of different cellular environments. Likewise, different binding modes of POU proteins to DNA could trigger alternative regulatory responses in the context of different genomic locations of the same cell. Prominent POU TFs such as Oct4, Brn2, Oct6 and Brn4 are not only essential regulators of development but have also been successfully employed to reprogram somatic cells to pluripotency and neural lineages. Here we review biochemical, structural, genomic and cellular reprogramming studies to examine how the ability of POU TFs to select regulatory DNA, alone or with partner factors, is tied to their capacity to epigenetically remodel chromatin and drive specific regulatory programs that give cells their identities.

Genomic and molecular control of cell type and cell type conversions

Organisms are made of a limited number of cell types that combine to form higher order tissues and organs. Cell types have traditionally been defined by their morphologies or biological activity, yet the underlying molecular controls of cell type remain unclear. The onset of single cell technologies, and more recently genomics (particularly single cell genomics), has substantially increased the understanding of the concept of cell type, but has also increased the complexity of this understanding. These new technologies have added a new genome wide molecular dimension to the description of cell type, with genome-wide expression and epigenetic data acting as a cell type ‘fingerprint’ to describe the cell state. Using these genomic fingerprints cell types are being increasingly defined based on specific genomic and molecular criteria, without necessarily a distinct biological function. In this review, we will discuss the molecular definitions of cell types and cell type control, and particularly how endogenous and exogenous transcription factors can control cell types and cell type conversions.

Coop-Seq Analysis Demonstrates that Sox2 Evokes Latent Specificities in the DNA Recognition by Pax6

Sox2 and Pax6 co-regulate genes in neural lineages and the lens by forming a ternary complex likely facilitated allosterically through DNA. We used the quantitative and scalable cooperativity-by-sequencing (Coop-seq) approach to interrogate Sox2/Pax6 dimerization on a DNA library where five positions of the Pax6 half-site were randomized yielding 1024 cooperativity factors. Consensus positions normally required for the high-affinity DNA binding by Pax6 need to be mutated for effective dimerization with Sox2. Out of the five randomized bases, a 5' thymidine is present in most of the top ranking elements. However, this thymidine maps to a region outside of the Pax half site and is not expected to directly interact with Pax6 in known binding modes suggesting structural reconfigurations. Re-analysis of ChIP-seq data identified several genomic regions where the cooperativity promoting sequence pattern is co-bound by Sox2 and Pax6. A highly conserved Sox2/Pax6 bound site near the Sprouty2 locus was verified to promote cooperative dimerization designating Sprouty2 as a potential target reliant on Sox2/Pax6 cooperativity in several neural cell types. Collectively, the functional interplay of Sox2 and Pax6 demands the relaxation of high-affinity binding sites and is enabled by alternative DNA sequences. We conclude that this binding mode evolved to warrant that a subset of target genes is only regulated in the presence of suitable partner factors.