Aarhus University Seal / Aarhus Universitets segl

Estimation of inbreeding using pedigree, 50k SNP chip genotypes and full sequence data in three cattle breeds

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review



  • Qianqian Zhang, Denmark
  • Mario P L Calus, Animal Breeding and Genomics Centre, Wageningen UR Livestock Research, Netherlands
  • Bernt Guldbrandtsen, Denmark
  • Mogens Sandø Lund
  • Goutam Sahana
Background: Levels of inbreeding in cattle populations have increased in the past due to the use of a limited number of bulls for artificial insemination. High levels of inbreeding lead to reduced genetic diversity and inbreeding
depression. Various estimators based on different sources, e.g., pedigree or genomic data, have been used to estimate inbreeding coefficients in cattle populations. However, the comparative advantage of using full sequence data to assess inbreeding is unknown. We used pedigree and genomic data at different densities from 50k to full sequence variants to compare how different methods performed for the estimation of inbreeding levels in three different cattle breeds.

Results: Five different estimates for inbreeding were calculated and compared in this study: pedigree based inbreeding coefficient (FPED); run of homozygosity (ROH)-based inbreeding coefficients (FROH); genomic relationship matrix (GRM)-based inbreeding coefficients (FGRM); inbreeding coefficients based on excess of homozygosity (FHOM) and correlation of uniting gametes (FUNI). Estimates using ROH provided the direct estimated levels of autozygosity in the current populations and are free effects of allele frequencies and incomplete pedigrees which may increase in inaccuracy in estimation of inbreeding. The highest correlations were observed between FROH estimated from the full sequence variants and the FROH estimated from 50k SNP (single nucleotide polymorphism) genotypes. The estimator based on the correlation between uniting gametes (FUNI) using full genome sequences was also
strongly correlated with FROH detected from sequence data.

Conclusions: Estimates based on ROH directly reflected levels of homozygosity and were not influenced by allele
frequencies, unlike the three other estimates evaluated (FGRM, FHOM and FUNI), which depended on estimated allele
frequencies. FPED suffered from limited pedigree depth. Marker density affects ROH estimation. Detecting ROH based
on 50k chip data was observed to give estimates similar to ROH from sequence data. In the absence of full sequence
data ROH based on 50k can be used to access homozygosity levels in individuals. However, genotypes denser than 50k
are required to accurately detect short ROH that are most likely identical by descent (IBD).
Original languageEnglish
JournalB M C Genetics
Pages (from-to)1-16
Number of pages16
Publication statusPublished - 22 Jul 2015

See relations at Aarhus University Citationformats

Download statistics

No data available

ID: 87141857