Improving palm oil quality through identification and mapping of the lipase gene causing oil deterioration

Abstract

The oil palm fruit mesocarp contains high lipase activity that increases free fatty acids and necessitates post-harvest inactivation by heat treatment of fruit bunches. Even before heat treatment the mesocarp lipase activity causes consequential oil losses and requires costly measures to limit free fatty acids quantities. Here we demonstrate that elite low-lipase lines yield oil with substantially less free fatty acids than standard genotypes, allowing more flexibility for post-harvest fruit processing and extended ripening for increased yields. We identify the lipase and its gene cosegregates with the low-/high-lipase trait, providing breeders a marker to rapidly identify potent elite genitors and introgress the trait into major cultivars. Overall, economic gains brought by wide adoption of this material could represent up to one billion dollars per year. Expected benefits concern all planters but are likely to be highest for African smallholders who would be more able to produce oil that meets international quality standards.

ARTICLE

Received 12 Dec 2012 | Accepted 17 Jun 2013 | Published 16 Jul 2013

DOI: 10.1038/ncomms3160 OPEN

Improving palm oil quality through identication and mapping of the lipase gene causing oil deterioration

F. Morcillo1,*, D. Cros2,*, N. Billotte2, G.-F. Ngando-Ebongue3, H. Domonhdo4, M. Pizot5, T. Cullar2, S. Espout2,R. Dhouib6, F. Bourgis7,8, S. Claverol9, T.J. Tranbarger5, B. Nouy10 & V. Arondel7,8

The oil palm fruit mesocarp contains high lipase activity that increases free fatty acids and necessitates post-harvest inactivation by heat treatment of fruit bunches. Even before heat treatment the mesocarp lipase activity causes consequential oil losses and requires costly measures to limit free fatty acids quantities. Here we demonstrate that elite low-lipase lines yield oil with substantially less free fatty acids than standard genotypes, allowing more exibility for post-harvest fruit processing and extended ripening for increased yields. We identify the lipase and its gene cosegregates with the low-/high-lipase trait, providing breeders a marker to rapidly identify potent elite genitors and introgress the trait into major cultivars. Overall, economic gains brought by wide adoption of this material could represent up to one billion dollars per year. Expected benets concern all planters but are likely to be highest for African smallholders who would be more able to produce oil that meets international quality standards.

1 CIRAD, UMR DIADE, F-34398 Montpellier, France. 2 CIRAD, UMR AGAP, F-34398 Montpellier, France. 3 IRAD, Centre de Recherches sur le Palmier Huile de la Dibamba (CEREPAH), BP243 Douala, Cameroon. 4 INRAB, CRA-PP, BP1, Pob, Benin. 5 IRD/UM2, UMR DIADE, F-34394 Montpellier, France. 6 CNRS, EIPL, UPR9025, 31 Chemin Joseph-Aiguier, F-13008 Marseille, France. 7 CNRS, LBM, UMR 5200, BP 71, F-33883 Villenave dOrnon, France. 8 University of Bordeaux, LBM, UMR 5200, BP 71, F-33883 Villenave dOrnon, France. 9 University of Bordeaux, Centre de Gnomique Fonctionnelle, 146 rue Lo Saignat, F-33000 Bordeaux, France. 10 PalmElit SAS, Parc Agropolis, 2214 Boulevard de la Lironde, F-34980 Montferrier-sur-Lez, France. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to V.A. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160

Oil extracted from the oil palm (Elaeis guineensis)1 fruit mesocarp is the largest source of vegetable oil, with a world market share of B32% (Fig. 1a,b). The share of

palm oil production by small, family-owned, estates is 30% worldwide and reaches 80% in Nigeria, Africas largest producer2. It is expected that the oil palm share of the market will need to increase to meet the doubled demand for vegetable oil predicted by 2030 (refs 36). One major constraint that weighs heavily on yield and harvest management is the rapid oil acidication due to the high lipase (HL) activity that releases free fatty acids (FFA) in the mesocarp of bruised ripe fruits710. Acidity has a strong impact on the quality of commodity oil because an FFA content of 45% is considered to be unt for human consumption11.

Lipase activity can hydrolyze 30% of the oil content within 5 min10 and large-scale exploitation of oil palm became feasible only after Fickendey12 found a century ago that the lipase could be inactivated by heat treatment of harvested bunches. However, the mesocarp lipase can be active prior to heat treatment and release FFA. While FFA level of crude palm oil from large estates is usually 34%13, oil from smallholders shows higher values especially in Africa where FFA level reaches 9% in Nigeria14 and can exceed 15%15. FFAs have to be removed through an energy-costly rening process that leads to oil losses13, so that most of the palm oil produced by African smallholders is not competitive on international markets14. To limit acidity, the bunches must be harvested before they are overripe, heat-sterilized rapidly (ideally within 1 day) to inactivate lipase activity and processed immediately. However, bunches also need to be harvested when they are as ripe as possible to maximize oil yield. These conicting requirements narrow the harvest window down to B1 week and oblige planters to harvest every 710 days, as oil palm continuously produces fruit bunches yearlong (Fig. 1a). As a result, more than a third of the plantation-operating costs are dedicated to harvest1. Palm oil is currently the least expensive oil on the market, due to high productivity and low labour costs in the producing countries, mainly Indonesia and Malaysia4. This will most likely change in the future as oil palm harvest is largely non-mechanized and labour-intensive, and wages are steadily increasing as Southeast Asia develops1,4. It is therefore important that management improves to keep harvest-associated costs as low as possible.

