Genetic parameters for carcase measurements and age at slaughter in commercial cattle.

Genetic parameters were estimated for cold carcase weight (CCW), carcase conformation (CON), carcase fat class (FAT), age at slaughter (AGE) and average daily carcase gain (ADCG) in 14 common UK breeds of cattle. These included crossbred animals but purebred datasets were also analysed for the most populous sire-breeds. Heritability estimates for beef breeds that were significant ranged from 0.24 to 0.44, 0.12 to 0.35, 0.12 to 0.36, 0.15 to 0.38 and 0.26 to 0.43 for CCW, CON, FAT, AGE and ADCG, respectively. For Holstein-Friesian, a dairy breed, heritability estimates were consistently lower than most beef breeds with estimates of 0.12, 0.13, 0.13, 0.06 and 0.15 for CCW, CON, FAT, AGE and ADCG, respectively. In all breed groups, genetic correlations were positive between CCW, CON and ADCG. In general, genetic correlations were moderate between CCW and CON (0.13 to 0.77), moderate to strong between CCW and ADCG (0.57 to 0.98) and weak or moderate between CON and ADCG (0.12 to 0.82). Genetic correlations for FAT with CCW (- 0.20 to - 0.42) and CON (- 0.16 to - 0.52) tended to be negative in the beef breed but were positive in the dairy breed, although not significant between CCW and FAT. For most beef breeds genetic correlations between AGE and carcase traits were not significant with the exceptions of AGE and CCW for Simmental (- 0.15) and Salers (- 0.24), AGE and CON for Limousin (0.15) and Simmental (0.14) and AGE and FAT from three sire-breeds (- 0.17 to - 0.35). However, the correlation between AGE and ADCG was negative and moderate to strong in magnitude (- 0.23 to - 0.67) in all beef breeds as expected since faster-growing animals reach slaughter age earlier. For Holstein-Friesian, all genetic correlations with AGE were negative and moderate to strong. Genetic correlations indicate that selection for increased carcase weight should simultaneously increase growth rate and improve conformation in all breeds and reduce carcase fatness in the majority of beef breeds. The results indicate that there is genetic variation in all five traits suitable for undertaking genetic improvement of carcase traits and age at slaughter; however, there are apparent breed differences. The use of abattoir-derived phenotypes for undertaking genetic improvement is an example where the supply chain can work together to share information to enable the cattle industry to move forward.

Predicting pregnancy status from mid-infrared spectroscopy in dairy cow milk using deep learning.

Accurately identifying pregnancy status is imperative for a profitable dairy enterprise. Mid-infrared (MIR) spectroscopy is routinely used to determine fat and protein concentrations in milk samples. Mid-infrared spectra have successfully been used to predict other economically important traits, including fatty acid content, mineral content, body energy status, lactoferrin, feed intake, and methane emissions. Machine learning has been used in a variety of fields to find patterns in vast quantities of data. This study aims to use deep learning, a sub-branch of machine learning, to establish pregnancy status from routinely collected milk MIR spectral data. Milk spectral data were obtained from National Milk Records (Chippenham, UK), who collect large volumes of data continuously on a monthly basis. Two approaches were followed: using genetic algorithms for feature selection and network design (model 1), and transfer learning with a pretrained DenseNet model (model 2). Feature selection in model 1 showed that the number of wave points in MIR data could be reduced from 1,060 to 196 wave points. The trained model converged after 162 epochs with validation accuracy and loss of 0.89 and 0.18, respectively. Although the accuracy was sufficiently high, the loss (in terms of predicting only 2 labels) was considered too high and suggested that the model would not be robust enough to apply to industry. Model 2 was trained in 2 stages of 100 epochs each with spectral data converted to gray-scale images and resulted in accuracy and loss of 0.97 and 0.08, respectively. Inspection on inference data showed prediction sensitivity of 0.89, specificity of 0.86, and prediction accuracy of 0.88. Results indicate that milk MIR data contains features relating to pregnancy status and the underlying metabolic changes in dairy cows, and such features can be identified by means of deep learning. Prediction equations from trained models can be used to alert farmers of nonviable pregnancies as well as to verify conception dates.

