Effects of cytoplasmic inheritance on production traits of dairy cattle

Pedigrees of 4461 cows were traced to the original female in a maternal line. Cytoplasmic origin was defined as the first female in the maternal lineage. There were 102 cytoplasmic lines. Most cows were at least 10 generations removed from the origin. After adjustment for sire, herd, calving year, calving month, and age, cytoplasmic effects accounted for 2.0, 1.8, 1.8, and 3.5% of total variation of milk yield, milk fat yield, 3.7% fat-corrected milk yield, and milk fat percentage in first lactation. Cytoplasmic effects were also in models that included adjustments for sires, maternal grand-sires, and dam's production. Correlations among independent subsets agreed with expectations. Cytoplasmic origin was a significant source of variation of production traits of dairy cattle.     

A further look at evidence for cytoplasmic inheritance of production traits in dairy cattle

A recent study interpreted significant effects of source of cytoplasm as evidence for cytoplasmic inheritance. An alternative interpretation that observed results reflected residual additive genetic effects unaccounted for by statistical models used was tested by computer simulation of records over 60 yr for milk yield (h2 = .30) and fat percentage (h2 = .60) under an additive genetic model with no cytoplasmic effects. First lactation records for the last 30 yr were analyzed by least squares according to models that included effects of cytoplasmic source, generations to cytoplasmic source, herd, and year as well as effects of sire (model 1), sire and record of dam (model 2), Predicted Difference of sire (model 3), or Predicted Difference of sire and 1/2 Predicted Difference of maternal grandsire (model 4). Ten replicates were performed. Significant cytoplasmic effects were observed for all models in at least 8 of 10 replicates for milk yield and in all replicates for fat percentage. Average F ratios for cytoplasmic effects were 1.43, 1.25, 1.46, and 1.43 for milk yield and 2.00, 1.29, 2.04, and 1.96 for fat percentage for models 1, 2, 3, and 4. Cytoplasmic variances were 1.4 and 3.2% of residual variances for milk yield and fat percentage for model 1. Unaccounted for additive genetic effects can produce spurious cytoplasmic effects.     

Genetic parameters for milk urea nitrogen in relation to milk production traits

The aim of this study was to estimate genetic parameters for test-day milk urea nitrogen (MUN) and its relationships with milk production traits. Three test-day morning milk samples were collected from 1,953 Holstein-Friesian heifers located on 398 commercial herds in the Netherlands. Each sample was analyzed for somatic cell count, net energy concentration, MUN, and the percentage of fat, protein, and lactose. Genetic parameters were estimated using an animal model with covariates for days in milk and age at first calving, fixed effects for season of calving and effect of test or proven bull, and random effects for herd-test day, animal, permanent environment, and error. Coefficient of variation for MUN was 33%. Estimated heritability for MUN was 0.14. Phenotypic correlation of MUN with each of the milk production traits was low. The genetic correlation was close to zero for MUN and lactose percentage (−0.09); was moderately positive for MUN and net energy concentration of milk (0.19), fat yield (0.41), protein yield (0.38), lactose yield (0.22), and milk yield (0.24), and percentage of fat (0.18), and percentage of protein (0.27); and was high for MUN and somatic cell score (0.85). Herd-test day explained 58% of the variation in MUN, which suggests that management adjustments at herd-level can reduce MUN. This study shows that it is possible to influence MUN by herd practice and by genetic selection.     

Lack of evidence of cytoplasmic inheritance in milk production traits of dairy cattle

