Neofunctionalization in an ancestral insect desaturase lineage led to rare Δ6 pheromone signals in the Chinese tussah silkworm

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.

CRISPR/Cas9 mediated mutagenesis of the major sex pheromone gene, acyl-CoA delta-9 desaturase (DES9) in Fall armyworm Spodoptera frugiperda (J. E. Smith) (Lepidoptera: Noctuidae)

The Fall armyworm, Spodoptera frugiperda is a significant global pest causing serious yield loss on several staple crops. In this regard, this pest defies several management approaches based on chemicals, Bt transgenics etc., requiring effective alternatives. Recently CRISPR/Cas9 mediated genome editing has opened up newer avenues to establish functions of various target genes before employing them for further application. The virgin female moths of S. frugiperda emit sex pheromones to draw conspecific males. Therefore, we have edited the key pheromone synthesis gene, fatty acyl-CoA Delta-9 desaturase (DES9) of the Indian population of S. frugiperda. In order to achieve a larger deletion of the DES9, we have designed two single guide RNA (sgRNA) in sense and antisense direction targeting the first exon instead of a single guide RNA. The sgRNA caused site-specific knockout with a larger deletion which impacted the mating. Crossing studies between wild male and mutant female resulted in no fecundity, while fecundity was normal when mutant male crossed with the wild female. This indicates that mating disruption is stronger in females where DES9 is mutated. The current work is the first of its kind to show that DES9 gene editing impacted the likelihood of mating in S. frugiperda.

Effect of Irradiation on Reproduction of Female Spodoptera litura (Fabr.) (Lepidoptera: Noctuidae) in Relation to the Inherited Sterility Technique

Radiobiological investigations on the reproductive behavior of female Spodoptera litura (Fabr.) were conducted with the aim of determining the suitable radio-sterilizing dose for females in order to release them along with sub-sterile males for effective implementation of the Inherited Sterility technique against this pest. Calling and copulation duration significantly increased, while mating success, oviposition, fertility and longevity significantly decreased with increasing radiation dose (100-200 Gy) compared to control. In view of the effect of irradiation on mating behavior and reproductive viability of female S. litura, 130 Gy was identified as a suitable radio-sterilization dose. Further molecular studies were conducted to corroborate this dose for female sterilization, along with a higher dose of 200 Gy in order to validate the gradational response of ionizing radiation. GC-MS analysis indicated decreased sex pheromone titer at 130 Gy, which was more pronounced at 200 Gy. Pheromone-associated genes, PBAN and PBAN-R showed decreased expression at 130 Gy, and were drastically reduced at 200 Gy. The fertility-related Vg gene also showed a negative correlation with radiation exposure. Based on these radiation responses of female S. litura, 130 Gy might be considered a suitable dose for complete female sterility and its inclusion in sterile insect programs against S. litura.

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.

Release of moth pheromone compounds from Nicotiana benthamiana upon transient expression of heterologous biosynthetic genes

BACKGROUND

Using genetically modified plants as natural dispensers of insect pheromones may eventually become part of a novel strategy for integrated pest management.

RESULTS

In the present study, we first characterized essential functional genes for sex pheromone biosynthesis in the rice stem borer Chilo suppressalis (Walker) by heterologous expression in Saccharomyces cerevisiae and Nicotiana benthamiana, including two desaturase genes CsupYPAQ and CsupKPSE and a reductase gene CsupFAR2. Subsequently, we co-expressed CsupYPAQ and CsupFAR2 together with the previously characterized moth desaturase Atr∆11 in N. benthamiana. This resulted in the production of (Z)-11-hexadecenol together with (Z)-11-hexadecenal, the major pheromone component of C. suppressalis. Both compounds were collected from the transformed N. benthamiana headspace volatiles using solid-phase microextraction. We finally added the expression of a yeast acetyltransferase gene ATF1 and could then confirm also (Z)-11-hexadecenyl acetate release from the plant.

CONCLUSIONS

Our results pave the way for stable transformation of plants to be used as biological pheromone sources in different pest control strategies.

Bioproduction of (Z,E)-9,12-tetradecadienyl acetate (ZETA), the major pheromone component of Plodia, Ephestia, and Spodoptera species in yeast

BACKGROUND

(Z,E)-9,12-tetradecadienyl acetate (ZETA, Z9,E12-14:OAc) is a major sex pheromone component for many stored-product moth species. This pheromone is used worldwide for mating disruption, detection, monitoring, and mass trapping in raw and processed food storage facilities. In this study, we demonstrate the biological production of ZETA pheromone by engineered yeast Saccharomyces cerevisiae.

RESULTS

We mined the pheromone gland transcriptome data of the almond moth, Ephestia (Cadra) cautella (Walker), to trace a novel E12 fatty acyl desaturase and expressed candidates heterologously in yeast and Sf9 systems. Furthermore, we demonstrated a tailor-made ZETA pheromone bioproduction in yeast through metabolic engineering using this E12 desaturase, in combination with three genes from various sources coding for a Z9 desaturase, a fatty acyl reductase, and an acetyltransferase, respectively. Electrophysiological assays (gas chromatography coupled to an electroantennographic detector) proved that the transgenic yeast-produced ZETA pheromone component elicits distinct antennal responses.

CONCLUSION

The reconstructed biosynthetic pathway in yeast efficiently produces ZETA pheromone, leaves an undetectable level of biosynthetic intermediates, and paves the way for the economically competitive high-demand ZETA pheromone's bioproduction technology for high-value storage pest control.

Evolution of the Sex Pheromone Communication System in Ostrinia Moths

Simple Summary