Predicting bovine tuberculosis status of dairy cows from mid-infrared spectral data of milk using deep learning.

Bovine tuberculosis (bTB) is a zoonotic disease in cattle that is transmissible to humans, distributed worldwide, and considered endemic throughout much of England and Wales. Mid-infrared (MIR) analysis of milk is used routinely to predict fat and protein concentration, and is also a robust predictor of several other economically important traits including individual fatty acids and body energy. This study predicted bTB status of UK dairy cows using their MIR spectral profiles collected as part of routine milk recording. Bovine tuberculosis data were collected as part of the national bTB testing program for Scotland, England, and Wales; these data provided information from over 40,500 bTB herd breakdowns. Corresponding individual cow life-history data were also available and provided information on births, movements, and deaths of all cows in the study. Data relating to single intradermal comparative cervical tuberculin (SICCT) skin-test results, culture, slaughter status, and presence of lesions were combined to create a binary bTB phenotype labeled 0 to represent nonresponders (i.e., healthy cows) and 1 to represent responders (i.e., bTB-affected cows). Contemporaneous individual milk MIR spectral data were collected as part of monthly routine milk recording and matched to bTB status of individual animals on the single intradermal comparative cervical tuberculin test date (±15 d). Deep learning, a sub-branch of machine learning, was used to train artificial neural networks and develop a prediction pipeline for subsequent use in national herds as part of routine milk recording. Spectra were first converted to 53 × 20-pixel PNG images, then used to train a deep convolutional neural network. Deep convolutional neural networks resulted in a bTB prediction accuracy (i.e., the number of correct predictions divided by the total number of predictions) of 71% after training for 278 epochs. This was accompanied by both a low validation loss (0.71) and moderate sensitivity and specificity (0.79 and 0.65, respectively). To balance data in each class, additional training data were synthesized using the synthetic minority over sampling technique. Accuracy was further increased to 95% (after 295 epochs), with corresponding validation loss minimized (0.26), when synthesized data were included during training of the network. Sensitivity and specificity also saw a 1.22- and 1.45-fold increase to 0.96 and 0.94, respectively, when synthesized data were included during training. We believe this study to be the first of its kind to predict bTB status from milk MIR spectral data. We also believe it to be the first study to use milk MIR spectral data to predict a disease phenotype, and posit that the automated prediction of bTB status at routine milk recording could provide farmers with a robust tool that enables them to make early management decisions on potential reactor cows, and thus help slow the spread of bTB.

Energy profiling of dairy cows from routine milk mid-infrared analysis.

The balance of body energy within and across lactations can have health and fertility consequences for the dairy cow. This study aimed to create a large calibration data set of dairy cow body energy traits across the cow's productive life, with concurrent milk mid-infrared (MIR) spectral data, to generate a prediction tool for use in commercial dairy herds. Detailed phenotypic data from 1,101 Holstein Friesian cows from the Langhill research herd (SRUC, Scotland) were used to generate energy balance (EB) and effective energy intake (EI), both in megajoules per day. Pretreatment of spectral data involved standardization to account for drift over time and machine. Body energy estimates were aligned with their spectral data to generate a prediction of these traits based on milk MIR spectroscopy. After data edits, partial least squares analysis generated prediction equations with a coefficient of determination from split sample 10-fold cross validation of 0.77 and 0.75 for EB and EI, respectively. These prediction equations were applied to national milk MIR spectra on over 11 million animal test dates (January 2013 to December 2016) from 4,453 farms. The predictions generated from these were subject to phenotypic analyses with a fixed regression model highlighting differences between the main dairy breeds in terms of energy traits. Genetic analyses generated heritability estimates for EB and EI ranging from 0.12 to 0.17 and 0.13 to 0.15, respectively. This study shows that MIR-based predictions from routinely collected national data can be used to generate predictions of dairy cow energy turnover profiles for both animal management and genetic improvement of such difficult and expensive-to-record traits.

The genetic and biological basis of feed efficiency in mid-lactation Holstein dairy cows.

