Email updates

Keep up to date with the latest news and content from BMC Biochemistry and BioMed Central.

Open Access Highly Accessed Research article

Identification of avian wax synthases

Eva-Maria Biester1*, Janine Hellenbrand1, Jens Gruber1, Mats Hamberg2 and Margrit Frentzen1

Author Affiliations

1 Institute of Biology I, RWTH Aachen University, (Worringer Weg 1), Aachen, (52074), Germany

2 Department of Medical Biochemistry and Biophysics, Karolinska Institute, (Scheeles Vag 2), Stockholm, (17176), Sweden

For all author emails, please log on.

BMC Biochemistry 2012, 13:4  doi:10.1186/1471-2091-13-4

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2091/13/4


Received:6 December 2011
Accepted:4 February 2012
Published:4 February 2012

© 2012 Biester et al; licensee 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.

Abstract

Background

Bird species show a high degree of variation in the composition of their preen gland waxes. For instance, galliform birds like chicken contain fatty acid esters of 2,3-alkanediols, while Anseriformes like goose or Strigiformes like barn owl contain wax monoesters in their preen gland secretions. The final biosynthetic step is catalyzed by wax synthases (WS) which have been identified in pro- and eukaryotic organisms.

Results

Sequence similarities enabled us to identify six cDNAs encoding putative wax synthesizing proteins in chicken and two from barn owl and goose. Expression studies in yeast under in vivo and in vitro conditions showed that three proteins from chicken performed WS activity while a sequence from chicken, goose and barn owl encoded a bifunctional enzyme catalyzing both wax ester and triacylglycerol synthesis. Mono- and bifunctional WS were found to differ in their substrate specificities especially with regard to branched-chain alcohols and acyl-CoA thioesters. According to the expression patterns of their transcripts and the properties of the enzymes, avian WS proteins might not be confined to preen glands.

Conclusions

We provide direct evidence that avian preen glands possess both monofunctional and bifunctional WS proteins which have different expression patterns and WS activities with different substrate specificities.

Background

Birds preen their feathers with a secretion produced by the uropygial gland, a holocrine bilobular gland located above their tail. The secretion consists of lipids, proteins and salts [1] and varies, for example, among species, age, season and sex [2-7]. These secretions confer different functions regarding sexual attraction, lubrication, waterproofing, antipathogenic effects and plumage maintenance [8-11]. Preen gland waxes show a high diversity of components; some species contain monoacyl esters, others diacyl esters or triacylglycerols (TAG). The distribution of fatty acids and alcohol residues is often unique, especially branched-chain, extremely long-chain or substituted fatty acids can be found here [12]. For instance, the preen gland secretions from chicken (Gallus gallus) consist of 50% wax diesters and 30% TAG [13]. Wax diesters contain erythro- and threo-alkane-2,3-diols with chain-lengths of 21 to 23 carbon atoms and saturated fatty acids of 12 to 20 carbon atoms [13,14]. Diester waxes are detected in other galliform birds as quail or pheasant as well as perching birds, pigeons, cranes or woodpeckers [12,15-17], although wax monoesters are the most frequently found components of preen gland secretions. Anseriformes like geese (Anser domesticus) contain wax monoesters in their preen gland secretion, in which di-, tri- or tetramethylated acyl groups are esterified with saturated straight-chain monoalcohols [11]. In goose 96% of the alcohol component is due to octadecanol, while the fatty acid residues consist of 76% 2,4,6,8-tetramethyldecanoic acid and 11% 2,4,6,8-tetramethylundecanoic acid [18]. In contrast to the wax monoesters of geese, those found in barn owl (Tyto alba) are rich in methyl-branched fatty alcohol and fatty acyl residues. About 60% of the components are monomethyl-branched, mainly 3-methyl-branched C13 or C17 acids and 2-, 3- or 4-methyl-branched fatty alcohols with 12 to 18 carbon atoms [19].

Production of wax esters has already been observed in preen gland membranes of chicken [20] and goose [21]in vitro, so it could be assumed, that genes essential for wax ester biosynthesis are expressed in preen gland tissue. The respective genes have not been identified in birds yet, but wax ester synthase sequences (acyl-CoA:alcohol acyltransferases, AWAT, WS) have already been described in other organisms including mammals [22,23], plants [24-27], bacteria [28-30] and protozoa [31]. Mammalian enzymes with wax synthase activity have been found within members of both DGAT1 and DGAT2 type acyltransferase families [22,32]. DGAT (acyl-CoA:diacylglycerol acyltransferases) catalyze the final step in storage lipid biosynthesis of TAG, but the human DGAT1 is capable of synthesizing wax monoesters, diesters and retinylesters as well [32,33]. Human wax synthases AWAT1 and AWAT2 belong to the DGAT2 type family. Like DGAT1, AWAT2 is a multifunctional acyltransferase which shows in vitro acyl-CoA:monoacylglycerol acyltransferase (MOGAT), DGAT, WS and retinylester synthase activities [33]. Bacterial wax synthases are at least bifunctional enzymes conferring WS activity next to DGAT and low MOGAT activity [28,29,34].

Wax esters are excellent lubricants because of their high stability under high temperature and pressure and high resistance to hydrolysis [35]. Unlike saturated long-chain monoesters, mono-unsaturated monoester or diester waxes combine good lubricity with good thermal and oxidative stability, high viscosity indices [36] and stability against lipases [37]. To achieve the renewable production of wax esters, currently attempts are made to identify new enzymes catalyzing respective esterification reactions [22,25,31,38]. Production of wax esters in oil crops [35,39] or microorganisms [38,40] might in future be able to surrogate fossil materials in technical industry.

Our studies identified WS genes of chicken, goose and barn owl, members of different bird families with distinct preen wax compositions. As the chicken genome was fully sequenced and assembled in March 2004 by the National Human Genome Research Institute [41,42], it served as a starting point for the identification of avian WS genes. Several sequences were successfully cloned and functionally analyzed in yeast cells.

Results

Identification of putative wax synthases from avian organisms

Sequence similarity based searches conducted with human AWAT1, AWAT2 and DGAT1 sequences as queries against the annotated Gallus gallus proteome resulted in five full-length sequences. Using mRNA isolated from preen glands as starting material, we succeeded in cloning the respective cDNAs of GgWS2 [NCBI: JQ031644], GgWS4 [NCBI: XP_419207.1], GgWS5 [NCBI: NP_0010261921.1] and GgDGAT1 [NCBI: JQ031642], while cDNAs of GgWS1 were synthesized because we failed to amplify the respective full-length sequence. Based on the WS sequences from chicken, we were able to clone WS4 and WS5 from preen glands of barn owl and goose.

The cloned 1056 bp GgDGAT1 sequence showed almost 90% identity with the annotated mRNA in databases [NCBI: XM_422267.2, ENSEMBL: ENSGALT00000006691], its central region, however, differed from both annotations. On protein level, GgDGAT1 displayed 98% identity to a recently identified sequence from turkey [NCBI: XP_003208594] but differed from both annotations of chicken (Additional File 1).

Additional file 1. Amino acid alignment of GgDGAT with database sequences. Amino acid alignment of GgDGAT1 [NCBI: JQ031642] and respective sequences from NCBI [XP_422267.2] and ENSEMBL [ENSGALT00000006691] database. The grey background highlights the different amino acids of the cloned protein and the predicted database sequences in the N-terminal region and the central part.

Format: PDF Size: 130KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The 1089 bp open reading frame (ORF) of GgWS2 differed from the annotations [NCBI: XM.426251.2, ENSEMBL: ENSGALT00000006967], especially in the 3'-terminal 60 nucleotides, while an alternation in the 5'-region causes a conservative substitution of V66I only (Additional File 2).