Heritability estimates, by year of freshening of daughter, were obtained from daughter-dam and granddaughter-granddam regressions using 61,482 triply matched first lactations of artificially sired Holstein cows obtained from the Northeast Dairy Records Processing Laboratory. After adjusting for herd-year-season effects, residual effects may include additive and other genetic effects of the animal, maternal effects, cytoplasmic effects, and other environmental effects. Analysis of residuals showed that cytoplasmic effects accounted for no variation in milk and fat yield and fat percent. Weighted yearly heritability estimates and standard errors from daughter on dam regressions were .35 +/- .01 for milk yield, .30 +/- .01 for milk fat yield, and .63 +/- .01 for milk fat percent and from daughter on granddam regressions were .34 +/- .03 for milk production, .28 +/- .03 for milk fat production, and .55 +/- .03 for milk fat percent. The differences between daughter-dam and daughter-granddam heritability estimates, which estimate twice the fraction of variance due to cytoplasmic effects, were negative and not statistically significant for milk fat yield and also were negative but highly significant for milk fat percent.     

Associations between milk protein polymorphisms and milk production traits.

Associations between milk protein genotypes and milk production traits were estimated from 6803 first lactation records. Exact tests of associated hypotheses and unbiased estimates of genotype effects were from an animal model. Milk protein genotype effects were estimated using a model in which each milk protein gene was analyzed separately (single-gene analysis) and a model in which all milk protein genes were analyzed simultaneously (multigene analysis). The results of the two models indicate that some effects ascribed to certain milk protein genes in the single-gene analysis are not effects of the milk protein gene itself but of linked genes. Results from this study and from literature indicate that the kappa-casein gene or a very closely linked gene affects protein percentage, and the beta-lactoglobulin gene or a very closely linked gene affects fat percentage. Furthermore, effects of beta-casein genotypes on milk production, fat percentage, and protein yield were significant, and beta-lactoglobulin genotypes had significant effects on milk production and protein yield. It is less clear whether those effects are due to effects of milk protein genes themselves or to effects of linked genes.     

Estimation of genetic parameters for milk fat depression in dairy cattle

The objective of this study was to apply reaction norm models to milk recording data to investigate genetic variation in and environmental sensitivity of susceptibility to milk fat depression (MFD). Data comprised 556,276 test-day records of 80,493 heifers in 1043 herds. Breeding values and genetic variances for fat percentage and fat yield were estimated by applying random regression models to average herd-test-day fat percentage. Genetic and permanent environmental correlations between fat yield expressed in different environments ranged, respectively, from 0.83 to 1.00 and from 0.29 to 1.00. Genetic and permanent environmental correlations between fat percentage expressed in different environments ranged, respectively, from 0.87 to 1.00 and from -0.05 to 0.99. Two traits were defined for MFD. The first trait reflected variation of milk fat percentage of animals within lactation after correction for year-season, herd-test-day, age-at-calving, and stage-of-lactation. This trait had an estimated heritability of about 5% and a genetic correlation between the fifth and 95th percentile of the data of 0.50. The second trait reflected the deviation of an animal's fat percentage on a test-day from its expected fat percentage based on fat percentage on the first test-day. This trait had an estimated heritability of about 4% and a genetic correlation between the fifth and 95th percentile of the data of 0.43. The correlation between estimated breeding values of sires for the 2 MFD traits was -0.3. Our results suggest that genetic variation in susceptibility to MFD is present and that selection for reduced susceptibility to MFD is possible.     

Estimating milk, fat, and protein lactation curves with a test day model.

Test day models were used to estimate lactation curves for milk, fat, protein, fat percentage, and protein percentage and to study the influence of age, season, and herd productivity on Holstein lactation curves. Random effects of lactation within herd and fixed effects of herd test date were absorbed. Fixed effects of cow's age on test day and either DIM (57 divisions) by parity (1, 2, greater than or equal to 3) class or season of calving (winter or summer) by DIM by parity class were estimated. Lactation curves for yield traits derived from DIM solutions were flatter for first versus later lactation, even without addition of age effects. Differences between lactation curves for the two seasons were slight, suggesting that most observed seasonal differences are caused by seasonal productivity accounted for by herd test date effects. At peak, winter calving cows yielded slightly more milk of similar fat percentage but of lower protein percentage than those calving in summer. Data were also partitioned into nine subsets based on rolling herd milk and fat percentage. Lactation curves for yield traits, but not percentage traits, varied with rolling herd milk. Lactation curves for fat yield and percentage varied with rolling herd fat percentage.     