The objective of this study was to identify genomic regions and candidate genes associated with feed efficiency in lactating Holstein cows. In total, 4,916 cows with actual or imputed genotypes for 60,671 single nucleotide polymorphisms having individual feed intake, milk yield, milk composition, and body weight records were used in this study. Cows were from research herds located in the United States, Canada, the Netherlands, and the United Kingdom. Feed efficiency, defined as residual feed intake (RFI), was calculated within location as the residual of the regression of dry matter intake (DMI) on milk energy (MilkE), metabolic body weight (MBW), change in body weight, and systematic effects. For RFI, DMI, MilkE, and MBW, bivariate analyses were performed considering each trait as a separate trait within parity group to estimate variance components and genetic correlations between them. Animal relationships were established using a genomic relationship matrix. Genome-wide association studies were performed separately by parity group for RFI, DMI, MilkE, and MBW using the Bayes B method with a prior assumption that 1% of single nucleotide polymorphisms have a nonzero effect. One-megabase windows with greatest percentage of the total genetic variation explained by the markers (TGVM) were identified, and adjacent windows with large proportion of the TGVM were combined and reanalyzed. Heritability estimates for RFI were 0.14 (±0.03; ±SE) in primiparous cows and 0.13 (±0.03) in multiparous cows. Genetic correlations between primiparous and multiparous cows were 0.76 for RFI, 0.78 for DMI, 0.92 for MBW, and 0.61 for MilkE. No single 1-Mb window explained a significant proportion of the TGVM for RFI; however, after combining windows, significance was met on Bos taurus autosome 27 in primiparous cows, and nearly reached on Bos taurus autosome 4 in multiparous cows. Among other genes, these regions contain ß-3 adrenergic receptor and the physiological candidate gene, leptin, respectively. Between the 2 parity groups, 3 of the 10 windows with the largest effects on DMI neighbored windows affecting RFI, but were not in the top 10 regions for MilkE or MBW. This result suggests a genetic basis for feed intake that is unrelated to energy consumption required for milk production or expected maintenance as determined by MBW. In conclusion, feed efficiency measured as RFI is a polygenic trait exhibiting a dynamic genetic basis and genetic variation distinct from that underlying expected maintenance requirements and milk energy output.

Genetic parameters of Visual Image Analysis primal cut carcass traits of commercial prime beef slaughter animals.

Visual Image analysis (VIA) of carcass traits provides the opportunity to estimate carcass primal cut yields on large numbers of slaughter animals. This allows carcases to be better differentiated and farmers to be paid based on the primal cut yields. It also creates more accurate genetic selection due to high volumes of data which enables breeders to breed cattle that better meet the abattoir specifications and market requirements. In order to implement genetic evaluations for VIA primal cut yields, genetic parameters must first be estimated and that was the aim of this study. Slaughter records from the UK prime slaughter population for VIA carcass traits was available from two processing plants. After edits, there were 17 765 VIA carcass records for six primal cut traits, carcass weight as well as the EUROP conformation and fat class grades. Heritability estimates after traits were adjusted for age ranged from 0.32 (0.03) for EUROP fat to 0.46 (0.03) for VIA Topside primal cut yield. Adjusting the VIA primal cut yields for carcass weight reduced the heritability estimates, with estimates of primal cut yields ranging from 0.23 (0.03) for Fillet to 0.29 (0.03) for Knuckle. Genetic correlations between VIA primal cut yields adjusted for carcass weight were very strong, ranging from 0.40 (0.06) between Fillet and Striploin to 0.92 (0.02) between Topside and Silverside. EUROP conformation was also positively correlated with the VIA primal cuts with genetic correlation estimates ranging from 0.59 to 0.84, whereas EUROP fat was estimated to have moderate negative correlations with primal cut yields, estimates ranged from -0.11 to -0.46. Based on these genetic parameter estimates, genetic evaluation of VIA primal cut yields can be undertaken to allow the UK beef industry to select carcases that better meet abattoir specification and market requirements.

Genetic evaluation for bovine tuberculosis resistance in dairy cattle.

Genetic evaluations for resistance to bovine tuberculosis (bTB) were calculated based on British national data including individual animal tuberculin skin test results, postmortem examination (presence of bTB lesions and bacteriological culture for Mycobacterium bovis), animal movement and location information, production history, and pedigree records. Holstein cows with identified sires in herds with bTB breakdowns (new herd incidents) occurring between the years 2000 and 2014 were considered. In the first instance, cows with a positive reaction to the skin test and a positive postmortem examination were defined as infected. Values of 0 and 1 were assigned to healthy and infected animal records, respectively. Data were analyzed with mixed models. Linear and logit function heritability estimates were 0.092 and 0.172, respectively. In subsequent analyses, breakdowns were split into 2-mo intervals to better model time of exposure and infection in the contemporary group. Intervals with at least one infected individual were retained and multiple intervals within the same breakdown were included. Healthy animal records were assigned values of 0, and infected records a value of 1 in the interval of infection and values reflecting a diminishing probability of infection in the preceding intervals. Heritability and repeatability estimates were 0.115 and 0.699, respectively. Reliabilities and across time stability of the genetic evaluation were improved with the interval model. Subsequently, 2 more definitions of "infected" were analyzed with the interval model: (1) all positive skin test reactors regardless of postmortem examination, and (2) all positive skin test reactors plus nonreactors with positive postmortem examination. Estimated heritability was 0.085 and 0.089, respectively; corresponding repeatability estimates were 0.701 and 0.697. Genetic evaluation reliabilities and across time stability did not change. Correlations of genetic evaluations for bTB with other traits in the current breeding goal were mostly not different from zero. Correlation with the UK Profitable Lifetime Index was moderate, significant, and favorable. Results demonstrated the feasibility of a national genetic evaluation for bTB resistance. Selection for enhanced resistance will have a positive effect on profitability and no antagonistic effects on current breeding goal traits. Official genetic evaluations are now based on the interval model and the last bTB trait definition.

Genetic analysis of carcass traits in beef cattle using random regression models.

Livestock mature at different rates depending, in part, on their genetic merit; therefore, the optimal age at slaughter for progeny of certain sires may differ. The objective of the present study was to examine sire-level genetic profiles for carcass weight, carcass conformation, and carcass fat in cattle of multiple beef and dairy breeds, including crossbreeds. Slaughter records from 126,214 heifers and 124,641 steers aged between 360 and 1,200 d and from 86,089 young bulls aged between 360 and 720 d were used in the analysis; animals were from 15,127 sires. Variance components for each trait across age at slaughter were generated using sire random regression models that included quadratic polynomials for fixed and random effects; heterogeneous residual variances were assumed across ages. Heritability estimates across genders ranged from 0.08 (±0.02) to 0.34 (±0.02) for carcass weight, from 0.24 (±0.02) to 0.42 (±0.01) for conformation, and from 0.16 (±0.03) to 0.40 (±0.02) for fat score. Genetic correlations within each trait across ages weakened as the interval between ages compared lengthened but were all >0.64, suggesting a similar genetic background for each trait across different ages. Eigenvalues and eigenfunctions of the additive genetic covariance matrix revealed genetic variability among animals in their growth profiles for carcass traits, although most of the genetic variability was associated with the height of the growth profile. At the same age, a positive genetic correlation (0.60 to 0.78; SE ranged from 0.01 to 0.04) existed between carcass weight and conformation, whereas negative genetic correlations existed between fatness and both conformation (-0.46 to 0.08; SE ranged from 0.02 to 0.09) and carcass weight (-0.48 to -0.16; SE ranged from 0.02 to 0.14) at the same age. The estimated genetic parameters in the present study indicate genetic variability in the growth trajectory in cattle, which can be exploited through breeding programs and used in decision support tools.

Hot topic: Definition and implementation of a breeding value for feed efficiency in dairy cows.

A new breeding value that combines the amount of feed saved through improved metabolic efficiency with predicted maintenance requirements is described. The breeding value includes a genomic component for residual feed intake (RFI) combined with maintenance requirements calculated from either a genomic or pedigree estimated breeding value (EBV) for body weight (BW) predicted using conformation traits. Residual feed intake is only available for genotyped Holsteins; however, BW is available for all breeds. The RFI component of the "feed saved" EBV has 2 parts: Australian calf RFI and Australian lactating cow RFI. Genomic breeding values for RFI were estimated from a reference population of 2,036 individuals in a multi-trait analysis including Australian calf RFI (n=843), Australian lactating cow RFI (n=234), and UK and Dutch lactating cow RFI (n=958). In all cases, the RFI phenotypes were deviations from a mean of 0, calculated by correcting dry matter intake for BW, growth, and milk yield (in the case of lactating cows). Single nucleotide polymorphism effects were calculated from the output of genomic BLUP and used to predict breeding values of 4,106 Holstein sires that were genotyped but did not have RFI phenotypes themselves. These bulls already had BW breeding values calculated from type traits, from which maintenance requirements in kilograms of feed per year were inferred. Finally, RFI and the feed required for maintenance (through BW) were used to calculate a feed saved breeding value and expressed as the predicted amount of feed saved per year. Animals that were 1 standard deviation above the mean were predicted to eat 66 kg dry matter less per year at the same level of milk production. In a data set of genotyped Holstein sires, the mean reliability of the feed saved breeding value was 0.37. For Holsteins that are not genotyped and for breeds other than Holsteins, feed saved is calculated using BW only. From April 2015, feed saved has been included as part of the Australian national selection index, the Balanced Performance Index (BPI). Selection on the BPI is expected to lead to modest gains in feed efficiency.

Transcutaneous bilirubin levels in newborns <35 weeks' gestation.

OBJECTIVE: In infants <35 weeks' gestation, we sought to define the transcutaneous bilirubin (TcB) levels at which a total serum bilirubin (TSB) level suggesting the need for phototherapy is unlikely to occur and a TSB measurement can, therefore, be avoided. STUDY DESIGN: Nursing staff performed 896 TcB measurements within 1 h of a TSB on 225 neonates 26 0/7-34 6/7 weeks' postmenstrual age (PMA). Generalized linear models were fit with generalized estimating equations (GEEs) to model the probability of having a TSB level at or above the phototherapy initiation cutpoint as a function of the TcB; these methods allow for multiple tests per infant. RESULTS: The mean difference between TcB and TSB measurements was <1 mg dl(-1) for each PMA category. When the TcB was at least 3 mg dl(-1) below the TSB cutpoint for phototherapy, there was a ⩾98% probability that the TSB was not at, or above, the recommended phototherapy level. The single exception to this was a phototherapy level of 6 mg dl(-1) for infants of 28 0/7-29 6/7 weeks' PMA, where a TcB of 4 mg dl(-1) below the phototherapy level (ie a TcB ⩽2 mg dl(-1)) was necessary to achieve ⩾98% probability. CONCLUSION: Our data support the use of routine TcB screening for infants 28-34 6/7 weeks' gestation. TcB screening in the neonatal intensive care unit can identify infants who require a TSB to confirm or exclude the need for phototherapy.

On the value of the phenotypes in the genomic era.

Genetic improvement programs around the world rely on the collection of accurate phenotypic data. These phenotypes have an inherent value that can be estimated as the contribution of an additional record to genetic gain. Here, the contribution of phenotypes to genetic gain was calculated using traditional progeny testing (PT) and 2 genomic selection (GS) strategies that, for simplicity, included either males or females in the reference population. A procedure to estimate the theoretical economic contribution of a phenotype to a breeding program is described for both GS and PT breeding programs through the increment in genetic gain per unit of increase in estimated breeding value reliability obtained when an additional phenotypic record is added. The main factors affecting the value of a phenotype were the economic value of the trait, the number of phenotypic records already available for the trait, and its heritability. Furthermore, the value of a phenotype was affected by several other factors, including the cost of establishing the breeding program and the cost of phenotyping and genotyping. The cost of achieving a reliability of 0.60 was assessed for different reference populations for GS. Genomic reference populations of more sires with small progeny group sizes (e.g., 20 equivalent daughters) had a lower cost than those reference populations with either large progeny group sizes for fewer genotyped sires, or female reference populations, unless the heritability was large and the cost of phenotyping exceeded a few hundred dollars; then, female reference populations were preferable from an economic perspective.

International genetic evaluations for feed intake in dairy cattle through the collation of data from multiple sources.

Feed represents a large proportion of the variable costs in dairy production systems. The omission of feed intake measures explicitly from national dairy cow breeding objectives is predominantly due to a lack of information from which to make selection decisions. However, individual cow feed intake data are available in different countries, mostly from research or nucleus herds. None of these data sets are sufficiently large enough on their own to generate accurate genetic evaluations. In the current study, we collate data from 10 populations in 9 countries and estimate genetic parameters for dry matter intake (DMI). A total of 224,174 test-day records from 10,068 parity 1 to 5 records of 6,957 cows were available, as well as records from 1,784 growing heifers. Random regression models were fit to the lactating cow test-day records and predicted feed intake at 70 d postcalving was extracted from these fitted profiles. The random regression model included a fixed polynomial regression for each lactation separately, as well as herd-year-season of calving and experimental treatment as fixed effects; random effects fit in the model included individual animal deviation from the fixed regression for each parity as well as mean herd-specific deviations from the fixed regression. Predicted DMI at 70 d postcalving was used as the phenotype for the subsequent genetic analyses undertaken using an animal repeatability model. Heritability estimates of predicted cow feed intake 70 d postcalving was 0.34 across the entire data set and varied, within population, from 0.08 to 0.52. Repeatability of feed intake across lactations was 0.66. Heritability of feed intake in the growing heifers was 0.20 to 0.34 in the 2 populations with heifer data. The genetic correlation between feed intake in lactating cows and growing heifers was 0.67. A combined pedigree and genomic relationship matrix was used to improve linkages between populations for the estimation of genetic correlations of DMI in lactating cows; genotype information was available on 5,429 of the animals. Populations were categorized as North America, grazing, other low input, and high input European Union. Albeit associated with large standard errors, genetic correlation estimates for DMI between populations varied from 0.14 to 0.84 but were stronger (0.76 to 0.84) between the populations representative of high-input production systems. Genetic correlations with the grazing populations were weak to moderate, varying from 0.14 to 0.57. Genetic evaluations for DMI can be undertaken using data collated from international populations; however, genotype-by-environment interactions with grazing production systems need to be considered.

Short communication: Validation of genomic breeding value predictions for feed intake and feed efficiency traits.

Validating genomic prediction equations in independent populations is an important part of evaluating genomic selection. Published genomic predictions from 2 studies on (1) residual feed intake and (2) dry matter intake (DMI) were validated in a cohort of 78 multiparous Holsteins from Australia. The mean realized accuracy of genomic prediction for residual feed intake was 0.27 when the reference population included phenotypes from 939 New Zealand and 843 Australian growing heifers (aged 5-8 mo) genotyped on high density (770k) single nucleotide polymorphism chips. The 90% bootstrapped confidence interval of this estimate was between 0.16 and 0.36. The mean realized accuracy was slightly lower (0.25) when the reference population comprised only Australian growing heifers. Higher realized accuracies were achieved for DMI in the same validation population and using a multicountry model that included 958 lactating cows from the Netherlands and United Kingdom in addition to 843 growing heifers from Australia. The multicountry analysis for DMI generated 3 sets of genomic predictions for validation animals, one on each country scale. The highest mean accuracy (0.72) was obtained when the genomic breeding values were expressed on the Dutch scale. Although the validation population used in this study was small (n=78), the results illustrate that genomic selection for DMI and residual feed intake is feasible. Multicountry collaboration in the area of dairy cow feed efficiency is the evident pathway to achieving reasonable genomic prediction accuracies for these valuable traits.

Direct and maternal genetic relationships between calving ease, gestation length, milk production, fertility, type, and lifespan of Holstein-Friesian primiparous cows.

