BACKGROUND/OBJECTIVES: Whole-genome sequencing (WGS) data provide valuable information about the genetic architecture of local livestock but have not yet been applied to Russian native goats, in particular, the Orenburg and Karachay breeds. A preliminary search for selection signatures based on single nucleotide polymorphism (SNP) genotype data in these breeds was not informative. Therefore, in this study, we aimed to address runs of homozygosity (ROHs) patterns and find the respective signatures of selection overlapping candidate genes in Orenburg and Karachay goats using the WGS approach.

METHODS

Paired-end libraries (150 bp reads) were constructed for each animal. Next-generation sequencing was performed using a NovaSeq 6000 sequencer (Illumina, Inc., San Diego, CA, USA), with ~20X genome coverage. ROHs were identified in sliding windows, and ROH segments shared by at least 50% of the samples were considered as ROH islands.

RESULTS

ROH islands were identified on chromosomes CHI3, CHI5, CHI7, CHI12, CHI13, and CHI15 in Karachay goats; and CHI3, CHI11, CHI12, CHI15, and CHI16 in Orenburg goats. Shared ROH islands were found on CHI12 (containing the and candidate genes) and on CHI15 (harboring and ). The Karachay breed had greater ROH length and higher ROH number compared to the Orenburg breed (134.13 Mb and 695 vs. 78.43 Mb and 438, respectively). The genomic inbreeding coefficient () varied from 0.032 in the Orenburg breed to 0.054 in the Karachay breed. Candidate genes associated with reproduction, milk production, immunity-related traits, embryogenesis, growth, and development were identified in ROH islands in the studied breeds.

CONCLUSIONS

Here, we present the first attempt of elucidating the ROH landscape and signatures of selection in Russian local goat breeds using WGS analysis. Our findings will pave the way for further insights into the genetic mechanisms underlying adaption and economically important traits in native goats.