A look under the hood of genomic-estimated breed compositions for brangus cattle: What have we learned?

The Brangus cattle were developed to utilize the superior traits of Angus and Brahman cattle. Their genetic compositions are expected to be stabilized at 3/8 Brahman and 5/8 Angus. Previous studies have shown more than expected Angus lineage with Brangus cattle, and the reasons are yet to be investigated. In this study, we revisited the breed compositions for 3,605 Brangus cattle from three perspectives: genome-wise (GBC), per chromosomes (CBC), and per chromosome segments (SBC). The former (GBC) depicted an overall picture of the "mosaic" genome of the Brangus attributable to their ancestors, whereas the latter two criteria (CBC and SBC) corresponded to local ancestral contributions. The average GBC for the 3,605 Brangus cattle were 70.2% Angus and 29.8% Brahman. The K-means clustering supported the postulation of the mixture of 1/2 Ultrablack (UB) animals in Brangus. For the non-UB Brangus animals, the average GBC were estimated to be 67.4% Angus and 32.6% Brahman. The 95% confidence intervals of their overall GBC were 60.4%-73.5% Angus and 26.5%-39.6% Brahman. Possibly, genetic selection and drifting have resulted in an approximately 5% average deviation toward Angus lineage. The estimated ancestral contributions by chromosomes were heavily distributed toward Angus, with 27 chromosomes having an average Angus CBC greater than 62.5% but only two chromosomes (5 and 20) having Brahman CBC greater than 37.5%. The chromosomal regions with high Angus breed proportions were prevalent, tending to form larger blocks on most chromosomes. In contrast, chromosome segments with high Brahman breed proportion were relatively few and isolated, presenting only on seven chromosomes. Hence, genomic hitchhiking effects were strong where Angus favorable alleles resided but weak where Brahman favorable alleles were present. The functions of genes identified in the chromosomal regions with high ( ≥ 75 % ) Angus compositions were diverse yet may were related to growth and body development. In contrast, the genes identified in the regions with high ( ≥ 37.5 % ) Brahman compositions were primarily responsible for disease resistance. In conclusion, we have addressed the questions concerning the Brangus genetic make-ups. The results can help form a dynamic picture of the Brangus breed formation and the genomic reshaping.

A Causality Perspective of Genomic Breed Composition for Composite Animals.

Genomic breed composition (GBC) of an individual animal refers to the partition of its genome according to the inheritance from its ancestors or ancestral breeds. For crossbred or composite animals, knowing their GBC is useful to estimate heterosis, to characterize their actual inheritance from foundation breeds, and to make management decisions for crossbreeding programs. Various statistical approaches have been proposed to estimate GBC in animals, but the interpretations of estimates have varied with these methods. In the present study, we proposed a causality interpretation of GBC based on path analysis. We applied this method to estimating GBC in two composite breeds of beef cattle, namely Brangus and Beefmaster. Three SNP panels were used to estimate GBC: a 10K SNP panel consisting of 10,226 common SNPs across three GeneSeek Genomic Profiler (GGP) bovine SNP arrays (GGP 30K, GGP 40K, and GGP 50K), and two subsets (1K and 5K) of uniformly distributed SNPs. The path analysis decomposed the relationships between the ancestors and the composite animals into direct and indirect path effects, and GBC was measured by the relative ratio of the coefficients of direct (D-GBC) and combined (C-GBC) effects from each ancestral breed to the progeny, respectively. Estimated GBC varied only slightly between different genotyping platforms and between the three SNP panels. In the Brangus cattle, because the two ancestral breeds had a very distant relationship, the estimated D-GBC and C-GBC were comparable to each other in the path analysis, and they corresponded roughly to the estimated GBC from the linear regression and the admixture model. In the Beefmaster, however, the strong relationship in allelic frequencies between Hereford and Shorthorn imposed a challenge for the linear regression and the admixture model to estimated GBC reliably. Instead, D-GBC by the path analysis included only direct ancestral effects, and it was robust to bias due to high genomic correlations between reference (ancestral) breeds.

Genomic mating as sustainable breeding for Chinese indigenous Ningxiang pigs.

