Authors: Anna Wirth, Jürgen Duda, Ottmar Distl
Categories: Article, German Brown cows, ancestral inbreeding, inbreeding depression, longevity
Source: Animals : an Open Access Journal from MDPI
Doi: 10.3390/ani13172765
Observed increasing levels of inbreeding in German brown cattle imply the need to evaluate the resulting effects. If increasing levels of inbreeding come along with a reduction in the phenotypic performance, this is referred to as inbreeding depression. However, not all inbreeding is assumed to be unfavorable. In the context of inbreeding in combination with selection, it is presumed that negative variants can be eliminated from the population over time, implying ancestral inbreeding to be less harmful than new inbreeding. This is referred to as purging. To evaluate the effect of purging, ancestral inbreeding coefficients have been developed to allow to distinguish between new and ancestral inbreeding. The aim of this study was to evaluate inbreeding depression and purging for lifetime, lifetime performance, survival, and culling rates due to different reasons in German brown cows, by calculating the effects of classical and ancestral inbreeding coefficients. All longevity traits under study were affected by inbreeding depression. As the effects of ancestral inbreeding coefficients were significantly negative for lifetime and lifetime performance, while new inbreeding had no significant effect, there was no evidence of purging in the population under study. Thus, considering inbreeding levels in future mating plans helps to avoid a further decline in longevity due to inbreeding.
A recent study on the population structure of the German Brown population found increasing levels of classical and ancestral inbreeding coefficients. Thus, the aim of this study was to determine the effects of inbreeding depression and purging on longevity traits using classical and ancestral inbreeding coefficients according to Kalinowski (2002) (Fa_Kal, FNew), Ballou (1997) (Fa_Bal), and Baumung (2015) (Ahc). For this purpose, uncensored data of 480,440 cows born between 1990 and 2001 were available. We analyzed 17 longevity traits, including herd life, length of productive life, number of calvings, lifetime and effective lifetime production for milk, fat, and protein yield, the survival to the 2nd, 4th, 6th, 8th, and 10th lactation number, and the culling frequencies due to infertility, or udder and foot and leg problems. Inbreeding depression was significant and negative for all traits but for culling due to udder and to foot and leg problems. When expressed in percentages of genetic standard deviations, inbreeding depression per 1% increase in inbreeding was −3.61 to −10.98%, −2.42 to −2.99%, −2.21 to −4.58%, and 5.13% for lifetime production traits, lifetime traits, survival rates, and culling due to infertility, respectively. Heterosis and recombination effects due to US Brown Swiss genes were positive and counteracted inbreeding depression. The effects of FNew were not significantly different from zero, while Fa_Kal had negative effects on lifetime and lifetime production traits. Similarly, the interaction of F with Fa_Bal was significantly negative. Thus, purging effects could not be shown for longevity traits in German Brown. A possible explanation may be seen in the breed history of the German Brown, that through the introgression of US Brown Swiss bulls ancestral inbreeding increased and longevity decreased. Our results show, that reducing a further increase in inbreeding in mating plans is advisable to prevent a further decline in longevity due to inbreeding depression, as purging effects were very unlikely in this population.
Keywords: inbreeding depression, ancestral inbreeding, longevity, German Brown cows
Targeted selection programs in dairy cattle populations led to an enormous genetic progress, especially in milk production traits, over recent decades, at the cost of increasing levels of inbreeding as a result of focusing on a low number of highly selected elite bulls. German Brown is characterized by high milk production with high quality and superiority in length of productive life compared to Holsteins and German Fleckvieh. Introgression of US Brown Swiss bulls since 1966 contributed to the improvement in milk production traits in German Brown and, on the other hand, has led to increasing levels of individual inbreeding in animals with higher proportions of US Brown Swiss genes [1]. A previous study on milk performance for the first three lactation numbers of German Brown from the birth years between 1980 and 1992 indicated an inbreeding depression by 10 to 13 kg milk yield, 0.36 to 0.51 kg fat yield, and 0.35 to 0.47 kg protein yield per 1% increase in the inbreeding coefficient [2]. A slight tendency to lower inbreeding coefficients was seen in cows surviving the first three lactations compared to cows that left the herd before completing three lactations [2]. In Swiss Brown cows, first lactation and milk and fat yield with first calving between 1971 and 1986 were reduced by 26 and 0.08 kg per 1% of the inbreeding coefficient [3]. In Austrian Brown cows with first parturition between 1979 and 1991, estimates for inbreeding depression by 1% increase in the inbreeding coefficient were 6.3 to 10.4 kg milk yield, but close to zero for fat and protein percentages [4]. Length of productive life and lifetime energy-corrected milk yield decreased by 4.3 days and 165 kg per 1% inbreeding, respectively [4]. In addition, juvenile mortality of calves and heifers of Austrian Brown, born in the years 2001 to 2007, have been shown to be negatively affected by 0.49% per 1% increase in inbreeding [5]. Next to the effect of inbreeding depression, there is an ongoing debate about the existence of positive effects of inbreeding caused by inbred ancestors and referred to as purging [6,7]. Purging assumes that intensive selection may lead to an elimination of deleterious alleles from a population and therefore counteracting inbreeding depression [8]. This means that after a few generations of inbreeding, the best-performing individuals survive and reach the performance level of the non-inbred or less inbred individuals, or even achieve higher performance, while the poorer-performing individuals carrying deleterious alleles die or do not reproduce, so that these alleles are eliminated from the population [9,10]. According to Dickerson [11], the response to selection is stronger in mating systems with complete sibling mating than in systems with random mating. To date, the presence of purging has been studied in various experiments, including self-fertilization in plants [12], sibling mating in animals [13,14,15,16], by mainly examining traits related to life history [17], or in small captive populations [18,19]. For example, in a previous breeding experiment with quails, it was shown that populations with inbreeding in previous generations were less prone to inbreeding depression than populations without inbreeding history. Furthermore, the reproductive performance of populations with intensive inbreeding was better after a few generations than that of populations with random mating [13,20]. However, the intensity of purging depends on the extent to which a trait is affected by deleterious alleles and how these alleles affect fertility and survival. Lethal or semi-lethal alleles are more easily eliminated than less deleterious ones [8,9,17]. To analyze the effect of purging based on pedigree data, ancestral inbreeding coefficients have to be calculated [18,21,22]. Only few studies reported ancestral coefficients for different Holstein dairy cattle populations. The effects of the ancestral inbreeding coefficients on production, fertility, birth traits, and survival traits were not consistent [6,7,23,24,25]. For example, in Irish and Dutch Holstein evidence of the presence of purging effects for production traits could be shown [6,23], whereas results for Canadian and Iranian Holstein were not significant or even contradictory [24,25]. As the extent of inbreeding depression as well as purging may depend on different factors such as the population structure, the selection history, and the considered trait [8,10], it is crucial to analyze and monitor these effects in every population to adopt appropriate measures. For the German Brown population, classical and ancestral inbreeding coefficients have recently been reported and increasing trends with birth years were shown for classical and ancestral inbreeding coefficients [26]. The breed proportion of US Brown Swiss is increasing, but length of life and productive life were not positively correlated with increasing breed proportions of US Brown Swiss; therefore, this follow-up study in German Brown has been initiated with the objective to analyze inbreeding depression and purging effects on longevity traits including herd life (HL), length of productive life (LPL), lifetime milk production, and survival and culling incidences. We used an uncensored data set of Brown German cows born in the years 1990 to 2001. In order to account for the US Brown Swiss breed percentage, heterosis and recombination effects were considered simultaneously in animal models.
The data were obtained from the official milk recording organization of Bavaria (Landeskuratorium der Erzeugerringe für tierische Veredelung in Bayern e.V., LKV). This data set comprised all German Brown with their available pedigrees and lifetime milk performance since 1990 [26]. For the analysis of longevity traits, all first-calving heifers born between 1990 and 2001 were regarded. All these animals had known parents as well as registered dates and reasons why they left the herd and, therefore, they had complete lifetime records. The final data set contained 480,440 cows with known parents and a corresponding pedigree file of 820,558 individuals. The data set was used to analyze measures of longevity, defined as functional, as they were corrected for the relative fat and protein yield of the single cow within their herd [27,28]. As suggested by Schuster et al. (2020) [29], herd life (HL) was defined as the time period from birth until culling. Length of productive life (LPL) included the time period from first calving to culling. Number of lactations (NC) was the sum of lactation numbers per lifetime which the respective cow reached. Lifetime milk yield (LMY), fat yield (LFY), and protein yield (LPY) were defined as total yield across all lactations. Lifetime efficiency for milk yield, (EffLMY), fat yield (EffLFY), and protein yield (EffLPY) were calculated as the ratio of lifetime production to herd life in days. Survival to lactation two (Surv1), four (Surv3), six (Surv5), eight (Surv7), and ten (Surv9) was calculated as the proportion of cows that reached the respective lactation number, e.g., cows that survived the first lactation and had at least a second calving date were encoded with Surv1 = 1 and all other cows with 0. Results are displayed as relative frequencies. The frequencies of the three most common culling reasons, namely culling due to foot and leg problems (CulCL), infertility (CulINF), and udder problems (CulUD), were analyzed as binary traits, indicating that the cow did leave the herd due to this respective specific reason. Classical and ancestral inbreeding coefficients were available from our previous study and merged with the present data set including 448, 440 cows [26].
To account for the effect of crossbreeding with US Brown Swiss bulls starting in 1966, the coefficients of heterosis (HET = PS (1 − PD) + PD (1 − PS)) and recombination loss (REC = PD (1 − PD) + PS (1 − PS)) were defined based on the maternal (PD) and paternal (PS) breed proportion of Brown Swiss genes [30].
The effect of the different inbreeding coefficients was calculated by regressing phenotypic values on the inbreeding coefficient and accounting for animal and environmental effects. The following linear univariate animal model, parameterized for non-genetic effects according to Punsmann, Duda and Distl [28], was employed
where Yijklmnopqr = dependent variate with HL, LPL, NC, LMY, LFY, LPY, EffLMY, EffLFY, EffLPY, Surv1, Surv3, Surv5, Surv7, Surv9, CulCL, CulINF, or CulUD. BYi is the ith class of birth year for i = 1 (1990–1993), 2 (1994–1997), and 3 (1998–2001), BMj is the jth class of birth month for j = 1–12, AFCk = age at first calving in months and b1 the corresponding linear regression coefficient on the age at first calving, WHPLl = within herd relative production level calculated as deviation of the sum of fat and protein production of each cow from the herd mean for fat and protein yield and expressed in standard deviation classes for l = 1–12, HPLm = herd production level of fat and protein yield expressed in standard deviations for m-classes from 1 to 8, HETn = coefficient of heterosis effects, RECo = coefficient of recombination effects, INBp = inbreeding coefficient F, Fa_Bal, or Ahc, b2, b3, and b4 are the linear regression coefficients of heterosis, recombination, and inbreeding effects, respectively, animal = random animal effect, herd = random herd effect, e = random error term.
The expected values for F, HET, and REC in the current population were calculated based on the observed estimates for the regression coefficients obtained from model 1 and the average coefficients of F, HET, and REC that were 0.018, 0.441, and 0.369, respectively. As Fa_Bal does not consider the probability that an individual is identical by descent (IBD) itself, we included the interaction of Fa_Bal with F as suggested by Ballou (1997) [18] in model 2
where fixed and random effects were defined as above and F is the coefficient of inbreeding, F×Fa_Bal = interaction of the coefficient of inbreeding with the ancestral inbreeding coefficient according to Ballou [18], and b4 and b5 are the corresponding linear regression coefficients.
Model 3 included the effects of either Ahc (ancestral history coefficient according to Baumung, Farkas, Boichard, Mészáros, Sölkner and Curik [22]) or Fa_Bal (ancestral inbreeding coefficient according to Ballou [18]) as FANC together with F
In model 4 the effect of the inbreeding coefficients according to Kalinowski, Hedrick and Miller [21] were analyzed by considering Fa_Kal and FNew simultaneously in the
Variance components were estimated using VCE 6.0.2 [31], and were used for estimation of random and fixed effects as well as linear regressions using PEST, version 4.2.6. Further statistical analyses were performed in SAS, version 9.4 (Statistical Analysis System, Cary, NC, USA, 2023). The model employed for the estimation of heritabilities included all fixed and random effects models as the above-mentioned models, and the proportion of US Brown Swiss in classes of 10% (<31%, 31–40%, 41–50%, 51–60%, 61–70%, 71–80%, >80%) [21]. Animal models with a binomial distribution function for 0/1-traits yielded very similar estimates for heritabilities, and thus, we employed linear animal models.
Phenotypic means of the longevity measures of German Brown cows born between 1990 and 2001 are presented in Table 1. German Brown cows reached an average age of 6.15 years and were productive for 56.2% (3.46 years) of their lifetime. The 95% quantiles for HL and LPL were 11.16–20.98 and 8.48–18.09 years with 8–17 calvings, respectively. Quartiles for LMY were from 7537 to 32,736 kg milk and the top 5% cows had an LMY of 58,065–191,049 kg milk. For LFY and LPY, the top 5% cows produced 2418–8045 kg fat and 2075–6787 kg protein, respectively.
Heritability estimates for HL, LPL, and NC were around 0.16 (Table 2). Effective lifetime production showed higher heritabilities (0.128–0.154) compared to lifetime production (0.126–0.135). For survival traits, heritability was highest for Surv5 (0.088) and lowest for Surv1 (0.025), whereas heritability for the three most common culling reasons ranged from 0.034 to 0.055. Heritability estimates with their corresponding standard errors for Surv1, Surv3, Surv5, Surv7, and Surv9 were 0.048 (0.005), 0.111 (0.006), 0.182 (0.007), 0.282 (0.011), and 0.319 (0.024), and for CulCL, CulINF, and CulUD 0.147 (0.008), 0.063 (0.005), and 0.099 (0.007), respectively, when transformed to the underlying scale with the formula of Dempster and Lerner [32].
Pedigree completeness was measured through the number of complete equivalent generations (GE). The GE was on average 5.78 with an increasing trend from 4.54 in 1990 to 6.13 in 2001. The average individual inbreeding coefficient for the reference population was 0.018 and increased over time from 0.013 in 1990 to 0.024 in 2001. Of the 480,440 cows, 88.5% were inbred, with the highest coefficient of individual inbreeding being 0.3175. The mean ancestral coefficients of inbreeding ranged from 0.001 to 0.012 and showed an increasing trend over time, with the highest increase observed for Fa_Bal and Ahc, which increased from about 0.013 in 1990 to 0.019 in 2001 (Supplementary Table S1). Correlation between Fa_Bal and Ahc was 0.99. The lowest correlations (0.07) were found between FNew and Fa_Bal as well as Ahc. Correlations between coefficients of heterosis and coefficients of inbreeding (F), FNew, Fa_Kal, Fa_Bal, and Ahc were −0.33, −0.07, −0.06, −0.11 and −0.11, respectively (Supplementary Table S2).
The results from model 1 revealed a significant unfavorable effect of F on all longevity traits studied, except for CulCL and CulUD (Table 3). All estimates from linear regression coefficients refer to 100% inbreeding and consider the difference between a non-inbred and a fully inbred animal. A 1% increase in inbreeding reduces HL and LPL by 7.3 and 7.7 days, respectively. Regarding survival, inbreeding depression was higher for cows with less lactation numbers than for cows with more lactation numbers and was highest for Surv3 (Table 3). When expressed in percentages of genetic standard deviations, lifetime production traits were more influenced by inbreeding than HL, LPL, NC, and survival to the following lactation number (Surv1 to Surv9), which had decreasing inbreeding depression when cows reached higher lactation numbers (Supplementary Table S3). Inbreeding significantly increased the frequency of CulINF, with an increase by 0.399 per 100% inbreeding. Significant positive heterosis effects were observed for most traits with 100% heterosis increasing HL by 161 days and LMY by 1700 kg, and increasing Surv3 and Surv5 by 6.8% and 5.9%, respectively. Recombination effects were also significantly positive for HL, LPL, NC, and Surv1 to Surv7. No significant heterosis and recombination effects were found for effective lifetime production traits (Table 3). Regarding the expected cumulated effects of inbreeding, heterosis, and recombination for the current population, the combined effects were favorable for all traits apart from effective lifetime production and culling due to infertility. Projecting these effects to the birth year 2014, based on coefficients from 2014, heterosis and recombination effects decrease, while inbreeding effects increase, leading to expected decreasing combined effects and even to more unfavorable negative effects for LPY and effective lifetime production (Supplementary Table S4).
Results for the effects of ancestral inbreeding coefficients are presented in Table 4 and Table 5. Using model 2, the regression coefficients of the interaction of F×Fa_Bal were significantly negative for HL, LPL, NC, lifetime, and effective lifetime production traits and the corresponding effects of F were slightly reduced compared with the results of model 1.
Significant negative effects of Fa_Kal were found for HL, LPL, NC, Surv3, and lifetime production traits using model 4, whereas the effect of FNew was close to zero and not significantly different from zero for these traits (Table 4).
Considering Ahc, negative regression coefficients were observed for HL, LPL, NC, lifetime production, and Surv3 and 5. Estimates for regression coefficients from model 4 were similar to the estimates obtained from model 1 (Table 5). Moreover, the regression coefficients for Ahc equaled those of Fa_Bal in the corresponding models and, therefore, only the results for Ahc are presented.
The expected phenotypes of highly (95% percentile) and lowly (5% percentile) inbred cows were calculated to further illustrate the effect of changing the inbreeding coefficients. For classical inbreeding, the differences corresponded to 33 days for HL and 873 kg for LMY. Differences in ancestral inbreeding levels were lower with 7 and 11 days of HL and 125 and 174 kg of LMY for Fa_Kal and Ahc, respectively (Supplementary Table S5).
Inbreeding depression was obvious for all traits analyzed but for CulCL and CulUD. On the other hand, all longevity traits and lifetime performance traits but lifetime efficiency showed significantly positive heterosis effects through crossbreeding with US Brown Swiss bulls. Recombination effects were also positive but smaller than heterosis effects and had importance for longevity traits only. Models regarding ancestral inbreeding coefficients revealed significant negative effects, and thus we were not able to demonstrate purging effects. Under purging, ancestral inbreeding would have been expected to exert positive effects on the traits.
The negative effect of inbreeding on LPL of 7.7 (7.70 = 2.11 years × 365 days/100) days per 1% inbreeding was slightly higher than estimates reported for Austrian Brown, Simmentals including Simmental x Red Friesian crosses, and pure Simmentals with −4.3, −5.7, and −4.9 days per 1% inbreeding, respectively [4]. In a US American study for 2.6 million Holsteins from 1983 to 1993, LPL was reduced by 5.96 days per 1% inbreeding [33].
Lifetime production traits were less negatively affected by inbreeding in Austrian Simmental, Simmental x Red Friesian, and Brown Swiss compared to this study. In the latter breeds, fat-corrected milk was reduced by −136.7, −109.9, and −165.0 kg per 1% inbreeding, respectively [4]. Similarly, US Holsteins showed a decrease due to 1% inbreeding of −177.17 kg, −6.01 kg, and −5.45 kg lifetime milk, and fat and protein yield, respectively [33]. Nevertheless, lifetime was limited in the Austrian study to 10 lactations [4] and in the US study to 84 months after first calving [33], whereas our study included all records until the cows left the herd.
The higher impact of inbreeding on lifetime performance compared to longevity traits, HL, LPL, and NC, becomes obvious when expressed in percentages of their phenotypic and genetic standard deviations. This may indicate that the reduced lifetime production in German Brown cattle in this study is not only due to the shortened LPL, but the production traits themselves may be negatively influenced by inbreeding. This is in agreement with previous studies on different dairy cattle populations that reported inbreeding depression for milk performance in the first three lactations [6,23,24,25]. The highest effects due to inbreeding depression, expressed in percentages of the genetic standard deviations, were demonstrated for LMY, Surv1, and EffLFY. Therefore, inbreeding depression may have the largest negative effects on the expected breeding progress of these traits.
In order to evaluate contributions from US Brown Swiss genes, we also employed a model with F and its interaction with classes of US Brown Swiss genes. These analyses showed, for all traits under analysis, that cows with less than 50% US Brown Swiss genes had less inbreeding depression and inbreeding depression increased up to 70–80% US Brown Swiss genes and then decreased with even higher proportions of US Brown Swiss genes. This may indicate that negative effects of inbreeding seem related with proportions of US Brown Swiss genes and the strategy of using bulls with different proportions of US Brown Swiss genes. Introducing US Brown Swiss bulls may be associated with less inbreeding depression, because these US Brown Swiss bulls may not be so closely related with the German Brown population in comparison to German Brown bulls with less than 100% US Brown Swiss genes.
Inbreeding depression may have caused a higher risk for younger cows to be culled. Cows which survived more lactations may be assumed to be prone to a lower risk to be culled due to the effects on inbreeding depression. Thus, cows surviving to higher lactation numbers are expected to show decreasing inbreeding coefficients (Supplementary Tables S5c and S6). The decreasing inbreeding coefficients in cows reaching a higher lactation number reduce the size of inbreeding depression in these cows. Selection of cows for following lactation numbers may also lower the increase in inbreeding and the amount of inbreeding depression in the next generation. Survivors to higher lactations may have more progeny than cows leaving the herd early and therefore, increase in inbreeding and inbreeding depression may be smaller compared to a situation when all cows have the same chance to produce progeny. The regression coefficients were highest for Surv3, implying that an inbred cow is at highest risk to leave the herd before the fourth lactation. As this is also the period of time when cows leave the herd on average, i.e., with 3.51 calvings and an LPL of 3.46 years in the studied population, this further supports the lifetime-limiting effect of inbreeding. Nevertheless, the effect of inbreeding on survival was small, as it would explain a reduction of 0.8% (−0.8% = (0.018) × (−0.457) ×100) for Surv3, in relation to the average degree of inbreeding in the population studied. With respect to survival to the second lactation, a comparable regression coefficient of −0.3 was reported for Irish Holstein based on 42,723 survival records [23]. The mean level of inbreeding was higher in Irish Holstein with F = 0.0268 than in German Brown with 0.018 though. A decrease in survival from the second through fifth and sixth lactation in comparison to lactation one was observed in American Holstein and Jersey, respectively, with increasing inbreeding [34,35]. When considering the stayability to 48 months of age, the effects of inbreeding in US Ayrshire, Guernsey, Holstein, Jersey, and Brown Swiss ranged from −0.011 to −0.002 in a linear sire model and were closer to zero than the estimates for Surv1 and Surv3 in this study [36]. Because of the high standard errors of 0.0017, the authors concluded that inbreeding had no negative effect on stayability to 48 months [36]. However, it has to be considered that the degree of inbreeding in this study was as low as 0.009. When longevity was expressed as relative culling risk, small but significant effects were found in Canadian Holstein, Jersey, and Ayrshire compared with non-inbred cows. The effect was greater for cows in inbreeding classes beyond 12.5% [37]. In US Jersey cows that calved for the first time between 1981 and 2000, a slightly higher culling risk was observed in animals with an inbreeding coefficient greater than 10% than in animals with an inbreeding coefficient less than 5% [38]. Regarding the reasons why cows left the herd, the increase in the frequency of culling due to infertility was related with inbreeding depression. Infertility is the most common reason in German Brown cows for leaving the herd and is one of the most common reasons in dairy cattle populations [39,40,41]. The proportion of cows leaving the herd due to infertility has already been shown to be highest in lower lactation numbers in German Brown [42] and Australian dairy cattle [41]. Cows surviving fewer lactations were more inbred in this study than cows surviving higher lactations; thus, it can be concluded that German Brown cows with higher levels of inbreeding are more likely to be affected by infertility, which is in agreement with the results on Spanish Holsteins, where cows with high inbreeding coefficients (6.25–12.5%) have a 1.68% lower pregnancy rate than non-inbred cows [43]. In general, negative effects of inbreeding on fertility have already been reported for different dairy cattle populations [6,23,43,44,45]. Previous studies reported positive correlations between fertility and longevity traits, further explaining the link between inbreeding, longevity, and fertility [46,47]. Thus, preventing inbreeding may improve fertility and reduce culling due to infertility, which would lead to an increasing HL and LPL [42].
We concurrently considered heterosis and recombination effects because the German Brown is a crossbred population of Original Brown and US Brown Swiss with an increasing trend of breed proportions of US Brown Swiss [26]. Considering the effects in relation to the mean heterosis and inbreeding coefficients of the population under study, heterosis effects are positive and are counteracting the negative impact of inbreeding. Using population averages, heterosis is expected to increase HL by 70 days, while inbreeding decreases HL by 13 days. Similar expectations were found for lifetime production traits. Heterosis estimates for Surv1-5 in this study range from 3.4% to 6.8% and are in the range of the results from studies of crossbred populations involving various dairy cattle breeds from New Zealand [48], Sweden [49], and Denmark [50]. The introduction of US Brown Swiss genes was performed with only a small number of sires that were mated frequently within the population. Thus, the high breed proportion of the today’s German Brown population stems from a small gene pool of US Brown Swiss genetics [26]. This means that crossing with US Brown Swiss sires led to increasing levels of inbreeding, as co-ancestries between cows rapidly increased in the first crossbred generations [26]. This may furthermore explain the positive and significant recombination effects for HL, LPL, NC, and survival. When summarizing the effects of heterosis, recombination, and inbreeding, the positive combined effects suggest a balanced contribution of inbreeding and the introduction of new genes of US Brown Swiss. When projecting these results to the 2014 birth cohort, based on the data of our previous study [26], the combined effects for all traits decrease, mainly due to the increasing trend in inbreeding, and the increasing breed proportion of US Brown Swiss is expected to lead to decreasing heterosis effects. Nevertheless, it has to be considered that estimates of inbreeding effects may change over time, e.g., due to purging; thus, the effects for future generations may deviate from the extrapolation of the current estimates.
The breeding history of the German Brown is also important for the interpretation of ancestral inbreeding effects and purging. So far, there are only a few studies having analyzed the effects of ancestral inbreeding coefficients in dairy cattle populations, whereby longevity traits have rarely been under study (Supplementary Table S7). In German Brown cattle, the effects of ancestral inbreeding coefficients were negative and significant for lifetime and lifetime performance traits. None of the models accounting for ancestral inbreeding coefficients revealed evidence of purging. Considering the definition of Ballou (1997) [18], a significant favorable regression coefficient of the interaction of F with Fa_Bal would provide evidence for purging, while this study revealed significant negative effects for HL, LPL, NC, and lifetime production. Moreover, the effects of Fa_Bal were also negative when considered in the model as a regression coefficient alone. In Irish Holsteins, a positive, but not significant interaction of F with Fa_Bal was observed for survival from first to second lactation and 305-day production traits, while the effect of the simple regression model with Fa_Bal was significantly positive only for milk and protein yield [23].
The approach with Ahc assumes that alleles that have met more frequently in the past are less likely to have a detrimental effect than alleles that have never or only a few times before been identical-by-descent (IBD) [22]. Since the inbreeding coefficient indicates how often a randomly selected allele was IBD during pedigree segregation, an increasing Ahc could be associated with a beneficial effect on phenotype [6,22]. In Dutch Holsteins, positive effects for Ahc were reported for first lactation protein yield [6], whereas the current study showed only negative effects for Ahc on lifetime and lifetime production traits. In addition, inclusion of Ahc alone or simultaneously with F resulted in similar estimates for the regression coefficients of Ahc. We also tested the interaction of Ahc and F using model 2, which revealed comparable results to the interaction of Fa_Bal and F. In general, results for Fa_Bal and Ahc were very similar in this study with slightly higher estimates of Fa_Bal for some of the traits when included in the different models, which may indicate equal meaning of these ancestral inbreeding measures, at least for the current pedigree structure of German Brown cows, where correlation of Fa_Bal and Ahc was close to one [26,51].
According to Kalinowski et al. (2000) [21], successful purging would be present if FNew had a more detrimental effect on the trait compared to Fa_Kal. In contrast, in this study the effects of FNew were close to zero and not significant, while Fa_Kal had significantly negative effects on HL, LPL, NC, and lifetime production. These results indicate that the detrimental effects of inbreeding for longevity traits are mainly due to ancestral generations. In Irish Holsteins, a significantly negative effect of Fa_Kal was reported for survival to second lactation, while FNew was not significant [23]. For 305-day production traits, FNew was responsible for greater losses than Fa_Kal in Irish Holsteins [23]. Also, in Dutch Holsteins, FNew was significantly negatively associated with first lactation yield [6]. No significant effects of Fa_Kal and FNew were observed in Canadian Holsteins for lactation yield, but for milk and protein yield the effects of FNew were unfavorable while those from Fa_Kal were favorable [24].
Our results for German Brown using different models to account for ancestral inbreeding were consistent and did not show favorable effects of inbreeding from more ancient generations for lifetime and lifetime performance traits. Therefore, we assume that purging effects are very unlikely. However, the informative value of the different models used strongly depends on the structure of the pedigree. Deeper pedigrees are likely to reveal more inbreeding, which could also affect the results in terms of inbreeding depression and purging. The depth of the pedigree in this study was comparable to other studies, as we have discussed previously [19].
There may be a few reasons why purging was not detected in the German Brown population for longevity traits, but instead negative genetic contributions of the ancestral generations were discovered. First, the efficiency of purging depends on the structure of inbreeding within the population and the rate at which inbreeding increases [8] implying differences in the success due to different breeding and selecting histories. The overall level of inbreeding and ancestral inbreeding in the current population was not high due to introgression of US Brown Swiss bulls; thus, purging might possibly occur in future generations, where a further increase in inbreeding has been observed [26]. Additional analysis of survival traits using birth year cohorts from 2002 to 2008 and records through 2018 for Fa_Kal and FNew gave estimates for Surv1 of 0.0998 ± 0.3412 and −0.1168 ± 0.0522 with p-Values of 0.7698 and 0.0251, respectively, and for Surv3 estimates of 0.0159 ± 0.4229 and −0.11515 ± 0.0646 with p-Values of 0.9699 and 0.07488, respectively. Mean values and standard deviations for Fa_Kal and FNew were 0.00183 ± 0.00264 and 0.01569 ± 0.01695, respectively. These preliminary data suggest that significant purging effects are not yet present, even though positive effects for Fa_Kal were observed. On the other hand, estimates for F (Surv1: −0.2790 ± 0.0438, p-Value < 0.001 and Surv3: −0.5821 ± 0.0548, p-Value < 0.001) and FNew indicated larger negative effects due to inbreeding depression. A similar picture emerged for effective lifetime performance for milk, fat, and protein yield in this data set with positive non-significant estimates for Fa_Kal and negative non-significant estimates for FNew. Purging may not even counterbalance the negative effects of new inbreeding in these data from 2002 to 2008. Other opportunities for including survival in more recent data would be multiple trait animal models that only regard survival up to the first four lactation numbers, but split the first or each lactation into three different periods to account for the distribution of reasons for culling in the different lactation periods. These multiple trait animal models need to be tested to see whether they can be extended to cows surviving more than seven or nine lactations if records from the younger birth cohorts are incomplete [52,53,54].
Increase in inbreeding in the ancestral generations was mainly due to multiple use of a limited number of bulls. In the German Brown population introgression with US Brown Swiss bulls started in 1966 resulting in higher milk production at the cost of a steeper increase in inbreeding [26] and reduced lifetime [28]. Thus, selecting towards higher milk production favored the more intense use of US Brown Swiss genetics. Sufficient purging is characterized through an elimination of deleterious genes. Through the introduction of US Brown Swiss, more genes were introduced that negatively influenced longevity and inbreeding depression did not decline through purging and elimination of unfavorable genes associated with longer lifetime in more inbred animals. Thus, it is likely that the negative ancestral effects on longevity traits might predominantly stem from ancestors representing a higher breed proportion of US Brown Swiss.
Furthermore, the greatest success of purging has been found under high selection pressure [9]. In the German Brown, the productivity traits fat and protein yield and protein percentage were economically weighted with 48% until 2015, resulting in a high selection pressure for these traits. In contrast, economical weights for functional traits such as longevity and fertility were only 16.1% and 8.6%, indicating less intensive selection. Estimation of breeding values for longevity started in 1995, so these traits had less time to express purging. In addition, purging depends on the genetic structure of the trait [7]. Lower heritabilities for functional traits compared to production traits make selection less successful than for productivity traits. To date, ancestral inbreeding effects on fertility were mostly not significant or at least less conclusive considering purging [6,23,24]. Thus, assuming that both fertility and longevity are functional traits subject to lower selection pressure, our results are consistent with these previous findings.
Significant inbreeding depression was shown in German Brown cows for lifetime, lifetime production, and effective lifetime production traits. However, compared to the current level of inbreeding, the negative effects through inbreeding depression are counterbalanced by positive heterosis and recombination effects, leading to positive combined effects. We were not able to demonstrate purging as all models accounting for ancestral inbreeding did not reveal positive effects on lifetime, lifetime production, survival, and main culling reasons for Fa_Kal, Fa_Bal, Ahc, and F×Fa_Bal. Inbreeding depression was predominantly caused by ancestors from more ancient generations than inbreeding in newer generations. The intensive use of US Brown Swiss bulls in the first generations after immigration started may also have contributed these outcomes. In addition, longevity traits have been under less selection pressure at this time, thus positive effects due to selection against inbreeding depression are not likely to appear in the data of the present study. Thus, it is worth considering the dynamics of the different ancestral inbreeding coefficients in future generations. As ancestral inbreeding is increasing more strongly than new inbreeding in German Brown cows to recent birth years, purging effects for longevity can possibly appear in future generations. Our results indicate that a further decline of lifetime and lifetime production can be prevented when measures can be implemented to slow down the increase in ancestral inbreeding and the further increase in new inbreeding through US Brown Swiss genes. At present, our data have some limitations in terms of the depth of the pedigree and the endpoint of the data (2001), future studies should possibly use censored data and survival analysis methods or multiple trait animal models for survival to different parities and periods within the respective lactations to gain insights into the actual development of longevity traits.
We thank the official milk recording organization of Bavaria (Landeskuratorium der Erzeugerringe für tierische Veredelung in Bayern e.V., LKV) for providing the data and particularly Anton Neuhuber for preparing the data sets.
The following supporting information can be downloaded https://www.mdpi.com/article/10.3390/ani13172765/s1, Table S1: Means, standard deviations (SD), minima, maxima of classical inbreeding (F), inbreeding coefficient according to Ballou (1997) (Fa_Bal), ancestral history coefficient according to Baumung (2005) (Ahc), and ancient (Fa_Kal) and new (FNew) inbreeding coefficient according to Kalinowski (2002) for German Brown born between 1990 and 2001; Table S2: Pearson correlation coefficients of inbreeding coefficients, heterosis, and recombination coefficients for German Brown cows born between 1990 and 2001; Table S3: Inbreeding depression per 1% inbreeding expressed in percentage of the respective trait mean (x¯), and phenotypic (σp) and genetic (σa) standard deviation; Table S4: Expected effects of heterosis, recombination, and inbreeding for longevity traits in German Brown cows analyzed (birth years 1990 to 2001) and projected for German Brown cows of the birth year 2014; Table S5a: Expected phenotypes and differences between cows with a low (5% percentile) and a high (95% percentile) degree of inbreeding and differences between fully inbred and non-inbred cows for the classical inbreeding coefficient F; Table S5b: Expected phenotypes of lowly (5% percentile) and highly (95% percentile) inbred cows for Fa_Kal, FNew, and Ahc; Table S5c: Observed inbreeding coefficients (F, Fa_Kal, FNew, and Ahc) for survival to 2nd (Surv1), 4th (Surv3), 6th (Surv5), 8th (Surv7), and 10th (Surv9) lactation number and culling rate due to foot and leg problems (CulCL), infertility (CulINF), and udder diseases (CulUD); Table S6: Inbreeding coefficients of cows that survived the respective lactation number (1, 3, 5, 7, 9) compared with non-survivors (Surv1, Surv3, Surv5, Surv7, Surv9); Table S7: Studies analyzing classical (F) and ancestral inbreeding coefficients (Fa_Bal, Fa_Kal, FNew, and Ahc) with their means and standard deviations (in brackets) in Holstein dairy cattle.
Conceptualization, O.D. and A.W.; methodology, O.D. and A.W.; software, O.D.; validation, O.D. and A.W.; formal analysis, O.D. and A.W.; investigation, A.W. and O.D.; resources, O.D. and J.D.; data curation, O.D. and J.D.; writing—original draft preparation, A.W.; writing—review and editing, O.D. and J.D.; visualization, A.W.; supervision, O.D.; project administration, O.D. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
Restrictions apply to the availability of these data. Data were obtained from the Landeskuratorium der Erzeugerringe für tierische Veredelung in Bayern e.V. and are available from the authors with the permission of the Landeskuratorium der Erzeugerringe für tierische Veredelung in Bayern e.V.
The authors declare no conflict of interest and there are no relevant financial or non-financial competing interests to report.
This research received no external funding. This Open Access publication was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—491094227 “Open Access Publication Funding” and the University of Veterinary Medicine Hannover, Foundation.
Restrictions apply to the availability of these data. Data were obtained from the Landeskuratorium der Erzeugerringe für tierische Veredelung in Bayern e.V. and are available from the authors with the permission of the Landeskuratorium der Erzeugerringe für tierische Veredelung in Bayern e.V.