Multi-omic Profiling Reveals Dynamics of the Phased Progression of Pluripotency

The transcription factor OCT6 promotes the dissolution of the naïve pluripotent state by repressing Nanog and activating a formative state gene regulatory network

In the mouse embryo, the transition from the preimplantation to the postimplantation epiblast is governed by changes in the gene regulatory network (GRN) that lead to transcriptional, epigenetic, and functional changes. This transition can be faithfully recapitulated in vitro by the differentiation of mouse embryonic stem cells (mESCs) to epiblast-like cells (EpiLCs), that reside in naïve and formative states of pluripotency, respectively. However, the GRN that drives this conversion is not fully elucidated. Here we demonstrate that the transcription factor OCT6 is a key driver of this process. Firstly, we show that Oct6 is not expressed in mESCs but is rapidly induced as cells exit the naïve pluripotent state. By deleting Oct6 in mESCs, we find that knockout cells fail to acquire the typical morphological changes associated with the formative state when induced to differentiate. Additionally, the key naïve pluripotency TFs Nanog, Klf2, Nr5a2, Prdm14, and Esrrb were expressed at higher levels than in wild-type cells, indicating an incomplete dismantling of the naïve pluripotency GRN. Conversely, premature expression of Oct6 in naïve cells triggered a rapid morphological transformation mirroring differentiation, that was accompanied by the upregulation of the endogenous Oct6 as well as the formative genes Sox3, Zic2/3, Foxp1, Dnmt3A and FGF5. Strikingly, we found that OCT6 represses Nanog in a bistable manner and that this regulation is at the transcriptional level. Moreover, our findings also reveal that Oct6 is repressed by NANOG. Collectively, our results establish OCT6 as a key TF in the dissolution of the naïve pluripotent state and support a model where Oct6 and Nanog form a double negative feedback loop which could act as an important toggle mediating the transition to the formative state.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

ZFP982 confers mouse embryonic stem cell characteristics by regulating expression of Nanog, Zfp42, and Dppa3

BACKGROUND

Understanding the genetic underpinnings of protein networks conferring stemness is of broad interest for basic and translational research.

METHODS

We used multi-omics analyses to identify and characterize stemness genes, and focused on the zinc finger protein 982 (Zfp982) that regulates stemness through the expression of Nanog, Zfp42, and Dppa3 in mouse embryonic stem cells (mESC).

RESULTS

Zfp982 was expressed in stem cells, and bound to chromatin through a GCAGAGKC motif, for example near the stemness genes Nanog, Zfp42, and Dppa3. Nanog and Zfp42 were direct targets of ZFP982 that decreased in expression upon knockdown and increased upon overexpression of Zfp982. We show that ZFP982 expression strongly correlated with stem cell characteristics, both on the transcriptional and morphological levels. Zfp982 expression decreased with progressive differentiation into ecto-, endo- and mesodermal cell lineages, and knockdown of Zfp982 correlated with morphological and transcriptional features of differentiated cells. Zfp982 showed transcriptional overlap with members of the Hippo signaling pathway, one of which was Yap1, the major co-activator of Hippo signaling. Despite the observation that ZFP982 and YAP1 interacted and localized predominantly to the cytoplasm upon differentiation, the localization of YAP1 was not influenced by ZFP982 localization.

CONCLUSIONS

Together, our study identified ZFP982 as a transcriptional regulator of early stemness genes, and since ZFP982 is under the control of the Hippo pathway, underscored the importance of the context-dependent Hippo signals for stem cell characteristics.

KDM3A and KDM3B Maintain Naïve Pluripotency Through the Regulation of Alternative Splicing

Histone modifying enzymes play a central role in maintaining cell identity by establishing a conducive chromatin environment for lineage specific transcription factor activity. Pluripotent embryonic stem cell (ESC) identity is characterized by a lower abundance of gene repression associated histone modifications that enables rapid response to differentiation cues. The KDM3 family of histone demethylases removes the repressive histone H3 lysine 9 dimethylation (H3K9me2). Here we uncover a surprising role for the KDM3 proteins in the maintenance of the pluripotent state through post-transcriptional regulation. We find that KDM3A and KDM3B interact with RNA processing factors such as EFTUD2 and PRMT5. Acute selective degradation of the endogenous KDM3A and KDM3B proteins resulted in altered splicing independent of H3K9me2 status or catalytic activity. These splicing changes partially resemble the splicing pattern of the more blastocyst-like ground state of pluripotency and occurred in important chromatin and transcription factors such as Dnmt3b, Tbx3 and Tcf12. Our findings reveal non-canonical roles of histone demethylating enzymes in splicing to regulate cell identity.

Impact of eIF2α phosphorylation on the translational landscape of mouse embryonic stem cells

The integrated stress response (ISR) is critical for cell survival under stress. In response to diverse environmental cues, eIF2α becomes phosphorylated, engendering a dramatic change in mRNA translation. The activation of ISR plays a pivotal role in the early embryogenesis, but the eIF2-dependent translational landscape in pluripotent embryonic stem cells (ESCs) is largely unexplored. We employ a multi-omics approach consisting of ribosome profiling, proteomics, and metabolomics in wild-type (eIF2α+/+) and phosphorylation-deficient mutant eIF2α (eIF2αA/A) mouse ESCs (mESCs) to investigate phosphorylated (p)-eIF2α-dependent translational control of naive pluripotency. We show a transient increase in p-eIF2α in the naive epiblast layer of E4.5 embryos. Absence of eIF2α phosphorylation engenders an exit from naive pluripotency following 2i (two chemical inhibitors of MEK1/2 and GSK3α/β) withdrawal. p-eIF2α controls translation of mRNAs encoding proteins that govern pluripotency, chromatin organization, and glutathione synthesis. Thus, p-eIF2α acts as a key regulator of the naive pluripotency gene regulatory network.

MLL3/MLL4 enzymatic activity shapes DNA replication timing

Mammalian genomes are replicated in a precise order during S phase, which is cell-type-specific1-3 and correlates with local transcriptional activity2,4-8, chromatin modifications9 and chromatin architecture1,10,11,12. However, the causal relationships between these features and the key regulators of DNA replication timing (RT) are largely unknown. Here, machine learning was applied to quantify chromatin features, including epigenetic marks, histone variants and chromatin architectural factors, best predicting local RT under steady-state and RT changes during early embryonic stem (ES) cell differentiation. About one-third of genome exhibited RT changes during the differentiation. Combined, chromatin features predicted steady-state RT and RT changes with high accuracy. Of these features, histone H3 lysine 4 monomethylation (H3K4me1) catalyzed by MLL3/4 (also known as KMT2C/D) emerged as a top predictor. Loss of Mll3/4 (but not Mll3 alone) or their enzymatic activity resulted in erasure of genome-wide RT dynamics during ES cell differentiation. Sites that normally gain H3K4me1 in a MLL3/4-dependent fashion during the transition failed to transition towards earlier RT, often with transcriptional activation unaffected. Further analysis revealed a requirement for MLL3/4 in promoting DNA replication initiation zones through MCM2 recruitment, providing a direct link for its role in regulating RT. Our results uncover MLL3/4-dependent H3K4me1 as a functional regulator of RT and highlight a causal relationship between the epigenome and RT that is largely uncoupled from transcription. These findings uncover a previously unknown role for MLL3/4-dependent chromatin functions which is likely relevant to the numerous diseases associated with MLL3/4 mutations.

