Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genotype-dependent changes of gene expression during somatic embryogenesis in oil palm hybrids (Elaeis oleifera x E. guineensis)

Abstract

To understand the molecular processes triggered during the different steps of somatic embryogenesis (SE) in oil palm, the expression of 19 genes associated to SE identified in proteomic and transcriptomic studies was investigated by qRT-PCR. To evaluate the differential expression of these genes, two interspecific hybrid genotypes (Elaeis oleifera x Elaeis guineensis) contrasting for the acquisition of embryogenic competence were used. Aclorophyllated leaves of both hybrids, one responsive (B351733) and the other non-responsive (B352933) to SE were submitted to callus induction and collected at different time points: 0 (before induction), 14, 30, 90 and 150 days of callus induction (doi). The results obtained showed that all evaluated genes were downregulated at 14 doi in the responsive genotype when compared to the non-responsive. It was also possible to observe that most of the genes changed their expression behavior at 30 doi and were upregulated thereafter until 150 doi, with the exception of the pathogenesis-related PRB1-3-like (PRB1-3) gene, which did not show differential expression at 30 doi and was downregulated at 90 and 150 doi when compared to the non-responsive hybrid. These results indicate that 30 doi is a turning point in gene expression, probably associated to embryogenic competence acquisition. We also show that the expression behavior of the responsive genotype is more stable than that of the non-responsive when the different induction time points are compared to 0 doi (before induction). Moreover, the results obtained in this study corroborate our hypothesis that the regulation of genes involved in the control of oxidative stress and energy metabolism are crucial for the acquisition of embryogenic competence in oil palm.

Introduction

Oil palm is an important crop in the tropics [1] that produces two oils of great economic importance, commonly known as palm oil and palm kernel oil, obtained from the mesocarp and endosperm, respectively [2]. The palm oil is mainly used for cooking, preparing margarine and also for non-food applications (fuel, soap, detergent, cosmetics, etc) [3]. The increasing demand for palm oil is due to the rapid population growth and economic development of several countries [4].

The American palm (Elaeis oleifera), known as caiaue or American oil palm, is found in humid tropical America, widely dispersed from Central to Southern America [57]. E. oleifera is considered a promising genetic resource for oil palm breeding programs, since it represents an important source of genetic variability. Among the beneficial agronomic traits are reduced trunk growth rate, reduced size and lower operating costs, resistance to fatal yellowing, as well as high content of unsaturated fatty acids, which gives greater fluidity to the oil under natural conditions, allowing its exploitation for biodiesel production [810]. The African oil palm (E. guineensis Jacq.) is originated from West Africa [11] and is also known as dendem, coconut-palm and palm of guinea. In the last years, due to the presence of diseases, plantations of interspecific hybrids, obtained from the cross between E. oleifera and E. guineensis (American and African oil palm, respectively), have increased substantially in Latin America due to their apparent partial resistance to the bud rot disease caused by Phytophthora palmivora [12]. E. guineensis Jacq and E. oleifera are two monocotyledonous species with a single apical meristem that does not present tillering, and therefore these species cannot be propagated vegetatively by conventional techniques. Therefore, in these species, micropropagation through somatic embryogenesis (SE) is one of the only alternatives for clonal propagation [13].

SE is the process by which cells differentiate to form embryos, reorganizing their epigenetic properties and cell cycle for tissue formation in morphological stages similar to those obtained in zygotic embryogenesis (ZE) [14]. Since the first reports of SE in E. guineensis by Staritsky [15] and Rabechault et al [16], several studies were performed focused on the induction, proliferation, regeneration and protein profiling in the different SE steps in oil palm [13, 1719]. The optimization of the protocols has been the central aim of the research, however studies on specific stages of the process such as multiplication, maturation and conversion of somatic embryos remain limited in many palm species [14].

Although millions of plants of different species are produced annually through SE, the molecular mechanisms governing the different steps of the process are not well understood [20]. Thus, the key issue in SE is to understand what triggers cells to change their fate and become embryogenic. Proteomic and gene expression studies, including transcriptomic analyses, during somatic embryogenesis in oil palm were performed and several genes potentially involved in SE were identified [2125]. Genes have also been described as specific markers of SE in plants, such as the well known SERK gene (Somatic Embryogenesis Receptor Kinase), first described in competent carrot cell cultures (DcSERK) [26]. SERK has been used as an SE marker to differentiate competent from non-competent cells [27]. Other examples of genes activated and differentially expressed during SE are those encoding glutathione S-transferases (GSTs), which is characteristic of the transition to the development of somatic embryos in in vitro culture systems [28], as well as a number of embryogenesis related genes, such as Late-Embryogenic Abundant (LEA), Baby Boom (BBM), Agamous-like 15 (AGL15), Leafy Cotyledon 2 (LEC2), Leafy Cotyledon 1 (LEC1), and Fusca3 (FUS3), expressed during both zygotic and somatic embryogenesis [29].

Although these genes have been identified as related to SE, a more detailed analysis of their expression patterns is necessary and could certainly help advance our knowledge on the crucial biological processes that trigger embryogenic competence acquisition. In this context, the objective of this work was to evaluate the expression of genes involved in the acquisition of embryogenic competence by qRT-PCR, in oil palm genotypes contrasting for the acquisition of embryogenic competence. Genes previously identified by our group were selected, as well genes described in the literature with possible involvement in SE. The results obtained in this study contribute to a better understanding of the genes involved in embryogenic competence acquisition and can help optimize the protocols by altering environmental factors that in turn modulate gene expression during the different steps of SE.

Materials and methods

Plant material

In this study, two interspecific oil palm F1 hybrids obtained from the cross between E. oleifera and E. guineensis Jacq. were used. These hybrids are unique individuals from the Active germplasm Bank and were selected for cloning propagation since they represent the most productive plants among 573 other plants analyzed for fruit production during 4 consecutive years [30]. However, studies from our group showed that one of them was responsive (B351733) to the SE process while the other (B352933) was not (Scherwinski-Pereira, personal communication). In order to study the differential expression of genes potentially involved in the SE process, a relative expression analysis comparing the responsive (R) and the non-responsive (NR) hybrids at 0, 14, 30, 90 and 150 days of callus induction (doi) was performed. Subsequently, a relative expression analysis of the genes throughout the steps (Time-course) for the R and NR genotypes, using 0 days of induction (doi) as the normalize, was also performed.

Induction of embryogenic callus and sampling

Aclorophyllated leaves of both oil palm hybrids (B351733 and B352933) were collected and submitted to disinfestation with 70% ethanol for 5 min and with 2.5% sodium hypochlorite for 30 min. The leaves were washed with sterile distilled water and then subjected to the callus induction process. For each replicate, 3 Petri dishes were used (10 x 90 mm) containing six explants, measuring about 1 cm2. Each plate contained 25 mL of MS culture medium [31], supplemented with 30 g Lˉ1 sucrose; 0,5 g Lˉ1 glutamine; 0,5 g Lˉ1 of hydrolyzed casein; 2,5 g Lˉ1 of activated charcoal; 2,5 g Lˉ1 Phytagel and 450 μM picloram, as described by Balzon et al [32], with modifications. During the callus induction period, the explants were stored in a growth room, at a temperature of 25 ± 2°C, in the dark. The collection of the plant material (leaf explants) for the analyses was carried out at 0 (before induction), 14, 30, 90 and 150 days of callus induction (doi), macerated in liquid nitrogen and stored at -80°C until further use for RNA extraction. Three biological replicates were analyzed for each hybrid at each time point, totalizing 30 samples.

RNA Purification and cDNA synthesis

Total RNA was extracted from 0.1 g of tissue using the Trizol method according to Simms et al [33], with modifications. The quantification was performed using a NanoDrop spectrophotometer (ND-1000 UV-Vis -Termo Fisher). The integrity of the isolated RNA was confirmed on agarose gel (1% agarose; TAE Buffer 1X), pre-stained with ethidium bromide (0.5 mg / mL) [34]. Before cDNA synthesis, RNA was treated with Turbo DNAse (Applied Biosystems/Ambion) to eliminate possible contaminations with genomic DNA. The cDNA was synthesized using 2 μg of RNA treated using the Next Generation MMLV RNAse H Minus First-Strand cDNA Synthesis (DNA express). The cDNA was stored at -20°C and used for qRT-PCR.

Gene selection, primer design and gene ontology (GO) analysis

Proteins potentially involved in the acquisition of embryogenic competence identified in proteomic studies performed by our research group [35] were selected and the corresponding gene was searched for primer design, including PGM, ACT1, ANN1, CAT2, OASA, ENO1, MDAR5, PFP-BETA, RUBA, PDIL1-4, PRB1, PRB1-3, EF1 and HSP81-1 (Table 1). Some genes were also selected from previous published studies, such as LEA, FIE2 and GBSS1 [36], BGLUC and SERk1[23] as well as the reference genes PD00380 [37] and ACT2 [38] (Table 1). The software Primer3 [39] was used to design all the primers, and the absence of amplification and nonspecific products were evaluated using the software OligoAnalyzer 3.1 [40].

thumbnail
Table 1. General information on the genes selected for gene expression analysis.

https://doi.org/10.1371/journal.pone.0209445.t001

Gene ontology (GO) analysis was performed to understand the biological processes in which these genes are involved using Blast2GO software [41]. The genes were grouped into four categories according their GO terms and according to the literature. Although the same gene can be classified in more than one biological process, in this study only one major biological process was assigned to each gene.

qRT-PCR experiments and data analysis

The qRT-PCR experiments were performed in the thermal cycler 7300 RealTime PCR System (Applied Biosystems). All reactions were composed of 5 μL of Fast SYBR Green Master Mix (Applied Biosystems), 0.2 μL of each primer at an initial concentration of 10 μM (forward and reverse) and 2 μL of single strand cDNA for each sample (diluted 20X). The PCR program used was a step at 95° C for 10 min to activate the Taq polymerase enzyme (hot start), followed by 40 cycles at 95° C for 15 sec, 60° C for 60 seconds. The denaturation curve (melting curve) was analyzed after the end of the amplification to verify the presence of primer dimers and nonspecific products. The program used was 95°C for15 s, 60°C for 60 s, and increase of 0.3°C a up to 95°C. All experiments were performed using three independent biological replicates and for each biological replicate, three technical replicates were analyzed.

The qRT-PCR reactions were performed in 96 well plates and in each plate a negative control without cDNA was placed to confirm the absence of contamination. The raw fluorescence data from all runs were imported into the Real-time PCR Miner software [42], to determine the Cq value and the PCR efficiency. The analyses of expression and statistics were performed using the REST software [43]. The normalization was performed using internal reference genes. To compare the differences in expression between groups, the t-test was used with 4 degrees of freedom, n = 9 (3 biological replicates and 3 techniques for each biological replicate) as determined by the REST software.

Results and discussion

Functional classification

To understand the biological processes in which the 19 genes selected for analysis are involved, a gene ontology (GO) search was performed and revealed important roles related especially to stress responses, defense responses, energy metabolism and development (Table 2). These genes appear to participate in the control of stress and adaptation of the genotypes to the in vitro condition during SE. Similar results were obtained by Tan et al [18] in proteomic analyses from mature leaves of oil palm with high or low proliferation rate. Genes related to defense responses and development were also reported by Lin et al [22].

thumbnail
Table 2. Gene ontology classification (biological process) of the genes potentially involved in embryogenic competence acquisition in oil palm.

https://doi.org/10.1371/journal.pone.0209445.t002

Gene expression analysis by qRT-PCR

In order to evaluate the expression profile of genes potentially involved in the acquisition of embryogenic competence, the relative expression of 19 genes was analyzed comparing the R and NR genotypes at different SE induction steps. We considered 0, 14 and 30 days of induction as initial steps, 90 days of induction as intermediate step and 150 days of induction as later step before the callus formation.

The results showed that most of the genes analyzed in this study were downregulated at 14 doi when the R genotype was compared with the NR genotype (Fig 1). It was interesting to observe that the relative expression of the majority of the genes analyzed in this comparison changed their expression profile at 30 doi and were upregulated until 150 doi. This upregulation at 30 doi indicates that at this time point the mechanisms of response to the acquisition of embryogenic competence were possibly activated in the R genotype. This increase can be explained by the intense cell division and stretching that occurs at the initiation of callus induction, with the need for synthesis of specific proteins that are responsible for morphological and biochemical alterations [44].

thumbnail
Fig 1. Relative expression analysis of the responsive (R) genotype compared to the non-responsive (NR) genotype at 0, 14, 30, 90 and 150 days of callus induction (doi).

A- Response to stress; B- Development; C- Energetic metabolism and D- Defense response. The red and blue bars indicate statistically significant (p≤ 0.05) up or downregulated genes, respectively. Gray bars indicate genes that did not show statistical significance.

https://doi.org/10.1371/journal.pone.0209445.g001

A time-course analysis of the relative expression was also performed for each genotype using the sampling point 0 doi (before induction) as the normalizer. The results showed that the relative expression of the genes throughout the steps in the R genotype is more stable than in the NR genotype (Fig 2). These results suggest that stress levels in the NR genotype are possibly higher in relation to the R genotype, indicating that stress modulation may be crucial for genotype adaptation in the SE induction medium. Indeed, stress adaptation has been reported to influence somatic embryo induction in several plants and promote dedifferentiation [45, 46].

thumbnail
Fig 2. Relative expression analysis of genes throughout the steps (Time-course) for the responsive (R) and non-responsive (NR) genotypes, using 0 days of induction (doi) as the normalizer.

A- Response to stress; B- Development; C- Energetic metabolism and D- Defense response.

https://doi.org/10.1371/journal.pone.0209445.g002

The results obtained indicate that the contrasting responses of acquisition of embryogenic competence between the R genotype and the NR genotype may be related to the increased expression of genes involved in response to stress, defense response, energy metabolism and development, since, at 30, 90 and 150 doi, the expression of the genes involved in these processes had a significant increase in the R genotype when compared to NR genotype. These results indicate that these genes may play an important role in the higher adaptation ability of the R genotype.

Role of stress response genes during somatic embryogenesis

The first days of callus induction during SE are characterized by the modulation of many genes related to stress, which lead to the assumption that acquisition of SE is an adaptation response of cells cultivated in vitro [46]. The control of oxidative stress at adequate levels has been reported as one of the main events for the development SE [47]. The ability to control oxidative stress and preserve protein structure seems to be the key to the success of embryogenic competence acquisition [19]. As part of the response to cellular oxidative stress is the regulation of the expression of genes encoding antioxidant enzymes, which alleviate the damage caused by reactive oxygen species (ROS) [48]. According to Libik et al [49] and Konieczny et al. [50], the increase in the activity of antioxidant enzymes is related to oxidative stress and this stress contributes to accelerate the process of embryogenic response in somatic tissues.

In this study, 6 of the 19 genes evaluated were associated with the response to oxidative stress: CAT2, MDAR5, ACT1, PDIL1-4 and LEA (Table 2). These genes were downregulated at 14 doi and upregulated at 30, 90 and 150 doi in the R genotype when compared to NR (Fig 1A). These results indicate that stress control at 30 doi in the R genotype may be a determining factor for the acquisition of embryogenic competence in oil palm. Therefore, 30 doi seems to be the key moment in which the R genotype adapts to the stress conditions during SE. On the other hand, the decreased expression of the CAT2, MDAR5, ACT1, LEA and PDIL1-4 genes in the NR genotype at 30, 90 and 150 doi may have contributed to the high levels of oxidation that may be one of the causes for the absence of callus formation in this genotype. According to Araldi et al. [51] and Zhou et al. [52], catalase is an antioxidant protein that acts on the regulation of ROS, and its reduced activity results in the accumulation of hydrogen peroxide (H2O2). The activity of catalase in scavenging H2O2 was associated with the embryogenic potential of the calli of Mesembryanthemum crystallinum [49] and Larix leptolepis [53] during SE. Similarly, it was reported in Lycium barbarum culture that in the first days of cell differentiation, catalase activity was decreased and its activity was increased with the division of embryogenic cells [54]. Actin has been associated with the process of programmed cell death of the plant cell [55] and, because of this, is considered an effective signaling marker during this process [56]. MDAR5 also has an antioxidant activity [53, 57]; in Crocus sativus, MDAR5 had a high activity from the nodular embryo to the mature embryo [58], proving to be of paramount importance in the SE process.

The LEA gene is also related to stress response [59], and was upregulated at 30, 90 and 150 doi in the R genotype when compared to NR. This high expression of the LEA gene at these stages is possibly due to the osmotic stress that occurs during the SE induction process. The products of this family of genes may be related to an osmoprotection function, such as defense of cellular structures in mature seed embryos during desiccation and prevention of early germination of zygotic embryos during seed development [60, 61].

In addition to genes encoding antioxidant enzymes, other genes also related to stress in SE, such as defense response genes, were analyzed, including PRB1, PRB1-3 and HSP801-1 (Table 2). These genes were downregulated at 14 doi when the R genotype was compared to NR (Fig 1D). PR proteins are a heterogeneous group of proteins encoded by genes expressed in response to pathogen attack, but also seem to present related roles in the response to abiotic stresses [62]. Therefore, it is possible that the activation of PRs is related to the adaptation of stress in the induction phase of SE in oil palm.

Energy metabolism

Energy metabolism proteins have been reported in other SE studies in E. guineensis showing that some cellular activities under in vitro conditions require high energy levels, which explains the relative abundance of these proteins during some stages of SE development [19, 63]. Many studies have been carried out showing the comparison between somatic embryos and zygotic embryos and it has been described that the energy metabolism is more active in somatic embryos [64, 65]. According to Rolland et al. [66], energy metabolism involves the biosynthesis of carbohydrates and lipids, which store energy and are then broken down and consumed in the form of ATP, providing energy for most cellular activities.

In the present study, some genes were related to energy metabolism such as ENO1, PGM, PFP-BETA, BGLUC, RUBA and GBSS1 (Fig 1C and Table 2). These genes were downregulated at 14 doi and upregulated at 30, 90 and 150 doi in the R genotype when compared to the NR genotype. The behavior of these genes may indicate a high energy expenditure in callus induction. Thus, high energy consumption requires a high production of ATP, and this may explain the increased expression of these genes at 30, 90 and 150 doi in the R genotype indicating a strong connection to the process of acquisition of embryogenic competence.

Lippert et al [67], also described enolase as involved in the SE process, presenting increased expression. According to Barkla et al [68], enolase has as main role in the conversion of 2-phosphoglycerate to phosphoenolpyruvate (PEP) in the glycolytic pathway, converting glucose into pyruvate, resulting in the final product NADH and ATP. Therefore, a high abundance of enolase is expected [69].

It is important to highlight the participation of the RUBA gene in SE. This gene was upregulated at 30, 90 and 150 doi in the R genotype when compared to the NR genotype. This gene has already been reported by Tan et al [18], also showing high abundance in young leaves. It is possible that this gene is related to the photorespiration of the tissue, which prevents the accumulation of toxic products such as 2-phosphoglycerate [70]. Thus, these results may also explain the low oxidative level of the tissues in the R genotype when compared to NR genotype, since RUBA could be acting in the control of oxidative stress.

The PFP-BETA and PGM genes have already been associated with the energy supply during callus induction in Nelumbo nucifera Gaertn. spp. baijianlian, in the formation and development of Coffea arabica torpedo embryos and also in E. guineensis cell suspensions in response to different concentrations of auxin [7173]. The GBSS1 gene is involved in amylose synthesis [74] and histological studies performed by Kanchanapoom and Domyoas [75] and Almeida et al. [35] showed the accumulation of starch granules during the induction of SE oil palm. These results may explain the increased transcription of GBSS1 at 30, 90 and 150 doi in the R genotype when compared to the NR genotype.

Genes associated with development

During SE, biochemical and morphological changes occur throughout the development of the induced tissues, which is closely related to changes in gene expression. Several genes are differentially expressed during induction of somatic embryogenesis, while others are expressed during differentiation and maturation of the embryo for plant development [76]. Among the genes associated with development analyzed in this study are ANN1, EF1, SERK, OASA and FIE2 (Table 2). The ANN1, EF1, SERK and FIE2 genes were downregulated at 14 doi and upregulated at 30, 90 and 150 doi when the R genotype was compared to NR (Fig 1B). It has been described that annexins are related to processes of plant development, such as pollen germination and cotton fiber elongation [77, 78]. Moreover, this gene family has been reported to act in the early steps of somatic embryogenesis from zygotic oil palm embryos, during the somatic embryogenesis of Coffea arabica L. and as a molecular marker of the stages of development of somatic embryos of Manihot esculenta Crantz [19, 73, 79].

Among the genes involved in the induction of SE, the SERK gene plays an important role. This gene was first isolated from embryogenic carrot cells and is considered as a molecular marker for SE [80, 81]. SERK is a leucine rich repeat (LRR) transmembrane protein kinase that increased the capacity of the apical meristem in Arabidopsis to form somatic embryos [82]. The expression of SERK is observed in competent cells until the globular stage of the somatic embryo, however, it is not detected in non-embryogenic stages. The expression of this gene was associated with SE in several species, including Dactylis glomerata [81], Arabidopsis thaliana [82], Medicago truncatula [63] and Helianthus annuus [83]. Interestingly, in our study, the SERK gene showed expression at all time points analyzed (0–150 doi). However, a downregulation was observed at 14 doi followed by an increase from 30–150 doi in the R genotype. These results reinforce our hypothesis that 30 doi is the crucial time point when SE competence is acquired.

Another gene evaluated in our study and related to protein processing was OASA (Table 2). This gene was downregulated at 0 doi and at 14 doi and upregulated at 30, 90 and 150 doi in the R genotype when compared to NR. It has already been described in some studies that in the presence of auxin, OASA also acts in the regulation of oxidative stress, being able to stimulate cell differentiation and promote the formation of somatic embryos [84, 85]. Thus, this may explain the increase of this gene in the R genotype indicating that this gene could be a potential marker of the acquisition of embryogenic competence.

Conclusion

The genes evaluated in this study have important roles related mainly to stress and defense response, energy metabolism and development, participating in the control of oxidation and adaptation of the genotypes to in vitro conditions during SE induction. The results showed that overall the evaluated genes were downregulated at 14 doi and upregulated from 30 doi to 150 doi when the R genotype was compared with the NR genotype. The results presented here suggest that 30 doi is the key moment when the R genotype changes gene expression leading to the development of somatic embryogenesis. In addition, all genes appear to play a role in the acquisition of embryogenic competence and can be explored as potential markers of SE in oil palm, even PRB1-3 gene that showed a different expression profile, with no modulation at 30 doi and downregulation at 90 and 150 doi in the responsive genotype. The results obtained in this study have an important impact in advancing the current knowledge regarding gene expression during somatic embryogenesis and particularly in understanding the differential response of genotypes to SE. This knowledge is of paramount importance for clonal propagation of elite oil palm genotypes.

References

  1. 1. Zheng Y, Li Y, Zhang Y, Ruan X, Zhang R. Purification, characterization, synthesis, in vitro ACE inhibition and in vivo antihypertensive activity of bioactive peptides derived from oil palm kernel glutelin-2 hydrolysates. J Funct Foods. 2017;28:48–58.
  2. 2. Singh R, Ong-Abdullah M, Low E-TL, Manaf MAA, Rosli R, Nookiah R, et al. Oil palm genome sequence reveals divergence of interfertile species in Old and New worlds. Nature. 2013;500(7462):335. pmid:23883927
  3. 3. Sanputawong S, Te-Chato S. Analysis of somaclonal variation of callus, somatic embryo and plant regeneration of in vitro oil palm (Elaeis guineensis Jacq.). J Agricult Technol. 2011;7(2):531–45.
  4. 4. Or KH, Putra A, Selamat MZ. Oil palm empty fruit bunch fibres as sustainable acoustic absorber. Appl Acoust. 2017;119:9–16.
  5. 5. Rajanaidu N. Elaeis oleifera collection in Central and South America. Proceedings of the International Workshop on Oil Palm Germplasm and Utilization, Palm Oil Research Institute of Malaysia (PORIM). 1986:84–94.
  6. 6. Balslev H, Henderson A. Elaeis oleifera (Palmae) encontrada en el Ecuador. Publ Mus Cienc Nat Ecuador. 1986;5:45–9.
  7. 7. Escobar R, editor Preliminary results of the collection and evaluation of the American oil palm Elaeis oleifera (HBK,[Humbaldt Bonbaldt, Kunith] Cortes) in Costa Rica [Panama and Colombia]. International Conference on Oil Palm in Agriculture in the Eighties Session-Breeding, genetics and selection II, Kuala Lumpur (Malaysia), 17–20 Jun 1981; 1982: PPP (ISP).
  8. 8. Hardon J, Tan G. Interspecific hybrids in the genus Elaeis I. crossability, cytogenetics and fertility of F1 hybrids of E. guineensis x E. oleifera. Euphytica. 1969;18(3):372–9.
  9. 9. Hormaza P, Fuquen EM, Romero HM. Phenology of the oil palm interspecific hybrid Elaeis oleifera× Elaeis guineensis. SCI AGR. 2012;69(4):275–80.
  10. 10. Moretzsohn MdC, Ferreira M, Amaral Z, Coelho PJdA, Grattapaglia D, Ferreira ME. Genetic diversity of Brazilian oil palm (Elaeis oleifera HBK) germplasm collected in the Amazon Forest. Euphytica. 2002;124(1):35.
  11. 11. Zeven A. On the origin of the oil palm (Elaeis guineensis Jacq.). Grana. 1964;5(1):121–3.
  12. 12. Torres G, Sarria G, Varon F, Coffey M, Elliott M, Martinez G. First report of bud rot caused by Phytophthora palmivora on African oil palm in Colombia. Plant Dis. 2010;94(9):1163.
  13. 13. Scherwinski-Pereira JE, da Guedes RS, Fermino PCP, Silva TL, Costa FHS. Somatic embryogenesis and plant regeneration in oil palm using the thin cell layer technique. In Vitro Cell Dev Biol Plant. 2010;46(4):378–85.
  14. 14. Ree JF, Guerra MP. Palm (Arecaceae) somatic embryogenesis. In Vitro Cellular & Developmental Biology-Plant. 2015;51(6):589–602.
  15. 15. Staritsky G. Tissue culture of the oil palm (Elaeis guineensisJacq.) as a tool for its vegetative propagation. Euphytica. 1970;19(3):288–92.
  16. 16. Rabechault H, Ahée J, Guenin G. Colonies cellulaires et formes embryoides obtenues in vitro a partir de cultures d'embryons de Palmier a huile (Elaeis guineensis Jacq. var. dura Becc.). Acad Sci Compt Rend Ser D. 1970;270(25):3067–70.
  17. 17. Konan KE, Durand-Gasselin T, Kouadio YJ, Flori A, Rival A, Duval Y, et al. In vitro conservation of oil palm somatic embryos for 20 years on a hormone-free culture medium: characteristics of the embryogenic cultures, derived plantlets and adult palms. Plant Cell Rep. 2010;29(1):1. pmid:19898823
  18. 18. Tan HS, Liddell S, Abdullah MO, Wong WC, Chin CF. Differential proteomic analysis of embryogenic lines in oil palm (Elaeis guineensis Jacq). J Proteomics. 2016;143:334–45. pmid:27130535
  19. 19. de Carvalho Silva R, Carmo LST, Luis ZG, Silva LP, Scherwinski-Pereira JE, Mehta A. Proteomic identification of differentially expressed proteins during the acquisition of somatic embryogenesis in oil palm (Elaeis guineensis Jacq.). J Proteomics. 2014;104:112–27. pmid:24675181
  20. 20. Gomes HT, Bartos PMC, Silva CO, do Amaral LIV, Scherwinski-Pereira JE. Comparative biochemical profiling during the stages of acquisition and development of somatic embryogenesis in African oil palm (Elaeis guineensis Jacq.). Plant Growth Regul. 2014;74(2):199–208.
  21. 21. Pattarapimol T, Thuzar M, Vanavichit A, Tragoonrung S, Roytrakul S, Jantasuriyarat C. Identification of genes involved in somatic embryogenesis development in oil palm (Elaeis guineensis Jacq.) using cDNA AFLP. J Oil Palm Res. 2015;27(1):1–11.
  22. 22. Lin H-C, Morcillo F, Dussert S, Tranchant-Dubreuil C, Tregear JW, Tranbarger TJ. Transcriptome analysis during somatic embryogenesis of the tropical monocot Elaeis guineensis: evidence for conserved gene functions in early development. Plant Mol Biol. 2009;70(1):173–92. pmid:19199047
  23. 23. Low E-TL, Alias H, Boon S-H, Shariff EM, Tan C-YA, Ooi LCL, et al. Oil palm (Elaeis guineensis Jacq.) tissue culture ESTs: Identifying genes associated with callogenesis and embryogenesis. BMC Plant Biol. 2008;8:62–. PMC2442076. pmid:18507865
  24. 24. Silva RdC, Carmo LST, Luis ZG, Silva LP, Scherwinski-Pereira JE, Mehta A. Proteomic identification of differentially expressed proteins during the acquisition of somatic embryogenesis in oil palm (Elaeis guineensis Jacq.). J Proteomics. 2014;104:112–27. https://doi.org/10.1016/j.jprot.2014.03.013. pmid:24675181
  25. 25. Tan HS, Liddell S, Ong Abdullah M, Wong WC, Chin CF. Differential proteomic analysis of embryogenic lines in oil palm (Elaeis guineensis Jacq). J Proteomics. 2016;143:334–45. https://doi.org/10.1016/j.jprot.2016.04.039. pmid:27130535
  26. 26. Schmidt ED, Guzzo F, Toonen MA, de Vries SC. A leucine-rich repeat containing receptor-like kinase marks somatic plant cells competent to form embryos. Development. 1997;124(10):2049–62. pmid:9169851
  27. 27. Chugh A, Khurana P. Gene expression during somatic embryogenesis—Recent advances. Current Science. 2002;83(6):715–30.
  28. 28. Low E-TL, Alias H, Boon S-H, Shariff EM, Tan C-YA, Ooi LC, et al. Oil palm (Elaeis guineensis Jacq.) tissue culture ESTs: identifying genes associated with callogenesis and embryogenesis. BMC Plant Biol. 2008;8:62–. pmid:18507865.
  29. 29. Ikeda M, Umehara M, Kamada H. Embryogenesis-related genes; Its expression and roles during somatic and zygotic embryogenesis in carrot and Arabidopsis2006. 153–61 p.
  30. 30. Lopes R, Cunha RNVd, Resende MDVd. Produção de cachos e parâmetros genéticos de híbridos de caiaué com dendezeiro. Pesquisa Agropecuária Brasileira. 2012;47:1496–503.
  31. 31. Murashige T, Skoog F. A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol Plant. 1962;15(3):473–97.
  32. 32. Balzon TA, Luis ZG, Scherwinski-Pereira JE. New approaches to improve the efficiency of somatic embryogenesis in oil palm (Elaeis guineensis Jacq.) from mature zygotic embryos. In Vitro Cell Dev Biol Plant. 2013;49(1):41–50.
  33. 33. Simms D, Cizdziel PE, Chomczynski P. TRIzol: A new reagent for optimal single-step isolation of RNA. Focus. 1993;15(4):532–5.
  34. 34. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual: Cold Spring Harb Protoc.; 1989.
  35. 35. Almeida RF. Clonagem de genitores pisifera e proteômica diferencial durante a aquisição de competência embriogênica em palma de óleo (Elaeis guineensis Jacq.): Universidade de Brasilia; 2017.
  36. 36. Chan PL, Ma L-S, Low L, Shariff E, Ooi L, Cheah SC, et al. Normalized embryoid cDNA library of oil palm (Elaeis guineensis)2010. 10–1 p.
  37. 37. Chan P-L, Rose RJ, Abdul Murad AM, Zainal Z, Leslie Low E-T, Ooi LC-L, et al. Evaluation of Reference Genes for Quantitative Real-Time PCR in Oil Palm Elite Planting Materials Propagated by Tissue Culture. PLoS ONE. 2014;9(6):e99774. PMC4057393. pmid:24927412
  38. 38. Xia W, Mason AS, Xiao Y, Liu Z, Yang Y, Lei X, et al. Analysis of multiple transcriptomes of the African oil palm (Elaeis guineensis) to identify reference genes for RT-qPCR. J Biotechnol. 2014;184:63–73. https://doi.org/10.1016/j.jbiotec.2014.05.008. pmid:24862192
  39. 39. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115–e. pmid:22730293
  40. 40. Owczarzy R, Tataurov AV, Wu Y, Manthey JA, McQuisten KA, Almabrazi HG, et al. IDT SciTools: a suite for analysis and design of nucleic acid oligomers. Nucleic acids research. 2008;36(suppl_2):W163–W9.
  41. 41. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6. pmid:16081474
  42. 42. Zhao S, Fernald RD. Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comput Biol. 2005;12(8):1047–64. pmid:16241897
  43. 43. Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic acids Res. 2002;30(9):e36–e. pmid:11972351
  44. 44. Dewitte W, Murray JA. The plant cell cycle. Annu Rev Plant Biol. 2003;54(1):235–64.
  45. 45. Feher A, Pasternak TP, Dudits D. Transition of somatic plant cells to an embryogenic state. Plant Cell Tissue Organ Cult. 2003;74(3):201–28.
  46. 46. Dudits D, Györgyey J, Bögre L, Bakó L. Molecular Biology of Somatic Embryogenesis. In: Thorpe TA, editor. In Vitro Embryogenesis in Plants. Dordrecht: Springer Netherlands; 1995. p. 267–308.
  47. 47. Sharifi G, Ebrahimzadeh H, Ghareyazie B, Gharechahi J, Vatankhah E. Identification of differentially accumulated proteins associated with embryogenic and non-embryogenic calli in saffron (Crocus sativus L.). Proteome Sci. 2012;10(1):3. pmid:22243837
  48. 48. Cyrne L, us Martins L, Fernandes L, Marinho HS. Regulation of antioxidant enzymes gene expression in the yeast Saccharomyces cerevisiae during stationary phase. Free Radic Biol Med. 2003;34(3):385–93. pmid:12543254
  49. 49. Libik M, Konieczny R, Pater B, Ślesak I, Miszalski Z. Differences in the activities of some antioxidant enzymes and in H 2 O 2 content during rhizogenesis and somatic embryogenesis in callus cultures of the ice plant. Plant Cell Rep. 2005;23(12):834–41. pmid:15517278
  50. 50. Konieczny R, Libik M, Tuleja M, Niewiadomska E, Miszalski Z. Oxidative events during in vitro regeneration of sunflower. Acta Physiol Plant. 2008;30(1):71–9.
  51. 51. Araldi CG, Coelho CMM, Gaziola SA, Azevedo RA. Storage elicits a fast antioxidant enzyme activity in Araucaria angustifolia embryos. Acta Physiol Plant. 2016;38(8):201.
  52. 52. Zhou T, Yang X, Guo K, Deng J, Xu J, Gao W, et al. ROS homeostasis regulates somatic embryogenesis via the regulation of auxin signaling in cotton. Mol Cell Proteomics. 2016:mcp. M115. 049338.
  53. 53. Zhang S-g, Han S-y, Yang W-h, Wei H-l, Zhang M, Qi L-w. Changes in H2O2 content and antioxidant enzyme gene expression during the somatic embryogenesis of Larix leptolepis. Plant Cell Tissue Organ Cult. 2010;100(1):21.
  54. 54. Kairong C, Gengsheng X, Xinmin L, Gengmei X, Yafu W. Effect of hydrogen peroxide on somatic embryogenesis of Lycium barbarum L. Plant Sci. 1999;146(1):9–16.
  55. 55. Smertenko A, Franklin-Tong V. Organisation and regulation of the cytoskeleton in plant programmed cell death. Cell Death Differ. 2011;18(8):1263. pmid:21566662
  56. 56. Gourlay CW, Ayscough KR. The actin cytoskeleton: a key regulator of apoptosis and ageing? at Rev Mol Cell Biol. 2005;6(7):583.
  57. 57. Díaz-Vivancos P, Majourhat K, Fernández JA, Hernández JA, Piqueras A. Study of the antioxidant enzymatic system during shoot development from cultured intercalar meristems of saffron. Plant Growth Regul. 2011;65(1):119–26.
  58. 58. Blazquez S, Olmos E, Hernández JA, Fernández-García N, Fernández JA, Piqueras A. Somatic embryogenesis in saffron (Crocus sativus L.). Histological differentiation and implication of some components of the antioxidant enzymatic system. Plant Cell Tissue Organ Cult. 2009;97(1):49–57.
  59. 59. Godoy JA, Lunar R, Torres-Schumann S, Moreno J, Rodrigo RM, Pintor-Toro JA. Expression, tissue distribution and subcellular localization of dehydrin TAS14 in salt-stressed tomato plants. Plant Mol Biol. 1994;26(6):1921–34. pmid:7858227
  60. 60. Dong J-Z, Dunstan D. Molecular biology of somatic embryogenesis in conifers. Springer2000. p. 51–87.
  61. 61. Wilde H, Seffens W, Thomas T. Gene expression in somatic embryos. Springer1995. p. 41–52.
  62. 62. Seo PJ, Lee A-K, Xiang F, Park C-M. Molecular and functional profiling of Arabidopsis pathogenesis-related genes: insights into their roles in salt response of seed germination. Plant and Cell Physiology. 2008;49(3):334–44. pmid:18203731
  63. 63. Nolan KE, Irwanto RR, Rose RJ. Auxin up-regulates MtSERK1 expression in both Medicago truncatula root-forming and embryogenic cultures. Plant Physiology. 2003;133(1):218–30. pmid:12970488
  64. 64. Morel A, Teyssier C, Trontin JF, Eliášová K, Pešek B, Beaufour M, et al. Early molecular events involved in Pinus pinaster Ait. somatic embryo development under reduced water availability: transcriptomic and proteomic analyses. Physiol Plant. 2014;152(1):184–201. pmid:24460664
  65. 65. Sghaier-Hammami B, Drira N, Jorrín-Novo JV. Comparative 2-DE proteomic analysis of date palm (Phoenix dactylifera L.) somatic and zygotic embryos. J proteomics. 2009;73(1):161–77. pmid:19619682
  66. 66. Rolland F, Moore B, Sheen J. Sugar sensing and signaling in plants. The plant cell. 2002;14(suppl 1):S185–S205.
  67. 67. Lippert D, Zhuang J, Ralph S, Ellis DE, Gilbert M, Olafson R, et al. Proteome analysis of early somatic embryogenesis in Picea glauca. Proteomics. 2005;5(2):461–73. pmid:15627954
  68. 68. Barkla BJ, Vera-Estrella R, Hernández-Coronado M, Pantoja O. Quantitative proteomics of the tonoplast reveals a role for glycolytic enzymes in salt tolerance. The plant cell. 2009;21(12):4044–58. pmid:20028841
  69. 69. Rode C, Gallien S, Heintz D, Van Dorsselaer A, Braun H-P, Winkelmann T. Enolases: storage compounds in seeds? Evidence from a proteomic comparison of zygotic and somatic embryos of Cyclamen persicum Mill. Plant Mol Biol. 2011;75(3):305–19. pmid:21249422
  70. 70. Bauwe H, Hagemann M, Kern R, Timm S. Photorespiration has a dual origin and manifold links to central metabolism. Curr Opin Plant Biol. 2012;15(3):269–75. pmid:22284850
  71. 71. Liu Y, Chaturvedi P, Fu J, Cai Q, Weckwerth W, Yang P. Induction and quantitative proteomic analysis of cell dedifferentiation during callus formation of lotus (Nelumbo nucifera Gaertn. spp. baijianlian). J proteomics. 2016;131:61–70. pmid:26459404
  72. 72. Roowi SH, Ho C-L, Alwee SSRS, Abdullah MO, Napis S. Isolation and characterization of differentially expressed transcripts from the suspension cells of oil palm (Elaeis guineensis Jacq.) in response to different concentration of auxins. Mol Biotechnol. 2010;46(1):1–19. pmid:20390382
  73. 73. Tonietto Â, Sato JH, Teixeira JB, de Souza EM, Pedrosa FO, Franco OL, et al. Proteomic analysis of developing somatic embryos of Coffea arabica. Plant Mol Biol Report. 2012;30(6):1393–9.
  74. 74. Miyazawa Y, Sakai A, Miyagishima S-y, Takano H, Kawano S, Kuroiwa T. Auxin and cytokinin have opposite effects on amyloplast development and the expression of starch synthesis genes in cultured bright yellow-2 tobacco cells. Plant physiol. 1999;121(2):461–70. pmid:10517837
  75. 75. Kanchanapoom K, Domyoas P. The origin and development of embryoids in oil palm (Elaeis guineensis Jacq) embryo culture. Sci Asia. 1999;25:195–202.
  76. 76. Santos MO, Aragão FJ. Role of SERK genes in plant environmental response. Plant Signal Behav. 2009;4(12):1111–3. pmid:20514223
  77. 77. Shao J, Wang T, Yan X. Proteomics of pollen development and germinaion. J Proteome Res. 2007;6(12):4556–63. pmid:17958392
  78. 78. Zhao P-M, Wang L-L, Han L-B, Wang J, Yao Y, Wang H-Y, et al. Proteomic identification of differentially expressed proteins in the Ligon lintless mutant of upland cotton (Gossypium hirsutum L.). J Proteome Res. 2009;9(2):1076–87.
  79. 79. Baba AI, Nogueira FC, Pinheiro CB, Brasil JN, Jereissati ES, Juca TL, et al. Proteome analysis of secondary somatic embryogenesis in cassava (Manihot esculenta). Plant Scie. 2008;175(5):717–23.
  80. 80. Schmidt E, Guzzo F, Toonen M, De Vries S. A leucine-rich repeat containing receptor-like kinase marks somatic plant cells competent to form embryos. Development. 1997;124(10):2049–62. pmid:9169851
  81. 81. Somleva M, Schmidt E, De Vries S. Embryogenic cells in Dactylis glomerata L.(Poaceae) explants identified by cell tracking and by SERK expression. Plant Cell Rep. 2000;19(7):718–26.
  82. 82. Hecht V, Vielle-Calzada J-P, Hartog MV, Schmidt ED, Boutilier K, Grossniklaus U, et al. The Arabidopsis SOMATIC EMBRYOGENESIS RECEPTOR KINASE 1 gene is expressed in developing ovules and embryos and enhances embryogenic competence in culture. Plant Physiol. 2001;127(3):803–16. pmid:11706164
  83. 83. Thomas C, Meyer D, Himber C, Steinmetz A. Spatial expression of a sunflower SERK gene during induction of somatic embryogenesis and shoot organogenesis. Plant Physiol Biochem. 2004;42(1):35–42. pmid:15061082
  84. 84. Pan Z, Zhu S, Guan R, Deng X. Identification of 2, 4-D-responsive proteins in embryogenic callus of Valencia sweet orange (Citrus sinensis Osbeck) following osmotic stress. Plant Cell Tissue Organ Cult. 2010;103(2):145–53.
  85. 85. Youssefian S, Nakamura M, Orudgev E, Kondo N. Increased cysteine biosynthesis capacity of transgenic tobacco overexpressing an O-acetylserine (thiol) lyase modifies plant responses to oxidative stress. Plant Physiol. 2001;126(3):1001–11. pmid:11457951