Diversification of reprogramming trajectories revealed by parallel single-cell transcriptome and chromatin accessibility sequencing

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.

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.

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.

Single-nucleus multiomic mapping of m6A methylomes and transcriptomes in native populations of cells with sn-m6A-CT

N6-methyladenosine (m6A) RNA modification plays important roles in the governance of gene expression and is temporally regulated in different cell states. In contrast to global m6A profiling in bulk sequencing, single-cell technologies for analyzing m6A heterogeneity are not extensively established. Here, we developed single-nucleus m6A-CUT&Tag (sn-m6A-CT) for simultaneous profiling of m6A methylomes and transcriptomes within a single nucleus using mouse embryonic stem cells (mESCs). m6A-CT is capable of enriching m6A-marked RNA molecules in situ, without isolating RNAs from cells. We adapted m6A-CT to the droplet-based single-cell omics platform and demonstrated high-throughput performance in analyzing nuclei isolated from thousands of cells from various cell types. We show that sn-m6A-CT profiling is sufficient to determine cell identity and allows the generation of cell-type-specific m6A methylome landscapes from heterogeneous populations. These indicate that sn-m6A-CT provides additional dimensions to multimodal datasets and insights into epitranscriptomic landscape in defining cell fate identity and states.

Transient naive reprogramming corrects hiPS cells functionally and epigenetically

Cells undergo a major epigenome reconfiguration when reprogrammed to human induced pluripotent stem cells (hiPS cells). However, the epigenomes of hiPS cells and human embryonic stem (hES) cells differ significantly, which affects hiPS cell function1-8. These differences include epigenetic memory and aberrations that emerge during reprogramming, for which the mechanisms remain unknown. Here we characterized the persistence and emergence of these epigenetic differences by performing genome-wide DNA methylation profiling throughout primed and naive reprogramming of human somatic cells to hiPS cells. We found that reprogramming-induced epigenetic aberrations emerge midway through primed reprogramming, whereas DNA demethylation begins early in naive reprogramming. Using this knowledge, we developed a transient-naive-treatment (TNT) reprogramming strategy that emulates the embryonic epigenetic reset. We show that the epigenetic memory in hiPS cells is concentrated in cell of origin-dependent repressive chromatin marked by H3K9me3, lamin-B1 and aberrant CpH methylation. TNT reprogramming reconfigures these domains to a hES cell-like state and does not disrupt genomic imprinting. Using an isogenic system, we demonstrate that TNT reprogramming can correct the transposable element overexpression and differential gene expression seen in conventional hiPS cells, and that TNT-reprogrammed hiPS and hES cells show similar differentiation efficiencies. Moreover, TNT reprogramming enhances the differentiation of hiPS cells derived from multiple cell types. Thus, TNT reprogramming corrects epigenetic memory and aberrations, producing hiPS cells that are molecularly and functionally more similar to hES cells than conventional hiPS cells. We foresee TNT reprogramming becoming a new standard for biomedical and therapeutic applications and providing a novel system for studying epigenetic memory.

Chemical-induced epigenome resetting for regeneration program activation in human cells

Human somatic cells can be reprogrammed to pluripotent stem cells by small molecules through an intermediate stage with a regeneration signature, but how this regeneration state is induced remains largely unknown. Here, through integrated single-cell analysis of transcriptome, we demonstrate that the pathway of human chemical reprogramming with regeneration state is distinct from that of transcription-factor-mediated reprogramming. Time-course construction of chromatin landscapes unveils hierarchical histone modification remodeling underlying the regeneration program, which involved sequential enhancer recommissioning and mirrored the reversal process of regeneration potential lost in organisms as they mature. In addition, LEF1 is identified as a key upstream regulator for regeneration gene program activation. Furthermore, we reveal that regeneration program activation requires sequential enhancer silencing of somatic and proinflammatory programs. Altogether, chemical reprogramming resets the epigenome through reversal of the loss of natural regeneration, representing a distinct concept for cellular reprogramming and advancing the development of regenerative therapeutic strategies.

Transfer learning enables predictions in network biology

