Dynamic modes of Notch transcription hubs conferring memory and stochastic activation revealed by live imaging the co-activator Mastermind.

Developmental programming involves the accurate conversion of signalling levels and dynamics to transcriptional outputs. The transcriptional relay in the Notch pathway relies on nuclear complexes containing the co-activator Mastermind (Mam). By tracking these complexes in real time, we reveal that they promote the formation of a dynamic transcription hub in Notch ON nuclei which concentrates key factors including the Mediator CDK module. The composition of the hub is labile and persists after Notch withdrawal conferring a memory that enables rapid reformation. Surprisingly, only a third of Notch ON hubs progress to a state with nascent transcription, which correlates with polymerase II and core Mediator recruitment. This probability is increased by a second signal. The discovery that target-gene transcription is probabilistic has far-reaching implications because it implies that stochastic differences in Notch pathway output can arise downstream of receptor activation.

To correctly give rise to future tissues, cells in an embryo must receive and respond to the right signals, at the right time, in the right way. This involves genes being switched on quickly, with cells often ensuring that a range of molecular actors physically come together at 'transcription hubs' in the nucleus ­ the compartment that houses genetic information. These hubs are thought to foster a microenvironment that facilitates the assembly of the machinery that will activate and copy the required genes into messenger RNA molecules. The resulting 'mRNAs' act as templates for producing the corresponding proteins, allowing cells to adequately respond to signals. For example, the activation at the cell surface of a molecule called Notch triggers a series of events that lead to important developmental genes being transcribed within minutes. This process involves a dedicated group of proteins, known as Notch nuclear complexes, quickly getting together in the nucleus and interacting with the transcriptional machinery. How they do this efficiently at the right gene locations is, however, still poorly understood. In particular, it remained unclear whether Notch nuclear complexes participate in the formation of transcription hubs, as well as how these influence mRNA production and the way cells 'remember' having been exposed to Notch activity. To investigate these questions, DeHaro-Arbona et al. genetically engineered fruit flies so that their Notch nuclear complexes and Notch target genes both carried visible tags that could be tracked in living cells in real time. Microscopy imaging of fly tissues revealed that, due to their characteristics, Notch complexes clustered with the transcription machinery and formed transcription hubs near their target genes. All cells exposed to Notch exhibited these hubs, but only a third produced the mRNAs associated with Notch target genes; adding a second signal (an insect hormone) significantly increased the proportion. This illustrates how 'chance' and collaboration influence the way the organism responds to Notch signalling. Finally, the experiments revealed that the hubs persisted for at least a day after removing the Notch signal. This 'molecular memory' led to cells responding faster when presented with Notch activity again. The work by DeHaro-Arbona sheds light on how individual cells respond to Notch signalling, and the factors that influence the activation of its target genes. This knowledge may prove useful when trying to better understand diseases in which this pathway is implicated, such as cancer.

A stochastic approach for co-evolution process of virus and human immune system.

Infectious diseases have long been a shaping force in human history, necessitating a comprehensive understanding of their dynamics. This study introduces a co-evolution model that integrates both epidemiological and evolutionary dynamics. Utilizing a system of differential equations, the model represents the interactions among susceptible, infected, and recovered populations for both ancestral and evolved viral strains. Methodologically rigorous, the model's existence and uniqueness have been verified, and it accommodates both deterministic and stochastic cases. A myriad of graphical techniques have been employed to elucidate the model's dynamics. Beyond its theoretical contributions, this model serves as a critical instrument for public health strategy, particularly predicting future outbreaks in scenarios where viral mutations compromise existing interventions.

A stochastic world model on gravity for stability inference.