Additional file 2. Protein alignment of GgWS2 with database sequences. Protein alignments of GgWS2 [NCBI: JQ031644] and respective protein sequences on NCBI [XP_426251.2] and ENSEMBL [ENSGALT00000006967] database. The grey background highlights the different amino acids between the cloned sequence and the predicted database sequences at the C-terminus and the I66V substitution.

Format: PDF Size: 82KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Unlike GgWS2, the GgWS4 ORF of 951 bp resembled that deposited at NCBI database [NCBI: XM_419207.2] except for two conservative nucleotide substitutions (T300C and T312C), which did not alter the amino acid sequence [NCBI: XP_419207.1]. Comparison with the respective sequences from goose [NCBI: JQ031643] and barn owl [NCBI: JQ031645] showed about 90% identity at the cDNA level and 94% sequence identity at the protein level (Additional File 3).

Additional file 3. Protein alignment of WS4 homologs. Protein sequence alignment of WS4 homologs from chicken [Gg, NCBI: XP_419207], goose [Ad, NCBI: JQ031643] and barn owl [Ta, NCBI: JQ031645] WS4. The grey background highlights the differences in the amino acid sequences of different avian WS4 proteins.

Format: PDF Size: 74KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

WS5 homologs from chicken [NCBI: NM_001031021], goose [NCBI: JQ031647] and barn owl [NCBI: JQ031646] were identical at the nucleotide level. Contrasting the NCBI sequence, we found two nucleotide substitutions, T299C and T317G, which led to the exchange of V100A and V124G on protein level [NCBI: NP_001026192.1].

In silico analysis of putative avian wax synthases

Comparison of the avian WS protein sequences with respective sequences from different organisms gave the results illustrated in Figure 1. It reveals that GgWS1 and GgWS2 possess the highest similarity to human acyltransferases of the DGAT2-type (up to 60%) while GgDGAT1 shows the highest identity to DGAT1 family members, namely 62% to human acyl-CoA:cholesterol acyltransferase (ACAT) HsACAT1, 43% to HsACAT2 and 16% to HsDGAT1. Avian WS4 and WS5 proteins share less than 15% sequence identity to both DGAT1 and DGAT2 family members and build an own branch of enzymes (Figure 1). WS4 and WS5 proteins comprise 55% identity to each other and WS5 proteins are most similar to transmembrane protein 68, a protein of unknown function found in various organisms (Figure 1).

thumbnailFigure 1. Phylogenetic analysis of putative avian wax synthases in comparison to different acyltransferases from human, plants, yeast, Tetrahymena and bacteria. The phylogram shows the relation of the analyzed avian proteins with acyltransferases from different organisms. The tree is based on the alignment of the following sequences (NCBI accession numbers are given): HsAWAT1 [Homo sapiens, NP_001013597.1], HsAWAT2 [H. sapiens, NP_001002254.1], HsMOGAT3 [H. sapiens, NP_835470], HsMOGAT2 [H. sapiens, NP_079374], HsMOGAT1 [H. sapiens, NP_477513.2], HsACAT1 [H. sapiens, NP_003092.4], HsACAT2 [H. sapiens, NP_003569.1], HsDGAT2 [H. sapiens, NP_115953.2], HsDGAT1 [H. sapiens, NP_036211.2], SchWS [Simmondsia chinensis, AF149919_1], EgWS [Euglena gracilis, ADI60058.1], PxhWS [Petunia x hybrida, AAZ08051.1], AcWS/DGAT [Acinetobacter sp. ADP1, YP_045555.1], GgWS1 [Gallus gallus, XP_424082.2], GgWS2 [G. gallus, JQ031643], GgWS4 [G. gallus, XP_419207.1], GgWS5 [G. gallus, NP_001026192.1], AdWS5 [Anser domesticus, JQ031647], TaWS5 [Tyto alba, JQ031646], AdWS4 [A. domesticus, JQ031643], TaWS4 [T. alba, JQ031645], HsTMEM68 [H. sapiens, Q96MH6.2], GgDGAT1 [G. gallus, JQ031642], Tetrahymena WS [Tetrahymena thermophila, XP_001027910, XP_001026090, XP_001008104, XP_001019739], MhWS1 [Marinobacter hydrocarbonoclasticus, ABO21021.1], Arabidopsis WSD [Arabidopsis thaliana, NP_568547.1, NP_177356.1, NP_850307.1, NP_200151.2] MrDGAT2B [Umbelopsis ramanniana, AAK84180.1], MrDGAT2A [U. ramanniana,AAK84179.1], ScDGA1p [Saccharomyces cerevisiae, NP_014888.1], ScARE1p [S. cerevisiae, CAA42296.1], ScARE2p [S. cerevisiae, CAA96298.1], Arabidopsis WS [A. thaliana, NP_200345.1, XP_002866091.1, NP_200349.1, NP_200346.1]. The scale corresponds to amino acid substitutions per site in the alignment of 41 sequences with a total of 188 positions. Numbers at the branches are bootstrap values indicating the probability of this relationship in %. Values above 95 can be regarded as correct. The dendrogram was created with ClustalX2 and MEGA5 software.

The relation to the different mammalian acyltransferase families was reflected in further characteristics of the avian proteins like the molecular mass, the transmembrane structure and acyltransferase motifs. Conserved domain search revealed that all DGAT1 family members possessed an MBOAT (membrane bound O-acyltransferase) superfamily motif (Pfam cl00738) and the FYxDWWN motif (Figure 2a), which is a predicted acyl-CoA binding site identified in mammalian DGAT1 family members [43,44]. Furthermore, GgDGAT1 and the human ACAT proteins contained the (H/Y)SF motif [43] which might play a role in sterol binding. Unlike DGAT1 family members, DGAT2 family members like GgWS1 and GgWS2 were found to contain an acyltransferase superfamily motif (Pfam cl00357), that was found in WS4 and WS5 as well. In addition, the HPHG motif which is typical for DGAT2 family members and likely comprises a part of the active site [45] was also found in GgWS1 and GgWS2, while in WS4 and WS5 homologs the motif was modified to YYHG and FYHG, respectively (Figure 2b). GgWS1 and GgWS2 have masses of about 40 kDa with one predicted N-terminal transmembrane domain (TMD) and their C-terminal parts stretched into the cytosol like other DGAT2 related proteins. On the other hand, GgDGAT1 is a 60 kDa protein with several predicted TMDs typical of DGAT1 related proteins. With regard to the transmembrane structure, the avian WS4 homologs resemble DGAT2 family members unlike WS5 homologs, which differed insofar as the TMD prediction indicated that the C-terminal part of the protein was not located in the cytosol.

thumbnailFigure 2. Alignment of typical motifs of DGAT1 (a) and DGAT2 (b) families. Figure 2a shows the alignment of the FYxDWWN motif which is a potential acyl-CoA binding motif in proteins of the DGAT1 family [44]. Figure 2b represents the partially modified HPHG motif representing a potential part of the active site in DGAT2 family proteins [45]. NCBI accession numbers of reference proteins: HsDGAT1 [Homo sapiens, NP_036211.2], HsACAT1 [H. sapiens, NP_003092.4], HsACAT2 [H. sapiens, NP_003569.1], HsDGAT2 [H. sapiens, NP_115953.2], HsAWAT1 [H. sapiens, NP_001013597.1], HsAWAT2 [H. sapiens, NP_001002254.1].

Expression profiles of wax synthase sequences in chicken

Expression profiles of four chicken sequences were analyzed with different chicken tissues to reveal whether the sequences were preferentially expressed in preen glands. As the use of oligo(dT) primers did not lead to the amplification of any wax synthase sequence, gene-specific primers were used to synthesize the respective cDNAs from 1 μg of RNA. The comparison of products from the different tissues indicated that GgWS1 and GgDGAT1 were almost exclusively expressed in the preen gland, while GgWS2 and GgWS4 were expressed in the four analyzed tissues (Figure 3).

thumbnailFigure 3. Expression profiles of avian WS sequences in different tissues of chicken. RNA was converted to cDNA with gene specific primers and used as template for partial transcript amplification. Reactions were performed with or without (+ and - marks) reverse transcriptase (C, positive control with plasmid DNA).

Functional expression in yeast

To analyze the identities of the avian proteins, the respective cDNAs were expressed in a yeast mutant strain lacking TAG synthesis. Transgenic yeast cultures were supplemented with fatty alcohols from 10 to 18 carbon atoms to enable wax ester production. Lipid analyses of the transgenic yeast cells suggest that most of the avian proteins were functionally expressed and caused accumulation of storage lipids in significant but different levels and patterns (Figure 4a). Expression of GgWS1 resulted in the highest levels of wax esters which comprised 2 μmol wax ester/g fresh weight, while TAG was formed in very low levels only (Additional file 4). GgDGAT1 expressing yeast cells produced almost 1 μmol wax esters/g but no TAG (Additional file 5), whereas all WS4 homologs catalyzed the synthesis of higher levels of TAG (250 to 500 nmol/g) than of wax esters (100 to 250 nmol/g). The produced wax esters contained mainly dodecanol (12:0-OH) and tetradecanol (14:0-OH) esterified with palmitoleic (16:1) and oleic (18:1) acid irrespective of the expressed sequence. In contrast, TAG produced by WS4 homologs consisted of almost equal amounts of saturated and unsaturated fatty acyl residues. GgWS2 expressing cells contained storage lipids in very low levels similar to those of the control cells (Additional file 4). Additional feeding of GgWS2 expressing yeast cells with myristic acid (14:0) gave a 5-fold increase in wax production while the wax level of control strains was not affected. Such stimulation was also observed in GgWS4 expressing cells (Additional File 6).

thumbnailFigure 4. Production of lipids in transgenic yeast cultures. Cultures expressing one of the respective WS sequences were supplemented with a mixture of 125 μM 10:0-OH, 12:0-OH, 14:0-OH, 16:0-OH and 18:0-OH and induced for 48 hours. Lipids of the harvested cells were extracted and analyzed. Mean values and standard deviations of two independent experiments are given (a). TLC analysis of lipid extracts from yeast cells expressing GgFAR1 or co-expressing GgFAR1 and GgWS4 from cultures supplemented with 14:0-fatty acid (b). (FA: fatty acids, FOH: fatty alcohols, TAG: triacylglycerols, WE: wax esters)

Additional file 4. GC analyses of transmethylated wax esters and triacylglycerols. Lipids were extracted from transgenic yeast cells expressing the empty vector (control) or one of the avian proteins under standard conditions. WE and TAG were reextracted from TLC plates, transmethylated and analyzed by GC. (1) 10:0-OH, (2) 12:0-OH, (3) 14:0-OH, (4) 16:1-FAME, (5)16:0-FAME, (6) 16:0-OH, (7) 18:1-FAME, (8) 18:1-FAME, (9) 22:1-ME (internal standard)

Format: PDF Size: 195KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 5. GC analyses of intact wax esters and transmethylated triacylglycerols from control yeast strains and yeast cultures expressing GgDGAT1. The yeast cultures expressing the empty vector (control) or GgDGAT1 were cultivated under standard conditions, lipids were extracted and separated by TLC. WE were extracted from TLC and analyzed as intact WE. TAG were extracted, transmethylated and analyzed as methylester-derivatives. (1) 22:1-ME (internal standard in WE analysis, 30 nmol), (2) 26:1-WE, (3) 26:0-WE, (4) 28:1-WE, (5) 28:0-WE, (6) 30:1-WE, (7) 30:0-WE, (8) 32:1-WE, (9) 32:0-WE, (10) 16:1-ME, (11) 16:0-ME, (12) 18:1-ME, (13) 18:0-ME, (14) 22:0-ME (internal standard in TAG analysis, 30 nmol)

Format: PDF Size: 122KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 6. Wax ester production of transgenic yeast cells expressing GgWS2 or GgWS4 under different conditions. Yeast cultures expressing the empty vector (control), GgWS2 or GgWS4 were induced for 48 hours in SD-medium containing 125 μM 10:0-, 12:0-, 14:0-, 16:0- and 18:0-alcohol (A), 500 μM 14:0 alcohol and fatty acid (B) or 500 μM 14:0 and 16:0 alcohol and fatty acid (C). The lipids were extracted and analyzed by GC, the total WE-amounts per gram fresh weight are given.

Format: PDF Size: 204KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

To reconstitute the avian wax biosynthesis in yeast cells, WS were co-expressed with the recently identified fatty acyl-CoA reductase GgFAR1 [46]. This enzyme has a specificity for 16:0-acyl-chains, but is also able to produce 14:0-OH in yeast cultures supplemented with 14:0-FA [46]. GgWS4 was chosen for co-expression because it performed the highest WS activities of all analyzed WS enzymes with long chain acyl-acceptors. Supplementation with 14:0-FA of cultures expressing both GgWS4 and GgFAR1 resulted in the production of fatty alcohols, 350 nmol TAG/g and 550 nmol wax esters per g fresh weight, while expression of GgFAR1 alone led only to the production of fatty alcohols (Figure 4b).

Properties of avian proteins

To analyze the substrate specificities of the different avian proteins, in vitro assays were performed with membranes of yeast cells expressing the respective sequences as enzyme source. WS activity was detectable in membranes, but not in soluble fractions, which is in line with the predicted transmembrane domains. Under standard assay conditions the WS activities were constant for at least 2 μg protein in an incubation time of 20 minutes at 35°C (data not shown). This holds true with regard to yeast membranes harboring GgWS1, GgWS2 and the WS4 proteins. Unlike GgWS5, which did not show acyltransferase activities under any conditions and with any substrates tested, GgDGAT1 was catalytically active in yeast cells, but not in isolated yeast membranes (Additional File 7). Addition of detergents, divalent cations or bovine serum albumin did not improve but partially inhibited incorporation rates of the labeled 16:0-acyl-groups into lipophilic reaction products (data not shown).

Additional file 7. TLC analysis of lipophilic reaction products from WS assays with yeast membranes. Assays were performed with 16:0-CoA and 10:0-OH under standard conditions using membranes of yeast cells expressing one of the respective sequences. Reaction products were extracted from the assays, separated by TLC and visualized using the FLA-3000 imaging system. The analysis is representative of several repetitions.

Format: PDF Size: 122KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

GgWS1, GgWS2 and the avian WS4 proteins showing WS activity were further analyzed concerning their reaction products and their acyl-donor and -acceptor specificities. As given in Figure 5, GgWS1 and GgWS2 formed wax esters as main reaction products, while WS4 homologs possessed both WS and DGAT activities. These results which were in line with yeast expression experiments (Figure 4) clearly demonstrate the mono- or bifunctionality of the respective avian WS (Additional File 7).

thumbnailFigure 5. Acyl-CoA specificities of avian enzymes. Relative WS (a and c) and DGAT (b) activities of isolated yeast membranes from yeast cells expressing an avian WS sequence. Assays were conducted with 10:0-OH and the given straight-chain (a and b) or branched-chain (c) acyl-CoA thioesters under standard assay conditions (a and b) or increased protein amount (10 μg) and incubation time (2 h) (c). The relative activities of each protein are given, 100% corresponds to the activities with 10:0-OH and 16:0-CoA under identical conditions. Values are mean values from at least two independent assay series. (GgWS1: 116 pmol*min-1*mg-1; GgWS2: 105 pmol*min-1*mg-1; GgWS4: 149 pmol*min-1*mg-1; AdWS4: 126 pmol*min-1*mg-1; TaWS4: 149 pmol*min-1*mg-1).

