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

Multi-Omics Approaches in Oil Palm Research: A Comprehensive Review of Metabolomics, Proteomics, and Transcriptomics Based on Low-Temperature Stress

Oil palm (Elaeis guineensis Jacq.) is a typical tropical oil crop with a temperature of 26-28 °C, providing approximately 35% of the total world's vegetable oil. Growth and productivity are significantly affected by low-temperature stress, resulting in inhibited growth and substantial yield losses. To comprehend the intricate molecular mechanisms underlying the response and acclimation of oil palm under low-temperature stress, multi-omics approaches, including metabolomics, proteomics, and transcriptomics, have emerged as powerful tools. This comprehensive review aims to provide an in-depth analysis of recent advancements in multi-omics studies on oil palm under low-temperature stress, including the key findings from omics-based research, highlighting changes in metabolite profiles, protein expression, and gene transcription, as well as including the potential of integrating multi-omics data to reveal novel insights into the molecular networks and regulatory pathways involved in the response to low-temperature stress. This review also emphasizes the challenges and prospects of multi-omics approaches in oil palm research, providing a roadmap for future investigations. Overall, a better understanding of the molecular basis of the response of oil palm to low-temperature stress will facilitate the development of effective breeding and biotechnological strategies to improve the crop's resilience and productivity in changing climate scenarios.

Towards DNA-free CRISPR/Cas9 genome editing for sustainable oil palm improvement

The CRISPR/Cas9 genome editing system has been in the spotlight compared to programmable nucleases such as ZFNs and TALENs due to its simplicity, versatility, and high efficiency. CRISPR/Cas9 has revolutionized plant genetic engineering and is broadly used to edit various plants' genomes, including those transformation-recalcitrant species such as oil palm. This review will comprehensively present the CRISPR-Cas9 system's brief history and underlying mechanisms. We then highlighted the establishment of the CRISPR/Cas9 system in plants with an emphasis on the strategies of highly efficient guide RNA design, the establishment of various CRISPR/Cas9 vector systems, approaches of multiplex editing, methods of transformation for stable and transient techniques, available methods for detecting and analyzing mutations, which have been applied and could be adopted for CRISPR/Cas9 genome editing in oil palm. In addition, we also provide insight into the strategy of DNA-free genome editing and its potential application in oil palm.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Crocus genome reveals the evolutionary origin of crocin biosynthesis

Crocus sativus (saffron) is a globally autumn-flowering plant, and its stigmas are the most expensive spice and valuable herb medicine. Crocus specialized metabolites, crocins, are biosynthesized in distant species, Gardenia (eudicot) and Crocus (monocot), and the evolution of crocin biosynthesis remains poorly understood. With the chromosome-level Crocus genome assembly, we revealed that two rounds of lineage-specific whole genome triplication occurred, contributing important roles in the production of carotenoids and apocarotenoids. According to the kingdom-wide identification, phylogenetic analysis, and functional assays of carotenoid cleavage dioxygenases (CCDs), we deduced that the duplication, site positive selection, and neofunctionalization of Crocus-specific CCD2 from CCD1 members are responsible for the crocin biosynthesis. In addition, site mutation of CsCCD2 revealed the key amino acids, including I143, L146, R161, E181, T259, and S292 related to the catalytic activity of zeaxanthin cleavage. Our study provides important insights into the origin and evolution of plant specialized metabolites, which are derived by duplication events of biosynthetic genes.

EgMADS3 directly regulates EgLPAAT to mediate medium-chain fatty acids (MCFA) anabolism in the mesocarp of oil palm

KEY MESSAGE

EgMADS3, a pivotal transcription factor, positively regulates MCFA accumulation via binding to the EgLPAAT promoter, advancing lipid content in mesocarp of oil palm. Lipids function as the structural components of cell membranes, which serve as permeable barriers to the external environment of cells. The medium-chain fatty acid in the stored lipids of plants is an important renewable energy. Most research on MCFA production in plant lipid synthesis is based on biochemical methods, and the importance of transcriptional regulation in MCFA synthesis and its incorporation into TAGs needs further research. Oil palm is the most productive oil crop in the world and has the highest productivity among the main oil crops. In this study, the MADS transcription factor (EgMADS3) in the mesocarp of oil palm was characterized. Through the VIGS-virus induced gene silencing, it was determined that the potential target gene of EgMADS3 was related to the biosynthesis of medium-chain fatty acid (MCFA). Transient transformation in protoplasts and qRT-PCR analysis showed that EgMADS3 positively regulated the expression of EgLPAAT. The results of the yeast one-hybrid assays and EMSA indicated the interaction between EgMADS3 and EgLPAAT promoter. Through genetic transformation and fatty acid analysis, it is concluded that EgMADS3 directly regulates the mid-chain fatty acid synthesis pathway of the potential target gene EgLPAAT, thus promotes the accumulation of MCFA and improves the total lipid content. This study is innovative in the functional analysis of the MADS family transcription factor in the metabolism of medium-chain fatty acids (MCFA) of oil palm, provides a certain research basis for improving the metabolic pathway of chain fatty acids in oil palm, and improves the synthesis of MCFA in plants. Our results will provide a reference direction for further research on improving the oil quality through biotechnology of oil palm.

De novo transcriptomic analysis of Doum Palm (Hyphaene compressa) revealed an insight into its potential drought tolerance

BACKGROUND

Doum palms (Hyphaene compressa) perform a crucial starring role in the lives of Kenya's arid and semi-arid people for empowerment and sustenance. Despite the crop's potential for economic gain, there is a lack of genetic resources and detailed information about its domestication at the molecular level. Given the doum palm's vast potential as a widely distributed plant in semi-arid and arid climates and a source of many applications, coupled with the current changing climate scenario, it is essential to understand the molecular processes that provide drought resistance to this plant.

RESULTS

Assembly of the first transcriptome of doum palms subjected to water stress generated about 39.97 Gb of RNA-Seq data. The assembled transcriptome revealed 193,167 unigenes with an average length of 1655 bp, with 128,708 (66.63%) successfully annotated in seven public databases. Unigenes exhibited significant differentially expressed genes (DEGs) in well-watered and stressed-treated plants, with 45071 and 42457 accounting for up-regulated and down-regulated DEGs, respectively. GO term, KEGG, and KOG analysis showed that DEGs were functionally enriched cellular processes, metabolic processes, cellular and catalytic activity, metabolism, genetic information processing, signal transduction mechanisms, and posttranslational modification pathways. Transcription factors (TF), such as the MYB, WRKY, NAC family, FAR1, B3, bHLH, and bZIP, were the prominent TF families identified as doum palm DEGs encoding drought stress tolerance.

CONCLUSIONS

This study provides a complete understanding of DEGs involved in drought stress at the transcriptome level in doum palms. This research is, therefore, the foundation for the characterization of potential genes, leading to a clear understanding of its drought stress responses and providing resources for improved genetic modification.

Genotyping by sequencing for the construction of oil palm (Elaeis guineensis Jacq.) genetic linkage map and mapping of yield related quantitative trait loci

Background

Oil palm (Elaeis guineensis Jacq.) is one of the major oil-producing crops. Improving the quality and increasing the production yield of oil palm have been the primary focuses of both conventional and modern breeding approaches. However, the conventional breeding approach for oil palm is very challenging due to its longevity, which results in a long breeding cycle. Thus, the establishment of marker assisted selection (MAS) for oil palm breeding programs would speed up the breeding pipeline by generating new oil palm varieties that possess high commercial traits. With the decreasing cost of sequencing, Genotyping-by-sequencing (GBS) is currently feasible to many researchers and it provides a platform to accelerate the discovery of single nucleotide polymorphism (SNP) as well as insertion and deletion (InDel) markers for the construction of a genetic linkage map. A genetic linkage map facilitates the identification of significant DNA regions associated with the trait of interest via quantitative trait loci (QTL) analysis.