As the emphasis in cattle breeding is shifting from traits that increase income toward traits that reduce costs, national breeding indices are expanding to include functional traits such as calving ease (CE). However, one issue is the lack of knowledge of genetic relationships between CE and other dairy traits. The same can be said about gestation length (GL), a potential novel selection trait with considerable heritabilities and possible genetic relationships with the calving process. This study aimed to estimate the genetic relationships between CE, GL, and other dairy traits of interest using a national data set of 31,053 primiparous cow performance records, as well as to separate direct and maternal genetic effects. Chosen dairy traits included fertility (calving interval, days to first service, nonreturn rate after 56 d, number of inseminations per conception), milk production (milk yield at d 110 in milk, accumulated 305-d milk yield, accumulated 305-d fat yield, accumulated 305-d protein yield), type (udder depth, chest width, rump width, rump angle, mammary composition, stature, body depth), and lifespan traits (functional days of productive life). To allow the separation of direct and maternal genetic effects, a random sire of the calf effect was included in the multi-trait linear trivariate sire models fitted using ASReml. Significant results showed that easily born individuals were genetically prone to high milk yield and reduced fertility in first lactation. Difficult calving primiparous cows were likely associated with being high-producing, wide and deep animals, with a reduced ability to subsequently conceive. Individuals that were born relatively early were associated with good genetic merit for milk production. Finally, individuals carrying their offspring longer were genetically associated with being wide and large animals that were themselves born relatively early. The study shows that it is feasible and valuable to separate direct and maternal effects when estimating genetic correlations between calving and other dairy traits. Furthermore, gestation length is best used as an indicator trait for lowly heritable calving traits, rather than as a novel selection trait. As estimated direct and maternal genetic correlations differ, we can conclude that genetic relationships between CE, GL, and traits of interest are present, but caution is required if these traits are implemented in national breeding indices.

The sensitivity of predicted financial and genetic gains in Holsteins to changes in the economic value of traits.

The objective of this study was to assess the impact of using different relative economic values (REVs) in selection indices on predicted financial and trait gains from selection of sires of cows and on the choice of leading Holstein bulls available in the UK dairy industry. Breeding objective traits were milk yield, fat yield, protein yield, lifespan, mastitis, non-return rate, calving interval and lameness. Relative importance of a trait, as estimated by a.h(2), was only moderately related to the rate of financial loss or total economic merit (ΔTEM) per percentage under- or overestimation of REV (r = 0.38 and 0.29, respectively) as a result of the variance-covariance structure of traits. The effects on TEM of under- or overestimating trait REVs were non-symmetrical. TEM was most sensitive to incorrect REVs for protein, fat, milk and lifespan and least sensitive to incorrect calving interval, lameness, non-return and mastitis REVs. A guide to deciding which dairy traits require the most rigorous analysis in the calculation of their REVs is given. Varying the REVs within a fairly wide range resulted in different bulls being selected by index and their differing predicted transmitting abilities would result in the herds moving in different directions in the long term (20 years). It is suggested that customized indices, where the farmer creates rankings of bulls tailored to their specific farm circumstances, can be worthwhile.

Genome-wide associations for feed utilisation complex in primiparous Holstein-Friesian dairy cows from experimental research herds in four European countries.

Genome-wide association studies for difficult-to-measure traits are generally limited by the sample size with accurate phenotypic data. The objective of this study was to utilise data on primiparous Holstein­Friesian cows from experimental farms in Ireland, the United Kingdom, the Netherlands and Sweden to identify genomic regions associated with the feed utilisation complex: fat and protein corrected milk yield (FPCM), dry matter intake (DMI), body condition score (BCS) and live-weight (LW). Phenotypic data and 37 590 single nucleotide polymorphisms (SNPs) were available on up to 1629 animals. Genetic parameters of the traits were estimated using a linear animal model with pedigree information, and univariate genome-wide association analyses were undertaken using Bayesian stochastic search variable selection performed using Gibbs sampling. The variation in the phenotypes explained by the SNPs on each chromosome was related to the size of the chromosome and was relatively consistent for each trait with the possible exceptions of BTA4 for BCS, BTA7, BTA13, BTA14, BTA18 for LW and BTA27 for DMI. For LW, BCS, DMI and FPCM, 266, 178, 206 and 254 SNPs had a Bayes factor .3, respectively. Olfactory genes and genes involved in the sensory smell process were overrepresented in a 500 kbp window around the significant SNPs. Potential candidate genes were involved with functions linked to insulin, epidermal growth factor and tryptophan.

