Predicting methane emission in Canadian Holstein dairy cattle using milk mid-infrared reflectance spectroscopy and other commonly available predictors via artificial neural networks

Interest in reducing eructed CH₄ is growing, but measuring CH₄ emissions is expensive and difficult in large populations. In this study, we investigated the effectiveness of milk mid-infrared spectroscopy (MIRS) data to predict CH₄ emission in lactating Canadian Holstein cows. A total of 181 weekly average CH₄ records from 158 Canadian cows and 217 records from 44 Danish cows were used. For each milk spectra record, the corresponding weekly average CH₄ emission (g/d), test-day milk yield (MY, kg/d), fat yield (FY, g/d), and protein yield (PY, g/d) were available. The weekly average CH₄ emission was predicted using various artificial neural networks (ANN), partial least squares regression, and different sets of predictors. The ANN architectures consisted of 3 training algorithms, 1 to 10 neurons with hyperbolic tangent activation function in the hidden layer, and 1 neuron with linear (purine) activation function in the hidden layer. Random cross-validation was used to compared the predictor sets: MY (set 1); FY (set 2); PY (set 3); MY and FY (set 4); MY and PY (set 5); MY, FY, and PY (set 6); MIRS (set 7); and MY, FY, PY, and MIRS (set 8). All predictor sets also included age at calving and days in milk, in addition to country, season of calving, and lactation number as categorical effects. Using only MY (set 1), the predictive accuracy (r) ranged from 0.245 to 0.457 and the root mean square error (RMSE) ranged from 87.28 to 99.39 across all prediction models and validation sets. Replacing MY with FY (set 2; r = 0.288–0.491; RMSE = 85.94–98.04) improved the predictive accuracy, but using PY (set 3; r = 0.260–0.468; RMSE = 86.95–98.47) did not. Adding FY to MY (set 4; r = 0.272–0.469; RMSE = 87.21–100.76) led to a negligible improvement compared with sets 1 and 3, but it slightly decreased accuracy compared with set 2. Adding PY to MY (set 5; r = 0.250–0.451; RMSE = 87.66–100.94) did not improve prediction ability. Combining MY, FY, and PY (set 6; r = 0.252–0.455; RMSE = 87.74–101.93) yielded accuracy slightly lower than sets 2 and 3. Using only MIRS data (set 7; r = 0.586–0.717; RMSE = 69.09–96.20) resulted in superior accuracy compared with all previous sets. Finally, the combination of MIRS data with MY, FY, and PY (set 8; r = 0.590–0.727; RMSE = 68.02–87.78) yielded similar accuracy to set 7. Overall, sets including the MIRS data yielded significantly better predictions than the other sets. To assess the predictive ability in a new unseen herd, a limited block cross-validation was performed using 20 cows in the same Canadian herd, which yielded r = 0.229 and RMSE = 154.44, which were clearly much worse than the average r = 0.704 and RMSE = 70.83 when predictions were made by random cross-validation. These results warrant further investigation when more data become available to allow for a more comprehensive block cross-validation before applying the calibrated models for large-scale prediction of CH₄ emissions.     

Genetic analysis of milk urea nitrogen and relationships with yield and fertility across lactation

The aim of this project was to investigate the relationship of milk urea nitrogen (MUN) with 3 milk production traits [milk yield (MY), fat yield (FY), protein yield (PY)] and 6 fertility measures (number of inseminations, calving interval, interval from calving to first insemination, interval from calving to last insemination, interval from first to last insemination, and pregnancy at first insemination). Data consisted of 635,289 test-day records of MY, FY, PY, and MUN on 76,959 first-lactation Swedish Holstein cows calving from 2001 to 2003, and corresponding lactation records for the fertility traits. Yields and MUN were analyzed with a random regression model followed by a multi-trait model in which the lactation was broken into 10 monthly periods. Heritability for MUN was stable across lactation (between 0.16 and 0.18), whereas MY, FY, and PY had low heritability at the beginning of lactation, which increased with time and stabilized after 100 d in milk, at 0.47, 0.36, and 0.44, respectively. Fertility traits had low heritabilities (0.02 to 0.05). Phenotypic correlations of MUN and milk production traits were between 0.13 (beginning of lactation) and 0.00 (end of lactation). Genetic correlations of MUN and MY, FY, and PY followed similar trends and were positive (0.22) at the beginning and negative (−0.15) at the end of lactation. Phenotypic correlations of MUN and fertility were close to zero. A surprising result was that genetic correlations of MUN and fertility traits suggest a positive relationship between the 2 traits for most of the lactation, indicating that animals with breeding values for increased MUN also had breeding values for improved fertility. This result was obtained with a random regression model as well as with a multi-trait model. The analyzed group of cows had a moderate level of MUN concentration. In such a population MUN concentration may increase slightly due to selection for improved fertility. Conversely, selection for increased MUN concentration may improve fertility slightly.     

