Genotype-by-gestational thermal environment interaction and its impact on the future performance of tropical composite beef cattle offspring.

With global warming, there are growing challenges for raising taurine and composite beef cattle populations in tropical regions, including elevated temperatures, limited forage availability, parasite infestation, and infectious diseases. These environmental factors can trigger specific physiological responses in the developing fetus, which may have long-term implications on its performance. Therefore, the main objective of this study was to assess the influence of naturally-induced thermal stress during the gestation period on the subsequent performance of tropical composite beef cattle progeny. Furthermore, we aimed to investigate the impact of genotype-by-gestational thermal environment interaction (G×Eg) on traits under selection pressure in the breeding population. A total of 157,414 animals from 58 farms located in various Brazilian states were recorded for birth weight (BW), pre-weaning weight gain (PWG), yearling weight (YW), hip height (HH), scrotal circumference (SC), and days to first calving (DFC). We first applied a linear regression model to the BW data, which revealed that the last 40 days of gestation were suitable for calculating the mean temperature humidity index (THIg). Subsequent regression analyses revealed that for every 10-unit increase in THIg, detrimental effects of approximately 1.13% to 16.34% are expected for all traits evaluated. Genetic parameters were estimated through a reaction norm model using THIg as the environmental descriptor. The posterior means of heritability estimates (standard deviation) were 0.35 (0.07), 0.25 (0.03), 0.31 (0.03), 0.37 (0.01), 0.29 (0.07), and 0.20 (0.09) for the direct effect of BW, PWG, YW, HH, SC, and DFC, respectively. These estimates varied along the range of THIg values, suggesting a variable response to selection depending on the thermal environment during gestation. Genetic correlation estimates between more divergent THIg values were low or negative for YW, PWG, and DFC, indicating that the best-performing individuals at low THIg values may not perform as well at high THIg values and vice versa. Overall, thermal stress during gestation impacts the future performance of beef cattle offspring. Our results indicate the need for developing effective breeding strategies that take into account G×Eg effects and the re-ranking of breeding animals along the THIg scale, particularly for traits such as DFC that are highly sensitive to thermal stress.

Transgenerational effects of the maternal gestational environment on the post-natal performance of beef cattle: A reaction norm approach.

In tropical beef cattle production systems, animals are commonly raised on pastures, exposing them to potential stressors. The end of gestation typically overlaps with a dry period characterized by limited food availability. Late gestation is pivotal for fetal development, making it an ideal scenario for inter- and transgenerational effects of the maternal gestational environment. Intergenerational effects occur due to exposure during gestation, impacting the development of the embryo and its future germline. Transgenerational effects, however, extend beyond direct exposure to the subsequent generations. The objective of the present study was to verify these effects on the post-natal performance of zebu beef cattle. We extended the use of a reaction norm model to identify genetic variation in the animals' responses to transgenerational effects. The inter- and transgenerational effects were predominantly positive (-0.09% to 19.74%) for growth and reproductive traits, indicating improved animal performance on the phenotypic scale in more favourable maternal gestational environments. Additionally, these effects were more pronounced in the reproductive performance of females. On average, the ratio of direct additive genetic variances of the slope and intercept of the reaction norm ranged from 1.23% to 3.60% for direct and from 10.17% to 11.42% for maternal effects. Despite its relatively modest magnitude, this variation proved sufficient to prompt modifications in parameter estimates. The average percentage variation of direct heritability estimates ranged from 19.3% for scrotal circumference to 33.2% for yearling weight across the environmental descriptors evaluated. Genetic correlations between distant environments for the studied traits were generally high for direct effects and far from unity for maternal effects. Changes in EBV rankings of sires across different gestational environments were also observed. Due to the multifaceted nature of inter- and transgenerational effects of the maternal gestational environment on various traits of beef cattle raised under tropical pasture conditions, they should not be overlooked by producers and breeders. There were differences in the specific response of beef cattle to variations in the quality of the maternal gestational environment, which can be partially explained by transgenerational epigenetic inheritance. Adopting a reaction norm model to capture a portion of the additive variance induced by inter- or transgenerational effects could be an alternative for future research and animal genetic evaluations.

