Chlamydomonas reinhardtii as a model system for pro-active herbicide resistance evolution research

Modern herbicides greatly contribute to world agricultural production but their sustainability is threatened by the widespread evolution of herbicide resistant weedy plant populations. Despite the commercial and scientific importance of resistance, there has not been an experimental model system for pro-actively evaluating the potential for herbicide resistance evolution. Here, utilizing the rapidly growing, unicellular photosynthetic microalgae Chlamydomona s reinhardtii (Dangeard), a ratchet protocol has been developed that solves the problem of maintaining both large populations and strong herbicide selection. The ratchet protocol is a progressive set of cycles, each cycle commencing with a population of approximately one million individuals apportioned amongst three herbicide doses for 14 days. Whenever the evolving population demonstrates growth across the three herbicide selection intensities, then the population ratchets to the next cycle of higher herbicide dose. Therefore, by always maintaining large populations under selection pressure, this system offers the opportunity for beneficial mutations to arise and be enriched. Using the well-characterized atrazine herbicide, the ratchet protocol resulted in rapid evolution of populations with different levels of resistance. This robust laboratory based Chlamydomonas system is proposed for application in establishing the respective propensity for resistance evolution to herbicides or other selecting agents.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 257–266.     

Sexual selection expedites the evolution of pesticide resistance

The evolution of insecticide resistance by crop pests and disease vectors causes serious problems for agriculture and health. Sexual selection can accelerate or hinder adaptation to abiotic challenges in a variety of ways, but the effect of sexual selection on resistance evolution is little studied. Here, we examine this question using experimental evolution in the pest insect Tribolium castaneum. The experimental removal of sexual selection slowed the evolution of resistance in populations treated with pyrethroid pesticide, and also reduced the rate at which resistance was lost from pesticide‐free populations. These results suggest that selection arising from variance in mating and fertilization success can augment natural selection on pesticide resistance, meaning that sexual selection should be considered when designing strategies to limit the evolution of pesticide resistance.     

Evolution of resistance against powdery mildew in winter wheat populations conducted under dynamic management. II. Adult plant resistance

The evolution of adult plant resistance towards powdery mildew (caused by Blumeria graminis f. sp. tritici) was investigated in 11 wheat populations cultivated for 10 years in a French network for dynamic management (DM) of wheat genetic resources. The aims of the study were to compare the evolution of resistance in sites submitted to different powdery mildew pressure and to investigate the implication of specific resistance gene action in adult plant resistance. For this, 7 of the 11 populations were characterized for their composition of specific resistance genes (results presented in a former paper). Even though no population differed significantly from the initial PA0 pool for mean adult plant resistance, divergence appeared among the final populations. The populations with the highest adult plant resistance level originated from sites where powdery mildew pressure is known to be high (Vervins, Le Rheu), whereas populations with the lowest adult plant resistance corresponded to areas with no, or very low, powdery mildew pressure (Toulouse, Montreuil-Bellay). A residual effect of defeated specific resistance genes was hypothesized, as lines accumulating at least two specific resistance genes appeared more resistant. Additional quantitative resistance seemed to be involved in adult plant resistance. DM lines appeared then as an interesting source of variability for resistance towards powdery mildew. Moreover, as these lines had been grown in mixed populations they may be appropriate as components of a composite cultivar. Received: 15 December 1999 / Accepted: 30 December 1999     

Comparison of Parameters from Rhythmometric Models with Multiple Components on Hybrid Data

