Optimal-Transport Analysis of Single-Cell Gene Expression Identifies Developmental Trajectories in Reprogramming

Employing core regulatory circuits to define cell identity

The interplay between extrinsic signaling and downstream gene networks controls the establishment of cell identity during development and its maintenance in adult life. Advances in next-generation sequencing and single-cell technologies have revealed additional layers of complexity in cell identity. Here, we review our current understanding of transcription factor (TF) networks as key determinants of cell identity. We discuss the concept of the core regulatory circuit as a set of TFs and interacting factors that together define the gene expression profile of the cell. We propose the core regulatory circuit as a comprehensive conceptual framework for defining cellular identity and discuss its connections to cell function in different contexts.

What are the applications of single-cell RNA sequencing in cancer research: a systematic review

Single-cell RNA sequencing (scRNA-seq) is a tool for studying gene expression at the single-cell level that has been widely used due to its unprecedented high resolution. In the present review, we outline the preparation process and sequencing platforms for the scRNA-seq analysis of solid tumor specimens and discuss the main steps and methods used during data analysis, including quality control, batch-effect correction, normalization, cell cycle phase assignment, clustering, cell trajectory and pseudo-time reconstruction, differential expression analysis and gene set enrichment analysis, as well as gene regulatory network inference. Traditional bulk RNA sequencing does not address the heterogeneity within and between tumors, and since the development of the first scRNA-seq technique, this approach has been widely used in cancer research to better understand cancer cell biology and pathogenetic mechanisms. ScRNA-seq has been of great significance for the development of targeted therapy and immunotherapy. In the second part of this review, we focus on the application of scRNA-seq in solid tumors, and summarize the findings and achievements in tumor research afforded by its use. ScRNA-seq holds promise for improving our understanding of the molecular characteristics of cancer, and potentially contributing to improved diagnosis, prognosis, and therapeutics.

Dynamics-based data science in biology

With the increasingly accumulated bio-data, dynamics-based data-science has been progressing as an efficient way to reveal mechanisms of dynamical biological processes. We review three applications on detecting the tipping-points of diseases, quantifying cell's potency, and predicting time-series, to show the importance of dynamics-based data-science.

A cellular and spatial map of the choroid plexus across brain ventricles and ages

The choroid plexus (ChP) in each brain ventricle produces cerebrospinal fluid (CSF) and forms the blood-CSF barrier. Here, we construct a single-cell and spatial atlas of each ChP in the developing, adult, and aged mouse brain. We delineate diverse cell types, subtypes, cell states, and expression programs in epithelial and mesenchymal cells across ages and ventricles. In the developing ChP, we predict a common progenitor pool for epithelial and neuronal cells, validated by lineage tracing. Epithelial and fibroblast cells show regionalized expression by ventricle, starting at embryonic stages and persisting with age, with a dramatic transcriptional shift with maturation, and a smaller shift in each aged cell type. With aging, epithelial cells upregulate host-defense programs, and resident macrophages upregulate interleukin-1β (IL-1β) signaling genes. Our atlas reveals cellular diversity, architecture and signaling across ventricles during development, maturation, and aging of the ChP-brain barrier.

A single-embryo, single-cell time-resolved model for mouse gastrulation

Mouse embryonic development is a canonical model system for studying mammalian cell fate acquisition. Recently, single-cell atlases comprehensively charted embryonic transcriptional landscapes, yet inference of the coordinated dynamics of cells over such atlases remains challenging. Here, we introduce a temporal model for mouse gastrulation, consisting of data from 153 individually sampled embryos spanning 36 h of molecular diversification. Using algorithms and precise timing, we infer differentiation flows and lineage specification dynamics over the embryonic transcriptional manifold. Rapid transcriptional bifurcations characterize the commitment of early specialized node and blood cells. However, for most lineages, we observe combinatorial multi-furcation dynamics rather than hierarchical transcriptional transitions. In the mesoderm, dozens of transcription factors combinatorially regulate multifurcations, as we exemplify using time-matched chimeric embryos of Foxc1/Foxc2 mutants. Our study rejects the notion of differentiation being governed by a series of binary choices, providing an alternative quantitative model for cell fate acquisition.

The Cellular and Physiological Basis for Lung Repair and Regeneration: Past, Present, and Future

The respiratory system, which includes the trachea, airways, and distal alveoli, is a complex multi-cellular organ that intimately links with the cardiovascular system to accomplish gas exchange. In this review and as members of the NIH/NHLBI-supported Progenitor Cell Translational Consortium, we discuss key aspects of lung repair and regeneration. We focus on the cellular compositions within functional niches, cell-cell signaling in homeostatic health, the responses to injury, and new methods to study lung repair and regeneration. We also provide future directions for an improved understanding of the cell biology of the respiratory system, as well as new therapeutic avenues.

Statistical mechanics meets single-cell biology

Single-cell omics is transforming our understanding of cell biology and disease, yet the systems-level analysis and interpretation of single-cell data faces many challenges. In this Perspective, we describe the impact that fundamental concepts from statistical mechanics, notably entropy, stochastic processes and critical phenomena, are having on single-cell data analysis. We further advocate the need for more bottom-up modelling of single-cell data and to embrace a statistical mechanics analysis paradigm to help attain a deeper understanding of single-cell systems biology.

Diverse cellular players orchestrate regeneration after wounding

Fibrosis is one of the largest sources of human morbidity. The skin is a complex organ where interplay between diverse cell types and signalling pathways is essential both in homeostasis and wound repair, which can result in fibrosis or regeneration. This makes skin a useful model to study fibrosis and regeneration. While fibrosis often occurs postinjury, both clinical and laboratory observations suggest skin regeneration, complete with reconstituted cell diversity and de novo hair follicles, is possible. Extensive research performed in pursuit of skin regeneration has elucidated the key players, both cellular and molecular. Interestingly, some cells known for their homeostatic function are not implicated in regeneration or wound-induced hair neogenesis (WIHN), suggesting regeneration harnesses separate functional pathways from embryogenesis or other non-homeostatic mechanisms. For example, classic bulge cells, noted for their role in normally cycling hair follicles, do not finally contribute to long-lived cells in the regenerated tissue. During healing, multiple populations of cells, among them specific epithelial lineages, mesenchymal cells, and immune cells promote regenerative outcomes in the wounded skin. Ultimately, targeting specific populations of cells will be essential in manipulating a postwound environment to favour regeneration in lieu of fibrosis.