While oil acidication concerns all planters, it is more acute for smallholders who may not have access to efcient milling

facilities close to their plantation and whose overripe bunches may be rejected by mills. Indeed, the problem extends from the whole industry down to the consumer in those African countries where crude, unrened, palm oil is a major commodity oil. The economic importance of palm oil acidity made the mesocarp lipase the most widely studied enzyme of oil palm16, as early as in 1910 (ref. 12). However, the source of this lipase activity has not yet been identied.

We previously identied trees with a low mesocarp lipase activity11. Here we show that the oil extracted from these trees, using different procedures that mimic standard processes, have signicantly lower acidity than that from the standard, HL genotypes currently cultivated. We also demonstrate that the low-/high-acidity trait is monogenic and identify the lipase responsible for oil acidity. Finally, we show that the lipase gene cosegregates with the low (LL)/HL trait, making it a useful DNA marker for marker-assisted selection of the trait, a key technique for fast breeding of crop trees17.

ResultsOil FFA content is lower in LL lines. To assess the performance of LL trees, we used either elite tenera material (Fig. 1) from a Deli (D1) La M (L1) cross already used for commercial oil

production or progenies of its dura Deli (D1) parent. The oil extraction ratio was 22.15% to 21.73% in HL and LL lines from D1, respectively, thus showing no inuence of the LL trait on oil yield (Table 1).

When ripe bunches were heat-sterilized a few hours after harvest, the oil extracted from standard fruits contained0.660.36% (s.d.) FFA, while the oil from LL fruits showed less than half FFA, down to 0.300.02% (Fig. 2a). These ideal conditions are rarely met in practice and the average acidity of crude palm oil from large estates is usually B34%13 and can exceed 10% in small plantations. Therefore, we further focused on two critical steps of fruit processing that increase oil acidity: delayed harvesting of bunches and post-harvest delayed processing. When bunches were left on trees up to 30 days after ripeness, the oil from standard trees contained 412% FFA, whereas that from LL trees showed only half that value (Fig. 2b). Therefore, it is clear from our data that the use of LL genotypes should allow signicant extended ripening without a negative

Oil palm mesocarp

Soybean Rapeseed Sunflower Groundnut

Oil palm kernel

50

45

40

35

30

25

20

10

15

5

0

Oil production (Mt)

1995

1990 2000 2005 2010Year

Figure 1 | Oil palm is a major oil-crop. (a) Palm tree with several bunches at different stages of ripeness; a palm tree continuously produces bunches yearlong (828 bunches per year, depending on age). While most cultivated trees are derived from Deli germplasm, several other germplasms such as La Me are signicantly used by breeders. (b) Open tenera fruits showing eshy mesocarp from which palm oil is extracted and kernel that contains an oil with shorter fatty acids. The amount of mesocarp per fruit is heavily inuenced by the Sh gene that governs shell thickness. Tenera trees, which are the most productive and most widely cultivated, are heterozygous Sh / . Dura trees are Sh / and contain a thicker shell and a thinner mesocarp,

thus producing less oil. (c) World vegetable oil production by crop, showing the increasing share of palm oil (FAOSTAT, http://faostat.org, Dec 2011).

2 NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160 ARTICLE

Table 1 | Yield is not affected in low-lipase lines.

Line Fruit to bunch Wet mesocarp to fruit Oil to dry mesocarp Mesocarp water content Oil to bunch High Lipase 59.270.58 68.180.77 79.260.28 31.620.31* 21.730.28

Low Lipase 59.410.89 70.121.18 79.630.43 33.180.47* 22.150.43

Ratios of main components of yield are in %s.e. Means are from 301 bunches from 21 high-lipase trees and 126 bunches from 9 low-lipase trees. A one-way analysis of variance (F test, factor group) was carried out after having tested the model s validity (normality, outliers).*indicates a statistically (P 0.009) signicant difference between low- and high-lipase lines. Only mesocarp water content appears to differ slightly between the two lines.

Oil acidity (%FFA)

1

0.75

0.5

0.25

0

15

10

5

0

Standard

Low lipase

Oil acidity (%FFA)

Standard

Low lipase

Genotypes

Control

Standard Low lipase

10 20 30

Days after normal harvest

Treated

Standard Low lipase

Oil acidity (%FFA)

7.5

5.0

2.5

0

Oil acidity (%FFA)

7.5

5.0

2.5

0

6 9 12

6 9 12

Days before processing

Days before processing

Figure 2 | LL genotype yields oil of lesser acidity than standard HL genotypes. (a) Ripe bunches (15 fallen fruits) were heat-sterilized a few hours after harvest. Oil was extracted immediately and its acidity titrated, expressed as %FFA). n 4 for LL fruits and 8 for HL (b) Overripe bunches were left

up to 30 days on the tree after ripeness (15 fallen fruits) before harvest. Bunches and fallen fruits were then sterilized and processed. (c,d) Fruits were sprayed with water (control, c) or with fungicide (d), sun-dried and stored up to 12 days before heat sterilization and oil extraction; n 35. Analyses

of variance were carried out according to the general linear model and differences between standard and LL materials were signicant at Po0.001 (ac) and Po0.012 (d) (Tukeys test). Bunches were harvested from the segregating D1 L1 cross (a,c and d) or from selfed D1 progenies (LL P4 and

standard P5, see plant material in Material and Methods section) (b). Error bars represent s.d.

impact on oil quality. This should increase oil yield as oil content of the bunches increases by B7% after rst fruit drop starts18,19.

Another critical step is the delay between harvest and bunch processing, especially in small estates distant from mills. When ripe bunches were processed 3 days after harvest, FFA content of oil from the standard genotype was 1.500.51% compared with only 0.350.09% from LL trees. The latter value is almost identical to that observed when bunches are processed immediately after harvest. When processing was further delayed, FFA content was still o5% at 12 days in LL lines (Fig. 2c). This result clearly indicates that LL lines will allow larger delays for processing bunches. This should strongly benet smallholders distant from mills in areas lacking good transportation infrastructures2, as is frequently the case in Africa and for those

who traditionally leave the bunches to ferment for several days before processing to allow fruit loosening20.

There is a controversy as to the respective contributions of endogenous and fungal lipases to oil acidity21. Treating bunches with fungicide upon harvest (Fig. 2d) led to reduced oil acidity, with the oil from LL fruits containing 3% FFA after 12 days. Again, the oil from LL trees contained 30% less FFA compared to that of standard trees. Our data show that, although both mesocarp and fungal lipases are contributors to oil acidity, the endogenous lipase is clearly the dominant contributor to the immediate post-harvest rise in FFA. To conclude, in conditions known to favour high levels of FFA, LL genotypes always showed a delayed development of acidity when compared to the standard genotype. Therefore, clearly the LL genotypes can produce an oil

NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160

of signicantly lower acidity and therefore higher quality than that from presently used genotypes.

Mesocarp lipase abundance correlates with lipase activity. We devised a functional proteomic approach to identify lipases from a crude protein extract. The approach is based on the use of an inhibitor, tetrahydrolipstatin22, that can bind covalently to the lipase catalytic serine residue. Thus, a crude protein extract from ripe mesocarp from the HL genotype was incubated with radiolabelled tetrahydrolipstatin. A single protein band of 55 kDa was subsequently detected by autoradiography (Fig. 3a). De novo sequencing (Supplementary Fig. S1A) showed signicant

homology to a protein with high similarity to castor bean acid lipase23. The MS/MS spectra allowed us to identify from palm sequences2426 a cDNA named EgLIP1 (GenBank accession number JX556215), coding for a protein with strong homology to castor bean acid lipase (34.6% identity, 51.8% similarity, Supplementary Fig. S1B). Antibodies raised against a specic peptide detected a 55-kDa protein in HL mesocarp, whereas no protein could be detected in the mesocarp from the LL line (Fig. 3b). Immunolocalization data (Fig. 3c) showed that the putative lipase is localized adjacent to the outer surface of the lipid droplets. Taken together, these data suggest that EgLIP1 encodes the mesocarp lipase.

a

SDS-PAGE

14C Label

b

116 kDa

66 kDa

45 kDa

35 kDa

25 kDa

18.4 kDa

HL LL

85 kDa

47 kDa

EgLip1

c

d

Lipase activity (U g1 DW)

800 90

75

60

45

30

15

0

600

400

200

0 30 70 100 120 140 160 DAP

Lipids (% DW)

Ld

Ld

Ld

85 kDa

47 kDa

70 100 120 140 160 DAP

EgLip1

Ld

Relative abundance

of EgLip1 transcript

400

300

200

100

Ld Ld

0 0 30 100 120 140 160 gs1 gs2

DAP Germinated seeds

Figure 3 | Isolation and characterization of the oil palm mesocarp lipase EgLIP1. (a) Mesocarp crude protein extract was incubated with radioactive tetrahydrolipstatin and loaded on a 10% SDS-PAGE gel (left panel). The gel was dried and autoradiographed (right panel). (b) Western analysisof the EgLIP1 protein in the fruit mesocarp from LL and HL lines (D1 progenies). (c) Immunolocalization of EgLIP1. Confocal analysis of mesocarp sections stained with Nile red (upper panel) to detect lipids in lipid droplets (Ld)26 labelled with anti-EgLIP1 and DAPI (middle panel) or pre-immune serum and DAPI as control (lower panel). Yellow arrow indicates anti-EgLIP1 signal (yellow spots) at border of Lds in mesocarp cell. Scale bar is 5 mm. (d) Lipase activity, lipid, EgLIP1 protein and transcript proles in mesocarp during fruit ripening. Fruit and early germinated seedlings were from the P1 progeny.

DAP: days after pollination. gs1, embryo from germinated seed; gs2, embryo from early post-germinated seed. Error bars represent s.d.