An important economic reason for the loss of local breeds is that they tend to be less productive, and hence having less market value than commercial breeds. Nevertheless, local breeds often have irreplaceable values, genetically and sociologically. In the breeding programs with local breeds, it is crucial to balance the selection for genetic gain and the maintaining of genetic diversity. These two objectives are often conflicting, and finding the optimal point of the trade-off has been a challenge for breeders. Genomic selection (GS) provides a revolutionary tool for the genetic improvement of farm animals. At the same time, it can increase inbreeding and produce a more rapid depletion of genetic variability of the selected traits in future generations. Optimum-contribution selection (OCS) represents an approach to maximize genetic gain while constraining inbreeding within a targeted range. In the present study, 515 Ningxiang pigs were genotyped with the Illumina Porcine SNP60 array or the GeneSeek Genomic Profiler Porcine 50K array. The Ningxiang pigs were found to be highly inbred at the genomic level. Average locus-wise inbreeding coefficients were 0.41 and 0.37 for the two SNP arrays used, whereas genomic inbreeding coefficients based on runs of homozygosity were 0.24 and 0.25, respectively. Simulated phenotypic data were used to assess the utility of genomic OCS (GOCS) in comparison with GS without inbreeding control. GOCS was conducted under two scenarios, selecting sires only (GOCS_S) or selecting sires and dams (GOCS_SD), while kinships were constrained on selected parents. The genetic gain for average daily body weight gain (ADG) per generation was between 18.99 and 20.55 g with GOCS_S, and between 23.20 and 28.92 with GOCS_SD, and it varied from 25.38 to 48.38 g under GS without controlling inbreeding. While the rate of genetic gain per generation obtained using GS was substantially larger than that obtained by the two scenarios of genomic OCS in the beginning generations of selection, the difference in the genetic gain of ADG between GS and GOCS reduced quickly in latter generations. At generation ten, the difference in the realized rates of genetic gain between GS and GOCS_SD diminished and ended up with even a slightly higher genetic gain with GOCS_SD, due to the rapid loss of genetic variance with GS and fixation of causative genes. The rate of inbreeding was mostly maintained below 5% per generation with genomic OCS, whereas it increased to between 10.5% and 15.3% per generation with GS. Therefore, genomic OCS appears to be a sustainable strategy for the genetic improvement of local breeds such as Ningxiang pigs, but keeping mind that a variety of GOCS methods exist and the optimal forms remain to be exploited further.

Estimation of Genomic Breed Composition for Purebred and Crossbred Animals Using Sparsely Regularized Admixture Models.

A variety of statistical methods, such as admixture models, have been used to estimate genomic breed composition (GBC). These methods, however, tend to produce non-zero components to reference breeds that shared some genomic similarity with a test animal. These non-essential GBC components, in turn, offset the estimated GBC for the breed to which it belongs. As a result, not all purebred animals have 100% GBC of their respective breeds, which statistically indicates an elevated false-negative rate in the identification of purebred animals with 100% GBC as the cutoff. Otherwise, a lower cutoff of estimated GBC will have to be used, which is arbitrary, and the results are less interpretable. In the present study, three admixture models with regularization were proposed, which produced sparse solutions through suppressing the noise in the estimated GBC due to genomic similarities. The regularization or penalty forms included the L1 norm penalty, minimax concave penalty (MCP), and smooth clipped absolute deviation (SCAD). The performances of these regularized admixture models on the estimation of GBC were examined in purebred and composite animals, respectively, and compared to that of the non-regularized admixture model as the baseline model. The results showed that, given optimal values for λ, the three sparsely regularized admixture models had higher power and thus reduced the false-negative rate for the breed identification of purebred animals than the non-regularized admixture model. Of the three regularized admixture models, the two with a non-convex penalty outperformed the one with L1 norm penalty. In the Brangus, a composite cattle breed, estimated GBC were roughly comparable among the four admixture models, but all the four models underestimated the GBC for these composite animals when non-ancestral breeds were included as the reference. In conclusion, the admixture models with sparse regularization gave more parsimonious, consistent and interpretable results of estimated GBC for purebred animals than the non-regularized admixture model. Nevertheless, the utility of regularized admixture models for estimating GBC in crossbred or composite animals needs to be taken with caution.

[Impacts of SNP genotyping call rate and SNP genotyping error rate on imputation accuracy inHolsteincattle].