Joint profiling of DNA and proteins in single cells to dissect genotype-phenotype associations in leukemia

Studies of acute myeloid leukemia rely on DNA sequencing and immunophenotyping by flow cytometry as primary tools for disease characterization. However, leukemia tumor heterogeneity complicates integration of DNA variants and immunophenotypes from separate measurements. Here we introduce DAb-seq, a technology for simultaneous capture of DNA genotype and cell surface phenotype from single cells at high throughput, enabling direct profiling of proteogenomic states in tens of thousands of cells. To demonstrate the approach, we analyze the disease of three patients with leukemia over multiple treatment timepoints and disease recurrences. We observe complex genotype-phenotype dynamics that illustrate the subtlety of the disease process and the degree of incongruity between blast cell genotype and phenotype in different clinical scenarios. Our results highlight the importance of combined single-cell DNA and protein measurements to fully characterize the heterogeneity of leukemia.

Revealing lineage-related signals in single-cell gene expression using random matrix theory

Gene expression profiles of a cellular population, generated by single-cell RNA sequencing, contains rich information about biological state, including cell type, cell cycle phase, gene regulatory patterns, and location within the tissue of origin. A major challenge is to disentangle information about these different biological states from each other, including distinguishing from cell lineage, since the correlation of cellular expression patterns is necessarily contaminated by ancestry. Here, we use a recent advance in random matrix theory, discovered in the context of protein phylogeny, to identify differentiation or ancestry-related processes in single-cell data. Qin and Colwell [C. Qin, L. J. Colwell, Proc. Natl. Acad. Sci. U.S.A. 115, 690-695 (2018)] showed that ancestral relationships in protein sequences create a power-law signature in the covariance eigenvalue distribution. We demonstrate the existence of such signatures in scRNA-seq data and that the genes driving them are indeed related to differentiation and developmental pathways. We predict the existence of similar power-law signatures for cells along linear trajectories and demonstrate this for linearly differentiating systems. Furthermore, we generalize to show that the same signatures can arise for cells along tissue-specific spatial trajectories. We illustrate these principles in diverse tissues and organisms, including the mammalian epidermis and lung, Drosophila whole-embryo, adult Hydra, dendritic cells, the intestinal epithelium, and cells undergoing induced pluripotent stem cells (iPSC) reprogramming. We show how these results can be used to interpret the gradual dynamics of lineage structure along iPSC reprogramming. Together, we provide a framework that can be used to identify signatures of specific biological processes in single-cell data without prior knowledge and identify candidate genes associated with these processes.

Identifying transposable element expression dynamics and heterogeneity during development at the single-cell level with a processing pipeline scTE

Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity in unclear ways. Single-cell sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE expression has not been addressed. Here, we develop a single-cell TE processing pipeline, scTE, and report the expression of TEs in single cells in a range of biological contexts. Specific TE types are expressed in subpopulations of embryonic stem cells and are dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs are expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we apply scTE to single-cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results classify the dynamic patterns of TEs in single cells and their contributions to cell heterogeneity.

How transposable elements (TE) contribute to cell fate changes is unclear. Here, the authors generate a pipeline to quantify TE expression from single cell data. They show the dynamic expression of TEs from gastrulation to somatic cell reprogramming and human disease

Taming human brain organoids one cell at a time

Human brain organoids are self-organizing three-dimensional structures that emerge from human pluripotent stem cells and mimic aspects of the cellular composition and functionality of the developing human brain. Despite their impressive self-organizing capacity, organoids lack the stereotypic structural anatomy of their in vivo counterpart, making conventional analysis techniques underpowered to assess cellular composition and gene network regulation in organoids. Advances in single cell transcriptomics have recently allowed characterization and improvement of organoid protocols, as they continue to evolve, by enabling identification of cell types and states along with their developmental origins. In this review, we summarize recent approaches, progresses and challenges in resolving brain organoid's complexity through single-cell transcriptomics. We then discuss emerging technologies that may complement single-cell RNA sequencing by providing additional readouts of cellular states to generate an organ-level view of developmental processes. Altogether, these integrative technologies will allow monitoring of global gene regulation in thousands of individual cells and will offer an unprecedented opportunity to investigate features of human brain development and disease across multiple cellular modalities and with cell-type resolution.

Physics approach to the variable-mass optimal-transport problem

Optimal transport (OT) has become a discipline by itself that offers solutions to a wide range of theoretical problems in probability and mathematics with applications in several applied fields such as imaging sciences, machine learning, and in data sciences in general. The traditional OT problem suffers from a severe limitation: its balance condition imposes that the two distributions to be compared be normalized and have the same total mass. However, it is important for many applications to be able to relax this constraint and allow for mass creation and/or destruction. This is true, for example, in all problems requiring partial matching. In this paper, we propose an approach to solving a generalized version of the OT problem, which we refer to as the discrete variable-mass optimal-transport (VMOT) problem, using techniques adapted from statistical physics. Our first contribution is to fully describe this formalism, including all the proofs of its main claims. In particular, we derive a strongly concave effective free-energy function that captures the constraints of the VMOT problem at a finite temperature. From its maximum we derive a weak distance (i.e., a divergence) between possibly unbalanced distribution functions. The temperature-dependent OT distance decreases monotonically to the standard variable-mass OT distance, providing a robust framework for temperature annealing. Our second contribution is to show that the implementation of this formalism has the same properties as the regularized OT algorithms in time complexity, making it a competitive approach to solving the VMOT problem. We illustrate applications of the framework to the problem of partial two- and three-dimensional shape-matching problems.

