Identification of lipases involved in PBAN stimulated pheromone production in Bombyx mori using the DGE and RNAi approaches

Lipases and carboxylesterases affect moth sex pheromone compounds involved in interspecific mate recognition

Moth sex pheromones are a classical model for studying sexual selection. Females typically produce a species-specific pheromone blend that attracts males. Revealing the enzymes involved in the interspecific variation in blend composition is key for understanding the evolution of these sexual communication systems. The nature of the enzymes involved in the variation of acetate esters, which are prominent compounds in moth pheromone blends, remains unclear. We identify enzymes involved in acetate degradation using two closely related moth species: Heliothis (Chloridea) subflexa and H. (C.) virescens, which have different quantities of acetate esters in their sex pheromone. Through comparative transcriptomic analyses and CRISPR/Cas9 knockouts, we show that two lipases and two esterases from H. virescens reduce the levels of pheromone acetate esters when expressed in H. subflexa females. Together, our results show that lipases and carboxylesterases are involved in tuning Lepidoptera pheromones composition.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

ATP Synthase Subunit α from Helicoverpa armigera Acts as a Receptor of Bacillus thuringiensis Cry1Ac and Synergizes Cry1Ac Toxicity

Insect resistance to Bacillus thuringiensis (Bt) toxins has led to an urgent need to explore the insecticidal mechanisms of Bt. Previous studies indicated that Helicoverpa armigera ATP synthase subunit α (HaATPs-α) is involved in Cry1Ac resistance. In this study, a real-time quantitative polymerase chain reaction (RT-PCR) confirmed that HaATPs-α expression was significantly reduced in the Cry1Ac-resistant strain (BtR). Cry1Ac feeding induced the downregulated expression of HaATPs-α in the susceptible strain, but not in the BtR strain. Furthermore, the interaction between HaATPs-α and Cry1Ac was verified by ligand blotting and homologous competition experiments. The in vitro gain and loss of function analyses showed HaATPs-α involved in Cry1Ac toxicity by expressing endogenous HaATPs-α and HaATPs-α double-stranded RNAs in Sf9 and midgut cells, respectively. Importantly, purified HaATPs-α synergized Cry1Ac toxicity to H. armigera larvae. These findings provide the first evidence that HaATPs-α is a potential receptor of Cry1Ac, it shows downregulated participation in Cry1Ac resistance, and it exhibits higher enhancement of Cry1Ac toxicity to H. armigera larvae.

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.

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

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

Pyruvate Kinase Is Required for Sex Pheromone Biosynthesis in Helicoverpa armigera

Pyruvate kinase (PYK) is a speed-limited enzyme of glycolysis that catalyzes the formation of pyruvate, and plays an important role in acetyl-CoA synthesis. The acetyl-CoA is the precursor of sex pheromone biosynthesis in Helicoverpa armigera. However, the role of PYK in sex pheromone biosynthesis remains elusive. Here, PYK in H. armigera (HaPYK) was found to be highly expressed in the pheromone glands (PGs). The developmental expression profile of HaPYK was consistent with the fluctuation of sex pheromone release. Function analysis revealed that the knockdown of HaPYK led to a decrease in the levels of pyruvic acid and acetyl-CoA in PGs, which in turn caused a significant decrease in cis-11-hexadecenal (Z11-16: Ald) production, female capability to attract males, and mating frequency. Further study demonstrated that sugar feeding (5% sugar) increased the transcription and enzyme activity of HaPYK, thereby facilitating sex pheromone biosynthesis. Moreover, pheromone biosynthesis activating neuropeptide (PBAN) upregulated HaPYK activity through protein kinase C (PKC), as shown by PKC-specific inhibitor analysis. Altogether, our results revealed that PBAN activated HaPYK by Ca²⁺/PKC, thereby regulating the synthesis of pyruvate and subsequent acetyl-CoA, ensuring the supply of sex pheromone precursor, and finally facilitating sex pheromone biosynthesis and mating behavior.

Supplemental Sugar Is Required for Sex Pheromone Biosynthesis in Mythimna separata

