Population structure and gene flow in the global pest, Helicoverpa armigera

Historical museum samples reveal signals of selection and drift in response to changing insecticide use in an agricultural pest moth

In response to environmental and human-imposed selective pressures, agroecosystem pests frequently undergo rapid evolution, with some species having a remarkable capacity to rapidly develop pesticide resistance. Temporal sampling of genomic data can comprehensively capture such adaptive changes over time, for example, by elucidating allele frequency shifts in pesticide resistance loci in response to different pesticides. Here, we leveraged museum specimens spanning over a century of collections to generate temporal contrasts between pre- and post-insecticide populations of an agricultural pest moth, Helicoverpa armigera. We used targeted exon sequencing of 254 samples collected across Australia from the pre-1950s (prior to insecticide introduction) to the 1990s, encompassing decades of changing insecticide use. Our sequencing approach focused on genes that are known to be involved in insecticide resistance, environmental sensation, and stress tolerance. We found an overall lack of spatial and temporal population structure change across Australia. In some decades (e.g., 1960s and 1970s), we found a moderate reduction of genetic diversity, implying stochasticity in evolutionary trajectories due to genetic drift. Temporal genome scans showed extensive evidence of selection following insecticide use, although the majority of selected variants were low impact. Finally, alternating trajectories of allele frequency change were suggestive of potential antagonistic pleiotropy. Our results provide new insights into recent evolutionary responses in an agricultural pest and show how temporal contrasts using museum specimens can improve mechanistic understanding of rapid evolution.

Lure specificity, phenology, and damage caused by Epiphyas moths (Lepidoptera: Tortricidae) in Western Australian apple orchards

Multiple Epiphyas species inhabit southwestern Western Australia, including Light Brown Apple Moth (LBAM) Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae), a globally significant, polyphagous pest. This study evaluated the efficacy and specificity of lures designed for 3 Epiphyas species: E. postvittana, Epiphyas pulla (Turner), and the undescribed Epiphyas sp. (1) (Common). Additionally, the study sought to determine the presence and distribution of Epiphyas species in 3 significant apple-growing localities. Trapping, together with partial sequencing of the mitochondrial COI gene, found LBAM to be restricted to the Perth Hills and E. pulla, to apple orchards near Manjimup and Pemberton. This geographic disjunction remains unexplained. Epiphyas sp. (1) was not recorded despite using a specifically designed lure. The E. pulla and LBAM traps demonstrated superior efficacy in capturing their target species, while the catch in Epiphyas sp. (1) traps did not significantly differ between the 2. Both E. pulla and LBAM exhibited peak abundance from late spring to the end of summer (October-February), with variations in timing and peak catch of male moths across species, locations, and years. Surveys conducted in April during the harvest period (February-May), when moth traps caught an average of 1-1.8 moths/trap/week, found no Epiphyas larvae or damage on 140,400 mature apples or on 26,000 leaves. While E. pulla and LBAM traps effectively monitor their target moths, genetic identification of trap catch would be necessary if they co-occurred. Encouragingly, the results indicate that both species become relatively rare as harvest season approaches, and neither inflicts significant damage to mature apples under existing management.

Chromosome-level genome assembly and comparative genomics shed light on Helicoverpa assulta ecology and pest management

BACKGROUND

The Oriental tobacco budworm, Helicoverpa assulta, a specialist herbivorous insect that exclusively feeds on plants of the Solanaceae family, causes considerable damage to crops, such as tobacco and hot pepper. The absence of a genome sequence for this species hinders further research on its pest management and ecological adaptation.

RESULTS

Here, we present a high-quality chromosome-level genome of a Korean strain of H. assulta (Pyeongchang strain, K18). The total assembly spans 424.4 Mb with an N50 length of 14.54 Mb and 37% GC content. The assembled genome (ASM2961881v1) comprises 31 chromosomes, similar to other congeneric generalist species including H. armigera and H. zea. In terms of genomic assembly quality, the complete BUSCOs and repeat content accounted for 98.3% and 33.01% of the genome, respectively. Based on this assembly, 19 485 protein-coding genes were predicted in the genome annotation. A comparative analysis was conducted using the identified number of protein-coding genes in H. armigera (24154) and H. zea (23696). Out of the 19 485 predicted genes, 137 genes in 15 orthogroups were found to have expanded significantly in H. assulta, while 149 genes in 95 orthogroups contracted rapidly.

CONCLUSION

This study revealed specific gene expansions and contractions in H. assulta compared to those in its close relatives, indicating potential adaptations related to its specialized feeding habits. Also, the comparative genome analysis provides valuable insights for the integrated pest management of H. assulta and other globally significant pests in the Heliothinae subfamily. © 2024 Society of Chemical Industry.

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Virulence of Beauveria sp. and Metarhizium sp. fungi towards fall armyworm (Spodoptera frugiperda)

The development of effective pest management strategies for Spodoptera frugiperda is a high priority for crop protection across its invasive ranges. Here, we examined six Beauveria and five Metarhizium fungal isolates against this pest. Two Beauveria isolates (B-0571, B-1311) induced high mortality toward 3rd and 6th instar caterpillars and adults. For B-0571 mortality was 82.81 ± 5.75%, 61.46 ± 6.83%, and 93.75 ± 3.61%, and 73.72 ± 2.51%, 71.88 ± 5.41%, and 97.92 ± 2.08% for B-1311, with deaths in caterpillars largely occurring under 24 h (3rd instar control 0.74 ± 0.33%, B-0571 73.96 ± 7.85% and B-1311 62.08 ± 3.67%; 6th instar control 0%, B-0571 66.67% ± 11.02% and B-1311 62.5% ± 9.55%). Infection from both Beauveria isolates fully prevented reproduction in surviving S. frugiperda females. In contrast, all five Metarhizium isolates tested and the remaining four Beauveria isolates exhibited lower virulence. The discovery of two highly virulent Beauveria fungal isolates to S. frugiperda opens avenues to develop novel biological control tools against this highly invasive pest.

Adaptive evolution to the natural and anthropogenic environment in a global invasive crop pest, the cotton bollworm

The cotton bollworm, Helicoverpa armigera, is set to become the most economically devastating crop pest in the world, threatening food security and biosafety as its range expands across the globe. Key to understanding the eco-evolutionary dynamics of H. armigera, and thus its management, is an understanding of population connectivity and the adaptations that allow the pest to establish in unique environments. We assembled a chromosome-scale reference genome and re-sequenced 503 individuals spanning the species range to delineate global patterns of connectivity, uncovering a previously cryptic population structure. Using a genome-wide association study (GWAS) and cell line expression of major effect loci, we show that adaptive changes in a temperature- and light-sensitive developmental pathway enable facultative diapause and that adaptation of trehalose synthesis and transport underlies cold tolerance in extreme environments. Incorporating extensive pesticide resistance monitoring, we also characterize a suite of novel pesticide and Bt resistance alleles under selection in East China. These findings offer avenues for more effective management strategies and provide insight into how insects adapt to variable climatic conditions and newly colonized environments.

The Biological Hierarchy, Time, and Temporal 'Omics in Evolutionary Biology: A Perspective

