Next Article in Journal
Seroprevalence of Swine Influenza A Virus (swIAV) Infections in Commercial Farrow-to-Finish Pig Farms in Greece
Next Article in Special Issue
Influence of Sex and a High-Fiber Diet on the Gut Microbiome of Alentejano Pigs Raised to Heavy Weights
Previous Article in Journal
Novel Prion Protein Gene Polymorphisms in Awassi Sheep in Three Regions of the Fertile Crescent
Previous Article in Special Issue
Cloning of the RNA m6A Methyltransferase 3 and Its Impact on the Proliferation and Differentiation of Quail Myoblasts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity of Kazakhstani Equus caballus (Linnaeus, 1758) Horse Breeds Inferred from Microsatellite Markers

by
Zarina Orazymbetova
1,2,†,
Daniya Ualiyeva
2,3,4,*,
Kairat Dossybayev
1,2,†,
Aibyn Torekhanov
1,
Dauren Sydykov
1,
Aizhan Mussayeva
2 and
Gabiden Baktybayev
1
1
Kazakh Research Institute of Livestock and Fodder Production, Almaty 050035, Kazakhstan
2
Institute of Genetics and Physiology, Committee of Science of the Ministry of Science and Higher Education of the Republic of Kazakhstan, Almaty 050060, Kazakhstan
3
Institute of Zoology, Committee of Science of the Ministry of Science and Higher Education of the Republic of Kazakhstan, Almaty 050060, Kazakhstan
4
Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu 610041, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Vet. Sci. 2023, 10(10), 598; https://doi.org/10.3390/vetsci10100598
Submission received: 25 August 2023 / Revised: 21 September 2023 / Accepted: 28 September 2023 / Published: 30 September 2023
(This article belongs to the Special Issue Genetic Analysis in Domestic Animal Breeds and Performance)

Abstract

:

Simple Summary

Traditional horse breeding has long developed several kinds and lineages of horse breeds in Kazakhstan. Among them, Kushum and Mugalzhar are the breeds most prolific and resistant to harsh climatic conditions. Microsatellite analysis is employed to examine present genetic variability and population structure. The subpopulation structure shows three regional groups indicating (I) purebred Kushum populations from western Kazakhstan, (II) the Kozhamberdy type of Mugalzhar interbreed in populations from central Kazakhstan, and (III) admixed Kushum–Mugalzhar populations from western and southwestern Kazakhstan. The majority of microsatellite markers utilized are informative, with seven of them being extremely variable. A high level of genetic admixture among Kazakhstani Equus caballus breeds is found, as well as a shallow level of genetic differentiation between the examined populations.

Abstract

Understanding the genetic diversity and structure of domesticated horse (Equus caballus) populations is critical for long-term herd management and breeding programs. This study examines 435 horses from Kazakhstan, covering seven groups in three geographic areas using 11 STR markers. Identified are 136 alleles, with the mean number of alleles per locus ranging from 9 to 19. VHL20 is the most variable locus across groups, while loci HTG4, AHT4, AHT5, HTG7, and HMS3 are variable in most populations. The locus AHT5 in the Emba population shows the highest frequency of rare alleles, while the lowest frequency, 0.005, is observed in the Kulandy population. All loci were highly informative for the Kazakhstani populations of E. caballus, with PIC values higher than 0.5. Pairwise variations in Wright’s FST distances show that the examined varieties have little genetic differentiation (0.05%), indicating a high degree of admixture and a continuing lineage sorting process. Phylogenetic and population structure analyses reveal three major clusters of Kazakh horses, representing (I) the Uralsk population of the Kushum breed and the monophyly of two groups: (II) the Kozhamberdy population of the Mugalzhar breed, and (III) the Mugalzhar–Kushum breed populations. Kazakhstani horse populations, while being regionally isolated, were recently in contact with each other.

Graphical Abstract

1. Introduction

The horse was the last domesticated animal which had a long-term significant and enduring impact on human civilizations, leading to advancements in transportation and trade, influencing people’s lifestyle, and revolutionizing the nature of warfare. In turn, human selection changed genetic diversity in horse populations, resulting in the variation seen among current horse characteristics and types. Domestic horses were subject to constant genetic restocking from the wild, mostly from females [1,2,3,4].
Kazakhstan, the largest Central Asian country, combines several climatic and geomorphologic features that offer an optimal environment for agricultural growth due to its geographical location at the intersection of the western and eastern Palearctic. The FAO database (Domestic Animal Diversity Information System, accessed on 5 December 2022), identifies Kazakhstan as having eight registered domesticated horse breeds (Kushum, Mugalzhar, Kostanay, Akhal-Teke, Aday, Kazakh, Jabe, and Karabair) [5]. Among them the Kushum breed was formed via purebred crossbreeding, whereas the Mugalzhar breed was established by traditional crossing methods, and these two breeds are dedicated to meat and milk production. The dairy production presented by fermented mare’s milk called “Kymyz” and fresh milk “Saumal” are national drink products of Kazakhstan. They are rich in minerals and vitamins, and benefit human health with anticancer and anti-inflammation effects [6]. The overall number of horses in Kazakhstan now surpasses 3 million, with 90% of them raised by herds in practically all areas of the country [7].
The Kushum horse was bred by the herd method of Ural and Aktobe farms, named after the Kushum River, flowing across western Kazakhstan. Initial approval as a separate breed came in 1976. The breed was established by a complex reproductive crossover procedure. Trotting, Thoroughbred riding, and Don breeds were crossed with local native mares [8]. The farmers chose the best animals in terms of growth and maintained them on semi-desert cereal–wormwood pastures throughout the year. As a result, Kushum horses are ideally suited for herding. Crosses and rigorous selection of stallions for desired traits were carried out in the middle of the twentieth century in terms of external features and adaptation to local environments. Later, the crossed horses were bred “in themselves” to reinforce breeding results [9].
Mugalzhar horses are an upgraded dual-purpose form of the Kazakh horse breed Jabe type, developed by Kazakhstani scientists between 1969 and 1998 [10]. They are distinguished by high-quality meat and milk production, as well as resilience to extreme weather conditions and year-round grassland farming [11].
Many studies have been conducted on a global scale in the horse breeding sector involving research of mtDNA analysis of modern horses [12,13,14,15,16] and fossils [17,18]. Microsatellite analysis has been widely utilized for parentage testing [19,20,21], as well as a single nucleotide polymorphism (SNP) analysis to determine the population structure and genetic variability of Equids [22,23,24].
Genetic distinctiveness can be clarified by population genetic structure, leading to breed preservation, including future breeding methods and management plans. Microsatellites have been a widely utilized genetic marker and effectively applied to studies of inter- and intrabreed variability in domestic and feral horse populations [25,26,27,28,29,30,31,32].
Here we determine the genetic diversity of the Kushum and Mugalzhar horse breeds within populations and between them. Applied are several different approaches to evaluate the distribution of molecular indices, genetic distances, and population structural variation, which can elucidate lineage sorting processes, and be helpful for organizing individual genetic feature panels of the type or line as well as to improve the inbreeding management strategy of Kazakh horse-breed populations.

