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.     

Association of single nucleotide polymorphisms in candidate genes previously related to genetic variation in fertility with phenotypic measurements of reproductive function in Holstein cows

Many genetic markers related to health or production traits are not evaluated in populations independent of the discovery population or related to phenotype. Here we evaluated 68 single nucleotide polymorphisms (SNP) in candidate genes previously associated with genetic merit for fertility and production traits for association with phenotypic measurements of fertility in a population of Holstein cows that was selected based on predicted transmitting ability (PTA) for daughter pregnancy rate (DPR; high, ≥1, n = 989; low, ≤ ?1.0, n = 1,285). Cows with a high PTA for DPR had higher pregnancy rate at first service, fewer services per conception, and fewer days open than cows with a low PTA for DPR. Of the 68 SNP, 11 were associated with pregnancy rate at first service, 16 with services per conception, and 19 with days open. Single nucleotide polymorphisms in 12 genes (BDH2, BSP3, CAST, CD2, CD14, FUT1, FYB, GCNT3, HSD17B7, IBSP, OCLN, and PCCB) had significant associations with 2 fertility traits, and SNP in 4 genes (CSPP1, FCER1G, PMM2, and TBC1D24) had significant associations with each of the 3 traits. Results from this experiment were compared with results from 2 earlier studies in which the SNP were associated with genetic estimates of fertility. One study involved the same animals as used here, and the other study was of an independent population of bulls. A total of 13 SNP associated with 1 or more phenotypic estimates of fertility were directionally associated with genetic estimates of fertility in the same cow population. Moreover, 14 SNP associated with reproductive phenotype were directionally associated with genetic estimates of fertility in the bull population. Nine SNP (located in BCAS, BSP3, CAST, FUT1, HSD17B7, OCLN, PCCB, PMM2, and TBC1D24) had a directional association with fertility in all 3 studies. Examination of the function of the genes with SNP associated with reproduction in more than one study indicates the importance of steroid hormones and immune function as determinants of reproductive function. All but 1 of the 68 evaluated SNP were variable in 11 breeds besides Holstein, indicating the potential effects of these SNP on reproductive function across breeds of cattle.     

Short communication: Iodine content in bovine milk is lowly heritable and shows limited genetic variation

Milk and dairy products are considered important sources of iodine in several countries. Despite this, there is a paucity of studies that have investigated sources of variation of milk iodine, especially on a large scale. So far, it is not clear if milk iodine content could be increased through breeding in dairy cattle. Recently, a mid-infrared spectroscopy prediction model has been developed for an indirect quantification of iodine content in cow milk, as it is a faster and less expensive method that allows the prediction at population level. The model has coefficient of determination and ratio of performance to deviation in external validation of 0.57 and 1.44, respectively, and it was used in the present study to predict the iodine content from historical milk spectral data to investigate phenotypic and genetic aspects in the Italian Holstein cattle. Based on the accuracy of the model, the prediction was interpreted as proxy for the real milk iodine concentration (IOD_P). The data set comprised 33,776 test-day records with IOD_P from 4,072 cows. Data of IOD_P were transformed through natural logarithm to achieve a normal distribution. The effect of parity, lactation stage, and month of sampling were investigated, and genetic parameters were estimated using a test-day repeatability animal model. Milk IOD_P decreased with parities and was the lowest in early lactation. Heritability of IOD_P was low (0.025) and it was positively genetically correlated with milk yield and negatively with fat content. Results suggested that it would be challenging to directly improve this trait through breeding strategies in dairy cattle, because IOD_P is mainly affected by temporary environmental factors and thus, cannot be easily improved through genetics. Although preliminary, findings of this study suggest that it would be more convenient to develop feeding and management strategies to drive milk iodine level than to put efforts and resources into breeding strategies. Further studies should validate IOD_P as an indicator trait of milk iodine content by improving reference data and estimating genetic correlation between predicted and measured values.     

Genetic parameters for milk coagulation properties in Estonian Holstein cows