Random regression analysis of milk yield and milk composition in the first and second lactations of Murciano-Granadina goats

Records from the milk recording scheme of Spanish Murciano-Granadina goats were studied to estimate genetic (co)variance components and breeding values throughout the first and second lactations. The data used consisted of 49,696 monthly test-day records of milk (MY), protein (PY), fat (FY), and dry matter (DMY) yields from 5,163 goats, distributed in 20 herds, offspring of 2,086 does and 206 bucks. These records were analyzed by 2-trait random regression models (RRM) and a repeatability test-day model (RTDM). At the middle of lactation, heritability estimates for MY, DMY, and FY obtained with RTDM were larger than those estimated with RRM, and the opposite was true for PY. The RRM estimates of heritability for MY, FY, and PY were very similar throughout the trajectories of both lactations. Heritability estimates for DMY decreased through the lactation period. The genetic correlations between the first and second lactation records estimated for all traits by RRM were positive and ranged from 0.43 to 0.80 throughout the lactation curves. The correlation between BV estimated with RTDM and RRM was 0.742 for MY and 0.664 for DMY. The RRM could be a useful alternative to RTDM for the prediction of BV in this breed.     

Comparison of methane production,intensity, and yield throughout lactation in Holstein cows

Genetic selection to reduce methane (CH₄) emissions from dairy cows is an attractive means of reducing the impact of agricultural production on climate change. In this study, we investigated the feasibility of such an approach by characterizing the interactions between CH₄ and several traits of interest in dairy cows. We measured CH₄, dry matter intake (DMI), fat- and protein-corrected milk (FPCM), body weight (BW), and body condition score (BCS) from 107 first- and second-parity Holstein cows from December 2019 to November 2021. Methane emissions were measured using a GreenFeed device and expressed in terms of production (MeP, in g/d), yield (MeY, in g/kg DMI), and intensity (MeI, in g/kg FPCM). Because of the limited number of cows, only animal parameters were estimated. Both MeP and MeI were moderately repeatable (>0.45), whereas MeY presented low repeatability, especially in early lactation. Mid lactation was the most stable and representative period of CH₄ emissions throughout lactation, with animal correlations above 0.9. The average animal correlations of MeP with DMI, FPCM, and BW were 0.62, 0.48, and 0.36, respectively. The MeI was negatively correlated with FCPM (<−0.5) and DMI (>−0.25), and positively correlated with BW and BCS. The MeY presented stable and weakly positive correlations with the 4 other traits throughout lactation, with the exception of slightly negative animal correlations with FPCM and DMI after the 35th week. The MeP, MeI, and MeY were positively correlated at all lactation stages and, assuming animal and genetic correlations do not strongly differ, selection on one trait should lead to improvements in all. Overall, selection for MeI is probably not optimal as its change would result more from CH₄ dilution in increased milk yield than from real decrease in methane emission. Instead, MeY is related to rumen function and is only weakly associated with DMI, FPCM, BW, and BCS; it thus appears to be the most promising CH₄ trait for selection, provided that this would not deteriorate feed efficiency and that a system of large-scale phenotyping is developed. The MeP is easier to measure and thus may represent an acceptable alternative, although care would need to be taken to avoid undesirable changes in FPCM and BW.     

Is rumination time an indicator of methane production in dairy cows?

As long as large-scale recording of expensive-to-measure and labor-consuming traits, such as dry matter intake (DMI) and CH₄ production (CH₄P), continues to be challenging in practical conditions, alternative traits that are already routinely recorded in dairy herds should be investigated. An ideal indicator trait must, in addition to expressing genetic variation, have a strong correlation with the trait of interest. Our aim was to estimate individual level and phenotypic correlations between rumination time (RT), CH₄P, and DMI to determine if RT could be used as an indicator trait for CH₄P and DMI. Data from 343 Danish Holstein cows were collected at the Danish Cattle Research Centre for a period of approximately 3 yr. The data set consisted of 14,890 records for DMI, 15,835 for RT, and 6,693 for CH₄P. Data were divided in primiparous cows only (PC) and all cows (MC), and then divided in lactation stage (early, mid, late, and whole lactation) to analyze the changes over lactation. Linear mixed models, including an animal effect but no pedigree, were used to estimate the correlations among traits. Phenotypic and individual level correlations between RT and both CH₄P and DMI were close to zero, regardless of lactation stage and data set (PC or MC). However, CH₄P and DMI were highly correlated, both across lactation stages and data sets. In conclusion, RT is unsuitable to be used as an indicator trait for either CH₄P or DMI. Our study failed to validate RT as a useful indicator trait for both CH₄P and DMI, but more studies with novel phenotypes can offer different approaches to select and incorporate important yet difficult to record traits into breeding goals and selection indexes.     

Estimation of genetic parameters for productive life,reproduction, and milk-production traits in US dairy goats

Heritabilities and correlations for milk yield (MY), fat yield (FY), protein yield (PY), combined fat and protein yield (FPY), fat percentage (F%), protein percentage (P%), age at first kidding (AFK), interval between the first and second kidding (KI), and real and functional productive life at 72 mo (FPL72) of 33,725 US dairy goats, were estimated using animal models. Productive life was defined as the total days in production until 72 mo of age (PL72) for goats having the opportunity to express the trait. Functional productive life was obtained by correcting PL72 for MY, FY, PY, and final type score (FS). Six selection indexes were used, including or excluding PL72, with 6 groups of different economic weights, to estimate the responses to selection considering MY, FY, PY, and PL72 as selection criteria. The main criteria that determined the culling of a goat from the herd were low FS, MY, and FY per lactation. Heritability estimates were 0.22, 0.17, 0.37, 0.37, 0.38, 0.39, 0.54, 0.64, 0.09, and 0.16 for PL72, FPL72, MY, FY, PY, FPY, F%, P%, KI, and AFK, respectively. Most genetic correlations between the evaluated traits and PL72 or FPL72 were positive, except for F% (−0.04 and −0.06, respectively), P% (−0.002 and −0.03, respectively), and AFK (−0.03 and −0.01, respectively). The highest genetic correlations were between FPL72 and MY (0.39) and between PL72 and MY (0.33). Most phenotypic correlations between the traits evaluated and FPL72 and PL72 were positive (>0.23 and >0.26, respectively), except for F% (−0.004 and −0.02, respectively), P% (−0.05 and −0.02), KI (−0.01 and −0.07), and AFK (−0.08 and −0.08). The direct selection for PL72 increased it by 102.28 d per generation. The use of MY, FY, PY, KI, or AFK as selection criteria increased PL72 by 39.21, 27.33, 35.90, −8.28, or 2.77 d per generation, respectively. The inclusion of PL72 as selection criterion increased the expected response per generation from 0.15 to 17.35% in all selection indices studied.     

Comparative accuracies of genetic values predicted for economically important milk traits,genome-wide association,and linkage disequilibrium patterns of Canadian Holstein cows

Genomic selection methodologies and genome-wide association studies use powerful statistical procedures that correlate large amounts of high-density SNP genotypes and phenotypic data. Actual 305-d milk (MY), fat (FY), and protein (PY) yield data on 695 cows and 76,355 genotyping-by-sequencing-generated SNP marker genotypes from Canadian Holstein dairy cows were used to characterize linkage disequilibrium (LD) structure of Canadian Holstein cows. Also, the comparison of pedigree-based BLUP, genomic BLUP (GBLUP), and Bayesian (BayesB) statistical methods in the genomic selection methodologies and the comparison of Bayesian ridge regression and BayesB statistical methods in the genome-wide association studies were carried out for MY, FY, and PY. Results from LD analysis revealed that as marker distance decreases, LD increases through chromosomes. However, unexpected high peaks in LD were observed between marker pairs with larger marker distances on all chromosomes. The GBLUP and BayesB models resulted in similar heritability estimates through 10-fold cross-validation for MY and PY; however, the GBLUP model resulted in higher heritability estimates than BayesB model for FY. The predictive ability of GBLUP model was significantly lower than that of BayesB for MY, FY, and PY. Association analyses indicated that 28 high-effect markers and markers on Bos taurus autosome 14 located within 6 genes (DOP1B, TONSL, CPSF1, ADCK5, PARP10, and GRINA) associated significantly with FY.     

Genetic parameters for fertility of dairy heifers and cows at different parities and relationships with production traits in first lactation

The objectives of this study were to estimate genetic parameters for fertility of Brown Swiss cattle, considering reproductive measures in different parities as different traits, and to estimate relationships between production traits of first lactation and fertility of heifers and first-parity and second-parity cows. Reproductive indicators were interval from parturition to first service, interval from first service to conception, interval from parturition to conception, number of inseminations to conception, conception rate at first service, and nonreturn rate at 56 d after first service. Production traits were peak milk yield (pMY), lactation milk yield, and lactation length (LL). Data included 37,546 records on heifers, and 24,098 and 15,653 records on first- and second-parity cows, respectively. Cows were reared in 2,035 herds, calved from 1999 to 2007, and were progeny of 527 AI bulls. Gibbs sampling was implemented to obtain (co)variance components using both univariate and bivariate threshold and censored linear sire models. Estimates of heritability for reproductive traits in heifers (0.016 to 0.026) were lower than those in first-parity (0.017 to 0.142) and second-parity (0.026 to 0.115) cows. Genetic correlations for fertility in first- and second-parity cows were very high (>0.920), whereas those between heifers and lactating cows were moderate (0.348 to 0.709). The latter result indicates that fertility in heifers is a different trait than fertility in lactating cows, and hence it cannot be used as robust indicator of cow fertility. Heifer fertility was not related to production traits in first lactation (genetic correlations between −0.215 and 0.251). Peak milk yield exerted a moderate and unfavorable effect on the interval from parturition to first service (genetic correlations of 0.414 and 0.353 after first and second calving, respectively), and a low and unfavorable effect on other fertility traits (genetic correlations between −0.281 and 0.295). Infertility after first calving caused a strong elongation of the lactation, and LL was negatively correlated with fertility of cows after second calving, so that LL can itself be regarded as a measure of fertility. Lactation milk yield depends on both pMY and LL, and, as such, is a cause and consequence of (in)fertility.     

