- Research article
Development and characterization of low α-linolenic acid Brassica oleracea lines bearing a novel mutation in a ‘class a’ FATTY ACID DESATURASE 3gene
BMC Geneticsvolume 15, Article number: 94 (2014)
Traditional canola (Brassica napus L.; AACC, 2n = 38) cultivars yield seed oil with a relatively high proportion of α-linolenic acid (ALA; C18:3cisΔ9,12,15), which is desirable from a health perspective. Unfortunately, due to the instability of this fatty acid, elevated levels also result in oils that exhibit a short shelf life and problems associated with use at high temperatures. As a result, the development of cultivars bearing reduced amounts of ALA in their seeds is becoming a priority. To date, several low ALA B. napus cultivars (~2-3% ALA of total fatty acids) have been developed and molecular analyses have revealed that the low ALA phenotype of lines tested thus far is a result of mutations within two ‘class b’ FATTY ACID DESATURASE 3 (FAD3) genes. Since B. napus possesses six FAD3 genes (two ‘class a’, two ‘class b’ and two ‘class c’) and ALA levels of approximately 2-3% remain in these low ALA lines, it is likely that the mutation of additional FAD3 genes could further decrease the content of this fatty acid.
In this study, we generated low ALA (≤2%) lines of B. oleracea, which is the C genome progenitor species of B. napus, via ethyl methanesulphonate (EMS) mutagenesis. We identified a novel nonsense mutation within the ‘class a’ FAD3 gene (BoFAD3-2) in these lines, which would result in the production of an encoded protein lacking 110 amino acids at its C terminus. When expressed in Saccharomyces cerevisiae, this mutant protein exhibited a drastic decline in its Δ-15 desaturase activity compared to the wild-type (wt) protein. Furthermore, we demonstrated that the expression of the mutant BoFAD3-2 gene was significantly reduced in developing seeds of low ALA lines when compared to expression in wt plants.
Given the additive nature of FAD3 mutations on ALA content and the ease with which B. napus can be re-synthesized from its progenitor species, the mutant isolated here has the potential to be used for the future development of B. napus cultivars exhibiting further reductions in ALA content.
Canola (Brassica napus L., AACC, 2n = 38) is an amphidiploid species bearing A and C genomes derived from the parental species B. rapa (A genome, n = 10) and B. oleracea (C genome, n = 9) ,, and is one of the most important oil crops worldwide. Typically, the seed oil of traditional canola cultivars is characterized by a low amount of saturated fatty acids (FAs; ~7%) and high contents of unsaturated FAs such as oleic acid (C18:1cisΔ9; ~55-60%), linoleic acid (LA, C18:2cisΔ9,12; ~20%) and α-linolenic acid (ALA, C18:3cisΔ9,12,15; ~10%). This relatively high level of ALA is desirable in terms of nutrition as it is an essential component of the human diet and is believed to provide various health benefits ,.
Indeed, it is the proportions of the three major unsaturated FAs (oleic acid, LA and ALA) that largely determine the usefulness of this oil from both food and industrial perspectives. While high levels of ALA may confer nutritional advantages, this same FA has been found to lead to rapid oxidation and instability during frying ,. It follows then that decreasing the level of ALA in canola seeds would be of value to enhance shelf life and stability of its oil at high temperatures [7-9]. To this effect, the development of low ALA (≤2-3%) cultivars has recently become one of the major breeding goals for the improvement of B. napus.
In oilseeds, ALA is mainly synthesized via the stepwise desaturation of oleic acid to LA, and LA to ALA, by endoplasmic reticulum (ER)-bound fatty acid desaturases (FADs), including the Δ-12 FAD2 and Δ-15 FAD3 enzymes, respectively ,. The B. napus genome contains 6 FAD3 genes, three of which are derived from the A genome and three from the C genome ,. These six genes (termed BnaA.FAD3.a, BnaA.FAD3.b, BnaA.FAD3.c, BnaC.FAD3.a, BnaC.FAD3.b and BnaC.FAD3.c) can be phylogenetically grouped into three classes (a, b and c), with each class comprising a single gene from the A and C genomes, respectively .
To date, several low ALA (approximately 2-3%) B. napus lines harbouring mutations within FAD3 genes have been characterized. Intriguingly, in every case, these mutations were apparent within only the ‘class b’ FAD3 genes, comprising two of the six FAD3 genes present in this species (BnaA.FAD3.b and BnaC.FAD3.b) ,,. Moreover, the mutations themselves were always identical, with BnaA.FAD3.b exhibiting a C to T substitution in the third position of the sixth codon in the seventh exon resulting in a missense mutation, and BnaC.FAD3.b bearing a G to A transition in the 5′ splice donor site of the sixth intron resulting in abnormal splicing. When combined, these two FAD3 mutations exhibited a cumulative effect on the reduction of ALA content , which insinuates that mutations in additional B. napus FAD3 genes will be key for obtaining further reductions in ALA within seed oil.
The alteration of ALA content to extremely low levels in B. napus through direct mutagenesis of this species would be a very complex endeavor as the trait is controlled by at least six gene loci (three from each of the A and C genomes). Since this would necessitate the use of an enormously large mutagenized population, we instead chose to carry out this study on one of its parental species for ease of manipulation due to the fact that it carries only half of the gene loci. Lines of the parental species exhibiting the desired mutant gene(s) can subsequently be utilized in B. napus breeding either through the re-synthesis of this species  or interspecific crossing with B. napus,.
We previously isolated two FAD3 transcripts from the seeds of B. oleracea and characterized low ALA lines bearing a missense mutation in the BoFAD3-1 coding region (corresponding to B. napus BnaC.FAD3.b) that resulted in a significant reduction in the desaturase activity of its encoded enzyme . While the mutation itself yielded a novel allele, it occurred in the gene corresponding to the previously identified ‘class b’ C genome FAD3 gene in which a mutation at a different position was identified in low ALA lines of B. napus,,. In this study, we isolated low ALA (approximately 2%) B. oleracea lines derived from ethyl methanesulphonate (EMS) mutagenesis and identified a nonsense mutation within the BoFAD3-2 gene (corresponding to B. napus BnaC.FAD3.a) that drastically reduced its enzymatic activity. To the best of our knowledge, this is the first instance in which a mutation in a Brassica ‘class a’ FAD3 gene has been linked to a low ALA phenotype. This novel mutant has the potential to be of use in the breeding of B. napus to produce seed oil bearing less than 2-3% ALA through its combination with previously identified mutant FAD3 alleles.
EMS treatment and generation of Mplants
A total of 600 seeds were treated with 0.5% EMS. While treated seeds only demonstrated a 3% reduction in viability when compared to wt B. oleracea, EMS treatment had a significant effect on M1 plants (Table 1). Indeed, approximately 30% of M1 plants failed to produce seed, with M1 plants producing an average of 29 seeds/plant and wt B. oleracea producing approximately 90 seeds/plant.
Selection of a low ALA phenotype in M2 to M7 generations
Of a total of 4,759 M2 plants transplanted to field plots, 97.5% produced seeds under open-pollination. Mature seeds harvested from 2,127 M2 plants (i.e. M3 seeds) were used for FA analysis. The proportion of ALA in this population varied between 3.89-15.92% (mean 11.48% ± 0.03 SE) of total FA content, while the seeds of wt B. oleracea contained 8.44-11.55% ALA (mean 10.15% ± 0.14 SE) (Figure 1, Additional file 1: Table S1). Of the 2,127 M3 seed families, 304 (14.3%) had a significantly lower ALA content than wt B. oleracea (confidence limit, CL0.05 = 9.78-10.52%). M3 generation plants were grown from 9 M3 seed families each bearing less than 7% ALA. The ALA content in seeds harvested from these M3 plants ranged from 2.28 to 11.13% (mean 5.55% ± 0.25 SE), which was significantly lower than wt B. oleracea (Figure 2, Additional file 1: Table S1). Further selection for low ALA content resulted in genetically stable low ALA mutant M7 lines with ≤ 2% ALA (Figure 2, Additional files 1 and 2).
When compared to wt B. oleracea grown under the same conditions, M7 lines not only exhibited a 6.5% reduction in ALA, but also a similar increase (6.9%) in the content of oleic acid. Conversely, the levels of LA in these lines were comparable to those in wt B. oleracea (Additional file 2: Table S2). No significant differences were noted between wt and mutant lines with respect to other FAs, such as C12:0, C14:0, C16:0, C18:0, C20:0, C20:1, C22:0, C22:1, C24:0, and C24:1 (Additional file 2: Table S2).
Low ALA mutant lines possess a single nonsense point mutation in the BoFAD3-2 coding region and reduced transcript levels compared to wt
In order to determine whether the low ALA phenotype in mutant lines was associated with one or more mutations within either the BoFAD3-1 or BoFAD3-2 coding regions, wt sequences were compared with those isolated from two low ALA EMS mutant lines. Interestingly, while the BoFAD3-1 sequences were identical in wt and mutant lines, BoFAD3-2 contained a single nucleotide substitution from G to A (Figure 3A) at position +822 (where +1 corresponds to the first nucleotide of the translational start codon) in both mutant lines compared to wt. This site lies within the putative seventh of eight exons of the BoFAD3-2 gene and results in the conversion of a TGG codon (encoding tryptophan) to a TGA stop codon. This mutation would correspond to premature translational termination and a deduced truncated protein product of 273 amino acids versus the deduced wt BoFAD3-2 protein of 383 amino acids (Figure 3B).
Semi-quantitative RT-PCR was carried out to determine whether any differences in expression levels or splice variants of BoFAD3-1 and/or BoFAD3-2 existed within the mutant lines compared to wt. While no variations in splicing were noted between wt and low ALA mutant lines in either gene, both low ALA mutant lines tested appeared to possess lower levels of BoFAD3-2, but not BoFAD3-1, transcripts than wt (Figure 4A). Subsequent quantitative real-time RT-PCR assays confirmed this apparent reduction in BoFAD3-2 transcript levels in low ALA mutant lines compared to wt (Figure 4B).
The mutant BoFAD3-2 protein exhibits reduced desaturase activity in yeast compared to wt
To establish whether the mutant BoFAD3-2 coding region from low ALA lines encoded a protein with reduced Δ-15 desaturase activity compared to the wt sequence, both wt and mutant BoFAD3-2 coding regions were expressed in S. cerevisiae. Supplementation of the transformed yeast cultures bearing wt BoFAD3-2 with LA led to the production of ALA at an average level of 5.76% ± 0.98 SD of the total FA content. In contrast, yeast containing the mutant BoFAD3-2 sequence produced only extremely low levels of ALA when supplemented with LA (0.00267% ± 0.00032 SD) while yeast bearing empty vector produced no ALA in any case (Figure 5A). The average conversion rate in LA-supplemented cultures expressing wt BoFAD3-2 was 52.25% ± 7.36 SD, compared to a 0.025% ± 0.003 SD conversion rate in cultures expressing the mutant BoFAD3-2 (Figure 6). Neither yeast transformed with empty vector, nor transformed yeast expressing either of the two BoFAD3 variants that were not supplied with exogenous LA, produced any detectable ALA (Figure 5B). Furthermore, yeast fed with exogenous LA and expressing the mutant BoFAD3-2 coding region contained significantly increased levels of LA (10.81% ± 0.31 SD) than yeast bearing the wt BoFAD3-2 sequence (5.23% ± 0.67 SD).
Interestingly, although cultures that were not supplemented with exogenous LA did not produce any detectable ALA, yeast containing the wt BoFAD3-2 not supplemented with LA generated significantly higher levels of endogenous C18:2 (0.85% ± 0.16 SD of total FA content) than non-fed cultures bearing either the mutant BoFAD3-2 or empty vector (0.04% ± 0.04 SD and 0.02% ± 0.03 SD of total FA content, respectively), which were not significantly different from one another (Figure 5B). This latter result was likely due to the Δ-15 desaturation of endogenous oleic acid (C18:1cisΔ9) to C18:2 cisΔ9,15 by the wt BoFAD3-2 enzyme, as has been reported to occur previously when both flax (Linum usitatissimum) and B. napus FAD3 enzymes were expressed in S. cerevisiae,.
As a means to enhance the stability and shelf life of canola oil, the development of B. napus cultivars bearing low levels of ALA in their seed oil is an important consideration for breeders. At present, several low ALA (approximately 2-3%) cultivars are available that are the result of mutations within FAD3 genes ,[22-24]. However, of the six FAD3 genes present within the B. napus genome, low ALA phenotypes have only been linked to mutations within the two ‘class b’ FAD3 genes as of yet (BnaA.FAD3.b and BnaC.FAD3.b) ,,. Since FAD3 mutations appear to have an additive effect on the reduction of ALA content in various plant species ,,, the isolation of low ALA lines resulting from mutations in additional classes of Brassica FAD3 genes would be valuable for the future production of lines bearing further reductions in this FA. In this study, we generated and analyzed low ALA (≤2%) EMS mutant lines of B. oleracea, which is one of the progenitor species of B. napus (C genome), and identified a nonsense mutation within the ‘class a’ BoFAD3-2 gene (corresponding to BnaC.FAD3.a). This is the first instance in which a mutation within a ‘class a’ Brassica FAD3 gene  has been linked to low ALA content.
Low ALA content in the seeds of mutant B. napus lines has previously been attributed to a C to T substitution resulting in a missense mutation in BnaA.FAD3.b and a G to A transition in the 5′ splice donor site of the sixth intron resulting in an impairment in splicing in BnaC.FAD3.b,,. Additionally, low ALA mutant B. oleracea lines have been generated that bore a single missense mutation within the putative third exon of the BoFAD3-1 gene (corresponding to BnaC.FAD3.b), resulting in a significant reduction in its desaturase activity . While no nonsense mutations have yet been found in low ALA Brassica mutants prior to this study, FAD3 genes exhibiting this type of mutation have been implicated in low ALA phenotypes in other plant species, such as soybean (Glycine max) [26-28] and flax ,.
In addition to the identified BoFAD3-2 mutation, analyses of the expression levels of BoFAD3-1 and BoFAD3-2 in the mutant low ALA B. oleracea lines indicated that while BoFAD3-1 transcript levels were apparently unchanged in mutant lines compared to wt, BoFAD3-2 transcripts were significantly reduced in the mutants compared to wt (Figure 4). A similar reduction in expression was noted previously in a mutant FAD3 gene isolated from a low ALA line of flax , and was suggested to be a consequence of nonsense-mediated mRNA decay (NMD), which serves as an mRNA quality control system that identifies and degrades mRNAs containing premature stop codons in eukaryotic organisms . In plants, both unusually long 3′ UTRs and an intron located at least 50-nt downstream of the stop codon often induce NMD ,. This is almost certainly also the case in this study, as the BoFAD3-2 mutation occurs in the seventh of eight exons, with the putative exon-exon junction occurring over 130-nt downstream of the mutation. However, an additional mutation within the BoFAD3-2 promoter region resulting in down-regulation of the gene cannot be ruled out at this point.
While one would expect a nonsense mutation occurring at a reasonable distance from the 3′ end of a gene to produce a truncated, and theoretically inactive, version of the encoded protein, it appears that the nonsense mutation within the BoFAD3-2 gene identified in this study may not have completely abolished the desaturation activity of this enzyme (Figures 5 and 6). Although the levels of ALA generated in yeast expressing a mutant copy of the BoFAD3-2 coding region were miniscule compared to yeast bearing a wt copy of the same gene (0.0027% ± 0.00032 SD vs. 5.76% ± 0.98 SD of the total FA content, respectively), yeast containing the empty vector control did not produce any ALA whatsoever. This suggests that the mutant BoFAD3-2 protein retained at least a very small portion of residual activity. FAD3 proteins are known to contain three highly conserved histidine boxes involved in the catalytic site of the enzyme that are essential for their function , as well as a C-terminal dilysine motif (KSKIN) required for localization to the ER and their enzymatic activity . The deduced mutant BoFAD3-2 protein identified here would lack 110 amino acids of its C terminal end (of 383 amino acids total; Figure 3B), including one of the three histidine boxes and the ER-retrieval motif, which indicates that a very small level of activity may be possible even in the absence of these two domains. This contrasts with a previous finding, whereby a flax FAD3 protein lacking only 101 amino acids of its C-terminal end was found to be completely inactive when expressed in S. cerevisiae. While it is possible that the flax and B. oleracea truncated FAD3 enzymes exhibited different activities, it is also feasible that the extremely low levels of activity observed here were simply below the limit of detection in the previous study.
Seeds from the low ALA B. oleracea mutant lines isolated in this study also retained at least some level of Δ-15 desaturase activity, as they generated approximately 2% ALA (Figure 2, Additional file 1: Table S1). While it is possible that the mutant BoFAD3-2 enzyme may have contributed a very small proportion of this ALA, it is likely that the vast majority was mainly attributable to the remaining wt BoFAD3 enzyme(s) present within these lines. While the B. oleracea genome contains three BoFAD3 genes, we were previously able to isolate only two of the corresponding transcripts (BoFAD3-1 and BoFAD3-2) from developing seeds , which implies that these may be the major Δ-15 desaturases within the seeds of this species. In line with this, both BnaC.FAD3.a (corresponding to BoFAD3-2) and BnaC.FAD3.b (corresponding to BoFAD3-1) transcripts have been amplified from B. napus previously; however, the same is not true for the third C genome FAD3 gene, BnaC.FAD3.c. It is also feasible that the activity of another enzyme, such as a plastidial desaturase, could be contributing to the low levels of ALA present in these lines. This appears to be the case in Arabidopsis lines bearing mutations in their single microsomal FAD3 gene, which still retain low levels of ALA in their seeds ,.
FAD3 mutations have not only been found to be associated with low ALA content in various plant species, but also a concomitant increase in LA content ,. Since LA acts as a substrate for FAD3, the inability of a mutant FAD3 gene to effectively convert LA to ALA results in the depletion of ALA and a buildup of LA. This phenomenon was noted in our yeast expression assays, where yeast supplemented with exogenous LA substrate and bearing the mutant BoFAD3-2 coding region contained significantly higher levels of LA (10.81% ± 0.31 SD) compared to yeast bearing the wt BoFAD3-2 sequence (5.23% ± 0.67 SD). Conversely, this was not the case in our low ALA B. olearacea mutant lines, which displayed no significant difference in LA content compared to wt lines (Additional file 1: Table S1). Instead, oleic acid content was significantly increased in the mutant plants. While certain single quantitative trait loci have previously been found to affect several FAs simultaneously in other plant species, such as B. rapa, soybean ,, and maize (Zea mays) , the fact that a similar increase in oleic acid was not apparent in our yeast expression assays (data not shown) indicate that a second uncharacterized mutation, possibly in a FAD2 gene, may also be present in our mutant B. oleracea lines.
In conclusion, we have isolated low ALA (≤2%) B. oleracea lines bearing a nonsense mutation in the BoFAD3-2 gene. This mutation resulted in a severe diminishment of its Δ-15 desaturase activity when expressed in S. cerevisiae and may also have elicited nonsense-mediated decay of the mRNA transcript within the mutant plants themselves. As of yet, this is the first mutation to be identified in a Brassica ‘class a’ FAD3 gene that has been linked to a low ALA phenotype. Due to the additive nature of FAD3 mutations on ALA content in plants, this novel mutant has the potential to be of great value in the future production of designer B. napus lines exhibiting further reductions in this FA than have been achieved previously.
A self-compatible variant of B. oleracea var. alboglabra (Chinese kale; hereafter referred to as B. oleracea) that exhibits a spring growth habit, and thus does not require vernalization for flowering, was used in this study (original seeds obtained from Lantmännen SW Seed, Sweden). EMS mutagenesis was carried out on a S6 generation inbred line developed through self-pollination of single plants. The ALA content in the seed oil of this particular line was approximately 9% (see Additional file 2: Table S2 for complete FA profile).
Seed mutagenesis was performed using EMS as described previously . In brief, seeds were immersed in 0.5% (v/v) EMS solution for 12 h at room temperature, and were subsequently rinsed and allowed to germinate on sterilized sand in a Petri dish. The viability of EMS-treated seeds was recorded 7-10 days after plating. M1 generation plants were transferred to soil (Sunshine Professional Growing Mix; Sunshine Horticulture, Bellevue, USA) and grown in a greenhouse (21°C/18°C ± 3°C day/night temperature and 16 h light with an intensity of 450 μMol M-2’s-1).
Generation of mutant population and selection of low ALA phenotype
M1 plants were found to produce a much smaller quantity of visible pollen than wild type (wt) B. oleracea. Therefore, in addition to being allowed to self-pollinate within isolation bags, all M1 plants were also self-pollinated manually in order to ensure they would produce a sufficient number of M2 seeds. Manual self-pollination was carried out by pollinating 10-20 newly opened flowers on each plant with pollen obtained from the same plant. M2 seeds were harvested at full maturity and were initially grown in a greenhouse (with conditions as described in the previous section) in 4°Cm ’ 4°Cm ’ 4.5°Cm (L ’ W ’ D) cells for approximately 45 days, after which time the seedlings were transplanted to the University of Alberta Research Farm in Edmonton, Alberta, Canada. This research field was isolated from other Brassica crops by a distance of approximately 300 M. Following transplantation, the field was irrigated once, and crop management practices recommended for Brassica species were followed.
Open-pollinated M3 seeds (i.e. seeds produced on M2 plants) were harvested at maturity on an individual plant basis and were utilized for FA analysis. Subsequent generations (M3 to M7) were grown either in a greenhouse (conditions as described in the previous section) or in a growth chamber (20°C/15°C day/night temperature and 16 h light with a photosynthetic photon flux density of 450 μMol M-2’s-1 at the plant level) and were self-pollinated using isolation bags. Single plants were selected for low ALA content in each generation. In every case, wt B. oleracea was grown alongside mutant plants as a control.
Isolation and sequence analysis of B. oleracea FAD3 (BoFAD3) cDNA
Total RNA was extracted from developing siliques 25-30 days after pollination (DAP) from wt B. oleracea, as well as two low ALA mutant B. oleracea lines, using the Qiagen Plant RNeasy kit according to the manufacturer’s instructions (Qiagen, Toronto, ON). Contaminating genomic DNA was eliminated using the TURBO DNA-free system (Ambion, Life Technologies Inc., Burlington, ON) and RNA concentrations were determined using a Nanodrop spectrophotometer (Nanodrop Products, Wilmington, DE). First-strand cDNA synthesis was conducted using the SuperScript’ III first-strand cDNA synthesis kit (Invitrogen, Life Technologies Inc.) with 1 μg total RNA and an oligo-dT primer. Subsequent PCR assays were carried out to amplify BoFAD3-1 [GenBank:JX866747] and BoFAD3-2 [GenBank:JX866748] transcripts, which had been found previously to be expressed in B. oleracea developing siliques , using Pfx50 high fidelity DNA polymerase (Invitrogen) along with primers BoFAD3-1F1Y and BoFAD3-1RY1  to amplify the full-length 1143-bp BoFAD3-1 coding region, and BoFAD3-2FY1 (5′ ’ GCA ATG GTT GTT GCT ATG TAC - 3′) and BoFAD3-2FY2 (5′ ’ AAG TTA ATT GAT TTT AGA TTT GTC - 3′) to amplify the full-length 1152-bp BoFAD3-2 coding region. PCR cycling parameters were as follows: 94°C for 5 Min, 33°Cycles of 94°C for 15’s, 54°C for 30 s, and 68°C for 2 Min, followed by a final extension at 68°C for 5 Min (for amplification of BoFAD3-1) or 94°C for 5 Min, 5°Cycles of 94°C for 15’s, 60°C for 30 s and 68°C for 2 Min, 6°Cycles of 94°C for 15’s, 57°C for 30 s and 68°C for 2 Min, 26°Cycles of 94°C for 15’s, 54°C for 30 s and 68°C for 2 Min, followed by a final extension of 68°C for 5 Min (for amplification of BoFAD3-2). The resulting amplification products were cloned into either the pCR 4-TOPO (Invitrogen) or pGEM-T easy vectors (Promega, Madison, WI) and sequenced. In the case of BoFAD3-1 sequences, 10 clones derived from wt plants and 7 clones derived from mutant plants (two separate lines) were sequenced. In the case of BoFAD3-2 sequences, 6 clones derived from wt plants and 17 clones derived from mutant plants (two separate lines) were sequenced. Sequences obtained were aligned with each other using CLUSTALW .
To ensure that the extreme 5′ and 3′ ends of the BoFAD3 coding regions did not exhibit any differences between low ALA and wt lines, 5′ and 3′ RACE were carried out as described previously . In each case, purified PCR products were cloned into the pGEM-T easy vector and sequenced, and full-length BoFAD3 coding sequences were obtained by combining the RT-PCR sequencing results with the RACE data.
Expression analysis of BoFADtranscripts in wt and low ALA mutant lines
To determine whether BoFAD3 transcripts from wt and low ALA mutant lines exhibited any variations in splicing or expression levels, semi-quantitative RT-PCR was carried out for both BoFAD3-1 and BoFAD3-2 transcripts as described previously  using 1 μL silique-derived cDNA (25-30 DAP) isolated from wt B. oleracea, as well as two low-ALA mutant lines, respectively, generated as described above. Primers BoFAD3F3c and BoFAD3R5  were utilized to amplify a 1108-nt BoFAD3-1-specific product, while BoFAD3F3b and BoFAD3R5  were utilized to amplify a 1110-nt BoFAD3-2-specific product. Additionally, a 677-nt BoPP2AA3 transcript [GenBank: DK472124] was amplified as an internal control using primers BoPP2AA3F1 and BoPP2AA3R1 . The conditions for PCR amplification of BoFAD3-specific products were 94°C for 2 Min, 28°Cycles of 94°C for 15’s, 55°C for 30 s, 68°C for 1.5 Min, with a final extension of 68°C for 7 Min. The same general parameters were used to amplify BoPP2AA3-specific fragments, with the exception of the annealing temperature, which was 58°C, and the extension time, which was 1 Min. PCR products were resolved on 1% agarose gels and visualized with SYBR Safe (Invitrogen).
To provide additional evidence of differences in BoFAD3-2 expression levels between wt and low ALA mutant lines, quantitative real-time RT-PCR was carried out on DNA-free total RNA isolated as described above from developing siliques (25-30 DAP) of two biological replicates of wt B. oleracea, as well as two separate low ALA lines. First-strand cDNA synthesis was conducted using 500 ng total RNA, along with Superscript III (Invitrogen) and an oligo-dT primer in a final volume of 10 μL. Quantitative PCR assays were performed in triplicate using 1 μL of a 1/20 dilution of each cDNA as template along with SYBR green PCR master mix in a final volume of 10 μL according to the manufacturer’s instructions (Applied Biosystems, Life Technologies Inc., Burlington, ON). Assays were carried out on an ABI 7900HT Fast Real-Time PCR System (Applied Biosystems) using primers BoFAD3F3b  and BoFAD3-2qPCRR7 (5′ ’ TCT CAA AGG ACT CTT CAC CCA G - 3′) to amplify a 115-nt fragment of the BoFAD3-2 transcript and GAPDH-LEFT and GAPDH-RIGHT  to amplify a 112-nt fragment of the constitutively expressed internal standard transcript, BoGAPDH [GenBank: EF123055]. Reactions lacking template cDNA and those without reverse transcriptase were included in each case as negative controls. Thermal parameters for amplification were 95°C for 2 Min, followed by 40°Cycles of 95°C for 15’s and 60°C for 1 Min. Dissociation curves were generated to ascertain that only a single product was produced in each assay. Relative levels of gene expression were attained using the standard curve method and SDS v2.4 software (Applied Biosystems), with all BoFAD3-2 expression data comprising mean values of the biological replicates normalized to those of BoGAPDH.
Heterologous expression of wt and mutant BoFAD3-2 in yeast
To determine the relative activities of wt and mutant BoFAD3-2 proteins, full-length coding regions were amplified from previously generated cDNA using the high-fidelity Pfx50 DNA polymerase along with primers BoFAD3-2F1Y and BoFAD3-2R1Y, which were designed to specifically amplify BoFAD3-2. The resulting PCR products were inserted downstream of the GAL1 promoter in the pYES2.1 TOPO®TA yeast expression vector (Invitrogen) according to the manufacturer’s instructions and were confirmed by sequencing in each case. Saccharomyces cerevisiae strain SCY62 (MAT a ADE2 can1-100 his3-11,15 leu2-3,112 trp1-1 ura3-1) was transformed with the resulting expression vectors containing the wt and mutant BoFAD3-2 coding regions, as well as empty vector, respectively, using the lithium acetate method .
To induce expression of the heterologous BoFAD3-2 sequences, three separate colonies bearing each vector, including the empty vector control, respectively, were initially utilized to produce yeast pre-cultures, which were incubated at 30°C for two days in minimal media supplemented with 2% raffinose and lacking uracil. These pre-cultures were then used to inoculate 2 ’ 80 ML of induction media (minimal media containing 2% galactose and 1% raffinose, but lacking uracil) to an OD600 of 0.2. In each case, one of the two cultures was supplied with 150 μM LA (in ethanol) along with 0.1% tyloxapol, while the other was supplemented with the same volume of ethanol and tyloxapol, without FA. The cultures were subsequently grown at 18°C with shaking at 225 rpm for 3 d.
FA analysis of B. oleracea seed oil was carried out via gas chromatography of FA methyl esters derived from mature bulk seed (0.1 ’ 0.25 g) harvested from M2 to M7 generation plants as described previously . Determination of the FA compositions of total lipids obtained from transformed yeast bearing wt BoFAD3-2, mutant BoFAD3-2, or empty vector, was accomplished using gas chromatography/mass spectrometry as described previously . In the case of yeast FA analyses, relative percentages of ALA were derived from peak areas and conversion rates of FA substrate to product were calculated by dividing the weight percent product by the sum of the weight percent substrate and product, then multiplying by 100.
Scattered diagrams and statistical analyses including mean, variance, standard error, confidence intervals and t-tests were calculated using the EXCEL program. For paired 2-tailed t-tests, differences were considered significant at P ≤ 0.05.
Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Jap J Bot. 1935, 7: 389-452.
Rana D, Boogaart T, O’Neill CM, Hynes L, Bent E, Macpherson L, Park JY, Lim YP, Bancroft I: Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives. Plant J. 2004, 40: 725-733. 10.1111/j.1365-313X.2004.02244.x.
Simopoulos AP: The importance of the omega-6/omega-3 fatty acid ratio in cardiovascular disease and other chronic diseases. Exp Biol Med. 2008, 233: 674-688. 10.3181/0711-MR-311.
Deshpande R, Mansara P, Suryavanshi S, Kaul-Ghanekar R: Alpha-linolenic acid regulates the growth of breast and cervical cancer cell lines through regulation of NO release and induction of lipid peroxidation. J Mol Biochem. 2013, 2: 6-17.
Barker GC, Larson RT, Graham IA, Lynn JR, Graham JK: Novel insights into seed fatty acid synthesis and modification pathways from genetic diversity and quantitative trait loci analysis of the Brassica C genome. Plant Physiol. 2007, 144: 1827-1842. 10.1104/pp.107.096172.
Scarth RT, Tang JH: Modification of oil using conventional and transgenic approaches. Crop Sci. 2006, 46: 1225-1236. 10.2135/cropsci2005.08-0245.
Warner K, Knowlton S: Frying quality and oxidative stability of high-oleic corn oils. J Am Oil Chem Soc. 1997, 74: 1317-1322. 10.1007/s11746-997-0063-7.
Matthäus B: Utilization of high-oleic rapeseed oil for deep-fat frying of French fries compared to other commonly used edible oils. Eur J Lipid Sci Technol. 2006, 108: 200-211. 10.1002/ejlt.200500249.
Merrill LI, Pike OA, Ogden LV, Dunn ML: Oxidative stability of conventional and high-oleic vegetable oils with added antioxidants. J Am Oil Chem Soc. 2008, 85: 771-776. 10.1007/s11746-008-1256-4.
Hu X, Sullivan-Gilbert M, Gupta M, Thompson SA: Mapping of the loci controlling oleic and linolenic acid contents and development of fad2 and fad3 allele-specific markers in canola (Brassica napus L.). Theor Appl Genet. 2006, 113: 497-507. 10.1007/s00122-006-0315-1.
Shanklin J, Cahoon EB: Desaturation and related modifications of fatty acids. Ann Rev Plant Physiol Plant Mol Biol. 1998, 49: 611-641. 10.1146/annurev.arplant.49.1.611.
Los DA, Murata N: Structure and expression of fatty acid desaturases. Biochimica Biophysica Acta. 1998, 1394: 3-15. 10.1016/S0005-2760(98)00091-5.
Smooker AM, Wells R, Morgan C, Beaudoin F, Cho K, Fraser F, Bancroft I: The identification and mapping of candidate genes and QTL involved in the fatty acid desaturation pathway in Brassica napus. Theor Appl Genet. 2011, 122: 1075-1090. 10.1007/s00122-010-1512-5.
Yang Q, Fan C, Guo Z, Qin J, Wu J, Li Q, Fu T, Zhou Y: Identification of FAD2 and FAD3 genes in Brassica napus genome and development of allele-specific markers for high oleic and low linolenic acid contents. Theor Appl Genet. 2012, 125: 715-729. 10.1007/s00122-012-1863-1.
Mikolajczyk K, Dabert M, Karlowski WM, Spasibionek S, Nowakowska J, Cegielska-Taras T, Bartkowiak-Broda I: Allele-specific SNP markers for the new linolenic mutant genotype of winter oilseed rape. Plant Breed. 2010, 129: 502-507.
Rahman MH: Fatty acid composition of resynthesized Brassica napus and trigenomic Brassica void of genes for erucic acid in their A genomes. Plant Breed. 2002, 121: 357-359. 10.1046/j.1439-0523.2002.00711.x.
Rahman MH, Bennett RA, Yang R-C, Thiagarajah MR: Exploitation of the late flowering species Brassica oleracea L. for the improvement of earliness in B. napus L. ’ an untraditional approach. Euphytica. 2011, 177: 365-374. 10.1007/s10681-010-0253-5.
Rahman H: Review: Breeding spring canola (Brassica napus L.) by the use of exotic germplasm. Can J Plant Sci. 2013, 93: 363-373. 10.4141/cjps2012-074.
Rahman H, Singer SD, Weselake RJ: Development of low-linolenic acid Brassica oleracea lines through seed mutagenesis and molecular characterization of mutants. Theor Appl Genet. 2013, 126: 1587-1598. 10.1007/s00122-013-2076-y.
Vrinten P, Hu Z, Munchinsky M-A, Rowland G, Qiu X: Two FAD3 desaturase genes control the level of linolenic acid in flax seed. Plant Physiol. 2005, 139: 79-87. 10.1104/pp.105.064451.
Reed DW, Schäfer UA, Covello PS: Characterization of the Brassica napus extraplastidial linoleate desaturase by expression in Saccharomyces cerevisiae. Plant Physiol. 2000, 122: 715-720. 10.1104/pp.122.3.715.
Röbbelen G, Nitsch A: Genetical and physiological investigations on mutants for polyenoic fatty acids in rapeseed, B. napus L. I. Selection and description of new mutants. Z Pflanzenz°Chtg. 1975, 75: 93-105.
Auld DL, Heikkinen MK, Erickson DA, Sernyk JL, Romero JE: Rapeseed mutants with reduced levels of polyunsaturated fatty acids and increased levels of oleic acid. Crop Sci. 1992, 32: 657-662. 10.2135/cropsci1992.0011183X003200030016x.
Spasibionek S: New mutants of winter rapeseed (Brassica napus L.) with changed fatty acid composition. Plant Breed. 2006, 125: 259-267. 10.1111/j.1439-0523.2006.01213.x.
Hu X, Sullivan-Gilbert M, Gupta M, Thompson SA: G-to-A mutation at the 5′ splice site of fad3c caused impaired splicing in a low linolenic mutant of canola (Brassica napus L.). Plant Biotech. 2007, 24: 397-400. 10.5511/plantbiotechnology.24.397.
Reinprecht Y, Luk-Labey S-Y, Larsen J, Poysa VW, Yu K, Rajcan I, Ablett GR, Pauls KP: Molecular basis of the low linolenic acid trait in soybean EMS mutant line RG10. Plant Breed. 2009, 128: 253-258. 10.1111/j.1439-0523.2008.01571.x.
Bilyeu K, Gillman JD, LeRoy AR: Novel FAD3 mutant allele combinations produce soybeans containing 1% linolenic acid in the seed oil. Crop Sci. 2011, 51: 259-264. 10.2135/cropsci2010.01.0044.
Chappell AS, Bilyeu KD: A GmFAD3A mutation in the low linolenic acid soybean mutant C1640. Plant Breed. 2006, 125: 535-536. 10.1111/j.1439-0523.2006.01271.x.
Thambugala D, Duguid S, Loewen E, Rowland G, Booker H, You FM, Cloutier S: Genetic variation of six desaturase genes in flax and their impact on fatty acid composition. Theor Appl Genet. 2013, 126: 2627-2641. 10.1007/s00122-013-2161-2.
Kert’sz S, Kerényi S, Mérai Z, Bartos I, Pélfy T, Barta E, Silhavy D: Both introns and long 3′-UTRs operate as cis-acting elements to trigger nonsense-mediated decay in plants. Nucleic Acids Res. 2006, 34: 6147-6157. 10.1093/nar/gkl737.
Nyikó T, Kerényi F, Szabadkai L, Benkovics AH, Major P, Sonkoly B, Mérai Z, Barta E, Niemiec E, Kufel J, Silhavy D: Plant nonsense-mediated mRNA decay is controlled by different autoregulatory circuits and can be induced by an EJC-like complex. Nucleic Acids Res. 2013, 41: 6715-6728. 10.1093/nar/gkt366.
McCartney AW, Dyer JM, Dhanoa PK, Kim PK, Andrews DW, McNew JA, Mullen RT: Membrane-bound fatty acid desaturases are inserted co-translationally into the ER and contain different ER retrieval motifs at their carboxy termini. Plant J. 2004, 37: 156-173. 10.1111/j.1365-313X.2004.01949.x.
James DW, Dooner JK: Isolation of EMS-induced mutants in Arabidopsis with alterations in seed lipid fatty acid composition. Theor Appl Genet. 1990, 80: 241-245. 10.1007/BF00224393.
Lemieux B, Miquel M, Somerville C, Browse J: Mutants of Arabidopsis with alterations in seed lipid fatty acid composition. Theor Appl Genet. 1990, 80: 234-240. 10.1007/BF00224392.
Cloutier S, Ragupathy R, Niu Z, Duguid S: SSR-based linkage map of flax (Linum usitatissimum L.) and mapping of QTLs underlying fatty acid composition traits. Mol Breeding. 2011, 28: 437-451. 10.1007/s11032-010-9494-1.
Tanhuanpaa PK, Vikki JP, Vikki JH: Mapping of a QTL for oleic acid concentration in spring turnip rape (Brassica rapa ssp oleifera). Theor Appl Genet. 1996, 92: 952-956. 10.1007/BF00224034.
Reinprecht Y, Poysa VW, Yu K, Rajcan I, Ablett GR, Pauls KP: Seed and agronomic QTL in low linolenic acid, lipoxygenase-free soybean (Glycine max L. Merrill) germplasm. Genome. 2006, 49: 1510-1527. 10.1139/g06-112.
Wassom JJ, Mikkelinem V, Bohn MO, Rocheford TR: QTL for fatty acid composition of maize kernel oil in Illinois high oil x B73 backcross-derived lines. Crop Sci. 2008, 48: 69-78. 10.2135/cropsci2007.04.0208.
Thurling M, Depittayanan V: EMS induction of early flowering mutants in spring rapeseed (Brassica napus). Plant Breed. 1992, 108: 177-184. 10.1111/j.1439-0523.1992.tb00118.x.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Broekgaarden C, Poelman EH, Steenhuis G, Voorrips RE, Dicke M, Vosman B: Genotypic variation in genome-wide transcription profiles induced by insect feeding: Brassica oleracea ’ Pieris rapae interactions. BMC Genomics. 2007, 8: 239-252. 10.1186/1471-2164-8-239.
Gietz RD, Woods RA: Transformation of yeast by lithium acetate/single-stranded carrier DNA/polyethylene glycol method. Methods Enzymol. 2002, 350: 87-96. 10.1016/S0076-6879(02)50957-5.
HR and RJW are grateful for the support provided by AVAC Ltd., the Canada Foundation for Innovation and the Research Capacity Program of Alberta Enterprise and Advanced Education. RJW is also grateful for the support provided by the Alberta Innovates Bio Solutions, Canada Research Chairs Program and the Natural Sciences and Engineering Research Council of Canada. The authors would also like to thank Dr. Mohan Thiagarajah for suggestions with regards to EMS treatment, Dr. Nidhi Sharma for collecting developing siliques, Ms. An Vo for FA analysis of the mutagenized poXpulations, and other laboratory staff for the technical assistance provided.
The authors declare that they have no competing interests.
SDS carried out molecular analyses, as well as enzyme assays, and wrote a portion of the manuscript. RJW investigated desaturase activity, contributed to the improvement of the manuscript, and secured funding for this research. HR developed the mutant lines and wrote a portion of the manuscript. All authors read and approved the final manuscript.