Supplemental nutrients of adult moths maximize moth fitness and contribute to the pollination of many plants. Previous reports have revealed that sugar feeding promotes to sex pheromone biosynthesis by increasing the haemolymph trehalose concentration in mating moths. Here, Mythimna separata adults were employed as a model to investigate the effect of sugar feeding on sex pheromone biosynthesis. Results showed that in virgin females, sugar feeding markedly increased the concentrations of trehalose, pyruvic acid, and acyl-CoA in pheromone glands (PGs), which in turn led to an increase in sex pheromone titer, female ability to attract males and successfully mating frequency in sugar-fed females. Consistently, sugar-fed females laid more eggs than water-fed females. Furthermore, the refeeding of starved females also caused significantly increase in the concentrations of trehalose, pyruvic acid, and acyl-CoA in PGs, thus facilitating a significant increase in sex pheromone production. Most importantly, RNAi-mediated knockdown of trehalase (leading to PG starvation) resulted in an increase in trehalose content, and decrease in the concentrations of pyruvic acid, and acyl-CoA in PGs, which in turn led to a decrease of sex pheromone titer, female ability to attract males and successful mating efficacy. Altogether, results revealed a mechanism by which sugar feeding contributed to trehalose utilization in PGs, promoted to significantly increased sex pheromone precursor by increasing the concentrations of pyruvic acid and acyl-CoA, and facilitated to sex pheromone biosynthesis and successful mating.

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.

Comprehensive History of CSP Genes: Evolution, Phylogenetic Distribution and Functions

In this review we present the developmental, histological, evolutionary and functional properties of insect chemosensory proteins (CSPs) in insect species. CSPs are small globular proteins folded like a prism and notoriously known for their complex and arguably obscure function(s), particularly in pheromone olfaction. Here, we focus on direct functional consequences on protein function depending on duplication, expression and RNA editing. The result of our analysis is important for understanding the significance of RNA-editing on functionality of CSP genes, particularly in the brain tissue.

Comparative transcriptomics of the pheromone glands provides new insights into the differentiation of sex pheromone between two host populations of Chilo suppressalis

Reproductive isolation between different host populations is often based on intraspecific sex pheromone differences. The mechanisms underlying these differences have not been thoroughly elucidated to date. Previous studies suggested that Chilo suppressalis has differentiated into rice and water-oat host populations, and these two populations manifest clear differences in sex pheromone titer and mating rhythm. Hence, this moth is an ideal model to investigate the endogenous mechanisms of intraspecific reproductive isolation. Here, we identified a series of putative genes associated with sex pheromone biosynthesis based on the C. suppressalis pheromone gland transcriptome data. Transcripts of most genes were at higher level in the rice population. Then we obtained 11 pivotal differentially expressed genes (DEGs). The expression levels of these DEGs exhibited a distinct increase in the rice population. Moreover, we also observed the expression rhythm of these DEGs is discrepant between two host populations. Our study offers a new understanding to elucidate the mechanisms of intraspecific reproductive isolation.

Calcineurin is required for male sex pheromone biosynthesis and female acceptance

Lepidoptera sex pheromone biosynthesis is regulated by pheromone biosynthesis activating neuropeptide (PBAN). PBAN regulates not only female sex pheromone biosynthesis but also male sex pheromone biosynthesis. Previous research has confirmed that PBAN regulates sex pheromone biosynthesis using Ca²⁺ as a secondary messenger in all examined species to date. However, the downstream signal of Ca²⁺ has remained elusive. In the present study, calcineurin A (CNA), a downstream signal of Ca²⁺ , was discovered in Helicoverpa armigera male hairpencil and named HaCNA. Sequence analysis demonstrated that the open reading frame of HaCNA contains 1488 nucleotides encoding 495 amino acid residues. A homology search revealed that HaCNA shares a high amino acid identity with the CNA of other insects. Developmental and spatial expression analyses revealed that the mRNA levels of HaCNA peaked at 24 h after emergence and that HaCNA expression was ubiquitous in all examined tissues. Activity analysis revealed that PBAN activates HaCNA, and a Ca²⁺ inhibitor, Lacl₃ , attenuated the effect of PBAN by decreasing HaCNA activity. Pharmacological inhibitor and RNA interference-mediated knockdown assays revealed that both activity inhibition and decreased mRNA levels of HaCNA led to a significant decrease in the production of the male sex pheromone components [octadecanol and (Z)-11 hexadecanol)] and in the efficacy of female mating acceptance. Our results demonstrate that HaCNA acts as downstream signal of PBAN/Ca²⁺ and plays an important role in PBAN-induced male sex pheromone biosynthesis and female mating acceptance.

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.

Regulation of pheromone biosynthesis in moths

