|
|
||||||||
a USDA-ARS, Plant Genetics Research Unit, 210 Waters Hall, Univ. of Missouri Columbia, MO 65211
b Dep. of Agronomy, Univ. of Missouri, Columbia, MO 65211
* Corresponding author (BilyeuK{at}missouri.edu)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: PCR, polymerase chain reaction SNP, single nucleotide polymorphism
| INTRODUCTION |
|---|
|
|
|---|
Breeding based on chemical phenotype for single seed composition traits has been successful in lowering linolenic acid from wild-type levels of 60 to 100 g kg1 oil to as low as 10 g kg1 oil for lines containing three independent mutations (Fehr et al., 1992; Rahman et al., 1998; Ross et al., 2000). Because linolenic acid content in soybean is controlled by multiple genes, a more efficient breeding strategy would be to use molecular markers for rapid introgression of the trait in a backcross breeding strategy to develop elite lines (Bilyeu et al., 2005). Molecular markers for the lowered linolenic acid trait would also allow the combination of multiple traits in one line where a backcrossing strategy could be utilized to capture mutant alleles of a low linolenic acid donor parent with a recurrent parent that has other desirable seed composition traits with or without associated molecular markers.
Genetics of the lowered linolenic acid trait have received considerable attention. After the original fan mutation corresponding to lower linolenic acid levels was detected in line C1640 (Wilcox and Cavins, 1985; Wilcox and Cavins, 1987), other independent loci (fan2, fan3, and fanx) were also identified in lines with further reductions in linolenic acid content (Fehr et al., 1992; Rahman et al., 1998; Ross et al., 2000). In an effort to provide a molecular understanding of the low linolenic acid trait, we identified and characterized the soybean homologs of the Arabidopsis omega-three fatty acid desaturase gene, FAD3 (Bilyeu et al., 2003; Yadav et al., 1993). At least three versions of this gene exist in the soybean genome, GmFAD3A, GmFAD3B, and GmFAD3C, and GmFAD3A was shown to be Fan (Bilyeu et al., 2003; Anai et al., 2005). Recently, we determined the molecular genetic basis for the low linolenic acid soybean line CX151244, which had mutant alleles of both GmFAD3A and GmFAD3C (Bilyeu et al., 2005). The GmFAD3A gene was shown to have a greater impact on seed linolenic acid levels than GmFAD3C, consistent with higher expression of GmFAD3A in developing seeds (Bilyeu et al., 2005; Bilyeu et al., 2003).
The lowest level of linolenic acid reported for a soybean line is 1% of the oil fraction (10 g kg1 oil) described for line A29 (Ross et al., 2000). A29, which is not related to other low linolenic acid lines with reported molecular information (CX151244, J18, M5, and M24 [Anai et al., 2005; Bilyeu et al., 2005]), was developed by combining three independent mutations: fan from line A5 (Hammond and Fehr, 1983), fan2 from A23, and fan3 from a mutagenized derivative of line A89144003 (Ross et al., 2000). The fan mutation in line A5 is caused by a GmFAD3A deletion (Bilyeu et al., 2003; Byrum et al., 1997). The very low linolenic acid lines A29, IA3017, and IA3018 are related by pedigree (W. Fehr, personal communication, 2005). The objective of this research was to investigate the FAD3 genes in soybean line A29 to identify possible mutations and determine if the mutant alleles were associated with the 1% linolenic acid trait.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Plant Growth
The low linolenic acid soybean lines A29 (10 g kg1 oil), IA3017 (10 g kg1 oil), and IA3018 (26 g kg1 oil) were chosen as candidates for potential FAD3 gene mutations. Low linolenic soybean lines A29, IA3017, and IA3018 were provided by Dr. Walter Fehr, Iowa State University (Ross et al., 2000). The low linolenic acid line 2721 (30 g kg1 oil) was an F3:F6 line produced from a cross between Pana and CX151244. Line 2721 was homozygous for the CX151244 mutant alleles of GmFAD3A and GmFAD3C (Bilyeu et al., 2005).
SS976946 is an experimental breeding line with normal fatty acid profiles (77 g kg1 oil), resistance to multiple soybean cyst nematode HG types, and appropriate maturity (D. Sleper, unpublished data, 2003). A cross of SS976946 x IA3017 produced F1 seeds in the summer of 2003 at the Bradford Research and Extension Center (BREC), located near Columbia, MO. F1 plants were grown in Guacima, Costa Rica, and F2 seed was chipped for fatty acid analysis with the remainder of the seed planted and grown in Nazareth, Costa Rica, during the fall of 2003 and spring of 2004. Chips were shipped to BREC and analyzed for linolenic acid concentration. Single-plant threshed F2:F3 seed of lines containing the lowest linolenic acid content (18 lines containing up to linolenic acid concentrations of 15 g kg1 oil) was returned to BREC for planting progeny rows in summer 2004. Because no lines were recovered which were homozygous for all three mutations, five individual aabbcc F2:F3 plants derived from an aabbCc F2 parent were used as pollen donors for a SS976946 x F2:F3 backcross (SS976946 x [SS976946 x IA3017]). The resulting three confirmed BC1F1 seeds were germinated in moist germination packets (Mega International, St. Paul, MN) and transferred to soil for growth in controlled chambers set at 27.5/23°C day/night with 14.5 h daylength at a light intensity of 750 µmol m2 s1. Eighty seeds were produced, and fatty acid and genotype data were collected for 61 individuals.
Seeds of a A29 x 2721 cross were produced in Costa Rica in 2003. F1 (1 seed), F2 (62 seeds), and F2:F3 seeds were germinated in moist germination packets and transferred to soil for growth in controlled chambers set at 27.5/23°C day/night with 14.5 h daylength at a light intensity of 750 µmol m2 s1. Twelve lines were identified in the F2 generation to have homozygous mutant alleles for GmFAD3B and GmFAD3C, and a subset (8) of those lines were followed to the F3 generation to test for segregation of GmFAD3A alleles.
Genotype Analysis
Detection of mutant alleles followed the procedure as described by Bilyeu et al. (2005) with PCR amplification of the genomic region encompassing the single nucleotide polymorphism (SNP) followed by restriction enzyme digestion of products. When primer sequences are indicated, the forward primer is listed first, and all primers are listed in the 5' to 3' orientation.
For the A29 GmFAD3B allele assay, amplification primers were B932: AGCCACAGAACTCACCATCAA and IABrev: TGGCACAGTGAATCTAATG. Following amplification conditions as described for the CX151244 alleles (Bilyeu et al., 2005), products were digested with 2.5 U HpaI (New England Biolabs, Beverly, MA) per reaction and resolved on agarose gels. Wild-type alleles produced fragments of 196 and 162 bp while A29-derived mutant alleles produced a single 358-bp fragment.
For the A29 GmFAD3C allele assay, amplification primers were IACleft: TTGGATCAACAACATTCACCA and IACright: CATCACATGTTTGTGGTCTTGA. Following amplification, products were digested with 5 U BccI (New England Biolabs) per reaction and resolved by McSNP analysis (Ye et al., 2002) or on agarose gels. Wild-type alleles produced fragments of 199 and 56 bp while A29-derived mutant alleles produced a single 250-bp fragment.
For the GmFAD3A deletion assay, a Taqman assay was used in a reaction with primers that detected both the GmFAD3A gene (3AD1/3AIX amplification primers [Bilyeu et al., 2005]; dual-labeled probe [FAM/BHQ1] WTprobe: TCGCGGCAAGGTAACAAAAA) and the PEPC16 gene (PEPC16f: TTCCTTTATCAGAAATAACGAGTTTAGCT, PEPC16r: TGTCTCATTTTGCGGCAGC; dual-labeled probe [VIC/BHQ1] PEC16probe: CCCTCCCCTGTACCCATGTTTCCATTATAA [Tuteja et al., 2004]). Reactions for PCR included 0.33 µM of each amplification primer, 0.21 µM of each dual labeled probe, and 1x QuantiTect Probe PCR mix (Qiagen Inc., Valencia, CA). PCR conditions were 95°C for 12 min, 44 cycles of 95°C for 20 s, 59°C for 20 s, and fluorescence read. In samples for which there was amplification of the PEPC16 target, fluorescence was scored for the amplification of the GmFAD3A target.
Templates for all genotype PCRs consisted of 2-mm washed FTA (Whatman, Clifton, NJ) card punches prepared from leaves according to the manufacturer's instructions. Genomic DNA was isolated for control and parental lines either using the DNeasy Plant Mini Kit (Qiagen, Inc.) or as described previously (Bilyeu et al., 2003) and used at 5 to 50 ng per PCR reaction.
Phenotype Analysis
Fatty acid analysis was conducted as described for seed chips or whole crushed seed (Beuselinck et al., 2006; Bilyeu et al., 2005). For the F2:F3 seed fatty acid determination, three or five F2:F3 seeds were crushed and analyzed individually. The concentration of linolenic acid in the seed sample was determined as a percentage of the total fatty acids of extracted oil by lipid gas chromatography of fatty acid methyl esters.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
When cDNA from line A29 was initially used as template for RT-PCR, the GmFAD3B gene reactions failed to produce a product. Therefore, genomic DNA from line A29 was used to amplify the GmFAD3B gene region, and the product was cloned, sequenced, and compared to the partial GmFAD3B genomic sequence present in GenBank (accession AX088031). Among several minor differences in intron sequences, a SNP was identified in the 3' splice site consensus (Brown et al., 1996) sequence preceding exon 4 (Fig. 1 ). Improper splicing was subsequently detected in mRNA from this region that would lead to a frameshift and premature termination of the protein product (data not shown). An assay was designed to distinguish wild-type and mutant alleles of GmFAD3B based on the presence or absence of a HpaI restriction enzyme site for wild-type and A29 derived mutant alleles, respectively (Fig. 1C).
|
|
Using Molecular Markers in a Backcross Breeding Strategy
Part of our interest in developing molecular markers for the low linolenic acid trait is to provide the tools to accelerate introgression of traits into elite soybean cultivars. The use of molecular markers supercedes the need for extensive chemical analysis and facilitates the capture of the trait in a backcrossing program that can rely on selection of unrelated traits. Toward this end, we initiated a cross between SS976946 (D. Sleper, unpublished data, 2003) and IA3017. Following the production of F1 seed, a typical breeding strategy was used in which F1 plants were grown, allowed to self-fertilize, and the F2 seed was chipped for fatty acid analysis. Agronomic selections and selections based on fatty acid analysis were made on F2 plants. As an initial genotyping screen, two individuals from each F2:F3 family were genotyped. No F2:F3 families were identified that contained all three homozygous mutant alleles of the FAD3 genes. However, individual plants were identified by genotype to be homozygous mutants, and these individual homozygous mutant F2:F3 plants were used as pollen donors in backcrosses to SS976946.
The association of genotype and phenotype was assessed in field grown plants derived from a single cross (SS976946 x IA3017). F2:F4 seeds produced on field grown selfed F2:F3 plants of selected plants with known genotypes were tested for fatty acid composition. While the recurrent parent germplasm contained an average of 7.7% linolenic acid (77.0 g kg1 oil), the five F2:F3 lines identified with the homozygous mutant genotype for all three genes produced F2:F4 seed with an average of 1.3% linolenic acid (13.0 g kg1 oil). Grown in the same field environment, IA3017 produced seeds with an average of 1.2% linolenic acid (12.0 g kg1 oil).
Over the subsequent cycle of backcross breeding (SS976946 x [SS976946 x IA3017]), plants were grown in controlled environment chambers. BC1F1 plants (derived from SS976946 x F2:F3) were confirmed using molecular markers. The BC1F2 seeds were chipped for fatty acid analysis and germinated for further growth and genotyping (Fig. 3 ). One BC1F2 line was recovered with the homozygous mutant genotype for all three GmFAD3 genes, and it contained the lowest concentration of linolenic acid (12 g kg1 oil). Because 27 genotypes are possible from the segregation of the three GmFAD3 loci, and the assay for the GmFAD3A deletion does not distinguish wild-type from heterozygous individuals, the association of genotype and phenotype is not entirely clear. In general, samples homozygous for the GmFAD3A deletion produce the lowest linolenic acid phenotype, regardless of the genotype at GmFAD3B and GmFAD3C. In addition, mutations in GmFAD3B and GmFAD3C also appear to additively contribute to lowered linolenic acid levels as demonstrated by the relatively low linolenic acid phenotype for samples with at least one copy of a wild-type GmFAD3A allele in combination with homozygous mutations in GmFAD3B and GmFAD3C (Fig. 3, samples 8, 12, 15, 17, 19, and 36, gray circles). Conversely, samples with at least one copy of a wild-type GmFAD3A allele in combination with homozygous wild-type GmFAD3B and GmFAD3C contained some of the highest linolenic acid levels.
|
|
The association data presented here supports a model in which all three GmFAD3 genes contribute to the total enzymatic capacity to produce linolenic acid in the seed oil. Consistent with our earlier work, the mutation in the GmFAD3A gene described here caused a greater reduction in linolenic acid concentration than mutations in GmFAD3B or GmFAD3C (Bilyeu et al., 2005; Bilyeu et al., 2003).
Recombination of Mutant Alleles
Together with the previously identified mutations in GmFAD3A and GmFAD3C present in the low linolenic acid lines C1640 and CX151244 (Bilyeu et al., 2005; Chappell and Bilyeu, 2006), there are now a total of six available mutant alleles with perfect (i.e., specific to the causative mutation) molecular markers at three soybean GmFAD3 loci (Table 1). The mutant GmFAD3 genotype in A29 can be designated aaA29bbA29ccA29 and distinguished from the CX151244 derived mutant genotype, designated aacxBBcccx. No mutations were identified in the CX151244 GmFAD3B gene (Bilyeu et al., 2005).
|
We initiated a cross between A29 (aaA29bbA29ccA29) and 2721 (aacxBBcccx, Bilyeu et al., 2005), and analyzed progeny for genotype and phenotype in the F2 and F2:F3 generations to select lines homozygous for different combinations of the three SNP-containing mutant alleles (aacxbbA29cccx or aacx bbA29ccA29). As expected for segregation of only mutant alleles at GmFAD3A and GmFAD3C and mutant and wild-type alleles of GmFAD3B, the range of linolenic acid concentrations in F2 seeds was narrow (12 to 51 g kg1 oil). Two F2 individuals (lines designated B152 and B162) were confirmed in the F3 generation to be homozygous for both CX151244-derived mutant alleles and homozygous for the A29-derived GmFAD3B allele, aacxbbA29cccx. One F2 individual (designated line B115) had the alternate genotype, aacx bbA29ccA29. The linolenic acid contents for these lines are listed in Table 2, along with a line with the aaA29bbA29cccx genotype. These results indicate that the available mutant alleles can be substituted for each other and in combinations that enable SNP detection for all three FAD3 genes as well as stable inheritance of the 1% linolenic acid phenotype.
|
| CONCLUSIONS |
|---|
|
|
|---|
The mutant alleles identified as part of this project were readily substituted with other mutant alleles of GmFAD3A and GmFAD3C that had been previously identified (Bilyeu et al., 2005). Breeding programs that have already incorporated a fan mutation (an allele of GmFAD3A) to lower linolenic acid concentration now have the option of accomplishing further reductions in linolenic acid by incorporating mutant alleles of GmFAD3B and GmFAD3C. In addition, the novel combinations of mutant GmFAD3 genes that impart a 1% (10.0 g kg1 oil) linolenic acid phenotype developed as part of this project provide new resources to the soybean breeding community.
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
Received for publication March 8, 2006.
| REFERENCES |
|---|
|
|
|---|
linolenic acid mutants and elucidation of their contribution to fatty acid biosynthesis in soybean seed. Plant Sci. 168:16151623.[CrossRef]This article has been cited by other articles:
![]() |
W. R. Fehr Breeding for Modified Fatty Acid Composition in Soybean Crop Sci., December 18, 2007; 47(Supplement_3): S-72 - S-87. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||