Direct Genomic Values Research Articles

The marker effects of a single-step random regression model for 4 test-day traits in German Holsteins

The single-step genomic model has become the golden standard for routine evaluation in livestock species, such as Holstein dairy cattle. The single-step genomic model with direct estimation of marker effects has been proven to be efficient in accurately accounting for millions of genotype records. For diverse applications including frequent genomic evaluation updates on a weekly basis, estimates of the marker effects from the single-step evaluations play a central role in genomic prediction. In this study we focused on exploring the marker effect estimates from the single-step evaluation. Phenotypic, genotypic, and pedigree data were taken from the official evaluation for German dairy breeds in April 2021. A multilactation random regression test-day model was applied to more than 242 million test-day records separately for 4 traits: milk, fat, and protein yields, and somatic cell scores (SCS). Approximately one million genotyped Holstein animals were considered in the single-step genomic evaluations including ∼21 million animals in pedigree. Deregressed multiple across-country breeding values of Holstein bulls having daughters outside Germany were integrated into the national test-day data to increase the reliability of genomic breeding values. To assess the stability and bias of the marker effects of the single-step model, test-day records of the last 4 yr were deleted, and the integrated bulls born in the last 4 yr were truncated from the complete phenotypic dataset. Estimates of the marker effects were shown to be highly correlated, with correlations ∼0.9, between the full and truncated evaluations. Regression slope values of the marker-effect estimates from the full on the truncated evaluations were all close to their expected value, being ∼1.03. Calculated using random regression coefficients of the marker effect estimates, drastically different shapes of the genetic lactation curve were seen for 2 markers on chromosome 14 for the 4 test-day traits. The contribution of individual chromosomes to the total additive genetic variances seemed to follow the polygenic inheritance mode for protein yield and SCS. However, chromosome 14 was found to make an exceptionally large contribution to the total additive genetic variance for milk and fat yields because of markers near the major gene DGAT1. For the first lactation test-day traits, we obtained ∼0 correlations of chromosomal direct genomic values between any pair of the chromosomes; no spurious correlations were found in our analysis, thanks to the large reference population. For trait milk yield, chromosomal direct genomic values appeared to have a large variation in the between-lactation correlations among the chromosomes, especially between first and second or third lactations. The optimal features of the random regression test-day model and the single-step marker model allowed us to track the differences in the shapes of genetic lactation curves down to the individual markers. Furthermore, the single-step random regression test-day model enabled us to better understand the inheritance mode of the yield traits and SCS (e.g., variable chromosomal contributions to the total additive genetic variance and to the genetic correlations between lactations).

Evaluating the potential of (epi)genotype-by-low pass nanopore sequencing in dairy cattle: a study on direct genomic value and methylation analysis

BackgroundGenotype-by-sequencing has been proposed as an alternative to SNP genotyping arrays in genomic selection to obtain a high density of markers along the genome. It requires a low sequencing depth to be cost effective, which may increase the error at the genotype assigment. Third generation nanopore sequencing technology offers low cost sequencing and the possibility to detect genome methylation, which provides added value to genotype-by-sequencing. The aim of this study was to evaluate the performance of genotype-by-low pass nanopore sequencing for estimating the direct genomic value in dairy cattle, and the possibility to obtain methylation marks simultaneously.ResultsLatest nanopore chemistry (LSK14 and Q20) achieved a modal base calling accuracy of 99.55%, whereas previous kit (LSK109) achieved slightly lower accuracy (99.1%). The direct genomic value accuracy from genotype-by-low pass sequencing ranged between 0.79 and 0.99, depending on the trait (milk, fat or protein yield), with a sequencing depth as low as 2 × and using the latest chemistry (LSK114). Lower sequencing depth led to biased estimates, yet with high rank correlations. The LSK109 and Q20 achieved lower accuracies (0.57–0.93). More than one million high reliable methylated sites were obtained, even at low sequencing depth, located mainly in distal intergenic (87%) and promoter (5%) regions.ConclusionsThis study showed that the latest nanopore technology in useful in a LowPass sequencing framework to estimate direct genomic values with high reliability. It may provide advantages in populations with no available SNP chip, or when a large density of markers with a wide range of allele frequencies is needed. In addition, low pass sequencing provided nucleotide methylation status of > 1 million nucleotides at ≥ 10 × , which is an added value for epigenetic studies.

Genomic prediction of fertility and calving traits in Holstein cattle based on models including epistatic genetic effects.

The goal of this study was to investigate whether the inclusion of genomic information and epistatic (additive by additive) genetic effects would increase the accuracy of predicting phenotypes adjusted for known environmental effects, reduce prediction bias and minimize the confounding between additive and additive by additive epistatic effects on fertility and calving traits in Holstein cattle. Phenotypic and genotypic records were available for 6090 cows. Eight cow traits were assessed including 56-day nonreturn rate (NRR), number of services (NS), calving to first insemination (CTFS), first insemination to conception (FSTC), gestation length (GL), calving ease (CE), stillbirth (SB) and calf size (CZ). Four scenarios were assessed for their ability to predict adjusted phenotypes, which included: (1) traditional pedigree-based Best Linear Unbiased Prediction (P-BLUP) for additive genetic effects (PA); (2) P-BLUP for additive and epistatic (additive by additive) genetic effects (PAE); (3) genomic BLUP (G-BLUP) for additive genetic effects (GA); and (4) G-BLUP for additive and epistatic genetic effects (GAEn, where n = 1-3 depending on the alternative ways to construct the epistatic genomic matrix used). Constructing epistatic relationship matrix as the Hadamard product of the additive genomic relationship matrix (GAE1), which is the usual method and implicitly assumes a model that fits all pairwise interactions between markers twice and includes the interactions of the markers with themselves (dominance). Two additional constructions of the epistatic genomic relationship matrix were compared to test whether removing the double counting of interactions and the interaction of the markers with themselves (GAE2), and removing double counting of interactions between markers, but including the interaction of the markers with themselves (GAE3) would had an impact on the prediction and estimation error correlation (i.e. confounding) between additive and epistatic genetic effects. Fitting epistatic genetic effects explained up to 5.7% of the variance for NRR (GAE3), 7.7% for NS (GAE1), 11.9% for CTFS (GAE3), 11.1% for FSTC (GAE2), 25.7% for GL (GAE1), 2.3% for CE (GAE1), 14.3% for SB (GAE3) and 15.2% for CZ (GAE1). Despite a substantial proportion of variance being explained by epistatic effects for some traits, the prediction accuracies were similar or lower for GAE models compared with pedigree models and genomic models without epistatic effects. Although the prediction accuracy of direct genomic values did not change significantly between the three variations of the epistatic genetic relationship matrix used, removing the interaction of the markers with themselves reduced the confounding between additive and additive by additive epistatic effects. These results suggest that epistatic genetic effects contribute to the variance of some fertility and calving traits in Holstein cattle. However, the inclusion of epistatic genetic effects in the genomic prediction of these traits is complex and warrant further investigation.

Association between body condition genomic values and feed intake, milk production, and body weight in French Holstein cows.

Body condition score (BCS) offers a good estimate of the amount of stored fat on the body, and its variations can be used as a proxy for energy balance. Many countries have implemented a genomic evaluation of BCS, including France, where estimated breeding values are based on an individual BCS determination during the first lactation. In this article, we investigate the degree to which this genomic estimated breeding value based on a single phenotype record per cow might reflect different profiles of body reserves throughout lactation and be used to predict, and perhaps limit, their mobilization during early lactation. We also investigate whether selection on BCS affects other traits. A data set including 686 lactations of 435 Holstein cows from 3 experimental farms not used in the reference population for genomic evaluation was used to estimate the effects of the BCS direct genomic value (iBCS) on BCS, body weight, feed intake, milk production, and fat and protein contents throughout the lactation period. For each trait, the model included different iBCS regressions and an effect of the direct genomic value of the trait itself when available. It thus appeared that cows with a positive iBCS always had a higher BCS than negative iBCS cows, whatever the lactation stage, and that this difference increased during the first 6 mo to reach a difference of 0.8 point. A similar effect was seen regarding body weight, but it was the opposite for milk production, with negative iBCS cows producing slightly more milk (difference of about 3% over lactation). Feed intake increased slightly faster at the beginning of lactation for cows with positive iBCS. Therefore, iBCS is a promising tool that could help to limit intense mobilization during early lactation. Should feed efficiency be included in the breeding goal, greater attention should be paid to BCS to avoid further body mobilization in early lactation.

Comparison of Marker Effects and Breeding Values at Two Levels at THI for Milk Yield and Quality Traits in Brazilian Holstein Cows.