Methods

A mapping population of 112 F1 individuals from a cross of Deli dura and Serdang pisifera was used in this study. GBS libraries were constructed using the double digestion method with HindIII and TaqI enzymes. Reduced representation libraries (RRL) of 112 F1 progeny and their parents were sequenced and the reads were mapped against the E. guineensis reference genome. To construct the oil palm genetic linkage map, informative SNP and InDel markers were used to discover significant DNA regions associated with the traits of interest. The nine traits of interest in this study were fresh fruit bunch (FFB) yield, oil yield (OY), oil to bunch ratio (O/B), oil to dry mesocarp ratio (O/DM) ratio, oil to wet mesocarp ratio (O/WM), mesocarp to fruit ratio (M/F), kernel to fruit ratio (K/F), shell to fruit ratio (S/F), and fruit to bunch ratio (F/B).

Results

A total of 2.5 million SNP and 153,547 InDel markers were identified. However, only a subset of 5,278 markers comprising of 4,838 SNPs and 440 InDels were informative for the construction of a genetic linkage map. Sixteen linkage groups were produced, spanning 2,737.6 cM for the maternal map and 4,571.6 cM for the paternal map, with average marker densities of one marker per 2.9 cM and one per 2.0 cM respectively, were produced. A QTL analysis was performed on nine traits; however, only QTL regions linked to M/F, K/F and S/F were declared to be significant. Of those QTLs were detected: two for M/F, four for K/F and one for S/F. These QTLs explained 18.1-25.6% of the phenotypic variance and were located near putative genes, such as casein kinase II and the zinc finger CCCH domain, which are involved in seed germination and growth. The identified QTL regions for M/F, K/F and S/F from this study could be applied in an oil palm breeding program and used to screen palms with desired traits via marker assisted selection (MAS).

Not just females and males: Unravelling the complex sex determinism of the hemp palm, Trachycarpus fortunei

PREMISE

The ornamental Asian palm Trachycarpus fortunei (Arecaceae: Coryphoideae) is widely planted in temperate regions. In Europe, it has spread outside of gardens, particularly on the southern side of the Alps. Sexual expression in the species is complex, varying from dioecy to polygamy. This study investigated (1) sexual floral development and (2) genetic markers implicated in sex determinism.

METHODS

The morphology and anatomy of floral organs at different developmental stages were studied using SEM observations and anatomical section. Sex determinism was explored using a genome-wide association study approach, searching for correlations between 31,000 single-nucleotide polymorphisms and sex affiliation of 122 palms from 21 wild populations.

RESULTS

We observed that sexual differentiation appears late in floral development of T. fortunei. Morpho-anatomical characters of flowers conducive to panmixia were observed, such as well-differentiated septal nectaries that are thought to promote cross-pollination. At the molecular level, homozygous and heterozygous allelic systems with closely linked regions were found for sex determinism in individuals with female and "dominant-male" phenotypes, respectively. Through our wide sampling in the southern Alps, the closely linked genetic regions in males suggest that at least fifteen percent of wild palms are the direct offspring of "males" that can also produce fertile pistillate flowers.

CONCLUSIONS

Trachycarpus fortunei is a further example of unstable sexual expression found in the family Arecaceae and represents an evolutionary path towards an XY genetic system. Our structural and genetic results may explain the high species dispersal ability in the southern Alps.

Identification and Analysis of MADS-Box Genes Expressed in the Mesocarp of Oil Palm Fruit (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis) is the most important tropical oil-bearing crop species worldwide. MADS-box proteins, which play crucial roles in plant growth and development and are involved in various physiological and biochemical processes, compose one of the largest families of plant transcription factors. In this study, 42 MADS-box genes were screened from the mesocarp transcriptome database of oil palm fruit, and their phylogenetic relationships with Arabidopsis thaliana MADS-box genes were analyzed. Based on the results, MADS-box genes from oil palm mesocarp were classified into four groups: MIKCc-type, MIKC*-type, Mα-type, and Mγ-type MADS-box genes. Members of the subfamilies were classified according to the presence of three specific protein motifs. To explore the differential expression of the MADS-box genes, the dynamic expression of all selected MADS-box genes in oil palm was measured by RNA-seq. The high expression of specific MADS-box genes in the mesocarp of oil palm during different developmental stages indicates that those genes may play important roles in the cell division of and metabolite accumulation in the fruit and could become important targets for fruit development and oil accumulation research in oil palm.

Molecular insights into the VIRESCENS amino acid sequence and its implication in anthocyanin production in red- and yellow-fruited cultivars of date palm

This study examined the amino acid sequence of the VIRESCENS gene (VIR), which regulates the production of anthocyanin in 12 cultivars of the date palm (Phoenix dactylifera L.), grown in Al-Madinah Al-Munawarah of the Kingdom of Saudi Arabia. The gene products were amplified via polymerase chain reactions, amplifying both exons and introns. The products were sequenced for the reconstruction of a phylogenetic tree, which used the associated amino acid sequences. The ripening stages of Khalal, Rutab, and Tamar varied among the cultivars. Regarding VIR genotype, the red date had the wild-type gene (VIR+), while the yellow date carried a dominant mutation (VIRIM), i.e., long terminal repeat retrotransposons (LTR-RTs). The DNA sequence of VIRIM revealed that the insertion length of the LTR-RTs ranged between 386 and 476 bp. The R2 and R3 motifs in both VIR+ and VIRIM were conserved. The C-terminus motifs S6A, S6B, and S6C were found in the VIR+ protein sequence. However, the amino acids at positions 123, 161, 166, and 168 differed between VIR+ and VIRIM, and were not included in the C-terminus motifs. Within the VIR+ allele, the lysine at position 187 in the C-terminus was located immediately after S6B, with a protein binding score of 0.3, which was unique to the dark, red-fruited cultivars Ajwah, Anbarah, and Safawi. In the lighter, red-fruited cultivars, the presence of glutamic acid at the same position suggested that the anthocyanin regulation of date palm might be outside the R2 and R3 domains in the N-terminus.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Molecular characterization of oleosin genes in Cyperus esculentus, a Cyperaceae plant producing oil in underground tubers

KEY MESSAGE

CeOLE genes exhibit a tuber-predominant expression pattern and their mRNA/protein abundances are positively correlated with oil accumulation during tuber development. Overexpression could significantly increase the oil content of tobacco leaves. Oleosins (OLEs) are abundant structural proteins of lipid droplets (LDs) that function in LD formation and stabilization in seeds of oil crops. However, little information is available on their roles in vegetative tissues. In this study, we present the first genome-wide characterization of the oleosin family in tigernut (Cyperus esculentus L., Cyperaceae), a rare example accumulating high amounts of oil in underground tubers. Six members identified represent three previously defined clades (i.e. U, SL and SH) or six out of seven orthogroups (i.e. U, SL1, SL2, and SH1-3) proposed in this study. Comparative genomics analysis reveals that lineage-specific expansion of Clades SL and SH was contributed by whole-genome duplication and dispersed duplication, respectively. Moreover, presence of SL2 and SH3 in Juncus effuses implies their appearance sometime before Cyperaceae-Juncaceae divergence, whereas SH2 appears to be Cyperaceae specific. Expression analysis showed that CeOLE genes exhibit a tuber-predominant expression pattern and transcript levels are considerably more abundant than homologs in the close relative Cyperus rotundus. Moreover, CeOLE mRNA and protein abundances were shown to positively correlate with oil accumulation during tuber development. Additionally, two dominant isoforms (i.e. CeOLE2 and -5) were shown to locate in LDs as well as the endoplasmic reticulum of tobacco (Nicotiana benthamiana) leaves, and are more likely to function in homo and heteromultimers. Furthermore, overexpression of CeOLE2 and -5 in tobacco leaves could significantly increase the oil content, supporting their roles in oil accumulation. These findings provide insights into lineage-specific family evolution and putative roles of CeOLE genes in oil accumulation of vegetative tissues, which facilitate further genetic improvement for tigernut.

