Strategies for genomic predictions of an indicine multi-breed population using single-step GBLUP.

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.

Reliabilities of estimated breeding values in models with metafounders.

BACKGROUND: Reliabilities of best linear unbiased predictions (BLUP) of breeding values are defined as the squared correlation between true and estimated breeding values and are helpful in assessing risk and genetic gain. Reliabilities can be computed from the prediction error variances for models with a single base population but are undefined for models that include several base populations and when unknown parent groups are modeled as fixed effects. In such a case, the use of metafounders in principle enables reliabilities to be derived. METHODS: We propose to compute the reliability of the contrast of an individual's estimated breeding value with that of a metafounder based on the prediction error variances of the individual and the metafounder, their prediction error covariance, and their genetic relationship. Computation of the required terms demands only little extra work once the sparse inverse of the mixed model equations is obtained, or they can be approximated. This also allows the reliabilities of the metafounders to be obtained. We studied the reliabilities for both BLUP and single-step genomic BLUP (ssGBLUP), using several definitions of reliability in a large dataset with 1,961,687 dairy sheep and rams, most of which had phenotypes and among which 27,000 rams were genotyped with a 50K single nucleotide polymorphism (SNP) chip. There were 23 metafounders with progeny sizes between 100,000 and 2000 individuals. RESULTS: In models with metafounders, directly using the prediction error variance instead of the contrast with a metafounder leads to artificially low reliabilities because they refer to a population with maximum heterozygosity. When only one metafounder is fitted in the model, the reliability of the contrast is shown to be equivalent to the reliability of the individual in a model without metafounders. When there are several metafounders in the model, using a contrast with the oldest metafounder yields reliabilities that are on a meaningful scale and very close to reliabilities obtained from models without metafounders. The reliabilities using contrasts with ssGBLUP also resulted in meaningful values. CONCLUSIONS: This work provides a general method to obtain reliabilities for both BLUP and ssGBLUP when several base populations are included through metafounders.

Lifting COVID-19 mitigation measures in Spain (May-June 2020).

Introduction: The state of alarm was declared in Spain due to the COVID-19 epidemic on March 14, 2020, and established population confinement measures. The objective is to describe the process of lifting these mitigation measures. Methods: The Plan for the Transition to a New Normality, approved on April 28, contained four sequential phases with progressive increase in socio-economic activities and population mobility. In parallel, a new strategy for early diagnosis, surveillance and control was implemented. A bilateral decision mechanism was established between the Spanish Government and the autonomous communities (AC), guided by a set of qualitative and quantitative indicators capturing the epidemiological situation and core capacities. The territorial units were established ad-hoc and could be from Basic Health Zones to entire AC. Results: The process run from May 4 to June 21, 2020. AC implemented plans for reinforcement of core capacities. Incidence decreased from a median (50% of territories) of 7.4 per 100,000 in 7 days at the beginning to 2.5 at the end. Median PCR testing increased from 53% to 89% of suspected cases and PCR total capacity from 4.5 to 9.8 per 1000 inhabitants weekly; positivity rate decreased from 3.5% to 1.8%. Median proportion of cases with traced contacts increased from 82% to 100%. Conclusion: Systematic data collection, analysis, and interterritorial dialogue allowed adequate process control. The epidemiological situation improved but, mostly, the process entailed a great reinforcement of core response capacities nation-wide, under common criteria. Maintaining and further reinforcing capacities remained crucial for responding to future waves.

