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 (rg) for MY, FY, FC, PC and CI were found to decrease as differences between environmental categories increased. The rg 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 rg 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.

Linear and nonlinear genetic relationships between type traits and productive life in US dairy goats.

Linear or nonlinear genetic relationships between productive life and functional productive life at 72 mo, with final score (SCO), stature, strength, dairyness (DAI), teat diameter, rear legs (side view), rump angle, rump width (RUW), fore udder attachment (FUA), rear udder height, rear udder arch, udder depth (UDD), suspensory ligament (SUS), and teat placement, as well as heritabilities and correlations were estimated from multibreed US dairy goat records. Productive life was defined as the total days in production until 72 mo of age (PL72) for goats having the opportunity to express the trait. Functional productive life (FPL72) was analyzed by incorporating first lactation milk yield, fat yield, protein yield, and SCO in the statistical model. Heritabilities and correlations were estimated using linear mixed models with pedigree additive genetic relationships and ASReml software. Nonlinearity of genetic relationships was assessed based on second-degree polynomial (quadratic) regression models, with the breeding values of PL72 or FPL72 as responses and the breeding values for each type trait (linear and quadratic) as predictor variables. Heritability estimates were 0.19, 0.14, 0.18, 0.20, 0.14, 0.07, 0.28, 0.20, 0.15, 0.13, 0.25, 0.18, 0.20, 0.21, 0.21, and 0.32 for PL72, FPL72, SCO, stature, strength, DAI, teat diameter, rear legs, rump angle, RUW, FUA, rear udder height, rear udder arch, UDD, SUS, and teat placement, respectively. The type traits SCO, RUW, and FUA were the most correlated with PL72 and FPL72, so these may be used as selection criteria to increase longevity in dairy goats. An increase in the coefficient of determination >1% for the second degree, compared with that for the linear model for either PL72 or FPL72, was taken as evidence of a nonlinear genetic relationship. Using this criterion, PL72 showed maximum values at intermediate scores in DAI, UDD, and RUW, and maximum values at extreme scores in FUA and SUS, whereas FPL72 showed maximum values at intermediate scores in DAI and UDD, and maximum values at extreme scores in FUA, RUW, and SUS. Selecting for increased SCO, RUW, and FUA will lead to an increase of FPL72 in goats. Consideration of nonlinear relationships between DAI, FUA, RUW, SUS, and UDD may help in the design of more efficient breeding programs for dairy goats using conformation traits.

Effect of reference population size and available ancestor genotypes on imputation of Mexican Holstein genotypes.

The effects of reference population size and the availability of information from genotyped ancestors on the accuracy of imputation of single nucleotide polymorphisms (SNP) were investigated for Mexican Holstein cattle. Three scenarios for reference population size were examined: (1) a local population of 2,011 genotyped Mexican Holsteins, (2) animals in scenario 1 plus 866 Holsteins in the US genotype database (GDB) with genotyped Mexican daughters, and (3) animals in scenario 1 and all US GDB Holsteins (338,073). Genotypes from 4 chip densities (2 low density, 1 mid density, and 1 high density) were imputed using findhap (version 3) to the 45,195 markers on the mid-density chip. Imputation success was determined by comparing the numbers of SNP with 1 or 2 alleles missing and the numbers of differently predicted SNP (conflicts) among the 3 scenarios. Imputation accuracy improved as chip density and numbers of genotyped ancestors increased, and the percentage of SNP with 1 missing allele was greater than that for 2 missing alleles for all scenarios. The largest numbers of conflicts were found between scenarios 1 and 3. The inclusion of information from direct ancestors (dam or sire) with US GDB genotypes in the imputation of Mexican Holstein genotypes increased imputation accuracy by 1 percentage point for low-density genotypes and by 0.5 percentage points for high-density genotypes, which was about half the gain found with information from all US GDB Holsteins. A larger reference population and the availability of genotyped ancestors improved imputation; animals with genotyped parents in a large reference population had higher imputation accuracy than those with no or few genotyped relatives in a small reference population. For small local populations, including genotypes from other related populations can aid in improving imputation accuracy.

