Ground rules of the pluripotency gene regulatory network

Single-cell profiling identifies LIN28A mRNA targets in the mouse pluripotent-to-2C-like transition and somatic cell reprogramming

RNA-binding proteins (RBPs) regulate totipotency, pluripotency maintenance, and induction. The intricacies of how they modulate these processes through their interaction with RNAs remain to be elucidated. Here we employed Targets of RBPs Identified By Editing (TRIBE) with single-cell resolution (scTRIBE) to profile the mRNA targets of the key pluripotency regulator LIN28A in mouse embryonic stem cells (ESCs), 2-cell embryo-like cells (2CLCs) and somatic cell reprogramming. LIN28A is known to act by controlling the maturation of the let-7 microRNA but, in addition, it binds to multiple mRNAs and influences their stability and translation efficiency. However, the mRNA targets of LIN28A in 2CLCs and reprogramming are unclear. Through quantitative single-cell analysis of the scTRIBE dataset, we observed a marked increase in the binding of LIN28A to mRNAs of ribosome biogenesis factors and a selected group of totipotency factors in 2CLCs within ESC cultures. Our results suggest that LIN28A extends the half-life of at least some of these mRNAs, providing new insights into its role in the totipotent state. We also uncovered the distinct trajectory-specific LIN28A-mRNA networks in reprogramming, helping explain how LIN28A facilitates the mesenchymal-to-epithelial transition and pluripotency acquisition. Our study not only clarifies the multifunctional role of LIN28A in these processes but also highlights the importance of decoding RNA-protein interactions at the single-cell level.

LMNA-Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-Derived iPSC Differentiation Support Cell Type and Lineage-Specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA-related dilated cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C (LMNA) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. The molecular mechanisms of the disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA-related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA-mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (four from Patients and eight from Controls) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for cardiac progenitors to cardiomyocytes (CMs) and epicardium-derived cells (EPDCs). Data integration and comparative analyses of Patient and Control cells found cell type and lineage-specific differentially expressed genes (DEGs) with enrichment, supporting pathway dysregulation. Top DEGs and enriched pathways included 10 ZNF genes and RNA polymerase II transcription in pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CMs; LMNA and epigenetic regulation, as well as DDIT4 and mTORC1 signaling in EPDCs. Top DEGs also included XIST and other X-linked genes, six imprinted genes (SNRPN, PWAR6, NDN, PEG10, MEG3, MEG8), and enriched gene sets related to metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs, as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

Poised PABP-RNA hubs implement signal-dependent mRNA decay in development

Signaling pathways drive cell fate transitions largely by changing gene expression. However, the mechanisms for rapid and selective transcriptome rewiring in response to signaling cues remain elusive. Here we use deep learning to deconvolve both the sequence determinants and the trans-acting regulators that trigger extracellular signal-regulated kinase (ERK)-mitogen-activated protein kinase kinase (MEK)-induced decay of the naive pluripotency mRNAs. Timing of decay is coupled to embryo implantation through ERK-MEK phosphorylation of LIN28A, which repositions pLIN28A to the highly A+U-rich 3' untranslated region (3'UTR) termini of naive pluripotency mRNAs. Interestingly, these A+U-rich 3'UTR termini serve as poly(A)-binding protein (PABP)-binding hubs, poised for signal-induced convergence with LIN28A. The multivalency of AUU motifs determines the efficacy of pLIN28A-PABP convergence, which enhances PABP 3'UTR binding, decreases the protection of poly(A) tails and activates mRNA decay to enable progression toward primed pluripotency. Thus, the signal-induced convergence of LIN28A with PABP-RNA hubs drives the rapid selection of naive mRNAs for decay, enabling the transcriptome remodeling that ensures swift developmental progression.

Gene regulatory network inference with covariance dynamics

Determining gene regulatory network (GRN) structure is a central problem in biology, with a variety of inference methods available for different types of data. For a widely prevalent and challenging use case, namely single-cell gene expression data measured after intervention at multiple time points with unknown joint distributions, there is only one known specifically developed method, which does not fully utilize the rich information contained in this data type. We develop an inference method for the GRN in this case, netWork infErence by covariaNce DYnamics, dubbed WENDY. The core idea of WENDY is to model the dynamics of the covariance matrix, and solve this dynamics as an optimization problem to determine the regulatory relationships. To evaluate its effectiveness, we compare WENDY with other inference methods using synthetic data and experimental data. Our results demonstrate that WENDY performs well across different data sets.

Canalizing kernel for cell fate determination

The tendency for cell fate to be robust to most perturbations, yet sensitive to certain perturbations raises intriguing questions about the existence of a key path within the underlying molecular network that critically determines distinct cell fates. Reprogramming and trans-differentiation clearly show examples of cell fate change by regulating only a few or even a single molecular switch. However, it is still unknown how to identify such a switch, called a master regulator, and how cell fate is determined by its regulation. Here, we present CAESAR, a computational framework that can systematically identify master regulators and unravel the resulting canalizing kernel, a key substructure of interconnected feedbacks that is critical for cell fate determination. We demonstrate that CAESAR can successfully predict reprogramming factors for de-differentiation into mouse embryonic stem cells and trans-differentiation of hematopoietic stem cells, while unveiling the underlying essential mechanism through the canalizing kernel. CAESAR provides a system-level understanding of how complex molecular networks determine cell fates.

LMNA -Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-derived iPSC Differentiation Support Cell type and Lineage-specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA -Related Dilated Cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C ( LMNA ) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. Molecular mechanisms of disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA -Related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (4 Patient and 8 Control) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for Cardiac Progenitors to Cardiomyocytes (CM) and Epicardium-Derived Cells (EPDC). Data integration and comparative analyses of Patient and Control cells found cell type and lineage differentially expressed genes (DEG) with enrichment to support pathway dysregulation. Top DEG and enriched pathways included: 10 ZNF genes and RNA polymerase II transcription in Pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CM; LMNA and epigenetic regulation and DDIT4 and mTORC1 signaling in EPDC. Top DEG also included: XIST and other X-linked genes, six imprinted genes: SNRPN , PWAR6 , NDN , PEG10 , MEG3 , MEG8 , and enriched gene sets in metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

A Novel Approach to Comparative RNA-Seq Does Not Support a Conserved Set of Orthologs Underlying Animal Regeneration

Molecular studies of animal regeneration typically focus on conserved genes and signaling pathways that underlie morphogenesis. To date, a holistic analysis of gene expression across animals has not been attempted, as it presents a suite of problems related to differences in experimental design and gene homology. By combining orthology analyses with a novel statistical method for testing gene enrichment across large data sets, we are able to test whether tissue regeneration across animals shares transcriptional regulation. We applied this method to a meta-analysis of six publicly available RNA-Seq data sets from diverse examples of animal regeneration. We recovered 160 conserved orthologous gene clusters, which are enriched in structural genes as opposed to those regulating morphogenesis. A breakdown of gene presence/absence provides limited support for the conservation of pathways typically implicated in regeneration, such as Wnt signaling and cell pluripotency pathways. Such pathways are only conserved if we permit large amounts of paralog switching through evolution. Overall, our analysis does not support the hypothesis that a shared set of ancestral genes underlie regeneration mechanisms in animals. After applying the same method to heat shock studies and getting similar results, we raise broader questions about the ability of comparative RNA-Seq to reveal conserved gene pathways across deep evolutionary relationships.

