Pedigree and genome-based patterns of homozygosity in the South African Ayrshire, Holstein, and Jersey breeds

Abstract

The erosion of genetic diversity limits long-term genetic gain and impedes the sustainability of livestock production. In the South African (SA) dairy industry, the major commercial dairy breeds have been applying estimated breeding values (EBVs) and/or have been participating in Multiple Across Country Evaluations (MACE). The transition to genomic estimated breeding values (GEBVs) in selection strategies requires monitoring of the genetic diversity and inbreeding of current genotyped animals, especially considering the comparatively small population sizes of global dairy breeds in SA. This study aimed to perform a homozygosity-based evaluation of the SA Ayrshire (AYR), Holstein (HST), and Jersey (JER) dairy cattle breeds. Three sources of information, namely 1) single nucleotide polymorphism (SNP) genotypes (3,199 animals genotyped for 35,572 SNPs) 2) pedigree records (7,885 AYR; 28,391 HST; 18,755 JER), and 3) identified runs of homozygosity (ROH) segments were used to quantify inbreeding related parameters. The lowest pedigree completeness was for the HST population reducing from a value of 0.990 to 0.186 for generation depths of one to six. Across all breeds, 46.7% of the detected ROH were between 4 megabase pairs (Mb) and 8 Mb in length. Two conserved homozygous haplotypes were identified in more than 70% of the JER population on Bos taurus autosome (BTA) 7. The JER breed displayed the highest level of inbreeding across all inbreeding coefficients. The mean (± standard deviation) pedigree-based inbreeding coefficient (FPED) ranged from 0.051 (±0.020) for AYR to 0.062 (±0.027) for JER, whereas SNP-based inbreeding coefficients (FSNP) ranged from 0.020 (HST) to 0.190 (JER) and ROH-based inbreeding coefficients, considering all ROH segment coverage (FROH), ranged from 0.053 (AYR) to 0.085 (JER). Within-breed Spearman correlations between pedigree-based and genome-based estimates ranged from weak (AYR: 0.132 between FPED and FROH calculated for ROH <4Mb in size) to moderate (HST: 0.584 between FPED and FSNP). Correlations strengthened between FPED and FROH as the ROH length category was considered lengthened, suggesting a dependency on breed-specific pedigree depth. The genomic homozygosity-based parameters studied proved useful in investigating the current inbreeding status of reference populations genotyped to implement genomic selection in the three most prominent South African dairy cattle breeds.

Genomic analysis revealed hotspots of genetic adaptation and risk of disappearance in the Brazilian goat populations

“…In general, variations in the number and distribution of ROH may reflect historical differences in genomic structures. Highly selected breeds often show a greater abundance of ROH and broader coverage compared to local breeds [15,45–47]. In our study, the Marota landrace is distinguished by its high number of ROH, contrasting with the Anglo-Nubian commercial breed.…”

Genetic diversity and runs of homozygosity (ROH): A portrait of the quantitative academic publication dynamic and scientific metadata

“…[44,46]. The R package detectRUNS has seen a signi cant surge in its utilization and is likely one of the most widely adopted programs for analyzing SNP genotypes and runs of homozygosity (ROH) and Heterozygosity-Rich Regions (HRR), across both animal populations [47][48][49][50]…”