Female moths release sex pheromones for attracting males from a distance. Most moths are nocturnal so there is a periodicity to the release of sex pheromone. The temporal release of sex pheromone in most moths is regulated by calling behavior and by the biosynthesis of sex pheromone. In most moths, biosynthesis occurs in the pheromone gland and is controlled by the neuropeptide PBAN (pheromone biosynthesis activating neuropeptide). PBAN is produced in the subesophageal ganglion and released into circulation where it travels to the pheromone gland to activate pheromone biosynthesis. The G-protein coupled receptor that binds PBAN has been identified as well as aspects of signal transduction to activate the biosynthetic pathway. This review will highlight recent advances in the study of regulation of pheromone biosynthesis in moths.

Calcineurin-mediated Dephosphorylation of Acetyl-coA Carboxylase is Required for Pheromone Biosynthesis Activating Neuropeptide (PBAN)-induced Sex Pheromone Biosynthesis in Helicoverpa armigera

Chemical signaling plays a critical role in the behavior and physiology of many animals. Female insects, as many other animals, release sex pheromones to attract males for mating. The evolutionary and ecological success of insects therefore hinges on their ability to precisely mediate (including initiation and termination) pheromone biosynthesis. Pheromone biosynthesis activating neuropeptide (PBAN) acts directly on pheromone glands to regulate sex pheromone production using Ca2+ and cyclic-AMP as secondary messengers in the majority of species. However, the molecular mechanism downstream of the secondary messengers has not yet been elucidated in heliothine species. The present study shows that calcineurin, protein kinase A (PKA) and acetyl-coA carboxylase (ACC) are key components involved in PBAN-induced sex pheromone biosynthesis in Helicoverpa armigera using PBAN-dependent phosphoproteomics in combination with transcriptomics. RNAi-mediated knockdown and inhibitor assay demonstrated that calcineurin A is required for PBAN-induced ACC activation and sex pheromone production. Calcineurin-dependent phosphoproteomics and in vitro calcineurin phosphorylation assay further revealed that calcineurin regulated ACC activity by dephosphorylating ser84 and ser92. In addition, PKA-dependent phosphoproteomics and activity analysis revealed that PKA reduces the activity of AMP-activated protein kinase (AMPK), a negative regulator of ACC by phosphorylating the conserved ser92. Taken together, our findings indicate that calcineurin acts as the downstream signal of PBAN/G-protein receptor/Ca2+ to activate ACC through dephosphorylation while inactivating AMPK via PKA to reduce ACC phosphorylation, thus facilitating calcineurin activation of ACC.

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.

Sex pheromone in the moth Heliothis virescens is produced as a mixture of two pools: de novo and via precursor storage in glycerolipids

Most species of moths use a female-produced volatile sex pheromone, typically produced via de novo fatty acid synthesis in a specialized gland, for communication among mates. While de novo biosynthesis of pheromone (DNP) is rapid, suggesting transient precursor acids, substantial amounts of pheromone precursor (and other) acids are stored, predominantly in triacylglycerols in the pheromone gland. Whether these stored acids are converted to pheromone later or not has been the subject of some debate. Using a tracer/tracee approach, in which we fed female Heliothis virescens U-¹³C-glucose, we were able to distinguish two pools of pheromone, in which precursors were temporally separated (after and before feeding on labeled glucose): DNP synthesized from a mixed tracer/tracee acetyl CoA pool after feeding, and pheromone made from precursor acids primarily synthesized before feeding, which we call recycled precursor fat pheromone (RPP). DNP titer varied from high (during scotophase) to low (photophase) and with presence/absence of pheromone biosynthesis activating neuropeptide (PBAN), in accord with native pheromone titer previously observed. By contrast, RPP was constant throughout the photoperiod and did not change with PBAN presence/absence. The amount of RPP (6.3-10.3 ng/female) was typically much lower than that of DNP, especially during the scotophase (peak DNP, 105 ng/female). We propose an integral role for stored fats in pheromone biosynthesis, in which they are hydrolyzed and re-esterified throughout the photoperiod, with a small proportion of liberated precursor acyl CoAs being converted to pheromone. During the sexually active period, release of PBAN results in increased flux of glucose (from trehalose) and hydrolyzed acids entering the mitochondria, producing acetyl CoA precursor for de novo fat and pheromone biosynthesis.

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.

Multicopper oxidase-1 is required for iron homeostasis in Malpighian tubules of Helicoverpa armigera

