Cell Reprogramming: The Many Roads to Success.

Cellular reprogramming experiments from somatic cell types have demonstrated the plasticity of terminally differentiated cell states. Recent efforts in understanding the mechanisms of cellular reprogramming have begun to elucidate the differentiation trajectories along the reprogramming processes. In this review, we focus mainly on direct reprogramming strategies by transcription factors and highlight the variables that contribute to cell fate conversion outcomes. We review key studies that shed light on the cellular and molecular mechanisms by investigating differentiation trajectories and alternative cell states as well as transcription factor regulatory activities during cell fate reprogramming. Finally, we highlight a few concepts that we believe require attention, particularly when measuring the success of cell reprogramming experiments.

An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming.

Alternative splicing (AS) is a key process underlying the expansion of proteomic diversity and the regulation of gene expression. Here, we identify an evolutionarily conserved embryonic stem cell (ESC)-specific AS event that changes the DNA-binding preference of the forkhead family transcription factor FOXP1. We show that the ESC-specific isoform of FOXP1 stimulates the expression of transcription factor genes required for pluripotency, including OCT4, NANOG, NR5A2, and GDF3, while concomitantly repressing genes required for ESC differentiation. This isoform also promotes the maintenance of ESC pluripotency and contributes to efficient reprogramming of somatic cells into induced pluripotent stem cells. These results reveal a pivotal role for an AS event in the regulation of pluripotency through the control of critical ESC-specific transcriptional programs.

CTCF-mediated topological boundaries during development foster appropriate gene regulation.

The genome is organized into repeating topologically associated domains (TADs), each of which is spatially isolated from its neighbor by poorly understood boundary elements thought to be conserved across cell types. Here, we show that deletion of CTCF (CCCTC-binding factor)-binding sites at TAD and sub-TAD topological boundaries that form within the HoxA and HoxC clusters during differentiation not only disturbs local chromatin domain organization and regulatory interactions but also results in homeotic transformations typical of Hox gene misregulation. Moreover, our data suggest that CTCF-dependent boundary function can be modulated by competing forces, such as the self-assembly of polycomb domains within the nucleus. Therefore, CTCF boundaries are not merely static structural components of the genome but instead are locally dynamic regulatory structures that control gene expression during development.

Differential abilities to engage inaccessible chromatin diversify vertebrate Hox binding patterns.

Although Hox genes encode for conserved transcription factors (TFs), they are further divided into anterior, central and posterior groups based on their DNA-binding domain similarity. The posterior Hox group expanded in the deuterostome clade and patterns caudal and distal structures. We aimed to address how similar Hox TFs diverge to induce different positional identities. We studied Hox TF DNA-binding and regulatory activity during an in vitro motor neuron differentiation system that recapitulates embryonic development. We found diversity in the genomic binding profiles of different Hox TFs, even among the posterior group paralogs that share similar DNA-binding domains. These differences in genomic binding were explained by differing abilities to bind to previously inaccessible sites. For example, the posterior group HOXC9 had a greater ability to bind occluded sites than the posterior HOXC10, producing different binding patterns and driving differential gene expression programs. From these results, we propose that the differential abilities of posterior Hox TFs to bind to previously inaccessible chromatin drive patterning diversification.This article has an associated 'The people behind the papers' interview.

Feedback from rhodopsin controls rhodopsin exclusion in Drosophila photoreceptors.

Sensory systems with high discriminatory power use neurons that express only one of several alternative sensory receptor proteins. This exclusive receptor gene expression restricts the sensitivity spectrum of neurons and is coordinated with the choice of their synaptic targets. However, little is known about how it is maintained throughout the life of a neuron. Here we show that the green-light sensing receptor rhodopsin 6 (Rh6) acts to exclude an alternative blue-sensitive rhodopsin 5 (Rh5) from a subset of Drosophila R8 photoreceptor neurons. Loss of Rh6 leads to a gradual expansion of Rh5 expression into all R8 photoreceptors of the ageing adult retina. The Rh6 feedback signal results in repression of the rh5 promoter and can be mimicked by other Drosophila rhodopsins; it is partly dependent on activation of rhodopsin by light, and relies on G(αq) activity, but not on the subsequent steps of the phototransduction cascade. Our observations reveal a thus far unappreciated spectral plasticity of R8 photoreceptors, and identify rhodopsin feedback as an exclusion mechanism.

An integrated model of multiple-condition ChIP-Seq data reveals predeterminants of Cdx2 binding.