Functional characterization of the MSP-C6 promoter as a potential tool for mesocarp-preferential expression of transgenes

Modification of lipid composition in the mesocarp tissue of oil palm involves genetic manipulation of multiple genes. More than one mesocarp-preferential promoter is necessary for the expression of individual transgenes in the same plant to obviate transcriptional gene silencing. This study aimed to identify genes that are preferentially expressed in the mesocarp tissue and characterize selected candidate mesocarp-preferential promoters. Ten transcripts that were preferentially expressed in the mesocarp tissue were identified from the analysis of 82 transcriptome datasets of 12 different oil palm tissues. The expression of two candidate genes, MSP-C1 and MSP-C6, was verified to be preferentially expressed in the mesocarp tissues and shown to have a low expression level in non-mesocarp tissues by reverse transcription quantitative real-time PCR (RT-qPCR). MSP-C6 promoter fragments of different lengths were transformed into tomato plants for further characterization. Both unripe and ripe fruits of transgenic tomato plants transformed with a construct harboring the MSP-C6-F1 (2014 bp) promoter were shown to have high beta-glucuronidase (GUS) activities. The findings of this study suggest the potential applications of the MSP-C6 promoter as a molecular tool for genetic engineering of novel traits in fruit crops.

Identification of oil palm cis-regulatory elements based on DNA free energy and single nucleotide polymorphism density

Transcription control through cis-regulatory elements (CREs) is one of important regulators of gene expression. This study aimed to identify the location of CREs in oil palm (Elaeis guineensis Jacq.) using the combination of DNA free energy and single nucleotide polymorphism (SNP) density approaches. Promoter region sequences were extracted oil palm genome spanning from 1500 nucleotides (nt) upstream to 1000 nt downstream of every annotated transcription start sites (TSS). Free energy profiles of each promoter region were calculated using PromPredict software. Raw reads from the deep sequencing of 59 oil palm origins were used to calculate SNP density of each promoter region. The result showed that the average free energy (AFE) on the upstream region of TSS is about 1.5 kcal/mol higher compared to the downstream region. Using DNA free energy method, 16,281 regions of CREs were predicted. Most of predicted CREs was located between 1 and 500 nt upstream of TSS. Anti-correlation pattern between free energy and SNP density was observed on the predicted regions of CREs. This anti-correlated pattern was also observed on an experimentally determined promoter of the oil palm metallothionein gene, EgMSP1. Considering the increasing use of promoter information on plant biotechnology, an easy and accurate promoter prediction using the combination of free energy and SNP density method could be recommended.

Interaction between phenylpropane metabolism and oil accumulation in the developing seed of Brassica napus revealed by high temporal-resolution transcriptomes

BACKGROUND

Brassica napus is an important oilseed crop providing high-quality vegetable oils for human consumption and non-food applications. However, the regulation between embryo and seed coat for the synthesis of oil and phenylpropanoid compounds remains largely unclear.

RESULTS

Here, we analyzed the transcriptomes in developing seeds at 2-day intervals from 14 days after flowering (DAF) to 64 DAF. The 26 high-resolution time-course transcriptomes are clearly clustered into five distinct groups from stage I to stage V. A total of 2217 genes including 136 transcription factors, are specifically expressed in the seed and show high temporal specificity by being expressed only at certain stages of seed development. Furthermore, we analyzed the co-expression networks during seed development, which mainly included master regulatory transcription factors, lipid, and phenylpropane metabolism genes. The results show that the phenylpropane pathway is prominent during seed development, and the key enzymes in the phenylpropane metabolic pathway, including TT5, BAN, and the transporter TT19, were directly or indirectly related to many key enzymes and transcription factors involved in oil accumulation. We identified candidate genes that may regulate seed oil content based on the co-expression network analysis combined with correlation analysis of the gene expression with seed oil content and seed coat content.

CONCLUSIONS

Overall, these results reveal the transcriptional regulation between lipid and phenylpropane accumulation during B. napus seed development. The established co-expression networks and predicted key factors provide important resources for future studies to reveal the genetic control of oil accumulation in B. napus seeds.

High-resolution genetic linkage map and height-related QTLs in an oil palm (Elaeis guineensis) family planted across multiple sites

A refined SNP array containing 92,459 probes was developed and applied for chromosome scanning, construction of a high-density genetic linkage map and QTL analysis in a selfed Nigerian oil palm family (T128). Genotyping of the T128 mapping family generated 76,447 good quality SNPs for detailed scanning of aberration and homozygosity in the individual pseudo-chromosomes. Of them, 25,364 polymorphic SNPs were used for linkage analysis resulting in an 84.4% mapping rate. A total of 21,413 SNPs were mapped into 16 linkage groups (LGs), covering a total map length of 1364.5 cM. This genetic map is 16X denser than the previous version used to establish pseudo-chromosomes of the oil palm reference genome published in 2013. The QTLs associated with height, height increment and rachis length were identified in LGs TT05, 06, 08, 15 and 16. The present QTLs as well as those published previously were tagged to the reference genome to determine their chromosomal locations. Almost all the QTLs identified in this study were either close to or co-located with those reported in other populations. Determining the QTL position on chromosomes was also helpful in mining for the underlying candidate genes. In total, 55 putative genes and transcription factors involved in the biosynthesis, conjugation and signalling of the major phytohormones, especially for gibberellins and cell wall morphogenesis were found to be present in the identified genomic QTL regions, and their potential roles in plant dwarfism are discussed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01360-2.

Molecular Interplay between Non-Host Resistance, Pathogens and Basal Immunity as a Background for Fatal Yellowing in Oil Palm (Elaeis guineensis Jacq.) Plants

An oil palm (Elaeis guineensis Jacq.) bud rod disorder of unknown etiology, named Fatal Yellowing (FY) disease, is regarded as one of the top constraints with respect to the growth of the palm oil industry in Brazil. FY etiology has been a challenge embraced by several research groups in plant pathology throughout the last 50 years in Brazil, with no success in completing Koch's postulates. Most recently, the hypothesis of having an abiotic stressor as the initial cause of FY has gained ground, and oxygen deficiency (hypoxia) damaging the root system has become a candidate for stress. Here, a comprehensive, large-scale, single- and multi-omics integration analysis of the metabolome and transcriptome profiles on the leaves of oil palm plants contrasting in terms of FY symptomatology-asymptomatic and symptomatic-and collected in two distinct seasons-dry and rainy-is reported. The changes observed in the physicochemical attributes of the soil and the chemical attributes and metabolome profiles of the leaves did not allow the discrimination of plants which were asymptomatic or symptomatic for this disease, not even in the rainy season, when the soil became waterlogged. However, the multi-omics integration analysis of enzymes and metabolites differentially expressed in asymptomatic and/or symptomatic plants in the rainy season compared to the dry season allowed the identification of the metabolic pathways most affected by the changes in the environment, opening an opportunity for additional characterization of the role of hypoxia in FY symptom intensification. Finally, the initial analysis of a set of 56 proteins/genes differentially expressed in symptomatic plants compared to the asymptomatic ones, independent of the season, has presented pieces of evidence suggesting that breaks in the non-host resistance to non-adapted pathogens and the basal immunity to adapted pathogens, caused by the anaerobic conditions experienced by the plants, might be linked to the onset of this disease. This set of genes might offer the opportunity to develop biomarkers for selecting oil palm plants resistant to this disease and to help pave the way to employing strategies to keep the safety barriers raised and strong.

