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

Accurate identification of Helicoverpa armigera-Helicoverpa zea hybrids using genome admixture analysis: implications for genomic surveillance

Helicoverpa armigera, the cotton bollworm moth, is one of the world's most important crop pests, and is spreading throughout the New World from its original range in the Old World. In Brazil, invasive H. armigera has been reported to hybridize with local populations of Helicoverpa zea. The correct identification of H. armigera-H. zea hybrids is important in understanding the origin, spread and future outlook for New World regions that are affected by outbreaks, given that hybridization can potentially facilitate H. zea pesticide resistance and host plant range via introgression of H. armigera genes. Here, we present a genome admixture analysis of high quality genome sequences generated from two H. armigera-H. zea F1 hybrids generated in two different labs. Our admixture pipeline predicts 48.8% and 48.9% H. armigera for the two F1 hybrids, confirming its accuracy. Genome sequences from five H. zea and one H. armigera that were generated as part of the study show no evidence of hybridization. Interestingly, we show that four H. zea genomes generated from a previous study are predicted to possess a proportion of H. armigera genetic material. Using unsupervised clustering to identify non-hybridized H. armigera and H. zea genomes, 8511 ancestry informative markers (AIMs) were identified. Their relative frequencies are consistent with a minor H. armigera component in the four genomes, however its origin remains to be established. We show that the size and quality of genomic reference datasets are critical for accurate hybridization prediction. Consequently, we discuss potential pitfalls in genome admixture analysis of H. armigera-H. zea hybrids, and suggest measures that will improve such analyses.

Field-based recombinase polymerase amplification and lab-based qPCR assays for detection of Helicoverpa armigera

Helicoverpa armigera (Hübner) is a major crop pest native to Europe, Asia, Australia, and Africa which has recently invaded South America and has caused billions of dollars in agricultural losses. Because of challenges in differentiating between H. armigera and Helicoverpa zea (Boddie), a closely related species native to North and South America, genetic tests have previously been developed to detect H. armigera DNA in pooled samples of moth legs. In this study, a field-based recombinase polymerase amplification (RPA) assay using a lateral flow strip and a qPCR melt curve assay were developed for specific detection of H. armigera DNA in pooled moth samples. In addition, a crude DNA extraction protocol for whole moths was developed to allow rapid preparation of DNA samples. The RPA field test was able to detect ≥ 10 pg of purified H. armigera DNA and the crude DNA of one H. armigera sample in a background of 999 H. zea equivalents. The qPCR assay was able to detect ≥ 100 fg of purified H. armigera DNA and the crude DNA of one H. armigera sample in a background of up to 99,999 H. zea equivalents. Both RPA and qPCR assays detected H. armigera in the crude DNA extracted in the field from a pool of one H. armigera moth and 999 H. zea moths. These newly developed molecular assays to detect H. armigera will contribute to large-scale surveillance programs of H. armigera.

Complex multiple introductions drive fall armyworm invasions into Asia and Australia

The fall armyworm (FAW) Spodoptera frugiperda is thought to have undergone a rapid 'west-to-east' spread since 2016 when it was first identified in western Africa. Between 2018 and 2020, it was recorded from South Asia (SA), Southeast Asia (SEA), East Asia (EA), and Pacific/Australia (PA). Population genomic analyses enabled the understanding of pathways, population sources, and gene flow in this notorious agricultural pest species. Using neutral single nucleotide polymorphic (SNP) DNA markers, we detected genome introgression that suggested most populations in this study were overwhelmingly C- and R-strain hybrids (n = 252/262). SNP and mitochondrial DNA markers identified multiple introductions that were most parsimoniously explained by anthropogenic-assisted spread, i.e., associated with international trade of live/fresh plants and plant products, and involved 'bridgehead populations' in countries to enable successful pest establishment in neighbouring countries. Distinct population genomic signatures between Myanmar and China do not support the 'African origin spread' nor the 'Myanmar source population to China' hypotheses. Significant genetic differentiation between populations from different Australian states supported multiple pathways involving distinct SEA populations. Our study identified Asia as a biosecurity hotspot and a FAW genetic melting pot, and demonstrated the use of genome analysis to disentangle preventable human-assisted pest introductions from unpreventable natural pest spread.

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.

Coocclusion of Helicoverpa armigera Single Nucleopolyhedrovirus (HearSNPV) and Helicoverpa armigera Multiple Nucleopolyhedrovirus (HearMNPV): Pathogenicity and Stability in Homologous and Heterologous Hosts

