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

Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression

Although more than half of all genes generate transcripts that differ in 3'UTR length, current analysis pipelines only quantify the amount but not the length of mRNA transcripts. 3'UTR length is determined by 3' end cleavage sites (CS). We map CS in more than 200 primary human and mouse cell types and increase CS annotations relative to the GENCODE database by 40%. Approximately half of all CS are used in few cell types, revealing that most genes only have one or two major 3' ends. We incorporate the CS annotations into a computational pipeline, called scUTRquant, for rapid, accurate, and simultaneous quantification of gene and 3'UTR isoform expression from single-cell RNA sequencing (scRNA-seq) data. When applying scUTRquant to data from 474 cell types and 2134 perturbations, we discover extensive 3'UTR length changes across cell types that are as widespread and coordinately regulated as gene expression changes but affect mostly different genes. Our data indicate that mRNA abundance and mRNA length are two largely independent axes of gene regulation that together determine the amount and spatial organization of protein synthesis.

Cell-Reprogramming-Inspired Dynamically Responsive Hydrogel Boosts the Induction of Pluripotency via Phase-Separated Biomolecular Condensates

Induced pluripotent stem cells (iPSCs) have wide applications in disease modeling, personalized medicine, and tissue engineering. The generation of iPSCs from somatic cells via transcriptional-factor- or chemical molecule-based approaches are time-consuming and inefficient. Here, a cell-reprogramming-inspired dynamically responsive hydrogel is fabricated via a synthetic-biology-based strategy. Human and mouse somatic cells (including senescent cells) are efficiently reprogrammed into iPSCs that exhibit key features of embryonic stem cells. The cell-reprogramming-responsive hydrogel possesses dynamic bioresponsiveness, and it faithfully senses metabolic remodeling and extracellular acidification during cell reprogramming, responding by changing its mechanical properties accordingly. Mechanistic study demonstrates that the autonomous change of the mechanical properties of the cell-reprogramming-responsive hydrogel elicits the formation of Yes-associated protein (YAP) biomolecular condensates with the appropriate timing during cell reprogramming, ensuring a faster and more efficient generation of iPSCs than conventional cell reprogramming approach. Taken together, this study reveals the robust induction of pluripotency by coordination of cell-reprogramming-inspired dynamically responsive hydrogel and phase-separated biomolecular condensates.

Induced pluripotent stem cells (iPSCs): molecular mechanisms of induction and applications

The induced pluripotent stem cell (iPSC) technology has transformed in vitro research and holds great promise to advance regenerative medicine. iPSCs have the capacity for an almost unlimited expansion, are amenable to genetic engineering, and can be differentiated into most somatic cell types. iPSCs have been widely applied to model human development and diseases, perform drug screening, and develop cell therapies. In this review, we outline key developments in the iPSC field and highlight the immense versatility of the iPSC technology for in vitro modeling and therapeutic applications. We begin by discussing the pivotal discoveries that revealed the potential of a somatic cell nucleus for reprogramming and led to successful generation of iPSCs. We consider the molecular mechanisms and dynamics of somatic cell reprogramming as well as the numerous methods available to induce pluripotency. Subsequently, we discuss various iPSC-based cellular models, from mono-cultures of a single cell type to complex three-dimensional organoids, and how these models can be applied to elucidate the mechanisms of human development and diseases. We use examples of neurological disorders, coronavirus disease 2019 (COVID-19), and cancer to highlight the diversity of disease-specific phenotypes that can be modeled using iPSC-derived cells. We also consider how iPSC-derived cellular models can be used in high-throughput drug screening and drug toxicity studies. Finally, we discuss the process of developing autologous and allogeneic iPSC-based cell therapies and their potential to alleviate human diseases.

Cellular plasticity in reprogramming, rejuvenation and tumorigenesis: a pioneer TF perspective

The multistep process of in vivo reprogramming, mediated by the transcription factors (TFs) Oct4, Sox2, Klf4, and c-Myc (OSKM), holds great promise for the development of rejuvenating and regenerative strategies. However, most of the approaches developed so far are accompanied by a persistent risk of tumorigenicity. Here, we review the groundbreaking effects of in vivo reprogramming with a particular focus on rejuvenation and regeneration. We discuss how the activity of pioneer TFs generates cellular plasticity that may be critical for inducing not only reprogramming and regeneration, but also cancer initiation. Finally, we highlight how a better understanding of the uncoupled control of cellular identity, plasticity, and aging during reprogramming might pave the way to the development of rejuvenating/regenerating strategies in a nontumorigenic manner.

Retrospective identification of cell-intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of OCT4, KLF4, SOX2, and MYC. Usually only a rare subset of cells reprogram, and the molecular characteristics of this subset remain unknown. We apply retrospective clone tracing to identify and characterize the rare human fibroblasts primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis increased the reprogramming efficiency. We provide evidence for a unified model in which cells can move into and out of the primed state over time, explaining how reprogramming appears deterministic at short timescales and stochastic at long timescales. Furthermore, inhibiting the activity of LSD1 enlarged the pool of cells that were primed for reprogramming. Thus, even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Metabolic control of induced pluripotency

Pluripotent stem cells of the mammalian epiblast and their cultured counterparts-embryonic stem cells (ESCs) and epiblast stem cells (EpiSCs)-have the capacity to differentiate in all cell types of adult organisms. An artificial process of reactivation of the pluripotency program in terminally differentiated cells was established in 2006, which allowed for the generation of induced pluripotent stem cells (iPSCs). This iPSC technology has become an invaluable tool in investigating the molecular mechanisms of human diseases and therapeutic drug development, and it also holds tremendous promise for iPSC applications in regenerative medicine. Since the process of induced reprogramming of differentiated cells to a pluripotent state was discovered, many questions about the molecular mechanisms involved in this process have been clarified. Studies conducted over the past 2 decades have established that metabolic pathways and retrograde mitochondrial signals are involved in the regulation of various aspects of stem cell biology, including differentiation, pluripotency acquisition, and maintenance. During the reprogramming process, cells undergo major transformations, progressing through three distinct stages that are regulated by different signaling pathways, transcription factor networks, and inputs from metabolic pathways. Among the main metabolic features of this process, representing a switch from the dominance of oxidative phosphorylation to aerobic glycolysis and anabolic processes, are many critical stage-specific metabolic signals that control the path of differentiated cells toward a pluripotent state. In this review, we discuss the achievements in the current understanding of the molecular mechanisms of processes controlled by metabolic pathways, and vice versa, during the reprogramming process.

Deciphering the decisive factors driving fate bifurcations in somatic cell reprogramming