Transcriptomic responses of oil palm (Elaeis guineensis) stem to waterlogging at plantation in relation to precipitation seasonality

Global warming-induced climate change causes significant agricultural problems by increasing the incidence of drought and flooding events. Waterlogging is an inevitable consequence of these changes but its effects on oil palms have received little attention and are poorly understood. Recent waterlogging studies have focused on oil palm seedlings, with particular emphasis on phenology. However, the transcriptomic waterlogging response of mature oil palms remains elusive in real environments. We therefore investigated transcriptomic changes over time in adult oil palms at plantations over a two-year period with pronounced seasonal variation in precipitation. A significant transcriptional waterlogging response was observed in the oil palm stem core but not in leaf samples when gene expression was correlated with cumulative precipitation over two-day periods. Pathways and processes upregulated or enriched in the stem core response included hypoxia, ethylene signaling, and carbon metabolism. Post-waterlogging recovery in oil palms was found to be associated with responses to heat stress and carotenoid biosynthesis. Nineteen transcription factors (TFs) potentially involved in the waterlogging response of mature oil palms were also identified. These data provide new insights into the transcriptomic responses of planted oil palms to waterlogging and offer valuable guidance on the sensitivity of oil palm plantations to future climate changes.

Transcriptomic analysis of mesocarp tissue during fruit development of the oil palm revealed specific isozymes related to starch metabolism that control oil yield

The oil palm (Elaeis guineensis Jacq.) produces a large amount of oil from the fruit. However, increasing the oil production in this fruit is still challenging. A recent study has shown that starch metabolism is essential for oil synthesis in fruit-producing species. Therefore, the transcriptomic analysis by RNA-seq was performed to observe gene expression alteration related to starch metabolism genes throughout the maturity stages of oil palm fruit with different oil yields. Gene expression profiles were examined with three different oil yields group (low, medium, and high) at six fruit development phases (4, 8, 12, 16, 20, and 22 weeks after pollination). We successfully identified and analyzed differentially expressed genes in oil palm mesocarps during development. The results showed that the transcriptome profile for each developmental phase was unique. Sucrose flux to the mesocarp tissue, rapid starch turnover, and high glycolytic activity have been identified as critical factors for oil production in oil palms. For starch metabolism and the glycolytic pathway, we identified specific gene expressions of enzyme isoforms (isozymes) that correlated with oil production, which may determine the oil content. This study provides valuable information for creating new high-oil-yielding palm varieties via breeding programs or genome editing approaches.

Ganoderma boninense: general characteristics of pathogenicity and methods of control

Ganoderma boninense (G. boninense) is a soil-borne fungus threatening oil palm at the present. It causes basal stem rot disease on oil palm. Within six months, this fungus can cause an oil palm plantation to suffer a significant 43% economic loss. The high persistence and nature of spread of G. boninense in soil make control of the disease challenging. Therefore, controlling the pathogen requires a thorough understanding of the mechanisms that underlie pathogenicity as well as its interactions with host plants. In this paper, we present the general characteristics, the pathogenic mechanisms, and the host's defensive system of G. boninense. We also review upcoming and most promising techniques for disease management that will have the least negative effects on the environment and natural resources.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

A Chromosome-level Reference Genome of African Oil Palm Provides Insights into Its Divergence and Stress Adaptation

The palm family (Arecaceae), consisting of ∼ 2600 species, is the third most economically important family of plants. The African oil palm (Elaeis guineensis) is one of the most important palms. However, the genome sequences of palms that are currently available are still limited and fragmented. Here, we report a high-quality chromosome-level reference genome of an oil palm, Dura, assembled by integrating long reads with ∼ 150× genome coverage. The assembled genome was 1.7 Gb in size, covering 94.5% of the estimated genome, of which 91.6% was assigned into 16 pseudochromosomes and 73.7% was repetitive sequences. Relying on the conserved synteny with oil palm, the existing draft genome sequences of both date palm and coconut were further assembled into chromosomal level. Transposon burst, particularly long terminal repeat retrotransposons, following the last whole-genome duplication, likely explains the genome size variation across palms. Sequence analysis of the VIRESCENS gene in palms suggests that DNA variations in this gene are related to fruit colors. Recent duplications of highly tandemly repeated pathogenesis-related proteins from the same tandem arrays play an important role in defense responses to Ganoderma. Whole-genome resequencing of both ancestral African and introduced oil palms in Southeast Asia reveals that genes under putative selection are notably associated with stress responses, suggesting adaptation to stresses in the new habitat. The genomic resources and insights gained in this study could be exploited for accelerating genetic improvement and understanding the evolution of palms.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

Genomics and Transcriptomics Reveal Genetic Contribution to Population Diversity and Specific Traits in Coconut

Coconut is an economically important palm species with a long history of human use. It has applications in various food, nutraceuticals, and cosmetic products, and there has been renewed interest in coconut in recent years due to its unique nutritional and medicinal properties. Unfortunately, the sustainable growth of the coconut industry has been hampered due to a shortage of good quality seedlings. Genetic improvement through the traditional breeding approach faced considerable obstacles due to its perennial nature, protracted juvenile period, and high heterozygosity. Molecular biotechnological tools, including molecular markers and next-generation sequencing (NGS), could expedite genetic improvement efforts in coconut. Researchers have employed various molecular markers to reveal genetic diversity among coconut populations and for the construction of a genetic map for exploitation in coconut breeding programs worldwide. Whole genome sequencing and transcriptomics on the different varieties have generated a massive amount of publicly accessible sequence data, substantially improving the ability to analyze and understand molecular mechanisms affecting crop performance. The production of high-yielding and disease-resilient coconuts and the deciphering of the complex coconut genome's structure can profit tremendously from these technologies. This paper aims to provide a comprehensive review of the progress of coconut research, using genomics, transcriptomics, and molecular markers initiatives.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Small RNAs and Karma methylation in Elaeis guineensis mother palms are linked to high clonal mantling

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.

Basal Stem Rot of Oil Palm: The Pathogen, Disease Incidence, and Control Methods

Basal stem rot of oil palm caused by Ganoderma boninense is the most serious disease of oil palm in Malaysia, Indonesia, and other oil-palm-producing countries. Economic losses caused by the disease can be up to USD500 million a year. For many years, basal stem rot was found to infect older palm trees of more than 25 to 30 years in age. Only in the 1950s, the disease began to appear in much younger palm trees, 10 to 15 years old, and, in the last decade or so, palm trees as young as 1 year were infected by the disease. The highest incidence occurs in coastal areas of Southeast Asia but the disease has now infected oil palm in inland areas, mainly oil palm planted in peat soils. Disease incidence is also high in areas previously growing coconut or forest. Basal stem rot infection and spread occur through root-to-root contact, and basidiospores that colonize the roots also play a role. In the early stages of infection by G. boninense, the pathogen behaves as a biotroph and later as a necrotroph, secreting cell-wall-degrading enzymes and triggering host defense responses. Genes, gene products, and metabolic pathways involved in oil palm defense mechanisms against G. boninense have been identified and these metabolites have the potential to be used as markers for early detection of the disease. Integrated disease management used to control basal stem rot includes cultural practices, chemical control, and application of biocontrol agents or fertilizers. Early detection tools have also been developed that could assist in management of basal stem rot infections. Development of resistant or tolerant oil palm is still at an early stage; therefore, the existing integrated disease management practices remain the most appropriate methods for managing basal stem rot of oil palm.

