Next Article in Journal
Behavioral and Cognitive Performance Following Exposure to Second-Hand Smoke (SHS) from Tobacco Products Associated with Oxidative-Stress-Induced DNA Damage and Repair and Disruption of the Gut Microbiome
Next Article in Special Issue
Comparative Mitogenome Analysis of Two Native Apple Snail Species (Ampullariidae, Pomacea) from Peruvian Amazon
Previous Article in Journal
Genomic Complexity and Complex Chromosomal Rearrangements in Genetic Diagnosis: Two Illustrative Cases on Chromosome 7
Previous Article in Special Issue
Comparison of Mitochondrial Genome Sequences between Two Palaemon Species of the Family Palaemonidae (Decapoda: Caridea): Gene Rearrangement and Phylogenetic Implications
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 9
10.3390/genes14091701
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences

by
Olga Koshkina
1,
Tatiana Deniskova
1,*,
Arsen Dotsev
1,
Elisabeth Kunz
2,
Marina Selionova
3,
Ivica Medugorac
2 and
Natalia Zinovieva
1
1
L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy 60, Podolsk Municipal District, Moscow 142132, Russia
2
Population Genomics Group, Department of Veterinary Sciences, Ludwig-Maximilians-University Munich, 82152 Munich, Germany
3
Timiryazev Agricultural Academy, Russian State Agrarian University-Moscow, Timiryazevskaya Street, 41, Moscow 127550, Russia
*
Author to whom correspondence should be addressed.
Genes 2023, 14(9), 1701; https://doi.org/10.3390/genes14091701
Submission received: 2 August 2023 / Revised: 22 August 2023 / Accepted: 25 August 2023 / Published: 27 August 2023
(This article belongs to the Special Issue Mitochondrial DNA Replication and Transcription)
Browse Figures
Review Reports Versions Notes

Abstract

:
Eurasia is represented by all climatic zones and various environments. A unique breed variety of farm animals has been developed in Russia, whose territory covers a large area of the continent. A total of 69 local breeds and types of dairy, wool, and meat sheep (Ovis aries) are maintained here. However, the genetic diversity and maternal origin of these local breeds have not been comprehensively investigated. In this study, we describe the diversity and phylogeny of Russian sheep breeds inhabiting different geographical regions based on the analysis of complete sequences of mitochondrial genomes (mtDNA). Complete mtDNA sequences of the studied sheep were obtained using next-generation sequencing technology (NGS). All investigated geographical groups of sheep were characterized by high haplotype (Hd = 0.9992) and nucleotide diversity (π = 0.00378). Analysis of the AMOVA results showed that genetic diversity was majorly determined by within-population differences (77.87%). We identified 128 haplotypes in all studied sheep. Haplotypes belonged to the following haplogroups: B (64.8%), A (28.9%), C (5.5%), and D (0.8%). Haplogroup B was predominant in the western part of Russia. A high level of mtDNA polymorphism in the studied groups of local sheep indicates the presence of a significant reserve of unique genotypes in Russia, which is to be explored.

1. Introduction

Domestic sheep (O. aries) is one of the first domesticated, adaptable, and widespread livestock species [1]. The Asian mouflon (O. orientalis) is proposed as the ancestor of domestic sheep. Several studies reported that most of the functional genetic diversity in sheep came from mouflon during progressive domestication [2,3].
According to FAO statistics, the population number of domestic sheep worldwide was 1.284 billion in 2021, and 1384 different breeds of sheep have been identified in 2022 [4]. The center of sheep domestication is the Fertile Crescent territory in the Middle East [3]. This was indirectly confirmed by the higher nucleotide diversity of mtDNA sequences of modern sheep in the Levant region [5,6].
The spreading of sheep from the domestication center, mediated by a variety of human activities, took place during the transition from the Mesolithic culture to the Neolithic era (about 8000–9000 years ago) [1]. Chessa et al. [1] used the integration of endogenous retroviruses in the sheep genome as genetic markers and postulated two major waves of sheep migration toward northwest Europe. The first wave began with the spreading of “primitive” sheep from the Fertile Crescent. The second wave of improved sheep migration occurred about 5000 years ago [7].
After approximately 1000 years of domestication, domestic animals and plants expanded and replaced the foraging lifestyle in larger parts of Southwest Asia [8]. Further, they subsequently spread to Europe along two routes: through the Mediterranean Sea and through the Danube Valley. Phylogenetic analysis of lentiviruses of small ruminants supported evidence of these events [9]. The Mediterranean route passed from Cyprus through the Balkan and Apennine peninsulas to Corsica and Sardinia. Later the settlement of Northern Italy and Southern France took place. Then, along the Danube route, sheep spread to Central and Northern Europe [10]. Analyzing mtDNA of Copper Age European sheep, Olivieri et al. [7] concluded that sheep ended up in the Alps a little more than 5000 years ago. In the Iberian Peninsula, the Neolithic settlements were established in Spain (7700–7600 years ago) [11], and then in Portugal (7400–7300 years ago) [5]. However, according to another hypothesis, the introduction of the Ovis genus had occurred repeatedly from Africa and Central Europe to the Iberian Peninsula. At the same time, there was a settlement of Northern Europe through the center of domestication through the Caucasus and Russia [12,13].
According to references, the colonization of Asia began from the Middle East to the Mongolian plateau and the Indian subcontinent, from where sheep migrated to the north and southwest of China [14,15]. Sheep distribution across Asia could have happened through several routes as it had happened in Europe [16].
The territory of Russia covering a large part of Eurasia is in the arctic, subarctic, temperate, and partly in the subtropical climatic zones. Thus, the country is characterized by a wide variety of climatic conditions. A total of 69 different sheep breeds have been registered on its territory [17], which produce agricultural products in various climatic conditions. The diversity of environments and long-term adaptation to these resulted in the great genetic diversity of Russian breeds.
The genetic diversity of farm animals plays an important role in the more complete performance of economically significant, fitness-related traits [18]. The conservation of local breeds of livestock is important because local breeds meet human needs in the face of impending climate change scenarios [19] and present a valuable pool of allelic variants for future unknown usage [20]. Moreover, the study of sheep genetics is relevant for understanding aspects of selection, breed diversity, and domestication history as well as for breed sustainable maintenance and conservation. Using genetic markers to address breeding problems in animal husbandry helps to preserve the existing gene pool of the breeds and to improve the economic profitability of the industry.
Mitochondrial DNA (mtDNA) is inherited only through the maternal lineages and has great variability within the same species. In this regard, mtDNA polymorphism analysis is widely used to investigate phylogenetic links and genetic structure of various livestock species [21,22,23,24,25]. Numerous studies including tracing the maternal lineage, and establishing phylogenetic relationships, structure, and diversity of the population have been carried out based on the D-loop mtDNA region and the CytB gene [26,27].
Based on the analysis of sheep mtDNA control region polymorphism, seven haplogroups were identified. Haplogroups F and G disappeared, and five others (A, B, C, D, E) are present in modern breeds [5,28]. The first identified haplogroups were A and B, which corresponded to Asian and European types of origin, respectively [28,29]. Later, the third recognized phylogenetic branch, haplogroup C, was discovered [30,31]. Haplogroup C is characterized by a larger genetic diversity than A and B [31]. This haplogroup with low frequency was detected in Portuguese native sheep [32], in individuals from the Caucasus, the Middle East, and Asia. Then, the fourth maternal lineage, haplogroup D, which is the closest to the mouflon (O. gmelini anatolica Valenciennes), was identified [12,33,34]. Using D-loop and CytB allowed for the detection of the fifth haplogroup E [5]. The latest discovered haplogroups are the rarest ones, which were identified only in animals from Turkey, the Caucasus [5,12,34], Europe [35,36], Tibet [37], and Iran [38]. These haplogroups were formed presumably from 5 to 35 thousand years ago [39].
Thus, a comprehensive study of mtDNA polymorphism allows unlocking the evolution events, which had occurred during the migration of domestic sheep populations throughout the history of their domestication. However, most sheep population studies are based on using the D-loop [40] or the CytB gene [41] sequences. The polymorphism of other mitochondrial genes has not been studied so precisely. The study of complete mitochondrial genomes will expand knowledge about the role of mtDNA in the livelihood of sheep.
In this research, for the first time, an investigation of the genetic diversity and structure of Russian sheep populations based on the polymorphism of complete mitochondrial genomes is presented.

2. Materials and Methods

2.1. Sample Collection

Tissue specimens (ear fragments) of Russian sheep breeds were used as biological material for this research. Sheep samples were provided by the Bioresource Collection of Farm Animals of the L.K. Ernst Federal Science Center for Animal Husbandry (deposited to the collection in the period from 2003 to 2018).
Based on pedigree records, the unrelated individuals were chosen. Additionally, all samples were genotyped using STR markers. The matrix of STR-genotypes was analyzed in ML-Relate software [https://www.montana.edu/kalinowski/software/mL-relate/index.html (accessed on 7 December 2022) [42]. Only unrelated animals were selected for mtDNA sequencing.
The final studied sample included 135 sheep inhabiting six territorial districts of Russia, corresponding to diverse geographical zones (Table 1).

2.2. DNA Extraction, PCR Amplification, and Sequencing

DNA was isolated from ear fragments using commercial «DNA-Extran-2» kits (CJSC Syntol, Moscow, Russia) in accordance with the manufacturer’s guidelines. An integrity of DNA was checked by gel electrophoresis in 1% agarose gel via the colorimetric gel documentation system Uvitec FireReader V10 imaging System (Cleaver Scientific Ltd, Rugby, Warwickshire, UK). The OD260/OD280 ratio and concentrations of DNA solutions were measured using a NanoDrop 8000 (Thermo Fisher Scientific, Waltham, MA, USA) and a Qubit 4.0 fluorimeter (Thermo Fisher Scientific, Waltham, MA, USA), respectively. The isolated DNA samples were preserved at −20 °C until further laboratory manipulations.
For the next generation sequencing of complete mitogenomes, six mtDNA fragments with overlapping regions were amplified. Primer sequences for the amplification of target overlapping mtDNA fragments were selected based on the reference sequence of complete mitochondrial genome of sheep, which were presented in the National Center for Biotechnology Information NCBI (GenBank accession number NC_001941.1, accessed on 10 January 2023) according to O. aries ARS-UI_Ramb_v2.0 sheep genome version using the online resource Basic Local Alignment Search Tool (BLAST) (https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 10 January 2023)) [43]. The nucleotide sequences of the used primers are shown in Table 2.
The final volume for PCR was 25 µL and contained 10 µL PCR-buffer (2.5½ HF Reaction buffer), 10.25 µL H2O, 2.5 µL dNTPs (1 mmol/L), 1 µL primer mix, 0.25 µL SmartTaq HF-FuZZ DNA polymerase (Dialat Ltd., Moscow, Russia), and 1 µL of DNA template.
Quality of PCR amplification was checked by gel electrophoresis in 1% agarose gel using the colorimetric gel documentation system Uvitec FireReader V10 imaging System (Cleaver Scientific Ltd, Rugby, Warwickshire, UK).
PCR products were purified using the «Cleanup Mini kit for the isolation and purification of DNA from agarose gel and reaction mixtures» (JSC Eurogene, Moscow, Russia) and were used to prepare sequencing libraries with NEBNext Ultra II DNA Library Prep Kit for Illumina (New England Biolabs Inc., Ipswich, UK) following the manufacturer’s recommendations. To check the quality, sequencing libraries were submitted to CJSC Syntol (Russia).
Sequencing of the prepared libraries was performed by 300-bp paired-end procedure on MiSeq (Illumina, Inc., San Diego, CA, USA).

2.3. Bioinformatics Data Processing

MtDNA sequences were aligned and edited using the MUSCLE algorithm [44] as implemented in MEGA 7.0.26 software (https://www.megasoftware.net/history (accessed on 12 May 2023)) [45].
Measures of genetic diversity parameters including number of polymorphic sites (S), average number of nucleotide differences (K), number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (π), and standard error mean (±SEM) were calculated in DNASP 6.12.01 software (https://www.bioinformaticshome.com/db/tool/DnaSP (accessed on 7 June 2023)) [46].
Haplogroups for Russian sheep were identified using the MitoToolPy software (http://www.mitotool.org/ (accessed on 15 June 2023)) [47].
The best evolution models were determined using the PartitionFinder 2 program [48] using the Akaike information corrected criterion (AICc) (AICc) [49]. The evolutionary model HKY + I + G was the most optimal.
The phylogenetic relationship was established by the Markov chain Monte Carlo search running with four chains for 10,000,000 generations, with trees being sampled every 500 generations (the first 25% of the trees were discarded as ‘burnin’) in the MrBayes 3.2.7 program (http://nbisweden.github.io/MrBayes/ (accessed on 20 June 2023)) [50] and further visualized in FigTree 1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 20 June 2023)) [51].
To analyze evolutionary relationships, a median-joining network was constructed in PopART 1.7 software (https://popart.maths.otago.ac.nz/ (accessed on 30 June 2023)) [52].
Wright’s F-statistics [53] and analysis of molecular variance (AMOVA) [54], implemented in Arlequin v 3.5 (http://cmpg.unibe.ch/software/arlequin35/Arl35Downloads.html (accessed on 7 July 2023)) [55], were used to assess population differentiation and to evaluate genetic variability among and within a population.
The demographic history of populations was assessed based on pairwise mismatch distribution analyses using the DnaSP 6.12.01 software [46], and tests for neutrality (Tajima D test, Fu’s Fs test) were performed in DnaSP 6.12.01 [46] and Arlequin v 3.5 [55].

3. Results

3.1. Genetic Diversity and Haplotype Network

To determine single nucleotide polymorphisms (SNPs) in the complete mitochondrial genomes of sheep, genetic diversity indices were calculated (Table 3).
A total of 804 polymorphic sites were identified in 135 studied animals. The largest number of polymorphic sites (S = 468) was determined in sheep of the North Caucasian region, and the smallest number (S = 49) was identified in the Volga region. The lowest values of nucleotide diversity and the smallest average number of nucleotide differences were recorded in sheep from the Volga region (π = 0.00123 ± 0.00017; K = 20.400). Sheep from the Far Eastern region were characterized by the highest values of these indicators (π = 0.00589 ± 0.00136; K = 97.893).
The results of Tajima’s D tests were not significant in all sheep groups which may indicate population stability. However, high negative Fu’s Fs values were found in the groups from the Southern and North Caucasian regions. Most likely, this indicates foreign gene flow due to spatial expansion. Perhaps this is due to the presence of representatives of rare haplogroups (C and D), which are characterized by high mtDNA polymorphism, in these geographical regions.
In all studied sheep, a total of 128 haplotypes were identified including 7 haplotypes in FAE, 5 in VOL, 46 in NCN, 17 in SIB, 38 in SOU, and 19 in CEN, respectively. The distribution of haplotypes of complete mitochondrial genomes of sheep is presented in Figure 1 as a median-joining network. There are seven haplotypes, each present in two animals: the haplotype Hap_42 was detected in two samples from the North Caucasian region, Hap_103 and Hap_105 in two samples from the Central, Hap_115 in two samples from Siberian, Hap_89 in two samples from Far Eastern region, Hap_6 in one animal from the North Caucasus and one from the Central region, and Hap_109 in one animal from the Central and one from the Siberian region. The remaining haplotypes were specific to the studied sample. The investigated animals were characterized by high haplotype diversity (Hd = 0.9992 ± 0.01). At the same time, the greatest haplotype diversity was determined in sheep from the Southern, Siberian and Volga regions (Hd = 1.000 ± 0.006; 1.000 ± 0.020 and 1.000 ± 0.126, respectively).
We checked all identified haplotypes with those published worldwide by comparing nucleotide sequences using the BLAST tool. We did not find the 100% nucleotide similarity (Table S1). In this regard, we may assume that the haplotypes found in Russian sheep breeds have not been described previously.

3.2. Haplogroup Assignment

In the studied groups of sheep four haplogroups were determined (A, B, C, and D), using the MitoToolPy program.
Analysis of median-joining network (Figure 1) and the Bayesian tree (Figure 2) revealed that the most frequent among the studied groups of sheep was haplogroup B, which included 83 haplotypes (64.8%). Haplogroup A was accounted for by 37 haplotypes (28.9%), haplogroup C was represented by seven haplotypes (5.5%), and only one haplotype was assigned to haplogroup D (0.8%). Haplogroups A and B were present in all regional groups of sheep. Haplogroup C was represented by the haplotypes from the North Caucasus and Far Eastern regions, while haplogroup D was found in a single animal from the southern regional group.

3.3. Genetic Differentiation between Russian Sheep from Diverse Geographical Regions

To evaluate genetic differentiation of sheep from different regional groups, pairwise FST values were calculated (Table 4).
The greatest genetic remoteness was showed between Far Eastern and Central geographical groups (FST = 0.43960) and between Volga and Far East geographical groups (FST = 0.39243), respectively (Table 4). Obviously, this is due to the territorial remoteness of breeding these groups of sheep. The pairs from the Central and Volga regions (FST = 0.00818), as well as the Southern and North Caucasian regions (FST = 0.00937) were characterized by the greatest genetic similarity.
Analysis of molecular variance (AMOVA) (Table 5) confirmed a significant genetic differentiation within the studied groups of sheep (77.87%), while differences among studied populations accounted for 22.13% of the genetic variability.

3.4. Population Dynamics and Mismatch Distribution of Russian Sheep from Different Regions

The mismatch distribution (distribution of the number of pairwise differences between sequences) for all groups of sheep was multimodal (Figure 3). All studied groups showed at least two peaks, which suggests that several independent events of gene influx from outside occurred in the past. This may be a consequence of the artificial movement of domestic sheep between different populations or of the high genetic diversity within and between haplogroups in originally Neolithic stock.

4. Discussion

The study of genetic diversity of sheep provides a better understanding of aspects of breeding and domestication processes. MtDNA is distinguished by its uniqueness and is one of the most informative genetic markers for phylogenetic studies, because mtDNA has a maternal inheritance pattern and allows determining introgression only in females. Here, we investigated complete sequences of mitochondrial genomes of Russian sheep originated from different geographic regions. The phylogeny of mitochondrial genome of the genus Ovis is formed from five major haplogroups (A, B, C, D, and E) [6].
Based on analysis of polymorphism of complete mtDNA nucleotide sequences, we identified four sheep haplogroups, including A, B, C, and D, in studied sheep populations.
According to previous studies [6,28,29,56], haplogroup B is typical for sheep of European type of origin. This is confirmed in our study as well, because the most frequent among the studied groups of sheep was haplogroup B, which was represented by 89 animals from all regions. The second largest group (n = 38) was haplogroup A, which was also found in all regions, but most of the animals were bred in the Siberian and Far Eastern geographical regions.
As in previous studies [6,57,58], haplogroup A prevailed in Asian sheep populations. Seven animals were assigned to haplogroup C. They were mainly from the North Caucasian region, but two representatives were from the Far Eastern geographical region as well. According to the review by Machová K. et al., the haplogroup C was found in Portugal, the Caucasus, the Middle East, and Asia [59]. In addition, Lv et al. in their study provides a map of haplogroup distribution in the world [14]. According to it, haplogroup C is found mainly in the Middle East, the Caspian Sea region, Northern China, and the Mongolian Plateau. This corresponds to the data obtained in our study. Haplogroups D and E were the rarest [14] as well as in our study only one animal from the Southern ecoregion was belonged to haplogroup D. Taylor et al. [58] previously described a similar distribution of sheep haplogroups. The presence of Haplogroup A and B in all regional groups may indirectly indicate on several ways of settlement from the main center of domestication.
Haplotype (Hd = 0.999 ± 0.010) and nucleotide (π = 0.00378 ± 0.00027) diversity of Russian sheep were comparable to the values obtained earlier in African sheep (Hd = 0.993 ± 0.002, π = 0.00254 ± 0.012) [60] and Indian sheep (Hd = 0.8502 ± 0.024, π = 0.00373 ± 0.012) [57].
The results of the AMOVA determined higher genetic differentiation within the studied groups of sheep (77.87%), rather than among populations (23.13%). Our data were corresponded to the results, which were obtained in Indian sheep (variability among sheep groups is accounted for 11.1%) [57]. Such a pattern of genetic differentiation may indicate on a certain commonality of the maternal origin of the studied groups of sheep from different geographical regions. As already mentioned by several authors e.g., [61], there is hardly any significant breed structure among the maternal lineages of the modern sheep groups, making it difficult to address their relationship to ancient sheep based on mtDNA alone.
Tajima’s D neutrality test results were close to zero and statistically insignificant in all sheep populations, suggesting that the evolution of maternal lineages was rather random. However, high negative Fu’s Fs (−10.572 **; −8.987 **) values were found in several groups and suggest population expansion. Similar results were obtained in Indian (Tajima’s D = −2.1; −0.97; Fu’s Fs = −26.12; −25.66) [57] and African sheep (Tajima’s D = −0.921; 3.715; Fu’s Fs = −24.003; −10.426) [61].
Thus, based on the results, we may summarize that there is a significant reserve of unique genotypes in Russian sheep populations. In future we would continue our investigations of genetic diversity of these populations for better understanding of their developmental history, and genetic relationships.

5. Conclusions

On the territory of Russia, there is a large population of domestic sheep with specific genotypes, the genetic potential of which has not been fully explored yet. In the present study, we analyzed the polymorphism of complete mtDNA nucleotide sequences of 135 Russian sheep from different habitats. A large number of haplotypes (H = 128) was identified, of which 94.5% (H = 121) were unique. High nucleotide and haplotype diversity was revealed that corresponded to high variability of mtDNA. At the same time, the median-joining network of haplotypes formed joined clusters among the populations, which may indicate on a common maternal origin of the studied sheep. The results obtained in this study will be useful for improving strategies for the conservation of sheep genetic resources. Continuing research on large samples will provide a more complete pattern of genetic diversity and demographic history of sheep inhabiting different geographical areas in Eurasia.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14091701/s1, Table S1: The results of comparative analysis of identified haplotypes with those published worldwide.

Author Contributions

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

Funding

This research was financially supported by the Ministry of Science and Higher Education of the Russian Federation, grant no. 075-15-2021-1037 (internal no. 15.БPK.21.0001) (sample collection and NGS of sheep breeds), and theme no. FGGN-2022-0002 (NGS of sheep breeds from Central region). Russian Science Foundation (RSF) within project no. 21-46-00001 financially supported NGS of Karachaev and Edilbai sheep breeds.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Commission on the Ethics of Animal Experiments of the L.K. Ernst Federal Science Center for Animal Husbandry (the protocol No. 2 (28 April 2022)).

Informed Consent Statement

Not applicable.

Data Availability Statement

The complete nucleotide sequences of mitochondrial genomes of the studied sheep were deposited as population sets according to their geographical origin to NCBI GenBank (accessed on 16 August 2023) under accession numbers OR459640-OR459774.

Acknowledgments

The equipment of the Center for Biological Resources and Bioengineering of Agricultural Animals (L.K. Ernst Federal Science Center for Animal Husbandry) was used in the research.

Conflicts of Interest

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

References

  1. Chessa, B.; Pereira, F.; Arnaud, F.; Amorim, A.; Goyache, F.; Mainland, I.; Kao, R.R.; Pemberton, J.M.; Beraldi, D.; Stear, M.J.; et al. Revealing the history of sheep domestication using retrovirus integrations. Science 2009, 324, 532–536. [Google Scholar] [CrossRef] [PubMed]
  2. Luigi-Sierra, M.G.; Mármol-Sánchez, E.; Amills, M. Comparing the diversity of the casein genes in the Asian mouflon and domestic sheep. Anim. Genet. 2020, 51, 470–475. [Google Scholar] [CrossRef] [PubMed]
  3. Zeder, M.A. Domestication and early agriculture in the Mediterranean Basin: Origins, diffusion, and impact. Proc. Natl. Acad. Sci. USA 2008, 105, 11597–11604. [Google Scholar] [CrossRef]
  4. Food and Agriculture Organization of the United Nations (FAO). Available online: https://www.fao.org/livestock-systems/global-distributions/sheep/en/. (accessed on 21 August 2023).
  5. Chen, Z.H.; Xu, Y.X.; Xie, X.L.; Wang, D.F.; Aguilar-Gómez, D.; Liu, G.J.; Li, X.; Esmailizadeh, A.; Rezaei, V.; Kantanen, J.; et al. Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep. Commun. Biol. 2021, 4, 1307. [Google Scholar] [CrossRef] [PubMed]
  6. Meadows, J.R.; Cemal, I.; Karaca, O.; Gootwine, E.; Kijas, J.W. Five ovine mitochondrial lineages identified from sheep breeds of the near East. Genetics 2007, 175, 1371–1379. [Google Scholar] [CrossRef] [PubMed]
  7. Olivieri, C.; Ermini, L.; Rizzi, E.; Corti, G.; Luciani, S.; Marota, I.; De Bellis, G.; Rollo, F. Phylogenetic position of a copper age sheep (Ovis aries) mitochondrial DNA. PLoS ONE 2012, 7, e33792. [Google Scholar] [CrossRef]
  8. Peters, J.; von den Driesch, A.; Helmer, D. The upper Euphrates-Tigris Basin: Cradle of agro-pastoralism? In The First Steps of Animal Domestication. New Archaeological Approaches; Vigne, J.D., Peters, J., Helmer, D., Eds.; Oxbow Books: Oxford, UK, 2005; pp. 96–124. [Google Scholar]
  9. Molaee, V.; Bazzucchi, M.; De Mia, G.M.; Otarod, V.; Abdollahi, D.; Rosati, S.; Lühken, G. Phylogenetic analysis of small ruminant lentiviruses in Germany and Iran suggests their expansion with domestic sheep. Sci. Rep. 2020, 10, 2243. [Google Scholar] [CrossRef]
  10. Larsson, M.N.A.; Miranda, P.M.; Pan, L.; Başak Vural, K.; Kaptan, D.; Rodrigues Soares, A.E.; Kivikero, H.; Kantanen, J.; Somel, M.; Özer, F.; et al. Ancient sheep genomes reveal four Millennia of North European short-tailed sheep in the Baltic Sea region. bioRxiv 2023, 544912. [Google Scholar] [CrossRef]
  11. Zeder, M.A. Domestication as a model system for the extended evolutionary synthesis. Interface Focus 2017, 7, 20160133. [Google Scholar] [CrossRef]
  12. Tapio, M.; Marzanov, N.; Ozerov, M.; Cinkulov, M.; Gonzarenko, G.; Kiselyova, T.; Murawski, M.; Viinalass, H.; Kantanen, J. Sheep mitochondrial DNA variation in European, Caucasian, and Central Asian areas. Mol. Biol. Evol. 2006, 23, 1776–1783. [Google Scholar] [CrossRef]
  13. Niemi, M.; Bläuer, A.; Iso-Touru, T.; Nyström, V.; Harjula, J.; Taavitsainen, J.P.; Storå, J.; Lidén, K.; Kantanen, J. Mitochondrial DNA and Y-chromosomal diversity in ancient populations of domestic sheep (Ovis aries) in Finland: Comparison with contemporary sheep breeds. Genet. Sel. Evol. 2013, 45, 2. [Google Scholar] [CrossRef] [PubMed]
  14. Lv, F.H.; Peng, W.F.; Yang, J.; Zhao, Y.X.; Li, W.R.; Liu, M.J.; Ma, Y.H.; Zhao, Q.J.; Yang, G.L.; Wang, F.; et al. Mitogenomic meta-analysis identifies two phases of migration in the history of Eastern Eurasian sheep. Mol. Biol. Evol. 2015, 32, 2515–2533. [Google Scholar] [CrossRef] [PubMed]
  15. Ganbold, O.; Lee, S.H.; Paek, W.K.; Munkhbayar, M.; Seo, D.; Manjula, P.; Khujuu, T.; Purevee, E.; Lee, J.H. Mitochondrial DNA variation and phylogeography of native Mongolian goats. Asian Australas J. Anim. Sci. 2020, 33, 902–912. [Google Scholar] [CrossRef] [PubMed]
  16. Singh, S.; Kumar, S., Jr.; Kolte, A.P.; Kumar, S. Extensive variation and sub-structuring in lineage A mtDNA in Indian sheep: Genetic evidence for domestication of sheep in India. PLoS ONE 2013, 8, e77858. [Google Scholar] [CrossRef]
  17. Ministry of Agriculture of the Russian Federation. Open Data. Available online: http://opendata.mcx.ru/opendata. (accessed on 17 April 2023).
  18. Zhong, Y.; Tang, Z.; Huang, L.; Wang, D.; Lu, Z. Genetic diversity of Procambarus clarkii populations based on mitochondrial DNA and microsatellite markers in different areas of Guangxi, China. Mitochondrial DNA Part A 2020, 31, 48–56. [Google Scholar] [CrossRef] [PubMed]
  19. Scherf, B.D. World Watch List for Domestic Animal Diversity, 3rd ed.; FAO: Rome, Italy, 2000; pp. 7–13. [Google Scholar]
  20. Ramljak, J.; Bunevski, G.; Bytyqi, H.; Marković, B.; Brka, M.; Ivanković, A.; Kume, K.; Stojanović, S.; Nikolov, V.; Simčič, M.; et al. Conservation of a domestic metapopulation structured into related and partly admixed strains. Mol. Ecol. 2018, 27, 1633–1650. [Google Scholar] [CrossRef]
  21. Mishra, A.K.; Ahlawat, S.; Sharma, R.; Arora, R.; Singh, S.; Jain, A. Assessment of genetic diversity of the fat-tailed Dumba sheep of India by mitochondrial and microsatellite markers. Anim. Biotechnol. 2023; advance online publication. [Google Scholar]
  22. Salim, B.; Alasmari, S.; Mohamed, N.S.; Ahmed, M.A.; Nakao, R.; Hanotte, O. Genetic variation and demographic history of Sudan desert sheep reveal two diversified lineages. BMC Genom. 2023, 24, 118. [Google Scholar] [CrossRef]
  23. Al-Araimi, N.A.; Al-Atiyat, R.M.; Gaafar, O.M.; Vasconcelos, R.; Luzuriaga-Neira, A.; Eisa, M.O.; Amir, N.; Benaissa, M.H.; Alfaris, A.A.; Aljumaah, R.S.; et al. Maternal genetic diversity and phylogeography of native Arabian goats. Livest. Sci. 2017, 206, 88–94. [Google Scholar] [CrossRef]
  24. Ibrahim, A.; Baliarti, E.; Budisatria, I.G.S.; Artama, W.T.; Widayanti, R.; Maharani, D.; Tavares, L.; Margawati, E.T. Genetic diversity and relationship among Indonesian local sheep breeds on Java Island based on mitochondrial cytochrome b gene sequences. J. Genet. Eng. Biotechnol. 2023, 21, 34. [Google Scholar] [CrossRef]
  25. Dymova, M.A.; Zadorozhny, A.V.; Mishukova, O.V.; Khrapov, E.A.; Druzhkova, A.S.; Trifonov, V.A.; Kichigin, I.G.; Tishkin, A.A.; Grushin, S.P.; Filipenko, M.L. Mitochondrial DNA analysis of ancient sheep from Altai. Anim. Genet. 2017, 48, 615–618. [Google Scholar] [CrossRef]
  26. Kamalakkannan, R.; Kumar, S.; Bhavana, K.; Prabhu, V.R.; Machado, C.B.; Singha, H.S.; Sureshgopi, D.; Vijay, V.; Nagarajan, M. Evidence for independent domestication of sheep mtDNA lineage A in India and introduction of lineage B through Arabian sea route. Sci. Rep. 2021, 11, 19733. [Google Scholar] [CrossRef] [PubMed]
  27. Mukhametzharova, I.; Islamov, Y.; Shauyenov, S.; Ibrayev, D.; Atavliyeva, S.; Tarlykov, P. Genetic characterization of Kazakh native sheep breeds using mitochondrial DNA. J. Biol. Sci. 2018, 18, 341–348. [Google Scholar] [CrossRef]
  28. Wood, N.J.; Phua, S.H. Variation in the control region sequence of the sheep mitochondrial genome. Anim. Genet. 1996, 27, 25–33. [Google Scholar] [CrossRef]
  29. Hiendleder, S.; Mainz, K.; Plante, Y.; Lewalski, H. Analysis of mitochondrial DNA indicates that domestic sheep are derived from two different ancestral maternal sources: No evidence for contributions from urial and argali sheep. J. Hered. 1998, 89, 113–120. [Google Scholar] [CrossRef] [PubMed]
  30. Guo, J.; Du, L.X.; Ma, Y.H.; Guan, W.J.; Li, H.B.; Zhao, Q.J.; Li, X.; Rao, S.Q. A novel maternal lineage revealed in sheep (Ovis aries). Anim. Genet. 2005, 36, 331–336. [Google Scholar] [CrossRef] [PubMed]
  31. Pedrosa, S.; Uzun, M.; Arranz, J.J.; Gutiérrez-Gil, B.; San Primitivo, F.; Bayón, Y. Evidence of three maternal lineages in Near Eastern sheep supporting multiple domestication events. Proc. Biol. Sci. 2005, 272, 2211–2217. [Google Scholar] [CrossRef] [PubMed]
  32. Pereira, F.; Davis, S.J.; Pereira, L.; McEvoy, B.; Bradley, D.G.; Amorim, A. Genetic signatures of a Mediterranean influence in Iberian Peninsula sheep husbandry. Mol. Biol. Evol. 2006, 23, 1420–1426. [Google Scholar] [CrossRef] [PubMed]
  33. Sanna, D.; Barbato, M.; Hadjisterkotis, E.; Cossu, P.; Decandia, L.; Trova, S.; Pirastru, M.; Leoni, G.G.; Naitana, S.; Francalacci, P.; et al. The first mitogenome of the Cyprus Mouflon (Ovis gmelini ophion): New insights into the phylogeny of the genus Ovis. PLoS ONE 2015, 10, e0144257. [Google Scholar] [CrossRef]
  34. Demirci, S.; Koban Baştanlar, E.; Dağtaş, N.D.; Pişkin, E.; Engin, A.; Ozer, F.; Yüncü, E.; Doğan, S.A.; Togan, I. Mitochondrial DNA diversity of modern, ancient and wild sheep (Ovis gmelinii anatolica) from Turkey: New insights on the evolutionary history of sheep. PLoS ONE 2013, 8, e81952. [Google Scholar] [CrossRef]
  35. Mariotti, M.; Valentini, A.; Marsan, P.A.; Pariset, L. Mitochondrial DNA of seven Italian sheep breeds shows faint signatures of domestication and suggests recent breed formation. Mitochondrial DNA 2013, 24, 577–583. [Google Scholar] [CrossRef]
  36. Gáspárdy, A.; Berger, B.; Zabavnik-Piano, J.; Kovács, E.; Annus, K.; Zenke, P.; Sáfár, L.; Maróti-Agóts, Á. Comparison of mtDNA control region among descendant breeds of the extinct Zaupel sheep revealed haplogroup C and D in Central Europe. Vet. Med. Sci. 2021, 7, 2330–2338. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, J.; Ding, X.; Zeng, Y.; Yue, Y.; Guo, X.; Guo, T.; Chu, M.; Wang, F.; Han, J.; Feng, R.; et al. Genetic diversity and phylogenetic evolution of Tibetan sheep based on mtDNA D-Loop sequences. PLoS ONE 2016, 11, e0159308. [Google Scholar]
  38. Rafia, P.; Tarang, A. Sequence variations of mitochondrial DNA displacement-loop in Iranian indigenous sheep breeds. J. Appl. Anim. Sci. 2016, 6, 363–368. [Google Scholar]
  39. Rezaei, H.R.; Naderi, S.; Chintauan-Marquier, I.C.; Jordan, S.; Taberlet, P.; Virk, A.T.; Naghash, H.R.; Rioux, D.; Kaboli, M.; Luikart, G.; et al. Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae). Mol. Phylogenet. Evol. 2010, 54, 315–326. [Google Scholar] [CrossRef] [PubMed]
  40. Barbato, M.; Masseti, M.; Pirastru, M.; Columbano, N.; Scali, M.; Vignani, R.; Mereu, P. Islands as time capsules for genetic diversity conservation: The case of the Giglio Island mouflon. Diversity 2022, 14, 609. [Google Scholar] [CrossRef]
  41. Liu, J.; Lu, Z.; Yuan, C.; Wang, F.; Yang, B. Phylogeography and phylogenetic evolution in Tibetan sheep based on mt-CyB sequences. Animals 2020, 10, 1177. [Google Scholar] [CrossRef]
  42. Kalinowski, S.T.; Wagner, A.P.; Taper, M.L. ML-Relate: A computer program for maximum likelihood estimation of relatedness and relationship. Mol. Ecol. Notes. 2006, 6, 576–579. [Google Scholar] [CrossRef]
  43. Basic Local Alignment Search Tool (BLAST). Available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 10 January 2023).
  44. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef]
  45. 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]
  46. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  47. Peng, M.S.; Fan, L.; Shi, N.N.; Ning, T.; Yao, Y.G.; Murphy, R.W.; Wang, W.; Zhang, Y.P. DomeTree: A canonical toolkit for mitochondrial DNA analyses in domesticated animals. Mol. Ecol. Resour. 2015, 15, 1238–1242. [Google Scholar] [CrossRef] [PubMed]
  48. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. PartitionFinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 2017, 34, 772–773. [Google Scholar] [CrossRef] [PubMed]
  49. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  50. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed]
  51. Molecular Evolution, Phylogenetics and Epidemiology. Available online: http://tree.bio.ed.ac.uk/soft-ware/figtree. (accessed on 1 April 2023).
  52. Leigh, J.W.; Bryant, D. Popart: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  53. Wright, S. Evolution in Mendelian populations. Genetics 1931, 16, 97–159. [Google Scholar] [CrossRef]
  54. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 1992, 131, 479–491. [Google Scholar] [CrossRef]
  55. Excoffier, L.; Lischer, H.E.L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour 2010, 10, 564–567. [Google Scholar] [CrossRef]
  56. Resende, A.; Gonçalves, J.; Muigai, A.W.; Pereira, F. Mitochondrial DNA variation of domestic sheep (Ovis aries) in Kenya. Anim. Genet. 2016, 47, 377–381. [Google Scholar] [CrossRef]
  57. Sharma, R.; Ahlawat, S.; Sharma, H.; Sharma, P.; Panchal, P.; Arora, R.; Tantia, M.S. Microsatellite and mitochondrial DNA analyses unveil the genetic structure of native sheep breeds from three major agro-ecological regions of India. Sci. Rep. 2020, 10, 20422. [Google Scholar] [CrossRef]
  58. Taylor, W.T.T.; Pruvost, M.; Posth, C.; Rendu, W.; Krajcarz, M.T.; Abdykanova, A.; Brancaleoni, G.; Spengler, R.; Hermes, T.; Schiavinato, S.; et al. Evidence for early dispersal of domestic sheep into Central Asia. Nat. Hum. Behav. 2021, 5, 1169–1179. [Google Scholar] [CrossRef]
  59. Machová, K.; Málková, A.; Vostrý, L. Sheep post-domestication expansion in the context of mitochondrial and Y chromosome haplogroups and haplotypes. Genes 2022, 13, 613. [Google Scholar] [CrossRef] [PubMed]
  60. Wanjala, G.; Bagi, Z.; Kusza, S. Meta-analysis of mitochondrial DNA control region diversity to shed light on phylogenetic relationship and demographic history of African sheep (Ovis aries) Breeds. Biology 2021, 10, 762. [Google Scholar] [CrossRef] [PubMed]
  61. Cai, D.; Tang, Z.; Yu, H.; Han, L.; Ren, X.; Zhao, X.; Zhu, H.; Zhou, H. Early history of Chinese domestic sheep indicated by ancient DNA analysis of Bronze Age Individuals. J. Archaeol. Sci. 2011, 38, 896–902. [Google Scholar] [CrossRef]
Genes 14 01701 g001
Figure 1. Median-joining network characterizing haplotype connections identified in Russian local sheep breeds based on analysis of the nucleotide sequence of complete mitogenomes. Note: the studied geographical groups of sheep: FAE: Far Eastern; VOL: Volga; NCN: North Caucasian; SIB: Siberian; SOU: Southern; CEN: Central. Colors within each haplogroup ring correspond to different geographical regions. Hap89, Hap42, Hap109, Hap115, Hap 103, and Hap 105 are shared haplotypes.
Figure 1. Median-joining network characterizing haplotype connections identified in Russian local sheep breeds based on analysis of the nucleotide sequence of complete mitogenomes. Note: the studied geographical groups of sheep: FAE: Far Eastern; VOL: Volga; NCN: North Caucasian; SIB: Siberian; SOU: Southern; CEN: Central. Colors within each haplogroup ring correspond to different geographical regions. Hap89, Hap42, Hap109, Hap115, Hap 103, and Hap 105 are shared haplotypes.
Genes 14 01701 g001
Genes 14 01701 g002
Figure 2. Bayesian phylogenetic tree reflecting genetic connections of Russian sheep breeds based on analysis of nucleotide sequences of complete mitochondrial genomes. Note: the studied geographical groups of sheep: FAE: Far Eastern; VOL: Volga; NCN: North Caucasian; SIB: Siberian; SOU: Southern; CEN: Central. Colors of branches correspond to different haplogroups (A, B, C, and D).
Figure 2. Bayesian phylogenetic tree reflecting genetic connections of Russian sheep breeds based on analysis of nucleotide sequences of complete mitochondrial genomes. Note: the studied geographical groups of sheep: FAE: Far Eastern; VOL: Volga; NCN: North Caucasian; SIB: Siberian; SOU: Southern; CEN: Central. Colors of branches correspond to different haplogroups (A, B, C, and D).
Genes 14 01701 g002
Genes 14 01701 g003
Figure 3. Graph of the mismatch distribution of Russian sheep from different geographical regions: (A) Southern; (B) North Caucasian; (C) Siberian; (D) Volga; (E) Far Eastern; (F) Central. The x and y axes represent the number of pairwise differences and the relative frequency of pairwise comparisons, respectively.
Figure 3. Graph of the mismatch distribution of Russian sheep from different geographical regions: (A) Southern; (B) North Caucasian; (C) Siberian; (D) Volga; (E) Far Eastern; (F) Central. The x and y axes represent the number of pairwise differences and the relative frequency of pairwise comparisons, respectively.
Genes 14 01701 g003
Table 1. Groups of studied sheep.
Table 1. Groups of studied sheep.
Geographical Regionn *Group Composition (Breed)
Far Eastern (FAE)8Buubei (n = 5)
Baikal fine-fleeced (n = 3)
Volga (VOL)5Tsigai (n = 5)
North Caucasian (NCN) 47Dagestan Mountain (n = 4)
Karachaev (n = 10)
Lezgin (n = 7)
Manych Merino (n = 5)
North Caucasian (n = 5)
Soviet Merino (n = 3)
Stavropol (n = 5)
Tushin (n = 8)
Siberian (SIB)17Altai Mountain (n = 5)
Buryat (n = 2)
Kulundin (n = 5)
Mongol (n = 1)
Tuva short-fat-tailed (n = 4)
Southern (SOU)38Volgograd (n = 8)
Groznenk (n = 7)
Kalmyk fat-rumped (n = 5)
Karakul (n = 6)
Salsk (n = 5)
Edilbai (n = 7)
Central (CEN)20Kuibyshev (n = 3)
Kuchugur (n = 3)
Romanov (n = 8)
Russian Longhaired (n = 6)
* n = sample size.
Table 2. Primer sequences selected for amplification of sequences of complete mitochondrial genome of sheep (O. aries).
Table 2. Primer sequences selected for amplification of sequences of complete mitochondrial genome of sheep (O. aries).
PrimerNucleotide SequencesLength of Amplification Product, b.p.
mtDNA_For15′-AGTACGGCGTAAAGCGTGTT-3′2987
mtDNA_Rev15′-AATGGTGCTCGGTTTGTTTC-3′
mtDNA_For25′-GAAAAGGCCCAAACGTTGTA-3′3294
mtDNA_Rev25′-GATATTATGGCTCATACTATTCCTATATA-3′
mtDNA_For35′-TCCTATATCAACACCTATTCTGATTCTT-3′3525
mtDNA_Rev35′-GGAAGTCAGAATGCGATGGT-3′
mtDNA_For45′-ACACCAAACCCACGCTTATC-3′3183
mtDNA_Rev45′-AAAATTGATTGCTGCGATGGGT-3′
mtDNA_For55′-TGAACGAGTTCACAGCCGAA-3′2810
mtDNA_Rev55′-ATTGTAAGTGGTGGGGTTGG-3′
mtDNA_For65′-AGCAATTCCCATAGCCTCCT-3′3517
mtDNA_Rev65′-GGCTGTTGCGGTAGTACTCT-3′
Table 3. Indices of genetic diversity and neutrality of six regional groups of Russian local breeds based on mitogenomes.
Table 3. Indices of genetic diversity and neutrality of six regional groups of Russian local breeds based on mitogenomes.
Group 1n 2S 3K 4H 5Hd 6
(±SEM) 7		π 8
(±SEM)		Tajima’s D 9Fu’s Fs 10
FAE824697,89370.964
±0.077		0.00589
±0.00136		0.17425 ns133.582 ns13
VOL54920,40051.000
±0.126		0.00123
±0.00017		−1.00076 ns130.544 ns13
NCN4746878,247460.999
±0.005		0.00471
±0.00052		−0.97712 ns13−8.987 *11
SIB1716656,191171.000
±0.020		0.00339
±0.00026		0.58964 ns13−2.097 **12
SOU3841455,092381.000
±0.006		0.00332
±0.00041		−1.67111 ns13−10.572 *11
CEN2021435,926190.995
±0.018		0.00216
±0.00059		−1.68045 ns13−2.62 **12
Overall13580462,6151280.9992
±0.010		0.00378
±0.00027		−1.91342 *11−31.631 *11
1 Groups of sheep: FAE: Far Eastern; VOL: Volga; NCN: North Caucasian; SIB: Siberian; SOU: Southern; CEN: Central; 2 n == sample number; 3 S = variable sites; 4 K = the average number of nucleotide differences; 5 H = the number of haplotypes; 6 Hd = haplotype diversity; 7 SEM = standard error mean; 8 π = nucleotide diversity; 9 Tajima’s D = Tajima’s neutrality test; 10 Fu’s Fs = neutrality test; *11 = indicates statistical significance p < 0.1, **12 = indicates statistical significance p < 0.05; 13ns = non-significant.
Table 4. Pairwise FST values calculated for the studied geographical groups of sheep.
Table 4. Pairwise FST values calculated for the studied geographical groups of sheep.
GroupsFAE 1VOL 2SIB 3NCN 4SOU 5CEN 6
FAE 10
VOL 20.392430
SIB 30.097070.299060
NCN 40.160990.048730.057810
SOU 50.298470.032480.108040.009370
CEN 60.439600.008180.280910.072880.040580
1 FAE: Far Eastern; 2 VOL: Volga; 3 NCN: North Caucasian; 4 SIB: Siberian; 5 SOU: Southern; 6 CEN: Central.
Table 5. AMOVA results for the studied groups of sheep.
Table 5. AMOVA results for the studied groups of sheep.
Source of
Variation		Freedom Degrees
d.f.		Sum of
Squares,
SS		Variance
Components,
VC		Percentage
of Variation
V %
Among
populations		2249.32210.2023422.13
Within
populations		27969.31235.9004477.87
Total291218.63346.10277
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

Koshkina, O.; Deniskova, T.; Dotsev, A.; Kunz, E.; Selionova, M.; Medugorac, I.; Zinovieva, N. Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences. Genes 2023, 14, 1701. https://doi.org/10.3390/genes14091701

AMA Style

Koshkina O, Deniskova T, Dotsev A, Kunz E, Selionova M, Medugorac I, Zinovieva N. Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences. Genes. 2023; 14(9):1701. https://doi.org/10.3390/genes14091701

Chicago/Turabian Style

Koshkina, Olga, Tatiana Deniskova, Arsen Dotsev, Elisabeth Kunz, Marina Selionova, Ivica Medugorac, and Natalia Zinovieva. 2023. "Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences" Genes 14, no. 9: 1701. https://doi.org/10.3390/genes14091701

APA Style

Koshkina, O., Deniskova, T., Dotsev, A., Kunz, E., Selionova, M., Medugorac, I., & Zinovieva, N. (2023). Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences. Genes, 14(9), 1701. https://doi.org/10.3390/genes14091701

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