Oil palm (Elaeis guineensis Jacq.) tissue culture ESTs: identifying genes associated with callogenesis and embryogenesis

GWAS determined genetic loci associated with callus induction in oil palm tissue culture

KEY MESSAGE

GWAS identified six loci at 25 kb downstream of WAK2, a crucial gene for cell wall and callus formation, enabling development of a SNP marker for enhanced callus induction potential. Efficient callus induction is vital for successful oil palm tissue culture, yet identifying genomic loci and markers for early detection of genotypes with high potential of callus induction remains unclear. In this study, immature male inflorescences from 198 oil palm accessions (dura, tenera and pisifera) were used as explants for tissue culture. Callus induction rates were collected at one-, two- and three-months after inoculation (C1, C2 and C3) as phenotypes. Resequencing generated 11,475,258 high quality single nucleotide polymorphisms (SNPs) as genotypes. GWAS was then performed, and correlation analysis revealed a positive association of C1 with both C2 (R = 0.81) and C3 (R = 0.50), indicating that C1 could be used as the major phenotype for callus induction rate. Therefore, only significant SNPs (P ≤ 0.05) in C1 were identified to develop markers for screening individuals with high potential of callus induction. Among 21 significant SNPs in C1, LD block analysis revealed six SNPs on chromosome 12 (Chr12) potentially linked to callus formation. Subsequently, 13 SNP markers were identified from these loci and electrophoresis results showed that marker C-12 at locus Chr12_12704856 can be used effectively to distinguish the GG allele, which showed the highest probability (69%) of callus induction. Furthermore, a rapid SNP variant detection method without electrophoresis was established via qPCR-based melting curve analysis. Our findings facilitated marker-assisted selection for specific palms with high potential of callus induction using immature male inflorescence as explant, aiding ortet palm selection in oil palm tissue culture.

Transcriptome profiling of high and low somatic embryogenesis rate of oil palm (Elaeis guineensis Jacq. var. Tenera)

Oil palm micropropagation through tissue culture is a technique to provide elite oil palms to meet the desired traits. This technique is commonly carried out through somatic embryogenesis. However, the oil palm's somatic embryogenesis rate is quite low. Several approaches have been made to overcome this problem, including transcriptome profiling through RNA-seq to identify key genes involved in oil palm somatic embryogenesis. RNA sequencing was applied in high- and low-embryogenic ortets of Tenera varieties based on the somatic embryoid rate at the callus, globular, scutellar, and coleoptilar embryoid stages. Cellular analysis of embryoid inductions and proliferations showed that high-embryogenic ortets resulted in higher embryoid proliferation and germinations than low-embryogenic ortets. Transcriptome profiling showed that there are a total of 1,911 differentially expressed genes (DEGs) between high- and low-embryogenic ortets. ABA signaling-related genes such as LEA, DDX28, and vicilin-like protein are upregulated in high-embryogenic ortets. Furthermore, DEGs associated with other hormone signaling, such as HD-ZIP associated with brassinosteroids and NPF associated with auxin, are upregulated in high-embryogenic ortets. This result suggests a physiological difference between high- and low-embryogenic ortets that is connected to their capacity for somatic embryogenesis. These DEGs will be used as potential biomarkers for high-embryogenic ortets and will be validated in further studies.

Leaf transcriptomic signatures for somatic embryogenesis potential of Elaeis guineensis

Potentially embryogenic oil palms can be identified through leaf transcriptomic signatures. Differential expression of genes involved in flowering time, and stress and light responses may associate with somatic embryogenesis potential. Clonal propagation is an attractive approach for the mass propagation of high yielding oil palms. A major issue hampering the effectiveness of oil palm tissue culture is the low somatic embryogenesis rate. Previous studies have identified numerous genes involved in oil palm somatic embryogenesis, but their association with embryogenic potential has not been determined. In this study, differential expression analysis of leaf transcriptomes from embryogenic and non-embryogenic mother palms revealed that transcriptome profiles from non- and poor embryogenic mother palms were more similar than highly embryogenic palms. A total of 171 genes exhibiting differential expression in non- and low embryogenesis groups could also discriminate high from poor embryogenesis groups of another tissue culture agency. Genes related to flowering time or transition such as FTIP, FRIGIDA-LIKE, and NF-YA were up-regulated in embryogenic ortets, suggesting that reproduction timing of the plant may associate with somatic embryogenesis potential. Several light response or photosynthesis-related genes were down-regulated in embryogenic ortets, suggesting a link between photosynthesis activity and embryogenic potential. As expression profiles of the differentially expressed genes are very similar between non- and low embryogenic groups, machine learning approaches with several candidate genes may generate a more sensitive model to better discriminate non-embryogenic from embryogenic ortets.

Early nodulin 93 protein gene: essential for induction of somatic embryogenesis in oil palm

Transcript profiling during the early induction phase of oil palm tissue culture and RNAi studies in a model somatic embryogenesis system showed that EgENOD93 expression is essential for somatic embryogenesis. Micropropagation of oil palm through tissue culture is vital for the generation of superior and uniform elite planting materials. Studies were carried out to identify genes to distinguish between leaf explants with the potential to develop into embryogenic or non-embryogenic callus. Oil palm cDNA microarrays were co-hybridized with cDNA probes of reference tissue, separately with embryo forming (media T527) and non-embryo (media T694) forming leaf explants sampled at Day 7, Day 14 and Day 21. Analysis of the normalized datasets has identified 77, 115 and 127 significantly differentially expressed genes at Day 7, Day 14, and Day 21, respectively. An early nodulin 93 protein gene (ENOD93), was highly expressed at Day 7, Day 14, and Day 21 and in callus (media T527), as assessed by RT-qPCR. Validation of EgENOD93 across tissue culture lines of different genetic background and media composition showed the potential of this gene as an embryogenic marker. In situ RNA hybridization and functional characterization in Medicago truncatula provided additional evidence that ENOD93 is essential for somatic embryogenesis. This study supports the suitability of EgENOD93 as a marker to predict the potential of leaf explants to produce embryogenic callus. Crosstalk among stresses, auxin, and Nod-factor like signalling molecules likely induces the expression of EgENOD93 for embryogenic callus formation.

Development of SSR markers based on transcriptome data and association mapping analysis for fruit shell thickness associated traits in oil palm (Elaeis guineensis Jacq.)

Present study mainly aimed to ascertain the distribution characteristics of gene-based microsatellite loci and to develop polymorphic SSR markers from the already available transcriptome data of Elaeis guineensis Jacq, an important oil crop. From this study, we identified the sum of 5791 SSRs across 51,425 unigenes from the transcripts of oil palm. We were able to evaluate 331primer pairs and characterized 183 polymorphic gene-based SSR markers. We identified a total of 506 alleles from the 183 polymorphic SSR loci, with an average of 2.77 alleles per locus. The characterized gene-based SSR markers from the transcriptome data of oil palm exhibited moderate levels of polymorphism with a significant level of heterozygosity ranges from 0.096 to 0.594 (mean = 0.336 ± 0.11). Among the identified SSR markers, sixty polymorphic markers were used to analyze genotypes of 55 oil palm accessions selected from three different provinces of China. Association mapping analysis provided the information of four markers that are associated with fruit shell thickness trait of oil palm. Among the four markers identified from association analysis, one SSR marker obtained from Unigene17150 is strictly associated with the oil palm fruit shell thickness trait.

