An investigation into effects of long-distance seed dispersal on organelle population genetic structure and colonization rate: a model analysis

A simulation-based modelling approach is used to examine the effects of stratified seed dispersal (representing the distribution of the majority of dispersal around the maternal parent and also rare long-distance dispersal) on the genetic structure of maternally inherited genomes and the colonization rate of expanding plant populations. The model is parameterized to approximate postglacial oak colonization in the UK, but is relevant to plant populations that exhibit stratified seed dispersal. The modelling approach considers the colonization of individual plants over a large area (three 500 km x 10 km rolled transects are used to approximate a 500 km x 300 km area). Our approach shows how the interaction of plant population dynamics with stratified dispersal can result in a spatially patchy haplotype structure. We show that while both colonization speeds and the resulting genetic structure are influenced by the characteristics of the dispersal kernel, they are robust to changes in the periodicity of long-distance events, provided the average number of long-distance dispersal events remains constant. We also consider the effects of additional physical and environmental mechanisms on plant colonization. Results show significant changes in genetic structure when the initial colonization of different haplotypes is staggered over time and when a barrier to colonization is introduced. Environmental influences on survivorship and fecundity affect both the genetic structure and the speed of colonization. The importance of these mechanisms in relation to the postglacial spread and genetic structure of oak in the UK is discussed.

Immunohistological analysis of Tannerella forsythia-induced lesions in a murine model

Tannerella forsythia has been implicated as a defined periodontal pathogen. In the present study a mouse model was used to determine the phenotype of leukocytes in the lesions induced by subcutaneous injections of either live (group A) or nonviable (group B) T. forsythia. Control mice (group C) received the vehicle only. Lesions were excised at days 1, 2, 4, and 7. An avidin-biotin immunoperoxidase method was used to stain infiltrating CD4(+) and CD8(+) T cells, CD14(+) macrophages, CD19(+) B cells, and neutrophils. Hematoxylin and eosin sections demonstrated lesions with central necrotic cores surrounded by neutrophils, macrophages and lymphocytes in both group A and group B mice. Lesions from control mice exhibited no or only occasional solitary leukocytes. In both groups A and B, neutrophils were the dominant leukocyte in the lesion 1 day after injection, the numbers decreasing over the 7-day experimental period. There was a relatively low mean percent of CD4(+) and CD8(+) T cells in the lesions and, whereas the percent of CD8(+) T cells remained constant, there was a significant increase in the percent of CD4(+) T cells at day 7. This increase was more evident in group A mice. The mean percent of CD14(+) macrophages and CD19(+) B cells remained low over the experimental period, although there was a significantly higher mean percent of CD19(+) B cells at day 1. In conclusion, the results showed that immunization of mice with live T. forsythia induced a stronger immune response than nonviable organisms. The inflammatory response presented as a nonspecific immune response with evidence of an adaptive (T-cell) response by day 7. Unlike Porphyromonas gingivalis, there was no inhibition of neutrophil migration.

An experimental and finite element poroelastic creep response analysis of an intervertebral hydrogel disc model in axial compression

A hydrogel intervertebral disc (lVD) model consisting of an inner nucleus core and an outer anulus ring was manufactured from 30 and 35% by weight Poly(vinyl alcohol) hydrogel (PVA-H) concentrations and subjected to axial compression in between saturated porous endplates at 200 N for 11 h, 30 min. Repeat experiments (n = 4) on different samples (N = 2) show good reproducibility of fluid loss and axial deformation. An axisymmetric nonlinear poroelastic finite element model with variable permeability was developed using commercial finite element software to compare axial deformation and predicted fluid loss with experimental data. The FE predictions indicate differential fluid loss similar to that of biological IVDs, with the nucleus losing more water than the anulus, and there is overall good agreement between experimental and finite element predicted fluid loss. The stress distribution pattern indicates important similarities with the biological lVD that includes stress transference from the nucleus to the anulus upon sustained loading and renders it suitable as a model that can be used in future studies to better understand the role of fluid and stress in biological IVDs. (C) 2005 Springer Science + Business Media, Inc.

Functional-structural plant modelling using a combination of architectural analysis, L-systems and a canonical model of function

This paper presents a new method for producing a functional-structural plant model that simulates response to different growth conditions, yet does not require detailed knowledge of underlying physiology. The example used to present this method is the modelling of the mountain birch tree. This new functional-structural modelling approach is based on linking an L-system representation of the dynamic structure of the plant with a canonical mathematical model of plant function. Growth indicated by the canonical model is allocated to the structural model according to probabilistic growth rules, such as rules for the placement and length of new shoots, which were derived from an analysis of architectural data. The main advantage of the approach is that it is relatively simple compared to the prevalent process-based functional-structural plant models and does not require a detailed understanding of underlying physiological processes, yet it is able to capture important aspects of plant function and adaptability, unlike simple empirical models. This approach, combining canonical modelling, architectural analysis and L-systems, thus fills the important role of providing an intermediate level of abstraction between the two extremes of deeply mechanistic process-based modelling and purely empirical modelling. We also investigated the relative importance of various aspects of this integrated modelling approach by analysing the sensitivity of the standard birch model to a number of variations in its parameters, functions and algorithms. The results show that using light as the sole factor determining the structural location of new growth gives satisfactory results. Including the influence of additional regulating factors made little difference to global characteristics of the emergent architecture. Changing the form of the probability functions and using alternative methods for choosing the sites of new growth also had little effect. (c) 2004 Elsevier B.V. All rights reserved.

Wave-generated flow on coral reefs - an analysis for two-dimensional horizontal reef-tops with steep faces

Waves breaking on the seaward rim of a coral reef generate a flow of water from the exposed side of the reef to the sheltered side and/or to either channels through the reef-rim or lower sections of the latter. This wave-generated flow is driven by the water surface gradient resulting from the wave set-up created by the breaking waves. This paper reviews previous approaches to modelling wave-generated flows across coral reefs and discusses the influence of reef morphology and roughness upon these flows. Laboratory measurements upon a two-dimensional horizontal reef platform with a steep reef face provide the basis for extending a previous theoretical analysis for wave set-up on a reef in the absence of a flow [Gourlay, M.R., 1996b. Wave set-up on coral reefs. 2. Set-up on reefs with various profiles. Coastal Engineering 28, 1755] to include the interaction between a unidirectional flow and the wave set-up. The laboratory model results are then used to demonstrate that there are two basic reef-top flow regimes-reef-top control and reef-rim control. Using open channel flow theory, analytical relationships are derived for the reef-top current velocity in terms of the offreef wave conditions, the reef-top water depth and the physical characteristics of the reef-top topography. The wave set-up and wave-generated flow relationships are found to predict experimental values with reasonable accuracy in most cases. The analytical relationships are used to investigate wave-generated flows into a boat harbour channel on Heron Reef in the southern Great Barrier Reef. (c) 2005 Elsevier B.V. All rights reserved.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

Leading coordinate analysis of reaction pathways in proton chain transfer: Application to a two-proton transfer model for the green fluorescent protein

The ‘leading coordinate’ approach to computing an approximate reaction pathway, with subsequent determination of the true minimum energy profile, is applied to a two-proton chain transfer model based on the chromophore and its surrounding moieties within the green fluorescent protein (GFP). Using an ab initio quantum chemical method, a number of different relaxed energy profiles are found for several plausible guesses at leading coordinates. The results obtained for different trial leading coordinates are rationalized through the calculation of a two-dimensional relaxed potential energy surface (PES) for the system. Analysis of the 2-D relaxed PES reveals that two of the trial pathways are entirely spurious, while two others contain useful information and can be used to furnish starting points for successful saddle-point searches. Implications for selection of trial leading coordinates in this class of proton chain transfer reactions are discussed, and a simple diagnostic function is proposed for revealing whether or not a relaxed pathway based on a trial leading coordinate is likely to furnish useful information.

Investigating the performance analysis of EASI algorithm and EKENS algorithm in nonlinear model

This paper investigates the performance of EASI algorithm and the proposed EKENS algorithm for linear and nonlinear mixtures. The proposed EKENS algorithm is based on the modified equivariant algorithm and kernel density estimation. Theory and characteristic of both the algorithms are discussed for blind source separation model. The separation structure of nonlinear mixtures is based on a nonlinear stage followed by a linear stage. Simulations with artificial and natural data demonstrate the feasibility and good performance of the proposed EKENS algorithm.

Motion analysis of right ventricular wall based on an electromechanical biventricular model

Based on our previously developed electrical heart model, an electromechanical biventricular model, which couples the electrical property and mechanical property of the heart, was constructed and the right ventricular wall motion and deformation was simulated using this model. The model was developed on the basis of composite material theory and finite element method. The excitation propagation was simulated by electrical heart model, and the resultant active forces were used to study the ventricular wall motion during systole. The simulation results show that: (1) The right ventricular free wall moves towards the septum, and at the same time, the base and middle of free wall move towards the apex, which reduce the volume of right ventricle; (2) The minimum principle strain (E3) is largest at the apex, then at the middle of free wall, and its direction is in the approximate direction of epicardial muscle fibers. These results are in good accordance with solutions obtained from MR tagging images. It suggests that such electromechanical biventricular model can be used to assess the mechanical function of two ventricles.

Formal analysis of human-computer interaction using model checking

Experiments with simulators allow psychologists to better understand the causes of human errors and build models of cognitive processes to be used in human reliability assessment (HRA). This paper investigates an approach to task failure analysis based on patterns of behaviour, by contrast to more traditional event-based approaches. It considers, as a case study, a formal model of an air traffic control (ATC) system which incorporates controller behaviour. The cognitive model is formalised in the CSP process algebra. Patterns of behaviour are expressed as temporal logic properties. Then a model-checking technique is used to verify whether the decomposition of the operator's behaviour into patterns is sound and complete with respect to the cognitive model. The decomposition is shown to be incomplete and a new behavioural pattern is identified, which appears to have been overlooked in the analysis of the data provided by the experiments with the simulator. This illustrates how formal analysis of operator models can yield fresh insights into how failures may arise in interactive systems.