Occurrence of triploids in oil palm and their origin

BACKGROUND AND AIMS

Oil palms showing exceptional vigour and dubbed as 'giant palms' were identified in some progeny during breeding. A panel of phenotypical traits were studied to characterize these trees. The hypothesis that gigantism and other anomalies might be linked to polyploidy was investigated.

METHODS

Twenty sib pairs of palms from different crosses, each comprising a giant and a normal oil palm, were studied by flow cytometry with rice 'Nipponbare' as standard reference. In parallel, palms were assessed in the field using 11 phenotypic traits. A principal component analysis (PCA) was conducted to define relationships between these phenotypical traits, and a linear discriminant analysis (LDA) to predict ploidy level and giant classification. Finally, a co-dominant molecular marker study was implemented to highlight the sexual process leading to the formation of 2n gametes.

KEY RESULTS

The first group of oil palms presented an oil palm/rice peak ratio of around 4.8 corresponding to diploid oil palms, whereas the second group presented a ratio of around 7, classifying these plants as triploid. The PCA enabled the classification of the plants in three classes: 21 were normal diploid palms; ten were giant diploid palms; while 11 were giant triploid palms. The LDA revealed three predictors for ploidy classification: phyllotaxy, petiole size and circumference of the plant, but surprisingly not height. The molecular study revealed that triploid palms arose from 2n gametes resulting from the second division restitution of meiosis in parents.

CONCLUSIONS

This study confirms and details the process of sexual polyploidization in oil palm. It also identifies three phenotypical traits to assess the ploidy level of the giant oil palms in the field. In practical terms, our results provide a cheap scientific method to identify polyploid palms in the field.

Integrated consensus genetic map and genomic scaffold re-ordering of oil palm (Elaeis guineensis) genome

A high-quality reference genome is an important resource that can help decipher the genetic basis of traits in combination with linkage or association analyses. The publicly available oil palm draft genome sequence of AVROS pisifera (EG5) accounts for 1.535 Gb of the 1.8 Gb oil palm genome. However, the assemblies are fragmented, and the earlier assembly only had 43% of the sequences placed on pseudo-chromosomes. By integrating a number of SNP and SSR-based genetic maps, a consensus map (AM_EG5.1), comprising of 828.243 Mb genomic scaffolds anchored to 16 pseudo-chromosomes, was generated. This accounted for 54% of the genome assembly, which is a significant improvement to the original assembly. The total length of N50 scaffolds anchored to the pseudo-chromosomes increased by ∼18% compared to the previous assembly. A total of 139 quantitative trait loci for agronomically important quantitative traits, sourced from literature, were successfully mapped on the new pseudo-chromosomes. The improved assembly could also be used as a reference to identify potential errors in placement of specific markers in the linkage groups of the genetic maps used to assemble the consensus map. The 3422 unique markers from five genetic maps, anchored to the pseudo-chromosomes of AM_EG5.1, are an important resource that can be used preferentially to either construct new maps or fill gaps in existing genetic maps. Synteny analysis further revealed that the AM_EG5.1 had high collinearity with the date palm genome cultivar 'Barhee BC4' and shared most of its segmental duplications. This improved chromosomal-level genome is a valuable resource for genetic research in oil palm.

Impact Comparison of El Niño and Ageing Crops on Malaysian Oil Palm Yield

Ageing oil palm crops show a significant correlation with the declining oil palm yield in Malaysia. Not only do aged crops result in lower production, but they are also more costly and difficult to harvest. The Malaysian oil palm yield recovered to the pre-El Niño level after the 1997/98 El Niño event. However, the oil palm yield failed to recover after the recent 2015/16 El Niño. Due to the accumulation of aged oil palm plantations in Malaysia, the financial losses from different magnitudes of El Niño events are increasing. Thirty-four years of monthly oil palm yield trends in Malaysia were compared with the El Niño-free yield dataset to show that the oil palm yield downtrend pattern is the same with or without El Niño events in Malaysia for the most recent 15 years (2005 to 2019). The performance of oil palm yield did not show any significant difference from 2000 to 2019. This study estimates that ageing oil palms would lead to a minimum opportunity loss of USD 431 million by December 2022. Without a proper replanting program, the total combined loss attributable to the ageing crops from 2009 to 2022 is estimated to be USD 3.94 billion, which is more profound than losses due to El Niño events within the same period. This study also concluded that a continuous 7-year replanting scheme of at least 115,000 hectares per year is needed to address the adverse impact of ageing crops on the Malaysian oil palm yield, which accounts for nearly 30% of the global palm oil production.

Multiplex CRISPR/Cas9 gene-editing platform in oil palm targeting mutations in EgFAD2 and EgPAT genes

BACKGROUND

CRISPR/Cas9 is the most powerful and versatile genome-editing tool that permits multiplexed-targeted gene modifications for the genetic enhancement of oil palm. Multiplex genome-editing has recently been developed for modifying multiple loci in a gene or multiple genes in a genome with high precision. This study focuses on the development of high-oleic oil palm, the primary target trait for healthy low-saturated oil. To achieve this, the fatty acid desaturase 2 (FAD2) and palmitoyl-acyl carrier protein thioesterase (PAT) genes, both of which are associated with fatty acid metabolism biosynthesis pathways in oil palm, need to be knocked out. The knockout of FAD2 and PAT leads to an accumulation of oleic acid content in oil palms.

RESULTS

A total of four single-guide RNAs (sgRNAs) were designed in silico based on the genomic sequences of EgFAD2 and EgPAT. Using robust plant CRISPR/Cas9 vector technology, multiple sgRNA expression cassettes were efficiently constructed into a single-binary CRISPR/Cas9 vector to edit the EgFAD2 and EgPAT genes. Each of the constructed transformation vectors was then delivered into oil palm embryogenic calli using the biolistic, Agrobacterium-mediated, and PEG-mediated protoplast transformation methods. Sequence analysis of PCR products from 15 samples confirmed that mutations were introduced at four target sites of the oil palm EgFAD2 and EgPAT genes. Single- and double-knockout mutants of both genes were generated, with large and small deletions within the targeted regions. Mutations found at EgFAD2 and EgPAT target sites indicate that the Cas9/sgRNA genome-editing system effectively knocked out both genes in oil palm.

CONCLUSION

This technology is the first in oil palm to use CRISPR/Cas9 genome-editing to target high-oleic-associated genes. These findings showed that multiplex genome-editing in oil palm could be achieved using multiple sgRNAs. Targeted mutations detected establish that the CRISPR/Cas9 technology offers a great potential for oil palm.

Proteomic changes associated with the development of açaí (Euterpe oleracea Mart.) seeds

Açaí palm (Euterpe oleracea Mart.) seeds are a rich source of mannans, which can be used to generate bioethanol or be converted to high-value D-mannose, in addition to being a source of polyphenols with beneficial health properties. Here, we present a quantitative proteome dataset of açaí seeds at four stages of development (S1, S2, S3, and S4 stages), in which 2465 high confidence proteins were identified and 524 of them show statistically different abundance profiles during development. Several enzymes involved in the biosynthesis of nucleotide-sugars were quantified, especially those dedicated to the formation of GDP-mannose, which showed an increase in abundance between stages S1 and S3. Our data suggest that linear mannans found abundantly in endosperm cell walls are initially deposited as galactomannans, and during development lose the galactosyl groups. Two isoforms of alpha-galactosidase enzymes showed significantly increased abundances in the S3 and S4 stages. Additionally, we quantified the enzymes participating in the central pathway of flavonoid biosynthesis responsible for the formation of catechin and epicatechin, which are subunits of procyanidins, the main class of polyphenols in the açaí seeds. These proteins showed the same pattern of deposition, in which higher abundances were seen in the S1 and S2 stages.

