Transcriptome comparison of the sex pheromone glands from two sibling Helicoverpa species with opposite sex pheromone components

Knockout of the delta11-desaturase SfruDES1 disrupts sex pheromone biosynthesis, mating and oviposition in the fall armyworm, Spodoptera frugiperda

Moth insects rely on sex pheromones for long distance attraction and searching for sex partners. The biosynthesis of moth sex pheromones involves the catalytic action of multiple enzymes, with desaturases playing a crucial role in the process of carbon chain desaturation. However, the specific desaturases involved in sex pheromone biosynthesis in fall armyworm (FAW), Spodoptera frugiperda, have not been clarified. In this study, a Δ11 desaturase (SfruDES1) gene in FAW was knocked out using the CRISPR/Cas9 genome editing system. A homozygous mutant of SfruDES1 was obtained through genetic crosses. The gas chromatography-mass spectrometry (GC-MS) analysis results showed that the three main sex pheromone components (Z7-12:Ac, Z9-14:Ac, and Z11-16:Ac) and the three minor components (Z9-14:Ald, E11-14:Ac and Z11-14:Ac) of FAW were not detected in homozygous mutant females compared to the wild type. Furthermore, behavioral assay demonstrated that the loss of SfruDES1 resulted in a significant reduction in the attractiveness of females to males, along with disruptions in mating behavior and oviposition. Additionally, in a heterologous expression system, recombinant SfruDES1 could introduce a cis double bond at the Δ11 position in palmitic acid, which resulted in the changes in components of the synthesized products. These findings suggest desaturase plays a key role in the biosynthesis of sex pheromones, and knockout of the SfruDES1 disrupts sex pheromone biosynthesis and mating behavior in FAW. The SfruDES1 could serve as tool to develop a control method for S. frugiperda.

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 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.

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 a unique three-component sex pheromone produced by the tea black tussock moth, Dasychira baibarana (Lepidoptera: Erebidae: Lymantriinae)

BACKGROUND

The tea black tussock moth Dasychira baibarana Matsumura is a devastating pest in tea plantations that causes substantial economic losses. Presently, there is no effective method to control this pest other than pesticide application. The identified sex pheromone of D. baibarana could be used for detecting and monitoring this pest.

RESULTS

Gas chromatography-electroantennogram detection showed that D. baibarana male moth antennae responded strongly to three components in the female sex pheromone gland, which were identified as: (3Z,6Z)-cis-9,10-epoxyhenicosa-3,6-diene (Z3,Z6,epo9-21:H) (I), (3Z,6Z,11E)-cis-9,10-epoxyhenicosa-3,6,11-triene (Z3,Z6,epo9,E11-21:H) (II) and (3Z,6Z)-henicosa-3,6-dien-11-one (Z3,Z6-21:11-one) (III). Dasychira baibarana uses a unique composition of an epoxydiene, epoxytriene, and dienone with the same 3,6-dienyl motif as its sex pheromone. The epoxytriene and dienone were not previously characterized as insect pheromone components. Electroantennogram analysis showed that each synthetic compound strongly stimulated male antennae, and compounds II and III elicited stronger responses than compound I. A wind tunnel bioassay and field trapping experiments proved that, ternary blends of compounds I-III attracted D. baibarana. Efficient attraction was achieved with a rubber septum baited with 500 μg of a mixture of compounds I-III at the ratio 25:20:55.

CONCLUSION

The three identified compounds elicited an electroantennogram response in D. baibarana male moth antennae, and a mixture of the three components at the ratio 25:20:55 attracted D. baibarana male moths in a wind tunnel assay and field trapping experiments. This blend could be used for integrated management of D. baibarana in tea plantations. © 2022 Society of Chemical Industry.

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.

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.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

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.

The whole body transcriptome of Coleophora obducta reveals important olfactory proteins

Background

The tiny casebearer moth Coleophora obducta, an important defoliator of Larix spp., is a major threat to ecological security in north China. Studies have shown that C. obducta is strongly specific to host plants; it is unable complete its life cycle without Larix spp. The sex pheromones of C. obducta Z5-10:OH have been elucidated; and eight types of antennae sensilla, have been detected, indicating that an exploration of its olfactory proteins is necessary, due to the general lack of information on this topic.

Methods

We investigated the whole body transcriptome of C. obducta, performed a phylogenetic analysis of its olfactory proteins and produced expression profiles of three pheromone-binding proteins (CobdPBPs) by qRT–PCR.

Results

We identified 16 odorant binding proteins, 14 chemosensory proteins, three sensory neuron membrane proteins, six odorant degrading enzymes, five antennal esterases, 13 odorant receptors, seven ionotropic receptors and 10 gustatory receptors, including three PBPs and one odorant co-receptor. Additionally, three putative pheromone receptors, two bitter gustatory receptors and five functional ionotropic receptors were found by phylogenetic analysis. The expression profiles of three PBPs in males and females showed that all of them exhibited male-specific expression and two were expressed at significantly higher levels in males. These data provide a molecular foundation from which to explore the olfactory recognition process and may be useful in the development of a new integrated pest management strategy targeting olfactory recognition of C. obducta.

Functional identification of fatty acyl reductases in female pheromone gland and tarsi of the corn earworm, Helicoverpa zea

Most moths utilize sex pheromones released by the female to attract a mate. Females produce the sex pheromone in the pheromone gland in a biosynthetic pathway which consists of several key enzymes. Fatty acyl-CoA reductase is one of the key enzymes, which catalyzes the conversion of fatty acyl-CoA to the corresponding alcohol, playing an important role in producing the final proportion of each pheromone component. In Helicoverpa zea, (Z)-11-hexadecenal is the major sex pheromone component in female pheromone glands and previously a large amount of hexadecanal was also found in female and male tarsi. In our previous study, we compared the transcriptome between pheromone glands and tarsi and found 20 fatty acyl-CoA reductases in both tissues. In this study, we functionally characterized four FARs which were expressed at high levels according to the transcriptome of pheromone glands and tarsi. Fatty acyl-CoA reductase 1 was homologous to other moth pheromone gland specific fatty acyl-CoA reductases, and it was also present in male tarsi. Functional expression in yeast cells indicates that only fatty acyl-CoA reductase 1 was able to produce fatty alcohols. In addition, a decreased mRNA level of fatty acyl-CoA reductase 1 in female pheromone glands and male tarsi by RNAi knockdown caused a significant decrease in the production of (Z)-11-hexadecenal in pheromone glands and hexadecanal in male tarsi. This study is the first to demonstrate the direct function of a fatty acyl-CoA reductase in male tarsi and also confirms its role in sex pheromone biosynthesis in H. zea.

Candidates for chemosensory genes identified in the Chinese citrus fly, Bactrocera minax, through a transcriptomic analysis

Background

The males of many Bactrocera species (Diptera: Tephritidae) respond strongly to plant-derived chemicals (male lures) and can be divided into cue lure/raspberry ketone (CL/RK) responders, methyl eugenol (ME) responders and non-responders. Representing a non-responders, Bactrocera minax display unique olfactory sensory characteristics compared with other Bactrocera species. The chemical senses of insects mediate behaviors that are associated with survival and reproduction. Here, we report the generation of transcriptomes from antennae and the rectal glands of both male and female adults of B. minax using Illumina sequencing technology, and annotated gene families potentially responsible for chemosensory.

Results

We developed four transcriptomes from different tissues of B. minax and identified a set of candidate genes potentially responsible for chemosensory by analyzing the transcriptomic data. The candidates included 40 unigenes coding for odorant receptors (ORs), 30 for ionotropic receptors (IRs), 17 for gustatory receptors (GRs), three for sensory neuron membrane proteins (SNMPs), 33 for odorant-binding proteins (OBPs), four for chemosensory proteins (CSPs). Sex- and tissue-specific expression profiles for candidate chemosensory genes were analyzed via transcriptomic data analyses, and expression profiles of all ORs and antennal IRs were investigated by real-time quantitative PCR (RT-qPCR). Phylogenetic analyses were also conducted on gene families and paralogs from other insect species together.

Conclusions

A large number of chemosensory genes were identified from transcriptomic data. Identification of these candidate genes and their expression profiles in various tissues provide useful information for future studies towards revealing their function in B. minax.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6022-5) contains supplementary material, which is available to authorized users.

Pheromone gland transcriptome of the pink bollworm moth, Pectinophora gossypiella: Comparison between a laboratory and field population

The pink bollworm, Pectinophora gossypiella, is a world-wide pest of cotton and in some parts of the cotton growing region is controlled by the mating disruption technique using synthetic sex pheromone. The sex pheromone consists of two compounds, (Z,Z)- and (Z,E)-7,11-hexadecadienyl acetates, in about a 50:50 ratio. However, recently, a population with sex pheromone compound ratios of about 62:38 were found in cotton fields that use mating disruption in Israel. To investigate how the change developed, we compared the pheromone gland transcriptomes between a reference laboratory population and a population obtained from an Israeli cotton field utilizing mating disruption. We analyzed four biological replicates from each population and found transcripts encoding 17 desaturases, 8 reductases, and 17 candidate acetyltransferases in both populations, which could be involved in sex pheromone biosynthesis. The expression abundance of some genes between the two populations was different. Some desaturases and candidate acetyltransferases were found to have mutated in one of the populations. The differentially expressed genes play potential roles in sex pheromone biosynthesis and could be involved in causing altered female sex pheromone ratios in the field population.

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.

Identification of Sex Pheromone of Miltochrista striata (Lepidoptera: Arctiidae)

A species of lichen moth, Miltochrista striata (Bremer & Grey, 1852), feeds on lichens in Chinese tea plantations (Camellia sinensis (L.) O. Kuntze (Ericales:Theaceae)). A previous sex attractant screening test showed that male moths of M. striata were attracted by a mixture of (Z,Z,Z)-3,6,9-octadecatriene (Z3,Z6,Z9-18:H), (Z,Z,Z)-3,6,9-nonadecatriene (Z3,Z6,Z9-19:H), and their monoepoxy derivatives. To determine which of the component is an effective sex attractant for M. striata, the sex pheromone glands of female moths were excised and extracted with n-hexane. By comparison with the retention time and mass spectra of synthetic chemicals, two compounds in the crude extracts were identified as Z3,Z6,Z9-18:H and (Z,Z)-3,6-cis-9,10-epoxy-octadecadiene (Z3,Z6,epo9-18:H) using gas chromatography-mass spectrometry. The results of electroantennographic tests showed that the electrophysiological activities of Z3,Z6,Z9-18:H and Z3,Z6,epo9-18:H were distinctly higher than those of (Z,Z)-6,9-cis-3,4-epoxy-octadecadiene, (Z,Z)-3,9-cis-6,7-epoxy-octadecadiene, (Z,Z,Z)-3,6,9-nonadecatriene, (Z,Z)-6,9-cis-3,4-epoxy-nonadecatriene, (Z,Z)-3,9-cis-6,7-epoxy-nonadecatriene, and (Z,Z)-3,6-cis-9,10-epoxy-nonadecatriene. Field trapping showed that only a mixture of Z3,Z6,Z9-18:H and Z3,Z6,epo9-18:H attracted male moths, and the optimal mixture of these compounds was the ratio of 4:6 at 1.0-mg dosage. The results represent the first determination of the sex pheromone of a lichen moth in a tea plantation and provide a scientific basis to develop an effective protocol using sex pheromone to monitor populations of M. striata.

Evidence of Premating Isolation Between Two Sibling Moths: Ectropis grisescens and Ectropis obliqua (Lepidoptera: Geometridae)

The sex pheromones of Ectropis grisescens Warren and Ectropis obliqua Prout were both reported to contain (Z,Z,Z)-3,6,9-octadecatriene (Z3,Z6,Z9-18:H) and (Z,Z)-3,9-cis-6,7-epoxy-octadecadiene (Z3,epo6,Z9-18:H). To clarify how these two sibling geometrids maintain premating isolation, the female sex pheromones of the two species were reexamined. Gas chromatography-electroantennographic detection (GC-EAD) and gas chromatography-mass spectrometry revealed two GC-EAD-active compounds, Z3,Z6,Z9-18:H and Z3,epo6,Z9-18:H, in E. grisescens female pheromone glands as well as an additional GC-EAD-active compound, (Z,Z)-3,9-cis-6,7-epoxy-nonadecadiene (Z3,epo6,Z9-19:H), in E. obliqua female pheromone glands. Synthesized Z3,Z6,Z9-18:H and Z3,epo6,Z9-18:H elicited dose-dependent electroantennogram (EAG) responses from male antennae of both E. grisescens and E. obliqua. However, Z3,epo6,Z9-19:H only elicited dose-dependent EAG responses from E. obliqua and limited EAG responses from E. grisescens at all doses. In wind-tunnel studies, lures that contained Z3,Z6,Z9-18:H and Z3,epo6,Z9-18:H attracted E. grisescens males and had no effect on E. obliqua males. The addition of Z3,epo6,Z9-19:H to the blend of Z3,Z6,Z9-18:H and Z3,epo6,Z9-18:H strongly attracted E. obliqua males but had a limited attraction for E. grisescens males. Thus, Z3,Z6,Z9-18:H and Z3,epo6,Z9-18:H were sex pheromone components of E. grisescens, whereas Z3,Z6,Z9-18:H, Z3,epo6,Z9-18:H and Z3,epo6,Z9-19:H were sex pheromone components of E. obliqua. The presence or absence of Z3,epo6,Z9-19:H played a central role in the premating isolation of these two sibling species.

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.

Transcriptome analysis of Helicoverpa armigera male hairpencils: Alcohol biosynthesis and requirement for mating success

Many female animals use different strategies to assess male quality to increase their own reproductive fitness. In moths, females usually use chemical signals (sex pheromones) to attract males from a distance. Once males approach a female, they release close range pheromones from hairpencils to facilitate female acceptance. However, detailed mechanisms involved in male sex pheromone biosynthesis and its action in promoting female acceptance have not yet been fully characterized. This study screened a series of candidate genes via a transcriptome analysis of the male hairpencil of Helicoverpa armigera. Using pharmacological inhibitor and RNAi-mediated knockdown assays, we demonstrated that Ca²⁺ and cyclic-AMP were involved in pheromone biosynthesis activating neuropeptide (PBAN)-induced male sex pheromone biosynthesis. The functional analysis of candidate enzymes involved in the male sex pheromone biosynthesis pathway demonstrated that a decreased mRNA levels of acetyl-CoA carboxylase, Δ11-desaturase, and fatty-acyl reductase 2 by RNAi-mediated knockdown led to a significant decrease in the production of fatty acyl alcohols and the efficacy of female acceptance. Our results demonstrated the important role of the fatty acyl alcohol biosynthetic pathway in a PBAN-induced male sex pheromone biosynthesis and the importance of hairpencil compounds in female mating acceptance.

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.

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.