Variances and correlations of milk production, fertility, longevity, and type traits over time in Australian Holstein cattle

Telomere Length in Neonatal Dairy Calves in Relation to Lifetime Parameters

Telomere length (TL) has gained attention as a biomarker for longevity and productivity in dairy cattle. This study explored the association between neonatal TL in Holstein calves and lifetime parameters (lifespan, milk production, and reproduction). Blood samples were collected from 210 calves (≤10d old) across four dairy farms in Flanders, Belgium. Telomere length was measured using qPCR and analyzed as a continuous variable and across three groups: the 10% shortest, the 10% longest, and the remaining 80%. Survival analyses showed no association between TL and lifespan (p = 0.1) or TL groups (p = 0.8). Similarly, TL showed no significant association with production traits. However, categorical analyses revealed that calves with the longest TL had lower lifetime fat (p = 0.01) and protein yields (p = 0.01) than those with the shortest TL. Reproductive analyses showed cows in the long TL group required fewer inseminations per lactation (p = 0.02) and exhibited longer calving intervals (p = 0.05). These findings suggest that while neonatal TL may not predict productive lifespan, it may provide insight into reproductive efficiency. Future studies should prioritize longitudinal assessments of TL dynamics to better understand their interactions with management practices and application in herd improvement.

Archetypal clustering reveals physiological mechanisms linking milk yield and fertility in dairy cattle

Fertility in dairy cattle has declined as an unintended consequence of single-trait selection for high milk yield. The unfavorable genetic correlation between milk yield and fertility is now well documented; however, the underlying physiological mechanisms are still uncertain. To understand the relationship between these traits, we developed a method that clusters variants with similar patterns of effects and, after the integration of gene expression data, identifies the genes through which they are likely to act. Biological processes that are enriched in the genes of each cluster were then identified. We identified several clusters with unique patterns of effects. One of the clusters included variants associated with increased milk yield and decreased fertility, where the "archetypal" variant (i.e., the one with the largest effect) was associated with the GC gene, whereas others were associated with TRIM32, LRRK2, and U6-associated snRNA. These genes have been linked to transcription and alternative splicing, suggesting that these processes are likely contributors to the unfavorable relationship between the 2 traits. Another cluster, with archetypal variant near DGAT1 and including variants associated with CDH2, BTRC, SFRP2, ZFHX3, and SLITRK5, appeared to affect milk yield but have little effect on fertility. These genes have been linked to insulin, adipose tissue, and energy metabolism. A third cluster with archetypal variant near ZNF613 and including variants associated with ROBO1, EFNA5, PALLD, GPC6, and PTPRT were associated with fertility but not milk yield. These genes have been linked to GnRH neuronal migration, embryonic development, or ovarian function. The use of archetypal clustering to group variants with similar patterns of effects may assist in identifying the biological processes underlying correlated traits. The method is hypothesis generating and requires experimental confirmation. However, we have uncovered several novel mechanisms potentially affecting milk production and fertility such as GnRH neuronal migration. We anticipate our method to be a starting point for experimental research into novel pathways, which have been previously unexplored within the context of dairy production.

Temporal dynamics of genetic parameters and SNP effects for performance and disorder traits in poultry undergoing genomic selection

