A reference gene set for sex pheromone biosynthesis and degradation genes from the diamondback moth, Plutella xylostella, based on genome and transcriptome digital gene expression analyses

Duplications and Losses of the Detoxification Enzyme Glycosyltransferase 1 Are Related to Insect Adaptations to Plant Feeding

Insects have developed sophisticated detoxification systems to protect them from plant secondary metabolites while feeding on plants to obtain necessary nutrients. As an important enzyme in the system, glycosyltransferase 1 (GT1) conjugates toxic compounds to mitigate their harm to insects. However, the evolutionary link between GT1s and insect plant feeding remains elusive. In this study, we explored the evolution of GT1s across different insect orders and feeding niches using publicly available insect genomes. GT1 is widely present in insect species; however, its gene number differs among insect orders. Notably, plant-sap-feeding species have the highest GT1 gene numbers, whereas blood-feeding species display the lowest. GT1s appear to be associated with insect adaptations to different plant substrates in different orders, while the shift to non-plant feeding is related to several losses of GT1s. Most large gene numbers are likely the consequence of tandem duplications showing variations in collinearity among insect orders. These results reveal the potential relationships between the evolution of GT1s and insect adaptation to plant feeding, facilitating our understanding of the molecular mechanisms underlying insect-plant interactions.

Predation evaluation of the green lacewing, Chrysopa pallens on the pink tea mite pest, Acaphylla theae (Watt) (Acarina: Eriophyidae)

A better understanding of predator-prey interactions is crucial for the development of biological control strategies. The green lacewing, Chrysopa pallens, is a well-known generalist predator and reportedly functions as one of the most important biological control agents of insect pests. However, information regarding C. pallens' predation on tea plant pests, particularly notorious tea mites, remains largely unknown. In this study, we focused on the predator-prey relationship between C. pallens and an important tea mite pest, Acaphylla theae. We designed species-specific primers for the detection of A. theae DNA and established a PCR-based DNA gut content analysis assay. These results demonstrated that the primers were A. theae-specific and suitable for its molecular identification. The laboratory feeding experiment showed that the detectability success (DS50) of A. theae DNA remaining in C. pallens' guts was 2.9 h. We then performed a molecular detection of field predation, and achieved a 23.53% positive detection rate of A. theae DNA in the guts of field-collected C. pallens. This, for the first time, provides direct evidence that C. pallens can prey on A. theae in tea plantations. Finally, we tested the prey preference and estimated the predation ability of C. pallens on different developmental stages of A. theae. The results revealed that C. pallens had no significant preference for different developmental stages of A. theae. The functional responses of C. pallens' predation on different densities of A. theae at different developmental stages followed a Type II Holling model. The initial attack rate (a') ranged from 0.735 to 0.858 and the handling time (Th) was approximately 0.01. This study is the first to demonstrate the trophic interactions between C. pallens and A. theae and provides evidence for the development of biological control strategies against A. theae using C. pallens as a candidate predator.

Knockout of tyramine receptor 1 results in a decrease of oviposition, mating, and sex pheromone biosynthesis in female Plutella xylostella

BACKGROUND

Mating and oviposition are essential and closely coordinated events in the reproduction of moths. Although tyramine, a biogenic amine, can affect insect reproduction by binding its receptors, the specific regulatory mechanism has not yet been fully elucidated.

RESULTS

Plutella xylostella mutant with tyramine receptor 1 (TAR1) knockout (homozygous mutant with 7-bp deletion, Mut7) was developed by the CRISPR/Cas9 system to investigate the effect of TAR1 knockout on the reproduction of the moth. Compared with wild-type (WT), the egg yield of Mut7 female (Mut7F ) was significantly lower, no significant difference was observed in the egg size and hatching ratio between the groups. Further analysis showed that TAR1 knockout adversely affected ovary development, characterized by shorter ovarioles and fewer mature oocyte. Additionally, TAR1 knockout significantly reduced the occurrence of mating, resulting in a decrease in egg yield in Mut7F . The amounts of sex pheromones were quantified using gas chromatography-mass spectrometry. Results showed that the amounts of sex pheromone released by Mut7F were significantly lower before mating. Correspondingly, the messenger RNA (mRNA) levels of sex pheromone biosynthesis enzymes, including acetyl-CoA carboxylase (ACC) and desaturase (DES), were significantly lower in the Mut7F pheromone gland. The decreased sex pheromone biosynthesis in Mut7F , especially before re-mating, may be related to the underexpression of pheromone biosynthesis-activated neuropeptide (PBAN).

CONCLUSION

Overall, this study investigated the effect of PxTAR1 on oviposition and mating of P. xylostella. We report for the first time that TAR1 knockout could reduce the sex pheromone biosynthesis. These findings provide insights for developing a novel integrated pest control strategy based on mating interference. © 2023 Society of Chemical Industry.

Functional Analysis of Pheromone Biosynthesis Activating Neuropeptide Receptor Isoforms in Maruca vitrata

Insect sex pheromones are volatile chemicals that induce mating behavior between conspecific individuals. In moths, sex pheromone biosynthesis is initiated when pheromone biosynthesis-activating neuropeptide (PBAN) synthesized in the suboesophageal ganglion binds to its receptor on the epithelial cell membrane of the pheromone gland. To investigate the function of PBAN receptor (PBANR), we identified two PBANR isoforms, MviPBANR-B and MviPBANR-C, in the pheromone glands of Maruca vitrata. These two genes belong to G protein-coupled receptors (GPCRs) and have differences in the C-terminus but share a 7-transmembrane region and GPCR family 1 signature. These isoforms were expressed in all developmental stages and adult tissues. MviPBANR-C had the highest expression level in pheromone glands among the examined tissues. Through in vitro heterologous expression in HeLa cell lines, only MviPBANR-C-transfected cells responded to MviPBAN (≥5 µM MviPBAN), inducing Ca²⁺ influx. Sex pheromone production and mating behavior were investigated using gas chromatography and a bioassay after MviPBANR-C suppression by RNA interference, which resulted in the major sex pheromone component, E10E12-16:Ald, being quantitatively reduced compared to the control, thereby decreasing the mating rate. Our findings indicate that MviPBANR-C is involved in the signal transduction of sex pheromone biosynthesis in M. vitrata and that the C-terminal tail plays an important role in its function.

Two Antenna-Enriched Carboxylesterases Mediate Olfactory Responses and Degradation of Ester Volatiles in the German Cockroach Blattella germanica

Insects have evolved an extremely sensitive olfactory system that is essential for a series of physiological and behavioral activities. Some carboxylesterases (CCEs) comprise a major subfamily of odorant-degrading enzymes (ODEs) playing a crucial role in odorant signal inactivation to maintain the odorant receptor sensitivity. In this study, 93 CCEs were annotated in the genome of the German cockroach Blattella germanica, a serious urban pest. Phylogenetic and digital tissue expression pattern analyses identified two antenna-enriched CCEs, BgerCCE021e3 and BgerCCE021d1, as candidate ODEs. RNA interference (RNAi)-mediated knockdown of BgerCCE021e3 and BgerCCE021d1 resulted in partial anosmia with experimental insects exhibiting reduced attraction to ester volatile resources and slower olfactory responses than controls. Furthermore, enzymatic conversion of geranyl acetate by crude male antennal extracts from BgerCCE021e3 and BgerCCE021d1 RNAi insects was also significantly reduced. Our results provide evidence for CCE function in German cockroach olfaction and provide a basis for further exploring behavioral inhibitors that target olfactory-related CCEs.

Dynamic Roles of Insect Carboxyl/Cholinesterases in Chemical Adaptation

Insects have evolved several intricate defense mechanisms to adapt to their chemical environment. Due to their versatile capabilities in hydrolytic biotransformation, insect carboxyl/cholinesterases (CCEs) play vital roles in the development of pesticide resistance, facilitating the adaptation of insects to their host plants, and manipulating insect behaviors through the olfaction system. CCEs confer insecticide resistance through the mechanisms of qualitative or quantitative changes of CCE-mediated enhanced metabolism or target-site insensitivity, and may contribute to the host plant adaptation. CCEs represent the first odorant-degrading enzymes (ODEs) discovered to degrade insect pheromones and plant odors and remain the most promising ODE candidates. Here, we summarize insect CCE classification, currently characterized insect CCE protein structure characteristics, and the dynamic roles of insect CCEs in chemical adaptation.

Transcriptomic Identification and Expression Profile Analysis of Odorant-Degrading Enzymes from the Asian Corn Borer Moth, Ostrinia furnacalis

The Asian corn borer moth Ostrinia furnacalis is an important lepidopteran pest of maize in Asia. Odorant-degrading enzymes (ODEs), including carboxylesterases (CCEs), glutathione S-transferases (GSTs), cytochrome P450s (CYPs), UDP-glycosyltransferases (UGTs), and aldehyde oxidases (AOXs), are responsible for rapid inactivation of odorant signals in the insect antennae. In this study, we performed a transcriptome assembly for the antennae of O. furnacalis to identify putative ODE genes. Transcriptome sequencing revealed 35,056 unigenes, and 21,012 (59.94%) of these were annotated by searching against the reference sequences in the NCBI non-redundant (NR) protein database. For functional classification, these unigenes were subjected to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. We identified 79 genes encoding putative ODEs: 19 CCEs, 17 GSTs, 24 CYPs, 13 UGTs, and 6 AOXs. BLASTX best hit results indicated that these genes shared quite high amino acid identities with their respective orthologs from other lepidopteran species. Reverse transcription-quantitative PCR showed that OfurCCE2, OfurCCE5, and OfurCCE18 were enriched in male antennae, while OfurCCE7 and OfurCCE10 were enriched in female antennae. OfurCCE14 and OfurCCE15 were expressed at near-equal amounts in the antennae of both sexes. Our findings establish a solid foundation for future studies aimed at understanding the olfactory functions of these genes in O. furnacalis.

Aldehyde oxidases mediate plant toxicant susceptibility and fecundity in the red flour beetle, Tribolium castaneum

Aldehyde oxidases (AOXs) are a group of metabolic enzymes that play critical roles in the degradation of xenobiotics and chemicals. However, the physiological function of this enzyme in insects remains poorly understood. In this study, three TcAOX genes (TcAOX1, TcAOX2, TcAOX3) were identified and characterized from Tribolium castaneum genome. Spatiotemporal expression profiling showed that TcAOX1 expression was most highly expressed at the early pupal stage and was predominantly expressed in the antennae of adults, indicating that TcAOX1 was involved in the degradation of chemical signals; TcAOX2 expression was most highly expressed at the late pupal stage and was mainly expressed in the fat body, epidermis of larvae and adults, respectively; and TcAOX3 expression was in all stages and was primarily expressed in the head of adults. Moreover, the transcripts of TcAOX2 and TcAOX3 were significantly induced after exposure to plant oil, and RNA interference (RNAi) targeting of each of them enhanced the susceptibility of beetles to this plant toxicant, suggesting that these two genes are associated with plant toxicant detoxification. Intriguingly, knockdown of the TcAOX1 led to reductions in female egg-laying but unchanged the hatchability and the development of genital organs, suggesting that this gene may mediate fecundity by effecting the inactivation of chemical signals in T. castaneum. Overall, these results shed new light on the function of AOX genes in insects, and could facilitate the development of research on pest control management.

Transcriptome analysis and identification of sex pheromone biosynthesis and transport related genes in Atrijuglans hetaohei (Lepidoptera: Gelechioidea)

Atrijuglans hetaohei Yang (Lepidoptera: Gelechioidea) is one of the major pests that can seriously damage the walnut tree, leading to harvest loss. Sex pheromones regulate mating communication and reproduction in insects and provide targets for developing a novel pest control strategy. In this study, by transcriptomic sequencing and analysis of the female pheromone gland (PG) and male genitalia of A. hetaohei, we identified 92 putative genes, of which 7 desaturases (Dess), 8 fatty acyl reductases (FARs), 4 fatty acid synthetases (FASs), 2 aldehyde oxidases (AOXs), 4 acetyltransferases (ACTs), 1 chemosensory protein (CSP), and 2 odorant-binding proteins (OBPs) were predominantly expressed in the female PG, while 5 Dess, 11 FARs, 7 FASs, 6 AOXs, 1 ACT, and 1 CSP showed more robust expression in the male genitalia. Moreover, phylogenetic analysis revealed that 7 Dess and 1 FAR were grouped with genes involved in pheromone synthesis in other Lepidoptera species. Thus, we proposed that these candidate genes are possibly involved in the sex pheromone biosynthetic pathway in A. hetaohei. Our findings will provide a solid genetic basis for further exploring the function of the tissue-biased genes and may be useful to screen potential targets for interfering chemical communication in A. hetaohei.