Introducción: El 14 de marzo de 2020 España declaró el estado de alarma por la pandemia por COVID-19 incluyendo medidas de confinamiento. El objetivo es describir el proceso de desescalada de estas medidas. Métodos: Un plan de transición hacia una nueva normalidad, del 28 de abril, incluía 4 fases secuenciales incrementando progresivamente las actividades socioeconómicas y la movilidad. Concomitantemente, se implementó una nueva estrategia de diagnóstico precoz, vigilancia y control. Se estableció un mecanismo de decisión bilateral entre Gobierno central y comunidades autónomas (CCAA), guiado por un panel de indicadores cualitativos y cuantitativos de la situación epidemiológica y las capacidades básicas. Las unidades territoriales evaluadas comprendían desde zonas básicas de salud hasta CCAA. Resultados: El proceso se extendió del 4 de mayo al 21 de junio y se asoció a planes de refuerzo de las capacidades en las CCAA. La incidencia disminuyó de una mediana inicial de 7,4 por 100.000 en 7 días a 2,5 al final del proceso. La mediana de pruebas PCR aumentó del 53% al 89% de los casos sospechosos, y la capacidad total de 4,5 a 9,8 pruebas semanales por 1.000 habitantes; la positividad disminuyó del 3,5% al 1,8%. La mediana de casos con contactos trazados aumentó del 82% al 100%. Conclusión: La recogida y análisis sistemático de información y el diálogo interterritorial logaron un adecuado control del proceso. La situación epidemiológica mejoró, pero sobre todo, se aumentaron las capacidades, en todo el país y con criterios comunes, cuyo mantenimiento y refuerzo fue clave en olas sucesivas.

Prediction ability of an alternative multi-trait genomic evaluation for residual feed intake.

Selection for feed efficiency is the goal for many genetic breeding programs in beef cattle. Residual feed intake has been included in genetic evaluations to reduce feed intake without compromising performance traits as liveweight, body gain or carcass traits. However, measuring feed intake is expensive, and only a small percentage of selection candidates are phenotyped. Genomic selection has become a very important tool to achieve effective genetic progress in these traits. Another effective strategy has been the implementation of multi-trait prediction using easily recordable predictor traits on both reference animals and candidates without phenotypes, and this could be another inexpensive way to increase accuracy. The objective of this work was to analyse and compare the prediction ability of two alternative different approaches to predict GEBVs for RFI. The population of inference was Hereford bulls in Uruguay that were genotyped candidates for to selection. The first model was the conventional univariate model for RFI and the second model was a multi-trait model which included a predictor trait (weaning weight, WW), in addition to the traits used in the first one (dry matter intake, metabolic mid test weight, average daily gain and ultrasound back fat) (DMI, MWT, ADG, UBF, respectively). GEBVs from the multi-trait model were combined using selection index theory to derive RFI values. All analyses were performed using ssGBLUP procedure. The prediction ability of both models was tested using two validation strategies (30 different replicates of random groups of animals and validation across 9 different feed intake tests). The prediction quality was assessed by the following parameters: bias, dispersion, ratio of accuracies and the relative increase in accuracy by adding phenotypic information. All parameters showed that the univariate model outperforms the multi-trait model, regardless of the validation strategy considered. These results indicate that including WW as a proxy trait in a multi-trait analysis does not improve the prediction ability when all animals to be predicted are genotyped.

Implication of the order of blending and tuning when computing the genomic relationship matrix in single-step GBLUP.

Single-step genomic BLUP (ssGBLUP) relies on the combination of the genomic ( G $$ \mathbf{G} $$ ) and pedigree relationship matrices for all ( A $$ \mathbf{A} $$ ) and genotyped ( A 22 $$ {\mathbf{A}}_{22} $$ ) animals. The procedure ensures G $$ \mathbf{G} $$ and A 22 $$ {\mathbf{A}}_{22} $$ are compatible so that both matrices refer to the same genetic base ('tuning'). Then G $$ \mathbf{G} $$ is combined with a proportion of A 22 $$ {\mathbf{A}}_{22} $$ ('blending') to avoid singularity problems and to account for the polygenic component not accounted for by markers. This computational procedure has been implemented in the reverse order (blending before tuning) following the sequential research developments. However, blending before tuning may result in less optimal tuning because the blended matrix already contains a proportion of A 22 $$ {\mathbf{A}}_{22} $$ . In this study, the impact of 'tuning before blending' was compared with 'blending before tuning' on genomic estimated breeding values (GEBV), single nucleotide polymorphism (SNP) effects and indirect predictions (IP) from ssGBLUP using American Angus Association and Holstein Association USA, Inc. data. Two slightly different tuning methods were used; one that adjusts the mean diagonals and off-diagonals of G $$ \mathbf{G} $$ to be similar to those in A 22 $$ {\mathbf{A}}_{22} $$ and another one that adjusts based on the average difference between all elements of G $$ \mathbf{G} $$ and A 22 $$ {\mathbf{A}}_{22} $$ . Over 6 million Angus growth records and 5.9 million Holstein udder depth records were available. Genomic information was available on 51,478 Angus and 105,116 Holstein animals. Average realized relationship estimates among groups of animals were similar across scenarios. Scatterplots show that GEBV, SNP effects and IP did not noticeably change for all animals in the evaluation regardless of the order of computations and when using blending parameter of 0.05. Formulas were derived to determine the blending parameter that maximizes changes in the genomic relationship matrix and GEBV when changing the order of blending and tuning. Algebraically, the change is maximized when the blending parameter is equal to 0.5. Overall, tuning G $$ \mathbf{G} $$ before blending, regardless of blending parameter used, had a negligible impact on genomic predictions and SNP effects in this study.