Helicoverpa armigera single nucleopolyhedrovirus (HearSNPV) is a virulent pathogen of lepidopterans in the genera Heliothis and Helicoverpa, whereas Helicoverpa armigera multiple nucleopolyhedrovirus (HearSNPV) is a different virus species with a broader host range. This study aimed to examine the consequences of coocclusion of HearSNPV and HearMNPV on the pathogenicity, stability and host range of mixed-virus occlusion bodies (OBs). HearSNPV OBs were approximately 6-fold more pathogenic than HearMNPV OBs, showed faster killing by approximately 13 h, and were approximately 45% more productive in terms of OB production per larva. For coocclusion, H. armigera larvae were first inoculated with HearMNPV OBs and subsequently inoculated with HearSNPV OBs at intervals of 0–72 h after the initial inoculation. When the interval between inoculations was 12–24 h, OBs collected from virus-killed insects were found to comprise 41–57% of HearSNPV genomes, but the prevalence of HearSNPV genomes was greatly reduced (3–4%) at later time points. Quantitative PCR (qPCR) analysis revealed the presence of HearSNPV genomes in a small fraction of multinucleocapsid ODVs representing 0.47–0.88% of the genomes quantified in ODV samples, indicating that both viruses had replicated in coinfected host cells. End-point dilution assays on ODVs from cooccluded mixed-virus OBs confirmed the presence of both viruses in 41.9–55.6% of wells that were predicted to have been infected by a single ODV. A control experiment indicated that this result was unlikely to be due to the adhesion of HearSNPV ODVs to HearMNPV ODVs or accidental contamination during ODV band extraction. Therefore, the disparity between the qPCR and end-point dilution estimates of the prevalence of mixed-virus ODVs likely reflected virus-specific differences in replication efficiency in cell culture and the higher infectivity of pseudotyped ODVs that were produced in coinfected parental cells. Bioassays on H. armigera, Spodoptera frugiperda and Mamestra brassicae larvae revealed that mixed-virus OBs were capable of infecting heterologous hosts, but relative potency values largely reflected the proportion of HearMNPV present in each mixed-virus preparation. The cooccluded mixtures were unstable in serial passage; HearSNPV rapidly dominated during passage in H. armigera whereas HearMNPV rapidly dominated during passage in the heterologous hosts. We conclude that mixed-virus coocclusion technology may be useful for producing precise mixtures of viruses with host range properties suitable for the control of complexes of lepidopteran pests in particular crops, although this requires validation by field testing.

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.

Fruit Herbivory Alters Plant Electrome: Evidence for Fruit-Shoot Long-Distance Electrical Signaling in Tomato Plants

The electrical activity of tomato plants subjected to fruit herbivory was investigated. The study aimed to test the hypothesis that tomato fruits transmit long-distance electrical signals to the shoot when subjected to herbivory. For such, time series classification by machine learning techniques and analyses related to the oxidative response were employed. Tomato plants (cv. “Micro-Tom”) were placed into a Faraday's cage and an electrode pair was inserted in the fruit's peduncle. Helicoverpa armigera caterpillars were placed on the fruit (either green and ripe) for 24 h. The time series were recorded before and after the fruit's exposure of the caterpillars. The plant material for chemical analyses was collected 24 and 48 h after the end of the acquisition of electrophysiological data. The time series were analyzed by the following techniques: Fast Fourier Transform (FFT), Wavelet Transform, Power Spectral Density (PSD), and Approximate Entropy. The following features from FFT, PSD, and Wavelet Transform were used for PCA (Principal Component Analysis): average, maximum and minimum value, variance, skewness, and kurtosis. Additionally, these features were used in Machine Learning (ML) analyses for looking for classifiable patterns between tomato plants before and after fruit herbivory. Also, we compared the electrome before and after herbivory in the green and ripe fruits. To evaluate an oxidative response in different organs, hydrogen peroxide, superoxide anion, catalase, ascorbate peroxidase, guaiacol peroxidase, and superoxide dismutase activity were evaluated in fruit and leaves. The results show with 90% of accuracy that the electrome registered in the fruit's peduncle before herbivory is different from the electrome during predation on the fruits. Interestingly, there was also a sharp difference in the electrome of the green and ripe fruits' peduncles before, but not during, the herbivory, which demonstrates that the signals generated by the herbivory stand over the others. Biochemical analysis showed that herbivory in the fruit triggered an oxidative response in other parts of the plant. Here, we demonstrate that the fruit perceives biotic stimuli and transmits electrical signals to the shoot of tomato plants. This study raises new possibilities for studies involving electrical signals in signaling and systemic response, as well as for the applicability of ML to classify electrophysiological data and its use in early diagnosis.

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.

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.

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.

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.