Sequencing data-genomics, transcriptomics, epigenomics, proteomics, and metabolomics-have revolutionized biological research, enabling a more detailed study of processes, ranging from subcellular to evolutionary, that drive biological organization. These processes, collectively, are responsible for generating patterns of phenotypic variation and can operate over dramatically different timescales (milliseconds to billions of years). While researchers often study phenotypic variation at specific levels of biological organization to isolate processes operating at that particular scale, the varying types of sequence data, or 'omics, can also provide complementary inferences to link molecular and phenotypic variation to produce an integrated view of evolutionary biology, ranging from molecular pathways to speciation. We briefly describe how 'omics has been used across biological levels and then demonstrate the utility of integrating different types of sequencing data across multiple biological levels within the same study to better understand biological phenomena. However, single-time-point studies cannot evaluate the temporal dynamics of these biological processes. Therefore, we put forward temporal 'omics as a framework that can better enable researchers to study the temporal dynamics of target processes. Temporal 'omics is not infallible, as the temporal sampling regime directly impacts inferential ability. Thus, we also discuss the role the temporal sampling regime plays in deriving inferences about the environmental conditions driving biological processes and provide examples that demonstrate the impact of the sampling regime on biological inference. Finally, we forecast the future of temporal 'omics by highlighting current methodological advancements that will enable temporal 'omics to be extended across species and timescales. We extend this discussion to using temporal multi-omics to integrate across the biological hierarchy to evaluate and link the temporal dynamics of processes that generate phenotypic variation.

Resistance Bioassays and Allele Characterization Inform Analysis of Spodoptera frugiperda (Lepidoptera: Noctuidae) Introduction Pathways in Asia and Australia

The fall armyworm (FAW) Spodoptera frugiperda (Smith; Lepidoptera: Noctuidae) is present in over 70 countries in Africa, Asia, and Oceania. Its rapid dispersal since 2016 when it was first reported in western Africa, and associated devastation to agricultural productivity, highlight the challenges posed by this pest. Currently, its management largely relies on insecticide sprays and transgenic Bacillus thuringiensis toxins, therefore understanding their responses to these agents and characteristics of any resistance genes enables adaptive strategies. In Australia, S. frugiperda was reported at the end of January 2020 in northern Queensland and by March 2020, also in northern Western Australia. As an urgent first response we undertook bioassays on two Australian populations, one each from these initial points of establishment. To assist with preliminary sensitivity assessment, two endemic noctuid pest species, Helicoverpa armigera (Hübner; Lepidoptera, Noctuidae) and Spodoptera litura (Fabricius; Lepidoptera, Noctuidae), were concurrently screened to obtain larval LC50 estimates against various insecticides. We characterized known resistance alleles from the VGSC, ACE-1, RyR, and ABCC2 genes to compare with published allele frequencies and bioassay responses from native and invasive S. frugiperda populations. An approximately 10× LC50 difference for indoxacarb was detected between Australian populations, which was approximately 28× higher than that reported from an Indian population. Characterization of ACE-1 and VGSC alleles provided further evidence of multiple introductions in Asia, and multiple pathways involving genetically distinct individuals in Australia. The preliminary bioassay results and resistance allele patterns from invasive S. frugiperda populations suggest multiple introductions have contributed to the pest's spread and challenge the axiom of its rapid 'west-to-east' spread.

Subspecific Synonym of Monochamus alternatus (Coleoptera: Cerambycidae): Population Genetics and Morphological Reassessment

Monochamus alternatus Hope, 1842, is a major forest pest that hosts the pathogenic pinewood nematode (PWN), Bursaphelenchus xylophilus (Steiner and Buhrer, 1934) Nickle 1970. Taxonomically, M. alternatus is currently divided into two subspecies, based on morphology and geography: Monochamus alternatus alternatus Hope, 1842 in China, Taiwan, Tibet, Vietnam, and Laos and Monochamus alternatus endai Makihara, 2004 in South Korea and Japan. Despite their economic importance, the subspecies taxonomy of M. alternatus has never been tested after the first description. In this study, we aimed to reassess the subspecies taxonomy of M. alternatus using molecular and morphological data. For morphological analysis, we examined three major morphological characters (pronotal longitudinal band, granulation on humeri, and elytral proximomedial spine) from 191 individuals from China, Korea, and Taiwan. Population genetic structures were examined using 85 de novo sequences and 82 public COI sequences from China, Korea, Japan, Malaysia, Taiwan, and a few intercepted specimens from the United States. All the genetic data were aligned as three different multiple sequence alignments. Individuals from each subspecies were morphologically and genetically scattered, not clustered according to subspecies in any of the analyses. Therefore, a new synonymy is proposed: Monochamus alternatus Hope, 1842 = Monochamus alternatus endai, syn. n. This study suggests a more robust classification of M. alternatus for the first time and ultimately will pose a substantial impact on implementing quarantine or forestry policies.

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.

Genetic Variation and Population Structure of the Old World Bollworm Helicoverpa armigera (Hübner, 1808) (Lepidoptera: Noctuidae) in Ethiopia

Helicoverpa armigera is one of the most destructive insect pests of economically valuable crops in the world. Despite its economic importance, the population genetic structure of this insect remains unexplored in Ethiopia. To investigate the genetic diversity and population structure of H. armigera, we sampled 170 individuals from 15 populations throughout Ethiopia. We sequenced a fragment of the mitochondrial cytochrome b (cyt b) gene and five exon-primed intron-crossing (EPIC) markers. Twenty cyt b haplotypes with low-to-moderate haplotype diversity (mean Hd = 0.537) and high nucleotide diversity (mean Pi = 0.00339) were identified. The most frequently observed and widely distributed cyt b haplotype was designated as Hap_1 (67.058%), which is identical to sequences found across the globe. Tajima's D and Fu's F for the cyt b data were negative, supporting a model of population expansion. Within populations, a mean of 2.493 alleles/locus was recorded across the five EPIC loci, ranging from 1.200 to 3.600 alleles/locus. The highest mean effective number of alleles/population was 2.369 and the lowest was 1.178. The mean observed heterozygosity (HO) of the five loci (0-0.289; mean 0.104 ± 0.020) was lower than the expected heterozygosity (HE) (0.095-0.523; mean 0.258 ± 0.028). AMOVA detected significant genetic structure with 61% of the total molecular genetic variation of EPIC genotypes occurring between populations, suggesting a considerable degree of differentiation among populations. STRUCTURE analyses clustered the H. armigera populations into three distinct population groups but very low isolation by distance (R2 = 0.0132, P < 0.05).

Slow and steady wins the race: spatial and stochastic processes and the failure of suppression gene drives

Gene drives that skew sex ratios offer a new management tool to suppress or eradicate pest populations. Early models and empirical work suggest that these suppression drives can completely eradicate well‐mixed populations, but models that incorporate stochasticity and space (i.e. drift and recolonization events) often result in loss or failure of the drive. We developed a stochastic model to examine these processes in a simple one‐dimensional space. This simple space allows us to map the events and outcomes that emerged and examine how properties of the drive's wave of invasion affect outcomes. Our simulations, across a biologically realistic section of parameter space, suggest that drive failure might be a common outcome in spatially explicit, stochastic systems, and that properties of the drive wave appear to mediate outcomes. Surprisingly, the drives that would be considered fittest in an aspatial model were strongly associated with failure in the spatial setting. The fittest drives cause relatively fast moving, and narrow waves that have a high chance of being penetrated by wild‐types (WTs) leading to WT recolonization, leading to failure. Our results also show that high rates of dispersal reduce the chance of failure because drive waves get disproportionately wider than WT waves as dispersal rates increase. Overall, wide, slow‐moving drive waves were much less prone to failure. Our results point to the complexity inherent in using a genetic system to effect demographic outcomes and speak to a clear need for ecological and evolutionary modelling to inform the drive design process.

Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda

Bean leaf beetle (BLB) (Ootheca mutabilis) has emerged as an important bean pest in Uganda, leading to devastating crop losses. There is limited information on the population genetic structure of BLB despite its importance. In this study, novel microsatellite DNA markers and the partial mitochondrial cytochrome oxidase subunit I (mtCOI) gene sequences were used to analyze the spatial population genetic structure, genetic differentiation and haplotype diversity of 86 O. mutabilis samples from 16 (districts) populations. We identified 19,356 simple sequence repeats (SSRs) (mono, di-, tri-, tetra-, penta-, and hexa-nucleotides) of which 81 di, tri and tetra-nucleotides were selected for primer synthesis. Five highly polymorphic SSR markers (4-21 alleles, heterozygosity 0.59-0.84, polymorphic information content (PIC) 50.13-83.14%) were used for this study. Analyses of the 16 O. mutabilis populations with these five novel SSRs found nearly all the genetic variation occurring within populations and there was no evidence of genetic differentiation detected for both types of markers. Also, there was no evidence of isolation by distance between geographical and genetic distances for SSR data and mtCOI data except in one agro-ecological zone for mtCOI data. Bayesian clustering identified a signature of admixture that suggests genetic contributions from two hypothetical ancestral genetic lineages for both types of markers, and the minimum-spanning haplotype network showed low differentiation in minor haplotypes from the most common haplotype with the most common haplotype occurring in all the 16 districts. A lack of genetic differentiation indicates unrestricted migrations between populations. This information will contribute to the design of BLB control strategies.

Genetic diversity and demographic history of the Old World Bollworm, Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae), in Ethiopia inferred from mitochondrial gene sequences

The Old World bollworm, Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae), is a globally distributed agricultural and horticultural insect pest. Despite the economic importance of this insect in Ethiopia, its genetic diversity and demographic history are poorly understood. We examined the nucleotide variation of the mitochondrial cytochrome c oxidase subunit I (COI) gene fragment of 74 H. armigera individuals from six collection sites in Ethiopia. We recorded 15 COI haplotypes in H. armigera, ten globally shared and five exclusive to Ethiopia (HaET15, HaET14, HaET10, HaET7, and HaET4). Haplotype HaET1 was the most widely geographically distributed and frequent (71.62%). Analysis of molecular variance (AMOVA) revealed a high and significant level of variation within H. armigera populations (θ_ST = −0.0135). Negative values of the neutrality test and nonsignificant index of mismatch distribution supported the demographic expansion of H. armigera populations in Ethiopia; furthermore, this was also supported by the nonsignificant values of the sum of squared deviations (SSD) and raggedness index (r). The high genetic variation and population expansion of H. armigera have immense implications for devising locally adapted management strategies in area‐wide integrated pest management IPM programs. However, a comprehensive study of H. armigera genetic diversity and population structure using various molecular markers is needed for future confirmation.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest’s global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

This population genomics study identifies the complex multiple introduction history of Spodoptera frugiperda (fall armyworm) from the Americas, into Africa and Asia. This provides new insight into the ‘east-to-west’ directionality of gene flow, and suggests ample genomic exchange at the nuclear level.

Helicoverpa armigera and Helicoverpa zea hybridization: constraints, heterosis, and implications for pest management

BACKGROUND

The invasion of Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) into the New World has made it possible for this pest to hybridize with a native American species, H. zea (Boddie), under natural conditions. We investigated the viability and development of hybrids of these two Helicoverpa species. We reared the parental species and evaluated crosses between H. armigera males and H. zea females and vice versa, two intercrosses between hybrids, and eight backcrosses between hybrids and parental species. We estimated the length of immature stages, fecundity, survival, sex ratio, and heterosis.

RESULTS

Although hybridization occcurred, with heterosis during the development of immatures, reproductive incompatibilities also were observed between the parental species and between hybrids from subsequent crosses. The interspecific crosses between hybrids and backcrosses confirmed the possibility of introgression events and their perpetuation in field populations. The results indicate that hybridization events are favored at high population levels, while at low population levels the 'species identities' will be maintained.

CONCLUSIONS

The possibility of interspecific gene flow and its perpetuation through successive crosses and backcrosses suggests several recommenations for management. Populations of both species should be maintained at an equilibrium level to reduce the chance of interspecific crosses, which are presumably more likely to occur during pest outbreaks. The existence of hybridization and resistance to different active pesticide ingredients should be monitored. All practices related to managing the resistance of these pests to chemical and biological insecticides should be systematized to reduce the chance of selecting for resistant individuals.

Population genetic structure of the carrot weevil ( Listronotus oregonensis ) in North America

Population genetic studies of insect pests enhance our ability to anticipate problems in agroecosystems, such as pest outbreaks, insecticide resistance, or expansions of the host range. This study focuses on geographic distance and host plant selection as potential determinants of genetic differentiation of the carrot weevil Listronotus oregonensis, a major pest of several apiaceous crops in North America. To undertake genetic studies on this species, we assembled the first complete genome sequence for L. oregonensis. Then, we used both haplotype discrimination with mitochondrial DNA (mtDNA) and a genotyping‐by‐sequencing (GBS) approach to characterize the genetic population structure. A total of 220 individuals were sampled from 17 localities in the provinces of Québec, Ontario, Nova Scotia (Canada), and the state of Ohio (USA). Our results showed significant genetic differences between distant populations across North America, indicating that geographic distance represents an important factor of differentiation for the carrot weevil. Furthermore, the GBS analysis revealed more different clusters than COI analysis between Québec and Nova Scotia populations, suggesting a recent differentiation in the latter province. In contrast, we found no clear evidence of population structure associated with the four cultivated apiaceous plants tested (carrot, parsley, celery, and celeriac) using populations from Québec. This first characterization of the genetic structure of the carrot weevil contributes to a better understanding of the gene flow of the species and helps to adapt local pest management measures to better control this agricultural pest.

Climate Change Increases the Expansion Risk of Helicoverpa zea in China According to Potential Geographical Distribution Estimation

Simple Summary

Helicoverpa zea is one of the most destructive lepidopteran agricultural pests in the world and can disperse long distances both with and without human transportation. It is listed in the catalog of quarantine pests for plants imported to the People’s Republic of China but has not yet been reported in China. On the basis of 1781 global distribution records of H. zea and eight bioclimatic variables, we predicted the potential geographical distributions (PGDs) of H. zea by using a calibrated MaxEnt model. The results showed that the PGDs of H. zea under the current climate are large in China. Future climate changes under shared socioeconomic pathways (SSP) 1-2.6, SSP2-4.5, and SSP5-8.5 for both the 2030s and 2050s will facilitate the expansion of PGDs for H. zea. Helicoverpa zea has a high capacity for colonization by introduced individuals in China. Customs ports should pay attention to the host plants of H. zea and containers harboring this pest.

Abstract

Helicoverpa zea, a well-documented and endemic pest throughout most of the Americas, affecting more than 100 species of host plants. It is a quarantine pest according to the Asia and Pacific Plant Protection Commission (APPPC) and the catalog of quarantine pests for plants imported to the People’s Republic of China. Based on 1781 global distribution records of H. zea and eight bioclimatic variables, the potential geographical distributions (PGDs) of H. zea were predicted by using a calibrated MaxEnt model. The contribution rate of bioclimatic variables and the jackknife method were integrated to assess the significant variables governing the PGDs. The response curves of bioclimatic variables were quantitatively determined to predict the PGDs of H. zea under climate change. The results showed that: (1) four out of the eight variables contributed the most to the model performance, namely, mean diurnal range (bio2), precipitation seasonality (bio15), precipitation of the driest quarter (bio17) and precipitation of the warmest quarter (bio18); (2) PGDs of H. zea under the current climate covered 418.15 × 10⁴ km², and were large in China; and (3) future climate change will facilitate the expansion of PGDs for H. zea under shared socioeconomic pathways (SSP) 1-2.6, SSP2-4.5, and SSP5-8.5 in both the 2030s and 2050s. The conversion of unsuitable to low suitability habitat and moderately to high suitability habitat increased by 8.43% and 2.35%, respectively. From the present day to the 2030s, under SSP1-2.6, SSP2-4.5 and SSP5-8.5, the centroid of the suitable habitats of H. zea showed a general tendency to move eastward; from 2030s to the 2050s, under SSP1-2.6 and SSP5-8.5, it moved southward, and it moved slightly northward under SSP2-4.5. According to bioclimatic conditions, H. zea has a high capacity for colonization by introduced individuals in China. Customs ports should pay attention to host plants and containers of H. zea and should exchange information to strengthen plant quarantine and pest monitoring, thus enhancing target management.

