'Particle genetics': treating every cell as unique

Cell-to-cell expression dispersion of B-cell surface proteins is linked to genetic variants in humans

Variability in gene expression across a population of homogeneous cells is known to influence various biological processes. In model organisms, natural genetic variants were found that modify expression dispersion (variability at a fixed mean) but very few studies have detected such effects in humans. Here, we analyzed single-cell expression of four proteins (CD23, CD55, CD63 and CD86) across cell lines derived from individuals of the Yoruba population. Using data from over 30 million cells, we found substantial inter-individual variation of dispersion. We demonstrate, via de novo cell line generation and subcloning experiments, that this variation exceeds the variation associated with cellular immortalization. We detected a genetic association between the expression dispersion of CD63 and the rs971 SNP. Our results show that human DNA variants can have inherently-probabilistic effects on gene expression. Such subtle genetic effects may participate to phenotypic variation and disease outcome.

Triqueneaux, Burny, Symmons et al. show association between gene expression noise and genotypes, using single-cell expression of four proteins across human-derived lymphoblastoid cell lines. This study suggests that very subtle regulatory effects of human DNA variants may contribute to phenotypic variation and disease outcome.

Poikilosis - pervasive biological variation

Biological systems are dynamic and display heterogeneity at all levels. Ubiquitous heterogeneity, here called for poikilosis, is an integral and important property of organisms and in molecules, systems and processes within them. Traditionally, heterogeneity in biology and experiments has been considered as unwanted noise, here poikilosis is shown to be the normal state. Acceptable variation ranges are called as lagom. Non-lagom, variations that are too extensive, have negative effects, which influence interconnected levels and once the variation is large enough cause a disease and can lead even to death. Poikilosis has numerous applications and consequences e.g. for how to design, analyze and report experiments, how to develop and apply prediction and modelling methods, and in diagnosis and treatment of diseases. Poikilosis-aware new and practical definitions are provided for life, death, senescence, disease, and lagom. Poikilosis is the first new unifying theory in biology since evolution and should be considered in every scientific study.

Poikilosis - pervasive biological variation

Biological systems are dynamic and display heterogeneity at all levels. Ubiquitous heterogeneity, here called for poikilosis, is an integral and important property of organisms and in molecules, systems and processes within them. Traditionally, heterogeneity in biology and experiments has been considered as unwanted noise, here poikilosis is shown to be the normal state. Acceptable variation ranges are called as lagom. Non-lagom, variations that are too extensive, have negative effects, which influence interconnected levels and once the variation is large enough cause a disease and can lead even to death. Poikilosis has numerous applications and consequences e.g. for how to design, analyze and report experiments, how to develop and apply prediction and modelling methods, and in diagnosis and treatment of diseases. Poikilosis-aware new and practical definitions are provided for life, death, senescence, disease, and lagom. Poikilosis is the first new unifying theory in biology since evolution and should be considered in every scientific study.

Transcriptional adaptation in Caenorhabditis elegans

Transcriptional adaptation is a recently described phenomenon by which a mutation in one gene leads to the transcriptional modulation of related genes, termed adapting genes. At the molecular level, it has been proposed that the mutant mRNA, rather than the loss of protein function, activates this response. While several examples of transcriptional adaptation have been reported in zebrafish embryos and in mouse cell lines, it is not known whether this phenomenon is observed across metazoans. Here we report transcriptional adaptation in C. elegans, and find that this process requires factors involved in mutant mRNA decay, as in zebrafish and mouse. We further uncover a requirement for Argonaute proteins and Dicer, factors involved in small RNA maturation and transport into the nucleus. Altogether, these results provide evidence for transcriptional adaptation in C. elegans, a powerful model to further investigate underlying molecular mechanisms.

WASABI: a dynamic iterative framework for gene regulatory network inference

BACKGROUND

Inference of gene regulatory networks from gene expression data has been a long-standing and notoriously difficult task in systems biology. Recently, single-cell transcriptomic data have been massively used for gene regulatory network inference, with both successes and limitations.

RESULTS

