Matching queried single-cell open-chromatin profiles to large pools of single-cell transcriptomes and epigenomes for reference supported analysis.

The true benefits of large single-cell transcriptome and epigenome data sets can be realized only with the development of new approaches and search tools for annotating individual cells. Matching a single-cell epigenome profile to a large pool of reference cells remains a major challenge. Here, we present scEpiSearch, which enables searching, comparison, and independent classification of single-cell open-chromatin profiles against a large reference of single-cell expression and open-chromatin data sets. Across performance benchmarks, scEpiSearch outperformed multiple methods in accuracy of search and low-dimensional coembedding of single-cell profiles, irrespective of platforms and species. Here we also demonstrate the unconventional utilities of scEpiSearch by applying it on single-cell epigenome profiles of K562 cells and samples from patients with acute leukaemia to reveal different aspects of their heterogeneity, multipotent behavior, and dedifferentiated states. Applying scEpiSearch on our single-cell open-chromatin profiles from embryonic stem cells (ESCs), we identified ESC subpopulations with more activity and poising for endoplasmic reticulum stress and unfolded protein response. Thus, scEpiSearch solves the nontrivial problem of amalgamating information from a large pool of single cells to identify and study the regulatory states of cells using their single-cell epigenomes.

TENET: gene network reconstruction using transfer entropy reveals key regulatory factors from single cell transcriptomic data.

Accurate prediction of gene regulatory rules is important towards understanding of cellular processes. Existing computational algorithms devised for bulk transcriptomics typically require a large number of time points to infer gene regulatory networks (GRNs), are applicable for a small number of genes and fail to detect potential causal relationships effectively. Here, we propose a novel approach 'TENET' to reconstruct GRNs from single cell RNA sequencing (scRNAseq) datasets. Employing transfer entropy (TE) to measure the amount of causal relationships between genes, TENET predicts large-scale gene regulatory cascades/relationships from scRNAseq data. TENET showed better performance than other GRN reconstructors, in identifying key regulators from public datasets. Specifically from scRNAseq, TENET identified key transcriptional factors in embryonic stem cells (ESCs) and during direct cardiomyocytes reprogramming, where other predictors failed. We further demonstrate that known target genes have significantly higher TE values, and TENET predicted higher TE genes were more influenced by the perturbation of their regulator. Using TENET, we identified and validated that Nme2 is a culture condition specific stem cell factor. These results indicate that TENET is uniquely capable of identifying key regulators from scRNAseq data.

Transcriptional Dynamics of Hepatic Sinusoid-Associated Cells After Liver Injury.

BACKGROUND AND AIMS: Hepatic sinusoidal cells are known actors in the fibrogenic response to injury. Activated hepatic stellate cells (HSCs), liver sinusoidal endothelial cells, and Kupffer cells are responsible for sinusoidal capillarization and perisinusoidal matrix deposition, impairing vascular exchange and heightening the risk of advanced fibrosis. While the overall pathogenesis is well understood, functional relations between cellular transitions during fibrogenesis are only beginning to be resolved. At single-cell resolution, we here explored the heterogeneity of individual cell types and dissected their transitions and crosstalk during fibrogenesis. APPROACH AND RESULTS: We applied single-cell transcriptomics to map the heterogeneity of sinusoid-associated cells in healthy and injured livers and reconstructed the single-lineage HSC trajectory from pericyte to myofibroblast. Stratifying each sinusoidal cell population by activation state, we projected shifts in sinusoidal communication upon injury. Weighted gene correlation network analysis of the HSC trajectory led to the identification of core genes whose expression proved highly predictive of advanced fibrosis in patients with nonalcoholic steatohepatitis (NASH). Among the core members of the injury-repressed gene module, we identified plasmalemma vesicle-associated protein (PLVAP) as a protein amply expressed by mouse and human HSCs. PLVAP expression was suppressed in activated HSCs upon injury and may hence define hitherto unknown roles for HSCs in the regulation of microcirculatory exchange and its breakdown in chronic liver disease. CONCLUSIONS: Our study offers a single-cell resolved account of drug-induced injury of the mammalian liver and identifies key genes that may serve important roles in sinusoidal integrity and as markers of advanced fibrosis in human NASH.

Power analysis of single-cell RNA-sequencing experiments.