The fact that objects without proper support will fall to the ground is not only a natural phenomenon, but also common sense in mind. Previous studies suggest that humans may infer objects' stability through a world model that performs mental simulations with a priori knowledge of gravity acting upon the objects. Here we measured participants' sensitivity to gravity to investigate how the world model works. We found that the world model on gravity was not a faithful replica of the physical laws, but instead encoded gravity's vertical direction as a Gaussian distribution. The world model with this stochastic feature fit nicely with participants' subjective sense of objects' stability and explained the illusion that taller objects are perceived as more likely to fall. Furthermore, a computational model with reinforcement learning revealed that the stochastic characteristic likely originated from experience-dependent comparisons between predictions formed by internal simulations and the realities observed in the external world, which illustrated the ecological advantage of stochastic representation in balancing accuracy and speed for efficient stability inference. The stochastic world model on gravity provides an example of how a priori knowledge of the physical world is implemented in mind that helps humans operate flexibly in open-ended environments.

DNA lesion bypass and the stochastic dynamics of transcription-coupled repair.

DNA base damage is a major source of oncogenic mutations and disruption to gene expression. The stalling of RNA polymerase II (RNAP) at sites of DNA damage and the subsequent triggering of repair processes have major roles in shaping the genome-wide distribution of mutations, clearing barriers to transcription, and minimizing the production of miscoded gene products. Despite its importance for genetic integrity, key mechanistic features of this transcription-coupled repair (TCR) process are controversial or unknown. Here, we exploited a well-powered in vivo mammalian model system to explore the mechanistic properties and parameters of TCR for alkylation damage at fine spatial resolution and with discrimination of the damaged DNA strand. For rigorous interpretation, a generalizable mathematical model of DNA damage and TCR was developed. Fitting experimental data to the model and simulation revealed that RNA polymerases frequently bypass lesions without triggering repair, indicating that small alkylation adducts are unlikely to be an efficient barrier to gene expression. Following a burst of damage, the efficiency of transcription-coupled repair gradually decays through gene bodies with implications for the occurrence and accurate inference of driver mutations in cancer. The reinitation of transcription from the repair site is not a general feature of transcription-coupled repair, and the observed data is consistent with reinitiation never taking place. Collectively, these results reveal how the directional but stochastic activity of TCR shapes the distribution of mutations following DNA damage.

On the validity of the state space correspondence strategy based on k-nearest neighbor cross-predictability in assessing directionality in stochastic systems: Application to cardiorespiratory coupling estimation.

We tested the validity of the state space correspondence (SSC) strategy based on k-nearest neighbor cross-predictability (KNNCP) to assess the directionality of coupling in stochastic nonlinear bivariate autoregressive (NBAR) processes. The approach was applied to assess closed-loop cardiorespiratory interactions between heart period (HP) variability and respiration (R) during a controlled respiration (CR) protocol in 19 healthy humans (aged from 27 to 35 yrs, 11 females) and during active standing (STAND) in 25 athletes (aged from 20 to 40 yrs, all men) and 25 non-athletes (aged from 20 to 40 yrs, all men). Over simulated NBAR processes, we found that (i) the SSC approach can detect the correct causal relationship as the direction leads to better KNNCP from the past of the driver to the future state of the target and (ii) simulations suggest that the ability of the method is preserved in any condition of complexity of the interacting series. Over CR and STAND protocols, we found that (a) slowing the breathing rate increases the strength of the causal relationship in both temporal directions in a balanced modality; (b) STAND is more powerful in modulating the coupling strength on the pathway from HP to R; (c) regardless of protocol and experimental condition, the strength of the link from HP to R is stronger than that from R to HP; (d) significant causal relationships in both temporal directions are found regardless of the level of complexity of HP variability and R. The SSC strategy is useful to disentangle closed-loop cardiorespiratory interactions.

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions.

Many imaging techniques for biological systems-like fixation of cells coupled with fluorescence microscopy-provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. These snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

Genetic, Environmental, and Stochastic Components of Lifespan Variability: The Drosophila Paradigm.