Accurate genetic parameters are crucial for predicting breeding values and selection responses in breeding programs. Genetic parameters change with selection, reducing additive genetic variance and changing genetic correlations. This study investigates the dynamic changes in genetic parameters for residual feed intake (RFI), gain (GAIN), breast percentage (BP), and femoral head necrosis (FHN) in a broiler population that undergoes selection, both with and without the use of genomic information. Changes in single nucleotide polymorphism (SNP) effects were also investigated when including genomic information. The dataset containing 200,093 phenotypes for RFI, 42,895 for BP, 203,060 for GAIN, and 63,349 for FHN was obtained from 55 mating groups. The pedigree included 1,252,619 purebred broilers, of which 154,318 were genotyped with a 60K Illumina Chicken SNP BeadChip. A Bayesian approach within the GIBBSF90 + software was applied to estimate the genetic parameters for single-, two-, and four-trait models with sliding time intervals. For all models, we used genomic-based (GEN) and pedigree-based approaches (PED), meaning with or without genotypes. For GEN (PED), heritability varied from 0.19 to 0.2 (0.31 to 0.21) for RFI, 0.18 to 0.11 (0.25 to 0.14) for GAIN, 0.45 to 0.38 (0.61 to 0.47) for BP, and 0.35 to 0.24 (0.53 to 0.28) for FHN, across the intervals. Changes in genetic correlations estimated by GEN (PED) were 0.32 to 0.33 (0.12 to 0.25) for RFI-GAIN, -0.04 to -0.27 (-0.18 to -0.27) for RFI-BP, -0.04 to -0.07 (-0.02 to -0.08) for RFI-FHN, -0.04 to 0.04 (0.06 to 0.2) for GAIN-BP, -0.17 to -0.06 (-0.02 to -0.01) for GAIN-FHN, and 0.02 to 0.07 (0.06 to 0.07) for BP-FHN. Heritabilities tended to decrease over time while genetic correlations showed both increases and decreases depending on the traits. Similar to heritabilities, correlations between SNP effects declined from 0.78 to 0.2 for RFI, 0.8 to 0.2 for GAIN, 0.73 to 0.16 for BP, and 0.71 to 0.14 for FHN over the eight intervals with genomic information, suggesting potential epistatic interactions affecting genetic trait architecture. Given rapid genetic architecture changes and differing estimates between genomic and pedigree-based approaches, using more recent data and genomic information to estimate variance components is recommended for populations undergoing genomic selection to avoid potential biases in genetic parameters.

Identification of genetic associations and functional SNPs of bovine KLF6 gene on milk production traits in Chinese holstein

BACKGROUND

Our previous research identified the Kruppel like factor 6 (KLF6) gene as a prospective candidate for milk production traits in dairy cattle. The expression of KLF6 in the livers of Holstein cows during the peak of lactation was significantly higher than that during the dry and early lactation periods. Notably, it plays an essential role in activating peroxisome proliferator-activated receptor α (PPARα) signaling pathways. The primary aim of this study was to further substantiate whether the KLF6 gene has significant genetic effects on milk traits in dairy cattle.

RESULTS

Through direct sequencing of PCR products with pooled DNA, we totally identified 12 single nucleotide polymorphisms (SNPs) within the KLF6 gene. The set of SNPs encompasses 7 located in 5' flanking region, 2 located in exon 2 and 3 located in 3' untranslated region (UTR). Of these, the g.44601035G > A is a missense mutation that resulting in the replacement of arginine (CGG) with glutamine (CAG), consequently leading to alterations in the secondary structure of the KLF6 protein, as predicted by SOPMA. The remaining 7 regulatory SNPs significantly impacted the transcriptional activity of KLF6 following mutation (P < 0.005), manifesting as changes in transcription factor binding sites. Additionally, 4 SNPs located in both the UTR and exons were predicted to influence the secondary structure of KLF6 mRNA using the RNAfold web server. Furthermore, we performed the genotype-phenotype association analysis using SAS 9.2 which found all the 12 SNPs were significantly correlated to milk yield, fat yield, fat percentage, protein yield and protein percentage within both the first and second lactations (P < 0.0001 ~ 0.0441). Also, with Haploview 4.2 software, we found the 12 SNPs linked closely and formed a haplotype block, which was strongly associated with five milk traits (P < 0.0001 ~ 0.0203).

CONCLUSIONS

In summary, our study represented the KLF6 gene has significant impacts on milk yield and composition traits in dairy cattle. Among the identified SNPs, 7 were implicated in modulating milk traits by impacting transcriptional activity, 4 by altering mRNA secondary structure, and 1 by affecting the protein secondary structure of KLF6. These findings provided valuable molecular insights for genomic selection program of dairy cattle.

Genetic Variance Estimation over Time in Broiler Breeding Programmes for Growth and Reproductive Traits

Monitoring the genetic variance of traits is a key priority to ensure the sustainability of breeding programmes in populations under directional selection, since directional selection can decrease genetic variation over time. Studies monitoring changes in genetic variation have typically used long-term data from small experimental populations selected for a handful of traits. Here, we used a large dataset from a commercial breeding line spread over a period of twenty-three years. A total of 2,059,869 records and 2,062,112 animals in the pedigree were used for the estimations of variance components for the traits: body weight (BWT; 2,059,869 records) and hen-housed egg production (HHP; 45,939 records). Data were analysed with three estimation approaches: sliding overlapping windows, under frequentist (restricted maximum likelihood (REML)) and Bayesian (Gibbs sampling) methods; expected variances using coefficients of the full relationship matrix; and a "double trait covariances" analysis by computing correlations and covariances between the same trait in two distinct consecutive windows. The genetic variance showed marginal fluctuations in its estimation over time. Whereas genetic, maternal permanent environmental, and residual variances were similar for BWT in both the REML and Gibbs methods, variance components when using the Gibbs method for HHP were smaller than the variances estimated when using REML. Large data amounts were needed to estimate variance components and detect their changes. For Gibbs (REML), the changes in genetic variance from 1999-2001 to 2020-2022 were 82.29 to 93.75 (82.84 to 93.68) for BWT and 76.68 to 95.67 (98.42 to 109.04) for HHP. Heritability presented a similar pattern as the genetic variance estimation, changing from 0.32 to 0.36 (0.32 to 0.36) for BWT and 0.16 to 0.15 (0.21 to 0.18) for HHP. On the whole, genetic parameters tended slightly to increase over time. The expected variance estimates were lower than the estimates when using overlapping windows. That indicates the low effect of the drift-selection process on the genetic variance, or likely, the presence of genetic variation sources compensating for the loss. Double trait covariance analysis confirmed the maintenance of variances over time, presenting genetic correlations >0.86 for BWT and >0.82 for HHP. Monitoring genetic variance in broiler breeding programmes is important to sustain genetic progress. Although the genetic variances of both traits fluctuated over time, in some windows, particularly between 2003 and 2020, increasing trends were observed, which warrants further research on the impact of other factors, such as novel mutations, operating on the dynamics of genetic variance.

Genetic analysis of longevity and their associations with fertility traits in Holstein cattle

The increase of longevity is intended to reduce involuntary culling rates, not extend the life span, and it reflects the ability of animals to successfully cope with the environment and disease during production. Sire model, animal model and repeatability animal models were used to estimate the (co) variance components of longevity and fertility traits. Six longevity and thirteen fertility traits were analysed, including herd life (HL), productive life (PL), number of days between first calving and the end of first lactation or culling (L1); number of days between first calving and the end of the second lactation or culling (L2); number of days between first calving and the end of the third lactation or culling (L3); number of days between first calving and the end of the fourth lactation or culling (L4); age at first service, age at first calving (AFC), the interval from first to last inseminations in heifer (IFLh), conception rate of first insemination in heifer, days open (DO), calving interval, gestation length, interval from calving to first insemination (ICF), interval from first to last inseminations in cow (IFLc), conception rate of first insemination in cow, calving ease (CE), birth weight, and calf survival. The estimated heritabilities (±SE) were 0.018 (±0.003), 0.015 (±0.003), 0.049 (±0.004), 0.025 (±0.003), 0.009 (±0.002) and 0.011 (±0.002) for HL, PL, L1, L2, L3 and L4, respectively. Strong correlations were appeared in HL and PL; the genetic and phenotypic correlation coefficients were 0.998 and 0.985, respectively. There were high genetic and phenotypic correlations which were observed in L1 and L2, L2 and L3, L3 and L4, respectively. All fertility traits of heifer showed medium to high heritability, while the cow showed low heritability. All heifer fertility traits had low genetic associations with longevity traits, ranging from -0.018 (L2 and IFLh) to 0.257 (L3 and AFC). Most of the fertility traits showed negative correlations with longevity traits in different parities, and we recommend DO, ICF, IFLc and CE as indirect indicators of longevity traits in dairy cows, but we also need to take into account the differences between parities.

Dairy cow longevity: Impact of animal health and farmers' investment decisions

A dairy farmer's decision to cull or keep dairy cows is likely a complex decision based on animal health and farm management practices. The present paper investigated the relationship between cow longevity and animal health, and between longevity and farm investments, while controlling for farm-specific characteristics and animal management practices, by using Swedish dairy farm and production data for the period 2009 to 2018. We used the ordinary least square and unconditional quantile regression model to perform mean-based and heterogeneous-based analysis, respectively. Findings from the study indicate that, on average, animal health has a negative but insignificant effect on dairy herd longevity. This implies that culling is predominantly done for other reasons than poor health status. Investment in farm infrastructure has a positive and significant effect on dairy herd longevity. The investment in farm infrastructure creates room for new or superior recruitment heifers without the need to cull existing dairy cows. Production variables that prolong dairy cow longevity include higher milk yield and an extended calving interval. Findings from this study imply that the relatively short longevity of dairy cows in Sweden compared with some dairy producing countries is not a result of problems with health and welfare. Rather, dairy cow longevity in Sweden hinges on the farmers' investment decisions, farm-specific characteristics and animal management practices.

Trend analysis and prediction of seasonal changes in milk composition from a pasture-based dairy research herd

The composition of seasonal pasture-produced milk is influenced by stage of lactation, animal genetics, and nutrition, which affects milk nutritional profile and processing characteristics. The objective was to study the effect of lactation stage (early, mid, and late lactation) and diet on milk composition in an Irish spring calving dairy research herd from 2012 to 2020 using principal component and predictive analytics. Crude protein, casein, fat, and solids increased from 2012 to 2020, whereas lactose concentration peaked in 2017, then decreased. Based on seasonal data from 2013 to 2016, forecasting models were successfully created to predict milk composition for 2017 to 2020. The diet of cows in this study is dependent upon grass growth rates across the milk production season, which in turn, are influenced by weather patterns, whereby extreme weather conditions (rainfall and temperature) were correlated with decreasing grass growth and increasing nonprotein nitrogen levels in milk. The study demonstrates a significant change in milk composition since 2012 and highlights the effect that seasonal changes such as weather and grass growth have on milk composition of pasture-based systems. The potential to forecast milk composition at different stages of lactation benefits processers by facilitating the optimization of in-process and supply logistics of dairy products.

The correlation of substitution effects across populations and generations in the presence of nonadditive functional gene action

Allele substitution effects at quantitative trait loci (QTL) are part of the basis of quantitative genetics theory and applications such as association analysis and genomic prediction. In the presence of nonadditive functional gene action, substitution effects are not constant across populations. We develop an original approach to model the difference in substitution effects across populations as a first order Taylor series expansion from a "focal" population. This expansion involves the difference in allele frequencies and second-order statistical effects (additive by additive and dominance). The change in allele frequencies is a function of relationships (or genetic distances) across populations. As a result, it is possible to estimate the correlation of substitution effects across two populations using three elements: magnitudes of additive, dominance, and additive by additive variances; relationships (Nei's minimum distances or Fst indexes); and assumed heterozygosities. Similarly, the theory applies as well to distinct generations in a population, in which case the distance across generations is a function of increase of inbreeding. Simulation results confirmed our derivations. Slight biases were observed, depending on the nonadditive mechanism and the reference allele. Our derivations are useful to understand and forecast the possibility of prediction across populations and the similarity of GWAS effects.

Analysis of Longevity Traits in Holstein Cattle: A Review

Dairy cow longevity is an essential economic trait that can supplement the breeding value of production traits, which is related to the herd time and lifetime milk yield of dairy cows. However, longevity is a relatively difficult trait to select for dairy cow breeding due to low heritability and numerous influence factors of the longevity in dairy cows. Longevity trait has been used as an important breeding target of a comprehensive selection index in many dairy developed countries; however, it has not been included in performance index in many developing countries. At present, cows in these countries are still in the primary stage of “large quantity, low quality, high cost, and low yield.” The average parity of dairy cows is less than 2.7, which is difficult to maintain the production efficiency to meet the demands of the dairy industry. Therefore, there is an urgent need to select and breed for the longevity of dairy cows. The various definitions and models (including linear, threshold, random regression, sire, and survival analysis) of longevity were reviewed and standardized. Survival analysis is the optimal model to evaluate longevity, and the longevity heritability is 0.01–0.30 by using different definitions and models. Additionally, the relationship between longevity and other traits was summarized, and found that longevity was regulated by multiple factors, and there were low or medium genetic correlations between them. Conformation traits, milk production traits, reproductive traits, and health traits may be used as indicators to select and breed the longevity of dairy cows. The genetic assessment methods, heritability, influencing factors, importance, breeding, and genetics of longevity were reviewed in the manuscript, which could provide a valuable reference for the selective breeding to extend the productive life of Holstein cattle.

Improving the accuracy of predictions for cow survival by multivariate evaluation model

Context Cow survival measures the ability of cows to survive from the current to subsequent lactation. In addition to economic gain, genetic selection for survival could improve animal welfare by increasing the adaptability and resilience of the cows to both environmental and health challenges. However, survival is a complex trait because it results from a diverse range of reasons for culling of cows from the herd. Consequently, the accuracy of genetic predictions of direct survival are often low. Aims Our aim was to increase the accuracy of predictions of survival in Holstein and Jersey sires by including important predictor traits in multi-trait evaluation models. Methods Phenotypic and genetic correlations between survival trait deviations (TDs) and 35 routinely measured traits (including milk yield, fertility and type traits) were estimated using bivariate sire models. Survival TDs for 538 394 Holstein and 63 839 Jersey cows were used in our study; these cows or their close relatives also had milk, fertility and type traits records between 2002 and 2019. These genetic parameters were required to assess the potential usefulness of predictor traits for the prediction of survival. Key results Survival was genetically correlated with milk, fat and protein yields, overall type, composite mammary system and fertility TDs in both Holstein and Jersey. Further, most of the type traits related to feet and legs, and rump, were also correlated with survival TDs in Jersey. For sires, the accuracy of predictions for survival increased by 0.05 for Holsteins (from 0.54 to 0.59) and for Jerseys (from 0.48 to 0.53) through the use of multivariate models compared with univariate models. Conclusions Survival was genetically associated with traits affecting voluntary and involuntary culling and when included in multi-trait genetic evaluation models, they moderately improved the accuracy of genetic prediction of survival. Implications Predictor traits can be used to increase the accuracy of predictions of survival through the use of multi-trait models. The inclusion of breed-specific predictor traits should be considered, especially for Jerseys in genetic evaluations of survival.

Impact of a multiple-test strategy on breeding index development for the Australian dairy industry

Context A high level of acceptance and use of breeding indices by farmers and breeding companies that target a National breeding objective is an effective strategy to achieve high rates of genetic gain. Indices require maintenance to ensure that they reflect current economic and genetic trends and farmer preferences. Often, indices are tested on an average herd on the basis of, for example, milk composition and calving pattern. However, this strategy does not differentiate the impact on breeds. Australian dairy farmers routinely make breeding decisions by using the balanced performance index (BPI) or the health weighted index, published by DataGene. Aims The aim of the present study was to test new selection indices on the most popular breeds to better understand the genetic progress that each breed is expected to make. Existing economic models were updated to reflect changing trends in input costs and milk income. Consultative processes identified opportunities to improve alignment between farmer preferences and Australia’s National Breeding Objective. In response, more than 20 selection index options were developed and options were discussed with industry. Methods Indices were evaluated on three breeds in the following three ways: (1) expected response to selection from the use of each index, (2) index and trait correlations, and (3) relative trait emphasis. Key results Farmer trait preferences varied by breed and this information was considered in the development of economic weights. The updated BPI has primary emphasis on production traits (44% in Holstein, 49% in Reds), secondary emphasis on health and fertility (35% in Holstein, 29% in Reds), tertiary emphasis on type, workability and feed saved. The equivalent index for Jerseys is similar, but following stakeholder feedback to multiple tests, it was decided to remove emphasis on the feed saved estimated breeding values, so that the percentage emphasis on trait groups in Jerseys is 51% production, 32% health and fertility and the remainder on type and workability. Implications Understanding trait preferences and testing indices on different breeds can change the decisions that are made during index development. Conclusions Developing a better understanding of the differences among breeds had a positive impact on farmer engagement and resulted in a modified BPI for the Jersey breed.

Analysis of culling reasons and age at culling in Australian dairy cattle

Context A thorough analysis of the reasons for culling was made to understand the phenotypic trend in herd life. In addition, identification of culling reasons could enable to develop a strategy for further evaluation of longevity in Australian dairy cows. Aims The aim of this study was to investigate the main causes of culling in Australian dairy herds and thereby to assess the trend of reason-specific culling over time. Methods Culling reasons in Australian dairy cattle were studied based on culling records from 1995 through 2016. A total of 2452124 individual cow culling observations were obtained from Datagene, Australia, of which 2140337 were Holstein and 311787 were from Jersey cows. A binary logistic regression model was used to estimate effects of breed and age and the trend of a particular culling reason over time. Key results The most important culling reasons identified over the 21-year period were infertility (17.0%), mastitis (12.9%), low production (9.3%), sold for dairy purpose (6.4%) and old age (6.2%), whereas 37.4% were ‘other reasons not reported’. The average age at culling was nearly the same for Holstein (6.75 years) and Jersey (6.73 years) cows. The estimated age at culling was slightly increased for Holstein cows (by 3.7 days) and somewhat decreased for Jersey cows (by 11 days) over the last two decades. The probability of culling cows for infertility and low production was high in early parities and consistently declined as age advanced, and culling due to mastitis was higher in older cows. The trend of main culling reasons over time was evaluated, indicating that the probability of culling due to infertility has progressively increased over the years in both breeds, and culling for mastitis in Jersey cows has also increased. Culling of cows due to low production sharply decreased from 2.5 to –8% for Holstein and from 73 to 60% for Jersey cows over the 21-year period. Conclusions Culling age has changed only little in both breeds whereas culling reasons have changed over the last two decades, with low production becoming a less important reason for culling and infertility becoming more important for Holstein and Jersey breeds. Implications Due to changes of culling reasons, there could be a change in the meaning of survival over time as well. As a result, genetic correlation with survival and other traits might be changed and accuracy and bias of genetic evaluations could be affected.

Beef trait genetic parameters based on old and recent data and its implications for genomic predictions in Italian Simmental cattle

This study aimed to evaluate the changes in variance components over time to identify a subset of data from the Italian Simmental (IS) population that would yield the most appropriate estimates of genetic parameters and breeding values for beef traits to select young bulls. Data from bulls raised between 1986 and 2017 were used to estimate genetic parameters and breeding values for four beef traits (average daily gain [ADG], body size [BS], muscularity [MUS], and feet and legs [FL]). The phenotypic mean increased during the years of the study for ADG, but it decreased for BS, MUS, and FL. The complete dataset (ALL) was divided into four generational subsets (Gen1, Gen2, Gen3, and Gen4). Additionally, ALL was divided into two larger subsets: the first one (OLD) combined data from Gen1 and Gen2 to represent the starting population, and the second one (CUR) combined data from Gen3 and Gen4 to represent a subpopulation with stronger ties to the current population. Genetic parameters were estimated with a four-trait genomic animal model using a single-step genomic average information restricted maximum likelihood algorithm. Heritability estimates from ALL were 0.26 ± 0.03 for ADG, 0.33 ± 0.04 for BS, 0.55 ± 0.03 for MUS, and 0.23 ± 0.03 for FL. Higher heritability estimates were obtained with OLD and ALL than with CUR. Considerable changes in heritability existed between Gen1 and Gen4 due to fluctuations in both additive genetic and residual variances. Genetic correlations also changed over time, with some values moving from positive to negative or even to zero. Genetic correlations from OLD were stronger than those from CUR. Changes in genetic parameters over time indicated that they should be updated regularly to avoid biases in genomic estimated breeding values (GEBV) and low selection accuracies. GEBV estimated using CUR variance components were less biased and more consistent than those estimated with OLD and ALL variance components. Validation results indicated that data from recent generations produced genetic parameters that more appropriately represent the structure of the current population, yielding accurate GEBV to select young animals and increasing the likelihood of higher genetic gains.

Using Random Regression Models to Genetically Evaluate Functional Longevity Traits in North American Angus Cattle

Simple Summary

Cattle longevity is usually defined as the duration of life of a cow from first calving to death. In addition to a longer lifespan, it is crucial that cows are productive throughout their lives. Incorporating optimal indicators of productive longevity in breeding schemes will directly improve the economic profitability of the beef cattle herd and long-term sustainability of the industry. Thus, the impact of different longevity indicators in the selection of North American Angus cattle was evaluated and optimal parameters were defined to perform the evaluations.

Abstract

This study aimed to propose novel longevity indicators by comparing genetic parameters for traditional (TL; i.e., the cow’s lifespan after the first calving) and functional (FL; i.e., how long the cow stayed in the herd while also calving; assuming no missing (FLa) or missing (FLb) records for unknown calving) longevity, considering different culling reasons (natural death, structural problems, disease, fertility, performance, and miscellaneous). Longevity definitions were evaluated from 2 to 15 years of age, using single- and multiple-trait Bayesian random regression models (RRM). The RRM fitting heterogenous residual variance and fourth order Legendre polynomials were considered as the optimal models for the majority of longevity indicators. The average heritability estimates over ages for FLb (from 0.08 to 0.25) were always higher than those for FLa (from 0.07 to 0.19), and higher or equal to the ones estimated for TL (from 0.07 to 0.23), considering the different culling reasons. The average genetic correlations estimated between ages were low to moderate (~0.40), for all longevity definitions and culling reasons. However, removing the extreme ages (i.e., 2 and >12 years) increased the average correlation between ages (from ~0.40 to >0.70). The genetic correlations estimated between culling reasons were low (0.12 and 0.20 on average, considering all ages and ages between 3 and 12 years old, respectively), indicating that longevity based on different culling reasons should be considered as different traits in the genetic evaluations. Higher average genetic correlations (estimated from 3 to 12 years old) were observed between TL and FLb (0.73) in comparison to TL and FLa (0.64), or FLa and FLb (0.65). Consequently, a higher average proportion of commonly-selected sires, for the top 1% sires, was also observed between TL and FLb (91.74%), compared to TL and FLa (59.68%), or FLa and FLb (61.01%). Higher prediction accuracies for the expected daughter performances (calculated based on the pedigree information) were obtained for FLb in comparison to TL and FLa. Our findings indicate that FLb is preferred for the genetic evaluation of longevity. In addition, it is recommended including multiple longevity traits based on different groups of culling reasons in a selection sub-index, as they are genetically-different traits. Genetic selection based on breeding values at the age of four years is expected to result in greater selection responses for increased longevity in North American Angus cattle.

Application of multivariate single-step SNP best linear unbiased predictor model and revised SNP list for genomic evaluation of dairy cattle in Australia

The objectives of this study were (1) to evaluate the computational feasibility of the multitrait test-day single-step SNP-BLUP (ssSNP-BLUP) model using phenotypic records of genotyped and nongenotyped animals, and (2) to compare accuracies (coefficient of determination; R²) and bias of genomic estimated breeding values (GEBV) and de-regressed proofs as response variables in 3 Australian dairy cattle breeds (i.e., Holstein, Jersey, and Red breeds). Additive genomic random regression coefficients for milk, fat, protein yield and somatic cell score were predicted in the first, second, and third lactation. The predicted coefficients were used to derive 305-d GEBV and were compared with the traditional parent averages obtained from a BLUP model without genomic information. Cow fertility traits were evaluated from the 5-trait repeatability model (i.e., calving interval, days from calving to first service, pregnancy diagnosis, first service nonreturn rate, and lactation length). The de-regressed proofs were only for calving interval. Our results showed that ssSNP-BLUP using multitrait test-day model increased reliability and reduced bias of breeding values of young animals when compared with parent average from traditional BLUP in Australian Holsten, Jersey, and Red breeds. The use of a custom selection of approximately 46,000 SNP (custom XT SNP list) increased the reliability of GEBV compared with the results obtained using the commercial Illumina 50K chip (Illumina, San Diego, CA). The use of the second preconditioner substantially improved the convergence rate of the preconditioned conjugate gradient method, but further work is needed to improve the efficiency of the computation of the Kronecker matrix product by vector. Application of ssSNP-BLUP to multitrait random regression models is computationally feasible.

Changes in genetic parameters for fitness and growth traits in pigs under genomic selection

Genomic selection increases accuracy and decreases generation interval, speeding up genetic changes in the populations. However, intensive changes caused by selection can reduce the genetic variation and can strengthen undesirable genetic correlations. The purpose of this study was to investigate changes in genetic parameters for fitness traits related with prolificacy (FT1) and litter survival (FT2 and FT3), and for growth (GT1 and GT2) traits in pigs over time. The data set contained 21,269 (FT1), 23,246 (FT2), 23,246 (FT3), 150,492 (GT1), and 150,493 (GT2) phenotypic records obtained from 2009 to 2018. The pedigree file included 369,776 animals born between 2001 and 2018, of which 39,103 were genotyped. Genetic parameters were estimated with bivariate models (FT1-GT1, FT1-GT2, FT2-GT1, FT2-GT2, FT3-GT1, and FT3-GT2) using 3-yr sliding subsets. With a Bayesian implementation using the GIBBS3F90 program computations were performed as genomic analysis (GEN) or pedigree-based analysis (PED), that is, with or without genotypes, respectively. For GEN (PED), the changes in heritability from the first to the last year interval, that is, from 2009–2011 to 2015–2018 were 8.6 to 5.6 (7.9 to 8.8) for FT1, 7.8 to 7.2 (7.7 to 10.8) for FT2, 11.4 to 7.6 (10.1 to 7.5) for FT3, 35.1 to 16.5 (32.5 to 23.7) for GT1, and 35.9 to 16.5 (32.6 to 24.1) for GT2. Differences were also observed for genetic correlations as they changed from −0.31 to −0.58 (−0.28 to −0.73) for FT1-GT1, −0.32 to −0.50 (−0.29 to −0.74) for FT1-GT2, −0.27 to −0.45 (−0.30 to −0.65) for FT2-GT1, −0.28 to −0.45 (−0.32 to −0.66) for FT2-GT2, 0.14 to 0.17 (0.11 to 0.04) for FT3-GT1, and 0.14 to 0.18 (0.11 to 0.05) for FT3-GT2. Strong selection in pigs reduced heritabilities and emphasized the antagonistic genetic relationships between fitness and growth traits. With genotypes considered, heritability estimates were smaller and genetic correlations were greater than estimates with only pedigree and phenotypes. When selection is based on genomic information, genetic parameters estimated without this information can be biased because preselection is not accounted for by the model.

Patterns of environmental variance across environments and traits in domestic cattle

The variance in phenotypic trait values is a product of environmental and genetic variation. The sensitivity of traits to environmental variation has a genetic component and is likely to be under selection. However, there are few studies investigating the evolution of this sensitivity, in part due to the challenges of estimating the environmental variance. The livestock literature provides a wealth of studies that accurately partition components of phenotypic variance, including the environmental variance, in well-defined environments. These studies involve breeds that have been under strong selection on mean phenotype in optimal environments for many generations, and therefore represent an opportunity to study the potential evolution of trait sensitivity to environmental conditions. Here, we use literature on domestic cattle to examine the evolution of micro-environmental variance (CVR-the coefficient of residual variance) by testing for differences in expression of CVR in animals from the same breed reared in different environments. Traits that have been under strong selection did not follow a null expectation of an increase in CVR in heterogenous environments (e.g., grazing), a pattern that may reflect evolution of increased uniformity in heterogeneous environments. When comparing CVR across environments of different levels of optimality, here measured by trait mean, we found a reduction in CVR in the more optimal environments for both life history and growth traits. Selection aimed at increasing trait means in livestock breeds typically occurs in the more optimal environments, and we therefore suspect that the decreased CVR is a consequence of evolution of the expression of micro-environmental variance in this environment. Our results highlight the heterogeneity in micro-environmental variance across environments and point to possible connections to the intensity of selection on trait means.

Prediction of random-regression coefficient for daily milk yield after 305 days in milk by using the regression-coefficient estimates from the first 305 days

Objective

Because lactation periods in dairy cows lengthen with increasing total milk production, it is important to predict individual productivities after 305 days in milk (DIM) to determine the optimal lactation period. We therefore examined whether the random regression (RR) coefficient from 306 to 450 DIM (M2) can be predicted from those during the first 305 DIM (M1) by using a RR model.

Methods

We analyzed test-day milk records from 85,690 Holstein cows in their first lactations and 131,727 cows in their later (second to fifth) lactations. Data in M1 and M2 were analyzed separately by using different single-trait RR animal models. We then performed a multiple regression analysis of the RR coefficients of M2 on those of M1 during the first and later lactations.

Results

The first-order Legendre polynomials were practical covariates of RR for the milk yields of M2. All RR coefficients for the additive genetic (AG) effect and the intercept for the permanent environmental (PE) effect of M2 had moderate to strong correlations with the intercept for the AG effect of M1. The coefficients of determination for multiple regression of the combined intercepts for the AG and PE effects of M2 on the coefficients for the AG effect of M1 were moderate to high. The daily milk yields of M2 predicted by using the RR coefficients for the AG effect of M1 were highly correlated with those obtained by using the coefficients of M2.

Conclusion

Milk production after 305 DIM can be predicted by using the RR coefficient estimates of the AG effect during the first 305 DIM.