Estimate Marker Effects Research Articles

Maximizing crossbred performance through purebred genomic selection.

BackgroundIn livestock production, many animals are crossbred, with two distinct advantages: heterosis and breed complementarity. Genomic selection (GS) can be used to select purebred parental lines for crossbred performance (CP). Dominance being the likely genetic basis of heterosis, explicitly including dominance in the GS model may be an advantage to select purebreds for CP. Estimated breeding values for CP can be calculated from additive and dominance effects of alleles that are estimated using pure line data. The objective of this simulation study was to investigate the benefits of applying GS to select purebred animals for CP, based on purebred phenotypic and genotypic information. A second objective was to compare the use of two separate pure line reference populations to that of a single reference population that combines both pure lines. These objectives were investigated under two conditions, i.e. either a low or a high correlation of linkage disequilibrium (LD) phase between the pure lines.ResultsThe results demonstrate that the gain in CP was higher when parental lines were selected for CP, rather than purebred performance, both with a low and a high correlation of LD phase. For a low correlation of LD phase between the pure lines, the use of two separate reference populations yielded a higher gain in CP than use of a single reference population that combines both pure lines. However, for a high correlation of LD phase, marker effects that were estimated using a single combined reference population increased the gain in CP.ConclusionsUnder the hypothesis that performance of crossbred animals differs from that of purebred animals due to dominance, a dominance model can be used for GS of purebred individuals for CP, without using crossbred data. Furthermore, if the correlation of LD phase between pure lines is high, accuracy of selection can be increased by combining the two pure lines into a single reference population to estimate marker effects.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-015-0099-3) contains supplementary material, which is available to authorized users.

Upweighting rare favourable alleles increases long-term genetic gain in genomic selection programs.

BackgroundThe short-term impact of using different genomic prediction (GP) models in genomic selection has been intensively studied, but their long-term impact is poorly understood. Furthermore, long-term genetic gain of genomic selection is expected to improve by using Jannink’s weighting (JW) method, in which rare favourable marker alleles are upweighted in the selection criterion. In this paper, we extend the JW method by including an additional parameter to decrease the emphasis on rare favourable alleles over the time horizon, with the purpose of further improving the long-term genetic gain. We call this new method dynamic weighting (DW). The paper explores the long-term impact of different GP models with or without weighting methods.MethodsDifferent selection criteria were tested by simulating a population of 500 animals with truncation selection of five males and 50 females. Selection criteria included unweighted and weighted genomic estimated breeding values using the JW or DW methods, for which ridge regression (RR) and Bayesian lasso (BL) were used to estimate marker effects. The impacts of these selection criteria were compared under three genetic architectures, i.e. varying numbers of QTL for the trait and for two time horizons of 15 (TH15) or 40 (TH40) generations.ResultsFor unweighted GP, BL resulted in up to 21.4% higher long-term genetic gain and 23.5% lower rate of inbreeding under TH40 than RR. For weighted GP, DW resulted in 1.3 to 5.5% higher long-term gain compared to unweighted GP. JW, however, showed a 6.8% lower long-term genetic gain relative to unweighted GP when BL was used to estimate the marker effects. Under TH40, both DW and JW obtained significantly higher genetic gain than unweighted GP. With DW, the long-term genetic gain was increased by up to 30.8% relative to unweighted GP, and also increased by 8% relative to JW, although at the expense of a lower short-term gain.ConclusionsIrrespective of the number of QTL simulated, BL is superior to RR in maintaining genetic variance and therefore results in higher long-term genetic gain. Moreover, DW is a promising method with which high long-term genetic gain can be expected within a fixed time frame.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-015-0101-0) contains supplementary material, which is available to authorized users.

Redesigning the exploitation of wheat genetic resources.

More than half a million wheat genetic resources are resting in gene banks worldwide. Unlocking their hidden favorable genetic diversity for breeding is pivotal for enhancing grain yield potential, and averting future food shortages. Here, we propose exploiting recent advances in hybrid wheat technology to uncover the masked breeding values of wheat genetic resources. The gathered phenotypic information will enable a targeted choice of accessions with high value for pre-breeding among this plethora of genetic resources. We intend to provoke a paradigm shift in pre-breeding strategies for grain yield, moving away from allele mining toward genome-wide selection to bridge the yield gap between genetic resources and elite breeding pools.

