Uncovering genetic parameters and environmental influences on fertility, milk production, and quality in autochthonous Reggiana cattle

Reggiana is a local cattle breed from northern Italy known for its rusticity and profitability, due to the production of branded Parmigiano Reggiano cheese. To ensure the persistence of such profitability in the long term, an adequate breeding program is required. To this aim, in the present study we estimate the genetic parameters of the main productive and reproductive traits, and we evaluate the effect of genotype by environment interaction (GxE) on these traits using 2 environmental covariates: (1) productivity and (2) temperature-humidity index (THI). Milk, fat, protein, and casein yield were considered as daily production traits, whereas protein, fat, casein percentage, casein index, and somatic cell score were considered as milk quality traits. Finally, reproductive traits such as the number of inseminations, days open, calving interval, and calving-to-first-insemination interval were evaluated. Reggiana cattle produce an average of 19 kg of milk per day with 3.7% fat and 3.4% protein content and have excellent fertility parameters. Compared with other breeds, they have slightly lower heritability for production and quality for production traits (e.g., 0.12 [0.09; 0.15] for milk yield), but similar heritability for fertility traits. Milk, protein, and fat daily yields are highly correlated but negatively correlated with the percentage of protein, fat, and casein, whereas fertility traits have an unfavorable genetic correlation with daily production traits. When considering productivity, a consistent amount of variability due to GxE was observed for all daily production traits, somatic cell count, and casein index. A modest amount of GxE was observed for fertility parameters, while the percentage of solid content showed almost no GxE effect. A similar situation occurred when considering the THI, but no GxE interaction was observed for reproduction traits. In conclusion, this study provides useful information for the implementation of accurate selection plans in this local breed, accounting for environmental plasticity measured through the consistent GxE interaction observed.

Genotype by Environment Interaction and Selection Response for Milk Yield Traits and Conformation in a Local Cattle Breed Using a Reaction Norm Approach

Local breeds are often reared in various environmental conditions (EC), suggesting that genotype by environment interaction (GxE) could influence genetic progress. This study aimed at investigating GxE and response to selection (R) in Rendena cattle under diverse EC. Traits included milk, fat, and protein yields, fat and protein percentage, and somatic cell score, three-factor scores and 24 linear type traits. The traits belonged to 11,085 cows (615 sires). Variance components were estimated in a two-step reaction norm model (RNM). A single trait animal model was run to obtain the solutions of herd-EC effect, then included in a random regression sire model. A multivariate response to selection (R) in different EC was computed for traits under selection including beef traits from a performance test. GxE accounted on average for 10% of phenotypic variance, and an average rank correlation of over 0.97 was found between bull estimated breeding values (EBVs) by either including or not including GxE, with changing top ranks. For various traits, significantly greater genetic components and R were observed in plain farms, loose housing rearing system, feeding total mixed ration, and without summer pasture. Conversely, for beef traits, a greater R was found for mountain farms, loose housing, hay-based feeding and summer pasture.

Derivation and genome-wide association study of a principal component-based measure of heat tolerance in dairy cattle

Heat stress represents a key factor that negatively affects the productive and reproductive performance of farm animals. In the present work, a new measure of tolerance to heat stress for dairy cattle was developed using principal component analysis. Data were from 590,174 test-day records for milk yield, fat and protein percentages, and somatic cell score of 39,261 Italian Holstein cows. Test-day records adjusted for main systematic factors were grouped into 11 temperature-humidity index (THI) classes. Daughter trait deviations (DTD) were calculated for 1,540 bulls as means of the adjusted test-day records for each THI class. Principal component analysis was performed on the DTD for each bull. The first 2 principal components (PC) explained 42 to 51% of the total variance of the system across the 4 traits. The first PC, a measure of the level at which the curve is located, was interpreted as a measure of the level at which the DTD curve was located. The second PC, which shows the slope of increasing or decreases DTD curves, synthesized the behavior of the DTD pattern. Heritability of the 2 component scores was moderate to high for level across all traits (range = 0.23-0.82) and low to moderate for slope (range = 0.16-0.28). For each trait, phenotypic and genetic correlations between level and slope were equal to zero. A genome-wide association analysis was carried out on a subsample of 423 bulls genotyped with the Illumina 50K bovine bead chip (Illumina, San Diego, CA). Two single nucleotide polymorphisms were significantly associated with slope for milk yield, 4 with level for fat percentage, and 2 with level and slope of protein percentage, respectively. The gene discovery was carried out considering windows of 0.5 Mb surrounding the significant markers and highlighted some interesting candidate genes. Some of them have been already associated with the mechanism of heat tolerance as the heat shock transcription factor (HSF1) and the malonyl-CoA-acyl carrier protein transacylase (MCAT). The 2 PC were able to describe the overall level and the slope of response of milk production traits across increasing levels of THI index. Moreover, they exhibited genetic variability and were genetically uncorrelated. These features suggest their use as measures of thermotolerance in dairy cattle breeding schemes.

Genotype by environment (climate) interaction improves genomic prediction for production traits in US Holstein cattle

Genotype by environment interaction (G × E) in dairy cattle productive traits has been shown to exist, but current genetic evaluation methods do not take this component into account. As several environmental descriptors (e.g., climate, farming system) are known to vary within the United States, not accounting for the G × E could lead to reranking of bulls and loss in genetic gain. Using test-day records on milk yield, somatic cell score, fat, and protein percentage from all over the United States, we computed within herd-year-season daughter yield deviations for 1,087 Holstein bulls and regressed them on genetic and environmental information to estimate variance components and to assess prediction accuracy. Genomic information was obtained from a 50k SNP marker panel. Environmental effect inputs included herd (160 levels), geographical region (7 levels), geographical location (2 variables), climate information (7 variables), and management conditions of the herds (16 total variables divided in 4 subgroups). For each set of environmental descriptors, environmental, genomic, and G × E components were sequentially fitted. Variance components estimates confirmed the presence of G × E on milk yield, with its effect being larger than main genetic effect and the environmental effect for some models. Conversely, G × E was moderate for somatic cell score and small for milk composition. Genotype by environment interaction, when included, partially eroded the genomic effect (as compared with the models where G × E was not included), suggesting that the genomic variance could at least in part be attributed to G × E not appropriately accounted for. Model predictive ability was assessed using 3 cross-validation schemes (new bulls, incomplete progeny test, and new environmental conditions), and performance was compared with a reference model including only the main genomic effect. In each scenario, at least 1 of the models including G × E was able to perform better than the reference model, although it was not possible to find the overall best-performing model that included the same set of environmental descriptors. In general, the methodology used is promising in accounting for G × E in genomic predictions, but challenges exist in identifying a unique set of covariates capable of describing the entire variety of environments.

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.

Genotype by environment interaction and breeding for robustness in livestock

The increasing size of the human population is projected to result in an increase in meat consumption. However, at the same time, the dominant position of meat as the center of meals is on the decline. Modern objections to the consumption of meat include public concerns with animal welfare in livestock production systems. Animal breeding practices have become part of the debate since it became recognized that animals in a population that have been selected for high production efficiency are more at risk for behavioral, physiological and immunological problems. As a solution, animal breeding practices need to include selection for robustness traits, which can be implemented through the use of reaction norms analysis, or though the direct inclusion of robustness traits in the breeding objective and in the selection index. This review gives an overview of genotype × environment interactions (the influence of the environment, reaction norms, phenotypic plasticity, canalization, and genetic homeostasis), reaction norms analysis in livestock production, options for selection for increased levels of production and against environmental sensitivity, and direct inclusion of robustness traits in the selection index. Ethical considerations of breeding for improved animal welfare are discussed. The discussion on animal breeding practices has been initiated and is very alive today. This positive trend is part of the sustainable food production movement that aims at feeding 9.15 billion people not just in the near future but also beyond.

Linear reaction norm models for genetic merit prediction of Angus cattle under genotype by environment interaction

