Mouse Embryonic Stem Cells Research Articles

Capture of Totipotency in Mouse Embryonic Stem Cells in the Absence of Pdzk1.

Totipotent cells can differentiate into three lineages: the epiblast, primitive endoderm, and trophectoderm. Naturally, only early fertilized embryos possess totipotency, and they lose this ability as they develop. The expansion of stem cell differentiation potential has been a hot topic in developmental biology for years, particularly with respect to the generation totipotent-like stem cells. Here, the study describes the establishment of totipotency in embryonic stem cells (ESCs) via the deletion of a single gene, Pdzk1. Pdzk1-knockout (KO) ESCs substantially contribute to the fetus, placenta, and yolk sac in chimera assays but can also self-organize to form standard blastocyst-like structures containing the three lineages efficiently; thus, they exhibit full developmental potential as early blastomeres. Single-cell transcriptome and bulk RNA-seq comprehensive analyses revealed that Pdzk1-KO activates several lineage inducers (C1qa, C1qb, Fgf5, and Cdx2) to break down barriers between embryonic and extraembryonic tissues, making these lineages switch smoothly and resulting in a totipotent-like state. This versatile and scalable system provides a robust experimental model for differentiation potency and cell fate studies.

Long-range regulation of transcription scales with genomic distance in a gene-specific manner.

Although critical for tuning the timing and level of transcription, enhancer communication with distal promoters is not well understood. Here, we bypass the need for sequence-specific transcription factors (TFs) and recruit activators directly using a chimeric array of gRNA oligos to target dCas9 fused to the activator VP64-p65-Rta (CARGO-VPR). We show that this approach achieves effective activator recruitment to arbitrary genomic sites, even those inaccessible when targeted with a single guide. We utilize CARGO-VPR across the Prdm8-Fgf5 locus in mouse embryonic stem cells (mESCs), where neither gene is expressed. Although activator recruitment to any tested region results in the transcriptional induction of at least one gene, the expression level strongly depends on the genomic distance between the promoter and activator recruitment site. However, the expression-distance relationship for each gene scales distinctly in a manner not attributable to differences in 3D contact frequency, promoter DNA sequence, or the presence of repressive chromatin marks at the locus.

Enhancer cooperativity can compensate for loss of activity over large genomic distances.

Enhancers are short DNA sequences that activate their target promoter from a distance; however, increasing the genomic distance between the enhancer and the promoter decreases expression levels. Many genes are controlled by combinations of multiple enhancers, yet the interaction and cooperation of individual enhancer elements are not well understood. Here, we developed a synthetic platform in mouse embryonic stem cells that allows building complex regulatory landscapes from the bottom up. We tested the system by integrating individual enhancers at different distances and confirmed that the strength of an enhancer contributes to how strongly it is affected by increased genomic distance. Furthermore, synergy between two enhancer elements depends on the distance at which the two elements are integrated: introducing a weak enhancer between a strong enhancer and the promoter strongly increases reporter gene expression, allowing enhancers to activate from increased genomic distances.

DEAD-box RNA helicase 10 is required for 18S rRNA maturation by controlling the release of U3 snoRNA from pre-rRNA in embryonic stem cells

Ribosome biogenesis plays a pivotal role in maintaining stem cell homeostasis, yet the precise regulatory mechanisms governing this process in mouse embryonic stem cells (mESCs) remain largely unknown. In this investigation, we ascertain that DEAD-box RNA helicase 10 (DDX10) is indispensable for upholding cellular homeostasis and the viability of mESCs. Positioned predominantly at the nucleolar dense fibrillar component (DFC) and granular component (GC), DDX10 predominantly binds to 45S ribosomal RNA (rRNA) and orchestrates ribosome biogenesis. Degradation of DDX10 prevents the release of U3 snoRNA from pre-rRNA, leading to perturbed pre-rRNA processing and compromised maturation of the 18S rRNA, thereby disrupting the biogenesis of the small ribosomal subunit. Moreover, DDX10 participates in the process of liquid-liquid phase separation (LLPS), which is necessary for efficient ribosome biogenesis. Notably, the NUP98-DDX10 fusion associated with acute myelocytic leukemia (AML) alters the cellular localization of DDX10 and results in loss of ability to regulate pre-rRNA processing. Collectively, this study reveals the critical role of DDX10 as a pivotal regulator of ribosome biogenesis in mESCs.

