Performance of multi-trait genomic selection for Eucalyptus robusta breeding program

In forest tree genetic improvement, multi-trait genomic selection (GS) may have advantages in improving the accuracy of the genotype estimation and shortening selection cycles. For the breeding of Eucalyptus robusta, one of the most exotic planted species in Madagascar, volume at 49 months (V49), total lignin (TL), and holo-cellulose (Holo) were considered. For GS, 2919 single nucleotide polymorphisms (SNP) were used with the genomic best linear unbiased predictor (GBLUP) method, which was as efficient as the reproducing kernel Hilbert space (RKHS) and elastic net methods (EN), but more adapted to multi-trait modeling. The efficiency of individual I model, including the genomic data, was much higher than the provenance effect P model. For example, with V49, mean goodness-of-fit was: rI_Full = 0.79, rP_Full = 0.37 for I and P, respectively. The prediction accuracies using the cross-validation procedure were lower for V49: rI = 0.29 rP = 0.28. The genetic gains resulting from the indexes associating (V49, TL) and (V49, Holo) were higher using I than for the P model; for V49, the relative genetic gain was 37 and 20%, respectively, with 5% of selection intensity. The single-trait approach was as efficient as the multi-trait approach given the weak correlations between V49 and TL or Holo. The I model also brings greater diversity: for V49 the number of provenances represented in a selected population was two and three with the P model, and 6 and 16 with the I model.