Agricultural insect pests as models for studying stress-induced evolutionary processes

Agricultural insect pests (AIPs) are widely successful in adapting to natural and anthropogenic stressors, repeatedly overcoming population bottlenecks and acquiring resistance to intensive management practices. Although they have been largely overlooked in evolutionary studies, AIPs are ideal systems for understanding rapid adaptation under novel environmental conditions. Researchers have identified several genomic mechanisms that likely contribute to adaptive stress responses, including positive selection on de novo mutations, polygenic selection on standing allelic variation and phenotypic plasticity (e.g., hormesis). However, new theory suggests that stress itself may induce epigenetic modifications, which may confer heritable physiological changes (i.e., stress-resistant phenotypes). In this perspective, we discuss how environmental stress from agricultural management generates the epigenetic and genetic modifications that are associated with rapid adaptation in AIPs. We summarise existing evidence for stress-induced evolutionary processes in the context of insecticide resistance. Ultimately, we propose that studying AIPs offers new opportunities and resources for advancing our knowledge of stress-induced evolution.

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.

The Transcriptomic Response of the Boll Weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), following Exposure to the Organophosphate Insecticide Malathion

Insecticide tolerance and resistance have evolved countless times in insect systems. Molecular drivers of resistance include mutations in the insecticide target site and/or gene duplication, and increased gene expression of detoxification enzymes. The boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), is a pest of commercial cotton and has developed resistance in the field to several insecticides; however, the current organophosphate insecticide, malathion, used by USA eradication programs remains effective despite its long-term use. Here, we present findings from an RNA-seq experiment documenting gene expression post-exposure to field-relevant concentrations of malathion, which was used to provide insight on the boll weevil's continued susceptibility to this insecticide. Additionally, we incorporated a large collection of boll weevil whole-genome resequencing data from nearly 200 individuals collected from three geographically distinct areas to determine SNP allele frequency of the malathion target site, as a proxy for directional selection in response to malathion exposure. No evidence was found in the gene expression data or SNP data consistent with a mechanism of enhanced tolerance or resistance adaptation to malathion in the boll weevil. Although this suggests continued effectiveness of malathion in the field, we identified important temporal and qualitative differences in gene expression between weevils exposed to two different concentrations of malathion. We also identified several tandem isoforms of the detoxifying esterase B1 and glutathione S-transferases, which are putatively associated with organophosphate resistance.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae)

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

Rapid evolution of insecticide resistance in the Colorado potato beetle, Leptinotarsa decemlineata

Despite considerable research, efforts to manage insecticide resistance continue to fail. The Colorado potato beetle (CPB), Leptinotarsa decemlineata Say (Coleoptera: Chrysomelidae), epitomizes this problem, as it has repeatedly and rapidly evolved resistance to>50 insecticides. The patterns of resistance evolution are intriguing, as they defy models where resistance evolves from rare mutations. Here, we synthesize recent research on insecticide resistance in CPB showing that polygenic resistance drawn from standing genetic diversity explains genomic patterns of insecticide resistance evolution. However, rapid gene regulatory evolution suggests that other mechanisms might also facilitate adaptive change. We explore the hypothesis that sublethal stress from insecticide exposure could alter heritable epigenetic modifications, and discuss the range of experimental approaches needed to fully understand insecticide resistance evolution in this super pest.

Museum Genomics of an Agricultural Super-Pest, the Colorado Potato Beetle, Leptinotarsa decemlineata (Chrysomelidae), Provides Evidence of Adaptation from Standing Variation

Despite extensive research on agricultural pests, our knowledge about their evolutionary history is often limited. A mechanistic understanding of the demographic changes and modes of adaptation remains an important goal, as it improves our understanding of organismal responses to environmental change and our ability to sustainably manage pest populations. Emerging genomic datasets now allow for characterization of demographic and adaptive processes, but face limits when they are drawn from contemporary samples, especially in the context of strong demographic change, repeated selection, or adaptation involving modest shifts in allele frequency at many loci. Temporal sampling, however, can improve our ability to reconstruct evolutionary events. Here, we leverage museum samples to examine whether population genomic diversity and structure has changed over time, and to identify genomic regions that appear to be under selection. We focus on the Colorado potato beetle (CPB), Leptinotarsa decemlineata (Say 1824; Coleoptera: Chrysomelidae), which is widely regarded as a super-pest due to its rapid, and repeated, evolution to insecticides. By combining whole genome resequencing data from 78 museum samples with modern sampling, we demonstrate that CPB expanded rapidly in the 19th century, leading to a reduction in diversity and limited genetic structure from the Midwest to Northeast United States. Temporal genome scans provide extensive evidence for selection acting in resistant field populations in Wisconsin and New York, including numerous known insecticide resistance genes. We also validate these results by showing that known selective sweeps in modern populations are identified by our genome scan. Perhaps most importantly, temporal analysis indicates selection on standing genetic variation, as we find evidence for parallel evolution in the two geographical regions. Parallel evolution involves a range of phenotypic traits not previously identified as under selection in CPB, such as reproductive and morphological functional pathways that might be important for adaptation to agricultural habitats.

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.

Elevated rates of positive selection drive the evolution of pestiferousness in the Colorado potato beetle (Leptinotarsa decemlineata, Say)

Contextualizing evolutionary history and identifying genomic features of an insect that might contribute to its pest status is important in developing early detection and control tactics. In order to understand the evolution of pestiferousness, which we define as the accumulation of traits that contribute to an insect population's success in an agroecosystem, we tested the importance of known genomic properties associated with rapid adaptation in the Colorado potato beetle (CPB), Leptinotarsa decemlineata Say. Within the leaf beetle genus Leptinotarsa, only CPB, and a few populations therein, has risen to pest status on cultivated nightshades, Solanum. Using whole genomes from ten closely related Leptinotarsa species native to the United States, we reconstructed a high-quality species tree and used this phylogenetic framework to assess evolutionary patterns in four genomic features of rapid adaptation: standing genetic variation, gene family expansion and contraction, transposable element abundance and location, and positive selection at protein-coding genes. Throughout approximately 20 million years of history, Leptinotarsa species show little evidence of gene family turnover and transposable element variation. However, there is a clear pattern of CPB experiencing higher rates of positive selection on protein-coding genes. We determine that these rates are associated with greater standing genetic variation due to larger effective population size, which supports the theory that the demographic history contributes to rates of protein evolution. Furthermore, we identify a suite of coding genes under positive selection that are putatively associated with pestiferousness in the Colorado potato beetle lineage. They are involved in the biological processes of xenobiotic detoxification, chemosensation and hormone function.