Single-cell 3D genome structure reveals distinct human pluripotent states

BACKGROUND

Pluripotent states of embryonic stem cells (ESCs) with distinct transcriptional profiles affect ESC differentiative capacity and therapeutic potential. Although single-cell RNA sequencing has revealed additional subpopulations and specific features of naive and primed human pluripotent stem cells (hPSCs), the underlying mechanisms that regulate their specific transcription and that control their pluripotent states remain elusive.

RESULTS

By single-cell analysis of high-resolution, three-dimensional (3D) genomic structure, we herein demonstrate that remodeling of genomic structure is highly associated with the pluripotent states of human ESCs (hESCs). The naive pluripotent state is featured with specialized 3D genomic structures and clear chromatin compartmentalization that is distinct from the primed state. The naive pluripotent state is achieved by remodeling the active euchromatin compartment and reducing chromatin interactions at the nuclear center. This unique genomic organization is linked to enhanced chromatin accessibility on enhancers and elevated expression levels of naive pluripotent genes localized to this region. In contradistinction, the primed state exhibits intermingled genomic organization. Moreover, active euchromatin and primed pluripotent genes are distributed at the nuclear periphery, while repressive heterochromatin is densely concentrated at the nuclear center, reducing chromatin accessibility and the transcription of naive genes.

CONCLUSIONS

Our data provide insights into the chromatin structure of ESCs in their naive and primed states, and we identify specific patterns of modifications in transcription and chromatin structure that might explain the genes that are differentially expressed between naive and primed hESCs. Thus, the inversion or relocation of heterochromatin to euchromatin via compartmentalization is related to the regulation of chromatin accessibility, thereby defining pluripotent states and cellular identity.

Multifaceted effects on even-skipped transcriptional dynamics upon Krüppel dosage changes

Although fluctuations in transcription factor (TF) dosage are often well tolerated, TF dosage modulation can change the target gene expression dynamics and result in significant non-lethal developmental phenotypes. Using MS2/MCP-mediated quantitative live imaging in early Drosophila embryos, we analyzed how changing levels of the gap gene Krüppel (Kr) affects transcriptional dynamics of the pair-rule gene even-skipped (eve). Halving the Kr dosage leads to a transient posterior expansion of the eve stripe 2 and an anterior shift of stripe 5. Surprisingly, the most significant changes are observed in eve stripes 3 and 4, the enhancers of which do not contain Kr-binding sites. In Kr heterozygous embryos, both stripes 3 and 4 display narrower widths, anteriorly shifted boundaries and reduced mRNA production levels. We show that Kr dosage indirectly affects stripe 3 and 4 dynamics by modulating other gap gene dynamics. We quantitatively correlate moderate body segment phenotypes of Kr heterozygotes with spatiotemporal changes in eve expression. Our results indicate that nonlinear relationships between TF dosage and phenotypes underlie direct TF-DNA and indirect TF-TF interactions.

Complex regulatory networks influence pluripotent cell state transitions in human iPSCs

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

Transcription factor co-expression mediates lineage priming for embryonic and extra-embryonic differentiation

In early mammalian development, cleavage stage blastomeres and inner cell mass (ICM) cells co-express embryonic and extra-embryonic transcriptional determinants. Using a protein-based double reporter we identify an embryonic stem cell (ESC) population that co-expresses the extra-embryonic factor GATA6 alongside the embryonic factor SOX2. Based on single cell transcriptomics, we find this population resembles the unsegregated ICM, exhibiting enhanced differentiation potential for endoderm while maintaining epiblast competence. To relate transcription factor binding in these cells to future fate, we describe a complete enhancer set in both ESCs and naive extra-embryonic endoderm stem cells and assess SOX2 and GATA6 binding at these elements in the ICM-like ESC sub-population. Both factors support cooperative recognition in these lineages, with GATA6 bound alongside SOX2 on a fraction of pluripotency enhancers and SOX2 alongside GATA6 more extensively on endoderm enhancers, suggesting that cooperative binding between these antagonistic factors both supports self-renewal and prepares progenitor cells for later differentiation.

Ubiquitination-specific protease 7 enhances stemness of hepatocellular carcinoma by stabilizing basic transcription factor 3

This study aims to explore the molecular regulation mechanism of ubiquitination-specific protease 7 (USP7) in facilitating the stemness properties of hepatocellular carcinoma (HCC). Gain-of-function and loss-of-function assays were conducted in SK-Hep1 and HepG2 cells transfected with USP7 overexpression/knockdown plasmids and USP7 inhibitor P22077. The proliferation, migration, invasion, and self-renewal capacity of hepatocellular carcinoma cells were detected by CCK-8, colony formation, Transwell, scratch, and tumor sphere formation, respectively. MS was performed to identify the potential substrate of USP7 following P22077 treatment. Co-IP assay was used to verify the interaction between USP7 and basic transcription factor 3 (BTF3) in HCC cells. The overexpression of USP7 could promote the proliferation, migration, invasion, and colony formation capacity of SK-Hep1 and HepG2 cells. Additionally, ectopic UPS7 enhanced the epithelial-mesenchymal transition (EMT) and stem-like characteristics of the HCC cells. In contrast, USP7 depletion by knockdown of USP7 or administrating inhibitor P22077 significantly inhibited these malignant phenotypes of SK-Hep1 and HepG2 cells. Following MS analysis, BTF3 was identified as a potential substrate for USP7. USP7 could interact with BTF3 and upregulate its protein level, while USP7 depletion significantly upregulated the ubiquitination levels. Overexpression of BTF3 partially rescue the inhibitory effects of USP7 depletion on the malignant phenotypes and stemness properties of SK-Hep1 and HepG2 cells. USP7 can promote the stemness and malignant phenotype of HCC by stabilizing BTF3.

Boldine promotes stemness of human urine-derived stem cells by activating the Wnt/β-catenin signaling pathway

Human urine-derived stem cells (hUSCs) process self-renewal and multilineage differentiation ability. Due to their non-invasive and easily available clinical source, hUSCs represent a promising alternative source of mesenchymal stem cells (MSCs) for application potential in cytotherapy. However, technical limitations, such as stemness property maintenance, have hindered hUSCs' clinical application. Certain some small molecules have been recognized with advantage in maintaining the stemness of stem cells. In this study, we identified stemness-regulated key targets of hUSCs based on the StemCellNet database, CMAP database and literature mining. Furthermore, we identified a small molecule compound, boldine, which may have the potential to promote the stemness of hUSCs. It promotes cell proliferation, multilineage differentiation and maintains stemness of hUSCs by cell viability assay, single-cell clone formation, osteogenic differentiation and stemness marker expression (OCT-4 and C-MYC). We identified that boldine may be a potential GSK-3β inhibitor by molecular docking and confirmed that it can upregulate the level of β-catenin and promote translocation of β-catenin into nucleus of hUSCs using Western blotting and immunofluorescence analysis. Our study indicates boldine activates the Wnt/β-catenin signaling pathway in hUSCs and provides an effective strategy for MSCs research and application of small molecules in maintaining the stemness of hUSCs.