Single nucleotide polymorphism (SNP) chips have been widely used in genetic studies and breeding applications in animal and plant species. The quality of SNP genotypes is of paramount importance. More often than not, there are situations in which a number of genotypes may fail, requiring them to be imputed. There are also situations in which ungenotyped loci need to be imputed between different chips, or high-density genotypes need to be imputed based on low-density genotypes. Under these circumstances, the validity and reliability of subsequent data analyses is subject to the accuracy of these imputed genotypes. For justifying a better understanding of factors affecting imputation accuracy, in the present study, the impacts of SNP genotyping call rate and SNP genotyping error rate on the accuracy of genotype imputation were investigated under two scenarios in 20 116 U.S. Holstein cattle, each genotyped with a GGP 50K SNP chip. When the two factors were not correlated in scenario 1, simulated genotyping call rate varied from 50% to 100% and simulated genotyping error rate changed from 0% to 50%, with both factors being independent of each other. In scenario 2, genotyping error rates were correlated with genotyping call rate, and the relationship was set up by fitting a linear regression model between the two variables on a real dataset. That is, the simulated SNP call rate varied from 100% to 50% whereas the SNP genotyping rate changed from 0% to 13.55%. Finally, a 5-fold cross-validation was used to assess the subsequent imputation accuracy. The results showed that when original SNP genotyping call rate were independent of SNP genotyping error rate, the imputation accuracy did not change significantly with the original genotyping call rate (P>0.05), but it decreased significantly as the genotyping error rate increased (P<0.01). However, when original genotyping call rate was negatively correlated with genotyping error rate, the imputation error increased with elevated original genotyping error rate. In both scenarios, genotyping call rate needs to be no less than 0.90 in order to obtain 98% or higher genotype imputation accuracy. The present results can provide guidance for establishing quality assurance criteria for SNP genotyping in practice.

Comparing SNP panels and statistical methods for estimating genomic breed composition of individual animals in ten cattle breeds.

BACKGROUND: SNPs are informative to estimate genomic breed composition (GBC) of individual animals, but selected SNPs for this purpose were not made available in the commercial bovine SNP chips prior to the present study. The primary objective of the present study was to select five common SNP panels for estimating GBC of individual animals initially involving 10 cattle breeds (two dairy breeds and eight beef breeds). The performance of the five common SNP panels was evaluated based on admixture model and linear regression model, respectively. Finally, the downstream implication of GBC on genomic prediction accuracies was investigated and discussed in a Santa Gertrudis cattle population. RESULTS: There were 15,708 common SNPs across five currently-available commercial bovine SNP chips. From this set, four subsets (1,000, 3,000, 5,000, and 10,000 SNPs) were selected by maximizing average Euclidean distance (AED) of SNP allelic frequencies among the ten cattle breeds. For 198 animals presented as Akaushi, estimated GBC of the Akaushi breed (GBCA) based on the admixture model agreed very well among the five SNP panels, identifying 166 animals with GBCA = 1. Using the same SNP panels, the linear regression approach reported fewer animals with GBCA = 1. Nevertheless, estimated GBCA using both models were highly correlated (r = 0.953 to 0.992). In the genomic prediction of a Santa Gertrudis population (and crosses), the results showed that the predictability of molecular breeding values using SNP effects obtained from 1,225 animals with no less than 0.90 GBC of Santa Gertrudis (GBCSG) decreased on crossbred animals with lower GBCSG. CONCLUSIONS: Of the two statistical models used to compute GBC, the admixture model gave more consistent results among the five selected SNP panels than the linear regression model. The availability of these common SNP panels facilitates identification and estimation of breed compositions using currently-available bovine SNP chips. In view of utility, the 1 K panel is the most cost effective and it is convenient to be included as add-on content in future development of bovine SNP chips, whereas the 10 K and 16 K SNP panels can be more resourceful if used independently for imputation to intermediate or high-density genotypes.

[Estimating genomic breed composition of individual animals using selected SNPs].

Natural and artificial selection, geographical segregation and genetic drift can result in differentiation of allelic frequencies of single nucleotide polymorphism (SNP) at many loci in the animal genome. For individuals whose ancestors originated from different populations, their genetic compositions exhibit multiple components correlated with the genotypes or allele frequencies of these breeds or populations. Therefore, by using an appropriate statistical method, one can estimate the genomic contribution of each breed (ancestor) to the genome of each individual animal, which is referred to as the genomic breed composition (GBC). This paper reviews the principles, statistical methods and steps for estimating GBC of individual animals using SNP genotype data. Based on a linear regression model and an admixture model respectively, the protocols were demonstrated by the breed characterization of 198 purported Akaushi cattle, which included selection of reference SNPs and reference individual animals, and computing GBC for animals to be evaluated. The reference populations consist of 36 574 cattle from five cattle breeds (Akaushi, Angus, Hereford, Holstein and Jersey), each genotyped on either a 40K or 50K SNP chip. Four common SNP panels scanned from commercial chips for estimating GBC of individual animals are optimally selected, thereby expanding the functionalities of the currently available commercial SNP chips. It remains to be explored in future studies as to how estimated GBC can be incorporated to improve the accuracies on genomic prediction in purebred animals and crossbreds as well.