Mapping gene networks requires large amounts of transcriptomic data to learn the connections between genes, which impedes discoveries in settings with limited data, including rare diseases and diseases affecting clinically inaccessible tissues. Recently, transfer learning has revolutionized fields such as natural language understanding1,2 and computer vision3 by leveraging deep learning models pretrained on large-scale general datasets that can then be fine-tuned towards a vast array of downstream tasks with limited task-specific data. Here, we developed a context-aware, attention-based deep learning model, Geneformer, pretrained on a large-scale corpus of about 30 million single-cell transcriptomes to enable context-specific predictions in settings with limited data in network biology. During pretraining, Geneformer gained a fundamental understanding of network dynamics, encoding network hierarchy in the attention weights of the model in a completely self-supervised manner. Fine-tuning towards a diverse panel of downstream tasks relevant to chromatin and network dynamics using limited task-specific data demonstrated that Geneformer consistently boosted predictive accuracy. Applied to disease modelling with limited patient data, Geneformer identified candidate therapeutic targets for cardiomyopathy. Overall, Geneformer represents a pretrained deep learning model from which fine-tuning towards a broad range of downstream applications can be pursued to accelerate discovery of key network regulators and candidate therapeutic targets.

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.

Manipulation of the nucleoscaffold potentiates cellular reprogramming kinetics

Somatic cell fate is an outcome set by the activities of specific transcription factors and the chromatin landscape and is maintained by gene silencing of alternate cell fates through physical interactions with the nuclear scaffold. Here, we evaluate the role of the nuclear scaffold as a guardian of cell fate in human fibroblasts by comparing the effects of transient loss (knockdown) and mutation (progeria) of functional Lamin A/C, a core component of the nuclear scaffold. We observed that Lamin A/C deficiency or mutation disrupts nuclear morphology, heterochromatin levels, and increases access to DNA in lamina-associated domains. Changes in Lamin A/C were also found to impact the mechanical properties of the nucleus when measured by a microfluidic cellular squeezing device. We also show that transient loss of Lamin A/C accelerates the kinetics of cellular reprogramming to pluripotency through opening of previously silenced heterochromatin domains while genetic mutation of Lamin A/C into progerin induces a senescent phenotype that inhibits the induction of reprogramming genes. Our results highlight the physical role of the nuclear scaffold in safeguarding cellular fate.

Ribosomal proteins regulate 2-cell-stage transcriptome in mouse embryonic stem cells

A rare sub-population of mouse embryonic stem cells (mESCs), the 2-cell-like cell, is defined by the expression of MERVL and 2-cell-stage-specific transcript (2C transcript). Here, we report that the ribosomal proteins (RPs) RPL14, RPL18, and RPL23 maintain the identity of mESCs and regulate the expression of 2C transcripts. Disregulation of the RPs induces DUX-dependent expression of 2C transcripts and alters the chromatin landscape. Mechanically, knockdown (KD) of RPs triggers the binding of RPL11 to MDM2, an interaction known to prevent P53 protein degradation. Increased P53 protein upon RP KD further activates its downstream pathways, including DUX. Our study delineates the critical roles of RPs in 2C transcript activation, ascribing a novel function to these essential proteins.

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.

Gene regulatory network reconfiguration in direct lineage reprogramming

In direct lineage conversion, transcription factor (TF) overexpression reconfigures gene regulatory networks (GRNs) to reprogram cell identity. We previously developed CellOracle, a computational method to infer GRNs from single-cell transcriptome and epigenome data. Using inferred GRNs, CellOracle simulates gene expression changes in response to TF perturbation, enabling in silico interrogation of network reconfiguration. Here, we combine CellOracle analysis with lineage tracing of fibroblast to induced endoderm progenitor (iEP) conversion, a prototypical direct reprogramming paradigm. By linking early network state to reprogramming outcome, we reveal distinct network configurations underlying successful and failed fate conversion. Via in silico simulation of TF perturbation, we identify new factors to coax cells into successfully converting their identity, uncovering a central role for the AP-1 subunit Fos with the Hippo signaling effector, Yap1. Together, these results demonstrate the efficacy of CellOracle to infer and interpret cell-type-specific GRN configurations, providing new mechanistic insights into lineage reprogramming.

Latent disease similarities and therapeutic repurposing possibilities uncovered by multi-modal generative topic modeling of human diseases

Motivation

Human diseases are characterized by multiple features such as their pathophysiological, molecular and genetic changes. The rapid expansion of such multi-modal disease-omics space provides an opportunity to re-classify diverse human diseases and to uncover their latent molecular similarities, which could be exploited to repurpose a therapeutic-target for one disease to another.

Results

Herein, we probe this underexplored space by soft-clustering 6955 human diseases by multi-modal generative topic modeling. Focusing on chronic kidney disease and myocardial infarction, two most life-threatening diseases, unveiled are their previously underrecognized molecular similarities to neoplasia and mental/neurological-disorders, and 69 repurposable therapeutic-targets for these diseases. Using an edit-distance-based pathway-classifier, we also find molecular pathways by which these targets could elicit their clinical effects. Importantly, for the 17 targets, the evidence for their therapeutic usefulness is retrospectively found in the pre-clinical and clinical space, illustrating the effectiveness of the method, and suggesting its broader applications across diverse human diseases.

Availability and implementation

The code reported in this article is available at: https://github.com/skozawa170301ktx/MultiModalDiseaseModeling.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Single-cell RNA sequencing reveals cellular and molecular reprograming landscape of gliomas and lung cancer brain metastases

BACKGROUND

Brain malignancies encompass gliomas and brain metastases originating from extracranial tumours including lung cancer. Approximately 50% of patients with lung adenocarcinoma (LUAD) will eventually develop brain metastases. However, the specific characteristics of gliomas and lung-to-brain metastases (LC) are largely unknown.

METHODS

We applied single-cell RNA sequencing to profile immune and nonimmune cells in 4 glioma and 10 LC samples.

RESULTS

Our analysis revealed that tumour microenvironment (TME) cells are present in heterogeneous subpopulations. LC reprogramed cells into immune suppressed state, including microglia, macrophages, endothelial cells, and CD8⁺ T cells, with unique cell proportions and gene signatures. Particularly, we identified that a subset of macrophages was associated with poor prognosis. ROS (reactive oxygen species)-producing neutrophils was found to participant in angiogenesis. Furthermore, endothelial cells participated in active communication with fibroblasts. Metastatic epithelial cells exhibited high heterogeneity in chromosomal instability (CIN) and cell population.

CONCLUSIONS

Our findings provide a comprehensive understanding of the heterogenicity of the tumor microenvironment and tumour cells and it will be crucial for successful immunotherapy development for brain metastasis of lung cancer.

Transcriptome Analysis Reveals Critical Factors For Survival After Adenovirus Serotype 4 Infection

Fowl adenovirus serotype-4 (FAdV-4) is highly lethal to poultry, making it one of the leading causes of economic losses in the poultry industry. However, a small proportion of poultry can survive after FAdV-4 infection. It is unclear whether there are genetic factors that protect chickens from FAdV-4 infection. Therefore, the livers from chickens uninfected with FAdV-4 (Normal), dead after FAdV-4 infection (Dead) or surviving after FAdV-4 infection (Survivor) were collected for RNA-seq, and 2,649 differentially expressed genes (DEGs) were identified. Among these, many immune-related cytokines and chemokines were significantly upregulated in the Dead group compared with the Survivor group, which might indicate that death is related to an excessive inflammatory immune response (cytokine storm). Subsequently, the KEGG results for DEGs specifically expressed in each comparison group indicated that cell cycle and apoptosis-related DEGs were upregulated and metabolism-related DEGs were downregulated in the Dead group, which also validated the reliability of the samples. Furthermore, GO and KEGG results showed DEGs expressed in all three groups were mainly associated with cell cycle. Among them, BRCA1, CDK1, ODC1, and MCM3 were screened as factors that might influence FAdV-4 infection. The qPCR results demonstrated that these 4 factors were not only upregulated in the Dead group but also significantly upregulated in the LMH cells after 24 h infection by FAdV-4. Moreover, interfering with BRCA1, CDK1, ODC1, and MCM3 significantly attenuated viral replication of FAdV-4. And interfering of BRCA1, CDK1, and MCM3 had more substantial hindering effects. These results provided novel insights into the molecular changes following FAdV-4 infection but also shed light on potential factors driving the survival of FAdV-4 infection in chickens.

SETDB1 acts as a topological accessory to Cohesin via an H3K9me3-independent, genomic shunt for regulating cell fates

