What's Luck Got to Do with It: Single Cells, Multiple Fates, and Biological Nondeterminism

Transcriptional silencing in Saccharomyces cerevisiae: known unknowns

Transcriptional silencing in Saccharomyces cerevisiae is a persistent and highly stable form of gene repression. It involves DNA silencers and repressor proteins that bind nucleosomes. The silenced state is influenced by numerous factors including the concentration of repressors, nature of activators, architecture of regulatory elements, modifying enzymes and the dynamics of chromatin.Silencers function to increase the residence time of repressor Sir proteins at silenced domains while clustering of silenced domains enables increased concentrations of repressors and helps facilitate long-range interactions. The presence of an accessible NDR at the regulatory regions of silenced genes, the cycling of chromatin configurations at regulatory sites, the mobility of Sir proteins, and the non-uniform distribution of the Sir proteins across the silenced domain, all result in silenced chromatin that only stably silences weak promoters and enhancers via changes in transcription burst duration and frequency.These data collectively suggest that silencing is probabilistic and the robustness of silencing is achieved through sub-optimization of many different nodes of action such that a stable expression state is generated and maintained even though individual constituents are in constant flux.

Global transcription regulation revealed from dynamical correlations in time-resolved single-cell RNA sequencing

Single-cell transcriptomics reveals significant variations in transcriptional activity across cells. Yet, it remains challenging to identify mechanisms of transcription dynamics from static snapshots. It is thus still unknown what drives global transcription dynamics in single cells. We present a stochastic model of gene expression with cell size- and cell cycle-dependent rates in growing and dividing cells that harnesses temporal dimensions of single-cell RNA sequencing through metabolic labeling protocols and cel lcycle reporters. We develop a parallel and highly scalable approximate Bayesian computation method that corrects for technical variation and accurately quantifies absolute burst frequency, burst size, and degradation rate along the cell cycle at a transcriptome-wide scale. Using Bayesian model selection, we reveal scaling between transcription rates and cell size and unveil waves of gene regulation across the cell cycle-dependent transcriptome. Our study shows that stochastic modeling of dynamical correlations identifies global mechanisms of transcription regulation. A record of this paper's transparent peer review process is included in the supplemental information.

Exploring the Surface: Sampling of Potential Skin Cancer Biomarkers Kynurenine and Tryptophan, Studied on 3D Melanocyte and Melanoma Models

Early detection of cancer via biomarkers is vital for improving patient survival rates. In the case of skin cancers, low-molecular-weight biomarkers can penetrate the skin barrier, enabling non-invasive sampling at an early stage. This study focuses on detecting tryptophan (Trp) and kynurenine (Kyn) on the surface of reconstructed 3D melanoma and melanocyte models. This is examined in connection with IDO-1 and IL-6 expression in response to IFN-γ or UVB stimulation, both crucial factors of the melanoma tumor microenvironment (TME). Using a polystyrene scaffold, full-thickness human skin equivalents containing fibroblasts, keratinocytes, and melanocytes or melanoma cells were developed. The samples were stimulated with IFN-γ or UVB, and Trp and Kyn secretion was measured using HPLC-PDA and HPLC-MS. The expression of IDO-1 and IL-6 was measured using RT-qPCR. Increased Trp catabolism to Kyn was observed in IFN-γ-stimulated melanoma and melanocyte models, along with higher IDO-1 expression. UVB exposure led to significant changes in Kyn levels but only in the melanoma model. This study demonstrates the potential of skin surface Trp and Kyn monitoring to capture TME metabolic changes. It also lays the groundwork for future in vivo studies, aiding in understanding and monitoring skin cancer progression.

Single-nucleus proteomics identifies regulators of protein transport

The physiological response of a cell to stimulation depends on its proteome configuration. Therefore, the abundance variation of regulatory proteins across unstimulated single cells can be associatively linked with their response to stimulation. Here we developed an approach that leverages this association across individual cells and nuclei to systematically identify potential regulators of biological processes, followed by targeted validation. Specifically, we applied this approach to identify regulators of nucleocytoplasmic protein transport in macrophages stimulated with lipopolysaccharide (LPS). To this end, we quantified the proteomes of 3,412 individual nuclei, sampling the dynamic response to LPS treatment, and linking functional variability to proteomic variability. Minutes after the stimulation, the protein transport in individual nuclei correlated strongly with the abundance of known protein transport regulators, thus revealing the impact of natural protein variability on functional cellular response. We found that simple biophysical constraints, such as the quantity of nuclear pores, partially explain the variability in LPS-induced nucleocytoplasmic transport. Among the many proteins newly identified to be associated with the response, we selected 16 for targeted validation by knockdown. The knockdown phenotypes confirmed the inferences derived from natural protein and functional variation of single nuclei, thus demonstrating the potential of (sub-)single-cell proteomics to infer functional regulation. We expect this approach to generalize to broad applications and enhance the functional interpretability of single-cell omics data.

Gram-positive bacteria are primed for surviving lethal doses of antibiotics and chemical stress

Antibiotic resistance kills millions worldwide yearly. However, a major contributor to recurrent infections lies in a small fraction of bacterial cells, known as persisters. These cells are not inherently antibiotic-resistant, yet they lead to increased antibiotic usage, raising the risk of developing resistant progenies. In a bacterial population, individual cells exhibit considerable fluctuations in their gene expression levels despite being cultivated under identical, stable conditions. This variability in cell-to-cell characteristics (phenotypic diversity) within an isogenic population enables persister cells to withstand antibiotic exposure by entering a non-dividing state. We recently showed the existence of "primed cells" in E. coli. Primed cells are dividing cells prepared for antibiotic stress before encountering it and are more prone to form persisters. They also pass their "prepared state" down for several generations through epigenetic memory. Here, we show that primed cells are common among distant bacterial lineages, allowing for survival against antibiotics and other chemical stress, and form in different growth phases. They are also responsible for increased persister levels in transition and stationary phases compared to the log phase. We tested and showed that the Gram-positive bacterium Bacillus megaterium, evolutionarily very distant from E. coli, forms primed cells and has a transient epigenetic memory that is maintained for 7 generations or more. We showed this using ciprofloxacin and the non-antibiotic chemical stress fluoride. It is well established that persister levels are higher in the stationary phase than in the log phase, and B. megaterium persisters levels are nearly identical from the early to late-log phase but are ~2-fold and ~4-fold higher in the transition and stationary phase, respectively. It was previously proposed that there are two distinct types of persisters: Type II forms in the log phase, while Type I forms in the stationary phase. However, we show that primed cells lead to increased persisters in the transition and stationary phase and found no evidence of Type I or II persisters with distant phenotypes. Overall, we have provided substantial evidence of the importance of primed cells and their transitory epigenetic memories to surviving stress.