Multicopper oxidases (MCOs) are enzymes that contain 10 conserved histidine residues and 1 cysteine residue. MCO1 has been extensively investigated in the midgut because this MCO is implicated in ascorbate oxidation, iron homeostasis and immune responses. However, information regarding the action of MCO1 in Malpighian tubules is limited. In this study, Helicoverpa armigera was used as a model to investigate the function of MCO1 in Malpighian tubules. Sequence analysis results revealed that HaMCO1 exhibits typical MCO characteristics, with 10 histidine and 1 cysteine residues for copper ion binding. HaMCO1 was also found to be highly abundant in Malpighian tubules. Temporal expression patterns indicated that HaMCO1 is mainly expressed during larval molting stages. Hormone treatments [the molting hormone 20-hydroxyecdysone (20E) and juvenile hormone (JH)] revealed that 20E inhibits HaMCO1 transcript expression via its heterodimer receptor, which consists of ecdysone receptor (EcR) and ultraspiracle (USP), and that JH counteracts the action of 20E to activate HaMCO1 transcript expression via its intracellular receptor methoprene-tolerant (Met). HaMCO1 knockdown caused a significant decrease in iron accumulation and also significantly reduced transferrin and ferritin transcript expression. Therefore, HaMCO1 is coordinately regulated by 20E and JH and is required for iron homeostasis in Malpighian tubules.

Transcriptome analysis of sexually dimorphic Chinese white wax scale insects reveals key differences in developmental programs and transcription factor expression

The Chinese white wax scale insect, Ericerus pela, represents one of the most dramatic examples of sexual dimorphism in any insect species. In this study, we showed that although E. pela males display complete metamorphosis similar to holometabolous insects, the species forms the sister group to Acyrthosiphon pisum and cluster with hemimetabolous insects. The gene expression profile and Gene Ontology (GO) analyses revealed that the two sexes engaged in distinct developmental programs. In particular, female development appeared to prioritize the expression of genes related to cellular, metabolic, and developmental processes and to anatomical structure formation in nymphs. By contrast, male nymphal development is characterized by the significant down-regulation of genes involved in chitin, the respiratory system, and neurons. The wing and appendage morphogenesis, anatomical and tissue structure morphogenesis programs activated after male nymphal development. Transcription factors (that convey juvenile hormone or ecdysone signals, and Hox genes) and DNA methyltransferase were also differentially expressed between females and males. These results may indicate the roles that these differentially expressed genes play in regulating sexual dimorphism through orchestrating complex genetic programs. This differential expression was particularly prominent for processes linked to female development and wing development in males.

Glycerol-3-phosphate O-acyltransferase is required for PBAN-induced sex pheromone biosynthesis in Bombyx mori

Female moths employ their own pheromone blends as a communicational medium in mating behavior. The biosynthesis and release of sex pheromone in female moths are regulated by pheromone biosynthesis activating neuropeptide (PBAN) and the corresponding action of PBAN has been well elucidated in Bombyx mori. However, very little is known about the molecular mechanism regarding the biosynthesis of sex pheromone precursor. In this study, quantitative proteomics was utilized to comprehensively elucidate the expression dynamics of pheromone glands (PGs) during development. Proteomic analysis revealed a serial of differentially expressed sex pheromone biosynthesis-associated proteins at the different time points of B. mori development. Most interestingly B. mori glycerol-3-phosphate O-acyltransferase (BmGPAT) was found to be expressed during the key periods of sex pheromone biosynthesis. RNAi knockdown of BmGPAT confirmed the important function of this protein in the biosynthesis of sex pheromone precursor, triacylglcerol (TAG), and subsequently PBAN-induced production of sex pheromone, bombykol. Behavioral analysis showed that RNAi knockdown of GPAT significantly impaired the ability of females to attract males. Our findings indicate that GPAT acts to regulate the biosynthesis of sex pheromone precursor, TAG, thus influencing PBAN-induced sex pheromone production and subsequent mating behavior.

Identification of differentially expressed genes in the pheromone glands of mated and virgin Bombyx mori by digital gene expression profiling

Background

Mating decreases female receptivity and terminates sex pheromone production in moths. Although significant progress has been made in elucidating the mating-regulated inactivation of pheromone biosynthesis-activating neuropeptide (PBAN) secretion, little is known about the mating induced gene expression profiles in pheromone glands (PGs). In this study, the associated genes involved in Bombyx mori mating were identified through digital gene expression (DGE) profiling and subsequent RNA interference (RNAi) to elucidate the molecular mechanisms underlying the mating-regulated gene expression in PGs.

Results

Eight DGE libraries were constructed from the PGs of mated and virgin females: 1 h mating (M1)/virgin (V1) PGs, 3 h mating (M3)/virgin (V3) PGs, 24 h mating (M24)/virgin (V24) PGs and 48 h mating (M48)/virgin (V48) PGs (M48 and V48). These libraries were used to investigate the gene expression profiles affected by mating. DGE profiling revealed a series of genes showing differential expression in each set of mated and virgin female samples, including immune-associated genes, sex pheromone synthesis-associated genes, juvenile hormone (JH) signal-associated genes, etc. Most interestingly, JH signal was found to be activated by mating. Application of the JH mimics, methoprene to the newly-emerged virgin females leaded to the significant reduction of sex pheromone production. RNAi-mediated knockdown of putative JH receptor gene, Methoprene tolerant 1 (Met1), in female pupa resulted in a significant decrease in sex pheromone production in mature females, suggesting the importance of JH in sex pheromone synthesis.

Conclusion

A series of differentially expressed genes in PGs in response to mating was identified. This study improves our understanding of the role of JH signaling on the mating-elicited termination of sex pheromone production.

Molecular identification of a pancreatic lipase-like gene involved in sex pheromone biosynthesis of Bombyx mori

Cytoplasmic lipid droplet (LD) lipolysis is regulated by pheromone biosynthesis activating neuropeptide (PBAN) in Bombyx mori. To elucidate the molecular mechanism of cytoplasm LD lipolysis, the pancreatic lipase-like gene in B. mori pheromone glands (PGs), designated as B. mori pancreatic lipase-like gene (BmPLLG), was identified in this study. Spatial expression analysis revealed that BmPLLG is a ubiquitous gene present in all studied tissues, such as PGs, brain, epidermis, egg, midgut, flight muscle and fat body. Temporal expression analysis showed that the BmPLLG transcript begins to express 96 h before eclosion (-96 h), continues to increase, peaks in newly emerged females and steadily decreases after eclosion. Translational expression analysis of BmPLLG using a prepared antiserum demonstrated that BmPLLG was expressed in an age-dependent pattern at different development stages in B. mori. This finding was similar to the transcript expression pattern. Further RNA interference-mediated knockdown of BmPLLG significantly inhibited bombykol production. Overall, these results demonstrated that BmPLLG is involved in PBAN-induced sex pheromone biosynthesis and release.

Expression analysis of several antiviral related genes to BmNPV in different resistant strains of silkworm, Bombyx mori

Bombyx mori L. (Lepidoptera: Bombycidae) nucleopolyhedrovirus (BmNPV) is a highly pathogenic virus in the sericultural industry, often causing severe damage leading to large economic losses. The immune mechanisms of B. mori against this virus remain obscure. Previous studies had demonstrated Bmlipase-1, BmNox and Bmserine protease-2 showing antiviral activity in vitro, but data on the transcription levels of these proteins in different resistant strains were not reported. In order to determine the resistance level of the four different strains (P50, A35, A40, A53) and gain a better understanding of the mechanism of resistance to BmNPV in B. mori, the relative expression level of the genes coding the three antiviral proteins in larval haemolymph and midgut of different B. mori strains resistant to BmNPV was determined. The results showed that these genes expressed significantly higher in the resistant strains compared to the susceptible strain, and the differential expression levels were consistent with the LC50 values in different strains. The transcription level of the target genes almost all up-regulated in the larvae midgut and down-regulated in the haemolymph. The results indicate the correlation of these genes to BmNPV resistance in B. mori.

Transcriptome analysis of silkworm, Bombyx mori, during early response to Beauveria bassiana challenges

Host–pathogen interactions are complex processes and it is a central challenge to reveal these interactions. Fungal infection of silkworm, Bombyx mori, may induce a variety of responsive reaction. However, little is known about the molecular mechanism of silkworm immune response against the fungal infection. To obtain an overview of the interaction between silkworm and an entomopathogenic fungus Beauveria bassiana, Digital Gene Expression profiling, a tag based high-throughput transcriptome sequencing method, was employed to screen and identify differentially expressed genes (DEGs, FDR≤0.001, ∣log₂ratio∣≥1) of silkworm larvae during early response against B. bassiana infection. Total 1430 DEGs including 960 up-regulated and 470 down-regulated ones were identified, of which 627 DEGs can be classified into GO categories by Gene Ontology (GO) analysis. KEGG pathways analysis of these DEGs suggested that many biological processes, such as defense and response, signal transduction, phagocytosis, regulation of gene expression, RNA splicing, biosynthesis and metabolism, protein transport etc. were involved in the interaction between the silkworm and B. bassiana. A number of differentially expressed fungal genes were also identified by mapping the sequencing tags to B. bassiana genome. These results provided new insights to the molecular mechanism of silkworm immune response to B. bassiana infection.