Population multiple components is a statistical tool useful for the analysis of time-dependent hybrid data. With a small number of parameters, it is possible to model and to predict the periodic behavior of a population. In this article, we propose two methods to compare among populations rhythmometric parameters obtained by multiple component analysis. The first is a parametric method based in the usual statistical techniques for comparison of mean vectors in multivariate normal populations. The method, through MANOVA analysis, allows comparison of the MESOR and amplitude-acrophase pair of each component among two or more populations. The second is a nonparametric method, based in bootstrap techniques, to compare parameters from two populations. This test allows one to compare the MESOR, the amplitude, and the acrophase of each fitted component, as well as the global amplitude, orthophase, and bathyphase estimated when all fitted components are harmonics of a fundamental period. The idea is to calculate a confidence interval for the difference of the parameters of interest. If this interval does not contain zero, it can be concluded that the parameters from the two models are different with high probability. An estimation of p-value for the corresponding test can also be calculated. Both methods are illustrated with an example, based on clinical data. The nonparametric test can also be applied to paired data, a special situation of great interest in practice. By the use of similar bootstrap techniques, we illustrate how to construct confidence intervals for any rhythmometric parameter estimated from population multiple components models, including the orthophase, bathyphase, and global amplitude. These tests for comparison of parameters among populations are a needed tool when modeling the nonsinusoidal rhythmic behavior of hybrid data by population multiple component analysis.     

Comparison of parameters from rhythmometric models with multiple components on hybrid data

Population multiple components is a statistical tool useful for the analysis of time-dependent hybrid data. With a small number of parameters, it is possible to model and to predict the periodic behavior of a population. In this article, we propose two methods to compare among populations rhythmometric parameters obtained by multiple component analysis. The first is a parametric method based in the usual statistical techniques for comparison of mean vectors in multivariate normal populations. The method, through MANOVA analysis, allows comparison of the MESOR and amplitude-acrophase pair of each component among two or more populations. The second is a nonparametric method, based in bootstrap techniques, to compare parameters from two populations. This test allows one to compare the MESOR, the amplitude, and the acrophase of each fitted component, as well as the global amplitude, orthophase, and bathyphase estimated when all fitted components are harmonics of a fundamental period. The idea is to calculate a confidence interval for the difference of the parameters of interest. If this interval does not contain zero, it can be concluded that the parameters from the two models are different with high probability. An estimation of p-value for the corresponding test can also be calculated. Both methods are illustrated with an example, based on clinical data. The nonparametric test can also be applied to paired data, a special situation of great interest in practice. By the use of similar bootstrap techniques, we illustrate how to construct confidence intervals for any rhythmometric parameter estimated from population multiple components models, including the orthophase, bathyphase, and global amplitude. These tests for comparison of parameters among populations are a needed tool when modeling the nonsinusoidal rhythmic behavior of hybrid data by population multiple component analysis.     

Residual contact vial method for the rapid on-site detection of insecticide resistance in Thrips palmi

A residual contact vial plus water (RCVpW) bioassay method, in which water was supplemented to minimize control mortality, was established to monitor insecticide resistance in field populations of the melon thrips, Thrips palmi. In the RCVpW, median lethal doses (LD₅₀) of six insecticides commonly used in T. palmi control, were determined at 8 h post-treatment, using a susceptible RDA strain according to the RCVpW protocol. Diagnostic doses for on-site resistance monitoring of the six insecticides, which were determined as doses two-fold higher than required to achieve LD₉₀ in the RDA strain, were in the range of 0.299 to 164.3 μg^?1 cm². Insecticide resistance levels in five field populations of T. palmi were evaluated to test the applicability of RCVpW in monitoring the pest. Although the RDA strain exhibited 100% mortality to diagnostic doses, field populations showed a reduced mortality in response to all test insecticides, indicating different degrees of resistance. In particular, all test field populations exhibited a significantly low mortality in response to spinosad, suggesting a wide distribution of spinosad resistance. Synergistic bioassay revealed that cytochrome P450-mediated metabolic factor is involved in spinosad resistance in the Korean population. Interestingly, an apparently reduced mortality to emamectin benzoate and chlofenapyr was observed in some field populations, perhaps suggesting uneven distribution of resistance to these insecticides in field populations. Our study showed that the RCVpW protocol can be employed both as an on-site resistance monitoring method for major thrip species, and in the selection of appropriate insecticides for their control.     