4 NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160 ARTICLE

To obtain more evidence for a role of EgLIP1 in the mesocarp, we studied its pattern of expression, protein accumulation and activity. Lipase activity assays revealed a 4100-fold increase during the maturation and ripening stages26 of fruit development (Fig. 3d), which coincided with the onset of lipid accumulation. The EgLIP1 transcript peak accumulation and the increase in the EgLIP1 protein (Fig. 3d, Supplementary Fig. S2) occurred concomitant to lipase activity. Quantitative RT-PCR analysis of transcripts in select palm tissues indicated that the EgLIP1 gene is almost exclusively expressed in the mesocarp during fruit ripening. EgLIP1 transcript was not detected in the embryo, kernel, owers, roots, leaves or the vegetative shoot apex (Supplementary Fig. S3). In early post-germinated seedlings, the EgLIP1 transcript was present, but at 34-fold less than that observed in the mesocarp (Fig. 3d, Supplementary Fig. S3C). To examine the relationships between EgLIP1 and the LL trait, lipase activity, transcript accumulation and the amounts of lipase protein were analyzed in the P2 segregating progeny (Fig. 4a, Supplementary Fig. S2, Supplementary Fig. S4). Low oil acidity was associated with LL activity, and barely detectable amounts of EgLIP1 protein and corresponding transcript, whereas high levels of acidity correlated with high levels of lipase activity, EgLIP1 protein and transcripts. These results provide strong evidence that oil acidity is directly related to the amount of transcript that encodes the EgLIP1 lipase.

The acidity trait cosegregates with mesocarp lipase gene. To determine the inheritance of the low-acidity trait, we analyzed the D1 self-crossed progeny (S). Two groups were observed, one with an average FFA content of 11.6% and the other with 48.5% (Fig. 4b). The low-acidity group represented 23.8% of the progeny, which is consistent with a 3:1 Mendelian segregation and a recessive low-acidity trait. We further analyzed three progenies derived from two successive D1 self-crosses (Fig. 4b). One contained only HL individuals (P1), another only LL individuals (P3) and the third progeny (P2) showed a 3:1 segregation pattern similar to that of S. These results conrmed that the oil acidity trait is monogenic with the low acidity being recessive, and that D1 is indeed heterozygous for the trait. Notably, similar results were also obtained from D1 L1 material (Supplementary Fig.

S5), which further indicated that the LL trait of D1 and L1 trees belonged to the same complementation group.

Analysis of genomic DNA from S and P2 revealed the existence of a (TG)8 (TA)16 microsatellite marker locus (mEgCIR_LIP03)

between 744 bp and 792 bp upstream of the EgLIP1 initiating ATG codon, with two alleles that showed a few base differences. A bulk segregant analysis (BSA), carried out in the P2 progeny showed a cosegregation of mEgCIR_LIP03 with the oil acidity trait. Two marker alleles (267 bp and 269 bp) were present in the bulk for HL activity and one single marker allele (267 bp) in the bulk for LL activity (Fig. 4c). Our BSA with 10 genotypes per bulk showed a very low probability (2 10 6) that this marker locus

Parent of P2 progeny

Signal intensity

Lipase activity (U g1 DW)

500

400

300

200

100

0 85 kDa

47 kDa

2

(arbitrary units)

267 269

BSA pool of high lipase P2

BSA pool of low lipase P2

Signal intensity

(arbitrary units)

2

1

1

1 4 6

3 5 7 8 9 Trees of P2 progeny

EgLip1

0

High Low

267 269

Acidity (% FFA)

80

60

40

20

0

2

0

Signal intensity

(arbitrary units)

2

1

0

267 269 Fragment sizes in base pair (bp)

S P1 P2 P3 Progeny

Figure 4 | The low-acidity trait is monogenic and cosegregates with EgLIP1 lipase gene. (a) Lipase activity and EgLIP1 protein in P2 progeny. Error bars represent s.d. (b) Acidity measurements from the S, P1, P2 and P3 progenies. Selng D1 led to a progeny of 21 trees (S), 3 of which were further self-crossed to give rise to progenies P1, P2 and P3 (see plant material in Material and Methods section). For each progeny, width of Boxplots is proportional to the number of trees per class. S and P2 showed a 3:1 ratio consistent with a monogenic recessive trait (S: w2 0.016, P 0.90, d.f. 1; P2: w2 0.205,

P 0.65, d.f. 1) (c) The low-acidity trait is linked to the EgLIP1 locus. The microsatellite marker EgCIR_LIP03 is located in the promoter region of EgLIP1

and shows two alleles that are visualized on electrophoregrams as 267 and 269-bp-long amplicons. Bulk segregation analysis was carried out on P2 progeny. Each genomic DNA pool (HL or LL) was represented by 10 randomly chosen trees in each acidity class.

NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160

is polymorphic without a physical link with the major gene controlling the oil acidity. The cosegregation provides further evidence that EgLIP1 is the gene that controls oil acidity in the pulp of mature fruits, or is at least closely linked to another gene controlling the activity of EgLIP1 itself. In either case, the mEgCirLIP03 marker or EgLIP1 itself provide useful markers for introgression of the LL trait into other elite oil palm genotypes.