Single-cell studies have demonstrated that somatic cell reprogramming is a continuous process of cell fates transition. Only partial reprogramming intermediates can overcome the molecular bottlenecks to acquire pluripotency. To decipher the underlying decisive factors driving cell fate, we identified induced pluripotent stem cells or stromal-like cells (iPSCs/SLCs) and iPSCs or trophoblast-like cells (iPSCs/TLCs) fate bifurcations by reconstructing cellular trajectory. The mesenchymal-epithelial transition and the activation of pluripotency networks are the main molecular series in successful reprogramming. Correspondingly, intermediates diverge into SLCs accompanied by the inhibition of cell cycle genes and the activation of extracellular matrix genes, whereas the TLCs fate is characterized by the up-regulation of placenta development genes. Combining putative gene regulatory networks, seven (Taf7, Ezh2, Klf2, etc.) and three key factors (Cdc5l, Klf4, and Nanog) were individually identified as drivers of the successful reprogramming by triggering downstream pluripotent networks during iPSCs/SLCs and iPSCs/TLCs fate bifurcation. Conversely, 11 factors (Cebpb, Sox4, Junb, etc.) and four factors (Gata2, Jund, Ctnnb1, etc.) drive SLCs fate and TLCs fate, respectively. Our study sheds new light on the understanding of decisive factors driving cell fate, which is helpful for improving reprogramming efficiency through manipulating cell fates to avoid alternative fates.

Complex haploinsufficiency in pluripotent cells yields somatic cells with DNA methylation abnormalities and pluripotency induction defects

A complete knockout of a single key pluripotency gene may drastically affect embryonic stem cell function and epigenetic reprogramming. In contrast, elimination of only one allele of a single pluripotency gene is mostly considered harmless to the cell. To understand whether complex haploinsufficiency exists in pluripotent cells, we simultaneously eliminated a single allele in different combinations of two pluripotency genes (i.e., Nanog+/-;Sall4+/-, Nanog+/-;Utf1+/-, Nanog+/-;Esrrb+/- and Sox2+/-;Sall4+/-). Although these double heterozygous mutant lines similarly contribute to chimeras, fibroblasts derived from these systems show a significant decrease in their ability to induce pluripotency. Tracing the stochastic expression of Sall4 and Nanog at early phases of reprogramming could not explain the seen delay or blockage. Further exploration identifies abnormal methylation around pluripotent and developmental genes in the double heterozygous mutant fibroblasts, which could be rescued by hypomethylating agent or high OSKM levels. This study emphasizes the importance of maintaining two intact alleles for pluripotency induction.

Cyclic Stretch Promotes Cellular Reprogramming Process through Cytoskeletal-Nuclear Mechano-Coupling and Epigenetic Modification

Advancing the technologies for cellular reprogramming with high efficiency has significant impact on regenerative therapy, disease modeling, and drug discovery. Biophysical cues can tune the cell fate, yet the precise role of external physical forces during reprogramming remains elusive. Here the authors show that temporal cyclic-stretching of fibroblasts significantly enhances the efficiency of induced pluripotent stem cell (iPSC) production. Generated iPSCs are proven to express pluripotency markers and exhibit in vivo functionality. Bulk RNA-sequencing reveales that cyclic-stretching enhances biological characteristics required for pluripotency acquisition, including increased cell division and mesenchymal-epithelial transition. Of note, cyclic-stretching activates key mechanosensitive molecules (integrins, perinuclear actins, nesprin-2, and YAP), across the cytoskeletal-to-nuclear space. Furthermore, stretch-mediated cytoskeletal-nuclear mechano-coupling leads to altered epigenetic modifications, mainly downregulation in H3K9 methylation, and its global gene occupancy change, as revealed by genome-wide ChIP-sequencing and pharmacological inhibition tests. Single cell RNA-sequencing further identifies subcluster of mechano-responsive iPSCs and key epigenetic modifier in stretched cells. Collectively, cyclic-stretching activates iPSC reprogramming through mechanotransduction process and epigenetic changes accompanied by altered occupancy of mechanosensitive genes. This study highlights the strong link between external physical forces with subsequent mechanotransduction process and the epigenetic changes with expression of related genes in cellular reprogramming, holding substantial implications in the field of cell biology, tissue engineering, and regenerative medicine.

Transcription factor stoichiometry, motif affinity and syntax regulate single-cell chromatin dynamics during fibroblast reprogramming to pluripotency

Ectopic expression of OCT4, SOX2, KLF4 and MYC (OSKM) transforms differentiated cells into induced pluripotent stem cells. To refine our mechanistic understanding of reprogramming, especially during the early stages, we profiled chromatin accessibility and gene expression at single-cell resolution across a densely sampled time course of human fibroblast reprogramming. Using neural networks that map DNA sequence to ATAC-seq profiles at base-resolution, we annotated cell-state-specific predictive transcription factor (TF) motif syntax in regulatory elements, inferred affinity- and concentration-dependent dynamics of Tn5-bias corrected TF footprints, linked peaks to putative target genes, and elucidated rewiring of TF-to-gene cis-regulatory networks. Our models reveal that early in reprogramming, OSK, at supraphysiological concentrations, rapidly open transient regulatory elements by occupying non-canonical low-affinity binding sites. As OSK concentration falls, the accessibility of these transient elements decays as a function of motif affinity. We find that these OSK-dependent transient elements sequester the somatic TF AP-1. This redistribution is strongly associated with the silencing of fibroblast-specific genes within individual nuclei. Together, our integrated single-cell resource and models reveal insights into the cis-regulatory code of reprogramming at unprecedented resolution, connect TF stoichiometry and motif syntax to diversification of cell fate trajectories, and provide new perspectives on the dynamics and role of transient regulatory elements in somatic silencing.

PC3T: a signature-driven predictor of chemical compounds for cellular transition

Cellular transitions hold great promise in translational medicine research. However, therapeutic applications are limited by the low efficiency and safety concerns of using transcription factors. Small molecules provide a temporal and highly tunable approach to overcome these issues. Here, we present PC3T, a computational framework to enrich molecules that induce desired cellular transitions, and PC3T was able to consistently enrich small molecules that had been experimentally validated in both bulk and single-cell datasets. We then predicted small molecule reprogramming of fibroblasts into hepatic progenitor-like cells (HPLCs). The converted cells exhibited epithelial cell-like morphology and HPLC-like gene expression pattern. Hepatic functions were also observed, such as glycogen storage and lipid accumulation. Finally, we collected and manually curated a cell state transition resource containing 224 time-course gene expression datasets and 153 cell types. Our framework, together with the data resource, is freely available at http://pc3t.idrug.net.cn/ . We believe that PC3T is a powerful tool to promote chemical-induced cell state transitions.

Chemical reprogramming for cell fate manipulation: Methods, applications, and perspectives

Chemical reprogramming offers an unprecedented opportunity to control somatic cell fate and generate desired cell types including pluripotent stem cells for applications in biomedicine in a precise, flexible, and controllable manner. Recent success in the chemical reprogramming of human somatic cells by activating a regeneration-like program provides an alternative way of producing stem cells for clinical translation. Likewise, chemical manipulation enables the capture of multiple (stem) cell states, ranging from totipotency to the stabilization of somatic fates in vitro. Here, we review progress in using chemical approaches for cell fate manipulation in addition to future opportunities in this promising field.

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.

Profiling A-to-I RNA editing during mouse somatic reprogramming at the single-cell level

