HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Figures Only Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (16) Citing Articles via Google Scholar Google Scholar Articles by Moreau, L. Articles by Gallais, A. Search for Related Content PubMed Articles by Moreau, L. Articles by Gallais, A. Agricola Articles by Moreau, L. Articles by Gallais, A.

CROP BREEDING, GENETICS & CYTOLOGY

Economic Efficiency of One Cycle of Marker-Assisted Selection

^a I.N.R.A.-U.P.S.-I.N.A.-P.G., Station de Génétique Végétale, Ferme du Moulon, F-91190 Gif-sur-Yvette, France
^b I.N.R.A., S.E.R.D. Département d'Economie et de Sociologie Rurales, Université Pierre Mendès France, BP 47, F-38040 Grenoble cedex 9, France

moreau{at}moulon.inra.fr

	ABSTRACT

TOP ABSTRACT INTRODUCTION Methods Results Discussion REFERENCES

 
The efficiency of molecular markers to improve genetic prediction has been proved by many studies. Nevertheless, the additional cost due to marker genotyping is seldom considered in the comparison between marker-assisted selection (MAS) and phenotypic selection. In the context of plant breeding, the relative cost efficiency of MAS in the first cycle of selection is evaluated through an analytical approach taking into account the effect of the experimental design (population size, number of trials, and replications per trial) on quantitative trait loci (QTL) detection. The same global cost is assumed for both methods. In a first step, the optimal allocation of the experimental resources is studied for each method before comparing them at their optimum. For traits sensitive to genotype x environment interactions, unreplicated trials are optimal for both methods but the optimal number of trials is different. It increases with the investment for phenotypic selection while it becomes nearly equal to one for MAS. The loss of efficiency due to non-optimal designs is evaluated. The expected economic return of MAS compared with phenotypic selection decreases with the cost of genotyping. When this cost is high, MAS interest is limited to traits with a low heritability, provided that the investment is high enough to evaluate a large population size, which is the necessary condition to explain with markers a large part of genetic variation. The maximal genotyping cost that is acceptable for MAS to be efficient is given for different values of investment and trait heritability.

Abbreviations: DH, doubled haploids • MAS, marker-assisted selection • QTL, quantitative trait locus

	INTRODUCTION

TOP ABSTRACT INTRODUCTION Methods Results Discussion REFERENCES

 
FOR THE LAST DECADE, rapid progress has been made in molecular marker techniques. For many species, the number of polymorphic markers is now sufficient to build dense maps and to localize individual QTL. Once marker-QTL associations have been detected, they can be used to improve the prediction of genetic values. This use of markers in selection schemes has received extensive interest in the recent past. Lande and Thompson (1990) proposed to identify marker-QTL associations by multiple regression of phenotype on marker type and to select individuals on an index combining phenotype and molecular value predicted with markers. The efficiency of this method of marker-assisted selection (MAS) relative to purely phenotypic selection has been widely studied, mainly in the case of populations derived from the cross of two inbred lines, through analytical approaches (Lande and Thompson, 1990; Luo et al., 1997; Moreau et al., 1998) and simulations (Zhang and Smith, 1992, 1993; Gimelfarb and Lande, 1994, 1995; Whittaker et al., 1995, 1997; Hospital et al., 1997). All these theoretical studies concluded that in many situations, MAS could be more efficient than phenotypic selection.

Despite these promising results, to our knowledge, few MAS experiments have been conducted with the approach proposed by Lande and Thompson (1990). This is certainly partly due to the cost of MAS experiments, which require both phenotypic and molecular evaluation. In all the theoretical studies mentioned above, the problem of the additional cost due to marker evaluation is not taken into account. Both methods of selection are assumed to be conducted in the same agricultural conditions, that is to say, assuming equal population sizes and identical experimental designs. In this situation, the additional cost associated with the marker evaluation is ignored, which leads to overestimating the interest of MAS compared with the phenotypic selection. Recently, Xie and Xu (1998) evaluated analytically the relative efficiency of different strategies of MAS, and incorporated the costs associated with measuring phenotype and scoring marker loci into the objective function, to maximize the gain per cost unit. Nevertheless, when estimating the genetic gain associated with MAS, the percentage of variance associated with the detected QTLs was assumed to be a known and fixed parameter. This would only be the case if the QTL detection and the estimation of the QTL effects were done in an infinite population. The same assumption was made by Knapp (1998) who compared the number of progenies that must be evaluated in phenotypic selection or MAS to retain a given proportion of superior genotypes. It is now well established that using finite populations generates important uncertainties with respect to QTL detection and effect estimation (Beavis, 1994). Thus, the size of the population under selection and also the experimental design, which affects the heritability of the trait, must be taken into account since they are deciding parameters of QTL detection and consequently MAS efficiency (Lande and Thompson, 1990; Zhang and Smith, 1992, 1993; Gimelfarb and Lande, 1994; Whittaker et al., 1995, 1997; Hospital et al., 1997; Moreau et al., 1998).

