Selective Sweeps in a Nutshell: The Genomic Footprint of Rapid Insecticide Resistance Evolution in the Almond Agroecosystem

Evolution of parasites in the Anthropocene: new pressures, new adaptive directions

The Anthropocene is seeing the human footprint rapidly spreading to all of Earth's ecosystems. The fast-changing biotic and abiotic conditions experienced by all organisms are exerting new and strong selective pressures, and there is a growing list of examples of human-induced evolution in response to anthropogenic impacts. No organism is exempt from these novel selective pressures. Here, we synthesise current knowledge on human-induced evolution in eukaryotic parasites of animals, and present a multidisciplinary framework for its study and monitoring. Parasites generally have short generation times and huge fecundity, features that predispose them for rapid evolution. We begin by reviewing evidence that parasites often have substantial standing genetic variation, and examples of their rapid evolution both under conditions of livestock production and in serial passage experiments. We then present a two-step conceptual overview of the causal chain linking anthropogenic impacts to parasite evolution. First, we review the major anthropogenic factors impacting parasites, and identify the selective pressures they exert on parasites through increased mortality of either infective stages or adult parasites, or through changes in host density, quality or immunity. Second, we discuss what new phenotypic traits are likely to be favoured by the new selective pressures resulting from altered parasite mortality or host changes; we focus mostly on parasite virulence and basic life-history traits, as these most directly influence the transmission success of parasites and the pathology they induce. To illustrate the kinds of evolutionary changes in parasites anticipated in the Anthropocene, we present a few scenarios, either already documented or hypothetical but plausible, involving parasite taxa in livestock, aquaculture and natural systems. Finally, we offer several approaches for investigations and real-time monitoring of rapid, human-induced evolution in parasites, ranging from controlled experiments to the use of state-of-the-art genomic tools. The implications of fast-evolving parasites in the Anthropocene for disease emergence and the dynamics of infections in domestic animals and wildlife are concerning. Broader recognition that it is not only the conditions for parasite transmission that are changing, but the parasites themselves, is needed to meet better the challenges ahead.

Demographic fluctuations and selection during host-parasite co-evolution interactively increase genetic diversity

Host-parasite interactions can cause strong demographic fluctuations accompanied by selective sweeps of resistance/infectivity alleles. Both demographic bottlenecks and frequent sweeps are expected to reduce the amount of segregating genetic variation and therefore might constrain adaptation during co-evolution. Recent studies, however, suggest that the interaction of demographic and selective processes is a key component of co-evolutionary dynamics and may rather positively affect levels of genetic diversity available for adaptation. Here, we provide direct experimental testing of this hypothesis by disentangling the effects of demography, selection and their interaction in an experimental host-parasite system. We grew 12 populations of a unicellular, asexually reproducing algae (Chlorella variabilis) that experienced either growth followed by constant population sizes (three populations), demographic fluctuations (three populations), selection induced by exposure to a virus (three populations), or demographic fluctuations together with virus-induced selection (three populations). After 50 days (~50 generations), we conducted whole-genome sequencing of each algal host population. We observed more genetic diversity in populations that jointly experienced selection and demographic fluctuations than in populations where these processes were experimentally separated. In addition, in those three populations that jointly experienced selection and demographic fluctuations, experimentally measured diversity exceeds expected values of diversity that account for the cultures' population sizes. Our results suggest that eco-evolutionary feedbacks can positively affect genetic diversity and provide the necessary empirical measures to guide further improvements of theoretical models of adaptation during host-parasite co-evolution.

Evolutionary origin and distribution of amino acid mutations associated with resistance to sodium channel modulators in onion thrips, Thrips tabaci

In onion thrips Thrips tabaci, reduced sensitivity of the sodium channel caused by several sodium channel mutations have been correlated with pyrethroid resistance. For this study, using mitochondrial cytochrome c oxidase subunit I gene sequences, we examined the phylogenetic relation among a total of 52 thelytokous and arrhenotokous strains with different genotypes of the sodium channel mutations. Then, we used flow cytometry to estimate their ploidy. Results showed that the strains are divisible into three groups: diploid thelytoky, triploid thelytoky, and diploid arrhenotoky. Using 23 whole genome resequencing data obtained from 20 strains out of 52, we examined their genetic relation further using principal component analysis, admixture analysis, and a fixation index. Results showed that diploid and triploid thelytokous groups are further classifiable into two based on the sodium channel mutations harbored by the respective group members (strains). The greatest genetic divergence was observed between thelytokous and arrhenotokous groups with a pair of T929I and K1774N. Nevertheless, they shared a genomic region with virtually no polymorphism around the sodium channel gene loci, suggesting a hard selective sweep. Based on these findings, we discuss the evolutionary origin and distribution of the sodium channel mutations in T. tabaci.

