Genotype by environment interaction in organic and conventional production systems and their consequences for breeding objectives in Austrian Fleckvieh cattle

Heterogeneity of variance and genetic parameters for milk production in cattle, using Bayesian inference

The goal of this study was to verify the effect of heterogeneity of variance (HV) on milk production in up to 305 days of lactation (L305) of daughters of Girolando, Gir and Holstein sires, as well as in the genetic evaluation of these sires and their progenies. in Brazil. The model included contemporary groups (consisting of herd, year and calving season) as a fixed effect, cow age at calving (linear and quadratic effects) and heterozygosity (linear effect) as covariates, in addition to the random effects of direct additive genetic and environmental, permanent and residual. The first analysis consisted of the single-trait animal model, with L305 records (disregarding HV). The second considered classes of standard deviations (SD): two-trait model including low and high classes (considering HV), according to the standardized means of L305 for herd-year of calving. The low SD class was composed of herds with SD equal to or less than zero and the high class with positive SD values. Estimates of (co)variance components and breeding values were obtained separately for each scenario using Bayesian inference via Gibbs sampling. Different heritability was estimated. Higher for the high DP class in the Gir (0.20) and Holstein (0.15) breeds, not occurring the same in the Girolando breed, with a lower value among the classes for the high DP (0.10). High values of genetic correlations were also found between low and high SD classes (0.88; 0.85 and 0.79) for the Girolando, Gir and Holstein breeds, respectively. Like the order correlations (Spearman) which were also high for the three breeds analyzed (equal to or above 0.92). Thus, the presence of HV had a smaller impact for L305 and did not affect the genetic evaluation of sires.

Assessing various environmental descriptors with respect to genotype x environment interaction for milk production traits in Bavarian Fleckvieh cattle

Genotype x environment (GxE) interaction for production traits in Fleckvieh cattle was assessed by means of various environmental descriptors (EDs). It was also of particular interest to search for EDs useful for studying differing robustness or resilience of individuals which implies reasonable GxE interaction. The set of studied EDs included farm/herd environment (e.g. herd size, housing/feeding regimes, herd production level), geographical location (e.g. height above sea level), temperature humidity index and fat-to-protein ratio. Milk, fat and protein yield deviations for the first test day, the first 100 days and the 305-day-yield of the first lactation were used as phenotypes. Genetic correlations were estimated with bi- and multivariate sire models in case of categorical EDs and reaction norm sire models in case of continuous EDs. Further, rank correlations of bulls' estimated breeding values were calculated for different environmental levels/classes. Almost all estimated genetic correlations were significantly different from 1, but not <0.93. There were some exceptions for extremely different classes/levels of some EDs (e.g. average herd-year production, organic vs conventional farm systems, way of concentrated feed supply). Rank correlations were rarely below 0.95. In summary, no substantial GxE interactions for milk production traits were found with the studied EDs.

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.

Proofs for genotype by environment interactions considering pedigree and genomic data from organic and conventional cow reference populations

The aim of the present study was to prove genotype by environment interactions (G × E) for production, longevity, and health traits considering conventional and organic German Holstein dairy cattle subpopulations. The full data set included 141,778 Holstein cows from 57 conventional herds and 7,915 cows from 9 organic herds. The analyzed traits were first-lactation milk yield and fat percentage (FP), the length of productive life (LPL) and the health traits mastitis, ovarian cycle disorders, and digital dermatitis in first lactation. A subset of phenotyped cows was genotyped and used for the implementation of separate cow reference populations. After SNP quality controls, the cow reference sets considered 40,830 SNP from 19,700 conventional cows and the same 40,830 SNP from 1,282 organic cows. The proof of possible G × E was made via multiple-trait model applications, considering same traits from the conventional and organic population as different traits. In this regard, pedigree (A), genomic (G) and combined relationship (H) matrices were constructed. For the production traits, heritabilities were very similar in both organic and conventional populations (i.e., close to 0.70 for FP and close to 0.40 for milk yield). For low heritability health traits and LPL, stronger heritability fluctuations were observed, especially for digital dermatitis with 0.05 ± 0.01 (organic, A matrix) to 0.33 ± 0.04 (conventional, G matrix). Quite large genetic correlations between same traits from the 2 environments were estimated for production traits, especially for high heritability FP. For LPL, the genetic correlation was 0.67 (A matrix) and 0.66 (H matrix). The genetic correlation between LPL organic with LPL conventional was 0.94 when considering the G matrix, but only 213 genotyped cows were included. For health traits, genetic correlations were throughout lower than 0.80, indicating possible G × E. Genetic correlations from the different matrices A, G, and H for health and production traits followed the same pattern, but the estimates from G for health traits were associated with quite large standard errors. In genome-wide association studies, significantly associated SNP for production traits overlapped in the conventional and organic population. In contrast, for low heritability LPL and health traits, significantly associated SNP and annotated potential candidate genes differed in both populations. In this regard, significantly associated SNP for mastitis from conventional cows were located on Bos taurus autosomes 6 and 19, but on Bos taurus autosomes 1, 10, and 22 in the organic population. For the remaining health traits and LPL, different potential candidate genes were annotated, but the different genes reflect similar physiological pathways. We found evidence of G × E for low heritability functional traits, suggesting different breeding approaches in organic and conventional populations. Nevertheless, for a verification of results and implementation of alternative breeding strategies, it is imperative to increase the organic cow reference population.