Systematic analysis of olfactory protein-protein interactions network of fruitfly, Drosophila melanogaster

Olfaction is one of the physiological traits of insect behavior. Insects have evolved a sophisticated olfactory system and use a combined coding strategy to process general odor. Drosophila melanogaster is a powerful model to reveal the molecular and cellular mechanisms of odor detection. Identifying new olfactory targets through complex interactions will contribute to a better understanding of the functions, interactions, and signaling pathways of olfactory proteins. However, the mechanism of D. melanogaster olfaction is still unclear, and more olfactory proteins are required to be discovered. In this study, we tried to explore essential proteins in the olfactory system of D. melanogaster and conduct protein-protein interactions (PPIs) analysis. We constructed the PPIs network of the olfactory system of D. melanogaster, consisting of 863 proteins and 18,959 interactions. Various methods were used to perform functional enrichment analysis, topological analysis and cluster analysis. Our results confirmed that Class B scavenger receptors (SR-Bs), glutathione S-transferases (GSTs), and UDP-glycosyltransferases (UGTs) play an essential role in olfaction of D. melanogaster. The proteins obtained in this study can be used for subsequent functional identification in D. melanogaster olfactory.

Identification and Characterization of UDP-Glycosyltransferase Genes in a Cerambycid Beetle, Pharsalia antennata Gahan, 1894 (Coleoptera: Cerambycidae)

The cerambycid beetle, Pharsalia antennata Gahan, 1894 (Coleoptera: Cerambycidae), is a wood-boring pest that spends most of its life cycle in the trunks or under the bark of trees. These distinctive biological characteristics make it likely that this beetle will encounter a number of plant defensive compounds, coupled with a broad range of host plants, possibly resulting in the overexpression or expansion of uridine diphosphate (UDP)-glycosyltransferase (UGT) genes. Here, we identified and characterized the UGT gene family in P. antennata through transcriptome data, sequence and phylogenetic analyses, and PCR and homology modeling approaches. In total, 59 transcripts encoding UGTs were identified, 34 of which harbored full-length sequences and shared high conservation with the UGTs of Anoplophora glabripennis. Of the 34 PantUGTs, only 31.78% amino acid identity was observed on average, but catalytic and sugar binding residues were highly conserved. Phylogenetic analyses revealed four Cerambycidae-specific clades, including 30 members from P. antennata. Combining the transcriptome and PCR data showed that PantUGTs had a wide tissue expression, and the majority of the genes were presented mainly in antennae or abdomens, suggesting their putative roles in olfaction and detoxification. This study provides, for the first time, information on the molecular and genetic basis of P. antennata, greatly enhancing our knowledge of the detoxification-related UGT gene family.

Genome-Wide Identification of Aldehyde Oxidase Genes in Moths and Butterflies Suggests New Insights Into Their Function as Odorant-Degrading Enzymes

Aldehyde oxidases (AOXs) are common detoxifying enzymes in several organisms. In insects, AOXs act in xenobiotic metabolism and as odorant-degrading enzymes (ODEs). These last appear as crucial enzymes in the life cycle of insects, helping to reset their olfactory system, particularly in lepidopterans, which fulfill important ecological roles (e.g., pollination or destructive life cycles). A comprehensive understanding of their olfactory system has provided opportunities to study key chemosensory proteins. However, no significant advance has been made around lepidopteran AOXs research, and even less around butterflies, a recently evolved lineage. In this study we have identified novel AOX gene families in moths and butterflies in order to understand their role as ODEs. Eighteen genomes from both moths and butterflies were used for phylogenetics, molecular evolution and sequence analyses. We identified 164 AOXs, from which 91 are new. Their phylogeny showed two main clades that are potentially related to odorant-degrading function, where both moths and butterflies have AOXs. A first ODE-related clade seems to have a non-ditrysian origin, likely related to plant volatiles. A second ODE-related clade could be more pheromone-biased. Molecular evolution analysis suggests a slight purifying selection process, though a number of sites appeared under positive selection. ODE-related AOXs have changed a phenylalanine residue by proline in the active site. Finally, this study could serve as a reference for further evolutionary and functional studies around Lepidopteran AOXs.

Exploring the Terminal Pathway of Sex Pheromone Biosynthesis and Metabolism in the Silkworm

Simple Summary

Insect sex pheromone biosynthesis has received widespread attention, while the terminal pathway related to aldehyde synthesis and metabolism is still poorly understood at a molecular level. Previous studies found that the silkworm, Bombyx mori (Lepidoptera, Bombycidae), has two pheromone compounds, bombykol and bombykal, with a ratio of 11:1, while its closest wild relative, B. mandarina, only uses bombykol as a pheromone. In this study, sex pheromone gland transcriptomes were compared between the domestic and wild silkworms. All the candidate gene families were identified. Then we used the differentially expressed information, tissue and developmental expression profiles, and phylogenetic analysis to identify the putative causal genes involved in the terminal pathway. Our findings provide insights into the aldehyde synthesis and metabolism pathways and evolutionary conservation in moths.

Abstract

Sex pheromones are vital to sexual communication and reproduction in insects. Although some key enzymes in pheromone production have been well studied, information on genes involved in the terminal pathway is limited. The domestic silkworm employs a pheromone blend containing (E,Z)-10,12-hexadecadienol (bombykol) and analogous (E,Z)-10,12-hexadecadienal (bombykal); whereas, its wild ancestor B. mandarina uses only bombykol. The two closely related moths might be a good model for exploring the genes involved in aldehyde pheromone synthesis and metabolism. By deep sequencing and analyzing the sex pheromone gland (PG) transcriptomes; we identified 116 candidate genes that may be related to pheromone biosynthesis, metabolism, and chemoreception. Spatiotemporal expression profiles and differentially expressed analysis revealed that four alcohol oxidases (BmorAO1; 2; 3; and 4); one aldehyde reductase (BmorAR1); and one aldehyde oxidase (BmorAOX5) might be involved in the terminal pathway. Phylogenetic analysis showed that, except for BmorAO3 and MsexAO3, AOs did not show a conversed orthologous relationship among moths; whereas, ARs and AOXs were phylogenetically conserved. This study provides crucial candidates for further functional elucidation, and which may be utilized as potential targets to disrupt sexual communication in other moth pests.

Functional analysis of the odorant receptor coreceptor in odor detection in Grapholita molesta (lepidoptera: Tortricidae)

The olfactory system must detect and discriminate various semiochemicals in the environment. In response to such diversity, insects have evolved a family of odorant-gated ion channels composed of a common receptor (coreceptor, Orco) and a ligand-binding tuning odorant receptor (OR) that confers odour specificity. This study aims to examine the expression pattern of Orco gene of Grapholita molesta (GmolOrco) and to elucidate the role of GmolOrco in detecting G. molesta sex pheromone and green leaf volatiles by using gene silencing via RNA interference (RNAi) coupled antennal electrophysiological (EAG). Multiple sequence alignment showed that GmolOrco shared high sequence similarities with the Orco ortholog of lepidopterans. The results of real-time quantitative PCR detection demonstrated that GmolOrco was predominantly expressed in adult antennae and had the highest expression quantity in adult period among the different developmental stages. Compared with the noninjected controls, GmolOrco expression in GmolOrcodouble-stranded RNA (dsRNA)-injected males was reduced to 39.92% and that in females was reduced to 40.43%. EAG assays showed that the responses of GmolOrco-dsRNA injected males to sex pheromones (Z)-8-dodecenyl acetate (Z8-12:OAc) and (Z)-8-dodecenyl alcohol (Z8-12:OH) were significantly reduced, and the GmolOrco-dsRNA-injected female to green leaf volatile (Z)-3-hexenyl acetate also significantly declined. We inferred that Orco-mediated olfaction was different in male and female G. molesta adults and was mainly involved in recognizing the sex pheromones released by female moths.

Genes involved in the Type I pheromone biosynthesis pathway and chemoreception from the sex pheromone gland transcriptome of Dioryctria abietella

Dioryctria abietella is a coniferous seed orchard pest that can damage a series of host plants and cause huge losses to the forest economy. Sex pheromones play an important role in lepidopteran sex communication for reproduction and can be used as biological control agents to monitor and trap pests. However, the genes involved in the biosynthesis, transportation, and degradation of D. abietella sex pheromones have not been studied extensively. Transcriptome analysis of female D. abietella sex pheromone glands (PGs) revealed that 210 candidate genes might be involved in sex pheromone biosynthesis (139 genes) and chemoreception systems (71 genes). The gene expression patterns exhibited four desaturase genes (DabiDES4-7) and one fatty acid reductase gene (DabiFAR6), which were more highly expressed in sex pheromone glands than in other tissues, suggesting that these enzymes play an important role in D. abietella sex pheromone synthesis. In addition, most DabiOBPs showed high expression in antennae, but only DabiOBP4 exhibited specific expression in sex pheromone glands, suggesting that they may play many physiological roles in D. abietella. We put forth a reasonable hypothesis about type I pheromone biosynthesis pathways based on these genes identified in the D. abietella sex pheromone gland transcriptome. Our findings lay a foundation for population monitoring, mating disruption, mass trapping, and the development of ecologically acceptable management strategies.

Two carboxylesterase genes in Plutella xylostella associated with sex pheromones and plant volatiles degradation

BACKGROUND

Carboxyl/cholinesterases (CCEs) are thought to play a pivotal role in the degradation of sex pheromones and plant-derived odorants in insects, but their exact biochemistry and physiological functions remain unclear.

RESULTS

In this study, two paralogous antennae-enriched CCEs from Plutella xylostella (PxylCCE16a and 16c) were identified and functionally characterized. High-purity protein preparations of active recombinant PxylCCE16a and 16c have been obtained from Sf9 insect cells by Ni²⁺ affinity purification. Our results revealed that the purified recombinant PxylCCE016c is able to degrade two sex pheromone components Z9-14:Ac and Z11-16:Ac at 27.64 ± 0.79% and 24.40 ± 3.07%, respectively, while PxylCCE016a presented relatively lower activity. Additionally, a similar difference in activity was measured in plant-derived odorants. Furthermore, both CCEs displayed obvious preferences for the two sex pheromone components, especially on Z11-16:Ac (K_m values are in the range 7.82-45.06 μmol L^-1 ) which much lower than plant odorants (K_m values are in the range 1290-4030 μmol L^-1 ). Furthermore, the activity of the two newly identified CCEs is pH-dependent. The activity at pH 6.5 is obviously higher than that at pH 5.0. Interestingly, only PxylCCE016c can be inhibited by a common esterase inhibitor triphenyl phosphate (TPP) with LC₅₀ of 1570 ± 520 μmol L^-1 .

CONCLUSION

PxylCCE16c plays a more essential role in odorant degradation than PxylCCE16a. Moreover, the current study provides novel potential pesticide targets for the notorious moth Plutella xylostella. © 2021 Society of Chemical Industry.

Glycosyltransferases: the multifaceted enzymatic regulator in insects

Glycosyltransferases (GTs) catalyse the reaction of glyco-conjugation of various biomolecules by transferring the saccharide moieties from an activated nucleotide sugar to nucleophilic glycosyl acceptor. In insects, GTs show diverse temporal and site-specific expression patterns and thus play significant roles in forming the complex biomolecular structures that are necessary for insect survival, growth and development. Several insects exhibit GT-mediated detoxification as a key defence strategy against plant allelochemicals and xenobiotic compounds, as well as a mechanism for pesticide cross-resistance. Also, these enzymes act as crucial effectors and modulators in various developmental processes of insects such as eye development, UV shielding, cuticle formation, epithelial development and other specialized functions. Furthermore, many of the known insect GTs have been shown to play a fundamental role in other physiological processes like body pigmentation, cuticular tanning, chemosensation and stress response. This review provides a detailed overview of the multifaceted functionality of insect GTs and summarizes numerous case studies associated with it.

An Overview of Antennal Esterases in Lepidoptera

Lepidoptera are used as a model for the study of insect olfactory proteins. Among them, odorant degrading enzymes (ODEs), that degrade odorant molecules to maintain the sensitivity of antennae, have received less attention. In particular, antennal esterases (AEs; responsible for ester degradation) are crucial for intraspecific communication in Lepidoptera. Currently, transcriptomic and genomic studies have provided AEs in several species. However, efforts in gene annotation, classification, and functional assignment are still lacking. Therefore, we propose to combine evidence at evolutionary, structural, and functional level to update ODEs as well as key information into an easier classification, particularly of AEs. Finally, the kinetic parameters for putative inhibition of ODEs are discussed in terms of its role in future integrated pest management (IPM) strategies.

A candidate aldehyde oxidase in the antennae of the diamondback moth, Plutella xylostella (L.), is potentially involved in the degradation of pheromones, plant-derived volatiles and the detoxification of xenobiotics