Current stewardship practices in invasion biology limit the value and secondary use of genomic data

Invasive species threaten native biota, putting fragile ecosystems at risk and having a large-scale impact on primary industries. Growing trade networks and the popularity of personal travel make incursions a more frequent risk, one only compounded by global climate change. With increasing publication of whole-genome sequences lies an opportunity for cross-species assessment of invasive potential. However, the degree to which published sequences are accompanied by satisfactory spatiotemporal data is unclear. We assessed the metadata associated with 199 whole-genome assemblies of 89 invasive terrestrial invertebrate species and found that only 38% of these were derived from field-collected samples. Seventy-six assemblies (38%) reported an 'undescribed' sample origin and, while further examination of associated literature closed this gap to 23.6%, an absence of spatial data remained for 47 of the total assemblies. Of the 76 assemblies that were ultimately determined to be field-collected, associated metadata relevant for invasion studies was predominantly lacking: only 35% (27 assemblies) provided granular location data, and 33% (n = 25) lacked sufficient collection date information. Our results support recent calls for standardized metadata in genome sequencing data submissions, highlighting the impact of missing metadata on current research in invasion biology (and likely other fields). Notably, large-scale consortia tended to provide the most complete metadata submissions in our analysis-such cross-institutional collaborations can foster a culture of increased adherence to improved metadata submission standards and a standard of metadata stewardship that enables reuse of genomes in invasion science.

A chromosome scale assembly of the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), genome

OBJECTIVE

The tarnished plant bug (TPB), Lygus lineolaris (Palisot de Beauvois) (Hemiptera: Miridae), is a pest damaging many cultivated crops in North America. Although partial transcriptome data are available for this pest, a genome assembly was not available for this species. This assembly of a high-quality chromosome-length genome of TPB is aimed to develop the genetic resources that can provide the foundation required for advancing research on this species.

RESULTS

The initial genome of TPB assembled with paired-end nucleotide sequences generated with Illumina technology was scaffolded with Illumina HiseqX reads generated from a proximity ligated (HiC) library to obtain a high-quality genome assembly. The final assembly contained 3963 scaffolds longer than 1 kbp to yield a genome of 599.96 Mbp. The N50 of the TPB genome assembly was 35.64 Mbp and 98.68% of the genome was assembled into 17 scaffolds larger than 1 Mbp. This megabase scaffold number is the same as the number of chromosomes observed in karyotyping of this insect. The TPB genome is known to have high repetitive DNA content, and the reduced assembled genome size compared to flowcytometric estimates of approximately 860 Mbp may be due to the collapsed assembly of highly similar regions.

A major 6 Mb superlocus is involved in pyrethroid resistance in the common bed bug Cimex lectularius

In the last few years, the bed bug Cimex lectularius has been an increasing problem worldwide, mainly due to the development of insecticide resistance to pyrethroids. The characterization of resistance alleles is a prerequisite to improve surveillance and resistance management. To identify genomic variants associated with pyrethroid resistance in Cimex lectularius, we compared the genetic composition of two recent and resistant populations with that of two ancient-susceptible strains using a genome-wide pool-seq design. We identified a large 6 Mb "superlocus" showing particularly high genetic differentiation and association with the resistance phenotype. This superlocus contained several clustered resistance genes and was also characterized by a high density of structural variants (inversions, duplications). The possibility that this superlocus constitutes a resistance "supergene" that evolved after the clustering of alleles adapted to insecticide and after reduction in recombination is discussed.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

Despite the huge human and economic costs of invasive insects, which are the main group of invasive species, their environmental impacts through various mechanisms remain inadequately explained in databases and much of the invasion biology literature. High-throughput sequencing technology, especially whole-genome sequencing, has been used as a powerful method to study the mechanisms through which insects achieve invasion. In this study, we reviewed whole-genome sequencing-based advances in revealing several important invasion mechanisms of invasive insects, including (1) the rapid genetic variation and evolution of invasive populations, (2) invasion history and dispersal paths, (3) rapid adaptation to different host plant ranges, (4) strong environmental adaptation, (5) the development of insecticide resistance, and (6) the synergistic damage caused by invasive insects and endosymbiotic bacteria. We also discussed prevention and control technologies based on whole-genome sequencing and their prospects.