Population genomics provides insights into lineage divergence and local adaptation within the cotton bollworm

The cotton bollworm Helicoverpa armigera is a cosmopolitan pest and its diverse habitats plausibly contribute to the formation of diverse lineages. Despite the significant threat it poses to economic crops worldwide, its evolutionary history and genetic basis of local adaptation are poorly understood. In this study, we de novo assembled a high-quality chromosome-level reference genome of H. a. armigera (contig N50 = 7.34 Mb), with 99.13% of the HaSCD2 assembly assigned into 31 chromosomes (Z-chromosome + 30 autosomes). We constructed an ultra-dense variation map across 14 cotton bollworm populations and identified a novel lineage in northwestern China. Historical inference showed that effective population size changes coincided with global temperature fluctuation. We identified nine differentiated genes in the three H. armigera lineages (H. a. armigera, H. a. conferta, and the new northwestern Chinese lineage), of which per and clk genes are involved in circadian rhythm. Selective sweep analyses identified a series of GO categories related to climate adaptation, feeding behavior and insecticide tolerance. Our findings reveal fundamental knowledge of the local adaptation of different cotton bollworm lineages and will guide the formulation of cotton bollworm management measures at different scales.

Determinants of Insecticide Resistance Evolution: Comparative Analysis Among Heliothines

It is increasingly clear that pest species vary widely in their propensities to develop insecticide resistance. This review uses a comparative approach to analyze the key pest management practices and ecological and biochemical or genetic characteristics of the target that contribute to this variation. We focus on six heliothine species, three of which, Helicoverpa armigera, Heliothis virescens, and Helicoverpa zea, have developed resistances to many pesticide classes. The three others, Helicoverpa punctigera, Helicoverpa assulta, and Helicoverpa gelotopoeon, also significant pests, have developed resistance to very few pesticide classes. We find that host range and movement between alternate hosts are key ecological traits that influence effective selection intensities for resistance. Operational issues are also critical; area-wide, cross-pesticide management practices that account for these ecological factors are key to reducing selection intensity. Without such management, treatment using broad-spectrum chemicals serves to multiply the effects of host plant preference, preadaptive detoxification ability, and high genetic diversity to create a pesticide treadmill for the three high-propensity species.Without rigorous ongoing management, such a treadmill could still develop for newer, more selective chemistries and insecticidal transgenic crops.

Anatomical Comparison of Antennal Lobes in Two Sibling Ectropis Moths: Emphasis on the Macroglomerular Complex

Ectropis obliqua and Ectropis grisescens are two sibling moth species of tea plantations in China. The male antennae of both species can detect shared and specific sex pheromone components. Thus, the primary olfactory center, i.e., the antennal lobe (AL), plays a vital role in distinguishing the sex pheromones. To provide evidence for the possible mechanism allowing this distinction, in this study, we compared the macroglomerular complex (MGC) of the AL between the males of the two species by immunostaining using presynaptic antibody and propidium iodide (PI) with antennal backfills, and confocal imaging and digital 3D-reconstruction. The results showed that MGC of both E. obliqua and E. grisescens contained five glomeruli at invariant positions between the species. However, the volumes of the anterior-lateral glomerulus (ALG) and posterior-ventral (PV) glomerulus differed between the species, possibly related to differences in sensing sex pheromone compounds and their ratios between E. obliqua and E. grisescens. Our results provide an important basis for the mechanism of mating isolation between these sibling moth species.

Whole-genome sequencing to detect mutations associated with resistance to insecticides and Bt proteins in Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a major pest native to the Americas that has recently invaded the Old World. Point mutations in the target-site proteins acetylcholinesterase-1 (ace-1), voltage-gated sodium channel (VGSC) and ryanodine receptor (RyR) have been identified in S. frugiperda as major resistance mechanisms to organophosphate, pyrethroid and diamide insecticides respectively. Mutations in the adenosine triphosphate-binding cassette transporter C2 gene (ABCC2) have also been identified to confer resistance to Cry1F protein. In this study, we applied a whole-genome sequencing (WGS) approach to identify point mutations in the target-site genes in 150 FAW individuals collected from China, Malawi, Uganda and Brazil. This approach revealed three amino acid substitutions (A201S, G227A and F290V) of S. frugiperda ace-1, which are known to be associated with organophosphate resistance. The Brazilian population had all three ace-1 point mutations and the 227A allele (mean frequency = 0.54) was the most common. Populations from China, Malawi and Uganda harbored two of the three ace-1 point mutations (A201S and F290V) with the 290V allele (0.47-0.58) as the dominant allele. Point mutations in VGSC (T929I, L932F and L1014F) and RyR (I4790M and G4946E) were not detected in any of the 150 individuals. A novel 12-bp insertion mutation in exon 15 of the ABCC2 gene was identified in some of the Brazilian individuals but absent in the invasive populations. Our results not only demonstrate robustness of the WGS-based genomic approach for detection of resistance mutations, but also provide insights for improvement of resistance management tactics in S. frugiperda.

On species delimitation, hybridization and population structure of cassava whitefly in Africa

The Bemisia cassava whitefly complex includes species that cause severe crop damage through vectoring cassava viruses in eastern Africa. Currently, this whitefly complex is divided into species and subgroups (SG) based on very limited molecular markers that do not allow clear definition of species and population structure. Based on 14,358 genome-wide SNPs from 62 Bemisia cassava whitefly individuals belonging to sub-Saharan African species (SSA1, SSA2 and SSA4), and using a well-curated mtCOI gene database, we show clear incongruities in previous taxonomic approaches underpinned by effects from pseudogenes. We show that the SSA4 species is nested within SSA2, and that populations of the SSA1 species comprise well-defined south-eastern (Madagascar, Tanzania) and north-western (Nigeria, Democratic Republic of Congo, Burundi) putative sub-species. Signatures of allopatric incipient speciation, and the presence of a 'hybrid zone' separating the two putative sub-species were also detected. These findings provide insights into the evolution and molecular ecology of a highly cryptic hemipteran insect complex in African, and allow the systematic use of genomic data to be incorporated in the development of management strategies for this cassava pest.

Migration trajectories of the diamondback moth Plutella xylostella in China inferred from population genomic variation

BACKGROUND

The diamondback moth (DBM), Plutella xylostella (Lepidoptera: Plutellidae), is a notorious pest of cruciferous plants. In temperate areas, annual populations of DBM originate from adult migrants. However, the source populations and migration trajectories of immigrants remain unclear. Here, we investigated migration trajectories of DBM in China using genome-wide single nucleotide polymorphisms (SNPs) genotyped using double-digest RAD (ddRAD) sequencing. We first analyzed patterns of spatial and temporal genetic structure among southern source and northern recipient populations, then inferred migration trajectories into northern regions using discriminant analysis of principal components (DAPC), assignment tests, and spatial kinship patterns.

RESULTS

Temporal genetic differentiation among populations was low, indicating that sources of recipient populations and migration trajectories are stable. Spatial genetic structure indicated three genetic clusters in the southern source populations. Assignment tests linked northern populations to the Sichuan cluster, and central-eastern populations to the southern and Yunnan clusters, indicating that Sichuan populations are sources of northern immigrants and southern and Yunnan populations are sources of central-eastern populations. First-order (full-sib) and second-order (half-sib) kin pairs were always found within populations, but ~ 35-40% of third-order (cousin) pairs were found in different populations. Closely related individuals in different populations were found at distances of 900-1500 km in ~ 35-40% of cases, while some were separated by > 2000 km.

CONCLUSION

This study unravels seasonal migration patterns in the DBM. We demonstrate how careful sampling and population genomic analyses can be combined to help understand cryptic migration patterns in insects. © 2020 Society of Chemical Industry.

Helicoverpa armigera Harm 1 Haplotype Predominates in the Heliothinae (Lepidoptera: Noctuidae) Complex Infesting Tomato Crops in Brazil

Tomato (Solanum lycopersicum L.) is a natural host for the Helicoverpa-Chloridea (Lepidoptera: Noctuidae: Heliothinae) pest complex. The species Helicoverpa armigera (Hübner) was responsible for significant yield losses in several crops after its detection in Brazil. The morphology of its larval stage resembles common Heliothinae species, making pest control decisions difficult. The overall lack of studies on the Heliothinae associated with tomatoes in Brazil and the establishment of H. armigera in the country plus their recent outbreaks supported our investigation about the relative importance of the insects from the Helicoverpa-Chloridea complex in this vegetable crop. A nationwide survey was carried out across fresh-market and processing tomato fields. Molecular analyses targeting a segment of the mitochondrial cytochrome oxidase subunit I (mtCOI) gene and their sequence analyses indicated the presence of a pest complex, comprising the introduced species H. armigera and the indigenous species, Helicoverpa zea (Boddie), and Chloridea virescens (Fabricius). The Harm 1 haplotype of H. armigera was identified as the predominant Heliothinae pest infesting fresh-market tomatoes. The New World species Chloridea subflexa (Guenée) as well as the exotic Solanaceae-specific species Helicoverpa assulta (Guenée) were not found in our survey. Additional larvae surveys in processing tomato fields during 2013/2014 in Central Brazil also indicated H. armigera as the most abundant Heliothinae species (95%) together with H. zea (4.75%) and C. virescens (0.25%). The occurrence of distinct Helicoverpa species (which are potentially capable of interbreeding) indicates that novel crop management strategies will be necessary in order to minimize damages caused by this pest complex in tomatoes.

Rapid parallel adaptation despite gene flow in silent crickets

Gene flow is predicted to impede parallel adaptation via de novo mutation, because it can introduce pre-existing adaptive alleles from population to population. We test this using Hawaiian crickets (Teleogryllus oceanicus) in which 'flatwing' males that lack sound-producing wing structures recently arose and spread under selection from an acoustically-orienting parasitoid. Morphometric and genetic comparisons identify distinct flatwing phenotypes in populations on three islands, localized to different loci. Nevertheless, we detect strong, recent and ongoing gene flow among the populations. Using genome scans and gene expression analysis we find that parallel evolution of flatwing on different islands is associated with shared genomic hotspots of adaptation that contain the gene doublesex, but the form of selection differs among islands and corresponds to known flatwing demographics in the wild. We thus show how parallel adaptation can occur on contemporary timescales despite gene flow, indicating that it could be less constrained than previously appreciated.

Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. Rhododendron is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of R. cyanocarpum, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in R. cyanocarpum, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, R. cyanocarpum had a higher genetic diversity (π = 0.0702 ± 0.0017, H_e = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (F_st from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (N_e) of R. cyanocarpum from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of N_e over a short period, and a subsequent gradual decline in N_e to date. Implications for conserving R. cyanocarpum based on these main results are then discussed.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

Multiple invasions of a generalist herbivore-Secondary contact between two divergent lineages of Nezara viridula Linnaeus in Australia

The presence of distinct evolutionary lineages within herbivorous pest insect taxa requires close attention. Scientific understanding, biosecurity planning and practice, and pest management decision-making each suffer when such situations remain poorly understood. The pest bug Nezara viridula Linnaeus has been recorded from numerous host plants and has two globally distributed mitochondrial (mtDNA) lineages. These mtDNA lineages co-occur in few locations globally, and the consequences of their divergence and recent secondary contact have not been assessed. We present evidence that both mtDNA lineages of N. viridula are present in Australia and their haplotype groups have a mostly separate distribution from one another. The north-western population has only Asian mtDNA haplotypes, and the population with an eastern distribution is characterized mostly by European mtDNA haplotypes. Haplotypes of both lineages were detected together at only one site in the north of eastern Australia, and microsatellite data indicate that this secondary contact has resulted in mating across the lineages. Admixture and the movement of mtDNA haplotypes outside of this limited area of overlap has not, however, been extensive. Some degree of mating incompatibility or differences in the climatic requirements and tolerances of the two lineages, and perhaps a combination of these influences, might limit introgression and the movement of individuals, but this needs to be tested. This work provides the foundation for further ecological investigation of the lineages of N. viridula, particularly the consequences of admixture on the ecology of this widespread pest. We propose that for now, the Asian and European lineages of N. viridula would best be investigated as subspecies, so that "pure" and admixed populations of this bug can each be considered directly with respect to management and research priorities.

Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control

Cryptic taxa have often been observed in the form of host‐associated species that diverged as the result of adaptation to alternate host plants. Untangling cryptic diversity in species complexes that encompass invasive species is a mandatory task for pest management. Moreover, investigating the evolutionary history of a species complex may help to understand the drivers of their diversification. The mealybug Hypogeococcus pungens was believed to be a polyphagous species from South America and has been reported as a pest devastating native cacti in Puerto Rico, also threatening cactus diversity in the Caribbean and North America. There is neither certainty about the identity of the pest nor the source population from South America. Recent studies pointed to substantial genetic differentiation among local populations, suggesting that H. pungens is a species complex. In this study, we used a combination of genome‐wide SNPs and mtDNA variation to investigate species diversity within H. pungens sensu lato to establish host plant ranges of each one of the putative members of the complex, to evaluate whether the pattern of host plant association drove diversification in the species complex, and to determine the source population of the Puerto Rican cactus pest. Our results suggested that H. pungens comprises at least five different species, each one strongly associated with specific host plants. We also established that the Puerto Rican cactus pest derives from southeastern Brazilian mealybugs. This is an important achievement because it will help to design reliable strategies for biological control using natural enemies of the pest from its native range.

Genetic Differentiation of a New World Screwworm Fly Population from Uruguay Detected by SNPs, Mitochondrial DNA and Microsatellites in Two Consecutive Years

The New World screwworm (NWS) fly, Cochliomyia hominivorax (Diptera: Calliphoridae), is an economically important ectoparasite currently distributed in South America and in the Caribbean basin. The successful eradication of this species in USA, Mexico and continental Central America was achieved by a control program based on the sterile insect technique (SIT). In order to implement a genetic control strategy over the NWS fly's current area of occurrence, first, it is necessary to understand the species dynamics and population structure. In order to address this objective, the spatial genetic structure of the NWS fly was previously reported in South America based on different genetic markers; however, to date, no study has investigated temporal changes in the genetic composition of its populations. In the current study, the temporal genetic structure of a NWS fly population from Uruguay was investigated through two consecutive samplings from the same locality over an interval of approximately 18 generations. The genetic structure was accessed with neutral and under selection SNPs obtained with genotyping-by-sequencing. The results gathered with these data were compared to estimates achieved with mitochondrial DNA sequences and eight microsatellite markers. Temporal changes in the genetic composition were revealed by all three molecular markers, which may be attributed to seasonal changes in the NWS fly's southern distribution. SNPs were employed for the first time for estimating the genetic structure in a NWS fly population; these results provide new clues and perspectives on its population genetic structure. This approach could have significant implications for the planning and implementation of management programs.

Arthropod Invasions Versus Soybean Production in Brazil: A Review

Soybean production in Brazil has been markedly affected by invasions of non-native arthropod species that feed on the crop, severely impacting biodiversity, food security, health, and economic development. Data on soybean production losses and increase in insecticide usage over the last two decades have not been explored in association with past invasion events, and the dynamics underlying the recent blitz of invasive species into Brazil remain largely unclear. We provide a review of arthropod invasions in the Brazilian soybean agroecosystem since 1990, indicating that the introductions of Bemisia tabaci (Gennadius) MEAM1 (Hemiptera: Aleyrodidae), Tetranychus urticae (Koch) (Acari: Tetranychidae), and Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) are likely correlated with periods of increase in insecticide usage for soybean production. Using these three cases as examples, we review factors that could lead to increased likelihood of future invasions by particular pests, outlining four possible criteria to evaluate potential invasiveness of non-native arthropods: likelihood of entry, likelihood of establishment, biological features of the species, and availability of control measures. Spodoptera litura (F.) (Lepidoptera: Noctuidae) and Aphis glycines (Matsumura) (Hemiptera: Sternorrhynca) are examples of highly damaging soybean pests, related to one or more of these factors, that could be introduced into Brazil over the next years and which could lead to problematic scenarios. Melanagromyza sojae (Zehnter) (Diptera: Agromyzidae) also meets these criteria and has successfully invaded and colonized Brazilian soybean fields in recent years. Our review identifies current issues within soybean pest management in Brazil and highlights the need to adopt management measures to offset future costs and minimize lost revenue.

Large-scale genetic admixture suggests high dispersal in an insect pest, the apple fruit moth

Knowledge about population genetic structure and dispersal capabilities is important for the development of targeted management strategies for agricultural pest species. The apple fruit moth, Argyresthia conjugella (Lepidoptera, Yponomeutidae), is a pre-dispersal seed predator. Larvae feed on rowanberries (Sorbus aucuparia), and when rowanberry seed production is low (i.e., inter-masting), the moth switches from laying eggs in rowanberries to apples (Malus domestica), resulting in devastating losses in apple crops. Using genetic methods, we investigated if this small moth expresses any local genetic structure, or alternatively if gene flow may be high within the Scandinavian Peninsula (~850.000 km², 55^o - 69^o N). Genetic diversity was found to be high (n = 669, mean He = 0.71). For three out of ten tetranucleotide STRs, we detected heterozygote deficiency caused by null alleles, but tests showed little impact on the overall results. Genetic differentiation between the 28 sampling locations was very low (average FST = 0.016, P < 0.000). Surprisingly, we found that all individuals could be assigned to one of two non-geographic genetic clusters, and that a third, geographic cluster was found to be associated with 30% of the sampling locations, with weak but significant signals of isolation-by-distance. Conclusively, our findings suggest wind-aided dispersal and spatial synchrony of both sexes of the apple fruit moth over large areas and across very different climatic zones. We speculate that the species may recently have had two separate genetic origins caused by a genetic bottleneck after inter-masting, followed by rapid dispersal and homogenization of the gene pool across the landscape. We suggest further investigations of spatial genetic similarities and differences of the apple fruit moth at larger geographical scales, through life-stages, across inter-masting, and during attacks by the parasitoid wasp (Microgaster politus).

Comparative Toxicity of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) to Selected Insecticides

