Short Communication: Changes in the Association Between Milk Yield and Age at First Calving in Holstein Cows with Herd Environment Level for Milk Yield

Genotype by climate zone interactions for fertility, somatic cell score, and production in Iranian Holsteins

The objective of this study was to investigate genetic variation and genotype by environment (G × E) interactions for fertility (including age at first calving and calving interval), somatic cell score (SCS), and milk production traits for Iranian Holsteins. Different environments were defined based on the climatic zones (cold, semi-cold, and moderate) and considering the herd location. Data were collected between 2003 and 2018 by the National Animal Breeding Center of Iran (Karaj). Variance and covariance components and genetic correlations were estimated using 2 different models, which were analyzed using Bayesian methods. For both models, performance of traits in each climate were considered as different traits. Fertility traits were analyzed using a trivariate model. Furthermore, SCS and production traits were analyzed using trivariate random regression models (records in different climate zones considered as different traits). For the fertility traits, the largest estimates of heritability were observed in cold climate. Fertility performance was always better in cold environment. Genetic correlations between climatic zones ranged from 0.85 to 0.94. For daily measurements of SCS and production traits, heritability ranged from 0.031 to 0.037 and 0.069 to 0.209, respectively. Genetic variances were the highest in the semi-cold and moderate climates for the SCS and production traits, respectively. Furthermore, across the studied climates, 305-d genetic correlation ranged from 0.756 to 0.884 for SCS and from 0.925 to 0.957 for the production traits. The structure of genetic correlation within each climate indicated a negative correlation between early and late lactation for SCS, especially in the cold climate and for milk production in the moderate climate. For fat percentage, in all climatic zones, the lowest genetic correlations were observed between early and mid-lactation. In addition, for protein production in the cold climate, a negative correlation was observed between early and late lactation. Results indicated heterogeneous variance components for all the studied traits across various climatic zones. Estimated genetic correlations for SCS revealed that the genetic expression of animals may vary by climatic zone. Results indicated the existence of G × E interaction due to the climatic condition, only for SCS. Therefore, in Iranian Holsteins, the effect of G × E interactions should not be neglected, especially for SCS, as different sires might be optimal for use in different climatic zones.

Genotype × region and genotype × production level interactions in Holstein cows

The identification of the presence of genotype by environment interaction effects on important traits in Holstein cattle allows for the use of international genetic evaluations and a more efficient design of regional genetic evaluation programmes. The aim of this study was to determine the genotype × environment interaction effects in Chilean Holstein dairy cattle through the analysis of records corresponding to calvings between 1998 and 2015. Herds were classified in the central and southern regions of Chile based on herd location as well as by high and low levels of production environments based on the fat plus protein yield averages per herd within each region. The central region has a Mediterranean climate and a confined production system while the southern region has a humid temperate climate and a production system based on grazing with supplementation. Traits studied were milk yield (MY), fat yield (FY), protein yield (PY), fat content (FC) and protein content (PC) by lactation, age at first calving (AFC) and calving interval (CI). Several four-trait mixed animal models were applied to environmental category data as different traits, which included herd-year-calving season (herd-year-birth season for AFC) and lactation number as fixed effects, and animal additive genetic, sire-herd, permanent environment and residual effects as random effects. Genetic correlations (r_g) for MY, FY, FC, PC and CI were found to decrease as differences between environmental categories increased. The r_g between the most extreme environmental categories considered in this study for AFC (0.26) was the only one found statistically lower than 0.60. Genetic correlation values statistically lower than 0.80 (P < 0.05) were observed for AFC, CI, MY, FY and PY between some environmental categories. If separate genetic evaluations are adopted as practical criteria when the value of r_g is lower than 0.60, the consequence of improving a multi-trait economic breeding objective in this population is likely to be small unless extreme environmental categories are considered. However, a moderate decrease in selection response and re-ranking of selection candidates is expected for AFC, CI and yield traits when selection is performed in different environmental conditions. Genotype × environment interaction effects involving production systems in a Mediterranean climate and confinement vs. Temperate Oceanic climate and grazing with supplementation, and between two fat plus protein yield level categories within each environment, were at most moderate for the studied traits.

Genomic Tools for the Conservation and Genetic Improvement of a Highly Fragmented Breed-The Ramo Grande Cattle from the Azores

Simple Summary

Inbreeding control is a key concern in managing local endangered breeds, which often have developed unique adaptation features. Ramo Grande is a local cattle breed raised in the Azores archipelago under very harsh conditions, with a census of about 1300 cows dispersed by various islands. This fragmentation is a challenge when the goal is to keep inbreeding under control. Currently, panels of genetic markers are available which enable the assessment of inbreeding and the occurrence of previous bottlenecks in a population. These panels also allow the identification of genes associated with specific production traits, if reliable phenotypic information is available. We used a panel of genetic markers and estimated that the degree of inbreeding was approaching a level of concern, while some exotic gene inflow may have occurred in the past. We were able to identify genetic markers significantly associated with longevity, which reflects the ability of these cattle to remain productive under severe environmental conditions. Genetic markers were also identified as significantly associated with age at first calving and calf growth rate. The results indicate that genomic information can be used to control inbreeding and to implement genomic selection in Ramo Grande cattle to enhance adaptation and production traits.

Abstract

Ramo Grande is a local cattle breed raised in the archipelago of Azores, with a small and dispersed census, where inbreeding control is of utmost importance. A single nucleotide polymorphism (SNP) Beadchip array was used to assess inbreeding, by analysis of genomic regions harboring contiguous homozygous genotypes named runs of homozygosity (ROH), and to estimate past effective population size by analysis of linkage disequilibrium (LD). Genetic markers associated with production traits were also investigated, exploiting the unique genetic and adaptation features of this breed. A total of 639 ROH with length >4 Mb were identified, with mean length of 14.96 Mb. The mean genomic inbreeding was 0.09, and long segments of ROH were common, indicating recent inbred matings. The LD pattern indicates a large effective population size, suggesting the inflow of exotic germplasm in the past. The genome-wide association study identified novel markers significantly affecting longevity, age at first calving and direct genetic effects on calf weight. These results provide the first evidence of the association of longevity with genes related with DNA recognition and repair, and the association of age at first calving with aquaporin proteins, which are known to have a crucial role in reproduction.

Somatic cell counts, reproduction indicators, and technical efficiency in milk production: A stochastic frontier analysis for Spanish dairy farms

Health and reproductive conditions affect milk yields in dairy cows and may lead to differences in technical efficiency across farms. To investigate this, we created a novel panel data set of 197 dairy farms in northern Spain observed over the period 2006-2014 by combining information from 2 different sources and including data on production variables, somatic cell count, and reproduction indicators, as well as a genetic index. We used these data to estimate a stochastic production frontier where the somatic cell count, age at first calving, and calving interval are included as determinants of technical efficiency. Higher somatic cell count, greater age at first calving, and longer calving intervals were found to decrease technical efficiency, translating into significant losses in profits which we quantify through a simulation exercise. We also quantified the effect of genetic selection on profits, taking into account not only its direct effect on productivity but also its indirect effect through reduced technical efficiency due to the effect of genetic selection on the health and reproduction variables.

Phenotypic and genetic relationships between age at first calving, its component traits, and survival of heifers up to second calving

The aim of this study was to answer the question whether models for genetic evaluations of longevity should include a correction for age at first calving (AFC). For this purpose, phenotypic and genetic relationships between AFC, its component traits age at first insemination (AFI) and interval from first to last insemination (FLI), and survival of different periods of the first lactation (S1: 0 to 49 d, S2: 50 to 249 d, S3: 250 d to second calving) were investigated. Data of 721,919 German Holstein heifers, being inseminated for the first time during the years from 2003 to 2012, were used for the analyses. Phenotypic correlations of AFI, FLI, and AFC to S1 to S3 were negative. Mean estimated heritabilities were 0.239 (AFI), 0.007 (FLI), and 0.103 (AFC) and 0.023 (S1), 0.016 (S2), and 0.028 (S3) on the observed scale. The genetic correlation between AFI and FLI was close to zero. Genetic correlations between AFI and the survival traits were -0.08 (S1), -0.02 (S2), and -0.10 (S3); those between FLI and the survival traits were -0.14 (S1), -0.20 (S2), and -0.44 (S3); and those between AFC and the survival traits were -0.09 (S1), -0.06 (S2), and -0.20 (S3). Some of these genetic correlations were different from zero, which suggests that correcting for AFC in genetic evaluations for longevity in dairy cows might remove functional genetic variance and should be reconsidered.

Breed differences over time and heritability estimates for production and reproduction traits of dairy goats in the United States

To aid in improvement of breeding programs for production and reproduction traits of US dairy goats, breed differences over time were documented and genetic parameters were estimated. Data were from herds with ≥2 breeds (Alpine, LaMancha, Nubian, Oberhasli, Saanen, or Toggenburg), but only purebred data were analyzed. Three kidding periods were examined: 1976 through 1984, 1985 through 1994, and 1995 through 2005. Univariate repeatability mixed models were used to estimate least squares means by kidding period-breed and genetic parameters for milk, fat, and protein yields, combined fat and protein yield, fat and protein percentages, protein:fat ratio, age at first kidding, and kidding interval. Trends across kidding periods were favorable for most yield traits for all breeds but generally unfavorable for reproduction traits. Saanens had the highest milk (1,063 to 1,125 kg) and protein yields (31 to 33 kg). Nubians had the highest fat yields (37 to 40 kg) and lowest milk yields (791 to 851 kg). Oberhaslis had the lowest fat (31 to 33 kg) and protein (23 to 27 kg) yields. Alpines had the largest increase in milk yield (7.4%); Oberhaslis had the largest increase in protein (17.4%) and combined fat and protein (13.2%) yields. Combined fat and protein yield was higher for Nubians, Saanens, and Alpines (65 to 72 kg) than for LaManchas, Toggenburgs, and Oberhaslis (53 to 67 kg). Nubians had the highest fat (4.7 to 4.8%) and protein (3.6 to 3.8%) percentages. Only Nubians increased in fat percentage (2.1%); protein percentage increased most for Toggenburgs (7.4%) and Alpines (7.1%). Protein:fat ratio was highest for Toggenburgs (0.84 to 0.89) and lowest for Nubians (0.76 to 0.81), but Nubians had the largest increase in protein:fat ratio (6.6%). Saanens were oldest at first kidding (509 to 589 d), and Toggenburgs and LaManchas generally were youngest (435 to 545 d); age at first kidding increased most for Alpines (21.8%) and LaManchas (21.6%). Kidding intervals generally were shorter for Oberhaslis, LaManchas, and Nubians (350 to 377 d) than for Toggenburgs, Alpines, and Saanens (373 to 387 d). Kidding interval increased most for Nubians (3.9%) and Saanens (3.8%) and decreased only for Oberhaslis (5.4%). Heritability estimates across breeds were 0.35 for milk and fat yields, 0.37 for protein yield and protein:fat ratio, 0.36 for combined fat and protein yield, 0.52 for fat percentage, 0.54 for protein percentage, 0.23 for age at first kidding, and 0.05 for kidding interval. Genetic selection within breed is feasible for production and reproduction traits of US dairy goats.