Environmental regulation of sex determination in oil palm: current knowledge and insights from other species

Nonlinear impacts of climate anomalies on oil palm productivity

Oil palm contributes to various global needs as one of the most productive oil crops, but there exist ongoing concerns regarding its yield reductions and associated environmental impacts resulting from land conversion. This is the first detailed report investigating the nonlinear threats to estate-level oil palm yields posed by El Niño Southern Oscillation (ENSO) in the equatorial Pacific Ocean, a major driver of climate variability. Using the Malaysian Palm Oil Board administrative records on monthly performances reported by oil palm estates through the e-submissions portal spanning from January 2015 to June 2023, we focused on elucidating the impacts of ENSO on fresh fruit bunch yield, oil extraction rate, and oil yield. We found that both El Niño and La Niña conditions, characterized by extreme levels of ENSO indices cumulated over lags of 0-23 months prior to harvest, were associated with statistically significant reductions in yields. Lag association patterns unveiled that production risks were linked to pre-harvest exposure to extreme ENSO indices in various time windows. Subgroup analyses further revealed that the effects were pronounced in labor-intensive estates and those lacking fertilizer investments. This study underscores the necessity for adaptation strategies in response to future climate anomalies.

Development and validation of an oil palm model for a wide range of planting densities and soil textures in Malaysian growing conditions

A semi-mechanistic oil palm growth and yield model called Sawit.jl was developed to account for a wide range of planting densities and soil textures under Malaysia's climate conditions. The model comprises components related to meteorology, photosynthesis, energy balance, soil water content, and crop growth. The model simulates instantaneous meteorological properties using daily weather data, calculates simultaneous evaporation from crop and soil with the Shuttleworth-Wallace model, determines soil water content through Darcy's law, and adapts a biochemical C3 model for photosynthesis. The model is also parameterized using updated measurements from the newer tenera oil palm, including temperature-dependent Rubisco kinetics, specific leaf area, and the partitioning of nutrients and dry matter between various tree parts. Sawit.jl was validated using historical field measurement data from seven Malaysian oil palm sites, encompassing palm ages spanning 1-23 years. These seven sites differed in soil type (Inceptisols and Ultisols), planting density (82-299 palms ha-1), soil texture (27-59 % clay and 7-67 % sand), and rainfall (1800-2800 mm yr-1). The model showed overall good accuracy in simulating oil palm parameters (except for trunk weight) across diverse conditions, with model agreement metrics ranging from 6 to 27 % for model absolute errors, -22 to +17 % for model bias, and 0.38 to 0.98 for the Kling-Gupta Efficiency index. The model also predicted the response of oil palm yield to abrupt rainfall changes, such as those during El Niño and La Niña events, while accounting for how soil texture, rainfall, and other meteorological factors influence water deficits and crop photosynthesis. However, model accuracy varied by site, planting density, and oil palm parameter. Model accuracy can be increased by more accurately representing the oil palm microclimate, incorporating fruiting activity, and refining the dry matter partitioning mechanism for the trunk.

Genome-Wide Analysis of SPL/miR156 Module and Its Expression Analysis in Vegetative and Reproductive Organs of Oil Palm (Elaeis guineensis)

The SPL (SQUAMOSA-promoter binding protein-like) gene family is one of the largest plant transcription factors and is known to be involved in the regulation of plant growth, development, and stress responses. The genome-wide analysis of SPL gene members in a diverse range of crops has been elucidated. However, none of the genome-wide studies on the SPL gene family have been carried out for oil palm, an important oil-yielding plant. In this research, a total of 24 EgSPL genes were identified via a genome-wide approach. Phylogenetic analysis revealed that most of the EgSPLs are closely related to the Arabidopsis and rice SPL gene members. EgSPL genes were mapped onto the only nine chromosomes of the oil palm genome. Motif analysis revealed conservation of the SBP domain and the occurrence of 1-10 motifs in EgSPL gene members. Gene duplication analysis demonstrated the tandem duplication of SPL members in the oil palm genome. Heatmap analysis indicated the significant expression of SPL genes in shoot and flower organs of oil palm plants. Among the identified EgSPL genes, a total 14 EgSPLs were shown to be targets of miR156. Real-time PCR analysis of 14 SPL genes showed that most of the EgSPL genes were more highly expressed in female and male inflorescences of oil palm plants than in vegetative tissues. Altogether, the present study revealed the significant role of EgSPL genes in inflorescence development.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

Variability of Flowering Sex and Its Effect on Agronomic Trait Expression in White Guinea Yam