A single-step genomic model with direct estimation of marker effects

Compared with the currently widely used multi-step genomic models for genomic evaluation, single-step genomic models can provide more accurate genomic evaluation by jointly analyzing phenotypes and genotypes of all animals and can properly correct for the effect of genomic preselection on genetic evaluations. The objectives of this study were to introduce a single-step genomic model, allowing a direct estimation of single nucleotide polymorphism (SNP) effects, and to develop efficient computing algorithms for solving equations of the single-step SNP model. We proposed an alternative to the current single-step genomic model based on the genomic relationship matrix by including an additional step for estimating the effects of SNP markers. Our single-step SNP model allowed flexible modeling of SNP effects in terms of the number and variance of SNP markers. Moreover, our single-step SNP model included a residual polygenic effect with trait-specific variance for reducing inflation in genomic prediction. A kernel calculation of the SNP model involved repeated multiplications of the inverse of the pedigree relationship matrix of genotyped animals with a vector, for which numerical methods such as preconditioned conjugate gradients can be used. For estimating SNP effects, a special updating algorithm was proposed to separate residual polygenic effects from the SNP effects. We extended our single-step SNP model to general multiple-trait cases. By taking advantage of a block-diagonal (co)variance matrix of SNP effects, we showed how to estimate multivariate SNP effects in an efficient way. A general prediction formula was derived for candidates without phenotypes, which can be used for frequent, interim genomic evaluations without running the whole genomic evaluation process. We discussed various issues related to implementation of the single-step SNP model in Holstein populations with an across-country genomic reference population.

DAIRRy-BLUP: A High-Performance Computing Approach to Genomic Prediction

In genomic prediction, common analysis methods rely on a linear mixed-model framework to estimate SNP marker effects and breeding values of animals or plants. Ridge regression-best linear unbiased prediction (RR-BLUP) is based on the assumptions that SNP marker effects are normally distributed, are uncorrelated, and have equal variances. We propose DAIRRy-BLUP, a parallel, Distributed-memory RR-BLUP implementation, based on single-trait observations ( Y: ), that uses the Average Information algorithm for restricted maximum-likelihood estimation of the variance components. The goal of DAIRRy-BLUP is to enable the analysis of large-scale data sets to provide more accurate estimates of marker effects and breeding values. A distributed-memory framework is required since the dimensionality of the problem, determined by the number of SNP markers, can become too large to be analyzed by a single computing node. Initial results show that DAIRRy-BLUP enables the analysis of very large-scale data sets (up to 1,000,000 individuals and 360,000 SNPs) and indicate that increasing the number of phenotypic and genotypic records has a more significant effect on the prediction accuracy than increasing the density of SNP arrays.

Genomic prediction for rust resistance in diverse wheat landraces.

We have demonstrated that genomic selection in diverse wheat landraces for resistance to leaf, stem and strip rust is possible, as genomic breeding values were moderately accurate. Markers with large effects in the Bayesian analysis confirmed many known genes, while also discovering many previously uncharacterised genome regions associated with rust scores. Genomic selection, where selection decisions are based on genomic estimated breeding values (GEBVs) derived from genome-wide DNA markers, could accelerate genetic progress in plant breeding. In this study, we assessed the accuracy of GEBVs for rust resistance in 206 hexaploid wheat (Triticum aestivum) landraces from the Watkins collection of phenotypically diverse wheat genotypes from 32 countries. The landraces were genotyped for 5,568 SNPs using an Illumina iSelect 9 K bead chip assay and phenotyped for field-based leaf rust (Lr), stem rust (Sr) and stripe rust (Yr) responses across multiple years. Genomic Best Linear Unbiased Prediction (GBLUP) and a Bayesian Regression method (BayesR) were used to predict GEBVs. Based on fivefold cross-validation, the accuracy of genomic prediction averaged across years was 0.35, 0.27 and 0.44 for Lr, Sr and Yr using GBLUP and 0.33, 0.38 and 0.30 for Lr, Sr and Yr using BayesR, respectively. Inclusion of PCR-predicted genotypes for known rust resistance genes increased accuracy more substantially when the marker was diagnostic (Lr34/Sr57/Yr18) for the presence-absence of the gene rather than just linked (Sr2). Investigation of the impact of genetic relatedness between validation and reference lines on accuracy of genomic prediction showed that accuracy will be higher when each validation line had at least one close relationship to the reference lines. Overall, the prediction accuracies achieved in this study are encouraging, and confirm the feasibility of genomic selection in wheat. In several instances, estimated marker effects were confirmed by published literature and results of mapping experiments using Watkins accessions.