Mouse somatic cells can be reprogrammed into induced pluripotent stem cells through a highly heterogeneous process regulated by numerous biological factors, including adenosine-to-inosine (A-to-I) RNA editing. In this study, we analyzed A-to-I RNA editing sites using a single-cell RNA sequencing (scRNA-seq) dataset with high-depth and full-length coverage. Our method revealed that A-to-I RNA editing frequency varied widely at the single-cell level and underwent dynamic changes. We also found that A-to-I RNA editing level was correlated with the expression of the RNA editing enzyme ADAR1. The analysis combined with gene ontology (GO) enrichment revealed that ADAR1-dependent A-to-I editing may downregulate the expression levels of Igtp, Irgm2, Mndal, Ifi202b, and Tapbp in the early stage, to inhibit the pathways of cellular response to interferon-beta and regulation of protein complex stability to promote mesenchymal-epithelial transition (MET). Notably, we identified a negative correlation between A-to-I RNA editing frequency and the expression of certain genes, such as Nras, Ube2l6, Zfp987, and Adsl.

Inferring cell diversity in single cell data using consortium-scale epigenetic data as a biological anchor for cell identity

Methods for cell clustering and gene expression from single-cell RNA sequencing (scRNA-seq) data are essential for biological interpretation of cell processes. Here, we present TRIAGE-Cluster which uses genome-wide epigenetic data from diverse bio-samples to identify genes demarcating cell diversity in scRNA-seq data. By integrating patterns of repressive chromatin deposited across diverse cell types with weighted density estimation, TRIAGE-Cluster determines cell type clusters in a 2D UMAP space. We then present TRIAGE-ParseR, a machine learning method which evaluates gene expression rank lists to define gene groups governing the identity and function of cell types. We demonstrate the utility of this two-step approach using atlases of in vivo and in vitro cell diversification and organogenesis. We also provide a web accessible dashboard for analysis and download of data and software. Collectively, genome-wide epigenetic repression provides a versatile strategy to define cell diversity and study gene regulation of scRNA-seq data.

Harnessing endogenous transcription factors directly by small molecules for chemically induced pluripotency inception

Chemistry-alone approach has recently been applied for incepting pluripotency in somatic cells, representing a breakthrough in biology. However, chemical reprogramming is hampered by low efficiency, and the underlying molecular mechanisms remain unclear. Particularly, chemical compounds do not have specific DNA-recognition domains or transcription regulatory domains, and then how do small molecules work as a driving force for reinstating pluripotency in somatic cells? Furthermore, how to efficiently clear materials and structures of an old cell to prepare the rebuilding of a new one? Here, we show that small molecule CD3254 activates endogenous existing transcription factor RXRα to significantly promote mouse chemical reprogramming. Mechanistically, CD3254-RXRα axis can directly activate all the 11 RNA exosome component genes (Exosc1-10 and Dis3) at transcriptional level. Unexpectedly, rather than degrading mRNAs as its substrates, RNA exosome mainly modulates the degradation of transposable element (TE)-associated RNAs, particularly MMVL30, which is identified as a new barrier for cell-fate determination. In turn, MMVL30-mediated inflammation (IFN-γ and TNF-α pathways) is reduced, contributing to the promotion of successful reprogramming. Collectively, our study provides conceptual advances for translating environmental cues into pluripotency inception, particularly, identifies that CD3254-RXRα-RNA exosome axis can promote chemical reprogramming, and suggests modulation of TE-mediated inflammation via CD3254-inducible RNA exosome as important opportunities for controlling cell fates and regenerative medicine.

A scalable sparse neural network framework for rare cell type annotation of single-cell transcriptome data

Automatic cell type annotation methods are increasingly used in single-cell RNA sequencing (scRNA-seq) analysis due to their fast and precise advantages. However, current methods often fail to account for the imbalance of scRNA-seq datasets and ignore information from smaller populations, leading to significant biological analysis errors. Here, we introduce scBalance, an integrated sparse neural network framework that incorporates adaptive weight sampling and dropout techniques for auto-annotation tasks. Using 20 scRNA-seq datasets with varying scales and degrees of imbalance, we demonstrate that scBalance outperforms current methods in both intra- and inter-dataset annotation tasks. Additionally, scBalance displays impressive scalability in identifying rare cell types in million-level datasets, as shown in the bronchoalveolar cell landscape. scBalance is also significantly faster than commonly used tools and comes in a user-friendly format, making it a superior tool for scRNA-seq analysis on the Python-based platform.

Cellular population dynamics shape the route to human pluripotency

Human cellular reprogramming to induced pluripotency is still an inefficient process, which has hindered studying the role of critical intermediate stages. Here we take advantage of high efficiency reprogramming in microfluidics and temporal multi-omics to identify and resolve distinct sub-populations and their interactions. We perform secretome analysis and single-cell transcriptomics to show functional extrinsic pathways of protein communication between reprogramming sub-populations and the re-shaping of a permissive extracellular environment. We pinpoint the HGF/MET/STAT3 axis as a potent enhancer of reprogramming, which acts via HGF accumulation within the confined system of microfluidics, and in conventional dishes needs to be supplied exogenously to enhance efficiency. Our data suggest that human cellular reprogramming is a transcription factor-driven process that it is deeply dependent on extracellular context and cell population determinants.

Characterizing Cellular Differentiation Potency and Waddington Landscape via Energy Indicator

The precise characterization of cellular differentiation potency remains an open question, which is fundamentally important for deciphering the dynamics mechanism related to cell fate transition. We quantitatively evaluated the differentiation potency of different stem cells based on the Hopfield neural network (HNN). The results emphasized that cellular differentiation potency can be approximated by Hopfield energy values. We then profiled the Waddington energy landscape of embryogenesis and cell reprogramming processes. The energy landscape at single-cell resolution further confirmed that cell fate decision is progressively specified in a continuous process. Moreover, the transition of cells from one steady state to another in embryogenesis and cell reprogramming processes was dynamically simulated on the energy ladder. These two processes can be metaphorized as the motion of descending and ascending ladders, respectively. We further deciphered the dynamics of the gene regulatory network (GRN) for driving cell fate transition. Our study proposes a new energy indicator to quantitatively characterize cellular differentiation potency without prior knowledge, facilitating the further exploration of the potential mechanism of cellular plasticity.

Sexually Dimorphic Alterations in the Transcriptome and Behavior with Loss of Histone Demethylase KDM5C

Chromatin dysregulation has emerged as a major hallmark of neurodevelopmental disorders such as intellectual disability (ID) and autism spectrum disorders (ASD). The prevalence of ID and ASD is higher in males compared to females, with unknown mechanisms. Intellectual developmental disorder, X-linked syndromic, Claes-Jensen type (MRXSCJ), is caused by loss-of-function mutations of lysine demethylase 5C (KDM5C), a histone H3K4 demethylase gene. KDM5C escapes X-inactivation, thereby presenting at a higher level in females. Initially, MRXSCJ was exclusively reported in males, while it is increasingly evident that females with heterozygous KDM5C mutations can show cognitive deficits. The mouse model of MRXSCJ, male Kdm5c-hemizygous knockout animals, recapitulates key features of human male patients. However, the behavioral and molecular traits of Kdm5c-heterozygous female mice remain incompletely characterized. Here, we report that gene expression and behavioral abnormalities are readily detectable in Kdm5c-heterozygous female mice, demonstrating the requirement for a higher KDM5C dose in females. Furthermore, we found both shared and sex-specific consequences of a reduced KDM5C dose in social behavior, gene expression, and genetic interaction with the counteracting enzyme KMT2A. These observations provide an essential insight into the sex-biased manifestation of neurodevelopmental disorders and sex chromosome evolution.