Moths typically rely on sex pheromone communication to find a mate. This involves the production of species-specific sex pheromones by females (the signaller) and the corresponding selective detection by conspecific males (the receiver). A key question in the evolution of the pheromone communication system is how the female signals can diversify and still be tracked by the receivers over the process of speciation. The genus Ostrinia, which comprises 20 species worldwide including several well-recognised agricultural pests, is an attractive model in the study of the evolution of sex pheromone communication, as the closely related species and strains provide an ideal example of ongoing speciation. This review presents a comprehensive overview of the research on pheromone communication in different Ostrinia species over the past four decades, from the identity and biosynthesis of pheromones in the females to the molecular and neuronal basis of the pheromone perception in males. The evolutionary insights from these studies are discussed and the directions for future research are outlined.

Abstract

It remains a conundrum in the evolution of sexual communication how the signals and responses can co-ordinate the changes during speciation. The genus Ostrinia contains several closely related species as well as distinctive strains with pheromone polymorphism and represents an example of ongoing speciation. Extensive studies in the genus, especially in the species the European corn borer O. nubilalis (ECB), the Asian corn borer O. furnacalis (ACB) and the adzuki bean borer O. scapulalis (ABB), have provided valuable insights into the evolution of sex pheromone communication. This review presents a comprehensive overview of the research on pheromone communication in different Ostrinia species over the past four decades, including pheromone identification and biosynthesis, the ligand profiles of pheromone receptor (PR) genes, the physiology of peripheral olfactory sensory neurons (OSNs) and the projection pattern to the antennal lobe. By integrating and comparing the closely related Ostrinia species and strains, it provides an evolutionary perspective on the sex pheromone communication in moths in general and also outlines the outstanding questions that await to be elucidated by future studies.

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.

CRISPR/Cas9 mutagenesis against sex pheromone biosynthesis leads to loss of female attractiveness in Spodoptera exigua, an insect pestt

Virgin female moths are known to release sex pheromones to attract conspecific males. Accurate sex pheromones are required for their chemical communication. Sex pheromones of Spodoptera exigua, a lepidopteran insect, contain unsaturated fatty acid derivatives having a double bond at the 12^th carbon position. A desaturase of S. exigua (SexiDES5) was proposed to have dual functions by forming double bonds at the 11^th and 12^th carbons to synthesize Z9,E12-tetradecedienoic acid, which could be acetylated to be a main sex pheromone component Z9,E12-tetradecenoic acetate (Z9E12-14:Ac). A deletion of SexiDES5 using CRISPR/Cas9 was generated and inbred to obtain homozygotes. Mutant females could not produce Z9E12-14:Ac along with Z9-14:Ac and Z11-14:Ac. Subsequently, pheromone extract of mutant females did not induce a sensory signal in male antennae. They failed to induce male mating behavior including hair pencil erection and orientation. In the field, these mutant females did not attract any males while control females attracted males. These results indicate that SexiDES5 can catalyze the desaturation at the 11th and 12th positions to produce sex pheromone components in S. exigua. This study also suggests an application of the genome editing technology to insect pest control by generating non-attractive female moths.

Developments in Fatty Acid-Derived Insect Pheromone Production Using Engineered Yeasts

The use of traditional chemical insecticides for pest control often leads to environmental pollution and a decrease in biodiversity. Recently, insect sex pheromones were applied for sustainable biocontrol of pests in fields, due to their limited adverse impacts on biodiversity and food safety compared to that of other conventional insecticides. However, the structures of insect pheromones are complex, and their chemical synthesis is not commercially feasible. As yeasts have been widely used for fatty acid-derived pheromone production in the past few years, using engineered yeasts may be promising and sustainable for the low-cost production of fatty acid-derived pheromones. The primary fatty acids produced by Saccharomyces cerevisiae and other yeasts are C16 and C18, and it is also possible to rewire/reprogram the metabolic flux for other fatty acids or fatty acid derivatives. This review summarizes the fatty acid biosynthetic pathway in S. cerevisiae and recent progress in yeast engineering in terms of metabolic engineering and synthetic biology strategies to produce insect pheromones. In the future, insect pheromones produced by yeasts might provide an eco-friendly pest control method in agricultural fields.

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

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

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.

Evolution of the codling moth pheromone via an ancient gene duplication

BACKGROUND

Defining the origin of genetic novelty is central to our understanding of the evolution of novel traits. Diversification among fatty acid desaturase (FAD) genes has played a fundamental role in the introduction of structural variation in fatty acyl derivatives. Because of its central role in generating diversity in insect semiochemicals, the FAD gene family has become a model to study how gene family expansions can contribute to the evolution of lineage-specific innovations. Here we used the codling moth (Cydia pomonella) as a study system to decipher the proximate mechanism underlying the production of the ∆8∆10 signature structure of olethreutine moths. Biosynthesis of the codling moth sex pheromone, (E8,E10)-dodecadienol (codlemone), involves two consecutive desaturation steps, the first of which is unusual in that it generates an E9 unsaturation. The second step is also atypical: it generates a conjugated diene system from the E9 monoene C12 intermediate via 1,4-desaturation.

RESULTS

Here we describe the characterization of the FAD gene acting in codlemone biosynthesis. We identify 27 FAD genes corresponding to the various functional classes identified in insects and Lepidoptera. These genes are distributed across the C. pomonella genome in tandem arrays or isolated genes, indicating that the FAD repertoire consists of both ancient and recent duplications and expansions. Using transcriptomics, we show large divergence in expression domains: some genes appear ubiquitously expressed across tissue and developmental stages; others appear more restricted in their expression pattern. Functional assays using heterologous expression systems reveal that one gene, Cpo_CPRQ, which is prominently and exclusively expressed in the female pheromone gland, encodes an FAD that possesses both E9 and ∆8∆10 desaturation activities. Phylogenetically, Cpo_CPRQ clusters within the Lepidoptera-specific ∆10/∆11 clade of FADs, a classic reservoir of unusual desaturase activities in moths.

CONCLUSIONS

Our integrative approach shows that the evolution of the signature pheromone structure of olethreutine moths relied on a gene belonging to an ancient gene expansion. Members of other expanded FAD subfamilies do not appear to play a role in chemical communication. This advises for caution when postulating the consequences of lineage-specific expansions based on genomics alone.

Biosynthesis of the Sex Pheromone Component (E,Z)-7,9-Dodecadienyl Acetate in the European Grapevine Moth, Lobesia botrana, Involving ∆11 Desaturation and an Elusive ∆7 Desaturase

The European grapevine moth, Lobesia botrana, uses (E,Z)-7,9-dodecadienyl acetate as its major sex pheromone component. Through in vivo labeling experiments we demonstrated that the doubly unsaturated pheromone component is produced by ∆11 desaturation of tetradecanoic acid, followed by chain shortening of (Z)-11-tetradecenoic acid to (Z)-9-dodecenoic acid, and subsequently introduction of the second double bond by an unknown ∆7 desaturase, before final reduction and acetylation. By sequencing and analyzing the transcriptome of female pheromone glands of L. botrana, we obtained 41 candidate genes that may be involved in sex pheromone production, including the genes encoding 17 fatty acyl desaturases, 13 fatty acyl reductases, 1 fatty acid synthase, 3 acyl-CoA oxidases, 1 acetyl-CoA carboxylase, 4 fatty acid transport proteins and 2 acyl-CoA binding proteins. A functional assay of desaturase and acyl-CoA oxidase gene candidates in yeast and insect cell (Sf9) heterologous expression systems revealed that Lbo_PPTQ encodes a ∆11 desaturase producing (Z)-11-tetradecenoic acid from tetradecanoic acid. Further, Lbo_31670 and Lbo_49602 encode two acyl-CoA oxidases that may produce (Z)-9-dodecenoic acid by chain shortening (Z)-11-tetradecenoic acid. The gene encoding the enzyme introducing the E7 double bond into (Z)-9-dodecenoic acid remains elusive even though we assayed 17 candidate desaturases in the two heterologous systems.

Transcriptome Analysis of Ostrinia furnacalis Female Pheromone Gland: Esters Biosynthesis and Requirement for Mating Success

Female moths use sex pheromones to attract males, and corresponding regulatory mechanism underlying sex pheromone biosynthesis is species-dependent. However, the detailed mechanism involved in sex pheromone biosynthesis in Ostrinia furnacalis has not yet been fully addressed. In the present study, transcriptome sequencing of O. furnacalis pheromone glands screened a serials of candidate genes involved in sex pheromone biosynthesis. Our analysis showed that sex pheromone release in O. furnacalis females arrives its peak at the 2^nd scotophase, consistent with its mating behavior. Pheromone biosynthesis-activating neuropeptide (PBAN) was confirmed to regulate sex pheromone biosynthesis, and Ca²⁺ is the secondary messenger of PBAN signaling in O. furnacalis. The functional analysis of candidate genes demonstrated that the decreased mRNA levels or activities of calcineurin (CaN) and acetyl-CoA carboxylase (ACC) led to significant decrease in sex pheromone production and female capability to attract males, as demonstrated by RNAi-mediated knockdown and pharmacological inhibitor assay. Most importantly, the activities of CaN and ACC depend on the activation of PBAN/PBANR/Ca²⁺. Furthermore, fatty-acyl reductase 14 was involved in PBAN-mediated sex pheromone biosynthesis. Altogether, our results demonstrated that PBAN regulates sex pheromone biosynthesis through PBANR/Ca²⁺/CaN/ACC pathway to promote sex pheromone biosynthesis in O. furnacalis and provided a reference for non-model organism to study neuropeptide signal transduction.

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.

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

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

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.

Identification of Δ6-unsaturated, monoenyl type I pheromone compounds from the cashew stem borer Anthistarcha binocularis (Lepidoptera: Gelechiidae)

BACKGROUND

The cashew stem borer Anthistarcha binocularis (Meyrick) is a major pest in cashew orchards in Brazil. The damage caused by the larvae results in economic losses, and the available chemical control is not suitable because of the endophytic nature of A. binocularis. The identification of the A. binocularis sex pheromone could provide novel applications for the detection and monitoring of this species.

RESULTS

Two compounds from female sex gland extracts elicited electrophysiological responses on male antennae. They were identified as dodec-6-en-1-ol and dodec-6-en-1-yl acetate by gas chromatography-mass spectrometry. E/Z stereoisomers of both compounds were synthesized, and the electroantennograms for the synthetic compounds showed the strongest responses for the (E)-stereoisomers of the alcohol and acetate. In a field trial, the E6-12:OH/E6-12:OAc mixture attracted male A. binocularis, whereas the Z6-12:OH/Z6-12:OAc mixture attracted no specimens.

CONCLUSIONS

The bioactive compounds from the sex pheromone of A. binocularis have been identified as a mixture of two previously unidentified pheromone compounds: E6-12:OH and E6-12:OAc. The mixture of both compounds was attractive to males in preliminary field experiments, and this study is the first report of Δ6-unsaturated monoenyl pheromone compounds in Lepidoptera. © 2019 Society of Chemical Industry.

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.

Carbon-carbon double bond position elucidation in fatty acids using ozone-coupled direct analysis in real time mass spectrometry

The carbon-carbon double bond positions of unsaturated fatty acids can have markedly different effects on biological function and also serve as biomarkers of disease pathology, dietary history, and species identity. As such, there is great interest in developing methods for the facile determination of double bond position for natural product chemistry, the pharmaceutical industry, and forensics. We paired ozonolysis with direct analysis in real time mass spectrometry (DART MS) to cleave and rapidly identify carbon-carbon double bond position in fatty acids, fatty alcohols, wax esters, and crude fatty acid extracts. In addition, ozone exposure time and DART ion source temperature were investigated to identify optimal conditions. Our results reveal that brief, offline exposure to ozone-generated aldehyde and carboxylate products that are indicative of carbon-carbon double bond position. The relative abundance of diagnostic fragments quantitatively reflects the ratios of isobaric fatty acid positional isomers in a mixture with a correlation coefficient of 0.99. Lastly, the unsaturation profile generated from unfractionated, fatty acid extracts can be used to differentiate insect species and populations. The ability to rapidly elucidate lipid double bond position by combining ozonolysis with DART MS will be useful for lipid structural elucidation, assessing isobaric purity, and potentially distinguishing between animals fed on different diets or belonging to different ecological populations.

A Δ9 desaturase (SlitDes11) is associated with the biosynthesis of ester sex pheromone components in Spodoptera litura

Sex pheromone biosynthesis in moths relies on the activity of multiple enzymes, including Δ9 desaturase, which plays an important role in catalyzing desaturation at the Δ9 position of the carbon chain. However, the physiological function of moth Δ9 desaturase has not been elucidated in vivo. In this study, we used the CRISPR/Cas9 system to knockout the Δ9 desaturase gene (SlitDes11) of Spodoptera litura to analyze its role in sex pheromone biosynthesis. First, through the direct injection of SlitDes11-single guide RNA (sgRNA)/Cas9 messenger RNA into newly laid eggs, gene editing was induced in around 30% of eggs 24 h after injection and was induced in 20.8% of the resulting adult moths. Second, using a sibling-crossing strategy, insects with mutant SlitDes11 (bearing a premature stop codon) were selected, and homozygous mutants were obtained in the G5 generation. Third, pheromone gland extracts of adult female homozygous SlitDes11 mutants were analyzed using Gas chromatography (GC). The results showed that titers of all three ester sex pheromone components; Z9, E11-14:Ac, Z9,E12-14:Ac, and Z9-14:Ac; were reduced by 62.40%, 78.50%, and 72.50%, respectively. This study provides the first direct evidence for the role of SlitDes11 in sex pheromone biosynthesis in S. litura, and indicates the gene could be as potential target to disrupt sexual communication in S. litura for developing a new pollution-free insecticide.

Multi-Functional Desaturases in Two Spodoptera Moths with ∆11 and ∆12 Desaturation Activities

The beet armyworm, Spodoptera exigua, uses (Z,E)-9,12-tetradecadienyl acetate as the major component of its sex pheromone. Previous isotope-labeling experiments demonstrated an unusual ∆12 desaturase activity involved in the biosynthesis of this compound; however, the putative ∆12 desaturase gene has not been reported to date. In the present study, we confirmed this ∆12 desaturation pathway by in vivo labeling experiments, and characterized candidate desaturase genes in a yeast heterologous expression system. We demonstrated that a pheromone gland-specific desaturase, SexiDes5, uses palmitic acid and the subsequently chain-shortened product (Z)-9-tetradecenoic acid as substrates to produce (Z)-11-hexadecenoic and (Z,E)-9,12-tetradecadienoic acids, respectively. In addition, the homologous desaturase SlitDes5 from the congeneric Spodoptera litura had analogous functions.

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.

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.

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.

Sex Attractant Pheromone of the Luna Moth, Actias luna (Linnaeus)

Giant silk moths (Lepidoptera: Saturniidae) typically are not well represented as larvae or adults in community level inventories of Lepidoptera, and as a result, little is known about their population dynamics. Furthermore, in recent years, many species of silk moths appear to have experienced population declines. Volatile sex pheromones are powerful sampling tools that can be used in operational conservation and monitoring programs for insects. Here, we describe the identification of the sex attractant pheromone of a giant silk moth, the luna moth Actias luna. Coupled gas chromatography-electroantennographic detection and gas chromatography-mass spectrometric analyses of extracts from pheromone glands of female luna moths supported the identification of (6E,11Z)-6,11-octadecadienal (E6,Z11-18:Ald), (6E)-6-octadecenal (E6-18:Ald), and (11Z)-11-octadecenal (Z11-18:Ald) as the compounds in extracts that elicited responses from antennae of male moths. These identifications were confirmed by synthesis, followed by testing of blends of the synthetic compounds in field trials in Ontario, Canada, and Kentucky, USA. Male moths were attracted to synthetic E6,Z11-18:Ald as a single component. Attraction appeared to be enhanced by addition of E6-18:Ald but not Z11-18:Ald, suggesting that the luna moth pheromone consists of a blend of E6,Z11-18:Ald and E6-18:Ald.

Sequence variation determining stereochemistry of a Δ11 desaturase active in moth sex pheromone biosynthesis

A Δ11 desaturase from the oblique banded leaf roller moth Choristoneura rosaceana takes the saturated myristic acid and produces a mixture of (E)-11-tetradecenoate and (Z)-11-tetradecenoate with an excess of the Z isomer (35:65). A desaturase from the spotted fireworm moth Choristoneura parallela also operates on myristic acid substrate but produces almost pure (E)-11-tetradecenoate. The two desaturases share 92% amino acid identity and 97% amino acid similarity. There are 24 amino acids differing between these two desaturases. We constructed mutations at all of these positions to pinpoint the sites that determine the product stereochemistry. We demonstrated with a yeast functional assay that one amino acid at the cytosolic carboxyl terminus of the protein (258E) is critical for the Z activity of the C. rosaceana desaturase. Mutating the glutamic acid (E) into aspartic acid (D) transforms the C. rosaceana enzyme into a desaturase with C. parallela-like activity, whereas the reciprocal mutation of the C. parallela desaturase transformed it into an enzyme producing an intermediate 64:36 E/Z product ratio. We discuss the causal link between this amino acid change and the stereochemical properties of the desaturase and the role of desaturase mutations in pheromone evolution.

Evolution of moth sex pheromone composition by a single amino acid substitution in a fatty acid desaturase

For sexual communication, moths primarily use blends of fatty acid derivatives containing one or more double bonds in various positions and configurations, called sex pheromones (SPs). To study the molecular basis of novel SP component (SPC) acquisition, we used the tobacco hornworm (Manduca sexta), which uses a blend of mono-, di-, and uncommon triunsaturated fatty acid (3UFA) derivatives as SP. We identified pheromone-biosynthetic fatty acid desaturases (FADs) MsexD3, MsexD5, and MsexD6 abundantly expressed in the M. sexta female pheromone gland. Their functional characterization and in vivo application of FAD substrates indicated that MsexD3 and MsexD5 biosynthesize 3UFAs via E/Z14 desaturation from diunsaturated fatty acids produced by previously characterized Z11-desaturase/conjugase MsexD2. Site-directed mutagenesis of sequentially highly similar MsexD3 and MsexD2 demonstrated that swapping of a single amino acid in the fatty acyl substrate binding tunnel introduces E/Z14-desaturase specificity to mutated MsexD2. Reconstruction of FAD gene phylogeny indicates that MsexD3 was recruited for biosynthesis of 3UFA SPCs in M. sexta lineage via gene duplication and neofunctionalization, whereas MsexD5 representing an alternative 3UFA-producing FAD has been acquired via activation of a presumably inactive ancestral MsexD5. Our results demonstrate that a change as small as a single amino acid substitution in a FAD enzyme might result in the acquisition of new SP compounds.

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.

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.

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.

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

Differences in sex pheromone component can lead to reproductive isolation. The sibling noctuid species, Helicoverpa armigera and Helicoverpa assulta, share the same two sex pheromone components, Z9-16:Ald and Z11-16:Ald, but in opposite ratios, providing an typical example of such reproductive isolation. To investigate how the ratios of the pheromone components are differently regulated in the two species, we sequenced cDNA libraries from the pheromone glands of H. armigera and H. assulta. After assembly and annotation, we identified 108 and 93 transcripts putatively involved in pheromone biosynthesis, transport, and degradation in H. armigera and H. assulta, respectively. Semi-quantitative RT-PCR, qRT-PCR, phylogenetic, and mRNA abundance analyses suggested that some of these transcripts involved in the sex pheromone biosynthesis pathways perform. Based on these results, we postulate that the regulation of desaturases, KPSE and LPAQ, might be key factor regulating the opposite component ratios in the two sibling moths. In addition, our study has yielded large-scale sequence information for further studies and can be used to identify potential targets for the bio-control of these species by disrupting their sexual communication.

Large number of putative chemoreception and pheromone biosynthesis genes revealed by analyzing transcriptome from ovipositor-pheromone glands of Chilo suppressalis

The chemoreception role of moth ovipositor has long been suggested, but its molecular mechanism is mostly unknown. By transcriptomic analysis of the female ovipositor-pheromone glands (OV-PG) of Chilo suppressalis, we obtained 31 putative chemoreception genes (9 OBPs, 10 CSPs, 2 ORs, 1 SNMP, 8 CXEs and 1 AOX), in addition to 32 genes related to sex pheromone biosynthesis (1 FAS, 6 Dess, 10 FARs, 2 ACOs, 1 ACC, 4 FATPs, 3 ACBPs and 5 ELOs). Tissue expression profiles further revealed that CsupCSP2 and CsupCSP10 were OV-PG biased, while most chemoreception genes were highly and preferably expressed in antennae. This suggests that OV-PG employs mostly the same chemoreception proteins as in antennae, although the physiological roles of these proteins might be different in OV-PG. Of the 32 pheromone biosynthesis related genes, CsupDes4, CsupDes5 and CsupFAR2 are strongly OV-PG biased, and clustered with functionally validated genes from other moths, strongly indicating their involvement in specific step of the pheromone biosynthesis. Our study for the first time identified a large number of putative chemoreception genes, and provided an important basis for exploring the chemoreception mechanisms of OV-PG in C. suppressalis, as well as other moth species.

Female sex pheromones of two Japanese saturniid species, Rhodinia fugax and Loepa sakaei: identification, synthesis, and field evaluation

While 11 species in the family Saturniidae are found in Japan, no sex pheromones of the native species had been investigated previously. We collected larvae of Rhodinia fugax in Nagano and Tottori Prefecture, and of Loepa sakaei in Okinawa Prefecture, and extracted sex pheromones of these two species from virgin female moths. In gas chromatography-electroantennogram detection (GC-EAD) analyses, male antennae of each species responded to one component in the respective pheromone extracts of conspecific females. Chemical analyses of the extracts by GC/mass spectrometry revealed that the EAD-active compounds of R. fugax and L. sakaei were a hexadecadienal and a tetradecadienyl acetate, respectively. The two species belong to the subfamily Saturniinae, and the mass spectra of both were similar to that of the 6,11-hexadecadienyl acetate identified from Antheraea polyphemus, classified in the same subfamily, suggesting the same 6,11-dienyl structure for the C16 aldehyde and a 4,9-dienyl structure for the C14 acetate. Based on this assumption, four geometrical isomers of each dienyl compound were stereoselectively synthesized via acetylene intermediates, compared to the natural products, and tested in the field. Male catches confirmed the pheromone structures of the two Japanese saturniid species as (6E,11Z)-6,11-hexadecadienal for R. fugax and (4E,9Z)-4,9-tetradecadienyl acetate for L. sakaei. The compounds have a characteristic 1,6-dienyl motif common to the pheromones of Saturniinae species.

Unexpected functional diversity in the fatty acid desaturases of the flour beetle Tribolium castaneum and identification of key residues determining activity

Desaturases catalyse modifications to fatty acids which are essential to homeostasis and for pheromone and defensive chemical production. All desaturases of the flour beetle Tribolium castaneum were investigated via query of the sequenced genome which yielded 15 putative acyl-Coenzyme A genes. Eleven desaturase mRNA were obtained in full length and functionally expressed in yeast. Phylogenetic analysis separated the desaturases into 4 distinct clades; one clade contained conserved beetle Δ9 desaturases, second clade was Tribolium-specific having diverse activities including Δ5, Δ9 and Δ12 desaturation and the other 2 clades had mixed insect representatives. Three members of this clade contained unusual inserted sequences of ∼20 residues in the C-terminal region and were related to desaturases that all contained similar inserts. Deletion of the entirety of the insert in the flour beetle Δ12 desaturase abolished its activity but this was partially restored by the reintroduction of two histidine residues, suggesting the histidine(s) are required for activity but the full length insert is not. Five new desaturase activities were discovered: Δ9 desaturation of C12:0-C16:0 substrates; two unprecedented Δ5 enzymes acting on C18:0 and C16:0; Δ9 activity exclusively on C16:0 and a further stearate Δ9 desaturase. qPCR analysis ruled out a role in sex pheromone synthesis for the Δ5 and Δ9/C16:0 desaturases. The flour beetle genome has underpinned an examination of all transcribed desaturases in the organism and revealed a diversity of novel and unusual activities, an improved understanding of the evolutionary relationships among insect desaturases and sequence determinants of activity.

Sex pheromone biosynthetic pathways are conserved between moths and the butterfly Bicyclus anynana

Although phylogenetically nested within the moths, butterflies have diverged extensively in a number of life history traits. Whereas moths rely greatly on chemical signals, visual advertisement is the hallmark of mate finding in butterflies. In the context of courtship, however, male chemical signals are widespread in both groups although they likely have multiple evolutionary origins. Here, we report that in males of the butterfly Bicyclus anynana, courtship scents are produced de novo via biosynthetic pathways shared with females of many moth species. We show that two of the pheromone components that play a major role in mate choice, namely the (Z)-9-tetradecenol and hexadecanal, are produced through the activity of a fatty acyl Δ11-desaturase and two specialized alcohol-forming fatty acyl reductases. Our study provides the first evidence of conservation and sharing of ancestral genetic modules for the production of FA-derived pheromones over a long evolutionary timeframe thereby reconciling mate communication in moths and butterflies.

Little is known about the evolutionary origins of the genes involved in butterfly pheromone synthesis. Here, Liénard et al. show that the biosynthetic pathways involved in the production of male courtship scents of the butterfly, Bicyclus anynana, are shared with females of many moth species.

Putative pathway of sex pheromone biosynthesis and degradation by expression patterns of genes identified from female pheromone gland and adult antenna of Sesamia inferens (Walker)

The general pathway of biosynthesis and degradation for Type-I sex pheromones in moths is well established, but some genes involved in this pathway remain to be characterized. The purple stem borer, Sesamia inferens, employs a pheromone blend containing components with three different terminal functional groups (Z11-16:OAc, Z11-16:OH, and Z11-16:Ald) of Type-I sex pheromones. Thus, it provides a good model to study the diversity of genes involved in pheromone biosynthesis and degradation pathways. By analyzing previously obtained transcriptomic data of the sex pheromone glands and antennae, we identified 73 novel genes that are possibly related to pheromone biosynthesis (46 genes) or degradation (27 genes). Gene expression patterns and phylogenetic analysis revealed that one desaturase (SinfDes4), one fatty acid reductase (SinfFAR2), and one fatty acid xtransport protein (SinfFATP1) genes were predominantly expressed in pheromone glands, and clustered with genes involved in pheromone synthesis in other moth species. Ten genes including five carboxylesterases (SinfCXE10, 13, 14, 18, and 20), three aldehyde oxidases (SinfAOX1, 2 and 3), and two alcohol dehydrogenases (SinfAD1 and 3) were expressed specifically or predominantly in antennae, and could be candidate genes involved in pheromone degradation. SinfAD1 and 3 are the first reported alcohol dehydrogenase genes with antennae-biased expression. Based on these results we propose a pathway involving these potential enzyme-encoding gene candidates in sex pheromone biosynthesis and degradation in S. inferens. This study provides robust background information for further elucidation of the genetic basis of sex pheromone biosynthesis and degradation, and ultimately provides potential targets to disrupt sexual communication in S. inferens for control purposes.

A novel fatty acyl desaturase from the pheromone glands of Ctenopseustis obliquana and C. herana with specific Z5-desaturase activity on myristic acid

Sexual communication in the Lepidoptera typically involves a female-produced sex pheromone that attracts males of the same species. The most common type of moth sex pheromone comprises individual or blends of fatty acyl derivatives that are synthesized by a specific enzymatic pathway in the female’s pheromone gland, often including a desaturation step. This reaction is catalyzed by fatty acyl desaturases that introduce double bonds at specific locations in the fatty acid precursor backbone. The two tortricid moths, Ctenopseustis obliquana and C. herana (brown-headed leafrollers), which are endemic in New Zealand, both use (Z)-5-tetradecenyl acetate as part of their sex pheromone. In C. herana, (Z)-5-tetradecenyl acetate is the sole component of the pheromone. Labeling experiments have revealed that this compound is produced via an unusual Δ5-desaturation of myristic acid. Previously six desaturases were identified from the pheromone glands of Ctenopseustis and its sibling genus Planotortrix, with one differentially regulated to produce the distinct blends used by individual species. However, none were able to conduct the Δ5-desaturation observed in C. herana, and presumably C. obliquana. We have now identified an additional desaturase gene, desat7, expressed in the pheromone glands of both Ctenopseustis species, which is not closely related to any previously described moth pheromone desaturase. The encoded enzyme displays Δ5-desaturase activity on myristic acid when heterologously expressed in yeast, but is not able to desaturate any other fatty acid (C8–C16). We conclude that desat7 represents a new group of desaturases that has evolved a role in the biosynthesis of sex pheromones in moths.

Functional evolution of a multigene family: orthologous and paralogous pheromone receptor genes in the turnip moth, Agrotis segetum

Lepidopteran pheromone receptors (PRs), for which orthologies are evident among closely related species, provide an intriguing example of gene family evolution in terms of how new functions may arise. However, only a limited number of PRs have been functionally characterized so far and thus evolutionary scenarios suffer from elements of speculation. In this study we investigated the turnip moth Agrotis segetum, in which female moths produce a mixture of chemically related pheromone components that elicit specific responses from receptor cells on male antennae. We cloned nine A. segetum PR genes and the Orco gene by degenerate primer based RT-PCR. The nine PR genes, named as AsegOR1 and AsegOR3-10, fall into four distinct orthologous clusters of known lepidopteran PRs, of which one contains six paralogues. The paralogues are under relaxed selective pressure, contrasting with the purifying selection on other clusters. We identified the receptors AsegOR9, AsegOR4 and AsegOR5, specific for the respective homologous pheromone components (Z)-5-decenyl, (Z)-7-dodecenyl and (Z)-9-tetradecenyl acetates, by two-electrode voltage clamp recording from Xenopus laevis oocytes co-expressing Orco and each PR candidate. These receptors occur in three different orthologous clusters. We also found that the six paralogues with high sequence similarity vary dramatically in ligand selectivity and sensitivity. Different from AsegOR9, AsegOR6 showed a relatively large response to the behavioural antagonist (Z)-5-decenol, and a small response to (Z)-5-decenyl acetate. AsegOR1 was broadly tuned, but most responsive to (Z)-5-decenyl acetate, (Z)-7-dodecenyl acetate and the behavioural antagonist (Z)-8-dodecenyl acetate. AsegOR8 and AsegOR7, which differ from AsegOR6 and AsegOR1 by 7 and 10 aa respectively, showed much lower sensitivities. AsegOR10 showed only small responses to all the tested compounds. These results suggest that new receptors arise through gene duplication, and relaxed evolutionary constraints or positive selection among paralogues allow functional divergence to occur in spite of purifying selection being the norm.

Cloning and tissue distribution of a fatty acyl Δ6-desaturase-like gene and effects of dietary lipid levels on its expression in the hepatopancreas of Chinese mitten crab (Eriocheir sinensis)

Fatty acyl Δ6-desaturase is the rate-limiting enzyme in the biosynthetic pathway of highly unsaturated fatty acids (HUFAs) in vertebrates. In this report, a fatty acyl Δ6-desaturase-like cDNA was cloned from the hepatopancreas of Eriocheir sinensis (Chinese mitten crab) and characterized by performing rapid-amplification of cDNA ends. The 2278-bp long full-length cDNA encodes a polypeptide with 442 amino acids. Gene expression analysis via real-time quantitative polymerase chain reaction revealed that the fatty acyl Δ6-desaturase-like transcripts are widely distributed in various tissues, with high expression levels in the hepatopancreas and cranial ganglia. This study focuses on the nutritional regulation of genes involved in the HUFA biosynthetic pathway in Chinese mitten crab. A feeding trial was performed whereby crablets were fed for 238 d with four different diets: control diet without oil lipids (added with 3% basic lipid of the fundamental diets); fish oil diet (FO; added with 3% of the fundamental diets); soybean oil diet (SO; added with 3% of the fundamental diets); and FO/SO diet (1:1; added with 3% of the fundamental diets). The hepatopancreas of crabs sampled at 168 d and 238 d to determine the effects on fatty acyl Δ6-desaturase-like mRNA expression. The results show that the expression of fatty acyl Δ6-desaturase-like is higher in the hepatopancreas of crabs fed with SO diet than those fed with FO diet. Furthermore, gene expression increased by 2.45-fold in the hepatopancreas of crabs fed with SO after 238 d than those fed after 168 d but remained steady for those fed with FO after 238 d.

Semi-selective fatty acyl reductases from four heliothine moths influence the specific pheromone composition

Background

Sex pheromones are essential in moth mate communication. Information on pheromone biosynthetic genes and enzymes is needed to comprehend the mechanisms that contribute to specificity of pheromone signals. Most heliothine moths use sex pheromones with (Z)–11–hexadecenal as the major component in combination with minor fatty aldehydes and alcohols. In this study we focus on four closely related species, Heliothis virescens, Heliothis subflexa, Helicoverpa armigera and Helicoverpa assulta, which use (Z)–11–hexadecenal, (Z)–9–tetradecanal, and (Z)–9–hexadecenal in different ratios in their pheromone blend. The components are produced from saturated fatty acid precursors by desaturation, β–oxidation, reduction and oxidation.

Results

We analyzed the composition of fatty acyl pheromone precursors and correlated it to the pheromone composition. Next, we investigated whether the downstream fatty–acyl reduction step modulates the ratio of alcohol intermediates before the final oxidation step. By isolating and functionally characterizing the Fatty Acyl Reductase (pgFAR) from each species we found that the pgFARs were active on a broad set of C8 to C16 fatty acyl substrates including the key pheromone precursors, Z9–14, Z9–16 and Z11–16:acyls. When presenting the three precursors in equal ratios to yeast cultures expressing any of the four pgFARs, all reduced (Z)–9–tetradecenoate preferentially over (Z)–11–hexadecenoate, and the latter over (Z)–9–hexadecenoate. Finally, when manipulating the precursor ratios in vitro, we found that the pgFARs display small differences in the biochemical activity on various substrates.

Conclusions

We conclude that a pgFAR with broad specificity is involved in heliothine moth pheromone biosynthesis, functioning as a semi–selective funnel that produces species–specific alcohol product ratios depending on the fatty–acyl precursor ratio in the pheromone gland. This study further supports the key role of these in pheromone biosynthesis and emphasizes the interplay between the pheromone fatty acyl precursors and the Lepidoptera specific pgFARs in shaping the pheromone composition.

Sex pheromone evolution is associated with differential regulation of the same desaturase gene in two genera of leafroller moths

Chemical signals are prevalent in sexual communication systems. Mate recognition has been extensively studied within the Lepidoptera, where the production and recognition of species-specific sex pheromone signals are typically the defining character. While the specific blend of compounds that makes up the sex pheromones of many species has been characterized, the molecular mechanisms underpinning the evolution of pheromone-based mate recognition systems remain largely unknown. We have focused on two sets of sibling species within the leafroller moth genera Ctenopseustis and Planotortrix that have rapidly evolved the use of distinct sex pheromone blends. The compounds within these blends differ almost exclusively in the relative position of double bonds that are introduced by desaturase enzymes. Of the six desaturase orthologs isolated from all four species, functional analyses in yeast and gene expression in pheromone glands implicate three in pheromone biosynthesis, two Δ9-desaturases, and a Δ10-desaturase, while the remaining three desaturases include a Δ6-desaturase, a terminal desaturase, and a non-functional desaturase. Comparative quantitative real-time PCR reveals that the Δ10-desaturase is differentially expressed in the pheromone glands of the two sets of sibling species, consistent with differences in the pheromone blend in both species pairs. In the pheromone glands of species that utilize (Z)-8-tetradecenyl acetate as sex pheromone component (Ctenopseustis obliquana and Planotortrix octo), the expression levels of the Δ10-desaturase are significantly higher than in the pheromone glands of their respective sibling species (C. herana and P. excessana). Our results demonstrate that interspecific sex pheromone differences are associated with differential regulation of the same desaturase gene in two genera of moths. We suggest that differential gene regulation among members of a multigene family may be an important mechanism of molecular innovation in sex pheromone evolution and speciation.

Chemical signals are prevalent in sexual communication systems, especially within the Lepidoptera where sex pheromones are typically one of the defining characteristics of species. We have isolated six desaturases from two groups of sibling species of leafroller moths belonging to the genera Ctenopseustis and Planotortrix. Functional analyses in yeast and quantitative RT–PCR indicate that three of the desaturases are involved in the biosynthesis of sex pheromone components in these species. One of three enzymes is a Δ10-desaturase that is differentially expressed in the pheromone glands of the two sets of sibling species, consistent with differences in the pheromone blend in both species pairs. In the pheromone glands of species that utilize (Z)-8-tetradecenyl acetate as sex pheromone component (C. obliquana and P. octo), the expression levels of the Δ10-desaturase are significantly higher than pheromone gland expression levels in their sibling species (C. herana and P. excessana). Our results demonstrate that interspecific sex pheromone differences are associated with differential regulation of the same desaturase gene in these two genera of moths. Based on these findings differential gene regulation among members of a multigene family may be an important mechanism of molecular innovation in sex pheromone evolution and speciation.

Terminal fatty-acyl-CoA desaturase involved in sex pheromone biosynthesis in the winter moth (Operophtera brumata)

The winter moth (Operophtera brumata L., Lepidoptera: Geometridae) utilizes a single hydrocarbon, 1,Z3,Z6,Z9-nonadecatetraene, as its sex pheromone. We tested the hypothesis that a fatty acid precursor, Z11,Z14,Z17,19-nonadecanoic acid, is biosynthesized from α-linolenic acid, through chain elongation by one 2-carbon unit, and subsequent methyl-terminus desaturation. Our results show that labeled α-linolenic acid is indeed incorporated into the pheromone component in vivo. A fatty-acyl-CoA desaturase gene that we found to be expressed in the abdominal epidermal tissue, the presumed site of biosynthesis for type II pheromones, was characterized and expressed heterologously in a yeast system. The transgenic yeast expressing this insect derived gene could convert Z11,Z14,Z17-eicosatrienoic acid into Z11,Z14,Z17,19-eicosatetraenoic acid. These results provide evidence that a terminal desaturation step is involved in the winter moth pheromone biosynthesis, prior to the decarboxylation.