Estimating autozygosity from high-throughput information: effects of SNP density and genotyping errors

How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species

Background

PLINK is probably the most used program for analyzing SNP genotypes and runs of homozygosity (ROH), both in human and in animal populations. The last decade, ROH analyses have become the state-of-the-art method for inbreeding assessment. In PLINK, the --homozyg function is used to perform ROH analyses and relies on several input settings. These settings can have a large impact on the outcome and default values are not always appropriate for medium density SNP array data. Guidelines for a robust and uniform ROH analysis in PLINK using medium density data are lacking, albeit these guidelines are vital for comparing different ROH studies. In this study, 8 populations of different livestock and pet species are used to demonstrate the importance of PLINK input settings. Moreover, the effects of pruning SNPs for low minor allele frequencies and linkage disequilibrium on ROH detection are shown.

Results

We introduce the genome coverage parameter to appropriately estimate F_ROH and to check the validity of ROH analyses. The effect of pruning for linkage disequilibrium and low minor allele frequencies on ROH analyses is highly population dependent and such pruning may result in missed ROH. PLINK’s minimal density requirement is crucial for medium density genotypes and if set too low, genome coverage of the ROH analysis is limited. Finally, we provide recommendations for the maximal gap, scanning window length and threshold settings.

Conclusions

In this study, we present guidelines for an adequate and robust ROH analysis in PLINK on medium density SNP data. Furthermore, we advise to report parameter settings in publications, and to validate them prior to analysis. Moreover, we encourage authors to report genome coverage to reflect the ROH analysis’ validity. Implementing these guidelines will substantially improve the overall quality and uniformity of ROH analyses.

Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon

Genomic approaches permit direct estimation of inbreeding and its effect on fitness. We used genomic-based estimates of inbreeding to investigate their relationship with eight adult traits in a captive-reared Pacific salmonid that is released into the wild. Estimates were also used to determine whether alternative broodstock management approaches reduced risks of inbreeding. Specifically, 1,100 unlinked restriction-site associated (RAD) loci were used to compare pairwise relatedness, derived from a relationship matrix, and individual inbreeding, estimated by comparing observed and expected homozygosity, across four generations in two hatchery lines of Chinook salmon that were derived from the same source. The lines are managed as "integrated" with the founding wild stock, with ongoing gene flow, and as "segregated" with no gene flow. While relatedness and inbreeding increased in the first generation of both lines, possibly due to population subdivision caused by hatchery initiation, the integrated line had significantly lower levels in some subsequent generations (relatedness: F₂ -F₄ ; inbreeding F₂ ). Generally, inbreeding was similar between the lines despite large differences in effective numbers of breeders. Inbreeding did not affect fecundity, reproductive effort, return timing, fork length, weight, condition factor, and daily growth coefficient. However, it delayed spawn timing by 1.75 days per one standard deviation increase in F (~0.16). The results indicate that integrated management may reduce inbreeding but also suggest that it is relatively low in a small, segregated hatchery population that maximized number of breeders. Our findings demonstrate the utility of genomics to monitor inbreeding under alternative management strategies in captive breeding programs.

Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (FROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (FHOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (FPED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. FHOM was highly correlated (0.671 to 0.985) with FROH measures in all breeds. Apulo-Calabrese and Casertana had the highest FROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (FROH1) to 16 Mbp (FROH16)), whereas the lowest FROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). FROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all FROH measures and FPED (from 0.514 to 0.523) and between FHOM and FPED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between FROH and FPED. This might be due to the imperfect measure of FPED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.

Unraveling the genetic diversity of Belgian Milk Sheep using medium-density SNP genotypes

The present study focuses on the Belgian Milk Sheep in Flanders (Belgium) and compares its genetic diversity and relationship with the Flemish Sheep, the Friesian Milk Sheep, the French Lacaune dairy sheep and other Northern European breeds. For this study, 94 Belgian Milk Sheep, 23 Flemish Sheep and 22 Friesian Milk Sheep were genotyped with the OvineSNP50 array. In addition, 29 unregistered animals phenotypically similar to Belgian Milk Sheep were genotyped using the 15K ISGC chip. Both Belgian and Friesian Milk Sheep as well as the East Friesian Sheep were found to be less diverse than the other seven breeds included in this study. Genomic inbreeding coefficients based on runs of homozygosity (ROH) were estimated at 14.5, 12.4 and 10.2% for Belgian Milk Sheep, Flemish Sheep and Friesian Milk Sheep respectively. Out of 29 unregistered Belgian Milk Sheep, 28 mapped in the registered Belgian Milk Sheep population. Ancestry analysis, PCA and F_ST calculations showed that Belgian Milk Sheep are more related to Friesian Milk Sheep than to Flemish Sheep, which was contrary to the breeders' expectations. Consequently, breeders may prefer to crossbreed Belgian Milk Sheep with Friesian sheep populations (Friesian Milk Sheep or East Friesian Sheep) in order to increase diversity. This research underlines the usefulness of SNP chip genotyping and ROH analyses for monitoring genetic diversity and studying genetic links in small livestock populations, profiting from internationally available genotypes. As assessment of genetic diversity is vital for long-term breed survival, these results will aid flockbooks to preserve genetic diversity.

SNP genotypes reveal breed substructure, selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population‐specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.

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

Simple Summary

Runs of homozygosity (ROH) regions are known to be common in the genomes of cattle and have become a subject of interest of various research in recent years. ROH can be used as a valuable tool to estimate inbreeding, which needs to be controlled in livestock populations. Moreover, analysis of ROH is considered to be an effective method of identifying genome regions that are a subject of selection pressure, which may help in understanding the genetic aspects of production traits under selection. In this study, we analyzed ROH characteristics of 11 cattle breeds, both commercial and native, maintained in Poland. We presented distinct differences in the length, quantity and frequency of ROH between the analyzed breeds as well as in the levels of genomic inbreeding. Higher levels of inbreeding were characteristic for commercial breeds, especially beef breeds. In addition, within ROH islands, we observed a number of genes with a confirmed influence on the level of production traits. The presented results and identified genes can be a basis for further research focused on the identification of genes and markers essential in the determination of the most important production traits in cattle.

Abstract

In the presented research, BovineSNP50 microarrays (Illumina) were applied to determine runs of homozygosity in the genomes of 11 cattle breeds maintained in Poland. These cattle breeds represent three basic utility types: milk, meat and dual purpose. Analysis of runs of homozygosity allowed the evaluation of the level of autozygosity within each breed in order to calculate the genomic inbreeding coefficient (F_ROH), as well as to identify regions of the genome with a high frequency of ROH occurrence, which may reflect traces of directional selectin left in their genomes. Visible differences in the length and distribution of runs of homozygosity in the genomes of the analyzed cattle breeds have been observed. The highest mean number and mean sums of lengths of runs of homozygosity were characteristic for Hereford cattle and intermediate for the Holstein-Friesian Black-and-White variety, Holstein-Friesian Red-and-White variety, Simmental, Limousin, Montbeliarde and Charolais breeds. However, lower values were observed for cattle of conserved breeds. Moreover, the selected livestock differed in the level of inbreeding estimated using the F_ROH coefficient. In regions of the genome with a high frequency of ROH occurrence, which may reflect the impact of directional selection, a number of genes were observed that can be potentially related to the production traits which are under selection pressure for specific production types. The most important detected genes were GHR, MSTN, DGAT1, FABP4, and TRH, with a known influence on the milk and meat traits of the studied cattle breeds.

Genomic Patterns of Homozygosity in Chinese Local Cattle

Genome-wide single nucleotide polymorphism (SNP) arrays can be used to explore homozygosity segments, where two haplotypes inherited from the parents are identical. In this study, we identified a total of 27,358 runs of homozygosity (ROH) with an average of 153 ROH events per animal in Chinese local cattle. The sizes of ROH events varied considerably ranging from 0.5 to 66 Mb, with an average length of 1.22 Mb. The highest average proportion of the genome covered by ROH (~11.54% of the cattle genome) was found in Nanda cattle (NDC) from South China, whereas the lowest average proportion (~3.1%) was observed in Yanhuang cattle (YHC). The average estimated F_ROH ranged from 0.03 in YHC to 0.12 in NDC. For each of three ROH classes with different sizes (Small 0.5-1 Mb, Medium 1-5 Mb and Large >5 Mb), the numbers and total lengths of ROH per individual showed considerable differences across breeds. Moreover, we obtained 993 to 3603 ROH hotspots (which were defined where ROH frequency at a SNP within each breed exceeded the 1% threshold) among eight cattle breeds. Our results also revealed several candidate genes embedded with ROH hotspots which may be related to environmental conditions and local adaptation. In conclusion, we generated baselines for homozygosity patterns in diverse Chinese cattle breeds. Our results suggested that selection has, at least partially, played a role with other factors in shaping the genomic patterns of ROH in Chinese local cattle and might provide valuable insights for understanding the genetic basis of economic and adaptive traits.

A Population Genomics Analysis of the Native Irish Galway Sheep Breed

The Galway sheep population is the only native Irish sheep breed and this livestock genetic resource is currently categorised as 'at-risk'. In the present study, comparative population genomics analyses of Galway sheep and other sheep populations of European origin were used to investigate the microevolution and recent genetic history of the breed. These analyses support the hypothesis that British Leicester sheep were used in the formation of the Galway. When compared to conventional and endangered breeds, the Galway breed was intermediate in effective population size, genomic inbreeding and runs of homozygosity. This indicates that, although the Galway breed is declining, it is still relatively genetically diverse and that conservation and management plans informed by genomic information may aid its recovery. The Galway breed also exhibited distinct genomic signatures of artificial or natural selection when compared to other breeds, which highlighted candidate genes that may be involved in production and health traits.

Inbreeding depression due to recent and ancient inbreeding in Dutch Holstein-Friesian dairy cattle

Background

Inbreeding decreases animal performance (inbreeding depression), but not all inbreeding is expected to be equally harmful. Recent inbreeding is expected to be more harmful than ancient inbreeding, because selection decreases the frequency of deleterious alleles over time. Selection efficiency is increased by inbreeding, a process called purging. Our objective was to investigate effects of recent and ancient inbreeding on yield, fertility and udder health traits in Dutch Holstein–Friesian cows.

Methods

In total, 38,792 first-parity cows were included. Pedigree inbreeding (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{PED}$$\end{document}FPED) was computed and 75 k genotype data were used to compute genomic inbreeding, among others based on regions of homozygosity (ROH) in the genome (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{ROH}$$\end{document}FROH).

Results

Inbreeding depression was observed, e.g. a 1% increase in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{ROH}$$\end{document}FROH was associated with a 36.3 kg (SE = 2.4) decrease in 305-day milk yield, a 0.48 day (SE = 0.15) increase in calving interval and a 0.86 unit (SE = 0.28) increase in somatic cell score for day 150 through to 400. These effects equalled − 0.45, 0.12 and 0.05% of the trait means, respectively. When \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{PED}$$\end{document}FPED was split into generation-based components, inbreeding on recent generations was more harmful than inbreeding on more distant generations for yield traits. When \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{PED}$$\end{document}FPED was split into new and ancestral components, based on whether alleles were identical-by-descent for the first time or not, new inbreeding was more harmful than ancestral inbreeding, especially for yield traits. For example, a 1% increase in new inbreeding was associated with a 2.42 kg (SE = 0.41) decrease in 305-day fat yield, compared to a 0.03 kg (SE = 0.71) increase for ancestral inbreeding. There were no clear differences between effects of long ROH (recent inbreeding) and short ROH (ancient inbreeding).

Conclusions

Inbreeding depression was observed for yield, fertility and udder health traits. For yield traits and based on pedigree, inbreeding on recent generations was more harmful than inbreeding on distant generations and there was evidence of purging. Across all traits, long and short ROH contributed to inbreeding depression. In future work, inbreeding depression and purging should be assessed in more detail at the genomic level, using higher density information and genomic time series.

The impact of genomic selection on genetic diversity and genetic gain in three French dairy cattle breeds

Background

In France, implementation of genomic evaluations in dairy cattle breeds started in 2009 and this has modified the breeding schemes drastically. In this context, the goal of our study was to understand the impact of genomic selection on the genetic diversity of bulls from three French dairy cattle breeds born between 2005 and 2015 (Montbéliarde, Normande and Holstein) and the factors that are involved.

Methods

We compared annual genetic gains, inbreeding rates based on runs of homozygosity (ROH) and pedigree data, and mean ROH length within breeds, before and after the implementation of genomic selection.

Results

Genomic selection induced an increase in mean annual genetic gains of 50, 71 and 33% for Montbéliarde, Normande and Holstein bulls, respectively, and in parallel, the generation intervals were reduced by a factor of 1.7, 1.9 and 2, respectively. We found no significant change in inbreeding rate for the two national breeds, Montbéliarde and Normande, and a significant increase in inbreeding rate for the Holstein international breed, which is now as high as 0.55% per year based on ROH and 0.49% per year based on pedigree data (equivalent to a rate of 1.36 and 1.39% per generation, respectively). The mean ROH length was longer for bulls from the Holstein breed than for those from the other two breeds.

Conclusions

With the implementation of genomic selection, the annual genetic gain increased for bulls from the three major French dairy cattle breeds. At the same time, the annual loss of genetic diversity increased for Holstein bulls, possibly because of the massive use of a few elite bulls in this breed, but not for Montbéliarde and Normande bulls. The increase in mean ROH length in Holstein may reflect the occurrence of recent inbreeding. New strategies in breeding schemes, such as female donor stations and embryo transfer, and recent implementation of genomic evaluations in small regional breeds should be studied carefully in order to ensure the sustainability of breeding schemes in the future.

Key Genetic Parameters for Population Management

Population management has the primary task of maximizing the long-term competitiveness of a breed. Breeds compete with each other for being able to supply consumer demands at low costs and also for funds from conservation programs. The competition for consumer preference is won by breeds with high genetic gain for total merit who maintained a sufficiently high genetic diversity, whereas the competition for funds is won by breeds with high conservation value. The conservation value of a breed could be improved by increasing its contribution to the gene pool of the species. This may include the recovery of its original genetic background and the maintenance of a high genetic diversity at native haplotype segments. The primary objective of a breeding program depends on the genetic state of the population and its intended usage. In this paper, we review the key genetic parameters that are relevant for population management, compare the methods for estimating them, derive the formulas for predicting their value at a future time, and clarify their usage in various types of breeding programs that differ in their main objectives. These key parameters are kinships, native kinships, breeding values, Mendelian sampling variances, native contributions, and mutational effects. Population management currently experiences a transition from using pedigree-based estimates to marker-based estimates, which improves the accuracies of these estimates and thereby increases response to selection. In addition, improved measures of the factors that determine the competitiveness of a breed and utilize auxiliary parameters, such as Mendelian sampling variances, mutational effects, and native kinships, enable to improve further upon historic recommendations for genetic population management.

A mating advice system in dairy cattle incorporating genomic information

This study investigated the effects of alternative mating programs that incorporate genomic information on expected progeny herd performance and inbreeding, as well as methods to include un-genotyped animals in such mating programs. A total of 54,535 Holstein-Friesian cattle with imputed high-density genotypes (547,650 SNP after edits) were available. First, to quantify the accuracy of imputing un-genotyped animals (often an issue in populations), a sub-population of 729 genotyped animals had their genotypes masked, and their allele dosages were imputed, using linear regression exploiting information on genotyped relatives. The reference population for imputation included all genotyped animals, excluding the 729 selected animals and their sires, dams, and grandsires, and had either (1) their sires' genotypes, (2) their dams' genotypes (3) both their sires' and their dams' genotypes, or (4) both their sires' and maternal grandsires' genotypes introduced into the reference population. The correlations between true genotypes and the imputed allele dosages ranged from 0.58 (sire only) to 0.68 (both sire and dam). A herd of 100 cows was then simulated (1,000 replicates) from the sub-population of 729 imputed animals. The top 10 bulls from the genotyped population, based on their total genetic merit index (TMI) were selected to be used as sires. Three mating allotment methods were investigated: (1) random mating, (2) sequential mating based on maximizing only the expected TMI of the progeny, and (3) linear programming to maximize a generated index constructed to maximize genetic merit and minimize expected progeny inbreeding as well as intra- and inter-progeny variability in genetic merit. Relationships among candidate parents were calculated using either the pedigree relationship matrix or the genomic relationship matrix; the latter was constructed using either the true genotypes of both parents or the true genotypes of the sire plus the imputed allele dosages of the dam. Using the genomic co-ancestry estimates resulted in lower average herd expected genomic inbreeding levels compared with using the pedigree-based co-ancestry estimates. Additionally, if the dams were not genotyped, using their imputed allele dosages also resulted in lower average herd expected inbreeding levels compared with using the pedigree co-ancestry estimates. The inter-progeny coefficient of variation for selected traits, milk and fertility, estimated breeding values were reduced by 12 to 65% using the linear programing method compared with sequential mating.

Timing and Extent of Inbreeding in African Goats

Genetic characterization of African goats is one of the current priorities in the improvement of goats in the continent. This study contributes to the characterization effort by determining the levels and number of generations to common ancestors ("age") associated with inbreeding in African goat breeds and identifies regions that contain copy number variation mistyped as being homozygous. Illumina 50k single nucleotide polymorphism genotype data for 608 goats from 31 breeds were used to compute the level and age of inbreeding at both local (marker) and global levels (F_G) using a model-based approach based on a hidden Markov model. Runs of homozygosity (ROH) segments detected using the Viterbi algorithm led to ROH-based inbreeding coefficients for all ROH (F_ROH) and for ROH longer than 2 Mb (F_{ROH > 2Mb}). Some of the genomic regions identified as having ROH are likely to be hemizygous regions (copy number deletions) mistyped as homozygous regions. Although the proportion of these miscalled ROH is small and does not substantially affect estimates of levels of inbreeding for individual animals, the inbreeding metrics were adjusted by removing these regions from the ROH. All the inbreeding metrics varied widely across breeds, with overall means of 0.0408, 0.0370, and 0.0691 and medians of 0.0125, 0.0098, and 0.0366 for F_ROH, F_{ROH > 2Mb}, and F_G, respectively. Several breeds (including Menabe and Sofia from Madagascar) had high proportions of recent inbreeding, while Small East African, Ethiopian, and most of the West African breeds (including West African Dwarf) had more ancient inbreeding.

Genomic measures of inbreeding in the Norwegian-Swedish Coldblooded Trotter and their associations with known QTL for reproduction and health traits

Background

Since the 1950s, the Norwegian–Swedish Coldblooded trotter (NSCT) has been intensively selected for harness racing performance. As a result, the racing performance of the NSCT has improved remarkably; however, this improved racing performance has also been accompanied by a gradual increase in inbreeding level. Inbreeding in NSCT has historically been monitored by using traditional methods that are based on pedigree analysis, but with recent advancements in genomics, the NSCT industry has shown interest in adopting molecular approaches for the selection and maintenance of this breed. Consequently, the aims of the current study were to estimate genomic-based inbreeding coefficients, i.e. the proportion of runs of homozygosity (ROH), for a sample of NSCT individuals using high-density genotyping array data, and subsequently to compare the resulting rate of genomic-based F (F_ROH) to that of pedigree-based F (F_PED) coefficients within the breed.

Results

A total of 566 raced NSCT were available for analyses. Average F_ROH ranged from 1.78 to 13.95%. Correlations between F_ROH and F_PED were significant (P < 0.001) and ranged from 0.27 to 0.56, with F_PED and F_ROH from 2000 to 2009 increasing by 1.48 and 3.15%, respectively. Comparisons of ROH between individuals yielded 1403 regions that were present in at least 95% of the sampled horses. The average percentage of a single chromosome covered in ROH ranged from 9.84 to 18.82% with chromosome 31 and 18 showing, respectively, the largest and smallest amount of homozygosity.

Conclusions

Genomic inbreeding coefficients were higher than pedigree inbreeding coefficients with both methods showing a gradual increase in inbreeding level in the NSCT breed between 2000 and 2009. Opportunities exist for the NSCT industry to develop programs that provide breeders with easily interpretable feedback on regions of the genome that are suboptimal from the perspective of genetic merit or that are sensitive to inbreeding within the population. The use of molecular data to identify genomic regions that may contribute to inbreeding depression in the NSCT will likely prove to be a valuable tool for the preservation of its genetic diversity in the long term.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0465-7) contains supplementary material, which is available to authorized users.