Short communication: Heritability of methane production and genetic correlations with milk yield and body weight in Holstein-Friesian dairy cows

Greenhouse gases originating from the dairy sector, including methane (CH₄), contribute to global warming. A possible strategy to reduce CH₄ production is to use genetic selection. This requires genetic parameters for CH₄ production and correlations with production traits. Data were available on 184 Holstein-Friesian cows. Methane production was measured in the milking robot during milking from December 2009 to April 2010. In total 2,456 observations for CH₄ production were available. Milk yield (MY) and body weight (BW) were obtained at every milking from November 2008 to October 2010. In total 4,567 observations for milk yield and 4,570 observations for BW were available. Restricted maximum likelihood, using random regression models, was used to analyze the data. Heritability (standard error given in parentheses) for CH₄ production ranged from 0.12 (0.16) to 0.45 (0.11), and genetic correlations with MY ranged from 0.49 (0.12) to 0.54 (0.26). The positive genetic correlation between CH₄ production and milk yield indicates that care needs to be taken when genetically selecting for lower CH₄ production, to avoid a decrease in MY at the animal level. However, this study shows that CH₄ production is moderately heritable and therefore progress through genetic selection is possible.     

Genetic parameters for buffalo milk yield and milk quality traits using Bayesian inference

The availability of accurate genetic parameters for important economic traits in milking buffaloes is critical for implementation of a genetic evaluation program. In the present study, heritabilities and genetic correlations for fat (FY305), protein (PY305), and milk (MY305) yields, milk fat (%F) and protein (%P) percentages, and SCS were estimated using Bayesian methodology. A total of 4,907 lactations from 1,985 cows were used. The (co)variance components were estimated using multiple-trait analysis by Bayesian inference method, applying an animal model, through Gibbs sampling. The model included the fixed effects of contemporary groups (herd-year and calving season), number of milking (2 levels), and age of cow at calving as (co)variable (quadratic and linear effect). The additive genetic, permanent environmental, and residual effects were included as random effects in the model. The posterior means of heritability distributions for MY305, FY305, PY305, %F, P%, and SCS were 0.22, 0.21, 0.23, 0.33, 0.39, and 0.26, respectively. The genetic correlation estimates ranged from −0.13 (between %P and SCS) to 0.94 (between MY305 and PY305). The permanent environmental correlation estimates ranged from −0.38 (between MY305 and %P) to 0.97 (between MY305 and PY305). Residual and phenotypic correlation estimates ranged from −0.26 (between PY305 and SCS) to 0.97 (between MY305 and PY305) and from −0.26 (between MY305 and SCS) to 0.97 (between MY305 and PY305), respectively. Milk yield, milk components, and milk somatic cells counts have enough genetic variation for selection purposes. The genetic correlation estimates suggest that milk components and milk somatic cell counts would be only slightly affected if increasing milk yield were the selection goal. Selecting to increase FY305 or PY305 will also increase MY305, %P, and %F.     

Breeding for reduced methane emission and feed-efficient Holstein cows: An international response

Selecting for lower methane (CH₄) emitting animals is one of the best approaches to reduce CH₄ given that genetic progress is permanent and cumulative over generations. As genetic selection requires a large number of animals with records and few countries actively record CH₄, combining data from different countries could help to expedite accurate genetic parameters for CH₄ traits and build a future genomic reference population. Additionally, if we want to include CH₄ in the breeding goal, it is important to know the genetic correlations of CH₄ traits with other economically important traits. Therefore, the aim of this study was first to estimate genetic parameters of 7 suggested methane traits, as well as genetic correlations between methane traits and production, maintenance, and efficiency traits using a multicountry database. The second aim was to estimate genetic correlations within parities and stages of lactation for CH₄. The third aim was to evaluate the expected response of economically important traits by including CH₄ traits in the breeding goal. A total of 15,320 methane production (MeP, g/d) records from 2,990 cows belonging to 4 countries (Canada, Australia, Switzerland, and Denmark) were analyzed. Records on dry matter intake (DMI), body weight (BW), body condition score, and milk yield (MY) were also available. Additional traits such as methane yield (MeY; g/kg DMI), methane intensity (MeI; g/kg energy-corrected milk), a genetic standardized methane production, and 3 definitions of residual methane production (g/d), residual feed intake, metabolic BW (MBW), BW change, and energy-corrected milk were calculated. The estimated heritability of MeP was 0.21, whereas heritability estimates for MeY and MeI were 0.30 and 0.38, and for the residual methane traits heritability ranged from 0.13 to 0.16. Genetic correlations between different methane traits were moderate to high (0.41 to 0.97). Genetic correlations between MeP and economically important traits ranged from 0.29 (MY) to 0.65 (BW and MBW), being 0.41 for DMI. Selection index calculations showed that residual methane had the most potential for inclusion in the breeding goal when compared with MeP, MeY, and MeI, as residual methane allows for selection of low methane emitting animals without compromising other economically important traits. Inclusion of residual feed intake in the breeding goal could further reduce methane, as the correlation with residual methane is moderate and elicits a favorable correlated response. Adding a negative economic value for methane could facilitate a substantial reduction in methane emissions while maintaining an increase in milk production.     

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire,Holstein, and Jersey dairy cattle

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.     