Genomic evaluation of commercial herds with different pedigree structures using the single-step genomic BLUP in Nelore cattle.

The aim of this work was to evaluate the impact of applying genomic information in pedigree uncertainty situations on genetic evaluations for growth- and cow productivity-related traits in Nelore commercial herds. Records for accumulated cow productivity (ACP) and adjusted weight at 450 days of age (W450) were used, as well as genotypes of registered and commercial herd animals, genotyped with the Clarifide Nelore 3.1 panel (~29,000 SNPs). The genetic values for commercial and registered populations were estimated using different approaches that included (ssGBLUP) or did not include genomic information (BLUP), with different pedigree structures. Different scenarios were tested, varying the proportion of young animals with unknown sires (0, 25, 50, 75, and 100%), and unknown maternal grandsires (0, 25, 50, 75, and 100%). The prediction accuracies and abilities were calculated. The estimated breeding value accuracies decreased as the proportion of unknown sires and maternal grandsires increased. The genomic estimated breeding value accuracy using the ssGBLUP was higher in scenarios with a lower proportion of known pedigree when compared to the BLUP methodology. The results obtained with the ssGBLUP showed that it is possible to obtain reliable direct and indirect predictions for young animals from commercial herds without pedigree structure.

Theoretical accuracy for indirect predictions based on SNP effects from single-step GBLUP.

BACKGROUND: Although single-step GBLUP (ssGBLUP) is an animal model, SNP effects can be backsolved from genomic estimated breeding values (GEBV). Predicted SNP effects allow to compute indirect prediction (IP) per individual as the sum of the SNP effects multiplied by its gene content, which is helpful when the number of genotyped animals is large, for genotyped animals not in the official evaluations, and when interim evaluations are needed. Typically, IP are obtained for new batches of genotyped individuals, all of them young and without phenotypes. Individual (theoretical) accuracies for IP are rarely reported, but they are nevertheless of interest. Our first objective was to present equations to compute individual accuracy of IP, based on prediction error covariance (PEC) of SNP effects, and in turn, are obtained from PEC of GEBV in ssGBLUP. The second objective was to test the algorithm for proven and young (APY) in PEC computations. With large datasets, it is impossible to handle the full PEC matrix, thus the third objective was to examine the minimum number of genotyped animals needed in PEC computations to achieve IP accuracies that are equivalent to GEBV accuracies. RESULTS: Correlations between GEBV and IP for the validation animals using SNP effects from ssGBLUP evaluations were ≥ 0.99. When all available genotyped animals were used for PEC computations, correlations between GEBV and IP accuracy were ≥ 0.99. In addition, IP accuracies were compatible with GEBV accuracies either with direct inversion of the genomic relationship matrix (G) or using the algorithm for proven and young (APY) to obtain the inverse of G. As the number of genotyped animals included in the PEC computations decreased from around 55,000 to 15,000, correlations were still ≥ 0.96, but IP accuracies were biased downwards. CONCLUSIONS: Theoretical accuracy of indirect prediction can be successfully obtained by computing SNP PEC out of GEBV PEC from ssGBLUP equations using direct or APY G inverse. It is possible to reduce the number of genotyped animals in PEC computations, but accuracies may be underestimated. Further research is needed to approximate SNP PEC from ssGBLUP to limit the computational requirements with many genotyped animals.

Certification program of Addiction Centres for hepatitis C virus elimination in Spain. HepCelentes Project. / Programa de certificación de Centros de Adicciones para la eliminación del virus de la hepatitis C en España. Proyecto HepCelentes.

Microelimination strategies for the hepatitis C virus (HCV) in vulnerable populations, such as users of Addiction Centres (AC), are key for the elimination of hepatitis C. The aim of the HepCelentes project was to design a certification program for AC from the generation of a guide with the criteria to favour the prevention, diagnosis, control, and treatment of HCV in Spain. The project was structured in 4 phases: normalisation, implementation, certification, and communication. In the first phase, developed between July and December 2020, a Steering Committee was created (formed by representatives of scientific societies, healthcare professionals from AC, primary care centres and hospital units, and patient associations) that, from of an exhaustive bibliographic review, generated by consensus an accreditation guide for AC. The guide consists of 22 criteria (15 mandatory and 7 recommended) structured based on the requirements to be met by AC, justification for the selection, level of action (management, prevention, diagnosis and treatment/follow-up), measurement of the indicator, objective level to be achieved, evidence of compliance, clarifications to improve understanding, and mandatory / recommendation (depending on their relevance to achieve HCV elimination and its feasibility for implementation in real practice). The development of a certification system for the AC, based on consensus and coordination of multidisciplinary teams, is intended to favour the management of hepatitis C and its elimination in AC users, supporting the international, national, and regional elimination strategies.

Las estrategias de microeliminación del virus de la hepatitis C (VHC) en poblaciones vulnerables, como los usuarios de los centros de adicciones (CA), son fundamentales para lograr la eliminación de la hepatitis C. El objetivo del proyecto HepCelentes fue diseñar un programa de certificación para los CA, a partir de la generación de una guía con los criterios para favorecer la prevención, diagnóstico, control y tratamiento del VHC en España. El proyecto se estructuró en 4 fases: normalización, implementación, certificación y comunicación. En la primera fase, desarrollada entre julio y diciembre de 2020, se creó un Comité de Normalización (formado por representantes de sociedades científicas, profesionales sanitarios de CA, centros de atención primaria, unidades hospitalarias, y asociaciones de pacientes) que, a partir de una revisión bibliográfica exhaustiva, generó por consenso una guía de certificación de los CA. La guía consta de 22 criterios (15 obligatorios y 7 recomendados) estructurados en base a la definición del criterio, justificación de su selección, nivel de actuación (gestión, prevención, diagnóstico y tratamiento/seguimiento), fórmula de medición, nivel objetivo a alcanzar, evidencias de su cumplimiento, aclaraciones para mejorar su comprensión y obligatoriedad/recomendación (en función de la relevancia en la eliminación y capacidad de implementación). El desarrollo de un sistema de certificación para los CA, a partir del consenso y la coordinación de equipos multidisciplinares, pretende favorecer el manejo de la hepatitis C y su eliminación en los usuarios de los CA, apoyando las estrategias de eliminación internacionales, nacionales y autonómicas.

Assessing the accuracy of imputation in the Gyr breed using different SNP panels.

