Genetic correlations among female fertility, 305-day milk yield and persistency during the first three lactations of Japanese Holstein cows

Genetic parameters for milk yield and reproductive traits in Honduran Holstein cattle

The aim of this study was to estimate the (co)variance components and genetic parameters for milk yield adjusted to 305d (MY305), calving-to-conception interval (CCI), number of services per conception (NSC) and calving interval (CI) of Honduran Holstein cows, by fitting a bivariate animal model using Maximum Restricted Likelihood procedures. Model included the fixed effects of calving number, the contemporary calving group (farm-season-year of calving and the cow age as covariate). The estimated means and standard deviations for MY, CCI, NSC and CI were, 5098.60 ± 1564.32 kg, 168.27 ± 104.71 days, 2.46 ± 1.69 services, and 448.73 ± 109.16 days, respectively; and their estimated heritabilities were 0.21 ± 0.05, 0.03 ± 0.028, 0.02 ± 0.024 and 0.06 ± 0.04, respectively. The genetic correlations between MY305 and CCI, NSC and CI were positive and antagonist, with values of 0.64 ± 0.52, 0.99 ± 0.56, and 0.32 ± 0.24 respectively. Even though moderate to low heritability was estimated for MY305, systematic selection for milk yield, with a reduction in reproductive efficiency, if considered as the only selection criterion is important to be considered. By including reproductive traits and considering permanent environment effects into the breeding program, might yield a slow, but constant and permanent improvement over time.

Effect of temperature-humidity index on the evolution of trade-offs between fertility and production in dairy cattle

BACKGROUND

In the current context of climate change, livestock production faces many challenges to improve the sustainability of systems. Dairy farming, in particular, must find ways to select animals that will be able to achieve sufficient overall production while maintaining their reproductive ability in environments with increasing temperatures. With future forecasted climate conditions in mind, this study used data from Holstein and Montbeliarde dairy cattle to: (1) estimate the genetic-by-temperature-humidity index (THI) interactions for female fertility, and (2) evaluate the production-fertility trade-off with increasing values of THI.

RESULTS

Two-trait random regression models were fitted for conception rate (fertility) and test-day protein yield (production). For fertility, genetic correlations between different THI values were generally above 0.75, suggesting weak genotype-by-THI interactions for conception rate in both breeds. However, the genetic correlations between the conception rate breeding values at the current average THI (THI = 50, corresponding to a 24-h average temperature of 8 °C at 50% relative humidity) and their slopes (i.e., potential reranking) for heat stress scenarios (THI > 70), were different for each breed. For Montbeliarde, this correlation tended to be positive (i.e., overall the best reproducers are less affected by heat stress), whereas for Holstein it was approximately zero. Finally, our results indicated a weak antagonism between production and fertility, although for Montbeliarde this antagonism intensified with increasing THI.

CONCLUSIONS

Within the range of weather conditions studied, increasing temperatures are not expected to exacerbate the fertility-production trade-off. However, our results indicated that the animals with the best breeding values for production today will be the most affected by temperature increases, both in terms of fertility and production. Nonetheless, these animals should remain among the most productive ones during heat waves. For Montbeliarde, the current selection program for fertility seems to be adequate for ensuring the adaptation of fertility traits to temperature increases, without adverse effects on production. Such a conclusion cannot be drawn for Holstein. In the future, the incorporation of a heat tolerance index into dairy cattle breeding programs would be valuable to promote the selection of animals adapted to future climate conditions.

Genomic Predictions in Korean Hanwoo Cows: A Comparative Analysis of Genomic BLUP and Bayesian Methods for Reproductive Traits

This study aimed to predict the accuracy of genomic estimated breeding values (GEBVs) for reproductive traits in Hanwoo cows using the GBLUP, BayesB, BayesLASSO, and BayesR methods. Accuracy estimates of GEBVs for reproductive traits were derived through fivefold cross-validation, analyzing a dataset comprising 11,348 animals and employing an Illumina Bovine 50K SNP chip. GBLUP showed an accuracy of 0.26 for AFC, while BayesB, BayesLASSO, and BayesR demonstrated values of 0.28, 0.29, and 0.29, respectively. For CI, GBLUP attained an accuracy of 0.19, whereas BayesB, BayesLASSO, and BayesR scored 0.21, 0.24, and 0.25, respectively. The accuracy for GL was uniform across GBLUP, BayesB, and BayesR at 0.31, whereas BayesLASSO showed a slightly higher accuracy of 0.33. For NAIPC, GBLUP showed an accuracy of 0.24, while BayesB, BayesLASSO, and BayesR recorded 0.22, 0.27, and 0.30, respectively. The variation in genomic prediction accuracy among methods indicated Bayesian approaches slightly outperformed GBLUP. The findings suggest that Bayesian methods, notably BayesLASSO and BayesR, offer improved predictive capabilities for reproductive traits. Future research may explore more advanced genomic approaches to enhance predictive accuracy and genetic gains in Hanwoo cattle breeding programs.

Identification of Genomic Instability in Cows Infected with BVD Virus

An important factor for dairy cattle farmers is the profitability of cattle rearing, which is influenced by the animals' health and reproductive parameters, as well as their genomic stability and integrity. Bovine viral diarrhea (BVD) negatively affects the health of dairy cattle and causes reproductive problems. The aim of the study was to identify genomic instability in cows with reproductive disorders following infection with the BVD virus. The material for analysis was peripheral blood from Holstein-Friesian cows with reproductive problems, which had tested positive for BVD, and from healthy cows with no reproductive problems, which had tested negative for BVD. Three cytogenetic tests were used: the sister chromatid exchange assay, fragile sites assay, and comet assay. Statistically significant differences were noted between the groups and between the individual cows in the average frequency of damage. The assays were good biomarkers of genomic stability and enabled the identification of individuals with an increased frequency of damage to genetic material that posed a negative impact on their health. The assays can be used to prevent disease during its course and evaluate the genetic resistance of animals. This is especially important for the breeder, both for economic and breeding reasons. Of the three assays, the comet assay proved to be the most sensitive for identifying DNA damage in the animals.

Prediction of persistency for day 305 of lactation at the moment of the insemination decision

When deciding on the voluntary waiting period of an individual cow, it might be useful to have insight into the persistency for the remainder of that lactation at the moment of the insemination decision, especially for farmers who consider persistency in their reproduction management. Currently, breeding values for persistency are calculated for dairy cows but, to our knowledge, prediction models to accurately predict persistency at different moments of insemination are lacking. This study aimed to predict lactation persistency for DIM 305 at different insemination moments (DIM 50, 75, 100, and 125). Available cow and herd level data from 2005 to 2022 were collected for a total of 20,508 cows from 85 herds located in the Netherlands and Belgium. Lactation curve characteristics were estimated for every daily record using the data up to and including that day. Persistency was defined as the number of days it takes for the milk production to decrease by half during the declining stage of lactation, and calculated from the estimated lactation curve characteristic 'decay'. Four linear regression models for each of the selected insemination moment were built separately to predict decay at DIM 305 (decay-305). Independent variables included the lactation curve characteristics at the selected insemination moment, daily milk yield, age, calving season, parity group and other herd variables. The average decay-305 of primiparous cows was lower than that of multiparous cows (1.55 *10-3 vs. 2.41*10-3, equivalent to a persistency of 447 vs. 288 days, respectively). Results showed that our models had limitations in accurately predicting persistency, although predictions improved slightly at later insemination moments, with R2 values ranging between 0.27 and 0.41. It can thus be concluded that, based only on cow and herd milk production information, accurate prediction of persistency for DIM 305 is not feasible.

Integrating genome-wide methylation and transcriptome-wide analyses to reveal the genetic mechanism of milk traits in Kazakh horses

Horse Milk has important quantitative characteristics and high economic value. However, the DNA methylation regulators involved in horse milk traits have not been clarified. To explore the important role of genome-wide DNA methylation in regulating equine milk yield, this study systematically investigated the genome-wide DNA methylation profiles of Kazakh horse blood by comparing a high-production group (HP, average daily milk yield of 7.5 kg) and low-production group (LP, average daily milk yield of 3.2 kg) using deep whole-genome bisulfite sequencing. First, both groups showed similar proportions of methylation at CpG sites. Subsequently, we identified 26,677 differential methylated regions (DMRs) of CG, 15 DMRs of CHG, 480 DMRs of CHH and 8268 DMR-related genes (DMGs). GO and KEGG analyses revealed that some DMGs were involved in regulating milk and milk component formation. By combining the WGBS-seq and the previous RNA-seq data, a total of 94 overlapping genes were obtained. Finally, we found that 9 DMGs are likely involved in milk production by Kazakh horses.

