Abstract

To identify genetic variation for fatty acid composition in mature soybean seeds, 4566 M3 generation seed samples from a chemically mutagenized population were subjected to fatty acid profiling by gas chromatography. In the population, a wide range of variation in the content for each of the five major fatty acids was observed. Seventy-nine lines were identified which contained significantly high or low levels of one of the five major soybean fatty acids. These lines were advanced to the subsequent generation. Of the 79 lines showing a variant fatty acid profile in the M3, 52 showed clear heritability for the oil composition in the seeds of the subsequent generation. These lines are likely to represent 52 distinct genetic mutations. These mutants may represent new loci involved in the determination of soybean seed oil content or could be new isolates or alleles of previously identified genetic variants for soybean oil composition.

1. Introduction

Soybean oil is composed of five major fatty acids, which are synthesized in the seed during development [1]. In general, wild-type soybeans contain 10–12% palmitic acid (16 : 0), 3-4% stearic acid (18 : 0), 20–25% oleic acid (18 : 1), 50–55% linoleic acid (18 : 2), and 8–10% linolenic acid (18 : 3). Soybean oil has industrial uses and is also a source of vegetable oil for human consumption. For these divergent applications, a modified fatty acid composition profile is sometimes desirable. For example, oil low in the saturated fatty acid palmitate has health benefits in food oils, and reduction of linolenic acid results in increased oxidative stability without the need for hydrogenation [2]. However, for some industrial applications, increased levels of saturated fatty acids (stearate and palmitate) could be desirable [3]. Additional alleles providing new ways to incorporate new and existing oil composition traits could be of use to breeders for many reasons.

Genetic variation is an inexpensive and nontransgenic way to achieve alterations in oil content, which does not require postharvest processing. Soybean mutagenesis has previously been used for the successful improvement of seed oil composition [46]. In addition, novel alleles have been identified in germplasm collections of natural accessions, and the combination of these two sources of genetic variation has resulted in the identification of a number of soybean genotypes with modified oil content [2, 3]. In some cases genetic lesions underlying these traits have been identified at the molecular level. In all published instances for soybean, the mutations are found in genes encoding enzymes that function within the soybean oil biosynthetic pathway. A number of mutations in one of the major enzymes required for fatty acid biosynthesis in the developing soybean seed have been identified in current composition variants. Mutation in the gene encoding the beta-keto-acyl-ACP synthase (KASIIa) results in seeds containing high levels of palmitic acid [7]. The low linolenic acid-containing mutants of soybean carry mutations in the fatty acid desaturase3 (FAD3) genes [810]. Lines carrying mutations in the omega-6 fatty acid desaturase (FAD2-1a) gene results in increased levels of oleic acid [11]. The A6 and FAM94-41 lines carry mutations in the gene encoding the delta-9-stearoyl-ACP-desaturase (SACPD-C), which results in high levels of stearic acid [12].

While many variants with significantly improved oil quality for various applications have been identified and used in breeding programs, previous efforts have focused on screening a population for a specific, favorable altered composition trait. In this project we characterize and describe a large number of mutants with composition distinct from that of the parental genotype. In some cases variants with an unfavorable oil quality profile (one that is unsuitable for any current downstream applications) may add to the present body of knowledge about other genes required for normal oil deposition and soybean seed development and offer new avenues for the genetic improvement of oil content. At the time of this writing, all of the known mutants in oil composition traits have implicated biosynthetic enzymes, and it is likely that other genes encoding proteins that affect the activity or localization of these enzymes could impact oil quality. While such mutations may affect quality to a lesser degree, they may offer scope for more dramatic improvements using biotechnological approaches or in combination with other alleles. Here we describe the screening of a large mutant population for oil composition and the identification of a number of heritable mutations that affect the determination of soybean oil quality.

2. Methods

2.1. Plant Materials and Growth Conditions

The mutant population consists of N-nitroso-N-methylurea (NMU) mutagenized Glycine max cv. Williams-82 (W82) and was described previously [13, 14]. To screen for new, and as-yet unidentified, mutations that contribute to the fatty acid composition of mature soybean seeds, a forward genetic screen was conducted for seed oil composition. A total of 4566 mutagenized lines were subjected to fatty acid profiling by gas chromatography (GC). The fatty acid composition was initially screened by profiling samples of M3 seeds (each sample consisting of three seeds from an individual M2 plant) produced in the field at West Lafayette, IN, USA during the 2005 growing season and comparing these to the W82 wild type grown in parallel. The heritability of the identified variants was subsequently confirmed in the next generation (5 M4 seeds from 3 to 25 individual M3 plants) produced in the field in West Lafayette, IN during the 2008, 2009, or 2010 growing seasons. To estimate the level of environmentally induced variation and to calculate statistical significance for plant-to-plant differences in oil content, ten W82 plants were harvested individually from the field (2009 growing season) and their seed was analyzed for oil content. Level of significance was calculated by Welch’s t-test comparing values from M4 samples to these ten wild-type seed samples.