A Casz1-NuRD complex regulates temporal identity transitions in neural progenitors

Neural progenitor cells undergo identity transitions during development to ensure the generation different types of neurons and glia in the correct sequence and proportions. A number of temporal identity factors that control these transitions in progenitor competence have been identified, but the molecular mechanisms underlying their function remain unclear. Here, we asked how Casz1, the mammalian orthologue of Drosophila castor, regulates competence during retinal development. We show that Casz1 is required to control the transition between neurogenesis and gliogenesis. Using BioID proteomics, we reveal that Casz1 interacts with the nucleosome remodeling and deacetylase (NuRD) complex in retinal cells. Finally, we show that both the NuRD and the polycomb repressor complexes are required for Casz1 to promote the rod fate and suppress gliogenesis. As additional temporal identity factors have been found to interact with the NuRD complex in other contexts, we propose that these factors might act through this common biochemical process to regulate neurogenesis.

Deciphering cell-cell interactions and communication from gene expression

Cell-cell interactions orchestrate organismal development, homeostasis and single-cell functions. When cells do not properly interact or improperly decode molecular messages, disease ensues. Thus, the identification and quantification of intercellular signalling pathways has become a common analysis performed across diverse disciplines. The expansion of protein-protein interaction databases and recent advances in RNA sequencing technologies have enabled routine analyses of intercellular signalling from gene expression measurements of bulk and single-cell data sets. In particular, ligand-receptor pairs can be used to infer intercellular communication from the coordinated expression of their cognate genes. In this Review, we highlight discoveries enabled by analyses of cell-cell interactions from transcriptomic data and review the methods and tools used in this context.

What machine learning can do for developmental biology

Developmental biology has grown into a data intensive science with the development of high-throughput imaging and multi-omics approaches. Machine learning is a versatile set of techniques that can help make sense of these large datasets with minimal human intervention, through tasks such as image segmentation, super-resolution microscopy and cell clustering. In this Spotlight, I introduce the key concepts, advantages and limitations of machine learning, and discuss how these methods are being applied to problems in developmental biology. Specifically, I focus on how machine learning is improving microscopy and single-cell 'omics' techniques and data analysis. Finally, I provide an outlook for the futures of these fields and suggest ways to foster new interdisciplinary developments.

Computational Stem Cell Biology: Open Questions and Guiding Principles

Computational biology is enabling an explosive growth in our understanding of stem cells and our ability to use them for disease modeling, regenerative medicine, and drug discovery. We discuss four topics that exemplify applications of computation to stem cell biology: cell typing, lineage tracing, trajectory inference, and regulatory networks. We use these examples to articulate principles that have guided computational biology broadly and call for renewed attention to these principles as computation becomes increasingly important in stem cell biology. We also discuss important challenges for this field with the hope that it will inspire more to join this exciting area.

Dppa2/4 as a trigger of signaling pathways to promote zygote genome activation by binding to CG-rich region

Developmental pluripotency-associated 2 (Dppa2) and developmental pluripotency-associated 4 (Dppa4) as positive drivers were helpful for transcriptional regulation of zygotic genome activation (ZGA). Here, we systematically assessed the cooperative interplay of Dppa2 and Dppa4 in regulating cell pluripotency and found that simultaneous overexpression of Dppa2/4 can make induced pluripotent stem cells closer to embryonic stem cells (ESCs). Compared with other pluripotency transcription factors, Dppa2/4 can regulate majorities of signaling pathways by binding on CG-rich region of proximal promoter (0-500 bp), of which 85% and 77% signaling pathways were significantly activated by Dppa2 and Dppa4, respectively. Notably, Dppa2/4 also can dramatically trigger the decisive signaling pathways for facilitating ZGA, including Hippo, MAPK and TGF-beta signaling pathways and so on. At last, we found alkaline phosphatase, placental-like 2 (Alppl2) was completely silenced when Dppa2 and 4 single- or double-knockout in ESC, which is consistent with Dux. Moreover, Alppl2 was significantly activated in mouse 2-cell embryos and 4-8 cells stage of human embryos, further predicted that Alppl2 was directly regulated by Dppa2/4 as a ZGA candidate driver to facilitate pre-embryonic development.

D-EE: Distributed software for visualizing intrinsic structure of large-scale single-cell data

BACKGROUND

Dimensionality reduction and visualization play vital roles in single-cell RNA sequencing (scRNA-seq) data analysis. While they have been extensively studied, state-of-the-art dimensionality reduction algorithms are often unable to preserve the global structures underlying data. Elastic embedding (EE), a nonlinear dimensionality reduction method, has shown promise in revealing low-dimensional intrinsic local and global data structure. However, the current implementation of the EE algorithm lacks scalability to large-scale scRNA-seq data.

RESULTS

We present a distributed optimization implementation of the EE algorithm, termed distributed elastic embedding (D-EE). D-EE reveals the low-dimensional intrinsic structures of data with accuracy equal to that of elastic embedding, and it is scalable to large-scale scRNA-seq data. It leverages distributed storage and distributed computation, achieving memory efficiency and high-performance computing simultaneously. In addition, an extended version of D-EE, termed distributed optimization implementation of time-series elastic embedding (D-TSEE), enables the user to visualize large-scale time-series scRNA-seq data by incorporating experimentally temporal information. Results with large-scale scRNA-seq data indicate that D-TSEE can uncover oscillatory gene expression patterns by using experimentally temporal information.

CONCLUSIONS

D-EE is a distributed dimensionality reduction and visualization tool. Its distributed storage and distributed computation technique allow us to efficiently analyze large-scale single-cell data at the cost of constant time speedup. The source code for D-EE algorithm based on C and MPI tailored to a high-performance computing cluster is available at https://github.com/ShaokunAn/D-EE.

A Synthesis Concerning Conservation and Divergence of Cell Types across Epithelia