With regard to the acyl-CoA specificities, avian enzymes were most active with saturated acyl-CoA thioesters of 14 to 18 carbon atoms, while the WS and DGAT activities were rather low with shorter or unsaturated thioesters. GgWS1 and GgWS2 were almost exclusively active with 16:0-CoA while the avian WS4 homologs produced the highest wax amounts with 16:0- or 18:0-CoA and the highest TAG levels with 14:0-CoA (Figure 5a and 5b). Hence bifunctional WS, especially those of goose and barn owl, displayed different acyl-CoA specificities regarding their WS and DGAT activities.

In addition to straight-chain acyl-donors, we assayed 2-methyl-branched acyl-CoA thioesters of 14 to 18 carbon atoms as especially barn owls contain many methyl-branched fatty acids in their preen gland secretion [19]. To detect WS activity with these substrates we had to increase the protein amount and incubation time. Under these conditions GgWS1 and GgWS2 showed appreciable activities with branched-chain acyl-CoA in comparison to 16:0-CoA, while activity of WS4 proteins were very low in comparison to 16:0-CoA (Figure 5c). GgWS1 displayed the highest activity with 2-methyl-branched 16:0-CoA that comprised 70% of the activity determined with 16:0-CoA under the same conditions. GgWS2 performed lower WS activities than GgWS1 but the substrate specificities of both enzymes were very similar to each other (Figure 5c).

Analysis of acyl-acceptor specificities showed that all proteins performed the highest WS activities with decanol (10:0-OH) or undecanol (11:0-OH) in combination with 16:0-CoA, whereas no activities could be observed with alcohols shorter than 9 carbon atoms (Figure 6). While the enzymes showed very similar activity patterns with saturated alcohols of 9 to 16 carbon atoms, assays with unsaturated or branched-chain alcohols revealed differences in the substrate specificities of the enzymes. Unsaturated alcohols were no suitable substrates for GgWS1 or GgWS2 while the WS4 proteins used 16:1-OH and 18:1-OH with higher activities than saturated fatty alcohols of the same chain length. With regard to branched-chain substrates, we analyzed 3,7-dimethyl-octanol (3,7-diMe-8:0-OH) and isoprenols as substrates. Although octanol (8:0-OH) was not esterified by any enzyme, the 3,7-diMe-8:0-OH was esterified with high activities by all enzymes except GgWS4 and AdWS4. Isoprenols were converted to prenyl esters by AdWS4 and TaWS4 proteins, while chicken proteins showed low prenyl ester synthase activities only. Among the WS4 enzymes, geranylgeraniol was the preferred alcohol of TaWS4 and AdWS4 in contrast to the respective chicken protein (Figure 6). As diesters are the naturally occurring waxes in chicken preen glands [13,14], we analyzed 1,2- and 1,12-isomers of dodecanediol as substrates. None of the avian enzymes was capable of esterifying these diols, although activity could be detected in assays with membranes isolated from chicken preen glands (Figure 7).

thumbnailFigure 6. Acyl-acceptor specificities of avian proteins. Relative WS activities of membrane fractions of yeast cells harboring an avian enzyme with 16:0-CoA and the given acyl-acceptors under otherwise standard conditions. 100% corresponds to the activity with 16:0-CoA and 10:0-OH (GgWS1: 116 pmol*min-1*mg-1; GgWS2: 105 pmol*min-1*mg-1; GgWS4: 149 pmol*min-1*mg-1; AdWS4: 126 pmol*min-1*mg-1; TaWS4: 149 pmol*min-1*mg-1), values are mean values from at least two independent assay series.

thumbnailFigure 7. Properties of WS activities of preen gland membranes of chicken, goose and barn owl. Relative WS activities of preen gland membranes with 12:0-OH and different acyl-CoA thioesters (a) or 16:0-CoA and different alcohols (b). 100% corresponds to the turnover of 12:0-OH with 16:0-CoA. Standard assays were used apart of higher protein amounts (50 μg) and incubation times (2 h).

In summary, we have identified both monofunctional and bifunctional WS enzymes from birds, which differ in their substrate specificities, especially with regard to branched-chain substrates.

In vitro analysis of WS activities in preen gland membranes

To analyze the in vitro WS activities of preen glands, assays were performed with isolated membrane preparations. Strong WS activity could be detected in membranes of chicken preen glands, but comparably low activities were obtained in skin or liver preparations, suggesting that WS activity is restricted to preen gland tissue.

Figure 7a displays the acyl-donor specificities of avian WS using preen gland membranes of chicken, goose and barn owl as enzyme source. The highest WS activities were obtained with 14:0-CoA and almost no activity with 18:1-CoA. In contrast to chicken and goose, barn owl performed high activities with 2-methyl-branched acyl-CoA thioesters as well (Figure 7a).

The comparison of different acyl-acceptors proved that membranes of preen glands of all tested species were able to catalyze the esterification of straight-chain alcohols, especially 12:0-OH, or branched-chain alcohols like geranylgeraniol or 3,7-diMe-8:0-OH (Figure 7b). In addition to WS, membranes of chicken preen glands exhibited strong DGAT activity and also high diester synthase activity with 1,2-dodecanediol.

Discussion and Conclusion

In this study we present the identification and characterization of the first avian wax synthase sequences. Although WS activity has already been demonstrated in vitro for cell free preparations from preen glands of chicken and turkey [47], neither a specific enzyme catalyzing this reaction nor the respective gene has been identified yet. Six of nine proteins identified by sequence homologies catalyzed wax ester syntheses in vivo or in vitro (GgDGAT1, GgWS1, GgWS2, GgWS4, AdWS4 and TaWS4) and WS4 homologs additionally showed DGAT activity.

With most avian proteins, in vivo and in vitro experiments gave consistent results. WS activity of GgDGAT1, however, was clearly detectable in yeast cells (Figure 4), but not in enzymatic assays with yeast membranes, in spite of enzymatic assays being more sensitive than feeding experiments. Human ACAT1, to which GgDGAT1 displays the highest sequence identity, was functionally expressed in yeast cells and showed in vitro ACAT activity with cholesterol and oleate [43], but GgDGAT1 did not show ACAT activity. Perhaps we have not tested GgDGAT1 under conditions suitable for this enzyme so far, or the chicken enzyme has a very low protein stability under in vitro conditions and is rapidly degraded by yeast proteases. This still remains to be determined.

The protein termed WS5 is highly conserved among vertebrates and identical in different bird species. It contains a lysophospholipid acyltransferase motif (cl00357 and cd07987), but the enzymatic activities have not been studied in any organism to date. Our experiments could not confirm any acyltransferase activity in vivo or in vitro.

All avian WS are found to be most active with saturated medium-chain alcohols (10:0-OH to 12:0-OH) and saturated long-chain acyl-CoA thioesters (14:0 to 18:0), but they differ in their reaction products and their activities with certain substrates. GgWS4, like the respective homologs from goose and barn owl, are bifunctional enzymes which catalyze both WE and TAG biosynthesis. They are more active with unsaturated than with saturated long-chain alcohols and effectively utilize branched-chain alcohols like isoprenols as acyl-acceptors, especially AdWS4 and TaWS4, but they show hardly any activities with branched-chain acyl-CoA. In contrast to the WS4 homologs, GgWS1 and GgWS2 are monofunctional enzymes which are inactive with saturated and unsaturated alcohols of more than 14 carbon atoms (Figure 6), but are active with methyl-branched alcohols and branched-chain acyl-CoA thioesters (Figure 5c). The ability of the chicken enzyme to utilize branched-chain substrates was surprising, as these components are typical for wax esters of barn owl and goose preen glands, but not for those of chicken. Perhaps WS1 and WS2 homologs from barn owl and goose might be even more active with branched-chain substrates than the chicken enzymes. These homologs likely caused the relatively high WE formation rates determined in assays with preen gland membranes of barn owl (Figure 7), while the relatively low activities of preen gland membranes of goose might be due to the wax esters containing tetramethylated but not monomethylated acyl-groups like barn owl glands [18,19].

