Differentiation is accompanied by a progressive loss in transcriptional memory.

BACKGROUND: Cell differentiation requires the integration of two opposite processes, a stabilizing cellular memory, especially at the transcriptional scale, and a burst of gene expression variability which follows the differentiation induction. Therefore, the actual capacity of a cell to undergo phenotypic change during a differentiation process relies upon a modification in this balance which favors change-inducing gene expression variability. However, there are no experimental data providing insight on how fast the transcriptomes of identical cells would diverge on the scale of the very first two cell divisions during the differentiation process. RESULTS: In order to quantitatively address this question, we developed different experimental methods to recover the transcriptomes of related cells, after one and two divisions, while preserving the information about their lineage at the scale of a single cell division. We analyzed the transcriptomes of related cells from two differentiation biological systems (human CD34+ cells and T2EC chicken primary erythrocytic progenitors) using two different single-cell transcriptomics technologies (scRT-qPCR and scRNA-seq). CONCLUSIONS: We identified that the gene transcription profiles of differentiating sister cells are more similar to each other than to those of non-related cells of the same type, sharing the same environment and undergoing similar biological processes. More importantly, we observed greater discrepancies between differentiating sister cells than between self-renewing sister cells. Furthermore, a progressive increase in this divergence from first generation to second generation was observed when comparing differentiating cousin cells to self renewing cousin cells. Our results are in favor of a gradual erasure of transcriptional memory during the differentiation process.

TopoDoE: a design of experiment strategy for selection and refinement in ensembles of executable gene regulatory networks.

BACKGROUND: Inference of Gene Regulatory Networks (GRNs) is a difficult and long-standing question in Systems Biology. Numerous approaches have been proposed with the latest methods exploring the richness of single-cell data. One of the current difficulties lies in the fact that many methods of GRN inference do not result in one proposed GRN but in a collection of plausible networks that need to be further refined. In this work, we present a Design of Experiment strategy to use as a second stage after the inference process. It is specifically fitted for identifying the next most informative experiment to perform for deciding between multiple network topologies, in the case where proposed GRNs are executable models. This strategy first performs a topological analysis to reduce the number of perturbations that need to be tested, then predicts the outcome of the retained perturbations by simulation of the GRNs and finally compares predictions with novel experimental data. RESULTS: We apply this method to the results of our divide-and-conquer algorithm called WASABI, adapt its gene expression model to produce perturbations and compare our predictions with experimental results. We show that our networks were able to produce in silico predictions on the outcome of a gene knock-out, which were qualitatively validated for 48 out of 49 genes. Finally, we eliminate as many as two thirds of the candidate networks for which we could identify an incorrect topology, thus greatly improving the accuracy of our predictions. CONCLUSION: These results both confirm the inference accuracy of WASABI and show how executable gene expression models can be leveraged to further refine the topology of inferred GRNs. We hope this strategy will help systems biologists further explore their data and encourage the development of more executable GRN models.

Evidence for close molecular proximity between reverting and undifferentiated cells.

BACKGROUND: According to Waddington's epigenetic landscape concept, the differentiation process can be illustrated by a cell akin to a ball rolling down from the top of a hill (proliferation state) and crossing furrows before stopping in basins or "attractor states" to reach its stable differentiated state. However, it is now clear that some committed cells can retain a certain degree of plasticity and reacquire phenotypical characteristics of a more pluripotent cell state. In line with this dynamic model, we have previously shown that differentiating cells (chicken erythrocytic progenitors (T2EC)) retain for 24 h the ability to self-renew when transferred back in self-renewal conditions. Despite those intriguing and promising results, the underlying molecular state of those "reverting" cells remains unexplored. The aim of the present study was therefore to molecularly characterize the T2EC reversion process by combining advanced statistical tools to make the most of single-cell transcriptomic data. For this purpose, T2EC, initially maintained in a self-renewal medium (0H), were induced to differentiate for 24H (24H differentiating cells); then, a part of these cells was transferred back to the self-renewal medium (48H reverting cells) and the other part was maintained in the differentiation medium for another 24H (48H differentiating cells). For each time point, cell transcriptomes were generated using scRT-qPCR and scRNAseq. RESULTS: Our results showed a strong overlap between 0H and 48H reverting cells when applying dimensional reduction. Moreover, the statistical comparison of cell distributions and differential expression analysis indicated no significant differences between these two cell groups. Interestingly, gene pattern distributions highlighted that, while 48H reverting cells have gene expression pattern more similar to 0H cells, they are not completely identical, which suggest that for some genes a longer delay may be required for the cells to fully recover. Finally, sparse PLS (sparse partial least square) analysis showed that only the expression of 3 genes discriminates 48H reverting and 0H cells. CONCLUSIONS: Altogether, we show that reverting cells return to an earlier molecular state almost identical to undifferentiated cells and demonstrate a previously undocumented physiological and molecular plasticity during the differentiation process, which most likely results from the dynamic behavior of the underlying molecular network.