Advances in single-cell RNA-seq (scRNA-seq) and computational analysis have enabled the systematic interrogation of the cellular composition of tissues. Combined with tools from developmental biology, cell biology, and genetics, these approaches are revealing fundamental aspects of tissue geometry and physiology, including the distribution, origins, and inferred functions of specialized cell types, and the dynamics of cellular turnover and differentiation. By comparing different tissues, such studies can delineate shared and specialized features of cell types and their lineage. Here, we compare two developmentally related murine epithelia, the airway and the small intestinal epithelia, which are both derived from the embryonic endodermal gut tube. We examine how airway and intestine generate and functionalize common archetypal cell types to fulfill similar shared physiologic functionalities. We point to cases in which similar cell types are repurposed to accommodate each tissue's unique physiologic role, and highlight tissue-specific cells whose specializations contribute to the distinct functional roles of each organ. We discuss how archetypal and unique cell types are incorporated within a cellular lineage, and how the regulation of the proportions of these cell types enables tissue-level organization to meet functional demands and maintain homeostasis.

Single-cell transcriptomics uncover distinct innate and adaptive cell subsets during tissue homeostasis and regeneration

Recently, immune cell-mediated tissue repair and regeneration has been an emerging paradigm of regenerative medicine. Immune cells form an essential part of the wound as induction of inflammation is a necessary step to elicit tissue healing. Rapid progress in transcriptomic analyses by high-throughput next-generation sequencing has been developed to study gene regulatory network and establish molecular signatures of immune cells that could potentially predict their functional roles in tissue repair and regeneration. However, the identification of cellular heterogeneity especially on the rare cell subsets has been limited in transcriptomic analyses of bulk cell populations. Therefore, genome-wide, single-cell RNA sequencing (scRNA-Seq) has offered an unprecedented approach to unravel cellular diversity and to study novel immune cell populations involved in tissue repair and regeneration through unsupervised sampling of individual cells without the need to rely on prior knowledge about cell-specific markers. The analysis of gene expression patterns at a single-cell resolution also holds promises to uncover the mechanisms and therefore the development of therapeutic strategy promoting immunoregenerative medicine. In this review, we will discuss how scRNA-Seq facilitates the characterization of immune cells, including macrophages, innate lymphoid cells and T and B lymphocytes, discovery of immune cell heterogeneity, identification of novel subsets, and tracking of developmental trajectories of distinct immune cells during tissue homeostasis, repair, and regeneration.

Batch-Corrected Distance Mitigates Temporal and Spatial Variability for Clustering and Visualization of Single-Cell Gene Expression Data

Clustering and visualization are essential parts of single-cell gene expression data analysis. The Euclidean distance used in most distance-based methods is not optimal. Batch effect, i.e., the variability among samples gathered from different times, tissues, and patients, introduces large between-group distance and obscures the true identities of cells. To solve this problem, we introduce Batch-Corrected Distance (BCD), a metric using temporal/spatial locality of the batch effect to control for such factors. We validate BCD on a simulated data as well as applied it to a mouse retina development dataset and a lung dataset. We also found the utility of our approach in understanding the progression of the Coronavirus Disease 2019 (COVID-19). BCD achieves more accurate clusters and better visualizations than state-of-the-art batch correction methods on longitudinal datasets. BCD can be directly integrated with most clustering and visualization methods to enable more scientific findings.

JMJD3 acts in tandem with KLF4 to facilitate reprogramming to pluripotency

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3’s catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription.

Previous work suggested that histone demethylase JMJD3 is detrimental to somatic cell reprogramming. Here, the authors show that while JMJD3 has a context-independent detrimental effect on early stages of reprogramming, during late stages it activates epithelial and pluripotency genes together with Klf4.

The Ground Zero of Organismal Life and Aging

Cells may naturally proceed or be forced to transition to a state with a radically lower biological age, that is, be rejuvenated. Examples are the conversion of somatic cells to induced pluripotent stem cells and rejuvenation of the germline with each generation. We posit that these processes converge to the same 'ground zero', the mid-embryonic state characterized by the lowest biological age where both organismal life and aging begin. It may also be related to the phylotypic state. The ground zero model clarifies the relationship between aging, development, rejuvenation, and de-differentiation, which are distinct throughout life. By extending the rejuvenation phase during early embryogenesis and editing the genome, it may be possible to achieve the biological age at the ground zero lower than that achieved naturally.

Methodologies for Following EMT In Vivo at Single Cell Resolution

An epithelial-mesenchymal transition (EMT) occurs in almost every metazoan embryo at the time mesoderm begins to differentiate. Several embryos have a long record as models for studying an EMT given that a known population of cells enters the EMT at a known time thereby enabling a detailed study of the process. Often, however, it is difficult to learn the molecular details of these model EMT systems because the transitioning cells are a minority of the population of cells in the embryo and in most cases there is an inability to isolate that population. Here we provide a method that enables an examination of genes expressed before, during, and after the EMT with a focus on just the cells that undergo the transition. Single cell RNA-seq (scRNA-seq) has advanced as a technology making it feasible to study the trajectory of gene expression specifically in the cells of interest, in vivo, and without the background noise of other cell populations. The sea urchin skeletogenic cells constitute only 5% of the total number of cells in the embryo yet with scRNA-seq it is possible to study the genes expressed by these cells without background noise. This approach, though not perfect, adds a new tool for uncovering the mechanism of EMT in this cell type.

Single-cell multiomics: technologies and data analysis methods