Selection signatures in four German warmblood horse breeds: Tracing breeding history in the modern sport horse

The study of selection signatures helps to find genomic regions that have been under selective pressure and might host genes or variants that modulate important phenotypes. Such knowledge improves our understanding of how breeding programmes have shaped the genomes of livestock. In this study, 942 stallions were included from four, exemplarily chosen, German warmblood breeds with divergent historical and recent selection focus and different crossbreeding policies: Trakehner (N = 44), Holsteiner (N = 358), Hanoverian (N = 319) and Oldenburger (N = 221). Those breeds are nowadays bred for athletic performance and aptitude for show-jumping, dressage or eventing, with a particular focus of Holsteiner on the first discipline. Blood samples were collected during the health exams of the stallion preselections before licensing and were genotyped with the Illumina EquineSNP50 BeadChip. Autosomal markers were used for a multi-method search for signals of positive selection. Analyses within and across breeds were conducted by using the integrated Haplotype Score (iHS), cross-population Extended Haplotype Homozygosity (xpEHH) and Runs of Homozygosity (ROH). Oldenburger and Hanoverian showed very similar iHS signatures, but breed specificities were detected on multiple chromosomes with the xpEHH. The Trakehner clustered as a distinct group in a principal component analysis and also showed the highest number of ROHs, which reflects their historical bottleneck. Beside breed specific differences, we found shared selection signals in an across breed iHS analysis on chromosomes 1, 4 and 7. After investigation of these iHS signals and shared ROH for potential functional candidate genes and affected pathways including enrichment analyses, we suggest that genes affecting muscle functionality (TPM1, TMOD2-3, MYO5A, MYO5C), energy metabolism and growth (AEBP1, RALGAPA2, IGFBP1, IGFBP3-4), embryonic development (HOXB-complex) and fertility (THEGL, ZPBP1-2, TEX14, ZP1, SUN3 and CFAP61) have been targeted by selection in all breeds. Our findings also indicate selection pressure on KITLG, which is well-documented for influencing pigmentation.

Linkage disequilibrium and within-breed genetic diversity in Iranian Zandi sheep

This research aimed to measure the extent of linkage disequilibrium (LD), effective population size (Ne), and runs of homozygosity (ROHs) in one of the major Iranian sheep breeds (Zandi) using 96 samples genotyped with Illumina Ovine SNP50 BeadChip. The amount of LD (r2) for single-nucleotide polymorphism (SNP) pairs in short distances (10–20 kb) was 0.21±0.25 but rapidly decreased to 0.10±0.16 by increasing the distance between SNP pairs (40–60 kb). The Ne of Zandi sheep in past (approximately 3500 generations ago) and recent (five generations ago) populations was estimated to be 6475 and 122, respectively. The ROH-based inbreeding was 0.023. We found 558 ROH regions, of which 37 % were relatively long (>10 Mb). Compared with the rate of LD reduction in other species (e.g., cattle and pigs), in Zandi, it was reduced more rapidly by increasing the distance between SNP pairs. According to the LD pattern and high genetic diversity of Zandi sheep, we need to use an SNP panel with a higher density than Illumina Ovine SNP50 BeadChip for genomic selection and genome-wide association studies in this breed.

Assessment of Autozygosity Derived From Runs of Homozygosity in Jinhua Pigs Disclosed by Sequencing Data

Jinhua pig, a well-known Chinese indigenous breed, has evolved as a pig breed with excellent meat quality, greater disease resistance, and higher prolificacy. The reduction in the number of Jinhua pigs over the past years has raised concerns about inbreeding. Runs of homozygosity (ROH) along the genome have been applied to quantify individual autozygosity to improve the understanding of inbreeding depression and identify genes associated with traits of interest. Here, we investigated the occurrence and distribution of ROH using next-generation sequencing data to characterize autozygosity in 202 Jinhua pigs, as well as to identify the genomic regions with high ROH frequencies within individuals. The average inbreeding coefficient, based on ROH longer than 1 Mb, was 0.168 ± 0.052. In total, 18,690 ROH were identified in all individuals, among which shorter segments (1-5 Mb) predominated. Individual ROH autosome coverage ranged from 5.32 to 29.14% in the Jinhua population. On average, approximately 16.8% of the whole genome was covered by ROH segments, with the lowest coverage on SSC11 and the highest coverage on SSC17. A total of 824 SNPs (about 0.5%) and 11 ROH island regions were identified (occurring in over 45% of the samples). Genes associated with reproduction (HOXA3, HOXA7, HOXA10, and HOXA11), meat quality (MYOD1, LPIN3, and CTNNBL1), appetite (NUCB2) and disease resistance traits (MUC4, MUC13, MUC20, LMLN, ITGB5, HEG1, SLC12A8, and MYLK) were identified in ROH islands. Moreover, several quantitative trait loci for ham weight and ham fat thickness were detected. Genes in ROH islands suggested, at least partially, a selection for economic traits and environmental adaptation, and should be subject of future investigation. These findings contribute to the understanding of the effects of environmental and artificial selection in shaping the distribution of functional variants in the pig genome.

Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates

Before availability of dense SNP data, genetic diversity was characterized and managed with pedigree-based information. Besides this classical approach, 2 methodologies have been proposed in recent years to characterize and manage diversity from dense SNP data: the SNP-by-SNP approach and the alternative based on runs of homozygosity (ROH). The establishment of criteria to identify ROH is a current constraint in the literature dealing with ROH. The objective of this study was, using a medium-density SNP chip, to quantify by 3 methods (pedigree, SNP-by-SNP, and ROH) the genetic diversity on 5 selected French dairy sheep subpopulations and breeds and to assess the effect of the definition of ROH on these estimates. The data set available included individuals from the breeds Basco-Béarnaise, Manech Tête Noire, Manech Tête Rousse, and 2 subpopulations of Lacaune: Lacaune Confederation and Lacaune Ovitest. Animals were genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). After filtering, the genomic data included 38,287 autosomal SNP and 8,700 individuals, which comprised 72,803 animals in the pedigree. The results indicated that no significant differences were observed in effective population size estimates obtained from pedigree or genomic (SNP-by-SNP or ROH) information. In general, estimates of effective population size were above 200 in Lacaune Confederation and Lacaune Ovitest subpopulations and below 200 in Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds. The minimum length that constituted a ROH, the minimum number of SNP that constituted a ROH, as well as the minimum density and the maximum distance allowed between 2 homozygous SNP are ROH-defining factors with important implications in the estimation of the rate of inbreeding. The ROH-based rates of inbreeding in concordance with those obtained from pedigree information require a specific set of values. This particular set of values is different from that identified to obtain ROH-based rates of inbreeding similar to those obtained on a SNP-by-SNP basis. Factors to define ROH do not change the results much unless extreme values are considered, although further research on ROH-based inbreeding is still required.

