Special Issue on Spatial Moment Techniques for Modelling Biological Processes

Over the last two decades, there has been an increasing awareness of, and interest in, the use of spatial moment techniques to provide insight into a range of biological and ecological processes. Models that incorporate spatial moments can be viewed as extensions of mean-field models. These mean-field models often consist of systems of classical ordinary differential equations and partial differential equations, whose derivation, at some point, hinges on the simplifying assumption that individuals in the underlying stochastic process encounter each other at a rate that is proportional to the average abundance of individuals. This assumption has several implications, the most striking of which is that mean-field models essentially neglect any impact of the spatial structure of individuals in the system. Moment dynamics models extend traditional mean-field descriptions by accounting for the dynamics of pairs, triples and higher n-tuples of individuals. This means that moment dynamics models can, to some extent, account for how the spatial structure affects the dynamics of the system in question.

Simulating the evolution of glyphosate resistance in grains farming in northern Australia.

Background and Aims: The evolution of resistance to herbicides is a substantial problem in contemporary agriculture. Solutions to this problem generally consist of the use of practices to control the resistant population once it evolves, and/or to institute preventative measures before populations become resistant. Herbicide resistance evolves in populations over years or decades, so predicting the effectiveness of preventative strategies in particular relies on computational modelling approaches. While models of herbicide resistance already exist, none deals with the complex regional variability in the northern Australian sub-tropical grains farming region. For this reason, a new computer model was developed. Methods: The model consists of an age- and stage-structured population model of weeds, with an existing crop model used to simulate plant growth and competition, and extensions to the crop model added to simulate seed bank ecology and population genetics factors. Using awnless barnyard grass (Echinochloa colona) as a test case, the model was used to investigate the likely rate of evolution under conditions expected to produce high selection pressure. Key Results: Simulating continuous summer fallows with glyphosate used as the only means of weed control resulted in predicted resistant weed populations after approx. 15 years. Validation of the model against the paddock history for the first real-world glyphosate-resistant awnless barnyard grass population shows that the model predicted resistance evolution to within a few years of the real situation. Conclusions: This validation work shows that empirical validation of herbicide resistance models is problematic. However, the model simulates the complexities of sub-tropical grains farming in Australia well, and can be used to investigate, generate and improve glyphosate resistance prevention strategies.

Molecular dynamics simulations of sliding in an Fe-Cu tribopair system

Experimental evidence suggests that high strain rates, stresses, strains and temperatures are experienced near sliding interfaces. The associated microstructural changes are due to several dynamic an interacting phenomena. 3D non-equilibrium molecular dynamics (MD) simulations of sliding were conducted with the aim of understanding the dynamic processes taking place in crystalline tribopairs, with a focus on plastic deformation and microstructural evolution. Embedded atom potentials were employed for simulating sliding of an Fe-Cu tribopair. Sliding velocity, crystal orientation and presence of lattice defects were some of the variables in these simulations. Extensive plastic deformation involving dislocation and twin activity, dynamic recrystallization, amorphization and/or nanocrystallization, mechanical mixing and material transfer were observed. Mechanical mixing in the vicinity of the sliding interface was observed even in the Fe-Cu system, which would cluster under equilibrium conditions, hinting at the ballistic nature of the process. Flow localization was observed at high velocities implying the possible role of adiabatic heating. The presence of preexisting defects (such as dislocations and interfaces) played a pivotal role in determining friction and microstructural evolution. The study also shed light on the relationship between adhesion and plastic deformation, and friction. Comparisons with experiments suggest that such simulations can indeed provide valuable insights that are difficult to obtain from experiments.

Complete and interpretable conformance checking of business processes

This article presents a method for checking the conformance between an event log capturing the actual execution of a business process, and a model capturing its expected or normative execution. Given a business process model and an event log, the method returns a set of statements in natural language describing the behavior allowed by the process model but not observed in the log and vice versa. The method relies on a unified representation of process models and event logs based on a well-known model of concurrency, namely event structures. Specifically, the problem of conformance checking is approached by folding the input event log into an event structure, unfolding the process model into another event structure, and comparing the two event structures via an error-correcting synchronized product. Each behavioral difference detected in the synchronized product is then verbalized as a natural language statement. An empirical evaluation shows that the proposed method scales up to real-life datasets while producing more concise and higher-level difference descriptions than state-of-the-art conformance checking methods.

Impacts of improved grazing land management on sediment yields, Part 1: Hillslope processes.

Summary Poor land condition resulting from unsustainable grazing practices can reduce enterprise profitability and increase water, sediment and associated nutrient runoff from properties and catchments. This paper presents the results of a 6 year field study that used a series of hillslope flume experiments to evaluate the impact of improved grazing land management (GLM) on hillslope runoff and sediment yields. The study was carried out on a commercial grazing property in a catchment draining to the Burdekin River in northern Australia. During this study average ground cover on hillslopes increased from ~35% to ~75%, although average biomass and litter levels are still relatively low for this landscape type (~60 increasing to 1100 kg of dry matter per hectare). Pasture recovery was greatest on the upper and middle parts of hillslopes. Areas that did not respond to the improved grazing management had <10% cover and were on the lower slopes associated with the location of sodic soil and the initiation of gullies. Comparison of ground cover changes and soil conditions with adjacent properties suggest that grazing management, and not just improved rainfall conditions, were responsible for the improvements in ground cover in this study. The ground cover improvements resulted in progressively lower runoff coefficients for the first event in each wet season, however, runoff coefficients were not reduced at the annual time scale. The hillslope annual sediment yields declined by ~70% on two out of three hillslopes, although where bare patches (with <10% cover) were connected to gullies and streams, annual sediment yields increased in response to higher rainfall in latter years of the study. It appears that bare patches are the primary source areas for both runoff and erosion on these hillslopes. Achieving further reductions in runoff and erosion in these landscapes may require management practices that improve ground cover and biomass in bare areas, particularly when they are located adjacent to concentrated drainage lines.

Characterization of the mitochondrial active-site serine protein LACTB: filaments in the mitochondrial intermembrane space

Mitochondria have evolved from endosymbiotic alpha-proteobacteria. During the endosymbiotic process early eukaryotes dumped the major component of the bacterial cell wall, the peptidoglycan layer. Peptidoglycan is synthesized and maintained by active-site serine enzymes belonging to the penicillin-binding protein and the β-lactamase superfamily. Mammals harbor a protein named LACTB that shares sequence similarity with bacterial penicillin-binding proteins and β-lactamases. Since eukaryotes lack the synthesis machinery for peptidoglycan, the physiological role of LACTB is intriguing. Recently, LACTB has been validated in vivo to be causative for obesity, suggesting that LACTB is implicated in metabolic processes. The aim of this study was to investigate the phylogeny, structure, biochemistry and cell biology of LACTB in order to elucidate its physiological function. Phylogenetic analysis revealed that LACTB has evolved from penicillin binding-proteins present in the bacterial periplasmic space. A structural model of LACTB indicates that LACTB shares characteristic features common to all penicillin-binding proteins and β-lactamases. Recombinat LACTB protein expressed in E. coli was recovered in significant quantities. Biochemical and cell biology studies showed that LACTB is a soluble protein localized in the mitochondrial intermembrane space. Further analysis showed that LACTB preprotein underwent proteolytic processing disclosing an N-terminal tetrapeptide motif also found in a set of cell death-inducing proteins. Electron microscopy structural studies revealed that LACTB can polymerize to form stable filaments with lengths ranging from twenty to several hundred nanometers. These data suggest that LACTB filaments define a distinct microdomain in the intermembrane space. A possible role of LACTB filaments is proposed in the intramitochondrial membrane organization and microcompartmentation. The implications of these findings offer novel insight into the evolution of mitochondria. Further studies of the LACTB function might provide a tool to treat mitochondria-related metabolic diseases.

Simple software processes and tests improve the reliability and usefulness of a model.

Models are abstractions of reality that have predetermined limits (often not consciously thought through) on what problem domains the models can be used to explore. These limits are determined by the range of observed data used to construct and validate the model. However, it is important to remember that operating the model beyond these limits, one of the reasons for building the model in the first place, potentially brings unwanted behaviour and thus reduces the usefulness of the model. Our experience with the Agricultural Production Systems Simulator (APSIM), a farming systems model, has led us to adapt techniques from the disciplines of modelling and software development to create a model development process. This process is simple, easy to follow, and brings a much higher level of stability to the development effort, which then delivers a much more useful model. A major part of the process relies on having a range of detailed model tests (unit, simulation, sensibility, validation) that exercise a model at various levels (sub-model, model and simulation). To underline the usefulness of testing, we examine several case studies where simulated output can be compared with simple relationships. For example, output is compared with crop water use efficiency relationships gleaned from the literature to check that the model reproduces the expected function. Similarly, another case study attempts to reproduce generalised hydrological relationships found in the literature. This paper then describes a simple model development process (using version control, automated testing and differencing tools), that will enhance the reliability and usefulness of a model.

Structural origin of set-reset processes in Ge15Te83Si2 glass investigated using in situ Raman scattering and transmission electron microscopy

We report here that the structural origin of an easily reversible Ge15Te83Si2 glass can be a promising candidate for phase change random access memories. In situ Raman scattering studies on Ge15Te83Si2 sample, undertaken during the amorphous set and reset processes, indicate that the degree of disorder in the glass is reduced from off to set state. It is also found that the local structure of the sample under reset condition is similar to that in the amorphous off state. Electron microscopic studies on switched samples indicate the formation of nanometric sized particles of c-SiTe2 structure. ©2009 American Institute of Physics

Occurrence, divergence and evolution of intrinsic terminators across Eubacteria

In Escherichia coli, the canonical intrinsic terminator of transcription includes a palindrome followed by a U-trail on the transcript. The apparent underrepresentation of such terminators in eubacterial genomes led us to develop a rapid and accurate algorithm, GeSTer, to predict putative intrinsic terminators. Now, we have analyzed 378 genome sequences with an improved version of GeSTer. Our results indicate that the canonical E. coli type terminators are not overwhelmingly abundant in eubacteria. The atypical structures, having stem-loop structures but lacking ‘U’ trail, occur downstream of genes in all the analyzed genomes but different phyla show conserved preference for different types of terminators. This propensity correlates with genomic GC content and presence of the factor, Rho. 60–70% of identified terminators in all the genomes show “optimized” stem-length and ΔG. These results provide evidence that eubacteria extensively rely on the mechanism of intrinsic termination, with a considerable divergence in their structure, positioning and prevalence. The software and detailed results for individual genomes are freely available on request

An Online Implementable Differential Evolution Tuned Optimal Guidance Law

This paper proposes a novel application of differential evolution to solve a difficult dynamic optimisation or optimal control problem. The miss distance in a missile-target engagement is minimised using differential evolution. The difficulty of solving it by existing conventional techniques in optimal control theory is caused by the nonlinearity of the dynamic constraint equation, inequality constraint on the control input and inequality constraint on another parameter that enters problem indirectly. The optimal control problem of finding the minimum miss distance has an analytical solution subject to several simplifying assumptions. In the approach proposed in this paper, the initial population is generated around the seed value given by this analytical solution. Thereafter, the algorithm progresses to an acceptable final solution within a few generations, satisfying the constraints at every iteration. Since this solution or the control input has to be obtained in real time to be of any use in practice, the feasibility of online implementation is also illustrated.

Microstructural evolution and grain boundary sliding in a superplastic magnesium AZ31 alloy

Tensile experiments on a fine-grained single-phase Mg–Zn–Al alloy (AZ31) at 673 K revealed superplastic behavior with an elongation to failure of 475% at 1 × 10−4 s−1 and non-superplastic behavior with an elongation to failure of 160% at 1 × 10−2 s−1; the corresponding strain rate sensitivities under these conditions were 0.5 and 0.2, respectively. Measurements indicated that the grain boundary sliding (GBS) contribution to strain ξ was 30% under non-superplastic conditions; there was also a significant sharpening in texture during such deformation. Under superplastic conditions, ξ was 50% at both low and high elongations of 20% and 120%; the initial texture became more random under such conditions. In non-superplastic conditions, deformation occurred under steady-state conditions without grain growth before significant flow localization whereas, under superplastic conditions, there was grain growth during the early stages of deformation, leading to strain hardening. The grains retained equiaxed shapes under all experimental conditions. Superplastic deformation is attributed to GBS, while non-superplastic deformation is attributed to intragranular dislocation creep with some contribution from GBS. The retention of equiaxed grain shapes during dislocation creep is consistent with a model based on local recovery related to the disturbance of triple junctions.

Transient pulse evolution of active-mode locking in an intra-cavity frequency-doubled laser

The theory of transient mode locking for an active modulator in an intracavity frequency-doubled laser is presented. The theory is applied to mode-locked and intracavity frequency-doubled Nd:YAG laser and the mode-locked pulse width is plotted as a function of number of round trips inside the cavity. It is found that the pulse compression is faster and the system takes a very short time to approach the steady state in the presence of a second harmonic generating crystal inside the laser cavity. The effect of modulation depth and the second harmonic conversion efficiency on the temporal behavior of the pulse width is discussed.

Stanniocalcin-1 in cell stress and differentiation

Stanniocalcin-1 (STC-1) is a 56 kD homodimeric protein which was originally identified in bony fish, where it regulates calcium/phosphate homeostasis and protects against toxic hypercalcemia. STC-1 was considered unique to fish until the cloning of cDNA for human STC-1 in 1995 and mouse Stc-1 in 1996. STC-1 is conserved through evolution with human and salmon STC-1 sharing 60% identity and 80% similarity. The surprisingly high homology between mammalian and fish STC-1 and the protective actions of STC-1 in terminally differentiated neurons, originally reported by my colleagues, prompted me to further study the role of STC-1 in cell stress and differentiation. One purpose was to determine whether there is an inter-relationship between terminally differentiated cells and STC-1 expression. The study revealed an accumulation of STC-1 in mature megakaryocytes and adipocytes, i.e. postmitotic cells with limited or lost proliferative capacity. Still proliferating uninduced cells were negative for STC-1 mRNA and protein, whereas differentiating cells accumulated STC-1 in their cytoplasm. Interestingly, in liposarcomas the grade inversely correlated with STC-1 expression. Another aim was to study how STC-1 gene expression is regulated. Given that IL-6 is a cytokine with neuroprotective actions, by unknown mechanisms, we examined whether IL-6 regulates STC-1 gene expression. Treatment of human neural Paju cells with IL-6 induced a dose-dependent upregulation of STC-1 mRNA levels. This induction of STC-1 expression by IL-6 occurred mainly through the MAPK signaling pathway. Furthermore, I studied the role of IL-6-mediated STC-1 expression as a mechanism of cytoprotection conferred by hypoxic preconditioning (HOPC) in brain and heart. My findings show that Stc-1 was upregulated in brain after hypoxia treatment. In the brain of IL-6 deficient mice, however, no upregulation of Stc-1 expression was evident. After induced brain injury the STC-1 response in brains of IL-6 transgenic mice, with IL-6 overexpression in astroglial cells, was stronger than in brains of WT mice. These results indicate that IL-6-mediated expression of STC-1 is one molecular mechanism of HOPC-induced tolerance to brain ischemia. The protection conferred by HOPC in heart occurs during a bimodal time course comprising early and delayed preconditioning. Interestingly, my results showed that the expression of Stc-1 in heart was upregulated in a biphasic manner during HOPC. IL-6 deficient mice did not, however, show a similar biphasic manner of Stc-1 upregulation as did WT mice. Instead, only an early upregulation of Stc-1 expression was evident. The results suggest that the upregulation of Stc-1 during the delayed preconditioning is IL-6-dependent. The upregulated expression of Stc-1 during the early preconditioning, however, is only partly IL-6-dependent and possibly also directly mediated by HIF-1. These findings suggest that STC-1 is a pro-survival protein for terminally differentiated cells and that STC-1 expression may in fact be regulated by stress. In addition, I show that STC-1 gene upregulation, mediated in part by IL-6, is a new mechanism of protection conferred by HOPC in brain and heart. Because of its importance for fundamental biological processes, such as differentiation and cytoprotection, STC-1 may have therapeutic implications for management of stroke, neurodegenerative diseases, cancer, and obesity.

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.