SnapKin: a snapshot deep learning ensemble for kinase-substrate prediction from phosphoproteomics data

A major challenge in mass spectrometry-based phosphoproteomics lies in identifying the substrates of kinases, as currently only a small fraction of substrates identified can be confidently linked with a known kinase. Machine learning techniques are promising approaches for leveraging large-scale phosphoproteomics data to computationally predict substrates of kinases. However, the small number of experimentally validated kinase substrates (true positive) and the high data noise in many phosphoproteomics datasets together limit their applicability and utility. Here, we aim to develop advanced kinase-substrate prediction methods to address these challenges. Using a collection of seven large phosphoproteomics datasets, and both traditional and deep learning models, we first demonstrate that a 'pseudo-positive' learning strategy for alleviating small sample size is effective at improving model predictive performance. We next show that a data resampling-based ensemble learning strategy is useful for improving model stability while further enhancing prediction. Lastly, we introduce an ensemble deep learning model ('SnapKin') by incorporating the above two learning strategies into a 'snapshot' ensemble learning algorithm. We propose SnapKin, an ensemble deep learning method, for predicting substrates of kinases from large-scale phosphoproteomics data. We demonstrate that SnapKin consistently outperforms existing methods in kinase-substrate prediction. SnapKin is freely available at https://github.com/PYangLab/SnapKin.

AKT1 induces Nanog promoter in a SUMOylation-dependent manner in different pluripotent contexts

AKT/PKB is a kinase crucial for pluripotency maintenance in pluripotent stem cells. Multiple post-translational modifications modulate its activity. We have previously demonstrated that AKT1 induces the expression of the pluripotency transcription factor Nanog in a SUMOylation-dependent manner in mouse embryonic stem cells. Here, we studied different cellular contexts and main candidates that could mediate this induction. Our results strongly suggest the pluripotency transcription factors OCT4 and SOX2 are not essential mediators. Additionally, we concluded that this induction takes place in different pluripotent contexts but not in terminally differentiated cells. Finally, the cross-matching analysis of ESCs, iPSCs and MEFs transcriptomes and AKT1 phosphorylation targets provided new clues about possible factors that could be involved in the SUMOylation-dependent Nanog induction by AKT.

Homocysteine impedes neurite outgrowth recovery after intracerebral haemorrhage by downregulating pCAMK2A

Hyperhomocysteinemia (HHcy) is independently associated with poorer long-term prognosis in patients with intracerebral haemorrhage (ICH); however, the effect and mechanisms of HHcy on ICH are still unclear. Here, we evaluated neurite outgrowth and neurological functional recovery using simulated models of ICH with HHcy in vitro and in vivo. We found that the neurite outgrowth velocity and motor functional recovery in the ICH plus HHcy group were significantly slower than that in the control group, indicating that homocysteine (Hcy) significantly impedes the neurite outgrowth recovery after ICH. Furthermore, phosphoproteomic data and signalome analysis of perihematomal brain tissues suggested that calmodulin-dependent protein kinases 2 (CAMK2A) kinase substrate pairs were significantly downregulated in ICH with HHcy compared with autologous blood injection only, both western blot and immunofluorescence staining confirmed this finding. Additionally, upregulation of pCAMK2A significantly increased neurite outgrowth recovery in ICH with HHcy. Collectively, we clarify the mechanism of HHcy-hindered neurite outgrowth recovery, and pCAMK2A may serve as a therapeutic strategy for promoting neurological recovery after ICH.

Parallel Genome-Wide CRISPR Screens to Identify State-Dependent Self-Renewal Regulators of Mouse Embryonic Stem Cells

The pluripotency of embryonic stem cells (ESCs) is more accurately viewed as a continuous developmental process rather than a fixed state. However, the factors that play general or state-specific roles in regulating self-renewal in different pluripotency states remain poorly defined. In this study, parallel genome-wide CRISPR/Cas9 knockout (KO) screens were applied in ESCs cultured in the serum plus LIF (SL) and in the 2i plus LIF (2iL) conditions. The candidate genes were classified into seven groups based on their positive or negative effects on self-renewal, and whether this effect was general or state-specific for ESCs under SL and 2iL culture conditions. We characterized the expression and function of genes in these seven groups. The loss of function of novel pluripotent candidate genes Usp28, Zfp598, and Zfp296 was further evaluated in mouse ESCs. Consistent with our screen, the knockout of Usp28 promotes the proliferation of SL-ESCs and 2iL-ESCs, whereas Zfp598 is indispensable for the self-renewal of ESCs under both culture conditions. The cell phenotypes of Zfp296 KO ESCs under SL and 2iL culture conditions were different. Our work provided a valuable resource for dissecting the molecular regulation of ESC self-renewal in different pluripotency states.

Stable maternal proteins underlie distinct transcriptome, translatome, and proteome reprogramming during mouse oocyte-to-embryo transition

BACKGROUND

The oocyte-to-embryo transition (OET) converts terminally differentiated gametes into a totipotent embryo and is critically controlled by maternal mRNAs and proteins, while the genome is silent until zygotic genome activation. How the transcriptome, translatome, and proteome are coordinated during this critical developmental window remains poorly understood.

RESULTS

Utilizing a highly sensitive and quantitative mass spectrometry approach, we obtain high-quality proteome data spanning seven mouse stages, from full-grown oocyte (FGO) to blastocyst, using 100 oocytes/embryos at each stage. Integrative analyses reveal distinct proteome reprogramming compared to that of the transcriptome or translatome. FGO to 8-cell proteomes are dominated by FGO-stockpiled proteins, while the transcriptome and translatome are more dynamic. FGO-originated proteins frequently persist to blastocyst while corresponding transcripts are already downregulated or decayed. Improved concordance between protein and translation or transcription is observed for genes starting translation upon meiotic resumption, as well as those transcribed and translated only in embryos. Concordance between protein and transcription/translation is also observed for proteins with short half-lives. We built a kinetic model that predicts protein dynamics by incorporating both initial protein abundance in FGOs and translation kinetics across developmental stages.

CONCLUSIONS

Through integrative analyses of datasets generated by ultrasensitive methods, our study reveals that the proteome shows distinct dynamics compared to the translatome and transcriptome during mouse OET. We propose that the remarkably stable oocyte-originated proteome may help save resources to accommodate the demanding needs of growing embryos. This study will advance our understanding of mammalian OET and the fundamental principles governing gene expression.

Human amnion-derived mesenchymal stem cells attenuate acute lung injury in two different acute lung injury mice models