A methodology based on insecticide impregnated filter paper for monitoring resistance to deltamethrin in Triatoma infestans field populations

The domiciliary presence of Triatoma infestans (Klug) (Hemiptera: Reduviidae) after control interventions was reported in recent years. Toxicological studies showed high levels of resistance to pyrethroids suggesting resistance as one of the main causes of deficient control. The aim of the present study was to develop a protocol to test resistance to deltamethrin in T. infestans collected from the field by discriminate concentration. To evaluate field insects, the effect of age (early vs. later) and nutritional state (starved vs. fed) on the deltamethrin susceptibility of each developmental stage was studied. Topical and insecticide impregnated paper bioassays were used. Using the impregnated paper, the susceptibility to deltamethrin was not affected by the age of the stadium and the nutritional states, and varied with the post‐exposure time and with the different developmental stages. A discriminant concentration of deltamethrin (0.36% w/v) impregnated in filter paper was established for all developmental stages. Finally, the methodology and the discriminant concentration were evaluated in the laboratory showing high sensitivity in the discrimination of resistance. The present study developed a methodology of exposure to insecticide impregnated papers and proposes a protocol to test T. infestans in field populations with the aim to detect early evolution of resistance to deltamethrin.     

Weed response to herbicides: regional-scale distribution of herbicide resistance alleles in the grass weed Alopecurus myosuroides

Effective herbicide resistance management requires an assessment of the range of spatial dispersion of resistance genes among weed populations and identification of the vectors of this dispersion. In the grass weed Alopecurus myosuroides (black-grass), seven alleles of the acetyl-CoA carboxylase (ACCase) gene are known to confer herbicide resistance. Here, we assessed their respective frequencies and spatial distribution on two nested geographical scales (the whole of France and the French administrative district of C?te d'Or) by genotyping 13 151 plants originating from 243 fields. Genetic variation in ACCase was structured in local populations at both geographical scales. No spatial structure in the distribution of resistant ACCase alleles and no isolation by distance were detected at either geographical scale investigated. These data, together with ACCase sequencing and data from the literature, suggest that evolution of A. myosuroides resistance to herbicides occurred at the level of the field or group of adjacent fields by multiple, independent appearances of mutant ACCase alleles that seem to have rather restricted spatial propagation. Seed transportation by farm machinery seems the most likely vector for resistance gene dispersal in A. myosuroides.     

Testing some hypotheses on Human Evolution and sexual dimorphism

Research on human evolution and sexual dimorphism motivates an interesting test problem. In studying hominid phylogeny it is of interest to test whether parallel evolution plays a role. With regard to sexual dimorphism it is of interest to know whether the directions of sexual dimorphism in the populations being compared are the same. We show that testing these two problems gives rise to the same type of hypothesis testing, viz. the problem of testing the hypothesis that the means of independent, normally distributed random vectors with unit covariance matrices are situated on a straight line through the origin. A test is proposed and applied to study the sexual dimorphism of 20 recent skull populations. In this example the hypothesis of equal directions of sexual dimorphism is rejected. The classical theory of constructing multiple discriminant functions (canonical variates) is adapted to the problem of comparing sexual dimorphisms.     

Revisiting the role of introgression vs shared ancestral polymorphisms as key processes shaping genetic diversity in the recently separated sibling species of the Anopheles gambiae complex

The role of interspecific hybridisation in the evolution of pest species is poorly understood. In mosquito disease vectors this is of particular importance due to the evolution of insecticide resistance and the proposed release of transgenic strains that are refractory to the malaria parasite. In this study, we apply population genetic methods in a novel manner to determine whether mitochondrial DNA sequences have introgressed between the closely related African malaria vectors Anopheles gambiae and A. arabiensis. Our results suggest that speciation was geologically recent and ancestral haplotypes at the ND5 locus are retained in both species. In addition, comparing haplotype frequencies in allopatric and sympatric populations, suggest locale specific unidirectional introgression of mitochondria from A. arabiensis into A. gambiae.     

Testing some hypotheses on human evolution and sexual dimorphism

Research on human evolution and sexual dimorphism motivates an interesting test problem. In studying hominid phylogeny it is of interest to test whether parallel evolution plays a role. With regard to sexual dimorphism it is of interest to known whether the directions of sexual dimorphism in the populations being compared are the same. We show that testing these two problems gives rise to the same type of hypothesis testing, viz. the problem of testing the hypothesis that the means of independent, normally distributed random vectors with unit covariance matrices are situated on a straight line through the origin. A test is proposed and applied to study the sexual dimorphism of 20 recent skull populations. In this example the hypothesis of equal directions of sexual dimorphism is rejected. The classical theory of constructing multiple discriminant functions (canonical variates) is adapted to the problem of comparing sexual dimorphisms.     

Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so‐called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral‐like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low‐dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.     

Least and most powerful phylogenetic tests to elucidate the origin of the seed plants in the presence of conflicting signals under misspecified models

Several tests of molecular phylogenies have been proposed over the last decades, but most of them lead to strikingly different P-values. I propose that such discrepancies are principally due to different forms of null hypotheses. To support this hypothesis, two new tests are described. Both consider the composite null hypothesis that all the topologies are equidistant from the true but unknown topology. This composite hypothesis can either be reduced to the simple hypothesis at the least favorable distribution (frequentist significance test [FST]) or to the maximum likelihood topology (frequentist hypothesis test [FHT]). In both cases, the reduced null hypothesis is tested against each topology included in the analysis. The tests proposed have an information-theoretic justification, and the distribution of their test statistic is estimated by a nonparametric bootstrap, adjusting P-values for multiple comparisons. I applied the new tests to the reanalysis of two chloroplast genes, psaA and psbB, and compared the results with those of previously described tests. As expected, the FST and the FHT behaved approximately like the Shimodaira-Hasegawa test and the bootstrap, respectively. Although the tests give overconfidence in a wrong tree when an overly simple nucleotide substitution model is assumed, more complex models incorporating heterogeneity among codon positions resolve some conflicts. To further investigate the influence of the null hypothesis, a power study was conducted. Simulations showed that FST and the Shimodaira-Hasegawa test are the least powerful and FHT is the most powerful across the parameter space. Although the size of all the tests is affected by misspecification, the two new tests appear more robust against misspecification of the model of evolution and consistently supported the hypothesis that the Gnetales are nested within gymnosperms.     

Detection and evolution of resistance to the pyrethroid cypermethrin in Helicoverpa zea (Lepidoptera: Noctuidae) populations in Texas

The bollworm, Helicoverpa zea (Boddie), is a key pest of cotton in Texas. Bollworm populations are widely controlled with pyrethroid insecticides in cotton and exposed to pyrethroids in other major crops such as grain sorghum, corn, and soybeans. A statewide program that evaluated cypermethrin resistance in male bollworm populations using an adult vial test was conducted from 2003 to 2006 in the major cotton production regions of Texas. Estimated parameters from the most susceptible field population currently available (Burleson County, September 2005) were used to calculate resistance ratios and their statistical significance. Populations from several counties had statistically significant (P < or = 0.05) resistance ratios for the LC(50), indicating that bollworm-resistant populations are widespread in Texas. The highest resistance ratios for the LC(50) were observed for populations in Burleson County in 2000 and 2003, Nueces County in 2004, and Williamson and Uvalde Counties in 2005. These findings explain the observed pyrethroid control failures in various counties in Texas. Based on the assumption that resistance is caused by a single gene, the Hardy-Weinberg equilibrium formula was used for estimation of frequencies for the putative resistant allele (q) using 3 and 10 microg/vial as discriminatory dosages for susceptible and heterozygote resistant insects, respectively. The influence of migration on local levels of resistance was estimated by analysis of wind trajectories, which partially clarifies the rapid evolution of resistance to cypermethrin in bollworm populations. This approach could be used in evaluating resistance evolution in other migratory pests.     

Tailored HIV-1 Vectors for Genetic Modification of Primary Human Dendritic Cells and Monocytes

Monocyte-derived dendritic cells (MDDCs) play a key role in the regulation of the immune system and are the target of numerous gene therapy applications. The genetic modification of MDDCs is possible with human immunodeficiency virus type 1 (HIV-1)-derived lentiviral vectors (LVs) but requires high viral doses to bypass their natural resistance to viral infection, and this in turn affects their physiological properties. To date, a single viral protein is able to counter this restrictive phenotype, Vpx, a protein derived from members of the HIV-2/simian immunodeficiency virus SM lineage that counters at least two restriction factors present in myeloid cells. By tagging Vpx with a short heterologous membrane-targeting domain, we have obtained HIV-1 LVs incorporating high levels of this protein (HIV-1-Src-Vpx). These vectors efficiently transduce differentiated MDDCs and monocytes either as previously purified populations or as populations within unsorted peripheral blood mononuclear cells (PBMCs). In addition, these vectors can be efficiently pseudotyped with receptor-specific envelopes, further restricting their cellular tropism almost uniquely to MDDCs. Compared to conventional HIV-1 LVs, these novel vectors allow for an efficient genetic modification of MDDCs and, more importantly, do not cause their maturation or affect their survival, which are unwanted side effects of the transduction process. This study describes HIV-1-Src-Vpx LVs as a novel potent tool for the genetic modification of differentiated MDDCs and of circulating monocyte precursors with strong potential for a wide range of gene therapy applications.     

Reverse evolution: selection against costly resistance in disease-free microcosm populations of Paramecium caudatum

Evolutionary costs of parasite resistance arise if genes conferring resistance reduce fitness in the absence of parasites. Thus, parasite-mediated selection may lead to increased resistance and a correlated decrease in fitness, whereas relaxed parasite-mediated selection may lead to reverse evolution of increased fitness and a correlated decrease in resistance. We tested this idea in experimental populations of the protozoan Paramecium caudatum and the parasitic bacterium Holospora undulata. After eight years, resistance to infection and asexual reproduction were compared among paramecia from (1) "infected" populations, (2) uninfected "naive" populations, and (3) previously infected, parasite-free "recovered" populations. Paramecia from "infected" populations were more resistant (+12%), but had lower reproduction (-15%) than "naive" paramecia, indicating an evolutionary trade-off between resistance and fitness. Recovered populations showed similar reproduction to naive populations; however, resistance of recently (<3 years) recovered populations was similar to paramecia from infected populations, whereas longer (>3 years) recovered populations were as susceptible as naive populations. This suggests a weak, convex trade-off between resistance and fitness, allowing recovery of fitness, without complete loss of resistance, favoring the maintenance of a generalist strategy of intermediate fitness and resistance. Our results indicate that (co)evolution with parasites can leave a genetic signature in disease-free populations.     

Evolution of resistance against powdery mildew in winter wheat populations conducted under dynamic management. I – Is specific seedling resistance selected?

Dynamic management has been proposed as a complementary strategy to gene banks for the conservation of genetic resources. The evolution of frequencies of genes for specific resistance towards powdery mildew (caused by Blumeria graminis f. sp. tritici) in populations of a French network for dynamic management of bread wheat genetic resources was investigated after 10 years of multiplication without human selection. The objective was to determine whether specific resistance gene diversity was maintained in the populations and whether any changes could be attributed to selection due to pathogen pressure. Seven populations, originating from four of the network sites, were characterized and compared to the initial population for six specific resistance gene frequencies detected by nine Blumeria graminis f. sp. tritici isolates. Diversity decreased at the population level, but because of a strong differentiation between the populations, this diversity was maintained at the network level. The comparison of Fst parameters estimated on neutral markers (RFLP) and on resistance gene data revealed that in two of the populations specific resistance genes had been selected by pathogen pressure, whereas evolution in two other populations seemed to be the result of genetic drift. For the three last populations, conclusions were less clear, as one had probably experienced a strong bottleneck and the other two presented intermediate Fst values. A dynamic management network with sites contrasted for pathogen pressure, allowing genetic drift in some populations and selection in others, appeared, at least on the short term, to be a good tool for maintaining the diversity of genes for specific resistance to powdery mildew. Received: 15 December 1999 / Accepted: 30 December 1999     

Using experimental evolution to study the physiological mechanisms of desiccation resistance in Drosophila melanogaster

Data from populations undergoing experimental evolution can be used to make comparisons between physiologically differentiated populations and to determine evolutionary trajectories. Comparisons of long-established laboratory populations of Drosophila melanogaster that are strongly differentiated with respect to desiccation resistance are used to test alternative hypotheses concerning the mechanisms that fruit flies use to survive bouts of extreme desiccation. This comparative study supports the hypothesis that, in at least one case, D. melanogaster can evolve increased resistance to desiccation by decreasing water loss rates and by increasing bulk water content but not by increasing metabolic water content or dehydration tolerance. While glycogen was involved in water storage, its primary role was in water binding, not the production of metabolic water. Measurement of the trajectories of these component mechanisms during selection for desiccation resistance is used to demonstrate that water loss rate quickly plateaus in response to selection, while water content continues to improve. This disparity reveals the value of studying evolutionary trajectories and the need for longer-term selection studies in evolutionary physiology.     

Immune selection suppresses the emergence of drug resistance in malaria parasites but facilitates its spread

Although drug resistance in Plasmodium falciparum typically evolves in regions of low transmission, resistance spreads readily following introduction to regions with a heavier disease burden. This suggests that the origin and the spread of resistance are governed by different processes, and that high transmission intensity specifically impedes the origin. Factors associated with high transmission, such as highly immune hosts and competition within genetically diverse infections, are associated with suppression of resistant lineages within hosts. However, interactions between these factors have rarely been investigated and the specific relationship between adaptive immunity and selection for resistance has not been explored. Here, we developed a multiscale, agent-based model of Plasmodium parasites, hosts, and vectors to examine how host and parasite dynamics shape the evolution of resistance in populations with different transmission intensities. We found that selection for antigenic novelty (“immune selection”) suppressed the evolution of resistance in high transmission settings. We show that high levels of population immunity increased the strength of immune selection relative to selection for resistance. As a result, immune selection delayed the evolution of resistance in high transmission populations by allowing novel, sensitive lineages to remain in circulation at the expense of the spread of a resistant lineage.In contrast, in low transmission settings, we observed that resistant strains were able to sweep to high population prevalence without interference. Additionally, we found that the relationship between immune selection and resistance changed when resistance was widespread. Once resistance was common enough to be found on many antigenic backgrounds, immune selection stably maintained resistant parasites in the population by allowing them to proliferate, even in untreated hosts, when resistance was linked to a novel epitope. Our results suggest that immune selection plays a role in the global pattern of resistance evolution.     

A TEST OF THE ROLE OF EPISTASIS IN DIVERGENCE UNDER UNIFORM SELECTION

Five populations of Drosophila melanogaster have previously been shown to be replicably different in their responses to artificial selection for knockdown resistance to ethanol fumes (Cohan and Hoffmann, 1986). The present study tests whether this divergence could be attributed to the epistatic mechanism assumed by Wright's shifting-balance model of evolution, in which alleles favored in the genetic background of one population are not favored in that of another. If this were the mechanism of divergence, crosses between selected lines from different populations would be expected to yield an epistatic loss of the selected phenotype. However, all such crosses showed a good fit to an additive model with dominance. Divergence by an epistatic mechanism may also be associated with epistatic variance within populations, but no evidence for such epistasis was found. The populations therefore appear to have responded in different ways to selection not because of epistasis but because knockdown-resistance alleles that were common in some populations were absent (or at least less common) in others.