Milk protein composition in purebred Holsteins and in first/second-generation crossbred cows from Swedish Red,Montbeliarde and Brown Swiss bulls

The aim of this study was to analyze milk protein composition in purebred and crossbred dairy cattle and estimate the effects of individual sources of variation on the investigated traits. Milk samples were collected from 505 cows from three commercial farms located in Northern Italy, some of which had originated from crossbreeding programs, although most were purebred Holsteins (HO). The basic crossbreeding scheme was a three-breed rotational system using Swedish Red (SR) semen on HO cows (SR×HO), Montbeliarde (MO) semen on SR×HO cows (MO×(SR×HO)) and HO semen again on MO×(SR×HO) cows. A smaller number of purebred HO from each of the herds were mated inverting the breed order (MO×HO and SR×(MO×HO)) or using Brown Swiss (BS) bulls (BS×HO) then MO bulls (MO×(BS×HO)). Milk samples were analyzed by reverse-phase HPLC to obtain protein fraction amounts (g/l) and proportions (% of total true protein). Traits were analyzed using a linear model, which included the fixed effects of herd-test-day (HTD), parity, days in milk and breed combination. Results showed that milk protein fractions were influenced by HTD, stage of lactation, parity and breed combination. The increase in protein concentration during lactation was due in particular to β-casein (β-CN), α_S1-CN and β-lactoglobulin (β-LG). The higher protein content of primiparous milk was mainly due to higher concentrations of all casein fractions. The milk from crossbred cows had higher contents and proportions of κ-CN and α-lactalbumin (α-LA), lower proportions of β-LG and greater proportion of caseins/smaller in whey proteins on milk true protein than purebred HO. The three-way crossbreds differed from two-way crossbreds only in having greater proportions of α-LA in their milk. Of the three-way crossbreds, the SR sired cows yielded milk with a smaller content and proportion of β-LG than the MO sired cows, and, consequently, a higher proportion of caseins than whey proteins. Results from this study support the feasibility of using crossbreeding programs to alter milk protein profiles with the aim of improving milk quality and cheese-making properties.     

The effects of crossbreeding with Norwegian Red dairy cattle on common postpartum diseases,fertility and body condition score

Norwegian Red bulls, selected in Norway, have been used for crossbreeding with Israeli Holstein on commercial farms. The aim of this project was to investigate Norwegian Red×Israeli Holstein (NRX) performance to see how the daughters perform in a different environment than the one their sires were selected in. This was done by comparing health and fertility of NRX with their Israeli Holstein (HO) counterparts. The data consisted of 71 911 HO records and 10 595 NRX records from 33 855 cows in 23 Israeli dairy herds. Calving events took place between 2006 and 2016. Five postpartum disorders (mean frequency in HO v. NRX, %) recorded by veterinarians were analyzed: anestrus (37.4 v. 41.2), metritis (40.1 v. 28.6), ketosis (11.9 v. 7.1), lameness (7.1 v. 3.1) and retained placenta (6.2 v. 4.0). The incidence of abortions was also analyzed; HO had a mean frequency of 9.9% and NRX 8.2%. These traits were defined as binary traits, with ‘1’ indicating that the disorder was present and a treatment took place at least once, or ‘0’ if the cow did not show signs of that disorder. Days open (i.e. the number of days from calving to conception), body condition score (BCS) recorded on a 1 to 5 scale and changes in BCS from calving to peak lactation were also analyzed. A logistic model was used for the health traits, while days open and BCS were analyzed with linear models. The model included breed group, herd-year of calving, birth year and parity as fixed effects. There was a significantly higher risk (odds ratio for HO v. NRX in parentheses) of ketosis (1.46), metritis (1.78), lameness (2.07), retained placenta (1.41) and abortion (1.13) in HO compared with NRX. Israeli Holstein heifers and cows in parity 3 to 6 had fewer cases of anestrus than NRX but no differences were found between the groups in parities 1 and 2. Body condition score was higher for NRX than HO and there was less change in BCS from calving to peak lactation in NRX compared with HO. Likewise, NRX had fewer days open than HO. Results indicate that crossbreeding can produce cows with better fertility that are less susceptible to postpartum disorders.     

Quantitative trait loci for clinical mastitis on chromosomes 2, 6, 14 and 20 in Norwegian Red cattle

