A Comprehensive Analysis of Runs of Homozygosity of Eleven Cattle Breeds Representing Different Production Types

Simple Summary Runs of homozygosity (ROH) regions are known to be common in the genomes of cattle and have become a subject of interest of various research in recent years. ROH can be used as a valuable tool to estimate inbreeding, which needs to be controlled in livestock populations. Moreover, analysis of ROH is considered to be an effective method of identifying genome regions that are a subject of selection pressure, which may help in understanding the genetic aspects of production traits under selection. In this study, we analyzed ROH characteristics of 11 cattle breeds, both commercial and native, maintained in Poland. We presented distinct differences in the length, quantity and frequency of ROH between the analyzed breeds as well as in the levels of genomic inbreeding. Higher levels of inbreeding were characteristic for commercial breeds, especially beef breeds. In addition, within ROH islands, we observed a number of genes with a confirmed influence on the level of production traits. The presented results and identified genes can be a basis for further research focused on the identification of genes and markers essential in the determination of the most important production traits in cattle. Abstract In the presented research, BovineSNP50 microarrays (Illumina) were applied to determine runs of homozygosity in the genomes of 11 cattle breeds maintained in Poland. These cattle breeds represent three basic utility types: milk, meat and dual purpose. Analysis of runs of homozygosity allowed the evaluation of the level of autozygosity within each breed in order to calculate the genomic inbreeding coefficient (FROH), as well as to identify regions of the genome with a high frequency of ROH occurrence, which may reflect traces of directional selectin left in their genomes. Visible differences in the length and distribution of runs of homozygosity in the genomes of the analyzed cattle breeds have been observed. The highest mean number and mean sums of lengths of runs of homozygosity were characteristic for Hereford cattle and intermediate for the Holstein-Friesian Black-and-White variety, Holstein-Friesian Red-and-White variety, Simmental, Limousin, Montbeliarde and Charolais breeds. However, lower values were observed for cattle of conserved breeds. Moreover, the selected livestock differed in the level of inbreeding estimated using the FROH coefficient. In regions of the genome with a high frequency of ROH occurrence, which may reflect the impact of directional selection, a number of genes were observed that can be potentially related to the production traits which are under selection pressure for specific production types. The most important detected genes were GHR, MSTN, DGAT1, FABP4, and TRH, with a known influence on the milk and meat traits of the studied cattle breeds.


Introduction
Runs of homozygosity (ROH) can be defined as long homozygous chromosomal regions within which both haplotypes inherited from parents are identical. The identity of large chromosome segments inherited from both parents by the offspring may be the result of their descent from a common ancestor and point to a certain level of parental relatedness [1]. Over the past 20 years, there has been a rapid development of high-throughput genomic analysis methods, including next-generation sequencing (NGS) and genotype-based microarrays (SNPs). Both of these methods allow for effective identification of runs of homozygosity [2]. However, SNP microarrays remain the most popular tool for determining ROH, mainly due to their lower costs and less time-consuming analysis compared to NGS.
ROH can be used as a valuable tool to estimate inbreeding, which is one of the methods often used in animal husbandry utilized to preserve characteristics of outstanding individuals in the population [3]. The inbreeding coefficient determined using ROH (called F ROH ) is calculated as the ratio of the total length of ROH for each individual in the selected ROH length category to the total length of the autosomal chromosomes covered by SNPs [4]. The determination of the ROH-based inbreeding coefficient shows several advantages compared to the classical inbreeding coefficient calculated on the basis of pedigree data (F PED ). F ROH more accurately predicts the actual level of autozygosity of the genome and it can be estimated in any animal with genotypic data-even if information about its pedigree is not available [1,3,5].
An increase in the homozygosity of certain regions of the genome may also take place as a result of intense artificial selection, which leads to the increase of frequency of beneficial alleles in the population in a process similar to genetic drift [6][7][8]. Identification of such regions that show a reduction or elimination of polymorphisms (known as selective sweeps) may indicate the occurrence of directional selection and help in understanding the genetic aspects of selected production traits [9,10]. Analysis of the frequency of alleles and genomic homozygosity are considered to be effective methods of identifying genome regions that are a subject of selection pressure [10,11]. Therefore, it is assumed that the size and location of ROH may be correlated with the effect of selection pressure and that ROH should occur non-randomly in the genome [12]. Moreover, it has been shown that ROH can create unique patterns in the genomes of different breeds, probably under the influence of natural or artificial selection [13,14]. This information can be potentially used in evolutionary research as well as in gene mapping related to production traits [15]. In the research of Nothnagel et al. [13] and Pemberton et al. [14], it has been shown that ROH are not evenly distributed throughout the genome and form so-called "ROH islands". These regions are widely observed in populations and form specific patterns within the genome. These patterns can be used as a useful tool to identify the phenomenon called "selective sweeps" and genome regions subject to selection pressure [10][11][12]16].
The purpose of this study was to characterize runs of homozygosity and identify ROH patterns in the genomes of 11 selected cattle breeds (Holstein-Friesian Black-and-White variety (HO), Holstein-Friesian Red-and-White variety (RW), Simmental (SM), Limusin (LM), Hereford (HH), Charolais (CH), Montbeliarde (MO), White-Backed (BG), Polish Red (RP), Polish Red-and-White (ZR) and Polish black and white (ZB)) maintained in Poland. As part of the study, genome regions were identified that are characterized by a high incidence of ROH, i.e., potentially under the influence of directional selection. Moreover, the structure of homozygosity sequences and their patterns was determined in the examined cattle breeds, which will contribute to a better understanding of genetic differences between the breeds, their breeding history and changes in the genome under the influence of intensive selection. In addition, the aim of this study was to assess the inbred of the studied animal populations based on the inbreeding coefficient calculated based on homozygosity sequences (F ROH ). This is particularly important in the case of animals belonging to herds of conserved breeds, for which pedigree information is often shallow and the level of inbreeding is crucial for ongoing programs for the protection of genetic resources.

Research Material, Genotyping and Data Filtering
The research material was genomic DNA obtained from 1931 randomly selected cows and bulls belonging to 11 breeds of cattle kept in Poland. For meat and milk production type breeds, the material was randomly sampled from Biological Material Bank of National Institute of Animal Production which included samples shipped by breeders for parentage verification. For conserved breeds, animals were randomly sampled from small individual farms located in various regions of Poland among the animals included in conservation program. In both cases, animals were controlled not to include closely related individuals. The detailed characteristics of cattle regarding the use type and the number of individuals of a given breed are presented in Table 1 and additional information regarding the conserved breeds is presented in Supplementary Material S1. Some of the animals from four of the selected breeds (HO, RP, LM and SM) were previously analyzed in our earlier publication [17]. However, the results described in our previous publication were maintained on a different SNP dataset (different filtration of microarrays) and the ROH identification was maintained in a different way. Moreover, all the animal procedures were approved by the Local Animal Care Ethics Committee No. II in Kraków-permission number 1293/2016 in accordance with EU regulations.
DNA was isolated from semen, ear tissue or whole blood using Sherlock AX (A & A Biotechnology, Gdynia, Poland) or QuickGene DNA whole blood kit S (KURABO, Kurashiki, Japan) kits. The purity and concentration of the obtained DNA was determined using a NanoDrop2000 spectrophotometer. High-quality DNA was normalized to the required concentration (50 ng/µL) and analyzed using Illumina BovineSNP50 v2 BeadChip microarrays (Illumina Inc., San Diego, CA, USA). Finally, the microarrays were scanned using the HiScanSQ system (Illumina). All the procedures were carried out in accordance with the manufacturer's protocol.
Filtering of the genotypic data obtained after scanning was performed together for all the tested breeds while ROH identification was performed separately for each examined breed. The dataset contained 54,609 SNPs before filtering. Only animals with more than 95% genotypes (CallRate) and SNP with GenCall quality coefficients above 0.7 and GenTrain above 0.4 were used for further analysis. Next, SNPs on the X, Y chromosome and those without a fixed genomic position (mapped to contigs) were removed. The final marker panel included probes for 42717 SNP markers localized in the genome with a mean distance of 55.4 kb (±44.4 kb).

Identification of ROH
The SNP panel obtained after filtering was used to identify ROH for each animal individually using cgaTOH software [18]. The following parameters were used for the identification of ROH: a minimum number of 30 consecutive homozygous SNPs in ROH with a maximum distance between SNPs equal to 1 Mb. The identified ROHs were assigned to five length categories: 1-2 Mb, 2-4 Mb, 4-8 Mb, 8-16 Mb and above 16 Mb. To calculate the number of heterozygotes to be included in the identification for each ROH length category, ROH was initially determined without allowing any heterozygotes. This enabled the calculation of the number of SNPs for each category of ROH length. Based on this data, assuming a 0.2% genotyping error for the Illumina microarrays [19], 0 heterozygotes were assigned in the first four ROH length categories (1-2 Mb, 2-4 Mb, 4-8 Mb, 8-16 Mb) and 1 heterozygote for the last ROH category with a size above 16 Mb. The number of missing genotypes that were allowed in ROH was calculated according to the methodology proposed by Ferenčaković et al. [20] and the following criteria were applied: in the categories of ROH 1-2 Mb and 2-4 Mb, no missing SNP genotypes were allowed, in the category 4-8 Mb, one missing SNP, in the category 8-16 Mb, two missing SNPs, while in the category over 16 Mb, four missing SNPs were allowed. The average sums of ROH in selected categories were calculated by summing all ROHs identified for each animal in each category and averaging the results within the breed.

Determination of Inbreeding Coefficient Based on ROH
The F ROH inbreeding coefficient was calculated according to the methodology proposed by McQuillan et al. [4] by dividing the total ROH length for each individual in the selected ROH length category by the total length of autosomal chromosomes covered by SNPs (2510.6 Mb). The ROH length categories were as follows: >1 Mb, >2 Mb, >4 Mb, >8 Mb and >16 Mb.
The literature data show that F ROH calculated for ROH above 1 Mb is the most frequently used measure of genomic inbreeding and best describes recent animal relatedness. However, due to some errors in the identification of short ROH observed in the case of 54 K microarrays [1], F ROH for ROH with length above 4 Mb was also calculated to compare the obtained coefficients and eliminate errors resulting from the use of short ROH segments for calculations.

Identification of Genome Regions with a High Frequency of ROH Occurrence
To identify genome regions characterized by a high frequency of ROH occurrence, result files determining how many times each SNP appeared in a given population in ROH were generated for each population. Subsequently, 1% of the highest occurrence values were chosen as the threshold value above which the identified genome regions were classified as "ROH islands" in a given population. These regions were analyzed for overlapped genes using the UCSC Genome Browser tool [21] based on the UMD3.1 bovine genome assembly and later, using the Panther Classification System [22] to identify their molecular functions and related biological processes. In addition, the size, range and frequency of the occurrence of ROH between the studied breeds was compared and the genomic ROH patterns were checked for similarities between the individual breeds and separate production types.

Characteristics of Runs of Homozygosity
One of the most important parameters characterizing the structure of ROH is the number and average sum of lengths of ROH per animal. The highest mean number of ROH per animal in individual populations were observed for the Hereford breed (80.6), while the correspondingly lower mean values were characteristic for conserved breeds: White-Backed (23.9), Polish Red (23.3) and Polish Red-and-White (21.8). The Holstein-Friesian Black-and-White variety, Holstein-Friesian Red-and-White variety, Simmental, Limousin, Montbeliarde and Charolais were characterized by a slightly higher number of ROH relative to conserved breeds with an average number of ROH in the range of 38.8 to 53.3 per animal.
The same trend was observed in the case of average sums of ROH lengths for individual animals. The highest average ROH lengths were observed for Hereford breed (378.0 Mb) and the lowest for The highest sum of ROH lengths was observed for the Simmental individual (874 Mb). Animals with a high total length of ROH were also observed in the following breeds: Charolais (759 Mb), Limousin (757.9 Mb) and the Holstein-Friesian Black-and-White variety (745 Mb). Basic statistics on the length and number of ROH are described in Table 2 and are shown in Figures 1 and 2. The highest sum of ROH lengths was observed for the Simmental individual (874 Mb). Animals with a high total length of ROH were also observed in the following breeds: Charolais (759 Mb), Limousin (757.9 Mb) and the Holstein-Friesian Black-and-White variety (745 Mb). Basic statistics on the length and number of ROH are described in Table 2 and are shown in Figures 1 and 2.    The highest sum of ROH lengths was observed for the Simmental individual (874 Mb). Animals with a high total length of ROH were also observed in the following breeds: Charolais (759 Mb), Limousin (757.9 Mb) and the Holstein-Friesian Black-and-White variety (745 Mb). Basic statistics on the length and number of ROH are described in Table 2 and are shown in Figures 1 and 2.

Analysis of the F ROH Inbreeding Coefficient
One of the most reliable inbreeding coefficients describing the phenomenon of inbreeding resulting from both past and recent relationship is F ROH coefficient calculated for all identified ROHs, i.e., those with a length above 1 Mb.
On the one hand, the highest mean F ROH values calculated for ROH above 1 Mb were observed for the Hereford breed  Table 3 and are shown in Figure 3.

Characteristics of Genomic Regions with a High Frequency of ROH Occurrence
In order to identify genomic regions with the highest frequency of ROH and thus potentially under the influence of selection, ROH was analyzed on all 29 autosomes of 11 breeds of cattle maintained in Poland. Next, we calculated the frequency of individual SNPs that were present in ROH and selected the top 1% of the markers most commonly occurring in ROH. The neighboring markers with the highest frequency of ROH occurrence pointed to regions of the genome in which there are likely to be haplotypes that are a subject of selective pressure. The achieved results are presented in detail in Table 4 and shown schematically in Figure 4.
The analyzes allowed for the identification of 5 to 17 separate regions of the genome with a high frequency of ROH occurrence for particular breeds. The highest number of such regions was detected in the Polish Black-and-White breed (17) and the lowest in the Charolais breed (5). These regions were located on 5 (Charolaise) to 14 (Polish Red-and-White) autosomal chromosomes and had the length in range of 77.6 kb (Polish Red-and-White, three consecutive SNPs) to 1258.3 kb (Limousin; 206 consecutive SNPs). The average length of the region calculated for all breeds was 24,111.2 kb, while the average number of SNP per region was 38.5.

Characteristics of Genomic Regions with a High Frequency of ROH Occurrence
In order to identify genomic regions with the highest frequency of ROH and thus potentially under the influence of selection, ROH was analyzed on all 29 autosomes of 11 breeds of cattle maintained in Poland. Next, we calculated the frequency of individual SNPs that were present in ROH and selected the top 1% of the markers most commonly occurring in ROH. The neighboring markers with the highest frequency of ROH occurrence pointed to regions of the genome in which there are likely to be haplotypes that are a subject of selective pressure. The achieved results are presented in detail in Table 4 and shown schematically in Figure 4.     The analyzes allowed for the identification of 5 to 17 separate regions of the genome with a high frequency of ROH occurrence for particular breeds. The highest number of such regions was detected in the Polish Black-and-White breed (17) and the lowest in the Charolais breed (5). These regions were located on 5 (Charolaise) to 14 (Polish Red-and-White) autosomal chromosomes and had the length in range of 77.6 kb (Polish Red-and-White, three consecutive SNPs) to 1258.3 kb (Limousin; 206 consecutive SNPs). The average length of the region calculated for all breeds was 24,111.2 kb, while the average number of SNP per region was 38.5.
In the identified regions with a high ROH occurrence, from 98 (LM) to 250 (RW) genes were detected depending on the breed. The number of genes located within ROH islands and detected in multiple breeds is shown in Figure 5 and the list of genes is presented in Supplementary Material S2.
The highest number of genes located within ROH islands and detected in multiple breeds was observed for MO and SM (60), HO and RW (49) as well as BG and ZB (51) breeds. Within dairy and meat production types ( Figure 5B), no such genes for all breeds were observed, while genes located within ROH islands and detected in multiple breeds with the most similar constitutional type (HO and RW, MO and SM) were observed. In the case of conserved breeds ( Figure 5C), there were 13 such genes between the four studied cattle breeds.
Genes located in regions with the highest frequency of ROH in individual populations were associated with numerous biological processes. The results of gene classification without enrichment analysis of individual processes are presented in Table 5. In the studied breeds, the largest number of genes were involved in cellular processes (from 42 to 109 genes) and metabolic processes (from 33 to 76 genes). Characteristic for conserved breeds, compared to other breeds, was a higher proportion of genes involved in processes related to the functioning of the immune system, reproduction, cell regulation and the organization of cellular components or biogenesis.  Venn diagram showing the number of genes located within ROH islands and detected between studied cattle breeds with division into breeds and cattle production types. (A) milk production type, (B) meat production type, (C) dual purpose production type. Table 5. Biological processes most often associated with genes identified within genome regions with a high frequency of ROH occurrence. Among the identified genes, a significant part were involved in the pathways presented in Table 6. In most breeds, genes associated with the pathways responsible for integrin signaling were identified. Characteristic pathways for dairy breeds were those associated with angiogenesis and signaling of thyrotropin and gonadotropin releasing receptor. In the case of conserved dual-purpose breeds, pathways characterized by the highest number of identified genes were associated with the activation of B and T cells and inflammatory processes involving the cytokine and chemokine signaling pathway.

HO RW SM LM HH CH MO BG RP ZB ZR
Within the regions of the genome with a high frequency of ROH occurrence, a number of genes with a confirmed influence on the level of production features have also been identified, including DGAT1, MSTN, FABP4, ERBB3I, STAT1, GHR, SUFU, BTRC or CHUK.

Characteristics of Runs of Homozygosity
The characteristics of the ROH detected in this study corresponds well to the results obtained by other authors. Purfield et al. [1] identified ROH for various breeds of cattle, including European breeds: Holstein-Friesian, Limousin, Simmental, Hereford and Charolais. The authors showed that the average sum of ROH length for segments longer than 5 Mb and identified using BovineHD BeadChip microarrays was the highest for the Hereford (145 Mb) and Holstein-Friesian (115 Mb) breeds. However, the Limousin, Simmental and Charolais breeds were characterized by much lower values, respectively 45 Mb, 55 Mb and 50 Mb. In the case of all ROH (>1 Mb), the average sum of the ROH lengths was respectively higher: Hereford-245 Mb; Holstein-Frisian-195 Mb; Limousin-68 Mb; Simmental-85 Mb; and Charolais-85 Mb. Similar results were also obtained by Iacolina et al. [23], who compared European cattle breeds to the European bison using Illumina BovineHD BeadChip microarrays. The authors observed the highest average ROH lengths for Angus and Hereford breeds for the 1-5 Mb category (approximately 200 Mb). In the case of the Holstein-Friesian, Charolais, Simmental and Limousin breeds, these values were lower-in the range of approximately 100 Mb. In turn, in the research of Peripolli et al. [24] performed on Gyr cattle (Bos indicus), it was shown that the average sum of ROH lengths was in the range of 100 Mb per animal and the longest region was 108.9 Mb. These results are comparable with the results described in the present study for ROH with lengths >1 Mb and >4 Mb. In this study, the highest average ROH length was observed for the Hereford breed (378 Mb for ROH above 1 Mb and 253.4 Mb above 4 Mb), followed by the Holstein-Friesian Black-and-White variety (295.1 and 220.8 MB), and later, respectively, the following breeds: Charolais (163.2 Mb and 97.2 Mb), Simmental (169.1 and 87.4 Mb) and Limousin (147.4 Mb for ROH above 1 Mb and 80.3 Mb above 4 Mb). Apart from analyzing the ROH of high-production cattle breeds kept in Poland (HO, RW, SM, LM, MO, CH and HH), which are largely genetically close to their European counterparts, Polish indigenous breeds of cattle were also examined (BG, RP, ZR and ZB). These breeds are included in the genetic resources conservation programs. An analysis of ROH showed notable differences in the length and amount of ROH between native breeds and production breeds. Regarding the average number of ROH per animal for native breeds, it was shown to be in the range from 21 to 30, while for production breeds, it was in the range from 39 to 81. A similar trend can be observed for the average sum of ROH lengths-in the case of native breeds, it was lower and ranged from 105 to 135 Mb, while in the case of production breeds, it was in the range from 147 to 378 Mb.
According to the results obtained, we hypothesized that much higher mean length values and ROH numbers obtained for highly selected breeds in comparison to native breeds may be associated with strong directional selection and widespread and intensive use of artificial insemination. It has been proven that the use of biotechnics and intensive artificial selection significantly reduce genetic variability and increase the degree of relationship between animals [25]. In addition, native breeds were characterized by low, compared to commercial breeds, average amounts of ROH with a length of more than 8 Mb (2.2-3.8 in the case of conserved breeds and 3.1-10.1 in production breeds). Due to the fact that the occurrence of long segments of ROH is related to the recent relationship within the parental population [26], small amounts of long ROH segments indicate a low degree of the close relationship between animals from conserved herds. Similar differences between native breeds and production breeds were also observed by other authors [1]. The implementation of the genetic resources conservation program itself (Program for the Conservation of Genetic Resources of Farm Animals), which focuses on minimizing inbreeding within protected populations, seems to be of great importance for the obtained results. Not without significance is the admission in the protected populations of a certain share of the blood of other breeds, which increases the genetic variation and diversity of segregating haplotypes.
Differences in the ROH statistics between highly selected and native breeds may also result from the intensity of selection used within these populations. Kim et al. [27] studied Holstein cattle from populations with varying intensity of selection and observed that in animals undergoing intensive selection, the average length of ROH per animal was about 6.67 Mb, while for extensively selected animals, it was 6.26 Mb. In the case of an average amount of ROH per animal, the values were as follows: 40.4 for intensively selected animals and 31.1 for extensively selected animals. The same trend can be seen in the results of the present research by comparing intensively the selected production breeds and native breeds covered by the genetic resources conservation program.

Runs of Homozygosity as a Tool of Inbreeding Estimation
F ROH (inbreeding coefficient calculated on the basis of ROH) is widely regarded as a reliable measure of individual autozygosity and provides information on the degree of inbreeding of individual herds of animals, taking into account both past and recent relatedness of individuals [3,28]. The literature data also indicated that F ROH may be a better estimator of individual autozygosity than those based on pedigree data (F PED ) [28]. However, it should be noted that the use of medium density SNP arrays, such as the one used in this research, leads to limited identification of short ROH that contribute significantly to ancient inbreeding [17].
In many studies on cattle, a strong or moderate correlation between F PED and F ROH values was observed: 0.73 for ROH longer than 1 Mb and 0.70 for ROH longer than 10 Mb [1] and from 0.485 to 0.715 for ROH longer than 1 Mb [28]. In our previous study [29] performed on Polish Holstein-Friesian cattle, slightly lower values of the correlation coefficients between F ROH and F PED were observed (in the range from 0.308 to 0.505 for ROH with a length of more than 1 Mb). This is mainly the result of using Spearman's rank correlation coefficients, more appropriate for data sets with a large number of outliers and those not showing a normal distribution [30], but often giving lower correlation values.
It is worth mentioning that the correlation between F ROH and F PED coefficients rises with the increase in the length of ROH segments used for calculations. According to Marras et al. [26], this is due to the fact that ROH reflects both past and recent animal relatedness, while F PED coefficients are based on pedigree records that may not contain information from many past generations. Therefore, F ROH coefficient calculated on the basis of long ROH better reflects the recent relationship, and the F ROH -F PED correlation is usually higher considering only long ROH segments in the calculation. In studies conducted on pigs, Saura et al. [31] showed that the mean value of F ROH coefficient calculated for ROH with lengths >5 Mb was close to F PED , while the mean for F ROH calculated for ROH longer than 5 Mb was about four times lower than the average F PED . Another study was presented by Scraggs et al. [32], which suggested that F PED does not determine true kinship within the Wagyu cattle population, as there were clear differences between F ROH and F PED coefficients. The authors showed much lower values of F PED coefficients compared to F ROH . These results are consistent with data obtained by other authors for cattle [26,33] and pigs [31], in which the F ROH coefficient was characterized by higher values than F PED , suggesting that F PED may underestimate the degree of inbreeding of the studied populations.
Due to the information presented above, in this research, thes F ROH coefficient was used to assess the level of inbreeding of the selected cattle population maintained in Poland. The highest mean values of inbreeding coefficient calculated for ROH with a length above 1 Mb were observed in the Hereford (0.151), Holstein-Friesian Black-and-White variety (0.118) and Montbeliarde (0.108). The lowest F ROH values, as expected, were noted for native breeds: White-Backed (0.051), Polish Red (0.053), Polish Red-and-White (0.042) and Polish Black-and-White (0.054). In addition, mean F ROH coefficients calculated for ROH longer than 8 and 16 Mb were the highest in Montbeliarde, Holstein-Friesian Black-and-White and Hereford breeds, which suggests a relatively high degree of relatedness between animals and a reduction in genetic variation.
The inbreeding level estimated in this work for individual cattle breeds generally corresponds to the results obtained by other authors. In the study of Ferenčaković et al. [28], F ROH inbreeding coefficients for ROH with lengths above 1 Mb was calculated and high average autozygosity levels (0.151) were found for strongly selected Brown Swiss cattle population and significantly lower for native cattle (0.052 for Pinzgauer cattle and 0.066 for Tyrol Gray). In addition, in another study, Ferenčaković et al. [20] observed varying F ROH values, depending on the cattle breed: 0.156 for brown Swiss cattle; 0.088 for Fleckvieh cattle; 0.099 for Norwegian Red cattle; 0.087 for Tyrol Gray cattle and 0.09 for Simmental cattle. In the case of Holstein cattle maintained in the USA, F ROH calculated for ROH with a length of more than 5 Mb was low and amounted to 0.038 [34]. In Holstein cattle kept in Europe, higher F ROH coefficients were found in the range from 0.081 for ROH > 1 Mb to 0.046 for ROH > 5 Mb [1]. F ROH values for the Polish Holstein-Friesian cattle described in this paper was slightly higher than that observed for other populations of Holstein-Friesian cattle in Europe; however, similar to the other high-production European dairy breeds [1].