DiscussionOur results indicate the EgLIP1 transcript is almost exclusively present in the mesocarp, suggesting that EgLIP1 function is restricted to this tissue (Supplementary Fig. S3). The oil contained in the mesocarp of oil palm fruit serves to attract disseminating animals and does not participate in the metabolism of other tissues and organs of the plant. One can speculate that the lipase function may enhance the palatability of fat, given this is due in part to receptors to fatty acids and not to triacylglycerol27, therefore increasing fruit attractiveness towards disseminating animals. As the main components of yield are not affected in LL trees (Table 1), the LL allele was evidently not selected against by breeders. As Deli genotypes originated from a narrow (4 trees) genetic basis, the LL allele is likely present in a signicant number of elite trees currently used for seed production. Indeed, a previous study allowed us to identify 4 LL trees out of 37 (ref. 11). We are currently screening a wide range of material for the LL allele to assess whether the microsatellite marker can be used in genotypes other than Deli and to increase our mapping population to achieve ne mapping of the trait. Preliminary mapping data conrm BSA results (N. Billotte, unpublished data).

The availability of a molecular marker for the trait will allow breeders to rapidly identify potent genitors that carry a LL allele, and to introgress the trait into main elite genotypes during the next breeding cycle for seed production. The low-lipase genotypes should allow more exibility in the management of harvest and bunch processing. Delaying harvest by 10 days in large estates should still reduce losses in the form of FFA, provide an additional 23% oil yield due to extended ripening and reduce harvesting frequency by one-third, thus lowering labour costs. Taken together, these impacts should represent an economic gain of B1 billion dollars per year, based on todays production.

Substantially higher gains are anticipated provided the LL trait is associated to non- or low-shedding fruits28. However, the most important benets are likely to be for African smallholders whose improved oil will gain access to international markets.

Methods

Plant material. All plant material were derived either from a cross between a tree (D1) from the Deli population and a tree (L1) from the La Me population or from progenies of D1 (Supplementary Fig. S6). The offspring of D1 L1 cross

corresponds to elite material planted worldwide. This cross was used to assess the performances of LL trees, unless otherwise indicated. Analyses were conducted using progenies from D1. S progeny resulted from self-crossing of D1 and contains 21 trees, 5 of which were further self-crossed to give rise to progenies P1 (15 trees), P2 (104 trees), P3 (17 trees), P4 (75 trees) and P5 (75 trees). These trees are cultivated both at CRA-PP (Pob, Bnin) and CEREPAH (La Dibamba, Cameroon). A bunch was considered to be ripe when one to ve fruits detached. Each experiment was performed 35 times.

Identication of LL and HL trees. The amount of endogenous lipase was estimated by measuring the percentage of FFA released from a mesocarp homogenate for 1 hour after tissue disruption10. Only attached, undamaged fruits that showed no visual signs of fungal infestation were used. Mesocarp (10 g) was ground with mortar and pestle and stored at room temperature for 1 hour.

The mixture was extracted with hexane (10 ml g 1 of pulp) and debris ltered away. Hexane was evaporated and the oil dried in a desiccator for 1 hour. Oil acidity was titrated with 0.1 N KOH according to standard procedures11. FFA content was calculated based on the molecular mass of palmitic acid.

Oil extraction from bunches. Oil was extracted using procedures similar to those used at mills. Bunches were subdivided into spikelets, which were boiled for 12 h; the fruits detached manually and the oil was extracted using a screw-press. After decantation, the oil was skimmed, dried, and its acidity was determined by titration. Data were BoxCox transformed to achieve a normal distribution, ascertained by AndersonDarling test. Analyses of variance were performed according to the general linear model. Multiple mean comparisons were made by Tukeys test.

Fungicide treatment. After harvest, each bunch was divided; one half was sprayed with a fungicide solution of DITHANE M45 (mancozeb 80%) at 2 g l 1 and the other was used as a control, and sprayed with water. Both samples were sun-dried and kept in a closed room at room temperature for 6, 9 and 12 days before the oil was extracted. Multiple mean comparisons were made by Tukeys test.

Mesocarp constituent analyses. For each stage of development, three independent bunches were collected from three individuals of the same progeny. Fifteen spikelets were collected in the centre of each bunch, from which ve spikelets were randomly sampled. Three undamaged fruits were then randomly sampled from each of the ve spikelets. Their mesocarp was dissected and ash-frozen in liquid nitrogen. Total lipids were extracted from 1-g samples of freeze-dried powder using a modied Folch method as described previously29.

For all electrophoresis gel analyses, total mesocarp proteins were extracted using aqueous phenol30. Labelling the lipase with radioactive tetrahydrolipstatin was carried out as described30, except that THL was dissolved in hexane and added to the hexane extract of mesocarp without using detergent. The gel was dried and submitted to autoradiography. The labelled band was submitted to trypsin digestion and analyzed by LC-MS/MS at Bordeaux Proteomic Facility as described previously30,31. Lipase activity was assayed by continuous titration of released fatty acids with a pHSTAT8. Antibodies were raised by Biotem France against the peptides RDYVNGGRLAKITGKYS and DRKKGFPKELSRND. The peptides were conjugated to a protein carrier and injected with Freund complete adjuvant. Rabbits were boosted twice and exsanguinated at day 45.

Total RNA from oil palm tissues was extracted as previously described32,33 and 1 mg was used to synthesize cDNA using a rst-strand cDNA synthesis kit (ImProm-II Reverse Transcription System, Promega). qPCR was conducted on a LightCycler 480 (Roche) in 96-well plates in a volume of 10 ml containing 2 ml of cDNA diluted 1/100, 1.5 ml of each primer (2 mM; 50-ACTACGTCAACGGTG GAAGG-30, 50-CCGTTCCAAATCCAGTTCTT-30) and 5 ml SYBR Green Mastermix (Roche). PCR was initiated by denaturation at 95 C for 10 min, followed by 45 cycles of 95 C for 15 s, 60 C for 15 s and a nal extension at 70 C for 1 min. Transcript levels were normalized to that of EgEfa1 (GenBank accession number: AY550990, oligonucleotides 50-TATCAAAGGATGGGCAGACC-30 and 50-TCATCATACCTTGCCTTGGA-30) mRNA, and relative mRNA abundance was determined as previously reported34. Control using RNA matrices were also conducted to validate the absence of DNA in each sample. Each time point was replicated three times from 35 independent biological samples, and all amplied cDNA fragments were sequenced by Beckman-Cogenics. Gene abundance is expressed as mean and standard error bars are calculated from three technical replicates.

Immunolocalization of EgLIP1 lipase. Mesocarp tissue samples were xed in0.1 M PBS buffer containing 8% (w/v) paraformaldehyde and 4% (w/v) N-(3-Dimethylaminopropyl)-N0-ethylcarbodiimide hydrochloride for 16 h at 4 C.

Samples were subsequently washed in 0.1 M PBS buffer and cut into 150 mm sections using a vibratome. Sections were saturated in 0.1 M PBS containing 3% (w/v) foetal bovine serum for 24 h at 4 C and transferred into 0.1 M PBS buffer containing either anti-EgLIP1 antibodies or pre-immune serum (diluted 1/100) for 24 h at 4 C. Sections were then washed 4 times in 0.1 M PBS and transferred into a0.1 M PBS solution containing goat anti-rabbit IgG antibody Alexa Fluor 545 conjugated (dilution 1/500) for 3 h. Sections were washed four times in PBS (0.1 M) and transferred to PBS buffer containing DAPI (14.3 mM) to visualize nucleus and membrane of oil droplets. Observation was performed using a confocal microscope (Zeiss LM510 meta) at magnications 40/1.2 numeric aperture, at the Cirad

imagery facility (Montpellier, France). Excitation was performed at 543 nm for Alexa Fluor 545 and at 405 nm for DAPI.

Cloning and sequencing of the EgLIP1 genomic region. A bacterial articial chromosome (BAC) library (92,160 clones, B6 genome equivalents) made from the oil palm genotype L1 was screened using EgLIP1 cDNA as a probe. Three positive BAC clones were sequenced (with at least 50-fold coverage) using the Roche/454 GS FLX sequencing technique and assembled by Beckman-Cogenics (USA). The BAC contigs were assembled into one single haplotype sequence using CodonCodeAligner version 3.7.1 (CodonCode Corporation, USA), although the tree used as source of DNA proved to be heterozygous for the lipase trait. A search of microsatellite loci in the EgLIP1 genomic region and the design of PCR primers for microsatellite locus mEgCirLIP03 were performed using the SSR Analysis Tool (SAT) software35.

6 NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160 ARTICLE

Bulk segregant analysis. Genomic DNA was extracted using the commercial kit NucleoSpin 96 Plant according to the manufacturer instructions (Macherey-Nagel, Germany). Legitimacy of P2 genotypes was conrmed by ADNid (Montpellier, France) company using 14 E. guineensis microsatellite markers36. A BSA was performed on the progeny P2 showing a 3:1 Mendelian segregation pattern for the oil acidity trait. The mEgCIR_LIP03 marker locus was amplied in the parent of P2 and in two equimolar mixtures of 10 DNAs each and made from P2 individuals randomly chosen within the two extreme phenotypic classes of oil acidity: (i) of HL phenotypes (oil acidity 440%) supposed to be homozygous ( ) or

heterozygous ( ) for EgLIP1, and (ii) of LL genotypes (acidity o20%)

supposed homozygous ( ) for that gene. The DNAs were PCR amplied in

condition of uorescence36,37 using the mEgCirLIP03 primers (forward: 50-TCAAA GAAACTTGTAGCATATATCAAA-30; reverse: 50-CATCCAGTAAGCTAACA CACAAATG-30) in a Mastercycler ep 384 thermocycler (Eppendorf, Germany) on 25 ng of DNA in a 10 ml nal volume of buffer (10 mM TrisHCl (pH 8), 200 mM of dNTP, 0.5 mM MgCl2, 0.05% w/v gelatin) containing 0.08 mM of the M13-labelled primer, 0.1 mM of the other primer and 1 U of Taq DNA polymerase. The PCR programme was: initial denaturation at 94 C for 4 min; 10 )Touchdown* cycles each of 94 C for 45 s, 55 C ( 0.5 C per cycle) for 1 min, 72 C for 1 min; and

