Evolution of multicomponent pheromone signals in small ermine moths involves a single fatty-acyl reductase gene

Regulatory Changes in the Fatty Acid Elongase eloF Underlie the Evolution of Sex-specific Pheromone Profiles in Drosophila prolongata

Pheromones play a key role in regulating sexual behavior throughout the animal kingdom. In Drosophila and other insects, many cuticular hydrocarbons (CHCs) are sexually dimorphic, and some are known to perform pheromonal functions. However, the genetic control of sex-specific CHC production is not understood outside of the model species D. melanogaster. A recent evolutionary change is found in D. prolongata, which, compared to its closest relatives, shows greatly increased sexual dimorphism in both CHCs and the chemosensory system responsible for their perception. A key transition involves a male-specific increase in the proportion of long-chain CHCs. Perfuming D. prolongata females with the male-biased CHCs reduces copulation success, suggesting that these compounds function as sex pheromones. The evolutionary change in CHC profiles correlates with a male-specific increase in the expression of multiple genes involved in CHC biosynthesis, including fatty acid elongases and reductases and other key enzymes. In particular, elongase F, which is responsible for producing female-specific pheromones in D. melanogaster, is strongly upregulated in D. prolongata males compared both to females and to males of the sibling species. Induced mutations in eloF reduce the amount of long-chain CHCs, resulting in a partial feminization of pheromone profiles in D. prolongata males while having minimal effect in females. Transgenic experiments show that sex-biased expression of eloF is caused in part by a putative transposable element insertion in its regulatory region. These results reveal one of the genetic mechanisms responsible for a recent evolutionary change in sexual communication.

Fatty Acid Origin of Insect Pheromones

Pheromones are utilized to a great extent in insects. Many of these pheromones are biosynthesized through a pathway involving fatty acids. This chapter will provide examples where the biosynthetic pathways of fatty acid-derived pheromones have been studied in detail. These include pheromones from Lepidoptera, Coleoptera, and Hymenoptera. Many species of Lepidoptera utilize fatty acids as precursors to pheromones with a functional group that include aldehydes, alcohols, and acetate esters. In addition, the biosynthesis of hydrocarbons will be briefly examined because many insects utilize hydrocarbons or modified hydrocarbons as pheromones.

A novel domain-duplicated SlitFAR3 gene involved in sex pheromone biosynthesis in Spodoptera litura

Fatty acyl reductases (FARs) are key enzymes that participate in sex pheromone biosynthesis by reducing fatty acids to fatty alcohols. Lepidoptera typically harbor numerous FAR gene family members. Although FAR genes are involved in the biosynthesis of sex pheromones in moths, the key FAR gene of Spodoptera litura remains unclear. In this work, we predicted 30 FAR genes from the S. litura genome and identified a domain duplication within gene SlitFAR3, which exhibited high and preferential expression in the sexually mature female pheromone glands (PGs) and a rhythmic expression pattern during the scotophase of sex pheromone production. The molecular docking of SlitFAR3, as predicted using a 3D model, revealed a co-factor NADPH binding cavity and 2 substrate binding cavities. Functional expression in yeast cells combined with comprehensive gas chromatography indicated that the SlitFAR3 gene could produce fatty alcohol products. This study is the first to focus on the special phenomenon of FAR domain duplication, which will advance our understanding of biosynthesis-related genes from the perspective of evolutionary biology.

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

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

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

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

Cathepsin L Contributes to Reproductive Diapause by Regulating Lipid Storage and Survival of Coccinella septempunctata (Linnaeus)

Cathepsin L protease, which belongs to the papain-like cysteine proteases family, is an important player in many physiological and pathological processes. However, little was known about the role of cathepsin L in ladybird beetles (Coccinella septempuctata Linnaeus) during diapause. Here, we analyzed the characteristics of cathepsin L (CsCatL) in the females of C. septempunctata and its role during the diapause of the ladybeetle. CsCatL was cloned and identified from beetle specimens by rapid amplification of cDNA-ends (RACE). The cDNA sequence of CsCatL was 971 bp in length, including an 843 bp open reading frame encoding a protein of 280 amino acids. It was identified as the cathepsin L group by phylogenetic analysis. Knockdown of CsCatL by RNA interference led to decreased expression levels of fatty acid synthase 2 (fas 2) genes and suppressed lipid accumulation. Furthermore, silencing the CsCatL gene distinctly reduced diapause-related features and the survival of female C. spetempunctata under diapause-inducing conditions. The results suggested that the CsCatL gene was involved in fatty acid biosynthesis and played a crucial role in the survival of adult C. septempunctata during the diapause preparation stage.

Fatty acyl-CoA reductase influences wax biosynthesis in the cotton mealybug, Phenacoccus solenopsis Tinsley

Mealybugs are highly aggressive to a diversity of plants. The waxy layer covering the outermost part of the integument is an important protective defense of these pests. However, the molecular mechanisms underlying wax biosynthesis in mealybugs remain largely unknown. Here, we analyzed multi-omics data on wax biosynthesis by the cotton mealybug, Phenacoccus solenopsis Tinsley, and found that a fatty acyl-CoA reductase (PsFAR) gene, which was highly expressed in the fat bodies of female mealybugs, contributed to wax biosynthesis by regulating the production of the dominant chemical components of wax, cuticular hydrocarbons (CHCs). RNA interference (RNAi) against PsFAR by dsRNA microinjection and allowing mealybugs to feed on transgenic tobacco expressing target dsRNA resulted in a reduction of CHC contents in the waxy layer, and an increase in mealybug mortality under desiccation and deltamethrin treatments. In conclusion, PsFAR plays crucial roles in the wax biosynthesis of mealybugs, thereby contributing to their adaptation to water loss and insecticide stress.

The role of a fatty acyl-CoA reductase (PsFAR) in wax biosynthesis of cotton mealybug is investigated, RNAi against PsFAR resulted in insects with lower generation of waxy filaments and higher mortality under desiccation and deltamethrin treatments.

Natural variation at a single gene generates sexual antagonism across fitness components in Drosophila

Mutations with conflicting fitness effects in males and females accumulate in sexual populations, reducing their adaptive capacity.^1,2 Although quantitative genetic studies indicate that sexually antagonistic polymorphisms are common,^3-5 their molecular basis and population genetic properties remain poorly understood.^6,7 Here, we show in fruit flies how natural variation at a single gene generates sexual antagonism through phenotypic effects on cuticular hydrocarbon (CHC) traits that function as both mate signals and protectors against abiotic stress⁸ across a latitudinal gradient. Tropical populations of Drosophila serrata have polymorphic CHCs producing sexual antagonism through opposing but sex-limited effects on these two fitness-related functions. We dissected this polymorphism to a single fatty-acyl CoA reductase gene, DsFAR2-B, that is expressed in oenocyte cells where CHCs are synthesized. RNAi-mediated disruption of the DsFAR2-B ortholog in D. melanogaster oenocytes affected CHCs in a similar way to that seen in D. serrata. Population genomic analysis revealed that balancing selection likely operates at the DsFAR2-B locus in the wild. Our study provides insights into the genetic basis of sexual antagonism in nature and connects sexually varying antagonistic selection on phenotypes with balancing selection on genotypes that maintains molecular variation.

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.

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.

Cloning and reproductive regulation of a trypsin precursor gene in Adelphocorissuturalis

Adelphocoris suturalis is a major pest of cotton. Here, we identified a trypsin precursor gene (AsTryP) in A. suturali, which has an open reading frame region of 873 bp and belongs to the trypsin superfamily. The mRNA of the AsTryP gene was detectable in every life stage and different tissues of 8-day-old females, and the gene was highly expressed in fourth-instar nymphs and the thorax of 8-day-old females. Down-regulation of AsTryP by the injection of double-stranded RNA suppressed the ovarian development and female fertility. These results reveal that trypsin precursor is involved in the reproductive process of A. suturali, and may facilitate the development of new strategies for a better control of A. suturalis.

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.

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.

Integrating genomics and multivariate evolutionary quantitative genetics: a case study of constraints on sexual selection in Drosophila serrata

