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.

Genetic parameters and predicted response to selection for lean growth in an Australian commercial nucleus pig herd

Genetic parameters for performance traits in a pig population were estimated using a multi-trait derivative-free REML algorithm. The 2590 total data included 922 restrictively fed male and 1668 ad libitum fed female records. Estimates of heritability (standard error in parentheses) were 0.25 (0.03), 0.15 (0.03), and 0.30 (0.05) for lifetime daily gain, test daily gain, and P2-fat depth in males, respectively; and 0.27 (0.04) and 0.38 (0.05) for average daily gain and P2-fat depth in females, respectively. The genetic correlation between P2-fat depth and test daily gain in males was -0.17 (0.06) and between P2-fat and lifetime average daily gain in females 0.44 (0.09). Genetic correlations between sexes were 0.71 (0.11) for average daily gain and -0.30 (0.10) for P2-fat depth. Genetic response per standard deviation of selection on an index combining all traits was predicted at $AU120 per sow per year. Responses in daily gain and backfat were expected to be higher when using only male selection than when using only female selection. Selection for growth rate in males will improve growth rate and carcass leanness simultaneously.

A reassessment of genetic limits to evolutionary change

An absence of genetic variance in traits under selection is perhaps the oldest explanation for a limit to evolutionary change, but has also been the most easily dismissed. We review a range of theoretical and empirical results covering single traits to more complex multivariate systems, and show that an absence of genetic variance may be more common than is currently appreciated. From a single-trait perspective, we highlight that it is becoming clear that some trait types do not display significant levels of genetic variation, and we raise the possibility that species with restricted ranges may differ qualitatively from more widespread species in levels of genetic variance in ecologically important traits. A common misconception in many life-history studies is that a lack of genetic variance in single traits, and genetic constraints as a consequence of bivariate genetic correlations, are different causes of selection limits. We detail how interpretations of bivariate patterns are unlikely to demonstrate genetic limits to selection in many cases. We advocate a multivariate definition of genetic constraints that emphasizes the presence (or otherwise) of genetic variance in the multivariate direction of selection. For multitrait systems, recent results using longer term studies of organisms, in which more is understood concerning what traits may be under selection, have indicated that selection may exhaust genetic variance, resulting in a limit to the selection response.

Advantage of single-trial models for response to selection in wheat breeding multi-environment trials

An investigation was conducted to evaluate the impact of experimental designs and spatial analyses (single-trial models) of the response to selection for grain yield in the northern grains region of Australia (Queensland and northern New South Wales). Two sets of multi-environment experiments were considered. One set, based on 33 trials conducted from 1994 to 1996, was used to represent the testing system of the wheat breeding program and is referred to as the multi-environment trial (MET). The second set, based on 47 trials conducted from 1986 to 1993, sampled a more diverse set of years and management regimes and was used to represent the target population of environments (TPE). There were 18 genotypes in common between the MET and TPE sets of trials. From indirect selection theory, the phenotypic correlation coefficient between the MET and TPE single-trial adjusted genotype means [r(p(MT))] was used to determine the effect of the single-trial model on the expected indirect response to selection for grain yield in the TPE based on selection in the MET. Five single-trial models were considered: randomised complete block (RCB), incomplete block (IB), spatial analysis (SS), spatial analysis with a measurement error (SSM) and a combination of spatial analysis and experimental design information to identify the preferred (PF) model. Bootstrap-resampling methodology was used to construct multiple MET data sets, ranging in size from 2 to 20 environments per MET sample. The size and environmental composition of the MET and the single-trial model influenced the r(p(MT)). On average, the PF model resulted in a higher r(p(MT)) than the IB, SS and SSM models, which were in turn superior to the RCB model for MET sizes based on fewer than ten environments. For METs based on ten or more environments, the r(p(MT)) was similar for all single-trial models.

Yield response of rice (Oryza sativa L.) genotypes to drought under rainfed lowlands 4. Vegetative stage screening in the dry season

Screening for drought resistance of rainfed lowland rice using drought score (leaf death) as a selection index has a long history of use in breeding programs. Genotypic variation for drought score during the vegetative stage in two dry season screens was examined among 128 recombinant inbred lines from four biparental crosses. The genotypic variation detected for drought score in the dry season was used to examine the reliability of the dry season screening method to estimate relative grain yield of genotypes under different types of drought stress in the wet season. Large genotypic variation for drought score existed in two experiments (A and B). However, there was no relationship between the drought scores of genotypes determined in these two experiments. Different patterns of development and severity of drought stress in these two experiments, i.e. slow development and mild plant water deficit in experiment A and fast development and severe plant water deficit in experiment B, were identified as the major factors contributing to the genotypes responding differently. Larger drought score in the dry season experiments was associated with lower grain yield under specific drought stress conditions in the wet season, but the association was weak to moderate and significant only in particular drought conditions. In most cases, a significant phenotypic and moderate genetic correlation between drought score in the dry season and grain yield in the wet season existed only when both drought score and grain yield of genotypes were affected by similar patterns and severity of drought stress in their respective experimental environments. The dry season environments used to measure genotypic variation for drought score should be managed to correspond to relevant types of drought environment that are frequent in the wet season. The efficiency of using the drought score as an indirect selection criterion for improving grain yield for drought conditions was lower than the direct selection for grain yield, and hence wet season screening with grain yield as a selection criterion would be more efficient. However, using drought score as a selection index, a larger number of genotypes can be evaluated than for wet season grain yield. Therefore, it is possible to apply higher selection intensities using the drought score system, and the selected lines can be further tested for grain yield in the wet season. (C) 2004 Elsevier B.V. All rights reserved.

Experimental evidence for multivariate stabilizing sexual selection

Stabilizing selection is a fundamental concept in evolutionary biology. In the presence of a single intermediate optimum phenotype (fitness peak) on the fitness surface, stabilizing selection should cause the population to evolve toward such a peak. This prediction has seldom been tested, particularly for suites of correlated traits. The lack of tests for an evolutionary match between population means and adaptive peaks may be due, at least in part, to problems associated with empirically detecting multivariate stabilizing selection and with testing whether population means are at the peak of multivariate fitness surfaces. Here we show how canonical analysis of the fitness surface, combined with the estimation of confidence regions for stationary points on quadratic response surfaces, may be used to define multivariate stabilizing selection on a suite of traits and to establish whether natural populations reside on the multivariate peak. We manufactured artificial advertisement calls of the male cricket Teleogryllus commodus and played them back to females in laboratory phonotaxis trials to estimate the linear and nonlinear sexual selection that female phonotactic choice imposes on male call structure. Significant nonlinear selection on the major axes of the fitness surface was convex in nature and displayed an intermediate optimum, indicating multivariate stabilizing selection. The mean phenotypes of four independent samples of males, from the same population as the females used in phonotaxis trials, were within the 95% confidence region for the fitness peak. These experiments indicate that stabilizing sexual selection may play an important role in the evolution of male call properties in natural populations of T. commodus.

And the beak shall inherit - evolution in response to invasion

The increased demographic performance of biological invaders may often depend on their escape from specifically adapted enemies. Here we report that native taxa in colonized regions may swiftly evolve to exploit such emancipated exotic species because of selection caused by invaders. A native Australian true bug has expanded it host range to include a vine imported from tropical America that has become a serious environmental weed. Based on field comparisons and historical museum specimens, we show that over the past 30-40 years, seed feeding soapberry bugs have evolved 5-10% longer mouthparts, better suited to attack the forest-invading balloon vines, which have large fruits. Laboratory experiments show that these differences are genetically based, and result in a near-doubling of the rate at which seeds are attacked. Thus a native biota that initially permits invasion may rapidly respond in ways that ultimately facilitate control.

CTL recognition of a bulged viral peptide involves biased TCR selection

MHC class I molecules generally present peptides of 8-10 aa long, forming an extended coil in the HLA cleft. Although longer peptides can also bind to class I molecules, they tend to bulge from the cleft and it is not known whether the TCR repertoire has sufficient plasticity to recognize these determinants during the antiviral CTL response. In this study, we show that unrelated individuals infected with EBV generate a significant CTL response directed toward an HLA-B*3501-restricted, 11-mer epitope from the BZLF1 Ag. The 11-mer determinant adopts a highly bulged conformation with seven of the peptide side chains being solvent-exposed and available for TCR interaction. Such a complex potentially creates a structural challenge for TCR corecognition of both HLA-B*3501 and the peptide Ag. Surprisingly, unrelated B*3501 donors recognizing the 11-mer use identical or closely related alpha beta TCR sequences that share particular CDR3 motifs. Within the small number of dominant CTL clonotypes observed, each has discrete fine specificity for the exposed-side chain residues of the peptide. The data show that bulged viral peptides are indeed immunogenic but suggest that the highly constrained TCR repertoire reflects a limit to TCR diversity when responding to some unusual MHC peptide ligands.

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.

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.