Inbreeding depression in line 1 Hereford cattle population using pedigree and genomic information

This study aimed at assessing inbreeding and its effect on growth and fertility traits using the longtime closed line 1 Hereford cattle population. Inbreeding was estimated based on pedigree (FPED) and genomic information. For the latter, three estimates were derived based on the diagonal elements of the genomic relationship matrix using estimated (FGRM) or fixed (FGRM0.5) minor allele frequencies or runs of homozygosity (ROH) (FROH). A pedigree containing 10,186 animals was used to calculate FPED. Genomic inbreeding was evaluated using 785 animals genotyped for 30,810 SNP. Traits analyzed were birth weight (BWT), weaning weight (WWT), yearling weight (YWT), ADG, and age at first calving (AFC). The number of ROH per animal ranged between 6 and 119 segments with an average of 83. The shortest and longest segments were 1.36 and 64.86 Mb long, respectively, reflecting both ancient and recent inbreeding occurring in the last 30 to 40 generations. The average inbreeding was 29.2%, 16.1%, 30.2%, and 22.9% for FPED, FGRM, FGRM0.5, and FROH, respectively. FROH provided the highest correlations with FPED (r = 0.66). Across paternal half-sib families, with minimal variation in FPED, there were substantial variations in their genomic inbreeding. Inbreeding depression analyses showed that a 1% increase in an animal's FPED resulted in a decrease of 1.20 kg, 2.03 kg, and 0.004 kg/d in WWT, YWT, and ADG, respectively. Maternal inbreeding showed significantly negative effects on progeny growth performance. AFC increased by 1.4 and 0.8 d for each 1% increase in FPED of the cow and her dam, respectively. Using genomic inbreeding, similar impact on growth traits was observed although the magnitude of the effect varied between methods. Across all genomic measures, WWT, YWT, and ADG decreased by 0.21 to 0.53 kg, 0.46 to 1.13 kg, and 0.002 to 0.006 kg/d for each 1% increase in genomic inbreeding, respectively. Four chromosomes (9, 12, 17, and 27) were identified to have a significant association between their homozygosity (FROH-CHR) and growth traits. Variability in genomic inbreeding could be useful when deciding between full and half-sib selection candidates. Despite the high level of inbreeding in this study, its negative impact on growth performance was not as severe as expected, which may be attributed to the purging of the deleterious alleles due to natural or artificial selection over time.

Linkage disequilibrium and past effective population size in native Tunisian cattle

To carry out effective genome-wide association studies, information about linkage disequilibrium (LD) is essential. Here, we used medium-density SNP chips to provide estimates of LD in native Tunisian cattle. The two measures of LD that were used, mean r2 and D', decreased from 0.26 to 0.05 and from 0.73 to 0.40, respectively, when the distance between markers increased from less than 20 Kb to 200 Kb. The decay in LD over physical distance occurred at a faster rate than that reported for European and other indigenous breeds, and reached background levels at less than 500 Kb distance. This is consistent with the absence of strong selective pressure within the Tunisian population and suggests that, in order to be effective, any potential genome-wide association mapping studies will need to use chips with higher marker density. An analysis of effective population size (Ne) based on LD data showed a decline in past Ne, with a sudden drop starting about eight generations ago. This finding, combined with the high levels of recent inbreeding revealed by runs of homozygosity (ROH) analysis, indicate that this population is endangered and may be in urgent need of a conservation plan that includes a well-designed genetic management program.

Optimum contribution selection for animal breeding and conservation: the R package optiSel

Background

Selecting animals for breeding in the optimum way plays an essential role for the management of genetic resources and in selective breeding of livestock species. It requires to compute the optimum genetic contribution of each selection candidate to the next generation. Current software packages for optimum contribution selection (OCS) are not able to handle the main conflicting objectives of animal breeding programs simultaneously, which includes to increase genetic gain, to increase or to maintain genetic diversity, to recover the original genetic background of endangered breeds with historic introgression, and to maintain or increase genetic diversity at native alleles.

Results

The free R package optiSel offers functions for estimating the above mentioned parameters from pedigree and marker data, and for solving OCS problems. One parameter can be optimized, whereas the remaining ones can be constrained. The results reveal the optimum numbers of offspring of all selection candidates, and can subsequently be used for mate allocation. Different solvers can be used. Solver slsqp was superior when the genetic diversity at native alleles was to be maximized, whereas solvers cccp and cccp2 were superior for all other OCS problems.

Conclusion

Optimum contribution selection applied to local breeds requires special attention due to the conflicting objectives of their breeding programs. The free R package optiSel is an easy-to-use software taking these conflicting objectives into account.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2450-5) contains supplementary material, which is available to authorized users.

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances

Background

Runs of homozygosity (ROH) islands are stretches of homozygous sequence in the genome of a large proportion of individuals in a population. Algorithms for the detection of ROH depend on the similarity of haplotypes. Coverage gaps and copy number variants (CNV) may result in incorrect identification of such similarity, leading to the detection of ROH islands where none exists. Misidentified hemizygous regions will also appear as homozygous based on sequence variation alone. Our aim was to identify ROH islands influenced by marker coverage gaps or CNV, using Illumina BovineHD BeadChip (777 K) single nucleotide polymorphism (SNP) data for Austrian Brown Swiss, Tyrol Grey and Pinzgauer cattle.

Methods

ROH were detected using clustering, and ROH islands were determined from population inbreeding levels for each marker. CNV were detected using a multivariate copy number analysis method and a hidden Markov model. SNP coverage gaps were defined as genomic regions with intermarker distances on average longer than 9.24 kb. ROH islands that overlapped CNV regions (CNVR) or SNP coverage gaps were considered as potential artefacts. Permutation tests were used to determine if overlaps between CNVR with copy losses and ROH islands were due to chance. Diversity of the haplotypes in the ROH islands was assessed by haplotype analyses.

Results

In Brown Swiss, Tyrol Grey and Pinzgauer, we identified 13, 22, and 24 ROH islands covering 26.6, 389.0 and 35.8 Mb, respectively, and we detected 30, 50 and 71 CNVR derived from CNV by using both algorithms, respectively. Overlaps between ROH islands, CNVR or coverage gaps occurred for 7, 14 and 16 ROH islands, respectively. About 37, 44 and 52% of the ROH islands coverage in Brown Swiss, Tyrol Grey and Pinzgauer, respectively, were affected by copy loss. Intersections between ROH islands and CNVR were small, but significantly larger compared to ROH islands at random locations across the genome, implying an association between ROH islands and CNVR. Haplotype diversity for reliable ROH islands was lower than for ROH islands that intersected with copy loss CNVR.

Conclusions

Our findings show that a significant proportion of the ROH islands in the bovine genome are artefacts due to CNV or SNP coverage gaps.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0414-x) contains supplementary material, which is available to authorized users.

Genome-Wide Characterization of Selection Signatures and Runs of Homozygosity in Ugandan Goat Breeds

Both natural and artificial selection are among the main driving forces shaping genetic variation across the genome of livestock species. Selection typically leaves signatures in the genome, which are often characterized by high genetic differentiation across breeds and/or a strong reduction in genetic diversity in regions associated with traits under intense selection pressure. In this study, we evaluated selection signatures and genomic inbreeding coefficients, FROH, based on runs of homozygosity (ROH), in six Ugandan goat breeds: Boer (n = 13), and the indigenous breeds Karamojong (n = 15), Kigezi (n = 29), Mubende (n = 29), Small East African (n = 29), and Sebei (n = 29). After genotyping quality control, 45,294 autosomal single nucleotide polymorphisms (SNPs) remained for further analyses. A total of 394 and 6 breed-specific putative selection signatures were identified across all breeds, based on marker-specific fixation index (FST-values) and haplotype differentiation (hapFLK), respectively. These regions were enriched with genes involved in signaling pathways associated directly or indirectly with environmental adaptation, such as immune response (e.g., IL10RB and IL23A), growth and fatty acid composition (e.g., FGF9 and IGF1), and thermo-tolerance (e.g., MTOR and MAPK3). The study revealed little overlap between breeds in genomic regions under selection and generally did not display the typical classic selection signatures as expected due to the complex nature of the traits. In the Boer breed, candidate genes associated with production traits, such as body size and growth (e.g., GJB2 and GJA3) were also identified. Furthermore, analysis of ROH in indigenous goat breeds showed very low levels of genomic inbreeding (with the mean FROH per breed ranging from 0.8% to 2.4%), as compared to higher inbreeding in Boer (mean FROH = 13.8%). Short ROH were more frequent than long ROH, except in Karamojong, providing insight in the developmental history of these goat breeds. This study provides insights into the effects of long-term selection in Boer and indigenous Ugandan goat breeds, which are relevant for implementation of breeding programs and conservation of genetic resources, as well as their sustainable use and management.

Runs of homozygosity in a selected cattle population with extremely inbred bulls: Descriptive and functional analyses revealed highly variable patterns

The analysis of runs of homozygosity (ROH), using high throughput genomic data, has become a valuable and frequently used methodology to characterize the genomic and inbreeding variation of livestock and wildlife animal populations. However, this methodology has been scarcely used in highly inbred domestic animals. Here, we analyzed and characterized the occurrence of ROH fragments in highly inbred (HI; average pedigree-based inbreeding coefficient F_PED = 0.164; 0.103 to 0.306) and outbred Retinta bulls (LI; average F_PED = 0.008; 0 to 0.025). We studied the length of the fragments, their abundance, and genome distribution using high-density microarray data. The number of ROH was significantly higher in the HI group, especially for long fragments (>8Mb). In the LI group, the number of ROH continuously decreased with fragment length. Genome-wide distribution of ROH was highly variable between samples. Some chromosomes presented a larger number of fragments (BTA1, BTA19, BTA29), others had longer fragments (BTA4, BTA12, BTA17), while other ones showed an increased ROH accumulation over specific loci (BTA2, BTA7, BTA23, BTA29). Similar differences were observed in the analysis of 12 individuals produced by a similar inbred event (F_PED3 = 0.125). The correlation between the fraction of the genome covered by ROH (F_ROH) and F_PED was high (0.79), suggesting that ROH-based estimations are indicative of inbreeding levels. On the other hand, the correlation between F_PED and the microsatellite-based inbreeding coefficient (F_MIC) was only moderate (r = 0.44), suggesting that STR-based inbreeding estimations should be avoided. Similarly, we found a very low correlation (r = -0.0132) between recombination rate and ROH abundance across the genome. Finally, we performed functional annotation analyses of genome regions with significantly enriched ROH abundance. Results revealed gene clusters related to pregnancy-associated proteins and immune reaction. The same analysis performed for regions enriched with recently formed ROH (> 8 Mb) showed gene clusters related to flagellum assembly. In both cases, the processes were related to male and female reproductive functions, which may partially explain the reduced fertility associated with inbred populations.

High-density Genotyping reveals Genomic Characterization, Population Structure and Genetic Diversity of Indian Mithun (Bos frontalis)

The current study aimed at genomic characterization and improved understanding of genetic diversity of two Indian mithun populations (both farm, 48 animals and field, 24 animals) using genome wide genotype data generated with Illumina BovineHD BeadChip. Eight additional populations of taurine cattle (Holstein and NDama), indicine cattle (Gir) and other evolutionarily closely related species (Bali cattle, Yak, Bison, Gaur and wild buffalo) were also included in this analysis (N = 137) for comparative purposes. Our results show that the genetic background of mithun populations was uniform with few possible signs of indicine admixture. In general, observed and expected heterozygosities were quite similar in these two populations. We also observed increased frequencies of small-sized runs of homozygosity (ROH) in the farm population compared to field mithuns. On the other hand, longer ROH were more frequent in field mithuns, which suggests recent founder effects and subsequent genetic drift due to close breeding in farmer herds. This represents the first study providing genetic evidence about the population structure and genomic diversity of Indian mithun. The information generated will be utilized for devising suitable breeding and conservation programme for mithun, an endangered bovine species in India.

Founder-specific inbreeding depression affects racing performance in Thoroughbred horses

The Thoroughbred horse has played an important role in both sporting and economic aspects of society since the establishment of the breed in the 1700s. The extensive pedigree and phenotypic information available for the Thoroughbred horse population provides a unique opportunity to examine the effects of 300 years of selective breeding on genetic load. By analysing the relationship between inbreeding and racing performance of 135,572 individuals, we found that selective breeding has not efficiently alleviated the Australian Thoroughbred population of its genetic load. However, we found evidence for purging in the population that might have improved racing performance over time. Over 80% of inbreeding in the contemporary population is accounted for by a small number of ancestors from the foundation of the breed. Inbreeding to these ancestors has variable effects on fitness, demonstrating that an understanding of the distribution of genetic load is important in improving the phenotypic value of a population in the future. Our findings hold value not only for Thoroughbred and other domestic breeds, but also for small and endangered populations where such comprehensive information is not available.

A heuristic method to identify runs of homozygosity associated with reduced performance in livestock

Although, for the most part, genome-wide metrics are currently used in managing livestock inbreeding, genomic data offer, in principle, the ability to identify functional inbreeding. Here, we present a heuristic method to identify haplotypes contained within a run of homozygosity (ROH) associated with reduced performance. Results are presented for simulated and swine data. The algorithm comprises 3 steps. Step 1 scans the genome based on marker windows of decreasing size and identifies ROH genotypes associated with an unfavorable phenotype. Within this stage, multiple aggregation steps reduce the haplotype to the smallest possible length. In step 2, the resulting regions are formally tested for significance with the use of a linear mixed model. Lastly, step 3 removes nested windows. The effect of the unfavorable haplotypes identified and their associated haplotype probabilities for a progeny of a given mating pair or an individual can be used to generate an inbreeding load matrix (ILM). Diagonals of ILM characterize the functional individual inbreeding load (IIL). We estimated the accuracy of predicting the phenotype based on IIL. We further compared the significance of the regression coefficient for IIL on phenotypes with genome-wide inbreeding metrics. We tested the algorithm using simulated scenarios (12 scenarios), combining different levels of linkage disequilibrium (LD) and number of loci impacting a quantitative trait. Additionally, we investigated 9 traits from 2 maternal purebred swine lines. In simulated data, as the LD in the population increased, the algorithm identified a greater proportion of the true unfavorable ROH effects. For example, the proportion of highly unfavorable true ROH effects identified rose from 32 to 41% for the low- to the high-LD scenario. In both simulated and real data, the haplotypes identified were contained within a much larger ROH (9.12-12.1 Mb). The IIL prediction accuracy was greater than 0 across all scenarios for simulated data (mean of 0.49 [95% confidence interval 0.47-0.52] for the high-LD scenario) and for nearly all swine traits (mean of 0.17 [SD 0.10]). On average, across simulated and swine data sets, the IIL regression coefficient was more closely related to progeny performance than any genome-wide inbreeding metric. A heuristic method was developed that identified ROH genotypes with reduced performance and characterized the combined effects of ROH genotypes within and across individuals.

Trends in genome-wide and region-specific genetic diversity in the Dutch-Flemish Holstein-Friesian breeding program from 1986 to 2015

Background

In recent decades, Holstein–Friesian (HF) selection schemes have undergone profound changes, including the introduction of optimal contribution selection (OCS; around 2000), a major shift in breeding goal composition (around 2000) and the implementation of genomic selection (GS; around 2010). These changes are expected to have influenced genetic diversity trends. Our aim was to evaluate genome-wide and region-specific diversity in HF artificial insemination (AI) bulls in the Dutch-Flemish breeding program from 1986 to 2015.

Methods

Pedigree and genotype data (~ 75.5 k) of 6280 AI-bulls were used to estimate rates of genome-wide inbreeding and kinship and corresponding effective population sizes. Region-specific inbreeding trends were evaluated using regions of homozygosity (ROH). Changes in observed allele frequencies were compared to those expected under pure drift to identify putative regions under selection. We also investigated the direction of changes in allele frequency over time.

Results

Effective population size estimates for the 1986–2015 period ranged from 69 to 102. Two major breakpoints were observed in genome-wide inbreeding and kinship trends. Around 2000, inbreeding and kinship levels temporarily dropped. From 2010 onwards, they steeply increased, with pedigree-based, ROH-based and marker-based inbreeding rates as high as 1.8, 2.1 and 2.8% per generation, respectively. Accumulation of inbreeding varied substantially across the genome. A considerable fraction of markers showed changes in allele frequency that were greater than expected under pure drift. Putative selected regions harboured many quantitative trait loci (QTL) associated to a wide range of traits. In consecutive 5-year periods, allele frequencies changed more often in the same direction than in opposite directions, except when comparing the 1996–2000 and 2001–2005 periods.

Conclusions

Genome-wide and region-specific diversity trends reflect major changes in the Dutch-Flemish HF breeding program. Introduction of OCS and the shift in breeding goal were followed by a drop in inbreeding and kinship and a shift in the direction of changes in allele frequency. After introduction of GS, rates of inbreeding and kinship increased substantially while allele frequencies continued to change in the same direction as before GS. These results provide insight in the effect of breeding practices on genomic diversity and emphasize the need for efficient management of genetic diversity in GS schemes.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0385-y) contains supplementary material, which is available to authorized users.

Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Age-based partitioning of individual genomic inbreeding levels in Belgian Blue cattle

Background

Inbreeding coefficients can be estimated either from pedigree data or from genomic data, and with genomic data, they are either global or local (when the linkage map is used). Recently, we developed a new hidden Markov model (HMM) that estimates probabilities of homozygosity-by-descent (HBD) at each marker position and automatically partitions autozygosity in multiple age-related classes (based on the length of HBD segments). Our objectives were to: (1) characterize inbreeding with our model in an intensively selected population such as the Belgian Blue Beef (BBB) cattle breed; (2) compare the properties of the model at different marker densities; and (3) compare our model with other methods.

Results

When using 600 K single nucleotide polymorphisms (SNPs), the inbreeding coefficient (probability of sampling an HBD locus in an individual) was on average 0.303 (ranging from 0.258 to 0.375). HBD-classes associated to historical ancestors (with small segments ≤ 200 kb) accounted for 21.6% of the genome length (71.4% of the total length of the genome in HBD segments), whereas classes associated to more recent ancestors accounted for only 22.6% of the total length of the genome in HBD segments. However, these recent classes presented more individual variation than more ancient classes. Although inbreeding coefficients obtained with low SNP densities (7 and 32 K) were much lower (0.060 and 0.093), they were highly correlated with those obtained at higher density (r = 0.934 and 0.975, respectively), indicating that they captured most of the individual variation. At higher SNP density, smaller HBD segments are identified and, thus, more past generations can be explored. We observed very high correlations between our estimates and those based on homozygosity (r = 0.95) or on runs-of-homozygosity (r = 0.95). As expected, pedigree-based estimates were mainly correlated with recent HBD-classes (r = 0.56).

Conclusions

