Skip to main content

The nature, scope and impact of genomic prediction in beef cattle in the United States

Abstract

Artificial selection has proven to be effective at altering the performance of animal production systems. Nevertheless, selection based on assessment of the genetic superiority of candidates is suboptimal as a result of errors in the prediction of genetic merit. Conventional breeding programs may extend phenotypic measurements on selection candidates to include correlated indicator traits, or delay selection decisions well beyond puberty so that phenotypic performance can be observed on progeny or other relatives. Extending the generation interval to increase the accuracy of selection reduces annual rates of gain compared to accurate selection and use of parents of the next generation at the immediate time they reach breeding age. Genomic prediction aims at reducing prediction errors at breeding age by exploiting information on the transmission of chromosome fragments from parents to selection candidates, in conjunction with knowledge on the value of every chromosome fragment. For genomic prediction to influence beef cattle breeding programs and the rate or cost of genetic gains, training analyses must be undertaken, and genomic prediction tools made available for breeders and other industry stakeholders. This paper reviews the nature or kind of studies currently underway, the scope or extent of some of those studies, and comments on the likely predictive value of genomic information for beef cattle improvement.

Background

Genetic improvement results from selection of above-average candidates as parents of the next generation. In a competitive market, above-average candidates would be those that improve consumer satisfaction, influencing immediate eating quality, purchase cost, long-term health implications of consumption, care of the environment in the production and processing of the beef; and welfare of the animals. Satisfied consumers demand and pay more for desirable beef, and under perfect competition this will be reflected along the production chain by increased farm-gate prices for cow-calf producers. Seedstock suppliers that sell bulls to cow-calf producers would be expected to respond by developing and implementing breeding programs that provide successive crops of bulls that outperform their predecessors.

Inspection of genetic trends, e.g. [1, 2], shows that beef cattle selection has resulted in animals with increased merit for early growth and improved rib eye area and marbling scores. There is no evidence for genetic improvement in reproductive performance. Selection has resulted in animals with larger mature size [1] and greater cow maintenance requirements [2], which increase production costs, as cow maintenance requirements are a major determinant of the total feed required in the production system [3]. Beef cattle selection has therefore failed in practice to achieve balanced improvement across the spectrum of traits that contribute to breeding goals. One reason has been our inability to cost-effectively rank selection candidates for all the attributes of interest [4]. This is the case because reliably quantifying the merits of animals in terms of their breeding values has been totally reliant on recording pedigree and performance information, primarily on the selection candidates themselves, their parents and perhaps their offspring. This has led to improvement programs that have been phenotype driven, i.e. programs that are focused on easy to measure traits that are recorded at young ages, such as early growth and ultrasound assessment of carcass attributes, rather than being goal driven and focused on all the attributes that influence consumer satisfaction [5]. The fundamental reason for this failure is that mixed model predictions of merit using the relationship matrix and applied to young animals can, with sufficient historical data, reliably predict the parent average (PA) effects, but are unable to predict the Mendelian sampling effects without having phenotypic observations on the individual or its descendants [6]. Accordingly, with only ancestral records, there is little information to discriminate among paternal half-sibs other than based on the merit of the dams. In that setting, it is seldom possible to identify young selection candidates with merit superior to existing selected sires. In the beef cattle context, this has led to low selection accuracy for mature size, lifetime reproductive performance, stayability/longevity, and disease resistance. Other important traits such as tenderness of beef, other aspects of eating quality, and feed efficiency, have had no prospects for selection as there are no phenotypic measures that can be readily and cost-effectively obtained on large numbers of seedstock animals.

Molecular-based information has long held promise to improve the prediction of young animals by first using phenotypic markers, second using microsatellite markers, and most recently using ever-increasing densities of single nucleotide polymorphisms (SNP).

Phenotypic markers such as blood groups were found to characterize the inheritance of certain chromosomal regions, proving useful for selection if that region contained a major gene responsible for variation in a trait of interest [7]. Unfortunately, there are insufficient simply inherited phenotypic attributes to characterize the entire genome.

Highly polymorphic microsatellite markers provided new opportunities to find major genes or quantitative trait loci (QTL) that influence important traits [8]. These markers that can have many alleles at each locus, can be informative in much of the population, and are well distributed along the genome. The offspring of any heterozygous parent can be segregated on the basis of marker information, to distinguish the marker haplotype inherited from each parent in a particular genomic region. Microsatellite genotyping was and is expensive and consequently many experiments lacked sufficient power to characterize regions well, and therefore detected only the largest effects [9]. Relatively few QTL were found that were useful for beef cattle improvement [10], although many interesting scientific discoveries arose from these endeavors.

Following the sequencing of the bovine genome, which led to the discovery of millions of bi-allelic SNP, and the creation of subsets of SNP that can characterize the genome and be multiplexed for cheap and efficient genotyping [11], molecular-based studies to predict animal merit have been based on high-density SNP genotypes. This review documents the current status of whole-genome prediction of breeding merit in beef cattle and describes its implementation for the purposes of selection.

Breeding objective

The breeding objective comprises a list of traits that influence the breeding goal, along with their relative emphasis [12]. An ideal breeding objective would include all the traits that will in the future influence the breeding goal. A profit-based goal would motivate the list to include all attributes that will influence income or costs. For beef cattle, these clearly include: traits that influence productivity such as reproductive performance, growth rate and survival; traits that influence cost of production such as feed intake; and traits that influence product quality such as tenderness and taste. In recent times, the list of traits has been expanding to include attributes that have been externalities. These include traits that impact the long-term contribution of beef consumption on human healthfulness, such as factors that influence anemia, cancer, obesity, diabetes and heart disease; traits that influence the environment in its broadest context, comprising air quality, water quality, soil degradation, visual farm/feedlot appearance and competition with wildlife throughout the production, finishing and processing system; and welfare factors, both of the animals in terms of exhibiting natural behaviors and being free of disease, suffering, and mortality, and of the labor in terms of worker safety. In this context, the design of a beef cattle improvement program should holistically consider traits that influence production efficiency such as individual animal measures of inputs and outputs, traits that influence the quality of the eating experience, traits that influence animal health, and traits that influence the human healthfulness of the consumed beef.

The tools available to the animal breeder to improve consumer satisfaction from beef include: the choice of breed, the choice of mating plan to exploit complementarities and heterosis, and selection for within-breed improvement [12]. The main tools for selection for within-breed improvement are the estimated breeding values (EBV) and corresponding indexes that arise from national cattle evaluations (NCE), which are available in many countries and empowers genetic improvement within the seedstock sector [4]. In the absence of genotype-environment interactions that can occur when seedstock animals are managed in different and typically superior environments compared to those of commercial animals [13], those gains are passed on to the commercial cow-calf sector by the sale of improved bulls (or semen) to be used as sires.

The current focus of the use of genetic markers for genomic prediction is to improve within-breed selection, by increasing the accuracy of existing EBV by the time the selection candidate reaches puberty, or by providing new EBV for attributes that influence the breeding goal but have not been available from conventional performance recording. Other genomic analyses that will not be considered in this review include correct assignment of parents, identification of genetic diseases, detection of signatures of selection, prediction of breed composition of crossbred animals and identification of QTL.

Estimated breeding values from national cattle evaluations in the United States

National cattle evaluations (NCE) in beef cattle began with measures of weight traits, and now include birth, weaning and yearling weights, and to a lesser extent mature weights. Rather than reporting EBV, US breed associations typically report Expected Progeny Differences (EPD), that are one-half the EBV. A summary of the traits for which EPD are typically reported is in Table 1 for the 16 most prominent US beef cattle breeds. Calving ease has been added to most national evaluation systems and, like weaning weight, includes EPD that reflect direct and maternal contributions [14]. Carcass traits have typically been problematic to collect in seedstock herds, so most carcass information tends to come from ultrasound measures of rib-eye area (REA), intramuscular fat (IMF) and fat depth [4]. Not all breed associations provide carcass EPD. Eating quality is principally limited to tenderness, but this is difficult to measure in most processing plants. In the US, carcass marbling has been used as a surrogate for tenderness/eating quality. More recently, QTL in the region of the calpain and calpastatin genes have been exploited for marker-assisted selection, using SNP that vary among breeds, most notably between Bos indicus and Bos taurus breeds. Reproductive measures have been difficult to evaluate since most breed associations have not used inventory recording systems until relatively recently, so it is impossible to determine if a female not represented as a dam actually calved or not [5]. Reproductive EPD have therefore been limited to scrotal circumference, and more recently, heifer pregnancy. There are no routine measures of input traits on a significant scale, as feed intake is problematic to measure, especially in grazing circumstances. Maintenance energy requirements have been predicted from knowledge on mature weight, condition score and milk production potential [3].

Table 1 Traits reported in national cattle evaluation for the 16 most prominent beef cattle breeds in the US

Genomic prediction

The concept of using high-density SNP genotypes to predict genetic merit was popularized by the landmark publication of Meuwissen et al. [15]. Their approach involved the computation of EBV for individual chromosome fragments, characterized by SNP genotypes or haplotypes. Estimated breeding values of selection candidates are subsequently obtained by summing up the values of all inherited chromosome fragments. This estimate is referred to as a molecular breeding value (MBV). A variety of methods has been proposed to derive EBV of chromosome fragments [16], and these can be broadly categorized into methods that fit all SNP, and methods that use mixture models that assume that not all but a fraction of the SNP have effects on the trait. All methods can be reparameterized in terms of equations that fit animal genetic effects rather than SNP effects and obtain the MBV directly, using the inverse of a genomic-based rather than a pedigree-based relationship matrix in the mixed model equations [17]. The concept of genomic prediction using a genotype-based relationship matrix predates [15] by several years [18]. In practice, so-called genomic training populations that are used to derive prediction equations, may be of inadequate size for reliable prediction of all but the largest chromosome fragments [19], leading to predictions that account for just a fraction of the additive genetic variance [20]. In this circumstance, blending the MBV and the conventional PA will improve accuracy [21]. Given the genotypes, blending can be achieved in the same analysis as the genomic training, using an inverse relationship matrix constructed from pedigree information on non-genotyped individuals and genomic information on genotyped animals [22, 23]. In the absence of the genotypes, the blending can be achieved using MBV as a correlated trait [24]. That approach requires knowledge of the covariance components relating the MBV to the trait, typically represented in publications as the genetic correlation [25, 26].

Whereas microsatellite marker studies have typically failed to identify QTL and subsequently SNP that could apply equally well across a range of breeds, there was hope that the reduced cost and the increased density of multiplexed SNP panels would lead to discoveries that could be exploited across breeds. The reduced cost per genotype for panels of 50,000 or more multiplexed SNP compared to microsatellite markers allows for more animals to be used in analyses, increasing power. In both conventional QTL studies and in genomic prediction, detection of effects relies on an association between the segregating marker genotype and the segregating causal polymorphism. The strength of this association reflects the extent of linkage disequilibrium (LD), which can be represented by the squared correlation between genotypes at two loci. Microsatellite studies exploited linkage relationships to create LD between the flanking sparse markers and a QTL within families, even when the marker was in linkage equilibrium with the QTL from a population perspective. Genomic prediction does not require family structures but takes advantage of the higher density of SNP markers and the fact that physically close loci tend to have higher LD than distant loci. Provided the genome is saturated with SNP markers, any QTL should be near some genotyped SNP and hopefully at least one will be in sufficient LD with the QTL.

Research studies of genomic prediction in livestock populations began with the release by Illumina of a high-density bovine panel of some 54,001 SNP markers [27]. In any particular breed, a proportion of these SNP will not be segregating, so the genotypes will be described in this paper as coming from a 50k panel.

Beef cattle training populations

Training involves statistical analyses that exploit individuals with both high-density genotypes and recorded performance [28]. The amount of data required for training depends upon a number of factors, including the heritability of the trait [29]. One approach to training is to use sires whose genetic merit can be assessed more reliably using progeny performance than would be the case using only measurements on the individual sire itself [9]. This may be more problematic in beef cattle than dairy cattle, as the recorded population of even the largest beef cattle breed is much smaller than that of the Holstein breed. Further, artificial insemination (AI) is much less used in beef cattle seedstock herds than in dairy herds, collectively resulting in fewer highly reliable sires available for use in training.

Industry populations have advantages for genomic prediction. In the case of elite or widely used industry animals, the individuals included in the training data will be relevant to the commercial population. For AI sires, DNA is readily accessible despite the disparate ownership or physical location of the animals. The principal source of performance information comes as EPD from NCE and is well represented for growth traits, moderately well for ultrasound traits, poorly for behavior, reproduction and longevity traits, and typically with no information on many other traits such as disease resistance or eating quality. Since most recorded animals are purebred, training on crossbred data is seldom an option using NCE data and is limited to those few breed associations that collect crossbred data.

A US repository of DNA from over 3,000 Angus bulls born since 1948 was assembled by the University of Missouri [30]. These bulls are represented in American Angus Association pedigrees and have generally been widely used. Accordingly, these bulls have EPD and accuracies for production traits: calving ease (direct); birth weight; weaning weight; yearling weight; yearling height; scrotal circumference; maternal traits: maternal calving ease; milk; mature weight; mature height; carcass traits: carcass weight; marbling; rib eye area; fat depth; along with some newly released trait EPD: docility; and heifer pregnancy. The accuracies of EPD on old bulls are limited for some traits. Igenity, a genomic testing service owned by the animal health company Merial, has used the results from the analysis of this Angus population, along with other resource populations, to market a reduced panel comprised of a subset of informative SNP referred to as a 50k-derived product. It is marketed in the US in conjunction with the American Angus Association and costs $65 [31].

The US Meat Animal Research Center (US-MARC) at Clay Center Nebraska has worked with some breed associations to develop a repository of some 2,026 influential or upcoming bulls in 16 of the most prominent beef breeds in the US with EPD from NCE and includes: Angus, Beefmaster, Brahman, Brangus, Braunvieh, Charolais, Chiangus, Gelbvieh, Hereford, Limousin, Maine-Anjou, Red Angus, Salers, Santa Gertrudis, Shorthorn, and Simmental. Initial plans for the use of this repository were to provide genomic predictions of these bulls from training analyses based on a US-MARC crossbred population [32] and to carry out multi-breed training. These SNP genotypes have now been made available to the respective breed associations.

The alternative to training on widely-used sires is to train using phenotypes collected specifically for genomic analyses. This could be achieved using non seedstock field data, but in many cases the mating designs and contemporary group classifications are not entirely adequate for the purpose. Most field data comprise offspring from natural mating, so sires tend to be nested within rather than cross-classified by contemporary groups. In the case of carcass traits, animals tend to have their ownership transferred several times between weaning and harvest, making it difficult to ensure harvest cohorts were managed together throughout their entire lifetime. For reproductive traits, it is difficult to obtain sizeable cohorts of animals for comparison, particularly for phenotypic measurements obtained after first calving, as birth cohorts get subdivided according to sex of calf, age of dam, and whether or not yearlings became pregnant. These problems can be overcome by sourcing animals from large herds and by designing the study prior to the birth of the study animals, which may be several years prior to the collection of phenotypes.

The US carcass merit project (CMP) was one such long-term industry-funded semi-structured undertaking initiated in 1998 that collected carcass data, tenderness and sensory attributes on over 8,200 progeny. Some of the half-sib offspring of more than 70 sires across 13 breeds were DNA sampled. The sires were widely-used AI bulls from various breeds and dams were commercial cows [33]. The dataset has been valuable to validate early genomic tests being commercialized in the US. Validation of tests using these data has been undertaken by the National Beef Cattle Evaluation Consortium (NBCEC) and the details having been published on-line by Van Eenennaam et al. [34]. More recently, the CMP dataset has been genotyped using high-density SNP chips by at least two different organizations to identify genes and to apply whole-genome prediction, which will prevent this resource from being used for independent validation of future tests derived from that data.

Collecting data for more novel phenotypes requires the deliberate generation of suitable populations. Given the current dominant market position of the Angus breed in the US, it was an obvious candidate for any new studies to expand the scope of traits for genomic prediction. Two large studies have been undertaken, one at Iowa State University to investigate fatty acid and mineral content in beef as possible targets for improving the human healthfulness of beef, and another at Colorado State University to investigate feedlot health. The healthfulness study involved several cohorts representing 2,300 predominately Angus cattle assessed for carcass and meat quality attributes, including tenderness and sensory information, in addition to extensive phenotyping of traits that might influence the human healthfulness of beef. These healthy beef traits include mineral and fatty acid compositions of key muscles [35]. The feedlot health study used two annual crops of about 1,500 composite British and Continental steers from one ranch in Nebraska. The animals were extensively phenotyped for feedlot health, particularly respiratory disease and response to treatment. Sickness was assessed visually, by temperature profiles and by lung damage scores. Data includes temperament and immunological measures [36]. Both experiments included body weight and a number of carcass and meat quality phenotypes. These collective resources have been used, along with other populations, to develop an Angus 50k product for production and carcass traits that Pfizer Animal Genetics has marketed in the US for $124-$139, depending upon the number of animals tested [37], with predictions from this panel now incorporated in NCE undertaken for the American Angus Association.

Research herds with deep phenotyping are also candidates for studies of genomic prediction. The most comprehensive such resource is represented by the US-MARC germplasm evaluation studies, the recent cohorts being known as the Cycle VII and F1-squared populations. In addition to an across-breed training analysis for which single-SNP effects have been published for birth, weaning and yearling weights and their respective gains [38], this population was used to develop a low-density 196-SNP panel with markers believed to be informative for weaning weight. Such reduced panels comprised of only the most informative markers were believed to be more cost-effective and therefore more likely to be widely adopted by the beef industry. That panel was used in a project coordinated by the NBCEC to demonstrate the use of reduced panels in seedstock herds, and the incorporation of the resulting MBV into NCE [39].