Retrospective identification of intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of the four "Yamanaka" factors OCT4, KLF4, SOX2, and MYC. However, even in homogeneous conditions, usually only a rare subset of cells admit reprogramming, and the molecular characteristics of this subset remain unknown. Here, we apply retrospective clone tracing to identify and characterize the individual human fibroblast cells that are primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis led to increased reprogramming efficiency, identifying it as a barrier to reprogramming. Changing the frequency of reprogramming by inhibiting the activity of LSD1 led to an enlarging of the pool of cells that were primed for reprogramming. Our results show that even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

RSL24D1 sustains steady-state ribosome biogenesis and pluripotency translational programs in embryonic stem cells

Embryonic stem cell (ESC) fate decisions are regulated by a complex circuitry that coordinates gene expression at multiple levels from chromatin to mRNA processing. Recently, ribosome biogenesis and translation have emerged as key pathways that efficiently control stem cell homeostasis, yet the underlying molecular mechanisms remain largely unknown. Here, we identified RSL24D1 as highly expressed in both mouse and human pluripotent stem cells. RSL24D1 is associated with nuclear pre-ribosomes and is required for the biogenesis of 60S subunits in mouse ESCs. Interestingly, RSL24D1 depletion significantly impairs global translation, particularly of key pluripotency factors and of components from the Polycomb Repressive Complex 2 (PRC2). While having a moderate impact on differentiation, RSL24D1 depletion significantly alters ESC self-renewal and lineage commitment choices. Altogether, these results demonstrate that RSL24D1-dependant ribosome biogenesis is both required to sustain the expression of pluripotent transcriptional programs and to silence PRC2-regulated developmental programs, which concertedly dictate ESC homeostasis.

Comparative roadmaps of reprogramming and oncogenic transformation identify Bcl11b and Atoh8 as broad regulators of cellular plasticity

Coordinated changes of cellular plasticity and identity are critical for pluripotent reprogramming and oncogenic transformation. However, the sequences of events that orchestrate these intermingled modifications have never been comparatively dissected. Here, we deconvolute the cellular trajectories of reprogramming (via Oct4/Sox2/Klf4/c-Myc) and transformation (via Ras/c-Myc) at the single-cell resolution and reveal how the two processes intersect before they bifurcate. This approach led us to identify the transcription factor Bcl11b as a broad-range regulator of cell fate changes, as well as a pertinent marker to capture early cellular intermediates that emerge simultaneously during reprogramming and transformation. Multiomics characterization of these intermediates unveiled a c-Myc/Atoh8/Sfrp1 regulatory axis that constrains reprogramming, transformation and transdifferentiation. Mechanistically, we found that Atoh8 restrains cellular plasticity, independent of cellular identity, by binding a specific enhancer network. This study provides insights into the partitioned control of cellular plasticity and identity for both regenerative and cancer biology.

Huyghe, Furlan et al. compare pluripotent reprogramming with oncogenic transformation and identify Bcl11b and Atoh8 as regulators of cellular plasticity in both processes, thus offering a unifying theory on the factors constraining cell fate changes.

Mettl14-driven senescence-associated secretory phenotype facilitates somatic cell reprogramming

The METTL3-METTL14 complex, the “writer” of N⁶-methyladenosine (m⁶A), plays an important role in many biological processes. Previous studies have shown that Mettl3 overexpression can increase the level of m⁶A and promote somatic cell reprogramming. Here, we demonstrate that Mettl14, another component of the methyltransferase complex, can significantly enhance the generation of induced pluripotent stem cells (iPSCs) in an m⁶A-independent manner. In cooperation with Oct4, Sox2, Klf4, and c-Myc, overexpressed Mettl14 transiently promoted senescence-associated secretory phenotype (SASP) gene expression in non-reprogrammed cells in the late stage of reprogramming. Subsequently, we demonstrated that interleukin-6 (IL-6), a component of the SASP, significantly enhanced somatic cell reprogramming. In contrast, blocking the SASP using a senolytic agent or a nuclear factor κB (NF-κB) inhibitor impaired the effect of Mettl14 on reprogramming. Our results highlight the m⁶A-independent function of Mettl14 in reprogramming and provide new insight into the interplay between senescence and reprogramming in vitro.

•

Mettl14 can facilitate somatic cell reprogramming in an m⁶A-independent manner

•

Mettl14 transcriptionally drives the senescence-associated secretory phenotype (SASP)

•

Mettl14-driven SASPs are mainly secreted from non-reprogramming cells

•

Blocking of SASP impairs the effect of Mettl14 on reprogramming

Diverse partial reprogramming strategies restore youthful gene expression and transiently suppress cell identity

Partial pluripotent reprogramming can reverse features of aging in mammalian cells, but the impact on somatic identity and the necessity of individual reprogramming factors remain unknown. Here, we used single-cell genomics to map the identity trajectory induced by partial reprogramming in multiple murine cell types and dissected the influence of each factor by screening all Yamanaka Factor subsets with pooled single-cell screens. We found that partial reprogramming restored youthful expression in adipogenic and mesenchymal stem cells but also temporarily suppressed somatic identity programs. Our pooled screens revealed that many subsets of the Yamanaka Factors both restore youthful expression and suppress somatic identity, but these effects were not tightly entangled. We also found that a multipotent reprogramming strategy inspired by amphibian regeneration restored youthful expression in myogenic cells. Our results suggest that various sets of reprogramming factors can restore youthful expression with varying degrees of somatic identity suppression. A record of this paper's Transparent Peer Review process is included in the supplemental information.

Practical bioinformatics pipelines for single-cell RNA-seq data analysis

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool to explore cells. With an increasing number of scRNA-seq data analysis tools that have been developed, it is challenging for users to choose and compare their performance. Here, we present an overview of the workflow for computational analysis of scRNA-seq data. We detail the steps of a typical scRNA-seq analysis, including experimental design, pre-processing and quality control, feature selection, dimensionality reduction, cell clustering and annotation, and downstream analysis including batch correction, trajectory inference and cell-cell communication. We provide guidelines according to our best practice. This review will be helpful for the experimentalists interested in analyzing their data, and will aid the users seeking to update their analysis pipelines.

Comparative parallel multi-omics analysis during the induction of pluripotent and trophectoderm states

