Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species

Metabolism of Furanocoumarins by Three Recombinant CYP9A Proteins From the Polyphagous Cotton Bollworm Helicoverpa armigera

Furanocoumarins are a class of chemical compounds with phototoxic properties. For herbivores, efficient detoxification of such defense compounds is the prerequisite to feed successfully on furanocoumarin-containing plants. The cotton bollworm Helicoverpa armigera is a very important polyphagous pest in agriculture, but how it copes with toxic furanocoumarins in some of its host plants is not well understood. Given that cytochrome P450s are well known for their capacity in xenobiotic metabolism, this study attempted to explore the potential roles of cytochrome P450s in furanocoumarin transformation in this pest. Our data showed that two linear structures (psoralen and xanthotoxin) could be metabolized by three recombinant CYP9A enzymes, but no detectable depletion was observed for the linear one with the 8-dimethylallyloxy substituent on the coumarin moiety (imperatorin) and the angular furanocoumarin (angelicin). Initial epoxidation of the double bond connecting C2' and C3' of the furano ring following by cleavage of the epoxidated furan ring, leading to the formation of more soluble, less reactive and nonphotosensitizing metabolites, was identified as a common mechanism of linear furanocoumarin metabolism using a quadrupole/time-of-flight (Q-TOF) mass spectrometry interfaced with a high performance liquid chromatography (HPLC) system. Our data demonstrated that multiple P450s were involved in the detoxification of linear furanocoumarins in the cotton bollworm. These findings contribute to a better understanding of the biochemical basis of adaptation to plant defense chemicals in this economically important pest.

The effects of Cl- channel inhibitors and pyrethroid insecticides on calcium-activated chloride channels in neurons of Helicoverpa armigera

TMEM16A, a member of the Transmembrane protein 16 family, serves as the molecular basis for calcium activated chloride channels (CaCCs). We use RT-PCR to demonstrate the expression of TMEM16A in the neurons of Helicoverpa armigera, and record the CaCCs current of acute isolated neurons of H. armigera for the first time using patch clamp technology. In order to screen effective inhibitors of calcium-activated chloride channels, the inhibitory effects of four chloride channel inhibitors, CaCCinh-A01, NPPB, DIDS, and SITS, on CaCCs were compared. The inhibitory effects of the four inhibitors on the outward current of CaCCs were CaCCinh-A01 (10 μM, 56.31 %), NPPB (200 μM, 43.69 %), SITS (1 mM, 12.41 %) and DIDS (1 mM, 13.29 %). Among these inhibitors, CaCCinh-A01 demonstrated the highest efficacy as a blocker. To further explore whether calcium channel proteins can serve as potential targets of pyrethroids, we compared the effects of (type I) tefluthrin and (type II) deltamethrin on CaCCs. 10 μM and 100 μM tefluthrin can stimulate a large tail current in CaCCs, prolonging their deactivation time by 10.44 ms and 31.49 ms, and the V0.5 shifted in the hyperpolarization by 2-8 mV. Then, deltamethrin had no obvious effect on the deactivation and activation of CaCCs. Therefore, CaCCs of H. armigera can be used as a potential target of pyrethroids, but type I and type II pyrethroids have different effects on CaCCs.

Non-CG DNA methylation marks the transition from pupa to adult in Helicoverpa armigera

DNA methylation in insects is generally low in abundance, and its role is not well understood. It is often localised in protein coding regions and associated with the expression of 'housekeeping' genes. Few studies have explored DNA methylation dynamics during lifecycle stage transitions in holometabolous (metamorphosing) insects. Using targeted mass spectrometry, we have found a significant difference in global DNA methylation levels between larvae, pupae and adults of Helicoverpa armigera (Lepidoptera: Noctuidae) Hübner, a polyphagous pest of agricultural importance. Whole-genome bisulfite sequencing confirmed these observations and pointed to non-CG context being the primary explanation for the difference observed between pupa and adult. Non-CG methylation was enriched in genes specific to various signalling pathways (Hippo signalling, Hedgehog signalling and mitogen-activated protein kinase (MAPK) signalling) and ATP-dependent chromatin remodelling. Understanding the function of this epigenetic mark could be a target in future studies focusing on integrated pest 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.

Tyramine-Mediated Hyperactivity Modulates the Dietary Habits in Helicoverpa armigera

Helicoverpa armigera exhibits extensive variability in feeding habits and food selection. Neuronal regulation of H. armigera feeding behavior is primarily influenced by biogenic amines such as Tyramine (TA) and Octopamine (OA). The molecular responses of H. armigera to dietary challenges in the presence of TA or OA have yet to be studied. This investigation dissects the impact of OA and TA on H. armigera feeding choices and behaviors under non-host nutritional stress. It has been observed that feeding behavior remains unaltered during the exogenous administration of OA and TA through an artificial diet (AD). Ingestion of higher OA or TA concentrations leads to increased mortality. OA and TA treatment in combination with host and non-host diets results in the induction of feeding and higher locomotion toward food, particularly in the case of TA treatment. Increased expression of markers, prominin-like, and tachykinin-related peptide receptor-like transcripts further assessed increased locomotion activity. Insects subjected to a non-host diet with TA treatment exhibited increased feeding and overexpression of the feeding indicator, the Neuropeptide F receptor, and the feeding regulator, Sulfakinin, compared with other conditions. Expression of sensation and biogenic amine synthesis genesis elevated in insects fed a non-host diet in combination with OA or TA. Metabolomics analysis revealed a decreased concentration of the feeding behavior elicitor, dopamine, in insects fed a non-host diet containing TA. This work highlights the complex interplay between biogenic amine functions during dietary stress and suggests the role of tyramine in feeding promotion under stressed conditions.

Homology-based characterization of the cis-regulatory elements modulate flavone induction of CYP321A1 in Helicoverpa armigera

Xenobiotic response element (XRE) to flavone was the cis- regulatory elements that mediates the induction of the allelochemical-metabolizing CYP321A1 gene from Helicoverpa zea. However, it was unknown whether the XRE-Fla element existed in other species. Recently we have identified and cloned the CYP321A1 gene with promoter region in a related species, Helicoverpa armigera. Sequence similarity of two orthologous CYP321A1 genes was 97.27%, but the promoter sequence similarity was only 56.32%. Sequence alignment showed the XRE-Fla like element owns three mutations in H. armigera compared with H. zea. Progressive 5' deletions and internal mutation indicated that H. armigera XRE-Fla was the essential element of CYP321A1 gene in response to flavone. XRE-Fla mutations and EMSA analysis confirmed that the H. armigera XRE-Fla element binding factor was stronger than H. zea. The findings indicate the XRE element mutations mainly contribute to the differences between the flavone-induced expressions of two CYP321A1 genes, which improve the flexibility and adaptability for allelochemical response of H. armigera.

The role of aquaporins in osmotic cell lysis induced by Bacillus thuringiensis Cry1Ac toxin in Helicoverpa armigera