The core triacylglycerol toolbox in woody oil plants reveals targets for oil production bioengineering

Woody oil plants are the most productive oil-bearing species that produce seeds with high levels of valuable triacylglycerols (TAGs). TAGs and their derivatives are the raw materials for many macromolecular bio-based products, such as nylon precursors, and biomass-based diesel. Here, we identified 280 genes encoding seven distinct classes of enzymes (i.e., G3PAT, LPAAT, PAP, DGAT, PDCT, PDAT, and CPT) involved in TAGs-biosynthesis. Several multigene families are expanded by large-scale duplication events, such as G3PATs, and PAPs. RNA-seq was used to survey the expression profiles of these TAG pathway-related genes in different tissues or development, indicating functional redundancy for some duplicated genes originated from the large-scale duplication events, and neo-functionalization or sub-functionalization for some of them. Sixty-two genes showed strong, preferential expression during the period of rapid seed lipid synthesis, suggesting that their might represented the core TAG-toolbox. We also revealed for the first time that there is no PDCT pathway in Vernicia fordii and Xanthoceras sorbifolium. The identification of key genes involved in lipid biosynthesis will be the foundation to plan strategies to develop woody oil plant varieties with enhanced processing properties and high oil content.

EgNRT2.3 and EgNAR2 expression are controlled by nitrogen deprivation and encode proteins that function as a two-component nitrate uptake system in oil palm

Oil palm (Elaeis guineensis Jacq.) is an important crop for oil and biodiesel production. Oil palm plantations require extensive fertilizer additions to achieve a high yield. Fertilizer application decisions and management for oil palm farming rely on leaf tissue and soil nutrient analyses with little information available to describe the key players for nutrient uptake. A molecular understanding of how nutrients, especially nitrogen (N), are taken up in oil palm is very important to improve fertilizer use and formulation practice in oil palm plantations. In this work, two nitrate uptake genes in oil palm, EgNRT2.3 and EgNAR2, were cloned and characterized. Spatial expression analysis showed high expression of these two genes was mainly found in un-lignified young roots. Interestingly, EgNRT2.3 and EgNAR2 were up-regulated by N deprivation, but their expression pattern depended on the form of N source. Promoter analysis of these two genes confirmed the presence of regulatory elements that support these expression patterns. The Xenopus oocyte assay showed that EgNRT2.3 and EgNAR2 had to act together to take up nitrate. The results suggest that EgNRT2.3 and EgNAR2 act as a two-component nitrate uptake system in oil palm.

The genome of Areca catechu provides insights into sex determination of monoecious plants

The areca palm (Areca catechu) has a monoecious spadix, with male flowers on the apical side and females on the basal side. Here, we applied multiomics analysis to investigate sex determination and floral organ development in areca palms. We generated a chromosome-level reference genome of A. catechu with 16 pseudochromosomes, composed of 2.73 Gb and encoding 31 406 genes. Data from RNA-seq and ATAC-seq (assay for transposase accessible chromatin sequencing) suggested that jasmonic acid (JA) synthesis and signal transduction-related genes were differentially expressed between female and male flowers via epigenetic modifications. JA concentration in female flowers was c. 10 times than that in males on the same inflorescence, while JA concentration in hermaphroditic flowers of abnormal inflorescences was about twice that in male flowers of normal inflorescences. JA promotes the development of female flower organs by decreasing the expression of B-function genes, including AGL16, AP3, PIb and PIc. There is also a region on pseudochromosome 15 harboring sex-related genes, including CYP703, LOG, GPAT, AMS and BiP. Among them, CYP703, AMS and BiP were specifically expressed in male flowers.

Genomic selection for morphological and yield-related traits using genome-wide SNPs in oil palm

Oil palm is the most important oil crop worldwide. Colombia is the fourth largest producer, primarily relying on production from interspecific hybrids, derived from crosses between Elaeis oleifera and Elaeis guineensis (OxG). However, conventional breeding can take up to 20 years to generate a new variety. Therefore, reducing the breeding cycle while improving the genetic gain for complex traits is desirable. Genomic selection (GS) is an approach with the potential to achieve this goal. In this study, we evaluated 431 F1 interspecific hybrids (OxG) and 444 backcrosses (BC1) for morphological and yield-related traits. Genomic predictions were performed with the G-BLUP model using three different population datasets for training the model: the same population (TRN1), the other population (TRN2), and both populations (TRN1+2). Higher multi-family prediction accuracies were obtained for foliar area (0.3 in OxG) and trunk height (0.47 in BC1) when the model was trained with TRN1. Single-family prediction accuracies were lower in the OxG compared to BC1 families for traits such as trunk diameter, trunk height, bunch number, and yield using TRN1. Conversely, lower prediction accuracies were obtained for most traits when the model was trained using TRN2 (< 0.1). Multi-trait models showed a substantial increase of the predictions for traits such as yield (0.22 for OxG and 0.44 for BC1), because of the genetic correlations between traits. The results herein highlighted the potential of GS for parental selection in OxG and BC1 populations, but further studies are required to improve the models to select individuals by their genetic value.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01341-5.

A chromosome-level genome assembly for Dracaena cochinchinensis reveals the molecular basis of its longevity and formation of dragon's blood

Dracaena, a remarkably long-lived and slowly maturing species of plant, is world famous for its ability to produce dragon's blood, a precious traditional medicine used by different cultures since ancient times. However, there is no detailed and high-quality genome available for this species at present; thus, the molecular mechanisms that underlie its important traits are largely unknown. These factors seriously limit the protection and regeneration of this rare and endangered plant resource. Here, we sequenced and assembled the genome of Dracaena cochinchinensis at the chromosome level. The D. cochinchinensis genome covers 1.21 Gb with a scaffold N50 of 50.06 Mb and encodes 31 619 predicted protein-coding genes. Analysis showed that D. cochinchinensis has undergone two whole-genome duplications and two bursts of long terminal repeat insertions. The expansion of two gene classes, cis-zeatin O-glucosyltransferase and small auxin upregulated RNA, were found to account for its longevity and slow growth. Two transcription factors (bHLH and MYB) were found to be core regulators of the flavonoid biosynthesis pathway, and reactive oxygen species were identified as the specific signaling molecules responsible for the injury-induced formation of dragon's blood. Our study provides high-quality genomic information relating to D. cochinchinensis and significant insight into the molecular mechanisms responsible for its longevity and formation of dragon's blood. These findings will facilitate resource protection and sustainable utilization of Dracaena.

Structural and functional analysis of stress-inducible genes and their promoters selected from young oil palm (Elaeis guineensis) under salt stress

Background

Soil salinity is a problem in more than 100 countries across all continents. It is one of the abiotic stress that threatens agriculture the most, negatively affecting crops and reducing productivity. Transcriptomics is a technology applied to characterize the transcriptome in a cell, tissue, or organism at a given time via RNA-Seq, also known as full-transcriptome shotgun sequencing. This technology allows the identification of most genes expressed at a particular stage, and different isoforms are separated and transcript expression levels measured. Once determined by this technology, the expression profile of a gene must undergo validation by another, such as quantitative real-time PCR (qRT-PCR). This study aimed to select, annotate, and validate stress-inducible genes—and their promoters—differentially expressed in the leaves of oil palm (Elaeis guineensis) plants under saline stress.

Results

The transcriptome analysis led to the selection of 14 genes that underwent structural and functional annotation, besides having their expression validated using the qRT-PCR technique. When compared, the RNA-Seq and qRT-PCR profiles of those genes resulted in some inconsistencies. The structural and functional annotation analysis of proteins coded by the selected genes showed that some of them are orthologs of genes reported as conferring resistance to salinity in other species. There were those coding for proteins related to the transport of salt into and out of cells, transcriptional regulatory activity, and opening and closing of stomata. The annotation analysis performed on the promoter sequence revealed 22 distinct types of cis-acting elements, and 14 of them are known to be involved in abiotic stress.

Conclusion

This study has helped validate the process of an accurate selection of genes responsive to salt stress with a specific and predefined expression profile and their promoter sequence. Its results also can be used in molecular-genetics-assisted breeding programs. In addition, using the identified genes is a window of opportunity for strategies trying to relieve the damages arising from the salt stress in many glycophyte crops with economic importance.

Insights from a Multi-Omics Integration (MOI) Study in Oil Palm (Elaeis guineensis Jacq.) Response to Abiotic Stresses: Part Two-Drought

Drought and salinity are two of the most severe abiotic stresses affecting agriculture worldwide and bear some similarities regarding the responses of plants to them. The first is also known as osmotic stress and shows similarities mainly with the osmotic effect, the first phase of salinity stress. Multi-Omics Integration (MOI) offers a new opportunity for the non-trivial challenge of unraveling the mechanisms behind multigenic traits, such as drought and salinity resistance. The current study carried out a comprehensive, large-scale, single-omics analysis (SOA) and MOI studies on the leaves of young oil palm plants submitted to water deprivation. After performing SOA, 1955 DE enzymes from transcriptomics analysis, 131 DE enzymes from proteomics analysis, and 269 DE metabolites underwent MOI analysis, revealing several pathways affected by this stress, with at least one DE molecule in all three omics platforms used. Moreover, the similarities and dissimilarities in the molecular response of those plants to those two abiotic stresses underwent mapping. Cysteine and methionine metabolism (map00270) was the most affected pathway in all scenarios evaluated. The correlation analysis revealed that 91.55% of those enzymes expressed under both stresses had similar qualitative profiles, corroborating the already known fact that plant responses to drought and salinity show several similarities. At last, the results shed light on some candidate genes for engineering crop species resilient to both abiotic stresses.

The heterologous expression of conserved Glycine max (soybean) mitogen activated protein kinase 3 (MAPK3) paralogs suppresses Meloidogyne incognita parasitism in Gossypium hirsutum (upland cotton)

Two conserved Glycine max (soybean) mitogen activated protein kinase 3 (MAPK3) paralogs function in defense to the parasitic soybean cyst nematode Heterodera glycines. Gene Ontology analyses of RNA seq data obtained from MAPK3-1-overexpressing (OE) and MAPK3-2-OE roots compared to their control, as well as MAPK3-1-RNA interference (RNAi) and MAPK3-2-RNAi compared to their control, hierarchically orders the induced and suppressed genes, strengthening the hypothesis that their heterologous expression in Gossypium hirsutum (upland cotton) would impair parasitism by the root knot nematode (RKN) Meloidogyne incognita. MAPK3-1 expression (E) in G. hirsutum suppresses the production of M. incognita root galls, egg masses, and second stage juveniles (J2s) by 80.32%, 82.37%, and 88.21%, respectfully. Unexpectedly, egg number increases by 28.99% but J2s are inviable. MAPK3-2-E effects are identical, statistically. MAPK3-1-E and MAPK3-2-E decreases root mass 1.49-fold and 1.55-fold, respectively, as compared to the pRAP15-ccdB-E control. The reproductive factor (RF) of M. incognita for G. hirsutum roots expressing MAPK3-1-E or MAPK3-2-E decreases 60.39% and 50.46%, respectively, compared to controls. The results are consistent with upstream pathogen activated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) functioning in defense to H. glycines. The experiments showcase the feasibility of employing MAPK3, through heterologous expression, to combat M. incognita parasitism, possibly overcoming impediments otherwise making G. hirsutum's defense platform deficient. MAPK homologs are identified in other important crop species for future functional analyses.

The early response of oil palm (Elaeis guineensis Jacq.) plants to water deprivation: Expression analysis of miRNAs and their putative target genes, and similarities with the response to salinity stress

Oil palm (Elaeis guineensis Jacq.) is a oilseed crop of great economic importance drastically affected by abiotic stresses. MicroRNAs (miRNAs) play crucial roles in transcription and post-transcription regulation of gene expression, being essential molecules in the response of plants to abiotic stress. To better understand the molecular mechanisms behind the response of young oil palm plants to drought stress, this study reports on the prediction and characterization of miRNAs and their putative target genes in the apical leaf of plants subjected to 14 days of water deprivation. Then, the data from this study were compared to the data from a similar study that focused on salinity stress. Both, the drought-and salt-responsive miRNAs and their putative target genes underwent correlation analysis to identify similarities and dissimilarities among them. Among the 81 identified miRNAs, 29 are specific for oil palm, including two (egu-miR28ds and egu-miR29ds) new ones - described for the first time. As for the expression profile, 62 miRNAs were significantly differentially expressed under drought stress, being five up-regulated (miR396e, miR159b, miR529b, egu-miR19sds, and egu-miR29ds) and 57 down-regulated. Transcription factors, such as MYBs, HOXs, and NF-Ys, were predicted as putative miRNA-target genes in oil palm under water deprivation; making them the most predominant group of such genes. Finally, the correlation analysis study revealed a group of putative target genes with similar behavior under salt and drought stresses. Those genes that are upregulated by these two abiotic stresses encode lncRNAs and proteins linked to stress tolerance, stress memory, modulation of ROS signaling, and defense response regulation to abiotic and biotic stresses. In summary, this study provides molecular evidence for the possible involvement of miRNAs in the drought stress response in oil palm. Besides, it shows that, at the molecular level, there are many similarities in the response of young oil palm plants to these two abiotic stresses.

Characterization and functional analysis of the MADS-box EgAGL9 transcription factor from the mesocarp of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is one of the most important oil crops in the world, and compared to all oil crops, it has the highest productive efficiency. In the present study, a MADS-box transcription factor of the AGAMOUS class, named EgAGL9, was identified by expression profile analysis in the different developmental stages of oil palm mesocarp. Real-time quantitative PCR results confirmed that the expression of EgAGL9 increased rapidly during the last stages of oil palm mesocarp development. Then, three downstream genes, including EgSAD (Stearoyl-ACP desaturase), EgTSA (Tryptophan synthase) and EgSDH (Succinate dehydrogenase), were screened by ChIP-Seq and data analysis. EMSA analysis verified that EgAGL9 interacted with the promoter regions of EgSAD, EgTSA and EgSDH. Moreover, the expression levels of EgSAD, EgTSA and EgSDH were downregulated in EgAGL9-overexpressing protoplasts and calli of oil palm. Compared to WT, the total lipid content and ratio of unsaturated fatty acids in transgenic calli (including oleic acid, linoleic acid and linolenic acid) were significantly decreased. Together, these results revealed that these three EgAGL9-regulated genes are involved in regulatory pathways in the oil palm mesocarp. Compared with previous studies, the present study provides a new research strategy for understanding of the molecular regulatory pathways of lipid metabolism in mesocarp of oil palm. The obtained results will bring a new perspective for a comprehensive understanding of the regulation of the metabolic accumulation in the oil palm mesocarp.