The aim of this study was to evaluate the accuracy of imputation in a Gyr population using two medium-density panels (Bos taurus - Bos indicus) and to test whether the inclusion of the Nellore breed increases the imputation accuracy in the Gyr population. The database consisted of 289 Gyr females from Brazil genotyped with the GGP Bovine LDv4 chip containing 30 000 SNPs and 158 Gyr females from Colombia genotyped with the GGP indicus chip containing 35 000 SNPs. A customized chip was created that contained the information of 9109 SNPs (9K) to test the imputation accuracy in Gyr populations; 604 Nellore animals with information of LD SNPs tested in the scenarios were included in the reference population. Four scenarios were tested: LD9K_30KGIR, LD9K_35INDGIR, LD9K_30KGIR_NEL, and LD9K_35INDGIR_NEL. Principal component analysis (PCA) was computed for the genomic matrix and sample-specific imputation accuracies were calculated using Pearson's correlation (CS) and the concordance rate (CR) for imputed genotypes. The results of PCA of the Colombian and Brazilian Gyr populations demonstrated the genomic relationship between the two populations. The CS and CR ranged from 0.88 to 0.94 and from 0.93 to 0.96, respectively. Among the scenarios tested, the highest CS (0.94) was observed for the LD9K_30KGIR scenario. The present results highlight the importance of the choice of chip for imputation in the Gyr breed. However, the variation in SNPs may reduce the imputation accuracy even when the chip of the Bos indicus subspecies is used.

Reproductive performance of Angus, Hereford, Salers and Nellore crossbred females: Additive and non-additive effects.

Reproductive traits in breeding herds can be improved through crossbreeding, which results in breed differences, heterosis and breed complementarity. The aim of this study was to estimate group additive genetic and dominance effects for three reproductive traits; probability of artificial insemination (AIP); calving success (CS); and days to calving (DC) for Hereford (H), Angus (A), Nellore (N) and Salers (S) breeds under grazing conditions. Data were obtained from an experiment carried out during 1992-2002 by the Faculty of Agronomy, Universidad de la Republica (UdelaR), Uruguay and Caja Notarial de Seguridad Social. The data set contained reproductive information of 1,164 females from 11 different genetic groups (GG) consisting of crosses between H, N, S and A. AIP, CS and DC were examined in first-calf heifers, while CS and DC were examined in second-calf and 3- to 7-year-old cows. Least square means for each GG and group additive genetic and dominance effects were estimated for each trait. F1 crossbreed females performed better for artificial insemination probability than purebred females. Crossbred A/H heifers had the highest AIPs and CS rates, while crossbred N/H 3- to 7-year-old cows recorded the highest averages for CS and DC. Estimates of group additive genetic effects did not differ amongst A, S, N and H; however, dominance increased the AIP and CS of the heifers.

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.

Bias and accuracy of dairy sheep evaluations using BLUP and SSGBLUP with metafounders and unknown parent groups.

BACKGROUND: Bias has been reported in genetic or genomic evaluations of several species. Common biases are systematic differences between averages of estimated and true breeding values, and their over- or under-dispersion. In addition, comparing accuracies of pedigree versus genomic predictions is a difficult task. This work proposes to analyse biases and accuracies in the genetic evaluation of milk yield in Manech Tête Rousse dairy sheep, over several years, by testing five models and using the estimators of the linear regression method. We tested models with and without genomic information [best linear unbiased prediction (BLUP) and single-step genomic BLUP (SSGBLUP)] and using three strategies to handle missing pedigree [unknown parent groups (UPG), UPG with QP transformation in the [Formula: see text] matrix (EUPG) and metafounders (MF)]. METHODS: We compared estimated breeding values (EBV) of selected rams at birth with the EBV of the same rams obtained each year from the first daughters with phenotypes up to 2017. We compared within and across models. Finally, we compared EBV at birth of the rams with and without genomic information. RESULTS: Within models, bias and over-dispersion were small (bias: 0.20 to 0.40 genetic standard deviations; slope of the dispersion: 0.95 to 0.99) except for model SSGBLUP-EUPG that presented an important over-dispersion (0.87). The estimates of accuracies confirm that the addition of genomic information increases the accuracy of EBV in young rams. The smallest bias was observed with BLUP-MF and SSGBLUP-MF. When we estimated dispersion by comparing a model with no markers to models with markers, SSGBLUP-MF showed a value close to 1, indicating that there was no problem in dispersion, whereas SSGBLUP-EUPG and SSGBLUP-UPG showed a significant under-dispersion. Another important observation was the heterogeneous behaviour of the estimates over time, which suggests that a single check could be insufficient to make a good analysis of genetic/genomic evaluations. CONCLUSIONS: The addition of genomic information increases the accuracy of EBV of young rams in Manech Tête Rousse. In this population that has missing pedigrees, the use of UPG and EUPG in SSGBLUP produced bias, whereas MF yielded unbiased estimates, and we recommend its use. We also recommend assessing biases and accuracies using multiple truncation points, since these statistics are subject to random variation across years.

Effects of ignoring inbreeding in model-based accuracy for BLUP and SSGBLUP.

Model-based accuracy, defined as the theoretical correlation between true and estimated breeding value, can be obtained for each individual as a function of its prediction error variance (PEV) and inbreeding coefficient F, in BLUP, GBLUP and SSGBLUP genetic evaluations. However, for computational convenience, inbreeding is often ignored in two places. First, in the computation of reliability = 1-PEV/(1 + F). Second, in the set-up, using Henderson's rules, of the inverse of the pedigree-based relationship matrix A. Both approximations have an effect in the computation of model-based accuracy and result in wrong values. In this work, first we present a reminder of the theory and extend it to SSGBLUP. Second, we quantify the error of ignoring inbreeding with real data in three scenarios: BLUP evaluation and SSGBLUP in Uruguayan dairy cattle, and BLUP evaluations in a line of rabbit closed for >40 generations with steady increase of inbreeding up to an average of 0.30. We show that ignoring inbreeding in the set-up of the A-inverse is equivalent to assume that non-inbred animals are actually inbred. This results in an increase of apparent PEV that is negligible for dairy cattle but considerable for rabbit. Ignoring inbreeding in reliability = 1-PEV/(1 + F) leads to underestimation of reliability for BLUP evaluations, and this underestimation is very large for rabbit. For SSGBLUP in dairy cattle, it leads to both underestimation and overestimation of reliability, both for genotyped and non-genotyped animals. We strongly recommend to include inbreeding both in the set-up of A-inverse and in the computation of reliability from PEVs.

Autoregressive and random regression test-day models for multiple lactations in genetic evaluation of Brazilian Holstein cattle.

Autoregressive (AR) and random regression (RR) models were fitted to test-day records from the first three lactations of Brazilian Holstein cattle with the objective of comparing their efficiency for national genetic evaluations. The data comprised 4,142,740 records of milk yield (MY) and somatic cell score (SCS) from 274,335 cows belonging to 2,322 herds. Although heritabilities were similar between models and traits, additive genetic variance estimates using AR were 7.0 (MY) and 22.2% (SCS) higher than those obtained from RR model. On the other hand, residual variances were lower in both traits when estimated through AR model. The rank correlation between EBV obtained from AR and RR models was 0.96 and 0.94 (MY) and 0.97 and 0.95 (SCS), respectively, for bulls (with 10 or more daughters) and cows. Estimated annual genetic gains for bulls (cows) obtained using AR were 46.11 (49.50) kg for MY and -0.019 (-0.025) score for SCS; whereas using RR these values were 47.70 (55.56) kg and -0.022 (-0.028) score. Akaike information criterion was lower for AR in both traits. Although AR model is more parsimonious, RR model assumes genetic correlations different from the unity within and across lactations. Thus, when these correlations are relatively high, these models tend to yield to similar predictions; otherwise, they will differ more and RR model would be theoretically sounder.

Frequentist p-values for large-scale-single step genome-wide association, with an application to birth weight in American Angus cattle.

BACKGROUND: Single-step genomic best linear unbiased prediction (SSGBLUP) is a comprehensive method for genomic prediction. Point estimates of marker effects from SSGBLUP are often used for genome-wide association studies (GWAS) without a formal framework of hypothesis testing. Our objective was to implement p-values for single-marker GWAS studies within the single-step GWAS (SSGWAS) framework by deriving computational algorithms and procedures, and by applying these to a large beef cattle population. METHODS: P-values were obtained based on the prediction error (co)variances for single nucleotide polymorphisms (SNPs), which were obtained from the prediction error (co)variances of genomic predictions based on the inverse of the coefficient matrix and formulas to estimate SNP effects. RESULTS: Computation of p-values took a negligible time for a dataset with almost 2 million animals in the pedigree and 1424 genotyped sires, and no inflation of statistics was observed. The SNPs that passed the Bonferroni threshold of 10-5.9 were the same as those that explained the highest proportion of additive genetic variance, but even at the same significance levels and effects, some of them explained less genetic variance due to lower allele frequency. CONCLUSIONS: The use of a p-value for SSGWAS is a very general and efficient strategy to identify quantitative trait loci (QTL). It can be used for complex datasets such as those used in animal breeding, where only a proportion of the pedigreed animals are genotyped.

Resource allocation optimization with multi-trait genomic prediction for bread wheat (Triticum aestivum L.) baking quality.

KEY MESSAGE: Multi-trait genomic prediction models are useful to allocate available resources in breeding programs by targeted phenotyping of correlated traits when predicting expensive and labor-intensive quality parameters. Multi-trait genomic prediction models can be used to predict labor-intensive or expensive correlated traits where phenotyping depth of correlated traits could be larger than phenotyping depth of targeted traits, reducing resources and improving prediction accuracy. This is particularly important in the context of allocating phenotyping resource in plant breeding programs. The objective of this work was to evaluate multi-trait models predictive ability with different depth of phenotypic information from correlated traits. We evaluated 495 wheat advanced breeding lines for eight baking quality traits which were genotyped with genotyping-by-sequencing. Through different approaches for cross-validation, we evaluated the predictive ability of a single-trait model and a multi-trait model. Moreover, we evaluated different sizes of the training population (from 50 to 396 individuals) for the trait of interest, different depth of phenotypic information for correlated traits (50 and 100%) and the number of correlated traits to be used (one to three). There was no loss in the predictive ability by reducing the training population up to a 30% (149 individuals) when using correlated traits. A multi-trait model with one highly correlated trait phenotyped for both the training and testing sets was the best model considering phenotyping resources and the gain in predictive ability. The inclusion of correlated traits in the training and testing lines is a strategic approach to replace phenotyping of labor-intensive and high cost traits in a breeding program.

Genome-wide association between single nucleotide polymorphisms with beef fatty acid profile in Nellore cattle using the single step procedure.

BACKGROUND: Saturated fatty acids can be detrimental to human health and have received considerable attention in recent years. Several studies using taurine breeds showed the existence of genetic variability and thus the possibility of genetic improvement of the fatty acid profile in beef. This study identified the regions of the genome associated with saturated, mono- and polyunsaturated fatty acids, and n-6 to n-3 ratios in the Longissimus thoracis of Nellore finished in feedlot, using the single-step method. RESULTS: The results showed that 115 windows explain more than 1 % of the additive genetic variance for the 22 studied fatty acids. Thirty-one genomic regions that explain more than 1 % of the additive genetic variance were observed for total saturated fatty acids, C12:0, C14:0, C16:0 and C18:0. Nineteen genomic regions, distributed in sixteen different chromosomes accounted for more than 1 % of the additive genetic variance for the monounsaturated fatty acids, such as the sum of monounsaturated fatty acids, C14:1 cis-9, C18:1 trans-11, C18:1 cis-9, and C18:1 trans-9. Forty genomic regions explained more than 1 % of the additive variance for the polyunsaturated fatty acids group, which are related to the total polyunsaturated fatty acids, C20:4 n-6, C18:2 cis-9 cis12 n-6, C18:3 n-3, C18:3 n-6, C22:6 n-3 and C20:3 n-6 cis-8 cis-11 cis-14. Twenty-one genomic regions accounted for more than 1 % of the genetic variance for the group of omega-3, omega-6 and the n-6:n-3 ratio. CONCLUSIONS: The identification of such regions and the respective candidate genes, such as ELOVL5, ESSRG, PCYT1A and genes of the ABC group (ABC5, ABC6 and ABC10), should contribute to form a genetic basis of the fatty acid profile of Nellore (Bos indicus) beef, contributing to better selection of the traits associated with improving human health.