Regulatory proteins can bind to different sets of genomic targets in various cell types or conditions. To reliably characterize such condition-specific regulatory binding we introduce MultiGPS, an integrated machine learning approach for the analysis of multiple related ChIP-seq experiments. MultiGPS is based on a generalized Expectation Maximization framework that shares information across multiple experiments for binding event discovery. We demonstrate that our framework enables the simultaneous modeling of sparse condition-specific binding changes, sequence dependence, and replicate-specific noise sources. MultiGPS encourages consistency in reported binding event locations across multiple-condition ChIP-seq datasets and provides accurate estimation of ChIP enrichment levels at each event. MultiGPS's multi-experiment modeling approach thus provides a reliable platform for detecting differential binding enrichment across experimental conditions. We demonstrate the advantages of MultiGPS with an analysis of Cdx2 binding in three distinct developmental contexts. By accurately characterizing condition-specific Cdx2 binding, MultiGPS enables novel insight into the mechanistic basis of Cdx2 site selectivity. Specifically, the condition-specific Cdx2 sites characterized by MultiGPS are highly associated with pre-existing genomic context, suggesting that such sites are pre-determined by cell-specific regulatory architecture. However, MultiGPS-defined condition-independent sites are not predicted by pre-existing regulatory signals, suggesting that Cdx2 can bind to a subset of locations regardless of genomic environment. A summary of this paper appears in the proceedings of the RECOMB 2014 conference, April 2-5.

Embryonic stem cell-based mapping of developmental transcriptional programs.

The study of developmentally regulated transcription factors by chromatin immunoprecipitation and deep sequencing (ChIP-seq) faces two major obstacles: availability of ChIP-grade antibodies and access to sufficient number of cells. We describe versatile genome-wide analysis of transcription-factor binding sites by combining directed differentiation of embryonic stem cells and inducible expression of tagged proteins. We demonstrate its utility by mapping DNA-binding sites of transcription factors involved in motor neuron specification.

Genome-wide screening identifies Trim33 as an essential regulator of dendritic cell differentiation.

The development of dendritic cells (DCs), including antigen-presenting conventional DCs (cDCs) and cytokine-producing plasmacytoid DCs (pDCs), is controlled by the growth factor Flt3 ligand (Flt3L) and its receptor Flt3. We genetically dissected Flt3L-driven DC differentiation using CRISPR-Cas9-based screening. Genome-wide screening identified multiple regulators of DC differentiation including subunits of TSC and GATOR1 complexes, which restricted progenitor growth but enabled DC differentiation by inhibiting mTOR signaling. An orthogonal screen identified the transcriptional repressor Trim33 (TIF-1γ) as a regulator of DC differentiation. Conditional targeting in vivo revealed an essential role of Trim33 in the development of all DCs, but not of monocytes or granulocytes. In particular, deletion of Trim33 caused rapid loss of DC progenitors, pDCs, and the cross-presenting cDC1 subset. Trim33-deficient Flt3+ progenitors up-regulated pro-inflammatory and macrophage-specific genes but failed to induce the DC differentiation program. Collectively, these data elucidate mechanisms that control Flt3L-driven differentiation of the entire DC lineage and identify Trim33 as its essential regulator.

Lineage-based identification of cellular states and expression programs.

We present a method, LineageProgram, that uses the developmental lineage relationship of observed gene expression measurements to improve the learning of developmentally relevant cellular states and expression programs. We find that incorporating lineage information allows us to significantly improve both the predictive power and interpretability of expression programs that are derived from expression measurements from in vitro differentiation experiments. The lineage tree of a differentiation experiment is a tree graph whose nodes describe all of the unique expression states in the input expression measurements, and edges describe the experimental perturbations applied to cells. Our method, LineageProgram, is based on a log-linear model with parameters that reflect changes along the lineage tree. Regularization with L(1) that based methods controls the parameters in three distinct ways: the number of genes change between two cellular states, the number of unique cellular states, and the number of underlying factors responsible for changes in cell state. The model is estimated with proximal operators to quickly discover a small number of key cell states and gene sets. Comparisons with existing factorization, techniques, such as singular value decomposition and non-negative matrix factorization show that our method provides higher predictive power in held, out tests while inducing sparse and biologically relevant gene sets.

Essential transcription factors for induced neuron differentiation.

Neurogenins are proneural transcription factors required to specify neuronal identity. Their overexpression in human pluripotent stem cells rapidly produces cortical-like neurons with spiking activity and, because of this, they have been widely adopted for human neuron disease models. However, we do not fully understand the key downstream regulatory effectors responsible for driving neural differentiation. Here, using inducible expression of NEUROG1 and NEUROG2, we identify transcription factors (TFs) required for directed neuronal differentiation by combining expression and chromatin accessibility analyses with a pooled in vitro CRISPR-Cas9 screen targeting all ~1900 TFs in the human genome. The loss of one of these essential TFs (ZBTB18) yields few MAP2-positive neurons. Differentiated ZBTB18-null cells have radically altered gene expression, leading to cytoskeletal defects and stunted neurites and spines. In addition to identifying key downstream TFs for neuronal differentiation, our work develops an integrative multi-omics and TFome-wide perturbation platform to rapidly characterize essential TFs for the differentiation of any human cell type.

ISSCR Education Committee syllabus and learning guide for enhancing stem cell literacy.

The rapidly evolving stem cell field puts much stress on developing educational resources. The ISSCR Education Committee has created a flexible stem cell syllabus rooted in core concepts to facilitate stem cell literacy. The free syllabus will be updated regularly to maintain accuracy and relevance.

Distinct visual pathways mediate Drosophila larval light avoidance and circadian clock entrainment.

Visual organs perceive environmental stimuli required for rapid initiation of behaviors and can also entrain the circadian clock. The larval eye of Drosophila is capable of both functions. Each eye contains only 12 photoreceptors (PRs), which can be subdivided into two subtypes. Four PRs express blue-sensitive rhodopsin5 (rh5) and eight express green-sensitive rhodopsin6 (rh6). We found that either PR-subtype is sufficient to entrain the molecular clock by light, while only the Rh5-PR subtype is essential for light avoidance. Acetylcholine released from PRs confers both functions. Both subtypes of larval PRs innervate the main circadian pacemaker neurons of the larva, the neuropeptide PDF (pigment-dispersing factor)-expressing lateral neurons (LNs), providing sensory input to control circadian rhythms. However, we show that PDF-expressing LNs are dispensable for light avoidance, and a distinct set of three clock neurons is required. Thus we have identified distinct sensory and central circuitry regulating light avoidance behavior and clock entrainment. Our findings provide insights into the coding of sensory information for distinct behavioral functions and the underlying molecular and neuronal circuitry.

Stochastic spineless expression creates the retinal mosaic for colour vision.

Drosophila colour vision is achieved by R7 and R8 photoreceptor cells present in every ommatidium. The fly retina contains two types of ommatidia, called 'pale' and 'yellow', defined by different rhodopsin pairs expressed in R7 and R8 cells. Similar to the human cone photoreceptors, these ommatidial subtypes are distributed stochastically in the retina. The choice between pale versus yellow ommatidia is made in R7 cells, which then impose their fate onto R8. Here we report that the Drosophila dioxin receptor Spineless is both necessary and sufficient for the formation of the ommatidial mosaic. A short burst of spineless expression at mid-pupation in a large subset of R7 cells precedes rhodopsin expression. In spineless mutants, all R7 and most R8 cells adopt the pale fate, whereas overexpression of spineless is sufficient to induce the yellow R7 fate. Therefore, this study suggests that the entire retinal mosaic required for colour vision is defined by the stochastic expression of a single transcription factor, Spineless.

Foxa2 and Pet1 Direct and Indirect Synergy Drive Serotonergic Neuronal Differentiation.

Neuronal programming by forced expression of transcription factors (TFs) holds promise for clinical applications of regenerative medicine. However, the mechanisms by which TFs coordinate their activities on the genome and control distinct neuronal fates remain obscure. Using direct neuronal programming of embryonic stem cells, we dissected the contribution of a series of TFs to specific neuronal regulatory programs. We deconstructed the Ascl1-Lmx1b-Foxa2-Pet1 TF combination that has been shown to generate serotonergic neurons and found that stepwise addition of TFs to Ascl1 canalizes the neuronal fate into a diffuse monoaminergic fate. The addition of pioneer factor Foxa2 represses Phox2b to induce serotonergic fate, similar to in vivo regulatory networks. Foxa2 and Pet1 appear to act synergistically to upregulate serotonergic fate. Foxa2 and Pet1 co-bind to a small fraction of genomic regions but mostly bind to different regulatory sites. In contrast to the combinatorial binding activities of other programming TFs, Pet1 does not strictly follow the Foxa2 pioneer. These findings highlight the challenges in formulating generalizable rules for describing the behavior of TF combinations that program distinct neuronal subtypes.

PRC1 sustains the integrity of neural fate in the absence of PRC2 function.