Oil Palm Fruits Ripeness Classification Based on the Characteristics of Protein, Lipid, Carotene, and Guanine/Cytosine from the Raman Spectra

The capacity of palm oil production is directly affected by the ripeness of the fresh fruit bunches (FFB) upon harvesting. Conventional harvesting standards rely on rigid harvesting scheduling as well as the number of fruitlets that have loosened from the bunch. Harvesting is usually done every 10 to 14 days, and an FFB is deemed ready to be harvested if there are around 5 to 10 empty sockets on the fruit bunch. Technology aided by imaging techniques relies heavily on the color of the fruit bunch, which is highly dependent on the surrounding light intensities. In this study, Raman spectroscopy is used for ripeness classification of oil palm fruits, based on the molecular assignments extracted from the Raman bands between 1240 cm^-1 and 1360 cm^-1. The Raman spectra of 52 oil palm fruit samples which contain the fingerprints of different organic compounds were collected. Signal processing was applied to perform baseline correction and to reduce background noises. Characteristic data of the organic compounds were extracted through deconvolution and curve fitting processes. Subsequently, a correlation study between organic compounds was developed and eight hidden Raman peaks including protein, beta carotene, carotene, lipid, guanine/cytosine, chlorophyll-a, and tryptophan were successfully located. Through ANOVA statistical analysis, a total of six peak intensities from proteins through Amide III (β-sheet), beta-carotene, carotene, lipid, guanine/cytosine, and carotene and one peak location from lipid were found to be significant. An automated oil palm fruit ripeness classification system deployed with artificial neural network (ANN) using the seven signification features showed an overall performance of 97.9% accuracy. An efficient and accurate ripeness classification model which uses seven significant Raman peak features from the correlation analysis between organic compounds was successfully developed.

Genome-wide analysis of the CAD gene family reveals two bona fide CAD genes in oil palm

Cinnamyl alcohol dehydrogenase (CAD) is the key enzyme for lignin biosynthesis in plants. In this study, genome-wide analysis was performed to identify CAD genes in oil palm (Elaeis guineensis). Phylogenetic analysis was then conducted to select the bona fide EgCADs. The bona fide EgCAD genes and their respective 5' flanking regions were cloned and analysed. Their expression profiles were evaluated in various organs using RT-PCR. Seven EgCAD genes (EgCAD1-7) were identified and divided into four phylogenetic groups. EgCAD1 and EgCAD2 display high sequence similarities with other bona fide CADs and possess all the signature motifs of the bona fide CAD. They also display similar 3D protein structures. Gene expression analysis showed that EgCAD1 was expressed most abundantly in the root tissues, while EgCAD2 was expressed constitutively in all the tissues studied. EgCAD1 possesses only one transcription start site, while EgCAD2 has five. Interestingly, a TC microsatellite was found in the 5' flanking region of EgCAD2. The 5' flanking regions of EgCAD1 and EgCAD2 contain lignin-associated regulatory elements i.e. AC-elements, and other defence-related motifs, including W-box, GT-1 motif and CGTCA-motif. Altogether, these results imply that EgCAD1 and EgCAD2 are bona fide CAD involved in lignin biosynthesis during the normal development of oil palm and in response to stresses. Our findings shed some light on the roles of the bona fide CAD genes in oil palm and pave the way for manipulating lignin content in oil palm through a genetic approach.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03208-0.

The slow-evolving Acorus tatarinowii genome sheds light on ancestral monocot evolution

Monocots are one of the most diverse groups of flowering plants, and tracing the evolution of their ancestral genome into modern species is essential for understanding their evolutionary success. Here, we report a high-quality assembly of the Acorus tatarinowii genome, a species that diverged early from all the other monocots. Genome-wide comparisons with a range of representative monocots characterized Acorus as a slowly evolved genome with one whole-genome duplication. Our inference of the ancestral monocot karyotypes provides new insights into the chromosomal evolutionary history assigned to modern species and reveals the probable molecular functions and processes related to the early adaptation of monocots to wetland or aquatic habitats (that is, low levels of inorganic phosphate, parallel leaf venation and ephemeral primary roots). The evolution of ancestral gene order in monocots is constrained by gene structural and functional features. The newly obtained Acorus genome offers crucial evidence for delineating the origin and diversification of monocots, including grasses.

Genome Assembly and Analysis of the Flavonoid and Phenylpropanoid Biosynthetic Pathways in Fingerroot Ginger (Boesenbergia rotunda)

Boesenbergia rotunda (Zingiberaceae), is a high-value culinary and ethno-medicinal plant of Southeast Asia. The rhizomes of this herb have a high flavanone and chalcone content. Here we report the genome analysis of B. rotunda together with a complete genome sequence as a hybrid assembly. B. rotunda has an estimated genome size of 2.4 Gb which is assembled as 27,491 contigs with an N50 size of 12.386 Mb. The highly heterozygous genome encodes 71,072 protein-coding genes and has a 72% repeat content, with class I TEs occupying ~67% of the assembled genome. Fluorescence in situ hybridization of the 18 chromosome pairs at the metaphase showed six sites of 45S rDNA and two sites of 5S rDNA. An SSR analysis identified 238,441 gSSRs and 4604 EST-SSRs with 49 SSR markers common among related species. Genome-wide methylation percentages ranged from 73% CpG, 36% CHG and 34% CHH in the leaf to 53% CpG, 18% CHG and 25% CHH in the embryogenic callus. Panduratin A biosynthetic unigenes were most highly expressed in the watery callus. B rotunda has a relatively large genome with a high heterozygosity and TE content. This assembly and data (PRJNA71294) comprise a source for further research on the functional genomics of B. rotunda, the evolution of the ginger plant family and the potential genetic selection or improvement of gingers.

Insights from a Multi-Omics Integration (MOI) Study in Oil Palm (Elaeis guineensis Jacq.) Response to Abiotic Stresses: Part One—Salinity

Oil palm (Elaeis guineensis Jacq.) is the number one source of consumed vegetable oil nowadays. It is cultivated in areas of tropical rainforest, where it meets its natural condition of high rainfall throughout the year. The palm oil industry faces criticism due to a series of practices that was considered not environmentally sustainable, and it finds itself under pressure to adopt new and innovative procedures to reverse this negative public perception. Cultivating this oilseed crop outside the rainforest zone is only possible using artificial irrigation. Close to 30% of the world’s irrigated agricultural lands also face problems due to salinity stress. Consequently, the research community must consider drought and salinity together when studying to empower breeding programs in order to develop superior genotypes adapted to those potential new areas for oil palm cultivation. Multi-Omics Integration (MOI) offers a new window of opportunity for the non-trivial challenge of unraveling the mechanisms behind multigenic traits, such as drought and salinity tolerance. The current study carried out a comprehensive, large-scale, single-omics analysis (SOA), and MOI study on the leaves of young oil palm plants submitted to very high salinity stress. Taken together, a total of 1239 proteins were positively regulated, and 1660 were negatively regulated in transcriptomics and proteomics analyses. Meanwhile, the metabolomics analysis revealed 37 metabolites that were upregulated and 92 that were downregulated. After performing SOA, 436 differentially expressed (DE) full-length transcripts, 74 DE proteins, and 19 DE metabolites underwent MOI analysis, revealing several pathways affected by this stress, with at least one DE molecule in all three omics platforms used. The Cysteine and methionine metabolism (map00270) and Glycolysis/Gluconeogenesis (map00010) pathways were the most affected ones, each one with 20 DE molecules.