Heteroscedastic Ridge Regression Approaches for Genome-Wide Prediction With a Focus on Computational Efficiency and Accurate Effect Estimation

Ridge regression with heteroscedastic marker variances provides an alternative to Bayesian genome-wide prediction methods. Our objectives were to suggest new methods to determine marker-specific shrinkage factors for heteroscedastic ridge regression and to investigate their properties with respect to computational efficiency and accuracy of estimated effects. We analyzed published data sets of maize, wheat, and sugar beet as well as simulated data with the new methods. Ridge regression with shrinkage factors that were proportional to single-marker analysis of variance estimates of variance components (i.e., RRWA) was the fastest method. It required computation times of less than 1 sec for medium-sized data sets, which have dimensions that are common in plant breeding. A modification of the expectation-maximization algorithm that yields heteroscedastic marker variances (i.e., RMLV) resulted in the most accurate marker effect estimates. It outperformed the homoscedastic ridge regression approach for best linear unbiased prediction in particular for situations with high marker density and strong linkage disequilibrium along the chromosomes, a situation that occurs often in plant breeding populations. We conclude that the RRWA and RMLV approaches provide alternatives to the commonly used Bayesian methods, in particular for applications in which computational feasibility or accuracy of effect estimates are important, such as detection or functional analysis of genes or planning crosses.

Increasing long-term response by selecting for favorable minor alleles.

Long-term response of genomic selection can be improved by considering allele frequencies of selected markers or quantitative trait loci (QTLs). A previous formula to weight allele frequency of favorable minor alleles was tested, and 2 new formulas were developed. The previous formula used nonlinear weights based on square root of frequency of the favorable allele. The new formulas included a parameter δ to balance long- and short-term progress; one used square root and the other used simple linear weights. The formulas were tested by simulation of 20 generations (population size of 3,000 for each generation) with direct selection on 3,000 QTLs (100 per chromosome). A QTL distribution with normally distributed allele effects and a heavy-tailed distribution were tested. Optimum δ from simulation was applied to data from Holstein, Jersey and Brown Swiss dairy cattle to compare differences of adjusted and official genomic evaluations. From simulation, optimum δ was 0.4 for the heavy-tailed QTL distribution but only 0.1 or 0.2 for a normal distribution. The previous formula had slower response than unweighted selection in early generations and did not recover by generation 20. Long-term response was slightly greater with the new formulas than with unweighted selection; the linear formula may be best for routine use because of more progress in early generations compared to nonlinear formula. Official and adjusted U.S. evaluations based on actual genotypes and estimated marker effects were correlated by 0.994 for Holsteins and Jerseys and 0.989 for Brown Swiss using linear weighting of allele frequency, which was higher than nonlinear weighting. The difference between adjusted and official evaluations was highly correlated negatively with an animal’s average genomic relationship to the population. Thus, strategies to reduce genomic inbreeding may achieve almost as much long-term progress as selection of favorable minor alleles.

Accuracy of genomic prediction within and across populations for nematode resistance and body weight traits in sheep

Genomic prediction utilizes single nucleotide polymorphism (SNP) chip data to predict animal genetic merit. It has the advantage of potentially capturing the effects of the majority of loci that contribute to genetic variation in a trait, even when the effects of the individual loci are very small. To implement genomic prediction, marker effects are estimated with a training set, including individuals with marker genotypes and trait phenotypes; subsequently, genomic estimated breeding values (GEBV) for any genotyped individual in the population can be calculated using the estimated marker effects. In this study, we aimed to: (i) evaluate the potential of genomic prediction to predict GEBV for nematode resistance traits and BW in sheep, within and across populations; (ii) evaluate the accuracy of these predictions through within-population cross-validation; and (iii) explore the impact of population structure on the accuracy of prediction. Four data sets comprising 752 lambs from a Scottish Blackface population, 2371 from a Sarda×Lacaune backcross population, 1000 from a Martinik Black-Belly×Romane backcross population and 64 from a British Texel population were used in this study. Traits available for the analysis were faecal egg count for Nematodirus and Strongyles and BW at different ages or as average effect, depending on the population. Moreover, immunoglobulin A was also available for the Scottish Blackface population. Results show that GEBV had moderate to good within-population predictive accuracy, whereas across-population predictions had accuracies close to zero. This can be explained by our finding that in most cases the accuracy estimates were mostly because of additive genetic relatedness between animals, rather than linkage disequilibrium between SNP and quantitative trait loci. Therefore, our results suggest that genomic prediction for nematode resistance and BW may be of value in closely related animals, but that with the current SNP chip genomic predictions are unlikely to work across breeds.

Molecular approaches for wheat improvement under drought and heat stress

Novel molecular breeding selection strategies were implemented in breeding programs in India for wheat improvement under drought and heat stress. Elite Indian varieties HD2733 and GW322 were targeted for improvement to abiotic stresses through marker assisted backcross breeding approach. Backcross populations were advanced to BC1/BC2 F2 after tracking QTLs for foreground selection covering traits such as canopy temperature (CT), chlorophyll content, normalized difference vegetation index (NDVI), grain filling duration (GFD), thousand kernel weight, grain yield and a range of polymorphic microsatellite markers covering entire genome (4–5 SSR markers per chromosome) being used for tracking the recurrent parent genome. Marker assisted recurrent selection (MARS) involved estimation of marker effects for several small effect/major QTLs followed by two recombination cycles. Designed intermating among selected F5 families were carried out after conducting ANOVA and AMMI analysis on multi-location yield trials and polymorphic markers among the parents in two biparental base populations (F5). Inter-family intermating among the best identified families carried out in different combinations to accumulate and recombine 4–8 QTLs per intermated progeny and confirmed in selfed homozygous lines.

Accuracy of genomic predictions in Bos indicus (Nellore) cattle

BackgroundNellore cattle play an important role in beef production in tropical systems and there is great interest in determining if genomic selection can contribute to accelerate genetic improvement of production and fertility in this breed. We present the first results of the implementation of genomic prediction in a Bos indicus (Nellore) population.MethodsInfluential bulls were genotyped with the Illumina Bovine HD chip in order to assess genomic predictive ability for weight and carcass traits, gestation length, scrotal circumference and two selection indices. 685 samples and 320 238 single nucleotide polymorphisms (SNPs) were used in the analyses. A forward-prediction scheme was adopted to predict the genomic breeding values (DGV). In the training step, the estimated breeding values (EBV) of bulls were deregressed (dEBV) and used as pseudo-phenotypes to estimate marker effects using four methods: genomic BLUP with or without a residual polygenic effect (GBLUP20 and GBLUP0, respectively), a mixture model (Bayes C) and Bayesian LASSO (BLASSO). Empirical accuracies of the resulting genomic predictions were assessed based on the correlation between DGV and dEBV for the testing group.ResultsAccuracies of genomic predictions ranged from 0.17 (navel at weaning) to 0.74 (finishing precocity). Across traits, Bayesian regression models (Bayes C and BLASSO) were more accurate than GBLUP. The average empirical accuracies were 0.39 (GBLUP0), 0.40 (GBLUP20) and 0.44 (Bayes C and BLASSO). Bayes C and BLASSO tended to produce deflated predictions (i.e. slope of the regression of dEBV on DGV greater than 1). Further analyses suggested that higher-than-expected accuracies were observed for traits for which EBV means differed significantly between two breeding subgroups that were identified in a principal component analysis based on genomic relationships.ConclusionsBayesian regression models are of interest for future applications of genomic selection in this population, but further improvements are needed to reduce deflation of their predictions. Recurrent updates of the training population would be required to enable accurate prediction of the genetic merit of young animals. The technical feasibility of applying genomic prediction in a Bos indicus (Nellore) population was demonstrated. Further research is needed to permit cost-effective selection decisions using genomic information.

Pedigree and genomic analyses of feed consumption and residual feed intake in laying hens

