Cell population-based framework of genetic epidemiology in the single-cell omics era

Genetic epidemiology is a rapidly advancing field due to the recent availability of large amounts of omics data. In recent years, it has become possible to obtain omics information at the single-cell level, so genetic epidemiological models need to be updated to integrate with single-cell expression data. In this perspective paper, we propose a cell population-based framework for genetic epidemiology in the single-cell era. In this framework, genetic diversity influences phenotypic diversity through the diversity of cell population profiles, which are defined as high-dimensional probability distributions of the state spaces of biomolecules of each omics layer. We discuss how biomolecular experimental measurement data can capture the different properties of this distribution. In particular, single-cell data constitute a sample from this population distribution where only some coordinate values are observable. From a data analysis standpoint, we introduce methodology for feature extraction from cell population profiles. Finally, we discuss how this framework can be applied not only to genetic epidemiology but also to systems biology.     

Scalable and flexible inference framework for stochastic dynamic single-cell models

Understanding the inherited nature of how biological processes dynamically change over time and exhibit intra- and inter-individual variability, due to the different responses to environmental stimuli and when interacting with other processes, has been a major focus of systems biology. The rise of single-cell fluorescent microscopy has enabled the study of those phenomena. The analysis of single-cell data with mechanistic models offers an invaluable tool to describe dynamic cellular processes and to rationalise cell-to-cell variability within the population. However, extracting mechanistic information from single-cell data has proven difficult. This requires statistical methods to infer unknown model parameters from dynamic, multi-individual data accounting for heterogeneity caused by both intrinsic (e.g. variations in chemical reactions) and extrinsic (e.g. variability in protein concentrations) noise. Although several inference methods exist, the availability of efficient, general and accessible methods that facilitate modelling of single-cell data, remains lacking. Here we present a scalable and flexible framework for Bayesian inference in state-space mixed-effects single-cell models with stochastic dynamic. Our approach infers model parameters when intrinsic noise is modelled by either exact or approximate stochastic simulators, and when extrinsic noise is modelled by either time-varying, or time-constant parameters that vary between cells. We demonstrate the relevance of our approach by studying how cell-to-cell variation in carbon source utilisation affects heterogeneity in the budding yeast Saccharomyces cerevisiae SNF1 nutrient sensing pathway. We identify hexokinase activity as a source of extrinsic noise and deduce that sugar availability dictates cell-to-cell variability.     

Integrating multiple references for single-cell assignment

Efficient single-cell assignment is essential for single-cell sequencing data analysis. With the explosive growth of single-cell sequencing data, multiple single-cell sequencing data sources are available for the same kind of tissue, which can be integrated to further improve single-cell assignment; however, an efficient integration strategy is still lacking due to the great challenges of data heterogeneity existing in multiple references. To this end, we present mtSC, a flexible single-cell assignment framework that integrates multiple references based on multitask deep metric learning designed specifically for cell type identification within tissues with multiple single-cell sequencing data as references. We evaluated mtSC on a comprehensive set of publicly available benchmark datasets and demonstrated its state-of-the-art effectiveness for integrative single-cell assignment with multiple references.     

Quantifying Stochastic Noise in Cultured Circadian Reporter Cells

Stochastic noise at the cellular level has been shown to play a fundamental role in circadian oscillations, influencing how groups of cells entrain to external cues and likely serving as the mechanism by which cell-autonomous rhythms are generated. Despite this importance, few studies have investigated how clock perturbations affect stochastic noise—even as increasing numbers of high-throughput screens categorize how gene knockdowns or small molecules can change clock period and amplitude. This absence is likely due to the difficulty associated with measuring cell-autonomous stochastic noise directly, which currently requires the careful collection and processing of single-cell data. In this study, we show that the damping rate of population-level bioluminescence recordings can serve as an accurate measure of overall stochastic noise, and one that can be applied to future and existing high-throughput circadian screens. Using cell-autonomous fibroblast data, we first show directly that higher noise at the single-cell results in faster damping at the population level. Next, we show that the damping rate of cultured cells can be changed in a dose-dependent fashion by small molecule modulators, and confirm that such a change can be explained by single-cell noise using a mathematical model. We further demonstrate the insights that can be gained by applying our method to a genome-wide siRNA screen, revealing that stochastic noise is altered independently from period, amplitude, and phase. Finally, we hypothesize that the unperturbed clock is highly optimized for robust rhythms, as very few gene perturbations are capable of simultaneously increasing amplitude and lowering stochastic noise. Ultimately, this study demonstrates the importance of considering the effect of circadian perturbations on stochastic noise, particularly with regard to the development of small-molecule circadian therapeutics.     

Optimal transport analysis reveals trajectories in steady-state systems

Understanding how cells change their identity and behaviour in living systems is an important question in many fields of biology. The problem of inferring cell trajectories from single-cell measurements has been a major topic in the single-cell analysis community, with different methods developed for equilibrium and non-equilibrium systems (e.g. haematopoeisis vs. embryonic development). We show that optimal transport analysis, a technique originally designed for analysing time-courses, may also be applied to infer cellular trajectories from a single snapshot of a population in equilibrium. Therefore, optimal transport provides a unified approach to inferring trajectories that is applicable to both stationary and non-stationary systems. Our method, StationaryOT, is mathematically motivated in a natural way from the hypothesis of a Waddington’s epigenetic landscape. We implement StationaryOT as a software package and demonstrate its efficacy in applications to simulated data as well as single-cell data from Arabidopsis thaliana root development.     

Flexible copula model for integrating correlated multi-omics data from single-cell experiments

With recent advances in technologies to profile multi-omics data at the single-cell level, integrative multi-omics data analysis has been increasingly popular. It is increasingly common that information such as methylation changes, chromatin accessibility, and gene expression are jointly collected in a single-cell experiment. In biomedical studies, it is often of interest to study the associations between various data types and to examine how these associations might change according to other factors such as cell types and gene regulatory components. However, since each data type usually has a distinct marginal distribution, joint analysis of these changes of associations using multi-omics data is statistically challenging. In this paper, we propose a flexible copula-based framework to model covariate-dependent correlation structures independent of their marginals. In addition, the proposed approach could jointly combine a wide variety of univariate marginal distributions, either discrete or continuous, including the class of zero-inflated distributions. The performance of the proposed framework is demonstrated through a series of simulation studies. Finally, it is applied to a set of experimental data to investigate the dynamic relationship between single-cell RNA sequencing, chromatin accessibility, and DNA methylation at different germ layers during mouse gastrulation.     

SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing

The lion's share of bacteria in various environments cannot be cloned in the laboratory and thus cannot be sequenced using existing technologies. A major goal of single-cell genomics is to complement gene-centric metagenomic data with whole-genome assemblies of uncultivated organisms. Assembly of single-cell data is challenging because of highly non-uniform read coverage as well as elevated levels of sequencing errors and chimeric reads. We describe SPAdes, a new assembler for both single-cell and standard (multicell) assembly, and demonstrate that it improves on the recently released E+V-SC assembler (specialized for single-cell data) and on popular assemblers Velvet and SoapDeNovo (for multicell data). SPAdes generates single-cell assemblies, providing information about genomes of uncultivatable bacteria that vastly exceeds what may be obtained via traditional metagenomics studies. SPAdes is available online ( http://bioinf.spbau.ru/spades ). It is distributed as open source software.     

scLM: Automatic Detection of Consensus Gene Clusters Across Multiple Single-cell Datasets

In gene expression profiling studies, including single-cell RNA sequencing(sc RNA-seq)analyses, the identification and characterization of co-expressed genes provides critical information on cell identity and function. Gene co-expression clustering in sc RNA-seq data presents certain challenges. We show that commonly used methods for single-cell data are not capable of identifying co-expressed genes accurately, and produce results that substantially limit biological expectations of co-expressed genes. Herein, we present single-cell Latent-variable Model(sc LM), a gene coclustering algorithm tailored to single-cell data that performs well at detecting gene clusters with significant biologic context. Importantly, sc LM can simultaneously cluster multiple single-cell datasets, i.e., consensus clustering, enabling users to leverage single-cell data from multiple sources for novel comparative analysis. sc LM takes raw count data as input and preserves biological variation without being influenced by batch effects from multiple datasets. Results from both simulation data and experimental data demonstrate that sc LM outperforms the existing methods with considerably improved accuracy. To illustrate the biological insights of sc LM, we apply it to our in-house and public experimental sc RNA-seq datasets. sc LM identifies novel functional gene modules and refines cell states, which facilitates mechanism discovery and understanding of complex biosystems such as cancers. A user-friendly R package with all the key features of the sc LM method is available at https://github.com/QSong-github/sc LM.     

Single-cell copy number variation detection

Detection of chromosomal aberrations from a single cell by array comparative genomic hybridization (single-cell array CGH), instead of from a population of cells, is an emerging technique. However, such detection is challenging because of the genome artifacts and the DNA amplification process inherent to the single cell approach. Current normalization algorithms result in inaccurate aberration detection for single-cell data. We propose a normalization method based on channel, genome composition and recurrent genome artifact corrections. We demonstrate that the proposed channel clone normalization significantly improves the copy number variation detection in both simulated and real single-cell array CGH data.     

Recent progress in mass spectrometry for single-cell metabolomics

Metabolites are the end products of cellular vital activities and can reflect the state of cellular to a certain extent. Rapid change of metabolites and the low abundance of signature metabolites cause difficulties in single-cell detection, which is a great challenge in single-cell metabolomics analysis. Mass spectrometry (MS) is a powerful tool that uniquely suited to detect intracellular small-molecule metabolites and has shown good application in single-cell metabolite analysis. In this mini-review, we describe three types of emerging technologies for MS-based single-cell metabolic analysis in recent years, including nano-ESI-MS based single-cell metabolomics analysis, high-throughput analysis via flow cytometry, and cellular metabolic imaging analysis. These techniques provide a large amount of single-cell metabolic data, allowing the potential of MS in single-cell metabolic analysis is gradually being explored and is of great importance in disease and life science research.     

Applications of single-cell technology on bacterial analysis

Background: Traditionally, scientists studied microbiology through the manner of batch cultures, to conclude the dynamics or outputs by averaging all individuals. However, as the researches go further, the heterogeneities among the individuals have been proven to be crucial for the population dynamics and fates. Results: Due to the limit of technology, single-cell analysis methods were not widely used to decipher the inherent connections between individual cells and populations. Since the early decades of this century, the rapid development of microfluidics, fluorescent labelling, next-generation sequencing, and high-resolution microscopy have speeded up the development of single-cell technologies and further facilitated the applications of these technologies on bacterial analysis. Conclusions: In this review, we summarized the recent processes of single-cell technologies applied in bacterial analysis in terms of intracellular characteristics, cell physiology dynamics, and group behaviors, and discussed how single-cell technologies could be more applicable for future bacterial researches.     

Different levels of Notch signaling regulate quiescence, renewal and differentiation in pancreatic endocrine progenitors

Genetic studies have implicated Notch signaling in the maintenance of pancreatic progenitors. However, how Notch signaling regulates the quiescent, proliferative or differentiation behaviors of pancreatic progenitors at the single-cell level remains unclear. Here, using single-cell genetic analyses and a new transgenic system that allows dynamic assessment of Notch signaling, we address how discrete levels of Notch signaling regulate the behavior of endocrine progenitors in the zebrafish intrapancreatic duct. We find that these progenitors experience different levels of Notch signaling, which in turn regulate distinct cellular outcomes. High levels of Notch signaling induce quiescence, whereas lower levels promote progenitor amplification. The sustained downregulation of Notch signaling triggers a multistep process that includes cell cycle entry and progenitor amplification prior to endocrine differentiation. Importantly, progenitor amplification and differentiation can be uncoupled by modulating the duration and/or extent of Notch signaling downregulation, indicating that these processes are triggered by distinct levels of Notch signaling. These data show that different levels of Notch signaling drive distinct behaviors in a progenitor population.     

Lipid biosynthesis monitored at the single-cell level in Saccharomyces cerevisiae

There is increasing interest in bioengineering of lipids for use in functional foods, pharmaceuticals, and biofuels. Saccharomyces cerevisiae is a widely utilized cell factory for biotechnological production, thus a tempting alternative. Herein, we show how its neutral lipid accumulation varies throughout metabolic phases under nutritional conditions relevant for large-scale fermentation. Population-averaged metabolic data were correlated with lipid storage at the single-cell level monitored at submicron resolution by label-free coherent anti-Stokes Raman scattering (CARS) microscopy. While lipid droplet sizes are fairly constant, the number of droplets is a dynamic parameter determined by glucose and ethanol levels. The lowest number of lipid droplets is observed in the transition phase between glucose and ethanol fermentation. It is followed by a buildup during the ethanol phase. The surplus of accumulated lipids is then mobilized at concurrent glucose and ethanol starvation in the subsequent stationary phase. Thus, the highest amount of lipids is found in the ethanol phase, which is about 0.3 fL/cell. Our results indicate that the budding yeast, S. cerevisiae, can be used for the biosynthesis of lipids and demonstrate the strength of CARS microscopy for monitoring the dynamics of lipid metabolism at the single-cell level of importance for optimized lipid production.