White Guinea yam (Dioscorea rotundata) is mainly a dioecious tuberous crop that produces flowers of varying sex phenotypes. Agronomic traits in Guinea yam differ according to the sex phenotype, but the precise interaction between the traits and sex phenotype is not clearly understood. This might be due to the high heterozygosity of yam where cultivars with different flowering sex have different genetic backgrounds, which mask the sole effect of sex phenotype on the agronomic traits. This study used F₁-derived clonal progenies from a bi-parental cross to minimize the impact of different genetic backgrounds among the plants with different sex phenotypes. The impact of plant sex on agronomic traits, specifically tuber yield, was evaluated through field trials conducted for four years. The results showed that only plants with a female genotype exhibited varying sex phenotypes even within the clones of same accession grown in the same experimental field. The significant effects of sex genotype and phenotype on agronomic traits were detected. Our results revealed that the flowering date was delayed in the plants with female genotypes compared to male genotypes, even when compared only among the plants with male phenotypes. The flowering date is the most important reason for the sexual differences in tuber yield. A high tuber yield was obtained when plants with the female phenotype flowered before tuber enlargement. This result can be attributed to the fact that the low flowering intensity in female plants increases the availability of carbon resources for leaf development. Female plants also showed a large negative effect of late flowering on tuber yield owing to resource competition between flowering and tuber enlargement. These findings demonstrate the feasibility of yield improvement by controlling the flowering time, with a higher effectiveness achieved in female than in male plants.

Transcriptomic dataset for early inflorescence stages of oil palm in response to defoliation stress

Oil palm breeding and seed development have been hindered due to the male parent's incapacity to produce male inflorescence as a source of pollen under normal conditions. On the other hand, a young oil palm plantation has a low pollination rate due to a lack of male flowers. These are the common problem of sex ratio in the oil palm industry. Nevertheless, the regulation of sex ratio in oil palm plants is a complex mechanism and remains an open question until now. Researchers have previously used complete defoliation to induce male inflorescences, but the biological and molecular mechanisms underlying this morphological change have yet to be discovered. Here, we present an RNA-seq dataset from three early stages of an oil palm inflorescence under normal conditions and complete defoliation stress. This transcriptomic dataset is a valuable resource to improve our understanding of sex determination mechanisms in oil palm inflorescence.

Micro-RNA-Regulated SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Expression and Cytokinin Accumulation Distinguish Early-Developing Male and Female Inflorescences in Oil Palm (Elaeis guineensis)

Sexual differentiation of inflorescences and flowers is important for reproduction and affects crop plant productivity. We report here on a molecular study of the process of sexual differentiation in the immature inflorescence of oil palm (Elaeis guineensis). This species is monoecious and exhibits gender diphasy, producing male and female inflorescences separately on the same plant in alternation. Three main approaches were used: small RNA-seq to characterise and study the expression of miRNA genes; RNA-seq to monitor mRNA accumulation patterns; hormone quantification to assess the role of cytokinins and auxins in inflorescence differentiation. Our study allowed the characterisation of 30 previously unreported palm MIRNA genes. In differential gene and miRNA expression studies, we identified a number of key developmental genes and miRNA-mRNA target modules previously described in relation to their developmental regulatory role in the cereal panicle, notably the miR156/529/535-SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) gene regulatory module. Gene enrichment analysis highlighted the importance of hormone-related genes, and this observation was corroborated by the detection of much higher levels of cytokinins in the female inflorescence. Our data illustrate the importance of branching regulation within the developmental window studied, during which the female inflorescence, unlike its male counterpart, produces flower clusters on new successive axes by sympodial growth.

Inferring the Genetic Basis of Sex Determination from the Genome of a Dioecious Nightshade

Dissecting the genetic mechanisms underlying dioecy (i.e., separate female and male individuals) is critical for understanding the evolution of this pervasive reproductive strategy. Nonetheless, the genetic basis of sex determination remains unclear in many cases, especially in systems where dioecy has arisen recently. Within the economically important plant genus Solanum (∼2,000 species), dioecy is thought to have evolved independently at least 4 times across roughly 20 species. Here, we generate the first genome sequence of a dioecious Solanum and use it to ascertain the genetic basis of sex determination in this species. We de novo assembled and annotated the genome of Solanum appendiculatum (assembly size: ∼750 Mb scaffold N50: 0.92 Mb; ∼35,000 genes), identified sex-specific sequences and their locations in the genome, and inferred that males in this species are the heterogametic sex. We also analyzed gene expression patterns in floral tissues of males and females, finding approximately 100 genes that are differentially expressed between the sexes. These analyses, together with observed patterns of gene-family evolution specific to S. appendiculatum, consistently implicate a suite of genes from the regulatory network controlling pectin degradation and modification in the expression of sex. Furthermore, the genome of a species with a relatively young sex-determination system provides the foundational resources for future studies on the independent evolution of dioecy in this clade.

Genomic diversity and genome-wide association analysis related to yield and fatty acid composition of wild American oil palm

Existing Elaeis guineensis cultivars lack sufficient genetic diversity due to extensive breeding. Harnessing variation in wild crop relatives is necessary to expand the breadth of agronomically valuable traits. Using RAD sequencing, we examine the natural diversity of wild American oil palm populations (Elaeis oleifera), a sister species of the cultivated Elaeis guineensis oil palm. We genotyped 192 wild E. oleifera palms collected from seven Latin American countries along with four cultivated E. guineensis palms. Honduras, Costa Rica, Panama and Colombia palms are panmictic and genetically similar. Genomic patterns of diversity suggest that these populations likely originated from the Amazon Basin. Despite evidence of a genetic bottleneck and high inbreeding observed in these populations, there is considerable genetic and phenotypic variation for agronomically valuable traits. Genome-wide association revealed several candidate genes associated with fatty acid composition along with vegetative and yield-related traits. These observations provide valuable insight into the geographic distribution of diversity, phenotypic variation and its genetic architecture that will guide choices of wild genotypes for crop improvement.

Comparative proteomic analysis of oil palm (Elaeis guineensis Jacq.) during early fruit development

To gain insights on protein changes in fruit setting and growth in oil palm, a comparative proteomic approach was undertaken to study proteome changes during its early development. The variations in the proteome at five early developmental stages were investigated via a gel-based proteomic technique. A total of 129 variant proteins were determined using mass spectrometric analysis, resulting in 80 identifications. The majority of the identified protein species were classified as energy and metabolism, stress response/defence and cell structure during early oil palm development representing potential candidates for the control of final fruit size and composition. Seven prominent protein species were then characterised using real-time polymerase chain reaction to validate the mRNA expression against the protein abundant profiles. Transcript and protein profiles were parallel across the developmental stages, but divergent expression was observed in one protein spot, indicative of possible post-transcriptional events. Our results revealed protein changes in early oil palm fruit development provide valuable information in the understanding of fruit growth and metabolism during early stages that may contribute towards improving agronomic traits. BIOLOGICAL SIGNIFICANCE: Two-dimensional gel electrophoresis coupled with mass spectrometry approach was used in this study to identify differentially expressed proteins during early oil palm fruit development. A total of 80 protein spots with significant change in abundance were successfully identified and selected genes were analysed using real time PCR to validate their expression. The dynamic changes in oil palm fruit proteome during early development were mostly active in primary and energy metabolism, stress responses, cell structure and protein metabolism. This study reveals the physiological processes during early oil palm fruit development and provides a reference proteome for further improvements in fruit quality traits.

A Review of Southeast Asian Oil Palm and Its CO2 Fluxes

Palm oil production is a key industry in tropical regions, driven by the demand for affordable vegetable oil. Palm oil production has been increasing by 9% every year, mostly due to expanding biofuel markets. However, the oil palm industry has been associated with key environmental issues, such as deforestation, peatland exploitation and biomass burning that release carbon dioxide (CO2) into the atmosphere, leading to climate change. This review therefore aims to discuss the characteristics of oil palm plantations and their impacts, especially CO2 emissions in the Southeast Asian region. The tropical climate and soil in Southeast Asian countries, such as Malaysia and Indonesia, are very suitable for growing oil palm trees. However, due to the scarcity of available plantation areas deforestation occurs, especially in peat swamp areas. Total carbon losses from both biomass and peat due to the conversion of tropical virgin peat swamp forest into oil palm plantations are estimated to be around 427.2 ± 90.7 t C ha−1 and 17.1 ± 3.6 t C ha−1 year−1, respectively. Even though measured CO2 fluxes have shown that overall, oil palm plantation CO2 emissions are about one to two times higher than other major crops, the ability of oil palms to absorb CO2 (a net of 64 tons of CO2 per hectare each year) and produce around 18 tons of oxygen per hectare per year is one of the main advantages of this crop. Since the oil palm industry plays a crucial role in the socio-economic development of Southeast Asian countries, sustainable and environmentally friendly practices would provide economic benefits while minimizing environmental impacts. A comprehensive review of all existing oil plantation procedures is needed to ensure that this high yielding crop has highly competitive environmental benefits.

Environmental and trophic determinism of fruit abscission and outlook with climate change in tropical regions

Fruit abscission facilitates the optimal conditions and timing of seed dispersal. Environmental regulation of tropical fruit abscission has received little attention, even though climate change may have its strongest impacts in tropical regions. In this study, oil palm fruit abscission was monitored during multiple years in the Benin Republic to take advantage of the climatic seasonality and the continuous fruit production by this species. An innovative multivariable statistical method was used to identify the best predictors of fruit abscission among a set of climate and ecophysiological variables, and the stage of inflorescence and fruit bunch development when the variables are perceived. The effects of climate scenarios on fruit abscission were then predicted based on the calibrated model. We found complex regulation takes place at specific stages of inflorescence and bunch development, even long before the fruit abscission zone is competent to execute abscission. Among the predictors selected, temperature variations during inflorescence and fruit bunch development are major determinants of the fruit abscission process. Furthermore, the timing of ripe fruit drop is determined by temperature in combination with the trophic status. Finally, climate simulations revealed that the abscission process is robust and is more affected by seasonal variations than by extreme scenarios. Our investigations highlighted the central function of the abscission zone as the sensor of environmental signals during reproductive development. Coupling ecophysiological and statistical modeling was an efficient approach to disentangle this complex environmental regulation.

Developmental basis for flower sex determination and effects of cytokinin on sex determination in Plukenetia volubilis (Euphorbiaceae)

KEY MESSAGE