Genome-Wide Characterization, Expression Profile Analysis of WRKY Family Genes in Santalum album and Functional Identification of Their Role in Abiotic Stress

WRKY proteins are a large superfamily of transcription factors that are involved in diverse biological processes including development, as well as biotic and abiotic stress responses in plants. WRKY family proteins have been extensively characterized and analyzed in many plant species, including Arabidopsis, rice, and poplar. However, knowledge on WRKY transcription factors in Santalum album is scarce. Based on S. album genome and transcriptome data, 64 SaWRKY genes were identified in this study. A phylogenetic analysis based on the structures of WRKY protein sequences divided these genes into three major groups (I, II, III) together with WRKY protein sequences from Arabidopsis. Tissue-specific expression patterns showed that 37 SaWRKY genes were expressed in at least one of five tissues (leaves, roots, heartwood, sapwood, or the transition zone), while the remaining four genes weakly expressed in all of these tissues. Analysis of the expression profiles of the 42 SaWRKY genes after callus was initiated by salicylic acid (SA) and methyl jasmonate (MeJA) revealed that 25 and 24 SaWRKY genes, respectively, were significantly induced. The function of SaWRKY1, which was significantly up-regulated by SA and MeJA, was analyzed. SaWRKY1 was localized in the nucleus and its overexpression improved salt tolerance in transgenic Arabidopsis. Our study provides important information to further identify the functions of SaWRKY genes and to understand the roles of SaWRKY family genes involved in the development and in SA- and MeJA-mediated stress responses.

Tissue Culture of Oil Palm: Finding the Balance Between Mass Propagation and Somaclonal Variation

The oil palm (Elaeis guineensis Jacq.) is typically propagated in vitro by indirect somatic embryogenesis, a process in which somatic cells of an explant of choice are, via an intermediate phase of callus growth, induced to differentiate into somatic embryos. The architecture of the oil palm, lacking axillary shoots, does not allow for vegetative propagation. Therefore, somatic embryogenesis is the only alternative to seed propagation, which is hampered by long germination times and low germination rates, for the production of planting material. The current oil palm somatic embryogenesis procedure is associated with several difficulties, which are described in this review. The limited availability of explants, combined with low somatic embryo initiation and regeneration rates, necessitate the proliferation of embryogenic structures, increasing the risk for somaclonal variants such as the mantled phenotype. Several ways to improve the efficiency of the tissue culture method and to reduce the risk of somaclonal variation are described. These include the use of alternative explants and propagation techniques, the introduction of specific embryo maturation treatments and the detection of the mantled abnormality in an early stage. These methods have not yet been fully explored and provide interesting research field for the future. The development of an efficient oil palm micropropagation protocol is needed to keep up with the increasing demand for palm oil in a sustainable way. Mass production of selected, high-yielding palms by tissue culture could raise yields on existing plantations, reducing the need for further expansion of the cultivated area, which is often associated with negative environmental impacts.

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

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.

Mining and comparative survey of EST-SSR markers among members of Euphorbiaceae family

Euphorbiaceae represents flowering plants family of tropical and sub-tropical region rich in secondary metabolites of economic importance. To understand and assess the genetic makeup among the members, this study was undertaken to characterize and compare SSR markers from publicly available ESTs and GSSs of nine selected species of the family. Mining of SSRs was performed by MISA, primer designing by Primer3, while functional annotation, gene ontology (GO) and enrichment analysis were performed by Blast2GO. A total 12,878 number of SSRs were detected from 101,701 number of EST sequences. SSR density ranged from 1 SSR/3.22 kb to 1 SSR/15.65 kb. A total of 1873 primer pairs were designed for the annotated SSR-Contigs. About 77.07% SSR-ESTs could be assigned a significant match to the protein database. 3037 unique SSR-FDM were assigned and IPR003657 (WRKY Domain) was found to be the most dominant FDM among the members. 1810 unique GO terms obtained were further subjected to enrichment analysis to obtain 513 statistically significant GO terms mapped to the SSR containing ESTs. Most frequent enriched GO terms were, GO:0003824 for molecular function, GO:0006350 for biological process and GO:0005886 for cellular component, justifying the richness of defensive secondary metabolites and phytomedicine within the family. The results from this study provides tangible insight to genetic make-up and distribution of SSRs. Functional annotation corresponded many genes of unknown functions which may be considered as novel genes or genes responsible for stress specific secondary metabolites. Further studies are required to understand stress specific genes accountable for leveraging the synthesis of secondary metabolites.

Codon usage and codon pair patterns in non-grass monocot genomes

BACKGROUND AND AIMS

Studies on codon usage in monocots have focused on grasses, and observed patterns of this taxon were generalized to all monocot species. Here, non-grass monocot species were analysed to investigate the differences between grass and non-grass monocots.

METHODS

First, studies of codon usage in monocots were reviewed. The current information was then extended regarding codon usage, as well as codon-pair context bias, using four completely sequenced non-grass monocot genomes (Musa acuminata, Musa balbisiana, Phoenix dactylifera and Spirodela polyrhiza) for which comparable transcriptome datasets are available. Measurements were taken regarding relative synonymous codon usage, effective number of codons, derived optimal codon and GC content and then the relationships investigated to infer the underlying evolutionary forces.

KEY RESULTS

The research identified optimal codons, rare codons and preferred codon-pair context in the non-grass monocot species studied. In contrast to the bimodal distribution of GC3 (GC content in third codon position) in grasses, non-grass monocots showed a unimodal distribution. Disproportionate use of G and C (and of A and T) in two- and four-codon amino acids detected in the analysis rules out the mutational bias hypothesis as an explanation of genomic variation in GC content. There was found to be a positive relationship between CAI (codon adaptation index; predicts the level of expression of a gene) and GC3. In addition, a strong correlation was observed between coding and genomic GC content and negative correlation of GC3 with gene length, indicating a strong impact of GC-biased gene conversion (gBGC) in shaping codon usage and nucleotide composition in non-grass monocots.

CONCLUSION

Optimal codons in these non-grass monocots show a preference for G/C in the third codon position. These results support the concept that codon usage and nucleotide composition in non-grass monocots are mainly driven by gBGC.

Genome-Wide Identification and Transferability of Microsatellite Markers between Palmae Species

The Palmae family contains 202 genera and approximately 2800 species. Except for Elaeis guineensis and Phoenix dactylifera, almost no genetic and genomic information is available for Palmae species. Therefore, this is an obstacle to the conservation and genetic assessment of Palmae species, especially those that are currently endangered. The study was performed to develop a large number of microsatellite markers which can be used for genetic analysis in different Palmae species. Based on the assembled genome of E. guineensis and P. dactylifera, a total of 814 383 and 371 629 microsatellites were identified. Among these microsatellites identified in E. guineensis, 734 509 primer pairs could be designed from the flanking sequences of these microsatellites. The majority (618 762) of these designed primer pairs had in silico products in the genome of E. guineensis. These 618 762 primer pairs were subsequently used to in silico amplify the genome of P. dactylifera. A total of 7 265 conserved microsatellites were identified between E. guineensis and P. dactylifera. One hundred and thirty-five primer pairs flanking the conserved SSRs were stochastically selected and validated to have high cross-genera transferability, varying from 16.7 to 93.3% with an average of 73.7%. These genome-wide conserved microsatellite markers will provide a useful tool for genetic assessment and conservation of different Palmae species in the future.