Identification of genes expressed in the sex pheromone gland of the black cutworm Agrotis ipsilon with putative roles in sex pheromone biosynthesis and transport

Background

One of the challenges in insect chemical ecology is to understand how insect pheromones are synthesised, detected and degraded. Genome wide survey by comparative sequencing and gene specific expression profiling provide rich resources for this challenge. A. ipsilon is a destructive pest of many crops and further characterization of the genes involved in pheromone biosynthesis and transport could offer potential targets for disruption of their chemical communication and for crop protection.

Results

Here we report 454 next-generation sequencing of the A. ipsilon pheromone gland transcriptome, identification and expression profiling of genes putatively involved in pheromone production, transport and degradation. A total of 23473 unigenes were obtained from the transcriptome analysis, 86% of which were A. ipsilon specific. 42 transcripts encoded enzymes putatively involved in pheromone biosynthesis, of which 15 were specifically, or mainly, expressed in the pheromone glands at 5 to 120-fold higher levels than in the body. Two transcripts encoding for a fatty acid synthase and a desaturase were highly abundant in the transcriptome and expressed more than 40-fold higher in the glands than in the body. The transcripts encoding for 2 acetyl-CoA carboxylases, 1 fatty acid synthase, 2 desaturases, 3 acyl-CoA reductases, 2 alcohol oxidases, 2 aldehyde reductases and 3 acetyltransferases were expressed at a significantly higher level in the pheromone glands than in the body. 17 esterase transcripts were not gland-specific and 7 of these were expressed highly in the antennae. Seven transcripts encoding odorant binding proteins (OBPs) and 8 encoding chemosensory proteins (CSPs) were identified. Two CSP transcripts (AipsCSP2, AipsCSP8) were highly abundant in the pheromone gland transcriptome and this was confirmed by qRT-PCR. One OBP (AipsOBP6) were pheromone gland-enriched and three OBPs (AipsOBP1, AipsOBP2 and AipsOBP4) were antennal-enriched. Based on these studies we proposed possible A. ipsilon biosynthesis pathways for major and minor sex pheromone components.

Conclusions

Our study identified genes potentially involved in sex pheromone biosynthesis and transport in A. ipsilon. The identified genes are likely to play essential roles in sex pheromone production, transport and degradation and could serve as targets to interfere with pheromone release. The identification of highly expressed CSPs and OBPs in the pheromone gland suggests that they may play a role in the binding, transport and release of sex pheromones during sex pheromone production in A. ipsilon and other Lepidoptera insects.

Synthetic rates of key stored fatty acids in the biosynthesis of sex pheromone in the moth Heliothis virescens

Using a tracer-tracee approach, we fed 1-d-old virgin Heliothis virescens U-(13)C-glucose and analyzed the key labeled fatty acids, (Z)-11-hexadecenoate, hexadecanoate and octadecanoate, known to be intermediates in pheromone biosynthesis, by mass isotopomer distribution analysis. This method allowed determination of enrichment, and fractional (FSR) and absolute (ASR) synthetic rates. As expected, FSRs and ASRs for all three moieties were greater in the scotophase than photophase. However, in whole gland extracts, FSRs and ASRs of (Z)-11-hexadecenoate and hexadecanoate were much lower than those of the major pheromone component, (Z)-11-hexadecenal, determined previously. Since pheromone is made via these acids, we postulated that pheromone was produced directly and very rapidly via a small pool of acyl CoA thioesters of these acids and that the pool of acids we analyzed in our whole gland extract was largely a 'dead end' pool of excess acids (i.e., not converted directly to pheromone) stored in glycerolipids. We tested this by fractionating the whole glandular extract and analyzing the glycerolipid fraction. FSRs and ASRs for the two acids in the glycerolipid fraction were similar to those for the whole gland extract, confirming our postulate. Thus, most acetate produced in the pheromone gland is converted rapidly and directly to pheromone, while excess fatty acids are stored in glycerolipids and remain relatively inaccessible for pheromone production, at least over the two periods studied. Precursor enrichment of octadecanoate was substantially lower than that determined for the two 16-carbon acids and pheromone component. This suggests that hexadecanoate is the principal product of the multi-enzyme complex fatty acid synthase in the gland, and that octadecanoate is formed by subsequent chain elongation of hexadecanoate.