The collection of feed intake on large numbers of animals is still problematic from a practical viewpoint, and to date, such data has been limited to measuring relatively small disparate groups of animals during finishing, with findings focused on QTL detection rather than genomic prediction. Other datasets of limited size have been collected on a range of traits, including reproductive performance and tick resistance but have not yet had any findings published from a genomic prediction perspective.

Funding for genotyping training populations

Costs for conventional pedigree and performance recording and for NCE have been met by producer funds in the US. Public funds have been used for the development of NCE methodology. Public funds were not immediately available for extensive genotyping of training populations, and neither seedstock breeders nor breed associations had funds to adopt this technology beforehand given the uncertain nature of its value. Fortunately, applications of this approach in beef cattle improvement were considered as business opportunities by commercial companies such as Merial Igenity and Pfizer Animal Genetics to invest in the training phase, presumably with expectations of recouping returns on that investment through future sale of genomic tests. However, this situation has changed industry dynamics, introducing competitive partners into the process of ranking animals, and has increased the proprietary nature of performance information, genotypes and analytical approaches. This is one reason for the dearth of refereed publications on the accuracy of genomic prediction in beef cattle, in contrast to the dairy cattle situation.

Predictive ability of whole-genome findings

Confidence in genomic predictions can only be provided by validation in a group of animals that are not included in the training population. Close relationships between animals in training and validation populations tend to lead to better predictive ability than when the groups are more distantly related [40]. Analysis of simulated data suggests that methods based on mixture models provide better predictive ability than methods that assume all the SNP have predictive value [15], while analysis of field data tends to demonstrate relatively little difference between alternative methods, and some inconsistencies appear from trait to trait as to which is the most predictive method [41, 42]. There appears to be more variation in predictive ability according to the choice of validation population than there is between methods.

Within-breed 50k predictions

One of the few reports on accuracy of genomic predictions in beef cattle analysed deregressed EPD [43] from NCE to quantify cross-validation results from 2,100 Angus AI bulls [44]. The data were partitioned into three subsets, with training animals in two groups and validation animals in the third. Subsets were created so that no sire had sons in both the training and validation groups. Genomic predictions were obtained from the training data using method Bayes C [41]. Predictive ability was quantified as correlations between 50k predictions and realized (deregressed) performance (Table 2). The general conclusion is that correlations between genomic predictions from 50k SNP and deregressed EPD in independent datasets of related animals are 0.5-0.7. It is not possible from these correlations to readily derive the genetic correlation between genomic prediction and the true BV, because of heterogeneity of variance among the deregressed EPD. This heterogeneity does not impact the expectation of the estimated covariance between genomic predictions and deregressed EPD, but it does impact the estimated variance of the deregressed EPD. Furthermore, the genotyped animals represent AI sires, and these represent highly selected individuals, so their genetic variance is not likely to be representative of the population genetic variance. Also, correlations between genomic prediction and EPD do not provide expectation on the genetic correlation, due to the varying degrees of shrinkage influencing EPD, which vary in their information content. Accordingly, correlations between genomic prediction and EPD or deregressed EPD provide a guide to accuracy, but cannot be interpreted as quantifying the proportion of variation accounted for by the genomic prediction applied to new animals. This would not be the case for correlations between genomic prediction and homogeneous information such as individual phenotypic observations.

Table 2 Correlations of 50k or 600 SNP predictions with deregressed EPD for various traits using cross-validation with three subsets of the data

Other numerically important breeds tend to have fewer registrations than Angus and it will be difficult to collect comparable sized training populations of AI sires. In contrast to the dairy industry, most bulls are used solely in commercial herds that do not record parentage or individual performance and therefore do not obtain progeny information for training or validation. The American Hereford Association has increased the 50k genotypes provided by US-MARC to develop a training population of 800 animals, but no results have been published yet. The other US breeds have even fewer animals ready for training.

Genomic prediction for beef cattle healthfulness has shown varying levels of predictive ability, as determined by the proportion of variation accounted for by markers [35]. Using samples from the Longissimus dorsi, iron concentration of beef could be readily predicted, whereas magnesium, manganese, phosphorus and zinc concentrations appeared to be under less genetic control. For other minerals such as calcium, copper, potassium and sodium, concentrations could not be predicted. Prediction of the fatty acid's concentrations showed similar trends to that of the minerals' concentration. For the predominant even-numbered saturated fatty acids C14:0, C16:0 and C18:0, monounsaturated C18:1 and polyunsaturated C18:2, prediction was good, while for C18:3 and conjugated linoleic acid (CLA) concentrations, predictions were not conclusive. These results look promising to develop tools capable of modifying the concentration of saturated fatty acids, or the relative proportions of saturated and unsaturated fatty acids. For these traits, the challenge will consist in developing a market for beef with modified fatty acid composition.

Using the same dataset as for beef healthfulness, it has been shown that carcass and beef quality traits can be predicted [35]. Hot carcass weight, calculated yield grade, marbling score and fat thickness had 40-50% of phenotypic variance explained by the 50k markers, whereas markers accounted for less than 30% of the variation for dressing percentage, loin eye area and tenderness assessed by Warner-Bratzler shear force. Cross validation results were not reported.

Within-breed reduced panels

Reduced SNP panels can be produced either to be highly informative for a particular trait or for several traits by including the most strongly associated SNP, or to be informative for high-density genomic prediction after imputing the high-density panel from a reduced set of evenly spaced SNP with high minor allele frequency [45]. To date, the beef industry focus has been on subsets of markers chosen to be informative for a subset of traits that are believed to have the most economic relevance and greatest market opportunity.

