Rapid evolution in insect pests: the importance of space and time in population genomics studies

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

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.

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.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Population genomic insights into invasion success in a polyphagous agricultural pest, Halyomorpha halys

Invasive species are increasingly threatening ecosystems and agriculture by rapidly expanding their range and adapting to environmental and human-imposed selective pressures. The genomic mechanisms that underlie such rapid changes remain unclear, especially for agriculturally important pests. Here, we used genome-wide polymorphisms derived from native, invasive, and intercepted samples and populations of the brown marmorated stink bug (BMSB), Halyomorpha halys, to gain insights into population genomics processes that have promoted the successful global invasion of this polyphagous pest. Our analysis demonstrated that BMSB exhibits spatial structure but admixture rates are high among introduced populations, resulting in similar levels of genomic diversity across native and introduced populations. These spatial genomic patterns suggest a complex invasion scenario, potentially with multiple bridgehead events, posing a challenge for accurately assigning BMSB incursions to their source using reduced-representation genomic data. By associating allele frequencies with the invasion status of BMSB populations, we found significantly differentiated single nucleotide polymorphisms (SNPs) located in close proximity to genes for insecticide resistance and olfaction. Comparing variations in allele frequencies among populations for outlier SNPs suggests that BMSB invasion success has probably evolved from standing genetic variation. In addition to being a major nuisance of households, BMSB has caused significant economic losses to agriculture in recent years and continues to expand its range. Despite no record of BMSB insecticide resistance to date, our results show high capacity for potential evolution of such traits, highlighting the need for future sustainable and targeted management strategies.

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.

Advances in data-collection tools and analytics for crop pest and disease management

Innovative methods in data collection and analytics for pest and disease management are advancing together with computational efficiency. Tools, such as the open-data kit, research electronic data capture, fall armyworm monitoring, and early warning- system application and remote sensing have aided the efficiency of all types of data collection, including text, location, images, audio, video, and others. Concurrently, data analytics have also evolved with the application of artificial intelligence and machine learning (ML) for early warning and decision-support systems. ML has repeatedly been used for the detection, diagnosis, modeling, and prediction of crop pests and diseases. This paper thus highlights the innovations, implications, and future progression of these technologies for sustainability.

Polymorphisms in the hypervariable control region of the mitochondrial DNA differentiate BPH populations

The brown planthopper (BPH; Nilaparvata lugens) is one of India's most destructive pests of rice. BPH, a monophagous migratory insect, reported from all major rice-growing ecosystems of the country, is capable of traversing large distances and causing massive crop loss. A crucial step for developing viable management strategies is understanding its population dynamics. Very few reliable markers are currently available to screen BPH populations for their diversity. In the current investigation, we developed a combinatorial approach using the polymorphism present within the mitochondrial Control Region of BPH and in the nuclear genome (genomic simple sequence repeats; gSSRs) to unravel the diversity present in BPH populations collected from various rice-growing regions of India. Using two specific primer pairs, the complete Control Region (1112 to 2612 bp) was PCR amplified as two overlapping fragments, cloned and sequenced from BPH individuals representing nine different populations. Results revealed extensive polymorphism within this region due to a variable number of tandem repeats. The three selected gSSR markers also exhibited population-specific amplification patterns. Overall genetic diversity between the nine populations was high (>5%). Further, in silico double-digestion of the consensus sequences of the Control Region, with HpyCH4IV and Tsp45I restriction enzymes, revealed unique restriction fragment length polymorphisms (digital-RFLPs; dRFLPs) that differentiated all the nine BPH populations. To the best of our knowledge, this is the first report of markers developed from the Control Region of the BPH mitogenome that can differentiate populations. Eventually, such reliable and rapid marker-based identification of BPH populations will pave the way for an efficient pest management strategy.

Unfolding the mitochondrial genome structure of green semilooper (Chrysodeixis acuta Walker): An emerging pest of onion (Allium cepa L.)

Onion is the most important crop challenged by a diverse group of insect pests in the agricultural ecosystem. The green semilooper (Chrysodeixis acuta Walker), a widespread tomato and soybean pest, has lately been described as an emergent onion crop pest in India. C. acuta whole mitochondrial genome was sequenced in this work. The circular genome of C. acuta measured 15,743 base pairs (bp) in length. Thirteen protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and one control region were found in the 37 sequence elements. With an average 395 bp gene length, the maximum and minimum gene length observed was 1749 bp and 63 bp of nad5 and trnR, respectively. Nine of the thirteen PCGs have (ATN) as a stop codon, while the other four have a single (T) as a stop codon. Except for trnS1, all of the tRNAs were capable of producing a conventional clover leaf structure. Conserved ATAGA motif sequences and poly-T stretch were identified at the start of the control region. Six overlapping areas and 18 intergenic spacer regions were found, with sizes ranged from 1 to 20 bp and 1 to 111 bp correspondingly. Phylogenetically, C. acuta belongs to the Plusiinae subfamily of the Noctuidae superfamily, and is closely linked to Trichoplusia ni species from the same subfamily. In the present study, the emerging onion pest C. acuta has its complete mitochondrial genome sequenced for the first time.

Geographic dispersion of invasive crop pests: the role of basal, plastic climate stress tolerance and other complementary traits in the tropics

Global pest invasions have significantly increased in recent years. These invasions together with climate warming directly impact agriculture. Tropical climates feature extreme weather events, including high temperatures and seasonal droughts. Thus, successful invasive pests in tropics have to adapt to these extreme climate features. The intrinsic factors relevant to tropical invasion of insects have been explored in many studies, but the knowledge is rather dispersed in contemporary literature. Here, we reviewed the potential biophysical characters of successful invasive pests' adaption to tropical environments including [1] inherent high basal stress tolerance and advanced life-history performances [2], phenotypic plasticity [3], rapid evolution to environmental stress, polyphagy, diverse reproductive strategies and high fecundity. We summarised how these traits and their interactive effects enhance pest invasions in the tropics. Comprehensive understanding of how these characters facilitate invasion improves models for predicting ecological consequences of climate change on invasive pest species for improved pest management.

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.

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

Genome evolution in an agricultural pest following adoption of transgenic crops

Evolution of resistance to management approaches in agricultural landscapes is common and results in economic losses. Early detection of pest resistance prior to significant crop damage would benefit the agricultural community. It has been hypothesized that new genomic approaches could track molecular signals of emerging resistance and trigger efforts to preempt widespread damage. We tested this hypothesis by quantifying genomic changes in the pest Helicoverpa zea over a 15-y period concurrent with commercialization of transgenic Bacillus thuringiensis–expressing crops and their subsequent loss of efficacy. Our results demonstrate the complex nature of evolution in agricultural ecosystems and provide insight into the potential and pitfalls of using genomic approaches for resistance monitoring.

Replacing synthetic insecticides with transgenic crops for pest management has been economically and environmentally beneficial, but these benefits erode as pests evolve resistance. It has been proposed that novel genomic approaches could track molecular signals of emerging resistance to aid in resistance management. To test this, we quantified patterns of genomic change in Helicoverpa zea, a major lepidopteran pest and target of transgenic Bacillus thuringiensis (Bt) crops, between 2002 and 2017 as both Bt crop adoption and resistance increased in North America. Genomic scans of wild H. zea were paired with quantitative trait locus (QTL) analyses and showed the genomic architecture of field-evolved Cry1Ab resistance was polygenic, likely arising from standing genetic variation. Resistance to pyramided Cry1A.105 and Cry2Ab2 toxins was controlled by fewer loci. Of the 11 previously described Bt resistance genes, 9 showed no significant change over time or major effects on resistance. We were unable to rule out a contribution of aminopeptidases (apns), as a cluster of apn genes were found within a Cry-associated QTL. Molecular signals of emerging Bt resistance were detectable as early as 2012 in our samples, and we discuss the potential and pitfalls of whole-genome analysis for resistance monitoring based on our findings. This first study of Bt resistance evolution using whole-genome analysis of field-collected specimens demonstrates the need for a more holistic approach to examining rapid adaptation to novel selection pressures in agricultural ecosystems.

Understanding Organismal Capacity to Respond to Anthropogenic Change: Barriers and Solutions