Efficiency of production is increasingly important with the current escalation of feed costs and demands to minimize the environmental footprint. The objectives of this study were 1) to estimate heritabilities for daily feed consumption and residual feed intake and their genetic correlations with production and egg-quality traits; 2) to evaluate accuracies of estimated breeding values from pedigree- and marker-based prediction models; and 3) to localize genomic regions associated with feed efficiency in a brown egg layer line. Individual feed intake data collected over 2-wk trial periods were available for approximately 6,000 birds from 8 generations. Genetic parameters were estimated with a multitrait animal model; methods BayesB and BayesCπ were used to estimate marker effects and find genomic regions associated with feed efficiency. Using pedigree information, feed efficiency was found to be moderately heritable (h2 = 0.46 for daily feed consumption and 0.47 for residual feed intake). Hens that consumed more feed and had greater residual feed intake (lower efficiency) had a genetic tendency to lay slightly more eggs with greater yolk weights and albumen heights. Regions on chromosomes 1, 2, 4, 7, 13, and Z were found to be associated with feed intake and efficiency. The accuracy from genomic prediction was higher and more persistent (better maintained across generations) than that from pedigree-based prediction. These results indicate that genomic selection can be used to improve feed efficiency in layers.

Genomic Selection: Status in Different Species and Challenges for Breeding

Technical advances and development in the market for genomic tools have facilitated access to whole-genome data across species. Building-up on the acquired knowledge of the genome sequences, large-scale genotyping has been optimized for broad use, so genotype information can be routinely used to predict genetic merit. Genomic selection (GS) refers to the use of aggregates of estimated marker effects as predictors which allow improved individual differentiation at young age. Realizable benefits of GS are influenced by several factors and vary in quantity and quality between species. General characteristics and challenges of GS in implementation and routine application are described, followed by an overview over the current status of its use, prospects and challenges in important animal species. Genetic gain for a particular trait can be enhanced by shortening of the generation interval, increased selection accuracy and increased selection intensity, with species- and breed-specific relevance of the determinants. Reliable predictions based on genetic marker effects require assembly of a reference for linking of phenotype and genotype data to allow estimation and regular re-estimation. Experiences from dairy breeding have shown that international collaboration can set the course for fast and successful implementation of innovative selection tools, so genomics may significantly impact the structures of future breeding and breeding programmes. Traits of great and increasing importance, which were difficult to improve in the conventional systems, could be emphasized, if continuous availability of high-quality phenotype data can be assured. Equally elaborate strategies for genotyping and phenotyping will allow tailored approaches to balance efficient animal production, sustainability, animal health and welfare in future.

Prediction of Complex Human Traits Using the Genomic Best Linear Unbiased Predictor

Despite important advances from Genome Wide Association Studies (GWAS), for most complex human traits and diseases, a sizable proportion of genetic variance remains unexplained and prediction accuracy (PA) is usually low. Evidence suggests that PA can be improved using Whole-Genome Regression (WGR) models where phenotypes are regressed on hundreds of thousands of variants simultaneously. The Genomic Best Linear Unbiased Prediction (G-BLUP, a ridge-regression type method) is a commonly used WGR method and has shown good predictive performance when applied to plant and animal breeding populations. However, breeding and human populations differ greatly in a number of factors that can affect the predictive performance of G-BLUP. Using theory, simulations, and real data analysis, we study the performance of G-BLUP when applied to data from related and unrelated human subjects. Under perfect linkage disequilibrium (LD) between markers and QTL, the prediction R-squared (R2) of G-BLUP reaches trait-heritability, asymptotically. However, under imperfect LD between markers and QTL, prediction R2 based on G-BLUP has a much lower upper bound. We show that the minimum decrease in prediction accuracy caused by imperfect LD between markers and QTL is given by (1−b)2, where b is the regression of marker-derived genomic relationships on those realized at causal loci. For pairs of related individuals, due to within-family disequilibrium, the patterns of realized genomic similarity are similar across the genome; therefore b is close to one inducing small decrease in R2. However, with distantly related individuals b reaches very low values imposing a very low upper bound on prediction R2. Our simulations suggest that for the analysis of data from unrelated individuals, the asymptotic upper bound on R2 may be of the order of 20% of the trait heritability. We show how PA can be enhanced with use of variable selection or differential shrinkage of estimates of marker effects.

Accelerating improvement of livestock with genomic selection.

