Simulating the effects of dominance and epistasis on selection response in the CIMMYT wheat breeding program using QuCim

Plant breeders use many different breeding methods to develop superior cultivars. However, it is difficult, cumbersome, and expensive to evaluate the performance of a breeding method or to compare the efficiencies of different breeding methods within an ongoing breeding program. To facilitate comparisons, we developed a QU-GENE module called QuCim that can simulate a large number of breeding strategies for self-pollinated species. The wheat breeding strategy Selected Bulk used by CIMMYT's wheat breeding program was defined in QuCim as an example of how this is done. This selection method was simulated in QuCim to investigate the effects of deviations from the additive genetic model, in the form of dominance and epistasis, on selection outcomes. The simulation results indicate that the partial dominance model does not greatly influence genetic advance compared with the pure additive model. Genetic advance in genetic systems with overdominance and epistasis are slower than when gene effects are purely additive or partially dominant. The additive gene effect is an appropriate indicator of the change in gene frequency following selection when epistasis is absent. In the absence of epistasis, the additive variance decreases rapidly with selection. However, after several cycles of selection it remains relatively fixed when epistasis is present. The variance from partial dominance is relatively small and therefore hard to detect by the covariance among half sibs and the covariance among full sibs. The dominance variance from the overdominance model can be identified successfully, but it does not change significantly, which confirms that overdominance cannot be utilized by an inbred breeding program. QuCim is an effective tool to compare selection strategies and to validate some theories in quantitative genetics.

The specification of vector autoregressive moving average models

In this paper we propose a new identification method based on the residual white noise autoregressive criterion (Pukkila et al. , 1990) to select the order of VARMA structures. Results from extensive simulation experiments based on different model structures with varying number of observations and number of component series are used to demonstrate the performance of this new procedure. We also use economic and business data to compare the model structures selected by this order selection method with those identified in other published studies.

Use of drought response index for identification of drought tolerant genotypes in rainfed lowland rice

Drought is a major constraint for rice production in the rainfed lowlands in Southeast Asia and Eastern India. The breeding programs for tainted lowland rice in these regions focus on adaptation to a range of drought conditions. However, a method of selection of drought tolerant genotypes has not been established and is considered to be one of the constraints faced by rice breeders. Drought response index (DRI) is based on grain yield adjusted for variation in potential yield and flowering date, and has been used recently, but its consistency among drought environments and hence its usefulness is not certain. In order to establish a selection method and subsequently to identify donor parents for drought resistance breeding, a series of experiments with 15 contrasting genotypes was conducted under well-watered and managed drought conditions at two sites for 5 years in Cambodia. Water level in the field was recorded and used to estimate the relative water level (WLREL) around flowering as an index of the severity of water deficit at the time of flowering for each entry. This was used to determine if DRI or yield reduction was due to drought tolerance or related to the amount of available water at flowering, i.e. drought escape. Grain yield reduction due to drought ranged from 12 to 46%. The drought occurred mainly during the reproductive phase, while four experiments had water stress from the early vegetative stage. There was significant variation for water availability around flowering among the nine experiments and this was associated with variation in mean yield reduction. Genotypic variation in DRI was consistent among most experiments, and genotypic mean DRI ranged from -0.54 to 0.47 (LSD 5% = 0.47). Genotypic variation in DRI was not related to WLREL around flowering in the nine environments. It is concluded that selection for DRI under drought conditions would allow breeders to identify donor lines with high drought tolerance as an important component of breeding better adapted varieties for the rainfed lowlands; two genotypes were identified with high DRI and low yield reduction and were subsequently used in the breeding program in Cambodia. (c) 2006 Elsevier B.V. All rights reserved.

Predictors of osteoarthritis(OA) improvement in patients treated with Hylan G-F for knee OA.

