Over the past 100 yr, the range of traits considered for genetic selection in dairy cattle populations has progressed to meet the demands of both industry and society. At the turn of the 20th century, dairy farmers were interested in increasing milk production; however, a systematic strategy for selection was not available. Organized milk performance recording took shape, followed quickly by conformation scoring. Methodological advances in both genetic theory and statistics around the middle of the century, together with technological innovations in computing, paved the way for powerful multitrait analyses. As more sophisticated analytical techniques for traits were developed and incorporated into selection programs, production began to increase rapidly, and the wheels of genetic progress began to turn. By the end of the century, the focus of selection had moved away from being purely production oriented toward a more balanced breeding goal. This shift occurred partly due to increasing health and fertility issues and partly due to societal pressure and welfare concerns. Traits encompassing longevity, fertility, calving, health, and workability have now been integrated into selection indices. Current research focuses on fitness, health, welfare, milk quality, and environmental sustainability, underlying the concentrated emphasis on a more comprehensive breeding goal. In the future, on-farm sensors, data loggers, precision measurement techniques, and other technological aids will provide even more data for use in selection, and the difficulty will lie not in measuring phenotypes but rather in choosing which traits to select for.

Traditionally, pedigree-based relationship coefficients have been used to manage the inbreeding and degree of inbreeding depression that exists within a population. The widespread incorporation of genomic information in dairy cattle genetic evaluations allows for the opportunity to develop and implement methods to manage populations at the genomic level. As a result, the realized proportion of the genome that 2 individuals share can be more accurately estimated instead of using pedigree information to estimate the expected proportion of shared alleles. Furthermore, genomic information allows genome-wide relationship or inbreeding estimates to be augmented to characterize relationships for specific regions of the genome. Region-specific stretches can be used to more effectively manage areas of low genetic diversity or areas that, when homozygous, result in reduced performance across economically important traits. The use of region-specific metrics should allow breeders to more precisely manage the trade-off between the genetic value of the progeny and undesirable side effects associated with inbreeding. Methods tailored toward more effectively identifying regions affected by inbreeding and their associated use to manage the genome at the herd level, however, still need to be developed. We have reviewed topics related to inbreeding, measures of relatedness, genetic diversity and methods to manage populations at the genomic level, and we discuss future challenges related to managing populations through implementing genomic methods at the herd and population levels.

While autozygosity as a consequence of selection is well understood, there is limited information on the ability of different methods to measure true inbreeding. In the present study, a gene dropping simulation was performed and inbreeding estimates based on runs of homozygosity (ROH), pedigree, and the genomic relationship matrix were compared to true inbreeding. Inbreeding based on ROH was estimated using SNP1101, PLINK, and BCFtools software with different threshold parameters. The effects of different selection methods on ROH patterns were also compared. Furthermore, inbreeding coefficients were estimated in a sample of genotyped North American Holstein animals born from 1990 to 2016 using 50 k chip data and ROH patterns were assessed before and after genomic selection. Using ROH with a minimum window size of 20 to 50 using SNP1101 provided the closest estimates to true inbreeding in simulation study. Pedigree inbreeding tended to underestimate true inbreeding, and results for genomic inbreeding varied depending on assumptions about base allele frequencies. Using an ROH approach also made it possible to assess the effect of population structure and selection on distribution of runs of autozygosity across the genome. In the simulation, the longest individual ROH and the largest average length of ROH were observed when selection was based on best linear unbiased prediction (BLUP), whereas genomic selection showed the largest number of small ROH compared to BLUP estimated breeding values (BLUP-EBV). In North American Holsteins, the average number of ROH segments of 1 Mb or more per individual increased from 57 in 1990 to 82 in 2016. The rate of increase in the last 5 years was almost double that of previous 5 year periods. Genomic selection results in less autozygosity per generation, but more per year given the reduced generation interval. This study shows that existing software based on the measurement of ROH can accurately identify autozygosity across the genome, provided appropriate threshold parameters are used. Our results show how different selection strategies affect the distribution of ROH, and how the distribution of ROH has changed in the North American dairy cattle population over the last 25 years.

