Microsatellite-based high density linkage map in oil palm (Elaeis guineensis Jacq.)

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).

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.

Development of Intron Polymorphism Markers and Their Association With Fatty Acid Component Variation in Oil Palm

Oil palm (Elaeis guineensis Jacq.) is a tropical woody oil crop of the palm family and is known as "the oil king of the world," but its palm oil contains about 50% palmitic acid, which is considered unhealthy for humans. Intron polymorphisms (IP) are highly efficient and easily examined molecular markers located adjacent to exon regions of functional genes, thus may be associated with targeted trait variation. In order to speed up the breeding of oil palm fatty acid composition, the current study identified a total of 310 introns located within 52 candidate genes involved in fatty acid biosynthesis in the oil palm genome. Based on the intron sequences, 205 primer pairs were designed, 64 of which showed polymorphism among 70 oil palm individuals. Phenotypic variation of fatty acid content in the 70 oil palm individuals was also investigated. Association analysis revealed that 13 IP markers were significantly associated with fatty acid content variation, and these IP markers were located on chromosomes 2, 5, 6, 8, 9, and 10 of oil palm. The development of such IP markers may be useful for the genetic improvement of fatty acid composition in oil palm.

Oil Palm Breeding in the Modern Era: Challenges and Opportunities

Oil palm, a cross-pollinated crop with long generation time, poses a lot of challenges in achieving sustainable oil palm with high yield and quality. The African oil palm (Elaeis guineensis Jacq.) is the most productive and versatile oil-yielding crop in the world, producing more than any other oil-yielding crop. Despite recent challenges, such as stress tolerance, superior oil quality, disease tolerance, and the need for new market niches, there is a growing need to explore and develop new varieties with high yield potential and the genetic diversity required to maintain oil palm yield stability. Breeding is an indispensable part of producing high-quality planting materials to increase oil palm yield. Biotechnological technologies have transformed conventional plant breeding approaches by introducing novel genotypes for breeding. Innovative pre-breeding and breeding approaches, such as identifying candidate genes in wild or land races using genomics tools, can pave the way for genetic improvement in oil palm. In this review, we highlighted the modern breeding tools, including genomics, marker-assisted breeding, genetic engineering, and genome editing techniques in oil palm crops, and we explored certain concerns connected to the techniques and their applications in practical breeding.

In silico QTL mapping in an oil palm breeding program reveals a quantitative and complex genetic resistance to Ganoderma boninense

Basal stem rot caused by Ganoderma boninense is the major threat to oil palm cultivation in Southeast Asia, which accounts for 80% of palm oil production worldwide, and this disease is increasing in Africa. The use of resistant planting material as part of an integrated pest management of this disease is one sustainable solution. However, breeding for Ganoderma resistance requires long-term and costly research, which could greatly benefit from marker-assisted selection (MAS). In this study, we evaluated the effectiveness of an in silico genetic mapping approach that took advantage of extensive data recorded in an ongoing breeding program. A pedigree-based QTL mapping approach applied to more than 10 years' worth of data collected during pre-nursery tests revealed the quantitative nature of Ganoderma resistance and identified underlying loci segregating in genetic diversity that is directly relevant for the breeding program supporting the study. To assess the consistency of QTL effects between pre-nursery and field environments, information was collected on the disease status of the genitors planted in genealogical gardens and modeled with pre-nursery-based QTL genotypes. In the field, individuals were less likely to be infected with Ganoderma when they carried more favorable alleles at the pre-nursery QTL. Our results pave the way for a MAS of Ganoderma resistant and high yielding planting material, and the provided proof-of-concept of this efficient and cost-effective approach could motivate similar studies based on diverse breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01246-9.

Assessment of genetic diversity and population structure of oil palm (Elaeis guineensis Jacq.) field genebank: A step towards molecular-assisted germplasm conservation

Oil palm (Elaeis guineensis) germplasm is exclusively maintained as ex situ living collections in the field for genetic conservation and evaluation. However, this is not for long term and the maintenance of field genebanks is expensive and challenging. Large area of land is required and the germplasms are exposed to extreme weather conditions and casualty from pests and diseases. By using 107 SSR markers, this study aimed to examine the genetic diversity and relatedness of 186 palms from a Nigerian-based oil palm germplasm and to identify core collection for conservation. On average, 8.67 alleles per SSR locus were scored with average effective number of alleles per population ranging from 1.96 to 3.34 and private alleles were detected in all populations. Mean expected heterozygosity was 0.576 ranging from 0.437 to 0.661 and the Wright's fixation index calculated was -0.110. Overall moderate genetic differentiation among populations was detected (mean pairwise population FST = 0.120, gene flow Nm = 1.117 and Nei's genetic distance = 0.466) and this was further confirmed by AMOVA analysis. UPGMA dendogram and Bayesian structure analysis concomitantly clustered the 12 populations into eight genetic groups. The best core collection assembled by Core Hunter ver. 3.2.1 consisted of 58 palms accounting for 31.2% of the original population, which was a smaller core set than using PowerCore 1.0. This core set attained perfect allelic coverage with good representation, high genetic distance between entries, and maintained genetic diversity and structure of the germplasm. This study reported the first molecular characterization and validation of core collections for oil palm field genebank. The established core collection via molecular approach, which captures maximum genetic diversity with minimum redundancy, would allow effective use of genetic resources for introgression and for sustainable oil palm germplasm conservation. The way forward to efficiently conserve the field genebanks into next generation without losing their diversity was further discussed.

Optimal set of microsatellite markers required to detect illegitimate progenies in selected oil palm (Elaeis guineensis Jacq.) breeding crosses

Oil palm is continually being improved via controlled crossing of selected palms to ensure sustainable yields and productivity. As such, correct parental assignment is important as the presence of illegitimates will compromise the progress of improvement. In the present study, we determined the optimal number of microsatellite (SSR) markers for detection of illegitimates in selected oil palm crosses with high confidence. Determining the optimal number of markers to assign parentage will ensure that the DNA fingerprinting will be cost effective for routine use as a quality control tool in oil palm improvement programs. Here, we evaluated a wide range of crosses that included a cross derived from wild germplasm palm. The results revealed that markers with high PIC are informative and detect most of the alleles present in a cross, including those exhibited by the illegitimates. A larger number of optimum sets of markers are needed to detect all illegitimates for crosses with higher levels of genetic diversity. The optimal number of polymorphic SSR markers determined in the present study can ensure that appropriate quality control is implemented for oil palm improvement programs.

Analysis of genetic diversity and discrimination of Oil Palm DxP populations based on the origins of pisifera elite parents

A total of 251 Dura cross Pisifera (DxP) hybrid palms from six populations descending from six parental African Pisifera origins and involving 12 progenies were analyzed with 19 selected Simple Sequence Repeats (SSR) markers. A total of 110 alleles were produced, ranging from three to eight per SSR, with a mean of 5.8 alleles per SSR locus. Of these, 68.5% were considered shared alleles by more than one population and the remaining 31.5% were population specific alleles. They generated between six and 21 haplotypes in all populations, and depending on the SSR marker, between one and 10 haplotypes within populations. Various parameters for analyzing genetic variability, differentiation and genetic structure were computed using GenAlEx, Structure and Darwin software. The obtained results confirmed microsatellites as a robust, feasible and trustful method for obtaining DNA fingerprints, tracing the source of oil palm samples. With respect to the authenticity of materials or for solving legitimacy issues, accession belonging to each population by SSR markers could be distinguished, but additional SSR should be screened for improving this process.