Advances in single-cell isolation and barcoding technologies offer unprecedented opportunities to profile DNA, mRNA, and proteins at a single-cell resolution. Recently, bulk multiomics analyses, such as multidimensional genomic and proteogenomic analyses, have proven beneficial for obtaining a comprehensive understanding of cellular events. This benefit has facilitated the development of single-cell multiomics analysis, which enables cell type-specific gene regulation to be examined. The cardinal features of single-cell multiomics analysis include (1) technologies for single-cell isolation, barcoding, and sequencing to measure multiple types of molecules from individual cells and (2) the integrative analysis of molecules to characterize cell types and their functions regarding pathophysiological processes based on molecular signatures. Here, we summarize the technologies for single-cell multiomics analyses (mRNA-genome, mRNA-DNA methylation, mRNA-chromatin accessibility, and mRNA-protein) as well as the methods for the integrative analysis of single-cell multiomics data.

The expansion of single-cell profiling technologies will provide unprecedented insights into the molecular mechanisms inherent in disease. Novel technologies known collectively as ‘single-cell multiomics’ enable systematic, high-resolution profiling of DNA, RNA and proteins in individual cells. This provides valuable data about gene regulation and molecular populations, and cellular processes during disease development and progression. Daehee Hwang and co-workers at Seoul National University, Seoul, South Korea, reviewed existing single-cell multiomics technologies and highlighted ways to integrate the data generated. Analytical features of multiomics allow scientists to isolate, sequence and label (or ‘barcode’) multiple molecules in single cells. Different sequencing techniques can be used for different purposes, such as exploring gene mutation coverage or measuring RNA transcripts. Combining these sequencing data will help identify links between significant features during disease.

Tempora: Cell trajectory inference using time-series single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) can map cell types, states and transitions during dynamic biological processes such as tissue development and regeneration. Many trajectory inference methods have been developed to order cells by their progression through a dynamic process. However, when time series data is available, most of these methods do not consider the available time information when ordering cells and are instead designed to work only on a single scRNA-seq data snapshot. We present Tempora, a novel cell trajectory inference method that orders cells using time information from time-series scRNA-seq data. In performance comparison tests, Tempora inferred known developmental lineages from three diverse tissue development time series data sets, beating state of the art methods in accuracy and speed. Tempora works at the level of cell clusters (types) and uses biological pathway information to help identify cell type relationships. This approach increases gene expression signal from single cells, processing speed, and interpretability of the inferred trajectory. Our results demonstrate the utility of a combination of time and pathway information to supervise trajectory inference for scRNA-seq based analysis.

Single-cell RNA sequencing (scRNA-seq) enables an unparalleled ability to map the heterogeneity of dynamic multicellular processes, such as tissue development, tumor growth, wound response and repair, and inflammation. Multiple methods have been developed to order cells along a pseudotime axis that represents a trajectory through such processes using the concept that cells that are closely related in a lineage will have similar transcriptomes. However, time series experiments provide another useful information source to order cells, from earlier to later time point. By introducing a novel use of biological pathway prior information, our Tempora algorithm improves the accuracy and speed of cell trajectory inference from time-series scRNA-seq data as measured by reconstructing known developmental trajectories from three diverse data sets. By analyzing scRNA-seq data at the cluster (cell type) level instead of at the single-cell level and by using known pathway information, Tempora amplifies gene expression signals from one cell using similar cells in a cluster and similar genes within a pathway. This approach also reduces computational time and resources needed to analyze large data sets because it works with a relatively small number of clusters instead of a potentially large number of cells. Finally, it eases interpretation, via operating on a relatively small number of clusters which usually represent known cell types, as well as by identifying time-dependent pathways. Tempora is useful for finding novel insights in dynamic processes.

scIGANs: single-cell RNA-seq imputation using generative adversarial networks

Single-cell RNA-sequencing (scRNA-seq) enables the characterization of transcriptomic profiles at the single-cell resolution with increasingly high throughput. However, it suffers from many sources of technical noises, including insufficient mRNA molecules that lead to excess false zero values, termed dropouts. Computational approaches have been proposed to recover the biologically meaningful expression by borrowing information from similar cells in the observed dataset. However, these methods suffer from oversmoothing and removal of natural cell-to-cell stochasticity in gene expression. Here, we propose the generative adversarial networks (GANs) for scRNA-seq imputation (scIGANs), which uses generated cells rather than observed cells to avoid these limitations and balances the performance between major and rare cell populations. Evaluations based on a variety of simulated and real scRNA-seq datasets show that scIGANs is effective for dropout imputation and enhances various downstream analysis. ScIGANs is robust to small datasets that have very few genes with low expression and/or cell-to-cell variance. ScIGANs works equally well on datasets from different scRNA-seq protocols and is scalable to datasets with over 100 000 cells. We demonstrated in many ways with compelling evidence that scIGANs is not only an application of GANs in omics data but also represents a competing imputation method for the scRNA-seq data.

Characterization of DNA Methylation Patterns and Mining of Epigenetic Markers During Genomic Reprogramming in SCNT Embryos

Somatic cell nuclear transfer (SCNT), also known as somatic cell cloning, is a commonly used technique to study epigenetic reprogramming. Although SCNT has the advantages of being safe and able to obtain pluripotent cells, early developmental arrest happens in most SCNT embryos. Overcoming epigenetic barriers is currently the primary strategy for improving reprogramming efficiency and improving developmental rate in SCNT embryos. In this study, we analyzed DNA methylation profiles of in vivo fertilized embryos and SCNT embryos with different developmental fates. Overall DNA methylation level was higher in SCNT embryos during global de-methylation process compared to in vivo fertilized embryos. In addition, promoter region, first intron and 3′UTR were found to be the major genomic regions that were hyper-methylated in SCNT embryos. Surprisingly, we found the length of re-methylated region was directly related to the change of methylation level. Furthermore, a number of genes including Dppa2 and Dppa4 which are important for early zygotic genome activation (ZGA) were not properly activated in SCNT embryos. This study comprehensively analyzed genome-wide DNA methylation patterns in SCNT embryos and provided candidate target genes for improving efficiency of genomic reprogramming in SCNT embryos. Since SCNT technology has been widely used in agricultural and pastoral production, protection of endangered animals, and therapeutic cloning, the findings of this study have significant importance for all these fields.