Mixture models such as Bayes B and Bayes C [41] assume that some fraction of the SNP have zero effect on the trait. The posterior frequency with which any particular SNP was fitted in an MCMC analysis reflects the informativeness of particular SNP and can be used for SNP selection. Subsets of 600 SNP markers created by selecting the 20 markers on each bovine chromosome with the highest model frequency, from Bayes C analyses with 90% of 50k SNP assumed to have zero effect, demonstrated relatively little loss of predictive ability compared to 50k predictions [43]. Cheaper genotyping can be achieved by reducing the number of markers to a single set of 384 SNP, chosen for predictive ability across the portfolio of traits of interest. However, reducing the number of SNP below 600 reduces predictive ability. For example, the correlation reported in [43] for sets of the best 50, 100, 150 or 200 SNP chosen to predict marbling in Angus were 0.28, 0.29, 0.39, and 0.43, well below the 0.67 achieved with 600 SNP. A single set of 384 markers chosen from the above analysis for predictive ability across a range of traits was validated in a new population of 275 Angus bulls [43]. The correlations from that analysis were 0.59 for marbling, 0.32 for backfat, 0.58 for rib eye, 0.44 for carcass weight, 0.39 for heifer pregnancy and 0.35 for yearling weight.

In the study on beef healthfulness [35], subsets of as few as 10 markers retained more than half of the predictive ability of the 50k SNP chip when used to predict the even-numbered saturated fatty acids C14:0 and C16:0. The genomic architecture of mineral and fatty acid concentrations is likely to be much simpler, as the biochemical pathways and enzymes involved in metabolizing and catabolizing these compounds have been identified and seem to be somewhat straightforward, in contrast to traits such as growth rate, which are the collective result of genes influencing bone growth, muscle growth, fat accumulation, visceral weight among other factors.

The development of reduced panels for any quantitative trait in breeds other than Angus is currently limited by the lack of training populations. In contrast to the dairy industry, where reduced panels are being used for imputation of 50k markers for genomic prediction [46], target populations in beef cattle are diverse in terms of species (Bos indicus and Bos taurus) and breeds. Furthermore, many pre-pubertal selection candidates are offspring of natural mating rather than of AI sires. Collectively, these facts increase the genetic distance between the training and target populations.

Across-breed panels

Prediction across breeds is more problematic because different breeds may exhibit different QTL, dominance or epistasis can occur, and allele frequencies may vary between populations. Linkage disequilibrium (LD) is not very consistent across breeds and therefore training in one beef cattle breed using 50k genotypes will not be very effective to predict a different breed [47]. Simulated data using actual 50k genotypes from the CMP and an Angus dataset as if they were causal genes and adding a random environmental effect to represent a trait with 50% heritability, demonstrated that predictive ability varied according to the number of simulated QTL. The best results were achieved for the smallest number of QTL, since in that scenario the average size of the QTL was larger than when more QTL were simulated. The across-breed predicted correlation from the simulation [47] varied from a high of 0.4 for 50 QTL down to 0.2-0.3 for 500 QTL. These correlations account for up to 18% of genetic variance for 50 genes and less than 10% of variance for 500 genes. Unpublished data predicting the merit of Hereford bulls using training results from Angus bulls always resulted in positive correlations, but typically less than 0.10, with the best correlation being 0.18 for birth weight and slightly less for yearling weight. Genomic prediction in beef cattle based only on 50k genotypes will therefore require training individuals from every target breed, confirming findings from simulations [48].

Recently released next generation Illumina HD or Affymetrix Bos-1 panels, with more than a 10-fold increase in SNP density beyond the 50k, will allow imputation of missing SNP genotypes in animals already genotyped for 50k panels [45, 46]. It is hoped that the 10-fold increased SNP density will improve across-breed prediction, avoiding the need for large training populations of every target breed, but this has yet to be demonstrated in practice.

Genomic prediction across-breed using reduced panels will be inferior to 50k based predictions. A subset of 192 SNP markers was chosen from the US-MARC association analysis for weaning weight reported in [38] and applied to predict merit for weaning weight and post-weaning gain in purebred calves representing seven of the breeds represented as crossbreds in the US-MARC training data. The genetic correlation estimated between the MBV and direct effects for weaning weight was slightly negative (-0.05) in one breed, 0.0 in another, and ranged from 0.10-0.28 in the remaining breeds [39]. These results are disappointingly low.

Incorporation of genomic information in US national cattle evaluation

Both predictions from Merial Igenity and Pfizer Animal Genetics are currently used in the American Angus Association (AAA) NCE by including them as correlated traits. The estimated genetic correlations for the Merial Igenity MBV are 0.54 for carcass weight, 0.58 for REA, 0.50 for fat and 0.65 for marbling [25]. Corresponding values have not yet been reported for the Pfizer Animal Genetics MBV. Procedurally, breeders send DNA samples to AAA, where they are anonymously recoded and forwarded to the relevant genomics company. The MBV are reported back to AAA to be provided to the breeders and included in NCE. In this circumstance, retraining to improve the accuracy of genomic prediction is not an option as no party has access to both the genotypes and EPD or phenotypic performance of the genotyped individuals.

Future hopes

Predictive ability is influenced by effective population size, heritability, and the number of animals in the training data, among other factors [20, 29]. Increasing the number of genotyped animals should increase predictive ability. Ideally, the training data should accumulate as the seedstock producers genotype individuals for selection purposes. Unlike for the dairy industry, this is not occurring yet in the beef industry, since genomics companies are marketing predictions without the genotypes going into the national databases administered by the breed associations. Research populations may therefore be critical to the accumulation of training animals in the near term. In Australia, industry has actively promoted an information nucleus for this very purpose [49]. The presence of such populations will inevitably place strain on the relationship between genomics companies that want to keep information of a proprietary nature and public/industry funding efforts. Pooling training populations across countries provides an opportunity to increase training data size, but may add complications. Different countries sometimes define traits in different ways (e.g. age-adjusted or weight-adjusted), and have different harvest end-points (e.g. weight-constant or fat-constant), resulting in imperfect relationships between the traits in different countries. Further, genotype by environment interactions can also be important because production conditions tend to be more diverse in beef cattle than in dairy production. Pooling training data across breeds provides an appealing alternative to increase predictive power but will require the use or imputation of new higher-density SNP panels. The use of haplotypes [50] may also provide additional power, although this has yet to be demonstrated in beef cattle with field rather than simulated data. Cost-effective use of the technology will likely result in approaches that exploit genotype imputation, and use mixed densities of genotyping on individual animals. This will likely include the DNA sequencing of individual animals [51], such as widely-used AI sires, and the imputation of sequences. However, additional SNP information alone may reduce predictive ability [47] unless the size of the training populations increases. Exploiting bioinformatics, such as from expression analyses and knowledge of the location of genes known to influence traits in beef cattle or other species, may help to increase predictive ability by allowing focusing on additional SNP only in the regions that lack sufficient LD. New analytical methods, such as approaches that explicitly fit QTL effects [52] rather than SNP effects (such as methods that jointly account for LD and linkage information [53]) may also help.

