Comparison of wholegenome prediction models for traits with contrasting genetic architecture in a diversity panel of maize inbred lines
 Christian Riedelsheimer^{1},
 Frank Technow^{1} and
 Albrecht E Melchinger^{1}Email author
DOI: 10.1186/1471216413452
© Riedelsheimer et al.; licensee BioMed Central Ltd. 2012
Received: 27 February 2012
Accepted: 14 August 2012
Published: 4 September 2012
Abstract
Background
There is increasing empirical evidence that wholegenome prediction (WGP) is a powerful tool for predicting line and hybrid performance in maize. However, there is a lack of knowledge about the sensitivity of WGP models towards the genetic architecture of the trait. Whereas previous studies exclusively focused on highly polygenic traits, important agronomic traits such as disease resistances, nutrifunctional or climate adaptational traits have a genetic architecture which is either much less complex or unknown. For such cases, information about model robustness and guidelines for model selection are lacking. Here, we compared five WGP models with different assumptions about the distribution of the underlying genetic effects. As contrasting model traits, we chose three highly polygenic agronomic traits and three metabolites each with a major QTL explaining 22 to 30% of the genetic variance in a panel of 289 diverse maize inbred lines genotyped with 56,110 SNPs.
Results
We found the five WGP models to be remarkable robust towards trait architecture with the largest differences in prediction accuracies ranging between 0.05 and 0.14 for the same trait, most likely as the result of the high level of linkage disequilibrium prevailing in elite maize germplasm. Whereas RRBLUP performed best for the agronomic traits, it was inferior to LASSO or elastic net for the three metabolites. We found the approach of genome partitioning of genetic variance, first applied in human genetics, as useful in guiding the breeder which model to choose, if prior knowledge of the trait architecture is lacking.
Conclusions
Our results suggest that in diverse germplasm of elite maize inbred lines with a high level of LD, WGP models differ only slightly in their accuracies, irrespective of the number and effects of QTL found in previous linkage or association mapping studies. However, small gains in prediction accuracies can be achieved if the WGP model is selected according to the genetic architecture of the trait. If the trait architecture is unknown e.g. for novel traits which only recently received attention in breeding, we suggest to inspect the distribution of the genetic variance explained by each chromosome for guiding model selection in WGP.
Keywords
Genomic selection Wholegenome prediction Genetic architecture Complex traits Zea maysBackground
Wholegenome prediction (WGP) is expected to reshape plant breeding fundametally in the near future [1–3]. Whereas the approach has been initially proposed [4] and rapidly implemented in animal breeding [5], recent empirical studies demonstrated also its potential in hybrid maize breeding [6–8]. Recently, we showed that WGP allows a reliable screening of large germplasm collections of diverse maize inbred lines for their potential to create superior hybrids [9]. However, these studies exclusively focused on predicting highly polygenic traits such as grain yield or biomass accumulation with genetic architectures close to the infinitesimal genetic model [10].
In maize, several economically important traits are genetically less complex with few quantitative trait loci (QTL) explaining a large proportion of the genetic variance. Examples include pest and disease resistances or nutrifunctional compounds such as bioavailable minerals [11] or βcarotene [12]. In addition, disease resistances are often found to be controlled by a combination of racespecific resistance loci with large effects involved in pathogen recognition, and a large number of loci with small effects involved in basal resistance. Such a mixed QTL effect distribution can be found in maize e.g. for rust [13], Giberella ear rot [14, 15] or to a lesser extent for Northern corn leaf blight [16, 17].
For such traits, the assumption of normally distributed SNP effects underlying ridge regression, the most commonly applied WGP model, is severely violated. Heslot et al.[18] found for polygenic traits in several plant species only minor differences between ridge regression and models with different assumptions of the underlying distribution of SNP effects. However, these differences are expected to be much larger for traits controlled by only a few QTL. Recently, Clark et al.[19] simulated this situation under the assumption of the historical population structure of Holstein cattle. They found that under the assumption of either few common or few rare quantitative trait loci, a Bayesian variable selection model (BayesB) outperforms ridge regression by far. For Holstein cattle, Hayes et al.[20] found also the BayesA model to be superior to ridge regression in the case of coat color or milkfat percentage.
Cattle differs greatly in its population structure and LD level from elite maize germplasm, which has faced severe genetic bottlenecks during domestication and the creation of genetically distinct heterotic pools to maximize exploitation of heterosis in hybrid breeding [21, 22]. Hence, results from cattle might not be directly transferable to maize, for which little is known about WGP for traits with a simpler genetic architecture. Moreover, the genetic architecture of a trait is often unclear in crops. Especially if the trait has not yet been extensively dissected by linkage or association mapping, which might be the case for traits which gained only recently in importance such as nutritional properties, nutrient acquisition traits or traits related to climate change adaptation.

