The E(NK) model: extending the NK model to incorporate gene-by-environment interactions and epistasis for diploid genomes

The haploid NK model developed by Kauffman can be extended to diploid genomes and to incorporate gene-by-environment interaction effects in combination with epistasis. To provide the flexibility to include a wide range of forms of gene-by-environment interactions, a target population of environment types (TPE) is defined. The TPE consists of a set of E different environment types, each with their own frequency of occurrence. Each environment type conditions a different NK gene network structure or series of gene effects for a given network structure, providing the framework for defining gene-by-environment interactions. Thus, different NK models can be partially or completely nested within the E environment types of a TPE, giving rise to the E(NK) model for a biological system. With this model it is possible to examine how populations of genotypes evolve in context with properties of the environment that influence the contributions of genes to the fitness values of genotypes. We are using the E(NK) model to investigate how both epistasis and gene-by-environment interactions influence the genetic improvement of quantitative traits by plant breeding strategies applied to agricultural systems. © 2002 Wiley Periodicals, Inc.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

Control of Aedes vectors of dengue in three provinces of Vietnam by use of Mesocyclops (Copepoda) and community-based methods validated by entomologic, clinical, and serological surveillance

We describe remarkable success in controlling dengue vectors, Aedes aegypti (L.) and Aedes albopictus (Skuse), in 6 communes with 11,675 households and 49,647 people in the northern provinces of Haiphong, Hung Yen, and Nam Dinh in Vietnam. The communes were selected for high-frequency use of large outdoor concrete tanks and wells. These were found to be the source of 49.6-98.4% of Ae. aegypti larvae, which were amenable to treatment with local Mesocyclops, mainly M. woutersi Van der Velde, M. aspericornis (Daday) and M. thermocyclopoides Harada. Knowledge, attitude, and practice surveys were performed to determine whether the communities viewed dengue and dengue hemorrhagic fever as a serious health threat; to determine their knowledge of the etiology, attitudes, and practices regarding control methods including Mesocyclops; and to determine their receptivity to various information methods. On the basis of the knowledge, attitude, and practice data, the community-based dengue control program comprised a system of local leaders, health volunteer teachers, and schoolchildren, supported by health professionals. Recycling of discards for economic gain was enhanced, where appropriate, and this, plus 37 clean-up campaigns, removed small containers unsuitable for Mesocyclops treatment. A previously successful eradication at Phan Boi village (Hung Yen province) was extended to 7 other villages forming Di Su commune (1,750 households) in the current study. Complete control was also achieved in Nghia Hiep (Hung Yen province) and in Xuan Phong (Nam Dinh province); control efficacy was greater than or equal to 99.7% in the other 3 communes (Lac Vien in Haiphong, Nghia Dong, and Xuan Kien in Nam Dinh). Although tanks and wells were the key container types of Ae. aegypti productivity, discarded materials were the source of 51% of the standing crop of Ae. albopictus. Aedes albopictus larvae were eliminated from the 3 Nam Dinh communes, and 86-98% control was achieved in the other 3 communes. Variable dengue attack rates made the clinical and serological comparison of control and untreated communes problematic, but these data indicate that clinical surveillance by itself is inadequate to monitor dengue transmission.

Present and future pharacotherapy for heart failure

The pharmacotherapy currently recommended by the American College of Cardiology and the American Heart Association for heart failure (HF) is a diuretic, an angiotensin-converting enzyme inhibitor (ACEI), a β-adrenoceptor antagonist and (usually) digitalis. This current treatment of HF may be improved by optimising the dose of ACEI used, as increasing the dose of lisinopril increases its benefits in HF. Selective angiotensin receptor-1 (AT1) antagonists are effective alternatives for those who cannot tolerate ACEIs. AT1 antagonists may also be used in combination with ACEIs, as some studies have shown cumulative benefits for the combination. In addition to being used in Stage IV HF patients, in whom it has a marked benefit, spironolactone should be studied in less severe HF and in the presence of β-blockers. The use of carvedilol, extended-release metoprolol and bisoprolol should be extended to severe HF patients as these agents have been shown to decrease mortality in this group. The ancillary properties of carvedilol, particularly antagonism at prejunctional β-adrenoceptors, may give it additional benefits to selective β1-adrenoceptor antagonists. Celiprolol and bucindolol are not the β-blockers of choice in HF, as they do not decrease mortality. Although digitalis does not reduce mortality, it remains the only option for a long-term positive inotropic effect, as the long-term use of the phosphodiesterase inhibitors is associated with increased mortality. The calcium sensitising drug levosimendan may be useful in the hospital treatment of decompensated HF to increase cardiac output and improve dyspnoea and fatigue. The antiarrhythmic drug amiodarone should probably be used in patients at high risk of arrhythmic or sudden death, although this treatment may soon be superseded by the more expensive implanted cardioverter defibrillators, which are probably more effective and have fewer side effects. The natriuretic peptide nesiritide has recently been introduced for the hospital treatment of decompensated HF. Novel drugs that may be beneficial in the treatment of HF include the vasopeptidase inhibitors and the selective endothelin-A receptor antagonists but these require much more investigation. However, disappointing results have been obtained in a large clinical trial of the tumour necrosis factor α antagonist etanercept, where no likelihood of a difference between placebo and etanercept was observed. Small clinical trials with recombinant growth hormone to thicken ventricles in dilated cardiomyopathy have given variable results.

Distribution of two splice variants of the glutamate transporter GLT-1 in the developing rat retina

The distributions of a carboxyl terminal splice variant of the glutamate transporter GLT-1, referred to as GLT-1B, and the carboxyl terminus of the originally described variant of GLT-1, referred to hereafter as GLT-1alpha, were examined using specific antisera. GLT-1B was present in the retina at very early developmental stages. Labelling was demonstrable at embryonic day 14, and strong labelling was evident by embryonic day 18. Such labelling was initially restricted to populations of cone photoreceptors, the processes of which extended through the entire thickness of the retina and appeared to make contact with the retinal ganglion cells. During postnatal development the GLT-1B-positive photoreceptor processes retracted to form the outer plexiform layer, and around postnatal day 7, GLT-1B-immunoreactive bipolar cells appeared. The pattern of labelling of bipolar cell processes within the inner plexiform layer changed during postnatal development. Two strata of strongly immunoreactive terminals were initially evident in the inner plexiform layer, but by adulthood these two bands were no longer evident and labelling was restricted to the somata and processes (but not synaptic terminals) of the bipolar cells, as well as the somata, processes, and terminals of cone photoreceptors. By contrast, GLT-1alpha appeared late in postnatal development and was restricted mainly to a population of amacrine cells, although transient labelling was also associated with punctate elements in the outer plexiform layer, which may represent photoreceptor terminals, (C) 2002 Wiley-Liss, Inc.

Product form approximations for highly linear loss networks with trunk reservation

Admission controls, such as trunk reservation, are often used in loss networks to optimise their performance. Since the numerical evaluation of performance measures is complex, much attention has been given to finding approximation methods. The Erlang Fixed-Point (EFP) approximation, which is based on an independent blocking assumption, has been used for networks both with and without controls. Several more elaborate approximation methods which account for dependencies in blocking behaviour have been developed for the uncontrolled setting. This paper is an exploratory investigation of extensions and synthesis of these methods to systems with controls, in particular, trunk reservation. In order to isolate the dependency factor, we restrict our attention to a highly linear network. We will compare the performance of the resulting approximations against the benchmark of the EFP approximation extended to the trunk reservation setting. By doing this, we seek to gain insight into the critical factors in constructing an effective approximation. (C) 2003 Elsevier Ltd. All rights reserved.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

Predicting cave initiation and propagation

The most widely used method for predicting the onset of continuous caving is Laubscher's caving chart. A detailed examination of this method was undertaken which concluded that it had limitations which may impact on results, particularly when dealing with stronger rock masses that are outside current experience. These limitations relate to inadequate guidelines for adjustment factors to rock mass rating (RMR), concerns about the position on the chart of critical case history data, undocumented changes to the method and an inadequate number of data points to be confident of stability boundaries. A review was undertaken on the application and reliability of a numerical method of assessing cavability. The review highlighted a number of issues, which at this stage, make numerical continuum methods problematic for predicting cavability. This is in particular reference to sensitivity to input parameters that are difficult to determine accurately and mesh dependency. An extended version of the Mathews method for open stope design was developed as an alternative method of predicting the onset of continuous caving. A number of caving case histories were collected and analyzed and a caving boundary delineated statistically on the Mathews stability graph. The definition of the caving boundary was aided by the existence of a large and wide-ranging stability database from non-caving mines. A caving rate model was extrapolated from the extended Mathews stability graph but could only be partially validated due to a lack of reliable data.

The variety of variables in automated real-time refinement

The refinement calculus is a well-established theory for deriving program code from specifications. Recent research has extended the theory to handle timing requirements, as well as functional ones, and we have developed an interactive programming tool based on these extensions. Through a number of case studies completed using the tool, this paper explains how the tool helps the programmer by supporting the many forms of variables needed in the theory. These include simple state variables as in the untimed calculus, trace variables that model the evolution of properties over time, auxiliary variables that exist only to support formal reasoning, subroutine parameters, and variables shared between parallel processes.

Numerical study of hydrogenic effective mass theory for an impurity P donor in Si in the presence of an electric field and interfaces

In this paper we examine the effects of varying several experimental parameters in the Kane quantum computer architecture: A-gate voltage, the qubit depth below the silicon oxide barrier, and the back gate depth to explore how these variables affect the electron density of the donor electron. In particular, we calculate the resonance frequency of the donor nuclei as a function of these parameters. To do this we calculated the donor electron wave function variationally using an effective-mass Hamiltonian approach, using a basis of deformed hydrogenic orbitals. This approach was then extended to include the electric-field Hamiltonian and the silicon host geometry. We found that the phosphorous donor electron wave function was very sensitive to all the experimental variables studied in our work, and thus to optimize the operation of these devices it is necessary to control all parameters varied in this paper.

Caenorhabditis elegans mutants resistant to phosphine toxicity show increased longevity and cross-resistance to the synergistic action of oxygen

Phosphine (hydrogen phosphide, PH3) is the fumigant most widely used to protect stored products from pest infestation. Despite the importance of this chemical, little is known about its mode of action. We have created three phosphine-resistant lines (pre-1, pre-7, pre-33) in the model organism C. elegans, with LC50 values 2, 5, and 9 times greater than the fully susceptible parental strain. Molecular oxygen was shown to be an extremely effective synergist with phosphine as, under hyperoxic conditions, 100% mortality was observed in wild-type nematodes exposed to 0.1 mg/l phosphine, a nonlethal concentration in air. All three mutants were resistant to the synergistic effects of oxygen in proportion to their resistance to phosphine with one mutant, pre-33, showing complete resistance to this synergism. We take the proportionality of cross-resistance between phosphine and the synergistic effect of oxygen to imply that all three mutants circumvent a mechanism of phosphine toxicity that is directly coupled to oxygen metabolism. Compared with the wild-type strain, all three mutants have an extended average life expectancy of from 12.5 to 25.3%. This is consistent with the proposed involvement of oxidative stress in both phosphine toxicity and ageing. Because the wild-type and mutant nematodes develop at the same rate, the longevity is unlikely to be caused by a clk-type reduction in oxidative metabolism, a potential alternative mechanism of phosphine resistance.

Early pregnancy factor treatment suppresses the inflammatory response and adhesion molecule expression in the spinal cord of SJL/J mice with experimental autoimmune encephalomyelitis and the delayed-type hypersensitivity reaction to trinitrochlorobenzene in normal BALB/c mice

Early pregnancy factor (EPF) is a secreted protein, present in serum during early pregnancy and essential for maintaining viability of the embryo. It is a homologue of chaperonin 10 (Cpn10) but, unlike Cpn10, it has an extracellular role. EPF has immunosuppressive and growth regulatory properties. Previously we have reported the preparation of recombinant EPF (rEPF) and shown that treatment with rEPF will suppress clinical signs of MBP-EAE in Lewis rats and PLP-EAE in SJL/J mice. In the present study, these findings have been extended to investigate possible mechanisms involved in the action of EPF. Following treatment of mice with rEPF from the day of inoculation, there were fewer infiltrating CD3+ and CD4+ cells in the parenchyma of the spinal cord during the onset of disease and after the initial episode, compared with mice treated with vehicle. Expression of the integrins LFA-1, VLA-4 and Mac-1 and of members of the immunoglobulin superfamily of adhesion molecules ICAM-1 and VCAM-1 was suppressed in the central nervous system (CNS) following rEPF treatment. The expression of PECAM-1 was not affected. To determine if rEPF suppressed T cell activation in the periphery, the delayed-type hypersensitivity (DTH) reaction of normal BALB/c mice to trinitrochlorobenzene (TNCB) following treatment with rEPF was studied. The results showed that treatment with rEPF suppressed the DTH reaction, demonstrating the ability of EPF to downregulate the cell-mediated immune response. These results indicate that suppression of immunological mechanisms by rEPF plays a major role in the reduction of clinical signs of disease in experimental autoimmune encephalomyelitis (EAE). (C) 2003 Elsevier Science B.V. All rights reserved.

Narrating the guillotine: Punishment technology as myth and symbol

The work of Michel Foucault sees modern penal technology its ann expression of power that operates through and is motivated by a dry instrumental reason. This article draws upon Durkheim and Bakhtin to advance a radically alternative approach. It is suggested that such technology is invested with sacred and profane symbolism and is understood via emotionally charged, dramatically compelling narrative frames. Tensions between official and unauthorized discourses can be understood through a center/periphery model of culture. In an extended case study of the guillotine, it is shown dial the apparatus was initially legitimated as an expression of a sacred revolutionary code. Such a discourse was subsequently destabilized by popular medical debates that raised the specter of pain after decapitation. While inconclusive, these new motifs mobilized Gothic and grotesque themes that confronted the rationalist aesthetics of the guillotine. A situation of Bakhtinian hetoroglossia eventuated. Uncertainty, the uncanny and fable entered a discursive field of increasing complexity.

Quadriceps activation in closed and in open kinetic chain exercise

Purpose: For treatment of various knee disorders, muscles are trained in open or closed kinetic chain tasks. Coordination between the heads of the quadriceps muscle is important for stability and optimal joint loading for both the tibiofemoral and the patellofemoral joint. The aim of this study was to examine whether the quadriceps femoris muscles are activated differently in open versus closed kinetic chain tasks. Methods: Ten healthy men and women (mean age 28.5 +/- 0.7) extended the knees isometrically in open and closed kinetic chain tasks in a reaction time paradigm using moderate force. Surface electromyography (EMG) recordings were made from four different parts of the quadriceps muscle. The onset and amplitude of EMG and force data were measured. Results: In closed chain knee extension, the onset of EMG activity of the four different muscle portions of the quadriceps was more simultaneous than in the open chain. In open chain, rectus femoris (RF) had the earliest EMG onset while vastus medialis obliquus was activated last (7 +/- 13 ms after RF EMG onset) and with smaller amplitude (40 +/- 30% of maximal voluntary contraction (MVC)) than in closed chain (46 +/- 43% MVC). Conclusions: Exercise in closed kinetic chain promotes more balanced initial quadriceps activation than does exercise in open kinetic chain. This may be of importance in designing training programs aimed toward control of the patellofemoral joint.

Winding patterns for actively shielded shim coils with asymmetric target-fields

A method is presented for calculating the currents and winding patterns required to design independent zonal and tesseral shim coils for magnetic resonance imaging. Both actively shielded and unshielded configurations are considered, and the region of interest can be located asymmetrically with respect to the coil's length. Streamline, target-field and Fourier-series methods are utilized. The desired target-field is specified at two cylindrical radii, on and inside a circular conducting cylinder of length 2L and radius a. The specification is over some asymmetric portion pL < z < qL of the coil's length (-1 < p < q < 1). Arbitrary functions are used in the outer sections, -L < z < pL and qL < z < L, to ensure continuity of the magnetic field across the entire length of the coil. The entire field is then periodically extended as a half-range cosine Fourier series about either end of the coil. The resultant Fourier coefficients are then substituted into the Fourier-series expressions for the internal and external magnetic fields, and current densities and stream functions on both the primary coil and shield. A contour plot of the stream function directly gives the required coil winding patterns. Spherical harmonic analysis and shielding analysis on field calculations from a ZX shim coil indicate that example designs and theory are well matched.