Assessment of Genetic and Health Management of Tunisian Holstein Dairy Herds with a Focus on Longevity

In Tunisia, the recognition of the possibility of including longevity and disease resistance in dairy cattle selection objectives has been hypothesized as a useful strategy by both researchers and producers. However, in this paper, the state of the art, with a focus on health and longevity, is reviewed. Along the same lines, the heritability for the milk traits, fertility traits, and longevity of Tunisian Holstein dairy cows complies with the literature. Therefore, the influence of genetics on some diseases of the dairy cow was investigated. In addition, a decreasing efficiency in cow fertility has been observed over the last few years. The results showed that the risk of culling increased with common diseases. When analyzed with the Weibull model, functional lifespan was strongly influenced by milk yield; therefore, the risk increased with a reduced milk yield. In her first three lactations, the relative risk of selection increased gradually with lactation. Thus, the risk of thinning is highest at the beginning and end of the first feeding and the end of her second feeding. In conclusion, the risk of culling was reduced in parity. The factors that influence the life of the herd, such as health, husbandry, environmental conditions, and management, are often ignored when evaluating longevity.

Estimation of Genetic Parameters for Heifer and Cow Fertility Traits Derived from On-Farm AI Service Records of South African Holstein Cattle

This study aimed to derive additional fertility traits from service data and estimate their genetic parameters for their possible inclusion in the South African Holstein cattle breeding programs. Service records (n = 64,464) were collected from 18 South African Holstein herds using on-farm automated milk recording systems. Using a multivariate model, the data were used to estimate heritabilities and correlations among several fertility traits. The pedigree data consisted of information on 18,592 animals born between 1981 and 2013. Heritability estimates observed were low to moderate (0.02 ± 0.00 to 0.24 ± 0.00), indicating that there is some genetic basis for the explored fertility traits to warrant selection. The genetic correlations observed between fertility traits were generally favorable, with some high correlations between age at first service (AFS) and services per conception for heifers (SPCh) (0.73 ± 0.00) and between days from calving to first service (CFS) and services per conception for cows (SPC) (0.90 ± 0.01). Positive genetic correlations indicate that improvement in one trait is coupled with a correlated genetic increase in another trait. The studied fertility traits could be used in addition to AFC and CI to serve as a basis for the selection of reproduction in dairy cattle to minimize selection bias.

Random Regression Model for Genetic Evaluation and Early Selection in the Iranian Holstein Population

Simple Summary

The objective of this study was to use a model to predict breeding values for sires and cows at an early stage of the first lactation of cows and progeny groups in the Iranian Holstein population to support the early selection of sires. Our results show that we can select sires according to their daughters’ early lactation performance before they finish first lactation. Cross-validation results show that early selection accuracy can be high, and such an early selection can decrease the generation interval and lead to an increased genetic gain in the Iranian Holstein population.

Abstract

The objective of this study was to use a model to predict breeding values for sires and cows at an early stage of the first lactation of cows and progeny groups in the Iranian Holstein population to enable the early selection of sires. An additional objective was to estimate genetic and phenotypic parameters associated with this model. The accuracy of predicted breeding values was investigated using cross-validation based on sequential genetic evaluations emulating yearly evaluation runs. The data consisted of 2,166,925 test-day records from 456,712 cows calving between 1990 and 2015. (Co)-variance components and breeding values were estimated using a random regression test-day model and the average information (AI) restricted maximum likelihood method (REML). Legendre polynomial functions of order three were chosen to fit the additive genetic and permanent environmental effects, and a homogeneous residual variance was assumed throughout lactation. The lowest heritability of daily milk yield was estimated to be just under 0.14 in early lactation, and the highest heritability of daily milk yield was estimated to be 0.18 in mid-lactation. Cross-validation showed a highly positive correlation of predicted breeding values between consecutive yearly evaluations for both cows and sires. Correlation between predicted breeding values based only on records of early lactation (5–90 days) and records including late lactation (181–305 days) were 0.77–0.87 for cows and 0.81–0.94 for sires. These results show that we can select sires according to their daughters’ early lactation information before they finish the first lactation. This can be used to decrease generation interval and to increase genetic gain in the Iranian Holstein population.