2.2. Fatty Acid Composition Profiling

Fatty acid profiling was performed at the National Center for Agricultural Utilization Research (Peoria, IL) using the following protocol: fatty acids were extracted into CHCl3 : hexane : MeOH (8 : 5 : 2) for 4 hours at room temperature from three beans cracked with a small hammer. A volume of 0.1 mL 0.35 M methoxide was added to the samples. Samples were analyzed on an Agilent 6890 GC with an Agilent J&W GC column (DB225, 30 m × 0.25 mm id × 0.25 μm film thickness). Data was normalized such that the palmitic, stearic, oleic, linoleic, and linolenic fractions total is 100, and content is expressed as a percentage. Data was analyzed and plots were generated in Microsoft Excel and R (http://cran.r-project.org/). Product names are necessary to report factually on available data. However, the USDA neither guarantees nor warrants the standard of the product and the use of the names implies no approval of the product to the exclusion of others that may also be suitable.

3. Results

Screens for seed fatty acid content were conducted over four years. For each fatty acid, a normal distribution was observed among the 4566 lines profiled. Each distribution showed a number of outliers that potentially represent mutants with extreme levels of one of the major fatty acids (Figure 1). In the population as a whole, a broad range of composition of the five major fatty acids was observed. Specifically, palmitate levels ranged from 5 to 16%, stearate levels ranged from 2% to 15.3%, oleate levels ranged from 14.6 to 46.1%, linoleate levels ranged from 34.6% to 62.6%, and linolenate levels ranged from 4% to 12.9% (Figure 1). For each fatty acid, upper and lower cut-offs were chosen based on the distribution: nine lines with seed palmitate levels above 14% or below 7% were selected for further characterization. Eight lines with stearate content less than 3.5% or greater than 8% were selected, 39 lines with oleate content less than 19% or greater than 30% were selected, 5 lines with less than 40% or more than 60% linoleate were selected, and 18 lines with less than 5% or more than 10% linolenate were selected for further characterization. Mutant M3 seeds of these 79 selected lines were planted in the field during the 2008–2010 growing seasons, and M4 seed samples derived from self-pollinated individual plants from these lines were profiled for oil composition. Data from the M4 seed samples allows determination if the phenotype was transmitted to the following generation, and if the M3 lines were heterozygous or homozygous. Figure 2 shows the reproducibility of fatty acid levels in the M3 and corresponding M4 progeny seed samples for 52 lines that were selected for propagation and crosses on the basis of reproducibility of the phenotype in the M4 samples, taking into consideration that the mutation may not have been homozygous in the M3 generation. In particular, in three lines (lines 8, 34, and 36) evidence of both aberrant and wild-type fatty acid levels in the M4 progeny is apparent, which may indicate that the original M3 isolate was heterozygous. While a range of variation in the levels of fatty acids are observed in the M4 progeny of several of the mutant lines, in contrast in the W82 control plants, variation between individual plants for the level of each fatty acid was less than 0.5%, and minimal year to year variation was observed (Figure 2). For the majority of selected lines, fatty acid values from mutant M4 differed significantly from the wild type ( ).

Complete fatty acid profiles for the M3 lines that were selected for further study are listed in Table 1, along with statistical significance data for oil content differences in the M4 seed samples (see Section 2). A total of three (lines 1–3) low-palmitic acid-containing mutants were identified, two of these contained normal levels of oleic acid. Six (lines 4–9) high-palmitic acid-containing mutants were identified (see Section 4). Four low-stearic acid-containing mutants (lines 10–13) and four mutants with high stearate levels (lines 14–17) were identified. Six low-oleate mutants (lines 18–23) had increased levels of linoleate, and five of these six had elevated levels of linolenate (in all cases these mutants would also have been identified as outliers in the population on the basis of linoleate or linolenate content). The high linoleate mutants (lines 40–42) had marginally reduced levels of oleate. All 16 high-oleate mutants (lines 24–39) had lower than average levels of linoleate. A decrease in linolenate occurred along with an increase in linoleate or oleate in two low-linolenate mutants (lines 43 and 44). All of the mutants with high levels of linolenate (lines 45–52) had reduced levels of oleate.

4. Discussion

Here we have described 52 newly isolated soybean seed composition mutants. These mutants may represent new alleles of previously identified loci that are known to be involved in the determination of seed composition, which will be determined by the sequencing of candidate genes encoding biosynthetic enzymes. Based on the mutagen used to generate these lines, the lesions in these mutants are predicted to be single-base changes. Alternatively, some of these mutants may carry lesions in genes not previously implicated in the determination of soybean seed oil content.

In this study, three mutants (lines 1–3 in Table 1) containing 5-6% palmitate were identified. Other mutants with low levels of palmitic acid (ranging from 5% to 8%) have been identified in previous studies, and at least three nonallelic loci are thought to be involved: fap*, fap1 (C1726), and fap3 (A22) [1519]. The allele fapnc, which is allelic to fap3, has been associated with mutation of the FAT1B gene [19, 20]. Lines 4–9 contain high levels of palmitate, ranging from 12.4% to 16.4%. Mutants containing similarly elevated levels of palmitic acid have been described previously. The allele fap2 has been associated with mutation of the KASIIa gene [7]. Other loci exist which confer a high palmitic acid phenotype which include fap4, fap5, fap6, and fap7. Therefore, there are likely to be other genes present that when mutated influence palmitic acid levels in the seed [16, 2123].

Four mutants containing low levels of stearate (lines 10–13) were isolated in this study, however no mutants with this phenotype are characterized at the molecular level. Four mutants (lines 14–17) with reproducibly high levels of stearate ranging from 7 to 15% were identified in this study. Mutations in the SACPD-C gene are known to result in increased levels of stearic acid ranging from 9% to 30% [5, 12, 24]. There are thought to be at least two genetic loci that contribute to the high stearate phenotype [25].

Sixteen mutants (lines 24–39) were identified in this study that contain high levels of oleic acid. Plants carrying deletions in the FAD2-1a gene contain high levels (up to 50%) of oleic acid in the seeds [11]. A number of other genes affecting oleate levels are known, but have not been characterized at the molecular level at the time of writing. One possible reason for this is that the dependence of oleic acid levels in soybean seeds upon environmental conditions complicates the molecular mapping of these genes [2628]. RNA interference of FAD2-1 (which probably silences multiple FAD transcripts in the seed) results in even higher oleate levels in the range of 70–80%, therefore it is thought that other FAD genes may contribute to the conversion of 18 : 1 to 18 : 2 [29, 30]. Six mutants containing low levels of oleic acid (lines 18–23) were also identified in this study.

Two mutants containing low levels of linolenate (line 43 and 44) were identified in this screen. Mutants with low linolenate (designated fan) have been identified previously, and two independent low linolenic mutants have been shown to contain lesions in the FAD3A gene [5, 10, 31, 32]. When all three FAD3 homologs are downregulated with RNAi in seeds, linolenic acid levels can be further reduced, therefore, it is possible that mutation in FAD3B or FAD3C may also result in low levels of linolenic acid [33]. Eight lines containing 10.5–12.9% linolenic acid were isolated in this study (lines 45–52). One mutant with comparable levels of linolenic acid (12.6%) is known, but the nature of the molecular lesion that causes this phenotype has not been determined [34].

5. Conclusions

Candidate gene sequencing and complementation tests or genetic mapping will be necessary to determine if the lines identified in this study represent new alleles of the known enzymes that control seed oil composition or represent mutations in other genes not yet implicated in the determination of soybean oil quality. These mutations have been isolated from a heavily mutagenized population created for the reverse genetics TILLING approach and carry a number of base changes at additional sites [13]. As with all variants created using a mutation breeding strategy, it will therefore be necessary to introgress these loci into elite cultivars for several generations to eliminate potentially deleterious mutations at other loci. Lines developed in this way will likely be of substantial utility for further development of new soybean varieties with both improved oil content and competitive agronomic characteristics.

Acknowledgments

The authors are grateful to Scott Taylor in the Functional Foods Research Unit at the National Center for Agricultural Utilization Research for sample profiling and to Matthew Hudson for editing the paper. Carrie Anderson, Megan Comerford, Jean Galbraith, Tim Galos, and Peter Kilanowski provided technical assistance. Funding for this work was provided through the USDA-ARS Current Research Information System 3602-21000-004-00D.