Non-linear growth models for bullfrog tadpoles

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Development and Application of Predictive Models for Survival, Growth, and Death of Enteric Pathogens in Leafy Greens Supply Chain

Leafy greens are essential part of a healthy diet. Because of their health benefits, production and consumption of leafy greens has increased considerably in the U.S. in the last few decades. However, leafy greens are also associated with a large number of foodborne disease outbreaks in the last few years. The overall goal of this dissertation was to use the current knowledge of predictive models and available data to understand the growth, survival, and death of enteric pathogens in leafy greens at pre- and post-harvest levels. Temperature plays a major role in the growth and death of bacteria in foods. A growth-death model was developed for Salmonella and Listeria monocytogenes in leafy greens for varying temperature conditions typically encountered during supply chain. The developed growth-death models were validated using experimental dynamic time-temperature profiles available in the literature. Furthermore, these growth-death models for Salmonella and Listeria monocytogenes and a similar model for E. coli O157:H7 were used to predict the growth of these pathogens in leafy greens during transportation without temperature control. Refrigeration of leafy greens meets the purposes of increasing their shelf-life and mitigating the bacterial growth, but at the same time, storage of foods at lower temperature increases the storage cost. Nonlinear programming was used to optimize the storage temperature of leafy greens during supply chain while minimizing the storage cost and maintaining the desired levels of sensory quality and microbial safety. Most of the outbreaks associated with consumption of leafy greens contaminated with E. coli O157:H7 have occurred during July-November in the U.S. A dynamic system model consisting of subsystems and inputs (soil, irrigation, cattle, wildlife, and rainfall) simulating a farm in a major leafy greens producing area in California was developed. The model was simulated incorporating the events of planting, irrigation, harvesting, ground preparation for the new crop, contamination of soil and plants, and survival of E. coli O157:H7. The predictions of this system model are in agreement with the seasonality of outbreaks. This dissertation utilized the growth, survival, and death models of enteric pathogens in leafy greens during production and supply chain.

Effect of female sex hormones on Chlamydia trachomatis growth and gene expression

Transmissible diseases are re-emerging as a global problem, with Sexually Transmitted Diseases (STDs) becoming endemic. Chlamydia trachomatis is the leading cause of bacterially-acquired STD worldwide, with the Australian cost of infection estimated at $90 - $160 million annually. Studies using animal models of genital tract Chlamydia infection suggested that the hormonal status of the genital tract epithelium at the time of exposure may influence the outcome of infection. Oral contraceptive use also increased the risk of contracting chlamydial infections compared to women not using contraception. Generally it was suggested that the outcome of chlamydial infection is determined in part by the hormonal status of the epithelium at the time of exposure. Using the human endolmetrial cell line ECC-1 this study investigated the effects of C. trachomatis serovar D infection, in conjunction with the female sex hormones, 17β-estradiol and progesterone, on chlamydial gene expression. While previous studies have examined the host response, this is the first study to examine C.trachomatis gene expression under different hormonal conditions. We have highlighted a basic model of C. trachomatis gene regulation in the presence of steroid hormones by identifying 60 genes that were regulated by addition of estradiol and/or progesterone. In addition, the third chapter of this thesis discussed and compared the significance of the current findings in the context of data from other research groups to improve our understanding of the molecular basis of chlamydial persistence under hormonal different conditions. In addition, this study analysed the effects of these female sex hormones and C. trachomatis Serovar D infection, on host susceptibility and bacterial growth. Our results clearly demonstrated that addition of steroid hormones not only had a great impact on the level of infectivity of epithelial cells with C.trachomatis serovar D, but also the morphology of chlamydial inclusions was affected by hormone supplementation.

In vivo performance of melimine as an antimicrobial coating for contact lenses in models of CLARE and CLPU

Purpose: One strategy to minimize bacteria-associated adverse responses such as microbial keratitis, contact lens–induced acute red eye (CLARE), and contact lens induced peripheral ulcers (CLPUs) that occur with contact lens wear is the development of an antimicrobial or antiadhesive contact lens. Cationic peptides represent a novel approach for the development of antimicrobial lenses.---------- Methods: A novel cationic peptide, melimine, was covalently incorporated into silicone hydrogel lenses. Confirmation tests to determine the presence of peptide and anti-microbial activity were performed. Cationic lenses were then tested for their ability to prevent CLPU in the Staphylococcus aureus rabbit model and CLARE in the Pseudomonas aeruginosa guinea pig model. ---------- Results: In the rabbit model of CLPU, melimine-coated lenses resulted in significant reductions in ocular symptom scores and in the extent of corneal infiltration (P < 0.05). Evaluation of the performance of melimine lenses in the CLARE model showed significant improvement in all ocular response parameters measured, including the percentage of eyes with corneal infiltrates, compared with those observed in the eyes fitted with the control lens (P ≤ 0.05). ---------- Conclusions: Cationic coating of contact lenses with the peptide melimine may represent a novel method of prevention of bacterial growth on contact lenses and consequently result in reduction of the incidence and severity of adverse responses due to Gram-positive and -negative bacteria during lens wear.

Morphology-based model for predicting osteocyte-like cell line (MLO-Y4) growth process

This work has led to the development of empirical mathematical models to quantitatively predicate the changes of morphology in osteocyte-like cell lines (MLO-Y4) in culture. MLO-Y4 cells were cultured at low density and the changes in morphology recorded over 11 hours. Cell area and three dimensional shape features including aspect ratio, circularity and solidity were then determined using widely accepted image analysis software (ImageJTM). Based on the data obtained from the imaging analysis, mathematical models were developed using the non-linear regression method. The developed mathematical models accurately predict the morphology of MLO-Y4 cells for different culture times and can, therefore, be used as a reference model for analyzing MLO-Y4 cell morphology changes within various biological/mechanical studies, as necessary.

Quantifying uncertainty in parameter estimates for stochastic models of collective cell spreading using approximate Bayesian computation

Wound healing and tumour growth involve collective cell spreading, which is driven by individual motility and proliferation events within a population of cells. Mathematical models are often used to interpret experimental data and to estimate the parameters so that predictions can be made. Existing methods for parameter estimation typically assume that these parameters are constants and often ignore any uncertainty in the estimated values. We use approximate Bayesian computation (ABC) to estimate the cell diffusivity, D, and the cell proliferation rate, λ, from a discrete model of collective cell spreading, and we quantify the uncertainty associated with these estimates using Bayesian inference. We use a detailed experimental data set describing the collective cell spreading of 3T3 fibroblast cells. The ABC analysis is conducted for different combinations of initial cell densities and experimental times in two separate scenarios: (i) where collective cell spreading is driven by cell motility alone, and (ii) where collective cell spreading is driven by combined cell motility and cell proliferation. We find that D can be estimated precisely, with a small coefficient of variation (CV) of 2–6%. Our results indicate that D appears to depend on the experimental time, which is a feature that has been previously overlooked. Assuming that the values of D are the same in both experimental scenarios, we use the information about D from the first experimental scenario to obtain reasonably precise estimates of λ, with a CV between 4 and 12%. Our estimates of D and λ are consistent with previously reported values; however, our method is based on a straightforward measurement of the position of the leading edge whereas previous approaches have involved expensive cell counting techniques. Additional insights gained using a fully Bayesian approach justify the computational cost, especially since it allows us to accommodate information from different experiments in a principled way.

Mathematical modelling of human corneal oxygenation and cell kinetics in corneal epithelial wound healing

Healthy transparent cornea depends upon the regulation of fluid, nutrient and oxygen transport through the tissue to sustain cell metabolism and other critical processes for normal functioning. This research considers the corneal geometry and investigates oxygen distribution using a two-dimensional Monod kinetic model, showing that previous studies make assumptions that lead to predictions of near-anoxic levels of oxygen tension in the limbal regions of the cornea. It also considers the comparison of experimental spatial and temporal data with the predictions of novel mathematical models with respect to distributed mitotic rates during corneal epithelial wound healing.

Modelling within-host spatiotemporal dynamics of invasive bacterial disease

Mechanistic determinants of bacterial growth, death, and spread within mammalian hosts cannot be fully resolved studying a single bacterial population. They are also currently poorly understood. Here, we report on the application of sophisticated experimental approaches to map spatiotemporal population dynamics of bacteria during an infection. We analyzed heterogeneous traits of simultaneous infections with tagged Salmonella enterica populations (wild-type isogenic tagged strains [WITS]) in wild-type and gene-targeted mice. WITS are phenotypically identical but can be distinguished and enumerated by quantitative PCR, making it possible, using probabilistic models, to estimate bacterial death rate based on the disappearance of strains through time. This multidisciplinary approach allowed us to establish the timing, relative occurrence, and immune control of key infection parameters in a true host-pathogen combination. Our analyses support a model in which shortly after infection, concomitant death and rapid bacterial replication lead to the establishment of independent bacterial subpopulations in different organs, a process controlled by host antimicrobial mechanisms. Later, decreased microbial mortality leads to an exponential increase in the number of bacteria that spread locally, with subsequent mixing of bacteria between organs via bacteraemia and further stochastic selection. This approach provides us with an unprecedented outlook on the pathogenesis of S. enterica infections, illustrating the complex spatial and stochastic effects that drive an infectious disease. The application of the novel method that we present in appropriate and diverse host-pathogen combinations, together with modelling of the data that result, will facilitate a comprehensive view of the spatial and stochastic nature of within-host dynamics. © 2008 Grant et al.

Interception of host angiogenic signalling limits mycobacterial growth.

Pathogenic mycobacteria induce the formation of complex cellular aggregates called granulomas that are the hallmark of tuberculosis. Here we examine the development and consequences of vascularization of the tuberculous granuloma in the zebrafish-Mycobacterium marinum infection model, which is characterized by organized granulomas with necrotic cores that bear striking resemblance to those of human tuberculosis. Using intravital microscopy in the transparent larval zebrafish, we show that granuloma formation is intimately associated with angiogenesis. The initiation of angiogenesis in turn coincides with the generation of local hypoxia and transcriptional induction of the canonical pro-angiogenic molecule Vegfaa. Pharmacological inhibition of the Vegf pathway suppresses granuloma-associated angiogenesis, reduces infection burden and limits dissemination. Moreover, anti-angiogenic therapies synergize with the first-line anti-tubercular antibiotic rifampicin, as well as with the antibiotic metronidazole, which targets hypoxic bacterial populations. Our data indicate that mycobacteria induce granuloma-associated angiogenesis, which promotes mycobacterial growth and increases spread of infection to new tissue sites. We propose the use of anti-angiogenic agents, now being used in cancer regimens, as a host-targeting tuberculosis therapy, particularly in extensively drug-resistant disease for which current antibiotic regimens are largely ineffective.

The inoculum effect and band-pass bacterial response to periodic antibiotic treatment.

The inoculum effect (IE) refers to the decreasing efficacy of an antibiotic with increasing bacterial density. It represents a unique strategy of antibiotic tolerance and it can complicate design of effective antibiotic treatment of bacterial infections. To gain insight into this phenomenon, we have analyzed responses of a lab strain of Escherichia coli to antibiotics that target the ribosome. We show that the IE can be explained by bistable inhibition of bacterial growth. A critical requirement for this bistability is sufficiently fast degradation of ribosomes, which can result from antibiotic-induced heat-shock response. Furthermore, antibiotics that elicit the IE can lead to 'band-pass' response of bacterial growth to periodic antibiotic treatment: the treatment efficacy drastically diminishes at intermediate frequencies of treatment. Our proposed mechanism for the IE may be generally applicable to other bacterial species treated with antibiotics targeting the ribosomes.

Ocean acidification with (de)eutrophication will alter future phytoplankton growth and succession.

Human activity causes ocean acidification (OA) though the dissolution of anthropogenically generated CO2 into seawater, and eutrophication through the addition of inorganic nutrients. Eutrophication increases the phytoplankton biomass that can be supported during a bloom, and the resultant uptake of dissolved inorganic carbon during photosynthesis increases water-column pH (bloom-induced basification). This increased pH can adversely affect plankton growth. With OA, basification commences at a lower pH. Using experimental analyses of the growth of three contrasting phytoplankton under different pH scenarios, coupled with mathematical models describing growth and death as functions of pH and nutrient status, we show how different conditions of pH modify the scope for competitive interactions between phytoplankton species. We then use the models previously configured against experimental data to explore how the commencement of bloom-induced basification at lower pH with OA, and operating against a background of changing patterns in nutrient loads, may modify phytoplankton growth and competition. We conclude that OA and changed nutrient supply into shelf seas with eutrophication or de-eutrophication (the latter owing to pollution control) has clear scope to alter phytoplankton succession, thus affecting future trophic dynamics and impacting both biogeochemical cycling and fisheries.

Design of neuro-fuzzy models by evolutionary and gradient-based algorithms

All systems found in nature exhibit, with different degrees, a nonlinear behavior. To emulate this behavior, classical systems identification techniques use, typically, linear models, for mathematical simplicity. Models inspired by biological principles (artificial neural networks) and linguistically motivated (fuzzy systems), due to their universal approximation property, are becoming alternatives to classical mathematical models. In systems identification, the design of this type of models is an iterative process, requiring, among other steps, the need to identify the model structure, as well as the estimation of the model parameters. This thesis addresses the applicability of gradient-basis algorithms for the parameter estimation phase, and the use of evolutionary algorithms for model structure selection, for the design of neuro-fuzzy systems, i.e., models that offer the transparency property found in fuzzy systems, but use, for their design, algorithms introduced in the context of neural networks. A new methodology, based on the minimization of the integral of the error, and exploiting the parameter separability property typically found in neuro-fuzzy systems, is proposed for parameter estimation. A recent evolutionary technique (bacterial algorithms), based on the natural phenomenon of microbial evolution, is combined with genetic programming, and the resulting algorithm, bacterial programming, advocated for structure determination. Different versions of this evolutionary technique are combined with gradient-based algorithms, solving problems found in fuzzy and neuro-fuzzy design, namely incorporation of a-priori knowledge, gradient algorithms initialization and model complexity reduction.

Numerical model of solid phase transformations governed by nucleation and growth: microstructure development during isothermal crystallization

A simple numerical model which calculates the kinetics of crystallization involving randomly distributed nucleation and isotropic growth is presented. The model can be applied to different thermal histories and no restrictions are imposed on the time and the temperature dependences of the nucleation and growth rates. We also develop an algorithm which evaluates the corresponding emerging grain-size distribution. The algorithm is easy to implement and particularly flexible, making it possible to simulate several experimental conditions. Its simplicity and minimal computer requirements allow high accuracy for two- and three-dimensional growth simulations. The algorithm is applied to explore the grain morphology development during isothermal treatments for several nucleation regimes. In particular, thermal nucleation, preexisting nuclei, and the combination of both nucleation mechanisms are analyzed. For the first two cases, the universal grain-size distribution is obtained. The high accuracy of the model is stated from its comparison to analytical predictions. Finally, the validity of the Kolmogorov-Johnson-Mehl-Avrami model SSSR, is verified for all the cases studied

A review and classification of the existing models of cyanobacteria

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling of cyanobacteria in freshwaters is an important tool for understanding their population dynamics and predicting bloom occurrence in lakes and rivers. In this paper existing key models of cyanobacteria are reviewed, evaluated and classified. Two major groups emerge: deterministic mathematical and artificial neural network models. Mathematical models can be further subcategorized into those models concerned with impounded water bodies and those concerned with rivers. Most existing models focus on a single aspect such as the growth of transport mechanisms, but there are a few models which couple both.