Genetic (co-)variation of methane emissions,efficiency, and production traits in Danish Holstein cattle along and across lactations

Methane emissions in ruminant livestock has become a hot topic, given the pressure to reduce greenhouse gas emissions drastically in the European Union over the next 10 to 30 yr. During the 2021 United Nations Climate Change conference, countries also made collective commitments to curb methane emissions by 2050. Genetic selection for low-methane-emitting animals, particularly dairy cows, is one possible strategy for mitigation. However, it is essential to understand how methane emissions in lactating animals vary along lactation and across lactations. This understanding is useful when making decisions for future phenotyping strategies, such as the frequency and duration of phenotyping within and across lactations. Therefore, the objectives of this study were to estimate (1) genetic parameters for 2 methane traits: methane concentration (MeC) and methane production (MeP) at 2 parity levels in Danish Holstein cows across the entire lactation using random regression models; (2) genetic correlations within and between methane traits across the entire lactation; and (3) genetic correlations between the methane traits and economically important traits throughout first lactation. Methane concentration (n = 19,639) records of 575 Danish Holstein cows from a research farm measured between 2013 and 2020 were available. Subsequently, CH₄ production in grams/day (MeP; n = 13,866) was calculated; MeP and MeC for first and second lactation (L1 and L2) were analyzed as separate traits: MeC_L1, MeP_L1, MeC_L2, and MeP_L2. Heritabilities, variance components, and genetic correlations within and between the 4 CH₄ traits were estimated using random regression models with Legendre polynomials. The additive genetic and permanent environmental effects were modeled using second-order Legendre polynomial for lactation weeks. Estimated heritabilities for MeP_L1 ranged between 0.11 and 0.49, for MeC_L1 between 0.10 and 0.28, for MeP_L2 between 0.14 and 0.36, and for MeC_L2 between 0.13 and 0.29. In general, heritability estimates of MeC traits were lower and more stable throughout lactation and were similar between lactations compared with MeP. Genetic correlations (within trait) at different lactation weeks were generally highly positive (0.7) for most of the first lactation, except for the correlation of early lactation (<10 wk) with late lactation (>40 wk) where the correlation was the lowest (<0.5). Genetic correlations between methane traits were moderate to highly correlated during early and mid lactation. Finally, MeP_L1 has stronger genetic correlations with energy-corrected milk and dry matter intake compared with MeC_L1. In conclusion, both traits are different along (and across) lactation(s) and they correlated differently with production, maintenance, and intake traits, which is important to consider when including one of them in a future breeding objective.     

Thin and fat cows,and the nonlinear genetic relationship between body condition score and fertility

Thin and fat cows are often credited for low fertility, but body condition score (BCS) has been traditionally treated as a linear trait when genetic correlations with reproductive performance have been estimated. The aims of this study were to assess genetic parameters for fertility, production, and body condition traits in the Brown Swiss population reared in the Alps (Bolzano-Bozen Province, Italy), and to investigate the possible nonlinearity among BCS and other traits by analyzing fat and thin cows. Records of BCS measured on a 5-point scale were preadjusted for year-season and days in milk at scoring, and were considered positive (1) for fat cows if they exceeded the value of 1 residual standard deviation or null (0) otherwise, whereas positive values for thin cows were imputed to records below −1 residual standard deviation. Fertility indicators measured on first- and second-parity cows were interval from parturition to first service, interval from first service to conception, interval from parturition to conception, number of inseminations to conception, conception at first service, and nonreturn rate at 56 d after first service. Production traits were peak milk yield, lactation milk yield, and lactation length. Data were from 1,413 herds and included 16,324 records of BCS, fertility, and production for first-parity, and 10,086 fertility records for second-parity cows. Animals calved from 2002 to 2007 and were progeny of 420 artificial insemination bulls. Genetic parameters for the aforementioned traits were obtained under univariate and bivariate threshold and censored linear sire models implemented in a Bayesian framework. Posterior means of heritabilities for BCS, fat cows, and thin cows were 0.141, 0.122, and 0.115, respectively. Genetic correlations of body condition traits with contemporary production were moderate to high and were between −0.556 and 0.623. Body condition score was moderately related to fertility in first (−0.280 to 0.497) and second (−0.392 to 0.248) lactation. The fat cow trait was scarcely related to fertility, particularly in first-parity cows (−0.203 to 0.281). Finally, the genetic relationships between thin cows and fertility were higher than those between BCS and fertility, both in first (−0.456 to 0.431) and second (−0.335 to 0.524) lactation. Body condition score can be considered a predictor of fertility, and it could be included in evaluation either as linear measure or as thin cow. In the second case, the genetic relationship with fertility was stronger, exacerbating the poorest body condition and considering the possible nonlinearity between fertility and energy reserves of the cow.     