Genomic tools can help in the selection of animals genetically resistant to heat stress, especially the genome-wide association studies (GWAS). The objective of this study was to compare the variance explained by SNPs and direct genomic breeding values (DGVs) at two levels of a temperature and humidity index (THI). Records of milk yield (MY), somatic cell score (SCS), and percentages of casein (CAS), saturated fatty acids (SFA), and unsaturated fatty acids (UFA) in milk from 1157 Holstein cows were used. Traditional breeding values (EBV) were determined in a previous study and used as pseudo-phenotypes. Two levels of THI (heat comfort zone and heat stress zone) were used as environments and were treated as "traits" in a bi-trait model. The GWAS was performed using the genomic best linear unbiased prediction (GBLUP) method. Considering the top 50 SNPs, a total of 36 SNPs were not common between environments, eight of which were located in gene regions related to the evaluated traits. Even for those SNPs that had differences in their explained variances between the two environments, the differences were very small. The animals showed virtually no rank order, with rank correlation values of 0.90, 0.88, 1.00, 0.88, and 0.97 for MY, CAS, SCS, SFA, and UFA, respectively. The small difference between the environments studied can be attributed to the small difference in the pseudo-phenotypes used between the environments, on-farm acclimation, the polygenic nature of the traits, and the THI values studied near the threshold between comfort and heat stress. It is recommended that future studies be conducted with a larger number of animals and at more extreme THI levels.

On the equivalence between marker effect models and breeding value models and direct genomic values with the Algorithm for Proven and Young

BackgroundSingle-step genomic predictions obtained from a breeding value model require calculating the inverse of the genomic relationship matrix ({mathbf{G}}^{-1}). The Algorithm for Proven and Young (APY) creates a sparse representation of {mathbf{G}}^{-1} with a low computational cost. APY consists of selecting a group of core animals and expressing the breeding values of the remaining animals as a linear combination of those from the core animals plus an error term. The objectives of this study were to: (1) extend APY to marker effects models; (2) derive equations for marker effect estimates when APY is used for breeding value models, and (3) show the implication of selecting a specific group of core animals in terms of a marker effects model.ResultsWe derived a family of marker effects models called APY-SNP-BLUP. It differs from the classic marker effects model in that the row space of the genotype matrix is reduced and an error term is fitted for non-core animals. We derived formulas for marker effect estimates that take this error term in account. The prediction error variance (PEV) of the marker effect estimates depends on the PEV for core animals but not directly on the PEV of the non-core animals. We extended the APY-SNP-BLUP to include a residual polygenic effect and accommodate non-genotyped animals. We show that selecting a specific group of core animals is equivalent to select a subspace of the row space of the genotype matrix. As the number of core animals increases, subspaces corresponding to different sets of core animals tend to overlap, showing that random selection of core animals is algebraically justified.ConclusionsThe APY-(ss)GBLUP models can be expressed in terms of marker effect models. When the number of core animals is equal to the rank of the genotype matrix, APY-SNP-BLUP is identical to the classic marker effects model. If the number of core animals is less than the rank of the genotype matrix, genotypes for non-core animals are imputed as a linear combination of the genotypes of the core animals. For estimating SNP effects, only relationships and estimated breeding values for core animals are needed.

Direct genomic value daughter pregnancy rate and services per conception are associated with characteristics of day-16 conceptuses and hormone signaling for maternal recognition of pregnancy in lactating Holstein cows