SETDB1 is a key regulator of lineage-specific genes and endogenous retroviral elements (ERVs) through its deposition of repressive H3K9me3 mark. Apart from its H3K9me3 regulatory role, SETDB1 has seldom been studied in terms of its other potential regulatory roles. To investigate this, a genomic survey of SETDB1 binding in mouse embryonic stem cells across multiple libraries was conducted, leading to the unexpected discovery of regions bereft of common repressive histone marks (H3K9me3, H3K27me3). These regions were enriched with the CTCF motif that is often associated with the topological regulator Cohesin. Further profiling of these non-H3K9me3 regions led to the discovery of a cluster of non-repeat loci that were co-bound by SETDB1 and Cohesin. These regions, which we named DiSCs (domains involving SETDB1 and Cohesin) were seen to be proximal to the gene promoters involved in embryonic stem cell pluripotency and lineage development. Importantly, it was found that SETDB1-Cohesin co-regulate target gene expression and genome topology at these DiSCs. Depletion of SETDB1 led to localized dysregulation of Cohesin binding thereby locally disrupting topological structures. Dysregulated gene expression trends revealed the importance of this cluster in ES cell maintenance as well as at gene 'islands' that drive differentiation to other lineages. The 'unearthing' of the DiSCs thus unravels a unique topological and transcriptional axis of control regulated chiefly by SETDB1.

A reference single-cell regulomic and transcriptomic map of cynomolgus monkeys

Non-human primates are attractive laboratory animal models that accurately reflect both developmental and pathological features of humans. Here we present a compendium of cell types across multiple organs in cynomolgus monkeys (Macaca fascicularis) using both single-cell chromatin accessibility and RNA sequencing data. The integrated cell map enables in-depth dissection and comparison of molecular dynamics, cell-type compositions and cellular heterogeneity across multiple tissues and organs. Using single-cell transcriptomic data, we infer pseudotime cell trajectories and cell-cell communications to uncover key molecular signatures underlying their cellular processes. Furthermore, we identify various cell-specific cis-regulatory elements and construct organ-specific gene regulatory networks at the single-cell level. Finally, we perform comparative analyses of single-cell landscapes among mouse, monkey and human. We show that cynomolgus monkey has strikingly higher degree of similarities in terms of immune-associated gene expression patterns and cellular communications to human than mouse. Taken together, our study provides a valuable resource for non-human primate cell biology.

Non-human primates are attractive laboratory animal models that can accurately reflect some developmental and pathological features of humans. Here the authors chart a reference cell map of cynomolgus monkeys using both scATAC-seq and scRNA-seq data across multiple organs, providing insights into the molecular dynamics and cellular heterogeneity of this organism.

Delineating chromatin accessibility re-patterning at single cell level during early stage of direct cardiac reprogramming

Direct conversion of cardiac fibroblast into induced cardiomyocytes (iCMs) by forced expression of cardiac transcription factors, such as Mef2c, Gata4, and Tbx5 (MGT), holds great promise for regenerative medicine. The process of cardiac reprogramming consists of waves of transcriptome remodelling events. However, how this transcriptome remodelling is driven by the upstream chromatin landscape alteration is still unclear. In this study, we performed single-cell ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) on early reprogramming iCMs given the known epigenetic changes as early as day 3. This approach unveiled networks of transcription factors (TFs) involved in the early shift of chromatin accessibility during cardiac reprogramming. Combining our analysis with functional assays, we identified Smad3 to be a bimodal TF in cardiac reprogramming, a barrier in the initiation of reprogramming and a facilitator during the intermediate stage of reprogramming. Moreover, integrative analysis of scATAC-seq with scRNA-seq data led to the identification of active TFs important for iCM conversion. Finally, we discovered a global rewiring of cis-regulatory interactions of cardiac genes along the reprogramming trajectory. Collectively, our scATAC-seq study and the integrative analysis with scRNA-seq data provided valuable resources to understand the epigenomic heterogeneity and its alteration in relation to transcription changes during early stage of cardiac reprogramming.

MYCL promotes iPSC-like colony formation via MYC Box 0 and 2 domains

