Key biosynthetic gene subfamily recruited for pheromone production prior to the extensive radiation of Lepidoptera
© Liénard et al; licensee BioMed Central Ltd. 2008
Received: 04 February 2008
Accepted: 02 October 2008
Published: 02 October 2008
Moths have evolved highly successful mating systems, relying on species-specific mixtures of sex pheromone components for long-distance mate communication. Acyl-CoA desaturases are key enzymes in the biosynthesis of these compounds and to a large extent they account for the great diversity of pheromone structures in Lepidoptera. A novel desaturase gene subfamily that displays Δ11 catalytic activities has been highlighted to account for most of the unique pheromone signatures of the taxonomically advanced ditrysian species. To assess the mechanisms driving pheromone evolution, information is needed about the signalling machinery of primitive moths. The currant shoot borer, Lampronia capitella, is the sole reported primitive non-ditrysian moth known to use unsaturated fatty-acid derivatives as sex-pheromone. By combining biochemical and molecular approaches we elucidated the biosynthesis paths of its main pheromone component, the (Z,Z)-9,11-tetradecadien-1-ol and bring new insights into the time point of the recruitment of the key Δ11-desaturase gene subfamily in moth pheromone biosynthesis.
The reconstructed evolutionary tree of desaturases evidenced two ditrysian-specific lineages (the Δ11 and Δ9 (18C>16C)) to have orthologs in the primitive moth L. capitella despite being absent in Diptera and other insect genomes. Four acyl-CoA desaturase cDNAs were isolated from the pheromone gland, three of which are related to Δ9-desaturases whereas the fourth cDNA clusters with Δ11-desaturases. We demonstrated that this transcript (Lca-KPVQ) exclusively accounts for both steps of desaturation involved in pheromone biosynthesis. This enzyme possesses a Z11-desaturase activity that allows transforming the palmitate precursor (C16:0) into (Z)-11-hexadecenoic acid and the (Z)-9-tetradecenoic acid into the conjugated intermediate (Z,Z)-9,11-tetradecadienoic acid.
The involvement of a single Z11-desaturase in pheromone biosynthesis of a non-ditrysian moth species, supports that the duplication event leading to the origin of the Lepidoptera-specific Δ11-desaturase gene subfamily took place before radiation of ditrysian moths and their divergence from other heteroneuran lineages. Our findings uncover that this novel class of enzymes affords complex combinations of unique unsaturated fatty acyl-moieties of variable chain-lengths, regio- and stereo-specificities since early in moth history and contributes a notable innovation in the early evolution of moth-pheromones.
An important step in insect pheromone biosynthesis is the involvement of specific desaturases introducing double bonds at specific positions in the fatty-acyl chain by removing two hydrogen atoms [7–10]. Insect desaturases are homologous to the ancestral Δ9 acyl-CoA desaturases of plants, vertebrates and fungi and are functioning as part of a multienzyme complex residing in the endoplasmic reticulum (ER) [11–14]. Certain moth desaturase lineages are believed to have arisen subsequently to the divergence of Lepidoptera and Diptera that took place around 300 million years ago (MYA) [10, 15] and may have played a significant role in the evolution of sex-pheromones in Lepidoptera. Interestingly, the biosynthesis of many ditrysian oxygenated pheromones proceeds with various Δ11-desaturation reactions [7, 15, 16]. Up to this date desaturase-encoding genes have been investigated in a dozen ditrysian moths, depicting examples of elaborate evolution in which a minimal number of enzymes account for the chemical diversity and species-specificity of pheromone components found among different species [17–32].
Only a few primitive moth species have so far been investigated considering their chemical communication. Their pheromones are made of short-chain alcohols and ketones and release through sternal globular glands whose openings are located on the Vth abdominal segment [33–35]. The modern distal female abdominal gland and the typical ditrysian calling postures likely arose prior to divergence of the Ditrysia as a few examples are known in some monotrysian Heteroneura [reviewed in  and ]. However, the exact localization of distal gland producing-pheromone in non-ditrysian Heteroneura remains to be determined. Likewise, the origin of the ditrysian sex-pheromone, i.e., the emergence of mating-signals derived from saturated long-chain fatty acids and involving genes (e.g., desaturases) specific for pheromone production, remains obscure although it has been suggested to originate in a moth lineage prior to the divergence of Ditrysia and the other Heteroneura lineages [36, 38] (Fig. 1). Recently, long-chain unsaturated pheromone components were identified from the abdominal tip in a moth species outside the Ditrysia; the sex pheromone of the currant shoot borer, Lampronia capitella (Incurvarioidea: Prodoxidae), was characterized as a mixture of (Z,Z)-9,11-tetradecadienol, and the corresponding acetate and aldehyde . This monotrysian species constitutes the most primitive moth for which long-chain unsaturated pheromone components have been evidenced and provides a unique opening for investigating evolutionary aspects of the ditrysian pheromone mating-signals.
In this study we report on pheromone-gland precursor identification, molecular characterization and functional expression of desaturase transcripts from the pheromone gland of L. capitella, of which a functional Δ11-desaturase transcript that exhibits all of the required biochemical activities to biosynthesize the conjugated chemical structures of the major sex pheromone component. Our findings evidence that the novel desaturase lineages, (Δ9 (18C>16C) and Δ11) evolved before the split between the ditrysian and other heteroneuran moth lineages. By extension, our findings suggest that the specific biosynthetic functions inherited from early moth history and which contributed to structural variations in mating-signals may have played a role in the radiation of the higher Lepidoptera.
Fatty-acid pheromone precursor identification
The double-bond positions of the methyl tetradecadienoate were identified as 9,11 by analysis of its MTAD derivatives that exhibited the diagnostic ions at m/z 351 (M+), 194 and 322. The presence of a methyl 11,13-hexadecadienoate was also found (data not shown), exhibiting diagnostic MTAD adducts at m/z 379 (M+), 194 and 350, which can be rationalised to be formed by chain-elongation of the relatively abundant Z9,Z11–14:acyl.
When the glands were incubated with D9-Z11–16:acid for 24 hours in vivo, the label was incorporated into Z9–14:Me (data not shown), a potential intermediate in pheromone biosynthesis, thereby indicating that Z9–14:acyl could be produced by β-oxidation of Z11–16:acyl.
Characterization of desaturase transcripts
Total RNA was isolated from female abdominal tips of L. capitella. In PCR reactions using the corresponding cDNA as template and primers targeting conserved desaturase motifs , 550-bp DNA fragments were amplified that encompassed the central region of a desaturase gene. From the DNA information provided by the central region, four distinct transcripts were identified and their full-length cDNA sequences subsequently cloned.
The first desaturase transcript, Lca-QPAQ spans 1,427 bp and encompasses an ORF encoding a protein of 353 aa residues, which shares high aa sequence similarity with known Z9-desaturases (16C>18C), e.g., 81% with the Z9-desaturase from Choristoneura parallela [GenBank:AAN39701] and 80% with the Z9-desaturase from Epiphyas postvittana [GenBank:AAL35750].
The second desaturase transcript, Lca-SPVE spans 1,953 bp and contains an ORF encoding a 352-aa residue protein. This deduced protein shows 64% sequence similarity to Lca-QPAQ but higher similarity to other Z9-desaturases (18C>16C), e.g., 77% aa-sequence similarity to the Z9-desaturase from Trichoplusia ni [GenBank:AAB92583] and 76% to the Z9-desaturase from Spodoptera littoralis [GenBank:AAQ74257].
The third desaturase transcript, Lca-GATD spans 1,852 bp and encompasses an ORF encoding a 380-aa residue protein, which displays 51% sequence similarity to both Lca-QPAQ and Lca-SPVE and even higher similarity (73%) to the Z9-desaturase (C14-C26) from C. parallela [GenBank:AAQ12887].
Functional assay of Δ9-desaturases by complementation in mutant yeast using the expression vector YEpOLEX
A desaturase-deficient yeast cell line (ole1) was transformed with YEpOLEX plasmids containing the Lca-QPAQ, Lca-SPVE or Lca-GATD ORF, respectively. YEpOLEX-Lca-QPAQ transformants were able to grow on medium lacking supplemental unsaturated fatty-acids (UFAs), indicating that the Lca-QPAQ cDNA was encoding a desaturase that complemented the UFA auxotrophic ole1 strain. In contrast, YEpOLEX-Lca-SPVE and YEpOLEX-Lca-GATD transformants were unable to grow on media lacking UFAs. These results were confirmed using distinct verified gene constructs under identical experimental conditions. Chromatograms of methylated fatty-acid extracts from the Lca-QPAQ transformants showed three peaks with retention times corresponding to Z9–14:Me, Z9–16:Me and Z9–18:Me, in a 2:62:36 ratio. Fatty-acid methyl esters from yeast transformed with the Z9-desaturase gene of H. assulta (YEpOLEX-Hass-KPSE, ) were also prepared for use as a positive control. In this case the transformation resulted in the production of Z9–16:Me and Z9–18:Me, with a preference for palmitic acid, as shown in GC-MS analyses (data not presented). Thus Lca-QPAQ and Hass-KPSE are encoding Z9-desaturases displaying a substrate preference for palmitic acid, which is consistent as these desaturases are phylogenetically closely related (Δ9 16C>18C cluster) (Fig. 4). Double-bond position in the reaction products was confirmed by analyses of DMDS derivatives, which exhibited the characteristic fragment at m/z 217 at expected retention times (data not shown). No distinctive products were detected in FAME or MTAD analyses after incubation of yeasts transformed with any of the YEpOLEX-Lca-constructs in the presence of Z11–14:Me.
Functional assay of a Δ11-desaturase by complementation in mutant yeast using the expression vector pYEX-CHT
Desaturases as key factor for pheromone evolution
There is no single explanation for the evolutionary and ecological success of Lepidoptera but it is rather thought to be the result of a cascade of successful evolutionary innovations . The recruitment of certain lineages of desaturase genes to serve in pheromone biosynthesis might be one of these innovations that have played a significant role in the evolution of Lepidoptera, setting the scene for the development of a very effective mate-communication system [10, 38]. The radiation of the Lepidoptera took off about 125 MYA, approximately at the time when the Ditrysia diverged from other heteroneuran Lepidoptera . This evolutionary event appears to coincide with the emergence of a novel pheromone communication system, characterized by the involvement of novel desaturases in the production of long-chain unsaturated fatty-acid derivatives as pheromone components. In particular, Δ11-desaturases contributed significantly to the biosynthesis of a broad range of unsaturated pheromone components in ditrysian species and strongly increased the structural diversity of pheromone precursors [6–8]. In this study we demonstrate the involvement of a Δ11-desaturase in the pheromone biosynthesis in a non-ditrysian moth, Lampronia capitella (Incurvarioidea, Prodoxidae), which strongly supports the interpretation that the recruitment of the specialized desaturases to serve in the production of moth pheromone components occurred in the early evolution of Lepidoptera, before their extensive radiation.
Involvement of a Δ11-desaturase gene in mate signalling in the primitive monotrysian moth L. capitella
Four full-length desaturase encoding cDNAs were characterized from the L. capitella abdominal tip and their deduced amino acid sequences notably revealed the three histidine-rich motifs (HIS boxes) that are characteristic of membrane-bound desaturases . The Lca-QPAQ, Lca-SPVE and Lca-GATD transcripts shared high homologies with known Δ9-desaturases while the fourth transcript (Lca-KPVQ) showed a high amino acid sequence similarity with several Δ10, Δ11 and bifunctional desaturases from ditrysian moth species. Sequence analysis evidenced that the isolated transcript represented a typical fatty-acyl desaturase gene. Its predicted primary protein structure (Fig. 3) shared several key features with other known insect desaturases, like the position and the length of the transmembrane domains relative to the conserved HIS boxes and the xxxQ signature motif according to a desaturase-nomenclature previously proposed . A reconstruction of the evolutionary history of desaturase genes (Fig. 4) finally indicated this candidate desaturase to cluster with members of the Δ11-desaturase lineage.
Functional characterization of the Lca-KPVQ gene product in a desaturase- and elongase-deficient yeast strain was performed and analyses of the unsaturated fatty-acid profiles of pYEX-Lca-KPVQ-expressing transformants established that the gene product displayed a Z11-desaturase activity with relatively broad substrate specificity. Saturated substrates from C8 to C22 occur naturally in the ole 1 elo1 yeast strain and the gene-product catalysed the removal of two hydrogen atoms from C14, C16, C18 and C20, with a strong substrate preference for palmitic acid (C16) as shown in Fig. 5A. This enzyme is also the first insect Δ11-desaturase reported to act on the C20:acyl substrate. After complete GC-MS analysis of both yeast methanolyzed extracts and their DMDS derivatives, the Z configuration of the double bond was confirmed in all newly formed unsaturated substrates and no E isomers were detected. Two Δ11-desaturases (Ostrinia sp., GenBank acc nos. AAL32060 and AAL35331) were characterized that produced Z-isomers of C16 and C18 acyl precursors in addition to Z and E isomers of C14 acyl precursors that were later converted into the active pheromone compounds . Though trans-unsaturated pheromone components are not utilized by L. capitella, minor amounts of E11–14:Me were detected in DMDS analyses of FAMEs from female abdominal tips. This was not confirmed in vitro, which might be due to the limits of detection and does not alter the above evidences of an active Δ11-gene.
After supplementation with excess of Z9–14:Me, the yeast produced small amounts of the conjugated Z9,Z11-tetradecadienoate, which was detected upon GC-MS analyses of the cis MTAD-adduct (Fig. 5B). In presence of large amounts of both pheromone biosynthetic intermediates (i.e., the 16:acyl and the Z9–14:acyl), the Lca-KPVQ transcript always preferentially catalyzed the formation of the Z11–16:acyl monounsaturated intermediate.
Δ9-stearyl-CoA desaturases are primordial enzymes regulating the level of unsaturated fatty-acid biosynthesis and thereby allowing living organisms to maintain the physical structure and fluidity of membrane lipid bilayers [42–44]. When transforming yeasts with YEpOLEX-Lca-SPVE or YEpOLEX-Lca-GATD constructs, the Δ9 desaturase activity of the ole1 yeast strain could not be restored. Δ9-desaturases (18C>16C) exist that catalyse the formation of Z9-monoenes in the pheromone gland of some ditrysian moth species [19–21, 23, 25]; although in most cases pheromone biosynthesis does not involve Δ9 desaturation . It is also not uncommon to encounter inactive desaturase transcripts as exemplified by two Lca-GATD paralogs, the HassGATD transcript that was found to be non-functional in pheromone biosynthesis in H. assulta  or the CpaZ9GATD (C14-C26) transcript that produced a long series of monounsaturated fatty acids in C. parallela, an activity that could however not be linked to sex-pheromone biosynthesis . It might then be envisioned that both Lca-SPVE and Lca-GATD are catalytically inactive desaturases in this moth species and as previously suggested, could represent pseudogenes or encode other functions yet to be determined [10, 24, 28]. Analyses of esterified total-lipid extracts of YEpOLEX-Lca-QPAQ-transformed yeast revealed that the gene product was producing a series of Z9-monoenes with a preference for palmitic acid like other Δ9 16C>18C insect-desaturases (Fig. 4). The Z9–14:acyl, intermediate in the pheromone biosynthesis (pathway b, Fig. 6), accounted for 2% of the total amount of UFAs produced and could thus be formed by the Lca-Z9-QPAQ desaturase. However, topical application of labelled precursors indicated that the D9-Z11–16:acyl is also significantly incorporated into Z9–14:acyl. Knowing that the Z11–16:acyl is produced by Δ11-desaturation and that significant chain-shortening occurs in the gland, it is more likely that the Z9–14:acyl derives from β-oxidation of the Z11–16:acyl precursor. The Lca-Z9-QPAQ desaturase is most likely not involved in the pathway bof pheromone production. Besides, no conjugated Z9,Z11–14:acyl could be detected from Lca-Z9-QPAQ yeasts transformants grown in presence of Z11–14:Me and this transcript is therefore not involved in the Δ9-desaturation reaction proposed in the pathway c(Fig. 6).
In conclusion, our in vitro analyses indicate that the Lca-Z11-KPVQ transcript encodes a functional protein displaying enzymatic properties consistent with the proposed pheromone biosynthetic pathway of L. capitella (pathway a, Fig. 6). While the first Z11-desaturation step on palmitic acid is a prevalent reaction that has been characterized in several other ditrysian moth species (see e.g., [17, 26, 27]), the Δ11-desaturase activity on monoene precursors has been evidenced in only a few moth species up to date. In E. postvittana, an E11-desaturase catalyses the formation of E9,E11–14:acid from E9–14:acid  while in B. mori, Manduca sexta and S. littoralis bifunctional Δ11 and Δ10,12-desaturases catalyse the formation of conjugated Δ10,12 dienes from Z11-monounsaturated precursors [27, 29, 31].
Insights into the origin of the Δ11-desaturase subfamily
Any gene in any genome is selectively constrained and most mutations that change the fitness of an organism are expected to be deleterious [45–47]. Gene duplications represent major opportunities that contribute to functional novelties in all living organisms thereby playing a vital role during evolution [47, 48] and the predominant motive during the evolution of a novel gene function is to gain a selective advantage . The need for efficient mate signalling and species-specific recognition might have been a strong motive for establishment of new functions within pheromone communication as derived from duplicated genes. The pheromone-desaturase family most likely originated in a common ancestor of Diptera and Lepidoptera before their divergence in the early Carboniferous, i.e., between 330 and 350 MYA . Several extant flies hence use orthologs of the Δ9 (16C>18C) group (Fig. 4)-which represents the metabolic ancestral function of the gene family [10, 15]- for synthesizing cuticular-hydrocarbon sex-pheromones [9, 50]. Subsequently the lepidopteran desaturase gene family evolved under a birth-and-death evolutionary process and underwent several duplication events leading to five well-supported clades [10, 21, 51] (Fig. 4). Reconstructing the evolutionary history of desaturase genes shows that only the Δ11 and the Δ9 (18C>16C) lineages have no orthologous genes in any Dipteran genome examined to date (Fig. 4), which suggests that the genes were either lost in Diptera subsequent to the divergence of moths and flies or that both lineages may have been recruited following a duplication point succeeding the divergence between Diptera and Lepidoptera [10, 21]. We have examined the evolution of desaturases using predicted genes available from other insect orders including Hymenoptera, Hemiptera and Coleoptera (see additional File 1) and we found no evidence for orthologous genes of the Δ11 and Δ9 (18C>16C lineages) in these insect orders. Both desaturase lineages are thus most likely innovations of lepidopterans, of which the Δ11-desaturase lineage exclusively serves in the pheromone production. Our molecular and functional investigations of pheromone biosynthetic genes in L. capitella strongly indicate that the duplication events that gave rise to both desaturase-gene lineages took place before the divergence of Heteroneura lineages from Ditrysia, around 125 MYA. The alternative explanation that both duplications occurred independently in two different lineages is a less parsimonious explanation.
Clearly, the emergence and establishment of the Δ11-desaturase subfamily in moth genomes might have promoted structural diversity among pheromones early in the evolution of Lepidoptera. Whether the duplication events took place before or after the divergence of Lepidoptera and its sister group Trichoptera will be interesting to investigate. In this context it would be of interest to find out whether Δ9- or Δ11-desaturases (or another desaturase function yet unidentified) underlie the double-bond formation in the biosynthesis of (Z)-4-hepten-2-ol and (Z)-6-nonen-2-ol (and their corresponding ketones) as observed in Eriocraniidae moths and in Trichoptera. Both compounds could hypothetically be produced by Δ11-desaturation of tetradecanoate, followed by chain shortening by successive steps of β-oxidation and finally decarboxylation . Investigating pheromone biosynthesis in these more primitive non-ditrysian moth lineages will further highlight to what extent the recruitment of discrete pheromone production genes contributed to the emergence of the typical ditrysian moth signalling chemistry and by extension promoted the adaptive success of the Lepidoptera.
Likewise, the exact location of the distal pheromone-producing cells and their morphological features in Monotrysia should also be investigated. Such an investigation could resolve the intriguing issue whether the pheromone-producing desaturases in monotrysian moths are catalysing the formation of biosynthetic precursors in a pheromone gland structure alike those of most ditrysian species or within -a yet uncharacterized internal structure- of the abdominal tip.
The renowned sexual communication of Lepidoptera has likely evolved through a combination of mechanisms [2, 36] involving the diversification of female-produced chemical signals that mediate mate attraction and corresponding changes in male responses to these signals . Variation in pheromone-biosynthetic genes have been supported for promoting changes in the emitted signals thereby indirectly directing major adaptations in moth mate-recognition systems . A core component of the evolution of moth chemistry has been evidenced that imply gene duplications of biosynthetic genes and their maintenance in the insect genomes [10, 15]. We here brought evidence that Lepidoptera-specific pheromone-production desaturase genes evolved before radiation of the ditrysian lineages. Especially, a Δ11-desaturase (Lca-KPVQ) was shown catalyzing key reactions leading to mono- and di-unsaturated fatty acyl-moieties in the primitive moth Lampronia capitella, the sole non-ditrysian species identified to date that relies on long-chain di-unsaturated fatty-acid derivatives for mate-attraction. Like many other characterized ditrysian-moth desaturases, Lca-KPVQ catalyses a few more desaturation reactions in vitro than those leading to the biologically relevant components used by L. capitella. This supports the current view that subtle alterations in enzymatic activities at different levels in the biosynthesis may have lead to subtle shifts in pheromones and ultimately to new species-specific communication channels [15, 21]. Recruiting Δ11-desaturases in mate signalling prior to the extensive lepidopteran radiation may thus have contributed, as part of a complex biological framework, to the evolution of the typical moth pheromone chemistry and supports the hypothesis that a limited number of biosynthetic genes played a pivotal role in ditrysian moth evolution.
Immature larvae were collected in a black currant (Ribes nigrum) orchard in mid-May in Sörfors, northwest of Umeå (Västerbotten, Sweden). Infested branches were cut, brought to the laboratory and kept at 15°C until the larvae had pupated. Pupae were separated from branches or folded leaves and were maintained at 23°C in 60% relative humidity and in a 17:7 hours light:dark photoperiod.
Application of labelled precursors
Virgin 0 to 2 day-old female moths were used for topical application of labelled precursors onto the abdominal tip. In the early photophase, individual calling females were anaesthetized using carbon dioxide and the abdomen was then gently squeezed in order to completely expose the abdominal tip. D9-Z11–16:acid (4 μg in 0.2 μl DMSO) was topically applied onto the abdominal tip and the solution was allowed to be absorbed by the gland during a few minutes. Females were reintroduced into individual cages for 24 hours. Abdominal tips from 10 females (including ovipositor) were then carefully dissected under a stereomicroscope and extraction of pheromone components was performed for 1 hour in a glass capillary containing 80 μl hexane.
For analysis of fatty-acid methyl esters, a total lipid extraction was performed using chloroform:methanol (2:1 v/v). Base methanolysis was then utilized to convert fatty-acyl moieties into the corresponding methyl esters. Abdominal tip tissues were removed and the solvent was evaporated under a gentle N2 stream. Concentrated extracts were treated with 100 μl 0.5 M KOH/methanol and allowed to react for 1 hour at room temperature. For base methanolysis of total lipid yeast extracts, yeast residues were directly treated with 500 μl 0.5 M KOH/methanol after solvent evaporation . Abdominal tip samples or yeast extracts were acidified by addition of 100 or 500 μl 1.0 M HCl, respectively and the resulting fatty-acid methyl esters (FAME) were collected in hexane, and analysed by GC-MS analysis.
Determination of double bond positions by DMDS and MTAD
Methanolyzed extracts were converted into dimethyl-disulfide (DMDS) derivatives by addition of 50 μl DMDS and 5 μl iodine 5% in diethyl ether and then incubated over night at 40°C. Two hundred μl of hexane were added to the sample and the reaction was neutralized by addition of 20 to 50 μl Na2S2O3 5% . The organic phase was removed, concentrated under a gentle N2 stream and subjected to GC-MS analysis. The 4-methyl-1,2,4-triazoline-3,5-dione (MTAD) adducts were prepared by transferring 5 μl of the methyl ester extracts into a glass vial containing 10 μl CH2Cl2 and treating the resulting solution with 10 μl of a MTAD solution (2 μg/μl; in CH2Cl2). Reactions were incubated for 15 min at room temperature and 2 μl were subjected to GC-MS analysis .
Gas chromatography-mass spectrometry analyses
Samples were analysed on a gas chromatograph (GC) (Hewlett Packard HP 6890 GC system) equipped with an auto-injector (HP 7683) and coupled to a mass selective detector (HP 5973). The GC was equipped with an HP1-MS column (100% methyl siloxane; 30 m 3 0.25 mm, df: 0.25 μm) and helium was used as carrier gas (velocity: 32 cm/s). For FAME and DMDS analyses, the oven temperature was set to 55°C (or 80°C) for 2 min, then increased by 10°C/min up to 250°C, followed by a hold at 250°C for 10 min, and then increased by 20°C/min up to 300°C followed by a hold at 300°C for 5 min. For MTAD adduct analyses, the injector temperature was set to 300°C, the oven temperature was set to 100°C and then increased by 15°C/min up to 300°C followed by a hold at 300°C for 20 min according to procedures previously described .
(11Z)-Tetradeca-11-en-1-ol and (9Z,11E)- and (9Z,11Z)-tetradeca-9,11-dien-1-ols (12–115 mg), respectively, were oxidised at 0°C with Jones reagent according to  followed by acid-catalysed esterification of the resulting acid in methanol at room temperature.
Total RNA isolation and cDNA synthesis
Abdominal tips were carefully dissected from 50 virgin female moths on the first day after their emergence, immediately snap frozen in liquid nitrogen and stored at -80°C. Total RNA was isolated and purified from abdominal tips using the Trizol® reagent (Invitrogen™ Life technologies) according to recommended procedures. One μg of total RNA was reverse transcribed into cDNA using Stratascript (Stratagene) following the manufacturer's protocol.
PCR screening and cloning of pheromone desaturase-encoding cDNAs
Oligonucleotide primers used for desaturase gene characterization and functional expression
Gene specific b
Full-length amplification b
Primer sequence c
Primer sequence c
DNA sequences were analysed using Sequencher V.3.0 software (Gene Codes Corporation, MI, USA) and the Bioedit Sequence alignment editor program V. 5.0.9 .
Sequence analyses and Neighbour-Joining tree construction
Comparison of sequence information to publicly available DNA information was performed using BLAST . Multiple sequence alignments were performed using the Clustal W 18.2 algorithm  and edited in BOXSHADE . Prediction of the endoplasmic retention signal was performed using PSORT . Desaturase aa sequences used for phylogenetic reconstructions were retrieved from the GenBank non-redundant (nr) protein database  and the Neighbour-Joining tree was constructed using MEGA version 3.1 .
Functional assay of Δ9-desaturases by complementation in mutant yeast using the expression vector YEpOLEX
Gene-specific primers including XbaI and SacI restriction sites (Table 1) were designed based on the ORFs identified for each of the three Δ9-like desaturase transcripts. After PCR amplification from abdominal tip cDNA, PCR products were purified by standard procedures followed by trimming by double restriction digestion using XbaI and SacI. Digestion products were purified by gel electrophoresis followed by gel elution and then directionally inserted by replacement cloning in the linearized YEpOLEX plasmid . After verification by DNA sequencing, the final recombinant plasmids were designated YEpOLEX-Lca-QPAQ, YEpOLEX-Lca-SPVE and YEpOLEX-Lca-GATD and used for transformation of a desaturase-deficient yeast strain (MATα ole1::LEU2 leu2–3 leu2–112 trp1-1 ura3–52 his4) . The YEpOLEX-Hass-KPSE plasmid construct from Helicoverpa assulta  was used as a positive control. For selection of uracil prototrophs, transformed yeasts were plated on selective medium containing 0.7% YNB and a drop-out medium lacking uracil and leucine (ForMedium™), 2% glucose, 0.5 mM palmitoleic acid, 0.5 mM oleic acid (Larodan, Sweden), 1% tergitol (type Nonidet NP-40, Sigma) and 0.01% adenine (Sigma). Individual uracil prototroph colonies were transferred to YPAD plates for complementation of the UFA auxotrophy of the ole1 yeast strain used. UFA prototrophs were inoculated in 10 ml YPAD medium with or without Z11–14:Me and incubated at 30°C and 300 rpm for 48 hours. Cells were collected by centrifugation and cell pellets were washed with sterile water three times and then processed for FAME analyses. The double-bond positions of FAMEs were determined by GC-MS analysis of DMDS adducts.
Functional assay of a Δ11-desaturase by complementation in mutant yeast using the expression vector pYEX-CHT
The pYEX-CHT expression vector  was used for functional assay in combination with the double deficient ole1 elo1 S. cerevisiae strain (MATa elo1::HIS3 ole1::LEU2 ade2 his3 leu2 ura3), defective in both desaturase and elongase gene functions . This expression vector is derived from the pYEX plasmid (ClonTech), and is related to the vector used for assaying the biochemical activity of the Δ11-desaturase in S. littoralis and Thaumetopoea pityocampa [29, 30]. Two gene-specific primers, Lc4ORFs and Lc4ORFas (Table 1) encompassing the BamH1 and EcoR1 restriction sites were designed to amplify Lca-KPVQ ORF using the Advantage®2 PCR enzyme system (Clontech), which was then ligated into the TOPO® TA PCR 2.1 vector (Invitrogen) and the construct was transformed into TOPO 10 cells. The insert was amplified, purified, released from plasmid by BamH1 and EcoR1 restriction digestion and subcloned into the BamH1 and EcoR1 cloning sites of the linearized pYEX vector. Constructs were transformed into TOPO 10 cells and analysed by sequencing. The final construct designated pYEX-Lca-KPVQ was transformed for functional expression into the ole1 elo1 yeast. pYEX-Lca-KPVQ yeast transformants were grown on selective medium as described above for 4 days at 30°C. Individual transformants were then isolated and inoculated in 10 ml YPAD medium containing 1 mM CuSO4 and grown for 48 hours at 30°C and 300 rpm. Cells were harvested by centrifugation and then suspended in fresh YPAD medium containing 1 mM CuSO4 with or without addition of Z9–14:Me (Larodan, Sweden). The cells were subsequently incubated for 48 hours, collected by centrifugation and then washed with water. Fatty-acid methyl esters were extracted and analysed by GC-MS. The position of double bonds in monoenes and dienes was determined by DMDS and MTAD derivatizations, followed by GC-MS analysis.
The cDNA-sequence information reported in this study has been deposited and are available by the GenBank accession numbers EU152332 to EU152335.
gas chromatography coupled to mass spectrometry
fatty-acid methyl ester
million years ago
open reading frame
polymerase chain reaction
rapid amplification of cDNA ends
total ion current
unsaturated fatty acid
yeast nitrogen base
acid: (Z)-9-tetradecenoic acid
acid: (Z)-11-tetradecenoic acid
acid: (E)-11-tetradecenoic acid
acid: (Z)-9-hexadecenoic acid (palmitoleic acid)
acid: (Z)-11-hexadecenoic acid
acid: (Z)-9-octadecenoic acid (oleic acid)
acid: (Z)-11-octadecenoic acid
Z11–14:acid: (Z,Z)-9,11-tetradecadienoic acid
E11–14:acid: (Z,E)-9,11-tetradecadienoic acid (fatty acyls and FAMEs are named correspondingly).
We gratefully acknowledge Sven Hellqvist and Martin Barrlund for collecting Lampronia capitella larvae and Erling Jirle for technical assistance; Charles Martin for providing the ole1 yeast strain; Roger Schneiter for providing the ole1 elo1 yeast strain; Joanne Caine for providing the pYEX-CHT vector; Wendell Roelofs for stimulating discussions and Douglas Knipple for providing the YEpOLEX plasmid; Jure Piškur for supervision and valuable advice in heterologous expression; Aleš Svatoš for advice regarding the MTAD reaction; Jean-Marc Lassance for valuable advice in GC-MS analyses and comments on the manuscript and Fredrik Andersson and Anna Nilsson for contributing to methyl esters synthesis. Three anonymous referees are acknowledged for valuable comments on a previous version of this manuscript. This work was supported by the Swedish Research Council (VR) (for CL) and EU (Objective 1 the region of South Forest Countries) and Länsstyrelsen i Västernorrlands län (for EH).
- Scoble M: The Lepidoptera: Form, function and diversity. 1992, Oxford University Press, OxfordGoogle Scholar
- Grimaldi D, Engel MS: Evolution of the insects. 2005, Cambridge University PressGoogle Scholar
- Roelofs WL: Chemistry of sex attraction. Proc Natl Acad Sci USA. 1995, 92: 44-49. 10.1073/pnas.92.1.44.PubMed CentralView ArticlePubMedGoogle Scholar
- El-Sayed AM: The pherobase: Database of Insect pheromones and semiochemicals. 2007, [http://www.pherobase.com/]Google Scholar
- Percy-Cunningham JE, MacDonald JA: Biology and ultrastructure of sex pheromone producing glands. Pheromone biochemistry. 1987, 27-76.Google Scholar
- Bjostad LB, Roelofs WL: Sex pheromone biosynthesis in Trichoplusia ni: key steps involve delta-11 desaturation and chain-shortening. Science. 1983, 220: 1387-1389. 10.1126/science.220.4604.1387.View ArticlePubMedGoogle Scholar
- Bjostad LB, Wolf W, Roelofs WL: Pheromone biosynthesis in lepidopterans: desaturation and chain shortening. Pheromone Biochemistry. Edited by: Blomquist GJ, Prestwich GD. 1987, Academic Press, New-York, 77-120.Google Scholar
- Jurenka RA, Roelofs WL: Biosynthesis and endocrine regulation of fatty acid derived pheromones in moths. Insect lipids: chemistry, biochemistry and biology. Edited by: Stanley-Samuelson DW, Nelson DR. 1993, University of Nebraska Press, Lincoln, Nebraska, 353-388.Google Scholar
- Tillman JA, Seybold SJ, Jurenka RA, Blomquist GJ: Insect pheromones-an overview of biosynthesis and endocrine regulation. Insect Biochem Molec Biol. 1999, 29: 481-514. 10.1016/S0965-1748(99)00016-8.View ArticleGoogle Scholar
- Roelofs WL, Rooney AP: Molecular genetics and evolution of pheromone biosynthesis in Lepidoptera. Proc Natl Acad Sci USA. 2003, 100: 9179-9184. 10.1073/pnas.1233767100.PubMed CentralView ArticlePubMedGoogle Scholar
- Harwood JL: Plant acyl lipids: structure, distribution and analysis. The biochemistry of plants. Edited by: Stumpf PK, Conn EE. 1980, Academic Press, New-York, 1-55.Google Scholar
- Stuckey JE, McDonough VM, Martin CE: The OLE1 gene of Saccharomyces cerevisiae encodes the Δ9 fatty acid desaturase and can be functionally replaced by the rat stearoyl-CoA desaturase gene. J Biol Chem. 1990, 265: 20144-20149.Google Scholar
- Shanklin J, Whittle E, Fox BG: Eight histidine residues are catalytically essential in a membrane-associated iron enzyme, stearoyl-CoA desaturase, and are conserved in alkane hydroxylase and xylene monooxygenase. Biochemistry. 1994, 33: 12787-12794. 10.1021/bi00209a009.View ArticlePubMedGoogle Scholar
- Shanklin J, Cahoon EB: Desaturation and related modifications of fatty acids. Annu Rev Plant Physiol Mol Biol. 1998, 49: 611-641. 10.1146/annurev.arplant.49.1.611.View ArticleGoogle Scholar
- Knipple DC, Rosenfield C-L, Nielsen R, You KM, Jeong SE: Evolution of the integral membrane desaturase gene family in moths and flies. Genetics. 2002, 162: 1737-1752.PubMed CentralPubMedGoogle Scholar
- Roelofs WL, Bjostad L: Biosynthesis of Lepidopteran pheromones. Bioorg Chem. 1984, 12: 279-298. 10.1016/0045-2068(84)90011-7.View ArticleGoogle Scholar
- Knipple DC, Rosenfield C-L, Miller SJ, Liu W, Tang J, Ma PWK, Roelofs WL: Cloning and functional expression of a cDNA encoding a pheromone gland-specific acyl-CoA Delta11-desaturase of the cabbage looper moth, Trichoplusia ni. Proc Natl Acad Sci USA. 1998, 95: 15287-15292. 10.1073/pnas.95.26.15287.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenfield C-L, You KM, Knipple DC: Structural and functional conservation and divergence among acyl-CoA desaturases of two noctuid species, the corn earworm, Helicoverpa zea, and the cabbage looper, Trichoplusia ni. Insect Biochem Molec Biol. 2001, 31: 949-964. 10.1016/S0965-1748(01)00043-1.View ArticleGoogle Scholar
- 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 Molec Biol. 2001, 32: 961-966. 10.1016/S0965-1748(01)00176-X.View ArticleGoogle Scholar
- Hao G, O'Connor M, Liu W, Roelofs WL: Characterization of Z/E11- and Z9-desaturases from the obliquebanded leafroller moth, Choristoneura rosaceana. J Ins Sci. 2002, 2: 26-32.Google Scholar
- Roelofs WL, Liu W, Hao G, Jiao H, Rooney AP, Linn CEJr: Evolution of moth sex pheromones via ancestral genes. Proc Natl Acad Sci USA. 2002, 99: 13621-26. 10.1073/pnas.152445399.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu W, Ma PWK, Marsella-Herrick P, Rosenfield C-L, Knipple DC, Roelofs WL: Cloning and functional expression of a cDNA encoding a metabolic acyl-CoA Δ9-desaturase of the cabbage looper moth, Trichoplusia ni. Insect Biochem Molec Biol. 1999, 29: 435-443. 10.1016/S0965-1748(99)00020-X.View ArticleGoogle Scholar
- Liu W, Jiao H, O'Connor M, Roelofs WL: Moth desaturase characterized that produces both Z and E isomers of Δ11-tetradecenoic acids. Insect Biochem Molec Biol. 2002, 32: 1489-1495. 10.1016/S0965-1748(02)00069-3.View ArticleGoogle Scholar
- Liu W, Jiao H, Murray NC, O'Connor M, Roelofs WL: Gene characterized for membrane desaturase that produces (E)-11 isomers of mono- and diunsaturated fatty acids. Proc Natl Acad Sci USA. 2002, 99: 620-624. 10.1073/pnas.221601498.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu W, Rooney AP, Xue B, 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-311. 10.1016/j.gene.2004.08.017.View ArticlePubMedGoogle Scholar
- Jeong SE, Rosenfield C-L, Marsella-Herrick P, You KM, Knipple DC: Multiple acyl-CoA desaturase-encoding transcripts in pheromone glands of Helicoverpa assulta, the oriental tobacco budworm. Insect Biochem Molec Biol. 2003, 33: 609-622. 10.1016/S0965-1748(03)00043-2.View ArticleGoogle Scholar
- Moto K, Suzuki MG, Hull JJ, Kurata R, Takahashi S, Yamamoto M, Okano K, Imai K, Ando T, Matsumoto S: Involvement of a bifunctional fatty-acyl desaturase in the biosynthesis of the silkmoth, Bombyx mori, sex pheromone. Proc Natl Acad Sci USA. 2004, 101: 8631-8636. 10.1073/pnas.0402056101.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodríguez S, Hao G, Liu W, Piña B, Rooney AP, Camps F, Roelofs WL, Fabriàs G: Expression and evolution of Δ9 and Δ11 desaturase genes in the moth Spodoptera littoralis. Insect Biochem Molec Biol. 2004, 34: 1315-1328. 10.1016/j.ibmb.2004.09.003.View ArticleGoogle Scholar
- Serra M, Piña B, Bujons J, Camps F, Fabriàs G: Biosynthesis of 10,12-dienoic fatty acids by a bifunctional Δ11 desaturase in Spodoptera littoralis. Insect Biochem Molec Biol. 2006, 36: 634-641. 10.1016/j.ibmb.2006.05.005.View ArticleGoogle Scholar
- Serra M, Piña B, Abad JL, Camps F, Fabriàs G: A multifunctional desaturase involved in the biosynthesis of the processionary moth sex pheromone. Proc Natl Acad Sci USA. 2007, 104: 16444-16449. 10.1073/pnas.0705385104.PubMed CentralView ArticlePubMedGoogle Scholar
- Matoušková P, Pichová I, Svatoš A: Functional characterization of a desaturase from the tobacco hornworm moth (Manduca sexta) with bifunctional Z11 and 10,12-desaturase activity. Insect Biochem Molec Biol. 2007, 37: 601-610. 10.1016/j.ibmb.2007.03.004.View ArticleGoogle Scholar
- Park HY, Kim MS, Paek A, Jeong SE, Knipple DC: An abundant Acyl-CoA (Δ9) desaturase transcript in pheromone glands of the cabbage moth Mamestra brassicae, encodes a catalytically inactive protein. Insect Biochem Molec Biol. 2008, 38: 581-595. 10.1016/j.ibmb.2008.01.005.View ArticleGoogle Scholar
- Tóth M, Szöcs G, van Nieukerken EJ, Philipp P, Schmidt F, Francke W: Novel type of sex pheromone structure identified from Stigmella malella (Stainton) (Lepidoptera: Nepticulidae). J Chem Ecol. 1994, 21: 13-23. 10.1007/BF02033659.View ArticleGoogle Scholar
- Zhu J, Kozlov MV, Philipp P, Francke W, Löfstedt C: Identification of a novel moth sex-pheromone in Eriocrania cicatricella (Zett.) (Lepidoptera: Eriocraniidae) and its phylogenetic implications. J Chem Ecol. 1995, 21: 29-43. 10.1007/BF02033660.View ArticlePubMedGoogle Scholar
- Koslov MV, Zhu J, Philipp P, Francke W, Zvereva EL, Hansson BS, Löfstedt C: Pheromone specificity in Eriocrania semipurpurella (Stephens) and E. sangii (Wood) (Lepidoptera: Eriocraniidae) based on chirality of semiochemicals. J Chem Ecol. 1996, 22: 431-454. 10.1007/BF02033647.View ArticleGoogle Scholar
- Phelan PL: Evolution of mate-signalling in moths: phylogenetic considerations and predictions from the asymmetric tracking hypothesis. The evolution of mating systems in Insects and Arachnids. 1997, 240-256.View ArticleGoogle Scholar
- Löfstedt C, Zhu J, Kozlov MV, Buda V, Jirle EV, Hellqvist S, Löfqvist J, Plass E, Franke S, Francke W: Identification of the sex pheromone of the currant shoot borer Lampronia capitella. J Chem Ecol. 2004, 30: 643-657. 10.1023/B:JOEC.0000018635.40128.2e.View ArticlePubMedGoogle Scholar
- Löfstedt C, Kozlov M: A phylogenetic analysis of pheromone communication in primitive moths. Insect pheromone research. 1997, 473-489.View ArticleGoogle Scholar
- Patel O, Fernley R, Macreadie I: Saccharomyces cerevisiae expression vectors with thrombin-cleavable N- and C-terminal 6x(His) tags. Biotechnol Lett. 2003, 25: 331-334. 10.1023/A:1022384828795.View ArticlePubMedGoogle Scholar
- Schneiter R, Tatzer V, Gogg G, Leitner E, Kohlwein SD: Elo1-dependent carboxy-terminal elongation of C14:1Delta(9) to C16:1Delta(11) fatty acids in Saccharomyces cerevisiae. J Bacteriol. 2000, 182: 3655-3660. 10.1128/JB.182.13.3655-3660.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Reaney JTM, Liu YD, Taylor WG: Gas chromatographic analysis of Diels-Alder adducts of geometrical and positional isomers of conjugated linoleic acid. JAOCS. 2001, 78: 1083-1086. 10.1007/s11746-001-0393-x.Google Scholar
- Russel JN: Mechanisms of thermal adaptation in bacteria: blueprints for survival. Trends Biochem Sci. 1984, 9: 108-112. 10.1016/0968-0004(84)90106-3.View ArticleGoogle Scholar
- Maresca B, Cossins AR: Fatty feedback and fluidity. Nature. 1993, 365: 606-607. 10.1038/365606a0.View ArticlePubMedGoogle Scholar
- Los D, Horvath I, Vigh L, Murata N: The temperature-dependent expression of the desaturase gene desA in Synechocystis PCC6803. FEBS lett. 1993, 318: 57-60. 10.1016/0014-5793(93)81327-V.View ArticlePubMedGoogle Scholar
- Hurles M: Gene duplication: The genomic trade in spare parts. PLoS Biol. 2004, 2: 900-904. 10.1371/journal.pbio.0020206.View ArticleGoogle Scholar
- Ohno S: Evolution by Gene Duplication. 1970, Springer-Verlag, Berlin-Heidelberg-New YorkView ArticleGoogle Scholar
- Hughes AL: The evolution of functionally novel proteins after gene duplication. Proc R Soc Lond B Biol Sci. 1994, 256: 119-124. 10.1098/rspb.1994.0058.View ArticleGoogle Scholar
- Lynch M, O'Hely M, Walsh B, Force A: The probability of preservation of a newly arisen gene duplicate. Genetics. 2001, 159: 1789-1804.PubMed CentralPubMedGoogle Scholar
- Gaunt MW, Miles MA: An insect molecular clock dates the origin of the insects and accords with palaeontological and biogeographic landmarks. Mol Biol Evol. 2002, 19: 748-761.View ArticlePubMedGoogle Scholar
- Dallerac R, Labeur C, Jallon JM, Knipple DC, Roelofs WL, Wicker-Thomas C: A delta 9 desaturase gene with a different substrate specificity is responsible for the cuticular diene hydrocarbon polymorphism in Drosophila melanogaster. Proc Natl Acad Sci USA. 2000, 97: 9449-9454. 10.1073/pnas.150243997.PubMed CentralView ArticlePubMedGoogle Scholar
- Nei M, Rooney AP: Concerted birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.PubMed CentralView ArticlePubMedGoogle Scholar
- Cardé RT, Haynes KF: Structure of the pheromone communication channels in moths. Advances in insect chemical ecology. Edited by: Cardé RT, Millar JG. 2004, Cambridge University Press, 283-332.View ArticleGoogle Scholar
- Löfstedt C: Moth pheromone genetics and evolution. Phil Trans R Soc Lond B. 1993, 340: 167-177. 10.1098/rstb.1993.0055.View ArticleGoogle Scholar
- Buser H-R, Arn H, Guerin P, Rauscher S: Determination of double bond position in mono-unsaturated acetates by mass spectrometry of dimethyl disulfide adducts. Anal Chem. 1983, 55: 818-822. 10.1021/ac00257a003.View ArticleGoogle Scholar
- Marques FA, Millar JG, McElfresh S: Efficient method to locate double bond positions in conjugated trienes. J Chromatogr A. 2004, 1048: 59-65.View ArticlePubMedGoogle Scholar
- Berglund P, Holmquist M, Hedenström E, Hult K, Högberg H-E: 2-Methylalkanoic acids resolved by esterification catalysed by lipase from Candida rugosa: Alcohol chain length and enantioselectivity. Tetrahedron: Assymetry. 1993, 4: 1869-1878. 10.1016/S0957-4166(00)80427-9.View ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Higgins DJ, Thompson JD, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Printing and Shading of Multiple-Alignment files, BOXSHADE 3.21. [http://www.ch.embnet.org/index.html]
- Prediction of protein sorting signals and localization sites in amino acid sequences (PSORT). [http://psort.ims.u-tokyo.ac.jp/form2.html]
- National Centre for Biotechnology Information (NCBI). [http://www.ncbi.nlm.nih.gov]
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Briefings in Bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.