In vitro somatic embryogenesis and plantlet regeneration from immature male inflorescence of adult dura and tenera palms of Elaeis guineensis (Jacq.)

We report here a method for plant regeneration through somatic embryogenesis from explants collected from immature male inflorescence of adult oil palm cultivated in India. Callus induction was successful from tissues of immature male inflorescence collected from both dura and tenera varieties of oil palm. A modified Y3 (Eeuwens) media supplemented with several additives and activated charcoal (3%) were used for the experiments. Out of four different auxin treatments, 4-amino-3,5,6-trichloro-2-pyridinecarboxylic acid (picloram) produced maximum callus induction (82%) and it was not significantly different from 2,4-dichlorophenoxyacetic acid (2,4-D) and a combination of 2,4-D + picloram. The callus induction obtained with auxin α-naphthalene acetic acid was only 54% and it was significantly low as compared to the other treatments. Highest embryogenesis was obtained with a combination of 2,4-D + picloram (4.9%) followed by picloram (3.4%). Genotypic variation in response to the same auxins was observed both for callus induction and embryogenesis. Callus induction and embryogenesis ranged from 42 to 72% and 6.8 to 9.35%, respectively in tenera. The formation of embryogenic calli was marked by the appearance of white to yellowish globular or nodular structures which subsequently formed clear somatic embryos. Somatic embryogenesis was asynchronous and at one time we could find different stages of embryogenesis like the globular, torpedo and the cotyledonary stages. The somatic embryos when exposed to light in the same basal media along with 6-benzyladenine (18 µM), abscisic acid (3.78 µM) and gibberellic acid (5.78 µM) regenerated into plantlets. To the best of our knowledge this is the first report o f callus induction and somatic embryogenesis from immature male inflorescence of oil palm.

Molecular defense response of oil palm to Ganoderma infection

Basal stem rot (BSR) of oil palm roots is due to the invasion of fungal mycelia of Ganoderma species which spreads to the bole of the stem. In addition to root contact, BSR can also spread by airborne basidiospores. These fungi are able to break down cell wall components including lignin. BSR not only decreases oil yield, it also causes the stands to collapse thus causing severe economic loss to the oil palm industry. The transmission and mode of action of Ganoderma, its interactions with oil palm as a hemibiotroph, and the molecular defence responses of oil palm to the infection of Ganoderma boninense in BSR are reviewed, based on the transcript profiles of infected oil palms. The knowledge gaps that need to be filled in oil palm-Ganoderma molecular interactions i.e. the associations of hypersensitive reaction (HR)-induced cell death and reactive oxygen species (ROS) kinetics to the susceptibility of oil palm to Ganoderma spp., the interactions of phytohormones (salicylate, jasmonate and ethylene) at early and late stages of BSR, and cell wall strengthening through increased production of guaiacyl (G)-type lignin, are also discussed.

Global scale transcriptome analysis of Arabidopsis embryogenesis in vitro

Background

Somatic embryogenesis (SE) in plants is a process by which embryos are generated directly from somatic cells, rather than from the fused products of male and female gametes. Despite the detailed expression analysis of several somatic-to-embryonic marker genes, a comprehensive understanding of SE at a molecular level is still lacking. The present study was designed to generate high resolution transcriptome datasets for early SE providing the way for future research to understand the underlying molecular mechanisms that regulate this process. We sequenced Arabidopsis thaliana somatic embryos collected from three distinct developmental time-points (5, 10 and 15 d after in vitro culture) using the Illumina HiSeq 2000 platform.

Results

This study yielded a total of 426,001,826 sequence reads mapped to 26,520 genes in the A. thaliana reference genome. Analysis of embryonic cultures after 5 and 10 d showed differential expression of 1,195 genes; these included 778 genes that were more highly expressed after 5 d as compared to 10 d. Moreover, 1,718 genes were differentially expressed in embryonic cultures between 10 and 15 d. Our data also showed at least eight different expression patterns during early SE; the majority of genes are transcriptionally more active in embryos after 5 d. Comparison of transcriptomes derived from somatic embryos and leaf tissues revealed that at least 4,951 genes are transcriptionally more active in embryos than in the leaf; increased expression of genes involved in DNA cytosine methylation and histone deacetylation were noted in embryogenic tissues. In silico expression analysis based on microarray data found that approximately 5% of these genes are transcriptionally more active in somatic embryos than in actively dividing callus and non-dividing leaf tissues. Moreover, this identified 49 genes expressed at a higher level in somatic embryos than in other tissues. This included several genes with unknown function, as well as others related to oxidative and osmotic stress, and auxin signalling.

Conclusions

The transcriptome information provided here will form the foundation for future research on genetic and epigenetic control of plant embryogenesis at a molecular level. In follow-up studies, these data could be used to construct a regulatory network for SE; the genes more highly expressed in somatic embryos than in vegetative tissues can be considered as potential candidates to validate these networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1504-6) contains supplementary material, which is available to authorized users.

Exploiting transcriptome data for the development and characterization of gene-based SSR markers related to cold tolerance in oil palm (Elaeis guineensis)

BACKGROUND

The oil palm (Elaeis guineensis, 2n = 32) has the highest oil yield of any crop species, as well as comprising the richest dietary source of provitamin A. For the tropical species, the best mean growth temperature is about 27°C, with a minimal growth temperature of 15°C. Hence, the plantation area is limited into the geographical ranges of 10°N to 10°S. Enhancing cold tolerance capability will increase the total cultivation area and subsequently oil productivity of this tropical species. Developing molecular markers related to cold tolerance would be helpful for molecular breeding of cold tolerant Elaeis guineensis.

RESULTS

In total, 5791 gene-based SSRs were identified in 51,452 expressed sequences from Elaeis guineensis transcriptome data: approximately one SSR was detected per 10 expressed sequences. Of these 5791 gene-based SSRs, 916 were derived from expressed sequences up- or down-regulated at least two-fold in response to cold stress. A total of 182 polymorphic markers were developed and characterized from 442 primer pairs flanking these cold-responsive SSR repeats. The polymorphic information content (PIC) of these polymorphic SSR markers across 24 lines of Elaeis guineensis varied from 0.08 to 0.65 (mean = 0.31 ± 0.12). Using in-silico mapping, 137 (75.3%) of the 182 polymorphic SSR markers were located onto the 16 Elaeis guineensis chromosomes. Total coverage of 473 Mbp was achieved, with an average physical distance of 3.4 Mbp between adjacent markers (range 96 bp - 20.8 Mbp). Meanwhile, Comparative analysis of transcriptome under cold stress revealed that one ICE1 putative ortholog, five CBF putative orthologs, 19 NAC transcription factors and four cold-induced orhologs were up-regulated at least two fold in response to cold stress. Interestingly, 5' untranslated region of both Unigene21287 (ICE1) and CL2628.Contig1 (NAC) both contained an SSR markers.

CONCLUSIONS

In the present study, a series of SSR markers were developed based on sequences differentially expressed in response to cold stress. These EST-SSR markers would be particularly useful for gene mapping and population structure analysis in Elaeis guineensis. Meanwhile, the EST-SSR loci were inducible expressed in response to low temperature, which may have potential application in identifying trait-associated markers in oil palm in the future.

RNA-seq analysis of oil palm under cold stress reveals a different C-repeat binding factor (CBF) mediated gene expression pattern in Elaeis guineensis compared to other species