In evolutionary quantitative genetics, the genetic variance-covariance matrix, G, and the vector of directional selection gradients, β, are key parameters for predicting multivariate selection responses and genetic constraints. Historically, investigations of G and β have not overlapped with those dissecting the genetic basis of quantitative traits. Thus, it remains unknown whether these parameters reflect pleiotropic effects at individual loci. Here, we integrate multivariate genome-wide association study (GWAS) with G and β estimation in a well-studied system of multivariate constraint: sexual selection on male cuticular hydrocarbons (CHCs) in Drosophila serrata. In a panel of wild-derived re-sequenced lines, we augment genome-based restricted maximum likelihood to estimate G alongside multivariate single nucleotide polymorphism (SNP) effects, detecting 532 significant associations from 1 652 276 SNPs. Constraint was evident, with β lying in a direction of G with low evolvability. Interestingly, minor frequency alleles typically increased male CHC-attractiveness suggesting opposing natural selection on β. SNP effects were significantly misaligned with the major eigenvector of G, gmax, but well aligned to the second and third eigenvectors g2 and g3. We discuss potential factors leading to these varied results including multivariate stabilizing selection and mutational bias. Our framework may be useful as researchers increasingly access genomic methods to study multivariate selection responses in wild populations.

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.

Production of Long Chain Fatty Alcohols Found in Bumblebee Pheromones by Yarrowia lipolytica

Fatty alcohols (FA-OH) are aliphatic unbranched primary alcohols with a chain of four or more carbon atoms. Besides potential industrial applications, fatty alcohols have important biological functions as well. In nature, fatty alcohols are produced as a part of a mixture of pheromones in several insect species, such as moths, termites, bees, wasps, etc. In addition, FA-OHs have a potential for agricultural applications, for example, they may be used as a suitable substitute for commercial insecticides. The insecticides have several drawbacks associated with their preparation, and they exert a negative impact on the environment. Currently, pheromone components are prepared mainly through the catalytic hydrogenation of plant oils and petrochemicals, which is an unsustainable, ecologically unfriendly, and highly expensive process. The biotechnological production of the pheromone components using engineered microbial strains and through the expression of the enzymes participating in the biosynthesis of these components is a promising approach that ensures ecological sustenance as well. The present study was aimed at evaluating the production of FA-OHs in the oleaginous yeast, Yarrowia lipolytica, with different lengths of fatty-acyl chains by expressing the fatty acyl-CoA reductase (FAR) BlapFAR4 from B. lapidarius, producing C16:0-OH, C16:1Δ⁹-OH, and lower quantities of both C14:0-OH and C18:1Δ⁹-OH, and BlucFAR1 from B. lucorum, producing FA-OHs with a chain length of 18-26 carbon atoms, in this yeast. Among the different novel Y. lipolytica strains used in the present study, the best results were obtained with JMY7086, which carried several lipid metabolism modifications and expressed the BlucFAR1 gene under the control of a strong constitutive promoter 8UAS-pTEF. JMY7086 produced only saturated fatty alcohols with chain lengths from 18 to 24 carbon atoms. The highest titer and accumulation achieved were 166.6 mg/L and 15.6 mg/g DCW of fatty alcohols, respectively. Unlike JMY7086, the BlapFAR4-expressing strain JMY7090 produced only 16 carbon atom-long FA-OHs with a titer of 14.6 mg/L.

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.

3-Acetoxy-fatty acid isoprenyl esters from androconia of the ithomiine butterfly Ithomia salapia

Male ithomiine butterflies (Nymphalidae: Danainae) have hairpencils on the forewings (i.e., androconia) that disseminate semiochemicals during courtship. While most ithomiines are known to contain derivatives of pyrrolizidine alkaloids, dihydropyrrolizines, or γ-lactones in these androconia, here we report on a new class of fatty acid esters identified in two subspecies, Ithomia salapia aquinia and I. s. derasa. The major components were identified as isoprenyl (3-methyl-3-butenyl) (Z)-3-acetoxy-11-octadecenoate, isoprenyl (Z)-3-acetoxy-13-octadecenoate (12) and isoprenyl 3-acetoxyoctadecanoate (11) by GC/MS and GC/IR analyses, microderivatizations, and synthesis of representative compounds. The absolute configuration of 12 was determined to be R. The two subspecies differed not only in the composition of the ester bouquet, but also in the composition of more volatile androconial constituents. While some individuals of I. s. aquinia contained ithomiolide A (3), a pyrrolizidine alkaloid derived γ-lactone, I. s. derasa carried the sesquiterpene α-elemol (8) in the androconia. These differences might be important for the reproductive isolation of the two subspecies, in line with previously reported low gene exchange between the two species in regions where they co-occur. Furthermore, the occurrence of positional isomers of unsaturated fatty acid derivatives indicates activity of two different desaturases within these butterflies, Δ9 and Δ11, which has not been reported before in male Lepidoptera.

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.

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.

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.