In the present work we propose an iterative algorithm called WASABI, dedicated to inferring a causal dynamical network from time-stamped single-cell data, which tackles some of the limitations associated with current approaches. We first introduce the concept of waves, which posits that the information provided by an external stimulus will affect genes one-by-one through a cascade, like waves spreading through a network. This concept allows us to infer the network one gene at a time, after genes have been ordered regarding their time of regulation. We then demonstrate the ability of WASABI to correctly infer small networks, which have been simulated in silico using a mechanistic model consisting of coupled piecewise-deterministic Markov processes for the proper description of gene expression at the single-cell level. We finally apply WASABI on in vitro generated data on an avian model of erythroid differentiation. The structure of the resulting gene regulatory network sheds a new light on the molecular mechanisms controlling this process. In particular, we find no evidence for hub genes and a much more distributed network structure than expected. Interestingly, we find that a majority of genes are under the direct control of the differentiation-inducing stimulus.

CONCLUSIONS

Together, these results demonstrate WASABI versatility and ability to tackle some general gene regulatory networks inference issues. It is our hope that WASABI will prove useful in helping biologists to fully exploit the power of time-stamped single-cell data.

Resolving the Complex Genetic Basis of Phenotypic Variation and Variability of Cellular Growth

In all organisms, the majority of traits vary continuously between individuals. Explaining the genetic basis of quantitative trait variation requires comprehensively accounting for genetic and nongenetic factors as well as their interactions. The growth of microbial cells can be characterized by a lag duration, an exponential growth phase, and a stationary phase. Parameters that characterize these growth phases can vary among genotypes (phenotypic variation), environmental conditions (phenotypic plasticity), and among isogenic cells in a given environment (phenotypic variability). We used a high-throughput microscopy assay to map genetic loci determining variation in lag duration and exponential growth rate in growth rate-limiting and nonlimiting glucose concentrations, using segregants from a cross of two natural isolates of the budding yeast, Saccharomyces cerevisiae We find that some quantitative trait loci (QTL) are common between traits and environments whereas some are unique, exhibiting gene-by-environment interactions. Furthermore, whereas variation in the central tendency of growth rate or lag duration is explained by many additive loci, differences in phenotypic variability are primarily the result of genetic interactions. We used bulk segregant mapping to increase QTL resolution by performing whole-genome sequencing of complex mixtures of an advanced intercross mapping population grown in selective conditions using glucose-limited chemostats. We find that sequence variation in the high-affinity glucose transporter HXT7 contributes to variation in growth rate and lag duration. Allele replacements of the entire locus, as well as of a single polymorphic amino acid, reveal that the effect of variation in HXT7 depends on genetic, and allelic, background. Amplifications of HXT7 are frequently selected in experimental evolution in glucose-limited environments, but we find that HXT7 amplifications result in antagonistic pleiotropy that is absent in naturally occurring variants of HXT7 Our study highlights the complex nature of the genotype-to-phenotype map within and between environments.

Single-Cell-Based Analysis Highlights a Surge in Cell-to-Cell Molecular Variability Preceding Irreversible Commitment in a Differentiation Process

In some recent studies, a view emerged that stochastic dynamics governing the switching of cells from one differentiation state to another could be characterized by a peak in gene expression variability at the point of fate commitment. We have tested this hypothesis at the single-cell level by analyzing primary chicken erythroid progenitors through their differentiation process and measuring the expression of selected genes at six sequential time-points after induction of differentiation. In contrast to population-based expression data, single-cell gene expression data revealed a high cell-to-cell variability, which was masked by averaging. We were able to show that the correlation network was a very dynamical entity and that a subgroup of genes tend to follow the predictions from the dynamical network biomarker (DNB) theory. In addition, we also identified a small group of functionally related genes encoding proteins involved in sterol synthesis that could act as the initial drivers of the differentiation. In order to assess quantitatively the cell-to-cell variability in gene expression and its evolution in time, we used Shannon entropy as a measure of the heterogeneity. Entropy values showed a significant increase in the first 8 h of the differentiation process, reaching a peak between 8 and 24 h, before decreasing to significantly lower values. Moreover, we observed that the previous point of maximum entropy precedes two paramount key points: an irreversible commitment to differentiation between 24 and 48 h followed by a significant increase in cell size variability at 48 h. In conclusion, when analyzed at the single cell level, the differentiation process looks very different from its classical population average view. New observables (like entropy) can be computed, the behavior of which is fully compatible with the idea that differentiation is not a “simple” program that all cells execute identically but results from the dynamical behavior of the underlying molecular network.

A single-cell transcriptomics analysis offers a new dynamical view of the differentiation process, involving an increase in between-cell variability prior to commitment.