Extension of genomic predictions to the full range of traits that influence consumer satisfaction will further require a focus on the collection of reliable phenotypic information across the broad spectrum of traits. Collecting such information will likely rely on public funding efforts, but even then will be limited by the availability of meaningful phenotypes for some traits. New electronic technologies that facilitate the collection of phenotypes on large cohorts will also be invaluable.

Conclusion

Genomic prediction offers accuracies that exceed those of pedigree-based parent average of young selection candidates. The highest accuracies are achieved for offspring of the training population. Accuracies can be equivalent to progeny tests based on up to 10 or so offspring, providing a slightly higher predictive ability than a single phenotypic observation on the individual. These accuracies are not yet sufficiently high to warrant selection in the absence of phenotypic information, particularly as these accuracies tend to erode when assessed in validation populations that are more distant from the training population in terms of the number of meioses separating generations. Accuracies are expected to improve with further research, as the training population grows in terms of numbers of genotyped animals, and density of SNP genotypes per animal.

Phenotyping is now the principal limitation in expanding the series of traits beyond those routinely recorded for NCE. In the meantime, applying genomic prediction will influence traits that were easy to record in conventional improvement programs, rather than addressing the traits difficult and costly to measure.

Sharing of information among parties to the benefit of industry is still in its infancy, as is the incorporation of MBV into NCE. The latter activity will cause particular challenges for small breed associations which lack the funding or expertise to change their NCE systems. Whereas it had been hoped that genomic prediction would facilitate selection in small breed associations with fewer registered animals, the current need for within-breed training will serve only to increase the technology gap between the breeds and facilitate faster rates of change in those breeds that have a large market share.

Abbreviations

CMP:

(carcass merit project)

EBV:

(estimated breeding value)

EPD:

(expected progeny difference)

LD:

(linkage disequilibrium)

IMF:

(intramuscular fat)

MBV:

(molecular breeding value)

NBCEC:

(National Beef Cattle Evaluation Consortium)

NCBA:

