Simultaneous improvement of grain yield and protein content in durum wheat by different phenotypic indices and genomic selection

The additive relationship matrix plays an important role in mixed model prediction of breeding values. For genotype matrix X (loci in columns), the product XX′ is widely used as a realized relationship matrix, but the scaling of this matrix is ambiguous. Our first objective was to derive a proper scaling such that the mean diagonal element equals 1+f, where f is the inbreeding coefficient of the current population. The result is a formula involving the covariance matrix for sampling genomic loci, which must be estimated with markers. Our second objective was to investigate whether shrinkage estimation of this covariance matrix can improve the accuracy of breeding value (GEBV) predictions with low-density markers. Using an analytical formula for shrinkage intensity that is optimal with respect to mean-squared error, simulations revealed that shrinkage can significantly increase GEBV accuracy in unstructured populations, but only for phenotyped lines; there was no benefit for unphenotyped lines. The accuracy gain from shrinkage increased with heritability, but at high heritability (> 0.6) this benefit was irrelevant because phenotypic accuracy was comparable. These trends were confirmed in a commercial pig population with progeny-test-estimated breeding values. For an anonymous trait where phenotypic accuracy was 0.58, shrinkage increased the average GEBV accuracy from 0.56 to 0.62 (SE < 0.00) when using random sets of 384 markers from a 60K array. We conclude that when moderate-accuracy phenotypes and low-density markers are available for the candidates of genomic selection, shrinkage estimation of the relationship matrix can improve genetic gain.

Grain yield is a major goal for the improvement of durum wheat, particularly in drought-prone areas. In this study, the genetic basis of grain yield (GY), heading date (HD), and plant height (PH) was investigated in a durum wheat population of 249 recombinant inbred lines evaluated in 16 environments (10 rainfed and 6 irrigated) characterized by a broad range of water availability and GY (from 5.6 to 58.8 q ha(-1)). Among the 16 quantitative trait loci (QTL) that affected GY, two major QTL on chromosomes 2BL and 3BS showed significant effects in 8 and 7 environments, with R2 values of 21.5 and 13.8% (mean data of all 16 environments), respectively. In both cases, extensive overlap was observed between the LOD profiles of GY and PH, but not with those for HD. QTL specific for PH were identified on chromosomes 1BS, 3AL, and 7AS. Additionally, three major QTL for HD on chromosomes 2AS, 2BL, and 7BS showed limited or no effects on GY. For both PH and GY, notable epistasis between the chromosome 2BL and 3BS QTL was detected across several environments.