- Research article
- Open Access
Analysis of the agrotis segetum pheromone gland transcriptome in the light of Sex pheromone biosynthesis
BMC Genomics volume 16, Article number: 711 (2015)
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.
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.
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.
Moths rely heavily on sex pheromones for mate finding . Usually the females produce and emit the pheromone from a specialized structure, the sex pheromone gland located at the intersegmental membrane between the 8th and 9th abdominal segment, associated with the ovipositor at the end of the adult female abdomen . Over the last five decades, sex pheromones have been identified from more than 700 species [3, 4]. The pheromone compounds are mostly fatty acid derivatives, with carbon chain length C10-C18, with 0–4 double bonds , and an oxygenated functional group (alcohol, aldehyde, acetate ester) [3–5]. A substantial portion is made up by acetate esters . Most moths use a combination of two or more compounds in a specific ratio, constituting a more or less species-specific blend. The pheromone biosynthesis pathways have been studied extensively and are well documented in many moth species [6, 7]. Characterization of the enzymes involved in the process of pheromone biosynthesis not only helps to understand the evolution of sexual communication and speciation, but could ultimately also aid in pest control by allowing the design of drugs that block the biosynthetic machinery or by allowing the synthetic biologist to produce species-specific pheromones for mass trapping or mating disruption in biological systems like cell factories or genetically modified plants [8–10].
During the last two decades, desaturases introducing double bonds into the acyl chain in ∆6 , ∆9 [12, 13], ∆10 , ∆11 [12, 13, 15–17] and ∆14  position have been cloned from many moth species, and their functions have been characterized in various heterologous expression systems. Also an omega-desaturase that introduces a double bond in the methyl terminal carbon has been cloned and characterized . In some cases, two double bonds can be produced by one desaturase  or by the consecutive activities of two desaturases . These studies highlight desaturases as key players for the diversity of pheromone structures in moth species and their role in reproductive isolation and speciation. A variety of desaturases in combination with limited chain-shortening or chain-elongation account for the diversity of double bond isomerism observed among moth sex pheromone compounds.
After the double bonds are in place and the acyl chain length is adjusted, the carbonyl carbon is modified to form a functional group. Firstly it requires a step that converts the fatty-acyl precursors into fatty alcohols. Great progress has been made since the first fatty-acyl reductase gene was identified in Bombyx mori . Several Ostrinia spp. FARs have been characterized [21, 22]. The Ostrinia FARs are very essential to determine the final pheromone compositions , minimal changes in sequence can cause the pheromone component ratio to shift dramatically . On the other hand, there are also FARs that are very versatile in terms of substrate specificity [24, 25], involved in the biosynthesis of multicomponent pheromones.
Fatty alcohols serve as the actual pheromone components in a number of moth species , but mostly, fatty alcohols will be either oxidized to aldehydes [26–28] or esterified to form acetate esters [29–31]. Acetyltransferases are important enzymes since the acetate esters are very commonly occurring pheromone components among moth species. They have not been cloned from insects yet, but have been investigated in some cases by in vivo labeling studies [32–34]. Unlike FARs, these functional group modification enzymes have not been studied extensively.
The chain-shortening pathway has not been characterized at the enzymatic level in insects, but it was noted in a couple of cases that mutations in the β-oxidation pathway did affect the final pheromone compositions. The major pheromone component of cabbage looper moth Trichoplusia ni is Z7-12:OAc, whereas a mutant strain produced a greatly increased amount of Z9-14:OAc [35, 36].
The turnip moth, A. segetum uses a series of short chain acetate esters as sex pheromone components, including the homologues (Z)-5-decenyl, (Z)-7-dodecenyl, and (Z)-9 tetradecenyl acetate (Z5-10:OAc, Z7-12:OAc, and Z9-14:OAc) [37–39]. The biosynthesis of this type of pheromone involves desaturation of palmitic acid (16C), a product of the ubiquitous fatty acid synthase machinery. The unsaturated fatty acid undergoes chain shortening, reduction, and acetylation . Populations of A. segetum from different geographic areas differ in the ratio of their pheromone components [40, 41]. For instance the Swedish population has a ratio of Z9-14:OAc/Z7-12:OAc/Z5-10:OAc = 29/59/12 whereas the Zimbabwean population has a ratio of 2/20/78 . This shift in ratios could be due to differences in the chain-shortening, in the FAR activity or less likely in the activity of the acetyltransferase.
An EST-library was previously constructed from the A. segetum pheromone gland, revealing candidate genes involved in pheromone production  and a ∆11 desaturase and a FAR involved in pheromone production were already characterized . The EST library, however, contains only 2, 285 objects. We now report a more extensive database of candidate genes potentially involved in the A. segetum sex pheromone biosynthetic pathway, generated by next generation sequencing technology (NGS). We constructed transcriptome libraries from two tissues of A. segetum: pheromone gland (As_PG) and abdomen (As_AB). Furthermore, to test the function of several of the candidate genes we assayed them in a yeast heterologous expression system.
Results and discussion
Illumina sequencing, unigene assembly, and analysis of transcripts
In total 53 million raw clean reads were obtained from each (As_PG and As_AB) library, with a total of 4.8G clean nucleotides in each (Table 1). The clean reads from the two libraries were assembled into 62,165 consensus contigs (Table 2), including 22,633 distinct clusters (referred to as CL) and 39,532 distinct singletons (referred to as Unigene). These consensus contigs have a mean length of 733 nt, and N50  =1,150, total length 45 Mb. Size distributions of the unigenes can be seen in Fig. 1 and their differential expression in the two tissues is displayed in Fig. 2. This transcriptome shotgun assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GBCW00000000. The version described in this paper is the first version, GBCW01000000.
After alignment by blastx to protein databases NCBI Nr, Swiss-Prot, KEGG and COG (e-value < 0.00001), and alignment by blastn to nucleotide database NCBI Nt, annotations were retrieved from the highest sequence similarity with the given unigenes along with their protein functional annotations. About half of the unigenes have a hit in one or more of the databases (Table 3). 14,108 unigenes hit the first record in the Nr database with Bombyx mori (Table 4). A. ipsilon got 293 hits, but this does not mean the A. segetum is closer to B. mori than A. ipsilon, it is just because the B. mori got more records deposited in GenBank than A. ipsilon. In terms of clusters of orthologous groups (COG), 19,911 unigenes were classed in one or more of the 26 COG functional categories (Fig. 3).
By running the Blast2GO program  using Nr annotation, 11,246 (18 %) unigenes were assigned to one or more GO categories. The GO terms cellular process, binding, catalytic activity, and metabolic process were the most abundantly represented categories (over 5,000 unigenes, details see Fig. 4). These numbers and percentages are similar to the results that Gu and coauthors  presented from A. ipsilon.
The biosynthetic pathway (Fig. 5) leading to the pheromone components of A. segetum is similar to what has been reported for many other moth species [37, 45]. The key players among enzymes involved or postulated to be involved are desaturases, β-oxidation enzymes, fatty-acyl reductases, and acetyltransferases. Pheromone biosynthesis is reported to be under control by a pheromone biosynthesis activating neuropeptide (PBAN) . In the following we present candidate genes related to each step in the biosynthetic pathway, their expression levels and their functional assay in yeast heterologous expression systems.
Pheromone biosynthesis activating neuropeptide (PBAN) receptor
PBAN is released from the subesophageal ganglion (located near the brain) and is transported through hemolymph (or via the ventral nerve cord) to the pheromone gland. Upon binding to the PBAN receptor present on the pheromone gland cell membrane , it will induce the opening of calcium channels causing influx of extracellular calcium . Then the calcium binds to calmodulin, that stimulates phosphatase (and/or kinase), which subsequently activates the FAR in the case of Bombyx mori  or in other cases the acetyl-CoA carboxylase  or the acetyltransferase  maybe regulated. In As_PG, we found one gene, Ase_17579 [GenBank: KJ622075], which is 97 % identical to Helicoverpa zea PBAN receptor. Its expression level is similar in As_PG and As_AB, 15.6 FPKM and 16.9 FPKM, respectively.
The rate-limiting step in fatty acid biosynthesis  is the ATP-dependent carboxylation of acetyl-CoA to malonyl-CoA catalyzed by acetyl-CoA carboxylase (ACCase) , the first step in saturated long chain fatty acid biosynthesis. ACCase is a large protein with multiple catalytic activities, working coordinately and providing malonyl-CoA substrate for the biosynthesis of fatty acids . In A. segetum, we found one contig Ase_7442 (KJ622074), encoding the full-length of this protein. It shares 69 % amino acid identity with the ACCase of Tribolium castaneum and 88 % aa identity with B. mori ACCase. The expression level of Ase_7442 is not significantly different in As_AB (37 FPKM) and in As_PG (26 FPKM).
Fatty acid synthase
Fatty acid synthase (FAS) is a multifunctional protein  that produces saturated fatty acids using malonyl-CoA and acetyl-CoA as substrate and that requires NADPH as reducing agent, in a cyclic process in which an acetyl primer undergoes a series of decarboxylative condensations with several malonyl moieties . The resulting products are palmitic and stearic acid in insects, as proven by labeling studies [49, 53, 54]. We found six unigenes (KJ622068-KJ622073) that are homologous to the Aip_FAS_JX989151 . In total, these unigenes are seven times more expressed (significant) in As_PG than in As_AB (FPKM 67.1 / FPKM 9).
Fatty acid desaturases catalyze the introduction of double bonds into acyl chains with strict regioselectivity and stereoselectivity, and can be divided into four categories : 1) first desaturases, inserting a double bond into the saturated acyl chain; 2) front-end desaturases, introducing a double bond between an existing double bond and the carboxylic end; 3) omega desaturases, inserting a double bond between an existing double bond and the methyl end; 4) sphingolipid desaturases, introducing double bonds into sphingolipids which are important components of eukaryotic plasma membranes. We found 10 desaturase candidates from the two tissues, and they belong to the first desaturase (7/10), the front-end desaturase (2/10), and the sphingolipid desaturase (1/10) subfamilies. Among the first desaturase subfamily, several groups have been recognized based on phylogenetic and functional analysis and a four-letter “signature motif” has been suggested for each of them as shown in Fig. 6. These signature motifs strongly associate with the location of double bond that the desaturase is inserting in the fatty-acyl chain, with less emphasis on the chain length selectivity: ∆9_KPSE (C16 > C18) desaturases and ∆9_NPVE (C18 > C16) desaturases that are mostly involved in fatty acid metabolism , and ∆11,∆10 and bifunctional desaturases with the “xxxQ” motif (with a few exceptions having “xxxE” motif) exclusively involved in pheromone biosynthesis [13, 57] (Fig. 6). The Ase_1623 and Ase_4567 have the signature motif of KPSE and NPVE, respectively, that strongly suggest their roles in ordinary metabolic pathways as ∆9 desaturases. The Ase_21308 and Ase_5534 both have “xxxQ” motif, suggesting that they represent ∆11 desaturases. The Ase_5534 is the obvious candidate for pheromone biosynthesis since its expression level in As_PG is 1–2 magnitudes higher than any other candidate, and it has a very low expression level in As_AB. The functional assay confirms that it possess the ability of creating a ∆11-double bond mainly on palmitate (Fig. 7), corroborating that Ase_5534 [GenBank:KJ622051] is the desaturase responsible for pheromone production in A. segetum. This result is consistent with the recent study in which this gene was expressed in different yeast strains  and this desaturase was as a matter of fact found already in the EST analysis reported by Strandh et al. . This desaturase introduces a double bond with Z configuration on 16:Acyl. Besides Z11-16:Acyl, some other minor products were detected: ∆11-12:Acyl, Z11-14:Acyl, E11-14:Acyl, and Z11-15:Acyl (Fig. 7). The Ase_21308 displayed no activity in our yeast expression system. Gu et al.  found a desaturase that is very close to our Ase_5534 (as shown in Fig. 6) and which should be the one involved in pheromone biosynthesis in A. ipsilon. Compared to [42, 45, 58], our dataset includes a larger set of desaturases isolated from the moth pheromone gland and surrounding tissues.
After the palmitate has undergone ∆11 desaturation in the A. segetum pheromone gland, it is subject to limited chain shortening by β-oxidation, resulting in three homologous fatty-acyl pheromone precursors with 14C, 12C, and 10C chain length (Fig. 5). Chain-shortening by β-oxidation is the action of a series of enzymes, working sequentially and forming a reaction spiral.
In the first reaction, acyl-CoA is converted to E2-enoyl-CoA, by acyl-CoA oxidases (in peroxisomes) and acyl-CoA dehydrogenases (in mitochondria). Four acyl-CoA dehydrogenases with different chain length specificities cooperate to assure that the complete degradation of all fatty acids with different chain length. The names of the four dehydrogenases, short-chain, medium-chain, long-chain, and very long-chain acyl-CoA dehydrogenases, reflect their chain-length specificities. Short-chain acyl-CoA dehydrogenase only acts on short-chain substrates like butyryl-CoA and hexanoyl-CoA. Medium-chain acyl-CoA dehydrogenase is most active with substrates from hexanoyl-CoA to dodecanoyl-CoA, whereas long-chain acyl-CoA dehydrogenase preferentially acts on octanoyl-CoA and longer-chain substrates . Very-long-chain acyl-CoA dehydrogenase extends the activity spectrum to longer-chain substrates, including those having acyl chains of 22 and 24 carbon atoms . We found a full spectrum of dehydrogenases and oxidases in A. segetum (Table 5). It is noteworthy that unigene6715 was expressed 72 times higher (386FPKM) in As_PG than in As_AB (5FPKM). In addition, we found two unigenes of isovaleryl-CoA dehydrogenase, which is specific for metabolism of branched-chain fatty acids .
The second step of β-oxidation E2-enoyl-CoA is reversibly hydrated by enoyl-CoA hydratase to L-3-hydroxyacyl-CoA. Two categories of enoyl-CoA hydratases have been identified in mitochondria . One is specialized for crotonyl-CoA (4C). The second one is long-chain enoyl-CoA hydratase, effectively hydrates medium-chain and long-chain substrates. Long-chain enoyl-CoA hydratase is a component enzyme of the trifunctional β-oxidation complex, which additionally exhibits long-chain activities of L-3-hydroxyacyl-CoA dehydrogenase and 3-ketoacyl-CoA thiolase , and resides in the inner mitochondrial membrane. We found six enoyl-CoA hydratases from A. segetum (Table 5), with similar expression level in both As_PG and As_AB.
The third reaction is the reversible dehydrogenation of L-3-hydroxyacyl-CoA to 3-ketoacyl-CoA catalyzed by L-3-hydroxyacyl-CoA dehydrogenase. Four categories of L-3-hydroxyacyl-CoA dehydrogenases have been identified in mitochondria. Long chain L-3-hydroxyacyl-CoA dehydrogenase is a component enzyme of the trifunctional β-oxidation , which is most active with long-chain substrates. Medium-chain and short-chain L-3-hydroxyacyl-CoA dehydrogenase, are both soluble matrix enzymes, processing medium- and short-chain substrates. These three enzymes complement each other and thus assure high rates of dehydrogenation over the whole spectrum of β-oxidation intermediates . In A. segetum this group of enzymes is represented by five unigenes (Table 5), with similar expression level in both tissues.
The final step in which 3-ketoacyl-CoA is cleaved by thiolase between its α and β carbon atoms, makes the substrate two carbons shorter. Three classes of thiolase exist in motochondria: acetoacetyl-CoA thiolase or acetyl-CoA acetyltransferase (specific for acetoacetyl-CoA), 3-ketoacyl-CoA thiolase or acetyl-CoA acyltransferase (act on C4-C16), and long-chain 3-ketoacyl-CoA thiolase that is a component enzyme of the membrane-bond trifunctional β-oxidation complex, whereas the first two thiolases are soluble matrix enzymes . We found six unigenes of 3-ketoacyl-CoA thiolase from A. segetum (Table 5), and all of them expressed at similar level in As_PG and As_AB.
The degradation of unsaturated fatty acids requires auxiliary enzymes like ∆3,∆2-enoyl-CoA isomerase and 2,4-dienoyl-CoA reductase to modify the structure of double bonds during the β-oxidation process to ensure a continuous flow of intermediates through the β-oxidation spiral . We found two unigenes of the ∆3,∆2-enoyl-CoA isomerase, one mitochondrial type and one peroxisomal. In addition, we found a ∆3,5∆2,4-dienoyl-CoA isomerase that is specialized for processing odd-numbered double bonds .
Insects in general have the ability to shorten long-chain fatty acids to a specific shorter chain length . Jurenka [5, 66] suggests that this kind of limited chain shortening takes place in the peroxisomes. Unlike the β-oxidation in mitochondria in which substrates are thoroughly degraded into two-carbon units, the β-oxidation in peroxisomes ceases at the formation of medium-chain fatty-acyl-CoAs, because acyl-CoA oxidase is inactive toward substrates having acyl moieties of eight or fewer carbon atoms . Moth pheromone components are commonly 12C and 14C, suggesting that there is still something specific about the chain-shortening being involved in moth pheromone biosynthesis.
In the present study we found representatives of all the key players and auxiliary enzymes of β-oxidation, and some of them are very high in expression level and differentially expressed among the As_PG and As_AB (Table 5), forming promising candidates to tackle the role of β-oxidation on pheromone biosynthesis either by heterologous expression or by RNAi.
Chain-shortened fatty-acyl precursors are reduced to the corresponding alcohols by Fatty-Acyl Reductases (FARs) . FARs have been cloned and characterized from several moth species [21–23] since the first one was found in B. mori . We found ten full-length unigenes in both the As_PG and As_AB (Fig. 8) and two of them cluster within the pgFAR clade . We expressed the two ORFs that belong to the pheromone-producing clade in our yeast expression system. The results showed that it is the one named Ase_1929 that is responsible for reducing all the three fatty-acyl precursors (Z5-10:CoA, Z7-12:CoA, Z9-14:CoA) into their corresponding fatty alcohols (Fig. 9), and it is the most abundant transcript among all these ten unigenes (Fig. 8). This result is consistent with previous findings that a single FAR takes a wide range of the fatty-acyl substrates and convert them into fatty alcohols [24, 25, 67]. Our results complement the findings of Hagström et al. , who presented the Ase_1929 with longer chain acyl substrates (16C) and found it to be very active on them. But what appears in Fig. 9 as if the Ase_1929 when expressed in yeast is more active towards longer chain length substrates than towards shorter ones, does not necessarily represent its actual activity in the pheromone gland. This difference could partly be due to the volatility of the substrates and products, since shorter chain methyl esters and alcohols are more volatile and thus may escape from the yeast expression system during incubation. The other unigene Ase_20982 which also clustered with the pgFAR clade, did not show any function in our heterologous expression system. The FAR [GenBank: JX989146, protein ID: AGR49323] found by Gu et al.  clustered very close to our Ase_1929, suggesting that it is likely involved in pheromone biosynthesis in A. ipsilon.
The genes involved in acetylation of fatty pheromone alcohols have not been cloned from any insect species. Acetyltransferases probably belong to a huge family of acyl CoA-utilizing enzymes whose products include a variety of chemicals, such as neurotransmitters , plant volatile esters [69, 70], constitutive defense compounds, waxes , phytoalexins, lignin, phenolics, alkaloids, and anthocyanins [72–74], which makes it very difficult to make functional predictions from primary sequence information alone. Attempts have been made but ended up with getting other member of the family . This step is the last step in the A. segetum pheromone biosynthetic pathway. Most likely it does not significantly influence the ratio between the homologous acetate pheromone components [32, 33], although in several moth species the Z isomers are produced faster than the E isomers and yield more products, when the alcohol substrates were supplied to the pheromone gland homogenate . Previously, a number of acetyltransferases that were cloned from plant species were studied. For example, the kiwi alcohol acetyltransferase AT9 produced butyl acetate and butyl propionate with highest catalytic activities when tested with short chain alcohols . Another example is the apple alcohol acyltransferases (MpAAT1) that use coenzyme A (CoA) donors together with alcohol acceptors as substrates, which were cloned and characterized . The MpAAT1 recombinant enzyme can utilize a range of alcohol substrates from short to medium straight chain (C3-C10), branched chain, aromatic and terpene alcohols. The enzyme can also utilize a range of short to medium chain CoAs , but alcohols longer than C10 have not been tested. Similar work was done with enzymes forming volatiles esters from banana, melon, and strawberry. The recombinant enzymes were capable of producing esters from a wide range of alcohols and acyl-CoA [70, 77, 78]. Overall, these enzymes are members of the BAHD [73, 74], generally recognized by their active site motif (HXXXD) and a conserved region (DFGWG) with likely structural significance . A tBLASTn search with current published BAHDs as query against our A. segetum transcriptome got no hit, suggesting this moth may not express this gene family in its pheromone gland or they have undergone substantial evolutionary changes. But according to previously published putative acetyltransferase [45, 58] and the annotation results we found 34 candidates in our transcriptome (Table 6) and we heterologously expressed them in our yeast system. The ORF of each of these putative acetyltransferases was thus cloned into a yeast expression vector, pYES-DEST52, under the control of galactose inducible promoter. The ∆ATF1 yeast strain (knock out strain that lacking most of the acetylation activity) was transformed with an individual construct and incubated in liquid culture with the culture medium supplemented with a mixture of Z9-14:OH, Z7-12:OH and Z5-10:OH. After incubation for 2 days, the total lipids were extracted and analyzed by GC-MS for acetate esters. The results did not reveal any of insect candidate genes being capable of esterifying fatty alcohols into acetate esters (Fig. 10), whereas the strain overexpressing ATF1 was highly active (positive control).
We explored the data obtained from massive sequencing of pheromone producing tissue of the turnip moth and compared the expression levels of candidate unigenes with their expression levels in the abdominal epidermal tissue. This allowed identification of key parts involved in the pheromone biosynthesis pathway such as: β-oxidation enzymes, a desaturase, a fatty-acyl reductase, putative acetyltransferases, and other components involved in the fatty-acid metabolism, like acetyl-CoA carboxylase and fatty-acid synthetase. By phylogenetic analyses of desaturases and FARs, we found the most promising candidates for each gene and confirmed their function in pheromone biosynthesis. The ∆11-desaturase is a specialist interacting preferentially with 16:CoA and producing Z11-16:CoA. The FAR we found is a generalist that can reduce a broad range of saturated and unsaturated acyl substrates from 10C to 16C. We had specifically hoped to be able to clone and characterize an acetyltransferase involved in pheromone biosynthesis that has been postulated in many studies. We tested 34 genes that were annotated to be acetyltransferases in our yeast expression system (which can successfully express plant-derived acetyltransferase ) but it turned out that none of them was functional in comparison to the ATF1 positive control. The nature of the acetyltransferase involved in moth pheromone biosynthesis remains illusive. In addition, we generated a list of β-oxidation enzymes, ACC and FAS that can be further tested either by heterologous expression or by RNAi.
Insects and tissue collection
Agrotis segetum were obtained from a laboratory culture, continuously maintained for more than 20 years in Lund, but repeatedly rejuvenated by addition of field-collected insects. The larvae were reared on a semisynthetic bean-based diet  and kept at 25 °C under a 16 h:8 h light: dark cycle. Pupae were separated by sex and placed in different jars. The third day after they emerged, 30 pheromone glands [58, 80] and abdominal epidermal tissue  from 3 individuals of female moths were dissected 3–4 h into scotophase  and stored in −80 °C freezer until RNA extraction.
Total RNA from pheromone gland and abdominal tissue were extracted using the TRIzol reagent (Life Technologies, Lidingö, Sweden) according to the manufacturer’s instructions except for three additional ethanol washes before dissolving RNA in water. RNA concentration and purity were checked on NanoDrop2000 (Thermo Scientific, Saveen Werner, Malmö, Sweden).
Illumina sequencing and bioinformatic analysis
Twenty μg of total RNA from each sample were sent to BGI (Hong Kong Co., Ltd) for library construction, Illumina sequencing and subsequent bioinformatic analysis.
Reads were assembled using Trinity  into contigs. Then the reads are mapped back to contigs, get sequences without Ns and cannot be extended on either end. Such sequences are defined as unigenes. Then TGICL  is used to assemble all the unigenes from As_PG and As_AB to form a single set of non-redundant unigenes. Then gene family clustering was performed and unigenes were divided to two classes. One is clusters, which the prefix is CL and the cluster id is behind. In one cluster, there are several unigenes which similarity between them is more than 70 %. The other is singletons, which the prefix is unigene. Unigene sequences are aligned with blastdb using blastx (E-value < 0.00001). Sequence orientations are determined according to the best hit in the database. Bioinformatic data were viewed and further processed using Geneious version (6.1.6), created by Biomatters and available from http://www.geneious.com.
The calculation of unigene expression was performed using the FPKM method (Fragments Per kb per Million fragments) [84, 85], the formula is FPKM = 106C/NL/103. Set FPKM to be the expression of unigene A, and C to be number of fragments that uniquely aligned to unigene A, N to be total number of fragments that uniquely aligned to all unigenes, and L to be the base number in the CDS of unigene A. The FPKM method eliminates the influence of different gene lengths and sequencing levels on the calculation of gene expression. Therefore the calculated gene expression can be directly used for comparing the differences in gene expression between samples .
Functional annotations of unigenes were conducted, based on protein sequence similarity, towards the KEGG Pathway, COG  and Gene Ontoloty (GO) databases. Briefly, we search all unigene sequences against protein databases (NR, SwissProt, KEGG, COG) using blastx (E-value < 0.00001). Based on NR annotation, we use Blast2GO program  to get GO annotation of all unigenes. After getting GO annotation for every unigene, we use WEGO software  to do GO functional classification for all unigenes.
Sequences used for phylogenetic reconstructions were retrieved from the GenBank (http://www.ncbi.nlm.nih.gov) database. Neighbor-Joining trees were constructed using Mega version 4.0 . Briefly, multiple sequence alignments were run using the MAFFT (online version) and the output FASTA format data were plugged into Mega software. Genetic distance model was JTT, 1500 replicates, Neighbor-Joining as tree building method.
Fatty acid precursors
Z9-18:Me, 16:Me were purchased from Larodan Fine Chemicals AB (Malmö, Sweden). Z5-10:OH, Z7-12:OH, Z9-14:OH were purchased from Pherobank (Wageningen, The Netherlands, 98 % purity). Z9-14:Me, Z7-12:Me, Z5-10:Me were prepared from their corresponding alcohols by previously described methods . All FAMEs and OH were dissolved in 96 % ethanol in a 0.1 M stock solution. All alcohols used as reference compounds were from our laboratory collection of pheromone compounds.
Construction of expression vector for functional assay
For the construction of yeast expression vectors containing the candidate genes, specific primers with attB1 and attB2 sites incorporated were designed for amplifying the ORF of genes of interest. The PCR products were subjected to agarose gel electrophoresis and purified using the Wizard® SV Gel and PCR Clean up system (Promega Biotech AB, Nacka, Sweden). The ORFs were subcloned into the pDONR221 vector in presence of BP clonase (Life Technologies), after confirmation by sequencing, the correct entry clones were selected and do LR reaction with pYES2-DEST52 (for FARs and acetyltransferases) or pYEX-CHT-DEST vector (for desaturases), and resulting expression clones were analyzed by sequencing.
Functional assay in yeast
The resultant recombinant expression vectors harboring the candidate genes were introduced into the INVSc (MATa HIS3 LEU2 trp1-289 ura3-52) (for FARs), the ∆ATF1 knockout strain (for acetyltransferase), or the double deficient ole1 elo1 strain (MATa elo1::HIS3 ole1::LEU2 ade2 his3 leu2 ura3) (for desaturase) of the yeast Saccharomyces cerevisiae  using the S.c. easy yeast transformation kit (Life Technologies). For selection of uracil (and leucine) prototrophs, the transformed yeast was allowed to grow on SC plate containing 0.7 % YNB (w/o aa, with Ammonium sulfate) and a complete drop-out medium lacking uracil (and leucine) (Formedium™ LTD, Norwich, England), 2 % glucose, 1 % tergitol (type Nonidet NP-40, Sigma-Aldrich Sweden AB, Stockholm, Sweden), 0.01 % adenine (Sigma) and containing 0.5 mM oleic acid (Sigma) as extra fatty acid source. After 2 days (7 days for ole1 elo1 strain) at 30 °C, individual colonies were picked up to inoculate 10 mL selective medium at 30 °C and grown at 300 rpm for 48 h. Yeast cultures were diluted to an OD600 of 0.4 in 10 mL fresh selective medium containing 2 % galactose (2 mM CuSO4) with supplementation of a biosynthetic precursor. Each FAME or fatty alcohol precursor was prepared at a concentration of 100 mM in 96 % ethanol and added to reach a final concentration of 0.5 mM in the culture medium . In the acetyltransferase assay the yeasts were supplemented with a mixture of the three alcohols Z9-14:OH, Z7-12:OH and Z5-10:OH. Yeasts were cultured in 30 °C in a shaking incubator at 30 °C.
Fatty acid/alcohol/acetate analysis
After 48 h of incubation yeast cells were harvested by centrifugation at 3,000 rpm. For the analysis of desaturase products, total lipids were extracted using 3.75 mL of methanol/chloroform (2:1, v/v), in a glass tube. One mL of HAc (0.15 M) and 1.25 mL of water were added to the tube to wash the chloroform phase. Tubes were vortexed vigorously and centrifuged at 2000 rpm for 2 min. The bottom chloroform phase, about 1 mL, containing the total lipids, were transferred to a new glass tube. Fatty acid methylesters (FAMEs) were made from this total lipid extract. The solution of total lipids was evaporated to dryness under gentle nitrogen flow. One mL of sulfuric acid (2 % in methanol) was added to the tube, which was then vortexed vigorously and incubated at 90 °C for an hour. After incubation, 1 mL of water was added, mixed well, and then 1 mL of hexane was used to extract the FAMEs .
Fatty alcohols and acetates were extracted from cells using 800 μL of hexane plus sonication. After brief centrifugation, the supernatant was transferred to a new tube and subjected to GC-MS analysis.
Double bond positions were confirmed by dimethyl disulfide (DMDS) derivatization , followed by GC-MS analysis. FAMEs (50 μL) were transferred to a new tube and 50 μL DMDS was added and incubated at 40 °C overnight, in the presence of 5 μL of iodine (5 % in diethyl ether) as catalyst. Hexane (200 μL) was added to the sample and the reaction was neutralized by addition of 50-100 μL Na2S2O3 (5 % in water). The organic phase was recovered and concentrated under a gentle nitrogen stream to 40-50 μL.
Gas chromatography - mass spectrometry (GC-MS)
The methylesters, fatty alcohols and acetates were subjected to GC-MS analyses on a Hewlett Packard 6890 GC coupled to a mass spectrometer HP 5973. The GC was equipped with an INNOWax column (30 m × 0.25 mm i.d. × 0.25 μm film thickness, Agilent Technologies), and helium was used as carrier gas (average velocity: 33 cm/s). The MS was operated in electron impact mode (70 eV), scaning between m/z 30 and m/z 400, and the injector was configured in splitless mode at 220 °C. The oven temperature was set to 80 °C for 1 min, then increased at a rate of 10 °C/min up to 210 °C, followed by a hold at 210 °C for 15 min, and then increased at a rate of 10 °C/min up to 230 °C followed by a hold at 230 °C for 20 min.
DMDS derivatives were analyzed on an Agilent 6890 GC system equipped with HP-5MS capillary column (30 m × 0.25 mm i.d. × 0.25 μm film thickness, Agilent Technologies) coupled with an HP 5973 mass spectrometer. The oven temperature was set at 80 °C for 1 min, raised to 140 °C at a rate of 20 °C/min, then to 250 °C at a rate of 4 °C/min and held for 20 min .
Data were analyzed using the ChemStation software (Agilent, Technologies, USA).
β-oxidation enzymes: KJ622076–KJ622113
fatty acid synthase
gas chromatography mass spectrometry
polymerase chain reaction
(Z)-9-tetradecenoic methyl ester
(Z)-7-dodecenoic methyl ester
(Z)-5-decenoic methyl ester
Wyatt TD. Pheromones and Animal Behaviour. Communication by Smell and Taste. Cambridge: Cambridge University Press; 2003.
Percy-Cunningham JE, MacDonald JA. Biology and ultrastructure of sex pheromone-producing glands. In: Prestwich GD, Blomquist GJ, editors. Pheromone biochemistry. Orlando Fl: Academic; 1987. p. 27–69.
El-Sayed AM. The Pherobase: Database of Pheromones and Semiochemicals. 2011, [http://www.pherobase.net]
Jurenka R. Insect pheromone biosynthesis. Topics in Curr Chem. 2004;239:97–132.
Blomquist GJ, Jurenka R, Schal C, Tittiger C. Biochemistry and molecular biology of pheromone production. In: Gilbet LI, Latrou K, Gill SS, editors. Comprehensive molecular insect science. London: Elsevier; 2005. p. 705–51.
Tillman JA, Seybold SJ, Jurenka RA, Blomquist GJ. Insect pheromones - an overview of biosynthesis and endocrine regulation. Insect Biochem Mol Biol. 1999;29:481–541.
Cook SM, Khan ZR, Pickett JA. The use of push-pull strategies in integrated pest management. Annu Rev Entomol. 2007;52:375–400.
Hagström ÅK, Wang HL, Liénard MA, Lassance JM, Johansson T, Löfstedt C. A moth pheromone brewery: production of (Z)-11-hexadecenol by heterologous co-expression of two biosynthetic genes from a noctuid moth in a yeast cell factory. Microb Cell Fact. 2013;12:125.
Ding BJ, Hofvander P, Wang HL, Durrett TP, Stymne S, Löfstedt C. A plant factory for moth pheromone production. Nat Commun. 2014;5:3353.
Wang HL, Liénard MA, Zhao CH, Wang CZ, Löfstedt C. Neofunctionalization in an ancestral insect desaturase lineage led to rare ∆6 pheromone signals in the Chinese tussah silkworm. Insect Biochem Mol Biol. 2010;40:742–51.
Hao GX, O’Conner M, Liu WT, Roelofs WL. Chraracterization of Z/E11- and Z9-desaturases from the obliquebanded leafroller moth, Choristoneura rosaceana. J Insect Sci. 2002;2:26.
Liénard MA, Strandh M, Hedenström E, Johansson T, Löfstedt C. Key biosynthetic gene subfamily recruited for pheromone production prior to the extensive radiation of Lepidoptera. BMC Evol Biol. 2008;8:270.
Hao G, Liu W, O'Connor M, Roelofs WL. Acyl-CoA Z9 and Z10-desaturase genes from a New Zealand leafroller moth species, Planotortrix octo. Insect Biochem Mol Biol. 2002;32:961–6.
Liu WT, Jiao HM, O'Connor M, Roelofs WL. Moth desaturase characterized that produces both Z and E isomers of ∆11-tetradecenoic acids. Insect Biochem Mol Biol. 2002;32:1489–95.
Liu WT, Rooney AP, Xue BY, Roelofs WL. Desaturases from the spotted fireworm moth (Choristoneura parallela) shed light on the evolutionary origins of novel moth sex pheromone desaturases. Gene. 2004;342:303–11.
Moto K, Suzuki MG, Hull JJ, Kurata R, Takahashi S, Yamamoto M, et al. Involvement of a bifunctional fatty-acyl desaturase in the biosynthesis of the silkmoth, Bombyx mori, sex pheromone. Proc Natl Acad Sci U S A. 2004;101:8631–6.
Roelofs WL, Liu W, Hao G, Jiao H, Rooney AP, Linn Jr CE. Evolution of moth sex pheromones via ancestral genes. Proc Natl Acad Sci U S A. 2002;99(21):13621–6.
Ding BJ, Liénard MA, Wang HL, Zhao CH, Löfstedt C. Terminal fatty-acyl-CoA desaturase involved in sex pheromone biosynthesis in the winter moth (Operophtera brumata). Insect Bioch Mol Biol. 2011;41:715–22.
Moto K, Yoshiga T, Yamamoto M, Takahashi S, Okano K, Ando T, et al. Pheromone gland-specific fatty-acyl reductase of the silkmoth, Bombyx mori. Proc Natl Acad Sci U S A. 2003;100(16):9156–61.
Antony B, Fujii T, Moto K, Matsumoto S, Fukuzawa M, Nakano R, et al. Pheromone-gland-specific fatty-acyl reductase in the adzuki bean borer, Ostrinia scapulalis (Lepidoptera: Crambidae). Insect Biochem Mol Biol. 2009;39(2):90–5.
Lassance JM, Groot AT, Liénard MA, Antony B, Borgwardt C, Andersson F, et al. Allelic variation in a fatty-acyl reductase gene causes divergence in moth sex pheromones. Nature. 2010;466:486–9.
Lassance JM, Liénard MA, Antony B, Qian S, Fujii T, Tabata J, et al. Functional consequences of sequence variation in the pheromone biosynthetic gene pgFAR for Ostrinia moths. Proc Natl Acad Sci U S A. 2013;110:3967–72.
Liénard MA, Hagström ÅK, Lassance JM, Löfstedt C. Evolution of multicomponent pheromone signals in small ermine moths involves a single fatty-acyl reductase gene. Proc Natl Acad Sci U S A. 2010;107:10955–60.
Hagström ÅK, Liénard MA, Groot AT, Hedenström E, Löfstedt C. Semi-selective fatty acyl reductases from four heliothine moths influence the specific pheromone composition. PLoS ONE. 2012;7(5):e37230.
Teal PEA, Tumlinson JH. Properties of cuticular oxidases used for sex pheromone biosynthesis by Heliothis zea. J Chem Ecol. 1998;14:2131–45.
Fang N, Teal PEA, Tumlinson JH. Correlation between glycerolipids and pheromone aldehydes in the sex pheromone gland of female tobacco hornworm moths, Manduca sexta (L.). Arch Insect Biochem Physiol. 1995;4:321–36.
Hoskovec MA, Luxová A, Svatoš A, Boland W. Biosynthesis of sex pheromones in moths: stereochemistry of fatty alcohol oxidation in Manduca sexta. Tetrahedron. 2002;58:9193–201.
Bestmann HJ, Herrig M, Attygalle AB. Terminal acetylation in pheromone biosynthesis by Mamestra brassicae L. (Lepidoptera:Noctuidae). Experientia. 1987;43:1033–4.
Teal PEA, Tumlinson JH. The role of alcohols in pheromone biosynthesis by two noctuid moths that use acetate pheromone components. Arch Insect Biochem Physiol. 1987;4:261–9.
Zhu JW, Zhao CH, Bengtsson M, Löfstedt C. Reductase specificity and the ratio regulation of E/Z isomers in the pheromone biosynthesis of the European corn borer, Ostrinia nubilalis (Lepidoptera: Pyralidae). Insect Biochem Mol Biol. 1996;26:171–6.
Morse D, Meighen E. Biosynthesis of the acetate ester precursor of the spruce budworm sex pheromone by an acetyl CoA: Fatty alcohol acetyltransferase. Insect Biochem. 1987;17:53–9.
Zhao CH, Lu F, Bengtsson M, Löfstedt C. Substrate specificity of acetyltransferase and reductase enzyme systems used in pheromone biosynthesis by Asian corn borer, Ostrinia furnacalis. J Chem Ecology. 1995;21:1495–510.
Jurenka RA, Roelofs WL. Characterization of the acetyltransferase used in pheromone biosynthesis in moths: Specificity for the Z isomer in tortricidae. Insect Biochem. 1989;19(7):639–44.
Haynes KF, Hunt RE. A mutation in pheromonal communication system of cabbage looper moth, Trichoplusia ni. J Chem Ecol. 1990;16:1249–57.
Jurenka RA, Haynes KF, Adlof RO, Bengtsson M, Roelofs WL. Sex pheromone component ratio in the cabbage looper moth altered by a mutation affecting the fatty acid chain-shortening reactions in the pheromone biosynthetic pathway. Insect Biochem Mol Biol. 1994;24:373–81.
Wu WQ, Zhu JW, Millar J, Löfstedt C. A comparative study of sex pheromone biosynthesis in two strains of the turnip moth, Agrotis segetum, producing different ratios of sex pheromone components. Insect Biochem Mol Biol. 1998;28:895–900.
Löfstedt C, Elmfors A, Sjögren M, Wijk E. Confirmation of sex pheromone biosynthesis from (16-D3) palmitic acid in the turnip moth using capillary gas chromatography. Experientia. 1986;42:1059–61.
Löfstedt C, Van Der Pers JNC, Löfqvist J, Lanne BS, Appelgren M, Bergstrom G, et al. Sex pheromone components of the turnip moth, Agrotis segetum: chemical identification, electrophysiological evaluation and behavioral activity. J Chem Ecol. 1982;8:1305–21.
Löfstedt C, Löfquist J, Lanne BS, Van Der Peers JNC, Hansson BS. Pheromone dialects in European turnip moths Agrotis segetum. Oikos. 1986;46:250–7.
Tòth M, Löfstedt C, Blair BW, Cabello T, Farag AI, Hansson BS, et al. Attraction of male turnip moths Agrotis segetum (lepidoptera, noctuidae) to sex-pheromone components and their mixtures at 11 sites in Europe, Asia, and Africa. J Chem Ecol. 1992;18(8):1337–47.
Strandh M, Johansson T, Ahren D, Löfstedt C. Transcriptional analysis of the pheromone gland of the turnip moth, Agrotis segetum (Noctuidae), reveals candidate genes involved in pheromone production. Insect Mol Biol. 2008;17(1):73–85.
Miller JR, Koren S, Sutton G. Assembly algorithms for next-generation sequencing data. Genomics. 2010;95:315–27.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Gu SH, Wu KM, Guo YY, Pickett JA, Field LM, Zhou JJ, et al. Identification of genes expressed in the sex pheromone gland of the black cutworm Agrotis ipsilon with putative roles in sex pheromone biosynthesis and transport. BMC Genomics. 2013;14:636.
Zhu J, Millar J, Löfstedt C. Hormonal regulation of sex pheromone biosynthesis in the turnip moth, Agrotis segetum. Archs Insect Biochem Physiol. 1995;30:41–59.
Choi MY, Fuerst EJ, Rafaeli A, Jurenka R. Identification of a G protein-coupled receptor for pheromone biosynthesis activating neuropeptide from pheromone glands of the moth Helicoverpa zea. Proc Natl Acad Sci U S A. 2003;100:9721–6.
Matsumoto S, Ozawa R, Nagamine T, Kim GH, Uchiumi K, Shono T, et al. Intracellular transduction in the regulation of pheromone biosynthesis of the silkworm, Bombyx mori: suggested involvement of calmodulin and phosphoprotein phosphatase. Biosci Biotechnol Biochem. 1995;59:560–2.
Jurenka RA, Jacquin E, Roelofs WL. Stimulation of pheromone biosynthesis in the moth Helicoverpa zea: Action of a brain hormone on pheromone glands involves Ca2+ and cAMP as second messengers. Proc Natl Acad Sci U S A. 1991;88:8621–5.
Mas E, Llòria J, Quero C, Camps F, Fabriàs G. Control of the biosynthetic pathway of Sesamia nonagrioides sex pheromone by the pheromone biosynthesis activating neuropeptide. Insect Biochem Mol Biol. 2000;30:455–9.
Pape ME, Lopezcasillas F, Kim KH. Physiological regulation of acetyl-CoA carboxylase gene expression: effects of diet, diabetes, and lactation on acetyl-CoA carboxylase mRNA. Arch Biochem Biophys. 1988;267(1):104–9.
Volpe JJ, Vagelos PR. Saturated fatty acid biosynthesis and its regulation. Annu Rev Biochem. 1973;42:21–60.
Tang JD, Charlton RE, Jurenka RA, Wolf WA, Phelan PL, Sreng L, et al. Regulation of Pheromone Biosynthesis by a Brain Hormone in 2 Moth Species. Proc Natl Acad Sci U S A. 1989;86(6):1806–10.
Bjostad LB, Roelofs WL. Biosynthesis of sex pheromone components and glycerolipid precursors from sodium [1–14C] acetate in redbanded leafroller moth. J Chem Ecol. 1984;10:681–91.
Hashimoto K, Yoshizawa AC, Okuda S, Kuma K, Goto S, Kanehisa M. The repertoire of desaturases and elongases reveals fatty acid variations in 56 eukaryotic genomes. J Lipid Res. 2008;49:183–91.
Knipple DC, Rosenfield CL, Nielsen R, You KM, Jeong SE. Evolution of the integral membrane desaturase gene family in moths and flies. Genetics. 2002;162(4):1737–52.
Jeong SE, Rosenfield CL, Marsella-Herrick P, Man You K, Knipple DC. Multiple acyl-CoA desaturase-encoding transcripts in pheromone glands of Helicoverpa assulta, the oriental tobacco budworm. Insect Biochem Mol Biol. 2003;33:609–22.
Vogel H, Heidel AJ, Heckel DG, Groot AT. Transcriptome analysis of the sex pheromone gland of the noctuid moth Heliothis virescens. BMC Genomics. 2010;11:29.
Ikeda Y, Okamura-Ikeda K, Tanaka K. Purification and characterization of short-chain, medium-chain, and long-chain acyl-CoA dehydrogenases from rat liver mitochondria. Isolation of the holo- and apoenzymes and conversion of the apoenzyme to the holoenzyme. J Biol Chem. 1985;260:1311–25.
Izai K, Uchida Y, Orii T, Yamamoto S, Hashimoto T. Novel fatty acid beta-oxidation enzymes in rat liver mitochondria. I. Purification and properties of very-long-chain acyl-coenzyme A dehydrogenase. J Biol Chem. 1992;267:1027–33.
Schulz H. Oxidation of fatty acids in eukaryotes. In: Vance DE, Vance J, editors. Biochemistry of Lipids, Lipoproteins and Membranes. 5th ed. Amsterdam: Elsevier; 2008. p. 131–54.
Kunau WH, Dommes V, Schulz H. beta-oxidation of fatty acids in mitochondria, peroxisomes, and bacteria: a century of continued progress. Prog Lipid Res. 1995;34:267–342.
Uchida Y, Izai K, Orii T, Hashimoto T. Novel fatty acid beta-oxidation enzymes in rat liver mitochondria. II. Purification and properties of enoyl-coenzyme A (CoA) hydratase/3-hydroxyacyl-CoA dehydrogenase/3-ketoacyl-CoA thiolase trifunctional protein. J Biol Chem. 1992;267:1034–41.
Schreurs M, Kuipers F, van der Leij FR. Regulatory enzymes of mitochondrial beta-oxidation as targets for treatment of the metabolic syndrome. Obes Rev. 2010;11:380–8.
Stanley-Samuelson DW, Jurenka RA, Cripps C, Blomquist GJ, de Renobales M. Fatty acids in insects: composition, metabolism, and biological significance. Arch Insect Biochem Physiol. 1988;9:1–33.
Hashimoto T. Peroxisomal β-oxidation: enzymology and molecular biology. Ann N Y Acad Sci. 1996;804:86–98.
Liénard MA, Löfstedt C. Functional flexibility as a prelude to signal diversity?: role of a fatty acyl reductase in moth pheromone evolution. Commun Integr Biol. 2010;3(6):586–8.
Itoh N, Slemmon J, Hawke D. Cloning of Drosophila choline acetyltransferase cDNA. Proc Natl Acad Sci U S A. 1986;83:4081–5.
Günther CS, Chervin C, Marsh KB, Newcomb RD, Souleyre EJ. Characterisation of two alcohol acyltransferases from kiwifruit (Actinidia spp.) reveals distinct substrate preferences. Phytochem. 2011;72:700–10.
Beekwilder J, Alvarez-Huerta M, Neef E, Verstappen FWA, Bouwmeester HJ, Aharoni A. Functional characterization of enzymes forming volatile esters from strawberry and banana. Plant Physiol. 2004;135:1865–78.
Kalscheuer R, Stoveken T, Luftmann H, Malkus U, Reichelt R, Steinbüchel A. Neutral lipid biosynthesis in engineered Escherichia coli: jojoba oil-like wax esters and fatty acid butyl esters. Appl Environ Microbiol. 2006;72:1373–9.
Uthoff S, Stoveken T, Weber N, Vosmann K, Klein E, Kalscheuer R, et al. Thio wax ester biosynthesis utilizing the unspecific bifunctional wax ester synthase/acyl coenzyme A: diacylglycerol acyltransferase of Acinetobacter sp. strain ADP1. Appl Environ Microbiol. 2005;71:790–6.
St-Pierre B, De Luca V. Evolution of acyltransferase genes: origin and diversification of the BAHD superfamily of acyltransferases involved in secondary metabolism. In: Ibrahim R, Varin L, De Luca V, Romeo JT, editors. Recent Advances in Phytochemistry Evolution of Metabolic Pathways. 2000, 34: 285–315
D’Auria J. Acyltransferases in plants: a good time to be BAHD. Curr Opin Plant Biol. 2006;9:331–40.
Fujii T, Ito K, Katsuma S, Nakano R, Shimada T, Ishikawa Y. Molecular and functional characterization of an acetyl-CoA acetyltransferase from the adzuki bean borer moth Ostrinia scapulalis (Lepidoptera: Crambidae). Insect Biochem Mol Biol. 2009;40:74–8.
Souleyre EJF, Greenwood DR, Friel EN, Karunairetnam S, Newcomb RD. An alcohol acyl transferase from apple (cv. Royal Gala), MpAAT1, produces esters involved in apple fruit flavor. FEBS J. 2005;272:3132–44.
Yahyaoui FE, Wongs-Aree C, Latche A, Hackett R, Grierson D, Pech JC. Molecular and biochemical characteristics of a gene encoding an alcohol acyl-transferase involved in the generation of aroma volatile esters during melon ripening. Eur J Biochem. 2002;269:2359–66.
Aharoni A, Keizer LC, Bouwmeester HJ, Sun Z, Alvarez-Huerta M, Verhoeven HA, et al. Identification of the SAAT gene involved in strawberry flavor biogenesis by use of DNA microarrays. Plant Cell. 2000;12:647–62.
Zhu JW, Löfstedt C, Bengtsson BO. Genetic variation in the strongly canalized sex pheromone communication system of the European corn borer Ostrinia nubilalis Hübner (Lepidoptera; Pyralidae). Genetics. 1996;144:757–66.
Ma PWK, Ramaswamy SB. Biology and ultrastructure of sex pheromone-producing tissue. In: Blomquist GJ, Vogt RG, editors. Insect pheromone biochemistry and molecular biology. London: Elsevier; 2003. p. 19–51.
Rosén W. Endogenous control of circadian rhythms of pheromone production in the turnip moth, Agrotis segetum. Arch Insect Biochem Physiol. 2002;50:21–30.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, et al. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19:651–2.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.
Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science. 1997;278:631–7.
Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34:W293–297.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.
Bjostad BL, Roelofs LW. Sex pheromone biosynthetic precursors in Bombyx mori. Insect Biochem. 1984;14:275–8.
The study was supported by grants from Formas and VR to Christer Löfstedt and grants from the Royal Physiographic Society in Lund to Baojian Ding. We thank Åsa Hagström for advice on and assistance with the construction of gene trees, Tomas Johansson for helpful comments on the manuscript and Erling Jirle for technical support. We thank Fredrik Andersson and Erik Hedenström at Mid Sweden University for synthesizing the Z5-10:Me, Z7-12:Me and Z9-14:Me from corresponding alcohols.
The authors declare that they have no competing interests.
BJD and CL conceived the study, BJD and CL designed the study, BJD performed the research, BJD analyzed the data, BJD and CL wrote the manuscript. Both authors read and approved the final manuscript.
About this article
Cite this article
Ding, BJ., Löfstedt, C. Analysis of the agrotis segetum pheromone gland transcriptome in the light of Sex pheromone biosynthesis. BMC Genomics 16, 711 (2015). https://doi.org/10.1186/s12864-015-1909-2
- Agrotis segetum
- Pheromone biosynthesis
- Fatty-acyl reductase
- Functional assay
- Yeast expression