To be fair, MAS and phenotypic selection must be compared at their optimum for a given cost attributed to the selection process and for a given number of selected entries at the end of the cycle. For the same cost, phenotypic selection, which does not require marker evaluation, allows breeders to increase the population size and thus the intensity of selection compared with the one used for MAS. Alternatively, breeders can improve the phenotypic evaluation and thus the heritability of the trait by doing replications, as is usually done in plant breeding. This possibility of doing replications requires (i) either populations composed of fixed genotypes like inbred lines, or (ii) progeny tests, often used in QTL detection experiments (for instance, F2 population evaluated through the mean performance of derived F3 families). Instead of doing replications, other methods of prediction based on performances of relatives (classical full-sib or half-sib schemes) or on associated traits may be more adapted for improving the genetic gains for MAS and phenotypic selection (see Lande and Thompson, 1990) but will not be considered here. When it is possible to replicate individuals, several allocations of the experimental resources are possible (population size, number of locations, replications per location) with a fixed experimental cost. Before comparing the efficiency of MAS and phenotypic selection, the allocation of experimental resources, which leads to the highest genetic gain, must be determined for each method.

In the plant breeding context, the optimization of the number of replications and the population size in experiments of phenotypic selection has already been studied in several papers (Bos, 1983; Wricke and Weber, 1986; Gauch and Zobel, 1996). This has not yet been done in the MAS context. The determination of the optimal design requires one to be able to predict the genetic gain for a large range of situations in order to determine the one which gives the highest genetic progress. As simulations are too time consuming, the ideal is to use an analytical approach. In this paper, we use the approach developed by Moreau et al. (1998) which, unlike the formulation of Lande and Thompson (1990), fully takes into account the effect of population size. Moreover it was shown to give results consistent with those obtained by simulation with a more complex genetic model (Hospital et al., 1997).

Extending the study of Moreau et al. (1998), the aim of this paper is to investigate (i) the optimal allocation of experimental resources for MAS and phenotypic selection and, then (ii) to evaluate the domains where MAS is cost efficient. This work is restricted to the plant breeding context: we assume that the population is derived from a cross between two inbred lines and that the heritability of the trait can be improved by replicating each entry several times, in different field trials or within the same trial.

	Methods

TOP ABSTRACT INTRODUCTION Methods Results Discussion REFERENCES

 
Phenotypic and Marker-Assisted Selection
We consider that the aim of the selection is to improve the mean value of the population over several environments. For a pure phenotypic selection, the genetic value of an individual, i, of a population sample can be simply estimated from its average performance (P_i) within the whole experimental design. If each individual i is evaluated in L different trials, corresponding to different locations, and if its evaluation is replicated K times within each trial (divided, for instance, in K complete blocks), then:

(1)

where p_ilk is the kth performance of i evaluated within the trial l. The basic experimental unit used to evaluate p_ilk will be further referred to as a plot. Each plot is sown with only one genotype i but may contain one or several plants. For instance, in most species the yield performance of a genotype is evaluated through the bulked harvest of several plants.

Following Lande and Thompson (1990), in MAS, individuals are selected on an index including both phenotype and marker information:

(2)

where P_i is the average performance of i (see above), M_i is the sum of the effects of the markers associated to the trait and is called Molecular score, _p and _m are the weight coefficients of the index. , where q indicates markers significantly associated with the trait of interest, â_q is the estimated additive effect associated with one of the parental alleles at marker q, and _iq is a dummy variable that takes the value -1, 0, or 1 according to the number of copies of the given parental allele at marker q is 0, 1, or 2, respectively. With the aim of predicting the mean value of the progeny of a given individual in the next generation, only additive effects need to be evaluated. The weight coefficients that maximize the genetic gain are defined by:

(3)

where ²_m is the estimated percentage of phenotypic variance associated with markers and h²_b is the broad sense heritability associated with the average performance P_i in the experimental design.

(4)

²_G is the additive genetic variance common to all the environments, ²_GE is the variance of the genotype x environment interaction and ²_e is the residual variance due to experimental error. We can define another heritability that will be further mentioned as individual heritability and noted h²_i:

(5)

which represents the proportion of the phenotypic variation, observed at the level of an individual plot, due to genetic effects common to all the environments. This parameter does not depend on the experimental design and is specific to the trait.

