Computational methods for trajectory inference from single-cell transcriptomics

Optimal sequencing budget allocation for trajectory reconstruction of single cells

BACKGROUND

Charting cellular trajectories over gene expression is key to understanding dynamic cellular processes and their underlying mechanisms. While advances in single-cell RNA-sequencing technologies and computational methods have pushed forward the recovery of such trajectories, trajectory inference remains a challenge due to the noisy, sparse, and high-dimensional nature of single-cell data. This challenge can be alleviated by increasing either the number of cells sampled along the trajectory (breadth) or the sequencing depth, i.e. the number of reads captured per cell (depth). Generally, these two factors are coupled due to an inherent breadth-depth tradeoff that arises when the sequencing budget is constrained due to financial or technical limitations.

RESULTS

Here we study the optimal allocation of a fixed sequencing budget to optimize the recovery of trajectory attributes. Empirical results reveal that reconstruction accuracy of internal cell structure in expression space scales with the logarithm of either the breadth or depth of sequencing. We additionally observe a power law relationship between the optimal number of sampled cells and the corresponding sequencing budget. For linear trajectories, non-monotonicity in trajectory reconstruction across the breadth-depth tradeoff can impact downstream inference, such as expression pattern analysis along the trajectory. We demonstrate these results for five single-cell RNA-sequencing datasets encompassing differentiation of embryonic stem cells, pancreatic beta cells, hepatoblast and multipotent hematopoietic cells, as well as induced reprogramming of embryonic fibroblasts into neurons. By addressing the challenges of single-cell data, our study offers insights into maximizing the efficiency of cellular trajectory analysis through strategic allocation of sequencing resources.

Spatially Self-Organized Three-Dimensional Neural Concentroid as a Novel Reductionist Humanized Model to Study Neurovascular Development

Although human pluripotent stem cell (PSC)-derived brain organoids have enabled researchers to gain insight into human brain development and disease, these organoids contain solely ectodermal cells and are not vascularized as occurs during brain development. Here it is created less complex and more homogenous large neural constructs starting from PSC-derived neuroprogenitor cells (NPC), by fusing small NPC spheroids into so-called concentroids. Such concentroids consisted of a pro-angiogenic core, containing neuronal and outer radial glia cells, surrounded by an astroglia-dense outer layer. Incorporating PSC-derived endothelial cells (EC) around and/or in the concentroids promoted vascularization, accompanied by differential outgrowth and differentiation of neuronal and astroglia cells, as well as the development of ectodermal-derived pericyte-like mural cells co-localizing with EC networks. Single nucleus transcriptomic analysis revealed an enhanced neural cell subtype maturation and diversity in EC-containing concentroids, which better resemble the fetal human brain compared to classical organoids or NPC-only concentroids. This PSC-derived "vascularized" concentroid brain model will facilitate the study of neurovascular/blood-brain barrier development, neural cell migration, and the development of effective in vitro vascularization strategies of brain mimics.

Inferring Tree-Shaped Single-Cell Trajectories with Totem

Single-cell transcriptomics allows unbiased characterization of cell heterogeneity in a sample by profiling gene expression at single-cell level. These profiles capture snapshots of transient or steady states in dynamic processes, such as cell cycle, activation, or differentiation, which can be computationally ordered into a "flip-book" of cell development using trajectory inference methods. However, prediction of more complex topology structures, such as multifurcations or trees, remains challenging. In this chapter, we present two user-friendly protocols for inferring tree-shaped single-cell trajectories and pseudotime from single-cell transcriptomics data with Totem. Totem is a trajectory inference method that offers flexibility in inferring both nonlinear and linear trajectories and usability by avoiding the cumbersome fine-tuning of parameters. The QuickStart protocol provides a simple and practical example, whereas the GuidedStart protocol details the analysis step-by-step. Both protocols are demonstrated using a case study of human bone marrow CD34+ cells, allowing the study of the branching of three lineages: erythroid, lymphoid, and myeloid. All the analyses can be fully reproduced in Linux, macOS, and Windows operating systems (amd64 architecture) with >8 Gb of RAM using the provided docker image distributed with notebooks, scripts, and data in Docker Hub (elolab/repro-totem-ti). These materials are shared online under open-source license at https://elolab.github.io/Totem-protocol .

Interpretable trajectory inference with single-cell Linear Adaptive Negative-binomial Expression (scLANE) testing

The rapid proliferation of trajectory inference methods for single-cell RNA-seq data has allowed researchers to investigate complex biological processes by examining underlying gene expression dynamics. After estimating a latent cell ordering, statistical models are used to determine which genes exhibit changes in expression that are significantly associated with progression through the biological trajectory. While a few techniques for performing trajectory differential expression exist, most rely on the flexibility of generalized additive models in order to account for the inherent nonlinearity of changes in gene expression. As such, the results can be difficult to interpret, and biological conclusions often rest on subjective visual inspections of the most dynamic genes. To address this challenge, we propose scLANE testing, which is built around an interpretable generalized linear model and handles nonlinearity with basis splines chosen empirically for each gene. In addition, extensions to estimating equations and mixed models allow for reliable trajectory testing under complex experimental designs. After validating the accuracy of scLANE under several different simulation scenarios, we apply it to a set of diverse biological datasets and display its ability to provide novel biological information when used downstream of both pseudotime and RNA velocity estimation methods.

A manifold-based framework for studying the dynamics of the vaginal microbiome

The vaginal microbiome plays a crucial role in our health. The composition of this community can be classified into five community state types (CSTs), four of which are primarily consisted of Lactobacillus species and considered healthy, while the fifth features non-Lactobacillus populations and signifies a disease state termed Bacterial vaginosis (BV), which is associated with various symptoms and increased susceptibility to diseases. Importantly, however, the exact mechanisms and dynamics underlying BV development are not yet fully understood, including specifically possible routes from a healthy to a BV state. To address this gap, this study set out to characterize the progression from healthy- to BV-associated compositions by analyzing 8026 vaginal samples and using a manifold-detection framework. This approach, inspired by single-cell analysis, aims to identify low-dimensional trajectories in the high-dimensional composition space. It further orders samples along these trajectories and assigns a score (pseudo-time) to each analyzed or new sample based on its proximity to the BV state. Our results reveal distinct routes of progression between healthy and BV states for each CST, with pseudo-time scores correlating with community diversity and quantifying the health state of each sample. Several BV indicators can also be successfully predicted based on pseudo-time scores, and key taxa involved in BV development can be identified using this approach. Taken together, these findings demonstrate how manifold detection can be used to successfully characterize the progression from healthy Lactobacillus-dominant populations to BV and to accurately quantify the health condition of new samples along the route of BV development.

Elucidating dynamic cell lineages and gene networks in time-course single cell differentiation

