Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2025 May 24;16(6):631.
doi: 10.3390/genes16060631.

Examination of Runs of Homozygosity Distribution Patterns and Relevant Candidate Genes of Potential Economic Interest in Russian Goat Breeds Using Whole-Genome Sequencing

Affiliations

Examination of Runs of Homozygosity Distribution Patterns and Relevant Candidate Genes of Potential Economic Interest in Russian Goat Breeds Using Whole-Genome Sequencing

Tatiana E Deniskova et al. Genes (Basel). .

Abstract

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 PARP4 and MPHOSPH8 candidate genes) and on CHI15 (harboring STIM1 and RRM1). 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 (FROH) 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.

Keywords: Russian local breeds; candidate genes; goat (Capra hircus); runs of homozygosity (ROHs); signatures of selection; whole-genome sequencing (WGS).

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflicts of interest. The RSF had no role in the design of this study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Figures

Figure 1
Figure 1
The representative studied goat breeds of Russian origin: (a) a family of Orenburg goats (a black male and two white females); and (b) a herd of Karachay goats.
Figure 2
Figure 2
The map with the corresponding sampling sites of goats for the Orenburg (▼) and Karachay (▼) breeds.
Figure 3
Figure 3
Individual Neighbor-Net graph using IBS distances plotted for the Orenburg and Karachay breeds based on the WGS data analysis.
Figure 4
Figure 4
Principal component analysis performed for Orenburg and Karachay goats based on the WGS data analysis: (a) for the first two principal components (PC1 and PC2); and (b) for the first and third principal components (PC1 and PC3).
Figure 5
Figure 5
Distribution of runs of homozygosity (ROH) according to length classes in Orenburg and Karachay goats based on the WGS data analysis: (a) mean number of ROHs per breed; and (b) mean length of ROHs per breed.
Figure 6
Figure 6
ROH islands identified in the Orenburg and Karachay breeds on the following chromosomes: (a) CHI3, (b) CHI5, (c) CHI7, (d) CHI11, (e) CHI12, (f) CHI13, (g) CHI15, and (h) CHI16.

Similar articles

References

    1. Signer-Hasler H., Henkel J., Bangerter E., Bulut Z., Drögemüller C., Leeb T., Flury C. Runs of homozygosity in Swiss goats reveal genetic changes associated with domestication and modern selection. Genet. Sel. Evol. 2022;54:6. doi: 10.1186/s12711-022-00695-w. - DOI - PMC - PubMed
    1. Saif R., Henkel J., Mahmood T., Ejaz A., Ahmad F., Zia S. Detection of whole genome selection signatures of Pakistani Teddy goat. Mol. Biol. Rep. 2021;48:7273–7280. doi: 10.1007/s11033-021-06726-x. - DOI - PubMed
    1. Dan H., Zhong H., Akhatayeva Z., Lin K., Xu S. Whole-genome selective scans detect genes associated with cashmere traits and climatic adaptation in cashmere goats (Capra hircus) in China. Genes. 2025;16:292. doi: 10.3390/genes16030292. - DOI - PMC - PubMed
    1. Gudra D., Valdovska A., Jonkus D., Kairisa D., Galina D., Ustinova M., Viksne K., Fridmanis D., Kalnina I. Genetic characterization of the Latvian local goat breed and genetic traits associated with somatic cell count. Animal. 2024;18:101154. doi: 10.1016/j.animal.2024.101154. - DOI - PubMed
    1. Liu J., Dong S., Lv J., Li Y., Sun B., Guo Y., Deng M., Liu D., Liu G. Screening of SNP loci related to leg length trait in Leizhou goats based on whole-genome resequencing. Int. J. Mol. Sci. 2024;25:12450. doi: 10.3390/ijms252212450. - DOI - PMC - PubMed

LinkOut - more resources