Metagenomics has lowered the barrier to microbial discovery--enabling the identification of novel microbes without isolation--but cultures remain imperative for the deep study of microbes. Cultivation and isolation of non-model microbes remains a major challenge, despite advances in high-throughput culturomic methods. The quantity of simultaneous experimental variables is constrained by time and resources, but the list can be reduced using computational biology. Given an annotated genome, metabolic modelling can be used to predict source nutrients required for the growth of a microbe, which acts as an initial screen to inform culture and isolation experiments. This chapter provides an overview of metabolic networks and modelling and how they can be used to predict the nutrient requirements of a microorganism, followed by a sample protocol using a toy metabolic network, which is then expanded to a genome-scale metabolic network application. These methods can be applied to any metabolic network of interest--which in turn can be created from any genome of interest--and are a starting point for experimental validation of source nutrients required for microorganisms that remain uncultiv
A model for morphoelastic growth, that is, growth influenced by elastic stress, driven by the absorption of nutrients is considered. The model features a multiplicative decomposition of the deformation gradient into an elastic contribution and a growth tensor. While the evolution of the system is governed by an ordinary differential equation for the growth tensor on a suitable Banach space, which depends on the elastic stresses and the concentration of a nutrient field, the total deformation is given by the solution of a quasi-static equilibrium equation arising from the formal Euler-Lagrange equations of a hyperelastic variational integral. The nutrient concentration is determined by a linear elliptic reaction-diffusion equation which is formulated in Lagrangian coordinates and whose coefficients depend on the growth tensor as well as the deformation gradient accounting for the change of material properties due to elastic deformation and growth. Existence and uniqueness of solutions of this fully coupled system of differential equations is proven via a fixed-point argument.
Regularity properties of solutions for a class of quasi-stationary models in one spatial dimension for stress-modulated growth in the presence of a nutrient field are proven. At a given point in time the configuration of a body after pure growth is determined by means of a family of ordinary differential equations in every point in space. Subsequently, an elastic deformation, which is given by the minimizer of a hyperelastic variational integral, is applied in order to restore Dirichlet boundary conditions. While the ordinary differential equations governing the growth process depend on the elastic stress and the pullback of a nutrient concentration in the current configuration, the hyperelastic variational problem is solved on the intermediate configuration after pure growth. Additionally, the coefficients of the reaction-diffusion equation determining the nutrient concentration in the current configuration depend on the elastic deformation and the deformation due to pure growth.
In this paper, we study a tumor growth model with nutrients. The contact inhibition for the tumor cells, presented in the model, results in the evolution of a congested tumor patch. We study the regularity of the tumor patch as the nutrients' diffusion strength $D$ diminishes. In particular, we show that for small $D>0$ the boundary of the tumor patch stays in a small neighborhood of the smooth tumor patch boundary obtained with $D=0$, uniformly with respect to the Hausdorff distance.
In this paper, we establish a reaction-diffusion-advection partial differential equation model to describe the growth of algae depending on both nutrients and inorganic carbon in a poorly mixed water column. Nutrients from the water bottom and inorganic carbon from the water surface form an asymmetric resource supply mechanism on the algal growth. The existence and stability of semi-trivial steady state and coexistence steady state of the model are proved, and a threshold condition for the regime shift from extinction to survival of algae is established. The influence of environmental parameters on the vertical distribution of algae is investigated in the water column. It is shown that the vertical distribution of algae can exhibit many different profiles under the joint limitation of nutrients and inorganic carbon.
In this paper we study a tumor growth model with nutrients. The model presents dynamic patch solutions due to the contact inhibition among the tumor cells. We show that when the nutrients do not diffuse and the cells do not die, the tumor density exhibits regularizing dynamics. In particular, we provide contraction estimates, exponential rate of asymptotic convergence, and boundary regularity of the tumor patch. These results are in sharp contrast to the models either with nutrient diffusion or with death rate in tumor cells.
For tumor growth, the morphological instability provides a mechanism for invasion via tumor fingering and fragmentation. This work considers the asymptotic stability of a free boundary tumor model with a periodic supply of external nutrients. The model consists of two elliptic equations describing the concentration of nutrients and the distribution of the internal pressure in the tumor tissue, respectively. The effect of the parameter $μ$ representing a measure of mitosis on the morphological stability are taken into account. It was recently established in [25] that there exists a critical value $μ_\ast$ such that the unique spherical periodic positive solution is linearly stable for $0<μ<μ_\ast$ and linearly unstable for $μ>μ_\ast$. In this paper, we further prove that the spherical periodic positive solution is asymptotically stable for $0<μ<μ_\ast$ for the fully nonlinear problem.
Existence and uniqueness of solutions for a class of models for stress-modulated growth is proven in one spatial dimension. The model features the multiplicative decomposition of the deformation gradient $F$ into an elastic part $F_e$ and a growth-related part $G$. After the transformation due to the growth process, governed by $G$, an elastic deformation described by $F_e$ is applied in order to restore the Dirichlet boundary conditions and therefore the current configuration might be stressed with a stress tensor $S$. The growth of the material at each point in the reference configuration is given by an ordinary differential equation for which the right-hand side may depend on the stress $S$ and the pull-back of a nutrient concentration in the current configuration, leading to a coupled system of ordinary differential equations.
Microbial communities routinely have several possible species compositions or community states observed for the same environmental parameters. Changes in these parameters can trigger abrupt and persistent transitions (regime shifts) between such community states. Yet little is known about the main determinants and mechanisms of multistability in microbial communities. Here we introduce and study a resource-explicit model in which microbes compete for two types of essential nutrients. We adapt game-theoretical methods of the stable matching problem to identify all possible species compositions of a microbial community. We then classify them by their resilience against three types of perturbations: fluctuations in nutrient supply, invasions by new species, and small changes of abundances of existing ones. We observe multistability and explore an intricate network of regime shifts between stable states in our model. Our results suggest that multistability requires microbial species to have different stoichiometries of essential nutrients. We also find that balanced nutrient supply promote multistability and species diversity yet make individual community states less stable.
Microbial communities routinely have several alternative stable states observed for the same environmental parameters. Sudden and irreversible transitions between these states make external manipulation of these systems more complicated. To better understand the mechanisms and origins of multistability in microbial communities, we introduce and study a model of a microbial ecosystem colonized by multiple specialist species selected from a fixed pool. Growth of each species can be limited by essential nutrients of two types, e.g. carbon and nitrogen, each represented in the environment by multiple metabolites. We demonstrate that our model has an exponentially large number of potential stable states realized for different environmental parameters. Using game theoretical methods adapted from the stable marriage problem we predict all of these states based only on ranked lists of competitive abilities of species for each of the nutrients. We show that for every set of nutrient influxes, several mutually uninvadable stable states are generally feasible and we distinguish them based upon their dynamic stability. We further explore an intricate network of discontinuous transitions (regim
In this paper, we consider a multi-layer tumor model with a periodic provision of external nutrients. The domain occupied by tumor has a different shape (flat shape) than spherical shape which has been studied widely. The important parameters are periodic external nutrients $Φ(t)$ and threshold concentration for proliferation $\widetildeσ$. In this paper, we give a complete classification about $Φ(t)$ and $\widetildeσ$ according to global stability of zero equilibrium solution or global stability of the positive periodic solution. Precisely, if $\frac{1}{T} \int_{0}^{T} Φ(t)d t\leqslant\widetildeσ$, then the zero equilibrium solution is globally stable while if $\frac{1}{T} \int_{0}^{T} Φ(t)d t>\widetildeσ$, then there exists a unique positive T-periodic solution and it is globally stable.
One of the major characteristics of living organisms is metabolic rate, which is the amount of energy produced per unit of time. When the mass of organisms increases, the metabolic rate also increases (usually as a power function of mass), but usually slower than mass. This effect is called metabolic allometric scaling. Its causes are considered unknown. The effect has important implications for individual and population organismal development. It was shown in the first part of this study, presented in a separate paper, that in the case of multicellular organisms, this effect is a consequence of natural selection and optimization of nutrient distribution between the species of a food chain, sharing resources of a common habitat. Here, in the second part that studies unicellular organisms, we discover that the same principle of natural selection guided by optimization of nutrient distribution between the species of a food chain defines also metabolic allometric scaling of unicellular organisms. To find that, we consider the metabolic properties of Amoeba proteus, fission yeast Schizosaccharomyces pombe, Escherichia coli, Bacillus subtilis, Staphylococcus. The sharing of nutrients is
The present study investigates different strategies for the treatment of a mixture of digestate from an anaerobic digester diluted and secondary effluent from a high rate algal pond. To this aim, the performance of two photo-sequencing batch reactors (PSBRs) operated at high nutrients loading rates and different solids retention times (SRTs) were compared with a semi-continuous photobioreactor (SC). Performances were evaluated in terms of wastewater treatment, biomass composition and biopolymers accumulation during 30 days of operation. PSBRs were operated at a hydraulic retention time (HRT) of 2 days and SRTs of 10 and 5 days (PSBR2-10 and PSBR2-5, respectively), whereas the semi-continuous reactor was operated at a coupled HRT/SRT of 10 days (SC10-10). Results showed that PSBR2-5 achieved the highest removal rates in terms of TN (6.7 mg L-1 d-1), TP (0.31 mg L-1 d-1), TOC (29.32 mg L-1 d-1) and TIC (3.91 mg L-1 d-1). These results were in general 3-6 times higher than the removal rates obtained in the SC10-10 (TN 29.74 mg L-1 d-1, TP 0.96 mg L-1 d-1, TOC 29.32 mg L-1 d-1 and TIC 3.91 mg L-1 d-1). Furthermore, both PSBRs were able to produce biomass up to 0.09 g L-1 d-1, more than
Crop rotation impacts on soil nutrients are typically assessed using field-averaged or single-nutrient analyses that ignore spatial heterogeneity and multivariate interactions. We propose a multivariate lattice model treating soil as a 4D tensor (space, time, and N, P, K channels). Crop rotations are represented as force vectors, with soil buffering capacity ("stiffness") varying spatially with texture. Lateral nutrient movement is introduced via kernel smoothing. Cumulative impact is quantified by Euclidean distance in N-P-K space, with significance assessed via Cramer-von Mises permutation tests. Simulating a three-year corn-soybean-wheat rotation on a 20 x 20 heterogeneous grid shows mean stress of 0.63 after one cycle, with maximum 0.91 in sandy areas. Phosphorus depletion (17.9%) exceeds nitrogen (10.8%), dominating stress in 19% of cells - obscured by single-nutrient analyses. Continuous corn increases mean stress by 41%. Cramer-von Mises tests detect significant deviation (p <= 0.002), and Moran's I (0.29-0.30) confirms spatial autocorrelation. Our framework identifies risk zones and guides site-specific management, bridging geostatistics with mechanistic crop models.
The majority of therapeutic monoclonal antibodies (mAbs) on the market are produced using Chinese Hamster Ovary (CHO) cells cultured at scale in chemically defined cell culture medium. Because of the high costs associated with mammalian cell cultures, obtaining high cell densities to produce high product titers is desired. These bioprocesses require high concentrations of nutrients in the basal media and periodically adding concentrated feed media to sustain cell growth and therapeutic protein productivity. Unfortunately, the desired or optimal nutrient concentrations of the feed media are often solubility limited due to precipitation of chemical complexes that form in the solution. Experimentally screening the various cell culture media configurations which contain 50 to 100 compounds can be expensive and laborious. This article lays the foundation for utilizing computational tools to understand precipitation of nutrients in cell culture media by studying the pairwise interactions between amino acids in thermodynamic models. Activity coefficient data for one amino acid in water and amino acid solubility data of two amino acids in water have been used to determine a single set of U
Life as we know it requires the presence of liquid water and the availability of nutrients, which are mainly based on the elements C, H, N, O, P, and S (CHNOPS) and trace metal micronutrients. We aim to understand the presence of these nutrients within atmospheres that show the presence of water cloud condensates, potentially allowing the existence of aerial biospheres. In this paper we introduce a framework of nutrient availability levels based on the presence of water condensates and the chemical state of the CHNOPS elements. These nutrient availability levels are applied to a set of atmospheric models based on different planetary surface compositions resulting in a range of atmospheric compositions. The atmospheric model is a bottom-to-top equilibrium chemistry atmospheric model which includes the atmosphere-crust interaction and the element depletion due to the formation of clouds. While the reduced forms of CNS are present at the water cloud base for most atmospheric compositions, P and metals are lacking. This indicates the potential bio-availability of CNS, while P and metals are limiting factors for aerial biospheres.
Through switching contraction patterns driving digestive flows, the small intestine balances nutrient uptake and waste removal. Complex segmentation contraction patterns are associated with higher nutrient uptake, while peristaltic contractions primarily serve to flush out unused remnants. However, the impact of the Non-Newtonian behavior of digestive fluids on the efficacy of these contraction patterns remains unclear. Here, we present finite-element simulations that model nutrient transport, diffusion, and uptake within both segmentation and peristaltic contractions along the small intestine for Newtonian and Non-Newtonian fluids. Our simulations reveal that diffusion plays a key role in uptake, with nutrient absorption directly linked to fluid viscosity, which governs molecular diffusivity. Further, we present an analytical prediction for uptake in peristalsis as a function of molecular diffusivity, which aligns closely with our simulation results. Notably, we find that shear-thinning properties of Non-Newtonian fluids enhance nutrient uptake, particularly in segmentation contractions compared to peristalsis. Our results demonstrate the fluid dynamical principles underlying inte
Biofilm infections on medical implants are difficult to eradicate because insufficient nutrient availability promotes antibiotic-tolerant persister cells that survive treatment and reseed growth. Existing mathematical models usually omit nutrient-dependent phenotypic switching between proliferative and persister states. Without this mechanism, models cannot capture how environmental conditions control the balance between active growth and dormancy, which is central to biofilm persistence. We present a continuum model that couples nutrient transport with the dynamics of proliferative bacteria, persisters, dead cells, and extracellular polymeric substances. The switching rates between proliferative and persister phenotypes depend on local nutrient concentration through two thresholds, enabling adaptation across nutrient-poor, intermediate, and nutrient-rich regimes. Simulations show that nutrient limitation produces a high and sustained proportion of persister cells even when biomass is reduced, whereas nutrient-rich conditions support reversion to proliferative growth and lead to greater biomass. The model also predicts that persister populations peak at times that vary with nutrien
Sessile ciliates attach to substrates and generate feeding currents to capture passing particulates and dissolved nutrients. Optimal ciliary activity that maximizes nutrient flux at the cell surface while minimizing the rate of hydrodynamic energy dissipation is well characterized in uniform nutrient fields. However, it is unclear how ciliary motion should change when nutrients are non-uniform or patchy. To address this question, we modeled the sessile ciliate and feeding currents using the spherical envelope model, and used an unsteady advection-diffusion equation to describe the nutrient scalar field. In the absence of flows, we calculated the diffusive nutrient uptake analytically in linear nutrient gradients and found no advantage over uptake in uniform fields. With ciliary activity driving feeding currents, we used a spectral method to solve for the unsteady nutrient concentration. We found that, when the axis of symmetry of the ciliary motion is aligned with the nutrient gradient, nutrient uptake at the cell surface increases steadily over time, with highest uptake achieved by the treadmill ciliary motion which is optimal in uniform fields as well. The associated nutrient upt
An early, non-invasive, and on-site detection of nutrient deficiencies is critical to enable timely actions to prevent major losses of crops caused by lack of nutrients. While acquiring labeled data is very expensive, collecting images from multiple views of a crop is straightforward. Despite its relevance for practical applications, unsupervised domain adaptation where multiple views are available for the labeled source domain as well as the unlabeled target domain is an unexplored research area. In this work, we thus propose an approach that leverages multiple camera views in the source and target domain for unsupervised domain adaptation. We evaluate the proposed approach on two nutrient deficiency datasets. The proposed method achieves state-of-the-art results on both datasets compared to other unsupervised domain adaptation methods. The dataset and source code are available at https://github.com/jh-yi/MV-Match.