The insecticidal crystalline (Cry) and vegetative insecticidal (Vip) proteins derived from Bacillus thuringiensis (Bt) are used globally to manage insect pests, including the cotton bollworm, Helicoverpa armigera, one of the world's most damaging agricultural pests. Cry proteins bind to the ATP-binding cassette transporter C2 (ABCC2) receptor on the membrane surface of larval midgut cells, resulting in Cry toxin pores, and ultimately leading to cell swelling and/or lysis. Insect aquaporin (AQP) proteins within the membranes of larval midgut cells are proposed to allow the rapid influx of water into enterocytes following the osmotic imbalance triggered by the formation of Cry toxin pores. Here, we examined the involvement of H. armigera AQPs in Cry1Ac-induced osmotic cell swelling. We identified and characterized eight H. armigera AQPs and demonstrated that five are functional water channel proteins. Three of these (HaDrip1, HaPrip, and HaEglp1) were found to be expressed in the larval midgut. Xenopus laevis oocytes co-expressing the known Cry1Ac receptor HaABCC2 and each of the three HaAQPs displayed abnormal morphology and were lysed following exposure to Cry1Ac, suggesting a rapid influx of water was induced after Cry1Ac pore formation. In contrast, oocytes producing either HaABCC2 or HaAQP alone failed to swell or lyse after treatment with Cry1Ac, implying that both Cry1Ac pore formation and HaAQP function are needed for osmotic cell swelling. However, CRISPR/Cas9-mediated knockout of any one of the three HaAQP genes failed to cause significant changes in susceptibility to the Bt toxins Cry1Ac, Cry2Ab, or Vip3Aa. Our findings suggest that the multiple HaAQPs produced in larval midgut cells compensate for each other in allowing for the rapid influx of water in H. armigera midgut cells following Cry toxin pore formation, and that mutations affecting a single HaAQP are unlikely to confer resistance to Bt proteins.

The genome sequence of the Cotton Bollworm moth, Helicoverpa armigera (Hübner, 1808)

We present a genome assembly from an adult female Cotton Bollworm moth, Helicoverpa armigera (Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence has a total length of 362.20 megabases. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the W and Z sex chromosomes. The mitochondrial genome has also been assembled and is 15.36 kilobases in length.

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.

A Female-Biased Chemosensory Protein PxutCSP19 in the Antennae of Papilio xuthus Tuned to Host Volatiles and Insecticides

Chemosensory protein (CSP) genes significantly enriched in the female antennae are potential molecular candidates for mediating female oviposition behaviors. In this study, we presented the interaction mechanisms of a female-antenna-biased PxutCSP19 in Papilio xuthus to 47 host volatiles, four biopesticides and 24 synthetic insecticides. Using a bioinformatics-based homology search, 22 genes orthologous to PxutCSP19 were identified from 22 other Papilio butterflies with high sequence identities to each other (73.20~98.72%). Multiple alignment analyses revealed a particularly extended N-terminus of Papilio CSP19s (an average of 154 residues) compared to insects' typical CSPs (approximately 120 residues). The expression profiles indicated that PxutCSP19 was significantly enriched in the female antennae, with a 31.81-fold difference relative to the male antennae. In ligand-binding assays, PxutCSP19 could strongly bind six host odorants with high affinities, ranging from dissociation constant (Ki) values of 20.44 ± 0.64 μM to 22.71 ± 0.73 μM. Notably, this protein was tuned to a monoterpenoid alcohol, linalool, which generally existed in the Rutaceae plants and elicited electrophysiological and behavioral activities of the swallowtail butterfly. On the other hand, PxutCSP19 was also capable of binding eight insecticides with stronger binding abilities (Ki < 12 μM) compared to host odorants. When an extended N-terminal region of PxutCSP19 was truncated into two different proteins, they did not significantly affect the binding of PxutCSP19 to ligands with high affinities, suggesting that this extended N-terminal sequences were not involved in the specificity of ligand recognition. Altogether, our study sheds light on the putative roles of PxutCSP19 enriched in the female antennae of P. xuthus in the perception of host volatiles and the sequestering of insecticides, and it complements the knowledge of butterfly CSPs in olfaction and insecticide resistance.

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.

Comparative Genomics and the Salivary Transcriptome of the Redbanded Stink Bug Shed Light on Its High Damage Potential to Soybean

The redbanded stink bug, Piezodorus guildinii (Westwood) (Hemiptera: Pentatomidae), is a significant soybean pest in the Americas, which inflicts more physical damage on soybean than other native stink bugs. Studies suggest that its heightened impact is attributed to the aggressive digestive properties of its saliva. Despite its agricultural importance, the factors driving its greater ability to degrade plant tissues have remained unexplored in a genomic evolutionary context. In this study, we hypothesized that lineage-specific gene family expansions have increased the copy number of digestive genes expressed in the salivary glands. To investigate this, we annotated a previously published genome assembly of the redbanded stink bug, performed a comparative genomic analysis on 11 hemipteran species, and reconstructed patterns of gene duplication, gain, and loss in the redbanded stink bug. We also performed RNA-seq on the redbanded stink bug's salivary tissues, along with the rest of the body without salivary glands. We identified hundreds of differentially expressed salivary genes, including a subset lost in other stink bug lineages, but retained and expressed in the redbanded stink bug's salivary glands. These genes were significantly enriched with protein families involved in proteolysis, potentially explaining the redbanded stink bug's heightened damage to soybeans. Contrary to our hypothesis, we found no support for an enrichment of duplicated digestive genes that are also differentially expressed in the salivary glands of the redbanded stink bug. Nonetheless, these results provide insight into the evolution of this important crop pest, establishing a link between its genomic history and its agriculturally important physiology.

Genome-wide identification and phylogenetic analysis of the tetraspanin gene family in lepidopteran insects and expression profiling analysis in Helicoverpa armigera

The tetraspanin gene family encodes cell-surface proteins that span the membrane 4 times and play critical roles in a wide range of biological processes across numerous organisms. Recent findings highlight the involvement of a tetraspanin of the lepidopteran pest Helicoverpa armigera in resistance to Bacillus thuringiensis Cry insecticidal proteins, which are extensively used in transgenic crops. Thus, a better understanding of lepidopteran tetraspanins is urgently needed. In the current study, genome scanning in 10 lepidopteran species identified a total of 283 sequences encoding potential tetraspanins. Based on conserved cysteine patterns in the large extracellular loop and their phylogenetic relationships, these tetraspanins were classified into 8 subfamilies (TspA to TspH). Six ancestral introns were identified within lepidopteran tetraspanin genes. Tetraspanins in TspA, TspB, TspC, and TspD subfamilies exhibit highly similar gene organization, while tetraspanins in the remaining 4 subfamilies exhibited variation in intron loss and/or gain during evolution. Analysis of chromosomal distribution revealed a lepidopteran-specific cluster of 10 to 11 tetraspanins, likely formed by tandem duplication events. Selective pressure analysis indicated negative selection across all orthologous groups, with ω values ranging between 0.004 and 0.362. However, positive selection was identified at 18 sites within TspB5, TspC5, TspE3, and TspF10. Furthermore, spatiotemporal expression analysis of H. armigera tetraspanins demonstrated variable expression levels across different developmental stages and tissues, suggesting diverse functions of tetraspanin members in this globally important insect pest. Our findings establish a solid foundation for subsequent functional investigations of tetraspanins in lepidopteran species.

The larva and adult of Helicoverpa armigera use differential gustatory receptors to sense sucrose

Almost all herbivorous insects feed on plants and use sucrose as a feeding stimulant, but the molecular basis of their sucrose reception remains unclear. Helicoverpa armigera as a notorious crop pest worldwide mainly feeds on reproductive organs of many plant species in the larval stage, and its adult draws nectar. In this study, we determined that the sucrose sensory neurons located in the contact chemosensilla on larval maxillary galea were 100-1000 times more sensitive to sucrose than those on adult antennae, tarsi, and proboscis. Using the Xenopus expression system, we discovered that Gr10 highly expressed in the larval sensilla was specifically tuned to sucrose, while Gr6 highly expressed in the adult sensilla responded to fucose, sucrose and fructose. Moreover, using CRISPR/Cas9, we revealed that Gr10 was mainly used by larvae to detect lower sucrose, while Gr6 was primarily used by adults to detect higher sucrose and other saccharides, which results in differences in selectivity and sensitivity between larval and adult sugar sensory neurons. Our results demonstrate the sugar receptors in this moth are evolved to adapt toward the larval and adult foods with different types and amounts of sugar, and fill in a gap in sweet taste of animals.

Genomic identification and evolutionary analysis of chemosensory receptor gene families in two Phthorimaea pest species: insights into chemical ecology and host adaptation

BACKGROUND

Insects rely on sophisticated sensitive chemosensory systems to sense their complex chemical environment. This sensory process involves a combination of odorant receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs) in the chemosensory system. This study focused on the identification and characterization of these three types of chemosensory receptor genes in two closely related Phthorimaea pest species, Phthorimaea operculella (potato tuber moth) and Phthorimaea absoluta (tomato leaf miner).

