Resolving Cell Fate Decisions during Somatic Cell Reprogramming by Single-Cell RNA-Seq.

Somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs), which is a highly heterogeneous process. Here we report the cell fate continuum during somatic cell reprogramming at single-cell resolution. We first develop SOT to analyze cell fate continuum from Oct4/Sox2/Klf4- or OSK-mediated reprogramming and show that cells bifurcate into two categories, reprogramming potential (RP) or non-reprogramming (NR). We further show that Klf4 contributes to Cd34+/Fxyd5+/Psca+ keratinocyte-like NR fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the RP cells. We analyze more than 150,000 single cells from both OSK and chemical reprograming and identify additional NR/RP bifurcation points. Our work reveals a generic bifurcation model for cell fate decisions during somatic cell reprogramming that may be applicable to other systems and inspire further improvements for reprogramming.

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

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

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements.

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

A programmable system to methylate and demethylate N6-methyladenosine (m6A) on specific RNA transcripts in mammalian cells.

RNA N6-methyladenosine (m6A) is the most abundant internal mRNA modification and forms part of an epitranscriptomic system that modulates RNA function. m6A is reversibly catalyzed by specific enzymes, and those modifications can be recognized by RNA-binding proteins that in turn regulate biological processes. Although there are many reports demonstrating m6A participation in critical biological functions, this exploration has mainly been conducted through the global KO or knockdown of the writers, erasers, or readers of m6A. Consequently, there is a lack of information about the role of m6A on single transcripts in biological processes, posing a challenge in understanding the biological functions of m6A. Here, we demonstrate a CRISPR/dCas13a-based RNA m6A editors, which can target RNAs using a single or multiple CRISPR RNA array to methylate or demethylate m6A in human 293T cells and mouse embryonic stem cells. We systematically assay its capabilities to enable the targeted rewriting of m6A dynamics, including modulation of circular RNA translation and transcript half-life. Finally, we use the system to specifically modulate m6A levels on the noncoding XIST (X-inactive specific transcript) to modulate X chromosome silencing and activation. The editors described here can be used to explore the roles of m6A in biological processes.

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

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

Systematic screening of CTCF binding partners identifies that BHLHE40 regulates CTCF genome-wide distribution and long-range chromatin interactions.

CTCF plays a pivotal role in mediating chromatin interactions, but it does not do so alone. A number of factors have been reported to co-localize with CTCF and regulate CTCF loops, but no comprehensive analysis of binding partners has been performed. This prompted us to identify CTCF loop participants and regulators by co-localization analysis with CTCF. We screened all factors that had ChIP-seq data in humans by co-localization analysis with human super conserved CTCF (hscCTCF) binding sites, and identified many new factors that overlapped with hscCTCF binding sites. Combined with CTCF loop information, we observed that clustered factors could promote CTCF loops. After in-depth mining of each factor, we found that many factors might have the potential to promote CTCF loops. Our data further demonstrated that BHLHE40 affected CTCF loops by regulating CTCF binding. Together, this study revealed that many factors have the potential to participate in or regulate CTCF loops, and discovered a new role for BHLHE40 in modulating CTCF loop formation.

Capturing the interactome of newly transcribed RNA.

We combine the labeling of newly transcribed RNAs with 5-ethynyluridine with the characterization of bound proteins. This approach, named capture of the newly transcribed RNA interactome using click chemistry (RICK), systematically captures proteins bound to a wide range of RNAs, including nascent RNAs and traditionally neglected nonpolyadenylated RNAs. RICK has identified mitotic regulators amongst other novel RNA-binding proteins with preferential affinity for nonpolyadenylated RNAs, revealed a link between metabolic enzymes/factors and nascent RNAs, and expanded the known RNA-bound proteome of mouse embryonic stem cells. RICK will facilitate an in-depth interrogation of the total RNA-bound proteome in different cells and systems.

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

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

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

The current classification of cells in an organism is largely based on their anatomic and developmental origin. Cells types and tissues are traditionally classified into those that arise from the three embryonic germ layers, the ectoderm, mesoderm and endoderm, but this model does not take into account the organization of cell type-specific patterns of gene expression. Here, we present computational models for cell type and tissue specification derived from a collection of 921 RNA-sequencing samples from 272 distinct mouse cell types or tissues. In an unbiased fashion, this analysis accurately predicts the three known germ layers. Unexpectedly, this analysis also suggests that in total there are eight major domains of cell type-specification, corresponding to the neurectoderm, neural crest, surface ectoderm, endoderm, mesoderm, blood mesoderm, germ cells and the embryonic domain. Further, we identify putative genes responsible for specifying the domain and the cell type. This model has implications for understanding trans-lineage differentiation for stem cells, developmental cell biology and regenerative medicine.

