KEY MESSAGE: 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.
The mantled phenotype is an abnormal somaclonal variant arising from the oil palm cloning process and severe phenotypes lead to oil yield losses. Hypomethylation of the Karma retrotransposon within the B-type MADS-box EgDEF1 gene has been associated with this phenotype. While abnormal Karma-EgDEF1 hypomethylation was detected in mantled clones, we examined the methylation state of Karma in ortets that gave rise to high mantling rates in their clones. Small RNAs (sRNAs) were proposed to play a role in Karma hypomethylation as part of the RNA-directed DNA methylation process, hence differential expression analysis of sRNAs between the ortet groups was conducted. While no sRNA was differentially expressed at the Karma-EgDEF1 region, three sRNA clusters were differentially regulated in high-mantling ortets. The first two down-regulated clusters were possibly derived from long non-coding RNAs while the third up-regulated cluster was derived from the intron of a DnaJ chaperone gene. Several predicted mRNA targets for the first two sRNA clusters conversely displayed increased expression in high-mantling relative to low-mantling ortets. These predicted mRNA targets may be associated with defense or pathogenesis response. In addition, several differentially methylated regions (DMRs) were identified in Karma and its surrounding regions, mainly comprising subtle CHH hypomethylation in high-mantling ortets. Four of the 12 DMRs were located in a region corresponding to hypomethylated areas at the 3'end of Karma previously reported in mantled clones. Further investigations on these sRNAs and DMRs may indicate the predisposition of certain ortets towards mantled somaclonal variation.