Conformational state switching and pathways of chromosome dynamics in cell cycle

The cell cycle is a process and function of a cell with different phases essential for cell growth, proliferation, and replication. It depends on the structure and dynamics of the underlying DNA molecule, which underpins the genome function. A microscopic structural-level understanding of how a genome or its functional module chromosome performs the cell cycle in terms of large-scale conformational transformation between different phases, such as the interphase and the mitotic phase, is still challenging. Here, we develop a non-equilibrium, excitation-relaxation energy landscape-switching model to quantify the underlying chromosome conformational transitions through (de-)condensation for a complete microscopic understanding of the cell cycle. We show that the chromosome conformational transition mechanism from the interphase to the mitotic phase follows a two-stage scenario, in good agreement with the experiments. In contrast, the mitotic exit pathways show the existence of an over-expanded chromosome that recapitulates the chromosome in the experimentally identified intermediate state at the telophase. We find the conformational pathways are heterogeneous and irreversible as a result of the non-equilibrium dynamics of the cell cycle from both structural and kinetic perspectives. We suggest that the irreversibility is mainly due to the distinct participation of the ATP-dependent structural maintenance of chromosomal protein complexes during the cell cycle. Our findings provide crucial insights into the microscopic molecular structural and dynamical physical mechanism for the cell cycle beyond the previous more macroscopic descriptions. Our non-equilibrium landscape framework is general and applicable to study diverse non-equilibrium physical and biological processes such as active matter, differentiation/development, and cancer.

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

Cellular reprogramming suffers from low efficiency especially for the human cells. To deconstruct the heterogeneity and unravel the mechanisms for successful reprogramming, we adopted single-cell RNA sequencing (scRNA-Seq) and single-cell assay for transposase-accessible chromatin (scATAC-Seq) to profile reprogramming cells across various time points. Our analysis revealed that reprogramming cells proceed in an asynchronous trajectory and diversify into heterogeneous subpopulations. We identified fluorescent probes and surface markers to enrich for the early reprogrammed human cells. Furthermore, combinatory usage of the surface markers enabled the fine segregation of the early-intermediate cells with diverse reprogramming propensities. scATAC-Seq analysis further uncovered the genomic partitions and transcription factors responsible for the regulatory phasing of reprogramming process. Binary choice between a FOSL1 and a TEAD4-centric regulatory network determines the outcome of a successful reprogramming. Together, our study illuminates the multitude of diverse routes transversed by individual reprogramming cells and presents an integrative roadmap for identifying the mechanistic part list of the reprogramming machinery.

Tuning parameters of dimensionality reduction methods for single-cell RNA-seq analysis

BACKGROUND

Many computational methods have been developed recently to analyze single-cell RNA-seq (scRNA-seq) data. Several benchmark studies have compared these methods on their ability for dimensionality reduction, clustering, or differential analysis, often relying on default parameters. Yet, given the biological diversity of scRNA-seq datasets, parameter tuning might be essential for the optimal usage of methods, and determining how to tune parameters remains an unmet need.

RESULTS

Here, we propose a benchmark to assess the performance of five methods, systematically varying their tunable parameters, for dimension reduction of scRNA-seq data, a common first step to many downstream applications such as cell type identification or trajectory inference. We run a total of 1.5 million experiments to assess the influence of parameter changes on the performance of each method, and propose two strategies to automatically tune parameters for methods that need it.

CONCLUSIONS

We find that principal component analysis (PCA)-based methods like scran and Seurat are competitive with default parameters but do not benefit much from parameter tuning, while more complex models like ZinbWave, DCA, and scVI can reach better performance but after parameter tuning.

Defining Human Pluripotency

Human pluripotent stem cells harbor the capacity to differentiate into cells from the three embryonic germ layers, and this ability grants them a central role in modeling human disorders and in the field of regenerative medicine. Here, we review pluripotency in human cells with respect to four different aspects: (1) embryonic development, (2) transcriptomes of pluripotent cell stages, (3) genes and pathways that reprogram somatic cells into pluripotent stem cells, and finally (4) the recent identification of the human pluripotent stem cell essentialome. These four aspects of pluripotency collectively culminate in a broader understanding of what makes a cell pluripotent.

Emergence of a High-Plasticity Cell State during Lung Cancer Evolution

Tumor evolution from a single cell into a malignant, heterogeneous tissue remains poorly understood. Here, we profile single-cell transcriptomes of genetically engineered mouse lung tumors at seven stages, from pre-neoplastic hyperplasia to adenocarcinoma. The diversity of transcriptional states increases over time and is reproducible across tumors and mice. Cancer cells progressively adopt alternate lineage identities, computationally predicted to be mediated through a common transitional, high-plasticity cell state (HPCS). Accordingly, HPCS cells prospectively isolated from mouse tumors and human patient-derived xenografts display high capacity for differentiation and proliferation. The HPCS program is associated with poor survival across human cancers and demonstrates chemoresistance in mice. Our study reveals a central principle underpinning intra-tumoral heterogeneity and motivates therapeutic targeting of the HPCS.

A Rigorous Quantitative Approach to Analyzing Phagocytosis Assays

Studying monocytic cells in isolated systems in vitro contributes significantly to the understanding of innate immune physiology. Functional assays produce read outs which can be used to measure responses to selected stimuli, such as pathogen exposure, antigen loading, and cytokine stimulation. Integration of these results with high quality in vivo models allows for the development of therapeutics which target these cell populations. Current methodologies to quantify phagocytic function of monocytic cells in vitro either measure phagocytic activity of individual cells (average number of beads or particles/cell), or a population outcome (% cells that contain phagocytosed material). Here we address technical challenges and shortcomings of these methods and present a protocol for collecting and analyzing data derived from a functional assay which measures phagocytic activity of macrophage and macrophage-like cells. We apply this method to two different experimental conditions, and compare to existing work flows. We also provide an online tool for users to upload and analyze data using this method.

Cumulus provides cloud-based data analysis for large-scale single-cell and single-nucleus RNA-seq

