Heterogeneity in Pure Microbial Systems: Experimental Measurements and Modeling

Synthetic microbial ecology: engineering habitats for modular consortia

Microbiomes, the complex networks of micro-organisms and the molecules through which they interact, play a crucial role in health and ecology. Over at least the past two decades, engineering biology has made significant progress, impacting the bio-based industry, health, and environmental sectors; but has only recently begun to explore the engineering of microbial ecosystems. The creation of synthetic microbial communities presents opportunities to help us understand the dynamics of wild ecosystems, learn how to manipulate and interact with existing microbiomes for therapeutic and other purposes, and to create entirely new microbial communities capable of undertaking tasks for industrial biology. Here, we describe how synthetic ecosystems can be constructed and controlled, focusing on how the available methods and interaction mechanisms facilitate the regulation of community composition and output. While experimental decisions are dictated by intended applications, the vast number of tools available suggests great opportunity for researchers to develop a diverse array of novel microbial ecosystems.

bglG Regulates the Heterogeneity Driven by the Acid Tolerance Response in Lacticaseibacillus paracasei L9

The acid tolerance of lactic acid bacteria is crucial for their fermentation and probiotic functions. Acid adaption significantly enhances the acid tolerance of strains, and the phenotypic heterogeneity driven by the acid tolerance response (ATR) contributes to this process by providing a selective advantage in harsh environments. The mechanism of heterogeneity under the ATR is not yet clear, but individual gene expression differences are recognized as the cause. In this study, we observed four heterogeneous subpopulations (viable, injured, dead, and unstained) of Lacticaseibacillus paracasei L9 (L9) induced by acid adaption (pH 5.0, 40 min) using flow cytometry. The viable subpopulation represented a significantly superior acid tolerance to the injured subpopulation or total population. Different subpopulations were sorted and transcriptomic analysis was performed. Five genes were found to be upregulated in the viable subpopulation and downregulated in the injured subpopulation, and bglG (LPL9_RS14735) was identified as having a key role in this process. Using salicin (glucoside)-inducing gene expression and gene insertion mutagenesis, we verified that bglG regulated the heterogeneity of the acid stress response and that the relevant mechanisms might be related to activating hsp20. This study provides new evidence for the mechanism of the ATR and may contribute to the theoretical basis of improving the acid tolerance of Lacticaseibacillus paracasei L9.

Metabolic heterogeneity and cross-feeding within isogenic yeast populations captured by DILAC

Genetically identical cells are known to differ in many physiological parameters such as growth rate and drug tolerance. Metabolic specialization is believed to be a cause of such phenotypic heterogeneity, but detection of metabolically divergent subpopulations remains technically challenging. We developed a proteomics-based technology, termed differential isotope labelling by amino acids (DILAC), that can detect producer and consumer subpopulations of a particular amino acid within an isogenic cell population by monitoring peptides with multiple occurrences of the amino acid. We reveal that young, morphologically undifferentiated yeast colonies contain subpopulations of lysine producers and consumers that emerge due to nutrient gradients. Deconvoluting their proteomes using DILAC, we find evidence for in situ cross-feeding where rapidly growing cells ferment and provide the more slowly growing, respiring cells with ethanol. Finally, by combining DILAC with fluorescence-activated cell sorting, we show that the metabolic subpopulations diverge phenotypically, as exemplified by a different tolerance to the antifungal drug amphotericin B. Overall, DILAC captures previously unnoticed metabolic heterogeneity and provides experimental evidence for the role of metabolic specialization and cross-feeding interactions as a source of phenotypic heterogeneity in isogenic cell populations.

Diverse mechanisms of bioproduction heterogeneity in fermentation and their control strategies

Microbial bioproduction often faces challenges related to populational heterogeneity, where cells exhibit varying biosynthesis capabilities. Bioproduction heterogeneity can stem from genetic and non-genetic factors, resulting in decreased titer, yield, stability, and reproducibility. Consequently, understanding and controlling bioproduction heterogeneity are crucial for enhancing the economic competitiveness of large-scale biomanufacturing. In this review, we provide a comprehensive overview of current understandings of the various mechanisms underlying bioproduction heterogeneity. Additionally, we examine common strategies for controlling bioproduction heterogeneity based on these mechanisms. By implementing more robust measures to mitigate heterogeneity, we anticipate substantial enhancements in the scalability and stability of bioproduction processes.

ONE-SENTENCE SUMMARY

This review summarizes current understandings of different mechanisms of bioproduction heterogeneity and common control strategies based on these mechanisms.

Advances in automated real-time flow cytometry for monitoring of bioreactor processes

Flow cytometry and its technological possibilities have greatly advanced in the past decade as analysis tool for single cell properties and population distributions of different cell types in bioreactors. Along the way, some solutions for automated real-time flow cytometry (ART-FCM) were developed for monitoring of bioreactor processes without operator interference over extended periods with variable sampling frequency. However, there is still great potential for ART-FCM to evolve and possibly become a standard application in bioprocess monitoring and process control. This review first addresses different components of an ART-FCM, including the sampling device, the sample-processing unit, the unit for sample delivery to the flow cytometer and the settings for measurement of pre-processed samples. Also, available algorithms are presented for automated data analysis of multi-parameter fluorescence datasets derived from ART-FCM experiments. Furthermore, challenges are discussed for integration of fluorescence-activated cell sorting into an ART-FCM setup for isolation and separation of interesting subpopulations that can be further characterized by for instance omics-methods. As the application of ART-FCM is especially of interest for bioreactor process monitoring, including investigation of population heterogeneity and automated process control, a summary of already existing setups for these purposes is given. Additionally, the general future potential of ART-FCM is addressed.

Study of plasmid-based expression level heterogeneity under plasmid-curing like conditions in Cupriavidus necator

Plasmid expression level heterogeneity in Cupriavidus necator was studied in response to stringent culture conditions, supposed to enhance plasmid instability, through plasmid curing strategies. Two plasmid curing strategies were compared based on their efficiency at generating heterogeneity in batch: rifampicin addition and temperature increase. A temperature increase from 30° to 37 °C was the most efficient plasmid curing strategy. To generate a heterogeneous population in terms of plasmid expression levels, successive batches at supra-optimal culture temperature (i.e. 37 °C) were initially conducted. Three distinct fluorescent subpopulations P₀ (not fluorescent), P₁ (low fluorescence intensity, median = 1 10³) and P₂ (high fluorescence intensity, median = 6 10³) were obtained. From there, the chemostat culture was implemented to study the long-term stress response under well-controlled environment at defined dilution rates. For dilution rates comprised between 0.05 and 0.10 h^-1, the subpopulation P₂ (62% vs 90%) was favored compared to P₁ cells (54% vs 1%), especially when growth rate increased. Our biosensor was efficient at discriminating subpopulation presenting different expression levels under stringent culture conditions. Plus, we showed that controlling growth kinetics had a stabilizing impact on plasmid expression levels, even under heterogeneous expression conditions.

Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol^-1, and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.

The Inclusion of the Food Microstructural Influence in Predictive Microbiology: State-of-the-Art

Predictive microbiology has steadily evolved into one of the most important tools to assess and control the microbiological safety of food products. Predictive models were traditionally developed based on experiments in liquid laboratory media, meaning that food microstructural effects were not represented in these models. Since food microstructure is known to exert a significant effect on microbial growth and inactivation dynamics, the applicability of predictive models is limited if food microstructure is not taken into account. Over the last 10-20 years, researchers, therefore, developed a variety of models that do include certain food microstructural influences. This review provides an overview of the most notable microstructure-including models which were developed over the years, both for microbial growth and inactivation.

Recombinant Protein Production with Escherichia coli in Glucose and Glycerol Limited Chemostats

Recently, there has been a resurgence of interest in continuous bioprocessing as a cost-optimised production strategy, driven by a rising global requirement for recombinant proteins used as biological drugs. This strategy could provide several benefits over traditional batch processing, including smaller bioreactors, smaller facilities, and overall reduced plant footprints and investment costs. Continuous processes may also offer improved product quality and minimise heterogeneity, both in the culture and in the product. In this paper, a model protein, green fluorescent protein (GFP) mut3*, was used to test the recombinant protein expression in an Escherichia coli strain with industrial relevance grown in chemostat. An important factor in enabling stable productivity in continuous cultures is the carbon source. We have studied the viability and heterogeneity of the chemostat cultures using a chemically defined medium based on glucose or glycerol as the single carbon source. As a by-product of biodiesel production, glycerol is expected to become a sustainable alternative substrate to glucose. We have found that although glycerol gives a higher cell density, it also generates higher heterogeneity in the culture and a less stable recombinant protein production. We suggest that manipulating the balance between different subpopulations to increase the proportion of productive cells may be a possible solution for making glycerol a successful alternative to glucose.

Computational modeling of dynein motor proteins at work

Along with various experimental methods, a combination of theoretical and computational methods is essential to explore different length-scale and time-scale processes in the biological system. The functional mechanism of a dynein, an ATP-fueled motor protein, working in a multiprotein complex, involves a wide range of length/time-scale events. It generates mechanical force from chemical energy and moves on microtubules towards the minus end direction while performing a large number of biological processes including ciliary beating, intracellular material transport, and cell division. Like in the cases of other conventional motor proteins, a combination of experimental techniques including X-crystallography, cryo-electron microscopy, and single molecular assay have provided a wealth of information about the mechanochemical cycle of a dynein. Dyneins have a large and complex structural architecture and therefore, computational modeling of different aspects of a dynein is extremely challenging. As the process of dynein movement involves varying length and timescales, it demands, like in experiments, a combination of computational methods covering such a wide range of processes for the comprehensive investigation of the mechanochemical cycle. In this review article, we will summarize how the use of state-of-the-art computational methods can provide a detailed molecular understanding of the mechanochemical cycle of the dynein. We implemented all-atom molecular dynamics simulations and hybrid quantum-mechanics/molecular-mechanics simulations to explore the ATP hydrolysis mechanisms at the primary ATPase site (AAA1) of dynein. To investigate the large-scale conformational changes we employed coarse-grained structure-based molecular dynamics simulations to capture the domain motions. Here we explored the conformational changes upon binding of ATP at AAA1, nucleotide state-dependent regulation of the mechanochemical cycle, and inter-head coordination by inter-head tension. Additionally, implementing a phenomenological theoretical model we explore the force-dependent detachment rate of a motorhead from the microtubule and the principle of multi-dynein cooperation during cargo transport.

Engineering Pseudomonas putida for efficient aromatic conversion to bioproduct using high throughput screening in a bioreactor

Pseudomonas putida KT2440 is an emerging biomanufacturing host amenable for use with renewable carbon streams including aromatics such as para-coumarate. We used a pooled transposon library disrupting nearly all (4,778) non-essential genes to characterize this microbe under common stirred-tank bioreactor parameters with quantitative fitness assays. Assessing differential fitness values by monitoring changes in mutant strain abundance identified 33 gene mutants with improved fitness across multiple stirred-tank bioreactor formats. Twenty-one deletion strains from this subset were reconstructed, including GacA, a regulator, TtgB, an ABC transporter, and PP_0063, a lipid A acyltransferase. Thirteen deletion strains with roles in varying cellular functions were evaluated for conversion of para-coumarate, to a heterologous bioproduct, indigoidine. Several mutants, such as the ΔgacA strain improved fitness in a bioreactor by 35 fold and showed an 8-fold improvement in indigoidine production (4.5 g/L, 0.29 g/g, 23% of maximum theoretical yield) from para-coumarate as the carbon source.

Microbial single-cell omics: the crux of the matter

Abstract

Single-cell genomics and transcriptomics can provide reliable context for assembled genome fragments and gene expression activity on the level of individual prokaryotic genomes. These methods are rapidly emerging as an essential complement to cultivation-based, metagenomics, metatranscriptomics, and microbial community-focused research approaches by allowing direct access to information from individual microorganisms, even from deep-branching phylogenetic groups that currently lack cultured representatives. Their integration and binning with environmental ‘omics data already provides unprecedented insights into microbial diversity and metabolic potential, enabling us to provide information on individual organisms and the structure and dynamics of natural microbial populations in complex environments. This review highlights the pitfalls and recent advances in the field of single-cell omics and its importance in microbiological and biotechnological studies.

Key points

• Single-cell omics expands the tree of life through the discovery of novel organisms, genes, and metabolic pathways.