Our objective was to test the hypothesis that sorting Holstein dairy cows with high direct genomic value for daughter pregnancy rate (DGV-DPR) and low services per conception records improves interferon tau production and endometrial signaling for maternal recognition of pregnancy. Lactating Holstein cows (n = 86) were sorted by DGV-DPR into high- and low-fertility groups with the Clarifide ( https://www.zoetisus.com/animal-genetics/dairy/clarifide/clarifide.aspx# ) DNA test and low and high services per conception from farm records. High-fertility nonpregnant (NP) cows were not inseminated (n = 7) and served as controls. High-fertility pregnant (HP; n = 7) and low-fertility pregnant (LP; n = 6) cows received timed AI. Conceptuses, uterine flushings, endometrial biopsies, and peripheral blood mononuclear cells were collected on d 16 of pregnancy or the estrous cycle. Conceptuses from HP cows were longer ( P < 0.05) than those from LP conceptuses. Concentrations of interferon tau in uterine flushings were greater ( P < 0.05) in HP compared with LP and NP cows, positively correlated with DGV-DPR (r = +0.68; P < 0.05), and negatively correlated with services per conception (r = −0.59; P < 0.05). Endometrial ISG15 mRNA and protein were upregulated ( P < 0.05) in HP compared with NP. Low-fertility pregnant cows had intermediate endometrial levels of ISG15 mRNA between HP and NP, and ISG15 protein tended ( P < 0.10) to be greater in LP compared with NP cows. This tendency was also observed for increased ISG15 in HP compared with NP in peripheral blood mononuclear cells. Results supported the hypothesis and provided evidence to suggest that Holstein cows with greater DNA-based DGV-DPR values and fewer services per conception have increased conceptus length and enhanced signaling for maternal recognition of pregnancy.

Infrared Predictions Are a Valuable Alternative to Actual Measures of Dry-Cured Ham Weight Loss in the Training of Genome-Enabled Prediction Models.

Simple SummaryExcessive ham weight losses during dry-curing (WL) result in a loss of marketable product, hindering the quality of dry-cured hams. Genetic selection for reducing WL requires individual traceability of hams throughout the dry-curing process and the measurement of WL is expensive, time-consuming, and can be performed only at the end of seasoning, resulting in long generation intervals. Infrared spectroscopy provides early, cost-effective, high-throughput predictions of WL that are highly genetically correlated with the actual measures. This study focused on the accuracy of genomic prediction models for observed and infrared-predicted WL. Models were tested on crossbred pigs and their purebred sires in random cross-validation and in a leave-one-family-out training-validation scheme. Accuracy of prediction of sire genetic merit, estimated from crossbred training data for actual ham WL, was 0.38, slightly higher than the accuracy attainable by a model trained on infrared-predicted WL (0.32). While the accuracy of genomic predictions is satisfactory for both the observed and infrared-predicted WL, the use of infrared predictions results in considerably lower phenotyping costs, enabling the construction of larger reference populations.Selection to reduce ham weight losses during dry-curing (WL) requires individual traceability of hams throughout dry-curing, with high phenotyping costs and long generation intervals. Infrared spectroscopy enables cost-effective, high-throughput phenotyping for WL 24 h after slaughter. Direct genomic values (DGV) of crossbred pigs and their purebred sires were estimated, for observed (OB) and infrared-predicted WL (IR), through models developed from 640 and 956 crossbred pigs, respectively. Five Bayesian models and two pseudo-phenotypes (estimated breeding value, EBV, and adjusted phenotype) were tested in random cross-validation and leave-one-family-out validation. The use of EBV as pseudo-phenotypes resulted in the highest accuracies. Accuracies in leave-one-family-out validation were much lower than those obtained in random cross-validation but still satisfactory and very similar for both traits. For sires in the leave-one-family-out validation scenario, the correlation between the DGV for IR and EBV for OB was slightly lower (0.32) than the correlation between the DGV for OB and EBV for OB (0.38). While genomic prediction of OB and IR can be equally suggested to be incorporated in future selection programs aiming at reducing WL, the use of IR enables an early, cost-effective phenotyping, favoring the construction of larger reference populations, with accuracies comparable to those achievable using OB phenotype.

Genomic prediction in Nordic Red dairy cattle considering breed origin of alleles

This study investigated the reliability of genomic prediction (GP) using breed origin of alleles (BOA) approach in the Nordic Red (RDC) population, which has an admixed population structure. The RDC population consists of animals with varying degrees of genetic materials from the Danish Red (RDM), Swedish Red (SRB), Finnish Ayrshire (FAY), and Holstein (HOL) because bulls have been used across the breeds. The BOA approach was tested using 39,550 RDC animals in the reference population and 11,786 in the validation population. Deregressed proofs (DRP) of milk, fat and protein were used as response variable for GP. Direct genomic breeding values (DGV) for animals in the validation population were calculated with (BOA model) or without (joint model) considering breed origin of alleles. The joint model assumed homogeneous marker effects and a single set of marker effects were estimated, whereas BOA model assumed heterogeneous marker effects, and different sets of marker effects were estimated across the breeds. For the BOA approach, we tested scenarios assuming both correlated (BOA_cor) and uncorrelated (BOA_uncor) marker effects between the breeds. Additionally, we investigated GP using a standard Illumina 50K chip and including SNP selected from imputed whole-genome sequencing (50K+WGS). We also studied the effect of estimating (co)variances for genome regions of different sizes to exploit the information of the genome regions contributing to the (co)variance between the breeds. Region sizes were set as 1 SNP, a group of 30 or 100 adjacent SNP, or the whole genome. Reliability of DGV was measured as squared correlations between DGV and DRP divided by the reliability of DRP. Across the 3 traits, in general, RS30 and RS100 SNP yielded the highest reliabilities. Including WGS SNP improved reliabilities in almost all scenarios (0.297 on average for 50K and 0.307 on average for 50K+WGS). The BOA_uncor (0.233 on average) was inferior to the joint model (0.339 on average), but the reliabilities obtained using BOA_cor (0.334 on average) in most cases were not significantly different from those obtained using the joint model. The results indicate that both including additional whole-genome sequencing SNP and dividing the genome into fixed regions improve GP in the RDC. The BOA models have the potential to increase the reliability of GP, but the benefit is limited in populations with a high exchange of genetic material for a long time, as is the case for RDC.

Review on Animal Diseases Resistance and Adaptation Improvement through Molecular Genetics

This review is conducted to highlight the application of molecular genetics for improvement of diseases resistance in livestock. Diseases prevention is a serious issue to restrict economic damage due to economically important diseases outbreaks worldwide. Over the last decades, the excessive use of drugs has been criticised because of the possible development of drug-resistant zoonotic organisms and the potential dangers of drug residues in food animal products for human consumption. To prevent this problem following the ways for host resistance improvement to disease is a low cost and sustainable approach. So, this literature review was done on the theory of how to use molecular markers to select for quantitative trait loci (QTL) in genetic improvement programs, both within populations and for introgression of QTL from one population to another. If the effect of each marker is known, then an animal with no available phenotype for disease can be genotyped and its direct genomic value can be estimated based on its genotype only. The advantage of such a method is that accurate estimates of genetic merit can be achieved exploiting knowledge of the genotype of the animal even if the animal is very young. So, it is possible to put recommendation for further research on the generation of phenotypes that are resistant to diseases and development accurate bio-markers that can be readily measured in large numbers of animals at a relatively low cost.

Accuracy of genomic breeding values and predictive ability for postweaning liveweight and age at first calving in a Nellore cattle population with missing sire information.

The multiple sire system (MSS) is a common mating scheme in extensive beef production systems. However, MSS does not allow paternity identification and lead to inaccurate genetic predictions. The objective of this study was to investigate the implementation of single-step genomic BLUP (ssGBLUP) in different scenarios of uncertain paternity in the evaluation for 450-day adjusted liveweight (W450) and age at first calving (AFC) in a Nellore cattle population. To estimate the variance components using BLUP and ssGBLUP, the relationship matrix (A) with different proportions of animals with missing sires (MS) (scenarios 0, 25, 50, 75, and 100% of MS) was created. The genotyped animals with MS were randomly chosen, and ten replicates were performed for each scenario and trait. Five groups of animals were evaluated in each scenario: PHE, all animals with phenotypic records in the population; SIR, proven sires; GEN, genotyped animals; YNG, young animals without phenotypes and progeny; and YNGEN, young genotyped animals. The additive genetic variance decreased for both traits as the proportion of MS increased in the population when using the regular REML. When using the ssGBLUP, accuracies ranged from 0.13 to 0.47 for W450 and from 0.10 to 0.25 for AFC. For both traits, the prediction ability of the direct genomic value (DGV) decreased as the percentage of MS increased. These results emphasize that indirect prediction via DGV of young animals is more accurate when the SNP effects are derived from ssGBLUP with a reference population with known sires. The ssGBLUP could be applied in situations of uncertain paternity, especially when selecting young animals. This methodology is shown to be accurate, mainly in scenarios with a high percentage of MS.

Genome-wide association and genomic prediction for scrotal circumference in Hereford and Braford bulls

Scrotal circumference (SC) is widely used as a selection criterion for bulls in breeding programs, since it is easily assessed and correlated with several desirable reproductive traits. The objectives of this study were: to perform a genome-wide association study (GWAS) to identify genomic regions associated with SC adjusted for age (SCa) and for both age and weight (SCaw); to select Tag SNPs from GWAS to construct low-density panel for genomic prediction; and to compare the prediction accuracy of the SC through different methods for Braford and Hereford bulls from the same genetic breeding program. Data of SC from 18,172 bulls (30.4 ± 3.7 cm) and of genotypes from 131 sires and 3,545 animals were used. From GWAS, the top 1% of 1-Mb windows were observed on chromosome (BTA) 2, 20, 7, 8, 15, 3, 16, 27, 6 and 8 for SCa and on BTA 8, 15, 16, 21, 19, 2, 6, 5 and 10 for SCaw, representing 17.4% and 18.8% of the additive genetic variance of SCa and SCaw, respectively. The MeSH analysis was able to translate genomic information providing biological meanings of more specific gene functions related to the SCa and SCaw. The genomic enhancement methods, especially single step GBLUP, that combined phenotype and pedigree data with direct genomic values generated gains in accuracy in relation to pedigree BLUP, suggesting that genomic predictions should be applied to improve genetic gain and to narrow the generation interval compared to traditional methods. The proposed Tag-SNP panels may be useful for lower-cost commercial genomic prediction applications in the future, when the number of bulls in the reference population increases for SC in Hereford and Braford breeds.

Prediction of genomic breeding values of milk traits in Brazilian Saanen goats.

The study's objective was to compare the genomic prediction ability methods for the traits milk yield, milk composition and somatic cell count of Saanen Brazilian goats. Nine hundred forty goats, genotyped with an Axiom_OviCap (Caprine) panel, Affimetrix customized array with 62,557 single nucleotide polymorphisms (SNPs), were used for the genomic selection analyses. The genomic methods studied to estimate the effects of SNPs and direct genomic values (DGV) were as follows: (a) genomic BLUP (GBLUP), (b) Bayes Cπ and (c) Bayesian Lasso (BLASSO). Estimated breeding values (EBV) and deregressed estimated breeding values (dEBV) were used as response variables for the genomic predictions. The prediction ability was assessed by Pearson's correlation between DGV and response variables (EBV and dEBV). Regression coefficients of the response variables on the DGV were obtained to verify if the genomic predictions were biased. In addition, the mean square error of prediction (MSE) was used as a measure of verification of model fit to the data. The means of prediction accuracy, when EBV was used as a response variable, were 0.68, 0.68 and 0.67 for GBLUP, Bayes Cπ and BLASSO, respectively. With dEBV, the mean prediction accuracy was 0.50 for all models. The averages of the EBV regression coefficients on DGV were 1.08 for all models (GBLUP, Bayes Cπ and BLASSO), higher than those obtained for the regression coefficient of dEBV on DGV, which presented values of 1.05, 1.05 and 1.08 for GBLUP, Bayes Cπ and BLASSO, respectively. None of the methods stood out in terms of prediction ability; however, the GBLUP method was the most appropriate for estimating the DGV, in a slightly more reliable and less biased way, besides presenting the lowest computational cost. In the context of the present study, EBV was the preferred response variables considering the genomic prediction accuracy despite dEBV also presented lower bias.

Genetic parameters for methane emission traits in Australian dairy cows

Methane is a greenhouse gas of high interest to the dairy industry, with 57% of Australia's dairy emissions attributed to enteric methane. Enteric methane emissions also constitute a loss of approximately 6.5% of ingested energy. Genetic selection offers a unique mitigation strategy to decrease the methane emissions of dairy cattle, while simultaneously improving their energy efficiency. Breeding objectives should focus on improving the overall sustainability of dairy cattle by reducing methane emissions without negatively affecting important economic traits. Common definitions for methane production, methane yield, and methane intensity are widely accepted, but there is not yet consensus for the most appropriate method to calculate residual methane production, as the different methods have not been compared. In this study, we examined 9 definitions of residual methane production. Records of individual cow methane, dry matter intake (DMI), and energy corrected milk (ECM) were obtained from 379 animals and measured over a 5-d period from 12 batches across 5 yr using the SF6 tracer method and an electronic feed recording system, respectively. The 9 methods of calculating residual methane involved genetic and phenotypic regression of methane production on a combination of DMI and ECM corrected for days in milk, parity, and experimental batch using phenotypes or direct genomic values. As direct genomic values (DGV) for DMI are not routinely evaluated in Australia at this time, DGV for FeedSaved, which is derived from DGV for residual feed intake and estimated breeding value for bodyweight, were used. Heritability estimates were calculated using univariate models, and correlations were estimated using bivariate models corrected for the fixed effects of year-batch, days in milk, and lactation number, and fitted using a genomic relationship matrix. Residual methane production candidate traits had low to moderate heritability (0.10 ± 0.09 to 0.21 ± 0.10), with residual methane production corrected for ECM being the highest. All definitions of residual methane were highly correlated phenotypically (>0.87) and genetically (>0.79) with one another and moderately to highly with other methane candidate traits (>0.59), with high standard errors. The results suggest that direct selection for a residual methane production trait would result in indirect, favorable improvement in all other methane traits. The high standard errors highlight the importance of expanding data sets by measuring more animals for their methane emissions and DMI, or through exploration of proxy traits and combining data via international collaboration.

Association of genomically enhanced and parent average breeding values with cow performance in Nordic dairy cattle

This study compared the abilities of virgin heifer genomically enhanced breeding values (GEBV) and parent average breeding values (PA) to predict future cow performance. To increase confidence in genomic technology among farmers, a clear demonstration of the relationship between genomic predictions and future phenotypes is needed. We analyzed 12 different traits in first parity, including production, conformation, fertility, and other functional traits. Phenotype data were obtained from national milk recording schemes and breeding values from the Nordic Cattle Genetic Evaluation. Direct genomic breeding values were calculated using genomic BLUP and combined with traditional breeding values, using bivariate blending. The data covered 14,862 Red Dairy Cattle, 17,145 Holstein, and 7,330 Jersey genotyped virgin heifers born between 2013 and 2015 in Denmark, Finland, and Sweden. Phenotypes adjusted for systematic environmental effects were used as measures of cow performance. Two methods were used to compared virgin heifer GEBV and PA regarding their ability to predict future cow performance: (1) correlations between breeding values and adjusted phenotypes, (2) ranking cows into 4 quartiles for their virgin heifer GEBV or PA, and calculating actual cow performance for each quartile. We showed that virgin heifer GEBV predicted cow performance significantly better than PA for the vast majority of analyzed traits. The correlations with adjusted phenotypes were 38 to 136% higher for GEBV than for PA in Red Dairy Cattle, 42 to 194% higher for GEBV in Holstein, and 11 to 78% higher for GEBV in Jersey. The relative change between GEBV bottom and top quartiles compared with that between PA bottom and top quartiles ranged from 9 to 261% for RDC, 42 to 138% for Holstein, and 4 to 90% for Jersey. Hence, farmers in Denmark, Finland, and Sweden can have confidence in using genomic technology on their herds.

Current status of genomic evaluation.

Early application of genomic selection relied on SNP estimation with phenotypes or de-regressed proofs (DRP). Chips of 50k SNP seemed sufficient for an accurate estimation of SNP effects. Genomic estimated breeding values (GEBV) were composed of an index with parent average, direct genomic value, and deduction of a parental index to eliminate double counting. Use of SNP selection or weighting increased accuracy with small data sets but had minimal to no impact with large data sets. Efforts to include potentially causative SNP derived from sequence data or high-density chips showed limited or no gain in accuracy. After the implementation of genomic selection, EBV by BLUP became biased because of genomic preselection and DRP computed based on EBV required adjustments, and the creation of DRP for females is hard and subject to double counting. Genomic selection was greatly simplified by single-step genomic BLUP (ssGBLUP). This method based on combining genomic and pedigree relationships automatically creates an index with all sources of information, can use any combination of male and female genotypes, and accounts for preselection. To avoid biases, especially under strong selection, ssGBLUP requires that pedigree and genomic relationships are compatible. Because the inversion of the genomic relationship matrix (G) becomes costly with more than 100k genotyped animals, large data computations in ssGBLUP were solved by exploiting limited dimensionality of genomic data due to limited effective population size. With such dimensionality ranging from 4k in chickens to about 15k in cattle, the inverse of G can be created directly (e.g., by the algorithm for proven and young) at a linear cost. Due to its simplicity and accuracy, ssGBLUP is routinely used for genomic selection by the major chicken, pig, and beef industries. Single step can be used to derive SNP effects for indirect prediction and for genome-wide association studies, including computations of the P-values. Alternative single-step formulations exist that use SNP effects for genotyped or for all animals. Although genomics is the new standard in breeding and genetics, there are still some problems that need to be solved. This involves new validation procedures that are unaffected by selection, parameter estimation that accounts for all the genomic data used in selection, and strategies to address reduction in genetic variances after genomic selection was implemented.

Accuracies of direct genomic breeding values for birth and weaning weights of registered Charolais cattle in Mexico

Context Genomic prediction is now routinely used in many livestock species to rank individuals based on genomic breeding values (GEBV). Aims This study reports the first assessment aimed to evaluate the accuracy of direct GEBV for birth (BW) and weaning (WW) weights of registered Charolais cattle in Mexico. Methods The population assessed included 823 animals genotyped with an array of 77000 single nucleotide polymorphisms. Genomic prediction used genomic best linear unbiased prediction (GBLUP), Bayes C (BC), and single-step Bayesian regression (SSBR) methods in comparison with a pedigree-based BLUP method. Key results Our results show that the genomic prediction methods provided low and similar accuracies to BLUP. The prediction accuracy of GBLUP and BC were identical at 0.31 for BW and 0.29 for WW, similar to BLUP. Prediction accuracies of SSBR for BW and WW were up to 4% higher than those by BLUP. Conclusions Genomic prediction is feasible under current conditions, and provides a slight improvement using SSBR. Implications Some limitations on reference population size and structure were identified and need to be addressed to obtain more accurate predictions in liveweight traits under the prevalent cattle breeding conditions of Mexico.

Comparison of genomic prediction methods for evaluation of adaptation and productive efficiency traits in Braford and Hereford cattle

Various strategies to the use of genetic markers in genetic evaluations have been proposed and substantial genetic progress has been achieved in livestock breeding programs by incorporating genomic information. However, the implementation of genomic selection in beef cattle breeds, especially in composites, adapted to tropical and sub-tropical conditions, is well behind other cattle breeds. Alternative methods to perform genomic predictions for novel adaptation and productive efficiency traits in Brazilian Hereford and Braford need to be further investigated in order to improve the competitiveness of these beef cattle populations. Pedigree, phenotypes and a total of 3,680 genotypes (2,997 Braford and 683 Hereford) were used in this study to assess the predictive ability of genomic predictions for tick resistance, prepuce (navel) score, hair length, ocular pigmentation score, birth weight, pre- and post-weaning weight gain, scrotal circumference, body conformation, precocity, muscularity, and body size. Accuracies of Direct Genomic Values (DGVs) across all scenarios ranged from 0.08 to 0.66, and accuracies of Genomic Estimated Breeding Values (GEBVs) ranged from 0.08 to 0.81. In general, there were no statistical differences (P>0.05) between the DGVs and GEBVs accuracies, and the regression coefficients were not statistically different from 1 (i.e. not biased) for the large majority of traits and scenarios. Our findings indicate that accurate breeding values can be obtained at an early age for various adaptation and productive efficiency traits in Hereford and Braford. In addition, both DGVs and GEBVs predicted by single-step and two-step GBLUP methods produced similar accuracy and amount of bias. Considering the implementation's feasibility, the use the single-step GBLUP method for genetic evaluation of adaptation and productive efficiency traits in Braford and Hereford cattle is recommended.

Alternative SNP weighting for single-step genomic best linear unbiased predictor evaluation of stature in US Holsteins in the presence of selected sequence variants

Causal variants inferred from sequence data analysis are expected to increase accuracy of genomic selection. In this work we evaluated the gain in reliability of genomic predictions, for stature in US Holsteins, when adding selected sequence variants to a pre-existent SNP chip. Two prediction methods were tested: de-regressed proofs assuming heterogeneous (genomic BLUP; GBLUP) residual variances and by single-step GBLUP (ssGBLUP) using actual phenotypes. Phenotypic data included 3,999,631 records for stature on 3,027,304 Holstein cows. Genotypes on 54,087 SNP markers (54k) were available for 26,877 bulls. Additionally, 16,648 selected sequence variants were combined with the 54k markers, for a total of 70,735 (70k) markers. In all methods, SNP in the genomic relationship matrix (G) were unweighted or weighted iteratively, with weights derived either by SNP effects squared or by a nonlinear method that resembles BayesA (nonlinear A). Reliability of genomic predictions were obtained by cross validation. With unweighted G derived from 54k markers, the reliabilities (× 100) were 72.4 for GBLUP and 75.3 for ssGBLUP. With unweighted G derived from 70k markers, the reliabilities were 73.4 and 76.0, respectively. Weighting by nonlinear A changed reliabilities to 73.3, and 75.9, respectively. Addition of selected sequence variants had a small effect on reliabilities. Weighting by quadratic functions reduced reliabilities. Weighting by nonlinear A increased reliabilities for GBLUP but had only a small effect in ssGBLUP. Reliabilities for direct genomic values extracted from ssGBLUP using unweighted G with 54k were higher than reliabilities by any GBLUP. Thus, ssGBLUP seems to capture more information than GBLUP and there is less room for extra reliability. Improvements in GBLUP may be because the weights in G change the covariance structure, which can explain a proportion of the variance that is accounted for when a heterogeneous residual variance is assumed by considering a different number of daughters per bull.

Genomic selection improves the possibility of applying multiple breeding programs in different environments

One joint breeding program (BP) for different dairy cattle environments can be advantageous for genetic gain depending on the genetic correlation between environments (rg). The break-even correlation (rb) refers to the specific rg where genetic gain with 1 joint BP is equal to the genetic gain of 2 environment-specific BP. One joint BP has the highest genetic gain if rg is higher than rb, whereas 2 environment-specific BP have higher genetic gain if rg is lower than rb. Genetic gain in this context is evaluated from a breeding company's perspective that aims to improve genetic gain in both environments. With the implementation of genomic selection, 2 types of collaboration can be identified: exchanging breeding animals and exchanging genomic information. The aim of this study was to study genetic gain in multiple environments with different breeding strategies with genomic selection. The specific aims were (1) to find rb when applying genomic selection; (2) to assess how much genetic gain is lost when applying a suboptimal breeding strategy; (3) to study the effect of the reliability of direct genomic values, number of genotyped animals, and environments of different size on rb and genetic gain; and (4) to find rb from each environment's point of view. Three breeding strategies were simulated: 1 joint BP for both environments, 2 environment-specific BP with selection of bulls across environments, and 2 environment-specific BP with selection of bulls within environments. The rb was 0.65 and not different from rb with progeny-testing breeding programs when compared at the same selection intensity. The maximum loss in genetic gain in a suboptimal breeding strategy was 24%. A higher direct genomic value reliability and an increased number of genotyped selection candidates increased genetic gain, and the effect on rb was not large. A different size in 2 environments decreased rb by, at most, 0.10 points. From a large environment's point of view, 1 joint BP was the optimal breeding strategy in most scenarios. From a small environment's point of view, 1 joint BP was only the optimal breeding strategy at high rg. When the exchange of breeding animals between environments was restricted, genetic gain could still increase in each environment. This was due to the exchange of genomic information between environments, even when rg between environments were as low as 0.4. Thus, genomic selection improves the possibility of applying environment-specific BP.