|
|
||||||||
|
First published online August 19, 2005; 10.1104/pp.105.064451 Plant Physiology 139:79-87 (2005) © 2005 American Society of Plant Biologists Two FAD3 Desaturase Genes Control the Level of Linolenic Acid in Flax SeedBioriginal Food and Science Corporation, Saskatoon, Saskatchewan, Canada S7N 0W9 (P.V., Z.H., M.-A.M., X.Q.); and Crop Development Centre, University of Saskatchewan, Saskatoon, Saskatchewan, Canada S7N 0W9 (G.R.)
Industrial oil flax (Linum usitatissimum) and edible oil or solin flax differ markedly in seed linolenic acid levels. Despite the economic importance of low-linolenic-acid or solin flax, the molecular mechanism underlying this trait has not been established. Two independently inherited genes control the low-linolenic-acid trait in flax. Here, we identified two genes, LuFAD3A and LuFAD3B that encode microsomal desaturases capable of desaturating linoleic acid. The deduced proteins encoded by these genes shared 95.4% identity. In the low-linolenic-acid line solin 593-708, both LuFAD3A and LuFAD3B carry point mutations that produce premature stop codons. Expression of these genes in yeast (Saccharomyces cerevisiae) demonstrated that, while the wild-type proteins were capable of desaturating linoleic acid, the truncated proteins were inactive. Furthermore, the low-linolenic-acid phenotype in flax was complemented by transformation with a wild-type gene. Codominant DNA markers were developed to distinguish between null and wild-type alleles of both genes, and linolenic acid levels cosegregated with genotypes, providing further proof that LuFAD3A and LuFAD3B are the major genes controlling linolenic acid levels in flax. The level of LuFAD3 transcripts in seeds peaked at about 20 d after flowering, and transcripts were not detectable in leaf, root, or stem tissue. A dramatic reduction in transcript levels of both genes occurred in the low-linolenic-acid solin line, which was likely due to nonsense-mediated decay.
The seed oil of flax (Linum usitatissimum) is notable for its high level of linolenic acid, generally around 45% to 65%, which gives it a high drying quality, making it useful for industrial purposes. The development of low-linolenic-acid flax lines (Green, 1986
Linolenic acid is produced through the desaturation of linoleic acid by omega-3/delta-15 desaturases. In plants, this reaction occurs both in the plastids and in the endoplasmic reticulum (ER). Arabidopsis (Arabidopsis thaliana) carries two plastidial desaturases, FAD7 and FAD8 (Browse et al., 1986
Flax Has Two FAD3 Genes, LuFAD3A and LuFAD3B As a first step in obtaining FAD3 cDNAs from flax, degenerate primers that targeted the first and third His-rich motifs of FAD3-like proteins were designed. When these primers were used to amplify first-strand cDNA from developing Normandy (wild-type) seeds, a PCR fragment with high similarity to FAD3/omega-3 desaturase sequences from various plant species was produced. Screening of a Normandy cDNA library with this fragment led to the identification of a full-length cDNA clone, 1,475 bp in length and containing an open reading frame encoding a 392-amino acid protein, LuFAD3A (Fig. 1).
Several primer sets designed to amplify the intron regions of the LuFAD3A gene produced two distinct bands with Normandy genomic DNA. The sequences of these pairs of bands were very similar but not identical, and only one of the two bands had exon sequences matching those of LuFAD3A. This suggested that the additional band represented a second LuFAD3 gene. Primers specific for this second gene were used to amplify overlapping 3' and 5' regions from solin line 593-708 adaptor-ligated cDNA, and primers were designed from these sequences to amplify the full-length clones from Normandy and solin line 593-708 cDNA. The deduced amino acid sequences of LuFAD3A and LuFAD3B from wild-type (Normandy) flax are 95.4% identical. LuFAD3A carries two additional amino acids in the N-terminal region and lacks one amino acid in the C-terminal region. Besides these differences, only 15 amino acids vary between the LuFAD3A and LuFAD3B proteins (Fig. 1). The LuFAD3A and LuFAD3B cDNA sequences are deposited in GenBank as accession numbers DQ116424 and DQ116425, respectively.
The LuFAD3A cDNA from the solin line was isolated using the same adaptor-ligated cDNA library as that used for the solin LuFAD3B cDNA. Overlapping fragments were isolated using primers specific to the LuFAD3A gene, these fragments were aligned, and a reverse transcription (RT)-PCR employing primers from each end of the alignment was used to obtain the full-length clone. When compared to the Normandy LuFAD3A sequence, the solin LuFAD3A cDNA sequence contained a single nucleotide substitution 874 nucleotides from the translational start, converting an Arg codon (CGA) to a stop codon (TGA; Fig. 2). Therefore, the protein encoded by solin LuFAD3A was truncated at the C terminus and consisted of only 291 amino acids, as opposed to the wild-type 392-amino acid Normandy LuFAD3A protein. Comparison of the solin-type and wild-type (Normandy) LuFAD3B sequences showed that solin LuFAD3B also carried a nonsense mutation. This substitution was 162 nucleotides from the start site, converting a Trp codon (TGG) to a stop codon (TGA). As a result, the solin LuFAD3B cDNA encoded a very short protein of only 53 amino acids, as opposed to the 391-amino acid protein encoded by the wild-type LuFAD3B gene. The solin-type and wild-type LuFAD3A and LuFAD3B cDNAs were identical except for these point mutations, although the poly-A tail was attached at slightly different positions in the solin- and wild-type LuFAD3B cDNAs. To verify the occurrence of these point mutations, the regions containing the mutations were amplified from solin line 593-708 genomic DNA. The sequences of the genomic PCR products were identical to the corresponding regions of the respective cDNAs (data not shown).
Alignments of the cDNAs with a partial LuFAD3A genomic clone from Normandy indicated that the LuFAD3A mutation occurred in the fifth of six exons, 100 nucleotides upstream of the exon-intron junction, whereas the mutation in the LuFAD3B gene occurred in the first exon.
A partial wild-type LuFAD3A clone isolated from a Normandy genomic library included the region from the translation start site to approximately 1,000 bp upstream of this site. This region was capable of directing seed-specific expression in flax (Fig. 3). Histochemical staining of transgenic flax plants carrying this promoter region in front of the uidA (
When we attempted to amplify this promoter region from Normandy genomic DNA, we obtained two distinct bands. While the sequence of one of the PCR products matched exactly to that of the genomic clone, the other product was similar but not identical, and the alignment of the two products contained a number of gaps. The second product appeared to represent the LuFAD3B promoter. To confirm this, an internal sequence was reamplified using a primer from within the putative LuFAD3B promoter and a primer specific to the LuFAD3B coding region (Fig. 4). LuFAD3A and LuFAD3B upstream regions were also amplified from solin 593-708 genomic DNA, and, while the LuFAD3B upstream regions from each genotype were identical, the LuFAD3A promoter from the solin-type line carried a single substitution (C to T) 559 bp upstream from the translation start site (data not shown).
Wild-Type Flax FAD3 Genes Demonstrate Desaturase Activity when Expressed in Yeast To determine whether the wild-type (Normandy) cDNAs encoded proteins capable of delta-15 desaturase activity and whether the mutations carried by the solin-type cDNAs resulted in the production of inactive proteins, the LuFAD3A and LuFAD3B cDNAs from both genotypes were expressed in yeast (Saccharomyces cerevisiae; Fig. 5). Supplement of yeast cultures carrying either Normandy LuFAD3A or LuFAD3B with linoleic acid led to the production of linolenic acid, while desaturation of linoleic acid did not occur in yeast cultures carrying either of the solin 593-708 FAD3 cDNAs (Fig. 5) or in cultures carrying a vector control (data not shown). The average conversion rate (weight percent product/(weight percent substrate plus product) x100) in three cultures expressing the wild-type LuFAD3B was 7.79% ± 0.32% SD, which was slightly higher than that seen in cultures expressing LuFAD3A (5.75% ± 0.32% SD).
In control samples, endogenous 18:1n-9 was not desaturated. However, in samples carrying the Normandy LuFAD3 cDNAs without exogenous substrates, a small peak corresponding to the expected position of 18:2(9,15) was detected (data not shown), which likely resulted from the desaturation of endogenous 18:1n-9; this desaturation activity was also detected with a Brassica napus FAD3 gene (Reed et al., 2000
When cultures were fed
In untransformed solin-type flax, the level of linolenic acid is less than 2% of total seed fatty acids. However, in solin 593-708 flax transformed with the wild-type (Normandy) LuFAD3A gene under the control of the flax 2S storage protein conlinin promoter (Truksa et al., 2003
Although expression in yeast is a useful test for enzymatic activity, occasionally the activity of an enzyme when expressed in yeast will not exactly mirror its activity in plants. For example, Heilmann et al. (2004)
To confirm that the two genes we identified were the major genes responsible for controlling the level of desaturation of linoleic acid in flax, we designed markers to distinguish plant genotypes, allowing us to correlate genotypes with phenotypes. A comparison of the sequences of the wild-type and mutated genes indicated that the point mutation in the LuFAD3A gene resulted in the loss of a PvuI site, whereas a BsaJI site was lost due to the point mutation in the LuFAD3B gene. Primers were designed for each gene such that the amplification product would contain the mutation, and digestion of the product with the appropriate enzyme would produce fragments with distinguishable sizes. Using these codominant markers, it was possible to classify plants as homozygous mutant, homozygous wild type, or heterozygous at each of the two LuFAD3 loci (Fig. 6).
These markers were used to analyze progeny from a number of F7 plants originating from a cross between a wild-type flax line and a solin line (Fig. 7). Since flax is strictly self-pollinated, it is likely that the F7 plants were homozygous at most loci. Three seeds from each plant were analyzed by cutting seeds in half and performing gas chromatography analysis on one half-seed while germinating the other half and subjecting the seedling to marker analysis. The markers for LuFAD3A and LuFAD3B identified several lines that were homozygous mutant for one gene and homozygous wild type for the second gene (Fig. 7). Gas chromatography analysis of the seed from which these plants originated confirmed that the LuFAD3 markers segregated with plant phenotype. Plants with less than 2% of the total fatty acids in seed as linolenic acid carried markers for both mutated genes, and plants with linolenic acid contents above 50% carried markers for both wild-type genes (Table I). Plants carrying a homozygous mutation at a single LuFAD3 locus had intermediate levels of linolenic acid, and the lack of active LuFAD3B seemed to have a greater effect on linolenic acid content than did the lack of LuFAD3A. While the ratio of 18:2n-6/18:3n-3 was less than one in null LuFAD3A/wild-type LuFAD3B seeds, it was greater than one in null LuFAD3B/wild-type LuFAD3A seeds (Table I).
Transcript Levels of LuFAD3A and LuFAD3B Are Dramatically Reduced in Low-Linolenic-Acid Lines
Northern-blot analysis of RNA from various tissues of Normandy and solin 593-708 was performed using a probe amplified from the coding region of the LuFAD3A gene (Fig. 8). LuFAD3A and LuFAD3B share 95% identity at the nucleotide level in the coding region; therefore, this probe should hybridize well to both mRNA species. LuFAD3 expression in wild-type seeds peaked at around 20 d after flowering (DAF; Fig. 8A). While flax seed reaches physiological maturity at about 60 DAF, the exponential phase of linolenic acid accumulation begins at about 20 DAF, and by 32 DAF one-half of the total seed lipids have been synthesized (Stymne et al., 1992
Flax is unusual among crop plants in its very high level of seed linolenic acid. This appears to be due at least in part to the presence in this species of two genes encoding very similar enzymes capable of desaturating linoleic acid. These genes have additive effects, and lines carrying a mutation in a single gene have intermediate levels of linolenic acid (Table I). Lesions were identified in the coding regions of both of these genes in the low-linolenic-acid line solin 593-708. This low-linolenic-acid line was selected following EMS mutation of a wild-type line (Rowland, 1991
Even when both LuFAD3A and LuFAD3B were mutated, there was still a low level of linolenic acid (<2%) present in seeds. This suggests that one of the FAD3 genes may have some residual activity. However, although the mutation in the LuFAD3A gene is closer to the C-terminal end of the gene, we did not detect desaturase activity when this gene was expressed in yeast (Fig. 5). Therefore, it appears more likely that the activity of another enzyme, probably a plastidial desaturase, results in the accumulation of low levels of linolenic acid in the solin line. A second low-linolenic-acid line obtained independently by EMS, which also carries two mutations, had linolenic acid levels similar to those found in solin lines (Tonnet and Green, 1987
Both LuFAD3A and LuFAD3B transcript levels were dramatically lower in the low-linolenic-acid line solin 593-708 compared to the wild-type Normandy line. The presence of a premature stop codon in a transcript has been associated with a reduction in transcript levels in plants, yeast, and animals. This RNA degradation functions as an mRNA quality control system, known as nonsense-mediated mRNA decay (or NMD; for review, see Schell et al., 2002
While the deduced protein of the null LuFAD3B gene is missing 338 amino acids, that of the null LuFAD3A gene lacks only 101 amino acids from its C terminus. However, the C-terminal regions of ER-localized desaturases contain ER retrieval motifs (McCartney et al., 2004 The two LuFAD3 genes have very similar amino acid sequences, and both were capable of desaturating linoleic acid in yeast (Figs. 1 and 4). Furthermore, expression of the wild-type LuFAD3A gene under the control of the conlinin promoter was able to rescue the mutant phenotype. The association of intermediate levels of linolenic acid with single homozygous mutations and low linolenic acid with double null plants (Fig. 7; Table I) provides further evidence that LuFAD3A and LuFAD3B are the major genes responsible for desaturation of linoleic acid in flax. The codominant markers described here are also useful for plant breeding purposes, particularly for the rapid introgression of the low-linolenic-acid trait into adapted genotypes. Developing similar markers for other genes involved in fatty acid biosynthesis and combining them with the markers described here may accelerate the production of "designer" flax oil.
Plant Materials F8 plants from a cross between solin line E1747 and the fiber flax (Linum usitatissimum) variety Viking were used for marker analysis. The linseed flax cultivars CDC Normandy and solin line 593-708 were used for all other experiments. Flax plants were grown in a growth chamber at temperatures of 22°C/17°C (day/night) under a 16-h photoperiod.
Total RNA was extracted from flax embryos at 15 to 20 DAF using an RNeasy Plant mini kit (Qiagen), and mRNA was isolated using the Dynabeads Oligo(dT)25 system (Dynal). The library was constructed with the ZAP cDNA synthesis kit and cloned into the UniZap XR vector (Stratagene). As a first step in isolating the LuFAD3A cDNA, the degenerate primer set D15F [5'-AT(ACGT) T(GT)(ACGT) GG(AG) AA(ACGT) A(GA)(GA) TG(AG) TG-3'] and D15R [5'-(AG) T(AGCT) GG(AGCT) CA(TC) GA(TC) TG(TC) GG(AGCT) CA-3'] was used to obtain a PCR product that had high homology to FAD3 genes. This fragment was then cloned into pCR 2.1 (Invitrogen), and the plasmid insert was used to screen the Normandy flaxseed cDNA library. The 3' end of the wild-type (Normandy) LuFAD3B gene was isolated from this cDNA library via PCR using a T7 primer and a LuFAD3B primer based on sequences from the solin-type LuFAD3B sequence (see below), while 5' sequences were obtained by RT-PCR using primers designed based on the solin LuFAD3B sequence.
Seed from solin line 593-708 was collected at 15 to 20 DAF, removed from the pod, frozen in liquid nitrogen, and ground to a fine powder. RNA was extracted according to the method described by Carpenter and Simon (1998)
DNA was isolated from 15-d-old flax cv Normandy seedlings using Qiagen Genomic Tips (Qiagen). After partial digestion and size fractionation, fragments between 16 to 21 kb were cloned into the Lambda Fix II vector (Stratagene). The LuFAD3A cDNA from Normandy flax was used as a probe to screen the library.
Precultures were grown at 30°C overnight in minimal media supplemented with 2% Glc and lacking uracil. Precultures were then washed and used to inoculate 20 mL of induction media (minimal media lacking uracil and supplemented with 2% Gal and 1% raffinose) to an OD600 of 0.5. Cultures were induced for 6 h at 20°C, then supplemented with 150 µM linoleic acid (18:2n-6), and grown at 20°C for 3 d.
For analysis of fatty acids from yeast (Saccharomyces cerevisiae), cells were pelleted by centrifugation, washed once with induction media and once with water, suspended in 2 mL of 3 N methanolic HCl, and incubated at 80°C for 30 min. After the addition of 1 mL of 0.9% NaCl, the sample was extracted with 2.5 mL of hexane. The hexane extract was dried and resuspended in 200 µL of hexane and analyzed on a Hewlett-Packard 5890A gas chromatograph equipped with a DB-23 column (30 mx250 µmx0.23 µm). The temperature was programmed to hold at 160°C for 1 min, increase to 240°C at 4°C/min, and hold for 10 min. For analysis of fatty acids from plants, single seeds were cut in half, and one or both pieces were incubated in 1 mL of 3 N methanolic HCl at 80°C for 30 min, after which the extraction essentially followed the procedure described above.
Hypocotyls of 5- to 6-d-old seedlings of flax were used as explants. Agrobacterium tumefaciens strain GV3101 containing the helper plasmid pMP90 was used as the host for a binary vector containing either the LuFAD3A cDNA under the control of the flax 2S storage protein conlinin promoter (Truksa et al., 2003
Tissues were stained overnight in staining solution (0.5 mg/mL 5-bromo-4-chloro-3-indolyl-beta-D-glucuronide, 50 mM Na phosphate buffer, pH 7, 1 mM K3Fe(CN)6, 1 mM K4Fe(CN)6, 20 mM EDTA), and fixed after staining in 5% v/v formaldehyde, 5% v/v acetic acid, 20% v/v ethanol.
Genomic DNA was isolated from flax leaf discs using the method described by Edwards et al. (1991) The markers were tested against 11 F8 lines from a recombinant inbred population from the cross E1747/Viking. Seeds from each line were cut in half and fatty acid analysis was performed on one half-seed, while the second half-seed was germinated and DNA was extracted and analyzed as described above. Fatty acid and PCR analyses were performed on three seeds from each line.
The gene-specific primer sets RTAF2 (5'-CCACCAACGGCAGCACC-3')/RTAR2 (5'-AATCCAGTAGAGAGGCCAGACAG-3') and RTBF2 (5'-CCACCACCAACGGCAATG-3')/RTBR2 (5'-ACCCAGTAAAGAGGCCAGAAAG-3') were designed to amplify approximately 250-bp regions from LuFAD3A and LuFAD3B, respectively. To ensure that the measured levels would be specific to a particular gene, a PCR experiment was conducted under the same conditions as those used for the real-time PCR, with plasmids containing the LuFAD3 cDNAs as templates. The LuFAD3A primers were unable to amplify the LuFAD3B cDNA and vice versa.
Total RNA was treated with FPLCpure DNAse I (Pharmacia), and cDNA was synthesized using the Superscript III first-strand synthesis system for RT-PCR (Invitrogen) according to the manufacturer's instructions. Real-time PCR was performed on an iCycler (Bio-Rad). The PCR was performed with Platinum SYBR Green qPCR SuperMix UDG (Invitrogen) according to the manufacturer's instructions. Primers for a conlinin seed storage protein gene were used for an endogenous control for normalization. Three separate first-strand cDNA reactions were analyzed in duplicate for each sample, and expression levels were calculated as described by Livak and Schmittgen (2001) Sequence data from this article can be found in the GenBank/EMBL data libraries under accession numbers DQ116424 and DQ116425.
We thank Dr. Nagamani Datla for performing flax transformation experiments and Dr. Stephen Knapp for providing the recombinant inbred population of E1747/Viking. Received April 20, 2005; returned for revision July 8, 2005; accepted July 11, 2005.
1 Present address: Lineberger Comprehensive Cancer Center, CB# 7295, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599. Article, publication date, and citation information can be found at www.plantphysiol.org/cgi/doi/10.1104/pp.105.064451. * Corresponding author; e-mail xqiu{at}bioriginal.com; fax 3069754839.
Anderson P (1995) Mutagenesis. Methods Cell Biol 48: 3158[Web of Science][Medline] Baker KE, Parker R (2004) Nonsense-mediated mRNA decay: terminating erroneous gene expression. Curr Opin Cell Biol 16: 293299[CrossRef][Web of Science][Medline] Bhatty RS (1995) Nutrient composition of whole flaxseed and flaxseed meal. In SC Cunnane, LU Thompson, eds, Flaxseed in Human Nutrition. AOCS Press, Champaign, IL, pp 2242
Browse J, McConn M, James D Jr, Miquel M (1993) Mutants of Arabidopsis deficient in the synthesis of
Browse JA, McCourt PJ, Somerville CR (1986) A mutant of Arabidopsis deficient in C18:3 and C16:3 leaf lipids. Plant Physiol 81: 859864 Carpenter CD, Simon AE (1998) Preparation of RNA. In JM Martinez-Zapater, J Salinas, eds, Methods in Molecular Biology, Vol 82. Arabidopsis Protocols. Humana Press, Totowa, NJ, pp 8589
Edwards K, Johnstone C, Thompson C (1991) A simple and rapid method for the preparation of plant genomic DNA for PCR analysis. Nucleic Acids Res 19: 1349 Green AG (1986) A mutant genotype of flax (Linum usitatissimum L.) containing very low levels of linolenic acid in its seed oil. Can J Plant Sci 66: 499503 Green AG, Marshall DR (1984) Isolation of induced mutants in linseed (Linum usitatissimum) having reduced linolenic acid content. Euphytica 33: 321328
Heilmann I, Pidkowich MS, Girke T, Shanklin J (2004) Switching desaturase enzyme specificity by alternate subcellular targeting. Proc Natl Acad Sci USA 101: 1026610271 James DW Jr, Dooner HK (1990) Isolation of EMS-induced mutants in Arabidopsis altered in seed fatty acid composition. Theor Appl Genet 80: 241245 Lemieux B, Miquel M, Somerville C, Browse J (1990) Mutants of Arabidopsis with alterations in seed lipid fatty acid composition. Theor Appl Genet 80: 234240[Web of Science]
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2 McCartney AW, Dyer JM, Dhanoa PK, Kim PK, Andrews DW, McNew JA, Mullen RT (2004) Membrane-bound fatty acid desaturases are inserted co-translationally into the ER and contain different ER retrieval motifs at their carboxy termini. Plant J 37: 156173[CrossRef][Web of Science][Medline]
McConn M, Hugly S, Browse J, Somerville C (1994) A mutation at the fad8 locus of Arabidopsis identifies a second chloroplast
Pendle AD, Clark GP, Boon R, Lewandowska D, Lam YW, Andersen J, Mann M, Lamond AI, Brown JWS, Shaw PJ (2005) Proteomic analysis of the Arabidopsis nucleolus suggests novel nucleolar functions. Mol Biol Cell 16: 260269
Reed EW, Schafer UA, Covello PS (2000) Characterization of the Brassica napus extraplastidial linoleate desaturase by expression in Saccharomyces cerevisiae. Plant Physiol 122: 715720 Rowland GG (1991) An EMS-induced low-linolenic-acid mutant in McGregor flax (Linum usitatissimum L.). Can J Plant Sci 71: 393396 Schell T, Kulozik AE, Hentze MW (2002) Integration of splicing, transport and translation to achieve mRNA quality control by the nonsense-mediated decay pathway. Genome Biol 3: REVIEWS1006 Stymne S, Tonnet ML, Green AG (1992) Biosynthesis of linolenate in developing embryos and cell-free preparations of high-linolenate linseed (Linum usitatissimum) and low-linolenate mutants. Arch Biochem Biophys 294: 557563[Medline] Tonnet ML, Green AG (1987) Characterization of the seed and leaf lipids of high and low linolenic acid flax genotypes. Arch Biochem Biophys 252: 646654[Medline] Truksa M, MacKenzie SL, Qiu X (2003) Molecular analysis of flax 2S storage protein conlinin and seed specific activity of its promoter. Plant Physiol Biochem 41: 141147[CrossRef]
Wagner E, Lykke-Andersen J (2002) mRNA surveillance: The perfect persist. J Cell Sci 115: 30333038 This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ASPB Publications | PLANT PHYSIOLOGY® | THE PLANT CELL | |
|---|---|---|---|