Acute lung injury (ALI) is one of the most common clinical emergencies with limited effective pharmaceutical treatment in the clinic, especially when it progresses to acute respiratory distress syndrome (ARDS). Currently, mesenchymal stem cells (MSCs) exhibit specific superiority for ALI/ARDS treatment. However, stem cells from different sources may result in controversial effects on similar disease conditions. This study aimed to determine the effects of human amnion-derived mesenchymal stem cells (hAMSCs) on two different ALI mice model. The administered hAMSCs effectively accumulated in the lung tissues in all hAMSC-treated groups. Compared with the model and 1% human serum albumin (HSA) groups, high-dose hAMSCs (1.0 × 10⁶ cells) group significantly alleviated alveolar-capillary permeability, oxidative stress, inflammatory factors level and histopathological damage. In addition, the NF-κB signaling pathway is one of the key pathways activated during lipopolysaccharide (LPS) or paraquat (PQ)-induced lung injury. Our results indicated that hAMSCs (1.0 × 10⁶ cells) obviously inhibited the expression of p-IKKα/β, p-IκBα, and p-p65 in the lung tissue (p < 0.05). The high-dose (HD) hAMSC treatment exerted beneficial therapeutic effects on ALI mice models without detectable adverse reactions. The therapeutic effect of hAMSCs might involve NF-κB signaling pathway inhibition. hAMSC treatment is a potential candidate therapy for ALI.

PAD2: interactive exploration of transcription factor genomic colocalization using ChIP-seq data

Characterizing transcription factor (TF) genomic colocalization is essential for identifying cooperative binding of TFs in controlling gene expression. Here, we introduce a protocol for using PAD2, an interactive web application that enables the investigation of colocalization of various TFs and chromatin-regulating proteins from mouse embryonic stem cells at various functional genomic regions. We describe steps for accessing and searching the PAD2 database and selecting and submitting genomic regions. We then detail protein colocalization analysis using heatmap and ranked correlation plot. For complete details on the use and execution of this protocol, please refer to Kim et al. (2022).¹.

Optogenetic manipulation identifies the roles of ERK and AKT dynamics in controlling mouse embryonic stem cell exit from pluripotency

ERK and AKT signaling control pluripotent cell self-renewal versus differentiation. ERK pathway activity over time (i.e., dynamics) is heterogeneous between individual pluripotent cells, even in response to the same stimuli. To analyze potential functions of ERK and AKT dynamics in controlling mouse embryonic stem cell (ESC) fates, we developed ESC lines and experimental pipelines for the simultaneous long-term manipulation and quantification of ERK or AKT dynamics and cell fates. We show that ERK activity duration or amplitude or the type of ERK dynamics (e.g., transient, sustained, or oscillatory) alone does not influence exit from pluripotency, but the sum of activity over time does. Interestingly, cells retain memory of previous ERK pulses, with duration of memory retention dependent on duration of previous pulse length. FGF receptor/AKT dynamics counteract ERK-induced pluripotency exit. These findings improve our understanding of how cells integrate dynamics from multiple signaling pathways and translate them into cell fate cues.

The dynamical organization of the core pluripotency transcription factors responds to differentiation cues in early S-phase

DNA replication in stem cells is a major challenge for pluripotency preservation and cell fate decisions. This process involves massive changes in the chromatin architecture and the reorganization of many transcription-related molecules in different spatial and temporal scales. Pluripotency is controlled by the master transcription factors (TFs) OCT4, SOX2 and NANOG that partition into condensates in the nucleus of embryonic stem cells. These condensates are proposed to play relevant roles in the regulation of gene expression and the maintenance of pluripotency. Here, we asked whether the dynamical distribution of the pluripotency TFs changes during the cell cycle, particularly during DNA replication. Since the S phase is considered to be a window of opportunity for cell fate decisions, we explored if differentiation cues in G1 phase trigger changes in the distribution of these TFs during the subsequent S phase. Our results show a spatial redistribution of TFs condensates during DNA replication which was not directly related to chromatin compaction. Additionally, fluorescence fluctuation spectroscopy revealed TF-specific, subtle changes in the landscape of TF-chromatin interactions, consistent with their particularities as key players of the pluripotency network. Moreover, we found that differentiation stimuli in the preceding G1 phase triggered a relatively fast and massive reorganization of pluripotency TFs in early-S phase. Particularly, OCT4 and SOX2 condensates dissolved whereas the lifetimes of TF-chromatin interactions increased suggesting that the reorganization of condensates is accompanied with a change in the landscape of TF-chromatin interactions. Notably, NANOG showed impaired interactions with chromatin in stimulated early-S cells in line with its role as naïve pluripotency TF. Together, these findings provide new insights into the regulation of the core pluripotency TFs during DNA replication of embryonic stem cells and highlight their different roles at early differentiation stages.

Esrrb guides naive pluripotent cells through the formative transcriptional programme

During embryonic development, naive pluripotent epiblast cells transit to a formative state. The formative epiblast cells form a polarized epithelium, exhibit distinct transcriptional and epigenetic profiles and acquire competence to differentiate into all somatic and germline lineages. However, we have limited understanding of how the transition to a formative state is molecularly controlled. Here we used murine embryonic stem cell models to show that ESRRB is both required and sufficient to activate formative genes. Genetic inactivation of Esrrb leads to illegitimate expression of mesendoderm and extra-embryonic markers, impaired formative expression and failure to self-organize in 3D. Functionally, this results in impaired ability to generate formative stem cells and primordial germ cells in the absence of Esrrb. Computational modelling and genomic analyses revealed that ESRRB occupies key formative genes in naive cells and throughout the formative state. In so doing, ESRRB kickstarts the formative transition, leading to timely and unbiased capacity for multi-lineage differentiation.

PRC2.1- and PRC2.2-specific accessory proteins drive recruitment of different forms of canonical PRC1

Polycomb repressive complex 2 (PRC2) mediates H3K27me3 deposition, which is thought to recruit canonical PRC1 (cPRC1) via chromodomain-containing CBX proteins to promote stable repression of developmental genes. PRC2 forms two major subcomplexes, PRC2.1 and PRC2.2, but their specific roles remain unclear. Through genetic knockout (KO) and replacement of PRC2 subcomplex-specific subunits in naïve and primed pluripotent cells, we uncover distinct roles for PRC2.1 and PRC2.2 in mediating the recruitment of different forms of cPRC1. PRC2.1 catalyzes the majority of H3K27me3 at Polycomb target genes and is sufficient to promote recruitment of CBX2/4-cPRC1 but not CBX7-cPRC1. Conversely, while PRC2.2 is poor at catalyzing H3K27me3, we find that its accessory protein JARID2 is essential for recruitment of CBX7-cPRC1 and the consequent 3D chromatin interactions at Polycomb target genes. We therefore define distinct contributions of PRC2.1- and PRC2.2-specific accessory proteins to Polycomb-mediated repression and uncover a new mechanism for cPRC1 recruitment.

Loss of MLL3/4 decouples enhancer H3K4 monomethylation, H3K27 acetylation, and gene activation during embryonic stem cell differentiation

BACKGROUND