Complementing Cell Taxonomies with a Multicellular Analysis of Tissues

The application of single-cell molecular profiling coupled with spatial technologies has enabled charting of cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single-cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single-cell and spatial datasets to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.

Unraveling IFN-I response dynamics and TNF crosstalk in the pathophysiology of systemic lupus erythematosus

Introduction

The innate immune system serves the crucial first line of defense against a wide variety of potential threats, during which the production of pro-inflammatory cytokines IFN-I and TNFα are key. This astonishing power to fight invaders, however, comes at the cost of risking IFN-I-related pathologies, such as observed during autoimmune diseases, during which IFN-I and TNFα response dynamics are dysregulated. Therefore, these response dynamics must be tightly regulated, and precisely matched with the potential threat. This regulation is currently far from understood.

Methods

Using droplet-based microfluidics and ODE modeling, we studied the fundamentals of single-cell decision-making upon TLR signaling in human primary immune cells (n = 23). Next, using biologicals used for treating autoimmune diseases [i.e., anti-TNFα, and JAK inhibitors], we unraveled the crosstalk between IFN-I and TNFα signaling dynamics. Finally, we studied primary immune cells isolated from SLE patients (n = 8) to provide insights into SLE pathophysiology.

Results

single-cell IFN-I and TNFα response dynamics display remarkable differences, yet both being highly heterogeneous. Blocking TNFα signaling increases the percentage of IFN-I-producing cells, while blocking IFN-I signaling decreases the percentage of TNFα-producing cells. Single-cell decision-making in SLE patients is dysregulated, pointing towards a dysregulated crosstalk between IFN-I and TNFα response dynamics.

Discussion

We provide a solid droplet-based microfluidic platform to study inherent immune secretory behaviors, substantiated by ODE modeling, which can challenge the conceptualization within and between different immune signaling systems. These insights will build towards an improved fundamental understanding on single-cell decision-making in health and disease.

Protocol to study the inheritance and propagation of non-genetically encoded states using barcode decay lineage tracing

Here, we present a protocol to perform barcode decay lineage tracing followed by single-cell transcriptome analysis (BdLT-Seq). We describe steps for BdLT-Seq experimental design, building barcoded episome reporters, performing episome transfection, and barcode retrieval. We then describe procedures for sequencing library construction while providing options for sample multiplexing and data analysis. This BdLT-Seq technique enables the assessment of clonal evolution in a directional manner while preserving isogeneity, thus allowing the comparison of non-genetic molecular features between isogenic cell lineages. For complete details on the use and execution of this protocol, please refer to Shlyakhtina et al. (2023).1.

Variability in Neural Circuit Formation

The study of neural development is usually concerned with the question of how nervous systems get put together. Variation in these processes is usually of interest as a means of revealing these normative mechanisms. However, variation itself can be an object of study and is of interest from multiple angles. First, the nature of variation in both the processes and the outcomes of neural development is relevant to our understanding of how these processes and outcomes are encoded in the genome. Second, variation in the wiring of the brain in humans may underlie variation in all kinds of psychological and behavioral traits, as well as neurodevelopmental disorders. And third, genetic variation that affects circuit development provides the raw material for evolutionary change. Here, I examine these different aspects of variation in circuit development and consider what they may tell us about these larger questions.

Retrospective identification of cell-intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of OCT4, KLF4, SOX2, and MYC. Usually only a rare subset of cells reprogram, and the molecular characteristics of this subset remain unknown. We apply retrospective clone tracing to identify and characterize the rare human fibroblasts primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis increased the reprogramming efficiency. We provide evidence for a unified model in which cells can move into and out of the primed state over time, explaining how reprogramming appears deterministic at short timescales and stochastic at long timescales. Furthermore, inhibiting the activity of LSD1 enlarged the pool of cells that were primed for reprogramming. Thus, even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Enhancing Single-Cell Western Blotting Sensitivity Using Diffusive Analyte Blotting and Antibody Conjugate Amplification

While there are many techniques to achieve highly sensitive, multiplex detection of RNA and DNA from single cells, detecting protein content often suffers from low limits of detection and throughput. Miniaturized, high-sensitivity Western blots on single cells (scWesterns) are attractive because they do not require advanced instrumentation. By physically separating analytes, scWesterns also uniquely mitigate limitations to target protein multiplexing posed by the affinity reagent performance. However, a fundamental limitation of scWesterns is their limited sensitivity for detecting low-abundance proteins, which arises from transport barriers posed by the separation gel against detection species. Here we address the sensitivity by decoupling the electrophoretic separation medium from the detection medium. We transfer scWestern separations to a nitrocellulose blotting medium with distinct mass transfer advantages over traditional in-gel probing, yielding a 5.9-fold improvement in the limit of detection. We next amplify probing of blotted proteins with enzyme-antibody conjugates, which are incompatible with traditional in-gel probing to achieve further improvement in the limit of detection to 1000 molecules, a 120-fold improvement. This enables us to detect 100% of cells in an EGFP-expressing population using fluorescently tagged and enzyme-conjugated antibodies compared to 84.5% of cells using in-gel detection. These results suggest the compatibility of nitrocellulose-immobilized scWesterns with a variety of affinity reagents─not previously accessible for in-gel use─for further signal amplification and detection of low-abundance targets.

mtFociCounter for automated single-cell mitochondrial nucleoid quantification and reproducible foci analysis

