1. Introduction
German Brown cattle have been developed in the Alpine region in South Germany as a dual-purpose breed for butter and cheese production as well as for fattening of surplus calves to produce valuable carcasses. This breed is characterized by good quality legs and claws, so that the cows can withstand the harsh environmental conditions in this mountainous region. Their robust health and high fertility enable the cows to have high lifetime production and a long life [
1]. The breeding history of the German Brown and other European Brown populations is shaped by the introgression of US Brown Swiss bulls over the last 50 years [
2]. With increasing US Brown Swiss blood proportions, milk performance increased, but longevity and lifetime performance were highest with US Brown Swiss blood proportions of 41–50% [
3,
4]. The ancestral (inbred common ancestors in the pedigree) and new inbreeding (inbred animal in the pedigree for the first time) coefficients from pedigree data in the German Brown population show increasing trends, but the effects of inbreeding depression on longevity and lifetime production were counterbalanced by positive heterosis effects [
5,
6]. A further increase in inbreeding is expected to reduce the positive heterosis effects and negatively influence the lifetime production and cow survival rates [
6]. A decrease in genetic diversity and effective population size can have negative effects on adaptability to a changing environment [
5].
An evaluation of inbreeding and its effects on the genome architecture in the actual German Brown population can be obtained from genomic data. A popular method used in genomic studies is the analysis of runs of homozygosity (ROH) [
7,
8,
9,
10,
11,
12,
13,
14,
15]. They are assumed to be the result of the transmission of identical haplotypes from parents to offspring and are considered as an indicator of the degree of autozygosity. Since recombination during meiosis leads to shortening of ROH segments, it is assumed that short ROH are the result of inbreeding and selection [
16]. The distribution and frequency of ROH of different lengths is population-specific and depends on the selection intensity and direction to which the respective population was exposed in its breeding history [
7,
15]. Furthermore, genomic regions that contain a high percentage of ROH (ROH islands) are regions subject to strong selection, as previously shown by overlapping the results of the analysis of selection signatures and ROH islands [
16,
17], and they may harbor genes involved in phenotypic characteristics of a breed [
12,
18,
19]. Further analysis of these regions can therefore shed light on the nature of selection within a breed and reveal loci associated with economically important traits [
20].
For European Brown Swiss populations, genome wide association studies (GWAS) have identified regions associated with production, udder morphology, fertility, calving and birth, body conformation and carcass traits based on different numbers of animals and SNPs, using deregressed breeding values of bulls as phenotypes for the analysis (
Supplementary Material Tables S1 and S2) [
21,
22,
23,
24,
25,
26,
27]. Few studies have distinguished between Original Brown breeds and modern Brown Swiss breeds that contain a proportion of US Brown Swiss blood [
25,
27], but it has been shown that modern Brown populations are one of the most differentiated breeds, with a strong genetic divergence from the Original Brown breeds [
28]. It can therefore be expected that different selection objectives between Original Brown breeds and modern Brown breeds have targeted different genomic regions, resulting in different genomic localizations and the sizes of the ROH and ROH islands. Studies of ROH islands and selection signatures in Brown Swiss populations have been conducted on US Brown Swiss, Italian, German, Austrian and Swiss Brown and Original Brown, mostly based on medium density SNP data, with 27 K–62 K SNPs (
Supplementary Material Table S3).
Nevertheless, differences between Brown Swiss and Original Brown populations were evident between and within studies. In Brown Swiss populations, the region on BTA 6 between 80–95 Mb and the region on BTA 5 between 75 and 80 Mb were frequently identified as ROH islands or selection signatures (
Supplementary Material Table S4), harboring genes associated with milk production traits, udder health and longevity. The QTL around 90 Mb was identified as a target of ongoing selection in Brown Swiss [
23], confirming the strong focus of selection on production and udder conformation in Brown Swiss populations (
Supplementary Material Table S4). In contrast, the selection pressure in Original Brown populations was more on BTA 11, as ROH islands [
18] and selection signatures [
29,
30,
31] were detected between 65 and 73 Mb with genes related to fertility, fat deposition, meat quality, adaption and immune response, and on BTA 26 between 21 and 23 Mb with genes such as
FGF8, associated with meat and carcass quality (
Supplementary Material Table S5). A QTL for body size and leg conformation was already found in this region in Brown Swiss [
23].
To our knowledge, no study to date has analyzed the ROH structure in a large number of randomly selected German Brown cows and bulls with special attention to the breed proportion of US Brown Swiss and the association with survival as an important functional trait for this breed. Therefore, the aim of this study was to characterize the ROH structure and ROH-based inbreeding coefficients to evaluate the breeding history in the context of selection and inbreeding in 2364 randomly selected German Brown cattle and to compare these results with previous reports. In addition, we aimed to analyze common ROH and ROH islands to show the effects of selection on the genomic architecture in the German Brown. To better understand the influence of US Brown Swiss sires in German Brown cattle, we also performed these analyses for different classes of US Brown Swiss breed proportions and survival rates to higher lactation numbers.
4. Discussion
In this study, we aimed to characterize the patterns of ROH in German Brown cattle based on 50 K SNP data with special emphasis on the influence of the breed proportion of US Brown Swiss and survival rates for the 1st to 9th lactation in order to find genomic regions associated with survivability to high lactation numbers and whether the introgression of US Brown Swiss did effect survival rates.
The effective population size calculated in this study based on LD showed similar patterns to Italian Holsteins, where an increase from generation 5 to the youngest generation was reported [
10], but the effective population size was lower for Italian Holsteins with estimates of 96 and 120 for generations 6–9 and 120 in the youngest generation, respectively. A lower N
e was also found for Italian Brown bulls, namely 237.6 and 90.7 for generation 50 and 5, respectively [
50] and for Canadian Holsteins and Jerseys with N
e estimates of 58 and 120 for generation 5, respectively [
51]. In agreement with results for Italian Modeneses [
13] and US Holsteins and Jerseys [
51], N
e estimates using regression on time resulted in lower values. Even smaller N
e estimates were obtained from the realized rate of inbreeding. Nevertheless, the N
e estimates exceeded the critical value of 50, at which long term negative effects due to inbreeding depression may be expected [
5,
10,
13,
51].
When comparing different studies on ROH analyses, different ROH definitions complicated the interpretation. In previous studies analyzing ROH in Brown populations, the minimum number of SNPs ranged from 15–60, with a minimum length between 1–2 Mb, allowing either one or no heterozygous or 0–2 missing SNPs [
17,
18,
19,
29,
31,
49,
52]. In addition, in most studies, a pruning for MAF was performed. As this could lead to an underestimation of genomic inbreeding, this was not considered in this study, possibly leading to higher inbreeding estimates [
35,
37,
52]. The minimum number of SNPs and the minimum length in this study were calculated based on the number of SNPs, the number of animals, and the average heterozygosity, and are thus higher than most of the studies, being most consistent with those in the study on US Brown Swiss [
52]. In agreement with previous studies [
35,
37,
52], it may be expected that in our data, when applying a lower minimum number of SNPs, a lower minimum length of ROH, a larger gap size between consecutive ROH, a smaller window size and without pruning for MAF, more ROH will be detected, resulting in higher genomic inbreeding coefficients (
Supplementary Material Table S7).
In German Browns, almost half of the ROH were between 4 Mb and 8 Mb long, which are assumed to correspond to 12.5 and 6.25 ancestral generations ago, assuming 1 cM = 1 Mb. In Original and modern Brown populations, the most abundant ROH class was 4–8 Mb [
18], and in Swiss Brown populations 5–10 Mb [
29]. In Italian Brown Swiss [
19,
50] and US Brown Swiss [
52], a decreasing number of ROH was found with increasing ROH length class, so that most ROH were 1–2 Mb and 2–4 Mb respectively. Since the minimum ROH length in the present study was set to 3.191 Mb, ROH of 1 to 3 Mb were not considered in this study. Using the ROH settings of the US study [
52] but without MAF pruning, we observed in our data an increase of ROH segments with 3–4 Mb and a shift of the frequency from 17.04% to 26.90% for the ROH class < 4 Mb across all animals. The number of ROH < 4 Mb by US Brown Swiss classes increased from 5.9–6.5 to 10.6–11.5 in the present study.
Nevertheless, the average length of a ROH of 8.323 Mb is comparable to the average length for Italian Browns of 8.54 Mb [
18]. In US Brown Swiss, the average ROH length was slightly lower at 7.54 Mb, with a higher average number of ROH at 47.69 [
52]. Applying the definitions for ROH like in the US study [
52] but without MAF pruning, the average ROH length decreased to 7.657 Mb and the average number of ROH increased to 40.709 in our data. Therefore, the reason for this difference may be mainly due to the definition of ROH and, to a lesser extent, that more detected ROH in US Brown Swiss belong to shorter ROH length classes, accounting for inbreeding in former times. Therefore, short ROH < 3.191 Mb, which were not retrieved in our study, led to a larger number of ROH of smaller sizes in US Brown Swiss. Nevertheless, studies comparing ROH of Original Browns to modern brown breeds reported a lower average ROH length [
18] and a lower average number of ROH [
29] for the Original Browns. Although we did not study animals of different Brown cattle subpopulations, the evaluation according to different US Brown Swiss classes showed a similar trend for the distribution of ROH by length classes. A higher average number of ROH per animal was found for animals in BS 90–99% compared to animals of <60% and 70–89% breed proportion of US Brown Swiss. However, ROH islands of the 99th percentile threshold were not specific for US Brown Swiss classes. All four ROH islands of the 99th percentile threshold found for the entire sample were also present in each of the US Brown Swiss classes. Nevertheless, a greater differentiation of the genomic architecture between US Brown Swiss classes was found when searching for ROH islands of the 95th percentile threshold. When comparing the 18 ROH islands of the 95th percentile threshold identified in the entire sample, 7/18 ROH islands of the 95th percentile threshold were equally distributed across all US Brown Swiss classes and 5/18 across four US Brown Swiss classes. Within US Brown Swiss classes, a total number of 43 ROH islands was identified. There were 7/43 ROH islands common to all US Brown Swiss classes, 11/43 ROH islands common to BS <70% and 9/43 ROH islands frequent in BS >69%. In summary, the distribution of the ROH islands of the 99th percentile threshold indicates changes in the size of the ROH islands between US Brown Swiss classes, but not to such an extent that animals in BS 90–99% can be considered a distinct subpopulation. Even 67% of the 18 ROH islands of the 95th percentile threshold, as shown for the entire sample, are shared by all five or at least four US Brown Swiss classes. Nevertheless, ROH analyses within US Brown Swiss classes suggested that animals of the different US Brown Swiss classes differ in 10/43 ROH islands of the 95th percentile threshold and there may be trait loci, which may be specific for traits of the animals in these respective US Brown Swiss classes (
Supplementary Material Table S20). Most distinct ROH islands of the 95th percentile threshold were found in BS <60% (7/10) and only 1/10 in BS 90–99%.
The proportion of animals sired by US Brown Swiss bulls amounted to 4.47% in the German Brown birth cohorts from 1990–2018. German Brown sires are therefore mostly crosses with US Brown Swiss bulls. The 75% and 95% confidence intervals for the proportion of US Brown Swiss genes in the sires for the 1990–2018 birth cohorts of the German Brown population were 69–88% and 58–97%, respectively. This may explain why the ROH islands of the 95th percentile threshold were found in all five US Brown Swiss classes and 12/18 of the ROH islands of the 95th percentile threshold were found in at least four US Brown Swiss classes.
The average F
ROH in the present study was higher than what has been reported for Original Brown, Braunvieh (Brown cattle with incrossings of US Brown Swiss bulls) and US Brown Swiss sires in Switzerland with F
ROH of 0.029, 0.074 and 0.091, respectively [
29]. Higher values for F
ROH of 0.15 were found for US Brown Swiss in the US [
52]. Using the same settings for defining ROH as in the latter two studies, the same tendency for F
ROH was confirmed (
Supplementary Material Table S7).
The increasing trend in F
ROH and pedigree based inbreeding coefficients across birth years for the 2364 animals in this study is consistent with our recent study in German Browns based on pedigree data, in which the inbreeding coefficients increased from 0.013 to 0.036 between 1990 and 2014 [
5]. The average values for lifetime traits and pedigree-based inbreeding coefficients from the present data set also agree well with the entire German Brown population (
Supplementary Material Table S6). The animals in this study therefore represent a random sample from the entire population.
Based on the pedigree data, the relationship of F
New and F
PED shows that more than two thirds of the inbreeding comes from alleles that are IBD for the first time and the high correlations between F
New and F
PED underline the large influence of inbreeding from more recent generations. The correlation between F
New and F
ROH was rather similar to that of F
ROH>4 and F
ROH>8, while the correlations between F
New and F
ROH>16 or F
ROH>32 were considerably smaller. Considering the correlation between F
ROH>8 and F
ROH with an estimate of 0.936, which corresponds to the correlation coefficient between F
New and F
PED, it can be assumed that a large portion of new inbreeding is covered by ROH segments between 8–16 Mb, i.e., the ROH length class with the highest total length. This relates to 3.125 to 6.25 generations. Underlying a generation interval of 7 years, this corresponds to a period of 22 to 44 years and falls into the time after the beginning of the introgression of US Brown Swiss in 1966, i.e., 24 to 52 years ago in relation to the birth years 1990–2018 in this study. Along this line, correlation coefficients between the different F
ROH and proportion of US Brown Swiss genes decreased with longer ROH segments. Therefore, it is likely that ROH < 6.25 Mb originated from Original Brown ancestors prior to the introduction of US Brown Swiss. These results support the outcomes of our previous study: that the breeding scheme during that time increased the level of inbreeding due to the unbalanced use of a few top US Brown Swiss sires and contributed to the current level of inbreeding [
5]. Nevertheless, systematic mating of closely related individuals cannot have taken place on a large scale as the F
IS value is negative and close to zero.
The tendency toward higher F
ROH with increasing breed proportion of US Brown Swiss supports the impact of US Brown Swiss on inbreeding. However, the lack of differences of F
ROH>16 and F
ROH>32 between the US Brown Swiss classes and their correlations close to zero also indicate that in more recent generations the effect of using US Brown Swiss genetics has lost its impact on inbreeding levels in the German Brown population. The fact that not all animals showed ROH of more than 16 Mb could also explain the lack of significant differences. In studies comparing Original Browns with Brown Swiss, higher F
ROH values were also reported for the modern Brown Swiss populations [
28,
29].
Testing for significant differences is influenced by samples sizes (
n) since sample size is connected with the size of standard errors (
, with SD = standard deviation). This may have had an influence on the non-significant
p-values when comparing average numbers of ROH, F
IS, F
ROH and F
ROH>4 between BS <60%, with a small sample size of 70, and BS 70–79%. A larger sample size for BS <60% at >200 would have resulted in significant
p-values. Rather uneven sample sizes are also common in previous reports on ROH [
17,
18,
19,
20,
28,
29,
30].
It has already been shown that the proportion of US Brown Swiss genes is associated with longevity [
3], and that longevity is negatively associated with inbreeding [
6]. Considering survival to the following calving, cows that left the herd earlier had significantly higher F
ROH values compared to those surviving until the following calving. This means that cows that are more inbred have a higher risk of leaving the herd. Cows that survived to at least the 5th to the 9th lactation exhibited the lowest levels of inbreeding. This relation is consistent with the results based on pedigree data [
6]. When accounting for the effect of the interaction of US Brown Swiss classes with survival, the differences between the survivors and not survivors were largest for F
ROH and BS <60% but decreasing in the BS classes with more US Brown Swiss genetics (
Supplementary Material Table S13). We may propose that inbreeding may be more detrimental for survival to high lactation numbers in German Brown cows when the proportion of US Brown Swiss is decreasing, particularly <60%. Recent inbreeding should go to zero independent of the US Brown Swiss proportion in order to increase the possibility of long survival times of the cows. Genomic inbreeding was lowest for cows with BS <60% having survived 5, 7 and 9 lactations and went to zero for F
ROH>32.
The region of strongest selection was detected on BTA 6 between 73 and 91 Mb, which was the longest ROH island and the longest consensus ROH shared by approximately 30% of the animals, with the region between 85 and 88 Mb shared by approximately 50% of the animals. In this genomic region, the casein genes are located, which are central for dairy cattle breeding. The casein genes
CSN1S1,
CSN2,
CSN1S2 and
CSN3 have been associated with milk production [
23] and milk protein percentage [
53] in European Brown breeds and other dairy cattle populations [
23,
53,
54]. The
CSN3 gene was within the ROH shared by 50% of the animals and encodes kappa casein, which is particularly important for the coagulation properties required for cheese production [
55]. Also, a QTL for milking speed was identified in this region [
24].
Furthermore, in Brown cattle populations, genome wide association studies have found this region to be associated with mammary gland morphology [
23], particularly for udder length and teat diameter [
21]. The region between 85 and 88 Mb has also been associated with genes involved in udder health and morphology in different dairy breeds such as
SLC4A4,
NPFFR2, GC and
RASSF6 [
56,
57,
58,
59] and longevity (
NPFFR2 and
ADAMTS3) [
60,
61,
62], both important traits for German Brown cattle breeding. In agreement with our study, this region was also most frequently reported in studies analyzing ROH islands or selection signatures in modern Brown cattle populations. A selection signature in German Browns using the Cross Population Extended Haplotype Homozygosity (XP-EHH) was detected between 84 and 96 Mb [
30]. Also, in Italian Brown Swiss, where selection signatures were identified for SNPs with a ROH count per SNP >50%, a signature around 85 Mb was detected [
19]. In a study comparing Original Brown and modern Brown breeds in Italy, a ROH island, defined as the top 0.999 SNPs of the percentile distribution, was detected on BTA 6 between 86 and 87 Mb for modern Brown breeds [
18]. In US Brown Swiss, BTA 6 was also reported as a ROH island [
52].
Also, the region on BTA5 between 71 and 86 Mb was discovered in previous studies as a ROH island in US Brown Swiss and Italian Brown Swiss [
18,
29,
52]. This region harbors QTLs [
63,
64] and genes primarily associated with milk production such as
NCF4,
RAC2 and
CSF2RB [
65,
66], and is therefore also a region of particular interest in dairy cows.
The ROH island on BTA 16 in the German Brown overlapped the smaller ROH island found in US Brown Swiss [
52], suggesting that it originated from US Brown Swiss. This is also supported by the decline of this ROH island with decreasing breed proportion of US Brown Swiss, where for BS <60% this ROH island comprises only 5 SNPs, whereas for BS 90–99% this ROH island reaches its largest size. Moreover, this region overlaps with QTLs mainly associated with milk production and includes genes related to fertility (
TGFB2) [
67] and udder morphology (
SUSD4) [
68]. In addition, the region between 21 and 31 Mb was reported to harbor a selection signature in US Brown Swiss and Swiss Brown in Switzerland [
29]. The region on BTA 5 between 12 and 24 Mb was uniquely herein discovered in BS 60–69%, including e.g., the
TMTC2 gene, which was identified as a candidate gene for udder morphology in Brown cattle populations [
23].
Nevertheless, the aim of this study was to assess to which extent Original German Brown genes are still under selection pressure. In studies on Original Browns, the region on BTA 11 around 65 to 71 Mb was found to be under strong selection [
29,
30,
31]. This typical Original Brown signature was also detected in the present study as a consensus ROH shared by 30% of the animals at 64 to 69 Mb containing genes associated with meat quality such as
CAPN14 and
CAPN13 and fertility (
PROKR1) [
18,
30]. This may suggest that the Original Brown genes are still fixed to a lower extent within the actual German Brown population. Interestingly, when defining ROH islands as the 95th percentile to account for the relatively small number of animals in low US Brown Swiss classes, the region on BTA 11 was detected as a ROH island with nearly the same length across all US Brown Swiss classes. Thus, the homozygosity originating from the Original German Brown animals has been dissolved in about 70% of the actual German Brown population due to the incrossing of US Brown Swiss in a former time, which is consistent with the development from a dual-purpose cattle breed to a more milk-emphasized cattle breed in terms of the associated traits mentioned above. When comparing the ROH islands of the 95th percentile threshold between US Brown Swiss classes, we found that 11 ROH islands, which were present in BS <70%, disappeared in BS >69%, while 5 new ROH islands developed in BS >79%.
Expanding the ROH-definition to the 95th percentile revealed another ROH island on BTA 2 between 76 and 81 Mb which was only detected in BS <60%. This region harbors genes such as
GYPC and
CNTNAP5 (
Supplementary Material Table S16b–d) that have been associated with growth and carcass traits in Beninese indigenous cattle as well as meat quality in Hanwoo [
69,
70] and milk fat composition [
71], representing traits that are more important in beef and dual-purpose cattle breeds than in pure dairy breeds and are therefore likely to be of Original Brown origin. In addition, this region has been reported as a ROH island in an old Polish dual-purpose cattle breed, the Polish Black-and-White [
72].
Two further ROH islands on BTA 22 (12–31 Mb) and 28 (15–23 Mb) were detected at BS <60%, BS 60–69% and BS 70–79%, decreasing in length with increasing US Brown Swiss class and absent in animals of the highest US Brown Swiss classes (
Supplementary Material Table S16b–d). In Original Browns, SNPs around 24 Mb were associated with multiple births, but not in Brown Swiss [
27]. On the other hand, this ROH island was located near the
MITF gene, which was associated with a white spotted coat color in Brown Swiss [
73].
The ROH island on BTA 28 has also been found in one of the oldest Polish cattle breeds, the White-Backed breed [
72]. Furthermore, this region has been associated with various traits in dairy cattle. Traits such as kappa casein percentage [
74], claw health [
75,
76], and limb and claw conformation, are of particular interest in German Brown cattle breeding [
68,
77]. Interestingly, the regions on BTA 2, 22 and 28 have also been associated with longevity [
68,
78,
79]. In Chinese Holsteins, a genome-wide association study on nine longevity traits identified the region around 22.9 Mb on BTA 28, which includes
CTNNA3, as associated with a productive life span from first calving until the end of the first and second lactation [
78]. Thus, increasing breed proportions of US Brown Swiss also increase heterozygosity in a region that is likely to affect longevity.
Increasing heterozygosity was also observed in animals in higher lactation numbers, as the ROH islands based on the 99th percentile on BTA 6 and 16 were only found in animals surviving the third lactation, but not in cows surviving higher lactations.
The definition based on the 95th percentile revealed a ROH island on BTA 7 between 42 and 45 Mb only for cows surviving lactation 9. Since a SNP in this BTA 7 region at position 43,904,171 was associated with herd life in Chinese Holsteins [
78], this may indicate that this region may be important for longevity in German Browns. Interestingly, this ROH island on BTA 7 was detected in animals with BS <69% but not in animals with higher proportions of US Brown Swiss genes, supporting the previously reported negative correlation of breed proportion of US Brown Swiss and longevity [
3]. The ROH island on BTA 7 overlaps with previously reported QTLs for milk yield [
80] and milking speed in French dairy cattle [
81], stature and body conformation in Canadian Holstein bulls [
82], lean meat yield [
83] and somatic cell scores [
84].
Thus, further investigations of the regions on BTA 7 and 28 with regard to longevity traits in German Brown cows seem necessary to clarify the importance of these regions for longevity traits and to enable the adaption of breeding programs, to prevent further decline of longevity, as an outstanding trait in German Brown cows.