Cytokinin might be an important factor to regulate floral sex at the very early stage of flower development in sacha inchi. Sacha inchi (Plukenetia volubilis, Euphorbiaceae) is characterized by having female and male flowers in a thyrse with particular differences. The mechanisms involved in the development of unisexual flowers are very poorly understood. In this study, the inflorescence and flower development of P. volubilis were investigated using light microscopy and scanning electron microscopy. We also investigated the effects of cytokinin on flower sex determination by exogenous application of 6-benzyladenine (BA) in P. volubilis. The floral development of P. volubilis was divided into eight stages, and the first morphological divergence between the male and female flowers was found to occur at stage 3. Both female and male flowers can be structurally distinguished by differences in the shape and size of the flower apex after sepal primordia initiation. There are no traces of gynoecia in male flowers or of androecia in female flowers. Exogenous application of BA effectively induced gynoecium primordia initiation and female flower development, especially at the early flower developmental stages. We propose that flower sex is determined earlier and probably occurs before flower initiation, either prior to or at inflorescence development due to the difference in the position of the female and male primordia in the inflorescence and in the time of the female and male primordia being initiated. The influence of cytokinin on female primordia during flower development in P. volubilis strongly suggests a feminization role for cytokinin in sex determination. These results indicate that cytokinin could modify the fate of the apical meristem of male flower and promote the formation of carpel primordia in P. volubilis.

Identification of sex-linked markers in the sexually cryptic coco de mer: are males and females produced in equal proportions?

Lodoicea maldivica (coco de mer) is a long-lived dioecious palm in which male and female plants are visually indistinguishable when immature, only becoming sexually dimorphic as adults, which in natural forest can take as much as 50 years. Most adult populations in the Seychelles exhibit biased sex ratios, but it is unknown whether this is due to different proportions of male and female plants being produced or to differential mortality. In this study, we developed sex-linked markers in Lodoicea using ddRAD sequencing, enabling us to reliably determine the gender of immature individuals. We screened 589 immature individuals to explore sex ratios across life stages in Lodoicea. The two sex-specific markers resulted in the amplification of male-specific bands (Lm123977 at 405 bp and Lm435135 at 130 bp). Our study of four sub-populations of Lodoicea on the islands of Praslin and Curieuse revealed that the two sexes were produced in approximately equal numbers, with no significant deviation from a 1:1 ratio before the adult stage. We conclude that sex in Lodoicea is genetically determined, suggesting that Lodoicea has a chromosomal sex determination system in which males are the heterogametic sex (XY) and females are homogametic (XX). We discuss the potential causes for observed biased sex ratios in adult populations, and the implications of our results for the life history, ecology and conservation management of Lodoicea.

Small temperature variations are a key regulator of reproductive growth and assimilate storage in oil palm (Elaeis guineensis)

Oil palm is an important crop for global vegetable oil production, and is widely grown in the humid tropical regions of Southeast Asia. Projected future climate change may well threaten palm oil production. However, oil palm plantations currently produce large amounts of unutilised biological waste. Oil palm stems – which comprise two-thirds of the waste - are especially relevant because they can contain high levels of non-structural carbohydrates (NSC) that can serve as feedstock for biorefineries. The NSC in stem are also considered a potent buffer to source-sink imbalances. In the present study, we monitored stem NSC levels and female reproductive growth. We then applied convergent cross mapping (CCM) to assess the causal relationship between the time-series. Mutual causal relationships between female reproductive growth and the stem NSC were detected, with the exception of a relationship between female reproductive organ growth and starch levels. The NSC levels were also influenced by long-term cumulative temperature, with the relationship showing a seven-month time lag. The dynamic between NSC levels and long-term cumulative rainfall showed a shorter time lag. The lower temperatures and higher cumulative rainfall observed from October to December identify this as a period with maximum stem NSC stocks.

Soil moisture regime and palm height influence embolism resistance in oil palm

With the prospect of climate change and more frequent El Niño-related dry spells, the drought tolerance of oil palm (Elaeis guineensis Jacq.), one of the most important tropical crop species, is of major concern. We studied the influence of soil water availability and palm height on the plasticity of xylem anatomy of oil palm fronds and their embolism resistance at well-drained and seasonally flooded riparian sites in lowland Sumatra, Indonesia. We found overall mean P12 and P50 values, i.e., the xylem pressures at 12% or 50% loss of hydraulic conductance, of -1.05 and - 1.86 MPa, respectively, indicating a rather vulnerable frond xylem of oil palm. This matches diurnal courses of stomatal conductance, which in combination with the observed low xylem safety evidence a sensitive water loss regulation. While the xylem anatomical traits vessel diameter (Dh), vessel density and potential hydraulic conductivity (Kp) were not different between the sites, palms in the moister riparian plots had on average by 0.4 MPa higher P50 values than plants in the well-drained plots. This could largely be attributed to differences in palm height between systems. As a consequence, palms of equal height had 1.3 MPa less negative P50 values in the moister riparian plots than in the well-drained plots. While palm height was positively related to P50, Dh and Kp decreased with height. The high plasticity in embolism resistance may be an element of the drought response strategy of oil palm, which, as a monocot, has a relatively deterministic hydraulic architecture. We conclude that oil palm fronds develop a vulnerable water transport system, which may expose the palms to increasing drought stress in a warmer and drier climate. However, the risk of hydraulic failure may be reduced by considerable plasticity in the hydraulic system and the environmental control of embolism resistance, and a presumably large stem capacitance.