Lifespan is a complex quantitative trait involving genetic and non-genetic factors as well as the peculiarities of ontogenesis. As with all quantitative traits, lifespan shows considerable variation within populations and between individuals. Drosophila, a favourite object of geneticists, has greatly advanced our understanding of how different forms of variability affect lifespan. This review considers the role of heritable genetic variability, phenotypic plasticity and stochastic variability in controlling lifespan in Drosophila melanogaster. We discuss the major historical milestones in the development of the genetic approach to study lifespan, the breeding of long-lived lines, advances in lifespan QTL mapping, the environmental factors that have the greatest influence on lifespan in laboratory maintained flies, and the mechanisms, by which individual development affects longevity. The interplay between approaches to study ageing and lifespan limitation will also be discussed. Particular attention will be paid to the interaction of different types of variability in the control of lifespan.

Stochastic parabolic growth promotes coexistence and a relaxed error threshold in RNA-like replicator populations.

The RNA world hypothesis proposes that during the early evolution of life, primordial genomes of the first self-propagating evolutionary units existed in the form of RNA-like polymers. Autonomous, non-enzymatic, and sustained replication of such information carriers presents a problem, because product formation and hybridization between template and copy strands reduces replication speed. Kinetics of growth is then parabolic with the benefit of entailing competitive coexistence, thereby maintaining diversity. Here, we test the information-maintaining ability of parabolic growth in stochastic multispecies population models under the constraints of constant total population size and chemostat conditions. We find that large population sizes and small differences in the replication rates favor the stable coexistence of the vast majority of replicator species ('genes'), while the error threshold problem is alleviated relative to exponential amplification. In addition, sequence properties (GC content) and the strength of resource competition mediated by the rate of resource inflow determine the number of coexisting variants, suggesting that fluctuations in building block availability favored repeated cycles of exploration and exploitation. Stochastic parabolic growth could thus have played a pivotal role in preserving viable sequences generated by random abiotic synthesis and providing diverse genetic raw material to the early evolution of functional ribozymes.

All living things use molecules known as nucleic acids to store instructions on how to grow and maintain themselves and pass these instructions down to the next generation. However, it remains unclear how these systems may have evolved from simple molecules in the environment when life began over 3.6 billion years ago. One idea proposes that, before the first cells evolved, abiotic chemical processes gave rise to substantial building blocks of ribonucleic acids (or RNAs, for short). Over time, RNAs could have combined to form polymers of random sequences that started to copy themselves to make simple machines, only carrying the information required to make more of the same RNAs. Later on, these RNA molecules teamed up with proteins, fats and other molecules to make the first cells. When RNA replicates, the parent molecule is used as a template to assemble a new copy. While the new RNA molecule remains attached to its template it prevents the template being used to make more RNA. Therefore, it is thought that the speed at which a specific RNA machine copied itself may have varied in a pattern known as parabolic growth. Furthermore, when RNA replicates without the help of other biological molecules, the process is very prone to errors, which would have severely limited how much information the RNA machines were able to pass on to the next generation. Theoretical work suggested that under certain conditions, parabolic growth may favor the maintenance of a large amount of RNA sequence-coded information, but it is not clear if this is actually possible in nature. To address this question, Paczkó et al. developed mathematical models to investigate the effect of parabolic growth on the ability of RNA to replicate without other biological molecules. The models show that when large numbers of RNAs are present, small differences in how quickly different RNAs replicated favored the stable coexistence of different RNA sequences. Parabolic growth decreased the adverse effect of copying errors, allowing larger pieces of RNA to faithfully replicate themselves. This work suggests that parabolic growth may help to maintain different types of RNA (or similar replicating molecules) in a population and in turn, help new simple life forms to evolve. In the future, these findings may be used as a framework for laboratory experiments to better understand how early life forms may have evolved.

A Genuinely Hybrid, Multiscale 3D Cancer Invasion and Metastasis Modelling Framework.

