Genome-wide detection of positive and balancing signatures of selection shared by four domesticated rainbow trout populations (Oncorhynchus mykiss).

BACKGROUND: Evolutionary processes leave footprints along the genome over time. Highly homozygous regions may correspond to positive selection of favorable alleles, while maintenance of heterozygous regions may be due to balancing selection phenomena. We analyzed data from 176 fish from four disconnected domestic rainbow trout populations that were genotyped using a high-density Axiom Trout genotyping 665K single nucleotide polymorphism array, including 20 from the US and 156 from three French lines. Using methods based on runs of homozygosity and extended haplotype homozygosity, we detected signatures of selection in these four populations. RESULTS: Nine genomic regions that included 253 genes were identified as being under positive selection in all four populations Most were located on chromosome 2 but also on chromosomes 12, 15, 16, and 20. In addition, four heterozygous regions that contain 29 genes that are putatively under balancing selection were also shared by the four populations. These were located on chromosomes 10, 13, and 19. Regardless of the homozygous or heterozygous nature of the regions, in each region, we detected several genes that are highly conserved among vertebrates due to their critical roles in cellular and nuclear organization, embryonic development, or immunity. We identified new candidate genes involved in rainbow trout fitness, as well as 17 genes that were previously identified to be under positive selection, 10 of which in other fishes (auts2, atp1b3, zp4, znf135, igf-1α, brd2, col9a2, mrap2, pbx1, and emilin-3). CONCLUSIONS: Using material from disconnected populations of different origins allowed us to draw a genome-wide map of signatures of positive selection that are shared between these rainbow trout populations, and to identify several regions that are putatively under balancing selection. These results provide a valuable resource for future investigations of the dynamics of genetic diversity and genome evolution during domestication.

Reintroducing genetic diversity in populations from cryopreserved material: the case of Abondance, a French local dairy cattle breed.

BACKGROUND: Genetic diversity is a necessary condition for populations to evolve under natural adaptation, artificial selection, or both. However, genetic diversity is often threatened, in particular in domestic animal populations where artificial selection, genetic drift and inbreeding are strong. In this context, cryopreserved genetic resources are a promising option to reintroduce lost variants and to limit inbreeding. However, while the use of ancient genetic resources is more common in plant breeding, it is less documented in animals due to a longer generation interval, making it difficult to fill the gap in performance due to continuous selection. This study investigates one of the only concrete cases available in animals, for which cryopreserved semen from a bull born in 1977 in a lost lineage was introduced into the breeding scheme of a French local dairy cattle breed, the Abondance breed, more than 20 years later. RESULTS: We found that this re-introduced bull was genetically distinct with respect to the current population and thus allowed part of the genetic diversity lost over time to be restored. The expected negative gap in milk production due to continuous selection was absorbed in a few years by preferential mating with elite cows. Moreover, the re-use of this bull more than two decades later did not increase the level of inbreeding, and even tended to reduce it by avoiding mating with relatives. Finally, the reintroduction of a bull from a lost lineage in the breeding scheme allowed for improved performance for reproductive abilities, a trait that was less subject to selection in the past. CONCLUSIONS: The use of cryopreserved material is an efficient way to manage the genetic diversity of an animal population, by mitigating the effects of both inbreeding and strong selection. However, attention should be paid to mating of animals to limit the disadvantages associated with incorporating original genetic material, notably a discrepancy in the breeding values for selected traits or an increase in inbreeding. Therefore, careful characterization of the genetic resources available in cryobanks could help to ensure the sustainable management of populations, in particular local or small populations. These results could also be transferred to the conservation of wild threatened populations.

Managing genetic diversity in breeding programs of small populations: the case of French local chicken breeds.