Estimation of genetic parameters for productive life, reproduction, and milk-production traits in US dairy goats.

Heritabilities and correlations for milk yield (MY), fat yield (FY), protein yield (PY), combined fat and protein yield (FPY), fat percentage (F%), protein percentage (P%), age at first kidding (AFK), interval between the first and second kidding (KI), and real and functional productive life at 72mo (FPL72) of 33,725 US dairy goats, were estimated using animal models. Productive life was defined as the total days in production until 72mo of age (PL72) for goats having the opportunity to express the trait. Functional productive life was obtained by correcting PL72 for MY, FY, PY, and final type score (FS). Six selection indexes were used, including or excluding PL72, with 6 groups of different economic weights, to estimate the responses to selection considering MY, FY, PY, and PL72 as selection criteria. The main criteria that determined the culling of a goat from the herd were low FS, MY, and FY per lactation. Heritability estimates were 0.22, 0.17, 0.37, 0.37, 0.38, 0.39, 0.54, 0.64, 0.09, and 0.16 for PL72, FPL72, MY, FY, PY, FPY, F%, P%, KI, and AFK, respectively. Most genetic correlations between the evaluated traits and PL72 or FPL72 were positive, except for F% (-0.04 and -0.06, respectively), P% (-0.002 and -0.03, respectively), and AFK (-0.03 and -0.01, respectively). The highest genetic correlations were between FPL72 and MY (0.39) and between PL72 and MY (0.33). Most phenotypic correlations between the traits evaluated and FPL72 and PL72 were positive (>0.23 and >0.26, respectively), except for F% (-0.004 and -0.02, respectively), P% (-0.05 and -0.02), KI (-0.01 and -0.07), and AFK (-0.08 and -0.08). The direct selection for PL72 increased it by 102.28d per generation. The use of MY, FY, PY, KI, or AFK as selection criteria increased PL72 by 39.21, 27.33, 35.90, -8.28, or 2.77d per generation, respectively. The inclusion of PL72 as selection criterion increased the expected response per generation from 0.15 to 17.35% in all selection indices studied.

Evaluation of heat stress effects on production traits and somatic cell score of Holsteins in a temperate environment.

This study was aimed to evaluate the degree of thermal stress exhibited by Holsteins under a continental temperate climate. Milk, fat, protein, and somatic cell count test-day records collected between 2000 and 2011 from 23,963 cows in 604 herds were combined with meteorological data from 14 public weather stations in Luxembourg. Daily values of 6 different thermal indices (TI) weighted in term of temperature, relative humidity, solar radiation, and wind speed were calculated by averaging hourly TI over 24h. Heat stress thresholds were first identified by a broken-line regression model. Regression models were thereafter applied to quantify milk production losses due to heat stress. The tipping points at which milk and protein yields declined were effectively identified. For fat yield, no valid threshold was identified for any of the studied TI. Daily fat yields tended to decrease steadily with increasing values of TI. Daily somatic cell score patterns were marked by increased values at both lowest and highest TI ranges, with a more pronounced reaction to cold stress for apparent temperature indices. Thresholds differed between TI and traits. For production traits, they ranged from 62 (TI(1)) to 80 (TI(3)) for temperature-humidity indices (THI) and from 16 (TI(5)) to 20 (TI(6)) for apparent temperature indices. Corresponding somatic cell score thresholds were higher and ranged from 66 (TI(1)) to 82 (TI(3)) and from 20 (TI(5)) to 23 (TI(6)), respectively. The largest milk decline per unit of mild, moderate, and extreme heat stress levels of 0.164, 0.356, and 0.955 kg, respectively, was observed when using the conventional THI (TI(1)). The highest yearly milk, fat, and protein losses of 54, 5.7, and 4.2 kg, respectively, were detected by TI(2), the THI index that is adjusted for wind speed and solar radiation. The latter index could be considered as the best indicator of heat stress to be used for forecast and herd management in a first step in temperate regions under anticipated climate changes.

Breed and breed x environment interaction effects for growth traits and survival rate from birth to weaning in crossbred lambs.

This study was conducted to estimate the effects of sire breed (Charollais, Dorper, Dorset, Hampshire, and Suffolk) and genetic group of the dam [local whiteface (Criolla), grade Hampshire, grade Suffolk, F(1) Dorper × Pelibuey, Australian composite, hair sheep (Barbados Blackbelly or Pelibuey), Dorset, Hampshire and Suffolk] on birth weight, survival rate from birth to weaning, and weaning weight adjusted to 65 d of age. Data were recorded from 2004 to 2006 on crossbreed lambs obtained by artificial insemination from 114 flocks in Central Mexico. High (above the mean) or low (below the mean) environmental categories were assigned to each record from flock-year-season effects solutions obtained in fixed-effects linear model analyses. Birth weights of lambs sired by Charollais rams (3.94 kg) were heavier (P < 0.05) than those sired by Suffolk (3.69 kg) whereas Dorper-, Dorset-, and Hampshire-sired lambs were intermediate. For weaning weight, lambs sired by Charollais were the heaviest (19.16 kg) compared with lambs sired by Hampshire (17.86 kg), Suffolk (17.79 kg), and Dorper (17.28 kg) whereas Dorset sired were the lightest (16.77 kg; P < 0.05). Lambs sired by Dorset rams had a lower survival rate (77.1%) than lambs sired by Charollais, Dorper, Hampshire, or Suffolk (81.6 to 83.3%; P < 0.05). Lambs from Australian composite, Dorset, Hampshire, and Suffolk dams were heavier for birth and weaning compared with lambs from local whiteface (Criolla), hair sheep breeds (Barbados Blackbelly or Pelibuey), or crosses (F(1) Pelibuey × Dorper) dams (P < 0.05) whereas lambs from grade Hampshire and grade Suffolk were intermediate. The genetic group of the dam had no effect (P > 0.05) on survival rate. Sire breed × environmental category interaction effect was significant for birth weight (P < 0.01) and weaning weight (P < 0.05), but relatively small changes on the ranking of sire breeds were observed between environmental categories for weaning weight. Genetic group of the dam × environmental category interaction effect was significant for birth and weaning weights (P < 0.01). Weaning weight of lambs from hair sheep breeds (Barbados Blackbelly or Pelibuey) and crosses (F(1) Pelibuey × Dorper) as well as purebred Dorset, Hampshire, and Suffolk dams were more affected when changing from the high to the low environmental category compared with the other genetic groups. No breed of the sire or genetic group of the dam × environmental category interactions were (P > 0.05) observed for survival rate.

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.

Genetic and environmental parameters for milk production, udder health, and fertility traits in Mexican Holstein cows.

Genetic and phenotypic parameters for Mexican Holstein cows were estimated for first- to third-parity cows with records from 1998 to 2003 (n=2,971-15,927) for 305-d mature equivalent milk production (MEM), fat production (MEF), and protein production (MEP), somatic cell score (SCS), subsequent calving interval (CAI), and age at first calving (AFC). Genetic parameters were obtained by average information matrix-REML methodology using 6-trait (first-parity data) and 5-trait (second- and third-parity data) animal models. Heritability estimates for production traits were between 0.17+/-0.02 and 0.23+/-0.02 for first- and second-parity cows and between 0.12+/-0.03 and 0.13+/-0.03 for third-parity cows. Heritability estimates for SCS were similar for all parities (0.10+/-0.02 to 0.11+/-0.03). For CAI, estimates of heritability were 0.01+/-0.05 for third-parity cows and 0.02+/-0.02 for second-parity cows. The heritability for AFC was moderate (0.28+/-0.03). No unfavorable estimates of correlations were found among MEM, MEF, MEP, CAI, and SCS. Estimates of environmental and phenotypic correlations were large and positive among production traits; favorable between SCS and CAI; slightly favorable between MEM, MEF, and MEP and SCS, between AFC and SCS, and between SCS and CAI; and small but unfavorable between production traits and CAI. Estimates of genetic variation and heritability indicate that selection would result in genetic improvement of production traits, AFC, and SCS. Estimates of both heritability and genetic variation for CAI were small, which indicates that genetic improvement would be difficult.

Short communication: genetic and environmental relationships between milk yield and kidding interval in dairy goats.

First-parity 305-d milk yield and intervals between first and second kiddings from 1975 through 2005 were analyzed to estimate genetic and environmental parameters for United States Alpine, LaMancha, Nubian, Saanen, and Toggenburg dairy goats. The data set included 43,612 does with first-parity milk yield and 25,863 does with first kidding interval. Parameters were estimated by REML using bivariate mixed models within and across breeds. Fixed effects were month and herd-year of kidding within breed and month and breed-herd-year of kidding across breeds. Random effects were animal and residual. The heritability estimate for first-parity milk yield was 0.36+/-0.01 across breeds and ranged from 0.35 to 0.38 within breed; heritability for first kidding interval was 0.05+/-0.01 across breeds and ranged from 0.00 to 0.15 within breed. The estimate of the genetic correlation between first-parity milk yield and first kidding interval was positive (unfavorable) across breeds (0.35+/-0.09). Estimates of environmental correlations within and across breeds were positive (0.16-0.25). The presence of unfavorable genetic relationships between milk yield and kidding interval indicates a need to include reproductive performance as a selection criterion.

Selection response for milk production in conventional production systems in Mexico, using genetic evaluations of Holstein sires from Canada and the United States.

Polynomial regression models of the first, second, and third order were used to fit milk production deviations of daughters in Mexico on Canadian and US predicted transmitting ability values for 305-d mature-equivalent milk production (kg). For the pairs Canada-Mexico and Mexico-United States, 40 and 73 bulls with a minimum reliability of 0.75 were analyzed, respectively. Genetic correlations between pairs of countries were also estimated. The parameters were evaluated for all data, and for sires grouped according to the mean of the average phenotypic milk production (high and low) of their daughters' herd mates. Quadratic and cubic effects were not significant in any analysis. From linear regression models, slopes of Mexican daughter deviations on US and Canadian predicted transmitting abilities were 1.01 and 0.93, respectively. Slopes were greater but intercepts were smaller for the high versus low level of production of the sires' herd mates in Mexico. A greater difference between the genetic correlations was found for the high versus low environmental level than for the low level (0.79 vs. 0.70) for Mexico-US data compared with Canada-Mexico data (0.81 vs. 0.78). Genetic correlations between Mexico and the United States (0.74), and between Mexico and Canada (0.77), were smaller than the genetic correlation between the same Canadian and US sires (0.92), suggesting the presence of a moderate degree of genotype-environment interaction for milk production between Canada and the United States, and Mexico.

Short communication: Changes in the association between milk yield and age at first calving in holstein cows with herd environment level for milk yield.

The aim of this study was to evaluate the effect of herd environment class on the genetic and phenotypic relationships of mature equivalent milk yield (MY) on age at first calving (AFC). Data analyzed were 248,230 first parity records of Holstein cows, daughters of 588 sires in 3,042 herds in the United States. Heritability for AFC was 0.33 +/- 0.01 and 0.20 +/- 0.01 in high and low environment herds, respectively, and 0.47 +/- 0.01 in the complete data set. The correlation between AFC sires' predicted breeding values of low and high classes was 0.69. Genetic correlations between MY and AFC were -0.52 +/- 0.02 and -0.31 +/- 0.03 in high and low environment herds, respectively, and -0.44 +/- 0.02 in the complete data set representing intermediate environments. If selection is based on the whole data set, expected correlated responses for AFC estimated as a result of 1,000 kg of genetic gain in MY, for high and low herd environment herds were -26.1 and -15.3 d, respectively, and -32.6 for the complete data set; hence the highest reduction in AFC occurs in intermediate environment herds. Different estimates of the heritability of AFC, the correlation between AFC breeding values of low and high classes as well as changes in the genetic correlation between MY and AFC across environments indicate genotype x environment interaction. Caution in interpretation is warranted because genetic relationships are dynamic, especially in populations undergoing selection. Current relationships may differ from those during the time period of the present study (1987-1994). Notwithstanding this possibility, methods and findings from the present study provide insight about the complexity of genetic association and genotype x environment interactions between AFC and MY.