|
|
||||||||
9 Stearoyl-ACP Desaturases
a Microbiology Dep., North Carolina State Univ., Raleigh, NC 27695
b Crop Science Dep., North Carolina State Univ., Raleigh, NC 27695
c USDA-ARS Soybean and Nitrogen Fixation Unit and Plant Pathology Dep., North Carolina State Univ., Raleigh, NC 27695
* Corresponding author (greg_upchurch{at}ncsu.edu)
| ABSTRACT |
|---|
|
|
|---|
9 stearoyl acyl-carrier protein desaturase (SACPD) gene of soybean [Glycine max (L.) Merrill] encodes a soluble enzyme that converts stearic to oleic acid. Understanding the regulation of SACPD expression and enzyme activity are thus important steps toward developing soybean lines with altered stearic or oleic acid content. Using primers designed to a G. max SACPD cDNA sequence, a 3648-bp product was cloned and sequenced from the genome of cultivar Dare. Comparison of the third SACPD exon protein sequence with other available Glycine SACPD sequences revealed unique amino acid variability at positions 310 and 313. Sequence-specific primers were designed for Real-time RT-PCR (reverse transcriptase-polymerase chain reaction) for this region of exon 3. Diagnostic and specific products were recovered with these primers using Dare cDNA template and Dare genomic DNA. Sequencing of a second genomic clone from Dare confirmed that there were two SACPD genes, designated A and B, in this cultivar. Survey of the genomes of 51 soybean lines and cultivars with PCR and the gene-specific primers indicated that all 51 had both A and B. Differences between SACPD-A and -B transcript abundance in soybean tissues, while quantifiable, were not dramatic. SACPD-A and -B transcript accumulation for three seed developmental stages between R5 and R6 was essentially equal. Biochemical analysis of the proteins encoded by these two SACPD genes may reveal whether the amino acid variability uncovered in this study has any relation to enzyme activity.
Abbreviations: EST, expressed sequence tag FAD, fatty acid desaturase RT-PCR, reverse transcription-polymerase chain reaction SACPD, stearoyl acyl carrier protein desaturase
| INTRODUCTION |
|---|
|
|
|---|
Increasing the stearic acid content of soybean oil is desirable for certain food-processing applications since increased stearic acid content offers the potential for the production of solid fat products without hydrogenation (Spencer et al., 2003). This would reduce the content of transfat in food products and reduce a current health concern. High stearic acid oil could also be used as a replacement for tropical oils that are also high in palmitic acid, which poses a heart disease risk.
Stearic acid concentration in soybean is genetically determined by alterations at the Fas locus. In addition to FAM9441 (Pantalone et al., 2002), a high stearic acid line carrying a natural mutation (fasnc), five other soybean germplasm lines have been reported to carry mutated stearic acid alleles that increase stearic acid: A6, fasa (Hammond and Fehr, 1983); FA41545, fasb (Graef et al., 1985a); A81606085, fas (Graef et al., 1985b); KK-2, st1 (Rahman et al., 1997); and M25, st2 (Rahman et al., 1997). Fasa (30% stearic acid), fasb (15% stearic acid), and fas (19% stearic acid) are allelic and represent different mutations in the same gene (Burton et al., 2004; Spencer et al., 2003). One candidate gene, although unproven, for the Fas locus in soybean, is the
9 stearoyl-acyl carrier protein (ACP) desaturase (Wilson, 2004). The
9 stearoyl-ACP desaturase enzyme (SACPD), through the insertion of a double bond at C9, converts stearic to oleic acid. Thus, SACPD occupies a key position in C18 fatty acid biosynthesis since perturbation of SACPD gene expression and/or enzyme activity may modulate the relative level of both stearic and oleic acid in soybean oil. Down-regulated SACPD expression or enzyme activity could produce oil with greater stearic acid, while increased expression of SACPD or enzyme activity could produce oil with greater oleic acid content.
In this study, we determined that there are two SACPD genes in soybean, a situation not previously recognized. The objective of this research was to characterize SACPD gene structure, determine the distribution of the SACPD genes in Glycine, and analyze SACPD expression levels in stages of seed development.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Primers were also designed to a region of variability in exon 3 of SACPD. PCR with these primers produced products that discriminated between SACPD gene A (133 bp) and gene B (111 bp). These gene-specific primers were used to type soybean cultivars as well as to quantify steady state mRNA accumulation by real-time RT-PCR. Gene-specific primers for SACPD-A (primer set A) were 5'CCTGTTTGATAACTACTCTGCC3'-forward and 5'TCTTCCCTCACCTGAAAGTCCG3'-reverse. Gene-specific primers for SACPD-B (primer set B) were 5'CCTGTTTGATAGCTACTCTTCG3'-forward and 5'GTTAGCTGCTCCACCTCC3'-reverse. Primers for the soybean housekeeping gene actin were 5'GAGCTATGAATTGCCTGATGG3' forward and 5'CGTTTCATGAATTCCAGTAGC3' reverse derived from GenBank accession number U60500 (Moniz de Sa and Drouin 1996). Amplification of actin was done according to the protocol outlined above for SACPD. Amplification reactions (25 µL) for PCR with SACPD gene-specific primers contained 50 ng of EcoRI restricted genomic DNA, 2 µM each primer, 250 µM dNTPs, 1 unit Taq DNA polymerase (Fisher, Atlanta, GA), and 1x PCR buffer [300 mM Tris-HCl, 75 mM (NH4)2SO4, 7.5 mM MgCl2, pH 8.5]. Reaction conditions were 94°C for 2 min (1 cycle), then 94°C for 1 min, 66.8°C for 2 min, 72°C for 3 min (40 cycles), then 72°C for 7 min.
PCR amplification reactions were performed in a MJ Research PTC-100 thermocycler (Watertown, MA). PCR products were analyzed on ethidium bromide stained, 1% (w/v) agarose gels after electrophoresis with a 100-bp DNA ladder molecular weight markers (Invitrogen, Carlsbad, CA). Amplicons were cloned into the sequencing vector pCR 2.1 with the TOPO TA cloning kit supplied by Invitrogen (Carlsbad, CA). Both strands of each PCR clone were sequenced at the Iowa State University Biotechnology Center, Ames, IA. Contigs were assembled by Vector NTI Advance 9.1 (Invitrogen). Assignment of open reading frames, translation to amino acid sequences, and sequence comparisons were done with the NCBI tools for data mining (Gish and States, 1993; Altschul et al., 1990) and with Chromas (http://www.technelysium.com.au/; verified 18 November 2005). Sequences were aligned by Blossum62 Divide-and-Conquer Multiple Sequence Alignment v. 1.0 (BiBiServ, Bielefeld University Bioinformatics Server). Assignment of SACPD exons and introns was done by MacVector v. 7(IBI, New Haven, CT) and NetGene2 WWW Server (Center for Biological Sequence Analysis, The Technical University of Denmark).
Genomic DNA Gel Blot Analysis
Genomic DNA (6 µg) was digested with restriction enzymes, electrophoresed through a 0.8% (w/v) agarose gel, and transferred to Nytran Plus (Schleicher & Schuell, Keene, NH) membranes. The membranes were probed with [
-32P] dCTP-labeled DNA fragments that correspond to an exon 3 region for either SACPD-A or SACPD-B and exposed to X-ray film. Membrane hybridization and washing procedures were those described for Ultrahyb solution (Ambion, Austin, TX).
RNA Isolation and Processing
Glycine max Dare plants were grown under greenhouse conditions and plant tissues (pod, leaf, lateral root) harvested at 18 and 35 d after flowering (DAF). Dare plants were also grown under controlled day/night temperatures of 26/22°C and seeds were harvested at four stages (18, 23, 28, and 35 DAF) of development between R5 and R6. Tissues and seeds were quickly frozen in liquid nitrogen and stored at 80°C until RNA was extracted. Samples were pooled from three plants, and RNA was isolated from 100 mg of frozen powdered tissue with the Qiagen RNeasy Plant Mini kit (Valencia, CA) following the manufacturer's protocol. The RNA extraction was repeated on two other pooled samples. RNA samples were DNase treated with Ambion DNA-free (Austin, TX) according to the manufacturer's protocol. RNA concentrations were determined spectrophotometrically with absorbance at 260 nm. Samples were diluted to 50 ng/µL in sterile water and aliquots stored at 80°C until use. To verify RNA integrity, 500 ng of total RNA of each sample was examined on a 1% (w/v) agarose gel after electrophoresis and staining with ethidium bromide.
Real-Time Reverse Transcriptase PCR
Real-time reverse transcriptase PCR (Winer et al., 1998; Bustin, 2002) was performed with the iCycler iQ (Bio-Rad, Hercules, CA) using the QuantiTech SYBR Green RT-PCR kit (Qiagen, Valencia, CA). Each reaction contained 12.5 µL of 2x SYBR Green PCR master mix, 250 nM forward and reverse gene-specific primers for SACPD gene A or gene B, 0.25 µL MgCl2 (25 mM), 0.25 µL RT mix, 250 ng RNA, and nuclease-free water to a final volume of 25 µL. The reactions were performed in a 96-well plates (0.2-mL tube volume) sealed with optical tape. Reverse transcription was performed at 50°C for 30 min followed by 95°C for 15 min to inactivate the reverse transcriptases and to activate the HotStar Taq DNA polymerase. PCR amplification involved 45 cycles of 15 s at 94°C, 30 s at 60°C, and 30 s at 72°C followed by melt curve analysis over a 10°C temperature gradient at 0.05°C s1 from 78 to 88°C. Duplicate reactions were done for each sample. Steady state transcript levels for SACPD genes were mathematically determined by comparison of individual cycle threshold (Ct) values with a standard curve generated from serial dilutions of a PCR standard (Peirson et al., 2003) from soybean genomic DNA. PCR efficiency ranged from 93 to 98%, and negative control reactions did not produce any products. Targets (copy number), initially determined per microgram of total RNA isolated, were normalized as a percentage of soybean actin gene expression level.
Determination of Seed Fatty Acid Composition
Fatty acid methyl esters (FAME) of mature soybean seed samples were prepared by acid methanolysis. Seed tissue, ground to a powder, was heated at 85°C for 90 min in a (v/v) 5% HCl95% methanol solution. FAME was partitioned 2x into hexane and transferred to 2-mL vials for analysis. The FAMEs were separated by gas chromatography with an HP 6890 GC (Agilent Technologies, Inc., Wilmington, DE) equipped with a DB-23 30-m x 0.53-mm column (same source). Operating conditions were 1-µL injection volume, a 20-to-1 split ratio, and He carrier gas flow of 6 mL min1. Temperatures were 250, 200, and 275°C for the injector, oven and FID, respectively. Chromatograms were analyzed by HP ChemStation software.
| RESULTS |
|---|
|
|
|---|
9 Stearoyl-ACP Desaturase
9 stearoyl-ACP desaturase cDNA (accession no. L34346) from soybean. Sequencing and analysis of these fragments revealed the Glycine SACPD gene structure (Fig. 1A
) to be the following: a 111-bp exon 1 sequence encoding a putative 37 amino acid transit sequence; a 1763-bp intron 1; a 504-bp exon 2 sequence encoding 168 amino acids; a 423-bp intron 2; and a 618-bp exon 3 sequence encoding 206 amino acids. The nucleotide and amino acid sequence data of the Dare SACPDs can be found at the GenBank database as accession numbers AY885234 (SACPD-A) and AY885233 (SACPD-B).
|
Alignment of Other Soybean SACPD Exon 3 Sequences
To assess the significance of the amino acid variability found in the SACPD exon 3s of Dare, other Gm SACPD exon 3s were PCR cloned and sequenced. We sequenced an exon 3 from Gm cultivar Bragg and an exon 3 from N013544, a mid-oleic acid Gm line. To these sequences were added 5 soybean cDNA exon 3-derived sequences from the GenBank database. Nucleic acid alignments for bp positions 3183 to 3211 in exon 3, depicted in Fig. 2A
, show that the nine SACPD exon 3 sequences can be differentiated into two groups: Dare-A, Bragg, N013544, AW 755581, and BG 882246 (Williams) designated group A, and Dare-B, AI 941223 (Williams), BI 471396 (also Bragg), and BG 363272 designated group B. Amino acid alignments of the corresponding region of the exon 3s (position 307321) revealed that variation in amino acids at position 310 and 313 permit exactly the same group A and B differentiation. Moreover, both SACPD-A and -B sequences were represented for cultivars Bragg (our sequence and BI 471396) and Williams (BG 882246 and AI 941223). On the basis of these results, we hypothesized that soybean possesses at least two SACPD genes.
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
9 stearoyl-ACP desaturase of soybean, like all the soluble desaturases using acyl-ACP substrates, is localized to the stroma fraction of plastids in developing seeds (Murphy and Piffanelli, 1998). A short N-terminal transit peptide of 37 amino acids was identified for the soybean SACPDs that is presumably responsible for stroma targeting (Fig. 1A). Most of our knowledge about plant
9 stearoyl-ACP desaturases (SACPDs) comes from the study of the soluble enzyme from castor seed (Lindqvist et al., 1996). Structure analysis of the crystallized protein shows it to be a µ-oxo-bridged di-iron enzyme that belongs to the structural class I of large helix bundle proteins that catalyzes the NADPH and O2dependent insertion of a cis-double bond between C-9 and C-10 positions in stearoyl-ACP (Moche et al., 2003). The enzyme is a homodimer with each mature subunit of 41.6 kDa containing an independent binuclear iron cluster. At the core of the desaturase structure, two iron atoms are coordinated within a central four helix bundle in which the motif (D/E)-E-X-R-H is present in two of the four helices (Ohlrogge and Browse, 1995). During the desaturation reaction, the two-electron reduced, di-iron center binds oxygen and the high valent iron-oxygen complex formed abstracts hydrogen from the substrate CH bond (Whittle and Shanklin, 2001).
Inspection of G. max SACPD exon 3 protein sequences derived from cDNA data suggested that soybean contains two SACPD genes (Fig. 2A and B). On the basis of the amino acid variable region of exon 3, what we now call SACPD-A and -B were found in cultivar Williams. Further sequence analysis of exon 3 clones from cultivars Dare and Bragg indicated that they also contained SACPD-A and -B. Sequence analysis of Dare genomic clones confirmed that there are two different, soluble SACPD genes in this cultivar. Amino acid variability was found only in the Dare protein sequence of exon 3 and not in the transit peptide region (exon 1) or in exon 2. Previously, multiple plant genes were identified only for microsomal desaturases. For example, three different microsomal
12 oleate desaturase genes (FAD2s) were reported for sunflower (Martinez-Rivas et al., 2001), and two microsomal FAD2s were reported each for soybean (Heppard et al., 1996) and for cotton (Liu et al., 1997), while three different microsomal
-3 desaturase genes (FAD3s) were reported for soybean (Bilyeu et al., 2003).
Rather dramatic differences in microsomal desaturase gene expression levels and gene distribution have been found in soybean. Of the two soybean FAD2 genes, FAD21 was expressed specifically in seeds, and FAD22 was expressed in all tissues (Heppard et al., 1996). Only one of the three soybean FAD3 genes, FAD3A, was predominatly expressed in developing seeds (Bilyeu et al., 2003). In the same study, Bilyeu et al. (2003) found that the low linolenic acid breeding line A5 contained two of the FAD3 genes but lacked a third gene, FAD3A. In our study, we found that the differences between the transcript abundance of the soluble SACPD-A and -B in soybean tissues, while quantifiable, were not dramatic (Fig. 4 and 5). However, both soybean SACPDs were much more highly expressed in seeds than in roots.
We thought that Glycine linescultivars might be identified that lacked one of the SACPD genes; however, both A and B genes were found in the genomes of all 51 Glycine linescultivars examined (Table 1). Group I consisted of eight G. max cultivars and lines of varying maturity, including the high stearic acid mutant A6 and the midoleic acid line NO13544, group II consisted of eight maturity group V G. soja lines varying in oleic acid content, and group III consisted of the 35 soybean cultivars of Gizlice et al. (1994). Group III genotypes were chosen because they define the genetic base of North American soybean cultivars and represent 95% of the genes found in modern cultivars. Although GenBank accession L34346 may be indicative of a third soybean SACPD gene, efforts to find this sequence in cultivar Dare were unsuccessful. Since the primary sequence of L34346 fits the amino acid substitution pattern of neither SACPD-A nor -B exactly, L34346 may represent an allelic form of one of the two SACPD genes.
The manipulation of fatty acid desaturases to achieve a desired fatty acid composition in soybean oil has a strong rationale. Support for this rationale comes from research that has shown that downregulation of the
-12 fatty acid desaturase gene FAD21 elevates oleic acid content in the oil. Transgenic seeds with oleic acid content of approximately 75 to 80% of the total oil have been recovered after this gene was silenced in somatic embryos (Kinney, 1997) or the gene transcript was ribozyme-terminated (Buhr et al., 2002). As mentioned previously, the low linolenic acid breeding line A5 was found to lack the FAD3A gene at the Fan locus in soybean (Bilyeu et al., 2003). In addition, FAM9441 (Spencer et al., 2003), a high stearic acid line, was found to carry a natural mutation, fasnc, at the Fas locus. On the basis of our findings, we will continue to characterize the enzymatic activity of the two soybean
9 soluble desaturases from both the wild-type Dare cultivar and the high stearic acid mutant line A6. Results from these experiments may provide a means to achieve the stable production of high stearic acid soybean oil.
| ACKNOWLEDGMENTS |
|---|
Received for publication June 30, 2005.
| REFERENCES |
|---|
|
|
|---|
-6 desaturase genes in soybeans. Plant Physiol. 110:311319.[Abstract]
9stearoyl0acyl carrier protein desaturase from castor seed and its relationship to other di-iron proteins. EMBO J. 15:40814092.[Web of Science][Medline]
-6 desaturase genes in cotton (Gossypium hirsutum L.). p. 383385. In J.P. Williams et al. (ed.) Physiology, biochemistry and molecular biology of plant lipids. Kluwer Academic Publishers, Dordrecht, Netherlands.
9 stearoyl-acyl carrier protein desaturase. J. Biol. Chem. 278:2507225080.
916:0-Acyl Carrier Protein (ACP) desaturase specificity based on combinatorial saturation mutagenesis and logical redesign of the castor
918:0-ACP desaturase. J. Biol. Chem. 276:2150021505.This article has been cited by other articles:
![]() |
P. Zhang, J. W. Burton, R. G. Upchurch, E. Whittle, J. Shanklin, and R. E. Dewey Mutations in a {Delta}9-Stearoyl-ACP-Desaturase Gene Are Associated with Enhanced Stearic Acid Levels in Soybean Seeds Crop Sci., November 24, 2008; 48(6): 2305 - 2313. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. E. Byfield and R. G. Upchurch Effect of Temperature on Microsomal Omega-3 Linoleate Desaturase Gene Expression and Linolenic Acid Content in Developing Soybean Seeds Crop Sci., November 7, 2007; 47(6): 2445 - 2452. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. E. Byfield and R. G. Upchurch Effect of Temperature on Delta-9 Stearoyl-ACP and Microsomal Omega-6 Desaturase Gene Expression and Fatty Acid Content in Developing Soybean Seeds Crop Sci., July 30, 2007; 47(4): 1698 - 1704. [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 | |||