Frequency and effect of the bovine acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) K232A polymorphism in Swedish dairy cattle

Acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) is a key enzyme in triacylglycerol synthesis in the mammary gland, and the corresponding gene has emerged as a strong candidate for the variation in milk fat percentage. In this study, the allele frequencies and effects of the DGAT1 K232A variants in the Swedish dairy breeds Swedish Red and Swedish Holstein were investigated. A total of 239 cows, 143 of the Swedish Red breed and 96 of the Swedish Holstein breed, in the experimental herd at the Swedish University of Agricultural Sciences were genotyped for the DGAT1 polymorphism. The Swedish Red cows in the herd belonged to 1 of 2 selection lines with high or low milk fat percentage, respectively, but with similar high total milk energy production. The frequency of the K variant was found to be significantly greater in the high-fat line than in the low-fat line. The average frequency of the K variant in the 2 lines of the Swedish Red cows was 0.09 compared with 0.12 among the Swedish Holstein cows. Mixed model analysis was used to estimate the effect of the DGAT1 K232A polymorphism based on 16,866 test-day records for milk production traits. In accordance with previous studies, the most pronounced effects were found for fat and protein percentages and milk yield; and the K variant was associated with an increase in milk fat and protein percentages but less milk yield compared with the A variant. Less pronounced effects were found for yields of fat and protein for which the K variant was associated with greater fat yield but less protein yield.     

Genetic parameters of milk yield and composition and their relationships with alternative breeding goals

Milk records on 32,077 Holstein cows in first lactation from 1116 herds on the Quebec Dairy Herd Analysis Service between 1973 and 1983 were used to estimate heritabilities of and genetic correlations between milk, fat, protein, and carrier (water, lactose, ash) yields; fat, protein, and carrier percentages; and ratios of protein to fat and protein plus fat to carrier. The cows were daughters of 451 sires; each sire was represented by five or more daughters. Heritabilities were from .20 to .31 for yields, .59 to .61 for percentages, and .58 to .61 for ratio measures. The heritability of protein yield was considerably lower than other yield traits (.20 vs. .28 to .31). The expected effects of multiple trait selection under 10 alternative breeding goals, using different sets of relative economic values, were considered using the genetic parameter estimates. Selection under the existing pricing system (milk yield plus fat percentage differential) had similar expected genetic gains for fat and protein yields as would be achieved under goals that either gave weight to milk yield and fat yield or to milk yield and the ratio of protein plus fat to carrier. Inclusion of information on protein percentage in an index gave relative little additional genetic progress in most traits, except some increases in protein yield and protein percentage were expected to be greater.     

Genetic relationships and additive genetic variation of productive and reproductive traits in Guernsey dairy cattle

Data were first lactation production and reproduction records initiated from 1958 to 1981 in two experiment station Guernsey herds. Heritability estimates using paternal half sib groups were .24 +/- .12 for milk yield, .27 +/- .12 for fat yield, and .77 +/- .15 for fat percentage. Heritability estimates for reproductive traits ranged from .01 to .04 for number of services, service period, conception rate, and days open, but were higher for days in milk at first breeding (.12) and age at first calving (.13). Except for age at first calving, coefficients of additive genetic variation were larger for reproductive traits than for productive traits. Genetic correlations between measures of production and reproduction were moderate to large and antagonistic, except that the relationship between production and age at first calving was favorable. Breeding value estimates for milk yield and reproduction were negatively correlated for sires with above average breeding values for milk yield. Huge phenotypic variances for reproductive traits masked substantial additive genetic variation for these traits. When all things are considered it seems unwise to ignore reproductive performance in selection programs for dairy cattle.     

Prenatal maternal effects on body condition score, female fertility, and milk yield of dairy cows

In this study, maternal effects were described as age of dam at first and second calving, first-lactation body condition score (BCS) of the dam during gestation, and milk yield of the dam. The impact of these effects on first-lactation daughter BCS, fertility, and test-day milk yield was assessed. The effect of milk yield of dam on daughter 305-d yield in the latter's first 3 lactations was also investigated. The proportion of total phenotypic variance in daughter traits accounted for by maternal effects was calculated. Dams calving early for the first time (18 to 23 mo of age) had daughters that produced 4.5% more first-lactation daily milk, had 7% higher BCS, and had their first service 3 d earlier than cows whose dams calved late (30 to 36 mo). However, daughters of dams that calved early had difficulties conceiving as they needed 7% more inseminations and had a 7.5% higher return rate. Cows from second calvings of relatively young (36 to 41 mo) dams produced 6% more first-lactation daily milk, had 2% higher BCS, and showed a significantly better fertility profile than cows whose dams calved at a late age (47 to 55 mo). High maternal BCS during gestation had a favorable effect on daughter BCS, nonreturn rate, and number of inseminations per conception. However, it was also associated with a small decrease in daughter daily milk yield. Changes in dam BCS during gestation did not affect daughter performance significantly. Maternal effects of milk yield of the dam, expressed as her permanent environment during lactation, adversely affected daughter 305-d milk, fat, and protein yield. However, although the effect was significant, it was practically negligible (<0.3% of the mean). Finally, overall maternal effects accounted for a significant proportion of the total phenotypic variance of calving interval (1.4 ± 0.6%) and nonreturn rate (1.1 ± 0.5%).     

Additive and nonadditive genetic effects on milk production in dairy cattle: evidence for major individual heterosis

Coefficients for individual and maternal breed composition and the expected contributions of individual and maternal heterosis and breed source of cytoplasm were assigned to 42,554 primiparous Holstein-Friesian, Jersey, and crossbred cows. The individual additive genetic breed effect influenced all milk production traits. Highly significant maternal additive genetic breed effects equivalent to 3% of the mean were identified for milk yield and milk fat percentage. Individual heterosis was highly significant for milk yield and milk fat yield. A primiparous first cross cow produced 6.1% more milk and 7.2% more milk fat than the average of straightbred cows of both breeds. For milk fat yield, the individual heterosis effect was higher than the individual additive genetic breed difference between Jersey and Holstein-Friesian. A small negative maternal heterosis and a small effect of breed source of cytoplasm were estimated for milk fat percentage. Results suggest that individual heterosis is a major genetic effect for milk yield and milk fat yield. This heterosis could be utilized through a stratified breeding scheme in which high genetic merit nucleus herds maintain genetic progress in the two straightbred populations, and commercial dairy herds employ a rotational cross-breeding scheme to take advantage of both the additive genetic progress and nonadditive genetic effects.     

Short communication: Genetic parameters of individual fatty acids in milk of Canadian Holsteins

The objective of the present study was to estimate heritabilities of milk fatty acids (FA) and genetic and phenotypic correlations among milk FA and milk production traits in Canadian Holsteins. One morning milk sample was collected from each of 3,185 dairy cows between February and June 2010 from 52 commercial herds enrolled in Valacta (Ste-Anne-de-Bellevue, Quebec, Canada). Individual FA percentages (g/100 g of total FA) were determined for each sample by gas chromatography. After editing the data, 2,573 cows from 46 herds remained. Genetic parameters were estimated using multitrait animal models fitted under REML. The model included fixed effects of age at calving and stage of lactation each nested within parity and random effects of herd-year-season of calving, animal, and residual. The pedigree of animals with data was traced back 5 generations on both the male and female sides to account for relationships among animals. The estimates of heritability for individual FA ranged from 0.01 to 0.39, with standard errors ranging from 0.01 to 0.06. Generally, monounsaturated FA (MUFA) and saturated FA (SFA) showed higher heritability estimates than polyunsaturated FA (PUFA). Overall, SFA were negatively genetically correlated with MUFA and PUFA, whereas genetic correlations between MUFA and PUFA were positive. The SFA showed positive associations with fat yield and fat percentage, whereas unsaturated FA were negatively associated with fat yield and fat percentage. Bovine milk FA composition could be improved through genetic selection, and selection for MUFA or against SFA could alter the bovine milk fat profile in a desirable direction.     

Genetic and environmental relationships between body condition score and milk production traits in Canadian Holsteins

The objective of this research was to estimate genetic parameters of first-lactation body condition score (BCS), milk yield, fat percentage (Fat%), protein percentage (Prot%), somatic cell score (SCS), milk urea nitrogen (MUN), lactose percentage (Lact%), and fat to protein ratio (F:P) using multiple-trait random regression animal models. Changes in covariances between BCS and milk production traits on a daily basis have not been investigated before and could be useful for determining which BCS estimated breeding values (EBV) might be practical for selection in the future. Field staff from Valacta milk recording agency (Sainte-Anne-de-Bellevue, QC, Canada) collected BCS from Québec herds several times per cow throughout the lactation. Average daily heritabilities and genetic correlations among the various traits were similar to literature values. On an average daily basis, BCS was genetically unfavorably correlated with milk yield (i.e., increased milk yield was associated with lower body condition). The unfavorable genetic correlation between BCS and milk yield became stronger as lactation progressed, but was equivalent to zero for the first month of lactation. Favorable genetic correlations were found between BCS with Prot%, SCS, and Lact% (i.e., greater BCS was associated with greater Prot%, lower SCS, and greater Lact%). These correlations were strongest in early lactation. On an average daily basis, BCS was not genetically correlated with Fat% or MUN, but was negatively correlated with F:P. Furthermore, BCS at 5 and 50 d in milk (DIM) had the most favorable genetic correlations with milk production traits over the lactation (at 5, 50, 150, and 250 DIM). Thus, early lactation BCS EBV shows potential for selection. Regardless, this study showed that the level of association BCS has with milk production traits is not constant over the lactation. Simultaneous selection for both BCS and milk production traits should be considered, mainly due to the unfavorable genetic correlation between BCS with milk yield.     

Sodium sesquicarbonate for early lactation dairy cows fed corn silage-based diets

An experiment was conducted to evaluate the effects of dietary sodium sesquicarbonate on milk production and composition, feed intake, ruminal fermentation, and urine pH and composition in 32 multiparous Holstein cows in early lactation. Cows were fed high energy diets of 50% corn silage: 50% concentrate containing 0 or .75% sodium sesquicarbonate (total ration DM basis) for 84 d postpartum. Sodium sesquicarbonate increased milk fat percentage (3.46 vs. 3.27%), milk fat yield (1.39 vs. 1.26 kg), and 4% FCM yield (37.1 vs. 34.7 kg) but did not affect milk yield, milk protein percentage or yield, feed intake, ruminal fluid pH, or ruminal VFA profile. Urine composition was not affected by sodium sesquicarbonate, but urine pH (8.14 vs. 8.04) and total urine output as estimated by creatinine concentration were higher than that of controls. In summary, sodium sesquicarbonate, added as .75% of the ration DM, was effective in improving milk fat percentage, milk fat yield, and 4% FCM yield in early lactation dairy cows fed corn silage-based diets.     

Lactation curves of Sarda breed goats estimated with test day models

Test day records of milk yield (38,765), fat and protein contents (11,357) of Sarda goats (the most numerous Italian goat breed) were analysed with mixed linear models in order to estimate the effects of test date (month and year of kidding for fat and protein contents) parity, number of kids born, altitude of location of flocks (<200 m asl, 200-500 m asl, >500 m asl), flocks within altitude and lactation stage (eight days-in-milk intervals of 30 d each) on milk production. All factors considered in the models affected milk traits significantly. Milk yield was lower in first parity goats than in higher parities whereas fat and protein contents showed an opposite trend. Goats with two kids at parturition had a higher milk yield than goats with one kid and tended to have lower fat and protein percentages. Repeatability between test days within lactation was 0.34, 0.17 and 0.45 for milk yield, fat content and protein content, respectively. Lactation curves of goats farmed at different altitudes were clearly separated, especially for milk yield. Results of the present study highlight differences in milk production traits among the three subpopulations that have been previously identified within the Sarda breed on the basis of the morphological structure of animals and altitude of location of flocks.     

Chromosome substitution effects associated with kappa-casein and beta-lactoglobulin in Holstein cattle.

Allelic variants of kappa-casein and beta-lactoglobulin were used to estimate chromosome substitution effects on transmitting abilities for yield traits among sons of two Holstein sires. Sire 1 was a heterozygote for both loci, whereas sire 2 was a heterozygote for beta-lactoglobulin only. Transmitting abilities for fat percentage and the exchange of protein and fat (protein yield minus fat yield) were influenced by kappa-casein genotype among offspring of sire 1. The B allele of kappa-casein was associated with a decrease in fat percentage, whereas fat yield was not affected; the result was 4.5 kg more transmitting ability for protein than for fat. Favorable chromosome substitution effects on transmitting abilities for fat percentage, protein yield, and the exchange of protein and fat were also observed for the A allele of beta-lactoglobulin among offspring of sire 1; there were consistent but not significant estimates for sire 2. These results suggest that different estimates could be obtained from offspring of different parents, perhaps from linkage disequilibrium between families. The chromosome substitution effects appeared to be in the direction opposite to the average correlation of milk fat and protein content; marker-assisted selection for increased yield of protein would result in decreased fat yield.     

Modeling heat stress effect on Holstein cows under hot and dry conditions: Selection tools

Data from milk recording of Holstein-Friesian cows together with weather information from 2 regions in Southern Spain were used to define the models that can better describe heat stress response for production traits and somatic cell score (SCS). Two sets of analyses were performed, one aimed at defining the population phenotypic response and the other at studying the genetic components. The first involved 2,514,762 test-day records from up to 5 lactations of 128,112 cows. Two models, one fitting a comfort threshold for temperature and a slope of decay after the threshold, and the other a cubic Legendre polynomial (LP) model were tested. Average (T_AVE) and maximum daily temperatures were alternatively considered as covariates. The LP model using T_AVE as covariate showed the best goodness of fit for all traits. Estimated rates of decay from this model for production at 25 and 34°C were 36 and 170, 3.8 and 3.0, and 3.9 and 8.2 g/d per degree Celsius for milk, fat, and protein yield, respectively. In the second set of analyses, a sample of 280,958 test-day records from first lactations of 29,114 cows was used. Random regression models including quadratic or cubic LP regressions (TEM_) on T_AVE or a fixed threshold and an unknown slope (DUMMY), including or not cubic regressions on days in milk (DIM3_), were tested. For milk and SCS, the best models were the DIM3_ models. In contrast, for fat and protein yield, the best model was TEM3. The DIM3DUMMY models showed similar performance to DIM3TEM3. The estimated genetic correlations between the same trait under cold and hot temperatures (ρ) indicated the existence of a large genotype by environment interaction for fat (ρ = 0.53 for model TEM3) and protein yield (ρ around 0.6 for DIM3TEM3) and for SCS (ρ = 0.64 for model DIM3TEM3), and a small genotype by environment interaction for milk (ρ over 0.8). The eigendecomposition of the additive genetic covariance matrix from model TEM3 showed the existence of a dominant component, a constant term that is not affected by temperature, representing from 64% of the variation for SCS to 91% of the variation for milk. The second component, showing a flat pattern at intermediate temperatures and increasing or decreasing slopes for the extremes, gathered 15, 11, and 24% of the variation for fat and protein yield and SCS, respectively. This component could be further evaluated as a selection criterion for heat tolerance independently of the production level.     

Genetic analysis of milk urea nitrogen and lactose and their relationships with other production traits in Canadian Holstein cattle

The objective of this research was to estimate heritabilities of milk urea nitrogen (MUN) and lactose in the first 3 parities and their genetic relationships with milk, fat, protein, and SCS in Canadian Holsteins. Data were a random sample of complete herds (60,645 test day records of 5,022 cows from 91 herds) extracted from the edited data set, which included 892,039 test-day records of 144,622 Holstein cows from 4,570 herds. A test-day animal model with multiple-trait random regression and the Gibbs sampling method were used for parameter estimation. Regression curves were modeled using Legendre polynomials of order 4. A total of 6 separate 4-trait analyses, which included MUN, lactose, or both (yield or percentage) with different combinations of production traits (milk, fat and protein yield, fat and protein percentages, and somatic cell score) were performed. Average daily heritabilities were moderately high for MUN (from 0.384 to 0.414), lactose kilograms (from 0.466 to 0.539), and lactose percentage (from 0.478 to 0.508). Lactose yield was highly correlated with milk yield (0.979). Lactose percentage and MUN were not genetically correlated with milk yield. However, lactose percentage was significantly correlated with somatic cell score (−0.202). The MUN was correlated with fat (0.425) and protein percentages (0.20). Genetic correlations among parities were high for MUN, lactose percentage, and yield. Estimated breeding values (EBV) of bulls for MUN were correlated with fat percentage EBV (0.287) and EBV of lactose percentage were correlated with lactation persistency EBV (0.329). Correlations between lactose percentage and MUN with fertility traits were close to zero, thus diminishing the potential of using those traits as possible indicators of fertility.     

Effect of heat stress on production of Mediterranean dairy sheep

A study on heat stress in Mediterranean dairy sheep was undertaken with the objective to examine the relationship between milk production and heat stress, to estimate the additive genetic variances of milk production traits and heat tolerance, and to investigate the possibility of future selection for increased heat tolerance. Production data included 59,661 test-day records belonging to 6624 lactations of 4428 lactating ewes from 17 flocks collected from 1994 through 2003. The traits investigated were daily milk yield, fat and protein percentage, and daily yield of fat-plus-protein. The pedigree file consisted of 5306 animals; in addition to the 4428 animals with records, 188 male and 690 female ancestors were included. Heat stress was modeled by using data from a weather station. Apart from the effects of the weather conditions of the milk recording test-day, the effects of the preceding 1, 2, and 3 d were determined. Because longer periods of heat stress might have a more severe effect than shorter periods, 2-, 3-, and 4-d periods were also considered, by averaging the weather data measurements. Fixed regression analyses were based on models that included effects of flock nested within year of test-day, DIM (days in milk) class x parity class, and several types of weather indicators. The preferred model using the temperature-humidity index (THI) gave a smoother pattern than did the model with temperature x humidity interaction. Both daily milk and fat-plus-protein yield appeared to decrease at THI > or = 23, in all periods considered. Based on the 4-d period, yield decreased for each unit increase of THI above 23 [-62.8 g/unit (-4.2%) for daily milk yield and -8.9 g/unit (-4.9%) for daily fat-plus-protein yield]. Fat and protein percentages appeared to be unaffected by heat stress. A test-day repeatability model was applied for estimation of genetic parameters. The genetic correlations between the general additive effect and the additive effect of heat tolerance were negative (approximately -0.8) for both daily milk and fat-plus-protein yields in all periods considered. Therefore, milk yield is antagonistic with heat tolerance, and selection only for increased milk production will reduce heat tolerance.