RENGE infers gene regulatory networks using time-series single-cell RNA-seq data with CRISPR perturbations

Single-cell RNA-seq analysis coupled with CRISPR-based perturbation has enabled the inference of gene regulatory networks with causal relationships. However, a snapshot of single-cell CRISPR data may not lead to an accurate inference, since a gene knockout can influence multi-layered downstream over time. Here, we developed RENGE, a computational method that infers gene regulatory networks using a time-series single-cell CRISPR dataset. RENGE models the propagation process of the effects elicited by a gene knockout on its regulatory network. It can distinguish between direct and indirect regulations, which allows for the inference of regulations by genes that are not knocked out. RENGE therefore outperforms current methods in the accuracy of inferring gene regulatory networks. When used on a dataset we derived from human-induced pluripotent stem cells, RENGE yielded a network consistent with multiple databases and literature. Accurate inference of gene regulatory networks by RENGE would enable the identification of key factors for various biological systems.

Reactivation of embryonic genetic programs in tissue regeneration and disease

Embryonic genetic programs are reactivated in response to various types of tissue damage, providing cell plasticity for tissue regeneration or disease progression. In acute conditions, these programs remedy the damage and then halt to allow a return to homeostasis. In chronic situations, including inflammatory diseases, fibrosis and cancer, prolonged activation of embryonic programs leads to disease progression and tissue deterioration. Induction of progenitor identity and cell plasticity, for example, epithelial-mesenchymal plasticity, are critical outcomes of reactivated embryonic programs. In this Review, we describe molecular players governing reactivated embryonic genetic programs, their role during disease progression, their similarities and differences and lineage reversion in pathology and discuss associated therapeutics and drug-resistance mechanisms across many organs. We also discuss the diversity of reactivated programs in different disease contexts. A comprehensive overview of commonalities between development and disease will provide better understanding of the biology and therapeutic strategies.

METTL3-Dependent N6-Methyladenosine Modification Programs Human Neural Progenitor Cell Proliferation

METTL3, a methyltransferase responsible for N6-methyladenosine (m6A) modification, plays key regulatory roles in mammal central neural system (CNS) development. However, the specific epigenetic mechanisms governing human CNS development remain poorly elucidated. Here, we generated small-molecule-assisted shut-off (SMASh)-tagged hESC lines to reduce METTL3 protein levels, and found that METTL3 is not required for human neural progenitor cell (hNPC) formation and neuron differentiation. However, METTL3 deficiency inhibited hNPC proliferation by reducing SLIT2 expression. Mechanistic studies revealed that METTL3 degradation in hNPCs significantly decreased the enrichment of m6A in SLIT2 mRNA, consequently reducing its expression. Our findings reveal a novel functional target (SLIT2) for METTL3 in hNPCs and contribute to a better understanding of m6A-dependent mechanisms in hNPC proliferation.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Identification of the regulatory circuit governing corneal epithelial fate determination and disease

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity.

Ascl1 and Ngn2 convert mouse embryonic stem cells to neurons via functionally distinct paths

Ascl1 and Ngn2, closely related proneural transcription factors, are able to convert mouse embryonic stem cells into induced neurons. Despite their similarities, these factors elicit only partially overlapping transcriptional programs, and it remains unknown whether cells are converted via distinct mechanisms. Here we show that Ascl1 and Ngn2 induce mutually exclusive side populations by binding and activating distinct lineage drivers. Furthermore, Ascl1 rapidly dismantles the pluripotency network and installs neuronal and trophoblast cell fates, while Ngn2 generates a neural stem cell-like intermediate supported by incomplete shutdown of the pluripotency network. Using CRISPR-Cas9 knockout screening, we find that Ascl1 relies more on factors regulating pluripotency and the cell cycle, such as Tcf7l1. In the absence of Tcf7l1, Ascl1 still represses core pluripotency genes but fails to exit the cell cycle. However, overexpression of Cdkn1c induces cell cycle exit and restores the generation of neurons. These findings highlight that cell type conversion can occur through two distinct mechanistic paths, even when induced by closely related transcription factors.

MPP8 Governs the Activity of the LIF/STAT3 Pathway and Plays a Crucial Role in the Differentiation of Mouse Embryonic Stem Cells

Mouse embryonic stem cells (mESCs) possess the remarkable characteristics of unlimited self-renewal and pluripotency, which render them highly valuable for both fundamental research and clinical applications. A comprehensive understanding of the molecular mechanisms underlying mESC function is of the utmost importance. The Human Silence Hub (HUSH) complex, comprising FAM208A, MPP8, and periphilin, constitutes an epigenetic silencing complex involved in suppressing retroviruses and transposons during early embryonic development. However, its precise role in regulating mESC pluripotency and differentiation remains elusive. In this study, we generated homogenous miniIAA7-tagged Mpp8 mouse ES cell lines. Upon induction of MPP8 protein degradation, we observed the impaired proliferation and reduced colony formation ability of mESCs. Furthermore, this study unveils the involvement of MPP8 in regulating the activity of the LIF/STAT3 signaling pathway and Nanog expression in mESCs. Finally, we provide compelling evidence that degradation of the MPP8 protein impairs the differentiation of mESC.

USP7 represses lineage differentiation genes in mouse embryonic stem cells by both catalytic and noncatalytic activities

USP7, a ubiquitin-specific peptidase (USP), plays an important role in many cellular processes through its catalytic deubiquitination of various substrates. However, its nuclear function that shapes the transcriptional network in mouse embryonic stem cells (mESCs) remains poorly understood. We report that USP7 maintains mESC identity through both catalytic activity-dependent and -independent repression of lineage differentiation genes. Usp7 depletion attenuates SOX2 levels and derepresses lineage differentiation genes thereby compromising mESC pluripotency. Mechanistically, USP7 deubiquitinates and stabilizes SOX2 to repress mesoendodermal (ME) lineage genes. Moreover, USP7 assembles into RYBP-variant Polycomb repressive complex 1 and contributes to Polycomb chromatin-mediated repression of ME lineage genes in a catalytic activity-dependent manner. USP7 deficiency in its deubiquitination function is able to maintain RYBP binding to chromatin for repressing primitive endoderm-associated genes. Our study demonstrates that USP7 harbors both catalytic and noncatalytic activities to repress different lineage differentiation genes, thereby revealing a previously unrecognized role in controlling gene expression for maintaining mESC identity.

A multi-omics integrative analysis based on CRISPR screens re-defines the pluripotency regulatory network in ESCs

