Heat stress on breeding value prediction for milk yield and composition of a Brazilian Holstein cattle population.

Due to the high milk production of Holstein cows, many countries have chosen to import semen to improve local dairy herds. This strategy would be more effective if this semen was used in the same environment conditions in which the bulls were selected. If the effect of genotype by environment (G × E) interaction is not considered, the estimated breeding values (EBVs) may vary, potentially reducing the selection response. We evaluate the impact of heat stress on selection for milk yield and composition of Holstein cows using random regression models. To verify the interference of heat stress in milk yield (MY) and composition traits (fat, protein, total saturated, and total unsaturated fatty acids content in milk), temperature-humidity index (THI) on test-day milk records was used. The threshold value to divide the environments using test-day information from Brazilian Holstein cows was 72 units of THI, i.e., < 72 represented no heat stress and > 72 represented heat stress. Legendre polynomials of second-order (Leg 2) model and two lactation points (33 and 122 DIM) were used to estimate heritabilities and EBVs for five important dairy traits. The heritabilities of milk components and fatty acids were low (0.09-0.29), regardless of lactation period and degree of heat stress, with the exception of protein content (0.30-0.35). Fat content was the only milk component that was reduced according to the degree of heat stress and lactation period. The EBVs tended to decrease in heat stress conditions, thus animals with high genetic potential demonstrated evidence of G × E interaction. However, acclimatization of dairy cows to heat stress in the farm production systems may have been responsible for the low differences among genetic parameters and EBVs with and without heat stress found in this study.

NASA POWER satellite meteorological system is a good tool for obtaining estimates of the temperature-humidity index under Brazilian conditions compared to INMET weather stations data.

Heat stress negatively affects livestock, with undesirable effects on animals' production and reproduction. Temperature and humidity index (THI) is a climatic variable used worldwide to study the effect of heat stress on farm animals. Temperature and humidity data can be obtained in Brazil through the National Institute of Meteorology (INMET), but complete data may not be available due to temporary failures on weather stations. An alternative to obtaining meteorological data is the National Aeronautics and Space Administration Prediction of Worldwide Energy Resources (NASA POWER) satellite-based weather system. We aimed to compare THI estimates obtained from INMET weather stations and NASA POWER meteorological information sources using Pearson correlation and linear regression. After quality check, data from 489 INMET weather stations were used. The hourly, average daily and maximum daily THI were evaluated. We found greater correlations and better regression evaluation metrics when average daily THI values were considered, followed by maximum daily THI, and hourly THI. NASA POWER satellite-based weather system is a suitable tool for obtaining the average and maximum THI values using information collected from Brazil, showing high correlations with THI estimates from INMET and good regression evaluation metrics, and can assist studies that aim to analyze the impact of heat stress on livestock production in Brazil, providing additional data to complement the existing information available in the INMET database.

Genomic prediction in Brazilian Guzerá cattle: application of a single-step approach to productive and reproductive traits.

This study aimed to investigate the feasibility of genomic prediction for productive and reproductive traits in Guzerá cattle using single-step genomic best linear unbiased prediction (ssGBLUP). Evaluations included the 305-day cumulative yields (first lactation, in kg) of milk, lactose, protein, fat, and total solids; adjusted body weight (kg) at the ages of 450, 365, and 210 days; and age at first calving (in days), from a database containing 197,283 measurements from Guzerá males and females born between 1954 and 2018. The pedigree included 433,823 animals spanning up to 14 overlapping generations. A total of 1618 animals were genotyped. The analyses were performed using ssGBLUP and traditional BLUP methods. Predictive ability and bias were accessed using cross-validation: predictive ability was similar between the methods and ranged from 0.27 to 0.47 for the genomic-based model and from 0.30 to 0.45 for the pedigree-based model; the bias was also similar between the methods, ranging from 0.88 to 1.35 in the genomic-based model and from 0.96 to 1.41 in the pedigree-based model. The individual accuracies of breeding values were evidently increased in the genomic evaluation, with values ranging from 0.41 to 0.56 in the genomic-based model and from 0.26 to 0.54 in the pedigree-based model. Even based on a small number of genotyped animals and a small database for some traits, the results suggest that ssGBLUP is feasible and may be applied to national genetic evaluation of the breed to increase the accuracy of breeding values without greatly impacting predictive ability and bias.

Assessment of FecGE genotypes on reproductive traits in Brazilian Morada Nova and Santa Inês sheep.