Genotype-by-environment interactions at the trait level and total merit index level for milk production and functional traits in Brown Swiss cattle

The production environments of the German-Austrian Brown Swiss population show a wide range due to differences in topography, landscapes, local climates, and different farm management systems. Extensive production systems such as organic farming have become increasingly popular in recent decades because of interest in sustainability and consumer preferences. Compared with conventional farmers, organic farmers put more weight on fitness traits. Besides the official total merit index (TMI), a selection index applying relative economic weights (REWs) suitable for organic production systems is provided for Brown Swiss cattle in Germany. The aim of the study was to investigate genotype-by-environment interactions (GxE) for milk production traits and functional traits (including longevity, fertility traits, and calving traits) in a sample of the German-Austrian Brown Swiss population housed in Baden-Wuerttemberg (southern Germany) by applying bivariate and random regression sire models. For bivariate analyses, the production environment was binary classified by farm management system (organic and conventional) and altitude of farm location (above or below 800 m above sea level (ASL)). Milk energy yields (MEY) obtained from herd effects were used as continuously scaled environmental descriptor in the reaction norm approach. The TMIs for sires were calculated based on breeding values estimated with different models and environment-specific REWs to determine possible GxE at TMI levels and rerankings of sires. In bivariate analyses, genetic correlations at the trait level were high and ranged from r_g = 0.99 (calving to first insemination, cystic ovaries, and maternal stillbirth rate) to r_g = 0.79 (first insemination to conception for altitude). Except for the latter, no severe GxE were found at the trait level using the bivariate models. Fat yield was the only trait showing minor GxE in the reaction norm model approach. Investigating the environmental sensitivity at the TMI level revealed rank correlations between the different environment-specific TMIs that were close to unity, implying no severe reranking effects. The results show no need to account for different environments in Brown Swiss cattle breeding programs.

Genotype by environment interaction in response to cold stress in a composite beef cattle breed

Extreme weather conditions such as cold stress influence the productivity and survivability of beef cattle raised on pasture. The objective of this study was to identify and evaluate the extent of the impact of genotype by environment interaction due to cold stress on birth weight (BW) and weaning weight (WW) in a composite beef cattle population. The effect of cold stress was modelled as the accumulation of total cold load (TCL) calculated using the Comprehensive Climate Index units, considering three TCL classes defined based on temperature: less than -5°C (TCL5), -15°C (TCL15) and -25°C (TCL25). A total of 4221 and 4217 records for BW and WW, respectively, were used from a composite beef cattle population (50% Red Angus, 25% Charolais and 25% Tarentaise) between 2002 and 2015. For both BW and WW, a univariate model (ignoring cold stress) and a reaction norm model were implemented. As cold load increased, the direct heritability slightly increased in both BW and WW for TCL5 class; however, this heritability remained consistent across the cold load of TCL25 class. In contrast, the maternal heritability of BW was constant with cold load increase in all TCL classes, although a slight increase of maternal heritability was observed for TCL5 and TCL15. The direct and maternal genetic correlation for BW and maternal genetic correlation for WW across different cold loads between all TCL classes were high (r > 0.99), whereas the lowest direct genetic correlations observed for WW were 0.88 for TCL5 and 0.85 for TCL15. The Spearman rank correlation between the estimated breeding value of top bulls (n = 79) using univariate and reaction norm models across TCL classes showed some re-ranking in direct and maternal effects for both BW and WW particularly for TCL5 and TCL15. In general, cold stress did not have a big impact on direct and maternal genetic effects of BW and WW.

Genotype by environment interaction for yearling weight in Nellore cattle applying reaction norms models

In the present study, a possible existence of genotype × environment interaction was verified for yearling weight in Nellore cattle, utilising a reaction norms model. Therefore, possible changes in the breeding value were evaluated for 46 032 animals, from three distinct herds, according to the environmental gradient variation of the different contemporary groups. Under a Bayesian approach, analyses were carried out utilising INTERGEN software resulting in solutions of contemporary groups dispersed in the environmental gradient from –90 to +100 kg. The estimates of heritability coefficients ranged from 0.19 to 0.63 through the environmental gradient and the genetic correlation between intercept and slope of the reaction norms was 0.76. The genetic correlation considering all animals of the herds in the environmental gradient ranged from 0.83 to 1.0, and the correlation between breeding values of bulls in different environments ranged from 0.79 to 1.0. The results showed no effect of genotype × environment interaction on yearling weight in the herds of this study. However, it is important to verify a possible influence of the genotype × environment in the genetic evaluation of beef cattle, as different environments might cause interference in gene expression and consequently difference in phenotypic response.