RESULTS

Based on manual annotation of the genome, we identified a total of 349 chemoreceptor genes from the genome of P. operculella, including 93 OR, 206 GR and 50 IR genes, while for P. absoluta, we identified 72 OR, 122 GR and 46 IR genes. Through phylogenetic analysis, we observed minimal differences in the number and types of ORs and IRs between the potato tuber moth and tomato leaf miner. In addition, we found that compared with those of tomato leaf miners, the gustatory receptor branch of P. operculella has undergone a large expansion, which may be related to P. absoluta having a narrower host range than P. operculella. Through analysis of differentially expressed genes (DEGs) of male and female antennae, we uncovered 45 DEGs (including 32ORs, 9 GRs, and 4 IRs).

CONCLUSIONS

Our research provides a foundation for exploring the chemical ecology of these two pests and offers new insights into the dietary differentiation of lepidopteran insects, while simultaneously providing molecular targets for developing environmentally friendly pest control methods based on insect chemoreception.

Genome sequencing and assembly of Indian golden silkmoth, Antheraea assamensis Helfer (Saturniidae, Lepidoptera)

Muga silkworm (Antheraea assamensis), one of the economically important wild silkmoths, is unique among saturniid silkmoths. It is confined to the North-eastern part of India. Muga silk has the highest value among the other silks. Unlike other silkmoths, A. assamensis has a low chromosome number (n = 15), and ZZ/ZO sex chromosome system. Here, we report the first high-quality draft genome of A. assamensis, assembled by employing the Illumina and PacBio sequencing platforms. The assembled genome of A. assamensis is 501.18 Mb long, with 2697 scaffolds and an N50 of 683.23 Kb. The genome encompasses 18,385 protein-coding genes, 86.29% of which were functionally annotated. Phylogenetic analysis of A. assamensis revealed its divergence from other Antheraea species approximately 28.7 million years ago. Moreover, an investigation into detoxification-related gene families, CYP450, GST, and ABC-transporter, revealed a significant expansion in A. assamensis as compared to the Bombyx mori. This expansion is comparable to Spodoptera litura, suggesting adaptive responses linked to the polyphagous behavior observed in these insects. This study provides valuable insights into the molecular basis of evolutionary divergence and adaptations in muga silkmoth. The genome assembly reported in this study will significantly help in the functional genomics studies on A. assamensis and other Antheraea species along with comparative genomics analyses of Bombycoidea insects.

Cross-pollination in seed-blended refuge and selection for Vip3A resistance in a lepidopteran pest as detected by genomic monitoring

The evolution of pest resistance to management tools reduces productivity and results in economic losses in agricultural systems. To slow its emergence and spread, monitoring and prevention practices are implemented in resistance management programs. Recent work suggests that genomic approaches can identify signs of emerging resistance to aid in resistance management. Here, we empirically examined the sensitivity of genomic monitoring for resistance management in transgenic Bt crops, a globally important agricultural innovation. Whole genome resequencing of wild North American Helicoverpa zea collected from non-expressing refuge and plants expressing Cry1Ab confirmed that resistance-associated signatures of selection were detectable after a single generation of exposure. Upon demonstrating its sensitivity, we applied genomic monitoring to wild H. zea that survived Vip3A exposure resulting from cross-pollination of refuge plants in seed-blended plots. Refuge seed interplanted with transgenic seed exposed H. zea to sublethal doses of Vip3A protein in corn ears and was associated with allele frequency divergence across the genome. Some of the greatest allele frequency divergence occurred in genomic regions adjacent to a previously described candidate gene for Vip3A resistance. Our work highlights the power of genomic monitoring to sensitively detect heritable changes associated with field exposure to Bt toxins and suggests that seed-blended refuge will likely hasten the evolution of resistance to Vip3A in lepidopteran pests.

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.

Analysis of insecticide resistance and de novo transcriptome assembly of resistance associated genes in the European grapevine moth, Lobesia botrana (Lepidoptera: Tortricidae)

The European grapevine moth Lobesia botrana (Denis & Shiffermüller 1776) is an economically important pest of the vine-growing areas worldwide. Chemical insecticides have been used for its control; however, its resistance status is largely unknown in many regions. We monitored the susceptibility of several L. botrana populations from Greece and Turkey. In addition, based on RNAseq transcriptome analysis, we identified and phylogenetically classify the cytochrome P450 genes of L. botrana, as well as analysed target site sequences and looked for the presence of known resistance mutations. Resistance against chlorantraniliprole, alpha-cypermethrin, spinetoram, etofenprox, and acetamiprid was very low (below 2.5-fold in all cases, compared to a reference strain from Greece) in all populations from Greece that were included in the study. However, resistance against indoxacarb (4-30-fold), spinosad (5-59-fold), and deltamethrin (18-30 fold) was detected in the L. botrana populations from Turkey, compared to a reference population from Turkey. De novo transcriptome assembly and manual annotation, and subsequent PCR-based analysis of insecticide target sequences (i.e. voltage-gated sodium channel - VGSC: target of pyrethroids and oxadiazines; nicotinic acetylcholine receptor subunit a6 - nAChR_α6: target of spinosad; ryanodine receptor - RyR: target of diamides; glutamate-gated chloride channel - GluCl: target of avermectins and; acetylcholinesterase - AChE: target of organophosphates) showed the absence of known resistance mutations in all specimens from both countries. Finally, the L. botrana CYPome (116 genes) was manually analysed and phylogenetically characterised, to provide resources for future studies that will aim the analysis of metabolic resistance.

A female-specific odorant receptor mediates oviposition deterrence in the moth Helicoverpa armigera

Finding ideal oviposition sites is a task of vital importance for all female insects. To ensure optimal conditions for their progeny, females of herbivorous insects detect not only the odors of a relevant host plant but also chemicals released by eggs, named oviposition-deterring pheromones (ODPs). It is reported that such chemicals play critical roles in suppressing female oviposition behavior; however, the molecular mechanism underlying the detection of egg-derived ODPs remains elusive. Here, we have identified three specific fatty acid methyl esters from the surface of eggs of Helicoverpa armigera serving as ODPs-methyl oleate (C18:1ME), methyl palmitate (C16:0ME), and methyl stearate (C18:0ME). We demonstrated that these ODPs are detected by the receptor, HarmOR56, exclusively expressed in sensilla trichodea on female antennae. To assess the significance of this receptor, we disrupted HarmOR56 in H. armigera using CRISPR-Cas9 and found that mutant females did not respond to the ODPs, neither in behavioral nor in electrophysiological tests. We therefore conclude that HarmOR56 is indispensable for identifying the ODPs. This study explores, for the first time, how a female-specific odorant receptor detects chemicals from conspecific eggs. Our data elucidate the intriguing biological phenomenon of repulsion to conspecific eggs during oviposition and contribute new insight into a female-specific olfactory pathway linked to reproduction.