Mastitis is the most frequent and costly disease in dairy production and solutions leading to a reduction in the incidence of mastitis are highly demanded. Here a genome-wide association study was performed to identify polymorphisms affecting susceptibility to mastitis. Genotypes for 17 349 SNPs distributed across the 29 bovine autosomal chromosomes from a total of 2589 sires with 1 389 776 daughters with records on clinical mastitis were included in the analysis. Records of occurrence of clinical mastitis were divided into seven time periods in the first three lactations in order to identify quantitative trait loci affecting mastitis susceptibility in particular phases of lactation. The most convincing results from the association mapping were followed up and validated by a combined linkage disequilibrium and linkage analysis. The study revealed quantitative trait loci affecting occurrence of clinical mastitis in the periparturient period on chromosomes 2, 6 and 20 and a quantitative trait locus affecting occurrence of clinical mastitis in late lactation on chromosome 14. None of the quantitative trait loci for clinical mastitis detected in the study seemed to affect lactation average of somatic cell score. The SNPs highly associated with clinical mastitis lie near both the gene encoding interleukin 8 on chromosome 6 and the genes encoding the two interleukin 8 receptors on chromosome 2.     

Smallholder experiences with dairy cattle crossbreeding in the tropics: from introduction to impact

Crossbreeding of indigenous tropical and improved western dairy cattle breeds as tool to improve dairy cattle performance on smallholder farms has been widely advocated, criticised and yet applied. The government of Ethiopia supported this technology for decades but adoption rate is low. Constraints are documented but there is little information about farm level introduction and development of crossbreeding. A total 122 smallholders with mixed crop livestock farms and at least 8 years of successful crossbreeding were interviewed using a pre-tested questionnaire in two contexts in Amhara Regional state in north-western Ethiopia. Crossbreeding initiator was either uncoordinated government extension or a coordinated development project, also implemented with governmental support. Qualitative and quantitative data on farmers’ motivations, crossbreeding introduction, initiator support, breeding adaptation and impacts at farm level were analysed. Results show that even though motives vary between contexts the underlying reason to introduce crossbreeding was economic profit. To be able to introduce crossbreeding support of initiators (e.g. extension) and other farmers was essential. The crossbreeding introduction context had some influence. Governmental actors were the main source of support and supplier of exotic genetics but the farmer network acted as safety net filling gaps of government support. Breeding strategies focused on performance increase. A lack of basic understanding of crossbreeding has been identified. A surprising, probably biased, result was general satisfaction with initiator support and with breeding services. It was challenged by the high proportion of farmers unable to follow a breeding strategy due to insufficient bull and/or semen supply. Crossbreeding changed the smallholder production system to a high input – high output system. Except for crossbred adaptation problems, challenges were ranked context specific and influenced by the initiator. Farmers perceived crossbreeding as success and recommended it. We conclude that farmers can realize income increase with crossbreeding. The complexity of this technology, high initial investment and the need for support services and external production inputs are probable reasons why crossbreeding uptake is low. Improving the availability of semen and/or bulls must be the top priority for breeding service providers to enable farmers to follow a breeding strategy and reach a suitable and sustainable herd performance. Access to investment capital, input supply, strong technical support and market linkages are crucial for successful crossbreeding.     

Beef production potential of Norwegian Red and Holstein-Friesian bulls slaughtered at two ages

