The role of desaturases in the biosynthesis of marking pheromones in bumblebee males

The influence of first desaturase subfamily genes on fatty acid synthesis, desiccation tolerance and inter-caste nutrient transfer in the termite Coptotermes formosanus

Desaturase enzymes play an essential role in the biosynthesis of unsaturated fatty acids (UFAs). In this study, we identified seven "first desaturase" subfamily genes (Cfor-desatA1, Cfor-desatA2-a, Cfor-desatA2-b, Cfor-desatB-a, Cfor-desatB-b, Cfor-desatD and Cfor-desatE) from the Formosan subterranean termite Coptotermes formosanus. These desaturases were highly expressed in the cuticle and fat body of C. formosanus. Inhibition of either the Cfor-desatA2-a or Cfor-desatA2-b gene resulted in a significant decrease in the contents of fatty acids (C16:0, C18:0, C18:1 and C18:2) in worker castes. Moreover, we observed that inhibition of most of desaturase genes identified in this study had a negative impact on the survival rate and desiccation tolerance of workers. Interestingly, when normal soldiers were reared together with dsCfor-desatA2-b-treated workers, they exhibited higher mortality, suggesting that desaturase had an impact on trophallaxis among C. formosanus castes. Our findings shed light on the novel roles of desaturase family genes in the eusocial termite C. formosanus.

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.

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.

The Desaturase Gene Nlug-desatA2 Regulates the Performance of the Brown Planthopper Nilaparvata lugens and Its Relationship with Rice

Insect desaturases are known to play an important role in chemical communication between individuals. However, their roles in insect growth, development and fecundity, and in regulating interactions of insects with plants, remain largely unknown. In this study, we explored the functions of Nlug-desatA2, a desaturase gene of the brown planthopper (BPH), Nilaparvata lugens (Stål). The RNA interference-based knockdown of Nlug-desatA2 decreased the ratio of monounsaturated fatty acids to saturated fatty acids, and the level of fatty acids and triglycerides in BPH. Nlug-desatA2-knockdown also reduced the food intake, body mass and fecundity of female BPH adults, and led to abdomen atrophy and ovarian agenesis. Nlug-desatA2-knockdown suppressed the transcription of TOR (target of rapamycin), Lpp (Lipophorin) and AKHR (adipokinetic hormone receptor) in female adults. Moreover, the corrected survival rate of BPH with Nlug-desatA2-knockdown fed an artificial diet was higher than the survival rate of those fed on rice plants. Higher levels of salicylic acid in rice infested by Nlug-desatA2-knockdown female BPH adults than in rice infested by control BPH may be the reason. These findings demonstrate that Nlug-desatA2 has an essential role in lipid metabolism and is involved in the food intake, survival, development and fecundity of BPH. In addition, this gene is likely involved in regulating the responses of rice to BPH infestation.

Desaturase specificity is controlled by the physicochemical properties of a single amino acid residue in the substrate binding tunnel

Membrane fatty acyl desaturases (mFAD) are ubiquitous enzymes in eukaryotes. They introduce double bonds into fatty acids (FAs), producing structurally diverse unsaturated FAs which serve as membrane lipid components or precursors of signaling molecules. The mechanisms controlling enzymatic specificity and selectivity of desaturation are, however, poorly understood. We found that the physicochemical properties, particularly side chain volume, of a single amino acid (aa) residue in insect mFADs (Lepidoptera: Bombyx mori and Manduca sexta) control the desaturation products. Molecular dynamics simulations of systems comprising wild-type or mutant mFADs with fatty acyl-CoA substrates revealed that the single aa substitution likely directs the outcome of the desaturation reaction by modulating the distance between substrate fatty acyl carbon atoms and active center metal ions. These findings, as well as our methodology combining mFAD mutational screening with molecular dynamics simulations, will facilitate prediction of desaturation products and facilitate engineering of mFADs for biotechnological applications.

The Desaturase Gene Family is Crucially Required for Fatty Acid Metabolism and Survival of the Brown Planthopper, Nilaparvata lugens