We introduce in this paper substantial enhancements to a previously proposed hybrid multiscale cancer invasion modelling framework to better reflect the biological reality and dynamics of cancer. These model updates contribute to a more accurate representation of cancer dynamics, they provide deeper insights and enhance our predictive capabilities. Key updates include the integration of porous medium-like diffusion for the evolution of Epithelial-like Cancer Cells and other essential cellular constituents of the system, more realistic modelling of Epithelial-Mesenchymal Transition and Mesenchymal-Epithelial Transition models with the inclusion of Transforming Growth Factor beta within the tumour microenvironment, and the introduction of Compound Poisson Process in the Stochastic Differential Equations that describe the migration behaviour of the Mesenchymal-like Cancer Cells. Another innovative feature of the model is its extension into a multi-organ metastatic framework. This framework connects various organs through a circulatory network, enabling the study of how cancer cells spread to secondary sites.

Modeling single cell trajectory using forward-backward stochastic differential equations.

Recent advances in single-cell sequencing technology have provided opportunities for mathematical modeling of dynamic developmental processes at the single-cell level, such as inferring developmental trajectories. Optimal transport has emerged as a promising theoretical framework for this task by computing pairings between cells from different time points. However, optimal transport methods have limitations in capturing nonlinear trajectories, as they are static and can only infer linear paths between endpoints. In contrast, stochastic differential equations (SDEs) offer a dynamic and flexible approach that can model non-linear trajectories, including the shape of the path. Nevertheless, existing SDE methods often rely on numerical approximations that can lead to inaccurate inferences, deviating from true trajectories. To address this challenge, we propose a novel approach combining forward-backward stochastic differential equations (FBSDE) with a refined approximation procedure. Our FBSDE model integrates the forward and backward movements of two SDEs in time, aiming to capture the underlying dynamics of single-cell developmental trajectories. Through comprehensive benchmarking on multiple scRNA-seq datasets, we demonstrate the superior performance of FBSDE compared to other methods, highlighting its efficacy in accurately inferring developmental trajectories.

Stochastic portfolio optimization: A regret-based approach on volatility risk measures: An empirical evidence from The New York stock market.

Portfolio optimization involves finding the ideal combination of securities and shares to reduce risk and increase profit in an investment. To assess the impact of risk in portfolio optimization, we utilize a significant volatility risk measure series. Behavioral finance biases play a critical role in portfolio optimization and the efficient allocation of stocks. Regret, within the realm of behavioral finance, is the feeling of remorse that causes hesitation in making significant decisions and avoiding actions that could lead to poor investment choices. This behavior often leads investors to hold onto losing investments for extended periods, refusing to acknowledge mistakes and accept losses. Ironically, by evading regret, investors may miss out on potential opportunities. in this paper, our purpose is to compare investment scenarios in the decision-making process and calculate the amount of regret obtained in each scenario. To accomplish this, we consider volatility risk metrics and utilize stochastic optimization to identify the most suitable scenario that not only maximizes yield in the investment portfolio and minimizes risk, but also minimizes resulting regret. To convert each multi-objective model into a single objective, we employ the augmented epsilon constraint (AEC) method to establish the Pareto efficiency frontier. As a means of validating the solution of this method, we analyze data spanning 20, 50, and 100 weeks from 150 selected stocks in the New York market based on fundamental analysis. The results show that the selection of the mad risk measure in the time horizon of 100 weeks with a regret rate of 0.104 is the most appropriate research scenario. this article recommended that investors diversify their portfolios by investing in a variety of assets. This can help reduce risk and increase overall returns and improve financial literacy among investors.

Quantitative analysis and stochastic modeling of osteophyte formation and growth process on human vertebrae based on radiographs: a follow-up study.

Osteophytes are frequently observed in elderly people and most commonly appear at the anterior edge of the cervical and lumbar vertebrae body. The anterior osteophytes keep developing and will lead to neck/back pain over time. In clinical practice, the accurate measurement of the anterior osteophyte length and the understanding of the temporal progression of anterior osteophyte growth are of vital importance to clinicians for effective treatment planning. This study proposes a new measuring method using the osteophyte ratio index to quantify anterior osteophyte length based on lateral radiographs. Moreover, we develop a continuous stochastic degradation model with time-related functions to characterize the anterior osteophyte formation and growth process on cervical and lumbar vertebrae over time. Follow-up data of anterior osteophytes up to 9 years are obtained for measurement and model validation. The agreement test indicates excellent reproducibility for our measuring method. The proposed model accurately fits the osteophyte growth paths. The model predicts the mean time to onset of pain and obtained survival function of the degenerative vertebrae. This research opens the door to future quantification and mathematical modeling of the anterior osteophyte growth on human cervical and lumbar vertebrae. The measured follow-up data is shared for future studies.

Effects of mRNA conformational switching on translational noise in gene circuits.

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model.

The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: µECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets.

Gossip-based distributed stochastic mirror descent for constrained optimization.

This paper considers a distributed constrained optimization problem over a multi-agent network in the non-Euclidean sense. The gossip protocol is adopted to relieve the communication burden, which also adapts to the constantly changing topology of the network. Based on this idea, a gossip-based distributed stochastic mirror descent (GB-DSMD) algorithm is proposed to handle the problem under consideration. The performances of GB-DSMD algorithms with constant and diminishing step sizes are analyzed, respectively. When the step size is constant, the error bound between the optimal function value and the expected function value corresponding to the average iteration output of the algorithm is derived. While for the case of the diminishing step size, it is proved that the output of the algorithm uniformly approaches to the optimal value with probability 1. Finally, as a numerical example, the distributed logistic regression is reported to demonstrate the effectiveness of the GB-DSMD algorithm.

Solving stochastic gene-expression models using queueing theory: A tutorial review.

Stochastic models of gene expression are typically formulated using the chemical master equation, which can be solved exactly or approximately using a repertoire of analytical methods. Here, we provide a tutorial review of an alternative approach based on queueing theory that has rarely been used in the literature of gene expression. We discuss the interpretation of six types of infinite-server queues from the angle of stochastic single-cell biology and provide analytical expressions for the stationary and nonstationary distributions and/or moments of mRNA/protein numbers and bounds on the Fano factor. This approach may enable the solution of complex models that have hitherto evaded analytical solution.

Bacterial and fungal community structures in Hulun Lake are regulated by both stochastic processes and environmental factors.

Microorganisms are a crucial component of lake ecosystems and significant contributors to biogeochemical cycles. However, the understanding of how primary microorganism groups (e.g., bacteria and fungi) are distributed and constructed within different lake habitats is lacking. We investigated the bacterial and fungal communities of Hulun Lake using high-throughput sequencing techniques targeting 16S rRNA and Internal Transcribed Spacer 2 genes, including a range of ecological and statistical methodologies. Our findings reveal that environmental factors have high spatial and temporal variability. The composition and community structures vary significantly depending on differences in habitats. Variance partitioning analysis showed that environmental and geographical factors accounted for <20% of the community variation. Canonical correlation analysis showed that among the environmental factors, temperature, pH, and dissolved oxygen had strong control over microbial communities. However, the microbial communities (bacterial and fungal) were primarily controlled by the dispersal limitations of stochastic processes. This study offers fresh perspectives regarding the maintenance mechanism of bacterial and fungal biodiversity in lake ecosystems, especially regarding the responses of microbial communities under identical environmental stress.IMPORTANCELake ecosystems are an important part of the freshwater ecosystem. Lake microorganisms play an important role in material circulation and energy flow owing to their unique enzymatic and metabolic capacity. In this study, we observed that bacterial and fungal communities varied widely in the water and sediments of Hulun Lake. The primary factor affecting their formation was identified as dispersal limitation during stochastic processes. Environmental and geographical factors accounted for <20% of the variation in bacterial and fungal communities, with pH, temperature, and dissolved oxygen being important environmental factors. Our findings provide new insights into the responses of bacteria and fungi to the environment, shed light on the ecological processes of community building, and deepen our understanding of lake ecosystems. The results of this study provide a reference for lake management and conservation, particularly with respect to monitoring and understanding microbial communities in response to environmental changes.

Stability against fluctuations: a two-dimensional study of scaling, bifurcations and spontaneous symmetry breaking in stochastic models of synaptic plasticity.

Stochastic models of synaptic plasticity must confront the corrosive influence of fluctuations in synaptic strength on patterns of synaptic connectivity. To solve this problem, we have proposed that synapses act as filters, integrating plasticity induction signals and expressing changes in synaptic strength only upon reaching filter threshold. Our earlier analytical study calculated the lifetimes of quasi-stable patterns of synaptic connectivity with synaptic filtering. We showed that the plasticity step size in a stochastic model of spike-timing-dependent plasticity (STDP) acts as a temperature-like parameter, exhibiting a critical value below which neuronal structure formation occurs. The filter threshold scales this temperature-like parameter downwards, cooling the dynamics and enhancing stability. A key step in this calculation was a resetting approximation, essentially reducing the dynamics to one-dimensional processes. Here, we revisit our earlier study to examine this resetting approximation, with the aim of understanding in detail why it works so well by comparing it, and a simpler approximation, to the system's full dynamics consisting of various embedded two-dimensional processes without resetting. Comparing the full system to the simpler approximation, to our original resetting approximation, and to a one-afferent system, we show that their equilibrium distributions of synaptic strengths and critical plasticity step sizes are all qualitatively similar, and increasingly quantitatively similar as the filter threshold increases. This increasing similarity is due to the decorrelation in changes in synaptic strength between different afferents caused by our STDP model, and the amplification of this decorrelation with larger synaptic filters.

Neutral competition explains the clonal composition of neural organoids.

Neural organoids model the development of the human brain and are an indispensable tool for studying neurodevelopment. Whole-organoid lineage tracing has revealed the number of progenies arising from each initial stem cell to be highly diverse, with lineage sizes ranging from one to more than 20,000 cells. This high variability exceeds what can be explained by existing stochastic models of corticogenesis and indicates the existence of an additional source of stochasticity. To explain this variability, we introduce the SAN model which distinguishes Symmetrically diving, Asymmetrically dividing, and Non-proliferating cells. In the SAN model, the additional source of stochasticity is the survival time of a lineage's pool of symmetrically dividing cells. These survival times result from neutral competition within the sub-population of all symmetrically dividing cells. We demonstrate that our model explains the experimentally observed variability of lineage sizes and derive the quantitative relationship between survival time and lineage size. We also show that our model implies the existence of a regulatory mechanism which keeps the size of the symmetrically dividing cell population constant. Our results provide quantitative insight into the clonal composition of neural organoids and how it arises. This is relevant for many applications of neural organoids, and similar processes may occur in other developing tissues both in vitro and in vivo.

Scaling of stochastic growth and division dynamics: A comparative study of individual rod-shaped cells in the Mother Machine and SChemostat platforms.

Microfluidic platforms enable long-term quantification of stochastic behaviors of individual bacterial cells under precisely controlled growth conditions. Yet, quantitative comparisons of physiological parameters and cell behaviors of different microorganisms in different experimental and device modalities is not available due to experiment-specific details affecting cell physiology. To rigorously assess the effects of mechanical confinement, we designed, engineered, and performed side-by-side experiments under otherwise identical conditions in the Mother Machine (with confinement) and the SChemostat (without confinement), using the latter as the ideal comparator. We established a protocol to cultivate a suitably engineered rod-shaped mutant of Caulobacter crescentus in the Mother Machine and benchmarked the differences in stochastic growth and division dynamics with respect to the SChemostat. While the single-cell growth rate distributions are remarkably similar, the mechanically confined cells in the Mother Machine experience a substantial increase in interdivision times. However, we find that the division ratio distribution precisely compensates for this increase, which in turn reflects identical emergent simplicities governing stochastic intergenerational homeostasis of cell sizes across device and experimental configurations, provided the cell sizes are appropriately mean-rescaled in each condition. Our results provide insights into the nature of the robustness of the bacterial growth and division machinery.