The differentiation process has classically been seen as a stereotyped program leading from one progenitor toward a functional cell. This vision was based upon cell population-based analyses averaged over millions of cells. However, new methods have recently emerged that allow interrogation of the molecular content at the single-cell level, challenging this view with a new model suggesting that cell-to-cell gene expression stochasticity could play a key role in differentiation. We took advantage of a physiologically relevant avian cellular model to analyze the expression level of 92 genes in individual cells collected at several time-points during differentiation. We first observed that the process analyzed at the single-cell level is very different and much less well ordered than the population-based average view. Furthermore, we showed that cell-to-cell variability in gene expression peaks transiently before strongly decreasing. This rise in variability precedes two key events: an irreversible commitment to differentiation, followed by a significant increase in cell size variability. Altogether, our results support the idea that differentiation is not a “simple” series of well-ordered molecular events executed identically by all cells in a population but likely results from dynamical behavior of the underlying molecular network.

Exploiting Single-Cell Quantitative Data to Map Genetic Variants Having Probabilistic Effects

Despite the recent progress in sequencing technologies, genome-wide association studies (GWAS) remain limited by a statistical-power issue: many polymorphisms contribute little to common trait variation and therefore escape detection. The small contribution sometimes corresponds to incomplete penetrance, which may result from probabilistic effects on molecular regulations. In such cases, genetic mapping may benefit from the wealth of data produced by single-cell technologies. We present here the development of a novel genetic mapping method that allows to scan genomes for single-cell Probabilistic Trait Loci that modify the statistical properties of cellular-level quantitative traits. Phenotypic values are acquired on thousands of individual cells, and genetic association is obtained from a multivariate analysis of a matrix of Kantorovich distances. No prior assumption is required on the mode of action of the genetic loci involved and, by exploiting all single-cell values, the method can reveal non-deterministic effects. Using both simulations and yeast experimental datasets, we show that it can detect linkages that are missed by classical genetic mapping. A probabilistic effect of a single SNP on cell shape was detected and validated. The method also detected a novel locus associated with elevated gene expression noise of the yeast galactose regulon. Our results illustrate how single-cell technologies can be exploited to improve the genetic dissection of certain common traits. The method is available as an open source R package called ptlmapper.

Stochastic developmental variation, an epigenetic source of phenotypic diversity with far-reaching biological consequences

This article reviews the production of different phenotypes from the same genotype in the same environment by stochastic cellular events, nonlinear mechanisms during patterning and morphogenesis, and probabilistic self-reinforcing circuitries in the adult life. These aspects of phenotypic variation are summarized under the term 'stochastic developmental variation' (SDV) in the following. In the past, SDV has been viewed primarily as a nuisance, impairing laboratory experiments, pharmaceutical testing, and true-to-type breeding. This article also emphasizes the positive biological effects of SDV and discusses implications for genotype-to-phenotype mapping, biological individuation, ecology, evolution, and applied biology. There is strong evidence from experiments with genetically identical organisms performed in narrowly standardized laboratory set-ups that SDV is a source of phenotypic variation in its own right aside from genetic variation and environmental variation. It is obviously mediated by molecular and higher-order epigenetic mechanisms. Comparison of SDV in animals, plants, fungi, protists, bacteria, archaeans, and viruses suggests that it is a ubiquitous and phylogenetically old phenomenon. In animals, it is usually smallest for morphometric traits and highest for life history traits and behaviour. SDV is thought to contribute to phenotypic diversity in all populations but is particularly relevant for asexually reproducing and genetically impoverished populations, where it generates individuality despite genetic uniformity. In each generation, SDV produces a range of phenotypes around a well-adapted target phenotype, which is interpreted as a bet-hedging strategy to cope with the unpredictability of dynamic environments. At least some manifestations of SDV are heritable, adaptable, selectable, and evolvable, and therefore, SDV may be seen as a hitherto overlooked evolution factor. SDV is also relevant for husbandry, agriculture, and medicine because most pathogens are asexuals that exploit this third source of phenotypic variation to modify infectivity and resistance to antibiotics. Since SDV affects all types of organisms and almost all aspects of life, it urgently requires more intense research and a better integration into biological thinking.

Variation is function: Are single cell differences functionally important?: Testing the hypothesis that single cell variation is required for aggregate function