2. Materials and Methods

2.1. Sampling

DNA isolation of 435 samples from blood cohorts were taken from 7 populations across Kazakhstan, representing Mugalzhar and Kushum horse breeds (Table 1; Figure 1).

2.2. Laboratory Protocol

A panel of 11 highly polymorphic microsatellite loci (AHT4, AHT5, ASB2, ASB23, HMS3, HMS6, HMS7, HTG4, HTG7, LEX3, and VHL20) were chosen for typing using methods described by Zabek and Fornal (2009) [33].
DNA extractions were conducted using the “DNA-sorb-B” set (AmpliSens, Moscow, Russia). Biomaterial was processed by the solid-phase sorption method, which consists in adding a lysing solution, DNA sorption on a sorbent, repeated washing and resorption of DNA with a buffer solution, as a result of which a purified solution containing DNA was obtained. Furthermore, the spectrometric quantification was performed as well as amplification reactions carried out using the StockMarks Equine Kit (Applied Biosystems, Waltham, MA, USA) [34]. Separation and analysis of amplified fragments were carried out by capillary electrophoresis using a genetic analyzer. PCR amplifications were performed on Thermocycler 2730 (Applied Biosystems) following a touchdown cycling protocol with an initial denaturation at 95 °C for 15 min, followed by 30 cycles of: the first 4 cycles, 58 °C (30 s), 59 °C (120 s), 72 °C (75 s); the next 6 cycles, 94 °C (30 s), 59 °C (120 s), 72 °C (75 s); the next 20 cycles, 90 °C (30 s), 59 °C (120 s), 72 °C (75 s); ending in a 5 min extension at 68 °C, with a 4 °C hold temperature. Amplification product separation was performed by capillary electrophoresis on an automatic genetic analyzer AB 3130 (Applied Biosystems), using the GeneMapper™ v. 4.0 program. Amplified DNA fragments were interpreted using a control DNA profile with a known genotype and data from international comparative tests (Horse Comparison Tests) conducted by ISAG.

2.3. Population Genetic Structure

Allele frequencies and polymorphic information content (PIC) were calculated using Cervus 3.0 software [35,36]. Genetic diversity within and between breeds, as well as basic parameters, including total number of allele variants (NA), effective number of alleles (NE), estimation of observed (HO), expected (HE), and unbiased expected (UHE) heterozygosity, and Shannon’s information index (I) were measured using GenAlEx 6.5 software (New Brunswick, NJ, USA) [37]. Variance components of microsatellite diversity within and between populations for all pairs of populations were analyzed using analysis of molecular variance (AMOVA) with permutations set to 999 in the GenAlEx 6.5 [37]. Chi-square tests of Hardy–Weinberg equilibrium and rare alleles were calculated for each population using Microsatellite Analyzer v. 4.05 (MSA) [38]. Fixation indices (FIT, FIS, and FST) of Wright’s F-statistics were obtained using GenAlex 6.5 and Excel microsatellite toolkit (version 3.1) [37]. Neighbor joining of Saitou and Nei (1987) [39] was used to construct a phylogenetic tree based on Nei’s genetic distance in MEGA 7 [40]. Factorial correspondence analysis (FCA) was investigated based on the individual multilocus genotype using GENETIX version 4.03 [41]. Bayesian clustering analysis was implemented in Structure 2.3.4 [42] without prior structure information. All possibilities were considered by dividing 7 populations into 7 groups. An ad hoc quantity based on the second order rate of change in the likelihood function with respect to K (K) was used for estimating the number of clusters from structure analysis [43]. In addition, we also use ln(Pr(X|K) values in order to identify the k for which Pr(K = k) is highest (as described in STRUCTURE’s manual, Section 5.1. Twenty runs for K = 1 to 7 were analyzed under the admixture model, correlated allele frequencies, and a burn-in of 250,000 followed by 1,000,000 Markov chain Monte Carlo (MCMC) iterations. Structure Harvester 0.6.93 [44] was applied to choose the optimal K-value based on the Delta K method. The 20 replicates for the chosen K-value were merged using CLUMPP 1.1.2 [45] and the final plots were generated using DISTRUCT 1.1 [46].

3. Results

Microsatellite Genotyping, and Population Genetic Diversity and Structure

A total of 136 alleles at 11 STR loci from 435 genotyped individuals of two Kazakh horse breeds from seven populations were identified. All markers were found to be polymorphic (p ≥ 0.05) (Table 2). The mean number of alleles varied from 9 at loci AHT4, HMS6, and HMS7 to 19 at locus ASB23. The mean number of alleles (Na) per locus was 12.36, and the effective number of alleles (Ne) was 5.82. The expected heterozygosity (He), which is a widely accepted measure of genetic diversity in a population, ranged from 0.64 in locus HTG4 to 0.83 in locus VHL20, with an average He of 0.77 across the seven populations for the 11 microsatellite loci analyzed. The observed heterozygosity (Ho) fluctuated from 0.48 in locus LEX3 to 0.85 in locus VHL20, with a population mean of 0.68, indicating that all studied lineages are characterized by considerable genetic variability. The polymorphic information content (PIC) varied from 0.62 for the marker ASB2 to 0.82 for the AHT4 locus. The average PIC for the 11 microsatellite markers was 0.74 and there were no markers with a PIC of less than 0.5, indicating that all loci were found to be highly polymorphic. Shannon’s information (diversity) index (I), which is an indicator of the genetic variability of a population, ranged from 1.33 in locus HTG4 to 1.99 in the VHL20 marker. The average value of the I-index for all seven populations was equal to 1.73, which reflects the level of allele abundancy (Table 2). Further, FIS, FIT, and FST indices were calculated for each marker in whole populations. FIS ranged from −0.031 (VHL20) to 0.093 (AHT4) with an average value of 0.211 for all loci. FIT presented a mean value of 0.157 ranging from −0.010 for HTG7 to 0.439 for LEX3. The calculation of FIS was between 0.030 (AHT5) and 0.074 (ASB2) with a mean value of 0.041 in the total population.
AMOVA analysis performed on seven populations, suggests that the majority of the variation occurred within individuals—70% (Table 3). Fixation indices based on standard permutation demonstrated differences (p ≥ 0.001) indicating a reduction of heterozygosity, panmixia, and inbreeding processes which occurred in Kazakhstani populations of E. caballus.
Across Kazakhstan, horse breeds are identified with rare alleles that are typical for each population: VHL20, HTG4, HMS3, HMS6, HMS7, AHT4, AHT5, ASB2, ASB23, and LEX3. Among them, five unique alleles were observed in the Uralsk (Population 1) and Kulandy (Population 7). For the Uralsk population unique alleles were found in the HTG locus at 121 bp and 125 bp lengths, at HMS6 at 155 bp length, and HMS7 169 bp length. For the Kulandy population a unique allele was found in the locus LEX3 143 bp length.
The value for gene differentiation based on F-statistic (FST) distance over all loci between populations of the Kushum breed (Uralsk with Aktobe) was 4.5%, whereas between populations of the Mugalzhar breed it varied in the range of 0.008–2.8%. Genetic variability between the Uralsk population and Mugalzhar populations was from 3.3% to 4.7%, which indicates that 4.7% of the variability could be attributed to differences between breeds (Table 4). A chi-squared test observed statistically significant (p ≤ 0.05) results at all loci, rejecting the null hypothesis of random mating (Table S1).
Factorial correspondence analysis revealed three clusters of horse populations are distinct at three axes with variance of 35.98%, 25.27%, and 17.40%, respectively (Figure 2). An FCA plot demonstrated that the Uralsk population of the Kushum breed was clearly separated from other horses and thus the result is consistent with the phylogenetic tree and structure inferences.
An unrooted neighbor-joining tree for all samples was constructed using a pairwise population matrix of Nei’s genetic distances in order to represent relationships among seven populations of Kazakh horse breeds (Figure 3). Three main groups were recovered: Group I is the Uralsk population with a distance of 2.8% to monophyletic groups II and III. Group II is the Kozhamberdy type population representing three lines: IIa, Meiman; IIb, Maupas; and IIc, Mesker (1.4%). Group III clusters members of two breeds consisting of IIIa, Aktobe (Kushum breed); IIIb, Emba (Mugalzhar breed); and IIIc, Kulandy type populations (Mugalzhar breed) (1.6–1.7%).
Bayesian cluster analysis performed with STRUCTURE [33] showed that the independent runs from K = 2 to K = 7 produced consistent results, where the most likely K values were identified at K = 3 and 4 (ΔK = 12.953; 12.935), respectively (Figure 4); the subpopulations’ structure [32] using the median values of Ln Prob of data to calculate Prob(K = k) yielded the uppermost value of K = 7. A plot with the clustering of individuals is presented in Figure S1.

4. Discussion

A comprehensive genetic analysis of microsatellite markers conducted for seven populations of two main horse breeds of Kazakhstan revealed a high genetic diversity. Wright’s FST analysis revealed a genetic differentiation of 4.5% between subpopulations of the Kushum breed, while to the Kulandy population of Mugalzhar breed the genetic distance was 4.7%. Nonetheless, the level of genetic variability between the Uralsk and Aktobe populations was determined to be weak, with p ≤ 0.05. This might be explained by the attribution to the same breed and possible stallion interchange given their geographic proximity.
The mean number of alleles in this study was 9 at locus AHT4 which is slightly lower than in Danubian horses (11.2) [47]. The same locus was characterized by high PIC values (0.82), which is in concordance with Halla horses [48]. The average of expected heterozygosity across the loci was computed at 0.77 in Kazakh horses and Mugalzhar horses [27], and 0.8 in Mongolian horses, which was higher than in Thoroughbred horses (0.72) [48]. The minimum observed heterozygosity was identified at locus LEX3 (0.479), where notably the same finding was reported in the study of [27].
The subpopulation inbreeding coefficient relative to the whole population showed that genetic differentiation is higher between breeds, compared to subpopulations [49]. Comparable results were found in Algerian horse breeds, which are not isolated breeds [50]. In the present study, FST findings were considerably higher than those in Bhutan horse breeds (0.003–0.008) [51], Turkish Rahvan (0.0019) [52], and Guizhou and Luoping (0.001) horse breeds [53]. Phylogenetic relationships recovered three main groupings representing (I) the Uralsk population, (II) the Mugalzhar populations, and (III) admixed populations of Kushum and Mugalzhar. Purebred Mugalzhar populations (II) of Meiman (IIa), Maupas (IIb), and Mesker (IIc) have a demonstrated relatedness of 1.6–1.7%. In contrast, a monophyletic group of admixed (IIIa) Kushum and (IIIb, c) Mugalzhar breed populations has a genetic distance of 2.8%. The Emba (IIIb) and Kulandy (IIIc) populations together group with the Aktobe (IIIa) population of the Kushum breed and have a difference in similarity of 1.3–1.4%.
Cluster analysis results [42] reveal the substructure of populations based on Δ K variation, which suggest the existence of three groups, also supported by phylogenetic analysis and FCA as described in Section 3. However, another approach of using the median values of Ln Prob is to calculate Prob(K = k) [43] detected for the possible number of subpopulations equal to seven. This may be explained when alpha is close to zero, which means most individuals are essentially from one population or another; however, when alpha ≥ 1, it means that most individuals are admixed [54]. On the other hand, STRUCTURE analysis indicates differences within subpopulations of studied horses. At K = 4 and K = 7 it clearly represents the difference between Uralsk and Aktobe populations despite being from the same breed. This is not a real surprise. It might be the result of the short history of the breed and the different initial mare populations. These findings are consistent with the results of phylogenetic analysis.
Further, Kulandy-type horses (IIIc) revealed a different genetic background in all K, while almost the same ancestral component was distinguished in the three lines of Kozhamberdy horses (II). At all values of K, a more similar pattern was observed between Aktobe (IIIa) and Emba (IIIb) populations. In the present study, the horse breeds, types and lines are also phenotypically diverse, sampled from different locations and various pasture areas. For instance, depending on the natural, climatic, and forage conditions of breeding, the Mugalzhar breed is divided into three types: Kozhamberdy, Emba, and Kulandy. Kozhamberdy forms the type of horses established in the Karagandy region of central Kazakhstan, with Emba from the Aktobe region of western Kazakhstan, and the Kulandy type of horses produced in the Kyzylorda region of the southwest part of the country. Moreover, this type of horse has high meat and milk productivity, and is well adapted to harsh environmental conditions, showing characteristic features of main breeding traits and steadfastly passing them on to their offspring.
Taking into account that the biological unit in domesticated animals is usually breed, differences of genetic traits among individuals characterize the endurance and adaptability of an organism. The STR loci are a universal tool for the examination of genetic variability among horse breeds [33,55,56]. Further, to examine genetic variability within and between breeds, the observed and expected heterozygosity were estimated. Results indicate a higher degree of related mating of individuals in a subpopulation and are inconsistent with previous studies [27]. To the best of our knowledge, herd farmers, to improve herds, usually use the most productive stallions on a farm without control of their origin. This may cause a low level of genetic variation in studied populations.
In the present study among identified rare alleles, private alleles were detected in loci HTG4, HMS6, and HMS7 for the Uralsk (Population 1), and in the locus LEX3 for the Kulandy (Population 7), which is highlighting their favor in genetic differentiation from each other [57,58,59]. Highly variable alleles for all populations were found in the locus VHL20, as well as the HTG4 locus for Populations 1–5 and 7; AHT4 for Populations 3 and 7; AHT5 for Populations 1–2 and 5–7; HTG7 for Populations 1–3 and 5–7; and locus HMS3 for Populations 3–4. This is consistent with Zaitcev et al. (2021) [59], where loci VHL20 and HTG4 were highly variable for the Priobskaya Russian horse breed population. The mean of observed heterozygosity was (0.68) similar to identified values of Lipizzaner (0.68) and Thoroughbred (0.68) horses. However, the mean of expected heterozygosity (0.77) was higher than those of Lipizzaner (0.69) and Thoroughbred (0.69) [60,61,62]. As a result, the calculated observed and expected heterozygosity showed that the studied populations are under inbreeding. Moreover, these results were confirmed by the F-statistics findings that FIT, FIS, and FST reflected a reduction of heterozygosity and weak divergency (FST ≤ 0.05) in the populations [63]. Confirmation is acquired here that the STR loci are appropriate for assessing genetic diversity and variation within and among Central Asian horse populations.
Considering the fact that the SNP (single nucleotide polymorphisms) genotyping approach has been widely applied to the investigations of Equines’ population structure [64,65], yet has a low mutation rate and abundance in genome compared to the microsatellite markers, we recommend continuing the study in the perspective of genome-wide SNP analysis to fulfill the whole picture of the population structure and genetic diversity of the Kazakh horses.

5. Conclusions

In the present study, microsatellite analysis allowed us to identify a division of the Uralsk population from populations of Aktobe (Kushum breed) and Kulandy (Mugalzhar breed), the latter of each differentiated by 4.5% and 4.7%, respectively. Moreover, a significant degree of population admixture between breeds and within types and lines of Kazakhstani Equus caballus was detected. STR loci were shown to be particularly helpful for characterizing the allele pool and assessing the extent of genetic differentiation of studied populations. Assessment of genetic variation and F-statistics presented a reduction of heterozygosity in all populations. Conducting a comparative microsatellite typing is critical for monitoring genetic variability and tracing patterns in the selection process. Furthermore, data on genetic diversity computed within and across breed variants might give insight into future conservation initiatives and better horse breeding management.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/vetsci10100598/s1, Figure S1: Estimation of the population structure with different K values (K = 2 to K = 7). Table S1: Chi-squared test of Hardy–Weinberg equilibrium for the 11 STR loci.

Author Contributions

Conceptualization, Z.O. and K.D.; methodology, K.D. and D.U.; software, K.D.; validation, D.U.; formal analysis, D.S., A.M. and G.B.; investigation, Z.O.; resources, Z.O., D.S., G.B. and A.T.; data curation, K.D. and D.U.; writing—original draft preparation, D.U.; writing—review and editing, K.D., A.T., D.S. and G.B.; visualization, K.D. and D.U.; supervision, Z.O.; project administration, D.S., A.M. and G.B.; funding acquisition, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Agriculture of the Republic of Kazakhstan, grant number BR10764999. The APC was funded by “Kazakh Research Institute of Livestock and Fodder Production”.

Institutional Review Board Statement

The animal study protocol was approved by the Local Ethics Committee of RSE “Institute of Genetics and Physiology” of the CS MSHE of the Republic of Kazakhstan (No 15-70/10.03.2023).

Informed Consent Statement

Agreements No 53-1/2022 and No 53-2/2022 on project BR10764999 co-financing and cooperation in research activities between peasant farms and “Kazakh Research Institute of Livestock and Fodder Production”.

Data Availability Statement

Not applicable.

Acknowledgments

This article benefited greatly from insightful comments by the Editors and two anonymous reviewers. We are thankful to all of the participants of this project and to J. Robert Macey for help with English editing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lister, A.M.; Kadwell, M.; Kaagan, L.M.; Jordan, W.C.; Richards, M.B.; Stanley, H.F. Ancient and modern DNA in a study of horse domestication. Anc. Biomol. 1998, 2, 267–280. [Google Scholar]
  2. Vilà, C.; Leonard, J.A.; Götherström, A.; Marklund, S.; Sandberg, K.; Liden, K.; Wayne, R.K.; Ellegren, H. Widespread Origins of Domestic Horse Lineages. Science 2001, 291, 474–477. [Google Scholar] [CrossRef] [PubMed]
  3. Achilli, A.; Olivieri, A.; Soares, P.; Lancioni, H.; Hooshiar Kashani, B.; Perego, U.A.; Nergadze, S.G.; Carossa, V.; Santagostino, M.; Capomaccio, S.; et al. Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication. Proc. Natl. Acad. Sci. USA 2012, 109, 2449–2454. [Google Scholar] [CrossRef]
  4. Warmuth, V.; Eriksson, A.; Bower, M.A.; Barker, G.; Barrett, E.; Hanks, B.K.; Li, S.; Lomitashvili, D.; Ochir-Goryaeva, M.; Sizonov, G.V.; et al. Reconstructing the origin and spread of horse domestication in the Eurasian steppe. Proc. Natl. Acad. Sci. USA 2012, 109, 8202–8206. [Google Scholar] [CrossRef]
  5. The Second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture in Brief. Available online: https://www.fao.org (accessed on 23 May 2022).
  6. Horse Milk Has More Health Benefits Than Cow’s Milk—Researchers. Available online: https://guardian.ng/news/horse-milk (accessed on 12 March 2023).
  7. Horse Breeding in Kazakhstan: Development Prospects and Today’s Realities. Available online: https://world-nan.kz/blogs/innovatsionnyyi-granulyator-kormov-predstavili-na-eurotier-2022 (accessed on 10 March 2022). (In Russian).
  8. Dmitriez, N.; Ernst, L. KUSHUM (Kushumskaya), 330–331. In Animal Genetics Resources of the USSR; Dmitriez, N., Ernst, L., Eds.; Food and Agriculture Organization of the United Nations: Roma, Italy, 1989. [Google Scholar]
  9. Kushum Horse. Available online: https://world-animal.com/6333949-kushum-horse (accessed on 15 March 2022).
  10. Rzabayev, S.S. Mugalzhar Horse Breed; LLP Information and Printing Center—Kokzhiyek: Aktobe, Kazakhstan, 2007; p. 154. [Google Scholar]
  11. Satybaldin, A.A. Current state of horse breeding and horse sports in Kazakhstan. In Proceedings of the First International Conference, Kostanay, Kazakhstan, 2002. [Google Scholar]
  12. Moridi, M.; Masoudi, A.A.; Vaez Torshizi, R.; Hill, E.W. Mitochondrial DNA D-loop sequence variation in maternal lineages of Iranian native horses. Anim. Genet. 2013, 44, 209–213. [Google Scholar] [CrossRef] [PubMed]
  13. Bowling, A.T.; Del Valle, A.; Bowling, M. A pedigree-based study of mitochondrial D-loop DNA sequence variation among Arabian horses. Anim. Genet. 2000, 31, 1–7. [Google Scholar] [CrossRef] [PubMed]
  14. Csizmár, N.; Mihók, S.; Jávor, A.; Kusza, S. Genetic analysis of the Hungarian draft horse population using partial mitochondrial DNA D-loop sequencing. PeerJ 2018, 6, e4198. [Google Scholar] [CrossRef]
  15. Effa, K.; Rosenbom, S.; Han, J.; Dessie, T.; Beja-Pereira, A. Genetic Diversities and Historical Dynamics of Native Ethiopian Horse Populations (Equus caballus) Inferred from Mitochondrial DNA Polymorphisms. Genes 2021, 12, 155. [Google Scholar] [CrossRef]
  16. Almarzook, S.; Reissmann, M.; Brockmann, G.A. Diversity of mitochondrial DNA in three Arabian horse strains. J. Appl. Genet. 2017, 58, 273–276. [Google Scholar] [CrossRef]
  17. Hong, J.H.; Oh, C.H.; Kim, S.; Kang, I.U.; Shin, D.H. Genetic analysis of mitochondrial DNA from ancient Equus caballus bones found at archaeological site of Joseon Dynasty period capital area. Anim. Biosci. 2022, 35, 1141–1150. [Google Scholar] [CrossRef]
  18. McGahern, A.M.; Edwards, C.J.; Bower, M.A.; Heffernan, A.; Park, S.D.E.; Brophy, P.O.; Hill, E.W. Mitochondrial DNA sequence diversity in extant Irish horse populations and in ancient horses. Anim. Genet. 2006, 37, 498–502. [Google Scholar] [CrossRef] [PubMed]
  19. Yatkın, S.; Özdil, F.; Özkan Ünal, E.; Genç, S.; Kaplan, S.; Gürcan, E.K.; Arat, S.; Soysal, M.İ. Genetic Characterization of Native Donkey (Equus asinus) Populations of Turkey Using Microsatellite Markers. Animals 2020, 10, 1093. [Google Scholar] [CrossRef] [PubMed]
  20. Kim, S.M.; Yun, S.W.; Cho, G.J. Assessment of genetic diversity using microsatellite markers to compare donkeys (Equus asinus) with horses (Equus caballus). Anim. Biosci. 2021, 34, 1460–1465. [Google Scholar] [CrossRef] [PubMed]
  21. Tozaki, T.; Kakoi, H.; Mashima, S.; Hirota, K.-I.; Hasegawa, T.; Ishida, N.; Miura, N.; Choi-Miura, N.-H.; Tomita, M. Population study and validation of paternity testing for Thoroughbred Horses by 15 microsatellite loci. J. Vet. Med. Sci. 2001, 63, 1191–1197. [Google Scholar] [CrossRef] [PubMed]
  22. Grilz-Seger, G.; Druml, T.; Neuditschko, M.; Dobretsberger, M.; Horna, M.; Brem, G. High-resolution population structure and runs of homozygosity reveal the genetic architecture of complex traits in the Lipizzan horse. BMC Genom. 2019, 20, 174. [Google Scholar] [CrossRef]
  23. Zhang, C.; Ni, P.; Ahmad, H.I.; Gemingguli, M.; Baizilaitibei, A.; Gulibaheti, D.; Fang, Y.; Wang, H.; Asif, A.R.; Xiao, C.; et al. Detecting the Population Structure and Scanning for Signatures of Selection in Horses (Equus caballus) From Whole-Genome Sequencing Data. Evol. Bioinform. 2018, 14, 1176934318775106. [Google Scholar] [CrossRef] [PubMed]
  24. Ma, H.; Wang, S.; Zeng, G.; Guo, J.; Guo, M.; Dong, X.; Hua, G.; Liu, Y.; Wang, M.; Ling, Y.; et al. The Origin of a Coastal Indigenous Horse Breed in China Revealed by Genome-Wide SNP Data. Genes 2019, 10, 241. [Google Scholar] [CrossRef]
  25. Kargayeva, M.T.; Baimukhanov, D.A.; Nurbayev, S.D.; Baimukhanov, A.D.; Alikhanov, O.; Yusupbayev, Z. Identification of kazakh horses by microsatellite DNA using modern analytical methods. Bull. Natl. Acad. Sci. Repub. Kazakhstan 2020, 4, 55–61. [Google Scholar] [CrossRef]
  26. Gemingguli, M.; Iskhan, K.R.; Li, Y.; Qi, W.; Wunirifu, W.; Ding, L.Y.; Wumaierjiang, A. Genetic diversity and population structure of Kazakh horses (Equus caballus) inferred from mtDNA sequences. Genet. Mol. Res. 2016, 15, 15048618. [Google Scholar] [CrossRef]
  27. Seleuova, L.A.; Naimanov, D.K.; Jarowski, Z.; Aubakirov, M.Z.; Mustafin, B.M.; Safronova, O.S.; Baktybayev, G.T.; Turabaev, A.T.; Domatski, V.N. Population genetic characteristic of horses of Mugalzhar breed by STR-markers. Biomed. Res. 2018, 29, 3508–3511. [Google Scholar]
  28. Nguyen, T.B.; Paul, R.C.; Okuda, Y.; Le, T.N.A.; Pham, P.T.K.; Kaissar, K.J.; Kazhmurat, A.; Bibigul, S.; Bakhtin, M.; Kazymbet, P.; et al. Genetic characterization of Kushum horses in Kazakhstan based on haplotypes of mtDNA and Y chromosome, and genes associated with important traits of the horses. J. Equine. Sci. 2020, 31, 35–43. [Google Scholar] [CrossRef] [PubMed]
  29. Kalashnikov, V.; Koveshnikov, V.; Sergienko, S. Prospects for the development of horse breeding in Russia in the new millennium. In Problems and Prospects of Horse Breeding in Russia at the Turn of the Century; Moscow, Russia, 2000; pp. 5–7. [Google Scholar]
  30. Canon, J.; Checa, M.L.; Carleos, C.; Vega-Pla, J.L.; Vallejo, M.; Dunner, S. The genetic structure of Spanish Celtic horse breeds inferred from microsatellite data. Anim. Genet. 2000, 31, 39–48. [Google Scholar] [CrossRef] [PubMed]
  31. Bjornstad, G.; Gunby, E.; Roed, K.H. Genetic structure of Norwegian horse breeds. J. Anim. Breed. Genet. 2000, 117, 307–317. [Google Scholar] [CrossRef]
  32. Aberle, K.S.; Hamann, H.; Drogemuller, C.; Distl, O. Genetic diversity in German drought horse breeds compared with a group of primitive, riding and wild horses by means of microsatellite DNA markers. Anim. Genet. 2004, 35, 270–277. [Google Scholar] [CrossRef] [PubMed]
  33. Zabek, T.; Fornal, A. Evaluation of the 17-plex STR kit for parentage testing of Polish coldblood and Hucul horses. Ann. Anim. Sci. 2009, 9, 363–372. [Google Scholar]
  34. Dimsoski, P. Development of a 17-plex microsatellite polymerase chain reaction kit for genotyping horses. Croat. Med. J. 2003, 44, 332–335. [Google Scholar]
  35. Marshal, T.C.; Slate, J.; Kruuk, L.E.B.; Pemberton, J.M. Statistical confidence for likelihood-based paternity inference in natural populations. Mol. Ecol. 1998, 7, 639–655. [Google Scholar] [CrossRef]
  36. Kalinowski, S.T.; Taper, M.L.; Marshall, T.C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol. Ecol. 2003, 16, 1099–1106. [Google Scholar] [CrossRef]
  37. Peakall, R.; Smouse, P.E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 2006, 6, 288–295. [Google Scholar] [CrossRef]
  38. Dieringer, D.; Schlötterer, C. Microsatellite analyser (MSA): A platform independent analysis tool for large microsatellite data sets. Mol. Ecol. Notes 2003, 3, 167–169. [Google Scholar] [CrossRef]
  39. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [PubMed]
  40. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  41. Belkhir, K. GENETIX 4.05, Logiciel Sous Windows TM Pour la Génétique Des Populations. 2004. Available online: http://www.genetix.univ-montp2.fr/genetix/genetix.htm. (accessed on 17 December 2022).
  42. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef] [PubMed]
  43. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the numbers of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef]
  44. Earl, D.A.; von Holdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  45. Jakobsson, M.; Rosenberg Noah, A. CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 2007, 23, 1801–1806. [Google Scholar] [CrossRef]
  46. Rosenberg, N.A. Distruct: A program for the graphical display of population structure. Mol. Ecol. Notes 2004, 4, 137–138. [Google Scholar] [CrossRef]
  47. Yordanov, G.; Mehandjyiski, I.; Palova, N.; Atsenova, N.; Boyko, N.; Radoslavov, G.; Hristov, P. Genetic diversity and structure of the main Danubian horse paternal genealogical lineages based on microsatellite genotyping. Vet. Sci. 2022, 9, 333. [Google Scholar] [CrossRef]
  48. Jung, J.S.; Seong, J.; Lee, G.H.; Kim, Y.; An, J.H.; Yun, J.H.; Kong, H.S. Genetic diversity and relationship of Halla horse based on polymorphism in microsatellites. J. Anim. Reprod. Biotechnol. 2021, 36, 76–81. [Google Scholar] [CrossRef]
  49. Wright, S. The genetical structure of populations. Ann. Eugen. 1951, 15, 323–354. [Google Scholar] [CrossRef]
  50. Benhamadi, M.E.A.; Berber, N.; Benyarou, M.; Ameur, A.A.; Haddam, H.Y.; Piro, M.; Gaouar, S.B.S. Molecular Characterization of Eight Horse Breeds in Algeria Using Microsatellite Markers. Biodiversitas 2020, 21, d210923. [Google Scholar] [CrossRef]
  51. Dorji, J.; Tamang, S.; Tshewang, T.; Dorji, T.; Dorji, T.Y. Genetic diversity and population structure of three traditional horse breeds of Bhutan based on 29 DNA microsatellite markers. PLoS ONE 2018, 13, e199376. [Google Scholar] [CrossRef] [PubMed]
  52. Yüceer, B.; Erdoğan, M.; Yaralı, C.; Özarslan, B.; Özbeyaz, C. Genetic diversity between Rahvan (pacing) horses in Turkey. Ankara Üniv. Vet. Fak. Derg. 2016, 63, 201–210, (In Turkısh). [Google Scholar]
  53. Ling, Y.H.; Ma, Y.H.; Guan, W.J.; Cheng, Y.J.; Wang, Y.P.; Han, J.L.; Mang, L.; Zhao, Q.J.; He, X.H.; Pu, Y.B.; et al. Evaluation of the genetic diversity and population structure of Chinese indigenous horse breeds using 27 microsatellite markers. Anim. Genet. 2011, 42, 56–65. [Google Scholar] [CrossRef]
  54. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [CrossRef]
  55. Behl, R.; Behl, J.; Gupta, N.; Gupta, S.C. Genetic relationships of five Indian horse breeds using a microsatellite markers. Animal 2007, 1, 483–488. [Google Scholar] [CrossRef]
  56. Solis, A.; Jugo, B.M.; Meriaux, J.C.; Iriondo, M.; Mazon, L.I.; Aguirre, A.I.; Vicario, A.; Estomba, A. Genetic diversity within and among four south European Native Horse breeds based on microsatellite DNA analysis: Implications for conservation. J. Hered. 2005, 96, 670–678. [Google Scholar] [CrossRef]
  57. Voronkova, V.N.; Nikolaeva, E.A.; Piskunov, A.K.; Babayan, O.V.; Takasu, M.; Tozaki, T.; Svishcheva, G.R.; Stolpovsky, Y.A. Assessment of Genetic Diversity and Structure of Russian and Mongolian Autochthonous Horse Breeds Using Nuclear and Mitochondrial DNA Markers. Russ. J. Genet. 2022, 58, 927–943. [Google Scholar] [CrossRef]
  58. Blokhina, N.V.; Khrabrova, L.A.; Gavrilicheva, I.S. Application of modern technologies in identifying distinctive features in the subpopulation of Novoaltaiskaya hoses. IOP Conf. Ser. Earth Environ. Sci. 2021, 624, 012019. [Google Scholar] [CrossRef]
  59. Zaitcev, A.M.; Gavrilicheva, I.S.; Blohina, N.V.; Khrabrova, L.A.; Kokorina, N.B. Assessment of the population structure of horses of the Priobskaya breed based on modern technologies. International Conference on World Technological trends in Agribusiness. IOP Conf. Ser. Earth Environ. Sci. 2021, 624, 012032. [Google Scholar] [CrossRef]
  60. Iwańczyk, E.; Juras, R.; Cholewiński, G.; Cothran, E.G. Genetic structure and phylogenetic relationships of the Polish Heavy horse. J. Appl Genet. 2006, 47, 353–359. [Google Scholar] [CrossRef] [PubMed]
  61. Juras, R.; Cothran, E.G.; Klimas, R. Genetic Analysis of Three Lithuanian Native Horse Breeds. Acta. Agric. Scand. 2003, 53–54, 180–185. [Google Scholar] [CrossRef]
  62. Funk, S.M.; Guedaoura, S.; Juras, R.; Raziq, A.; Landolsi, F.; Luís, C.; Martínez, A.M.; Mayaki, A.M.; Mujica, F.; Oom, M.D.M.; et al. Major inconsistencies of inferred population genetic structure estimated in a large set of domestic horse breeds using microsatellites. Ecol. Evol. 2020, 10, 4261–4279. [Google Scholar] [CrossRef] [PubMed]
  63. Wright, S. Evolution and the genetics of populations, vol. 4. In Variability Within and among Natural Populations; University of Chicago Press: Chicago, IL, USA, 1978; p. 580. [Google Scholar]
  64. Petersen, J.L.; Mickelson, J.R.; Cothran, E.G.; Andersson, L.S.; Axelsson, J.; Bailey, E.; Bannasch, D.; Binns, M.M.; Borges, A.S.; Brama, P.; et al. Genetic diversity in the modern horse illustrated from genome-wide SNP data. PLoS ONE 2013, 8, e54997. [Google Scholar] [CrossRef] [PubMed]
  65. Marchiori, C.M.; Pereira, G.L.; Maiorano, A.M.; Rogatto, G.M.; Assoni, A.D.; Silva, J.A.V., II; Chardulo, L.A.L.; Curi, R.A. Linkage disequilibrium and population structure characterization in the cutting and racing lines of Quarter Horses bred in Brazil. Livest. Sci. 2019, 219, 45–51. [Google Scholar] [CrossRef]