Single-cell RNA sequencing (scRNA-seq) has become an established and powerful method to investigate transcriptomic cell-to-cell variation, thereby revealing new cell types and providing insights into developmental processes and transcriptional stochasticity. A key question is how the variety of available protocols compare in terms of their ability to detect and accurately quantify gene expression. Here, we assessed the protocol sensitivity and accuracy of many published data sets, on the basis of spike-in standards and uniform data processing. For our workflow, we developed a flexible tool for counting the number of unique molecular identifiers (https://github.com/vals/umis/). We compared 15 protocols computationally and 4 protocols experimentally for batch-matched cell populations, in addition to investigating the effects of spike-in molecular degradation. Our analysis provides an integrated framework for comparing scRNA-seq protocols.

SC3: consensus clustering of single-cell RNA-seq data.

Single-cell RNA-seq enables the quantitative characterization of cell types based on global transcriptome profiles. We present single-cell consensus clustering (SC3), a user-friendly tool for unsupervised clustering, which achieves high accuracy and robustness by combining multiple clustering solutions through a consensus approach (http://bioconductor.org/packages/SC3). We demonstrate that SC3 is capable of identifying subclones from the transcriptomes of neoplastic cells collected from patients.

BioModels: ten-year anniversary.

BioModels (http://www.ebi.ac.uk/biomodels/) is a repository of mathematical models of biological processes. A large set of models is curated to verify both correspondence to the biological process that the model seeks to represent, and reproducibility of the simulation results as described in the corresponding peer-reviewed publication. Many models submitted to the database are annotated, cross-referencing its components to external resources such as database records, and terms from controlled vocabularies and ontologies. BioModels comprises two main branches: one is composed of models derived from literature, while the second is generated through automated processes. BioModels currently hosts over 1200 models derived directly from the literature, as well as in excess of 140,000 models automatically generated from pathway resources. This represents an approximate 60-fold growth for literature-based model numbers alone, since BioModels' first release a decade ago. This article describes updates to the resource over this period, which include changes to the user interface, the annotation profiles of models in the curation pipeline, major infrastructure changes, ability to perform online simulations and the availability of model content in Linked Data form. We also outline planned improvements to cope with a diverse array of new challenges.

The systems biology format converter.

BACKGROUND: Interoperability between formats is a recurring problem in systems biology research. Many tools have been developed to convert computational models from one format to another. However, they have been developed independently, resulting in redundancy of efforts and lack of synergy. RESULTS: Here we present the System Biology Format Converter (SBFC), which provide a generic framework to potentially convert any format into another. The framework currently includes several converters translating between the following formats: SBML, BioPAX, SBGN-ML, Matlab, Octave, XPP, GPML, Dot, MDL and APM. This software is written in Java and can be used as a standalone executable or web service. CONCLUSIONS: The SBFC framework is an evolving software project. Existing converters can be used and improved, and new converters can be easily added, making SBFC useful to both modellers and developers. The source code and documentation of the framework are freely available from the project web site.

Computational assignment of cell-cycle stage from single-cell transcriptome data.

The transcriptome of single cells can reveal important information about cellular states and heterogeneity within populations of cells. Recently, single-cell RNA-sequencing has facilitated expression profiling of large numbers of single cells in parallel. To fully exploit these data, it is critical that suitable computational approaches are developed. One key challenge, especially pertinent when considering dividing populations of cells, is to understand the cell-cycle stage of each captured cell. Here we describe and compare five established supervised machine learning methods and a custom-built predictor for allocating cells to their cell-cycle stage on the basis of their transcriptome. In particular, we assess the impact of different normalisation strategies and the usage of prior knowledge on the predictive power of the classifiers. We tested the methods on previously published datasets and found that a PCA-based approach and the custom predictor performed best. Moreover, our analysis shows that the performance depends strongly on normalisation and the usage of prior knowledge. Only by leveraging prior knowledge in form of cell-cycle annotated genes and by preprocessing the data using a rank-based normalisation, is it possible to robustly capture the transcriptional cell-cycle signature across different cell types, organisms and experimental protocols.

Extensive folding variability between homologous chromosomes in mammalian cells.

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Oncogenic signaling pathways in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is the most common (∼90% cases) pancreatic neoplasm and one of the most lethal cancer among all malignances. PDAC harbor aberrant oncogenic signaling that may result from the multiple genetic and epigenetic alterations such as the mutation in driver genes (KRAS, CDKN2A, p53), genomic amplification of regulatory genes (MYC, IGF2BP2, ROIK3), deregulation of chromatin-modifying proteins (HDAC, WDR5) among others. A key event is the formation of Pancreatic Intraepithelial Neoplasia (PanIN) that often results from the activating mutation in KRAS. Mutated KRAS can direct a variety of signaling pathways and modulate downstream targets including MYC, which play an important role in cancer progression. In this review, we discuss recent literature shedding light on the origins of PDAC from the perspective of major oncogenic signaling pathways. We highlight how MYC directly and indirectly, with cooperation with KRAS, affect epigenetic reprogramming and metastasis. Additionally, we summarize the recent findings from single cell genomic approaches that highlight heterogeneity in PDAC and tumor microenvironment, and provide molecular avenues for PDAC treatment in the future.

A single-cell, time-resolved profiling of Xenopus mucociliary epithelium reveals nonhierarchical model of development.

The specialized cell types of the mucociliary epithelium (MCE) lining the respiratory tract enable continuous airway clearing, with its defects leading to chronic respiratory diseases. The molecular mechanisms driving cell fate acquisition and temporal specialization during mucociliary epithelial development remain largely unknown. Here, we profile the developing Xenopus MCE from pluripotent to mature stages by single-cell transcriptomics, identifying multipotent early epithelial progenitors that execute multilineage cues before specializing into late-stage ionocytes and goblet and basal cells. Combining in silico lineage inference, in situ hybridization, and single-cell multiplexed RNA imaging, we capture the initial bifurcation into early epithelial and multiciliated progenitors and chart cell type emergence and fate progression into specialized cell types. Comparative analysis of nine airway atlases reveals an evolutionary conserved transcriptional module in ciliated cells, whereas secretory and basal types execute distinct function-specific programs across vertebrates. We uncover a continuous nonhierarchical model of MCE development alongside a data resource for understanding respiratory biology.

CRISPR-Cas12a-integrated transgenes in genomic safe harbors retain high expression in human hematopoietic iPSC-derived lineages and primary cells.

Discovery of genomic safe harbor sites (SHSs) is fundamental for multiple transgene integrations, such as reporter genes, chimeric antigen receptors (CARs), and safety switches, which are required for safe cell products for regenerative cell therapies and immunotherapies. Here we identified and characterized potential SHS in human cells. Using the CRISPR-MAD7 system, we integrated transgenes at these sites in induced pluripotent stem cells (iPSCs), primary T and natural killer (NK) cells, and Jurkat cell line, and demonstrated efficient and stable expression at these loci. Subsequently, we validated the differentiation potential of engineered iPSC toward CD34+ hematopoietic stem and progenitor cells (HSPCs), lymphoid progenitor cells (LPCs), and NK cells and showed that transgene expression was perpetuated in these lineages. Finally, we demonstrated that engineered iPSC-derived NK cells retained expression of a non-virally integrated anti-CD19 CAR, suggesting that several of the investigated SHSs can be used to engineer cells for adoptive immunotherapies.

Exploration of cell state heterogeneity using single-cell proteomics through sensitivity-tailored data-independent acquisition.

Single-cell resolution analysis of complex biological tissues is fundamental to capture cell-state heterogeneity and distinct cellular signaling patterns that remain obscured with population-based techniques. The limited amount of material encapsulated in a single cell however, raises significant technical challenges to molecular profiling. Due to extensive optimization efforts, single-cell proteomics by Mass Spectrometry (scp-MS) has emerged as a powerful tool to facilitate proteome profiling from ultra-low amounts of input, although further development is needed to realize its full potential. To this end, we carry out comprehensive analysis of orbitrap-based data-independent acquisition (DIA) for limited material proteomics. Notably, we find a fundamental difference between optimal DIA methods for high- and low-load samples. We further improve our low-input DIA method by relying on high-resolution MS1 quantification, thus enhancing sensitivity by more efficiently utilizing available mass analyzer time. With our ultra-low input tailored DIA method, we are able to accommodate long injection times and high resolution, while keeping the scan cycle time low enough to ensure robust quantification. Finally, we demonstrate the capability of our approach by profiling mouse embryonic stem cell culture conditions, showcasing heterogeneity in global proteomes and highlighting distinct differences in key metabolic enzyme expression in distinct cell subclusters.

Cell Competition in Carcinogenesis.

The majority of human cancers evolve over time through the stepwise accumulation of somatic mutations followed by clonal selection akin to Darwinian evolution. However, the in-depth mechanisms that govern clonal dynamics and selection remain elusive, particularly during the earliest stages of tissue transformation. Cell competition (CC), often referred to as 'survival of the fittest' at the cellular level, results in the elimination of less fit cells by their more fit neighbors supporting optimal organism health and function. Alternatively, CC may allow an uncontrolled expansion of super-fit cancer cells to outcompete their less fit neighbors thereby fueling tumorigenesis. Recent research discussed herein highlights the various non-cell-autonomous principles, including interclonal competition and cancer microenvironment competition supporting the ability of a tumor to progress from the initial stages to tissue colonization. In addition, we extend current insights from CC-mediated clonal interactions and selection in normal tissues to better comprehend those factors that contribute to cancer development.

PPARγ lipodystrophy mutants reveal intermolecular interactions required for enhancer activation.

Peroxisome proliferator-activated receptor γ (PPARγ) is the master regulator of adipocyte differentiation, and mutations that interfere with its function cause lipodystrophy. PPARγ is a highly modular protein, and structural studies indicate that PPARγ domains engage in several intra- and inter-molecular interactions. How these interactions modulate PPARγ's ability to activate target genes in a cellular context is currently poorly understood. Here we take advantage of two previously uncharacterized lipodystrophy mutations, R212Q and E379K, that are predicted to interfere with the interaction of the hinge of PPARγ with DNA and with the interaction of PPARγ ligand binding domain (LBD) with the DNA-binding domain (DBD) of the retinoid X receptor, respectively. Using biochemical and genome-wide approaches we show that these mutations impair PPARγ function on an overlapping subset of target enhancers. The hinge region-DNA interaction appears mostly important for binding and remodelling of target enhancers in inaccessible chromatin, whereas the PPARγ-LBD:RXR-DBD interface stabilizes the PPARγ:RXR:DNA ternary complex. Our data demonstrate how in-depth analyses of lipodystrophy mutants can unravel molecular mechanisms of PPARγ function.

Predicting gene regulatory networks from cell atlases.

Recent single-cell RNA-sequencing atlases have surveyed and identified major cell types across different mouse tissues. Here, we computationally reconstruct gene regulatory networks from three major mouse cell atlases to capture functional regulators critical for cell identity, while accounting for a variety of technical differences, including sampled tissues, sequencing depth, and author assigned cell type labels. Extracting the regulatory crosstalk from mouse atlases, we identify and distinguish global regulons active in multiple cell types from specialised cell type-specific regulons. We demonstrate that regulon activities accurately distinguish individual cell types, despite differences between individual atlases. We generate an integrated network that further uncovers regulon modules with coordinated activities critical for cell types, and validate modules using available experimental data. Inferring regulatory networks during myeloid differentiation from wild-type and Irf8 KO cells, we uncover functional contribution of Irf8 regulon activity and composition towards monocyte lineage. Our analysis provides an avenue to further extract and integrate the regulatory crosstalk from single-cell expression data.

Publisher Correction: Single-cell RNA-sequencing of differentiating iPS cells reveals dynamic genetic effects on gene expression.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Single-cell transcriptomics identifies CD44 as a marker and regulator of endothelial to haematopoietic transition.

The endothelial to haematopoietic transition (EHT) is the process whereby haemogenic endothelium differentiates into haematopoietic stem and progenitor cells (HSPCs). The intermediary steps of this process are unclear, in particular the identity of endothelial cells that give rise to HSPCs is unknown. Using single-cell transcriptome analysis and antibody screening, we identify CD44 as a marker of EHT enabling us to isolate robustly the different stages of EHT in the aorta-gonad-mesonephros (AGM) region. This allows us to provide a detailed phenotypical and transcriptional profile of CD44-positive arterial endothelial cells from which HSPCs emerge. They are characterized with high expression of genes related to Notch signalling, TGFbeta/BMP antagonists, a downregulation of genes related to glycolysis and the TCA cycle, and a lower rate of cell cycle. Moreover, we demonstrate that by inhibiting the interaction between CD44 and its ligand hyaluronan, we can block EHT, identifying an additional regulator of HSPC development.

Single-cell RNA-sequencing of differentiating iPS cells reveals dynamic genetic effects on gene expression.

Recent developments in stem cell biology have enabled the study of cell fate decisions in early human development that are impossible to study in vivo. However, understanding how development varies across individuals and, in particular, the influence of common genetic variants during this process has not been characterised. Here, we exploit human iPS cell lines from 125 donors, a pooled experimental design, and single-cell RNA-sequencing to study population variation of endoderm differentiation. We identify molecular markers that are predictive of differentiation efficiency of individual lines, and utilise heterogeneity in the genetic background across individuals to map hundreds of expression quantitative trait loci that influence expression dynamically during differentiation and across cellular contexts.

Single-Cell Tagged Reverse Transcription (STRT-Seq).

Single-cell RNA sequencing (scRNA-seq) has become an established approach to profile entire transcriptomes of individual cells from different cell types, tissues, species, and organisms. Single-cell tagged reverse transcription sequencing (STRT-seq) is one of the early single-cell methods which utilize 5' tag counting of transcripts. STRT-seq performed on microfluidics Fluidigm C1 platform (STRT-C1) is a flexible scRNA-seq approach that allows for accurate, sensitive and importantly molecular counting of transcripts at single-cell level. Herein, I describe the STRT-C1 method and the steps involved in capturing 96 cells across C1 microfluidics chip, cDNA synthesis, and preparing single-cell libraries for Illumina short-read sequencing.