Insect antennae play a fundamental role in perceiving and recognizing a broad spectrum of conventional semiochemicals and host plant-derived odors. As such, genes that are tightly associated with the antennae are thought to have olfactory-related roles related to signal transduction mechanisms. Several mechanisms suggest that enzymatic inactivation could contribute to the signal termination process, such as odorant-degrading enzymes (ODEs). To date, a few ODEs have been identified and characterized in detail in insect herbivores, but little is known about aldehyde oxidases (AOXs); moreover, direct in vivo experimental evidence is needed. AOXs are a major family of metabolic enzymes that oxidize a variety of aromatic aldehydes, and they may also play a significant role in detoxification and degradation of environmental chemical cues. Here, we report on the identification and characterization of a novel cDNA encoding the putative odorant-degrading enzyme, PxylAOX3, from the antennae of the diamondback moth, (DBM), Plutella xylostella (L.) (Lepidoptera: Plutellidae). The purified recombinant protein showed a wide-range of substrate zymography oxidizing both sex pheromone compounds as well as plant-derived aldehydes with distinct activities. Our data suggest PxylAOX3 might be involved in the degradation of many structurally diverse aldehyde odorants. Furthermore, PxylAOX3 could participate in olfactory neuron protection by inactivation of redundant odorants and xenobiotic detoxification, making it a potential target for pesticide development as well.

Genomic divergences between the two polyphagous Spodoptera relatives provide cues for successful invasion of the fall armyworm

The fall armyworm Spodoptera frugiperda recently invaded China, ravaging crops in many provinces. Deciphering the possible genetic basics for its successful invasion is critical for innovative and specific control for this gluttonous pest. Here we generated comparative genomic analyses between S. frugiperda and its native relative, S. litura, which differs in host preference, locomotivity and production behavior. We demonstrated that S. frugiperda genes are enriched in taste sensory perception and nervous system, obviously different from those of S. litura. Potential host adaptation genes showed generally an elevated ratio of non-synonymous substitution rate to synonymous substitution rate, suggesting a faster evolution during the divergence of the two species. Focusing on these sets of genes, we identified 23 genes being under positive selection in S. frugiperda. Among them are two notable genes involved in sensory perception, gustatory receptor (GR) and an acetaldehyde oxidase, which are important for host detection in invasion and expansion processes. Another two genes are mitochondrial adenosine triphosphate synthase β subunit and ferritin heavy chain, which may be associated with the enhanced locomotivity and resistance, which fascinated long-distance migration needed for invasion and rapid expansion. Another interesting gene is chorion protein, in which positive selection sites in S. frugiperda were found and a replacement in one site is predicted to affect the protein function, which might be associated with competent reproductivity in S. frugiperda to ensure genetic resources for expansion.

Expression, Affinity, and Functional Characterization of the Specific Binding of Two Putative Pheromone-Binding Proteins in the Omnivorous German Cockroach Blattella germanica

The German cockroach Blattella germanica (L.) is an important pest in medical, veterinary, and public health. Studies on the olfaction mechanism of hemimetabolous insects have rarely been reported, especially in cockroaches. Pheromone-binding proteins (PBPs) play a vital role in insect sex pheromone recognition, which solubilize and carry the hydrophobic pheromonal compounds through the antennal lymph to receptors. In this study, two potential PBPs (BgerOBP26 and BgerOBP40) were identified on the basis of their biased expression in male antennae using tissue transcriptome data and verified by the quantitative real-time polymerase chain reaction approach. We then expressed and purified the two identified odorant-binding proteins (OBPs) using the Escherichia coli expression system and affinity purification. In vitro binding studies showed that the two OBPs display stronger binding affinities to the female volatile sex pheromone blattellaquinone than to its analogues and contact sex pheromone components. Finally, three-dimensional modeling of the two OBPs and dock conformation with sex pheromone molecules showed BgerOBP26 has a larger odorant cavity and more conservative active amino acid residues than BgerOBP40. These results illuminated the binding characteristics of potential PBPs of B. germanica, which could lay the groundwork for improved understanding of many aspects of the chemical ecology of B. germanica. Moreover, this information complements the understanding of the olfactory molecular mechanism in cockroaches and provides potential gene targets for B. germanica control.

Identification of genes involved in sex pheromone biosynthesis and metabolic pathway in the Chinese oak silkworm, Antheraea pernyi

The Chinese oak silkworm, Antheraea pernyi, has not only been semi-domesticated as an important economical insect but also used for genetic research. The female moths of A. pernyi employ a pheromone blend containing (E,Z)-6,11-hexadecadienal (E6,Z11-16:Ald), (E,Z)-6,11-hexadecadienyl acetate (E6,Z11-16:OAc), and (E,Z)-4,9-tetradecadienyl acetate (E4,Z9-14:OAc). While its biosynthesis pathway is largely unknown. By deep sequencing and de novo assembly of sex pheromone gland (PG) transcriptome, we identified 141 candidate genes that are putatively related to pheromone biosynthesis, degradation, and chemoreception in A. pernyi. Gene expression patterns and phylogenetic analysis revealed that two desaturases (AperDES1 and 2), two fatty acid reductase (AperFAR1 and 2), and three acetyltransferase genes (AperACT1, 2 and 3) showed PG-biased or specific expression and were phylogenetically related to genes known to be involved in pheromone synthesis in other species. Furthermore, two carboxylesterases (AperCOE6 and 11) and two chemosensory protein (AperCSP1 and 6) were also expressed specifically or predominantly in the PGs, which might be related to sex pheromone degradation and transportation, respectively. Based on these results, the sex pheromone biosynthesis and metabolic pathway was proposed in A. pernyi. This study provides some crucial candidates for further functional elucidation, and may be used for interfering sexual communication in other Saturniidae pests.

Alcohol dehydrogenase 5 of Helicoverpa armigera interacts with the CYP6B6 promoter in response to 2-tridecanone

Alcohol dehydrogenase 5 (ADH5) is a member of medium-chain dehydrogenase/reductase family and takes part in cellular formaldehyde and S-nitrosoglutathione metabolic network. 2-tridecanone (2-TD) is a toxic compound in many Solanaceae crops to defend against a variety of herbivory insects. In the broader context of insect development and pest control strategies, this study investigates how a new ADH5 from Helicoverpa armigera (HaADH5) regulates the expression of CYP6B6, a gene involved in molting and metamorphosis, in response to 2-TD treatment. Cloning of the HaADH5 complementary DNA sequence revealed that its 1002 bp open reading frame encodes 334 amino acids with a predicted molecular weight of 36.5 kD. HaADH5 protein was purified in the Escherichia coli Transetta (pET32a-HaADH5) strain using a prokaryotic expression system. The ability of HaADH5 protein to interact with the 2-TD responsive region within the promoter of CYP6B6 was confirmed by an in vitro electrophoretic mobility shift assay and transcription activity validation in yeast. Finally, the expression levels of both HaADH5 and CYP6B6 were found to be significantly decreased in the midgut of 6th instar larvae after 48 h of treatment with 10 mg/g 2-TD artificial diet. These results indicate that upon 2-TD treatment of cotton bollworm, HaADH5 regulates the expression of CYP6B6 by interacting with its promoter. As HaADH5 regulation of CYP6B6 expression may contribute to the larval xenobiotic detoxification, molting and metamorphosis, HaADH5 is a candidate target for controlling the growth and development of cotton bollworm.

Pheromone biosynthetic pathway and chemoreception proteins in sex pheromone gland of Eogystia hippophaecolus

The moth Eogystia hippophaecolus (Hua et al.) is a major threat to sea buckthorn plantations in China. Specific and highly efficient artificial sex pheromone traps have been developed and used to control this pest species. However, the biosynthesis of sex pheromones Z7-14: Ac and E3-14:Ac remains poorly understood. We investigated the female pheromone gland transcriptome of E. hippophaecolus and identified two pheromone biosynthesis-activating neuropeptides (PBANs), two pheromone biosynthesis-activating neuropeptide receptors (PBANrs), five acetyl-CoA carboxylases (ACCs), six fatty acid synthases (FASs), 16 Acyl-CoA desaturases (DESs), 26 reductases (REDs), 13 acetyltransferases (ACTs), one fatty acid transport protein (FATP), one acyl-CoA-binding protein (ACBP), and five elongation of very long-chain fatty acid proteins (ELOs) in pheromone biosynthesis pathways. Additionally, we identified 11 odorant-degrading enzymes (ODEs) and 16 odorant-binding proteins (OBPs), 14 chemosensory proteins (CSPs), two sensory neuron membrane proteins (SNMPs), three odorant receptors (ORs), seven ionotropic receptors (IRs), and six gustatory receptors (GRs). 77 unigenes involved in female pheromone biosynthesis, 31 chemoreception proteins and 11 odorant degradation enzymes were identified, which provided insight into the regulation of the pheromone components and pheromone recognition in the sex pheromone gland, and knowledge pertinent to new integrated pest management strategy of interference pheromone biosynthesis and recognition.

Overexpression of UDP-glycosyltransferase potentially involved in insecticide resistance in Aphis gossypii Glover collected from Bt cotton fields in China

BACKGROUND

The cotton aphid Aphis gossypii Glover is one of the most destructive insect pests. It has evolved resistance to numerous insecticides around the world due to the application of insecticides. Uridine diphosphate (UDP)-glycosyltransferases (UGTs) have been reported to potentially facilitate the detoxification process of imidacloprid and thiamethoxam in A. gossypii.

RESULTS

In this study, the field populations of A. gossypii developed different levels of resistance to multiple insecticides. A UGT inhibitor, 5-nitrouracil, dramatically increased the toxicity of acetamiprid in resistant populations, moderately increased the toxicity of sulfoxaflor in the imidacloprid susceptible (IMI_S) population, and populations from Yuncheng in Shanxi Province (SXYC) and Jingzhou in Hubei Province (HBJZ), and increased the toxicity of bifenthrin in the IMI_S and HBJZ populations, but there was no synergism on omethoate or carbosulfan. Quantitative real-time PCR analysis revealed that UGT344B4 and UGT344C7 were overexpressed in all field populations, and UGT344N4 was overexpressed in the SDBZ and HBZJ populations. Furthermore, the suppression of UGT344B4 or UGT344C7 by RNA interference significantly increased the susceptibility to bifenthrin in the IMI_S population and the susceptibility to sulfoxaflor in the SXYC population.

CONCLUSION

These results suggested that UGTs are potentially involved in the detoxification of neonicotinoid, sulfoximine, and pyrethroid insecticides in A. gossypii. Furthermore, the overexpression of UGTs could be associated with insecticide resistance in field populations of A. gossypii. The results might be helpful for the management of insecticide resistance in field populations of A. gossypii. © 2019 Society of Chemical Industry.

Comparative transcriptomics of the pheromone glands provides new insights into the differentiation of sex pheromone between two host populations of Chilo suppressalis

Reproductive isolation between different host populations is often based on intraspecific sex pheromone differences. The mechanisms underlying these differences have not been thoroughly elucidated to date. Previous studies suggested that Chilo suppressalis has differentiated into rice and water-oat host populations, and these two populations manifest clear differences in sex pheromone titer and mating rhythm. Hence, this moth is an ideal model to investigate the endogenous mechanisms of intraspecific reproductive isolation. Here, we identified a series of putative genes associated with sex pheromone biosynthesis based on the C. suppressalis pheromone gland transcriptome data. Transcripts of most genes were at higher level in the rice population. Then we obtained 11 pivotal differentially expressed genes (DEGs). The expression levels of these DEGs exhibited a distinct increase in the rice population. Moreover, we also observed the expression rhythm of these DEGs is discrepant between two host populations. Our study offers a new understanding to elucidate the mechanisms of intraspecific reproductive isolation.

Modulation of Sex Pheromone Discrimination by A UDP-Glycosyltransferase in Drosophila melanogaster

The detection and processing of chemical stimuli involve coordinated neuronal networks that process sensory information. This allows animals, such as the model species Drosophila melanogaster, to detect food sources and to choose a potential mate. In peripheral olfactory tissues, several classes of proteins are acting to modulate the detection of chemosensory signals. This includes odorant-binding proteins together with odorant-degrading enzymes (ODEs). These enzymes, which primarily act to eliminate toxic compounds from the whole organism also modulate chemodetection. ODEs are thought to neutralize the stimulus molecule concurrently to its detection, avoiding receptor saturation thus allowing chemosensory neurons to respond to the next stimulus. Here, we show that one UDP-glycosyltransferase (UGT36E1) expressed in D. melanogaster antennal olfactory sensory neurons (OSNs) is involved in sex pheromone discrimination. UGT36E1 overexpression caused by an insertion mutation affected male behavioral ability to discriminate sex pheromones while it increased OSN electrophysiological activity to male pheromones. Reciprocally, the decreased expression of UGT36E1, controlled by an RNAi transgene, improved male ability to discriminate sex pheromones whereas it decreased electrophysiological activity in the relevant OSNs. When we combined the two genotypes (mutation and RNAi), we restored wild-type-like levels both for the behavioral discrimination and UGT36E1 expression. Taken together, our results strongly suggest that this UGT plays a pivotal role in Drosophila pheromonal detection.

Transcriptional comparison between pheromone gland-ovipositor and tarsi in the corn earworm moth Helicoverpa zea

The corn earworm, Helicoverpa zea, utilizes (Z)-11-hexadecenal as the major sex pheromone component. The saturated fatty acid derivative hexadecanal is also found in the pheromone gland and recently a large amount (0.5-1.5 μg) was found in male tarsi with lower amounts (0.05-0.5 μg) in female tarsi. In this study, we compared the transcriptome between female pheromone glands (including the ovipositor) and female and male tarsi to identify differences between these tissues, particularly the genes involved in sex pheromone biosynthesis and chemosensation. We found transcripts encoding 9 fatty acyl-CoA desaturases, 20 fatty acyl-CoA reductases, 8 alcohol oxidases, some G protein-coupled receptors and many transcripts involved in signal transduction and pheromone transportation. Also we found gustatory and olfactory receptors associated with the tarsi and ovipositor. Differential expression analysis showed that there were many genes differentially expressed between tissues, including the candidate desaturases, fatty acyl-CoA reductases, and alcohol oxidases. We discuss how some of these genes produce proteins that could be involved in the biosynthesis of hexadecanal in tarsi and (Z)-11-hexadecenal in the pheromone gland and the possible role of proteins in chemosensation of the tarsi and ovipositor.

Synthesis and bioactivity studies of sex pheromone analogs of the diamond back moth, Plutella xylostella

BACKGROUND

Diamondback moth, Plutella xylostella L., is a very important pest of cruciferous vegetables causing excessive economic losses worldwide. Bioactivities of halo-, diazo-, and cyclopropane acetates of P. xylostella sex pheromone have been evaluated using electrophysiology and enzyme inhibition assays.

RESULTS

A total of 23 sex pheromone analogs of P. xylostella were designed and synthesized and the result shows that (11Z)-hexadec-11-en-1-yl 2,2,2-trifluoroacetate, (11Z)-hexadec-11-en-1-yl 2,2,3,3,3-pentafluoropropanoate, and (11Z)-hexadec-11-en-1-yl trifluoromethanesulfonate elicited potential inhibitory effects at all doses tested in the electrophysiology and enzyme inhibition assays. Interference of locating the sex pheromone source was found strongest when these three analogs were mixed with the sex pheromone at a 10:1 ratio. In addition, field test showed that the rate of mating disruption was over 90% when (11Z)-hexadec-11-en-1-yl 2,2,2-trifluoroacetate or (11Z)-hexadec-11-en-1-yl 2,2,3,3,3-pentafluoropropanoate was mixed with the sex pheromone at a 10:1 ratio.

CONCLUSION

Two sex pheromone antagonists were screen out by electrophysiology, enzyme inhibition assays, wind tunnel and field tests. We believe that these antagonists could be used to establish a novel eco-friendly measure to control P. xylostella and provide evidence for clarifying the specific functions and molecular mechanisms of sex pheromone antagonists. © 2018 Society of Chemical Industry.

Gypsy moth genome provides insights into flight capability and virus-host interactions

Since its accidental introduction to Massachusetts in the late 1800s, the European gypsy moth (EGM; Lymantria dispar dispar) has become a major defoliator in North American forests. However, in part because females are flightless, the spread of the EGM across the United States and Canada has been relatively slow over the past 150 years. In contrast, females of the Asian gypsy moth (AGM; Lymantria dispar asiatica) subspecies have fully developed wings and can fly, thereby posing a serious economic threat if populations are established in North America. To explore the genetic determinants of these phenotypic differences, we sequenced and annotated a draft genome of L. dispar and used it to identify genetic variation between EGM and AGM populations. The 865-Mb gypsy moth genome is the largest Lepidoptera genome sequenced to date and encodes ∼13,300 proteins. Gene ontology analyses of EGM and AGM samples revealed divergence between these populations in genes enriched for several gene ontology categories related to muscle adaptation, chemosensory communication, detoxification of food plant foliage, and immunity. These genetic differences likely contribute to variations in flight ability, chemical sensing, and pathogen interactions among EGM and AGM populations. Finally, we use our new genomic and transcriptomic tools to provide insights into genome-wide gene-expression changes of the gypsy moth after viral infection. Characterizing the immunological response of gypsy moths to virus infection may aid in the improvement of virus-based bioinsecticides currently used to control larval populations.

Glutathione-S-Transferases in the Olfactory Organ of the Noctuid Moth Spodoptera littoralis, Diversity and Conservation of Chemosensory Clades

Glutathione-S-transferases (GSTs) are conjugating enzymes involved in the detoxification of a wide range of xenobiotic compounds. The expression of GSTs as well as their activities have been also highlighted in the olfactory organs of several species, including insects, where they could play a role in the signal termination and in odorant clearance. Using a transcriptomic approach, we identified 33 putative GSTs expressed in the antennae of the cotton leafworm Spodoptera littoralis. We established their expression patterns and revealed four olfactory-enriched genes in adults. In order to investigate the evolution of antennal GST repertoires in moths, we re-annotated antennal transcripts corresponding to GSTs in two moth and one coleopteran species. We performed a large phylogenetic analysis that revealed an unsuspected structural—and potentially functional—diversity of GSTs within the olfactory organ of insects. This led us to identify a conserved clade containing most of the already identified antennal-specific and antennal-enriched GSTs from moths. In addition, for all the sequences from this clade, we were able to identify a signal peptide, which is an unusual structural feature for GSTs. Taken together, these data highlight the diversity and evolution of GSTs in the olfactory organ of a pest species and more generally in the olfactory system of moths, and also the conservation of putative extracellular members across multiple insect orders.

Expression Profiles and Functional Characterization of Two Odorant-Binding Proteins From the Apple Buprestid Beetle Agrilus mali (Coleoptera: Buprestidae)

The apple buprestid beetle, Agrilus mali Matsumura (Coleoptera: Buprestidae), can respond to various volatiles, but the underlying mechanism of odorant perception for this insect is poorly understood. Here, we cloned A. mali's odorant-binding proteins 3 (AmalOBP3) and 8 (AmalOBP8) and characterized their expression patterns and binding profiles. Sequence and phylogenetic analyses showed that AmalOBP3 and AmalOBP8 were distributed in the classic and minus-C OBP subfamily, respectively. AmalOBP3 was specifically and abundantly expressed in antennae of both sexes. AmalOBP8 displayed high transcript levels in antennae of both sexes, abdomens of males, and wings of both sexes. Both AmalOBPs exhibited much higher expression in male antennae than in female antennae, suggesting that they could be important in perception of male-specific olfactory cues (e.g., some sex pheromones). Out of the 40 odorant ligands tested, AmalOBP3 and AmalOBP8 bound to 15 and 21 different odorants, respectively, indicating a distinct and selective binding profile for them. Both AmalOBPs seemed to have very strong binding affinity to aliphatic alcohols and aldehydes with 12 to 15 carbon atoms. Alcohols, esters, and terpenoids were more likely to be good ligands for both AmalOBPs than aldehydes and alkanes. Together with its broad expression in different tissues, strong binding with higher numbers of putative ligands for AmalOPB8 means that this protein can have more extensive functional roles in chemosensation of A. mali. Our results provide insights into the molecular basis of chemosensation in A. mali, as well as a basis for developing detection, monitoring, and management tools for this serious pest.

Analysis of natural female post-mating responses of Anopheles gambiae and Anopheles coluzzii unravels similarities and differences in their reproductive ecology

Anopheles gambiae and An. coluzzii, the two most important malaria vectors in sub-Saharan Africa, are recently radiated sibling species that are reproductively isolated even in areas of sympatry. In females from these species, sexual transfer of male accessory gland products, including the steroid hormone 20-hydroxyecdysone (20E), induces vast behavioral, physiological, and transcriptional changes that profoundly shape their post-mating ecology, and that may have contributed to the insurgence of post-mating, prezygotic reproductive barriers. As these barriers can be detected by studying transcriptional changes induced by mating, we set out to analyze the post-mating response of An. gambiae and An. coluzzii females captured in natural mating swarms in Burkina Faso. While the molecular pathways shaping short- and long-term mating-induced changes are largely conserved in females from the two species, we unravel significant inter-specific differences that suggest divergent regulation of key reproductive processes such as egg development, processing of seminal secretion, and mating behavior, that may have played a role in reproductive isolation. Interestingly, a number of these changes occur in genes previously shown to be regulated by the sexual transfer of 20E and may be due to divergent utilization of this steroid hormone in the two species.

Identification of Chemosensory Genes Based on the Transcriptomic Analysis of Six Different Chemosensory Organs in Spodoptera exigua

Insects have a complex chemosensory system that accurately perceives external chemicals and plays a pivotal role in many insect life activities. Thus, the study of the chemosensory mechanism has become an important research topic in entomology. Spodoptera exigua Hübner (Lepidoptera: Noctuidae) is a major agricultural polyphagous pest that causes significant agricultural economic losses worldwide. However, except for a few genes that have been discovered, its olfactory and gustatory mechanisms remain uncertain. In the present study, we acquired 144,479 unigenes of S. exigua by assembling 65.81 giga base reads from 6 chemosensory organs (female and male antennae, female and male proboscises, and female and male labial palps), and identified many differentially expressed genes in the gustatory and olfactory organs. Analysis of the transcriptome data obtained 159 putative chemosensory genes, including 24 odorant binding proteins (OBPs; 3 were new), 19 chemosensory proteins (4 were new), 64 odorant receptors (57 were new), 22 ionotropic receptors (16 were new), and 30 new gustatory receptors. Phylogenetic analyses of all genes and SexiGRs expression patterns using quantitative real-time polymerase chain reactions were investigated. Our results found that several of these genes had differential expression features in the olfactory organs compared to the gustatory organs that might play crucial roles in the chemosensory system of S. exigua, and could be utilized as targets for future functional studies to assist in the interpretation of the molecular mechanism of the system. They could also be used for developing novel behavioral disturbance agents to control the population of the moths in the future.

Gene Identification of Pheromone Gland Genes Involved in Type II Sex Pheromone Biosynthesis and Transportation in Female Tea Pest Ectropis grisescens

Moths can biosynthesize sex pheromones in the female sex pheromone glands (PGs) and can distinguish species-specific sex pheromones using their antennae. However, the biosynthesis and transportation mechanism for Type II sex pheromone components has rarely been documented in moths. In this study, we constructed a massive PG transcriptome database (14.72 Gb) from a moth species, Ectropis grisescens, which uses type II sex pheromones and is a major tea pest in China. We further identified putative sex pheromone biosynthesis and transportation-related unigenes: 111 cytochrome P450 monooxygenases (CYPs), 25 odorant-binding proteins (OBPs), and 20 chemosensory proteins (CSPs). Tissue expression and phylogenetic tree analyses showed that one CYP (EgriCYP341-fragment3), one OBP (EgriOBP4), and one CSP (EgriCSP10) gene displayed an enriched expression in the PGs, and that EgriOBP2, 3, and 25 are clustered in the moth pheromone-binding protein clade. We considered these our candidate genes. Our results yielded large-scale PG sequence information for further functional studies.

A Transcriptome Survey Spanning Life Stages and Sexes of the Harlequin Bug, Murgantia histrionica

The harlequin bug, Murgantia histrionica (Hahn), is an agricultural pest in the continental United States, particularly in southern states. Reliable gene sequence data are especially useful to the development of species-specific, environmentally friendly molecular biopesticides and effective biolures for this insect. Here, mRNAs were sampled from whole insects at the 2nd and 4th nymphal instars, as well as sexed adults, and sequenced using Illumina RNA-Seq technology. A global assembly of these data identified 72,540 putative unique transcripts bearing high levels of similarity to transcripts identified in other taxa, with over 99% of conserved single-copy orthologs among insects being detected. Gene ontology and protein family analyses were conducted to explore the functional potential of the harlequin bug's gene repertoire, and phylogenetic analyses were conducted on gene families germane to xenobiotic detoxification, including glutathione S-transferases, carboxylesterases and cytochrome P450s. Genic content in harlequin bug was compared with that of the closely related invasive pest, the brown marmorated stink bug, Halyomorpha halys (Stål). Quantitative analyses of harlequin bug gene expression levels, experimentally validated using quantitative real-time PCR, identified genes differentially expressed between life stages and/or sexes.