The quiescent fraction of chronic myeloid leukemic stem cells depends on BMPR1B, Stat3 and BMP4-niche signals to persist in patients in remission.

Chronic myelogenous leukemia arises from the transformation of hematopoietic stem cells by the BCR-ABL oncogene. Though transformed cells are predominantly BCR-ABL-dependent and sensitive to tyrosine kinase inhibitor treatment, some BMPR1B+ leukemic stem cells are treatment-insensitive and rely, among others, on the bone morphogenetic protein (BMP) pathway for their survival via a BMP4 autocrine loop. Here, we further studied the involvement of BMP signaling in favoring residual leukemic stem cell persistence in the bone marrow of patients having achieved remission under treatment. We demonstrate by single-cell RNA-Seq analysis that a sub-fraction of surviving BMPR1B+ leukemic stem cells are co-enriched in BMP signaling, quiescence and stem cell signatures, without modulation of the canonical BMP target genes, but enrichment in actors of the Jak2/Stat3 signaling pathway. Indeed, based on a new model of persisting CD34+CD38- leukemic stem cells, we show that BMPR1B+ cells display co-activated Smad1/5/8 and Stat3 pathways. Interestingly, we reveal that only the BMPR1B+ cells adhering to stromal cells display a quiescent status. Surprisingly, this quiescence is induced by treatment, while non-adherent BMPR1B+ cells treated with tyrosine kinase inhibitors continued to proliferate. The subsequent targeting of BMPR1B and Jak2 pathways decreased quiescent leukemic stem cells by promoting their cell cycle re-entry and differentiation. Moreover, while Jak2-inhibitors alone increased BMP4 production by mesenchymal cells, the addition of the newly described BMPR1B inhibitor (E6201) impaired BMP4-mediated production by stromal cells. Altogether, our data demonstrate that targeting both BMPR1B and Jak2/Stat3 efficiently impacts persisting and dormant leukemic stem cells hidden in their bone marrow microenvironment.

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.

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.

An image-guided microfluidic system for single-cell lineage tracking.

Cell lineage tracking is a long-standing and unresolved problem in biology. Microfluidic technologies have the potential to address this problem, by virtue of their ability to manipulate and process single-cells in a rapid, controllable and efficient manner. Indeed, when coupled with traditional imaging approaches, microfluidic systems allow the experimentalist to follow single-cell divisions over time. Herein, we present a valve-based microfluidic system able to probe the decision-making processes of single-cells, by tracking their lineage over multiple generations. The system operates by trapping single-cells within growth chambers, allowing the trapped cells to grow and divide, isolating sister cells after a user-defined number of divisions and finally extracting them for downstream transcriptome analysis. The platform incorporates multiple cell manipulation operations, image processing-based automation for cell loading and growth monitoring, reagent addition and device washing. To demonstrate the efficacy of the microfluidic workflow, 6C2 (chicken erythroleukemia) and T2EC (primary chicken erythrocytic progenitors) cells are tracked inside the microfluidic device over two generations, with a cell viability rate in excess of 90%. Sister cells are successfully isolated after division and extracted within a 500 nL volume, which was demonstrated to be compatible with downstream single-cell RNA sequencing analysis.

Drugs modulating stochastic gene expression affect the erythroid differentiation process.

To better understand the mechanisms behind cells decision-making to differentiate, we assessed the influence of stochastic gene expression (SGE) modulation on the erythroid differentiation process. It has been suggested that stochastic gene expression has a role in cell fate decision-making which is revealed by single-cell analyses but studies dedicated to demonstrate the consistency of this link are still lacking. Recent observations showed that SGE significantly increased during differentiation and a few showed that an increase of the level of SGE is accompanied by an increase in the differentiation process. However, a consistent relation in both increasing and decreasing directions has never been shown in the same cellular system. Such demonstration would require to be able to experimentally manipulate simultaneously the level of SGE and cell differentiation in order to observe if cell behavior matches with the current theory. We identified three drugs that modulate SGE in primary erythroid progenitor cells. Both Artemisinin and Indomethacin decreased SGE and reduced the amount of differentiated cells. On the contrary, a third component called MB-3 simultaneously increased the level of SGE and the amount of differentiated cells. We then used a dynamical modelling approach which confirmed that differentiation rates were indeed affected by the drug treatment. Using single-cell analysis and modeling tools, we provide experimental evidence that, in a physiologically relevant cellular system, SGE is linked to differentiation.

Erythroid differentiation displays a peak of energy consumption concomitant with glycolytic metabolism rearrangements.