Three recent breakthroughs have resulted in the current widespread use of DNA information: the genomic selection (GS) methodology, which is a form of marker-assisted selection on a genome-wide scale, and the discovery of large numbers of single-nucleotide markers and cost effective methods to genotype them. GS estimates the effect of thousands of DNA markers simultaneously. Nonlinear estimation methods yield higher accuracy, especially for traits with major genes. The marker effects are estimated in a genotyped and phenotyped training population and are used for the estimation of breeding values of selection candidates by combining their genotypes with the estimated marker effects. The benefits of GS are greatest when selection is for traits that are not themselves recorded on the selection candidates before they can be selected. In the future, genome sequence data may replace SNP genotypes as markers. This could increase GS accuracy because the causative mutations should be included in the data.

Accuracy of direct genomic breeding values for nationally evaluated traits in US Limousin and Simmental beef cattle

BackgroundIn national evaluations, direct genomic breeding values can be considered as correlated traits to those for which phenotypes are available for traditional estimation of breeding values. For this purpose, estimates of the accuracy of direct genomic breeding values expressed as genetic correlations between traits and their respective direct genomic breeding values are required.MethodsWe derived direct genomic breeding values for 2239 registered Limousin and 2703 registered Simmental beef cattle genotyped with either the Illumina BovineSNP50 BeadChip or the Illumina BovineHD BeadChip. For the 264 Simmental animals that were genotyped with the BovineHD BeadChip, genotypes for markers present on the BovineSNP50 BeadChip were extracted. Deregressed estimated breeding values were used as observations in weighted analyses that estimated marker effects to derive direct genomic breeding values for each breed. For each breed, genotyped individuals were clustered into five groups using K-means clustering, with the aim of increasing within-group and decreasing between-group pedigree relationships. Cross-validation was performed five times for each breed, using four groups for training and the fifth group for validation. For each trait, we then applied a weighted bivariate analysis of the direct genomic breeding values of genotyped animals from all five validation sets and their corresponding deregressed estimated breeding values to estimate variance and covariance components.ResultsAfter minimizing relationships between training and validation groups, estimated genetic correlations between each trait and its direct genomic breeding values ranged from 0.39 to 0.76 in Limousin and from 0.29 to 0.65 in Simmental. The efficiency of selection based on direct genomic breeding values relative to selection based on parent average information ranged from 0.68 to 1.28 in genotyped Limousin and from 0.51 to 1.44 in genotyped Simmental animals. The efficiencies were higher for 323 non-genotyped young Simmental animals, born after January 2012, and ranged from 0.60 to 2.04.ConclusionsDirect genomic breeding values show promise for routine use by Limousin and Simmental breeders to improve the accuracy of predicted genetic merit of their animals at a young age and increase response to selection. Benefits from selecting on direct genomic breeding values are greater for breeders who use natural mating sires in their herds than for those who use artificial insemination sires. Producers with unregistered commercial Limousin and Simmental cattle could also benefit from being able to identify genetically superior animals in their herds, an opportunity that has in the past been limited to seed stock animals.

Estimates of marker effects for measures of milk flow in the Italian brown Swiss dairy cattle population

BackgroundMilkability is a complex trait that is characterized by milk flow traits including average milk flow rate, maximum milk flow rate and total milking time. Milkability has long been recognized as an economically important trait that can be improved through selection. By improving milkability, management costs of milking decrease through reduced labor and improved efficiency of the automatic milking system, which has been identified as an important factor affecting net profit. The objective of this study was to identify markers associated with electronically measured milk flow traits, in the Italian Brown Swiss population that could potentially improve selection based on genomic predictions.ResultsSires (n = 1351) of cows with milk flow information were genotyped for 33,074 single nucleotide polymorphism (SNP) markers distributed across 29 Bos taurus autosomes (BTA). Among the six milk flow traits collected, ascending time, time of plateau, descending time, total milking time, maximum milk flow and average milk flow, there were 6,929 (time of plateau) to 14,585 (maximum milk flow) significant SNP markers identified for each trait across all BTA. Unique regions were found for each of the 6 traits providing evidence that each individual milk flow trait offers distinct genetic information about milk flow. This study was also successful in identifying functional processes and genes associated with SNPs that influences milk flow.ConclusionsIn addition to verifying the presence of previously identified milking speed quantitative trait loci (QTL) within the Italian Brown Swiss population, this study revealed a number of genomic regions associated with milk flow traits that have never been reported as milking speed QTL. While several of these regions were not associated with a known gene or QTL, a number of regions were associated with QTL that have been formerly reported as regions associated with somatic cell count, somatic cell score and udder morphometrics. This provides further evidence of the complexity of milk flow traits and the underlying relationship it has with other economically important traits for dairy cattle. Improved understanding of the overall milking pattern will aid in identification of cows with lower management costs and improved udder health.