Desaturases are essentially required for unsaturated fatty acid (UFA) biosynthesis. We identified 10 genes encoding putative desaturases in the transcriptome database of the brown planthopper (BPH), Nilaparvata lugens. These include eight First Desaturase family genes, one cytochrome b5 fused desaturase gene (Nlug-Cytb5r) and one Sphingolipid Desaturase gene (Nlug-ifc). Transcript level profiling revealed significant variation in the expression patterns of these genes across tissues and developmental stages, which occur in a gene-specific manner. Interestingly, their expression was also modulated by the insect food source: the mRNA levels of Nlug-desatC and Nlug-Cytb5r were down-regulated, but the expression level of Nlug-desatA1-b and Nlug-desatA1-c were elevated in the BPH fed on the resistant rice variety Babawee as compared to the non-resistant variety Taichun Native 1 (TN1). Silencing Nlug-desatA1-b, Nlug-desatA1-c, or Nlug-Ifc reduced fatty acid composition and abundance in female BPH 1-d-old-adults compared to controls. Whereas, single knockdown of all ten desaturase genes significantly increased mortality of BPH nymphs compared with controls. Of the ten desaturase genes, knockdown of Nlug-desatA1-b and Nlug-desatA2 caused the highest mortality in BPH (91% and 97%, respectively). Our findings offer a base for expression and functional characterization of newly identified desaturase genes in BPH, and may contribute to RNA interference-based pest management strategies.

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.

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.

Fatty acid composition in native bees: Associations with thermal and feeding ecology

Fatty acid (FA) composition of lipids plays a crucial role in the functioning of lipid-containing structures in organisms and may be affected by the temperature an organism experiences, as well as its diet. We compared FA composition among four bee genera: Andrena, Bombus, Megachile, and Osmia which differ in their thermal ecology and diet. Fatty acid methyl esters (FAME) were prepared by direct transesterification with KOH and analyzed using gas-liquid chromatography with a flame ionization detector. Sixteen total FAs ranging in chain length from eight to 22 carbon atoms were identified. Linear discriminant analysis separated the bees based on their FA composition. Andrena was characterized by relatively high concentrations of polyunsaturated FAs, Bombus by high monounsaturated FAs and Megachilids (Megachile and Osmia) by relatively high amounts of saturated FAs. These differences in FA composition may in part be explained by variation in the diets of these bees. Because tongue (proboscis) length may be used as a proxy for the types of flowers bees may visit for nectar and pollen, we compared FA composition among Bombus that differed in proboscis length (but have similar thermal ecology). A clear separation in FA composition within Bombus with varying proboscis lengths was found using linear discriminant analysis. Further, comparing the relationship between each genus by cluster analysis revealed aggregations by genus that were not completely separated, suggesting potential overlap in dietary acquisition of FAs.

Exploring complex pheromone biosynthetic processes in the bumblebee male labial gland by RNA sequencing

Male marking pheromones (MPs) are used by the majority of bumblebee species (Hymenoptera: Apidae), including a commercially important greenhouse pollinator, the buff-tailed bumblebee (Bombus terrestris), to attract conspecific females. MP biosynthetic processes in the cephalic part of the bumblebee male labial gland (LG) are of extraordinary complexity, involving enzymes of fatty acid and isoprenoid biosynthesis, which jointly produce more than 50 compounds. We employed a differential transcriptomic approach to identify candidate genes involved in MP biosynthesis by sequencing Bombus terrestris LG and fat body (FB) transcriptomes. We identified 12 454 abundantly expressed gene products (reads per kilobase of exon model per million mapped reads value > 1) that had significant hits in the GenBank nonredundant database. Of these, 876 were upregulated in the LG (> 4-fold difference). We identified more than 140 candidate genes potentially involved in MP biosynthesis, including esterases, fatty acid reductases, lipases, enzymes involved in limited fatty acid chain shortening, neuropeptide receptors and enzymes involved in biosynthesis of triacylglycerols, isoprenoids and fatty acids. For selected candidates, we confirmed their abundant expression in LG using quantitative real-time reverse transcription-PCR (qRT-PCR). Our study shows that the Bombus terrestris LG transcriptome reflects both fatty acid and isoprenoid MP biosynthetic processes and identifies rational gene targets for future studies to disentangle the molecular basis of MP biosynthesis. Additionally, LG and FB transcriptomes enrich the available transcriptomic resources for Bombus terrestris.