There is a paucity of data on the beef production potential of Norwegian Red (NOR) compared with 'modern' Holstein-Friesian (HF) cattle. The present study used a total of 64 bulls in a 2 × 2 factorial design study encompassing two breeds (HF and NOR) and two slaughter ages (485; E, and 610; L, days). The mean initial age and live weight of the HF bulls were 179 (s.d. 47.1) days and 203 (s.d. 64.0) kg, while the corresponding data for the NOR bulls were 176 (s.d. 39.7) days and 185 (s.d. 63.6) kg, respectively. Bulls were offered a 50 : 50 mixture (dry matter (DM) basis) of grass silage and concentrates. No breed × slaughter group interactions were recorded for any parameters evaluated (P > 0.05). HF bulls had higher (P < 0.001) DM intake and poorer (P < 0.01) efficiency of conversion of food to carcass gain than NOR bulls. HF bulls tended (P = 0.07) to have a higher rate of live-weight gain and were heavier (P < 0.001) at slaughter than NOR bulls, though both carcass weight and rate of carcass gain did not differ between the breeds (P > 0.05). NOR bulls had higher (P < 0.001) dressing proportion and carcass conformation score than HF bulls, while breed of bull had no influence (P > 0.05) on carcass fat classification, depth of subcutaneous fat, marbling score or on the weight of fat in the internal depots. Daily food intakes did not differ (P > 0.05) across the two slaughter age groups, though efficiency of conversion of food to carcass gain was poorer (P < 0.05) in the L compared with E bulls. Rate of live-weight gain was lower (P < 0.01) for L bulls, although rate of carcass gain did not differ (P > 0.05) between the E and L bulls. Increasing age at slaughter increased (P < 0.01 or greater) dressing proportion, carcass fat class, depth of subcutaneous fat, marbling score and internal fat depots, but had no effect (P > 0.05) on the carcass conformation score. Instrumental measures of meat quality indicated that meat from NOR bulls was tougher (P < 0.01) than meat from HF bulls, while delaying slaughter increased (P < 0.001) a* and C*ab, and decreased (P < 0.01) h0, indicating improved redness. It is concluded that NOR bulls have higher food efficiency and produce more highly conformed carcasses than HF bulls, but HF bulls produce more tender meat.     

Managing the transition from purebred to rotational crossbred dairy cattle herds: three technical pathways from a retrospective case-study analysis

The growing interest in rotational crossbreeding in Western countries is due to its potential to improve reproductive and health performances of cows. Although a large amount of research focuses on assessing crossbred cows’ performances, how to manage the transition from purebred to rotational crossbred herds is under-explored. Based on a retrospective analysis of French dairy herd case studies, we aimed to identify and characterise technical pathways to make such a transition. In 2018, we performed semi-directive interviews on 26 commercial dairy farms. Data were collected to describe changes in breeding, replacement and culling management practices from the first crossbred mating with purebred cows to the management of a mainly crossbred herd in 2018. Based on a multivariate analysis, we identified two main guidelines structuring technical pathways to move towards rotational crossbred herds: (i) the depth and scale of change (i.e. farm v. herd) associated with the introduction of rotational crossbreeding in the whole-farm dynamics and (ii) the changes in herd replacement and breeding practices to adapt to the evolution of herd demographics induced by the evolution of the dairy crossbred mating rate over time (high from the beginning v. distributed over time). Hierarchical clustering discriminated three groups of farmers differing in their technical pathway to move towards a rotational crossbred herd. In pathway 1, farmers customised one or several rotational crossbreeding schemes to support whole-farm transition towards an organic or grass-based system. Once the scheme stabilised, they quickly implemented it and had to readjust replacement and culling practices to regulate imbalance in herd demographics induced by the improvement in cow fertility. In pathway 2, farmers also customised one or several rotational crossbreeding schemes to support whole-farm redesign but they implemented it more gradually in the herd, which induced no major imbalance in herd demographics. In pathway 3, farmers predefined a relatively well-known rotational crossbreeding scheme to correct fertility issues of purebred cows without any changes at the farm level. They implemented it quickly from the beginning and had to adapt herd replacement and culling to regulate imbalance in herd demographics induced by the improvement in cow fertility. These first empirical evidences on how dairy farmers manage the transition from a purebred to rotational crossbred herd provide original scientific and operational contributions.     

A study on associaton between mastitis and serologically defined class I bovine lymphocyte antigens (BOLA-A) in Norwegian cows

A total of 102 cows was tested for class I antigens of the bovine major histocompatibility complex. Half of the animals (51) had completed at least four lactations without any veterinary treatment for mastitis. The distribution of BoLA-A antigens among these relative mastitis-resistant cows was compared to that in the other half of the material (51), which comprised animals with at least one recorded treatment for mastitis. There were no statistically significant differences in BoLA-A antigen frequency between cows with mastitis and cows without mastitis. The two most common antigens were A2 and w16. The frequency of these two antigens deviated from earlier estimates within the Norwegian cattle (NRF) population, the difference for w16 being statistically significant.     

Comparison of artificial insemination and natural service cost effectiveness in dairy cattle

Reproductive efficiency in the dairy herd is the most important factor for its economic success and a major concern for dairy farmers when using artificial insemination (AI) or natural service (NS). Our objectives were to estimate, compare and analyse the costs associated with breeding cattle by do-it-yourself (DIY) AI and NS and identify the factors that influence them, under typical dairy farming conditions in Greece. A simulation study was designed based on data from 120 dairy cattle farms that differed in size (range 40 to 285 cows) and milk production level (4000 to 9300 kg per cow per year). Different scenarios were employed to estimate costs associated directly with AI and NS as well as potentially extended calving intervals (ECI) due to AI. Results showed that bull maintenance costs for NS were €1440 to €1670 per year ($1,820 to $2,111). Direct AI costs were higher than those for NS for farms with more than 30 cows and ECI constituted a considerable additional burden. In fact, amongst the factors that affected the amount of milk needed to cover total extra AI costs, number of days open was the dominant one. Semen, feed and heifer prices had a very small effect. When, hypothetically, use of NS bulls results in a calving interval of 12 months, AI daughters with a calving interval of 13.5 months have to produce about 705 kg of additional milk in order to cover the extra cost. Their actual milk production, however, exceeds this limit by more than 25%. When real calving intervals are considered (13.0 v. 13.7 months for NS and AI, respectively) AI daughters turn out to produce more than twice the additional amount of milk needed. It was concluded that even under less than average management conditions, AI is more profitable than the best NS scenario. The efficient communication of this message should be a primary concern of the AI industry.     

Bayesian single-step genomic evaluations combining local and foreign information in Walloon Holsteins

Most dairy cattle populations found in different countries around the world are small to medium sized and use many artificial insemination bulls imported from different foreign countries. The Walloon population in the southern part of Belgium is a good example for such a small-scale population. Wallonia has also a very active community of Holstein breeders requesting high level genetic evaluation services. Single-step Genomic BLUP (ssGBLUP) methods allow the simultaneous use of genomic, pedigree and phenotypic information and could reduce potential biases in the estimation of genomically enhanced breeding values (GEBV). Therefore, in the context of implementing a Walloon genomic evaluation system for Holsteins, it was considered as the best option. However, in contrast to multi-step genomic predictions, natively ssGBLUP will only use local phenotypic information and is unable to use directly important other sources of information coming from abroad, for example Multiple Across Country Evaluation (MACE) results as provided by the Interbull Center (Uppsala, Sweden). Therefore, we developed and implemented single-step Genomic Bayesian Prediction (ssGBayes), as an alternative method for the Walloon genomic evaluations. The ssGBayes method approximated the correct system of equations directly using estimated breeding values (EBV) and associated reliabilities (REL) without any explicit deregression step. In the Walloon genomic evaluation, local information refers to Walloon EBV and REL and foreign information refers to MACE EBV and associated REL. Combining simultaneously all available genotypes, pedigree, local and foreign information in an evaluation can be achieved but adding contributions to left-hand and right-hand sides subtracting double-counted contributions. Correct propagation of external information avoiding double counting of contributions due to relationships and due to records can be achieved. This ssGBayes method computed more accurate predictions for all types of animals. For example, for genotyped animals with low Walloon REL (<0.25) without MACE results but sired by genotyped bulls with MACE results, the average increase of REL for the studied traits was 0.38 points of which 0.08 points could be traced to the inclusion of MACE information. For other categories of genotyped animals, the contribution by MACE information was also high. The Walloon genomic evaluation system passed for the first time the Interbull GEBV tests for several traits in July 2013. Recent experiences reported here refer to its use in April 2016 for the routine genomic evaluations of milk production, udder health and type traits. Results showed that the proposed methodology should also be of interest for other, similar, populations.     

High levels of linkage disequilibria between serologically defined class I bovine lymphocyte antigens (BOLA-A) and class II DQB restriction fragment length polymorphism (RFLP) in Norwegian cows

Serum defined BoLA-A antigens, together with BoLA-DQB RFLP patterns, were determined in 87 almost unrelated Norwegian cattle. Statistical analysis revealed strong linkage disequilibria between these loci at the population level. A total of 13 haplotypes were found to be present at frequencies significantly greater than those predicted on the basis of their component gene frequencies. Among these, the subgroups 1A and 1B of the DQ1 haplotype were found to be closely associated with the class I antigens A11 and w16, respectively. The association between A11 and DQ1A is of particular interest, as two independent studies, one employing class I serology, and the other RFLP analysis of the class II locus DQ, have previously indicated that A11 and DQ1A confer relative susceptibility to mastitis.     