Integrative omics analysis on phytohormones involved in oil palm seed germination

BACKGROUND

Heat treatment is widely used to break dormancy for seed germination and phytohormones could be deeply involved. However, effect of heat treatment on phytohormone related genes/proteins/metabolites and possible relationship with dormancy release remains unclear in oil palm. In this study, oil palm seeds were heat-treated at 39 °C for 60 days according to the method for commercial production. The embryos of seeds during heat treatment (0 d, 15 d, 30 d, 45 d and 60 d) and of germinated seeds (70 d and 75 d) were selected to discover the mechanisms involved in oil palm seed germination. RNA-seq and iTRAQ were applied to investigate DEGs and DEPs related to seed germination; qPCR and western blot were used as validation accordingly; endogenous phytohormones were determined by LC-MS/MS and exogenous phytohormones were also applied to validate their effects on seed germination.

RESULTS

RNA-seq results showed that plant hormone signal transduction was one of the most important pathways and eight phytohormones involved, while six of them (ABA, GA, ET, CTK, IAA and JA) were also identified by iTRAQ. Both RNA-seq and iTRAQ results showed that the expression of ABA decreased after heat treatment, which was further validated by qPCR and western blot. Furthermore, changes in endogenous phytohormones showed that ABA decreased rapidly to about 9% of the control at 30 d and then stayed at very low levels until germination; GA and CTK increased while IAA was not affected by heat treatment. Besides, exogenous ABA treatments (10, 100, 1000 mg/L) showed that the germination rate decreased to 63, 42 and 16% of the control, respectively, suggesting that ABA suppress seed germination and the inhibition effect increase with higher concentration; while the germination rates of exogenous GA and IAA treatments barely changed among different concentrations.

CONCLUSIONS

Phytohormones are deeply involved in oil palm seed germination and ABA acts as an inhibitor. Heat treatment can eliminate endogenous ABA and break dormancy, while GA and CTK may also involve in dormancy release. At least 30 days of heat treatment might be necessary. This study provided informative perspectives on oil palm seed germination, which could be also applicable in other palm species.

Communities of oil palm flower-visiting insects: investigating the covariation of Elaeidobius kamerunicus and other dominant species

Insects visit flowers not only to forage for nectar or pollen but also to search for hosts or prey, and to look for suitable habitats for breeding sites. In oil palm flowers, it has been documented that not all flower-visiting insects are pollinators, but some insects are recognized as predators, parasitoids or saprophages, which may affect the abundance and persistence of the weevil pollinating oil palm, Elaeidobius kamerunicus. We studied the community of oil palm flower-visiting insects and investigated the covariation between the abundance E. kamerunicus and that of other dominant species. Ecological research was conducted in oil palm plantations with different tree ages in Central Borneo. Our results found that tree age and flower type of oil palm did not influence the abundance and species richness of flower-visiting insects, but significantly affected their species composition. There was a significant positive relationship between the abundance of E. kamerunicus and the fly Scaptodrosophila sp, indicating that these species covariate in oil palm flowers. These findings suggest that understanding the covariation between E. kamerunicus and Scaptodrosophila sp may help develop the conservation strategies for E. kamerunicus to support the sustainable production of oil palm.

Development of High-Density SNP Markers and Their Application in Evaluating Genetic Diversity and Population Structure in Elaeis guineensis

High-density single nucleotide polymorphisms (SNPs) are used as highly favored makers to analyze genetic diversity and population structure, to construct high-density genetic maps and provide genotypes for genome-wide association analysis. In order to develop genome-wide SNP markers in oil palm (Elaeis guineensis), single locus amplified fragment sequencing (SLAF-seq) technology was performed in a diversity panel of 200 oil palm individuals and 1,261,501 SNPs were identified with minor allele frequency > 0.05 and integrity > 1. Among them, only 17.81% can be mapped within the genic region and the remaining was located into the intergenic region. A positive correlation was detected between the distribution of SNP markers and retrotransposons [transposable elements (TEs)]. Population structure analysis showed that the 200 individuals of oil palm can be divided into five subgroups based on cross-validation errors. However, the subpopulations divided for the 200 oil palm individuals based on the SNP markers were not accurately related to their geographical origins and 80 oil palm individuals from Malaysia showed highest genetic diversity. In addition, the physical distance of linkage disequilibrium (LD) decay in the analyzed oil palm population was 14.516 kb when r² = 0.1. The LD decay distances for different chromosomes varied from 3.324 (chromosome 15) to 19.983 kb (chromosome 7). Our research provides genome-wide SNPs for future targeted breeding in palm oil.

Transcriptome analysis of oil palm inflorescences revealed candidate genes for an auxin signaling pathway involved in parthenocarpy

Oil palm parthenocarpic fruits, which are produced without fertilization, can be targeted to increase oil content because the majority of the fruit is occupied by mesocarp, the part in which palm oil is stored. Consequently, gaining an understanding of the parthenocarpic mechanism would be instrumental for producing parthenocarpic oil palm. This study aims to determine effects of auxin treatment and analyze differentially expressed genes in oil palm pistils at the pollination/anthesis stage, using an RNA sequencing (RNA seq) approach. The auxin treatment caused 100% parthenocarpy when auxin was sprayed before stigmas opened. The parthenocarpy decreased to 55%, 8% and 5% when the auxin was sprayed 1, 2 and 3 days after the opening of stigmas, respectively. Oil palm plants used for RNA seq were plants untreated with auxin as controls and auxin-treated plants on the day before pollination and 1 day after pollination. The number of raw reads ranged from 8,425,859 to 11,811,166 reads, with an average size ranging from 99 to 137 base pairs (bp). When compared with the oil palm transcriptome, the mapped reads ranged from 8,179,948 to 11,320,799 reads, representing 95.85–98.01% of the oil palm matching. Based on five comparisons between RNA seq of treatments and controls, and confirmation using reverse transcription polymerase chain reaction and quantitative real-time RT-PCR expression, five candidate genes, including probable indole-3-acetic acid (IAA)-amido synthetase GH3.8 (EgGH3.8), IAA-amido synthetase GH3.1 (EgGH3.1), IAA induced ARG7 like (EgARG7), tryptophan amino transferase-related protein 3-like (EgTAA3) and flavin-containing monooxygenase 1 (EgFMO1), were differentially expressed between auxin-treated and untreated samples. This evidence suggests a pathway of parthenocarpic fruit development at the beginning of fruit development. However, more research is needed to identify which genes are definitely involved in parthenocarpy.

Climate Based Predictability of Oil Palm Tree Yield in Malaysia

The influence of local conditions and remote climate modes on the interannual variability of oil palm fresh fruit bunches (FFB) total yields in Malaysia and two major regions (Peninsular Malaysia and Sabah/Sarawak) is explored. On a country scale, the state of sea-surface temperatures (SST) in the tropical Pacific Ocean during the previous boreal winter is found to influence the regional climate. When El Niño occurs in the Pacific Ocean, rainfall in Malaysia reduces but air temperature increases, generating a high level of water stress for palm trees. As a result, the yearly production of FFB becomes lower than that of a normal year since the water stress during the boreal spring has an important impact on the total annual yields of FFB. Conversely, La Niña sets favorable conditions for palm trees to produce more FFB by reducing chances of water stress risk. The region of the Leeuwin current also seems to play a secondary role through the Ningaloo Niño/ Niña in the interannual variability of FFB yields. Based on these findings, a linear model is constructed and its ability to reproduce the interannual signal is assessed. This model has shown some skills in predicting the total FFB yield.

Comparative de novo transcriptome analysis of male and female Sea buckthorn

Sea buckthorn is a dioecious medicinal plant found at high altitude. The plant has both male and female reproductive organs in separate individuals. In this article, whole transcriptome de novo assemblies of male and female flower bud samples were carried out using Illumina NextSeq 500 platform to determine the role of the genes involved in sex determination. Moreover, genes with differential expression in male and female transcriptomes were identified to understand the underlying sex determination mechanism. The current study showed 63,904 and 62,272 coding sequences (CDS) in female and male transcriptome data sets, respectively. 16,831 common CDS were screened out from both transcriptomes, out of which 625 were upregulated and 491 were found to be downregulated. To understand the potential regulatory roles of differentially expressed genes in metabolic networks and biosynthetic pathways: KEGG mapping, gene ontology, and co-expression network analysis were performed. Comparison with Flowering Interactive Database (FLOR-ID) resulted in eight differentially expressed genes viz. CHD3-type chromatin-remodeling factor PICKLE (PKL), phytochrome-associated serine/threonine-protein phosphatase (FYPP), protein TOPLESS (TPL), sensitive to freezing 6 (SFR6), lysine-specific histone demethylase 1 homolog 1 (LDL1), pre-mRNA-processing-splicing factor 8A (PRP8A), sucrose synthase 4 (SUS4), ubiquitin carboxyl-terminal hydrolase 12 (UBP12), known to be broadly involved in flowering, photoperiodism, embryo development, and cold response pathways. Male and female flower bud transcriptome data of Sea buckthorn may provide comprehensive information at genomic level for the identification of genetic regulation involved in sex determination.

The association of changes in DNA methylation with temperature-dependent sex determination in cucumber