The interplay of DNA repair context with target sequence predictably biases Cas9-generated mutations

Repair of double-stranded breaks generated by CRISPR/Cas9 is highly dependent on the flanking DNA sequence. To learn about interactions between DNA repair and target sequence, we measure frequencies of over 236,000 distinct Cas9-generated mutational outcomes at over 2800 synthetic target sequences in 18 DNA repair deficient mouse embryonic stem cells lines. We classify the outcomes in an unbiased way, finding a specialised role for Prkdc (DNA-PKcs protein) and Polm in creating 1 bp insertions matching the nucleotide on the protospacer-adjacent motif side of the break, a variable involvement of Nbn and Polq in the creation of different deletion outcomes, and uni-directional deletions dependent on both end-protection and end-resection. Using our dataset, we build predictive models of the mutagenic outcomes of Cas9 scission that outperform the current standards. This work improves our understanding of DNA repair gene function, and provides avenues for more precise modulation of Cas9-generated mutations.

Phosphorylation-mediated disassembly of C-terminal binding protein 2 tetramer impedes epigenetic silencing of pluripotency in mouse embryonic stem cells.

Cells need to overcome both intrinsic and extrinsic threats. Although pluripotency is associated with damage responses, how stem cells respond to DNA damage remains controversial. Here, we elucidate that DNA damage activates Chk2, leading to the phosphorylation of serine 164 on C-terminal binding protein 2 (Ctbp2). The phosphorylation of Ctbp2 induces the disruption of Ctbp2 tetramer, weakening interactions with zinc finger proteins, leading to the dissociation of phosphorylated Ctbp2 from chromatin. This transition to a monomeric state results in the separation of histone deacetylase 1 from Ctbp2, consequently slowing the rate of H3K27 deacetylation. In contrast to the nucleosome remodeling and deacetylase complex, phosphorylated Ctbp2 increased binding affinity to polycomb repressive complex (PRC)2, interacting through the N-terminal domain of Suz12. Through this domain, Ctbp2 competes with Jarid2, inhibiting the function of PRC2. Thus, the phosphorylation of Ctbp2 under stress conditions represents a precise mechanism aimed at preserving stemness traits by inhibiting permanent transcriptional shutdown.

CXXC5 stabilizes DNA methylation patterns in mouse embryonic stem cells.

Mammalian genomes encode 12 proteins that contain a CXXC zinc finger domain. Most members of this family are large multi-domain proteins that function in the control of DNA methylation and histone methylation patterns. CXXC5 is a smaller member of the family, along with its closest homologue CXXC4. These two proteins lack known catalytic domains. Here, we have characterized CXXC5 in mouse embryonic stem (ES) cells. We used gene knockouts, RNA sequencing, and DNA methylation analysis by whole-genome bisulfite sequencing. We show that CXXC5 is a nuclear protein that interacts with 5-methylcytosine oxidases (TET proteins). Removal of CXXC5 from ES cells leads to very few changes in gene expression. CXXC5 extensively colocalizes with TET1 and TET2 at CpG islands. CXXC5 inactivation leads to a substantial reduction of DNA methylation levels that affects all genomic compartments including genic and intergenic regions and CpG island shores. We propose a model in which CXXC5 serves as an anchor for TET proteins at CpG islands. In the absence of CXXC5, the 5-methylcytosine oxidases become dislodged from CpG islands and are enabled to induce genome-scale DNA demethylation. Thus, CXXC5 serves as a stabilizer of DNA methylation patterns.

Gastruloids are competent to specify both cardiac and skeletal muscle lineages

Cardiopharyngeal mesoderm contributes to the formation of the heart and head muscles. However, the mechanisms governing cardiopharyngeal mesoderm specification remain unclear. Here, we reproduce cardiopharyngeal mesoderm specification towards cardiac and skeletal muscle lineages with gastruloids from mouse embryonic stem cells. By conducting a comprehensive temporal analysis of cardiopharyngeal mesoderm development and differentiation in gastruloids compared to mouse embryos, we present the evidence for skeletal myogenesis in gastruloids. We identify different subpopulations of cardiomyocytes and skeletal muscles, the latter of which most likely correspond to different states of myogenesis with “head-like” and “trunk-like” skeletal myoblasts. In this work, we unveil the potential of gastruloids to undergo specification into both cardiac and skeletal muscle lineages, allowing the investigation of the mechanisms of cardiopharyngeal mesoderm differentiation in development and how this could be affected in congenital diseases.