Definition of Economic Parameters
Assuming that in phenotypic selection, the experimental costs are directly proportional to the number of plots, the total experimental resources I^P required to achieve one cycle of selection are: I^P = N^P L^P K^P C, where N^P is the total number of the genotypes of the population sample evaluated, C is the individual cost needed to evaluate the phenotypic value of one genotype in one plot (i.e., p_ilk), L^P and K^P are the number of trials and the number of replications per trial, defined above. The superscript P refers to the phenotypic selection.

For MAS, the investment ( I^M) needed to achieve one cycle is I^M = N^M (L^M K^M C + C_mark), where N^M, L^M and K^M are the N, L, and K parameters for MAS (indicated with the superscript M). C_mark is the cost associated with the genotyping of one individual. This quantity can be expressed as an equivalent number of plots that could be done with the same experimental resources: C_mark = C R_C, where R_C is the relative cost of genotyping one individual compared with the cost of evaluating one individual in one plot. Obviously, the cost of genotyping depends on the number of markers and on the individual cost of one marker-data point. In the first cycle of selection, marker-QTL associations have to be determined, which requires enough markers to have good coverage of the genome. Because this number depends on the genome length, for sake of simplicity, we did not evaluate the effect of this number and only considered in our formulations the global genotyping cost. With these notations, all the costs can be expressed in field plot unit. . The comparison of the two methods for given fixed experimental resources ( I₀) leads to the relationship:

(6)

As R_C 0, the above relation implies that the two methods are conducted with different sets of L, K, and N values. For a given I₀ and R_C, the problem is to determine the values of N, L, and K that give the highest genetic gain for each method.

Evaluation of the Genetic Gains
Assuming that the selected trait is normally distributed in the population, the expected genetic gain after one cycle of phenotypic selection is equal to (Falconer, 1981) G_P = ih_bG, where i is the intensity of selection in standard unit, h_b is the square root of the broad sense heritability defined above, and _G is the genetic standard deviation. In an infinite population, i is equal to z/p, where z is the height of the ordinate at the point of truncation of a normal distribution and p is the proportion of selected individuals. The broad sense heritability h²_b depends on the experimental design (L^P and K^P) and on the proportions , and , which are specific to the trait.

For MAS, the analytical expression of the expected genetic gain is more complex and implies more parameters. Among them, the population sample size has a major effect. It strongly affects the power of QTL detection and the precision of the estimation. In finite population samples, when QTL detection and genetic prediction are made on the same sample, the selection of the markers included in the index (Eq. [2]) leads to an overestimation of the effects associated with markers (Lande and Thompson, 1990; Beavis, 1994). This bias of selection is difficult to correct (Whittaker et al., 1997; Melchinger et al., 1998) and was taken into account to evaluate MAS efficiency (Moreau et al., 1998). Other parameters influence the MAS genetic gain: the genetic determinism of the trait (the number of QTLs involved, their localization, the distribution of their effects), the molecular map (the number of markers and their position relative to QTLs), the method of QTL detection used, and the rule for the selection of markers in the index.

To take these parameters into account, we used the approach of Moreau et al. (1998) (see this paper for a full description and discussion of the model). We consider that the genetic variability of the trait is explained by x_l independent QTLs of equal effect, with no epistasis. Two genetic models are considered: x_l = 5 or 10 QTLs. Each QTL is linked with a single marker. Other markers (30 - x_l), unlinked with QTL, are added to the model to simulate possible false detections. The total number of 30 markers was chosen to roughly correspond to 10 chromosomes and three nearly independent markers located on each. Marker-QTL associations are detected by a simple regression, taking a type I risk level equal to 5%, chosen following the results of Moreau et al. (1998). The rate of recombination (noted r) between each QTL and its linked marker is assumed to be the same for each QTL. This parameter determines m², the maximum percentage of genetic variance that can be explained by markers. In a population composed of doubled haploids (DH), m² is equal to (1 - 2r)². In this study, we do not address the problem of the optimization of the map density. Thus, we consider a fixed marker-QTL distance, corresponding to an m² value equal to 0.8. In a DH population, this value is obtained with an average density of one marker every 20 centimorgans, which was shown to be nearly optimal (Gimelfarb and Lande, 1994; Darvarsi et al., 1994). Moreau et al. (1998) highlighted that a small number of QTLs with equal effects provides a good approximation of the results obtained with more QTLs and a geometric distribution of their effects, corresponding to the same equivalent number of QTL (as defined by Lande, 1981, and further noted x_e). If QTLs have variable contributions to the genetic variance, the equivalent number of QTLs is expected to be rather small (between 3 and 20) (Chevalet, 1994). In many studies addressing the potential benefit of MAS by means of simulations, an equivalent number of 10 QTLs is assumed (Gimelfarb and Lande, 1994; Hospital et al., 1997). This model was shown to give consistent results with such simulations conducted with a more realistic genetic model and marker allocation, and a multiple regression for the QTL detection. Thus, the model considered here is certainly too simple, but rapidly yields results that can reasonably be extrapolated to more realistic situations, at least in the first cycle.