Until recently, the Old World bollworm (OWB) Helicoverpa armigera (Hübner) and the corn earworm Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) were geographically isolated. Both species are major pests of agricultural commodities that are known to develop insecticide resistance, and they now coexist in areas where H. armigera invaded the Americas. This is the first study to compare the susceptibility of the two species to conventional insecticides. The susceptibility of third instar H. armigera and H. zea larvae to indoxacarb, methomyl, spinetoram, and spinosad was determined using a diet-overlay bioassay in a quarantine laboratory in Puerto Rico. Mortality was assessed at 48 h after exposure for up to eight concentrations per insecticide. Spinetoram exhibited the highest acute toxicity against H. armigera, with a median lethal concentration (LC50) of 0.11 µg a.i./cm2, followed by indoxacarb and spinosad (0.17 µg a.i./cm2 for both) and methomyl (0.32 µg a.i./cm2). Spinetoram was also the most toxic to H. zea (LC50 of 0.08 µg a.i./cm2), followed by spinosad (0.17 µg a.i./cm2) and methomyl (0.18 µg a.i./cm2). Indoxacarb was the least toxic to H. zea, with an LC50 of 0.21 µg a.i./cm2. These findings could serve as a comparative reference for monitoring the susceptibility of H. armigera and H. zea to indoxacarb, methomyl, spinetoram, and spinosad in Puerto Rico, and may facilitate the detection of field-selected resistance for these two species and their potential hybrids in areas recently invaded by H. armigera.

Oxidative stress delays development and alters gene expression in the agricultural pest moth, Helicoverpa armigera

Stress is a widespread phenomenon that all organisms must endure. Common in nature is oxidative stress, which can interrupt cell homeostasis to cause cell damage and may be derived from respiration or from environmental exposure through diet. As a result of the routine exposure from respiration, many organisms can mitigate the effects of oxidative stress, but less is known about responses to oxidative stress from other sources. Helicoverpa armigera is a major agricultural pest moth that causes significant damage to crops worldwide. Here, we examined the effects of oxidative stress on H. armigera by chronically exposing individuals to paraquat-a free radical producer-and measuring changes in development (weight, developmental rate, lifespan), and gene expression. We found that oxidative stress strongly affected development in H. armigera, with stressed samples spending more time as caterpillars than control samples (>24 vs. ~15 days, respectively) and therefore living longer overall. We found 1,618 up- and 761 down-regulated genes, respectively, in stressed versus control samples. In the up-regulated gene set, was an over-representation of biological processes related to cuticle and chitin development, glycine metabolism, and oxidation-reduction. Oxidative stress clearly impacts physiology and biochemistry in H. armigera and the interesting finding of an extended lifespan in stressed individuals could demonstrate hormesis, the phenomenon whereby toxic compounds can actually be beneficial at low doses. Collectively, our findings provide new insights into physiological and gene expression responses to oxidative stress in invertebrates.

Hybridization and introgression between Helicoverpa armigera and H. zea: an adaptational bridge

BACKGROUND

Invasion of organisms into new ecosystems is increasingly common, due to the global trade in commodities. One of the most complex post-invasion scenarios occurs when an invasive species is related to a native pest, and even more so when they can hybridize and produce fertile progeny. The global pest Helicoverpa armigera was first detected in Brazil in 2013 and generated a wave of speculations about the possibility of hybridization with the native sister taxon Helicoverpa zea. In the present study, we used genome-wide single nucleotide polymorphisms from field-collected individuals to estimate hybridization between H. armigera and H. zea in different Brazilian agricultural landscapes.

RESULTS

The frequency of hybridization varied from 15 to 30% depending on the statistical analyses. These methods showed more congruence in estimating that hybrids contained approximately 10% mixed ancestry (i.e. introgression) from either species. Hybridization also varied considerably depending on the geographic locations where the sample was collected, forming a 'mosaic' hybrid zone where introgression may be facilitated by environmental and landscape variables. Both landscape composition and bioclimatic variables indicated that maize and soybean cropland are the main factors responsible for high levels of introgression in agricultural landscapes. The impact of multiple H. armigera incursions is reflected in the structured and inbred pattern of genetic diversity.

CONCLUSIONS

Our data showed that the landscape composition and bioclimatic variables influence the introgression rate between H. armigera and H. zea in agricultural areas. Continuous monitoring of the hybridization process in the field is necessary, since agricultural expansion, climatic fluctuations, changing composition of crop species and varieties, and dynamic planting seasons are some factors in South America that could cause a sudden alteration in the introgression rate between Helicoverpa species. Introgression between invasive and native pests can dramatically impact the evolution of host ranges and resistance management.

Global gene expression changes induced by knockout of a protease gene cluster in Helicoverpa armigera with CRISPR/Cas9