Transcriptome Analysis of Female and Male Conopomorpha sinensis (Lepidoptera: Gracilariidae) Adults With a Focus on Hormone and Reproduction

Conopomorpha sinensis Bradley is the dominant borer pest of litchi and longan in the Asian-pacific area. Reduction or interference of reproduction and mating of adult moths is one of the most used strategies to control C. sinensis. Insect reproduction is a critical biological process closely related to endocrine control. Conopomorpha sinensis genome and transcriptome information is limited, hampering both our understanding of the molecular mechanisms underlying hormone activity and reproduction and the development of control strategies for this borer pest. To explore the sex differences in gene expression profiles influencing these biological processes, de novo transcriptomes were assembled from female and male adult C. sinensis specimens. This analysis yielded 184,422 unigenes with an average length of 903 bp and 405,961 transcripts after sequencing and assembly. About 45.06, 22.41, 19.53, 34.05, 35.82, 36.42, and 19.85% of the unigenes had significant matches in seven public databases. Subsequently, gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis revealed comprehensive information about the function of each gene and identified enriched categories and pathways that were associated with the 2,890 female-biased genes and 2,964 male-biased genes. In addition, we identified some important unigenes related to hormone activity and reproduction among the sex-differentially expressed genes (DEGs), including unigenes coding for ecdysone-induced protein 78C, juvenile hormone (JH)-regulated gene fatty acyl-CoA reductase, vitellogenin, etc. Our findings provide a more comprehensive portrait of the sex differences involved in the relationship of two important physiological features-hormone activity and reproduction in C. sinensis and members of the family Gracillariidae.

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.

Candidate genes involved in spiroacetal biosynthesis in the oriental fruit fly, Bactrocera dorsalis

Spiroacetals are widespread in nature as components of volatile semiochemical secretions from many insect species. The general pathway for spiroacetal biosynthesis in Bactrocera sp. is preliminarily established, but many genes involved in this pathway remain to be characterized. By analyzing transcriptomes of the rectal glands (RGs) from immature and mature females of the oriental fruit fly, Bactrocera dorsalis, a set of genes encoding two acetyl-CoA carboxylases (ACCs), two fatty acid synthases (FASs), eight desaturases (DESs), twelve fatty acyl-CoA reductases (FARs), seventy-two cytochrome P450 enzymes (CYPs), and twenty-three odorant binding proteins (OBPs) were identified. We investigated the expression of candidate genes in immature and mature stages based on the RNA-seq data and Real-time quantitative PCR. Expression profiling revealed that some of these genes were primarily expressed in female rectal glands among different tissues, and were up-regulated in mature females. Semi-quantitative RT-PCR assays were also adapted to examine tissue-specific expression of selected candidate genes. Additionally, their putative functions in spiroacetal synthesis and transportation are proposed. Our study provided large-scale sequence information for further functional studies on spiroacetal biosynthetic pathways.

FAR gene enables the brown planthopper to walk and jump on water in paddy field

Many insects can live on water and survive being caught in the rain. Current research has shown that insect cuticular hydrocarbons (CHC) confer desiccation resistance to maintain water balance. In this study, we identified a fatty acyl-CoA reductase gene (NlFAR) of the rice brown planthopper, Nilaparvata lugens that is essential for the production of CHCs, and found that NlFAR is essential for N. lugens to walk and jump on water when moving from one rice plant to another in paddy fields. NlFAR was mainly expressed in the integument at the beginning of each molt. Cuticular surface analysis by scanning electron microscopy and characterization of CHC extracts indicated that N. lugens with knockdown of NlFAR using RNA inference (RNAi) had a neater epicuticle layer and a significant decrease in CHC contents. Knockdown of NlFAR did not influence the desiccation resistance of N. lugens, but the dsNlFAR-treated insects were easily adhered and moistened by water droplets or their own secreted honeydew and unable to walk or jump on water. These results suggested that NlFAR is a crucial enzyme for CHC biosynthesis and cuticle waterproofing, but not for water retention of N. lugens, which may provide a potential strategy for pest management.

Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera

Fatty acyl reductases (FARs) are involved in the biosynthesis of fatty alcohols that serve a range of biological roles. Insects typically harbor numerous FAR gene family members. While some FARs are involved in pheromone biosynthesis, the biological significance of the large number of FARs in insect genomes remains unclear.

Using bumble bee (Bombini) FAR expression analysis and functional characterization, hymenopteran FAR gene tree reconstruction, and inspection of transposable elements (TEs) in the genomic environment of FARs, we uncovered a massive expansion of the FAR gene family in Hymenoptera, presumably facilitated by TEs. The expansion occurred in the common ancestor of bumble bees and stingless bees (Meliponini). We found that bumble bee FARs from the expanded FAR-A ortholog group contribute to the species-specific pheromone composition. Our results indicate that expansion and functional diversification of the FAR gene family played a key role in the evolution of pheromone communication in Hymenoptera.

Mountain pine beetle (Dendroctonus ponderosae) CYP4Gs convert long and short chain alcohols and aldehydes to hydrocarbons

Hydrocarbon biosynthesis in insects involves the elongation of fatty acyl-CoAs to very-long chain fatty acyl-CoAs that are then reduced and converted to hydrocarbon, with the last step involving the oxidative decarbonylation of an aldehyde to hydrocarbon and carbon dioxide. Cytochromes P450 in the 4G family decarbonylate aldehydes to hydrocarbon. All insect acyl-CoA reductases studied to date reduce fatty acyl-CoAs to alcohols. The results of the work reported herein demonstrate that CYP4G55 and CYP4G56 from the mountain pine beetle, Dendroctonus ponderosae, expressed as fusion proteins with house fly cytochrome P450 reductase (CPR), convert both long chain aldehydes and long chain alcohols to hydrocarbons. CYP4G55 and CYP4G56 appear to prefer primary alcohols to aldehydes as substrates. These data strongly suggest that hydrocarbon biosynthesis in insects occurs by the two-step reduction of very long chain fatty acyl-CoAs to alcohols, which are then oxidized to aldehydes and then oxidatively decarbonylated to hydrocarbon by CYP4G enzymes. In addition, both CYP4G55 and CYP4G56 fusion proteins convert C10 alcohols and aldehydes to hydrocarbons, including the conversion of (Z)-7-decenal, a putative intermediate in the exo-brevicomin pheromone biosynthetic pathway, to (Z)-3-nonene. These data demonstrate that the highly conserved CYP4G enzymes accept a broad range of carbon chain lengths, including C10 and C18, and have evolved to function in cuticular hydrocarbon biosynthesis and pheromone production.

BmSUC1 is essential for glycometabolism modulation in the silkworm, Bombyx mori

Sucrose is the most commonly transported sugar in plants and is easily assimilated by insects to fulfill the requirement of physiological metabolism. BmSuc1 is a novel animal β-fructofuranosidase (β-FFase, EC 3.2.1.26)-encoding gene that was firstly cloned and identified in silkworm, Bombyx mori. BmSUC1 was presumed to play an important role in the silkworm-mulberry enzymatic adaptation system by effectively hydrolyzing sucrose absorbed from mulberry leaves. However, this has not been proved with direct evidence thus far. In this study, we investigated sucrose hydrolysis activity in the larval midgut of B. mori by inhibition tests and found that sucrase activity mainly stemmed from β-FFase, not α-glucosidase. Next, we performed shRNA-mediated transgenic RNAi to analyze the growth characteristics of silkworm larvae and variations in glycometabolism in vivo in transgenic silkworms. The results showed that in the RNAi-BmSuc1 transgenic line, larval development was delayed, and their body size was markedly reduced. Finally, the activity of several disaccharidases alone in the midgut and the sugar distribution, total sugar and glycogen in the midgut, hemolymph and fat body were then determined and compared. Our results demonstrated that silencing BmSuc1 significantly reduced glucose and apparently activated maltase and trehalase in the midgut. Together with a clear decrease in both glycogen and trehalose in the fat body, we conclude that BmSUC1 acts as an essential sucrase by directly modulating the degree of sucrose hydrolysis in the silkworm larval midgut, and insufficient sugar storage in the fat body may be responsible for larval malnutrition and abnormal petite phenotypes.

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.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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

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

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.