Comparing strategies for selection of low-density SNPs for imputation-mediated genomic prediction in U. S. Holsteins.

SNP chips are commonly used for genotyping animals in genomic selection but strategies for selecting low-density (LD) SNPs for imputation-mediated genomic selection have not been addressed adequately. The main purpose of the present study was to compare the performance of eight LD (6K) SNP panels, each selected by a different strategy exploiting a combination of three major factors: evenly-spaced SNPs, increased minor allele frequencies, and SNP-trait associations either for single traits independently or for all the three traits jointly. The imputation accuracies from 6K to 80K SNP genotypes were between 96.2 and 98.2%. Genomic prediction accuracies obtained using imputed 80K genotypes were between 0.817 and 0.821 for daughter pregnancy rate, between 0.838 and 0.844 for fat yield, and between 0.850 and 0.863 for milk yield. The two SNP panels optimized on the three major factors had the highest genomic prediction accuracy (0.821-0.863), and these accuracies were very close to those obtained using observed 80K genotypes (0.825-0.868). Further exploration of the underlying relationships showed that genomic prediction accuracies did not respond linearly to imputation accuracies, but were significantly affected by genotype (imputation) errors of SNPs in association with the traits to be predicted. SNPs optimal for map coverage and MAF were favorable for obtaining accurate imputation of genotypes whereas trait-associated SNPs improved genomic prediction accuracies. Thus, optimal LD SNP panels were the ones that combined both strengths. The present results have practical implications on the design of LD SNP chips for imputation-enabled genomic prediction.

Optimal Design of Low-Density SNP Arrays for Genomic Prediction: Algorithm and Applications.

Low-density (LD) single nucleotide polymorphism (SNP) arrays provide a cost-effective solution for genomic prediction and selection, but algorithms and computational tools are needed for the optimal design of LD SNP chips. A multiple-objective, local optimization (MOLO) algorithm was developed for design of optimal LD SNP chips that can be imputed accurately to medium-density (MD) or high-density (HD) SNP genotypes for genomic prediction. The objective function facilitates maximization of non-gap map length and system information for the SNP chip, and the latter is computed either as locus-averaged (LASE) or haplotype-averaged Shannon entropy (HASE) and adjusted for uniformity of the SNP distribution. HASE performed better than LASE with ≤1,000 SNPs, but required considerably more computing time. Nevertheless, the differences diminished when >5,000 SNPs were selected. Optimization was accomplished conditionally on the presence of SNPs that were obligated to each chromosome. The frame location of SNPs on a chip can be either uniform (evenly spaced) or non-uniform. For the latter design, a tunable empirical Beta distribution was used to guide location distribution of frame SNPs such that both ends of each chromosome were enriched with SNPs. The SNP distribution on each chromosome was finalized through the objective function that was locally and empirically maximized. This MOLO algorithm was capable of selecting a set of approximately evenly-spaced and highly-informative SNPs, which in turn led to increased imputation accuracy compared with selection solely of evenly-spaced SNPs. Imputation accuracy increased with LD chip size, and imputation error rate was extremely low for chips with ≥3,000 SNPs. Assuming that genotyping or imputation error occurs at random, imputation error rate can be viewed as the upper limit for genomic prediction error. Our results show that about 25% of imputation error rate was propagated to genomic prediction in an Angus population. The utility of this MOLO algorithm was also demonstrated in a real application, in which a 6K SNP panel was optimized conditional on 5,260 obligatory SNP selected based on SNP-trait association in U.S. Holstein animals. With this MOLO algorithm, both imputation error rate and genomic prediction error rate were minimal.

Predicting expected progeny difference for marbling score in Angus cattle using artificial neural networks and Bayesian regression models.

BACKGROUND: Artificial neural networks (ANN) mimic the function of the human brain and are capable of performing massively parallel computations for data processing and knowledge representation. ANN can capture nonlinear relationships between predictors and responses and can adaptively learn complex functional forms, in particular, for situations where conventional regression models are ineffective. In a previous study, ANN with Bayesian regularization outperformed a benchmark linear model when predicting milk yield in dairy cattle or grain yield of wheat. Although breeding values rely on the assumption of additive inheritance, the predictive capabilities of ANN are of interest from the perspective of their potential to increase the accuracy of prediction of molecular breeding values used for genomic selection. This motivated the present study, in which the aim was to investigate the accuracy of ANN when predicting the expected progeny difference (EPD) of marbling score in Angus cattle. Various ANN architectures were explored, which involved two training algorithms, two types of activation functions, and from 1 to 4 neurons in hidden layers. For comparison, BayesCπ models were used to select a subset of optimal markers (referred to as feature selection), under the assumption of additive inheritance, and then the marker effects were estimated using BayesCπ with π set equal to zero. This procedure is referred to as BayesCpC and was implemented on a high-throughput computing cluster. RESULTS: The ANN with Bayesian regularization method performed equally well for prediction of EPD as BayesCpC, based on prediction accuracy and sum of squared errors. With the 3K-SNP panel, for example, prediction accuracy was 0.776 using BayesCpC, and ranged from 0.776 to 0.807 using BRANN. With the selected 700-SNP panel, prediction accuracy was 0.863 for BayesCpC and ranged from 0.842 to 0.858 for BRANN. However, prediction accuracy for the ANN with scaled conjugate gradient back-propagation was lower, ranging from 0.653 to 0.689 with the 3K-SNP panel, and from 0.743 to 0.793 with the selected 700-SNP panel. CONCLUSIONS: ANN with Bayesian regularization performed as well as linear Bayesian regression models in predicting additive genetic values, supporting the idea that ANN are useful as universal approximators of functions of interest in breeding contexts.

A novel analytical method, Birth Date Selection Mapping, detects response of the Angus (Bos taurus) genome to selection on complex traits.

BACKGROUND: Several methods have recently been developed to identify regions of the genome that have been exposed to strong selection. However, recent theoretical and empirical work suggests that polygenic models are required to identify the genomic regions that are more moderately responding to ongoing selection on complex traits. We examine the effects of multi-trait selection on the genome of a population of US registered Angus beef cattle born over a 50-year period representing approximately 10 generations of selection. We present results from the application of a quantitative genetic model, called Birth Date Selection Mapping, to identify signatures of recent ongoing selection. RESULTS: We show that US Angus cattle have been systematically selected to alter their mean additive genetic merit for most of the 16 production traits routinely recorded by breeders. Using Birth Date Selection Mapping, we estimate the time-dependency of allele frequency for 44,817 SNP loci using genomic best linear unbiased prediction, generalized least squares, and BayesCπ analyses. Finally, we reconstruct the primary phenotypes that have historically been exposed to selection from a genome-wide analysis of the 16 production traits and gene ontology enrichment analysis. CONCLUSIONS: We demonstrate that Birth Date Selection Mapping utilizing mixed models corrects for time-dependent pedigree sampling effects that lead to spurious SNP associations and reveals genomic signatures of ongoing selection on complex traits. Because multiple traits have historically been selected in concert and most quantitative trait loci have small effects, selection has incrementally altered allele frequencies throughout the genome. Two quantitative trait loci of large effect were not the most strongly selected of the loci due to their antagonistic pleiotropic effects on strongly selected phenotypes. Birth Date Selection Mapping may readily be extended to temporally-stratified human or model organism populations.

An ensemble-based approach to imputation of moderate-density genotypes for genomic selection with application to Angus cattle.

