Original Research ARTICLE
Characterization of a stearoyl-acyl carrier protein desaturase gene family from chocolate tree, Theobroma cacao L
- 1Huck Institutes of the Life Sciences, The Pennsylvania State University, University Park, PA, USA
- 2Department of Plant Science, The Pennsylvania State University, University Park, PA, USA
In plants, the conversion of stearoyl-ACP to oleoyol-ACP is catalyzed by a plastid-localized soluble stearoyl-acyl carrier protein (ACP) desaturase (SAD). The activity of SAD significantly impacts the ratio of saturated and unsaturated fatty acids, and is thus a major determinant of fatty acid composition. The cacao genome contains eight putative SAD isoforms with high amino acid sequence similarities and functional domain conservation with SAD genes from other species. Sequence variation in known functional domains between different SAD family members suggested that these eight SAD isoforms might have distinct functions in plant development, a hypothesis supported by their diverse expression patterns in various cacao tissues. Notably, TcSAD1 is universally expressed across all the tissues, and its expression pattern in seeds is highly correlated with the dramatic change in fatty acid composition during seed maturation. Interestingly, TcSAD3 and TcSAD4 appear to be exclusively and highly expressed in flowers, functions of which remain unknown. To test the function of TcSAD1 in vivo, transgenic complementation of the Arabidopsis ssi2 mutant was performed, demonstrating that TcSAD1 successfully rescued all AtSSI2 related phenotypes further supporting the functional orthology between these two genes. The identification of the major SAD gene responsible for cocoa butter biosynthesis provides new strategies for screening for novel genotypes with desirable fatty acid compositions, and for use in breeding programs to help pyramid genes for quality and other traits such as disease resistance.
Theobroma cacao L. (cacao) is an understory tropical tree domesticated in the Amazon basin and today widely cultivated in West Africa, Central and South America, and Southeast Asia (Argout et al., 2011). Cacao pods, containing around 40 seeds, are harvested approximate 20 weeks after pollination, and the seeds contain about 50% total lipids (cocoa butter), which provides a main raw ingredient for chocolate manufacturing as well as ingredients for pharmaceutical and cosmetic products (Liendo et al., 1997). Notably, cocoa butter composition consists of almost equal amount of palmitic acid (16:0), stearic acid (18:0), and oleic acid (18:1n−9), the composition of which determines its unique melting temperature very close to human body temperature, thus providing the smoothness and mouth feel of chocolate, as well as the creamy texture of cosmetics on skin (Liendo et al., 1997). Remarkably, during the 20 weeks of cacao seed development and maturation, dramatic transitions of fatty acid profiles occur within the only 20 days when a polyunsaturated fatty acid (linoleic acid (18:2n−6) and α-linolenic acid (18:3n−3))-rich membrane-like profile transitions to a final fatty acid content rich in saturated and monounsaturated-storage lipids (Patel et al., 1994). Moreover, as the dominant form of unsaturated fatty acid, 18:1n−9 gradually accumulates in developing cacao seeds, resulting in a significant impact on the physical properties of cocoa butter owing to its much lower melting temperature (16°C) compared to the other saturated fatty acids (16:0 – 62°C; 18:0 – 68°C) (Kachroo et al., 2007).
In plants, stearoyl-acyl carrier protein (ACP) desaturase (SAD, EC 220.127.116.11) is the nuclear-encoded, plastid-localized soluble desaturase that introduces the first Δ9 double bond into the saturated fatty acid resulting in the conversion of 18:0-ACP into 18:1n−9-ACP (Fox et al., 1993). As the dominant form of monounsaturated fatty acid exported from the plastid, 18:1n−9 can be further desaturated into polyunsaturated fatty acid derivatives, such as18:2n−6 and 18:3n−3, both of which serve as major components of cell membrane systems in the form of phospholipids (Ohlrogge and Browse, 1995). In this respect, given the fact that most plants lack other desaturases that utilize 18:0 as substrate, the activity of SAD is of particular interest because of its significant effects on the ratio of saturated fatty acids to unsaturated fatty acids (Lindqvist et al., 1996), and the effects on the fluidity and rigidity of membrane system and the relationship of this to the adaption of plants to various environmental conditions. For example, sunflower, flax, and castor grown in lower temperature conditions contained higher proportions of polyunsaturated fatty acid compared to those in high temperature growing regions (Green, 1986; Garces et al., 1992).
Given the functional importance of fatty acid saturation in plant development and industrial application, SAD genes from many plant species have been identified and characterized and a high correlations between the activities of SADs and levels of 18:0 and 18:1n−9 have been widely observed (Knutzon et al., 1992; Nishida et al., 1992; Slocombe et al., 1992; Gibson, 1993; Tong et al., 2006; Schluter et al., 2011; Shilman et al., 2011). For instance, a mutation of the Arabidopsis AtSSI2 gene, resulting in a single amino acid substitution resulted in an elevated level of 18:0 in Arabidopsis leaves (Lightner et al., 1994). Likewise, antisense-mediated reduction of SAD enzymatic activity in Brassica led to a dramatically increased level of 18:0 in the mature seed oil, from 2% up to 40% (Knutzon et al., 1992). Likewise, in soybean, a high 18:0 phenotype was strongly correlated with the activity of one of the SAD isozymes SACPD-C, an observation useful for screening of soybean varieties for unique seed oil profiles (Zhang et al., 2008).
In Arabidopsis, the family of acyl-ACP desaturases exhibits a high degree of amino acid sequence similarity (>70%) in spite of their minor functional differences (Kachroo et al., 2007). The three-dimensional crystal structure of homodimeric SAD protein from castor seed revealed that each SAD monomer comprises 11 conserved α-helices, which are crucial for ligand binding affinity, substrate chain-length selectivity, and double bond insertion position (Fox et al., 1993; Lindqvist et al., 1996; Cahoon et al., 1997). More impressively, replacements of several key amino acid residues in those conserved functional domains were sufficient to significantly alter the preferential mode of SAD. For example, replacement of five amino acids (T181A/F200A/N205S/T206L/A207G) of a castor Δ9 SAD resulted in a new specificity as a Δ6-16:0-ACP desaturase (Cahoon et al., 1997). Similarly, mutation of three amino acids (T117R/G188L/D280K) in castor Δ9 SAD enables it convert stearoyl-ACP to the allylic alcohol (E)-10-18:1-9-OH instead of 18:1n−9 (Whittle et al., 2008). Therefore, it is feasible to precisely manipulate activity of SAD taking advantage of the in-depth knowledge of its functional mechanism and advanced genetic approaches.
Interestingly, activities of SADs and levels of 18:1 are also involved in regulation of the defense response in many plant species via the fatty acid (FA)-derived signaling pathway (Kachroo and Kachroo, 2009). In Arabidopsis, decreased levels of 18:1n−9 in the ssi2 mutant induced the SA mediated defense pathway, resulting in constitutive expression of PR genes, activation of the hypersensitive response and enhancement of broad-spectrum resistance to bacterial and oomycete pathogens (Shah et al., 2001; Kachroo et al., 2003, 2004). Similar defense responses were also observed in soybean and rice (Kachroo et al., 2008; Jiang et al., 2009), showing that silenced SAD activities were sufficient to induce the same type of defense-signaling pathway and enhance resistance to multiple pathogens. Together, it appears that 18:1 derived defense signaling pathway is conserved among many plant species. We recently demonstrated reduction of 18:1n−9 levels in cacao leaves, induced by application of exogenous glycerol, can induce a hypersensitive-like response and enhance the resistance against the oomycete pathogen Phytophthora capsici (Zhang et al., 2014b).
The cacao genome was previously shown to contain eight putative SAD genes (Argout et al., 2011). In this study, we examined the gene family in detail, and explored the expression patterns of each SAD gene in various cacao tissues. A single gene primarily involved in the synthesis of 18:1 pools in developing cacao seeds was identified and functionally characterized in Arabidopsis ssi2 mutant. This information can be used to develop biomarkers for screening and breeding of new cacao varieties with novel fatty acid compositions of cocoa butter.
Material and Methods
Gene Identification, Phylogenetic Analysis, and Protein Prediction
SAD isoforms in Theobroma cacao were identified by blastp (Altschul et al., 1990) using full-length amino acid sequence of Arabidopsis AtSSI2 as the query (E-value cut-off 1e−5). Multiple protein sequence alignment was performed by MUSCLE (Edgar, 2004). The phylogenetic tree was constructed by MEGA4.1 using neighbor-joining algorithm with Poisson correction model and pairwise deletion (Tamura et al., 2007). Bootstrap values represent 1000 replicates. The phylogenetic tree was rooted using the amino acid sequence of PpSAD from Physcomitrella patens. Molecular mass and isoelectric point of SAD isozymes were predicted on ExPASy server (http://web.expasy.org/compute_pi/) (Gasteiger et al., 2005). Transient signal peptides were predicted using ChloroP 1.1 server (Emanuelsson et al., 1999) and TargetP 1.1 server (Emanuelsson et al., 2000).
Genbank accessions of the genes used in phylogenetic analysis:
Fatty Acid Profiling by GC-MS
Plant tissue from glycerol treated and control leaves (four biological replicates) were ground in liquid nitrogen and fatty acid methyl esters (FAME) were prepared from each sample using approximately 30 mg of tissue per sample. Briefly, 1 ml of a MeOH/fuming HCl/Dichloromethane (10:1:1 v/v) solution was added to each tissue sample and incubated without shaking at 80°C for 2 h. Fatty acid methyl esters were re-extracted in 1 ml buffer H2O/Hexane/Dichloromethane (5/4/1, v/v) with vortexing for 1 min. The hexane (upper phase) was separated by centrifugation at 1500 g for 5 min, transferred to glass GC vials (Agilent) and evaporated to dryness under vacuum. The FAMEs were then dissolved in 500 μl hexane for GC-MS analysis. Pentadecanoic acid (C15:0) (Sigma, Cat. P6125) was used as the internal standard added prior to the extraction and methyl nonadecanoate (C19:0-methyl ester) (Sigma, Cat. N5377) was used as the spike control, added into the sample prior to the GC injection. Samples were analyzed on an Agilent 6890N gas chromatograph coupled to a Waters GCT time of flight mass spectrometer. Mass spectra were acquired in electron ionization mode (70 eV) from 45 to 500 Da at a rate of 1 scan/s. The samples were separated on an Omegawax® 250 Capillary GC column (30 m × 0.25 mm 0.25 uM phase thickness, Sigma, Cat. 24136) using helium at a constant flow of 1.0 ml/min. The initial oven temperature was 100°C held for 1 min then increased at 15°C/min to a temperature of 150°C, and then increased at 4°C/min to a final temperature of 280°C. Samples (1 ul) were injected onto the column using a split/splitless injector maintained at 240°C with a split ratio of 50/1.
RNA Extraction and RT-qPCR Analysis of Gene Expression
Plant tissues were first ground in liquid nitrogen. Three biological replicates for each tissue were subjected for the gene expression analysis. Total RNA was extracted using Plant RNA Purification Reagent (Life Technologies, Cat. 12322-012, following the manufactures instructions). The concentration of RNA was measured using a Nanodrop 2000c (Thermo Scientific). Five hundred nanograms of RNA was further treated with RQ1 RNase-free DNase (Promega, Cat. M6101) at 37°C for 30 min to remove potential genomic DNA contamination (following the manufacturer's protocol). The treated RNA was reverse-transcribed by M-MuLV Reverse Transcriptase (New England Biolabs) with oligo-dT15 primers to obtain cDNA. RT-qPCR was performed in total reaction volume of 10 μl containing 4 μl diluted-cDNA (1:50), 5 μl SYBR Green PCR Master Mix (Takara), 0.2 μl Rox, and 0.4 μl each 5 μM primers. Each reaction was performed in duplicates using Roche Applied Biosystem Step One Plus Realtime PCR System under the following program: 15 min at 94°C, 40 cycle of 15 s at 94°C, 20 s at 60°C, and 40 s at 72°C. The specificity of the primer pair was verified by PCR visualized on a 2% agarose gel and analysis of the qPCR dissociation curve. A tubulin gene (Tc06g000360, TcTUB1, TcTUB1-5′: GGAGGAGTCTCTATAAGCTTGCAGTTGG and TcTUB1-3′: ACATAAGCATAGCCAGCTAGAGCCAG) and a gene encoding an acyl-carrier protein (Tc01g039970, TcACP1, TcACP1-5′: GGAAAGCAAGGGTGTCTCGTTGAA and TcACP1-3′: GCGAGTTGAAATCTGCTGTTGTTTGG) were used as reference genes.
Characterization of Arabidopsis Mutants and Arabidopsis Transformation
All Arabidopsis plants were grown in a Conviron growth chamber at 22°C with 16 h light/8 h dark cycle. Arabidopsis ssi2 mutant was kindly provided by Dr. Kachroo (University of Kentucky), and homozygous fab2 mutant (SALK_036854) was obtained from the Arabidopsis Biological Resource Center at Ohio State University (https://abrc.osu.edu/). The base pair mutation of AtSSI2 in ssi2 mutant was confirmed by sequencing the PCR product amplified by the following primers (LP1: TGAAGAAACCATTTACGCCAC; p3: CGTGTTGACATGAGGCAGATCG). The presence and homozygosity of T-DNA insertion in the fab2 mutant line were confirmed by genotyping using the following primers (LB: CTTTGACGTTGGAGTCCAC; Up1F: TGAAACAGGTGCTAGTCCTACTTCA; Dn1R: CACCTGAAAGCCCGGTTAAGTC) (Schluter et al., 2011). Expression levels of AtSSI2 in Col-0, ssi2, and fab2 were examined by semi-quantitative RT-PCR using the following intron spanning primers (AtSSI2-5′: GGCCCCAAGGAGGTTGAGAG; AtSSI2-3′: ATCTGGAATGGATCCGCGGAC).
The full length coding sequence of TcSAD1 was first amplified from leaf cDNA with the following primers (TcSAD1-5′-XbaI: GCTCTAGAATGGCTCTGAAATTGAATCCCAT and TcSAD1-3′-XbaI: TCTAGACTAGAGCTTCACTTCTCTATCAAAAATCCAAC). The resulting fragment was cloned into intermediate vector pE2113 using XbaI, driven under enhanced E12-Ω CaMV 35S promoter (Mitsuhara et al., 1996). The whole gene cassette was then cloned into the binary plasmid pCAMBIA-1300 (Hajdukiewicz et al., 1994) using EcoRI and HindIII sites, resulting in plasmid pGZ14.1021 (GenBank Accession KP100426) and was further transformed into A. tumefaciens strain AGL1 by electroporation. Col-0, ssi2, and fab2 Arabidopsis genotypes were transformed using the floral dip method (Clough and Bent, 1998). Transformants were selected on hygromycin B containing MS media (Harrison et al., 2006). Expression levels of TcSAD1 in transgenic plants were confirmed by RT-PCR using the following primers: (TcSAD1-RT-5′: GCTCTGAAATTGAATCCCATCACTTCTCAA; TcSAD1-RT-3′: TGGCTCCAAAGAGGTTGAGAATGTC). Arabidopsis Ubiquitin (AtUbi-5′: ACCGGCAAGACCATCACTCT; AtUbi-3′: AGGCCTCAACTGGTTGCTGT) was used as a reference gene.
Results and Discussion
The T. cacao Contains Eight Highly Conserved SAD Isoforms in the Genome
To identify the orthologous genes encoding SAD(s) in cacao genome, the full-length amino acid sequence of AtSSI2 (At2g43710) was blasted against the predicted cacao proteome of the Belizean Criollo genotype (B97-61/B2) (http://cocoagendb.cirad.fr/ Argout et al., 2011) using blastp algorithm with E-value cut-off 1e−5 (Altschul et al., 1990). Eight putative SAD genes were identified in cacao genome with reliable sequencing data and predicted gene structures (Table 1). Similar results were found with Blast analysis of the predicted cacao proteome of the genotype Matina1-6 (http://www.cacaogenomedb.org; Motamayor et al., 2013). The reciprocal best hits were further identified from both cacao genomes through blast and phylogenetic analysis (Supplemental File 1). Interestingly, of these eight SAD isoforms, four genes are located on chromosomes 4, and three of them (Tc04g017510, Tc04g017520, Tc04g017540) are clustered except for one predicted transposable element gene in between, which is similar to in the Arabidopsis genome where three out of seven SAD isoforms (At3g02610, At3g02620, and At3g02630) are located in tandem on chromosome 3 (Kachroo et al., 2007). Multiple amino acid sequence comparison and phylogenetic analysis including SAD isoforms from cacao, castor, and Arabidopsis revealed that Tc04g017510 (designated as TcSAD1) is closely clustered with RcSAD1 and AtSSI2 (Figure 1), both of which have significant impacts on oleic acid contents and fatty acid profiles in Arabidopsis (Lightner et al., 1994) and castor (Lindqvist et al., 1996), respectively, suggesting a potential significant role of TcSAD1 in cacao seed oil biosynthesis. Moreover, Tc05g012840 (designated as TcSAD2), Tc04g017520 (designated as TcSAD3), and Tc04g017540 (designated as TcSAD4) also exhibit high sequence similarities to RcSAD1, implying that they may also contribute to the activity of SAD in cacao.
Figure 1. Neighbor-joining bootstrap tree phylogeny based on full-length amino acid sequences of stearoyl-acyl carrier protein-desaturase (SAD) gene families in Arabidopsis (AtSAD), castor (RcSAD), and cacao (TcSAD). The scale bar represents 0.05 estimated substitutions per residue and values next to nodes indicated bootstrap values from 1000 replicates. The phylogenetic tree was rooted by using the amino acid sequence of PpSAD from Physcomitrella patens.
The metabolic pathways in which SADs participate in are well-known to occur in the lumen of plastids (Cahoon et al., 1997). Thus, we performed bioinformatics analysis of transit peptides (TP) in the N-terminal sequences of the predicted cacao SAD proteins using the TargetP (Emanuelsson et al., 2000) and ChloroP (Emanuelsson et al., 1999) algorithms. AtSSI2 (At2g43710) in Arabidopsis and RcSAD1 (XP_002531889) in castor were also included in this analysis since the structures and the functions of these two proteins have been well-studied (Lindqvist et al., 1996; Cahoon et al., 1997; Kachroo et al., 2007; Whittle et al., 2008). High-confidence chloroplast transient peptides were successfully detected in AtSSI2 and RcSAD1, consistent with previous analysis (Kachroo et al., 2007; Schluter et al., 2011). Of eight cacao SAD isoforms, chloroplast transient peptides were identified in TcSAD1, TcSAD2, TcSAD5, and TcSAD7. A mitochondrial targeting peptide was predicted in TcSAD3 and a signal peptide related with the secretory pathway was detected in TcSAD4, suggesting that these isoforms may exhibit distinct functions in mitochondria. No transient peptides were detected in TcSAD6 and TcSAD8 (Table 2).
SAD Isozymes in T. cacao Share Conserved Overall Protein Structures with Noteworthy Differences
The three-dimensional protein structure of RcSAD1 from castor seed has been extensively studied in detail as a model of Δ9 desaturase in plants (Schneider et al., 1992; Fox et al., 1993; Lindqvist et al., 1996). To provide further insights into potential catalytic activities of SAD isozymes in T. cacao, a multiple amino acid sequence alignment was conducted including RcSAD1, AtSSI2, and eight cacao SADs, and the determinative amino acid residues that are highly associated with substrate binding activities and chain length specificities of the SADs were evaluated (Figure 2). As mentioned above, a high degree of variation was observed within putative transient signal peptides in N-terminal sequences. Consistent with crystal structure of RcSAD1, cacao SADs also consist of 11 highly conserved α-helices, except that putative TcSAD8 lacks of the region from α3b- to α5-helices, which might greatly compromise its normal functionality. Notably, a further detailed investigation of the determinative residues revealed that several prominent divergences exist within cacao SADs, implying that they might have distinct substrate preferences. For example, amino acid residues at positions of 117, 118, 189, and 206 (positions assigned according to castor RcSAD1 amino acid sequence, indicated in Figure 2) were reported to be crucial for the function of RcSAD1 as a Δ9-18:0-ACP desaturase, and replacements of these residues converted RcSAD1 into an enzyme that either functioned as a Δ6-16:0-ACP desaturase or converted 18:0-ACP into the allylic alcohol trans-isomer (E)-10-18:1-9-OH (Cahoon et al., 1997; Whittle et al., 2008). In this respect, TcSAD3 and TcSAD4 share the same amino acid residues at those positions which differ from RcSAD1, AtSSI2, and other cacao SADs, suggesting that TcSAD3 and TcSAD4 may exert different functions from the other cacao SADs. Likewise, at positions of 179 and 181, which are crucial for the substrate specificity of RcSAD1 (Cahoon et al., 1997), TcSAD5, TcSAD6, and TcSAD8 contain unique varied amino acid residues compared to RcSAD1, AtSSI2, and the other cacao SADs.
Figure 2. Multiple amino acid sequence alignment of TcSADs from cacao with SAD from castor bean (RcSAD1) and Arabidopsis (AtSSI2). Secondary protein structures of SAD were compared, deduced, and annotated according to the crystal structure of RcSAD1. Shadings were performed by Genedoc software (Nicholas et al., 1997) in a conservation mode. Residues in black represent 100% conservation of all sequences; residues in dark gray represent >75% identity of all sequences; residues in light gray >50% identity of all sequences. Functional determinative amino acids were marked in rectangles and numbered according to the start codon of RcSAD1.
Identification and Validation of Housekeeping Genes Across T. cacao Tissues
In order to accurately characterize the expression pattern of SAD isoforms, suitable reference genes for qRT-PCR normalization were first identified and validated in various cacao tissues including: leaves at developmental stages A, C, and E (defined in Mejia et al., 2012), unopened flowers, open flowers, roots, and zygotic seeds at 14, 16, 18, and 20 weeks after pollination (WAP). Ten candidate genes primarily involved in cellular structures and central metabolisms were selected and evaluated including: tubulin (TUB), actin (ACT), acyl carrier protein (ACP), elongation factor (EF1α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and malate dehydrogenase (MDH). All the genes were successfully amplified in all the tested cacao tissues using gene specific primers that preferentially bind near the 3′ end of the cDNA (Supplemental Table 1). The specificities of the primers were confirmed by melting curve analysis (Supplemental File 2). The reference gene stability analysis was conducted applying two statistical approaches: geNorm (Vandesompele et al., 2002), which calculates M-values to suggest the mean expression stability of given genes, and NormFinder (Andersen et al., 2004), which estimates not only the overall variation of candidate reference genes but also the variation within biological sample subgroups and experimental conditions. According to both geNorm and NormFinder algorithms, TUB1 (Tc06g000360) has the lowest average expression stability values (Figures 3A,B), suggesting TUB1 is the most stable gene in the tested gene set across all these cacao tissues. Therefore, TUB1 was selected as the reference gene to characterize the gene expression pattern of SAD isoforms in various cacao tissues. Additionally, given central roles of ACPs in transporting and stabilizing the growing fatty acid carbon chains throughout the fatty acid biosynthetic pathway (Crosby and Crump, 2012), ACP1 was selected as the reference gene to normalize gene expression level during cacao seed development since ACP1 was also recognized as the second most stable reference gene from NormFinder analysis (Figure 3B).
Figure 3. Stability values of candidate reference genes. (A) Average expression stability M-values of expression levels of candidate reference genes in various cacao tissues calculated by geNorm software. (B) Stability value of candidate reference genes in various cacao tissues calculated by NormFinder software. Tissues include: leaves stage A, C, and E [early, mid, and late development, see (LA, LC, LE), developing flowers (UF) and mature flowers (OF), root (Rt), and zygotic cotyledons from seeds collected at 14, 16, 18, and 20 weeks after pollination (14, 16, 18, and 20 Cot)].
Tissue Specific Expression Pattern of SAD Isoforms in T. cacao
To further characterize functions of the cacao SAD family, the expression profiles were examined by qRT-PCR in various cacao tissues including: leaves at developmental stages A, C, and E (defined in Mejia et al., 2012), unopened flowers, open flowers, roots, and zygotic seeds at 14, 16, 18, and 20 WAP. Since SAD isoforms are highly conserved in their functional domains (α-helices and β-sheets), isoform specific qRT-PCR primers were designed to preferably target at predicted 5′ and 3′ untranslated regions, and signal transient peptide regions, all of which exhibited low sequence similarities among SAD isoforms in cacao (Supplemental Table 2). The specificities of each set of primers were determined by examining amplification products on the agarose gel and further confirmed by melting point analysis (Supplemental File 3). Expression levels of SAD isoforms were normalized to TcTUB1 (Table 3). As shown in Figure 4, transcripts of TcSAD1 and TcSAD2 were detected in all the examined tissues, in which the expression levels of TcSAD1 were higher than those of TcSAD2. Interestingly, TcSAD3 and TcSAD4 were predominantly expressed in unopened and open flowers, where the expression levels of TcSAD3 and TcSAD4 were significantly higher than any other SAD isoforms; however, the transcripts of both were barely detectable in all the other tissues, indicating that they are flower specific genes. The expression of TcSAD5 was primarily detected in roots; however, at levels lower than those of TcSAD1, TcSAD2, and TcSAD7 in roots. TcSAD7 constitutively expresses in all examined tissues with higher expression in roots and zygotic seeds at all developmental stage (Figure 4). Overall, on the basis of gene expression data, both TcSAD1 and TcSAD7 contribute to the major proportions to the total SAD transcripts in developing zygotic embryos, with TcSAD7 being higher expressed over all the stages except in zygotic seeds at 20 WAP, implying that the activities of both TcSAD1 and TcSAD7 are possibly involved in the synthesis and accumulation of 18:1 in cacao embryos (Figure 4). Notably, the transcripts of TcSAD6 and TcSAD8 were barely detected in any of the tissues analyzed, suggesting that they might be pseudogenes or expressed in a tissue, stage of development or induction condition not tested.
Figure 4. TcSAD expression patterns in various cacao tissues. Bar charts illustrate the relative expression levels of TcSAD isoforms in various cacao tissues normalized to that of TcTUB1 measured by qRT-PCR. The same set of cacao tissues as indicated in Figure 3 was included in this analysis. Three biological samples for each tissue were included in the analysis.
Profoundly, TcSAD3 and TcSAD4 appear to be exclusively and highly expressed in floral tissues. These two proteins are unique within the cacao SAD gene family in containing predicted mictochondrial targeting and secretory peptides respectively (Table 2). Interestingly, a previous study in orchid revealed that the activities of SAD isoforms were involved in the biosynthesis of alkene species with different double bond positions, which result in the major odor differences among orchid species attracting differential pollinators (Schluter et al., 2011). Moreover, a more recent study in tobacco showed that the normal function of one isoform of stearoyl-ACP desaturase was critical to maintain membrane lipid composition during ovule development and thus affected female fertility (Zhang et al., 2014a). Therefore, even though functions of TcSAD3 and TcSAD4 in cacao floral tissues remain largely unknown, it is reasonable to speculate that their activity might be associated with the synthesis of fatty acid derived metabolites that are vital to the pollination process and the development of cacao reproductive tissues.
Oleic Acid Accumulation and SAD Isoform Expression during Cacao Seed Development
Notably, the maturation of cacao seeds initiates from the axial ends and gradually toward the cotyledonary ends. Seeds during this rapid transition are ideal to evaluate the correlation between oleic acid contents and the activities of TcSADs since biological variations are greatly minimized. To explore the role of each TcSAD isoform in seed lipid biosynthesis, maturing cacao seeds were isolated from fruit at 12 WAP (Figure 5A), and the entire cacao seed was transversely dissected into six segments along the developmental gradient, with the least developed segment (segment 1, color coded as dark green, Figure 5A) at the cotyledonary end and the most developed segment (segment 6, color coded as dark red, Figure 5A) at the axial end. The analysis of fatty acid composition in each segment indicated that the percentage of 18:1n−9 gradually increased from 12% in segment 1 to 35% in segment 6 (Figure 5B), which is consistent with previous reports of fatty acid composition changes during cacao seed development over time (Patel et al., 1994). Thereafter, the expression of two major TcSAD isoforms in seeds (TcSAD1 and TcSAD7, Figure 4) was examined in each segment by qRT-PCR, and the correlation between the expression levels of TcSADs and fatty acid composition was evaluated by linear regression (Figure 5C). Expression levels of SAD isoforms were normalized to TcACP1. In general, the levels of SAD gene expression shown in Figure 5C suggested that TcSAD1 expression (left panel) was highly correlated with gradual alterations of fatty acid composition during development. However, expression of TcSAD7 did not correlate well with fatty acid composition changes (right panel). Specifically, TcSAD1 expression was positively correlated with the molar percentage of 18:1n−9 (r = 0.799) and 18:0 (r = 0.522), and negatively correlated with the molar percentage of 18:2n−6 (r = −0.638); however, poor correlations were observed between the expression of TcSAD7 and these fatty acid species. Additionally, the expression of TcSAD1 was also positively correlated with the level of total molar percentage of unsaturated fatty acids (18:1n−9 + 18:2n−6 + 18:3n−3) (r = 0.438) and the ratio of unsaturated fatty acid (as above) over saturated fatty acid (16:0 + 18:0) (r = 0.464). Again, the expression of TcSAD7 was poorly correlated with these two aspects.
Figure 5. TcSAD relative expression and fatty acid composition in maturing cacao seeds. (A) (Upper) Mid-stage cacao seed sliced longitudinally, within ovule and exocarp tissue (pod shell) 12 weeks after pollination. The embryo root axial meristem is to the right, expanding cotyledonary tissues to the left. (Lower) Six radial serial thick sections dissected from a single cacao seed oriented as in the top image, then rotated 90 degrees to display the cross sections. Each section is sequentially numbered from left to right and color coded by green, light green, yellow, orange, light red, red, respectively, from youngest cotyledonary end (left) to the most mature axial end (right). Scale bar = 1 cm. (B) Molar percentages of fatty acid compositions in each seed section. (n = 3, mean ± SE). (C) The correlation of TcSAD1 and TcSAD7 gene expression and molar percentage of fatty acid composition in various cross sections of developing cacao seeds. Left column, TcSAD1, right column, TcSAD7. Color-coded data points represent different seed sections as labeled in (A). Three biological replicates from each seed section were analyzed. Each independent measurement is shown as a separate data point. Correlation coefficient r was calculated based on the linear regression model by Minitab (Ryan et al., 2004).
Taken together, the evidence suggests that TcSAD1 is the dominant isoform for the synthesis of seed storage fatty acids in cacao since the increasing expression levels of TcSAD1 were highly correlated with the gradual changes of fatty acid compositions. It further suggests that TcSAD7 may play a minor role in this process and could play additional roles in seed development unrelated to storage lipid accumulation. Regarding this, a recent characterization of AtSAD6, a TcSAD7 homolog in Arabidopsis, revealed that AtSAD6 may be involved in the fatty acid desaturation and lipid metabolism in crown galls under hypoxia and drought stress conditions (Klinkenberg et al., 2014).
TcSAD1 is Functional Equivalent to AtSSI2
To further investigate in vivo TcSAD1 and TcSAD7 protein functions, TcSAD1 and TcSAD7 were overexpressed in Arabidopsis ssi2 mutant background, driven by Cauliflower mosaic virus 35S promoter. Owing to the poor transformation rate and mortality rate of the ssi2 mutant (Schluter et al., 2011), only six transgenic lines overexpressing TcSAD1 were identified; however, no transgenic plant overexpressing TcSAD7 was obtained from three Arabidopsis transformation attempts. Notably, two out of six 35S::TcSAD1 (ssi2) transgenic lines were morphologically indistinguishable from wild type Col-0 plants (Figure 6A), suggesting that the function of TcSAD1 successfully rescued the dwarf phenotype of ssi2 mutant. RT-PCR analysis confirmed that transcripts of TcSAD1 were only detected in the transgenic lines, but not in wild type or ssi2 mutant plants (Figure 6B). Fatty acid composition analysis on wild type Col-0, ssi2 mutants and the 35S::TcSAD1 (ssi2) transgenic lines revealed that overexpression of TcSAD1 restored the wild type-like level of 18:0 (Figure 6C). Another dramatic phenotype of the Arabidopsis ssi2 is the constitutive expression of PR genes and spontaneous cell death (Shah et al., 2001). In this respect, 35S::TcSAD1 (ssi2) plants exhibited no visible cell death on the leaves and contained only basal level of PR-1 expression (Figure 6B). Taken together, we concluded that the overexpression of TcSAD1 was sufficient to functional complement the reduced activity of SSI2 in the ssi2 Arabidopsis mutant, demonstrating that TcSAD1 is a functional ortholog of AtSSI2.
Figure 6. Transgenic complementation of the Arabidopsis ssi2 mutation with the cacao TcSAD1 gene. Comparison of morphological, molecular, and fatty acid composition phenotypes of ssi2 and 35S::TcSAD1 (ssi2) plants. (A) Morphological comparison of wild-type Col-0, ssi2, and 35S::TcSAD1 (ssi2) at 4-week-old plants. Overexpression of TcSAD1 in Arabidopsis ssi2 successfully rescued dwarf phenotype and eliminated hypersensitive responses. Scale bar = 1 cm. (B) Overexpression of TcSAD1 eliminated constitutive overexpression of AtPR1 (At2g14610) in ssi2, but did not affect the level of AtSSI2. Arabidopsis Ubiquitin was used as the reference gene in RT-PCR analysis. (C) Fatty acid compositions of wild-type Col-0, ssi2, and 35S::TcSAD1 (ssi2) (n = 6 for Col-0 and ssi2, n = 2 for 35S::TcSAD1, mean ± SD).
In conclusion, the eight TcSAD isoforms in cacao genome share high amino acid sequence conservation but they also have specific differences in key determinative amino acid residues and distinct tissue specific expression patterns. Among them, TcSAD1 is phylogentically most closely clustered with RcSAD1 and AtSAD1, and is functional equivalent to AtSAD1. Notably, the expression level of TcSAD1 was also highly positively correlated with the level of 18:1n−9 and therefore affected final fatty acid profiles in maturing cacao seeds, which makes it as an ideal biomarker to screen for desirable fatty acid compositions of cocoa butter. TcSAD1 successfully complimented the Arabidopsis ssi2 mutant. Together, our results strongly support the conclusion that TcSAD1 is the functional ortholog of AtSSI2, and plays a major role in determining the fatty acid composition of cacao seeds.
YZ conceived and performed all the experiments and drafted the manuscript. SM was involved in designing and directing the experiments, and revising the manuscript. MG contributed to the conception of the study, gave advice on experiments, and finalized the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Lena Sheaffer and Brian Rutkowski for the technical assistance in maintenance of our greenhouse. This work was supported in part by The Pennsylvania State University, College of Agricultural Sciences, The Huck Institutes of Life Sciences, and the American Research Institute Penn State Endowed Program in the Molecular Biology of Cacao.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fpls.2015.00239/abstract
Supplemental File 1. Neighbor-joining bootstrap tree phylogeny based on full-length amino acid sequences of stearoyl-acyl carrier protein-desaturase (SAD) gene families in Criollo and Matina1-6 genotypes of cacao.
Supplemental File 2. Melt curve analysis of 10 candidate reference genes in cacao.
Supplemental File 3: Melt curve analysis of SAD isoforms in cacao.
Andersen, C. L., Jensen, J. L., and Ørntoft, T.F. (2004). Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 64, 5245–5250. doi: 10.1158/0008-5472.CAN-04-0496
Cahoon, E. B., Lindqvist, Y., Schneider, G., and Shanklin, J. (1997). Redesign of soluble fatty acid desaturases from plants for altered substrate specificity and double bond position. Proc. Natl. Acad. Sci. U.S.A. 94, 4872–4877. doi: 10.1073/pnas.94.10.4872
Emanuelsson, O., Nielsen, H., Brunak, S., and Von Heijne, G. (2000). Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J. Mol. Biol. 300, 1005–1016. doi: 10.1006/jmbi.2000.3903
Emanuelsson, O., Nielsen, H., and Von Heijne, G. (1999). ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci. 8, 978–984. doi: 10.1110/ps.8.5.978
Fox, B. G., Shanklin, J., Somerville, C., and Munck, E. (1993). Stearoyl-acyl carrier protein delta 9 desaturase from Ricinus communis is a diiron-oxo protein. Proc. Natl. Acad. Sci. U.S.A. 90, 2486–2490. doi: 10.1073/pnas.90.6.2486
Gasteiger, E., Hoogland, C., Gattiker, A., Duvaud, S. E., Wilkins, M., Appel, R., et al. (2005). “Protein identification and analysis tools on the ExPASy server,” in The Proteomics Protocols Handbook, ed J. Walker (Totawa, NJ: Humana Press), 571–607.
Harrison, S. J., Mott, E. K., Parsley, K., Aspinall, S., Gray, J. C., and Cottage, A. (2006). A rapid and robust method of identifying transformed Arabidopsis thaliana seedlings following floral dip transformation. Plant Methods 2:19. doi: 10.1186/1746-4811-2-19
Jiang, C. J., Shimono, M., Maeda, S., Inoue, H., Mori, M., Hasegawa, M., et al. (2009). Suppression of the rice fatty-acid desaturase gene OsSSI2 enhances resistance to blast and leaf blight diseases in rice. Mol. Plant Microbe. Interact. 22, 820–829. doi: 10.1094/MPMI-22-7-0820
Kachroo, A., Fu, D. Q., Havens, W., Navarre, D., Kachroo, P., and Ghabrial, S. A. (2008). An oleic acid-mediated pathway induces constitutive defense signaling and enhanced resistance to multiple pathogens in soybean. Mol. Plant Microbe Interact. 21, 564–575. doi: 10.1094/MPMI-21-5-0564
Kachroo, A., Lapchyk, L., Fukushige, H., Hildebrand, D., Klessig, D., and Kachroo, P. (2003). Plastidial fatty acid signaling modulates salicylic acid- and jasmonic acid-mediated defense pathways in the Arabidopsis ssi2 mutant. Plant Cell 15, 2952–2965. doi: 10.1105/tpc.017301
Kachroo, A., Shanklin, J., Whittle, E., Lapchyk, L., Hildebrand, D., and Kachroo, P. (2007). The Arabidopsis stearoyl-acyl carrier protein-desaturase family and the contribution of leaf isoforms to oleic acid synthesis. Plant Mol. Biol. 63, 257–271. doi: 10.1007/s11103-006-9086-y
Kachroo, A., Venugopal, S. C., Lapchyk, L., Falcone, D., Hildebrand, D., and Kachroo, P. (2004). Oleic acid levels regulated by glycerolipid metabolism modulate defense gene expression in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 101, 5152–5157. doi: 10.1073/pnas.0401315101
Klinkenberg, J., Faist, H., Saupe, S., Lambertz, S., Krischke, M., Stingl, N., et al. (2014). Two fatty acid desaturases, STEAROYL-ACYL CARRIER PROTEIN Δ9-DESATURASE6 and FATTY ACID DESATURASE3, are involved in drought and hypoxia stress signaling in Arabidopsis crown galls. Plant Physiol. 164, 570–583. doi: 10.1104/pp.113.230326
Knutzon, D. S., Thompson, G. A., Radke, S. E., Johnson, W. B., Knauf, V. C., and Kridl, J. C. (1992). Modification of Brassica seed oil by antisense expression of a stearoyl-acyl carrier protein desaturase gene. Proc. Natl. Acad. Sci. U.S.A. 89, 2624–2628. doi: 10.1073/pnas.89.7.2624
Liendo, R., Padilla, F. C., and Quintana, A. (1997). Characterization of cocoa butter extracted from Criollo cultivars of Theobroma cacao L. Food Res. Int. 30, 727–731. doi: 10.1016/S0963-9969(98)00025-8
Lindqvist, Y., Huang, W., Schneider, G., and Shanklin, J. (1996). Crystal structure of delta9 stearoyl-acyl carrier protein desaturase from castor seed and its relationship to other di-iron proteins. EMBO J. 15, 4081–4092.
Mejia, L. C., Guiltinan, M. J., Shi, Z., Landherr, L., and Maximova, S. N. (2012). “Expression of designed antimicrobial peptides in Theobroma cacao L. trees reduces leaf necrosis caused by Phytophthora spp.,” in Small Wonders: Peptides for Disease Control, Vol. 1095 (Washington, DC: American Chemical Society), 379–395. doi: 10.1021/bk-2012-1095.ch018
Mitsuhara, I., Ugaki, M., Hirochika, H., Ohshima, M., Murakami, T., Gotoh, Y., et al. (1996). Efficient promoter cassettes for enhanced expression of foreign genes in dicotyledonous and monocotyledonous plants. Plant Cell Physiol. 37, 49–59. doi: 10.1093/oxfordjournals.pcp.a028913
Motamayor, J. C., Mockaitis, K., Schmutz, J., Haiminen, N., Donald Iii, L., Cornejo, O., et al. (2013). The genome sequence of the most widely cultivated cacao type and its use to identify candidate genes regulating pod color. Genome Biol. 14:r53. doi: 10.1186/gb-2013-14-6-r53
Nishida, I., Beppu, T., Matsuo, T., and Murata, N. (1992). Nucleotide sequence of a cDNA clone encoding a precursor to stearoyl-(acyl-carrier-protein) desaturase from spinach, Spinacia oleracea. Plant Mol. Biol. 19, 711–713. doi: 10.1007/BF00026799
Patel, V. K., Shanklin, J., and Furtek, D. B. (1994). Changes in fatty-acid composition and stearoyl-acyl carrier protein desaturase expression in developing Theobroma-cacao L embryos. Planta 193, 83–88. doi: 10.1007/BF00191610
Ryan, T. A., Joiner, B. L., and Ryan, B. F. (2004). “Minitab™,” in Encyclopedia of Statistical Sciences, eds S. Kotz, C. B. Read, N. Balakrishnan and B. Vidakovic (Hoboken, NJ: John Wiley & Sons, Inc.). doi: 10.1002/0471667196.ess1647.pub2
Schluter, P. M., Xu, S., Gagliardini, V., Whittle, E., Shanklin, J., Grossniklaus, U., et al. (2011). Stearoyl-acyl carrier protein desaturases are associated with floral isolation in sexually deceptive orchids. Proc. Natl. Acad. Sci. U.S.A. 108, 5696–5701. doi: 10.1073/pnas.1013313108
Schneider, G., Lindqvist, Y., Shanklin, J., and Somerville, C. (1992). Preliminary crystallographic data for stearoyl-acyl carrier protein desaturase from castor seed. J. Mol. Biol. 225, 561–564. doi: 10.1016/0022-2836(92)90941-C
Shah, J., Kachroo, P., Nandi, A., and Klessig, D. F. (2001). A recessive mutation in the Arabidopsis SSI2 gene confers SA- and NPR1-independent expression of PR genes and resistance against bacterial and oomycete pathogens. Plant J. 25, 563–574. doi: 10.1046/j.1365-313x.2001.00992.x
Shilman, F., Brand, Y., Brand, A., Hedvat, I., and Hovav, R. (2011). Identification and molecular characterization of homeologous delta 9-stearoyl acyl carrier protein desaturase 3 genes from the Allotetraploid Peanut (Arachis hypogaea). Plant Mol. Biol. Rep. 29, 232–241. doi: 10.1007/s11105-010-0226-9
Slocombe, S. P., Cummins, I., Jarvis, R. P., and Murphy, D. J. (1992). Nucleotide sequence and temporal regulation of a seed-specific Brassica napus cDNA encoding a stearoyl-acyl carrier protein (ACP) desaturase. Plant Mol. Biol. 20, 151–155. doi: 10.1007/BF00029157
Tong, L., Shu-Ming, P., Wu-Yuan, D., Dan-Wei, M., Ying, X., Meng, X., et al. (2006). Characterization of a new stearoyl-acyl carrier protein desaturase gene from Jatropha curcas. Biotechnol. Lett. 28, 657–662. doi: 10.1007/s10529-006-0034-3
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:research0034. doi: 10.1186/gb-2002-3-7-research0034
Whittle, E. J., Tremblay, A. E., Buist, P. H., and Shanklin, J. (2008). Revealing the catalytic potential of an acyl-ACP desaturase: tandem selective oxidation of saturated fatty acids. Proc. Natl. Acad. Sci. U.S.A. 105, 14738–14743. doi: 10.1073/pnas.0805645105
Zhang, J., Li, J., Garcia-Ruiz, H., Bates, P. D., Mirkov, T. E., and Wang, X. (2014a). A stearoyl-acyl carrier protein desaturase, NbSACPD-C, is critical for ovule development in Nicotiana benthamiana. Plant J. 80, 489–502. doi: 10.1111/tpj.12649
Zhang, P., Burton, J. W., Upchurch, R. G., Whittle, E., Shanklin, J., and Dewey, R. E. (2008). Mutations in a delta(9)-stearoyl-ACP-desaturase gene are associated with enhanced stearic acid levels in soybean seeds. Crop Sci. 48, 2305–2313. doi: 10.2135/cropsci2008.02.0084
Zhang, Y., Smith, P., Maximova, S. N., and Guiltinan, M. J. (2014b). Application of glycerol as a foliar spray activates the defence response and enhances disease resistance of Theobroma cacao. Mol. Plant Pathol. 16, 27–37. doi: 10.1111/mpp.12158
Keywords: Theobroma cacao, stearoyl-ACP-desaturase, oleic acid, ssi2, cocoa butter, fatty acid composition
Citation: Zhang Y, Maximova SN and Guiltinan MJ (2015) Characterization of a stearoyl-acyl carrier protein desaturase gene family from chocolate tree, Theobroma cacao L. Front. Plant Sci. 6:239. doi: 10.3389/fpls.2015.00239
Received: 29 January 2015; Accepted: 25 March 2015;
Published: 14 April 2015.
Edited by:Ramin Yadegari, The University of Arizona, USA
Reviewed by:Chaofu Lu, Montana State University, USA
Karen Hudson, Agricultural Research Service - United States Department of Agriculture, USA
Gary Norman Drews, University of Utah, USA
Copyright © 2015 Zhang, Maximova and Guiltinan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mark J. Guiltinan, Huck Institutes of the Life Sciences, Department of Plant Science, The Pennsylvania State University, University Park, 422 Life Sciences Building, PA 16802, USA email@example.com