SCITO-seq: single-cell combinatorial indexed cytometry sequencing.

The development of DNA-barcoded antibodies to tag cell surface molecules has enabled the use of droplet-based single-cell sequencing (dsc-seq) to profile protein abundances from thousands of cells simultaneously. As compared to flow and mass cytometry, the high per cell cost of current dsc-seq-based workflows precludes their use in clinical applications and large-scale pooled screens. Here, we introduce SCITO-seq, a workflow that uses splint oligonucleotides (oligos) to enable combinatorially indexed dsc-seq of DNA-barcoded antibodies from over 105 cells per reaction using commercial microfluidics. By encoding sample barcodes into splint oligos, we demonstrate that multiplexed SCITO-seq produces reproducible estimates of cellular composition and surface protein expression comparable to those from mass cytometry. We further demonstrate two modified splint oligo designs that extend SCITO-seq to achieve compatibility with commercial DNA-barcoded antibodies and simultaneous expression profiling of the transcriptome and surface proteins from the same cell. These results demonstrate SCITO-seq as a flexible and ultra-high-throughput platform for sequencing-based single-cell protein and multimodal profiling.

Scalable, multimodal profiling of chromatin accessibility, gene expression and protein levels in single cells.

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary T cells.

Soluble α-synuclein-antibody complexes activate the NLRP3 inflammasome in hiPSC-derived microglia.

Parkinson's disease is characterized by accumulation of α-synuclein (αSyn). Release of oligomeric/fibrillar αSyn from damaged neurons may potentiate neuronal death in part via microglial activation. Heretofore, it remained unknown if oligomeric/fibrillar αSyn could activate the nucleotide-binding oligomerization domain (NOD)-like receptor (NLR) family pyrin domain-containing 3 (NLRP3) inflammasome in human microglia and whether anti-αSyn antibodies could prevent this effect. Here, we show that αSyn activates the NLRP3 inflammasome in human induced pluripotent stem cell (hiPSC)-derived microglia (hiMG) via dual stimulation involving Toll-like receptor 2 (TLR2) engagement and mitochondrial damage. In vitro, hiMG can be activated by mutant (A53T) αSyn secreted from hiPSC-derived A9-dopaminergic neurons. Surprisingly, αSyn-antibody complexes enhanced rather than suppressed inflammasome-mediated interleukin-1ß (IL-1ß) secretion, indicating these complexes are neuroinflammatory in a human context. A further increase in inflammation was observed with addition of oligomerized amyloid-ß peptide (Aß) and its cognate antibody. In vivo, engraftment of hiMG with αSyn in humanized mouse brain resulted in caspase-1 activation and neurotoxicity, which was exacerbated by αSyn antibody. These findings may have important implications for antibody therapies aimed at depleting misfolded/aggregated proteins from the human brain, as they may paradoxically trigger inflammation in human microglia.

Joint probabilistic modeling of single-cell multi-omic data with totalVI.

The paired measurement of RNA and surface proteins in single cells with cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) is a promising approach to connect transcriptional variation with cell phenotypes and functions. However, combining these paired views into a unified representation of cell state is made challenging by the unique technical characteristics of each measurement. Here we present Total Variational Inference (totalVI; https://scvi-tools.org ), a framework for end-to-end joint analysis of CITE-seq data that probabilistically represents the data as a composite of biological and technical factors, including protein background and batch effects. To evaluate totalVI's performance, we profiled immune cells from murine spleen and lymph nodes with CITE-seq, measuring over 100 surface proteins. We demonstrate that totalVI provides a cohesive solution for common analysis tasks such as dimensionality reduction, the integration of datasets with different measured proteins, estimation of correlations between molecules and differential expression testing.

Mouse models of neutropenia reveal progenitor-stage-specific defects.

Advances in genetics and sequencing have identified a plethora of disease-associated and disease-causing genetic alterations. To determine causality between genetics and disease, accurate models for molecular dissection are required; however, the rapid expansion of transcriptional populations identified through single-cell analyses presents a major challenge for accurate comparisons between mutant and wild-type cells. Here we generate mouse models of human severe congenital neutropenia (SCN) using patient-derived mutations in the GFI1 transcription factor. To determine the effects of SCN mutations, we generated single-cell references for granulopoietic genomic states with linked epitopes1, aligned mutant cells to their wild-type equivalents and identified differentially expressed genes and epigenetic loci. We find that GFI1-target genes are altered sequentially, as cells go through successive states of differentiation. These insights facilitated the genetic rescue of granulocytic specification but not post-commitment defects in innate immune effector function, and underscore the importance of evaluating the effects of mutations and therapy within each relevant cell state.

Single-Cell Virtual Cytometer allows user-friendly and versatile analysis and visualization of multimodal single cell RNAseq datasets.

The development of single-cell transcriptomic technologies yields large datasets comprising multimodal informations, such as transcriptomes and immunophenotypes. Despite the current explosion of methods for pre-processing and integrating multimodal single-cell data, there is currently no user-friendly software to display easily and simultaneously both immunophenotype and transcriptome-based UMAP/t-SNE plots from the pre-processed data. Here, we introduce Single-Cell Virtual Cytometer, an open-source software for flow cytometry-like visualization and exploration of pre-processed multi-omics single cell datasets. Using an original CITE-seq dataset of PBMC from an healthy donor, we illustrate its use for the integrated analysis of transcriptomes and epitopes of functional maturation in human peripheral T lymphocytes. So this free and open-source algorithm constitutes a unique resource for biologists seeking for a user-friendly analytic tool for multimodal single cell datasets.

Human Pluripotent Stem Cell-Derived Multipotent Vascular Progenitors of the Mesothelium Lineage Have Utility in Tissue Engineering and Repair.

In this report we describe a human pluripotent stem cell-derived vascular progenitor (MesoT) cell of the mesothelium lineage. MesoT cells are multipotent and generate smooth muscle cells, endothelial cells, and pericytes and self-assemble into vessel-like networks in vitro. MesoT cells transplanted into mechanically damaged neonatal mouse heart migrate into the injured tissue and contribute to nascent coronary vessels in the repair zone. When seeded onto decellularized vascular scaffolds, MesoT cells differentiate into the major vascular lineages and self-assemble into vasculature capable of supporting peripheral blood flow following transplantation. These findings demonstrate in vivo functionality and the potential utility of MesoT cells in vascular engineering applications.

Molecular analyses of neurogenic defects in a human pluripotent stem cell model of fragile X syndrome.

New research suggests that common pathways are altered in many neurodevelopmental disorders including autism spectrum disorder; however, little is known about early molecular events that contribute to the pathology of these diseases. The study of monogenic, neurodevelopmental disorders with a high incidence of autistic behaviours, such as fragile X syndrome, has the potential to identify genes and pathways that are dysregulated in autism spectrum disorder as well as fragile X syndrome. In vitro generation of human disease-relevant cell types provides the ability to investigate aspects of disease that are impossible to study in patients or animal models. Differentiation of human pluripotent stem cells recapitulates development of the neocortex, an area affected in both fragile X syndrome and autism spectrum disorder. We have generated induced human pluripotent stem cells from several individuals clinically diagnosed with fragile X syndrome and autism spectrum disorder. When differentiated to dorsal forebrain cell fates, our fragile X syndrome human pluripotent stem cell lines exhibited reproducible aberrant neurogenic phenotypes. Using global gene expression and DNA methylation profiling, we have analysed the early stages of neurogenesis in fragile X syndrome human pluripotent stem cells. We discovered aberrant DNA methylation patterns at specific genomic regions in fragile X syndrome cells, and identified dysregulated gene- and network-level correlates of fragile X syndrome that are associated with developmental signalling, cell migration, and neuronal maturation. Integration of our gene expression and epigenetic analysis identified altered epigenetic-mediated transcriptional regulation of a distinct set of genes in fragile X syndrome. These fragile X syndrome-aberrant networks are significantly enriched for genes associated with autism spectrum disorder, giving support to the idea that underlying similarities exist among these neurodevelopmental diseases.

Mission critical: the need for proteomics in the era of next-generation sequencing and precision medicine.

Next generation sequencing (NGS) has ignited an unprecedented pace of discovery in the biomedical sciences that is fundamentally transforming the way that we understand, diagnose and treat disease, and has motivated the belief that true precision medicine - medicine that is tailored to an individual's genetic, biochemical and exposure profile - will be a reality in the near term. With minimal sample requirement, NGS can enable the concurrent genome-wide study of genetic variations, transcriptomes, and certain epigenetic modifications. However, interrogating proteins as efficiently as DNA and RNA can be interrogated with NGS is lacking and this hampers more comprehensive views of molecular physiology and limits advances in biomedical science and precision medicine. The fact is that innovations in proteomic technologies pale in comparison to the advances in NGS, with current methodologies suffering from issues related to reproducibility, sensitivity, sample requirements, and limited multiplexing capacity. The development of proteomic technologies to overcome these limitations would fill the void in systems biology research, catalyze clinical innovations, and expedite the realization of precision medicine.

Whole-genome mutational burden analysis of three pluripotency induction methods.

There is concern that the stresses of inducing pluripotency may lead to deleterious DNA mutations in induced pluripotent stem cell (iPSC) lines, which would compromise their use for cell therapies. Here we report comparative genomic analysis of nine isogenic iPSC lines generated using three reprogramming methods: integrating retroviral vectors, non-integrating Sendai virus and synthetic mRNAs. We used whole-genome sequencing and de novo genome mapping to identify single-nucleotide variants, insertions and deletions, and structural variants. Our results show a moderate number of variants in the iPSCs that were not evident in the parental fibroblasts, which may result from reprogramming. There were only small differences in the total numbers and types of variants among different reprogramming methods. Most importantly, a thorough genomic analysis showed that the variants were generally benign. We conclude that the process of reprogramming is unlikely to introduce variants that would make the cells inappropriate for therapy.

Human stem cells from single blastomeres reveal pathways of embryonic or trophoblast fate specification.

Mechanisms of initial cell fate decisions differ among species. To gain insights into lineage allocation in humans, we derived ten human embryonic stem cell lines (designated UCSFB1-10) from single blastomeres of four 8-cell embryos and one 12-cell embryo from a single couple. Compared with numerous conventional lines from blastocysts, they had unique gene expression and DNA methylation patterns that were, in part, indicative of trophoblast competence. At a transcriptional level, UCSFB lines from different embryos were often more closely related than those from the same embryo. As predicted by the transcriptomic data, immunolocalization of EOMES, T brachyury, GDF15 and active ß-catenin revealed differential expression among blastomeres of 8- to 10-cell human embryos. The UCSFB lines formed derivatives of the three germ layers and CDX2-positive progeny, from which we derived the first human trophoblast stem cell line. Our data suggest heterogeneity among early-stage blastomeres and that the UCSFB lines have unique properties, indicative of a more immature state than conventional lines.

Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells.

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400-800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.

DNA methylation fingerprint of neuroblastoma reveals new biological and clinical insights.

AIM: To define the DNA methylation landscape of neuroblastoma and its clinicopathological impact. MATERIALS & METHODS: Microarray DNA methylation data were analyzed and associated with functional/regulatory genome annotation data, transcriptional profiles and clinicobiological parameters. RESULTS: DNA methylation changes in neuroblastoma affect not only promoters but also intragenic and intergenic regions at cytosine-phosphate-guanine (CpG) and non-CpG sites, and target functional chromatin domains of development and cancer-related genes such as CCND1. Tumors with diverse clinical risk showed differences affecting CpG and, remarkably, non-CpG sites. Non-CpG methylation observed essentially in clinically favorable cases was associated with the differentiation status of neuroblastoma and expression of key genes such as ALK. CONCLUSION: This epigenetic fingerprint of neuroblastoma provides new insights into the pathogenesis and clinical behavior of this pediatric tumor.

Increased risk of genetic and epigenetic instability in human embryonic stem cells associated with specific culture conditions.

The self-renewal and differentiation capacities of human pluripotent stem cells (hPSCs) make them a promising source of material for cell transplantation therapy, drug development, and studies of cellular differentiation and development. However, the large numbers of cells necessary for many of these applications require extensive expansion of hPSC cultures, a process that has been associated with genetic and epigenetic alterations. We have performed a combinatorial study on both hESCs and hiPSCs to compare the effects of enzymatic vs. mechanical passaging, and feeder-free vs. mouse embryonic fibroblast feeder substrate, on the genetic and epigenetic stability and the phenotypic characteristics of hPSCs. In extensive experiments involving over 100 continuous passages, we observed that both enzymatic passaging and feeder-free culture were associated with genetic instability, higher rates of cell proliferation, and persistence of OCT4/POU5F1-positive cells in teratomas, with enzymatic passaging having the stronger effect. In all combinations of culture conditions except for mechanical passaging on feeder layers, we noted recurrent deletions in the genomic region containing the tumor suppressor gene TP53, which was associated with decreased mRNA expression of TP53, as well as alterations in the expression of several downstream genes consistent with a decrease in the activity of the TP53 pathway. Among the hESC cultures, we also observed culture-associated variations in global gene expression and DNA methylation. The effects of enzymatic passaging and feeder-free conditions were also observed in hiPSC cultures. Our results highlight the need for careful assessment of the effects of culture conditions on cells intended for clinical therapies.

Application of a low cost array-based technique - TAB-Array - for quantifying and mapping both 5mC and 5hmC at single base resolution in human pluripotent stem cells.

5-hydroxymethylcytosine (5hmC), an oxidized derivative of 5-methylcytosine (5mC), has been implicated as an important epigenetic regulator of mammalian development. Current procedures use DNA sequencing methods to discriminate 5hmC from 5mC, limiting their accessibility to the scientific community. Here we report a method that combines TET-assisted bisulfite conversion with Illumina 450K DNA methylation arrays for a low-cost high-throughput approach that distinguishes 5hmC and 5mC signals at base resolution. Implementing this approach, termed "TAB-array", we assessed DNA methylation dynamics in the differentiation of human pluripotent stem cells into cardiovascular progenitors and neural precursor cells. With the ability to discriminate 5mC and 5hmC, we identified a large number of novel dynamically methylated genomic regions that are implicated in the development of these lineages. The increased resolution and accuracy afforded by this approach provides a powerful means to investigate the distinct contributions of 5mC and 5hmC in human development and disease.

Epigenetic therapy for Friedreich ataxia.

OBJECTIVE: To investigate whether a histone deacetylase inhibitor (HDACi) would be effective in an in vitro model for the neurodegenerative disease Friedreich ataxia (FRDA) and to evaluate safety and surrogate markers of efficacy in a phase I clinical trial in patients. METHODS: We used a human FRDA neuronal cell model, derived from patient induced pluripotent stem cells, to determine the efficacy of a 2-aminobenzamide HDACi (109) as a modulator of FXN gene expression and chromatin histone modifications. FRDA patients were dosed in 4 cohorts, ranging from 30mg/day to 240mg/day of the formulated drug product of HDACi 109, RG2833. Patients were monitored for adverse effects as well as for increases in FXN mRNA, frataxin protein, and chromatin modification in blood cells. RESULTS: In the neuronal cell model, HDACi 109/RG2833 increases FXN mRNA levels and frataxin protein, with concomitant changes in the epigenetic state of the gene. Chromatin signatures indicate that histone H3 lysine 9 is a key residue for gene silencing through methylation and reactivation through acetylation, mediated by the HDACi. Drug treatment in FRDA patients demonstrated increased FXN mRNA and H3 lysine 9 acetylation in peripheral blood mononuclear cells. No safety issues were encountered. INTERPRETATION: Drug exposure inducing epigenetic changes in neurons in vitro is comparable to the exposure required in patients to see epigenetic changes in circulating lymphoid cells and increases in gene expression. These findings provide a proof of concept for the development of an epigenetic therapy for this fatal neurological disease.

Epigenetic regulation of pluripotency and differentiation.

The precise, temporal order of gene expression during development is critical to ensure proper lineage commitment, cell fate determination, and ultimately, organogenesis. Epigenetic regulation of chromatin structure is fundamental to the activation or repression of genes during embryonic development. In recent years, there has been an explosion of research relating to various modes of epigenetic regulation, such as DNA methylation, post-translational histone tail modifications, noncoding RNA control of chromatin structure, and nucleosome remodeling. Technological advances in genome-wide epigenetic profiling and pluripotent stem cell differentiation have been primary drivers for elucidating the epigenetic control of cellular identity during development and nuclear reprogramming. Not only do epigenetic mechanisms regulate transcriptional states in a cell-type-specific manner but also they establish higher order genomic topology and nuclear architecture. Here, we review the epigenetic control of pluripotency and changes associated with pluripotent stem cell differentiation. We focus on DNA methylation, DNA demethylation, and common histone tail modifications. Finally, we briefly discuss epigenetic heterogeneity among pluripotent stem cell lines and the influence of epigenetic patterns on genome topology.

BMP4-directed trophoblast differentiation of human embryonic stem cells is mediated through a ΔNp63+ cytotrophoblast stem cell state.

The placenta is a transient organ that is necessary for proper fetal development. Its main functional component is the trophoblast, which is derived from extra-embryonic ectoderm. Little is known about early trophoblast differentiation in the human embryo, owing to lack of a proper in vitro model system. Human embryonic stem cells (hESCs) differentiate into functional trophoblast following BMP4 treatment in the presence of feeder-conditioned media; however, this model has not been widely accepted, in part owing to a lack of proof for a trophoblast progenitor population. We have previously shown that p63, a member of the p53 family of nuclear proteins, is expressed in proliferative cytotrophoblast (CTB), precursors to terminally differentiated syncytiotrophoblast (STB) in chorionic villi and extravillous trophoblast (EVT) at the implantation site. Here, we show that BMP4-treated hESCs differentiate into bona fide CTB by direct comparison with primary human placental tissues and isolated CTB through gene expression profiling. We show that, in primary CTB, p63 levels are reduced as cells differentiate into STB, and that forced expression of p63 maintains cyclin B1 and inhibits STB differentiation. We also establish that, similar to in vivo events, hESC differentiation into trophoblast is characterized by a p63(+)/KRT7(+) CTB stem cell state, followed by formation of functional KLF4(+) STB and HLA-G(+) EVT. Finally, we illustrate that downregulation of p63 by shRNA inhibits differentiation of hESCs into functional trophoblast. Taken together, our results establish that BMP4-treated hESCs are an excellent model of human trophoblast differentiation, closely mimicking the in vivo progression from p63(+) CTB stem cells to terminally differentiated trophoblast subtypes.

Matched miRNA and mRNA signatures from an hESC-based in vitro model of pancreatic differentiation reveal novel regulatory interactions.

The differentiation of human pluripotent stem cells (hPSCs) to insulin-expressing beta islet-like cells is a promising in vitro model system for studying the molecular signaling pathways underlying beta cell differentiation, as well as a potential source of cells for the treatment of type 1 diabetes. MicroRNAs (miRNAs) are a class of small non-coding RNAs that regulate many biological processes, including cellular differentiation. We studied the miRNA and mRNA expression profiles of hPSCs at five stages of in vitro differentiation along the pancreatic beta cell lineage (definitive endoderm, primitive gut tube, posterior foregut, pancreatic progenitor and hormone-expressing endocrine cells) in the context of samples of primary human fetal pancreas and purified adult islet cells using microarray analysis. Bioinformatic analysis of the resulting data identified a unique miRNA signature in differentiated beta islet cells, and predicted the effects of key miRNAs on mRNA expression. Many of the predicted miRNA-mRNA interactions involved mRNAs known to play key roles in the epithelial-mesenchymal transition process and pancreatic differentiation. We validated a subset of the predictions using qRT-PCR, luciferase reporter assays and western blotting, including the known interaction between miR-200 and ZEB2 (involved in epithelial-mesenchymal transition) and the novel interaction between miR-200 and SOX17 (a key transcription factor in specification of definitive endoderm). In addition, we found that miR-30d and let-7e, two miRNAs induced during differentiation, regulated the expression of RFX6, a transcription factor that directs pancreatic islet formation. These findings suggest that precise control of target mRNA expression by miRNAs ensures proper lineage specification during pancreatic development.