Dataset Information

A Note on the Conditioning of the H^-1 Matrix Used in Single-Step GBLUP.

ABSTRACT: The single-step genomic BLUP (ssGBLUP) is used worldwide for the simultaneous genetic evaluation of genotyped and non-genotyped animals. It is easily extendible to all BLUP models by replacing the pedigree-based additive genetic relationship matrix (A) with an augmented pedigree-genomic relationship matrix (H). Theoretically, H does not introduce any artificially inflated variance. However, inflated genetic variances have been observed due to the incomparability between the genomic relationship matrix (G) and A used in H. Usually, G is blended and tuned with A22 (the block of A for genotyped animals) to improve its numerical condition and compatibility. If deflation/inflation is still needed, a common approach is weighting G-1-A22-1 in the form of τG-1-ωA22-1, added to A-1 to form H-1. In some situations, this can violate the conditional properties upon which H is built. Different ways of weighting the H-1 components (A-1, G-1, A22-1, and H-1 itself) were studied to avoid/minimise the violations of the conditional properties of H. Data were simulated on ten populations and twenty generations. Responses to weighting different components of H-1 were measured in terms of the regression of phenotypes on the estimated breeding values (the lower the slope, the higher the inflation) and the correlation between phenotypes and the estimated breeding values (predictive ability). Increasing the weight on H-1 increased the inflation. The responses to weighting G-1 were similar to those for H-1. Increasing the weight on A-1 (together with A22-1) was not influential and slightly increased the inflation. Predictive ability is a direct function of the slope of the regression line and followed similar trends. Responses to weighting G-1-A22-1 depend on the inflation/deflation of evaluations from A-1 to H-1 and the compatibility of the two matrices with the heritability used in the model. One possibility is a combination of weighting G-1-A22-1 and weighting H-1. Given recent advances in ssGBLUP, conditioning H-1 might become an interim solution from the past and then not be needed in the future.

SUBMITTER: Nilforooshan MA

PROVIDER: S-EPMC9686757 | biostudies-literature | 2022 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A Note on the Conditioning of the H-1 Matrix Used in Single-Step GBLUP.

Nilforooshan Mohammad Ali MA

Animals : an open access journal from MDPI 20221119 22

The single-step genomic BLUP (ssGBLUP) is used worldwide for the simultaneous genetic evaluation of genotyped and non-genotyped animals. It is easily extendible to all BLUP models by replacing the pedigree-based additive genetic relationship matrix (A) with an augmented pedigree-genomic relationship matrix (H). Theoretically, H does not introduce any artificially inflated variance. However, inflated genetic variances have been observed due to the incomparability between the ...[more]

PMID: 36428435

Similar Datasets

Project description:BackgroundThe preconditioned conjugate gradient (PCG) method is the current method of choice for iterative solving of genetic evaluations. The relative difference between two successive iterates and the relative residual of the system of equations are usually chosen as a termination criterion for the PCG method in animal breeding. However, our initial analyses showed that these two commonly used termination criteria may report that a PCG method applied to a single-step single nucleotide polymorphism best linear unbiased prediction (ssSNPBLUP) is not converged yet, whereas the solutions are accurate enough for practical use. Therefore, the aim of this study was to propose two termination criteria that have been (partly) developed in other fields, but are new in animal breeding, and to compare their behavior to that of the two termination criteria widely used in animal breeding for the PCG method applied to ssSNPBLUP. The convergence patterns of ssSNPBLUP were also compared to the convergence patterns of single-step genomic BLUP (ssGBLUP).ResultsBuilding upon previous work, we propose two termination criteria that take the properties of the system of equations into account. These two termination criteria are directly related to the relative error of the iterates with respect to the true solutions. Based on pig and dairy cattle datasets, we show that the preconditioned coefficient matrices of ssSNPBLUP and ssGBLUP have similar properties when using a second-level preconditioner for ssSNPBLUP. Therefore, the PCG method applied to ssSNPBLUP and ssGBLUP converged similarly based on the relative error of the iterates with respect to the true solutions. This similar convergence behavior between ssSNPBLUP and ssGBLUP was observed for both proposed termination criteria. This was, however, not the case for the termination criterion defined as the relative residual when applied to the dairy cattle evaluations.ConclusionOur results showed that the PCG method can converge similarly when applied to ssSNPBLUP and to ssGBLUP. The two proposed termination criteria always depicted these similar convergence behaviors, and we recommend them for comparing convergence properties of different models and for routine evaluations.