Population structure and interspecific hybridisation of two invasive blowflies (Diptera: Calliphoridae) following replicated incursions into New Zealand

Rates of biological invasion are increasing globally, with associated negative effects on native biodiversity and ecosystem services. Among other genetic processes, hybridisation can facilitate invasion by producing new combinations of genetic variation that increase adaptive potential and associated population fitness. Yet the role of hybridisation (and resulting gene flow) in biological invasion in invertebrate species is under-studied. Calliphora hilli and Calliphora stygia are blowflies proposed to have invaded New Zealand separately from Australia between 1779 and 1841, and are now widespread throughout the country. Here, we analysed genome-wide single nucleotide polymorphisms (SNPs), generating genotyping-by-sequencing data for 154 individuals collected from 24 populations across New Zealand and Australia to assess the extent of gene flow and hybridisation occurring within and between these blowflies and to better understand their overall population structure. We found that New Zealand populations of both species had weak genetic structure, suggesting high gene flow and an absence of dispersal limitations across the country. We also found evidence that interspecific hybridisation is occurring in the wild between C. hilli and C. stygia in both the native and invasive ranges, and that intraspecific admixture is occurring among populations at appreciable rates. Collectively, these findings provide new insights into the population structure of these two invasive invertebrates and highlight the potential importance of hybridisation and gene flow in biological invasion.

Discovery of Insect Attractants Based on the Functional Analyses of Female-Biased Odorant Receptors and Their Orthologs in Two Closely Related Species

Olfaction plays an instrumental role in host plant selection by phytophagous insects. Helicoverpa assulta and Helicoverpa armigera are two closely related moth species with different host plant ranges. In this study, we first comparatively analyzed the function of 11 female-biased odorant receptors (ORs) and their orthologs in the two species by the Drosophila T1 neuron expression system and then examined the electroantennography responses of the two species to the most effective OR ligands. Behavioral assays using a Y-tube olfactometer indicate that guaiene, the primary ligand of HassOR21-2 and HarmOR21-2, only attracts the females, while benzyl acetone, the main ligand of HassOR35 and HarmOR35, attracts both sexes of the two species. Oviposition preference experiments further confirm that guaiene and benzyl acetone are potent oviposition attractants for the mated females of both species. These findings deepen our understanding of the olfactory coding mechanisms of host plant selection in herbivorous insects and provide valuable attractants for managing pest populations.

The most polyphagous insect herbivore? Host plant associations of the Meadow spittlebug, Philaenus spumarius (L.)

A comprehensive list of all known host plant species utilised by the Meadow Spittlebug (Philaenus spumarius (L.)) is presented, compiled from published and unpublished sources. P. spumarius feeds on 1311 host plants in 631 genera and 117 families. This appears, by a large margin, to be the greatest number of host species exploited by any herbivorous insect. The Asteraceae (222 species) and Rosaceae (110) together account for 25% of all host species. The Fabaceae (76) and Poaceae (73), are nearly tied for third and fourth place and these four families, combined with the Lamiaceae (62), Apiaceae (50), Brassicaceae (43) and Caprifoliaceae (34), comprise about half of all host species. Hosts are concentrated among herbaceous dicots but range from ferns and grasses to shrubs and trees. Philaenus spumarius is an "extreme polyphage", which appears to have evolved from a monophage ancestor in the past 3.7 to 7.9 million years. It is also the primary European vector of the emerging plant pathogen Xylella fastidiosa. Its vast host range suggests that it has the potential to spread X. fastidiosa among multiple hosts in any environment in which both the spittlebug and bacterium are present. Fully 47.9% of all known hosts were recorded in the Xylella-inspired BRIGIT citizen science P. spumarius host survey, including 358 hosts new to the documentary record, 27.3% of the 1311 total. This is a strong demonstration of the power of organized amateur observers to contribute to scientific knowledge.

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.

Absence of known knockdown resistance mutations but fixation of CYP337B3 was detected in field populations of Helicoverpa armigera across China

The cotton bollworm (Helicoverpa armigera) is a worldwide agricultural pest that infests many important crops. Pyrethroids targeting the voltage-gated sodium channel (VGSC) have been long used in the control of the cotton bollworm. Two amino acid substitutions (D1561V and E1565G) in H. armigera VGSC (HaVGSC) and the presence of a chimeric P450 gene (CYP337B3) have been documented to be associated with pyrethroid resistance. To understand the current occurrence of kdr mutations and the CYP337B3 gene in Chinese H. armigera populations, high-throughput amplicon sequencing was adopted to detect potential nucleotide variations in three fragments of the VGSC gene that cover 10 reported knockdown resistance (kdr) sites in insects, and gene-specific PCR was performed to examine the presence of CYP337B3 gene in H. armigera samples collected across China. The nucleotide variation analysis revealed a wealth of nucleotide variations in not only exons but also introns in the VGSC gene in Chinese H. armigera populations. However, neither previously reported kdr-conferring amino acid replacements nor other non-synonymous mutations were observed in a total of 1439 examined individuals. Population genetic analysis suggested that the H. armigera population in Nanchang, Jiangxi Province (JNC) had a moderate genetic differentiation from other populations, while no significant divergence was observed in other populations in northern and northwestern China. The CYP337B3 was present in all the examined individuals, indicating that CYP337B3 is extensively fixed in H. armigera populations across China. These results support that point mutations in VGSC are not a major factor involved in the current pyrethroid resistance in H. armigera. Instead, CYP337B3 plays a prevalent role in the development of resistance to pyrethroids in H. armigera.

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.

Evolutionary shifts in pheromone receptors contribute to speciation in four Helicoverpa species

Male moths utilize their pheromone communication systems to distinguish potential mates from other sympatric species, which contributes to maintaining reproductive isolation and even drives speciation. The molecular mechanisms underlying the evolution of pheromone communication systems are usually studied between closely-related moth species for their similar but divergent traits associated with pheromone production, detection, and/or processing. In this study, we first identified the functional differentiation in two orthologous pheromone receptors, OR14b, and OR16, in four Helicoverpa species, Helicoverpa armigera, H. assulta, H. zea, and H. gelotopoeon. To understand the substrate response specificity of these two PRs, we performed all-atom molecular dynamics simulations of OR14b and OR16 based on AlphaFold2 structural prediction, and molecular docking, allowing us to predict a few key amino acids involved in substrate binding. These candidate residues were further tested and validated by site-directed mutagenesis and functional analysis. These results together identified two hydrophobic amino acids at positions 164 and 232 are the determinants of the response specificity of HarmOR14b and HzeaOR14b to Z9-14:Ald and Z9-16:Ald by directly interacting with the substrates. Interestingly, in OR16 orthologs, we found that position 66 alone determines the specific binding of Z11-16:OH, likely via allosteric interactions. Overall, we have developed an effective integrated method to identify the critical residues for substrate selectivity of ORs and elucidated the molecular mechanism of the diversification of pheromone recognition systems.

Global patterns of genomic and phenotypic variation in the invasive harlequin ladybird

BACKGROUND

The harlequin ladybird Harmonia axyridis (Coleoptera: Coccinellidae), native to Asia, has been introduced to other major continents where it has caused serious negative impacts on local biodiversity. Though notable advances to understand its invasion success have been made during the past decade, especially with then newer molecular tools, the conclusions reached remain to be confirmed with more advanced genomic analyses and especially using more samples from larger geographical regions across the native range. Furthermore, although H. axyridis is one of the best studied invasive insect species with respect to life history traits (often comparing invasive and native populations), the traits responsible for its colonization success in non-native areas warrant more research.

