Methods, Challenges and Potentials of Single Cell RNA-seq

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Transcriptomic Approaches to Cardiomyocyte-Biomaterial Interactions: A Review

Biomaterials, essential for supporting, enhancing, and repairing damaged tissues, play a critical role in various medical applications. This Review focuses on the interaction of biomaterials and cardiomyocytes, emphasizing the unique significance of transcriptomic approaches in understanding their interactions, which are pivotal in cardiac bioengineering and regenerative medicine. Transcriptomic approaches serve as powerful tools to investigate how cardiomyocytes respond to biomaterials, shedding light on the gene expression patterns, regulatory pathways, and cellular processes involved in these interactions. Emerging technologies such as bulk RNA-seq, single-cell RNA-seq, single-nucleus RNA-seq, and spatial transcriptomics offer promising avenues for more precise and in-depth investigations. Longitudinal studies, pathway analyses, and machine learning techniques further improve the ability to explore the complex regulatory mechanisms involved. This review also discusses the challenges and opportunities of utilizing transcriptomic techniques in cardiomyocyte-biomaterial research. Although there are ongoing challenges such as costs, cell size limitation, sample differences, and complex analytical process, there exist exciting prospects in comprehensive gene expression analyses, biomaterial design, cardiac disease treatment, and drug testing. These multimodal methodologies have the capacity to deepen our understanding of the intricate interaction network between cardiomyocytes and biomaterials, potentially revolutionizing cardiac research with the aim of promoting heart health, and they are also promising for studying interactions between biomaterials and other cell types.

scVIC: deep generative modeling of heterogeneity for scRNA-seq data

Motivation

Single-cell RNA sequencing (scRNA-seq) has become a valuable tool for studying cellular heterogeneity. However, the analysis of scRNA-seq data is challenging because of inherent noise and technical variability. Existing methods often struggle to simultaneously explore heterogeneity across cells, handle dropout events, and account for batch effects. These drawbacks call for a robust and comprehensive method that can address these challenges and provide accurate insights into heterogeneity at the single-cell level.

Results

In this study, we introduce scVIC, an algorithm designed to account for variational inference, while simultaneously handling biological heterogeneity and batch effects at the single-cell level. scVIC explicitly models both biological heterogeneity and technical variability to learn cellular heterogeneity in a manner free from dropout events and the bias of batch effects. By leveraging variational inference, we provide a robust framework for inferring the parameters of scVIC. To test the performance of scVIC, we employed both simulated and biological scRNA-seq datasets, either including, or not, batch effects. scVIC was found to outperform other approaches because of its superior clustering ability and circumvention of the batch effects problem.

Availability and implementation

The code of scVIC and replication for this study are available at https://github.com/HiBearME/scVIC/tree/v1.0.

Research progress of single-cell transcriptome sequencing in autoimmune diseases and autoinflammatory disease: A review

Autoimmunity refers to the phenomenon that the body's immune system produces antibodies or sensitized lymphocytes to its own tissues to cause an immune response. Immune disorders caused by autoimmunity can mediate autoimmune diseases. Autoimmune diseases have complicated pathogenesis due to the many types of cells involved, and the mechanism is still unclear. The emergence of single-cell research technology can solve the problem that ordinary transcriptome technology cannot be accurate to cell type. It provides unbiased results through independent analysis of cells in tissues and provides more mRNA information for identifying cell subpopulations, which provides a novel approach to study disruption of immune tolerance and disturbance of pro-inflammatory pathways on a cellular basis. It may fundamentally change the understanding of molecular pathways in the pathogenesis of autoimmune diseases and develop targeted drugs. Single-cell transcriptome sequencing (scRNA-seq) has been widely applied in autoimmune diseases, which provides a powerful tool for demonstrating the cellular heterogeneity of tissues involved in various immune inflammations, identifying pathogenic cell populations, and revealing the mechanism of disease occurrence and development. This review describes the principles of scRNA-seq, introduces common sequencing platforms and practical procedures, and focuses on the progress of scRNA-seq in 41 autoimmune diseases, which include 9 systemic autoimmune diseases and autoinflammatory diseases (rheumatoid arthritis, systemic lupus erythematosus, etc.) and 32 organ-specific autoimmune diseases (5 Skin diseases, 3 Nervous system diseases, 4 Eye diseases, 2 Respiratory system diseases, 2 Circulatory system diseases, 6 Liver, Gallbladder and Pancreas diseases, 2 Gastrointestinal system diseases, 3 Muscle, Bones and joint diseases, 3 Urinary system diseases, 2 Reproductive system diseases). This review also prospects the molecular mechanism targets of autoimmune diseases from the multi-molecular level and multi-dimensional analysis combined with single-cell multi-omics sequencing technology (such as scRNA-seq, Single cell ATAC-seq and single cell immune group library sequencing), which provides a reference for further exploring the pathogenesis and marker screening of autoimmune diseases and autoimmune inflammatory diseases in the future.

Decoding brain memory formation by single-cell RNA sequencing