A stepwise mode of TGFβ-SMAD signaling and DNA methylation regulates naïve-to-primed pluripotency and differentiation

The formation of transcription regulatory complexes by the association of Smad4 with Smad2 and Smad3 (Smad2/3) is crucial in the canonical TGFβ pathway. Although the central requirement of Smad4 as a common mediator is emphasized in regulating TGFβ signaling, it is not obligatory for all responses. The role of Smad2/3 independently of Smad4 remains understudied. Here, we introduce a stepwise paradigm in which Smad2/3 regulate the lineage priming and differentiation of mouse embryonic stem cells (mESCs) by collaboration with different effectors. During the naïve-to-primed transition, Smad2/3 upregulate DNA methyltransferase 3b (Dnmt3b), which establishes the proper DNA methylation patterns and, in turn, enables Smad2/3 binding to the hypomethylated centers of promoters and enhancers of epiblast marker genes. Consequently, in the absence of Smad2/3, Smad4 alone cannot initiate epiblast-specific gene transcription. When primed epiblast cells begin to differentiate, Dnmt3b becomes less actively engaged in global genome methylation, and Smad4 takes over the baton in this relay race, forming a complex with Smad2/3 to support mesendoderm induction. Thus, mESCs lacking Smad4 can undergo the priming process but struggle with the downstream differentiation. This work sheds light on the intricate mechanisms underlying TGFβ signaling and its role in cellular processes.

ZBTB24 is a conserved multifaceted transcription factor at genes and centromeres that governs the DNA methylation state and expression of satellite repeats.

Since its discovery as a causative gene of the Immunodeficiency with Centromeric instability and Facial anomalies syndrome, ZBTB24 has emerged as a key player in DNA methylation, immunity and development. By extensively analyzing ZBTB24 genomic functions in ICF-relevant mouse and human cellular models, we document here its multiple facets as a transcription factor, with key roles in immune response-related genes expression and also in early embryonic development. Using a constitutive Zbtb24 ICF-like mutant and an auxin-inducible degron system in mouse embryonic stem cells, we showed that ZBTB24 is recruited to centromeric satellite DNA where it is required to establish and maintain the correct DNA methylation patterns through the recruitment of DNMT3B. The ability of ZBTB24 to occupy centromeric satellite DNA is conserved in human cells. Together, our results unveiled an essential and underappreciated role for ZBTB24 at mouse and human centromeric satellite repeat arrays by controlling their DNA methylation and transcription status.

RADIP technology comprehensively identifies H3K27me3-associated RNA-chromatin interactions.

Many RNAs associate with chromatin, either directly or indirectly. Several technologies for mapping regions where RNAs interact across the genome have been developed to investigate the function of these RNAs. Obtaining information on the proteins involved in these RNA-chromatin interactions is critical for further analysis. Here, we developed RADIP [RNA and DNA interacting complexes ligated and sequenced (RADICL-seq) with immunoprecipitation], a novel technology that combines RADICL-seq technology with chromatin immunoprecipitation to characterize RNA-chromatin interactions mediated by individual proteins. Building upon the foundational principles of RADICL-seq, RADIP extends its advantages by increasing genomic coverage and unique mapping rate efficiency compared to existing methods. To demonstrate its effectiveness, we applied an anti-H3K27me3 antibody to the RADIP technology and generated libraries from mouse embryonic stem cells (mESCs). We identified a multitude of RNAs, including RNAs from protein-coding genes and non-coding RNAs, that are associated with chromatin via H3K27me3 and that likely facilitate the spread of Polycomb repressive complexes over broad regions of the mammalian genome, thereby affecting gene expression, chromatin structures and pluripotency of mESCs. Our study demonstrates the applicability of RADIP to investigations of the functions of chromatin-associated RNAs.

Lrif1 modulates Trim28-mediated repression of the Dux locus in mouse embryonic stem cells.