Human induced pluripotent stem cells (hiPSCs) can differentiate into cells of the three germ layers and are promising cell sources for regenerative medicine therapies. However, current protocols generate hiPSCs with low efficiency, and the generated iPSCs have variable differentiation capacity among different clones. Our previous study reported that MYC proteins (c-MYC and MYCL) are essential for reprogramming and germline transmission but that MYCL can generate hiPSC colonies more efficiently than c-MYC. The molecular underpinnings for the different reprogramming efficiencies between c-MYC and MYCL, however, are unknown. In this study, we found that MYC Box 0 (MB0) and MB2, two functional domains conserved in the MYC protein family, contribute to the phenotypic differences and promote hiPSC generation in MYCL-induced reprogramming. Proteome analyses suggested that in MYCL-induced reprogramming, cell adhesion-related cytoskeletal proteins are regulated by the MB0 domain, while the MB2 domain regulates RNA processes. These findings provide a molecular explanation for why MYCL has higher reprogramming efficiency than c-MYC.

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.

Panobinostat Effectively Increases Histone Acetylation and Alters Chromatin Accessibility Landscape in Canine Embryonic Fibroblasts but Does Not Enhance Cellular Reprogramming

Robust and reproducible protocols to efficiently reprogram adult canine cells to induced pluripotent stem cells are still elusive. Somatic cell reprogramming requires global chromatin remodeling that is finely orchestrated spatially and temporally. Histone acetylation and deacetylation are key regulators of chromatin condensation, mediated by histone acetyltransferases and histone deacetylases (HDACs), respectively. HDAC inhibitors have been used to increase histone acetylation, chromatin accessibility, and somatic cell reprogramming in human and mice cells. We hypothesized that inhibition of HDACs in canine fibroblasts would increase their reprogramming efficiency by altering the epigenomic landscape and enabling greater chromatin accessibility. We report that a combined treatment of panobinostat (LBH589) and vitamin C effectively inhibits HDAC function and increases histone acetylation in canine embryonic fibroblasts in vitro, with no significant cytotoxic effects. We further determined the effect of this treatment on global chromatin accessibility via Assay for Transposase-Accessible Chromatin using sequencing. Finally, the treatment did not induce any significant increase in cellular reprogramming efficiency. Although our data demonstrate that the unique epigenetic landscape of canine cells does not make them amenable to cellular reprogramming through the proposed treatment, it provides a rationale for a targeted, canine-specific, reprogramming approach by enhancing the expression of transcription factors such as CEBP.

Multi-species single-cell transcriptomic analysis of ocular compartment regulons

The retina is a widely profiled tissue in multiple species by single-cell RNA sequencing studies. However, integrative research of the retina across species is lacking. Here, we construct the first single-cell atlas of the human and porcine ocular compartments and study inter-species differences in the retina. In addition to that, we identify putative adult stem cells present in the iris tissue. We also create a disease map of genes involved in eye disorders across compartments of the eye. Furthermore, we probe the regulons of different cell populations, which include transcription factors and receptor-ligand interactions and reveal unique directional signalling between ocular cell types. In addition, we study conservation of regulons across vertebrates and zebrafish to identify common core factors. Here, we show perturbation of KLF7 gene expression during retinal ganglion cells differentiation and conclude that it plays a significant role in the maturation of retinal ganglion cells.

A comprehensive analysis of the ocular networks among various tissues is necessary to understand eye physiology in health and disease. Here the authors present a multi-species single-cell transcriptomic atlas consisting of cells of the cornea, iris, ciliary body, neural retina, retinal pigmented epithelium, and choroid.

Identification of New Transcription Factors that Can Promote Pluripotent Reprogramming

BACKGROUND

Four transcription factors, Oct4, Sox2, Klf4, and c-Myc (the Yamanka factors), can reprogram somatic cells to induced pluripotent stem cells (iPSCs). Many studies have provided a number of alternative combinations to the non-Yamanaka factors. However, it is clear that many additional transcription factors that can generate iPSCs remain to be discovered.

METHODS

The chromatin accessibility and transcriptional level of human embryonic stem cells and human urine cells were compared by Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing (RNA-seq) to identify potential reprogramming factors. Selected transcription factors were employed to reprogram urine cells, and the reprogramming efficiency was measured. Urine-derived iPSCs were detected for pluripotency by Immunofluorescence, quantitative polymerase chain reaction, RNA sequencing and teratoma formation test. Finally, we assessed the differentiation potential of the new iPSCs to cardiomyocytes in vitro.