Cucumber (Cucumis sativus L.) is characterized by its diverse and flexible sexual types. Here, we evaluated the effect of low temperature (LT) exposure on cucumber femaleness under short-day conditions. Shoot apices were subjected to whole-genome bisulfate sequencing (WGBS), mRNA-seq, and sRNA-seq. The results showed that temperature had a substantial and global impact on transposable element (TE)-related small RNA-directed DNA methylation (RdDM) mechanisms, resulting in large amounts of CHH-type cytosine demethylation. In the cucumber genome, TEs are common in regions near genes that are also subject to DNA demethylation. TE-gene interactions showed very strong reactions to LT treatment, as nearly 80% of the differentially methylated regions (DMRs) were distributed in genic regions. Demethylation near genes led to the co-ordinated expression of genes and TEs. More importantly, genome-wide de novo methylation changes also resulted in small amounts of CG- and CHG-type DMRs. Methylation changes in CG-DMRs located <600 bp from the transcription start and end sites (TSSs/TESs) negatively correlated with transcription changes in differentially expressed genes (DEGs), probably indicating epiregulation. Ethylene is called the 'sex hormone' of cucumbers. We observed the up-regulation of ethylene biosynthesis-related CsACO3 and the down-regulation of an Arabidopsis RAP2.4-like ethylene-responsive (AP2/ERF) transcription factor, demonstrating the inferred epiregulation. Our study characterized the response of the apex methylome to LT and predicted the possible epiregulation of temperature-dependent sex determination (TSD) in cucumber.

Diel Patterns of Activity for Insect Pollinators of Two Oil Palm Species (Arecales : Arecaceae)

The pollination of two oil palm species, Elaeis guineensis Jacquin and Elaeis oleifera Cortés (Arecales: Arecaceae), depends on a mutualistic relation with insects, which use male inflorescences as a brood site, and visits female inflorescences lured by the emitted odor, which is similar to that of males. Although the activity of visiting the inflorescences by these insects is critical for the adequate natural pollination of the host plant, their activity is poorly documented. In the present study, we determine the diel activity of two specialized pollinator weevils (Coleoptera: Curculionidae) on inflorescences of their respective host-palm: Elaeidobius kamerunicus Faust specialized on E. guineensis, and Grasidius hybridus O'Brien and Beserra specialized on E. oleifera. The average timing of activity was studied by using passive interception traps. Then the pattern and the duration were refined by using aspiration trapping within the active period for each insect species at the male and female inflorescences. All the experiments were conducted in an Ecuadorian oil palm plantation, located close to Amazonian forest. El. kamerunicus and G. hybridus were found to be the pollinators of E. guineensis and E. oleifera, respectively. The two species differed in their diel pattern of activity: E. kamerunicus was active in the morning and G. hybridus during a short period at dusk. For both palm species, insect visits were synchronous on both male and female inflorescences. The synchronicity is discussed as a strategy to maintain the relation mutualistic between partners. These findings increase our understanding of the oil palm pollination system.

A Century of Sex Determination in Flowering Plants

Plants have evolved a diverse array of strategies for sexual reproduction, particularly through the modification of male and female organs at distinct points in development. The immense variation in sexual systems across the land plants provides a unique opportunity to study the genetic, epigenetic, phylogenetic, and ecological underpinnings of sex determination. Here, we reflect on more than a century of research into flowering plant sex determination, placing a particular focus on the foundational genetic and cytogenetic observations, experiments, and hypotheses. Building on the seminal work on the genetics of plant sex, modern comparative genomic analyses now allow us to address longstanding questions about sex determination and the origins of sex chromosomes.

The Quest for Molecular Regulation Underlying Unisexual Flower Development

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled "Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber" published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.

Comparative expression profiling of three early inflorescence stages of oil palm indicates that vegetative to reproductive phase transition of meristem is regulated by sugar balance

Breeding and seed production activities in oil palm have been hampered because of the inability of the male parent Pisifera to produce male inflorescence as source of pollen under normal conditions. Researchers are using complete defoliation to induce male inflorescences, but the biological and molecular processes responsible for this morphological change are yet to be revealed. To understand the underlying network of genes that initiate and control this phenotypically documented activity, we initiated a study aimed at identifying differentially expressed genes (DEGs) in three stages of an oil palm inflorescence under complete defoliation stress using RNA-seq. Sequencing on an Illumina platform produced 82631476 reads consisting of 8345779076 bases. A total of 60700 genes were obtained after transcript filtering and normalisation and 54% of them were downregulated. Differences in gene expression levels were significant between tissues under stress. The farther the distance between tissues, the more DEGs recorded. Comparison between stage 2 and stage 1 induced 3893 DEGs whereas 10136 DEGs were induced between stage 3 and stage 1. Stress response genes and flower development genes were among the highly expressed genes. This study suggests a link between complete defoliation and meristem differentiation from vegetative to reproductive phase in oil palm.

Oil palm natural diversity and the potential for yield improvement

African oil palm has the highest productivity amongst cultivated oleaginous crops. Species can constitute a single crop capable to fulfill the growing global demand for vegetable oils, which is estimated to reach 240 million tons by 2050. Two types of vegetable oil are extracted from the palm fruit on commercial scale. The crude palm oil and kernel palm oil have different fatty acid profiles, which increases versatility of the crop in industrial applications. Plantations of the current varieties have economic life-span around 25-30 years and produce fruits around the year. Thus, predictable annual palm oil supply enables marketing plans and adjustments in line with the economic forecasts. Oil palm cultivation is one of the most profitable land uses in the humid tropics. Oil palm fruits are the richest plant source of pro-vitamin A and vitamin E. Hence, crop both alleviates poverty, and could provide a simple practical solution to eliminate global pro-vitamin A deficiency. Oil palm is a perennial, evergreen tree adapted to cultivation in biodiversity rich equatorial land areas. The growing demand for the palm oil threatens the future of the rain forests and has a large negative impact on biodiversity. Plant science faces three major challenges to make oil palm the key element of building the future sustainable world. The global average yield of 3.5 tons of oil per hectare (t) should be raised to the full yield potential estimated at 11-18t. The tree architecture must be changed to lower labor intensity and improve mechanization of the harvest. Oil composition should be tailored to the evolving needs of the food, oleochemical and fuel industries. The release of the oil palm reference genome sequence in 2013 was the key step toward this goal. The molecular bases of agronomically important traits can be and are beginning to be understood at the single base pair resolution, enabling gene-centered breeding and engineering of this remarkable crop.

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.

Sex determination in flowering plants: papaya as a model system

Unisexuality in flowering plants evolved from a hermaphrodite ancestor. Transition from hermaphrodite to unisexual flowers has occurred multiple times across the different lineages of the angiosperms. Sexuality in plants is regulated by genetic, epigenetic and physiological mechanisms. The most specialized mechanism of sex determination is sex chromosomes. The sex chromosomes ensure the stable segregation of sexual phenotypes by preventing the recombination of sex determining genes. Despite continuous efforts, sex determining genes of dioecious plants have not yet been cloned. Concerted efforts with various model systems are necessary to understand the complex mechanism of sex determination in plants. Papaya (Carica papaya L.) is a tropical fruit tree with three sex forms, male, hermaphrodite, and female. Sexuality in papaya is determined by an XY chromosome system that is in an early evolutionary stage. The male and hermaphrodite of papaya are controlled by two different types of Y chromosomes: Y and Y(h). Large amounts of information in the area of genetics, genomics, and epigenetics of papaya have been accumulated over the last few decades. Relatively short lifecycle, small genome size, and readily available genetic and genomic resources render papaya an excellent model system to study sex determination and sex chromosomes in flowering plants.

Hormonal interactions and gene regulation can link monoecy and environmental plasticity to the evolution of dioecy in plants

Most models for dioecy in flowering plants assume that dioecy arises directly from hermaphroditism through a series of independent feminizing and masculinizing mutations that become chromosomally linked. However, dioecy appears to evolve most frequently through monoecious grades. The major genetic models do not explain the evolution of unisexual flowers in monoecious and submonoecious populations, nor do they account for environmentally induced sexual plasticity. In this review, we explore the roles of environmental stress and hormones on sex determination, and propose a model that can explain the evolution of dioecy through monoecy, and the mechanisms of environmental sex determination. Environmental stresses elicit hormones that allow plants to mediate the negative effects of the stresses. Many of these same hormones are involved in the regulation of floral developmental genes. Recent studies have elucidated the mechanisms whereby these hormones interact and can act as switchpoints in regulatory pathways. Consequently, differential concentrations of plant hormones can regulate whole developmental pathways, providing a mechanism for differential development within isogenic individuals such as seen in monoecious plants. Sex-determining genes in such systems will evolve to generate clusters of coexpressed suites. Coexpression rather than coinheritance of gender-specific genes will define the sexual developmental fate. Therefore, selection for gender type will drive evolution of the regulatory sequences of such genes rather than their synteny. Subsequent mutations to hyper- or hyposensitive alleles within the hormone response pathway can result in segregating dioecious populations. Simultaneously, such developmental systems will remain sensitive to external stimuli that modify hormone responses.

Effect of source/sink ratios on yield components, growth dynamics and structural characteristics of oil palm (Elaeis guineensis) bunches

Source/sink ratios are known to be one of the main determinants of oil palm growth and development. A long-term experiment (9 years) was conducted in Indonesia on mature oil palms subjected to continuous bunch ablation and partial defoliation treatments to artificially modify source/sink ratios. During the experiment, all harvested bunches were dissected and phenological measurements were carried out to analyse the effect of source/sink ratios on yield components explaining variations in bunch number, the number of fruits per bunch and oil dry weight per fruit. An integrative variable (supply/demand ratio) describing the ratio between the assimilate supply from sources and the growing organ demand for carbohydrate was computed for each plant on a daily basis from observations of the number of developing organs and their sink strength, and of climate variables. Defoliation and bunch ablation affected the bunch number and the fruit number per bunch. Variations in bunch number per month were mainly due to variations in the fraction of aborted inflorescence and in the ratio between female and male inflorescences. Under fluctuating trophic conditions, variations in fruit number per bunch resulted both from changes in fruit-set and in the number of branches (rachillae) per inflorescence. For defoliated plants, the decrease in the number of developing reproductive sinks appeared to be sufficient to maintain fruit weight and oil concentration at the control level, without any major decrease in the concentration of non-structural carbohydrate reserves. Computation of the supply/demand ratio revealed that each yield component had a specific phase of sensitivity to supply/demand ratios during inflorescence development. Establishing quantitative relationships between supply/demand ratios, competition and yield components is the first step towards a functional model for oil palm.