Our previous single-cell based gene expression analysis pointed out significant variations of LDHA level during erythroid differentiation. Deeper investigations highlighted that a metabolic switch occurred along differentiation of erythroid cells. More precisely we showed that self-renewing progenitors relied mostly upon lactate-productive glycolysis, and required LDHA activity, whereas differentiating cells, mainly involved mitochondrial oxidative phosphorylation (OXPHOS). These metabolic rearrangements were coming along with a particular temporary event, occurring within the first 24h of erythroid differentiation. The activity of glycolytic metabolism and OXPHOS rose jointly with oxgene consumption dedicated to ATP production at 12-24h of the differentiation process before lactate-productive glycolysis sharply fall down and energy needs decline. Finally, we demonstrated that the metabolic switch mediated through LDHA drop and OXPHOS upkeep might be necessary for erythroid differentiation. We also discuss the possibility that metabolism, gene expression and epigenetics could act together in a circular manner as a driving force for differentiation.

Automated cell cycle and cell size measurements for single-cell gene expression studies.

OBJECTIVES: Recent rise of single-cell studies revealed the importance of understanding the role of cell-to-cell variability, especially at the transcriptomic level. One of the numerous sources of cell-to-cell variation in gene expression is the heterogeneity in cell proliferation state. In order to identify how cell cycle and cell size influences gene expression variability at the single-cell level, we provide an universal and automatic toxic-free label method, compatible with single-cell high-throughput RT-qPCR. The method consists of isolating cells after a double-stained, analyzing their morphological parameters and performing a transcriptomic analysis on the same identified cells. RESULTS: This led to an unbiased gene expression analysis and could be also used for improving single-cell tracking and imaging when combined with cell isolation. As an application for this technique, we showed that cell-to-cell variability in chicken erythroid progenitors was negligibly influenced by cell size nor cell cycle.

Large-scale analysis by SAGE reveals new mechanisms of v-erbA oncogene action.

BACKGROUND: The v-erbA oncogene, carried by the Avian Erythroblastosis Virus, derives from the c-erbAalpha proto-oncogene that encodes the nuclear receptor for triiodothyronine (T3R). v-ErbA transforms erythroid progenitors in vitro by blocking their differentiation, supposedly by interference with T3R and RAR (Retinoic Acid Receptor). However, v-ErbA target genes involved in its transforming activity still remain to be identified. RESULTS: By using Serial Analysis of Gene Expression (SAGE), we identified 110 genes deregulated by v-ErbA and potentially implicated in the transformation process. Bioinformatic analysis of promoter sequence and transcriptional assays point out a potential role of c-Myb in the v-ErbA effect. Furthermore, grouping of newly identified target genes by function revealed both expected (chromatin/transcription) and unexpected (protein metabolism) functions potentially deregulated by v-ErbA. We then focused our study on 15 of the new v-ErbA target genes and demonstrated by real time PCR that in majority their expression was activated neither by T3, nor RA, nor during differentiation. This was unexpected based upon the previously known role of v-ErbA. CONCLUSION: This paper suggests the involvement of a wealth of new unanticipated mechanisms of v-ErbA action.

Global transcription analysis of immature avian erythrocytic progenitors: from self-renewal to differentiation.

The molecular mechanisms regulating the cell fate decision between self-renewal and differentiation/apoptosis in stem and progenitor cells are poorly understood. Here, we report the first comprehensive identification of genes potentially involved in the switch from self-renewal toward differentiation of primary, non-immortalized erythroid avian progenitor cells (T2EC cells). We used the Serial Analysis of Gene Expression (SAGE) technique in order to identify and quantify the genome fraction functionally active in a self-renewing versus a differentiating cell population. We generated two SAGE libraries and sequenced a total of 37,589 tags, thereby obtaining the first transcriptional profile characterization of a chicken cell. Tag identification was performed using a new relational database (Identitag) developed in the laboratory, which allowed a highly satisfactory level of identification. Among 123 differentially expressed genes, 11 were investigated further and for nine of them the differential expression was subsequently confirmed by real-time PCR. The comparison of tag abundance between the two libraries revealed that only a small fraction of transcripts was differentially expressed. The analysis of their functions argue against a prominent role for a master switch in T2EC cells decision-making, but are in favor of a critical role for coordinated small variations in a relatively small number of genes that can lead to essential cellular identity changes.

Identification of human, rat and chicken ribosomal proteins by a combination of two-dimensional polyacrylamide gel electrophoresis and mass spectrometry.

To identify the exact spot position of human, rat and chicken ribosomal proteins (RP) separated by two-dimensional polyacrylamide gel electrophoresis (2-DE), a 2-DE system was designed to separate RP with a pI>8.6 according to their charge in the first dimension and to their molecular mass in the second dimension. Individual proteins were excised from the gels and identified by mass spectrometry after digestion by trypsin. In addition, a mixture of purified RP from these three species was also analyzed by tandem mass tag spectrometry. By combining those two methods 74 RP from human, 76 from rat and 67 from chicken were identified according to the nomenclature initially defined for rat liver RP and by using the Swiss-Prot/trEMBL databases. Whereas human and rat RP were well described, most of RP from chicken were not characterized in databases, since 35 out of 67 chicken RP identified in this study were not listed yet. We propose here the first comprehensive description of chicken RP and their comparison to those from human and rat.