BACKGROUND: On-going climate change will drastically modify agriculture in the future, with a need for more sustainable systems, in particular regarding animal production. In this context, genetic diversity is a key factor for adaptation to new conditions: local breeds likely harbor unique adaptive features and represent a key component of diversity to reach resilience. However, local breeds often suffer from small population sizes, which puts these valuable resources at risk of extinction. In chickens, population management programs were initiated a few decades ago in France, relying on a particular niche market that aims at promoting and protecting local breeds. We conducted a unique comprehensive study of 22 French local breeds, along with four commercial lines, to evaluate their genetic conservation status and the efficiency of the population management programs. RESULTS: Using a 57K single nucleotide polymorphism (SNP) chip, we demonstrated that both the between- and within-breed genetic diversity levels are high in the French local chicken populations. Diversity is mainly structured according to the breeds' selection and history. Nevertheless, we observed a prominent sub-structuring of breeds according to farmers' practices in terms of exchange, leading to more or less isolated flocks. By analysing demographic parameters and molecular information, we showed that consistent management programs are efficient in conserving genetic diversity, since breeds that integrated such programs earlier had older inbreeding. CONCLUSIONS: Management programs of French local chicken breeds have maintained their genetic diversity at a good level. We recommend that future programs sample as many individuals as possible, with emphasis on both males and females from the start, and focus on a quick and strong increase of population size while conserving as many families as possible. We also stress the usefulness of molecular tools to monitor small populations for which pedigrees are not always available. Finally, the breed appears to be an appropriate operational unit for the conservation of genetic diversity, even for local breeds, for which varieties, if present, could also be taken into account.

Small freshwater ecosystems with dissimilar microbial communities exhibit similar temporal patterns.

Despite small freshwater ecosystems being biodiversity reservoirs and contributing significantly to greenhouse fluxes, their microbial communities remain largely understudied. Yet, microorganisms intervene in biogeochemical cycling and impact water quality. Because of their small size, these ecosystems are in principle more sensitive to disturbances, seasonal variation and pluri-annual climate change. However, how microbial community composition varies over space and time, and whether archaeal, bacterial and microbial eukaryote communities behave similarly remain unanswered. Here, we aim to unravel the composition and intra/interannual temporal dynamic patterns for archaea, bacteria and microbial eukaryotes in a set of small freshwater ecosystems. We monitored archaeal and bacterial community composition during 24 consecutive months in four ponds and one brook from northwestern France by 16S rRNA gene amplicon sequencing (microbial eukaryotes were previously investigated for the same systems). Unexpectedly for oxic environments, bacterial Candidate Phyla Radiation (CPR) were highly diverse and locally abundant. Our results suggest that microbial community structure is mainly driven by environmental conditions acting over space (ecosystems) and time (seasons). A low proportion of operational taxonomic units (OTUs) (<1%) was shared by the five ecosystems despite their geographical proximity (2-9 km away), making microbial communities almost unique in each ecosystem and highlighting the strong selective influence of local environmental conditions. Marked and similar seasonality patterns were observed for archaea, bacteria and microbial eukaryotes in all ecosystems despite strong turnovers of rare OTUs. Over the 2-year survey, microbial community composition varied despite relatively stable environmental parameters. This suggests that biotic associations play an important role in interannual community assembly.

Intensified Use of Reproductive Technologies and Reduced Dimensions of Breeding Schemes Put Genetic Diversity at Risk in Dairy Cattle Breeds.

In the management of dairy cattle breeds, two recent trends have arisen that pose potential threats to genetic diversity: the use of reproductive technologies (RT) and a reduction in the number of bulls in breeding schemes. The expected outcome of these changes, in terms of both genetic gain and genetic diversity, is not trivial to predict. Here, we simulated 15 breeding schemes similar to those carried out in large French dairy cattle breeds; breeding schemes differed with respect to their dimensions, the intensity of RT use, and the type of RT involved. We found that intensive use of RT resulted in improved genetic gain, but deteriorated genetic diversity. Specifically, a reduction in the interval between generations through the use of ovum pick-up and in vitro fertilization (OPU-IVF) resulted in a large increase in the inbreeding rate both per year and per generation, suggesting that OPU-IVF could have severe adverse effects on genetic diversity. To achieve a given level of genetic gain, the scenarios that best maintained genetic diversity were those with a higher number of sires/bulls and a medium intensity of RT use or those with a higher number of female donors to compensate for the increased intensity of RT.

The impact of genomic selection on genetic diversity and genetic gain in three French dairy cattle breeds.

BACKGROUND: In France, implementation of genomic evaluations in dairy cattle breeds started in 2009 and this has modified the breeding schemes drastically. In this context, the goal of our study was to understand the impact of genomic selection on the genetic diversity of bulls from three French dairy cattle breeds born between 2005 and 2015 (Montbéliarde, Normande and Holstein) and the factors that are involved. METHODS: We compared annual genetic gains, inbreeding rates based on runs of homozygosity (ROH) and pedigree data, and mean ROH length within breeds, before and after the implementation of genomic selection. RESULTS: Genomic selection induced an increase in mean annual genetic gains of 50, 71 and 33% for Montbéliarde, Normande and Holstein bulls, respectively, and in parallel, the generation intervals were reduced by a factor of 1.7, 1.9 and 2, respectively. We found no significant change in inbreeding rate for the two national breeds, Montbéliarde and Normande, and a significant increase in inbreeding rate for the Holstein international breed, which is now as high as 0.55% per year based on ROH and 0.49% per year based on pedigree data (equivalent to a rate of 1.36 and 1.39% per generation, respectively). The mean ROH length was longer for bulls from the Holstein breed than for those from the other two breeds. CONCLUSIONS: With the implementation of genomic selection, the annual genetic gain increased for bulls from the three major French dairy cattle breeds. At the same time, the annual loss of genetic diversity increased for Holstein bulls, possibly because of the massive use of a few elite bulls in this breed, but not for Montbéliarde and Normande bulls. The increase in mean ROH length in Holstein may reflect the occurrence of recent inbreeding. New strategies in breeding schemes, such as female donor stations and embryo transfer, and recent implementation of genomic evaluations in small regional breeds should be studied carefully in order to ensure the sustainability of breeding schemes in the future.

Revealing the selection history of adaptive loci using genome-wide scans for selection: an example from domestic sheep.

BACKGROUND: One of the approaches to detect genetics variants affecting fitness traits is to identify their surrounding genomic signatures of past selection. With established methods for detecting selection signatures and the current and future availability of large datasets, such studies should have the power to not only detect these signatures but also to infer their selective histories. Domesticated animals offer a powerful model for these approaches as they adapted rapidly to environmental and human-mediated constraints in a relatively short time. We investigated this question by studying a large dataset of 542 individuals from 27 domestic sheep populations raised in France, genotyped for more than 500,000 SNPs. RESULTS: Population structure analysis revealed that this set of populations harbour a large part of European sheep diversity in a small geographical area, offering a powerful model for the study of adaptation. Identification of extreme SNP and haplotype frequency differences between populations listed 126 genomic regions likely affected by selection. These signatures revealed selection at loci commonly identified as selection targets in many species ("selection hotspots") including ABCG2, LCORL/NCAPG, MSTN, and coat colour genes such as ASIP, MC1R, MITF, and TYRP1. For one of these regions (ABCG2, LCORL/NCAPG), we could propose a historical scenario leading to the introgression of an adaptive allele into a new genetic background. Among selection signatures, we found clear evidence for parallel selection events in different genetic backgrounds, most likely for different mutations. We confirmed this allelic heterogeneity in one case by resequencing the MC1R gene in three black-faced breeds. CONCLUSIONS: Our study illustrates how dense genetic data in multiple populations allows the deciphering of evolutionary history of populations and of their adaptive mutations.

Which Individuals To Choose To Update the Reference Population? Minimizing the Loss of Genetic Diversity in Animal Genomic Selection Programs.

Genomic selection (GS) is commonly used in livestock and increasingly in plant breeding. Relying on phenotypes and genotypes of a reference population, GS allows performance prediction for young individuals having only genotypes. This is expected to achieve fast high genetic gain but with a potential loss of genetic diversity. Existing methods to conserve genetic diversity depend mostly on the choice of the breeding individuals. In this study, we propose a modification of the reference population composition to mitigate diversity loss. Since the high cost of phenotyping is the limiting factor for GS, our findings are of major economic interest. This study aims to answer the following questions: how would decisions on the reference population affect the breeding population, and how to best select individuals to update the reference population and balance maximizing genetic gain and minimizing loss of genetic diversity? We investigated three updating strategies for the reference population: random, truncation, and optimal contribution (OC) strategies. OC maximizes genetic merit for a fixed loss of genetic diversity. A French Montbéliarde dairy cattle population with 50K SNP chip genotypes and simulations over 10 generations were used to compare these different strategies using milk production as the trait of interest. Candidates were selected to update the reference population. Prediction bias and both genetic merit and diversity were measured. Changes in the reference population composition slightly affected the breeding population. Optimal contribution strategy appeared to be an acceptable compromise to maintain both genetic gain and diversity in the reference and the breeding populations.

Individual and mixture acute toxicity of model pesticides chlordecone and pyriproxyfen in the estuarine copepod Eurytemora affinis.

Due to the increase in the use of phytosanitary products during the last few decades, the importance to study the effect of pesticide mixtures has been established. In this study, we investigated the acute toxicity of two model insecticides, chlordecone (CLD) and pyriproxyfen (PXF), alone and in mixtures, in the estuarine copepod Eurytemora affinis. After 48 h of exposure, the relative LC50 were 73.24 and 131.61 µg/L for PXF and CLD, respectively. The lower concentration tested (10 µg/L) did not affect the mortality of E. affinis whatever the considered chemical compound. To understand the interaction between compounds in mixture, the results were fitted to the concentration addition, Vølund, and Hewlett models. The best fit was obtained with the Hewlett model, suggesting a synergistic effect of the mixture.

Resilience of Freshwater Communities of Small Microbial Eukaryotes Undergoing Severe Drought Events.

Small and shallow aquatic ecosystems such as ponds and streams constitute a significant proportion of continental surface waters, especially in temperate zones. In comparison with bigger lakes and rivers, they harbor higher biodiversity but they also exhibit reduced buffering capacity face to environmental shifts, such that climate global change can affect them in a more drastic way. For instance, many temperate areas are predicted to undergo droughts with increasing frequency in the near future, which may lead to the temporal desiccation of streams and ponds. In this work, we monitored temporal dynamics of planktonic communities of microbial eukaryotes (cell size range: 0.2-5 µm) in one brook and one pond that experienced recurrent droughts from 1 to 5 consecutive months during a temporal survey carried out monthly for 2 years based on high-throughput 18S rDNA metabarcoding. During drought-induced desiccation events, protist communities present in the remaining dry sediment, though highly diverse, differed radically from their planktonic counterparts. However, after water refill, the aquatic protist assemblages recovered their original structure within a month. This rapid recovery indicates that these eukaryotic communities are resilient to droughts, most likely via the entrance in dormancy. This property is essential for the long-term survival and functional stability of small freshwater ecosystems.

Marked seasonality and high spatial variability of protist communities in shallow freshwater systems.

Small eukaryotes have key roles in aquatic ecosystems, influencing their local environment, global biogeochemical cycles and climate. Their impact depends on community structure, which varies along time. However, very few studies take into account temporal variation. This is especially true for small, shallow freshwater systems, which remain largely understudied despite their wide variety, global surface and intense microbial activity. We have monthly followed changes in the community structure of small microbial eukaryotes (0.2-5 µm cell diameter) for 2 years in four ponds and one brook located in North-Western France based on massive 18S rDNA amplicon 454 pyrosequencing. We detected a total of 3742 stringently defined operational taxonomic units (OTUs) encompassing all recognized eukaryotic supergroups and lineages of uncertain affiliation. Although geographically close, protist communities in the five ecosystems were contrasting, with very few shared OTUs, suggesting that environmental selection mainly drives community structure. The temporal dynamics of different high-rank taxa appeared complex and rapid at monthly scales. Despite this, a clear and reproducible seasonality was observed. As expected, low-abundance OTUs dominated the community. Although some of them appeared sporadically or remained at low frequencies during the survey, others occasionally reached relatively high abundances, sometimes recurrently. This shows that at least a fraction of low-abundance eukaryotes constitutes a seed bank. The annual proportion of primary producers, free-living heterotrophs and parasites appeared remarkably constant among the different ecosystems, suggesting underlying trends of ecosystem carrying capacity for these functional groups.

Efficiently tracking selection in a multiparental population: the case of earliness in wheat.

Multiparental populations are innovative tools for fine mapping large numbers of loci. Here we explored the application of a wheat Multiparent Advanced Generation Inter-Cross (MAGIC) population for QTL mapping. This population was created by 12 generations of free recombination among 60 founder lines, following modification of the mating system from strict selfing to strict outcrossing using the ms1b nuclear male sterility gene. Available parents and a subset of 380 SSD lines of the resulting MAGIC population were phenotyped for earliness and genotyped with the 9K i-Select SNP array and additional markers in candidate genes controlling heading date. We demonstrated that 12 generations of strict outcrossing rapidly and drastically reduced linkage disequilibrium to very low levels even at short map distances and also greatly reduced the population structure exhibited among the parents. We developed a Bayesian method, based on allelic frequency, to estimate the contribution of each parent in the evolved population. To detect loci under selection and estimate selective pressure, we also developed a new method comparing shifts in allelic frequency between the initial and the evolved populations due to both selection and genetic drift with expectations under drift only. This evolutionary approach allowed us to identify 26 genomic areas under selection. Using association tests between flowering time and polymorphisms, 6 of these genomic areas appeared to carry flowering time QTL, 1 of which corresponds to Ppd-D1, a major gene involved in the photoperiod sensitivity. Frequency shifts at 4 of 6 areas were consistent with earlier flowering of the evolved population relative to the initial population. The use of this new outcrossing wheat population, mixing numerous initial parental lines through multiple generations of panmixia, is discussed in terms of power to detect genes under selection and association mapping. Furthermore we provide new statistical methods for use in future analyses of multiparental populations.

Complex communities of small protists and unexpected occurrence of typical marine lineages in shallow freshwater systems.

Although inland water bodies are more heterogeneous and sensitive to environmental variation than oceans, the diversity of small protists in these ecosystems is much less well known. Some molecular surveys of lakes exist, but little information is available from smaller, shallower and often ephemeral freshwater systems, despite their global distribution and ecological importance. We carried out a comparative study based on massive pyrosequencing of amplified 18S rRNA gene fragments of protists in the 0.2-5 µm size range in one brook and four shallow ponds located in the Natural Regional Park of the Chevreuse Valley, France. Our study revealed a wide diversity of small protists, with 812 stringently defined operational taxonomic units (OTUs) belonging to the recognized eukaryotic supergroups (SAR--Stramenopiles, Alveolata, Rhizaria--Archaeplastida, Excavata, Amoebozoa, Opisthokonta) and to groups of unresolved phylogenetic position (Cryptophyta, Haptophyta, Centrohelida, Katablepharida, Telonemida, Apusozoa). Some OTUs represented deep-branching lineages (Cryptomycota, Aphelida, Colpodellida, Tremulida, clade-10 Cercozoa, HAP-1 Haptophyta). We identified several lineages previously thought to be marine including, in addition to MAST-2 and MAST-12, already detected in freshwater, MAST-3 and possibly MAST-6. Protist community structures were different in the five ecosystems. These differences did not correlate with geographical distances, but seemed to be influenced by environmental parameters.

Toxicity of sediment-bound pollutants in the Seine estuary, France, using a Eurytemora affinis larval bioassay.

Coastal urbanisation exposes surrounding estuarine environments to urban-related contaminants such as polycyclic aromatic hydrocarbons (PAHs), polychlorobiphenyls (PCBs) and pesticide mixtures. Hydrophobic contaminants can adsorb on estuarine sediments. They can subsequently be released on a massive scale in the aquatic environment due to artificial or natural phenomena (e.g. dredging, tides), thereby threatening living organisms. The contamination of sediment is a significant ecological issue in the Seine estuary, France. However, few relevant methods have been developed to assess sediment toxicity and its ecological impacts in a cost-effective way. In this context, we aimed to assess the toxicity of natural sediments from the Seine estuary on the development of the calanoid copepod Eurytemora affinis using a previously developed larval bioassay. This assay involves direct exposure of nauplii to elutriates of sediments for six days. Sediments were collected along the Seine estuary from six polluted sites and one reference site. Pollutants in this estuary included PAHs, PCBs and OCPs (organochlorine pesticides). Nauplius survival was significantly more affected by exposure to all contaminated sediment elutriates, than by exposure to sediment from Yville-sur-Seine (the reference site), whereas nauplius growth was significantly reduced after exposure to contaminated sediment elutriates from four of the six contaminated sites. We identified two distinct site clusters, one including both the sand-rich and the least polluted sediments (Oissel, Quillebeuf-sur-Seine, Caudebec-en-Caux) and the other including both the clay- and silt-rich, and the most polluted sediments (La Bouille, Poses, Pont de Normandie). As expected, survival was significantly more impacted after exposure to elutriates from the second cluster than from the first. This work enables (i) assessment of the toxicity of natural sediments in the Seine estuary and (ii) validation of the larval bioassay previously developed using sorbed sediment with model molecules.

Variability of individual genetic load: consequences for the detection of inbreeding depression.

