Chromatin profiling has emerged as a powerful means of genome annotation and detection of regulatory activity. The approach is especially well suited to the characterization of non-coding portions of the genome, which critically contribute to cellular phenotypes yet remain largely uncharted. Here we map nine chromatin marks across nine cell types to systematically characterize regulatory elements, their cell-type specificities and their functional interactions. Focusing on cell-type-specific patterns of promoters and enhancers, we define multicell activity profiles for chromatin state, gene expression, regulatory motif enrichment and regulator expression. We use correlations between these profiles to link enhancers to putative target genes, and predict the cell-type-specific activators and repressors that modulate them. The resulting annotations and regulatory predictions have implications for the interpretation of genome-wide association studies. Top-scoring disease single nucleotide polymorphisms are frequently positioned within enhancer elements specifically active in relevant cell types, and in some cases affect a motif instance for a predicted regulator, thus suggesting a mechanism for the association. Our study presents a general framework for deciphering cis-regulatory connections and their roles in disease.

Despite rapid progress in characterizing transcription factor-driven reprogramming of somatic cells to an induced pluripotent stem cell (iPSC) state, many mechanistic questions still remain. To gain insight into the earliest events in the reprogramming process, we systematically analyzed the transcriptional and epigenetic changes that occur during early factor induction after discrete numbers of divisions. We observed rapid, genome-wide changes in the euchromatic histone modification, H3K4me2, at more than a thousand loci including large subsets of pluripotency-related or developmentally regulated gene promoters and enhancers. In contrast, patterns of the repressive H3K27me3 modification remained largely unchanged except for focused depletion specifically at positions where H3K4 methylation is gained. These chromatin regulatory events precede transcriptional changes within the corresponding loci. Our data provide evidence for an early, organized, and population-wide epigenetic response to ectopic reprogramming factors that clarify the temporal order through which somatic identity is reset during reprogramming.

In embryonic stem (ES) cells, bivalent chromatin domains with overlapping repressive (H3 lysine 27 tri-methylation) and activating (H3 lysine 4 tri-methylation) histone modifications mark the promoters of more than 2,000 genes. To gain insight into the structure and function of bivalent domains, we mapped key histone modifications and subunits of Polycomb-repressive complexes 1 and 2 (PRC1 and PRC2) genomewide in human and mouse ES cells by chromatin immunoprecipitation, followed by ultra high-throughput sequencing. We find that bivalent domains can be segregated into two classes -- the first occupied by both PRC2 and PRC1 (PRC1-positive) and the second specifically bound by PRC2 (PRC2-only). PRC1-positive bivalent domains appear functionally distinct as they more efficiently retain lysine 27 tri-methylation upon differentiation, show stringent conservation of chromatin state, and associate with an overwhelming number of developmental regulator gene promoters. We also used computational genomics to search for sequence determinants of Polycomb binding. This analysis revealed that the genomewide locations of PRC2 and PRC1 can be largely predicted from the locations, sizes, and underlying motif contents of CpG islands. We propose that large CpG islands depleted of activating motifs confer epigenetic memory by recruiting the full repertoire of Polycomb complexes in pluripotent cells.

Gene-editing nucleases enable targeted modification of DNA sequences in living cells, thereby facilitating efficient knockout and precise editing of endogenous loci. Engineered nucleases also have the potential to introduce mutations at off-target sites of action. Such unintended alterations can confound interpretation of experiments and can have implications for development of therapeutic applications. Recently, two improved methods for identifying the off-target effects of zinc finger nucleases (ZFNs) were described-one using an in vitro cleavage site selection method and the other exploiting the insertion of integration-defective lentiviruses into nuclease-induced double-stranded DNA breaks. However, application of these two methods to a ZFN pair targeted to the human CCR5 gene led to identification of largely non-overlapping off-target sites, raising the possibility that additional off-target sites might exist. Here, we show that in silico abstraction of ZFN cleavage profiles obtained from in vitro cleavage site selections can greatly enhance the ability to identify potential off-target sites in human cells. Our improved method should enable more comprehensive profiling of ZFN specificities.

Autoimmune diabetes is a complex multifactorial disease with genetic and environmental factors playing pivotal roles. While many genes associated with the risk of diabetes have been identified to date, the mechanisms by which external triggers contribute to the genetic predisposition remain unclear. Here, we derived embryonic stem (ES) cell lines from diabetes-prone non-obese diabetic (NOD) and healthy C57BL/6 (B6) mice. While overall pluripotency markers were indistinguishable between newly derived NOD and B6 ES cells, we discovered several differentially expressed genes that normally are not expressed in ES cells. Several genes that reside in previously identified insulin-dependent diabetics (Idd) genomic regions were up-regulated in NOD ES cells. Gene set enrichment analysis showed that different groups of genes associated with immune functions are differentially expressed in NOD. Transcriptomic analysis of NOD blastocysts validated several differentially overexpressed Idd genes compared to B6. Genome-wide mapping of active histone modifications using ChIP-Seq supports active expression as the promoters and enhancers of activated genes are also marked by active histone modifications. We have also found that NOD ES cells secrete more inflammatory cytokines. Our data suggest that the known genetic predisposition of NOD to autoimmune diabetes leads to epigenetic instability of several Idd regions.

Direct conversion of human somatic fibroblasts into induced neurons (iNs) allows for the generation of functional neurons while bypassing any stem cell intermediary stages. Although iN technology has an enormous potential for modeling age-related diseases, as well as therapeutic approaches, the technology faces limitations due to variable conversion efficiencies and a lack of thorough understanding of the signaling pathways directing iN conversion. Here, we introduce a new all-in-one inducible lentiviral system that simplifies fibroblast transgenesis for the two pioneer transcription factors, Ngn2 and Ascl1, and markedly improves iN yields. Further, our timeline RNA-Seq data across the course of conversion has identified signaling pathways that become transcriptionally enriched during iN conversion. Small molecular modulators were identified for four signaling pathways that reliably increase the yield of iNs. Taken together, these advances provide an improved toolkit for iN technology and new insight into the mechanisms influencing direct iN conversion.

While B cells play a significant role in the onset of type-1 diabetes (T1D), little is know about their role in those early stages. Thus, to gain new insights into the role of B cells in T1D, we converted a physiological early pancreas-infiltrating B cell into a novel BCR mouse model using Somatic Cell Nuclear Transfer (SCNT). Strikingly, SCNT-derived B1411 model displayed neither developmental block nor anergy. Instead, B1411 underwent spontaneous germinal center reactions. Without T cell help, B1411-Rag1-/- was capable of forming peri-/intra-pancreatic lymph nodes, and undergoing class-switching. RNA-Seq analysis identified 93 differentially expressed genes in B1411 compared to WT B cells, including Irf7, Usp18, and Mda5 that had been linked to a potential viral etiology of T1D. We also found various members of the oligoadenylate synthase (OAS) family to be enriched in B1411, such as Oas1, which had recently also been linked to T1D. Strikingly, when challenged with glucose B1411-Rag1-/- mice displayed impaired glucose tolerance.

Sporadic Alzheimer's disease (AD) exclusively affects elderly people. Using direct conversion of AD patient fibroblasts into induced neurons (iNs), we generated an age-equivalent neuronal model. AD patient-derived iNs exhibit strong neuronal transcriptome signatures characterized by downregulation of mature neuronal properties and upregulation of immature and progenitor-like signaling pathways. Mapping iNs to longitudinal neuronal differentiation trajectory data demonstrated that AD iNs reflect a hypo-mature neuronal identity characterized by markers of stress, cell cycle, and de-differentiation. Epigenetic landscape profiling revealed an underlying aberrant neuronal state that shares similarities with malignant transformation and age-dependent epigenetic erosion. To probe for the involvement of aging, we generated rejuvenated iPSC-derived neurons that showed no significant disease-related transcriptome signatures, a feature that is consistent with epigenetic clock and brain ontogenesis mapping, which indicate that fibroblast-derived iNs more closely reflect old adult brain stages. Our findings identify AD-related neuronal changes as age-dependent cellular programs that impair neuronal identity.

We have previously reported an antigen-specific protocol to induce transplant tolerance and linked suppression to human embryonic stem cell (hESC)-derived tissues in immunocompetent mice through coreceptor and costimulation blockade. However, the exact mechanisms of acquired immune tolerance in this model have remained unclear.

The diverse functions of WASP, the deficiency of which causes Wiskott-Aldrich syndrome (WAS), remain poorly defined. We generated three isogenic WAS models using patient induced pluripotent stem cells and genome editing. These models recapitulated WAS phenotypes and revealed that WASP deficiency causes an upregulation of numerous RNA splicing factors and widespread altered splicing. Loss of WASP binding to splicing factor gene promoters frequently leads to aberrant epigenetic activation. WASP interacts with dozens of nuclear speckle constituents and constrains SRSF2 mobility. Using an optogenetic system, we showed that WASP forms phase-separated condensates that encompasses SRSF2, nascent RNA and active Pol II. The role of WASP in gene body condensates is corroborated by ChIPseq and RIPseq. Together our data reveal that WASP is a nexus regulator of RNA splicing that controls the transcription of splicing factors epigenetically and the dynamics of the splicing machinery through liquid-liquid phase separation.

Single-cell transcriptomics of neocortical neurons have revealed more than 100 clusters corresponding to putative cell types. For inhibitory and subcortical projection neurons (SCPNs), there is a strong concordance between clusters and anatomical descriptions of cell types. In contrast, cortico-cortical projection neurons (CCPNs) separate into surprisingly few transcriptomic clusters, despite their diverse anatomical projection types. We used projection-dependent single-cell transcriptomic analyses and monosynaptic rabies tracing to compare mouse primary visual cortex CCPNs projecting to different higher visual areas. We find that layer 2/3 CCPNs with different anatomical projections differ systematically in their gene expressions, despite forming only a single genetic cluster. Furthermore, these neurons receive feedback selectively from the same areas to which they project. These findings demonstrate that gene-expression analysis in isolation is insufficient to identify neuron types and have important implications for understanding the functional role of cortical feedback circuits.

Aging is a major risk factor for many human diseases, and in vitro generation of human neurons is an attractive approach for modeling aging-related brain disorders. However, modeling aging in differentiated human neurons has proved challenging. We generated neurons from human donors across a broad range of ages, either by iPSC-based reprogramming and differentiation or by direct conversion into induced neurons (iNs). While iPSCs and derived neurons did not retain aging-associated gene signatures, iNs displayed age-specific transcriptional profiles and revealed age-associated decreases in the nuclear transport receptor RanBP17. We detected an age-dependent loss of nucleocytoplasmic compartmentalization (NCC) in donor fibroblasts and corresponding iNs and found that reduced RanBP17 impaired NCC in young cells, while iPSC rejuvenation restored NCC in aged cells. These results show that iNs retain important aging-related signatures, thus allowing modeling of the aging process in vitro, and they identify impaired NCC as an important factor in human aging.

Polycomb proteins are epigenetic regulators that localize to developmental loci in the early embryo where they mediate lineage-specific gene repression. In Drosophila, these repressors are recruited to sequence elements by DNA binding proteins associated with Polycomb repressive complex 2 (PRC2). However, the sequences that recruit PRC2 in mammalian cells have remained obscure. To address this, we integrated a series of engineered bacterial artificial chromosomes into embryonic stem (ES) cells and examined their chromatin. We found that a 44 kb region corresponding to the Zfpm2 locus initiates de novo recruitment of PRC2. We then pinpointed a CpG island within this locus as both necessary and sufficient for PRC2 recruitment. Based on this causal demonstration and prior genomic analyses, we hypothesized that large GC-rich elements depleted of activating transcription factor motifs mediate PRC2 recruitment in mammals. We validated this model in two ways. First, we showed that a constitutively active CpG island is able to recruit PRC2 after excision of a cluster of activating motifs. Second, we showed that two 1 kb sequence intervals from the Escherichia coli genome with GC-contents comparable to a mammalian CpG island are both capable of recruiting PRC2 when integrated into the ES cell genome. Our findings demonstrate a causal role for GC-rich sequences in PRC2 recruitment and implicate a specific subset of CpG islands depleted of activating motifs as instrumental for the initial localization of this key regulator in mammalian genomes.

Despite widespread interest in using human induced pluripotent stem cells (hiPSCs) in neurological disease modeling, a suitable model system to study human neuronal connectivity is lacking. Here, we report a comprehensive and efficient differentiation paradigm for hiPSCs that generate multiple CA3 pyramidal neuron subtypes as detected by single-cell RNA sequencing (RNA-seq). This differentiation paradigm exhibits characteristics of neuronal network maturation, and rabies virus tracing revealed synaptic connections between stem cell-derived dentate gyrus (DG) and CA3 neurons in vitro recapitulating the neuronal connectivity within the hippocampus. Because hippocampal dysfunction has been implicated in schizophrenia, we applied DG and CA3 differentiation paradigms to schizophrenia-patient-derived hiPSCs. We detected reduced activity in DG-CA3 co-culture and deficits in spontaneous and evoked activity in CA3 neurons from schizophrenia-patient-derived hiPSCs. Our approach offers critical insights into the network activity aspects of schizophrenia and may serve as a promising tool for modeling diseases with hippocampal vulnerability. VIDEO ABSTRACT.

To date, it remains unclear if there are specific cell-surface markers for purifying glucose-responsive pancreatic β-like cells derived from human pluripotent stem cells (hPSCs). In searching for this, we generated an efficient protocol for differentiating β-like cells from human embryonic stem cells. We performed single-cell RNA sequencing and found that CD9 is a negative cell-surface marker of β-like cells, as most INS+ cells are CD9-. We purified β-like cells for spontaneous formation of islet-like organoids against CD9, and found significantly more NKX6.1+MAFA+C-PEPTIDE+ β-like cells in the CD9- than in the CD9+ population. CD9- cells also demonstrate better glucose responsiveness than CD9+ cells. In humans, we observe more CD9+C-PEPTIDE+ β cells in the fetal than in the adult cadaveric islets and more Ki67+ proliferating cells among CD9+ fetal β cells. Taken together, our experiments show that CD9 is a cell-surface marker for negative enrichment of glucose-responsive β-like cells differentiated from hPSCs.

The histone variant H2A.Z has been implicated in nucleosome exchange, transcriptional activation and Polycomb repression. However, the relationships among these seemingly disparate functions remain obscure.

Astrocyte dysfunction and neuroinflammation are detrimental features in multiple pathologies of the CNS. Therefore, the development of methods that produce functional human astrocytes represents an advance in the study of neurological diseases. Here we report an efficient method for inflammation-responsive astrocyte generation from induced pluripotent stem cells (iPSCs) and embryonic stem cells. This protocol uses an intermediate glial progenitor stage and generates functional astrocytes that show levels of glutamate uptake and calcium activation comparable with those observed in human primary astrocytes. Stimulation of stem cell-derived astrocytes with interleukin-1β or tumor necrosis factor α elicits a strong and rapid pro-inflammatory response. RNA-sequencing transcriptome profiling confirmed that similar gene expression changes occurred in iPSC-derived and primary astrocytes upon stimulation with interleukin-1β. This protocol represents an important tool for modeling in-a-dish neurological diseases with an inflammatory component, allowing for the investigation of the role of diseased astrocytes in neuronal degeneration.

No abstract available

Rhabdomyosarcomas (RMS) are phenotypically and functionally heterogeneous. Both primary human RMS cultures and low-passage Myf6Cre,Pax3:Foxo1,p53 mouse RMS cell lines, which express the fusion oncoprotein Pax3:Foxo1 and lack the tumor suppressor Tp53 (Myf6Cre,Pax3:Foxo1,p53), exhibit marked heterogeneity in PAX3:FOXO1 (P3F) expression at the single cell level. In mouse RMS cells, P3F expression is directed by the Pax3 promoter and coupled to eYFP YFPlow/P3Flow mouse RMS cells included 87% G0/G1 cells and reorganized their actin cytoskeleton to produce a cellular phenotype characterized by more efficient adhesion and migration. This translated into higher tumor-propagating cell frequencies of YFPlow/P3Flow compared with YFPhigh/P3Fhigh cells. Both YFPlow/P3Flow and YFPhigh/P3Fhigh cells gave rise to mixed clones in vitro, consistent with fluctuations in P3F expression over time. Exposure to the anti-tropomyosin compound TR100 disrupted the cytoskeleton and reversed enhanced migration and adhesion of YFPlow/P3Flow RMS cells. Heterogeneous expression of PAX3:FOXO1 at the single cell level may provide a critical advantage during tumor progression.

Glioblastomas are malignant tumors of the central nervous system hallmarked by subclonal diversity and dynamic adaptation amid developmental hierarchies. The source of dynamic reorganization within the spatial context of these tumors remains elusive. Here, we characterized glioblastomas by spatially resolved transcriptomics, metabolomics, and proteomics. By deciphering regionally shared transcriptional programs across patients, we infer that glioblastoma is organized by spatial segregation of lineage states and adapts to inflammatory and/or metabolic stimuli, reminiscent of the reactive transformation in mature astrocytes. Integration of metabolic imaging and imaging mass cytometry uncovered locoregional tumor-host interdependence, resulting in spatially exclusive adaptive transcriptional programs. Inferring copy-number alterations emphasizes a spatially cohesive organization of subclones associated with reactive transcriptional programs, confirming that environmental stress gives rise to selection pressure. A model of glioblastoma stem cells implanted into human and rodent neocortical tissue mimicking various environments confirmed that transcriptional states originate from dynamic adaptation to various environments.