Transcriptome analysis of the sex pheromone gland of the noctuid moth Heliothis virescens

Identification of attractants for adult Spodoptera litura based on the interaction between odorant-binding protein 34 and host volatiles

Odorant-binding proteins (OBPs) play key roles in host plant location by insects, and can accordingly serve as important targets for the development of attractants. In this study, we detected the high expression of SlitOBP34 in male antennae of Spodoptera litura. Subsequently, the fluorescence competitive binding experiments displayed that the SlitOBP34 protein has binding affinity for different ligands. Then, protein-ligand interaction analyses found the presence of six amino acid residues may serve as key recognition sites. Further electroantennographic and biobehavioral assessments revealed that the electrophysiological responses of male antennae were evoked in response to stimulation with the six identified host volatiles, and that these volatiles attracted male moths to varying extents. Notably, low concentrations of benzaldehyde, 1-hexanol, and cis-3-hexenyl acetate were found to have significant attractant effects on male moths, thereby identifying these three host volatiles as potential candidates for the development of male attractants. These findings advance our current understanding of the olfactory-encoded mechanisms of host plants selection in S. litura and have enabled us to develop novel adult attractants for controlling the pest in the future.

Transcriptome-informed identification and characterization of Planococcus citri cis- and trans-isoprenyl diphosphate synthase genes

Insect physiology and reproduction depend on several terpenoid compounds, whose biosynthesis is mainly unknown. One enigmatic group of insect monoterpenoids are mealybug sex pheromones, presumably resulting from the irregular coupling activity of unidentified isoprenyl diphosphate synthases (IDSs). Here, we performed a comprehensive search for IDS coding sequences of the pest mealybug Planococcus citri. We queried the available genomic and newly generated short- and long-read P. citri transcriptomic data and identified 18 putative IDS genes, whose phylogenetic analysis indicates several gene family expansion events. In vitro testing confirmed regular short-chain coupling activity with five gene products. With the candidate with highest IDS activity, we also detected low amounts of irregular coupling products, and determined amino acid residues important for chain-length preference and irregular coupling activity. This work therefore provides an important foundation for deciphering terpenoid biosynthesis in mealybugs, including the sex pheromone biosynthesis in P. citri.

Functional Characterization of Pheromone Receptors in the Beet Webworm, Loxostege sticticalis (Lepidoptera: Pyralidae)

Lepidopteran insects mainly rely on sex pheromones to complete sexual communications. Pheromone receptors (PRs) are expressed on the olfactory receptor neurons (ORNs) of the sensilla trichodea and play an essential role in sexual communication. Despite extensive investigations into the mechanisms of peripheral recognition of sex pheromones in Lepidoptera, knowledge about these mechanisms in L. sticticalis remains limited. In this study, five candidate LstiPRs were analyzed in a phylogenetic tree with those of other Lepidopteran insects. Electroantennography (EAG) assays showed that the major sex pheromone component E11-14:OAc elicited a stronger antennal response than other compounds in male moths. Moreover, two types of neurons in sensilla trichodea were classified by single sensillum recordings, of which the "a" neuron specifically responded to E11-14:OAc. Five candidate PRs were functionally assayed by the heterologous expression system of Xenopus oocytes, and LstiPR2 responded to the major sex pheromone E11-14:OAc. Our findings suggest that LstiPR2 is a PR sensitive to L. sticticalis's major sex pheromone compound, E11-14:OAc. Furthermore, this study offers valuable insights into the sexual communication behavior of L. sticticalis, forming a foundation for further analysis of the species' central nervous system.

Identification and Expression Profiles of Candidate Sex Pheromone Biosynthesis Genes by the Transcriptome Analysis of Sex Pheromone Glands in Spodoptera litura and Spodoptera exigua

Like many insects, females of the Noctuid moth Spodoptera litura and Spodoptera exigua release chemical signals to attract males from a long distance for successful mating. In this study, 98 and 86 genes related to the sex pheromone biosynthesis of S. litura and S. exigua were identified. The tissue expression profiles of highly expressed genes in sex pheromone glands (PGs) were further examined by real-time quantitative polymerase chain reaction. The results displayed that only SlitDes5 and SexiDes5 gene were specifically and significantly overexpressed in the PGs of S. litura and S. exigua. The functional study of SlitDes5 gene showed that RNA interference reduced its expression level by 49.42%. In addition, the content of the sex pheromones of S. litura, Z9E11-14:OAc, Z9E12-14:OAc, E11-14:OAc, and Z9-14:OAc, decreased by 41.98% on average. Our findings provide a basis for better understanding the key genes that affect the biosynthesis of sex pheromones and for determining potential gene targets for pest control strategies.

Genome-wide identification of FAR gene family and functional analysis of NlFAR10 during embryogenesis in the brown planthopper Nilaparvata lugens

Fatty acyl-CoA reductases (FARs) catalyze the synthesis of fatty alcohols from corresponding fatty acid precursors in organisms. However, the function of FARs in insect fecundity and embryogenesis remains largely unclear. Here, a total of 22 putative FAR proteins were identified in the brown planthopper Nilaparvata lugens, a hemipteran insect pest of rice, and most of them were highly expressed in embryonic stages. Among them, NlFAR10 was specifically and highly expressed in the later embryogenesis, but was promiscuously expressed in tissues of adults. The heterologously expressed NlFAR10 was able to produce the intermediate fatty acid alcohols from the corresponding acyl-CoA precursors. When NlFAR10 was silenced through RNAi in vivo, the embryogenesis was obviously inhibited, resulting in low hatching rates. Moreover, the metabolome analyses indicated that loss of NlFAR10 affected lipid metabolism and purine metabolism during embryogenesis. To the best of our knowledge, this is the first report of a FAR member affecting insect embryogenesis, thus providing a new target for future pest management.

Genome-Wide Identification and Expression Profiling of Candidate Sex Pheromone Biosynthesis Genes in the Fall Armyworm (Spodoptera frugiperda)

Spodoptera frugiperda is an agricultural pest causing substantial damage and losses to commercial crops. Sex pheromones are critical for successful mating in Lepidoptera and have been used for monitoring and control of many pest species. The sex pheromone of S. frugiperda is known, but the genes involved in its biosynthesis have not been identified. We systematically studied 99 candidate sex pheromone genes in the genome of S. frugiperda including 1 acetyl-CoA carboxylase (ACC), 11 fatty acid synthases (FASs), 17 desaturases (DESs), 4 fatty acid transport proteins (FATPs), 29 fatty acyl-CoA reductases (FARs), 17 acetyl-CoA acetyltransferases (ACTs), 5 acyl-CoA dehydrogenase (ACDs), 3 enoyl-CoA hydratases (ECHs), 3 hydroxyacyl-CoA dehydrogenases (HCDs), 6 ethyl-CoA thiolases (KCTs), and 3 acyl-CoA-binding proteins (ACBPs). Based on the comparative transcriptome results, we found 22 candidate sex pheromone biosynthesis genes predominately expressed in pheromone glands (PGs) than abdomens without PGs including SfruFAS4, SfruFATP3, SfruACD5, SfruKCT3, SfruDES2, SfruDES5, SfruDES11, SfruDES13, SfruFAR1, SfruFAR2, SfruFAR3, SfruFAR6, SfruFAR7, SfruFAR8, SfruFAR9, SfruFAR10, SfruFAR11, SfruFAR14, SfruFAR16, SfruFAR29, SfruACT6, and SfruACT10. A combination of phylogenetic and tissue-specific transcriptomic analyses indicated that SfruDES5, SfruDES11, SfruFAR2, SfruFAR3, and SfruFAR9 may be key genes involved in the sex pheromone synthesis of S. frugiperda. Our results could provide a theoretical basis for understanding the molecular mechanisms of sex pheromone biosynthesis in S. frugiperda, and also provide new targets for developing novel pest control methods based on disrupting sexual communication.

Suppression of pheromone biosynthesis and mating behavior by RNA interference of pheromone gland-specific fatty acyl reductase in Maruca vitrata