Enhancers are essential in defining cell fates through the control of cell-type-specific gene expression. Enhancer activation is a multi-step process involving chromatin remodelers and histone modifiers including the monomethylation of H3K4 (H3K4me1) by MLL3 (KMT2C) and MLL4 (KMT2D). MLL3/4 are thought to be critical for enhancer activation and cognate gene expression including through the recruitment of acetyltransferases for H3K27.

RESULTS

Here we test this model by evaluating the impact of MLL3/4 loss on chromatin and transcription during early differentiation of mouse embryonic stem cells. We find that MLL3/4 activity is required at most if not all sites that gain or lose H3K4me1 but is largely dispensable at sites that remain stably methylated during this transition. This requirement extends to H3K27 acetylation (H3K27ac) at most transitional sites. However, many sites gain H3K27ac independent of MLL3/4 or H3K4me1 including enhancers regulating key factors in early differentiation. Furthermore, despite the failure to gain active histone marks at thousands of enhancers, transcriptional activation of nearby genes is largely unaffected, thus uncoupling the regulation of these chromatin events from transcriptional changes during this transition. These data challenge current models of enhancer activation and imply distinct mechanisms between stable and dynamically changing enhancers.

CONCLUSIONS

Collectively, our study highlights gaps in knowledge about the steps and epistatic relationships of enzymes necessary for enhancer activation and cognate gene transcription.

Phosphoproteomics reveals rewiring of the insulin signaling network and multi-nodal defects in insulin resistance

The failure of metabolic tissues to appropriately respond to insulin ("insulin resistance") is an early marker in the pathogenesis of type 2 diabetes. Protein phosphorylation is central to the adipocyte insulin response, but how adipocyte signaling networks are dysregulated upon insulin resistance is unknown. Here we employ phosphoproteomics to delineate insulin signal transduction in adipocyte cells and adipose tissue. Across a range of insults causing insulin resistance, we observe a marked rewiring of the insulin signaling network. This includes both attenuated insulin-responsive phosphorylation, and the emergence of phosphorylation uniquely insulin-regulated in insulin resistance. Identifying dysregulated phosphosites common to multiple insults reveals subnetworks containing non-canonical regulators of insulin action, such as MARK2/3, and causal drivers of insulin resistance. The presence of several bona fide GSK3 substrates among these phosphosites led us to establish a pipeline for identifying context-specific kinase substrates, revealing widespread dysregulation of GSK3 signaling. Pharmacological inhibition of GSK3 partially reverses insulin resistance in cells and tissue explants. These data highlight that insulin resistance is a multi-nodal signaling defect that includes dysregulated MARK2/3 and GSK3 activity.

Computational systems approach towards phosphoproteomics and their downstream regulation

Protein phosphorylation plays an essential role in modulating cell signalling and its downstream transcriptional and translational regulations. Until recently, protein phosphorylation has been studied mostly using low-throughput biochemical assays. The advancement of mass spectrometry (MS)-based phosphoproteomics transformed the field by enabling measurement of proteome-wide phosphorylation events, where tens of thousands of phosphosites are routinely identified and quantified in an experiment. This has brought a significant challenge in analysing large-scale phosphoproteomic data, making computational methods and systems approaches integral parts of phosphoproteomics. Previous works have primarily focused on reviewing the experimental techniques in MS-based phosphoproteomics, yet a systematic survey of the computational landscape in this field is still missing. Here, we review computational methods and tools, and systems approaches that have been developed for phosphoproteomics data analysis. We categorise them into four aspects including data processing, functional analysis, phosphoproteome annotation and their integration with other omics, and in each aspect, we discuss the key methods and example studies. Lastly, we highlight some of the potential research directions on which future work would make a significant contribution to this fast-growing field. We hope this review provides a useful snapshot of the field of computational systems phosphoproteomics and stimulates new research that drives future development.

bulkAnalyseR: an accessible, interactive pipeline for analysing and sharing bulk multi-modal sequencing data

Bulk sequencing experiments (single- and multi-omics) are essential for exploring wide-ranging biological questions. To facilitate interactive, exploratory tasks, coupled with the sharing of easily accessible information, we present bulkAnalyseR, a package integrating state-of-the-art approaches using an expression matrix as the starting point (pre-processing functions are available as part of the package). Static summary images are replaced with interactive panels illustrating quality-checking, differential expression analysis (with noise detection) and biological interpretation (enrichment analyses, identification of expression patterns, followed by inference and comparison of regulatory interactions). bulkAnalyseR can handle different modalities, facilitating robust integration and comparison of cis-, trans- and customised regulatory networks.

Molecular versatility during pluripotency progression

A challenge during development is to ensure lineage segregation while preserving plasticity. Using pluripotency progression as a paradigm, we review how developmental transitions are coordinated by redeployments, rather than global resettings, of cellular components. We highlight how changes in response to extrinsic cues (FGF, WNT, Activin/Nodal, Netrin-1), context- and stoichiometry-dependent action of transcription factors (Oct4, Nanog) and reconfigurations of epigenetic regulators (enhancers, promoters, TrxG, PRC) may confer robustness to naïve to primed pluripotency transition. We propose the notion of Molecular Versatility to regroup mechanisms by which molecules are repurposed to exert different, sometimes opposite, functions in close stem cell configurations.

A high-content endogenous GLUT4 trafficking assay reveals new aspects of adipocyte biology

Insulin-induced GLUT4 translocation to the plasma membrane in muscle and adipocytes is crucial for whole-body glucose homeostasis. Currently, GLUT4 trafficking assays rely on overexpression of tagged GLUT4. Here we describe a high-content imaging platform for studying endogenous GLUT4 translocation in intact adipocytes. This method enables high fidelity analysis of GLUT4 responses to specific perturbations, multiplexing of other trafficking proteins and other features including lipid droplet morphology. Using this multiplexed approach we showed that Vps45 and Rab14 are selective regulators of GLUT4, but Trarg1, Stx6, Stx16, Tbc1d4 and Rab10 knockdown affected both GLUT4 and TfR translocation. Thus, GLUT4 and TfR translocation machinery likely have some overlap upon insulin-stimulation. In addition, we identified Kif13A, a Rab10 binding molecular motor, as a novel regulator of GLUT4 traffic. Finally, comparison of endogenous to overexpressed GLUT4 highlights that the endogenous GLUT4 methodology has an enhanced sensitivity to genetic perturbations and emphasises the advantage of studying endogenous protein trafficking for drug discovery and genetic analysis of insulin action in relevant cell types.

Multi-omics profiling visualizes dynamics of cardiac development and functions

Cardiogenesis is a tightly regulated dynamic process through a continuum of differentiation and proliferation events. Key factors and pathways governing this process remain incompletely understood. Here, we investigate mice hearts from embryonic day 10.5 to postnatal week 8 and dissect developmental changes in phosphoproteome-, proteome-, metabolome-, and transcriptome-encompassing cardiogenesis and cardiac maturation. We identify mitogen-activated protein kinases as core kinases involved in transcriptional regulation by mediating the phosphorylation of chromatin remodeling proteins during early cardiogenesis. We construct the reciprocal regulatory network of transcription factors (TFs) and identify a series of TFs controlling early cardiogenesis involved in cycling-dependent proliferation. After birth, we identify cardiac resident macrophages with high arachidonic acid metabolism activities likely involved in the clearance of injured apoptotic cardiomyocytes. Together, our comprehensive multi-omics data offer a panoramic view of cardiac development and maturation that provides a resource for further in-depth functional exploration.

