Advertisement
Research| Volume 102, ISSUE 12, P11207-11216, December 2019

Genetic parameters of endocrine fertility traits based on in-line milk progesterone profiles in Swedish Red and Holstein dairy cows

Open ArchivePublished:October 09, 2019DOI:https://doi.org/10.3168/jds.2019-16691

      ABSTRACT

      Evaluating fertility traits based on endocrine progesterone profiles is becoming a promising option to improve dairy cow fertility. Several studies have been conducted on endocrine fertility traits, mainly in the Holstein breed. In this study, focusing also on the Swedish Red (SR) breed, genetic parameters were estimated for classical and endocrine fertility traits, the latter based on in-line milk progesterone records obtained for 14 Swedish herds using DeLaval Herd Navigator (DeLaval International, Tumba, Sweden). A total of 210,403 observations from 3,437 lactations of 1,107 SR and 1,538 Holstein cows were used. Mixed linear animal models were used for estimation of genetic parameters. Least squares means analysis showed that Holstein cows had a 2.5-d-shorter interval from calving to commencement of luteal activity (C-LA) and longer length of first inter-ovulatory interval (IOI) than SR cows. The highest mean interval for C-LA, IOI, and first luteal phase length (LPL) was observed in the fourth parity. The incidence of short (<18 d), normal, (18–24 d), and long (>24 d) IOI was 29.3, 40.7, and 30%, respectively. Genetic analysis indicated moderate heritability (h2) for C-LA (h2 = 0.24), luteal activity during the first 60 d in milk (LA60, h2 = 0.15), proportion of samples with luteal activity (PLA, h2 = 0.13), and calving to first heat (CFH, h2 = 0.18), and low heritability estimates for LPL (h2 = 0.08) and IOI (h2 = 0.03) in the combined data set for both breeds. Similar heritability estimates were obtained for each breed separately except for IOI and LPL in SR cows, for which heritability was estimated to be zero. Swedish Red cows had 0.01 to 0.06 higher heritability estimates for C-LA, LA60, and PLA than did Holstein cows. Calving interval had moderate heritability among the classical traits for Holstein and the combined data set, but h2 was zero for SR. Commencement of luteal activity had a strong genetic correlation with LA60 (mean ± SE; −0.88 ± 0.06), PLA (−0.72 ± 0.11), and CFH (0.90 ± 0.04). Similarly, CFH had a strong genetic correlation with IOI (0.98 ± 0.20). Number of inseminations per series showed a weak genetic correlation with all endocrine traits except IOI. Overall, endocrine traits had higher heritability estimates than classical traits in both breeds, and may have a better potential to explain the actual reproductive status of dairy cows than classical traits. This might favor inclusion of some endocrine fertility traits—especially those related to commencement of luteal activity—as selection criteria and breeding goal traits if recording becomes more common in herds. Further studies on genetic and genomic evaluations for endocrine fertility traits may help to provide firm conclusions. A prerequisite is that the data from automatic devices be made available to recording and breeding organizations in the future and included in a central database.

      Key words

      INTRODUCTION

      Dairy cows have suffered a decline in fertility, reflected by an extended anestrous period, less pronounced estrous behavior (
      • Leroy J.L.M.R.
      • Opsomer G.
      • Van Soom A.
      • Goovaerts I.G.F.
      • Bols P.E.J.
      Reduced fertility in high-yielding dairy cows: are the oocyte and embryo in danger? Part I. The importance of negative energy balance and altered corpus luteum function to the reduction of oocyte and embryo quality in high-yielding dairy cows.
      ), and increased prevalence of atypical progesterone (P4) profiles (
      • Royal M.D.
      • Darwash A.O.
      • Flint A.P.F.
      • Webb R.
      • Woolliams J.A.
      • Lamming G.E.
      Declining fertility in dairy cattle: Changes in traditional and endocrine parameters of fertility.
      ;
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ). Understanding the underlying biology that contributes to low fertility, developing strategies to improve fertility, and utilizing automated tools to reveal cow fertility status remain major challenges in research on dairy cattle (
      • Walsh S.W.
      • Williams E.J.
      • Evans A.C.O.
      A review of the causes of poor fertility in high milk producing dairy cows.
      ). Most fertility improvement interventions to date have focused on classical traits based on insemination records. However, these traits are highly influenced by environmental factors and management decisions (
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). As an alternative, the use of endocrine traits, based on in-line milk measurements, to improve dairy cow fertility has been suggested. This is because they are less affected by environment and management decisions (
      • Royal M.
      • Mann G.E.
      • Flint A.P.F.
      Strategies for reversing the trend towards subfertility in dairy cattle.
      ;
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ) and provide objective and accurate measurements of ovarian activity (
      • Petersson K.J.
      • Gustafsson H.
      • Strandberg E.
      • Berglund B.
      Atypical progesterone profiles and fertility in Swedish dairy cows.
      ). Endocrine traits better reflect the physiology and reproductive status of the cow, which often leads to higher heritability estimates compared with classical traits (
      • Friggens N.C.
      • Chagunda M.G.G.
      Prediction of the reproductive status of cattle on the basis of milk progesterone measures: Model description.
      ;
      • Martin O.
      • Friggens N.C.
      • Dupont J.
      • Salvetti P.
      • Freret S.
      • Rame C.
      • Elis S.
      • Gatien J.
      • Disenhaus C.
      • Blanc F.
      Data-derived reference profiles with core presentation of progesterone, estradiol, LH, and FSH dynamics during the bovine estrous cycle.
      ;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). Recent studies on endocrine fertility traits show encouraging heritability and genetic correlation estimates and these traits can be useful indicators of fertility to consider in the selection criteria and breeding goal (
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ;
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ). Therefore, the aim of this study was to perform a comprehensive evaluation of endocrine fertility traits in the 2 most common breeds in Swedish dairy herds, using progesterone profiles derived from in-line milk progesterone measurements.
      In-line milk progesterone profiles are generated by the automated management system Herd Navigator (DeLaval International, Tumba, Sweden), a decision support tool developed for monitoring energy balance, mastitis, and fertility in dairy cows (
      • Friggens N.C.
      • Chagunda M.G.G.
      Prediction of the reproductive status of cattle on the basis of milk progesterone measures: Model description.
      ;
      • DeLaval International
      Herd Navigator and reproduction management.
      ). It enables monitoring of fertility status in dairy cows and reduces time-consuming visual heat detection (
      • Saint-Dizier M.
      • Chastant-Maillard S.
      Towards an automated detection of oestrus in dairy cattle.
      ;
      • Adriaens I.
      • Huybrechts T.
      • Geerinckx K.
      • Daems D.
      • Lammertyn J.
      • De Ketelaere B.
      • Saeys W.
      • Aernouts B.
      Mathematical characterization of the milk progesterone profile as a leg up to individualised monitoring of reproduction status in dairy cows.
      ). Moreover, it helps to identify the onset and end of luteal phases and to determine estrus and pregnancy status (
      • Adriaens I.
      • Saeys W.
      • Huybrechts T.
      • Lamberigts C.
      • François L.
      • Geerinckx K.
      • Leroy J.
      • Ketelaere B.D.
      • Aernouts B.
      A novel system for on-farm fertility monitoring based on milk progesterone.
      ;
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ). Overall, this tool can enable timely identification of a cow's reproductive status, helping to minimize fertility-related losses on dairy farms (e.g.,
      • Adriaens I.
      • Saeys W.
      • Huybrechts T.
      • Lamberigts C.
      • François L.
      • Geerinckx K.
      • Leroy J.
      • Ketelaere B.D.
      • Aernouts B.
      A novel system for on-farm fertility monitoring based on milk progesterone.
      ). In recent years, Herd Navigator has been installed on several dairy farms in Sweden. This provides the opportunity to evaluate endocrine fertility traits of Swedish Red cows in a large data set. Apart from a recent study in Finland (
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ), most progesterone-based fertility studies to date have focused on Holstein cows (
      • Berry D.P.
      • Bastiaansen J.W.M.
      • Veerkamp R.F.
      • Wijga S.
      • Wall E.
      • Berglund B.
      • Calus M.P.L.
      Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
      ;
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ;
      • Adriaens I.
      • Huybrechts T.
      • Geerinckx K.
      • Daems D.
      • Lammertyn J.
      • De Ketelaere B.
      • Saeys W.
      • Aernouts B.
      Mathematical characterization of the milk progesterone profile as a leg up to individualised monitoring of reproduction status in dairy cows.
      ,
      • Adriaens I.
      • Saeys W.
      • Huybrechts T.
      • Lamberigts C.
      • François L.
      • Geerinckx K.
      • Leroy J.
      • Ketelaere B.D.
      • Aernouts B.
      A novel system for on-farm fertility monitoring based on milk progesterone.
      ;
      • Bruinjé T.C.
      • Colazo M.G.
      • Ribeiro E.S.
      • Gobikrushanth M.D.
      • Ambrose J.
      Using in-line milk progesterone data to characterize parameters of luteal activity and their association with fertility in Holstein cows.
      ). In this study, we estimated genetic parameters for classical and endocrine traits, the latter based on in-line milk progesterone records, in both Swedish Red and Holstein dairy cows.

      MATERIALS AND METHODS

      A total of 14 Swedish dairy herds comprising Swedish Red (SR) and Holstein cows were included in the study. In-line milk P4 samples were automatically collected and analyzed by Herd Navigator. In the initial data set, 83.2% of the records were the only sample taken on that day; however, some cows had more than one P4 value for a given day. We used records from October 2015 to October 2017. The endocrine traits studied were interval from calving to commencement of luteal activity (C-LA), length of first inter-ovulatory interval (IOI), length of inter-luteal interval (ILI), first luteal phase length (LPL), luteal activity during the first 60 DIM (LA60; categorical 0 or 1), proportion of samples with luteal activity during the first 60 DIM (PLA; Figure 1), and interval from calving to first heat (CFH) and the occurrence of atypical profiles. The classical traits included calving interval (CInt), interval from calving to first insemination (CFI), interval from calving to last insemination (CLI), interval from first to last insemination (FLI), number of inseminations per series (NINS), and the mixed (endocrine and classical) trait interval from commencement of luteal activity to first service (CLAFS). Endocrine fertility traits were defined from progesterone profile data provided by Lattec I/S (Hillerød, Denmark) and based on a recent report by
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      . A progesterone concentration of 5 ng/mL was used as the threshold to define progesterone profiles. Pedigree, calving, and insemination data used to assess classical fertility traits were obtained from Växa Sverige (Stockholm, Sweden).
      Figure thumbnail gr1
      Figure 1Milk progesterone curve for endocrine fertility traits, generated from a cow used in the study. C-LA = commencement of luteal activity; LPL = first luteal phase length; IOI = length of first inter-ovulatory interval; ILI = length of first inter-luteal interval; CLAFS = interval from commencement of luteal activity to first service. A threshold of 5 ng/mL was used to determine whether the cow had luteal activity.
      Individual progesterone profiles were plotted against DIM using SAS software (ver. 9.2;
      • SAS Institute
      SAS/STAT® 12.1 User's Guide.
      ; Figure 1). The progesterone profiles were manually classified into normal and atypical (delayed cyclicity: P4 concentrations below the predefined P4 threshold for more than 56 d postpartum; prolonged luteal phase: high P4 levels >5 ng/mL for at least 20 d; cessation of cyclicity: cyclicity interrupted for at least 14 d with P4 concentrations below the P4 threshold value) progesterone profiles. Details of the atypical profile classification were reported by
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      . The incidence of IOI was classified as short (IOI <18 d), normal (18 ≤ IOI ≤24 d), or long (IOI >24 d) cycle (
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ).

      Filtering Criteria

      The data were filtered at the herd level and the lactation level. Herd-level filtering criteria were as follows: (1) lactation was retained if interval from calving to start of herd P4 recording was ≤20 d; and (2) lactation was retained if interval from calving to end of herd recording was ≥60 d. Lactation-level filtering criteria were as follows: (1) P4 samples taken before 20 DIM in lactation were excluded; and (2) lactations in which P4 recording started after 35 DIM were excluded. The numbers of samples, cows, and lactations retained after applying the filtering criteria are shown in Table 1.
      Table 1Number (N) of progesterone (P4) samples, lactations, and cows retained after applying filtering criteria at the herd and lactation levels
      CriteriaP4 samplesLactationCows
      N%N%N%
      Start data286,3211005,3711003,653100
      Herd level
      Herd-level filtering excluded lactations with interval from calving to start of P4 recording >20 d and lactations with interval from calving to end of P4 recording <60 d.
      243,072854,15477.33,09484.7
      Lactation level
      Lactation-level filtering excluded P4 records taken before 20 DIM and lactations with start of P4 recording >35 DIM.
      210,403733,437642,64572.4
      1 Herd-level filtering excluded lactations with interval from calving to start of P4 recording >20 d and lactations with interval from calving to end of P4 recording <60 d.
      2 Lactation-level filtering excluded P4 records taken before 20 DIM and lactations with start of P4 recording >35 DIM.
      For specific traits at the cow lactation level, the following restrictions were applied if the cow's lactation had no P4 records for more than 7 d: all traits were retained if the gap occurred one or more days before C-LA and the gap duration was less than 15 d, or if the gap occurred more than 24 d after C-LA. To handle herd-specific gaps, the following restrictions were applied to that herd: all traits were retained if a gap of at least 3 d occurred one or more days before C-LA. The following restrictions were applied if a gap of 3 or more days occurred after C-LA: (1) LPL was retained if the gap occurred more than 24 d after C-LA; (2) ILI was retained if the gap occurred at least 1 d before the start of inter-luteal interval; (3) IOI was retained if the gap occurred at least 1 d before the start of next luteal cycle. Insemination events that occurred before the start of recorded luteal activity were excluded. For CInt, CLI, and FLI, records above 730, 300, and 300 d, respectively, were excluded.

      Statistical Analysis

      Effects of breed, parity, and season on endocrine and classical fertility traits were estimated using a general linear model in SAS software (ver. 9.2;
      • SAS Institute
      SAS/STAT® 12.1 User's Guide.
      ):
      yijkl = μ + Bi + Pj + Sk + eijkl,
      [1]


      where yijkl is the analyzed trait; μ is the overall mean; Bi is the fixed effect of the ith breed (SR or Holstein); Pj is the fixed effect of jth parity (lactation) (j = 1, …, 5; parities >5 were merged into the fifth parity); Sk is the fixed effect of kth calving season (k = March–May, June–August, September–November, and December–February); and eijkl is the random error. Because 2 herds had only SH cows, it was not possible to estimate a herd effect free of the breed effect. Therefore, the herd effect was not included. Heritabilities and genetic correlations were estimated using a mixed linear animal model in the DMU package (
      • Madsen P.
      • Jensen J.
      A user's guide to DMU. A package for analyzing multivariate mixed models. Version 6, release 5.2.
      ). Breed and parity were used as fixed effects in the model for combined analysis of SR and Holstein cows, whereas parity was the only fixed effect used in separate analysis of the breeds. Herd-year-season (hys) combinations and cow were fitted as random effects in the model. Univariate analysis was used for all classical and endocrine fertility traits. Bivariate analysis was used to estimate the correlation between endocrine and classical fertility traits. All traits except LA60, PLA, and NINS were (natural) log-transformed. The model used for both univariate and bivariate analyses of the joint data set of Swedish Red and Holstein cows was
      yijklm = μ + Bi + Pj + hysk + al + eijklm,
      [2]


      where yijklm is the analyzed trait; μ, Bi, and Pj are as in model [1]; hysk = is the random effect of herd-year-season combination k (14 herds, years 2015, 2016, and 2017, and the 4 seasons described in model [1]), ~IND(0,σHYS2); where IND is independent, normally distributed, σHYS2 is the random additive herd year season variance, al is the random additive genetic effect of cow l with var (a), ~ND(0,Aσa2), where ND is normally distributed, σa2 is the additive genetic variance, A is the additive relationship matrix using pedigree information traced back 5 generations for 33,120 animals; and eijklm is the random residual with var (e) ~IND(0,σe2), where σe2 is the residual variance. Heritability was estimated as σˆa2/(σˆa2+σˆe2) and standard error (SE) as in
      • Edwards S.M.
      Standard errors of heritability estimates.
      .

      RESULTS

      Descriptive Analysis of Endocrine and Classical Fertility Traits

      The endocrine and classical traits of the herds in the study sample are summarized in Table 2. After applying the filtering criteria, the lowest number of observations remained for ILI (n = 1,412) among the endocrine traits and for CInt (n = 1,084) among the classical traits. The number of observations for atypical profiles was low, particularly because they were generated only from the first estrous cycle. The average values obtained for endocrine traits were higher in Holstein than in SR cows, except for C-LA and CFH. This was also the case for classical traits. The proportion of samples with luteal activity (PLA) was 6 percentage points lower for SR cows than for Holstein cows. The incidence of short, normal, and long IOI was 29.3, 40.7, and 30%, respectively. From the total of 2,793 P4 profiles generated from 2,283 cows, 70% of the lactations obtained from 75.1% of cows (data not shown) had a normal profile, and 19.2% of the cows had atypical profiles that constituted 16.5% (n = 460) of the total P4 profiles. The remaining 13.5% of the profiles obtained from 5.7% of the cows were categorized as neither normal nor atypical, and might be explained by early pregnancy. Prolonged luteal phase had an incidence of 12.9% (from total LPL) and cessation of cyclicity had an incidence of 4.5% (from total ILI).
      Table 2Summary statistics [number (N), mean, and SD] on endocrine and classical fertility traits in the 14 Swedish Herd Navigator (DeLaval International, Tumba, Sweden) dairy herds studied
      Trait
      C-LA = commencement of luteal activity; normal profile = a first increase in progesterone (P4) concentration before d 56 postpartum, followed by 2 wk of high P4 concentration and then 1 wk of low P4 concentration; LPL = first luteal phase length; ILI = length of first inter-luteal interval; IOI = length of first inter-ovulatory interval; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; CLAFS = commencement of luteal activity to first service; CFH = calving to first heat; Delayed = delayed cyclicity P4 level <5 ng/mL for >56 d postpartum; Cessation = cessation of cyclicity (cyclicity interrupted for at ≥14 d with P4 level <5 ng/mL); Prolonged = prolonged cyclicity (P4 level >5 ng/mL for ≥20 d); CInt = calving interval; CFI = calving to first insemination; CLI = calving to last insemination; FLI = first to last insemination; NINS = number of inseminations per series.
      Swedish Red (SR)HolsteinSR and Holstein
      NMeanSDNMeanSDMinimumMaximum
      C-LA (d)
       All1,17335.217.21,62032.717.120234
       Delayed8480.925.511978.125.857234
       Normal86931.709.51,09429.409.22055
      LPL (d)
       All61710.914.292513.516.21252
       Prolonged5342.533.514939.727.420252
      ILI (d)
       All5567.785.78568.06.1182
       Cessation2527.414. 53928.915.11482
      IOI (d)62321.812.797524.615.71219
      LA60 (0, 1)1,1730.940.241,6200.940.2401
      PLA (%)1,3380.590.191,8710.650.200.041
      CLAFS (d)64753.129.81,01062.933.84403
      CFH (d)1,37646.621.51,94045.819.421326
      CInt (d)46038046.362438850.8265705
      CFI (d)1,23387.329.61,76691.927.234297
      CLI (d)1,110116.546.81,612133.353.239299
      FLI (d)1,11329.340.71,63743.051.60275
      NINS (N)1,2341.901.11,7722.201.418
      1 C-LA = commencement of luteal activity; normal profile = a first increase in progesterone (P4) concentration before d 56 postpartum, followed by 2 wk of high P4 concentration and then 1 wk of low P4 concentration; LPL = first luteal phase length; ILI = length of first inter-luteal interval; IOI = length of first inter-ovulatory interval; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; CLAFS = commencement of luteal activity to first service; CFH = calving to first heat; Delayed = delayed cyclicity P4 level <5 ng/mL for >56 d postpartum; Cessation = cessation of cyclicity (cyclicity interrupted for at ≥14 d with P4 level <5 ng/mL); Prolonged = prolonged cyclicity (P4 level >5 ng/mL for ≥20 d); CInt = calving interval; CFI = calving to first insemination; CLI = calving to last insemination; FLI = first to last insemination; NINS = number of inseminations per series.
      Holstein cows started luteal activity 2.5 d earlier and had a luteal phase 2.4 d longer and IOI 3.0 d longer than SR cows (Table 3). Cows in their first parity started luteal activity 4.5 d earlier than cows in their fourth parity. Similarly, cows in their first parity had first estrus 6.8 d earlier, LPL 3.5 d shorter, ILI 2.1 d shorter, and ovulatory interval 5.3 d shorter than cows in fourth parity. Calving season had a significant influence on C-LA; the highest estimate was obtained in spring and the lowest in summer (Table 3). Similarly, calving season influenced LA60, PLA, NINS, and CInt. Cows calving during the summer had a C-LA 4 d earlier than cows calving during spring. Similarly, cows calving during summer had an LA60 6 percentage points higher and PLA 4 percentage points higher than spring calvers. However, we observed no association of season (P > 0.05) with LPL, ILI, IOI, or LA60. Number of inseminations per series (NINS) was 30 percentage points higher for Holstein than for SR cows and 27 percentage points higher in autumn than in spring.
      Table 3Least squares means and standard errors (LSM ± SE) of endocrine and classical fertility traits
      C-LA = commencement of luteal activity; LPL = first luteal phase length; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; ILI = length of first inter-luteal interval; IOI = length of first inter-ovulatory interval; NINS = number of inseminations per series; CInt = calving interval; N = number of observations.
      EffectC-LA (d)LPL (d)LA60PLA
      NLSM ± SENLSM ± SENLSM ± SENLSM ± SE
      Breed
      P < 0.01
      P < 0.01
      NS
      P < 0.01
       Swedish Red1,17336.2 ± 0.56
      LSM with different superscript letters differ significantly.
      63111.3 ± 0.57
      LSM with different superscript letters differ significantly.
      1,1730.93 ± 0.011,3380.58 ± 0.00
      LSM with different superscript letters differ significantly.
       Holstein1,62033.7 ± 0.49
      LSM with different superscript letters differ significantly.
      94514.1 ± 0.57
      LSM with different superscript letters differ significantly.
      1,6200.93 ± 0.001,8710.65 ± 0.01
      LSM with different superscript letters differ significantly.
      Parity
      P < 0.001.
      P < 0.05
      NS
      P < 0.05
       11,07532.7 ± 0.52
      LSM with different superscript letters differ significantly.
      51210.2 ± 0.69
      LSM with different superscript letters differ significantly.
      1,0750.94 ± 0.011,2260.63 ± 0.01
      LSM with different superscript letters differ significantly.
       282433.9 ± 0.59
      LSM with different superscript letters differ significantly.
      45612.6 ± 0.73
      LSM with different superscript letters differ significantly.
      8240.94 ± 0.019560.61 ± 0.01
      LSM with different superscript letters differ significantly.
       347234.1 ± 0.78
      LSM with different superscript letters differ significantly.
      33613.3 ± 0.84
      LSM with different superscript letters differ significantly.
      4720.95 ± 0.015530.61 ± 0.02
      LSM with different superscript letters differ significantly.
       427337.1 ± 1.0
      LSM with different superscript letters differ significantly.
      17714.0 ± 1.2
      LSM with different superscript letters differ significantly.
      2730.91 ± 0.013010.62 ± 0.01
      LSM with different superscript letters differ significantly.
       514936.8 ± 1.3
      LSM with different superscript letters differ significantly.
      9513.3 ± 1.6
      LSM with different superscript letters differ significantly.
      1490.92 ± 0.021730.59 ± 0.01
      LSM with different superscript letters differ significantly.
      Season
      Season = calving season; spring = March–May; summer = June–August; autumn = September–November; winter = December–February.
      P < 0.001.
      NS
      P < 0.001.
      P < 0.01
       Winter78235.4 ± 0.66
      LSM with different superscript letters differ significantly.
      41912.90 ± 0.797820.93 ± 0.00
      LSM with different superscript letters differ significantly.
      8880.61 ± 0.01
      LSM with different superscript letters differ significantly.
       Spring67737.8 ± 0.70
      LSM with different superscript letters differ significantly.
      37413.77 ± 0.846770.90 ± 0.00
      LSM with different superscript letters differ significantly.
      6920.60 ± 0.01
      LSM with different superscript letters differ significantly.
       Summer65832.9 ± 0.71
      LSM with different superscript letters differ significantly.
      43212.22 ± 0.806580.96 ± 0.01
      LSM with different superscript letters differ significantly.
      7520.64 ± 0.01
      LSM with different superscript letters differ significantly.
       Autumn67633.6 ± 0.71
      LSM with different superscript letters differ significantly.
      35111.87 ± 0.886760.94 ± 0.01
      LSM with different superscript letters differ significantly.
      8770.61 ± 0.01
      LSM with different superscript letters differ significantly.
      ILI (d)IOI (d)NINS (n)CInt (d)
      NLSM ± SENLSM ± SENLSM ± SENLSM ± SE
      BreedNS
      P < 0.001.
      P < 0.001.
      P < 0.01
       Swedish Red6607.98 ± 0.2862322.3 ± 0.651,2341.90 ± 0.04
      LSM with different superscript letters differ significantly.
      460382 ± 2.8
      LSM with different superscript letters differ significantly.
       Holstein1,0038.34 ± 0.2397425.2 ± 0.531,7722.20 ± 0.04
      LSM with different superscript letters differ significantly.
      624390 ± 2.5
      LSM with different superscript letters differ significantly.
      Parity
      P < 0.001.
      P < 0.001.
      NS
      P < 0.01
       14576.67 ± 0.28
      LSM with different superscript letters differ significantly.
      53620.8 ± 0.64
      LSM with different superscript letters differ significantly.
      1,1842.00 ± 0.04496378 ± 2.2
      LSM with different superscript letters differ significantly.
       24298.19 ± 0.29
      LSM with different superscript letters differ significantly.
      46724.2 ± 0.68
      LSM with different superscript letters differ significantly.
      9082.10 ± 0.04305385 ± 2.8
      LSM with different superscript letters differ significantly.
       32968.58 ± 0.34
      LSM with different superscript letters differ significantly.
      33524.0 ± 0.80
      LSM with different superscript letters differ significantly.
      5242.04 ± 0.06161389 ± 3.8
      LSM with different superscript letters differ significantly.
       41519.20 ± 0.48
      LSM with different superscript letters differ significantly.
      16626.1 ± 1.1
      LSM with different superscript letters differ significantly.
      2582.08 ± 0.0882390 ± 5.4
      LSM with different superscript letters differ significantly.
       5798.57 ± 0.67
      LSM with different superscript letters differ significantly.
      9323.6 ± 1.5
      LSM with different superscript letters differ significantly.
      1322.03 ± 0.1140389 ± 7.7
      LSM with different superscript letters differ significantly.
      Season
      Season = calving season; spring = March–May; summer = June–August; autumn = September–November; winter = December–February.
      NSNS
      P < 0.001.
      P < 0.01
       Winter4398.23 ± 0.3242524.1 ± 0.768321.99 ± 0.05
      LSM with different superscript letters differ significantly.
      304394 ± 3.1
      LSM with different superscript letters differ significantly.
       Spring4098.61 ± 0.3439024.3 ± 0.786621.95 ± 0.05
      LSM with different superscript letters differ significantly.
      234393 ± 3.6
      LSM with different superscript letters differ significantly.
       Summer4517.85 ± 0.3342422.8 ± 0.766862.05 ± 0.05
      LSM with different superscript letters differ significantly.
      213371 ± 3.7
      LSM with different superscript letters differ significantly.
       Autumn3647.94 ± 0.3535823.9 ± 0.828262.22 ± 0.05
      LSM with different superscript letters differ significantly.
      333386 ± 3.1
      LSM with different superscript letters differ significantly.
      a–c LSM with different superscript letters differ significantly.
      1 C-LA = commencement of luteal activity; LPL = first luteal phase length; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; ILI = length of first inter-luteal interval; IOI = length of first inter-ovulatory interval; NINS = number of inseminations per series; CInt = calving interval; N = number of observations.
      2 Season = calving season; spring = March–May; summer = June–August; autumn = September–November; winter = December–February.
      * P < 0.05
      ** P < 0.01
      *** P < 0.001.

      Heritability Estimates

      Among the endocrine fertility traits, C-LA had the highest estimates of heritability 0.24 to 0.26, whereas LA60, PLA, and CFH had moderate heritability estimates for SR, Holstein, and the combined data set (Table 4). The heritabilities for LA60 and PLA in SR cows were 0.20 and 0.14, respectively, which was higher than the heritability estimated for Holstein (0.14 and 0.12, respectively), albeit not significantly higher. However, unlike in Holstein, zero heritabilities were obtained for IOI, LPL, and ILI among the endocrine fertility traits, for CLAFS, and for the classical fertility traits CFI, CLI, and CInt in SR cows. Prolonged luteal phase and delayed cyclicity had moderate heritability estimates, whereas cessation of cyclicity had zero heritability. Among the classical traits, CInt had moderate heritability (0.15) in the combined data set and in Holstein cows (0.16). The estimated heritability for NINS was 0.02 to 0.06.
      Table 4Estimates of heritability (h2), standard error (SE), and additive genetic variance σ2a obtained for endocrine and classical fertility traits in Holstein and Swedish Red (SR) dairy cows
      Trait
      C-LA = commencement of luteal activity (all C-LA profiles); LPL = first luteal phase length (all LPL profiles); ILI = length of first inter-luteal interval (all ILI profiles); IOI = length of first inter-ovulatory interval; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; CLAFS = commencement of luteal activity to first service; CFH = calving to first heat; CFI = calving to first insemination in days; CLI = calving to last insemination in days; FLI = first to last insemination in days; NINS = number of inseminations per series; CInt = calving interval in days.
      Holstein and SRHolsteinSR
      h2SEσ2ah2SEσ2ah2SEσ2a
      C-LA0.240.040.070.250.230.0070.260.060.007
      LPL0.080.040.010.120.060.030.000.060.00
      ILI0.000.030.000.000.040.000.000.060.00
      IOI0.030.040.0020.090.060.0050.000.060.00
      LA600.150.040.010.140.040.010.200.060.01
      PLA0.130.030.0050.120.040.0050.140.050.005
      CLAFS0.030.040.0010.070.060.0030.000.060.00
      CFH0.180.050.0040.200.040.0050.140.050.003
      CFI0.030.0010.000.040.000.000.000.000.00
      CLI0.060.0010.0020.070.000.0020.000.030.00
      FLI0.040.020.030.000.050.000.020.030.01
      NINS0.020.060.080.060.030.100.040.040.05
      CInt0.150.080.000.160.120.000.000.090.00
      1 C-LA = commencement of luteal activity (all C-LA profiles); LPL = first luteal phase length (all LPL profiles); ILI = length of first inter-luteal interval (all ILI profiles); IOI = length of first inter-ovulatory interval; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; CLAFS = commencement of luteal activity to first service; CFH = calving to first heat; CFI = calving to first insemination in days; CLI = calving to last insemination in days; FLI = first to last insemination in days; NINS = number of inseminations per series; CInt = calving interval in days.

      Genetic Correlations

      Due to zero or very low heritability estimates and small sample size, genetic correlation analysis was carried out only for selected traits in the combined data set for both breeds (Table 5). The results revealed that C-LA had strong genetic correlations with CFH (0.90), LA60 (−0.88), and PLA (−0.72), and CLAFS was strongly correlated with PLA (0.99). Similarly, CFH had strong genetic correlations with IOI (0.98) and LA60 (0.88). In contrast, NINS had weak genetic correlations with all endocrine traits except IOI (0.61). The standard error for most variables was high (range 0.04–0.92).
      Table 5Estimates of genetic correlations with endocrine and classical fertility traits (above diagonal) and standard errors (below diagonal)
      Trait
      C-LA = commencement of luteal activity (all C-LA profiles); LPL = first luteal phase length (all LPL profiles); IOI = length of first inter-ovulatory interval; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; CLAFS = commencement of luteal activity to first service; CFH = calving to first heat; NINS = number of inseminations per series.
      C-LALPLIOICLAFSLA60PLACFHNINS
      C-LA0.080.29−0.53−0.88−0.720.900.21
      LPL0.27NC
      NC = not converged.
      0.29−0.040.540.490.19
      IOI0.39NC0.010.160.420.980.61
      CLAFS0.350.400.580.250.99−0.200.11
      LA600.060.420.390.630.390.880.17
      PLA0.110.300.250.370.17−0.31−0.11
      CFH0.040.200.200.350.040.140.40
      NINS0.250.920.460.190.290.240.23
      1 C-LA = commencement of luteal activity (all C-LA profiles); LPL = first luteal phase length (all LPL profiles); IOI = length of first inter-ovulatory interval; LA60 = luteal activity during the first 60 DIM; PLA = proportion of samples with luteal activity; CLAFS = commencement of luteal activity to first service; CFH = calving to first heat; NINS = number of inseminations per series.
      2 NC = not converged.

      DISCUSSION

      Endocrine fertility traits are suggested as traits of choice to improve fertility in dairy cows. The ability of these traits to explain the actual reproductive status of the cows with little influence of management makes them ideal for use in breeding (
      • Friggens N.C.
      • Chagunda M.G.G.
      Prediction of the reproductive status of cattle on the basis of milk progesterone measures: Model description.
      ;
      • Martin O.
      • Friggens N.C.
      • Dupont J.
      • Salvetti P.
      • Freret S.
      • Rame C.
      • Elis S.
      • Gatien J.
      • Disenhaus C.
      • Blanc F.
      Data-derived reference profiles with core presentation of progesterone, estradiol, LH, and FSH dynamics during the bovine estrous cycle.
      ). Moreover, the possibility for automated collection and analysis of milk progesterone is now resulting in availability of large quantities of data. In the present study, we evaluated fertility traits of SR and Holstein cows from 14 herds, based on their progesterone records obtained from Herd Navigator in-line milk progesterone analysis.

      Descriptive Evaluation of Endocrine and Classical Fertility Traits

      The average C-LA was 32.7 d for Holsteins and 35.2 d for SR cows (Table 2).
      • Bulman D.C.
      • Lamming G.E.
      Milk progesterone levels in relation to conception, repeat breeding and factors influencing acyclicity in dairy cows.
      reported that ovulation occurs on average 5 d before C-LA, which would imply that ovulation occurred approximately on d 28 to 30 postpartum. These results for C-LA are in agreement with previous findings for SR and Holstein cows manually sampled twice weekly until cyclicity, and then once weekly until insemination (C-LA = 34 d;
      • Petersson K.J.
      • Strandberg E.
      • Gustafsson H.
      • Berglund B.
      Environmental effects on progesterone profile measures of dairy cow fertility.
      ). However, they are smaller than values reported for Dutch Holstein-Friesian herds (C-LA = 38.1 d;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ) and for Nordic Red cows in Finland (C-LA = 39.8, 37.7, and 40.6 d for parities 1, 2 and 3, respectively;
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ). However, our C-LA estimate for Holsteins is 2.5 d longer than the overall mean reported for 4 Holstein-Friesian experimental herds in Sweden, the Netherlands, the UK, and Ireland (C-LA = 31.1;
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ). The latter value was based on manual sampling twice weekly in Sweden and the Netherlands and thrice weekly in Ireland and the UK. The average C-LA for SR in the present study was 2 to 5 d shorter than that reported for Nordic Red cows in Finland, despite the fact that this breed is closely related to SR in Sweden. Differences in data recording and editing between the studies may explain the difference. These reasons may also explain the difference in C-LA observed between Holsteins in the Dutch and Swedish dairy herds. The selection pressure put on milk yield may vary between countries, which could influence C-LA (
      • Lopez H.
      • Satter L.D.
      • Wiltbank M.C.
      Relationship between level of milk production and estrous behavior of lactating dairy cows.
      ;
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ;
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ). The filtering criteria used in the present study could also have contributed to the observed differences. For instance, when defining C-LA, we used a 5-d-earlier limit (at 20 DIM) than used in
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      .
      The intervals observed for all endocrine and classical traits (except C-LA, CFH, and atypical profiles) were longer for Holstein cows than for SR cows (Table 2). For instance, SR cows had an estimated 21-d cycle of IOI and an IOI up to 4.4 d shorter than Holsteins or Holstein-Friesians in the Netherlands (
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ) and the aggregated data set of Holstein-Friesian herds from 4 countries studied in
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      . Similarly, both breeds in the present study had shorter C-LA, LPL, and ILI than Dutch Holstein-Friesians (
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). A similar trend was observed for classical traits. However, although endocrine fertility traits were similar, classical traits of SR cows in the present study had slightly higher CFI, CLI, and FLI compared with the Nordic Red cows studied by
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      . Variations in herd management and in data editing could have contributed to these differences. The first heat was observed, on average, 11.4 d after C-LA for SR, which is 2.85 d longer than reported for Nordic Red cows in Finland (
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ). The incidence of short (29.3%), normal (40.7%), and long (30.0%) IOI determined in the present study was similar to that reported for Holstein-Friesian herds in Sweden, the Netherlands, the UK, and Ireland (
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ).
      To obtain a 1-yr calving interval, cow cyclicity should resume within 60 d postpartum (
      • Walsh S.W.
      • Williams E.J.
      • Evans A.C.O.
      A review of the causes of poor fertility in high milk producing dairy cows.
      ). This is why PLA was considered an important trait, indicating the start of luteal activity and also a cow's cyclicity after ovulation (
      • Petersson K.J.
      • Gustafsson H.
      • Strandberg E.
      • Berglund B.
      Atypical progesterone profiles and fertility in Swedish dairy cows.
      ;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). We obtained a PLA of 0.63 for the combined data set, 0.59 for SR, and 0.65 for Holstein. These values are higher than those reported for Swedish, UK, and Dutch Holstein-Friesian herds (
      • Petersson K.J.
      • Berglund B.
      • Strandberg E.
      • Gustafsson H.
      • Flint A.P.F.
      • Woolliams J.A.
      • Royal M.D.
      Genetic analysis of postpartum measures of luteal activity in dairy cows.
      ;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). This could be due to differences in the number of samples taken in a fixed period and persistence of the corpus luteum, which would inflate PLA values.
      We detected differences between the breeds and across parities for C-LA, LPL, IOI, and PLA. Holstein had 2.5 d earlier C-LA and longer LPL and IOI than SR. The highest mean intervals for C-LA, LPL, and IOI were observed in fourth parity and the lowest in first parity. These results are in agreement with the study by
      • Blavy P.
      • Derks M.
      • Martin O.
      • Höglund J.K.
      • Friggens N.C.
      Overview of progesterone profiles in dairy cows.
      , who found longer LPL for Holstein than for Danish Red and longer LPL by parity number. Cows calving in spring had 4.9-d-longer luteal activity than those calving in summer. In contrast, cows calving in summer (longer daylight than darkness) had LA60 6 percentage units higher than spring calvers. It has been found that photoperiod or temperature can affect luteal activity more strongly than diet or management in SR and Holstein cows (
      • Petersson K.J.
      • Gustafsson H.
      • Strandberg E.
      • Berglund B.
      Atypical progesterone profiles and fertility in Swedish dairy cows.
      ).
      In general, low values for C-LA, CFH, CLAFS, CInt, CFI, CLI, FLI, and NINS and high values of PLA and LA60 are indicators of better fertility. Values of LPL, ILI, and IOI that are too long indicate fertility problems, but there is an optimal lower limit of these intervals. A shorter C-LA (CFH) gives the farmer an opportunity to start to inseminate early and thus achieve a short calving interval. A short CLAFS indicates that the cow tended to show heat and was inseminated soon after first luteal activity.

      Estimates of Heritability and Genetic Correlation

      The genetic analysis produced higher heritability estimates for endocrine fertility traits (e.g., h2 = 0.26, maximum heritability value of C-LA for SR cows) than for the classical fertility traits (Table 4), indicating a comparative advantage of some endocrine fertility traits over classical traits for fertility evaluation in dairy cows. The low heritability estimates for the classical fertility traits could be attributed to a multitude of factors, including, for example, that the traits are truly under poor additive genetic control relative to nonadditive and nongenetic factors (
      • Berry D.P.
      • Bastiaansen J.W.M.
      • Veerkamp R.F.
      • Wijga S.
      • Wall E.
      • Berglund B.
      • Calus M.P.L.
      Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
      ). It is also claimed that low heritability estimates together with an unfavorable relationship between production and fertility traits complicates genetic improvement (
      • Nyman S.
      • Duchemin S.I.
      • de Koning D.J.
      • Berglund B.
      Genome-wide association study of normal and atypical progesterone profiles in Holstein-Friesian dairy cows.
      ).
      • Daetwyler H.D.
      • Villanueva B.
      • Woolliams J.A.
      Accuracy of predicting the genetic risk of disease using a genome-wide approach.
      pointed out that the level of heritability for a trait influences the success of genome-wide association studies. Thus, the use of more heritable hormonal profiles related to fertility performance should help increase genetic improvement (
      • Berry D.P.
      • Bastiaansen J.W.M.
      • Veerkamp R.F.
      • Wijga S.
      • Wall E.
      • Berglund B.
      • Calus M.P.L.
      Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
      ). The heritability estimates obtained for endocrine traits of SR and Holstein cows in the current study were higher than those reported for Holstein-Friesian herds in Ireland, the UK, the Netherlands, and Sweden (e.g., C-LA: h2 = 0.13;
      • Berry D.P.
      • Bastiaansen J.W.M.
      • Veerkamp R.F.
      • Wijga S.
      • Wall E.
      • Berglund B.
      • Calus M.P.L.
      Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
      ) and for Holstein-Friesian herds in the Netherlands based on in-line milk P4 information using Herd Navigator (h2 = 0.14;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). However, the values were similar to that reported for Nordic Red cows in Finland (h2 = 0.24, using univariate model;
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      ).
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      reported heritability estimates of 0.18, 0.08, 0.08, and 0.03 for C-LA, ILI, LPL, and IOI, respectively, for Swedish, Dutch, UK, and Irish Holstein-Friesian herds using manually collected progesterone information.
      Our heritability estimates for PLA (0.12 for Holstein and 0.14 for SR) are in agreement with that (0.12) reported for Dutch Holstein-Friesian herds (
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ).
      • Petersson K.J.
      • Berglund B.
      • Strandberg E.
      • Gustafsson H.
      • Flint A.P.F.
      • Woolliams J.A.
      • Royal M.D.
      Genetic analysis of postpartum measures of luteal activity in dairy cows.
      found lower estimates of heritability for PLA with more infrequent sampling (0.25, 0.20, and 0.14 for weekly, twice-monthly, and monthly sampling, respectively). Moderate heritability (0.16) for CInt was obtained for Holstein cows in the present study, but zero heritability for SR cows. The former was higher than the heritability of other classical traits in our study and compared with previous findings for Dutch herds (0.05;
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      ). However, it is important to note that we obtained zero heritability for several endocrine and classical fertility traits in SR cows. Similar findings were made in the study in Finland on Nordic Red cows (except for C-LA and CFH); therefore, these variables were excluded from their paper (Jarmo Juga, University of Helsinki, Helsinki, Finland, personal communication, 2018). This suggests that further studies are necessary on larger groups of Red cows in the Nordic countries to obtain reliable genetic parameter estimates for this breed. Length of first ILI consistently showed zero heritability estimates for SR, Holstein, and both breeds combined, suggesting that this may not be a useful endocrine trait to include in the breeding goal. Most of the correlation estimates had high standard errors, because of the small sample size (
      • Berry D.P.
      • Bastiaansen J.W.M.
      • Veerkamp R.F.
      • Wijga S.
      • Wall E.
      • Berglund B.
      • Calus M.P.L.
      Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
      ). The strong correlations of C-LA with LA60 (−0.88), PLA (−0.72), and CFH (0.90) suggest that these traits reflect almost the same parameter. Similar strong correlations for C-LA with LA60, PLA, and CFH have been reported by
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      ,
      • Tenghe A.M.
      • Bouwman A.C.
      • Berglund B.
      • Strandberg E.
      • Blom J.Y.
      • Veerkamp R.F.
      Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
      , and
      • Häggman J.
      • Christensen J.M.
      • Mäntysaari E.A.
      • Juga J.
      Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
      . In contrast to our findings,
      • Berry D.P.
      • Bastiaansen J.W.M.
      • Veerkamp R.F.
      • Wijga S.
      • Wall E.
      • Berglund B.
      • Calus M.P.L.
      Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
      obtained a very weak correlation (0.29) between C-LA and CFH. Consistent with
      • Nyman S.
      • Johansson K.
      • de Koning D.J.
      • Berry D.P.
      • Veerkamp R.F.
      • Wall E.
      • Berglund B.
      Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
      , IOI had a strong correlation with CFH (0.88) in the present study.
      • Petersson K.J.
      • Berglund B.
      • Strandberg E.
      • Gustafsson H.
      • Flint A.P.F.
      • Woolliams J.A.
      • Royal M.D.
      Genetic analysis of postpartum measures of luteal activity in dairy cows.
      suggest that PLA has an intermediate optimum value: too high a value for PLA is unfavorable, because it is associated with persistence of corpus luteum type 1, but too low a value could also be unfavorable, because it is associated with delayed ovulation type 1.
      Overall, we believe that endocrine traits have the potential to better reflect the biological fertility status of the dairy cows than classical traits do. Therefore, they should be considered as part of the selection criteria and the breeding goal. A prerequisite is that automated progesterone scoring or recording resources are made available for routine recording as part of national databases.

      CONCLUSIONS

      Most of the endocrine traits studied, including C-LA, LA60, PLA, and CFH, were found to be moderately heritable and strongly correlated with each other. Classical traits had lower heritability estimates and weak genetic correlations with endocrine traits. Our results confirm previous findings that endocrine traits might have the potential to better describe the actual reproductive status of dairy cows than classical traits, and suggests that some endocrine fertility traits, especially those related to commencement of luteal activity, may increase the rate of genetic improvement if used as additional selection criteria and in the breeding goal. Swedish Red cows appeared superior to Holstein in some fertility traits, but showed zero heritability for LPL, IOI, CLAFS, and CFH, which warrants further studies on a larger sample of Red cattle. It is important that endocrine progesterone data from automatic devices be made available to recording and breeding organizations in the future and included in a central database.

      ACKNOWLEDGMENTS

      We acknowledge the herd owners in Sweden who contributed data on their herds. This study was performed within the Swedish part of the Nordic research project NorFert (Improving Nordic dairy cow fertility through genetics). The main participant is the Swedish University of Agricultural Sciences (SLU). Industry partners are Växa Sverige (Stockholm, Sweden), VikingGenetics (Randers, Denmark), Nordic Cattle Genetic Evaluation Ltd. (NAV, Aarhus, Denmark), and Lattec I/S (Hillerod, Denmark). NorFert project partners are the Natural Resource Institute of Finland (Luke, Jokioinen, Finland) and Helsinki University, Helsinki, Finland, Aarhus University, Aarhus, Denmark, and NMBU, Aas, Norway. We are grateful for the funding for the project provided by the Swedish Farmers' Foundation for Agricultural Research (Project ID: O-15-20-587). The authors declare no conflicts of interest.

      REFERENCES

        • Adriaens I.
        • Huybrechts T.
        • Geerinckx K.
        • Daems D.
        • Lammertyn J.
        • De Ketelaere B.
        • Saeys W.
        • Aernouts B.
        Mathematical characterization of the milk progesterone profile as a leg up to individualised monitoring of reproduction status in dairy cows.
        Theriogenology. 2017; 103 (28779608): 44-51
        • Adriaens I.
        • Saeys W.
        • Huybrechts T.
        • Lamberigts C.
        • François L.
        • Geerinckx K.
        • Leroy J.
        • Ketelaere B.D.
        • Aernouts B.
        A novel system for on-farm fertility monitoring based on milk progesterone.
        J. Dairy Sci. 2018; 101 (29935821): 8369-8382
        • Berry D.P.
        • Bastiaansen J.W.M.
        • Veerkamp R.F.
        • Wijga S.
        • Wall E.
        • Berglund B.
        • Calus M.P.L.
        Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries.
        Animal. 2012; 6 (23217223): 1206-1215
        • Blavy P.
        • Derks M.
        • Martin O.
        • Höglund J.K.
        • Friggens N.C.
        Overview of progesterone profiles in dairy cows.
        Theriogenology. 2016; 86 (27177962): 1061-1071
        • Bruinjé T.C.
        • Colazo M.G.
        • Ribeiro E.S.
        • Gobikrushanth M.D.
        • Ambrose J.
        Using in-line milk progesterone data to characterize parameters of luteal activity and their association with fertility in Holstein cows.
        J. Dairy Sci. 2019; 102 (30343919): 780-798
        • Bulman D.C.
        • Lamming G.E.
        Milk progesterone levels in relation to conception, repeat breeding and factors influencing acyclicity in dairy cows.
        J. Reprod. Fertil. 1978; 54 (722692): 447-458
        • Daetwyler H.D.
        • Villanueva B.
        • Woolliams J.A.
        Accuracy of predicting the genetic risk of disease using a genome-wide approach.
        PLoS One. 2008; 3 (18852893)e3395
        • DeLaval International
        Herd Navigator and reproduction management.
        • Edwards S.M.
        Standard errors of heritability estimates.
        http://www.iysik.com/dmu/heritabilities
        Date: 2017
        Date accessed: February 11, 2019
        • Friggens N.C.
        • Chagunda M.G.G.
        Prediction of the reproductive status of cattle on the basis of milk progesterone measures: Model description.
        Theriogenology. 2005; 64 (15935851): 155-190
        • Häggman J.
        • Christensen J.M.
        • Mäntysaari E.A.
        • Juga J.
        Genetic parameters for endocrine and traditional fertility traits, hyperketonemia and milk yield in dairy cattle.
        Animal. 2019; 13 (29954471): 248-255
        • Leroy J.L.M.R.
        • Opsomer G.
        • Van Soom A.
        • Goovaerts I.G.F.
        • Bols P.E.J.
        Reduced fertility in high-yielding dairy cows: are the oocyte and embryo in danger? Part I. The importance of negative energy balance and altered corpus luteum function to the reduction of oocyte and embryo quality in high-yielding dairy cows.
        Reprod. Domest. Anim. 2008; 43 (18384499): 612-622
        • Lopez H.
        • Satter L.D.
        • Wiltbank M.C.
        Relationship between level of milk production and estrous behavior of lactating dairy cows.
        Anim. Reprod. Sci. 2004; 81 (14998648): 209-223
        • Madsen P.
        • Jensen J.
        A user's guide to DMU. A package for analyzing multivariate mixed models. Version 6, release 5.2.
        Department of Genetics and Biotechnology, Faculty of Life Sciences, University of Aarhus, Research Centre Foulum, Tjele, Denmark2013
        • Martin O.
        • Friggens N.C.
        • Dupont J.
        • Salvetti P.
        • Freret S.
        • Rame C.
        • Elis S.
        • Gatien J.
        • Disenhaus C.
        • Blanc F.
        Data-derived reference profiles with core presentation of progesterone, estradiol, LH, and FSH dynamics during the bovine estrous cycle.
        Theriogenology. 2013; 79 (23149413): 331-343
        • Nyman S.
        • Duchemin S.I.
        • de Koning D.J.
        • Berglund B.
        Genome-wide association study of normal and atypical progesterone profiles in Holstein-Friesian dairy cows.
        J. Dairy Sci. 2019; 102 (30799107): 3204-3215
        • Nyman S.
        • Johansson K.
        • de Koning D.J.
        • Berry D.P.
        • Veerkamp R.F.
        • Wall E.
        • Berglund B.
        Genetic analysis of atypical progesterone profiles in Holstein-Friesian cows from experimental research herds.
        J. Dairy Sci. 2014; 97 (25151875): 7230-7239
        • Petersson K.J.
        • Berglund B.
        • Strandberg E.
        • Gustafsson H.
        • Flint A.P.F.
        • Woolliams J.A.
        • Royal M.D.
        Genetic analysis of postpartum measures of luteal activity in dairy cows.
        J. Dairy Sci. 2007; 90 (17183111): 427-434
        • Petersson K.J.
        • Gustafsson H.
        • Strandberg E.
        • Berglund B.
        Atypical progesterone profiles and fertility in Swedish dairy cows.
        J. Dairy Sci. 2006; 89 (16772571): 2529-2538
        • Petersson K.J.
        • Strandberg E.
        • Gustafsson H.
        • Berglund B.
        Environmental effects on progesterone profile measures of dairy cow fertility.
        Anim. Reprod. Sci. 2006; 91 (15950407): 201-214
        • Royal M.D.
        • Darwash A.O.
        • Flint A.P.F.
        • Webb R.
        • Woolliams J.A.
        • Lamming G.E.
        Declining fertility in dairy cattle: Changes in traditional and endocrine parameters of fertility.
        Anim. Sci. 2000; 70: 487-501
        • Royal M.
        • Mann G.E.
        • Flint A.P.F.
        Strategies for reversing the trend towards subfertility in dairy cattle.
        Vet. J. 2000; 160 (10950135): 53-60
        • Saint-Dizier M.
        • Chastant-Maillard S.
        Towards an automated detection of oestrus in dairy cattle.
        Reprod. Domest. Anim. 2012; 47 (22214367): 1056-1061
        • SAS Institute
        SAS/STAT® 12.1 User's Guide.
        SAS Institute Inc., Cary, NC2012
        • Tenghe A.M.
        • Bouwman A.C.
        • Berglund B.
        • Strandberg E.
        • Blom J.Y.
        • Veerkamp R.F.
        Estimating genetic parameters for fertility in dairy cows from in-line milk progesterone profiles.
        J. Dairy Sci. 2015; 98 (26004838): 5763-5773
        • Walsh S.W.
        • Williams E.J.
        • Evans A.C.O.
        A review of the causes of poor fertility in high milk producing dairy cows.
        Anim. Reprod. Sci. 2011; 123 (21255947): 127-138