An assessment of linkage disequilibrium in Holstein cattle using a Bayesian network

Linkage disequilibrium (LD) is defined as a non-random association of the distributions of alleles at different loci within a population. This association between loci is valuable in prediction of quantitative traits in animals and plants and in genome-wide association studies. A question that arises is whether standard metrics such as D' and r(2) reflect complex associations in a genetic system properly. It seems reasonable to take the view that loci associate and interact together as a system or network, as opposed to in a simple pairwise manner. We used a Bayesian network (BN) as a representation of choice for an LD network. A BN is a graphical depiction of a probability distribution and can represent sets of conditional independencies. Moreover, it provides a visual display of the joint distribution of the set of random variables in question. The usefulness of BN for linkage disequilibrium was explored and illustrated using genetic marker loci found to have the strongest effects on milk protein in Holstein cattle based on three strategies for ranking marker effect estimates: posterior means, standardized posterior means and additive genetic variance. Two different algorithms, Tabu search (a local score-based algorithm) and incremental association Markov blanket (a constraint-based algorithm), coupled with the chi-square test, were used for learning the structure of the BN and were compared with the reference r(2) metric represented as an LD heat map. The BN captured several genetic markers associated as clusters, implying that markers are inter-related in a complicated manner. Further, the BN detected conditionally dependent markers. The results confirm that LD relationships are of a multivariate nature and that r(2) gives an incomplete description and understanding of LD. Use of an LD Bayesian network enables inferring associations between loci in a systems framework and provides a more accurate picture of LD than that resulting from the use of pairwise metrics.

Potential and Optimization of Genomic Selection for Fusarium Head Blight Resistance in Six‐Row Barley

ABSTRACTFusarium head blight (FHB),is a devastating disease of barley (Hordeum vulgare L.), causing reductions in yield and quality. Marker‐based selection for resistance to FHB and lowered deoxynivalenol (DON) grain concentration would save considerable costs and time associated with phenotyping. A marker‐based selection approach called genomic selection (GS) uses genomewide marker information to predict genetic value. We used a cross‐validation approach that separated training sets from validation sets by both entry and environment. We used this framework to test the potential of GS for genetic improvement of FHB and DON as well as test the effect of different factors on prediction accuracy. Prediction accuracy for FHB was found to be as high as 0.72 and that for DON was found to be as high as 0.68. Little difference was found between marker effect estimation methods in terms of prediction of entry genetic value. The extensive linkage disequilibrium (LD) present in this population allowed the marker set to be reduced to 384 markers and training population (TP) size to be reduced 200 with little effect on prediction accuracy. We found little to no advantage to combining subpopulations that correspond to neighboring breeding programs to increase TP size. Apparently, little genetic information is shared between subpopulations, either because of different marker–quantitative trait loci (QTL) linkage phases, different segregating QTL, or nonadditive gene action.

Impact of selective genotyping in the training population on accuracy and bias of genomic selection

Estimating marker effects based on routinely generated phenotypic data of breeding programs is a cost-effective strategy to implement genomic selection. Truncation selection in breeding populations, however, could have a strong impact on the accuracy to predict genomic breeding values. The main objective of our study was to investigate the influence of phenotypic selection on the accuracy and bias of genomic selection. We used experimental data of 788 testcross progenies from an elite maize breeding program. The testcross progenies were evaluated in unreplicated field trials in ten environments and fingerprinted with 857 SNP markers. Random regression best linear unbiased prediction method was used in combination with fivefold cross-validation based on genotypic sampling. We observed a substantial loss in the accuracy to predict genomic breeding values in unidirectional selected populations. In contrast, estimating marker effects based on bidirectional selected populations led to only a marginal decrease in the prediction accuracy of genomic breeding values. We concluded that bidirectional selection is a valuable approach to efficiently implement genomic selection in applied plant breeding programs.