Molecular mechanism of substrate speciﬁ city for delta 6 desaturase from Mortierella alpina and Micromonas pusilla

The (cid:2) 6 and (cid:2) 3 pathways are two major pathways in the biosynthesis of PUFAs. In both of these, delta 6 desaturase (FADS6) is a key bifunctional enzyme desaturating linoleic acid or (cid:3) -linolenic acid. Microbial species have different propensity for accumulating (cid:2) 6- or (cid:2) 3-series PUFAs, which may be determined by the substrate preference of FADS6 enzyme. In the present study, we analyzed the molecular mechanism of FADS6 substrate speciﬁ city. FADS6 cDNAs were cloned from Mortierella alpina (ATCC 32222) and Micromonas pusilla (CCMP1545) that synthesized high levels of arachidonic acid and EPA, respectively. M. alpina FADS6 (MaFADS6-I) showed substrate preference for LA; whereas, M. pusilla FADS6 (MpFADS6) preferred ALA. To understand the structural basis of substrate specificity, MaFADS6-I and MpFADS6 sequences were divided into ﬁ ve sections and a domain swapping approach was used to ex-amine the role of each section in substrate preference. Our results showed that sequences between the histidine boxes I and II played a pivotal role in substrate preference. Based on our domain swapping results, nine amino acid (aa) residues were targeted for further analysis by site-directed mutagenesis. substitutions substrate recognition, corresponding aa

level of GLA or AA is high ( 10 ). This suggests that their FADS6s may have preference for LA. The genus Mortierella has been extensively studied for its production of GLA ( 13 ) or AA (14)(15)(16). Some species of this genus are now used for the commercial production of single cell oil that is rich in AA ( 17 ). Among them, Mortierella alpina (ATCC 32222), whose genome has previously been characterized in our lab ( 18 ), has the ability to synthesize a wide range of PUFAs (up to 50% of its cell dry weight), including AA as high as 40% of total fatty acid. By contrast, its -3 PUFA level is low ( 18 ). Therefore, M. alpina was used as the other source of FADS6.
Although M. alpina has a disproportionate ratio of -6 to -3 PUFAs ( 19 ), that does not necessarily mean that its FADS6 (MaFADS6) has substrate preference for LA. Previous studies reported that M. alpina had low levels of ALA, SDA, and EPA ( 19 ), and expressed no delta 9 elongase ( Fig. 1 ), as previously demonstrated by M. alpina genome determination ( 18 ). Therefore, even if MaFADS6 had equivalent catalytic effi ciency toward LA and ALA, it would preference of FADS6 for LA and ALA is likely to be different in species which produce greatly divergent levels of -3 and -6 PUFAs. For example, in a small number of species, especially in marine microbes, EPA or DHA production is much higher than production of -6 PUFAs, probably because their FADS6s preferentially use ALA ( 12 ). The catalytic efficiency of FADS6s from various species with LA and ALA substrates is summarized in supplementary Table 1. From the results, we found that Micromonas pusilla , a type of marine microalgae whose FADS6 (MpFADS6) appears to function as an acyl-CoA desaturase, has the highest preference for the -3 substrate (ALA) compared with other plants and fungi. Its conversion effi ciency was found to be 4.9% in the -6 pool and 63% in the -3 pool ( 12 ). Mp-FADS6 has the highest conversion ratio of ALA/LA by far. Therefore, M. pusilla was used as one of the sources of FADS6 in our study.
In contrast, the EPA or DHA production is fairly low or even absent in numerous oleaginous species, while the  ( 34 ) including microorganisms, plants, and animals . The various routes for biosynthesis of AA, EPA, and DHA are shown, as mediated by the consecutive action of desaturases and elongases. The predominant delta 6 pathway is shown by black boxes with white bold font, as is the alternative D8 pathway. Three routes for DHA biosynthesis are shown: microbial D4 pathway, mammalian "Sprecher" pathway, and polyketide synthase (PKS) pathway [PUFAs are synthesized by polyketide synthase controlling chain length ( 35 )]. FADS12, FADS15, and FADS17 are found in plants, but not in mammals. FADS4 is also only found in plants and in marine teleost fi sh species ( 36 ). There is no alternative delta 8 pathway or microbial delta 4 pathway in M. alpina . ELOVL, elongation of long-chain fatty acids; FADS9, delta 9 desaturase; ELO9, delta 9 elongase; FADS12, delta 12 desaturase; FADS15, delta 15  GenBank accession number ADAG00000000 ([fi rst version)] and MpFADS6 gene sequences . The forward primers were FMa and FMp and the reverse primers were RMa and RMp (all primers are listed in Table 1 ). PCR was carried out in a total volume of 50 l. After initial denaturation at 94°C for 4 min, amplifi cation was performed in 30 cycles of 40 s at 94°C, 40 s at 55°C, and 2 min at 68°C, followed by a fi nal extension at 68°C for another 5 min. Amplifi cation products were fractionated on 1.0% agarose gels and subcloned into the pYES2 vector (Invitrogen) downstream of the GAL1 promoter to generate recombinant plasmids designated pYES2-MaFADS6 and pYES2-MpFADS6. Then the constructs were transformed into competent E. coli cells and positive clones were selected by colony PCR with T7 and pYES2.R primers. Positive clones were sequenced in both directions.