Genome-wide association study (GWAS) of major QTLs for bunch and oil yield related traits in Elaeis guineensis L

Oil palm (Elaeis guineensis Jacq.) is a long breeding cycle perennial crop with a genome size of 1.8 Gb. This is the first report of GWAS on large number of 310 African germplasm using 400 SSR markers till date. Highly significant correlation was found between leaf area (LA) and rachis length (RL) (0.75) followed by bunch weight (BW) and bunch index (BI) (0.65), whereas negative correlation was observed between bunch number (BN) and average bunch weight (ABW). First two principal component analysis (PCA) together explained maximum amount of variation (84.5 %). The PCA1 revealed that group 2 (Guinea Bissau and Cameroon) and group 4 (Zambia and Cameroon) genotypes are best suitable for BN, BI and BW traits. GWAS of six bunch yield and seven bunch oil yield traits with SSRs resulted in the identification 43 significant quantitative trait loci (QTLs) by mixed linear model (MLM) approach. Seven SSR loci were found to be linked to oil to dry mesocarp (ODM) on chromosomes 1,4,7,10,12 and 15. The SSR locus mEgCIR1753 for ODM was significantly linked at a p of ≤0.05 which explained 34.6 % of phenotypic variance. The important parameters like ODM, OWM and OB were located on 4, 10, 11 and 15 chromosomes. The leaf area and ODM were associated with candidate genes representing of low-temperature-induced 65 kDa proteins. The identified markers can be effectively used for marker assisted selection of high yielding oil palm genotypes.

Phenomics, genomics of oil palm (Elaeis guineensis Jacq.): way forward for making sustainable and high yielding quality oil palm

Oil palm (Elaeis guineensis Jacq.) is a heterogeneous, perennial crop having long breeding cycle with a genome size of 1.8 Gb. The demand for vegetable oil is steadily increasing, and expected that nearly 240-250 million tons of vegetable oil may be required by 2050. Genomics and next generation technologies plays crucial role in achieving the sustainable availability of oil palm with good yield and high quality. A successful breeding programme in oil palm depends on the availability of diverse gene pool, ex-situ conservation and their proper utilization for generating elite planting material. The major breeding methods adopted in oil palm are either modified recurrent selection or the modified reciprocal recurrent selection method. The QTLs of yield and related traits are chiefly located on chromosome 4, 10, 12 and 15 which is discussed in the current review. The probable chromosomal regions influencing the less height increment is observed to be on chromosomes 4, 10, 14 and 15. Advanced genomic approaches together with bioinformatics tools were discussed thoroughly for achieving sustainable oil palm where more efforts are needed. Major emphasis is given on oil palm crop improvement using holistic approaches of various genomic tools. Also a road map given on the milestones in the genomics and way forward for making oil palm to high yielding quality oil palm.

Molecular approaches for improving oil palm for oil

The oil palm, originating from Africa, is the most productive oil crop species. Palm oil is an important source of edible oil. Its current global plantation area is over 23 million ha. The theoretical oil yield potential of the oil palm is 18.2 tons/ha/year. However, current average oil yield is only 3.8 tons/ha/year. In the past 100 years, conventional breeding and improvement of field management played important roles in increasing oil yield. However, conventional breeding for trait improvement was limited by its very long (10-20 years) phenotypic selection cycle, although it improved oil yield by ~10-20% per generation. Molecular breeding using novel molecular technologies will accelerate genetic improvement and may reduce the need to deforest and to use arable land for expanding oil palm plantations, which in turn makes palm oil more sustainable. Here, we comprehensively synthesize information from relevant literature of the technologies, achievements, and challenges of molecular approaches, including tissue culture, haploid breeding, mutation breeding, marker-assisted selection (MAS), genomic selection (GS), and genome editing (GE). We propose the characteristics of ideal palms and suggest a road map to breed ideal palms for sustainable palm oil.

Transferability, development of simple sequence repeat (SSR) markers and application to the analysis of genetic diversity and population structure of the African fan palm (Borassus aethiopum Mart.) in Benin

Background

In Sub-Saharan Africa, Borassus aethiopum Mart. (African fan palm) is an important non-timber forest product-providing palm that faces multiple anthropogenic threats to its genetic diversity. However, this species is so far under-studied, which prevents its sustainable development as a resource. The present work is a first attempt at characterizing the genetic diversity and population structure of B. aethiopum across nine collection sites spanning the three climatic regions of Benin, West Africa, through the use of microsatellite markers.

Results

During a first phase we relied on the reported transferability of primers developed in other palm species. We find that, in disagreement with previously published results, only 22.5% of the markers tested enable amplification of B. aethiopum DNA and polymorphism detection is very low.

In a second phase, we generated a B. aethiopum-specific genomic dataset through high-throughput sequencing and used it for the de novo detection of microsatellite loci. Among the primer pairs targeting these, 11 detected polymorphisms and were further used for analyzing genetic diversity. Across the nine sites, expected heterozygosity (He) ranges from 0.263 to 0.451 with an overall average of 0.354, showing a low genetic diversity. Analysis of molecular variance (AMOVA) shows that within-site variation accounts for 53% of the genetic variation. Accordingly, the low number of migrants and positive values of the fixation index (F) in sites from both the Central (Sudano-Guinean) and the Southern (Guinean) climatic regions suggest limited gene flow between sites. The global correlation between genetic and geographic distances is weak; however, our clustering analyses indicate that B. aethiopum palms from Savè (Center) are genetically more similar to those from the North than to samples from other Central sites.

Conclusions

In the light of our results, we discuss the use of inter-species transfer vs. de novo development of microsatellite markers in genetic diversity analyses targeting under-studied species, and suggest future applications for our molecular resources. We propose that, while prominent short-range pollen and seed dispersal in Benin explain most of our results, gene flux between the Central and Northern regions, as a result of animal and/or human migrations, might underlie the Savè discrepancy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00955-y.

An Improved Oil Palm Genome Assembly as a Valuable Resource for Crop Improvement and Comparative Genomics in the Arecoideae Subfamily

Oil palm (Elaeis guineensis Jacq.) is the most traded crop among the economically important palm species. Here, we report an extended version genome of E. guineensis that is 1.2 Gb in length, an improvement of the physical genome coverage to 79% from the previous 43%. The improvement was made by assigning an additional 1968 originally unplaced scaffolds that were available publicly into the physical genome. By integrating three ultra-dense linkage maps and using them to place genomic scaffolds, the 16 pseudomolecules were extended. As we show, the improved genome has enhanced the mapping resolution for genome-wide association studies (GWAS) and permitted further identification of candidate genes/protein-coding regions (CDSs) and any non-coding RNA that may be associated with them for further studies. We then employed the new physical map in a comparative genomics study against two other agriculturally and economically important palm species—date palm (Phoenix dactylifera L.) and coconut palm (Cocos nucifera L.)—confirming the high level of conserved synteny among these palm species. We also used the improved oil palm genome assembly version as a palm genome reference to extend the date palm physical map. The improved genome of oil palm will enable molecular breeding approaches to expedite crop improvement, especially in the largest subfamily of Arecoideae, which consists of 107 species belonging to Arecaceae.

Genome-wide association study (GWAS) for morphological and yield-related traits in an oil palm hybrid (Elaeis oleifera x Elaeis guineensis) population

BACKGROUND

