Evolved polygenic herbicide resistance in Lolium rigidum by low-dose herbicide selection within standing genetic variation

The interaction between environment and genetic traits under selection is the basis of evolution. In this study, we have investigated the genetic basis of herbicide resistance in a highly characterized initially herbicide-susceptible Lolium rigidum population recurrently selected with low (below recommended label) doses of the herbicide diclofop-methyl. We report the variability in herbicide resistance levels observed in F1 families and the segregation of resistance observed in F2 and back-cross (BC) families. The selected herbicide resistance phenotypic trait(s) appear to be under complex polygenic control. The estimation of the effective minimum number of genes (NE), depending on the herbicide dose used, reveals at least three resistance genes had been enriched. A joint scaling test indicates that an additive-dominance model best explains gene interactions in parental, F1, F2 and BC families. The Mendelian study of six F2 and two BC segregating families confirmed involvement of more than one resistance gene. Cross-pollinated L. rigidum under selection at low herbicide dose can rapidly evolve polygenic broad-spectrum herbicide resistance by quantitative accumulation of additive genes of small effect. This can be minimized by using herbicides at the recommended dose which causes high mortality acting outside the normal range of phenotypic variation for herbicide susceptibility.


Introduction
The great advantage of herbicides is their ability to selectively remove weed plants from crop fields and for this reason herbicides are used worldwide (Oerke 2006). However, recurrent and persistent herbicide use has resulted in the evolution of resistance in many weed species (Heap 2012). The speed of resistance evolution is influenced by specifics of the herbicide selection (use history, dose applied, associated agronomic practices), the biology of the plant species under selection (population growth rate, genetic diversity, reproductive mode, etc.) and population genetic factors (Jasieniuk et al. 1996). Of the range of gene traits that can endow herbicide resistance often there is target-site resistance resulting from a point mutation in a gene encoding for a specific herbicide target enzyme (Powles and Yu 2010). Target-site resistance is thus usually single gene inherited (Darmency 1994). Conversely, resistance can be nontarget-site based. For example, cytochrome P450 mono-oxygenases or glutathione S-transferases can metabolize certain herbicides and evolved resistance can be due to enhanced rates of herbicide metabolism endowed by these enzymes (Yuan et al. 2007). In evolved herbicideresistant populations of the important cross-pollinated grass weeds Alopecurus myosuroides (Huds.) and Lolium rigidum (Gaud.), enhanced rates of herbicide metabolism have been shown (Preston et al. 1996;Reade et al. 2004). In these cases of nontarget-site herbicide resistance mediated by herbicide metabolism, individual genes can endow resistance to specific herbicides (Preston 2003) and/or there can be more complex genetic linkages (Busi et al. 2011b). Diverse patterns of herbicide resistance can be evident at both the individual and population level (Petit et al. 2010a,b).
Here, we are concerned with herbicide resistance evolution from selection at low herbicide doses. Our studies conducted with L. rigidum have shown that recurrent low herbicide dose selection, within the range of quantitative variation for herbicide response evident within a small population, can rapidly lead to herbicide resistance evolution (Neve and Powles 2005;Powles 2009, 2011;Manalil et al. 2011). In small populations, a few generations of recurrent selection can lead to significant phenotypic shifts depending on the extent of intra-population genetic variation and heritability of traits. This has been shown in plants and insects and interpreted as incremental stacking of several genes of minor effect (Ffrench-Constant et al. 2004). In addition to selection, founder effects and genetic drift may influence the dynamics of resistance evolution (Falconer 1981). Here, we investigated the genetic basis of herbicide resistance as selected by low herbicide doses in cross-pollinated, genetically diverse L. rigidum. We report the herbicide resistance observed in the both herbicide-resistant R and initially susceptible S parental population, the variability of F 1 families, the polygenic segregation found in F 2 and back-cross (BC) families and discuss the management of low-dose-selected polygenic herbicide resistance.

Parental lines
The well-characterized herbicide-susceptible L. rigidum parental population VLR1 (hereinafter referred to as S) was exposed to three cycles of recurrent selection with low (below the recommended label dose) doses of the ACCase-inhibiting herbicide diclofop-methyl (selecting agent) as described in (Neve and Powles 2005). The three-time selected line VLR1 (0.1 0.5 2.0) which exhibited the highest level of phenotypic resistance was chosen for this genetic study (hereinafter referred to as R). The coefficient 0.1, 0.5 and 2.0 represents the proportion of the recommended label dose of diclofop-methyl (375 g ha À1 ) applied to plants at the first, second and third cycles of recurrent selection, respectively. Plant survival at 0.1, 0.5 and 2.0 dose was 36%, 33% and 44%, respectively.
Generation of F 1 families R plants as described above were treated with a single herbicide dose (188 g diclofop-methyl ha -1 ) to confirm resistance, and six randomly chosen resistant R plants were each paircrossed to one plant of the susceptible original unselected S parental to produce a total of 6 F 1 pair crosses. Seed progeny Table 1. Statistical differences between survival proportion was also assessed by chi-square tests performed by using the statistical software R with the command prop.test. A pooled chi-square value was calculated considering the sum of all the survivors in F 1 families and a heterogeneity chisquare test was performed to compare the segregation frequencies obtained in each family (Sokal and Rohlf 1969 Generation of F 2 families from F 1 cloned resistant phenotypes Three types of F 2 families (L, M or H) were generated by pair-crossing two cloned plants within the same F 1 family that survived the specific L, M or H diclofop-methyl dose. Type L F 2 families were generated by pair-crossing cloned plants from F 1 families that were only able to survive 188 g diclofop-methyl ha À1 (cloned counterparts were killed at 375 and 1500 g ha À1 ). Type M F 2 families were generated by pair-crossing cloned plants surviving 188 and 375 g diclofop-methyl ha À1 but killed at the highest dose (1500 g ha À1 ). Type H F 2 families were generated by paircrossing cloned plants that were able to survive the high dose of 1500 g diclofop-methyl ha À1 .

Generation of back-cross families
Plants from ♀ R F 1 families were treated with 188 g diclofop-methyl ha À1 at the two-leaf stage to eliminate susceptible individuals. The mean plant survival ratio in F 1 families was 75%, significantly different from survival in the S (v 2 = 59.1; P < 0.01), but not significantly different to survival observed in the R plants (v 2 = 0.68; P = 0.40). This suggested the phenotypic herbicide resistance trait(s) were endowed by dominant gene(s). As suggested by Tabashnik (2001), F 1 survivors of families 5 and 6 were pair back-crossed to plants of the original parental S population. Seeds were collected from both plants (♀, ♂ S).

Data analysis
Effective number of resistance genes and scaling test analysis The effective (minimum) number of genes (N E ) for resistance to diclofop-methyl in parental, F 1 and F 2 families was estimated as suggested by Lande (1981) as: l P2 -l P1 represent the difference in the mean phenotypic value (aboveground plant biomass) of parental populations (R and S) and r 2 is the difference in variance for the same phenotypic trait between F 2 and F 1 families.

Scaling test
Additive and dominance effects of the gene(s) endowing diclofop-methyl resistance were assessed by a scaling test as described in Mather and Jinks (1982). Observed mean phenotypic values (aboveground biomass) following herbicide treatment of parental R and S, F 1 , F 2 and BC families were tested for conformity to an additive-dominance model. The quantities A and C and their variances were calculated as described by Mather's equations (Mather and Jinks 1982) and test the adequacy of an additive-dominance model to the observed data in this study. If the model is adequate, A and C are equal to zero, and a Student's t-test was used to assess values for A and C were significantly different from zero. Pooled aboveground dry biomass data of F 1 F 2 and BC were used to calculate the parameters: Maternal effects and herbicide dose-response of parental, F 1 , F 2 and BC genetic families The resistant parental R, the unselected susceptible parent S and four F 1 families (no. 5 ♀, ♂ R1 F 1 , and no. 6 ♀, ♂ R2 F 1 ), six F 2 families (no. 5 and no. 6 for each family L, M and H) and two ♀ BC families (no. 5 and 6) were tested in final diclofop-methyl doseresponse studies conducted at different times but under identical glasshouse conditions. Seeds were germinated on 0.6% (w/v) agar. Seedlings were transplanted and grown in 18-cm diameter pots containing a potting mix (50% peatmoss, 25% river sand and 25% pine bark) and maintained in glasshouse conditions. Plants at the two-leaf stage were treated with 0, 46, 188, 375 or 1500 g diclofop-methyl ha À1 (recommended label dose 375 g ha À1 ). In each experiment, there were three replicates per treatment with a minimum of 37 plants treated at each herbicide dose depending on seed availability and germination. Plant survival and aboveground dry biomass were evaluated 21 days after herbicide treatment.
As described in Busi et al. (2011b), maternal effects of diclofop-methyl resistance genes were assessed by nonlinear regression analysis. Data from the two experiments were pooled and plant survival data were expressed as percentages. Data sets were fitted to a three-parameter log-logistic model: where the parameter d is the upper limit, b is the slope of the curve, x is the herbicide dose, and e is the dose producing a 50% reduction in response. Regression assumptions were held under square root data transformation (i.e. Box-Cox transformation lambda k = 0.5) (Onofri et al. 2010). Statistical differences in plant survival between the parental and F 1 families over a range of the herbicide diclofopmethyl doses were assessed by a lack-of-fit F-test applied to data sets fitted with the above nonlinear logistic model.

Resistance segregation in F 2 and BC families
The segregation of genetic traits endowing diclofop-methyl resistance was assessed along a gradient of doses: 188 (low), 375 (medium) or 1500 (high) g diclofop-methyl ha À1 . Diclofop-methyl resistance was assumed to be endowed by at least one incompletely dominant additive resistance allele. As described by Tabashnik (2001) and Preston (2003), the segregation analysis in F 2 and BC families was based on survival as resistant (alive) or susceptible (dead) plants compared with the expected survival/mortality H o ratios (Table  S1). The expected survival ratio in F 2 and BC families was corrected by weighting according to the observed performances of R and S parents and F 1 families (Table S1). Thus, the expected F 2 or BC survival values as number of plants are calculated with the total number of plants herbicidetreated multiplied by the theoretical one, two or three segregation ratios, respectively, (e.g. for one-gene model that ratio is 0.25R:0.5F1:0.25S) multiplied by the observed survival (%) in R, F1 and S at that specific dose (see example in Table S1). As described by Busi et al. (2011b), for each segregating F 2 and BC family, a goodness of fit chi-square (v 2 ) test was used to compare the observed plant survival with the expected calculated values according to one-, two-or three-resistance gene segregation models (Tables 1-3; Table  S1). P-values were obtained indicating the probability of type II error in rejecting the null hypothesis (H o = the F 2 family segregates as one resistance-endowing gene in a 0.25R:0.5 F 1 :0.25S ratio or two resistance genes in a 0.3125R: 0.625 F 1 :0.0625S ratio or three genes in a 0.375R:0.609375 F 1 :0.015625S ratio; BC segregates as 0.5 F 1 :0.5S (1 gene), 0.75 F 1 :0.25S (2 genes) or 0.875 F 1 :0.125S (3 genes) ( Table  S1). The significance level was a = 0.05 (two-sided). Statistical differences between survival proportion pairs or multiple comparisons (heterogeneity test) were also assessed by v 2 tests performed by using the statistical software R with the command prop.test.
Herbicide cross-resistance study In a separate study, a sub-sample of 100 seedlings of each parental R and S, F 1 and F 2 families was also treated with a discriminating (label) dose of the ACCase-inhibiting herbicides sethoxydim (186 g ha À1 ) or diclofop-methyl (375 g ha À1 ) or the very dissimilar ALS-inhibiting herbicides, sulfometuron (30 g ha À1 ) or chlorsulfuron (30 g ha À1 ). Survivors of chlorsulfuron and diclofop-methyl were treated with 2 kg malathion ha À1 followed eight hours later by another treatment with the same herbicide at the same dose. In L. rigidum, malathion can overcome chlorsulfuron resistance by inhibiting cytochrome P450 enzymatic activity (Christopher et al. 1994;Preston et al. 1996). Thus, malathion was used to test for a common genetic basis for cross-resistance between the very dissimilar herbicides diclofop-methyl and chlorsulfuron in R, F 1 , F 2 and S plants. We tested whether common P450 enzymes (i.e. genes) for resistance to both herbicides would be equally inhibited by the P450 inhibitor malathion. Sethoxydim and sulfometuron were used to identify the presence of resistant individuals carrying mutations in the ACCase or ALS gene, respectively (Busi et al. 2011a). Metabolism-based resistance has never been reported for both sethoxydim and sulfometuron, and thus, survival to these herbicides is likely due to occurrence of mutant alleles in ACCase or ALS genes, respectively (Christopher et al. 1991;Tardif et al. 1993).

Minimum gene number and additive-dominance model
The estimation of the minimum number of resistance genes (N E ) indicated that, at the population level, at least three resistance genes were present in the original parental herbicide-susceptible (S) L. rigidum population. Most probably, all three genes were not all initially stacked within S individuals; however, over three cycles of recurrent low-dose diclofop-methyl selection, these three resistance genes became accumulated within surviving individuals through crosspollination. Once the three resistance genes were present within individuals, then they collectively endowed relatively high-level diclofop-methyl resistance. Pooled data obtained from parental families (R and S) and respective F 1 , F 2 and BC crosses at the three diclofop-methyl doses tested were assessed by a joint scaling test. Student's t-test revealed no significant deviation from zero of quantities A and C calculated in the scaling test (Table 4). Thus, these results indicate that gene interactions are well explained by an additivedominance model at each diclofop-methyl dose used in the Mendelian segregation analysis.
Herbicide dose-response and genetic analysis of F 1 families As expected, the herbicide dose-response study confirmed the R parent as diclofop-methyl resistant and the original parental S population as susceptible (Fig. 1). The comparison between dose-response curves of paternal versus  Survivors expected in F 2 is the calculated number of plants treated multiplied by the theoretical one, two or three segregation ratios, respectively, (e.g. for one-gene model that ratio is 0.25R:0.5F1:0.25S) multiplied by the observed survival (%) in R, F1 and S at that specific dose. maternal F 1 families confirmed pollen-transmitted, nuclear inheritance of herbicide resistance traits without maternal effects (P > 0.48) (Fig. 1A,B). Overall, the diclofop-methyl resistance level observed in F 1 families was lower than the respective R parent (P < 0.01).
F 1 cloned plants were confirmed to be diclofop-methyl resistant and cloned S confirmed as susceptible (Fig. 2). Overall, F 1 plants exhibited variability in response to herbicide treatment. Using cloned plants, it was possible to show the variability evident within and between F 1 families ( Fig. 2). At the highest herbicide dose tested (H), significant heterogeneity between F 1 families was evident, with three different response patterns. Higher survival than expected was observed in F 1 family no.1 and 2, versus lower than expected survival in family no. 3 and 4 and with no significant deviations from the pooled mean survival in families no. 5 and 6 ( Table 1). This was envisaged as we expected genetic variability between parental R plants used to generate F 1 families. F 1 families no. 5 and 6 were chosen to generate F 2 plants because there was a significant segregation (P = 0.005) between plant response at treatments with L versus H dose (188 vs 1500 g diclofopmethyl ha À1 ). The other F 1 families were not studied further because either this difference in herbicide response was not significant (P = 0.79) or there were no survivors at the H dose (Fig. 2).
Herbicide dose-response and genetic analysis with F 2 and BC families F 2 families Genetic analysis indicated two different phenotypic responses in families L, M and H of F 2 no. 5 (Fig. 1C). Lack-of-fit F-test revealed no significant difference between nonlinear models fitted to L and M F 2 data (F = 1.17; P = 0.35). Similarly, the chi-square test revealed that a one-gene model best fitted the data obtained with families L and M. For these two families, any other genetic model involving a greater number of genes did not fit the data ( Table 2). The family H exhibited a significant greater level of phenotypic resistance than families L and M (Fig. 1C) (F = 12.5; P = 0.002). Consistently, the one-gene model poorly fitted the segregation data in the F 2 no. 5 family H (P < 0.02), whereas two-gene and three-gene models were found to be appropriate with no significant statistical deviations (P > 0.13) ( Table 2). F 2 no. 6 families exhibited three different levels of diclofop-methyl resistance (Fig. 1D). Family L had the lowest resistance level, although the lack-of-fit test did not reveal a significant difference between curves fitted to families L and M (F = 1.49; P = 0.24). The observed survival in family L was significantly lower than expected by the one-gene model at doses 375 and 1500 g diclofop-methyl ha À1 (Table 3). Family M showed greater level diclofop-methyl resistance and a one-gene model best fitted the observed plant survival at the three different diclofop-methyl doses (Table 3). Two-and three-gene models did not fit the data obtained with the M family (P < 0.01). Family H exhibited the highest resistance level (F = 7.47; P = 0.001) (Fig. 2D). Either a two-or a three-gene model had the best fit with the data (P > 0.08) ( Table 3).

Back-cross families
As expected, diclofop-methyl resistance was evident in two BC families tested (Fig. 1F,G). The overall resistance level observed in BC no. 5 did not significantly decrease compared with the respective F 1 no. 5 lines (P = 0.40), whereas BC no. 6 exhibited a significantly lower resistance level than F 1 s no. 6 (P < 0.001) (Fig. 1E,F). The segregation analysis suggests two genes were associated with resistance in BC family no. 5, similar to that observed in the F 2 family H. A one-gene model was found to best fit the BC family no. 6 as in the F 2 family M (

Genetic basis of herbicide cross-resistance
Importantly, recurrent selection at low diclofop-methyl dose resulted not only in diclofop-methyl resistance but in a level of cross-resistance to a dissimilar herbicide (chlorsulfuron) of different mode of action. Here, parental R, F 1 and F 2 (H type) displayed diclofop-methyl resistance and a variable degree of cross-resistance to the dissimilar ALSinhibiting herbicide chlorsulfuron (Fig. 3A). It is emphasized that before recurrent diclofop-methyl low-dose selection, the parental S population was chlorsulfuron sensitive. Therefore, the diclofop-methyl resistance resulting from recurrent low-dose diclofop-methyl selection concomitantly resulted in chlorsulfuron resistance. A common linkage between these two otherwise very dissimilar herbicides is that both are metabolized by P450 enzymes (Busi et al. 2011b) and it has been documented in L. rigidum that malathion synergises chlorsulfuron because malathion can inhibit cytochrome P450 enzymes able to metabolize chlorsulfuron (Christopher et al. 1994). Therefore, if the exact same P450 enzyme(s) metabolize both diclofopmethyl and chlorsulfuron, then malathion should be able to synergize both herbicides. This was tested by first treating with either diclofop-methyl or chlorsulfuron, observing plant survival (resistant plants) then treating with diclofopmethyl or chlorsulfuron plus malathion. As expected, plants surviving chlorsulfuron alone (57% mean plant survival) were killed by a subsequent treatment with chlorsulfuron plus malathion (7% survival). However, malathion could not synergize diclofop-methyl as survivors from diclofop-methyl treatment (83% survival) had 77% survival (v 2 = 2.03; P = 0.15) by a subsequent diclofop-methyl plus malathion treatment (Fig. 3B). This indicates that the genes endowing resistance to diclofop-methyl are likely dif-  Figure 3 Cross-resistance between diclofop-methyl and chlorsulfuron observed in parental R and F 1 , H F 2 families following treatments with the recommended dose of diclofop-methyl (375 g ha À1 ; white bars) or chlorsulfuron (30 g ha À1 ; black bars). ferent to those endowing resistance to chlorsulfuron. Thus, low-dose selection with diclofop-methyl, enriched genes that confer diclofop-methyl resistance and other genes that confer chlorsulfuron resistance.
It is important to note that the diclofop-methyl resistance and the chlorsulfuron resistance were not ACCase or ALS target-site gene based. Sulfometuron (30 g ha À1 ) resulted in 100% mortality in R, S, F 1 and F 2 lines and likewise sethoxydim (186 g ha À1 ) killed all the plants treated as they do not have a resistance-endowing mutation in the ALS or ACCase gene, respectively (Neve and Powles 2005;Yu et al. 2008).

Discussion
A highly characterized herbicide-susceptible L. rigidum population when subjected to directional selection for three generations with low (below-label) dose of the herbicide diclofop-methyl evolved a high level of phenotypic resistance ( Fig. 1) (Neve and Powles 2005). Neve and Powles (2005) hypothesized that when herbicide selection occurs at low dose (within genetic variation at putatively many loci for herbicide response in a S population) then polygenic herbicide resistance can evolve, especially in a crosspollinated genus like Lolium. By contrast, the majority of the evolved cases of herbicide resistance in L. rigidum and other weed species is endowed by genes of major effect, typically a mutation in the herbicide target-site (Délye 2005;Yu et al. 2007;Powles and Yu 2010). Monogenic responses are also common in other human-mediated selective systems (Roush and McKenzie 1987;Palumbi 2001). It is argued that selection of pre-existing genetic variation may result in more rapid evolution because it is immediately available in response to selection, whereas there would be no time to allow new beneficial mutation to arise (Gomulkiewicz et al. 2010). Similarly, the results of this study suggest the relative greater importance of standing genetic variation versus novel mutations when herbicide selection occurs at low dose on a relatively limited number of individuals and highly resistant phenotypes are very rare. It is emphasized that the L. rigidum S population used in this study was well characterized and before selection did not contain major gene(s) endowing herbicide resistance (Yu et al. 2007). Most probably, single major and minor genes can additively and/or synergistically interact to endow a level of resistance allowing plants to survive herbicide field applications. As suggested by Neve et al. (2009), it is possible that herbicide resistance evolution occurs through a 'step-by-step' pathway involving mainly polygenes in the early phase of selection before resistant phenotypes are evident in the field. A similar pathway was also suggested as genetic mechanism for evolution of mimicry in insects (Orr and Coyne 1992). Our study indicates that a few quantitative genes of small effect present in a small unselected herbicide-susceptible L. rigidum population were enriched by directional low-dose herbicide selection and these genes accumulated through crosspollination among survivors to ultimately endow a substantial level of resistance. By contrast, herbicide resistance evolution because of gene stacking seems unlikely to occur in a self-fertilized species as shown in Arabidopsis thaliana L. following seven cycles of recurrent selection at low doses of glyphosate (Brotherton et al. 2007). Genetic drift can also cause significant shift in small populations (Falconer 1981). However, in the L. rigidum population analysed here, the persistent herbicide selection to produce resistant progeny likely prevented random loss of resistance alleles and resulted in high-level phenotypic resistance to diclofop-methyl (Neve and Powles 2005). Notably, recurrent selection with low diclofop-methyl dose resulted not only in diclofop-methyl resistance but cross-resistance to the ALS herbicide chlorsulfuron, a chemically unrelated and dissimilar herbicide mode of action. This diclofop-methyl and chlorsulfuron cross-resistance was not because of target-site mutations in either the ACCase or ALS gene as also established by previous extensive molecular analysis and biochemical work (Neve and Powles 2005;Yu et al. 2007). Lolium rigidum can detoxify diclofop-methyl and chlorsulfuron through P450-mediated metabolism (Christopher et al. 1991;Preston et al. 1996;Preston and Powles 1998), and we have evidence of enhanced rates of diclofop-methyl metabolism in these low-dose-selected resistant plants (Yu et al., unpublished manuscript). However, other resistance-endowing mechanism(s) cannot be ruled out. For example, gene amplification was documented as a major force for the evolution of resistance in bacteria, insects and plants (Sandegren and Andersson 2009;Gaines et al. 2010;Bass and Field 2011).
Our study here adds to only a small number of studies that have reported polygenic control of evolved herbicide resistance in weeds under herbicide selection (Letouze and Gasquez 2001;Preston 2003;Petit et al. 2010a,b;Busi et al. 2011b). In this study, the L. rigidum R plants exhibited high-level diclofop-methyl resistance, and cloned individuals from F 1 families treated at different doses of diclofopmethyl showed genetic complexity by segregation of different resistant phenotypes (Fig. 2). We tested whether the three different levels of resistance to diclofop-methyl in cloned F 1 plants would have resulted in gene segregation in F 2 families with different levels of inherited resistance. Within each F 1 family, we confirmed that the variability in response to diclofop-methyl dose was heritable with a more resistant F 2 progeny generated from F 1 phenotypes able to survive an incrementally higher diclofop-methyl dose (Fig. 1). The observed resistance levels obtained in segregating F 2 and BC families are consistent with the hypothesis of polygenic resistance selected by enrichment (accumulation) of quantitative genes (Neve and Powles 2005;Busi and Powles 2009;Manalil et al. 2011, Manalil et al. 2012. Genetic analyses conducted in this study have clearly shown that there is involvement of more than one gene, and probably, a minimum of three effective resistance genes were enriched from the original parental S population. A population genetic analysis conducted with the polygenic modelling platform qu-gene (Podlich and Cooper 1998) also indicates polygenic resistance. A theoretical population shift was simulated by qu-gene under a range of parameters and the assumption of additive genes (see Data S1 for methodology). The results obtained by modelling simulations were compared with the observed population shift reported by Neve and Powles (2005) and similar evolutionary dynamics (observed versus simulated) were obtained under the assumption of two or three genes involved in diclofop-methyl resistance. Conversely, a greater number of resistance genes would have not resulted in similar and rapid population shifts (Figs S1-S5). Busi et al. (2011b) reported evidence of genetic linkage between P450 genes for diclofop-methyl and chlorsulfuron resistance in a field-evolved resistant L. rigidum population. Here, in the S population, diclofop-methyl recurrent selection has increased the gene frequency for resistance to both ACCase-and ALS-inhibiting herbicides (diclofop-methyl and chlorsulfuron) in the R, F 1 and F 2 lines. Thus, we hypothesize that these are P450 genes also because diclofop-methyl and chlorsulfuron are known to be P450-metabolized in wheat (Werck-Reichhart et al. 2000). However, the different plant response observed with malathion (absence or presence of synergism between malathion and diclofop-methyl or chlorsulfuron, respectively) suggests that efficient gene(s) endowing diclofop-methyl resistance may be different to gene(s) endowing chlorsulfuron resistance. This supports the hypothesis of polygenic control of resistance in this population and presence of complex gene diversity to endow resistance across different herbicide classes (Manalil et al. 2011).
Herbicide use at low dose allows the selection of all genetic trait(s) that minimize herbicide damage and thus endow plant survival. On the contrary, use of high (recommended) herbicide doses may eliminate those weakly endowing resistance gene traits. Herbicide resistance evolution in weed species is a global challenge in modern agroecosystems and increasingly evolutionary concerns are raised for the sustainability of world food production in a changing climate (Neve et al. 2009;Powles and Yu 2010;Thrall et al. 2011). It appears that relatively small L. rigidum populations have great evolutionary potential and can quickly adapt to harsh selective environments. An immediate applied conclusion from this research is that avoiding low-dose herbicide usage could be one component of helping achieve herbicide sustainability. Herbicide use at low doses and/or under conditions that result in substantial plant survival increases the risks of enrichment of minor gene traits leading to herbicide resistance. As more information on the complex genetic basis of herbicide resistance is unravelled, it becomes evident that herbicides should be used at doses that cause very high weed mortality, and with herbicide diversity and use of nonchemical weed control strategies to minimize the likelihood of resistance genes being passed to the next generation.

Data archiving statement
Data for this submission available as supplementary material hosted by Wiley Online.

Supporting Information
Additional Supporting Information may be found in the online version of this article: Figure S1. Plant survival observed at 375 g diclofop ha À1 in the susceptible VLR1 parent population and the progenies 1 (VLR1 0.1), 2 (VLR1 0.1 0.5) and 3 (VLR1 0.1 0.5 2.0) subjected to recurrent selection at the reduced doses of diclofop-methyl. Figure S2. Response to selection simulated with QU-GENE model under the hypothesis of completely additive model (Podlich and Cooper 1998). Three generations of recurrent selection, two additive genes, population size of 100 plants, initial gene frequency of 0.06 (population mean as plant survival in the unselected population VLR1 at 375 g diclofop-methyl ha À1 ), selection proportion 0.2, heritability 50% and three generations of recurrent selection. Figure S3. Response to selection simulated with QU-GENE model under the hypothesis of completely additive model (Podlich and Cooper 1998). Three generations of recurrent selection, three additive genes, population size of 100 plants, initial gene frequency of 0.06 (population mean as plant survival in the unselected population VLR1 at 375 g diclofop-methyl ha À1 ), selection proportion 0.2, heritability 50% and three generations of recurrent selection. Figure S4. Response to selection simulated with QU-GENE model under the hypothesis of completely additive model (Podlich and Cooper 1998). Three generations of recurrent selection, five additive genes, population size of 100 plants, initial gene frequency of 0.06 (population mean as plant survival in the unselected population VLR1 at 375 g diclofop-methyl ha À1 ), selection proportion 0.2, heritability 50% and three generations of recurrent selection. Figure S5. Response to selection simulated with QU-GENE model under the hypothesis of completely additive model (Podlich and Cooper 1998). Three generations of recurrent selection, 10 additive genes, population size of 100 plants, initial gene frequency of 0.06 (population mean as plant survival in the unselected population VLR1 at 375 g diclofopmethyl ha À1 ), selection proportion 0.2, heritability 50% and three generations of recurrent selection. Table S1. Genotypic segregation pattern and expected phenotypic survival in F 2 and BC families according to one-gene (A), two-gene (B) and three-gene (C) resistance models.
Data S1. Metadata sets. Data S2. Methodology. Please note: Wiley-Blackwell are not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.