The objectives of this work were to assess alternative linear reaction norm (RN) models for genetic evaluation of Angus cattle in Brazil. That is, we investigated the interaction between genotypes and continuous descriptors of the environmental variation to examine evidence of genotype by environment interaction (G×E) in post-weaning BW gain (PWG) and to compare the environmental sensitivity of national and imported Angus sires. Data were collected by the Brazilian Angus Improvement Program from 1974 to 2005 and consisted of 63,098 records and a pedigree file with 95,896 animals. Six models were implemented using Bayesian inference and compared using the Deviance Information Criterion (DIC). The simplest model was M(1), a traditional animal model, which showed the largest DIC and hence the poorest fit when compared with the 4 alternative RN specifications accounting for G×E. In M(2), a 2-step procedure was implemented using the contemporary group posterior means of M(1) as the environmental gradient, ranging from -92.6 to +265.5 kg. Moreover, the benefits of jointly estimating all parameters in a 1-step approach were demonstrated by M(3). Additionally, we extended M(3) to allow for residual heteroskedasticity using an exponential function (M(4)) and the best fitting (smallest DIC) environmental classification model (M(5)) specification. Finally, M(6) added just heteroskedastic residual variance to M(1). Heritabilities were less at harsh environments and increased with the improvement of production conditions for all RN models. Rank correlations among genetic merit predictions obtained by M(1) and by the best fitting RN models M(3) (homoskedastic) and M(5) (heteroskedastic) at different environmental levels ranged from 0.79 and 0.81, suggesting biological importance of G×E in Brazilian Angus PWG. These results suggest that selection progress could be optimized by adopting environment-specific genetic merit predictions. The PWG environmental sensitivity of imported North American origin bulls (0.046 ± 0.009) was significantly larger (P < 0.05) than that of local sires (0.012 ± 0.013). Moreover, PWG of progeny of imported sires exceeded that of native sires in medium and superior production levels. On the other hand, Angus cattle locally selected in Brazil tended to be more robust to environmental changes and hence be more suitable when production environments for potential progeny is uncertain.

Genotype × environment interaction for long-yearling weight in Canchim cattle quantified by reaction norm analysis

The objective of this study was to investigate the presence of genotype × environment interactions (G×E) for long-yearling weight in Canchim cattle (5/8 Charolais + 3/8 zebu) in Brazil using reaction norms (RN). The hierarchical RN model included the fixed effect of age of the animal (linear coefficient) and random effects of contemporary groups and additive animal genetic intercept and slope of the RN and contemporary group effects as random effects. Contemporary groups as the most elemental representation of management conditions in beef cattle were chosen to represent the environmental covariate of the RN. The deviance information criteria demonstrated that a homoskedastic residual RN model provided a better data fit compared with a heteroskedastic counterpart and with a traditional animal model, which had the worst fit. The environmental gradient for long-yearling weight based on contemporary group effects ranged from -105 to 150 kg. The additive direct variance and heritability estimates increased with increasing environmental gradient from 74.33 ± 22.32 to 1,922.59 ± 258.99 kg(2) and from 0.08 ± 0.02 to 0.68 ± 0.03, respectively. The high genetic correlation (0.90 ± 0.03) between the intercept and the slope of the RN shows that animals with the greatest breeding values best responded to environmental improvement, characterizing scale effect as the source of G×E for long-yearling weight. The phenotypic plasticity demonstrated by the slope of the RN of the animal indicates the possibility to change genotype expression along the environmental gradient through selection. The results demonstrate the importance of accounting for G×E in the genetic evaluation of this population.

Efficiency of alternative MCMC strategies illustrated using the reaction norm model

The Markov chain Monte Carlo (MCMC) strategy provides remarkable flexibility for fitting complex hierarchical models. However, when parameters are highly correlated in their posterior distributions and their number is large, a particular MCMC algorithm may perform poorly and the resulting inferences may be affected. The objective of this study was to compare the efficiency (in terms of the asymptotic variance of features of posterior distributions of chosen parameters, and in terms of computing cost) of six MCMC strategies to sample parameters using simulated data generated with a reaction norm model with unknown covariates as an example. The six strategies are single-site Gibbs updates (SG), single-site Gibbs sampler for updating transformed (a priori independent) additive genetic values (TSG), pairwise Gibbs updates (PG), blocked (all location parameters are updated jointly) Gibbs updates (BG), Langevin-Hastings (LH) proposals, and finally Langevin-Hastings proposals for updating transformed additive genetic values (TLH). The ranking of the methods in terms of asymptotic variance is affected by the degree of the correlation structure of the data and by the true values of the parameters, and no method comes out as an overall winner across all parameters. TSG and BG show very good performance in terms of asymptotic variance especially when the posterior correlation between genetic effects is high. In terms of computing cost, TSG performs best except for dispersion parameters in the low correlation scenario where SG was the best strategy. The two LH proposals could not compete with any of the Gibbs sampling algorithms. In this study it was not possible to find an MCMC strategy that performs optimally across the range of target distributions and across all possible values of parameters. However, when the posterior correlation between parameters is high, TSG, BG and even PG show better mixing than SG.