RESULTS

ATAC-seq and RNA-seq datasets predicted TEAD2, TEAD4 and ZIC3 as potential factors involved in urine cell reprogramming. Transfection of TEAD2, TEAD4 and ZIC3 (in the presence of Yamanaka factors) significantly improved the reprogramming efficiency of urine cells. We confirmed that the newly generated iPSCs possessed pluripotency characteristics similar to normal H1 embryonic stem cells. We also confirmed that the new iPSCs could differentiate to functional cardiomyocytes.

CONCLUSIONS

In conclusion, TEAD2, TEAD4 and ZIC3 can increase the efficiency of reprogramming human urine cells into iPSCs, and provides a new stem cell sources for the clinical application and modeling of cardiovascular disease.

A plate-based single-cell ATAC-seq workflow for fast and robust profiling of chromatin accessibility

Profiling chromatin accessibility at the single-cell level provides critical information about cell type composition and cell-to-cell variation within a complex tissue. Emerging techniques for the interrogation of chromatin accessibility in individual cells allow investigation of the fundamental mechanisms that lead to the variability of different cells. This protocol describes a fast and robust method for single-cell chromatin accessibility profiling based on the assay for transposase-accessible chromatin using sequencing (ATAC-seq). The method combines up-front bulk Tn5 tagging of chromatin with flow cytometry to isolate single nuclei or cells. Reagents required to generate sequencing libraries are added to the same well in the plate where cells are sorted. The protocol described here generates data of high complexity and excellent signal-to-noise ratio and can be combined with index sorting for in-depth characterization of cell types. The whole experimental procedure can be finished within 1 or 2 d with a throughput of hundreds to thousands of nuclei, and the data can be processed by the provided computational pipeline. The execution of the protocol only requires basic techniques and equipment in a molecular biology laboratory with flow cytometry support.

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.

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.

Transcription Factor Reprogramming in the Inner Ear: Turning on Cell Fate Switches to Regenerate Sensory Hair Cells

Non-mammalian vertebrates can restore their auditory and vestibular hair cells naturally by triggering the regeneration of adjacent supporting cells. The transcription factor ATOH1 is a key regulator of hair cell development and regeneration in the inner ear. Following the death of hair cells, supporting cells upregulate ATOH1 and give rise to new hair cells. However, in the mature mammalian cochlea, such natural regeneration of hair cells is largely absent. Transcription factor reprogramming has been used in many tissues to convert one cell type into another, with the long-term hope of achieving tissue regeneration. Reprogramming transcription factors work by altering the transcriptomic and epigenetic landscapes in a target cell, resulting in a fate change to the desired cell type. Several studies have shown that ATOH1 is capable of reprogramming cochlear non-sensory tissue into cells resembling hair cells in young animals. However, the reprogramming ability of ATOH1 is lost with age, implying that the potency of individual hair cell-specific transcription factors may be reduced or lost over time by mechanisms that are still not clear. To circumvent this, combinations of key hair cell transcription factors have been used to promote hair cell regeneration in older animals. In this review, we summarize recent findings that have identified and studied these reprogramming factor combinations for hair cell regeneration. Finally, we discuss the important questions that emerge from these findings, particularly the feasibility of therapeutic strategies using reprogramming factors to restore human hearing in the future.

Machine Learning and Systems Biology Approaches to Characterize Dosage-Based Gene Dependencies in Cancer Cells

Mapping of cancer survivability factors allows for the identification of novel biological insights for drug targeting. Using genomic editing techniques, gene dependencies can be extracted in a high-throughput and quantitative manner. Dependencies have been predicted using machine learning techniques on -omics data, but the biological consequences of dependency predictor pairs has not been explored. In this work we devised a framework to explore gene dependency using an ensemble of machine learning methods, and our learned models captured meaningful biological information beyond just gene dependency prediction. We show that dosage-based dependent predictors (DDPs) primarily belonged to transcriptional regulation ontologies. We also found that anti-sense RNAs and long- noncoding RNA transcripts display DDPs. Network analyses revealed that SOX10, HLA-J, and ZEB2 act as a triad of network hubs in the dependent-predictor network. Collectively, we demonstrate the powerful combination of machine learning and systems biology approach can illuminate new insights in understanding gene dependency and guide novel targeting avenues.