Mitochondrial DNA (mtDNA) encodes the core subunits for OXPHOS, essential in near-all eukaryotes. Packed into distinct foci (nucleoids) inside mitochondria, the number of mtDNA copies differs between cell-types and is affected in several human diseases. Currently, common protocols estimate per-cell mtDNA-molecule numbers by sequencing or qPCR from bulk samples. However, this does not allow insight into cell-to-cell heterogeneity and can mask phenotypical sub-populations. Here, we present mtFociCounter, a single-cell image analysis tool for reproducible quantification of nucleoids and other foci. mtFociCounter is a light-weight, open-source freeware and overcomes current limitations to reproducible single-cell analysis of mitochondrial foci. We demonstrate its use by analysing 2165 single fibroblasts, and observe a large cell-to-cell heterogeneity in nucleoid numbers. In addition, mtFociCounter quantifies mitochondrial content and our results show good correlation (R = 0.90) between nucleoid number and mitochondrial area, and we find nucleoid density is less variable than nucleoid numbers in wild-type cells. Finally, we demonstrate mtFociCounter readily detects differences in foci-numbers upon sample treatment, and applies to Mitochondrial RNA Granules and superresolution microscopy. mtFociCounter provides a versatile solution to reproducibly quantify cellular foci in single cells and our results highlight the importance of accounting for cell-to-cell variance and mitochondrial context in mitochondrial foci analysis.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Escherichia coli cells are primed for survival before lethal antibiotic stress

Non-genetic factors can cause significant fluctuations in gene expression levels. Regardless of growing in a stable environment, this fluctuation leads to cell-to-cell variability in an isogenic population. This phenotypic heterogeneity allows a tiny subset of bacterial cells in a population called persister cells to tolerate long-term lethal antibiotic effects by entering into a non-dividing, metabolically repressed state. We occasionally noticed a high variation in persister levels, and to explore this, we tested clonal populations starting from a single cell using a modified Luria-Delbrück fluctuation test. Although we kept the conditions same, the diversity in persistence level among clones was relatively consistent: varying from ~60- to 100- and ~40- to 70-fold for ampicillin and apramycin, respectively. Then, we divided and diluted each clone to observe whether the same clone had comparable persister levels for more than one generation. Replicates had similar persister levels even when clones were divided, diluted by 1:20, and allowed to grow for approximately five generations. This result explicitly shows a cellular memory passed on for generations and eventually lost when cells are diluted to 1:100 and regrown (>seven generations). Our result demonstrates (1) the existence of a small population prepared for stress ("primed cells") resulting in higher persister numbers; (2) the primed memory state is reproducible and transient, passed down for generations but eventually lost; and (3) a heterogeneous persister population is a result of a transiently primed reversible cell state and not due to a pre-existing genetic mutation. IMPORTANCE Antibiotics have been highly effective in treating lethal infectious diseases for almost a century. However, the increasing threat of antibiotic resistance is again causing these diseases to become life-threatening. The longer a bacteria can survive antibiotics, the more likely it is to develop resistance. Complicating matters is that non-genetic factors can allow bacterial cells with identical DNA to gain transient resistance (also known as persistence). Here, we show that a small fraction of the bacterial population called primed cells can pass down non-genetic information ("memory") to their offspring, enabling them to survive lethal antibiotics for a long time. However, this memory is eventually lost. These results demonstrate how bacteria can leverage differences among genetically identical cells formed through non-genetic factors to form primed cells with a selective advantage to survive antibiotics.

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.

Diverse clonal fates emerge upon drug treatment of homogeneous cancer cells

Even among genetically identical cancer cells, resistance to therapy frequently emerges from a small subset of those cells1-7. Molecular differences in rare individual cells in the initial population enable certain cells to become resistant to therapy7-9; however, comparatively little is known about the variability in the resistance outcomes. Here we develop and apply FateMap, a framework that combines DNA barcoding with single-cell RNA sequencing, to reveal the fates of hundreds of thousands of clones exposed to anti-cancer therapies. We show that resistant clones emerging from single-cell-derived cancer cells adopt molecularly, morphologically and functionally distinct resistant types. These resistant types are largely predetermined by molecular differences between cells before drug addition and not by extrinsic factors. Changes in the dose and type of drug can switch the resistant type of an initial cell, resulting in the generation and elimination of certain resistant types. Samples from patients show evidence for the existence of these resistant types in a clinical context. We observed diversity in resistant types across several single-cell-derived cancer cell lines and cell types treated with a variety of drugs. The diversity of resistant types as a result of the variability in intrinsic cell states may be a generic feature of responses to external cues.

Enhancing single-cell western blotting sensitivity using diffusive analyte blotting and antibody conjugate amplification

While there are many techniques to achieve highly sensitive, multiplex detection of RNA and DNA from single cells, detecting protein contents often suffers from low limits of detection and throughput. Miniaturized, high-sensitivity western blots on single cells (scWesterns) are attractive since they do not require advanced instrumentation. By physically separating analytes, scWesterns also uniquely mitigate limitations to target protein multiplexing posed by affinity reagent performance. However, a fundamental limitation of scWesterns is their limited sensitivity for detecting low-abundance proteins, which arises from transport barriers posed by the separation gel against detection species. Here we address sensitivity by decoupling the electrophoretic separation medium from the detection medium. We transfer scWestern separations to a nitrocellulose blotting medium with distinct mass transfer advantages over traditional in-gel probing, yielding a 5.9-fold improvement in limit of detection. We next amplify probing of blotted proteins with enzyme-antibody conjugates which are incompatible with traditional in-gel probing to achieve further improvement in the limit of detection to 103 molecules, a 520-fold improvement. This enables us to detect 85% and 100% of cells in an EGFP-expressing population using fluorescently tagged and enzyme-conjugated antibodies respectively, compared to 47% of cells using in-gel detection. These results suggest compatibility of nitrocellulose-immobilized scWesterns with a variety of affinity reagents - not previously accessible for in-gel use - for further signal amplification and detection of low abundance targets.

A microfluidic device for simultaneous detection of enzyme secretion and elongation of a single hypha