Inbreeding depression is a key factor affecting the persistence of natural populations, particularly when they are fragmented. In species with mixed mating systems, inbreeding depression can be estimated at the population level by regressing the average progeny fitness by the selfing rate of their mothers. We applied this method using simulated populations to investigate how population genetic parameters can affect the detection power of inbreeding depression. We simulated individual selfing rates and genetic loads from which we computed fitness values. The regression method yielded high statistical power, inbreeding depression being detected as significant (5 % level) in 92 % of the simulations. High individual variation in selfing rate and high mean genetic load led to better detection of inbreeding depression while high among-individual variation in genetic load made it more difficult to detect inbreeding depression. For a constant sampling effort, increasing the number of progenies while decreasing the number of individuals per progeny enhanced the detection power of inbreeding depression. We discuss the implication of among-mother variability of genetic load and selfing rate on inbreeding depression studies.

Sunlight-exposed biofilm microbial communities are naturally resistant to chernobyl ionizing-radiation levels.

BACKGROUND: The Chernobyl accident represents a long-term experiment on the effects of exposure to ionizing radiation at the ecosystem level. Though studies of these effects on plants and animals are abundant, the study of how Chernobyl radiation levels affect prokaryotic and eukaryotic microbial communities is practically non-existent, except for a few reports on human pathogens or soil microorganisms. Environments enduring extreme desiccation and UV radiation, such as sunlight exposed biofilms could in principle select for organisms highly resistant to ionizing radiation as well. METHODOLOGY/PRINCIPAL FINDINGS: To test this hypothesis, we explored the diversity of microorganisms belonging to the three domains of life by cultivation-independent approaches in biofilms developing on concrete walls or pillars in the Chernobyl area exposed to different levels of radiation, and we compared them with a similar biofilm from a non-irradiated site in Northern Ireland. Actinobacteria, Alphaproteobacteria, Bacteroidetes, Acidobacteria and Deinococcales were the most consistently detected bacterial groups, whereas green algae (Chlorophyta) and ascomycete fungi (Ascomycota) dominated within the eukaryotes. Close relatives to the most radio-resistant organisms known, including Rubrobacter species, Deinococcales and melanized ascomycete fungi were always detected. The diversity of bacteria and eukaryotes found in the most highly irradiated samples was comparable to that of less irradiated Chernobyl sites and Northern Ireland. However, the study of mutation frequencies in non-coding ITS regions versus SSU rRNA genes in members of a same actinobacterial operational taxonomic unit (OTU) present in Chernobyl samples and Northern Ireland showed a positive correlation between increased radiation and mutation rates. CONCLUSIONS/SIGNIFICANCE: Our results show that biofilm microbial communities in the most irradiated samples are comparable to non-irradiated samples in terms of general diversity patterns, despite increased mutation levels at the single-OTU level. Therefore, biofilm communities growing in sunlight exposed substrates are capable of coping with increased mutation rates and appear pre-adapted to levels of ionizing radiation in Chernobyl due to their natural adaptation to periodical desiccation and ambient UV radiation.

Anthropogenic effects on population genetics of phytophagous insects associated with domesticated plants.

The hypothesis of isolation by distance (IBD) predicts that genetic differentiation between populations increases with geographic distance. However, gene flow is governed by numerous factors and the correlation between genetic differentiation and geographic distance is never simply linear. In this study, we analyze the interaction between the effects of geographic distance and of wild or domesticated status of the host plant on genetic differentiation in the bean beetle Acanthoscelides obvelatus. Geographic distance explained most of the among-population genetic differentiation. However, IBD varied depending on the kind of population pairs for which the correlation between genetic differentiation and geographic distance was examined. Whereas pairs of beetle populations associated with wild beans showed significant IBD (P < 10(-4)), no IBD was found when pairs of beetle populations on domesticated beans were examined (P= 0.2992). This latter result can be explained by long-distance migrations of beetles on domesticated plants resulting from human exchanges of bean seeds. Beetle populations associated with wild beans were also significantly more likely than those on domesticated plants to contain rare alleles. However, at the population level, beetles on cultivated beans were similar in allelic richness to those on wild beans. This similarity in allelic richness combined with differences in other aspects of the genetic diversity (i.e., IBD, allelic diversity) is compatible with strongly contrasting effects of migration and drift. This novel indirect effect of human actions on gene flow of a serious pest of a domesticated plant has important implications for the spread of new adaptations such as resistance to pesticides.