The objective of this study was to estimate heritabilities and repeatabilities for milk coagulation traits [milk coagulation time (RCT) and curd firmness (E₃₀)] and genetic and phenotypic correlations between milk yield and composition traits (milk fat percentage and protein percentage, urea, somatic cell count, pH) in first-lactation Estonian Holstein dairy cattle. A total of 17,577 test-day records from 4,191 Estonian Holstein cows in 73 herds across the country were collected during routine milk recordings. Measurements of RCT and E₃₀ determined with the Optigraph (Ysebaert, Frepillon, France) are based on an optical signal in the near-infrared region. The cows had at least 3 measurements taken during the period from April 2005 to January 2009. Data were analyzed using a repeatability animal model. There was substantial variation in milk coagulation traits with a coefficient of variation of 27% for E₃₀ and 9% for the log-transformed RCT. The percentage of variation explained by herd was 3% for E₃₀ and 4% for RCT, suggesting that milk coagulation traits are not strongly affected by herd conditions (e.g., feeding). Heritability was 0.28 for RCT and 0.41 for E₃₀, and repeatability estimates were 0.45 and 0.50, respectively. Genetic correlation between both milk coagulation traits was negligible, suggesting that RCT and E₃₀ have genetically different foundations. Milk coagulation time had a moderately high positive genetic (0.69) and phenotypic (0.61) correlation with milk pH indicating that a high pH is related to a less favorable RCT. Curd firmness had a moderate positive genetic (0.48) and phenotypic (0.45) correlation with the protein percentage. Therefore, a high protein percentage is associated with favorable curd firmness. All reported genetic parameters were statistically significantly different from zero. Additional univariate random regression analysis for milk coagulation traits yielded slightly higher average heritabilities of 0.38 and 0.47 for RCT and E₃₀ compared with the heritabilities of the repeatability model.     

A quantitative trait locus on Bos taurus autosome 17 explains a large proportion of the genetic variation in de novo synthesized milk fatty acids

A genomic region associated with milk fatty acid (FA) composition has been detected on Bos taurus autosome (BTA)17 based on 50,000 (50K) single nucleotide polymorphism (SNP) genotypes. The aim of our study was to fine-map BTA17 with imputed 777,000 (777K) SNP genotypes to identify candidate genes associated with milk FA composition. Phenotypes consisted of gas chromatography measurements of 14 FA based on winter and summer milk samples. Phenotypes and genotypes were available on 1,640 animals in winter milk, and on 1,581 animals in summer milk samples. Single-SNP analyses showed that several SNP in a region located between 29.0 and 34.0 Mbp were in strong association with C6:0, C8:0, and C10:0. This region was further characterized based on haplotypes. In summer milk samples, for example, these haplotypes explained almost 10% of the genetic variance in C6:0, 9% in C8:0, 3.5% in C10:0, 1.8% in C12:0, and 0.9% in C14:0. Two groups of haplotypes with distinct predicted effects could be defined, suggesting the presence of one causal variant. Predicted haplotype effects tended to increase from C6:0 to C14:0; however, the proportion of genetic variance explained by the haplotypes tended to decrease from C6:0 to C14:0. This is an indication that the quantitative trait locus (QTL) region is involved either in the elongation process or in early termination of de novo synthesized FA. Although many genes are present in this QTL region, most of these genes on BTA17 have not been characterized yet. The strongest association was found close to the progesterone receptor membrane component 2 (PGRMC2) gene, which has not yet been associated with milk FA composition. Therefore, no clear candidate gene associated with milk FA composition could be identified for this QTL.     

Statistical modeling of candidate gene effects on milk production traits in dairy cattle