There is a growing appreciation of the extent of transcriptome variation across individual cells of the same cell type. While expression variation may be a byproduct of, for example, dynamic or homeostatic processes, here we consider whether single-cell molecular variation per se might be crucial for population-level function. Under this hypothesis, molecular variation indicates a diversity of hidden functional capacities within an ensemble of identical cells, and this functional diversity facilitates collective behavior that would be inaccessible to a homogenous population. In reviewing this topic, we explore possible functions that might be carried by a heterogeneous ensemble of cells; however, this question has proven difficult to test, both because methods to manipulate molecular variation are limited and because it is complicated to define, and measure, population-level function. We consider several possible methods to further pursue the hypothesis that variation is function through the use of comparative analysis and novel experimental techniques.

Behavioral idiosyncrasy reveals genetic control of phenotypic variability

Quantitative genetics has primarily focused on describing genetic effects on trait means and largely ignored the effect of alternative alleles on trait variability, potentially missing an important axis of genetic variation contributing to phenotypic differences among individuals. To study the genetic effects on individual-to-individual phenotypic variability (or intragenotypic variability), we used Drosophila inbred lines and measured the spontaneous locomotor behavior of flies walking individually in Y-shaped mazes, focusing on variability in locomotor handedness, an assay optimized to measure variability. We discovered that some lines had consistently high levels of intragenotypic variability among individuals, whereas lines with low variability behaved as although they tossed a coin at each left/right turn decision. We demonstrate that the degree of variability is itself heritable. Using a genome-wide association study (GWAS) for the degree of intragenotypic variability as the phenotype across lines, we identified several genes expressed in the brain that affect variability in handedness without affecting the mean. One of these genes, Ten-a, implicates a neuropil in the central complex of the fly brain as influencing the magnitude of behavioral variability, a brain region involved in sensory integration and locomotor coordination. We validated these results using genetic deficiencies, null alleles, and inducible RNAi transgenes. Our study reveals the constellation of phenotypes that can arise from a single genotype and shows that different genetic backgrounds differ dramatically in their propensity for phenotypic variabililty. Because traditional mean-focused GWASs ignore the contribution of variability to overall phenotypic variation, current methods may miss important links between genotype and phenotype.

How does evolution tune biological noise?

Part of molecular and phenotypic differences between individual cells, between body parts, or between individuals can result from biological noise. This source of variation is becoming more and more apparent thanks to the recent advances in dynamic imaging and single-cell analysis. Some of these studies showed that the link between genotype and phenotype is not strictly deterministic. Mutations can change various statistical properties of a biochemical reaction, and thereby the probability of a trait outcome. The fact that they can modulate phenotypic noise brings up an intriguing question: how may selection act on these mutations? In this review, we approach this question by first covering the evidence that biological noise is under genetic control and therefore a substrate for evolution. We then sequentially inspect the possibilities of negative, neutral, and positive selection for mutations increasing biological noise. Finally, we hypothesize on the specific case of H2A.Z, which was shown to both buffer phenotypic noise and modulate transcriptional efficiency.

Changing genetic paradigms: creating next-generation genetic databases as tools to understand the emerging complexities of genotype/phenotype relationships

Understanding genotype/phenotype relationships has become more complicated as increasing amounts of inter- and intra-tissue genetic heterogeneity have been revealed through next-generation sequencing and evidence showing that factors such as epigenetic modifications, non-coding RNAs and RNA editing can play an important role in determining phenotype. Such findings have challenged a number of classic genetic assumptions including (i) analysis of genomic sequence obtained from blood is an accurate reflection of the genotype responsible for phenotype expression in an individual; (ii) that significant genetic alterations will be found only in diseased individuals, in germline tissues in inherited diseases, or in specific diseased tissues in somatic diseases such as cancer; and (iii) that mutation rates in putative disease-associated genes solely determine disease phenotypes. With the breakdown of our traditional understanding of genotype to phenotype relationships, it is becoming increasingly apparent that new analytical tools will be required to determine the relationship between genotype and phenotypic expression. To this end, we are proposing that next-generation genetic database (NGDB) platforms be created that include new bioinformatics tools based on algorithms that can evaluate genetic heterogeneity, as well as powerful systems biology analysis tools to actively process and evaluate the vast amounts of both genomic and genomic-modifying information required to reveal the true relationships between genotype and phenotype.