Following fertilization, it is only at the 32-64-cell stage when a clear segregation between cells of the inner cell mass and trophectoderm is observed, suggesting a 'T'-shaped model of specification. Here, we examine whether the acquisition of these two states in vitro, by nuclear reprogramming, share similar dynamics/trajectories. Using a comparative parallel multi-omics analysis (i.e., bulk RNA-seq, scRNA-seq, ATAC-seq, ChIP-seq, RRBS and CNVs) on cells undergoing reprogramming to pluripotency and TSC state we show that each reprogramming system exhibits specific trajectories from the onset of the process, suggesting 'V'-shaped model. We describe in detail the various trajectories toward the two states and illuminate reprogramming stage-specific markers, blockers, facilitators and TSC subpopulations. Finally, we show that while the acquisition of the TSC state involves the silencing of embryonic programs by DNA methylation, during the acquisition of pluripotency these regions are initially defined but retain inactive by the elimination of H3K27ac.

RNA helicase DDX5 enables STAT1 mRNA translation and interferon signalling in hepatitis B virus replicating hepatocytes

OBJECTIVE

RNA helicase DDX5 is downregulated during HBV replication and poor prognosis HBV-related hepatocellular carcinoma (HCC). The objective of this study is to investigate the role of DDX5 in interferon (IFN) signalling. We provide evidence of a novel mechanism involving DDX5 that enables translation of transcription factor STAT1 mediating the IFN response.

DESIGN AND RESULTS

Molecular, pharmacological and biophysical assays were used together with cellular models of HBV replication, HCC cell lines and liver tumours. We demonstrate that DDX5 regulates STAT1 mRNA translation by resolving a G-quadruplex (rG4) RNA structure, proximal to the 5' end of STAT1 5'UTR. We employed luciferase reporter assays comparing wild type (WT) versus mutant rG4 sequence, rG4-stabilising compounds, CRISPR/Cas9 editing of the STAT1-rG4 sequence and circular dichroism determination of the rG4 structure. STAT1-rG4 edited cell lines were resistant to the effect of rG4-stabilising compounds in response to IFN-α, while HCC cell lines expressing low DDX5 exhibited reduced IFN response. Ribonucleoprotein and electrophoretic mobility assays demonstrated direct and selective binding of RNA helicase-active DDX5 to the WT STAT1-rG4 sequence. Immunohistochemistry of normal liver and liver tumours demonstrated that absence of DDX5 corresponded to absence of STAT1. Significantly, knockdown of DDX5 in HBV infected HepaRG cells reduced the anti-viral effect of IFN-α.

CONCLUSION

RNA helicase DDX5 resolves a G-quadruplex structure in 5'UTR of STAT1 mRNA, enabling STAT1 translation. We propose that DDX5 is a key regulator of the dynamic range of IFN response during innate immunity and adjuvant IFN-α therapy.

Epigenetics, Enhancer Function and 3D Chromatin Organization in Reprogramming to Pluripotency

Genome architecture, epigenetics and enhancer function control the fate and identity of cells. Reprogramming to induced pluripotent stem cells (iPSCs) changes the transcriptional profile and chromatin landscape of the starting somatic cell to that of the pluripotent cell in a stepwise manner. Changes in the regulatory networks are tightly regulated during normal embryonic development to determine cell fate, and similarly need to function in cell fate control during reprogramming. Switching off the somatic program and turning on the pluripotent program involves a dynamic reorganization of the epigenetic landscape, enhancer function, chromatin accessibility and 3D chromatin topology. Within this context, we will review here the current knowledge on the processes that control the establishment and maintenance of pluripotency during somatic cell reprogramming.

Mass Cytometry Study on Hepatic Fibrosis and Its Drug-Induced Recovery Using Mouse Peripheral Blood Mononuclear Cells

The number of patients with liver diseases has increased significantly with the progress of global industrialization. Hepatic fibrosis, one of the most common liver diseases diagnosed in many developed countries, occurs in response to chronic liver injury and is primarily driven by the development of inflammation. Earlier immunological studies have been focused on the importance of the innate immune response in the pathophysiology of steatohepatitis and fibrosis, but recently, it has also been reported that adaptive immunity, particularly B cells, plays an essential role in hepatic inflammation and fibrosis. However, despite recent data showing the importance of adaptive immunity, relatively little is known about the role of B cells in the pathogenesis of steatohepatitis fibrosis. In this study, a single-cell-based, high-dimensional mass cytometric investigation of the peripheral blood mononuclear cells collected from mice belonging to three groups [normal chow (NC), thioacetamide (TAA), and 11beta-HSD inhibitor drug] was conducted to further understand the pathogenesis of liver fibrosis through reliable noninvasive biomarkers. Firstly, major immune cell types and their population changes were qualitatively analyzed using UMAP dimensionality reduction and two-dimensional visualization technique combined with a conventional manual gating strategy. The population of B cells displayed a twofold increase in the TAA group compared to that in the NC group, which was recovered slightly after treatment with the 11beta-HSD inhibitor drug. In contrast, the populations of NK cells, effector CD4⁺ T cells, and memory CD8⁺ T cells were significantly reduced in the TAA group compared with those in the NC group. Further identification and quantification of the major immune cell types and their subsets were conducted based on automated clustering approaches [PhenoGraph (PG) and FlowSOM]. The B-cell subset corresponding to PhenoGraph cluster PG#2 (CD62L^highCD44^highLy6c^high B cells) and PG#3 (CD62L^highCD44^highLy6c^low B cell) appears to play a major role in both the development of hepatic fibrosis and recovery via treatment, whereas PG#1 (CD62L^lowCD44^highLy6c^low B cell) seems to play a dominant role in the development of hepatic fibrosis. These findings provide insights into the roles of cellular subsets of B cells during the progression of, and recovery from, hepatic fibrosis.

Human induced-T-to-natural killer cells have potent anti-tumour activities

BACKGROUND

Adoptive cell therapy (ACT) is a particularly promising area of cancer immunotherapy, engineered T and NK cells that express chimeric antigen receptors (CAR) are being explored for treating hematopoietic malignancies but exhibit limited clinical benefits for solid tumour patients, successful cellular immunotherapy of solid tumors demands new strategies.

METHODS

Inactivation of BCL11B were performed by CRISPR/Cas9 in human T cells. Immunophenotypic and transcriptional profiles of sgBCL11B T cells were characterized by cytometer and transcriptomics, respectively. sgBCL11B T cells are further engineered with chimeric antigen receptor. Anti-tumor activity of ITNK or CAR-ITNK cells were evaluated in preclinical and clinical studies.

RESULTS

We report that inactivation of BCL11B in human CD8⁺ and CD4⁺ T cells induced their reprogramming into induced T-to-natural killer cells (ITNKs). ITNKs contained a diverse TCR repertoire; downregulated T cell-associated genes such as TCF7 and LEF1; and expressed high levels of NK cell lineage-associated genes. ITNKs and chimeric antigen receptor (CAR)-transduced ITNKs selectively lysed a variety of cancer cells in culture and suppressed the growth of solid tumors in xenograft models. In a preliminary clinical study, autologous administration of ITNKs in patients with advanced solid tumors was well tolerated, and tumor stabilization was seen in six out nine patients, with one partial remission.