A genome map of divergent artificial selection between Bos taurus dairy cattle and Bos taurus beef cattle

A number of cattle breeds have become highly specialized for milk or beef production, following strong artificial selection for these traits. In this paper, we compare allele frequencies from 9323 single nucleotide polymorphism (SNP) markers genotyped in dairy and beef cattle breeds averaged in sliding windows across the genome, with the aim of identifying divergently selected regions of the genome between the production types. The value of the method for identifying selection signatures was validated by four sources of evidence. First, differences in allele frequencies between dairy and beef cattle at individual SNPs were correlated with the effects of those SNPs on production traits. Secondly, large differences in allele frequencies generally occurred in the same location for two independent data sets (correlation 0.45) between sliding window averages. Thirdly, the largest differences in sliding window average difference in allele frequencies were found on chromosome 20 in the region of the growth hormone receptor gene, which carries a mutation known to have an effect on milk production traits in a number of dairy populations. Finally, for the chromosome tested, the location of selection signatures between dairy and beef cattle was correlated with the location of selection signatures within dairy cattle.     

Review: How to improve genomic predictions in small dairy cattle populations

This paper reviews strategies and methods to improve accuracies of genomic predictions from the perspective of a numerically small population. Improvements are realized by influencing one or both of the main factors: (1) improve or increase genomic connections to phenotypic records in training data. (2) Models and strategies to focus genomic predictions on markers closer to the causative variants. Combining populations into a joint reference population results in high improvements when combining populations of the same breed and diminishes as the genetic distance between populations increases. For distantly related breeds sophisticated Bayesian variable selection models in combination with denser markers sets or functional subsets of markers is needed. This is expected to be further improved by the efficient use of sequence information. In addition predictions can be improved by the use of phenotypes of genotyped and non-genotyped cows directly. For a small population the optimal approach will combine the above components.     

Identification of a QTL on BTA20 affecting susceptibility to Mycobacterium avium ssp. paratuberculosis infection in US Holsteins

The objective of this study was to identify QTL affecting susceptibility to Mycobacterium paratuberculosis infection in US Holsteins. Twelve paternal half-sib families were selected for the study based on large numbers of daughters in production and limited relationships among sires. Serum and faecal samples from 4350 daughters of these 12 sires were obtained for disease testing. Case definition for an infected cow was an ELISA sample-to-positive ratio >/=0.25, a positive faecal culture or both. Three families were selected for genotyping based on a high apparent prevalence (6.8-10.4% infected cows), high faecal culture prevalence (46.2-52.9% positive faecal cultures) and large numbers of daughters tested for disease (264-585). DNA pooling was used to genotype cows, with an average of 159 microsatellites within each sire family. Infected cows (the positive pool) were matched with two of their non-infected herdmates in the same lactation (the negative pool) to control for herd and age effects. Eight chromosomal regions putatively linked with susceptibility to M. paratuberculosis infection were identified using a Z-test (P < 0.01). Significant results were more rigorously tested by individually genotyping cows with three to five informative microsatellites within 15 cM of the significant markers identified with the DNA pools. Probability of infection based on both diagnostic tests was estimated for each individual and used as the dependent variable for interval mapping. Based on this analysis, evidence for the presence of a QTL segregating within families on BTA20 was found (chromosome-wide P-value = 0.0319).     

Estimation of genetic parameters for young stock survival in beef x dairy crossbred calves

Young stock survival is a trait of crucial importance in cattle breeding as calf mortality leads to economic losses and represents an animal welfare issue. The aim of this study was to estimate genetic parameters and sire breeding values for young stock survival in beef x dairy crossbred calves. Two traits were analysed with a univariate animal model: young stock survival between 1 to 30 days and 31 to 200 days after birth. Breed combinations with Belgian Blue sires outperformed all other sire breeds. The lowest survival rates were found for breed combinations with Jersey dams or Blonde d’Aquitaine sires. The results showed low but significant heritabilities (0.045 to 0.075) for both survival traits. Differences in breeding values between sires ranged from −2.5% to 3.5% and from −5.4% to 4.7% survival from 1 to 30 days and 31 to 200 days, respectively. Based on these findings, we concluded that it is feasible to breed for improved young stock survival in beef x dairy crossbred calves. This will hopefully contribute to increasing the survival rate of the calves and reduce economic losses for the farmers.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