A major objective of dairy cattle genomic research is to identify genes underlying the variability of milk production traits that could be useful in breeding programs. The candidate gene approach provides tools for searching for causative polymorphisms affecting quantitative traits. Genes with a possible effect on milk traits in cattle can be involved in different physiological pathways, such as triglyceride synthesis [acyl-CoA:diacylglycerol acyltransferase 1 gene (DGAT1)], fat secretion from the mammary epithelial tissue (butyrophilin), or entire-body energy homeostasis regulation (leptin and leptin receptor). In this study, based on data from 252 Black and White bulls from the active Polish dairy population, effects and potential interactions of 9 single nucleotide polymorphisms in the butyrophilin, DGAT1, leptin, and leptin receptor genes were investigated. Additionally, the effect of the number of additive, dominance, and epistatic genetic effects fitted into the model on the estimates of model parameters and model selection was illustrated. Phenotypic records were daughter yield deviations for milk, fat, and protein yields, obtained from a routine national genetic evaluation. Out of all the analyzed polymorphisms, DGAT1 K232A had a much larger effect on milk traits than the other single nucleotide polymorphisms considered. Estimates of the additive genetic effect of K232A expressed as half of the difference between Lys- and Ala-encoding variants were −107.4 kg of milk, 5.4 kg of fat, and −1.6 kg of protein at first parity, as well as −120 kg of milk and 6.8 kg of fat at second parity. In terms of model selection, it was demonstrated that the modified version of Bayesian information criterion selects models with the parameterization reflecting the genetic background of the analyzed trait, while the Bayesian information criterion chooses models that are too highly parameterized.     

A 100-Year Review: Methods and impact of genetic selection in dairy cattle—From daughter–dam comparisons to deep learning algorithms

In the early 1900s, breed society herdbooks had been established and milk-recording programs were in their infancy. Farmers wanted to improve the productivity of their cattle, but the foundations of population genetics, quantitative genetics, and animal breeding had not been laid. Early animal breeders struggled to identify genetically superior families using performance records that were influenced by local environmental conditions and herd-specific management practices. Daughter–dam comparisons were used for more than 30 yr and, although genetic progress was minimal, the attention given to performance recording, genetic theory, and statistical methods paid off in future years. Contemporary (herdmate) comparison methods allowed more accurate accounting for environmental factors and genetic progress began to accelerate when these methods were coupled with artificial insemination and progeny testing. Advances in computing facilitated the implementation of mixed linear models that used pedigree and performance data optimally and enabled accurate selection decisions. Sequencing of the bovine genome led to a revolution in dairy cattle breeding, and the pace of scientific discovery and genetic progress accelerated rapidly. Pedigree-based models have given way to whole-genome prediction, and Bayesian regression models and machine learning algorithms have joined mixed linear models in the toolbox of modern animal breeders. Future developments will likely include elucidation of the mechanisms of genetic inheritance and epigenetic modification in key biological pathways, and genomic data will be used with data from on-farm sensors to facilitate precision management on modern dairy farms.     

Genome-wide association study and in silico functional analysis of the number of embryos produced by Holstein donors

Superovulation or ovum pick-up and in vitro fertilization are technologies used to produce an increased number of embryos from elite females. Embryo production traits have been shown to be heritable, but the genes that cause this variability have not yet been assessed. The main objectives of this study were to perform a genome-wide association study (GWAS) to find single nucleotide polymorphisms (SNP) associated with embryo production traits and to identify candidate genes affecting the number of embryos produced by Holstein donors in Canada that may provide insight into the regulation of embryo production. Breeding values were estimated and de-regressed for all donors and sires using a data set of 150,971 records of superovulation or ovum pick-up and in vitro fertilization. A total of 11,607 animals were genotyped, but of that number only 5,118 were genotyped with at least a 50K SNP panel and had a de-regressed estimated breeding value reliability of at least 10%. For the GWAS, 606,406 imputed SNP on 29 autosomal chromosomes were considered after applying quality control measures. A single-SNP univariate mixed linear animal model was used to perform the GWAS, and a 5% false discovery rate was applied to adjust for multiple testing. We found 36 and 14 significant SNP associated with the total number of embryos and the number of viable embryos, respectively, with most of them located on chromosome 11. Using these significant SNP, positional genes located within 10,000 bp upstream and downstream of the SNP were retrieved. Thirteen genes were harboring or near the significant SNP for the total number of embryos, 4 of them also being near the significant SNP for viable embryos. Some of these genes (CRB2, DENND1A, MAD1L1, NDUFA8, PTGS1) could be considered as potential positional candidate genes related to the number of embryos produced by a donor. This list will need to be validated in an independent population to confirm the role of the genes for embryo production.     

Quantitative trait loci mapping of functional traits in the German Holstein cattle population