Short communication: Heterosis and breed effects for milk production and udder health traits in crosses between Danish Holstein,Danish Red,and Danish Jersey

During the last decade, the use of systematic crossbreeding in dairy cattle herds has increased in several countries of the world. The aim of this study was to estimate the effect of breed proportion and heterosis on milk production traits and udder health traits in dairy cattle. The study was based on records on milk yield (MY), protein yield (PY), fat yield (FY), somatic cell score (SCS), and mastitis (MAST) from 73,695 first-lactation dairy cows in 130 Danish herds applying systematic crossbreeding programs. Around 45% of the cows were crosses between Danish Holstein (DH), Danish Red (DR), or Danish Jersey (DJ), and the remaining were purebred DH, DR, or DJ. The statistical model included the fixed effects of herd-year, calving month, and calving age and an effect representing the lactation status of the cow. In addition, the model included a regression on calving interval from first to second lactation, a regression on the proportion of DH, DR, and DJ genes, and a regression on the degree of heterozygosity between DH and DR, DH and DJ, and DR and DJ. Random effects were the genetic effect of the cow and a residual. The effect of breed proportions was estimated relatively to DH. For MY, a pure DR yielded 461 kg milk less than DH, whereas a pure DJ yielded 2,259 kg milk less than a pure DH. Compared with DH, PY was 41.7 kg less for DJ, whereas PY for DR was 4.0 kg less than for DH. For FY, a DR yielded 10.6 kg less than DH, whereas there was no significant effect of breed proportion between DJ and DH. A DR cow had lower SCS (0.13) than DH, whereas DJ had higher SCS (0.14) than DH. There was no significant effect of breed proportion on MAST between the 3 breeds. Heterosis was significant in all combinations of breeds for MY, FY, and PY. Heterosis for crosses between DH and DR was 257 kg (3.2%), 11.9 kg (3.2%), and 8.9 kg (3.2%) for MY, PY, and FY, respectively. Corresponding figures for crosses between DH and DJ were 314 kg (4.4%), 14.3 kg (4.4%), and 10.4 kg (4.0%), whereas heterosis between DR and DJ was 462 kg (6.7%), 19.6 kg (6.7%), and 13.9 kg (5.4%) for MY, PY, and FY, respectively. Heterosis was only significant for SCS in the crosses between DH and DR. Heterosis effects for MAST were nonsignificant for all the crosses. The results obtained in this study demonstrate that in first lactation cows, there is a positive effect of heterosis on milk production traits, but limited effect on udder health traits.     

Accuracy of genomic predictions in Gyr (Bos indicus) dairy cattle

Genomic selection may accelerate genetic progress in breeding programs of indicine breeds when compared with traditional selection methods. We present results of genomic predictions in Gyr (Bos indicus) dairy cattle of Brazil for milk yield (MY), fat yield (FY), protein yield (PY), and age at first calving using information from bulls and cows. Four different single nucleotide polymorphism (SNP) chips were studied. Additionally, the effect of the use of imputed data on genomic prediction accuracy was studied. A total of 474 bulls and 1,688 cows were genotyped with the Illumina BovineHD (HD; San Diego, CA) and BovineSNP50 (50K) chip, respectively. Genotypes of cows were imputed to HD using FImpute v2.2. After quality check of data, 496,606 markers remained. The HD markers present on the GeneSeek SGGP-20Ki (15,727; Lincoln, NE), 50K (22,152), and GeneSeek GGP-75Ki (65,018) were subset and used to assess the effect of lower SNP density on accuracy of prediction. Deregressed breeding values were used as pseudophenotypes for model training. Data were split into reference and validation to mimic a forward prediction scheme. The reference population consisted of animals whose birth year was ≤2004 and consisted of either only bulls (TR1) or a combination of bulls and dams (TR2), whereas the validation set consisted of younger bulls (born after 2004). Genomic BLUP was used to estimate genomic breeding values (GEBV) and reliability of GEBV (R²_PEV) was based on the prediction error variance approach. Reliability of GEBV ranged from ∼0.46 (FY and PY) to 0.56 (MY) with TR1 and from 0.51 (PY) to 0.65 (MY) with TR2. When averaged across all traits, R²_PEV were substantially higher (R²_PEV of TR1 = 0.50 and TR2 = 0.57) compared with reliabilities of parent averages (0.35) computed from pedigree data and based on diagonals of the coefficient matrix (prediction error variance approach). Reliability was similar for all the 4 marker panels using either TR1 or TR2, except that imputed HD cow data set led to an inflation of reliability. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information. A reduced panel of ∼15K markers resulted in reliabilities similar to using HD markers. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information.     

Variations in methane yield and microbial community profiles in the rumen of dairy cows as they pass through stages of first lactation

Considerable interest exists both from an environmental and economic perspective in reducing methane emissions from agriculture. In ruminants, CH₄ is produced by a complex community of microorganisms that is established in early life but can be influenced by external factors such as feed. Although CH₄ emissions were thought to be constant once an animal reached maturity, recent studies have shown that CH₄ yield significantly increases from early to late lactation in dairy cows. The aim of this study was to test the hypothesis that increases in CH₄ yield over the lactation cycle are related to changes in rumen microbial community structure. Nine cows were monitored throughout their first lactation cycle. Methane and dry matter intake were measured to calculate CH₄ per dry matter intake (CH₄ yield) and ruminal fluid was collected during early, mid, and late lactation. A significant difference in bacterial and archaeal community structure during early and late lactation was observed. Furthermore, when ruminal short-chain fatty acid concentrations were measured, the ratio of acetate and butyrate to propionate was significantly higher in late lactation compared with early lactation. Propionate concentrations were higher in cows with low CH₄ yield during late lactation, but no differences were observed in bacterial or archaeal community structures. Prevotella dominated the rumen of cows followed by Succinclasticum; Treponema, Fibrobacter, Ruminococcus, and Bifidobacterium were also in high abundance relative to other bacterial genera. In general, positive correlations were stronger between the most relatively abundant bacterial genera and acetate and butyrate concentrations in the cows with high CH₄ and weaker between these genera and propionate concentration. This study indicates that increased CH₄ yield in late lactation is reflected in significant changes in microbial community structure.     

Climate sensitivity of milk production traits and milk fatty acids in genotyped Holstein dairy cows

The aim of this study was the evaluation of climate sensitivity via genomic reaction norm models [i.e., to infer cow milk production and milk fatty acid (FA) responses on temperature-humidity index (THI) alterations]. Test-day milk traits were recorded between 2010 and 2016 from 5,257 first-lactation genotyped Holstein dairy cows. The cows were kept in 16 large-scale cooperator herds, being daughters of 344 genotyped sires. The longitudinal data consisted of 47,789 test-day records for the production traits milk yield (MY), fat yield (FY), and protein yield (PY), and of 20,742 test-day records for 6 FA including C16:0, C18:0, saturated fatty acids (SFA), unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA). After quality control of the genotypic data, 41,057 SNP markers remained for genomic analyses. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI. Genomic reaction norm models were applied to estimate genetic parameters in a single-step approach for production traits and FA in dependency of THI at different lactation stages, and to evaluate the model stability. In a first evaluation strategy (New_sire), all phenotypic records from daughters of genotyped sires born after 2010 were masked, to mimic a validation population. In the second strategy (New_env), only daughter records of the new sires recorded in the most extreme THI classes were masked, aiming at predicting sire genomic estimated breeding values (GEBV) under heat stress conditions. Model stability was the correlation between GEBV of the new sires in the reduced data set with respective GEBV estimated from all phenotypic data. Among all test-day production traits, PY responded as the most sensitive to heat stress. As observed for the remaining production traits, genetic variances were quite stable across THI, but genetic correlations between PY from temperate climates with PY from extreme THI classes dropped to 0.68. Genetic variances in dependency of THI were very similar for C16:0 and SFA, indicating marginal climatic sensitivity. In the early lactation stage, genetic variances for C18:0, MUFA, PUFA, and UFA were significantly larger in the extreme THI classes compared with the estimates under thermoneutral conditions. For C18:0 and MUFA, PUFA, and UFA in the middle THI classes, genetic correlations in same traits from the early and the later lactation stages were lower than 0.50, indicating strong days in milk influence. Interestingly, within lactation stages, genetic correlations for C18:0 and UFA recorded at low and high THI were quite large, indicating similar genetic mechanisms under stress conditions. The model stability was improved when applying the New_env instead of New_sire strategy, especially for FA in the first stage of lactation. Results indicate moderately accurate genomic predictions for milk traits in extreme THI classes when considering phenotypic data from a broad range of remaining THI. Phenotypically, thermal stress conditions contributed to an increase of UFA, suggesting value as a heat stress biomarker. Furthermore, the quite large genetic variances for UFA at high THI suggest the consideration of UFA in selection strategies for improved heat stress resistance.     

Predicting dry matter intake in Canadian Holstein dairy cattle using milk mid-infrared reflectance spectroscopy and other commonly available predictors via artificial neural networks

Dry matter intake (DMI) is a fundamental component of the animal's feed efficiency, but measuring DMI of individual cows is expensive. Mid-infrared reflectance spectroscopy (MIRS) on milk samples could be an inexpensive alternative to predict DMI. The objectives of this study were (1) to assess if milk MIRS data could improve DMI predictions of Canadian Holstein cows using artificial neural networks (ANN); (2) to investigate the ability of different ANN architectures to predict unobserved DMI; and (3) to validate the robustness of developed prediction models. A total of 7,398 milk samples from 509 dairy cows distributed over Canada, Denmark, and the United States were analyzed. Data from Denmark and the United States were used to increase the training data size and variability to improve the generalization of the prediction models over the lactation. For each milk spectra record, the corresponding weekly average DMI (kg/d), test-day milk yield (MY, kg/d), fat yield (FY, g/d), and protein yield (PY, g/d), metabolic body weight (MBW), age at calving, year of calving, season of calving, days in milk, lactation number, country, and herd were available. The weekly average DMI was predicted with various ANN architectures using 7 predictor sets, which were created by different combinations MY, FY, PY, MBW, and MIRS data. All predictor sets also included age of calving and days in milk. In addition, the classification effects of season of calving, country, and lactation number were included in all models. The explored ANN architectures consisted of 3 training algorithms (Bayesian regularization, Levenberg-Marquardt, and scaled conjugate gradient), 2 types of activation functions (hyperbolic tangent and linear), and from 1 to 10 neurons in hidden layers). In addition, partial least squares regression was also applied to predict the DMI. Models were compared using cross-validation based on leaving out 10% of records (validation A) and leaving out 10% of cows (validation B). Superior fitting statistics of models comprising MIRS information compared with the models fitting milk, fat and protein yields suggest that other unknown milk components may help explain variation in weekly average DMI. For instance, using MY, FY, PY, and MBW as predictor variables produced a predictive accuracy (r) ranging from 0.510 to 0.652 across ANN models and validation sets. Using MIRS together with MY, FY, PY, and MBW as predictors resulted in improved fitting (r = 0.679–0.777). Including MIRS data improved the weekly average DMI prediction of Canadian Holstein cows, but it seems that MIRS predicts DMI mostly through its association with milk production traits and its utility to estimate a measure of feed efficiency that accounts for the level of production, such as residual feed intake, might be limited and needs further investigation. The better predictive ability of nonlinear ANN compared with linear ANN and partial least squares regression indicated possible nonlinear relationships between weekly average DMI and the predictor variables. In general, ANN using Bayesian regularization and scaled conjugate gradient training algorithms yielded slightly better weekly average DMI predictions compared with ANN using the Levenberg-Marquardt training algorithm.     

Genetic parameters for predicted methane production and potential for reducing enteric emissions through genomic selection

Mitigation of enteric methane (CH₄) emission in ruminants has become an important area of research because accumulation of CH₄ is linked to global warming. Nutritional and microbial opportunities to reduce CH₄ emissions have been extensively researched, but little is known about using natural variation to breed animals with lower CH₄ yield. Measuring CH₄ emission rates directly from animals is difficult and hinders direct selection on reduced CH₄ emission. However, improvements can be made through selection on associated traits (e.g., residual feed intake, RFI) or through selection on CH₄ predicted from feed intake and diet composition. The objective was to establish phenotypic and genetic variation in predicted CH₄ output, and to determine the potential of genetics to reduce methane emissions in dairy cattle. Experimental data were used and records on daily feed intake, weekly body weights, and weekly milk production were available from 548 heifers. Residual feed intake (MJ/d) is the difference between net energy intake and calculated net energy requirements for maintenance as a function of body weight and for fat- and protein-corrected milk production. Predicted methane emission (PME; g/d) is 6% of gross energy intake (Intergovernmental Panel on Climate Change methodology) corrected for energy content of methane (55.65 kJ/g). The estimated heritabilities for PME and RFI were 0.35 and 0.40, respectively. The positive genetic correlation between RFI and PME indicated that cows with lower RFI have lower PME (estimates ranging from 0.18 to 0.84). Hence, it is possible to decrease the methane production of a cow by selecting more-efficient cows, and the genetic variation suggests that reductions in the order of 11 to 26% in 10 yr are theoretically possible, and could be even higher in a genomic selection program. However, several uncertainties are discussed; for example, the lack of true methane measurements (and the key assumption that methane produced per unit feed is not affected by RFI level), as well as the limitations of predicting the biological consequences of selection. To overcome these limitations, an international effort is required to bring together data on feed intake and methane emissions of dairy cows.