Next Article in Journal
Prioritizing Plants around the Cross-Border Area of Greece and the Republic of North Macedonia: Integrated Conservation Actions and Sustainable Exploitation Potential
Next Article in Special Issue
The Spread of the Invasive Locust Digitate Leafminer Parectopa robiniella Clemens, 1863 (Lepidoptera: Gracillariidae) in Europe, with Special Reference to Ukraine
Previous Article in Journal
Prey Selectivity in Juvenile Red King Crabs from the Coastal Barents Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Genetic Diversity of Dagestan Mountain Cattle Using STR-Markers

by
Valeria V. Volkova
1,*,
Alexandra S. Abdelmanova
1,
Tatiana E. Deniskova
1,
Olga S. Romanenkova
1,
Abdusalam A. Khozhokov
2,
Alimsoltan A. Ozdemirov
2,
Alexander A. Sermyagin
1 and
Natalia A. Zinovieva
1
1
L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy 60, Podolsk Municipal District, 142132 Moscow, Moscow Region, Russia
2
Federal State Budgetary Scientific Institution “Federal Agrarian Scientific Center of the Republic of Dagestan”, Akushinsky Ave., Scientific Town, 367014 Makhachkala, Republic of Dagestan, Russia
*
Author to whom correspondence should be addressed.
Submission received: 17 June 2022 / Revised: 14 July 2022 / Accepted: 14 July 2022 / Published: 16 July 2022

Abstract

:
The Dagestan Mountain (DM) cattle breed was established to provide cheese and milk products to mountain dwellers in the specific conditions of the Republic of Dagestan in Southern Russia. Only 650 head of DM cattle were registered in 2020, and their “breed purity” is questionable. We aimed to assess the genetic diversity and population structure of modern DM cattle using short tandem repeat (STR) markers. The sample included 150 animals collected from private owners in Dagestan during a scientific expedition (n = 32) and provided by a gene pool farm (n = 118). An additional 166 samples from other cattle breeds distributed in the breeding zone of DM cattle were used as the comparison groups. The genotypes for the 11 STR loci recommended by ISAG were obtained using a genetic analyzer. We identified a high level of genetic diversity as revealed by allelic richness calculations (Ar = 6.82 vs. 4.38–5.82 in other cattle breeds) and observed heterozygosity indices (uHe = 0.76 vs. 0.65–0.72). Based on the STRUCTURE results, animals with low levels of admixture with other breeds were found within the DM cattle, which can be considered as candidates for use in germplasm conservation programs.

1. Introduction

Dagestan is located in the southernmost part of the Russian Federation on the north-eastern slopes of the Greater Caucasus and Caspian lowlands. The territory of Dagestan comprises several landscape zones, including plains, foothills, mountains, mountain valleys, and highlands. The mountainous zone occupies almost half of the territory (48%) [1]. A study of the labor activity of the ethnic population of Dagestan showed that the source of the highlanders’ livelihood was manual labor, which was complicated by unfavorable natural conditions, such as low atmospheric pressure, rarefied air, lack of oxygen, sudden temperature changes (day–night), and difficult terrain [2]. Simultaneously, large areas of mountain pastures, cheap pasture fodder, and a long pasture season favor cattle breeding. The inhabitants of the plain and foothill zones prefer cattle breeding, whereas mountain dwellers raise both sheep and cattle. Livestock management systems differ depending on the landscape zone. Combined pasture and stall-cattle keeping are traditional management systems in the Dagestan flatlands. From spring to late autumn, the cattle are kept on rural pastures with no feed supplements. Residents of some communities move cattle from rural to mountain pastures from spring to autumn, when raw milk is processed into butter, cheese, and cottage cheese. A transhumance system is used in the mountains [2].
The Republic of Dagestan has a long history of livestock husbandry. Great Caucasian and Lesser Caucasian cattle were raised in the highlands and flatlands until the 1930s; currently, these cattle are reared in several regions. The Dagestan Mountain (DM) breed was obtained by crossings of local and brown Swiss cattle with further improvement by the Kostroma and Lebedinsky breeds [3] (Figure 1).
DM cattle represent a valuable gene pool, and are well adapted to breeding in extremal highlands in the climate of southern Russia. Owing to their comparably low body weight, DM cattle effortlessly move around the highlands and consume vegetation in places that are inaccessible to other cattle breeds. In addition to adapting to environmental conditions, the most noticeable features of mountain cattle are their obedience and their ability to herd without human supervision. They are released from the byre in the morning to find their way through the outskirts of villages and mountain slopes as they forage for food. In the evening, they return to the byres [4]. These features facilitate the obtention of milk and meat from these animals without spending any resources on care and feeding. These valuable qualities allowed highlanders to survive harsh conditions for many centuries.
The human population has been growing steadily since the 20th century [5]. Extensive livestock farming can no longer meet the demand for animal proteins. Highland cattle cannot compete with highly productive commercial breeds owing to their low milk yield and fattening ability. Therefore, the breeding of these animals has not received proper attention in the Republic of Dagestan. Poor fodder resources and insufficient selection are the main reasons for the low productivity of DM cattle [6]. The intensification of cattle breeding is unsuccessful because of the specific natural geographic conditions in highland areas. The breeding of specialized commercial cattle in mountainous terrain leads to a considerable decrease in milk yield and early culling for various reasons [7]. In this regard, conservation of the DM cattle gene pool is an important priority for livestock husbandry in mountainous terrains.
Only 650 head of DM cattle were registered in 2020 and kept at a single gene pool farm in Dagestan. Most cattle lack complete pedigree records; therefore, their “breed purity” is questionable. We aimed to characterize the current gene pool, evaluate genetic diversity, and analyze the population structure of DM cattle in comparison with other cattle breeds reared in Russia based on microsatellite analysis.

2. Materials and Methods

2.1. Sample Collection

Tissue samples of DM cattle (n = 32) were collected from private owners in the mountain villages of Dagestan during a scientific expedition in 2019. In 2021, the gene pool farm involved in the conservation of DM cattle provided 118 additional samples (Supplementary Table S1 and Supplementary Figure S1). The final sample consisted of 150 specimens. The short tandem repeat (STR) genotypes of five other cattle breeds (n = 166) bred in this breeding zone, which could have contributed to the gene pool of DM cattle, were obtained from the Bioresource Collection of Farm Animals of the L.K. Ernst Federal Science Center for Animal Husbandry, supported by the Russian Ministry of Science and Higher Education [8] (Table 1).

2.2. DNA Extraction and STR-Genotyping

The genomic DNA was extracted using the DNA-Extran-2 and S-Sorb Kits (SyntolTM, Moscow, Russia) according to the manufacturer’s instructions.
In this study, we used eleven microsatellite loci, which were recommended by the International Society for Animal Genetics (ISAG) [9] and included BM1818, BM2113, BM1824, [10], ETH10 [11], ETH225 [12], INRA023 [13], SPS115 [14] and TGLA53, TGLA122, TGLA126, TGLA227 [15]. All microsatellite loci were autosomal. The 5′ fluorescently labeled primer sequences and chromosomal localization of each microsatellite locus are provided in Supplementary Table S2.
PCR amplification was performed according to standard protocol [16]. The total reaction volume of 9 µL comprised 1 µL 10× PCR-buffer, 1 µL dNTP (1 mmol/L), 0.5 µL MgCl2, 0.34 µL of primer master mix, 0.1 µL of Smart Taq DNA polymerase (Dialat ltd, Moscow, Russia), and 1 µL of DNA template. The PCR amplification included initial denaturation (95 °C for 10 min), 40 cycles (95 °C for 30 s, 60 °C for 40 s, and 72 °C for 1 min), and a final extension (72 °C for 10 min).
Microsatellite analysis was performed using ABI 3130xl Genetic Analyzer (Applied Biosystems, Waltham, MA, USA). Raw allele sizes were determined using the GeneMapper v.4.0 software (Applied Biosystems, Waltham, MA, USA).

2.3. Data Analysis

Allelic and genetic diversity parameters were calculated using GenAIEx 6.503 software [17] and the R package diversity [18].
To study the genetic differentiation between DM and other cattle populations, we calculated pairwise FST values [19] and Jost’s index of population differentiation (Jost’s D) [20] using the R Project for Statistical Computing software environment [21].
The Neighbor-net graph, which is based on Jost’s index of population differentiation, was constructed using SplitsTree 4.14.5 software [22].
Principal component analysis (PCA) was performed using the R package adegenet [23] and visualized using the R package ggplot2 [24]. The map with sampling sites was created using R package maps [25]. Data files were prepared using R 3.5.0 [21].
The genetic structure of the populations was estimated in Structure 2.3.4 software [26] using the following parameters: burn-in period—10,000; number of Markov chain Monte Carlo simulations (MCMC)—100,000 for each run. Ten iterations were performed for each K. CLUMPAK software [27,28] was used to visualize and to determine the most probable number of clusters in the studied sample based on the ΔK values according to the method proposed by Evanno et al. [29]. To analyze the population structure, we determined the average similarity scores produced in CLUMPAK for several independent runs with the same K-value.

3. Results

3.1. Genetic Variation among and within Breeds

The statistics of the analyzed loci are summarized in Table 2. The total number of alleles per locus ranged from 26 (BM1824) to 61 (TGLA53) alleles.
In total, 442 alleles across 11 loci and 34 private alleles across 10 loci were identified in the studied populations. No private alleles were found at the BM1818 locus. The DM population was characterized by the highest number of alleles (118 alleles, including 25 private alleles).
Table 3 presents the main statistical indicators used to estimate the current allele pool and the level of genetic diversity of the studied cattle populations.
We observed the highest level of genetic diversity in DM cattle compared to other breeds, as revealed by allelic richness (Ar = 6.827 ± 0.654 vs. 4.376–5.818) and unbiased expected heterozygosity (uHe = 0.764 ± 0.031 vs. 0.645–0.718). The breeds raised in the Caucasus showed a higher level of genetic diversity (Ar values varied from 5.278 ± 0.45 to 6.827 ± 0.654) compared to transboundary breeds (Ar = 4.376 ± 0.396 to 5.307 ± 0.373), which can be the consequence of lower selection pressure. We detected a larger heterozygote deficiency in DM cattle (uFis = 0.052), which can be caused by the low population size. The SIM and BS_G groups revealed an excess of heterozygotes (uFis = −0.073 and uFis = −0.058, respectively), while the heterozygote number of other cattle breeds did not show a significant deviation from the Hardy–Weinberg equilibrium.

3.2. Genetic Differentiation

Table 4 shows Jost’s D genetic distances and FST values, which characterize the degree of differentiation between the studied cattle populations, calculated for each pair of groups.
According to the classification [30], most pairwise FST values between the studied cattle populations corresponded to moderate genetic differentiation (0.053–0.138). We established that the BS_G, HOL, and SIM groups had the greatest genetic distance from the other populations.
The DM cattle demonstrated moderate differentiation in all groups, with the least genetic distance from the BS_D group (FST = 0.023), indicating the participation of DM cattle in the development of the Caucasian type of Brown Swiss cattle.
PCA (Figure 2a) showed that the first principal component (PC1), responsible for 5.62% of the genetic variability, clearly separated the Holstein breed from the other breeds. The second principal component (PC2) accounted for 3.83% of the genetic variability and separated the Simmental breed from the group of breeds, which were reared mainly in the Caucasus (DM, BS, and RS). The HOL and BS_G breeds formed more compact clusters, which may indicate higher selection pressure in these groups. Brown and red cattle populations, which are bred in the Caucasus, are characterized by greater genetic diversity than the transboundary breeds.
Analysis of the Neighbor-net graph (Figure 2b) showed the presence of genetic links between the DM, BS_D, and RS groups and provided evidence of the common origin of the BS_D and BS_G groups. The localization of the DM cattle at the edge of the Neighbor-net graph suggests its mixed origin. The HOL and SIM breeds form separate branches.
An analysis of the population structure (Figure 3) was performed using Structure for the number of clusters K from two to six. According to the algorithm proposed by Evanno et al. [29], the most probable number of clusters was three. At K = 2 and K = 3 (Supplementary Figure S2), we identified clusters corresponding to HOL and SIM breeds. With a further increase in K, at K = 6, BS, RS, HOL, and SIM groups formed clusters, whereas DM cattle revealed a mixed origin from several ancestral populations. We identified DM antagonists in the Dagestan population of Brown Swiss cattle, which reflects its contribution to the development of the allele pool BS_D cattle. Most of the DM animals did not show common ancestors with other studied breeds, while in a small number of animals, the admixture with BS and RS ancestors was visible. As a result of Structure analysis, DM samples were divided into two clusters. One was linked to the samples collected in the remote mountainous area, and the other was associated with the gene pool population of DM cattle. Some animals showed mixed origins.

4. Discussion

Economic significance, the high demand for dairy products and beef, and diverse environmental conditions have promoted the development of various local cattle breeds in Russia [31]. The breeds that display high productivity and satisfactory climate resilience settled widely outside the regions of their origin, while the other breeds are reared in their historical locations and are in a critical condition as a result of their small population sizes [3]. In addition, endemic breeds inhabit regions with extremely harsh environments and cannot be substituted by other local or commercial breeds without a significant increase in feeding and maintenance costs. DM cattle belong to the last group and dwell predominantly at high altitudes.
The Republic of Dagestan is a mountainous area in the Caucasus, with an abundance of natural highland and flatland pastures. Human tribes have inhabited the Caucasus since the Early Upper Paleolithic Age and most likely imported their cattle [32,33]. However, scientific data on the origin and history of DM cattle are limited and ambiguous. DM cattle belong to a branch of the Greater Caucasian group, which originated from the old dwarf cattle of Egypt [34]. The Greater Caucasian group (or Great Caucasus) included several local types, which were raised in different climatic zones within the mountainous area in the former USSR (Russia, Georgia, and Armenia) [35]. However, another study showed that local cattle, which were improved by the Brown Swiss, Kostroma, and Lebedinsky breeds, were transformed into the Lesser Caucasus cattle type [35], which corresponds to the official breed description [3].
Most local Russian cattle breeds have been investigated using STR [36,37] and single nucleotide polymorphism (SNP) markers [38,39,40,41]. However, little is known about the genetic diversity and population structure of DM cattle inhabiting the high-altitude warmer area of the Republic of Dagestan in southern Russia. Here, we present the results of the first genetic assessment of DM in cattle based on 11 microsatellite markers.
The level of genetic diversity of the DM cattle was the highest among the studied cattle populations and was compatible with the estimates obtained for other aboriginal cattle, including Indian cattle (Ho = 0.681–0.721, He = 0.702–0.751) [42], Creole (Ho = 0.719, He = 0.739) [43], and Taiwanese local cattle breeds (Ho = 0.530–0.607, He = 0.697–0.719) [44]. The DM group showed a moderate level of heterozygote deficiency. Such a pattern is frequent in native breeds [45,46,47,48], which may correspond to a lack of sires [42].
According to the classification of Hartl [30], most of the FST values among the studied cattle populations exhibit a moderate level of genetic differentiation, which agrees with the patterns of genetic distances between local breeds in other countries [44,48]. Pairwise FST values between Indian local cattle populations were lower (FST= 0.008–0.044) [42] than those obtained in our study between DM and other Russian cattle populations, which probably reflects higher selection pressure.
In our study, we found that the DM cattle group had the lowest differentiation evaluated by FST and Jost’s D-values with the Brown Swiss population raised in the same breeding zone. The definite genetic closeness between these groups was demonstrated using the Neighbor-net graph (Figure 2b). In addition, shared genetic components were identified in the DM and Brown Swiss groups using structural analysis (Table 3). Several animals within the DM population («endemic groups») demonstrated the presence of a genetic background, which was absent or found in barely noticeable traces in other studied cattle groups (Figure 3). These individuals can be considered valuable national genetic resources that retain most of their ancestral genetic components. Further studies using more powerful molecular genetic instruments at the genome-wide level are necessary to determine the source of this ancestry more precisely.

5. Conclusions

Our study provides data on the genetic diversity, population structure, and genetic differentiation of modern DM cattle for the first time. Conservation of valuable gene pools of DM cattle is necessary for the further development of animal husbandry in the Dagestan Mountains, as well as for the breeding and selection of animals with high productivity and adaptation to harsh conditions in the mountain areas. Our studies will be continued with a larger sample size using SNP genotyping arrays and whole-genome sequencing.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d14070569/s1, Figure S1: The map with sampling sites of Dagestan Mountain cattle; Figure S2: Population structure of the studied cattle breeds for K from 2 to 6; Table S1: STR-genotypes for the animals of Dagestan Mountain cattle; Table S2: The 5′ fluorescently labeled primer sequences and chromosomal localization of each microsatellite locus used in this study.

Author Contributions

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

Funding

This research was funded by the Ministry of Science and Higher Education of the Russian Federation (theme No. FGGN-2022-0002).

Institutional Review Board Statement

The animal study protocol was approved by the Ethics Committee of the L.K. Ernst Federal Research Center for Animal Husbandry (protocol № 1B of 21 January 22).

Data Availability Statement

The DNA genotypes of Dagestan Mountain cattle used in this study are presented in Supplementary Table S1. The DNA genotypes for other cattle breeds are available on request to the corresponding author.

Acknowledgments

To conduct this study, we used the equipment of the Center for Collective Use “Bioresources and Bioengineering of Farm Animals” of L.K. Ernst Federal Research Center for Animal Husbandry.

Conflicts of Interest

The authors declare no conflict of interest. The Ministry of Science and Higher Education of the Russian Federation had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ibragimov, R.E.; Dzhalalov, A.P.; Alilov, M.M. Outlook for the development to mountain beef cattle breeding in Dagestan. Collect. Sci. Pap. North Cauc. Res. Inst. Anim. Husb. 2013, 2, 68–74. [Google Scholar]
  2. Nuratinov, D.A. Traditional Forms of Management in Dagestan. Issues Struct. Econ. 2008, 4, 22–29. [Google Scholar]
  3. Dunin, I.M.; Dankvert, A.G. Directory of Breeds and Types of Farm Animals Bred in the Russian Federation; VNIIPLEM: Moscow, Russia, 2013; 560p. [Google Scholar]
  4. Petherbridge, G.; Ismailov, M.M.; Ismailov, S.M.; Rabadanov, M.K.; Gadzhiev, A.A.; Teymurov, A.A.; Rabadanov, M.R.; Daudova, M.G.; Abdulaev, A.; Verkhnee Gakvari, G.M. The contribution of adat, religious beliefs and public education to collective environmental management in an agro-pastoral community in the Dagestan high Caucasus. South Russ. Ecol. Dev. 2021, 16, 142–179. [Google Scholar] [CrossRef]
  5. Population United Nations. Available online: https://www.un.org/en/global-issues/population (accessed on 6 June 2022).
  6. Aligazieva, P.A.; Aligaziev, A.M.; Alidibirov, A.T. Influence of Factors on Milk Productivity of First-Heifer Heifers of Mountain Cattle. In The Current State and Main Directions of Development of Seed Production in the Republic of Dagestan: Materials of the All-Russian Scientific and Practical Conference; Publisher Federal State Budgetary Educational Institution of Higher Education “Dagestan State Agrarian University named after M.M. Dzhambulatov”: Makhachkala, Russia, 2019; pp. 119–125. [Google Scholar]
  7. Sadykov, M.M. Ways to increase the meat productivity of mountain cattle. Mt. Agric. 2016, 3, 167–170. [Google Scholar]
  8. UNU. Available online: https://www.vij.ru/infrastruktura/unu (accessed on 20 February 2020).
  9. FAO. Molecular Genetic Characterization of Animal Genetic Resources; FAO: Rome, Italy, 2011; pp. 68–69. Available online: http://www.fao.org/3/i2413e/i2413e00.pdf (accessed on 6 June 2022).
  10. Bishop, M.D.; Kappes, S.M.; Keele, J.W.; Stone, R.T.; Sunden, S.L.; Hawkins, G.A.; Toldo, S.S.; Fries, R.; Grosz, M.D.; Yoo, J.; et al. A genetic linkage map for cattle. Genetics 1994, 136, 619–639. [Google Scholar] [CrossRef]
  11. Toldo, S.S.; Fries, R.; Steffen, P.; Neiberg, H.L.; Barendse, W.; Womack, J.E.; Hetzel, D.J.S.; Stranzinger, G. Physically mapped, cosmid-derived microsatellite markers as anchor loci on bovine chromosomes. Mamm. Genome 1993, 4, 720–727. [Google Scholar] [CrossRef]
  12. Steffen, P.; Eggen, A.; Stranzinger, G.; Fries, R.; Dietz, A.B.; Womack, J.E. Isolation and mapping of polymorphic microsatellites in cattle. Anim. Genet. 2009, 24, 121–124. [Google Scholar] [CrossRef]
  13. Vaiman, D.; Mercier, D.; Moazami-Goudarzi, K.; Eggen, A.; Ciampolini, R.; Lepingle, A.; Velmala, R.; Kaukinen, J.; Varvio, S.L.; Martin, P.; et al. A set of 99 cattle microsatellites: Characterization, synteny mapping, and polymorphism. Mamm. Genome 1994, 5, 288–297. [Google Scholar] [CrossRef]
  14. Moore, S.S.; Byrne, K.; Berger, K.T.; Barendse, W.; McCarthy, F.; Womack, J.E.; Hetzel, D.J.S. Characterization of 65 bovine microsatellites. Mamm. Genome 1994, 5, 84–90. [Google Scholar] [CrossRef]
  15. Barendse, W.; Armitage, S.M.; Kossarek, L.M.; Shalom, A.; Kirkpatrick, B.W.; Ryan, A.M.; Clayton, D.; Li, L.; Neibergs, H.L.; Zhang, N.; et al. A genetic linkage map of the bovine genome. Nat. Genet. 1994, 6, 227–235. [Google Scholar] [CrossRef]
  16. Zinovieva, N.A.; Popov, A.N.; Ernst, L.K. Guidelines for the Use of the Polymerase Chain Reaction Method in Animal Husbandry; VIJ: Dubrovitsy, Russia, 1998; 47p. [Google Scholar]
  17. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Keenan, K.; McGinnity, P.; Cross, T.F.; Crozier, W.W.; Prodöhl, P.A. diveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol. Evol. 2013, 4, 782–788. [Google Scholar] [CrossRef] [Green Version]
  19. Weir, B.S.; Cockerham, C.C. Estimating F-Statistics for the Analysis of Population Structure. Evolution 1984, 38, 1358–1370. [Google Scholar] [CrossRef] [PubMed]
  20. Jost, L. GST and its relatives do not measure differentiation. Mol. Ecol. 2008, 17, 4015. [Google Scholar] [CrossRef]
  21. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: https://www.R-project.org/ (accessed on 23 July 2020).
  22. Huson, D.H.; Bryant, D. Application of Phylogenetic Networks in Evolutionary Studies. Mol. Biol. Evol. 2006, 23, 254–267. [Google Scholar] [CrossRef]
  23. Jombart, T. Adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics 2008, 24, 1403–1405. [Google Scholar] [CrossRef] [Green Version]
  24. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis (Use R!); Springer: New York, NY, USA, 2009; ISBN 978-0-387-98140-6. [Google Scholar]
  25. Becker, R.A.; Wilks, A.R.; Brownrigg, R.; Minka, T.P.; Deckmyn, A. Maps: Draw Geographical Maps. R Package Version 3.3.0. 2018. Available online: https://CRAN.R-project.org/package=maps (accessed on 22 December 2020).
  26. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  27. Kopelman, N.M.; Mayzel, J.; Jakobsson, M.; Rosenberg, N.A.; Mayrose, I. CLUMPAK: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 2015, 15, 1179–1191. [Google Scholar] [CrossRef] [Green Version]
  28. CLUMPAK—Clustering Markov Packager across K. Available online: http://clumpak.tau.ac.il (accessed on 23 May 2022).
  29. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  30. Hartl, D.L.; Clark, A.G. Principles of Population Genetics; Sinauer Associates: Sunderland, UK, 1997. [Google Scholar]
  31. Dmitriev, N.G.; Ernst, L.K. Animal Genetic Resources of the USSR; FAO: Rome, Italy, 1989; Available online: http://www.fao.org/3/ah759e/AH759E07.htm (accessed on 25 May 2022).
  32. Wang, C.; Reinhold, S.; Kalmykov, A.; Wissgott, A.; Brandt, G.; Jeong, C.; Cheronet, O.; Ferry, M.; Harney, E.; Keating, D.; et al. Ancient human genome-wide data from a 3000-year interval in the Caucasus corresponds with eco-geographic regions. Nat. Commun. 2019, 10, 590. [Google Scholar] [CrossRef] [Green Version]
  33. Yunusbayev, B.; Metspalu, M.; Järve, M.; Kutuev, I.; Rootsi, S.; Metspalu, E.; Behar, D.M.; Varendi, K.; Sahakyan, H.; Khusainova, R.; et al. The Caucasus as an asymmetric semipermeable barrier to ancient human migrations. Mol. Biol. Evol. 2012, 29, 359–365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Caucasian Cattle. (n.d.) The Great Soviet Encyclopedia. 3rd Edition. (1970–1979). Available online: https://encyclopedia2.thefreedictionary.com/Caucasian+Cattle (accessed on 23 July 2020).
  35. Porter, V.; Alderson, L.; Hall, S.J.; Sponenberg, D.P. Mason’s World Encyclopedia of Livestock Breeds and Breeding; CABI: Wallingford, UK, 2016; Volume 1–2. [Google Scholar]
  36. Abdelmanova, A.S.; Kharzinova, V.R.; Volkova, V.V.; Mishina, A.I.; Dotsev, A.V.; Sermyagin, A.A.; Boronetskaya, O.I.; Petrikeeva, L.V.; Chinarov, R.Y.; Brem, G.; et al. Genetic Diversity of Historical and Modern Populations of Russian Cattle Breeds Revealed by Microsatellite Analysis. Genes 2020, 11, 940. [Google Scholar] [CrossRef] [PubMed]
  37. Abdelmanova, A.S.; Kharzinova, V.R.; Volkova, V.V.; Dotsev, A.V.; Sermyagin, A.A.; Boronetskaya, O.I.; Chinarov, R.Y.; Lutshikhina, E.M.; Sölkner, J.; Brem, G.; et al. Comparative study of the genetic diversity of local steppe cattle breeds from Russia, Kazakhstan and Kyrgyzstan by microsatellite analysis of museum and modern samples. Diversity 2021, 13, 351. [Google Scholar] [CrossRef]
  38. Sermyagin, A.A.; Dotsev, A.V.; Gladyr, E.A.; Traspov, A.A.; Deniskova, T.E.; Kostyunina, O.V.; Reyer, H.; Wimmers, K.; Barbato, M.; Paronyan, I.A.; et al. Whole-genome SNP analysis elucidates the genetic structure of Russian cattle and its relationship with eurasian taurinebreeds. Genet. Sel. Evol. 2018, 50, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Yurchenko, A.; Yudin, N.; Aitnazarov, R.; Plyusnina, A.; Brukhin, V.; Soloshenko, V.; Lhasaranov, B.; Popov, R.; Paronyan, I.A.; Plemyashov, K.V.; et al. Genome-wide genotyping uncovers genetic profiles and history of the Russian cattle breeds. Heredity 2018, 120, 125–137. [Google Scholar] [CrossRef]
  40. Zinovieva, N.A.; Dotsev, A.V.; Sermyagin, A.A.; Deniskova, T.E.; Abdelmanova, A.S.; Kharzinova, V.R.; Sölkner, J.; Reyer, H.; Wimmers, K.; Brem, G. Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis. PLoS ONE 2020, 15, e0242200. [Google Scholar] [CrossRef]
  41. Zinovieva, N.A.; Sheiko, I.P.; Dotsev, A.V.; Sheiko, R.I.; Mikhailova, M.E.; Sermyagin, A.A.; Abdelmanova, A.S.; Kharzinova, V.R.; Reyer, H.; Wimmers, K.; et al. Genome-wide SNP analysis clearly distinguished the Belarusian Red cattle from other European cattle breeds. Anim. Genet. 2021, 52, 720–724. [Google Scholar] [CrossRef]
  42. Sharma, R.; Maitra, A.; Singh, P.K.; Tantia, M.S. Genetic diversity and relationship of cattle populations of East India: Distinguishing lesser known cattle populations and established breeds based on STR markers. SpringerPlus 2013, 2, 359. [Google Scholar] [CrossRef] [Green Version]
  43. Ginja, C.; Gama, L.T.; Cortes, Ó.; Delgado, J.V.; Dunner, S.; García, D.; Landi, V.; Martín-Burriel, I.; Martínez-Martínez, A.; Penedo, M.C.T.; et al. Analysis of conservation priorities of Iberoamerican cattle based on autosomal microsatellite markers. Genet. Sel. Evol. 2013, 45, 35. [Google Scholar] [CrossRef] [Green Version]
  44. Tu, P.A.; Lin, D.Y.; Li, G.F.; Huang, J.C.; Wang, D.C.; Wang, P.H. Characterization of the genetic diversity and population structure for the yellow cattle in Taiwan based on microsatellite markers. Anim. Biotechnol. 2014, 25, 234–249. [Google Scholar] [CrossRef]
  45. Ginja, C.; Gama, L.T.; Penedo, M.C. Analysis of STR markers reveals high genetic structure in Portuguese native cattle. J. Hered. 2010, 101, 201–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Laoun, A.; Harkat, S.; Lafri, M.; Bechir, S.; Gaouar, S.; Belabdi, I.; Ciani, E.; De Groot, M.; Blanquet, V.; Leroy, G.; et al. Inference of Breed Structure in Farm Animals: Empirical Comparison between SNP and Microsatellite Performance. Genes 2020, 11, 57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Martín-Burriel, I.; Rodellar, C.; Cañón, J.; Cortés, O.; Dunner, S.; Landi, V.; Martínez-Martínez, A.; Gama, L.T.; Ginja, C.; Penedo, M.C.; et al. Genetic diversity, structure, and breed relationships in Iberian cattle. J. Anim. Sci. 2011, 89, 893–906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Ndiaye, N.P.; Sow, A.; Dayo, G.K.; Ndiaye, S.; Sawadogo, G.J.; Sembène, M. Genetic diversity and phylogenetic relationships in local cattle breeds of Senegal based on autosomal microsatellite markers. Vet. World 2015, 8, 994–1005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Typical Dagestan Mountain cow (Republic of Dagestan). Photo provided by Alexander A. Sermyagin.
Figure 1. Typical Dagestan Mountain cow (Republic of Dagestan). Photo provided by Alexander A. Sermyagin.
Diversity 14 00569 g001
Figure 2. Genetic relationships between the studied cattle breeds: (a) Principal component analysis, X-axis: first component (PC1). Y-axis: second component (PC2); (b) dendrogram based on pairwise genetic distances (DJost), plotted using the Neighbor-net algorithm.
Figure 2. Genetic relationships between the studied cattle breeds: (a) Principal component analysis, X-axis: first component (PC1). Y-axis: second component (PC2); (b) dendrogram based on pairwise genetic distances (DJost), plotted using the Neighbor-net algorithm.
Diversity 14 00569 g002
Figure 3. Population structure of the studied cattle breeds. DM—Dagestan Mountain cattle; BS_D—sBrown Swiss, Caucasian branch; BS_G—Brown Swiss, Germany; RS—Red Steppe; SIM—Simmental; HOL—Holstein.
Figure 3. Population structure of the studied cattle breeds. DM—Dagestan Mountain cattle; BS_D—sBrown Swiss, Caucasian branch; BS_G—Brown Swiss, Germany; RS—Red Steppe; SIM—Simmental; HOL—Holstein.
Diversity 14 00569 g003
Table 1. Sampling information for studied cattle breeds.
Table 1. Sampling information for studied cattle breeds.
BreedCodenBreeding Region
Dagestan MountainDM150The Republic of Dagestan
Brown Swiss, Caucasian branchBS_D13The Republic of Dagestan
Brown Swiss, GermanyBS_G27Germany
Red SteppeRS26The Republic of Dagestan
SimmentalSIM50Oryol region
Holstein 1HOL50Holstein Association USA, Inc.
1 STR-genotypes of Holstein cattle breed were obtained from the Holstein Association USA, Inc. (Brattleboro, VT, USA).
Table 2. Allelic variability in studied breeds.
Table 2. Allelic variability in studied breeds.
LocusPopulationIn Total 1
DMBS_DBS_GRSHOLSIM
TGLA22712/17697/1849/2
BM211311/46765540/4
TGLA5317/299871161/2
ETH108/15/1456533/2
SPS1159/24564634/2
TGLA12217/668/27/17/1954/10
INRA2311/175848/143/2
TGLA1268/343446/129/4
BM181885634632/0
ETH22511/47/1575641/5
BM18246/14443526/1
Total118/2564/262/267/156/275/2442/34
1 Number of alleles/number of private alleles.
Table 3. Genetic diversity parameters estimated for 11 microsatellite markers.
Table 3. Genetic diversity parameters estimated for 11 microsatellite markers.
Populationn 1Ar 2 (M ± SE) 7Ho 3 (M ± SE)uHe 4 (M ± SE)uFis 5 [CI] 6
DM1506.827 ± 0.6540.723 ± 0.0320.764 ± 0.0310.052 [0.005; 0.099]
BS_D135.818 ± 0.4830.692 ± 0.0690.713 ± 0.0550.023 [−0.093; 0.139]
BS_G275.146 ± 0.4440.741 ± 0.0450.701 ± 0.039−0.058 [−0.111; −0.005]
RS265.278 ± 0.4500.692 ± 0.0540.718 ± 0.0260.046 [−0.059; 0.151]
HOL504.376 ± 0.3960.669 ± 0.0440.645 ± 0.039−0.038 [−0.110; 0.034]
SIM505.307 ± 0.3730.729 ± 0.0330.682 ± 0.031−0.073 [−0.115; −0.031]
1n, number of individuals; 2 Ar, rarefied allele richness; 3 Ho, observed heterozygosity; 4 uHe, unbiased expected heterozygosity; 5 uFis, unbiased inbreeding coefficient; 6 CI, variation range of coefficient of uFis at a confidence interval of 95%; 7 M ± SE, mean value and standard error.
Table 4. Genetic differentiation between the DM and other studied cattle breeds.
Table 4. Genetic differentiation between the DM and other studied cattle breeds.
PopulationDMBS_DBS_GRSHOLSIM
DM00.023 10.053 10.036 10.113 10.058 1
BS_D0.049 200.047 10.049 10.135 10.093 1
BS_G0.130 20.054 200.081 10.134 10.116 1
RS0.066 20.043 20.150 200.105 10.084 1
HOL0.271 20.261 20.250 20.189 200.138 1
SIM0.114 20.195 20.210 20.132 20.234 20
1 FST values are presented above the diagonal; 2 Jost’s D values are presented below the diagonal.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Volkova, V.V.; Abdelmanova, A.S.; Deniskova, T.E.; Romanenkova, O.S.; Khozhokov, A.A.; Ozdemirov, A.A.; Sermyagin, A.A.; Zinovieva, N.A. Investigation of the Genetic Diversity of Dagestan Mountain Cattle Using STR-Markers. Diversity 2022, 14, 569. https://0-doi-org.brum.beds.ac.uk/10.3390/d14070569

AMA Style

Volkova VV, Abdelmanova AS, Deniskova TE, Romanenkova OS, Khozhokov AA, Ozdemirov AA, Sermyagin AA, Zinovieva NA. Investigation of the Genetic Diversity of Dagestan Mountain Cattle Using STR-Markers. Diversity. 2022; 14(7):569. https://0-doi-org.brum.beds.ac.uk/10.3390/d14070569

Chicago/Turabian Style

Volkova, Valeria V., Alexandra S. Abdelmanova, Tatiana E. Deniskova, Olga S. Romanenkova, Abdusalam A. Khozhokov, Alimsoltan A. Ozdemirov, Alexander A. Sermyagin, and Natalia A. Zinovieva. 2022. "Investigation of the Genetic Diversity of Dagestan Mountain Cattle Using STR-Markers" Diversity 14, no. 7: 569. https://0-doi-org.brum.beds.ac.uk/10.3390/d14070569

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