In all this work, for both methods of selection, we consider that a fixed number of individuals (n_s) is selected (i.e., the same effective genetic size). Following Eq. [6], the proportion of selected individuals is p^P = n_s/N^P = K^P L^P (n_s/I₀) for phenotypic selection, and p^M = (K^M L^M + R_C) (n_s/I₀) for MAS. The results that will be presented are obtained with a fixed number of selected individuals n_s equal to 10. This number represents a compromise between short term selection efficiency and the risk of genetic drift in the long term. As the population sample size is limited when the investment is low, the intensity of selection i corresponding to the p, found in the table given by Falconer (1981) for an infinite population, was corrected by the approximate formulation given by Burrows (1972) to take the population sample size into account: i_N = i - (N - n_s)/[2n_s (N + 1)i], where i and i_N are the intensities of selection for an infinite population and for population sample size of N, respectively.

Optimization of the Experimental Parameters—Relative Efficiency of MAS
For a given trait (characterized by h²_i and ) and given values of I₀ and R_C, we first determine the combination of the three parameters N, K (between 1 and 25), and L (between 1 and 25), which maximizes the genetic gain. Because a strict analytical solution of this problem requires complex mathematical derivations, the determination of the optimal allocation of the experimental resources is made numerically by comparing the genetic gain obtained with different experimental designs. The values (L, K) leading to the highest genetic gain are considered as optimal for a given method. The genetic gains (G_M and G_P) obtained for each method with an optimal allocation of the same investment are later compared. The relative cost efficiency of MAS (further noted RE) is defined by the ratio:

(7)

The program was implemented in FORTRAN. For MAS, the genetic gain evaluation was performed by the statistical facilities of SAS (1988). Because these two programs were not easy to connect, MAS genetic gains were first evaluated for a large range of values of N and h²_b and stocked in a data file. The genetic gain for a particular set of values (N, h²_b) was estimated by linear approximations between the genetic gains found in the data file that were evaluated with the nearest flanking values of these parameters.

	Results

TOP ABSTRACT INTRODUCTION Methods Results Discussion REFERENCES

 
Optimal Allocation of the Experimental Resources
For MAS and phenotypic selection, when non-zero genotype x environment interactions are considered in the model, the highest genetic gain is always obtained when the genotypes are evaluated only once (K = 1) in each field trial (results not shown). This conclusion already drawn by Wricke and Weber (1986) for a phenotypic selection appears to be also valid for MAS. When genotype x environment interactions are considered as null, there is no need to evaluate the population in several locations. For practical reasons, the best strategy is to put all the genotypes in the same trial (L = 1) and, if needed, to replicate them within the field trial.

As a consequence, for both methods of selection, all results presented can be interpreted either as (i) the optimal number of trials (L) with one replication per trial if the trait is affected by genotype x environment interactions or (ii) the optimal number of replications (K) within one trial if the trait is not affected by such interactions. For sake of simplicity, we will henceforth always use the term replication without mentioning whether it corresponds to replications in different trials or within the same one. The other consequence is that the results do not depend on the ratio , the only effective parameter being h²_i.

The comparison of the results obtained for each method shows that the optimal number of replications is quite different for MAS and phenotypic selection (Fig. 1) . For phenotypic selection (see graph a), when the individual heritability is high (above 0.7), doing replications does not sufficiently increase the broad sense heritability to compensate for the decrease in selection intensity associated with a smaller population sample size. For lower heritabilities, replications can increase genetic gain provided that the total investment is high enough. Thus, the optimal number of replications increases with the investment and decreases with the individual heritability. So, for very low individual heritabilities (h²_i < 0.10), it is optimal to do a high number of replications (above 6) even with a medium investment (equivalent to 600 plots).

View larger version (26K): [in this window] [in a new window]   Fig. 1 Optimal number of replications for different investments, (I0, abscissa) and different individual heritabilities, (h2i, ordinate), for phenotypic selection (graph a) and marker-assisted selection with RC value equal to 1 (graph b), 4 (graph c), and 8 (graph d). We assumed that 10 individuals are selected and that the determinism of the trait implies an equivalent number of 5 QTLs

The fact that two replications instead of one are optimal when R_C = 8, even for high heritabilities, is due to the small impact of doing one or two replications on the population size in this case. When I₀ = 4000, doing one replication instead of two only increases the population sample size from 400 to 444 (+11%) when R_C = 8, while it increases the population sample size from 600 to 800 (+20%) when R_C = 4, and from 1333 to 2000 (+50%) when R_C = 1. Thus, for high R_C values, the gain in the selection intensity and in the power of QTL detection is small. Moreover, when h²_i is high, the phenotype is a good predictor of breeding values. The additional information that can be provided by markers is low and a moderate increase in the power of QTL detection has only a low impact on breeding value prediction and does not compensate for the broad sense heritability decrease.

The results obtained with x_e = 10 QTLs are not presented but shows similar trends to those observed with x_e = 5 QTLs, except that the optimal number of replications is slightly higher. With x_e = 10 QTLs, individual QTL effects are smaller. The power of QTL detection can less easily be improved by increasing the population sample size. This enhances the interest of doing replications.

For practical reasons, it is not always possible to achieve the optimal experimental design. It is thus important to evaluate the loss of efficiency associated with a non optimal number of replications. For I₀ = 2000, it appears that when the number of replications is close to the optimum, the loss of efficiency is generally small for both methods, especially for MAS when R_C > 1 and h²_i is high (Table 1) . When the marker genotyping is cheap (R_C = 1), it is always optimal to do one replication. In this case, adding one replication reduces the population sample size (from 1000 to 667) and significantly reduces MAS efficiency. For higher R_C values, the number of replications has a weaker impact on the population sample size and its effect on MAS efficiency becomes negligible near the optimum. Thus, the experimental design can be adapted to specific constraints without reducing the MAS efficiency.

View this table: [in this window] [in a new window]   Table 1 Genetic gains (G), in genetic standard deviation unit, and population sample sizes (N) corresponding to a given investment I0 = 2000 but different numbers of replications (rep) with phenotypic selection and MAS, for different individual heritabilities (h2i) and different relative costs of genotyping (RC). For MAS, an equivalent number of 5 QTLs is assumed

Relative Cost Efficiency of MAS
Once the optimal number of replications for each method of selection is determined for a given global cost, the ratio of the genetic gains at optimum leads to the RE of MAS. For both genetic models (x_e = 5 or 10 QTLs), the general trend is towards an increase of MAS relative efficiency with I₀ and a decrease with h²_i (Fig. 2) . When h²_i values are high, RE is always below 1. In such a case, the phenotype is a good predictor of the breeding value, the gain in accuracy provided by markers is too low to compensate for the additional experimental cost associated with the marker-genotyping. Thus, MAS is only valuable for low or medium h²_i, when markers can improve the accuracy of prediction. Nevertheless, as mentioned by Moreau et al. (1998), a small decrease in MAS RE can be observed for very low heritabilities because of the low power of QTL detection in this situation. MAS RE for low or medium heritabilities highly depends on the experimental resources (I₀), especially when the relative cost of marker genotyping is high. With x_e = 5 QTLs and R_C = 1, MAS is justified in nearly all the situations (except when h²_i > 0.75 and I₀ < 200). When R_C = 8, the domain of MAS efficiency is restricted to very high I₀ values (>1700) and low heritabilities. Thus, even for a low heritability, MAS could be cost inefficient if I₀ is too small. For instance, when I₀ = 400 and h²_i < 0.15, RE is smaller than 0.6. The population sample size allowed with such an investment is too limited for the markers to improve sufficiently the accuracy of genetic value prediction. With the same investment, phenotypic selection allows breeders to do more replications, which significantly increases the broad sense heritability and thus the efficiency of the phenotypic selection.

View larger version (26K): [in this window] [in a new window]   Fig. 2 Domains of cost efficiency of MAS for different values of the relative cost of genotyping (RC = 1, 4, or 8) and two genetic models implying an effective number of 5 or 10 QTLs, as a function of the total investment (I0, abscissa) and the individual heritability of the trait (h2i, ordinate)

Genetic determinism affects the domains of MAS efficiency (see the upper and lower graphs in Fig. 2). MAS appears less efficient when x_e = 10 QTLs than when x_e = 5 QTLs. This result is not surprising since when more QTLs are implied in the genetic determinism, their individual effects decrease and consequently their probability of being detected. Despite the reduction in domains of MAS efficiency, the general trends observed when x_e = 5 QTLs are consistent with those observed when x_e = 10 QTLs. The differences between the two genetic situations are more striking for low h²_i values, since, in this case, the weight given to the markers in the index (Eq. [2]) is important and the power of QTL detection is a deciding parameter. When , MAS requires an investment approximately twice as high when x_e = 10 QTLs than when x_e = 5 QTLs to become cost efficient.

	Discussion