then 25 cycles of 94 C for 45 s, 50 C for 60 s and 72 C for 60 s; and a nal elongation step at 72 C for 30 min. The PCR products were revealed by capillary electrophoresis on an Applied Biosystems 3500xL Genetic Analyzer (Life Technologies, USA) and the chromatogram allelic patterns were analyzed under the GeneMapper Software version 4.1 (Life Technologies, USA).

Determination of yield components. Ripe bunches were collected and analyzed as previously described38. Oil was extracted and quantied as described above. Results were analyzed using SAS v9.2 software (MIXED Procedure, SAS Institute, http://www.sas.com

Web End =http://www.sas.com ): a one-way analysis of variance (factor group) was carried out and the models residuals were studied to test the models validity (normality, outliers). The tree/bunch effect was tested with a likelihood ratio test. If signicant, it was kept in the model as a random effect. The results were expressed using a 95% condence interval and a margin of error of 5%.

References

1. Corley, R. H. V. & Tinker, P. B. The Oil Palm 4th edn (Blackwell Science, 2003).2. Vermeulen, S. & Goad, N. In Natural Resource Issues Series No. 5 (International Institute for Environment and Development, 2006).

3. Carlsson, A. S., Yilmaz, J. L., Green, A. G., Stymne, S. & Hofvander, P. Replacing fossil oil with fresh oilwith what and what for? Eur. J. Lipid Sci. Tech. 113, 812831 (2011).

4. Carter, C., Finley, W., Fry, J., Jackson, D. & Willis, L. Palm oil markets and future supply. Eur. J. Lipid Sci. Tech. 109, 307314 (2007).

5. Corley, R. H. V. How much palm oil do we need? Environ. Sci. Policy 12, 134139 (2009).

6. Murphy, D. J. Oil palm: future prospects for yield and quality improvements. Lipid Technology 21, 257260 (2009).

7. Henderson, J. & Osborne, D. J. Lipase activity in ripening and mature fruit of the oil palmstability in vivo and in vitro. Phytochemistry 30, 10731078 (1991).

8. Ngando-Ebongue, G. F., Dhouib, R., Carriere, F., Amvam Zollo, P. H. & Arondel, V. Assaying lipase activity from oil palm fruit (Elaeis guineensis Jacq.) mesocarp. Plant Physiol. Biochem. 44, 611617 (2006).

9. Sambanthamurthi, R., Oo, K. C. & Parman, S. H. Factors affecting lipase activity in Elaeis guineensis mesocarp. Plant Physiol. Biochem. 33, 353359 (1995).

10. Desassis, A. Palm oil acidication. Olagineux 12, 525534 (1957).11. Ngando Ebongue, G. F. et al. Identication of oil palm breeding genotypes producing an oil with low acid value. Eur. J. Lipid Sci. Tech. 110, 506509 (2008).

12. Fickendey, E. Uber die Verwendungsfahigkeit des Palmls als Speisefett. Tropenpanzer 14, 566570 (1910).

13. Gibon, V., De Greyt, W. & Kellens, M. Palm oil rening. Eur. J. Lipid Sci. Tech. 109, 315335 (2007).

14. Ohimain, E. I., Daokoru-Okukole, C., Isah, S. C. & Alaka, E. E. Assement of the quality of crude palm oil produced by smallholder processors in rivers state, Nigeria. Niger. J. Agri. Food Environ. 8, 2834 (2012).

15. Osei-Amponsah, C. et al. Processing practices of small-scale palm oil producers in the Kwaebibirem district, Ghana: a diagnostic study. NJAS-Wagen J. Life. Sci. 60-63, 4956 (2012).

16. Sambanthamurthi, R., Sundram, K. & Tan, Y. A. Chemistry and biochemistry of palm oil. Prog. Lipid Res. 39, 507558 (2000).

17. Mayes, S. et al. In Genomics of Tropical Crop Plants. (eds Moore, P. H. & Ming, R.) 371404 (Springer, 2008).

18. Rao, V., Gomez, M., Chayawat, N., Chia, C. C. & Lim, C. C. In Proceedings of the PIPOC 2001 International Palm Oil Congress (Agriculture) 167191.19. Corley, R. H. V. & Law, I. H. Ripening, Harvesting and Oil Extraction. The Planter 77, 507524 (2001).

20. Zu, K. S. A., Adjei-Nsiah, S. & Bani, R. J. Effect of processing equipment and duration of storage of palm fruits on palm oil yield and quality in the Kwaebibrem District, Ghana. Agric. Res. Rev. 1, 1825 (2012).

21. Tombs, M. P. & Stubbs, J. M. The absence of endogenous lipase from oil palm mesocarp. J. Sci. Food. Agr. 33, 892897 (1982).

22. Hadvary, P., Lengsfeld, H. & Wolfer, H. Inhibition of pancreatic lipasein vitro by the covalent inhibitor tetrahydrolipstatin. Biochem. J. 256, 357361 (1988).

23. Eastmond, P. J. Cloning and characterization of the acid lipase from castor beans. J. Biol. Chem. 279, 4554045545 (2004).

24. Bourgis, F. et al. Comparative transcriptome and metabolite analysis of oil palm and date palm mesocarp that differ dramatically in carbon partitioning. Proc. Natl Acad. Sci. USA 108, 1252712532 (2011).

25. Nurniwalis, A. W., Suhaimi, N., Siti Nor Akmar, A., Aminah, S. & Mohamad Arif, M. A. Gene discovery via expressed sequence tags from oil palm (Elaeis guineensis Jacq.). J. Oil Palm Res. 20, 8796 (2008).

26. Tranbarger, T. J. et al. Regulatory mechanisms underlying oil palm fruit mesocarp maturation, ripening, and functional specialization in lipid and carotenoid metabolism. Plant Physiol. 156, 564584 (2011).

27. Khan, N. A. & Besnard, P. Oro-sensory perception of dietary lipids: new insights into the fat taste transduction. Biochim. Biophys. Acta 1791, 149155 (2009).

28. Corley, R. H. V. Progress with non-shedding oil palms. Planter 82, 799807 (2006).

29. Laffargue, A., de Kochko, A. & Dussert, S. Development of solid-phase extraction and methylation procedures to analyse free fatty acids in lipid-rich seeds. Plant Physiol. Biochem. 45, 250257 (2007).

30. Dhouib, R. et al. Identication of a putative triacylglycerol lipase from papaya latex by functional proteomics. FEBS J. 278, 97110 (2011).

31. Deytieux, C. et al. Proteome analysis of grape skins during ripening. J. Exp. Bot. 58, 18511862 (2007).

32. Argout, X. et al. Towards the understanding of the cocoa transcriptome: production and analysis of an exhaustive dataset of ESTs of Theobroma cacao L. generated from various tissues and under various conditions. BMC Genomics 9, 512 (2008).

33. Morcillo, F. et al. Somaclonal variation in micropropagated oil palm. Characterization of two novel genes with enhanced expression in epigenetically abnormal cell lines and in response to auxin. Tree Physiol. 26, 585594 (2006).

34. Pfaf, M. W. A new mathematical model for relative quantication in real-time RT-PCR. Nucleic Acids Res. 29, e45 (2001).

35. Dereeper, A., Argout, X., Billot, C., Rami, J. F. & Ruiz, M. SAT, a exible and optimized Web application for SSR marker development. BMC Bioinformatics 8, 465 (2007).

36. Roy, R., Steffens, D. L., Gartside, B., Jang, G. Y. & Brumbaugh, J. A. Producing STR locus patterns from bloodstains and other forensic samples using an infrared uorescent automated DNA sequencer. J. Forensic Sci. 41, 418424 (1996).

37. Billotte, N. et al. Microsatellite-based high density linkage map in oil palm (Elaeis guineensis Jacq.). Theor. Appl. Genet. 110, 754765 (2005).

38. Blaak, G., Spamaaij, L. D. & Mendez, T. Breeding and inheritance in oil palm(E. guineensis Jacq.). II. Methods of bunch quality analysis. J. W. Afr. Inst. Oil Palm Res. 4, 145146 (1963).

Acknowledgements

This study was supported in part by PalmElit SAS, the Aquitaine region and by recurrent funding from IRD, Cirad, INRAB, IRAD, CNRS and University of Bordeaux. H.D. andG.-F.N.-E. were recipients of SCAC French embassy studentships, G.-F.N.-E. of an EMBO fellowship and F.B. of a post-doctoral fellowship from the Aquitaine region. We thank G. Piombo for helpful technical assistance and advice on oil acidity measurements and V. Poulain for statistical analysis of data. We especially thank S. Dussert for the lipid content determination, and A. Crabos and G. Bourdel for their help in RNA analysis. Radioactive tetrahydrolipstatin was obtained from Drs P. Hadvary and H. Lengsfeld (Roche) and the perl script for de novo sequencing from B. Valot (INRA, Ferme du Moulon). We are greatly indebted to J. Ohlrogge for helpful comments on the manuscript.

Author contributions

F.M., D.C., N.B., G.-F.N.-E., B.N. and V.A. designed the research. F.M., T.J.T., D.C., H.D., F.B., R.D., G.-F.N.-E., S.C., T.C., S.E., N.B., M.P. and V.A. performed research. V.A., F.M., D.C., N.B. and T.J.T. analyzed the data and wrote the paper. All authors discussed the results, and read and approved the nal submitted manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160

Competing nancial interests: The study was funded in part by PalmElit, a breeding company that might market low-lipase lines of oil palm in the future. Bruno Nouy is a PalmElit employee. All other authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Morcillo, F. et al. Improving palm oil quality through identication and mapping of the lipase gene causing oil deterioration. Nat. Commun. 4:2160 doi: 10.1038/ncomms3160 (2013).

This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/

Web End =http://creativecommons.org/licenses/by-nc-sa/3.0/

8 NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.