Polycomb repressive complexes (PRCs) 1 and 2 maintain stable cellular memories of early fate decisions by establishing heritable patterns of gene repression. PRCs repress transcription through histone modifications and chromatin compaction, but their roles in neuronal subtype diversification are poorly defined. We found that PRC1 is essential for the specification of segmentally restricted spinal motor neuron (MN) subtypes, while PRC2 activity is dispensable to maintain MN positional identities during terminal differentiation. Mutation of the core PRC1 component Ring1 in mice leads to increased chromatin accessibility and ectopic expression of a broad variety of fates determinants, including Hox transcription factors, while neuronal class-specific features are maintained. Loss of MN subtype identities in Ring1 mutants is due to the suppression of Hox-dependent specification programs by derepressed Hox13 paralogs (Hoxa13, Hoxb13, Hoxc13, Hoxd13). These results indicate that PRC1 can function in the absence of de novo PRC2-dependent histone methylation to maintain chromatin topology and postmitotic neuronal fate.

Capybara: A computational tool to measure cell identity and fate transitions.

Measuring cell identity in development, disease, and reprogramming is challenging as cell types and states are in continual transition. Here, we present Capybara, a computational tool to classify discrete cell identity and intermediate "hybrid" cell states, supporting a metric to quantify cell fate transition dynamics. We validate hybrid cells using experimental lineage tracing data to demonstrate the multi-lineage potential of these intermediate cell states. We apply Capybara to diagnose shortcomings in several cell engineering protocols, identifying hybrid states in cardiac reprogramming and off-target identities in motor neuron programming, which we alleviate by adding exogenous signaling factors. Further, we establish a putative in vivo correlate for induced endoderm progenitors. Together, these results showcase the utility of Capybara to dissect cell identity and fate transitions, prioritizing interventions to enhance the efficiency and fidelity of stem cell engineering.

CRISPR and biochemical screens identify MAZ as a cofactor in CTCF-mediated insulation at Hox clusters.

CCCTC-binding factor (CTCF) is critical to three-dimensional genome organization. Upon differentiation, CTCF insulates active and repressed genes within Hox gene clusters. We conducted a genome-wide CRISPR knockout (KO) screen to identify genes required for CTCF-boundary activity at the HoxA cluster, complemented by biochemical approaches. Among the candidates, we identified Myc-associated zinc-finger protein (MAZ) as a cofactor in CTCF insulation. MAZ colocalizes with CTCF at chromatin borders and, similar to CTCF, interacts with the cohesin subunit RAD21. MAZ KO disrupts gene expression and local contacts within topologically associating domains. Similar to CTCF motif deletions, MAZ motif deletions lead to derepression of posterior Hox genes immediately after CTCF boundaries upon differentiation, giving rise to homeotic transformations in mouse. Thus, MAZ is a factor contributing to appropriate insulation, gene expression and genomic architecture during development.

An expansion of the non-coding genome and its regulatory potential underlies vertebrate neuronal diversity.

Proper assembly and function of the nervous system requires the generation of a uniquely diverse population of neurons expressing a cell-type-specific combination of effector genes that collectively define neuronal morphology, connectivity, and function. How countless partially overlapping but cell-type-specific patterns of gene expression are controlled at the genomic level remains poorly understood. Here we show that neuronal genes are associated with highly complex gene regulatory systems composed of independent cell-type- and cell-stage-specific regulatory elements that reside in expanded non-coding genomic domains. Mapping enhancer-promoter interactions revealed that motor neuron enhancers are broadly distributed across the large chromatin domains. This distributed regulatory architecture is not a unique property of motor neurons but is employed throughout the nervous system. The number of regulatory elements increased dramatically during the transition from invertebrates to vertebrates, suggesting that acquisition of new enhancers might be a fundamental process underlying the evolutionary increase in cellular complexity.

The BTB transcription factors ZBTB11 and ZFP131 maintain pluripotency by repressing pro-differentiation genes.

In pluripotent cells, a delicate activation-repression balance maintains pro-differentiation genes ready for rapid activation. The identity of transcription factors (TFs) that specifically repress pro-differentiation genes remains obscure. By targeting ∼1,700 TFs with CRISPR loss-of-function screen, we found that ZBTB11 and ZFP131 are required for embryonic stem cell (ESC) pluripotency. ESCs without ZBTB11 or ZFP131 lose colony morphology, reduce proliferation rate, and upregulate transcription of genes associated with three germ layers. ZBTB11 and ZFP131 bind proximally to pro-differentiation genes. ZBTB11 or ZFP131 loss leads to an increase in H3K4me3, negative elongation factor (NELF) complex release, and concomitant transcription at associated genes. Together, our results suggest that ZBTB11 and ZFP131 maintain pluripotency by preventing premature expression of pro-differentiation genes and present a generalizable framework to maintain cellular potency.