The Sum of Two Halves May Be Different from the Whole-Effects of Splitting Sequencing Samples Across Lanes

The advances in high-throughput sequencing (HTS) have enabled the characterisation of biological processes at an unprecedented level of detail; most hypotheses in molecular biology rely on analyses of HTS data. However, achieving increased robustness and reproducibility of results remains a main challenge. Although variability in results may be introduced at various stages, e.g., alignment, summarisation or detection of differential expression, one source of variability was systematically omitted: the sequencing design, which propagates through analyses and may introduce an additional layer of technical variation. We illustrate qualitative and quantitative differences arising from splitting samples across lanes on bulk and single-cell sequencing. For bulk mRNAseq data, we focus on differential expression and enrichment analyses; for bulk ChIPseq data, we investigate the effect on peak calling and the peaks' properties. At the single-cell level, we concentrate on identifying cell subpopulations. We rely on markers used for assigning cell identities; both smartSeq and 10× data are presented. The observed reduction in the number of unique sequenced fragments limits the level of detail on which the different prediction approaches depend. Furthermore, the sequencing stochasticity adds in a weighting bias corroborated with variable sequencing depths and (yet unexplained) sequencing bias. Subsequently, we observe an overall reduction in sequencing complexity and a distortion in the biological signal across technologies, experimental contexts, organisms and tissues.

Progression of the pluripotent epiblast depends upon the NMD factor UPF2

Nonsense-mediated RNA decay (NMD) is a highly conserved RNA turnover pathway that degrades RNAs harboring in-frame stop codons in specific contexts. Loss of NMD factors leads to embryonic lethality in organisms spanning the phylogenetic scale, but the mechanism remains unknown. Here, we report that the core NMD factor, UPF2, is required for expansion of epiblast cells within the inner cell mass of mice in vivo. We identify NMD target mRNAs in mouse blastocysts - both canonical and alternatively processed mRNAs - including those encoding cell cycle arrest and apoptosis factors, raising the possibility that NMD is essential for embryonic cell proliferation and survival. In support, the inner cell mass of Upf2-null blastocysts rapidly regresses with outgrowth and is incompetent for embryonic stem cell derivation in vitro. In addition, we uncovered concordant temporal- and lineage-specific regulation of NMD factors and mRNA targets, indicative of a shift in NMD magnitude during peri-implantation development. Together, our results reveal developmental and molecular functions of the NMD pathway in the early embryo.

Integrative analysis reveals histone demethylase LSD1 promotes RNA polymerase II pausing

Lysine-specific demethylase 1 (LSD1) is well-known for its role in decommissioning enhancers during mouse embryonic stem cell (ESC) differentiation. Its role in gene promoters remains poorly understood despite its widespread presence at these sites. Here, we report that LSD1 promotes RNA polymerase II (RNAPII) pausing, a rate-limiting step in transcription regulation, in ESCs. We found the knockdown of LSD1 preferentially affects genes with higher RNAPII pausing. Next, we demonstrate that the co-localization sites of LSD1 and MYC, a factor known to regulate pause-release, are enriched for other RNAPII pausing factors. We show that LSD1 and MYC directly interact and MYC recruitment to genes co-regulated with LSD1 is dependent on LSD1 but not vice versa. The co-regulated gene set is significantly enriched for housekeeping processes and depleted of transcription factors compared to those bound by LSD1 alone. Collectively, our integrative analysis reveals a pleiotropic role of LSD1 in promoting RNAPII pausing.

•

LSD1 promotes RNA polymerase II pausing in mouse embryonic stem cells

•

LSD1 knockdown causes global reduction of RNAPII pausing

•

Co-localized sites of LSD1 and MYC are enriched for RNAPII pausing and releasing factors

•

MYC recruitment to co-regulated genes is dependent on LSD1 but not vice versa

H3K18 lactylation marks tissue-specific active enhancers

Background

Histone lactylation has been recently described as a novel histone post-translational modification linking cellular metabolism to epigenetic regulation.

Results

Given the expected relevance of this modification and current limited knowledge of its function, we generate genome-wide datasets of H3K18la distribution in various in vitro and in vivo samples, including mouse embryonic stem cells, macrophages, adipocytes, and mouse and human skeletal muscle. We compare them to profiles of well-established histone modifications and gene expression patterns. Supervised and unsupervised bioinformatics analysis shows that global H3K18la distribution resembles H3K27ac, although we also find notable differences. H3K18la marks active CpG island-containing promoters of highly expressed genes across most tissues assessed, including many housekeeping genes, and positively correlates with H3K27ac and H3K4me3 as well as with gene expression. In addition, H3K18la is enriched at active enhancers that lie in proximity to genes that are functionally important for the respective tissue.

Conclusions

Overall, our data suggests that H3K18la is not only a marker for active promoters, but also a mark of tissue specific active enhancers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02775-y.

RNA Pol II pausing facilitates phased pluripotency transitions by buffering transcription

Promoter-proximal RNA Pol II pausing is a critical step in transcriptional control. Pol II pausing has been predominantly studied in tissue culture systems. While Pol II pausing has been shown to be required for mammalian development, the phenotypic and mechanistic details of this requirement are unknown. Here, we found that loss of Pol II pausing stalls pluripotent state transitions within the epiblast of the early mouse embryo. Using Nelfb -/- mice and a NELFB degron mouse pluripotent stem cell model, we show that embryonic stem cells (ESCs) representing the naïve state of pluripotency successfully initiate a transition program but fail to balance levels of induced and repressed genes and enhancers in the absence of NELF. We found an increase in chromatin-associated NELF during transition from the naïve to later pluripotent states. Overall, our work defines the acute and long-term molecular consequences of NELF loss and reveals a role for Pol II pausing in the pluripotency continuum as a modulator of cell state transitions.

Epigenetics as "conductor" in "orchestra" of pluripotent states

Pluripotent character is described as the potency of cells to differentiate into all three germ layers. The best example to reinstate the term lies in the context of embryonic stem cells (ESCs). Pluripotent ESC describes the in vitro status of those cells that originate during the complex process of embryogenesis. Pre-implantation to post-implantation development of embryo embrace cells with different levels of stemness. Currently, four states of pluripotency have been recognized, in the progressing order of "naïve," "poised," "formative," and "primed." Epigenetics act as the "conductor" in this "orchestra" of transition in pluripotent states. With a distinguishable gene expression profile, these four states associate with different epigenetic signatures, sometimes distinct while otherwise overlapping. The present review focuses on how epigenetic factors, including DNA methylation, bivalent chromatin, chromatin remodelers, chromatin/nuclear architecture, and microRNA, could dictate pluripotent states and their transition among themselves.

The human amniotic epithelium confers a bias to differentiate toward the neuroectoderm lineage in human embryonic stem cells

Human embryonic stem cells (hESCs) derive from the epiblast and have pluripotent potential. To maintain the conventional conditions of the pluripotent potential in an undifferentiated state, inactivated mouse embryonic fibroblast (iMEF) is used as a feeder layer. However, it has been suggested that hESC under this conventional condition (hESC-iMEF) is an artifact that does not correspond to the in vitro counterpart of the human epiblast. Our previous studies demonstrated the use of an alternative feeder layer of human amniotic epithelial cells (hAECs) to derive and maintain hESC. We wondered if the hESC-hAEC culture could represent a different pluripotent stage than that of naïve or primed conventional conditions, simulating the stage in which the amniotic epithelium derives from the epiblast during peri-implantation. Like the conventional primed hESC-iMEF, hESC-hAEC has the same levels of expression as the ‘pluripotency core’ and does not express markers of naïve pluripotency. However, it presents a downregulation of HOX genes and genes associated with the endoderm and mesoderm, and it exhibits an increase in the expression of ectoderm lineage genes, specifically in the anterior neuroectoderm. Transcriptome analysis showed in hESC-hAEC an upregulated signature of genes coding for transcription factors involved in neural induction and forebrain development, and the ability to differentiate into a neural lineage was superior in comparison with conventional hESC-iMEF. We propose that the interaction of hESC with hAEC confers hESC a biased potential that resembles the anteriorized epiblast, which is predisposed to form the neural ectoderm.

Trafficking regulator of GLUT4-1 (TRARG1) is a GSK3 substrate

Trafficking regulator of GLUT4-1, TRARG1, positively regulates insulin-stimulated GLUT4 trafficking and insulin sensitivity. However, the mechanism(s) by which this occurs remain(s) unclear. Using biochemical and mass spectrometry analyses we found that TRARG1 is dephosphorylated in response to insulin in a PI3K/Akt-dependent manner and is a novel substrate for GSK3. Priming phosphorylation of murine TRARG1 at serine 84 allows for GSK3-directed phosphorylation at serines 72, 76 and 80. A similar pattern of phosphorylation was observed in human TRARG1, suggesting that our findings are translatable to human TRARG1. Pharmacological inhibition of GSK3 increased cell surface GLUT4 in cells stimulated with a submaximal insulin dose, and this was impaired following Trarg1 knockdown, suggesting that TRARG1 acts as a GSK3-mediated regulator in GLUT4 trafficking. These data place TRARG1 within the insulin signaling network and provide insights into how GSK3 regulates GLUT4 trafficking in adipocytes.

Time-resolved Phosphoproteome and Proteome Analysis Reveals Kinase Signalling on Master Transcription Factors During Myogenesis

Myogenesis is governed by signaling networks that are tightly regulated in a time-dependent manner. Although different protein kinases have been identified, knowledge of the global signaling networks and their downstream substrates during myogenesis remains incomplete. Here, we map the myogenic differentiation of C2C12 cells using phosphoproteomics and proteomics. From these data, we infer global kinase activity and predict the substrates that are involved in myogenesis. We found that multiple mitogen-activated protein kinases (MAPKs) mark the initial wave of signaling cascades. Further phosphoproteomic and proteomic profiling with MAPK1/3 and MAPK8/9 specific inhibitions unveil their shared and distinctive roles in myogenesis. Lastly, we identified and validated the transcription factor nuclear factor 1 X-type (NFIX) as a novel MAPK1/3 substrate and demonstrated the functional impact of NFIX phosphorylation on myogenesis. Altogether, these data characterize the dynamics, interactions, and downstream control of kinase signaling networks during myogenesis on a global scale.

•

Phosphoproteomic and proteomic maps of myogenic differentiation of C2C12 cells

•

Myogenic kinome activity and kinase-substrates prediction using machine learning

•

MAPK1/3 and MAPK8/9 inhibition unveil shared and distinctive effects on myogenesis

•

Validation of NFIX phosphorylation by MAPK1/3 and its impact on myogenesis

Personalized phosphoproteomics identifies functional signaling

Protein phosphorylation dynamically integrates environmental and cellular information to control biological processes. Identifying functional phosphorylation amongst the thousands of phosphosites regulated by a perturbation at a global scale is a major challenge. Here we introduce 'personalized phosphoproteomics', a combination of experimental and computational analyses to link signaling with biological function by utilizing human phenotypic variance. We measure individual subject phosphoproteome responses to interventions with corresponding phenotypes measured in parallel. Applying this approach to investigate how exercise potentiates insulin signaling in human skeletal muscle, we identify both known and previously unidentified phosphosites on proteins involved in glucose metabolism. This includes a cooperative relationship between mTOR and AMPK whereby the former directly phosphorylates the latter on S377, for which we find a role in metabolic regulation. These results establish personalized phosphoproteomics as a general approach for investigating the signal transduction underlying complex biology.

High temporal resolution proteome and phosphoproteome profiling of stem cell-derived hepatocyte development

Primary human hepatocytes are widely used to evaluate liver toxicity of drugs, but they are scarce and demanding to culture. Stem cell-derived hepatocytes are increasingly discussed as alternatives. To obtain a better appreciation of the molecular processes during the differentiation of induced pluripotent stem cells into hepatocytes, we employ a quantitative proteomic approach to follow the expression of 9,000 proteins, 12,000 phosphorylation sites, and 800 acetylation sites over time. The analysis reveals stage-specific markers, a major molecular switch between hepatic endoderm versus immature hepatocyte-like cells impacting, e.g., metabolism, the cell cycle, kinase activity, and the expression of drug transporters. Comparing the proteomes of two- (2D) and three-dimensional (3D)-derived hepatocytes with fetal and adult liver indicates a fetal-like status of the in vitro models and lower expression of important ADME/Tox proteins. The collective data enable constructing a molecular roadmap of hepatocyte development that serves as a valuable resource for future research.

Functional analysis of the stable phosphoproteome reveals cancer vulnerabilities

Motivation

The advance of mass spectrometry-based technologies enabled the profiling of the phosphoproteomes of a multitude of cell and tissue types. However, current research primarily focused on investigating the phosphorylation dynamics in specific cell types and experimental conditions, whereas the phosphorylation events that are common across cell/tissue types and stable regardless of experimental conditions are, so far, mostly ignored.

Results

Here, we developed a statistical framework to identify the stable phosphoproteome across 53 human phosphoproteomics datasets, covering 40 cell/tissue types and 194 conditions/treatments. We demonstrate that the stably phosphorylated sites (SPSs) identified from our statistical framework are evolutionarily conserved, functionally important and enriched in a range of core signaling and gene pathways. Particularly, we show that SPSs are highly enriched in the RNA splicing pathway, an essential cellular process in mammalian cells, and frequently disrupted by cancer mutations, suggesting a link between the dysregulation of RNA splicing and cancer development through mutations on SPSs.

Availability and implementation

The source code for data analysis in this study is available from Github repository https://github.com/PYangLab/SPSs under the open-source license of GPL-3. The data used in this study are publicly available (see Section 2.8).

Supplementary information

Supplementary data are available at Bioinformatics online.

RNA polymerase II pausing in development: orchestrating transcription

The coordinated regulation of transcriptional networks underpins cellular identity and developmental progression. RNA polymerase II promoter-proximal pausing (Pol II pausing) is a prevalent mechanism by which cells can control and synchronize transcription. Pol II pausing regulates the productive elongation step of transcription at key genes downstream of a variety of signalling pathways, such as FGF and Nodal. Recent advances in our understanding of the Pol II pausing machinery and its role in transcription call for an assessment of these findings within the context of development. In this review, we discuss our current understanding of the molecular basis of Pol II pausing and its function during organismal development. By critically assessing the tools used to study this process we conclude that combining recently developed genomics approaches with refined perturbation systems has the potential to expand our understanding of Pol II pausing mechanistically and functionally in the context of development and beyond.

Measuring Endocytosis and Endosomal Uptake at Single Cell Resolution

Endocytosis impacts many cell biological functions, including in embryonic stem cells (ESCs). It has been shown that endocytosis is necessary for adequate FGF-signaling within the preimplantation ESC to post-implantation epiblast (EpiLC) pluripotency continuum and is required for proper levels of ERK activation. Quantitative methods at single cell resolution are needed to study endocytosis as well as its regulation and roles in these transitioning populations. The methods in this chapter provide an easily adaptable, multiplexable platform to monitor and quantify endosomal uptake at single cell resolution in live cells following receptor-mediated and non-receptor-mediated endocytosis, including nonspecific mechanisms such as pinocytosis.

Decoding the function of bivalent chromatin in development and cancer

Bivalent chromatin is characterized by the simultaneous presence of H3K4me3 and H3K27me3, histone modifications generally associated with transcriptionally active and repressed chromatin, respectively. Prevalent in embryonic stem cells (ESCs), bivalency is postulated to poise/prime lineage-controlling developmental genes for rapid activation during embryogenesis while maintaining a transcriptionally repressed state in the absence of activation cues; however, this hypothesis remains to be directly tested. Most gene promoters DNA hypermethylated in adult human cancers are bivalently marked in ESCs, and it was speculated that bivalency predisposes them for aberrant de novo DNA methylation and irreversible silencing in cancer, but evidence supporting this model is largely lacking. Here, we show that bivalent chromatin does not poise genes for rapid activation but protects promoters from de novo DNA methylation. Genome-wide studies in differentiating ESCs reveal that activation of bivalent genes is no more rapid than that of other transcriptionally silent genes, challenging the premise that H3K4me3 is instructive for transcription. H3K4me3 at bivalent promoters-a product of the underlying DNA sequence-persists in nearly all cell types irrespective of gene expression and confers protection from de novo DNA methylation. Bivalent genes in ESCs that are frequent targets of aberrant hypermethylation in cancer are particularly strongly associated with loss of H3K4me3/bivalency in cancer. Altogether, our findings suggest that bivalency protects reversibly repressed genes from irreversible silencing and that loss of H3K4me3 may make them more susceptible to aberrant DNA methylation in diseases such as cancer. Bivalency may thus represent a distinct regulatory mechanism for maintaining epigenetic plasticity.

Regulation of NANOG and SOX2 expression by activin A and a canonical WNT agonist in bovine embryonic stem cells and blastocysts

Bovine embryonic stem cells (ESC) have features associated with the primed pluripotent state including low expression of one of the core pluripotency transcription factors, NANOG. It has been reported that NANOG expression can be upregulated in porcine ESC by treatment with activin A and the WNT agonist CHIR99021. Accordingly, it was tested whether expression of NANOG and another pluripotency factor SOX2 could be stimulated by activin A and the WNT agonist CHIR99021. Immunoreactive NANOG and SOX2 were analyzed for bovine ESC lines derived under conditions in which activin A and CHIR99021 were added singly or in combination. Activin A enhanced NANOG expression but also reduced SOX2 expression. CHIR99021 depressed expression of both NANOG and SOX2. In a second experiment, activin A enhanced blastocyst development while CHIR99021 treatment impaired blastocyst formation and reduced number of blastomeres. Activin A treatment decreased blastomeres in the blastocyst that were positive for either NANOG or SOX2 but increased those that were CDX2+ and that were GATA6+ outside the inner cell mass. CHIR99021 reduced SOX2+ and NANOG+ blastomeres without affecting the number or percent of blastomeres that were CDX2+ and GATA6+. Results indicate activation of activin A signaling stimulates NANOG expression during self-renewal of bovine ESC but suppresses cells expressing pluripotency markers in the blastocyst and increases cells expressing CDX2. Actions of activin A to promote blastocyst development may involve its role in promoting trophectoderm formation. Furthermore, results demonstrate the negative role of canonical WNT signaling in cattle for pluripotency marker expression in ESC and in formation of the inner cell mass and epiblast during embryonic development. This article has an associated First Person interview with the first author of the paper.

The ETS transcription factor ERF controls the exit from the naïve pluripotent state in a MAPK-dependent manner

The naïve epiblast transitions to a pluripotent primed state during embryo implantation. Despite the relevance of the FGF pathway during this period, little is known about the downstream effectors regulating this signaling. Here, we examined the molecular mechanisms coordinating the naïve to primed transition by using inducible ESC to genetically eliminate all RAS proteins. We show that differentiated RAS^KO ESC remain trapped in an intermediate state of pluripotency with naïve-associated features. Elimination of the transcription factor ERF overcomes the developmental blockage of RAS-deficient cells by naïve enhancer decommissioning. Mechanistically, ERF regulates NANOG expression and ensures naïve pluripotency by strengthening naïve transcription factor binding at ESC enhancers. Moreover, ERF negatively regulates the expression of the methyltransferase DNMT3B, which participates in the extinction of the naïve transcriptional program. Collectively, we demonstrated an essential role for ERF controlling the exit from naïve pluripotency in a MAPK-dependent manner during the progression to primed pluripotency.

Enhancer-associated H3K4 methylation safeguards in vitro germline competence

Germline specification in mammals occurs through an inductive process whereby competent cells in the post-implantation epiblast differentiate into primordial germ cells (PGC). The intrinsic factors that endow epiblast cells with the competence to respond to germline inductive signals remain unknown. Single-cell RNA sequencing across multiple stages of an in vitro PGC-like cells (PGCLC) differentiation system shows that PGCLC genes initially expressed in the naïve pluripotent stage become homogeneously dismantled in germline competent epiblast like-cells (EpiLC). In contrast, the decommissioning of enhancers associated with these germline genes is incomplete. Namely, a subset of these enhancers partly retain H3K4me1, accumulate less heterochromatic marks and remain accessible and responsive to transcriptional activators. Subsequently, as in vitro germline competence is lost, these enhancers get further decommissioned and lose their responsiveness to transcriptional activators. Importantly, using H3K4me1-deficient cells, we show that the loss of this histone modification reduces the germline competence of EpiLC and decreases PGCLC differentiation efficiency. Our work suggests that, although H3K4me1 might not be essential for enhancer function, it can facilitate the (re)activation of enhancers and the establishment of gene expression programs during specific developmental transitions.

An extensive and dynamic trans-omic network illustrating prominent regulatory mechanisms in response to insulin in the liver

An effective combination of multi-omic datasets can enhance our understanding of complex biological phenomena. To build a context-dependent network with multiple omic layers, i.e., a trans-omic network, we perform phosphoproteomics, transcriptomics, proteomics, and metabolomics of murine liver for 4 h after insulin administration and integrate the resulting time series. Structural characteristics and dynamic nature of the network are analyzed to elucidate the impact of insulin. Early and prominent changes in protein phosphorylation and persistent and asynchronous changes in mRNA and protein levels through non-transcriptional mechanisms indicate enhanced crosstalk between phosphorylation-mediated signaling and protein expression regulation. Metabolic response shows different temporal regulation with transient increases at early time points across categories and enhanced response in the amino acid and nucleotide categories at later time points as a result of process convergence. This extensive and dynamic view of the trans-omic network elucidates prominent regulatory mechanisms that drive insulin responses through intricate interlayer coordination.

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.

Phosphorylation of the HP1β hinge region sequesters KAP1 in heterochromatin and promotes the exit from naïve pluripotency

Heterochromatin binding protein HP1β plays an important role in chromatin organization and cell differentiation, however the underlying mechanisms remain unclear. Here, we generated HP1β^−/− embryonic stem cells and observed reduced heterochromatin clustering and impaired differentiation. We found that during stem cell differentiation, HP1β is phosphorylated at serine 89 by CK2, which creates a binding site for the pluripotency regulator KAP1. This phosphorylation dependent sequestration of KAP1 in heterochromatin compartments causes a downregulation of pluripotency factors and triggers pluripotency exit. Accordingly, HP1β^−/− and phospho-mutant cells exhibited impaired differentiation, while ubiquitination-deficient KAP1^−/− cells had the opposite phenotype with enhanced differentiation. These results suggest that KAP1 regulates pluripotency via its ubiquitination activity. We propose that the formation of subnuclear membraneless heterochromatin compartments may serve as a dynamic reservoir to trap or release cellular factors. The sequestration of essential regulators defines a novel and active role of heterochromatin in gene regulation and represents a dynamic mode of remote control to regulate cellular processes like cell fate decisions.

INO80 promotes H2A.Z occupancy to regulate cell fate transition in pluripotent stem cells

The INO80 chromatin remodeler is involved in many chromatin-dependent cellular functions. However, its role in pluripotency and cell fate transition is not fully defined. We examined the impact of Ino80 deletion in the naïve and primed pluripotent stem cells. We found that Ino80 deletion had minimal effect on self-renewal and gene expression in the naïve state, but led to cellular differentiation and de-repression of developmental genes in the transition toward and maintenance of the primed state. In the naïve state, INO80 pre-marked gene promoters that would adopt bivalent histone modifications by H3K4me3 and H3K27me3 upon transition into the primed state. In the primed state, in contrast to its known role in H2A.Z exchange, INO80 promoted H2A.Z occupancy at these bivalent promoters and facilitated H3K27me3 installation and maintenance as well as downstream gene repression. Together, our results identified an unexpected function of INO80 in H2A.Z deposition and gene regulation. We showed that INO80-dependent H2A.Z occupancy is a critical licensing step for the bivalent domains, and thereby uncovered an epigenetic mechanism by which chromatin remodeling, histone variant deposition and histone modification coordinately control cell fate.

Three-dimensional geometry controls division symmetry in stem cell colonies

Proper control of division orientation and symmetry, largely determined by spindle positioning, is essential to development and homeostasis. Spindle positioning has been extensively studied in cells dividing in two-dimensional (2D) environments and in epithelial tissues, where proteins such as NuMA (also known as NUMA1) orient division along the interphase long axis of the cell. However, little is known about how cells control spindle positioning in three-dimensional (3D) environments, such as early mammalian embryos and a variety of adult tissues. Here, we use mouse embryonic stem cells (ESCs), which grow in 3D colonies, as a model to investigate division in 3D. We observe that, at the periphery of 3D colonies, ESCs display high spindle mobility and divide asymmetrically. Our data suggest that enhanced spindle movements are due to unequal distribution of the cell-cell junction protein E-cadherin between future daughter cells. Interestingly, when cells progress towards differentiation, division becomes more symmetric, with more elongated shapes in metaphase and enhanced cortical NuMA recruitment in anaphase. Altogether, this study suggests that in 3D contexts, the geometry of the cell and its contacts with neighbors control division orientation and symmetry. This article has an associated First Person interview with the first author of the paper.

A proteome signature for acute incisional pain in dorsal root ganglia of mice

ABSTRACT

After surgery, acute pain is still managed insufficiently and may lead to short-term and long-term complications including chronic postsurgical pain and an increased prescription of opioids. Thus, identifying new targets specifically implicated in postoperative pain is of utmost importance to develop effective and nonaddictive analgesics. Here, we used an integrated and multimethod workflow to reveal unprecedented insights into proteome dynamics in dorsal root ganglia (DRG) of mice after plantar incision (INC). Based on a detailed characterization of INC-associated pain-related behavior profiles, including a novel paradigm for nonevoked pain, we performed quantitative mass-spectrometry-based proteomics in DRG 1 day after INC. Our data revealed a hitherto unknown INC-regulated protein signature in DRG with changes in distinct proteins and cellular signaling pathways. In particular, we show the differential regulation of 44 protein candidates, many of which are annotated with pathways related to immune and inflammatory responses such as MAPK/extracellular signal-regulated kinases signaling. Subsequent orthogonal assays comprised multiplex Western blotting, bioinformatic protein network analysis, and immunolabeling in independent mouse cohorts to validate (1) the INC-induced regulation of immune/inflammatory pathways and (2) the high priority candidate Annexin A1. Taken together, our results propose novel potential targets in the context of incision and, therefore, represent a highly valuable resource for further mechanistic and translational studies of postoperative pain.

Unraveling the Spatiotemporal Human Pluripotency in Embryonic Development

There have been significant advances in understanding human embryogenesis using human pluripotent stem cells (hPSCs) in conventional monolayer and 3D self-organized cultures. Thus, in vitro models have contributed to elucidate the molecular mechanisms for specification and differentiation during development. However, the molecular and functional spectrum of human pluripotency (i.e., intermediate states, pluripotency subtypes and regionalization) is still not fully understood. This review describes the mechanisms that establish and maintain pluripotency in human embryos and their differences with mouse embryos. Further, it describes a new pluripotent state representing a transition between naïve and primed pluripotency. This review also presents the data that divide pluripotency into substates expressing epiblast regionalization and amnion specification as well as primordial germ cells in primates. Finally, this work analyzes the amnion's relevance as an "signaling center" for regionalization before the onset of gastrulation.