Project description:BackgroundSingle-step genomic best linear unbiased prediction (GBLUP) involves a joint analysis of individuals with genotype information, and their ancestors, descendants, or contemporaries, without recorded genotypes. Livestock applications typically represent populations with fewer individuals with genotypes relative to the number not genotyped. Most breeding programmes are structured, consisting of a nucleus tier in which selection drives genetic gains that are propagated through descendants that represent parents in multiplier and commercial tiers. In some cases, the genotypes in the nucleus tier are proprietary to a breeding company, and not publicly available for a whole industry analysis. Bayesian inference involves combining a defined description of prior information with new information to generate a posterior distribution that contains all available information on parameters of interest. A natural extension of Bayesian analysis would be to use information from the posterior distribution to define the prior distribution in a subsequent analysis.MethodsWe derive the mixed model equations for inference on breeding values for non genotyped individuals in that subset of the population that is of current interest, using only data on the performance of current individuals and their immediate pedigree, along with prior information defined from the posterior distribution of an external BLUP or single-step GBLUP analysis of the ancestors of the current population.DiscussionIdentical estimates of breeding values and their prediction error covariances for current animals of interest in the multiplier or commercial tier can be obtained without requiring neither the genomic relationship matrix nor genotypes of any of their ancestors in the nucleus tier, as can be obtained from a single analysis using pedigree, performance, and genomic information from all tiers. The Bayesian analysis of the current population does not require explicit information on unselected genotyped animals in the external population.

Project description:Brazilian livestock breeding programmes strive to enhance the genetics of beef cattle, with a strong emphasis on the Nellore breed, which has an extensive database and has achieved significant genetic progress in the last years. There are other indicine breeds that are economically important in Brazil; however, these breeds have more modest sets of phenotypes, pedigree and genotypes, slowing down their genetic progress as their predictions are less accurate. Combining several breeds in a multi-breed evaluation could help enhance predictions for those breeds with less information available. This study aimed to evaluate the feasibility of multi-breed, single-step genomic best linear unbiased predictor genomic evaluations for Nellore, Brahman, Guzerat and Tabapua. Multi-breed evaluations were contrasted to the single-breed ones. Data were sourced from the National Association of Breeders and Researchers of Brazil and included pedigree (4,207,516), phenotypic (328,748), and genomic (63,492) information across all breeds. Phenotypes were available for adjusted weight at 210 and 450 days of age, and scrotal circumference at 365 days of age. Various scenarios were evaluated to ensure pedigree and genomic information compatibility when combining different breeds, including metafounders (MF) or building the genomic relationship matrix with breed-specific allele frequencies. Scenarios were compared using the linear regression method for bias, dispersion and accuracy. The results showed that using multi-breed evaluations significantly improved accuracy, especially for smaller breeds like Guzerat and Tabapua. The validation statistics indicated that the MF approach provided accurate predictions, albeit with some bias. While single-breed evaluations tended to have lower accuracy, merging all breeds in multi-breed evaluations increased accuracy and reduced dispersion. This study demonstrates that multi-breed genomic evaluations are proper for indicine beef cattle breeds. The MF approach may be particularly beneficial for less-represented breeds, addressing limitations related to small reference populations and incompatibilities between G and A22. By leveraging genomic information across breeds, breeders and producers can make more informed selection decisions, ultimately improving genetic gain in these cattle populations.