RESULTS

Our analyses of genome-wide nuclear population structure indicated that an eastern Chinese population could be the source of all non-native populations and revealed several putatively adaptive candidate genomic loci involved in body color variation, visual perception, and hemolymph synthesis. Our estimates of evolutionary history indicate (1) asymmetric migration with varying population sizes across its native and non-native range, (2) a recent admixture between eastern Chinese and American populations in Europe, (3) signatures of a large progressive, historical bottleneck in the common ancestors of both populations and smaller effective sizes of the non-native population, and (4) the southwest origin and subsequent dispersal routes within its native range in China. In addition, we found that while two mitochondrial haplotypes-Hap1 and Hap2 were dominant in the native range, Hap1 was the only dominant haplotype in the non-native range. Our laboratory observations in both China and USA found statistical yet slight differences between Hap1 and Hap2 in some of life history traits.

CONCLUSIONS

Our study on H. axyridis provides new insights into its invasion processes into other major continents from its native Asian range, reconstructs a geographic range evolution across its native region China, and tentatively suggests that its invasiveness may differ between mitochondrial haplotypes.

Mutation of the clock gene timeless disturbs diapause induction and adult emergence rhythm in Helicoverpa armigera

BACKGROUND

Circadian rhythms are physical and behavioral changes that follow the 24-h cycle of Earth's light and temperature and are regulated by clock genes. Timeless (Tim) has been identified as a canonical clock gene in some insects, however, its functions have been little studied in lepidopteran pests.

RESULTS

To investigate Tim (HaTim) gene function in Helicoverpa armigera, an important lepidopteran pest, we obtained the HaTim mutant using the CRISPR/Cas9 gene editing system. Our results showed that the transcript levels of HaTim rhythmically peaked at night in heads of the wild larvae and adult, and the diel expression of HaTim was sensitive to photoperiod and temperature. The expression rhythms of other clock genes, such as HaPer, HaCry1, HaCry2 and HaCwo, were disturbed in the HaTim mutant larvae, as that stage is a sensitivity period for diapause induction. Fifth-instar wild-type larvae could be induced to pupate in diapause under a short-day photoperiod and low temperature, however, fifth-instar HaTim mutant larvae could not be induced under the same conditions. In addition, the emergence of wild-type adults peaked early at night, but the rhythm was disturbed in the HaTim mutant with arrhythmic expression of some clock genes, such as HaPer, HaCry1 and HaCwo in adults.

CONCLUSION

Our results suggest that the clock gene Tim is involved in diapause induction and adult emergence in H. armigera, and is a potential target gene for controlling pest. © 2023 Society of Chemical Industry.

Population Genomics of Nonrecessive Resistance to Bt Toxin Cry1Ac in Helicoverpa armigera From Northern China

Transgenic crops that produce insecticidal proteins from Bacillus thuringiensis (Bt) have provided control of some key pests since 1996. However, the evolution of resistance by pests reduces the benefits of Bt crops. Resistance to Bt crops that is not recessively inherited is especially challenging to manage. Here we analyzed nonrecessive resistance to Bt toxin Cry1Ac in eight field populations of Helicoverpa armigera sampled in 2018 from northern China, where this global pest has been exposed to Cry1Ac in Bt cotton since 1997. Bioassays revealed 7.5% of field-derived larvae were resistant to Cry1Ac of which 87% had at least one allele conferring nonrecessive resistance. To analyze this nonrecessive resistance, we developed and applied a variant of a genomic mapping approach called quantitative trait locus (QTL)-seq. This analysis identified a region on chromosome 10 associated with nonrecessive resistance to Cry1Ac in all 21 backcross families derived from field-collected moths. Individual sequencing revealed that all 21 field-collected resistant grandparents of the backcross families had a previously identified dominant point mutation in the tetraspanin gene HaTSPAN1 that occurs in the region of chromosome 10 identified by QTL-seq. QTL-seq also revealed a region on chromosome 26 associated with nonrecessive resistance in at most 14% of the backcross families. Overall, the results imply the point mutation in HaTSPAN1 is the primary genetic basis of nonrecessive resistance to Cry1Ac in field populations of H. armigera from northern China. Moreover, because nonrecessive resistance is predominant, tracking the frequency of this point mutation could facilitate resistance monitoring in the region.

The plant specialized metabolite epicatechin- 3-gallate (EC3G) perturbs lipid metabolism and attenuates fat accumulation in pigeonpea pod borer, Helicoverpa armigera

Control of pod borer Helicoverpa armigera, a notorious polyphagous pest requires paramount attention with focus on environment-friendly management approaches. Overproduction of catechins (epigallocatechin-EGC and epicatechin-3-gallate-EC3G) in the pod borer-resistant pigeonpea wild relative, Cajanus platycarpus during continued herbivory prodded us to assess their underlying molecular effect on H. armigera. Significant reduction in larval and pupal growth parameters was observed when reared on artificial diet incorporated with 100 ppm EC3G vis a vis 100 ppm EGC and EGC + EC3G. Comparative RNAseq analyses of larvae that fed on normal and EC3G-incorporated diet revealed 62 differentially expressed genes dominated by detoxification and lipid metabolism. While lipase and fatty acid-binding protein 2-like were up-regulated, delta9-FADS-like involved in fatty acid synthesis was downregulated, indicating effect of EC3G on fat metabolism. Validation of RNAseq data by qPCR; midgut glutathione-S-transferase and esterase assays depicted increased lipolysis and reduced lipogenesis in EC3G-fed larvae. Additionally, differential accumulation of stearic acid and oleic acid in EC3G-fed and control larvae/adults ascertained perturbation in lipogenesis. Supported by modelling, molecular docking and simulations, we demonstrate the possible involvement of the insect adipokinetic hormone receptor (AKHR) in the EC3G-mediated response. The study demonstrates plant specialized metabolite EC3G as a potential candidate for H. armigera control.

Identification and characterization of ABC proteins in an important rice insect pest, Cnaphalocrocis medinalis unveil their response to Cry1C toxin

Rice leaffolder (Cnaphalocrocis medinalis) is an important insect pest in paddy fields. Due to their essential role in the physiology and insecticidal resistance, ATP-binding cassette (ABC) proteins were studied in many insects. In this study, we identified the ABC proteins in C. medinalis through genomic data and analyzed their molecular characteristics. A total of 37 sequences with nucleotide-binding domain (NBD) were identified as ABC proteins and belonged to eight families (ABCA-ABCH). Four structure styles of ABC proteins were found in C. medinalis, including full structure, half structure, single structure, and ABC2 structure. In addition to these structures, TMD-NBD-TMD, NBD-TMD-NBD, and NBD-TMD-NBD-NBD were found in C. medinalis ABC proteins. Docking studies suggested that in addition to the soluble ABC proteins, other ABC proteins including ABCC4, ABCH1, ABCG3, ABCB5, ABCG1, ABCC7, ABCB3, ABCA3, and ABCC5 binding with Cry1C had higher weighted scores. The upregulation of ABCB1 and downregulation of ABCB3, ABCC1, ABCC7, ABCG1, ABCG3, and ABCG6 were associated with the C. medinalis response to Cry1C toxin. Collectively, these results help elucidate the molecular characteristics of C. medinalis ABC proteins, pave the way for further functional studies of C. medinalis ABC proteins, including their interaction with Cry1C toxin, and provide potential insecticide targets.

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.

Removal of an Aminopeptidase N From Midgut Brush Border Does Not Affect Susceptibility of Spodoptera litura (Lepidoptera: Noctuidae) Larvae to Four Insecticidal Proteins of Bacillus thuringiensis (Bacillales: Bacillaceae)

Spodoptera litura is one of the most destructive lepidopteran insects of cabbages and cauliflowers in the world. Cry1 and Vip3 toxins from Bacillus thuringiensis have been reported to show toxicity in multiple lepidopteran insects. Binding of toxic molecules to specific receptors on the midgut epithelial cells is known to be a key step in the action mode of Bt toxins. Aminopeptidase N (APN) -like proteins have been reported to be binding sites of multiple Cry toxins in the midgut of Cry susceptible insects. In the present study, we identified six midgut APNs by analysis of the genome and midgut transcriptome of S. litura. CRISPR/Cas9 mediated gene-knockout system was utilized to mutate the GPI-anchor signal peptide at the C terminus of SlAPN1. SlAPN1 was verified to be removed from the midgut brush border membrane vesicles of a homozygous knockout strain of S. litura (SlAPN1-KO). Bioassay results indicated that susceptibility of the SlAPN1-KO strain to Cry1Aa, Cry1Ac, Cry1Ca, and Vip3Aa toxins was close to that of the wild-type strain of S. litura. RT-qPCR results showed that the transcriptional level of SlAPN2-6 was not up-regulated after knockout of the SlAPN1. Results in this study indicated that the SlAPN1 did not play a critical role in the pathway of toxicity of Cry1Aa, Cry1Ac, Cry1Ca, and Vip3Aa toxins in S. litura.

Down-regulation of HaABCC3, potentially mediated by a cis-regulatory mechanism, is involved in resistance to Cry1Ac in the cotton bollworm, Helicoverpa armigera

Evolution of resistance to Cry proteins in multiple pest insects has been threatening the sustainable use of Bacillus thuringiensis (Bt)-transgenic crops. Better understanding about the mechanism of resistance to Cry proteins in insects is needed. Our preliminary study reported that the transcription of HaABCC3 was significantly decreased in a near-isogenic line (LFC2) of a Cry1Ac-resistant strain (LF60) of the global pest Helicoverpa armigera. However, the causality between HaABCC3 downregulation and resistance to Cry1Ac remains to be verified, and the regulatory mechanism underlying the HaABCC3 downregulation is still unclear. In this study, our data showed that both HaABCC3 and HaABCC3 downregulation were genetically linked to resistance to Cry1Ac in LF60. However, no InDels were observed in the coding sequence of HaABCC3 from LF60. Furthermore, F₁ offspring from the cross of LF60 and a HaABCC2/3-knockout mutant exhibited moderate resistance to Cry1Ac toxin; this indicated that the high resistance to Cry1Ac toxin in LF60 may have resulted from multiple genetic factors, including HaABCC2 mis-splicing and HaABCC3 downregulation. Results from luciferase reporter assays showed that promoter activity of HaABCC3 in LF60 was significantly lower than that in the susceptible strain, which indicated that HaABCC3 downregulation was likely mediated by promoter variation. Consistently, multiple variations of the GATA- or FoxA-binding sites in the promoter region of HaABCC3 were identified. Collectively, all results in this study suggested that the downregulation of HaABCC3 observed in the H. armigera LF60 strain, which is resistant to Cry1Ac, may be mediated by a cis-regulatory mechanism.

Chromosome-level genome of black cutworm provides novel insights into polyphagy and seasonal migration in insects

BACKGROUND

The black cutworm, Agrotis ipsilon, is a serious global underground pest. Its distinct phenotypic traits, especially its polyphagy and ability to migrate long distances, contribute to its widening distribution and increasing difficulty of control. However, knowledge about these traits is still limited.

RESULTS

We generated a high-quality chromosome-level assembly of A. ipsilon using PacBio and Hi-C technology with a contig N50 length of ~ 6.7 Mb. Comparative genomic and transcriptomic analyses showed that detoxification-associated gene families were highly expanded and induced after insects fed on specific host plants. Knockout of genes that encoded two induced ABC transporters using CRISPR/Cas9 significantly reduced larval growth rate, consistent with their contribution to host adaptation. A comparative transcriptomic analysis between tethered-flight moths and migrating moths showed expression changes in the circadian rhythm gene AiCry2 involved in sensing photoperiod variations and may receipt magnetic fields accompanied by MagR and in genes that regulate the juvenile hormone pathway and energy metabolism, all involved in migration processes.

CONCLUSIONS

This study provides valuable genomic resources for elucidating the mechanisms involved in moth migration and developing innovative control strategies.

The oriental armyworm genome yields insights into the long-distance migration of noctuid moths

The oriental armyworm, Mythimna separata, is known for its long-distance seasonal migration and environment-dependent phase polymorphisms. Here, we present a chromosome-level genome reference and integrate multi-omics, functional genetics, and behavioral assays to explore the genetic bases of the hallmark traits of M. separata migration. Gene family comparisons show expansion of gustatory receptor genes in this cereal crop pest. Functional investigation of magnetoreception-related genes and associated flight behaviors suggest that M. separata may use the geomagnetic field to guide orientation in its nocturnal flight. Comparative transcriptome characterizes a suite of genes that may confer the observed plasticity between phases, including genes involved in protein processing, hormone regulation, and dopamine metabolism. We further report molecular signatures that underlie the dynamic regulation of a migratory syndrome coordinating reproduction and flight. Our study yields insights into environment-dependent developmental plasticity in moths and advances our understanding of long-distance migration in nocturnal insect pests.

A General Signal Pathway to Regulate Multiple Detoxification Genes Drives the Evolution of Helicoverpa armigera Adaptation to Xenobiotics

The study of insect adaptation to the defensive metabolites of host plants and various kinds of insecticides in order to acquire resistance is a hot topic in the pest-control field, but the mechanism is still unclear. In our study, we found that a general signal pathway exists in H. armigera which can regulate multiple P450s, GSTs and UGTs genes to help insects decrease their susceptibility to xenobiotics. Knockdown of HaNrf2 and HaAhR expression could significantly increase the toxicity of xenobiotics to H. armigera, and simultaneously decrease the gene expression of P450s, GSTs and UGTs which are related to the xenobiotic metabolism and synthesis of insect hormone pathways. Then, we used EMSA and dual luciferase assay to verify that a crosstalk exists between AhR and Nrf2 to regulate multiple P450s, GSTs and UGTs genes to mediate H. armigera susceptibility to plant allelochemicals and insecticides. The detoxification genes' expression network which can be regulated by Nrf2 and AhR is still unknown, and there were also no reports about the crosstalk between AhR and Nrf2 that exist in insects and can regulate multiple detoxification genes' expression. Our results provide a new general signaling pathway to reveal the adaptive mechanism of insects to xenobiotics and provides further insight into designing effective pest-management strategies to avoid the overuse of insecticides.

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.

Improved genome assembly provides new insights into the environmental adaptation of the American cockroach, Periplaneta americana

