Temporal Variation in Genetic Composition of Migratory Helicoverpa Zea in Peripheral Populations

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

A Chromosome-Scale Genome Assembly of a Helicoverpa zea Strain Resistant to Bacillus thuringiensis Cry1Ac Insecticidal Protein

Helicoverpa zea (Lepidoptera: Noctuidae) is an insect pest of major cultivated crops in North and South America. The species has adapted to different host plants and developed resistance to several insecticidal agents, including Bacillus thuringiensis (Bt) insecticidal proteins in transgenic cotton and maize. Helicoverpa zea populations persist year-round in tropical and subtropical regions, but seasonal migrations into temperate zones increase the geographic range of associated crop damage. To better understand the genetic basis of these physiological and ecological characteristics, we generated a high-quality chromosome-level assembly for a single H. zea male from Bt-resistant strain, HzStark_Cry1AcR. Hi-C data were used to scaffold an initial 375.2 Mb contig assembly into 30 autosomes and the Z sex chromosome (scaffold N50 = 12.8 Mb and L50 = 14). The scaffolded assembly was error-corrected with a novel pipeline, polishCLR. The mitochondrial genome was assembled through an improved pipeline and annotated. Assessment of this genome assembly indicated 98.8% of the Lepidopteran Benchmark Universal Single-Copy Ortholog set were complete (98.5% as complete single copy). Repetitive elements comprised approximately 29.5% of the assembly with the plurality (11.2%) classified as retroelements. This chromosome-scale reference assembly for H. zea, ilHelZeax1.1, will facilitate future research to evaluate and enhance sustainable crop production practices.

Temporal sampling and network analysis reveal rapid population turnover and dynamic migration pattern in overwintering regions of a cosmopolitan pest

Genetic makeup of insect pest is informative for source-sink dynamics, spreading of insecticide resistant genes, and effective management. However, collecting samples from field populations without considering temporal resolution and calculating parameters related to historical gene flow may not capture contemporary genetic pattern and metapopulation dynamics of highly dispersive pests. Plutella xylostella (L.), the most widely distributed Lepidopteran pest that developed resistance to almost all current insecticides, migrates heterogeneously across space and time. To investigate its real-time genetic pattern and dynamics, we executed four samplings over two consecutive years across Southern China and Southeast Asia, and constructed population network based on contemporary gene flow. Across 48 populations, genetic structure analysis identified two differentiated insect swarms, of which the one with higher genetic variation was replaced by the other over time. We further inferred gene flow by estimation of kinship relationship and constructed migration network in each sampling time. Interestingly, we found mean migration distance at around 1,000 km. Such distance might have contributed to the formation of step-stone migration and migration circuit over large geographical scale. Probing network clustering across sampling times, we found a dynamic P. xylostella metapopulation with more active migration in spring than in winter, and identified a consistent pattern that some regions are sources (e.g., Yunnan in China, Myanmar and Vietnam) while several others are sinks (e.g., Guangdong and Fujian in China) over 2 years. Rapid turnover of insect swarms and highly dynamic metapopulation highlight the importance of temporal sampling and network analysis in investigation of source-sink relationships and thus effective pest management of P. xylostella, and other highly dispersive insect pests.

Landscape-level variation in Bt crops predict Helicoverpa zea (Lepidoptera: Noctuidae) resistance in cotton agroecosystems

BACKGROUND

Helicoverpa zea (Boddie) damage to Bt cotton and maize has increased as a result of widespread Bt resistance across the USA Cotton Belt. Our objective was to link Bt crop production patterns to cotton damage through a series of spatial and temporal surveys of commercial fields to understand how Bt crop production relates to greater than expected H. zea damage to Bt cotton. To do this, we assembled longitudinal cotton damage data that spanned the Bt adoption period, collected cotton damage data since Bt resistance has been detected, and estimated local population susceptibility using replicated on-farm studies that included all Bt pyramids marketed in cotton.

RESULTS

Significant year effects of H. zea damage frequency in commercial cotton were observed throughout the Bt adoption period, with a recent damage increase after 2012. Landscape-level Bt crop production intensity over time was positively associated with the risk of H. zea damage in two- and three-toxin pyramided Bt cotton. Helicoverpa zea damage also varied across Bt toxin types in spatially replicated on-farm studies.

CONCLUSIONS

Landscape-level predictors of H. zea damage in Bt cotton can be used to identify heightened Bt resistance risk areas and serves as a model to understand factors that drive pest resistance evolution to Bt toxins in the southeastern United States. These results provide a framework for more effective insect resistance management strategies to be used in combination with conventional pest management practices that improve Bt trait durability while minimizing the environmental footprint of row crop agriculture. © 2021 Society of Chemical Industry. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Sweet Corn Sentinel Monitoring for Lepidopteran Field-Evolved Resistance to Bt Toxins

As part of an insect resistance management plan to preserve Bt transgenic technology, annual monitoring of target pests is mandated to detect susceptibility changes to Bt toxins. Currently Helicoverpa zea (Boddie) monitoring involves investigating unexpected injury in Bt crop fields and collecting larvae from non-Bt host plants for laboratory diet bioassays to determine mortality responses to diagnostic concentrations of Bt toxins. To date, this monitoring approach has not detected any significant change from the known range of baseline susceptibility to Bt toxins, yet practical field-evolved resistance in H. zea populations and numerous occurrences of unexpected injury occur in Bt crops. In this study, we implemented a network of 73 sentinel sweet corn trials, spanning 16 U.S. states and 4 Canadian provinces, for monitoring changes in H. zea susceptibility to Cry and Vip3A toxins by measuring differences in ear damage and larval infestations between isogenic pairs of non-Bt and Bt hybrids over three years. This approach can monitor susceptibility changes and regional differences in other ear-feeding lepidopteran pests. Temporal changes in the field efficacy of each toxin were evidenced by comparing our current results with earlier published studies, including baseline data for each Bt trait when first commercialized. Changes in amount of ear damage showed significant increases in H. zea resistance to Cry toxins and possibly lower susceptibility to Vip3a. Our findings demonstrate that the sentinel plot approach as an in-field screen can effectively monitor phenotypic resistance and document field-evolved resistance in target pest populations, improving resistance monitoring for Bt crops.