To what extent do distinct WGP models differ in their prediction accuracies for a diversity panel of maize inbred lines if the genetic architecture of the trait changes dramatically?

Are there guidelines for plants concerning the choice of the most promising WGP model?
Methods
Genetic material
The genetic material consisted of 289 maize inbred lines which were previously described in great detail [9, 23–25]. The population constituted a global sample of elite breeding material from worldwide sources with a focus on North America and Europe and encompassed 285 lines from the Dent heterotic pool (StiffStalk and nonStiffStalk) and 4 from the European Flint pool, which served as check genotypes.
Genotyping
The population was genotyped with the Illumina SNP chip MaizeSNP50 containing 56,110 SNPs [26]. Quality control preprocessing of SNPs was performed by eliminating SNPs that did not match the following criteria: (i) less than 10% missing values, (ii) minor allele frequency of greater than 2.5%, (iii) no more than three heterozygous genotypes, and (iv) unique allele assignment for the 22 replicated checks of genotype B73. A total of 38,019 SNPs remained and were used for further analysis. Linkage disequilibrium (LD) declined to r^{2}=0.1 at approximately 500 kb with a mean LD between adjacent SNPs of 0.34 [9].
Field trials
The population was phenotyped in six environments (three agroecologically diverse locations in the years 2008 and 2009) in Germany [25]. Briefly, the population was split into three maturity groups based on prior knowledge of their flowering time. In the trials of each of the three maturity groups, 100 genotypes, including five common check genotypes, were randomized in a 20 × 5 αlattice design with two replications and were planted in 2row plots. Plots were thinned to a final plant density of 100,000 plants/ha. The common check genotypes were used to adjust for potential differences in the soil fertility among trials in each environment.
Metabolites
Phenotypic correlations among traits
  Plant height  Lignin content  Dopamine  Ribitol  719700204 

Dry matter yield  0.62  0.32  0.28  0.02  0.24 
Plant height    0.50  0.17  0.07  0.12 
Lignin content      0.20  0.11  0.08 
Dopamine        0.10  0.08 
Ribitol          0.03 
Agronomic traits
Dry matter yield of wholeplant biomass (t/ha) and plant height (m) were measured per field plot of the inbred lines. Lignin content was measured as acid detergent lignin (ADL) in the harvested plant material of the inbred lines using calibrated nearinfrared spectroscopy (NIRS). The NIRS calibration model was built using phenotypic data from 20 inbred lines, 32 testcrosses and 3 hybrids grown in the same environments as the population of inbred lines analyzed in this study [24]. Heritability estimates and genotypic means were obtained using a onestep linear mixed model analysis as described previously [25]. Using a 1% Bonferroni corrected significance threshold, we could not find any significant SNPtrait association signal using the the same GWA model as for metabolites. Since population size, marker density, and heritabilites were sufficiently high for detecting QTL with large effects, the absence of any significant traitSNP associations suggest a highly polygenic genetic architecture for the agronomic traits with no major QTL.
Genome partitioning of the genetic variance
To further characterize the genetic architectures of the investigated traits irrespective of the significance thresholds for SNPtrait associations, we compared how the ten chromosomes contributed to the total genetic variance. Later on, we will use these results as a guideline for model selection based on trait architecture.
where y is a vector with n trait values, 1 is a vector of 1’s, Q is a matrix of size n×10 containing the first 10 principal components calculated from SNP data with β containing the corresponding regression coefficients, S is an incidence matrix allocating components of y to components of g_{ c }, which is a vector of length n with random genotypic effects attributable to chromosome c with $\left(\right)close="">{\mathbf{g}}_{c}\sim N(0,{\mathbf{G}}_{c}{\sigma}_{\mathrm{gc}}^{2})$ and $\left(\right)close="">{\mathbf{G}}_{c}={\mathbf{Z}}_{c}{\mathbf{Z}}_{c}^{T}/{p}_{c}$ where Z_{ c } is a matrix of size n×p_{ c }with standardized levels of SNP alleles on chromosome c. Vector e contains normally distributed residuals with $\left(\right)close="">\mathbf{e}\sim N(0,\mathbf{I}{\sigma}_{e}^{2})$. The genetic variance contributed by chromosome c was then estimated as $\left(\right)close="">{\sigma}_{\mathrm{gc}}^{2}/\left(\sum _{c=1}^{10}\right({\sigma}_{\mathrm{gc}}^{2})+{\sigma}_{e}^{2})$.
Variance components were estimated by restricted maximum likelihood (REML) using ASRemlR 3 [29]. Since matrices G_{ c }were often found to be singular, we used the algorithm of Higham [30] implemented in the function nearPD of the Rpackage Matrix [31], to approximate the nearest positive definite matrices.
WGP models
We investigated five WGP models that have been recently advocated in the literature for this purpose [4, 32–34].
where y is a vector with n trait values, μ is the overall mean, 1 is a vector with 1’s, Z is the n×p matrix of standardized values of SNP alleles, u is a vector with SNP effects, and e is a vector of residuals with $\left(\right)close="">\mathbf{e}\sim N(0,\mathbf{I}{\sigma}_{e}^{2})$. Depending on the trait, a combination of genotypic and phenotypic information was available for 276 to 280 genotypes which were used for WGP.
RRBLUP
The Lagrangian multiplier λ_{RR}is a regularization parameter which controls the amount of shrinkage. It can be estimated as $\left(\right)close="">{\lambda}_{\mathrm{RR}}={\sigma}_{e}^{2}/{\sigma}_{u}^{2}$ by regarding u as random effects with $\left(\right)close="">\mathbf{u}\sim N(0,\mathbf{I}{\sigma}_{u}^{2})$ with ${\sigma}_{u}^{2}$ being the SNP effect variance estimated by REML. In this setting, $\left(\right)close="">{\widehat{\mathbf{u}}}_{\mathrm{RR}}$ is equivalent to the best linear unbiased predictor (BLUP) of u[35, 36].
with g being a vector of random genotype effects with var(g) = $\left(\right)close="">\mathbf{G}{\sigma}_{g}^{2}$ and wholegenome relationship matrix G=ZZ^{ T }/p. The solution vector of SNP effects can then be obtained as $\left(\right)close="">{\widehat{\mathbf{u}}}_{\mathrm{RR}}={\mathbf{Z}}^{T}{\mathbf{G}}^{1}\widehat{\mathbf{g}}$[37]. Here, G is an innerproduct kernel which allows to perform all computations in the space of n genotypes instead of p SNPs, a shortcut which is well established in the field of kernelbased machine learning [38].
LASSO
which bounds the Manhattan (L_{1}) norm of u to a constraint: $\left(\right)close="">\left\right\mathbf{u}{}_{1}=\sum _{i=1}^{p}{u}_{i}{c}_{\mathrm{L}}$. The LASSO penalty is a diamond shaped constraint which allows not only to shrink coefficients towards zero but to set some coefficients to exactly zero (Figure 1). Unlike RR, LASSO, cannot be ’kernelized’, i.e., it is not possible to transform the LASSO estimator into an equivalent kernel regression problem in the space of n genotypes [38]. Hence, LASSO regression has to be carried out with the full set of SNPs. Here, we used the R package glmnet, a fast implementation using cyclic coordinate descent to compute the complete LASSO path solution [42].
Elastic net
and is a weighted mixture between the RR penalty (α=0) and the LASSO penalty (α=1) [33]. While the RR penalty encourages highly correlated variables to be averaged, the LASSO penalty encourages a sparse solution [38]. We again used the implementation in glmnet and performed a grid search to find the combination of α and λ_{EN}, which yielded the lowest mean squared error in the training population.
Reproducing kernel Hilbert space (RKHS) regression
The theory of RKHS regression is rooted in the field of kernelbased machinelearning [38] and has recently been advocated for wholegenome prediction [34]. The approach uses equation 4 but replaces the innerproduct matrix G with a kernel matrix K. The motivation behind RKHS regression lies in the ability to effectively perform nonlinear regression in a higherdimensional feature space so it might capture nonadditive genetic effects, if present. Here, we used a Gaussian kernel on genetic distances with ${K}_{\mathrm{ij}}=exp({\mathrm{GD}}_{\mathrm{ij}}/{\theta}^{2})$, where GD_{ ij } is the modified Rogers’ genetic distance (Euclidean distance scaled to fall between 0 and 1) between genotype i and j, and θ is a smoothing parameter which controls the rate of decay of K_{ ij } with increasing genetic distance. The optimum value for θ was chosen from a sequence from 0.1 to 100 at which the maximum likelihood was obtained.
BayesB
Priors used for BayesB
Parameter  Prior 

u _{ i }  $N(0,{\sigma}_{\mathrm{ui}}^{2})$ 
$\left(\right)close="">{\sigma}_{\mathrm{ui}}^{2}{v}_{u},{S}_{u}^{2}$  $\left\{\begin{array}{l}0\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{2.36043pt}{0ex}}\phantom{\rule{2.36043pt}{0ex}}\text{with probability}\phantom{\rule{2.36043pt}{0ex}}{\Pi}_{u},\\ {\chi}^{2}({v}_{u},{S}_{u}^{2})\phantom{\rule{2.36043pt}{0ex}}\phantom{\rule{2.36043pt}{0ex}}\text{with probability}\phantom{\rule{2.36043pt}{0ex}}(1{\Pi}_{u})\end{array}\right.$ 
v _{ u }  Gamma(k=5,θ=2) 
$\left(\right)close="">{S}_{u}^{2}$  Gamma(k=0.1,θ=10) 
Π _{ u }  Beta(α=7,β=3) 
$\left(\right)close="">{\sigma}_{e}^{2}$  $\left(\right)close="">{\chi}^{2}({v}_{e},{S}_{e}^{2}=\widehat{{\sigma}_{e}^{2}}({v}_{e}2)/{v}_{e})$ 
v_{ e }=4.001, $\left(\right)close="">\widehat{{\sigma}_{e}^{2}}$ estimated with REML 
Validation
A fivefold crossvalidation scheme was applied and repeated 20 times. In each repetition, the dataset was divided into 5 disjoint subsets of genotypes whereas one subset served as the validation set and the other four subsets served as the training population to estimate the model parameters for predicting the leftout genotypes in the validation set. In each of the five rounds, the Pearson correlation between the observed and predicted phenotypic values was calculated. The procedure was repeated twenty times to yield 100 crossvalidation runs. The predictive ability was then calculated as the Pearson correlation ($\left(\right)close="">{r}_{(\mathbf{y},\widehat{\mathbf{y}})}$) between the observed (y) and predicted ($\widehat{\mathbf{y}}$) phenotypic values. The ’prediction accuracy’ estimates the correlation (${r}_{(\mathbf{g},\widehat{\mathbf{g}})}$) between the predicted ($\widehat{\mathbf{g}}$) and unobserved true genetic values (g) and was calculated by ${r}_{(\mathbf{g},\widehat{\mathbf{g}})}={r}_{(\mathbf{y},\widehat{\mathbf{y}})}/h$ where h is the square root of the heritability on a linemean basis for the agronomic traits. For metabolites, the square root of the estimated repeatability was used.
Results
For the agronomic traits, the total genetic variance was largely uniformly distributed over all chromosomes. Using the same GWA model as for the metabolites [23], we found that for dry matter yield and plant height, the chromosomes which captured the largest portion of genetic variance contain the strongest GWA signals. However, in no instance was the 1% Bonferroni corrected significance threshold surpassed (dry matter yield: chr. 6, P=4.06×10^{−6}, position 139,284,469, explained genetic variance 8.8%; plant height: chr. 3, P=8.6×10^{−6}, position 163,617,228, explained genetic variance 8.0%).
When excluding chromosomes containing either these two association signals or major mQTL, we observed a tendency that longer chromosomes captured more genetic variance than shorter ones (Figure 2B). This trend was significant (P<0.10) for lignin content (r=0.88,P=7.0×10^{−4}) and dry matter yield (r=0.60,P=0.09). The grey area in Figure 2B was therefore regarded as the range in chromosomal genetic variance explainable by the length of the chromosome. On the other side, the genetic variance contributed by the leftout chromosomes was highly correlated (r=0.98,P=0.003) with the explained genetic variance of the individual SNPs found by GWA mapping (Figure 2C).
Prediction accuracies $\mathbf{\left(}{\mathbf{r}}_{\mathbf{(}\mathbf{g}\mathbf{,}\widehat{\mathbf{g}}\mathbf{)}}\mathbf{\right)}$ and their standard deviations (s.d.) for different WGP models
Trait  h ^{ 2 }  RRBLUP  LASSO  Elastic net  RKHS  BayesB  

${\mathbf{r}}_{\mathbf{(}\mathbf{g}\mathbf{,}\widehat{\mathbf{g}}\mathbf{)}}$  s.d.  ${\mathbf{r}}_{\mathbf{(}\mathbf{g}\mathbf{,}\widehat{\mathbf{g}}\mathbf{)}}$  s.d.  ${\mathbf{r}}_{\mathbf{(}\mathbf{g}\mathbf{,}\widehat{\mathbf{g}}\mathbf{)}}$  s.d.  ${\mathbf{r}}_{\mathbf{(}\mathbf{g}\mathbf{,}\widehat{\mathbf{g}}\mathbf{)}}$  s.d.  ${\mathbf{r}}_{\mathbf{(}\mathbf{g}\mathbf{,}\widehat{\mathbf{g}}\mathbf{)}}$  s.d.  
Dry matter yield  0.93  0.61  0.07  0.51  0.11  0.56  0.08  0.61  0.07  0.59  0.08 
Plant height  0.97  0.57  0.09  0.45  0.11  0.48  0.11  0.57  0.09  0.56  0.08 
Lignin content  0.88  0.69  0.07  0.60  0.08  0.60  0.10  0.68  0.07  0.58  0.09 
Dopamine  0.97  0.74  0.06  0.79  0.06  0.79  0.06  0.74  0.07  0.75  0.06 
Ribitol  0.95  0.49  0.12  0.61  0.10  0.63  0.10  0.50  0.10  0.50  0.11 
719700204  0.96  0.79  0.06  0.82  0.05  0.82  0.05  0.80  0.05  0.80  0.08 
Prediction accuracies of WGP ranged between 0.45 and 0.82 with standard deviations between 0.05 to 0.12 across traits and models (Table 3). The largest differences in accuracies between models ranged from 0.05 to 0.14 for the same trait. Between RRBLUP and RKHS, we found no difference in the prediction accuracies above 0.01 for any trait.
For agronomic traits, prediction accuracies were highest for RRBLUP with a drop of 0.09 to 0.12 if LASSO or elastic net was used and with a drop of 0.01 to 0.11 if BayesB was used.
Discussion
We found in a diverse panel of elite maize inbred lines that prediction accuracies obtained with five different WGP models were remarkable similar, even for traits with drastically deviating genetic architecture. Our results suggest that small gains in accuracies (up to 0.14) can be gained if the WGP model is selected according to the genetic architecture underlying the trait.
Recently, Heslot et al.[18] reported similar small differences for seven parametric WGP models when comparing them for different presumable highly polygenic agronomic traits over eight datasets of barley, Arabidopsis thaliana, maize, and wheat. For the metabolites, however, our results differ from those obtained from Clark et al.[19], who investigated the influence of genetic architecture on prediction accuracies achieved by RRBLUP or BayesB. Whereas these authors found only slight differences for simulated traits with a genetic architecture close to the infinitesimal genetic model, BayesB outperformed RRBLUP by an increase in prediction accuracy of ≈0.4 if the trait is controlled by either a few common or a few rare QTL. Simulation also predicted a drop in prediction accuracy in case of RRBLUP for traits controlled by a small number of QTL [44]. Although LASSO, elastic net, and BayesB showed higher accuracies compared to RRBLUP for metabolites, we found the differences to be remarkable small in case of LASSO or elastic net and negliable in the case of BayesB.
One major reason of the minor differences in prediction accuracies among the different models lies in the high level of LD found in elite breeding germplasm of maize. Our results suggest that with this level of LD (r^{2}=0.1 at ≈ 500 kb), accuracies are quite similar irrespective whether the effect of large QTL are precisely captured (as in the case of LASSO, elastic net, or BayesB) or spread over a larger region (as in the case of RRBLUP and RKHS). Since our population was highly diverse for elite maize germplasm in Europe, it is unlikely that breeders are confronted with lower levels of LD unless they work with highly exotic germplasm for which LD has been reported to decline within 510 kb [45].
Moreover, the high similarity of RKHS and RRBLUP suggest that either (i) nonadditive, epistatic genetic effects are not present, (ii) these are so small that they are negligible in WGP for the investigated traits, or (iii) RKHS regression is unable to capture them. In either case, for prediction purposes RKHS does not seem to yield any advancements over RRBLUP for situations comparable to our germplasm and traits. Dominance, as another source of nonadditive genetic, effects cannot be present in the inbred lines investigated in this study. For predicting heterozygeous F_{1} maize hybrids, however, it has been shown that modeling dominance effects can result in higher prediction accuracies [8].
Although BayesB reached for 5 of the 6 traits a higher prediction accuracy than the worst model, we cannot recommend it because of the excessively larger computation time and the negliable differences in prediction accuracies compared with RRBLUP in case of the metabolites as the result of probably only sampling error.
We found the approach to partition genetic variance over chromosomes useful for guiding the breeder which WGP model to prefer in the case of little or no prior knowledge on the genetic architecture. Whereas for the agronomic traits an approximately linear increase of cumulative explained genetic variance matched with a superiority of the L_{2} penalty (RRBLUP), the L_{1}penalty (LASSO) or a mixture of both penalties (elastic net) performed better in the case of the metabolites with a strong convex curve curvature (Figure 2A). Although for dry matter yield and plant height, barely significant association signals with a proportion of explained genetic variance <9% led to a chromosomal genetic variance slightly above the range expected from length of the chromosome (Figure 2B), these effects were too small to justify the use of the elastic net or LASSO.
As an alternative to this approach, Hayes et al.[20] estimated successively the genetic variance explained by each chromosome segment and compared it with the genetic variance captured by the remaining part of the genome. To correct for the nonindependence of neighbouring segments, they applied a bias correction using an expectation maximization (EM) algorithm. Such a correction is not necessary if the variance components for all chromosomes are estimated simultaneously as applied in this study; this is a further advantage besides its straightforward implementation using standard mixed model software packages such as ASReml.
Conclusions
Our empirical data of WGP in a large panel of diverse maize inbred lines suggest that (i) different WGP models differ only slightly in their prediction accuracies, irrespective of the number and effects of QTL found in association analysis, (ii) small gains in prediction accuracies can be obtained if the WGP model is selected according to the genetic architecture of the trait, (iii) genome partitioning of genetic variance offers a straightforward approach for model selection if the genetic architecture is unknown. The question of which WGP model to choose is therefore not expected to hamper implementation of WGP in maize breeding.
Declarations
Acknowledgements
We thank the staff of the experimental stations of the University of Hohenheim for conducting the field experiments. We thank the groups of Mark Stitt and Lothar Willmitzer of the Max Planck Institute of Molecular Plant Physiology for performing the metabolic profiling. This research was funded by the German Federal Ministry of Education and Research (BMBF) within the project GABIEnergy (FKZ: 0315045) and the AgroClustEr ’Synbreed  Synergistic plant and animal breeding’ (FKZ: 0315528D).
Authors’ Affiliations
References
 Morrell PL, Buckler ES, RossIbarra J: Crop genomics: advances and applications. Nat Rev Genet. 2011, 13: 8596.PubMed
 Heffner E, Sorrells M, Jannink JL: Genomic Selection for Crop Improvement. Crop Sci. 2009, 49: 112. 10.2135/cropsci2008.08.0512.View Article
 Jannink JL, Lorenz AJ, Iwata H: Genomic selection in plant breeding: from theory to practice. Brief Funct Genomics. 2010, 9: 166177. 10.1093/bfgp/elq001.View ArticlePubMed
 Meuwissen THE, Hayes BJ, Goddard ME: Prediction of total genetic value using genomewide dense marker maps. Genetics. 2001, 157: 18191829.PubMed CentralPubMed
 Hayes BJ, Bowman PJ, Chamberlain AJ, Goddard ME: Invited review: Genomic selection in dairy cattle: progress and challenges. J. Dairy Sci. 2009, 92: 433443. 10.3168/jds.20081646.View ArticlePubMed
 Albrecht T, Wimmer V, Auinger HJ, Erbe M, Knaak C, Ouzunova M, Simianer H, Schön CC: Genomebased prediction of testcross values in maize. Theor Appl Genet. 2011, 123: 339350. 10.1007/s0012201115877.View ArticlePubMed
 Zhao Y, Gowda M, Liu W, Würschrum T, Maurer HP, Longin FH, Ranc N, Reif JC: Accuracy of genomic selection in European maize elite breeding populations. Theor Appl Genet. 2012, 124: 769776. 10.1007/s001220111745y.View ArticlePubMed
 Technow F, Riedelsheimer C, Schrag TA, Melchinger AE: Genomic prediction of hybrid performance in maize with models incorporating dominance and population specific marker effects. Theor Appl Genet. 2012, 125: 11811194. 10.1007/s0012201219058.View ArticlePubMed
 Riedelsheimer C, CzedikEysenberg A, Grieder C, Lisec J, Technow F, Sulpice R, Altmann T, Stitt M, Willmitzer L, Melchinger AE: Genomic and metabolic prediction of complex heterotic traits in hybrid maize. Nat Genet. 2012, 44: 217220. 10.1038/ng.1033.View ArticlePubMed
 Fisher R: The correlation between relatives on the supposition of Mendelian inheritance. Trans Roy Soc Edinburgh. 1918, 52: 399433.View Article
 Simic D, Mladenovic Drinic S, Zdunic Z, Jambrovic A, Ledencan T, Brkic J, Brkic A, Brkic I: Quantitative trait Loci for biofortification traits in maize grain. J Hered. 2012, 103: 4754. 10.1093/jhered/esr122.View ArticlePubMed
 Wong JC, Lambert RJ, Wurtzel ET, Rocheford TR: QTL and candidate genes phytoene synthase and zetacarotene desaturase associated with the accumulation of carotenoids in maize. Theor Appl Genet. 2004, 108: 349359. 10.1007/s0012200314364.View ArticlePubMed
 Wisser RJ, BalintKurti PJ, Nelson RJ: The genetic architecture of disease resistance in maize: a synthesis of published studies. Phytopathology. 2006, 96: 120129. 10.1094/PHYTO960120.View ArticlePubMed
 Martin M, Miedaner T, Dhillon BS, Ufermann U, Kessel B, Ouzunova M, Schipprack W, Melchinger AE: Colocalization of QTL for Gibberella Ear Rot Resistance and Low Mycotoxin Contamination in Early European Maize. Crop Sci. 2011, 51: 19351945. 10.2135/cropsci2010.11.0664.View Article
 Martin M, Miedaner T, Schwegler DD, Kessel B, Ouzunova M, Dhillon BS, Schipprack W, Utz HF, Melchinger AE: Comparative Quantitative Trait Loci Mapping for Gibberella Ear Rot Resistance and Reduced Deoxynivalenol Contamination across Connected Maize Populations. Crop Sci. 2012, 52: 3243. 10.2135/cropsci2011.04.0214.View Article
 Welz H, Geiger H: Genes for resistance to northern corn leaf blight in diverse maize populations. Plant Breeding. 2000, 119: 114. 10.1046/j.14390523.2000.00462.x.View Article
 Poland JA, Bradbury PJ, Buckler E, Nelson RJ: Genomewide nested association mapping of quantitative resistance to northern leaf blight in maize. Proc Natl Acad Sci USA. 2011, 108: 68936898. 10.1073/pnas.1010894108.PubMed CentralView ArticlePubMed
 Heslot N, Yang HP, Sorrells ME, Jannink JL: Genomic Selection in Plant Breeding: A Comparison of Models. Crop Sci. 2012, 52: 146160.View Article
 Clark SA, Hickey JM, van der Werf JH: Different models of genetic variation and their effect on genomic evaluation. Genet Sel Evol. 2011, 43: 1810.1186/129796864318.PubMed CentralView ArticlePubMed
 Hayes BJ, Pryce J, Chamberlain AJ, Bowman PJ, Goddard ME: Genetic architecture of complex traits and accuracy of genomic prediction: coat colour, milkfat percentage, and type in Holstein cattle as contrasting model traits. PLoS Genet. 2010, 6: e100113910.1371/journal.pgen.1001139.PubMed CentralView ArticlePubMed
 EyreWalker A, Gaut RL, Hilton H, Feldman DL, Gaut BS: Investigation of the bottleneck leading to the domestication of maize. Proc Natl Acad Sci USA. 1998, 95: 44414446. 10.1073/pnas.95.8.4441.PubMed CentralView ArticlePubMed
 Duvick D, Smith J, Cooper M: LongTerm Selection in a Commercial Hybrid Maize Breeding Program. Plant Breed Rev. 2004, 24: 109151.
 Riedelsheimer C, Lisec J, CzedikEysenberg A, Sulpice R, Flis A, Grieder C, Altmann T, Stitt M, Willmitzer L, Melchinger AE: Genomewide association mapping of leaf metabolic profiles for dissecting complex traits in maize. Proc Natl Acad Sci USA. 2012, 109: 88728877. 10.1073/pnas.1120813109.PubMed CentralView ArticlePubMed
 Grieder C, Mittweg G, Dhillon BS, Montes JM, Orsini E, Melchinger AE: Determination of methane fermentation yield and its kinetics by near infrared spectroscopy and chemical composition in maize. JNIRS. 2011, 19: 463477.
 Grieder C, Dhillon B, Schipprack W, Melchinger A: Breeding maize as biogas substrate in Central Europe: II. Quantitativegenetic parameters for inbred lines and correlations with testcross performance. Theor Appl Genet. 2012, 124: 981988. 10.1007/s001220111762x.View ArticlePubMed
 Ganal MW, Durstewitz G, Polley A, Bérard A, Buckler ES, Charcosset A, Clarke JD, Graner EM, Hansen M, Joets J, Le Paslier MC, McMullen MD, Montalent P, Rose M, Schön CC, Sun Q, Walter H, Martin OC, Falque M: A Large Maize (Zea mays L.) SNP Genotyping Array: Development and Germplasm Genotyping, and Genetic Mapping to Compare with the B73 Reference Genome. PLoS ONE. 2011, 6: e2833410.1371/journal.pone.0028334.PubMed CentralView ArticlePubMed
 Lisec J, Schauer N, Kopka J, Willmitzer L, Fernie AR: Gas chromatography mass spectrometrybased metabolite profiling in plants. Nat Prot. 2006, 1: 387396. 10.1038/nprot.2006.59.View Article
 Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham JM, de Andrade M, Feenstra B, Feingold E, Hayes MG, Hill WG, Landi MT, Alonso A, Lettre G, Lin P, Ling H, Lowe W, Mathias RA, Melbye M, Pugh E, Cornelis MC, Weir BS, Goddard ME, Visscher PM: Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011, 43: 519525. 10.1038/ng.823.PubMed CentralView ArticlePubMed
 Butler D, Cullis B, Gilmour A, Gogel B: ASRemlR reference manual, version 3. 2009, Queensland, Australia: The Department of Primary Industries and Fisheries
 Higham N: Computing the nearest correlation matrix — a problem from finance. IMA J Numer Anal. 2002, 22: 329343. 10.1093/imanum/22.3.329.View Article
 Bates D, Maechler M: Matrix: Sparse and Dense Matrix Classes and Methods. 2012, [http://cran.rproject.org/package=Matrix]
 de los Campos G, Naya H, Gianola D, Crossa J, Legarra A, Manfredi E, Weigel K, Cotes JM: Predicting quantitative traits with regression models for dense molecular markers and pedigree. Genetics. 2009, 182: 375385. 10.1534/genetics.109.101501.PubMed CentralView ArticlePubMed
 Zou H, Hastie T: Regularization and variable selection via the elastic net. J R Statist Soc B. 2005, 67: 301320. 10.1111/j.14679868.2005.00503.x.View Article
 Gianola D, van Kaam JBCHM: Reproducing kernel hilbert spaces regression methods for genomic assisted prediction of quantitative traits. Genetics. 2008, 178: 22892303. 10.1534/genetics.107.084285.PubMed CentralView ArticlePubMed
 Ruppert D, Wand M, Carroll R: Semiparametric Regression. 2003, New York, NY, USA: Cambridge University PressView Article
 Piepho H: Ridge regression and extensions for genomewide selection in maize. Crop Sci. 2009, 49: 11651176. 10.2135/cropsci2008.10.0595.View Article
 Yang J, Lee SH, Goddard ME, Visscher PM: GCTA: a tool for genomewide complex trait analysis. Am J Hum Genet. 2011, 88: 7682. 10.1016/j.ajhg.2010.11.011.PubMed CentralView ArticlePubMed
 Hastie T, Tibshirani R, Friedman J: The Elements of Statistical Learning. 2009, New York, NY, USA: Springer Science+Business MediaView Article
 Tibshirani R: Regression shrinkage and selection via the lasso. J R Statist Soc B. 1996, 58: 267288.
 Li Z, Sillanpää MJ: Overview of LASSOrelated penalized regression methods for quantitative trait mapping and genomic selection. Theor Appl Genet. 2012, 125: 419435. 10.1007/s0012201218929.View ArticlePubMed
 Ogutu JO, SchulzStreeck T, Piepho HP: Genomic selection using regularized linear regression models: ridge regression, lasso, elastic net and their extensions. BMC Proceedings. 2012, 6 (Suppl 2): S1010.1186/175365616S2S10.PubMed CentralView ArticlePubMed
 Friedman J, Hastie T, Tibshirani R: Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010, 33: 122.PubMed CentralView ArticlePubMed
 Yang W, Tempelman RJ: A Bayesian Antedependence Model for Whole Genome Prediction. Genetics. 2011, 190: 14911501.View ArticlePubMed
 Habier D, Fernando RL, Dekkers JCM: The impact of genetic relationship information on genomeassisted breeding values. Genetics. 2007, 177: 23892397.PubMed CentralPubMed
 Yan J, Shah T, Warburton M, Buckler E, McMullen M, Crouch J: Genetic characterization and linkage disequilibrium estimation of a global maize collection using SNP markers. PLoS ONE. 2009, 4: e845110.1371/journal.pone.0008451.PubMed CentralView ArticlePubMed
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.