The genus Elaeis has two species of economic importance for the oil palm agroindustry: Elaeis oleifera (O), native to the Americas, and Elaeis guineensis (G), native to Africa. This work provides to our knowledge, the first association mapping study in an interspecific OxG oil palm population, which shows tolerance to pests and diseases, high oil quality, and acceptable fruit bunch production.

RESULTS

Using genotyping-by-sequencing (GBS), we identified a total of 3776 single nucleotide polymorphisms (SNPs) that were used to perform a genome-wide association analysis (GWAS) in 378 OxG hybrid population for 10 agronomic traits. Twelve genomic regions (SNPs) were located near candidate genes implicated in multiple functional categories, such as tissue growth, cellular trafficking, and physiological processes.

CONCLUSIONS

We provide new insights on genomic regions that mapped on candidate genes involved in plant architecture and yield. These potential candidate genes need to be confirmed for future targeted functional analyses. Associated markers to the traits of interest may be valuable resources for the development of marker-assisted selection in oil palm breeding.

Linkage-based genome assembly improvement of oil palm (Elaeis guineensis)

Meiotic crossovers in outbred species, such as oil palm (Elaeis guineensis Jacq., 2n = 32) contribute to allelic re-assortment in the genome. Such genetic variation is usually exploited in breeding to combine positive alleles for trait superiority. A good quality reference genome is essential for identifying the genetic factors underlying traits of interest through linkage or association studies. At the moment, an AVROS pisifera genome is publicly available for oil palm. Distribution and frequency of crossovers throughout chromosomes in different origins of oil palm are still unclear. Hence, an ultrahigh-density genomic linkage map of a commercial Deli dura x AVROS pisifera family was constructed using the OP200K SNP array, to evaluate the genetic alignment with the genome assembly. A total of 27,890 linked SNP markers generated a total map length of 1,151.7 cM and an average mapping interval of 0.04 cM. Nineteen linkage groups represented 16 pseudo-chromosomes of oil palm, with 61.7% of the mapped SNPs present in the published genome. Meanwhile, the physical map was also successfully extended from 658 Mb to 969 Mb by assigning unplaced scaffolds to the pseudo-chromosomes. A genic linkage map with major representation of sugar and lipid biosynthesis pathways was subsequently built for future studies on oil related quantitative trait loci (QTL). This study improves the current physical genome of the commercial oil palm, and provides important insights into its recombination landscape, eventually unlocking the full potential genome sequence-enabled biology for oil palm.

Development and validation of whole genome-wide and genic microsatellite markers in oil palm (Elaeis guineensis Jacq.): First microsatellite database (OpSatdb)

The availability of large expressed sequence tag (EST) and whole genome databases of oil palm enabled the development of a data base of microsatellite markers. For this purpose, an EST database consisting of 40,979 EST sequences spanning 27 Mb and a chromosome-wise whole genome databases were downloaded. A total of 3,950 primer pairs were identified and developed from EST sequences. The tri and tetra nucleotide repeat motifs were most prevalent (each 24.75%) followed by di-nucleotide repeat motifs. Whole genome-wide analysis found a total of 245,654 SSR repeats across the 16 chromosomes of oil palm, of which 38,717 were compound microsatellite repeats. A web application, OpSatdb, the first microsatellite database of oil palm, was developed using the PHP and MySQL database ( https://ssr.icar.gov.in/index.php ). It is a simple and systematic web-based search engine for searching SSRs based on repeat motif type, repeat type, and primer details. High synteny was observed between oil palm and rice genomes. The mapping of ESTs having SSRs by Blast2GO resulted in the identification of 19.2% sequences with gene ontology (GO) annotations. Randomly, a set of ten genic SSRs and five genomic SSRs were used for validation and genetic diversity on 100 genotypes belonging to the world oil palm genetic resources. The grouping pattern was observed to be broadly in accordance with the geographical origin of the genotypes. The identified genic and genome-wide SSRs can be effectively useful for various genomic applications of oil palm, such as genetic diversity, linkage map construction, mapping of QTLs, marker-assisted selection, and comparative population studies.

Cloning and characterization of EgGDSL, a gene associated with oil content in oil palm

Oil palm (Elaeis guineensis, Jacq.) is a key tropical oil crop, which provides over one third of the global vegetable oil production, but few genes related to oil yield have been characterized. In this study, a GDSL esterase/lipase gene, which was significantly associated with oil content, was isolated from oil palm and designated as EgGDSL. Its functional characterization was carried out through ectopic expression in Arabidopsis ecotype Col-0. It was shown that expression of EgGDSL in Arabidopsis led to the increased total fatty acid content by 9.5% compared with the wild type. Further analysis of the fatty acid composition revealed that stearic acid (18:0) increased in the seeds of the transgenic lines, but the levels of linoleic acid (18:2) plus 11-eicosenoic acid drastically declined. Quantitative real-time PCR (qPCR) revealed that in oil palm, EgGDSL was highly expressed in mesocarp followed by leaf, and the expression level was very low in the root. The expression level of EgGDSL gene began to increase at two months after flowering (MAF) and reached its peak by four MAF, then declined rapidly, and reached its lowest level during the mature period (6 MAF). The EgGDSL gene was more highly expressed in oil palm trees with high oil content than that with low oil content, demonstrating that the transcription level of EgGDSL correlated with the amount of oil accumulation. The gene may be valuable for engineering fatty acid metabolism in crop improvement programmes and for marker-assisted breeding.

Developing genome-wide SNPs and constructing an ultrahigh-density linkage map in oil palm

Oil palm (Elaeis guineensis Jacq.) is the leading oil-producing crops and the most important edible oil resource worldwide. DNA markers and genetic linkage maps are essential resources for marker-assisted selection to accelerate genetic improvement. We conducted RAD-seq on an Illumina NextSeq500 to discover genome-wide SNPs, and used the SNPs to construct a linkage map for an oil palm (Tenera) population derived from a cross between a Deli Dura and an AVROS Pisifera. The RAD-seq produced 1,076 million single-end reads across the breeding population containing 155 trees. Mining this dataset detected 510,251 loci. After filtering out loci with low accuracy and more than 20% missing data, 11,394 SNPs were retained. Using these SNPs, in combination with 188 anchor SNPs and 123 microsatellites, we constructed a linkage map containing 10,023 markers covering 16 chromosomes. The map length is 2,938.2 cM with an average marker space of 0.29 cM. The large number of SNPs will supply ample choices of DNA markers in analysing the genetic diversity, population structure and evolution of oil palm. This high-density linkage map will contribute to mapping quantitative trait loci (QTL) for important traits, thus accelerating oil palm genetic improvement.

High density SNP and DArT-based genetic linkage maps of two closely related oil palm populations

Oil palm (Elaeis guineensis Jacq.) is an outbreeding perennial tree crop with long breeding cycles, typically 12 years. Molecular marker technologies can greatly improve the breeding efficiency of oil palm. This study reports the first use of the DArTseq platform to genotype two closely related self-pollinated oil palm populations, namely AA0768 and AA0769 with 48 and 58 progeny respectively. Genetic maps were constructed using the DArT and SNP markers generated in combination with anchor SSR markers. Both maps consisted of 16 major independent linkage groups (2n = 2× = 32) with 1399 and 1466 mapped markers for the AA0768 and AA0769 populations, respectively, including the morphological trait "shell-thickness" (Sh). The map lengths were 1873.7 and 1720.6 cM with an average marker density of 1.34 and 1.17 cM, respectively. The integrated map was 1803.1 cM long with 2066 mapped markers and average marker density of 0.87 cM. A total of 82% of the DArTseq marker sequence tags identified a single site in the published genome sequence, suggesting preferential targeting of gene-rich regions by DArTseq markers. Map integration of higher density focused around the Sh region identified closely linked markers to the Sh, with D.15322 marker 0.24 cM away from the morphological trait and 5071 bp from the transcriptional start of the published SHELL gene. Identification of the Sh marker demonstrates the robustness of using the DArTseq platform to generate high density genetic maps of oil palm with good genome coverage. Both genetic maps and integrated maps will be useful for quantitative trait loci analysis of important yield traits as well as potentially assisting the anchoring of genetic maps to genomic sequences.

Genetic evidence of multiple invasions and a small number of founders of Asian Palmyra palm (Borassus flabellifer) in Thailand

BACKGROUND

Borassus flabellifer or Asian Palmyra palm is an important crop for local economies in the South and Southeast Asia for its fruit and palm sugar production. Archeological and historical evidence indicated the presence of this species in Southeast Asia dating back at least 1500 years. B. flabellifer is believed to be originated in Africa, spread to South Asia and introduced into Southeast Asia through commercial routes and dissemination of cultures, however, the nature of its invasion and settlement in Thailand is unclear.

RESULTS

Here, we analyzed genetic data of 230 B. flabellifer accessions across Thailand using 17 EST-SSR and 12 gSSR polymorphic markers. Clustering analysis revealed that the population consisted of two genetic clusters (STRUCTURE K = 2). Cluster I is found mainly in southern Thailand, while Cluster II is found mainly in the northeastern. Those found in the central are of an extensive mix between the two. These two clusters are in moderate differentiation (F _ST = 0.066 and N _M = 3.532) and have low genetic diversity (H_O = 0.371 and 0.416; A_R = 2.99 and 3.19, for the cluster I and II respectively). The minimum numbers of founders for each genetic group varies from 3 to 4 individuals, based on simulation using different allele frequency assumptions. These numbers coincide with that B. flabellifer is dioecious, and a number of seeds had to be simultaneously introduced for obtaining both male and female founders.

CONCLUSIONS

From these data and geographical and historical evidence, we hypothesize that there were at least two different invasive events of B. flabellifer in Thailand. B. flabellifer was likely brought through the Straits of Malacca to be propagated in the southern Thailand as one of the invasive events before spreading to the central Thailand. The second event likely occurred in Khmer Empire, currently Cambodia, before spreading to the northeastern Thailand.

Identification of Ganoderma Disease Resistance Loci Using Natural Field Infection of an Oil Palm Multiparental Population

Multi-parental populations are promising tools for identifying quantitative disease resistance loci. Stem rot caused by Ganoderma boninense is a major threat to palm oil production, with yield losses of up to 80% prompting premature replantation of palms. There is evidence of genetic resistance sources, but the genetic architecture of Ganoderma resistance has not yet been investigated. This study aimed to identify Ganoderma resistance loci using an oil palm multi-parental population derived from nine major founders of ongoing breeding programs. A total of 1200 palm trees of the multi-parental population was planted in plots naturally infected by Ganoderma, and their health status was assessed biannually over 25 yr. The data were treated as survival data, and modeled using the Cox regression model, including a spatial effect to take the spatial component in the spread of Ganoderma into account. Based on the genotypes of 757 palm trees out of the 1200 planted, and on pedigree information, resistance loci were identified using a random effect with identity-by-descent kinship matrices as covariance matrices in the Cox model. Four Ganoderma resistance loci were identified, two controlling the occurrence of the first Ganoderma symptoms, and two the death of palm trees, while favorable haplotypes were identified among a major gene pool for ongoing breeding programs. This study implemented an efficient and flexible QTL mapping approach, and generated unique valuable information for the selection of oil palm varieties resistant to Ganoderma disease.

Genome-wide identification of markers for selecting higher oil content in oil palm

BACKGROUND

Oil palm (Elaeis guineensis, Jacq.) is the most important source of edible oil. The improvement of oil yield is currently slow in conventional breeding programs due to long generation intervals. Marker-assisted selection (MAS) has the potential to accelerate genetic improvement. To identify DNA markers associated with oil content traits for MAS, we performed quantitative trait loci (QTL) mapping using genotyping by sequencing (GBS) in a breeding population derived from a cross between Deli Dura and Ghana Pisifera, containing 153 F₁ trees.

RESULTS

We constructed a high-density linkage map containing 1357 SNPs and 123 microsatellites. The 16 linkage groups (LGs) spanned 1527 cM, with an average marker space of 1.03 cM. One significant and three suggestive QTL for oil to bunch (O/B) and oil to dry mesocarp (O/DM) were mapped on LG1, LG8, and LG10 in a F₁ breeding population, respectively. These QTL explained 7.6-13.3% of phenotypic variance. DNA markers associated with oil content in these QTL were identified. Trees with beneficial genotypes at two QTL for O/B showed an average O/B of 30.97%, significantly (P < 0.01) higher than that of trees without any beneficial QTL genotypes (average O/B of 28.24%). QTL combinations showed that the higher the number of QTL with beneficial genotypes, the higher the resulting average O/B in the breeding population.

CONCLUSIONS

A linkage map with 1480 DNA markers was constructed and used to identify QTL for oil content traits. Pyramiding the identified QTL with beneficial genotypes associated with oil content traits using DNA markers has the potential to accelerate genetic improvement for oil yield in the breeding population of oil palm.

Genetic diversity of Elaeis oleifera (HBK) Cortes populations using cross species SSRs: implication's for germplasm utilization and conservation

Background

The Elaeis oleifera genetic materials were assembled from its center of diversity in South and Central America. These materials are currently being preserved in Malaysia as ex situ living collections. Maintaining such collections is expensive and requires sizable land. Information on the genetic diversity of these collections can help achieve efficient conservation via maintenance of core collection. For this purpose, we have applied fourteen unlinked microsatellite markers to evaluate 532 E. oleifera palms representing 19 populations distributed across Honduras, Costa Rica, Panama and Colombia.

Results

In general, the genetic diversity decreased from Costa Rica towards the north (Honduras) and south-east (Colombia). Principle coordinate analysis (PCoA) showed a single cluster indicating low divergence among palms. The phylogenetic tree and STRUCTURE analysis revealed clusters based on country of origin, indicating considerable gene flow among populations within countries. Based on the values of the genetic diversity parameters, some genetically diverse populations could be identified. Further, a total of 34 individual palms that collectively captured maximum allelic diversity with reduced redundancy were also identified. High pairwise genetic differentiation (Fst > 0.250) among populations was evident, particularly between the Colombian populations and those from Honduras, Panama and Costa Rica. Crossing selected palms from highly differentiated populations could generate off-springs that retain more genetic diversity.

Conclusion

The results attained are useful for selecting palms and populations for core collection. The selected materials can also be included into crossing scheme to generate offsprings that capture greater genetic diversity for selection gain in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0505-7) contains supplementary material, which is available to authorized users.

Development, identification and validation of CAPS marker for SHELL trait which governs dura, pisifera and tenera fruit forms in oil palm (Elaeis guineensis Jacq.)

The oil palm fruit forms (dura, pisifera and tenera) governed by the shell thickness gene (Sh) plays a major role in identification of fruit type and also influences palm oil yield. Identification of desired fruit type is a major asset to the breeders and oil palm workers for applications in breeding, seed certification and to reduce time, space and money spent on identification of fruit form. In the present study, we developed Sh gene specific primer pairs and bulk segregant analysis was done using 300 genomic and 8 genic SSR markers. We identified one cleaved amplified polymorphic site (CAPS) marker for differentiation of oil palm fruit type which produced two alleles (280 and 250bp) in dura genotypes, three alleles in tenera genotypes (550, 280, and 250bp) and one allele in pisifera genotypes (550bp). The shell allele sequencing results showed that two SNPs were present, of which SNP2 contributed for variation of fruit forms. The nucleotide ‘A’ was present in only dura genotypes, where as ‘T’ was present only in pisifera genotypes, which in turn led to the change of amino acid lysine to aspargine. The identified CAPS marker was validated on 300 dura, 25 pisifera and 80 tenera genotypes, 80 dura/ pisifera cross progenies and 60 lines of tenera/ tenera cross progeny. Association mapping of marker data with phenotypic data of eight oil yield related traits resulted in identification of seven significant QTLs by GLM approach, four by MLM approach at a significant threshold (P) level of 0.001. Significant QTLs were identified for fruit to bunch and oil to bunch traits, which explained R² of 12.9% and 11.5% respectively. The CAPS marker used in the present study facilitate selection and timely distribution of desirable high yielding tenera sprouts to the farmers instead of waiting for 4–5 years. This saves a lot of land, time and money which will be a major breakthrough to the oil palm community.

Draft genome sequence of an elite Dura palm and whole-genome patterns of DNA variation in oil palm

Oil palm is the world’s leading source of vegetable oil and fat. Dura, Pisifera and Tenera are three forms of oil palm. The genome sequence of Pisifera is available whereas the Dura form has not been sequenced yet. We sequenced the genome of one elite Dura palm, and re-sequenced 17 palm genomes. The assemble genome sequence of the elite Dura tree contained 10,971 scaffolds and was 1.701 Gb in length, covering 94.49% of the oil palm genome. 36,105 genes were predicted. Re-sequencing of 17 additional palm trees identified 18.1 million SNPs. We found high genetic variation among palms from different geographical regions, but lower variation among Southeast Asian Dura and Pisifera palms. We mapped 10,000 SNPs on the linkage map of oil palm. In addition, high linkage disequilibrium (LD) was detected in the oil palms used in breeding populations of Southeast Asia, suggesting that LD mapping is likely to be practical in this important oil crop. Our data provide a valuable resource for accelerating genetic improvement and studying the mechanism underlying phenotypic variations of important oil palm traits.

Mating System and Genetic Composition of the Macaw Palm (Acrocomia aculeata): Implications for Breeding and Genetic Conservation Programs

Acrocomia aculeata (Arecaceae), a palm endemic to South and Central America, is a potential oil crop. Knowledge of the mating system of this species is limited to its reproductive biology and to studies using molecular markers. The present study analyzed genetic diversity between its developmental stages and determined its prevailing mating system in order to support genetic conservation and breeding programs. We tested 9 microsatellite markers in 27 mother trees (adult plants) and 157 offspring (juvenile plants) from the southeastern region of Brazil. Heterozygosity levels differed between the 2 studied life stages, as indicated by the fixation index of adult and juvenile trees, suggesting that selection against homozygotes occurs during the plant life cycle. The mating system parameters analyzed indicate that A. aculeata is predominantly outcrossing (allogamous). However, its low levels of selfing suggest that there is individual variation with regard to self-incompatibility, which can be a survival strategy in isolated or fragmented habitats. Deviations in variance effective size were detected because of high mating rates among relatives and correlated matings. These findings indicate that the main source of inbreeding results from biparental inbreeding in the population and that the progenies are predominantly composed of full-sibs. The information provided by this study on the ecology and reproduction dynamics of A. aculeata should be useful to both breeding and genetic conservation programs, allowing the development of more precise mathematical models and the estimation of the appropriate number of mother trees for seed collection.

Mixed model approach for IBD-based QTL mapping in a complex oil palm pedigree

Background

Elaeis guineensis is the world’s leading source of vegetable oil, and the demand is still increasing. Oil palm breeding would benefit from marker-assisted selection but genetic studies are scarce and inconclusive. This study aims to identify genetic bases of oil palm production using a pedigree-based approach that is innovative in plant genetics.

Results

A quantitative trait locus (QTL) mapping approach involving two-step variance component analysis was employed using phenotypic data on 30852 palms from crosses between more than 300 genotyped parents of two heterotic groups. Genome scans were performed at parental level by modeling QTL effects as random terms in linear mixed models with identity-by-descent (IBD) kinship matrices. Eighteen QTL regions controlling production traits were identified among a large genetically diversified sample from breeding program. QTL patterns depended on the genetic origin, with only one region shared between heterotic groups. Contrasting effects of QTLs on bunch number and weights reflected the close negative correlation between the two traits.

Conclusions

The pedigree-based approach using data from ongoing breeding programs is a powerful, relevant and economic approach to map QTLs. Genetic determinisms contributing to heterotic effects have been identified and provide valuable information for orienting oil palm breeding strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1985-3) contains supplementary material, which is available to authorized users.

Genomic selection prediction accuracy in a perennial crop: case study of oil palm (Elaeis guineensis Jacq.)

Genomic selection empirically appeared valuable for reciprocal recurrent selection in oil palm as it could account for family effects and Mendelian sampling terms, despite small populations and low marker density. Genomic selection (GS) can increase the genetic gain in plants. In perennial crops, this is expected mainly through shortened breeding cycles and increased selection intensity, which requires sufficient GS accuracy in selection candidates, despite often small training populations. Our objective was to obtain the first empirical estimate of GS accuracy in oil palm (Elaeis guineensis), the major world oil crop. We used two parental populations involved in conventional reciprocal recurrent selection (Deli and Group B) with 131 individuals each, genotyped with 265 SSR. We estimated within-population GS accuracies when predicting breeding values of non-progeny-tested individuals for eight yield traits. We used three methods to sample training sets and five statistical methods to estimate genomic breeding values. The results showed that GS could account for family effects and Mendelian sampling terms in Group B but only for family effects in Deli. Presumably, this difference between populations originated from their contrasting breeding history. The GS accuracy ranged from -0.41 to 0.94 and was positively correlated with the relationship between training and test sets. Training sets optimized with the so-called CDmean criterion gave the highest accuracies, ranging from 0.49 (pulp to fruit ratio in Group B) to 0.94 (fruit weight in Group B). The statistical methods did not affect the accuracy. Finally, Group B could be preselected for progeny tests by applying GS to key yield traits, therefore increasing the selection intensity. Our results should be valuable for breeding programs with small populations, long breeding cycles, or reduced effective size.

Genome-wide SNP discovery and identification of QTL associated with agronomic traits in oil palm using genotyping-by-sequencing (GBS)

Oil palm has become one of the most important oil crops in the world. Marker-assisted selections have played a pivotal role in oil palm breeding programs. Here, we report the use of genotyping-by-sequencing (GBS) approach for a large-scale SNP discovery and genotyping of a mapping population. Reduced representation libraries of 108 F2 progeny were sequenced and a total of 524 million reads were obtained. We detected 21,471 single nucleotide substitutions, most of which (62.6%) represented transition events. Of 3417 fully informative SNP markers, we were able to place 1085 on a linkage map, which spanned 1429.6 cM and had an average of one marker every 1.26 cM. Three QTL affecting trunk height were detected on LG 10, 14 and 15, whereas a single QTL associated with fruit bunch weight was identified on LG 3. The use of GBS approach proved to be rapid, cost-effective and highly reproducible in this species.

A consensus linkage map of oil palm and a major QTL for stem height

Oil palm (Elaeis guinensis Jacquin) is the most important source of vegetable oil and fat. Several linkage maps had been constructed using dominant and co-dominant markers to facilitate mapping of QTL. However, dominant markers are not easily transferable among different laboratories. We constructed a consensus linkage map for oil palm using co-dominant markers (i.e. microsatellite and SNPs) and two F1 breeding populations generated by crossing Dura and Pisifera individuals. Four hundreds and forty-four microsatellites and 36 SNPs were mapped onto 16 linkage groups. The map length was 1565.6 cM, with an average marker space of 3.72 cM. A genome-wide scan of QTL identified a major QTL for stem height on the linkage group 5, which explained 51% of the phenotypic variation. Genes in the QTL were predicted using the palm genome sequence and bioinformatic tools. The linkage map supplies a base for mapping QTL for accelerating the genetic improvement, and will be also useful in the improvement of the assembly of the genome sequences. Markers linked to the QTL may be used in selecting dwarf trees. Genes within the QTL will be characterized to understand the mechanisms underlying dwarfing.

Exploiting transcriptome data for the development and characterization of gene-based SSR markers related to cold tolerance in oil palm (Elaeis guineensis)

BACKGROUND

The oil palm (Elaeis guineensis, 2n = 32) has the highest oil yield of any crop species, as well as comprising the richest dietary source of provitamin A. For the tropical species, the best mean growth temperature is about 27°C, with a minimal growth temperature of 15°C. Hence, the plantation area is limited into the geographical ranges of 10°N to 10°S. Enhancing cold tolerance capability will increase the total cultivation area and subsequently oil productivity of this tropical species. Developing molecular markers related to cold tolerance would be helpful for molecular breeding of cold tolerant Elaeis guineensis.

RESULTS

In total, 5791 gene-based SSRs were identified in 51,452 expressed sequences from Elaeis guineensis transcriptome data: approximately one SSR was detected per 10 expressed sequences. Of these 5791 gene-based SSRs, 916 were derived from expressed sequences up- or down-regulated at least two-fold in response to cold stress. A total of 182 polymorphic markers were developed and characterized from 442 primer pairs flanking these cold-responsive SSR repeats. The polymorphic information content (PIC) of these polymorphic SSR markers across 24 lines of Elaeis guineensis varied from 0.08 to 0.65 (mean = 0.31 ± 0.12). Using in-silico mapping, 137 (75.3%) of the 182 polymorphic SSR markers were located onto the 16 Elaeis guineensis chromosomes. Total coverage of 473 Mbp was achieved, with an average physical distance of 3.4 Mbp between adjacent markers (range 96 bp - 20.8 Mbp). Meanwhile, Comparative analysis of transcriptome under cold stress revealed that one ICE1 putative ortholog, five CBF putative orthologs, 19 NAC transcription factors and four cold-induced orhologs were up-regulated at least two fold in response to cold stress. Interestingly, 5' untranslated region of both Unigene21287 (ICE1) and CL2628.Contig1 (NAC) both contained an SSR markers.

CONCLUSIONS

In the present study, a series of SSR markers were developed based on sequences differentially expressed in response to cold stress. These EST-SSR markers would be particularly useful for gene mapping and population structure analysis in Elaeis guineensis. Meanwhile, the EST-SSR loci were inducible expressed in response to low temperature, which may have potential application in identifying trait-associated markers in oil palm in the future.

Illegitimacy and sibship assignments in oil palm (Elaeis guineensis Jacq.) half-sib families using single locus DNA microsatellite markers

Oil palm breeding has been progressing very well in Southeast Asia, especially in Malaysia and Indonesia. Despite this progress, there are still problems due to the difficulty of controlled crossing in oil palm. Contaminated/illegitimate progeny has appeared in some breeding programs; late and failure of detection by the traditional method causes a waste of time and labor. The use of molecular markers improves the integrity of breeding programs in perennial crops such as oil palm. Four half-sib families with a total of 200 progeny were used in this study. Thirty polymorphic single locus DNA microsatellites markers were typed to identify the illegitimate individuals and to obtain the correct parental and progeny assignments by using the CERVUS and COLONY programs. Three illegitimate palms (1.5%) were found, and 16 loci proved to be sufficient for sibship assignments without parental genotypes by using the COLONY program. The pairwise-likelihood score (PLS) method was better for half-sib family assignments than the full likelihood (FL) method.

Genetic architecture of palm oil fatty acid composition in cultivated oil palm (Elaeis guineensis Jacq.) compared to its wild relative E. oleifera (H.B.K) Cortés

We searched for quantitative trait loci (QTL) associated with the palm oil fatty acid composition of mature fruits of the oil palm E. guineensis Jacq. in comparison with its wild relative E. oleifera (H.B.K) Cortés. The oil palm cross LM2T x DA10D between two heterozygous parents was considered in our experiment as an intraspecific representative of E. guineensis. Its QTLs were compared to QTLs published for the same traits in an interspecific Elaeis pseudo-backcross used as an indirect representative of E. oleifera. Few correlations were found in E. guineensis between pulp fatty acid proportions and yield traits, allowing for the rather independent selection of both types of traits. Sixteen QTLs affecting palm oil fatty acid proportions and iodine value were identified in oil palm. The phenotypic variation explained by the detected QTLs was low to medium in E. guineensis, ranging between 10% and 36%. The explained cumulative variation was 29% for palmitic acid C16:0 (one QTL), 68% for stearic acid C18:0 (two QTLs), 50% for oleic acid C18:1 (three QTLs), 25% for linoleic acid C18:2 (one QTL), and 40% (two QTLs) for the iodine value. Good marker co-linearity was observed between the intraspecific and interspecific Simple Sequence Repeat (SSR) linkage maps. Specific QTL regions for several traits were found in each mapping population. Our comparative QTL results in both E. guineensis and interspecific materials strongly suggest that, apart from two common QTL zones, there are two specific QTL regions with major effects, which might be one in E. guineensis, the other in E. oleifera, which are independent of each other and harbor QTLs for several traits, indicating either pleiotropic effects or linkage. Using QTL maps connected by highly transferable SSR markers, our study established a good basis to decipher in the future such hypothesis at the Elaeis genus level.

High density SNP and SSR-based genetic maps of two independent oil palm hybrids

BACKGROUND

Oil palm is an important perennial oil crop with an extremely long selection cycle of 10 to 12 years. As such, any tool that speeds up its genetic improvement process, such as marker-assisted breeding is invaluable. Previously, genetic linkage maps based on AFLP, RFLP and SSR markers were developed and QTLs for fatty acid composition and yield components identified. High density genetic maps of crosses of different genetic backgrounds are indispensable tools for investigating oil palm genetics. They are also useful for comparative mapping analyses to identify markers closely linked to traits of interest.

RESULTS

A 4.5 K customized oil palm SNP array was developed using the Illumina Infinium platform. The SNPs and 252 SSRs were genotyped on two mapping populations, an intraspecific cross with 87 palms and an interspecific cross with 108 palms. Parental maps with 16 linkage groups (LGs), were constructed for the three fruit forms of E. guineensis (dura, pisifera and tenera). Map resolution was further increased by integrating the dura and pisifera maps into an intraspecific integrated map with 1,331 markers spanning 1,867 cM. We also report the first map of a Colombian E. oleifera, comprising 10 LGs with 65 markers spanning 471 cM. Although not very dense due to the high level of homozygosity in E. oleifera, the LGs were successfully integrated with the LGs of the tenera map. Direct comparison between the parental maps identified 603 transferable markers polymorphic in at least two of the parents. Further analysis revealed a high degree of marker transferability covering 1,075 cM, between the intra- and interspecific integrated maps. The interspecific cross displayed higher segregation distortion than the intraspecific cross. However, inclusion of distorted markers in the genetic maps did not disrupt the marker order and no map expansion was observed.

CONCLUSIONS

The high density SNP and SSR-based genetic maps reported in this paper have greatly improved marker density and genome coverage in comparison with the first reference map based on AFLP and SSR markers. Therefore, it is foreseen that they will be more useful for fine mapping of QTLs and whole genome association mapping studies in oil palm.

A first genetic map of date palm (Phoenix dactylifera) reveals long-range genome structure conservation in the palms

Background

The date palm is one of the oldest cultivated fruit trees. It is critical in many ways to cultures in arid lands by providing highly nutritious fruit while surviving extreme heat and environmental conditions. Despite its importance from antiquity, few genetic resources are available for improving the productivity and development of the dioecious date palm. To date there has been no genetic map and no sex chromosome has been identified.

Results

Here we present the first genetic map for date palm and identify the putative date palm sex chromosome. We placed ~4000 markers on the map using nearly 1200 framework markers spanning a total of 1293 cM. We have integrated the genetic map, derived from the Khalas cultivar, with the draft genome and placed up to 19% of the draft genome sequence scaffolds onto linkage groups for the first time. This analysis revealed approximately ~1.9 cM/Mb on the map. Comparison of the date palm linkage groups revealed significant long-range synteny to oil palm. Analysis of the date palm sex-determination region suggests it is telomeric on linkage group 12 and recombination is not suppressed in the full chromosome.

Conclusions

Based on a modified gentoyping-by-sequencing approach we have overcome challenges due to lack of genetic resources and provide the first genetic map for date palm. Combined with the recent draft genome sequence of the same cultivar, this resource offers a critical new tool for date palm biotechnology, palm comparative genomics and a better understanding of sex chromosome development in the palms.

Estimation of genealogical coancestry in plant species using a pedigree reconstruction algorithm and application to an oil palm breeding population

Explicit pedigree reconstruction by simulated annealing gave reliable estimates of genealogical coancestry in plant species, especially when selfing rate was lower than 0.6, using a realistic number of markers. Genealogical coancestry information is crucial in plant breeding to estimate genetic parameters and breeding values. The approach of Fernández and Toro (Mol Ecol 15:1657-1667, 2006) to estimate genealogical coancestries from molecular data through pedigree reconstruction was limited to species with separate sexes. In this study it was extended to plants, allowing hermaphroditism and monoecy, with possible selfing. Moreover, some improvements were made to take previous knowledge on the population demographic history into account. The new method was validated using simulated and real datasets. Simulations showed that accuracy of estimates was high with 30 microsatellites, with the best results obtained for selfing rates below 0.6. In these conditions, the root mean square error (RMSE) between the true and estimated genealogical coancestry was small (<0.07), although the number of ancestors was overestimated and the selfing rate could be biased. Simulations also showed that linkage disequilibrium between markers and departure from the Hardy-Weinberg equilibrium in the founder population did not affect the efficiency of the method. Real oil palm data confirmed the simulation results, with a high correlation between the true and estimated genealogical coancestry (>0.9) and a low RMSE (<0.08) using 38 markers. The method was applied to the Deli oil palm population for which pedigree data were scarce. The estimated genealogical coancestries were highly correlated (>0.9) with the molecular coancestries using 100 markers. Reconstructed pedigrees were used to estimate effective population sizes. In conclusion, this method gave reliable genealogical coancestry estimates. The strategy was implemented in the software MOLCOANC 3.0.

Analyses of hypomethylated oil palm gene space

Demand for palm oil has been increasing by an average of ∼8% the past decade and currently accounts for about 59% of the world's vegetable oil market. This drives the need to increase palm oil production. Nevertheless, due to the increasing need for sustainable production, it is imperative to increase productivity rather than the area cultivated. Studies on the oil palm genome are essential to help identify genes or markers that are associated with important processes or traits, such as flowering, yield and disease resistance. To achieve this, 294,115 and 150,744 sequences from the hypomethylated or gene-rich regions of Elaeis guineensis and E. oleifera genome were sequenced and assembled into contigs. An additional 16,427 shot-gun sequences and 176 bacterial artificial chromosomes (BAC) were also generated to check the quality of libraries constructed. Comparison of these sequences revealed that although the methylation-filtered libraries were sequenced at low coverage, they still tagged at least 66% of the RefSeq supported genes in the BAC and had a filtration power of at least 2.0. A total 33,752 microsatellites and 40,820 high-quality single nucleotide polymorphism (SNP) markers were identified. These represent the most comprehensive collection of microsatellites and SNPs to date and would be an important resource for genetic mapping and association studies. The gene models predicted from the assembled contigs were mined for genes of interest, and 242, 65 and 14 oil palm transcription factors, resistance genes and miRNAs were identified respectively. Examples of the transcriptional factors tagged include those associated with floral development and tissue culture, such as homeodomain proteins, MADS, Squamosa and Apetala2. The E. guineensis and E. oleifera hypomethylated sequences provide an important resource to understand the molecular mechanisms associated with important agronomic traits in oil palm.

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

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.

Lipase in the mesocarp of oil palm fruits reduces the quality of oil and leads to decreased yield. In this study, the authors identify the gene encoding the mesocarp lipase, raising the possibility of breeding palm genotypes that produce better quality oil.

The oil palm SHELL gene controls oil yield and encodes a homologue of SEEDSTICK

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 pisifera^1–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 Asia⁵. 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 conservation⁶.

Identification of QTLs associated with callogenesis and embryogenesis in oil palm using genetic linkage maps improved with SSR markers

Clonal reproduction of oil palm by means of tissue culture is a very inefficient process. Tissue culturability is known to be genotype dependent with some genotypes being more amenable to tissue culture than others. In this study, genetic linkage maps enriched with simple sequence repeat (SSR) markers were developed for dura (ENL48) and pisifera (ML161), the two fruit forms of oil palm, Elaeis guineensis. The SSR markers were mapped onto earlier reported parental maps based on amplified fragment length polymorphism (AFLP) and restriction fragment length polymorphism (RFLP) markers. The new linkage map of ENL48 contains 148 markers (33 AFLPs, 38 RFLPs and 77 SSRs) in 23 linkage groups (LGs), covering a total map length of 798.0 cM. The ML161 map contains 240 markers (50 AFLPs, 71 RFLPs and 119 SSRs) in 24 LGs covering a total of 1,328.1 cM. Using the improved maps, two quantitative trait loci (QTLs) associated with tissue culturability were identified each for callusing rate and embryogenesis rate. A QTL for callogenesis was identified in LGD4b of ENL48 and explained 17.5% of the phenotypic variation. For embryogenesis rate, a QTL was detected on LGP16b in ML161 and explained 20.1% of the variation. This study is the first attempt to identify QTL associated with tissue culture amenity in oil palm which is an important step towards understanding the molecular processes underlying clonal regeneration of oil palm.

Elaeis oleifera genomic-SSR markers: exploitation in oil palm germplasm diversity and cross-amplification in arecaceae

Species-specific simple sequence repeat (SSR) markers are favored for genetic studies and marker-assisted selection (MAS) breeding for oil palm genetic improvement. This report characterizes 20 SSR markers from an Elaeis oleifera genomic library (gSSR). Characterization of the repeat type in 2000 sequences revealed a high percentage of di-nucleotides (63.6%), followed by tri-nucleotides (24.2%). Primer pairs were successfully designed for 394 of the E. oleifera gSSRs. Subsequent analysis showed the ability of the 20 selected E. oleifera gSSR markers to reveal genetic diversity in the genus Elaeis. The average Polymorphism Information Content (PIC) value for the SSRs was 0.402, with the tri-repeats showing the highest average PIC (0.626). Low values of observed heterozygosity (H_o) (0.164) and highly positive fixation indices (F_is) in the E. oleifera germplasm collection, compared to the E. guineensis, indicated an excess of homozygosity in E. oleifera. The transferability of the markers to closely related palms, Elaeis guineensis, Cocos nucifera and ornamental palms is also reported. Sequencing the amplicons of three selected E. oleifera gSSRs across both species and palm taxa revealed variations in the repeat-units. The study showed the potential of E. oleifera gSSR markers to reveal genetic diversity in the genus Elaeis. The markers are also a valuable genetic resource for studying E. oleifera and other genus in the Arecaceae family.

SSR markers in transcripts of genes linked to post-transcriptional and transcriptional regulatory functions during vegetative and reproductive development of Elaeis guineensis

BACKGROUND

The oil palm (Elaeis guineensis Jacq.) is a perennial monocotyledonous tropical crop species that is now the world's number one source of edible vegetable oil, and the richest dietary source of provitamin A. While new elite genotypes from traditional breeding programs provide steady yield increases, the long selection cycle (10-12 years) and the large areas required to cultivate oil palm make genetic improvement slow and labor intensive. Molecular breeding programs have the potential to make significant impacts on the rate of genetic improvement but the limited molecular resources, in particular the lack of molecular markers for agronomic traits of interest, restrict the application of molecular breeding schemes for oil palm.

RESULTS

In the current study, 6,103 non-redundant ESTs derived from cDNA libraries of developing vegetative and reproductive tissues were annotated and searched for simple sequence repeats (SSRs). Primer pairs from sequences flanking 289 EST-SSRs were tested to detect polymorphisms in elite breeding parents and their crosses. 230 of these amplified PCR products, 88 of which were polymorphic within the breeding material tested. A detailed analysis and annotation of the EST-SSRs revealed the locations of the polymorphisms within the transcripts, and that the main functional category was related to transcription and post-transcriptional regulation. Indeed, SSR polymorphisms were found in sequences encoding AP2-like, bZIP, zinc finger, MADS-box, and NAC-like transcription factors in addition to other transcriptional regulatory proteins and several RNA interacting proteins.

CONCLUSIONS

The identification of new EST-SSRs that detect polymorphisms in elite breeding material provides tools for molecular breeding strategies. The identification of SSRs within transcripts, in particular those that encode proteins involved in transcriptional and post-transcriptional regulation, will allow insight into the functional roles of these proteins by studying the phenotypic traits that cosegregate with these markers. Finally, the oil palm EST-SSRs derived from vegetative and reproductive development will be useful for studies on the evolution of the functional diversity within the palm family.

Genetic linkage map of a high yielding FELDA deli×yangambi oil palm cross

Enroute to mapping QTLs for yield components in oil palm, we constructed the linkage map of a FELDA high yielding oil palm (Elaeis guineensis), hybrid cross. The parents of the mapping population are a Deli dura and a pisifera of Yangambi origin. The cross out-yielded the average by 8–21% in four trials all of which yielded comparably to the best current commercial planting materials. The higher yield derived from a higher fruit oil content. SSR markers in the public domain - from CIRAD and MPOB, as well as some developed in FELDA - were used for the mapping, augmented by locally-designed AFLP markers. The female parent linkage map comprised 317 marker loci and the male parent map 331 loci, both in 16 linkage groups each. The number of markers per group ranged from 8–47 in the former and 12–40 in the latter. The integrated map was 2,247.5 cM long and included 479 markers and 168 anchor points. The number of markers per linkage group was 15–57, the average being 29, and the average map density 4.7 cM. The linkage groups ranged in length from 77.5 cM to 223.7 cM, with an average of 137 cM. The map is currently being validated against a closely related population and also being expanded to include yield related QTLs.

Development of SSR and gene-targeted markers for construction of a framework linkage map of Catharanthus roseus

BACKGROUND AND AIMS

Catharanthus roseus is a plant of great medicinal importance, yet inadequate knowledge of its genome structure and the unavailability of genomic resources have been major impediments in the development of improved varieties. The aims of this study were to develop co-dominant sequence-tagged microsatellite sites (STMS) and gene-targeted markers (GTMs) and utilize them for the construction of a framework intraspecific linkage map of C. roseus.

METHODS

For simple sequence repeat (SSR) isolation, a genomic library enriched for (GA)(n) repeats was constructed from C. roseus 'Nirmal' (CrN1). In addition, GTMs were also designed from 12 genes of the TIA (terpenoid indole alkaloid) pathway - the medicinally most significant pathway in C. roseus. An F(2) mapping population was also generated by crossing two diverse accessions of C. roseus CrN1 (Nirmal)×CrN82 (Kew).

KEY RESULTS

A new set of 314 STMS markers and 64 GTMs were developed in this study. A segregating F(2) mapping population consisting of 111 F(2) individuals was generated. For generating the linkage map, a set of 423 co-dominant markers (378 newly developed and 45 published earlier) were screened for polymorphism between the parental genotypes, of which 134 were identified to be polymorphic. A total of 114 markers were mapped on eight linkage groups that spanned a 632·7 cM region of the genome with an average marker distance of 5·55 cM. Further, the mechanism of hypervariability at the gene-targeted loci was investigated at the sequence level.

CONCLUSIONS

For the first time, a large array of STMS markers and GTMs was generated in the model medicinal plant C. roseus. Moreover, the first microsatellite marker-based linkage map was described in this study. Together, these will serve as a foundation for future genomics studies related to quantitative trait loci analysis and molecular breeding in C. roseus.

Genetic variability of oil palm parental genotypes and performance of its' progenies as revealed by molecular markers and quantitative traits

Studies were conducted to assess the genetic relationships between the parental palms (dura and pisifera) and performance of their progenies based on nine microsatellite markers and 29 quantitative traits. Correlation analyses between genetic distances and hybrids performance were estimated. The coefficients of correlation values of genetic distances with hybrid performance were non-significant, except for mean nut weight and leaf number. However, the correlation coefficient of genetic distances with these characters was low to be used as predicted value. These results indicated that genetic distances based on the microsatellite markers may not be useful for predicting hybrid performance. The genetic distance analysis using UPGMA clustering system generated 5 genetic clusters with coefficient of 1.26 based on quantitative traits of progenies. The genotypes, DP16, DP14, DP4, DP13, DP12, DP15, DP8, DP1 and DP2 belonging to distant clusters and greater genetic distances could be selected for further breeding programs.