In moths, various enzymes, such as fatty acid synthases, fatty acyl desaturases, and fatty acyl reductases (FARs), are involved in pheromone biosynthesis. In particular, pheromone gland-specific FAR (pgFAR) plays an important role in converting the functional group from carboxylic to alcohol during pheromone biosynthesis. A novel pgFAR of Maruca vitrata, Mvi-pgFAR, was identified through transcriptome sequencing of its pheromone gland. To investigate the involvement of Mvi-pgFAR in pheromone biosynthesis, Mvi-pgFAR was cloned from the pheromone gland and suppressed by RNA interference (RNAi). Mvi-pgFAR harbored several conserved motifs related to NAD(P)H-binding, N-glycosylation, and adenosine / guanosine triphosphate binding. Phylogenetic analysis revealed that Mvi-pgFAR with other lepidopteran pgFARs formed an independent clade. Mvi-pgFAR was specifically expressed only in the pheromone gland. Quantitative real-time polymerase chain reaction showed that the diurnal expression levels of Mvi-pgFAR in the pheromone gland were the highest at 2 h before the scotophase. After primarily confirming Mvi-pgFAR suppression by RNAi, (E,E)-10,12-hexadecadienal (E10E12-16:Ald), a major sex pheromone component, was quantified by gas chromatography. When Mvi-pgFAR was successfully suppressed, E10E12-16:Ald production was reduced by up to half of that of the control, and the mating rate was subsequently decreased. Our results demonstrate that Mvi-pgFAR downregulation can suppress mating behavior by changing the relative sex pheromone component ratio, suggesting that Mvi-pgFAR can be used as a novel control target.

Differential expression profiling of Oxycarenus laetus Kirby (Hemiptera: Lygaeidae) upon exposure to gossypol

BACKGROUND

Gossypium hirsutum seeds are rich in gossypol. In addition to its diverse beneficial properties, it is a known anti-fertility inducing agent in humans. Oxycarenus laetus feeds on the cottonseeds and yet its courtship, mating and reproduction is unaffected.

METHODS AND RESULTS

In this study, we performed a transcriptome profiling of O. laetus fed on Abutilon indicum (AB-no gossypol), G. hirsutum (GH-natural gossypol) and 1400 ppm commercial gossypol-soaked GH seeds (GHGO). Illumina NextSeq-500 paired-end 75 bp reads were generated and de novo assembled (48,214 genes) to identify the differentially expressed transcripts (DET) between the samples. Gene enrichment, KEGG pathway and cluster profiling of the DETs resulted in the identification of vital genes involved in the detoxification, pheromone biosynthesis, cuticle protein in the GHGO sample. Cyp4C1, Cyp6a13, Cyp6a14, Cyp4g15, Cyp4em8, Cyp303a1 were the detoxification related genes identified. Similarly, SDR dehydrogenase family 11 and fatty acid synthase in pheromone biosynthesis and cuticle proteins (RR1 and RR2) coding transcripts were found to be differentially expressed.

CONCLUSION

This is the first study to report the expression of genes induced by gossypol in O. laetus. Based on the findings from the DET analysis, we conclude that the detoxification related genes of gossypol treated samples were affected.

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.

Identification of Odorant-Binding and Chemosensory Protein Genes in Mythimna separata Adult Brains Using Transcriptome Analyses

Large numbers of chemosensory genes have been identified in the peripheral sensory organs of the pest Mythimna separata (Walker) to increase our understanding of chemoreception-related molecular mechanisms and to identify molecular targets for pest control. Chemosensory-related genes are expressed in various tissues, including non-sensory organs, and they play diverse roles. To better understand the functions of chemosensory-related genes in non-sensory organs, transcriptomic analyses of M. separata brains were performed. In total, 29 odorant-binding proteins (OBPs) and 16 chemosensory proteins (CSPs) putative genes were identified in the transcriptomic data set. The further examination of sex- and tissue-specific expression using RT-PCR suggested that eight OBPs (OBP5, -7, -11, -13, -16, -18, -21, and -24) and eight CSPs (CSP2–4, -8, CSP10–12, and -15) genes were expressed in the brain. Furthermore, bands representing most OBPs and CSPs could be detected in antennae, except for a few that underwent sex-biased expression in abdomens, legs, or wings. An RT-qPCR analysis of the expression profiles of six OBPs (OBP3–5, -9, -10, and -16) and two CSPs (CSP3 and CSP4) in different tissues and sexes indicated that OBP16 was highly expressed in male brain, and CSP3 and CSP4 were female-biased and highly expressed in brain. The expression levels of OBP5 and OBP10 in brain were not significantly different between the sexes. The findings expand our current understanding of the expression patterns of OBPs and CSPs in M. separata sensory and non-sensory tissues. These results provide valuable reference data for exploring novel functions of OBPs and CSPs in M. separata and may help in developing effective biological control strategies for managing this pest by exploring novel molecular targets.

Transcriptome analysis identifies candidate genes in the biosynthetic pathway of sex pheromones from a zygaenid moth, Achelura yunnanensis (Lepidoptera: Zygaenidae)

In most moth species, sex pheromones responsible for mating and communication of both sexes are primarily produced by the pheromone glands (PGs) of female moths. Although the PG transcriptomes and pheromone production related genes from 24 moth species have been characterized, studies on the related information remain unknown in the Zygaenidae family. Here, we sequenced the PG transcriptome of a zygaenid moth, Achelura yunnanensis. Such the sequencing resulted in the yields of 47,632,610 clean reads that were assembled into 54,297 unigenes, coupled with RNA sequencing data from 12 other tissues. Based on the transcriptome, a total of 191 genes encoding pheromone biosynthesis and degradation enzymes were identified, 161 of which were predicted to have full-length sequences. A comparative analysis among 24 moth species of nine families indicated that the numbers of the genes were variable, ranging from 14 in two Grapholita species to 191 in A. yunnanensis. Phylogenetic analysis in parallel with the expression data highlighted some key genes, including three △9 and four △11 desaturases, four fatty acyl-CoA reductases (FARs) clustering in the pgFAR clade, and three significantly antennae-enriched aldehyde oxidases. An extensive tissue- and sex- expression profile revealed a broad distribution of the genes, in which 128 relatives were detected in the PGs and 127 in the antennae. This study reports, for the first time, the gene repertoires associated with the pheromone production in Zygaenidae, and provides a valuable resource for exploring putative roles of the PG-enriched genes in A. yunnanensis.

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.

Latest Developments in Insect Sex Pheromone Research and Its Application in Agricultural Pest Management

Since the first identification of the silkworm moth sex pheromone in 1959, significant research has been reported on identifying and unravelling the sex pheromone mechanisms of hundreds of insect species. In the past two decades, the number of research studies on new insect pheromones, pheromone biosynthesis, mode of action, peripheral olfactory and neural mechanisms, and their practical applications in Integrated Pest Management has increased dramatically. An interdisciplinary approach that uses the advances and new techniques in analytical chemistry, chemical ecology, neurophysiology, genetics, and evolutionary and molecular biology has helped us to better understand the pheromone perception mechanisms and its practical application in agricultural pest management. In this review, we present the most recent developments in pheromone research and its application in the past two decades.

Isolation, Identification, and Bioinformatic Analysis of Antibacterial Proteins and Peptides from Immunized Hemolymph of Red Palm Weevil Rhynchophorus ferrugineus

Red palm weevil (Rhynchophorus ferrugineus Olivier, 1791, Coleoptera: Curculionidae) is a destructive pest of palms, rapidly extending its native geographical range and causing large economic losses worldwide. The present work describes isolation, identification, and bioinformatic analysis of antibacterial proteins and peptides from the immunized hemolymph of this beetle. In total, 17 different bactericidal or bacteriostatic compounds were isolated via a series of high-pressure liquid chromatography steps, and their partial amino acid sequences were determined by N-terminal sequencing or by mass spectrometry. The bioinformatic analysis of the results facilitated identification and description of corresponding nucleotide coding sequences for each peptide and protein, based on the recently published R. ferrugineus transcriptome database. The identified compounds are represented by several well-known bactericidal factors: two peptides similar to defensins, one cecropin-A1-like peptide, and one attacin-B-like protein. Interestingly, we have also identified some unexpected compounds comprising five isoforms of pheromone-binding proteins as well as seven isoforms of odorant-binding proteins. The particular role of these factors in insect response to bacterial infection needs further investigation.

