Mathematical modeling of colony formation in algal blooms: phenotypic plasticity in cyanobacteria

In this paper, we analyzed a mathematical model of algal-grazer dynamics, including the effect of colony formation, which is an example of phenotypic plasticity. The model consists of three variables, which correspond to the biomasses of unicellular algae, colonial algae, and herbivorous zooplankton. Among these organisms, colonial algae are the main components of algal blooms. This aquatic system has two stable attractors, which can be identified as a zooplankton-dominated (ZD) state and an algal-dominated (AD) state, respectively. Assuming that the handling time of zooplankton on colonial algae increases with the colonial algae biomass, we discovered that bistability can occur within the model system. The applicability of alternative stable states in algae-grazer dynamics as a framework for explaining the algal blooms in real lake ecosystems, thus, seems to depend on whether the assumption mentioned above is met in natural circumstances.

Mathematical Modeling to Assess the Drivers of the Recent Emergence of Typhoid Fever in Blantyre, Malawi

BACKGROUND: Multiyear epidemics of Salmonella enterica serovar Typhi have been reported from countries across eastern and southern Africa in recent years. In Blantyre, Malawi, a dramatic increase in typhoid fever cases has recently occurred, and may be linked to the emergence of the H58 haplotype. Strains belonging to the H58 haplotype often exhibit multidrug resistance and may have a fitness advantage relative to other Salmonella Typhi strains.

METHODS: To explore hypotheses for the increased number of typhoid fever cases in Blantyre, we fit a mathematical model to culture-confirmed cases of Salmonella enterica infections at Queen Elizabeth Central Hospital, Blantyre. We explored 4 hypotheses: (1) an increase in the basic reproductive number (R0) in response to increasing population density; (2) a decrease in the incidence of cross-immunizing infection with Salmonella Enteritidis; (3) an increase in the duration of infectiousness due to failure to respond to first-line antibiotics; and (4) an increase in the transmission rate following the emergence of the H58 haplotype.

RESULTS: Increasing population density or decreasing cross-immunity could not fully explain the observed pattern of typhoid emergence in Blantyre, whereas models allowing for an increase in the duration of infectiousness and/or the transmission rate of typhoid following the emergence of the H58 haplotype provided a good fit to the data.

CONCLUSIONS: Our results suggest that an increase in the transmissibility of typhoid due to the emergence of drug resistance associated with the H58 haplotype may help to explain recent outbreaks of typhoid in Malawi and similar settings in Africa.

Genetics of Healthy Aging in Europe.The EU-Integrated Project GEHA (GEnetics of Healthy Aging)

The aim of the 5-year European Union (EU)-Integrated Project GEnetics of Healthy Aging (GEHA), constituted by 25 partners (24 from Europe plus the Beijing Genomics Institute from China), is to identify genes involved in healthy aging and longevity, which allow individuals to survive to advanced old age in good cognitive and physical function and in the absence of major age-related diseases. To achieve this aim a coherent, tightly integrated program of research that unites demographers, geriatricians, geneticists, genetic epidemiologists, molecular biologists, bioinfomaticians, and statisticians has been set up. The working plan is to: (a) collect DNA and information on the health status from an unprecedented number of long-lived 90+ sibpairs (n = 2650) and of younger ethnically matched controls (n = 2650) from 11 European countries; (b) perform a genome-wide linkage scannning in all the sibpairs (a total of 5300 individuals); this investigation will be followed by linkage disequilibrium mapping (LD mapping) of the candidate chromosomal regions; (c) study in cases (i.e., the 2650 probands of the sibpairs) and controls (2650 younger people), genomic regions (chromosome 4, D4S1564, chromosome 11, 11.p15.5) which were identified in previous studies as possible candidates to harbor longevity genes; (d) genotype all recruited subjects for apoE polymorphisms; and (e) genotype all recruited subjects for inherited as well as epigenetic variability of the mitochondrial DNA (mtDNA). The genetic analysis will be performed by 9 high-throughput platforms, within the framework of centralized databases for phenotypic, genetic, and mtDNA data. Additional advanced approaches (bioinformatics, advanced statistics, mathematical modeling, functional genomics and proteomics, molecular biology, molecular genetics) are envisaged to identify the gene variant(s) of interest. The experimental design will also allow (a) to identify gender-specific genes involved in healthy aging and longevity in women and men stratified for ethnic and geographic origin and apoE genotype; (b) to perform a longitudinal survival study to assess the impact of the identified genetic loci on 90+ people mortality; and (c) to develop mathematical and statistical models capable of combining genetic data with demographic characteristics, health status, socioeconomic factors, lifestyle habits.

The Origin of Phragmoplast Asymmetry

The phragmoplast coordinates cytokinesis in plants [1]. It directs vesicles to the midzone, the site where they coalesce to form the new cell plate. Failure in phragmoplast function results in aborted or incomplete cytokinesis leading to embryo lethality, morphological defects, or multinucleate cells [2, 3]. The asymmetry of vesicular traffic is regulated by microtubules [1, 4, 5, 6], and the current model suggests that this asymmetry is established and maintained through treadmilling of parallel microtubules. However, we have analyzed the behavior of microtubules in the phragmoplast using live-cell imaging coupled with mathematical modeling and dynamic simulations and report that microtubules initiate randomly in the phragmoplast and that the majority exhibit dynamic instability with higher turnover rates nearer to the midzone. The directional transport of vesicles is possible because the majority of the microtubules polymerize toward the midzone. Here, we propose the first inclusive model where microtubule dynamics and phragmoplast asymmetry are consistent with the localization and activity of proteins known to regulate microtubule assembly and disassembly.

Contemporary environmental issues of landfill leachate: assessment & remedies

Landfills are the primary option for waste disposal all over the world. Most of the landfill sites across the world are old and are not engineered to prevent contamination of the underlying soil and groundwater by the toxic leachate. The pollutants from landfill leachate have accumulative and detrimental effect on the ecology and food chains leading to carcinogenic effects, acute toxicity and genotoxicity among human beings. Management of this highly toxic leachate presents a challenging problem to the regulatory authorities who have set specific regulations regarding maximum limits of contaminants in treated leachate prior to disposal into the environment to ensure minimal environmental impact. There are different stages of leachate management such as monitoring of its formation and flow into the environment, identification of hazards associated with it and its treatment prior to disposal into the environment. This review focuses on: (i) leachate composition, (ii) Plume migration, (iii) Contaminant fate, (iv) Leachate plume monitoring techniques, (v) Risk assessment techniques, Hazard rating methods, mathematical modeling, and (vi) Recent innovations in leachate treatment technologies. However, due to seasonal fluctuations in leachate composition, flow rate and leachate volume, the management approaches cannot be stereotyped. Every scenario is unique and the strategy will vary accordingly. This paper lays out the choices for making an educated guess leading to the best management option.

Quantitative regulation of FLC via coordinated transcriptional initiation and elongation

The basis of quantitative regulation of gene expression is still poorly understood. In Arabidopsis thaliana, quantitative variation in expression of FLOWERING LOCUS C (FLC) influences the timing of flowering. In ambient temperatures, FLC expression is quantitatively modulated by a chromatin silencing mechanism involving alternative polyadenylation of antisense transcripts. Investigation of this mechanism unexpectedly showed that RNA polymerase II (Pol II) occupancy changes at FLC did not reflect RNA fold changes. Mathematical modeling of these transcriptional dynamics predicted a tight coordination of transcriptional initiation and elongation. This prediction was validated by detailed measurements of total and chromatin-bound FLC intronic RNA, a methodology appropriate for analyzing elongation rate changes in a range of organisms. Transcription initiation was found to vary ∼ 25-fold with elongation rate varying ∼ 8- to 12-fold. Premature sense transcript termination contributed very little to expression differences. This quantitative variation in transcription was coincident with variation in H3K36me3 and H3K4me2 over the FLC gene body. We propose different chromatin states coordinately influence transcriptional initiation and elongation rates and that this coordination is likely to be a general feature of quantitative gene regulation in a chromatin context.

Nonisothermal moisture transport in hygroscopic building materials: modeling for the determination of moisture transport coefficients

A mathematical model for calculating the nonisothermal moisture transfer in building materials is presented in the article. The coupled heat and moisture transfer problem was modeled. Vapor content and temperature were chosen as principal driving potentials. The coupled equations were solved by an analytical method, which consists of applying the Laplace transform technique and the Transfer Function Method. A new experimental methodology for determining the temperature gradient coefficient for building materials was also proposed. Both the moisture diffusion coefficient and the temperature gradient coefficient for building material were experimentally evaluated. Using the measured moisture transport coefficients, the temperature and vapor content distribution inside building materials were predicted by the new model. The results were compared with experimental data. A good agreement was obtained.

Modeling the geometric characteristics of five-dimensionally woven composites

An analytical modeling approach for the prediction of the geometric characteristics of five-dimensional (5D) woven composites has been formulated. The model is driven by readily available data including the weaving parameters and constituent material properties. The new model calculates the individual proportions of fiber in each direction, areal density, overall fiber volume fraction, and laminate thickness. This information is useful for the engineer in the design and manufacture of 5D woven composites. In addition the present model outputs the mathematical definition of the 5D woven composite unit cell, which could be implemented as the geometric input for a downstream analytical model that is capable of predicting the elastic stiffness of 5D woven composites. Input parameters have been sourced from existing published work and the subsequent predictions made by the model are compared with the available experimental data on 5D woven composites.

Modeling colloid deposition on a protein layer adsorbed to iron-oxide-coated sand

Our recent study reported that conformation change of granule-associated Bovine Serum Albumin (BSA) may influence the role of the protein controlling colloid deposition in porous media (Flynn et al., 2012). The present study conceptualized the observed phenomena with an ellipsoid morphology model, describing BSA as an ellipsoid taking a side-on or end-on conformation on granular surface, and identified the following processes: (1) at low adsorbed concentrations, BSA exhibited a side-on conformation blocking colloid deposition; (2) at high adsorbed concentrations, BSA adapted to an end-on conformation promoted colloid deposition; and (3) colloid deposition on the BSA layer may progressively generate end-on molecules (sites) by conformation change of side-on BSA, resulting in sustained increasing deposition rates. Generally, the protein layer lowered colloid attenuation by the porous medium, suggesting the overall effect of BSA was inhibitory at the experimental time scale. A mathematical model was developed to interpret the ripening curves. Modeling analysis identified the site generation efficiency of colloid as a control on the ripening rate (declining rate in colloid concentrations), and this efficiency was higher for BSA adsorbed from a more dilute BSA solution. © 2012 Elsevier B.V.

A Lagrangian approach for modeling road collisions using second-order models of traffic flow

In [M. Herty, A. Klein, S. Moutari, V. Schleper, and G. Steinaur, IMA J. Appl. Math., 78(5), 1087–1108, 2013] and [M. Herty and V. Schleper, ZAMM J. Appl. Math. Mech., 91, 763–776, 2011], a macroscopic approach, derived from fluid-dynamics models, has been introduced to infer traffic conditions prone to road traffic collisions along highways’ sections. In these studies, the governing equations are coupled within an Eulerian framework, which assumes fixed interfaces between the models. A coupling in Lagrangian coordinates would enable us to get rid of this (not very realistic) assumption. In this paper, we investigate the well-posedness and the suitability of the coupling of the governing equations within the Lagrangian framework. Further, we illustrate some features of the proposed approach through some numerical simulations.

Bdellovibrio predation in the presence of decoys:Three-way bacterial interactions revealed by mathematical and experimental analyses

Bdellovibrio bacteriovorus is a small, gram-negative, motile bacterium that preys upon other gram-negative bacteria, including several known human pathogens. Its predation efficiency is usually studied in pure cultures containing solely B. bacteriovorus and a suitable prey. However, in natural environments, as well as in any possible biomedical uses as an antimicrobial, Bdellovibrio is predatory in the presence of diverse decoys, including live nonsusceptible bacteria, eukaryotic cells, and cell debris. Here we gathered and mathematically modeled data from three-member cultures containing predator, prey, and nonsusceptible bacterial decoys. Specifically, we studied the rate of predation of planktonic late-log-phase Escherichia coli S17-1 prey by B. bacteriovorus HD100, both in the presence and in the absence of Bacillus subtilis nonsporulating strain 671, which acted as a live bacterial decoy. Interestingly, we found that although addition of the live Bacillus decoy did decrease the rate of Bdellovibrio predation in liquid cultures, this addition also resulted in a partially compensatory enhancement of the availability of prey for predation. This effect resulted in a higher final yield of Bdellovibrio than would be predicted for a simple inert decoy. Our mathematical model accounts for both negative and positive effects of predator-prey-decoy interactions in the closed batch environment. In addition, it informs considerations for predator dosing in any future therapeutic applications and sheds some light on considerations for modeling the massively complex interactions of real mixed bacterial populations in nature.