TOP ABSTRACT INTRODUCTION Methods Results Discussion REFERENCES

 
In our study, we used simple genetic and economic models. The assumptions made for the genetic model have been discussed above and a full discussion can be found in Moreau et al. (1998). Other simplifications concern the economic part of the model. The costs of genotyping and phenotyping one individual are assumed to be independent of the population size. This is not realistic since the fixed costs linked to the necessary investment in machines, staff, etc., are the same whatever the population sample size. Thus, the cost for one data point (phenotype or marker-genotype) decreases when more individuals are analyzed simultaneously. For the agronomic evaluation, doing several trials instead of a large one is more costly, and the fixed cost associated with a trial depends on the location and the experimental conditions (for instance irrigation or not). Since all these economic aspects are difficult to quantify, they were not included in our model.

Despite the simplifications made, the results obtained with our model can be reasonably extended to real situations and highlight the major parameters that must be taken into account before undertaking a MAS experiment. When genotype x environment interactions are included in the model, it is always optimal to perform only one replication per trial. When the investment is high enough, it appears optimal to do only a small number of trials (one or two), even when genotype x environment interactions are important.

Nevertheless, such a design presents several drawbacks. Experimental problems during sowing or harvest may occur and lead to the elimination of a trial from the analysis. Traits like lodging in maize (Zea mays L.) or disease tolerances in natural conditions strongly depend on uncontrolled environmental effects and may be totally unexpressed in a given trial, in a given year. Performing only one large unreplicated trial could then be very risky. Moreover, spatial heterogeneity is more likely to occur in large field trials and without replication, it is difficult to estimate the accuracy of trials.

For these reasons, it may be less risky to increase the number of trials compared with the optimum, since it only slightly decreases the efficiency of MAS (Table 1). It may also be useful to use checks and/or to replicate a small subset of the population sample within each trial. This permits to estimate the trial accuracies and can also be used to control trial heterogeneity using spatial models, as was shown by Cullis et al. (1989) for phenotypic selection and by Moreau et al. (1999) in MAS context.

Even if a nearly optimal design may be preferred to reduce experimental risks, it still remains that the optimum for MAS is to use a reduced number of trials when compared with phenotypic selection. This strategy favors the population sample size, which is the determining parameter of MAS efficiency. This is consistent with the results found in the context of QTL detection experiments by Knapp and Bridges (1990), who showed that it is more efficient to increase the population sample size rather than the accuracy of phenotypic evaluation.

Our results show that the relative efficiency of MAS is reduced when cost is taken into account. One can further add that for low heritabilities, doing replications is not the best strategy for phenotypic selection. Breeding schemes based on performances of relatives or associated traits are more efficient. Such schemes were not considered in this study but we can infer that they would further reduce the interest of MAS (as shown by Lande and Thompson, 1990).

According to our study, the benefit of MAS in the first generation is rather small compared with the conclusions drawn by Xie and Xu (1998) or Knapp (1998). These authors do not consider the same criterion as ours to assess the cost efficiency of MAS. Roughly speaking, Knapp (1998) considers the ratio of the cost needed for MAS and phenotypic selection to achieve the same genetic gain, while Xie and Xu (1998) consider the genetic gain per cost unit of each method, both conducted with a same population size and experimental design. Moreover, Knapp (1998) evaluated the effect of the objective of selection in the comparison between MAS and phenotypic selection that was not taken into account in the other studies. Nevertheless, the difference between the results of these studies and ours mainly comes from the fact that these authors neglect the impact of the experimental design on QTL detection and the uncertainty of the estimated QTL effects. Even if new methods of detection based on interval mapping with a multilocus analysis (Jansen and Stam, 1994; Zeng, 1994) increase the resolution and the power of QTL detection, the assumption of known QTL effects is still optimistic. The true nature of the QTLs is still undetermined, even if genes of known function are often found to be located near the detected QTLs and provide good candidates. One can expect that progress in molecular biology or in QTL detection will help to determine more precisely the QTL positions and their effects but until now, QTL detection is still necessary. Because of sampling problems, the QTL effects estimated in one experiment are known to be biased and can hardly be inferred from QTL effects previously found in another one, even made in the same genetic background (Beavis, 1994). To limit the impact of sampling error when estimating MAS efficiency, Melchinger et al. (1998) advised the use of QTL parameters obtained from a validation experiment rather than from the original discovery experiment. This strategy is efficient for taking into account sampling bias but requires two population samples (one for QTL detection and one for validation). Thus, its benefit on the economic efficiency of MAS is difficult to evaluate. For these reasons, we find it to be more realistic to estimate the QTL effects for the first cycle in the population under selection. For the same reasons we did not consider in our comparison the case of a selection on markers only as was done by Zhang and Smith (1992) or Xie and Xu (1998). When the QTL positions are unknown, both phenotypic evaluation and complete marker genotyping of the population sample are required in the first generation and the more efficient strategy is to select individuals on an index combining markers and phenotype (Zhang and Smith, 1992).

The most important limit to the cost efficiency of MAS comes from the costs of genotyping. In maize selection, considering that 80 markers are necessary, one can estimate that the relative cost with RFLP markers (restricted fragment length polymorphism) is near 8 for usual traits like grain yield in maize ($2 for each marker data point and $20 for individual phenotypic value). With such marker costs, an R_C value equal to 1 can only be achieved if the phenotypic evaluation is very expensive, around 2 x 80 = $160 per plot. Thus, for most traits, small values of R_C are certainly unrealistic and MAS interest in the first cycle is restricted to very high investments. New marker techniques, for instance microsatellites, can be automated (Schwemgel et al., 1994; Perlin et al., 1995). One can expect that the improvement of marker techniques will decrease the genotyping costs and consequently increase the domains of MAS efficiency. To assess whether it is worth doing a MAS experiment, one can determine for each set of parameters h²_i and I₀, the highest R_C value that is acceptable for MAS to be cost efficient (Table 2) . The maximum acceptable value to R_C increases when the individual heritability becomes lower, when the experimental resources become higher, and when the genetic determinism of the trait becomes simpler. Even with a R_C value of 8 and an investment equivalent to 3500 field plots, MAS is efficient for a complex trait when the individual heritability is equal to 0.15. Such heritability may seem very small but it is not unrealistic for traits which are influenced by genotype x environment interactions (like stalk lodging or even yield in maize). With these parameters, the optimal design consists of evaluating 292 individuals in five trials, which is rather close to designs that are now used in QTL detection experiments.

Received for publication January 19, 1999.

	REFERENCES

TOP ABSTRACT INTRODUCTION Methods Results Discussion REFERENCES

 

• Beavis, W.D. 1994. The power and deceit of QTL experiments: Lessons from comparative QTL studies. p. 250–266. In D. Wilkinson, (ed.) Proc. Annu. Corn Sorghum Indus. Res. Conf., 49th, Chicago, IL. Dec 1994. Am. Seed Trade Assoc., Washington, DC.
• Bos I. The optimum number of replications when testing lines or families on a fixed number of plots. Euphityca 1983;32:311-318.
• Burrows P.M. Expected selection differentials for directional selection. Biometrics 1972;22:1091-1100.
• Chevalet C. An approximate theory of selection assuming a finite number of quantitative trait loci. Gen. Sel. Evol. 1994;26:379-400.
• Cullis B.R., Lill W.J., Fisher J.A., Read J.B., Gleeson A.C. A new procedure for the analysis of early generation variety trials. Appl. Statistics 1989;38:361-375.
• Darvarsi A., Soller M. Optimum spacing of genetic markers for determining linkage between marker loci and quantitative trait loci. Theor. Appl. Genet. 1994;89:351-357.
• Falconer D.S. Introduction to quantitative genetics, 2nd ed New York: Longman, 1981.
• Gauch H.G., Zobel R.W. Optimal replication is selection experiments. Crop Sci. 1996;36:838-843.[Abstract/Free Full Text]
• Gimelfarb A., Lande R. Simulation of marker-assisted selection in hybrid populations. Genetical Res. 1994;63:39-47.
• Gimelfarb A., Lande R. Marker-assisted selection and marker-QTL associations in hybrid population. Theor. Appl. Gen. 1995;91:522-528.
• Hospital F., Moreau L., Charcosset A., Gallais A. More on the efficiency of marker assisted selection. Theor. Appl. Gen. 1997;95:1181-1189.
• Jansen R.C., Stam P. High resolution of quantitative traits into multiple loci via interval mapping. Genetics 1994;136:1447-1455.[Abstract]
• Knapp S.J. Marker-Assisted selection as a strategy for increasing the probability of selecting superior genotypes. Crop Sci. 1998;38:1164-1174.[Abstract/Free Full Text]
• Knapp S.J., Bridges W.C. Using molecular markers to estimate quantitative trait locus parameters: power and genetic variances for unreplicated and replicated progeny. Genetics 1990;126:769-777.[Abstract]
• Lande R. The minimum number of genes contributing to quantitative variation between and within populations. Genetics 1981;99:541-553.[Abstract/Free Full Text]
• Lande R., Thompson R. Efficiency of marker-assisted selection in the improvement of quantitative traits. Genetics 1990;124:743-756.[Abstract]
• Luo Z.W., Thompson R., Woolliams J.A. A population genetics model of marker-assisted selection. Genetics 1997;146:1173-1183.[Abstract]
• Melchinger A.E., Utz H.F., Schön C.C. Mapping using different testers and independent population samples in maize reveals low power of QTL detection and large biais in estimates of QTL effects. Genetics 1998;149:383-403.[Abstract/Free Full Text]
• Moreau L., Charcosset A., Hospital F., Gallais A. Marker-assisted selection efficiency in populations of finite size. Genetics 1998;148:1353-1365.[Abstract/Free Full Text]
• Moreau L., Monod H., Charcosset A., Gallais A. Marker-assisted selection with spatial analysis of unreplicated field trials. Theor. Appl. Genet. 1999;98:234-242.
• Perlin M.W., Lancia G., See-Kiong N. Toward fully automated genotyping: Genotyping microsattelites markers by deconvolution. Am. J. Hum. Genet. 1995;57:1199-1210.[Web of Science][Medline]
• SAS Institute Inc. 1988. SAS language: Reference. Version 6. First ed. SAS Institute Inc., Cary, NC.
• Schwemgel D.A., Jedlick A.E., Nanthakumar E.J., Weber J.L., Levitt R.C. Comparison of fluorescence-based semi-automated genotyping of multiple microsattelite loci with autoradiographic techniques. Genomic 1994;22:46-54.[Web of Science][Medline]
• Wricke G., Weber W.E. Quantitative genetics and selection in plant breeding. Berlin: Walter de Gruyter, 1986.
• Whittaker J.C., Curnow R.N., Haley C.S., Thompson R. Using marker-maps in maker-assisted selection. Genetical Res. 1995;66:255-265.
• Whittaker J.C., Haley C.S., Thompson R. Optimal weighting on information in marker-assisted selection. Genetical. Res. 1997;69:137-144.
• Xie C., Xu S. Efficiency of multistage marker-assisted selection in the improvement of multiple quantitative traits. Heredity 1998;80:489-498.
• Zeng Z.B. Precision mapping of quantitative trait loci. Genetics 1994;136:1457-1468.[Abstract]
• Zhang W., Smith C. Computer simulation of marker-assisted selection utilizing linkage disequilibrium. Theor. Appl. Genet. 1992;83:813-820.[Web of Science]
• Zhang W., Smith C. Simulation of marker-assisted selection utilizing linkage disequilibrium: the effect of several additional factors. Theor. Appl. Genet. 1993;86:492-496.[Web of Science]

This article has been cited by other articles:

				G. Abalo, P. Tongoona, J. Derera, and R. Edema A Comparative Analysis of Conventional and Marker-Assisted Selection Methods in Breeding Maize Streak Virus Resistance in Maize Crop Sci., March 17, 2009; 49(2): 509 - 520. [Abstract] [Full Text] [PDF]

				J. J. Ceron-Rojas, F. Castillo-Gonzalez, J. Sahagun-Castellanos, A. Santacruz-Varela, I. Benitez-Riquelme, and J. Crossa A Molecular Selection Index Method Based on Eigenanalysis Genetics, September 1, 2008; 180(1): 547 - 557. [Abstract] [Full Text] [PDF]

				C. Zhang, B. Tar'an, T. Warkentin', A. Tullu, K. E. Bett, B. Vandenberg, and D. J. Somers Selection for Lodging Resistance in Early Generations of Field Pea by Molecular Markers Crop Sci., January 24, 2006; 46(1): 321 - 329. [Abstract] [Full Text] [PDF]

				L. Flagel, J. R. Christensen, C. D. Gustus, K. P. Smith, P. M. Olhoft, D. A. Somers, and P. D. Matthews Inexpensive, High Throughput Microplate Format for Plant Nucleic Acid Extraction: Suitable for Multiplex Southern Analyses of Transgenes Crop Sci., August 26, 2005; 45(5): 1985 - 1989. [Abstract] [Full Text] [PDF]

				S. Crepieux, C. Lebreton, B. Servin, and G. Charmet Quantitative Trait Loci (QTL) Detection in Multicross Inbred Designs: Recovering QTL Identical-by-Descent Status Information From Marker Data Genetics, November 1, 2004; 168(3): 1737 - 1749. [Abstract] [Full Text] [PDF]

This Article Abstract Figures Only Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (16) Citing Articles via Google Scholar Google Scholar Articles by Moreau, L. Articles by Gallais, A. Search for Related Content PubMed Articles by Moreau, L. Articles by Gallais, A. Agricola Articles by Moreau, L. Articles by Gallais, A.