Genome-wide identification and analysis of the thiolase family in insects

Thiolases are important enzymes involved in lipid metabolism in both prokaryotes and eukaryotes, and are essential for a range of metabolic pathways, while, little is known for this important family in insects. To shed light on the evolutionary models and functional diversities of the thiolase family, 137 thiolase genes were identified in 20 representative insect genomes. They were mainly classified into five classes, namely cytosolic thiolase (CT-thiolase), T1-thiolase, T2-thiolase, trifunctional enzyme thiolase (TFE-thiolase), and sterol carrier protein 2 thiolase (SCP2-thiolase). The intron number and exon/intron structures of the thiolase genes reserve large diversification. Subcellular localization prediction indicated that all the thiolase proteins were mitochondrial, cytosolic, or peroxisomal enzymes. Four highly conserved sequence fingerprints were found in the insect thiolase proteins, including CxS-, NEAF-, GHP-, and CxGGGxG-motifs. Homology modeling indicated that insect thiolases share similar 3D structures with mammals, fishes, and microorganisms. In Bombyx mori, microarray data and reverse transcription-polymerase chain reaction (RT-PCR) analysis suggested that some thiolases might be involved in steroid metabolism, juvenile hormone (JH), and sex pheromone biosynthesis pathways. In general, sequence and structural characteristics were relatively conserved among insects, bacteria and vertebrates, while different classes of thiolases might have differentiation in specific functions and physiological processes. These results will provide an important foundation for future functional validation of insect thiolases.

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.

Chasing the Apomictic Factors in the Ranunculus auricomus Complex: Exploring Gene Expression Patterns in Microdissected Sexual and Apomictic Ovules

Apomixis, the asexual reproduction via seeds, is associated to polyploidy and hybridization. To identify possible signatures of apomixis, and possible candidate genes underlying the shift from sex to apomixis, microarray-based gene expression patterns of live microdissected ovules at four different developmental stages were compared between apomictic and sexual individuals of the Ranunculus auricomus complex. Following predictions from previous work on mechanisms underlying apomixis penetrance and expressivity in the genus, gene expression patterns were classified into three categories based on their relative expression in apomicts compared to their sexual parental ancestors. We found evidence of misregulation and differential gene expression between apomicts and sexuals, with the highest number of differences detected during meiosis progression and emergence of aposporous initial (AI) cells, a key developmental stage in the ovule of apomicts where a decision between divergent reproductive pathways takes place. While most of the differentially expressed genes (DEGs) could not be annotated, gene expression was classified into transgressive, parent of origin and ploidy effects. Genes related to gametogenesis and meiosis demonstrated patterns reflective of transgressive and genome dosage effects, which support the hypothesis of a dominant factor controlling apomixis in Ranunculus and modulated by secondary modifiers. Three genes with probable functions in sporogenesis and gametogenesis development are identified and characterized for future studies.

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.

Molecular cloning and comparative analysis of transcripts encoding chemosensory proteins from two plant bugs, Lygus lineolaris and Lygus hesperus

Chemosensory proteins (CSPs) are soluble carrier proteins typically characterized by a six-helix bundle structure joined by two disulfide bridges and a conserved Cys spacing pattern (C1-X_6-8 -C2-X_16-21 -C3-X₂ -C4). CSPs are functionally diverse with reported roles in chemosensation, immunity, development, and resistance. To expand our molecular understanding of CSP function in plant bugs, we used recently developed transcriptomic resources for Lygus lineolaris and Lygus hesperus to identify 17 and 14 CSP-like sequences, respectively. The Lygus CSPs are orthologous and share significant sequence identity with previously annotated CSPs. Three of the CSPs are predicted to deviate from the typical CSP structure with either five or seven helical segments rather than six. The seven helix CSP is further differentiated by an atypical C3-X₃ -C4 Cys spacing motif. Reverse transcriptase PCR-based profiling of CSP transcript abundance in adult L. lineolaris tissues revealed broad expression for most of the CSPs with antenna specific expression limited to a subset of the CSPs. Comparative sequence analyses and homology modeling suggest that variations in the amino acids that comprise the Lygus CSP binding pockets affect the size and nature of the ligands accommodated.

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.

Molecular evolution of gland cell types and chemical interactions in animals

Across the Metazoa, the emergence of new ecological interactions has been enabled by the repeated evolution of exocrine glands. Specialized glands have arisen recurrently and with great frequency, even in single genera or species, transforming how animals interact with their environment through trophic resource exploitation, pheromonal communication, chemical defense and parental care. The widespread convergent evolution of animal glands implies that exocrine secretory cells are a hotspot of metazoan cell type innovation. Each evolutionary origin of a novel gland involves a process of 'gland cell type assembly': the stitching together of unique biosynthesis pathways; coordinated changes in secretory systems to enable efficient chemical release; and transcriptional deployment of these machineries into cells constituting the gland. This molecular evolutionary process influences what types of compound a given species is capable of secreting, and, consequently, the kinds of ecological interactions that species can display. Here, we discuss what is known about the evolutionary assembly of gland cell types and propose a framework for how it may happen. We posit the existence of 'terminal selector' transcription factors that program gland function via regulatory recruitment of biosynthetic enzymes and secretory proteins. We suggest ancestral enzymes are initially co-opted into the novel gland, fostering pleiotropic conflict that drives enzyme duplication. This process has yielded the observed pattern of modular, gland-specific biosynthesis pathways optimized for manufacturing specific secretions. We anticipate that single-cell technologies and gene editing methods applicable in diverse species will transform the study of animal chemical interactions, revealing how gland cell types are assembled and functionally configured at a molecular level.

Identification of putative Type-I sex pheromone biosynthesis-related genes expressed in the female pheromone gland of Streltzoviella insularis

Species-specific sex pheromones play key roles in moth sexual communication. Although the general pathway of Type-I sex pheromone biosynthesis is well established, only a handful of genes encoding enzymes involved in this pathway have been characterized. Streltzoviella insularis is a destructive wood-boring pest of many street trees in China, and the female sex pheromone of this species comprises a blend of (Z)-3-tetradecenyl acetate, (E)-3-tetradecenyl acetate, and (Z)-5-dodecenyl acetate. This organism therefore provides an excellent model for research on the diversity of genes and molecular mechanisms involved in pheromone production. Herein, we assembled the pheromone gland transcriptome of S. insularis by next-generation sequencing and identified 74 genes encoding candidate key enzymes involved in the fatty acid biosynthesis, β-oxidation, and functional group modification. In addition, tissue expression patterns further showed that an acetyl-CoA carboxylase and two desaturases were highly expressed in the pheromone glands compared with the other tissues, indicating possible roles in S. insularis sex pheromone biosynthesis. Finally, we proposed putative S. insularis biosynthetic pathways for sex pheromone components and highlighted candidate genes. Our findings lay a solid foundation for understanding the molecular mechanisms underpinning S. insularis sex pheromone biosynthesis, and provide potential targets for disrupting chemical communication that could assist the development of novel pest control methods.

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.

Characterization and expression profiling of odorant-binding proteins in Anoplophora glabripennis Motsch

In insects, olfaction plays a critical role in locating hosts, recognizing mates, and selecting oviposition sites. The Asian long-horned beetle (Anoplophora glabripennis Motschulsky) feeds on 43 species of trees in 15 families, but its chemosensory mechanisms are poorly understood. Herein, genes encoding 61 odorant-binding proteins (OBPs) were identified from the published genome and our previous A. glabripennis transcriptomic data. To investigate their physiological functions, we performed expression profiling of all AglaOBPs in the antennae, legs, and maxillary palps of both sexes. Phylogenetic analysis clustered A. glabripennis OBPs into four subgroups, comprising 29 Minus-C OBPs, 15 Antennae-binding proteins (ABPIIs), 10 Classic OBPs, and one Plus-C OBP. 12 AglaOBP genes were expressed specifically in antennae, and AglaOBP3, AglaOBP18, AglaOBP21, AglaOBP33, AglaOBP41, AglaOBP45, and AglaOBP47 were particularly highly expressed in male antennae. These proteins may function in the detection of female sex pheromones. AglaOBP23 and AglaOBP44 were preferentially expressed in maxillary palps. Expression profiling suggests that many OBPs may be involved in olfaction and gustation, in addition to carrying hydrophobic molecules. The AglaOBPs family has acquired functional diversity concurrently with functional constraints, and further investigation could provide insight into the roles of OBPs in chemoreception.

Carabidae Semiochemistry: Current and Future Directions

Ground beetles (Carabidae) are recognized for their diverse, chemically-mediated defensive behaviors. Produced using a pair of pygidial glands, over 250 chemical constituents have been characterized across the family thus far, many of which are considered allomones. Over the past century, our knowledge of Carabidae exocrine chemistry has increased substantially, yet the role of these defensive compounds in mediating behavior other than repelling predators is largely unknown. It is also unclear whether non-defensive compounds produced by ground beetles mediate conspecific and heterospecific interactions, such as sex-aggregation pheromones or kairomones, respectively. Here we review the current state of non-exocrine Carabidae semiochemistry and behavioral research, discuss the importance of semiochemical research including but not limited to allomones, and describe next-generation methods for elucidating the underlying genetics and evolution of chemically-mediated behavior.

Sensilla Morphology and Complex Expression Pattern of Odorant Binding Proteins in the Vetch Aphid Megoura viciae (Hemiptera: Aphididae)

Chemoreception in insects is mediated by several components interacting at different levels and including odorant-binding proteins (OBPs). Although recent studies demonstrate that the function of OBPs cannot be restricted to an exclusively olfactory role, and that OBPs have been found also in organs generally not related to chemoreception, their feature of binding molecules remains undisputed. Studying the vetch aphid Megoura viciae (Buckton), we used a transcriptomic approach to identify ten OBPs in the antennae and we examined the ultrastructural morphology of sensilla and their distribution on the antennae, legs, mouthparts and cauda of wingless and winged adults by scanning electron microscopy (SEM). Three types of sensilla, trichoid, coeloconic and placoid, differently localized and distributed on antennae, mouthparts, legs and cauda, were described. The expression analysis of the ten OBPs was performed by RT-qPCR in the antennae and other body parts of the wingless adults and at different developmental stages and morphs. Five of the ten OBPs (MvicOBP1, MvicOBP3, MvicOBP6, MvicOBP7, and MvicOBP8), whose antibodies were already available, were selected for experiments of whole-mount immunolocalization on antennae, mouthparts, cornicles and cauda of adult aphids. Most of the ten OBPs were more expressed in antennae than in other body parts; MvicOBP1, MvicOBP3, MvicOBP6, MvicOBP7 were also immunolocalized in the sensilla on the antennae, suggesting a possible involvement of these proteins in chemoreception. MvicOBP6, MvicOBP7, MvicOBP8, MvicOBP9 were highly expressed in the heads and three of them (MvicOBP6, MvicOBP7, MvicOBP8) were immunolocalized in the sensilla on the mouthparts, supporting the hypothesis that also mouthparts may be involved in chemoreception. MvicOBP2, MvicOBP3, MvicOBP5, MvicOBP8 were highly expressed in the cornicles-cauda and two of them (MvicOBP3, MvicOBP8) were immunolocalized in cornicles and in cauda, suggesting a possible new function not related to chemoreception. Moreover, the response of M. viciae to different components of the alarm pheromone was assessed by behavioral assays on wingless adult morph; (-)-α-pinene and (+)-limonene were found to be the components mainly eliciting an alarm response. Taken together, our results represent a road map for subsequent in-depth analyses of the OBPs involved in several physiological functions in M. viciae, including chemoreception.

Pheromone Binding Protein EhipPBP1 Is Highly Enriched in the Male Antennae of the Seabuckthorn Carpenterworm and Is Binding to Sex Pheromone Components

The seabuckthorn carpenterworm moth Eogystia hippophaecolus is a major threat to seabuckthorn plantations, causing considerable ecological and economic losses in China. Transcriptomic analysis of E. hippophaecolus previously identified 137 olfactory proteins, including three pheromone-binding proteins (PBPs). We investigated the function of E. hippophaecolus PBP1 by studying its mRNA and protein expression profiles and its binding ability with different compounds. The highest levels of expression were in the antennae, particularly in males, with much lower levels of expression in the legs and external genitals. Recombinant PBP1 showed strong binding to sex-pheromone components, suggesting that antennal EhipPBP1 is involved in binding sex-pheromone components during pheromone communication.

Biotechnological potential of insect fatty acid-modifying enzymes

There are more than one million described insect species. This species richness is reflected in the diversity of insect metabolic processes. In particular, biosynthesis of secondary metabolites, such as defensive compounds and chemical signals, encompasses an extraordinarily wide range of chemicals that are generally unparalleled among natural products from other organisms. Insect genomes, transcriptomes and proteomes thus offer a valuable resource for discovery of novel enzymes with potential for biotechnological applications. Here, we focus on fatty acid (FA) metabolism-related enzymes, notably the fatty acyl desaturases and fatty acyl reductases involved in the biosynthesis of FA-derived pheromones. Research on insect pheromone-biosynthetic enzymes, which exhibit diverse enzymatic properties, has the potential to broaden the understanding of enzyme specificity determinants and contribute to engineering of enzymes with desired properties for biotechnological production of FA derivatives. Additionally, the application of such pheromone-biosynthetic enzymes represents an environmentally friendly and economic alternative to the chemical synthesis of pheromones that are used in insect pest management strategies.

Key genes of the sex pheromone biosynthesis pathway in female moths are required for pheromone quality and possibly mediate olfactory plasticity in conspecific male moths in Spodoptera litura

Ninety sex pheromone biosynthesis genes in Spodoptera litura were identified in transcriptome data and were investigated and classified into the following five groups: fatty acid synthase, beta oxidase, fatty acyl-coenzyme A (CoA) reductase, desaturase and acetyl-CoA acetyltransferase. Fourteen female-specific genes were identified through semi-quantitative PCR, and 15 additional genes had expression levels that were 3- to 10-fold higher in females than in males. The majority of the genes had higher expression levels in the sex pheromone glands. Injection of double-stranded RNA (dsRNA) against nine selected genes showed that down-regulation of Desaturase 3 (Des3), Des5 or fatty acyl coenzyme A reductase 17 (FAR17) significantly changed the ratio of the four sex pheromone components (Z,E) -9,11-tetradecadienyl acetate (Z9E11-14:Ac), (Z,E)-9,12-Tetradecadienyl acetate(Z9E12-14:Ac), (Z)-9-tetradecenyl acetate (Z9-14:Ac), (E)-11-Tetradecenyl acetate(E11-14:Ac). These key genes were differentially expressed in female moths collected from different geographical regions. Furthermore, field bioassays demonstrated geographical variation in the olfactory profile of male moths in response to the different sex pheromone mixtures, which therefore indicates that a significant variation in the sex pheromone components exists in the natural population. Our results suggest that a change in the expression of these key genes, Des3, Des5 and FAR17, in the sex pheromone biosynthesis pathway could change the ratio of the sex pheromone components. We surmise that the differential expression levels of the key genes of the sex pheromone biosynthesis pathway may lead to differential ratios of the sex pheromones in the field. Our field trapping experiment suggested that the change of the ratio of the sex pheromone components may have been adapted by the olfactory system and possibly mediate olfactory plasticity in conspecific male moths.

Transcriptome analysis in the beet webworm, Spoladea recurvalis (Lepidoptera: Crambidae)

The beet webworm, Spoladea recurvalis Fabricius, is a destructive pest on vegetable crops in tropics and subtropics; its main host plant is amaranth. It has become imperative to develop non-chemical methods to control S. recurvalis on amaranth. However, the lack of molecular information about this species has hindered the development of novel pest management strategies. In this study, high-throughput RNA sequencing covering de novo sequence assemblies, functional annotation of transcripts, gene function classification and enrichment was performed on S. recurvalis. Illumina sequencing generated a total of 120 435 transcript contigs ranging from 201 to 22 729 bases with a mean length of 688 bases. The assembled transcripts were subjected to Basic Local Alignment Search Tool-X (BLASTX) to obtain the annotations against non-redundant, Swiss-Prot, Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) protein databases. A subset of 58 225 transcript sequences returned hits from known proteins in the National Center for Biotechnology Information database, and the majority of the transcript sequences had the highest number of hits for Danaus plexippus (50.43%). A total of 1217 Gene Ontology-level 3 annotations were assigned to 51 805 transcripts, while 39 650 transcripts were predicted as functional protein-coding genes in the COG database and 20 037 transcripts were enriched to KEGG pathways. We identified 40 putative genes related to pheromone production and reception in S. recurvalis, with the expression of one gene between 0.29 and 1141.79 fragments per kilo base per million (FPKM) reads. The transcriptome sequence of S. recurvalis is a first step toward offering a comprehensive genomic resource which would enable better understanding of molecular mechanisms to enable development of effective pest management practices for this species.

Differential gene expression in the evolution of sex pheromone communication in New Zealand's endemic leafroller moths of the genera Ctenopseustis and Planotortrix

Background

Sex pheromone communication in moths has attracted the attention of evolutionary biologists due to the vast array of pheromone compounds used, addressing questions of how this diversity arose and how male reception has evolved in step with the female signal. Here we examine the role of changing gene expression in the evolution of mate recognition systems in leafroller moths, particularly focusing on genes involved in the biosynthetic pathways of sex pheromones in female pheromone glands and the peripheral reception repertoire in the antennae of males. From tissue-specific transcriptomes we mined and compared a database of genes expressed in the pheromone glands and antennae of males and females of four closely related species of leafroller moths endemic to New Zealand, Ctenopseutis herana and C. obliquana, and Planotortrix excessana and P. octo. The peculiarity of this group, compared to other Lepidoptera, is the use of (Z)-5-tetradecenyl acetate, (Z)-7-tetradecenyl acetate, and (Z)-8-tetradecenyl acetate as sex pheromone components.

Results

We identify orthologues of candidate genes from the pheromone biosynthesis pathway, degradation and transport, as well as genes of the periphery olfactory repertoire, including large families of binding proteins, receptors and odorant degrading enzymes. The production of distinct pheromone blends in the sibling species is associated with the differential expression of two desaturase genes, deast5 and desat7, in the pheromone glands. In male antennae, three odorant receptors, OR74, OR76a and OR30 are over-expressed, but their expression could not be clearly associated with the detection of species-specific pheromones components. In addition these species contain duplications of all three pheromone binding proteins (PBPs) that are also differentially expressed among species.

Conclusions

While in females differences in the expression of desaturases may be sufficient to explain pheromone blend differences among these New Zealand leafroller species, in males differential expression of several genes, including pheromone binding proteins, may underpin differences in the response by males to changing pheromone components among the species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4451-1) contains supplementary material, which is available to authorized users.

Identification of the pheromone biosynthesis genes from the sex pheromone gland transcriptome of the diamondback moth, Plutella xylostella

The diamondback moth was estimated to increase costs to the global agricultural economy as the global area increase of Brassica vegetable crops and oilseed rape. Sex pheromones traps are outstanding tools available in Integrated Pest Management for many years and provides an effective approach for DBM population monitoring and control. The ratio of two major sex pheromone compounds shows geographical variations. However, the limitation of our information in the DBM pheromone biosynthesis dampens our understanding of the ratio diversity of pheromone compounds. Here, we constructed a transcriptomic library from the DBM pheromone gland and identified genes putatively involved in the fatty acid biosynthesis, pheromones functional group transfer, and β-oxidation enzymes. In addition, odorant binding protein, chemosensory protein and pheromone binding protein genes encoded in the pheromone gland transcriptome, suggest that female DBM moths may receive odors or pheromone compounds via their pheromone gland and ovipositor system. Tissue expression profiles further revealed that two ALR, three DES and one FAR5 genes were pheromone gland tissue biased, while some chemoreception genes expressed extensively in PG, pupa, antenna and legs tissues. Finally, the candidate genes from large-scale transcriptome information may be useful for characterizing a presumed biosynthetic pathway of the DBM sex pheromone.

Expressional divergences of two desaturase genes determine the opposite ratios of two sex pheromone components in Helicoverpa armigera and Helicoverpa assulta

The sympatric closely related species Helicoverpa armigera and Helicoverpa assulta use 97:3 and 7:93 of (Z)-11-hexadecenal and (Z)-9-hexadecenal, respectively, as their sex pheromone to find/locate correct sex mates. Moreover, (Z)-11-hexadecenyl alcohol and (Z)-9-hexadecenyl alcohol are more abundant in the pheromone gland of H. assulta than in that of H. armigera. To clarify the molecular basis of these differences, we sequenced the pheromone gland transcriptomes of the two species and compared the expression patterns of the candidate enzyme genes involved in the pheromone biosynthetic pathways by FPKM values and quantitative RT-PCR analysis. We found that the desaturase gene LPAQ expressed about 70 times higher in H. armigera than in H. assulta, whereas another desaturase gene NPVE expressed about 60 times higher in H. assulta than in H. armigera. We also observed significantly higher expression of the fatty acyl reductase (FAR) gene FAR1 and the aldehyde reductase (AR) gene AR3 in H. assulta than in H. armigera. Examination of the pheromone glands of the backcross offspring of their hybrids to H. assulta showed a positive linear correlation between the expression level of LPAQ and the amount of Z11-16:Ald and between the expression level of NPVE and the amount of Z9-16:Ald in the pheromone glands. Taken together, these data demonstrate that the expressional divergences of LPAQ and NPVE determine the opposite sex pheromone component ratios in the two species and the divergent expression of FAR1 and AR3 may account for the greater accumulation of alcohols in the pheromone gland of H. assulta.

Identification and Expression Profiling of Chemosensory Genes in Dendrolimus punctatus Walker

Dendrolimus punctatus Walker is a serious pest affecting conifers in southern China. As extensive pesticide spraying is currently required to control D. punctatus, new control strategies are urgently needed. Chemosensory genes represent potential molecular targets for development of alternative pest control strategies, and the expression characteristics of these genes provide an indication of their function. To date, little information is available regarding chemosensory genes in D. punctatus or their expression profiles at different development stages and in various tissues. Here, we assembled and analyzed the transcriptomes of D. punctatus collected at different developmental stages and in a range of organs, using next-generation sequencing. A total of 171 putative chemosensory genes were identified, encoding 53 odorant binding proteins, 26 chemosensory proteins, 60 odorant receptors (OR), 12 gustatory receptors (GR), 18 ionotropic receptors (IR), and 2 sensory neuron membrane proteins (SNMPs). Expression analysis indicated that the antennae possess the largest number of highly expressed olfactory genes and that olfactory gene expression patterns in the eggs, larvae, and head were similar to one another, with each having moderate numbers of highly expressed olfactory genes. Fat body, ovary, midgut, and testis tissues also had similar olfactory gene expression patterns, including few highly expressed olfactory genes. Of particular note, we identified only two pheromone binding proteins and no pheromone receptors in D. punctatus, similar to our previous findings in Dendrolimus houi and Dendrolimus kikuchii, suggesting that insects of the Dendrolimus genus have different pheromone recognition characteristics to other Lepidopteran insects. Overall, this extensive expression profile analysis provides a clear map of D. punctatus chemosensory genes, and will facilitate functional studies and the development of new pest control methods in the future.

Antennal transcriptome analysis of the piercing moth Oraesia emarginata (Lepidoptera: Noctuidae)

The piercing fruit moth Oraesia emarginata is an economically significant pest; however, our understanding of its olfactory mechanisms in infestation is limited. The present study conducted antennal transcriptome analysis of olfactory genes using real-time quantitative reverse transcription PCR analysis (RT-qPCR). We identified a total of 104 candidate chemosensory genes from several gene families, including 35 olfactory receptors (ORs), 41 odorant-binding proteins, 20 chemosensory proteins, 6 ionotropic receptors, and 2 sensory neuron membrane proteins. Seven candidate pheromone receptors (PRs) and 3 candidate pheromone-binding proteins (PBPs) for sex pheromone recognition were found. OemaOR29 and OemaPBP1 had the highest fragments per kb per million fragments (FPKM) values in all ORs and OBPs, respectively. Eighteen olfactory genes were upregulated in females, including 5 candidate PRs, and 20 olfactory genes were upregulated in males, including 2 candidate PRs (OemaOR29 and 4) and 2 PBPs (OemaPBP1 and 3). These genes may have roles in mediating sex-specific behaviors. Most candidate olfactory genes of sex pheromone recognition (except OemaOR29 and OemaPBP3) in O. emarginata were not clustered with those of studied noctuid species (type I pheromone). In addition, OemaOR29 was belonged to cluster PRIII, which comprise proteins that recognize type II pheromones instead of type I pheromones. The structure and function of olfactory genes that encode sex pheromones in O. emarginata might thus differ from those of other studied noctuids. The findings of the present study may help explain the molecular mechanism underlying olfaction and the evolution of olfactory genes encoding sex pheromones in O. emarginata.

Molecular identification of differential expression genes associated with sex pheromone biosynthesis in Spodoptera exigua

Species-specific sex pheromone is biosynthesized and released in most female moths as a chemical cue in mating communication. However, information on genes involved in this pathway is limited. The beet armyworm, Spodoptera exigua, is a cosmopolitan agricultural pest that causes severe economic losses to many crops. In China, the female sex pheromones in sex pheromone glands (PGs) of S. exigua have been measured which comprises (Z,E)-9,12-tetradecadienyl acetate, (Z)-9-tetradecen-l-ol, (Z)-9-tetradecenyl acetate, and (Z,E)-9,12-tetradecadien-1-ol in a ratio of 47:18:18:17. Fifty-nine putative genes related to sex pheromone biosynthesis were identified in the present study by sequencing and analyzing the sex pheromone gland (PG) transcriptome of S. exigua. Expression profiles revealed that two desaturase (SexiDes5 and SexiDes11) and three fatty acyl reductase (SexiFAR2, 3, and 9) genes had PG-specific expression, and phylogenetic analysis demonstrated that they clustered with genes known to be involved in pheromone synthesis in other moth species. Our results provide crucial background information that could facilitate the elucidation of sex pheromone biosynthesis pathway of S. exigua as well as other Spodoptera species and help identify potential targets for disrupting sexual communication in S. exigua for developing novel environment-friendly pesticides.

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

Background

Female moths synthesize species-specific sex pheromone components and release them to attract male moths, which depend on precise sex pheromone chemosensory system to locate females. Two types of genes involved in the sex pheromone biosynthesis and degradation pathways play essential roles in this important moth behavior. To understand the function of genes in the sex pheromone pathway, this study investigated the genome-wide and digital gene expression of sex pheromone biosynthesis and degradation genes in various adult tissues in the diamondback moth (DBM), Plutella xylostella, which is a notorious vegetable pest worldwide.

Results

A massive transcriptome data (at least 39.04 Gb) was generated by sequencing 6 adult tissues including male antennae, female antennae, heads, legs, abdomen and female pheromone glands from DBM by using Illumina 4000 next-generation sequencing and mapping to a published DBM genome. Bioinformatics analysis yielded a total of 89,332 unigenes among which 87 transcripts were putatively related to seven gene families in the sex pheromone biosynthesis pathway. Among these, seven [two desaturases (DES), three fatty acyl-CoA reductases (FAR) one acetyltransferase (ACT) and one alcohol dehydrogenase (AD)] were mainly expressed in the pheromone glands with likely function in the three essential sex pheromone biosynthesis steps: desaturation, reduction, and esterification. We also identified 210 odorant-degradation related genes (including sex pheromone-degradation related genes) from seven major enzyme groups. Among these genes, 100 genes are new identified and two aldehyde oxidases (AOXs), one aldehyde dehydrogenase (ALDH), five carboxyl/cholinesterases (CCEs), five UDP-glycosyltransferases (UGTs), eight cytochrome P450 (CYP) and three glutathione S-transferases (GSTs) displayed more robust expression in the antennae, and thus are proposed to participate in the degradation of sex pheromone components and plant volatiles.

Conclusions

To date, this is the most comprehensive gene data set of sex pheromone biosynthesis and degradation enzyme related genes in DBM created by genome- and transcriptome-wide identification, characterization and expression profiling. Our findings provide a basis to better understand the function of genes with tissue enriched expression. The results also provide information on the genes involved in sex pheromone biosynthesis and degradation, and may be useful to identify potential gene targets for pest control strategies by disrupting the insect-insect communication using pheromone-based behavioral antagonists.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3592-y) contains supplementary material, which is available to authorized users.

Ultrastructure of sensilla of antennae and ovipositor of Sitotroga cerealella (Lepidoptera: Gelechiidae), and location of female sex pheromone gland

The Angoumois grain moth, Sitotroga cerealella, is a serious pest of stored grains worldwide. Presently, the best effective control against the moth is to disrupt the sexual communication between sexes. Sexual communication in moths includes two processes in which females produce and release pheromones from the sex pheromone gland and males detect and respond to them with a relatively sophisticated olfactory system in their antennae. To better understand these processes, we studied the ultrastructure of antennal and ovipositor sensilla of S. cerealella and determined the location of the female sex pheromone gland. Seven types of antennal sensilla were identified on both sexes: sensilla trichodea, sensilla chaetica, sensilla coeloconica, sensilla styloconica, sensilla auricillica, sensilla squamiformia and Bӧhm bristles. Of these sensilla, the sensilla trichodea were significantly more abundant on male antennae than on those of females, suggesting that these sensilla may detect the sex pheromones. On the ovipositor, only sensilla chaetica of various lengths were found. The sexual gland was an eversible sac of glandular epithelium that was situated dorsally in the intersegmental membrane between the 8^th and 9^th abdominal segments. These results will lead to a better understanding of mate finding with sex pheromones for this worldwide pest species.

Identification and Expression Analysis of Candidate Odorant-Binding Protein and Chemosensory Protein Genes by Antennal Transcriptome of Sitobion avenae

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) of aphids are thought to be responsible for the initial molecular interactions during olfaction that mediate detection of chemical signals. Analysis of the diversity of proteins involved comprises critical basic research work that will facilitate the development of sustainable pest control strategies. To help us better understand differences in the olfactory system between winged and wingless grain aphids, we constructed an antennal transcriptome from winged and wingless Sitobion avenae (Fabricius), one of the most serious pests of cereal fields worldwide. Among the 133,331 unigenes in the antennal assembly, 13 OBP and 5 CSP putative transcripts were identified with 6 OBP and 3 CSP sequences representing new S. avenae annotations. We used qPCR to examine the expression profile of these genes sets across S. avenae development and in various tissues. We found 7 SaveOBPs and 1 SaveCSP were specifically or significantly elevated in antennae compared with other tissues, and that some transcripts (SaveOBP8, SaveCSP2 and SaveCSP5) were abundantly expressed in the legs of winged or wingless aphids. The expression levels of the SaveOBPs and SaveCSPs varied depending on the developmental stage. Possible physiological functions of these genes are discussed. Further molecular and functional studies of these olfactory related genes will explore their potential as novel targets for controlling S. avenae.

Antennal transcriptome and differential expression of olfactory genes in the yellow peach moth, Conogethes punctiferalis (Lepidoptera: Crambidae)

The yellow peach moth (YPM), Conogethes punctiferalis (Guenée), is a multivoltine insect pest of crops and fruits. Antennal-expressed receptors are important for insects to detect olfactory cues for host finding, mate attraction and oviposition site selection. However, few olfactory related genes were reported in YPM until now. In the present study, we sequenced and characterized the antennal transcriptomes of male and female YPM. In total, 15 putative odorant binding proteins (OBPs), 46 putative odorant receptors (ORs) and 7 putative ionotropic receptors (IRs) were annotated and identified as olfactory-related genes of C. punctiferalis. Further analysis of RT-qPCR revealed that all these olfactory genes are primarily or uniquely expressed in male and female antennae. Among which, 3 OBPs (OBP4, OBP8 and PBP2) and 4 ORs (OR22, OR26, OR44 and OR46) were specially expressed in male antennae, whereas 4 ORs (OR5, OR16, OR25 and OR42) were primarily expressed in female antennae. The predicted protein sequences were compared with homologs in other lepidopteran species and model insects, which showed high sequence homologies between C. punctiferalis and O. furnacalis. Our work allows for further functional studies of pheromone and general odorant detection genes, which might be meaningful targets for pest management.

The Genetic Basis of Pheromone Evolution in Moths

Moth sexual pheromones are widely studied as a fine-tuned system of intraspecific sexual communication that reinforces interspecific reproductive isolation. However, their evolution poses a dilemma: How can the female pheromone and male preference simultaneously change to create a new pattern of species-specific attraction? Solving this puzzle requires us to identify the genes underlying intraspecific variation in signals and responses and to understand the evolutionary mechanisms responsible for their interspecific divergence. Candidate gene approaches and functional analyses have yielded insights into large families of biosynthetic enzymes and pheromone receptors, although the factors controlling their expression remain largely unexplored. Intra- and interspecific crosses have provided tantalizing evidence of regulatory genes, although, to date, mapping resolution has been insufficient to identify them. Recent advances in high-throughput genome and transcriptome sequencing, together with established techniques, have great potential to help scientists identify the specific genetic changes underlying divergence and resolve the mystery of how moth sexual communication systems evolve.

Identification and Expression Profiles of Sex Pheromone Biosynthesis and Transport Related Genes in Spodoptera litura

Although the general pathway of sex pheromone synthesis in moth species has been established, the molecular mechanisms remain poorly understood. The common cutworm Spodoptera litura is an important agricultural pest worldwide and causes huge economic losses annually. The female sex pheromone of S. litura comprises Z9,E11-14:OAc, Z9,E12-14:OAc, Z9-14:OAc, and E11-14:OAc. By sequencing and analyzing the transcriptomic data of the sex pheromone glands, we identified 94 candidate genes related to pheromone biosynthesis (55 genes) or chemoreception (39 genes). Gene expression patterns and phylogenetic analysis revealed that two desaturase genes (SlitDes5 and SlitDes11) and one fatty acyl reductase gene (SlitFAR3) showed pheromone gland (PG) biased or specific expression, and clustered with genes known to be involved in pheromone synthesis in other moth species. Furthermore, 4 chemoreception related genes (SlitOBP6, SlitOBP11, SlitCSP3, and SlitCSP14) also showed higher expression in the PG, and could be additional candidate genes involved in sex pheromone transport. This study provides the first solid background information that should facilitate further elucidation of sex pheromone biosynthesis and transport, and indicates potential targets to disrupt sexual communication in S. litura for a novel pest management strategy.

Diapause hormone in the Helicoverpa/Heliothis complex: A review of gene expression, peptide structure and activity, analog and antagonist development, and the receptor

This review summarizes recent studies focusing on diapause hormone (DH) in the Helicoverpa/Heliothis complex of agricultural pests. Moths in this complex overwinter in pupal diapause, a form of developmental arrest used to circumvent unfavorable seasons. DH was originally reported in the silkmoth Bombyx mori, a species that relies on DH to induce an embryonic diapause. But, in the case of Helicoverpa/Heliothis, levels of dh transcripts and DH peptides are more abundant in nondiapausing pupae than in diapausing individuals, and DH effectively terminates diapause within a specific temperature range. A structure activity relationship study indicated that the active core of DH is the C-terminal hepta-peptide, LWFGPRLa. We designed and synthesized a first generation of DH agonists and identified two agonists (PK-2Abf and PK-Etz) that were nearly 50- and 13-fold more potent than the native hormone. These studies revealed two structural characteristics of DH and its agonists that are essential for interaction with the receptor: a trans-Pro configuration to form a type I β-turn and a hydrophobic moiety involved in ligand binding. Modification of DH at the active core yielded a potent DH antagonist (DH-Jo, acetyl-GLWA[Jo]RLa) as well as an agonist (DH-2Abf-K). Three compounds (Decyl-1963, Dodecyl-1967, Heptyl-1965) were identified as agents capable of penetrating the cuticle of young pupae and thereby preventing the entry into diapause. DH receptor cDNA was cloned and an effective in vitro high throughput screen system was established for future use. This work sets the stage for further development of DH analogs and antagonists that have the potential to disrupt insect diapause as a tool for pest management.

Analysis of the Agrotis segetum pheromone gland transcriptome in the light of sex pheromone biosynthesis

BACKGROUND

Moths rely heavily on pheromone communication for mate finding. The pheromone components of most moths are modified from the products of normal fatty acid metabolism by a set of tissue-specific enzymes. The turnip moth, Agrotis segetum uses a series of homologous fatty-alcohol acetate esters ((Z)-5-decenyl, (Z)-7-dodecenyl, and (Z)-9 tetradecenyl acetate) as its sex pheromone components. The ratio of the components differs between populations, making this species an interesting subject for studies of the enzymes involved in the biosynthetic pathway and their influence on sex pheromone variation.

RESULTS

Illumina sequencing and comparative analysis of the transcriptomes of the pheromone gland and abdominal epidermal tissue, enabled us to identify genes coding for putative key enzymes involved in the pheromone biosynthetic pathway, such as fatty acid synthase, β-oxidation enzymes, fatty-acyl desaturases (FAD), fatty-acyl reductases (FAR), and acetyltransferases. We functionally assayed the previously identified ∆11-desaturase [GenBank: ES583599, JX679209] and FAR [GenBank: JX679210] and candidate acetyltransferases (34 genes) by heterologous expression in yeast. The functional assay confirmed that the ∆11-desaturase interacts with palmitate and produces (Z)-11-hexadecenoate, which is the common unsaturated precursor of three homologous pheromone component acetates produced by subsequent chain-shortening, reduction and acetylation. Much lower, but still visible, activity on 14C and 12C saturated acids may account for minor pheromone compounds previously observed in the pheromone gland. The FAR characterized can operate on various unsaturated fatty acids that are the immediate acyl precursors of the different A. segetum pheromone components. None of the putative acetyltransferases that we expressed heterologously did acetylate any of the fatty alcohols tested as substrates.

CONCLUSIONS

The massive sequencing technology generates enormous amounts of candidate genes potentially involved in pheromone biosynthesis but testing their function by heterologous expression or gene silencing is a bottleneck. We confirmed the function of a previously identified desaturase gene and a fatty-acyl reductase gene by heterologous expression, but the acetyltransferase postulated to be involved in pheromone biosynthesis remains illusive, in spite of 34 candidates being assayed. We also generated lists of gene candidates that may be useful for characterizing the acetyl-CoA carboxylase, fatty acid synthetase and β-oxidation enzymes.

The regulation of Δ11-desaturase gene expression in the pheromone gland of Mamestra brassicae (Lepidoptera; Noctuidae) during pheromonogenesis

Cabbage moth (Mamestra brassicae) females produce sex pheromones to attract conspecific males. In our M. brassicae colony, the pheromone blend is composed of Z11-hexadecenyl acetate (Z11-16Ac) and hexadecyl acetate (16Ac) in a 93:7 ratio. A fatty acyl Δ11-desaturase is involved in the production of the main pheromone component. The release of Pheromone Biosynthesis Activating Neuropeptide (PBAN) regulates the pheromone production in the pheromone gland (PG). We cloned a cDNA encoding the MambrΔ11-desaturase and analyzed its expression profile over time in M. brassicae tissues. Transcript levels of the Δ11-desaturase in larvae, pupal PGs, fat body, brain and muscle tissues were <0.1% of that in female PGs, whereas expression in male genitalia was 2%. In the PGs of virgin females the expression level increased continuously from eclosion to the end of the 1st day when it reached a plateau without further significant fluctuation up to the 8th day. In contrast, we recorded a characteristic daily rhythmicity in pheromone production with a maximum around 200 ng Z11-16Ac/PG. In some experiments, females were decapitated to prevent PBAN release and thereby inhibit pheromone production, which remarkably increased after treatment with Mambr-Pheromonotropin. Further experiments revealed that mating resulted in a significant suppression of pheromone production. However, expression of the Δ11-desaturase was not affected by any of these interventions, suggesting that it's not regulated by PBAN. Fluorescent microscopy was used to study the potential role of lipid droplets during pheromone production, however, no lipid droplets were identified indicating that pheromonogenesis is regulated via de novo fatty acid synthesis.