The role of structural variants in pest adaptation and genome evolution of the Colorado potato beetle, Leptinotarsa decemlineata (Say)

Structural variation has been associated with genetic diversity and adaptation. Despite these observations, it is not clear what their relative importance is for evolution, especially in rapidly adapting species. Here, we examine the significance of structural polymorphisms in pesticide resistance evolution of the agricultural super-pest, the Colorado potato beetle, Leptinotarsa decemlineata. By employing a parent offspring trio sequencing procedure, we develop highly contiguous reference genomes to characterize structural variation. These updated assemblies represent >100-fold improvement of contiguity and include derived pest and ancestral nonpest individuals. We identify >200,000 structural variations, which appear to be nonrandomly distributed across the genome as they co-occur with transposable elements and genes. Structural variations intersect with exons in a large proportion of gene annotations (~20%) that are associated with insecticide resistance (including cytochrome P450s), development, and transcription. To understand the role structural variations play in adaptation, we measure their allele frequencies among an additional 57 individuals using whole genome resequencing data, which represents pest and nonpest populations of North America. Incorporating multiple independent tests to detect the signature of natural selection using SNP data, we identify 14 genes that are probably under positive selection, include structural variations, and SNPs of elevated frequency within the pest lineages. Among these, three are associated with insecticide resistance based on previous research. One of these genes, CYP4g15, is coinduced during insecticide exposure with glycosyltransferase-13, which is a duplicated gene enclosed within a structural variant adjacent to the CYP4g15 genic region. These results demonstrate the significance of structural variations as a genomic feature to describe species history, genetic diversity, and adaptation.

Pest risk assessment of Amyelois transitella for the European Union

Following a request from the European Commission, the EFSA Panel on Plant Health performed a pest risk assessment of Amyelois transitella (Lepidoptera: Pyralidae), the navel orangeworm, for the EU. The quantitative assessment considered two scenarios: (i) current practices and (ii) a requirement for chilled transport. The assessment focused on pathways of introduction, climatic conditions and cultivation of hosts allowing establishment, spread and impact. A. transitella is a common pest of almonds, pistachios and walnuts in California, which is the main source for these nuts imported into the EU. Based on size of the trade and infestation at origin, importation of walnuts and almonds from the USA was identified as the most important pathways for entry of A. transitella. Using expert knowledge elicitation (EKE) and pathway modelling, a median estimate of 2,630 infested nuts is expected to enter the EU each year over the next 5 years (90% certainty range (CR) from 338 to 26,000 infested nuts per year). However, due to estimated small likelihoods of transfer to a host, mating upon transfer and survival of founder populations, the number of populations that establish was estimated to be 0.000698 year-1 (median, 90% CR: 0.0000126-0.0364 year-1). Accordingly, the expected period between founding events is 1,430 years (median, 90% CR: 27.5-79,400 year). The likelihood of entry resulting in establishment is therefore considered very small. However, this estimate has high uncertainty, mainly concerning the processes of transfer of the insect to hosts and the establishment of founder populations by those that successfully transfer. Climate matching and CLIMEX modelling indicate that conditions are most suitable for establishment in the southern EU, especially around the Mediterranean basin. The median rate of natural spread was estimated to be 5.6 km/year (median, 90% CR 0.8-19.3 km/year), after an initial lag period of 3.1 year (mean, 90% CR 1.7-6.2 year) following the establishment of a founder population. If A. transitella did establish, estimated median yield losses in nuts were estimated to be in the order of 1-2% depending on the nut species and production system. A scenario requiring imports of nuts to be transported under chilled conditions was shown to provide potential to further reduce the likelihood of entry.

Utility and challenges of using whole-genome resequencing to detect emerging insect and mite resistance in agroecosystems

Arthropods that invade agricultural ecosystems systematically evolve resistance to the control measures used against them, and this remains a significant and ongoing challenge for sustainable food production systems. Early detection of resistance evolution could prompt remedial action to slow the spread of resistance alleles in the landscape. Historical approaches used to detect emerging resistance included phenotypic monitoring of agricultural pest populations, as well as monitoring of allele frequency changes at one or a few candidate pesticide resistance genes. In this article, I discuss the successes and limitations of these traditional monitoring approaches and then consider whether whole-genome scanning could be applied to samples collected from agroecosystems over time for resistance monitoring. I examine the qualities of agroecosystems that could impact application of this approach to pesticide resistance monitoring and describe a recent retrospective analysis where genome scanning successfully detected an oligogenic response to selection by pesticides years prior to pest management failure. I conclude by considering areas of further study that will shed light on the feasibility of applying whole-genome scanning for resistance risk monitoring in agricultural pest species.

Evidence of hard-selective sweeps suggests independent adaptation to insecticides in Colorado potato beetle (Coleoptera: Chrysomelidae) populations

Pesticide resistance provides one of the best examples of rapid evolution to environmental change. The Colorado potato beetle (CPB) has a long and noteworthy history as a super-pest due to its ability to repeatedly develop resistance to novel insecticides and rapidly expand its geographic and host plant range. Here, we investigate regional differences in demography, recombination, and selection using whole-genome resequencing data from two highly resistant CPB populations in the United States (Hancock, Wisconsin and Long Island, New York). Demographic reconstruction corroborates historical records for a single pest origin during the colonization of the Midwestern and Eastern United States in the mid- to late-19th century and suggests that the effective population size might be higher in Long Island, NY than Hancock, WI despite contemporary potato acreage of Wisconsin being far greater. Population-based recombination maps show similar background recombination rates between these populations, as well as overlapping regions of low recombination that intersect with important metabolic detoxification genes. In both populations, we find compelling evidence for hard selective sweeps linked to insecticide resistance with multiple sweeps involving genes associated with xenobiotic metabolism, stress response, and defensive chemistry. Notably, only two candidate insecticide resistance genes are shared among both populations, but both appear to be independent hard selective sweep events. This suggests that repeated, rapid, and independent evolution of genes may underlie CPB's pest status among geographically distinct populations.

Loss and fixation of strongly favoured new variants: Understanding and extending Haldane's result via the Wright-Fisher model

The implications of Haldane's analysis for the fixation of beneficial alleles lies at the heart of much of 'population genetic thinking' and underlies many approaches that have been tailored to the detection of positive selection. Within the framework of a branching process, Haldane gave an approximation for the probability that fixation ultimately occurs when the selective advantage of a beneficial allele is small (≪1). Here, we make no use of branching processes. Rather, we work solely within a finite-population Wright-Fisher framework. We use this framework to analyse where Haldane's result applies, and extend Haldane's analysis. In particular, we present results for: (i) the domain of applicability of Haldane's analysis; (ii) the probability that loss occurs up to a given time; (iii) the probabilities that loss and fixation ultimately occur; (iv) an analytic approximation associated with the probability of loss and fixation ultimately occurring; (v) quantification of the crossover from weak to strong selection; (vi) determination of the number of invasive alleles that have a significant probability (>0.95) of invading a novel population. We note that the results obtained for (ii), (iii) and (iv) hold for an arbitrary initial number of mutations, and for selection that can be arbitrarily strong. Our results have fundamental implications for population and conservation genetics, and open up new avenues to identify traces of historically beneficial alleles through comparative genomics.

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.

RNAseq-based gene expression profiling of the Anopheles funestus pyrethroid-resistant strain FUMOZ highlights the predominant role of the duplicated CYP6P9a/b cytochrome P450s