To understand how distinct memories are formed and stored in the brain is an important and fundamental question in neuroscience and computational biology. A population of neurons, termed engram cells, represents the physiological manifestation of a specific memory trace and is characterized by dynamic changes in gene expression, which in turn alters the synaptic connectivity and excitability of these cells. Recent applications of single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) are promising approaches for delineating the dynamic expression profiles in these subsets of neurons, and thus understanding memory-specific genes, their combinatorial patterns and regulatory networks. The aim of this article is to review and discuss the experimental and computational procedures of sc/snRNA-seq, new studies of molecular mechanisms of memory aided by sc/snRNA-seq in human brain diseases and related mouse models, and computational challenges in understanding the regulatory mechanisms underlying long-term memory formation.

Frontiers in single cell analysis: multimodal technologies and their clinical perspectives

Single cell multimodal analysis is at the frontier of single cell research: it defines the roles and functions of distinct cell types through simultaneous analysis to provide unprecedented insight into cellular processes. Current single cell approaches are rapidly moving toward multimodal characterizations. It replaces one-dimensional single cell analysis, for example by allowing for simultaneous measurement of transcription and post-transcriptional regulation, epigenetic modifications and/or surface protein expression. By providing deeper insights into single cell processes, multimodal single cell analyses paves the way to new understandings in various cellular processes such as cell fate decisions, physiological heterogeneity or genotype-phenotype linkages. At the forefront of this, microfluidics is key for high-throughput single cell analysis. Here, we present an overview of the recent multimodal microfluidic platforms having a potential in biomedical research, with a specific focus on their potential clinical applications.

An optimized protocol for single nuclei isolation from clinical biopsies for RNA-seq

Single nuclei RNA sequencing (snRNA-seq) has evolved as a powerful tool to study complex human diseases. Single cell resolution enables the study of novel cell types, biological processes, cell trajectories, and cell-cell signaling pathways. snRNA-seq largely relies on the dissociation of intact nuclei from human tissues. However, the study of complex tissues using small core biopsies presents many technical challenges. Here, an optimized protocol for single nuclei isolation is presented for frozen and RNAlater preserved human kidney biopsies. The described protocol is fast, low cost, and time effective due to the elimination of cell sorting and ultra-centrifugation. Samples can be processed in 90 min or less. This method is effective for obtaining normal nuclei morphology without signs of structural damage. Using snRNA-seq, 16 distinct kidney cell clusters were recovered from normal and peri-transplant acute kidney injury allograft samples, including immune cell clusters. Quality control measurements demonstrated that these optimizations eliminated cellular debris and allowed for a high yield of high-quality nuclei and RNA for library preparation and sequencing. Cellular disassociation did not induce cellular stress responses, which recapitulated transcriptional patterns associated with standardized methods of nuclei isolation. Future applications of this protocol will allow for thorough investigations of small biobank biopsies, identifying cell-specific injury pathways and driving the discovery of novel diagnostics and therapeutic targets.

Anti-bias training for (sc)RNA-seq: experimental and computational approaches to improve precision

RNA-seq, including single cell RNA-seq (scRNA-seq), is plagued by insufficient sensitivity and lack of precision. As a result, the full potential of (sc)RNA-seq is limited. Major factors in this respect are the presence of global bias in most datasets, which affects detection and quantitation of RNA in a length-dependent fashion. In particular, scRNA-seq is affected by technical noise and a high rate of dropouts, where the vast majority of original transcripts is not converted into sequencing reads. We discuss these biases origins and implications, bioinformatics approaches to correct for them, and how biases can be exploited to infer characteristics of the sample preparation process, which in turn can be used to improve library preparation.

Imputation of spatially-resolved transcriptomes by graph-regularized tensor completion

High-throughput spatial-transcriptomics RNA sequencing (sptRNA-seq) based on in-situ capturing technologies has recently been developed to spatially resolve transcriptome-wide mRNA expressions mapped to the captured locations in a tissue sample. Due to the low RNA capture efficiency by in-situ capturing and the complication of tissue section preparation, sptRNA-seq data often only provides an incomplete profiling of the gene expressions over the spatial regions of the tissue. In this paper, we introduce a graph-regularized tensor completion model for imputing the missing mRNA expressions in sptRNA-seq data, namely FIST, Fast Imputation of Spatially-resolved transcriptomes by graph-regularized Tensor completion. We first model sptRNA-seq data as a 3-way sparse tensor in genes (p-mode) and the (x, y) spatial coordinates (x-mode and y-mode) of the observed gene expressions, and then consider the imputation of the unobserved entries or fibers as a tensor completion problem in Canonical Polyadic Decomposition (CPD) form. To improve the imputation of highly sparse sptRNA-seq data, we also introduce a protein-protein interaction network to add prior knowledge of gene functions, and a spatial graph to capture the the spatial relations among the capture spots. The tensor completion model is then regularized by a Cartesian product graph of protein-protein interaction network and the spatial graph to capture the high-order relations in the tensor. In the experiments, FIST was tested on ten 10x Genomics Visium spatial transcriptomic datasets of different tissue sections with cross-validation among the known entries in the imputation. FIST significantly outperformed the state-of-the-art methods for single-cell RNAseq data imputation. We also demonstrate that both the spatial graph and PPI network play an important role in improving the imputation. In a case study, we further analyzed the gene clusters obtained from the imputed gene expressions to show that the imputations by FIST indeed capture the spatial characteristics in the gene expressions and reveal functions that are highly relevant to three different kinds of tissues in mouse kidney.

Biological tissues are composed of different types of structurally organized cell units playing distinct functional roles. The exciting new spatial gene expression profiling methods have enabled the analysis of spatially resolved transcriptomes to understand the spatial and functional characteristics of these cells in the context of eco-environment of tissue. Due to the technical limitations, spatial transcriptomics data suffers from only sparsely measured mRNAs by in-situ capture and possibly missing spots in tissue regions that entirely failed fixing and permeabilizing RNAs. Our method, FIST (Fast Imputation of Spatially-resolved transcriptomes by graph-regularized Tensor completion), focuses on the spatial and high-sparsity nature of spatial transcriptomics data by modeling the data as a 3-way gene-by-(x, y)-location tensor and a product graph of a spatial graph and a protein-protein interaction network. Our comprehensive evaluation of FIST on ten 10x Genomics Visium spatial genomics datasets and comparison with the methods for single-cell RNA sequencing data imputation demonstrate that FIST is a better method more suitable for spatial gene expression imputation. Overall, we found FIST a useful new method for analyzing spatially resolved gene expressions based on novel modeling of spatial and functional information.

Deep soft K-means clustering with self-training for single-cell RNA sequence data

Single-cell RNA sequencing (scRNA-seq) allows researchers to study cell heterogeneity at the cellular level. A crucial step in analyzing scRNA-seq data is to cluster cells into subpopulations to facilitate subsequent downstream analysis. However, frequent dropout events and increasing size of scRNA-seq data make clustering such high-dimensional, sparse and massive transcriptional expression profiles challenging. Although some existing deep learning-based clustering algorithms for single cells combine dimensionality reduction with clustering, they either ignore the distance and affinity constraints between similar cells or make some additional latent space assumptions like mixture Gaussian distribution, failing to learn cluster-friendly low-dimensional space. Therefore, in this paper, we combine the deep learning technique with the use of a denoising autoencoder to characterize scRNA-seq data while propose a soft self-training K-means algorithm to cluster the cell population in the learned latent space. The self-training procedure can effectively aggregate the similar cells and pursue more cluster-friendly latent space. Our method, called ‘scziDesk’, alternately performs data compression, data reconstruction and soft clustering iteratively, and the results exhibit excellent compatibility and robustness in both simulated and real data. Moreover, our proposed method has perfect scalability in line with cell size on large-scale datasets.

Single-Cell Transcriptome Data Clustering via Multinomial Modeling and Adaptive Fuzzy K-Means Algorithm

Single-cell RNA sequencing technologies have enabled us to study tissue heterogeneity at cellular resolution. Fast-developing sequencing platforms like droplet-based sequencing make it feasible to parallel process thousands of single cells effectively. Although a unique molecular identifier (UMI) can remove bias from amplification noise to a certain extent, clustering for such sparse and high-dimensional large-scale discrete data remains intractable and challenging. Most existing deep learning-based clustering methods utilize the mean square error or negative binomial distribution with or without zero inflation to denoise single-cell UMI count data, which may underfit or overfit the gene expression profiles. In addition, neglecting the molecule sampling mechanism and extracting representation by simple linear dimension reduction with a hard clustering algorithm may distort data structure and lead to spurious analytical results. In this paper, we combined the deep autoencoder technique with statistical modeling and developed a novel and effective clustering method, scDMFK, for single-cell transcriptome UMI count data. ScDMFK utilizes multinomial distribution to characterize data structure and draw support from neural network to facilitate model parameter estimation. In the learned low-dimensional latent space, we proposed an adaptive fuzzy k-means algorithm with entropy regularization to perform soft clustering. Various simulation scenarios and the analysis of 10 real datasets have shown that scDMFK outperforms other state-of-the-art methods with respect to data modeling and clustering algorithms. Besides, scDMFK has excellent scalability for large-scale single-cell datasets.

A Branch Point on Differentiation Trajectory is the Bifurcating Event Revealed by Dynamical Network Biomarker Analysis of Single-Cell Data

The advance in single-cell profiling technologies and the development in computational algorithms provide the opportunity to reconstruct pseudo temporal trajectory with branch point of cellular development. On the other hand, theories such as dynamical network biomarkers (DNB) theory have been recently proposed to characterize the pre-transition state in biological systems. Few studies have validated whether the branch point identified in pseudo time is the critical point in dynamical system. In this study, the dynamical behavior of the branch point on the pseudo trajectory has been investigated. We study the pseudo temporal trajectories reconstructed by Wishbone and diffusion pseudotime analysis (DPT) algorithms, as well as the simulated trajectory. DNB theory is applied to justify the bifurcating event on the pseudo trajectories. Our results demonstrate that the branch point recovered by Wishbone and DPT algorithms is confirmed as a transition state in cell differentiation process by DNB theory. Furthermore, we show that an appropriate DNB group will amplify the comprehensive index of critical event as defined in DNB theory. Our study provides biological insights on pseudo trajectory with branch point in a dynamical view and also indicates that DNB theory may serve as a benchmark to check the validity of branch point.

Normalization Methods on Single-Cell RNA-seq Data: An Empirical Survey

Data normalization is vital to single-cell sequencing, addressing limitations presented by low input material and various forms of bias or noise present in the sequencing process. Several such normalization methods exist, some of which rely on spike-in genes, molecules added in known quantities to serve as a basis for a normalization model. Depending on available information and the type of data, some methods may express certain advantages over others. We compare the effectiveness of seven available normalization methods designed specifically for single-cell sequencing using two real data sets containing spike-in genes and one simulation study. Additionally, we test those methods not dependent on spike-in genes using a real data set with three distinct cell-cycle states and a real data set under the 10X Genomics GemCode platform with multiple cell types represented. We demonstrate the differences in effectiveness for the featured methods using visualization and classification assessment and conclude which methods are preferable for normalizing a certain type of data for further downstream analysis, such as classification or differential analysis. The comparison in computational time for all methods is addressed as well.

A Single Cell but Many Different Transcripts: A Journey into the World of Long Non-Coding RNAs

In late 2012 it was evidenced that most of the human genome is transcribed but only a small percentage of the transcripts are translated. This observation supported the importance of non-coding RNAs and it was confirmed in several organisms. The most abundant non-translated transcripts are long non-coding RNAs (lncRNAs). In contrast to protein-coding RNAs, they show a more cell-specific expression. To understand the function of lncRNAs, it is fundamental to investigate in which cells they are preferentially expressed and to detect their subcellular localization. Recent improvements of techniques that localize single RNA molecules in tissues like single-cell RNA sequencing and fluorescence amplification methods have given a considerable boost in the knowledge of the lncRNA functions. In recent years, single-cell transcription variability was associated with non-coding RNA expression, revealing this class of RNAs as important transcripts in the cell lineage specification. The purpose of this review is to collect updated information about lncRNA classification and new findings on their function derived from single-cell analysis. We also retained useful for all researchers to describe the methods available for single-cell analysis and the databases collecting single-cell and lncRNA data. Tables are included to schematize, describe, and compare exposed concepts.

Early Transcriptional Landscapes of Chlamydia trachomatis-Infected Epithelial Cells at Single Cell Resolution

Chlamydia are Gram-negative obligate intracellular bacterial pathogens responsible for a variety of disease in humans and animals worldwide. Chlamydia trachomatis causes trachoma in disadvantaged populations, and is the most common bacterial sexually transmitted infection in humans, causing reproductive tract disease. Antibiotic therapy successfully treats diagnosed chlamydial infections, however asymptomatic infections are common. High-throughput transcriptomic approaches have explored chlamydial gene expression and infected host cell gene expression. However, these were performed on large cell populations, averaging gene expression profiles across all cells sampled and potentially obscuring biologically relevant subsets of cells. We generated a pilot dataset, applying single cell RNA-Seq (scRNA-Seq) to C. trachomatis infected and mock-infected epithelial cells to assess the utility, pitfalls and challenges of single cell approaches applied to chlamydial biology, and to potentially identify early host cell biomarkers of chlamydial infection. Two hundred sixty-four time-matched C. trachomatis-infected and mock-infected HEp-2 cells were collected and subjected to scRNA-Seq. After quality control, 200 cells were retained for analysis. Two distinct clusters distinguished 3-h cells from 6- and 12-h. Pseudotime analysis identified a possible infection-specific cellular trajectory for Chlamydia-infected cells, while differential expression analyses found temporal expression of metallothioneins and genes involved with cell cycle regulation, innate immune responses, cytoskeletal components, lipid biosynthesis and cellular stress. We find that changes to the host cell transcriptome at early times of C. trachomatis infection are readily discernible by scRNA-Seq, supporting the utility of single cell approaches to identify host cell biomarkers of chlamydial infection, and to further deconvolute the complex host response to infection.

Machine learning and statistical methods for clustering single-cell RNA-sequencing data

 

Single-cell RNAsequencing (scRNA-seq) technologies have enabled the large-scale whole-transcriptome profiling of each individual single cell in a cell population. A core analysis of the scRNA-seq transcriptome profiles is to cluster the single cells to reveal cell subtypes and infer cell lineages based on the relations among the cells. This article reviews the machine learning and statistical methods for clustering scRNA-seq transcriptomes developed in the past few years. The review focuses on how conventional clustering techniques such as hierarchical clustering, graph-based clustering, mixture models, $k$-means, ensemble learning, neural networks and density-based clustering are modified or customized to tackle the unique challenges in scRNA-seq data analysis, such as the dropout of low-expression genes, low and uneven read coverage of transcripts, highly variable total mRNAs from single cells and ambiguous cell markers in the presence of technical biases and irrelevant confounding biological variations. We review how cell-specific normalization, the imputation of dropouts and dimension reduction methods can be applied with new statistical or optimization strategies to improve the clustering of single cells. We will also introduce those more advanced approaches to cluster scRNA-seq transcriptomes in time series data and multiple cell populations and to detect rare cell types. Several software packages developed to support the cluster analysis of scRNA-seq data are also reviewed and experimentally compared to evaluate their performance and efficiency. Finally, we conclude with useful observations and possible future directions in scRNA-seq data analytics.

Availability

All the source code and data are available at https://github.com/kuanglab/single-cell-review.

Predictive approaches for drug combination discovery in cancer

Drug combinations have been proposed as a promising therapeutic strategy to overcome drug resistance and improve efficacy of monotherapy regimens in cancer. This strategy aims at targeting multiple components of this complex disease. Despite the increasing number of drug combinations in use, many of them were empirically found in the clinic, and the molecular mechanisms underlying these drug combinations are often unclear. These challenges call for rational, systematic approaches for drug combination discovery. Although high-throughput screening of single-agent therapeutics has been successfully implemented, it is not feasible to test all possible drug combinations, even for a reduced subset of anticancer drugs. Hence, in vitro and in vivo screening of a large number of drug combinations are not practical. Therefore, devising computational methods to efficiently explore the space of drug combinations and to discover efficacious combinations has attracted a lot of attention from the scientific community in the past few years. Nevertheless, in the absence of consensus regarding the computational approaches used to predict efficacious drug combinations, a plethora of methods, techniques and hypotheses have been developed to date, while the research field lacks an elaborate categorization of the existing computational methods and the available data sources. In this manuscript, we review and categorize the state-of-the-art computational approaches for drug combination prediction, and elaborate on the limitations of these methods and the existing challenges. We also discuss about the recent pan-cancer drug combination data sets and their importance in revising the available methods or developing more performant approaches.

Technical Advances in Single-Cell RNA Sequencing and Applications in Normal and Malignant Hematopoiesis

Single-cell RNA sequencing (scRNA-seq) has been tremendously developed in the past decade owing to overcoming challenges associated with isolation of massive quantities of single cells. Previously, cell heterogeneity and low quantities of available biological material posed significant difficulties to scRNA-seq. Cell-to-cell variation and heterogeneity are fundamental and intrinsic characteristics of normal and malignant hematopoietic cells; this heterogeneity has often been ignored in omics studies. The application of scRNA-seq has profoundly changed our comprehension of many biological phenomena, including organ development and carcinogenesis. Hematopoiesis, is actually a maturation process for more than ten distinct blood and immune cells, and is thought to be critically involved in hematological homeostasis and in sustaining the physiological functions. However, aberrant hematopoiesis directly leads to hematological malignancy, and a deeper understanding of malignant hematopoiesis will provide deeper insights into diagnosis and prognosis for patients with hematological malignancies. Here, we aim to review the recent technical progress and future prospects for scRNA-seq, as applied in physiological and malignant hematopoiesis, in efforts to further understand the hematopoietic hierarchy and to illuminate personalized therapy and precision medicine approaches used in the clinical treatment of hematological malignancies.

A multitask clustering approach for single-cell RNA-seq analysis in Recessive Dystrophic Epidermolysis Bullosa

Single-cell RNA sequencing (scRNA-seq) has been widely applied to discover new cell types by detecting sub-populations in a heterogeneous group of cells. Since scRNA-seq experiments have lower read coverage/tag counts and introduce more technical biases compared to bulk RNA-seq experiments, the limited number of sampled cells combined with the experimental biases and other dataset specific variations presents a challenge to cross-dataset analysis and discovery of relevant biological variations across multiple cell populations. In this paper, we introduce a method of variance-driven multitask clustering of single-cell RNA-seq data (scVDMC) that utilizes multiple single-cell populations from biological replicates or different samples. scVDMC clusters single cells in multiple scRNA-seq experiments of similar cell types and markers but varying expression patterns such that the scRNA-seq data are better integrated than typical pooled analyses which only increase the sample size. By controlling the variance among the cell clusters within each dataset and across all the datasets, scVDMC detects cell sub-populations in each individual experiment with shared cell-type markers but varying cluster centers among all the experiments. Applied to two real scRNA-seq datasets with several replicates and one large-scale droplet-based dataset on three patient samples, scVDMC more accurately detected cell populations and known cell markers than pooled clustering and other recently proposed scRNA-seq clustering methods. In the case study applied to in-house Recessive Dystrophic Epidermolysis Bullosa (RDEB) scRNA-seq data, scVDMC revealed several new cell types and unknown markers validated by flow cytometry. MATLAB/Octave code available at https://github.com/kuanglab/scVDMC.

scRNA-seq enables detailed profiling of heterogeneous cell populations and can be used to reveal lineage relationships or discover new cell types. In the literature, there has been little effort directed towards developing computational methods for cross-population transcriptome analysis of multiple single-cell populations. The cross-cell-population clustering problem is different from the traditional clustering problem because single-cell populations can be collected from different patients, different samples of a tissue, or different experimental replicates. The accompanying biological and technical variation tends to dominate the signals for clustering the pooled single cells from the multiple populations. In this work, we have developed a multitask clustering method to address the cross-population clustering problem. The method simultaneously clusters each individual cell population and controls variance among the cell-type cluster centers within each cell population and across the cell populations. We demonstrate that our multitask clustering method significantly improves clustering accuracy and marker discovery in three public scRNA-seq datasets and also apply the method to an in-house Recessive Dystrophic Epidermolysis Bullosa (RDEB) dataset. Our results make it evident that multitask clustering is a promising new approach for cross-population analysis of scRNA-seq data.

How to design a single-cell RNA-sequencing experiment: pitfalls, challenges and perspectives

The sequencing of the transcriptome of single cells, or single-cell RNA-sequencing, has now become the dominant technology for the identification of novel cell types in heterogeneous cell populations or for the study of stochastic gene expression. In recent years, various experimental methods and computational tools for analysing single-cell RNA-sequencing data have been proposed. However, most of them are tailored to different experimental designs or biological questions, and in many cases, their performance has not been benchmarked yet, thus increasing the difficulty for a researcher to choose the optimal single-cell transcriptome sequencing (scRNA-seq) experiment and analysis workflow. In this review, we aim to provide an overview of the current available experimental and computational methods developed to handle single-cell RNA-sequencing data and, based on their peculiarities, we suggest possible analysis frameworks depending on specific experimental designs. Together, we propose an evaluation of challenges and open questions and future perspectives in the field. In particular, we go through the different steps of scRNA-seq experimental protocols such as cell isolation, messenger RNA capture, reverse transcription, amplification and use of quantitative standards such as spike-ins and Unique Molecular Identifiers (UMIs). We then analyse the current methodological challenges related to preprocessing, alignment, quantification, normalization, batch effect correction and methods to control for confounding effects.

Spatial transcriptomic analysis of cryosectioned tissue samples with Geo-seq

Conventional gene expression studies analyze multiple cells simultaneously or single cells, for which the exact in vivo or in situ position is unknown. Although cellular heterogeneity can be discerned when analyzing single cells, any spatially defined attributes that underpin the heterogeneous nature of the cells cannot be identified. Here, we describe how to use Geo-seq, a method that combines laser capture microdissection (LCM) and single-cell RNA-seq technology. The combination of these two methods enables the elucidation of cellular heterogeneity and spatial variance simultaneously. The Geo-seq protocol allows the profiling of transcriptome information from only a small number cells and retains their native spatial information. This protocol has wide potential applications to address biological and pathological questions of cellular properties such as prospective cell fates, biological function and the gene regulatory network. Geo-seq has been applied to investigate the spatial transcriptome of mouse early embryo, mouse brain, and pathological liver and sperm tissues. The entire protocol from tissue collection and microdissection to sequencing requires ∼5 d, Data analysis takes another 1 or 2 weeks, depending on the amount of data and the speed of the processor.

Modeling Enzyme Processivity Reveals that RNA-Seq Libraries Are Biased in Characteristic and Correctable Ways

Experimental procedures for preparing RNA-seq and single-cell (sc) RNA-seq libraries are based on assumptions regarding their underlying enzymatic reactions. Here, we show that the fairness of these assumptions varies within libraries: coverage by sequencing reads along and between transcripts exhibits characteristic, protocol-dependent biases. To understand the mechanistic basis of this bias, we present an integrated modeling framework that infers the relationship between enzyme reactions during library preparation and the characteristic coverage patterns observed for different protocols. Analysis of new and existing (sc)RNA-seq data from six different library preparation protocols reveals that polymerase processivity is the mechanistic origin of coverage biases. We apply our framework to demonstrate that lowering incubation temperature increases processivity, yield, and (sc)RNA-seq sensitivity in all protocols. We also provide correction factors based on our model for increasing accuracy of transcript quantification in existing samples prepared at standard temperatures. In total, our findings improve our ability to accurately reflect in vivo transcript abundances in (sc)RNA-seq libraries.

•

Characterization of global RNA-seq biases specific to library preparation protocols

•

Mathematical framework to reverse engineer enzyme reactions that cause bias

•

Insights from reverse engineering allow optimization of RNA-seq protocols

•

Lowered incubation temperatures during library preparation improve sensitivity

Single Cell Isolation and Analysis

Individual cell heterogeneity within a population can be critical to its peculiar function and fate. Subpopulations studies with mixed mutants and wild types may not be as informative regarding which cell responds to which drugs or clinical treatments. Cell to cell differences in RNA transcripts and protein expression can be key to answering questions in cancer, neurobiology, stem cell biology, immunology, and developmental biology. Conventional cell-based assays mainly analyze the average responses from a population of cells, without regarding individual cell phenotypes. To better understand the variations from cell to cell, scientists need to use single cell analyses to provide more detailed information for therapeutic decision making in precision medicine. In this review, we focus on the recent developments in single cell isolation and analysis, which include technologies, analyses and main applications. Here, we summarize the historical background, limitations, applications, and potential of single cell isolation technologies.

Undiscovered Physiology of Transcript and Protein Networks

The past two decades have witnessed a rapid evolution in our ability to measure RNA and protein from biological systems. As a result, new principles have arisen regarding how information is processed in cells, how decisions are made, and the role of networks in biology. This essay examines this technological evolution, reviewing (and critiquing) the conceptual framework that has emerged to explain how RNA and protein networks control cellular function. We identify how future investigations into transcriptomes, proteomes, and other cellular networks will enable development of more robust, quantitative models of cellular behavior whilst also providing new avenues to use knowledge of biological networks to improve human health. © 2016 American Physiological Society. Compr Physiol 6:1851-1872, 2016.

Single-cell Transcriptome Study as Big Data

The rapid growth of single-cell RNA-seq studies (scRNA-seq) demands efficient data storage, processing, and analysis. Big-data technology provides a framework that facilitates the comprehensive discovery of biological signals from inter-institutional scRNA-seq datasets. The strategies to solve the stochastic and heterogeneous single-cell transcriptome signal are discussed in this article. After extensively reviewing the available big-data applications of next-generation sequencing (NGS)-based studies, we propose a workflow that accounts for the unique characteristics of scRNA-seq data and primary objectives of single-cell studies.

Implications of Epigenetic Variability within a Cell Population for "Cell Type" Classification

Here, we propose a new approach to defining nerve “cell types” in reaction to recent advances in single cell analysis. Among cells previously thought to be equivalent, considerable differences in global gene expression and biased tendencies among differing developmental fates have been demonstrated within multiple lineages. The model of classifying cells into distinct types thus has to be revised to account for this intrinsic variability. A “cell type” could be a group of cells that possess similar, but not necessarily identical properties, variable within a spectrum of epigenetic adjustments that permit its developmental path toward a specific function to be achieved. Thus, the definition of a cell type is becoming more similar to the definition of a species: sharing essential properties with other members of its group, but permitting a certain amount of deviation in aspects that do not seriously impact function. This approach accommodates, even embraces the spectrum of natural variation found in various cell populations and consequently avoids the fallacy of false equivalence. For example, developing neurons will react to their microenvironments with epigenetic changes resulting in slight changes in gene expression and morphology. Addressing the new questions implied here will have significant implications for developmental neurobiology.

Recent advances and current issues in single-cell sequencing of tumors

Intratumoral heterogeneity is a recently recognized but important feature of cancer that underlies the various biocharacteristics of cancer tissues. The advent of next-generation sequencing technologies has facilitated large scale capture of genomic data, while the recent development of single-cell sequencing has allowed for more in-depth studies into the complex molecular mechanisms of intratumoral heterogeneity. In this review, the recent advances and current challenges in single-cell sequencing methodologies are discussed, highlighting the potential power of these data to provide insights into oncological processes, from tumorigenesis through progression to metastasis and therapy resistance.

Parallel RNA extraction using magnetic beads and a droplet array

Nucleic acid extraction is a necessary step for most genomic/transcriptomic analyses, but it often requires complicated mechanisms to be integrated into a lab-on-a-chip device. Here, we present a simple, effective configuration for rapidly obtaining purified RNA from low concentration cell medium. This Total RNA Extraction Droplet Array (TREDA) utilizes an array of surface-adhering droplets to facilitate the transportation of magnetic purification beads seamlessly through individual buffer solutions without solid structures. The fabrication of TREDA chips is rapid and does not require a microfabrication facility or expertise. The process takes less than 5 minutes. When purifying mRNA from bulk marine diatom samples, its repeatability and extraction efficiency are comparable to conventional tube-based operations. We demonstrate that TREDA can extract the total mRNA of about 10 marine diatom cells, indicating that the sensitivity of TREDA approaches single-digit cell numbers.

A new toolbox for assessing single cells

Unprecedented access to the biology of single cells is now feasible, enabled by recent technological advancements that allow us to manipulate and measure sparse samples and achieve a new level of resolution in space and time. This review focuses on advances in tools to study single cells for specific areas of biology. We examine both mature and nascent techniques to study single cells at the genomics, transcriptomics, and proteomics level. In addition, we provide an overview of tools that are well suited for following biological responses to defined perturbations with single-cell resolution. Techniques to analyze and manipulate single cells through soluble and chemical ligands, the microenvironment, and cell-cell interactions are provided. For each of these topics, we highlight the biological motivation, applications, methods, recent advances, and opportunities for improvement. The toolbox presented in this review can function as a starting point for the design of single-cell experiments.

Single-cell analysis of the transcriptome and its application in the characterization of stem cells and early embryos

Cellular heterogeneity within a cell population is a common phenomenon in multicellular organisms, tissues, cultured cells, and even FACS-sorted subpopulations. Important information may be masked if the cells are studied as a mass. Transcriptome profiling is a parameter that has been intensively studied, and relatively easier to address than protein composition. To understand the basis and importance of heterogeneity and stochastic aspects of the cell function and its mechanisms, it is essential to examine transcriptomes of a panel of single cells. High-throughput technologies, starting from microarrays and now RNA-seq, provide a full view of the expression of transcriptomes but are limited by the amount of RNA for analysis. Recently, several new approaches for amplification and sequencing the transcriptome of single cells or a limited low number of cells have been developed and applied. In this review, we summarize these major strategies, such as PCR-based methods, IVT-based methods, phi29-DNA polymerase-based methods, and several other methods, including their principles, characteristics, advantages, and limitations, with representative applications in cancer stem cells, early development, and embryonic stem cells. The prospects for development of future technology and application of transcriptome analysis in a single cell are also discussed.

Current challenges in the bioinformatics of single cell genomics

Single cell genomics is a rapidly growing field with many new techniques emerging in the past few years. However, few bioinformatics tools specific for single cell genomics analysis are available. Single cell DNA/RNA sequencing data usually have low genome coverage and high amplification bias, which makes bioinformatics analysis challenging. Many current bioinformatics tools developed for bulk cell sequencing do not work well with single cell sequencing data. Here, we summarize current challenges in the bioinformatics analysis of single cell genomic DNA sequencing and single cell transcriptomes. These challenges include calling copy number variations, identifying mutated genes in tumor samples, reconstructing cell lineages, recovering low abundant transcripts, and improving the accuracy of quantitative analysis of transcripts. Development in single cell genomics bioinformatics analysis will promote the application of this technology to basic biology and medical research.

Laser-assisted microdissection applied to floral tissues

Cellular context can be crucial when studying developmental processes as well as responses to environmental variation. Several different tools have been developed in recent years to isolate specific tissues or cell types. Laser-assisted microdissection (LAM) allows for the isolation of such specific tissue or single cell-types purely based on morphology and cytology. This has the advantage that (1) cell types that are rare can be isolated from heterogeneous tissue, (2) no marker line with cell type-specific expression needs to be established, and (3) the method can be applied to non-model species and species that are difficult to genetically transform. The rapid development of next-generation sequencing (NGS) approaches has greatly advanced the possibilities to perform molecular analyses in diverse organisms. However, there is a mismatch between currently available cell isolation tools and their applicability to non-model organisms. Therefore, LAM will become increasingly popular in the study of diverse agriculturally or ecologically relevant plant species. Here, we describe a protocol that has been successfully used for LAM to isolate either whole floral organs or even single cell types in plants, e.g., Arabidopsis thaliana, Boechera spp., or tomato.

High-throughput, multiparameter analysis of single cells

Heterogeneity of cell populations in various biological systems has been widely recognized, and the highly heterogeneous nature of cancer cells has been emerging with clinical relevance. Single-cell analysis using a combination of high-throughput and multiparameter approaches is capable of reflecting cell-to-cell variability, and at the same time of unraveling the complexity and interdependence of cellular processes in the individual cells of a heterogeneous population. In this review, analytical methods and microfluidic tools commonly used for high-throughput, multiparameter single-cell analysis of DNA, RNA, and proteins are discussed. Applications and limitations of currently available technologies for cancer research and diagnostics are reviewed in the light of the ultimate goal to establish clinically applicable assays.

Genomics and proteomics in solving brain complexity

The human brain is extraordinarily complex, composed of billions of neurons and trillions of synaptic connections. Neurons are organized into circuit assemblies that are modulated by specific interneurons and non-neuronal cells, such as glia and astrocytes. Data on human genome sequences predicts that each of these cells in the human brain has the potential of expressing ∼20 000 protein coding genes and tens of thousands of noncoding RNAs. A major challenge in neuroscience is to determine (1) how individual neurons and circuitry utilize this potential during development and maturation of the nervous system, and for higher brain functions such as cognition, and (2) how this potential is altered in neurological and psychiatric disorders. In this review, we will discuss how recent advances in next generation sequencing, proteomics and bioinformatics have transformed our understanding of gene expression and the functions of neural circuitry, memory storage, and disorders of cognition.

The role of single-cell analyses in understanding cell lineage commitment

The study of cell lineage commitment is critical for improving our understanding of tissue development and regeneration, and for realizing stem cell-based therapies and engineered tissue replacements. Recently, the discovery of an unanticipated degree of variability in fundamental biological processes, including divergent responses of genetically identical cells to various stimuli, has provided mechanistic insight into cellular decision making and the collective behavior of cell populations. Therefore, the study of lineage commitment with single-cell resolution could provide greater knowledge of cellular differentiation mechanisms and the influence of noise on cellular processes. This will require the adoption of new technologies for single-cell analysis as traditional methods typically measure average values of bulk population behavior. This review discusses the recent developments in methods for analyzing the behavior of individual cells, and how these approaches are leading to a deeper understanding and better control of cellular decision making.