Skip to main content
  • Research Article
  • Open access
  • Published:

Runs of homozygosity in the Italian goat breeds: impact of management practices in low-input systems

Abstract

Background

Climate and farming systems, several of which are considered as low-input agricultural systems, vary between goat populations from Northern and Southern Italy and have led to different management practices. These processes have impacted genome shaping in terms of inbreeding and regions under selection and resulted in differences between the northern and southern populations. Both inbreeding and signatures of selection can be pinpointed by the analysis of runs of homozygosity (ROH), which provides useful information to assist the management of this species in different rural areas.

Results

We analyzed the ROH distribution and inbreeding (FROH) in 902 goats from the Italian Goat Consortium2 dataset. We evaluated the differences in individual ROH number and length between goat breeds from Northern (NRD) and Central-southern (CSD) Italy. Then, we identified the signatures of selection that differentiate these two groups using three methods: ROH, ΔROH, and averaged FST. ROH analyses showed that some Italian goat breeds have a lower inbreeding coefficient, which is attributable to their management and history. ROH are longer in breeds that are undergoing non-optimal management or with small population size. In several small breeds, the ROH length classes are balanced, reflecting more accurate mating planning. The differences in climate and management between the NRD and CSD groups have resulted in different ROH lengths and numbers: the NRD populations bred in isolated valleys present more and shorter ROH segments, while the CSD populations have fewer and longer ROH, which is likely due to the fact that they have undergone more admixture events during the horizontal transhumance practice followed by a more recent standardization. We identified four genes within signatures of selection on chromosome 11 related to fertility in the NRD group, and 23 genes on chromosomes 5 and 6 related to growth in the CSD group. Finally, we identified 17 genes on chromosome 12 related to environmental adaptation and body size with high homozygosity in both groups.

Conclusions

These results show how different management practices have impacted the level of genomic inbreeding in two Italian goat groups and could be useful to assist management in a low-input system while safeguarding the diversity of small populations.

Background

Today, in light of the ongoing climate change, the management and conservation of livestock biodiversity are becoming an increasingly important goal at the global level [1]. To face this challenge, it is crucial to draw a precise picture of the genetic structure of the indigenous breeds and populations of farmed animals in different countries. It is necessary to understand the genetic basis of their adaptation, not only to the natural environment, but also to the breeding conditions and management strategies to which they have been subjected [2]. From this point of view, Italy provides a good model because it is characterized by a rich biodiversity in all domesticated species thanks to its varied history, environment, climate, and farming traditions [3, 4]. Goats, in particular, represent one of the greatest expressions of Italian biodiversity with more than 30 autochthonous breeds and populations reared under very diverse climates and farming conditions, several of which are considered as low-input agricultural systems [5].

In the Northern regions, goats are mainly bred in the Alps, where two diametrically opposed farming systems coexist. On the one hand, in the valleys and hilly regions, modern intensive and semi-intensive farming systems are present, which are specifically suited to milk and cheese production and usually exploit cosmopolitan dairy goat breeds, particularly Saanen and Alpine [6, 7]. In these systems, medium-to-large flocks are mostly kept indoor with controlled feeding and limited grazing, which is generally conducted in fenced pastures near the farm [8]. On the other hand, the traditional extensive farms, which can be considered as low-input/low-output systems, are mainly located in the mountainous areas and depend highly on natural grazing. On these farms, small flocks of local breeds are kept indoor during the winter and in pasture for the rest of the year because of extreme variations in climate and weather conditions, especially during the winter. Some farmers still practice the traditional vertical transhumance (alpeggio), which consists in transferring the animals to alpine pastures during the summer only [6,7,8]. The animals of this farming system are particularly influenced by the climate conditions.

Central-southern Italy and the islands count the largest number of goat farms and heads [9]. In these regions, which are characterized by a hotter and dryer climate [5] than in Northern Italy, the traditional extensive or semi-extensive farms of autochthonous goat breeds prevail and are generally located in marginal mountainous areas [10, 11]. The vertical transhumance was usually combined with a horizontal transhumance; for example, the shepherds transferred all their animals—cattle, sheep, goats, and shepherd dogs—in the mild Apulian plains during the winter and returned to the Abruzzo mountains in the summer [12].

These differences between Northern and Southern Italian goat populations in terms of animal nutrition, housing, and mating management, may have contributed to the genetic makeup of the Italian caprine diversity. Among the genomic tools and approaches that are now proposed to characterize animal biodiversity, the analysis of runs of homozygosity (ROH) is certainly one of the most useful [13]. ROH are long stretches of homozygous genotypes in the genome of an individual, which compose a pair of identical haplotypes. They are considered as a standard approach for the calculation of genomic inbreeding values (FROH) and for the detection of signatures of selection [14]. The length of a ROH can also be a useful indicator of the time of the inbreeding event with which it is associated, i.e. long ROH are associated with recent events of inbreeding in the history of a breed or of a single individual, whereas short ROH indicate a more ancient event [15]. The presence of several ROH in a particular region of the genome of a species or a population, regardless of their length, constitutes a so-called ROH island. The analysis of ROH islands can be a very effective tool to identify the regions of a genome that have been under selective pressure because they can contain variants that are shared between the individuals of a specific population [16]. For all these reasons, the analysis of ROH from genomic data and the derived inbreeding value (FROH) are increasingly used as a starting point to develop new management systems of animal populations, together with the more traditional pedigree information [17].

In this work, our aim was to characterize ROH in 902 goats from the Northern and Southern Italian groups, estimate their level of inbreeding, and analyze how it has evolved across generations according to management practices.

Methods

Dataset and quality control

In this work, we used the same Italian Goat Consortium2 (IGC2) dataset as described in Cortellari et al. [5]. Among the 34 populations present in that dataset, we decided to exclude the Bezoar, which in the previous work was used as an outgroup, the Maltese × Sarda crosses and the two Montecristo populations due to their unique history of isolation (feral) or farming (mainland). All the animals were genotyped with the Illumina Goat single nucleotide polymorphism (SNP)50 BeadChip. SNPs that had a missing genotype frequency higher than 0.2, or that were in unplaced scaffolds or on the X chromosome were excluded from the analysis, but we did not apply a threshold for minor allele frequency (MAF) to better identify ROH [18]. Individuals with a call rate lower than 95% were removed. All quality control procedures were carried out with the software PLINK 1.9 [19]. After the initial quality check on the 986 individuals and the 52,538 SNPs taken from the original IGC2 dataset, 902 animals grouped in 30 breeds and 46,995 SNPs were retained. Population structure of the goats included in the final dataset was investigated by multidimensional scaling analysis (MDS) and by building a phylogeny tree based on Reynolds genetic distances.

Expected heterozygosity (genetic diversity)

For each breed, the PLINK 1.9 software (-hardy option) was used to calculate the expected heterozygosity (HE), observed heterozygosity (HO), and Wright’s fixation index (FIS), which is defined as the correlation between the homologous alleles within an individual relative to the local population to which that individual belongs [20].

Definition of runs of homozygosity

In order to minimize the discovery of false positives within regions of low marker density, we selected rather stringent criteria [18]. ROH were calculated separately for each individual using the software PLINK 1.9 by applying a sliding window of 20 SNPs. A ROH was called if the following parameters were fulfilled: (i) no heterozygous genotypes, (ii) less than two missing genotypes, (iii) a minimum number of SNPs within a ROH ≥ 20, (iv) a minimum ROH length of more than 1 Mb, (v) a minimum SNP density of two SNPs per Mb, and (vi) a maximum gap of 500 kb between consecutive homozygous SNPs.

ROH distribution and genomic inbreeding

To characterize the ROH distribution, for each breed we estimated: the number of individuals without ROH, the mean number of ROH per individual, the mean total length of ROH per individual, the mean length of a ROH per individual, and the genomic inbreeding coefficient (FROH) for each individual. The FROH for each breed was computed following the method proposed by McQuillan [21]:

$${F}_{ROH, i} = \frac{{L}_{ROH} }{{L}_{AUTO}},$$

where \({L}_{ROH}\) is the sum of the total length of ROH in individual \(i\) and \({L}_{AUTO}\) is the total length of the autosomes covered by SNPs. In addition, we categorized the ROH for each breed into five length classes (1–2 Mb, 2–4 Mb, 4–8 Mb, 8–16 Mb, and > 16 Mb) to compare the distribution of the FROH across these categories between the considered breeds [22]. We focused on these length classes with the intent to investigate the percentage and the impact of ancient and more recent inbreeding events that occurred in the Italian goat breeds.

Identification of the groups of populations

In order to better disentangle the genetic differences between the Italian goat populations analyzed due to climatic conditions and breeding management techniques, we divided them into two large groups according to their geographical distribution: a group of ten populations from the Northern Italy breeds (NRD) and a group of 20 populations from the Central-southern Italy breeds (CSD), which also includes the two Maltese populations (MAL and SAM).

Statistical analysis

We performed two linear mixed models to evaluate the statistical significance of the difference between the ROH parameters identified for each population group [23] using the statistical software JMP 16 [24]. We modelized two variables (\(\mathrm{Y}\)): (i) a standardized ratio between the length of single ROH in each individual and the length of the corresponding chromosome \(\bigg(std\left(\frac{ROH\,length}{CHR\,length}\right)\bigg)\), and (ii) a standardized ratio between the number of ROH on each chromosome of each individual and the length of the chromosome \(\bigg(std\left(\frac{ROH\,number}{CHR\,length}\right)\bigg)\); both of these models included the same factors:

$$\mathrm{Y}=\upmu +\mathrm{CHR}+\mathrm{POPGROUP}*\mathrm{CHR}+\mathrm{BREED}\left[\mathrm{POPGROUP}\right]+\mathrm{POPGROUP}+\mathrm{id}+\mathrm{e},$$

where \(\upmu\) is the mean, \(\mathrm{CHR}\) is the fixed effect of the autosome (chromosomes 1 to 29), \(\mathrm{POPGROUP}\) is the fixed effect of the population group (CSD vs NRD), \(\mathrm{BREED}[\mathrm{POPGROUP}]\) is the fixed effect of the breed nested within the population groups (n = 20 in CSD and n = 10 in NRD), \(\mathrm{POPGROUP}*\mathrm{CHR}\) is the interaction between population group and autosome, \(\mathrm{id}\) is the random effect of the animal, and \(\mathrm{e}\) is the random residual. The covariance between the animals was assumed to be equal to 0.

Signatures of selection

In this work, we investigated the signatures of selection by using three methods: ROH, ΔROH islands and the Wright’s fixation index (FST).

For the first method, a homozygosity score (H-score) ranging from 0 (0%) to 1 (100%) was obtained for each SNP by counting how many times it appeared in a ROH and dividing the result by the number of the animals. This approach was applied to each population separately, and the top 1% H-scores were considered. SNPs that were within regions of 0.25 Mb were joined together, and regions with more than 15 SNPs were considered as ROH islands. Then, the identified ROH islands were investigated to list the annotated genes they contain based on the reference genome (ARS1) and the associated functions and pathways.

The ΔROH score was defined as the difference between the H-scores for the CSD and the NRD groups at a specific SNP. The regions of maximum difference in homozygosity were defined analogously to the ROH islands (top 1% values, SNPs within a region of 0.25 Mb combined together, and regions encompassing > 15 SNPs), thus resulting in ΔROH islands.

The Wright’s fixation index (FST) was calculated using the PLINK software, averaging the value of each SNP with the values of five adjacent SNPs in each flanking region to minimize the impact of outlier scores [22]. The top 1% averaged FST values were considered and investigated for annotated genes in the reference genome (ARS1) within a region spanning ± 0.25 Mb from each SNP.

Finally, the genes identified by these three methods were analyzed to detect the shared genes.

Results

Dataset composition and quality control

The dataset used for the analyses consisted of 902 goats belonging to 30 populations. The number of analyzed animals per breed are in Table 1. The geographic distributions, the MDS plot, and the phylogeny tree of the studied goat populations are reported in Additional file 1: Fig. S1.

Table 1 Dataset composition and mean ROH and genetic parameters for the studied Italian goat breeds

ROH description and genetic diversity

In total, 28,383 ROH were identified in the 902 individuals considered; five animals displayed no ROH: one Garganica individual, one Rossa Mediterranea individual, one Grigia Ciociara individual, and two Capra di Teramo individuals. In terms of average number of ROH per animal, the breed with the largest number was Orobica (92.8 ROH per individual), followed by Vallesana (76.6), both of these breeds being raised in the northern regions of Italy. The breeds with the smallest average number of ROH were two Sicilian breeds, Messinese and Argentata dell’Etna with 9.3 and 9.1 ROH per individual, respectively. Argentata dell’Etna and Messinese were also the breeds that, together with the Val di Livo (or Lariana) had the lowest average value of the total ROH length per individual (31.2, 36.8, and 56.8 Mb, respectively), while the two highest values were found in Vallesana (364.6 Mb) and Maltese bred in Sicily (347.5 Mb). When the average length of ROH per individual in each breed was considered, the Capra di Livo-Lariana, Argentata dell'Etna, and Orobica were the breeds with the lowest values, whereas the Capra di Teramo and Roccaverano were those with the highest values (Table 1).

The ROH-based inbreeding values (FROH) showed that the two breeds with the highest level of inbreeding were Maltese and Vallesana and those with the lowest levels were Messinese, Argentata dell’Etna, and Capra di Livo-Lariana. However, the distribution of the individual FROH within each population varied among breeds (see the boxplot in Fig. 1): some breeds such as Maltese, Capra di Teramo, Vallesana, and Girgentana showed a wide dispersion of the individual inbreeding values, while other breeds such as Saanen, Rossa Mediterranea, Capra di Livo-Lariana, and Messinese showed a more compact distribution. The genomic diversity parameters were similar across all the breeds, with the lowest HE and HO respectively in the Orobica (0.35) and Vallesana (0.35), and the highest ones in Roccaverano (0.41) and Saanen (0.41). The FIS highest values were found for Roccaverano and Derivata di Siria, and the lowest for the Jonica and Rossa Mediterranea breeds (Table 1).

Fig. 1
figure 1

Distribution of FROH values per individual and mean values per breed. Boxplot (a) of the single individual FROH distribution in each population (Northern breeds: blue boxplots, Central-southern breeds: red boxplots) with matching barplot (b) of the mean FROH values (each color representative of the different ROH length classes)

The ROH identified in all the populations were classified into the five length classes. The largest numbers of identified ROH belonged to the 1–2 Mb (11,294) and the 2–4 Mb (7525) classes. Nevertheless, the length classes that contributed most to the calculated inbreeding value (FROH) in the different populations were 8–16 Mb and > 16 Mb. Analysis of the distribution of the ROH categories among all the populations and of their proportional weight in the definition of the mean FROH value revealed three categories: (1) one for which the influence of the longest ROH on the total FROH value appeared to predominate, such as for the Capra di Teramo, Roccaverano, and Montefalcone breeds; (2) one for which the different ROH length classes were well balanced, such as for the Girgentana, Bianca Monticellana, and Nera di Verzasca breeds; and (3) one for which the short ROH were more important, such as for the Orobica, Capra di Livo-Lariana, and Rossa Mediterranea breeds (Fig. 1).

Statistical analysis

The two groups of Italian goat populations included 560 individuals belonging to 20 breeds for the Central-southern (CSD) and 342 individuals belonging to ten breeds for the Northern (NRD) groups. The statistical models showed a significant difference in the number and length of ROH between the two groups. Particularly, the first model (r2 = 0.18) showed that ROH length was significantly affected by the group (p = 0.003), the breeds within each group (p < 0.0001), and the chromosome (p < 0.0001). ROH were longer for the CSD than the NRD group (LSmean ± SE = 0.06 ± 0.02 vs − 0.02 ± 0.02). Figure 2a shows the mean standardized ratio between the length of ROH and the length of the corresponding chromosome \(std\left(\frac{ROH\,length}{CHR\,length}\right)\) for each chromosome in the two groups (CSD and NRD): the largest differences were found for chromosomes 3, 13, 25, and 29 and the smallest differences for chromosomes 8, 18, and 20. Interestingly, the smaller chromosomes presented relatively longer ROH in both groups. Previous studies on plants, yeasts, and humans have shown that recombination rates are inversely correlated with chromosome length, which could be due to the lower frequency of multiple crossovers within a chromosome [25,26,27]; this has also been reported in goats [28] and other species [29].

Fig. 2
figure 2

Comparison of mean standardized length and number of ROH in the two groups of Italian goat populations. Graphic representation of the mean standardized length (a) and number (b) of ROH divided by the corresponding chromosome length in the two groups of Italian goat populations

When the individual number of ROH per chromosome was modelled (r2 = 0.40), all the selected factors (group, breed within each group, chromosome, and the interaction between chromosome and group) were significant (p < 0.0001). ROH number was, on average (± SE), larger for the NRD group (0.19 ± 0.03) than the CSD group (− 0.07 ± 0.02). Figure 2b shows the mean standardized ratio between the number of ROH per chromosome and the corresponding chromosome length \(std\left(\frac{ROH\,number}{CHR\,length}\right)\) in the two groups (CSD and NRD): the largest differences were found for chromosome 13 and the smallest for chromosomes 8, 18, 20, and 27. It is worth mentioning that the number of ROH for chromosome 12 was large in both groups, particularly in the CSD group.

Signatures of selection

For both the CSD and NRD groups, we identified the genomic regions with the highest level of homozygosity, corresponding to the top 1% of SNPs (H-score value > 0.107 for the CSD and > 0.116 for the NRD group) and found ten regions distributed on seven chromosomes for the CSD group and 15 regions distributed on nine chromosomes for the NRD group. Among these regions, six were partially or totally shared because they were highly homozygous in both groups, while the remaining 13 were specific to only one of the two groups. Matching the positions of these regions with those of the genes annotated in the goat genome version ARS1 and excluding genes for which a symbol or orthologs were not available (i.e. beginning with “LOC”) and transfer RNA gene sequences (TRNA), we identified 133 genes specific to the NRD group, 47 genes specific to the CSD group, and 111 genes common to the two groups.

Then, we identified the regions that showed the largest difference in homozygosity between the two groups (ΔH-score values > 0.06) and found nine regions that were distributed on seven chromosomes and harbored 80 genes of the 291 previously identified genes. These regions were both highly homozygous within a group and capable of differentiating the two NRD and CSD groups.

Finally, when the results of the previous analyses were cross-referenced with the top 1% mean FST values (> 0.09), we identified 44 genes that were shared among all the genes detected by the three methods (Fig. 3). In particular, two groups of genes were specific to the CSD group, i.e. one on chromosome 6 and one on chromosome 5, and one group specific to the NRD group on chromosome 11; finally, a gene cluster on chromosome 12 that was revealed by the FST analysis was highly homozygous in both groups and had a high ΔROH score (Fig. 3 and see Additional file 2: Fig. S2). The complete list of the identified genes is in Additional file 3: Table S1.

Fig. 3
figure 3

Circos plot of the analysis of signatures of selection with a Venn diagram of the results. Circos plot of the signatures of selection in the two groups of Italian goat populations (external blue tracks), of the ΔROH (middle yellow track) and averaged FST (inner red track). The Venn diagram shows the number of regions and genes shared across methods. CSD Central-southern population group, NRD Northern population group

Discussion

Italy is characterized by a wide variety of breeding environments, managements, and traditions; the effect of this variability is particularly evident in goats, which have been traditionally bred in low-input systems, which are strongly affected by climatic conditions. In our study, we used the well-established genomic tool of runs of homozygosity to shed light on the impact of different management practices on Italian goat homozygosity.

Our results show that some populations present an extremely small number of ROH per individual and, consequently, a level of genomic inbreeding near zero. Among these, the two Sicilian breeds Argentata and Messinese are known to have been crossbred, mainly due to the typical management of traditional extensive farms that share common pastures [30]. Another interesting breed is the Capra di Livo-Lariana, with a very low level of inbreeding that can be explained by historical events, including introgression of many unknown individuals from the surrounding valleys.

Another fundamental aspect that emerges from our work is the possibility of monitoring the inbreeding management in the populations through the evaluation of the relationship among the various ROH length classes. Indeed, regardless of the absolute value of FROH, the populations that have a large preponderance of long ROH (> 16 Mb) are more likely to have been under a non-optimal management with frequent mating between closely related individuals, a possible consequence of the reduced number of individuals in the population. One example is the Capra di Teramo breed: the earthquakes that hit the regions where this breed is reared had catastrophic consequences on this already endangered population. On the contrary, Orobica has a more balanced ratio between the different ROH length categories and a smaller number of long ROH. In fact, Orobica is one of the first Italian breeds to have been standardized and reflects a long-term efficient management and the particular attention paid by shepherds to the mating plans.

Moreover, the statistical models performed on ROH number, length, and distribution revealed significant differences between the two NRD and CSD groups. In particular, for the NRD group the number of ROH per individual was larger than for the CSD group, whereas for the CSD group the ROH were longer than for the NRD group. The populations from the NRD group have always been bred in isolated valleys, with natural barriers that prevent the exchange of animals; for this reason, a large number of short ROH, indicating ancient inbreeding events, is expected. On the contrary, the breeds from the CSD group might have undergone, in the past, more admixture events due to the horizontal transhumance practice, the sharing of common pastures, and the presence of multi-breed farms; the more recent standardization is represented by longer ROH.

We also identified signatures of selection that characterize Italian goat populations according to their geographic location. For the NRD group, only four genes on chromosome 11 were found across all the analyses. Among these, the most interesting one is DENND1A, a fertility-related gene [31] that is involved in embryogenesis in cattle. Furthermore, other genes worthy of attention belong to highly homozygous regions but were not found by the ΔROH analysis. In particular, two genes on chromosome 11, HSPA5 and NR5A1, are linked to the production of the anti-Müllerian hormone in grazing cows [32] and are located in a large region on this chromosome that is related to milk production in European, American, and Asian goats [33]. Another region of interest is located on chromosome 13 and hosts genes that are important for pigmentation such as ASIP and RALY [34].

We found a particularly interesting group of genes for the CSD group in a region of chromosome 6 that distinguishes it from the NRD group. This region harbors different genes related to animal growth and development, such as LCORL [35], which has been shown to regulate body size in goats and several other mammals, HERC6 [36], and FAM184B [37]. A part of this homozygous region and another region that we also found on chromosome 5 for the CSD group were previously described by [30], who analyzed only the cluster of Sicilian goats.

Finally, in both NRD and CSD groups, we identified a highly homozygous region on chromosome 12 that was detected by all three analyses and contains 17 genes, which are mostly related to environmental adaptation, for example to hot and arid climates, and body size, including GJB2, GJA3 [38], and PSPC1 [39].

Conclusions

Our findings show that the analysis of ROH is a useful tool not only to identify regions under selection in different breeds, but also to evaluate how their management has evolved over generations. However, this is possible only if a representative recent sample of the specific populations is available, with the potential of expanding the study to historical samples to understand the evolution of a breed’s inbreeding and signatures of selection. ROH assessment can be adopted as a ‘checkpoint’ to assess whether selection in a population is leading to an increase in its average homozygosity and inbreeding, therefore indicating whether a fine-tuning of the breeding scheme is necessary. For these reasons, we recommend the implementation of this tool in the routine evaluation of biodiversity and, consequently, the management of autochthonous populations that are bred in a low-input system as typical of marginal rural areas.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available on Mendeley Data (https://doi.org/10.17632/hnd59x6gmg.1; URL: https://data.mendeley.com/datasets/hnd59x6gmg/1).

References

  1. Bruford MW, Ginja C, Hoffmann I, Joost S, Orozco-terWengel P, Alberto FJ, et al. Prospects and challenges for the conservation of farm animal genomic resources, 2015–2025. Front Genet. 2015;6:314.

    Article  Google Scholar 

  2. Meuwissen THE, Sonesson AK, Gebregiwergis G, Woolliams JA. Management of genetic diversity in the era of genomics. Front Genet. 2020;11:880.

    Article  Google Scholar 

  3. Senczuk G, Mastrangelo S, Ciani E, Battaglini L, Cendron F, Ciampolini R, et al. The genetic heritage of Alpine local cattle breeds using genomic SNP data. Genet Sel Evol. 2020;52:40.

    Article  CAS  Google Scholar 

  4. Talenti A, Dreger DL, Frattini S, Polli M, Marelli S, Harris AC, et al. Studies of modern Italian dog populations reveal multiple patterns for domestic breed evolution. Ecol Evol. 2018;8:2911–25.

    Article  Google Scholar 

  5. Cortellari M, Barbato M, Talenti A, Bionda A, Carta A, Ciampolini R, et al. The climatic and genetic heritage of Italian goat breeds with genomic SNP data. Sci Rep. 2021;11:10986.

    Article  CAS  Google Scholar 

  6. Manfredi MT, Di Cerbo AR, Zanzani S, Stradiotto K. Breeding management in goat farms of Lombardy, northern Italy: risk factors connected to gastrointestinal parasites. Small Rumin Res. 2010;88:113–8.

    Article  Google Scholar 

  7. Sandrucci A, Bava L, Tamburini A, Gislon G, Zucali M. Management practices and milk quality in dairy goat farms in Northern Italy. Ital J Anim Sci. 2019;18:1–12.

    Article  CAS  Google Scholar 

  8. Crepaldi P, Corti M, Cicogna M. Factors affecting milk production and prolificacy of Alpine goats in Lombardy (Italy). Small Rumin Res. 1999;32:83–8.

    Article  Google Scholar 

  9. ISTAT. 6° Censimento agricoltura. ISTAT; 2010. https://www.istat.it/it/archivio/66591/. Accessed 10 Nov 2021.

  10. Di Trana A, Sepe L, Di Gregorio P, Di Napoli MA, Giorgio D, Caputo AR, et al. The role of local sheep and goat breeds and their products as a tool for sustainability and safeguard of the Mediterranean environment. In: Vastola A, editor., et al., The sustainability of agro-food and natural resource systems in the Mediterranean basin. New York: Springer Publishing Company; 2015. p. 77–112.

    Google Scholar 

  11. Paschino P, Stocco G, Dettori ML, Pazzola M, Marongiu ML, Pilo CE, et al. Characterization of milk composition, coagulation properties, and cheese-making ability of goats reared in extensive farms. J Dairy Sci. 2020;103:5830–43.

    Article  CAS  Google Scholar 

  12. Nannini D, Calistri P, Giovannini A, Di Ventura M, Cafiero MA, Ferrari G, et al. Health management of large transhumant animal populations and risk of bluetongue spread to disease-free areas. Vet Ital. 2004;40:707–12.

    CAS  PubMed  Google Scholar 

  13. Eusebi PG, Martinez A, Cortes O. Genomic tools for effective conservation of livestock breed diversity. Diversity. 2019;12:8.

    Article  Google Scholar 

  14. Gorssen W, Meyermans R, Janssens S, Buys N. A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species. Genet Sel Evol. 2021;53:2.

    Article  CAS  Google Scholar 

  15. Kirin M, McQuillan R, Franklin CS, Campbell H, McKeigue PM. Genomic runs of homozygosity record population history and consanguinity. PLoS One. 2010;5:e13996.

    Article  Google Scholar 

  16. Nothnagel M, Lu TT, Kayser M, Krawczak M. Genomic and geographic distribution of SNP-defined runs of homozygosity in Europeans. Hum Mol Genet. 2010;19:2927–35.

    Article  CAS  Google Scholar 

  17. Rodríguez-Ramilo ST, Elsen JM, Legarra A. Inbreeding and effective population size in French dairy sheep: comparison between genomic and pedigree estimates. J Dairy Sci. 2019;102:4227–37.

    Article  Google Scholar 

  18. Meyermans R, Gorssen W, Buys N, Janssens S. How to study runs of homozygosity using plink? a guide for analyzing medium density SNP data in livestock and pet species. BMC Genomics. 2020;21:94.

    Article  CAS  Google Scholar 

  19. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015;4:7.

    Article  Google Scholar 

  20. Wright S. The genetical structure of populations. Ann Eugen. 1951;15:323–54.

    Article  CAS  Google Scholar 

  21. McQuillan R, Leutenegger A-L, Abdel-Rahman R, Franklin CS, Pericic M, Barac-Lauc L, et al. Runs of homozygosity in European populations. Am J Hum Genet. 2008;83:359–72.

    Article  CAS  Google Scholar 

  22. Onzima RB, Upadhyay MR, Doekes HP, Brito LF, Bosse M, Kanis E, et al. Genome-wide characterization of selection signatures and runs of homozygosity in Ugandan goat breeds. Front Genet. 2018;9:318.

    Article  Google Scholar 

  23. Macciotta NPP, Colli L, Cesarani A, Ajmone-Marsan P, Low WY, Tearle R, et al. The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events. Genet Sel Evol. 2021;53:20.

    Article  CAS  Google Scholar 

  24. JMP®, Version 16: SAS Institute Inc. 1989–2021. https://www.jmp.com/support/notes/35/282.html. Accessed 8th Oct 2021.

  25. Copenhaver GP, Browne WE, Preuss D. Assaying genome-wide recombination and centromere functions with Arabidopsis tetrads. Proc Natl Acad Sci USA. 1998;95:247–52.

    Article  CAS  Google Scholar 

  26. Kaback DB. Chromosome-size dependent control of meiotic recombination in humans. Nat Genet. 1996;13:20–1.

    Article  CAS  Google Scholar 

  27. Kaback DB, Barber D, Mahon J, Lamb J, You J. Chromosome size-dependent control of meiotic reciprocal recombination in Saccharomyces cerevisiae: the role of crossover interference. Genetics. 1999;152:1475–86.

    Article  CAS  Google Scholar 

  28. Mastrangelo S, Tolone M, Sardina MT, Sottile G, Sutera AM, Di Gerlando R, et al. Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep. Genet Sel Evol. 2017;49:84.

    Article  Google Scholar 

  29. He S, Di J, Han B, Chen L, Liu M, Li W. Genome-wide scan for runs of homozygosity identifies candidate genes related to economically important traits in Chinese Merino. Animals (Basel). 2020;10:524.

    Article  Google Scholar 

  30. Mastrangelo S, Di Gerlando R, Sardina MT, Sutera AM, Moscarelli A, Tolone M, et al. Genome-wide patterns of homozygosity reveal the conservation status in five Italian goat populations. Animals (Basel). 2021;11:1510.

    Article  Google Scholar 

  31. Jaton C. Genome-wide association study and in silico functional analysis of the number of embryos produced by Holstein donors. J Dairy Sci. 2018;101:7248–57.

    Article  CAS  Google Scholar 

  32. Gobikrushanth M, Purfield DC, Canadas ER, Herlihy MM, Kenneally J, Murray M, et al. Anti-Müllerian hormone in grazing dairy cows: Identification of factors affecting plasma concentration, relationship with phenotypic fertility, and genome-wide associations. J Dairy Sci. 2019;102:11622–35.

    Article  CAS  Google Scholar 

  33. Bertolini F, Cardoso TF, Marras G, Nicolazzi EL, Rothschild MF, Amills M, et al. Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats. Genet Sel Evol. 2018;50:59.

    Article  Google Scholar 

  34. Guo J, Tao H, Li P, Li L, Zhong T, Wang L, et al. Whole-genome sequencing reveals selection signatures associated with important traits in six goat breeds. Sci Rep. 2018;8:10405.

    Article  Google Scholar 

  35. Saif R, Henkel J, Jagannathan V, Drögemüller C, Flury C, Leeb T. The LCORL locus is under selection in large-sized Pakistani goat breeds. Genes (Basel). 2020;11:168.

    Article  CAS  Google Scholar 

  36. Cheng J, Cao X, Hanif Q, Pi L, Hu L, Huang Y, et al. Integrating genome-wide CNVs into QTLs and high confidence GWAS core regions identified positional candidates for sheep economic traits. Front Genet. 2020;11:569.

    Article  CAS  Google Scholar 

  37. An B, Xia J, Chang T, Wang X, Miao J, Xu L, et al. Genome-wide association study identifies loci and candidate genes for internal organ weights in Simmental beef cattle. Physiol Genomics. 2018;50:523–31.

    Article  CAS  Google Scholar 

  38. Kim ES, Elbeltagy AR, Aboul-Naga AM, Rischkowsky B, Sayre B, Mwacharo JM, et al. Multiple genomic signatures of selection in goats and sheep indigenous to a hot arid environment. Heredity (Edinb). 2016;116:255–64.

    Article  CAS  Google Scholar 

  39. Edea Z, Dadi H, Dessie T, Uzzaman MR, Rothschild MF, Kim ES, et al. Genome-wide scan reveals divergent selection among taurine and zebu cattle populations from different regions. Anim Genet. 2018;49:550–63.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors thank all the researchers involved in the Italian Goat Consortium and the Italian National Breeders Association of Sheep and Goat (Asso.Na.Pa). The authors acknowledge support from the University of Milan through the APC initiative.

Author information

Authors and Affiliations

Authors

Contributions

PC and AT conceived the study. MC, AB, and ES performed the analyses. AB, AN, and SF performed lab work. DM, EL, FMS, FP, LC, PAM, RC, SM, LL, and PC provided samples. MC, AB, AT, and PC contributed to data interpretation and drafted the manuscript. All authors contributed to the article with generation of data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Arianna Bionda.

Ethics declarations

Ethics approval and consent to participate

All authors declare that animal samples were obtained in compliance with local/national laws in force at the time of sampling. Genotyping data exchange was in accordance with national and international regulations. The management and handling of the animals involved in the study were performed according to the Italian and European legislation on animal welfare (D.lgs n. 146/2001, Council Directive 98/58/CE) and adhering to the ARRIVE ESSENTIAL 10 guidelines, where applicable. Blood samples were taken by veterinary surgeons according to the recommendations of the European directive 2010/63 and no actual experimental research on animals was performed. The Ethical committee of the Department of Veterinary Science of the University of Messina approved the experimental protocol (code 046/2020).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Geographic distribution (a), phylogeny tree (b), and multidimensional scaling analysis (c) of all the Italian goat breeds included in the study.

Additional file 2: Figure S2.

Manhattan plots representing the signals of signatures of selection in the two population groups (CSD and NRD), of the ΔROH, and averaged FST. CSD = Central-southern population group; NRD = Northern population group.

Additional file 3: Table S1.

List of the genes identified by the analyses of the signatures of selection: top 1% homozygosity score in CSD and NRD groups, ΔROH, and averaged FST. CSD = Central-southern population group; NRD = Northern population group.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cortellari, M., Bionda, A., Negro, A. et al. Runs of homozygosity in the Italian goat breeds: impact of management practices in low-input systems. Genet Sel Evol 53, 92 (2021). https://doi.org/10.1186/s12711-021-00685-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-021-00685-4