Although we observed high levels of autozygosity associated with small HBD segments in BBB cattle, recent inbreeding accounted for most of the individual variation. Recent autozygosity can be captured efficiently with low-density SNP arrays and relatively simple models (e.g., two HBD classes). The HMM framework provides local HBD probabilities that are still useful at lower SNP densities.

Electronic supplementary material

The online version of this article (10.1186/s12711-017-0370-x) contains supplementary material, which is available to authorized users.

Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep

Background

Because very large numbers of single nucleotide polymorphisms (SNPs) are now available throughout the genome, they are particularly suitable for the detection of genomic regions where a reduction in heterozygosity has occurred and they offer new opportunities to improve the accuracy of inbreeding (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F) estimates. Runs of homozygosity (ROH) are contiguous lengths of homozygous segments of the genome where the two haplotypes inherited from the parents are identical. Here, we investigated the occurrence and distribution of ROH using a medium-dense SNP panel to characterize autozygosity in 516 Valle del Belice sheep and to identify the genomic regions with high ROH frequencies.

Results

We identified 11,629 ROH and all individuals displayed at least one ROH longer than 1 Mb. The mean value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F estimated from ROH longer than1 Mb was 0.084 ± 0.061. ROH that were shorter than 10 Mb predominated. The highest and lowest coverages of Ovis aries chromosomes (OAR) by ROH were on OAR24 and OAR1, respectively. The number of ROH per chromosome length displayed a specific pattern, with higher values for the first three chromosomes. Both number of ROH and length of the genome covered by ROH varied considerably between animals. Two hundred and thirty-nine SNPs were considered as candidate markers that may be under directional selection and we identified 107 potential candidate genes. Six genomic regions located on six chromosomes, corresponding to ROH islands, are presented as hotspots of autozygosity, which frequently coincided with regions of medium recombination rate. According to the KEGG database, most of these genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. A genome scan revealed the presence of ROH islands in genomic regions that harbor candidate genes for selection in response to environmental stress and which underlie local adaptation.

Conclusions

These results suggest that natural selection has, at least partially, a role in shaping the genome of Valle del Belice sheep and that ROH in the ovine genome may help to detect genomic regions involved in the determinism of traits under selection.

Electronic supplementary material

The online version of this article (10.1186/s12711-017-0360-z) contains supplementary material, which is available to authorized users.

Population structure and genomic inbreeding in nine Swiss dairy cattle populations

BACKGROUND

Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF).

RESULTS

The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26.

CONCLUSIONS

The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.

Estimation of inbreeding depression on female fertility in the Finnish Ayrshire population

Single nucleotide polymorphism (SNP) data enable the estimation of inbreeding at the genome level. In this study, we estimated inbreeding levels for 19,075 Finnish Ayrshire cows genotyped with a low-density SNP panel (8K). The genotypes were imputed to 50K density, and after quality control, 39,144 SNPs remained for the analysis. Inbreeding coefficients were estimated for each animal based on the percentage of homozygous SNPs (F_PH ), runs of homozygosity (F_ROH ) and pedigree (F_PED ). Phenotypic records were available for 13,712 animals including non-return rate (NRR), number of inseminations (AIS) and interval from first to last insemination (IFL) for heifers and up to three parities for cows, as well as interval from calving to first insemination (ICF) for cows. Average F_PED was 0.02, F_ROH 0.06 and F_PH 0.63. A correlation of 0.71 was found between F_PED and F_ROH , 0.66 between F_PED and F_PH and 0.94 between F_ROH and F_PH . Pedigree-based inbreeding coefficients did not show inbreeding depression in any of the traits. However, when F_ROH or F_PH was used as a covariate, significant inbreeding depression was observed; a 10% increase in F_ROH was associated with 5 days longer IFL0 and IFL1, 2 weeks longer IFL3 and 3 days longer ICF2 compared to non-inbred cows.

Genome-wide mapping and estimation of inbreeding depression of semen quality traits in a cattle population

Inbreeding depression is known to affect quantitative traits such as male fertility and sperm quality, but the genetic basis for these associations is poorly understood. Most studies have been limited to examining how pedigree- or marker-derived genome-wide autozygosity is associated with quantitative phenotypes. In this study, we analyzed possible associations of genetic features of inbreeding depression with percentage of live spermatozoa and total number of spermatozoa in 19,720 ejaculates obtained from 554 Austrian Fleckvieh bulls during routine artificial insemination programs. Genome-wide inbreeding depression was estimated and genomic regions contributing to inbreeding depression were mapped. Inbreeding depression did affect total number of spermatozoa, and such depression was predicted by pedigree-based inbreeding levels and genome-wide inbreeding levels based on runs of homozygosity (ROH). Genome-wide inbreeding depression did not seem to affect percentage of live spermatozoa. A model incorporating genetic effects of the bull, environmental factors, and additive genetic and ROH status effects of individual single-nucleotide polymorphisms revealed genomic regions significantly associated with ROH status for total number of spermatozoa (4 regions) or percentage of live spermatozoa (5 regions). All but one region contains genes related to spermatogenesis and sperm morphology. These genomic regions contain genes affecting sperm morphogenesis and efficacy. The results highlight that next-generation sequencing may help explain some of the genetic factors contributing to inbreeding depression of sperm quality traits in Fleckvieh bulls.

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH ) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED ). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome-wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.

Characterization and management of long runs of homozygosity in parental nucleus lines and their associated crossbred progeny

BACKGROUND

In nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome.

METHODS

We calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated.

RESULTS

We detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (-0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions.

CONCLUSIONS

We show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions.

Novel Graphical Analyses of Runs of Homozygosity among Species and Livestock Breeds

Runs of homozygosity (ROH), uninterrupted stretches of homozygous genotypes resulting from parents transmitting identical haplotypes to their offspring, have emerged as informative genome-wide estimates of autozygosity (inbreeding). We used genomic profiles based on 698 K single nucleotide polymorphisms (SNPs) from nine breeds of domestic cattle (Bos taurus) and the European bison (Bison bonasus) to investigate how ROH distributions can be compared within and among species. We focused on two length classes: 0.5-15 Mb to investigate ancient events and >15 Mb to address recent events (approximately three generations). For each length class, we chose a few chromosomes with a high number of ROH, calculated the percentage of times a SNP appeared in a ROH, and plotted the results. We selected areas with distinct patterns including regions where (1) all groups revealed an increase or decrease of ROH, (2) bison differed from cattle, (3) one cattle breed or groups of breeds differed (e.g., dairy versus meat cattle). Examination of these regions in the cattle genome showed genes potentially important for natural and human-induced selection, concerning, for example, meat and milk quality, metabolism, growth, and immune function. The comparative methodology presented here permits visual identification of regions of interest for selection, breeding programs, and conservation.