(National Cattlemen's Beef Association)

NCE:

(national cattle evaluation)

PA:

(parent average)

QTL:

(quantitative trait locus)

REA:

(rib-eye area)

SNP:

(single nucleotide polymorphism)

US-MARC:

(United States Meat Animal Research Center).

References

  1. American Angus Association® - The Business Breed. (Accessed 18 February 2011), [http://www.angus.org/nce/genetictrends.aspx]

  2. Red Angus EDP Trends. (Accessed 18 February 2011), [http://redangus.org/genetics/epd-trends]

  3. Williams JL, Garrick DJ, Speidel SE: Reducing bias in maintenance energy expected progeny difference by accounting for selection on weaning and yearling weights. J Anim Sci. 2009, 87: 1628-1637. 10.2527/jas.2008-1442.

    Article  CAS  PubMed  Google Scholar 

  4. Golden BL, Garrick DJ, Benyshek LL: Milestones in beef cattle genetic evaluation. J Anim Sci. 2009, 87: E3-E10. 10.2527/jas.2008-1430.

    Article  CAS  PubMed  Google Scholar 

  5. Garrick DJ, Golden BL: Producing and using genetic evaluations in the United States beef industry of today. J Anim Sci. 2009, 87: E11-E18. 10.2527/jas.2008-1431.

    Article  CAS  PubMed  Google Scholar 

  6. VanRaden PM, Wiggans GR: Derivation, calculation, and use of national animal model information. J Dairy Sci. 1991, 74: 2737-2746. 10.3168/jds.S0022-0302(91)78453-1.

    Article  CAS  PubMed  Google Scholar 

  7. Hansen MP, Law GRJ, Van Zandt JN: Differences in susceptibility to Marek's disease in chickens carrying two different B locus blood group alleles. Poult Sci. 1967, 46: 1268-1268.

    Google Scholar 

  8. Georges M, Nielsen D, Mackinnon M, Mishra A, Okimoto R, Pasquino AT, Sargeant LS, Sorensen A, Steele MR, Zhao X, Womack JE, Hoeschele I: Mapping quantitative trait loci controlling milk production in dairy cattle by exploiting progeny testing. Genetics. 1995, 139: 907-920.

    PubMed Central  CAS  PubMed  Google Scholar 

  9. Weller JI, Kashi Y, Soller M: Power of daughter and granddaughter designs for determining linkage between marker loci and quantitative trait loci in dairy cattle. J Dairy Sci. 1990, 73: 2525-2537. 10.3168/jds.S0022-0302(90)78938-2.

    Article  CAS  PubMed  Google Scholar 

  10. Dekkers JCM: Commercial application of marker- and gene-assisted selection in livestock: strategies and lessons. J Anim Sci. 2004, 82: E313-E328.

    PubMed  Google Scholar 

  11. Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, C'Connell J, Moore SS, Smith TPL, Sonstegard TS, Van Tassell CP: Development and characterization of a high density SNP genotyping assay for cattle. PLoS ONE. 2009, 4: e5350-10.1371/journal.pone.0005350.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Harris DL, Stewart TS, Arboleda CR: Animal breeding programs: a systematic approach to their design. Advances in agricultural technology. 1984, Agricultural Research Service, USDA-ARS Bull, 1-14.

    Google Scholar 

  13. Morris CA, Baker RL, Hickey SM, Johnson DL, Cullen NG, Wilson JA: Evidence of genotype by environment interaction for reproductive and maternal traits in beef cattle. Anim Prod. 1993, 56: 69-83. 10.1017/S0003356100006176.

    Article  Google Scholar 

  14. Willham RL: The covariance between relatives for characters composed of components contributed by related individuals. Biometrics. 1963, 19: 18-10.2307/2527570.

    Article  Google Scholar 

  15. Meuwissen THE, Hayes B, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001, 157: 1819-1829.

    PubMed Central  CAS  PubMed  Google Scholar 

  16. Gianola D, de los Campos G, Hill WG, Manfredi E, Fernando R: Additive genetic variability and the Bayesian alphabet. Genetics. 2009, 183: 347-363. 10.1534/genetics.109.103952.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Strandén I, Garrick DJ: Technical note: Derivation of equivalent computing algorithms for genomic predictions and reliabilities of animal merit. J Dairy Sci. 2009, 92: 2971-2975. 10.3168/jds.2008-1929.

    Article  PubMed  Google Scholar 

  18. Nejati-Javaremi A, Smith C, Gibson JP: Effect of total allelic relationship on accuracy of evaluation and response to selection. J Anim Sci. 1997, 75: 1738-1745.

    CAS  PubMed  Google Scholar 

  19. MacLeod IM, Hayes BJ, Savin KW, Chamberlain AJ, McPartlan HC, Goddard ME: Power of a genome scan to detect and locate quantitative trait loci in cattle using dense single nucleotide polymorphisms. J Anim Breed Genet. 2009, 127: 133-142.

    Article  Google Scholar 

  20. Goddard ME: Genomic selection: prediction of accuracy and maximization of long term response. Genetica. 2009, 136: 245-257. 10.1007/s10709-008-9308-0.

    Article  PubMed  Google Scholar 

  21. VanRaden PM, Van Tassell CP, Wiggans GR, Sonstegard TS, Schnabel RD, Taylor JF, Schenkel FS: Invited review: Reliability of genomic predictions for North American Holstein bulls. J Dairy Sci. 2009, 92: 16-24. 10.3168/jds.2008-1514.

    Article  CAS  PubMed  Google Scholar 

  22. Legarra A, Aguilar I, Misztal I: A relationship matrix including full pedigree and genomic information. J Dairy Sci. 2009, 92: 4656-4663. 10.3168/jds.2009-2061.

    Article  CAS  PubMed  Google Scholar 

  23. Aguilar I, Misztal I, Johnson DL, Legarra A, Tsuruta S, Lawlor TJ: Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J Dairy Sci. 2010, 93: 743-752. 10.3168/jds.2009-2730.

    Article  CAS  PubMed  Google Scholar 

  24. MacNeil MD, Nkrumah JD, Woodward BW, Northcutt SL: Genetic evaluation of Angus cattle for carcass marbling using ultrasound and genomic indicators. J Anim Sci. 2010, 88: 517-522. 10.2527/jas.2009-2022.

    Article  CAS  PubMed  Google Scholar 

  25. MacNeil MD, Northcutt SL, Schnabel RD, Garrick DJ, Woodward BW, Taylor JF: Genetic correlations between carcass traits and molecular breeding values in Angus cattle. Proceedings of the Ninth World Congress on Genetics applied to Livestock Production: 1-6 August 2010; Leipzig. 2010, 482-

    Google Scholar 

  26. Hanford K, Thallman M, Kachman S, Kuehn L, Quaas R, Tempelman R, Fernando R, Pollak J: Estimation of the proportion of variation accounted for by DNA tests I. Genetic variance. Proceedings of the Ninth World Congress on Genetics applied to Livestock Production: 1-6 August 2010; Leipzig. 2010, 822-

    Google Scholar 

  27. Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O'Connell J, Moore SS, Smith TPL, Sonstegard TS, Van Tassell CP: Development and characterization of a high-density SNP genotyping assay for cattle. PLoS ONE. 2009, 4: e5350-10.1371/journal.pone.0005350.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Meuwissen T: Genomic selection: Marker-assisted selection on a genome wide scale. J Anim Breed Genet. 2007, 124: 321-322. 10.1111/j.1439-0388.2007.00708.x.

    Article  PubMed  Google Scholar 

  29. Daetwyler HD, Villanueva B, Woolliams JA: Accuracy of predicting the genetic risk of disease using a genome-wide approach. PLoS ONE. 2008, 3: e3395-10.1371/journal.pone.0003395.

    Article  PubMed Central  PubMed  Google Scholar 

  30. McClure MC, Morsci N, Schnabel RD, Kim JW, Yao P, Rolf MM, McKay SD, Gregg SJ, Chapple RH, Northcutt SL, Taylor JF: A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus. Anim Genet. 2010, 41: 597-607. 10.1111/j.1365-2052.2010.02063.x.

    Article  CAS  PubMed  Google Scholar 

  31. Igenity® profile for Angus goes HD. (Accessed 25 November 2010), [http://www.angusbeefbulletin.com/extra/2010/02feb10/0210mg_agi_merial.html]

  32. Thallman RM: Whole genome selection. (Accessed 25 November 2010), [http://animalscience.ucdavis.edu/animalbiotech/Outreach/Whole_Genome_Selection.pdf]

  33. Thallman RM, Moser DW, Dressler EW, Totir LR, Fernando RL, Kachman SD, Rumph JM, Dikeman ME, Pollak EJ: Carcass merit project: DNA marker validation. (Accessed 25 November 2010), [http://www.beefimprovement.org/proceedings/genetic-prediction-workshop/GPW-CarcassMeritProject-Final.pdf]

  34. National beef cattle evaluation consortium: commercial genetic test validation. (Accessed 27 April 2011), [http://www.nbcec.org/validation.html]

  35. Reecy JM, Tait RG, VanOverbeke DL, Garmyn AJ, Mateescu RG, Van Eenennaam AL, Duan Q, Liu Q, Schoonmaker JP, Drewnoski ME, Beitz DC, Kizilkaya K, Fernando RL, Garrick DJ: Use of genomics to improve healthfulness and quality of meat. Proceedings of the Ninth World Congress on Genetics applied to Livestock Production: 1-6 August 2010; Leipzig. 2010, 53-

    Google Scholar 

  36. Brigham BW, McCallister CM, Enns RM: Genetics of disease resistance. (Accessed 25 November 2010), [http://www.rangebeefcow.com/2009/documents/BrighamEnns2009RBCS_pp.pdf]

  37. Cooperative resources international: Genex to market Pfizer HD 50k. (Accessed 25 November 2010), [http://cri.crinet.com/news1001/GenexToMarketPfizerHD50K]

  38. Snelling WM, Allan MF, Keele JW, Kuehn LA, McDaneld T, Smith TPL, Sonstegard TS, Thallman RM, Bennett GL: Genome-wide association study of growth in crossbred beef cattle. J Anim Sci. 2010, 88: 837-848. 10.2527/jas.2009-2257.

    Article  CAS  PubMed  Google Scholar 

  39. Spangler ML, Kachman S, Hanford K, Bennett G, Snelling W, Kuehn L, Pollak J: Integration of DNA marker information into breeding value predictions. Nebraska Beef Cattle Report. 2011, Lincoln, 11-12. (Accessed 2 March 2011), [http://www.nlc.state.ne.us/epubs/U2300/T001-2011.pdf]

    Google Scholar 

  40. Habier D, Tetens J, Seefried FR, Lichtner P, Thaller G: The impact of genetic relationship information on genomic breeding values in German Holstein cattle. Genet Sel Evol. 2010, 42: 5-10.1186/1297-9686-42-5.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Habier D, Fernando RL, Kizilkaya K, Garrick DJ: Extension of the Bayesian alphabet for genomic selection. Proceedings of the Ninth World Congress on Genetics applied to Livestock Production: 1-6 August 2010; Leipzig. 2010, 468-

    Google Scholar 

  42. Wolc A, Stricker C, Arango J, Settar P, Fulton JE, O'Sullivan NP, Habier D, Fernando R, Garrick DJ, Lamont SJ, Dekkers JCM: Breeding value prediction for production traits in layers using pedigree or genomic relationships in a reduced animal model. Genet Sel Evol. 2011, 43: 5-10.1186/1297-9686-43-5.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Garrick DJ, Taylor JF, Fernando RL: Deregressing estimated breeding values and weighting information for genomic regression analyses. Genet Sel Evol. 2009, 41: 55-10.1186/1297-9686-41-55.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Garrick DJ: The nature and scope of some whole genome analyses in US beef cattle. Proceedings of the Beef Improvement Federation 41st Annual Research Symposium: 30 April-3May 2009; Sacramento. 2009, 41: 92-102.

    Google Scholar 

  45. Habier D, Fernando RL, Dekkers JC: Genomic selection using low-density marker panels. Genetics. 2009, 182: 343-353. 10.1534/genetics.108.100289.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Zhang Z, Druet T: Marker imputation with low-density marker panels in Dutch Holstein cattle. J Dairy Sci. 2010, 93: 5487-5494. 10.3168/jds.2010-3501.

    Article  CAS  PubMed  Google Scholar 

  47. Kizilkaya K, Fernando RL, Garrick DJ: Genomic prediction of simulated multibreed and purebred performance using observed fifty thousand single nucleotide polymorphism genotypes. J Anim Sci. 2010, 88: 544-551. 10.2527/jas.2009-2064.

    Article  CAS  PubMed  Google Scholar 

  48. Toosi A, Fernando RL, Dekkers JCM: Genomic selection in admixed and crossbred populations. J Anim Sci. 2010, 88: 32-46. 10.2527/jas.2009-1975.

    Article  CAS  PubMed  Google Scholar 

  49. van der Werf J, Banks RG: A genomic information nucleus to accelerate rates of genetic improvement in sheep. Proceedings of the Ninth World Congress on Genetics applied to Livestock Production: 1-6 August 2010; Leipzig. 2010, 46-

    Google Scholar 

  50. Hayes BJ, Chamberlain AJ, McPartlan H, MacLeod I, Sethuraman L, Goddard ME: Accuracy of marker-assisted selection with single markers and marker haplotypes in cattle. Genet Res (Camb). 2007, 89: 215-220.

    Article  CAS  Google Scholar 

  51. Meuwissen T, Goddard ME: Accurate prediction of genetic value for complex traits by whole-genome resequencing. Genetics. 2010, 185: 623-631. 10.1534/genetics.110.116590.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Fernando RL, Grossman M: Marker-assisted selection using best linear unbiased prediction. Genet Sel Evol. 1989, 21: 467-477. 10.1186/1297-9686-21-4-467.

    Article  PubMed Central  Google Scholar 

  53. Meuwissen THE, Karlsen A, Lien S, Olsaker I, Goddard ME: Fine mapping of a quantitative trait locus for twinning rate using combined linkage and linkage disequilibrium mapping. Genetics. 2002, 161: 373-379.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

DJG is supported by the United States Department of Agriculture, National Research Initiative grant USDA-NRI-2009-03924, Agriculture and Food Research Initiative competitive grant 2009-35205-05100 from the National Institute of Food and Agriculture Animal Genome Program, and by Hatch and State of Iowa funds through the Iowa Agricultural and Home Economic Experiment Station, Ames, IA. An anonymous referee is acknowledged for providing constructive comments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Dorian J Garrick.

Additional information

Competing interests

The author declares that they have no competing interests.

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Garrick, D.J. The nature, scope and impact of genomic prediction in beef cattle in the United States. Genet Sel Evol 43, 17 (2011). https://doi.org/10.1186/1297-9686-43-17

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1297-9686-43-17

Keywords