Genetic diversity in livestock populations is a significant contributor to the sustainability of animal production. Also, genetic diversity allows animal production to become more responsive to environmental changes and market demands. The loss of genetic diversity can result in a plateau in production and may also result in loss of fitness or viability in animal production. In this study, we investigated the rate of inbreeding (ΔF), rate of coancestry (Δf), and effective population size (Ne) as important quantitative indicators of genetic diversity and evaluated the effect of the recent implementation of genomic selection on the loss of genetic diversity in North American Holstein and Jersey dairy cattle. To estimate the rate of inbreeding and coancestry, inbreeding and coancestry coefficients were calculated using the traditional pedigree method and genomic methods estimated from segment- and marker-based approaches. Furthermore, we estimated Ne from the rate of inbreeding and coancestry and extent of linkage disequilibrium. A total of 205,755 and 89,238 pedigreed and genotyped animals born between 1990 and 2018 inclusively were available for Holsteins and Jerseys, respectively. The estimated average pedigree inbreeding coefficients were 7.74 and 7.20% for Holsteins and Jerseys, respectively. The corresponding values for the segment and marker-by-marker genomic inbreeding coefficients were 13.61, 15.64, and 31.40% for Holsteins and 21.16, 22.54, and 42.62% for Jerseys, respectively. The average coancestry coefficients were 8.33 and 15.84% for Holsteins and 9.23 and 23.46% for Jerseys with pedigree and genomic measures, respectively. Generation interval for the whole 29-yr time period averaged approximately 5 yr for all selection pathways combined. The ΔF per generation based on pedigree, segment, and marker-by-marker genomic measures for the entire 29-yr period was estimated to be 0.75, 1.10, 1.16, and 1.02% for Holstein animals and 0.67, 0.62, 0.63, and 0.59% for Jersey animals, respectively. The Δf was estimated to be 0.98 and 0.98% for Holsteins and 0.73 and 0.78% for Jerseys with pedigree and genomic measures, respectively. These ΔF and Δf translated to an Ne that ranged from 43 to 66 animals for Holsteins and 64 to 85 animals for Jerseys. In addition, the Ne based on linkage disequilibrium was 58 and 120 for Holsteins and Jerseys, respectively. The 10-yr period that involved the application of genomic selection resulted in an increased ΔF per generation with ranges from 1.19 to 2.06% for pedigree and genomic measures in Holsteins. Given the rate at which inbreeding is increasing after the implementation of genomic selection, there is a need to implement measures and means for controlling the rate of inbreeding per year, which will help to manage and maintain farm animal genetic resources.

Inbreeding depression is a growing concern in livestock because it can detrimentally affect animal fitness, health, and production levels. Genomic information can be used to more effectively capture variance in Mendelian sampling, thereby enabling more accurate estimation of inbreeding, but further progress is still required. The calculation of inbreeding for herd management purposes is largely still done using pedigree information only, although inbreeding coefficients calculated in this manner have been shown to be less accurate than genomic inbreeding measures. Continuous stretches of homozygous genotypes, so called runs of homozygosity, have been shown to provide a better estimate of autozygosity at the genomic level than conventional measures based on inbreeding coefficients calculated through conventional pedigree information or even genomic relationship matrices. For improved and targeted management of genomic inbreeding at the population level, the development of methods that incorporate genomic information in mate selection programs may provide a more precise tool for reducing the detrimental effects of inbreeding in dairy herds. Additionally, a better understanding of the genomic architecture of inbreeding and incorporating that knowledge into breeding programs could significantly refine current practices. Opportunities to maintain high levels of genetic progress in traits of interest while managing homozygosity and sustaining acceptable levels of heterozygosity in highly selected dairy populations exist and should be examined more closely for continued sustainability of both the dairy cattle population as well as the dairy industry. The inclusion of precise genomic measures of inbreeding, such as runs of homozygosity, inbreeding, and mating programs, may provide a path forward. In this symposium review article, we describe traditional measures of inbreeding and the recent developments made toward more precise measures of homozygosity using genomic information. The effects of homozygosity resulting from inbreeding on phenotypes, the identification and mapping of detrimental homozygosity haplotypes, management of inbreeding with genomic data, and areas in need of further research are discussed.

The success and sustainability of a breeding program incorporating genomic information is largely dependent on the accuracy of predictions. For low heritability traits, large training populations are required to achieve high accuracies of genomic estimated breeding values (GEBV). By including genotyped and nongenotyped animals simultaneously in the evaluation, the single-step genomic BLUP (ssGBLUP) approach has the potential to deliver more accurate and less biased genomic evaluations. The aim of this study was to compare the accuracy and bias of genomic predictions for various traits in Canadian Holstein cattle using ssGBLUP and multi-step genomic BLUP (msGBLUP) under different strategies, such as (1) adding genomic information of cows in the analysis, (2) testing different adjustments of the genomic relationship matrix, and (3) using a blending approach to obtain GEBV from msGBLUP. The following genomic predictions were evaluated regarding accuracy and bias: (1) GEBV estimated by ssGBLUP; (2) direct genomic value estimated by msGBLUP with polygenic effects of 5 and 20%; and (3) GEBV calculated by a blending approach of direct genomic value with estimated breeding values using polygenic effects of 5 and 20%. The effect of adding genomic information of cows in the evaluation was also assessed for each approach. When genomic information was included in the analyses, the average improvement in observed reliability of predictions was observed to be 7 and 13 percentage points for reproductive and workability traits, respectively, compared with traditional BLUP. Absolute deviation from 1 of the regression coefficient of the linear regression of de-regressed estimated breeding values on genomic predictions went from 0.19 when using traditional BLUP to 0.22 when using the msGBLUP method, and to 0.14 when using the ssGBLUP method. The use of polygenic weight of 20% in the msGBLUP slightly improved the reliability of predictions, while reducing the bias. A similar trend was observed when a blending approach was used. Adding genomic information of cows increased reliabilities, while decreasing bias of genomic predictions when using the ssGBLUP method. Differences between using a training population with cows and bulls or with only bulls for the msGBLUP method were small, likely due to the small number of cows included in the analysis. Predictions for lowly heritable traits benefit greatly from genomic information, especially when all phenotypes, pedigrees, and genotypes are used in a single-step approach.