Global environmental changes induced by human activities are forcing organisms to respond at an unprecedented pace. At present we have only a limited understanding of why some species possess the capacity to respond to these changes while others do not. We introduce the concept of multidimensional phenospace as an organizing construct to understanding organismal evolutionary responses to environmental change. We then describe five barriers that currently challenge our ability to understand these responses: 1) Understanding the parameters of environmental change and their fitness effects, 2) Mapping and integrating phenotypic and genotypic variation, 3) Understanding whether changes in phenospace are heritable, 4) Predicting consistency of genotype to phenotype patterns across space and time, and 5) Determining which traits should be prioritized to understand organismal response to environmental change. For each we suggest one or more solutions that would help us surmount the barrier and improve our ability to predict, and eventually manipulate, organismal capacity to respond to anthropogenic change. Additionally, we provide examples of target species that could be useful to examine interactions between phenotypic plasticity and adaptive evolution in changing phenospace.

Population genomics and phylogeography of the boll weevil, Anthonomus grandis Boheman (Coleoptera: Curculionidae), in the United States, northern Mexico, and Argentina

The boll weevil, Anthonomus grandis Boheman (Coleoptera: Curculionidae), is an important pest of commercial cotton across the Americas. In the United States, eradication of this species is complicated by re-infestations of areas where eradication has been previously successful and by the existence of morphologically similar variants that can confound identification efforts. To date, no study has applied a high-throughput sequencing approach to better understand the population genetic structure of the boll weevil. Furthermore, only a single study has investigated genetic relationships between populations in North and South America. We used double digest restriction site-associated DNA sequencing (ddRADseq) to resolve the population genomic structure of the boll weevil in the southern United States, northern Mexico, and Argentina. Additionally, we assembled the first complete mitochondrial genome for this species and generated a preliminary whole genome assembly, both of which were used to improve the identification of informative loci. Downstream analyses revealed two main lineages-one consisting of populations found geographically west of the Sierra Madre Occidental mountain range and the second consisting of populations found to the east-were revealed, and both were sub-structured. Population geographic structure was consistent with the isolation by distance model, indicating that geogrpahic distance is likely a primary mechanism driving divergence in this species. Boll weevil populations from Argentina were found to be more closely related to the eastern lineage, suggesting a recent colonization of South America by the eastern lineage, but additional sampling across Mexico, Central America and South America is needed to further clarify their origin. Finally, we uncovered an instance of population turnover or replacement, highlighting the temporal instability of population structure.

Population structure, adaptation and divergence of the meadow spittlebug, Philaenus spumarius (Hemiptera, Aphrophoridae), revealed by genomic and morphological data

Understanding patterns of population differentiation and gene flow in insect vectors of plant diseases is crucial for the implementation of management programs of disease. We investigated morphological and genome-wide variation across the distribution range of the spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera, Auchenorrhyncha, Aphrophoridae), presently the most important vector of the plant pathogenic bacterium Xylella fastidiosa Wells et al., 1987 in Europe. We found genome-wide divergence between P. spumarius and a very closely related species, P. tesselatus Melichar, 1899, at RAD sequencing markers. The two species may be identified by the morphology of male genitalia but are not differentiated at mitochondrial COI, making DNA barcoding with this gene ineffective. This highlights the importance of using integrative approaches in taxonomy. We detected admixture between P. tesselatus from Morocco and P. spumarius from the Iberian Peninsula, suggesting gene-flow between them. Within P. spumarius, we found a pattern of isolation-by-distance in European populations, likely acting alongside other factors restricting gene flow. Varying levels of co-occurrence of different lineages, showing heterogeneous levels of admixture, suggest other isolation mechanisms. The transatlantic populations of North America and Azores were genetically closer to the British population analyzed here, suggesting an origin from North-Western Europe, as already detected with mitochondrial DNA. Nevertheless, these may have been produced through different colonization events. We detected SNPs with signatures of positive selection associated with environmental variables, especially related to extremes and range variation in temperature and precipitation. The population genomics approach provided new insights into the patterns of divergence, gene flow and adaptation in these spittlebugs and led to several hypotheses that require further local investigation.

Spatial population genomics of a recent mosquito invasion

Population genomic approaches can characterize dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG) and four incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232-577. Close kin dyads revealed recent movement between islands 31-203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and nonadjacent islands. Private alleles and a co-ancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system.

Disjunction between canola distribution and the genetic structure of its recently described pest, the canola flower midge (Contarinia brassicola)

Population genomics is a useful tool to support integrated pest management as it can elucidate population dynamics, demography, and histories of invasion. Here, we use a restriction site-associated DNA sequencing approach combined with whole-genome amplification (WGA) to assess genomic population structure of a newly described pest of canola, the diminutive canola flower midge, Contarinia brassicola. Clustering analyses recovered little geographic structure across the main canola production region but differentiated several geographically disparate populations at edges of the agricultural zone. Given a lack of alternative hypotheses for this pattern, we suggest these data support alternative hosts for this species and thus our canola-centric view of this midge as a pest has limited our understanding of its biology. These results speak to the need for increased surveying efforts across multiple habitats and other potential hosts within Brassicaceae to improve both our ecological and evolutionary knowledge of this species and contribute to effective management strategies. We additionally found that use of WGA prior to library preparation was an effective method for increasing DNA quantity of these small insects prior to restriction site-associated DNA sequencing and had no discernible impact on genotyping consistency for population genetic analysis; WGA is therefore likely to be tractable for other similar studies that seek to randomly sample markers across the genome in small organisms.

What Is the Spatial Extent of a Bemisia tabaci Population?

Effective pest management depends on basic knowledge about insect dispersal patterns and gene flow in agroecosystems. The globally invasive sweet potato whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is considered a weak flier whose life history nonetheless predisposes it to frequent dispersal, but the scale over which populations exchange migrants, and should therefore be managed, is uncertain. In this review, we synthesize the emergent literature on B. tabaci population genetics to address the question: What spatial scales define B. tabaci populations? We find that within-species genetic differentiation among sites is often low, and evidence of population structuring by host plant or geography is rare. Heterozygote deficits prevail among populations, indicating that migrants from divergent populations are frequently sampled together. Overall, these results suggest that there is high ongoing gene flow over large spatial extents. However, genetic homogeneity typical of recently invading populations could obscure power to detect real isolation among populations. Genome-wide data collected systematically across space and time could distinguish signatures of invasion history from those of ongoing gene flow. Characterizing the spatial extent of B. tabaci populations could reveal whether insecticide rotations can be tailored to specific commodities or if coordination across linked commodities and regions is justified.

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.

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

Among the most familiar forms of human-driven evolution on ecological time scales is the rapid acquisition of resistance to pesticides by insects. Since the widespread adoption of synthetic organic insecticides in the mid-twentieth century, over 500 arthropod species have evolved resistance to at least one insecticide. Efforts to determine the genetic bases of insecticide resistance have historically focused on individual loci, but the availability of genomic tools has facilitated the screening of genome-wide characteristics. We resequenced three contemporary populations of the navel orangeworm (Amyelois transitella), the principal pest of almond orchards in California, differing in bifenthrin resistance status to examine insecticide-induced changes in the population genomic landscape of this species. We detected an exceptionally large region with virtually no polymorphisms, extending to up to 1.3 Mb in the resistant population. This selective sweep includes genes associated with pyrethroid and DDT resistance, including a cytochrome P450 gene cluster and the gene encoding the voltage-gated sodium channel para. Moreover, the sequence along the sweep is nearly identical in the genome assembled from a population founded in 1966, suggesting that the foundation for insecticide resistance may date back a half-century, when California’s Central Valley experienced massive area-wide applications of DDT for pest control.

Variable resistance to spinetoram in populations of Thrips palmi across a small area unconnected to genetic similarity

The melon thrips, Thrips palmi, is an increasingly important pest of vegetables in northern China. Some populations have developed resistance in the field to the insecticide spinetoram. Understanding the origin and dispersal of insecticide-resistant populations can shed light on resistance management strategies. In this study, we tested susceptibility of seven greenhouse populations of T. palmi to spinetoram collected from a small area of about 300 km2 in Shandong Province and examined population genetic structure across the area based on a segment of mitochondrial cox1 gene and 22 microsatellite loci to infer the possible origin and dispersal of insecticide resistance. Levels of resistance to spinetoram differed among seven populations, which included one population with high resistance (LC50 = 759.34 mg/L), three populations with medium resistance (LC50 ranged from 28.69 to 34.79 mg/L), and three populations with low resistance (LC50 ranged from 7.61 to 8.97 mg/L). The populations were genetically differentiated into two groups unrelated to both levels of resistance and geographic distance. The molecular data indicated high levels of gene flow between populations with different levels of resistance to spinetoram and low gene flow among populations with the same level of resistance, pointing to a likely separate history of resistance evolution. Resistance levels of two tested populations to spinetoram decreased 23 and 4.6 times after five generations without any exposure to the pesticide. We therefore suspect that resistance of T. palmi most likely evolved in response to local applications of the insecticide. Our study suggests that the development of resistance could be avoided or resistance even reversed by reducing usage of spinetoram.

Contrasting population structure and demographic history of cereal aphids in different environmental and agricultural landscapes

Genetic diversity of populations has important ecological and evolutionary consequences, whose understanding is fundamental to improve the sustainability of agricultural production. Studies of how differences in agricultural management and environment influence the population structure of insect pests are central to predict outbreaks and optimize control programs. Here, we have studied the population genetic diversity and evolution of Sitobion avenae and Sitobion miscanthi, which are among the most relevant aphid pests of cereals across Europe and Asia, respectively. We have used genotyping by sequencing (GBS) to identify genome-wide single nucleotide polymorphisms (SNPs) to infer the geographic structure and migration patterns. In the present study, we show that the population structure in present-day populations is different from that described in previous studies, which suggest that they have evolved recently possibly as a response to human-induced changes in agriculture. This study shows that S. avenae in England is predominantly parthenogenetic and there has been a demographic and spatial expansion of a single genetic cluster, which could correspond with the insecticide resistance superclone identified in previous studies. Conversely, in China, S. miscanthi populations are mostly cyclical parthenogenetic, with one sexual stage in autumn to produce overwintering eggs, and there are six genetically differentiated subpopulations and high genetic differentiation between geographic locations, which suggests that further taxonomical research is needed. Unlike S. avenae in England, there is no evidence for insecticide resistance and there is no predominance of a single lineage in S. miscanthi in China.

Genome-wide analysis of diamondback moth, Plutella xylostella L., from Brassica crops and wild host plants reveals no genetic structure in Australia

Molecular studies of population structure can reveal insight into the movement patterns of mobile insect pests in agricultural landscapes. The diamondback moth, Plutella xylostella L., a destructive pest of Brassica vegetable and oilseed crops worldwide, seasonally colonizes winter canola crops in southern Australia from alternative host plant sources. To investigate movement, we collected 59 P. xylostella populations from canola crops, Brassica vegetable and forage crops and brassicaceous wild host plants throughout southern Australia in 2014 and 2015 and genotyped 833 individuals using RAD-seq for genome-wide analysis. Despite a geographic sampling scale > 3,000 km and a statistically powerful set of 1,032 SNP markers, there was no genetic differentiation among P. xylostella populations irrespective of geographic location, host plant or sampling year, and no evidence for isolation-by-distance. Hierarchical STRUCTURE analysis at K = 2–5 showed nearly uniform ancestry in both years. Cluster analysis showed divergence of a small number of individuals at several locations, possibly reflecting an artefact of sampling related individuals. It is likely that genetic homogeneity within Australian P. xylostella largely reflects the recent colonization history of this species but is maintained through some level of present gene flow. Use of genome-wide neutral markers was uninformative for revealing the seasonal movements of P. xylostella within Australia, but may provide more insight in other global regions where the species has higher genetic diversity.

Genetic and phenotypic consequences of early domestication in black soldier flies (Hermetia illucens)

The black soldier fly, Hermetia illucens, is an emerging biotechnological agent with its larvae being effective converters of organic waste into usable bio-products including protein and lipids. To date, most operations use unimproved commercial populations produced by mass rearing, without cognisance of specific breeding strategies. The genetic and phenotypic consequences of these commercial practices remain unknown and could have a significant impact on long-term population viability and productivity. The aim of this study was thus to assess the genetic and phenotypic changes during the early phases of colony establishment and domestication in the black soldier fly. An experimental colony was established from wild founder flies and a new microsatellite marker panel was developed to assess population genetic parameters along with the phenotypic characteristics of each generational cohort under captive breeding. The experimental colony was characterised by a small effective population size, subsequent loss of genetic diversity and rapid genetic and phenotypic differentiation between the generational cohorts. Ultimately, the population collapsed by the fifth generation, most likely owing to the adverse effect of inbreeding depression following the fixation of deleterious alleles. Species with r-selected life history characteristics (e.g. short life-span, high fecundity and low larval survival) are known to pose particular challenges for genetic management. The current study suggests that sufficient genetic and phenotypic variations exist in the wild population and that domestication and strain development could be achieved with careful population augmentation and selection during the early stages of colony establishment.

A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues

High-throughput sequencing technologies are a proposed solution for accessing the molecular data in historical specimens. However, degraded DNA combined with the computational demands of short-read assemblies has posed significant laboratory and bioinformatics challenges for de novo genome assembly. Linked-read or "synthetic long-read" sequencing technologies, such as 10× Genomics, may provide a cost-effective alternative solution to assemble higher quality de novo genomes from degraded tissue samples. Here, we compare assembly quality (e.g., genome contiguity and completeness, presence of orthogroups) between four new deer mouse (Peromyscus spp.) genomes assembled using linked-read technology and four published genomes assembled from a single shotgun library. At a similar price-point, these approaches produce vastly different assemblies, with linked-read assemblies having overall higher contiguity and completeness, measured by larger N50 values and greater number of genes assembled, respectively. As a proof-of-concept, we used annotated genes from the four Peromyscus linked-read assemblies and eight additional rodent taxa to generate a phylogeny, which reconstructed the expected relationships among species with 100% support. Although not without caveats, our results suggest that linked-read sequencing approaches are a viable option to build de novo genomes from degraded tissues, which may prove particularly valuable for taxa that are extinct, rare or difficult to collect.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

Evolutionary epidemiology in the field: a proactive approach for identifying herbicide resistance in problematic crop weeds

This article is a Commentary on Comont et al., 223: 1584–1594.

Comparative Proteomics Analysis Between the Short-Term Stress and Long-Term Adaptation of the Blattella germanica (Blattodea: Blattellidae) in Response to Beta-Cypermethrin

A proteomic method combining two-dimensional polyacrylamide gel electrophoresis and tandem mass spectrometry was used to compare the hemolymph expression profiles of a beta-cypermethrin-resistant Blattella germanica L. strain (R) and a susceptible strain (S) after 24 h of beta-cypermethrin induction. The results showed that there were 42 differentially expressed proteins after induction of the R strain: 4 proteins were upregulated and 38 proteins were downregulated. One hundred one hemolymph proteins were differentially expressed after induction of the S strain: 53 proteins were upregulated and 48 proteins were downregulated. The identified proteins were mainly classified into the following categories: energy metabolism proteins such as arginine kinase and triose phosphate isomerase, detoxification-related proteins such as glutathione S-transferases (GSTs), signal molecule-regulated proteins such as nitric oxide synthase (NOS), and other proteins such as kinetic-related proteins and gene expression-related proteins. Several proteins show significant differences in response to short-term stress and long-term adaptation, and differential expression of these proteins reflects an overall change in cellular structure and metabolism associated with resistance to pyrethroid insecticides. In summary, our research has improved the understanding of the molecular mechanisms of beta-cypermethrin resistance in German cockroaches, which will facilitate the development of rational methods to improve the management of this pest.

Patterns of genetic differentiation in Colorado potato beetle correlate with contemporary, not historic, potato land cover

Changing landscape heterogeneity can influence connectivity and alter genetic variation in local populations, but there can be a lag between ecological change and evolutionary responses. Temporal lag effects might be acute in agroecosystems, where land cover has changed substantially in the last two centuries. Here, we evaluate how patterns of an insect pest's genetic differentiation are related to past and present agricultural land cover change over a 150-year period. We quantified change in the amount of potato, Solanum tuberosum L., land cover since 1850 using county-level agricultural census reports, obtained allele frequency data from 7,408 single-nucleotide polymorphism loci, and compared effects of historic and contemporary landscape connectivity on genetic differentiation of Colorado potato beetle, Leptinotarsa decemlineata Say, in two agricultural landscapes in the United States. We found that potato land cover peaked in Wisconsin in the early 1900s, followed by rapid decline and spatial concentration, whereas it increased in amount and extent in the Columbia Basin of Oregon and Washington beginning in the 1960s. In both landscapes, we found small effect sizes of landscape resistance on genetic differentiation, but a 20× to 1,000× larger effect of contemporary relative to historic landscape resistances. Demographic analyses suggest population size trajectories were largely consistent among regions and therefore are not likely to have differentially impacted the observed patterns of population structure in each region. Weak landscape genetic associations might instead be related to the coarse resolution of our historical land cover data. Despite rapid changes in agricultural landscapes over the last two centuries, genetic differentiation among L. decemlineata populations appears to reflect ongoing landscape change. The historical landscape genetic framework employed in this study is broadly applicable to other agricultural pests and might reveal general responses of pests to agricultural land-use change.

Beetles as Model Organisms in Physiological, Biomedical and Environmental Studies - A Review

Model organisms are often used in biological, medical and environmental research. Among insects, Drosophila melanogaster, Galleria mellonella, Apis mellifera, Bombyx mori, Periplaneta americana, and Locusta migratoria are often used. However, new model organisms still appear. In recent years, an increasing number of insect species has been suggested as model organisms in life sciences research due to their worldwide distribution and environmental significance, the possibility of extrapolating research studies to vertebrates and the relatively low cost of rearing. Beetles are the largest insect order, with their representative - Tribolium castaneum - being the first species with a completely sequenced genome, and seem to be emerging as new potential candidates for model organisms in various studies. Apart from T. castaneum, additional species representing various Coleoptera families, such as Nicrophorus vespilloides, Leptinotarsa decemlineata, Coccinella septempunctata, Poecilus cupreus, Tenebrio molitor and many others, have been used. They are increasingly often included in two major research aspects: biomedical and environmental studies. Biomedical studies focus mainly on unraveling mechanisms of basic life processes, such as feeding, neurotransmission or activity of the immune system, as well as on elucidating the mechanism of different diseases (neurodegenerative, cardiovascular, metabolic, or immunological) using beetles as models. Furthermore, pharmacological bioassays for testing novel biologically active substances in beetles have also been developed. It should be emphasized that beetles are a source of compounds with potential antimicrobial and anticancer activity. Environmental-based studies focus mainly on the development and testing of new potential pesticides of both chemical and natural origin. Additionally, beetles are used as food or for their valuable supplements. Different beetle families are also used as bioindicators. Another important research area using beetles as models is behavioral ecology studies, for instance, parental care. In this paper, we review the current knowledge regarding beetles as model organisms and their practical application in various fields of life science.

Evolutionary insights from large scale resequencing datasets in Drosophila melanogaster

Drosophila melanogaster has long been used as an evolutionary model system. Its small genome size, well-annotated genome, and ease of sampling, also makes it a choice species for genome resequencing studies. Hundreds of genomic samples from populations worldwide are available and are currently being used to tackle a wide range of evolutionary questions. In this review, we focused on three insights that have increased our understanding of the evolutionary history of this species, and that have implications for the study of evolutionary processes in other species as well. Because of technical limitations, most of the studies so far have focused on SNP variants. However, long-read sequencing techniques should allow us in the near future to include other type of genomic variants that also influence genome evolution.

Geographic Variation in Larval Metabolic Rate Between Northern and Southern Populations of the Invasive Gypsy Moth

Thermal regimes can diverge considerably across the geographic range of a species, and accordingly, populations can vary in their response to changing environmental conditions. Both local adaptation and acclimatization are important mechanisms for ectotherms to maintain homeostasis as environments become thermally stressful, which organisms often experience at their geographic range limits. The spatial spread of the gypsy moth (Lymantria dispar L.) (Lepidoptera: Erebidae) after introduction to North America provides an exemplary system for studying population variation in physiological traits given the gradient of climates encompassed by its current invasive range. This study quantifies differences in resting metabolic rate (RMR) across temperature for four populations of gypsy moth, two from the northern and two from southern regions of their introduced range in North America. Gypsy moth larvae were reared at high and low thermal regimes, and then metabolic activity was monitored at four temperatures using stop-flow respirometry to test for an acclimation response. For all populations, there was a significant increase in RMR as respirometry test temperature increased. Contrary to our expectations, we did not find evidence for metabolic adaptation to colder environments based on our comparisons between northern and southern populations. We also found no evidence for an acclimation response of RMR to rearing temperature for three of the four pairwise comparisons examined. Understanding the thermal sensitivity of metabolic rate in gypsy moth, and understanding the potential for changes in physiology at range extremes, is critical for estimating continued spatial spread of this invasive species both under current and potential future climatic constraints.