A whole-genome scan to detect quantitative trait loci (QTL) for functional traits was performed in the German Holstein cattle population. For this purpose, 263 genetic markers across all autosomes and the pseudoautosomal region of the sex chromosomes were genotyped in 16 granddaughter-design families with 872 sons. The traits investigated were deregressed breedingvalues for maternal and direct effects on dystocia (DYSm, DYSd) and stillbirth (STIm, STId) as well as maternal and paternal effects on nonreturn rates of 90 d (NR90m, NR90p). Furthermore, deregressed breeding values for functional herd life (FHL) and daughter yield deviation for somatic cell count (SCC) were investigated. Weighted multimarker regression analyses across families and permutation tests were applied for the detection of QTL and the calculation of statistical significance. A ten percent genomewise significant QTL was localized for DYSm on chromosome 8 and for SCC on chromosome 18. A further 24 putative QTL exceeding the 5% chromosomewise threshold were detected. On chromosomes 7, 8, 10, 18, and X/Yps, coincidence of QTL for several traits was observed. Our results suggest that loci with influence on udder health may also contribute to genetic variance of longevity. Prior to implementation of these QTL in marker assisted selection programs for functional traits, information about direct and correlated effects of these QTL as well as fine mapping of their chromosomal positions is required.     

Inbreeding trends and application of optimized selection in the UK Holstein population

Important increases in the rates of inbreeding have recently been observed in dairy cattle populations, and methods have been proposed to address these increases. The aims of this study were to estimate the current level and rates of inbreeding in the UK Holstein population and to investigate the potential of applying optimized selection to manage the rates of inbreeding. Inbreeding coefficients were calculated for the entire UK Holstein population using 1940 as the base year. Rates of inbreeding were obtained for 3 time periods by regressing mean inbreeding coefficients on the year of birth of the animals. The expected average pedigree index and expected inbreeding of offspring using optimized contributions for a given set of selection candidates was compared to the expected pedigree index and inbreeding of offspring for the same set of selection candidates using observed contributions. The rate of inbreeding in the UK Holstein population has increased substantially since 1990 when compared to previous time periods. This increase is most likely due to the large influence of a few related sires on the breed in the mid- to late 1980s. The introduction of the individual animal model in the early 1990s may also have contributed to increased inbreeding. Optimized selection appears to represent a promising selection tool, not only to manage rates of inbreeding, but also to increase genetic gain at the same rate of inbreeding.     

Effects of recombinant bovine interleukin-8 (rbIL-8) treatment on health,metabolism, and lactation performance in Holstein cattle I: Production and functional characterization of rbIL-8 in vitro and in vivo

In the present study, we standardized processes of cloning and purification of recombinant bovine interleukin-8 (rbIL-8) from bacterial culture and assessed its biological activity in Holstein cattle. Plasmid containing a subclone of bovine IL-8 was expressed using Escherichia coli BL21 and cell lysate was purified by chromatography. The presence of rbIL-8 was assessed by Western blot analyses and function was confirmed in vitro using a chemotaxis chamber. Based on optical density values, chemoattractant properties of rbIL-8 were 10-fold greater compared with control wells. Two in vivo studies were conducted to assess the biological activity of rbIL8. For study 1, one-year-old Holstein heifers (n = 20) were randomly allocated to receive a single intravaginal administration containing 1,125 µg of rbIL-8 diluted in 20 mL of saline solution (rbIL-8, n = 10) or a single intravaginal administration of 20 mL of saline solution (control, n = 10). For study 2, nonpregnant lactating Holstein cows (n = 31) were randomly allocated to receive an intrauterine administration with 1,125 µg of rbIL-8 diluted in 20 mL of saline solution (rbIL-8, n = 11), a positive control consisting of resin-purified lysate of E. coli BL21 not transfected with the plasmid coding for rbIL-8 diluted in 20 mL of saline solution (E. coli, n = 10), and a negative control administered with 20 mL of saline solution (control, n = 10). An increase in vaginal neutrophils was observed in heifers treated with rbIL-8 within 3 h of treatment, but not in control heifers. Additionally, intrauterine administration of rbIL-8 increased the proportion of PMN cells in uterine cytological samples from 3.5% before treatment to 75.8% 24 h later—an increase that was not observed in the negative control group and cows treated with resin-purified lysate of E. coli. To further evaluate the effect of local and systemic rbIL-8 stimulation on the dynamics of circulating white blood cells, a third study was conducted. In study 3, nonpregnant 8-mo-old Holstein heifers (n = 30) were randomly allocated into 1 of 3 treatment groups: intravenous rbIL-8 (1,125 µg of rbIL-8 diluted in 5 mL of saline solution, n = 10); intravaginal rbIL-8 (1,125 µg of rbIL-8 diluted in 20 mL of saline solution; n = 10); or intravaginal saline (20 mL of saline solution, n = 10). Intravenous injection of rbIL-8 resulted in a transient increase in rectal temperature, which was greater at 2 h after treatment compared with cows treated intravaginally with rbIL-8 or heifers treated with saline solution. Heifers treated with rbIL-8 intravenously displayed a marked reduction in neutrophils, basophils, lymphocytes, and monocytes within the first 4 h posttreatment compared with heifers treated intravaginally. However, at 6 h after treatment, heifers treated with rbIL-8 intravenously displayed a rebound in white blood cell counts caused by an increase in neutrophil counts. These results show that the presented purification method is effective and results in biologically active rbIL-8 that can be used safely to modulate immune responses in cattle.     

Mapping quantitative trait Loci for left-sided displacement of the abomasum in German Holstein dairy cows

A whole-genome scan using an affected paternal half-sib design was utilized to detect quantitative trait loci (QTL) for left-sided displaced abomasum (LDA) in German Holsteins. A total of 360 animals from 14 paternal half-sib families were genotyped, for a total of 306 polymorphic microsatellites. For a whole-genome scan, 221 markers were equally distributed over all 29 bovine autosomes, with an average distance of 13.7 cM. For fine-mapping, a total of 85 additional microsatellites were used. We identified genome-wide significant QTL on Bos taurus autosome (BTA) 1 (54.6 to 58.3 cM) and on BTA3 (5.9 cM). Furthermore, 3 chromosome-wide significant QTL were located on bovine chromosomes 21, 23, and 24. In addition, we found 11 QTL that cosegregated in grandsire families but that were not significant in the across-family analysis. These QTL were located on BTA5, 6, 10, 12, 15, 16, 17, 19, 23, and 26. This study is the first report on QTL for LDA and is a first step toward identifying single nucleotide polymorphisms for LDA-QTL.     

Genetic progress in multistage dairy cattle breeding schemes using genetic markers

The aim of this paper was to explore general characteristics of multistage breeding schemes and to evaluate multistage dairy cattle breeding schemes that use information on quantitative trait loci (QTL). Evaluation was either for additional genetic response or for reduction in number of progeny-tested bulls while maintaining the same response. The reduction in response in multistage breeding schemes relative to comparable single-stage breeding schemes (i.e., with the same overall selection intensity and the same amount of information in the final stage of selection) depended on the overall selection intensity, the selection intensity in the various stages of the breeding scheme, and the ratio of the accuracies of selection in the various stages of the breeding scheme. When overall selection intensity was constant, reduction in response increased with increasing selection intensity in the first stage. The decrease in response was highest in schemes with lower overall selection intensity. Reduction in response was limited in schemes with low to average emphasis on first-stage selection, especially if the accuracy of selection in the first stage was relatively high compared with the accuracy in the final stage.Closed nucleus breeding schemes in dairy cattle that use information on QTL were evaluated by deterministic simulation. In the base scheme, the selection index consisted of pedigree information and own performance (dams), or pedigree information and performance of 100 daughters (sires). In alternative breeding schemes, information on a QTL was accounted for by simulating an additional index trait. The fraction of the variance explained by the QTL determined the correlation between the additional index trait and the breeding goal trait. Response in progeny test schemes relative to a base breeding scheme without QTL information ranged from +4.5% (QTL explaining 5% of the additive genetic variance) to +21.2% (QTL explaining 50% of the additive genetic variance). A QTL explaining 5% of the additive genetic variance allowed a 35% reduction in the number of progeny tested bulls, while maintaining genetic response at the level of the base scheme. Genetic progress was up to 31.3% higher for schemes with increased embryo production and selection of embryos based on QTL information. The challenge for breeding organizations is to find the optimum breeding program with regard to additional genetic progress and additional (or reduced) cost.     

Factors affecting incorrect paternity assignment in the Israeli Holstein population

A total of 6040 Israeli Holstein cows from 181 Kibbutz herds listed as progeny of 11 sires were genotyped for 104 microsatellites. Seventeen markers were deleted due to a frequency of erroneous genotypes >1%, leaving 160,470 valid genotypes. Conflicts between the putative sire and daughter in at least 2 markers and for at least 10% of the markers genotyped per cow were required to reject paternity. Cows that did not meet the requirements for paternity confirmation or rejection were deleted from further analysis. The frequency of rejected paternity was 11.7%. The effects of recorded sire, birth year, geographical region, herd, and inseminator on the frequency of paternity rejection were analyzed with linear and nonlinear models. Only the effects of inseminator and recorded sire were significant in all models tested that included these effects. The main causes of incorrect paternity recording appear to be inseminator recording mistakes, and possibly mistakes with respect to semen labeling at the AI institutes. Incorrect paternity recording due to multiple inseminations by different sires could explain, at most, 20% of the paternity mistakes. Instituting a system of quality control, especially at the level of the inseminator, should reduce paternity errors to no more than 8%, and increase genetic progress by at least 1%.     

烟草株高和叶数性状QTL定位及候选基因预测

为明确烟草株高和叶数性状的遗传规律,发掘控制相关性状的主效位点,以台烟8号（TT8）为母本（P1）,NC82为父本（P2）,配置杂交组合,以TT8为轮回亲本回交,构建了包含200个单株的BC1F1群体。在此基础上,分别在四川、山东两个环境点种植群体材料,获得株高和叶数表型,进而利用高密度遗传图谱对株高和叶数性状进行QTL定位分析。结果表明,在两个环境条件下共定位到4个与株高和叶数相关的主效QTL,每个QTL均可以解释相应性状10%~20%的表型变异。两个环境点定位到2个叶数性状QTL,均位于23号连锁群上,非等位;定位到3个株高性状主效QTL,四川点1个位于23号连锁群上,山东点2个分别位于21号和23号连锁群上,其中山东点23号连锁群上的主效位点与控制叶数的主效QTL遗传位置一致。相关结果为进一步克隆控制株高和叶数性状的主效基因及烟草重要农艺性状分子改良奠定了基础。     

Genetic correlations among production, body size, udder, and productive life traits over time in Holsteins

Genetic correlations among milk, fat, and protein yields; body size composite (BSC); udder composite (UDC); and productive life (PL) in Holsteins were investigated over time. The data set contained 25,280 records of cows born in Wisconsin between 1979 and 1993. The multiple trait random regression (MT-RR) animal model included registration status, herd-year, age group, and stage of lactation as fixed effects; additive genetic effects with random regressions (RR) on year of birth using the first-order Legendre polynomial; and residual effects. Heterogeneous residual variances were considered in the model. Estimates of variance components and genetic correlations among traits from MT-RR were compared with those estimated with a multiple trait interval (MT-I) model, which assumed that every 3-yr interval was a separate trait and included the same effects as in the MT-RR model except for the RR. Genetic correlations estimated with MT-RR and MT-I models over time among all traits were compared with correlations among breeding values predicted with the single trait (ST) model without RR. Correlations among breeding values predicted with MT-RR, ST, and MT models were also calculated. Additive genetic and residual variances for all traits except PL increased over time; those for PL were constant. As a result, heritability estimates had no significant changes during the 15 yr. Genetic correlations of PL with milk, fat, protein, and BSC declined to zero or negative; those with UDC remained positive. Correlations among breeding values predicted with ST, MT, and MT-RR models were relatively high for all traits except PL. Genetic correlations between PL and other traits varied over time, with some correlations changing sign. For accurate indirect prediction of PL from other traits, the genetic correlations among the traits need to be re-estimated periodically.