Inhibition of protein kinase C increases Prdm14 level to promote self-renewal of embryonic stem cells through reducing Suv39h-induced H3K9 methylation

Inhibition of protein kinase C (PKC) efficiently promoted the self-renewal of embryonic stem cells (ESCs). However, information about the function of PKC inhibition remains lacking. Here, RNA-sequencing showed that the addition of Go6983 significantly inhibited the expression of de novo methyltransferases (Dnmt3a and Dnmt3b) and their regulator Dnmt3l, resulting in global hypomethylation of DNA in mouse ESCs. Mechanistically, PR domain-containing 14 (Prdm14), a site-specific transcriptional activator, partially contributed to Go6983-mediated repression of Dnmt3 genes. Administration of Go6983 increased Prdm14 expression mainly through the inhibition of PKCδ. High constitutive expression of Prdm14 phenocopied the ability of Go6983 to maintain` mouse ESC stemness in the absence of self-renewal-promoting cytokines. In contrast, the knockdown of Prdm14 eliminated the response to PKC inhibition and substantially impaired the Go6983-induced resistance of mouse ESCs to differentiation. Furthermore, liquid chromatography-mass spectrometry profiling and Western blotting revealed low levels of Suv39h1 and Suv39h2 in Go6983-treated mouse ESCs. Suv39h enzymes are histone methyltransferases that recognize dimethylated and trimethylated histone H3K9 specifically and usually function as transcriptional repressors. Consistently, the inhibition of Suv39h1 by RNA interference or the addition of the selective inhibitor chaetocin increased Prdm14 expression. Moreover, chromatin immunoprecipitation assay showed that Go6983 treatment led to decreased enrichment of dimethylation and trimethylation of H3K9 at the Prdm14 promoter but increased RNA polymerase Ⅱ binding affinity. Together, our results provide novel insights into the pivotal association between PKC inhibition-mediated self-renewal and epigenetic changes, which will help us better understand the regulatory network of stem cell pluripotency.

CNOT3 interacts with the Aurora B and MAPK/ERK kinases to promote survival of differentiating mesendodermal progenitor cells

Mesendoderm cells are key intermediate progenitors that form at the early primitive streak (PrS) and give rise to mesoderm and endoderm in the gastrulating embryo. We have identified an interaction between CNOT3 and the cell cycle kinase Aurora B that requires sequences in the NOT box domain of CNOT3 and regulates MAPK/ERK signaling during mesendoderm differentiation. Aurora B phosphorylates CNOT3 at two sites located close to a nuclear localization signal and promotes localization of CNOT3 to the nuclei of mouse embryonic stem cells (ESCs) and metastatic lung cancer cells. ESCs that have both sites mutated give rise to embryoid bodies that are largely devoid of mesoderm and endoderm and are composed mainly of cells with ectodermal characteristics. The mutant ESCs are also compromised in their ability to differentiate into mesendoderm in response to FGF2, BMP4, and Wnt3 due to reduced survival and proliferation of differentiating mesendoderm cells. We also show that the double mutation alters the balance of interaction of CNOT3 with Aurora B and with ERK and reduces phosphorylation of ERK in response to FGF2. Our results identify a potential adaptor function for CNOT3 that regulates the Ras/MEK/ERK pathway during embryogenesis.

Inhibition of LINE-1 retrotransposition represses telomere reprogramming during mouse 2-cell embryo development

PURPOSE

To investigate whether inhibition of LINE-1 affects telomere reprogramming during 2-cell embryo development.

METHODS

Mouse zygotes were cultured with or without 1 µM azidothymidine (AZT) for up to 15 h (early 2-cell, G1/S) or 24 h (late 2-cell, S/G2). Gene expression and DNA copy number were determined by RT-qPCR and qPCR respectively. Immunostaining and telomeric PNA-FISH were performed for co-localization between telomeres and ZSCAN4 or LINE-1-Orf1p.

RESULTS

LINE-1 copy number was remarkably reduced in later 2-cell embryos by exposure to 1 µM AZT, and telomere lengths in late 2-cell embryos with AZT were significantly shorter compared to control embryos (P = 0.0002). Additionally, in the absence of LINE-1 inhibition, Dux, Zscan4, and LINE-1 were highly transcribed in early 2-cell embryos, as compared to late 2-cell embryos (P < 0.0001), suggesting that these 2-cell genes are activated at the early 2-cell stage. However, in early 2-cell embryos with AZT treatment, mRNA levels of Dux, Zscan4, and LINE-1 were significantly decreased. Furthermore, both Zscan4 and LINE-1 encoded proteins localized to telomere regions in 2-cell embryos, but this co-localization was dramatically reduced after AZT treatment (P < 0.001).

CONCLUSIONS

Upon inhibition of LINE-1 retrotransposition in mouse 2-cell embryos, Dux, Zscan4, and LINE-1 were significantly downregulated, and telomere elongation was blocked. ZSCAN4 foci and their co-localization with telomeres were also significantly decreased, indicating that ZSCAN4 is an essential component of the telomere reprogramming that occurs in mice at the 2-cell stage. Our findings also suggest that LINE-1 may directly contribute to telomere reprogramming in addition to regulating gene expression.

PRC1 drives Polycomb-mediated gene repression by controlling transcription initiation and burst frequency

The Polycomb repressive system plays a fundamental role in controlling gene expression during mammalian development. To achieve this, Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) bind target genes and use histone modification-dependent feedback mechanisms to form Polycomb chromatin domains and repress transcription. The inter-relatedness of PRC1 and PRC2 activity at these sites has made it difficult to discover the specific components of Polycomb chromatin domains that drive gene repression and to understand mechanistically how this is achieved. Here, by exploiting rapid degron-based approaches and time-resolved genomics, we kinetically dissect Polycomb-mediated repression and discover that PRC1 functions independently of PRC2 to counteract RNA polymerase II binding and transcription initiation. Using single-cell gene expression analysis, we reveal that PRC1 acts uniformly within the cell population and that repression is achieved by controlling transcriptional burst frequency. These important new discoveries provide a mechanistic and conceptual framework for Polycomb-dependent transcriptional control.

p53 convergently activates Dux/DUX4 in embryonic stem cells and in facioscapulohumeral muscular dystrophy cell models

In mammalian embryos, proper zygotic genome activation (ZGA) underlies totipotent development. Double homeobox (DUX)-family factors participate in ZGA, and mouse Dux is required for forming cultured two-cell (2C)-like cells. Remarkably, in mouse embryonic stem cells, Dux is activated by the tumor suppressor p53, and Dux expression promotes differentiation into expanded-fate cell types. Long-read sequencing and assembly of the mouse Dux locus reveals its complex chromatin regulation including putative positive and negative feedback loops. We show that the p53-DUX/DUX4 regulatory axis is conserved in humans. Furthermore, we demonstrate that cells derived from patients with facioscapulohumeral muscular dystrophy (FSHD) activate human DUX4 during p53 signaling via a p53-binding site in a primate-specific subtelomeric long terminal repeat (LTR)10C element. In summary, our work shows that p53 activation convergently evolved to couple p53 to Dux/DUX4 activation in embryonic stem cells, embryos and cells from patients with FSHD, potentially uniting the developmental and disease regulation of DUX-family factors and identifying evidence-based therapeutic opportunities for FSHD.

Transcriptional profiling of mESC-derived tendon and fibrocartilage cell fate switch

The transcriptional regulators underlying induction and differentiation of dense connective tissues such as tendon and related fibrocartilaginous tissues (meniscus and annulus fibrosus) remain largely unknown. Using an iterative approach informed by developmental cues and single cell RNA sequencing (scRNA-seq), we establish directed differentiation models to generate tendon and fibrocartilage cells from mouse embryonic stem cells (mESCs) by activation of TGFβ and hedgehog pathways, achieving 90% induction efficiency. Transcriptional signatures of the mESC-derived cells recapitulate embryonic tendon and fibrocartilage signatures from the mouse tail. scRNA-seq further identify retinoic acid signaling as a critical regulator of cell fate switch between TGFβ-induced tendon and fibrocartilage lineages. Trajectory analysis by RNA sequencing define transcriptional modules underlying tendon and fibrocartilage fate induction and identify molecules associated with lineage-specific differentiation. Finally, we successfully generate 3-dimensional engineered tissues using these differentiation protocols and show activation of mechanotransduction markers with dynamic tensile loading. These findings provide a serum-free approach to generate tendon and fibrocartilage cells and tissues at high efficiency for modeling development and disease.

Epithelial cell plasticity drives endoderm formation during gastrulation

It is generally accepted that epiblast cells ingress into the primitive streak by epithelial-to-mesenchymal transition (EMT) to give rise to the mesoderm; however, it is less clear how the endoderm acquires an epithelial fate. Here, we used embryonic stem cell and mouse embryo knock-in reporter systems to combine time-resolved lineage labelling with high-resolution single-cell transcriptomics. This allowed us to resolve the morphogenetic programs that segregate the mesoderm from the endoderm germ layer. Strikingly, while the mesoderm is formed by classical EMT, the endoderm is formed independent of the key EMT transcription factor Snail1 by mechanisms of epithelial cell plasticity. Importantly, forkhead box transcription factor A2 (Foxa2) acts as an epithelial gatekeeper and EMT suppressor to shield the endoderm from undergoing a mesenchymal transition. Altogether, these results not only establish the morphogenetic details of germ layer formation, but also have broader implications for stem cell differentiation and cancer metastasis.

DNMTs Play an Important Role in Maintaining the Pluripotency of Leukemia Inhibitory Factor-Dependent Embryonic Stem Cells

Naive pluripotency can be maintained in medium with two inhibitors plus leukemia inhibitory factor (2i/LIF) supplementation, which primarily affects canonical WNT, FGF/ERK, and JAK/STAT3 signaling. However, whether one of these three supplements alone is sufficient to maintain naive self-renewal remains unclear. Here we show that LIF alone in medium is sufficient for adaptation of 2i/L-ESCs to embryonic stem cells (ESCs) in a hypermethylated state (L-ESCs). Global transcriptomic analysis shows that L-ESCs are close to 2i/L-ESCs and in a stable state between naive and primed pluripotency. Notably, our results demonstrate that DNA methyltransferases (DNMTs) play an important role in LIF-dependent mouse ESC adaptation and self-renewal. LIF-dependent ESC adaptation efficiency is significantly increased in serum treatment and reduced in Dnmt3a or Dnmt3l knockout ESCs. Importantly, unlike epiblast stem cells, L-ESCs contribute to somatic tissues and germ cells in chimeras. L-ESCs cultured under such simple conditions as in this study would provide a more conducive platform to clarify the molecular mechanism of ESCs in in vitro culture.

•

LIF alone supports ESC self-renewal and pluripotency in chemically defined media

•

L-ESCs re-establish the epigenetic state in LIF adaptation

•

DNMTs are important for LIF adaptation and L-ESC self-renewal

•

L-ESCs contribute to somatic tissues and germ cells in chimeras

Mode of action assessment of the genotoxic properties of antimony and its compounds evaluated in the ToxTracker assay

Antimony (Sb) and its compounds are negative in gene mutation assays in bacteria and cultured mammalian cells but positive in some assays for clastogenicity and/or DNA damage. In order to better understand the modes of action for antimony genotoxicity, we assessed reporter gene activation by antimony and antimony compounds in the new expanded ToxTracker assay. ToxTracker evaluates the activation of biomarkers for different cellular defense mechanisms using a series of green fluorescent protein reporters inserted into mouse embryonic stem cell lines. The assay responds to: 1) DNA damage and inhibition of DNA replication; 2) oxidative stress; 3) unfolded protein response (protein damage); and 4) p53-dependent cellular stress. Sb metal powder, six trivalent (Sb(III)) compounds, and five pentavalent antimony (Sb(V)) compounds were assessed. Sb powder and all six Sb(III) compounds activated oxidative stress ToxTracker reporters at non-toxic doses. Of the five Sb(V) compounds, antimony pentachloride and potassium hexahydroantimonate induced a robust oxidative stress response while sodium antimonate induced only a weak oxidative stress response. At higher concentrations (up to either 75 % toxicity or the highest dissolved concentration tested), Sb powder and all Sb(III) compounds except for antimony trichloride induced the unfolded protein response. Of the five Sb(V) compounds tested, only potassium hexahydroantimonate induced weak activation of the unfolded protein response and was also the only pentavalent compound to yield modest (30 %) cytotoxicity. None of the compounds tested activated the DNA damage/inhibition of DNA replication reporters, nor did they activate the p53-dependent response. All Sb(III) compounds, Sb powder, and three of the five Sb(V) compounds activated the oxidative stress reporters, but there was no activation of reporters associated with DNA damage and repair or p53-dependent cellular stress. The consistent activation of reporters for oxidative stress suggests this mode of action may underlie genotoxicity responses for antimony and its compounds.

Understanding the role of Beclin1 in mouse embryonic stem cell differentiation through CRISPR-Cas9-mediated gene editing

Autophagy is a vacuolar pathway for the regulated degradation and recycling of cellular components. Beclin1, a Bcl2-interacting protein, is a well-studied autophagy regulator. Homozygous loss of Beclin1 in mice leads to early embryonic lethality. However, the role of Beclin1 in regulating the pluripotency of embryonic stem cells and their differentiation remains poorly explored. To study this, we generated Beclin1-Knockout (KO) mouse embryonic stem cells (mESCs) using the CRISPR-Cas9 genome-editing tool. Interestingly, Beclin1-KO mESCs did not show any change in the expression of pluripotency marker genes. Beclin1-KO mESCs also displayed active autophagy, suggesting the presence of Beclin1-independent autophagy in mESCs. However, loss of Beclin1 resulted in compromised differentiation of mESCs in vitro and in vivo due to misregulated expression of transcription factors. Our results suggest that Beclin1 may play an autophagy-independent role in regulating the differentiation of mESCs.

Primordial germ cells can be differentiated by retinoic acid and progesterone induction from embryonic stem cells

This study aimed to examine the expression of the genes associated with different development stages of primordial germ cells (PGCs) in differentiating mouse embryonic stem cells (mESCs). The cells were cultured in three groups of control, 10^-8 M of all-trans retinoic acid and the combination of 10^-7 M of Progesterone and retinoic acid for 7, 12, 17, and 22 days. Immunofluorescent and Quantitative RT-PCR were used to evaluate the effect of progesterone on the differentiation of mESCs into primordial germ cells. RA-treated cells exhibited increased expression of Fragilis, Stella, Dazl, Stra8, Sycp3, and Gdf9 genes and decreased expression of Oct4, Mvh genes compared to the non-treated controls. Furthermore, RA in combination with progesterone (RA?P) led to increased expression of Oct4, Fragilis, Stella, Dazl, Sycp3, Gdf9 and decreased expression of Mvh, and Stra8 genes compared to the RA-treated scenario. Immunofluorescence detection of Stella and Mvh showed that the expression levels of the cells treated with RA+P are much higher than those of the other groups. Our project showed that under the influence of the induced factors, mESCs can spontaneously differentiate into germ cells. Also, the combination of RA+P can enhance and accelerate the differentiation of mESCs into germ cells.

Mouse embryonic stem cells self-organize into trunk-like structures with neural tube and somites

Post-implantation embryogenesis is a highly dynamic process comprising multiple lineage decisions and morphogenetic changes that are inaccessible to deep analysis in vivo. We found that pluripotent mouse embryonic stem cells (mESCs) form aggregates that upon embedding in an extracellular matrix compound induce the formation of highly organized "trunk-like structures" (TLSs) comprising the neural tube and somites. Comparative single-cell RNA sequencing analysis confirmed that this process is highly analogous to mouse development and follows the same stepwise gene-regulatory program. Tbx6 knockout TLSs developed additional neural tubes mirroring the embryonic mutant phenotype, and chemical modulation could induce excess somite formation. TLSs thus reveal an advanced level of self-organization and provide a powerful platform for investigating post-implantation embryogenesis in a dish.

Direct analysis of brain phenotypes via neural blastocyst complementation

We provide a protocol for generating forebrain structures in vivo from mouse embryonic stem cells (ESCs) via neural blastocyst complementation (NBC). We developed this protocol for studies of development and function of specific forebrain regions, including the cerebral cortex and hippocampus. We describe a complete workflow, from methods for modifying a given genomic locus in ESCs via CRISPR-Cas9-mediated editing to the generation of mouse chimeras with ESC-reconstituted forebrain regions that can be directly analyzed. The procedure begins with genetic editing of mouse ESCs via CRISPR-Cas9, which can be accomplished in ~4-8 weeks. We provide protocols to achieve fluorescent labeling of ESCs in ~2-3 weeks, which allows tracing of the injected, ESC-derived donor cells in chimeras generated via NBC. Once modified ESCs are ready, NBC chimeras are generated in ~3 weeks via injection of ESCs into genetically programmed blastocysts that are subsequently transferred into pseudo-pregnant fosters. Our in vivo brain organogenesis platform is efficient, allowing functional and systematic analysis of genes and other genomic factors in as little as 3 months, in the context of a whole organism.

Intracellular Ca2+ Homeostasis and Nuclear Export Mediate Exit from Naive Pluripotency

Progression through states of pluripotency is required for cells in early mammalian embryos to transition away from heightened self-renewal and toward competency for lineage specification. Here, we use a CRISPR mutagenesis screen in mouse embryonic stem cells (ESCs) to identify unexpected roles for nuclear export and intracellular Ca2+ homeostasis during the exit out of the naive state of pluripotency. Mutation of a plasma membrane Ca2+ pump encoded by Atp2b1 increased intracellular Ca2+ such that it overcame effects of intracellular Ca2+ reduction, which is required for naive exit. Persistent self-renewal of ESCs was supported both in Atp2b1-/-Tcf7l1-/- double-knockout ESCs passaged in defined media alone (no LIF or inhibitors) and in wild-type cells passaged in media containing only calcitonin and a GSK3 inhibitor. These new findings suggest a central role for intracellular Ca2+ in safeguarding naive pluripotency.

Derivation and maintenance of mouse haploid embryonic stem cells

Ploidy represents the number of chromosome sets in a cell. Although gametes have a haploid genome (n), most mammalian cells have diploid genomes (2n). The diploid status of most cells correlates with the number of probable alleles for each autosomal gene and makes it difficult to target these genes via mutagenesis techniques. Here, we describe a 7-week protocol for the derivation of mouse haploid embryonic stem cells (hESCs) from female gametes that also outlines how to maintain the cells once derived. We detail additional procedures that can be used with cell lines obtained from the mouse Haplobank, a biobank of >100,000 individual mouse hESC lines with targeted mutations in 16,970 genes. hESCs can spontaneously diploidize and can be maintained in both haploid and diploid states. Mouse hESCs are genomically and karyotypically stable, are innately immortal and isogenic, and can be derived in an array of differentiated cell types; they are thus highly amenable to genetic screens and to defining molecular connectivity pathways.

Compromised Chondrocyte Differentiation Capacity in TERC Knockout Mouse Embryonic Stem Cells Derived by Somatic Cell Nuclear Transfer

Mammalian telomere lengths are primarily regulated by telomerase, consisting of a reverse transcriptase protein (TERT) and an RNA subunit (TERC). We previously reported the generation of mouse Terc^+/- and Terc^-/- embryonic stem cells (ntESCs) by somatic cell nuclear transfer. In the present work, we investigated the germ layer development competence of Terc^-/-, Terc^+/- and wild-type (Terc^+/+) ntESCs. The telomere lengths are longest in wild-type but shortest in Terc^-/- ntESCs, and correlate reversely with the population doubling time. Interestingly, while in vitro embryoid body (EB) differentiation assay reveals EB size difference among ntESCs of different genotypes, the more stringent in vivo teratoma assay demonstrates that Terc^-/- ntESCs are severely defective in differentiating into the mesodermal lineage cartilage. Consistently, in a directed in vitro chondrocyte differentiation assay, the Terc^-/- cells failed in forming Collagen II expressing cells. These findings underscore the significance in maintaining proper telomere lengths in stem cells and their derivatives for regenerative medicine.

Promoter-Enhancer Communication Occurs Primarily within Insulated Neighborhoods

Metazoan chromosomes are sequentially partitioned into topologically associating domains (TADs) and then into smaller sub-domains. One class of sub-domains, insulated neighborhoods, are proposed to spatially sequester and insulate the enclosed genes through self-association and chromatin looping. However, it has not been determined functionally whether promoter-enhancer interactions and gene regulation are broadly restricted to within these loops. Here, we employed published datasets from murine embryonic stem cells (mESCs) to identify insulated neighborhoods that confine promoter-enhancer interactions and demarcate gene regulatory regions. To directly address the functionality of these regions, we depleted estrogen-related receptor β (Esrrb), which binds the Mediator co-activator complex, to impair enhancers of genes within 222 insulated neighborhoods without causing mESC differentiation. Esrrb depletion reduces Mediator binding, promoter-enhancer looping, and expression of both nascent RNA and mRNA within the insulated neighborhoods without significantly affecting the flanking genes. Our data indicate that insulated neighborhoods represent functional regulons in mammalian genomes.

Analysis of Technical and Biological Variability in Single-Cell RNA Sequencing

Profiling the transcriptomes of individual cells with single-cell RNA sequencing (scRNA-seq) has been widely applied to provide a detailed molecular characterization of cellular heterogeneity within a population of cells. Despite recent technological advances of scRNA-seq, technical variability of gene expression in scRNA-seq is still much higher than that in bulk RNA-seq. Accounting for technical variability is therefore a prerequisite for correctly analyzing single-cell data. This chapter describes a computational pipeline for detecting highly variable genes exhibiting higher cell-to-cell variability than expected by technical noise. The basic pipeline using the scater and scran R/Bioconductor packages includes deconvolution-based normalization, fitting the mean-variance trend, testing for nonzero biological variability, and visualization with highly variable genes. An outline of the underlying theory of detecting highly variable genes is also presented. We illustrate how the pipeline works by using two case studies, one from mouse embryonic stem cells with external RNA spike-ins, and the other from mouse dentate gyrus cells without spike-ins.

Heterochromatin delays CRISPR-Cas9 mutagenesis but does not influence the outcome of mutagenic DNA repair

Genome editing occurs in the context of chromatin, which is heterogeneous in structure and function across the genome. Chromatin heterogeneity is thought to affect genome editing efficiency, but this has been challenging to quantify due to the presence of confounding variables. Here, we develop a method that exploits the allele-specific chromatin status of imprinted genes in order to address this problem in cycling mouse embryonic stem cells (mESCs). Because maternal and paternal alleles of imprinted genes have identical DNA sequence and are situated in the same nucleus, allele-specific differences in the frequency and spectrum of mutations induced by CRISPR-Cas9 can be unequivocally attributed to epigenetic mechanisms. We found that heterochromatin can impede mutagenesis, but to a degree that depends on other key experimental parameters. Mutagenesis was impeded by up to 7-fold when Cas9 exposure was brief and when intracellular Cas9 expression was low. In contrast, the outcome of mutagenic DNA repair was unaffected by chromatin state, with similar efficiencies of homology-directed repair (HDR) and deletion spectra on maternal and paternal chromosomes. Combined, our data show that heterochromatin imposes a permeable barrier that influences the kinetics, but not the endpoint, of CRISPR-Cas9 genome editing and suggest that therapeutic applications involving low-level Cas9 exposure will be particularly affected by chromatin status.

Gain of CTCF-Anchored Chromatin Loops Marks the Exit from Naive Pluripotency

The genome of pluripotent stem cells adopts a unique three-dimensional architecture featuring weakly condensed heterochromatin and large nucleosome-free regions. Yet, it is unknown whether structural loops and contact domains display characteristics that distinguish embryonic stem cells (ESCs) from differentiated cell types. We used genome-wide chromosome conformation capture and super-resolution imaging to determine nuclear organization in mouse ESC and neural stem cell (NSC) derivatives. We found that loss of pluripotency is accompanied by widespread gain of structural loops. This general architectural change correlates with enhanced binding of CTCF and cohesins and more pronounced insulation of contacts across chromatin boundaries in lineage-committed cells. Reprogramming NSCs to pluripotency restores the unique features of ESC domain topology. Domains defined by the anchors of loops established upon differentiation are enriched for developmental genes. Chromatin loop formation is a pervasive structural alteration to the genome that accompanies exit from pluripotency and delineates the spatial segregation of developmentally regulated genes.

Modulation of STAT3 phosphorylation by PTPN2 inhibits naïve pluripotency of embryonic stem cells

STAT3 phosphorylation at tyrosine 705 (STAT3^pY705 ), triggered by the addition of the leukemia inhibitory factor (LIF), can maintain mouse embryonic stem cell (mESC) self-renewal and reprogram mouse epiblast stem cells (EpiSCs) to enter a naïve pluripotent state. The activation of STAT3^pY705 occurs mainly through Janus kinases. However, it remains unclear how STAT3^pY705 levels are decreased in mESCs. Our study shows that upregulation of the protein tyrosine phosphatase (PTPN2) inhibits STAT3 activity by reducing its phosphorylation level and promotes mESC differentiation, whereas PTPN2 knockout by CRISPR/CAS9 delays mESC differentiation. Consistently, PTPN2 knockdown facilitates the generation of mESC-like colonies in STAT3-overexpressing EpiSCs. PTPN2-mediated STAT3 activity, thus, contributes to the exit of ESCs from the pluripotent ground state. These findings expand the current understanding of the regulatory network of naïve pluripotency.

RAS Regulates the Transition from Naive to Primed Pluripotent Stem Cells

The transition from naive to primed state of pluripotent stem cells is hallmarked by epithelial-mesenchymal transition, metabolic switch from oxidative phosphorylation to aerobic glycolysis, and changes in the epigenetic landscape. Since these changes are also seen as putative hallmarks of neoplastic cell transformation, we hypothesized that oncogenic pathways may be involved in this process. We report that the activity of RAS is repressed in the naive state of mouse embryonic stem cells (ESCs) and that all three RAS isoforms are significantly activated upon early differentiation induced by LIF withdrawal, embryoid body formation, or transition to the primed state. Forced expression of active RAS and RAS inhibition have shown that RAS regulates glycolysis, CADHERIN expression, and the expression of repressive epigenetic marks in pluripotent stem cells. Altogether, this study indicates that RAS is located at a key junction of early ESC differentiation controlling key processes in priming of naive cells.

Neurons Generated by Mouse ESCs with Hippocampal or Cortical Identity Display Distinct Projection Patterns When Co-transplanted in the Adult Brain

The capability of generating neural precursor cells with distinct types of regional identity in vitro has recently opened new opportunities for cell replacement in animal models of neurodegenerative diseases. By manipulating Wnt and BMP signaling, we steered the differentiation of mouse embryonic stem cells (ESCs) toward isocortical or hippocampal molecular identity. These two types of cells showed different degrees of axonal outgrowth and targeted different regions when co-transplanted in healthy or lesioned isocortex or in hippocampus. In hippocampus, only precursor cells with hippocampal molecular identity were able to extend projections, contacting CA3. Conversely, isocortical-like cells were capable of extending long-range axonal projections only when transplanted in motor cortex, sending fibers toward both intra- and extra-cortical targets. Ischemic damage induced by photothrombosis greatly enhanced the capability of isocortical-like cells to extend far-reaching projections. Our results indicate that neural precursors generated by ESCs carry intrinsic signals specifying axonal extension in different environments.

•

Wnt signaling induces hippocampal fate in neuralized mouse ESCs

•

Transplanted cortical and hippocampal neurons target distinct regions in adult brain

•

Photothrombotic lesion favors neurite elongation of cortical transplanted cells

•

Cortical cell transplantation improves the motor performance after ischemic damage

MGA, L3MBTL2 and E2F6 determine genomic binding of the non-canonical Polycomb repressive complex PRC1.6

Diverse Polycomb repressive complexes 1 (PRC1) play essential roles in gene regulation, differentiation and development. Six major groups of PRC1 complexes that differ in their subunit composition have been identified in mammals. How the different PRC1 complexes are recruited to specific genomic sites is poorly understood. The Polycomb Ring finger protein PCGF6, the transcription factors MGA and E2F6, and the histone-binding protein L3MBTL2 are specific components of the non-canonical PRC1.6 complex. In this study, we have investigated their role in genomic targeting of PRC1.6. ChIP-seq analysis revealed colocalization of MGA, L3MBTL2, E2F6 and PCGF6 genome-wide. Ablation of MGA in a human cell line by CRISPR/Cas resulted in complete loss of PRC1.6 binding. Rescue experiments revealed that MGA recruits PRC1.6 to specific loci both by DNA binding-dependent and by DNA binding-independent mechanisms. Depletion of L3MBTL2 and E2F6 but not of PCGF6 resulted in differential, locus-specific loss of PRC1.6 binding illustrating that different subunits mediate PRC1.6 loading to distinct sets of promoters. Mga, L3mbtl2 and Pcgf6 colocalize also in mouse embryonic stem cells, where PRC1.6 has been linked to repression of germ cell-related genes. Our findings unveil strikingly different genomic recruitment mechanisms of the non-canonical PRC1.6 complex, which specify its cell type- and context-specific regulatory functions.

Polycomb group proteins assemble in two major repressive multi-subunit complexes (PRC1 and PRC2), which play important roles in many physiological processes, including stem cell maintenance, differentiation, cell cycle control and cancer. In mammals, six different groups of PRC1 complexes exist (PRC1.1 to PRC1.6), which differ in their subunit composition. The mechanisms that target the different PRC1 complexes to specific genomic sites appear diverse and are poorly understood. In this study, we have investigated the genomic targeting mechanisms of the non-canonical PRC1.6 complex. In PRC1.6, the defining subunit PCGF6 is specifically associated with several proteins including the transcription factors MGA and E2F6, and the histone-binding protein L3MBTL2. We found that MGA is absolutely essential for targeting PRC1.6. MGA executes recruitment of PRC1.6 to its target sites through two distinct functions. On the one hand it acts as a sequence-specific DNA-binding factor; on the other hand it has a scaffolding function, which is independent of its DNA binding capacity. E2F6 and L3MBTL2 are also important in genomic targeting of PRC1.6 as they promote binding of PRC1.6 to different sets of genes associated with distinct functions. Our finding that different components specify loading of PRC1.6 to distinct sets of genes could establish a paradigm for other chromatin-associated complexes.

Label-free and real-time monitoring of single cell attachment on template-stripped plasmonic nano-holes

Leveraging microfluidics and nano-plasmonics, we present in this paper a new method employing a micro-nano-device that is capable of monitoring the dynamic cell-substrate attachment process at single cell level in real time without labeling. The micro-nano-device essentially has a gold thin film as the substrate perforated with periodic, near-cm2-area, template-stripped nano-holes, which generate plasmonic extraordinary optical transmission (EOT) with a high sensitivity to refractive index changes at the metal-dielectric interface. Using this device, we successfully demonstrated label-free and real-time monitoring of the dynamic cell attachment process for single mouse embryonic stem cell (C3H10) and human tumor cell (HeLa) by collecting EOT spectrum data during 3-hour on-chip culture. We further collected the EOT spectral shift data at the start and end points of measurement during 3-hour on-chip culture for 50 C3H10 and 50 HeLa cells, respectively. The experiment results show that the single cell attachment process of both HeLa and C3H10 cells follow the logistic retarded growth model, but with different kinetic parameters. Variations in spectral shift during the same culture period across single cells present new evidence for cell heterogeneity. The micro-nano-device provides a new, label-free, real-time, and sensitive, platform to investigate the cell adhesion kinetics at single cell level.

A conceptual and computational framework for modelling and understanding the non-equilibrium gene regulatory networks of mouse embryonic stem cells

The capacity of pluripotent embryonic stem cells to differentiate into any cell type in the body makes them invaluable in the field of regenerative medicine. However, because of the complexity of both the core pluripotency network and the process of cell fate computation it is not yet possible to control the fate of stem cells. We present a theoretical model of stem cell fate computation that is based on Halley and Winkler’s Branching Process Theory (BPT) and on Greaves et al.’s agent-based computer simulation derived from that theoretical model. BPT abstracts the complex production and action of a Transcription Factor (TF) into a single critical branching process that may dissipate, maintain, or become supercritical. Here we take the single TF model and extend it to multiple interacting TFs, and build an agent-based simulation of multiple TFs to investigate the dynamics of such coupled systems. We have developed the simulation and the theoretical model together, in an iterative manner, with the aim of obtaining a deeper understanding of stem cell fate computation, in order to influence experimental efforts, which may in turn influence the outcome of cellular differentiation. The model used is an example of self-organization and could be more widely applicable to the modelling of other complex systems. The simulation based on this model, though currently limited in scope in terms of the biology it represents, supports the utility of the Halley and Winkler branching process model in describing the behaviour of stem cell gene regulatory networks. Our simulation demonstrates three key features: (i) the existence of a critical value of the branching process parameter, dependent on the details of the cistrome in question; (ii) the ability of an active cistrome to “ignite” an otherwise fully dissipated cistrome, and drive it to criticality; (iii) how coupling cistromes together can reduce their critical branching parameter values needed to drive them to criticality.

Pluripotent stem cells possess the capacity both to renew themselves indefinitely and to differentiate to any cell type in the body. Thus the ability to direct stem cell differentiation would have immense potential in regenerative medicine. There is a massive amount of biological data relevant to stem cells; here we exploit data relating to stem cell differentiation to help understand cell behaviour and complexity. These cells contain a dynamic, non-equilibrium network of genes regulated in part by transcription factors expressed by the network itself. Here we take an existing theoretical framework, Transcription Factor Branching Processes, which explains how these genetic networks can have critical behaviour, and can tip between low and full expression. We use this theory as the basis for the design and implementation of a computational simulation platform, which we then use to run a variety of simulation experiments, to gain a better understanding how these various transcription factors can combine, interact, and influence each other. The simulation parameters are derived from experimental data relating to the core factors in pluripotent stem cell differentiation. The simulation results determine the critical values of branching process parameters, and how these are modulated by the various interacting transcription factors.

A lncRNA fine tunes the dynamics of a cell state transition involving Lin28, let-7 and de novo DNA methylation

Execution of pluripotency requires progression from the naïve status represented by mouse embryonic stem cells (ESCs) to a state capacitated for lineage specification. This transition is coordinated at multiple levels. Non-coding RNAs may contribute to this regulatory orchestra. We identified a rodent-specific long non-coding RNA (lncRNA) linc1281, hereafter Ephemeron (Eprn), that modulates the dynamics of exit from naïve pluripotency. Eprn deletion delays the extinction of ESC identity, an effect associated with perduring Nanog expression. In the absence of Eprn, Lin28a expression is reduced which results in persistence of let-7 microRNAs, and the up-regulation of de novo methyltransferases Dnmt3a/b is delayed. Dnmt3a/b deletion retards ES cell transition, correlating with delayed Nanog promoter methylation and phenocopying loss of Eprn or Lin28a. The connection from lncRNA to miRNA and DNA methylation facilitates the acute extinction of naïve pluripotency, a pre-requisite for rapid progression from preimplantation epiblast to gastrulation in rodents. Eprn illustrates how lncRNAs may introduce species-specific network modulations.

Mobilization of LINE-1 retrotransposons is restricted by Tex19.1 in mouse embryonic stem cells

Mobilization of retrotransposons to new genomic locations is a significant driver of mammalian genome evolution, but these mutagenic events can also cause genetic disorders. In humans, retrotransposon mobilization is mediated primarily by proteins encoded by LINE-1 (L1) retrotransposons, which mobilize in pluripotent cells early in development. Here we show that TEX19.1, which is induced by developmentally programmed DNA hypomethylation, can directly interact with the L1-encoded protein L1-ORF1p, stimulate its polyubiquitylation and degradation, and restrict L1 mobilization. We also show that TEX19.1 likely acts, at least in part, through promoting the activity of the E3 ubiquitin ligase UBR2 towards L1-ORF1p. Moreover, loss of Tex19.1 increases L1-ORF1p levels and L1 mobilization in pluripotent mouse embryonic stem cells, implying that Tex19.1 prevents de novo retrotransposition in the pluripotent phase of the germline cycle. These data show that post-translational regulation of L1 retrotransposons plays a key role in maintaining trans-generational genome stability in mammals.

AIRE is a critical spindle-associated protein in embryonic stem cells

Embryonic stem (ES) cells go though embryo-like cell cycles regulated by specialized molecular mechanisms. However, it is not known whether there are ES cell-specific mechanisms regulating mitotic fidelity. Here we showed that Autoimmune Regulator (Aire), a transcription coordinator involved in immune tolerance processes, is a critical spindle-associated protein in mouse ES(mES) cells. BioID analysis showed that AIRE associates with spindle-associated proteins in mES cells. Loss of function analysis revealed that Aire was important for centrosome number regulation and spindle pole integrity specifically in mES cells. We also identified the c-terminal LESLL motif as a critical motif for AIRE's mitotic function. Combined maternal and zygotic knockout further revealed Aire's critical functions for spindle assembly in preimplantation embryos. These results uncovered a previously unappreciated function for Aire and provide new insights into the biology of stem cell proliferation and potential new angles to understand fertility defects in humans carrying Aire mutations.

Differentiating Mouse Embryonic Stem Cells into Embryoid Bodies in AggreWell Plates

Embryonic stem (ES) cells can develop into many types of differentiated tissues if they are placed into a differentiating environment. This can occur in vivo when the ES cells are injected into or aggregated with an embryo, or in vitro if their culture conditions are modified to induce differentiation. Many times ES cell differentiation proceeds through an intermediate stage called the embryoid body (EB). EBs are round structures composed of ES cells that have undergone some of the initial stages of differentiation. EBs can then be manipulated further to generate more specific cell types. The method described here makes use of commercially available AggreWell 400 plates with prefabricated indentations that cradle each EB. Although these plates are relatively expensive, they may be suitable for some high-throughput experiments. This protocol describes the preparation of embryoid bodies of defined size and shape in a large-scale format (approximately 1200 embryoid bodies per preparation).

Rbm46 regulates mouse embryonic stem cell differentiation by targeting β-Catenin mRNA for degradation

Embryonic stem cells (ESCs) are pluripotent cells and have the capability for differentiation into any of the three embryonic germ layers. The Wnt/β-Catenin pathway has been shown to play an essential role in ESC differentiation regulation. Activation of β-Catenin by post-translational modification has been extensively studied. However, mechanism(s) of post-transcriptional regulation of β-Catenin are not well defined. In this study, we report an RNA recognition motif-containing protein (RNA binding motif protein 46, RBM46) which regulates the degradation of β-Catenin mRNA. Our results show that Rbm46 is distributed primarily in the cytoplasm of mouse ESCs (mESCs) and is elevated during the process of ESC differentiation. In addition, overexpression of Rbm46 results in differentiation of mESCs into trophectoderm, while knock-down of Rbm46 leads to mESC differentiation into endoderm. β-Catenin, a key effector in the Wnt pathway which has been reported to play a significant role in the regulation of ESC differentiation, is post-transcriptionally regulated by Rbm46. Our study reveals Rbm46 plays a novel role in the regulation of ESC differentiation.

Involvement of Wnt pathway in ethanol-induced inhibition of mouse embryonic stem cell differentiation

Ethanol has been reported to have toxicity on embryonic stem cells (ESCs). The present study aims to address the teratogenic effects of ethanol on the growth and cardiac differentiation of ESCs. Mouse embryonic stem D3 cells were employed. 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) and lactate dehydrogenase (LDH) assays were used to determine cytotoxicity. Quantitative real time polymerase chain reaction (qRT-PCR) and Western blotting were used to analyze the expressions of cardiac differentiation-related and Wnt signaling factors. The beating profile of cardiomyocytes was recorded to assess cardiac differentiation. Ethanol induced growth inhibition in both undifferentiated and differentiated ESCs after 5 days of exposure. Ethanol inhibited the loss of pluripotent gene expressions including Nanog, Sox2 and Oct4. The expressions of cardiac markers, Nkx2.5, Mef2c, Tbx5, dHand, αMHC, Cx43 and troponin C1, were suppressed by ethanol treatment. Furthermore, ethanol delayed cardiac differentiation of ESCs till 11 days of differentiation. The expressions of Wnt-related regulators, β-catenin and its target cyclin D1, were downregulated by ethanol. Wnt pathway agonist wnt3a could greatly rescue ethanol-induced inhibition of cardiac differentiation and Wnt-pathway-related protein expressions. These finding suggested that ethanol suppresses mouse ESC differentiation largely by inhibiting Wnt signaling pathway.

Single Cell RNA-Sequencing of Pluripotent States Unlocks Modular Transcriptional Variation

Embryonic stem cell (ESC) culture conditions are important for maintaining long-term self-renewal, and they influence cellular pluripotency state. Here, we report single cell RNA-sequencing of mESCs cultured in three different conditions: serum, 2i, and the alternative ground state a2i. We find that the cellular transcriptomes of cells grown in these conditions are distinct, with 2i being the most similar to blastocyst cells and including a subpopulation resembling the two-cell embryo state. Overall levels of intercellular gene expression heterogeneity are comparable across the three conditions. However, this masks variable expression of pluripotency genes in serum cells and homogeneous expression in 2i and a2i cells. Additionally, genes related to the cell cycle are more variably expressed in the 2i and a2i conditions. Mining of our dataset for correlations in gene expression allowed us to identify additional components of the pluripotency network, including Ptma and Zfp640, illustrating its value as a resource for future discovery.

PAPST, a User Friendly and Powerful Java Platform for ChIP-Seq Peak Co-Localization Analysis and Beyond

Comparative co-localization analysis of transcription factors (TFs) and epigenetic marks (EMs) in specific biological contexts is one of the most critical areas of ChIP-Seq data analysis beyond peak calling. Yet there is a significant lack of user-friendly and powerful tools geared towards co-localization analysis based exploratory research. Most tools currently used for co-localization analysis are command line only and require extensive installation procedures and Linux expertise. Online tools partially address the usability issues of command line tools, but slow response times and few customization features make them unsuitable for rapid data-driven interactive exploratory research. We have developed PAPST: Peak Assignment and Profile Search Tool, a user-friendly yet powerful platform with a unique design, which integrates both gene-centric and peak-centric co-localization analysis into a single package. Most of PAPST’s functions can be completed in less than five seconds, allowing quick cycles of data-driven hypothesis generation and testing. With PAPST, a researcher with or without computational expertise can perform sophisticated co-localization pattern analysis of multiple TFs and EMs, either against all known genes or a set of genomic regions obtained from public repositories or prior analysis. PAPST is a versatile, efficient, and customizable tool for genome-wide data-driven exploratory research. Creatively used, PAPST can be quickly applied to any genomic data analysis that involves a comparison of two or more sets of genomic coordinate intervals, making it a powerful tool for a wide range of exploratory genomic research. We first present PAPST’s general purpose features then apply it to several public ChIP-Seq data sets to demonstrate its rapid execution and potential for cutting-edge research with a case study in enhancer analysis. To our knowledge, PAPST is the first software of its kind to provide efficient and sophisticated post peak-calling ChIP-Seq data analysis as an easy-to-use interactive application. PAPST is available at https://github.com/paulbible/papst and is a public domain work.

Oxygen transport and stem cell aggregation in stirred-suspension bioreactor cultures

Stirred-suspension bioreactors are a promising modality for large-scale culture of 3D aggregates of pluripotent stem cells and their progeny. Yet, cells within these clusters experience limitations in the transfer of factors and particularly O₂ which is characterized by low solubility in aqueous media. Cultured stem cells under different O₂ levels may exhibit significantly different proliferation, viability and differentiation potential. Here, a transient diffusion-reaction model was built encompassing the size distribution and ultrastructural characteristics of embryonic stem cell (ESC) aggregates. The model was coupled to experimental data from bioreactor and static cultures for extracting the effective diffusivity and kinetics of consumption of O₂ within mouse (mESC) and human ESC (hESC) clusters. Under agitation, mESC aggregates exhibited a higher maximum consumption rate than hESC aggregates. Moreover, the reaction-diffusion model was integrated with a population balance equation (PBE) for the temporal distribution of ESC clusters changing due to aggregation and cell proliferation. Hypoxia was found to be negligible for ESCs with a smaller radius than 100 µm but became appreciable for aggregates larger than 300 µm. The integrated model not only captured the O₂ profile both in the bioreactor bulk and inside ESC aggregates but also led to the calculation of the duration that fractions of cells experience a certain range of O₂ concentrations. The approach described in this study can be employed for gaining a deeper understanding of the effects of O₂ on the physiology of stem cells organized in 3D structures. Such frameworks can be extended to encompass the spatial and temporal availability of nutrients and differentiation factors and facilitate the design and control of relevant bioprocesses for the production of stem cell therapeutics.