Figure 1. Distribution map of two main Kazakh horse breeds, Kushum and Mugalzhar.
Figure 1. Distribution map of two main Kazakh horse breeds, Kushum and Mugalzhar.
Vetsci 10 00598 g001
Figure 2. Factorial correspondence analysis of 7 horse populations studied on the basis of 11 STR loci. Dashed lines representing the three clusters with the following colors: yellow—Uralsk, Kulandy—black, and waterloo color—Aktobe, Emba, Maupas, Meiman, Mesker.
Figure 2. Factorial correspondence analysis of 7 horse populations studied on the basis of 11 STR loci. Dashed lines representing the three clusters with the following colors: yellow—Uralsk, Kulandy—black, and waterloo color—Aktobe, Emba, Maupas, Meiman, Mesker.
Vetsci 10 00598 g002
Figure 3. Neighbor-joining dendrogram showing the relationships of seven Kazakh horse populations with Nei’s genetic distances plotted. Outline colors of the tree branches representing the breeds according to Figure 1; red = Kushum breed, green = Mugalzhar breed.
Figure 3. Neighbor-joining dendrogram showing the relationships of seven Kazakh horse populations with Nei’s genetic distances plotted. Outline colors of the tree branches representing the breeds according to Figure 1; red = Kushum breed, green = Mugalzhar breed.
Vetsci 10 00598 g003
Figure 4. The best K results of the hierarchical STRUCTURE clustering. (a) Method of [33]; (b) method of [32]; (c) variance of K = 3; (d) K = 4; (e) K = 7. The green color representing the Uralsk population, blue color − Aktobe and Emba populations, orange color − Kozhamberdy type (Maupas, Meiman, Mesker) populations, blue color − Kulandy population.
Figure 4. The best K results of the hierarchical STRUCTURE clustering. (a) Method of [33]; (b) method of [32]; (c) variance of K = 3; (d) K = 4; (e) K = 7. The green color representing the Uralsk population, blue color − Aktobe and Emba populations, orange color − Kozhamberdy type (Maupas, Meiman, Mesker) populations, blue color − Kulandy population.
Vetsci 10 00598 g004
Table 1. Sampling information of studied populations of E. caballus.
Table 1. Sampling information of studied populations of E. caballus.
NPopulation NsLocation
Kushum breed
1Uralsk type76Zhangali district, W KZ
2Aktobe type74Aktobe region, W KZ
Mugalzhar breed
6Emba type93Aktobe region, W KZ
7Kulandy type101Kyzylorda region, SW KZ
Kozhamberdy type
3Maupas interbreed line29Karaganda region, C KZ
4Meiman interbreed line28Karaganda region, C KZ
5Mesker interbreed line34Karaganda region, C KZ
N—consecutive number of populations; Ns—number of samples; W—western; SW—southwestern; C—central; KZ—Kazakhstan.
Table 2. Summary statistics of mean genetic diversity at 11 microsatellite loci in 435 individuals of E. caballus from Kazakhstan.
Table 2. Summary statistics of mean genetic diversity at 11 microsatellite loci in 435 individuals of E. caballus from Kazakhstan.
LocusNANEIHOHEUHEFISFITFSTPIC
VHL2010.576.0261.9910.8550.8300.839−0.0310.0040.0340.813
HTG46.4232.9001.3380.6340.6420.6500.0130.0440.0320.629
AHT48.4295.7011.8670.7460.8230.8320.0930.1250.0350.822
HMS76.5713.9111.4960.6670.7280.7360.0830.1120.0310.734
AHT57.0004.9081.6950.7790.7960.8050.0210.0510.0300.784
HMS67.5714.6601.7170.6520.7830.7930.1680.2060.0460.788
ASB239.8575.6781.9180.5020.8220.8370.3890.4150.0420.713
ASB210.006.3341.9270.8420.8060.896−0.0450.0330.0740.627
HTG77.0003.5681.4530.7510.7120.720−0.054−0.0100.0410.686
HMS38.1434.6761.7110.5660.7840.7930.2780.3060.0400.786
LEX39.5715.7281.9120.4790.8140.8250.4220.4390.0470.768
Mean 8.2864.9171.7300.6790.7760.7930.1210.1570.0410.741
Numbers of observed alleles (NA), number of effective alleles (NE), Shannon’s information index (I), observed heterozygosity (HO), expected heterozygosity (HE), unbiased expected heterozygosity (UHE), inbreeding coefficient (FIS), fixation index (FIT), population differentiation statistic (FST), p ≤ 0.001.
Table 3. Analysis of molecular variance.
Table 3. Analysis of molecular variance.
Source of Variationd.f.SSMSEst. Var.% of Variation
Among populations6240.41140.0690.2906%
Among individuals4282270.1835.3041.07424%
Within individuals4351372.5003.1553.15570%
Total8693883.094 4.520100%
FST = 0.064 (p ≥ 0.001); FIS = 0.254 (p ≥ 0.001); FIT = 0.302 (p ≥ 0.001)
d.f.—degrees of freedom; SS—sum of squares; MS—mean of squares; Est. Var.—estimated variance.
Table 4. Genetic distances of studied horse populations.
Table 4. Genetic distances of studied horse populations.
UralskAktobeMaupasMeimanMeskerEmbaKulandy
Uralsk
Aktobe0.045
Maupas0.0370.021
Meiman0.0350.0200.008
Mesker0.0410.0210.0100.014
Emba0.0330.0140.0140.0140.021
Kulandy0.0470.0250.0220.0210.0280.021
The Wright’s FST fixation indices given below diagonal.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Orazymbetova, Z.; Ualiyeva, D.; Dossybayev, K.; Torekhanov, A.; Sydykov, D.; Mussayeva, A.; Baktybayev, G. Genetic Diversity of Kazakhstani Equus caballus (Linnaeus, 1758) Horse Breeds Inferred from Microsatellite Markers. Vet. Sci. 2023, 10, 598. https://doi.org/10.3390/vetsci10100598

AMA Style

Orazymbetova Z, Ualiyeva D, Dossybayev K, Torekhanov A, Sydykov D, Mussayeva A, Baktybayev G. Genetic Diversity of Kazakhstani Equus caballus (Linnaeus, 1758) Horse Breeds Inferred from Microsatellite Markers. Veterinary Sciences. 2023; 10(10):598. https://doi.org/10.3390/vetsci10100598

Chicago/Turabian Style

Orazymbetova, Zarina, Daniya Ualiyeva, Kairat Dossybayev, Aibyn Torekhanov, Dauren Sydykov, Aizhan Mussayeva, and Gabiden Baktybayev. 2023. "Genetic Diversity of Kazakhstani Equus caballus (Linnaeus, 1758) Horse Breeds Inferred from Microsatellite Markers" Veterinary Sciences 10, no. 10: 598. https://doi.org/10.3390/vetsci10100598

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop