Single-cell variability in multicellular organisms

DMF-Bimol: Counting mRNA and Protein Molecules in Single Cells with Digital Microfluidics

Analyzing single-cell protein and mRNA levels yields invaluable insights into cellular functions and the intricacies of biologically heterogeneous systems. Current joint mRNAs and protein analysis methodologies suffer from relative quantification, low sensitivity, possible background interference, and tedious manual manipulation. Therefore, we propose DMF-Bimol that leverages addressable digital microfluidics to automate digital counting of single-cell mRNA and protein based on proximity ligation assay (PLA) and one-step RT-droplet digital PCR (RT-ddPCR). Through an engineered hydrophilic-hydrophobic interface, DMF-Bimol enables efficient single-cell isolation and lossless protein and nucleic acid processing. The closed droplet reaction system enhances the protein concentration and isolates exogenous contaminants, thereby dramatically improving the efficiency of the PLA reaction. The limit of detection of this approach achieves 3313 protein copies, marking a significant 17-fold enhancement in sensitivity over traditional benchtop PLA. This heightened sensitivity also uncovers a lower correlation between mRNA and protein levels in individual cells (Spearman r = 0.255) than bulk results, reflecting the complex relationship in heterogeneous cells. Using DMF-Bimol, we observed a significant upsurge of CD147 protein in CD138+ myeloma cells but consistent levels of CD147 mRNAs compared with normal leukocytes. This discovery indicates a possible consequence of CD147 oncogenic activation that tends to harness protein translation to bolster tumor cell survival and enhance invasiveness, highlighting the potential of DMF-Bimol in unveiling intricate dynamics in translation processes at the single-cell level.

Efficient and scalable prediction of stochastic reaction-diffusion processes using graph neural networks

The dynamics of locally interacting particles that are distributed in space give rise to a multitude of complex behaviours. However the simulation of reaction-diffusion processes which model such systems is highly computationally expensive, the cost increasing rapidly with the size of space. Here, we devise a graph neural network based approach that uses cheap Monte Carlo simulations of reaction-diffusion processes in a small space to cast predictions of the dynamics of the same processes in a much larger and complex space, including spaces modelled by networks with heterogeneous topology. By applying the method to two biological examples, we show that it leads to accurate results in a small fraction of the computation time of standard stochastic simulation methods. The scalability and accuracy of the method suggest it is a promising approach for studying reaction-diffusion processes in complex spatial domains such as those modelling biochemical reactions, population evolution and epidemic spreading.

Proto-neural networks from thermal proteins

Proteinoids are synthetic polymers that have structural similarities to natural proteins, and their formation is achieved through the application of heat to amino acid combinations in a dehydrated environment. The thermal proteins, initially synthesised by Sidney Fox during the 1960s, has the ability to undergo self-assembly, resulting in the formation of microspheres that resemble cells. These microspheres have fascinating biomimetic characteristics. In recent studies, substantial advancements have been made in elucidating the electrical signalling phenomena shown by proteinoids, hence showcasing their promising prospects in the field of neuro-inspired computing. This study demonstrates the advancement of experimental prototypes that employ proteinoids in the construction of fundamental neural network structures. The article provides an overview of significant achievements in proteinoid systems, such as the demonstration of electrical excitability, emulation of synaptic functions, capabilities in pattern recognition, and adaptability of network structures. This study examines the similarities and differences between proteinoid networks and spontaneous neural computation. We examine the persistent challenges associated with deciphering the underlying mechanisms of emergent proteinoid-based intelligence. Additionally, we explore the potential for developing bio-inspired computing systems using synthetic thermal proteins in forthcoming times. The results of this study offer a theoretical foundation for the advancement of adaptive, self-assembling electronic systems that operate using artificial bio-neural principles.

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.

Resource allocation in mammalian systems

Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.

Studying temporal dynamics of single cells: expression, lineage and regulatory networks

Learning how multicellular organs are developed from single cells to different cell types is a fundamental problem in biology. With the high-throughput scRNA-seq technology, computational methods have been developed to reveal the temporal dynamics of single cells from transcriptomic data, from phenomena on cell trajectories to the underlying mechanism that formed the trajectory. There are several distinct families of computational methods including Trajectory Inference (TI), Lineage Tracing (LT), and Gene Regulatory Network (GRN) Inference which are involved in such studies. This review summarizes these computational approaches which use scRNA-seq data to study cell differentiation and cell fate specification as well as the advantages and limitations of different methods. We further discuss how GRNs can potentially affect cell fate decisions and trajectory structures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-023-01090-5.