Analysis of Genomic Regions with a High Frequency of ROH Occurrence
The examination of runs of homozygosity distribution across the genome shows that their dispersal in particular regions of the genome is unique and they can form characteristic patterns depending on a given population. These genomic regions with the highest frequency of ROH occurrence were called by Nothnagel et al. [13] and Pemberton et al. [14] with a term "ROH islands". In studies of many authors, it is suggested that these regions are shaped by the influence of strong selection pressure on variants located at a given locus [10][11][12]14]. In addition, Zhang et al. [25] confirmed a significant correlation between short regions with a high frequency of occurrence of ROH and genomic regions considered to be under the influence of selection using methods based on F ST coefficient and iHS (integrated haplotype score).
An overview of the available literature indicates that many authors have observed and described the occurrence of ROH islands and their relationship with QTL for important production traits in cattle. Purfield et al. [1] observed that the genomic regions located on BTA 7, BTA 14, BTA 16 and BTA 18 were characterized by a high frequency of ROH occurrence and included important genes related to immune traits, muscularity and ease of calving. In particular, BTA 5 and BTA 9 were characterized by an increased number of long ROH above 20 Mb, within which there were numerous QTL regions associated with the production of milk fat and growth characteristics of cattle [1].
In this study, several genomic regions with extremely high levels of autozygosity and frequency of ROH were observed in individual breeds of cattle. An example of such a region may be the initial sequence of chromosome 2 in the Limousin breed, within which there is a myostatin gene locus (MSTN), which is a strong QTL for muscle traits in this breed [35]. A similar region on BTA 2 was found by Marras et al. [26] by studying Piedmontese cattle, where almost 90% of all animals were characterized by the occurrence of ROH in the proximal part of the second chromosome. Other regions of interest, potentially influenced by the directional selection, may be two regions located on chromosome 5 and 6 in the Charolais breed, within which genes related to growth factors and coat color are located. One of these genes is ERBB3-encoding the epidermal growth factor receptor. The literature showed that ERRB3 controls the proliferation and myogenic diversity of muscle stem cells [36]. Then, there are the INHBC and INHBE genes that are part of the TGFB signaling pathway-a transforming growth factor beta, which is an essential regulator of fibroblast proliferation, collagen synthesis in muscle tendons and in the case of muscular tissue can affect muscle atrophy [37]. Therefore, these genes may potentially be associated with meat traits of Charolais cattle that are subject to selection in this breed. Moreover, in various studies [38] within this ROH island, the PMEL gene was proposed as the one associated with the white coat color of Charolais cattle. Another region with high homozygosity was observed in the proximal part of chromosome 14. This region was characterized by a high frequency of occurrence of ROH, in the range of 60% to 80% in the population of LM and CH breeds (respectively); however, elevated levels of homozygosity were also observed for this region in other breeds of cattle: RW, BG, RP, ZB and ZR. The identified region included the DGAT1 gene locus, with known influence on fat percentage in milk and FABP4 gene associated in numerous previous studies with transport, regulation and lipid metabolism [39]. In addition, the FABP4 gene was presented as affecting the qualitative characteristics of meat in cattle and also affecting fat content in meat [40]. In turn, in the studies of Zhou et al. [41], the effect of the FABP4 gene on the amount of milk produced and protein content in milk was also observed.
When comparing the results of this research to literature data, a large similarity of the location of ROH islands identified in different breeds of cattle can be observed. In studies of Mészáros et al. [42], a ROH island was observed in Tyrol Gray cattle located on chromosome 6 in the 36-41 Mb position. In the same chromosomal region, in this study, regions with an increased incidence of ROH were identified for the Montbeliarde, Simmental and Charolais breeds. The research performed on domestic cattle allowed the detection of four ROH islands, two of which were also identified in this study: the first on BTA 6 in the location of 38.2-39.4 Mb, which coincides with the ROH islands observed in the MO, SM and CH breeds; the second on BTA 16 with the location of 43.8-45.0 Mb, which coincides with the regions with high ROH frequency for the MO, SM, HH, RP and ZB breeds [43]. In the studies of Marras et al. [26], the authors observed characteristic ROH islands on chromosomes 2 and 6. The ROH island on chromosome 2 was characteristic to Piemontese cattle and was associated with the MSTN gene. A similar effect resulting from strong selection pressure on the MSTN gene was observed in this research in Limousin cattle. In contrast, ROH island located on chromosome 6 in Holstein and Italian Brown dairy cattle included ABCG2 and FAM13A1 genes, which are known to be associated with milk traits [44,45]. The same genes were identified in regions with a high incidence of ROH in this study in the following cattle breeds: SM, CH, MO (ABCG2 gene) and HO, RW, SM, CH, MO and ZB (FAM13A1 gene).
In the present study, in regions with a high frequency of ROH occurrence, a number of other genes associated with cattle production traits were also observed. An example of such a gene identified in the ROH island on chromosome 2 is the STAT1 gene observed in ZB cattle. Cobanoglu et al. [46]'s study demonstrated the important role of this gene in the regulation process of transcription of other genes involved in the synthesis of milk proteins and fat metabolism. In dairy cattle breeds HO and RW, within ROH islands, the GHR gene was observed, which has a significant impact on the quantity and composition of milk [47,48]. In addition, the same breeds have demonstrated the potential effect of selection pressure on the TRH gene belonging to the thyrotropin secretion signaling pathway, which has a significant effect on growth hormone secretion and prolactin [49]. Literature data confirm that exogenous stimulation of these hormones promotes a significant increase in milk yield [50,51]. Further examples of genes identified in ROH islands in the RW, BG, ZR and ZB breeds are SUFU and BTRC, that regulate the proliferation of lactate stem cells. In humans, the SUFU gene acts as a negative regulator of the hedgehog signaling pathway, whereas the BTRC gene controls the self-renewal process of the mammary gland stem cells [52]. In BG, ZR and ZB breeds, in the ROH islands, we also detected the CHUK gene, which, in the literature, is described as one of the genes responsible for the development of mammary gland and lactation [53]. In addition, genes participating in the Wnt signaling pathway were identified in RW, SM, BG and RP cattle breeds. The Wnt signaling pathway is responsible for the regulation of cells proliferation, their morphology, apoptosis processes and differentiation as well as the development of the mammary gland [53]. Another group of genes identified in the ROH islands are genes related to the signaling pathway of the angiogenesis process observed in the following breeds: HO, SM, MO, RP and HH, which may also be associated with the development of the mammary gland and affect its blood supply [54]. In addition, in the case of conserved breeds, the highest number of genes identified in the ROH islands belonged to the metabolic pathways associated with the activation of B, T cells and inflammatory processes involving cytokine and the chemokine signaling pathway, which may be related to the higher resistance to diseases characteristic of these breeds.
When analyzing genes occurring in ROH islands and observed in more than one breed, 13 genes were found between BG, RP, ZB and ZR cattle breeds. These genes have been described previously as affecting growth in cattle and humans (TGS1, LYN, CHCHD7, SDR16C5, TMEM68) [55][56][57][58], ease of calving (RPS20) [59] and traits related to puberty (FAM110B) [60]. These traits are also subject of selection in cattle of different breeds and may have been fixed during a long-term breeding process.

Conclusions
To conclude, in the present study we examined the characteristics of ROH identified on the basis of genomes of 11 cattle breeds. The whole genome approach allowed the presentation of visible differences in the length, quantity and frequency of ROH between the analyzed breeds. These differences were particularly pronounced between commercial and conserved breeds. Moreover, we calculated the levels of genomic inbreeding within individual cattle breeds, which showed a higher level of inbreeding characteristic to commercial breeds, especially meat breeds. In addition, regions of the genome that were most commonly associated with ROH were identified, which may reveal signatures of recent selection and describe ROH patterns of the analysed cattle breeds. Within these genomic regions, a number of genes with a confirmed influence on the level of production traits were observed, including DGAT1, MSTN, FABP4, ERBB3, STAT1, GHR, SUFU, BTRC or CHUK. The presented results and identified genes can be a basis to a further research focused on the identification of genes and markers essential in determining the most important production traits in cattle.