Construction of chimeric genes by reciprocal section swapping
According to the conserved regions of FADS6 (HPGG and three His boxes), MaFADS6 and MpFADS6 genes were divided into fi ve regions ( Fig. 3A ). The corresponding regions of both enzymes were systematically exchanged to construct recombinant swap genes in order to determine enzymatic specifi city of each fragment. Recombinant swap genes were generated by overlap extension PCR with the primers listed in Table 1 . In the fi rst step, each fragment was amplifi ed individually from the Ma-FADS6 and MpFADS6 cDNA template using primers specifi ed in Table 1 . In the second step, two adjacent fragments that had homologous arms were fused using upstream and downstream primers to produce double fragments if the recombinant swap gene had three fragments. Similarly, double fragments and a complementary fragment were fused by the same method (with the longer one generally selected). The amplifi cation procedure included 10 min at 94°C, 2 min at 60°C, and 1 min per kilobase at 68°C, followed by 35 cycles of 50 s at 94°C, 50s at 60°C, and 1 min per kilobase at 68°C, and a terminal extension step of 5 min at 68°C. Amplifi cation products were fused in a subsequent PCR reaction, and the protocol used was the same as for the fi rst step. The hybrid genes were digested with Eco R I/ Not I and ligated into the pYES2 plasmid (Invitrogen) with a His tag sequence, and transformed into competent Top10 cells .

Site-directed mutagenesis of MpFADS6
For mutagenesis, oligonucleotide primers were used to introduce nucleotide substitutions into MpFADS6 using the Fast Site-Directed Mutagenesis kit (TIANGEN) according to the manufacturer's instructions. Seven mutants derived from the MpFADS6 gene were constructed, where aas were substituted with the corresponding residues in the MaFADS6 gene (V189L/Q190A, G194L, S197Q, Q209G, E222S, M227K, and V399I/I400E). Mutation primers are listed in Table 1 . All mutants were verifi ed by colony PCR with T7 and pYES2.R primers. Positive mutants were sequenced in both directions.
Yeast transformation, heterologous expression in S. cerevisiae , and determination of substrate preference for LA or ALA by measuring chimeric desaturase activity Constructs pYES2-chimera 1 ‫ف‬ 10, pYES2-MaFADS6, and pYES2-MpFADS6 were transformed into S. cerevisiae using the lithium acetate transformation method ( 27 ). After SC-U plate selection, the recombinant yeasts were selected on uracil-defi cient medium containing 1.0% (w/v) raffi nose as the single carbon source. MaFADS6 and MpFADS6 genes were induced under transcriptional control of the GAL1 promoter containing 2% galactose for 12 h at 28°C.
The pellets were collected by centrifuging 2 ml culture at 6,000 g for 1 min. Yeast pellets were resuspended in 80 l breaking not result in an equal level of its products, GLA and SDA (or longer-chain PUFAs), because initial levels of ALA are low. So the high levels of GLA and AA in M. alpina do not automatically imply that MaFADS6 has a strong substrate preference for LA.
Recently, some genes encoding FADS6 have been cloned and sequenced from many organisms ( 13,(20)(21)(22)(23)(24)(25)(26), but little is known about their molecular mechanism of FADS6 substrate preference. This is mainly due to the lack of crystal structure information. In this work, we cloned the MaFADS6 (MaFADS6 is equivalent to MaFADS6-I) and MpFADS6 genes, and characterized their substrate preference in Saccharomyces cerevisiae . Furthermore, domains were swapped with each other and the recombinant chimeras were expressed in S. cerevisiae to determine key areas responsible for the substrate preference. Finally, site-directed mutagenesis within key areas was used to identify sites that are important for substrate recognition.

Strains and plasmids
Wild-type M. alpina (ATCC 32222) was from our laboratory. The MpFADS6 gene was synthesized by Shanghai Sunny Biotechnology Company, Ltd. and Escherichia coli Top 10 was preserved in our laboratory. INVSc1 yeast strain (Invitrogen) was used for heterogeneous expression, including substrate preference determination and expression of chimeras. Plasmid pYES2/NT C (Invitrogen) was used for FADS6 expression.

RNA isolation and gene synthesis
Approximately 1 g (fresh weight) of M. alpina protonemal tissue was ground to a fi ne powder under liquid nitrogen using a precooled mortar and pestle. Total RNA was isolated from cooled powder using the HYQspin TM Total RNA Kit (HG403-03), and 1 g RNA was reverse-transcribed with the QuantScript RT kit (TIANGEN) according to the manufacturer's instructions. The cDNA transcribed was used as a template for MaFADS6 amplifi cation with primers. In this work, MaFADS6 is equivalent to MaFADS6-I, because the RNA transcriptional level of the MaFADS6-II isoform is quite low ( 16 ).
MpFADS6 gene was synthesized by Shanghai Sunny Biotechnology Company, Ltd. and subcloned into the vector pUC57 and transformed into DH5 ␣ . 0.25 mM cis -ALA, respectively) and 1% Tergitol Nonidet P-40 for the solubilization of fatty acids. Subsequently, cultures were grown for 12 h at 28°C. Cells were harvested by centrifugation, washed with distilled water, and freeze-dried. The dried cells were used to determine fatty acid composition. For the substrate concentration experiment, 0.125, 0.25, or 0.5 mM cis -LA/ALA was added in the induced medium.