Elaeis guineensis as a tropical oil-crop is particularly sensitive to low temperature. Improvement of cold-tolerance may significantly increase the total cultivation area of this tropical oil-crop worldwide. We sequenced cold-treated and control (untreated) samples of Elaeis guineensis. De novo assembly generated 51,452 unigenes with an average length of 703 bp. Subsequently, these expressed sequences were functionally annotated. In the K category (transcription factors) of COG (Cluster of Orthologous Group) annotation, the largest proportion of genes induced and repressed at least two-fold under cold stress were from the AP2/ERE family, indicating that C-repeat binding factor, (CBFs, members of the AP2/ERE family) may play a central role in cold tolerance in Elaeis guineensis. Subsequently, the CBF-mediated signal transduction pathway was reconstructed based on transcriptome data and the gene expression profile involving the pathway was examined using real-time quantitative RT-PCR (qRT-PCR). CBFs reached maximum transcript level both at medium (4 h) and long period time points (7 days), contrary to the expression pattern of CBFs in Arabidopsis and rice. Moreover, the promoters of downstream Cold Responsive gene (CORs) regulated by CBFs were analyzed. Conservation, mutation and absence of the DRE core motif were detected in the promoters of six CORs. These mutations in DRE motifs suggest that CORs may not be induced via cold stress in Elaeis guineensis.

De novo transcriptome sequence assembly from coconut leaves and seeds with a focus on factors involved in RNA-directed DNA methylation

Coconut palm (Cocos nucifera) is a symbol of the tropics and a source of numerous edible and nonedible products of economic value. Despite its nutritional and industrial significance, coconut remains under-represented in public repositories for genomic and transcriptomic data. We report de novo transcript assembly from RNA-seq data and analysis of gene expression in seed tissues (embryo and endosperm) and leaves of a dwarf coconut variety. Assembly of 10 GB sequencing data for each tissue resulted in 58,211 total unigenes in embryo, 61,152 in endosperm, and 33,446 in leaf. Within each unigene pool, 24,857 could be annotated in embryo, 29,731 could be annotated in endosperm, and 26,064 could be annotated in leaf. A KEGG analysis identified 138, 138, and 139 pathways, respectively, in transcriptomes of embryo, endosperm, and leaf tissues. Given the extraordinarily large size of coconut seeds and the importance of small RNA-mediated epigenetic regulation during seed development in model plants, we used homology searches to identify putative homologs of factors required for RNA-directed DNA methylation in coconut. The findings suggest that RNA-directed DNA methylation is important during coconut seed development, particularly in maturing endosperm. This dataset will expand the genomics resources available for coconut and provide a foundation for more detailed analyses that may assist molecular breeding strategies aimed at improving this major tropical crop.

The use of 2D-DIGE to understand the regeneration of somatic embryos in avocado

Avocado embryogenic cell cultures can be classified into two groups based on their morphology when cultured on a medium containing auxin: somatic embryo (SE) and proembryonic masses (PEM) type cultures. The calli of SE-type cell lines are able to go through the maturation process, whereas the calli of PEM cell lines rarely mature. We have investigated four independent avocado cell cultures (two SE and two PEM). The aim of this study was to link the differential regeneration capacity of the four cell cultures to a proteomic pattern and to gain insight into the regeneration capacity. A 2D-DIGE analysis followed by a blind multivariate analysis was able to separate the two SE lines from the PEM lines indicating that the protein profiles of SE and PEM calli are different. Based on the variable importance, that is, the differential protein pattern, we hypothesize that the regeneration capacity in avocado is correlated to the ability to overcome the physicochemical stress stimuli associated with the in vitro culture. Our identical culture conditions do not seem to trigger an appropriate response in PEM lines.

Identification and computational annotation of genes differentially expressed in pulp development of Cocos nucifera L. by suppression subtractive hybridization

BACKGROUND

Coconut (Cocos nucifera L.) is one of the world's most versatile, economically important tropical crops. Little is known about the physiological and molecular basis of coconut pulp (endosperm) development and only a few coconut genes and gene product sequences are available in public databases. This study identified genes that were differentially expressed during development of coconut pulp and functionally annotated these identified genes using bioinformatics analysis.

RESULTS

Pulp from three different coconut developmental stages was collected. Four suppression subtractive hybridization (SSH) libraries were constructed (forward and reverse libraries A and B between stages 1 and 2, and C and D between stages 2 and 3), and identified sequences were computationally annotated using Blast2GO software. A total of 1272 clones were obtained for analysis from four SSH libraries with 63% showing similarity to known proteins. Pairwise comparing of stage-specific gene ontology ids from libraries B-D, A-C, B-C and A-D showed that 32 genes were continuously upregulated and seven downregulated; 28 were transiently upregulated and 23 downregulated. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis showed that 1-acyl-sn-glycerol-3-phosphate acyltransferase (LPAAT), phospholipase D, acetyl-CoA carboxylase carboxyltransferase beta subunit, 3-hydroxyisobutyryl-CoA hydrolase-like and pyruvate dehydrogenase E1 β subunit were associated with fatty acid biosynthesis or metabolism. Triose phosphate isomerase, cellulose synthase and glucan 1,3-β-glucosidase were related to carbohydrate metabolism, and phosphoenolpyruvate carboxylase was related to both fatty acid and carbohydrate metabolism. Of 737 unigenes, 103 encoded enzymes were involved in fatty acid and carbohydrate biosynthesis and metabolism, and a number of transcription factors and other interesting genes with stage-specific expression were confirmed by real-time PCR, with validation of the SSH results as high as 66.6%. Based on determination of coconut endosperm fatty acids content by gas chromatography-mass spectrometry, a number of candidate genes in fatty acid anabolism were selected for further study.

CONCLUSION

Functional annotation of genes differentially expressed in coconut pulp development helped determine the molecular basis of coconut endosperm development. The SSH method identified genes related to fatty acids, carbohydrate and secondary metabolites. The results will be important for understanding gene functions and regulatory networks in coconut fruit.

Evaluation of reference genes for quantitative real-time PCR in oil palm elite planting materials propagated by tissue culture

BACKGROUND

The somatic embryogenesis tissue culture process has been utilized to propagate high yielding oil palm. Due to the low callogenesis and embryogenesis rates, molecular studies were initiated to identify genes regulating the process, and their expression levels are usually quantified using reverse transcription quantitative real-time PCR (RT-qPCR). With the recent release of oil palm genome sequences, it is crucial to establish a proper strategy for gene analysis using RT-qPCR. Selection of the most suitable reference genes should be performed for accurate quantification of gene expression levels.

RESULTS

In this study, eight candidate reference genes selected from cDNA microarray study and literature review were evaluated comprehensively across 26 tissue culture samples using RT-qPCR. These samples were collected from two tissue culture lines and media treatments, which consisted of leaf explants cultures, callus and embryoids from consecutive developmental stages. Three statistical algorithms (geNorm, NormFinder and BestKeeper) confirmed that the expression stability of novel reference genes (pOP-EA01332, PD00380 and PD00569) outperformed classical housekeeping genes (GAPDH, NAD5, TUBULIN, UBIQUITIN and ACTIN). PD00380 and PD00569 were identified as the most stably expressed genes in total samples, MA2 and MA8 tissue culture lines. Their applicability to validate the expression profiles of a putative ethylene-responsive transcription factor 3-like gene demonstrated the importance of using the geometric mean of two genes for normalization.

CONCLUSIONS

Systematic selection of the most stably expressed reference genes for RT-qPCR was established in oil palm tissue culture samples. PD00380 and PD00569 were selected for accurate and reliable normalization of gene expression data from RT-qPCR. These data will be valuable to the research associated with the tissue culture process. Also, the method described here will facilitate the selection of appropriate reference genes in other oil palm tissues and in the expression profiling of genes relating to yield, biotic and abiotic stresses.

Analysis of multiple transcriptomes of the African oil palm (Elaeis guineensis) to identify reference genes for RT-qPCR

The African oil palm (Elaeis guineensis), which is grown in tropical and subtropical regions, is a highly productive oil-bearing crop. For gene expression-based analyses such as reverse transcription-quantitative real time PCR (RT-qPCR), reference genes are essential to provide a baseline with which to quantify relative gene expression. Normalization using reliable reference genes is critical in correctly interpreting expression data from RT-qPCR. In order to identify suitable reference genes in African oil palm, 17 transcriptomes of different tissues obtained from NCBI were systematically assessed for gene expression variation. In total, 53 putative candidate reference genes with coefficient of variation values <3.0 were identified: 18 in reproductive tissue and 35 in vegetative tissue. Analysis for enriched functions showed that approximately 90% of identified genes were clustered in cell component gene functions, and 12 out of 53 genes were traditional housekeeping genes. We selected and validated 16 reference genes chosen from leaf tissue transcriptomes by using RT-qPCR in sets of cold, drought and high salinity treated samples, and ranked expression stability using statistical algorithms geNorm, Normfinder and Bestkeeper. Genes encoding actin, adenine phosphoribosyltransferase and eukaryotic initiation factor 4A genes were the most stable genes over the cold, drought and high salinity stresses. Identification of stably expressed genes as reference gene candidates from multiple transcriptome datasets was found to be reliable and efficient, and some traditional housekeeping genes were more stably expressed than others. We provide a useful molecular genetic resource for future gene expression studies in African oil palm, facilitating molecular genetics approaches for crop improvement in this species.

Insights from computational analysis of full-length β-ketoacyl-[ACP] synthase-II cDNA isolated from American and African oil palms

Background:

Palm oil derived from fruits (mesocarp) of African oil palm (Elaeis guineensis Jacq. Tenera) and American oil palm (E. oleifera) is important for food industry. Due to high yield, Elaeis guineensis (Tenera) is cultivated on commercial scale, though its oil contains high (~54%) level of saturated fatty acids. The rate-limiting activity of beta-ketoacyl-[ACP] synthase-II (KAS-II) is considered mainly responsible for the high (44%) level of palmitic acid (C_16:0) in the oil obtained from E. guineensis.

Objective:

The objective of this study was to annotate KAS-II cDNA isolated from American and African oil palms.

Materials and Methods:

The full-length E. oleifera KAS-II (EoKAS-II) cDNA clone was isolated using random method of gene isolation. Whereas, the E. guineensis KAS-II (EgTKAS-II) cDNA was isolated using reverse transcriptase polymerase chain reaction (RT-PCR) technique; and missing ends were obtained by employing 5’and 3’ RACE technique.

Results:

The results show that EoKAS-II and EgTKAS-II open reading frames (ORFs) are of 1689 and 1721 bp in length, respectively. Further analysis of the both EoKAS-II and EgTKAS-II predicted protein illustrates that they contains conserved domains for ‘KAS-I and II’, ‘elongating’ condensing enzymes, ‘condensing enzymes super-family’, and ‘3-oxoacyl-[ACP] synthase II’. The predicted protein sequences shows 95% similarity with each other. Consecutively, the three active sites (Cys, His, and His) were identified in both proteins. However, difference in positions of two active Histidine (His) residues was noticed.

Conclusion:

These insights may serve as the foundation in understanding the variable activity of KAS-II in American and African oil palms; and cDNA clones could be useful in the genetic engineering of oil palms.

Analyses of hypomethylated oil palm gene space

Demand for palm oil has been increasing by an average of ∼8% the past decade and currently accounts for about 59% of the world's vegetable oil market. This drives the need to increase palm oil production. Nevertheless, due to the increasing need for sustainable production, it is imperative to increase productivity rather than the area cultivated. Studies on the oil palm genome are essential to help identify genes or markers that are associated with important processes or traits, such as flowering, yield and disease resistance. To achieve this, 294,115 and 150,744 sequences from the hypomethylated or gene-rich regions of Elaeis guineensis and E. oleifera genome were sequenced and assembled into contigs. An additional 16,427 shot-gun sequences and 176 bacterial artificial chromosomes (BAC) were also generated to check the quality of libraries constructed. Comparison of these sequences revealed that although the methylation-filtered libraries were sequenced at low coverage, they still tagged at least 66% of the RefSeq supported genes in the BAC and had a filtration power of at least 2.0. A total 33,752 microsatellites and 40,820 high-quality single nucleotide polymorphism (SNP) markers were identified. These represent the most comprehensive collection of microsatellites and SNPs to date and would be an important resource for genetic mapping and association studies. The gene models predicted from the assembled contigs were mined for genes of interest, and 242, 65 and 14 oil palm transcription factors, resistance genes and miRNAs were identified respectively. Examples of the transcriptional factors tagged include those associated with floral development and tissue culture, such as homeodomain proteins, MADS, Squamosa and Apetala2. The E. guineensis and E. oleifera hypomethylated sequences provide an important resource to understand the molecular mechanisms associated with important agronomic traits in oil palm.

Oil palm genome sequence reveals divergence of interfertile species in Old and New worlds

Oil palm is the most productive oil-bearing crop. Planted on only 5% of the total vegetable oil acreage, palm oil accounts for 33% of vegetable oil, and 45% of edible oil worldwide, but increased cultivation competes with dwindling rainforest reserves. We report the 1.8 gigabase (Gb) genome sequence of the African oil palm Elaeis guineensis, the predominant source of worldwide oil production. 1.535 Gb of assembled sequence and transcriptome data from 30 tissue types were used to predict at least 34,802 genes, including oil biosynthesis genes and homologues of WRINKLED1 (WRI1), and other transcriptional regulators¹, which are highly expressed in the kernel. We also report the draft sequence of the S. American oil palm Elaeis oleifera, which has the same number of chromosomes (2n=32) and produces fertile interspecific hybrids with E. guineensis², but appears to have diverged in the new world. Segmental duplications of chromosome arms define the palaeotetraploid origin of palm trees. The oil palm sequence enables the discovery of genes for important traits as well as somaclonal epigenetic alterations which restrict the use of clones in commercial plantings³, and thus helps achieve sustainability for biofuels and edible oils, reducing the rainforest footprint of this tropical plantation crop.

Extensive modulation of the transcription factor transcriptome during somatic embryogenesis in Arabidopsis thaliana

Molecular mechanisms controlling plant totipotency are largely unknown and studies on somatic embryogenesis (SE), the process through which already differentiated cells reverse their developmental program and become embryogenic, provide a unique means for deciphering molecular mechanisms controlling developmental plasticity of somatic cells. Among various factors essential for embryogenic transition of somatic cells transcription factors (TFs), crucial regulators of genetic programs, are believed to play a central role. Herein, we used quantitative real-time polymerase chain reaction (qRT-PCR) to identify TF genes affected during SE induced by in vitro culture in Arabidopsis thaliana. Expression profiles of 1,880 TFs were evaluated in the highly embryogenic Col-0 accession and the non-embryogenic tanmei/emb2757 mutant. Our study revealed 729 TFs whose expression changes during the 10-days incubation period of SE; 141 TFs displayed distinct differences in expression patterns in embryogenic versus non-embryogenic cultures. The embryo-induction stage of SE occurring during the first 5 days of culture was associated with a robust and dramatic change of the TF transcriptome characterized by the drastic up-regulation of the expression of a great majority (over 80%) of the TFs active during embryogenic culture. In contrast to SE induction, the advanced stage of embryo formation showed attenuation and stabilization of transcript levels of many TFs. In total, 519 of the SE-modulated TFs were functionally annotated and transcripts related with plant development, phytohormones and stress responses were found to be most abundant. The involvement of selected TFs in SE was verified using T-DNA insertion lines and a significantly reduced embryogenic response was found for the majority of them. This study provides comprehensive data focused on the expression of TF genes during SE and suggests directions for further research on functional genomics of SE.

Differential metabolite profiles during fruit development in high-yielding oil palm mesocarp

To better understand lipid biosynthesis in oil palm mesocarp, in particular the differences in gene regulation leading to and including de novo fatty acid biosynthesis, a multi-platform metabolomics technology was used to profile mesocarp metabolites during six critical stages of fruit development in comparatively high- and low-yielding oil palm populations. Significantly higher amino acid levels preceding lipid biosynthesis and nucleosides during lipid biosynthesis were observed in a higher yielding commercial palm population. Levels of metabolites involved in glycolysis revealed interesting divergence of flux towards glycerol-3-phosphate, while carbon utilization differences in the TCA cycle were proven by an increase in malic acid/citric acid ratio. Apart from insights into the regulation of enhanced lipid production in oil palm, these results provide potentially useful metabolite yield markers and genes of interest for use in breeding programmes.

Identification of QTLs associated with callogenesis and embryogenesis in oil palm using genetic linkage maps improved with SSR markers

Clonal reproduction of oil palm by means of tissue culture is a very inefficient process. Tissue culturability is known to be genotype dependent with some genotypes being more amenable to tissue culture than others. In this study, genetic linkage maps enriched with simple sequence repeat (SSR) markers were developed for dura (ENL48) and pisifera (ML161), the two fruit forms of oil palm, Elaeis guineensis. The SSR markers were mapped onto earlier reported parental maps based on amplified fragment length polymorphism (AFLP) and restriction fragment length polymorphism (RFLP) markers. The new linkage map of ENL48 contains 148 markers (33 AFLPs, 38 RFLPs and 77 SSRs) in 23 linkage groups (LGs), covering a total map length of 798.0 cM. The ML161 map contains 240 markers (50 AFLPs, 71 RFLPs and 119 SSRs) in 24 LGs covering a total of 1,328.1 cM. Using the improved maps, two quantitative trait loci (QTLs) associated with tissue culturability were identified each for callusing rate and embryogenesis rate. A QTL for callogenesis was identified in LGD4b of ENL48 and explained 17.5% of the phenotypic variation. For embryogenesis rate, a QTL was detected on LGP16b in ML161 and explained 20.1% of the variation. This study is the first attempt to identify QTL associated with tissue culture amenity in oil palm which is an important step towards understanding the molecular processes underlying clonal regeneration of oil palm.

Arabinogalactan proteins and the extracellular matrix surface network during peach palm somatic embryogenesis

Somatic embryogenesis has been described in peach palm as a reliable method for its in vitro multiplication and conservation. In this study, we evaluated the possible role of arabinogalactan proteins (AGPs) during this morphogenetic pathway. The presence of Yariv reagent, a synthesized chemical antibody that specifically binds AGP molecules, affected somatic embryos and callus development rate, but no effect was observed on fresh weight increment. This substance also had profound effects on embryo morphology: somatic embryos presented loose cells in the protoderm and no signs of polarization could be observed. To better evaluate the role of AGPs, analyses of specific monoclonal antibodies (MAbs) against different AGP epitopes revealed a specific pattern of distribution for each epitope. MAb JIM13 had differential expression and showed intense signal on the embryogenic sector and some immediately adjacent layers. MAb JIM7 against pectin recognized cell walls and a specific layer over the developing somatic embryo, as well as over the shoot meristem region of mature somatic embryos. This corresponds to an extracellular matrix surface network (ECMSN) associated with the development of somatic embryos and closely related to the expression of MAb JIM13. Scanning electron microscopy confirmed the presence of an ECMSN covering a specific group of cells and ultra-structural analyses revealed that the ECMSN had lipophilic substances.

Overexpression of the oil palm (Elaeis guineensis Jacq.) TAPETUM DEVELOPMENT1-like Eg707 in rice affects cell division and differentiation and reduces fertility

The functional analysis of the TAPETUM DEVELOPMENT1-like analog Eg707 of oil palm was carried out in rice by over-expressing Eg707 under the control of a double cauliflower mosaic virus 35S promoter. Ectopic expression of Eg707 in rice induced dark green and matured compact brownish calli compared to pale wild type and negative control calli. Regenerated transgenic rice plants exhibited a reduction in organ size and plant height, rolled, erect leaves, less tillers, increased chlorophyll content, and reduced fertility with smaller green seeds. At the molecular level Eg707 overexpression caused an increase in the transcription of SAPK9, a SnRK2 protein kinase family member that is activated by ABA and hyperosmotic stress. Together, the results show that ectopic Eg707 expression influences cell division and differentiation, presumably via altered hormone homeostasis.

Transcriptional profiling of genes involved in embryogenic, non-embryogenic calluses and somatic embryogenesis of Valencia sweet orange by SSH-based microarray

Somatic embryogenesis (SE) is a most promising technology that is used for in vitro germplasm conservation and genetic improvement via biotechnological approaches in citrus. Herein, three suppression subtractive hybridization (SSH) libraries were constructed using calluses of Citrus sinensis cv. 'Valencia' to explore the molecular mechanisms that underlie the SE in citrus. A total of 880 unisequences were identified by microarray screening based on these three SSH libraries. Gene ontology analysis of the differentially expressed genes indicated that nucleolus associated regulation and biogenesis processes, hormone signal transduction, and stress factors might be involved in SE. Transcription factors might also play an important role. LEC1/B3 domain regulatory network genes (LEC1, L1L, FUS3, ABI3, and ABI5) were isolated in citrus SE. Some new transcription factors associated with citrus SE, like a B3 domain containing gene and HB4, were identified. To understand the influence of these isolated genes on SE competence, their expression profiles were compared among callus lines of seven citrus cultivars with different SE competence. The expression dynamics suggested that these genes could be necessary for the SE initiation and might play a role in embryogenic competence maintenance in different cultivars. On the basis of gene expression profiles, an overview of major physiological and biosynthesis processes at different developmental stages during citrus SE is presented. For the first time, these data provide a global resource for transcriptional events important for SE in citrus, and the specific genes offer new information for further investigation on citrus SE maintenance and development.

Elaeis oleifera genomic-SSR markers: exploitation in oil palm germplasm diversity and cross-amplification in arecaceae

Species-specific simple sequence repeat (SSR) markers are favored for genetic studies and marker-assisted selection (MAS) breeding for oil palm genetic improvement. This report characterizes 20 SSR markers from an Elaeis oleifera genomic library (gSSR). Characterization of the repeat type in 2000 sequences revealed a high percentage of di-nucleotides (63.6%), followed by tri-nucleotides (24.2%). Primer pairs were successfully designed for 394 of the E. oleifera gSSRs. Subsequent analysis showed the ability of the 20 selected E. oleifera gSSR markers to reveal genetic diversity in the genus Elaeis. The average Polymorphism Information Content (PIC) value for the SSRs was 0.402, with the tri-repeats showing the highest average PIC (0.626). Low values of observed heterozygosity (H_o) (0.164) and highly positive fixation indices (F_is) in the E. oleifera germplasm collection, compared to the E. guineensis, indicated an excess of homozygosity in E. oleifera. The transferability of the markers to closely related palms, Elaeis guineensis, Cocos nucifera and ornamental palms is also reported. Sequencing the amplicons of three selected E. oleifera gSSRs across both species and palm taxa revealed variations in the repeat-units. The study showed the potential of E. oleifera gSSR markers to reveal genetic diversity in the genus Elaeis. The markers are also a valuable genetic resource for studying E. oleifera and other genus in the Arecaceae family.

Expressed sequence tags (ESTs) from young leaves of Metroxylon sagu

Sago palm, or Metroxylon sagu, is a hardy and versatile plant that is able to tolerate many stresses, biotic and abiotic, during its growth. It is one of the plants that are able to grow in waterlogged area where others could not. Apart from that sago palm is also a source of starch, contributes economically to the people and an important export for the state of Sarawak. Despite the importance of sago palm especially in the production of starch and its ability to withstand stresses, so far, not many molecular studies have been reported on sago palm. To study the characters in sago palm, transcriptome analysis was conducted where it would give a better understanding of the plant development through gene expression. Here, we report the construction of a cDNA library and preliminary expressed sequence tags analysis from the young leaves of sago palm. A total of 434 clones were sequenced with inserts ranging from 1,000 to 3,000 bps with primary and amplified titers of 8 × 10⁵ and 1.0 × 10⁹ pfu/ml, respectively. Clustering of these sequences resulted in a set of 372 tentative unigenes comprising 340 singletons and 32 contigs. The database was also annotated with BLAST2GO which showed that majority of the transcripts were involved in primary metabolism and stress tolerance.

SSR markers in transcripts of genes linked to post-transcriptional and transcriptional regulatory functions during vegetative and reproductive development of Elaeis guineensis

BACKGROUND

The oil palm (Elaeis guineensis Jacq.) is a perennial monocotyledonous tropical crop species that is now the world's number one source of edible vegetable oil, and the richest dietary source of provitamin A. While new elite genotypes from traditional breeding programs provide steady yield increases, the long selection cycle (10-12 years) and the large areas required to cultivate oil palm make genetic improvement slow and labor intensive. Molecular breeding programs have the potential to make significant impacts on the rate of genetic improvement but the limited molecular resources, in particular the lack of molecular markers for agronomic traits of interest, restrict the application of molecular breeding schemes for oil palm.

RESULTS

In the current study, 6,103 non-redundant ESTs derived from cDNA libraries of developing vegetative and reproductive tissues were annotated and searched for simple sequence repeats (SSRs). Primer pairs from sequences flanking 289 EST-SSRs were tested to detect polymorphisms in elite breeding parents and their crosses. 230 of these amplified PCR products, 88 of which were polymorphic within the breeding material tested. A detailed analysis and annotation of the EST-SSRs revealed the locations of the polymorphisms within the transcripts, and that the main functional category was related to transcription and post-transcriptional regulation. Indeed, SSR polymorphisms were found in sequences encoding AP2-like, bZIP, zinc finger, MADS-box, and NAC-like transcription factors in addition to other transcriptional regulatory proteins and several RNA interacting proteins.

CONCLUSIONS

The identification of new EST-SSRs that detect polymorphisms in elite breeding material provides tools for molecular breeding strategies. The identification of SSRs within transcripts, in particular those that encode proteins involved in transcriptional and post-transcriptional regulation, will allow insight into the functional roles of these proteins by studying the phenotypic traits that cosegregate with these markers. Finally, the oil palm EST-SSRs derived from vegetative and reproductive development will be useful for studies on the evolution of the functional diversity within the palm family.

Epigenetic imbalance and the floral developmental abnormality of the in vitro-regenerated oil palm Elaeis guineensis

BACKGROUND

The large-scale clonal propagation of oil palm (Elaeis guineensis) is being stalled by the occurrence of the mantled somaclonal variation. Indeed, this abnormality which presents a homeotic-like conversion of male floral organs into carpelloid structures, hampers oil production since the supernumerary female organs are either sterile or produce fruits with poor oil yields.

SCOPE

In the last 15 years, the prevailing point of view on the origin of the mantled floral phenotype has evolved from a random mutation event triggered by in vitro culture to a hormone-dependent dysfunction of gene regulation processes. In this review, we retrace the history of the research on the mantled variation in the light of the parallel advances made in the understanding of plant development regulation in model systems and more specifically in the role of epigenetic mechanisms. An overview of the current state of oil palm genomic and transcriptomic resources, which are key to any comparison with model organisms, is given. We show that, while displaying original characteristics, the mantled phenotype of oil palm is morphologically, and possibly molecularly, related to MADS-box genes mutants described in model plants. We also discuss the occurrence of comparable floral phenotypes in other palm species.

CONCLUSIONS

Beyond its primary interest in the search for discriminating markers against an economically crippling phenotype, the study of the mantled abnormality also provides a unique opportunity to investigate the regulation of reproductive development in a perennial tropical palm. On the basis of recent results, we propose that future efforts should concentrate on the epigenetic regulation targeting MADS-box genes and transposable elements of oil palm, since both types of sequences are most likely to be involved in the mantled variant phenotype.

Genetic linkage map of a high yielding FELDA deli×yangambi oil palm cross

Enroute to mapping QTLs for yield components in oil palm, we constructed the linkage map of a FELDA high yielding oil palm (Elaeis guineensis), hybrid cross. The parents of the mapping population are a Deli dura and a pisifera of Yangambi origin. The cross out-yielded the average by 8–21% in four trials all of which yielded comparably to the best current commercial planting materials. The higher yield derived from a higher fruit oil content. SSR markers in the public domain - from CIRAD and MPOB, as well as some developed in FELDA - were used for the mapping, augmented by locally-designed AFLP markers. The female parent linkage map comprised 317 marker loci and the male parent map 331 loci, both in 16 linkage groups each. The number of markers per group ranged from 8–47 in the former and 12–40 in the latter. The integrated map was 2,247.5 cM long and included 479 markers and 168 anchor points. The number of markers per linkage group was 15–57, the average being 29, and the average map density 4.7 cM. The linkage groups ranged in length from 77.5 cM to 223.7 cM, with an average of 137 cM. The map is currently being validated against a closely related population and also being expanded to include yield related QTLs.

A novel transcript of oil palm (Elaeis guineensis Jacq.), Eg707, is specifically upregulated in tissues related to totipotency

In this study, we report the molecular characterization of clone Eg707 isolated from cell suspension culture of the oil palm. The deduced polypeptide of clone Eg707 is highly similar to an unknown protein from Arabidopsis thaliana. The presence of an Ald-Xan-dh-C2 superfamily domain in the deduced protein sequence suggested that Eg707 protein might be involved in abscisic acid biosynthesis. Eg707 might be present as a single copy gene in the oil palm genome. This gene is highly expressed in tissue cultured materials compared to vegetative and reproductive tissues, suggesting a role of this gene during oil palm somatic embryogenesis or at the early stages of embryo development. Expression analysis of Eg707 by RNA in situ hybridization showed that Eg707 transcripts were present throughout somatic embryo development starting from proembryo formation at the embryogenic callus stages till the maturing embryo stages. Since proembryo formation within the embryogenic callus is one of the first key factors in oil palm somatic embryo development, it is suggested that Eg707 could be used as a reliable molecular marker for detecting early stage of oil palm somatic embryogenesis.

Exploiting EST databases for the development and characterization of EST-SSR markers in castor bean (Ricinus communis L.)

BACKGROUND

The castor bean (Ricinus communis L.), a monotypic species in the spurge family (Euphorbiaceae, 2n = 20), is an important non-edible oilseed crop widely cultivated in tropical, sub-tropical and temperate countries for its high economic value. Because of the high level of ricinoleic acid (over 85%) in its seed oil, the castor bean seed derivatives are often used in aviation oil, lubricants, nylon, dyes, inks, soaps, adhesive and biodiesel. Due to lack of efficient molecular markers, little is known about the population genetic diversity and the genetic relationships among castor bean germplasm. Efficient and robust molecular markers are increasingly needed for breeding and improving varieties in castor bean. The advent of modern genomics has produced large amounts of publicly available DNA sequence data. In particular, expressed sequence tags (ESTs) provide valuable resources to develop gene-associated SSR markers.

RESULTS

In total, 18,928 publicly available non-redundant castor bean EST sequences, representing approximately 17.03 Mb, were evaluated and 7732 SSR sites in 5,122 ESTs were identified by data mining. Castor bean exhibited considerably high frequency of EST-SSRs. We developed and characterized 118 polymorphic EST-SSR markers from 379 primer pairs flanking repeats by screening 24 castor bean samples collected from different countries. A total of 350 alleles were identified from 118 polymorphic SSR loci, ranging from 2-6 per locus (A) with an average of 2.97. The EST-SSR markers developed displayed moderate gene diversity (He) with an average of 0.41. Genetic relationships among 24 germplasms were investigated using the genotypes of 350 alleles, showing geographic pattern of genotypes across genetic diversity centers of castor bean.

CONCLUSION

Castor bean EST sequences exhibited considerably high frequency of SSR sites, and were rich resources for developing EST-SSR markers. These EST-SSR markers would be particularly useful for both genetic mapping and population structure analysis, facilitating breeding and crop improvement of castor bean.

Development of expressed sequence tag resources for Vanda Mimi Palmer and data mining for EST-SSR

Vanda Mimi Palmer (VMP) is a highly sought as fragrant-orchid hybrid in Malaysia. It is economically important in cosmetic and beauty industries and also a famous potted ornamental plant. To date, no work on fragrance-related genes of vandaceous orchids has been reported from other research groups although the analysis of floral fragrance or volatiles have been extensively studied. An expressed sequence tag (EST) resource was developed for VMP principally to mine any potential fragrance-related expressed sequence tag-simple sequence repeat (EST-SSR) for future development as markers in the identification of fragrant vandaceous orchids endemic to Malaysia. Clustering, annotation and assembling of the ESTs identified 1,196 unigenes which defined 966 singletons and 230 contigs. The VMP dbEST was functionally classified by gene ontology (GO) into three groups: molecular functions (51.2%), cellular components (16.4%) and biological processes (24.6%) while the remaining 7.8% showed no hits with GO identifier. A total of 112 EST-SSR (9.4%) was mined on which at least five units of di-, tri-, tetra-, penta-, or hexa-nucleotide repeats were predicted. The di-nucleotide motif repeats appeared to be the most frequent repeats among the detected SSRs with the AT/TA types as the most abundant among the dimerics, while AAG/TTC, AGA/TCT-type were the most frequent trimerics. The mined EST-SSR is believed to be useful in the development of EST-SSR markers that is applicable in the screening and characterization of fragrance-related transcripts in closely related species.

Mapping quantitative trait loci (QTLs) for fatty acid composition in an interspecific cross of oil palm

BACKGROUND

Marker Assisted Selection (MAS) is well suited to a perennial crop like oil palm, in which the economic products are not produced until several years after planting. The use of DNA markers for selection in such crops can greatly reduce the number of breeding cycles needed. With the use of DNA markers, informed decisions can be made at the nursery stage, regarding which individuals should be retained as breeding stock, which are satisfactory for agricultural production, and which should be culled. The trait associated with oil quality, measured in terms of its fatty acid composition, is an important agronomic trait that can eventually be tracked using molecular markers. This will speed up the production of new and improved oil palm planting materials.

RESULTS

A map was constructed using AFLP, RFLP and SSR markers for an interspecific cross involving a Colombian Elaeis oleifera (UP1026) and a Nigerian E. guinneensis (T128). A framework map was generated for the male parent, T128, using Joinmap ver. 4.0. In the paternal (E. guineensis) map, 252 markers (199 AFLP, 38 RFLP and 15 SSR) could be ordered in 21 linkage groups (1815 cM). Interval mapping and multiple-QTL model (MQM) mapping (also known as composite interval mapping, CIM) were used to detect quantitative trait loci (QTLs) controlling oil quality (measured in terms of iodine value and fatty acid composition). At a 5% genome-wide significance threshold level, QTLs associated with iodine value (IV), myristic acid (C14:0), palmitic acid (C16:0), palmitoleic acid (C16:1), stearic acid (C18:0), oleic acid (C18:1) and linoleic acid (C18:2) content were detected. One genomic region on Group 1 appears to be influencing IV, C14:0, C16:0, C18:0 and C18:1 content. Significant QTL for C14:0, C16:1, C18:0 and C18:1 content was detected around the same locus on Group 15, thus revealing another major locus influencing fatty acid composition in oil palm. Additional QTL for C18:0 was detected on Group 3. A minor QTL for C18:2 was detected on Group 2.

CONCLUSION

This study describes the first successful detection of QTLs for fatty acid composition in oil palm. These QTLs constitute useful tools for application in breeding programmes.

Transcriptome analysis during somatic embryogenesis of the tropical monocot Elaeis guineensis: evidence for conserved gene functions in early development

With the aim of understanding the molecular mechanisms underlying somatic embryogenesis (SE) in oil palm, we examined transcriptome changes that occur when embryogenic suspension cells are initiated to develop somatic embryos. Two reciprocal suppression subtractive hybridization (SSH) libraries were constructed from oil palm embryogenic cell suspensions: one in which embryo development was blocked by the presence of the synthetic auxin analogue 2,4-dichlorophenoxyacetic acid (2,4-D: ) in the medium (proliferation library); and another in which cells were stimulated to form embryos by the removal of 2,4-D: from the medium (initiation library). A total of 1867 Expressed Sequence Tags (ESTs) consisting of 1567 potential unigenes were assembled from the two libraries. Functional annotation indicated that 928 of the ESTs correspond to proteins that have either no similarity to sequences in public databases or are of unknown function. Gene Ontology (GO) terms assigned to the two EST populations give clues to the underlying molecular functions, biological processes and cellular components involved in the initiation of embryo development. Macroarrays were used for transcript profiling the ESTs during SE. Hierarchical cluster analysis of differential transcript accumulation revealed 4 distinct profiles containing a total of 192 statistically significant developmentally regulated transcripts. Similarities and differences between the global results obtained with in vitro systems from dicots, monocots and gymnosperms will be discussed.