Helicoverpa armigera is one of the most serious agricultural insect pests of global importance. It is highly polyphagous and depends on digestive serine proteases to degrade proteins to peptides and to amino acids. H. armigera has evolved adaptive ability to compensate for the inhibition of plant defensive protease inhibitors (PIs) in its diet by overproduction of digestive enzymes. As far as we know, compensation for deletion of serine protease genes has not yet been studied in any herbivorous insect. In this study, we used CRISPR/Cas9 to knock out a cluster of 18 trypsin-like genes in H. armigera. Compared with the wild type SCD strain, activities of the total proteases, trypsins and chymotrypsins were not significantly changed in the gene cluster knockout strain (Tryp-KO). RNA-seq data showed 1492 upregulated and 461 downregulated DEGs in Try-KO. GO function classification and KEGG pathway analyses revealed these differentially expressed genes were enriched for terms related to binding, catalytic activity, metabolic process and signal transduction. In regard to serine protease genes, 35 were upregulated and 12 downregulated in Tryp-KO strain. Our study indicated that H. armigera can compensate for the deleted protease genes by overexpression of other trypsin and chymotrypsin genes in order to maintain its genetic and metabolic robustness. It also suggests that genetic perturbations created by genome editing tools can induce global gene expression changes.

Multiple incursion pathways for Helicoverpa armigera in Brazil show its genetic diversity spreading in a connected world

The Old World cotton bollworm Helicoverpa armigera was first detected in Brazil with subsequent reports from Paraguay, Argentina, Bolivia, and Uruguay. This pattern suggests that the H. armigera spread across the South American continent following incursions into northern/central Brazil, however, this hypothesis has not been tested. Here we compare northern and central Brazilian H. armigera mtDNA COI haplotypes with those from southern Brazil, Uruguay, Argentina, and Paraguay. We infer spatial genetic and gene flow patterns of this dispersive pest in the agricultural landscape of South America. We show that the spatial distribution of H. armigera mtDNA haplotypes and its inferred gene flow patterns in the southwestern region of South America exhibited signatures inconsistent with a single incursion hypothesis. Simulations on spatial distribution patterns show that the detection of rare and/or the absence of dominant mtDNA haplotypes in southern H. armigera populations are inconsistent with genetic signatures observed in northern and central Brazil. Incursions of H. armigera into the New World are therefore likely to have involved independent events in northern/central Brazil, and southern Brazil/Uruguay-Argentina-Paraguay. This study demonstrates the significant biosecurity challenges facing the South American continent, and highlights alternate pathways for introductions of alien species into the New World.

Use of DNA Markers for Grape Phylloxera Population and Evolutionary Genetics: From RAPDs to SSRs and Beyond

Grape phylloxera (Daktulosphaira vitifoliae Fitch) is a major pest of cultivated grapevines (Vitis spp.), occurring in virtually all viticultural regions around the world. Different grape phylloxera strains can be found at varying levels on leaves and roots on both own-rooted plants and in plants grafted onto partially resistant rootstocks. Considering its relevance for the adequate management of the pest in infested vineyards, the analysis of its genetic diversity has received considerable attention from the scientific community in the last decades. Here, we review 25 years of DNA-based molecular markers applied to the analysis of the genetic structure and the reproductive mode of grape phylloxera in its native range and in different introduced regions. The use given to RAPD, AFLP, mtDNA sequencing and microsatellite (SSR) genetic markers for the analysis of grape phylloxera diversity is discussed, and an overview of the main findings obtained after their application to different populations collected in diverse regions all around the world is shown. Lastly, we explore how recent advancements in molecular biology and in modern high throughput genotyping technologies may be applied to better understand grape phylloxera natural diversity at a molecular level.

Invasion origin, rapid population expansion, and the lack of genetic structure of cotton bollworm (Helicoverpa armigera) in the Americas

In 2013, Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) was officially declared as present in Brazil and, after two years, the species was detected in the Caribbean and North America. Information on genetic features and accurate distribution of pests is the basis for agricultural protection policies. Furthermore, such knowledge is imperative to develop control strategies, understand the geographical range, and genetic patterns of this species in the Americas. Here, we carried out the widest sampling of H. armigera in the South American continent and Puerto Rico, after we estimated the diversity, demographic parameters, and genetic structure. The Internal Transcribed Spacer 1 (ITS1) nuclear marker was used to investigate the presence of putative hybrids between H. armigera and H. zea, and they were observed at a frequency of 1.5%. An ABC analysis, based in COI gene fragment, suggested Europe as the origin of South America specimens of H. armigeraand following a movement northward through the Caribbean. Three mtDNA genes and three nDNA markers revealed high genetic diversity distributed without the defined population structure of H. armigera in South America. Most of the genetic variation is within populations with a multidirectional expansion of H. armigera among morphoclimatic regions. High genetic diversity, rapid population expansion, and hybridization have implications for pest management since they suggest that adaptive alleles are spread through wide areas in South America that favor rapid local adaptation of H. armigera to new and disturbed environments (e.g., in agricultural areas).

Mitochondrial DNA genomes of five major Helicoverpa pest species from the Old and New Worlds (Lepidoptera: Noctuidae)

Five species of noctuid moths, Helicoverpa armigera, H. punctigera, H. assulta, H. zea, and H. gelotopoeon, are major agricultural pests inhabiting various and often overlapping global distributions. Visual identification of these species requires a great deal of expertise and misidentification can have repercussions for pest management and agricultural biosecurity. Here, we report on the complete mitochondrial genomes of H. assulta assulta and H. assulta afra, H. gelotopoeon, H. punctigera, H. zea, and H. armigera armigera and H. armigera conferta' assembled from high-throughput sequencing data. This study significantly increases the mitogenome resources for these five agricultural pests with sequences assembled from across different continents, including an H. armigera individual collected from an invasive population in Brazil. We infer the phylogenetic relationships of these five Helicoverpa species based on the 13 mitochondrial DNA protein-coding genes (PCG's) and show that two publicly available mitogenomes of H. assulta (KP015198 and KR149448) have been misidentified or incorrectly assembled. We further consolidate existing PCR-RFLP methods to cover all five Helicoverpa pest species, providing an updated method that will contribute to species differentiation and to future monitoring efforts of Helicoverpa pest species across different continents. We discuss the value of Helicoverpa mitogenomes to assist with species identification in view of the context of the rapid spread of H. armigera in the New World. With this work, we provide the molecular resources necessary for future studies of the evolutionary history and ecology of these species.

Going global - genomic insights into insect invasions

The spread of invasive insect pests is becoming an increasing problem for agriculture globally. We discuss a number of invasive insects, already of major economic significance that have recently expanded their range to become truly global threats. These include the noctuid moths Helicoverpa and Spodoptera, whose caterpillars have long been among the worst pests in their native Old and New World habitats, respectively, and the whitefly Bemisia, a major vector of plant virus diseases. Importantly, genomic resources for these species have recently become available, allowing research to move beyond the restrictions imposed by earlier approaches limited to a single or few mitochondrial and nuclear markers, to employ genome-wide genotyping and resequencing protocols. These studies have shown hybridisation within the various species complexes, identified regions under selection in agricultural environments, and enable monitoring of genes important as biosecurity risks through introgression into established populations free of the genes. In all cases studied, global trade has emerged as the probable cause of insect spread, making it ever more important that biosecurity protocols and agencies work with researchers to make the most effective use of emerging genomic resources and tools.

Movement Ecology of Pest Helicoverpa: Implications for Ongoing Spread

The recent introduction and spread of Helicoverpa armigera throughout South America highlight the invasiveness and adaptability of moths in the Helicoverpa genus. Long-range movement in three key members, H. armigera, H. zea, and H. punctigera, occurs by migration and international trade. These movements facilitate high population admixture and genetic diversity, with important economic, biosecurity, and control implications in today's agricultural landscape. This is particularly true for the spread of resistance alleles to transgenic crops expressing Bacillus thuringiensis (Bt) toxins that are planted over vast areas to suppress Helicoverpa spp. The ability to track long-distance movement through radar technology, population genetic markers, and/or long-distance dispersal modeling has advanced in recent years, yet we still know relatively little about the population trajectories or migratory routes in Helicoverpa spp. Here, we consider how experimental and theoretical approaches can be integrated to fill key knowledge gaps and assist management practices.

Population differentiation between Australian and Chinese Helicoverpa armigera occurs in distinct blocks on the Z-chromosome

Over the last 40 years, many types of population genetic markers have been used to assess the population structure of the pest moth species Helicoverpa armigera. While this species is highly vagile, there is evidence of inter-continental population structure. Here, we examine Z-chromosome molecular markers within and between Chinese and Australian populations. Using 1352 polymorphic sites from 40 Z-linked loci, we compared two Chinese populations of moths separated by 700 km and found virtually no population structure (n = 41 and n = 54, with <1% of variation discriminating between populations). The levels of nucleotide diversity within these populations were consistent with previous estimates from introns in Z-linked genes of Australian samples (π = 0.028 vs. 0.03). Furthermore, all loci surveyed in these Chinese populations showed a skew toward rare variants, with ten loci having a significant Tajima's D statistic, suggesting that this species could have undergone a population expansion. Eight of the 40 loci had been examined in a previous study of Australian moths, of which six revealed very little inter-continental population structure. However, the two markers associated with the Cyp303a1 locus that has previously been proposed to be a target of a selective sweep, exhibited allele structuring between countries. Using a separate dataset of 19 Australian and four Chinese moths, we scanned the molecular variation distributed across the entire Z-chromosome and found distinct blocks of differentiation that include the region containing Cyp303a1. We recommend some of these loci join those associated with insecticide resistance to form a set of genes best suited to analyzing population structure in this global pest.

Multiple recombination events between two cytochrome P450 loci contribute to global pyrethroid resistance in Helicoverpa armigera

The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most serious insect pest species to evolve resistance against many insecticides from different chemical classes. This species has evolved resistance to the pyrethroid insecticides across its native range and is becoming a truly global pest after establishing in South America and having been recently recorded in North America. A chimeric cytochrome P450 gene, CYP337B3, has been identified as a resistance mechanism for resistance to fenvalerate and cypermethrin. Here we show that this resistance mechanism is common around the world with at least eight different alleles. It is present in South America and has probably introgressed into its closely related native sibling species, Helicoverpa zea. The different alleles of CYP337B3 are likely to have arisen independently in different geographic locations from selection on existing diversity. The alleles found in Brazil are those most commonly found in Asia, suggesting a potential origin for the incursion of H. armigera into the Americas.

Biotic Potential and Life Table of Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) from Three Brazilian Regions

This study aimed to evaluate the biotic potential and life table of individuals of Helicoverpa armigera (Hübner) from different host plants (citrus, corn, and cotton) and Brazilian states (São Paulo, Distrito Federal, and Bahia) in artificial diet, under laboratory conditions (25 ± 1°C, 70 ± 10% RH, 14 h photophase). The longevity, pre-, post- and oviposition periods, fecundity, and fertility of 15 mating pairs per origin were evaluated. We also compared the reproductive parameters of each group of insects (São Paulo (SP), Distrito Federal (DF), and Bahia (BA)), including the net reproductive rate (Ro), mean generation time (T), intrinsic rate of increase (r _m), and finite rate of increase (λ). Microsatellite analysis from individuals collected in different locations and host plants did not show differences among the parental insects. It was verified that parental progeny collected in cotton fields from Bahia had a higher biotic potential, a higher reproductive rate (Ro), and a better fecundity compared to the insects from remaining regions. The life table charts indicate that the highest values for the reproductive parameters of the Bahia progeny are associated with higher specific fertility, particularly in early adulthood. The greatest biotic potential of the Bahia progeny may be due to increased selection pressure from the insecticide used (organophosphate and pyrethroid) on cotton crops compared to that of other crops, as well due to the massive adoption of Bt cotton-producing areas of that state from 2013 outbreaks.

Epistasis confers resistance to Bt toxin Cry1Ac in the cotton bollworm

Evolution of resistance by insect pests reduces the benefits of extensively cultivated transgenic crops that produce insecticidal proteins from Bacillus thuringiensis (Bt). Previous work showed that resistance to Bt toxin Cry1Ac, which is produced by transgenic cotton, can be conferred by mutations disrupting a cadherin protein that binds this Bt toxin in the larval midgut. However, the potential for epistatic interactions between the cadherin gene and other genes has received little attention. Here, we report evidence of epistasis conferring resistance to Cry1Ac in the cotton bollworm, Helicoverpa armigera, one of the world's most devastating crop pests. Resistance to Cry1Ac in strain LF256 originated from a field-captured male and was autosomal, recessive, and 220-fold relative to susceptible strain SCD. We conducted complementation tests for allelism by crossing LF256 with a strain in which resistance to Cry1Ac is conferred by a recessive allele at the cadherin locus HaCad. The resulting F1 offspring were resistant, suggesting that resistance to Cry1Ac in LF256 is also conferred by resistance alleles at this locus. However, the HaCad amino acid sequence in LF256 lacked insertions and deletions, and did not differ consistently between LF256 and a susceptible strain. In addition, most of the cadherin alleles in LF256 were not derived from the field-captured male. Moreover, Cry1Ac resistance was not genetically linked with the HaCad locus in LF256. Furthermore, LF256 and the susceptible strain were similar in levels of HaCad transcript, cadherin protein, and binding of Cry1Ac to cadherin. Overall, the results imply that epistasis between HaCad and an unknown second locus in LF256 yielded the observed resistance in the F1 progeny from the complementation test. The observed epistasis has important implications for interpreting results of the F1 screen used widely to monitor and analyze resistance, as well as the potential to accelerate evolution of resistance.

Detection of sister-species in invasive populations of the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) from Uganda

The fall armyworm (FAW) Spodoptera frugiperda (J. E. Smith) is a species native to the Americas. This polyphagous lepidopteran pest was first reported in Nigeria and the Democratic Republic of São Tomé and Principe in 2016, but its presence in eastern Africa has not been confirmed via molecular characterisation. In this study, FAW specimens from western and central Uganda were identified based on the partial mtDNA COI gene sequences, with mtDNA COI haplotypes matching those identified in Nigeria and São Tomé. In this study, we sequence an additional partial mtDNA Cyt b gene and also the partial mtDNA COIII gene in Ugandan FAW samples. We detected identical mitochondrial DNA haplotypes for both the mtDNA Cyt b and COI partial genes, while combining the mtDNA COI/Cyt b haplotypes and mtDNA COIII haplotypes enabled a new maternal lineage in the Ugandan corn-preferred FAW samples to be identified. Our results suggested that the African incursions of S. frugiperda involved at least three maternal lineages. Recent full genome, phylogenetic and microsatellite analyses provided evidence to support S. frugiperda as likely consisted of two sympatric sister species known as the corn-preferred and rice-preferred strains. In our Ugandan FAW populations, we identified the presence of mtDNA haplotypes representative of both sister species. It is not known if both FAW sister species were originally introduced together or separately, and whether they have since spread as a single population. Further analyses of additional specimens originally collected from São Tomé, Nigeria and throughout Africa would be required to clarify this issue. Importantly, our finding showed that the genetic diversity of the African corn-preferred FAW species is higher than previously reported. This potentially contributed to the success of FAW establishment in Africa. Furthermore, with the additional maternal lineages detected, there is likely an increase in paternal lineages, thereby increasing the diversity of the African FAW population. Knowledge of the FAW genetic diversity will be needed to assess the risks of introducing Bt-resistance traits and to understand the FAW incursion pathways into the Old World and its potential onward spread. The agricultural implications of the presence of two evolutionary divergent FAW lineages (the corn and the rice lineage) in the African continent are further considered and discussed.