Genotype by environment interaction for tick resistance of Hereford and Braford beef cattle using reaction norm models.

BACKGROUND: The cattle tick is a parasite that adversely affects livestock performance in tropical areas. Although countries such as Australia and Brazil have developed genetic evaluations for tick resistance, these evaluations have not considered genotype by environment (G*E) interactions. Genetic gains could be adversely affected, since breedstock comparisons are environmentally dependent on the presence of G*E interactions, particularly if residual variability is also heterogeneous across environments. The objective of this study was to infer upon the existence of G*E interactions for tick resistance of cattle based on various models with different assumptions of genetic and residual variability. METHODS: Data were collected by the Delta G Connection Improvement program and included 10,673 records of tick counts on 4363 animals. Twelve models, including three traditional animal models (AM) and nine different hierarchical Bayesian reaction norm models (HBRNM), were investigated. One-step models that jointly estimate environmental covariates and reaction norms and two-step models based on previously estimated environmental covariates were used to infer upon G*E interactions. Model choice was based on the deviance criterion information. RESULTS: The best-fitting model specified heterogeneous residual variances across 10 subclasses that were bounded by every decile of the contemporary group (CG) estimates of tick count effects. One-step models generally had the highest estimated genetic variances. Heritability estimates were normally higher for HBRNM than for AM. One-step models based on heterogeneous residual variances also usually led to higher heritability estimates. Estimates of repeatability varied along the environmental gradient (ranging from 0.18 to 0.45), which implies that the relative importance of additive and permanent environmental effects for tick resistance is influenced by the environment. Estimated genetic correlations decreased as the tick infestation level increased, with negative correlations between extreme environmental levels, i.e., between more favorable (low infestation) and harsh environments (high infestation). CONCLUSIONS: HBRNM can be used to describe the presence of G*E interactions for tick resistance in Hereford and Braford beef cattle. The preferred model for the genetic evaluation of this population for tick counts in Brazilian climates was a one-step model that considered heteroscedastic residual variance. Reaction norm models are a powerful tool to identify and quantify G*E interactions and represent a promising alternative for genetic evaluation of tick resistance, since they are expected to lead to greater selection efficiency and genetic progress.

Accuracy of estimated breeding values with genomic information on males, females, or both: an example on broiler chicken.

BACKGROUND: As more and more genotypes become available, accuracy of genomic evaluations can potentially increase. However, the impact of genotype data on accuracy depends on the structure of the genotyped cohort. For populations such as dairy cattle, the greatest benefit has come from genotyping sires with high accuracy, whereas the benefit due to adding genotypes from cows was smaller. In broiler chicken breeding programs, males have less progeny than dairy bulls, females have more progeny than dairy cows, and most production traits are recorded for both sexes. Consequently, genotyping both sexes in broiler chickens may be more advantageous than in dairy cattle. METHODS: We studied the contribution of genotypes from males and females using a real dataset with genotypes on 15 723 broiler chickens. Genomic evaluations used three training sets that included only males (4648), only females (8100), and both sexes (12 748). Realized accuracies of genomic estimated breeding values (GEBV) were used to evaluate the benefit of including genotypes for different training populations on genomic predictions of young genotyped chickens. RESULTS: Using genotypes on males, the average increase in accuracy of GEBV over pedigree-based EBV for males and females was 12 and 1 percentage points, respectively. Using female genotypes, this increase was 1 and 18 percentage points, respectively. Using genotypes of both sexes increased accuracies by 19 points for males and 20 points for females. For two traits with similar heritabilities and amounts of information, realized accuracies from cross-validation were lower for the trait that was under strong selection. CONCLUSIONS: Overall, genotyping males and females improves predictions of all young genotyped chickens, regardless of sex. Therefore, when males and females both contribute to genetic progress of the population, genotyping both sexes may be the best option.