The effect of lactation length on greenhouse gas emissions from the national dairy herd.

Many governments have signed up to greenhouse gas emission (GHGE) reduction programmes under their national climate change obligations. Recently, it has been suggested that the use of extended lactations in dairy herds could result in reduced GHGE. Dairy GHGE were modelled on a national basis and the model was used to compare emissions from lactations of three different lengths (305, 370 and 440 days), and a current 'base' scenario on the basis of maintaining current milk production levels. In addition to comparing GHGE from the average 'National Herd' under these scenarios, results were used to investigate how accounting for lactations of different lengths might alter the estimation of emissions calculated from the National Inventory methodology currently recommended by Intergovernmental Panel on Climate Change. Data for the three lactation length scenarios were derived from nationally recorded dairy performance information and used in the GHGE model. Long lactations required fewer milking cows and replacements to maintain current milk yield levels than short ones, but GHGEs were found to rise from 1214 t of CO2 equivalent (CE)/farm per year for lactations of 305 days to 1371 t CE/farm per year for 440-day lactations. This apparent anomaly can be explained by the less efficient milk production (kg milk produced per kg cow weight) found in later lactation, a more pronounced effect in longer lactations. The sensitivity of the model to changes in replacement rate, persistency and level of milk yield was investigated. Changes in the replacement rate from 25% to 20% and in persistency by −10% to +20% resulted in very small changes in GHGE. Differences in GHGE due to the level of milk yield were much more dramatic with animals in the top 10% for yield, producing about 25% less GHGE/year than the average animal. National Inventory results were investigated using a more realistic spread of lactation lengths than recommended for such calculations using emissions calculated in the first part of the study. Current UK emission calculations based on the National Inventory were 329 Gg of methane per year from the dairy herd. Using the national distribution of lactation lengths, this was found to be an underestimate by about 10%. This work showed that the current rise in lactation length or a move towards calving every 18 months would increase GHGE by 7% to 14% compared with the current scenario, assuming the same milk yield in all models. Increased milk yield would have a much greater effect on reducing GHGE than changes to lactation length, replacement rate or persistency. National Inventory methodology appears to underestimate GHGE when the distribution of lactation lengths is considered and may need revising to provide more realistic figures.

Merging and characterising phenotypic data on conventional and rare traits from dairy cattle experimental resources in three countries.

This study set out to demonstrate the feasibility of merging data from different experimental resource dairy populations for joint genetic analyses. Data from four experimental herds located in three different countries (Scotland, Ireland and the Netherlands) were used for this purpose. Animals were first lactation Holstein cows that participated in ongoing or previously completed selection and feeding experiments. Data included a total of 60 058 weekly records from 1630 cows across the four herds; number of cows per herd ranged from 90 to 563. Weekly records were extracted from the individual herd databases and included seven traits: milk, fat and protein yield, milk somatic cell count, liveweight, dry matter intake and energy intake. Missing records were predicted with the use of random regression models, so that at the end there were 44 weekly records, corresponding to the typical 305-day lactation, for each cow. A total of 23 different lactation traits were derived from these records: total milk, fat and protein yield, average fat and protein percentage, average fat-to-protein ratio, total dry matter and energy intake and average dry matter intake-to-milk yield ratio in lactation weeks 1 to 44 and 1 to 15; average milk somatic cell count in lactation weeks 1 to 15 and 16 to 44; average liveweight in lactation weeks 1 to 44; and average energy balance in lactation weeks 1 to 44 and 1 to 15. Data were subsequently merged across the four herds into a single dataset, which was analysed with mixed linear models. Genetic variance and heritability estimates were greater (P < 0.05) than zero for all traits except for average milk somatic cell count in weeks 16 to 44. Proportion of total phenotypic variance due to genotype-by-environment (sire-by-herd) interaction was not different (P > 0.05) from zero. When estimable, the genetic correlation between herds ranged from 0.85 to 0.99. Results suggested that merging experimental herd data into a single dataset is both feasible and sensible, despite potential differences in management and recording of the animals in the four herds. Merging experimental data will increase power of detection in a genetic analysis and augment the potential reference population in genome-wide association studies, especially of difficult-to-record traits.