Single cell RNA sequencing (scRNA-seq) technologies provide researchers with an unprecedented opportunity to exploit cell heterogeneity. For example, the sequenced cells belong to various cell lineages, which may have different cell fates in stem and progenitor cells. Those cells may differentiate into various mature cell types in a cell differentiation process. To trace the behavior of cell differentiation, researchers reconstruct cell lineages and predict cell fates by ordering cells chronologically into a trajectory with a pseudo-time. However, in scRNA-seq experiments, there are no cell-to-cell correspondences along with the time to reconstruct the cell lineages, which creates a significant challenge for cell lineage tracing and cell fate prediction. Therefore, methods that can accurately reconstruct the dynamic cell lineages and predict cell fates are highly desirable. In this article, we develop an innovative machine-learning framework called Cell Smoothing Transformation (CellST) to elucidate the dynamic cell fate paths and construct gene networks in cell differentiation processes. Unlike the existing methods that construct one single bulk cell trajectory, CellST builds cell trajectories and tracks behaviors for each individual cell. Additionally, CellST can predict cell fates even for less frequent cell types. Based on the individual cell fate trajectories, CellST can further construct dynamic gene networks to model gene-gene relationships along the cell differentiation process and discover critical genes that potentially regulate cells into various mature cell types.

Detection of continuous hierarchical heterogeneity by single-cell surface antigen analysis in the prognosis evaluation of acute myeloid leukaemia

BACKGROUND

Acute myeloid leukaemia (AML) is characterised by the malignant accumulation of myeloid progenitors with a high recurrence rate after chemotherapy. Blasts (leukaemia cells) exhibit a complete myeloid differentiation hierarchy hiding a wide range of temporal information from initial to mature clones, including genesis, phenotypic transformation, and cell fate decisions, which might contribute to relapse in AML patients.

METHODS

Based on the landscape of AML surface antigens generated by mass cytometry (CyTOF), we combined manifold analysis and principal curve-based trajectory inference algorithm to align myelocytes on a single-linear evolution axis by considering their phenotype continuum that correlated with differentiation order. Backtracking the trajectory from mature clusters located automatically at the terminal, we recurred the molecular dynamics during AML progression and confirmed the evolution stage of single cells. We also designed a 'dispersive antigens in neighbouring clusters exhibition (DANCE)' feature selection method to simplify and unify trajectories, which enabled the exploration and comparison of relapse-related traits among 43 paediatric AML bone marrow specimens.

RESULTS

The feasibility of the proposed trajectory analysis method was verified with public datasets. After aligning single cells on the pseudotime axis, primitive clones were recognized precisely from AML blasts, and the expression of the inner molecules before and after drug stimulation was accurately plotted on the trajectory. Applying DANCE to 43 clinical samples with different responses for chemotherapy, we selected 12 antigens as a general panel for myeloblast differentiation performance, and obtain trajectories to those patients. For the trajectories with unified molecular dynamics, CD11c overexpression in the primitive stage indicated a good chemotherapy outcome. Moreover, a later initial peak of stemness heterogeneity tended to be associated with a higher risk of relapse compared with complete remission.

CONCLUSIONS

In this study, pseudotime was generated as a new single-cell feature. Minute differences in temporal traits among samples could be exhibited on a trajectory, thus providing a new strategy for predicting AML relapse and monitoring drug responses over time scale.

A statistical framework for differential pseudotime analysis with multiple single-cell RNA-seq samples

Pseudotime analysis with single-cell RNA-sequencing (scRNA-seq) data has been widely used to study dynamic gene regulatory programs along continuous biological processes. While many methods have been developed to infer the pseudotemporal trajectories of cells within a biological sample, it remains a challenge to compare pseudotemporal patterns with multiple samples (or replicates) across different experimental conditions. Here, we introduce Lamian, a comprehensive and statistically-rigorous computational framework for differential multi-sample pseudotime analysis. Lamian can be used to identify changes in a biological process associated with sample covariates, such as different biological conditions while adjusting for batch effects, and to detect changes in gene expression, cell density, and topology of a pseudotemporal trajectory. Unlike existing methods that ignore sample variability, Lamian draws statistical inference after accounting for cross-sample variability and hence substantially reduces sample-specific false discoveries that are not generalizable to new samples. Using both real scRNA-seq and simulation data, including an analysis of differential immune response programs between COVID-19 patients with different disease severity levels, we demonstrate the advantages of Lamian in decoding cellular gene expression programs in continuous biological processes.

What cannot be seen correctly in 2D visualizations of single-cell 'omics data?

A common strategy for exploring single-cell 'omics data is visualizing 2D nonlinear projections that aim to preserve high-dimensional data properties such as neighborhoods. Alternatively, mathematical theory and other computational tools can directly describe data geometry, while also showing that neighborhoods and other properties cannot be well-preserved in any 2D projection.

A Manifold-Based Framework for Studying the Dynamics of the Vaginal Microbiome

The vaginal bacterial community plays a crucial role in preventing infections. The composition of this community can be classified into five main groups, termed community state types (CSTs). Four of these CSTs, which are primarily consisted of Lactobacillus species, are considered healthy, while the fifth, which is composed of non-Lactobacillus populations, is considered less protective. This latter CST is often considered to represent a state termed Bacterial vaginosis (BV) - a common disease condition associated with unpleasant symptoms and increased susceptibility to sexually transmitted diseases. However, the exact mechanisms underlying BV development are not yet fully understood, including specifically, the dynamics of the vaginal microbiome in BV, and the possible routes it may take from a healthy to a BV state. This study aims to identify the progression from healthy Lactobacillus-dominant populations to symptomatic BV by analyzing 8,026 vaginal samples and using a manifold-detection framework. This approach is inspired by single-cell analysis and aims to identify low-dimensional trajectories in the high-dimensional composition space. This framework further order samples along these trajectories and assign a score (pseudo-time) to each sample based on its proximity to the BV state. Our results reveal distinct routes of progression between healthy and BV state for each CST, with pseudo-time scores correlating with community diversity and quantifying the health state of each sample. BV indicators, including Nugent score, positive Amsel's test, and several Amsel's criteria, can also be successfully predicted based on pseudo-time scores. Additionally, Gardnerella vaginalis can be identified as a key taxon in BV development using this approach, with increased abundance in samples with high pseudo-time, indicating an unhealthier state across all BV-development routes on the manifold. Taken together, these findings demonstrate how manifold detection can be used to successfully characterizes the progression from healthy Lactobacillus-dominant populations to BV and to accurately quantify the health condition of new samples along the route of BV development.

GCSTI: A Single-Cell Pseudotemporal Trajectory Inference Method Based on Graph Compression

The single-cell pseudotemporal trajectory inference is an important way to explore the process of developmental changes within a cell. Due to the uneven rate of cell growth, changes in gene expression depend less on the time of data collection and more on a cell's "internal clock". To overcome the challenges of gene analysis, and replicate biological developmental processes, several strategies have been put forth. However, due to the size of single-cell datasets, locating relevant signposts usually necessitate clustering analysis or a sizable amount of priori information. To this end, we propose a novel single-cell pseudotemporal trajectory inference technique: GCSTI method, which is based on graph compression and doesn't rely on a priori knowledge or clustering procedures, can handle the trajectory inference problem for a large network in a stable and efficient manner. Additionally, we simultaneously improve the pseudotime defining method currently employed in this study in order to obtain more trustworthy and beneficial outcomes for trajectory inference. Finally, we validate the efficacy and stability of the GCSTI method using datasets from human skeletal muscle myogenic cells and four simulated datasets.

Single-cell analysis of lizard blastema fibroblasts reveals phagocyte-dependent activation of Hedgehog-responsive chondrogenesis

Lizards cannot naturally regenerate limbs but are the closest known relatives of mammals capable of epimorphic tail regrowth. However, the mechanisms regulating lizard blastema formation and chondrogenesis remain unclear. Here, single-cell RNA sequencing analysis of regenerating lizard tails identifies fibroblast and phagocyte populations linked to cartilage formation. Pseudotime trajectory analyses suggest spp1+-activated fibroblasts as blastema cell sources, with subsets exhibiting sulf1 expression and chondrogenic potential. Tail blastema, but not limb, fibroblasts express sulf1 and form cartilage under Hedgehog signaling regulation. Depletion of phagocytes inhibits blastema formation, but treatment with pericytic phagocyte-conditioned media rescues blastema chondrogenesis and cartilage formation in amputated limbs. The results indicate a hierarchy of phagocyte-induced fibroblast gene activations during lizard blastema formation, culminating in sulf1+ pro-chondrogenic populations singularly responsive to Hedgehog signaling. These properties distinguish lizard blastema cells from homeostatic and injury-stimulated fibroblasts and indicate potential actionable targets for inducing regeneration in other species, including humans.

Single-Cell Multiomics

Single-cell RNA sequencing methods have led to improved understanding of the heterogeneity and transcriptomic states present in complex biological systems. Recently, the development of novel single-cell technologies for assaying additional modalities, specifically genomic, epigenomic, proteomic, and spatial data, allows for unprecedented insight into cellular biology. While certain technologies collect multiple measurements from the same cells simultaneously, even when modalities are separately assayed in different cells, we can apply novel computational methods to integrate these data. The application of computational integration methods to multimodal paired and unpaired data results in rich information about the identities of the cells present and the interactions between different levels of biology, such as between genetic variation and transcription. In this review, we both discuss the single-cell technologies for measuring these modalities and describe and characterize a variety of computational integration methods for combining the resulting data to leverage multimodal information toward greater biological insight.

HDSCC: A robust clustering approach for Single Cell RNA-seq data using Hyperdimensional Encoding

Significant improvement in Single Cell technologies has given a hand to researchers to measure RNA expression of considerable number of Single Cells simultaneously, resulting in noticeable progress in our knowledge of cellular structure. Microfluidics-based sequencing protocols employing unique molecular identifiers (UMIs) lead to not only high-quality processing but also screening of thousands of cells. However, analysis of said data has caused challenges when it comes to processing time and computational resources as well as analyzing noisy and highly sparse data. Addressing these issues, we proposed a new method to cluster large RNA-seq datasets effectively. In our proposed approach, for having a noise robust clustering, we employed Hyper Dimensional Computing (HDC) approach to analyze Single Cell RNA sequencing data for the first time as best of our knowledge. We compared our results with state-of-the-art works on single-cell clustering and it shows promising performance and robustness in comparison to them. We performed our experiments with a 3.2-GHz CPU and 32 GB of RAM laptop.

Single-cell transcriptomes in the heart: when every epigenome counts

The response of an organ to stimuli emerges from the actions of individual cells. Recent cardiac single-cell RNA-sequencing studies of development, injury, and reprogramming have uncovered heterogeneous populations even among previously well-defined cell types, raising questions about what level of experimental resolution corresponds to disease-relevant, tissue-level phenotypes. In this review, we explore the biological meaning behind this cellular heterogeneity by undertaking an exhaustive analysis of single-cell transcriptomics in the heart (including a comprehensive, annotated compendium of studies published to date) and evaluating new models for the cardiac function that have emerged from these studies (including discussion and schematics that depict new hypotheses in the field). We evaluate the evidence to support the biological actions of newly identified cell populations and debate questions related to the role of cell-to-cell variability in development and disease. Finally, we present emerging epigenomic approaches that, when combined with single-cell RNA-sequencing, can resolve basic mechanisms of gene regulation and variability in cell phenotype.

A Clustering Method Unifying Cell-Type Recognition and Subtype Identification for Tumor Heterogeneity Analysis

The rapid development of single-cell technology has opened up a whole new perspective for identifying cell types in multicellular organisms and understanding the relationships between them. Distinguishing different cell types and subtypes can identify the components of different immune cells and different tumor clones in the tumor microenvironment, which is the basic work of tumor heterogeneity analysis and can help researchers understand the mechanism of tumor immune escape. Existing algorithms treat both cell types and subtypes as populations of cells with specific gene expression patterns, which is not conducive to accurate cell typing. For that, we proposed a cell similarity metric that unifies cell type recognition and subtype identification (UCRSI), with the assumption that selectively expressed genes represent differences in underlying cell type with on/off manner, while differences in expression level represent different cell subtype with more/less manner. Our method calculates these two kinds of differences separately, and then combines them using a consensus adjacency matrix, and finally cell typing is completed using spectral clustering algorithm. The results show that UCRSI can reconstruct expert annotation of single-cell RNA sequencing datasets more robustly than existing methods. And, UCRSI is useful for analyzing tumor heterogeneity and improving visualization of large-scale cell clustering.

siVAE: interpretable deep generative models for single-cell transcriptomes

Neural networks such as variational autoencoders (VAE) perform dimensionality reduction for the visualization and analysis of genomic data, but are limited in their interpretability: it is unknown which data features are represented by each embedding dimension. We present siVAE, a VAE that is interpretable by design, thereby enhancing downstream analysis tasks. Through interpretation, siVAE also identifies gene modules and hubs without explicit gene network inference. We use siVAE to identify gene modules whose connectivity is associated with diverse phenotypes such as iPSC neuronal differentiation efficiency and dementia, showcasing the wide applicability of interpretable generative models for genomic data analysis.

Dynamic cellular and molecular characteristics of spermatogenesis in the viviparous marine teleost Sebastes schlegelii†

Spermatogenesis is a dynamic cell developmental process that is essential for reproductive success. Vertebrates utilize a variety of reproductive strategies, including sperm diversity, and internal and external fertilization. Research on the cellular and molecular dynamic changes involved in viviparous teleost spermatogenesis, however, is currently lacking. Here, we combined cytohistology, 10 × genomic single-cell RNA-seq, and transcriptome technology to determine the dynamic development characteristics of the spermatogenesis of Sebastes schlegelii. The expressions of lhcgr (Luteinizing hormone/Choriogonadotropin receptor), fshr (follicle-stimulating hormone receptor), ar (androgen receptor), pgr (progesterone receptor), and cox (cyclo-oxygen-ase), as well as the prostaglandin E and F levels peaked during the maturation period, indicating that they were important for sperm maturation and mating. Fifteen clusters were identified based on the 10 × genomic single-cell results. The cell markers of the sub-cluster were identified by their upregulation; piwil, dazl, and dmrt1 were upregulated and identified as spermatogonium markers, and sycp1/3 and spo11 were identified as spermatocyte markers. For S. schlegelii, the sperm head nucleus was elongated (spherical to streamlined in shape), which is a typical characteristic for sperm involved in internal fertilization. We also identified a series of crucial genes associated with spermiogenesis, such as spata6, spag16, kif20a, trip10, and klf10, while kif2c, kifap3, fez2, and spaca6 were found to be involved in nucleus elongation. The results of this study will enrich our cellular and molecular knowledge of spermatogenesis and spermiogenesis in fish that undergo internal fertilization.

Reconstruction of Single-Cell Trajectories Using Stochastic Tree Search

The recent advancement in single-cell RNA sequencing technologies enables the understanding of dynamic cellular processes at the single-cell level. Using trajectory inference methods, pseudotimes can be estimated based on reconstructed single-cell trajectories which can be further used to gain biological knowledge. Existing methods for modeling cell trajectories, such as minimal spanning tree or k-nearest neighbor graph, often lead to locally optimal solutions. In this paper, we propose a penalized likelihood-based framework and introduce a stochastic tree search (STS) algorithm aiming at the global solution in a large and non-convex tree space. Both simulated and real data experiments show that our approach is more accurate and robust than other existing methods in terms of cell ordering and pseudotime estimation.

Bulk and Single-Cell RNA-Seq Analyses for Studies of Spermatogonia

Robust methods have been developed that leverage next-generation sequencing (NGS) to measure abundance of all mRNAs (RNA-seq) in samples as small as individual cells in order to study the testicular transcriptome in mammals. In this chapter, we present robust options for implementing bioinformatics workflows for the analysis of bulk RNA-seq from aggregate samples of hundreds to millions of cells and single-cell RNA-seq from individual cells. We also provide detailed protocols for using the R packages DESeq2 and Seurat, important parameters for successful implementation, and considerations for drawing conclusions from the results.

From pseudo to real-time dynamics of T cell thymic differentiation

Numerous methods have recently emerged for ordering single cells along developmental trajectories. However, accurate depiction of developmental dynamics can only be achieved after rescaling the trajectory according to the relative time spent at each developmental point. We formulate a model which estimates local cell densities and fluxes, and incorporates cell division and apoptosis rates, to infer the real-time dimension of the developmental trajectory. We validate the model using mathematical simulations and apply it to experimental high dimensional cytometry data obtained from the mouse thymus to construct the true time profile of the thymocyte developmental process. Our method can easily be implemented in any of the existing tools for trajectory inference.

Unified epigenomic, transcriptomic, proteomic, and metabolomic taxonomy of Alzheimer's disease progression and heterogeneity

Alzheimer's disease (AD) is a heterogeneous disorder with abnormalities in multiple biological domains. In an advanced machine learning analysis of postmortem brain and in vivo blood multi-omics molecular data (N = 1863), we integrated epigenomic, transcriptomic, proteomic, and metabolomic profiles into a multilevel biological AD taxonomy. We obtained a personalized multilevel molecular index of AD dementia progression that predicts severity of neuropathologies, and identified three robust molecular-based subtypes that explain much of the pathologic and clinical heterogeneity of AD. These subtypes present distinct patterns of alteration in DNA methylation, RNA, proteins, and metabolites, identifiable in the brain and subsequently in blood. In addition, the genetic variations that predispose to the various AD subtypes in brain predict distinct spatial patterns of alteration in cell types, suggesting a unique influence of each putative AD variant on neuropathological mechanisms. These observations support that an individually tailored multi-omics molecular taxonomy of AD may represent distinct targets for preventive or treatment interventions.

Cascade Forest-Based Model for Prediction of RNA Velocity

In recent years, single-cell RNA sequencing technology (scRNA-seq) has developed rapidly and has been widely used in biological and medical research, such as in expression heterogeneity and transcriptome dynamics of single cells. The investigation of RNA velocity is a new topic in the study of cellular dynamics using single-cell RNA sequencing data. It can recover directional dynamic information from single-cell transcriptomics by linking measurements to the underlying dynamics of gene expression. Predicting the RNA velocity vector of each cell based on its gene expression data and formulating RNA velocity prediction as a classification problem is a new research direction. In this paper, we develop a cascade forest model to predict RNA velocity. Compared with other popular ensemble classifiers, such as XGBoost, RandomForest, LightGBM, NGBoost, and TabNet, it performs better in predicting RNA velocity. This paper provides guidance for researchers in selecting and applying appropriate classification tools in their analytical work and suggests some possible directions for future improvement of classification tools.

Single-cell and single-nuclei RNA sequencing as powerful tools to decipher cellular heterogeneity and dysregulation in neurodegenerative diseases

Neurodegenerative diseases affect millions of people worldwide and there are currently no cures. Two types of common neurodegenerative diseases are Alzheimer's (AD) and Parkinson's disease (PD). Single-cell and single-nuclei RNA sequencing (scRNA-seq and snRNA-seq) have become powerful tools to elucidate the inherent complexity and dynamics of the central nervous system at cellular resolution. This technology has allowed the identification of cell types and states, providing new insights into cellular susceptibilities and molecular mechanisms underlying neurodegenerative conditions. Exciting research using high throughput scRNA-seq and snRNA-seq technologies to study AD and PD is emerging. Herein we review the recent progress in understanding these neurodegenerative diseases using these state-of-the-art technologies. We discuss the fundamental principles and implications of single-cell sequencing of the human brain. Moreover, we review some examples of the computational and analytical tools required to interpret the extensive amount of data generated from these assays. We conclude by highlighting challenges and limitations in the application of these technologies in the study of AD and PD.

Alignment of single-cell trajectory trees with CAPITAL

Global alignment of complex pseudotime trajectories between different single-cell RNA-seq datasets is challenging, as existing tools mainly focus on linear alignment of single-cell trajectories. Here we present CAPITAL (comparative analysis of pseudotime trajectory inference with tree alignment), a method for comparing single-cell trajectories with tree alignment whereby branching trajectories can be automatically compared. Computational tests on synthetic datasets and authentic bone marrow cells datasets indicate that CAPITAL has achieved accurate and robust alignments of trajectory trees, revealing various gene expression dynamics including gene–gene correlation conservation between different species.

Global alignment of complex cell state trajectories between single-cell datasets remains challenging. Here, the authors present a computational method called CAPITAL to compare branching trajectories, and demonstrate that this method achieves accurate and robust alignments.

Principles of human and mouse nephron development

The mechanisms underlying kidney development in mice and humans is an area of intense study. Insights into kidney organogenesis have the potential to guide our understanding of the origin of congenital anomalies and enable the assembly of genetic diagnostic tools. A number of studies have delineated signalling nodes that regulate positional identities and cell fates of nephron progenitor and precursor cells, whereas cross-species comparisons have markedly enhanced our understanding of conserved and divergent features of mammalian kidney organogenesis. Greater insights into the complex cellular movements that occur as the proximal-distal axis is established have challenged our understanding of nephron patterning and provided important clues to the elaborate developmental context in which human kidney diseases can arise. Studies of kidney development in vivo have also facilitated efforts to recapitulate nephrogenesis in kidney organoids in vitro, by providing a detailed blueprint of signalling events, cell movements and patterning mechanisms that are required for the formation of correctly patterned nephrons and maturation of physiologically functional apparatus that are responsible for maintaining human health.

Prognostic Significance of Pseudotime from Texture Parameters of FDG PET/CT in Locally Advanced Non-Small-Cell Lung Cancer with Tri-Modality Therapy

Simple Summary

Although texture parameters of F-18 fluorodeoxyglucose positron emission tomography/computed tomography images were known to associate tumor biology and clinical features, the types and implications of parameters are too various and complicated. To overcome the limitation of texture parameter, we attempted to produce a new simplified parameter from texture parameters of F-18 fluorodeoxyglucose positron emission tomography/computed tomography images in lung cancer patients using pseudotime analysis. Pseudotime analysis is a recently developed method to explore changes in cell or tissue characteristics based on transcriptomic expression. It is the first study to apply pseudotime analysis into radiomics dataset other than transcriptomics data. Herein, we demonstrated that pseudotime can be successfully estimated from texture parameters. In the aspect of prognostic prediction, pseudotime was an independent prognostic factor for overall survival in contrast to conventional parameters such as metabolic tumor volume and total lesion glycolysis. This study showed possibility of integrating various texture parameters into single parameter which reflects disease progression status. Pseudotime, as a concrete value of disease progression, is expected to be used in clinical field to evaluate disease and predict prognosis.

Abstract

Texture analysis provides image parameters from F-18 fluorodeoxyglucose positron emission tomography/computed tomography (FDG PET/CT). Although some parameters are associated with tumor biology and clinical features, the types and implications of these parameters are complicated. We applied pseudotime analysis, which has recently been used to estimate changes in individual sample characteristics, to texture parameters from FDG PET/CT images of locally advanced non-small-cell lung cancer (NSCLC) patients undergoing neoadjuvant concurrent chemoradiation therapy (CCRT) followed by surgery. Our subjects were 303 NSCLC patients who underwent pretherapeutic FDG PET/CT and tri-modality therapy. Texture parameters of the primary tumor were calculated from FDG PET/CT images acquired before neoadjuvant CCRT. Pseudotime analysis was performed using the PhenoPath tool. Clinicopathologic features including survival data were collected and survival analysis was performed to compare the prognostic significances of pseudotime parameters with those of conventional PET parameters. Pseudotime was successfully estimated from texture parameters. Normalized co-occurrence homogeneity, normalized co-occurrence inverse difference moment, and black–white symmetry showed positive correlations with pseudotime, short run emphasis, normalized co-occurrence dissimilarity, and short zone emphasis negative correlation. The maximum standardized uptake value (SUV) and mean SUV were not associated with overall survival. Pseudotime, metabolic tumor volume (MTV), and total lesion glycolysis (TLG) showed significant associations with overall survival. In contrast to MTV and TLG, pseudotime was an independent prognostic factor for overall survival. Various metabolic texture parameters can be integrated into a single parameter using pseudotime analysis. Pseudotime of the primary tumor, estimated from FDG PET/CT images, better predicts overall survival in locally advanced NSCLC patients treated with tri-modality therapy than conventional PET parameters.

Discovering Immune-Mediated Mechanisms of Gastric Carcinogenesis Through Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNAseq) technology is still relatively new in the field of gastric cancer immunology but gaining significant traction. This technology now provides unprecedented insights into the intratumoral and intertumoral heterogeneities at the immunological, cellular, and molecular levels. Within the last few years, a volume of publications reported the usefulness of scRNAseq technology in identifying thus far elusive immunological mechanisms that may promote and impede gastric cancer development. These studies analyzed datasets generated from primary human gastric cancer tissues, metastatic ascites fluid from gastric cancer patients, and laboratory-generated data from in vitro and in vivo models of gastric diseases. In this review, we overview the exciting findings from scRNAseq datasets that uncovered the role of critical immune cells, including T cells, B cells, myeloid cells, mast cells, ILC2s, and other inflammatory stromal cells, like fibroblasts and endothelial cells. In addition, we also provide a synopsis of the initial scRNAseq findings on the interesting epithelial cell responses to inflammation. In summary, these new studies have implicated roles for T and B cells and subsets like NKT cells in tumor development and progression. The current studies identified diverse subsets of macrophages and mast cells in the tumor microenvironment, however, additional studies to determine their roles in promoting cancer growth are needed. Some groups specifically focus on the less prevalent ILC2 cell type that may contribute to early cancer development. ScRNAseq analysis also reveals that stromal cells, e.g., fibroblasts and endothelial cells, regulate inflammation and promote metastasis, making them key targets for future investigations. While evaluating the outcomes, we also highlight the gaps in the current findings and provide an assessment of what this technology holds for gastric cancer research in the coming years. With scRNAseq technology expanding rapidly, we stress the need for periodic review of the findings and assess the available scRNAseq analytical tools to guide future work on immunological mechanisms of gastric carcinogenesis.  

Energy landscape decomposition for cell differentiation with proliferation effect

Complex interactions between genes determine the development and differentiation of cells. We establish a landscape theory for cell differentiation with proliferation effect, in which the developmental process is modeled as a stochastic dynamical system with a birth-death term. We find that two different energy landscapes, denoted U and V, collectively contribute to the establishment of non-equilibrium steady differentiation. The potential U is known as the energy landscape leading to the steady distribution, whose metastable states stand for cell types, while V indicates the differentiation direction from pluripotent to differentiated cells. This interpretation of cell differentiation is different from the previous landscape theory without the proliferation effect. We propose feasible numerical methods and a mean-field approximation for constructing landscapes U and V. Successful applications to typical biological models demonstrate the energy landscape decomposition's validity and reveal biological insights into the considered processes.

Investigating the Clinico-Molecular and Immunological Evolution of Lung Adenocarcinoma Using Pseudotime Analysis

Introduction

As the molecular features of lung adenocarcinoma (LUAD) have been evaluated as a cross-sectional study, the course of tumor characteristics has not been modeled. The temporal evolution of the tumor immune microenvironment (TIME), as well as the clinico-molecular features of LUAD, could provide a precise strategy for immunotherapy and surrogate biomarkers for the course of LUAD.

Methods

A pseudotime trajectory was constructed in patients with LUAD from the Cancer Genome Atlas and non-small cell lung cancer radiogenomics datasets. Correlation analyses were performed between clinical features and pseudotime. Genes associated with pseudotime were selected, and gene ontology analysis was performed. F-18 fluorodeoxyglucose positron emission tomography images of subjects were collected, and imaging parameters, including standardized uptake value (SUV), were obtained. Correlation analyses were performed between imaging parameters and pseudotime. Correlation analyses were performed between the enrichment scores of various immune cell types and pseudotime. In addition, correlation analyses were performed between the expression of PD-L1, tumor mutation burden, and pseudotime.

Results

Pseudotime trajectories of LUAD corresponded to clinical stages. Molecular profiles related to cell division and natural killer cell activity were changed along the pseudotime. The maximal SUV of LUAD tumors showed a positive correlation with pseudotime. Type 1 helper T (Th1) cells showed a positive correlation, whereas M2 macrophages showed a negative correlation with pseudotime. PD-L1 expression showed a negative correlation, whereas tumor mutation burden showed a positive correlation with pseudotime.

Conclusion

The estimated pseudotime associated with the stage suggested that it could reflect the clinico-molecular evolution of LUAD. Specific immune cell types in the TIME as well as cell division and glucose metabolism were dynamically changed according to the progression of the pseudotime. As a molecular progression of LUAD, different cellular targets should be considered for immunotherapy.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.

Pseudotime Analysis Reveals Exponential Trends in DNA Methylation Aging with Mortality Associated Timescales

The epigenetic trajectory of DNA methylation profiles has a nonlinear relationship with time, reflecting rapid changes in DNA methylation early in life that progressively slow with age. In this study, we use pseudotime analysis to determine the functional form of these trajectories. Unlike epigenetic clocks that constrain the functional form of methylation changes with time, pseudotime analysis orders samples along a path, based on similarities in a latent dimension, to provide an unbiased trajectory. We show that pseudotime analysis can be applied to DNA methylation in human blood and brain tissue and find that it is highly correlated with the epigenetic states described by the Epigenetic Pacemaker. Moreover, we show that the pseudotime trajectory can be modeled with respect to time, using a sum of two exponentials, with coefficients that are close to the timescales of human age-associated mortality. Thus, for the first time, we can identify age-associated molecular changes that appear to track the exponential dynamics of mortality risk.

Revealing new biology from multiplexed, metal-isotope-tagged, single-cell readouts

Mass cytometry (MC) is a recent technology that pairs plasma-based ionization of cells in suspension with time-of-flight (TOF) mass spectrometry to sensitively quantify the single-cell abundance of metal-isotope-tagged affinity reagents to key proteins, RNA, and peptides. Given the ability to multiplex readouts (~50 per cell) and capture millions of cells per experiment, MC offers a robust way to assay rare, transitional cell states that are pertinent to human development and disease. Here, we review MC approaches that let us probe the dynamics of cellular regulation across multiple conditions and sample types in a single experiment. Additionally, we discuss current limitations and future extensions of MC as well as computational tools commonly used to extract biological insight from single-cell proteomic datasets.

A Bayesian factorization method to recover single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) offers opportunities to study gene expression of tens of thousands of single cells simultaneously, to investigate cell-to-cell variation, and to reconstruct cell-type-specific gene regulatory networks. Recovering dropout events in a sparse gene expression matrix for scRNA-seq data is a long-standing matrix completion problem. In this article, we introduce Bfimpute, a Bayesian factorization imputation algorithm that reconstructs two latent gene and cell matrices to impute the final gene expression matrix within each cell group, with or without the aid of cell type labels or bulk data. Bfimpute achieves better accuracy than ten other publicly notable scRNA-seq imputation methods on simulated and real scRNA-seq data, as measured by several different evaluation metrics. Bfimpute can also flexibly integrate any gene- or cell-related information that users provide to increase performance.

Dictionary learning allows model-free pseudotime estimation of transcriptomic data

Background

Pseudotime estimation from dynamic single-cell transcriptomic data enables characterisation and understanding of the underlying processes, for example developmental processes. Various pseudotime estimation methods have been proposed during the last years. Typically, these methods start with a dimension reduction step because the low-dimensional representation is usually easier to analyse. Approaches such as PCA, ICA or t-SNE belong to the most widely used methods for dimension reduction in pseudotime estimation methods. However, these methods usually make assumptions on the derived dimensions, which can result in important dataset properties being missed. In this paper, we suggest a new dictionary learning based approach, dynDLT, for dimension reduction and pseudotime estimation of dynamic transcriptomic data. Dictionary learning is a matrix factorisation approach that does not restrict the dependence of the derived dimensions. To evaluate the performance, we conduct a large simulation study and analyse 8 real-world datasets.

Results

The simulation studies reveal that firstly, dynDLT preserves the simulated patterns in low-dimension and the pseudotimes can be derived from the low-dimensional representation. Secondly, the results show that dynDLT is suitable for the detection of genes exhibiting the simulated dynamic patterns, thereby facilitating the interpretation of the compressed representation and thus the dynamic processes. For the real-world data analysis, we select datasets with samples that are taken at different time points throughout an experiment. The pseudotimes found by dynDLT have high correlations with the experimental times. We compare the results to other approaches used in pseudotime estimation, or those that are method-wise closely connected to dictionary learning: ICA, NMF, PCA, t-SNE, and UMAP. DynDLT has the best overall performance for the simulated and real-world datasets.

Conclusions

We introduce dynDLT, a method that is suitable for pseudotime estimation. Its main advantages are: (1) It presents a model-free approach, meaning that it does not restrict the dependence of the derived dimensions; (2) Genes that are relevant in the detected dynamic processes can be identified from the dictionary matrix; (3) By a restriction of the dictionary entries to positive values, the dictionary atoms are highly interpretable.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-08276-9).

Acute Infectious Gastroenteritis: The Causative Agents, Omics-Based Detection of Antigens and Novel Biomarkers

Acute infectious gastroenteritis (AGE) is among the leading causes of mortality in children less than 5 years of age worldwide. There are many causative agents that lead to this infection, with rotavirus being the commonest pathogen in the past decade. However, this trend is now being progressively replaced by another agent, which is the norovirus. Apart from the viruses, bacteria such as Salmonella and Escherichia coli and parasites such as Entamoeba histolytica also contribute to AGE. These agents can be recognised by their respective biological markers, which are mainly the specific antigens or genes to determine the causative pathogen. In conjunction to that, omics technologies are currently providing crucial insights into the diagnosis of acute infectious gastroenteritis at the molecular level. Recent advancement in omics technologies could be an important tool to further elucidate the potential causative agents for AGE. This review will explore the current available biomarkers and antigens available for the diagnosis and management of the different causative agents of AGE. Despite the high-priced multi-omics approaches, the idea for utilization of these technologies is to allow more robust discovery of novel antigens and biomarkers related to management AGE, which eventually can be developed using easier and cheaper detection methods for future clinical setting. Thus, prediction of prognosis, virulence and drug susceptibility for active infections can be obtained. Case management, risk prediction for hospital-acquired infections, outbreak detection, and antimicrobial accountability are aimed for further improvement by integrating these capabilities into a new clinical workflow.

Matrix factorization for biomedical link prediction and scRNA-seq data imputation: an empirical survey

Advances in high-throughput experimental technologies promote the accumulation of vast number of biomedical data. Biomedical link prediction and single-cell RNA-sequencing (scRNA-seq) data imputation are two essential tasks in biomedical data analyses, which can facilitate various downstream studies and gain insights into the mechanisms of complex diseases. Both tasks can be transformed into matrix completion problems. For a variety of matrix completion tasks, matrix factorization has shown promising performance. However, the sparseness and high dimensionality of biomedical networks and scRNA-seq data have raised new challenges. To resolve these issues, various matrix factorization methods have emerged recently. In this paper, we present a comprehensive review on such matrix factorization methods and their usage in biomedical link prediction and scRNA-seq data imputation. Moreover, we select representative matrix factorization methods and conduct a systematic empirical comparison on 15 real data sets to evaluate their performance under different scenarios. By summarizing the experimental results, we provide general guidelines for selecting matrix factorization methods for different biomedical matrix completion tasks and point out some future directions to further improve the performance for biomedical link prediction and scRNA-seq data imputation.

Molecular and Cellular Dynamics of Aortic Aneurysms Revealed by Single-Cell Transcriptomics

The aorta is highly heterogeneous, containing many different types of cells that perform sophisticated functions to maintain aortic homeostasis. Recently, single-cell RNA sequencing studies have provided substantial new insight into the heterogeneity of vascular cell types, the comprehensive molecular features of each cell type, and the phenotypic interrelationship between these cell populations. This new information has significantly improved our understanding of aortic biology and aneurysms at the molecular and cellular level. Here, we summarize these findings, with a focus on what single-cell RNA sequencing analysis has revealed about cellular heterogeneity, cellular transitions, communications among cell populations, and critical transcription factors in the vascular wall. We also review the information learned from single-cell RNA sequencing that has contributed to our understanding of the pathogenesis of vascular disease, such as the identification of cell types in which aneurysm-related genes and genetic variants function. Finally, we discuss the challenges and future directions of single-cell RNA sequencing applications in studies of aortic biology and diseases.

The San Diego Nathan Shock Center: tackling the heterogeneity of aging

Understanding basic mechanisms of aging holds great promise for developing interventions that prevent or delay many age-related declines and diseases simultaneously to increase human healthspan. However, a major confounding factor in aging research is the heterogeneity of the aging process itself. At the organismal level, it is clear that chronological age does not always predict biological age or susceptibility to frailty or pathology. While genetics and environment are major factors driving variable rates of aging, additional complexity arises because different organs, tissues, and cell types are intrinsically heterogeneous and exhibit different aging trajectories normally or in response to the stresses of the aging process (e.g., damage accumulation). Tackling the heterogeneity of aging requires new and specialized tools (e.g., single-cell analyses, mass spectrometry-based approaches, and advanced imaging) to identify novel signatures of aging across scales. Cutting-edge computational approaches are then needed to integrate these disparate datasets and elucidate network interactions between known aging hallmarks. There is also a need for improved, human cell-based models of aging to ensure that basic research findings are relevant to human aging and healthspan interventions. The San Diego Nathan Shock Center (SD-NSC) provides access to cutting-edge scientific resources to facilitate the study of the heterogeneity of aging in general and to promote the use of novel human cell models of aging. The center also has a robust Research Development Core that funds pilot projects on the heterogeneity of aging and organizes innovative training activities, including workshops and a personalized mentoring program, to help investigators new to the aging field succeed. Finally, the SD-NSC participates in outreach activities to educate the general community about the importance of aging research and promote the need for basic biology of aging research in particular.

Spatially distinct reprogramming of the tumor microenvironment based on tumor invasion in diffuse-type gastric cancers

PURPOSE

Histological features of diffuse-type gastric cancer (GC) indicate that the tumor microenvironment (TME) may substantially impact tumor invasiveness. However, cellular components and molecular features associated with cancer invasiveness in the TME of diffuse-type GCs are poorly understood.

EXPERIMENTAL DESIGN

We performed single-cell RNA-sequencing (scRNA-seq) using tissue samples from superficial and deep invasive layers of cancerous and paired normal tissues freshly harvested from five patients with diffuse-type GC. The scRNA-seq results were validated by immunohistochemistry and duplex in situ hybridization (ISH) in formalin-fixed paraffin-embedded tissues.

RESULTS

Seven major cell types were identified. Fibroblasts, endothelial cells, and myeloid cells were categorised as being enriched in the deep layers. Cell type-specific clustering further revealed that the superficial-to-deep layer transition is associated with enrichment in inflammatory endothelial cells and fibroblasts with upregulated CCL2 transcripts. Immunohistochemistry and duplex ISH revealed the distribution of the major cell types and CCL2-expressing endothelial cells and fibroblasts, indicating tumor invasion. Elevation of CCL2 levels along the superficial-to-deep layer axis revealed the immunosuppressive immune cell sub-types that may contribute to tumor cell aggressiveness in the deep invasive layers of diffuse-type GC. The analyses of public datasets revealed the high-level co-expression of stromal cell-specific genes and that CCL2 correlated with poor survival outcomes in GC patients.

CONCLUSIONS

This study reveals the spatial reprogramming of the TME that may underlie invasive tumor potential in diffuse-type GC. This TME profiling across tumor layers suggests new targets, such as CCL2, that can modify the TME to inhibit tumor progression in diffuse-type GC.

NovoSpaRc: flexible spatial reconstruction of single-cell gene expression with optimal transport

Single-cell RNA-sequencing (scRNA-seq) technologies have revolutionized modern biomedical sciences. A fundamental challenge is to incorporate spatial information to study tissue organization and spatial gene expression patterns. Here, we describe a detailed protocol for using novoSpaRc, a computational framework that probabilistically assigns cells to tissue locations. At the core of this framework lies a structural correspondence hypothesis, that cells in physical proximity share similar gene expression profiles. Given scRNA-seq data, novoSpaRc spatially reconstructs tissues based on this hypothesis, and optionally, by including a reference atlas of marker genes to improve reconstruction. We describe the novoSpaRc algorithm, and its implementation in an open-source Python package ( https://pypi.org/project/novosparc ). NovoSpaRc maps a scRNA-seq dataset of 10,000 cells onto 1,000 locations in <5 min. We describe results obtained using novoSpaRc to reconstruct the mouse organ of Corti de novo based on the structural correspondence assumption and human osteosarcoma cultured cells based on marker gene information, and provide a step-by-step guide to Drosophila embryo reconstruction in the Procedure to demonstrate how these two strategies can be combined.

The Trifecta of Single-Cell, Systems-Biology, and Machine-Learning Approaches

Together, single-cell technologies and systems biology have been used to investigate previously unanswerable questions in biomedicine with unparalleled detail. Despite these advances, gaps in analytical capacity remain. Machine learning, which has revolutionized biomedical imaging analysis, drug discovery, and systems biology, is an ideal strategy to fill these gaps in single-cell studies. Machine learning additionally has proven to be remarkably synergistic with single-cell data because it remedies unique challenges while capitalizing on the positive aspects of single-cell data. In this review, we describe how systems-biology algorithms have layered machine learning with biological components to provide systems level analyses of single-cell omics data, thus elucidating complex biological mechanisms. Accordingly, we highlight the trifecta of single-cell, systems-biology, and machine-learning approaches and illustrate how this trifecta can significantly contribute to five key areas of scientific research: cell trajectory and identity, individualized medicine, pharmacology, spatial omics, and multi-omics. Given its success to date, the systems-biology, single-cell omics, and machine-learning trifecta has proven to be a potent combination that will further advance biomedical research.

scGPS: Determining Cell States and Global Fate Potential of Subpopulations

Finding cell states and their transcriptional relatedness is a main outcome from analysing single-cell data. In developmental biology, determining whether cells are related in a differentiation lineage remains a major challenge. A seamless analysis pipeline from cell clustering to estimating the probability of transitions between cell clusters is lacking. Here, we present Single Cell Global fate Potential of Subpopulations (scGPS) to characterise transcriptional relationship between cell states. scGPS decomposes mixed cell populations in one or more samples into clusters (SCORE algorithm) and estimates pairwise transitioning potential (scGPS algorithm) of any pair of clusters. SCORE allows for the assessment and selection of stable clustering results, a major challenge in clustering analysis. scGPS implements a novel approach, with machine learning classification, to flexibly construct trajectory connections between clusters. scGPS also has a feature selection functionality by network and modelling approaches to find biological processes and driver genes that connect cell populations. We applied scGPS in diverse developmental contexts and show superior results compared to a range of clustering and trajectory analysis methods. scGPS is able to identify the dynamics of cellular plasticity in a user-friendly workflow, that is fast and memory efficient. scGPS is implemented in R with optimised functions using C++ and is publicly available in Bioconductor.

A statistical framework for differential pseudotime analysis with multiple single-cell RNA-seq samples

Pseudotime analysis with single-cell RNA-sequencing (scRNA-seq) data has been widely used to study dynamic gene regulatory programs along continuous biological processes. While many computational methods have been developed to infer the pseudo-temporal trajectories of cells within a biological sample, methods that compare pseudo-temporal patterns with multiple samples (or replicates) across different experimental conditions are lacking. Lamian is a comprehensive and statistically-rigorous computational framework for differential multi-sample pseudotime analysis. It can be used to identify changes in a biological process associated with sample covariates, such as different biological conditions, and also to detect changes in gene expression, cell density, and topology of a pseudotemporal trajectory. Unlike existing methods that ignore sample variability, Lamian draws statistical inference after accounting for cross-sample variability and hence substantially reduces sample-specific false discoveries that are not generalizable to new samples. Using both simulations and real scRNA-seq data, including an analysis of differential immune response programs between COVID-19 patients with different disease severity levels, we demonstrate the advantages of Lamian in decoding cellular gene expression programs in continuous biological processes.

PseudoGA: cell pseudotime reconstruction based on genetic algorithm

Dynamic regulation of gene expression is often governed by progression through transient cell states. Bulk RNA-seq analysis can only detect average change in expression levels and is unable to identify this dynamics. Single cell RNA-seq presents an unprecedented opportunity that helps in placing the cells on a hypothetical time trajectory that reflects gradual transition of their transcriptomes. This continuum trajectory or ‘pseudotime’, may reveal the developmental pathway and provide us with information on dynamic transcriptomic changes and other biological processes. Existing approaches to build pseudotime heavily depend on reducing huge dimension to extremely low dimensional subspaces and may lead to loss of information. We propose PseudoGA, a genetic algorithm based approach to order cells assuming that gene expressions vary according to a smooth curve along the pseudotime trajectory. We observe superior accuracy of our method in simulated as well as benchmarking real datasets. Generality of the assumption behind PseudoGA and no dependence on dimensionality reduction technique make it a robust choice for pseudotime estimation from single cell transcriptome data. PseudoGA is also time efficient when applied to a large single cell RNA-seq data and adaptable to parallel computing. R code for PseudoGA is freely available at https://github.com/indranillab/pseudoga.

Nasal ciliated cells are primary targets for SARS-CoV-2 replication in the early stage of COVID-19

The upper respiratory tract is compromised in the early period of COVID-19, but SARS-CoV-2 tropism at the cellular level is not fully defined. Unlike recent single-cell RNA-Seq analyses indicating uniformly low mRNA expression of SARS-CoV-2 entry-related host molecules in all nasal epithelial cells, we show that the protein levels are relatively high and that their localizations are restricted to the apical side of multiciliated epithelial cells. In addition, we provide evidence in patients with COVID-19 that SARS-CoV-2 is massively detected and replicated within the multiciliated cells. We observed these findings during the early stage of COVID-19, when infected ciliated cells were rapidly replaced by differentiating precursor cells. Moreover, our analyses revealed that SARS-CoV-2 cellular tropism was restricted to the nasal ciliated versus oral squamous epithelium. These results imply that targeting ciliated cells of the nasal epithelium during the early stage of COVID-19 could be an ideal strategy to prevent SARS-CoV-2 propagation.

Integrating molecular, histopathological, neuroimaging and clinical neuroscience data with NeuroPM-box

Understanding and treating heterogeneous brain disorders requires specialized techniques spanning genetics, proteomics, and neuroimaging. Designed to meet this need, NeuroPM-box is a user-friendly, open-access, multi-tool cross-platform software capable of characterizing multiscale and multifactorial neuropathological mechanisms. Using advanced analytical modeling for molecular, histopathological, brain-imaging and/or clinical evaluations, this framework has multiple applications, validated here with synthetic (N > 2900), in-vivo (N = 911) and post-mortem (N = 736) neurodegenerative data, and including the ability to characterize: (i) the series of sequential states (genetic, histopathological, imaging or clinical alterations) covering decades of disease progression, (ii) concurrent intra-brain spreading of pathological factors (e.g., amyloid, tau and alpha-synuclein proteins), (iii) synergistic interactions between multiple biological factors (e.g., toxic tau effects on brain atrophy), and (iv) biologically-defined patient stratification based on disease heterogeneity and/or therapeutic needs. This freely available toolbox ( neuropm-lab.com/neuropm-box.html ) could contribute significantly to a better understanding of complex brain processes and accelerating the implementation of Precision Medicine in Neurology.

Endocrine Pancreas Development and Dysfunction Through the Lens of Single-Cell RNA-Sequencing

A chronic inability to maintain blood glucose homeostasis leads to diabetes, which can damage multiple organs. The pancreatic islets regulate blood glucose levels through the coordinated action of islet cell-secreted hormones, with the insulin released by β-cells playing a crucial role in this process. Diabetes is caused by insufficient insulin secretion due to β-cell loss, or a pancreatic dysfunction. The restoration of a functional β-cell mass might, therefore, offer a cure. To this end, major efforts are underway to generate human β-cells de novo, in vitro, or in vivo. The efficient generation of functional β-cells requires a comprehensive knowledge of pancreas development, including the mechanisms driving cell fate decisions or endocrine cell maturation. Rapid progress in single-cell RNA sequencing (scRNA-Seq) technologies has brought a new dimension to pancreas development research. These methods can capture the transcriptomes of thousands of individual cells, including rare cell types, subtypes, and transient states. With such massive datasets, it is possible to infer the developmental trajectories of cell transitions and gene regulatory pathways. Here, we summarize recent advances in our understanding of endocrine pancreas development and function from scRNA-Seq studies on developing and adult pancreas and human endocrine differentiation models. We also discuss recent scRNA-Seq findings for the pathological pancreas in diabetes, and their implications for better treatment.

The Role of Cell Tracing and Fate Mapping Experiments in Cardiac Outflow Tract Development, New Opportunities through Emerging Technologies

Whilst knowledge regarding the pathophysiology of congenital heart disease (CHDs) has advanced greatly in recent years, the underlying developmental processes affecting the cardiac outflow tract (OFT) such as bicuspid aortic valve, tetralogy of Fallot and transposition of the great arteries remain poorly understood. Common among CHDs affecting the OFT, is a large variation in disease phenotypes. Even though the different cell lineages contributing to OFT development have been studied for many decades, it remains challenging to relate cell lineage dynamics to the morphologic variation observed in OFT pathologies. We postulate that the variation observed in cellular contribution in these congenital heart diseases might be related to underlying cell lineage dynamics of which little is known. We believe this gap in knowledge is mainly the result of technical limitations in experimental methods used for cell lineage analysis. The aim of this review is to provide an overview of historical fate mapping and cell tracing techniques used to study OFT development and introduce emerging technologies which provide new opportunities that will aid our understanding of the cellular dynamics underlying OFT pathology.