CONCLUSIONS

The novel ITNKs thus may be a promising novel cell source for cancer immunotherapy.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT03882840 . Registered 20 March 2019-Retrospectively registered.

Structurally-discovered KLF4 variants accelerate and stabilize reprogramming to pluripotency

Non-genetically modified somatic cells can only be inefficiently and stochastically reprogrammed to pluripotency by exogenous expression of reprogramming factors. Low competence of natural reprogramming factors may prevent the majority of cells to successfully and synchronously reprogram. Here we screened DNA-interacting amino acid residues in the zinc-finger domain of KLF4 for enhanced reprogramming efficiency using alanine-substitution scanning methods. Identified KLF4 L507A mutant accelerated and stabilized reprogramming to pluripotency in both mouse and human somatic cells. By testing all the variants of L507 position, variants with smaller amino acid residues in the KLF4 L507 position showed higher reprogramming efficiency. L507A bound more to promoters or enhancers of pluripotency genes, such as KLF5, and drove gene expression of these genes during reprogramming. Molecular dynamics simulations predicted that L507A formed additional interactions with DNA. Our study demonstrates how modifications in amino acid residues of DNA-binding domains enable next-generation reprogramming technology with engineered reprogramming factors.

•

KLF4 L507A variant accelerates and stabilizes reprogramming to pluripotency

•

KLF4 L507A has distinctive features of transcriptional binding and activation

•

KLF4 L507A may acquire a unique conformation with additional DNA interaction

•

Smaller amino acid residues in L507 position cause higher reprogramming efficiency

OUP accepted manuscript

Alternative polyadenylation increases transcript diversities at the 3’ end, regulating biological processes including cell differentiation, embryonic development and cancer progression. Here, we present a Bayesian method SCAPE, which enables de novo identification and quantification of polyadenylation (pA) sites at single-cell level by utilizing insert size information. We demonstrated its accuracy and robustness and identified 31 558 sites from 36 mouse organs, 43.8% (13 807) of which were novel. We illustrated that APA isoforms were associated with miRNAs binding and regulated in tissue-, cell type-and tumor-specific manners where no difference was found at gene expression level, providing an extra layer of information for cell clustering. Furthermore, we found genome-wide dynamic changes of APA usage during erythropoiesis and induced pluripotent stem cell (iPSC) differentiation, suggesting APA contributes to the functional flexibility and diversity of single cells. We expect SCAPE to aid the analyses of cellular dynamics and diversities in health and disease.

Enhanced chromatin accessibility contributes to X chromosome dosage compensation in mammals

BACKGROUND

Precise gene dosage of the X chromosomes is critical for normal development and cellular function. In mice, XX female somatic cells show transcriptional X chromosome upregulation of their single active X chromosome, while the other X chromosome is inactive. Moreover, the inactive X chromosome is reactivated during development in the inner cell mass and in germ cells through X chromosome reactivation, which can be studied in vitro by reprogramming of somatic cells to pluripotency. How chromatin processes and gene regulatory networks evolved to regulate X chromosome dosage in the somatic state and during X chromosome reactivation remains unclear.

RESULTS

Using genome-wide approaches, allele-specific ATAC-seq and single-cell RNA-seq, in female embryonic fibroblasts and during reprogramming to pluripotency, we show that chromatin accessibility on the upregulated mammalian active X chromosome is increased compared to autosomes. We further show that increased accessibility on the active X chromosome is erased by reprogramming, accompanied by erasure of transcriptional X chromosome upregulation and the loss of increased transcriptional burst frequency. In addition, we characterize gene regulatory networks during reprogramming and X chromosome reactivation, revealing changes in regulatory states. Our data show that ZFP42/REX1, a pluripotency-associated gene that evolved specifically in placental mammals, targets multiple X-linked genes, suggesting an evolutionary link between ZFP42/REX1, X chromosome reactivation, and pluripotency.

CONCLUSIONS

Our data reveal the existence of intrinsic compensatory mechanisms that involve modulation of chromatin accessibility to counteract X-to-Autosome gene dosage imbalances caused by evolutionary or in vitro X chromosome loss and X chromosome inactivation in mammalian cells.

Distinct human Langerhans cell subsets orchestrate reciprocal functions and require different developmental regulation

Langerhans cells (LCs) play a pivotal role in skin homeostasis, and the heterogeneity of LCs has long been considered. In this study, we have identified two steady-state (LC1 and LC2) and two activated LC subsets in the epidermis of human skin and in LCs derived from CD34⁺ hemopoietic stem cells (HSC-LCs) by utilizing single-cell RNA sequencing and mass cytometry. Analysis of HSC-LCs at multiple time-points during differentiation revealed that EGR1 and Notch signaling were among the top pathways regulating the bifurcation of LC1 and LC2. LC1 were characterized as classical LCs, mainly related to innate immunity and antigen processing. LC2 were similar to monocytes or myeloid dendritic cells, involving in immune responses and leukocyte activation. LC1 remained stable under inflammatory microenvironment, whereas LC2 were prone to being activated and demonstrated elevated expression of immuno-suppressive molecules. We revealed distinct human LC subsets that require different developmental regulation and orchestrate reciprocal functions.

Cell fate conversion prediction by group sparse optimization method utilizing single-cell and bulk OMICs data

Cell fate conversion by overexpressing defined factors is a powerful tool in regenerative medicine. However, identifying key factors for cell fate conversion requires laborious experimental efforts; thus, many of such conversions have not been achieved yet. Nevertheless, cell fate conversions found in many published studies were incomplete as the expression of important gene sets could not be manipulated thoroughly. Therefore, the identification of master transcription factors for complete and efficient conversion is crucial to render this technology more applicable clinically. In the past decade, systematic analyses on various single-cell and bulk OMICs data have uncovered numerous gene regulatory mechanisms, and made it possible to predict master gene regulators during cell fate conversion. By virtue of the sparse structure of master transcription factors and the group structure of their simultaneous regulatory effects on the cell fate conversion process, this study introduces a novel computational method predicting master transcription factors based on group sparse optimization technique integrating data from multi-OMICs levels, which can be applicable to both single-cell and bulk OMICs data with a high tolerance of data sparsity. When it is compared with current prediction methods by cross-referencing published and validated master transcription factors, it possesses superior performance. In short, this method facilitates fast identification of key regulators, give raise to the possibility of higher successful conversion rate and in the hope of reducing experimental cost.

Nuclear Transfer Arrest Embryos Show Massive Dysregulation of Genes Involved in Transcription Pathways

Somatic cell nuclear transfer (SCNT) technology can reprogram terminally differentiated cell nuclei into a totipotent state. However, the underlying molecular barriers of SCNT embryo development remain incompletely elucidated. Here, we observed that transcription-related pathways were incompletely activated in nuclear transfer arrest (NTA) embryos compared to normal SCNT embryos and in vivo fertilized (WT) embryos, which hinders the development of SCNT embryos. We further revealed the transcription pathway associated gene regulatory networks (GRNs) and found the aberrant transcription pathways can lead to the massive dysregulation of genes in NTA embryos. The predicted target genes of transcription pathways contain a series of crucial factors in WT embryos, which play an important role in catabolic process, pluripotency regulation, epigenetic modification and signal transduction. In NTA embryos, however, these genes were varying degrees of inhibition and show a defect in synergy. Overall, our research found that the incomplete activation of transcription pathways is another potential molecular barrier for SCNT embryos besides the incomplete reprogramming of epigenetic modifications, broadening the understanding of molecular mechanism of SCNT embryonic development.

The transcription factor code in iPSC reprogramming

Transcription factor (TF)-induced reprogramming of somatic cells across lineages and to induced pluripotent stem cells (iPSCs) has revealed a remarkable plasticity of differentiated cells and presents great opportunities for generating clinically relevant cell types for disease modeling and regenerative medicine. The understanding of iPSC reprogramming provides insights into the mechanisms that safeguard somatic cell identity, drive epigenetic reprogramming, and underlie cell fate specification in vivo. The combinatorial action of TFs has emerged as the key mechanism for the direct and indirect effects of reprogramming factors that induce the remodelling of the enhancer landscape. The interplay of TFs in iPSC reprogramming also yields trophectoderm- and extraembryonic endoderm-like cell populations, uncovering an intriguing plasticity of cell states and opening new avenues for exploring cell fate decisions during early embryogenesis.

Forkhead box family transcription factors as versatile regulators for cellular reprogramming to pluripotency

Forkhead box (Fox) transcription factors play important roles in mammalian development and disease. However, their function in mouse somatic cell reprogramming remains unclear. Here, we report that FoxD subfamily and FoxG1 accelerate induced pluripotent stem cells (iPSCs) generation from mouse fibroblasts as early as day4 while FoxA and FoxO subfamily impede this process obviously. More importantly, FoxD3, FoxD4 and FoxG1 can replace Oct4 respectively and generate iPSCs with germline transmission together with Sox2 and Klf4. On the contrary, FoxO6 almost totally blocks reprogramming through inhibiting cell proliferation, suppressing the expression of pluripotent genes and hindering the process of mesenchymal to epithelial transition (MET). Thus, our study uncovers unexpected roles of Fox transcription factors in reprogramming and offers new insights into cell fate transition.

Cattle In Vitro Induced Pluripotent Stem Cells Generated and Maintained in 5 or 20% Oxygen and Different Supplementation

The event of cellular reprogramming into pluripotency is influenced by several factors, such as in vitro culture conditions (e.g., culture medium and oxygen concentration). Herein, bovine iPSCs (biPSCs) were generated in different levels of oxygen tension (5% or 20% of oxygen) and supplementation (bFGF or bFGF + LIF + 2i-bFL2i) to evaluate the efficiency of pluripotency induction and maintenance in vitro. Initial reprogramming was observed in all groups and bFL2i supplementation initially resulted in a superior number of colonies. However, bFL2i supplementation in low oxygen led to a loss of self-renewal and pluripotency maintenance. All clonal lines were positive for alkaline phosphatase; they expressed endogenous pluripotency-related genes SOX2, OCT4 and STELLA. However, expression was decreased throughout the passages without the influence of oxygen tension. GLUT1 and GLUT3 were upregulated by low oxygen. The biPSCs were immunofluorescence-positive stained for OCT4 and SOX2 and they formed embryoid bodies which differentiated in ectoderm and mesoderm (all groups), as well as endoderm (one line from bFL2i in high oxygen). Our study is the first to compare high and low oxygen environments during and after induced reprogramming in cattle. In our conditions, a low oxygen environment did not favor the pluripotency maintenance of biPSCs.

Single-Cell Transcriptome Analysis as a Promising Tool to Study Pluripotent Stem Cell Reprogramming

Cells are the basic units of all organisms and are involved in all vital activities, such as proliferation, differentiation, senescence, and apoptosis. A human body consists of more than 30 trillion cells generated through repeated division and differentiation from a single-cell fertilized egg in a highly organized programmatic fashion. Since the recent formation of the Human Cell Atlas consortium, establishing the Human Cell Atlas at the single-cell level has been an ongoing activity with the goal of understanding the mechanisms underlying diseases and vital cellular activities at the level of the single cell. In particular, transcriptome analysis of embryonic stem cells at the single-cell level is of great importance, as these cells are responsible for determining cell fate. Here, we review single-cell analysis techniques that have been actively used in recent years, introduce the single-cell analysis studies currently in progress in pluripotent stem cells and reprogramming, and forecast future studies.

Reprogramming lineage identity through cell-cell fusion

The conversion of differentiated cells to a pluripotent state through somatic cell nuclear transfer provided the first unequivocal evidence that differentiation was reversible. In more recent times, introducing a combination of key transcription factors into terminally differentiated mammalian cells was shown to drive their conversion to induced pluripotent stem cells (iPSCs). These discoveries were transformative, but the relatively slow speed (2-3 weeks) and low efficiency of reprogramming (0.1-1%) made deciphering the underlying molecular mechanisms difficult and complex. Cell fusion provides an alternative reprogramming approach that is both efficient and tractable, particularly when combined with modern multi-omics analysis of individual cells. Here we review the history and the recent advances in cell-cell fusion that are enabling a better understanding cell fate conversion, and we discuss how this knowledge could be used to shape improved strategies for regenerative medicine.

Modular, efficient and constant-memory single-cell RNA-seq preprocessing

We describe a workflow for preprocessing of single-cell RNA-sequencing data that balances efficiency and accuracy. Our workflow is based on the kallisto and bustools programs, and is near optimal in speed with a constant memory requirement providing scalability for arbitrarily large datasets. The workflow is modular, and we demonstrate its flexibility by showing how it can be used for RNA velocity analyses.

Identifying transposable element expression dynamics and heterogeneity during development at the single-cell level with a processing pipeline scTE

Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity in unclear ways. Single-cell sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE expression has not been addressed. Here, we develop a single-cell TE processing pipeline, scTE, and report the expression of TEs in single cells in a range of biological contexts. Specific TE types are expressed in subpopulations of embryonic stem cells and are dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs are expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we apply scTE to single-cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results classify the dynamic patterns of TEs in single cells and their contributions to cell heterogeneity.

How transposable elements (TE) contribute to cell fate changes is unclear. Here, the authors generate a pipeline to quantify TE expression from single cell data. They show the dynamic expression of TEs from gastrulation to somatic cell reprogramming and human disease

Targeted expression profiling reveals distinct stages of early canine fibroblast reprogramming are regulated by 2-oxoglutarate hydroxylases

Background

Ectopic expression of a defined set of transcription factors allows the reprogramming of mammalian somatic cells to pluripotency. Despite continuous progress in primate and rodent reprogramming, limited attention has been paid to cell reprogramming in domestic and companion species. Previous studies attempting to reprogram canine cells have mostly assessed a small number of presumptive canine induced pluripotent stem cell (iPSC) lines for generic pluripotency attributes. However, why canine cell reprogramming remains extremely inefficient is poorly understood.