Massively parallel single-cell and single-nucleus RNA sequencing has opened the way to systematic tissue atlases in health and disease, but as the scale of data generation is growing, so is the need for computational pipelines for scaled analysis. Here we developed Cumulus-a cloud-based framework for analyzing large-scale single-cell and single-nucleus RNA sequencing datasets. Cumulus combines the power of cloud computing with improvements in algorithm and implementation to achieve high scalability, low cost, user-friendliness and integrated support for a comprehensive set of features. We benchmark Cumulus on the Human Cell Atlas Census of Immune Cells dataset of bone marrow cells and show that it substantially improves efficiency over conventional frameworks, while maintaining or improving the quality of results, enabling large-scale studies.

COMUNET: a tool to explore and visualize intercellular communication

MOTIVATION

Intercellular communication plays an essential role in multicellular organisms and several algorithms to analyze it from single-cell transcriptional data have been recently published, but the results are often hard to visualize and interpret.

RESULTS

We developed Cell cOmmunication exploration with MUltiplex NETworks (COMUNET), a tool that streamlines the interpretation of the results from cell-cell communication analyses. COMUNET uses multiplex networks to represent and cluster all potential communication patterns between cell types. The algorithm also enables the search for specific patterns of communication and can perform comparative analysis between two biological conditions. To exemplify its use, here we apply COMUNET to investigate cell communication patterns in single-cell transcriptomic datasets from mouse embryos and from an acute myeloid leukemia patient at diagnosis and after treatment.

AVAILABILITY AND IMPLEMENTATION

Our algorithm is implemented in an R package available from https://github.com/ScialdoneLab/COMUNET, along with all the code to perform the analyses reported here.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computational Methods for Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNA-seq) has provided a high-dimensional catalog of millions of cells across species and diseases. These data have spurred the development of hundreds of computational tools to derive novel biological insights. Here, we outline the components of scRNA-seq analytical pipelines and the computational methods that underlie these steps. We describe available methods, highlight well-executed benchmarking studies, and identify opportunities for additional benchmarking studies and computational methods. As the biochemical approaches for single-cell omics advance, we propose coupled development of robust analytical pipelines suited for the challenges that new data present and principled selection of analytical methods that are suited for the biological questions to be addressed.

Interneuron Types as Attractors and Controllers

Cortical interneurons display striking differences in shape, physiology, and other attributes, challenging us to appropriately classify them. We previously suggested that interneuron types should be defined by their role in cortical processing. Here, we revisit the question of how to codify their diversity based upon their division of labor and function as controllers of cortical information flow. We suggest that developmental trajectories provide a guide for appreciating interneuron diversity and argue that subtype identity is generated using a configurational (rather than combinatorial) code of transcription factors that produce attractor states in the underlying gene regulatory network. We present our updated three-stage model for interneuron specification: an initial cardinal step, allocating interneurons into a few major classes, followed by definitive refinement, creating subclasses upon settling within the cortex, and lastly, state determination, reflecting the incorporation of interneurons into functional circuit ensembles. We close by discussing findings indicating that major interneuron classes are both evolutionarily ancient and conserved. We propose that the complexity of cortical circuits is generated by phylogenetically old interneuron types, complemented by an evolutionary increase in principal neuron diversity. This suggests that a natural neurobiological definition of interneuron types might be derived from a match between their developmental origin and computational function.

Cell profiling of mouse acute kidney injury reveals conserved cellular responses to injury

After acute kidney injury (AKI), patients either recover or alternatively develop fibrosis and chronic kidney disease. Interactions between injured epithelia, stroma, and inflammatory cells determine whether kidneys repair or undergo fibrosis, but the molecular events that drive these processes are poorly understood. Here, we use single nucleus RNA sequencing of a mouse model of AKI to characterize cell states during repair from acute injury. We identify a distinct proinflammatory and profibrotic proximal tubule cell state that fails to repair. Deconvolution of bulk RNA-seq datasets indicates that this failed-repair proximal tubule cell (FR-PTC) state can be detected in other models of kidney injury, increasing during aging in rat kidney and over time in human kidney allografts. We also describe dynamic intercellular communication networks and discern transcriptional pathways driving successful vs. failed repair. Our study provides a detailed description of cellular responses after injury and suggests that the FR-PTC state may represent a therapeutic target to improve repair.

Single cell genomics and developmental biology: moving beyond the generation of cell type catalogues

Major developmental processes such as gastrulation and early embryogenesis rely on a complex network of cell-cell interactions, chromatin remodeling, and transcriptional regulators. This makes it challenging to study early development when using bulk populations of cells. Recent advances in single-cell technologies have allowed researchers to better understand the interactions between different molecular modalities and the heterogeneities within classically defined cell types. As new single-cell technologies mature, they have the potential of providing a step-change in our understanding of embryogenesis. In this review, we summarize recent advances in single-cell technologies with particular focus on those that lend insight into early organogenesis. We then discuss current pitfalls and implications for future research.

Quantile normalization of single-cell RNA-seq read counts without unique molecular identifiers

Single-cell RNA-seq (scRNA-seq) profiles gene expression of individual cells. Unique molecular identifiers (UMIs) remove duplicates in read counts resulting from polymerase chain reaction, a major source of noise. For scRNA-seq data lacking UMIs, we propose quasi-UMIs: quantile normalization of read counts to a compound Poisson distribution empirically derived from UMI datasets. When applied to ground-truth datasets having both reads and UMIs, quasi-UMI normalization has higher accuracy than competing methods. Using quasi-UMIs enables methods designed specifically for UMI data to be applied to non-UMI scRNA-seq datasets.

Comprehensive identification of mRNA isoforms reveals the diversity of neural cell-surface molecules with roles in retinal development and disease

Genes encoding cell-surface proteins control nervous system development and are implicated in neurological disorders. These genes produce alternative mRNA isoforms which remain poorly characterized, impeding understanding of how disease-associated mutations cause pathology. Here we introduce a strategy to define complete portfolios of full-length isoforms encoded by individual genes. Applying this approach to neural cell-surface molecules, we identify thousands of unannotated isoforms expressed in retina and brain. By mass spectrometry we confirm expression of newly-discovered proteins on the cell surface in vivo. Remarkably, we discover that the major isoform of a retinal degeneration gene, CRB1, was previously overlooked. This CRB1 isoform is the only one expressed by photoreceptors, the affected cells in CRB1 disease. Using mouse mutants, we identify a function for this isoform at photoreceptor-glial junctions and demonstrate that loss of this isoform accelerates photoreceptor death. Therefore, our isoform identification strategy enables discovery of new gene functions relevant to disease.

TrajectoryNet: A Dynamic Optimal Transport Network for Modeling Cellular Dynamics

It is increasingly common to encounter data from dynamic processes captured by static cross-sectional measurements over time, particularly in biomedical settings. Recent attempts to model individual trajectories from this data use optimal transport to create pairwise matchings between time points. However, these methods cannot model continuous dynamics and non-linear paths that entities can take in these systems. To address this issue, we establish a link between continuous normalizing flows and dynamic optimal transport, that allows us to model the expected paths of points over time. Continuous normalizing flows are generally under constrained, as they are allowed to take an arbitrary path from the source to the target distribution. We present TrajectoryNet, which controls the continuous paths taken between distributions to produce dynamic optimal transport. We show how this is particularly applicable for studying cellular dynamics in data from single-cell RNA sequencing (scRNA-seq) technologies, and that TrajectoryNet improves upon recently proposed static optimal transport-based models that can be used for interpolating cellular distributions.

Defining Reprogramming Checkpoints from Single-Cell Analyses of Induced Pluripotency

Elucidating the mechanism of reprogramming is confounded by heterogeneity due to the low efficiency and differential kinetics of obtaining induced pluripotent stem cells (iPSCs) from somatic cells. Therefore, we increased the efficiency with a combination of epigenomic modifiers and signaling molecules and profiled the transcriptomes of individual reprogramming cells. Contrary to the established temporal order, somatic gene inactivation and upregulation of cell cycle, epithelial, and early pluripotency genes can be triggered independently such that any combination of these events can occur in single cells. Sustained co-expression of Epcam, Nanog, and Sox2 with other genes is required to progress toward iPSCs. Ehf, Phlda2, and translation initiation factor Eif4a1 play functional roles in robust iPSC generation. Using regulatory network analysis, we identify a critical role for signaling inhibition by 2i in repressing somatic expression and synergy between the epigenomic modifiers ascorbic acid and a Dot1L inhibitor for pluripotency gene activation.

Advancing Stem Cell Research through Multimodal Single-Cell Analysis

Technological advances play a key role in furthering our understanding of stem cell biology, and advancing the prospects of regenerative therapies. Highly parallelized methods, developed in the last decade, can profile DNA, RNA, or proteins in thousands of cells and even capture data across two or more modalities (multiomics). This allows unbiased and precise definition of molecular cell states, thus allowing classification of cell types, tracking of differentiation trajectories, and discovery of underlying mechanisms. Despite being based on destructive techniques, novel experimental and bioinformatic approaches enable embedding and extraction of temporal information, which is essential for deconvolution of complex data and establishing cause and effect relationships. Here, we provide an overview of recent studies pertinent to stem cell biology, followed by an outlook on how further advances in single-cell molecular profiling and computational analysis have the potential to shape the future of both basic and translational research.

Perspectives on somatic reprogramming: spotlighting epigenetic regulation and cellular heterogeneity

Induction of pluripotency in somatic cells via ectopic expression of defined factors demonstrates the reversibility of developmental programming, thus serves as an ideal model to investigate the underlying principles of cell fate determination, which drives the differentiation of various cell types from an identical genome. Over the last decade, our understanding has grown considerably regarding how somatic reprogramming initiates and be regulated. Recent mechanistic investigations on chromatin architecture regulation and cell heterogeneity significantly contribute to the understanding of somatic reprogramming. Specifically, it is found that chromatin dynamics which are regulated by the orchestration of transcriptional factors and epigenetic regulation, play an important role during somatic reprogramming. Moreover, the widespread application of single-cell technologies reveals latent cell fate transition trajectories during the reprogramming process. Here, we highlighted several mechanistic studies on somatic reprogramming, particularly from epigenetic regulation and cell heterogeneity perspectives, and summarized their contribution to our current understanding of cell fate determination.

TFAP2C facilitates somatic cell reprogramming by inhibiting c-Myc-dependent apoptosis and promoting mesenchymal-to-epithelial transition

Transcription factors are known to mediate the conversion of somatic cells to induced pluripotent stem cells (iPSCs). Transcription factor TFAP2C plays important roles in the regulation of embryonic development and carcinogenesis; however, the roles of Tfap2c in regulating somatic cell reprogramming are not well understood. Here we demonstrate Tfap2c is induced during the generation of iPSCs from mouse fibroblasts and acts as a facilitator for iPSCs formation. Mechanistically, the c-Myc-dependent apoptosis, which is a roadblock to reprogramming, can be significantly mitigated by Tfap2c overexpression. Meanwhile, Tfap2c can greatly promote mesenchymal-to-epithelial transition (MET) at initiation stage of OSKM-induced reprogramming. Further analysis of gene expression and targets of Tfap2c during reprogramming by RNA-sequencing (RNA-seq) and ChIP-qPCR indicates that TFAP2C can promote epithelial gene expression by binding to their promoters directly. Finally, knockdown of E-cadherin (Cdh1), an important downstream target of TFAP2C and a critical regulator of MET antagonizes Tfap2c-mediated reprogramming. Taken together, we conclude that Tfap2c serves as a strong activator for somatic cell reprogramming through promoting the MET and inhibiting c-Myc-dependent apoptosis.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.