The secretion of chicken preen glands is rich in 2,3-diester waxes [14,20] and preen gland membranes of chicken, unlike those of goose and barn owl, effectively catalyze the esterification of 1,2-dodecanediol (Figure 7). On the other hand, diester synthase activities with 1,2-dodecanediol were not detected with transgenic yeast membranes harboring one of the various avian enzymes even when WS assays with yeast membranes were run under conditions identical to those with preen gland membranes. These data suggest that diester waxes might be formed by a WS in chicken which is unrelated to known WS classes or by a related protein, which could not be detected in yeast membranes so far, like GgDGAT1. Anyhow, chicken preen gland membranes displayed high monoester synthase activities (Figure 7), which is in line with previous findings [47] and supports our results of the WS enzymes expressed in yeast (Figures 5 and 6).

Methods

Identification of putative wax synthases

BLASTp (protein-to-protein Basic Local Alignment Search Tool) [48] studies were undertaken to investigate the predicted Gallus gallus proteome for putative wax synthases using the NCBI server [49]. Human wax synthase sequences HsAWAT1 [NCBI: NP_001013597] and HsAWAT2 [NCBI: NP_001002254] were used as query sequences for DGAT2 family members and human HsDGAT1 [NCBI: NP_036211] was used as query for DGAT1 like proteins. The obtained results were drawn back on to the originating mRNA sequences, which were used for generating cloning PCR primers.

Vector construction

Preparation of mRNA, cDNA synthesis and cloning was performed as published previously [46]. Specific primer sequences were deduced from NCBI sequences (DGAT1-For, 5'-ATGGCAGGTGAAGACTGTGTAAG-3', DGAT1-Rev, 5'-TTACATCTGCACGTGACATGACCAC-3', WS2-For, 5'-CACCATGAAAACAATCATTGCAGCCTG-3', WS2-Rev, 5'-TCAACTCCTGGGCCATGTGG-3', WS4-For, 5'-CACCATGACCTACCTAAGCTACTTTGC-3', WS4-Rev, 5'-CTAATCACTTTTACAATGTTTATC-3', WS5-For, 5'-CACCATGATAGGTAGCAATGAATCC-3', WS5-Rev, 5'-TTAGTCTTCTTTCTGTCGTGTTTG-3'). GgWS1 was synthesized (GeneArt) and derived in Gateway® compatible entry vectors. All sequences were cloned into yeast expression vector pYES-DEST52 (Invitrogen) and for co-expression of wax synthases with the chicken fatty acyl-CoA reductase GgFAR1 [46] the respective sequences were cloned into pAG425GAL vectors (Addgene plasmid 14203, Susan Lindquist). Cotransformation of pYES-DEST52 and pAG425GAL enabled the selection on medium lacking uracil and leucine. The respective destination vectors were transformed into Saccharomyces cerevisiae BY4741 Δlro1 Δdga1 (MATa, his3Δ1, leu2Δ0, met15Δ0, ura3Δ0, lro1-Δ::kanMX4, dga1-Δ::natMX4) [50].

Analysis of WS expression profiles

Expression profiles were analyzed in tissues of chicken pectoral muscle, skin, liver and preen gland. 1 μg of total RNA was digested with DNase (Fermentas) to remove genomic DNA contamination, cDNA synthesis was conducted as described previously [46]. The following reverse primers were used for first strand synthesis: DGAT1-Rev, 5'-TTACATCTGCACGTGACATGACCAC-3', WS1-Rev, 5'-CTATATAAATTCGAGGTGACTGTCTTCTG -3', WS2-Rev, 5'-TCAACTCCTGGGCCATGTGG-3', WS4-Rev, 5'-CTAATCACTTTTACAATGTTTATC-3'. Negative controls without reverse transcriptase were performed and treated likewise. The following PCR was conducted using SupraTherm Taq polymerase (Genecraft). The amplified internal fragments were chosen to span introns, in case of residual chromosomal DNA. (GgDGAT1-exon6 for, 5'-TGCATGTTCTGTGCCACGGTT-3', GgDGAT1-exon12 rev, 5'-AAGTGAACCAAGCACCTGTGCA-3', 435 bp fragment length, GgWS1-exon4 for, 5'-TCTTCCACGCAAGAGGTATC-3', GgWS1-exon5 rev 5'-GTTGGCTTGTGCTTCTTCAG-3', GgWS2-exon4 for, 5'-TTGCCGTGCCTGAGGAGATG-3' 112 bp fragment length, GgWS2-exon5 rev, 5'-CTGGGAACAGCTCGGAGAAG-3', 190 bp fragment length, GgWS4-exon3 for, 5'-GCTCTAGATGCGAAAGTGCCCA-3', GgWS4-exon5 rev, 5'-GGAGGAGTGAATGGCAATCG-3'; 129 bp fragment length). Plasmid DNA was used as positive control. The PCR products were analyzed on 2% agarose gels in TAE buffer.

Phylogenetic analysis and structure prediction

Sequence analyses were carried out using ClustalX2 [51] and GeneDoc [52] software. Phylograms were computed with MEGA5 [53] and neighbor-joining method [54] with 1000 bootstrap replicates using p-distance method. All gaps were deleted for computation of evolutionary distances.

Molecular mass and isoelectric points were calculated using ProtParam [55] on the ExPASy Server [56]. Transmembrane helices of avian proteins and mammalian wax synthases were predicted using TMHMM software [57-60]. Predictions were compared to Kyte Doolittle plots [61] with window parameters of 19 which revealed similar results. Acyltransferase superfamily motifs and putative acyl-acceptor binding pockets were discovered by NCBI conserved domain search [62,63].

Functional analysis in yeast

The S. cerevisiae BY4741 Δlro1 Δdga1 strain deficient in TAG synthesis [50] was used for expression experiments. Transgenic yeast cells expressing an avian sequence in pYES-DEST52 vectors (Invitrogen) were cultivated in SD minimal medium containing 0.17% (w/v) yeast nitrogen base (MP Biomedicals), 0.068% (w/v) complete supplement medium without uracil and leucine, 0.5% NH4SO4 (w/v), 0.01% (w/v) leucine and 2% glucose for 24 h at 28°C. The cells of 50 ml cultures were induced with 2% galactose for 48 h and were supplemented with 125 μM decanol, dodecanol, tetradecanol, hexadecanol and octadecanol (Sigma Aldrich). Transgenic yeast cells harboring both pYES-DEST52-GgFAR1 and pAG425GAL-GgWS4 were cultivated in SD medium without uracil and leucine and supplemented with 250 μM myristic acid. Cells were harvested, washed and stored at -20°C.

Extraction of yeast cells was performed according to Bligh and Dyer [64]. The lipid extracts were separated by TLC on preparative TLC plates (Silica Gel 60 plates 0.5 mm thickness, Merck) in heptane/diethyl ether/acetic acid (90/30/1, v/v/v) and visualized under UV light after spraying with dichlorofluorescein (0.3% (w/v) dissolved in isopropanol) [65]. Myristoyl-dodecanoate (Sigma Aldrich) and TAG isolated from sunflower oil were used as standards.

GC analysis of WE and TAG

Bands co-chromatographing with the WE and TAG standards were transmethylated in 0.5 M sulfuric acid and 3% dimethoxypropane in methanol for 1 hour at 80°C together with 250 nmol docosanoic acid as internal standard. Fatty acid methyl esters (FAMEs) and fatty alcohols were extracted with heptane, concentrated and analyzed via gas chromatography (GC) with flame ionization detection (FID). For quantification of WE the total amount of fatty alcohols in the fractions was summarized, for quantification of TAG the sum of FAMEs was divided by 3.

GC-FID analysis was carried out using the HP6890 gas chromatograph equipped with an OPTIMA225 column (Macherey & Nagel) (25 m length, 0.25 mm diameter, 0.25 μm film thickness). 1 μl of the extract was analyzed in splitless injection with N2 as carrier gas (constant flow, 0.9 bar pressure, total column flow 1 ml/min) and inlet and detector temperatures of 260°C. A temperature program was carried out starting at 120°C, 8°C/min to 144°C, 4°C/min to 240°C. Peaks were identified by comparison of the respective retention times with those of standard substances of different fatty alcohols and FAMEs (Sigma Aldrich).

Synthesis of [1-14C]-labeled branched-chain fatty acids

[1-14C]2-Methyltetradecanoic, 2-methylhexadecanoic and 2-methyloctadecanoic acids were prepared by α-methylation of the corresponding [1-14C]-labeled fatty acids via the sequence carboxylic acid → acyl chloride → diazoketone → chloroketone → 2-methylcarboxylic acid essentially as described [66]. Purification by reversed-phase HPLC (solvent system, acetonitrile/water/acetic acid 85:15:0.01, v/v/v) afforded > 98% pure materials having a specific radioactivity of 0.622 GBq/mmol.

Preparation of yeast membranes and in vitro wax synthase assay

Membrane preparation and WS assays were performed as described previously [67]. Briefly, transgenic yeast cells were harvested, washed in Tris-H2SO4 (50 mM, pH 7.6), frozen and disrupted. Cell supernatants were combined and sonicated, cell debris was sedimented (2,500 × g, 15 min and 4°C) and the membranes were isolated from the supernatant by high speed centrifugation (1 h, 140,000 × g, 4°C). The sedimented membranes were resuspended in Tris-H2SO4 buffer and stored in aliquots at -80°C. The protein concentration was determined [68].

WS activity was measured with unlabeled acyl-acceptors and labeled acyl-CoA-thioesters as outlined before [63]. The following acyl-CoA thioesters were used: [1-14C]-myristoyl-CoA and [1-C14]-stearoyl-CoA (Biotrend), specific activity 2.03 Bq/pmol, [1-14C]-palmitoyl-CoA and [1-14C]-oleoyl-CoA (Perkin Elmer), specific activity 2.22 Bq/pmol and 2.03 Bq/pmol, [1-14C]-decanoyl-CoA (0.09 Bq/pmol) and [1-14C]-dodecanoyl-CoA (0.7 Bq/pmol) and 2-methyl-acyl-CoA (14:0, 16:0 and 18:0) (0.62 Bq/pmol) kindly provided by Prof. Sten Stymne and members of his laboratory, SLU Alnarp, Sweden. The reaction mixture of standard assays consisted of 10 mM BIS-Tris-propane buffer (pH 9), 13 μM [1-14C]-labeled acyl-CoA and 300 μM acyl-acceptor (Sigma Aldrich) using 2 to 4 μg protein of total yeast membrane fractions as enzyme source. Acyl-acceptors were dissolved in heptane and evaporated to dryness in the reaction tubes before addition of further assay components. Incubation was carried out at 35°C for 20 minutes. Lipids were extracted, applied to TLC silica gel plates (Merck) and chromatographed in heptane/diethyl ether/acetic acid (90/20/1 v/v/v). The bands were visualized with the FLA-3000 bioimager system (Fujifilm) and quantified in a liquid scintillation counter LS 6500 (Beckman Coulter). DGAT activity was measured as by-product of WS activity without addition of DAG, endogenous substrates were used for TAG synthesis.

Before calculation of relative activities, background activities of membranes from yeast control strains were subtracted.

WS assay with isolated preen gland membranes

Membranes from dissected preen glands were prepared as described before [21]. Briefly, the secretions were removed, the tissues were cut into small pieces, homogenized in NaK-buffer (0.1 M Na-K-phosphate buffer, pH 7.6, 0.5 mM DTT, 1 mM MgCl2) using an Ultra-Turrax three times for 15 seconds. Residual fragments were removed by centrifugation at 2,500 × g for 15 min and membranes were sedimented by ultracentrifugation at 140,000 × g for one hour. The membranes were resuspended in an appropriate volume of NaK-buffer. The protein content of the prepared membrane fractions was determined [68] and about 50 μg protein was used in the WS assay as outlined above, but incubation time was extended to 2 hours.

Abbreviations

ACAT: acyl-CoA:cholesterol acyltransferase, retinyl ester synthase; CoA: Coenzyme A; DGAT: acyl-CoA:diacylglycerol acyltransferase; MOGAT: acyl-CoA:monoacylglycerol acyltransferase; ORF: open reading frame; TAG: triacylglycerols; TMD: transmembrane domain; WE: wax esters; WS: wax synthase

Authors' contributions

EMB designed and carried out experiments, analyzed all data and drafted the journal. JH identified and cloned GgFAR1, participated in the design of experiments and data analyses and revised the manuscript. JG cloned GgDGAT1 and participated in lipid analyses. MH synthesized methyl-branched fatty acids and wrote the corresponding method section. MF coordinated the study, was involved in drafting the manuscript and revised the manuscript. All authors read and approved the manuscript.

Acknowledgements

This work was funded in part by the European Commission through the FP7 ICON project. We would like to thank Professor Sten Stymne (SLU Alnarp, Sweden) for coordination of the ICON project and him and his research team for doing and teaching us the acyl-CoA synthesis. Furthermore we would like to thank Professor Hermann Wagner (RWTH Aachen, Germany) and his team for the preparation of preen glands from barn owls and Mr. Peter Ritte (Wegberg-Rickelrath, Germany) for the supply with tissues from chicken and goose.