Methods

To better characterize the initial steps of pluripotency induction in canine somatic cells, we optimized an experimental system where canine fetal fibroblasts (cFFs) are transduced with the Yamanaka reprogramming factors by Sendai virus vectors. We use quantitative PCR arrays to measure the expression of 80 target genes at various stages of canine cell reprogramming. We ask how cFF reprogramming is influenced by small molecules affecting the epigenomic modification 5-hydroxymethylcytosine, specifically L-ascorbic acid and retinoic acid (AA/RA).

Results

We found that the expression and catalytic output of a class of 2-oxoglutarate-dependent (2-OG) hydroxylases, known as ten-eleven translocation (TET) enzymes, can be modulated in canine cells treated with AA/RA. We further show that AA/RA treatment induces TET1 expression and facilitates early canine reprogramming, evidenced by upregulation of epithelial and pluripotency markers. Using a chemical inhibitor of 2-OG hydroxylases, we demonstrate that 2-OG hydroxylase activity regulates the expression of a subset of genes involved in mesenchymal-to-epithelial transition (MET) and pluripotency in early canine reprogramming. We identify a set of transcription factors depleted in maturing reprogramming intermediates compared to pluripotent canine embryonic stem cells.

Conclusions

Our findings highlight 2-OG hydroxylases have evolutionarily conserved and divergent functions regulating the early reprogramming of canine somatic cells and show reprogramming conditions can be rationally optimized for the generation of maturing canine iPSC.

IRF-1 expressed in the inner cell mass of the porcine early blastocyst enhances the pluripotency of induced pluripotent stem cells

BACKGROUND

Despite years of research, porcine-induced pluripotent stem cells (piPSCs) with germline chimeric capacity have not been established. Furthermore, the key transcription factors (TFs) defining the naïve state in piPSCs also remain elusive, even though TFs in the inner cell mass (ICM) are believed to be key molecular determinants of naïve pluripotency. In this study, interferon regulatory factor 1 (IRF-1) was screened to express higher in ICM than trophectoderm (TE). But the impact of IRF-1 on maintenance of pluripotency in piPSCs was not determined.

METHODS

Transcriptome profiles of the early ICM were analyzed to determine highly interconnected TFs. Cells carrying these TFs' reporter were used to as donor cells for somatic cell nuclear transfer to detect expression patterns in blastocysts. Next, IRF1-Flag was overexpressed in DOX-hLIF-2i piPSCs and AP staining, qRT-PCR, and RNA-seq were conducted to examine the effect of IRF-1 on pluripotency. Then, the expression of IRF-1 in DOX-hLIF-2i piPSCs was labeled by GFP and qRT-PCR was conducted to determine the difference between GFP-positive and GFP-negative cells. Next, ChIP-Seq was conducted to identify genes target by IRF-1. Treatment with IL7 in wild-type piPSCs and STAT3 phosphorylation inhibitor in IRF-1 overexpressing piPSCs was conducted to confirm the roles of JAK-STAT3 signaling pathway in IRF-1's regulation of pluripotency. Moreover, during reprogramming, IRF-1 was overexpressed and knocked down to determine the change of reprogramming efficiency.

RESULTS

IRF-1 was screened to be expressed higher in porcine ICM than TE of d6~7 SCNT blastocysts. First, overexpression of IRF-1 in the piPSCs was observed to promote the morphology, AP staining, and expression profiles of pluripotency genes as would be expected when cells approach the naïve state. Genes, KEGG pathways, and GO terms related to the process of differentiation were also downregulated. Next, in the wild-type piPSCs, high-level fluorescence activated by the IRF-1 promoter was associated with higher expression of naïve related genes in piPSCs. Analysis by ChIP-Seq indicated that genes related to the JAK-STAT pathway, and expression of IL7 and STAT3 were activated by IRF-1. The inhibitor of STAT3 phosphorylation was observed could revert the expression of primed genes in IRF-1 overexpressing cells, but the addition of IL7 in culture medium had no apparent change in the cell morphology, AP staining results, or expression of pluripotency related genes. In addition, knockdown of IRF-1 during reprogramming appeared to reduce reprogramming efficiency, whereas overexpression exerted the converse effect.

CONCLUSION

The IRF-1 expressed in the ICM of pigs' early blastocyst enhances the pluripotency of piPSCs, in part through promoting the JAK-STAT pathway.

Ancestry-dependent gene expression correlates with reprogramming to pluripotency and multiple dynamic biological processes

Induced pluripotent stem cells (iPSCs) can be derived from differentiated cells, enabling the generation of personalized disease models by differentiating patient-derived iPSCs into disease-relevant cell lines. While genetic variability between different iPSC lines affects differentiation potential, how this variability in somatic cells affects pluripotent potential is less understood. We generated and compared transcriptomic data from 72 dermal fibroblast-iPSC pairs with consistent variation in reprogramming efficiency. By considering equal numbers of samples from self-reported African Americans and White Americans, we identified both ancestry-dependent and ancestry-independent transcripts associated with reprogramming efficiency, suggesting that transcriptomic heterogeneity can substantially affect reprogramming. Moreover, reprogramming efficiency-associated genes are involved in diverse dynamic biological processes, including cancer and wound healing, and are predictive of 5-year breast cancer survival in an independent cohort. Candidate genes may provide insight into mechanisms of ancestry-dependent regulation of cell fate transitions and motivate additional studies for improvement of reprogramming.

JMJD3 acts in tandem with KLF4 to facilitate reprogramming to pluripotency

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3’s catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription.

Previous work suggested that histone demethylase JMJD3 is detrimental to somatic cell reprogramming. Here, the authors show that while JMJD3 has a context-independent detrimental effect on early stages of reprogramming, during late stages it activates epithelial and pluripotency genes together with Klf4.

Dynamical gene regulatory networks are tuned by transcriptional autoregulation with microRNA feedback

Concepts from dynamical systems theory, including multi-stability, oscillations, robustness and stochasticity, are critical for understanding gene regulation during cell fate decisions, inflammation and stem cell heterogeneity. However, the prevalence of the structures within gene networks that drive these dynamical behaviours, such as autoregulation or feedback by microRNAs, is unknown. We integrate transcription factor binding site (TFBS) and microRNA target data to generate a gene interaction network across 28 human tissues. This network was analysed for motifs capable of driving dynamical gene expression, including oscillations. Identified autoregulatory motifs involve 56% of transcription factors (TFs) studied. TFs that autoregulate have more interactions with microRNAs than non-autoregulatory genes and 89% of autoregulatory TFs were found in dual feedback motifs with a microRNA. Both autoregulatory and dual feedback motifs were enriched in the network. TFs that autoregulate were highly conserved between tissues. Dual feedback motifs with microRNAs were also conserved between tissues, but less so, and TFs regulate different combinations of microRNAs in a tissue-dependent manner. The study of these motifs highlights ever more genes that have complex regulatory dynamics. These data provide a resource for the identification of TFs which regulate the dynamical properties of human gene expression.