Transcriptome profiling identifies ABA mediated regulatory changes towards storage filling in developing seeds of castor bean (Ricinus communis L.)
© Chandrasekaran et al.; licensee BioMed Central Ltd. 2014
Received: 12 May 2014
Accepted: 12 June 2014
Published: 30 June 2014
The potential biodiesel plant castor bean (Ricinus communis) has been in the limelight for bioenergy research due to the availability of its genome which raises the bar for genome-wide studies claiming advances that impact the “genome-phenome challenge”. Here we report the application of phytohormone ABA as an exogenous factor for the improvement of storage reserve accumulation with a focus on the complex interaction of pathways associated with seed filling.
After the application of exogenous ABA treatments, we measured an increased ABA levels in the developing seeds cultured in vitro using the ELISA technique and quantified the content of major biomolecules (including total lipids, sugars and protein) in treated seeds. Exogenous ABA (10 μM) enhanced the accumulation of soluble sugar content (6.3%) followed by deposition of total lipid content (4.9 %). To elucidate the possible ABA signal transduction pathways towards overall seed filling, we studied the differential gene expression analysis using Illumina RNA-Sequencing technology, resulting in 2568 (1507-up/1061-down regulated) differentially expressed genes were identified. These genes were involved in sugar metabolism (such as glucose-6-phosphate, fructose 1,6 bis-phosphate, glycerol-3-phosphate, pyruvate kinase), lipid biosynthesis (such as ACS, ACBP, GPAT2, GPAT3, FAD2, FAD3, SAD1 and DGAT1), storage proteins synthesis (such as SGP1, zinc finger protein, RING H2 protein, nodulin 55 and cytochrome P450), and ABA biosynthesis (such as NCED1, NCED3 and beta carotene). Further, we confirmed the validation of RNA-Sequencing data by Semi-quantitative RT-PCR analysis.
Taken together, metabolite measurements supported by genes and pathway expression results indicated in this study provide new insights to understand the ABA signaling mechanism towards seed storage filling and also contribute useful information for facilitating oilseed crop functional genomics on an aim for utilizing castor bean agricultural and bioenergy use.
The plant hormone abscisic acid (ABA) is broadly involved in various stress-related responses and developmental regulation in higher plants . Many studies have demonstrated that ABA participates in regulating organ development such as flower, fruit, root and seed in diverse plants [2–4]. In particular, ABA regulates several important aspects of seed development such as embryogenesis , endosperm development  and the biosynthesis of storage reserves (lipids, proteins and starch) during seed filling . Endogenous ABA levels are closely associated to storage reserve accumulation in developing seeds [8, 9]. The pattern of endogenous ABA accumulation in developing seeds appeared to be varied in different plants; however, most cases observed showed that endogenous ABA accumulation was tightly associated with seed growth and storage material accumulation in developing seeds [4, 8]. On the other hand, exogenous ABA application enhanced the levels of TAG’s (nine fold increase) in somatic embryos of Picea glauca and rapeseed cultures [7, 10]. In addition, application of exogenous ABA increased the levels of total sugar content as well as reducing sugars in cultured seedling of Phaseolus vulgaris and in cell cultures of oilseed rape [11, 12]. A 19 kDa oil body protein encoding oleosin and a desaturase lipid gene were found to be up regulated, correlating to the physiological changes after exogenous ABA treatment in B. napus embryo cultures . In Arabidopsis, ABA responsive genes LEC2 (CCAAT family), FUS3, ABI3 and a maize VP1 (viviparous 1) were identified to be critical for the embryogenesis . In maize, the ABA responsive gene VP8 was demonstrated to have an essential role in coordinating embryo and endosperm development . Besides, ABI3 and FUS3 genes were found to be critical in mediating the accumulation of storage materials in Arabidopsis seeds [13, 14]. More than forty ABA response genes involved in embryogenesis, endosperm development and the accumulation of storage materials were identified in developing seeds of rice . Although diverse ABA responsive genes involved in seed development and storage material accumulation have been identified in different plants, the molecular basis for ABA dependent physiological responses in developing seeds is far from understanding.
Investigating the transcriptional changes of responsive genes to ABA signal by applying exogenous ABA as stimulant is of great help to identify critical ABA regulators in developing seeds. Many such efforts have been made and identified diverse genes involved in many aspects of seed development in different plants . In particular, high-throughout sequencing technology has recently become a powerful tool that allows to profile a genomic expression pattern and to measure modest changes among different samples at unprecedented perspectives [16, 17]. However limited information is known about the holistic transcriptional changes of responsive genes to ABA signals in developing seeds.
Castor bean (Ricinus communis, Euphorbiaceae) is one of the most important non-edible oilseed crops. Its seed is a highly specialized storage organ for accumulating lipids, proteins, and carbohydrates in endosperm . Its seed oils comprised of the unusual hydroxy-fatty acid (i.e. ricinoleic acid) is broadly used in industry and is considered as an ideal and unique feedstock for biodiesel production . Its seed storage protein ricin is extremely toxic and able to be reportedly used as a biochemical weapon or a specific immunotoxin for therapeutic purposes in different cancer treatments [20, 21]. Due to economical importance, castor bean has been widely accepted as an agricultural solution for all subtropical and tropical regions that addresses the need for commercial crops with low impute costs and at the same time provides traditional farming with a viable income from current non- productive lands. With increasing demand for production of castor bean seed oils in many countries, breeding and improvement of varieties are drawing great attention from breeders. Further efforts should be made to elucidate the molecular mechanism underlying the regulation of growth and development . Study on the molecular basis of ABA physiological response in developing seeds is of significance for understanding the mechanism controlling seed development and crop improvement in castor bean.
In this study, we inspected the effect of exogenous ABA signaling on the storage material accumulation and preformed high-throughput sequencing to identify the ABA responsive genes in castor bean seed filling. Results obtained here provide new insights into the understanding of molecular mechanisms underlying ABA regulation towards seed development and storage material accumulation in castor developing seeds.
Effects of exogenous ABA on seed development
Effect of exogenous ABA on seed reserve parameters
Protein content (mg/seed)
Lipid content (mg/seed)
Sugar content (mg/seed)
Dry Weight (mg/seed)
5.9 ± 1.3
1.78 ± 0.4
12.58 ± 2.5
25.9 ± 1.2
5.2 ± 1.1
3.23 ± 0.8*
15.01 ± 2.8*
27.4 ± 1.6
In addition, compartmentalization of lipid classes (neutral lipid-NL) mediated by exogenous ABA (10 μM) based on the TLC plate analysis was explored (Additional file 1: Figure S1). Triricinolein (TR3), the most prominent component in castor bean, is the high polar triglyceride. Triglycerides containing two ricinoleates and one unhydroxylated acyl moiety are less polar (TR2). Triglycerides of two or three unhydroxylated acyl moieties (TR1) are least polar and are found close to the solvent front. TR3 castor lipid class was found to be highly responsive by ABA treatments compared to the control samples followed by TR2 the second major lipid class of castor lipids was also responsive but not as high as that of TR3 comparing to the control sample density (Additional file 1: Figure S1). Minimal increase was observed in TR1 the least triglyceride density compared to control spot percent (Additional file 1: Figure S1). There were no any visual differences in seed morphology observed after ABA treatments during the course of the experiments. To compare the prolific changes caused by ABA signaling towards storage reserves at this stage, lipid content, sugar content and protein content were also measured in developing seeds at different developmental stages (7-63 DAP). From our results, protein synthesis was significantly high during the early-mid seed filling stages (7-35 DAP; 0.3-27.7 mg−1seed) comparing to lipid deposition which increased rapidly from 21-35 DAP (0.95-30 mg−1seed) although minimal amounts detected during early stages of seed filling (7-14 DAP) (Additional file 2: Figure S2). Synthesis of lipid and protein components declined on a slow rate during the maturation stages (42-63 DAP). On the other hand, high amount of sugar levels were noted only during the early stages (7-21 DAP; 0.8-6.7 mg−1seed) after which a rapid decline during the mid-developmental as well as the maturation stages was observed (Additional file 2: Figure S2). This analysis of metabolite levels in developing seeds (in vivo) is of fundamental importance to understand the changes caused by ABA treatment at 21 DAP (in vitro cultured). The advancement of lipid deposition by ABA in developing seeds depends on the percentage of intrusion in to the seed tissues. To make sure the elicitation of exogenous ABA, we measured the ABA content both in culture medium and treated seeds. The culture medium analyses were performed at 3 different time intervals 6, 24, 72 h after treatment. After 6 h of treatment about 79% of free ABA was detected in the medium which was highly stable until 24 h treatments (24%) after which only minimal amounts (1.5%) were detected in 72 h culture medium. Preference was given to the initial stages (24 h) in treated seeds to limit the catabolism of free form ABA in to its metabolites. Relatively, higher content of ABA was noted in seeds treated (24 h) with exogenous ABA than the control samples (Additional file 3: Figure S3). The data presented (means of five replicates) were statistically significant as determined by student’s t-test. The significant changes of storage reserves in castor seed filling after exogenous ABA treatment suggested that ABA play a critical role in mediating storage material metabolism in seed filling.
Global gene expression profile
Categorization and summary of abundance of RNA-Seq tags in two libraries (control and ABA treated)
Total tags (raw tags)
Tags mapped to gene
Unique tags mapped to gene
Tags mapped to genome
ABA mediates differential gene expression in castor developing seeds
ABA responsive genes associated with storage material metabolisms
Contrarily, some crucial genes participating in the sucrose metabolic process were also down-regulated. These include pyruvate dehydrogenase [28883.m000717], 3-phopshoglycerate dehydrogenase [3-PGA] [29801.m003106], NADP malic enzyme [29794.m003406], NADH dehyrogenase [29726.m003961], UDP-glucosyltransferase [29970.m000993], cellulose synthase [29863.m001055], UDP glucose 6-dehydrogenase [30174.m009132] and PEP carboxynase (ATP) [28179.m000470]. As a result of the mediation of crucial sucrose metabolic enzymes; rate limiting enzymes participating in lipid biosynthesis as well as metabolic pathway (Figure 4) were highly triggered. Of the positive responsive transcripts due to the effect of exogenous ABA, several rate limiting enzymes involved in acyl CoA dependent lipid synthesis pathway like acyl coA synthase (ACS) [29844.m003365], Acyl coA binding protein (ACBP) [30170.m013760], type 1 diacylglycerol acyl transferase (DGAT1) [29912.m005373], ER glycerol-3-phosphate acyltransferase (GPAT 2 and 3) [30174.m008615, 29969.m000267], and desaturase genes participating in acyl coA independent lipid synthesis pathway like acyl carrier protein desaturase (SAD1) [29929.m004514], omega-6-desaturase (FAD2) [29696.m000105], fatty acid desaturase (FAD3) [30174.m008679] genes were observed to be highly modulated with a high fold change (≥1) in castor bean developing seeds at this stage (Additional file 4: Table S1). Lipid breakdown genes like long chain fatty acid ligase (LCFA) [30076.m004616], phospholipase D (PL) [29784.m000369], monoglyceride lipase (MAGL) [29986.m001619], diacylglycerol synthase (DAGS) [28726.m000069] and triacylglycerol lipase (TAGL) [30193.m000713, 29158.m000191, 29726.m004101, 29904.m003044, 29620.m000558] were also influenced by ABA signaling. Some significant enzymes related to fatty acid β-oxidation like 3 ketoacyl CoA thiolase B (KAT2) [29904.m002984], 3-OH hydroxy acyl CoA dehydrogenease (29668.m000313) were also triggered by exogenous ABA activity (Additional file 4: Table S1).
In addition, cytochrome b5 [30213.m000673] a rate limiting enzyme (electron supply) for the synthesis of ricinoleic acid in castor bean seeds showed high expression comparing to control levels (Additional file 4: Table S1). On the other hand, critical fatty acid synthesis genes like ketoacyl synthase (KAS) [29693.m002034], keto acyl reductase (KAR) [29929.m004732], biotin carboxylase subunit of Het-ACCase (BACC, BCCP) [30185.m000954, 29630.m000809] and fatty acid desaturase [29794.m003308] were found to be down regulated by the exogenous ABA signaling (Additional file 4: Table S1). Although the total protein content was slightly decreased (Table 1), we examined the gene models encoding protein showing differential response for exogenous ABA as oil (45%) and protein (34%) comprise the major proportion of the storage reserve complex in castor bean seeds. 7S globulin [29993.m001039, 29844.m00331], cytochrome P450 [29970.m001003 (CYP81D1), 30170.m013965 (CYP2), 30170.m013873 (CYP707A4), 29813.m001518 (CYP72B1), 30115.m001196 (CYP707A4), 30174.m008617 (CYP734A1), 30076.m004534 (CYP72B1)], late embryogenesis abundant (LEA5) [29601.m000448, 29726.m004064], early nodulin [29983.m003167], zinc finger protein [28152.m000908, 28637.m000203, 29589.m001255, 30170.m013751, 29589.m001278], RING H2 finger protein [29810.m000139, 29986.m001641, 29048.m000065, 29983.m003157, 30078.m002255] and SGP1 G-protein [29666.m001485], were among the protein encoding genes playing vital role in seed development positively regulated with fold change ≥1 (Additional file 4: Table S1). Several other critical genes encoding storage and developmental proteins were down regulated under ABA signaling (Additional file 4: Table S1). These include cytochrome P450 [30094.m000683 (CYP711A1), 30128.m008568 (CYP51)], ring H2 finger protein [30068.m002520, 29603.m000511], glycine rich cell wall protein [30010.m000662], proline rich cell wall protein [29724.m000836] and cell division control protein [30171.m000404].
The intensity of ABA signaling pathway is well known to be controlled by multiple factors. Several physiological responses of seeds are closely regulated by the network of transcription factors (TFs). In relation to this, several seed specific transcription factors analyzed showed differential response for exogenous ABA. Among such MYB factor [30190.m011160, 30055.m001589], r2r3 MYB factor [30066.m000714; 28154.m000040; 29686.m000890; 30170.m013698], TINY [29938.m000611], F-box and wd 40 domain protein factor [29851.m002476], CCAAT binding factor (NF-YB6) [29656.m000483; 29912.m005288], WRKY factor [30147.m014474, 30174.m008663], binding factors related to ABI3/VP1 (RAV1) [29738.m001050], long chain fatty acid cofactor [27446.m000490], NAC domain factor [30076.m004490, 30076.m004486], bZIP [29634.m002076, 29646.m001108, 29832.m000310, 30076.m004659, 30131.m007163, 30170.m013868] and Bhlh factor [29827.m002528] showed up-regulated expression (Additional file 4: Table S1). These TFs identified belong to several families including bZIP, CCAAT, PHD, ARF, bHLH, NAC, MYB, WRKY and AP2/EREBP and of these families, number of TFs (homologous) identified have been previously reported to regulate seed development  (Additional file 4: Table S1). Very few transcription factors belonging to CCAAT, NAC domain, MYB, PHD family were down regulated by the influence of exogenous ABA (Additional file 4: Table S1). As shown in Additional file 4: Table S1, some of the genes directly participating in pentose phosphate pathway were also strongly influenced by exogenous ABA. Genes like ribose 5 phosphate isomerase [29932.m000590], Ribose-phosphate pyrophosphokinase [30054.m000794], 1-deoxyxylulose 5 phosphate synthase [29726.m003963], trehalose 6 phosphate [28863.m000088, 30147.m013978], which mediate support to sucrose assimilation process in plants oil seeds were also responsive to exogenous ABA signaling. The response of positively modulated metabolic genes were in high ratio as of lipid genes when compared to negatively regulated genes (Figure 1).
Hormone metabolism is triggered by exogenous ABA signaling
Cross talk between ABA and other hormone signaling pathways were analysed, in particular genes participating in gibberellin signaling, ethylene-mediated signaling, auxin-mediated signaling to exogenous ABA signal. Genes like gibberallin receptor GID1 [30170.m013907], indole 3 acetic acid synthase GH3.6 [27533.m000079], chitin inducible giberallin responsive protein [29648.m001919], ethylene receptor [29603.m000534], brassinosteroid regulated protein BRU1 [30065.m001165], auxin induced protein 5NG4 [29783.m000310], auxin responsive protein IAA1 [30169.m006326, 29841.m002749], auxin repressed protein [30190.m011325], auxin carrier protein [29816.m000677], ethylene Insensitive 3 [29708.m000185] showed high induction fold towards exogenous ABA (Additional file 4: Table S1). In addition, 9-cis-epoxycarotenoid dioxygenase (NCED1, 3) [30202.m000262, 29693.m001992] gene along with beta carotene hydrolase [30169.m006624] which function as key enzyme in ABA biosynthesis were highly expressed indicating the up regulation of ABA levels in the castor developing seeds as a result of exogenous ABA activity (Additional file 4: Table S1). NCED1 was the predominant ABA biosynthesis gene with an induction fold of 8.11 towards exogenous ABA signaling. Many protein phosphatase 2C [29794.m003349, 29912.m005442, 30169.m006515, 28211.m000131, 30189.m001658] genes which are negative regulators (i.e., low expression levels in seed) for ABA in seed were highly up-regulated. Contrarily, several genes related to hormone signaling involved in seed development like gibberallin 2 oxidase [29851.m002508], Indole 3 acetic acid ARG7 [30131.m007151], Brassinosteroid receptor kinase [27985.m000860], Brassinosteroid insensitive 1 [29631.m001053], Chitin inducible giberralin protein [29949.m000124] were down regulated by ABA signaling (Additional file 4: Table S1).
sqRT-PCR analysis confirms DGE validation
Seed filling, a phase defined by morphological, cellular, and metabolic changes in the endosperm and embryo that coincide with rapidly increasing storage reserves, determine seed development and storage reserve accumulation . Although diverse functioning ABA responsive factors were identified in plants, relatively limited information is known on ABA responsive genes involved in storage reserve accumulation during seed filling. Thus, the present study is focused on analyzing the effect of exogenous ABA signaling on storage reserve accumulation during seed filling at both physiological as well as at transcriptional level.
Since most metabolic pathways of storage reserve biosynthesis are branched, and there are alternation routes, futile cycles, and product turnover in developing seeds, the molecular mechanisms determining the differential partitioning of seed reserves into the major storage components remain largely unknown . Lipid and protein are the major storage reserves in castor bean seeds . Earlier studies reported that exogenous ABA promoted lipid accumulation in developing seeds of Picea glauca, developing wheat embryos  and embryo cultures of rapeseed  under in vitro conditions. In support to this, our current study demonstrated that exogenous ABA enhanced the sugar accumulation (6.3% increases) followed by a significant storage lipid enhancement (4.9% increases) in castor developing seeds. Triacylglycerols (TAGs) are the major storage lipids in diverse oilseeds . Storage TAGs in the developing seeds of castor bean are highly enriched with unusual fatty acids comparing to phospholipids . Thus our study particularly focused on neutral lipids containing different class distribution of TAGs (Additional file 1: Figure S1). As a highly reduced form of carbon, TAGs serve as an important energy reserve in plant seeds, providing for subsequent germination and seedling development. Marked effect of exogenous ABA towards TAG lipid classes, in particular TR2 and TR3 denoted the inductive action favoring unusual fatty acid synthesis, similar to the effect observed in isolated developing wheat embryos . A major metabolic activity in developing oilseeds is the conversion of sucrose via glycolysis to lipids , indicating that exogenous ABA might have resulted in an increased conversion of sucrose via glycolysis to lipids in the cultured seeds. Previous observations found that exogenous ABA enhanced the accumulation of total sugars in cultured seeds of soybean and Arabidopsis because of the high uptake of sucrose from nutrient media as a result of alteration in the osmotic potential caused by exogenous ABA signals [31–33]. In support to this, sugar levels in particular were measured in the present study as they are the major transport carbohydrate in oilseeds and also work as an important signaling molecule that regulates genes involved in photosynthesis, metabolism, and developmental processes . During seed germination and seedling development, storage TAGs are cleaved by lipases from their glycerol backbone in the oil body, and convert to carbohydrates major as sucrose for transport to the root and shoot axes . From the prolific comparison of metabolite levels after ABA treatment in developing seeds; it is noteworthy to mention that high uptake of sucrose (after ABA treatment) from the culture medium as well as the high proportion of sucrose synthesis (in vivo at 21 DAP) at this stage of development might have added to the pushing of carbon fluxes towards TAG accumulation along with a strong interaction between ABA and other carbohydrate metabolites (Additional file 2: Figure S2). This is in agreement with Schussler et al.  who elaborated the ABA stimulation on high uptake of sucrose in excised soybean embryos cultured under in vitro conditions . Based on our results and other reports it seems reasonable to think that the substantial increase in carbon assimilation in the form of sugars towards the TAG accumulation mediated by exogenous ABA might be the possible reason for the significant increase in lipid deposition without causing much difference in the dry weight proportion in developing castor seeds.
The de novo synthesis of fatty acids in plants occurs in the plastids through the activity of fatty acid synthetase. The synthesis of the malonyl-coA enzyme that is required for acyl-chain elongation requires the import of metabolites (participating in the sucrose breakdown) from the cytosol and their subsequent metabolism. The biosynthesis of lipids in plant seeds is dependent on the supply of many cytosolic metabolites, such as sucrose, acetate, malic enzymes, glucose 6-phosphate, pyruvate, glycerol 3-phosphate . In the current study, enhancement of sugar sink strength may directly push the carbon fluxes toward lipid biosynthesis, or indirectly enhance sugar metabolism, and correspondingly provide much more precursors for seed lipid biosynthesis. In addition, our current study showed that exogenous ABA application slightly decreased the protein content in cultured castor seeds. Exogenous ABA transiently decreased the protein content at the early stages (21 DPA) of grain filling period in wheat . Similar effect of exogenous ABA on enhancing lipid deposition with decrease in protein content has been previously reported in the embryos of Brassica napus L . Decrease in the enzyme activity of ammonia assimilation and substrate concentration of protein synthesis due to exogenous ABA action might have caused the reduction in protein concentration in castor developing seeds in agreement with the previous observation during wheat grain filling . In addition, ABA induced a decline in total nitrogen fixation resulting in a decrease in protein content after 5d of treatment independent of sucrose synthase (SS) activity (enzymes participating in starch synthesis) in seeds of pea . These clearly state that exogenous ABA induces a switch of the conversion from carbon fluxes toward lipids and proteins. Contrarily, some studies showed that the exogenous ABA enhanced accumulation of seed storage proteins in soybean embryos , Arabidopsis and date palm  independent of other storage metabolic pathway. From our results and previous reports, it is evident that exogenous ABA effect on protein accumulation in in vitro cultured seed is stage specific and species dependent.
Through DGE analysis, we obtained a total of 2568 differentially expressed transcripts, in which 1507 transcripts were expressed significantly higher and 1061 transcripts had a down regulated expression in ABA treated library (Figure 1). This discrepancy is a reflection of the effective signaling of ABA in to the seed tissues, in which genes associated with seed metabolic pathways are likely highly regulated. Among these differentially expressed transcripts, more number of cytosolic metabolites participating in sucrose breakdown was highly expressed in ABA treated seeds, respectively (Additional file 4: Table S1). Positive utilization of sucrose metabolites like NADPH, malic enzyme (ME), Glutamate-to-γ-Aminobutyrate to improve storage reserves have been reported in both plant and microbial systems . Analysis on differential expression of genes encoding enzymes involved in fatty acid and TAG biosynthesis may provide clues to identify the key genes that play pivotal roles in the limiting steps of storage lipid biosynthesis which is consistent with the previous observation of ABA influence towards triglyceride content and lipid biosynthesis genes in seed cultures of Brassica napus[7, 43]. Several rate limiting transcripts related to lipid biosynthesis and metabolism were found at higher expression levels in ABA treated seeds (Additional file 5: Table S2). Associated with lipid synthesis and metabolism, transcripts of acyl coA synthase (RcACS), diacylglycerol acyltrasnferase (RcDGAT), ER glycerol-3-phosphate (RcGPAT), fatty acid desaturase (RcFAD), acyl carrier desaturase (RcSAD), triacylglycerol lipase and long-chain-fatty-acid-CoA ligase (LCFA) were shown significantly higher expression levels After ABA treatment. Furthermore, six highly expressed bZIP transcription factors identified may be involved in regulating storage reservoir accumulation during seed filling in castor bean. In addition, differentially expressed transcripts were significantly enriched in seven GO terms, suggesting that ABA signaling may be complex and participate in diverse biological processes during seed filling.
The sqRT-PCR technique provided a sensitive and specific method to analyze multigenic expressions, enabling quantification of the weakly expressed transcripts in the picogram range of RNA allowing diverse studies of gene expression . Fully replicated sqRT-PCR analysis, to certain extent, complements the lack of DGE replications. All eight lipid genes subjected to the sqRT-PCR analysis in this study showed the similar trends with the DGE results, suggesting that our results be confirmed.
In conclusion, the present study showed that the application of exogenous ABA enhanced the storage reserve accumulation together with the global regulatory changes at transcriptional level towards seed filling. These findings provide a comprehensive coverage to discover many known genes of several major metabolic pathways substantially responsive to exogenous ABA signaling, which also serve as a potential contribution to available sequence resources for castor bean. From our results and diverse assessment, it was conclusive that ABA exhibits a stage specific and species dependent functioning when applied as an external supplement during seed filling process under in vitro conditions. These results serve as basic information to dissect the physiological mechanisms underlying ABA signaling towards seed filling and storage reserve accumulation in castor bean seeds.
Materials and methods
Seeds of castor bean var. ZB306 (kindly provided by Zibo Academy of Agricultural Sciences, Shandong, China) were germinated in Xishuangbanna Tropical Botanical Garden (21°560 N, 101°150 E, 600 m asl) under natural conditions. Mature female flowers were hand pollinated and tagged as 0 days after pollination (DAP). Developing seeds with capsulate at ca. 21 DAP were collected from the field. Young seeds were separated from capsulates and were stored at room temperature (18-25°C) for several hours. Decapsulated seeds were inoculated in MS media (Murashige and Skoog) supplemented with exogenous ABA (Sigma, China), as described previously , at 1, 10, 50 and 100 μM concentrations to sort out the optimal ABA concentration for further test. Same treatments without ABA in MS media were used as control. Once the optimal ABA concentration was determined, two decapsulated seeds from the same capsulate were considered as a pair for in vitro seed culture, one was cultured with ABA in MS media, and other as control (without ABA). At least five individual experiments were replicated. Tested seeds were dried overnight in Eppendorf tubes under vacuum centrifugation. The dried seeds were weighted (DW, dry weight). The contents of storage materials including total lipids, proteins and total sugars were measured for each cultured seed after 96 h incubation in dark at 25-28°C as well as in developing seeds. For the preparation of RNA samples, three seeds cultured with ABA for 24 h and the other three seeds from same capsulate cultured without ABA were collected as tissues for extracting RNA.
Measurement of storage materials
For lipid extraction, five replicate samples of cultured seeds (for cultured ABA seeds or controls, respectively) were weighted and frozen in liquid nitrogen for a few minutes, and ground to fine powder in a mortar. The fine powder was transferred into 10 ml glass tubes, and the total lipid was extracted using hexane–isopropanol method described previously . For total protein measurement, as described by Borek et al. , five replicate samples of cultured seeds for ABA or controls were weighed and homogenized (Fluko, China) using 0.1 M TRIS-HCl buffer, pH 7.8. The homogenate was centrifuged for 20 min at 22 000 g. Protein concentration in the supernatant was determined according to , using bovine serum albumin as a standard . For sugar estimation, a previously described dinitrosalicylic acid (DNS) colorimetric method was applied with minor modification . Briefly, seeds of 0.5 g in weight were homogenized and extracted twice with 90% ethanol. The supernatant (ethanol extracts) were collected and pooled to a final volume of 10 ml with double distilled water. A suitable aliquot (1 ml) was taken from the extract and 20 μl of concentrated hydrochloric acid was added and allowed to hydrolysis for 5 min at 90°C after which 0.05 ml of 5N KOH was added to neutralize the acid (alkaline condition for color development). Finally, 0.5 ml of DNS solution (0.1%) was added to the solution and kept for 10 min at 95°C boiling water. The solution was then cooled for color development and the absorbance of this solution was measured at 540 nm using a UV-VIS spectrophotometer (Model DU 640B, Beckman, USA). Student’s t-test was statistically used for all significance tests.
High-throughput RNA-tag sequencing and analyses
The total RNA was isolated from three replicated seeds cultured with ABA for 24 h and three control samples, respectively. RNA extraction was performed according to the manufacturer’s instructions of Trizol reagent (Invitrogen, USA) and purified using an RNeasy Mini Kit (Qiaden, Hilden, Germany). The total RNA was checked for quality and quantity using a Nano-Drop ND-1000 spectrophotometer, and the high quality RNA was used to construct tag libraries respectively for deep-sequencing. Briefly, total RNA (about 6 μg) was enriched by oligo (dT) magnetic beads and oligo (dT) beads were used as primer to synthesize the first and second strand cDNA. The cDNA was digested by two types of endonuclease: Nla III or Dpn II, acquiring 17 bp tags with different adaptors of both ends to form a tag library. After 15 cycles of linear PCR amplification, 105 bp fragments were purified by 6% PAGE gel electrophoresis. After denaturation, the single-chain molecules were fixed on to the Illumina Sequencing Chip. Each molecule was then allowed to grow into a single-molecule cluster sequencing template through in situ amplification. Four colored nucleotides were added for sequencing using the method of sequencing by synthesis (SBS). Millions of raw reads were generated with a sequencing length of 49 bp. Sequencing was performed in BGI Shenzhen (China). The raw data from the ABA treatment and control tag libraries of castor bean were preprocessed to filter out low quality reads and clipped adapter sequences. After that, all clean reads were mapped to the castor bean genome (http://castorbean.jcvi.org/index.php) to obtain unique reads and reads abundance using SOAP2 software . To compare the differential expression of genes between ABA treatment and control libraries, the expression level of each gene in the libraries was normalized to the number of transcripts per million (TPM). Genes with significantly different expression were determined by P < 0.05 and fold-change ≥ (log2 (sample1/sample2)) in two samples. To dissect the potential functions of genes, GO enrichment and KEGG pathway analyses were performed. For enrichment analysis of the function and pathway, we mapped all differentially expressed genes to terms in GO and KEGG database.
Semi-quantitative RT-PCR analysis
Total RNA from the same samples used for high-throughput Tag-Seq sequencing was used for sqRT-PCR experiments. A first-strand cDNA fragment was synthesized from total RNA using Superscript II reverse transcriptase (Invitrogen, USA). Gene-specific primers were designed according to the gene sequences using the Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA). RT-PCR amplifications were performed in a GeneAmp® PCR System (Applied Biosystems). The optimized PCR reactions (25 μl) contained 200 ng of first-strand cDNA template, 2.5 μl 10 × PCR buffer, 0.5 μM of each primer, 0.25 mM of dNTPs, and 1.2 μM of Trans-start™ HiFi DNA polymerase (Transgen, Beijing, China). The following PCR conditions were used: initial denaturation at 95°C for 1 min, 30 cycles of denaturation at 94°C for 30 sec, annealing at 55-58°C for 20 sec and extension at 72°C for 5 min. The castor bean Actin gene was used as internal control to normalize the relative amount of mRNAs for all samples. The amplified products were analyzed by electrophoresis system in 1% agarose gels and visualized by ethidium bromide staining. The primers used in this study are listed in Additional file 5: Table S2.
Extraction and quantitation of endogenous ABA levels in treated seeds and culture medium
Pre-weighed seed homogenates were prepared from harvested seeds as described by Holbrook et al.  and adjusted to give an FW equivalent of 250 mg FW/mL homogenate. Homogenates equivalent to 1g FW were lyophilized for 5 to 8 h 1% butylated hydroxytoluen (BHT). The residue was stirred at room temperature overnight, centrifuged at 2500 g for 20 min to pellet undissolved solids, and the supernatant, containing the free ABA and metabolites, was collected. In addition, Samples of culture media (24 h treated) were filtered through 2-μm nylon filters directly into sample vials . Because of the liability of the free ABA, samples were analyzed immediately or frozen at -80°C overnight and analyzed using a simplified ELISA technique (free of internal matrix) .
The authors thank Mr. Tao faqing, SW China Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences for his timely assistance in conducting ELISA experiments.
- Fitzpatrick AH, Shrestha N, Bhandari J, Crowell DN: Roles of farnesol and ABA in Arabidopsis flower development. Plant Signal Behav. 2011, 6: 1189-1191.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng WH, Endo A, Zhou L, Penny J, Chen HC, Arroya A, Leon P, Nambara E, Asami T, Seo M, Koshiba T, Sheen J: A unique short-chain dehydrogenase/reductase in Arabidopsis glucose signaling and abscisic acid biosynthesis and functions. Plant Cell. 2002, 14: 2723-2743.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Adams IP, Ratledge C: Malic enzyme: the controlling activity for lipid production? Overexpression of malic enzyme in Mucor circinelloides leads to a 2.5-fold increase in lipid accumulation. Microbiol. 2007, 153: 2013-2025.View ArticleGoogle Scholar
- Kanno Y, Jikumaru Y, Hanada A, Nambara E, Abrams SR, Kamiya Y, Seo M: Comprehensive hormone profiling in developing seeds of Arabidopsis seeds: Examination of the site of ABA biosynthesis, ABA transport and hormone interactions. Plant Cell Physiol. 2010, 51: 1988-2001.View ArticlePubMedGoogle Scholar
- Vahdati K, Bayat S, Ebrahimzadeh H, Jariteh M, Mirmasoumi M: Effect of exogenous ABA on somatic embryo maturation and germination in persian walnut (Juglans regia L.). Plant Cell Tissue Organ Cult. 2008, 93: 163-171.View ArticleGoogle Scholar
- Chen Z, Huang J, Muttucumaru N, Powers SJ, Halford NG: Expression analysis of abscisic acid (ABA) and metabolic signaling factors in developing endosperm and embryo of barley. J Cereal Sci. 2013, 58: 255-262.PubMed CentralView ArticlePubMedGoogle Scholar
- Zou J, Abrams GD, Barton DL, Taylor DC, Pomeroy MK, Abrams SR: Induction of lipid and oleosin biosynthesis by (+)-abscisic acid and its metabolites in microspore-derived embryos of Brassica napus 1. cv Reston. Plant Physiol. 1995, 108: 563-571.PubMed CentralPubMedGoogle Scholar
- Carrier DJ, Kendall EJ, Bock CA, Cunningham JE, Dunstan DI: Water content, lipid deposition, and (+)-abscisic acid content in developing white spruce seeds. J Exp Bot. 1999, 337: 1359-1364.View ArticleGoogle Scholar
- Chandrasekaran U, Liu A: Endogenous abscisic acid signaling towards storage reserve filling in developing seed tissues of castor bean (Ricinus communis L.). Plant Growth Regul. 2014, 72: 203-207.View ArticleGoogle Scholar
- Attree SM, Pomeroy MK, Fowke LC: Manipulation of conditions for the culture of somatic embryos of white spruce for improved triacylglycerol biosynthesis and desiccation tolerance. Planta. 1992, 187: 395-404.View ArticlePubMedGoogle Scholar
- Karmoker JL, van Steveninck RFM: The effect of abscisic acid on sugar levels in seedlings of Phaseolus vulgaris L. cv. redland pioneer. Planta. 1979, 146: 25-30.View ArticlePubMedGoogle Scholar
- Hauser C, Kwiatkowski J, Jung J, Grossmann K: Accumulation of abscisic acid in cell suspension cultures of oilseed rape treated with the growth retardant BAS 111.W: Effects on osmotic potential and potassium, water and sugar Contents. J Plant Physiol. 1992, 6: 747-753.View ArticleGoogle Scholar
- Suzuki M, Latshaw S, Sato Y, Settles AM, Koch KE, Hannah LC, Kojima M, Sakakibara H, McCarty DR: The maize viviparous8 locus, encoding a putative ALTERED MERISTEM PROGRAM1-Like peptidase, regulates abscisic acid accumulation and coordinates embryo and endosperm development 1[W]. Plant Physiol. 2008, 146: 1193-1206.PubMed CentralView ArticlePubMedGoogle Scholar
- Finkelstein RR, Gampala SSL, Rock CD: Abscisic acid signaling in seeds and seedlings. Plant Cell. 2002, 14: 15-45.Google Scholar
- Xue LJ, Zhang JJ, Xue HW: Genome wide analysis of the complex transcriptional networks of rice developing seeds. PLoS One. 2012, 7: e31081.PubMed CentralView ArticlePubMedGoogle Scholar
- Qi XH, Xu XW, Lin XJ, Zhang WJ, Chen XH: Identification of differentially expressed genes in cucumber (Cucumis sativus L.) root under water-logging stress by digital gene expression profile. Genomics. 2012, 99: 160-168.View ArticlePubMedGoogle Scholar
- Wu P, Qin Z, Zhao W, Zhou X, Wu T, Xin M, Guo Q: Transcriptome analysis reveals differentially expressed genes associated with propamocarb response in cucumber (Cucumis sativus L.) fruit. Acta Physiol Plant. 2013, 35: 2393-2406.View ArticleGoogle Scholar
- Chen GQ, Turner C, He X, Nguyen T, McKeon TA, Laudencia-Chingcuanco D: Expression profiles of genes involved in fatty acid and triacylglycerol synthesis in castor bean ( Ricinus communis L). Lipids. 2007, 42: 263-274.View ArticlePubMedGoogle Scholar
- Akpan UG, Jimoh A, Mohammed AD: Extraction, characterization and modification of castor seed oil. Leonardo J Sci. 2006, 8: 43-52.Google Scholar
- Schnell R, Staak O, Borchmann P, Schwatz C, Matthey B, Hansen H, Schindler J, Ghetie V, Vietta ES, Diehl V, Engert A: A Phase I study with an anti-CD30 ricin A-chain immunotoxin (Ki-4.dgA) in patients with refractory CD30+ Hodgkin’s and non-Hodgkin’s lymphoma. Clin Cancer Res. 2002, 8: 1779-1786.PubMedGoogle Scholar
- Fidias P, Grossbard M, Lynch TJ: A phase II study of the immunotoxin N901-blocked ricin in small-cell lung cancer. Clin Lung Cancer. 2002, 3: 219-222.View ArticlePubMedGoogle Scholar
- Qiu L, Chun Y, Bo T, Junbo Y, Ai-Zhong L: Exploiting EST databases for the development and characterization of EST-SSR markers in castor (Ricinus communis L.). BMC Plant Biol. 2010, 10: 278.PubMed CentralView ArticlePubMedGoogle Scholar
- Greenwood JS, Bewley JD: Seed development in Ricinus communis (castor bean). I, descriptive morphology. Can J Bot. 1981, 60: 1751-1760.View ArticleGoogle Scholar
- Kharenko OA, Zaharia LI, Giblin M, Eric VC, Taylor DC, Palmer CD, Abrams SR, Loewen MC: Abscisic acid metabolism and lipid accumulation of a cell suspension culture of Lesquerella fendleri. Plant Cell Tissue Organ Cult. 2011, 105: 415-422.View ArticleGoogle Scholar
- Houston NL, Hajduch M, Thelen JJ: Quantitative proteomics of seed filling in castor: comparison with soybean and rapeseed reveals differences between photosynthetic and non-photosynthetic seed metabolism. Plant Physiol. 2009, 151: 857-868.PubMed CentralView ArticlePubMedGoogle Scholar
- Borek S, Pukacka S, Michalski K, Ratajczak L: Lipid and protein accumulation in developing seeds of three lupine species: Lupinus luteus L., Lupinus albus L., and Lupinus mutabilis sweet. J Exp Bot. 2009, 60: 3453-3466.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez-Sotres R, Black M: Osmotic potential and abscisic acid regulate triacylglycerol synthesis in developing wheat embryos. Planta. 1994, 192: 9-15.Google Scholar
- Cagliari A, Margis-Pinheiro M, Loss G, Mastroberti AA, Mariath JEA, Margis R: Identification and expression analysis of castor bean (Ricinus communis L.) genes encoding enzymes from the triacylglycerol biosynthesis pathway. Plant Sci. 2010, 179: 499-509.View ArticlePubMedGoogle Scholar
- Brown AP, Kroon JTM, Swarbreck D, Febrer M, Larson TR, Graham IA, Caccamo M, Slabas AR: Tissue-specific whole transcriptome sequencing in castor, directed at understanding triacylglycerol lipid biosynthetic pathways. PLoS One. 2011, 7: e30100.View ArticleGoogle Scholar
- Baud S, Dubreucq B, Miquel M, Rochat CH, Lepiniec L: The Arabidopsis Book. Storage Reserve Accumulation in Arabidopsis: Metabolic and Developmental Control of Seed Filling. 10.1199/tab.0113 http://www.aspb.org/publications/arabidopsis/2008, Rockville MD: American Society of Plant Biologists,
- Schussler JR, Brenner ML, Brun WA: Abscisic acid and its relationship to seed filling in soybeans. Plant Physiol. 1984, 76: 301-306.PubMed CentralView ArticlePubMedGoogle Scholar
- Xing X, Zheng G, Ping D, Medrano H: Effect of exogenous abscisic acid and water stress on the growth response of subterranean clover of different genotypes. Acta Bot Sin. 2002, 44: 1425-1431.Google Scholar
- Verslues PE, Bray EA: Role of abscisic acid (ABA) and Arabidopsis thaliana ABA-insensitive loci in low water potential-induced ABA and proline accumulation. J Exp Bot. 2006, 57: 201-212.View ArticlePubMedGoogle Scholar
- Wind J, Smeekens S, Hanson J: Sucrose: Metabolite and signaling molecule. Phytoche. 2010, 71: 1610-1614.View ArticleGoogle Scholar
- Eastmond PJ: SUGAR-DEPENDENT1 encodes a patatin domain triacylglycerol lipase that initiate storage oil breakdown in germinating Arabidopsis seeds. Plant Cell. 2006, 18: 665-675.PubMed CentralView ArticlePubMedGoogle Scholar
- Rawsthorne S: Carbon flux and fatty acid synthesis in plants. Progr Lipid Res. 2002, 41: 182-196.View ArticleGoogle Scholar
- Yang W, Cai T, Ni Y, Li Y, Guo J, Peng D, Yang D, Yin Y, Wang Z: Effects of exogenous abscisic acid and gibberellic acid on filling process and nitrogen metabolism characteristics in wheat grains. Aus J Crop Sci. 2013, 7: 58-65.Google Scholar
- Gonzalez EM, Gailvez L, Arrese-Igor C: Abscisic acid induces a decline in nitrogen fixation that involves leghaemoglobin, but is independent of sucrose synthase activity. J Exp Bot. 2001, 52: 285-293.View ArticlePubMedGoogle Scholar
- Eisenberg AJ, Mascarenhas JP: Abscisic acid and the regulation of synthesis of specific seed proteins and their messenger RNAs during culture of soybean embryos. Planta. 1985, 166: 505-514.View ArticlePubMedGoogle Scholar
- Kagaya Y, Okuda R, Ban A, Toyoshima R, Tsutsumida K, Usui H, Yamamoto A, Hattori T: Indirect ABA-dependent regulation of seed storage protein genes by FUSCA3 transcription factor in Arabidopsis. Plant Cell Physiol. 2005, 46: 300-311.View ArticlePubMedGoogle Scholar
- Sghaier-Hammami B, Jorrín-Novo JV, Gargouri-Bouzid R, Drira N: Abscisic acid and sucrose increase the protein content in date palm somatic embryos, causing changes in 2-DE profile. Phytochem. 2010, 71: 1223-1236.View ArticleGoogle Scholar
- Fait A, Nesi AN, Angelovici R, Lehmann M, Pham PA, Song L, Haslam RP, Napier JA, Galli G, Fernie AR: Targeted enhancement of glutamate-to-γ-aminobutyrate conversion in Arabidopsis seeds affects carbon-nitrogen balance and storage reserves in a development-dependent manner. Plant Physiol. 2011, 157: 1026-1042.PubMed CentralView ArticlePubMedGoogle Scholar
- Jadhav AS, Taylor DC, Giblin M, Ferrie AM, Ambrose SJ, Rose AR, Nelson KM, Zaharia I, Sharma N, Anderson M, Fobert PR, Abrams SR: Hormonal regulation of oil accumulation in Brassica seeds: metabolism and biological activity of ABA, 7’-, 8’- and 9’-hydroxy ABA in microspore derived embryos of B. napus. Phytochem. 2008, 69: 2678-2688.View ArticleGoogle Scholar
- Pray L: The biotechnology revolution: PCR and the use of reverse transcriptase to clone expressed genes. Nat Educ. 2008, 1: 1.Google Scholar
- Xu R, Wang R, Liu A: Expression profiles of genes involved in fatty acid and triacylglycerol synthesis in developing seeds of Jatropha (Jatropha curcas L.). Biomass Bioenerg. 2011, 35: 1683-1692.View ArticleGoogle Scholar
- Bradford MM: A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein dye binding. Anal Biochem. 1976, 72: 248-254.View ArticlePubMedGoogle Scholar
- Negrulescu A, Patrulea V, Mincea MM, Ionascu C, Vlad-Oros BA, Ostafe V: Adapting the reducing sugars method with dinitrosalicylic acid to microtiter plates and microwave heating. J Braz Chem Soc. 2012, 23: 2176-2182.View ArticleGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP short oligonucleotide alignment program. Bioinformatics. 2008, 24: 713-714.View ArticlePubMedGoogle Scholar
- Holbrook LA, Magus JR, Taylor DC: Abscisic acid induction of elongase activity, biosynthesis and accumulation of very long chain monounsaturated fatty acids and oil body proteins in microspore-derived embryos of Brassica nupus L. cv Reston. Plant Sci. 1992, 84: 99-115.View ArticleGoogle Scholar
- Zhao J, Li G, Yi G, Wang B, Deng A, Nan T, Li Z, Li Q: Comparison between conventional indirect competitive enzyme-linked immunosorbent assay (icELISA) and simplified icELISA for small molecules. Anal Chim Acta. 2006, 571: 79-85.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.