Filamentous fungi grow through elongation of their apical region by exocytosis and secrete enzymes that can be of commercial or industrial importance. Their hyphae exhibit extensive branching, making it difficult to control hyphal growth for observation and analysis. Therefore, although hyphal morphology and productivity are closely related, the relationship between the two has not yet been clarified. Conventional morphology and productivity studies have only compared the results of macro imaging of fungal pellets cultured in bulk with the averaged products in the culture medium. Filamentous fungi are multicellular and their expression differs between different hyphae. To truly understand the relationship between morphology and productivity, it is necessary to compare the morphology and productivity of individual hyphae. To achieve this, we developed a microfluidic system that confines hyphae to individual channels for observation and investigated the relationship between their growth, morphology, and enzyme productivity. Furthermore, using Trichoderma reesei, a potent cellulase-producing fungus, as a model, we developed a cellulase detection assay with 4-MUC substrate to detect hyphal growth and enzyme secretion in a microfluidic device in real time. Using a strain that expresses cellobiohydrolase I (CBH I) fused with AcGFP1, we compared fluorescence from the detection assay with GFP fluorescence intensity, which showed a strong correlation between the two. These results indicate that extracellular enzymes can be easily detected in the microfluidic device in real time because the production of cellulase is synchronized in T. reesei. This microfluidic system enables real-time visualization of the dynamics of hypha and enzymes during carbon source exchange and the quantitative dynamics of gene expression. This technology can be applied to many biosystems from bioenergy production to human health.

BdLT-Seq as a barcode decay-based method to unravel lineage-linked transcriptome plasticity

Cell plasticity is a core biological process underlying a myriad of molecular and cellular events taking place throughout organismal development and evolution. It has been postulated that cellular systems thrive to balance the organization of meta-stable states underlying this phenomenon, thereby maintaining a degree of populational homeostasis compatible with an ever-changing environment and, thus, life. Notably, albeit circumstantial evidence has been gathered in favour of the latter conceptual framework, a direct observation of meta-state dynamics and the biological consequences of such a process in generating non-genetic clonal diversity and divergent phenotypic output remains largely unexplored. To fill this void, here we develop a lineage-tracing technology termed Barcode decay Lineage Tracing-Seq. BdLT-Seq is based on episome-encoded molecular identifiers that, supported by the dynamic decay of the tracing information upon cell division, ascribe directionality to a cell lineage tree whilst directly coupling non-genetic molecular features to phenotypes in comparable genomic landscapes. We show that cell transcriptome states are both inherited, and dynamically reshaped following constrained rules encoded within the cell lineage in basal growth conditions, upon oncogene activation and throughout the process of reversible resistance to therapeutic cues thus adjusting phenotypic output leading to intra-clonal non-genetic diversity.

Retrospective identification of intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of the four "Yamanaka" factors OCT4, KLF4, SOX2, and MYC. However, even in homogeneous conditions, usually only a rare subset of cells admit reprogramming, and the molecular characteristics of this subset remain unknown. Here, we apply retrospective clone tracing to identify and characterize the individual human fibroblast cells that are primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis led to increased reprogramming efficiency, identifying it as a barrier to reprogramming. Changing the frequency of reprogramming by inhibiting the activity of LSD1 led to an enlarging of the pool of cells that were primed for reprogramming. Our results show that even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Gene expression model inference from snapshot RNA data using Bayesian non-parametrics

Gene expression models, which are key towards understanding cellular regulatory response, underlie observations of single-cell transcriptional dynamics. Although RNA expression data encode information on gene expression models, existing computational frameworks do not perform simultaneous Bayesian inference of gene expression models and parameters from such data. Rather, gene expression models-composed of gene states, their connectivities and associated parameters-are currently deduced by pre-specifying gene state numbers and connectivity before learning associated rate parameters. Here we propose a method to learn full distributions over gene states, state connectivities and associated rate parameters, simultaneously and self-consistently from single-molecule RNA counts. We propagate noise from fluctuating RNA counts over models by treating models themselves as random variables. We achieve this within a Bayesian non-parametric paradigm. We demonstrate our method on the Escherichia coli lacZ pathway and the Saccharomyces cerevisiae STL1 pathway, and verify its robustness on synthetic data.

Transiently heritable fates and quorum sensing drive early IFN-I response dynamics

Type I interferon (IFN-I)-mediated antiviral responses are central to host defense against viral infections. Crucial is the tight and well-orchestrated control of cellular decision-making leading to the production of IFN-Is. Innovative single-cell approaches revealed that the initiation of IFN-I production is limited to only fractions of 1-3% of the total population, both found in vitro, in vivo, and across cell types, which were thought to be stochastically regulated. To challenge this dogma, we addressed the influence of various stochastic and deterministic host-intrinsic factors on dictating early IFN-I responses, using a murine fibroblast reporter model. Epigenetic drugs influenced the percentage of responding cells. Next, with the classical Luria-Delbrück fluctuation test, we provided evidence for transient heritability driving responder fates, which was verified with mathematical modeling. Finally, while studying varying cell densities, we substantiated an important role for cell density in dictating responsiveness, similar to the phenomenon of quorum sensing. Together, this systems immunology approach opens up new avenues to progress the fundamental understanding on cellular decision-making during early IFN-I responses, which can be translated to other (immune) signaling systems.

Allelic correlation is a marker of trade-offs between barriers to transmission of expression variability and signal responsiveness in genetic networks

Genetic networks should respond to signals but prevent the transmission of spontaneous fluctuations. Limited data from mammalian cells suggest that noise transmission is uncommon, but systematic claims about noise transmission have been limited by the inability to directly measure it. Here, we build a mathematical framework modeling allelic correlation and noise transmission, showing that allelic correlation and noise transmission correspond across model parameters and network architectures. Limiting noise transmission comes with the trade-off of being unresponsive to signals, and within responsive regimes, there is a further trade-off between response time and basal noise transmission. Analysis of allele-specific single-cell RNA-sequencing data revealed that genes encoding upstream factors in signaling pathways and cell-type-specific factors have higher allelic correlation than downstream factors, suggesting they are more subject to regulation. Overall, our findings suggest that some noise transmission must result from signal responsiveness, but it can be minimized by trading off for a slower response. A record of this paper's transparent peer review process is included in the supplemental information.

Novel biological insights revealed from the investigation of multiscale genome architecture

Gene expression and cell fate determination require precise and coordinated epigenetic regulation. The complex three-dimensional (3D) genome organization plays a critical role in transcription in myriad biological processes. A wide range of architectural features of the 3D genome, including chromatin loops, topologically associated domains (TADs), chromatin compartments, and phase separation, together regulate the chromatin state and transcriptional activity at multiple levels. With the help of 3D genome informatics, recent biochemistry and imaging approaches based on different strategies have revealed functional interactions among biomacromolecules, even at the single-cell level. Here, we review the occurrence, mechanistic basis, and functional implications of dynamic genome organization, and outline recent experimental and computational approaches for profiling multiscale genome architecture to provide robust tools for studying the 3D genome.

Single-Cell Proteomics Preparation for Mass Spectrometry Analysis Using Freeze-Heat Lysis and an Isobaric Carrier

Single-cell proteomics analysis requires sensitive, quantitatively accurate, widely accessible, and robust methods. To meet these requirements, the Single-Cell ProtEomics (SCoPE2) protocol was developed as a second-generation method for quantifying hundreds to thousands of proteins from limited samples, down to the level of a single cell. Experiments using this method have achieved quantifying over 3,000 proteins across 1,500 single mammalian cells (500-1,000 proteins per cell) in 10 days of mass spectrometer instrument time. SCoPE2 leverages a freeze-heat cycle for cell lysis, obviating the need for clean-up of single cells and consequently reducing sample losses, while expediting sample preparation and simplifying its automation. Additionally, the method uses an isobaric carrier, which aids protein identification and reduces sample losses. This video protocol provides detailed guidance to enable the adoption of automated single-cell protein analysis using only equipment and reagents that are widely accessible. We demonstrate critical steps in the procedure of preparing single cells for proteomic analysis, from harvesting up to injection to liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. Additionally, viewers are guided through the principles of experimental design with the isobaric carrier, quality control for both isobaric carrier and single-cell preparations, and representative results with a discussion of limitations of the approach.

Revising immune cell coordination: Origins and importance of single-cell variation

Moving from the optimalization of single-cell technologies to the interpretation of the multi-complex single-cell data, the field of immunoengineering is granted with numerous important insights into the coordination of immune cell activation and how to modulate it for therapeutic purposes. However, insights come with additional follow-up questions that challenge our perception on how immune responses are generated and fine-tuned to fight a wide array of pathogens in ever-changing and often unpredictable microenvironments. Are immune responses really either being tightly regulated by molecular determinants, or highly flexible attributed to stochasticity? What exactly makes up the basic rules by which single cells cooperate to establish tissue-level immunity? Taking the type I IFN system and its newest insights as a main example throughout this review, we revise the basic concepts of (single) immune cell coordination, redefine the concepts of noise, stochasticity and determinism, and highlight the importance of single-cell variation in immunology and beyond.

Multiple parameters shape the 3D chromatin structure of single nuclei at the doc locus in Drosophila

The spatial organization of chromatin at the scale of topologically associating domains (TADs) and below displays large cell-to-cell variations. Up until now, how this heterogeneity in chromatin conformation is shaped by chromatin condensation, TAD insulation, and transcription has remained mostly elusive. Here, we used Hi-M, a multiplexed DNA-FISH imaging technique providing developmental timing and transcriptional status, to show that the emergence of TADs at the ensemble level partially segregates the conformational space explored by single nuclei during the early development of Drosophila embryos. Surprisingly, a substantial fraction of nuclei display strong insulation even before TADs emerge. Moreover, active transcription within a TAD leads to minor changes to the local inter- and intra-TAD chromatin conformation in single nuclei and only weakly affects insulation to the neighboring TAD. Overall, our results indicate that multiple parameters contribute to shaping the chromatin architecture of single nuclei at the TAD scale.

AP-1 transcription factor network explains diverse patterns of cellular plasticity in melanoma cells

Cellular plasticity associated with fluctuations in transcriptional programs allows individual cells in a tumor to adopt heterogeneous differentiation states and switch phenotype during their adaptive responses to therapies. Despite increasing knowledge of such transcriptional programs, the molecular basis of cellular plasticity remains poorly understood. Here, we combine multiplexed transcriptional and protein measurements at population and single-cell levels with multivariate statistical modeling to show that the state of AP-1 transcription factor network plays a unifying role in explaining diverse patterns of plasticity in melanoma. We find that a regulated balance among AP-1 factors cJUN, JUND, FRA2, FRA1, and cFOS determines the intrinsic diversity of differentiation states and adaptive responses to MAPK inhibitors in melanoma cells. Perturbing this balance through genetic depletion of specific AP-1 proteins, or by MAPK inhibitors, shifts cellular heterogeneity in a predictable fashion. Thus, AP-1 may serve as a critical node for manipulating cellular plasticity with potential therapeutic implications.

Activation domains can decouple the mean and noise of gene expression

Regulatory mechanisms set a gene's average level of expression, but a gene's expression constantly fluctuates around that average. These stochastic fluctuations, or expression noise, play a role in cell-fate transitions, bet hedging in microbes, and the development of chemotherapeutic resistance in cancer. An outstanding question is what regulatory mechanisms contribute to noise. Here, we demonstrate that, for a fixed mean level of expression, strong activation domains (ADs) at low abundance produce high expression noise, while weak ADs at high abundance generate lower expression noise. We conclude that differences in noise can be explained by the interplay between a TF's nuclear concentration and the strength of its AD's effect on mean expression, without invoking differences between classes of ADs. These results raise the possibility of engineering gene expression noise independently of mean levels in synthetic biology contexts and provide a potential mechanism for natural selection to tune the noisiness of gene expression.

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.

Sphingolipids control dermal fibroblast heterogeneity

Human cells produce thousands of lipids that change during cell differentiation and can vary across individual cells of the same type. However, we are only starting to characterize the function of these cell-to-cell differences in lipid composition. Here, we measured the lipidomes and transcriptomes of individual human dermal fibroblasts by coupling high-resolution mass spectrometry imaging with single-cell transcriptomics. We found that the cell-to-cell variations of specific lipid metabolic pathways contribute to the establishment of cell states involved in the organization of skin architecture. Sphingolipid composition is shown to define fibroblast subpopulations, with sphingolipid metabolic rewiring driving cell-state transitions. Therefore, cell-to-cell lipid heterogeneity affects the determination of cell states, adding a new regulatory component to the self-organization of multicellular systems.

Hematopoietic differentiation is characterized by a transient peak of entropy at a single-cell level

Background

Mature blood cells arise from hematopoietic stem cells in the bone marrow by a process of differentiation along one of several different lineage trajectories. This is often represented as a series of discrete steps of increasing progenitor cell commitment to a given lineage, but as for differentiation in general, whether the process is instructive or stochastic remains controversial. Here, we examine this question by analyzing single-cell transcriptomic data from human bone marrow cells, assessing cell-to-cell variability along the trajectories of hematopoietic differentiation into four different types of mature blood cells. The instructive model predicts that cells will be following the same sequence of instructions and that there will be minimal variability of gene expression between them throughout the process, while the stochastic model predicts a role for cell-to-cell variability when lineage commitments are being made.

Results

Applying Shannon entropy to measure cell-to-cell variability among human hematopoietic bone marrow cells at the same stage of differentiation, we observed a transient peak of gene expression variability occurring at characteristic points in all hematopoietic differentiation pathways. Strikingly, the genes whose cell-to-cell variation of expression fluctuated the most over the course of a given differentiation trajectory are pathway-specific genes, whereas genes which showed the greatest variation of mean expression are common to all pathways. Finally, we showed that the level of cell-to-cell variation is increased in the most immature compartment of hematopoiesis in myelodysplastic syndromes.

Conclusions

These data suggest that human hematopoietic differentiation could be better conceptualized as a dynamical stochastic process with a transient stage of cellular indetermination, and strongly support the stochastic view of differentiation. They also highlight the need to consider the role of stochastic gene expression in complex physiological processes and pathologies such as cancers, paving the way for possible noise-based therapies through epigenetic regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01264-9.

Phenotypic variability, not noise, accounts for most of the cell-to-cell heterogeneity in IFN-γ and oncostatin M signaling responses

Cellular signaling responses show substantial cell-to-cell heterogeneity, which is often ascribed to the inherent randomness of biochemical reactions, termed molecular noise, wherein high noise implies low signaling fidelity. Alternatively, heterogeneity could arise from differences in molecular content between cells, termed molecular phenotypic variability, which does not necessarily imply imprecise signaling. The contribution of these two processes to signaling heterogeneity is unclear. Here, we fused fibroblasts to produce binuclear syncytia to distinguish noise from phenotypic variability in the analysis of cytokine signaling. We reasoned that the responses of the two nuclei within one syncytium could approximate the signaling outcomes of two cells with the same molecular content, thereby disclosing noise contribution, whereas comparison of different syncytia should reveal contribution of phenotypic variability. We found that ~90% of the variance in the primary response (which was the abundance of phosphorylated, nuclear STAT) to stimulation with the cytokines interferon-γ and oncostatin M resulted from differences in the molecular content of individual cells. Thus, our data reveal that cytokine signaling in the system used here operates in a reproducible, high-fidelity manner.

Olfactory receptor choice: a case study for gene regulation in a multi-enhancer system

The mammalian genome possesses >2000 olfactory receptor (OR) alleles regulated by 63 known OR-Enhancer elements, yet each olfactory sensory neuron (OSN) expresses only a single OR allele. Choreographed changes to OSN nuclear architecture are evidently necessary for OR expression. Additionally, the insulated organization of OR-enhancers around an OR allele is a hallmark of the chosen OR. However, the biology guiding OR choice itself is unclear. Innovations in single-cell and biophysics-based analysis of nuclear architecture are revising previous models of the nucleus to include its dynamic and probabilistic nature. In this review, we ground current knowledge of OR gene regulation in these emerging theories to speculate on mechanisms that may give rise to diverse and singular OR expression.

Multiplexed single-cell proteomics using SCoPE2

Many biological systems are composed of diverse single cells. This diversity necessitates functional and molecular single-cell analysis. Single-cell protein analysis has long relied on affinity reagents, but emerging mass-spectrometry methods (either label-free or multiplexed) have enabled quantifying >1,000 proteins per cell while simultaneously increasing the specificity of protein quantification. Here we describe the Single Cell ProtEomics (SCoPE2) protocol, which uses an isobaric carrier to enhance peptide sequence identification. Single cells are isolated by FACS or CellenONE into multiwell plates and lysed by Minimal ProteOmic sample Preparation (mPOP), and their peptides labeled by isobaric mass tags (TMT or TMTpro) for multiplexed analysis. SCoPE2 affords a cost-effective single-cell protein quantification that can be fully automated using widely available equipment and scaled to thousands of single cells. SCoPE2 uses inexpensive reagents and is applicable to any sample that can be processed to a single-cell suspension. The SCoPE2 workflow allows analyzing ~200 single cells per 24 h using only standard commercial equipment. We emphasize experimental steps and benchmarks required for achieving quantitative protein analysis.

Transcription Factor Dynamics

To predict transcription, one needs a mechanistic understanding of how the numerous required transcription factors (TFs) explore the nuclear space to find their target genes, assemble, cooperate, and compete with one another. Advances in fluorescence microscopy have made it possible to visualize real-time TF dynamics in living cells, leading to two intriguing observations: first, most TFs contact chromatin only transiently; and second, TFs can assemble into clusters through their intrinsically disordered regions. These findings suggest that highly dynamic events and spatially structured nuclear microenvironments might play key roles in transcription regulation that are not yet fully understood. The emerging model is that while some promoters directly convert TF-binding events into on/off cycles of transcription, many others apply complex regulatory layers that ultimately lead to diverse phenotypic outputs. Cracking this kinetic code is an ongoing and challenging task that is made possible by combining innovative imaging approaches with biophysical models.

Spatial Organization of Chromatin: Emergence of Chromatin Structure During Development

Nuclei are central hubs for information processing in eukaryotic cells. The need to fit large genomes into small nuclei imposes severe restrictions on genome organization and the mechanisms that drive genome-wide regulatory processes. How a disordered polymer such as chromatin, which has vast heterogeneity in its DNA and histone modification profiles, folds into discernibly consistent patterns is a fundamental question in biology. Outstanding questions include how genomes are spatially and temporally organized to regulate cellular processes with high precision and whether genome organization is causally linked to transcription regulation. The advent of next-generation sequencing, super-resolution imaging, multiplexed fluorescent in situ hybridization, and single-molecule imaging in individual living cells has caused a resurgence in efforts to understand the spatiotemporal organization of the genome. In this review, we discuss structural and mechanistic properties of genome organization at different length scales and examine changes in higher-order chromatin organization during important developmental transitions.

Inducible Population Quality Control of Engineered Bacillus subtilis for Improved N-Acetylneuraminic Acid Biosynthesis

Biosynthesis by microorganisms using renewable feedstocks is an important approach for realizing sustainable chemical manufacturing. However, cell-to-cell variation in biosynthesis capability during fermentation restricts the robustness and efficiency of bioproduction, hampering the industrialization of biosynthesis. Herein, we developed an inducible population quality control system (iPopQC) for dynamically modulating the producing and nonproducing subpopulations of engineered Bacillus subtilis, which was constructed via inducible promoter- and metabolite-responsive biosensor-based genetic circuit for regulating essential genes. Moreover, iPopQC achieved a 1.97-fold increase in N-acetylneuraminic acid (NeuAc) titer by enriching producing cell subpopulation during cultivation, representing 52% higher than that of previous PopQC. Strains with double-output iPopQC cocoupling the expression of double essential genes with NeuAc production improved production robustness further, retaining NeuAc production throughout 96 h of fermentation, upon which the strains cocoupling one essential gene expression with NeuAc production abolished the production ability.

The mechanics and dynamics of cancer cells sensing noisy 3D contact guidance

Contact guidance is a major physical cue that modulates cancer cell morphology and motility, and is directly linked to the prognosis of cancer patients. Under physiological conditions, particularly in the three-dimensional (3D) extracellular matrix (ECM), the disordered assembly of fibers presents a complex directional bias to the cells. It is unclear how cancer cells respond to these noncoherent contact guidance cues. Here we combine quantitative experiments, theoretical analysis, and computational modeling to study the morphological and migrational responses of breast cancer cells to 3D collagen ECM with varying degrees of fiber alignment. We quantify the strength of contact guidance using directional coherence of ECM fibers, and find that stronger contact guidance causes cells to polarize more strongly along the principal direction of the fibers. Interestingly, sensitivity to contact guidance is positively correlated with cell aspect ratio, with elongated cells responding more strongly to ECM alignment than rounded cells. Both experiments and simulations show that cell-ECM adhesions and actomyosin contractility modulate cell responses to contact guidance by inducing a population shift between rounded and elongated cells. We also find that cells rapidly change their morphology when navigating the ECM, and that ECM fiber coherence modulates cell transition rates between different morphological phenotypes. Taken together, we find that subcellular processes that integrate conflicting mechanical cues determine cell morphology, which predicts the polarization and migration dynamics of cancer cells in 3D ECM.

A mathematical modelling framework for the regulation of intra-cellular OCT4 in human pluripotent stem cells

Human pluripotent stem cells (hPSCs) have the potential to differentiate into all cell types, a property known as pluripotency. A deeper understanding of how pluripotency is regulated is required to assist in controlling pluripotency and differentiation trajectories experimentally. Mathematical modelling provides a non-invasive tool through which to explore, characterise and replicate the regulation of pluripotency and the consequences on cell fate. Here we use experimental data of the expression of the pluripotency transcription factor OCT4 in a growing hPSC colony to develop and evaluate mathematical models for temporal pluripotency regulation. We consider fractional Brownian motion and the stochastic logistic equation and explore the effects of both additive and multiplicative noise. We illustrate the use of time-dependent carrying capacities and the introduction of Allee effects to the stochastic logistic equation to describe cell differentiation. We conclude both methods adequately capture the decline in OCT4 upon differentiation, but the Allee effect model has the advantage of allowing differentiation to occur stochastically in a sub-set of cells. This mathematical framework for describing intra-cellular OCT4 regulation can be extended to other transcription factors and developed into predictive models.

Fractional response analysis reveals logarithmic cytokine responses in cellular populations

Although we can now measure single-cell signaling responses with multivariate, high-throughput techniques our ability to interpret such measurements is still limited. Even interpretation of dose–response based on single-cell data is not straightforward: signaling responses can differ significantly between cells, encompass multiple signaling effectors, and have dynamic character. Here, we use probabilistic modeling and information-theory to introduce fractional response analysis (FRA), which quantifies changes in fractions of cells with given response levels. FRA can be universally performed for heterogeneous, multivariate, and dynamic measurements and, as we demonstrate, quantifies otherwise hidden patterns in single-cell data. In particular, we show that fractional responses to type I interferon in human peripheral blood mononuclear cells are very similar across different cell types, despite significant differences in mean or median responses and degrees of cell-to-cell heterogeneity. Further, we demonstrate that fractional responses to cytokines scale linearly with the log of the cytokine dose, which uncovers that heterogeneous cellular populations are sensitive to fold-changes in the dose, as opposed to additive changes.

Our ability to interpret single-cell multivariate signaling responses is still limited. Here the authors introduce fractional response analysis (FRA), involving fractional cell counting, capable of deconvoluting heterogeneous multivariate responses of cellular populations.

Cell-to-Cell Heterogeneity in Trypanosomes

Recent developments in single-cell and single-molecule techniques have revealed surprising levels of heterogeneity among isogenic cells. These advances have transformed the study of cell-to-cell heterogeneity into a major area of biomedical research, revealing that it can confer essential advantages, such as priming populations of unicellular organisms for future environmental stresses. Protozoan parasites, such as trypanosomes, face multiple and often hostile environments, and to survive, they undergo multiple changes, including changes in morphology, gene expression, and metabolism. But why does only a subset of proliferative cells differentiate to the next life cycle stage? Why do only some bloodstream parasites undergo antigenic switching while others stably express one variant surface glycoprotein? And why do some parasites invade an organ while others remain in the bloodstream? Building on extensive research performed in bacteria, here we suggest that biological noise can contribute to the fitness of eukaryotic pathogens and discuss the importance of cell-to-cell heterogeneity in trypanosome infections. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Variability within rare cell states enables multiple paths toward drug resistance

Molecular differences between individual cells can lead to dramatic differences in cell fate, such as death versus survival of cancer cells upon drug treatment. These originating differences remain largely hidden due to difficulties in determining precisely what variable molecular features lead to which cellular fates. Thus, we developed Rewind, a methodology that combines genetic barcoding with RNA fluorescence in situ hybridization to directly capture rare cells that give rise to cellular behaviors of interest. Applying Rewind to BRAFV600E melanoma, we trace drug-resistant cell fates back to single-cell gene expression differences in their drug-naive precursors (initial frequency of ~1:1,000-1:10,000 cells) and relative persistence of MAP kinase signaling soon after drug treatment. Within this rare subpopulation, we uncover a rich substructure in which molecular differences among several distinct subpopulations predict future differences in phenotypic behavior, such as proliferative capacity of distinct resistant clones after drug treatment. Our results reveal hidden, rare-cell variability that underlies a range of latent phenotypic outcomes upon drug exposure.

Gene duplication and subsequent diversification strongly affect phenotypic evolvability and robustness

We study the effects of non-determinism and gene duplication on the structure of genotype-phenotype (GP) maps by introducing a non-deterministic version of the Polyomino self-assembly model. This model has previously been used in a variety of contexts to model the assembly and evolution of protein quaternary structure. Firstly, we show the limit of the current deterministic paradigm which leads to built-in anti-correlation between evolvability and robustness at the genotypic level. We develop a set of metrics to measure structural properties of GP maps in a non-deterministic setting and use them to evaluate the effects of gene duplication and subsequent diversification. Our generalized versions of evolvability and robustness exhibit positive correlation for a subset of genotypes. This positive correlation is only possible because non-deterministic phenotypes can contribute to both robustness and evolvability. Secondly, we show that duplication increases robustness and reduces evolvability initially, but that the subsequent diversification that duplication enables has a stronger, inverse effect, greatly increasing evolvability and reducing robustness relative to their original values.

Forecasting cellular states: from descriptive to predictive biology via single-cell multiomics

As the single cell field races to characterize each cell type, state, and behavior, the complexity of the computational analysis approaches the complexity of the biological systems. Single cell and imaging technologies now enable unprecedented measurements of state transitions in biological systems, providing high-throughput data that capture tens-of-thousands of measurements on hundreds-of-thousands of samples. Thus, the definition of cell type and state is evolving to encompass the broad range of biological questions now attainable. To answer these questions requires the development of computational tools for integrated multi-omics analysis. Merged with mathematical models, these algorithms will be able to forecast future states of biological systems, going from statistical inferences of phenotypes to time course predictions of the biological systems with dynamic maps analogous to weather systems. Thus, systems biology for forecasting biological system dynamics from multi-omic data represents the future of cell biology empowering a new generation of technology-driven predictive medicine.

From Bench to Keyboard and Back Again: A Brief History of Lambda Phage Modeling

Cellular decision making is the process whereby cells choose one developmental pathway from multiple possible ones, either spontaneously or due to environmental stimuli. Examples in various cell types suggest an almost inexhaustible plethora of underlying molecular mechanisms. In general, cellular decisions rely on the gene regulatory network, which integrates external signals to drive cell fate choice. The search for general principles of such a process benefits from appropriate biological model systems that reveal how and why certain gene regulatory mechanisms drive specific cellular decisions according to ecological context and evolutionary outcomes. In this article, we review the historical and ongoing development of the phage lambda lysis-lysogeny decision as a model system to investigate all aspects of cellular decision making. The unique generality, simplicity, and richness of phage lambda decision making render it a constant source ofmathematical modeling-aided inspiration across all of biology. We discuss the origins and progress of quantitative phage lambda modeling from the 1950s until today, as well as its possible future directions. We provide examples of how modeling enabled methods and theory development, leading to new biological insights by revealing gaps in the theory and pinpointing areas requiring further experimental investigation. Overall, we highlight the utility of theoretical approaches both as predictive tools, to forecast the outcome of novel experiments, and as explanatory tools, to elucidate the natural processes underlying experimental data.

A role for the fusogen eff-1 in epidermal stem cell number robustness in Caenorhabditis elegans

Developmental patterning in Caenorhabditis elegans is known to proceed in a highly stereotypical manner, which raises the question of how developmental robustness is achieved despite the inevitable stochastic noise. We focus here on a population of epidermal cells, the seam cells, which show stem cell-like behaviour and divide symmetrically and asymmetrically over post-embryonic development to generate epidermal and neuronal tissues. We have conducted a mutagenesis screen to identify mutants that introduce phenotypic variability in the normally invariant seam cell population. We report here that a null mutation in the fusogen eff-1 increases seam cell number variability. Using time-lapse microscopy and single molecule fluorescence hybridisation, we find that seam cell division and differentiation patterns are mostly unperturbed in eff-1 mutants, indicating that cell fusion is uncoupled from the cell differentiation programme. Nevertheless, seam cell losses due to the inappropriate differentiation of both daughter cells following division, as well as seam cell gains through symmetric divisions towards the seam cell fate were observed at low frequency. We show that these stochastic errors likely arise through accumulation of defects interrupting the continuity of the seam and changing seam cell shape, highlighting the role of tissue homeostasis in suppressing phenotypic variability during development.

Genome-wide impact of hydrogen peroxide on maintenance DNA methylation in replicating cells

Background

Environmental factors, such as oxidative stress, have the potential to modify the epigenetic landscape of cells. We have previously shown that DNA methyltransferase (DNMT) activity can be inhibited by sublethal doses of hydrogen peroxide (H₂O₂). However, site-specific changes in DNA methylation and the reversibility of any changes have not been explored. Using bead chip array technology, differential methylation was assessed in Jurkat T-lymphoma cells following exposure to H₂O₂.

Results

Sublethal H₂O₂ exposure was associated with an initial genome-wide decrease in DNA methylation in replicating cells, which was largely corrected 72 h later. However, some alterations were conserved through subsequent cycles of cell division. Significant changes to the variability of DNA methylation were also observed both globally and at the site-specific level.

Conclusions

This research indicates that increased exposure to H₂O₂ can result in long-term alterations to DNA methylation patterns, providing a mechanism for environmental factors to have prolonged impact on gene expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00388-6.