Genes involved in sex pheromone biosynthesis of Ephestia cautella, an important food storage pest, are determined by transcriptome sequencing

BACKGROUND

Insects use pheromones, chemical signals that underlie all animal behaviors, for communication and for attracting mates. Synthetic pheromones are widely used in pest control strategies because they are environmentally safe. The production of insect pheromones in transgenic plants, which could be more economical and effective in producing isomerically pure compounds, has recently been successfully demonstrated. This research requires information regarding the pheromone biosynthetic pathways and the characterization of pheromone biosynthetic enzymes (PBEs). We used Illumina sequencing to characterize the pheromone gland (PG) transcriptome of the Pyralid moth, Ephestia cautella, a destructive storage pest, to reveal putative candidate genes involved in pheromone biosynthesis, release, transport and degradation.

RESULTS

We isolated the E. cautella pheromone compound as (Z,E)-9,12-tetradecadienyl acetate, and the major pheromone precursors 16:acyl, 14:acyl, E14-16:acyl, E12-14:acyl and Z9,E12-14:acyl. Based on the abundance of precursors, two possible pheromone biosynthetic pathways are proposed. Both pathways initiate from C16:acyl-CoA, with one involving ∆14 and ∆9 desaturation to generate Z9,E12-14:acyl, and the other involving the chain shortening of C16:acyl-CoA to C14:acyl-CoA, followed by ∆12 and ∆9 desaturation to generate Z9,E12-14:acyl-CoA. Then, a final reduction and acetylation generates Z9,E12-14:OAc. Illumina sequencing yielded 83,792 transcripts, and we obtained a PG transcriptome of ~49.5 Mb. A total of 191 PBE transcripts, which included pheromone biosynthesis activating neuropeptides, fatty acid transport proteins, acetyl-CoA carboxylases, fatty acid synthases, desaturases, β-oxidation enzymes, fatty acyl-CoA reductases (FARs) and fatty acetyltransferases (FATs), were selected from the dataset. A comparison of the E. cautella transcriptome data with three other Lepidoptera PG datasets revealed that 45% of the sequences were shared. Phylogenetic trees were constructed for desaturases, FARs and FATs, and transcripts that clustered with the ∆14, ∆12 and ∆9 desaturases, PG-specific FARs and potential candidate FATs, respectively, were identified. Transcripts encoding putative pheromone degrading enzymes, and candidate pheromone carrier and receptor proteins expressed in the E. cautella PG, were also identified.

CONCLUSIONS

Our study provides important background information on the enzymes involved in pheromone biosynthesis. This information will be useful for the in vitro production of E. cautella sex pheromones and may provide potential targets for disrupting the pheromone-based communication system of E. cautella to prevent infestations.

Identification and Characterization of Candidate Chemosensory Gene Families from Spodoptera exigua Developmental Transcriptomes

Insect chemosensory genes have been considered as potential molecular targets to develop alternative strategies for pest control. However, in Spodoptera exigua, a seriously polyphagous agricultural pest, only a small part of such genes have been identified and characterized to date. Here, using a bioinformatics screen a total of 79 chemosensory genes were identified from a public transcriptomic data of different developmental stages (eggs, 1st to 5th instar larvae, pupae, female and male adults), including 34 odorant binding proteins (OBPs), 20 chemosensory proteins (CSPs), 22 chemosensory receptors (10 odorant receptors (ORs), six gustatory receptors (GRs) and six ionotropic receptors (IRs)) and three sensory neuron membrane proteins (SNMPs). Notably, a new group of lepidopteran SNMPs (SNMP3 group) was found for the first time in S. exigua, and confirmed in four other moth species. Further, reverse transcription PCR (RT-PCR) and quantitative real time PCR (qPCR) were employed respectively to validate the sequences and determine the expression patterns of 69 identified chemosensory genes regarding to sexes, tissues and stages. Results showed that 67 of these genes could be detected and reconstructed in at least one tissue tested. Further, 60 chemosensory genes were expressed in adult antennae and 52 in larval heads with the antennae, whereas over half of the genes were also detected in non-olfactory tissues like egg and thorax. Particularly, S. exigua OBP2 showed a predominantly larval head-biased expression, and functional studies further indicated its potentially olfactory roles in guiding food searching of larvae. This work suggests functional diversities of S. exigua chemosensory genes and could greatly facilitate the understanding of olfactory system in S. exigua and other lepidopteran species.

Generation of a Transcriptome in a Model Lepidopteran Pest, Heliothis virescens, Using Multiple Sequencing Strategies for Profiling Midgut Gene Expression

Heliothine pests such as the tobacco budworm, Heliothis virescens (F.), pose a significant threat to production of a variety of crops and ornamental plants and are models for developmental and physiological studies. The efforts to develop new control measures for H. virescens, as well as its use as a relevant biological model, are hampered by a lack of molecular resources. The present work demonstrates the utility of next-generation sequencing technologies for rapid molecular resource generation from this species for which lacks a sequenced genome. In order to amass a de novo transcriptome for this moth, transcript sequences generated from Illumina, Roche 454, and Sanger sequencing platforms were merged into a single de novo transcriptome assembly. This pooling strategy allowed a thorough sampling of transcripts produced under diverse environmental conditions, developmental stages, tissues, and infections with entomopathogens used for biological control, to provide the most complete transcriptome to date for this species. Over 138 million reads from the three platforms were assembled into the final set of 63,648 contigs. Of these, 29,978 had significant BLAST scores indicating orthologous relationships to transcripts of other insect species, with the top-hit species being the monarch butterfly (Danaus plexippus) and silkworm (Bombyx mori). Among identified H. virescens orthologs were immune effectors, signal transduction pathways, olfactory receptors, hormone biosynthetic pathways, peptide hormones and their receptors, digestive enzymes, and insecticide resistance enzymes. As an example, we demonstrate the utility of this transcriptomic resource to study gene expression profiling of larval midguts and detect transcripts of putative Bacillus thuringiensis (Bt) Cry toxin receptors. The substantial molecular resources described in this study will facilitate development of H. virescens as a relevant biological model for functional genomics and for new biological experimentation needed to develop efficient control efforts for this and related Noctuid pest moths.

Transcriptome and expression profiling analysis link patterns of gene expression to antennal responses in Spodoptera litura

Background

The study of olfaction is key to understanding the interaction of insects with their environment and provides opportunities to develop novel tactics for control of pest species. Recent developments in transcriptomic approaches enable the molecular basis of olfaction to be studied even in species with limited genomic information. Here we use transcriptome and expression profiling analysis to characterize the antennal transcriptome of the noctuid moth and polyphagous pest Spodoptera litura.

Results

We identify 74 candidate genes involved in odor detection and recognition, encoding 26 ORs, 21 OBPs, 18 CSPs and 9 IRs. We examine their expression levels in both sexes and seek evidence for their function by relating their expression with levels of EAG response in male and female antennae to 58 host and non-host plant volatiles and sex pheromone components. The majority of olfactory genes showed sex-biased expression, usually male-biased in ORs. A link between OR gene expression and antennal responses to odors was evident, a third of the compounds tested evoking a sex-biased response, in every case also male-biased. Two candidate pheromone receptors, OR14 and OR23 were especially strongly expressed and male-biased and we suggest that these may respond to the two female sex pheromone components of S. litura, Z9E11-14:OAc and Z9E12-14:OAc, which evoked strongly male-biased EAG responses.

Conclusions

Our results provide the molecular basis for elucidating the olfactory profile of moths and the sexual divergence of their behavior and could enable the targeting of particular genes, and behaviors for pest management.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1375-x) contains supplementary material, which is available to authorized users.