A comprehensive and precise definition of the pluripotency gene regulatory network (PGRN) is crucial for clarifying the regulatory mechanisms in embryonic stem cells (ESCs). Here, after a CRISPR/Cas9-based functional genomics screen and integrative analysis with other functional genomes, transcriptomes, proteomes and epigenome data, an expanded pluripotency-associated gene set is obtained, and a new PGRN with nine sub-classes is constructed. By integrating the DNA binding, epigenetic modification, chromatin conformation, and RNA expression profiles, the PGRN is resolved to six functionally independent transcriptional modules (CORE, MYC, PAF, PRC, PCGF and TBX). Spatiotemporal transcriptomics reveal activated CORE/MYC/PAF module activity and repressed PRC/PCGF/TBX module activity in both mouse ESCs (mESCs) and pluripotent cells of early embryos. Moreover, this module activity pattern is found to be shared by human ESCs (hESCs) and cancers. Thus, our results provide novel insights into elucidating the molecular basis of ESC pluripotency.

Long range inter-chromosomal interaction of Oct4 distal enhancer loci regulates ESCs pluripotency

Nuclear architecture underlies the transcriptional programs within the cell to establish cell identity. As previously demonstrated, long-range chromatin interactions of the Oct4 distal enhancer (DE) are correlated with active transcription in naïve state embryonic stem cells. Here, we identify and characterize extreme long-range interactions of the Oct4 DE through a novel CRISPR labeling technique we developed and chromosome conformation capture to identify lethal giant larvae 2 (Llgl2) and growth factor receptor-bound protein 7 (Grb7) as putative functional interacting target genes in different chromosomes. We show that the Oct4 DE directly regulates expression of Llgl2 and Grb7 in addition to Oct4. Expression of Llgl2 and Grb7 closely correlates with the pluripotent state, where knock down of either result in loss of pluripotency, and overexpression enhances somatic cell reprogramming. We demonstrated that biologically important interactions of the Oct4 DE can occur at extreme distances that are necessary for the maintenance of the pluripotent state.

Role of the Paf1 complex in the maintenance of stem cell pluripotency and development

Cell identity is determined by the transcriptional regulation of a cell-type-specific gene group. The Paf1 complex (Paf1C), an RNA polymerase II-associating factor, is an important transcriptional regulator that not only participates in transcription elongation and termination but also affects transcription-coupled histone modifications and chromatin organisation. Recent studies have shown that Paf1C is involved in the expression of genes required for self-renewal and pluripotency in stem cells and tumorigenesis. In this review, we focused on the role of Paf1C as a critical transcriptional regulator in cell fate decisions. Paf1C affects the pluripotency of stem cells by regulating the expression of core transcription factors such as Oct4 and Nanog. In addition, Paf1C directly binds to the promoters or distant elements of target genes, thereby maintaining the pluripotency in embryonic stem cells derived from an early stage of the mammalian embryo. Paf1C is upregulated in cancer stem cells, as compared with that in cancer cells, suggesting that Paf1C may be a target for cancer therapy. Interestingly, Paf1C is involved in multiple developmental stages in Drosophila, zebrafish, mice and even humans, thereby displaying a trend for the correlation between Paf1C and cell fate. Thus, we propose that Paf1C is a critical contributor to cell differentiation, cell specification and its characteristics and could be employed as a therapeutic target in developmental diseases.

Control of RNA degradation in cell fate decision

Cell fate is shaped by a unique gene expression program, which reflects the concerted action of multilayered precise regulation. Substantial research attention has been paid to the contribution of RNA biogenesis to cell fate decisions. However, increasing evidence shows that RNA degradation, well known for its function in RNA processing and the surveillance of aberrant transcripts, is broadly engaged in cell fate decisions, such as maternal-to-zygotic transition (MZT), stem cell differentiation, or somatic cell reprogramming. In this review, we first look at the diverse RNA degradation pathways in the cytoplasm and nucleus. Then, we summarize how selective transcript clearance is regulated and integrated into the gene expression regulation network for the establishment, maintenance, and exit from a special cellular state.

The role of Wnt signaling in the development of the epiblast and axial progenitors

Understanding how the body plan is established during embryogenesis remains a fundamental biological question. The Wnt/β-catenin signaling pathway plays a crucial and highly conserved role in body plan formation, functioning to polarize the primary anterior-posterior (AP) or head-to-tail body axis in most metazoans. In this chapter, we focus on the roles that the mammalian Wnt/β-catenin pathway plays to prepare the pluripotent epiblast for gastrulation, and to elicit the emergence of multipotent axial progenitors from the caudal epiblast. Interactions between Wnt and retinoic acid (RA), another powerful family of developmental signaling molecules, in axial progenitors will also be discussed. Gastrulation movements and somitogenesis result in the anterior displacement of the RA source (the rostral somites and lateral plate mesoderm (LPM)), from the posterior Wnt source (the primitive streak (PS)), leading to the establishment of antiparallel gradients of RA and Wnt that control the self-renewal and successive differentiation of neck, trunk and tail progenitors.

RNA degradation eliminates developmental transcripts during murine embryonic stem cell differentiation via CAPRIN1-XRN2

Embryonic stem cells (ESCs) are self-renewing and pluripotent. In recent years, factors that control pluripotency, mostly nuclear, have been identified. To identify non-nuclear regulators of ESCs, we screened an endogenously labeled fluorescent fusion-protein library in mouse ESCs. One of the more compelling hits was the cell-cycle-associated protein 1 (CAPRIN1). CAPRIN1 knockout had little effect in ESCs, but it significantly altered differentiation and gene expression programs. Using RIP-seq and SLAM-seq, we found that CAPRIN1 associates with, and promotes the degradation of, thousands of RNA transcripts. CAPRIN1 interactome identified XRN2 as the likely ribonuclease. Upon early ESC differentiation, XRN2 is located in the nucleus and colocalizes with CAPRIN1 in small RNA granules in a CAPRIN1-dependent manner. We propose that CAPRIN1 regulates an RNA degradation pathway operating during early ESC differentiation, thus eliminating undesired spuriously transcribed transcripts in ESCs.

Nucleosome breathing facilitates cooperative binding of pluripotency factors Sox2 and Oct4 to DNA

Critical lineage commitment events are staged by multiple transcription factors (TFs) binding to their cognate motifs, often positioned at nucleosome-enriched regions of chromatin. The underlying mechanism remains elusive due to difficulty in disentangling the heterogeneity in chromatin states. Using a novel coarse-grained model and molecular dynamics simulations, here we probe the association of Sox2 and Oct4 proteins that show clustered binding at the entry-exit region of a nucleosome. The model captures the conformational heterogeneity of nucleosome breathing dynamics that features repeated wrap-unwrap transitions of a DNA segment from one end of the nucleosome. During the dynamics, DNA forms bulges that diffuse stochastically and may regulate the target search dynamics of a protein by nonspecifically interacting with it. The overall search kinetics of the TF pair follows a "dissociation-compensated-association" mechanism, where Oct4 binding is facilitated by the association of Sox2. The cooperativity stems from a change in entropy caused by an alteration in the nucleosome dynamics upon TF binding. The binding pattern is consistent with a live-cell single-particle tracking experiment, suggesting the mechanism observed for clustered binding of a TF pair, which is a hallmark of cis-regulatory elements, has broader implications in understanding gene regulation in a complex chromatin environment.

KLF17 promotes human naive pluripotency through repressing MAPK3 and ZIC2

The pluripotent state of embryonic stem cells (ESCs) is regulated by a sophisticated network of transcription factors. High expression of KLF17 has recently been identified as a hallmark of naive state of human ESCs (hESCs). However, the functional role of KLF17 in naive state is not clear. Here, by employing various gain and loss-of-function approaches, we demonstrate that KLF17 is essential for the maintenance of naive state and promotes the primed to naive state transition in hESCs. Mechanistically, we identify MAPK3 and ZIC2 as two direct targets repressed by KLF17. Overexpression of MAPK3 or ZIC2 partially blocks KLF17 from promoting the naive pluripotency. Furthermore, we find that human and mouse homologs of KLF17 retain conserved functions in promoting naive pluripotency of both species. Finally, we show that Klf17 may be essential for early embryo development in mouse. These findings demonstrate the important and conserved function of KLF17 in promoting naive pluripotency and reveal two essential transcriptional targets of KLF17 that underlie its function.

Oct4:Sox2 binding is essential for establishing but not maintaining active and silent states of dynamically regulated genes in pluripotent cells

Much has been learned about the mechanisms of action of pluripotency factors Oct4 and Sox2. However, as with other regulators of cell identity, little is known about the impact of disrupting their binding motifs in a native environment or the characteristics of genes they regulate. By quantitatively examining dynamic ranges of gene expression instead of focusing on conventional measures of differential expression, we found that Oct4 and Sox2 enhancer binding is strongly enriched near genes subject to large dynamic ranges of expression among cell types, with binding sites near these genes usually within superenhancers. Mutagenesis of representative Oct4:Sox2 motifs near such active, dynamically regulated genes revealed critical roles in transcriptional activation during reprogramming, with more limited roles in transcriptional maintenance in the pluripotent state. Furthermore, representative motifs near silent genes were critical for establishing but not maintaining the fully silent state, while genes whose transcript levels varied by smaller magnitudes among cell types were unaffected by nearby Oct4:Sox2 motifs. These results suggest that Oct4 and Sox2 directly establish both active and silent transcriptional states in pluripotent cells at a large number of genes subject to dynamic regulation during mammalian development, but are less important than expected for maintaining transcriptional states.

Dynamic regulation of P-TEFb by 7SK snRNP is integral to the DNA damage response to regulate chemotherapy sensitivity

Testicular germ cell tumors and closely related embryonal stem cells are exquisitely sensitive to cisplatin, a feature thought to be linked to their pluripotent state and p53 status. It remains unclear whether and how cellular state is coordinated with p53 to confer cisplatin sensitivity. Here, we report that positive transcription elongation factor b (P-TEFb) determines cell fate upon DNA damage. We find that cisplatin rapidly activates P-TEFb by releasing it from inhibitory 7SK small nuclear ribonucleoprotein complex. P-TEFb directly phosphorylates pluripotency factor estrogen-related receptor beta (ESRRB), and induces its proteasomal degradation to enhance pro-survival glycolysis. On the other hand, P-TEFb is required for the transcription of a substantial portion of p53 target genes, triggering cell death during prolonged cisplatin treatment. These results reveal previously underappreciated roles of P-TEFb to coordinate the DNA damage response. We discuss the implications for using P-TEFb inhibitors to treat cancer and ameliorate cisplatin-induced ototoxicity.

•

P-TEFb regulates pro-survival and pro-death pathways during DNA damage response

•

P-TEFb promotes ESRRB proteasomal degradation to enhance pro-survival glycolysis

•

P-TEFb induces a substantial portion of p53 target genes to trigger cell death

•

Chemical inhibitors of P-TEFb blocks cisplatin- or UV-induced cell death

New developments in the biology of fibroblast growth factors

The fibroblast growth factor (FGF) family is composed of 18 secreted signaling proteins consisting of canonical FGFs and endocrine FGFs that activate four receptor tyrosine kinases (FGFRs 1-4) and four intracellular proteins (intracellular FGFs or iFGFs) that primarily function to regulate the activity of voltage-gated sodium channels and other molecules. The canonical FGFs, endocrine FGFs, and iFGFs have been reviewed extensively by us and others. In this review, we briefly summarize past reviews and then focus on new developments in the FGF field since our last review in 2015. Some of the highlights in the past 6 years include the use of optogenetic tools, viral vectors, and inducible transgenes to experimentally modulate FGF signaling, the clinical use of small molecule FGFR inhibitors, an expanded understanding of endocrine FGF signaling, functions for FGF signaling in stem cell pluripotency and differentiation, roles for FGF signaling in tissue homeostasis and regeneration, a continuing elaboration of mechanisms of FGF signaling in development, and an expanding appreciation of roles for FGF signaling in neuropsychiatric diseases. This article is categorized under: Cardiovascular Diseases > Molecular and Cellular Physiology Neurological Diseases > Molecular and Cellular Physiology Congenital Diseases > Stem Cells and Development Cancer > Stem Cells and Development.

Transcription factor HESX1 enhances mesendodermal commitment of human embryonic stem cells by modulating ERK1/2 signaling

Transcription factors are key determinants of lineage commitment during mammalian development. However, the function and molecular mechanism for the transcription factors in the formation of three primary germ layers during human embryonic development are not fully elucidated. Here, we report that homeobox-containing transcription factor HESX1 plays a critical role in mesendodermal (ME) commitment of human embryonic stem cells (hESCs). Our results show that expression of HESX1 in hESCs is regulated by OCT4 and NANOG, and that its expression level changes with hESC differentiation. We find that knockdown of HESX1 does not disrupt the undifferentiated state of hESCs, in terms of cell morphology and expression levels of pluripotency-associated genes. However, HESX1 deficiency in hESCs impairs their ME commitment, whereas forced expression of HESX1 significantly enhances ME marker expression during ME commitment. Interestingly, HESX1 knockdown in hESCs represses ERK1/2 signaling activated by ME induction, while overexpression of HESX1 markedly enhances ERK1/2 activity during ME commitment of hESCs. Of note, MEK inhibitor PD0325901 weakens or even eliminates HESX1 overexpression-mediated promotive effects on ME induction in a dosage-dependent manner. Together, this study identifies a novel role of HESX1 in hESC commitment to ME cells and establishes the functional link between a transcription factor and lineage-associated signaling. These findings would help to better understand early human development and develop more efficient protocols to induce hESC differentiation to desired lineages.

Cnot8 eliminates naïve regulation networks and is essential for naïve-to-formative pluripotency transition

Mammalian early epiblasts at different phases are characterized by naïve, formative, and primed pluripotency states, involving extensive transcriptome changes. Here, we report that deadenylase Cnot8 of Ccr4-Not complex plays essential roles during the transition from naïve to formative state. Knock out (KO) Cnot8 resulted in early embryonic lethality in mice, but Cnot8 KO embryonic stem cells (ESCs) could be established. Compared with the cells differentiated from normal ESCs, Cnot8 KO cells highly expressed a great many genes during their differentiation into the formative state, including several hundred naïve-like genes enriched in lipid metabolic process and gene expression regulation that may form the naïve regulation networks. Knockdown expression of the selected genes of naïve regulation networks partially rescued the differentiation defects of Cnot8 KO ESCs. Cnot8 depletion led to the deadenylation defects of its targets, increasing their poly(A) tail lengths and half-life, eventually elevating their expression levels. We further found that Cnot8 was involved in the clearance of targets through its deadenylase activity and the binding of Ccr4-Not complex, as well as the interacting with Tob1 and Pabpc1. Our results suggest that Cnot8 eliminates naïve regulation networks through mRNA clearance, and is essential for naïve-to-formative pluripotency transition.

Research Progress of Totipotent Stem Cells

Totipotent stem cells (TSCs), can develop into complete organisms, used in biological fields such as regenerative medicine, mammalian breeding, and conservation. However, cells from early-stage embryos cultured are hard to self-renew and maintain developmental totipotency, which becomes a key factor limiting the research of TSCs. Fortunately, a break-through in the study of induced pluripotent stem cells returning to their totipotent state has been made, resulting in the establishment of multiple TSCs and igniting a new wave of stem cell research. Furthermore, the blastocyst-like structures can be generated by the established TSCs, which lays a foundation for synthetic embryos in vitro. In this review, we summarize the totipotent stage of the early embryos, the establishment and cultivation of TSCs, and the developmental ability exploration of TSCs to promote further research of TSCs.

Generation and gene expression profiles of Grevy's zebra induced Pluripotent Stem Cells

Induced pluripotent stem cells (iPSCs) can serve as a biological resource for functional and conservation research for various species. This realisation has led to the generation of iPSCs from many species, including those identified as endangered. However, the understanding of species variation in mammalian iPSCs remains largely unknown. To gain insight into species variation in iPSCs, we generated iPSCs from a new species Grevy's zebra (Equus grevyi; gz-iPSCs), which has been listed as endangered in the IUCN (International Union for Conservation of Nature) Red List. We isolated primary fibroblast cells from an individual and successfully reprogrammed them into iPSCs. The generated gz-iPSCs continued to grow under primed-type culture condition and showed pluripotency and differentiation potential. To describe the molecular characteristics of gz-iPSCs, we performed RNA sequencing analysis. The gz-iPSC transcriptome showed robust expression of pluripotency associated genes reported in human and mouse, suggesting evolutionary conservation among the species. This study provides insight into the iPSCs from a rare species and helps the understanding of the gene expression basis underlying mammalian PSCs.

A balanced Oct4 interactome is crucial for maintaining pluripotency

Oct4 collaborates primarily with other transcriptional factors or coregulators to maintain pluripotency. However, how Oct4 exerts its function is still unclear. Here, we show that the Oct4 linker interface mediates competing yet balanced Oct4 protein interactions that are crucial for maintaining pluripotency. Oct4 linker mutant embryonic stem cells (ESCs) show decreased expression of self-renewal genes and increased expression of differentiation genes, resulting in impaired ESC self-renewal and early embryonic development. The linker mutation interrupts the balanced Oct4 interactome. In mutant ESCs, the interaction between Oct4 and Klf5 is decreased. In contrast, interactions between Oct4 and Cbx1, Ctr9, and Cdc73 are increased, disrupting the epigenetic state of ESCs. Control of the expression level of Klf5, Cbx1, or Cdc73 rebalances the Oct4 interactome and rescues the pluripotency of linker mutant ESCs, indicating that such factors interact with Oct4 competitively. Thus, we provide previously unidentified molecular insights into how Oct4 maintains pluripotency.

Defining the Pluripotent Marker Genes for Identification of Teleost Fish Cell Pluripotency During Reprogramming

Pluripotency is a transient state in early embryos, which is regulated by an interconnected network of pluripotency-related genes. The pluripotent state itself seems to be highly dynamic, which leads to significant differences in the description of induced pluripotent stem cells from different species at the molecular level. With the application of cell reprogramming technology in fish, the establishment of a set of molecular standards for defining pluripotency will be important for the research and potential application of induced pluripotent stem cells in fish. In this study, by BLAST search and expression pattern analysis, we screen out four pluripotent genes (Oct4, Nanog, Tdgf1, and Gdf3) in zebrafish (Danio rerio) and crucian carp (Carassius). These genes were highly expressed in the short period of early embryonic development, but significantly down-regulated after differentiation. Moreover, three genes (Oct4, Nanog and Tdgf1) have been verified that are suitable for identifying the pluripotency of induced pluripotent stem cells in zebrafish and crucian carp. Our study expands the understanding of the pluripotent markers of induced pluripotent stem cells in fish.

The spatial self-organization within pluripotent stem cell colonies is continued in detaching aggregates

Colonies of induced pluripotent stem cells (iPSCs) reveal aspects of self-organization even under culture conditions that maintain pluripotency. To investigate the dynamics of this process under spatial confinement, we used either polydimethylsiloxane (PDMS) pillars or micro-contact printing of vitronectin. There was a progressive upregulation of OCT4, E-cadherin, and NANOG within 70 μm from the outer rim of iPSC colonies. Single-cell RNA-sequencing and spatial reconstruction of gene expression demonstrated that OCT4^high subsets, residing at the edge of the colony, have pronounced up-regulation of the TGF-β pathway, particularly of NODAL and its inhibitor LEFTY. Interestingly, after 5-7 days, iPSC colonies detached spontaneously from micro-contact printed substrates to form 3D aggregates. This new method allowed generation of embryoid bodies (EBs) of controlled size without enzymatic or mechanical treatment. Within the early 3D aggregates, radial organization and differential gene expression continued in analogy to the changes observed during self-organization of iPSC colonies. Early self-detached aggregates revealed up-regulated germline-specific gene expression patterns as compared to conventional EBs. However, there were no marked differences after further directed differentiation toward hematopoietic, mesenchymal, and neuronal lineages. Our results provide further insight into the gradual self-organization within iPSC colonies and at their transition into EBs.

Periodic synchronization of isolated network elements facilitates simulating and inferring gene regulatory networks including stochastic molecular kinetics

Background

The temporal progression of many fundamental processes in cells and organisms, including homeostasis, differentiation and development, are governed by gene regulatory networks (GRNs). GRNs balance fluctuations in the output of their genes, which trace back to the stochasticity of molecular interactions. Although highly desirable to understand life processes, predicting the temporal progression of gene products within a GRN is challenging when considering stochastic events such as transcription factor–DNA interactions or protein production and degradation.

Results

We report a method to simulate and infer GRNs including genes and biochemical reactions at molecular detail. In our approach, we consider each network element to be isolated from other elements during small time intervals, after which we synchronize molecule numbers across all network elements. Thereby, the temporal behaviour of network elements is decoupled and can be treated by local stochastic or deterministic solutions. We demonstrate the working principle of this modular approach with a repressive gene cascade comprising four genes. By considering a deterministic time evolution within each time interval for all elements, our method approaches the solution of the system of deterministic differential equations associated with the GRN. By allowing genes to stochastically switch between on and off states or by considering stochastic production of gene outputs, we are able to include increasing levels of stochastic detail and approximate the solution of a Gillespie simulation. Thereby, CaiNet is able to reproduce noise-induced bi-stability and oscillations in dynamically complex GRNs. Notably, our modular approach further allows for a simple consideration of deterministic delays. We further infer relevant regulatory connections and steady-state parameters of a GRN of up to ten genes from steady-state measurements by identifying each gene of the network with a single perceptron in an artificial neuronal network and using a gradient decent method originally designed to train recurrent neural networks. To facilitate setting up GRNs and using our simulation and inference method, we provide a fast computer-aided interactive network simulation environment, CaiNet.

Conclusion

We developed a method to simulate GRNs at molecular detail and to infer the topology and steady-state parameters of GRNs. Our method and associated user-friendly framework CaiNet should prove helpful to analyze or predict the temporal progression of reaction networks or GRNs in cellular and organismic biology. CaiNet is freely available at https://gitlab.com/GebhardtLab/CaiNet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04541-6.

Pluripotency transcription factors at the focus: the phase separation paradigm in stem cells

The transcription factors (TFs) OCT4, SOX2 and NANOG are key players of the gene regulatory network of pluripotent stem cells. Evidence accumulated in recent years shows that even small imbalances in the expression levels or relative concentrations of these TFs affect both, the maintenance of pluripotency and cell fate decisions. In addition, many components of the transcriptional machinery including RNA polymerases, cofactors and TFs such as those required for pluripotency, do not distribute homogeneously in the nucleus but concentrate in multiple foci influencing the delivery of these molecules to their DNA-targets. How cells control strict levels of available pluripotency TFs in this heterogeneous space and the biological role of these foci remain elusive. In recent years, a wealth of evidence led to propose that many of the nuclear compartments are formed through a liquid-liquid phase separation process. This new paradigm early penetrated the stem cells field since many key players of the pluripotency circuitry seem to phase-separate. Overall, the formation of liquid compartments may modulate the kinetics of biochemical reactions and consequently regulate many nuclear processes. Here, we review the state-of-the-art knowledge of compartmentalization in the cell nucleus and the relevance of this process for transcriptional regulation, particularly in pluripotent stem cells. We also highlight the recent advances and new ideas in the field showing how compartmentalization may affect pluripotency preservation and cell fate decisions.

Cardiac specification during gastrulation - The Yellow Brick Road leading to Tinman

The question of how the heart develops, and the genetic networks governing this process have become intense areas of research over the past several decades. This research is propelled by classical developmental studies and potential clinical applications to understand and treat congenital conditions in which cardiac development is disrupted. Discovery of the tinman gene in Drosophila, and examination of its vertebrate homolog Nkx2.5, along with other core cardiac transcription factors has revealed how cardiac progenitor differentiation and maturation drives heart development. Careful observation of cardiac morphogenesis along with lineage tracing approaches indicated that cardiac progenitors can be divided into two broad classes of cells, namely the first and second heart fields, that contribute to the heart in two distinct waves of differentiation. Ample evidence suggests that the fate of individual cardiac progenitors is restricted to distinct cardiac structures quite early in development, well before the expression of canonical cardiac progenitor markers like Nkx2.5. Here we review the initial specification of cardiac progenitors, discuss evidence for the early patterning of cardiac progenitors during gastrulation, and consider how early gene expression programs and epigenetic patterns can direct their development. A complete understanding of when and how the developmental potential of cardiac progenitors is determined, and their potential plasticity, is of great interest developmentally and also has important implications for both the study of congenital heart disease and therapeutic approaches based on cardiac stem cell programming.

SOX transcription factors and glioma stem cells: Choosing between stemness and differentiation

Glioblastoma (GBM) is the most common, most aggressive and deadliest brain tumor. Recently, remarkable progress has been made towards understanding the cellular and molecular biology of gliomas. GBM tumor initiation, progression and relapse as well as resistance to treatments are associated with glioma stem cells (GSCs). GSCs exhibit a high proliferation rate and self-renewal capacity and the ability to differentiate into diverse cell types, generating a range of distinct cell types within the tumor, leading to cellular heterogeneity. GBM tumors may contain different subsets of GSCs, and some of them may adopt a quiescent state that protects them against chemotherapy and radiotherapy. GSCs enriched in recurrent gliomas acquire more aggressive and therapy-resistant properties, making them more malignant, able to rapidly spread. The impact of SOX transcription factors (TFs) on brain tumors has been extensively studied in the last decade. Almost all SOX genes are expressed in GBM, and their expression levels are associated with patient prognosis and survival. Numerous SOX TFs are involved in the maintenance of the stemness of GSCs or play a role in the initiation of GSC differentiation. The fine-tuning of SOX gene expression levels controls the balance between cell stemness and differentiation. Therefore, innovative therapies targeting SOX TFs are emerging as promising tools for combatting GBM. Combatting GBM has been a demanding and challenging goal for decades. The current therapeutic strategies have not yet provided a cure for GBM and have only resulted in a slight improvement in patient survival. Novel approaches will require the fine adjustment of multimodal therapeutic strategies that simultaneously target numerous hallmarks of cancer cells to win the battle against GBM.

noisyR: enhancing biological signal in sequencing datasets by characterizing random technical noise

High-throughput sequencing enables an unprecedented resolution in transcript quantification, at the cost of magnifying the impact of technical noise. The consistent reduction of random background noise to capture functionally meaningful biological signals is still challenging. Intrinsic sequencing variability introducing low-level expression variations can obscure patterns in downstream analyses. We introduce noisyR, a comprehensive noise filter to assess the variation in signal distribution and achieve an optimal information-consistency across replicates and samples; this selection also facilitates meaningful pattern recognition outside the background-noise range. noisyR is applicable to count matrices and sequencing data; it outputs sample-specific signal/noise thresholds and filtered expression matrices. We exemplify the effects of minimizing technical noise on several datasets, across various sequencing assays: coding, non-coding RNAs and interactions, at bulk and single-cell level. An immediate consequence of filtering out noise is the convergence of predictions (differential-expression calls, enrichment analyses and inference of gene regulatory networks) across different approaches.

microRNA regulation of pluripotent state transition

microRNAs (miRNAs) play essential roles in mouse embryonic stem cells (ESCs) and early embryo development. The exact mechanism by which miRNAs regulate cell fate transition during embryo development is still not clear. Recent studies have identified and captured various pluripotent stem cells (PSCs) that share similar characteristics with cells from different stages of pre- and post-implantation embryos. These PSCs provide valuable models to understand miRNA functions in early mammalian development. In this short review, we will summarize recent work towards understanding the function and mechanism of miRNAs in regulating the transition or conversion between different pluripotent states. In addition, we will highlight unresolved questions and key future directions related to miRNAs in pluripotent state transition. Studies in these areas will further our understanding of miRNA functions in early embryo development, and may lead to practical means to control human PSCs for clinical applications in regenerative medicine.

Pluripotency-State-Dependent Role of Dax1 in Embryonic Stem Cells Self-Renewal

Dax1(also known as Nr0b1) is regarded as an important component of the transcription factor network in mouse embryonic stem cells (ESCs). However, the role and the molecular mechanism of Dax1 in the maintenance of different pluripotency states are poorly understood. Here, we constructed a stable Dax1 knockout (KO) cell line using the CRISPR/Cas9 system to analyze the precise function of Dax1. We reported that 2i/LIF-ESCs had significantly lower Dax1 expression than LIF/serum-ESCs. Dax1KO ES cell lines could be established in 2i/LIF and their pluripotency was confirmed. In contrast, Dax1-null ESCs could not be continuously passaged in LIF/serum due to severe differentiation and apoptosis. In LIF/serum, the activities of the Core module and Myc module were significantly reduced, while the PRC2 module was activated after Dax1KO. The expression of most proapoptotic genes and lineage-commitment genes were drastically increased, while the downregulated expression of antiapoptotic genes and many pluripotency genes was observed. Our research on the pluripotent state-dependent role of Dax1 provides clues to understand the molecular regulation mechanism at different stages of early embryonic development.

ZFP451-mediated SUMOylation of SATB2 drives embryonic stem cell differentiation

Here, Urrutia et al. set out to study the mechanism that regulates the choice between pluripotency and differentiation in embryonic stem cells (ESCs). Using biochemical and genomic analyses, the authors identify SUMO2 modification of Satb2 by the E3 ligase Zfp451 as a driver of ESC differentiation.

The establishment of cell fates involves alterations of transcription factor repertoires and repurposing of transcription factors by post-translational modifications. In embryonic stem cells (ESCs), the chromatin organizers SATB2 and SATB1 balance pluripotency and differentiation by activating and repressing pluripotency genes, respectively. Here, we show that conditional Satb2 gene inactivation weakens ESC pluripotency, and we identify SUMO2 modification of SATB2 by the E3 ligase ZFP451 as a potential driver of ESC differentiation. Mutations of two SUMO-acceptor lysines of Satb2 (Satb2^{K →R}) or knockout of Zfp451 impair the ability of ESCs to silence pluripotency genes and activate differentiation-associated genes in response to retinoic acid (RA) treatment. Notably, the forced expression of a SUMO2-SATB2 fusion protein in either Satb2^{K →R} or Zfp451^−/− ESCs rescues, in part, their impaired differentiation potential and enhances the down-regulation of Nanog. The differentiation defect of Satb2^{K →R} ESCs correlates with altered higher-order chromatin interactions relative to Satb2^wt ESCs. Upon RA treatment of Satb2^wt ESCs, SATB2 interacts with ZFP451 and the LSD1/CoREST complex and gains binding at differentiation genes, which is not observed in RA-treated Satb2^{K →R} cells. Thus, SATB2 SUMOylation may contribute to the rewiring of transcriptional networks and the chromatin interactome of ESCs in the transition of pluripotency to differentiation.

Continuous expression of reprogramming factors induces and maintains mouse pluripotency without specific growth factors and signaling inhibitors

Objectives

Derivation and maintenance of pluripotent stem cells (PSCs) generally require optimized and complex culture media, which hinders the derivation of PSCs from various species. Expression of Oct4, Sox2, Klf4, and c‐Myc (OSKM) can reprogram somatic cells into induced PSCs (iPSCs), even for species possessing no optimal culture condition. Herein, we explored whether expression of OSKM could induce and maintain pluripotency without PSC‐specific growth factors and signaling inhibitors.

Methods

The culture medium of Tet‐On‐OSKM/Oct4‐GFP mouse embryonic stem cells (ESCs) was switched from N2B27 with MEK inhibitor, GSK3β inhibitor, and leukemia inhibitory factor (LIF) (2iL) to N2B27 with doxycycline. Tet‐On‐OSKM mouse embryonic fibroblast (MEF) cells were reprogrammed in N2B27 with doxycycline. Cell proliferation was traced. Pluripotency was assessed by expression of ESC marker genes, teratoma, and chimera formation. RNA‐Seq was conducted to analyze gene expression.

Results

Via continuous expression of OSKM, mouse ESCs (OSKM‐ESCs) and the resulting iPSCs (OSKM‐iPSCs) reprogrammed from MEF cells propagated stably, expressed pluripotency marker genes, and formed three germ layers in teratomas. Transcriptional landscapes of OSKM‐iPSCs resembled those of ESCs cultured in 2iL and were more similar to those of ESCs cultured in serum/LIF. Furthermore, OSKM‐iPSCs contributed to germline transmission.

Conclusions

Expression of OSKM could induce and maintain mouse pluripotency without specific culturing factors. Importantly, OSKM‐iPSCs could produce gene‐modified animals through germline transmission, with potential applications in other species.

Direct Reprogramming of Cardiac Fibroblasts to Repair the Injured Heart

Coronary heart disease is a leading cause of mortality and morbidity. Those that survive acute myocardial infarction are at significant risk of subsequent heart failure due to fibrotic remodelling of the infarcted myocardium. By applying knowledge from the study of embryonic cardiovascular development, modern medicine offers hope for treatment of this condition through regeneration of the myocardium by direct reprogramming of fibrotic scar tissue. Here, we will review mechanisms of cell fate specification leading to the generation of cardiovascular cell types in the embryo and use this as a framework in which to understand direct reprogramming. Driving expression of a network of transcription factors, micro RNA or small molecule epigenetic modifiers can reverse epigenetic silencing, reverting differentiated cells to a state of induced pluripotency. The pluripotent state can be bypassed by direct reprogramming in which one differentiated cell type can be transdifferentiated into another. Transdifferentiating cardiac fibroblasts to cardiomyocytes requires a network of transcription factors similar to that observed in embryonic multipotent cardiac progenitors. There is some flexibility in the composition of this network. These studies raise the possibility that the failing heart could one day be regenerated by directly reprogramming cardiac fibroblasts within post-infarct scar tissue.

Uncovering the RNA-binding protein landscape in the pluripotency network of human embryonic stem cells

Embryonic stem cell (ESC) self-renewal and cell fate decisions are driven by a broad array of molecular signals. While transcriptional regulators have been extensively studied in human ESCs (hESCs), the extent to which RNA-binding proteins (RBPs) contribute to human pluripotency remains unclear. Here, we carry out a proteome-wide screen and identify 810 proteins that bind RNA in hESCs. We reveal that RBPs are preferentially expressed in hESCs and dynamically regulated during early stem cell differentiation. Notably, many RBPs are affected by knockdown of OCT4, a master regulator of pluripotency, several dozen of which are directly targeted by this factor. Using cross-linking and immunoprecipitation (CLIP-seq), we find that the pluripotency-associated STAT3 and OCT4 transcription factors interact with RNA in hESCs and confirm the binding of STAT3 to the conserved NORAD long-noncoding RNA. Our findings indicate that RBPs have a more widespread role in human pluripotency than previously appreciated.