References

  1. Cater DB, Lawrie NR: Some Histochemical and Biochemical Observations on the Preen Gland.

    J Physiol-London 1950, 111:231-243. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Haahti E, Lagerspetz K, Nikkari T, Fales HM: Lipids of the Uropygial Gland of Birds.

    Comparative Biochemistry and Physiology 1964, 12:435-436. PubMed Abstract | Publisher Full Text OpenURL

  3. Jacob J, Balthazart J, Schoffeniels E: Sex Differences in the Chemical Composition of Uropygial Gland Waxes in Domestic Ducks.

    Biochemical Systematics and Ecology 1979, 7:149-153. Publisher Full Text OpenURL

  4. Piersma T, Dekker M, Damste JSS: An Avian Equivalent of Make-Up?

    Ecol Lett 1999, 2:201-203. Publisher Full Text OpenURL

  5. Whittaker DJ, Soini HA, Atwell JW, Hollars C, Novotny MV, Ketterson ED: Songbird Chemosignals: Volatile Compounds in Preen Gland Secretions Vary among Individuals, Sexes, and Populations.

    Behavioral Ecology 2010, 21:608-614. Publisher Full Text OpenURL

  6. Jacob J, Poltz J: Composition of Uropygial Gland Secretions of Birds of Prey.

    Lipids 1975, 10:1-8. PubMed Abstract | Publisher Full Text OpenURL

  7. Kolattukudy PE, Sawaya WN: Age Dependent Structural Changes in the Diol Esters of Uropygial Glands of Chicken.

    Lipids 1974, 9:290-292. PubMed Abstract | Publisher Full Text OpenURL

  8. Haribal M, Dhondt A, Rosane D, Rodriguez E: Chemistry of Preen Gland Secretions of Passerines: Different Pathways to Same Goal? Why?

    Chemoecology 2005, 15:251-260. Publisher Full Text OpenURL

  9. Salibian A, Montalti D: Physiological and Biochemical Aspects of the Avian Uropygial Gland.

    Brazilian Journal of Biology 2009, 69:437-446. OpenURL

  10. Hirao A, Aoyama M, Sugita S: The Role of Uropygial Gland on Sexual Behavior in Domestic Chicken Gallus gallus domesticus.

    Behavioural Processes 2009, 80:115-120. PubMed Abstract | Publisher Full Text OpenURL

  11. Morr M, Wray V, Fortkamp J, Schmid RD: (2r,4r,6r,8r)-2,4,6,8-Tetramethyldecan- und -Undecansäure aus dem Bürzeldrüsenwachs der Hausgans, Anser a. f. domesticus: Isolierung, Synthese einiger Derivate sowie der Rac-2,4,6,8-Tetramethyldecansäure.

    Liebigs Annalen der Chemie 1992, 1992:433-439. Publisher Full Text OpenURL

  12. Gamo M, Saito K: The Distribution of Diol Waxes in Preen Glands of Some Birds.

    Comparative Biochemistry and Physiology Part B: Comparative Biochemistry 1971, 39:151-157. Publisher Full Text OpenURL

  13. Wertz PW, Stover PM, Downing DT: A Survey of Polar and Nonpolar Lipids from Epidermis and Epidermal Appendages of the Chicken (Gallus domesticus).

    Comp Biochem Physiol B 1986, 84:203-206. PubMed Abstract | Publisher Full Text OpenURL

  14. Hansen IA, Tang BK, Edkins E: Erythro-Diols of Wax from the Uropygial Gland of the Turkey.

    J Lipid Res 1969, 10:267-270. PubMed Abstract | Publisher Full Text OpenURL

  15. Gamo M: The Distribution of Diol Waxes in Preen Glands of Some Birds. II.

    International Journal of Biochemistry 1971, 2:574-580. Publisher Full Text OpenURL

  16. Saito K, Gamo M: The Distribution of Diol Waxes in Preen Glands of Some Birds--III. The Occurrence of 1,2-Diols.

    Comparative Biochemistry and Physiology Part B: Comparative Biochemistry 1973, 45:603-606. Publisher Full Text OpenURL

  17. Saito K, Gamo M: The Occurrence of Diesters of 2,3-Dihydroxyoctadecane in the Preen Gland of Green Pheasant (Phasianus colchicus).

    J Biochem 1970, 67:841-849. PubMed Abstract | Publisher Full Text OpenURL

  18. Murray KE: Studies in Waxes XXI. Branched-Chain Acids of Preen Gland Wax of Goose.

    Aust J Chem 1962, 15:510-520. Publisher Full Text OpenURL

  19. Jacob J, Poltz J: Chemical Composition of Uropygial Gland Secretions of Owls.

    J Lipid Res 1974, 15:243-248. PubMed Abstract | Publisher Full Text OpenURL

  20. Tang BY, Hansen IA: Lipogenesis in Chicken Uropygial Glands.

    Eur J Biochem 1972, 31:372-377. PubMed Abstract | Publisher Full Text OpenURL

  21. Buckner JS, Kolattukudy PE: Lipid Biosynthesis in the Sebaceous Glands. Synthesis of Multibranched Fatty Acids from Methylmalonyl-Coenzyme A in Cell-Free Preparations from the Uropygial Gland of Goose.

    Biochemistry 1975, 14:1774-1782. PubMed Abstract | Publisher Full Text OpenURL

  22. Cheng JB, Russell DW: Mammalian Wax Biosynthesis II. Expression Cloning of Wax Synthase cDNAs Encoding a Member of the Acyltransferase Enzyme Family.

    J Biol Chem 2004, 279:37798-37807. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Turkish AR, Henneberry AL, Cromley D, Padamsee M, Oelkers P, Bazzi H, Christiano AM, Billheimer JT, Sturley SL: Identification of Two Novel Human Acyl-CoA Wax Alcohol Acyltransferases: Members of the Diacylglycerol Acyltransferase 2 (DGAT2) Gene Superfamily.

    J Biol Chem 2005, 280:14755-14764. PubMed Abstract | Publisher Full Text OpenURL

  24. Lassner MW, Lardizabal K, Metz JG: Producing Wax Esters in Transgenic Plants by Expression of Genes Derived from Jojoba. In Perspectives on New Crops and New Uses. Edited by Alexandria JJ. VA: ASHS Press; 1999:220-224. OpenURL

  25. Lardizabal KD, Metz JG, Sakamoto T, Hutton WC, Pollard MR, Lassner MW: Purification of a Jojoba Embryo Wax Synthase, Cloning of its cDNA, and Production of High Levels of Wax in Seeds of Transgenic Arabidopsis.

    Plant Physiol 2000, 122:645-656. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. King A, Nam JW, Han J, Hilliard J, Jaworski JG: Cuticular Wax Biosynthesis in Petunia Petals: Cloning and Characterization of an Alcohol-Acyltransferase that Synthesizes Wax-Esters.

    Planta 2007, 226:381-394. PubMed Abstract | Publisher Full Text OpenURL

  27. Li F, Wu X, Lam P, Bird D, Zheng H, Samuels L, Jetter R, Kunst L: Identification of the Wax Ester Synthase/Acyl-Coenzyme A: Diacylglycerol Acyltransferase WSD1 Required for Stem Wax Ester Biosynthesis in Arabidopsis.

    Plant Physiol 2008, 148:97-107. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Kalscheuer R, Steinbüchel A: A Novel Bifunctional Wax Ester Synthase/Acyl-CoA:Diacylglycerol Acyltransferase Mediates Wax Ester and Triacylglycerol Biosynthesis in Acinetobacter calcoaceticus ADP1.

    J Biol Chem 2003, 278:8075-8082. PubMed Abstract | Publisher Full Text OpenURL

  29. Stöveken T, Kalscheuer R, Malkus U, Reichelt R, Steinbüchel A: The Wax Ester Synthase/Acyl Coenzyme A:Diacylglycerol Acyltransferase from Acinetobacter Sp. Strain ADP1: Characterization of a Novel Type of Acyltransferase.

    J Bacteriol 2005, 187:1369-1376. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Holtzapple E, Schmidt-Dannert C: Biosynthesis of Isoprenoid Wax Ester in Marinobacter hydrocarbonoclasticus DSM 8798: Identification and Characterization of Isoprenoid Coenzyme A Synthetase and Wax Ester Synthases.

    J Bacteriol 2007, 189:3804-3812. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Teerawanichpan P, Qiu X: Fatty Acyl-CoA Reductase and Wax Synthase from Euglena gracilis in the Biosynthesis of Medium-Chain Wax Esters.

    Lipids 2010, 45:263-273. PubMed Abstract | Publisher Full Text OpenURL

  32. Yen CL, Monetti M, Burri BJ, Farese RV Jr: The Triacylglycerol Synthesis Enzyme DGAT1 also Catalyzes the Synthesis of Diacylglycerols, Waxes, and Retinyl Esters.

    J Lipid Res 2005, 46:1502-1511. PubMed Abstract | Publisher Full Text OpenURL

  33. Yen CL, Brown CHt, Monetti M, Farese RV Jr: A Human Skin Multifunctional O-Acyltransferase That Catalyzes the Synthesis of Acylglycerols, Waxes, and Retinyl Esters.

    J Lipid Res 2005, 46:2388-2397. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Kalscheuer R, Luftmann H, Steinbüchel A: Synthesis of Novel Lipids in Saccharomyces cerevisiae by Heterologous Expression of an Unspecific Bacterial Acyltransferase.

    Appl Environ Microbiol 2004, 70:7119-7125. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Carlsson AS, Lindberg Yilmaz J, Green AG, Stymne S, Hofvander P: Replacing Fossil Oil with Fresh Oil - with What and for What?

    European Journal of Lipid Science and Technology 2011, 113:812-831. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Kohashi H: Application of Fatty Acid Esters for Lubricating Oil. In World Conference on Oleochemicals into the 21st Century: Proceedings. Edited by Applewhite TH. Institiut Penyelidikan Minyak Kelapa Sawit Malaysia: American Oil Chemists' Society, ASEAN Oleochemical Manufacturers Group; 1991:243-255. OpenURL

  37. Gamo M, Saito K: Enzymatic Hydrolysis of Diol Waxes by Rat Liver Microsome.

    J Biochem 1970, 67:607-608. PubMed Abstract | Publisher Full Text OpenURL

  38. Kalscheuer R, Stöveken 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-1379. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Metzger J, Bornscheuer U: Lipids as Renewable Resources: Current State of Chemical and Biotechnological Conversion and Diversification.

    Applied Microbiology and Biotechnology 2006, 71:13-22. PubMed Abstract | Publisher Full Text OpenURL

  40. Alvarez HM: Biotechnological Production and Significance of Triacylglycerols and Wax Esters. In Handbook of Hydrocarbon and Lipid Microbiology. Edited by Timmis KN. Springer Berlin Heidelberg; 2010:2995-3002. OpenURL

  41. International Chicken Genome Sequencing Consortium: Sequence and Comparative Analysis of the Chicken Genome Provide Unique Perspectives on Vertebrate Evolution.

    Nature 2004, 432:695-716. PubMed Abstract | Publisher Full Text OpenURL

  42. International Chicken Polymorphism Map Consortium: A Genetic Variation Map for Chicken with 2.8 Million Single-Nucleotide Polymorphisms.

    Nature 2004, 432:717-722. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Oelkers P, Behari A, Cromley D, Billheimer JT, Sturley SL: Characterization of Two Human Genes Encoding Acyl Coenzyme A:Cholesterol Acyltransferase-Related Enzymes.

    J Biol Chem 1998, 273:26765-26771. PubMed Abstract | Publisher Full Text OpenURL

  44. Holmes RS: Comparative Genomics and Proteomics of Vertebrate Diacylglycerol Acyltransferase (DGAT), Acyl CoA Wax Alcohol Acyltransferase (AWAT) and Monoacylglycerol Acyltransferase (MGAT).

    Comp Biochem Physiol Part D Genomics Proteomics 2010, 5:45-54. PubMed Abstract | Publisher Full Text OpenURL

  45. Stone SJ, Levin MC, Farese RV Jr: Membrane Topology and Identification of Key Functional Amino Acid Residues of Murine Acyl-CoA:Diacylglycerol Acyltransferase-2.

    J Biol Chem 2006, 281:40273-40282. PubMed Abstract | Publisher Full Text OpenURL

  46. Hellenbrand J, Biester E-M, Gruber J, Hamberg M, Frentzen M: Fatty Acyl-CoA Reductases of Birds.

    BMC Biochemistry 2011, 12:64. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  47. Tang BY, Hansen IA: Wax Ester Synthesis in the Uropygial Glands of the Chicken and Turkey.

    Comp Biochem Physiol B 1976, 53:451-453. PubMed Abstract | Publisher Full Text OpenURL

  48. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic Local Alignment Search Tool.

    J Mol Biol 1990, 215:403-410. PubMed Abstract OpenURL

  49. National Center for Biotechnology Information [http://blast.ncbi.nlm.nih.gov/] webcite

  50. Wagner M, Hoppe K, Czabany T, Heilmann M, Daum G, Feussner I, Fulda M: Identification and Characterization of an Acyl-CoA:Diacylglycerol Acyltransferase 2 (DGAT2) Gene from the Microalga O. tauri.

    Plant Physiol Biochem 2010, 48:407-416. PubMed Abstract | Publisher Full Text OpenURL

  51. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al.: Clustal W and Clustal X Version 2.0.

    Bioinformatics 2007, 23:2947-2948. PubMed Abstract | Publisher Full Text OpenURL

  52. Nicholas KB, Nicholas HB, Deerfield DW: GeneDoc: Analysis and Visualization of Genetic Variation.

    EMBNEW News 1997, 4:14. OpenURL

  53. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods.

    Molecular Biology and Evolution 2011, 28:2731-2739. PubMed Abstract | Publisher Full Text OpenURL

  54. Saitou N, Nei M: The Neighbor-Joining Method: A New Method for Reconstructing Phylogenetic Trees.

    Molecular Biology and Evolution 1987, 4:406-425. PubMed Abstract | Publisher Full Text OpenURL

  55. Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein Identification and Analysis Tools on the Expasy Server. In The Proteomics Protocols Handbook. Edited by Walker JM. Humana Press; 2005:571-607. OpenURL

  56. ExPASy Bioinformatics Resource Portal [http://web.expasy.org/protparam/] webcite

  57. Sonnhammer EL, von Heijne G, Krogh A: A Hidden Markov Model for Predicting Transmembrane Helices in Protein Sequences. In Proceedings of the Sixth International Conference on Intelligent Systems for Molecular Biology. Edited by Glasgow JLT, Major F, Lathrop R, Sankoff D, Sensen C. AAAI Press, Menlo Park, CA; 1998:175-182. OpenURL

  58. Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting Transmembrane Protein Topology with a Hidden Markov Model: Application to Complete Genomes.

    J Mol Biol 2001, 305:567-580. PubMed Abstract | Publisher Full Text OpenURL

  59. Möller S, Croning MDR, Apweiler R: Evaluation of Methods for the Prediction of Membrane Spanning Regions.

    Bioinformatics 2001, 17:646-653. PubMed Abstract | Publisher Full Text OpenURL

  60. TMHMM Server v 2.0 [http://www.cbs.dtu.dk/services/TMHMM/] webcite

  61. Kyte J, Doolittle RF: A Simple Method for Displaying the Hydropathic Character of a Protein.

    J Mol Biol 1982, 157:105-132. PubMed Abstract | Publisher Full Text OpenURL

  62. NCBI Conserved Domain Search [http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd] webcite

  63. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, et al.: CDD: A Conserved Domain Database for the Functional Annotation of Proteins.

    Nucleic Acids Res 2011, 39:D225-D229. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Bligh EG, Dyer WJ: A Rapid Method of Total Lipid Extraction and Purification.

    Can J Biochem Physiol 1959, 37:911-917. PubMed Abstract | Publisher Full Text OpenURL

  65. Dunphy PJ, Whittle KJ, Pennock JF: On the Use of Fluorescein and Dichlorofluorescein as Non-Destructive Stains for Lipids.

    Chem Ind 1965, 27:1217-1218. PubMed Abstract OpenURL

  66. Gerson T, Schlenk H: 2-Methyl Substituted and 2-Enoic Fatty Acids by Favorsky Rearrangement of Haloketones.

    Chem Phys Lipids 1968, 2:213-219. PubMed Abstract | Publisher Full Text OpenURL

  67. Biester E-M, Hellenbrand J, Frentzen M: Multifunctional Acyltransferases from Tetrahymena thermophila.

    Lipids 2011.

    DOI 10.1007/s11745-011-3642-1

    OpenURL

  68. Bradford MM: A Rapid and Sensitive Method for the Quantitation of Microgram Quantities of Protein Utilizing the Principle of Protein-Dye Binding.

    Anal Biochem 1976, 72:248-254. PubMed Abstract | Publisher Full Text OpenURL