Inferring delays in partially observed gene regulation processes

MOTIVATION

Cell function is regulated by gene regulatory networks (GRNs) defined by protein-mediated interaction between constituent genes. Despite advances in experimental techniques, we can still measure only a fraction of the processes that govern GRN dynamics. To infer the properties of GRNs using partial observation, unobserved sequential processes can be replaced with distributed time delays, yielding non-Markovian models. Inference methods based on the resulting model suffer from the curse of dimensionality.

RESULTS

We develop a simulation-based Bayesian MCMC method employing an approximate likelihood for the efficient and accurate inference of GRN parameters when only some of their products are observed. We illustrate our approach using a two-step activation model: an activation signal leads to the accumulation of an unobserved regulatory protein, which triggers the expression of observed fluorescent proteins. With prior information about observed fluorescent protein synthesis, our method successfully infers the dynamics of the unobserved regulatory protein. We can estimate the delay and kinetic parameters characterizing target regulation including transcription, translation, and target searching of an unobserved protein from experimental measurements of the products of its target gene. Our method is scalable and can be used to analyze non-Markovian models with hidden components.

AVAILABILITY AND IMPLEMENTATION

Our code is implemented in R and is freely available with a simple example data at https://github.com/Mathbiomed/SimMCMC.

Concentration fluctuations in growing and dividing cells: Insights into the emergence of concentration homeostasis

Intracellular reaction rates depend on concentrations and hence their levels are often regulated. However classical models of stochastic gene expression lack a cell size description and cannot be used to predict noise in concentrations. Here, we construct a model of gene product dynamics that includes a description of cell growth, cell division, size-dependent gene expression, gene dosage compensation, and size control mechanisms that can vary with the cell cycle phase. We obtain expressions for the approximate distributions and power spectra of concentration fluctuations which lead to insight into the emergence of concentration homeostasis. We find that (i) the conditions necessary to suppress cell division-induced concentration oscillations are difficult to achieve; (ii) mRNA concentration and number distributions can have different number of modes; (iii) two-layer size control strategies such as sizer-timer or adder-timer are ideal because they maintain constant mean concentrations whilst minimising concentration noise; (iv) accurate concentration homeostasis requires a fine tuning of dosage compensation, replication timing, and size-dependent gene expression; (v) deviations from perfect concentration homeostasis show up as deviations of the concentration distribution from a gamma distribution. Some of these predictions are confirmed using data for E. coli, fission yeast, and budding yeast.

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.

Noise reduction by upstream open reading frames

Gene expression is prone to burst production, making it a highly noisy process that requires additional controls. Upstream open reading frames (uORFs) are widely present in the 5' leader sequences of 30-50% of eukaryotic messenger RNAs^1-3. The translation of uORFs can repress the translation efficiency of the downstream main coding sequences. Whether the low translation efficiency leads to a different variation, or noise, in gene expression has not been investigated, nor has the direct biological impact of uORF-repressed translation. Here we show that uORFs achieve low but precise protein production in plant cells, possibly by reducing the protein production rate. We also demonstrate that, by buffering a stable TIMING OF CAB EXPRESSION 1 (TOC1) protein production level, uORFs contribute to the robust operation of the plant circadian clock. Our results provide both an action model and the biological impact of uORFs in translational control to mitigate transcriptional noise for precise protein production.

The macroscopic limit to synchronization of cellular clocks in single cells of Neurospora crassa

We determined the macroscopic limit for phase synchronization of cellular clocks in an artificial tissue created by a “big chamber” microfluidic device to be about 150,000 cells or less. The dimensions of the microfluidic chamber allowed us to calculate an upper limit on the radius of a hypothesized quorum sensing signal molecule of 13.05 nm using a diffusion approximation for signal travel within the device. The use of a second microwell microfluidic device allowed the refinement of the macroscopic limit to a cell density of 2166 cells per fixed area of the device for phase synchronization. The measurement of averages over single cell trajectories in the microwell device supported a deterministic quorum sensing model identified by ensemble methods for clock phase synchronization. A strong inference framework was used to test the communication mechanism in phase synchronization of quorum sensing versus cell-to-cell contact, suggesting support for quorum sensing. Further evidence came from showing phase synchronization was density-dependent.

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.

A single-cell resolved cell-cell communication model explains lineage commitment in hematopoiesis

Cells do not make fate decisions independently. Arguably, every cell-fate decision occurs in response to environmental signals. In many cases, cell-cell communication alters the dynamics of the internal gene regulatory network of a cell to initiate cell-fate transitions, yet models rarely take this into account. Here, we have developed a multiscale perspective to study the granulocyte-monocyte versus megakaryocyte-erythrocyte fate decisions. This transition is dictated by the GATA1-PU.1 network: a classical example of a bistable cell-fate system. We show that, for a wide range of cell communication topologies, even subtle changes in signaling can have pronounced effects on cell-fate decisions. We go on to show how cell-cell coupling through signaling can spontaneously break the symmetry of a homogenous cell population. Noise, both intrinsic and extrinsic, shapes the decision landscape profoundly, and affects the transcriptional dynamics underlying this important hematopoietic cell-fate decision-making system.

This article has an associated ‘The people behind the papers’ interview.

Summary: Through theory and computational modeling, cell-cell communication is revealed to be a crucial and under-appreciated determinant of cell-fate decision-making during hematopoiesis.

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

Plant 3D genomics: the exploration and application of chromatin organization

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

Identifiability analysis for stochastic differential equation models in systems biology

Mathematical models are routinely calibrated to experimental data, with goals ranging from building predictive models to quantifying parameters that cannot be measured. Whether or not reliable parameter estimates are obtainable from the available data can easily be overlooked. Such issues of parameter identifiability have important ramifications for both the predictive power of a model, and the mechanistic insight that can be obtained. Identifiability analysis is well-established for deterministic, ordinary differential equation (ODE) models, but there are no commonly adopted methods for analysing identifiability in stochastic models. We provide an accessible introduction to identifiability analysis and demonstrate how existing ideas for analysis of ODE models can be applied to stochastic differential equation (SDE) models through four practical case studies. To assess structural identifiability, we study ODEs that describe the statistical moments of the stochastic process using open-source software tools. Using practically motivated synthetic data and Markov chain Monte Carlo methods, we assess parameter identifiability in the context of available data. Our analysis shows that SDE models can often extract more information about parameters than deterministic descriptions. All code used to perform the analysis is available on Github.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

Stochastic reaction networks in dynamic compartment populations

Many biochemical processes in living systems take place in compartmentalized environments, where individual compartments can interact with each other and undergo dynamic remodeling. Studying such processes through mathematical models poses formidable challenges because the underlying dynamics involve a large number of states, which evolve stochastically with time. Here we propose a mathematical framework to study stochastic biochemical networks in compartmentalized environments. We develop a generic population model, which tracks individual compartments and their molecular composition. We then show how the time evolution of this system can be studied effectively through a small number of differential equations, which track the statistics of the population. Our approach is versatile and renders an important class of biological systems computationally accessible.

Compartmentalization of biochemical processes underlies all biological systems, from the organelle to the tissue scale. Theoretical models to study the interplay between noisy reaction dynamics and compartmentalization are sparse, and typically very challenging to analyze computationally. Recent studies have made progress toward addressing this problem in the context of specific biological systems, but a general and sufficiently effective approach remains lacking. In this work, we propose a mathematical framework based on counting processes that allows us to study dynamic compartment populations with arbitrary interactions and internal biochemistry. We derive an efficient description of the dynamics in terms of differential equations which capture the statistics of the population. We demonstrate the relevance of our approach by analyzing models inspired by different biological processes, including subcellular compartmentalization and tissue homeostasis.

Peptide Binding Sites of Connexin Proteins

Intercellular gap junction (GJ) contacts formed by the coupling of connexin (Cx) hemichannels (HCs) embedded into the plasma membranes of neighboring cells play significant role in the development, signaling and malfunctions of mammalian tissues. Understanding and targeting GJ functions, however, calls for finding valid Cx subtype-specific inhibitors. We conjecture the lack of information about binding interactions between the GJ interface forming extracellular EL1 and EL2 loops and peptide mimetics designed to specifically inhibit Cx43HC coupling to Cx43GJ. Here, we explore active spots at the GJ interface using known peptide inhibitors that mimic various segments of EL1 and EL2. Binding interactions of these peptide inhibitors and the non-peptide inhibitor quinine has been modelled in combination with the use of blind docking molecular mechanics (MM). The neuron-specific Cx36HC and astrocyte-specific Cx43HC subtypes were modelled with a template derived from the high-resolution structure of Cx26GJ. GJ-coupled and free Cx36HC and Cx43HC models were obtained by dissection of GJs (GJ-coupled) followed by 50 ns molecular dynamics (free). Molecular mechanics (MM) calculations were performed by the docking of inhibitors, explicitly the designed Cx43 EL1 or EL2 loop sequence mimetics (GAP26, P5 or P180–195, GAP27, Peptide5, respectively) and the Cx36 subtype-specific quinine into the model structures. In order to explore specific binding interactions between inhibitors and CxHC subtypes, MM/Generalized Born Surface Area (MM/GBSA) ΔGbind values for representative conformers of peptide mimetics and quinine were evaluated by mapping the binding surface of Cx36HC and Cx43HC for all inhibitors. Quinine specifically contacts Cx36 EL1 residues V54-C55-N56-T57-L58, P60 and N63. Blocking the vestibule by the side of Cx36HC entry, quinine explicitly interacts with the non-conserved V54, L58, N63 residues of Cx36 EL1. In addition, our work challenges the predicted specificity of peptide mimetics, showing that the docking site of peptides is unrelated to the location of the sequence they mimic. Binding features, such as unaffected EL2 residues and the lack of Cx43 subtype-specificity of peptide mimetics, suggest critical roles for peptide stringency and dimension, possibly pertaining to the Cx subtype-specificity of peptide inhibitors.

Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art

Stochasticity is a key characteristic of intracellular processes such as gene regulation and chemical signalling. Therefore, characterizing stochastic effects in biochemical systems is essential to understand the complex dynamics of living things. Mathematical idealizations of biochemically reacting systems must be able to capture stochastic phenomena. While robust theory exists to describe such stochastic models, the computational challenges in exploring these models can be a significant burden in practice since realistic models are analytically intractable. Determining the expected behaviour and variability of a stochastic biochemical reaction network requires many probabilistic simulations of its evolution. Using a biochemical reaction network model to assist in the interpretation of time-course data from a biological experiment is an even greater challenge due to the intractability of the likelihood function for determining observation probabilities. These computational challenges have been subjects of active research for over four decades. In this review, we present an accessible discussion of the major historical developments and state-of-the-art computational techniques relevant to simulation and inference problems for stochastic biochemical reaction network models. Detailed algorithms for particularly important methods are described and complemented with Matlab^® implementations. As a result, this review provides a practical and accessible introduction to computational methods for stochastic models within the life sciences community.

Coupled membrane transporters reduce noise

Molecular systems are inherently probabilistic and operate in a noisy environment, yet, despite all these uncertainties, molecular functions are surprisingly reliable and robust. The principles used by natural systems to deal with noise are still not well understood, especially in a nonhomogeneous environment where molecules can diffuse across different compartments. In this paper we show that membrane transport mechanisms have very effective properties of noise reduction. In particular, we show that active transport mechanisms (those that can transport against a gradient of concentration by using energy or by means of the concentration gradient of other substances), such as symporters and antiporters, have surprising efficiency in noise reduction, which outperforms passive diffusion mechanisms and are well below Poisson levels. We link our results to the coupled transport of potassium, sodium, and glucose to show that the noise in internal glucose level can be greatly reduced. Our results show that compartmentalization can be a highly effective mechanism of noise reduction and suggests that membrane transport could give this extra benefit, contributing to the emergence of complex compartmentalization in eukaryotes.

Learning unsupervised feature representations for single cell microscopy images with paired cell inpainting

Cellular microscopy images contain rich insights about biology. To extract this information, researchers use features, or measurements of the patterns of interest in the images. Here, we introduce a convolutional neural network (CNN) to automatically design features for fluorescence microscopy. We use a self-supervised method to learn feature representations of single cells in microscopy images without labelled training data. We train CNNs on a simple task that leverages the inherent structure of microscopy images and controls for variation in cell morphology and imaging: given one cell from an image, the CNN is asked to predict the fluorescence pattern in a second different cell from the same image. We show that our method learns high-quality features that describe protein expression patterns in single cells both yeast and human microscopy datasets. Moreover, we demonstrate that our features are useful for exploratory biological analysis, by capturing high-resolution cellular components in a proteome-wide cluster analysis of human proteins, and by quantifying multi-localized proteins and single-cell variability. We believe paired cell inpainting is a generalizable method to obtain feature representations of single cells in multichannel microscopy images.

To understand the cell biology captured by microscopy images, researchers use features, or measurements of relevant properties of cells, such as the shape or size of cells, or the intensity of fluorescent markers. Features are the starting point of most image analysis pipelines, so their quality in representing cells is fundamental to the success of an analysis. Classically, researchers have relied on features manually defined by imaging experts. In contrast, deep learning techniques based on convolutional neural networks (CNNs) automatically learn features, which can outperform manually-defined features at image analysis tasks. However, most CNN methods require large manually-annotated training datasets to learn useful features, limiting their practical application. Here, we developed a new CNN method that learns high-quality features for single cells in microscopy images, without the need for any labeled training data. We show that our features surpass other comparable features in identifying protein localization from images, and that our method can generalize to diverse datasets. By exploiting our method, researchers will be able to automatically obtain high-quality features customized to their own image datasets, facilitating many downstream analyses, as we highlight by demonstrating many possible use cases of our features in this study.

Analyzing the stability of gene expression using a simple reaction-diffusion model in an early Drosophila embryo

In all complex organisms, the precise levels and timing of gene expression controls vital biological processes. In higher eukaryotes, including the fruit fly Drosophila melanogaster, the complex molecular control of transcription (the synthesis of RNA from DNA) and translation (the synthesis of proteins from RNA) events driving this gene expression are not fully understood. In particular, for Drosophila melanogaster, there is a plethora of experimental data, including quantitative measurements of both RNA and protein concentrations, but the precise mechanisms that control the dynamics of gene expression during early development and the processes which lead to steady-state levels of certain proteins remain elusive. This study analyzes a current mathematical modeling approach in an attempt to better understand the long-term behavior of gene regulation. The model is a modified reaction-diffusion equation which has been previously employed in predicting gene expression levels and studying the relative contributions of transcription and translation events to protein abundance [10,11,24]. Here, we use Matrix Algebra and Analysis techniques to study the stability of the gene expression system and analyze equilibria, using very general assumptions regarding the parameter values incorporated into the model. We prove that, given realistic biological parameter values, the system will result in a unique, stable equilibrium solution. Additionally, we give an example of this long-term behavior using the model alongside actual experimental data obtained from Drosophila embryos.

Spatial Stochastic Intracellular Kinetics: A Review of Modelling Approaches

Models of chemical kinetics that incorporate both stochasticity and diffusion are an increasingly common tool for studying biology. The variety of competing models is vast, but two stand out by virtue of their popularity: the reaction-diffusion master equation and Brownian dynamics. In this review, we critically address a number of open questions surrounding these models: How can they be justified physically? How do they relate to each other? How do they fit into the wider landscape of chemical models, ranging from the rate equations to molecular dynamics? This review assumes no prior knowledge of modelling chemical kinetics and should be accessible to a wide range of readers.

Bayesian inference of distributed time delay in transcriptional and translational regulation

MOTIVATION

Advances in experimental and imaging techniques have allowed for unprecedented insights into the dynamical processes within individual cells. However, many facets of intracellular dynamics remain hidden, or can be measured only indirectly. This makes it challenging to reconstruct the regulatory networks that govern the biochemical processes underlying various cell functions. Current estimation techniques for inferring reaction rates frequently rely on marginalization over unobserved processes and states. Even in simple systems this approach can be computationally challenging, and can lead to large uncertainties and lack of robustness in parameter estimates. Therefore we will require alternative approaches to efficiently uncover the interactions in complex biochemical networks.

RESULTS

We propose a Bayesian inference framework based on replacing uninteresting or unobserved reactions with time delays. Although the resulting models are non-Markovian, recent results on stochastic systems with random delays allow us to rigorously obtain expressions for the likelihoods of model parameters. In turn, this allows us to extend MCMC methods to efficiently estimate reaction rates, and delay distribution parameters, from single-cell assays. We illustrate the advantages, and potential pitfalls, of the approach using a birth-death model with both synthetic and experimental data, and show that we can robustly infer model parameters using a relatively small number of measurements. We demonstrate how to do so even when only the relative molecule count within the cell is measured, as in the case of fluorescence microscopy.

AVAILABILITY AND IMPLEMENTATION

Accompanying code in R is available at https://github.com/cbskust/DDE_BD.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Application of the Transcriptional Disease Signature (TDSs) to Screen Melanoma-Effective Compounds in a Small Fish Model

Cell culture and protein target-based compound screening strategies, though broadly utilized in selecting candidate compounds, often fail to eliminate candidate compounds with non-target effects and/or safety concerns until late in the drug developmental process. Phenotype screening using intact research animals is attractive because it can help identify small molecule candidate compounds that have a high probability of proceeding to clinical use. Most FDA approved, first-in-class small molecules were identified from phenotypic screening. However, phenotypic screening using rodent models is labor intensive, low-throughput, and very expensive. As a novel alternative for small molecule screening, we have been developing gene expression disease profiles, termed the Transcriptional Disease Signature (TDS), as readout of small molecule screens for therapeutic molecules. In this concept, compounds that can reverse, or otherwise affect known disease-associated gene expression patterns in whole animals may be rapidly identified for more detailed downstream direct testing of their efficacy and mode of action. To establish proof of concept for this screening strategy, we employed a transgenic strain of a small aquarium fish, medaka (Oryzias latipes), that overexpresses the malignant melanoma driver gene xmrk, a mutant egfr gene, that is driven by a pigment cell-specific mitf promoter. In this model, melanoma develops with 100% penetrance. Using the transgenic medaka malignant melanoma model, we established a screening system that employs the NanoString nCounter platform to quantify gene expression within custom sets of TDS gene targets that we had previously shown to exhibit differential transcription among xmrk-transgenic and wild-type medaka. Compound-modulated gene expression was identified using an internet-accessible custom-built data processing pipeline. The effect of a given drug on the entire TDS profile was estimated by comparing compound-modulated genes in the TDS using an activation Z-score and Kolmogorov-Smirnov statistics. TDS gene probes were designed that target common signaling pathways that include proliferation, development, toxicity, immune function, metabolism and detoxification. These pathways may be utilized to evaluate candidate compounds for potential favorable, or unfavorable, effects on melanoma-associated gene expression. Here we present the logistics of using medaka to screen compounds, as well as, the development of a user-friendly NanoString data analysis pipeline to support feasibility of this novel TDS drug-screening strategy.

Diffusion-dynamics laws in stochastic reaction networks

Many biological activities are induced by cellular chemical reactions of diffusing reactants. The dynamics of such systems can be captured by stochastic reaction networks. A recent numerical study has shown that diffusion can significantly enhance the fluctuations in gene regulatory networks. However, the universal relation between diffusion and stochastic system dynamics remains veiled. Within the approximation of reaction-diffusion master equation (RDME), we find general relation that the steady-state distribution in complex balanced networks is diffusion-independent. Here, complex balance is the nonequilibrium generalization of detailed balance. We also find that for a diffusion-included network with a Poisson-like steady-state distribution, the diffusion can be ignored at steady state. We then derive a necessary and sufficient condition for networks holding such steady-state distributions. Moreover, we show that for linear reaction networks the RDME reduces to the chemical master equation, which implies that the stochastic dynamics of networks is unaffected by diffusion at any arbitrary time. Our findings shed light on the fundamental question of when diffusion can be neglected, or (if nonnegligible) its effects on the stochastic dynamics of the reaction network.

Recent Advances and New Perspectives in Capillary Electrophoresis-Mass Spectrometry for Single Cell "Omics"

Accurate clinical therapeutics rely on understanding the metabolic responses of individual cells. However, the high level of heterogeneity between cells means that simply sampling from large populations of cells is not necessarily a reliable approximation of an individual cell's response. As a result, there have been numerous developments in the field of single-cell analysis to address this lack of knowledge. Many of these developments have focused on the coupling of capillary electrophoresis (CE), a separation technique with low sample consumption and high resolving power, and mass spectrometry (MS), a sensitive detection method for interrogating all ions in a sample in a single analysis. In recent years, there have been many notable advancements at each step of the single-cell CE-MS analysis workflow, including sampling, manipulation, separation, and MS analysis. In each of these areas, the combined improvements in analytical instrumentation and achievements of numerous researchers have served to drive the field forward to new frontiers. Consequently, notable biological discoveries have been made possible by the implementation of these methods. Although there is still room in the field for numerous further advances, researchers have effectively minimized various limitations in detection of analytes, and it is expected that there will be many more developments in the near future.

Role of noise and parametric variation in the dynamics of gene regulatory circuits

Stochasticity in gene expression impacts the dynamics and functions of gene regulatory circuits. Intrinsic noises, including those that are caused by low copy number of molecules and transcriptional bursting, are usually studied by stochastic simulations. However, the role of extrinsic factors, such as cell-to-cell variability and heterogeneity in the microenvironment, is still elusive. To evaluate the effects of both the intrinsic and extrinsic noises, we develop a method, named sRACIPE, by integrating stochastic analysis with random circuit perturbation (RACIPE) method. RACIPE uniquely generates and analyzes an ensemble of models with random kinetic parameters. Previously, we have shown that the gene expression from random models form robust and functionally related clusters. In sRACIPE we further develop two stochastic simulation schemes, aiming to reduce the computational cost without sacrificing the convergence of statistics. One scheme uses constant noise to capture the basins of attraction, and the other one uses simulated annealing to detect the stability of states. By testing the methods on several synthetic gene regulatory circuits and an epithelial-mesenchymal transition network in squamous cell carcinoma, we demonstrate that sRACIPE can interpret the experimental observations from single-cell gene expression data. We observe that parametric variation (the spread of parameters around a median value) increases the spread of the gene expression clusters, whereas high noise merges the states. Our approach quantifies the robustness of a gene circuit in the presence of noise and sheds light on a new mechanism of noise-induced hybrid states. We have implemented sRACIPE as an R package.

Systems, variation, individuality and plant hormones

Inter-individual variation in plants and particularly in hormone content, figures strongly in evolution and behaviour. Homo sapiens and Arabidopsis exhibit similar and substantial phenotypic and molecular variation. Whereas there is a very substantial degree of hormone variation in mankind, reports of inter-individual variation in plant hormone content are virtually absent but are likely to be as large if not larger than that in mankind. Reasons for this absence are discussed. Using an example of inter-individual variation in ethylene content in ripening, the article shows how biological time is compressed by hormones. It further resolves an old issue of very wide hormone dose response that result directly from negative regulation in hormone (and light) transduction. Negative regulation is used because of inter-individual variability in hormone synthesis, receptors and ancillary proteins, a consequence of substantial genomic and environmental variation. Somatic mosaics have been reported for several plant tissues and these too contribute to tissue variation and wide variation in hormone response. The article concludes by examining what variation exists in gravitropic responses. There are multiple sensing systems of gravity vectors and multiple routes towards curvature. These are an aspect of the need for reliability in both inter-individual variation and unpredictable environments. Plant hormone inter-individuality is a new area for research and is likely to change appreciation of the mechanisms that underpin individual behaviour.

Measurement of Cell-Matrix Adhesion at Single-Cell Resolution for Revealing the Functions of Biomaterials for Adherent Cell Culture

Cell adhesion is essential for a cell to maintain its functions, and biomaterials acting as the extracellular matrix (ECM) play a vital role. However, conventional methods for evaluating the functions of biomaterials become insufficient and sometimes incorrect when we give a deeper insight into single-cell research. In this work, we reported a novel methodology for the measurement of cell-matrix adhesion at single-cell resolution that could precisely evaluate the functions of biomaterials for adherent cell culture. A microfludic device, a live single-cell extractor (LSCE), was used for cell extraction. We applied this method to evaluate various modified biomaterials. The results indicated that poly(l-polylysine) (PLL)-coated glass and fibronection (FN)-coated glass slides showed the best biocompatibility for adherent cell culture following by the (3-aminopropyl)triethoxysilane (APTES)-coated glass, while piranha solution treated glass slide and octadecyltrichlorosilane (OTS)-coated glass showed weak biocompatibilities. Furthermore, APTES, PLL, and FN modifications enhanced the cell heterogeneity, while the OTS modification weakened the cell heterogeneity compare to the initial piranha solution treated glass. The method not only clarified the cell-matrix adhesion strength at single-cell resolution but also revealed the influences of biomaterials on cell-matrix adhesion and heterogeneity of cell-matrix adhesion for adherent cell culture. It might be a general strategy for precise evaluation of biomaterials.