Summary Imputation of moderate-density genotypes from low-density panels is of increasing interest in genomic selection, because it can dramatically reduce genotyping costs. Several imputation software packages have been developed, but they vary in imputation accuracy, and imputed genotypes may be inconsistent among methods. An AdaBoost-like approach is proposed to combine imputation results from several independent software packages, i.e. Beagle(v3.3), IMPUTE(v2.0), fastPHASE(v1.4), AlphaImpute, findhap(v2) and Fimpute(v2), with each package serving as a basic classifier in an ensemble-based system. The ensemble-based method computes weights sequentially for all classifiers, and combines results from component methods via weighted majority 'voting' to determine unknown genotypes. The data included 3078 registered Angus cattle, each genotyped with the Illumina BovineSNP50 BeadChip. SNP genotypes on three chromosomes (BTA1, BTA16 and BTA28) were used to compare imputation accuracy among methods, and the application involved the imputation of 50K genotypes covering 29 chromosomes based on a set of 5K genotypes. Beagle and Fimpute had the greatest accuracy among the six imputation packages, which ranged from 0·8677 to 0·9858. The proposed ensemble method was better than any of these packages, but the sequence of independent classifiers in the voting scheme affected imputation accuracy. The ensemble systems yielding the best imputation accuracies were those that had Beagle as first classifier, followed by one or two methods that utilized pedigree information. A salient feature of the proposed ensemble method is that it can solve imputation inconsistencies among different imputation methods, hence leading to a more reliable system for imputing genotypes relative to independent methods.

Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation.

BACKGROUND: Genomic selection is a recently developed technology that is beginning to revolutionize animal breeding. The objective of this study was to estimate marker effects to derive prediction equations for direct genomic values for 16 routinely recorded traits of American Angus beef cattle and quantify corresponding accuracies of prediction. METHODS: Deregressed estimated breeding values were used as observations in a weighted analysis to derive direct genomic values for 3570 sires genotyped using the Illumina BovineSNP50 BeadChip. These bulls were clustered into five groups using K-means clustering on pedigree estimates of additive genetic relationships between animals, with the aim of increasing within-group and decreasing between-group relationships. All five combinations of four groups were used for model training, with cross-validation performed in the group not used in training. Bivariate animal models were used for each trait to estimate the genetic correlation between deregressed estimated breeding values and direct genomic values. RESULTS: Accuracies of direct genomic values ranged from 0.22 to 0.69 for the studied traits, with an average of 0.44. Predictions were more accurate when animals within the validation group were more closely related to animals in the training set. When training and validation sets were formed by random allocation, the accuracies of direct genomic values ranged from 0.38 to 0.85, with an average of 0.65, reflecting the greater relationship between animals in training and validation. The accuracies of direct genomic values obtained from training on older animals and validating in younger animals were intermediate to the accuracies obtained from K-means clustering and random clustering for most traits. The genetic correlation between deregressed estimated breeding values and direct genomic values ranged from 0.15 to 0.80 for the traits studied. CONCLUSIONS: These results suggest that genomic estimates of genetic merit can be produced in beef cattle at a young age but the recurrent inclusion of genotyped sires in retraining analyses will be necessary to routinely produce for the industry the direct genomic values with the highest accuracy.

A primer on high-throughput computing for genomic selection.

High-throughput computing (HTC) uses computer clusters to solve advanced computational problems, with the goal of accomplishing high-throughput over relatively long periods of time. In genomic selection, for example, a set of markers covering the entire genome is used to train a model based on known data, and the resulting model is used to predict the genetic merit of selection candidates. Sophisticated models are very computationally demanding and, with several traits to be evaluated sequentially, computing time is long, and output is low. In this paper, we present scenarios and basic principles of how HTC can be used in genomic selection, implemented using various techniques from simple batch processing to pipelining in distributed computer clusters. Various scripting languages, such as shell scripting, Perl, and R, are also very useful to devise pipelines. By pipelining, we can reduce total computing time and consequently increase throughput. In comparison to the traditional data processing pipeline residing on the central processors, performing general-purpose computation on a graphics processing unit provide a new-generation approach to massive parallel computing in genomic selection. While the concept of HTC may still be new to many researchers in animal breeding, plant breeding, and genetics, HTC infrastructures have already been built in many institutions, such as the University of Wisconsin-Madison, which can be leveraged for genomic selection, in terms of central processing unit capacity, network connectivity, storage availability, and middleware connectivity. Exploring existing HTC infrastructures as well as general-purpose computing environments will further expand our capability to meet increasing computing demands posed by unprecedented genomic data that we have today. We anticipate that HTC will impact genomic selection via better statistical models, faster solutions, and more competitive products (e.g., from design of marker panels to realized genetic gain). Eventually, HTC may change our view of data analysis as well as decision-making in the post-genomic era of selection programs in animals and plants, or in the study of complex diseases in humans.

A comparison of prophylactic efficacy of tilmicosin and a new formulation of oxytetracycline in feedlot calves.

Two replicated-pen field studies were performed under commercial feedlot conditions in western Canada to compare the administration of long-acting oxytetracycline at 30 mg/kg body weight (BW) versus tilmicosin at 10 mg/kg BW to feedlot calves upon arrival at the feedlot. Ten thousand nine hundred and eighty-nine, recently weaned, auction market derived, crossbred beef steer and bull calves were randomly allocated upon arrival at the feedlot to one of 2 experimental groups as follows: oxytetracycline, which received intramuscular long-acting oxytetracycline (300 mg/mL formulation) at a rate of 30 mg/kg BW; or tilmicosin, which received subcutaneous tilmicosin (300 mg/mL formulation) at a rate of 10 mg/kg BW. There were 20 pens in each experimental group. In Study 1 and in the combined analysis, the initial undifferentiated fever (UF) treatment rate was significantly (P < 0.05) higher in the oxytetracycline group as compared with the tilmicosin group. There were no significant (P > or = 0.05) differences in first UF relapse, second UF relapse, third UF relapse, overall chronicity, overall rail, overall mortality, bovine respiratory disease (BRD) mortality, hemophilosis mortality, arthritis mortality, or miscellaneous mortality rates between the experimental groups in either study or in the combined analysis. In addition, there were no significant (P > or = 0.05) differences in initial weight, final weight, weight gain, days on feed, daily dry matter intake, average daily gain, or the dry matter intake to gain ratio between the experimental groups in either study or in the combined analyses. In the economic analysis, there was a net economic advantage of $5.22 CDN per animal in the oxytetracycline group, due to a lower prophylactic cost, even though the UF therapeutic cost was higher.

An evaluation of the relative efficacy of a new formulation of oxytetracycline for the treatment of undifferentiated fever in feedlot calves in western Canada.

A field trial was performed under commercial feedlot conditions in western Canada to compare the efficacy of a new formulation of long-acting oxytetracycline (LA 30) to a standard long-acting oxytetracycline formulation (LA 20) and florfenicol (FLOR) for the treatment of undifferentiated fever (UF) in calves that received metaphylactic tilmicosin upon arrival at the feed-lot. Seven hundred and ninety-seven recently weaned, auction market derived, crossbred, beef calves suffering from UF were allocated to 1 of 3 experimental groups as follows: LA 30, which received intramuscular long-acting oxytetracycline (300 mg/mL formulation) at the rate of 30 mg/kg body weight (BW) at the time of allocation; LA 20, which received intramuscular long-acting oxytetracycline (200 mg/mL formulation) at the rate of 20 mg/kg BW at the time of allocation; or FLOR, which received intramuscular florfenicol administered at the rate of 20 mg/kg BW at the time of allocation and again 48 hours later. Two hundred and sixty-six animals were allocated to the LA 30 group, 265 animals were allocated to the LA 20 group, and 266 animals were allocated to the FLOR group. The relative efficacy of the LA 30 group, as compared with the LA 20 and FLOR groups, was assessed by comparing relapse, chronicity, wastage, and mortality rates. The overall mortality (RR = 0.50) rate in the LA 30 group was significantly (P < 0.05) lower than in the LA 20 group. However, the overall chronicity (RR = 2.56) and overall wastage (RR = 6.97) rates of the LA 30 group were significantly (P < 0.05) higher than in the LA 20 group. There were no significant (P > or = 0.05) differences in UF relapse rates or cause specific mortality rates between the LA 30 and LA 20 groups. In the economic analysis, there was an advantage of $28.59 CDN per animal in the LA 30 group compared with the LA 20 group. The overall chronicity (RR = 2.25) and overall wastage (RR = 2.80) rates of the LA 30 group were significantly (P < 0.05) higher than the FLOR group. There were no significant (P > or = 0.05) differences in UF relapse rates, overall mortality rates, or cause specific mortality rates between the LA 30 and FLOR groups. In the economic analysis, there was an advantage of $12.90 CDN per animal in the LA 30 group compared with the FLOR group. In summary, the results of this study indicate that it is more cost-effective to use a new formulation of long-acting oxytetracycline (300 mg/mL formulation administered at a rate of 30 mg/kg BW) than a standard long-acting oxytetracycline formulation (200 mg/mL formulation administered at a rate of 20 mg/kg BW) or florfenicol for the treatment of UF in feedlot calves that have previously received metaphylactic tilmicosin upon arrival at the feedlot.