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

Somaclonal variation arises in plants and animals when differentiated somatic cells are induced into a pluripotent state, but the resulting clones differ from each other and from their parents. In agriculture, somaclonal variation has hindered micropropagation of elite hybrids and genetically modified crops, but the mechanism remains a mystery1. The oil palm fruit abnormality, mantled, is a somaclonal variant arising from tissue culture that drastically reduces yield, and has largely halted efforts to clone elite hybrids for oil production2–4. Widely regarded as epigenetic5, mantling has defied explanation, but here we identify the MANTLED gene using Epigenome Wide Association Studies. DNA hypomethylation of a LINE retrotransposon related to rice Karma, in the intron of the homeotic gene DEFICIENS, is common to all mantled clones and is associated with alternative splicing and premature termination. Dense methylation near the Karma splice site (the Good Karma epiallele) predicts normal fruit set, while hypomethylation (the Bad Karma epiallele) predicts homeotic transformation, parthenocarpy and dramatic loss of yield. Loss of Karma methylation and small RNA in tissue culture contributes to the origin of mantled, while restoration in spontaneous revertants accounts for non-Mendelian inheritance. The ability to predict and cull mantling at the plantlet stage will facilitate the introduction of higher performing clones and optimize environmentally sensitive land resources.

A microsatellite-based high-density linkage map for oil palm (Elaeis guinensis Jacq.) was constructed from a cross between two heterozygous parents, a tenera palm from the La Me population (LM2T) and a dura palm from the Deli population (DA10D). A set of 390 simple sequence repeat (SSR) markers was developed in oil palm from microsatellite-enriched libraries and evaluated for polymorphism along with 21 coconut SSRs. A dense and genome-wide microsatellite framework as well as saturating amplified fragments length polymorphisms (AFLPs) allowed the construction of a linkage map consisting of 255 microsatellites, 688 AFLPs and the locus of the Sh gene, which controls the presence or absence of a shell in the oil palm fruit. An AFLP marker E-Agg/M-CAA132 was mapped at 4.7 cM from the Sh locus. The 944 genetic markers were distributed on 16 linkage groups (LGs) and covered 1,743 cM. Our linkage map is the first in oil palm to have 16 independent linkage groups corresponding to the plant's 16 homologous chromosome pairs. It is also the only high-density linkage map with as many microsatellite markers in an Arecaceae species and represents an important step towards quantitative trait loci analysis and physical mapping in the E. guineensis species.

A key event in the domestication and breeding of the oil palm, Elaeis guineensis, was loss of the thick coconut-like shell surrounding the kernel. Modern E. guineensis has three fruit forms, dura (thick-shelled), pisifera (shell-less) and tenera (thin-shelled), a hybrid between dura and pisifera1–4. The pisifera palm is usually female-sterile but the tenera yields far more oil than dura, and is the basis for commercial palm oil production in all of Southeast Asia5. Here, we describe the mapping and identification of the Shell gene responsible for the different fruit forms. Using homozygosity mapping by sequencing we found two independent mutations in the DNA binding domain of a homologue of the MADS-box gene SEEDSTICK (STK) which controls ovule identity and seed development in Arabidopsis. The Shell gene is responsible for the tenera phenotype in both cultivated and wild palms from sub-Saharan Africa, and our findings provide a genetic explanation for the single gene heterosis attributed to Shell, via heterodimerization. This gene mutation explains the single most important economic trait in oil palm, and has implications for the competing interests of global edible oil production, biofuels and rainforest conservation6.

A quantitative trait locus (QTL) analysis designed for a multi-parent population was carried out and tested in oil palm (Elaeis guineensis Jacq.), which is a diploid cross-fertilising perennial species. A new extension of the MCQTL package was especially designed for crosses between heterozygous parents. The algorithm, which is now available for any allogamous species, was used to perform and compare two types of QTL search for small size families, within-family analysis and across-family analysis, using data from a 2 × 2 complete factorial mating experiment involving four parents from three selected gene pools. A consensus genetic map of the factorial design was produced using 251 microsatellite loci, the locus of the Sh major gene controlling fruit shell presence, and an AFLP marker of that gene. A set of 76 QTLs involved in 24 quantitative phenotypic traits was identified. A comparison of the QTL detection results showed that the across-family analysis proved to be efficient due to the interconnected families, but the family size issue is just partially solved. The identification of QTL markers for small progeny numbers and for marker-assisted selection strategies is discussed.Electronic supplementary materialThe online version of this article (doi:10.1007/s00122-010-1284-y) contains supplementary material, which is available to authorized users.