• Disadvantages of metagenome-assembled genomes are overcome by single-cell omics.

• Functional analysis of single cells explores the heterogeneity of gene expression.

• Technical challenges still limit this field, thus prompting new method developments.

Toward Engineering Biosystems With Emergent Collective Functions

Many complex behaviors in biological systems emerge from large populations of interacting molecules or cells, generating functions that go beyond the capabilities of the individual parts. Such collective phenomena are of great interest to bioengineers due to their robustness and scalability. However, engineering emergent collective functions is difficult because they arise as a consequence of complex multi-level feedback, which often spans many length-scales. Here, we present a perspective on how some of these challenges could be overcome by using multi-agent modeling as a design framework within synthetic biology. Using case studies covering the construction of synthetic ecologies to biological computation and synthetic cellularity, we show how multi-agent modeling can capture the core features of complex multi-scale systems and provide novel insights into the underlying mechanisms which guide emergent functionalities across scales. The ability to unravel design rules underpinning these behaviors offers a means to take synthetic biology beyond single molecules or cells and toward the creation of systems with functions that can only emerge from collectives at multiple scales.

Plasmid expression level heterogeneity monitoring via heterologous eGFP production at the single-cell level in Cupriavidus necator

A methodology for plasmid expression level monitoring of eGFP expression suitable for dynamic processes was assessed during fermentation. This technique was based on the expression of a fluorescent biosensor (eGFP) encoded on a recombinant plasmid coupled to single-cell analysis. Fluorescence intensity at single-cell level was measured by flow cytometry. We demonstrated that promoter evaluation based on single-cell analysis versus classic global analysis brings valuable insights. Single-cell analysis pointed out the fact that intrinsic fluorescence increased with the strength of the promoter up to a threshold. Beyond that, cell permeability increases to excrete the fluorescent protein in the medium. The metabolic load due to the increase in the eGFP production in the case of strong constitutive promoters leads to slower growth kinetics compared with plasmid-free cells. With the strain Cupriavidus necator Re2133, growth rate losses were measured from 3% with the weak constitutive promoter P_lac to 56% with the strong constitutive promoter P_j5. Through this work, it seems crucial to find a compromise between the fluorescence intensity in single cells and the metabolic load; in our conditions, the best compromise found was the weak promoter P_lac. The plasmid expression level monitoring method was tested in the presence of a heterogeneous population induced by plasmid-curing methods. For all the identified subpopulations, the plasmid expression level heterogeneity was significantly detected at the level of fluorescence intensity in single cells. After cell sorting, growth rate and cultivability were assessed for each subpopulation. In conclusion, this eGFP biosensor makes it possible to follow the variations in the level of plasmid expression under conditions of population heterogeneity.Key Points•Development of a plasmid expression level monitoring method at the single-cell level by flow cytometry.•Promoter evaluation by single-cell analysis: cell heterogeneity and strain robustness.•Reporter system optimization for efficient subpopulation detection in pure cultures.

INDISIM-Denitrification, an individual-based model for study the denitrification process

Denitrification is one of the key processes of the global nitrogen (N) cycle driven by bacteria. It has been widely known for more than 100 years as a process by which the biogeochemical N-cycle is balanced. To study this process, we develop an individual-based model called INDISIM-Denitrification. The model embeds a thermodynamic model for bacterial yield prediction inside the individual-based model INDISIM and is designed to simulate in aerobic and anaerobic conditions the cell growth kinetics of denitrifying bacteria. INDISIM-Denitrification simulates a bioreactor that contains a culture medium with succinate as a carbon source, ammonium as nitrogen source and various electron acceptors. To implement INDISIM-Denitrification, the individual-based model INDISIM was used to give sub-models for nutrient uptake, stirring and reproduction cycle. Using a thermodynamic approach, the denitrification pathway, cellular maintenance and individual mass degradation were modeled using microbial metabolic reactions. These equations are the basis of the sub-models for metabolic maintenance, individual mass synthesis and reducing internal cytotoxic products. The model was implemented in the open-access platform NetLogo. INDISIM-Denitrification is validated using a set of experimental data of two denitrifying bacteria in two different experimental conditions. This provides an interactive tool to study the denitrification process carried out by any denitrifying bacterium since INDISIM-Denitrification allows changes in the microbial empirical formula and in the energy-transfer-efficiency used to represent the metabolic pathways involved in the denitrification process. The simulator can be obtained from the authors on request.

Quantitative Flow Cytometry to Understand Population Heterogeneity in Response to Changes in Substrate Availability in Escherichia coli and Saccharomyces cerevisiae Chemostats

Microbial cells in bioprocesses are usually described with averaged parameters. But in fact, single cells within populations vary greatly in characteristics such as stress resistance, especially in response to carbon source gradients. Our aim was to introduce tools to quantify population heterogeneity in bioprocesses using a combination of reporter strains, flow cytometry, and easily comprehensible parameters. We calculated mean, mode, peak width, and coefficient of variance to describe distribution characteristics and temporal shifts in fluorescence intensity. The skewness and the slope of cumulative distribution function plots illustrated differences in distribution shape. These parameters are person-independent and precise. We demonstrated this by quantifying growth-related population heterogeneity of Saccharomyces cerevisiae and Escherichia coli reporter strains in steady-state of aerobic glucose-limited chemostat cultures at different dilution rates and in response to glucose pulses. Generally, slow-growing cells showed stronger responses to glucose excess than fast-growing cells. Cell robustness, measured as membrane integrity after exposure to freeze-thaw treatment, of fast-growing cells was strongly affected in subpopulations of low membrane robustness. Glucose pulses protected subpopulations of fast-growing but not slower-growing yeast cells against membrane damage. Our parameters could successfully describe population heterogeneity, thereby revealing physiological characteristics that might have been overlooked during traditional averaged analysis.

Coculturing Bacteria Leads to Reduced Phenotypic Heterogeneities

Isogenic bacterial populations are known to exhibit phenotypic heterogeneity at the single-cell level. Because of difficulties in assessing the phenotypic heterogeneity of a single taxon in a mixed community, the importance of this deeper level of organization remains relatively unknown for natural communities. In this study, we have used membrane-based microcosms that allow the probing of the phenotypic heterogeneity of a single taxon while interacting with a synthetic or natural community. Individual taxa were studied under axenic conditions, as members of a coculture with physical separation, and as a mixed culture. Phenotypic heterogeneity was assessed through both flow cytometry and Raman spectroscopy. Using this setup, we investigated the effect of microbial interactions on the individual phenotypic heterogeneities of two interacting drinking water isolates. Through flow cytometry we have demonstrated that interactions between these bacteria lead to a reduction of their individual phenotypic diversities and that this adjustment is conditional on the bacterial taxon. Single-cell Raman spectroscopy confirmed a taxon-dependent phenotypic shift due to the interaction. In conclusion, our data suggest that bacterial interactions may be a general driver of phenotypic heterogeneity in mixed microbial populations.IMPORTANCE Laboratory studies have shown the impact of phenotypic heterogeneity on the survival and functionality of isogenic populations. Because phenotypic heterogeneity plays an important role in pathogenicity and virulence, antibiotic resistance, biotechnological applications, and ecosystem properties, it is crucial to understand its influencing factors. An unanswered question is whether bacteria in mixed communities influence the phenotypic heterogeneity of their community partners. We found that coculturing bacteria leads to a reduction in their individual phenotypic heterogeneities, which led us to the hypothesis that the individual phenotypic diversity of a taxon is dependent on the community composition.

Maximum entropy and population heterogeneity in continuous cell cultures

Continuous cultures of mammalian cells are complex systems displaying hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis, as well as sharp transitions between different metabolic states. In this context mathematical models may suggest control strategies to steer the system towards desired states. Although even clonal populations are known to exhibit cell-to-cell variability, most of the currently studied models assume that the population is homogeneous. To overcome this limitation, we use the maximum entropy principle to model the phenotypic distribution of cells in a chemostat as a function of the dilution rate. We consider the coupling between cell metabolism and extracellular variables describing the state of the bioreactor and take into account the impact of toxic byproduct accumulation on cell viability. We present a formal solution for the stationary state of the chemostat and show how to apply it in two examples. First, a simplified model of cell metabolism where the exact solution is tractable, and then a genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line. Along the way we discuss several consequences of heterogeneity, such as: qualitative changes in the dynamical landscape of the system, increasing concentrations of byproducts that vanish in the homogeneous case, and larger population sizes.

Defining Cell Identity with Single-Cell Omics

Cells are a fundamental unit of life, and the ability to study the phenotypes and behaviors of individual cells is crucial to understanding the workings of complex biological systems. Cell phenotypes (epigenomic, transcriptomic, proteomic, and metabolomic) exhibit dramatic heterogeneity between and within the different cell types and states underlying cellular functional diversity. Cell genotypes can also display heterogeneity throughout an organism, in the form of somatic genetic variation-most notably in the emergence and evolution of tumors. Recent technical advances in single-cell isolation and the development of omics approaches sensitive enough to reveal these aspects of cell identity have enabled a revolution in the study of multicellular systems. In this review, we discuss the technologies available to resolve the genomes, epigenomes, transcriptomes, proteomes, and metabolomes of single cells from a wide variety of living systems.

Population heterogeneity in microbial bioprocesses: origin, analysis, mechanisms, and future perspectives

Population heterogeneity is omnipresent in all bioprocesses even in homogenous environments. Its origin, however, is only so well understood that potential strategies like bet-hedging, noise in gene expression and division of labour that lead to population heterogeneity can be derived from experimental studies simulating the dynamics in industrial scale bioprocesses. This review aims at summarizing the current state of the different parts of single cell studies in bioprocesses. This includes setups to visualize different phenotypes of single cells, computational approaches connecting single cell physiology with environmental influence and special cultivation setups like scale-down reactors that have been proven to be useful to simulate large-scale conditions. A step in between investigation of populations and single cells is studying subpopulations with distinct properties that differ from the rest of the population with sub-omics methods which are also presented here. Moreover, the current knowledge about population heterogeneity in bioprocesses is summarized for relevant industrial production hosts and mixed cultures, as they provide the unique opportunity to distribute metabolic burden and optimize production processes in a way that is impossible in traditional monocultures. In the end, approaches to explain the underlying mechanism of population heterogeneity and the evidences found to support each hypothesis are presented. For instance, population heterogeneity serving as a bet-hedging strategy that is used as coordinated action against bioprocess-related stresses while at the same time spreading the risk between individual cells as it ensures the survival of least a part of the population in any environment the cells encounter.

Digital Image Analysis of Yeast Single Cells Growing in Two Different Oxygen Concentrations to Analyze the Population Growth and to Assist Individual-Based Modeling

Nowadays control of the growth of Saccharomyces to obtain biomass or cellular wall components is crucial for specific industrial applications. The general aim of this contribution is to deal with experimental data obtained from yeast cells and from yeast cultures to attempt the integration of the two levels of information, individual and population, to progress in the control of yeast biotechnological processes by means of the overall analysis of this set of experimental data, and to assist in the improvement of an individual-based model, namely, INDISIM-Saccha. Populations of S. cerevisiae growing in liquid batch culture, in aerobic and microaerophilic conditions, were studied. A set of digital images was taken during the population growth, and a protocol for the treatment and analyses of the images obtained was established. The piecewise linear model of Buchanan was adjusted to the temporal evolutions of the yeast populations to determine the kinetic parameters and changes of growth phases. In parallel, for all the yeast cells analyzed, values of direct morphological parameters, such as area, perimeter, major diameter, minor diameter, and derived ones, such as circularity and elongation, were obtained. Graphical and numerical methods from descriptive statistics were applied to these data to characterize the growth phases and the budding state of the yeast cells in both experimental conditions, and inferential statistical methods were used to compare the diverse groups of data achieved. Oxidative metabolism of yeast in a medium with oxygen available and low initial sugar concentration can be taken into account in order to obtain a greater number of cells or larger cells. Morphological parameters were analyzed statistically to identify which were the most useful for the discrimination of the different states, according to budding and/or growth phase, in aerobic and microaerophilic conditions. The use of the experimental data for subsequent modeling work was then discussed and compared to simulation results generated with INDISIM-Saccha, which allowed us to advance in the development of this yeast model, and illustrated the utility of data at different levels of observation and the needs and logic behind the development of a microbial individual-based model.