Regulation of Isoprenoid Pheromone Biosynthesis in Bumblebee Males

Males of the closely related species Bombus terrestris and Bombus lucorum attract conspecific females by completely different marking pheromones. MP of B. terrestris and B. lucorum pheromones contain mainly isoprenoid (ISP) compounds and fatty acid derivatives, respectively. Here, we studied the regulation of ISP biosynthesis in both bumblebees. RNA-seq and qRT-PCR analyses indicated that acetoacetyl-CoA thiolase (AACT), 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), and farnesyl diphosphate synthase (FPPS) transcripts are abundant in the B. terrestris labial gland. Maximal abundance of these transcripts correlated well with AACT enzymatic activity detected in the LG extracts. In contrast, transcript abundances of AACT, HMGR, and FPPS in B. lucorum were low, and AACT activity was not detected in LGs. These results suggest that transcriptional regulation plays a key role in the control of ISP biosynthetic gene expression and ISP pheromone biosynthesis in bumblebee males.

Biosynthetic Studies of the Male Marking Pheromone in Bumblebees by Using Labelled Fatty Acids and Two-Dimensional Gas Chromatography with Mass Detection

Two-dimensional comprehensive gas chromatography (GC×GC) coupled with mass detection was used as a tool for biosynthetic studies of bumblebee pheromones. Prior to biosynthetic experiments, the chromatographic behaviour of isotopically modified esters in the GC×GC system as well as their behaviour in mass detection was studied. The male marking pheromones of Bombus lucorum, Bombus lapidarius and Bombus terrestris were investigated. Main pheromonal components are ethyl tetradec-9-enoate (53 %) and ethyl dodecanoate (6 %) in B. lucorum, hexadec-9-en-1-ol (52 %) and hexadecan-1-ol (31 %) in B. lapidarius, and 2,3-dihydrofarnesol (58 %) and ethyl dodecanoate (15 %) in B. terrestris. The research strategy was based on 1) in vivo incubation of isotopically (² H, ¹³ C) modified fatty acids (FAs) and analysis of their metabolites and 2) feeding experiments with ² H- and ¹³ C-labelled FAs mixed with food. It was observed that labelled FAs were modified into the most abundant aliphatic compounds present in labial gland secretions. In feeding experiments, the labelled FAs were transformed into pheromone components. Transport of the FA precursors from the fat body through haemolymph was confirmed. The results show that FAs, stored in the form of triacylglycerols in the fat body, are likely to participate in the biosynthesis of some aliphatic pheromone components.

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.

Exoskeleton formation in Apis mellifera: cuticular hydrocarbons profiles and expression of desaturase and elongase genes during pupal and adult development

Cuticular hydrocarbons (CHCs) are abundant in the superficial cuticular layer (envelope) of insects where they play roles as structural, anti-desiccation and semiochemical compounds. Many studies have investigated the CHC composition in the adult insects. However, studies on the profiles of these compounds during cuticle formation and differentiation are scarce and restrict to specific stages of a few insect species. We characterized the CHCs developmental profiles in the honeybee workers during an entire molting cycle (from pupal-to-adult ecdyses) and in mature adults (forager bees). Gas chromatography/mass spectrometry (GC/MS) analysis revealed remarkable differences in the relative quantities of CHCs, thus discriminating pupae, developing and newly-ecdysed adults, and foragers from each other. In parallel, the honeybee genome database was searched for predicted gene models using known amino acid sequences of insect enzymes catalyzing lipid desaturation (desaturases) or elongation (elongases) as queries in BLASTP analysis. The expression levels of six desaturase genes and ten elongase genes potentially involved in CHC biosynthesis were determined by reverse transcription and real time polymerase chain reaction (RT-qPCR) in the developing integument (cuticle and subjacent epidermis). Aiming to predict roles for these genes in CHC biosynthesis, the developmental profiles of CHCs and desaturase/elongase transcript levels were evaluated using Spearman correlation coefficient. This analysis pointed to differential roles for these gene products in the biosynthesis of certain CHC classes. Based on the assumption that homologous proteins may share a similar function, phylogenetic trees were reconstructed as an additional strategy to predict functions and evolutionary relationships of the honeybee desaturases and elongases. Together, these approaches highlighted the molecular complexity underlying the formation of the lesser known layer of the cuticular exoskeleton, the envelope.

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.

Δ12-Fatty acid desaturase from Candida parapsilosis is a multifunctional desaturase producing a range of polyunsaturated and hydroxylated fatty acids

Numerous Δ12-, Δ15- and multifunctional membrane fatty acid desaturases (FADs) have been identified in fungi, revealing great variability in the enzymatic specificities of FADs involved in biosynthesis of polyunsaturated fatty acids (PUFAs). Here, we report gene isolation and characterization of novel Δ12/Δ15- and Δ15-FADs named CpFad2 and CpFad3, respectively, from the opportunistic pathogenic yeast Candida parapsilosis. Overexpression of CpFad3 in Saccharomyces cerevisiae strains supplemented with linoleic acid (Δ9,Δ12-18:2) and hexadecadienoic acid (Δ9,Δ12-16:2) leads to accumulation of Δ15-PUFAs, i.e., α-linolenic acid (Δ9,Δ12,Δ15-18:3) and hexadecatrienoic acid with an unusual terminal double bond (Δ9,Δ12,Δ15-16:3). CpFad2 produces a range of Δ12- and Δ15-PUFAs. The major products of CpFad2 are linoleic and hexadecadienoic acid (Δ9,Δ12-16:2), accompanied by α-linolenic acid and hexadecatrienoic acid (Δ9,Δ12,Δ15-16:3). Using GC/MS analysis of trimethylsilyl derivatives, we identified ricinoleic acid (12-hydroxy-9-octadecenoic acid) as an additional product of CpFad2. These results demonstrate that CpFAD2 is a multifunctional FAD and indicate that detailed analysis of fatty acid derivatives might uncover a range of enzymatic selectivities in other Δ12-FADs from budding yeasts (Ascomycota: Saccharomycotina).

Characterization of neutral lipase BT-1 isolated from the labial gland of Bombus terrestris males

Background

In addition to their general role in the hydrolysis of storage lipids, bumblebee lipases can participate in the biosynthesis of fatty acids that serve as precursors of pheromones used for sexual communication.

Results

We studied the temporal dynamics of lipolytic activity in crude extracts from the cephalic part of Bombus terrestris labial glands. Extracts from 3-day-old males displayed the highest lipolytic activity. The highest lipase gene expression level was observed in freshly emerged bumblebees, and both gene expression and lipase activity were lower in bumblebees older than 3 days. Lipase was purified from labial glands, further characterized and named as BT-1. The B. terrestris orthologue shares 88% sequence identity with B. impatiens lipase HA. The molecular weight of B. terrestris lipase BT-1 was approximately 30 kDa, the pH optimum was 8.3, and the temperature optimum was 50°C. Lipase BT-1 showed a notable preference for C8-C10 p-nitrophenyl esters, with the highest activity toward p-nitrophenyl caprylate (C8). The Michaelis constant (K_m) and maximum reaction rate (V_max) for p-nitrophenyl laurate hydrolysis were K_m = 0.0011 mM and V_max = 0.15 U/mg.

Conclusion

This is the first report describing neutral lipase from the labial gland of B. terrestris. Our findings help increase understanding of its possible function in the labial gland.