Genomewide Association Analyses of Lactation Persistency and Milk Production Traits in Holstein Cattle Based on Imputed Whole-Genome Sequence Data

Lactation persistency and milk production are among the most economically important traits in the dairy industry. In this study, we explored the association of over 6.1 million imputed whole-genome sequence variants with lactation persistency (LP), milk yield (MILK), fat yield (FAT), fat percentage (FAT%), protein yield (PROT), and protein percentage (PROT%) in North American Holstein cattle. We identified 49, 3991, 2607, 4459, 805, and 5519 SNPs significantly associated with LP, MILK, FAT, FAT%, PROT, and PROT%, respectively. Various known associations were confirmed while several novel candidate genes were also revealed, including ARHGAP35, NPAS1, TMEM160, ZC3H4, SAE1, ZMIZ1, PPIF, LDB2, ABI3, SERPINB6, and SERPINB9 for LP; NIM1K, ZNF131, GABRG1, GABRA2, DCHS1, and SPIDR for MILK; NR6A1, OLFML2A, EXT2, POLD1, GOT1, and ETV6 for FAT; DPP6, LRRC26, and the KCN gene family for FAT%; CDC14A, RTCA, HSTN, and ODAM for PROT; and HERC3, HERC5, LALBA, CCL28, and NEURL1 for PROT%. Most of these genes are involved in relevant gene ontology (GO) terms such as fatty acid homeostasis, transporter regulator activity, response to progesterone and estradiol, response to steroid hormones, and lactation. The significant genomic regions found contribute to a better understanding of the molecular mechanisms related to LP and milk production in North American Holstein cattle.

Genetic parameters for test-day milk yield, lactation persistency, and fertility in low-, medium-, and high-production systems in Kenya

Genetic parameters for test-day milk yield, lactation persistency, and age at first calving (as a fertility trait) were estimated for the first 4 lactations in multiple-breed dairy cows in low-, medium-, and high-production systems in Kenya. Data included 223,285 test-day milk yield records from 11,450 cows calving from 1990 to 2015 in 148 herds. A multivariate random regression model was used to estimate variance and covariance components. The fixed effects in the model included herd, year, and test month, and age as a covariate. The lactation profile over days in milk (DIM) was fitted as a cubic smoothing spline. Random effects included herd, year, and test month interaction effects, genetic group effects, and additive genetic and permanent environmental effects modeled with a cubic Legendre polynomial function. The residual variance was heterogeneous with 11 classes. Consequently, the variance components were varied over the lactation and with the production system. The estimated heritability for milk yield was lower in the low-production system (0.04-0.48) than in the medium- (0.22-0.59) and high-production (0.21-0 60) systems. The genetic correlations estimated between different DIM within lactations decreased as the time interval increased, becoming negative between the ends of the lactations in the low- and medium-production systems. Low (0.05) to medium (0.60) genetic correlations were estimated among first lactation test-day milk yields across the 3 production systems. Genetic correlations between the first lactation test-day milk yield and age at first calving ranged from 0.27 to 0.49, 0 to 0.81, and -0.08 to 0.27 in the low-, medium-, and high-production systems, respectively. Medium to high heritabilities (0.17-0.44) were estimated for persistency, with moderate to high (0.30-0.87) genetic correlations between 305-d milk yield and persistency. This indicates that genetic improvement in persistency would lead to increased milk yield. The low to medium genetic correlations between test-day milk yield between production systems indicate that sires may be re-ranked between production systems. Therefore, we conclude that sires should be selected based on a genetic evaluation within the target production system.

The economic value of days open in Holstein cows in Japan based on simulated changes in conception rate

Increasing herd-life length and culling parity with increasing conception rate (CR) is necessary to increase lifetime profit from dairy cow production. Economic values of days open (DO) were determined by calculating changes in fertility traits, herd-life length, and milk yield when the simulated CR were changed in increments of 1% from -5% to 5% from the basal levels, which were obtained for Hokkaido and regions other than Hokkaido separately. When CR increased, number of artificial inseminations, DO, and milking length decreased. Furthermore, culling parity, and annual milk yield increased. Herd-life length increased in Hokkaido and decreased in the other regions. The economic values of CR were 1,623.8 to 946.8 yen and of DO were 857.4 to 399.0 yen. Relative economic values of milk yield to days open per genetic standard deviation were higher in the other regions than in Hokkaido where the economic effect of selection for DO was expected to be equal to selection for 305-day milk yield and herd life. If the survival rate of multipara cows in the other regions increases, then the economic value of DO would similarly increase in Hokkaido.

Genetic parameters for conception rate and milk production traits within and across Holstein herds with different housing types and feeding systems during the first 3 lactations

The housing types (HST) in which dairy cows are kept and the feeding systems (FDS) used differ among farmers in Japan. Here, we investigated the genetic relationships among conception rate at first insemination (CR) and milk production traits (PROD) during the first 3 lactations of Holstein cows by using a multiple-trait model that considered the trait values of herds with different HST [tiestall (TSL) barn, freestall (FS) barn, or grazing (GZ)] and FDS as separate traits. Milk production and conception records of Holstein cows in the Hokkaido region of Japan (283,611 records for first lactation, 253,902 for second, and 181,197 for third) were analyzed. We categorized herds with TSL or FS into 2 types of FDS for cows: separate feeding (SF) of roughage plus concentrate or feeding of total mixed ration, in which roughage and concentrates were mixed before feeding. The PROD analyzed were cumulative milk, fat, and protein yields within 305 d and lactation persistency, which we defined as the difference between milk yields at 240 and 60 d in milk. We estimated the heritabilities for CR or PROD within each HST or HST × FDS group and the genetic correlations between these traits within each group or across different groups within each lactation by using a 3-HST (TSL, FS, and GZ) × 2-trait (CR and each PROD) or 2-HST (TSL and FS) × 2-FDS × 2-trait animal model. Heritability estimates for CR in GZ were higher than those in TSL or FS, and genetic correlations for CR between GZ and TSL or FS barns were weaker than those between TSL and FS barns. In addition, genetic correlations between CR and PROD in GZ were weaker than those in TSL and FS barns. In the comparison among the 4 HST × FDS except GZ, heritability estimates for CR in FS × SF were higher than those in the others, and genetic correlations for CR between FS × SF and the other systems were relatively weak. These results indicated that differences in the production system for Holstein cows influence genotypic effects in terms of the cows' ability to conceive and the genetic relationships between fertility traits and milk production traits.

Different aspects of lactation persistency in dairy cows

Lactation persistency (cow’s ability to maintain milk production after reaching its peak) is a very important economic characteristic in the dairy cattle production system. Different definition and functions for describing and measuring of this trait were proposed by researchers. The random regression model using Legendre polynomial was one of the common and effective methodologies for evaluation of persistency in the last decade. Several factors affecting persistency such as different characteristics of lactation curve, environment factors, reproduction traits and health status of the dairy cow. Based on different studies the heritability of this trait was low to medium and negative or positive amount of genetic correlation between persistency and total milk yield in dairy cattle is attributed to persistency measures and method of data analysis. Persistency is related with low and later peak yield and selecting cows for peak yield will improve persistency and lactation curve traits. Analysis of relationships between persistency and other functional traits show signs that genetic improvement for persistency is possible and favorable. Different aspects and relationships of persistency with various lactation and other functional traits in dairy cows are reviewed in this article.

Genetic correlations between yield traits or days open measured in cows and semen production traits measured in bulls

We used a bivariate animal model to investigate the genetic correlations between yield traits or days open (DO) as characters measured in cows and semen production traits as characters measured in bulls. Lactation records of 305-day milk, fat, and protein yields, and DO, from 386 809 first-lactation Holstein cows in Hokkaido, Japan, that calved between 2008 and 2014 were used. Semen production records were collected between 2005 and 2014 and included volume per ejaculate (VOL), sperm concentration (CON), number of sperm per ejaculate (NUM), progressive motility index of sperm (MOT), and MOT after freeze-thawing (A-MOT). Number of sperm per ejaculate was log-transformed into a NUM score (NUMS). A total of 30 373 semen production records from 1196 bulls were obtained. The pedigree file used for analysing the records was involving 885 345 animals. Heritability was estimated for VOL (0.42), CON (0.12), NUMS (0.37), MOT (0.08), and A-MOT (0.11). Weak and negative genetic correlations were recorded between yield traits measured in cows and VOL, CON or NUMS measured in bulls. Moderate and negative genetic correlations were obtained between DO and MOT (-0.42) or A-MOT (-0.43). Selection focused on MOT or A-MOT measured in bulls may therefore improve DO measured in cows.

