Metabolic adaptation pilots the differentiation of human hematopoietic cells.

A continuous supply of energy is an essential prerequisite for survival and represents the highest priority for the cell. We hypothesize that cell differentiation is a process of optimization of energy flow in a changing environment through phenotypic adaptation. The mechanistic basis of this hypothesis is provided by the established link between core energy metabolism and epigenetic covalent modifications of chromatin. This theory predicts that early metabolic perturbations impact subsequent differentiation. To test this, we induced transient metabolic perturbations in undifferentiated human hematopoietic cells using pharmacological inhibitors targeting key metabolic reactions. We recorded changes in chromatin structure and gene expression, as well as phenotypic alterations by single-cell ATAC and RNA sequencing, time-lapse microscopy, and flow cytometry. Our observations suggest that these metabolic perturbations are shortly followed by alterations in chromatin structure, leading to changes in gene expression. We also show that these transient fluctuations alter the differentiation potential of the cells.

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.

Patient-derived organoids identify tailored therapeutic options and determinants of plasticity in sarcomatoid urothelial bladder cancer.

Sarcomatoid Urothelial Bladder Cancer (SARC) is a rare and aggressive histological subtype of bladder cancer for which therapeutic options are limited and experimental models are lacking. Here, we report the establishment of a long-term 3D organoid-like model derived from a SARC patient (SarBC-01). SarBC-01 emulates aggressive morphological, phenotypical, and transcriptional features of SARC and harbors somatic mutations in genes frequently altered in sarcomatoid tumors such as TP53 (p53) and RB1 (pRB). High-throughput drug screening, using a library comprising 1567 compounds in SarBC-01 and conventional urothelial carcinoma (UroCa) organoids, identified drug candidates active against SARC cells exclusively, or UroCa cells exclusively, or both. Among those, standard-of-care chemotherapeutic drugs inhibited both SARC and UroCa cells, while a subset of targeted drugs was specifically effective in SARC cells, including agents targeting the Glucocorticoid Receptor (GR) pathway. In two independent patient cohorts and in organoid models, GR and its encoding gene NR3C1 were found to be significantly more expressed in SARC as compared to UroCa, suggesting that high GR expression is a hallmark of SARC tumors. Further, glucocorticoid treatment impaired the mesenchymal morphology, abrogated the invasive ability of SARC cells, and led to transcriptomic changes associated with reversion of epithelial-to-mesenchymal transition, at single-cell level. Altogether, our study highlights the power of organoids for precision oncology and for providing key insights into factors driving rare tumor entities.

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.

Global genome decompaction leads to stochastic activation of gene expression as a first step toward fate commitment in human hematopoietic cells.

When human cord blood-derived CD34+ cells are induced to differentiate, they undergo rapid and dynamic morphological and molecular transformations that are critical for fate commitment. In particular, the cells pass through a transitory phase known as "multilineage-primed" state. These cells are characterized by a mixed gene expression profile, different in each cell, with the coexpression of many genes characteristic for concurrent cell lineages. The aim of our study is to understand the mechanisms of the establishment and the exit from this transitory state. We investigated this issue using single-cell RNA sequencing and ATAC-seq. Two phases were detected. The first phase is a rapid and global chromatin decompaction that makes most of the gene promoters in the genome accessible for transcription. It results 24 h later in enhanced and pervasive transcription of the genome leading to the concomitant increase in the cell-to-cell variability of transcriptional profiles. The second phase is the exit from the multilineage-primed phase marked by a slow chromatin closure and a subsequent overall down-regulation of gene transcription. This process is selective and results in the emergence of coherent expression profiles corresponding to distinct cell subpopulations. The typical time scale of these events spans 48 to 72 h. These observations suggest that the nonspecificity of genome decompaction is the condition for the generation of a highly variable multilineage expression profile. The nonspecific phase is followed by specific regulatory actions that stabilize and maintain the activity of key genes, while the rest of the genome becomes repressed again by the chromatin recompaction. Thus, the initiation of differentiation is reminiscent of a constrained optimization process that associates the spontaneous generation of gene expression diversity to subsequent regulatory actions that maintain the activity of some genes, while the rest of the genome sinks back to the repressive closed chromatin state.

Integrated time-lapse and single-cell transcription studies highlight the variable and dynamic nature of human hematopoietic cell fate commitment.

Individual cells take lineage commitment decisions in a way that is not necessarily uniform. We address this issue by characterising transcriptional changes in cord blood-derived CD34+ cells at the single-cell level and integrating data with cell division history and morphological changes determined by time-lapse microscopy. We show that major transcriptional changes leading to a multilineage-primed gene expression state occur very rapidly during the first cell cycle. One of the 2 stable lineage-primed patterns emerges gradually in each cell with variable timing. Some cells reach a stable morphology and molecular phenotype by the end of the first cell cycle and transmit it clonally. Others fluctuate between the 2 phenotypes over several cell cycles. Our analysis highlights the dynamic nature and variable timing of cell fate commitment in hematopoietic cells, links the gene expression pattern to cell morphology, and identifies a new category of cells with fluctuating phenotypic characteristics, demonstrating the complexity of the fate decision process (which is different from a simple binary switch between 2 options, as it is usually envisioned).