Aim of study: Different criteria for treatment response were explored to identify predictors of OA improvement. Analyses were based on data from a previously reported 1-year randomized controlled trial of appropriate care with or without hylan G-F 20 in patients with knee OA. Methods: Five definitions of ‘‘patient responder’’ from baseline to month 12 were examined: at least 20% reduction in WOMAC pain score; at least 20% reduction in WOMAC pain score and at least 20% reduction in either the WOMAC stiffness or function score; OARSI responder criteria (Propositions A and B) for intra-articular drugs; and OMERACT-OARSI responder criteria (Proposition D). As an a posteriori analysis, multivariable logistic regression models for each definition of patient responder were developed using a forward selection method. The following variables were defined prior to modeling and considered in the model along with two-way interactions: age (O65 years), BMI, gender, X-ray grade (0, I, II vs III, IV), co-morbidity (1 or 2 conditions vs 3 or more), duration of OA in study knee (years), previous surgery of study knee, hylan G-F 20 injection technique, WOMAC pain, stiffness and function, and treatment group. Results: Hylan G-F 20 was a predictor of improvement for all patient responder definitions P ! 0.001; odds of improvement were 2.7 or higher for patients in the hylan G-F 20 group compared to appropriate care without hylan G-F 20. For three of the five patient responder definitions, X-ray grade was a predictor of improvement (P ! 0.10; lower X-ray grade increased the odds of improvement). For four of the five patient responder definitions, duration of OA was a predictor of improvement (P ! 0.10; shorter duration of OA increased the odds of improvement). Conclusion: Analyses showed that appropriate care with hylan G-F 20 is the dominant predictor of patient improvement. While high grade structural damage does not preclude a response, patients who are targeted early in the disease process when less structural damage has occurred, may have a greater chance of improvement.

Is the method of signal analysis and test selection important for measuring standing balance in subjects with persistent whiplash?

Dizziness and or unsteadiness, associated with episodes of loss of balance, are frequent complaints in those suffering from persistent problems following a whiplash injury. Research has been inconclusive with respect to possible aetiology, discriminative tests and analyses used. The aim of this pilot research was to identify the test conditions and the most appropriate method for the analysis of sway that may differentiate subjects with persistent whiplash associated disorders (WAD) from healthy controls. The six conditions of the Clinical Test for Sensory Interaction in Balance was performed in both comfortable and tandem stance in 20 subjects with persistent WAD compared to 20 control subjects. The analyses were carried out using a traditional method of measurement, total sway distance, to results obtained from the use of wavelet analysis. Subjects with WAD were significantly less able to complete the tandem stance tests on a firm surface than controls. In comfortable stance, using wavelet analysis, significant differences between subjects with WAD and the control group were evident in total energy of the trace for all test conditions apart from eyes open on the firm surface. In contrast, the results of the analysis using total sway distance revealed no significant differences between groups across all six conditions. Wavelet analysis may be more appropriate for detecting disturbances in balance in whiplash subjects because the technique allows separation of the noise from the underlying systematic effect of sway. These findings will be used to direct future studies on the aeitiology of balance disturbances in WAD. (c) 2004 Elsevier B.V. All rights reserved.

Method for generation of homogeneous multicellular tumor spheroids applicable to a wide variety of cell types

Multicellular tumor spheroids (MCTS) are used as organotypic models of normal and solid tumor tissue. Traditional techniques for generating MCTS, such as growth on nonadherent surfaces, in suspension, or on scaffolds, have a number of drawbacks, including the need for manual selection to achieve a homogeneous population and the use of nonphysiological matrix compounds. In this study we describe a mild method for the generation of MCTS, in which individual spheroids form in hanging drops suspended from a microtiter plate. The method has been successfully applied to a broad range of cell lines and shows nearly 100% efficiency (i.e., one spheroid per drop). Using the hepatoma cell line, HepG2, the hanging drop method generated well-rounded MCTS with a narrow size distribution (coefficient of variation [CV] 10% to 15%, compared with 40% to 60% for growth on nonadherent surfaces). Structural analysis of HepG2 and a mammary gland adenocarcinoma cell line, MCF-7, composed spheroids, revealed highly organized, three-dimensional, tissue-like structures with an extensive extracellular matrix. The hanging drop method represents an attractive alternative for MCTS production, because it is mild, can be applied to a wide variety of cell lines, and can produce spheroids of a homogeneous size without the need for sieving or manual selection. The method has applications for basic studies of physiology and metabolism, tumor biology, toxicology, cellular organization, and the development of bioartificial tissue. (C) 2003 Wiley Periodicals, Inc.

Guideline to reference gene selection for quantitative real-time PCR

Today, quantitative real-time PCR is the method of choice for rapid and reliable quantification of mRNA transcription. However, for an exact comparison of mRNA transcription in different samples or tissues it is crucial to choose the appropriate reference gene. Recently glyceraldehyde 3-phosphate dehydrogenase and P-actin have been used for that purpose. However, it has been reported that these genes as well as alternatives, like rRNA genes, are unsuitable references, because their transcription is significantly regulated in various experimental settings and variable in different tissues. Therefore, quantitative real-time PCR was used to determine the mRNA transcription profiles of 13 putative reference genes, comparing their transcription in 16 different tissues and in CCRF-HSB-2 cells stimulated with 12-O-tetradecanoylphorbol-13-acetate and ionomycin. Our results show that Classical reference genes are indeed unsuitable, whereas the RNA polymerase II gene was the gene with the most constant expression in different tissues and following stimulation in CCRF-HSB-2 cells. (C) 2003 Elsevier Inc. All rights reserved.

The selection of a model microalgal species as biomaterial for a novel aquatic phytotoxicity assay

A phytotoxicity assay based on the ToxY-PAM dual-channel yield analyser has been developed and successfully incorporated into field assessments for the detection of phytotoxicants in water. As a means of further exploring the scope of the assay application and of selecting a model biomaterial to complement the instrument design, nine algal species were exposed to four chemical substances deemed of priority for water quality monitoring purposes (chlorpyrifos, copper, diuron and nonylphenol ethoxylate). Inter-species differences in sensitivity to the four toxicants varied by a factor of 1.9-100. Measurements of photosystem-II quantum yield using these nine single-celled microalgae as biomaterial corroborated previous studies which have shown that the ToxY-PAM dual-channel yield analyser is a highly sensitive method for the detection of PS-II impacting herbicides. Besides Phaeodactylum tricornutum, the previously applied biomaterial, three other species consistently performed well (Nitzschia closterium, Chlorella vulgaris and Dunaliella tertiolecta) and will be used in further test optimisation experiments. In addition to sensitivity, response time was evaluated and revealed a high degree of variation between species and toxicants. While most species displayed relatively weak and slow responses to copper, C. vulgaris demonstrated an IC10 of 51 μ g L-1, with maximum response measured within 25 minutes and inhibition being accompanied by a large decrease in fluorescence yield. The potential for this C vulgaris-based bioassay to be used for the detection of copper is discussed. There was no evidence that the standard ToxY-PAM protocol, using these unicellular algae species, could be used for the detection of chlorpyrifos or nonylphenol ethoxylate at environmentally relevant levels. © 2005 Elsevier B.V. All rights reserved.

Composite auction method for suppressing unreasonable electricity price spikes in a competitive electricity market

The existence of undesirable electricity price spikes in a competitive electricity market requires an efficient auction mechanism. However, many of the existing auction mechanism have difficulties in suppressing such unreasonable price spikes effectively. A new auction mechanism is proposed to suppress effectively unreasonable price spikes in a competitive electricity market. It optimally combines system marginal price auction and pay as bid auction mechanisms. A threshold value is determined to activate the switching between the marginal price auction and the proposed composite auction. Basically when the system marginal price is higher than the threshold value, the composite auction for high price electricity market is activated. The winning electricity sellers will sell their electricity at the system marginal price or their own bid prices, depending on their rights of being paid at the system marginal price and their offers' impact on suppressing undesirable price spikes. Such economic stimuli discourage sellers from practising economic and physical withholdings. Multiple price caps are proposed to regulate strong market power. We also compare other auction mechanisms to highlight the characteristics of the proposed one. Numerical simulation using the proposed auction mechanism is given to illustrate the procedure of this new auction mechanism.

A new method for conservation planning for the persistence of multiple species

Although the aim of conservation planning is the persistence of biodiversity, current methods trade-off ecological realism at a species level in favour of including multiple species and landscape features. For conservation planning to be relevant, the impact of landscape configuration on population processes and the viability of species needs to be considered. We present a novel method for selecting reserve systems that maximize persistence across multiple species, subject to a conservation budget. We use a spatially explicit metapopulation model to estimate extinction risk, a function of the ecology of the species and the amount, quality and configuration of habitat. We compare our new method with more traditional, area-based reserve selection methods, using a ten-species case study, and find that the expected loss of species is reduced 20-fold. Unlike previous methods, we avoid designating arbitrary weightings between reserve size and configuration; rather, our method is based on population processes and is grounded in ecological theory.

Numerical method for determination of the NMR frequency of the single-qubit operation in a silicon-based solid-state quantum computer

A numerical method is introduced to determine the nuclear magnetic resonance frequency of a donor (P-31) doped inside a silicon substrate under the influence of an applied electric field. This phosphorus donor has been suggested for operation as a qubit for the realization of a solid-state scalable quantum computer. The operation of the qubit is achieved by a combination of the rotation of the phosphorus nuclear spin through a globally applied magnetic field and the selection of the phosphorus nucleus through a locally applied electric field. To realize the selection function, it is required to know the relationship between the applied electric field and the change of the nuclear magnetic resonance frequency of phosphorus. In this study, based on the wave functions obtained by the effective-mass theory, we introduce an empirical correction factor to the wave functions at the donor nucleus. Using the corrected wave functions, we formulate a first-order perturbation theory for the perturbed system under the influence of an electric field. In order to calculate the potential distributions inside the silicon and the silicon dioxide layers due to the applied electric field, we use the multilayered Green's functions and solve an integral equation by the moment method. This enables us to consider more realistic, arbitrary shape, and three-dimensional qubit structures. With the calculation of the potential distributions, we have investigated the effects of the thicknesses of silicon and silicon dioxide layers, the relative position of the donor, and the applied electric field on the nuclear magnetic resonance frequency of the donor.

Conservation planning with irreplaceability: does the method matter?

A number of systematic conservation planning tools are available to aid in making land use decisions. Given the increasing worldwide use and application of reserve design tools, including measures of site irreplaceability, it is essential that methodological differences and their potential effect on conservation planning outcomes are understood. We compared the irreplaceability of sites for protecting ecosystems within the Brigalow Belt Bioregion, Queensland, Australia, using two alternative reserve system design tools, Marxan and C-Plan. We set Marxan to generate multiple reserve systems that met targets with minimal area; the first scenario ignored spatial objectives, while the second selected compact groups of areas. Marxan calculates the irreplaceability of each site as the proportion of solutions in which it occurs for each of these set scenarios. In contrast, C-Plan uses a statistical estimate of irreplaceability as the likelihood that each site is needed in all combinations of sites that satisfy the targets. We found that sites containing rare ecosystems are almost always irreplaceable regardless of the method. Importantly, Marxan and C-Plan gave similar outcomes when spatial objectives were ignored. Marxan with a compactness objective defined twice as much area as irreplaceable, including many sites with relatively common ecosystems. However, targets for all ecosystems were met using a similar amount of area in C-Plan and Marxan, even with compactness. The importance of differences in the outcomes of using the two methods will depend on the question being addressed; in general, the use of two or more complementary tools is beneficial.

A bayesian artificial neural network method to characterise laminar defects using dynamic measurements

This paper reports on the development of an artificial neural network (ANN) method to detect laminar defects following the pattern matching approach utilizing dynamic measurement. Although structural health monitoring (SHM) using ANN has attracted much attention in the last decade, the problem of how to select the optimal class of ANN models has not been investigated in great depth. It turns out that the lack of a rigorous ANN design methodology is one of the main reasons for the delay in the successful application of the promising technique in SHM. In this paper, a Bayesian method is applied in the selection of the optimal class of ANN models for a given set of input/target training data. The ANN design method is demonstrated for the case of the detection and characterisation of laminar defects in carbon fibre-reinforced beams using flexural vibration data for beams with and without non-symmetric delamination damage.

Australia: new screening method for cold tolerance during the reproducitve stage in rice

Low temperature, particularly during the reproductive stage of the development of rice, limits productivity in the Riverina region of New South Wales (NSW). This study primarily examined genotypic differences in cold damage that are associated with low temperature during reproductive development. Results from experiments in temperature-controlled rooms and the cold water facility were combined with four years of field experiments, which used natural exposure to low temperature to examine the response of over 50 cultivars from diverse origins. Plants were exposed to day/night air temperatures of 27°/13°C in temperature-controlled rooms and to a constant temperature of 19°C in the cold water facility. Low temperature treatments were imposed from panicle initiation (PI) to 50% heading. In field experiments several techniques were used to increase the likelihood of inducing cold damage such as sequential sowing dates (five to eight sowing dates each year), shallow water depths (5cm) and high nitrogen rates (e.g. 300kgN ha-1). Several cultivars were identified that were more cold tolerant than Australia’s commercial cultivars.