The aim of this study was to evaluate the effect of the polymorphic FecGE allele on reproductive traits in Santa Inês and Morada Nova ewes. The traits evaluated were as follows: total progeny weights at birth (PWB) and weaning (PWW) and progeny survival rates at birth (PSRB) and weaning (PSRW). A total of 389 animals, belonging to two Santa Inês herds and one Morada Nova herd, were genotyped. There was a difference between the averages for all the traits studied regarding type of parturition, herd/breed, genotype/herd, and genotype/type of parturition. For each additional progeny, if the female was FecGE/E, the PWB decreased by 1.02 kg and the PWW by 3.16 kg, also with a 0.04% reduction in PSRB and no change in PSRW. If the female was FecGE/+, the reduction in PWB was 0.24 kg, with an increase in PSRW by 0.11%, but no change in PWW and PSRB. In general, these results demonstrate that FecG+/+ females have a better ability to increase their number of progenies without reducing PWB and PWW (also similar to FecGE/+). Thus, it is suggested that further studies on the association between the traits of interest and candidate genes in sheep should be carried out so that the regions which have the greatest effect on the expression of these traits are actually identified. It was not possible to verify the effect of the FecGE allele on the PWB, PWW, PSRB, and PSRW in these Morada Nova and Santa Inês herds.

Searching for causal relationships among latent variables concerning performance, carcass, and meat quality traits in broilers.

In causal relationship studies, the latent variables may summarize the phenotypes in theoretical traits according to their phenotypic correlations, improving the understanding of causal relationships between broilers phenotypes. In this study, we aimed to investigate potential causal relationships among latent variables in broilers using a structural equation model in the context of genetic analysis. The data used in this study comprised 14 traits in broilers with 2,017 records each, and 104,154 animals in pedigree. Four latent variables (WEIGHT, LOSSES, COLOUR, and VISCERA) were defined and validated using Bayesian Confirmatory Factor Analysis. Subsequently, a search for causal linkage structures was performed, obtaining a single causal link structure between the latent variables. Then, this information was used to fit the structural equation model (SEM). The results from the SEM indicated positive causal effects of the variables WEIGHT and LOSSES on the variables VISCERA and COLOUR, respectively, with structural coefficient estimates of 1.006 and 0.040, respectively. On the other hand, an antagonist causal effect of the variable WEIGHT on the variable LOSSES was verified, with a structural coefficient estimate of -4.333. These results highlight the causal relationship between performance and meat quality traits, which may be associated with the natural processes involved in the conversion of muscle into meat and the structural changes in muscle tissues due to intense selection for high growth rates in broilers.

Genetic parameters for milk yield and quality traits of Brazilian Holstein cows as a function of temperature and humidity index.

Measurements of milk yield (MY), somatic cell score (SCS), percentage of fat (FP), protein (PP), lactose (LP), casein (CP) and percentage of palmitic (C16:0), stearic (C18:0), oleic (C18:1), total saturated (SFA), unsaturated (UFA), monounsaturated (MUFA) and polyunsaturated (PUFA) fatty acids in milk from 5,224 Holstein cows were evaluated as a function of a temperature and humidity index (THI). Legendre orthogonal polynomials from second to seventh order were tested. The best fit order for MY, PP and C18:0 was the third, whereas the second for all other traits. The heritability estimates decreased for MY (0.31 to 0.14), FP (0.28 to 0.16), LP (0.43 to 0.30), SCS (0.14 to 0.09), SFA (0.33 to 0.22) and C16:0 (0.31 to 0.26), whereas increased for CP (0.32 to 0.42), MUFA (0.08 to 0.13), UFA (0.07 to 0.11) and C18:1 (0.07 to 0.11) as the THI level increased. For PP, heritabilities (0.26 to 0.39) presented larger values in intermediate THI. For PUFA and C18:0, heritabilities were approximately constant (0.13 to 0.14 and 0.15, respectively). However, the greatest variations may have been the result of the limitations of Legendre polynomials at the extreme points of the curve, and the pattern of heritabilities curves was approximately constant for the evaluated traits. Spearman's rank correlations between breeding values in extreme THI levels were greater than 0.80 for all traits considering all animals, only cows and only bulls. When considering the top 1% and the top 50% animals (only cows, only bulls and all), Spearman correlations smaller than 0.70 were found, suggesting reranking of the animals. Although there was little variation in the variance components over THI, it is possible that there is no heat stress in the animals studied, because, on average, there was no great impact of the thermal load on the traits. One possible explanation is the use of herds with little climatic difference among herds, as well as the use of fans and sprinklers into the barns. However, the THI levels may be important factors in the selection process, as reranking of animals was verified.

Random-effect meta-analysis of genetic parameter estimates for carcass and meat quality traits in beef cattle.

Considerable variability of genetic parameter estimates is observed among different studies for the same trait, which is associated with the distinct effects included in the statistical model, population breed, and sample sizes. The random-effect meta-analysis summarizes genetic parameters considering the heterogeneity among studies. Therefore, the aim of this study was to perform a random-effect meta-analysis of heritability and genetic correlation estimates for carcass and meat quality traits in beef cattle. A total of 152 estimates of heritability and 83 genetic correlations for longissimus muscle area (LMA), back fat thickness (BFT), and marbling score (MRB) were used. High heterogeneity among published studies was observed for all traits, indicating the need of a random-effects model to perform the analysis. Estimates of heritability through the meta-analysis using the random-effects model were high (0.30 to 0.34), indicating that fast genetic progress can be obtained for these traits. However, genetic correlations had low magnitude (lower than 0.25), which suggested that all three traits should be included in the selection scheme.

Genetic parameters for milk, growth, and reproductive traits in Guzerá cattle under tropical conditions.

This study aimed to estimate the genetic parameters of milk (305-day milk yield (MY305)), growth (weaning weight (WW), yearling weight (YW), and weight at 550 days (W550)), and reproductive (age at first calving (AFC)) traits in Guzerá cattle by using Bayesian multiple-trait models. Systematic effects included sex and age at calving for the growth and milk traits, respectively. The additive genetic and contemporary groups (herd and year and season of birth) were included as random effects. Additionally, maternal genetic and permanent effects were also included as random effects for the WW. The heritability estimates were 0.29 (MY305), 0.42 (WW), 0.49 (YW), 0.56 (W550), and 0.25 (AFC). The genetic correlations among the growth traits were higher than 0.83; between the MY305 and WW, MY305 and YW, and MY305 and W550, the genetic correlations were 0.25, 0.32, and 0.36, respectively. The AFC was negatively correlated with the milk and growth traits. These results suggest the viability and potential of the joint selection for milk, beef, and reproductive traits in Guzerá cattle.

The contribution of a MOET nucleus scheme for the improvement of Guzerá (Bos indicus) cattle for milk traits in Brazil.

The Guzerá breed evolved from the introduction of breeds from India, mainly the Kankrej breed, into Brazilian livestock at the end of the 19th century. Guzerá adapted well to the climatic conditions of Brazil, where it is considered a dual-purpose breed and has been used for pasture-based beef, milk or dual-purpose production systems with the use of low-medium inputs. The importance of this genetic resource for milk production in tropical regions moved breeders to implement the National Breeding Program for the Improvement of Guzerá in 1994, based on both progeny testing and MOET nucleus schemes. We sought to evaluate the role of the MOET nucleus scheme in the phenotypic and genetic progress for milk traits in this breed. The initial database used in the present study consisted of 6,513 cows, daughters of 761 bulls. We performed genetic evaluations with different datasets using a linear mixed model in a single trait analysis, including the relationship matrix, in order to estimate breeding values. Inbreeding coefficients were also calculated using the relationship of descent between two parents. Annual phenotypic, genetic and inbreeding trends were obtained for each dataset, considering the genetic pathways of both the bull and the cow. The low genetic progress found for milk yield in the whole population (5.27 ± 0.30 kg/year) partially accounted for the dual-purpose selection goal, despite the higher genetic progress in the MOET nucleus (9.39 ± 0.79 kg/year). The inbreeding coefficient was minimized at the beginning of the breeding program based on the use of new lineages. Posteriorly, it started increasing again from 0.002 in 1991 to 0.008 in 2019. The results provided evidence of the significant contribution of the MOET nucleus scheme for the phenotypic and genetic progress of Guzerá breed for milk traits, as well as of the impact of the breeding program on the inbreeding coefficient rate in the early years. New strategies need to be designed for the Guzerá breed, to allow for greater improvement of milk traits and minimizing the rate of the inbreeding coefficient.

Genetic study of quantitative traits supports the use of Guzerá as dual-purpose cattle.

OBJECTIVE: The aim of this study was to estimate genetic parameters for 305-day cumulative milk yield and components, growth, and reproductive traits in Guzerá cattle. METHODS: The evaluated traits were 305-day first-lactation cumulative yields (kg) of milk (MY305), fat (FY305), protein (PY305), lactose (LY305), and total solids (SY305); age at first calving (AFC) in days; adjusted scrotal perimeter (cm) at the ages of 365 (SP365) and 450 (SP450) days; and adjusted body weight (kg) at the ages of 210 (W210), 365 (W365), and 450 (W450) days. The (co)variance components were estimated using the restricted maximum likelihood method for single-trait, bi-trait and tri-trait analyses. Contemporary groups and additive genetic effects were included in the general mixed model. Maternal genetic and permanent environmental effects were also included for W210. RESULTS: The direct heritability estimates ranged from 0.16 (W210) to 0.32 (MY305). The maternal heritability estimate for W210 was 0.03. Genetic correlation estimates among milk production traits and growth traits ranged from 0.92 to 0.99 and from 0.92 to 0.99, respectively. For milk production and growth traits, the genetic correlations ranged from 0.33 to 0.56. The genetic correlations among AFC and all other traits were negative (-0.43 to -0.27). Scrotal perimeter traits and body weights showed genetic correlations ranging from 0.41 to 0.46, and scrotal perimeter and milk production traits showed genetic correlations ranging from 0.11 to 0.30. The phenotypic correlations were similar in direction (same sign) and lower than the corresponding genetic correlations. CONCLUSION: These results suggest the viability and potential of joint selection for dairy and beef traits in Guzerá cattle, taking into account reproductive traits.

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.