Comparison of piecewise Weibull baseline survival models for estimation of true and functional longevity in Brown cattle raised in small herds

Piecewise Weibull proportional hazard models were used to investigate the effect of genetic and nongenetic factors on functional and true longevity traits of the Slovenian Brown cattle breed. Records of 37 908 Brown cows from 2401 Slovenian herds were used. As these herds were characterised by a relatively small average herd size starting from 6.7 in 1999 and increasing to 8.7 Brown cows per herd in 2008, milk yield classification was made within different herd size groups. The hazard rate was the lowest in the first part of each lactation and was increasing for later stages. Culling risk was lower for cows from herds increasing in size, for cows with higher milk production and for cows from a region with smaller herd sizes and tougher conditions for cattle breeding. The latter result is surprising and may be related to better attention to maintain the animals, despite their lower milk production. The introduction of the milk quota system and drought was found to have an important effect on culling policy between the last seasons of the years 2001 and 2003. Seasonal effects were not related to the milk quota year (from April to March), but to the effect of shortage in fodder during the winter time. The effect of age at first calving and the interaction between year and milk yield class were not found to be significant. Heritability for functional and for true longevity were similar at around 10% each. Inclusion of a correction for class of milk yield to approximate functional longevity increased the herd-year random effect variance by 53%, whereas the sire variance increased by only 14%. The correlation coefficient between ranks of breeding values for functional and true longevity was high (0.91), whereas genetic trends were not found to be significant. To assess their predictive ability, models were compared looking at the survival rate of 4212 second-crop daughters not included in the initial models. The average correlation between estimated breeding values and survival at different stages was 0.39 for true longevity and 0.43 for functional longevity. Results showed that ranking milk yield at population level is appropriate to correct for voluntary culling on low production in small herds.     

Genetic association between body energy measured throughout lactation and fertility in dairy cattle

The objective of this study was to quantify the genetic association of body energy assessed throughout lactation with a cow's fertility. Nine direct and indirect body energy traits were defined at different stages of lactation. Four were daily records of energy balance, energy content, cumulative effective energy (CEE) and body condition score (BCS) calculated between lactation days 4 and 311. The other five traits included duration of negative energy balance (DNEB), rate of recovery during DNEB (RNEB), sum of negative energy balance (SNEB), nadir of energy content (NEC) and number of days from calving to NEC. Of these traits, energy balance, DNEB, RNEB and SNEB were primarily based on individual cow feed intake and milk yield, and considered direct measures of body energy. The other traits were calculated from body lipid and protein changes, predicted from BCS and live weight profiles, and were considered indirect measures of body energy. Fertility was defined by number of days between calving and commencement of luteal activity (DLA), first observed oestrus (DH) and conception (DC), and number of services per conception. A total of 957 cows in their first four lactations were considered in the study. Genetic models fitted cubic splines to define longitudinal traits (energy balance, energy content, CEE and BCS) and calculate heritability and genetic correlation with fertility. Daily heritability estimate ranges were 0.10 to 0.34, 0.35 to 0.61, 0.32 to 0.53 and 0.24 to 0.56 for energy balance, energy content, CEE and BCS, respectively, and, in most cases, tended to increase towards the middle of lactation and remain relatively stable thereafter. Of the other body energy traits, heritability of NEC (0.44) was the most notable. Statistically significant (P < 0.05) genetic correlations of DH with daily energy balance, energy content, CEE and BCS ranged from -0.16 to -0.28, -0.35 to -0.48, -0.16 to -0.26 and -0.37 to -0.44, respectively. For DC, respective estimates were -0.28 to -0.64, -0.37 to -0.60, -0.30 to -0.48 and -0.29 to -0.53. For DLA, they ranged from -0.47 to -0.56 with energy content and from -0.50 to -0.74 with BCS. Of special interest was the genetic correlation of NEC with DH (-0.54) and DC (-0.48). Results suggest that indirect measures of body energy have the strongest genetic association with cow fertility. NEC and early lactation (circa day 50) BCS and energy content are the most useful traits for selection in terms of the correlated improvement in a cow's capacity to resume her reproductive activity post partum.