The synanthropic pest and a model organism for entomological research, American cockroach, Periplaneta americana (Linnaeus), can survive in unfavorable environments for humans. To investigate the genetic mechanisms of success in environmental adaptation of P. americana, we de novo reassembled its whole genome based on next-generation sequencing and PacBio sequencing. The final genome reassembly consisted of approximately 3.34 Gb with scaffold N50 of 465.51 Kb. The completeness (95.4%) of the complete genome was evaluated with single-copy orthologous genes using BUSCO. We identified 18,618 protein-coding genes, 16,443 (88.32%) of which were well supported by public protein databases. We identified 482.04 Mb (approximately 14.45%) repeat elements, 1,385,093 perfect microsatellites simple sequence repeats in P. americana genome, which was higher than other four Blattaria insects. Comparative genomics analysis revealed obvious expansion in the gene families associated with chemoreception (olfactory receptors, gustatory receptors, ionotropic glutamate receptors, chemosensory protein, and sensory neuron membrane protein), which provided the necessary information for functional characterization of the chemosensory receptors of P. americana, with potential for new or refined applications of semiochemicals-based control of this pest insect. Similarly, gene families (cytochrome P450s, carboxyl/choline esterases, and UDP-glycosyl-transferases) encoding receptors for bitter or toxic substances and detoxification enzymes were obviously expanded in P. americana, enabling its ability to detect and detoxify many toxins. Enrichment analysis of positively selected genes in P. americana revealed items associated with metabolic process and catalytic activity, which possibly contributed to the pesticide resistance of P. americana. We also analyzed the homologs to antimicrobial peptide genes reported in the Drosophila genome, and identified two attacins and seven defensins in P. americana. Our data and findings will substantially facilitate molecular studies in P. americana, including elucidation of detoxification mechanisms of xenobiotic, as well as development of new pest management strategies for the control of pests like P. americana.

Metabonomics reveals that entomopathogenic nematodes mediate tryptophan metabolites that kill host insects

The entomopathogenic nematode (EPN) Steinernema feltiae, which carries the symbiotic bacterium Xenorhabdus bovienii in its gut, is an important biocontrol agent. This EPN could produce a suite of complex metabolites and toxin proteins and lead to the death of host insects within 24–48 h. However, few studies have been performed on the key biomarkers released by EPNs to kill host insects. The objective of this study was to examine what substances produced by EPNs cause the death of host insects. We found that all densities of nematode suspensions exhibited insecticidal activities after hemocoelic injection into Galleria mellonella larvae. EPN infection 9 h later led to immunosuppression by activating insect esterase activity, but eventually, the host insect darkened and died. Before insect immunity was activated, we applied a high-resolution mass spectrometry-based metabolomics approach to determine the hemolymph of the wax moth G. mellonella infected by EPNs. The results indicated that the tryptophan (Trp) pathway of G. mellonella was significantly activated, and the contents of kynurenine (Kyn) and 3-hydroxyanthranilic acid (3-HAA) were markedly increased. Additionally, 3-HAA was highly toxic to G. mellonella and resulted in corrected mortalities of 62.50%. Tryptophan metabolites produced by EPNs are a potential marker to kill insects, opening up a novel line of inquiry into exploring the infestation mechanism of EPNs.

Comparative genomic analysis of ABC transporter genes in Tenebrio molitor and four other tenebrionid beetles (Coleoptera: Tenebrionidea)

ATP-binding cassette (ABC) transporters, one of the largest transmembrane protein families, transport a diverse number of substate across membranes. Details of their diverse physiological functions have not been established. Here, we identified 87 ABC transporter genes in the genomes of Tenebrio molitor along with those from Asbolus verrucosus (104), Hycleus cichorii (65), and Hycleus phaleratus (80). Combining these genes (336 in total) with genes reported in Tribolium castaneum (73), we analyzed the phylogeny of ABC transporter genes in all five Tenebrionids. They are assigned into eight subfamilies (ABCA-H). In comparison to other species, the ABCC subfamily in this group of beetles appears expanded. The expression profiles of the T. molitor genes at different life stages and in various tissues were also investigated using transcriptomic analysis. Most of them display developmental specific expression patterns, suggesting to us their possible roles in development. Most of them are highly expressed in detoxification-related tissues including gut and Malpighian tubule, from which we infer their roles in insecticide resistance. We detected specific or abundant expressions of many ABC transporter genes in various tissues such as salivary gland, ovary, testis, and antenna. This new information helps generate new hypotheses on their biological significance within tissues.

Functional analysis of a bitter gustatory receptor highly expressed in the larval maxillary galea of Helicoverpa armigera

Many plant secondary substances are feeding deterrents for insects and play a key role in the selection of host plants. The taste sensilla of phytophagous insects contain gustatory sensory neurons sensitive to deterrents but the molecular basis of deterrent chemoreception remains unknown. We investigated the function of Gr180, the most highly expressed bitter gustatory receptor in the maxillary galea of Helicoverpa armigera larvae. Functional analyses using the Xenopus oocyte expression system and two-electrode voltage clamp revealed that the oocytes expressing Gr180 responded to coumarin. Tip recording results showed that the medial sensilla styloconica of the maxilla of fifth instar larvae exhibited electrophysiological responses to coumarin. Two-choice feeding bioassays confirmed that coumarin inhibited larval feeding. A homozygous mutant strain of H. armigera with truncated Gr180 proteins (Gr180-/-) was established using the CRISPR-Cas9 system. The responses of the medial sensilla styloconica in Gr180-/- to coumarin were almost abolished, and the responses to sinigrin and strychnine were also significantly decreased. Knockout of Gr180 alleviated the feeding deterrent effects of coumarin, sinigrin, and strychnine. Thus, we conclude that Gr180 is a receptor responding to coumarin,and also participates in sensing sinigrin and strychnine. These results enhance our understanding of the gustatory sensing mechanisms of phytophagous insects to deterrents.

Mutagenesis of the odorant receptor co-receptor (Orco) reveals severe olfactory defects in the crop pest moth Helicoverpa armigera

BACKGROUND

Odorant receptors (ORs) as odorant-gated ion channels play a crucial role in insect olfaction. They are formed by a heteromultimeric complex of the odorant receptor co-receptor (Orco) and a ligand-selective Or. Other types of olfactory receptor proteins, such as ionotropic receptors (IRs) and some gustatory receptors (GRs), are also involved in the olfactory system of insects. Orco as an obligatory subunit of ORs is highly conserved, providing an opportunity to systematically evaluate OR-dependent olfactory responses.

RESULTS

Herein, we successfully established a homozygous mutant (Orco^-/-) of Helicoverpa armigera, a notorious crop pest, using the CRISPR/Cas9 gene editing technique. We then compared the olfactory response characteristics of wild type (WT) and Orco^-/- adults and larvae. Orco^-/- males were infertile, while Orco^-/- females were fertile. The lifespan of Orco^-/- females was longer than that of WT females. The expressions of most Ors, Irs, and other olfaction-related genes in adult antennae of Orco^-/- moths were not obviously affected, but some of them were up- or down-regulated. In addition, there was no change in the neuroanatomical phenotype of Orco^-/- moths at the level of the antennal lobe (including the macroglomerular complex region of the male). Using EAG and SSR techniques, we discovered that electrophysiological responses of Orco^-/- moths to sex pheromone components and many host plant odorants were absent. The upwind flight behaviors toward sex pheromones of Orco^-/- males were severely reduced in a wind tunnel experiment. The oviposition selectivity of Orco^-/- females to the host plant (green pepper) has completely disappeared, and the chemotaxis toward green pepper was also lost in Orco^-/- larvae.

CONCLUSIONS

Our study indicates that OR-mediated olfaction is essential for pheromone communication, oviposition selection, and larval chemotaxis of H. armigera, suggesting a strategy in which mate searching and host-seeking behaviors of moth pests could be disrupted by inhibiting or silencing Orco expression.

Discovery of ZQ-, a Novel Starting Point To Develop Inhibitors against the Potent Molecular Target Chitinase

In a previous study, we found that compound ZQ-8 had a strong impact on the growth and development of Helicoverpa armigera. However, the mechanism and target of ZQ-8 are not clear, which makes it difficult to optimize the structure of this compound. In this study, the preliminary mode of action of ZQ-8 was studied through RNA sequencing and molecular docking. We also analyzed the underlying mechanisms from the aspect of ZQ-8 with respect to chitinase. The results showed that ZQ-8 mainly affects chitinase activity in the epidermis of H. armigera. ZQ-8 can competitively combine with chitinase 2 and endochitinase to form a relatively stable complex of ZQ-8, resulting in the failure of chitinase to degrade chitin. These findings indicate that the epidermis of H. armigera was identified as the action site of ZQ-8 and chitinase 2 and endochitinase were potential targets.

Identification of Candidate Chemosensory Gene Families by Head Transcriptomes Analysis in the Mexican Fruit Fly, Anastrepha ludens Loew (Diptera: Tephritidae)

Insect chemosensory systems, such as smell and taste, are mediated by chemosensory receptor and non-receptor protein families. In the last decade, many studies have focused on discovering these families in Tephritidae species of agricultural importance. However, to date, there is no information on the Mexican fruit fly Anastrepha ludens Loew, a priority pest of quarantine importance in Mexico and other countries. This work represents the first effort to identify, classify and characterize the six chemosensory gene families by analyzing two head transcriptomes of sexually immature and mature adults of A. ludens from laboratory-reared and wild populations, respectively. We identified 120 chemosensory genes encoding 31 Odorant-Binding Proteins (OBPs), 5 Chemosensory Proteins (CSPs), 2 Sensory Neuron Membrane Proteins (SNMPs), 42 Odorant Receptors (ORs), 17 Ionotropic Receptors (IRs), and 23 Gustatory Receptors (GRs). The 120 described chemosensory proteins of the Mexican fruit fly significantly contribute to the genetic databases of insects, particularly dipterans. Except for some OBPs, this work reports for the first time the repertoire of olfactory proteins for one species of the genus Anastrepha, which provides a further basis for studying the olfactory system in the family Tephritidae, one of the most important for its economic and social impact worldwide.

The Genetic Basis of Gene Expression Divergence in Antennae of Two Closely Related Moth Species, Helicoverpa armigera and Helicoverpa assulta

The closely related species Helicoverpa armigera (H. armigera) and Helicoverpa assulta (H. assulta) have different host plant ranges and share two principal components of sex pheromones but with reversed ratios. The antennae are the main olfactory organ of insects and play a crucial role in host plant selection and mate seeking. However, the genetic basis for gene expression divergence in the antennae of the two species is unclear. We performed an allele-specific expression (ASE) analysis in the antennal transcriptomes of the two species and their F₁ hybrids, examining the connection between gene expression divergence and phenotypic differences. The results show that the proportion of genes classified as all cis was higher than that of all trans in males and reversed in females. The contribution of regulatory patterns to gene expression divergence in males was less than that in females, which explained the functional differentiation of male and female antennae. Among the five groups of F₁ hybrids, the fertile males from the cross of H. armigera female and H. assulta male had the lowest proportion of misexpressed genes, and the inferred regulatory patterns were more accurate. By using this group of F₁ hybrids, we discovered that cis-related regulations play a crucial role in gene expression divergence of sex pheromone perception-related proteins. These results are helpful for understanding how specific changes in the gene expression of olfactory-related genes can contribute to rapid evolutionary changes in important olfactory traits in closely related moths.

Family dinner: Transcriptional plasticity of five Noctuidae (Lepidoptera) feeding on three host plant species

Polyphagous insects often show specialization in feeding on different host plants in terms of survival and growth and, therefore, can be considered minor or major pests of particular hosts. Whether polyphagous insects employ a common transcriptional response to cope with defenses from diverse host plants is under-studied. We focused on patterns of transcriptional plasticity in polyphagous moths (Noctuidae), of which many species are notorious pests, in relation to herbivore performance on different host plants. We compared the transcriptional plasticity of five polyphagous moth species feeding and developing on three different host plant species. Using a comparative phylogenetic framework, we evaluated if successful herbivory, as measured by larval performance, is determined by a shared or lineage-specific transcriptional response. The upregulated transcriptional activity, or gene expression pattern, of larvae feeding on the different host plants and artificial control diet was highly plastic and moth species-specific. Specialization, defined as high herbivore success for specific host plants, was not generally linked to a lower number of induced genes. Moths that were more distantly related and showing high herbivore success for certain host plants showed shared expression of multiple homologous genes, indicating convergence. We further observed specific transcriptional responses within phylogenetic lineages. These expression patterns for specific host plant species are likely caused by shared evolutionary histories, for example, symplesiomorphic patterns, and could therefore not be associated with herbivore success alone. Multiple gene families, with roles in plant digestion and detoxification, were widely expressed in response to host plant feeding but again showed highly moth species-specific. Consequently, high herbivore success for specific host plants is also driven by species-specific transcriptional plasticity. Thus, potential pest moths display a complex and species-specific transcriptional plasticity.

Genome-Wide Analysis of Odorant and Gustatory Receptors in Six Papilio Butterflies (Lepidoptera: Papilionidae)

The chemical interactions of insects and host plants are shaping the evolution of chemosensory receptor gene families. However, the correlation between host range and chemoreceptor gene repertoire sizes is still elusive in Papilionidae. Here, we addressed the issue of whether host plant diversities are correlated with the expansions of odorant (ORs) or gustatory (GRs) receptors in six Papilio butterflies. By combining genomics, transcriptomics and bioinformatics approaches, 381 ORs and 328 GRs were annotated in the genomes of a generalist P. glaucus and five specialists, P. xuthus, P. polytes, P. memnon, P. machaon and P. dardanus. Orthologous ORs or GRs in Papilio had highly conserved gene structure. Five Papilio specialists exhibited a similar frequency of intron lengths for ORs or GRs, but which was different from those in the generalist. Phylogenetic analysis revealed 60 orthologous OR groups, 45 of which shared one-to-one relationships. Such a single gene in each butterfly also occurred in 26 GR groups. Intriguingly, bitter GRs had fewer introns than other GRs and clustered into a large clade. Focusing on the two chemoreceptor gene families in P. xuthus, most PxutORs (52/58) were expressed in antennae and 31 genes in reproductive tissues. Eleven out of 28 foretarsus-expressed PxutGRs were female-biased genes, as strong candidates for sensing oviposition stimulants. These results indicate that the host range may not shape the large-scale expansions of ORs and GRs in Papilio butterflies and identify important molecular targets involved in olfaction, oviposition or reproduction in P. xuthus.

Genomic Variation in the Tea Leafhopper Reveals the Basis of Adaptive Evolution

Tea green leafhopper (TGL), Empoasca onukii, is of biological and economic interest. Despite numerous studies, the mechanisms underlying its adaptation and evolution remain enigmatic. Here, we use previously untapped genome and population genetics approaches to examine how the pest adapted to different environmental variables and thus has expanded geographically. We complete a chromosome-level assembly and annotation of the E. onukii genome, showing notable expansions of gene families associated with adaptation to chemoreception and detoxification. Genomic signals indicating balancing selection highlight metabolic pathways involved in adaptation to a wide range of tea varieties grown across ecologically diverse regions. Patterns of genetic variations among 54 E. onukii samples unveil the population structure and evolutionary history across different tea-growing regions in China. Our results demonstrate that the genomic changes in key pathways, including those linked to metabolism, circadian rhythms, and immune system functions, may underlie the successful spread and adaptation of E. onukii. This work highlights the genetic and molecular bases underlying the evolutionary success of a species with broad economic impacts, and provides insights into insect adaptation to host plants, which will ultimately facilitate more sustainable pest management.