Genotype-environment interaction for milk production of Gyr cattle in Brazil and Colombia.

BACKGROUND: Genotype by environment interactions (G × E) can play an important role in cattle populations and should be included in breeding programs in order to select the best animals for different environments. OBJECTIVE: The aim of this study was to investigate the G × E for milk production of Gyr cattle in Brazil and Colombia by applying a reaction norm model used genomics information, and to identify genomic regions associated with milk production in the two countries. METHODS: The Brazilian and Colombian database included 464 animals (273 cows and 33 sires from Brazil and 158 cows from Colombia) and 27,505 SNPs. A two-trait animal model was used for milk yield adjusted to 305 days in Brazil and Colombia as a function of country of origin, which included genomic information obtained with a single-step genomic reaction norm model. The GIBBS3F90 and POSTGSf90 programs were used. RESULTS: The results obtained indicate G × E based on the reranking of bulls between Brazil and Colombia, demonstrating environmental differences between the two countries. The findings highlight the importance of considering the environment when choosing breeding animals in order to ensure the adequate performance of their progeny. Within this context, the reranking of bulls and the different SNPs associated with milk production in the two countries suggest that G × E is an important effect that should be included in the genetic evaluation of Dairy Gyr cattle in Brazil and Colombia. CONCLUSION: The Gyr breeding program can be optimized by choosing a selection environment that will allow maximum genetic progress in milk production in different environments within and between countries.

Status of the genetic diversity and population structure of the Pêga donkey.

Pedigree analysis was extended to the Pêga donkey population in order to evaluate the status of genetic diversity and population structure. All parameters were computed for three reference populations of animals born between 2004 and 2014: all animals, animals with mouse gray coat color, and animals with roan coat color. Although the average inbreeding coefficient observed was low (about 3 %), highly inbred animals are present in the current population. The effective population size based on the individual inbreeding rate was 35, while the effective population size based on the individual coancestry rate was about three times higher. The number of equivalent subpopulations was at around three, indicating that the Pêga donkey population is highly structured. There is no evidence of differentiation between subpopulations based on the coat color of the animals (Nei's minimum distance 0.10 %). The breeding policy of Pêga donkeys is predominantly intra-herd. The loss of genetic diversity since the founder generations can be considered small in Pêga donkeys (1.25 %). The excessive contribution of few ancestors to the gene pool may lead to narrower bottlenecks in the pedigree of this population in the future. The long generation interval in Pêga donkeys (10.7 years) may be considered an advantage to reduce the increase in inbreeding and to maintain the genetic diversity of these animals.

Reaction norm model to describe environmental sensitivity across first lactation in dairy cattle under tropical conditions.

Reaction norm models have been widely used to study genotype by environment interaction (G × E) in animal breeding. The objective of this study was to describe environmental sensitivity across first lactation in Brazilian Holstein cows using a reaction norm approach. A total of 50,168 individual monthly test day (TD) milk yields (10 test days) from 7476 complete first lactations of Holstein cattle were analyzed. The statistical models for all traits (10 TDs and for 305-day milk yield) included the fixed effects of contemporary group, age of cow (linear and quadratic effects), and days in milk (linear effect), except for 305-day milk yield. A hierarchical reaction norm model (HRNM) based on the unknown covariate was used. The present study showed the presence of G × E in milk yield across first lactation of Holstein cows. The variation in the heritability estimates implies differences in the response to selection depending on the environment where the animals of this population are evaluated. In the average environment, the heritabilities for all traits were rather similar, in range from 0.02 to 0.63. The scaling effect of G × E predominated throughout most of lactation. Particularly during the first 2 months of lactation, G × E caused reranking of breeding values. It is therefore important to include the environmental sensitivity of animals according to the phase of lactation in the genetic evaluations of Holstein cattle in tropical environments.

Multi-trait linear reaction norm model to describe the pattern of phenotypic expression of some economic traits in beef cattle across a range of environments.

The multi-trait reaction norm (MTRN) model was extended to beef cattle reared under tropical conditions with the following objectives: to compare multi-trait (MT) and MTRN models regarding the genetic parameters obtained; and to characterize G × E, the pattern of phenotypic expression, and the environmental sensitivity of animals for postweaning weight gain (PWG), scrotal circumference (SC), and annual average productivity of the cow (PRODAM). There was divergence in the estimates between the MT and MTRN models when the posterior probability intervals of additive genetic variances and heritability coefficients of PWG and PRODAM were analyzed. The MTRN model indicated an increase in heritability for PWG and PRODAM with improvement of the environmental conditions. For SC, heritability was practically the same, irrespective of the environmental conditions. The genetic correlations between the traits studied were low but varied over environments by the MTRN model. Considering genetic correlations obtained by the MTRN model for the same trait, lower estimates were obtained between extreme favorable and unfavorable environments. This finding suggest re-ranking of breeding values in different environments mainly for PWG and PRODAM. Thus, G × E is more important for PWG and PRODAM than for SC and should be included in the genetic evaluation of these traits. The traits PWG and PRODAM can be considered plastic traits, whereas SC is poorly plastic. The genetic trends in individual animal slopes indicate that the population is moving towards greater plasticity. This could be a matter of concern for breeders since greater plasticity seems to limit heritability and, consequently, the responses to selection.

Inclusion of weaning management group as a random effect in the genetic evaluation of postweaning traits in Nellore cattle.

The objective of this study was to evaluate the impact of including weaning management group (WMG) as an uncorrelated random effect in the genetic evaluation of postweaning traits. Data from 186,507 Nellore animals' sons to 1,734 sires and 75,016 cows were analyzed. Three single-trait models were studied. These models included the contemporary group (CG) as a fixed effect and age of animal at measurement and age of dam at calving as covariates, in addition to the direct additive breeding value as random effect. The CGs for postweaning traits varied between models which included or not WMG as part of the concatenation of fixed effects. In the model in which WMG was not part of the CG, the trait was included as an uncorrelated random effect. The results suggest that although no significant effects were observed on genetic parameter estimates or animal ranking, the inclusion of WMG as an uncorrelated random effect increased the number of observations per CG and contributed to maintaining animals in the analysis that would be discarded because they were in a CG with a small number of observations. This model could therefore be recommended for the genetic evaluation of this Nellore population.

Genetic association between milk yield, stayability, and mastitis in Holstein cows under tropical conditions.

The objective of this study was to estimate genetic parameters for milk yield, stayability, and the occurrence of clinical mastitis in Holstein cows, as well as studying the genetic relationship between them, in order to provide subsidies for the genetic evaluation of these traits. Records from 5,090 Holstein cows with calving varying from 1991 to 2010, were used in the analysis. Two standard multivariate analyses were carried out, one containing the trait of accumulated 305-day milk yields in the first lactation (MY1), stayability (STAY) until the third lactation, and clinical mastitis (CM), as well as the other traits, considering accumulated 305-day milk yields (Y305), STAY, and CM, including the first three lactations as repeated measures for Y305 and CM. The covariance components were obtained by a Bayesian approach. The heritability estimates obtained by multivariate analysis with MY1 were 0.19, 0.28, and 0.13 for MY1, STAY, and CM, respectively, whereas using the multivariate analysis with the Y305, the estimates were 0.19, 0.31, and 0.14, respectively. The genetic correlations between MY1 and STAY, MY1 and CM, and STAY and CM, respectively, were 0.38, 0.12, and -0.49. The genetic correlations between Y305 and STAY, Y305 and CM, and STAY and CM, respectively, were 0.66, -0.25, and -0.52.