A transgenic strategy for controlling plant bugs (Adelphocoris suturalis) through expression of double-stranded RNA homologous to fatty acyl-coenzyme A reductase in cotton

Plant bugs (Miridae species), which are sap-sucking insects, have emerged as major pests of cotton in China. Most Miridae species are not sensitive to commercial Bacillus thuringiensis (Bt) cotton, resulting in significant economic losses and an increased application of insecticide, which eventually may compromise the future of Bt cotton. We demonstrate that FATTY ACYL-COA REDUCTASE (AsFAR) plays an essential role in the reproduction of the bug Adelphocoris suturalis. Down-regulation of AsFAR expression by injection of double-stranded RNA suppresses ovarian development and female fertility, resulting in females producing few viable offspring. To determine the viability of an RNA interference approach to limit FAR expression and reproductive ability in A. suturalis, a dsRNA targeting the AsFAR gene (dsAsFAR) of A. suturalis was expressed in transgenic cotton plants. AsFAR transcription levels were significantly downregulated in A. suturalis feeding on the transgenic plants. In contained field trials, the transgenic cotton lines significantly suppressed the development of A. suturalis populations and were resistant to damage caused by plant bug infestation. These results suggest a new strategy for the management of plant bug pests of cotton.

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.

Epimerisation of chiral hydroxylactones by short-chain dehydrogenases/reductases accounts for sex pheromone evolution in Nasonia

Males of all species of the parasitic wasp genus Nasonia use (4R,5S)-5-hydroxy-4-decanolide (RS) as component of their sex pheromone while only N. vitripennis (Nv), employs additionally (4R,5R)-5-hydroxy-4-decanolide (RR). Three genes coding for the NAD⁺-dependent short-chain dehydrogenases/reductases (SDRs) NV10127, NV10128, and NV10129 are linked to the ability of Nv to produce RR. Here we show by assaying recombinant enzymes that SDRs from both Nv and N. giraulti (Ng), the latter a species with only RS in the pheromone, epimerise RS into RR and vice versa with (4R)-5-oxo-4-decanolide as an intermediate. Nv-derived SDR orthologues generally had higher epimerisation rates, which were also influenced by NAD⁺ availability. Semiquantitative protein analyses of the pheromone glands by tandem mass spectrometry revealed that NV10127 as well as NV10128 and/or NV10129 were more abundant in Nv compared to Ng. We conclude that the interplay of differential expression patterns and SDR epimerisation rates on the ancestral pheromone component RS accounts for the evolution of a novel pheromone phenotype in Nv.

Two fatty acyl reductases involved in moth pheromone biosynthesis

Fatty acyl reductases (FARs) constitute an evolutionarily conserved gene family found in all kingdoms of life. Members of the FAR gene family play diverse roles, including seed oil synthesis, insect pheromone biosynthesis, and mammalian wax biosynthesis. In insects, FAR genes dedicated to sex pheromone biosynthesis (pheromone-gland-specific fatty acyl reductase, pgFAR) form a unique clade that exhibits substantial modifications in gene structure and possesses unique specificity and selectivity for fatty acyl substrates. Highly selective and semi-selective ‘single pgFARs’ produce single and multicomponent pheromone signals in bombycid, pyralid, yponomeutid and noctuid moths. An intriguing question is how a ‘single reductase’ can direct the synthesis of several fatty alcohols of various chain lengths and isomeric forms. Here, we report two active pgFARs in the pheromone gland of Spodoptera, namely a semi-selective, C14:acyl-specific pgFAR and a highly selective, C16:acyl-specific pgFAR, and demonstrate that these pgFARs play a pivotal role in the formation of species-specific signals, a finding that is strongly supported by functional gene expression data. The study envisages a new area of research for disclosing evolutionary changes associated with C₁₄- and C₁₆-specific FARs in moth pheromone biosynthesis.

Molecular Characterization of Two Fatty Acyl-CoA Reductase Genes From Phenacoccus solenopsis (Hemiptera: Pseudococcidae)

Fatty acyl-CoA reductases (FARs) are key enzymes involved in fatty alcohol synthesis. Here, we cloned and characterized full-length cDNAs of two FAR genes from the cotton mealybug, Phenacoccus solenopsis. The results showed PsFAR I and PsFAR II cDNAs were 1,584 bp and 1,515 bp in length respectively. Both PsFAR I and PsFAR II were predicted to be located in the endoplasmic reticulum by Euk-mPLoc 2.0 approach. Both of them had a Rossmann folding region and a FAR_C region. Two conservative motifs were discovered in Rossmann folding region by sequence alignment including a NADPH combining motif, TGXXGG, and an active site motif, YXXXK. A phylogenetic tree made using MEGA 6.06 indicated that PsFAR I and PsFAR II were placed in two different branches. Gene expression analysis performed at different developmental stages showed that the expression of PsFar I is significantly higher than that of PsFar II in first and second instar nymphs and in male adults. Spirotetramat treatment at 125 mg/liter significantly increased the expression of PsFar I in third instar nymphs, but there was no effect in the expression of PsFar II Our results indicated these two FAR genes showed different expression patterns during insect development and after pesticide treatment, suggesting they play different roles in insect development and detoxification against pesticides.

The Yeast ATF1 Acetyltransferase Efficiently Acetylates Insect Pheromone Alcohols: Implications for the Biological Production of Moth Pheromones

Many moth pheromones are composed of mixtures of acetates of long-chain (≥10 carbon) fatty alcohols. Moth pheromone precursors such as fatty acids and fatty alcohols can be produced in yeast by the heterologous expression of genes involved in insect pheromone production. Acetyltransferases that subsequently catalyze the formation of acetates by transfer of the acetate unit from acetyl-CoA to a fatty alcohol have been postulated in pheromone biosynthesis. However, so far no fatty alcohol acetyltransferases responsible for the production of straight chain alkyl acetate pheromone components in insects have been identified. In search for a non-insect acetyltransferase alternative, we expressed a plant-derived diacylglycerol acetyltransferase (EaDAcT) (EC 2.3.1.20) cloned from the seed of the burning bush (Euonymus alatus) in a yeast system. EaDAcT transformed various fatty alcohol insect pheromone precursors into acetates but we also found high background acetylation activities. Only one enzyme in yeast was shown to be responsible for the majority of that background activity, the acetyltransferase ATF1 (EC 2.3.1.84). We further investigated the usefulness of ATF1 for the conversion of moth pheromone alcohols into acetates in comparison with EaDAcT. Overexpression of ATF1 revealed that it was capable of acetylating these fatty alcohols with chain lengths from 10 to 18 carbons with up to 27- and 10-fold higher in vivo and in vitro efficiency, respectively, compared to EaDAcT. The ATF1 enzyme thus has the potential to serve as the missing enzyme in the reconstruction of the biosynthetic pathway of insect acetate pheromones from precursor fatty acids in yeast.

The Genetic Basis of Pheromone Evolution in Moths

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

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

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

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.

Identification and characterization of a fatty acyl reductase from a Spodoptera littoralis female gland involved in pheromone biosynthesis

Fatty acyl-CoA reductases (FARs), the enzymes that catalyse reduction of a fatty acyl-CoA to the corresponding alcohol in insect pheromone biosynthesis, are postulated to play an important role in determining the proportion of each component in the pheromone blend. For the first time, we have isolated and characterized from the Egyptian cotton leaf worm Spodoptera littoralis (Lepidoptera: Noctuidae) a FAR cDNA (Slit-FAR1), which appeared to be expressed only in the pheromone gland and was undetectable in other female tissues, such as fat body, ovaries, wings, legs or thorax. The encoded protein has been successfully expressed in a recombinant system, and the recombinant enzyme is able to produce the intermediate fatty acid alcohols of the pheromone biosynthesis of S. littoralis from the corresponding acyl-CoA precursors. The kinetic variables Km and Vmax, which have been calculated for each acyl-CoA pheromone precursor, suggest that in S. littoralis pheromone biosynthesis other biosynthetic enzymes (e.g. desaturases, acetyl transferase) should also contribute to the final ratio of components of the pheromone blend. In a phylogenetic analysis, Slit-FAR1 appeared grouped in a cluster of other FARs involved in the pheromone biosynthesis of other insects, with little or non-specificity for the natural pheromone precursors.