Oct4 switches partnering from Sox2 to Sox17 to reinterpret the enhancer code and specify endoderm.

How regulatory information is encoded in the genome is poorly understood and poses a challenge when studying biological processes. We demonstrate here that genomic redistribution of Oct4 by alternative partnering with Sox2 and Sox17 is a fundamental regulatory event of endodermal specification. We show that Sox17 partners with Oct4 and binds to a unique 'compressed' Sox/Oct motif that earmarks endodermal genes. This is in contrast to the pluripotent state where Oct4 selectively partners with Sox2 at 'canonical' binding sites. The distinct selection of binding sites by alternative Sox/Oct partnering is underscored by our demonstration that rationally point-mutated Sox17 partners with Oct4 on pluripotency genes earmarked by the canonical Sox/Oct motif. In an endodermal differentiation assay, we demonstrate that the compressed motif is required for proper expression of endodermal genes. Evidently, Oct4 drives alternative developmental programs by switching Sox partners that affects enhancer selection, leading to either an endodermal or pluripotent cell fate. This work provides insights in understanding cell fate transcriptional regulation by highlighting the direct link between the DNA sequence of an enhancer and a developmental outcome.

PTP-central: a comprehensive resource of protein tyrosine phosphatases in eukaryotic genomes.

Reversible tyrosine phosphorylation is a fundamental signaling mechanism controlling a diversity of cellular processes. Whereas protein tyrosine kinases have long been implicated in many diseases, aberrant protein tyrosine phosphatase (PTP) activity is also increasingly being associated with a wide spectrum of conditions. PTPs are now regarded as key regulators of biochemical processes instead of simple "off" switches operating in tyrosine kinase signaling pathways. Despite the central importance that PTPs play in the cell's biochemistry, the tyrosine phosphatomes of most species remain uncharted. Here we present a highly sensitive and specific sequence-based method for the automatic classification of PTPs. As proof of principle we re-annotated the human tyrosine phosphatome, and discovered four new PTP genes that had not been reported before. Our method and the predicted tyrosine phosphatomes of 65 eukaryotic genomes are accessible online through the user-friendly PTP-central resource (http://www.PTP-central.org/), where users can also submit their own sequences for prediction. PTP-central is a comprehensive and continually developing resource that currently integrates the predicted tyrosine phosphatomes with structural data and genetic association disease studies, as well as homology relationships. PTP-central thus fills an important void for the systematic study of PTPs, both in model organisms and from an evolutionary perspective.

Characterization of PTPN2 and its use as a biomarker.

For years, the two main isoforms of PTPN2 have been an interesting yet academic topic of debate for researchers working on this phosphatase. In recent years, several studies were published in which these isoforms were attributed specific functions. Most importantly, differences in their stoichiometry have been reported to be associated with certain diseases such as inflammatory bowel diseases (IBDs). Hence, understanding the evolutionary ontogeny of the main transcripts and the physiological consequences of their expression have now become clinically relevant issues. Herein we describe the genomic controls placed upon PTPN2, the identified splice variants, the encoded PTPN2 proteins, and both the known and putative post-translational modifications that have been reported. Moreover, we examine the expression of PTPN2 isoforms in specific tissues as well as in a disease setting. PTPN2 is an important negative regulator of inflammation. Therefore, the following protocols are effective approaches for its adequate monitoring in inflammatory diseases' progression and outcome.

VGLL1 cooperates with TEAD4 to control human trophectoderm lineage specification.

In contrast to rodents, the mechanisms underlying human trophectoderm and early placenta specification are understudied due to ethical barriers and the scarcity of embryos. Recent reports have shown that human pluripotent stem cells (PSCs) can differentiate into trophectoderm (TE)-like cells (TELCs) and trophoblast stem cells (TSCs), offering a valuable in vitro model to study early placenta specification. Here, we demonstrate that the VGLL1 (vestigial-like family member 1), which is highly expressed during human and non-human primate TE specification in vivo but is negligibly expressed in mouse, is a critical regulator of cell fate determination and self-renewal in human TELCs and TSCs derived from naïve PSCs. Mechanistically, VGLL1 partners with the transcription factor TEAD4 (TEA domain transcription factor 4) to regulate chromatin accessibility at target gene loci through histone acetylation and acts in cooperation with GATA3 and TFAP2C. Our work is relevant to understand primate early embryogenesis and how it differs from other mammalian species.

Single-cell RNA sequencing reveals immune cell dysfunction in the peripheral blood of patients with highly aggressive gastric cancer.

Highly aggressive gastric cancer (HAGC) is a gastric cancer characterized by bone marrow metastasis and disseminated intravascular coagulation (DIC). Information about the disease is limited. Here we employed single-cell RNA sequencing to investigate peripheral blood mononuclear cells (PBMCs), aiming to unravel the immune response of patients toward HAGC. PBMCs from seven HAGC patients, six normal advanced gastric cancer (NAGC) patients, and five healthy individuals were analysed by single-cell RNA sequencing. The expression of genes of interest was validated by bulk RNA-sequencing and ELISA. We found a massive expansion of neutrophils in PBMCs of HAGC. These neutrophils are activated, but immature. Besides, mononuclear phagocytes exhibited an M2-like signature and T cells were suppressed and reduced in number. Analysis of cell-cell crosstalk revealed that several signalling pathways involved in neutrophil to T-cell suppression including APP-CD74, MIF-(CD74+CXCR2), and MIF-(CD74+CD44) pathways were increased in HAGC. NETosis-associated genes S100A8 and S100A9 as well as VEGF, PDGF, FGF, and NOTCH signalling that contribute to DIC development were upregulated in HAGC too. This study reveals significant changes in the distribution and interactions of the PBMC subsets and provides valuable insight into the immune response in patients with HAGC. S100A8 and S100A9 are highly expressed in HAGC neutrophils, suggesting their potential to be used as novel diagnostic and therapeutic targets for HAGC.

A spatiotemporal atlas of mouse liver homeostasis and regeneration.

The mechanism by which mammalian liver cell responses are coordinated during tissue homeostasis and perturbation is poorly understood, representing a major obstacle in our understanding of many diseases. This knowledge gap is caused by the difficulty involved with studying multiple cell types in different states and locations, particularly when these are transient. We have combined Stereo-seq (spatiotemporal enhanced resolution omics-sequencing) with single-cell transcriptomic profiling of 473,290 cells to generate a high-definition spatiotemporal atlas of mouse liver homeostasis and regeneration at the whole-lobe scale. Our integrative study dissects in detail the molecular gradients controlling liver cell function, systematically defining how gene networks are dynamically modulated through intercellular communication to promote regeneration. Among other important regulators, we identified the transcriptional cofactor TBL1XR1 as a rheostat linking inflammation to Wnt/ß-catenin signaling for facilitating hepatocyte proliferation. Our data and analytical pipelines lay the foundation for future high-definition tissue-scale atlases of organ physiology and malfunction.

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

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

c-JUN is a barrier in hESC to cardiomyocyte transition.

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

c-Jun as a one-way valve at the naive to primed interface.

BACKGROUND: c-Jun is a proto-oncogene functioning as a transcription factor to activate gene expression under many physiological and pathological conditions, particularly in somatic cells. However, its role in early embryonic development remains unknown. RESULTS: Here, we show that c-Jun acts as a one-way valve to preserve the primed state and impair reversion to the naïve state. c-Jun is induced during the naive to primed transition, and it works to stabilize the chromatin structure and inhibit the reverse transition. Loss of c-Jun has surprisingly little effect on the naïve to primed transition, and no phenotypic effect on primed cells, however, in primed cells the loss of c-Jun leads to a failure to correctly close naïve-specific enhancers. When the primed cells are induced to reprogram to a naïve state, these enhancers are more rapidly activated when c-Jun is lost or impaired, and the conversion is more efficient. CONCLUSIONS: The results of this study indicate that c-Jun can function as a chromatin stabilizer in primed EpiSCs, to maintain the epigenetic cell type state and act as a one-way valve for cell fate conversions.

Disruption of neuronal RHEB signaling impairs oligodendrocyte differentiation and myelination through mTORC1-DLK1 axis.

How neuronal signaling affects brain myelination remains poorly understood. We show dysregulated neuronal RHEB-mTORC1-DLK1 axis impairs brain myelination. Neuronal Rheb cKO impairs oligodendrocyte differentiation/myelination, with activated neuronal expression of the imprinted gene Dlk1. Neuronal Dlk1 cKO ameliorates myelination deficit in neuronal Rheb cKO mice, indicating that activated neuronal Dlk1 expression contributes to impaired myelination caused by Rheb cKO. The effect of Rheb cKO on Dlk1 expression is mediated by mTORC1; neuronal mTor cKO and Raptor cKO and pharmacological inhibition of mTORC1 recapitulate elevated neuronal Dlk1 expression. We demonstrate that both a secreted form of DLK1 and a membrane-bound DLK1 inhibit the differentiation of cultured oligodendrocyte precursor cells into oligodendrocytes expressing myelin proteins. Finally, neuronal expression of Dlk1 in transgenic mice reduces the formation of mature oligodendrocytes and myelination. This study identifies Dlk1 as an inhibitor of oligodendrocyte myelination and a mechanism linking altered neuronal signaling with oligodendrocyte dysfunction.

In vitro generation of mouse morula-like cells.

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