Project description:As crossbreeding is extensively used in some livestock species, we aimed to evaluate the performance of single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) methods to predict Genomic Estimated Breeding Values (GEBVs) of crossbred animals. Different training population scenarios were evaluated: (SC1) ssGBLUP based on a single-trait model considering purebred and crossbred animals in a joint training population; (SC2) ssGBLUP based on a multiple-trait model to enable considering phenotypes recorded in purebred and crossbred training animals as different traits; (SC3) WssGBLUP based on a single-trait model considering purebred and crossbred animals jointly in the training population (both populations were used for SNP weights' estimation); (SC4) WssGBLUP based on a single-trait model considering only purebred animals in the training population (crossbred population only used for SNP weights' estimation); (SC5) WssGBLUP based on a single-trait model and the training population characterized by purebred animals (purebred population used for SNP weights' estimation). A complex trait was simulated assuming alternative genetic architectures. Different scaling factors to blend the inverse of the genomic (G -1) and pedigree ( A 22 - 1 ) relationship matrices were also tested. The predictive performance of each scenario was evaluated based on the validation accuracy and regression coefficient. The genetic correlations across simulated populations in the different scenarios ranged from moderate to high (0.71-0.99). The scenario mimicking a completely polygenic trait ( h Q T L 2 = 0) yielded the lowest validation accuracy (0.12; for SC3 and SC4). The simulated scenarios assuming 4,500 QTLs affecting the trait and h Q T L 2 = h 2 resulted in the greatest GEBV accuracies (0.47; for SC1 and SC2). The regression coefficients ranged from 0.28 (for SC3 assuming polygenic effect) to 1.27 (for SC2 considering 4,500 QTLs). In general, SC3 and SC5 resulted in inflated GEBVs, whereas other scenarios yielded deflated GEBVs. The scaling factors used to combine G -1 and A 22 - 1 had a small influence on the validation accuracies, but a greater effect on the regression coefficients. Due to the complexity of multiple-trait models and WssGBLUP analyses, and a similar predictive performance across the methods evaluated, SC1 is recommended for genomic evaluation in crossbred populations with similar genetic structures [moderate-to-high (0.71-0.99) genetic correlations between purebred and crossbred populations].

Project description:There has been a growing interest in the genetic improvement of carcass traits as an important and primary breeding goal in the beef cattle industry over the last few decades. The use of correlated traits and molecular information can aid in obtaining more accurate estimates of breeding values. This study aimed to assess the improvement in the accuracy of genetic predictions for carcass traits by using ultrasound measurements and yearling weight along with genomic information in Hanwoo beef cattle by comparing four evaluation models using the estimators of the recently developed linear regression method. We compared the performance of single-trait pedigree best linear unbiased prediction [ST-BLUP and single-step genomic (ST-ssGBLUP)], as well as multi-trait (MT-BLUP and MT-ssGBLUP) models for the studied traits at birth and yearling date of steers. The data comprised of 15,796 phenotypic records for yearling weight and ultrasound traits as well as 5,622 records for carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score), resulting in 43,949 single-nucleotide polymorphisms from 4,284 steers and 2,332 bulls. Our results indicated that averaged across all traits, the accuracy of ssGBLUP models (0.52) was higher than that of pedigree-based BLUP (0.34), regardless of the use of single- or multi-trait models. On average, the accuracy of prediction can be further improved by implementing yearling weight and ultrasound data in the MT-ssGBLUP model (0.56) for the corresponding carcass traits compared to the ST-ssGBLUP model (0.49). Moreover, this study has shown the impact of genomic information and correlated traits on predictions at the yearling date (0.61) using MT-ssGBLUP models, which was advantageous compared to predictions at birth date (0.51) in terms of accuracy. Thus, using genomic information and high genetically correlated traits in the multi-trait model is a promising approach for practical genomic selection in Hanwoo cattle, especially for traits that are difficult to measure.

Dataset Information

A Note on the Conditioning of the H^-1 Matrix Used in Single-Step GBLUP.

Publications

A Note on the Conditioning of the <b>H</b><sup>-1</sup> Matrix Used in Single-Step GBLUP.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets