Abstract
Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (FROH). The average pedigree inbreeding was 0.07, while the average FROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for Fped and FROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.
Introduction
The Holstein dairy cattle have been bred in a pioneering effort to increase milk yield over the last century while, more recently, the emphasis on functional, longevity, and fertility traits has grown. This breed is counted in over 150 countries, and it is the most common dairy cattle breed worldwide (). Despite its abundance, concerns on the decreasing level of genetic diversity of the Holstein populations have been recently raised in several countries (–). These current alarms stem from the intense directional selection practiced in the last century, which has caused loss of genetic variability. The nineteenth century harbored the adoption of selection theories, such as the selection index theory (), and advanced statistical methods, such as the best linear unbiased prediction (BLUP) (), promoting a remarkable improvement in the genetics of dairy cattle breeds. In addition, the implementation of artificial insemination (AI) over 80 years ago boosted the impact of elite sires worldwide, thus allowing a superior genetic gain per generation (). In the last 20 years, the advances in high-throughput genotyping procedures allowed the development of single-nucleotide polymorphism (SNP) chips at reasonable price, thereby determining the application of genomic selection (GS) based on SNP arrays () in several Holstein breeding programs in different countries (, –). A recent study evaluated the impact of GS implementation on the US Holstein and showed an outstanding increase in the annual genetic gain rate, ranging from 50 to 100% for yield traits and from 200 to 300% for fitness traits (). The increase in the reliability of genomic breeding values (GEBVs) over traditional estimated breeding values (EBVs) for young bulls is quite remarkable, reaching up to a 20% increase for some traits (). Since young bulls can now be selected as sires based on their GEBV at a very early stage, generation intervals (GIs) have been shortened significantly (, , , ). However, there is increasing evidence that the maximization of genetic gain is counterbalanced by a reduction of genetic diversity within the breed. Moreover, it has been shown that the annual inbreeding rate has increased in several dairy cattle populations after the implementation of GS (, , ). This may be partly attributed to the fact that, while the overall number of sires of bulls has increased since the introduction of GS, the number of popular bulls (siring half of the young bulls entering AI) has remained fairly stable ().
Managing genetic diversity determines the long-term sustainability of the livestock production sector. The impact of climate change on livestock, market demand fluctuations, and the increase in human population urgently require a sufficient reservoir of genetic diversity (). It is therefore crucial to evaluate the overall pool of genetic diversity in both commercial and local breeds to preserve biodiversity. Traditionally, the study of genetic diversity relies on pedigree information, where inbreeding can be estimated as the probability of an individual to have two identical alleles by descent (). Analyses at the pedigree level are particularly effective to evaluate the state of genetic diversity in small and underdeveloped populations with limited financial resources that cause the unavailability of more advanced technologies such as genomic data (–). However, those estimates highly depend on quality and depth of the pedigree information and rely on the assumption that no relationship exists among founder animals; hence, pedigree-based inbreeding estimates tend to underestimate actual inbreeding coefficients (, ). The advent of genomics allowed researchers to gain insights into genetic diversity by using genotype data (). One of the most well-established methods to detect within-breed loss of genotypic diversity is the run-of-homozygosity (ROH) detection (). The ROHs are long consecutive homozygous segments distributed across the genome, which arise from identical-by-descendent haplotype (, ). Hence, ROHs have been commonly used to estimate genomic inbreeding (FROH) in several species such as cattle (, –), horses (–), pigs (, ), sheep (, ), and goats (, ). In contrast to pedigree-based inbreeding estimates, FROH can capture the variation due to Mendelian sampling and linkage during gamete formation ().
In Italy, the most reared dairy cattle breed is the Italian Holstein, counting for more than 1,000,000 live animals and about 9,500 breeders, with an average of 10,386 kg of milk produced per lactation/cow in 2020 (). However, few studies have evaluated the level of genetic diversity in this breed by studying SNP data. More precisely, previous reports on genomic inbreeding in the Italian Holstein were based on 50K SNP data in 2,093 bulls () and in a set of 96 animals (). Moreover, a recent study that aimed to evaluate the presence of genomic divergence in Italian Holstein cows bred for different production chains calculated the inbreeding based on ROH in 1,000 Italian Holstein cows (). Nevertheless, to the best of the authors' knowledge, the genetic variability in the Italian Holstein has not been fully explored, especially using a large database on the female side and with deep historical data. Thus, in this study, we aimed to investigate the genetic diversity using data from Italian Holstein cows born between 2002 and 2020. The specific aims of the study were to (i) calculate the inbreeding based on pedigree and genotype data, (ii) evaluate changes in the effective population size throughout generations, and (iii) test the effect of GS on genetic diversity and GI.
Materials and Methods
Records used in this study were obtained from archived data provided by the Italian National Association of Holstein, Brown Swiss, and Jersey Breeders (ANAFIBJ), and as such, no approval was required for animal experimental purposes from the Animal Care Committee unit of the University of Parma. The consent for data use was obtained by ANAFIBJ.
Pedigree Data
Pedigree records for the genotyped animals were provided by ANAFIBJ. Pedigree information consisted of 393,607 individuals born between 1898 and 2020 with 26,226 males and 367,381 females over 24 generations of pedigree depth. To evaluate the role of pedigree depth in the inbreeding estimates, the complete generation equivalent (CGE) was calculated using the optiSel package () in R software ().
Genotype Data
A total of 95,497 genotyped Italian Holstein cows born between 2002 and 2020 were available for this study. Cows were genotyped with a variety of SNP panels, ranging from low- to high-density panels. The animals genotyped with low-density panels were imputed to medium density (85K) using PedImpute (). To guarantee high accuracy during the imputation pipeline, females were retained for this study only when both sire and sire of the dam were (i) genotyped and (ii) used in the imputation pipeline. Quality control (QC) excluding poorly genotyped and faulty data was performed on the 29 autosomal chromosomes by using PLINK v1.90 (). The QC was based on the following criteria: call rate of <95%, parent–offspring SNP mismatch of <0.01, minor allele (<0.01) and genotype (<0.001) frequencies, and extreme deviation from the Hardy–Weinberg equilibrium (P < 0.005).
Pedigree and SNP-Based Inbreeding Coefficients
The pedigree-based inbreeding coefficient (Fped) was defined as the probability of an individual presenting two identical alleles by descent (). The coefficients were computed by using the optiSel () package in R (). The SNP-based inbreeding coefficient was derived by means of ROH assessments, a segment-based approach. The ROH segments were detected by using the detectRUNS package () in R () and defined as follows: (i) at least 15 SNPs in a run, (ii) a minimum length of a run equal to 1 Mb, (iii) a maximum distance between consecutive SNPs in a window of 500 kb, (iv) a lower density limit of one SNP per 100 kb, and (v) a maximum of one missing and one heterozygous SNP being allowed in a run. The genomic inbreeding coefficient (FROH) was calculated as follows ():
with LROH being the sum of the length of ROHs per cow and LAUTO the total length of the autosomal genome covered by SNPs (in this study 2.48 Gbp). The correlation between Fped and FROH was calculated by means of Pearson's product–moment correlation (r). To compare the variability of Fped and FROH, we calculated the coefficient of variation of these two measurements as the ratio between the standard deviation of inbreeding and its overall mean. Moreover, the rFped−FROH was estimated per year. To evaluate the effect of CGE on the relationship between Fped and FROH, the database was divided into cows with CGE ≤ 10 (N = 21,028) and cows with CGE > 10 (N = 53,457). Based on the hypothesis that ROH length reflects the chronological time points when inbreeding happened, the genomic inbreeding was expressed separately for six length ROH categories to differentiate old and recent inbreeding (1 > ROH ≤ 2, 2 > ROH ≤ 4, 4 > ROH ≤ 8, 8 > ROH ≤ 16, 16 > ROH ≤ 32, and ROH > 32 Mbp) by dividing the database into four birth-year classes, considering an average GI of 5 years as found from the pedigree data.
Effective Population Size
The effective population size (Ne) of an actual population can be defined as the size of a hypothetical ideal population resulting in the same amount of genetic drift as is present in the real population (). In this study, we estimated the historical and recent Ne based on both pedigree and SNP data. The optiSel () package in R () was used to estimate the Ne based on pedigree data from 1960 to 2020. The SNeP v.1.1 software was used to estimate the trends of historical Ne based on linkage disequilibrium (LD) on SNP data () by using animals born in the latest year (677 cows in 2020). The Ne was estimated for the latest 30 generations as the first introduction of the Dutch Friesian in Italy dates back to 1,870 followed by North American Holstein Friesians from 1923. Since different methods for Ne estimation based on LD are available in the literature, two analyses were performed by (a) using default settings, except for the recombination rate, which was inferred following the Sved and Feldman method (), and the mutation rate in cattle [α = 2.2 ()]; and (b) including a restriction on the maximum distance used to calculate LD. The latter parameterization, together with the Sved and Feldman mutation rate modifier, allowed us to estimate Ne for the most recent generations (, ).
The Role of GS on Genetic Diversity and GI
The rate of inbreeding (ΔF) per year was calculated as the inverse of the slope of the regression of on the year of birth, where x was equal to the average of the parameter each year (Fped and FROH) (). The annual rate was multiplied by the GI (in the Italian Holstein being equal on average to 5 years as shown in this study) to obtain the rate per generation (ΔFgen). To assess the effect of different selection strategies on genetic diversity in the Italian Holstein breed, i.e., classical progeny testing (PTS) vs. GS, we divided the database into two 5 year birth cohorts by considering all animals born before (2006 to 2010) and after the introduction of GS (2015 to 2019). Then, we tested the equality of PTS and GS means. Moreover, to evaluate the impact of the introduction of GS on genetic diversity, we used the following linear model, using the R function lm ():
where Yi is the variable of interest for each cow i (Fped and FROH), xi is the birth year of each cow i, and βPTS is the associated coefficient of regression if cow i was born in the PTS cohort or βPTS + δ if cow i was born in the GS cohort. The impact of GS on the inbreeding rate was measured with the δ coefficient. Analysis of variance was used to test the significance of the δ value. The relative change (RC) of the slopes of regression before and after GS was computed as RC = δ/β1, with β1 being the slope in the second evaluated period. The assessment of the RC allowed the evaluation of slope value through time (). The effect of the introduction of GS was also evaluated via the length of the GI. The GI, defined as the average age of parents when their offspring were born, was calculated for all the four selection pathways (sire of bulls, dam of bulls, sire of cows, and dam of cows) using the 393,607 individuals present in the pedigree. The GI was estimated using the Pedig software ().
Results
Pedigree and SNP-Based Inbreeding Coefficients
A total of 84,443 SNPs and 74,485 cows were kept after QC. The average genotyping call rate was 0.99, and the average pedigree depth based on CGE was equal to 10.67 (SD = 1.12). The cows descended from 3,058 sires and 59,377 dams. In total, over 50% of the cows was born between 2016 and 2020 (last 5 years), whereas 34,286 cows were born between 2002 and 2015. The mean Fped was equal to 0.07 (SD = 0.02) ranging between 0.01 and 0.32 (CV = 0.29). The FROH mean value was more than doubled (0.17; SD = 0.03), with a minimum and maximum of 0.05 and 0.50, respectively (CV = 0.20). A total of 839 cows showed an Fped exceeding the mean + 3 SD (corresponding to Fped ≥ 0.13), which were defined as highly inbred females. In the case of the FROH, 507 cows presented an inbreeding higher than the mean + 3 SD (FROH ≥ 0.27) (Figure 1). Approximately 23% of the highly inbred cows were in common from the comparison between the two methodologies. The Pearson correlation between Fped and FROH was equal to 0.68 (confidence interval: 0.676–0.683), P-value < 2.2e−16) as shown in Supplementary Figure 1. The average CGE in the entire database was equal to 10.06, showing an increase throughout the studied period from 2002 (average CGE = 7.5) to 2019 (average CGE = 11.9) (Figure 2A). The correlation ranged from 0.44 in 2005 (N = 136) to 0.89 in 2003 (N = 16); however, those estimates might be an artifact produced by the limited number of cows for those years; thus, they should be considered with caution. Since 2010, the number of born and genotyped cows per year has been above 1,000, and the rFped−FROH was steadier, ranging between 0.58 and 0.69. The rFped−FROH was found higher in cows with CGE > 10 compared to those with CGE ≤ 10 (0.68 and 0.59, respectively; Figure 2B).
Figure 1
Figure 2
The FROH estimates and average number of ROH per cow based upon the six ROH length classes (to differentiate old and recent inbreeding) of each of the four birth-year cohorts are presented in Table 1. In the case of ROH length classes of 1–2, 2–4, and 4–8 Mbp, all the cows, regardless of birth-year cohort, exhibited some degree of inbreeding with a mean value ranging from 0.02 in the 2–4 and 4–8 Mbp in the two oldest birth-year cohorts (2002–2005 and 2006–2010) to 0.04 in the shortest length class for all the birth classes. From the 8–16 Mbp class and above, not all the cows exhibited ROH of such lengths. Only 39% of the cows exhibited ROH in the longest length class (>32 Mbp) with most of them belonging to the latest two birth-year cohorts (2011–2015 and 2016–2020) (data not shown).
Table 1
| 2002–2005 (n. 248)a | 2006–2010 (n. 3,883)a | 2011–2015 (n. 30,156)a | 2016–2020 (n. 40,198)a | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Length Class (Mbp) | Mean | SD | N.b | Mean | SD | N.b | Mean | SD | N.b | Mean | SD | N.b |
| Inbreeding coefficients based on ROH (FROH) | ||||||||||||
| 1> ROH ≤ 2 | 0.04 | 0.01 | 69.4 | 0.04 | 0.01 | 68.7 | 0.04 | 0.01 | 69.8 | 0.04 | 0.01 | 70.4 |
| 2> ROH ≤ 4 | 0.02 | 0.01 | 20.1 | 0.02 | 0.01 | 20.9 | 0.03 | 0.01 | 22.8 | 0.03 | 0.01 | 24.0 |
| 4> ROH ≤ 8 | 0.03 | 0.01 | 11.4 | 0.03 | 0.01 | 12.6 | 0.03 | 0.01 | 14.1 | 0.04 | 0.01 | 15.6 |
| 8> ROH ≤ 16 | 0.03 | 0.01 | 6.10 | 0.04 | 0.01 | 6.98 | 0.04 | 0.01 | 7.88 | 0.04 | 0.01 | 9.07 |
| 16> ROH ≤ 32 | 0.02 | 0.01 | 2.12 | 0.03 | 0.02 | 2.67 | 0.03 | 0.02 | 2.96 | 0.03 | 0.02 | 3.38 |
| >32 | 0.02 | 0.01 | 0.37 | 0.02 | 0.01 | 0.47 | 0.02 | 0.01 | 0.52 | 0.02 | 0.01 | 0.60 |
| Total | 0.15 | 0.03 | 109.5 | 0.17 | 0.03 | 112.5 | 0.18 | 0.03 | 118.1 | 0.19 | 0.03 | 123.2 |
Descriptive statistics of inbreeding based on runs of homozygosity (ROH) divided by six length classes per each of the four birth year cohorts.
n., number of animals per each birth year cohort;
N., average number of ROH per animal.
Effective Population Size Based on Genotype and Pedigree Data
The two settings for the Ne estimation based on LD were comparable with an Ne reduction from generation 30 till generation 13 (Figure 3A). The Ne was equal to 140 animals at generation 30 and to 96 at generation 13 based on both settings. The second parameterization (Figure 3B) allowed the evaluation of the Ne trend at more recent generations and showed a sharp increase (almost doubled) within the last five generations, reaching to 120 in the most recent years.
Figure 3
The patterns of the Ne based on pedigree and SNP data were similar to each other, with differences in scale (Supplementary Figure 2). The Ne based on pedigree data was calculated from 1960 to 1964 (generation 12) to 2016–2020 (generation 1) (Supplementary Figure 2). This was equal to 87 and 55 animals in the oldest and earliest generations, respectively. A decrease in the Ne was found from generation 12 to generation 6 (1991–1995), with an Ne of 43 in generation 6. In contrast, from generation 5 (1996–2000), a general increase was observed. Likewise, for the Ne based on SNP data, a decrease was observed since generation 9 (1976–1980), followed by a plateau, with Ne equal to 89 animals from generation 8 to generation 6 (1981–1995). An increase in the Ne was observed with a value of 120 animals in the latest generation.
The Role of GS on Genetic Diversity and GI
GIs based on pedigree information were calculated from 1960 to 2018 for all four pathways of selection by using all animals in the pedigree file (Figure 4). For sire of bulls and cows, an increase in GI was observed till 1984, with a GI equal to 11.08 years for bulls and 8.74 years for cows. In contrast, a tendency to decrease was found from 1985 onwards for both pathways. A noticeable drop occurred in both the sire of bulls and sire of cows from 2011 to 2018, with the lowest GI in 2017 for the former (2.34 years) and in 2018 for the latter (3.6 years). The dam-of-bulls pathway decreased from 1962 (8.64 years) to 1992 (3.8 years) with some oscillations. The GI in this pathway remained approximately steady from 1993 to 2011, with a minimum of 3.68 years in 2010 and a maximum of 4.24 years in 1996, whereas it declined from 2012 to 2018 (2.43 years). In the case of the dam-of-cows pathway, a decrease was visible from 1960 (8.22 years) to 1992 (4.03 years) with few fluctuations. From 1993 onwards, the GI remained stable with an average value of 3.68 years, which varied between 3.01 and 4.0 years.
Figure 4
The overall inbreeding rate per year was equal to +0.27% and +0.44% for Fped and FROH throughout the studied period, which corresponds to about +1.35% and +2.2% ΔFgen, respectively. The mean difference in Fped and FROH based on a two-sample t-test between the two periods (PTS and GS) was significant (P < 0.001) for both inbreeding estimates (Figure 5). The average Fped was equal to 0.05 in the PTS and 0.07 in the GS. The average FROH was equal to 0.14 and 0.17 in the PTS and GS, respectively. The overall inbreeding rate per year in the PTS was equal to 0.14% and 0.32% based on Fped and FROH, whereas it increased up to 0.47% (based on Fped) and 0.70% (based on FROH) in the GS. The RC in inbreeding comparing the two periods (PTS and GS) was equal to 2.36 and 1.19 from the Fped and FROH (Table 2 and Figure 5). The overall GI (calculated as the average among the four pathways) decreased in the GS by a factor of 1.8 compared to the PTS period (Figure 5). Since this latter reduction, the ΔFgen also changed between GS and PTS, being equal to +0.75% (based on Fped) and +1.72% (based on FROH) in the PTS and 1.41 and 2.1% in the GS period, based on Fped and FROH, respectively.
Figure 5
Table 2
| Parameter | Fped | FROH |
|---|---|---|
| bPTS (±SE)a | 0.14% (2.0 × 10−04) | 0.32% (4.5 × 10 −04) |
| bGS (±SE)b | 0.47% (7.7 × 10−05) | 0.70% (1.4 × 10 −04) |
| δc | 0.0033 | 0.0031 |
| p-value of δ- | <0.0001 | <0.0001 |
| RCd | 2.36 | 1.19 |
Parameters used to estimate the differences in the inbreeding trend based on pedigree and genotype data between the PTS and GS periods.
bPTS is the slope in percentage of regression per Fped and FROH for cows born between 2006 and 2010 (PTS selection);
bGS - is the slope in percentage of regression per Fped and FROH for cows born between 2015 and 2019 (GS);
δ is the difference between the slopes of regression of each inbreeding measurement depending on the two 5-year birth class;
RC is the relative change equals to with β1 the slope of the second evaluated period.
Discussion
Pedigree and Genotype-Based Inbreeding Coefficients
In this study, we evaluated the genetic diversity in the Italian Holstein breed by using pedigree and SNP data from cows genotyped throughout a period of 19 years. The primary goal was to estimate variation in inbreeding, effective population size, and GI and to compare those aspects prior to and after the introduction of GS in the breed. In line with the development and spread of the genotyping technology in the dairy sector (), from 2010, the number of genotyped Italian Holstein cows per year rapidly increased, and it reached over 15,000 cows in 2016. This amount of information provided new opportunities to explore genetic diversity at the genomic level. The most traditional source of information to evaluate genetic diversity is pedigree information (, ). However, the accuracy of genetic diversity estimated using pedigree data highly relies on the quality and depth of the recorded data (). In this study, the pedigree depth showed high values with an average CGE of 10.67, which agrees with previous studies on the Dutch Holstein breed (CGE between ~10 and ~14) (, ). Nevertheless, the CGE influenced the correlation between Fped and FROH. Indeed, the correlation dropped to 0.59 when considering cows with a CGE lower than 10, and it reached a value of 0.68 for animals with CGE above 10. The reduction of the correlation values for animals with lower CGE was expected, as the genealogical approach strongly depends on pedigree completeness and quality. Nevertheless, the obtained correlations are in line with the values found in previous studies, regardless of the CGE assessment. In the Spanish Holstein population, a correlation of 0.57 was shown (), and in four different cattle breeds (Brown Swiss, Fleckvieh, Norwegian Red, and Tyrol Gray), correlations were reported to be between 0.50 and 0.72 (). A unity correlation is not expected between Fped and FROH, since the former does not account for the Mendelian sampling variation, as the latter does (, ).
Effective Population Size Based on Genotype and Pedigree Data
Prior to herdbook formation, the Ne was large for most cattle populations, i.e., in the order of tens of thousands (, ). However, the advent of closed policy in the breeding management from the nineteenth century onwards drastically reduced the gene flow among populations, causing loss of genetic diversity. The major changes in breeding programs in the last century have exacerbated this reduction. It has been shown that the current Ne in most modern and commercial dairy cattle breeds is close to 100 animals (). Estimates of Ne found in this study based on pedigree and SNP data are in the range of those formerly published in other Holstein populations. Pedigree-based Ne estimates ranged from 39 in the US Holstein in 2000 () to 114 in the Canadian Holstein population between 2000 and 2007 (). In this study, the Ne based on pedigree in 2000 was equal to 47 animals, thus slightly higher compared to the US Holstein population, yet comparable. In contrast, an over two-time lower Ne between 2000 and 2007 (Ne = 49) was found in our study compared to that reported in the Canadian Holstein population. Since 2000, there was a reduction in the inbreeding rate per generation in the Canadian Holstein () due to the introduction of the average relationship values in the breeding program. Since in the Italian Holstein breed this type of control was not implemented at that time, we suspect that this might be the reason behind those differences among the two Holstein populations in terms of Ne. SNP-based Ne estimates ranged from 69 in the Dutch Holstein () to roughly 135 in the New Zealand Holstein cattle (). In contrast to pedigree-based inbreeding estimates, where the whole genealogical data are generally used, the SNP-based estimates are highly dependent on the genotyped animals and SNP data used for the estimation. For this reason, the comparison with other studies might be less straightforward. In addition, several methods for the Ne estimations are available based on SNP data (as an example, based on homozygous segments, marker-by-marker homozygosity, and similarity and based on LD), leading to different results. The Ne values found in the Dutch Holstein population were calculated using several SNP-based estimates (ranging between 69 and 100) (), and most of them were lower than those found for the Italian Holstein. A possible explanation of this difference might be due to the use of different sexes in the two studies: in the Dutch Holstein (), proven sires were genotyped, whereas in this work only females were analyzed. Thus, due to differences in the intensity of selection, we expect to have higher Ne based on female data, when compared to proven sires. New Zealand and Australian Holstein cattle () showed analogous estimates to the Italian Holstein in the latest generations in our study. The present study and the one by de Roos et al. () were based on LD-based Ne estimates, making the comparison more forthright. Nevertheless, a similar Ne across countries is expected for highly specialized worldwide breeds, due to the exchange of genetic material and similar applied breeding strategies. The Holstein fully represents this scenario as it has been bred in over 150 countries, and it currently dominates commercial dairy production worldwide (). Although the overall breeding goal is not always uniform among countries, high selection intensities and a small number of sires are commonly used worldwide (). In the Italian dairy industry, over 80% of bovine milk is produced by this breed, which has been used over time for a multitude of purposes, from drinkable milk to highly specialized consortia for PDO cheese production (). This latter aspect might have caused a peculiar genomic architecture in the Italian Holstein which should be investigated further in comparison with other international Holstein populations. Nevertheless, in the Italian Holstein population, roughly 49% of the sires is local (Italian nationality), while the remaining ones are from other countries (i.e., 23% from the United States and 7% from Canada), highlighting that what we found in the Italian Holstein might resemble the status quo of other countries as well. Unfortunately, this implies that, despite the census of millions of animals, the Holstein breed is exposed to the same genetic drift and accumulation of inbreeding as a population of roughly 100 individuals. The use of a small number of sires propagated worldwide through AI and other advanced reproductive technologies, as well as mating of close relatives, has boosted this issue (). Analysis on genetic distance and variability of different Holstein populations worldwide may provide additional information on their current relatedness, which can further serve to quantify the genetic variability in this breed at a global scale.
For the latest generation considered in this work, the Ne was equal to 55 and to 120 animals from pedigree and SNP-based estimates, respectively. The Food and Agriculture Organization of the United Nations (FAO) set the Ne critical value at 50 animals, from which the population is expected to lose fitness and variability in the long term (). Since the Italian Holstein population seems to be very close to the critical value from pedigree analysis, appropriate strategies to reduce the loss of genetic diversity are needed to preserve genetic diversity and sustain fitness of the breed in the long term. The reduction in the Ne from 1960 to 1995 found in this study agrees with previous studies in the Canadian, Dutch, French, and US Holstein populations (, ). Interestingly, a similar Ne increase was observed from 2000 in the Canadian Holstein (). This increase was even more evident in our study when looking at SNP Ne estimates, which showed an increase of about 20 animals in the latest three generations. However, the latter result seems counterintuitive if evaluated with the general increase of inbreeding over time found in the Italian Holstein. A potential reason might be the more diverse genetic pool in terms of selected candidates thanks to GS (), leading to a reduction in LD over greater recombinant distances. With GS, we can quantify Mendelian sampling. This serves to differentiate, e.g., between full sibs, hence potentially reducing the co-selection of sibs. Thus, we expect a positive effect on Ne, which was indeed shown in this study based on both pedigree and SNP estimates. Nevertheless, we suspect that the more visible increase in Ne based on SNP compared to pedigree is due to the distinct methodologies employed based on the source of data available. In the case of Ne based on pedigree, the estimate depends on the quality and completeness of the genealogical data, whereas, in the case of SNP-based Ne, the estimate relies on the relationship between the LD variance and effective population size. Therefore, we believe that the latter estimate might be more sensible for the detection of Ne variation compared to pedigree data.
Impact of Genomics and Future Breeding Strategies
The average GI in the Italian Holstein has decreased since 1985, with two sharp declines from 1985 to 1990 and from 2009 onwards. The initial reduction in GI might be attributed to the implementation and use of the BLUP evaluation (). Nevertheless, from 1992 to GS, an almost constant GI was found which might be due to the general tendency of using proven bull sires. The remarkable reduction in the GI since 2009 might be attributed to GS, as predicted in a previous study (). The severe drop in the GI found in this study is in line with reported literature. After 2009, in the Canadian Holstein, a drop of 38% in the GI was registered, as well as in the US Holstein (37%) and in the Dutch Holstein (~35%) (, , ). The ΔFgen also fluctuated throughout the studied period, and it increased from 0.14 to 0.47% (pedigree based) and from 0.32 to 0.70% based on ROH comparing the PTS–GS periods. Similar rates of inbreeding were found in a previous study where over 4,000 Holstein bulls were analyzed, highlighting the effect of GS both on the annual genetic gain increase and on the inbreeding rate (). The observed ΔF in the Italian Holstein is higher than the 1% per generation suggested by FAO guidelines as the critical value for the maintenance of genetic variability in the long term ().
Currently, a web interface is available for breeders to evaluate and manage inbreeding within their herds, where advice on mating strategies is provided by specialists from the breeding association (). However, the routine use of this tool by breeders is not fully known. Thus, further activities to raise awareness on inbreeding control among breeders are advised to prevent further loss of genetic variability. This latter aspect is extremely important if evaluated together with the current and coming expansion of international demand for dairy products due to emerging economies, the need for high quality milk proteins in developing countries, and world population expansion (). In this study, we revealed loss of genetic diversity in the Italian Holstein, which might in turn cause inbreeding depression. So far, in the dairy sector, the management of genetic diversity has been kept within individual herds where ultimately breeding decisions are made. Unfortunately, this is a drawback for the implementation of one of the most effective methods to manage genetic variability, which is the optimal contribution selection (OCS) (). Nevertheless, thanks to the genomic era, the dairy genetic industry is slowly changing, moving toward scenarios where tighter control of the population as a whole is becoming possible (i.e., vertically integrated industries) (77). Thus, the implementation of genomic OCS might prevent any further loss of genetic variability in this breed. Finally, we suggest using inbreeding estimates based on genotype data, as we have shown that they can provide a more precise evaluation of the available variability.
Conclusion
The presented study in the Italian Holstein cattle showed the urgent matter of controlling the loss of genetic diversity in a highly specialized breed and that this loss is not enclosed to small and local breeds only. The implementation of OCS may help to alleviate this problem and prevent any further loss of genetic variability. The preservation of genetic resources is key to ensure the long-term sustainability of this breed which represents one of the most important players in the Italian dairy chain as well as to guarantee future market demands.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: Data supporting this paper were obtained from ANAFIBJ. The genotype data are available only upon agreement with ANAFIBJ. Requests to access these datasets should be directed to jtkaam@anafi.it.
Ethics statement
Ethical review and approval was not required for the animal study because records used in this study were obtained from archived data from the Italian National Association of Holstein, Brown Swiss and Jersey Breeders (ANAFIBJ), and as such, no approval was required for animal experimental purposes from the Animal Care Committee unit of the University of Parma.
Author contributions
ASu, CC-G, and MA conceived the idea and formulated the objectives of this study. CC-G, J-TK, and GS helped in data preparation. MA conducted the analysis and wrote the first draft of the paper. ASu and ASa supervised the project. CD and GS contributed on the data visualization. ASu, CC-G, ASa, GS, RF, and CD critically reviewed the text. All authors read and approved the final manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.773985/full#supplementary-material
Supplementary Figure 1Pearson correlation (above diagonal), scatterplot (below diagonal) and density (diagonal) of inbreeding coefficients measured by ROH (FROH) and pedigree data (Fped) in the Italian Holstein dairy cows.
Supplementary Figure 2Effective population size based on pedigree and SNP data (using the second method applied in the study for the Ne calculation) from 1960 till 2020.
References
1.
BuchananDS. Breeds of dairy cattle (major Bos taurus breeds). In: FuguayJW, editor. Reference Module in Food Science Encyclopedia of Dairy Sciences, 2nd Edn. (Oklahoma State University, OK: Elsevier). p. 1–9.
2.
DoekesHPVeerkampRFBijmaPHiemstraSJWindigJJ. Trends in genome-wide and region-specific genetic diversity in the Dutch-Flemish Holstein–Friesian breeding program from 1986 to 2015. Genet Sel Evol. (2018) 50:15. 10.1186/s12711-018-0385-y
3.
ForutanMAnsari MahyariSBaesCMelzerNSchenkelFSSargolzaeiM. Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle. BMC Genomics. (2018) 19:98. 10.1186/s12864-018-4453-z
4.
MakanjuolaBOMigliorFAbdallaEAMalteccaCSchenkelFSBaesCF. Effect of genomic selection on rate of inbreeding and coancestry and effective population size of Holstein and Jersey cattle populations. J Dairy Sci. (2020) 103:5183–99. 10.3168/jds.2019-18013
5.
HazelLNLushJL. The efficiency of three methods of selection. J Hered. (1942) 33:393–9. 10.1093/oxfordjournals.jhered.a105102
6.
HendersonCR. Best linear unbiased estimation and prediction under a selection model. Biometrics. (1975) 31:423. 10.2307/2529430
7.
BrotherstoneSGoddardM. Artificial selection and maintenance of genetic variance in the global dairy cow population. Philos Trans R Soc B Biol Sci. (2005) 360:1479–88. 10.1098/rstb.2005.1668
8.
MeuwissenTHEHayesBJGoddardME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. (2001) 157:1819–29. 10.1093/genetics/157.4.1819
9.
WiggansGRVanRadenPMCooperTA. The genomic evaluation system in the United States: past, present, future. J Dairy Sci. (2011) 94:3202–11. 10.3168/jds.2010-3866
10.
LundMSde RoosAPde VriesAGDruetTDucrocqVFritzSet al. A common reference population from four European Holstein populations increases reliability of genomic predictions. Genet Sel Evol. (2011) 43:43. 10.1186/1297-9686-43-43
11.
HayesBJBowmanPJChamberlainAJGoddardME. Invited review: genomic selection in dairy cattle: progress and challenges. J Dairy Sci. (2009) 92:433–43. 10.3168/jds.2008-1646
12.
García-RuizAColeJBVanRadenPMWiggansGRRuiz-LópezFJVan TassellCP. Changes in genetic selection differentials and generation intervals in US Holstein dairy cattle as a result of genomic selection. Proc Natl Acad Sci. (2016) 113:E3995–4004. 10.1073/pnas.1519061113
13.
DoubletACCroiseauPFritzSMichenetAHozéCDanchin-BurgeCet al. The impact of genomic selection on genetic diversity and genetic gain in three French dairy cattle breeds. Genet Sel Evol. (2019) 51:1–13. 10.1186/s12711-019-0495-1
14.
Rodríguez-RamiloSTFernándezJToroMAHernándezDVillanuevaB. Genome-Wide estimates of coancestry, inbreeding and effective population size in the spanish holstein population. PLoS ONE. (2015) 10:e0124157. 10.1371/journal.pone.0124157
15.
MigliorFBeaversL. Genetic Diversity Inbreeding: Before After Genomics. Progressive Dairyman (2014). Available online at: https://www.progressivedairy.com/topics/a-i-breeding/genetic-diversity-and-inbreeding-before-and-after-genomics (accessed July 28, 2021).
16.
Rojas-DowningMMNejadhashemiAPHarriganTWoznickiSA. Climate change and livestock: impacts, adaptation, and mitigation. Clim Risk Manag. (2017) 16:145–63. 10.1016/j.crm.2017.02.001
17.
WrightS. Coefficients of Inbreeding and Relationship. Am Nat. (1922) 56:330–8. 10.1086/279872
18.
MarianiESummerAAblondiMSabbioniA. Genetic variability and management in nero di parma swine breed to preserve local diversity. Animals. (2020) 10:538. 10.3390/ani10030538
19.
AblondiMVasiniMBerettiVSuperchiPSabbioniA. Exploring genetic diversity in an Italian horse native breed to develop strategies for preservation and management. J Anim Breed Genet. (2018) 135:450–9. 10.1111/jbg.12357
20.
FabbriMCde RezendeMPGDadousisCBiffaniSNegriniRCarneiroPLSet al. Population structure and genetic diversity of Italian beef breeds as a tool for planning conservation and selection strategies. Animals. (2019) 9:880. 10.3390/ani9110880
21.
CavaniLde SilvaRMOCarreñoLODOnoRKBertipagliaTSFarahMMet al. Genetic diversity of Brazilian Brahman cattle by pedigree analysis. Pesqui Agropecuária Bras. (2018) 53:74–9. 10.1590/s0100-204x2018000100008
22.
SchurinkAArtsDJGDucroBJ. Genetic diversity in the Dutch harness horse population using pedigree analysis. Livest Sci. (2012) 143:270–7. 10.1016/j.livsci.2011.10.005
23.
HowardJTPryceJEBaesCMalteccaC. Invited review: inbreeding in the genomics era: inbreeding, inbreeding depression, and management of genomic variability. J Dairy Sci. (2017) 100:6009–24. 10.3168/jds.2017-12787
24.
PeripolliEMunariDPSilvaMVGBLimaALFIrgangRBaldiF. Runs of homozygosity: current knowledge and applications in livestock. Anim Genet. (2017) 48:255–71. 10.1111/age.12526
25.
CeballosFCJoshiPKClarkDWRamsayMWilsonJF. Runs of homozygosity: windows into population history and trait architecture. Nat Rev Genet. (2018) 19:220–34. 10.1038/nrg.2017.109
26.
CurikIFerenčakovićMSölknerJ. Inbreeding and runs of homozygosity: a possible solution to an old problem. Livest Sci. (2014) 166:26–34. 10.1016/j.livsci.2014.05.034
27.
MakanjuolaBOMalteccaCMigliorFSchenkelFSBaesCF. Effect of recent and ancient inbreeding on production and fertility traits in Canadian Holsteins. BMC Genomics. (2020) 21:1–15. 10.1186/s12864-020-07031-w
28.
KimE-SColeJBHusonHWiggansGRVan TassellCPCrookerBAet al. Effect of artificial selection on runs of homozygosity in U.S. Holstein cattle. PLoS ONE. (2013) 8:e80813. 10.1371/journal.pone.0080813
29.
PurfieldDCBerryDPMcParlandSBradleyDG. Runs of homozygosity and population history in cattle. BMC Genet. (2012) 13:70. 10.1186/1471-2156-13-70
30.
Grilz-SegerGMesaričMCotmanMNeuditschkoMDrumlTBremG. Runs of homozygosity and population history of three horse breeds with small population size. J Equine Vet Sci. (2018) 71:27–34. 10.1016/j.jevs.2018.09.004
31.
SadeghiRMoradi-ShahrbabakMAshtianiSRMSchlampFCosgroveEJAntczakDF. Genetic diversity of Persian Arabian horses and their relationship to other native Iranian horse breeds. J Hered. (2019) 110:173–82. 10.1093/jhered/esy061
32.
MancinEAblondiMMantovaniRPigozziGSabbioniASartoriC. Genetic variability in the Italian heavy draught horse from pedigree data and genomic information. Animals. (2020) 10:1310. 10.3390/ani10081310
33.
AblondiMDadousisCVasiniMErikssonSMikkoSSabbioniA. Genetic diversity and signatures of selection in a native Italian horse breed based on SNP data. Animals. (2020) 10:1–15. 10.3390/ani10061005
34.
Gomez-RayaLRodríguezCBarragánCSilióL. Genomic inbreeding coefficients based on the distribution of the length of runs of homozygosity in a closed line of Iberian pigs. Genet Sel Evol. (2015) 47:81. 10.1186/s12711-015-0153-1
35.
SchiavoGBovoSBertoliniFTinarelliSDall'OlioSNanni CostaLet al. Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds. Animal. (2020) 14:910–920. 10.1017/S175173111900332X
36.
AddoSKlingelSThallerGHinrichsD. Genetic diversity and the application of runs of homozygosity-based methods for inbreeding estimation in German White-headed Mutton sheep. PLoS ONE. (2021) 16:e0250608. 10.1371/journal.pone.0250608
37.
Rodríguez-RamiloSTElsenJMLegarraA. Inbreeding and effective population size in French dairy sheep: comparison between genomic and pedigree estimates. J Dairy Sci. (2019) 102:4227–37. 10.3168/jds.2018-15405
38.
DadousisCCecchiFAblondiMFabbriMCStellaABozziR. Keep Garfagnina alive. An integrated study on patterns of homozygosity, genomic inbreeding, admixture and breed traceability of the Italian Garfagnina goat breed. PLoS ONE. (2021) 16:e0232436. 10.1371/journal.pone.0232436
39.
BurrenANeuditschkoMSigner-HaslerHFrischknechtMReberIMenziFet al. Genetic diversity analyses reveal first insights into breed-specific selection signatures within Swiss goat breeds. Anim Genet. (2016) 47:727–39. 10.1111/age.12476
40.
Allevatori AIA,. Italy: Milk Recording Activity - Official Statistics - Year 2020. Roma (2020). Available online at: http://bollettino.aia.it/Contenuti.aspx?CD_GruppoStampe=RS&CD_Specie=C4 (accessed June 15, 2021).
41.
MarrasGGaspaGSorboliniSDimauroCAjmone-MarsanPValentiniAet al. Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy. Anim Genet. (2015) 46:110–21. 10.1111/age.12259
42.
MastrangeloSSardinaMTToloneMDi GerlandoRSuteraAMFontanesiLet al. Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds. Animal. (2018) 12:2480–8. 10.1017/S1751731118000629
43.
AblondiMMalacarneMCipolat-GotetCvan KaamJ-TSabbioniASummerA. Genome-wide scan reveals genetic divergence in Italian Holstein cows bred within PDO cheese production chains. Sci Rep. (2021) 11:12601. 10.1038/s41598-021-92168-1
44.
WellmannR. Optimum contribution selection for animal breeding and conservation: the R package optiSel. BMC Bioinformatics. (2019) 20:25. 10.1186/s12859-018-2450-5
45.
R Development Core Team. R: A Language and Environment for Statistical Computing (2011).
46.
NicolazziELBiffaniSJansenG. Short communication: Imputing genotypes using PedImpute fast algorithm combining pedigree and population information. J Dairy Sci. (2013) 96:2649–53. 10.3168/jds.2012-6062
47.
PurcellSNealeBTodd-BrownKThomasLFerreiraMARRBenderDet al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:559–75. 10.1086/519795
48.
WrightS. Evolution and the Genetics of Populations. Vol. 4: Variability Within and Among Natural Populations. Chicago, IL: University of Chicago Press (1978).
49.
BiscariniFCozziPGaspaGMarrasG. detectRUNS: an R Package to Detect Runs of Homozygosity Heterozygosity in Diploid Genomes. (2019). Available online at: https://cran.r-project.org/web/packages/detectRUNS/vignettes/detectRUNS.vignette.html#references (accessed January 1, 2021).
50.
McQuillanRLeuteneggerA-LAbdel-RahmanRFranklinCSPericicMBarac-LaucLet al. Runs of homozygosity in European populations. Am J Hum Genet. (2008) 83:359–72. 10.1016/j.ajhg.2008.08.007
51.
BarbatoMOrozco-terWengelPTapioMBrufordMW. SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front Genet. (2015) 6:109. 10.3389/fgene.2015.00109
52.
SvedJAFeldmanMW. Correlation and probability methods for one and two loci. Theor Popul Biol. (1973) 4:129–32. 10.1016/0040-5809(73)90008-7
53.
CorbinLJLiuAYHBishopSCWoolliamsJA. Estimation of historical effective population size using linkage disequilibria with marker data. J Anim Breed Genet. (2012) 129:257–70. 10.1111/j.1439-0388.2012.01003.x
54.
BarbatoMHailerFUpadhyayMDel CorvoMColliLNegriniRet al. Adaptive introgression from indicine cattle into white cattle breeds from Central Italy. Sci Rep. (2020) 10:1279. 10.1038/s41598-020-57880-4
55.
BoichardD. PEDIG : a fortran package for pedigree analysis suited for large populations. In: 7th World Congress on Genetics Applied to Livestock Production, August 19-23.Montpellier (2002). p. 13–28.
56.
WiggansGRColeJBHubbardSMSonstegardTS. Genomic selection in dairy cattle: the USDA experience. Annu Rev Anim Biosci. (2017) 5:309–27. 10.1146/annurev-animal-021815-111422
57.
MeuwissenTLuoZ. Computing inbreeding coefficients in large populations. Genet Sel Evol. (1992) 24:305. 10.1186/1297-9686-24-4-305
58.
OliehoekPABijmaP. Effects of pedigree errors on the efficiency of conservation decisions. Genet Sel Evol. (2009) 41:9. 10.1186/1297-9686-41-9
59.
DoekesHPVeerkampRFBijmaPde JongGHiemstraSJWindigJJ. Inbreeding depression due to recent and ancient inbreeding in Dutch Holstein–Friesian dairy cattle. Genet Sel Evol. (2019) 51:54. 10.1186/s12711-019-0497-z
60.
DoekesHPVeerkampRFBijmaPHiemstraSJWindigJ. Value of the Dutch Holstein Friesian germplasm collection to increase genetic variability and improve genetic merit. J Dairy Sci. (2018) 101:10022–33. 10.3168/jds.2018-15217
61.
FerenčakovićMHamzićEGredlerBSolbergTRKlemetsdalGCurikIet al. Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations. J Anim Breed Genet. (2013) 130:286–93. 10.1111/jbg.12012
62.
BaesCFMakanjuolaBOMigliorFMarrasGHowardJTFlemingAet al. Symposium review: the genomic architecture of inbreeding: how homozygosity affects health and performance. J Dairy Sci. (2019) 102:2807–17. 10.3168/jds.2018-15520
63.
HillWG. Understanding and using quantitative genetic variation. Philos Trans R Soc B Biol Sci. (2010) 365:73–85. 10.1098/rstb.2009.0203
64.
WeigelKA. Controlling inbreeding in modern breeding programs. J Dairy Sci. (2001) 84:E177–84. 10.3168/jds.S0022-0302(01)70213-5
65.
StachowiczKSargolzaeiMMigliorFSchenkelFS. Rates of inbreeding and genetic diversity in Canadian Holstein and Jersey cattle. J Dairy Sci. (2011) 94:5160–75. 10.3168/jds.2010-3308
66.
de RoosAPWHayesBJSpelmanRJGoddardME. Linkage disequilibrium and persistence of phase in Holstein–Friesian, Jersey and Angus Cattle. Genetics. (2008) 179:1503–12. 10.1534/genetics.107.084301
67.
PowellRLNormanHDSandersAH. Progeny testing and selection intensity for holstein bulls in different countries. J Dairy Sci. (2003) 86:3386–93. 10.3168/jds.S0022-0302(03)73942-3
68.
SherfB. The Second Report on the State of the World's Animal Genetic Resources for Food and Agriculture. FAO (2015).
69.
Danchin-BurgeCHiemstraSJBlackburnH. Ex situ conservation of Holstein-Friesian cattle: comparing the Dutch, French, and US germplasm collections. J Dairy Sci. (2011) 94:4100–8. 10.3168/jds.2010-3957
70.
DaetwylerHDVillanuevaBBijmaPWoolliamsJA. Inbreeding in genome-wide selection. J Anim Breed Genet. (2007) 124:369–76. 10.1111/j.1439-0388.2007.00693.x
71.
GianolaDRosaGJM. One hundred years of statistical developments in animal breeding. Annu Rev Anim Biosci. (2015) 3:19–56. 10.1146/annurev-animal-022114-110733
72.
SchaefferLR. Strategy for applying genome-wide selection in dairy cattle. J Anim Breed Genet. (2006) 123:218–23. 10.1111/j.1439-0388.2006.00595.x
73.
ScherfBDPillingD. The Second Report on the State of the World's Animal Genetic Resources for Food Agriculture. Rome (2015). Available online at: http://www.fao.org/3/a-i4787e/index.html (accessed May 5, 2021).
74.
Web Anafi Mate - ANAFIBJ. Available online at: http://www.anafi.it/it/servizi/web-anafi-mate (accessed November 20, 2021).
75.
OECD-FAO Agricultural Outlook 2019-2028. OECD (2019).
76.
MeuwissenTHSonessonAK. Maximizing the response of selection with a predefined rate of inbreeding: overlapping generations. J Anim Sci. (1998) 76:2575. 10.2527/1998.76102575x
77.
MalteccaCTiezziFColeJBBaesC. Symposium review: exploiting homozygosity in the era of genomics—Selection, inbreeding, and mating programs. J Dairy Sci. (2020) 103:5302–13. 10.3168/jds.2019-17846
Summary
Keywords
inbreeding, runs of homozygosity, effective population size, cattle, genomic selection, sustainability
Citation
Ablondi M, Sabbioni A, Stocco G, Cipolat-Gotet C, Dadousis C, Kaam J-T, Finocchiaro R and Summer A (2022) Genetic Diversity in the Italian Holstein Dairy Cattle Based on Pedigree and SNP Data Prior and After Genomic Selection. Front. Vet. Sci. 8:773985. doi: 10.3389/fvets.2021.773985
Received
10 September 2021
Accepted
30 November 2021
Published
13 January 2022
Volume
8 - 2021
Edited by
Eveline M. Ibeagha-Awemu, Agriculture and Agri-Food Canada (AAFC), Canada
Reviewed by
Hugo H. Montaldo, National Autonomous University of Mexico, Mexico; Carina Visser, University of Pretoria, South Africa; Este van Marle-Köster, University of Pretoria, South Africa
Updates
Copyright
© 2022 Ablondi, Sabbioni, Stocco, Cipolat-Gotet, Dadousis, Kaam, Finocchiaro and Summer.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Claudio Cipolat-Gotet claudio.cipolatgotet@unipr.it
This article was submitted to Livestock Genomics, a section of the journal Frontiers in Veterinary Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.