Genetic relationship of lactation persistency with milk yield, somatic cell score, reproductive traits, and longevity in Slovak Holstein cattle

Abstract. The objective of this study was to estimate the breeding values (BVs) of lactation persistency, the test day of milk yield, the somatic cell score, reproductive traits (calving interval, days open), longevity in Slovak Holstein dairy cattle. BVs were used for the detection of relationships among the persistency of lactation and other selected traits. Data for the estimation of BVs of milk production and somatic cell score were collected from 855 240 cows. BVs for reproductive traits were estimated for 352 712 cows and for longevity for 528 362 cows. The highest correlations were confirmed between the BV of persistency and the BV of test day milk yield at 100, 200, and 305 days (−0.88, −0.65, and −0.61). Correlations between the BV of lactation persistency and the BV of somatic cell score at day 305 or the BV of somatic cell score persistency were favorable: −0.05 and −0.12, respectively. The relationship between the BV of persistency and the BV of the calving interval or the BV of days open was 0.11 and 0.10 respectively. The selection for the persistency of lactation may not improve longevity because there is no relation between the BV of persistency and the BV of longevity (rg = 0.06).

Genetic parameters of reproductive traits in Tunisian Holsteins

Abstract. Multi-trait Bayesian procedure was used to estimate genetic parameters for reproductive traits in Tunisian Holstein cows. A total of 31 348 lactations of the calving years 2005 to 2012 were analyzed. Fertility traits were the calving interval (CI), days open (DO), days to first insemination (DFI), days from first insemination to conception (FIC), and number of inseminations per conception (NI). Posterior means of heritabilities of CI, DO, DFI, FIC, and NI were 0.047, 0.03, 0.025, 0.024, and 0.069, respectively. Posterior means of repeatabilities of the same respective traits were 0.106, 0.094, 0.051, 0.036, and 0.17. Genetic correlations among female fertility traits were also computed. Calving interval and DO had the highest genetic correlation estimate (0.85) because they have overlapping genetic meanings. The lowest genetic correlation estimate (−0.25) was found between DFI and NI. Genetic parameter estimates are low and are even lower than those reported in most literature, implying that more focus should be put upon improving the management of reproduction in dairy cattle herds in Tunisia.

Effects of stage of pregnancy on variance components, daily milk yields and 305-day milk yield in Holstein cows, as estimated by using a test-day model

Pregnancy and calving are elements indispensable for dairy production, but the daily milk yield of cows decline as pregnancy progresses, especially during the late stages. Therefore, the effect of stage of pregnancy on daily milk yield must be clarified to accurately estimate the breeding values and lifetime productivity of cows. To improve the genetic evaluation model for daily milk yield and determine the effect of the timing of pregnancy on productivity, we used a test-day model to assess the effects of stage of pregnancy on variance component estimates, daily milk yields and 305-day milk yield during the first three lactations of Holstein cows. Data were 10 646 333 test-day records for the first lactation; 8 222 661 records for the second; and 5 513 039 records for the third. The data were analyzed within each lactation by using three single-trait random regression animal models: one model that did not account for the stage of pregnancy effect and two models that did. The effect of stage of pregnancy on test-day milk yield was included in the model by applying a regression on days pregnant or fitting a separate lactation curve for each days open (days from calving to pregnancy) class (eight levels). Stage of pregnancy did not affect the heritability estimates of daily milk yield, although the additive genetic and permanent environmental variances in late lactation were decreased by accounting for the stage of pregnancy effect. The effects of days pregnant on daily milk yield during late lactation were larger in the second and third lactations than in the first lactation. The rates of reduction of the 305-day milk yield of cows that conceived fewer than 90 days after the second or third calving were significantly (P<0.05) greater than that after the first calving. Therefore, we conclude that differences between the negative effects of early pregnancy in the first, compared with later, lactations should be included when determining the optimal number of days open to maximize lifetime productivity in dairy cows.