Improving the accuracy of genomic predictions in an outcrossing species with hybrid cultivars between heterozygote parents: a case study of oil palm (Elaeis guineensis Jacq.).

Genomic selection (GS) is a method of marker-assisted selection revolutionizing crop improvement, but it can still be optimized. For hybrid breeding between heterozygote parents of different populations or species, specific aspects can be considered to increase GS accuracy: (1) training population genotyping, i.e., only genotyping the hybrid parents or also a sample of hybrid individuals, and (2) marker effects modeling, i.e., using population-specific effects of single nucleotide polymorphism alleles model (PSAM) or across-population SNP genotype model (ASGM). Here, this was investigated empirically for the prediction of the performances of oil palm hybrids for yield traits. The GS model was trained on 352 hybrid crosses and validated on 213 independent hybrid crosses. The training and validation hybrid parents and 399 training hybrid individuals were genotyping by sequencing. Despite the small proportion of hybrid individuals genotyped and low parental heterozygosity, GS prediction accuracy increased on average by 5% (range 1.4-31.3%, depending on trait and model) when training was done using genomic data on hybrids and parents compared with only parental genomic data. With ASGM, GS prediction accuracy increased on average by 3% (- 10.2 to 40%, depending on trait and genotyping strategy) compared with PSAM. We conclude that the best GS strategy for oil palm is to aggregate genomic data of parents and hybrid individuals and to ignore the parental origin of marker alleles (ASGM). To gain a better insight into these results, future studies should examine the respective effect of capturing genetic variability within crosses and taking segregation distortion into account when genotyping hybrid individuals, and investigate the factors controlling the relative performances of ASGM and PSAM in hybrid crops.

Genomic preselection with genotyping-by-sequencing increases performance of commercial oil palm hybrid crosses.

BACKGROUND: There is great potential for the genetic improvement of oil palm yield. Traditional progeny tests allow accurate selection but limit the number of individuals evaluated. Genomic selection (GS) could overcome this constraint. We estimated the accuracy of GS prediction of seven oil yield components using A × B hybrid progeny tests with almost 500 crosses for training and 200 crosses for independent validation. Genotyping-by-sequencing (GBS) yielded +5000 single nucleotide polymorphisms (SNPs) on the parents of the crosses. The genomic best linear unbiased prediction method gave genomic predictions using the SNPs of the training and validation sets and the phenotypes of the training crosses. The practical impact was illustrated by quantifying the additional bunch production of the crosses selected in the validation experiment if genomic preselection had been applied in the parental populations before progeny tests. RESULTS: We found that prediction accuracies for cross values plateaued at 500 to 2000 SNPs, with high (0.73) or low (0.28) values depending on traits. Similar results were obtained when parental breeding values were predicted. GS was able to capture genetic differences within parental families, requiring at least 2000 SNPs with less than 5% missing data, imputed using pedigrees. Genomic preselection could have increased the selected hybrids bunch production by more than 10%. CONCLUSIONS: Finally, preselection for yield components using GBS is the first possible application of GS in oil palm. This will increase selection intensity, thus improving the performance of commercial hybrids. Further research is required to increase the benefits from GS, which should revolutionize oil palm breeding.

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.

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

KEY MESSAGE: 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 properties of key oil palm (Elaeis guineensis Jacq.) breeding populations.

A good knowledge of the genome properties of the populations makes it possible to optimize breeding methods, in particular genomic selection (GS). In oil palm (Elaeis guineensis Jacq), the world's main source of vegetable oil, this would provide insight into the promising GS results obtained so far. The present study considered two complex breeding populations, Deli and La Mé, with 943 individuals and 7324 single-nucleotide polymorphisms (SNPs) from genotyping-by-sequencing. Linkage disequilibrium (LD), haplotype sharing, effective size (Ne), and fixation index (Fst) were investigated. A genetic linkage map spanning 1778.52 cM and with a recombination rate of 2.85 cM/Mbp was constructed. The LD at r2=0.3, considered the minimum to get reliable GS results, spanned over 1.05 cM/0.22 Mbp in Deli and 0.9 cM/0.21 Mbp in La Mé. The significant degree of differentiation existing between Deli and La Mé was confirmed by the high Fst value (0.53), the pattern of correlation of SNP heterozygosity and allele frequency among populations, and the decrease of persistence of LD and of haplotype sharing among populations with increasing SNP distance. However, the level of resemblance between the two populations over short genomic distances (correlation of r values between populations >0.6 for SNPs separated by <0.5 cM/1 kbp and percentage of common haplotypes >40% for haplotypes <3600 bp/0.20 cM) likely explains the superiority of GS models ignoring the parental origin of marker alleles over models taking this information into account. The two populations had low Ne (<5). Population-specific genetic maps and reference genomes are recommended for future studies.

Genomic predictions improve clonal selection in oil palm (Elaeis guineensis Jacq.) hybrids.

The prediction of clonal genetic value for yield is challenging in oil palm (Elaeis guineensis Jacq.). Currently, clonal selection involves two stages of phenotypic selection (PS): ortet preselection on traits with sufficient heritability among a small number of individuals in the best crosses in progeny tests, and final selection on performance in clonal trials. The present study evaluated the efficiency of genomic selection (GS) for clonal selection. The training set comprised almost 300 Deli × La Mé crosses phenotyped for eight palm oil yield components and the validation set 42 Deli × La Mé ortets. Genotyping-by-sequencing (GBS) revealed 15,054 single nucleotide polymorphisms (SNP). The effects of the SNP dataset (density and percentage of missing data) and two GS modeling approaches, ignoring (ASGM) and considering (PSAM) the parental origin of alleles, were assessed. The results showed prediction accuracies ranging from 0.08 to 0.70 for ortet candidates without data records, depending on trait, SNP dataset and modeling. ASGM was better (on average slightly more accurate, less sensitive to SNP dataset and simpler), although PSAM appeared interesting for a few traits. With ASGM, the number of SNPs had to reach 7,000, while the percentage of missing data per SNP was of secondary importance, and GS prediction accuracies were higher than those of PS for most of the traits. Finally, this makes possible two practical applications of GS, that will increase genetic progress by improving ortet preselection before clonal trials: (1) preselection at the mature stage on all yield components jointly using ortet genotypes and phenotypes, and (2) genomic preselection on more yield components than PS, among a large population of the best possible crosses at nursery stage.

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.