|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 5, 2996-3003, February 2, 2007
Co-transcribed Genes for Long Chain Polyunsaturated Fatty Acid Biosynthesis in the Protozoon Perkinsus marinus Include a Plant-like FAE1 3-Ketoacyl Coenzyme A Synthase*From the Crop Performance and Improvement Division, Rothamsted Research, Harpenden, Herts AL5 2JQ, United Kingdom
Received for publication, July 25, 2006 , and in revised form, November 7, 2006.
The marine parasitic protozoon Perkinus marinus synthesizes the polyunsaturated fatty acid arachidonic acid via the unusual alternative 8 pathway in which elongation of C18 fatty acids generates substrate for two sequential desaturations. Here we have shown that genes encoding the three P. marinus activities responsible for arachidonic acid biosynthesis (C18 9-elongating activity, C20 8 desaturase, C20 5 desaturase) are genomically clustered and co-transcribed as an operon. The acyl elongation reaction, which underpins this pathway, is catalyzed by a FAE1 (fatty acid elongation 1)-like 3-ketoacyl-CoA synthase class of condensing enzyme previously only reported in higher plants and algae. This is the first example of an elongating activity involved in the biosynthesis of a polyunsaturated fatty acid that is not a member of the ELO/SUR4 family. The P. marinus FAE1-like elongating activity is sensitive to the herbicide flufenacet, similar to some higher plant 3-ketoacyl-CoA synthases, but unable to rescue the yeast elo2 /elo3 mutant consistent with a role in the elongation of polyunsaturated fatty acids. P. marinus represents a key organism in the taxonomic separation of the single-celled eukaryotes collectively known as the alveolates, and our data imply a lineage in which ancestral acquisition of plant-like genes, such as FAE1-like 3-ketoacyl-CoA synthases, occurred via endosymbiosis. The P. marinus FAE1-like elongating activity is also indicative of the independent evolution of the alternative 8 pathway, distinct from ELO/SUR4-dependent examples.
The parasitic protozoon Perkinsus marinus is an alveolate (the collective name for ciliates, dinoflagellates, and apicomplexans) that causes severe mortality in populations of the American eastern oyster Crassostrea virginica (1). As a consequence of infection by P. marinus (commonly known as "Dermo" based on its earlier classification as a fungal species, Dermocystidium sp.) (2), commercial production of native oysters in Virginian sections of the Chesapeake Bay has declined drastically (1). Although attempts have been made to identify and breed Dermo-resistant strains of C. virginica, the basis on which P. marinus is such an effective and prevalent pathogen of the eastern oyster is unclear (3). One unusual aspect of P. marinus is its use of the so-called alternative 8 pathway to synthesize long chain polyunsaturated fatty acids (LC-PUFAs),3 such as arachidonic acid (ARA) (20:4 5,8,11,14n-6 (4, 5), which are well known precursors of the eicosanoid class of signaling molecules (such as prostaglandins, thromboxanes, and leukotrienes) (5). It has been hypothesized that fatty acids such as ARA play a key role in the free-living (meront) infectious stage of the P. marinus lifecycle (6). There is ongoing debate as to the prevalence of the alternative 8 pathway in which 18:2n-6 and 18:3n-3 are elongated to generate C20 fatty acid substrates for a specific 8 desaturase (from which the pathway derives its name) followed by 5 desaturation. This 8 pathway therefore differs from the predominant conventional pathway found in the vast majority of LC-PUFA-synthesizing organisms, where 18:2n-6 and 18:3n-3 undergo 6 desaturation followed by elongation prior to 5 desaturation (7). Currently, only three unrelated unicellular organisms have been demonstrated to contain genes encoding the alternative 8 pathway (Euglena gracilis, Acanthamoeba castellanii, and Isochrysis galbana) (8, 9). Metabolic evidence for this 8 pathway in mammals, in particular diseased or cancerous cells, has been presented but also disputed (reviewed in Ref. 10). More recently, the presence of the alternative pathway in P. marinus has also been demonstrated biochemically with the free-living meront stage of this protozoon synthesizing arachidonic acid from the C18 precursor oleic acid via elongation of 18:2n-6 (11). In addition to its commercial significance as a pathogen of shellfish, P. marinus represents a key species in the taxonomic division of the aveolates. In some taxonomies, P. marinus has been accorded its own branch within the aveolatae (the perkinsozoa) (12), but a sublineage more closely related to dinoflagellates has also been proposed on the basis of multiple molecular phylogenies (13). However, the evolutionary path of this organism is still unclear and the basis of ongoing debate.
A number of recent studies have attempted to produce LC-PUFAs such as ARA in transgenic plants with the aim of providing an alternative, sustainable source of these important fatty acids (14, 15). In particular, the use of the alternative
Growth and Harvesting of P. marinusP. marinus cells were obtained from the American Type Culture Collection (ATCC) (Manassas, VA). P. marinus meronts were grown in ATCC Medium 1886 and prepared as described by the distributor. Meronts were inoculated ( 1 x 106 ml-1) and cultivated in 10-ml aliquots of medium in T-10 tissue culture flasks at 25 °C. Meronts at exponential growth phase (7 or 9 days old) were harvested by centrifugation at 200 x g for 5 min. Nucleic Acid Manipulation and PCR-based CloningTotal RNA was extracted from cells using an RNeasy mini kit (Qiagen). First strand cDNA was synthesized from total RNA using the SMART RACE cDNA amplification kit (BD Biosciences) according to the manufacturer's instructions and used for PCR-based cloning. Single-stranded cDNAs were amplified with primers specific to each desaturase gene as follows. The reactions were heated to 95 °C for 2 min followed by 30 cycles at 94 °C for 30 s, 30 s at temperatures ranging from 55 to 72 °C according to the primer design and 72 °C for 2 min, and then a single step at 72 °C for 10 min. PCR amplification products were cloned into TOPO vector (Invitrogen) and verified by sequencing. cDNAs were amplified with specific primers to synthesize full-length copies of the putative desaturases. Gene-specific primers were designed to the 5' and 3' ends of the coding regions of the corresponding nucleotide sequences with restriction sites to facilitate cloning into the yeast expression vector. Forward primers for cloning into yeast expression vector pYES2 (Invitrogen) were designed to contain a G/A at position -3 and a G at position +4 to improve translation initiation in eukaryotic cells. Sequence of PCR Oligonucleotides Used in This WorkThe following gene-specific primers were used: Elo2For (5'-ATGCAAGTTCCCGCGGAGCATCACTCC-3'), Elo2Rev (5'-GTTACGCATCAATATTATGCATAGCCAACC-3'), Des1 A for (5'-ATGACTACTTCAACCACTACTGTGCA-3'), Des1Arev (5'-ACTAGGATCGTTCGTTAGAGAGCTAC-3'), Des2Afor (5'-ATGTCTTCTCTTACCCTCTACAGAGGC-3'), and PemaDes2Rev (5'-CTATTCCACTATGGCAACAGCAGCCA-3'). The following yeast expression primers were used (restriction sites are shown in bold): KpnElo2For (5'-TTGGTACCATGGGATTTCCTGCGGAG-3'), SacElo1Rev (5'-GGAGCTCTTACGCATCAATATTATGCATAGC-3'), KpnDes1For (5'-CCGGTACCATGGCTACTTCAACCACTACTGTGC-3'), SacDes1Rev (5'GGGAGCTCCTACCTAGCAAGCAATCTCTCGATGTTACGC-3'), KpnDes2For (5'CCGGTACCATGGCTTCTCTTACCCTCTACAGAGGC-3'), and SacDes2Rev (5'GGGAGCTCCTATTCCACTATGGCAACAGCAGCC-3'). The following operon amplification primers were used: FAEfor (5'-ATCCAAGCAAGCCGTAGTGT-3'), D8rev (5'-ACCCTCGCATCCTCCTTAGT-3'), D8for (5'-TAGCCCAGGAGCTAGGACTG-3'), and D5rev (5'-TCACTGCCAGCATGAAACTC-3'). Controls for gDNA contamination were: FAE225f (5'-TTCCGGCACCAGTGACCATACGG-3') and FAE225r (5'-CTGTTCTGGAGAAGACGTTGAGC-3'). Transcript analysis of the P. marinus OperonTotal RNA was isolated from free-living meronts and subjected to treatment with RNase-free DNaseI. RNA was repurified and used as a template for cDNA synthesis with an oligo(dT) primer. A control reaction was carried out in which no reverse transcriptase was added. Resulting cDNA was used as a template for PCR with two gene-specific primer pairs (Faefor, D8rev and D8for, D5rev). Each pair spanned the intervening sequence between two open reading frames (ORFs) (illustrated in Fig. 3). As an additional control, we used a second FAE1-like gene (present on TIGR P. marinus clone 22547), which was predicted to contain three small introns. This allowed us to confirm the absence of genomic contamination in cDNA derived from DNaseI-treated RNA. Functional Expression in YeastORFs encoding putative desaturation activities were introduced in Saccharomyces cerevisiae strain W303-1A by lithium acetate. Cultures were grown at 22 °C in the presence of 2% (v/v) raffinose, and expression of the transgenes was induced by the addition of galactose to 2% (w/v) in the presence of the corresponding fatty acid at 0.5 mM and 1% (w/v) tergitol-Nonidet P-40 (Sigma).
Flufenacet Inhibition of Elongation in YeastFlufenacet (20500 µM) (Sigma-Aldrich) was added to 25 ml of yeast culture transformed with pYES2PmFAE at the time of induction. After 48 h of incubation time, the yeast cells were processed as described above, and fatty acid composition was compared with cell-expressing PmFAE grown in the absence of herbicide. Data presented represent the mean of three independent experiments with two technical replicates for each concentration tested. Inhibition of the herbicide was expressed at percentage of inhibition of the elongation of C18
Yeast Mutant ComplementationThe P. marinus FAE ORF was cloned into pADH for expression and complementation assays. This construct was introduced in elo2 Fatty Acid AnalysisFatty acids were extracted and methylated using standard protocols. Total fatty acids were analyzed by gas chromatography of methyl ester derivatives.
The genome sequence of P. marinus was examined for the presence of genes likely to encode LC-PUFA biosynthetic activities. Searches of the shotgun-sequenced unassembled reads were carried out via the WU-BLAST2 tool on The Institute for Genomic Research (TIGR) P. marinus genome site using N-terminal cytochrome b5 fusion desaturases containing Pfam domains 00173 and 00487 and ELO/SUR4-elongating activity (containing Pfam domain 01151) as templates. These two components (desaturase, ELO-like elongating activity) are absolutely conserved in all known examples of the aerobic PUFA biosynthetic pathway (14, 18). This search identified one genomic fragment (TIGR designation 1047306867) of 18 kb, which reported significant similarity to N-terminal cytochrome b5 fusion desaturases associated with LC-PUFA biosynthesis. Closer inspection of the deduced amino acid sequence of this particular clone in fact revealed the presence of two cytochrome b5 fusion desaturases in close proximity to one another. The deduced ORFs of these two putative LC-PUFA desaturases were separated by 117 bp, and neither gene appeared to contain introns. Although clearly related to each other, the two putative desaturases were only 29% identical, 40% similar (Fig. 1A). However, they both contained the diagnostic histidine boxes associated with LC-PUFA desaturases (including the variant Gln His substitution in the third box) as well as an N-terminal cytochrome b5 domain containing the highly conserved heme binding motif His-Pro-Gly-Gly (Fig. 1A) (19). These two ORFs were placed in a phylogenetic tree of known cytochrome b5 fusion desaturases as a tool to infer biochemical function (Fig. 1B). This failed to reveal any obvious lineage but grouped the two P. marinus cytochrome b5 fusion desaturases with 5 and 6 desaturases from the parasitic protozoon Leishmania major (20) and the algae Thraustochytrium sp. and Ostreococcus tauri (21). In contrast, the higher plant cytochrome b5 fusion desaturases clearly form a distinct branch (19). As can also be seen, neither of the P. marinus desaturases showed any close phylogenetic association with the two known examples of the alternative 8 pathway C20 8 desaturases from E. gracilis and A. castellanii (9) (Fig. 1B).
Further analysis of the sequence of clone 1047306867 identified one additional ORF also in close proximity (150 bp upstream) to these two desaturase-like genes. This ORF showed similarity to condensing enzymes involved in the biosynthesis of long chain fatty acids, indicating the genomic clustering of the two distinct classes of ORFs involved in this biochemical process (desaturases, elongating activity) and their organization in a distinct linear manner (illustrated in Fig. 2A). Perhaps more surprisingly, this putative condensing enzyme was not a member of the ELO/SUR4 family normally associated with this pathway but rather a member of the FAE1-like (22) class of 3-ketoacyl-CoA synthase (KCS) activities (Pfam cd00831.2), which until now have only been functionally described in higher plants and photosynthetic algae (and as such, have been assumed to be specific to the plant kingdom) (23). The FAE1 gene of Arabidopsis thaliana has been demonstrated biochemically to be a KCS microsomal condensing enzyme distantly related to chalcone synthase, and genetic analysis has confirmed FAE1 is specifically involved in the synthesis of C20+-saturated and -monounsaturated fatty acids (21, 23). Using the P. marinus FAE1-like ORF to search genomic databases at NCBI, TIGR, and Sanger, we confirmed the absence of FAE1-like sequences in metazoan, yeast, and fungal genomes but identified several such sequences in human pathogens, such as the parasitic protozoon Entamoeba histolytica and the soil-living amoeba Dictyostelium discoideum. We also observed the presence of FAE1-like sequences in photosynthetic algae, such as Phaeodactylum tricornutum, in agreement with earlier reports of a FAE1-like condensing enzyme identified from the alga Dunaliella salina (24). A phylogenetic tree comparing the P. marinus FAE1-like sequence with the 21 KCS FAE1-like orthologs present in the Arabidopsis genome and also amoeba and algal FAEs is shown in Fig. 2B. This indicates that the closest association with the P. marinus FAE1-like sequence is shown by the P. tricornutum FAE.
Given the operon-like organization of the three P. marinus ORFs predicted to be involved in the synthesis of LC-PUFAs, it was determined whether these three genes were co-transcribed or instead generated individual transcripts. Polycistronic transcription has been documented to occur in parasitic protozoa such as trypanosomes (25). No information exists regarding P. marinus transcriptional systems, although the proximity of the three genes makes such a process likely. Therefore, two primer pairs were designed to (a) the 3' end of the FAE1-like ORF and 5' end of the first desaturase ORF and (b) the 3' end of the first desaturase ORF and the 5' end of the second desaturase ORF (Fig. 3). Thus, if these three genes were co-transcribed as a single transcript, it would be possible to generate two discrete PCR products (of 461 and 424 bp). As shown in Fig. 3, two such products were amplified with the identity of these PCR products confirmed by cloning and sequencing. In view of the possibility that these PCR products were derived from genomic DNA still present in the cDNA template used for amplification, control primers were designed to another (uncharacterized) FAE-like gene we identified on P. marinus TIGR clone 22547 (hence designated FAE225a). This FAE-like gene contained three small introns (totaling
To confirm the predicted function of the P. marinus ORFs in the biosynthesis of arachidonic acid, the individual coding sequences for the three genes were cloned into yeast expression vectors. As shown in Fig. 4, this functional characterization in yeast demonstrated LC-PUFA biosynthetic activities for these P. marinus ORFs. The first ORF, the FAE1-like KCS, was shown to encode an elongating activity for 9-desaturated C18 fatty acids, converting 18:2 9,12n-6 and 18:3 9,12,15n-3 into 20:2 11,14n-6 and 20:3 11,14,17n-3 (the latter substrate is shown in Fig. 4B). Such activity represents the first reaction in the alternative 8 pathway for PUFA biosynthesis (also detailed in Fig. 4), but the only previously identified C18 9-elongating activity was the structurally unrelated ELO/SUR4-like ORF from I. galbana (26). The P. marinus C18 9-elongating activity (hereafter designated PmFAE) was tested with other substrates such as -linolenic acid and stearidonic acid but found to have no activity toward C18 6-desaturated fatty acids. The first cytochrome b5 fusion desaturase PmD8, which lies adjacent to the 3' end of the PmFAE gene, was shown to encode a C20 8 desaturase, the activity associated with the alternative pathway of PUFA biosynthesis. This enzyme utilized 20:2 11,14n-6 and 20:3 11,14,17n-3 to yield the desaturation products 20:3 8,11,14n-6 and 20:4 8,11,14,17n-3. When this desaturase was co-expressed with PmFAE (generating the n-3 C20 substrate in vivo), this latter fatty acid (20:4n-3) was observed (Fig. 4C) and its identity confirmed by gas chromatography/mass spectrometry (data not shown). The second cytochrome b5 fusion desaturase PmD5 was identified as a C20 5 desaturase, (Fig. 4D) responsible for the conversion of the C20 products of the 8 desaturase to ARA (n-6) and eicosapentanoic acid (20:5 5,8,11,14,17n-3; eicosapentaenoic acid), respectively (the latter product is shown in Fig. 4D). Fig. 4D not only confirms the activity of PmD5 as a C20 5 desaturase but also demonstrates that co-expression of all three P. marinus ORFs (PmFAE, PmD8, and PmD5) resulted in the successful heterologous reconstitution of the alternative LC-PUFA biosynthetic 8 pathway. The substrate preference of these three activities was also tested by individual expression of these activities in yeast. As can be seen in Table 1, although PmFAE and PmD8 showed a slight preference for n-3 substrates, PmD5 showed preference for n-6 substrates.
Based on the above results, the three clustered P. marinus genes encode all of the primary biosynthetic activities required for the synthesis of ARA and eicosapentaenoic acid via the alternative 8 pathway. Because these three genes are co-transcribed (Fig. 3), they represent a PUFA biosynthetic operon. Moreover, the three genes are structurally positioned in the same order as their encoded biochemical reactions proceed in the alternative 8 pathway (i.e. elongation, 8 desaturation, and 5 desaturation); however, whether this represents some aspect of selection pressure is currently unclear.
The functional identification and characterization of this P. marinus PUFA biosynthetic operon raises a number of intriguing points. First, the involvement of a FAE1-like KCS in LC-PUFA biosynthesis is remarkable, given the taxonomic distribution of FAE1-like orthologs. Moreover, these FAE1-like genes normally found only in photosynthetic organisms are required for the synthesis of saturated and monounsaturated long chain fatty acids primarily used in storage lipids and cuticular waxes (27), with no previously characterized examples shown to be involved in the elongation of polyunsaturated fatty acids (23). This is therefore the first observation of the aerobic biosynthesis of LC-PUFAs using a FAE1-like KCS-elongating activity rather than the ELO/SUR4 system. Second, we confirmed the presence of the alternative
Equally intriguing is the genomic organization of the P. marinus alternative pathway genes as a structured operon, which could indicate the co-evolution of all three genes. The presence of the FAE1-like KCS points toward an ancestral gene transfer from a photosynthetic alga, also implied by the phylogenetic association with the P. tricornutum FAE. Whether such a genetic exchange would involve the nuclear transfer of the entire operon or just the FAE1-like ORF is unknown, although this particular gene has significantly diverged in function from other FAE1-like condensing enzymes. Because this is the first description of all three activities of the alternative
Further evidence for the unique nature of the P. marinus FAE1-like elongating activity was provided by the inability of PmFAE to complement yeast mutants lacking the ELO2p- and ELO3p-elongating activities and therefore defective in the synthesis of saturated very long chain fatty acids present in sphingolipids. It has been shown that higher plant FAE1-like KCS ORFs can rescue the otherwise lethal yeast elo2 /elo3 double mutant, even though the two enzymes are structurally unrelated (17). However, as can be seen in Fig. 5, the PmFAE was unable to substitute for either ELO2 or ELO3, in contrast to two Arabidopsis genes At2g16280 (CUT1) and At4g34520 (FAE1), both involved in the elongation of C18-saturated and -monounsaturated fatty acids. In that respect, the PmFAE 9-elongating activity behaves in a similar manner as that observed for the ELO/SUR4-like PUFA-elongating activity from Caenorhabditis elegans PEA1 (7), which also fails to rescue the elo2 /elo3 double mutant (Fig. 5). This indicates that the primary function of the P. marinus FAE is likely to be the biosynthesis of LC-PUFAs as opposed to sphingolipids and highlights its functional divergence from higher plant FAE-like genes.
The observation of PUFA biosynthesis catalyzed by a FAE1-like KCS not only indicated the presence of a plant-like elongation system in P. marinus but also suggested a possible target for chemical inhibition. Recent studies on the FAE1-like KCSs of Arabidopsis have shown that the K3 class of anilide herbicides can act as potent inhibitors of the very long chain fatty acid elongation reactions catalyzed by these condensing enzymes (29). Expression in yeast of several Arabidopsis FAE1-like KCS activities were inhibited by K3 herbicides, such as flufenacet (N-(4-fluorophenyl)-N-(1-methylethyl)-2-(5-trifluoromethyl-1,3,4-thiadiazol-2-yl oxy)acetamide), although this compound had no effect on the endogenous yeast ELO/SUR4 elongases (30). We therefore examined the sensitivity of PmFAE to flufenacet using similar methods, examining the inhibition of elongation of substrate fatty acids such as 18:2 9,12 and 18:3 9,12,15. In view of the divergent nature of PmFAE (both in terms of function and primary sequence), it was not obvious whether the effect of anilide herbicides on this protist enzyme would be similar to that observed for the higher plant enzymes. However, application of flufenacet (in the range 50500 µM) to yeast cultures expressing PmFAE resulted in a strong inhibition of the C18
As mentioned above, the alternative
In conclusion, we have identified a PUFA biosynthetic operon that encodes for all of the primary enzymatic components of the alternative
The nucleotide sequence(s) reported in this paper has been submitted to the GenBankTM/EBI Data Bank with accession number(s) DQ508730, DQ508731, and DQ508732.
* This work was supported by a grant-in-aid from the Biotechnology and Biological Sciences Research Council (UK). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 Partially supported by the European Union FP6 Project LIPGENE (FOOD-CT-2003-505944). 2 To whom correspondence should be addressed. Tel.: 44-1582-763133; Fax: 44-1582-763010; E-mail: johnathan.napier{at}bbsrc.ac.uk.
3 The abbreviations used are: LC-PUFA, long chain polyunsaturated fatty acid; ORF, open reading frame; ARA, arachidonic acid; KCS, 3-ketoacyl-CoA synthase; FAE, fatty acid elongating activity; gDNA, genomic DNA.
We thank John Pickett for critical reading of the manuscript and BASF Plant Sciences (Limburgerhof, Germany) for support. We thank Teresa Dunn for generously providing yeast strains used in this study. Preliminary sequence data were obtained from The Institute for Genomic Research. Sequencing of P. marinus was accomplished with support from National Science Foundation.
This article has been cited by other articles:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||