Germline mutations in SMCHD1, DNMT3B and LRIF1 can cause facioscapulohumeral muscular dystrophy type 2 (FSHD2). FSHD is an epigenetic skeletal muscle disorder in which partial failure in heterochromatinization of the D4Z4 macrosatellite repeat causes spurious expression of the repeat-embedded DUX4 gene in skeletal muscle, ultimately leading to muscle weakness and wasting. All three proteins play a role in chromatin organization and gene silencing; however, their functional relationship has not been fully elucidated. Here, we show that knockdown of Lrif1 , but not of the other two FSHD2 genes, in mouse embryonic stem cells leads to modest upregulation of the 2-cell cleavage stage transcriptional program driven by the transcription factor Dux, which is the mouse functional homologue of human DUX4. Furthermore, we show that Lrif1 interacts with Trim28, a known Dux repressor and that this interaction is independent of Cbx proteins and Smchd1. We uncover that modest Dux upregulation in Lrif1 knockdown mESCs coincides with decreased Trim28 occupancy at the Dux locus. Together, our results provide evidence for a conserved function of Lrif1 in repressing an early zygotic genome activator in mice and humans.

PRC1 and CTCF-Mediated Transition from Poised to Active Chromatin Loops Drives Bivalent Gene Activation.

Polycomb Repressive Complex 1 (PRC1) and CCCTC-binding factor (CTCF) are critical regulators of 3D chromatin architecture that influence cellular transcriptional programs. Spatial chromatin structures comprise conserved compartments, topologically associating domains (TADs), and dynamic, cell-type-specific chromatin loops. Although the role of CTCF in chromatin organization is well-known, the involvement of PRC1 is less understood. In this study, we identified an unexpected, essential role for the canonical Pcgf2-containing PRC1 complex (cPRC1.2), a known transcriptional repressor, in activating bivalent genes during differentiation. Our Hi-C analysis revealed that cPRC1.2 forms chromatin loops at bivalent promoters, rendering them silent yet poised for activation. Using mouse embryonic stem cells (ESCs) with CRISPR/Cas9-mediated gene editing, we found that the loss of Pcgf2, though not affecting the global level of H2AK119ub1, disrupts these cPRC1.2 loops in ESCs and impairs the transcriptional induction of crucial target genes necessary for neuronal differentiation. Furthermore, we identified CTCF enrichment at cPRC1.2 loop anchors and at Polycomb group (PcG) bodies, nuclear foci with concentrated PRC1 and its tethered chromatin domains, suggesting that PRC1 and CTCF cooperatively shape chromatin loop structures. Through virtual 4C and other genomic analyses, we discovered that establishing neuronal progenitor cell (NPC) identity involves a switch from cPRC1.2-mediated chromatin loops to CTCF-mediated active loops, enabling the expression of critical lineage-specific factors. This study uncovers a novel mechanism by which pre-formed PRC1 and CTCF loops at lineage-specific genes maintain a poised state for subsequent gene activation, advancing our understanding of the role of chromatin architecture in controlling cell fate transitions.

Subcellular mRNA kinetic modeling reveals nuclear retention as rate-limiting

Eukaryotic mRNAs are transcribed, processed, translated, and degraded in different subcellular compartments. Here, we measured mRNA flow rates between subcellular compartments in mouse embryonic stem cells. By combining metabolic RNA labeling, biochemical fractionation, mRNA sequencing, and mathematical modeling, we determined the half-lives of nuclear pre-, nuclear mature, cytosolic, and membrane-associated mRNAs from over 9000 genes. In addition, we estimated transcript elongation rates. Many matured mRNAs have long nuclear half-lives, indicating nuclear retention as the rate-limiting step in the flow of mRNAs. In contrast, mRNA transcripts coding for transcription factors show fast kinetic rates, and in particular short nuclear half-lives. Differentially localized mRNAs have distinct rate constant combinations, implying modular regulation. Membrane stability is high for membrane-localized mRNA and cytosolic stability is high for cytosol-localized mRNA. mRNAs encoding target signals for membranes have low cytosolic and high membrane half-lives with minor differences between signals. Transcripts of nuclear-encoded mitochondrial proteins have long nuclear retention and cytoplasmic kinetics that do not reflect co-translational targeting. Our data and analyses provide a useful resource to study spatiotemporal gene expression regulation.

Lactylation of Hdac1 regulated by Ldh prevents the pluripotent-to-2C state conversion

BackgroundCellular metabolism regulates the pluripotency of embryonic stem cells (ESCs). Yet, how metabolism regulates the transition among different pluripotent states remains elusive. It has been shown that protein lactylation, which uses lactate, a metabolic product of glycolysis, as a substrate, plays a critical role in various biological events. Here we focused on that glycolysis regulates the conversion between ESCs and 2-cell-like cells (2CLCs) through protein lactylation.MethodsRNA-seq revealed the activation of 2-cell (2C) genes by suppression of Ldh. Stable isotope labeling by amino acids in cell culture (SILAC) coupled with lactylated peptide enrichment and quantitative mass spectrometric analysis was carried out to investigate the mechanism how protein lactylation regulates the pluripotent-to-2C transition. And we focused on Hdac1. Lactylation of Hdac1 required for silencing 2C genes was proved by quantitative reverse-transcription PCR (qRT-PCR), immunofluorescence (IF), Western blot and chimeric embryos. Chromatin immunoprecipitation coupled with sequencing (ChIP-seq) and in vitro deacetylation assay confirmed lactylation of Hdac1 promoting its binding at 2C genes and enhancing its deacetylase activity, thereby facilitating the removal of H3K27ac and the silencing of 2C genes.ResultsWe found that inhibition or depletion of Ldha, the enzyme converting pyruvate to lactate, leads to the activation of 2C genes, as well as reduced global lactylation in ESCs. To investigate the mechanism how protein lactylation regulates the pluripotent-to-2C transition, quantitative lactylome analysis was performed, and 1716 lactylated proteins were identified. We then focused on Hdac1, a histone deacetylase involved in the silencing of 2C genes. Lactylation of Hdac1 promotes its binding at 2C genes and enhances its deacetylase activity, thus facilitating the removal of H3K27ac and the silencing of 2C genes.ConclusionsIn summary, our study reveals a mechanistic link between cellular metabolism and pluripotency regulation through protein lactylation. Our research is the first time to reveal that quantitative lactylome analysis in mouse ESCs. We found that lactylated Hdac1 promotes its binding at 2C genes and enhances its deacetylase activity, thus facilitating the removal of H3K27ac and the silencing of 2C genes.

Repeat expansion in a Fragile X model is independent of double strand break repair mediated by Pol θ, Rad52, Rad54l or Rad54b.

Microsatellite instability is responsible for the human Repeat Expansion Disorders. The mutation responsible differs from classical cancer-associated microsatellite instability (MSI) in that it requires the mismatch repair proteins that normally protect against MSI. LIG4, an enzyme essential for non-homologous end-joining (NHEJ), the major pathway for double-strand break repair (DSBR) in mammalian cells, protects against expansion in mouse models. Thus, NHEJ may compete with the expansion pathway for access to a common intermediate. This raises the possibility that expansion involves an NHEJ-independent form of DSBR. Pol θ, a polymerase involved in the theta-mediated end joining (TMEJ) DSBR pathway, has been proposed to play a role in repeat expansion. Here we examine the effect of the loss of Pol θ on expansion in FXD mouse embryonic stem cells (mESCs), along with the effects of mutations in Rad52 , Rad54l and Rad54b, genes important for multiple DSBR pathways. None of these mutations significantly affected repeat expansion. These observations put major constraints on what pathways are likely to drive expansion. Together with our previous demonstration of the protective effect of nucleases like EXO1 and FAN1, and the importance of Pol β, they suggest a plausible model for late steps in the expansion process.

Copper Oxide Nanoparticles Impair Mouse Preimplantation Embryonic Development through Disruption of Mitophagy-Mediated Metabolism.

Copper oxide nanoparticles (CuONPs) have been widely applied, posing potential risks to human health. Although the toxicity of CuONPs on the liver and spleen has been reported, their effects on reproductive health remain unexplored. In this study, we investigate the effects of CuONPs on embryonic development and their potential mechanisms. Our results demonstrate that CuONPs exposure impairs mouse preimplantation embryonic development, particularly affecting the morula-to-blastocyst transition. Additionally, CuONPs were found to reduce the pluripotency of the inner cell mass (ICM) and mouse embryonic stem cells (mESCs). Mechanistically, CuONPs block autophagic flux and impair mitophagy, leading to the accumulation of damaged mitochondria. This mitochondrial dysfunction leads to reduced tricarboxylic acid (TCA) cycle activity and decreased α-ketoglutarate (α-KG) production. Insufficient α-KG induces the failure of DNA demethylation, reducing corresponding chromatin accessibility and consequently inhibiting ICM-specific genes expressions. Similar reduced development and inhibitions of pluripotency gene expression were observed in CuONPs-treated human blastocysts. Moreover, in women undergoing assisted reproductive technology (ART), a negative correlation was found between urinary Cu ion concentrations and clinical outcomes. Collectively, our study elucidates the mitophagy-mediated metabolic mechanisms of CuONPs embryotoxicity, improving our understanding of the potential reproductive toxicity associated with it.

H3.3K122A results in a neomorphic phenotype in mouse embryonic stem cells

Canonical histone H3 and histone variant H3.3 are posttranslationally modified with the genomic distribution of these marks denoting different features and these modifications may influence transcription. While the majority of posttranslational modifications occur on histone tails, there are defined modifications within the globular domain, such as acetylation of H3K122/H3.3K122. To understand the function of the amino acid H3.3K122 in transcriptional regulation, we attempted to generate H3.3K122A mouse embryonic stem (mES) cells but were unsuccessful. Through multi-omic profiling of mutant cell lines harboring two or three of four H3.3 targeted alleles, we have uncovered that H3.3K122A is neomorphic and results in lethality. This is surprising as prior studies demonstrate H3.3-null mES cells are viable and pluripotent but exhibit a reduced differentiation capacity. Together, these studies have uncovered a novel dependence of a globular domain residue within H3.3 for viability and broadened our understanding of how histone variants contribute to transcription regulation and pluripotency in mES cells.

Differential modulation of polycomb-associated histone marks by cBAF, pBAF, and gBAF complexes.

Chromatin regulators alter the physical properties of chromatin to make it more or less permissive to transcription by modulating another protein's access to a specific DNA sequence through changes in nucleosome occupancy or histone modifications at a particular locus. Mammalian SWI/SNF complexes are a group of ATPase-dependent chromatin remodelers. In mouse embryonic stem cells, there are three primary forms of mSWI/SNF: canonical BAF (cBAF), polybromo-associated BAF (pBAF), and GLTSCR-associated BAF (gBAF). Nkx2-9 is bivalent, meaning nucleosomes at the locus have active and repressive modifications. In this study, we used unique BAF subunits to recruit each of the three complexes to Nkx2-9 using dCas9-mediated inducible recruitment (FIRE-Cas9). We show that recruitment of cBAF complexes leads to a significant loss of the polycomb repressive-2 H3K27me3 histone mark and polycomb repressive-1 and repressive-2 complex proteins, whereas gBAF and pBAF do not. Moreover, nucleosome occupancy alone cannot explain the loss of these marks. Our results demonstrate that cBAF has a unique role in the direct opposition of polycomb-associated histone modifications that gBAF and pBAF do not share.

STAG2 promotes naive-primed transition via activating Lin28a transcription in mouse embryonic stem cells

Mouse embryonic stem cells (mESCs) exist in two distinct pluripotent states: the naive and the primed. Mainly by inducing differentiation of mESCs in vitro, conducting RNA sequencing analyses, and specifying expression of the regulatory genes, we explored the regulatory mechanisms underlying the transition between the naive and primed states. We found that, under the defined differentiation-inducing conditions, the naive state of mESCs shifted to the primed state within two days of differentiation induction, during which the cell cycle- and differentiation-related proteins changes significantly. Specifically, we uncovered that the expression of STAG2, a subunit of the Cohesin complex, was upregulated. We further revealed that knockout of STAG2 resulted in upregulation of the naive gene sets and downregulation of the primed gene sets, indicating importance of STAG2 in regulating the naive-primed transition. More importantly, STAG2 knockout led to a reduction in number of the bivalent genes, a decrease in Lin28a transcription, and a reduced cytoplasmic localization of Lin28a. Overexpressing Lin28a or a Lin28a variant lacking the nucleolar localization signal (Lin28aΔNoLS) in STAG2 knockout cells rescued the downregulation of the primed marker genes Dnmt3a/3b. Collectively, we conclude that STAG2 facilitates the naive-primed transition of mESCs by activating Lin28a transcription and that this work may offer a new insight into the regulation of pluripotency in mESCs.