Insecticide-based interventions, notably long-lasting insecticidal nets, against mosquito vectors of malaria are currently threatened by pyrethroid resistance. Here, we contrasted RNAseq-based gene expression profiling of laboratory-resistant (FUMOZ) and susceptible (FANG) strains of the major malaria vector Anopheles funestus. Cytochrome P450 genes were the predominant over-expressed detoxification genes in FUMOZ, with high expression of the duplicated CYP6P9a (fold-change of 82.23 vs FANG) and CYP6P9b (FC 11.15). Other over-expressed P450s belonged to the same cluster of P450s corresponding to the resistance to pyrethroid 1 (rp1) quantitative trait loci (QTL) on chromosome 2R. Several Epsilon class glutathione S-transferases were also over-expressed in FUMOZ, as was the ATP-binding cassette transporter AFUN019220 (ABCA) which also exhibited between-strain alternative splicing events at exon 7. Significant differences in single-nucleotide polymorphism frequencies between strains occurred in resistance QTLs rp1 (CYP6P9a/b and CYP6AA1), rp2 on chromosome 2L (CYP6Z1, CYP6M7, and CYP6Z3), and rp3 on chromosome 3R (CYP9J5, CYP9J4, and CYP9J3). Differences were also detected in CYP4G17 and CYP4G16 genes on the X chromosome, both of which are associated with cuticular resistance in Anopheles gambiae. A close analysis of nonsynonymous diversity at the CYP6P9a/b loci revealed a drastic loss of diversity in FUMOZ with only a single polymorphism and 2 haplotypes vs 18 substitutions and 8 haplotypes in FANG. By contrast, a lowly expressed cytochrome P450 (CYP4C36) did not exhibit diversity differences between strains. We also detected the known pyrethroid resistance conferring amino acid change N384S in CYP6P9b. This study further elucidates the molecular bases of resistance in An. funestus, informing strategies to better manage widespread resistance across Africa.

Genome-Wide Selective Signature Analysis Revealed Insecticide Resistance Mechanisms in Cydia pomonella

The codling moth, Cydia pomonella L. (Lepidoptera, Tortricidae), is a serious invasive pest of pome fruits. Currently, C. pomonella management mainly relies on the application of insecticides, which have driven the development of resistance in the insect. Understanding the genetic mechanisms of insecticide resistance is of great significance for developing new pest resistance management techniques and formulating effective resistance management strategies. Using existing genome resequencing data, we performed selective sweep analysis by comparing two resistant strains and one susceptible strain of the insect pest and identified seven genes, among which, two (glycine receptor and glutamate receptor) were under strong insecticide selection, suggesting their functional importance in insecticide resistance. We also found that eight genes including CYP6B2, CYP307a1, 5-hydroxytryptamine receptor, cuticle protein, and acetylcholinesterase, are potentially involved in cross-resistance to azinphos-methyl and deltamethrin. Moreover, among several P450s identified as positively selected genes, CYP6B2, CYP4C1, and CYP4d2 showed the highest expression level in larva compared to other stages tested, and CYP6B2 also showed the highest expression level in midgut, supporting the roles they may play in insecticide metabolism. Our results provide several potential genes that can be studied further to advance understanding of complexity of insecticide resistance mechanisms in C. pomonella.

Differential regulation of cytochrome P450 genes associated with biosynthesis and detoxification in bifenthrin-resistant populations of navel orangewom (Amyelois transitella)

Pyrethroid resistance was first reported in 2013 for the navel orangeworm, Amyelois transitella, but the genetic underpinnings of pyrethroid resistance are unknown. We investigated the role of cytochrome P450 monooxygenases (P450s) belonging to the CYP3 and CYP4 clans using colonies derived from individuals collected in 2016 from almond orchards in two counties. One colony (ALM) originated from an almond orchard in Madera County with no reported pyrethroid resistance and the second colony (R347) originated from the same Kern County orchard where pyrethroid resistance was first reported. We used high-throughput quantitative real-time PCR (qRT-PCR) analyses of 65 P450s in the CYP3 and CYP4 clans of A. transitella to identify P450s induced by bifenthrin and associated with pyrethroid resistance. Nine P450s were constitutively overexpressed in R347 compared to ALM, including CYP6AE54 (11.7-fold), belonging to a subfamily associated with metabolic pesticide detoxification in Lepidoptera and CYP4G89 (33-fold) belonging to a subfamily associated with cuticular hydrocarbon (CHC) synthesis and resistance via reduced pesticide penetrance. Cuticular hydrocarbons analysis revealed that R347 produced twice as many total CHCs in the egg and adult stages as ALM. Topical toxicity bioassays for R347 determined that egg mortality was reduced at low bifenthrin concentrations and larval mortality was reduced at high concentrations of bifenthrin compared to ALM. Our discovery of both changes in metabolism and production of CHCs for R347 have implications for the possible decreased efficacy of other classes of insecticide used to control this insect. The threat of widespread pyrethroid resistance combined with the potential for cross-resistance to develop through the mechanism of reduced penetrance warrants developing management strategies that facilitate insecticide passage across the cuticle.