Lipid extraction and fatty acid analysis
Lipids from an equivalent weight of freeze-dried cells and supernatant were extracted and methyl esterifi ed as described previously ( 28 ). Fatty acid methyl esters (FAMEs) were analyzed by GC. Pentadecanoic acid (C15:0; NU-CHEK) and heneicosanoic acid (C21:0; NU-CHEK) were added to the biomass samples as internal standards to quantify the fatty acid content. GC analysis was performed with a GC-2010 (Shimadzu Company, Japan) equipped with a flame ionization detector and a capillary DBWAX column (30 m × 0.32 mm, 0.25 m; Agilent ). The samples were measured with a split of 20:1 with the injector temperature set to 240°C. The column temperature was 180°C. FAMEs were identified by comparing with buffer [50 mM sodium phosphate (pH 7.4), 1 mM EDTA, 5% glycerol, 1 mM PMSF] with an equal volume of 0.5 mm glass beads and vibrated for 30 min at 4°C. Yeast lysates were mixed with 20 l 5× SDS sample buffer and heated for 10 min in boiling water, then set on ice. A 20 l sample was loaded onto SDS polyacrylamide gel and was run for about 1.5 h at 120 V. Then, one SDS-PAGE gel was used for Coomassie blue staining and another one for Western blotting analysis.
For Western blotting analysis, gels were transferred onto a polyvinylidene fl uoride membrane by electroblotting (150 mA, 2 h). The membrane was blocked with 5% skim milk in TBST [150 mM NaCl, 10 mM Tris-Cl (pH 7.5), 0.05% Tween 20) and probed with 1:5,000 mouse THETM His tag antibody, followed by 1:2,000 HRP-conjugated goat anti-mouse IgG (H&L). Chemiluminescence reagent was prepared by mixing equal volumes of enhanced luminol reagent and oxidizing reagent. Membranes were then incubated with the chemiluminescence reagent for 1 min with gentle agitation and photos were taken using FluorChem FC3 (ProteinSimple).
After induction, cultures were supplemented with 0.5 mM cis -LA or 0.5 mM cis -ALA or both substrates (0.25 mM cis -LA and  R I  CGCCG GAATTC ATGGCTGCTGCTCCCAGTGTGAGGAC  FD  -CTGATCATCCCGGTGGAAGCGTGATATTCTA  RD  -ATGCAAAAAGTCGTGCTGGACCCAACCGCA  FF  -CACGACTTTTTGCATAACTCGCTCACGGGG  RF  -GTGGTGAGTGTTGTGCATCTGGTTCCACATC  FH  -CACAACACTCACCACGCGACGCCGCAGAA  RH  -GTGGTGCTCGATCTGGCAGTTCAGATATCCC  FJ  -CAGATCGAGCACCACCTGTTCCCGGACAT  RJ  Not I  CGAAGGAAAAAA GCGGCCGC TCAGTGCGCCTTCTCC  RK  -TAGAATATCACGCTTCCACCGGGATGATCAG  FN  -TCACAAACTTCAAACATCCCGGTGGAAGTGTG  RN  -CCCCGTGAGCGAGTTATGCAAAAAGTCGTG  FP  -TGCGGTTGGGTCCAGCACGACTTTTTGCAT  RP  -TTCTGCGGCGTCGCGTGGTGAGTGTTGTG  FR - T TCTGGTTCCACATCTCCCCG  F-V399I/I400E  -GATATCTGAACTGCCAG A TC GAA CATCACCTGTTCCCGGA  R-V399I/I400E  -TCCGGGAACAGGTGATG TTC GA T   order to determine the structural basis of -3/ -6 substrate preference, multiple sequence alignment was performed among several FADS6 genes, arranged according to their substrate preference for LA or ALA (the closer to the top, the stronger preference for ALA; the closer to the bottom, the stronger preference for LA). The result of multiple sequence alignment among these FADS6s revealed that there were three conserved His-rich motifs, HDFLHH (Ma172-177) and HEGGHN (Mp191-196), HDFLHH (Ma209-214) and HNKHH (Mp228-232), and QIEHH (Ma395-399) and QVIHH (Mp398-402) ( Fig. 2 ). In addition, a cytochrome b5-like domain, HPGG, was found near the N terminus, which is required as an electron donor for fatty acid desaturation. Alignment and analysis of MaFADS6 and MpFADS6 sequences with other FADS6s indicated that the homology occurs mainly in the cytochrome b5-like domain and in the three conserved His-rich motif areas. We expected that FADS6s with -3 substrate preference should have a homologous site or fragment not present in FADS6s with -3 substrate preference, but our multiple sequence alignment results failed to uncover such a site or fragment.

Construction and expression of fusion genes in yeast
To identify which structural elements were functionally involved in substrate preference of MaFADS6 and MpFADS6, both enzymes were divided into fi ve sections as follows: section 1, the N-terminal end region minus the HPGG domain; section 2, from the HPGG domain to His box I (including HPGG); section 3, from His box I to His box II (including His box I); section 4, from His box II to His box III (including His box II); and section 5, the C-terminal region including His box III ( Fig. 3A ). After expression in S. cerevisiae, 10 different chimeras were quantifi ed by Western blotting analysis. Our results showed that none of the chimeras were clearly seen by SDS-PAGE, probably because the expression level of membrane-bound proteins in S. cerevisiae was too low to be observed with the naked eye. However, our Western blotting analysis showed that all chimeras were expressed commercial FAME standards (SDA-ME and 37 component FAME mix; Supelco).

Cloning of MaFADS6 and MpFADS6 and characterization of their substrate specifi city
To characterize the substrate specifi city of MaFADS6 and MpFADS6, a 1,374 bp fragment was amplifi ed using FMa and RMa primers, and a 1,392 bp fragment was amplifi ed using FMp and RMp primers (primers listed in Table 1 ). Successful expression of pYES2-MaFADS6 and pYES2-MpFADS6 was confi rmed by Western blotting. They were incubated with 0.5 mM LA, 0.5 mM ALA, or 0.25 mM of each of them as fatty acid substrates for 12 h at 28°C, and the resulting fatty acid composition is listed in Table 2 . When a single substrate was added, LA and ALA could each be catalyzed by MaFADS6 and their conversion rates were 45.6 ± 2.1% and 19.6 ± 1.3%, respectively. When LA and ALA were added at the same time, the LA conversion rate of MaFADS6 reached 58.4 ± 1.4%, but ALA conversion decreased to 2.0 ± 0.9%. These results show that MaFADS6 is capable of catalyzing LA and ALA conversion to GLA and SDA, respectively. However, when LA and ALA were present at the same time, Ma-FADS6 preferentially catalyzed LA to GLA. Similarly, Mp-FADS6 had a preference for the ALA substrate with a 66.5 ± 2.8% conversion rate, which is consistent with the results of Petrie et al. ( 12 ).

Multiple sequence alignment among FADS6s
The deduced aa sequences of MaFADS6 showed 20.35% identity with the aa sequence of MpFADS6. In

Determination of substrate preference by mutating aas from MpFADS6
To further understand the molecular mechanism of MpFADS6 substrate preference, we constructed a series of mutants in those 37 aas within section 3 and a double mutant within section 5. Targeted mutagenesis of V189L/ Q190A, S197Q, and Q209G did not lead to major changes in substrate preference. However, the relative conversion rate of ALA was reduced almost by half in double mutant V399I/I400E and single mutants E222S and M227K (40.26, 31.42, and 31.61%, respectively) compared with MpFADS6 (71.37%). In addition, the relative conversion rate of ALA was signifi cantly reduced to 6.50% in the presence of a G194L substitution. The relative conversion rate for LA remained low in all mutants ( Fig. 5 ).

DISCUSSION
Our result indicates that the 37 aas between the His boxes I and II ( Fig. 6 ) are in part responsible for the substrate preference ( Fig. 3 ). Comparing sequences between the group with preference index (PI) greater than 1 and that with PI smaller than Ϫ 1 within this 37-aa segment, four positions show clear differences, namely 190Q, 194G, 197S, and 209Q/H in the PI >1 group, and 170A, 175L, 178Q, and 190G are the corresponding aa residues in the same position in the PI <1 group. Among the single mutants corresponding to these sites, G194L displayed the largest decrease in ALA conversion rate. In the topological model of the FADS6, 194G is located within His box I, suggesting that 194G is a key residue in His box I that infl uences the ALA conversion to SDA. In the single mutant G194L, steric bulk presented by the larger Leu residue may cause the substrate to bind more loosely. Interestingly, it has been reported that sequences between His boxes I and II are involved in the regioselectivity of -3 and -6 fatty acid desaturases from Pichia pastoris ( 29 ) and Aspergillus nidulans ( 30 ). It has also been shown that this segment is implicated in fatty acid carbon length preference in FADS6 from Mucor rouxii , and the 213S and 218K residues, which correspond to the 222E and 227M residues from MpFADS6 mutants , may be critical ( 31 ). According to the predicted structural model, residues 222 and 227 are located right upstream of the second conserved His box. These two positions were also found to be important for ALA catalysis, but to a lower extent than 194G. Thus, these two residues are unlikely to be located near the binding site of the substrate or to have direct interactions with the ALA substrate. None of and their molecular mass was as predicted ( Fig. 4A and data not shown). The relative expression level of every recombinant protein was fairly stable, with a range of 0.7-2.1 compared with MaFADS6, defined as 1 ( Fig. 4B ).

Determination of substrate preference in 12 types of fusion FADS6 constructs
For substrate preference determination, lipids from dried cells and supernatant were extracted to determine the conversion effi ciency of each chimera (no fatty acids were detected in supernatant, data not shown). To account for differences in chimera expression, substrate conversion rates of all chimeras were converted to relative substrate conversion rates ( Fig. 3C ). Our analysis showed that the catalytic effi ciency of chimera 3 (MaFADS6 aa172 ‫ف‬ 208 replaced by MpFADS6 aa191 ‫ف‬ 227) for ALA increased from 4.0 to 15.0%; whereas, the catalytic efficiency of chimera 3 for LA decreased from 73.9 to 24.0%; and, conversely, the effi ciency of chimera 8 (Mp-FADS6 aa191 ‫ف‬ 227 replaced by MaFADS6 aa172 ‫ف‬ 208) for ALA decreased from 75.2 to 9.5%, but for LA increased from 5.2 to 16.3% compared with their corresponding wild-types. These results suggest that these 37 aas in section 3 are in part responsible for the substrate preference. Replacement of aa1-49 (chimera 1), aa50-171 (chimera 2), aa209-394 (chimera 4), or aa395-458 (chimera 5) in MaFADS6 caused partial or complete loss of activity. Similarly, replacement of aa1-68 (chimera 6), aa228-397 (chimera 9), or aa398-464 (chimera 10) in MpFADS6 resulted in partial or complete loss of activity. Replacement of aa69-190 (chimera 7) in MpFADS6, however, seemed to signifi cantly decrease its activity toward ALA, while slightly increasing its activity toward LA ( Fig. 3C ).

The effect of substrate concentration on the conversion rate
To understand the effect of substrate concentration on the conversion rate, 0.125, 0.25, and 0.5 mM cis -LA/ALA were used in yeast cultures expressing active fatty acid desaturases, i.e., MaFADS6, MpFADS6, and chimeras 1, 3, and 5-9. Our results revealed that substrate concentration had little effect on the conversion rate in yeast expressing MaFADS6, MpFADS6, and chimera 5; whereas, it signifi cantly augmented ALA conversion from 11.0 ± 1.7% to 19.2 ± 3.2% in chimera 3 and LA conversion from 8.7 ± 0.8% to 19.2 ± 0.7% in chimera 8. However, the substrate selectivity remained unchanged at any fatty acid concentrations ( Table 3 ).  Table 1 . Substrate conversion effi ciency of each chimera was classifi ed into three groups: ND, not detected; +, 1-10%; ++, 11-30%; +++, 31-100%. C: The relative substrate conversion rate of each recombinant protein expressed in S. cerevisiae , determined by adding 0.25 mM LA and 0.25 mM ALA after induction with 2% galactose. Relative substrate conversion rate = 100 × [product/(product + substrate)]/Density WB /Density SDS-PAGE . WB, Western blot. Histograms were established as the means of three independent samples (three transformants in SC-U plates) and error bars represent standard deviations. The + signs above the bars represent relative conversion rate determined from this study, and are as shown in (B). Because all known membrane-bound fatty acid desaturases contain the characteristic three His-box motifs, sequences between His boxes I and II may be a key the single or double mutants increased LA conversion rate, which suggests that the substrate preference of MpFADS6 may be dependent on several aa residues. The substrate preference can be further analyzed in terms of the molecular structure of cis -LA and cis -ALA. Figure 6A illustrates the C=C double-bond structure in the -3 position of cis -ALA and the C-C bond in the same position of cis -LA. We reasoned that the molecular mechanism of substrate preference for MaFADS6 and MpFADS6 may be their binding with the corresponding position of the respective substrate. It is likely that the catalytic areas of two enzymes are in common because, except for the -3 position substrate, overall molecular structure is identical. We propose the following model ( Fig. 6B ) based on the structures of stearoyl-CoA desaturase reported recently ( 32,33 ): there are four membrane-spanning helices (1, 2, 9, and 10) that are connected through two short endoplasmic reticulum lumen loops, and two membrane-embedding helices ( 3,13 ). Key aas between His boxes I and II and the catalytic center of FADS6, comprised by three His boxes and two ferric ions, should interact with and desaturate the sixth and seventh carbon positions.

CONCLUSIONS
In this study, MaFADS6 and MpFADS6 were confi rmed to have a preference for LA and ALA substrates, respectively, when expressed in S. cerevisiae . Reciprocal domain swapping was used to seek critical fragments in MaFADS6 and MpFADS6 that are linked to substrate preference. Site-directed mutagenesis was then used to provide further understanding of the specifi city determinants. Our study presents a step forward to gain an in-depth understanding of the relationship between the structure and function of FADS6s. area interacting with fatty acids and dictating substrate preference.
Domain swapping between MpFADS6 and MaFADS6 reduces overall desaturase activity toward LA and ALA, suggesting fusion proteins may alter their tertiary structure. The relative conversion rate of chimera 7 for LA was higher than the wild-type MpFADS6, suggesting that section 2 of MaFADS6 plays a part in catalysis of LA. Based on this result, it was rationalized that section 2 of MpFADS6 plays a part in catalysis of ALA and, thus, the conversion rate of chimera 2 for ALA should be higher than the wild-type MaFADS6. In fact, we did not detect any SDA product by chimera 2. It is probably because the three-dimensional structure of chimera 2 was disrupted, which caused its inactivation. Therefore the section 2s of MaFADS6 and MpFADS6 may not be fully exchangeable. Similarly, chimera 9 retained some activity toward LA and ALA, whereas chimera 4 lost its activity. It appears that MpFADS6 is more accommodating in accepting swapped sequences in section 2 and section 4. On the contrary, MaFADS6 is less stringent in section 5, swapping according to the activity results chimeras 5 and 10 ( Fig. 3 ). Sections 4 and 5 from chimeras 4 and 10 are in common and no product was detected for these mutants. Our interpretation of this result is that the combination of several key sites at the junction of these two sections may disrupt their function. Therefore, a double-mutant was created by substituting residues 399V and 400I near the junction with I (isoleucine) and E (glutamic acid) (from the corresponding sites of MaFADS6). This decreased ALA conversion approximately by half, suggesting that these residues were not very relevant for enzymatic activity or were located far from the substrate binding site.  Supplemental Material can be found at: